Diagnostic Value Investigation and Bioinformatics Analysis of miR-31 in Patients with Lymph Node Metastasis of Colorectal Cancer

Zhang, Wu-wen; Ming, Xin-liang; Rong, Yuan; Huang, Chao-qun; Weng, Hong; Chen, Hao; Bian, Jun-mei; Wang, Fu-bing

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

Analytical Cellular Pathology

On this page

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

Research Article | Open Access

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

Diagnostic Value Investigation and Bioinformatics Analysis of miR-31 in Patients with Lymph Node Metastasis of Colorectal Cancer

Wu-wen Zhang,¹Xin-liang Ming,¹Yuan Rong,¹Chao-qun Huang,²Hong Weng,³Hao Chen,⁴Jun-mei Bian,⁵and Fu-bing Wang¹

Academic Editor: Silvia Cantara

Received14 Oct 2019

Revised18 Nov 2019

Accepted26 Nov 2019

Published16 Dec 2019

Abstract

Colorectal cancer (CRC) is one of the most frequent cancers occurring in developed countries. Distant CRC metastasis causes more than 90% of CRC-associated mortality. MicroRNAs (miRNAs) play a key role in regulating tumor metastasis and could be potential diagnostic biomarkers in CRC patients. This study is aimed at identifying miRNAs that can be used as diagnostic biomarkers for CRC metastasis. Towards this goal, we compared the expression of five miRNAs commonly associated with metastasis (i.e., miR-10b, miR-200c, miR-155, miR-21, and miR-31) between primary CRC (pCRC) tissues and corresponding metastatic lymph nodes (mCRC). Further, bioinformatics analysis of miR-31 was performed to predict target genes and related signaling pathways. Results showed that miR-31, miR-21, miR-10b, and miR-155 expression was increased to different extents, while miR-200c expression was lower in mCRC than that in pCRC. Moreover, we found that the level of both miR-31 and miR-21 was notably increased in pCRC when lymph node metastasis (LNM) was present, and the increase of miR-31 expression was more profound. Hence, upregulated miR-31 and miR-21 expression might be a miRNA signature in CRC metastasis. Moreover, we detected a higher miR-31 level in the plasma of CRC patients with LNM compared to patients without LNM or healthy individuals. With the bioinformatics analysis of miR-31, 121 putative target genes and transition of mitotic cell cycle and Wnt signaling pathway were identified to possibly play a role in CRC progression. We next identified seven hub genes via module analysis; of these, TNS1 was most likely to be the target of miR-31 and had significant prognostic value for CRC patients. In conclusion, miR-31 is significantly increased in the cancer tissues and plasma of CRC patients with LNM; thus, a high level of miR-31 in the plasma is a potential biomarker for the diagnosis of LNM of CRC.

1. Introduction

Colorectal cancer (CRC) is one of the most common malignancies and the fourth leading cause of cancer-related death worldwide [1]. However, although the development of early screening modalities and adjuvant chemotherapies has improved the prognosis of CRC, advanced CRC with lymph node metastasis (LNM) remains to have extremely poor prognosis. LNM is common among patients with advanced CRC, and given the poor prognosis, a more sensitive and specific method for LNM diagnosis can significantly benefit the therapeutic planning and clinical follow-up for CRC patients [2]. Although several gene products seem to contribute to the malignancy of CRC, accurate predictive factors of the prognosis and recurrence of CRC are yet to be identified.

Among the gene products and signaling molecules involved in the development of CRC, microRNAs (miRNAs), which are short noncoding RNAs measuring 18–25 nucleotides in length, might serve as molecular targets for both diagnosis and therapy. miRNAs have been demonstrated to contribute to a new mechanism of regulating gene expression and are involved in various biological processes of human cancers [3]. They can regulate gene expression posttranscriptionally, and bioinformatics analysis has suggested that miRNAs are capable of regulating the expression of numerous mammalian genes, among which tumor-promoting genes and tumor suppressor genes are included [4]. miRNAs were reported to be active in cancer development, acting as oncogenes (e.g., miR-155 and miR-21), tumor suppressors (e.g., miR-15a and miR-16-1), or metastasis promoters (e.g., miR-10b, miR-182, and miR-29a) [5]. Aberrant miRNA expression is associated with human cancers, and thus miRNA profiling, as one of the most modern modalities for molecular characterization of tumors, is used for cancer diagnosis and prognostic prediction [6, 7].

