A Circulating miRNA-Based Scoring System Established by WGCNA to Predict Colon Cancer

Qin, Da; Wei, Rui; Liu, Si; Zhu, Shengtao; Zhang, Shutian; Min, Li

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

Analytical Cellular Pathology

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

Tumor Microenvironment and Early Diagnosis of Cancers

View this Special Issue

Research Article | Open Access

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

A Circulating miRNA-Based Scoring System Established by WGCNA to Predict Colon Cancer

Da Qin,¹Rui Wei,¹Si Liu,¹Shengtao Zhu,¹Shutian Zhang,¹and Li Min¹

Academic Editor: Silvia Cantara

Received29 Aug 2019

Revised24 Oct 2019

Accepted29 Oct 2019

Published01 Dec 2019

Abstract

Introduction. Circulation microRNAs (miRNAs) perform as potential diagnostic biomarkers of many kinds of cancers. This study is aimed at identifying circulation miRNAs as diagnostic biomarkers in colon cancer. Methods. We conducted a weighted gene coexpression network analysis (WGCNA) in miRNAs to find out the expression pattern among circulation miRNAs by using a “WGCNA” package in R. Correlation analysis was performed to find cancer-related modules. Differentially expressed miRNAs (DEmiRs) in colon cancer were identified by a “limma” package in R. Hub gene analysis was conducted for these DEmiRs in the cancer-related modules by the “closeness” method in cytoscape software. Then, logistic regression was performed to identify the independent risk factors, and a scoring system was constructed based on these independent risk factors. Then, we use data from the GEO database to confirm the reliability of this scoring system. Results. A total of 9 independent coexpression modules were constructed based on the expression levels of 848 miRNAs by WGCNA. After correlation analysis, green (, ) and yellow (, ) modules were strongly correlated with cancer development. 20 hub genes were found after hub gene analysis in these DEmiRs by cytoscape. Among all these hub genes, hsa-miR-23a-3p (, ) and hsa-miR-663a (, ) were identified as an independent risk factor of colon cancer by multivariate regression. Furthermore, a scoring system was built to predict the probability of colon cancer based on both of these miRNAs, the area under the curve (AUC) of which was 0.828. Data from GSE106817 and GSE112264 was used to confirm this scoring system. And the AUC of them was 0.980 and 0.917, respectively. Conclusion. We built a scoring system based on circulation hub miRNAs found by WGCNA to predict the development of colon cancer.

1. Introduction

The incidence of colon cancer continues to grow these years. According to large-scale epidemiological studies [1, 2], colon cancer has become a major disease threatening patients’ health. Researchers predicted that, in 2019, 145,600 new colon cancer would be diagnosed and 51,020 persons will die because of this disease [1]. Colon cancer develops insidiously, which means that it will not easily exhibit specific symptoms until turning into extremely terminal stages. This will cause irreparable effect to patients. Screening methods which could find colon cancer patients from population were important.

Although guidelines [3, 4] emphasized the importance of noninvasive tools to screen colon cancer, screening of colon cancer still relies on colonoscopies for high-risk population, which would be difficult to popularize in China because of its pain and poor compliance of patients. Moreover, traditional tumor biomarkers ([5], p., 1), such as CEA, were characterized by poor sensitivity and specificity, which may lead to missed diagnosis and waste of medical resources. The establishment of a sensitive and specific serological screening method is particularly important to improve the detection rate of colon cancer.

Previous studies [6–11] have revealed that many kinds of miRNAs play crucial roles in tumorigenesis. miRNAs could be found in many kinds of biofluids, such as blood and urine [12]. So, there would be a potential way to develop circulating miRNAs as minimally invasive biomarkers of cancer.

WGCNA is a method to analyze gene expression patterns in multiple samples. Genes with similar expression patterns would be clustered, and correlation between modules and specific traits or phenotypes could be calculated. Therefore, WGCNA [13] is widely used in the study of cancer researches. In this study, WGCNA-based methods were used to identify gene modules correlated with the development of colon cancer, and the most representative miRNAs were identified as circulation markers to help in screening colon cancer.

2. Materials and Methods

