Next Article in Journal
Molecular Characterization of Esophageal Squamous Cell Carcinoma Using Quantitative Proteomics
Next Article in Special Issue
High OCT4 Expression Might Be Associated with an Aggressive Phenotype in Rectal Cancer
Previous Article in Journal
Al-Biruni Earth Radius Optimization with Transfer Learning Based Histopathological Image Analysis for Lung and Colon Cancer Detection
Previous Article in Special Issue
Knockdown of UBQLN1 Functions as a Strategy to Inhibit CRC Progression through the ERK-c-Myc Pathway
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of a Twelve-microRNA Signature with Prognostic Value in Stage II Microsatellite Stable Colon Cancer

by
Ferran Moratalla-Navarro
1,2,3,4,
Anna Díez-Villanueva
1,2,
Ainhoa Garcia-Serrano
5,
Adrià Closa
6,
David Cordero
1,2,
Xavier Solé
7,8,
Elisabet Guinó
1,2,
Rebeca Sanz-Pamplona
1,2,3,9,
Xavier Sanjuan
2,10,
Cristina Santos
2,11,12,
Sebastiano Biondo
2,4,13,
Ramón Salazar
2,4,11,12 and
Victor Moreno
1,2,3,4,*
1
Oncology Data Analytics Program, Catalan Institute of Oncology (ICO), 08908 Barcelona, Spain
2
Colorectal Cancer Group, Bellvitge Biomedical Research Institute (IDIBELL), 08908 Barcelona, Spain
3
Consortium for Biomedical Research in Epidemiology and Public Health (CIBERESP), 28029 Madrid, Spain
4
Department of Clinical Sciences, Faculty of Medicine, University of Barcelona (UB), 08907 Barcelona, Spain
5
Department of Clinical Science, Intervention and Technology (CLINTEC), Karolinska Institutet, 14186 Stockholm, Sweden
6
Department of Pathology, Netherlands Cancer Institute, 1066 CX Amsterdam, The Netherlands
7
Molecular Biology CORE, Center for Biomedical Diagnostics, Hospital Clinic de Barcelona, 08036 Barcelona, Spain
8
Translational Genomic and Targeted Therapeutics in Solid Tumors, August Pi i Sunyer Biomedical Research Institute (IDIBAPS), 08036 Barcelona, Spain
9
Lozano Blesa University Hospital, Aragon Health Research Institute (IISA), Aragon I+D Foundation (ARAID), Government of Aragon, 50009 Zaragoza, Spain
10
Department of Pathology, Bellvitge University Hospital, 08907 Barcelona, Spain
11
Oncology Service, Catalan Institute of Oncology (ICO), 08908 Barcelona, Spain
12
Consortium for Biomedical Research in Oncology (CIBERONC), 28029 Madrid, Spain
13
Department of General and Digestive Surgery, Bellvitge University Hospital, 08907 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(13), 3301; https://doi.org/10.3390/cancers15133301
Submission received: 14 May 2023 / Revised: 16 June 2023 / Accepted: 20 June 2023 / Published: 23 June 2023
(This article belongs to the Special Issue Molecular Biology of Colorectal Cancers)

Abstract

:

Simple Summary

Colorectal cancer (CRC) is one of the most prevalent cancers, and approximately a quarter of patients diagnosed at stage II exhibit a significant risk of recurrence. In this study, we successfully identified a microRNA (miRNA) signature allowing the recognition of patients at high recurrence risk. The validity of these findings has been confirmed through an entirely separate group of patients diagnosed with stage II microsatellite stability (MSS) colon adenocarcinoma (COAD). Most of the miRNAs present in the signature have demonstrated prognostic relevance in various other cancer types. Upon examining their gene targets, we discovered that some of these miRNAs are intricately involved in pivotal pathways of cancer progression.

Abstract

We aimed to identify and validate a set of miRNAs that could serve as a prognostic signature useful to determine the recurrence risk for patients with COAD. Small RNAs from tumors of 100 stage II, untreated, MSS colon cancer patients were sequenced for the discovery step. For this purpose, we built an miRNA score using an elastic net Cox regression model based on the disease-free survival status. Patients were grouped into high or low recurrence risk categories based on the median value of the score. We then validated these results in an independent sample of stage II microsatellite stable tumor tissues, with a hazard ratio of 3.24, (CI95% = 1.05–10.0) and a 10-year area under the receiver operating characteristic curve of 0.67. Functional analysis of the miRNAs present in the signature identified key pathways in cancer progression. In conclusion, the proposed signature of 12 miRNAs can contribute to improving the prediction of disease relapse in patients with stage II MSS colorectal cancer, and might be useful in deciding which patients may benefit from adjuvant chemotherapy.

1. Introduction

CRC is the third most newly diagnosed cancer type worldwide. Although systematic screening programs have reduced the incidence of CRC in Western countries [1], it is still the second leading cause of cancer-related deaths worldwide for both men and women [2], with millions of cases being reported each year. Currently, stage at diagnosis is the most relevant predictor of prognosis. It is known that about a quarter of patients with CRC are diagnosed in stage II, with localized disease and no evidence of regional lymph node invasion [3]. Nevertheless, disease will recur or progress to distant metastasis in about 20–25% of these patients. Clinical and pathological risk factors, such as the size and location of the tumor, are used to identify patients at high risk of recurrence, but they are not always reliable. As a result, there has been growing interest in the use of biomarkers to improve the accuracy of prognostic and predictive testing for CRC patients [4,5,6,7,8]. However, these approaches are limited by small populations or accuracy [9,10]. Molecular biomarkers could be used to identify patients who are at high risk of disease recurrence or improve stratification of patients who could benefit from adjuvant chemotherapy or immunotherapy.
miRNAs are short, double-stranded, non-coding RNA molecules, typically between 19 and 24 nucleotides in length, which play a critical role in the regulation of gene expression. MiRNAs are involved in post-transcriptional regulation of multiple protein coding genes, mainly by binding to the 3’ untranslated regions (UTRs) of target genes, leading to inhibition of messenger RNA (mRNA) transcription [11]. Changes in miRNA expression affect target genes regulation, and consequently, their deregulation can lead to irregular cell processes related to tumor development and progression [12]. To date, several miRNAs have been proposed to be either oncogenic or tumor suppressors [13,14,15], and there have been some miRNA signatures proposed as molecular biomarkers in CRC for both diagnosis and prognosis, as well as treatment decisions [16,17,18,19].
In this study, we aimed to identify and validate a signature of miRNAs with prognostic value in stage II COAD patients. We used next-generation sequencing (NGS) techniques to obtain miRNA expression values for a set of tumor samples, and we tried to validate the findings in an independent sample series.

2. Materials and Methods

2.1. Subjects and Samples