Studies have suggested miRNA profiles in distinguishing CRC tissues from normal colorectal mucosa. Luo et al. have observed that 164 miRNAs were aberrantly expressed in CRC [8]. Further, miR-31 and miR-20a were reported to be significantly elevated, whereas miR-145 and miR-143 were significantly downregulated in CRC tissues [9–11]. In addition, several studies have further described the correlations between dysregulated miRNA expression levels with tumor features. For example, the expressions of miR-21, miR-31, and miR-20a were positively correlated with the expression levels of histological markers (Ki-67 and CD34) in CRC tissues [12]. miR-145 was negatively correlated with the expression of K-ras gene, while miR-21 was positively correlated with the expression of K-ras gene in CRC [13]. Results of these studies indicated that dysregulated miRNAs were involved in cell proliferation and angiogenesis in CRC development. However, few studies have focused on the relationship between miRNAs and tumor stages. miRNAs might influence migration, invasion, and intravasation of tumor cells and act as metastasis promoters in breast cancer [14, 15], urothelial carcinomas, melanoma, and CRC [16]. Interestingly, upregulation of miR-21 and miR-10b correlates with LNM of breast cancer, while downregulation of miR-31, miR-335, and miR-126 is associated with tumor recurrence [16]. Eslamizadeh et al. have investigated the diagnostic value of miRNA profile in CRC. According to their study, the plasma levels of miR-21, miR-31, miR-20a, and miR-135b were rising with the higher stages of CRC. By contrast, miR-145, miR-let-7g, and miR-200c were decreasing with the higher stages of CRC. And the expression levels of plasma miR-21, miR-31, and miR-135b were significantly different between patients with stage II and III CRC [17].

miRNAs might promote or inhibit tumor metastasis by regulating the expression of target genes [18]. miR-31 was found to target NF-κB-inducing kinase (NIK) to negatively regulate the noncanonical NF-κB pathway, and loss of miR-31 therefore triggers oncogenic signaling in adult T cell leukemia [19]. Another study on head and neck carcinoma suggested that miR-31 could activate hypoxia-inducible factor to promote cancer progression by targeting factor-inhibiting hypoxia-inducible factor [20]. Besides, the tumor suppressor gene RhoBTB1 was also suggested to be regulated by miR-31, which was associated with the progression of human colon cancer [21]. However, conclusive studies on the target genes and pathways in CRC are relatively rare. To further explore the possible mechanism by which miR-31 regulates CRC, we performed bioinformatics analysis on miR-31 to predict target genes and signaling pathways. And a comprehensive approach was used to screen for the most possible target genes, including GO and KEGG enrichment analysis, PPI network construction, validation of expression levels, correlation analysis, and survival analysis of the target genes.

The aim of this study was to (1) identify miRNAs that can be biomarkers for CRC metastasis. Towards this goal, we compared the expression of five miRNAs commonly associated with metastasis (i.e., miR-10b, miR-200c, miR-155, miR-21, and miR-31) between primary CRC (pCRC) tissues and corresponding metastatic lymph nodes (mCRC). Further, we aimed to (2) determine the mechanism by which miR-31 regulates CRC metastasis and thus adopted bioinformatics analysis to explore putative target genes and related pathways involved in CRC. Finally, we sought to (3) identify the most possible key target genes of miR-31 in CRC using comprehensive bioinformatics methods.

2. Materials and Methods

2.1. Patients and Samples

All experimental procedures and pathological classification were approved by the Zhongnan Hospital Ethics Committee of Wuhan University and were performed following the International Union Against Cancer and American Joint Committee on Cancer TNM staging system for colon cancer established in 2003. Informed consent was obtained from each patient prior to sampling. Tumor tissues were collected from CRC patients after tumor resection at Zhongnan Hospital of Wuhan University. None of them received any preoperative treatment. We designed three sets of experiments. In set 1, 9 primary CRC (pCRC) tissue samples and corresponding mCRC of the same patients (stage III/IV CRC) were paired to identify metastasis-related miRNAs. Each lymph node was tested via hematoxylin and eosin staining to confirm the LNM. In set 2, pCRC tissue samples were collected from matched CRC patients with or without LNM to validate the candidate miRNAs identified in set 1. In set 3, before surgery or preoperative treatments, blood plasmas were collected from 28 CRC patients (stage III/IV) with LNM, 28 patients (stage I/II) without LNM, and 28 age- and sex-matched healthy volunteers. All tissues and plasmas were well preserved immediately after collection until RNA extraction.

2.2. RNA Extraction

Small RNAs from the plasma were extracted with the miRVana PARIS RNA isolation kit (Qiagen, Germany). Briefly, 250 μl of plasma was thawed on ice followed by centrifugation at 14,000 rpm for 5 minutes to remove cell debris and organelles. We then lysed 150 μl of supernatant with an equal volume of solution (Qiagen, Germany). To normalize the variations among the samples, 25 fmol of synthetic Caenorhabditis elegans miRNA cel-miR-39 (BioVendor, Europe) was added into each denatured sample during the RNA extraction [22]. Small RNAs were then purified following the manufacturer’s protocol except that 45 μl of nuclease-free water was used to elute small RNAs.

