Identification of key microRNAs and their targets in exosomes of pancreatic cancer using bioinformatics analysis

Pancreatic cancer (PC) is one of the most lethal tumors, due to late diagnosis and limited surgical strategies. It has been reported that serum exosomal microRNAs (S-Exo-miRNAs) play a pivotal role as signaling molecules and serve as noninvasive diagnosis methods for PC. The combination of S-Exo-miRNAs with the corresponding target also plays an important role in the tumor microenvironment. Here we investigated S-Exo-miRNAs involved in PC. The gene expression profile was downloaded from the Gene Expression Omnibus (GEO) database. The analysis was carried out using GEO2R. The targets of differentially expressed serum exosomal miRNAs (DE-S-Exo-miRNAs) were predicted by 4 bioinformatic algorithms (miRanda, miRDB, miRWalk, and Targetscan). Further analysis with gene ontology (GO) and Kyoto Encyclopedia of Genomes pathway (KEGG) enrichment analyses were performed with Cytoscape software version 3.4.0. Subsequently, the interaction regulatory network of target genes was performed with the Search Tool for the Retrieval of Interacting Genes (STRING) database (http://www.string-db.org/) and visualized using Cytoscape software. We downloaded the gene expression profile GSE50632, which was based on an Agilent microarray GPL17660 platform containing 4 eligible samples. In total 467 DE-S-Exo-miRNAs were obtained, including 7 overexpressed miRNAs (1.50%), and 460 remaining underexpressed miRNAs (98.50%). The databases miRWalk, miRDB, miRanda, and TargetScan were used to predict their potential targets, which were subsequently submitted to Cytoscape software version 3.4.0 (www.cytoscape.org). Next the functional and pathway enrichment analysis were used for the KEGG pathway and GO categories analysis. The enrichment analysis identified the genes involved in such processes as developmental and negative regulation of multicellular organismal processes, regulation of anatomical structure morphogenesis, regulation of cell death, apoptotic processes and mitogen-activated protein kinase (MAPK) signaling pathway, transforming growth factor - beta (TGF -β) signaling pathway, cyclic adenosine monophosphate (cAMP) signaling pathway, and the phosphatidylinositol-3 kinases/Akt (PI3K-Akt) signaling pathway. Subsequently according to the protein–protein interaction (PPI) network, the top 10 genes were obtained. The enrichment analyses of the genes involved in a significant module revealed that these genes were related to the TGF-β signaling pathway. After reviewing the literature, we identified the apoptosis genes, and their corresponding miRNAs that have a relationship with apoptosis of the tumor. This analysis provides a comprehensive understanding of the roles of S-Exo-miRNAs and the related targets in the development of PC. Additionally, the present study provides promising candidate targets for early diagnosis and therapeutic intervention. However, these predictions require further experimental validation in future studies.