Study on miRNAs in Pan-Cancer of the Digestive Tract Based on the Illumina HiSeq System Data Sequencing

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.

were estimated to occur in 2018. ey were colorectal cancer (6.1%), stomach cancer (5.7%), and esophagus cancer (3.2%) with 1,096,601, 1,033,701, and 572,034 incidences, respectively [2]. e prognosis of cancers of the digestive tract is poor because of the advanced stage at the time of being diagnosed.
us, early and timely diagnosis is of great importance for patients with cancers of the digestive tract which determines whether the patient can receive radical surgical treatment or not [3,4]. erefore, reliable biomarkers are urgently required for detection of the cancers of the digestive tract in primary stages to improve the adverse outcome of patients [5]. In this study, we aim to find potential biomarkers for early diagnosis and explore their molecular mechanisms in cancers of the digestive tract.
MicroRNA (miRNA) is an endogenous noncoding RNA with a length of about 22 nucleotides [6]. It plays an oncogene or tumor suppressor role in biological function by regulating target genes, which can inhibit the translation of protein via binding directly to specific sequences of target gene transcription through regulatory mechanisms including positive regulation of mRNA targets and induction of direct degradation of mRNA [7]. e dysregulated miRNA expression is correlated with human diseases, and some miRNAs also play a causal role in the biological processes of the tumor microenvironment [8]. All the above data indicated that miRNAs have the unique advantages as a biomarker because of their sensitive regulatory function and superior stability. More and more evidence has demonstrated that many key miRNAs participated in the initiation and progression of cancers of the digestive tract [9][10][11]. erefore, the purpose of our study was to find the diagnostic and prognostic value of key miRNAs, and the biological function of target genes in cancers of the digestive tract needs to be better elucidated.
In this study, miRNA and mRNA expression profiles associated with ESCA, GAC, COAD, and READ were selected from e Cancer Genome Atlas (TCGA) database. TCGA includes comprehensive cancer genome data such as mutation, copy number variation, mRNA expression, miRNA expression, and methylation data. After identifying differentially expressed miRNAs (DEMIs), we estimated their diagnostic and prognostic value. Furthermore, we take the intersection of differentially expressed genes (DEGs) and the target genes predicted by miRWalk 2.0 and established miRNA-mRNA regulatory networks. Meanwhile, the relationship between miRNAs and their target genes was explored by constructing the regulated network. Furthermore, we investigated the target genes through bioinformatics methods including miRNAs-gene regulatory networks, 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 molecular mechanisms of tumor development and give a new insight into cancers of the digestive tract molecular study.

Collection of DEMIs and DEMs in Cancers.
e cancer genome data collected by TCGA can be divided into three levels. Level 1 represents the original sequencing data (fasta, fastq, etc.), level 2 comprises compared bam files, and level 3 comprises processed and standardized data. e Illumina HiSeq system high-throughput miRNA and mRNA expression profile data (level 3) and corresponding survival time clinicopathological data were obtained from TCGA in 2018, which included ESCA, GAC, COAD, and READ tissue samples and corresponding normal tissue samples. Genes with |log 2 fold change (FC)| > 1 as well as P value <0.05 were considered to be DEMIs and differentially expressed mRNAs (DEMs). Since the data are publicly available, additional approval from the ethics committee was not needed in the present study.

Identification of DEMIs and DEMs in the Cancers of the Digestive Tract.
e edgeR package in Bioconductor was used to screen the DEMIs and DEMs with |logFC| > 1 and adjusted P value < 0.05. We filtered the miRNAs from the above with |logFC| > 3 and adjusted P value < 0.05, and then we took the intersection in the upregulated or downregulated miRNAs, respectively, in the cancers of the digestive tract.

Screening of the miRNAs with Prognostic Value.
e obtained DEMIs in ESCA, GAC, COAD, and READ were matched with a corresponding survival time and stage. After removing more than 10% of the missing values in above data, DEMIs were analyzed by the univariate Cox risk regression model to screen prognostic factors. Subsequently, we performed intersection of DEMIs associated with the prognosis of ESCA, GAC, COAD, and READ which were identified through the above univariate Cox risk regression model analysis, and DEMIs with significant prognostic value in cancers of the digestive tract were obtained. Kaplan-Meier analysis was used to demonstrate their prognostic significance.

