Svoboda | Graniru | BBC Russia | Golosameriki | Facebook
Next Article in Journal
The Circulating miR-107 as a Potential Biomarker Up-Regulated in Castration-Resistant Prostate Cancer
Previous Article in Journal
Predicting the Effect of miRNA on Gene Regulation to Foster Translational Multi-Omics Research—A Review on the Role of Super-Enhancers
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
ncRNA
Volume 10
Issue 5
10.3390/ncrna10050046
ncrna-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MicroRNA Profiling as a Predictive Indicator for Time to First Treatment in Chronic Lymphocytic Leukemia: Insights from the O-CLL1 Prospective Study

by
Ennio Nano
1,†,
Francesco Reggiani
2,†,
Adriana Agnese Amaro
2,
Paola Monti
3,
Monica Colombo
1,
Nadia Bertola
1,
Fabiana Ferrero
1,4,
Franco Fais
1,4,
Antonella Bruzzese
5,
Enrica Antonia Martino
5,
Ernesto Vigna
5,
Noemi Puccio
6,7,
Mariaelena Pistoni
6,
Federica Torricelli
6,
Graziella D’Arrigo
8,
Gianluigi Greco
9,
Giovanni Tripepi
8,
Carlo Adornetto
9,
Massimo Gentile
5,10,
Manlio Ferrarini
4,
Massimo Negrini
11,
Fortunato Morabito
12,*,‡,
Antonino Neri
13,*,‡ and
Giovanna Cutrona
1,‡add Show full author list remove Hide full author list
1
Molecular Pathology Unit, IRCCS Ospedale Policlinico San Martino, 16132 Genoa, Italy
2
SSD Gene Expression Regulation, IRCCS Ospedale Policlinico San Martino, 16132 Genova, Italy
3
Mutagenesis and Cancer Prevention Unit, IRCCS Ospedale Policlinico San Martino, 16132 Genoa, Italy
4
Department of Experimental Medicine, University of Genoa, 16132 Genoa, Italy
5
Hematology Unit, Department of Onco-Hematology, Azienda Ospedaliera Annunziata, 87100 Cosenza, Italy
6
Laboratory of Translational Research, Azienda USL-IRCCS di Reggio Emilia, 42122 Reggio Emilia, Italy
7
Clinical and Experimental Medicine PhD Program, University of Modena and Reggio Emilia, 41121 Modena, Italy
8
Institute of Clinical Physiology (IFC-CNR), Section of Reggio Calabria, 89124 Reggio Calabria, Italy
9
Department of Mathematics and Computer Science, University of Calabria, 87100 Cosenza, Italy
10
Department of Pharmacy, Health and Nutritional Science, University of Calabria, 87036 Rende, Italy
11
Department of Translational Medicine, University of Ferrara, 44121 Ferrara, Italy
12
Gruppo Amici Dell’Ematologia Foundation-GrADE, 42122 Reggio Emilia, Italy
13
Scientific Directorate, Azienda USL-IRCCS di Reggio Emilia, 42122 Reggio Emilia, Italy
 
add Show full affiliation list
 
remove Hide full affiliation list
*
Authors to whom correspondence should be addressed.
†
These authors have contributed equally to this work as first name.
‡
These authors have contributed equally to this work as last name.
Non-Coding RNA 2024, 10(5), 46; https://doi.org/10.3390/ncrna10050046
Submission received: 25 July 2024 / Revised: 13 August 2024 / Accepted: 22 August 2024 / Published: 23 August 2024
Browse Figures
Review Reports Versions Notes

Abstract

:
A “watch and wait” strategy, delaying treatment until active disease manifests, is adopted for most CLL cases; however, prognostic models incorporating biomarkers have shown to be useful to predict treatment requirement. In our prospective O-CLL1 study including 224 patients, we investigated the predictive role of 513 microRNAs (miRNAs) on time to first treatment (TTFT). In the context of this study, six well-established variables (i.e., Rai stage, beta-2-microglobulin levels, IGVH mutational status, del11q, del17p, and NOTCH1 mutations) maintained significant associations with TTFT in a basic multivariable model, collectively yielding a Harrell’s C-index of 75% and explaining 45.4% of the variance in the prediction of TTFT. Concerning miRNAs, 73 out of 513 were significantly associated with TTFT in a univariable model; of these, 16 retained an independent relationship with the outcome in a multivariable analysis. For 8 of these (i.e., miR-582-3p, miR-33a-3p, miR-516a-5p, miR-99a-5p, and miR-296-3p, miR-502-5p, miR-625-5p, and miR-29c-3p), a lower expression correlated with a shorter TTFT, whereas in the remaining eight (i.e., miR-150-5p, miR-148a-3p, miR-28-5p, miR-144-5p, miR-671-5p, miR-1-3p, miR-193a-3p, and miR-124-3p), the higher expression was associated with shorter TTFT. Integrating these miRNAs into the basic model significantly enhanced predictive accuracy, raising the Harrell’s C-index to 81.1% and the explained variation in TTFT to 63.3%. Moreover, the inclusion of the miRNA scores enhanced the integrated discrimination improvement (IDI) and the net reclassification index (NRI), underscoring the potential of miRNAs to refine CLL prognostic models and providing insights for clinical decision-making. In silico analyses on the differently expressed miRNAs revealed their potential regulatory functions of several pathways, including those involved in the therapeutic responses. To add a biological context to the clinical evidence, an miRNA–mRNA correlation analysis revealed at least one significant negative correlation between 15 of the identified miRNAs and a set of 50 artificial intelligence (AI)-selected genes, previously identified by us as relevant for TTFT prediction in the same cohort of CLL patients. In conclusion, the identification of specific miRNAs as predictors of TTFT holds promise for enhancing risk stratification in CLL to predict therapeutic needs. However, further validation studies and in-depth functional analyses are required to confirm the robustness of these observations and to facilitate their translation into meaningful clinical utility.

1. Introduction

Chronic lymphocytic leukemia (CLL) is a B-cell disorder characterized by the monoclonal accumulation of CD5/CD23-positive lymphocytes at multiple sites [1,2,3] and a clinical heterogeneity [4,5] which has been correlated with different cytogenetic and molecular features of leukemic cells [6,7,8,9].
Many CLL patients do not require immediate therapeutic intervention, and for the majority, therapy can be initiated after a variable period of disease progression. To avoid useless treatment and to limit treatment only to the progressing cases, a “watch and wait” strategy is generally adopted in which only patients who progress, according to defined clinical criteria, are treated [1]. Due to the aforementioned variability of CLL, prognostic studies aimed at defining the risk of disease progression, warranting therapeutic intervention, have a relevant clinical impact [10,11].
While single prognostic factors, although relevant, may have a limited prognostic significance, the evaluation of multiple biomarkers with varying prognostic power may lead to more definite prognostic models [12]. To this end, prognostic scores have been developed that incorporate different biomarkers into comprehensive models. The CLL-International Prognostic Index (CLL-IPI) has emerged as one of the well-established models, demonstrating its efficacy in predicting overall survival, as well as time to first treatment (TTFT) and progression-free survival prediction in the context of chemo-immunotherapy [13,14]. Furthermore, our group has documented the predictive value of CLL-IPI for TTFT in early-stage CLL. In the same setting of Binet stage A CLL, Condoluci et al. introduced the International Prognostic Score for Early CLL patients (IPS-E) [15]. Further validated by Smolej et al. [16], the IPS-E and its alternative version (AIPS-E) effectively predicted TTFT and therefore may provide valuable guidance for clinical decision-making in early-stage CLL patients. In our prospective O-CLL1 study, evaluating newly diagnosed Binet stage A patients, both AIPS-E and IPS-E accurately predicted the need for therapy in early-stage CLL patients [17].
The advent of next-generation sequencing (NGS) has facilitated the discovery of an unprecedented number of gene mutations with potential prognostic and predictive value [9]. A recent study found that mutations were detected in roughly 35% of CLL patients, with frequencies ranging from 2.3% to 9.8% of cases, and that NOTCH1 mutations were the most common [18]. When the predictive impact of nine recurrently mutated genes was tested, eight of them (BIRC3, EGR2, NFKBIE, NOTCH1, POT1, SF3B1, TP53, and XPO1) were independently associated with shorter TTFT [18].
In recent years, the role of microRNAs (miRNAs) in CLL has garnered increasing attention (reviewed in [19,20,21]). MiRNAs are evolutionarily conserved, single-stranded non-coding RNA molecules that can redundantly and simultaneously regulate the expression of multiple genes. Acting post-transcriptionally, miRNAs bind to the messenger RNA of target genes, affecting translation or inducing degradation [22,23].
Several miRNAs have been described to be associated with CLL prognosis [24,25,26,27,28,29,30] and pathogenesis [31,32]. In addition, miRNAs may contribute to the deregulation of apoptosis [31,32,33], BCR signaling [34,35,36], or metabolism [37] in CLL cells, and treatment of CLL patients with chemotherapy or BCR inhibitors was shown to affect the expression of miRNAs involved in these processes [38,39,40]. However, the pathways regulated by most miRNAs remain to be elucidated.
MiRNA expression profiling studies have disclosed correlations between certain miRNA signatures and cytogenetic features and/or IGHV gene mutational status [30,41], which are recognized prognostic markers in CLL. Finally, certain miRNA signatures are associated with disease progression and outcome [24,30,34,42,43] or with the onset of Richter’s transformation [44,45,46], a lethal condition characterized by the development of aggressive lymphoma in CLL patients [1,47].
Specifically, miR-15a and miR-16-1, which are located on 13q14 [31], act as tumor suppressors and were the first miRNAs used to predict outcomes in CLL [24]. Consistent with the notion that miR-15 and miR-16 regulate cell apoptosis and proliferation [33,48], our group showed that transfection of miR-15 and miR-16 mimics into del(13)(q14) CLL cells significantly hindered their growth in NOD/Shi-scid, γ(c)(null) (NSG) mice, resulting in substantial tumor regression [49]. Moreover, our previous research revealed that CLL has a miRNA expression profile that closely resembles that of antigen-experienced B cells; some miRNAs in this profile are likely to influence disease progression, as suggested by their correlation with the clinical outcomes [30,50].
In the context of our prospective O-CLL1 study, we herein explored the predictive role of 513 miRNAs on TTFT by (i) incorporating significant miRNAs into a basic prognostic model including IGHV mutational status, del11q and del17p, beta-2-microglobulin (B2M), NOTCH1 mutation, and Rai stage, which are known to predict TTFT in CLL [51,52,53,54], and (ii) analyzing their potential association with 50 genes involved in TTFT that we previously identified by an artificial intelligence (AI)-based model. Our approach identified 16 out of 513 miRNAs independently associated with TTFT. Moreover, enrichment analysis was performed to evaluate biological processes regulated by the selected miRNAs, while correlation analysis and multiple bioinformatics software were used to evaluate possible interactions between the 16 prioritized miRNAs and the 50 genes selected as a CLL TTFT predictive signature, defined by an AI model [54].