Total RNAs from the tumor tissues were extracted using TRIzol reagent (TaKaRa, Tokyo, Japan) following the manufacturer’s protocol. The concentration and purity of total RNAs were measured using the Smart Spec Plus spectrophotometer (Thermo Scientific Inc., USA).

2.3. Real-Time PCR and Expression Analysis

Plasma RNA or 100 ng tissue RNA was polyadenylated and reverse transcribed into cDNA using the miScript Reverse Transcription kit (Qiagen, Germany). RT-PCR for each sample was performed in duplicates using the miScript SYBR Green PCR kit (Qiagen, Hilden, Germany) on a DNA Engine Opticon II system (Bio-Rad Laboratories, Inc., Hercules, CA, USA). The miRNA-specific primers were designed based on the miRNA sequences obtained from the miRBase database. Each amplification reaction was performed in a final volume of 20 μl containing 1 μl of cDNA, 0.25 mM of each primer, and Green PCR Master mix. At the end of PCR, the melting curve analysis and electrophoresis of the PCR products on 3.0% agarose gels were performed to determine the specificity of amplification. U6 snRNA was used as the internal control.

The quantification of mature miRNA was performed using BioRad CFX Manager (Bio-Rad Laboratories, Inc., Hercules, CA, USA). The cycle threshold (Ct) was defined as the cycle number when the fluorescent signal crossed the threshold in PCR.

The relative miRNA expression was normalized to U6 levels and calculated through the 2^-ΔΔCt method [23]. The data were presented as the fold changes of expression relative to normal tissues.

2.4. Identification of Metastasis-Related miRNAs in pCRC and mCRC

To identify miRNAs specific for metastatic CRC, miRNAs in nine matched pCRC tissues and corresponding mCRC were analyzed via RT-PCR. We first tested the expression of five miRNAs, i.e., miR-10b, miR-200c, miR-155, miR-21, and miR-31, because they have been suggested to closely associate with metastasis in multiple tumor types.

2.5. Prediction of Putative Target Genes of miR-31

Target genes of miR-31 were predicted using three online databases frequently used for miRNA target prediction: TargetScan Release 7.2, miRDB, and DIANA-microT web server v5.0 [24–26]. The overlapping parts of the possible target genes were considered as the putative target genes.

2.6. GO and KEGG Clustering Analysis of Putative Target Genes

Putative target genes were uploaded into the function annotation portal of The Database for Annotation, Visualization, and Integrated Discovery (DAVID), online bioinformatics resources for investigating the biological meaning of large gene lists [27]. Gene ontology (GO) analysis was utilized to investigate the possible roles of target genes involved in biological processes (BP), cellular components (CC), and molecular functions (MFs). Kyoto Encyclopedia of Genes and Genomes (KEGG) clustering analysis were adopted to map genes to a related pathway. A value of <0.05 and a false discovery rate of <0.05 were applied to identify significant GO and KEGG items.

2.7. Construction of PPI Network and Module Analysis

To study the interaction network among functional proteins, we constructed the protein-protein interactive (PPI) network with the StringApp plugin in Cytoscape from the STRING database [28]. Cytoscape 3.51 is a public platform for establishing biomolecular interaction networks [29]. Further, the CytoHubba containing 12 predicted algorithms (Stress, Betweenness, DMNC, Degree, MNC, MCC, BottleNeck, EPC, Closeness, Radiality, EcCentricity, and ClusteringCoefficient) was used to identify hub genes that were most likely key target genes of miR-31.

2.8. Validation of the Expression of miR-31 and Hub Genes in Colorectal Cancer

To validate the expression of miR-31 and hub genes in CRC, we downloaded RNA expression data of GDC TCGA Colon Cancer and Rectal Cancer and controls from the UCSC Xena platform [30]. All the expression data have been normalized with log2 transformation. Then, the RNA counts of miR-31 and hub genes were extracted. Expression levels of miR-31 and hub genes were compared between colon cancer, rectal cancer, and controls with Student’s -test using GraphPad Prism 6.0 (GraphPad Software, Inc., La Jolla, CA, USA). The difference in miR-31 and hub gene expression was shown as a box plot. In addition, to reveal the relationship between miR-31 and the hub genes, Spearman’s correlation analysis and linear regression plots construction were performed with GraphPad Prism 6.0.

2.9. Prognostic Value Evaluation of Hub Genes

To evaluate the value of these hub genes in CRC patients’ prognosis, we used the GTEx survival data of colon and rectal cancer patients with the GEPIA online tool [31]. Overall survival (OS) and disease-free survival (DFS) of colon and rectal cancer patients were estimated by plotting Kaplan-Meier survival curves. Hazard ratio (HR) was calculated with the log-rank test to assess the effects of hub gene expression on patients’ survival.

2.10. Statistical Analysis