In the discovery series, we included Colonomics (CLX): 98 tumor tissue samples, MSS stage II patients with a new diagnosis of COAD at the University Hospital of Bellvitge in Barcelona (Spain) between January 1996 and December 2000. Patients were selected from those that had donated fresh tissue to the biobank and had undergone a complete surgical resection of the tumor, but had not received adjuvant chemotherapy. In addition, a minimum of 3 years of follow-up was required.
The validation series included public independent samples of 130 COAD patients (stage II) from The Cancer Genome Atlas (TCGA) study.
The study was performed in accordance with relevant ethics guidelines and regulations. The Clinical Research Ethics Committee of the Bellvitge Hospital approved the study protocol (PR178/11). Individuals provided written informed consent to participate and for genetic analysis to be carried out on their samples. Additional information about the study can be found at www.colonomics.org (accessed on 15 May 2023). This study carefully follows the recommendations for reporting proposed by the REMARK guidelines [20].

2.2. Sample Processing

Tumor samples were cut by the pathologist from the surgical specimen during the first hour after removal and kept frozen at −80 °C in the hospital’s tumor bank. Total RNA was isolated from tissue samples using the miRCURYTM RNA isolation kit (Exiqon, Vedbaek, Denmark) according to the manufacturer’s protocol, quantified using a NanoDrop® ND-1000 Spectrophotometer (Nanodrop technologies, Wilmington, DE, USA) and stored at −80 °C. The quality of these RNA samples was assessed with the RNA 6000 Nano Assay (Agilent Technologies, Santa Clara, CA, USA) following the manufacturer’s recommendations. The RNA integrity number (RIN) showed high quality values for all the samples (mean = 7.89, sd = 0.86). The RNA purity was measured with the ratio of absorbance at 260 nm and 280 nm (mean = 1.96, sd = 0.04). The quality control for the small RNA fraction was assessed with the Small RNA Assay in the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) following the manufacturer’s recommendations.

2.3. Small RNA-seq Analysis of the Discovery Series

The small RNA-seq was performed through the SOLiD platform. The PureLink miRNA isolation kit was used to construct the libraries of compatible fragments with SOLiD from an enriched fraction of small RNA. Sequencing microspheres were obtained by applying an emulsion PCR into an equimolar mixture of 48 libraries followed by an enrichment process before charging in the reaction chamber. Finally, the reaction to obtain the sequences (35 nucleotides + 10 nucleotides barcode) from the small RNA fraction was performed with the Applied Biosystems SOLiD 4 System. Samples were randomly distributed among the different sequencing slides to minimize batch effects. The data quality was estimated using the SOLiD Experimental Tracking System (SETS) software.

2.4. Expression Data of the Discovery Series

CLX is a multiomics experiment design with different high-throughput sequencing data. In addition to small RNA-seq, it has microarray expression data for the same subjects. Sample processing, quality control and normalization are described elsewhere [21].

2.5. Bioinformatics Analysis

Quality control of sequenced reads was ensured using specially designed bioinformatics framework for the SOLiD system [22]. Total number of reads, proportion of miscalled reads and low-average-quality-score proportion reads were evaluated. All samples passed the quality control criteria and were selected for further analysis. Next, quantification of specific miRNAs was performed by mapping reads to the reference of mature miRNA sequences annotated in miRBase release 22 [23], containing 2641 human mature miRNA sequences. The FASTX-toolkit [24] was used to preprocess miRNA data and provide compatible sequences for mapping with Bowtie aligner. Read adapters were trimmed with cutadapt [25], and finally, a table of counts was generated with SAMtools [26]. A principal component analysis (PCA) was computed to detect possible outliers. A filter based on low variability of miRNAs across all samples was performed to remove unwanted noisy data (standard deviation < 0.1). Data normalization was performed using DESeq2 package [27], after which it was transformed with a logarithmic function to reduce positive skewness. MiRNAs with normalized expression values not detected in more than 90% of samples were filtered out due to low expression. This criterion was mandatory for both discovery and validation datasets.

2.6. Statistical Analysis of Prognosis

For this study, disease-free survival (DFS) was assessed, and disease progression, defined as local tumor recurrence, metastasis or cancer-related death, was the event of interest. First, we wanted to inspect miRNA profiles based on possible sources of confounding variables. For this purpose, differential miRNA expression analysis (DEA) was carried out with the DESeq2 package for sex and tumor site. In addition, a proportional hazards assumption test was performed to assess possible sources of analysis bias from common covariates: age, sex and tumor site (left or right colon). Univariate Cox proportional hazard models were computed for each miRNA and adjusted for age, sex, tumor site and sub-stage. Kaplan–Meier survival curves were used to graph the results, which were split by median normalized expression values. Next, in order to identify an miRNA signature that could capture disease progression, a regularized Cox regression model was performed with an alpha parameter (α = 0.5) and leave-one-out cross-validation. Adjustment variables were included in all models, and we selected the model that minimized the cross-validation error. The coefficients for the miRNAs that were not shrunk in the optimal model were used to compute an miRNA prognosis risk score (RS) as follows:
i = 1 n = β i   e x p r i
where n is the total number of active miRNAs in the model, βi is the coefficient for each active miRNA, and expri is the expression value of each active miRNA. The score obtained was ranked and split into two equal groups by the median value. The performance of the model was assessed with Cox proportional hazard models (selecting miRNA (RS) and adjustment variables, as mentioned before) and Kaplan–Meier survival curves. Univariate Cox proportional hazard models and Kaplan–Meier estimates were obtained with the survival R package [28]. Regularized Cox regression models were computed with the glmnet R package [29]. All statistical analysis was performed in R version 3.5 [30].

2.7. Validation Analysis

For the validation dataset, TCGA COAD samples were filtered in order to obtain similar clinical characteristics and reproducible analysis. Sample exclusion criteria included no clinical information for disease-free survival status or microsatellite instability (MSI) subtype, the latter of which was assessed in cBioPortal [31,32]. Normalization of miRNA expression values and filtering was identical in discovery and validation series. Independent prognosis analysis for stage II was assessed. The same coefficients and cutoffs obtained in the training dataset were used for the validation.

2.8. Functional Characterization

We conducted two separate approaches to characterize the resulting miRNA signature and score. On the one hand, for the signature study, two different miRNA–mRNA interaction resources were interrogated to extract high-confidence gene targets for each miRNA present in our signature. Only common interactions present in mirDB (release 6.0) [33,34] and miRTarBase (release 9.0) [35] were included to define the functional role of miRNAs associated with prognosis. mirDB annotates predicted miRNA–target interactions (MTI) while miRTarBase captures experimentally validated MTIs from research articles. Network analysis was carried out with igraph [36] to identify hub miRNAs; next, we performed an enrichment analysis with the ReactomePA [37] R package based on the REACTOME pathway database, including direct targets for each selected miRNA. On the other hand, we wanted to study the relationship between miRNA prognosis RS and the abundance of tissue-infiltrating immune cell populations that could potentially play different roles in the tumor microenvironment. For this purpose, a deconvolution method [38] was used to estimate cell population proportions (immune and non-immune stromal) from average gene expression signals. Non-parametric Spearman’s rank-order correlations were computed to evaluate correlation patterns between each of the ten different cell types and miRNA prognosis RS.

3. Results

3.1. Study Population Characteristics and Quality Control of Samples