2. Results

2.1. TTFT Prediction by the Basic Model

The estimated median TTFT for the 224 CLL cases from our prospective O-CLL1 study was 105.2 months (95% CI, 94.2–116.1 months), with 4-year and 8-year cumulative first-treatment-free survival of 76.8% (95% CI, 71.1–82.5) and 57.4% (95% CI, 50.1–64.6), respectively (Figure 1). Only six out of nine univariable predictors (i.e., Rai stage, B2M, IGVH and NOTCH1 mutational status, 11q and 17p deletions) remained significantly associated with TTFT when tested in a Cox multivariable model adopting a backward elimination strategy. Univariable and multivariable Cox analyses, focusing on established TTFT predictors, are detailed in Table S1A,B. The six variables together gave a prognostic value for TTFT (i.e., a Harrell’s C-index) of 75% and an explained variation of 45.4%.

2.2. TTFT Prediction by miRNAs

Of the 513 available miRNAs tested, 73 were found to be significantly associated with TTFT in univariable Cox regression analyses (all p ≤ 0.05) (Table 1). These 73 miRNAs were then simultaneously included in the same multivariable Cox regression model, and those significantly associated with TTFT were identified by a backward elimination strategy.
Sixteen miRNAs retained an independent association with TTFT (Table 2). For eight miRNAs (i.e., miR-582-3p, miR-33a-3p, miR-516a-5p, miR-99a-5p, and miR-296-3p, miR-502-5p, miR-625-5p, and miR-29c-3p), lower expression levels were associated with a higher likelihood of treatment (i.e., shorter TTFT). Conversely, for the remaining eight miRNAs (i.e., miR-150-5p, miR-148a-3p, miR-28-5p, miR-144-5p, miR-671-5p, miR-1-3p, miR-193a-3p, and miR-124-3p), a higher expression was associated with shorter TTFT.
Using these 16 miRNAs, a risk prediction score (ranging from 0 to 100%) was derived (Figure 2).
Inclusion of the miRNA scores into the basic model (i.e., a Cox model including Rai stage, B2M, IGVH and NOTCH1 mutational status, 11q and 17p deletions) significantly increased both the Harrell’s C-index (from 75.0% to 81.1%) and the explained variation in TTFT (from 45.4% to 63.3%). Remarkably, the inclusion of the miRNA scores into the basic model also yielded an IDI and an NRI of +14.9% and +44.2%, respectively (Table 3).
Supplementary Figure S1 illustrates the survival curve derived from the multivariable Cox proportional hazards model, which closely aligns with the Kaplan–Meier curve presented in Figure 1, thereby confirming the robustness and calibration of the model for predicting TTFT.

2.3. Correlation and Interaction Analysis between miRNAs and Genes Found to Be Related to TTFT by an AI Model

Next, we analyzed the relationship between the 16 miRNAs identified as predictive of TTFT and the previous 50 AI model prioritized genes (AI genes) that we found significantly associated with TTFT in the same patients of the O-CLL1 cohort [54]. Fifteen out of sixteen miRNAs showed at least one significant negative correlation with the AI genes; in particular, eight miRNAs (i.e., miR-29c-3p, miR-625-5p, miR-150-5p, miR-144-5p, miR-28-5p, and miR-516a-5p) showed a relationship with one of the top ten genes selected by the neural network (i.e., CEACAM19, PIGP, FADD, FIBP, IGF1R, COL28A1, QTRT1, MKL1, GNE, SLC39A6) (Figure 3). Only hsa-miR-124-3p showed no correlation with any of the AI genes; in addition, miR-29c-3p had a significant negative correlation with the highest number of genes prioritized by the AI model.
Comparing the previous correlation analysis with the miRComb model, this analysis retained only five miRNA–mRNA annotated interactions, displaying a mild to moderate correlation index. miR-29c-3p showed again the highest number of interactions, three of which were previously predicted; specifically, miR-29c-3p’s interaction with PRICKLE1 was annotated in the microCosm database, while interactions with ANKRD52 and ZBTB34 were annotated in the TargetScan database. miR-29c-3p’s interaction with KDM5B was also predicted by TargetScan with experimental validation, although with less strong evidence. MiR-625-5p was found to be significantly anti-correlated with IGF1R, an interaction annotated by miRTarBase with a strong validation score (Table S2). Interestingly, the MiRComb regulation score suggested an association between IGHV mutational status and the interactions between miR-29c-3p with the previously identified AI model-selected genes KDM5B, ANKRD52, ZBTB34, PRICKLE1; this dysregulation was not evident when considering the interaction between miR-625-5p and IGF1R (Figure 4). Although the magnitude of the effect of IGHV mutational status on the former interaction is strong, this is largely due to an over-expression of miR-29c-3p in the IGHV mutated group compared to the unmutated one and, to a lesser extent, to a slightly not significant downregulation of the previous four AI-selected genes in the IGHV mutated group compared to the IGHV unmutated one.
The aforementioned miRNA–mRNA interactions represented by the miRComb model were found to be specific for CLL. Although the statistical power might be affected by the small sample size, miRNA and mRNA expression profiles of tonsil B cells from normal individuals did not reveal any significant interactions between any of the 16 miRNAs and 50 AI-selected genes.

2.4. Pathways Regulation by the 16 miRNAs Linked to TTFT

To further investigate the biological significance of the 16 identified miRNAs, an enrichment analysis was carried out, showing that selected miRNAs are involved in cancer and cellular processes associated with disease progression and drug resistance.
Specifically, Kegg term overrepresentation analysis revealed a significant pathway (i.e., hsa05206) corresponding to the “MicroRNAs in cancer” pathway, associated with the miRNA genes MIR1-2, MIR28, MIR29C, MIR99A, MIR124-3, MIR150, and MIR625 encoding for miR-1-3p, miR-28-5p, miR-29c-3p, miR-99a-5p, miR-124-3p, and miR-150-5p, respectively (Table 4).
Overrepresentation analysis using Wikipath revealed six significant terms associated with the miRNA genes MIR29C, MIR33A, and MIR150 (Table 4). The terms WP1601, WP2249, and WP1545, signatures of “Fluoropyrimidine activity”, “Metastatic brain tumor”, and “miRNAs involved in DNA damage response”, respectively, were enriched in the MIR29C gene product; the terms WP299 and WP430, signatures of “Nuclear receptors in lipid metabolism and toxicity” and “Statin inhibition of cholesterol production”, respectively, were enriched in the MIR33A gene product. Finally, the WP2023 term, the signature of the “Cell differentiation expanded index”, was enriched in the MIR150 gene product.

3. Discussion