All data were expressed as the . Student’s -test and Mann-Whitney test were applied to assess the statistical significance. To determine the extent to which the obtained -ΔCt efficiently distinguishes different clinical subsets, the receiver operating characteristic (ROC) analysis was performed, and the area under curve (AUC) was used as an indicator of the distinguished capability. All statistical analyses were performed using SPSS 13.0 software, and values below 0.05 were defined as statistically significant.

3. Results

3.1. Metastasis-Related miRNAs in pCRC and mCRC

The expressions of miR-10b, miR-155, miR-31, and miR-21 were higher, and the expression of miR-200c was lower in mCRC than that in corresponding pCRC tissues. These changes were observed in all nine cases tested. In particular, the expression of miR-21 and miR-31 in mCRC was both more than two-folds higher than those in pCRC tissues (Table 1). Therefore, we focused on miR-21 and miR-31 in further investigations to disclose their significance.

3.2. miR-21 and miR-31 Are Upregulated in CRC with LNM

To explore the value of miR-21 and miR-31 in the diagnosis of CRC metastasis, we assessed their expression in the pCRC tissues collected from patients with or without LNM via RT-PCR. Each group comprised 33 patients. The miR-21 level was significantly higher by 1.62-fold in average in pCRC tissues when LNM was present, as compared with pCRC tissues without LNM () (Figure 1). Similarly, the miR-31 level was also considerably higher in the pCRC tissues of patients with LNM (average, 2.36-fold, ).

3.3. Plasma miR-31 Distinguishes CRC Metastasis

To assess whether the upregulation of miR-31 can identify CRC or, more specifically, CRC metastasis, we determined the miR-31 level in the plasma of CRC patients and healthy donors. In total, 84 plasma samples were obtained from 28 stage I/II patients without LNM, 28 stage III/IV patients with LNM, and 28 healthy donors. The -ΔCt value was used to indicate miR-31 expression in the plasma. Compared with that of healthy donors, CRC patients had significantly elevated miR-31 in the plasma regardless of LNM (). Furthermore, patients with LNM had even higher plasma miR-31 expression than patients without LNM () (Figure 2(a)). The ROC curve revealed that miR-31 was potentially a valuable biomarker for discriminating CRC patients with LNM from CRC patients without LNM, as indicated by an AUC of 0.89 (95% CI: 0.81–0.97, ) (Figure 2(b)). In the ROC assay, a -ΔCt value of -8.6 (normalized) in patients with LNM was identified as a cut-off to discriminate metastatic CRC from nonmetastatic CRC. The optimal specificity and sensitivity were 86.2% and 78.5%, respectively (Figure 2(b)).

(a)

(b)

3.4. Putative Target Genes of miR-31

There were 477, 613, and 595 target genes of miR-31 predicted using the TargetScan, miRDB, and DIANA-microT web server, respectively. The 121 overlapping genes of the 3 databases were considered target genes of miR-31 (Figure 3).

3.5. GO and KEGG Enrichment Analysis

There were four, three, and six terms significantly clustered for BP, CC, and MF through GO analysis, respectively (Figure 4). The top three GO terms were cerebral cortex development, transition of mitotic cell cycle, and locomotion involved in locomotory behavior in BP; cytoplasm, focal adhesion, and cAMP-dependent protein kinase complex in CC; and protein heterodimerization activity, SH3 domain binding, and ion channel binding in MF. Seven pathways were identified via KEGG enrichment analysis (Figure 4), and the top three pathways were melanogenesis, ubiquitin-mediated proteolysis, and Wnt signaling pathway.

3.6. PPI Network and Hub Genes

The PPI network was constructed with the 121 putative target genes using StringApp (Figure 5(a)). Each of the 121 target genes may associate with others, which may constitute interactive modules involved in the progression of CRC. After module analysis using 12 algorithms, the results were sorted in a descending order, and seven hub genes (ELAVL1, PPP3CA, DICER1, CBL, GNA13, SSH1, and TNS1) were identified, which are probably the key target genes of miR-31 (Figure 5(b)).

(a)

(b)

3.7. Validation of the Expression of miR-31 and Hub Genes

To compare the expression of miR-31 between CRC tissues and controls, miRNA expression data of 615 tumor tissues and 11controls was downloaded. The miR-31 level was significantly higher in CRC tissues than that in controls (Figure 6(a)). As for the hub genes, mRNA expression data of 638 GDC TCGA tumor tissues (including colon cancer and rectal cancer) and 51 controls were downloaded (Figure 6(a)). Except for ELAVL1, all of the other six hub genes were significantly downregulated in colon and rectal tumor tissues compared to those in controls. Considering that miR-31 was significantly upregulated in CRC, these six hub genes are possibly regulated by miR-31 in CRC. Further, correlation analysis showed that TNS1 level was negatively correlated with miR-31 (Figure 6(b)); meanwhile, no statistically significant relationships between the other six hub genes and miR-31 were found.

(a)

(b)

Figure 6