The CLX small-RNA sequencing dataset comprised 100 tumor tissues, stage II MSS, and 2641 different miRNAs were initially identified. The first filter removed 153 miRNAs due to low variability. The second filter was applied to remove low-expression features. A total of 928 miRNAs passed the filtering criteria. Two samples were filtered out after PCA analysis (see Supplementary Figure S1). The same procedure was performed for the validation dataset, which lowered the number of different miRNAs from 2117 to 796. Finally, 605 miRNAs were determined to be present in both datasets and were selected for the next analysis. Table 1 summarizes the main characteristics of the patients included in both datasets. DEA comparing tumor location resulted in only 10 significant miRNAs with different profiles between the left and right sides (seven overexpressed on the right side, three overexpressed on the left side). All of them had an absolute log2-fold change greater than 0.5. When comparing expression profiles between sex, only two significant miRNAs were observed (Supplementary Table S1 and Figure S2). None of the inspected covariates violated the proportional hazards assumption of the Cox model (p-values > 0.05, Supplementary Table S2). Univariate Cox proportional hazard models, adjusted by sex, age, site and stage, identified 55 miRNAs associated with disease-free survival (p-value < 0.05). However, none of them passed the false discovery rate (Benjamini–Hochberg-adjusted p-value < 0.05) (Supplementary Table S3).

3.2. miRNA Signature and Score