2.1. Weighted Gene Coexpression Network Analysis

We conducted a weighted gene coexpression network analysis in miRNAs to find out the expression pattern among circulation miRNAs by using a “WGCNA” package [14] in R. Data was from GSE61741 [15] in Gene Expression Omnibus (GEO) database. miRNAs with a similar expression pattern were regarded as a specific module and marked by a unique color. Correlations between these modules were calculated, and heat map was performed to show the independence between each of these modules. Then, correlation analysis was performed to find cancer-related modules. Modules which are strongly correlated with cancer development were chosen to undergo further analysis.

2.2. Hub Gene Selection

Differentially expressed miRNAs in colon cancer were identified by a “limma” [16] package in R. Hub gene analysis was conducted for these DEmiRs in the cancer-related modules by the closeness method in cytoscape [17] software. miRNA-miRNA interaction graph was performed to visualize the relationship between these circulation miRNAs.

2.3. Scoring System Establishing

Logistic regression was performed to identify the independent risk factor, and the scoring system was constructed based on the independent risk factors. Furthermore, we visualized this scoring system by nomogram using an “rms” package in R. A receiver operating characteristic (ROC) curve was performed to evaluate the reliability of this scoring system using a “pROC” package [18] in R.

2.4. External Validation

To confirm the reliability of this scoring system, data from GSE106817 [19] and GSE112264 [20] database was used. Violin graphs were performed to describe gene expression levels. ROC curves were graphed to evaluate the reliability of this scoring system using a “pROC” package [18] in R.

2.5. Statistical Analysis

