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

Lai, Chun-hui; Liang, Xu-zhi; Liang, Xiu-yun; Ma, Sheng-jun; Li, Jun-guo; Shi, Ming-fang; Zhu, Xu; Lan, Hui-hua; Zeng, Jiang-hui

doi:https://doi.org/10.1155/2019/8016120

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Prevention and Early Detection of Gastrointestinal Cancers

View this Special Issue

Research Article | Open Access

Volume 2019 | Article ID 8016120 | https://doi.org/10.1155/2019/8016120

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

Chun-hui Lai,¹Xu-zhi Liang,²Xiu-yun Liang,¹Sheng-jun Ma,¹Jun-guo Li,¹Ming-fang Shi,¹Xu Zhu,¹Hui-hua Lan,³and Jiang-hui Zeng¹

Guest Editor: Guoku Hu

Received26 May 2019

Revised16 Aug 2019

Accepted06 Sept 2019

Published16 Oct 2019

Abstract

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.

1. Introduction

Cancer of the digestive tract is one of the most common malignancies in the world, and it is also a common high-risk malignancy in China, mainly including esophageal carcinoma (ESCA), gastric adenocarcinoma (GAC), colorectal cancer (CRC), and rectal adenocarcinoma (READ). According to the data of the National Cancer Center of China (NCCR) [1], the total number of newly diagnosed cancers of the digestive system in China was 2,089,500 cases and 1,565,500 deaths in 2015. Also, over 2.7 million new cases of cancers of the digestive tract and 1.8 million deaths were estimated to occur in 2018. They 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]. The prognosis of cancers of the digestive tract is poor because of the advanced stage at the time of being diagnosed. Thus, 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]. Therefore, 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]. The 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–11]. Therefore, 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 The 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.

2. Materials and Methods

2.1. Collection of DEMIs and DEMs in Cancers

The 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. The 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₂ fold change (FC)| > 1 as well as 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.

2.2. Identification of DEMIs and DEMs in the Cancers of the Digestive Tract

The edgeR package in Bioconductor was used to screen the DEMIs and DEMs with and adjusted value < 0.05. We filtered the miRNAs from the above with and adjusted value < 0.05, and then we took the intersection in the upregulated or downregulated miRNAs, respectively, in the cancers of the digestive tract.

2.3. Screening of the miRNAs with Prognostic Value

The 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.

2.4. Screening of the miRNAs with Diagnostic Value

DEMIs were extracted (, ) from the above screened differentially expressed miRNAs in ESCA, GAC, COAD, and READ. The diagnostic value of miRNAs was evaluated by ROC curve analysis performed via the MedCalc software.

2.5. Prediction of the Target Genes of miRNAs

The 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, miRMap, miRNAMap, PicTar2, PITA, RNA22, RNAhybrid, and TargetScan). The source version and prediction mechanism of these databases are listed in Table 1.

2.6. 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. The 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. The 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.

2.7. PPI Network Construction

The 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.

2.8. 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.

2.9. Validation by Gene Expression Omnibus DataSets

The 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. The ROC curve analysis of the miRNAs was performed by the MedCalc software.

3. Results

3.1. 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 (, ). The 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).

(a)

(b)

(c)

(d)

Figure 1

Differential expression of the miRNA heat map in four kinds of cancers of the digestive tract: (a) esophageal carcinoma (ESCA); (b) gastric adenocarcinoma (GAC); (c) colon adenocarcinoma (COAD); (d) rectal adenocarcinoma (READ). Among 157 DEMIs in ESCA tissues, 86 were upregulated and 71 were downregulated, while 252 and 86 were upregulated and downregulated, respectively, in 338 DEMIs in GAC tissues, and 205 were upregulated and 250 were downregulated among the 455 DEMIs in COAD tissues. Of the 339 READ DEMIs, 175 were upregulated and 164 were downregulated.

(a)

(b)

(c)

(d)

Figure 2