The regularized Cox regression model resulted in a 12-miRNA signature computed in the discovery dataset (Table 2). Of note, all miRNAs present in the signature were statistically significant according to the univariate models (p < 0.05), and all of them showed the same trend at the individual coefficient signs, suggesting low collinearity between all of them (Supplementary Table S4), which was confirmed with a Spearman’s correlation matrix (absolute Spearman’s r ≤ 0.43 for all pairwise comparisons (Supplementary Figure S3).
The RS formula obtained was:
miRNA RS = hsa-miR-1185-5p × (−0.185) + hsa-miR-16-5p × (−0.111) + hsa-miR-181a-2-3p × 0.181 + hsa-miR-204-5p × 0.003 + hsa-miR-2355-3p × 0.242 + hsa-miR-29b-2-5p × (−0.306) + hsa-miR-331-3p × 0.153 + hsa-miR-423-3p × (−0.355) + hsa-miR-432-5p × (−0.187) + hsa-miR-497-5p × (−0.183) + hsa-miR-656-3p × (−0.526) + hsa-miR-935 × (−0.136)
Kaplan–Meier curves demonstrated a good performance of the model, clearly differentiating low- and high-risk patient log-rank p-values = 1.62 × 10−6. Patients in the high-risk group were found to have higher recurrence rates (HR = 33.59, 4.34–244.8, p < 0.001) (Figure 1a,c). Three-, five- and ten-year disease-free survival (DFS) were selected to compute the area under the ROC curve (AUC). AUCs of 0.89, 0.92 and 0.94 were obtained for three-, five- and ten-year DFS, respectively (Figure 2a).

3.3. Validation

CLX miRNA signatures were tested on an independent COAD dataset from TCGA. Overall, the miRNA predictive capacity was poor, which is mainly explained by the poor performance on stage II MSI tumor samples. However, it improved substantially when only stage II MSS samples were analyzed (n = 100). A similar trend to that in CLX was observed; the miRNA risk group had an HR of 3.24, with a range of 1.05–10.0, and p = 0.041 (Figure 1d), and the low- versus high-risk patient log-rank p-value was 0.034, as assessed using KM curves (Figure 1b). The AUCs were 0.60, 0.59 and 0.67 for three-, five- and ten-year DFS, respectively (Figure 2b).

3.4. Functional Characterization

Gene target candidates for each miRNA present in our signature were retrieved from mirDB and mirTarBase (10,352 and 3673, respectively). Overall, common MTIs from both datasets revealed 1015 MTIs. Network analysis showed three hub miRNAs: hsa-miR-16-5p (483 gene interactions), hsa-miR-497-5p (252 gene interactions) and hsa-miR-204-5p (119 gene interactions). Of note, the first two miRNAs shared 242 common gene interactions (Figure 3). Next, Reactome pathway analysis identified relevant cancer pathways associated with both hsa-miR-16-5p and hsa-miR -497-5p, such as signaling by the TGF-beta receptor complex (p = 1.84 × 10−7, p = 1.76 × 10−5, respectively), regulation of RUNX1 expression and activity (p = 2.82 × 10−7, 1.47 × 10−7, respectively) and aberrant regulation of the mitotic G1/S transition in cancer due to RB1 defects (p = 2.82 × 10−7, 1.46 × 10−7, respectively). The complete results of the enrichment analysis are available in Supplementary Tables S5–S13.
MCP-counter was applied to gene expression data from the discovery series. Abundances of reported immune-infiltrating cell populations, as well as other non-immune cell types, are summarized in Supplementary Figure S4. B-cell infiltrates (Spearman’s r = −0.34, p-value = 6.24 × 10−4), myeloid dendritic cell infiltrates (Spearman’s r = −0.32 p-value = 1.51 × 10−3) and T-cell infiltrates (Spearman’s r = −0.29 p-value = 3.79 × 10−3) appeared with a moderate inverse association with miRNA RS (Supplementary Figure S5). It is worth mentioning that all tested cell population abundances were inversely correlated with miRNA RS, suggesting an overall increase in immune infiltration in tumors with lower values for the computed RS.

4. Discussion

A comprehensive analysis has been conducted in order to identify an miRNA signature with potential value for stratifying patients into different disease progression risk groups in stage II MSS colon cancer.
Similar to other prognostic signatures previously published [16,17,18,39], our signature demonstrated a significant association. However, its ability to accurately predict which patients will experience recurrence was limited [40]. Nevertheless, this does not imply that the signature is not valuable since it can be used to classify patients into distinct risk groups.
To the best of our knowledge, only 2 of the 12 miRNAs included in the proposed signature have been previously included in other CRC miRNA signatures [41,42]. However, several of them have been associated with CRC development and/or prognosis; some of them have also been associated with tumoral progression in other tissues. Our findings agree with current knowledge concerning the miRNAs present in our signature. Recently, it has been found that overexpression of hsa-miR-16-5p can inhibit CRC cell proliferation, migration, immune modulation and invasion [43,44]. Another recent publication pointed out the relationship between the down-regulation of hsa-miR-16-5p and hsa-miR-497-5p and the progression of endometrial cancer mediated by circular RNA hsa-circ-0011324 [45]. Both miRNAs appeared to have lower expression values in the CLX series high-risk group, as reported in univariate Cox proportional hazard models. hsa-miR-656-3p was also included in one miRNA classifier for tumor recurrence in stage II CRC [42]. Interestingly, it has also been identified as an inhibitor of CRC cell migration in vitro [46]. In contrast, hsa-miR-204-5p has been identified to be negatively associated with CRC progression and chemoresistance [47,48]. This specific miRNA goes in the opposite direction both in the discovery and validation series. Another miRNA present in the signature, hsa-miR-935, has been studied in different cancer types and seems to have different behaviors depending on the targeted tissue, inhibiting or promoting tumor development in glioblastoma, liver and gastric cancer [49,50,51]. Hsa-miR-423-3p was seen to be down-regulated in hepatocellular carcinoma compared to healthy and liver cirrhosis samples [52]. In a recent study conducted on bladder urothelial carcinoma patients, hsa-miR-432-5p was found to be a good biomarker for diagnosis, being under-expressed in tumoral samples [53]. Regarding the remnant miRNAs, one study suggested opposite effects from what we have reported for hsa-miR-181a-2-3p in glioblastoma [54], and two studies found hsa-miR-331-3p to be under-expressed in prostate cancers [55] and CRC [56] compared to healthy groups. Higher levels of hsa-29b-2-5p expression were associated with the staging of esophageal and gastric cancer [57] in TCGA. However, this association was not observed in TCGA COAD stage II samples. No previous studies were found for hsa-miR-2355-3p and hsa-miR-1185-5p related to cancer.
To assess the complex functional interrelationships between the miRNAs and their putative target mRNAs, we analyzed a network of validated miRNA–target gene associations. Our results show a significant enrichment of genes involved in cellular processes relevant for cancer progression, such as cell cycle regulation, interleukin signaling and cell migration [58,59]. Since the information on validated target genes for miRNA is still incomplete, these results might be biased because more well-known miRNAs share more interactions with important cancer genes. However, this behavior is expected to be minimized with the release of the latest updated versions. In addition, to further evaluate the possible effects of the miRNA signature, the immune profile analysis of the miRNA risk groups revealed an overall increased abundance of all types of immune cell populations as measured through deconvolution analyses, suggesting a protective effect of immune infiltration in tumors [60]. Although the role of immune cell infiltrates in cancer progression is complex and context-dependent [61], it is thought that in early CRC stages, immune cells could help to control the growth and spread of cancer cells. However, as the tumor progresses, these immune cells can adapt a pro-tumorigenic role, promoting tumor growth and metastasis [62,63].
Besides the considerable efforts to generate biomarkers for risk assessment based on miRNA expression levels [16,17,18,42], other emerging tools are being investigated as potential prognostic factors for stage II CRC. One such approach involves examining molecular characteristics, including mutations or expression profile alterations in BRAF, KRAS or PIK3CA. These molecular features can aid in determining more targeted treatment strategies for patients [64,65,66]. Another emerging tool, circulating tumor DNA (ctDNA) analysis, has shown promise in identifying minimal residual disease [67], which results in a higher risk of recurrence for patients and may require closer surveillance or additional treatment options. Furthermore, the study of the gut microbiome is also a recent topic of interest for the assessment of both tumor onset and recurrent disease [68,69]. Lastly, exploring the tumor microenvironment [40,70], such as TILs, tumor-associated fribroblasts and stromal characteristics, and understanding the interactions between tumor cells and their microenvironment might provide valuable prognostic information.
This study has several limitations. The most important is the low number of events in both the discovery and the validation series, which is probably related to the good prognosis for early COAD diagnosis. In addition, our results may have been underestimated due to the lack of information regarding the administration of adjuvant chemotherapy or radiation therapy in the validation series. The low number of events together with the limited sample size reduces the statistical power for this kind of analysis; as a result, we saw a trend in the validation cohort with borderline statistical significance. Another limitation of the study is the lack of microsatellite instability (MSI) patients. The proportion of MSI in our hospital was low, around 8%; thus, the CLX series was only composed of MSS patients.
Overall, based on the reported results, this signature could be valuable to stratify MSS stage II COAD patients and identify those that require adjuvant chemotherapy. Additionally, the individual miRNA prognostic data provided in the discovery series contribute to increasing the knowledge on these markers in CRC.

5. Conclusions

In summary, we have identified a panel of 12 miRNAs that can be used to stratify prognosis in MSS stage II COAD. These miRNAs have been described to regulate a large list of genes involved in relevant cancer pathways, which reinforces the validity of the panel. Further studies with larger samples sizes are needed to improve our ability to classify patients with recurrence risk in a more general way.

Supplementary Materials

The following supporting information can be downloaded at: www.mdpi.com/xxx/s1.Supplementary Figure S1: PCA representation for the first two components of CLX normalized miRNA. Supplementary Figure S2: Volcano plots for DEA miRNA for tumor side (a) and sex (b). Supplementary Figure S3: Correlation plot of the 12 miRNA presents in the signature. Supplementary Figure S4: Histogram of ten cell populations abundance detected by MPC counter method. Supplementary Figure S5: Scatter plots for significant correlation between miRNA prognostic RS and B-cells, T-cells and myeloid Dendritic cells detected by MCP-counter. Supplementary Table S1: Differential miRNA expression analysis. Supplementary Table S2: Cox Proportional Hazards assumption test. Supplementary Table S3: Univariate Cox regression models for discovery series, adjusted by age, sex, site and stage. Supplementary Table S4: Univariate Cox regression models for validation series adjusted by age, sex and site. Supplementary Table S5: hsa-miR-16-5p target gene Enrichment. Supplementary Table S6: hsa-miR-432-5p target gene Enrichment. Supplementary Table S7: hsa-miR-497-5p target gene Enrichment. Supplementary Table S8: hsa-miR-432-5p target gene Enrichment. Supplementary Table S9: hsa-miR-181-a-2-3p target gene Enrichment. Supplementary Table S10: hsa-miR-204-5p target gene Enrichment. Supplementary Table S11: hsa-miR-2355-3p target gene Enrichment. Supplementary Table S12: hsa-miR-656-3p target gene Enrichment. Table S13: hsa-miR-935 target gene Enrichment.

Author Contributions

Conceptualization and methodology V.M., F.M.-N. and A.C.; formal analysis, F.M.-N., A.C., A.D.-V., A.S., D.C., X.S. (Xavier Solé) and E.G.; resources C.S., R.S., S.B. and X.S. (Xavier Sanjuan); writing—original draft preparation F.M.-N. and A.C.; writing—review and editing, F.M.-N., A.C., V.M. and R.S.-P. All authors have read and agreed to the published version of the manuscript.

Funding

This study received funding from Instituto de Salud Carlos III and was co-funded by FEDER funds—a way to build Europe—under the grant PI17-00092; the Spanish Association Against Cancer (AECC) Scientific Foundation under the grant GCTRA18022MORE; Consortium for Biomedical Research in Epidemiology and Public Health (CIBERESP), action Genrisk. ADV received funding from the Catalan Government, under fellowship PERIS SLT017/20/00042.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Clinical Research Ethics Committee of Bellvitge Hospital (PR178/11).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Full normalized miRNA expression data matrix of the training set and its associated clinical data are available at the project site: https://www.colonomics.org/data (accessed on 15 May 2023). Raw small RNA seq data are available upon request at the European Genome-Phenome Archive under the study ID: EGAS00001002453; dataset ID: EGAD00001004827.

Acknowledgments

The authors would like to thank Isabel Padrol and Carmen Atencia for their technical assistance. The “Xarxa de Bancs de Tumors de Catalunya (XBTX)” sponsored by “Pla Director d’Oncologia de Catalunya” and the ICO Biobank and PLATAFORMA BIOBANCOS PT13/0010/0013 helped with simple collection.

Conflicts of Interest

VM, co-investigator, received grants from Aniling SL.

References

  1. Navarro, M.; Nicolas, A.; Ferrandez, A.; Lanas, A. Colorectal Cancer Population Screening Programs Worldwide in 2016: An Update. World J. Gastroenterol. 2017, 23, 3632. [Google Scholar] [CrossRef] [PubMed]
  2. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA. Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  3. AJCC Cancer Staging Manual|SpringerLink. Available online: https://link.springer.com/book/9783319406176 (accessed on 4 January 2023).
  4. Guo, W.; Cai, Y.; Liu, X.; Ji, Y.; Zhang, C.; Wang, L.; Liao, W.; Liu, Y.; Cui, N.; Xiang, J.; et al. Single-Exosome Profiling Identifies ITGB3+ and ITGAM+ Exosome Subpopulations as Promising Early Diagnostic Biomarkers and Therapeutic Targets for Colorectal Cancer. Research 2023, 6, 0041. [Google Scholar] [CrossRef] [PubMed]
  5. Levy, J.J.; Zavras, J.P.; Veziroglu, E.M.; Nasir-Moin, M.; Kolling, F.W.; Christensen, B.C.; Salas, L.A.; Barney, R.E.; Palisoul, S.M.; Ren, B.; et al. Identification of Spatial Proteomic Signatures of Colon Tumor Metastasis: A Digital Spatial Profiling Approach. Am. J. Pathol. 2023, 193, 778–795. [Google Scholar] [CrossRef]
  6. Ciocan, A.; Ciocan, R.A.; Al Hajjar, N.; Benea, A.M.; Pandrea, S.L.; Cătană, C.S.; Drugan, C.; Oprea, V.C.; Dîrzu, D.S.; Bolboacă, S.D. Exploratory Evaluation of Neopterin and Chitotriosidase as Potential Circulating Biomarkers for Colorectal Cancer. Biomedicines 2023, 11, 894. [Google Scholar] [CrossRef] [PubMed]
  7. Linke, C.; Hunger, R.; Reinwald, M.; Deckert, M.; Mantke, R. Quantification of Mitochondrial CfDNA Reveals New Perspectives for Early Diagnosis of Colorectal Cancer. BMC Cancer 2023, 23, 291. [Google Scholar] [CrossRef]
  8. Hosseini, F.A.; Rejali, L.; Zabihi, M.R.; Salehi, Z.; Daskar-Abkenar, E.; Taraz, T.; Fatemi, N.; Hashemi, M.; Asadzadeh-Aghdaei, H.; Nazemalhosseini-Mojarad, E. Long Non-coding RNA LINC00460 Contributes as a Potential Prognostic Biomarker through Its Oncogenic Role with ANXA2 in Colorectal Polyps. Mol. Biol. Rep. 2023, 50, 4505–4515. [Google Scholar] [CrossRef]
  9. Cortes-Ciriano, I.; Lee, S.; Park, W.-Y.; Kim, T.-M.; Park, P.J. A Molecular Portrait of Microsatellite Instability across Multiple Cancers. Nat. Commun. 2017, 8, 15180. [Google Scholar] [CrossRef] [Green Version]
  10. Gibney, G.T.; Weiner, L.M.; Atkins, M.B. Predictive Biomarkers for Checkpoint Inhibitor-Based Immunotherapy. Lancet Oncol. 2016, 17, e542–e551. [Google Scholar] [CrossRef] [Green Version]
  11. Bartel, D.P. MicroRNAs: Genomics, Biogenesis, Mechanism, and Function. Cell 2004, 116, 281–297. [Google Scholar] [CrossRef] [Green Version]
  12. MacFarlane, L.-A.; Murphy, P.R. MicroRNA: Biogenesis, Function and Role in Cancer. Curr. Genom. 2010, 11, 537–561. [Google Scholar] [CrossRef] [Green Version]
  13. Yu, W.; Liang, X.; Li, X.; Zhang, Y.; Sun, Z.; Liu, Y.; Wang, J. MicroRNA-195: A Review of Its Role in Cancers. OncoTargets Ther. 2018, 11, 7109–7123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Chen, D.; Sun, Y.; Yuan, Y.; Han, Z.; Zhang, P.; Zhang, J.; You, M.J.; Teruya-Feldstein, J.; Wang, M.; Gupta, S.; et al. MiR-100 Induces Epithelial-Mesenchymal Transition but Suppresses Tumorigenesis, Migration and Invasion. PLoS Genet. 2014, 10, e1004177. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, X.; Zhang, H.; Yang, H.; Bai, M.; Ning, T.; Deng, T.; Liu, R.; Fan, Q.; Zhu, K.; Li, J.; et al. Exosome-Delivered CircRNA Promotes Glycolysis to Induce Chemoresistance through the MiR-122-PKM2 Axis in Colorectal Cancer. Mol. Oncol. 2020, 14, 539–555. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, Y.; Huang, L.; Shan, N.; Ma, H.; Lu, S.; Chen, X.; Long, H. Establishing a Three-MiRNA Signature as a Prognostic Model for Colorectal Cancer through Bioinformatics Analysis. Aging 2021, 13, 19894–19907. [Google Scholar] [CrossRef] [PubMed]
  17. Ma, R.; Zhao, Y.; He, M.; Zhao, H.; Zhang, Y.; Zhou, S.; Gao, M.; Di, D.; Wang, J.; Ding, J.; et al. Identifying a Ten-MicroRNA Signature as a Superior Prognosis Biomarker in Colon Adenocarcinoma. Cancer Cell Int. 2019, 19, 360. [Google Scholar] [CrossRef] [Green Version]
  18. Zhang, J.-X.; Song, W.; Chen, Z.-H.; Wei, J.-H.; Liao, Y.-J.; Lei, J.; Hu, M.; Chen, G.-Z.; Liao, B.; Lu, J.; et al. Prognostic and Predictive Value of a MicroRNA Signature in Stage II Colon Cancer: A MicroRNA Expression Analysis. Lancet Oncol. 2013, 14, 1295–1306. [Google Scholar] [CrossRef]
  19. Slattery, M.L.; Herrick, J.S.; Mullany, L.E.; Valeri, N.; Stevens, J.; Caan, B.J.; Samowitz, W.; Wolff, R.K. An Evaluation and Replication of MiRNAs with Disease Stage and Colorectal Cancer-Specific Mortality. Int. J. Cancer 2015, 137, 428–438. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. McShane, L.M.; Altman, D.G.; Sauerbrei, W.; Taube, S.E.; Gion, M.; Clark, G.M.; Statistics Subcommittee of the NCI-EORTC Working Group on Cancer Diagnostics. REporting Recommendations for Tumour MARKer Prognostic Studies (REMARK). Br. J. Cancer 2005, 93, 387–391. [Google Scholar] [CrossRef] [Green Version]
  21. Díez-Villanueva, A.; Sanz-Pamplona, R.; Solé, X.; Cordero, D.; Crous-Bou, M.; Guinó, E.; Lopez-Doriga, A.; Berenguer, A.; Aussó, S.; Paré-Brunet, L.; et al. COLONOMICS—Integrative Omics Data of One Hundred Paired Normal-Tumoral Samples from Colon Cancer Patients. Sci. Data 2022, 9, 595. [Google Scholar] [CrossRef]
  22. Sasson, A.; Michael, T.P. Filtering Error from SOLiD Output. Bioinformatics 2010, 26, 849–850. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kozomara, A.; Birgaoanu, M.; Griffiths-Jones, S. MiRBase: From MicroRNA Sequences to Function. Nucleic Acids Res. 2019, 47, D155–D162. [Google Scholar] [CrossRef]
  24. FASTX-Toolkit. Available online: http://hannonlab.cshl.edu/fastx_toolkit/ (accessed on 11 April 2023).
  25. Martin, M. Cutadapt Removes Adapter Sequences from High-Throughput Sequencing Reads. EMBnet. J. 2011, 17, 10. [Google Scholar] [CrossRef]
  26. Danecek, P.; Bonfield, J.K.; Liddle, J.; Marshall, J.; Ohan, V.; Pollard, M.O.; Whitwham, A.; Keane, T.; McCarthy, S.A.; Davies, R.M.; et al. Twelve Years of SAMtools and BCFtools. GigaScience 2021, 10, giab008. [Google Scholar] [CrossRef] [PubMed]
  27. Love, M.I.; Huber, W.; Anders, S. Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Therneau, T. A Package for Survival Analysis in R. Available online: https://cran.r-project.org/web/packages/survival/vignettes/survival.pdf (accessed on 20 April 2023).
  29. Simon, N.; Friedman, J.H.; Hastie, T.; Tibshirani, R. Regularization Paths for Cox’s Proportional Hazards Model via Coordinate Descent. J. Stat. Softw. 2011, 39, 1–13. [Google Scholar] [CrossRef] [PubMed]
  30. R core team. R: The R Project for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 25 April 2023).
  31. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The CBio Cancer Genomics Portal: An Open Platform for Exploring Multidimensional Cancer Genomics Data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [Green Version]
  32. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the CBioPortal. Sci. Signal. 2013, 6, pl1. [Google Scholar] [CrossRef] [Green Version]
  33. Chen, Y.; Wang, X. MiRDB: An Online Database for Prediction of Functional MicroRNA Targets. Nucleic Acids Res. 2020, 48, D127–D131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Liu, W.; Wang, X. Prediction of Functional MicroRNA Targets by Integrative Modeling of MicroRNA Binding and Target Expression Data. Genome Biol. 2019, 20, 18. [Google Scholar] [CrossRef]
  35. Huang, H.-Y.; Lin, Y.-C.-D.; Cui, S.; Huang, Y.; Tang, Y.; Xu, J.; Bao, J.; Li, Y.; Wen, J.; Zuo, H.; et al. MiRTarBase Update 2022: An Informative Resource for Experimentally Validated MiRNA-Target Interactions. Nucleic Acids Res. 2022, 50, D222–D230. [Google Scholar] [CrossRef]
  36. Csardi, G.; Nepusz, T. The Igraph Software Package for Complex Network Research. InterJ. Complex Syst. 2005, 1695, 1–9. [Google Scholar]
  37. Yu, G.; He, Q.-Y. ReactomePA: An R/Bioconductor Package for Reactome Pathway Analysis and Visualization. Mol. Biosyst. 2016, 12, 477–479. [Google Scholar] [CrossRef]
  38. Becht, E.; Giraldo, N.A.; Lacroix, L.; Buttard, B.; Elarouci, N.; Petitprez, F.; Selves, J.; Laurent-Puig, P.; Sautès-Fridman, C.; Fridman, W.H.; et al. Estimating the Population Abundance of Tissue-Infiltrating Immune and Stromal Cell Populations Using Gene Expression. Genome Biol. 2016, 17, 218. [Google Scholar] [CrossRef] [PubMed]
  39. Jiang, S.; Xie, X.; Jiang, H. Establishment of a 7-MicroRNA Prognostic Signature and Identification of Hub Target Genes in Colorectal Carcinoma. Transl. Cancer Res. 2022, 11, 367. [Google Scholar] [CrossRef]
  40. Sanz-Pamplona, R.; Berenguer, A.; Cordero, D.; Riccadonna, S.; Solé, X.; Crous-Bou, M.; Guinó, E.; Sanjuan, X.; Biondo, S.; Soriano, A.; et al. Clinical Value of Prognosis Gene Expression Signatures in Colorectal Cancer: A Systematic Review. PLoS ONE 2012, 7, e48877. [Google Scholar] [CrossRef] [Green Version]
  41. Yang, Z.; Lu, S.; Wang, Y.; Tang, H.; Wang, B.; Sun, X.; Qu, J.; Rao, B. A Novel Defined Necroptosis-Related MiRNAs Signature for Predicting the Prognosis of Colon Cancer. Int. J. Gen. Med. 2022, 15, 555–565. [Google Scholar] [CrossRef]
  42. Jacob, H.; Stanisavljevic, L.; Storli, K.E.; Hestetun, K.E.; Dahl, O.; Myklebust, M.P. A Four-MicroRNA Classifier as a Novel Prognostic Marker for Tumor Recurrence in Stage II Colon Cancer. Sci. Rep. 2018, 8, 6157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Huang, X.; Xu, X.; Ke, H.; Pan, X.; Ai, J.; Xie, R.; Lan, G.; Hu, Y.; Wu, Y. MicroRNA-16-5p Suppresses Cell Proliferation and Angiogenesis in Colorectal Cancer by Negatively Regulating Forkhead Box K1 to Block the PI3K/Akt/MTOR Pathway. Eur. J. Histochem. EJH 2022, 66, 3333. [Google Scholar] [CrossRef] [PubMed]
  44. Zhu, T.; Lin, Z.; Han, S.; Wei, Y.; Lu, G.; Zhang, Y.; Xiao, W.; Wang, Z.; Jia, X.; Gong, W. Low MiR-16 Expression Induces Regulatory CD4+NKG2D+ T Cells Involved in Colorectal Cancer Progression. Am. J. Cancer Res. 2021, 11, 1540–1556. [Google Scholar] [PubMed]
  45. Liu, D.; Bi, X.; Yang, Y. Circular RNA Hsa_circ_0011324 Is Involved in Endometrial Cancer Progression and the Evolution of Its Mechanism. Bioengineered 2022, 13, 7485–7499. [Google Scholar] [CrossRef] [PubMed]
  46. Zhang, B.; Gao, S.; Bao, Z.; Pan, C.; Tian, Q.; Tang, Q. MicroRNA-656-3p Inhibits Colorectal Cancer Cell Migration, Invasion, and Chemo-Resistance by Targeting Sphingosine-1-Phosphate Phosphatase 1. Bioengineered 2022, 13, 3810–3826. [Google Scholar] [CrossRef] [PubMed]
  47. Yao, S.; Yin, Y.; Jin, G.; Li, D.; Li, M.; Hu, Y.; Feng, Y.; Liu, Y.; Bian, Z.; Wang, X.; et al. Exosome-Mediated Delivery of MiR-204-5p Inhibits Tumor Growth and Chemoresistance. Cancer Med. 2020, 9, 5989–5998. [Google Scholar] [CrossRef] [PubMed]
  48. Bian, Z.; Jin, L.; Zhang, J.; Yin, Y.; Quan, C.; Hu, Y.; Feng, Y.; Liu, H.; Fei, B.; Mao, Y.; et al. LncRNA—UCA1 Enhances Cell Proliferation and 5-Fluorouracil Resistance in Colorectal Cancer by Inhibiting MiR-204-5p. Sci. Rep. 2016, 6, 23892. [Google Scholar] [CrossRef] [Green Version]
  49. Huang, Y.; Xiao, W.; Jiang, X.; Li, H. MicroRNA-935 Acts as a Prognostic Marker and Promotes Cell Proliferation, Migration, and Invasion in Colorectal Cancer. Cancer Biomark. 2019, 26, 229–237. [Google Scholar] [CrossRef]
  50. Zhang, D.; Ma, S.; Zhang, C.; Li, P.; Mao, B.; Guan, X.; Zhou, W.; Peng, J.; Wang, X.; Li, S.; et al. MicroRNA-935 Directly Targets FZD6 to Inhibit the Proliferation of Human Glioblastoma and Correlate to Glioma Malignancy and Prognosis. Front. Oncol. 2021, 11, 566492. [Google Scholar] [CrossRef]
  51. Yang, M.; Cui, G.; Ding, M.; Yang, W.; Liu, Y.; Dai, D.; Chen, L. MiR-935 Promotes Gastric Cancer Cell Proliferation by Targeting SOX7. Biomed. Pharmacother. Biomed. Pharmacother. 2016, 79, 153–158. [Google Scholar] [CrossRef] [PubMed]
  52. de Oliveira, A.R.C.P.; Castanhole-Nunes, M.M.U.; Biselli-Chicote, P.M.; Pavarino, É.C.; da Silva, R.C.M.A.; da Silva, R.F.; Goloni-Bertollo, E.M. Differential Expression of Angiogenesis-Related MiRNAs and VEGFA in Cirrhosis and Hepatocellular Carcinoma. Arch. Med. Sci. AMS 2020, 16, 1150–1157. [Google Scholar] [CrossRef]
  53. Yerukala Sathipati, S.; Tsai, M.-J.; Shukla, S.K.; Ho, S.-Y.; Liu, Y.; Beheshti, A. MicroRNA Signature for Estimating the Survival Time in Patients with Bladder Urothelial Carcinoma. Sci. Rep. 2022, 12, 4141. [Google Scholar] [CrossRef] [PubMed]
  54. Plata-Bello, J.; Fariña-Jerónimo, H.; Betancor, I.; Salido, E. High Expression of FOXP2 Is Associated with Worse Prognosis in Glioblastoma. World Neurosurg. 2021, 150, e253–e278. [Google Scholar] [CrossRef]
  55. Luedemann, C.; Reinersmann, J.-L.; Klinger, C.; Degener, S.; Dreger, N.M.; Roth, S.; Kaufmann, M.; Savelsbergh, A. Prostate Cancer-Associated MiRNAs in Saliva: First Steps to an Easily Accessible and Reliable Screening Tool. Biomolecules 2022, 12, 1366. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, S.; Xiang, J.; Li, Z.; Lu, S.; Hu, J.; Gao, X.; Yu, L.; Wang, L.; Wang, J.; Wu, Y.; et al. A Plasma MicroRNA Panel for Early Detection of Colorectal Cancer. Int. J. Cancer 2015, 136, 152–161. [Google Scholar] [CrossRef]
  57. Chen, Z.; Shen, Z.; Zhang, Z.; Zhao, D.; Xu, L.; Zhang, L. RNA-Associated Co-Expression Network Identifies Novel Biomarkers for Digestive System Cancer. Front. Genet. 2021, 12, 659788. [Google Scholar] [CrossRef] [PubMed]
  58. Ruan, X.-J.; Ye, B.-L.; Zheng, Z.-H.; Li, S.-T.; Zheng, X.-F.; Zhang, S.-Z. TGFβ1I1 Suppressed Cell Migration and Invasion in Colorectal Cancer by Inhibiting the TGF-β Pathway and EMT Progress. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 7294–7302. [Google Scholar] [CrossRef]
  59. Amirkhah, R.; Schmitz, U.; Linnebacher, M.; Wolkenhauer, O.; Farazmand, A. MicroRNA-MRNA Interactions in Colorectal Cancer and Their Role in Tumor Progression. Genes Chromosomes Cancer 2015, 54, 129–141. [Google Scholar] [CrossRef] [PubMed]
  60. Wu, X.; Li, J.; Zhang, Y.; Cheng, Y.; Wu, Z.; Zhan, W.; Deng, Y. Identification of Immune Cell Infiltration Landscape for Predicting Prognosis of Colorectal Cancer. Gastroenterol. Rep. 2023, 11, goad014. [Google Scholar] [CrossRef] [PubMed]
  61. Xiong, Y.; Wang, K.; Zhou, H.; Peng, L.; You, W.; Fu, Z. Profiles of Immune Infiltration in Colorectal Cancer and Their Clinical Significant: A Gene Expression-Based Study. Cancer Med. 2018, 7, 4496–4508. [Google Scholar] [CrossRef] [Green Version]
  62. Fridman, W.H.; Pagès, F.; Sautès-Fridman, C.; Galon, J. The Immune Contexture in Human Tumours: Impact on Clinical Outcome. Nat. Rev. Cancer 2012, 12, 298–306. [Google Scholar] [CrossRef]
  63. Teng, M.W.L.; Ngiow, S.F.; Ribas, A.; Smyth, M.J. Classifying Cancers Based on T-Cell Infiltration and PD-L1. Cancer Res. 2015, 75, 2139–2145. [Google Scholar] [CrossRef] [Green Version]
  64. Potocki, P.M.; Wójcik, P.; Chmura, Ł.; Goc, B.; Fedewicz, M.; Bielańska, Z.; Swadźba, J.; Konopka, K.; Kwinta, Ł.; Wysocki, P.J. Clinical Characterization of Targetable Mutations (BRAF V600E and KRAS G12C) in Advanced Colorectal Cancer-A Nation-Wide Study. Int. J. Mol. Sci. 2023, 24, 9073. [Google Scholar] [CrossRef]
  65. Ma, Z.; Qi, Z.; Gu, C.; Yang, Z.; Ding, Y.; Zhou, Y.; Wang, W.; Zou, Q. BRAFV600E Mutation Promoted the Growth and Chemoresistance of Colorectal Cancer. Am. J. Cancer Res. 2023, 13, 1486–1497. [Google Scholar] [PubMed]
  66. Mirzapoor Abbasabadi, Z.; Hamedi Asl, D.; Rahmani, B.; Shahbadori, R.; Karami, S.; Peymani, A.; Taghizadeh, S.; Samiee Rad, F. KRAS, NRAS, BRAF, and PIK3CA Mutation Rates, Clinicopathological Association, and Their Prognostic Value in Iranian Colorectal Cancer Patients. J. Clin. Lab. Anal. 2023, 37, e24868. [Google Scholar] [CrossRef] [PubMed]
  67. Mo, S.; Ye, L.; Wang, D.; Han, L.; Zhou, S.; Wang, H.; Dai, W.; Wang, Y.; Luo, W.; Wang, R.; et al. Early Detection of Molecular Residual Disease and Risk Stratification for Stage I to III Colorectal Cancer via Circulating Tumor DNA Methylation. JAMA Oncol. 2023, 9, 770–778. [Google Scholar] [CrossRef] [PubMed]
  68. Debelius, J.W.; Engstrand, L.; Matussek, A.; Brusselaers, N.; Morton, J.T.; Stenmarker, M.; Olsen, R.S. The Local Tumor Microbiome Is Associated with Survival in Late-Stage Colorectal Cancer Patients. Microbiol. Spectr. 2023, 11, e0506622. [Google Scholar] [CrossRef] [PubMed]
  69. Xu, Y.; Zhao, J.; Ma, Y.; Liu, J.; Cui, Y.; Yuan, Y.; Xiang, C.; Ma, D.; Liu, H. The Microbiome Types of Colorectal Tissue Are Potentially Associated with the Prognosis of Patients with Colorectal Cancer. Front. Microbiol. 2023, 14, 1100873. [Google Scholar] [CrossRef]
  70. Kamal, Y.; Dwan, D.; Hoehn, H.J.; Sanz-Pamplona, R.; Alonso, M.H.; Moreno, V.; Cheng, C.; Schell, M.J.; Kim, Y.; Felder, S.I.; et al. Tumor Immune Infiltration Estimated from Gene Expression Profiles Predicts Colorectal Cancer Relapse. Oncoimmunology 2021, 10, 1862529. [Google Scholar] [CrossRef]