Validation of the expression of miR-31 and hub genes in CRC and controls tissues. (a) Expression of miR-31 and hub gene were detected in 615 CRC tissues and 11 control tissues. miR-31 was significantly upregulated in CRC tissues compared to control tissues. Expression of the hub genes was detected in 638 CRC tissues and 51 control tissues. Six hub genes—PPP3CA, DICER1, CBL, GNA13, SSH1, and TNS1—were significantly downregulated in CRC tissues compared to control tissues. No significant difference of ELAVL1 expression was observed between CRC tissues and control tissues. (b) Correlation analysis showed that TNS1 expression level was negatively correlated with miR-31 expression level in CRC and control tissues.

3.8. Prognostic Value Evaluation of Hub Genes

Data of 362 CRC patients were used in Kaplan-Meier survival analysis for seven hub genes. The results indicated that low TNS1 expression was significantly associated with improved OS (, ) and DFS (, ) (Figure 7). However, for the other six hub genes including ELAVL1, PPP3CA, DICER1, CBL, GNA13, and SSH1, no statistical significance was found in the survival analysis (Supplement Figures 1–6). Further studies with larger cohorts need to be carried out.

(a)

(b)

4. Discussion

Amplification, deletion, and rearrangement of miRNAs are frequently present in human cancers. Some altered miRNA expression can promote tumorigenesis, and some miRNAs act as tumor suppressors [32]. Recently, the term “metastamiR” was proposed to describe miRNAs that are associated with tumor metastasis. MetastamiRs can be prometastatic or antimetastatic [16]. For example, miR-10b was firstly reported to contribute to breast cancer metastasis [33]. Later, it was reported that miR-335 could inhibit the invasion of metastatic breast cancer cell. Increasing evidence has demonstrated altered miRNA expression between primary and metastatic tumors, implying the important role of miRNAs in tumor metastasis. Previous studies have reported that 5 upregulated miRNAs and 14 downregulated miRNAs are involved in the metastasis of CRC, but the association between miRNA levels and CRC metastasis remains unclear [34, 35]. Thus, it is necessary to investigate the differential miRNA levels between primary tumors and metastatic tissues. A matched comparison method is optimal because it can exclude endogenous differences of miRNA expression. In our research, the levels of five metastasis-related miRNAs were detected and most of them were increased in mCRC. The trend of miRNA expression in our results was similar to previous investigations [33, 36, 37]. In our study, the higher expression of miR-31, miR-21, miR-10b, and miR-155 in mCRC indicated a miRNA signature might predict CRC metastasis. Because miR-31 and miR-21 were the most elevated miRNAs among the five selected miRNAs, we then wondered if they can be potential metastatic biomarkers of pCRC. As expected, their expression was significantly elevated in the pCRC with LNM. Interestingly, the increase of miR-31 was more profound, suggesting that miR-31 might be a more sensitive biomarker to predict CRC metastasis.

Although the alterations in miRNA expression in pCRC with LNM might be helpful for diagnosis, it is difficult for clinical practitioners to collect CRC tissues from patients. A more convenient and less invasive detection approach, such as blood testing, will be substantially beneficial for the prediction or diagnosis of CRC metastasis. Nucleic acid levels in the circulation can be used for the diagnosis of CRC [7, 38]. Previously, plasma miRNAs levels were reported to be highly correlated with miRNA expression in tumor tissue from breast cancer patients [39]. miRNAs have been detected in the serum and plasma of CRC, ovarian cancer, and prostate cancer patients. The plasma miRNAs are more stable and consistent than other circulating nucleic acids. Hence, they could be optimal biomarkers for cancer diagnosis. For example, increased plasma miR-92 levels can accurately discriminate CRC from gastric cancer and benign disease [40]. Plasma miR-141 has been proposed in diagnosing metastatic colon cancer [41]. We also evaluated the plasma miR-31 levels in CRC patients with or without LNM. Compared to CRC patients without LNM, patients with LNM have significantly higher plasma miR-31 level. miR-31 yielded a ROC curve area of 0.89 with a sensitivity of 78.5% and specificity of 86.2% in distinguishing CRC with LNM from CRC without LNM, using a cut-off value of -8.6 (normalized).

As far as we know, plasma miR-31 has been investigated as a biomarker for oral cancer [42]. Eslamizadeh et al. have reported the plasma miR-31 level was rising with the higher stages of CRC [17]. Here, we come to the conclusion that plasma miR-31 is of diagnostic value for CRC with LNM. However, to further validate our study, future investigations should recruit a large sample cohort.

To explore the associations of miR-31 with CRC metastasis, we predicted putative target genes of miR-31 and then carried out enrichment analysis of functions and signaling pathways of the target genes with bioinformatics methods. Results of enrichment analysis show that the target genes of miR-31 are possibly involved in some cancer-related biological process and signaling pathways such as G1/S transition of mitotic cell cycle, ATP binding, and Wnt signaling pathway. PPI network construction and module analysis were further conducted, and seven hub genes were identified, which were more likely to be target genes of miR-31. Through validation of expression and correlation analysis, among the seven hub genes, the TNS1 level was found lower in CRC tissues compared to controls and was negatively correlated with the miR-31 level. Thus, TNS1 is most possibly to be regulated by miR-31. Through survival analysis, we found that the TNS1 level was significantly associated with OS and DFS of CRC patients. In addition, low expression of TNS1 was predictive of improved OS and DFS for CRC patients.

TNS1 is a 220 kD protein localized to focal adhesions and regions of the plasma membrane where the cell attaches to the extracellular matrix. TNS1 protein plays a role in regulating cell motility and is suggested to be involved in tumorigenesis [43]. Elevated TNS1 levels were associated with a poor overall survival in CRC patients. Therefore, we suspect miR-31 may target TNS1, contributing to improved outcomes for CRC patients. Further experimental studies need to be performed to validate whether TNS1 is actually targeted by miR-31 in CRC.

5. Conclusion

In summary, miR-31 is significantly elevated in tumor tissues and plasma of CRC patients with LNM. Plasma miR-31 may be utilized as a biomarker for CRC with LNM. In addition, elevated miR-31 may contribute to improved outcomes for CRC patients by targeting TNS1.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

No potential conflicts of interest were disclosed.

Authors’ Contributions

Wu-wen Zhang, Xin-liang Ming, Yuan Rong, Fu-bing Wang, and Jun-mei Bian designed and performed the experiments and performed the data analysis and presentation. All authors prepared and revised the manuscript. Wu-wen Zhang, Xin-liang Ming, and Yuan Rong contributed equally to this work.

Acknowledgments

The authors thank the members of their laboratory and their collaborators for their research work. This work was supported by the Sixth Training Program Foundation for Wuhan Medical Talents by Health Commission of Wuhan. This work was also supported by the Science and Technology Innovation Fostering Foundation of Zhongnan Hospital of Wuhan University (grant No. znpy2018027).

Supplementary Materials

Supplement Figure 1: Kaplan-Meier survival curves of 362 CRC cases according to ELAVL1 expression levels. Supplement Figure 2: Kaplan-Meier survival curves of 362 CRC cases according to PPP3CA expression levels. Supplement Figure 3: Kaplan-Meier survival curves of 361 CRC cases according to DICER1 expression levels. Supplement Figure 4: Kaplan-Meier survival curves of 360 CRC cases according to CBL expression levels. Supplement Figure 5: Kaplan-Meier survival curves of 362 CRC cases according to GNA13 expression levels. Supplement Figure 6: Kaplan-Meier survival curves of 362 CRC cases according to SSH1 expression levels. (Supplementary Materials)

References