Screening of the miRNAs with Diagnostic Value.
DEMIs were extracted (|logFC| > 3, P < 0.05) from the above screened differentially expressed miRNAs in ESCA, GAC, COAD, and READ. e diagnostic value of miRNAs was evaluated by ROC curve analysis performed via the MedCalc software.

Prediction of the Target Genes of miRNAs.
e potential target genes of miRNAs were predicted by miRWalk 2.0, which is a cross-prediction website that provides experimentally validated miRNA-target genes as well as target genes predicted by binding specific 3′-end-untranslated regions of specific mRNAs and other methods to predict target genes for miRNA. It contains 12 silico databases (miRWalk, Microt4, miRanda, miRBridge, miRDB, miR-Map, miRNAMap, PicTar2, PITA, RNA22, RNAhybrid, and TargetScan). e source version and prediction mechanism of these databases are listed in Table 1.

Construction of the miRNA-mRNA Regulatory Networks.
Only the genes predicted by 5 or more databases were regarded as potential target genes for each DEMI. e intersection between the predicted genes mentioned above and DEMs of TCGA was considered target genes of DEMIs. Based on the principle of miRNAs negatively regulating their target genes, the miRNA-mRNA regulatory networks were constructed using Cytoscape 3.4.0 which is a public source network visualization as well as visualization analysis tool. e core of the software is the network, and the nodes are genes, proteins, or molecules, and the connection is the interaction between these biological structures.

PPI Network Construction.
e STRING database (http://www.string-db.org/) is a large open-source database for searching the direct physical interaction and indirect functional correlation between known proteins and predictive proteins, with 2031 species, 9.6 million protein information, and 13.8 million protein interactions. Its connection lines use different color identifications to distinguish known interactions such as experimentally validated and predicted interactions and others. We constructed the PPI network of the target genes of DEMIs by using it.

GO Annotation, KEGG Pathway Analysis, and RIDA
Analysis. To further investigate the functional annotation as well as signal pathways related to the target genes of DEMIs, we conducted GO annotation analysis and KEGG pathway analysis using DAVID (https://david.ncifcrf.gov/), a database with integrated biological information and analytical tools. GO annotation analysis contains 3 categories, including molecular function (MF), biological process (BP), and cellular component (CC), and the hierarchical structure of functional information was described in every catalog. KEGG is second to none in the study of molecular signaling pathways, which includes the GENES database with complete genomic information and sequences as well as the PATHWAY database with more complete functional information. We also performed a gene group RIDA analysis which includes related diseases, core networks, and related transcription factors by applying GCBI (https://www.gcbi. com.cn/). GCBI is a platform which contains more than 1.2 million samples of 20 TB of information data, and it provides many services such as genetic analysis, data analysis, and molecular medical information solution.

Validation by Gene Expression Omnibus DataSets.
e clinical roles of the six miRNAs were validated by collecting the relevant miRNA microarrays from Gene Expression Omnibus (GEO) DataSets. Differences in miRNA expression levels between different groups were assessed using Student's t-test. e ROC curve analysis of the miRNAs was performed by the MedCalc software.

DEMIs and DEMs in Cancers of the Digestive Tract on
TCGA Database. After obtaining miRNA expression data of ESCA (90 cases and 13 normal samples), GAC (446 cases and 45 normal samples), COAD (457 cases and 8 normal samples), and READ (162 cases and 3 normal samples), the R language edgeR package was used for difference analysis and screening (|logFC| > 1, P < 0.05). e results are shown in Figure 1. In the same way, mRNA expression data of ESCA (80 cases and 11 normal samples), GAC (375 cases and 32 normal samples), COAD (480 cases and 41 normal samples), and READ (167 cases and 10 normal samples) were analyzed. Because of the large amount of data, the DEMs were shown in a more intuitive volcanic map instead of a list ( Figure 2).  Table 2). After taking the intersection, miR-7-3 was found in ESCA and GAC intersection, miR-328 in GAC and COAD intersection, and miR-323a in COAD and READ intersection ( Figure 3). We found that the expression level of these 3 miRNAs associated with prognosis in cancers of the digestive tract had significant difference compared with each control sample group. e expression of miR-7-3 was greatly higher than that in four normal sample groups, the expression of miR-328 was all obviously lower than that in four normal sample groups, and the expression of miR-323a was lower in COAD and READ than that in the normal sample group ( Figure 4). Subsequently, we performed survival analysis of miR-7-3, miR-328, and miR-323a by Kaplan-Meier analysis, and the results are shown in Figure 5. It was shown that the miR-7-3 high-expression group in ESCA had a lower survival rate than the low-expression group, while the highexpression group of miR-7-3 in GAC had a higher survival rate than the low one, which suggested miR-7-3 was a risk factor in ESCA and a protective factor in GAC. In COAD and GAC, the results indicated that the high-expression group of miR-328 had a lower survival rate than the low one, and miR-328 was a risk factor in COAD and GAC. e K-M survival curve of miR-323a in COAD and READ showed that the miR-323a high-expression group had a higher survival rate than the low one, and miR-323a was a protective factor in READ. e HR value was 1.2207 > 1, suggesting miR-323a as a risk factor in COAD. Meanwhile, the prognostic value of a diversity of clinicopathological parameters was also explored. It was demonstrated by univariate Cox regression analysis that gender, disease stage, N stage, M stage, lymphatic invasion, venous invasion, treatment outcome, neoplasm recurrence, residual tumor, and tumor status were associated with the prognosis of colon cancer (Supplementary Table 1).

Diagnostic Significance of the DEMIs.
A total of 31 miRNAs (24 upregulated and 7 downregulated) in ESCA, 71 miRNAs (65 upregulated and 6 downregulated) in GAC, 200 miRNAs (146 upregulated and 54 downregulated) in COAD, and 225 miRNAs (139 upregulated and 86 downregulated) in READ were considered DEMIs by the standard of |logFC| > 3 and P < 0.05 by the edgeR package. One upregulated miRNA (miR-552) and two downregulated miRNAs (miR-490 and miR-133a-2) were identified as DEMIs according to the intersection taken between upregulated and downregulated miRNAs in the cancers of the digestive tract ( Figure 6). We found that the expression level of these 3 miRNAs in cancers of the digestive tract had significant difference compared with each control sample group. e expression of miR-552 was greatly higher than that in the normal sample group, and the expression of miR-490 and miR-133a-2 was obviously lower than that in the normal sample group (Table 3; Figure 7). We evaluated the diagnostic value of the three DEMIs by ROC curve analysis. e AUC of the three miRNAs is shown in Figure 8. e results indicated that three miRNAs exhibited high diagnostic value in READ and moderate diagnostic value in ESCA, GAC, and COAD. Interestingly, the AUC value of the three miRNAs was all biggest in READ, and miR-552 had the highest AUC of 1.0 in READ, which indicated favorable diagnostic value of the three miRNAs in READ.

Biological Function Analysis of Overlapping Genes.
To better understand the biological functions of the 281 overlapping genes, we performed the GO annotation and KEGG pathway analyses. All significant results of GO function analysis with 281 overlapping genes are listed in Table 4.
ese target genes mainly participate in the biological process (BP) including positive regulation of cell proliferation, cell migration, positive regulation of gene expression, and positive regulation of cell migration. For the term "cellular component" (CC), these genes were involved in the plasma membrane, extracellular region, cell surface, and cell junction. Moreover, the target genes were mainly involved in heparin binding, actin binding, growth factor activity, and extracellular matrix structural constituent classified by molecular function (MF) analysis. Furthermore, 15 significant pathways among these 281 target genes were obtained by KEGG pathway analysis (P < 0.05), and the  results indicated the target genes were significantly enriched in the PI3K-Akt signaling pathway, Rap1 signaling pathway, Ras signaling pathway, focal adhesion, and pathways in cancer (Table 5; Figure 12). ese target genes were also performed for cluster radar analysis, including related diseases, core networks, and related transcription factors in RIDA analysis by GCBI. From the analysis result of related diseases, it was suggested that these genes were associated with tumor metastasis, adenocarcinoma, GAC, and colorectal cancer (CRC) (Figure 13). Core network analysis indicated the relationships of genes and gene proteins with interaction, activation, inhibition, phosphorylation, and ubiquitination ( Figure 14). Transcription factor analysis showed that the transcription factor regulation contains DNA binding, activation, and inhibition ( Figure 15).

Discussion
Cancer of the digestive tract is a kind of cancer with high mortality due to its poor prognosis. Increasing evidence has demonstrated that miRNAs played a significant role in diverse tumors and were involved in many biological   processes, for instance, proliferation, apoptosis, differentiation, and metabolism [7]. In addition, miRNAs participate in carcinogenesis and progression of malignant tumors by regulation of cell cycle progression, angiogenesis, cell migration, and invasion [12]. A large number of researches have showed that miRNAs play a vital role and could be used as an independent prognosis factor in cancers [13][14][15][16]. Some exosomal miRNAs are considered stable biomarkers of tumor disease and play a critical role in carcinogenesis and cancer progression. Exosomal miR-6803-5p could act as a biomarker for diagnosis and prognosis through detecting the level of exosomal miRNA in the serum of the CRC patient [17]. Furthermore, Fu et al. [18] indicated that the circulating miR-17-5p and miR-92a-3p were significantly related to the pathological stage and grade in patients with CRC. us, miRNAs have unique advantages to serve as novel biomarkers for cancer diagnosis because of their important regulatory functions, superior stability, and relatively small amounts.
In recent years, more and more studies had been reported from the perspective of pan-cancer [19,20]. Analysis of scientific problems across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in different tumors. Chen et al. presented the characterization of a large number of expressed enhancers in a genome-wide analysis of 8928 tumor samples of 33 cancer types [21]. It was rarely reported from the perspective of pan-cancer to analyze the prognostic and diagnostic value of miRNA in digestive tract tumors. We present the prognostic and diagnostic value of microRNAs in digestive tract tumors from the perspective of pan-cancer.
Our study found that miR-7-3, miR-328, and miR-323a were of significant prognostic value in four cancers of the digestive tract. Recently, Huo found that miR-7-3 was an independent prognostic factor in gastric adenocarcinoma (GAC) by analyzing the clinical data of 393 cases of TCGA and the expression data of 1881 miRNAs via Cox regression and univariate analysis. e results showed that the survival time of the miR-7-3 group with high expression was significantly longer than that of the miR-7-3 group with low one (P � 0.01) [22], which was consistent with our study. In addition, a study about cervical cancer showed that hsa-miR-7-3 was a potential prognostic biomarker by discovering that increased expression of hsa-miR-7-3 was significantly associated with shortened survival [23], and functional enrichment analysis showed poor prognosis of hsa-miR-7-3related target genes, which were mainly enriched in the AMPK pathway.
Also, circulating miRNAs can be applied to predict disease progression of primary CRC, including the serum (c) Figure 9: Venn map of predictive target genes of differential miRNAs and common differential genes. 49 target genes were obtained by miR-490-3p, 63 target genes were obtained by miR-133a-3p, and 179 target genes were obtained by miR-552-3p. e light blue color represents the potential target genes for the prediction, the light purple color represents the up-or downregulated common differential genes, and the dark purple color represents the overlapping genes. miR-328-3p marker, which was correlated with the overall survival (OS) and disease-free survival (DFS) of phase I-III CRC, and the curative effect of adjuvant chemotherapy was also related with the expression of miR-328-3p [24], which was consistent with our results in COAD. Furthermore, miR-328 overexpression inhibits the invasion of side population (SP) cells in the CRC and reverses drug resistance by targeting matrix metallopeptidase 16 (MMP16) and ATPbinding cassette subfamily G member 2 (ABCG2) [25,26]. On the contrary, Xue et al. discussed the role of miR-328 as a biomarker to predict tumor recurrence by studying the luciferase activity in the simulant of miR-328 with different concentrations. It was found that miR-328 could target CD44, and they showed a negative regulatory relationship between them, which indicated that detecting the expression level of miR-328 can be used as a means to predict the prognosis of gastric cancer patients [27]. miR-328 has been found to be associated with the prognosis of a variety of tumors in previous studies. In addition to CRC [28], it also involved in non-small-cell lung carcinoma (NSCLC), nasopharyngeal carcinoma, breast cancer, acute myeloid leukemia, glioma, etc. [29][30][31][32][33]. is is similar to the result of Zhang et al. [27], who studied the correlation between five miRNAs (miR-191, miR-28-3p, miR-145, miR-328, and miR-18a) in the serum of NSCLC patients and the three-year overall survival rate of the two groups. e result showed that miR-328 was differentially expressed in paired tumors and normal tissues (P < 0.0001), that the target gene of miR-328 was mainly enriched in inflammatory/immune response pathways and pathways related to chemoradiotherapy and/ or targeted treatment resistance, and that the miR-328 signal could play a role as a prognostic factor for patients with advanced NSCLC [29,34]. Also, Ulivi et al. found that PRKCA, il-1, and c-raf1 were involved in the migration of cancer cells through the mRNA target analysis of miR-328 and that the downregulation of hsa-miR-328-3p and its target genes and transcription factors played a key role in the carcinogenesis of NSCLC [35]. Furthermore, in vitro functional tests of Ma cell lines (U251 and U87) showed that high expression of miR-328 could significantly inhibit the  proliferation, invasion, and migration of two types of glial cells, and miR-328 could inhibit the invasion and proliferation of malignant cells, indicating a good prognosis for glioma [36]. Similarly, Wu et al. established the expression profile of mRNAs related to miR-328 with 60 GBM samples and identified the gene set closely related to cell mitosis. It was found that the ectopic expression of miR-328 blocked the spreading, adhesion, and proliferation of U87 cells and induced the arrest of the U87 cell division cycle, suggesting the decreased expression of miR-328 in high-level glioma associated with poor survival [33]. In a word, miR-328 can be applied to the diagnosis, prognosis, and treatment of CRC metastasis. Nevertheless, further studies are needed to verify its effectiveness, such as epidemiological and clinical researches.
As for the prognosis factor of miR-323a, there is still no research showing the intrinsic relationship between miR-323a and digestive tract tumors. Our study is the first one to suggest that miR-323a was a protective factor in CRC. But other cancer studies have shown that miR-323a can be used as a prognostic marker for acute coronary syndrome (ACS), cervical cancer, and glioblastoma [37][38][39]. Pilbrow et al. investigated the level of circulating miRNAs in patients with ACS and studied the relationship between circulating miRNAs and neurohormones, cardiac structure, and function as well as the 5-year follow-up period. It was found that miR-323-3p circulating in ACS patients was significantly increased at least one year later and could be a stable biomarker of ACS [37]. In addition, the increased expression level of miR-323-3p may contribute to the activation of metastatic potential of cervical cancer and may be involved in cell adhesion, cytoskeleton remodeling, and cell proliferation and migration by targeting metastasis-related genes to regulate the steps of the metastasis cascade [38]. In addition, as a prognostic indicator for glioblastoma patients, it has been found that the recurrence of glioblastoma relates to the genes concerning migration and adhesion targeted by miR-323, and the low level of miR-323 is associated with long OS. Furthermore, survival analysis showed that the expression level of miR-323 was increased in patients with nonrecurrence or long recurrence time [39].
Moreover, the diagnostic ability of miRNAs in cancers of the digestive tract has been widely reported. Zeng et al. [40] explored the diagnostic and prognostic value of miRNAs in ESCA based on the TCGA data and found the combined AUC of five miRNAs (miRNA-93, miRNA-21, miRNA-4746, miRNA-196a-1, and miRNA-196a-2) was 0.985. In addition, Wang et al. [41] identified four miRNAs were related to GAC pathogenesis based on 130 patients using qRT-PCR. e AUC values of miR-106a-5p and miR-19b-3p were 0.786 and 0.786, respectively. Similarly, Liu et al. [42] detected the expression of miR-21, miR-29a, miR-92a, miR-125b, and miR-223 in 85 patients with CRC utilizing qRT-PCR, evaluated their diagnostic value, and found the combined AUC of the miRNA was 0.952. Compared with above previous studies, the present study investigated the  diagnostic ability of 3 miRNAs (miR-490, miR-133a-2, and miR-552) in pan-cancer of the digestive tract together. miR-490 has been reported in many digestive tract tumors, such as GAC [43,44] and CRC [45][46][47]. Qu et al. [43] detected the expression level of miR-490-3p in GAC patients by qRT-PCR. e result showed that miR-490-3p was downregulated. e expression level of miR-490-3p in Helicobacter pylori-(HP-) positive GAC was significantly associated with lymph node metastasis and tumor differentiation. In addition, 5-year overall survival of GAC patients with low miR-490-3p expression was shortened in the HP-positive cohort, and multivariate analysis indicated that low expression of miR-490-3p was an independent risk factor in HP-positive GAC. Also, Tang et al. [48] found that miR-490 expression was downregulated in cancers of the digestive tract by analyzing 1,765 tumor samples in TCGA. e study of Zheng et al. [45] revealed that hypermethylation of the miR-490-3p promoter downregulated the expression of miR-490-3p in colorectal cancer, which revealed that miR-490-3p suppresses the proliferation of tumor cells by inducing apoptosis and inhibits the initiation of epithelial-tomesenchymal transition (EMT) in vitro and in vivo, indicating that its expression was continuously downregulated in CRC malignancy and CRC cell lines. e ability of migration and invasion in tumor cells was inhibited if miR-490-3p was oversuppressed, and the opposite effect was produced by knocking down miR-490-3p. Further studies found that it can inhibit CRC cell migration and invasion by targeting the TGF-β signaling pathway [46]. For miR-133a-2, however, the study of its function in tumors is rare in recent years, and only GAC [49], CRC [50], and bladder cancer [51] have been reported. Chen et al. [49] found that the AUC of miR-133a-2 in GAC was 0.905, and the sensitivity and specificity were 76.75% and 92.68% based on the TCGA data. miR-133a-2 showed significant correlation between race, tumor grade, tumor pathology, metastasis, and tumor stage (P < 0.05). In addition, miR-133a-2 was less expressed in the early stage (T1/T2 grade and I/II tumor stage) and was significantly higher in the late stage (P < 0.05), which suggested that miR-133a-2 may be a new potential predictor in the tumor stage. miR-552 function in CRC has been widely reported [52][53][54][55][56], and it has also been studied in lung tumor [57], breast tumor [58], and pancreatic cancer [59]. Wang and Liu [52] suggested that miR-552 was highly expressed in CRC through qRT-PCR experiments on 183 CRC patients and precancerous tissues, and high expression of miR-552 was significantly related to histological grade, lymph node metastasis, and TNM stage. Overall survival in patients with high expression levels of miR-552 was obviously lower than that in patients with low expression levels of miR-552 (p � 0.0036), which indicated high expression of miR-552 was an independent prognostic factor for CRC patients with poor prognosis through univariate and multivariate analysis. Cao et al. [53] found that miR-552 was upregulated in CRC tissues and cell lines by qRT-PCR experiments in 20 pairs of CRC patient tissues and adjacent control tissues and cell lines. miR-552 could be applied as a potential predictor and diagnostic biomarker for CRC. Wang et al. [54] also found that miR-552 was highly expressed in CRC by similar experiments and found that miR-552 promotes CRC metastasis by targeting ADAM28, which suggested miR-552 may serve as a new target for CRC therapy. miR-552 was highly expressed in CRC tissues in analysis of 25 cases of CRC and 10 cases of normal colorectal using qRT-PCR [56], and it was found to be related to the clinical stage and distant metastasis (p < 0.05).
In this study, 3 DEMIs (miR-490, miR-133a-2, and miR-552) exhibited high diagnostic ability in READ and moderate diagnostic ability in ESCA, GAC, and COAD. e expressions of miR-490 and miR-133a-2 were downregulated, and miR-552 was upregulated in cancers of the digestive tract. Our results were consistent with the former studies and validated through analysis of the Gene Expression Omnibus (GEO) DataSets. To summarize, the significant prognostic value of miR-490, miR-133a-2, and miR-552 in ESCA, GAC, COAD, and READ provided solid reasons for being therapeutic targets for patients with pancancer of the digestive tract, for which the novel prognosis and candidates may serve certain functions in process of tumor development, as predicted.
Increasing evidence has demonstrated that new diagnostic or therapeutic targets could be obtained, and new insight into the molecular mechanisms of tumor development may be provided through network regulation and signaling pathway analysis [40,45,53,60,61]. So we construct the miRNA-gene network and further explore the biological process and pathways in which the target genes could be participated. We perform the PPI network analysis, GO annotation and KEGG pathway analysis, and gene group RIDA analysis. 22 hub genes were identified by PPI network analysis, and G protein c subunit 7 (GNG7) was most active. GNG7 is the target gene of miR-552, but the relationship  between GNG7 and miR-552 is unclear. e team of surgeons at Kyushu University in Japan has explored the role of GNG7 in different tumors. ey found that GNG7 was downregulated in pancreatic cancer, gastric cancer, intrahepatic cholangiocarcinoma, esophageal cancer, and CRC [62][63][64][65]. ey also found that the expression of GNG7 was significantly repressed by miR-328 and the suppression of GNG7 was associated with great invasiveness and poor prognosis [65].
en, we supposed that the upregulated miR-552 might suppress the expression of GNG7 to promote the development of cancers of the digestive tract; however, this hypothesis needs to be confirmed with in vitro and in vivo experiments in future. In these hub genes, CXCL12 has been reported to be related to poor prognosis in ESCA [66]. We previously identified CXCL12 as one of the hub genes in ESCA with 187 cases based on TCGA data [40].
e GO annotation analysis of the 281 overlapping genes suggested that the target genes mainly participate in positive regulation of cell proliferation and play their molecular role by heparin binding. e KEGG pathways were also identified, such as PI3K-Akt signaling pathway, Rap1 signaling pathway, and pathways in cancer, which play a key role in the occurrence and development of cancers of the digestive tract. In the gene group RIDA analysis, the target genes were related to neoplasm, neoplasm metastasis, stomach neoplasm, and colorectal neoplasm. In addition, analysis of the core network indicated the relationship of activation, inhibition, and ubiquitination between these target genes. In the analysis of related transcription factors, the target genes were significantly enriched in DNA binding, activation, and inhibition. All the results of bioinformatics analysis demonstrated that these three miRNAs might influence the tumorigenesis and deterioration of tumors by targeting manifold genes and pathways. e target genes might play their role in the occurrence and development of cancers of the digestive tract through a number of biological processes, gene and gene interconnections, and the regulation of transcription factors. e results in this ERBB2 S10 0A 4 EL K1 H G F M G M T   present study gave us a direction to further explore the molecular mechanisms of the three miRNAs in the development of cancers of the digestive tract by using the experimental technique.

Conclusion
To summarize, we identified three DEMIs (miR-7-3, miR-328, and miR-323a) which are of significant prognostic value and three DEMIs (miR-490-3p, miR-133a-3p, and miR-552-3p) with high diagnostic ability in READ and moderate diagnostic ability in ESCA, GAC, and COAD by analyzing the miRNA expression profiles in a large number of samples from the TCGA database and GEO DataSet. A total of 281 target genes were also found, and they were used to perform a number of biological function analyses. e results indicated that variant genes and signaling pathways were targeted by these three DEMIs in cancers of the digestive tract. However, in vitro and in vivo experimentally functional verification was needed for target gene functional analysis in the present study. From the perspective of pancancer, our study may provide significant contribution to validation studies and bring the diagnostic and prognostic miRNAs or target genes to clinics in the future.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.