Figure 1. Comparison of the prognostic value with stratification analysis by miRNA risk group using Kaplan–Meier disease-free survival curves. (a) Discovery series (Colonomics). (b) Validation series (TCGA stage II MSS). Multivariate COX regression hazard models. (c) Discovery series (Colonomics). (d) Validation series (TCGA stage II MSS). Statistically significant HRs denoted with *.
Figure 1. Comparison of the prognostic value with stratification analysis by miRNA risk group using Kaplan–Meier disease-free survival curves. (a) Discovery series (Colonomics). (b) Validation series (TCGA stage II MSS). Multivariate COX regression hazard models. (c) Discovery series (Colonomics). (d) Validation series (TCGA stage II MSS). Statistically significant HRs denoted with *.
Cancers 15 03301 g001aCancers 15 03301 g001b
Figure 2. ROC curves at 3-, 5- and 10-year disease-free survival. (a) Discovery series (Colonomics). (b) Validation series (TCGA stage II MSS).
Figure 2. ROC curves at 3-, 5- and 10-year disease-free survival. (a) Discovery series (Colonomics). (b) Validation series (TCGA stage II MSS).
Cancers 15 03301 g002
Figure 3. mRNA–miRNA interaction network representation. Light-green nodes represent miRNAs present in the signature, and orange nodes represent their direct target mRNAs. Gray lines represent an interaction present between an miRNA and an mRNA (validated in miRTarBase and also predicted by mirDB). Node size is proportional to its degree centrality measure (number of direct interactors).
Figure 3. mRNA–miRNA interaction network representation. Light-green nodes represent miRNAs present in the signature, and orange nodes represent their direct target mRNAs. Gray lines represent an interaction present between an miRNA and an mRNA (validated in miRTarBase and also predicted by mirDB). Node size is proportional to its degree centrality measure (number of direct interactors).
Cancers 15 03301 g003
Table 1. Summary of characteristics of the patients included in the study.
Table 1. Summary of characteristics of the patients included in the study.
Colonomics n (%)TCGA n (%)
Number of Patients98130
Gender
Male70 (71.43%)69 (53.08%)
Female28 (28.57%)61 (46.92%)
Median Age (Years)7169
Tumor Site
Right38 (38.78%)75 (57.69%)
Left60 (61.22%)50 (38.46%)
Stage
II-A90 (91.84%)99 (76.15%)
II-B8 (8.16%)6 (4.62%)
Disease-Free Survival
No Event76 (77.55%)104 (80.00%)
Event22 (22.45%)26 (20.00%)
Microsatellite Instability
MSS98 (100%)101 (77.69%)
MSI0 (0%)20 (15.38%)
Median Metastatic Lymph Nodes0 (100%)0 (100%)
Median Isolated Lymph Nodes18.520.0
Lymphatic Invasion
Yes7 (0.07%)26 (20.00%)
No86 (87.76%)92 (70.77%)
Perineural Invasion
Yes2 (2.04%)13 (10.00%)
No83 (84.69%)38 (29.23%)
Table 2. List of miRNAs present in the signature and the coefficient extracted from the elastic net Cox regression model.
Table 2. List of miRNAs present in the signature and the coefficient extracted from the elastic net Cox regression model.
miRNACoefficient
hsa-miR-1185-5p−0.185
hsa-miR-16-5p−0.111
hsa-miR-181a-2-3p0.181
hsa-miR-204-5p0.003
hsa-miR-2355-3p0.242
hsa-miR-29b-2-5p−0.306
hsa-miR-331-3p0.153
hsa-miR-423-3p−0.355
hsa-miR-432-5p−0.187
hsa-miR-497-5p−0.183
hsa-miR-656-3p−0.526
hsa-miR-935−0.136
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Moratalla-Navarro, F.; Díez-Villanueva, A.; Garcia-Serrano, A.; Closa, A.; Cordero, D.; Solé, X.; Guinó, E.; Sanz-Pamplona, R.; Sanjuan, X.; Santos, C.; et al. Identification of a Twelve-microRNA Signature with Prognostic Value in Stage II Microsatellite Stable Colon Cancer. Cancers 2023, 15, 3301. https://doi.org/10.3390/cancers15133301