Differential expression of the mRNA volcano map in four kinds of cancers of the digestive tract: (a) esophageal carcinoma (ESCA); (b) gastric adenocarcinoma (GAC); (c) colon adenocarcinoma (COAD); (d) rectal adenocarcinoma (READ). Among 6926 DEMs in ESCA tissues, 2515 were upregulated and 4411 were downregulated, while 11558 and 7922 were upregulated and downregulated, respectively, in 3636 DEMs in GAC tissues, and 205 were upregulated and 250 were downregulated among the 455 DEMs in COAD tissues. Of the 9335 READ DEMs, 5144 were upregulated and 4191 were downregulated. The volcano map was drawn by the R language gplots package. The difference in expression of upregulated miRNAs (, ) was marked in red, and the difference in expression of downregulated miRNAs (, ) was marked in green.

3.2. Prognostic Significance of the DEMIs

We evaluated the prognostic value of three DEMIs by univariate Cox risk regression model analysis. The DEMIs above were further screened by matching with survival time and removing more than 10% missing values. 101 DEMIs related to ESCA, 140 DEMIs related to GAC, 296 DEMIs related to COAD, and 261 DEMIs related to READ were obtained. By identifying a number of prognostic values of miRNA molecular analysis by the univariate Cox risk regression model, DEMIs with prognostic values were obtained, among which 22 were related to ESCA, 22 to GAC, 22 to COAD, and 11 to READ (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. The 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 high-expression 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. The 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. The 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).

(a)

(b)

(c)

Figure 4

Boxplot of differential expression of (a) miR-7-3, (c) miR-328, and (b) miR-323a in gastrointestinal adenocarcinoma—colon adenocarcinoma (COAD), esophageal carcinoma (ESCA), rectal adenocarcinoma (READ), and gastric adenocarcinoma (GAC). In the boxplot, the abscissa represents the type of cancer and the ordinate represents the value of the expression. The brownish yellow color represents the tumor sample, and the green color represents the sample of the paracancerous control tissue.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5

K-M survival curve of miR-7-3, miR-328, and miR-323a in the cancers of the digestive tract—ESCA (esophageal adenocarcinoma), GAC (gastric adenocarcinoma), COAD (colon adenocarcinoma), and READ (rectal adenocarcinoma). K-M survival curve was drawn by the SPSS 19 software. The transverse coordinates represent the survival time (the number of days), and the ordinate represents the cumulative survival function. The blue color represents the low expression of the patient, and the green color represents the high expression of the patient. (a) ESCA miR-7-3. (b) GAC miR-7-3. (c) COAD miR-328. (d) GAC miR-328. (e) COAD miR-323a. (f) READ miR-323a.

3.3. 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 and 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. The 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. The AUC of the three miRNAs is shown in Figure 8. The 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.

(a)

(b)

(a)

(b)

(c)

Figure 7

Differential expressions of miR-490, miR-552, and miR-133a-2 in adenocarcinoma of the digestive tract—ESCA (esophageal adenocarcinoma), GAC (gastric adenocarcinoma), COAD (colon adenocarcinoma), and READ (rectal adenocarcinoma). In the boxplot, the abscissa represents the type of cancer and the ordinate represents the value of the expression. The brownish yellow color represents the tumor sample, and the green color represents the sample of the paracancerous control tissue. (a) miR-490. (b) miR-133a-2. (c) miR-552.

(a)

(b)

(c)

(d)

3.4. Prediction of Target Genes and Construction of miRNA-mRNA Regulatory Networks

We predicted potential target genes for miR-490-3p, miR-133a-3p, and miR-552-3p by miRWalk 2.0. Genes that can be predicted in 5 database resources at least were identified as potential target genes to intersect with DEMs obtained from TCGA. Since DEMIs and target genes are negatively correlated under sequence complementation, we intersected the predicted potential target genes of upregulated DEMIs with common downregulated DEMs obtained from TCGA and also the downregulated DEMIs’ predicted target genes with common upregulated DEMs. miR-490-3p and miR-133a-3p were intersected with upregulated DEMs. Meanwhile, miR-552-3p and downregulated DEMs were intersected. Finally, 281 overlapped genes in total were obtained (Figure 9). The miRNA-mRNA regulatory networks were established with miR-490-3p and its 49 target genes, miR-133a-3p and its 63 target genes, and miR-552-3p and its 179 target genes using the Cytoscape software (Figure 10). The PPI network was constructed with 281 target genes. According to the frequency of network connection, 22 hub genes (GNG7, CXCL12, GNAI1, LPAR1, PTGER3, APLN, CXCL11, GPER, CACNB2, COL11A1, COL1A1, COL4A1, COL5A2, COL5A3, GNAO1, SERPINH1, CDC6, F2RL2, KIT, P4HA3, PGF, and PTGFR) were obtained (Figure 11).

(a)

(b)

(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. The 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. (a) miR-490-3p. (b) miR-133a-3p. (c) miR-552-3p.

(a)

(b)

3.5. 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. These 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 (), 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). These 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).

3.6. Validation of These miRNAs Using Gene Expression Omnibus DataSets

The expression differences and expression trends of miR-7-5p, miR-323a, and miR-328 in CRC, ESCA, and GAC were verified by GSE89974, GSE114110, GSE93415, and GSE54397. Among them, miR-7-5p was a mature body of miR-7-3 (Figure 16). The expression differences and expression trends of miR-490-3p, miR-133a, and miR-552 in CRC, ESCA, and GAC were verified by GSE89974, GSE55856, GSE93415, GSE114110, GSE30070, and GSE66274 (Figure 17). Furthermore, their diagnostic values were verified by calculating the AUC value. The three areas under the ROC curves (AUCs) of miRNAs (miR-490-3p, miR-133a, and miR-552) were 0.758, 0.613, and 0.920 (CRC, ESCA, and GAC), 0.771, 1.000, and 0.761 (CRC, ESCA, and GAC), and 0.838, 0.809, and 0.715 (CRC, ESCA, and GAC) (Figure 18).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

Figure 18

Verification of the diagnostic value of the three differentially expressed miRNAs based on Gene Expression Omnibus (GEO) DataSets. The receiver-operating characteristic (ROC) curves were drawn using the MedCalc software. (a) CRC GSE89974. (b) ESCA GSE55856. (c) GAC GSE93415. (d) CRC GSE89974. (e) ESCA GSE114110. (f) GAC GSE30070. (g) CRC GSE89974. (h) ESCA GSE66274. (i) GAC GSE30070. AUC: area under the ROC curve; CRC: colorectal cancer; ESCA: esophageal carcinoma; GAC: gastric adenocarcinoma.

4. 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–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. Thus, 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. The 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 () [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-3-related 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 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 ATP-binding 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–33]. This 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. The result showed that miR-328 was differentially expressed in paired tumors and normal tissues (), 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–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. The 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–47]. Qu et al. [43] detected the expression level of miR-490-3p in GAC patients by qRT-PCR. The result showed that miR-490-3p was downregulated. The 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. The 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-to-mesenchymal transition (EMT) in vitro and in vivo, indicating that its expression was continuously downregulated in CRC malignancy and CRC cell lines. The 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 (). 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 (), 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–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 (), 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 ().

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. The 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 pan-cancer 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 γ subunit 7 (GNG7) was most active. GNG7 is the target gene of miR-552, but the relationship between GNG7 and miR-552 is unclear. The team of surgeons at Kyushu University in Japan has explored the role of GNG7 in different tumors. They found that GNG7 was downregulated in pancreatic cancer, gastric cancer, intrahepatic cholangiocarcinoma, esophageal cancer, and CRC [62–65]. They 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]. Then, 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]. The 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. The 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. The 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. The results in this 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.

5. 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. The 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 pan-cancer, 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

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by funding from the Nanning Scientific Research and Technical Development Project (20133182), Guangxi Health Department Research Project (Z2014458, Z2013401, and Z2013675), and Guangxi Natural Science Foundation Project (2018GXNSFAA050150).

Supplementary Materials