Predicting overall and progression-free survival in CLL is a dynamic field of research that adapts to advances in biology and therapeutics, providing insights to personalized patient planning and guiding clinicians in tailored treatment decisions. Conversely, models predicting TTFT, which is not influenced by the treatment choice, play a pivotal role in counseling, family planning, surveillance, and identification of high-risk candidates for potential early intervention in clinical trials [55,56]. Indeed, several predictive models have been published for this purpose [57,58,59]. Recently, the European Research Initiative on CLL (ERIC) introduced a score that categorizes patients based on IGHV mutation and somatic hypermutation status [18].
Nevertheless, most of these predictive models typically require a combination of clinically accessible factors and molecular factors, which are relatively straightforward to assess. In the O-CLL1 trial, multivariable analyses identified a foundational model comprising Rai stage, B2M, IGHV, and NOTCH1 mutational status, in addition to 17p and 11q deletions [51,52,53,54]. This model also served as the basis for validating the improved prognostic efficacy of lymphocyte doubling time (LDT) [53] and, more recently, for incorporating our previously discovered AI gene expression data into an enhanced model predicting TTFT [54].
The current study focused on analyzing the profile of miRNAs at the time of diagnosis among treatment-naive CLL cases enrolled in the prospective O-CLL1 trial. The primary aim was to identify specific miRNAs that could function as predictive markers to better identify patients in need of early treatment. The establishment of a risk prediction score based on 16 miRNAs identified as independently associated with TTFT emphasizes the significant impact of miRNA data in predicting TTFT in CLL. Specifically, lower expression of eight miRNAs (i.e., miR-582-3p, miR-33a-3p, miR-516a-5p, miR-99a-5p, miR-296-3p, miR-502-5p, miR-625-5p, and miR-29c-3p) was associated with a shorter TTFT, while for the remaining eight miRNAs (i.e., miR-150-5p, miR-148a-3p, miR-28-5p, miR-144-5p, miR-671-5p, miR-1-3p, miR-193a-3p, and miR-124-3p), the higher expression was associated with a shorter TTFT.
Importantly, the integration of these miRNAs into the basic model significantly improved predictive accuracy, as reflected by the improved Harrell’s C-index from 75.0% to 81.1% and explained variation in TTFT from 45.4% to 63.3%. Additionally, this integration yielded an IDI of +14.9% and an NRI of +44.2%. This suggests that miRNA expression data may provide meaningful insights, offering a promising opportunity for counseling, and a more reasonable follow-up scheduling through early identification of cases with a higher likelihood of requiring early therapy.
To date, miR-15a and miR-16-1, localized on chromosome 13q14 and functioning as bona fide tumor suppressors [31,33,48], have been utilized as pioneering prognostic indicators in CLL [24]. Subsequent investigations have highlighted the potential of miRNA profiling in enhancing the accuracy of CLL prognostication. Specifically, the expression levels of dysregulated miR-155, miR-181b, miR-29a/b/c, and miR-34a have been systematically correlated with established prognostic biomarkers, including IGHV and TP53 mutational status, as well as ZAP70 expression, thereby exerting a discernible impact on the clinical outcomes of CLL patients [25,28,41,43,60,61]. Indeed, our previous study unveiled several dysregulated miRNAs in CLL, indicating their potential role in the pathogenesis of the disease and their contribution to its progression, ultimately influencing the initiation of therapy [30]. In particular, 8 miRNAs (i.e., miR-146b-5p, miR-222-3p, miR-26a-5p, miR-29c-3p, miR-29c-5p, miR-502-3p, and miR-503-5p) out of 15 dysregulated miRNAs demonstrated a significant role in predicting TTFT, with miR-26a-5p and miR-532-3p remaining significantly associated with TTFT after data adjustment for confounders.
Notably, in our current study, only miR-29c-3p was found to be significantly associated with TTFT after applying the Cox multivariable regression model. However, it is essential to consider that, in contrast to our previous analyses and other published studies, our current investigation diverges from the conventional focus based on dysregulated miRNAs. Instead, we opted for a broader statistical strategy, encompassing the evaluation of all the assessable miRNAs to gain a more comprehensive understanding of the miRNA landscape in the context of TTFT prediction. In the present investigation, we have delineated a distinctive signature comprising 16 miRNAs that influence the risk of initiating therapy in CLL. Of note, elevated levels of eight of them correlated with higher treatment likelihood, suggesting their role in disease progression, while the remaining showed an inverse relationship, lowering treatment initiation likelihood, indicating their dual ability to act either as promoters of oncogenic processes or as regulators of tumor suppression.
We reasoned that the understanding of the regulatory networks involving the identified 16 miRNAs could be of interest to further elucidate the CLL pathogenesis; therefore, we adopted a comprehensive in silico approach employing correlation, interaction, and enrichment models to validate and refine miRNA–mRNA analyses in our O-CLL1 dataset. The miRNA–mRNA correlation analysis revealed at least one significant negative correlation between 15 of the identified miRNAs and the set of 50 AI model-based genes we previously identified, with miR-29c-3p being the one with the highest number of correlations as confirmed by multiple tools. miR-625-5p was found to be significantly anti-correlated only with IGF1R; interestingly, IGF-1R was functionally validated as a direct target of miR-625-5p also in melanoma cancer cells, being positively modulated by the long non-coding RNA LINC01291 via miR-625-5p sponging [62]. Of particular interest, miR-29c was previously found to exhibit differential expression based on the IGHV status [24,30], which correlated with the expression of some of its specific target genes [61,63]. The present study expanded this observation to four AI-generated genes (ANKRD52, KDM5B, PRICKLE1, and ZBTB34) previously described to be strongly associated with TTFT in CLL [54]. Regulation of KDM5B by miR-29c was described in endometrial carcinoma, showing a correlation between elevated levels of KDM5B and tumor grade and paclitaxel resistance [64]. Finally, the enrichment analysis revealed the multifaceted involvement of miRNAs in cellular activities, from tumor development and resistance to lipid metabolism regulation.
Our previous investigation demonstrated that IL-23R expression on CLL cells independently estimated TTFT in the O-CLL1 cohort, supporting the notion of an essential IL-23 autocrine loop driving CLL expansion [65]. Furthermore, we found that miR-146b-5p modulates IL-12Rβ1 expression, which influences TTFT, with lower levels associated with shorter TTFT duration [50]. Notably, in the current study, miR-148a-3p emerged as an independent prognostic factor in the multivariable model, showcasing a striking complementarity of 18 out of 21 bases with a sequence within the interleukin 23 subunit alpha (IL-23A) mRNA (miRWalk: refseq ID NM_016584, http://mirwalk.umm.uni-heidelberg.de/human/gene/51561/ accessed on 21 August 2024), thus expanding the exploration of IL23A gene regulation.
In conclusion, the identification of specific miRNAs as predictors of TTFT in CLL represents a promising avenue for refining risk stratification and predicting therapeutic needs. The integration of miRNA data into predictive models holds the potential to improve the accuracy of clinical decision-making in CLL management. However, further validation studies and comprehensive functional analyses are needed to confirm the reliability and robustness of these findings.

4. Materials and Methods

4.1. Patient Population and Study Design

In the observational O-CLL1 study (clinicaltrials.gov identifier NCT00917540), 224 newly diagnosed Binet A CLL cases from 40 Italian institutions were prospectively enrolled for miRNAs analysis [30]. The miRNA expression data are deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus repository (http://www.ncbi.nlm.nih.gov/geo/ accessed on 21 August 2024) and are accessible via GEO Series accession number GSE40533. The gene expression data are accessible through GEO Series accession number GSE40570.
All participants gave written informed consent, and the study was approved by the appropriate institutional review boards. Detailed inclusion and exclusion criteria have been previously outlined [54]. Specifically, the recruitment was limited to cases diagnosed within 12 months, age ≤ 70 years, and at Binet stage A [30].

4.2. Assessment of Biological Markers

The diagnosis was confirmed by flow cytometric analysis, which assessed the proportion of CD5/CD19/CD23 triple-positive B cells; monoclonal antibodies (mAbs) against CD19-FITC (BD Biosciences Pharmigen, San Jose, CA, USA), CD23-PE (BD Biosciences, San Jose, CA), and CD5-PC5 (Beckman Coulter Immunotech, Marseille, France) were used for this purpose. CD38-positive leukemic cells were quantified through triple staining with CD19 FITC (BD Biosciences), CD38 PE (BD Biosciences), and CD5 PC5 (Beckman Coulter Immunotech) mABs. CD38 positivity was defined with a cutoff of ≥20%, following previous reports [30].
ZAP-70 was detected by flow cytometry using a ZAP-70-FITC (clone 2F3.2, Millipore, Temecula, CA, USA) or an isotype control mAb (mouse IgG2a-FITC, BD Biosciences) as previously described [30]. Briefly, peripheral mononuclear cells purified from fresh heparinized CLL samples by Ficoll-Hypaque gradient were first incubated with CD3 PE-CY7, CD19 PE, and CD5 PC5 mAbs (BD Biosciences, San Jose, CA, USA), fixed, permeabilized with fix and perm reagents (Caltag Laboratories, Buckingham, UK), and exposed to ZAP-70 or the isotype control mAb. A cutoff value of more than 30% for ZAP-70 positivity was used, as previously reported [30] and calculated by receiver-operating characteristic analysis as the most suitable ZAP-70 cutoff value to discriminate IGHV-Unmutated (UM) from IGHV-Mutated (M) cases. All flow cytometric analyses were performed on a FACSCalibur flow cytometer (BD Biosciences).
Cytogenetic abnormalities, including chromosome deletions 11q23 and 17p13, were examined by fluorescence in situ hybridization (FISH) in a purified CD19-positive population, according to established protocols [30]. The mutational status of the IGHV gene was assessed on cDNA samples [30] by aligning sequences to the IMGT directory and by analyzing them using IMGT/V-QUEST software(version: 3.6.3 (30 January 2024). The NOTCH1 mutation hotspot was determined by next-generation deep sequencing, as previously outlined [51].

4.3. miRNAs Analysis

Total RNA was extracted from CD19-positive purified B-cell samples using TRizol reagent (Life Technologies, Carlsbad, CA, USA). Subsequently, RNA quality was evaluated through the Agilent 2100 Bioanalyzer (Agilent Technologies). Any RNA sample with poor quality, as indicated by an RNA integrity number < 7, was excluded from microarray analyses. After sample collection, processing was performed according to the guidelines of the Agilent manual for the Human miRNA Microarray V2 platform (Agilent Technologies, Palo Alto, CA, USA), which includes miRNAs from the Sanger miRBase (v10.1). Expression values of miRNAs were computed using Agilent Feature Extraction Software 10.1, followed by summarization and background subtraction. MiRNAs with low expression (all detection calls missing), non-human miRNAs, and miRNAs expired according to Sanger miRBase release 15 (April 2010) were excluded. The obtained raw data underwent quantile normalization, conversion to positive values with a minimum value of 1, and log2 transformation using the R-2.14 statistical environment (http://www.r-project.org/ accessed on 21 August 2024) [30]. Array old identifiers were converted to the last version with Mirbase (i.e., miR-33a* = miR-33a-3p, miR-miR-99a = 99a-5p, miR-625 = miR-625-5p, miR-29c = miR-29c-3p, miR-150 = miR-150-5p, miR-148a = miR-148a-3p, miR-144* = miR-144-5p, miR-1 = miR-1-3p, and miR-124 = miR-124-3p).

4.4. miRNA–mRNA Correlation, Interaction, and Enrichment Analyses

To elucidate potential relationships between the miRNAs significantly associated with TTFT and the previously identified 50 genes implicated in TTFT as detected by an AI-based model [54], a Spearman correlation analysis between miRNA and mRNA expression profiles was performed using the O-CLL1 database [30]. To retain miRNA–mRNA correlations implicating an miRNA–mRNA interaction, an additional Spearman correlation analysis between expression profiles was performed using the miRComb model [66], selecting significant correlations that were suggestive of an interaction annotated in at least one of the following prediction databases, namely microCosm (version 5.18) and TargetScan 6.2. The sensitivity of the previous model was improved by extending the annotation query to the miRTarBase database [67], which annotated experimentally proven miRNA–mRNA interactions accompanied by a qualitative validation score including the two values “Low evidence” and “Strong evidence” based on the experimental techniques involved in the validation protocol. At the same time, the specificity of the model was increased by selecting correlations with a corresponding p-value, adjusted for multiple comparisons inherent to the whole gene expression and miRNA profile analysis, below a threshold of 0.05. Additionally, the mirComb model was applied to the expression dataset to calculate an miRNA–mRNA regulation score, which allowed the assessment of the effect of IGHV mutational status on each miRNA–mRNA interaction retained; the regulation score was calculated as follows:
s c o r e = − 2 log 2 ⁡ r a t i o m R N A   ·   log 2 ⁡ r a t i o m i R N A
where log2 ratiomRNA and log2 ratiomiRNA equal to the log-transformed fold changes of mRNA and miRNA expressions of a given interaction pair, respectively, calculated by comparing the IGHV mutated group with the unmutated group. Specifically, the previous score measured the magnitude of interaction deregulation between an miRNA and an mRNA as a result of the effect mediated by a categorical variable: a positive score implied interaction deregulation, while a negative score indicated a preserved interaction. To ensure the specificity of the miRNA–mRNA interactions detected from the O-CLL1 cohort for CLL pathogenesis, a prior correlation analysis was also carried out on a control set of normal B-cell sample subpopulations collected from tonsils of healthy individuals. In particular, miRNA and mRNA expression profiles were collected for 3 samples of Naive B cells, 2 samples of Marginal Zone-like B cells, and 3 samples of memory cells (GSE51529).
To elucidate the biological significance of the identified miRNAs that are correlated with clinical outcomes, we conducted an over-representation analysis (ORA) of Kyoto Encyclopedia of Genes and Genomes (KEGG) and Wikipathways terms. This analysis was performed using the R package ClusterProfiler, a versatile tool designed for the interpretation of omics data [68]. ORA allows the identification of enriched biological pathways or processes associated with a given gene set, thereby providing insight into the potential functional roles of the miRNAs under investigation.

4.5. Statistical Analyses

TTFT was calculated during the watch-and-wait period, which lasted from the date of diagnosis to the start of therapy or the last follow-up. The prognostic impact of standard risk factors and miRNAs for TTFT was preliminarily investigated by univariable Cox regression analyses, with data presented as hazard ratios (HRs) with 95% confidence intervals (CIs). To obtain parsimonious models, all univariate predictors of TTFT (i.e., variables with p ≤ 0.05) were tested in Cox analyses with a backward elimination strategy to identify relevant prognostic variables among standard risk factors or miRNAs. Harrell’s C-index, the explained variation in TTFT (an index combining calibration and discrimination), the integrated discrimination improvement (IDI), and the net reclassification index (NRI) were used to assess the accuracy of prognostic models and to measure the gain in prognostic accuracy attributable to relevant miRNAs. To generate a risk prediction rule based on miRNAs, a logistic regression model was fitted with TTFT as the dependent variable and the miRNAs which remained significantly associated with the study outcome by the backward elimination strategy as independent variables. This analysis resulted in an miRNA score (ranging from 0 to 100%), which was then used for further analyses. Statistical calculations were performed using SPSS for Windows v.21 (IBM, Chicago, IL, USA) and Stata 16 (StataCorp, College Station, TX, USA).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ncrna10050046/s1, Table S1: Univariable (A) and Multivariable Cox Regression (B) analyses for TTFT with identification of significant Independent Predictors via backward elimination strategy. Table S2: Evidence of interactions between genes and miRNAs shown in Figure 3 and Figure 4 according to mirTarbase, MicroCosm, and TargetScan. Figure S1: survival curve derived from the multivariable Cox proportional hazards model incorporating the risk score.

Author Contributions

Conceptualization, E.N., F.R., A.N., F.M. and G.C.; methodology, E.N., F.R., A.A.A., F.M., G.T., G.D, G.G., C.A. and P.M.; software E.N., F.R., A.A.A., G.G., G.T., C.A., G.D. and F.M.; formal analysis, E.N., F.R., F.M., A.A.A., G.T., G.D., G.G. and C.A.; investigation, EN., F.R., M.N., A.N., M.C., N.B., F.F. (Fabiana Ferrero), F.T., N.P., M.P., G.C. and F.M.; resources, A.B., E.A.M., E.V., M.G., F.F., M.N., M.F., G.T. and C.A.; data curation, M.N., A.N., F.M. and E.N.; writing—original draft preparation, F.M., G.C., E.N. and F.R.; writing—review and editing, A.N., M.F., P.M., A.A.A., F.F. (Franco Fais), F.M. and G.C.; visualization, E.N., F.R., F.M. and G.C.; supervision, F.M, G.C., A.A.A., M.F. and A.N.; project administration, F.M., G.C., A.N. and M.F.; funding acquisition, F.M., M.F., A.N., M.N., F.F. (Franco Fais), N.B. and G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Italian Association for Cancer Research (AIRC) Grant 5 × mille ID.9980, (M.F., F.M., M.N. and A.N.); AIRC I.G. ID.1426 (to M.F.), ID.24365 (to A.N.), ID.15426 (to F.F.); AIRC and Fondazione CaRiCal co-financed Multi-Unit Regional Grant 2014 n.16695 (to F.M.); Associazione Italiana contro le Leucemie-Linfomi e Mieloma (AIL, Cosenza, to F.F.); Gilead fellowship program 2016 (M.C.) and 2017 (G.C.); Italian Ministry of Health 5 × 1000 funds 2014 (to G.C.), 2015 (to F.F.), 2016 (to F.F. and G.C.), 2020 (to G.C. and P.M.), 2021 (to N.B.); Italian Ministry of Health RF-2021-12374376 (to G.C. and A.N.) and by Italian Ministry of Health (Ricerca Corrente 2022–2024).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committees of each Institution participating in the O-CLL1 Multicenter observational trial “Prospective Collection of Biological Data of Prognostic Relevance in Patients with B-Cell Chronic Lymphocytic Leukemia” (clinicaltrials.gov identifier NCT00917540) as mentioned in the Acknowledgments section. Specifically, the Ethics committee of the National Institute for Cancer Research-IST of Genoa-Italy (now called IRCCS Ospedale Policlinico San Martino) approved the study on 04/16/2007 with no. OMC007.002 for the centralized laboratory that processed and stored all the CLL samples for the subsequent biological studies. Since July 2020, the study for the Regional Ethics Committee has been identified as N. CER Liguria id 10115.

Informed Consent Statement

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

Data Availability Statement

Datasets used in this study are publicly available and can be accessed from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus repository (http://www.ncbi.nlm.nih.gov/geo/) accessed on 21 August 2024 and are accessible via GEO Series accession number GSE40533, GSE51529 and GSE40570.

Acknowledgments

In addition to the listed authors, the following investigators participated in this study as part of the Gruppo Italiano Studio Linfomi (GISL): Gianni Quintana, Divisione di Ematologia, Presidio Ospedaliero “A.Perrino”, Brindisi; Giovanni Bertoldero, Dipartimento di Oncologia, Ospedale Civile, Noale, Venezia; Paolo Di Tonno, Dipartimento di Ematologia, Ospedale di Venere, Bari; Robin Foà and Francesca R Mauro, Divisione di Ematologia, Università La Sapienza, Roma; Nicola Di Renzo, Unità di Ematologia, Ospedale Vito Fazzi, Lecce; Maria Cristina Cox, Ematologia, A.O. Sant’Andrea, Università La Sapienza, Roma; Stefano Molica, Dipartimento di Oncologia ed Ematologia, Pugliese-Ciaccio Hospital, Catanzaro; Attilio Guarini, Unità di Ematologia e Trapianto di Cellule Staminali, Istituto di Oncologia “Giovanni Paolo II”, Bari; Antonio Abbadessa, U.O.C. di Oncoematologia Ospedale “S. Anna e S. Sebastiano”, Caserta; Francesco Iuliano, U.O.C. di Oncologia, Ospedale Giannettasio, Rossano Calabro, Cosenza; Omar Racchi, Ospedale Villa Scassi Sampierdarena, Genova; Mauro Spriano, Ematologia, A.O. San Martino, Genova; Felicetto Ferrara, Divisione di Ematologia, Ospedale Cardarelli, Napoli; Monica Crugnola, Ematologia, CTMO, Azienda Ospedaliera Universitaria di Parma; Alessandro Andriani, Dipartimento di Ematologia, Ospedale Nuovo Regina Margherita, Roma; Nicola Cascavilla, Unità di Ematologia e Trapianto di Cellule Staminali, IRCCS Ospedale Casa Sollievo della Sofferenza, San Giovanni Rotondo; Lucia Ciuffreda, Unità di Ematologia, Ospedale San Nicola Pellegrino, Trani; Graziella Pinotti, U.O. Oncologia Medica, Ospedale di Circolo Fondazione Macchi, Varese; Anna Pascarella, Unità Operativa di Ematologia, Ospedale dell’Angelo, Venezia-Mestre; Maria Grazia Lipari, Divisione di Ematologia, Ospedale Policlinico, Palermo, Francesco Merli, Unità Operativa di Ematologia, A.O.S. Maria Nuova, Reggio Emilia; Luca Baldini Istituto di Ricovero e Cura a Carattere Scientifico Cà Granda-Maggiore Policlinico, Milano; Caterina Musolino, Divisione di Ematologia, Università di Messina; Agostino Cortelezzi, Ematologia and CTMO, Foundation IRCCS Ca’ Granda Ospedale Maggiore Policlinico, Milano; Francesco Angrilli, Dipartimento di Ematologia, Ospedale Santo Spirito, Pescara; Ugo Consoli, U.O.S. di Emato-Oncologia, Ospedale Garibaldi-Nesima, Catania; Gianluca Festini, Centro di Riferimento Ematologico-Seconda Medicina, Azienda Ospedaliero-Universitaria, Ospedali Riuniti, Trieste; Giuseppe Longo, Unità di Ematologia, Ospedale San Vincenzo, Taormina; Daniele Vallisa and Annalisa Arcari, Unità di Ematologia, Dipartimento di Onco-Ematologia, Guglielmo da Saliceto Hospital, Piacenza; Francesco Di Raimondo and Annalisa Chiarenza, Divisione di Ematologia, Università di Catania Ospedale Ferrarotto, Catania; Iolanda Vincelli, Unità di Ematologia, A.O. of Reggio Calabria; Donato Mannina, Divisione di Ematologia, Ospedale Papardo, Messina, Italy. For Each Institution, the Ethics Review Committee that approved the protocol is listed below. Unità di Ematologia e Trapianto di Cellule Staminali, Istituto di Oncologia “Giovanni Paolo II”, Bari comitato etico Istituto di Oncologia “Giovanni Paolo II”, di Bari, Italy; Dipartimento di Ematologia, Ospedale di Venere, Bari comitato etico Ospedale di Venere, Bari, Italy; U.O.C. Ematologia e Trapianto, Ospedale “Mons. R. Dimiccoli” -Barletta comitato etico Ospedale “Mons. R. Dimiccoli” -Barletta, Italy; Divisione di Ematologia, Presidio Ospedaliero “A. Perrino”, Brindisi comitato etico Presidio Ospedaliero “A. Perrino”, Brindisi, Italy; U.O.C. di Oncoematologia Ospedale “S. Anna e S. Sebastiano”, Caserta comitato etico Ospedale “S. Anna e S. Sebastiano”, Caserta, Italy; Divisione di Ematologia, Università di Catania Ospedale Ferrarotto, Catania comitato etico Ospedale Ferrarotto, Catania, Italy; U.O.S. di Emato-Oncologia, Ospedale Garibaldi-Nesima, Catania comitato etico Ospedale Garibaldi-Nesima, Catania, Italy, Dipartimento di Oncologia ed Ematologia, Pugliese-Ciaccio Hospital, Catanzaro comitato etico Pugliese-Ciaccio Hospital, Catanzaro, Italy; U.O.C di Ematologia, A.O. Cosenza comitato etico A.O. Cosenza, Italy; U.O.C. di Oncologia, Ospedale Giannettasio, Rossano Calabro, Cosenza comitato etico Ospedale Giannettasio, Rossano Calabro, Italy; Divisione di Ematologia, Ospedale Policlinico, Palermo, comitato etico Ospedale Policlinico, Palermo, Italy; Ematologia ospedale Goretti, Latina comitato etico ospedale Goretti, Latina, Italy; Clinica Ematologica, DIMI, Genova, comitato etico Ospedale San Martino, Genova, Italy; Oncologia medica C IRCCS Ospedale Policlinico San Martino, Genoa, Italy comitato etico IRCCS Ospedale Policlinico San Martino, Genova, Italy; Ospedale Villa Scassi Sampierdarena, Genova comitato etico Ospedale Villa Scassi, Genova, Italy; Ematologia, A.O. San Martino, Genova comitato etico A.O. San Martino, Genova, Italy; Unità di Ematologia, Ospedale Vito Fazzi, Lecce comitato etico Ospedale Vito Fazzi, Lecce, Italy; UOC di Ematologia Ospedale di Matera comitato etico Ospedale di Matera, Italy; Divisione di Ematologia, Ospedale Papardo, Messina, Italy. comitato etico Ospedale Papardo, Messina, Italy; Divisione di Ematologia, Università di Messina comitato etico Ospedale di Messina, Italy; Ematologia and CTMO, Foundation IRCCS Ca’ Granda Ospedale Maggiore Policlinico, Milano comitato etico Ospedale Maggiore Policlinico di Milano, Italy; Dipartimento di Oncologia, Ospedale Civile, Noale, Venezia comitato etico Ospedale Civile, Noale, Italy; Oncoematologia Policlinico di Modena comitato etico provinciale di Modena, Italy; Divisione di Ematologia, Ospedale Cardarelli, Napoli comitato etico Ospedale Cardarelli, Napoli, Italy; Ematologia, CTMO, Azienda Ospedaliera Universitaria di Parma comitato etico Azienda Ospedaliera Universitaria di Parma, Italy; Dipartimento di Ematologia, Ospedale Santo Spirito, Pescara comitato etico Ospedale Santo Spirito, Pescara, Italy; Unità di Ematologia, Dipartimento di Onco-Ematologia, Guglielmo da Saliceto Hospital, Piacenza comitato etico Guglielmo da Saliceto Hospital, Piacenza, Italy; Unità di Ematologia, A.O. of Reggio Calabria comitato etico A.O. of Reggio Calabria, Italy; Unità Operativa di Ematologia, A.O.S. Maria Nuova, Reggio Emilia comitato etico A.O.S. Maria Nuova, Reggio Emilia, Italy; S.C. di Ematologia e Trapianto di Cellule Staminali IRCCS-CROB di Rionero in Vulture, Potenza comitato etico IRCCS-CROB di Rionero in Vulture, Italy; Dipartimento di Ematologia, Ospedale Nuovo Regina Margherita, Roma comitato etico Ospedale Nuovo Regina Margherita, Roma, Italy; Ematologia, A.O. Sant’Andrea, Università La Sapienza, Roma comitato etico A.O. Sant’Andrea, Roma, Italy; Divisione di Ematologia, Università La Sapienza, Roma comitato etico Università La Sapienza, Roma, Italy; Unità di Ematologia e Trapianto di Cellule Staminali, IRCCS Ospedale Casa Sollievo della Sofferenza, San Giovanni Rotondo comitato etico IRCCS Ospedale Casa Sollievo della Sofferenza, San Giovanni Rotondo, Italy; Unità di Ematologia, Ospedale San Vincenzo, Taormina, comitato etico Ospedale San Vincenzo, Taormina, Italy; Unità di Ematologia, Ospedale San Nicola Pellegrino, Trani comitato etico Ospedale San Nicola Pellegrino, Trani, Italy; Centro di Riferimento Ematologico-Seconda Medicina, Azienda Ospedaliero-Universitaria, Ospedali Riuniti, Trieste comitato etico Ospedali Riuniti, Trieste, Italy; U.O. Oncologia Medica, Ospedale di Circolo Fondazione Macchi, Varese comitato etico Ospedale di Circolo Fondazione Macchi, Varese, Italy; Unità Operativa di Ematologia, Ospedale dell’Angelo, Venezia-Mestre comitato etico Ospedale dell’Angelo, Venezia-Mestre, Italy; Istituto di Ricovero e Cura a Carattere Scientifico Cà Granda-Maggiore Policlinico, Milano comitato etico Policlinico, Milano, Italy.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Hallek, M.; Cheson, B.D.; Catovsky, D.; Caligaris-Cappio, F.; Dighiero, G.; Dohner, H.; Hillmen, P.; Keating, M.; Montserrat, E.; Chiorazzi, N.; et al. iwCLL guidelines for diagnosis, indications for treatment, response assessment, and supportive management of CLL. Blood 2018, 131, 2745–2760. [Google Scholar] [CrossRef]
  2. Hallek, M.; Shanafelt, T.D.; Eichhorst, B. Chronic lymphocytic leukaemia. Lancet 2018, 391, 1524–1537. [Google Scholar] [CrossRef] [PubMed]
  3. Kipps, T.J.; Stevenson, F.K.; Wu, C.J.; Croce, C.M.; Packham, G.; Wierda, W.G.; O’Brien, S.; Gribben, J.; Rai, K. Chronic lymphocytic leukaemia. Nat. Rev. Dis. Primers 2017, 3, 16096. [Google Scholar] [CrossRef] [PubMed]
  4. Gaidano, G.; Rossi, D. The mutational landscape of chronic lymphocytic leukemia and its impact on prognosis and treatment. Hematol. Am. Soc. Hematol. Educ. Program. 2017, 2017, 329–337. [Google Scholar] [CrossRef]
  5. Guieze, R.; Wu, C.J. Genomic and epigenomic heterogeneity in chronic lymphocytic leukemia. Blood 2015, 126, 445–453. [Google Scholar] [CrossRef]
  6. Dohner, H.; Stilgenbauer, S.; Benner, A.; Leupolt, E.; Krober, A.; Bullinger, L.; Dohner, K.; Bentz, M.; Lichter, P. Genomic aberrations and survival in chronic lymphocytic leukemia. N. Engl. J. Med. 2000, 343, 1910–1916. [Google Scholar] [CrossRef] [PubMed]
  7. Landau, D.A.; Tausch, E.; Taylor-Weiner, A.N.; Stewart, C.; Reiter, J.G.; Bahlo, J.; Kluth, S.; Bozic, I.; Lawrence, M.; Bottcher, S.; et al. Mutations driving CLL and their evolution in progression and relapse. Nature 2015, 526, 525–530. [Google Scholar] [CrossRef]
  8. Puente, X.S.; Bea, S.; Valdes-Mas, R.; Villamor, N.; Gutierrez-Abril, J.; Martin-Subero, J.I.; Munar, M.; Rubio-Perez, C.; Jares, P.; Aymerich, M.; et al. Non-coding recurrent mutations in chronic lymphocytic leukaemia. Nature 2015, 526, 519–524. [Google Scholar] [CrossRef]
  9. Puente, X.S.; Pinyol, M.; Quesada, V.; Conde, L.; Ordonez, G.R.; Villamor, N.; Escaramis, G.; Jares, P.; Bea, S.; Gonzalez-Diaz, M.; et al. Whole-genome sequencing identifies recurrent mutations in chronic lymphocytic leukaemia. Nature 2011, 475, 101–105. [Google Scholar] [CrossRef]
  10. Baliakas, P.; Mattsson, M.; Stamatopoulos, K.; Rosenquist, R. Prognostic indices in chronic lymphocytic leukaemia: Where do we stand how do we proceed? J. Intern. Med. 2016, 279, 347–357. [Google Scholar] [CrossRef]
  11. Eichhorst, B.; Robak, T.; Montserrat, E.; Ghia, P.; Niemann, C.U.; Kater, A.P.; Gregor, M.; Cymbalista, F.; Buske, C.; Hillmen, P.; et al. Chronic lymphocytic leukaemia: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann. Oncol. 2021, 32, 23–33. [Google Scholar] [CrossRef] [PubMed]
  12. Gonzalez-Gascon, Y.M.I.; Munoz-Novas, C.; Rodriguez-Vicente, A.E.; Quijada-Alamo, M.; Hernandez-Sanchez, M.; Perez-Carretero, C.; Ramos-Ascanio, V.; Hernandez-Rivas, J.A. From Biomarkers to Models in the Changing Landscape of Chronic Lymphocytic Leukemia: Evolve or Become Extinct. Cancers 2021, 13, 1782. [Google Scholar] [CrossRef] [PubMed]
  13. Gentile, M.; Shanafelt, T.D.; Mauro, F.R.; Reda, G.; Rossi, D.; Laurenti, L.; Del Principe, M.I.; Cutrona, G.; Angeletti, I.; Coscia, M.; et al. Predictive value of the CLL-IPI in CLL patients receiving chemo-immunotherapy as first-line treatment. Eur. J. Haematol. 2018, 101, 703–706. [Google Scholar] [CrossRef]
  14. International, C.L.L.I.P.I.w.g. An international prognostic index for patients with chronic lymphocytic leukaemia (CLL-IPI): A meta-analysis of individual patient data. Lancet Oncol. 2016, 17, 779–790. [Google Scholar] [CrossRef]
  15. Condoluci, A.; Terzi di Bergamo, L.; Langerbeins, P.; Hoechstetter, M.A.; Herling, C.D.; De Paoli, L.; Delgado, J.; Rabe, K.G.; Gentile, M.; Doubek, M.; et al. International prognostic score for asymptomatic early-stage chronic lymphocytic leukemia. Blood 2020, 135, 1859–1869. [Google Scholar] [CrossRef]
  16. Smolej, L.; Turcsanyi, P.; Kubova, Z.; Zuchnicka, J.; Mihalyova, J.; Simkovic, M.; Vodarek, P.; Krcmeryova, M.; Mocikova, H.; Brejcha, M.; et al. External validation of International Prognostic Score for asymptomatic early stage chronic lymphocytic leukaemia and proposal of an alternative score. Br. J. Haematol. 2021, 193, 133–137. [Google Scholar] [CrossRef]
  17. Morabito, F.; Tripepi, G.; Vigna, E.; Bossio, S.; D’Arrigo, G.; Martino, E.A.; Storino, F.; Recchia, A.G.; Fronza, G.; Di Raimondo, F.; et al. Validation of the Alternative International Prognostic Score-E (AIPS-E): Analysis of Binet stage A chronic lymphocytic leukemia patients enrolled into the O-CLL1-GISL protocol. Eur. J. Haematol. 2021, 106, 831–835. [Google Scholar] [CrossRef]
  18. Mansouri, L.; Thorvaldsdottir, B.; Sutton, L.A.; Karakatsoulis, G.; Meggendorfer, M.; Parker, H.; Nadeu, F.; Brieghel, C.; Laidou, S.; Moia, R.; et al. Different prognostic impact of recurrent gene mutations in chronic lymphocytic leukemia depending on IGHV gene somatic hypermutation status: A study by ERIC in HARMONY. Leukemia 2023, 37, 339–347. [Google Scholar] [CrossRef] [PubMed]
  19. Bayraktar, R.; Fontana, B.; Calin, G.A.; Nemeth, K. miRNA Biology in Chronic Lymphocytic Leukemia. Semin. Hematol. 2024, 61, 181–193. [Google Scholar] [CrossRef]
  20. Balatti, V.; Croce, C.M. MicroRNA dysregulation and multi-targeted therapy for cancer treatment. Adv. Biol. Regul. 2020, 75, 100669. [Google Scholar] [CrossRef]
  21. Autore, F.; Ramassone, A.; Stirparo, L.; Pagotto, S.; Fresa, A.; Innocenti, I.; Visone, R.; Laurenti, L. Role of microRNAs in Chronic Lymphocytic Leukemia. Int. J. Mol. Sci. 2023, 24, 12471. [Google Scholar] [CrossRef] [PubMed]
  22. Shang, R.; Lee, S.; Senavirathne, G.; Lai, E.C. microRNAs in action: Biogenesis, function and regulation. Nat. Rev. Genet. 2023, 24, 816–833. [Google Scholar] [CrossRef] [PubMed]
  23. Nemeth, K.; Bayraktar, R.; Ferracin, M.; Calin, G.A. Non-coding RNAs in disease: From mechanisms to therapeutics. Nat. Rev. Genet. 2024, 25, 211–232. [Google Scholar] [CrossRef] [PubMed]
  24. Calin, G.A.; Ferracin, M.; Cimmino, A.; Di Leva, G.; Shimizu, M.; Wojcik, S.E.; Iorio, M.V.; Visone, R.; Sever, N.I.; Fabbri, M.; et al. A MicroRNA signature associated with prognosis and progression in chronic lymphocytic leukemia. N. Engl. J. Med. 2005, 353, 1793–1801. [Google Scholar] [CrossRef] [PubMed]
  25. Stamatopoulos, B.; Meuleman, N.; Haibe-Kains, B.; Saussoy, P.; Van Den Neste, E.; Michaux, L.; Heimann, P.; Martiat, P.; Bron, D.; Lagneaux, L. microRNA-29c and microRNA-223 down-regulation has in vivo significance in chronic lymphocytic leukemia and improves disease risk stratification. Blood 2009, 113, 5237–5245. [Google Scholar] [CrossRef]
  26. Rossi, S.; Shimizu, M.; Barbarotto, E.; Nicoloso, M.S.; Dimitri, F.; Sampath, D.; Fabbri, M.; Lerner, S.; Barron, L.L.; Rassenti, L.Z.; et al. microRNA fingerprinting of CLL patients with chromosome 17p deletion identify a miR-21 score that stratifies early survival. Blood 2010, 116, 945–952. [Google Scholar] [CrossRef]
  27. Fulci, V.; Chiaretti, S.; Goldoni, M.; Azzalin, G.; Carucci, N.; Tavolaro, S.; Castellano, L.; Magrelli, A.; Citarella, F.; Messina, M.; et al. Quantitative technologies establish a novel microRNA profile of chronic lymphocytic leukemia. Blood 2007, 109, 4944–4951. [Google Scholar] [CrossRef]
  28. Mraz, M.; Malinova, K.; Kotaskova, J.; Pavlova, S.; Tichy, B.; Malcikova, J.; Stano Kozubik, K.; Smardova, J.; Brychtova, Y.; Doubek, M.; et al. miR-34a, miR-29c and miR-17-5p are downregulated in CLL patients with TP53 abnormalities. Leukemia 2009, 23, 1159–1163. [Google Scholar] [CrossRef]
  29. Mraz, M.; Dolezalova, D.; Plevova, K.; Stano Kozubik, K.; Mayerova, V.; Cerna, K.; Musilova, K.; Tichy, B.; Pavlova, S.; Borsky, M.; et al. MicroRNA-650 expression is influenced by immunoglobulin gene rearrangement and affects the biology of chronic lymphocytic leukemia. Blood 2012, 119, 2110–2113. [Google Scholar] [CrossRef]
  30. Negrini, M.; Cutrona, G.; Bassi, C.; Fabris, S.; Zagatti, B.; Colombo, M.; Ferracin, M.; D’Abundo, L.; Saccenti, E.; Matis, S.; et al. microRNAome expression in chronic lymphocytic leukemia: Comparison with normal B-cell subsets and correlations with prognostic and clinical parameters. Clin. Cancer Res. 2014, 20, 4141–4153. [Google Scholar] [CrossRef]
  31. Calin, G.A.; Dumitru, C.D.; Shimizu, M.; Bichi, R.; Zupo, S.; Noch, E.; Aldler, H.; Rattan, S.; Keating, M.; Rai, K.; et al. Frequent deletions and down-regulation of micro- RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc. Natl. Acad. Sci. USA 2002, 99, 15524–15529. [Google Scholar] [CrossRef]
  32. Klein, U.; Lia, M.; Crespo, M.; Siegel, R.; Shen, Q.; Mo, T.; Ambesi-Impiombato, A.; Califano, A.; Migliazza, A.; Bhagat, G.; et al. The DLEU2/miR-15a/16-1 cluster controls B cell proliferation and its deletion leads to chronic lymphocytic leukemia. Cancer Cell 2010, 17, 28–40. [Google Scholar] [CrossRef] [PubMed]
  33. Cimmino, A.; Calin, G.A.; Fabbri, M.; Iorio, M.V.; Ferracin, M.; Shimizu, M.; Wojcik, S.E.; Aqeilan, R.I.; Zupo, S.; Dono, M.; et al. miR-15 and miR-16 induce apoptosis by targeting BCL2. Proc. Natl. Acad. Sci. USA 2005, 102, 13944–13949. [Google Scholar] [CrossRef]
  34. Cui, B.; Chen, L.; Zhang, S.; Mraz, M.; Fecteau, J.F.; Yu, J.; Ghia, E.M.; Zhang, L.; Bao, L.; Rassenti, L.Z.; et al. MicroRNA-155 influences B-cell receptor signaling and associates with aggressive disease in chronic lymphocytic leukemia. Blood 2014, 124, 546–554. [Google Scholar] [CrossRef]
  35. Mraz, M.; Chen, L.; Rassenti, L.Z.; Ghia, E.M.; Li, H.; Jepsen, K.; Smith, E.N.; Messer, K.; Frazer, K.A.; Kipps, T.J. miR-150 influences B-cell receptor signaling in chronic lymphocytic leukemia by regulating expression of GAB1 and FOXP1. Blood 2014, 124, 84–95. [Google Scholar] [CrossRef] [PubMed]
  36. Palacios, F.; Abreu, C.; Prieto, D.; Morande, P.; Ruiz, S.; Fernandez-Calero, T.; Naya, H.; Libisch, G.; Robello, C.; Landoni, A.I.; et al. Activation of the PI3K/AKT pathway by microRNA-22 results in CLL B-cell proliferation. Leukemia 2015, 29, 115–125. [Google Scholar] [CrossRef] [PubMed]
  37. Tili, E.; Michaille, J.J.; Luo, Z.; Volinia, S.; Rassenti, L.Z.; Kipps, T.J.; Croce, C.M. The down-regulation of miR-125b in chronic lymphocytic leukemias leads to metabolic adaptation of cells to a transformed state. Blood 2012, 120, 2631–2638. [Google Scholar] [CrossRef]
  38. Cerna, K.; Oppelt, J.; Chochola, V.; Musilova, K.; Seda, V.; Pavlasova, G.; Radova, L.; Arigoni, M.; Calogero, R.A.; Benes, V.; et al. MicroRNA miR-34a downregulates FOXP1 during DNA damage response to limit BCR signalling in chronic lymphocytic leukaemia B cells. Leukemia 2019, 33, 403–414. [Google Scholar] [CrossRef] [PubMed]
  39. Guinn, D.; Ruppert, A.S.; Maddocks, K.; Jaglowski, S.; Gordon, A.; Lin, T.S.; Larson, R.; Marcucci, G.; Hertlein, E.; Woyach, J.; et al. miR-155 expression is associated with chemoimmunotherapy outcome and is modulated by Bruton’s tyrosine kinase inhibition with Ibrutinib. Leukemia 2015, 29, 1210–1213. [Google Scholar] [CrossRef]
  40. Guinn, D.; Lehman, A.; Fabian, C.; Yu, L.; Maddocks, K.; Andritsos, L.A.; Jones, J.A.; Flynn, J.M.; Jaglowski, S.M.; Woyach, J.A.; et al. The regulation of tumor-suppressive microRNA, miR-126, in chronic lymphocytic leukemia. Cancer Med. 2017, 6, 778–787. [Google Scholar] [CrossRef]
  41. Visone, R.; Rassenti, L.Z.; Veronese, A.; Taccioli, C.; Costinean, S.; Aguda, B.D.; Volinia, S.; Ferracin, M.; Palatini, J.; Balatti, V.; et al. Karyotype-specific microRNA signature in chronic lymphocytic leukemia. Blood 2009, 114, 3872–3879. [Google Scholar] [CrossRef]
  42. Van Roosbroeck, K.; Calin, G.A. MicroRNAs in chronic lymphocytic leukemia: miRacle or miRage for prognosis and targeted therapies? Semin. Oncol. 2016, 43, 209–214. [Google Scholar] [CrossRef] [PubMed]
  43. Papakonstantinou, N.; Ntoufa, S.; Chartomatsidou, E.; Papadopoulos, G.; Hatzigeorgiou, A.; Anagnostopoulos, A.; Chlichlia, K.; Ghia, P.; Muzio, M.; Belessi, C.; et al. Differential microRNA profiles and their functional implications in different immunogenetic subsets of chronic lymphocytic leukemia. Mol. Med. 2013, 19, 115–123. [Google Scholar] [CrossRef]
  44. Van Roosbroeck, K.; Bayraktar, R.; Calin, S.; Bloehdorn, J.; Dragomir, M.P.; Okubo, K.; Bertilaccio, M.T.S.; Zupo, S.; You, M.J.; Gaidano, G.; et al. The involvement of microRNA in the pathogenesis of Richter syndrome. Haematologica 2019, 104, 1004–1015. [Google Scholar] [CrossRef]
  45. Balatti, V.; Tomasello, L.; Rassenti, L.Z.; Veneziano, D.; Nigita, G.; Wang, H.Y.; Thorson, J.A.; Kipps, T.J.; Pekarsky, Y.; Croce, C.M. miR-125a and miR-34a expression predicts Richter syndrome in chronic lymphocytic leukemia patients. Blood 2018, 132, 2179–2182. [Google Scholar] [CrossRef]
  46. Scandurra, M.; Rossi, D.; Deambrogi, C.; Rancoita, P.M.; Chigrinova, E.; Mian, M.; Cerri, M.; Rasi, S.; Sozzi, E.; Forconi, F.; et al. Genomic profiling of Richter’s syndrome: Recurrent lesions and differences with de novo diffuse large B-cell lymphomas. Hematol. Oncol. 2010, 28, 62–67. [Google Scholar] [CrossRef]
  47. Condoluci, A.; Rossi, D. Richter Syndrome. Curr. Oncol. Rep. 2021, 23, 26. [Google Scholar] [CrossRef] [PubMed]
  48. Calin, G.A.; Cimmino, A.; Fabbri, M.; Ferracin, M.; Wojcik, S.E.; Shimizu, M.; Taccioli, C.; Zanesi, N.; Garzon, R.; Aqeilan, R.I.; et al. MiR-15a and miR-16-1 cluster functions in human leukemia. Proc. Natl. Acad. Sci. USA 2008, 105, 5166–5171. [Google Scholar] [CrossRef]
  49. Cutrona, G.; Matis, S.; Colombo, M.; Massucco, C.; Baio, G.; Valdora, F.; Emionite, L.; Fabris, S.; Recchia, A.G.; Gentile, M.; et al. Effects of miRNA-15 and miRNA-16 expression replacement in chronic lymphocytic leukemia: Implication for therapy. Leukemia 2017, 31, 1894–1904. [Google Scholar] [CrossRef]
  50. Matis, S.; Grazia Recchia, A.; Colombo, M.; Cardillo, M.; Fabbi, M.; Todoerti, K.; Bossio, S.; Fabris, S.; Cancila, V.; Massara, R.; et al. MiR-146b-5p regulates IL-23 receptor complex expression in chronic lymphocytic leukemia cells. Blood Adv. 2022, 6, 5593–5612. [Google Scholar] [CrossRef] [PubMed]
  51. Lionetti, M.; Barbieri, M.; Favasuli, V.; Taiana, E.; Fabris, S.; Favoino, C.; Ciceri, G.; Matis, S.; Colombo, M.; Massara, R.; et al. Frequency and clinical relevance of coding and noncoding NOTCH1 mutations in early stage Binet A chronic lymphocytic leukemia patients. Hematol. Oncol. 2020, 38, 406–408. [Google Scholar] [CrossRef] [PubMed]
  52. Monti, P.; Lionetti, M.; De Luca, G.; Menichini, P.; Recchia, A.G.; Matis, S.; Colombo, M.; Fabris, S.; Speciale, A.; Barbieri, M.; et al. Time to first treatment and P53 dysfunction in chronic lymphocytic leukaemia: Results of the O-CLL1 study in early stage patients. Sci. Rep. 2020, 10, 18427. [Google Scholar] [CrossRef] [PubMed]
  53. Morabito, F.; Tripepi, G.; Moia, R.; Recchia, A.G.; Boggione, P.; Mauro, F.R.; Bossio, S.; D’Arrigo, G.; Martino, E.A.; Vigna, E.; et al. Lymphocyte Doubling Time As A Key Prognostic Factor To Predict Time To First Treatment In Early-Stage Chronic Lymphocytic Leukemia. Front. Oncol. 2021, 11, 684621. [Google Scholar] [CrossRef] [PubMed]
  54. Morabito, F.; Adornetto, C.; Monti, P.; Amaro, A.; Reggiani, F.; Colombo, M.; Rodriguez-Aldana, Y.; Tripepi, G.; D’Arrigo, G.; Vener, C.; et al. Genes selection using deep learning and explainable artificial intelligence for chronic lymphocytic leukemia predicting the need and time to therapy. Front. Oncol. 2023, 13, 1198992. [Google Scholar] [CrossRef]
  55. Herling, C.D.; Cymbalista, F.; Gross-Ophoff-Muller, C.; Bahlo, J.; Robrecht, S.; Langerbeins, P.; Fink, A.M.; Al-Sawaf, O.; Busch, R.; Porcher, R.; et al. Early treatment with FCR versus watch and wait in patients with stage Binet A high-risk chronic lymphocytic leukemia (CLL): A randomized phase 3 trial. Leukemia 2020, 34, 2038–2050. [Google Scholar] [CrossRef]
  56. Muchtar, E.; Kay, N.E.; Parikh, S.A. Early intervention in asymptomatic chronic lymphocytic leukemia. Clin. Adv. Hematol. Oncol. 2021, 19, 92–103. [Google Scholar]
  57. Gentile, M.; Shanafelt, T.D.; Mauro, F.R.; Laurenti, L.; Rossi, D.; Molica, S.; Vincelli, I.; Cutrona, G.; Uccello, G.; Pepe, S.; et al. Comparison between the CLL-IPI and the Barcelona-Brno prognostic model: Analysis of 1299 newly diagnosed cases. Am. J. Hematol. 2018, 93, E35–E37. [Google Scholar] [CrossRef] [PubMed]
  58. Gentile, M.; Shanafelt, T.D.; Rossi, D.; Laurenti, L.; Mauro, F.R.; Molica, S.; Cutrona, G.; Uccello, G.; Campanelli, M.; Vigna, E.; et al. Validation of the CLL-IPI and comparison with the MDACC prognostic index in newly diagnosed patients. Blood 2016, 128, 2093–2095. [Google Scholar] [CrossRef]
  59. Gentile, M.; Shanafelt, T.D.; Cutrona, G.; Molica, S.; Tripepi, G.; Alvarez, I.; Mauro, F.R.; Di Renzo, N.; Di Raimondo, F.; Vincelli, I.; et al. A progression-risk score to predict treatment-free survival for early stage chronic lymphocytic leukemia patients. Leukemia 2016, 30, 1440–1443. [Google Scholar] [CrossRef]
  60. Fabbri, M.; Bottoni, A.; Shimizu, M.; Spizzo, R.; Nicoloso, M.S.; Rossi, S.; Barbarotto, E.; Cimmino, A.; Adair, B.; Wojcik, S.E.; et al. Association of a microRNA/TP53 feedback circuitry with pathogenesis and outcome of B-cell chronic lymphocytic leukemia. JAMA 2011, 305, 59–67. [Google Scholar] [CrossRef]
  61. Sharma, S.; Pavlasova, G.M.; Seda, V.; Cerna, K.A.; Vojackova, E.; Filip, D.; Ondrisova, L.; Sandova, V.; Kostalova, L.; Zeni, P.F.; et al. miR-29 modulates CD40 signaling in chronic lymphocytic leukemia by targeting TRAF4: An axis affected by BCR inhibitors. Blood 2021, 137, 2481–2494. [Google Scholar] [CrossRef] [PubMed]
  62. Wu, L.; Li, K.; Lin, W.; Liu, J.; Qi, Q.; Shen, G.; Chen, W.; He, W. Long noncoding RNA LINC01291 promotes the aggressive properties of melanoma by functioning as a competing endogenous RNA for microRNA-625-5p and subsequently increasing IGF-1R expression. Cancer Gene Ther. 2022, 29, 341–357. [Google Scholar] [CrossRef] [PubMed]
  63. Kumari, R.; Roy, U.; Desai, S.; Nilavar, N.M.; Van Nieuwenhuijze, A.; Paranjape, A.; Radha, G.; Bawa, P.; Srivastava, M.; Nambiar, M.; et al. MicroRNA miR-29c regulates RAG1 expression and modulates V(D)J recombination during B cell development. Cell Rep. 2021, 36, 109390. [Google Scholar] [CrossRef]
  64. Li, L.; Shou, H.; Wang, Q.; Liu, S. Investigation of the potential theranostic role of KDM5B/miR-29c signaling axis in paclitaxel resistant endometrial carcinoma. Gene 2019, 694, 76–82. [Google Scholar] [CrossRef] [PubMed]
  65. Cutrona, G.; Tripodo, C.; Matis, S.; Recchia, A.G.; Massucco, C.; Fabbi, M.; Colombo, M.; Emionite, L.; Sangaletti, S.; Gulino, A.; et al. Microenvironmental regulation of the IL-23R/IL-23 axis overrides chronic lymphocytic leukemia indolence. Sci. Transl. Med. 2018, 10, eaal1571. [Google Scholar] [CrossRef]
  66. Vila-Casadesus, M.; Gironella, M.; Lozano, J.J. MiRComb: An R Package to Analyse miRNA-mRNA Interactions. Examples across Five Digestive Cancers. PLoS ONE 2016, 11, e0151127. [Google Scholar] [CrossRef] [PubMed]
  67. Huang, H.Y.; Lin, Y.C.; 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]
  68. Wu, T.; Hu, E.; Xu, S.; Chen, M.; Guo, P.; Dai, Z.; Feng, T.; Zhou, L.; Tang, W.; Zhan, L.; et al. clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation 2021, 2, 100141. [Google Scholar] [CrossRef]
Ncrna 10 00046 g001
Figure 1. Kaplan–Meier curve of TTFT in the 224 CLL cases from our prospective O-CLL1 study.
Figure 1. Kaplan–Meier curve of TTFT in the 224 CLL cases from our prospective O-CLL1 study.
Ncrna 10 00046 g001
Ncrna 10 00046 g002
Figure 2. Risk score histogram distribution of 224 cases from O-CLL1 study based on the 16 miRNAs significantly associated with TTFT in multivariate Cox regression model.
Figure 2. Risk score histogram distribution of 224 cases from O-CLL1 study based on the 16 miRNAs significantly associated with TTFT in multivariate Cox regression model.
Ncrna 10 00046 g002
Ncrna 10 00046 g003
Figure 3. Dot plot displaying the magnitude of negative correlation between the 16 miRNAs significantly associated with TTFT in the 224 CLL cases from our prospective O-CLL1 study and the 50 AI genes from the same cohort. AI genes are sorted by significance magnitude from the bottom to the top. The first 10 genes represent the top AI-based selected genes that were significantly associated with TTFT. The greater intensity of the color in the circle shapes corresponds to a higher magnitude of the negative correlation measured by a Spearman test. rho (or Spearman’s correlation coefficient, measured using a Spearman test) is a non-parametric measure that evaluates the strength and direction of a monotonic relationship between two ordinal or quantitative variables using the ranks of the observations instead of their absolute values.
Figure 3. Dot plot displaying the magnitude of negative correlation between the 16 miRNAs significantly associated with TTFT in the 224 CLL cases from our prospective O-CLL1 study and the 50 AI genes from the same cohort. AI genes are sorted by significance magnitude from the bottom to the top. The first 10 genes represent the top AI-based selected genes that were significantly associated with TTFT. The greater intensity of the color in the circle shapes corresponds to a higher magnitude of the negative correlation measured by a Spearman test. rho (or Spearman’s correlation coefficient, measured using a Spearman test) is a non-parametric measure that evaluates the strength and direction of a monotonic relationship between two ordinal or quantitative variables using the ranks of the observations instead of their absolute values.
Ncrna 10 00046 g003
Ncrna 10 00046 g004
Figure 4. Interaction analysis between miRNAs and genes which were found to be associated with TTFT. Interaction networks between has-miR-29c-3p with 4 AI genes ANKRD52, KDM5B, PRICKLE1, and ZBTB34 or hsa-miR-625-5p with IGF1R are shown. Arrows, circles, and squares indicate the miRNA–mRNA interaction scores, miRNA, and mRNA log-transformed fold-change (Log FC), respectively. Log FC values for miRNA (◯) and mRNA (□) are indicated by a color gradient scale. Regulation score parameter values related to the IGHV mutational status (arrows) are highlighted by a scale color gradient from light gray that represents a midpoint value of 0 to red that indicates positive values (interaction score).
Figure 4. Interaction analysis between miRNAs and genes which were found to be associated with TTFT. Interaction networks between has-miR-29c-3p with 4 AI genes ANKRD52, KDM5B, PRICKLE1, and ZBTB34 or hsa-miR-625-5p with IGF1R are shown. Arrows, circles, and squares indicate the miRNA–mRNA interaction scores, miRNA, and mRNA log-transformed fold-change (Log FC), respectively. Log FC values for miRNA (◯) and mRNA (□) are indicated by a color gradient scale. Regulation score parameter values related to the IGHV mutational status (arrows) are highlighted by a scale color gradient from light gray that represents a midpoint value of 0 to red that indicates positive values (interaction score).
Ncrna 10 00046 g004
Table 1. Cox univariable analyses of miRNAs significantly associated with TTFT.
Table 1. Cox univariable analyses of miRNAs significantly associated with TTFT.
miRNA-IDUnits of
Increase 1		HR95%CI
Lower Limit		95% CI
Upper Limit		p-Value
1-3p11.1661.0061.3520.041
103a-3p10.9980.9961.0000.05
106a-5p11.4381.0032.0620.048
10b-3p11.2391.0551.4560.009
1224-5p11.0641.0241.1050.002
1225-5p1001.0621.0011.1260.046
124-3p11.4581.1781.805<0.001
125b-5p10.7440.560.9890.042
138-5p10.5990.40.8960.013
140-3p10.9910.9860.9960.001
144-3p11.0041.0021.0060.002
144-5p11.0241.0081.0410.003
146b-5p10.9910.9830.9980.019
148a-3p11.0061.0021.0100.007
150-5p10000.9250.8580.9970.041
150-3p11.0351.0091.0620.008
151-3p10.9710.9460.9960.026
151-5p10.9960.9930.9980.001
155-5p1001.0581.0141.1040.009
15a-5p1001.0741.0341.117<0.001
18411.4721.1261.9250.005
193a-3p11.1321.0261.2500.014
20a-3p11.0491.0071.0930.022
21-5p10001.1521.0561.2570.001
222-3p10.9460.9070.9880.012
223-5p10.8190.7120.9430.006
24-1-5p11.3121.0131.6980.04
26a-5p1000.9010.8130.9990.047
28-5p11.0051.0011.0090.012
296-3p10.5700.3480.9340.026
29811.3071.0221.6700.033
29c-3p1000.9510.9230.980<0.001
29c-5p10.9520.9240.9820.002
301a-3p11.0381.0051.0730.024
30c-5p1000.6210.4000.9620.033
323-3p10.6460.4260.9790.04
338-5p10.7850.6430.9600.018
339-3p10.6890.5210.9100.009
33a-3p10.3750.2200.639<0.001
361-3p10.9890.9780.9990.034
37011.0581.0281.088<0.001
371-5p11.0821.0161.1520.014
373-5p11.0711.0061.140.031
376b-3p11.5231.0462.2180.028
491-3p11.7661.3312.343<0.001
500-3p10.8170.7170.9310.002
502-3p10.8720.7960.9550.003
502-5p10.6730.5040.8980.007
513a-5p11.0071.0021.0120.008
518c-5p11.1481.0331.2760.01
520b11.2341.0671.4260.005
532-3p10.8980.8410.9590.001
532-5p10.9390.8910.9890.018
55211.5561.0142.3880.043
55711.1971.0991.303<0.001
56611.6161.1882.1970.002
574-3p11.0301.0061.0550.015
582-3p10.4650.2740.7890.005
584-5p11.1621.0221.320.022
59610.5970.3910.9130.017
60111.0691.011.1310.022
60311.5521.0232.3560.039
625-5p10.9600.9400.981<0.001
628-3p10.6300.4170.9520.028
63111.1801.0061.3850.042
64511.6041.0912.3580.016
659-3p11.1141.011.2280.03
66110.5790.3420.9810.042
66511.1451.0081.3000.037
671-5p11.0461.0141.0790.004
877-5p11.2451.0311.5030.023
9-3p11.0861.0151.1630.017
99a-5p10.6150.4210.8980.012
1 To provide a clinically meaningful magnitude of the effect of each miRNA on the study outcome, the units of increase (1, 100, or 1000) were chosen according to the data distribution of the variable.
Table 2. 16 miRNAs independently associated with TTFT after a multivariable Cox regression analysis using a backward elimination strategy.
Table 2. 16 miRNAs independently associated with TTFT after a multivariable Cox regression analysis using a backward elimination strategy.
miRNA-IDUnits of
Increase 1		HR95% CI
Lower Limit		95% CI
Upper Limit		p-Value
582-3p10.2780.1450.535<0.001
33a-3p10.3340.160.6970.003
516a-5p10.4900.2970.8100.005
99a-5p10.5120.3410.7690.001
296-3p10.5390.3010.9670.038
502-5p10.6230.430.9050.013
625-5p10.9580.9370.98<0.001
29c-3p1000.9360.9030.970<0.001
150-5p10001.1121.0051.2310.039
148a-3p11.0091.0041.014<0.001
28-5p11.011.0051.014<0.001
144-5p11.0491.0261.072<0.001
671-5p11.0751.0271.1250.002
1-3p11.2611.0471.5170.014
193a-3p11.3431.1861.52<0.001
124-3p11.5361.2331.913<0.001
1 To provide a clinically meaningful magnitude of the effect of each miRNA on the study outcome, the units of increase (1,100, or 1000) were chosen according to the data distribution of the variable.
Table 3. Prognostic performance of the basic and the expanded models.
Table 3. Prognostic performance of the basic and the expanded models.
Basic ModelExpanded Model
Harrell’s C-index75.0%81.1%
Explained variation in TTFT45.4%63.3%
IDI 1-14.9%, p < 0.001
NRI 2-44.2%, p < 0.001
1 IDI: integrated discrimination improvement; 2 NRI: net reclassification index.
Table 4. Summary of significantly enriched ontologies in the 16 analyzed miRNAs displaying ontology-associated miRNAs and the corresponding adjusted p-values and q-values.
Table 4. Summary of significantly enriched ontologies in the 16 analyzed miRNAs displaying ontology-associated miRNAs and the corresponding adjusted p-values and q-values.
EnrichmentIDTerm SpecificationAssociated
miRNA Genes		Adjusted
p-Value		q-Value
KEGG ORA summaryhsa05206MicroRNAs in cancerMIR1-2, MIR28, MIR29C, MIR99A, MIR124-3, MIR150, MIR6251.41 × 10−107.41 × 10−10
WikiPathways ORA summaryWP299Nuclear receptors in lipid metabolism and toxicityMIR33A0.030.009
WP430Statin inhibition of cholesterol productionMIR33A0.030.009
WP1545miRNAs involved in DNA damage responseMIR29C0.040.01
WP1601Fluoropyrimidine activityMIR29C0.030.009
WP2023Cell differentiation expanded indexMIR1500.040.01
WP2249Metastatic brain tumorMIR29C0.030.009
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

Nano, E.; Reggiani, F.; Amaro, A.A.; Monti, P.; Colombo, M.; Bertola, N.; Ferrero, F.; Fais, F.; Bruzzese, A.; Martino, E.A.; et al. MicroRNA Profiling as a Predictive Indicator for Time to First Treatment in Chronic Lymphocytic Leukemia: Insights from the O-CLL1 Prospective Study. Non-Coding RNA 2024, 10, 46. https://doi.org/10.3390/ncrna10050046

AMA Style

Nano E, Reggiani F, Amaro AA, Monti P, Colombo M, Bertola N, Ferrero F, Fais F, Bruzzese A, Martino EA, et al. MicroRNA Profiling as a Predictive Indicator for Time to First Treatment in Chronic Lymphocytic Leukemia: Insights from the O-CLL1 Prospective Study. Non-Coding RNA. 2024; 10(5):46. https://doi.org/10.3390/ncrna10050046

Chicago/Turabian Style

Nano, Ennio, Francesco Reggiani, Adriana Agnese Amaro, Paola Monti, Monica Colombo, Nadia Bertola, Fabiana Ferrero, Franco Fais, Antonella Bruzzese, Enrica Antonia Martino, and et al. 2024. "MicroRNA Profiling as a Predictive Indicator for Time to First Treatment in Chronic Lymphocytic Leukemia: Insights from the O-CLL1 Prospective Study" Non-Coding RNA 10, no. 5: 46. https://doi.org/10.3390/ncrna10050046

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