AMA Style

Moratalla-Navarro F, Díez-Villanueva A, Garcia-Serrano A, Closa A, Cordero D, Solé X, Guinó E, Sanz-Pamplona R, Sanjuan X, Santos C, et al. Identification of a Twelve-microRNA Signature with Prognostic Value in Stage II Microsatellite Stable Colon Cancer. Cancers. 2023; 15(13):3301. https://doi.org/10.3390/cancers15133301

Chicago/Turabian Style

Moratalla-Navarro, Ferran, Anna Díez-Villanueva, Ainhoa Garcia-Serrano, Adrià Closa, David Cordero, Xavier Solé, Elisabet Guinó, Rebeca Sanz-Pamplona, Xavier Sanjuan, Cristina Santos, and et al. 2023. "Identification of a Twelve-microRNA Signature with Prognostic Value in Stage II Microsatellite Stable Colon Cancer" Cancers 15, no. 13: 3301. https://doi.org/10.3390/cancers15133301

APA Style

Moratalla-Navarro, F., Díez-Villanueva, A., Garcia-Serrano, A., Closa, A., Cordero, D., Solé, X., Guinó, E., Sanz-Pamplona, R., Sanjuan, X., Santos, C., Biondo, S., Salazar, R., & Moreno, V. (2023). Identification of a Twelve-microRNA Signature with Prognostic Value in Stage II Microsatellite Stable Colon Cancer. Cancers, 15(13), 3301. https://doi.org/10.3390/cancers15133301

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop