Abstract
About 20 percent of TNM-stage II colon cancer patients who are treated by surgical resection develop recurrence, and adjuvant chemotherapy in this group is still debated among researchers and clinicians. Currently, adverse histopathological and clinical factors are used to select patients for adjuvant chemotherapy following surgery. However, additional biomarkers to classify patients at risk of recurrence are needed. We have conducted a study using fresh frozen tumor tissue from 54 TNM-stage II colon cancer patients and performed microRNA profiling using next-generation sequencing. For the selection of the prognostic microRNAs, a LASSO Cox Regression model was employed. For the validation, we used the publically available TCGA-COAD cohort (n = 122). A prognostic panel of four micorRNAs (hsa-miR-5010-3p, hsa-miR-5100, hsa-miR-656-3p and hsa-miR-671-3p) was identified in the study cohort and validated in the TCGA-COAD cohort. The four-microRNA classifier successfully identified high-risk patients in the study cohort (P < 0.001) and the validation cohort (P = 0.005). Additionally, a number of established risk factors and the four-miRNA classifier were used to construct a nomogram to evaluate risk of recurrence. We identified a four-microRNA classifier in patients with TNM-stage II colon cancer that can be used to discriminate between patients at low- and high risk of recurrence.
Similar content being viewed by others
Introduction
Colon cancer is one of the most common cancers worldwide and is the second most common cancer in developed countries1,2. Surgery and chemotherapy are the most common treatments for colon cancer patients, and the TNM (Tumor-Node -Metastasis) - staging system has been widely used for colon cancer management allocation clinically3. Adjuvant chemotherapy after surgery is a standard regimen for patients with TNM-stage III to improve survival and reduce the risk of tumor recurrence4. In contrast, the benefit of adjuvant chemotherapy for TNM-stage II patients is still debated, and a quarter of these patients die from tumor recurrence after surgery5,6,7,8. Currently, clinicopathological factors do not clearly guide patients‘ prognostic stratification, nor do they identify the patients with benefit from adjuvant chemotherapy in TNM-stage II colon cancer9,10.
One of the most recently established prognostic factors in early stages of colon cancer is a defective DNA mismatch repair (MMR) system in the tumor cells. Failure of the MMR system can be tested by immunohistochemistry for loss of MMR proteins, or by testing for microsatellite instability (MSI) by PCR. A study based on a large series of unselected cases found that immunohistochemistry and microsatellite instability are comparable methodologies for detecting defective MMR system11. The prevalence of MMR system deficiency is higher in early stages (20% in stage I-II) and it is associated with better survival in TNM-stage II colon cancer12,13,14,15,16,17,18. Importantly, it is known that patients in TNM-stage II with MSI does not benefit from adjuvant 5-flurouracil based chemotherapy19,20,21. Other established clinicopathological risk factors for TNM-stage II colon cancer patients are pT4 and low total number of lymph nodes removed22,23. Unfortunately, these factors are not able to successfully identify all patients at high risk of recurrence. Consequently, there is a need for biomarkers to add prognostic value to the traditional clinicopathological factors and to guide the clinicians to discriminate high-risk and low-risk patients. As a result of the modest therapeutic benefit of adjuvant chemotherapy in TNM-stage II colon cancers, prognostic biomarkers have a major relevance.
Recently, microRNAs (miRNAs) are identified as prognostic biomarkers in TNM-stage II colon cancer patients, such as miR-21, miR-29a, miR-106a, miR-55624,25,26,27. A study identified a classifier based on six-miRNAs as a prognostic biomarker in TNM-stage II colon cancer patients28. Another study aimed to reproduce this study, identified three of these miRNAs as the best prognostic marker29. In our previous study, we identified a prognostic 16-miRNA signature among 83 selected miRNAs in TNM-stage II and III colon cancer patients30. These studies have used microarrays and quantitative real-time PCR (RT-qPCR).
In order to identify a prognostic biomarker with the ability to classify TNM-stage II patients at high- and low- risk of recurrence, we here investigated global miRNA expression profiling using next-generation sequencing (NGS).
Results
Cohorts
A summary of the demographical features for the 54 patients in the study (HDH-CC) cohort and the 122 patients in the validation (TCGA-COAD) cohort is shown in the Table 1. The median age of the HDH-CC cohort was 73 years. The patients of the TCGA-COAD cohort were younger (median age 67.8 years). There was no difference between fraction of women and men in the HDH-CC cohort (52% women and 48% men) and in the TCGA-COAD cohort (57% men and 43% women). In the HDH-CC cohort, 65% of the patients had proficient MMR status, and 59% of the TCGA-COAD cohort had MSS. MSI-status is used instead of MMR-status in the validation cohort because immunohistochemistry for the MMR proteins was not available for all the patients. Eighty five percent of the patients had well or moderate differentiated tumors and 93% had adenocarcinoma in the HDH-CC cohort. Information regarding tumor differentiation status was not available for the TCGA-COAD cohort. A total of 13% and 19% of the patients had recurrence in the HDH-CC cohort and TCGA-COAD cohort, respectively. The association between recurrence and the four-miRNA classifier was statistically significant (P < 0.001 and P = 0.006, respectively). The relation between tumor stage and the four-miRNA classifier was statistically significant in the HDH-CC cohort (P = 0.04). There were no statistically significant associations between the four-miRNA classifier and the other clinicopathological variables (Table 1).
Identification and construction of the four-miRNA classifier
In order to find the prognostic miRNAs in the HDH-CC cohort, we assessed TNM-stage II colon cancer patients using a LASSO Cox regression analysis model. A panel of four miRNAs (hsa-miR-5010-3p, hsa-miR-5100, hsa-miR-656-3p and hsa-miR-671-3p) was found to be the best predictor of recurrence in TNM-stage II colon cancer patients. The dichotomized miRNA expression values weighted by the coefficients from the LASSO regression generated a prognostic index for each patient using the formula where E-miR-n was the dichotomized expression of miRNAn.
The prognostic index was applied to all TNM-stage II colon cancer patients in the HDH-CC cohort and TCGA-COAD cohort. ROC-curve analyses were performed in order to generate the best cut-point of the prognostic index with optimal sensitivity and specificity at which low-risk patients can be distinguished from high-risk patients. The area under curve (AUC) in the ROC analyses for the HDH-CC and TCGA-COAD cohorts were statistically significant (AUC = 0.901, 95% CI: 0.75–1.00; P = 0.001; Fig. 1A and AUC = 0.658, 95% CI: 0.53–0.78; P = 0.0.19; Fig. 1B, respectively).
Eighty percent of the patients in the HDH-CC cohort and 64% of the patients in the TCGA-COAD cohort have a low risk of recurrence, according to our four-miRNA classifier (Table 1). We did not observe any statically significant associations between the level of the classifier and other clinicopathological features in the HDH-CC and TCGA-COAD cohort (Table 1).
The prognostic value of the four-miRNA classifier in the HDH-CC cohort
Association between the four-miRNA classifier level and 5-year DFS were evaluated by performing univariate and multivariate Cox regression analyses. For the patients in the HDH-CC cohort, Kaplan-Meier analyses showed statistically significant better 5-year DFS for low versus high level of the four-miRNA classifier (97% vs 42%, P < 0.001; Fig. 1C). The sensitivity and specificity of the four-miRNA classifier for detecting recurrence were 86% and 90%, respectively. Multivariate Cox-regression analyses were not performed in the HDH-CC cohort because of the small number of events (recurrence, n = 7).
Validation of the prognostic value of the four-miRNA classifier in an independent dataset (TCGA-COAD)
The prognostic value of the four-miRNA classifier was tested in the publically available TCGA-COAD dataset. In agreement with the findings for the HDH-CC cohort, high-risk patients had statistically significant shorter 5-year DFS than low-risk patients. Thus, patients with a low-level of the four-miRNA classifier have better 5-year DFS than patients with high-level, 82% versus 56% (P = 0.005; Fig. 1D). The sensitivity of the four-miRNA classifier was 61% and the specificity was 72% in this cohort. The multivariate Cox regression analyses showed that the four-miRNA classifier remained a statistically significant prognostic factor when adjusted for age, gender, MSI-status, total number of lymph node (TLN), and T-stage (HR 2.88; 95% CI: 1.20–6.94; P = 0.023; Table 2).
Construction of nomogram based on the four-miRNA classifier
We constructed a nomogram to predict the risk of recurrence in TNM-stage II colon cancer patients and to provide the clinician with a better estimation of the prognosis of the individual patient. The nomogram integrated age, gender, MSI-status, TLN, T-stage and the four-miRNA classifier (Fig. 2A). The nomogram identified T-stage and the four-miRNA classifier as the largest contributors to prognosis, followed by gender, age, TLN, and MSI. We calculated the risk score by adding up the scores from each variable, calculated the Total Point score and drew a straight line down to determine the estimated risk of recurrence. The prognostic accuracy of the nomogram was assessed using ROC analyses (Fig. 2B). The cut-point for the risk of recurrence determined by ROC curve analyses resulted in a sensitivity and specificity of 96% and 46%, respectively in the TCGA-COAD cohort. The KM plot of the TCGA-COAD cohort showed that the stratification of the patients in low- and high risk groups is statistically significant (P = 0.001; Fig. 2C). Further, the nomogram was tested in the HDH-CC cohort, using the same cut-point as for the TCGA-COAD cohort. The results showed that the nomogram identified the patients experiencing recurrence with a sensitivity of 100% and a specificity of 85% (P < 0.001; Fig. 2D). Waterfall plots were generated in order to illustrate the risk score for every patient in both cohorts based on the four-miRNA classifier (Fig. 3A,B), and based on the nomogram (Fig. 3C,D).
Discussion
Clinical trials have shown that although 20% of TNM-stage II colon cancer patients are at high risk of recurrence, the benefit of current adjuvant chemotherapy treatment is small in TNM-stage II colon cancer patients as a whole. In the current study, we constructed and validated a four-miRNA classifier as a prognostic biomarker to predict recurrence in TNM-stage II colon cancer patients. The four-miRNA classifier was developed using LASSO cox regression in the HDH-CC cohort (n = 54) and validated in TCGA-COAD cohort (n = 122). The classifier successfully discriminated between patients at high- and low risk of recurrence. High-risk patients had statistically significant shorter 5-year DFS than low-risk patients. Accordingly, a low level of the four-miRNA classifier was associated with better survival in both HDH-CC and TCGA-COAD cohort. The sensitivity and specificity of the four-miRNA classifier appears to be better in the HDH-CC cohort (86% and 90%, respectively) than in the TCGA-COAD cohort (61% and 72%, respectively). The 5-year DFS in the TCGA-COAD is lower than in the HDH-CC cohort (72% versus 85%, respectively), which may have an impact on the differences. The classifier remained statistically significant in multivariate analyses in the TCGA-COAD cohort. Because of the relative small number of patients and few recurrences in the HDH-CC cohort, we did not perform multivariate analyses in this cohort. Nomograms have been developed to evaluate several significant clinicopathological risk factors. We constructed a nomogram including the four-miRNA classifier to develop a practical method for clinicians to predict the probability of TNM-stage II colon cancer recurrence. Compared to the classifier alone, the prognostic nomogram improved the sensitivity from 86% to 100% in the HDH-CC cohort, and from 61% to 96% in the TCGA-COAD cohort. On the other hand, the nomogram does not improve the specificity; from 90% to 85% in the HDH-CC and from 72% to 46% in the TCGA-COAD cohort.
Therefore, we argue that the four-miRNA classifier alone is able to predict the probability of recurrence; however, we can not exclude the possibility that the nomogram could be more robust than the classifier alone for identifying patients at risk for recurrence in other colon cancer patient cohorts.
In clinical practice, histopathological and clinical factors are used to determine prognosis and treatment in colon cancer23. However, histologic tumor staging has not been consistently associated with increased risk of recurrence in TNM-stage II colon cancer31. Further, MMR or MSI status are documented as the most validated and important prognostic biomarkers in TNM-stage II colon cancer patients12,13. Patients with MSI have low recurrence rates and good outcomes without adjuvant chemotherapy17,21. A previous study describes different miRNA expression pattern in deficient and proficient MMR32. In the present study, there was no statistically significant association between MMR status and the four-miRNA classifier (Table 1). In the HDH-CC cohort, 13 out of 15 patients (87%) with deficient MMR were classified as low risk, and 27 out of 35 patients (77%) with proficient MMR were also classified as low risk patients (Supplementary Table 1). Therefore, the good prognostic value of the four-miRNA classifier is not associated with MMR-status in this study. Further, the four-miRNA classifier was not able to classify patients into low- and high risk groups in TNM-stage III colon cancer patients in the two cohorts (Supplementary Fig. 3).
Therefore, the prognostic value of the four-miRNA classifier is specific for TNM-stage II colon cancer patients.
In addition, several miRNA studies have identified miRNA signatures as prognostic biomarkers in colon cancer cancer28,29,30,33. The results of the tissue-based miRNA signature studies are based on microarray technology or RT-qPCR. Earlier, we have identified a 16-miRNA signature as a prognostic biomarker for TNM-stage II and III colon cancer patients. The study was based on RT-qPCR and the panel we used included only 83 selected, cancer related miRNAs30. The current study is based on next-generation sequencing which has the advantage of global miRNA analysis, higher sensitivity and the ability to detect novel miRNAs34. Interestingly, the four miRNAs included in the classifier have not previously been documented as prognostic markers in colon cancer. However, there are several studies linking these four miRNAs to other cancer types or tumor models. One previous study found an increased level of miR-5100 in non-small-cell lung cancer and its association with poor prognosis was documented35. A functional study showed that this miRNA promotes tumor growth as a result of targeting Rab6, a protein involved in membrane traffic from the Golgi apparatus to the endoplasmatic reticulum and exocytosis36. Downregulation of miR-656 in cluster with miR-379 were found in multiple cancers in a genome-wide study37. The miR-379/miR-656 cluster was also identified as a tumor suppressor locus in adult medulloblastomas38. Another study identified an eleven-miRNA signature including miR-671 as a prognostic biomarker for a subtype of kidney cancer with wild type BRCA1 associated protein-1 (BAP-1)39. Regarding miR-5010, no functional or clinical studies have been published yet; however, a number of genes are predicted to be targets of this miRNA with high relevance to tumor progression, among others, member of RAS oncogene family (RAB9A), APC membrane recruitment protein 1 (AMER1), beta-catenin interacting protein 1 (CTNNBIP1), Ras-related GTP binding D (RRAGD), etc. (miRDB; www.mirdb.org and TargetMiner; www.isical.ac.in).
In summary, we constructed a classifier based on expression of four miRNAs in fresh frozen tumor tissue from TNM-stage II colon cancer patients. The classifier is able to distinguish between high- and low-risk patients, which enables the clinicians to identify the patients that need to be followed closely in order to detect any recurrence at an early stage and administrate chemotherapy to patients at high-risk of recurrence.
Patients and Methods
Patient cohorts
This study was approved by The Regional Ethics Committee West and the Data Inspectorate for National Registries. The patients signed the informed consent and the original clinical trial was registered at U.S National Institute of Health (ClinicalTrials.gov, NCT00963352). All methods were conducted in accordance with approved guidelines and regulations.
Three hundred and ninety six TNM stage I-IV colon cancer patients had undergone surgical resection for primary tumor between January 2007 and December 2011 at Haraldsplass Deaconess Hospital (HDH-CC; Supplementary Fig. 1)40. We performed NGS of miRNAs in 128 TNM-stage I-IV tumor samples and 10 adjacent normal colon samples with available fresh frozen tissue. Out of 128 patients, 54 patients were classified as TNM-stage II. None of these patients received radiotherapy or adjuvant chemotherapy41. The patients had a regular check-up, which included clinical examination, chest x-ray, or CT, abdominal ultrasound or CT and CEA blood testing until the fifth year of the follow-up after surgery. Colonoscopy was performed three years after the date of surgery. The four mismatch repair proteins were assessed by immunohistochemistry using MLH1 (1:60, DAKO, M3640, clone ES05), MSH-2 (1:300, Biocare Medical, CM219B, clone FE11), MSH6 (1:50, DAKO, M3646, clone EP49) and PMS2 (1:50, DAKO, M3647, clone EP51). Epitope retrieval (HIER) was performed in TE-buffer pH9, and the detection kit used was MACH3 HRP-Polymer (Biocare Medical).
For validation of our findings, we downloaded miRNA-sequencing data from the TCGA-COAD cohort from NIH National Cancer Institute Genomic Data Commons Data Portal (The Cancer Genome Atlas; https://cancergenome.nih.gov/). We included all 122 TNM-stage II colon cancer patients for which miRNA sequencing results and the requested clinicopathological and survival information were available (Supplementary Fig. 2). MSI status for the TCGA-COAD patients was downloaded from the TCGA Data Portal (NIH). According to information available from the TCGA consortium, a panel of four mononucleotide repeat loci (BAT25, BAT26, BAT40, and transforming growth factor receptor type II) and three dinucleotide repeat loci (CA repeats in D2S123, D5S346, and D17S250) were assessed42. The patients were classified as MSI-high if more than 40% of the markers were altered, and patients with lower than 40%-altered marker were classified as MSS or MSI-low.
RNA extraction
Total RNA from the fresh frozen tissue samples was extracted using the miRNeasy, Mini kit (QIAGEN) according to the manufacturer’s protocol. Approximately 30 mg of each sample was homogenized in 700 µl QIAzol Lysis Reagent using Tissuelyzer (QIAGEN) for 10 minutes at 25 Hz. The homogenate was processed according to the manufacturer’s instruction and the total RNA was purified using a protocol where DNase treatment was included. RNA concentrations were measured by NanoDrop and the quality of RNA was measured by Agilent RNA Bioanalyzer (Agilent RNA 6000 Nano Assay protocol-edition April 2007). The RNA samples were stored at −80 °C until use.
MicroRNA sequencing
The Illumina TruSeq RNA sample Prep Kit (Illumina, Inc., San Diego, CA, USA) was used to construct cDNA libraries for sequencing. In brief, one µg of the total RNA was used to sequentially ligate 3′ and 5′ adapters to the ends and reverse transcribe to generate cDNA. The cDNA was amplified using 15 PCR cycles. The amplified cDNA was purified and quantified on a 6% Novex TBE gel. The quality control was performed using Agilent Technologies 2100 Bioanalyzer. The resulting pooled libraries were sequenced on Illumina HighSeq4000 (Illumina, Inc., San Diego, CA, USA).
Processing of the sequencing data
The raw miRNA sequences were output in FASTQ format. The raw reads were initially filtered for reads containing adapter sequences and low-quality sequences. miRDeep2 were used to map the sequences to miRNAs and quantify expression of the specific miRNAs43. The counts were then processed by edgeR Bioconductor package using the Trimmed Mean of M (TMM) method for normalization44.
Statistical analysis
The R software (version 3.3.2) and the “glmnet” package were used to perform the LASSO regression model analysis in the HDH-CC cohort. The penalized Cox regression model was determined through 200-times cross validation45,46. Based on the optimal lambda value, a panel of four miRNAs with associated coefficients was identified as the best prognostic marker for recurrence among the TNM-stage II colon cancer patients in the HDH-CC cohort. The risk score for each patient was calculated based on the dichotomized expression of each prognostic miRNA using X-tile plots (software version 3.6.1, Yale University School of Medicine, New Haven, CT, USA) and its associated regression coefficient47. The patients in the HDH-CC and TCGA-COAD cohort were split into low-risk or high-risk groups according to the receiver operating characteristic (ROC) curve. The disease- free survival (DFS) time of patients were from date of surgery to the time of recurrence or colon cancer related death before documented recurrence. Deaths of other causes were censored and maximum follow-up time was 5-years. Survival differences among the risk groups were analyzed by the Kaplan-Meier method with log-rank test, and the multivariate survival analyses were performed using Cox-regression model. ROC-curves and survival analyses were performed using SPSS (IBM SPSS Statistics, version 24) and all P-values considered statistically significant if P < 0.05. We designed a nomogram to provide the overall probability of risk of recurrence using the “rms” package of R software (version 3.3.2)48. The nomogram was generated in the TCGA-COAD cohort because the number of patients and recurrent events in this cohort were larger than in the study cohort. ROC curves were used to assess the accuracy of the nomogram in the TCGA-COAD - and the HDH-CC cohort.
References
Ferlay, J. et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. International journal of cancer 136, E359–386, https://doi.org/10.1002/ijc.29210 (2015).
DeSantis, C. E. et al. Cancer treatment and survivorship statistics, 2014. CA: a cancer journal for clinicians 64, 252–271, https://doi.org/10.3322/caac.21235 (2014).
Kawaguchi, T. et al. Comparison of prognostic compatibility between seventh AJCC/TNM of the esophagus and 14th JCGC staging systems in Siewert type II adenocarcinoma. Anticancer research 33, 3461–3465 (2013).
Gill, S. et al. Pooled analysis of fluorouracil-based adjuvant therapy for stage II and III colon cancer: who benefits and by how much? Journal of clinical oncology: official journal of the American Society of Clinical Oncology 22, 1797–1806, https://doi.org/10.1200/jco.2004.09.059 (2004).
Edge, S. B. & Compton, C. C. The American Joint Committee on Cancer: the 7th edition of the AJCC cancer staging manual and the future of TNM. Annals of surgical oncology 17, 1471–1474, https://doi.org/10.1245/s10434-010-0985-4 (2010).
Benson, A. B. 3rd et al. American Society of Clinical Oncology recommendations on adjuvant chemotherapy for stage II colon cancer. Journal of clinical oncology: official journal of the American Society of Clinical Oncology 22, 3408–3419, https://doi.org/10.1200/jco.2004.05.063 (2004).
Figueredo, A., Coombes, M. E. & Mukherjee, S. Adjuvant therapy for completely resected stage II colon cancer. The Cochrane database of systematic reviews, Cd005390, https://doi.org/10.1002/14651858.CD005390.pub2 (2008).
Ratkin, G. A. Adjuvant therapy for colon and rectal cancer. American family physician 55(2487–2492), 2495–2486 (1997).
O’Connor, E. S. et al. Adjuvant chemotherapy for stage II colon cancer with poor prognostic features. Journal of clinical oncology: official journal of the American Society of Clinical Oncology 29, 3381–3388, https://doi.org/10.1200/jco.2010.34.3426 (2011).
Morris, E. J., Maughan, N. J., Forman, D. & Quirke, P. Who to treat with adjuvant therapy in Dukes B/stage II colorectal cancer? The need for high quality pathology. Gut 56, 1419–1425, https://doi.org/10.1136/gut.2006.116830 (2007).
Lindor, N. M. et al. Immunohistochemistry Versus Microsatellite Instability Testing in Phenotyping Colorectal Tumors. Journal of Clinical Oncology 20, 1043–1048, https://doi.org/10.1200/jco.2002.20.4.1043 (2002).
Kawakami, H., Zaanan, A. & Sinicrope, F. A. Microsatellite instability testing and its role in the management of colorectal cancer. Current treatment options in oncology 16, 30, https://doi.org/10.1007/s11864-015-0348-2 (2015).
Sargent, D. J. et al. Prognostic impact of deficient mismatch repair (dMMR) in 7,803 stage II/III colon cancer (CC) patients (pts): A pooled individual pt data analysis of 17 adjuvant trials in the ACCENT database. Journal of Clinical Oncology 32, 3507–3507, https://doi.org/10.1200/jco.2014.32.15_suppl.3507 (2014).
Hutchins, G. et al. Value of mismatch repair, KRAS, and BRAF mutations in predicting recurrence and benefits from chemotherapy in colorectal cancer. Journal of clinical oncology: official journal of the American Society of Clinical Oncology 29, 1261–1270, https://doi.org/10.1200/jco.2010.30.1366 (2011).
Roth, A. D. et al. Integrated analysis of molecular and clinical prognostic factors in stage II/III colon cancer. Journal of the National Cancer Institute 104, 1635–1646, https://doi.org/10.1093/jnci/djs427 (2012).
Sinicrope, F. A. et al. Prognostic impact of deficient DNA mismatch repair in patients with stage III colon cancer from a randomized trial of FOLFOX-based adjuvant chemotherapy. Journal of clinical oncology: official journal of the American Society of Clinical Oncology 31, 3664–3672, https://doi.org/10.1200/jco.2013.48.9591 (2013).
Gavin, P. G. et al. Mutation profiling and microsatellite instability in stage II and III colon cancer: an assessment of their prognostic and oxaliplatin predictive value. Clinical cancer research: an official journal of the American Association for Cancer Research 18, 6531–6541, https://doi.org/10.1158/1078-0432.ccr-12-0605 (2012).
Roth, A. D. et al. Prognostic role of KRAS and BRAF in stage II and III resected colon cancer: results of the translational study on the PETACC-3, EORTC 40993, SAKK 60-00 trial. Journal of clinical oncology: official journal of the American Society of Clinical Oncology 28, 466–474, https://doi.org/10.1200/jco.2009.23.3452 (2010).
Merok, M. A. et al. Microsatellite instability has a positive prognostic impact on stage II colorectal cancer after complete resection: results from a large, consecutive Norwegian series. Annals of oncology: official journal of the European Society for Medical Oncology/ESMO 24, 1274–1282, https://doi.org/10.1093/annonc/mds614 (2013).
Dienstmann, R., Salazar, R. & Tabernero, J. Personalizing colon cancer adjuvant therapy: selecting optimal treatments for individual patients. Journal of clinical oncology: official journal of the American Society of Clinical Oncology 33, 1787–1796, https://doi.org/10.1200/jco.2014.60.0213 (2015).
Ribic, C. M. et al. Tumor microsatellite-instability status as a predictor of benefit from fluorouracil-based adjuvant chemotherapy for colon cancer. The New England journal of medicine 349, 247–257, https://doi.org/10.1056/NEJMoa022289 (2003).
Snaebjornsson, P. et al. pT4 stage II and III colon cancers carry the worst prognosis in a nationwide survival analysis. Shepherd’s local peritoneal involvement revisited. International journal of cancer 135, 467–478, https://doi.org/10.1002/ijc.28676 (2014).
Gray, R. et al. Adjuvant chemotherapy versus observation in patients with colorectal cancer: a randomised study. Lancet (London, England) 370, 2020–2029, https://doi.org/10.1016/s0140-6736(07)61866-2 (2007).
Hansen, T. F. et al. Redefining high-risk patients with stage II colon cancer by risk index and microRNA-21: results from a population-based cohort. British journal of cancer 111, 1285–1292, https://doi.org/10.1038/bjc.2014.409 (2014).
Oue, N. et al. High miR-21 expression from FFPE tissues is associated with poor survival and response to adjuvant chemotherapy in colon cancer. International journal of cancer 134, 1926–1934, https://doi.org/10.1002/ijc.28522 (2014).
Bullock, M. D. et al. Stratifying risk of recurrence in stage II colorectal cancer using deregulated stromal and epithelial microRNAs. Oncotarget 6, 7262–7279, https://doi.org/10.18632/oncotarget.3225 (2015).
Weissmann-Brenner, A. et al. Tumor microRNA-29a expression and the risk of recurrence in stage II colon cancer. International journal of oncology 40, 2097–2103, https://doi.org/10.3892/ijo.2012.1403 (2012).
Zhang, J. X. et al. Prognostic and predictive value of a microRNA signature in stage II colon cancer: a microRNA expression analysis. The Lancet. Oncology 14, 1295–1306, https://doi.org/10.1016/s1470-2045(13)70491-1 (2013).
Caritg, O. et al. Identifying High-Risk Stage II Colon Cancer Patients: A Three-MicroRNA-Based Score as a Prognostic Biomarker. Clinical colorectal cancer 15, e175–e182, https://doi.org/10.1016/j.clcc.2016.04.008 (2016).
Jacob, H. et al. Identification of a sixteen-microRNA signature as prognostic biomarker for stage II and III colon cancer. Oncotarget 8, 10 (2017).
Gray, R. G., Stanisavljevic, L. S. K., Hestetun, K., Dahl, O. & Myklebust, M. P. Validation study of a quantitative multigene reverse transcriptase-polymerase chain reaction assay for assessment of recurrence risk in patients with stage II colon cancer. Journal of clinical oncology: official journal of the American Society of Clinical Oncology 29, 4611–4619, https://doi.org/10.1200/jco.2010.32.8732 (2011).
Balaguer, F. et al. Colorectal cancers with microsatellite instability display unique miRNA profiles. Clinical cancer research: an official journal of the American Association for Cancer Research 17, 6239–6249, https://doi.org/10.1158/1078-0432.ccr-11-1424 (2011).
Bobowicz, M. et al. Prognostic value of 5-microRNA based signature in T2-T3N0 colon cancer. Clinical & experimental metastasis 33, 765–773, https://doi.org/10.1007/s10585-016-9810-1 (2016).
Berezikov, E., Cuppen, E. & Plasterk, R. H. Approaches to microRNA discovery. Nature genetics 38(Suppl), S2–7, https://doi.org/10.1038/ng1794 (2006).
Wang, T., Liu, X., Tian, Q., Liang, T. & Chang, P. Increasing expression of miR-5100 in non-small-cell lung cancer and correlation with prognosis. European review for medical and pharmacological sciences 21, 3592–3597 (2017).
Huang, H. et al. miR-5100 promotes tumor growth in lung cancer by targeting Rab6. Cancer letters 362, 15–24, https://doi.org/10.1016/j.canlet.2015.03.004 (2015).
Laddha, S. V. et al. Genome-wide analysis reveals downregulation of miR-379/miR-656 cluster in human cancers. Biology direct 8, 10, https://doi.org/10.1186/1745-6150-8-10 (2013).
Kaur, K. et al. Clinicopathological characteristics, molecular subgrouping, and expression of miR-379/miR-656 cluster (C14MC) in adult medulloblastomas. Journal of neuro-oncology 130, 423–430, https://doi.org/10.1007/s11060-016-2250-6 (2016).
Ge, Y. Z. et al. A BAP1 Mutation-specific MicroRNA Signature Predicts Clinical Outcomes in Clear Cell Renal Cell Carcinoma Patients with Wild-type BAP1. Journal of Cancer 8, 2643–2652, https://doi.org/10.7150/jca.20234 (2017).
Storli, K. E., Sondenaa, K., Furnes, B. & Eide, G. E. Outcome after introduction of complete mesocolic excision for colon cancer is similar for open and laparoscopic surgical treatments. Digestive surgery 30, 317–327, https://doi.org/10.1159/000354580 (2013).
Sondenaa, K. et al. The rationale behind complete mesocolic excision (CME) and a central vascular ligation for colon cancer in open and laparoscopic surgery: proceedings of a consensus conference. International journal of colorectal disease 29, 419–428, https://doi.org/10.1007/s00384-013-1818-2 (2014).
TCGA. Microsatellite data; https://wiki.nci.nih.gov/display/TCGA/Microsatellite+data (2012).
Friedlander, M. R., Mackowiak, S. D., Li, N., Chen, W. & Rajewsky, N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic acids research 40, 37–52, https://doi.org/10.1093/nar/gkr688 (2012).
Robinson, M. D. & Oshlack, A. A scaling normalization method for differential expression analysis of RNA-seq data. Genome biology 11, R25, https://doi.org/10.1186/gb-2010-11-3-r25 (2010).
Tibshirani, R. Regression Shrinkage and Selection via the Lasso. Journal of the Royal statistical society 58, 267–288 (1996).
Goeman, J. J. L1 penalized estimation in the Cox proportional hazards model. Biometrical journal. Biometrische Zeitschrift 52, 70–84, https://doi.org/10.1002/bimj.200900028 (2010).
Camp, R. L., Dolled-Filhart, M. & Rimm, D. L. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clinical cancer research: an official journal of the American Association for Cancer Research 10, 7252–7259, https://doi.org/10.1158/1078-0432.ccr-04-0713 (2004).
Harrell, F. Regression modeling stratigies: with application to linear models, logistic and ordinal regression, and survival analysis. Springer (2015).
Acknowledgements
This work was supported by the University of Bergen, Norway, and Western Norway Regional Health Authority.
Author information
Authors and Affiliations
Contributions
H.J. and M.P.M. conceived the experiments, wrote the manuscript, developed the methodology and performed data analysis. O.D. developed the methodology and wrote the manuscript. L.S. and K.E.H. provided clinical data. K.E.S. provided clinical samples.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Jacob, H., Stanisavljevic, L., Storli, K.E. et al. A four-microRNA classifier as a novel prognostic marker for tumor recurrence in stage II colon cancer. Sci Rep 8, 6157 (2018). https://doi.org/10.1038/s41598-018-24519-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-24519-4
This article is cited by
-
A novel miRNA-based classification model of risks and stages for clear cell renal cell carcinoma patients
BMC Bioinformatics (2021)
-
A seven-gene signature model predicts overall survival in kidney renal clear cell carcinoma
Hereditas (2020)
-
Identification of an alternative splicing signature as an independent factor in colon cancer
BMC Cancer (2020)
-
Prognostic Value of MicroRNAs in Stage II Colorectal Cancer Patients: A Systematic Review and Meta-Analysis
Molecular Diagnosis & Therapy (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.