Supplementary Table 1: predictive values of related clinical parameters. The clinicopathological parameters from patients with colon cancer were obtained from the TCGA data portal (https://portal.gdc.cancer.gov/). These data were matched with miRNA expression values, and parameters of 431 colon cancer patients were obtained. The clinicopathological parameters include gender, age, disease stage, T stage, N stage, M stage, lymphatic invasion, venous invasion, perineural invasion, treatment outcome, neoplasm recurrence, residual tumor, BMI, radiation therapy, and tumor status. 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 Materials)

References

W. Chen, R. Zheng, P. D. Baade et al., “Cancer statistics in China, 2015,” CA: A Cancer Journal for Clinicians, vol. 66, no. 2, pp. 115–132, 2016.
View at: Publisher Site | Google Scholar
F. Bray, J. Ferlay, I. Soerjomataram, R. L. Siegel, L. A. Torre, and A. Jemal, “Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries,” CA: A Cancer Journal for Clinicians, vol. 68, no. 6, pp. 394–424, 2018.
View at: Publisher Site | Google Scholar
M.-M. He, F. Wang, Y. Jin et al., “Phase II clinical trial of S-1 plus nanoparticle albumin-bound paclitaxel in untreated patients with metastatic gastric cancer,” Cancer Science, vol. 109, no. 11, pp. 3575–3582, 2018.
View at: Publisher Site | Google Scholar
B. E. C. Fava, W. L. da Costa Jr., M. L. L. Medeiros et al., “Neoadjuvant intraperitoneal chemotherapy followed by radical surgery and HIPEC in patients with very advanced gastric cancer and peritoneal metastases: report of an initial experience in a western single center,” World Journal of Surgical Oncology, vol. 16, no. 1, p. 62, 2018.
View at: Publisher Site | Google Scholar
H. Tanaka, S. Hazama, M. Iida et al., “miR-125b-1 and miR-378a are predictive biomarkers for the efficacy of vaccine treatment against colorectal cancer,” Cancer Science, vol. 108, no. 11, pp. 2229–2238, 2017.
View at: Publisher Site | Google Scholar
M. F. Souza, H. Kuasne, M. C. Barros-Filho et al., “Circulating mRNAs and miRNAs as candidate markers for the diagnosis and prognosis of prostate cancer,” PLoS One, vol. 12, no. 9, Article ID e0184094, 2017.
View at: Publisher Site | Google Scholar
F. Lovat, N. Valeri, and C. M. Croce, “MicroRNAs in the pathogenesis of cancer,” Seminars in Oncology, vol. 38, no. 6, pp. 724–733, 2011.
View at: Publisher Site | Google Scholar
C. J. Cheng, R. Bahal, I. A. Babar et al., “MicroRNA silencing for cancer therapy targeted to the tumour microenvironment,” Nature, vol. 518, no. 7537, pp. 107–110, 2015.
View at: Publisher Site | Google Scholar
Z. Zhang and D.-Q. Dai, “MicroRNA-596 acts as a tumor suppressor in gastric cancer and is upregulated by promotor demethylation,” World Journal of Gastroenterology, vol. 25, no. 10, pp. 1224–1237, 2019.
View at: Publisher Site | Google Scholar
L. H. Yao, G. R. Wang, Y. Cai et al., “The expressions and diagnostic values of miR-18a and miR-21 in esophageal cancer,” Zhonghua Zhong Liu Za Zhi, vol. 41, no. 2, pp. 107–111, 2019.
View at: Google Scholar
X. Xi, M. Teng, L. Zhang, L. Xia, J. Chen, and Z. Cui, “MicroRNA-204-3p represses colon cancer cells proliferation, migration, and invasion by targeting HMGA2,” Journal of Cellular Physiology, 2019.
View at: Publisher Site | Google Scholar
L. Fabris, Y. Ceder, A. M. Chinnaiyan et al., “The potential of MicroRNAs as prostate cancer biomarkers,” European Urology, vol. 70, no. 2, pp. 312–322, 2016.
View at: Publisher Site | Google Scholar
P. Li, Q. Ou, T. A. Braciak, G. Chen, and F. S. Oduncu, “MicroRNA-192-5p is a predictive biomarker of survival for stage IIIB colon cancer patients,” Japanese Journal of Clinical Oncology, vol. 48, no. 7, pp. 619–624, 2018.
View at: Publisher Site | Google Scholar
T. Jiang, L. Ye, Z. Han et al., “miR-19b-3p promotes colon cancer proliferation and oxaliplatin-based chemoresistance by targeting SMAD4: validation by bioinformatics and experimental analyses,” Journal of Experimental & Clinical Cancer Research, vol. 36, no. 1, p. 131, 2017.
View at: Publisher Site | Google Scholar
J. Yang, D. Ma, A. Fesler et al., “Expression analysis of microRNA as prognostic biomarkers in colorectal cancer,” Oncotarget, vol. 8, no. 32, pp. 52403–52412, 2017.
View at: Publisher Site | Google Scholar
F. Tong, Y. Ying, H. Pan, W. Zhao, H. Li, and X. Zhan, “MicroRNA-466 (miR-466) functions as a tumor suppressor and prognostic factor in colorectal cancer (CRC),” Bosnian Journal of Basic Medical Sciences, vol. 18, no. 3, pp. 252–259, 2018.
View at: Publisher Site | Google Scholar
S. Yan, Y. Jiang, C. Liang et al., “Exosomal miR-6803-5p as potential diagnostic and prognostic marker in colorectal cancer,” Journal of Cellular Biochemistry, vol. 119, no. 5, pp. 4113–4119, 2018.
View at: Publisher Site | Google Scholar
F. Fu, W. Jiang, L. Zhou, and Z. Chen, “Circulating exosomal mir-17-5p and mir-92a-3p predict pathologic stage and grade of colorectal cancer,” Translational Oncology, vol. 11, no. 2, pp. 221–232, 2018.
View at: Publisher Site | Google Scholar
J. Liu, T. Lichtenberg, K. A. Hoadley et al., “An integrated TCGA pan-cancer clinical data resource to drive high-quality survival outcome analytics,” Cell, vol. 173, no. 2, pp. 400–416 e11, 2018.
View at: Google Scholar
X. Ma, Y. Liu, Y. Liu et al., “Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours,” Nature, vol. 555, no. 7696, pp. 371–376, 2018.
View at: Publisher Site | Google Scholar
H. Chen, C. Li, X. Peng et al., “A pan-cancer analysis of enhancer expression in nearly 9000 patient samples,” Cell, vol. 173, no. 2, pp. 386–399 e12, 2018.
View at: Publisher Site | Google Scholar
Q. Huo, “Analysis of expression profile of miRNA in stomach adenocarcinoma,” Journal of BUON, vol. 22, no. 5, pp. 1154–1159, 2017.
View at: Google Scholar
Y. Zeng, K.-X. Wang, H. Xu, and Y. Hong, “Integrative miRNA analysis identifies hsa-miR-3154, hsa-miR-7-3, and hsa-miR-600 as potential prognostic biomarker for cervical cancer,” Journal of Cellular Biochemistry, vol. 119, no. 2, pp. 1558–1566, 2018.
View at: Publisher Site | Google Scholar
D. Ji, M. Qiao, Y. Yao et al., “Serum-based microRNA signature predicts relapse and therapeutic outcome of adjuvant chemotherapy in colorectal cancer patients,” EBioMedicine, vol. 35, pp. 189–197, 2018.
View at: Publisher Site | Google Scholar
X. T. Xu, Q. Xu, J. L. Tong et al., “MicroRNA expression profiling identifies miR-328 regulates cancer stem cell-like SP cells in colorectal cancer,” British Journal of Cancer, vol. 106, no. 7, pp. 1320–1330, 2012.
View at: Publisher Site | Google Scholar
P. Tokarz and J. Blasiak, “The role of microRNA in metastatic colorectal cancer and its significance in cancer prognosis and treatment,” Acta Biochimica Polonica, vol. 59, no. 4, pp. 467–474, 2012.
View at: Publisher Site | Google Scholar
H.-G. Xue, A.-H. Yang, X.-G. Sun, Y.-Y. Lu, and Z.-B. Tian, “Expression of microRNA-328 functions as a biomarker for recurrence of early gastric cancer (EGC) after endoscopic submucosal dissection (ESD) by modulating CD44,” Medical Science Monitor, vol. 22, pp. 4779–4785, 2016.
View at: Publisher Site | Google Scholar
M. Xu, Y. Kuang, M. Wang, X. Han, and Q. Yang, “A microRNA expression signature as a predictor of survival for colon adenocarcinoma,” Neoplasma, vol. 64, no. 1, pp. 56–64, 2017.
View at: Publisher Site | Google Scholar
K. Wang, M. Chen, and W. Wu, “Analysis of microRNA (miRNA) expression profiles reveals 11 key biomarkers associated with non-small cell lung cancer,” World Journal of Surgical Oncology, vol. 15, no. 1, p. 175, 2017.
View at: Publisher Site | Google Scholar
S. Zhang, P. Li, L. Zhao, and L. Xu, “LINC00210 as a miR-328-5p sponge promotes nasopharyngeal carcinoma tumorigenesis by activating NOTCH3 pathway,” Bioscience Reports, vol. 38, no. 6, Article ID BSR20181168, 2018.
View at: Publisher Site | Google Scholar
N. Al-Othman, H. Hammad, and M. Ahram, “Dihydrotestosterone regulates expression of CD44 via miR-328-3p in triple-negative breast cancer cells,” Gene, vol. 675, pp. 128–135, 2018.
View at: Publisher Site | Google Scholar
L. Liu, R. Chen, Y. Zhang, W. Fan, F. Xiao, and X. Yan, “Low expression of circulating microRNA-328 is associated with poor prognosis in patients with acute myeloid leukemia,” Diagnostic Pathology, vol. 10, no. 1, p. 109, 2015.
View at: Publisher Site | Google Scholar
Z. Wu, L. Sun, H. Wang et al., “MiR-328 expression is decreased in high-grade gliomas and is associated with worse survival in primary glioblastoma,” PLoS One, vol. 7, no. 10, Article ID e47270, 2012.
View at: Publisher Site | Google Scholar
Y. Zhang, J. A. Roth, H. Yu et al., “A 5-MicroRNA signature identified from serum MicroRNA profiling predicts survival in patients with advanced stage non-small cell lung cancer,” Carcinogenesis, vol. 40, no. 5, pp. 643–650, 2018.
View at: Google Scholar
P. Ulivi, G. Foschi, M. Mengozzi et al., “Peripheral blood miR-328 expression as a potential biomarker for the early diagnosis of NSCLC,” International Journal of Molecular Sciences, vol. 14, no. 5, pp. 10332–10342, 2013.
View at: Publisher Site | Google Scholar
J. Yuan, Z. Zheng, Y. Zheng, X. Lu, L. Xu, and L. Lin, “microRNA-328 is a favorable prognostic marker in human glioma via suppressing invasive and proliferative phenotypes of malignant cells,” International Journal of Neuroscience, vol. 126, no. 2, pp. 145–153, 2016.
View at: Publisher Site | Google Scholar
A. P. Pilbrow, L. Cordeddu, V. A. Cameron et al., “Circulating miR-323-3p and miR-652: candidate markers for the presence and progression of acute coronary syndromes,” International Journal of Cardiology, vol. 176, no. 2, pp. 375–385, 2014.
View at: Publisher Site | Google Scholar
H. Ding, Y.-L. Wu, Y.-X. Wang, and F.-F. Zhu, “Characterization of the microRNA expression profile of cervical squamous cell carcinoma metastases,” Asian Pacific Journal of Cancer Prevention, vol. 15, no. 4, pp. 1675–1679, 2014.
View at: Publisher Site | Google Scholar
S. Qiu, S. Lin, D. Hu, Y. Feng, Y. Tan, and Y. Peng, “Interactions of miR-323/miR-326/miR-329 and miR-130a/miR-155/miR-210 as prognostic indicators for clinical outcome of glioblastoma patients,” Journal of Translational Medicine, vol. 11, no. 1, p. 10, 2013.
View at: Publisher Site | Google Scholar
J. H. Zeng, D. D. Xiong, Y. Y. Pang et al., “Identification of molecular targets for esophageal carcinoma diagnosis using miRNA-seq and RNA-seq data from the cancer genome atlas: a study of 187 cases,” Oncotarget, vol. 8, no. 22, pp. 35681–35699, 2017.
View at: Publisher Site | Google Scholar
N. Wang, L. Wang, Y. Yang, L. Gong, B. Xiao, and X. Liu, “A serum exosomal microRNA panel as a potential biomarker test for gastric cancer,” Biochemical and Biophysical Research Communications, vol. 493, no. 3, pp. 1322–1328, 2017.
View at: Publisher Site | Google Scholar
H.-N. Liu, T.-T. Liu, H. Wu et al., “Serum microRNA signatures and metabolomics have high diagnostic value in colorectal cancer using two novel methods,” Cancer Science, vol. 109, no. 4, pp. 1185–1194, 2018.
View at: Publisher Site | Google Scholar
M. Qu, L. Li, and W. C. Zheng, “Reduced miR-490-3p expression is associated with poor prognosis of Helicobacter pylori induced gastric cancer,” European Review for Medical and Pharmacological Sciences, vol. 21, no. 15, pp. 3384–3388, 2017.
View at: Google Scholar
B. Zhou, Y. Wang, J. Jiang et al., “The long noncoding RNA colon cancer-associated transcript-1/miR-490 axis regulates gastric cancer cell migration by targeting hnRNPA1,” IUBMB Life, vol. 68, no. 3, pp. 201–210, 2016.
View at: Publisher Site | Google Scholar
K. Zheng, X. Zhou, J. Yu et al., “Epigenetic silencing of miR-490-3p promotes development of an aggressive colorectal cancer phenotype through activation of the Wnt/β-catenin signaling pathway,” Cancer Letters, vol. 376, no. 1, pp. 178–187, 2016.
View at: Publisher Site | Google Scholar
X. Xu, R. Chen, Z. Li et al., “MicroRNA-490-3p inhibits colorectal cancer metastasis by targeting TGFbetaR1,” BMC Cancer, vol. 15, no. 1, p. 1023, 2015.
View at: Publisher Site | Google Scholar
X. Xu, X. Wu, S. Wu et al., “Study on miR-490-5p and miR-363 as novel biomarkers for the diagnosis of colorectal cancer,” Zhonghua Wei Zhonghua Wei Chang Wai Ke Za Zhi, vol. 17, no. 1, pp. 45–50, 2014.
View at: Google Scholar
S. Tang, W. K. Wu, X. Li et al., “Stratification of digestive cancers with different pathological features and survival outcomes by MicroRNA expression,” Scientific Reports, vol. 6, no. 1, Article ID 24466, 2016.
View at: Publisher Site | Google Scholar
S. Chen, J. Zhu, F. Yu, Y. Tian, S. Ma, and X. Liu, “Combination of miRNA and RNA functions as potential biomarkers for gastric cancer,” Tumor Biology, vol. 36, no. 12, pp. 9909–9918, 2015.
View at: Publisher Site | Google Scholar
W.-S. Chen, C.-M. Leung, H.-W. Pan et al., “Silencing of miR-1-1 and miR-133a-2 cluster expression by DNA hypermethylation in colorectal cancer,” Oncology Reports, vol. 28, no. 3, pp. 1069–1076, 2012.
View at: Publisher Site | Google Scholar
T. Yamasaki, H. Yoshino, H. Enokida et al., “Novel molecular targets regulated by tumor suppressors microRNA-1 and microRNA-133a in bladder cancer,” International Journal of Oncology, vol. 40, pp. 1821–1830, 2012.
View at: Publisher Site | Google Scholar
N. Wang and W. Liu, “Increased expression of miR-552 acts as a potential predictor biomarker for poor prognosis of colorectal cancer,” European Review for Medical and Pharmacological Sciences, vol. 22, no. 2, pp. 412–416, 2018.
View at: Google Scholar
J. Cao, X.-R. Yan, T. Liu et al., “MicroRNA-552 promotes tumor cell proliferation and migration by directly targeting DACH1 via the Wnt/β-catenin signaling pathway in colorectal cancer,” Oncology Letters, vol. 14, no. 3, pp. 3795–3802, 2017.
View at: Publisher Site | Google Scholar
J. Wang, H. Li, Y. Wang et al., “MicroRNA-552 enhances metastatic capacity of colorectal cancer cells by targeting a disintegrin and metalloprotease 28,” Oncotarget, vol. 7, no. 43, pp. 70194–70210, 2016.
View at: Publisher Site | Google Scholar
Z.-S. Xia, L. Wang, T. Yu et al., “MiR-5000-3p, miR-5009-3P and miR-552: potential microRNA biomarkers of side population cells in colon cancer,” Oncology Reports, vol. 32, no. 2, pp. 589–596, 2014.
View at: Publisher Site | Google Scholar
W. C. Chen, M. S. Lin, Y. L. Ye, H. J. Gao, Z. Y. Song, and X. Y. Shen, “microRNA expression pattern and its alteration following celecoxib intervention in human colorectal cancer,” Experimental and Therapeutic Medicine, vol. 3, no. 6, pp. 1039–1048, 2012.
View at: Publisher Site | Google Scholar
J. Kim, N. J. Lim, S. G. Jang, H. K. Kim, and G. K. Lee, “miR-592 and miR-552 can distinguish between primary lung adenocarcinoma and colorectal cancer metastases in the lung,” Anticancer Research, vol. 34, no. 5, pp. 2297–2302, 2014.
View at: Google Scholar
S.-K. Leivonen, K. K. Sahlberg, R. Mäkelä et al., “High-throughput screens identify microRNAs essential for HER2 positive breast cancer cell growth,” Molecular Oncology, vol. 8, no. 1, pp. 93–104, 2014.
View at: Publisher Site | Google Scholar
Y. G. Park, K. H. Lee, J. K. Lee et al., “MicroRNA expression pattern in intraductal papillary mucinous neoplasm,” The Korean Journal of Gastroenterology, vol. 58, no. 4, pp. 190–200, 2011.
View at: Publisher Site | Google Scholar
J. Liu, F. Liu, Y. Shi, H. Tan, and L. Zhou, “Identification of key miRNAs and genes associated with stomach adenocarcinoma from the cancer genome atlas database,” FEBS Open Bio, vol. 8, no. 2, pp. 279–294, 2018.
View at: Publisher Site | Google Scholar
S. Wu, F. Wu, and Z. Jiang, “Identification of hub genes, key miRNAs and potential molecular mechanisms of colorectal cancer,” Oncology Reports, vol. 38, no. 4, pp. 2043–2050, 2017.
View at: Publisher Site | Google Scholar
K. Shibata, M. Mori, S. Tanaka, S. Kitano, and T. Akiyoshi, “Identification and cloning of human G-protein γ7, down-regulated in pancreatic cancer,” Biochemical and Biophysical Research Communications, vol. 246, no. 1, pp. 205–209, 1998.
View at: Publisher Site | Google Scholar
T. Utsunomiya, H. Inoue, K. Taguchi, M. Shimada, K. Sugimachi, and M. Mori, “G protein γ7 expression as a new clinicopathological marker in patients with intrahepatic cholangiocarcinoma,” Archives of Surgery, vol. 137, no. 2, pp. 181–185, 2002.
View at: Publisher Site | Google Scholar
K. Shibata, S. Tanaka, T. Shiraishi, S. Kitano, and M. Mori, “G-protein gamma 7 is down-regulated in cancers and associated with p 27kip1-induced growth arrest,” Cancer Research, vol. 59, no. 5, pp. 1096–1101, 1999.
View at: Google Scholar
M. Ohta, K. Mimori, Y. Fukuyoshi et al., “Clinical significance of the reduced expression of G protein gamma 7 (GNG7) in oesophageal cancer,” British Journal of Cancer, vol. 98, no. 2, pp. 410–417, 2008.
View at: Publisher Site | Google Scholar
Y. Uchi, H. Takeuchi, S. Matsuda et al., “CXCL12 expression promotes esophageal squamous cell carcinoma proliferation and worsens the prognosis,” BMC Cancer, vol. 16, no. 1, p. 514, 2016.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Chun-hui Lai et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1135

Downloads

1484

Citations