M. Arnold, M. S. Sierra, M. Laversanne, I. Soerjomataram, A. Jemal, and F. Bray, “Global patterns and trends in colorectal cancer incidence and mortality,” Gut, vol. 66, no. 4, pp. 683–691, 2017.
View at: Publisher Site | Google Scholar
I. D. Nagtegaal and H. J. Schmoll, “Colorectal cancer: what is the role of lymph node metastases in the progression of colorectal cancer?” Nature Reviews Gastroenterology & Hepatology, vol. 14, no. 11, pp. 633-634, 2017.
View at: Publisher Site | Google Scholar
A. A. Svoronos, D. M. Engelman, and F. J. Slack, “OncomiR or tumor suppressor? The duplicity of microRNAs in cancer,” Cancer Research, vol. 76, no. 13, pp. 3666–3670, 2016.
View at: Publisher Site | Google Scholar
T. Schepeler, J. T. Reinert, M. S. Ostenfeld et al., “Diagnostic and prognostic microRNAs in stage II colon cancer,” Cancer Research, vol. 68, no. 15, pp. 6416–6424, 2008.
View at: Publisher Site | Google Scholar
J. M. Bouyssou, S. Manier, D. Huynh, S. Issa, A. M. Roccaro, and I. M. Ghobrial, “Regulation of microRNAs in cancer metastasis,” Biochimica et Biophysica Acta, vol. 1845, no. 2, pp. 255–265, 2014.
View at: Publisher Site | Google Scholar
K. Hur, Y. Toiyama, Y. Okugawa et al., “Circulating microRNA-203 predicts prognosis and metastasis in human colorectal cancer,” Gut, vol. 66, no. 4, pp. 654–665, 2017.
View at: Publisher Site | Google Scholar
H. Imaoka, Y. Toiyama, H. Fujikawa et al., “Circulating microRNA-1290 as a novel diagnostic and prognostic biomarker in human colorectal cancer,” Annals of Oncology, vol. 27, no. 10, pp. 1879–1886, 2016.
View at: Publisher Site | Google Scholar
X. Luo, B. Burwinkel, S. Tao, and H. Brenner, “MicroRNA signatures: novel biomarker for colorectal cancer?” Cancer Epidemiology, Biomarkers & Prevention, vol. 20, no. 7, pp. 1272–1286, 2011.
View at: Publisher Site | Google Scholar
V. Kulda, M. Pesta, O. Topolcan et al., “Relevance of miR-21 and miR-143 expression in tissue samples of colorectal carcinoma and its liver metastases,” Cancer Genetics and Cytogenetics, vol. 200, no. 2, pp. 154–160, 2010.
View at: Publisher Site | Google Scholar
K. Motoyama, H. Inoue, Y. Takatsuno et al., “Over- and under-expressed microRNAs in human colorectal cancer,” International Journal of Oncology, vol. 34, no. 4, pp. 1069–1075, 2009.
View at: Publisher Site | Google Scholar
O. Slaby, M. Svoboda, P. Fabian et al., “Altered expression of miR-21, miR-31, miR-143 and miR-145 is related to clinicopathologic features of colorectal cancer,” Oncology, vol. 72, no. 5-6, pp. 397–402, 2007.
View at: Publisher Site | Google Scholar
S. S. Emami, A. Akbari, A.-A. Zare et al., “MicroRNA expression levels and histopathological features of colorectal cancer,” Journal of Gastrointestinal Cancer, vol. 50, no. 2, pp. 276–284, 2019.
View at: Publisher Site | Google Scholar
F. Moghadamnia, P. Ghoraeian, S. Minaeian, A. Talebi, F. Farsi, and A. Akbari, “MicroRNA expression and correlation with mRNA levels of colorectal cancer-related gene,” Journal of Gastrointestinal Cancer, 2019.
View at: Publisher Site | Google Scholar
M. Y. Fong, W. Zhou, L. Liu et al., “Breast-cancer-secreted miR-122 reprograms glucose metabolism in premetastatic niche to promote metastasis,” Nature Cell Biology, vol. 17, no. 2, pp. 183–194, 2015.
View at: Publisher Site | Google Scholar
J. N. Goh, S. Y. Loo, A. Datta et al., “microRNAs in breast cancer: regulatory roles governing the hallmarks of cancer,” Biological Reviews, vol. 91, no. 2, pp. 409–428, 2016.
View at: Publisher Site | Google Scholar
N. M. A. White, E. Fatoohi, M. Metias, K. Jung, C. Stephan, and G. M. Yousef, “Metastamirs: a stepping stone towards improved cancer management,” Nature Reviews Clinical Oncology, vol. 8, no. 2, pp. 75–84, 2011.
View at: Publisher Site | Google Scholar
S. Eslamizadeh, M. Heidari, S. Agah et al., “The role of microRNA signature as diagnostic biomarkers in different clinical stages of colorectal cancer,” Cell Journal, vol. 20, no. 2, pp. 220–230, 2018.
View at: Publisher Site | Google Scholar
H. Zhang, Y. Li, and M. Lai, “The microRNA network and tumor metastasis,” Oncogene, vol. 29, no. 7, pp. 937–948, 2010.
View at: Publisher Site | Google Scholar
M. Yamagishi, K. Nakano, A. Miyake et al., “Polycomb-mediated loss of miR-31 activates NIK-dependent NF-κB pathway in adult T cell leukemia and other cancers,” Cancer Cell, vol. 21, no. 1, pp. 121–135, 2012.
View at: Publisher Site | Google Scholar
C. J. Liu, M. M. Tsai, P. S. Hung et al., “miR-31 ablates expression of the HIF regulatory factor FIH to activate the HIF pathway in head and neck carcinoma,” Cancer Research, vol. 70, no. 4, pp. 1635–1644, 2010.
View at: Publisher Site | Google Scholar
R. S. Xu, X. D. Wu, S. Q. Zhang et al., “The tumor suppressor gene RhoBTB1 is a novel target of miR-31 in human colon cancer,” International Journal of Oncology, vol. 42, no. 2, pp. 676–682, 2013.
View at: Publisher Site | Google Scholar
E. M. Kroh, R. K. Parkin, P. S. Mitchell, and M. Tewari, “Analysis of circulating microRNA biomarkers in plasma and serum using quantitative reverse transcription-PCR (qRT-PCR),” Methods, vol. 50, no. 4, pp. 298–301, 2010.
View at: Publisher Site | Google Scholar
K. J. Livak and T. D. Schmittgen, “Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method,” Methods, vol. 25, no. 4, pp. 402–408, 2001.
View at: Publisher Site | Google Scholar
V. Agarwal, G. W. Bell, J. W. Nam, and D. P. Bartel, “Predicting effective microRNA target sites in mammalian mRNAs,” eLife, vol. 4, 2015.
View at: Publisher Site | Google Scholar
W. Liu and X. Wang, “Prediction of functional microRNA targets by integrative modeling of microRNA binding and target expression data,” Genome Biology, vol. 20, no. 1, p. 18, 2019.
View at: Publisher Site | Google Scholar
M. Reczko, M. Maragkakis, P. Alexiou, I. Grosse, and A. G. Hatzigeorgiou, “Functional microRNA targets in protein coding sequences,” Bioinformatics, vol. 28, no. 6, pp. 771–776, 2012.
View at: Publisher Site | Google Scholar
D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.
View at: Publisher Site | Google Scholar
S. Mishra, T. Yadav, and V. Rani, “Exploring miRNA based approaches in cancer diagnostics and therapeutics,” Critical Reviews in Oncology/Hematology, vol. 98, pp. 12–23, 2016.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
M. Goldman, B. Craft, M. Hastie et al., “The UCSC Xena platform for public and private cancer genomics data visualization and interpretation,” bioRxiv, p. 326470, 2019.
View at: Publisher Site | Google Scholar
Z. Tang, C. Li, B. Kang, G. Gao, C. Li, and Z. Zhang, “GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses,” Nucleic Acids Research, vol. 45, no. W1, pp. W98–w102, 2017.
View at: Publisher Site | Google Scholar
A. J. Schetter, S. Y. Leung, J. J. Sohn et al., “MicroRNA expression profiles associated with prognosis and therapeutic outcome in colon adenocarcinoma,” JAMA, vol. 299, no. 4, pp. 425–436, 2008.
View at: Publisher Site | Google Scholar
L. Ma, J. Teruya-Feldstein, and R. A. Weinberg, “Tumour invasion and metastasis initiated by microRNA-10b in breast cancer,” Nature, vol. 449, no. 7163, pp. 682–688, 2007.
View at: Publisher Site | Google Scholar
X. F. Le, O. Merchant, R. C. Bast, and G. A. Calin, “The roles of microRNAs in the cancer invasion-metastasis cascade,” Cancer Microenvironment, vol. 3, no. 1, pp. 137–147, 2010.
View at: Publisher Site | Google Scholar
I. de Krijger, L. J. Mekenkamp, C. J. Punt, and I. D. Nagtegaal, “MicroRNAs in colorectal cancer metastasis,” The Journal of Pathology, vol. 224, no. 4, pp. 438–447, 2011.
View at: Publisher Site | Google Scholar
B. Qian, D. Katsaros, L. Lu et al., “High miR-21 expression in breast cancer associated with poor disease-free survival in early stage disease and high TGF-β1,” Breast Cancer Research and Treatment, vol. 117, no. 1, pp. 131–140, 2009.
View at: Publisher Site | Google Scholar
Y. Dong, W. K. Wu, C. W. Wu, J. J. Sung, J. Yu, and S. S. Ng, “MicroRNA dysregulation in colorectal cancer: a clinical perspective,” British Journal of Cancer, vol. 104, no. 6, pp. 893–898, 2011.
View at: Publisher Site | Google Scholar
S. Wang, J. Xiang, Z. Li et al., “A plasma microRNA panel for early detection of colorectal cancer,” International Journal of Cancer, vol. 136, no. 1, pp. 152–161, 2015.
View at: Publisher Site | Google Scholar
F. Wang, Z. Zheng, J. Guo, and X. Ding, “Correlation and quantitation of microRNA aberrant expression in tissues and sera from patients with breast tumor,” Gynecologic Oncology, vol. 119, no. 3, pp. 586–593, 2010.
View at: Publisher Site | Google Scholar
E. K. Ng, W. W. Chong, H. Jin et al., “Differential expression of microRNAs in plasma of patients with colorectal cancer: a potential marker for colorectal cancer screening,” Gut, vol. 58, no. 10, pp. 1375–1381, 2009.
View at: Publisher Site | Google Scholar
P. S. Mitchell, R. K. Parkin, E. M. Kroh et al., “Circulating microRNAs as stable blood-based markers for cancer detection,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 30, pp. 10513–10518, 2008.
View at: Publisher Site | Google Scholar
C. J. Liu, S. Y. Kao, H. F. Tu, M. M. Tsai, K. W. Chang, and S. C. Lin, “Increase of microRNA miR-31 level in plasma could be a potential marker of oral cancer,” Oral Diseases, vol. 16, no. 4, pp. 360–364, 2010.
View at: Publisher Site | Google Scholar
E. H. Hall, A. E. Daugherty, C. K. Choi, A. F. Horwitz, and D. L. Brautigan, “Tensin1 requires protein phosphatase-1alpha in addition to RhoGAP DLC-1 to control cell polarization, migration, and invasion,” The Journal of Biological Chemistry, vol. 284, no. 50, pp. 34713–34722, 2009.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Wu-wen Zhang 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

737

Downloads

1195

Citations