Integrated analysis identified novel miRNAs and mRNA in endometriosis.

Hepatocellular carcinoma (HCC) is the most common types of hepatic malignancies. This study aimed to better understand the pathogenesis of HCC and may help facilitate the improvement of the diagnostic of HCC. The mRNA and miRNA expression profiles of HCC, which was retrieved from The Cancer Genome Atlas database, and the circRNA expression profiles of HCC, which was retrieved from Gene Expression Omnibus database, were included in this study to perform an integrated analysis. The differentially expressed mRNAs (DEmRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed circRNAs (DEcircRNAs) were identified, and competing endogenous RNA (ceRNA) (DEcircRNA-DEmiRNA-DEmRNA) regulatory network was conducted. Functional annotation of host gene of DEcircRNAs and DEmRNAs in ceRNA regulatory network was performed. Quantitative real-time polymerase chain reaction validation of the expression of the selected DEmRNAs, DEmiRNAs, and DEcircRNAs was performed. A total of 2982 DEmRNAs, 144 DEmiRNAs, and 264 DEcircRNAs were obtained. The ceRNA network contained 61 circRNA-miRNA pairs and 1149 miRNA-mRNA pairs, including 48 circRNAs, 30 miRNAs, and 1149 mRNAs. Functional annotation of DEmRNAs in ceRNA regulatory network revealed that these DEmRNAs were significantly enriched in tryptophan metabolism, fatty acid metabolism, and pathways in cancer. Except for ARNT2 and hsa-miR-214-3p, expression of the others in the quantitative real-time polymerase chain reaction results was consistent with that in our integrated analysis, generally. We speculate that hsa_circRNA_104268/hsa-miR-214-3p/E2F2, hsa_circRNA_104168/hsa-miR-139-5p/HRAS, and hsa_circRNA_104769/hsa-miR-93-5p/JUN interaction pairs may play a vital role in HCC. This study expected to provide a novel insight into the pathogenesis and therapy of HCC from the circRNA-miRNA-mRNA network view.

In acute aristolochic acid nephropathy (AAN), aristolochic acid (AA) induces renal injury and tubulointerstitial fibrosis. However, the roles of microRNAs (miRNAs/miRs) and mRNAs involved in AAN are not clearly understood. The aim of the present study was to examine AA-induced genome-wide differentially expressed (DE) miRNAs and DE mRNAs using deep sequencing in mouse kidneys, and to analyze their regulatory networks. In the present self-controlled study, mice were treated with 5 mg/kg/day AA for 5 days, following unilateral nephrectomy. AA-induced renal injury and tubulointerstitial fibrosis were detected using hematoxylin and eosin staining and Masson's trichrome staining in the mouse kidneys. A total of 82 DE miRNAs and 4,605 DE mRNAs were identified between the AA-treated group and the self-control group. Of these DE miRNAs and mRNAs, some were validated using reverse transcription-quantitative PCR. Expression levels of the profibrotic miR-21, miR-433 and miR-132 families were significantly increased, whereas expression levels of the anti-fibrotic miR-122-5p and let-7a-1-3p were significantly decreased. Functions and signaling pathways associated with the DE miRNAs and mRNAs were analyzed using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG). A total of 767 DE pairs (in opposing directions) of miRNAs and their mRNA targets were identified. Among these, regulatory networks of miRNAs and mRNAs were analyzed using KEGG to identify enriched signaling pathways and extracellular matrix-associated pathways. In conclusion, the present study identified genome-wide DE miRNAs and mRNAs in the kidneys of AA-treated mice, as well as their regulatory pairs and signaling networks. The present results may improve the understanding of the role of DE miRNAs and their mRNA targets in the pathophysiology of acute AAN.

Objective. To develop a putative microRNA (miRNA) and messenger RNA (mRNA) regulatory network of Danggui Buxue decoction's (DGBXD) amelioration of idiopathic pulmonary fibrosis (IPF). Methods. The Gene Expression Omnibus (GEO) database was used to identify differentially expressed miRNAs (DE-miRNAs) and differentially expressed mRNAs (DE-mRNAs). Using miRNet, the predicted target genes of identified DE-miRNAs were estimated, and then the target genes of DE-miRNAs in IPF were comprehensively examined. The Enrichr database was used to conduct functional enrichment and pathway enrichment. Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) was employed to obtain the target genes of DGBXD as well as active compounds. A putative miRNA-mRNA regulatory network of DGBXD acting on IPF was developed by intersecting the target genes of DGBXD with the DE-miRNA target genes in IPF. A bleomycin-induced mouse model was established and used to perform histopathology as well as real-time quantitative polymerase chain reaction (qRT-PCR) analyses of some miRNA-mRNA pairs. Results. Fourteen upmodulated DE-miRNAs and six downmodulated DE-miRNAs were screened. The downstream target genes of upmodulated and downmodulated DE-miRNAs were predicted. Subsequently, 1160 upmodulated DE-mRNAs and 1427 downmodulated DE-mRNAs were identified. Then, target genes of DE-miRNAs comprising 49 downmodulated and 53 upmodulated target genes were further screened to perform functional enrichment and pathway enrichment analyses. Subsequently, 196 target genes of DGBXD were obtained from TCMSP, with six downregulated target genes and six upregulated target genes of DGBXD acting on IPF being identified. A promising miRNA-mRNA regulatory network of DGBXD acting on IPF was developed in this study. Moreover, mir-493 together with its target gene Olr1 and mir-338 together with Hif1a were further validated by qRT-PCR. Conclusion. This study proposed detailed possible processes of miRNA-mRNA modulatory axis in IPF and constructed a prospective IPF-related miRNA-mRNA modulatory network with the aim of alleviating IPF with DGBXD.

Lung adenocarcinoma (LUAD) is a highly malignant cancer. Although competing endogenous RNA (ceRNA)‐based profiling has been investigated in patients with LUAD, it has not been specifically used to study metastasis in LUAD. We found 130 differentially expressed (DE) lncRNAs, 32 DE miRNAs and 981 DE mRNAs from patients with LUAD in The Cancer Genome Atlas (TCGA) database. We analysed the functions and pathways of 981 DE mRNAs using the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Based on the target DE mRNAs and DE lncRNAs of DE miRNAs, we established an lncRNA‐miRNA‐mRNA ceRNA network, comprising 37 DE lncRNAs, 22 DE miRNAs and 212 DE mRNAs. Subsequently, we constructed a protein‐protein interaction network of DE mRNAs in the ceRNA network. Among all, DE RNAs, 5 DE lncRNAs, 5 DE miRNAs and 45 DE mRNAs were confirmed found to be associated with clinical prognosis. Moreover, 3 DE lncRNAs, 4 DE miRNAs and 9 DE mRNAs in the ceRNA network were associated with clinical prognosis. We further screened 3 DE lncRNAs, 3 DE miRNAs and 3 DE mRNAs using clinical samples. These DE lncRNAs, DE miRNAs and DE mRNAs in ceRNA network may serve as independent biomarkers of LUAD metastasis.

BackgroundThe role of circular RNAs (circRNAs) and microRNAs (miRNAs) in osteosarcoma (OS) development has not been fully elucidated. Further, the contribution of the immune response to OS progression is not well defined. However, it is known that circRNAs and miRNAs can serve as biomarkers for the diagnosis, prognosis, and therapy of many cancers. Thus, the aim of this study was to identify novel key serum biomarkers for the diagnosis and metastatic prediction of OS by analysis of immune cell infiltration and associated RNA molecules.MethodsHuman OS differentially expressed circRNAs (DEcircRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed mRNAs (DEmRNAs) were identified by analysis of microarray data downloaded from Gene Expression Omnibus (GEO) datasets. Further, characteristic patterns of OS-infiltrating immune cells were analyzed. On this basis, we identified statistically significant transcription factors. Moreover we performed pathway enrichment analysis, constructed protein–protein interaction networks, and devised competitive endogenous RNA (ceRNA) networks. Biological targets of the ceRNA networks were evaluated and potential OS biomarkers confirmed by RT-qPCR analysis of the patients’ serum.ResultsSeven differentially expressed circRNAs, 166 differentially expressed miRNAs, and 175 differentially expressed mRNAs were identified. An evaluation of cellular OS infiltration identified the highest level of infiltration by M0 macrophages, M2 macrophages, and CD8+ T cells, with M0 macrophages and CD8+ T cells as the most prominent. Significant patterns of tumor-infiltrating immune cells were identified by principal component analysis. Moreover, 185 statistically significant transcription factors were associated with OS. Further, in association with immune cell infiltration, hsa-circ-0010220, hsa-miR-326, hsa-miR-338-3p, and FAM98A were identified as potential novel biomarkers for OS diagnosis. Of these, FAM98A had the most promise as a diagnostic marker for OS and OS metastasis. Most importantly, a novel diagnostic model consisting of these four biomarkers (hsa-circ-0010220, hsa-miR-326, hsa-miR-338-3p, and FAM98A) was established with a 0.928 AUC value.ConclusionsIn summary, potential serum biomarkers for OS diagnosis and metastatic prediction were identified based on an analysis of immune cell infiltration. A novel diagnostic model consisting of these four promising serum biomarkers was established. Taken together, the results of this study provide a new perspective by which to understand immunotherapy of OS.

The aim of the present study was to elaborate the underlying pathogenesis of laryngeal squamous cell carcinoma (LSCC). Micro (mi) RNA and messenger (m) RNA expression profiling of patients with LSCC were downloaded from The Cancer Genome Atlas (TCGA) database. Differentially expressed miRNAs (DEMIs) and differentially expressed mRNAs (DEMs) were identified in LSCC compared to normal control tissues. The DEMs targeted by DEMIs were identified and the negative correlation between DEMs and DEMIs was subjected to visualization. The potential functions of DEMs targeted by DEMIs were annotated in Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database. A total of 663 dysregulated DEMs (449 upregulated and 214 downregulated) and 33 DEMIs (24 upregulated and 8 downregulated) were identified in LSCC compared with normal controls. 502 negative correlations between DEMIs and DEMs were identified and subjected to construct interaction network. In the network, hsa‑miR‑486, ‑34c, ‑206 and ‑182 had the highest connectivity with DEMs, and respectively regulated 39, 33, 28 and 27 DEMs. DEMs targeted by DEMIs were significantly enriched in signal transduction, actin binding and extracellular region of GO terms and focal adhesion and extracellular matrix‑receptor interaction of KEGG pathways. The present study may provide valuable information for understanding the potential oncogenesis mechanism in LSCC and provide the foundation work for diagnosis biomarkers and therapeutic targets for LSCC.

BACKGROUNDPsoriasis is a chronic inflammatory skin disease, the pathogenesis of which is more complicated and often requires long-term treatment. In particular, moderate to severe psoriasis usually requires systemic treatment. Psoriasis is also associated with many diseases, such as cardiometabolic diseases, malignant tumors, infections, and mood disorders. Psoriasis can appear at any age, and lead to a substantial burden for individuals and society. At present, psoriasis is still a treatable, but incurable, disease. Previous studies have found that microRNAs (miRNAs) play an important regulatory role in the progression of various diseases. Currently, miRNAs studies in psoriasis and dermatology are relatively new. Therefore, the identification of key miRNAs in psoriasis is helpful to elucidate the molecular mechanism of psoriasis.AIMTo identify key molecular markers and signaling pathways to provide potential basis for the treatment and management of psoriasis.METHODSThe miRNA and mRNA data were obtained from the Gene Expression Omnibus database. Then, differentially expressed mRNAs (DEmRNAs) and differentially expressed miRNAs (DEmiRNAs) were screened out by limma R package. Subsequently, DEmRNAs were analyzed for Gene Ontology and Kyoto Encyclopedia of Genes and Genomics functional enrichment. The “WGCNA” R package was used to analyze the co-expression network of all miRNAs. In addition, we constructed miRNA-mRNA regulatory networks based on identified hub miRNAs. Finally, in vitro validation was performed. All experimental procedures were approved by the ethics committee of Chinese PLA General Hospital (S2021-012-01). RESULTSA total of 639 DEmRNAs and 84 DEmiRNAs were identified. DEmRNAs screening criteria were adjusted P (adj. P) value < 0.01 and |logFoldChange| (|logFC|) > 1. DEmiRNAs screening criteria were adj. P value < 0.01 and |logFC| > 1.5. KEGG functional analysis demonstrated that DEmRNAs were significantly enriched in immune-related biological functions, for example, toll-like receptor signaling pathway, cytokine-cytokine receptor interaction, and chemokine signaling pathway. In weighted gene co-expression network analysis, turquoise module was the hub module. Moreover, 10 hub miRNAs were identified. Among these 10 hub miRNAs, only 8 hub miRNAs predicted the corresponding target mRNAs. 97 negatively regulated miRNA-mRNA pairs were involved in the miRNA-mRNA regulatory network, for example, hsa-miR-21-5p-claudin 8 (CLDN8), hsa-miR-30a-3p-interleukin-1B (IL-1B), and hsa-miR-181a-5p/hsa-miR-30c-2-3p-C-X-C motif chemokine ligand 9 (CXCL9). Real-time polymerase chain reaction results showed that IL-1B and CXCL9 were up-regulated and CLDN8 was down-regulated in psoriasis with statistically significant differences.CONCLUSIONThe identification of potential key molecular markers and signaling pathways provides potential research directions for further understanding the molecular mechanisms of psoriasis. This may also provide new research ideas for the prevention and treatment of psoriasis in the future.

Background Circular RNA (circRNA) is a noncoding RNA that forms a closed-loop structure, and its abnormal expression may cause disease. We aimed to find potential network for circRNA-related competitive endogenous RNA (ceRNA) in atrial fibrillation (AF). Methods The circRNA, miRNA, and mRNA expression profiles in the heart tissue from AF patients were retrieved from the Gene Expression Omnibus database and analyzed comprehensively. Differentially expressed circRNAs (DEcircRNAs), differentially expressed miRNAs (DEmiRNAs), and differentially expressed mRNAs (DEmRNAs) were identified, followed by the establishment of DEcircRNA-DEmiRNA-DEmRNA regulatory network. Functional annotation analysis of host gene of DEcircRNAs and DEmRNAs in ceRNA regulatory network was performed. In vitro experiment and electronic validation were used to validate the expression of DEcircRNAs, DEmiRNAs, and DEmRNAs. Results A total of 1611 DEcircRNAs, 51 DEmiRNAs, and 1250 DEmRNAs were identified in AF. The DEcircRNA-DEmiRNA-DEmRNA network contained 62 circRNAs, 14 miRNAs, and 728 mRNAs. Among which, two ceRNA regulatory pairs of hsa-circRNA-100053-hsa-miR-455-5p-TRPV1 and hsa-circRNA-005843-hsa-miR-188-5p-SPON1 were identified. In addition, six miRNA-mRNA regulatory pairs including hsa-miR-34c-5p-INMT, hsa-miR-1253-DDIT4L, hsa-miR-508-5p-SMOC2, hsa-miR-943-ACTA1, hsa-miR-338-3p-WIPI1, and hsa-miR-199a-3p-RAP1GAP2 were also obtained. MTOR was a significantly enriched signaling pathway of host gene of DEcircRNAs. In addition, arrhythmogenic right ventricular cardiomyopathy, dilated cardiomyopathy, and hypertrophic cardiomyopathy were remarkably enriched signaling pathways of DEmRNAs in DEcircRNA-DEmiRNA-DEmRNA regulatory network. The expression validation of hsa-circRNA-402565, hsa-miR-34c-5p, hsa-miR-188-5p, SPON1, DDIT4L, SMOC2, and WIPI1 was consistent with the integrated analysis. Conclusion We speculated that hsa-circRNA-100053-hsa-miR-455-5p-TRPV1 and hsa-circRNA-005843-hsa-miR-188-5p-SPON1 interaction pairs may be involved in AF.

BackgroundThe aim of our study was to use the differentially expressed mRNAs (DEmRNAs) and differentially expressed miRNAs (DEmiRNAs) to illustrate the underlying mechanism of hypoxia in liver cancer.MethodsIn this study, a cell model of hypoxia was established, and autophagy activity was measured with western blotting and transmission electron microscopy. The effect of hypoxia conditions on the invasion of liver cancer cell was evaluated. RNA sequencing was used to identify DEmRNAs and DEmiRNAs to explore the mechanism of hypoxia in liver cancer cells.ResultsWe found that autophagy activation was triggered by hypoxia stress and hypoxia might promote liver cancer cell invasion. In addition, a total of 407 shared DEmRNAs and 57 shared DEmiRNAs were identified in both HCCLM3 hypoxia group and SMMC-7721 hypoxia group compared with control group. Furthermore, 278 DEmRNAs and 24 DEmiRNAs were identified as cancer hypoxia-specific DEmRNAs and DEmiRNAs. Finally, we obtained 19 DEmiRNAs with high degree based on the DEmiRNA-DEmRNA interaction network. Among them, hsa-miR-483-5p, hsa-miR-4739, hsa-miR-214-3p and hsa-miR-296-5p may be potential gene signatures related to liver cancer hypoxia.ConclusionsOur study may help to understand the potential molecular mechanism of hypoxia in liver cancer.

Integrated analysis of miRNA and mRNA expression profiles associated with wooden breast myopathy in broiler chickens

BackgroundVentilator-induced diaphragmatic dysfunction (VIDD) is associated with weaning difficulties, intensive care unit hospitalization (ICU), infant mortality, and poor long-term clinical outcomes. The expression patterns of long noncoding RNAs (lncRNAs) and mRNAs in the diaphragm in a rat controlled mechanical ventilation (CMV) model, however, remain to be investigated.ResultsThe diaphragms of five male Wistar rats in a CMV group and five control Wistar rats were used to explore lncRNA and mRNA expression profiles by RNA-sequencing (RNA-seq). Muscle force measurements and immunofluorescence (IF) staining were used to verify the successful establishment of the CMV model. A total of 906 differentially expressed (DE) lncRNAs and 2,139 DE mRNAs were found in the CMV group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to determine the biological functions or pathways of these DE mRNAs. Our results revealed that these DE mRNAs were related mainly related to complement and coagulation cascades, the PPAR signaling pathway, cholesterol metabolism, cytokine-cytokine receptor interaction, and the AMPK signaling pathway. Some DE lncRNAs and DE mRNAs determined by RNA-seq were validated by quantitative real-time polymerase chain reaction (qRT-PCR), which exhibited trends similar to those observed by RNA-sEq. Co-expression network analysis indicated that three selected muscle atrophy-related mRNAs (Myog, Trim63, and Fbxo32) were coexpressed with relatively newly discovered DE lncRNAs.ConclusionsThis study provides a novel perspective on the molecular mechanism of DE lncRNAs and mRNAs in a CMV model, and indicates that the inflammatory signaling pathway and lipid metabolism may play important roles in the pathophysiological mechanism and progression of VIDD.

BackgroundA sight-threatening, cataract is a common degenerative disease of the ocular lens. This study aimed to explore the regulatory mechanism of age-related cataract (ARC) formation and progression.MethodsCataracts in Sprague Dawley rats were induced by adopting the method that injected selenite subcutaneously in the nape. We performed high-throughput RNA sequencing technology to identify the mRNA and microRNA(miRNA) expression profiles of the capsular membrane of the lens from Na2SeO3-induced and saline-injected Sprague Dawley rats. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were carried out to forecast the regulatory and functional role of mRNAs in cataracts by DAVID and Metascape. The protein–protein interaction(PPI) network of differentially expressed mRNA(DEmRNAs) was built via the STRING. Target miRNAs of hub genes were predicted by miRBD and TargetScan. Furthermore, differentially expressed miRNA(DEmiRNAs) were selected as hub genes’ targets, validated by quantitative real-time polymerase chain reaction(qRT-PCR), and a DEmiRNA-DEmRNA regulatory network was constructed via Cytoscape.ResultIn total, 329 DEmRNAs including 40 upregulated and 289 downregulated genes were identified. Forty seven DEmiRNAs including 29 upregulated and 18 downregulated miRNAs were detected. The DEmRNAs are involved in lens development, visual perception, and aging-related biological processes. A protein–protein interaction network including 274 node genes was constructed to explore the interactions of DEmRNAs. Furthermore, a DEmiRNA-DEmRNA regulatory network related to cataracts was constructed, including 8 hub DEmRNAs, and 8 key DEmiRNAs which were confirmed by qRT-PCR analysis.ConclusionWe identified several differentially expressed genes and established a miRNA-mRNA-regulated network in a Na2SeO3-induced Sprague Dawley rat cataract model. These results may provide novel insights into the clinical treatment of cataracts, and the hub DEmRNAs and key DEmiRNAs could be potential therapeutic targets for ARC.

Screening key genes and microRNAs in sepsis by RNA-sequencing

BackgroundBone is the most common site of metastatic breast cancer, and it is a leading cause of breast cancer-related death. This study aimed to explore bone metastasis-related long non-coding RNAs (lncRNAs) in breast cancer.MethodsFour mRNA datasets and two lncRNA datasets of bone metastasis, lung metastasis and liver metastasis of breast cancer were downloaded from Gene Expression Omnibus (GEO) database. The differentially expressed mRNAs (DEmRNAs) and lncRNAs (DElncRNAs) in group of bone metastasis vs lung metastasis and bone metastasis vs liver metastasis, as well as the overlap of the two groups, were identified. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and protein–protein interaction (PPI) network construction of DEmRNAs were conducted. The cis nearby-targeted DEmRNAs of DElncRNAs were obtained. Quantitative real-time polymerase chain reactions (qRT-PCR) was used to detect the expression levels of selected DEmRNAs and DElncRNAs. LOC641518-lymphoid enhancer-binding factor 1 (LEF1) pair was selected to verify its role in migration and invasion capability of breast cancer cells by wounding healing assay and transwell invasion assay.ResultsA total of 237 DEmRNAs were obtained in bone metastasis compared with both lung metastasis and liver metastasis. A total of three DElncRNAs in bone metastasis compared with both lung metastasis and liver metastasis were obtained. A total of seven DElncRNA-nearby-targeted DEmRNA pairs and 15 DElncRNA-nearby-targeted DEmRNA pairs in group of bone metastasis vs lung metastasis and bone metastasis vs liver metastasis, were detected, respectively. Four cis LncRNA-mRNA interaction pairs were identified, which are LOC641518-LEF1, FLJ35024-Very Low Density Lipoprotein Receptor (VLDLR), LOC285972-Retinoic Acid Receptor Responder 2 (RARRES2) and LOC254896-TNF receptor superfamily member 10c (TNFRSF10C). qRT-PCR using clinical samples from our hospital confirms the bioinformatics prediction. siRNA knocking down LOC641518 down-regulates LEF1 mRNA expression, and reduces the migration and invasion capability of breast cancer cells.ConclusionsWe concluded that four LncRNA-mRNA pairs, including LOC641518-LEF1, may play a central role in breast cancer bone metastasis.

Objective miRNA has gained attention as a therapeutic target in various malignancies. The proposal of this study was to investigate the biological functions of key miRNAs and target genes in cancers of the digestive tract which include esophageal carcinoma (ESCA), gastric adenocarcinoma (GAC), colon adenocarcinoma (COAD), and rectal adenocarcinoma (READ). Materials and Methods After screening differentially expressed miRNAs (DEMIs) and differentially expressed mRNAs (DEMs) in four digestive cancers from The Cancer Genome Atlas (TCGA) database, the diagnostic value of above DEMIs was evaluated by receiver-operating characteristic (ROC) curve analysis. Then, corresponding DEMIs' target genes were predicted by miRWalk 2.0. Intersection of predicted target genes and DEMs was taken as the target genes of DEMIs, and miRNA-mRNA regulatory networks between DEMIs and target genes were constructed. Meanwhile, the univariate Cox risk regression model was used to screen miRNAs with distinct prognostic value, and Kaplan–Meier analysis was used to determine their significance of prognosis. Furthermore, we performed bioinformatics methods including protein-protein interaction (PPI) networks, gene ontology (GO) annotation, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis, and gene group RIDA analysis by Gene-Cloud of Biotechnology Information (GCBI) to explore the function and molecular mechanisms of DEMIs and predicted target genes in tumor development. Results Eventually, 3 DEMIs (miR-7-3, miR-328, and miR-323a) with significant prognostic value were obtained. In addition, 3 DEMIs (miR-490-3p, miR-133a-3p, and miR-552-3p) and 281 target genes were identified, and the 3 DEMIs showed high diagnostic value in READ and moderate diagnostic value in ESCA, GAC, and COAD. Also, the miRNA-mRNA regulatory network with 3 DEMIs and 281 overlapping genes was successfully established. Functional enrichment analysis showed that 281 overlapping genes were mainly related to regulation of cell proliferation, cell migration, and PI3K-Akt signaling pathway. Conclusion The diagnostic value and prognostic value of significant DEMIs in cancers of the digestive tract were identified, which may provide a novel direction for treatment and prognosis improvement of cancers of the digestive tract.