All statistical analyses were performed using R software (version 3.5.2; https://www.r-project.org/). WGCNA were conducted by a “WGCNA” package (version 1.68). Differentially expressed miRNAs were identified by a “limma” package (version 3.38.3). Then, hub gene analysis was conducted by the closeness method in cytoscape software (version 3.7.1). And an “rms” package (version 5.1-3.1) was used to visualize the scoring system. Then, ROC curves were graphed using a “pROC” package (version 1.15.3). Statistical significance was defined as .

3. Results

3.1. Weighted Gene Coexpression Network Analysis

123 samples from GSE61741 were clustered based on their gene expression patterns (Supplementary Figure A). Then, a total of 9 coexpression modules were constructed based on the expression levels of 848 miRNAs by WGCNA (Figure 1(a)). Heat maps were performed to prove the independence between each gene (Supplementary Figure B). After correlation analysis, green (, ) and yellow (, ) modules were strongly correlated with cancer development (Figure 1(b), Supplementary Figure C). Scatter diagrams were graphed to prove the correlation between each gene in modules and cancer development (Figures 1(c) and 1(d)).

(a)

(b)

(c)

(d)

3.2. Hub Gene Selection

DEmiRs in these modules were found by a “limma” package, and after hub gene analysis in these DEmiRs by closeness software in cytoscape, 20 hub genes were identified (Table 1). The miRNA-miRNA interaction network is shown in Figure 2(a) and Supplementary Figure D.

(a)

(b)

(c)

(d)

(e)

3.3. Scoring System Establishing

hsa-miR-23a-3p (Figure 2(b)) and hsa-miR-663a (Figure 2(c)) levels were significantly increased in the blood of colon cancer, values of which were <0.0001 and 0.0388, respectively. Then, hsa-miR-23a-3p (, ) and hsa-miR-663a (, ) were identified as an independent risk factor of colon cancer by multivariate regression. Furthermore, a scoring system was built to predict the probability of colon cancer based on both of these miRNAs. The area under the curve (AUC) of this scoring system was 0.828 (Figure 2(d)). A nomogram was performed to visualize the scoring system (Figure 2(e)).

3.4. External Validation

Data from GSE112264 and GSE106817 was used to confirm this scoring system.

For GSE112264, a total of 41 blood samples of healthy control and 50 blood samples of colon cancer were analyzed. hsa-miR-23a-3p (Figure 3(a)) and hsa-miR-663a (Figure 3(b)) levels were significantly increased in the blood of colon cancer, values of which were both <0.0001. The receiver operating curve is shown as Figure 3(c), AUC of which was 0.917. When choosing a threshold of 0.7, the sensitivity of this scoring system was 0.97 and the specificity was 0.82.

(a)

(b)

(c)

(d)

(e)

(f)

For GSE106817, a total of 115 blood samples of healthy control and 115 blood samples of colon cancer were analyzed. hsa-miR-23a-3p (Figure 3(d)) and hsa-miR-663a (Figure 3(e)) levels were significantly increased in the blood of colon cancer, values of which were both <0.0001. ROC is shown as Figure 3(f), AUC of which was 0.980. When choosing a threshold of 0.5, the sensitivity of this scoring system was 0.965 and the specificity was 0.965.

4. Discussion

Colon cancer develops insidiously. Screening methods which could find colon cancer patients from population were important. Although guidelines [3] emphasized the importance of noninvasive tools, such as the fecal immunochemical test [21–24], high-sensitivity guaiac-based fecal occult blood test [25, 26], and computed tomography colonography, to screen colon cancer, screening of colon cancer in China still relies on colonoscopies for high-risk population. A new kind of noninvasive way to diagnose CC is urgently needed.

Circulation miRNAs have been regarded as an effective noninvasive screen method for many kinds of cancers [27–29]; some studies about circulation miRNAs have been performed in colon cancer. It has been proved that circulation miRNAs could hold stable under various storage conditions [30], which provided basic theoretical support for miRNAs to become tumor biomarkers.

Several researchers did great jobs in identifying a miRNA biomarker for colon cancer. Sabry et al. [31] focused on HIF-1α-VEGF signaling pathway-related miRNAs, such as miR-210 and miR-21. They found that miR-210 and miR-21 were associated with clinical staging of colon cancer and showing a high diagnostic value with sensitivity and specificity 88.6%, 90.1% and 91.4%, 95.0%, respectively. Bilegsaikhan et al. [32] studied miR-338-5p which has been reported to be upregulated in colon cancer tissues and found that the AUC of it was 0.923 in distinguishing CRC from healthy controls. Moreover, Zhang et al. [33] studied the miRNA profile in the serum of patients with colon cancer by microarray. And they identified 8 DE miRNAs including miR-4463, miR-5704, miR-371b-3p, miR-1247-5p, miR-1293, miR-548at-5p, miR-107, and miR-139-3p, to distinguish colon cancer patients from normal controls. However, no diagnostic tests and AUC analysis were performed in this study. Krawczyk et al. [34] focused on the combination of circulation miR-506 and miR-4316, which showed an AUC value of 0.751 to predict colon cancer.

Their findings were mostly based on the aberrantly expressed circulation miRNAs or miRNA which has been reported upregulating in colon cancer tissues; deep data mining and analysis were lacking. Our study was based on both differentially expressed miRNAs and WGCNA. Instead of only focusing on differentially expressed miRNAs like most of the other researches, WGCNA utilizes information from all of the most varied genes, to identify interested gene module and to analyze their association with phenotypes. By making full usage of gene expression information, we can find hub genes which really matter in the development of cancers. Hub genes we discovered might not have the smallest value in DE analysis; however, it would be more stable and more biologically meaningful than genes selected by values of DE analysis, which could be proved by our validation datasets. This article chose a more powerful dimension to select target genes, which would explain why miR-23 and miR-663 were not selected in the other studies.

We have identified hsa-miR-23a-3p (, ) and hsa-miR-663a (, ) as diagnostic biomarkers in colon cancer using data from GSE61741 and confirmed our scoring system in GSE112264 and GSE106817. The AUC of these two databases was 0.917 and 0.980, respectively. Diagnostic test results reflect that it is reliable and stable for our scoring system to distinguish cancer patients from population, which means that it is possible to regard hsa-miR-23a-3p and hsa-miR-663a as diagnostic biomarkers in colon cancer.

For hsa-miR-663a, Kuroda et al. [35] indicated that after treatment with antimicrobial peptides, miR-663a could regulate cancer development of CC by targeting the CXCR4-p21 pathway. And for hsa-miR-23a-3p, no previous researches were performed to find its function in colon cancer. In hepatocellular cancer (HCC), Liu et al. [36] found that HCC patients with a high body fat ratio have higher miR-23a levels in both serum exosomes and tumor tissues than HCC patients with a low body fat ratio. Their further studies indicated that miR-23a was transported into cancer cells via exosomes, thus promoting HCC cell growth and migration. Our results exhibited that miR-23a was upregulated stably in blood samples of colon cancer patients, and it was also one of the hub genes in cancer-related modules, which indicated that hsa-miR-23a-3p could exhibit similar properties as HCC in colon cancer.

This is the first WGCNA research to identify circulation miRNA in colon cancer. The diagnostic biomarkers we defined were stable and reliable. Our study was based on the GEO database; although advanced analysis methods and multidatabase validation were applied, it still has some limitations because no hospital’s samples were involved. The prospective validation using our hospital’s samples would be performed when our cohorts are strong enough to validate these findings. Besides, we found that two hub miRNAs could be used as a circulation biomarker for colon cancer; however, no functional mining experiments were performed. Our further studies would focus on the function of hsa-miR-23a-3p and hsa-miR-663a in colon cancer tissues and circulation system to find its roles in CC development.

5. Conclusion

In this study, we built a scoring system based on circulation hub miRNAs found by WGCNA to predict the development of colon cancer.

Data Availability

All the data used in this article could be downloaded in the GEO database.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

The authors thank all the persons who also participated in this study. This work was fully supported by the National Natural Science Foundation of China (81702314), the Funding Program for Excellent Talents of Beijing (2017000021469G212), the Beijing Municipal Administration of Hospitals Youth Program (QML20180108), and the Digestive Medical Coordinated Development Center of Beijing Municipal Administration of Hospitals (XXZ02, XXZ01).

Supplementary Materials

Supplementary Figure A: sample clustering. Supplementary Figure B: heat map proving the independence between each gene. Supplementary Figure C: correlation analysis between gene modules and cancer development. Supplementary Figure D: miRNA350-miRNA interaction network in yellow module. (Supplementary Materials)

References

K. D. Miller, L. Nogueira, A. B. Mariotto et al., “Cancer treatment and survivorship statistics, 2019,” CA: A Cancer Journal for Clinicians, vol. 69, no. 5, pp. 363–385, 2019.
View at: Publisher Site | Google Scholar
K. D. Miller, R. L. Siegel, C. C. Lin et al., “Cancer treatment and survivorship statistics, 2016,” CA: a Cancer Journal for Clinicians, vol. 66, no. 4, pp. 271–289, 2016.
View at: Publisher Site | Google Scholar
R. A. Smith, K. S. Andrews, D. Brooks et al., “Cancer screening in the United States, 2019: a review of current American Cancer Society guidelines and current issues in cancer screening,” CA: a Cancer Journal for Clinicians, vol. 69, no. 3, pp. 184–210, 2019.
View at: Publisher Site | Google Scholar
A. M. D. Wolf, E. T. H. Fontham, T. R. Church et al., “Colorectal cancer screening for average-risk adults: 2018 guideline update from the American Cancer Society,” CA: A Cancer Journal for Clinicians, vol. 68, no. 4, pp. 250–281, 2018.
View at: Publisher Site | Google Scholar
A. M. Attallah, M. el-Far, A. R. Ibrahim et al., “Clinical value of a diagnostic score for colon cancer based on serum CEA, CA19-9, cytokeratin-1 and mucin-1,” British Journal of Biomedical Science, vol. 75, no. 3, pp. 122–127, 2018.
View at: Publisher Site | Google Scholar
Q.-R. Huang and X.-B. Pan, “Prognostic lncRNAs, miRNAs, and mRNAs form a competing endogenous RNA network in colon cancer,” Frontiers in Oncology, vol. 9, p. 712, 2019.
View at: Publisher Site | Google Scholar
Q. Lou, R. Liu, X. Yang et al., “miR-448 targets IDO1 and regulates CD8⁺ T cell response in human colon cancer,” Journal for ImmunoTherapy of Cancer, vol. 7, no. 1, p. 210, 2019.
View at: Publisher Site | Google Scholar
W. Zhu, X. Luo, H. Fu, L. Liu, P. Sun, and Z. Wang, “miR-3653 inhibits the metastasis and epithelial-mesenchymal transition of colon cancer by targeting Zeb2,” Pathology - Research and Practice, vol. 215, no. 10, article 152577, 2019.
View at: Publisher Site | Google Scholar
J. Mukohyama, T. Isobe, Q. Hu et al., “miR-221 targets QKI to enhance the tumorigenic capacity of human colorectal cancer stem cells,” Cancer Research, vol. 79, no. 20, pp. 5151–5158, 2019.
View at: Publisher Site | Google Scholar
Q. Ju, Y. J. Zhao, Y. Dong et al., “Identification of a miRNA-mRNA network associated with lymph node metastasis in colorectal cancer,” Oncology Letters, vol. 18, no. 2, pp. 1179–1188, 2019.
View at: Publisher Site | Google Scholar
M. Li, C. Lei, S. Liu et al., “Loss of circulating exosomal miR-92b is a novel biomarker of colorectal cancer at early stage,” International Journal of Medical Sciences, vol. 16, no. 9, pp. 1231–1237, 2019.
View at: Publisher Site | Google Scholar
A. Etheridge, I. Lee, L. Hood, D. Galas, and K. Wang, “Extracellular microRNA: a new source of biomarkers,” Mutation Research, vol. 717, no. 1–2, pp. 85–90, 2011.
View at: Publisher Site | Google Scholar
X. Zhou, Y.-N. Shang, R. Lu, C.-W. Fan, and X.-M. Mo, “HighANKZF1expression is associated with poor overall survival and recurrence-free survival in colon cancer,” Future Oncology, vol. 15, no. 18, pp. 2093–2106, 2019.
View at: Publisher Site | Google Scholar
P. Langfelder and S. Horvath, “WGCNA: an R package for weighted correlation network analysis,” BMC Bioinformatics, vol. 9, no. 1, p. 559, 2008.
View at: Publisher Site | Google Scholar
A. Keller, P. Leidinger, B. Vogel et al., “miRNAs can be generally associated with human pathologies as exemplified for miR-144,” BMC Medicine, vol. 12, no. 1, p. 224, 2014.
View at: Publisher Site | Google Scholar
M. E. Ritchie, B. Phipson, D. Wu et al., “limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, no. 7, pp. e47–e47, 2015.
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
X. Robin, N. Turck, A. Hainard et al., “pROC: an open-source package for R and S+ to analyze and compare ROC curves,” BMC Bioinformatics, vol. 12, no. 1, article 77, 2011.
View at: Publisher Site | Google Scholar
A. Yokoi, J. Matsuzaki, Y. Yamamoto et al., “Integrated extracellular microRNA profiling for ovarian cancer screening,” Nature Communications, vol. 9, no. 1, article 4319, 2018.
View at: Publisher Site | Google Scholar
F. Urabe, J. Matsuzaki, Y. Yamamoto et al., “Large-scale circulating microRNA profiling for the liquid biopsy of prostate cancer,” Clinical Cancer Research, vol. 25, no. 10, pp. 3016–3025, 2019.
View at: Publisher Site | Google Scholar
H. Salimzadeh, F. Bishehsari, C. Sauvaget et al., “Feasibility of colon cancer screening by fecal immunochemical test in Iran,” Archives of Iranian Medicine, vol. 20, no. 12, pp. 726–733, 2017.
View at: Google Scholar
S. Y. Yuan, W. Wu, J. Fu et al., “Quantitative immunochemical fecal occult blood test for neoplasia in colon cancer screening,” Journal of Digestive Diseases, vol. 20, no. 2, pp. 78–82, 2019.
View at: Publisher Site | Google Scholar
J. M. Daly, B. T. Levy, M. L. Merchant, and J. Wilbur, “Mailed fecal-immunochemical test for colon cancer screening,” Journal of Community Health, vol. 35, no. 3, pp. 235–239, 2010.
View at: Publisher Site | Google Scholar
C. A. Doubeni, C. D. Jensen, S. A. Fedewa et al., “Fecal immunochemical test (FIT) for colon cancer screening: variable performance with ambient temperature,” Journal of the American Board of Family Medicine, vol. 29, no. 6, pp. 672–681, 2016.
View at: Publisher Site | Google Scholar
M. S. Cappell, “Reducing the incidence and mortality of colon cancer: mass screening and colonoscopic polypectomy,” Gastroenterology Clinics of North America, vol. 37, no. 1, pp. 129–160, 2008.
View at: Publisher Site | Google Scholar
P. Rozen, E. Ron, Z. Fireman et al., “The relative value of fecal occult blood tests and flexible sigmoidoscopy in screening for large bowel neoplasia,” Cancer, vol. 60, no. 10, pp. 2553–2558, 1987.
View at: Publisher Site | Google Scholar
J. R. Chevillet, T. D. Khokhlova, M. D. Giraldez et al., “Release of cell-free microRNA tumor biomarkers into the blood circulation with pulsed focused ultrasound: a noninvasive, anatomically localized, molecular liquid biopsy,” Radiology, vol. 283, no. 1, pp. 158–167, 2017.
View at: Publisher Site | Google Scholar
J. V. Fritz, A. Heintz-Buschart, A. Ghosal et al., “Sources and functions of extracellular small RNAs in human circulation,” Annual Review of Nutrition, vol. 36, no. 1, pp. 301–336, 2016.
View at: Publisher Site | Google Scholar
L. Min, S. Zhu, L. Chen et al., “Evaluation of circulating small extracellular vesicles derived miRNAs as biomarkers of early colon cancer: a comparison with plasma total miRNAs,” Journal of Extracellular Vesicles, vol. 8, no. 1, article 1643670, 2019.
View at: Publisher Site | Google Scholar
Q. Ge, Y. Zhou, J. Lu, Y. Bai, X. Xie, and Z. Lu, “miRNA in plasma exosome is stable under different storage conditions,” Molecules, vol. 19, no. 2, pp. 1568–1575, 2014.
View at: Publisher Site | Google Scholar
D. Sabry, S. E. M. el-Deek, M. Maher et al., “Role of miRNA-210, miRNA-21 and miRNA-126 as diagnostic biomarkers in colorectal carcinoma: impact of HIF-1α-VEGF signaling pathway,” Molecular and Cellular Biochemistry, vol. 454, no. 1–2, pp. 177–189, 2019.
View at: Publisher Site | Google Scholar
E. Bilegsaikhan, H. N. Liu, X. Z. Shen, and T. T. Liu, “Circulating miR-338-5p is a potential diagnostic biomarker in colorectal cancer,” Journal of Digestive Diseases, vol. 19, no. 7, pp. 404–410, 2018.
View at: Publisher Site | Google Scholar
Y. Zhang, M. Li, Y. Ding et al., “Serum microRNA profile in patients with colon adenomas or cancer,” BMC Med Genomics, vol. 10, no. 1, p. 23, 2017.
View at: Publisher Site | Google Scholar
P. Krawczyk, T. Powrózek, T. Olesiński et al., “Evaluation of miR-506 and miR-4316 expression in early and non-invasive diagnosis of colorectal cancer,” International Journal of Colorectal Disease, vol. 32, no. 7, pp. 1057–1060, 2017.
View at: Publisher Site | Google Scholar
K. Kuroda, T. Fukuda, M. Krstic-Demonacos et al., “miR-663a regulates growth of colon cancer cells, after administration of antimicrobial peptides, by targeting CXCR4-p21 pathway,” BMC Cancer, vol. 17, no. 1, article 33, 2017.
View at: Publisher Site | Google Scholar
Y. Liu, J. Tan, S. Ou, J. Chen, and L. Chen, “Adipose-derived exosomes deliver miR-23a/b to regulate tumor growth in hepatocellular cancer by targeting the VHL/HIF axis,” Journal of Physiology and Biochemistry, vol. 75, no. 3, pp. 391–401, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Da Qin 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

1285

Downloads

1271

Citations