Competing endogenous RNA networks related to prognosis in chronic lymphocytic leukemia: comprehensive analyses and construction of a novel risk score model
Biomarker Research volume 10, Article number: 75 (2022)
Chronic lymphocytic leukemia (CLL) is a heterogeneous B-cell malignancy that lacks specific biomarkers and drug targets. Competing endogenous RNAs (ceRNAs) play vital roles in oncogenesis and tumor progression by sponging microRNAs (miRNAs). Nevertheless, the regulatory mechanisms of survival-related ceRNA networks in CLL remain to be uncovered.
We included 865 de novo CLL patients to investigate RNA expression profiles and Illumina sequencing was performed on four CLL patients, two CLL cell lines and six healthy donors in our center. According to univariate Cox regression, LASSO regression as well as multivariate Cox regression analyses, we established a novel risk score model in CLL patients. Immune signatures were compared between the low- and high-risk groups with CIBERSORT and ESTIMATE program. Afterwards, we analyzed the relationship between differentially expressed miRNAs (DEmiRNAs) and IGHV mutational status, p53 mutation status and del17p. Based on the survival analyses and differentially expressed RNAs with targeting relationships, the lncRNA/circRNA-miRNA-mRNA ceRNA networks were constructed. In addition, the circRNA circ_0002078/miR-185-3p/TCF7L1 axis was verified and their interrelations were delineated by dual-luciferase reporter gene assay.
Totally, 57 differentially expressed mRNAs (DEmRNAs) and 335 DEmiRNAs were identified between CLL patient specimens and normal B cells. A novel risk score model consisting of HTN3, IL3RA and NCK1 was established and validated. The concordance indexes of the model were 0.825, 0.719 and 0.773 in the training, test and total sets, respectively. The high-risk group was related to del(13q14) as well as shorter overall survival (OS). Moreover, we identified DEmiRNAs that related to cytogenetic abnormality of CLL patients, which revealed that miR-324-3p was associated with IGHV mutation, p53 mutation and del17p. The survival-related lncRNA/circRNA-miRNA-mRNA ceRNA networks were constructed to further facilitate the development of potential predictive biomarkers. Besides, the expression of circ_0002078 and TCF7L1 were significantly elevated and miR-185-3p was obviously decreased in CLL patients. Circ_0002078 regulated TCF7L1 expression by competing with TCF7L1 for miR-185-3p.
The comprehensive analyses of RNA expression profiles provide pioneering insights into the molecular mechanisms of CLL. The novel risk score model and survival-related ceRNA networks promote the development of prognostic biomarkers and potential therapeutic vulnerabilities for CLL.
As the most prevalent type of leukemia among adults in the western world , chronic lymphocytic leukemia (CLL) is defined as a monoclonal lymphoproliferative disease featuring the gradual cumulation of immunologically dysfunctional mature B lymphocytes [2, 3]. Even if there have been great advances in the diagnosis and treatment of CLL, the underlying pathogenesis of CLL has not been completely understood until now. With the development of RNA-sequencing, a myriad of potential biomarkers in association with the prognosis of CLL have been discovered. However, we still need to explore risk score models with highly robust predictive power for survival to provide a basis for individualized treatment of CLL patients.
Lately, accumulative evidence revealed that non-coding RNAs (ncRNAs) played pivotal roles in regulating gene expression [4, 5] and can affect the cancer process by regulating alternative splicing [6, 7]. NcRNAs, which are a group of RNAs that do not have protein coding function, contain many categories. Among them, microRNA (miRNA), long noncoding RNA (lncRNA) and circular RNA (circRNA) are mainly involved in post-transcriptional regulation .
Competing endogenous RNAs (ceRNAs) hypothesis was first proposed by the team of Pier Paolo Pandolfi, which assumed that if lncRNAs or other types of RNA molecules have the same miRNA response elements (MREs) as mRNA, they can act as ceRNAs to compete with mRNA for binding to miRNAs, accordingly influencing the regulation of miRNAs on its target genes . This hypothesis has been widely accepted and emerging evidence has showed the roles of ceRNA networks in tumorigenesis, invasion and metastasis [10, 11]. For example, lncRNA TMPO-AS1/miR-126-5p/BRCC3 axis facilitates gastric carcinoma development and angiogenesis . LncRNA APCDD1L-AS1 can suppress autophagic degradation of EGFR through miR-1322/miR-1972/miR-324-3p-SIRT5 axis, thereby contributes to icotinib resistance . The inhibition of lncRNA LAMP5-AS1 dramatically decreases colony formation and enhances primary MLL leukemia CD34+ cells differentiation . Wang et al. found that in hepatocellular carcinoma, circ_0001588/miR-874/CDK4 aixs promoted cell proliferation, invasion, and migration . CircRNA circ_0025033/miR-184/LSM4 axis was shown to accelerate the development of ovarian cancer . In colorectal cancer, circHERC4 inactivates miR-556-5p, resulting in the activation of the CTBP2/E-cadherin pathway, which promotes cancer metastasis .
In CLL, the functions of certain lncRNA/circRNA axes in the development of CLL have been elucidated. Circ-CBFB was highly expressed in CLL and promoted FZD3 expression by inhibiting miR-607, resulting in the activation of the Wnt/β-catenin pathway and subsequent CLL progression . Circ_0132266 was found to accelerate apoptosis and inhibit proliferation by miR-337-3p/PML pathway . In addition, Li et al. found that the relative expression of circ-RPL15 in CLL patients was up-regulated and was positively associated with IGHV mutation status. Circ-RPL15 can sponge miR-146b-3p and block its suppression of the RAS/RAF1/MEK/ERK pathway . LncRNA CRNDE was down-regulated in CLL cell lines, and its protective effect in CLL was to impair proliferation through the miR-28/NDRG2 axis . Yet, the roles of the ceRNA networks in CLL and their regulatory mechanism as well as more meaningful biomarkers still require further exploration.
In the present study, distinct RNA expression signatures were identified in CLL patients and CLL cell lines. Differentially expressed genes (DEGs) in association with the overall survival (OS) were revealed, and a novel risk score model including HTN3, IL3RA and NCK1 was established. Besides, we constructed survival-related ceRNA networks and validated the circRNA circ_0002078/miR-185-3p/TCF7L1 axis. Our results provide an in-depth understanding of the underlying mechanism of CLL pathogenesis. Moreover, the current study is helpful for uncovering novel prognostic biomarkers and therapeutic vulnerabilities of CLL.
Materials and methods
In this study, 865 de novo CLL patients were included. The CLL patients for Illumina sequencing were collected from Shandong Provincial Hospital CLL (SPHCLL) database and were diagnosed on the basis of the revised International Workshop on Chronic Lymphocytic Leukemia criteria . B cells of these samples were separated according to the previously reported method . Ethics approval was acquired from the Medical Ethical Committee of Shandong Provincial Hospital, and all of the patients and healthy donors have given their informed consent.
RNA sequence analyses and differentially expressed analyses
The lncRNA, circRNA, miRNA and mRNA expression data of normal B cells, CLL patient specimens and CLL cells were available from Illumina sequencing. Compare CLL patient specimens and CLL cells with normal B cells, we identified the differentially expressed mRNAs (DEmRNAs), differentially expressed miRNAs (DEmiRNAs), differentially expressed lncRNAs (DElncRNAs) and differentially expressed circRNAs (DEcircRNAs) with the statistical threshold of |log2fold change|> 2 as well as p values < 0.05.
Functional enrichment analyses
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to analyze the molecular function (MF), cell component (CC), biological processes (BP) and signaling pathways of DEmRNAs. P values < 0.05 were regarded as statistically significant.
Construction of the novel risk score model
Univariate Cox regression analysis was applied to initially filter DEGs in association with OS. LASSO regression analysis was conducted on these DEGs to avoid overfitting. Multivariate Cox regression method was used to establish the risk score model. The calculation method of the risk scoring system is as follows: Risk score = coefficient1 * gene1 expression value + coefficient2 * gene2 expression value + coefficienti * genei expression. The above analyses were conducted by R software.
Exploration of immune signatures between low-risk and high-risk groups
In order to explore the distinct immune signatures between low-risk and high-risk groups, we performed the CIBERSORT algorithm. The ESTIMATE method was conducted to compute the scores of immune and stromal cells as well as ESTIMATE score in low-risk and high-risk groups.
Construction of lncRNA/circRNA-miRNA-mRNA ceRNA networks
We used miRanda and psRobot software to constructed ceRNA regulatory networks with DElncRNA/DEcircRNA as decoy, DEmiRNA as the core, and DEmRNA as the target according to the RNA sequence results. Cytoscape 3.7.1 was used to visualize and select the lncRNA/circRNA-miRNA-mRNA ceRNA networks on the basis of the DEGs that associated with OS.
Cell lines and cell culture
The cell lines for RNA-sequencing were MEC1 cell line and EHEB cell line. The MEC1 cells were obtained from University of California, San Diego. The EHEB cells and 293 T cells were purchased from American type culture collection. IMDM with 10% fetal bovine serum (FBS) was the medium for culturing the MEC1 cells. The EHEB cells and 293 T cells were cultured in RPMI1640 and DMEM respectively with FBS (10%). Cell culture conditions are 37 °C, 5% CO2 and saturated humidity.
Quantitative real-time polymerase chain reaction (qRT-PCR)
After the ceRNA networks were established, we validated three DEGs (TRIM34, SLC30A10, HOXD4), circ_0007675 and the circ_0002078/miR-185-3p/TCF7L1 axis. Peripheral blood mononuclear cells of patients with CLL and healthy donors were extracted following the previous method [24, 25]. Total RNA was obtained by RNAiso Plus (TaKaRa, Dalian, China). Then reverse transcription reagents (Accurate Biology, Dalian, China) were used for reverse transcription reaction to obtain cDNA. The qRT-PCR was performed with SyberGreen (Accurate Biology, Changsha, China) in LightCycler 480II (Roche, Basel, Swizerland). Internal controls were U6 and GAPDH. Biosune (Shanghai, China) synthesized all of the primers. Relevant primer sequences are shown in Supplementary table S9.
Dual-luciferase reporter gene assay
The interrelations between circ_0002078 and miR-185-3p, miR-185-3p and TCF7L1 were confirmed through the dual-luciferase reporter gene assay. PmirGLO-circ_0002078-WT/Mut and pmirGLO-3′UTR TCF7L1-WT/Mut plasmids (constructed by Biosune (Shanghai, China)) were co-transfected with miR-185-3p mimics and miR-NC respectively into the 293 T cells using lipofectamine 3000 (Invitrogen, USA). Dual-Luciferase Assay System Kit (E1910, Promega, USA) was applied to detect the relative luciferase activity after 48 h.
Cell proliferation assays
Cell proliferation was assessed with cell Counting Kit-8 (CCK-8; YEASEN, Shanghai, China). A total of 1 × 104 MEC1 or EHEB cells (each well) were seeded in a 96-well culture plate. After 24, 48, and 72 h of incubation, 10 μl of CCK-8 solution was added to each well. The absorbance of each well was measured at 450 nm.
Analysis of cell apoptosis and cell cycle
PE Annexin V Apoptosis Detection Kit I (BD Pharmingen, Bedford, MA, USA) was used for the analysis of cell apoptosis. After collecting treated cells, resuspended the cells and added 5 μl each of Annexin V-PE and 7AAD and incubated the mixture at room temperature for 15 min. Cells for cell cycle analysis were soaked in 70% ethanol overnight at 4 °C and then incubated with PI/ RNase Staining Buffer (BD Biosciences, Bedford, MA, USA) for 15 min. Navios Flow Cytometer (Beckman Coulter, CA, USA) was used for cell apoptosis and cell cycle analyses.
Kaplan–Meier (K-M) method was performed to generate survival curves. For the purpose of quantitatively evaluating the prognosis of CLL, we used the “rms” R package to establish the prognostic DEGs nomogram. The relevance between risk score model and cytogenetic abnormalities were explored by Chi-square test. The relationship between gene expression and clinical or other characteristics were identified by T-test. the relationship between circ_0002078 and TCF7L1 was explored by Pearson correlation analysis. P values < 0.05 were considered to be statistically significant.
Differentially expressed RNAs between CLL patients or CLL cells and normal B cells
The DElncRNAs, DEcircRNAs, DEmiRNAs and DEmRNAs between CLL patient specimens, CLL cells and normal B cells were acquired from the Illumina sequencing. Totally, we identified 57 DEmRNAs and 335 DEmiRNAs by comparing CLL patient specimens with normal B cells. At the same time, between CLL cells and normal B cells, 482 mRNAs and 302 miRNAs were showed differential expression. DEmRNAs and DEmiRNAs in the samples were shown in the heatmap (Fig. 1A-B). The top 10 up-regulation and down-regulation differentially expressed RNAs were presented in Supplementary table S1-8.
Functional enrichment analyses of DEGs
The potential functions of DEGs were explored by GO and KEGG analyses. GO analysis results demonstrated that changes in MF of DEGs between CLL patients and healthy controls were significantly focused on heterocyclic compound binding, organic cyclic compound binding and sequence-specific DNA binding (Fig. 1C). Meanwhile, the changes in MF of DEGs between MEC1 and EHEB cell lines and controls were significantly focused on sequence-specific DNA binding, oxidoreductase activity and cofactor binding (Supplementary Fig. 1A). Changes in CC of DEGs between CLL patient specimens and healthy donors were principally involved in organelle part, chromatin and intracellular part (Fig. 1C), meanwhile, between CLL cells and healthy controls they were involved in Golgi membrane, ribosomal subunit and Golgi apparatus part (Supplementary Fig. 1A). Changes in BP of DEGs between CLL patients and healthy donors primarily included regulation of gene expression, regulation of macromolecule metabolic process and regulation of cellular metabolic process (Fig. 1C), which between CLL cells and healthy controls primarily included cation transmembrane transport, lipid metabolic process and negative regulation of RNA metabolic process (Supplementary Fig. 1A).
For KEGG pathway analysis, the DEGs between CLL specimens and healthy donors were primarily relevant to Notch signaling pathway, JAK-STAT signaling pathway and cGMP-PKG signaling pathway (Fig. 1D). Additionally, between CLL cell lines and normal groups, DEGs were significantly associated with mTOR signaling pathway, NF-kappaB signaling pathway and cell cycle (Supplementary Fig. 1B).
Construction and validation of a novel risk score model
In this study, in order to reduce ethnic differences and the limitation of small samples, we studied the relationship between DEGs and OS in the genome microarray map GSE22762 to identify key genes in association with the prognosis of CLL patients. The dataset was divided into two groups: training set and test set. The training set was used for establishing novel risk score model. We used test and total sets to verify the stability and predictive power of the risk score model. According to the univariate Cox regression analysis, 8 genes were obviously in association with OS in the training set (Table 1). To avoid overfitting, we performed LASSO regression to further analyze these 8 survival-related genes (Fig. 2A, B). Multivariate Cox regression finally screened out 3 survival-related genes (HTN3, IL3RA and NCK1) for risk model establishment. Among them, IL3RA and NCK1 were independent predictors (Table 1; Fig. 2 C). Subsequently, we drew a nomogram consisting HTN3, IL3RA and NCK1 to quantitatively predict 1-, 2- and 3-year OS in CLL patients (Fig. 2D). The calculation method of the novel risk score model was as follows: 2.071 * HTN3 expression value + (-3.228) * IL3RA expression value + (-1.847) * NCK1 expression value. The risk score for each CLL patient in the training, test and total sets was calculated according to the above formular. In accordance with the median risk score of training set, we divided CLL patients into low- and high-risk group.
K-M survival curve indicated that OS in the low-risk group was longer than that in high-risk group in the training set (p = 0.010; Fig. 2H), test set (p = 0.026; Fig. 2M) and total set (p = 0.001; Fig. 2R). The effectiveness of the risk score were assessed by the receiver operating characteristic (ROC). The areas under the curve (AUC) of training set, test set and total set were 0.825, 0.719 and 0.773, respectively (Fig. 2I, N, S). The risk scores of CLL patients in the training set, test set and total set were ranked and their distribution were shown in Fig. 2E, J, O. The dot plots demonstrated the CLL patients’ survival status in the training, test and total sets (Fig. 2F, K, P). The expression of 3 survival-related DEGs between low- and high-risk groups were displayed in the heatmap (Fig. 2G, L, Q).
Immune signatures of risk groups and the correlation between risk groups of CLL patients with cytogenetic abnormality
We revealed the immune signatures of the low- and high-risk groups. Histogram showed that the proportion of neutrophils in the high-risk group is higher than that in the low-risk group (Fig. 3A, B). The correlation between immune cells in the two groups were shown in (Fig. 3C, D). The vioplot demonstrated the differential expression of 22 reported immune cell types between low- and high-risk groups. B cells naïve (p < 0.001), T cells CD8 (p = 0.005), CD4 memory resting T cells (p = 0.003), follicular helper T cells (p = 0.039), resting NK cells (p = 0.024), activated NK cells (p = 0.035), monocytes (p = 0.003), macrophages M0 (p = 0.003), macrophages M1 (p = 0.038), resting mast cells (p = 0.021), activated mast cells (p < 0.001) and eosinophils (p = 0.025) were obviously associated with risk group (Fig. 3E). The stromal score representing stroma feature of the tumor microenvironment in the high-risk group were lower than that of the low-risk group (p < 0.01). The immune score that representing immune cell signatures of high-risk group was lower compared to low-risk group, but there was no statistical significance (Fig. 3F). Moreover, distinct expression of immune checkpoints between the two groups indicated that BTLA (p < 0.001), CD200 (p < 0.01), CD27 (p < 0.001), CD40 (p < 0.05), LAG3 (p < 0.05) and TIGIT (p < 0.05) expression levels were evidently higher in high-risk group (Fig. 3G). Heatmap of immune function analysis demonstrated that the immune status of the high-risk group was more suppressed compared to the low-risk group (Fig. 3H). The effects of different immune functions on OS of CLL patients were shown in Supplementary Fig. 2. We subsequently investigated the relationship between the risk groups and cytogenetic abnormality by chi-square test in GSE25571. The results demonstrated that the high-risk group were related to del 13q14 and del 13q14 single (Table 2).
The correlation between DEmiRNAs and clinical characteristics
To explore the prognostic values of the DEmiRNAs, we analysis the relationship between DEmiRNAs and the mutational status of IGHV, p53 mutational status and del17p in the GSE40533 and GSE45328. In total, we identified 13 DEmiRNAs that related to IGHV mutational status. Among them, the high expression of miR-19b-1-5p (p < 0.05), miR-20a-3p (p < 0.05), miR-24-1-5p (p < 0.05), miR-127-3p (p < 0.05), miR-590-5p (p < 0.05), miR-340-3p (p < 0.001), miR-193a-3p (p < 0.001), miR-654-3p (p < 0.05) and miR-199a-5p (p < 0.05) were associated with IGHV unmutated status (Fig. 4A-I), at the same time, the high expression of miR-33a-3p (p < 0.05), miR-338-5p (p < 0.05), miR-139-3p (p < 0.05) and miR-324-3p (p < 0.01) were related to IGHV gene mutation (Fig. 4J-M). Additionally, the low expression of miR-30e-5p (p < 0.05), miR-324-3p (p < 0.05), miR-324-5p(p < 0.001) and miR-581 (p < 0.05) were related to p53 mutation (Fig. 4N-Q). We also demonstrated that the low expression of miR-30e-5p (p < 0.05), miR-324-5p (p < 0.001) and miR-324-3p (p < 0.01) were associated with del17p (Fig. 4R-T).
Prediction and construction of the ceRNA networks
We constructed the lncRNA-miRNA-mRNA ceRNA network according to the survival-related genes, DEmiRNAs and DElncRNA in CLL patients and healthy donors. Totally, 68 lncRNA nodes, 8 miRNA nodes and 3 mRNA nodes were selected to construct the ceRNA network (Fig. 5). To decipher the functions of DEcircRNAs in CLL patients, we constructed an abnormal circRNA‐miRNA‐mRNA ceRNA network. In total, 3 DEmRNAs, 8 DEmiRNAs and 101 DEcircRNAs were included in the ceRNA network (Fig. 6).
Validation of the circ_0002078/miR-185-3p/TCF7L1 axis
We verified the expression of three DEGs (TRIM34, SLC30A10, HOXD4) and circ_0007675 by qRT-PCR. The relative expression of three DEGs and circ_0007675 in CLL patients were remarkably up-regulated compared with normal B cells (p < 0.05; Supplementary Fig. 3A-D).
The expression of TCF7L1, miR-185-3p, and circ_0002078 in the CLL patients and healthy donors were detected respectively. The results showed that TCF7L1 as well as circ_0002078 were up-regulated (p < 0.05; Fig. 7A and C), whereas miR-185-3p was down-regulated in CLL patients compared to normal B cells (p < 0.05; Fig. 7B). Pearson correlation analysis revealed the strong linear correlation between TCF7L1 and circ_0002078 (r = 0.743, p < 0.001; Fig. 7D). K-M survival analysis in GSE39671 demonstrated that higher expression of TCF7L1 was correlated to shorter time to first treatment (TTFT) (p = 0.009; Fig. 7E). Based on the clinical information and expression levels in CLL patients, we performed the survival analysis of circ_0002078 and miR-185-3p. The high expression of circ_0002078 was associated with a shorter OS (p = 0.031; Fig. 7F). CLL patients with higher miR-185-3p expression had longer OS than patients with lower expression of miR-185-3p (p = 0.047; Fig. 7G).
The interrelationship among circ_0002078, miR-185-3p and TCF7L1 were detected by dual-luciferase reporter gene assay. When luciferase reporter plasmid containing circ_0002078-WT (or TCF7L1-WT) co-transfected with miR-185-3p mimics into 293 T cells, the luciferase activities were evidently decreased compared to that co-transfected with miR-NC. However, when the luciferase reporter plasmid containing circ_0002078-Mut (or TCF7L1-Mut) co-transfected with miR-185-3p mimics and miR-NC respectively into 293 T cells, there was no significant difference between the results of two groups (p < 0.01; Fig. 7H, I).
To explore the role of miR-185-3p on cell proliferation, we performed CCK-8 assay. We transferred miR-185-3p mimics and miR-NC into MEC1 cells. The expression of miR-185-3p was up-regulated in MEC1 cells (p < 0.01; Supplementary Fig. 4A). The OD values of miR-185-3p overexpression group were significantly lower than control group (Fig. 7J), which indicated that miR-185-3p inhibited cell proliferation of the CLL cells. We predicted the potential functions of target genes of circ_0002078 by GO and KEGG analyses. GO analysis results demonstrated that the target genes of circ_0002078 were enriched in positive regulation of gene expression, lymphocyte proliferation and regulation of B cell apoptotic process (Supplementary Fig. 4B). KEGG pathway analysis showed that the target genes were primarily relevant to JAK-STAT signaling pathway and cytokine-cytokine receptor interaction (Supplementary Fig. 4C). The results suggested that circ_0002078 may promote the occurrence and development of CLL in terms of cell proliferation and apoptosis. We transfected siRNA of circ_0002078 in EHEB cells and the expression of circ_0002078 was down-regulated (p < 0.001; Fig. 7K). CCK-8 assay and apoptosis analysis showed that the knockdown of circ_0002078 obviously inhibited cell proliferation (Fig. 7L) and promoted cell apoptosis (p < 0.001; Fig. 7M, N). EHEB cells were blocked in G2/M phase after the transfection of si-circ_0002078 (Fig. 7O-Q). These results suggested that targeting circ_0002078 is a promising direction for the future treatment of CLL.
The interaction mechanism among circ_0002078, miR-185-3p and TCF7L1 is that circ_0002078 can absorb miR-185-3p like a sponge, accordingly suppressing the repression of miR-185-3p on TCF7L1. In CLL, the expression of circ_0002078 is abnormally elevated, therefore the binding of miR-185-3p to TCF7L1 is correspondingly reduced, leading to the decrease of the TCF7L1 mRNA degradation, and eventually resulting in the up-regulation of TCF7L1, which contributes to the promotion of cell proliferation, metastasis and other life processes that can conducive to the development of CLL (Fig. 8).
In the present study, we analyzed ceRNA networks obtained from RNA sequencing in CLL patients and cell lines for the first time. The potential biological functions and prognostic prediction value of the DEGs were explored. A novel risk score model composed of HTN3, IL3RA and NCK1 was established to facilitate prognostic assessment. MiR-324-3p was identified in association with IGHV unmutated status, p53 mutation and del17p. The survival-related ceRNA networks provided fundamental clues for promising prognostic biomarkers and therapeutic vulnerabilities of CLL. Moreover, the validation of circ_0002078/miR-185-3p/TCF7L1 axis provided pioneering insights into the molecular underpinnings of CLL.
CLL is a B-cell malignancy occurs mainly in advanced age and have changeable response to treatments [26, 27]. CLL presents a variable disease course and typically involves chromosomal deletions [28, 29]. The age of onset, genetic profile and other characteristics between different races are heterogenous . Due to the RNA sequencing technology, a plethora of ncRNAs have been explored and confirmed to exert pivotal roles in different physiopathological processes of cells [31, 32]. CeRNA network is a crucial mechanism in post-transcriptional layer of gene translation regulation . These ceRNAs contain many kinds of RNAs, such as circRNAs, lncRNAs and pseudogenes, and they regulate gene expression by competing for the same MREs . Recently, accumulating studies have showed that the lncRNA/circRNA-miRNA-mRNA ceRNA networks played crucial roles in tumor onset and development [35, 36]. However, the roles of ceRNA networks in CLL remain uncovered. In the present investigation, we identified the DEmRNAs, DEmiRNAs, DElncRNAs and DEcircRNAs between CLL patients and healthy donors as well as CLL cells and normal controls. GO analysis results demonstrated that these DEGs between CLL specimens and control were significantly related to sequence-specific DNA binding, chromosomal part and cellular biosynthetic process. KEGG analysis indicated that these DEGs enriched in pathways that were correlated with oncogenesis and tumor development. Notch signaling is deregulated in multiple solid tumors and hematological malignancies . JAK-STAT signaling pathway played important roles in driving aggressive growth, invasion, therapeutic resistance and tumor-mediated immunosuppression . These analyses indicated that the identified DEGs are correlated with the tumorigenesis and progression of CLL, which is worthy of our further exploration.
Understanding of prognostic biomarkers and risk score systems contributes to the prognosis of CLL . In this study, we identified 8 DEGs that associated with the OS of CLL via survival analysis. Among them, IL3RA and NCK1 were independent prognostic indicators for the OS in CLL patients. Subsequently, we constructed a novel risk score model that showed highly robust predictive power for OS. The high-risk group was associated with del(13q14) as well as shorter OS. Distinct immune signatures in two groups indicated that the proportion of macrophages and activated CD4+ memory T cells were larger in high-risk group. Yucai Wang et al. discovered that Richter syndrome nodal tissue had higher infiltration of FOXP3+ T cells and CD163+ macrophages . The activated CLL signature was obviously related to macrophages M2 and activated CD4+ memory T cells . Collectively, these findings implied the intrinsic relations between the risk score model and immune pathways, which provide valuable clues for promising drug discovery in further investigations.
Recently, studies of miRNAs as biomarkers for diagnosis and treatment of disease have attracted great attention. For example, circulating microRNAs are important predictors of fasting blood glucose in prediabetic patients . MicroRNA let-7b downregulates the expression of oncogene AML1-ETO, making it a potential therapeutic target in t (8;21) AML . MiR-142-5p can differentiate between treatment-naïve CML patients who did not respond to imatinib during treatment and those who did . Nevertheless, the roles of miRNA in CLL are largely unexplored. The current study identified DEmiRNAs that related to IGHV mutational status, p53 mutational status and del17p. Intriguingly, the low expression of miR-324-3p was associated with IGHV unmutated status, p53 mutation and del17p, which implied that low expression of miR-324-3p was related to dismal prognosis. The potential functions of miR-324-3p in CLL deserve further discussion.
Based on the influence of DEGs on prognosis and DElncRNA/DEcircRNA in CLL patients, we established the lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA networks. These intricate and intertwined networks contributed to uncovering the regulatory mechanism of the ceRNAs in CLL. The key nodes in the networks presented an obvious relevance to OS, which imply that they play significant roles in survival prediction. Among the DEmiRNAs in our RNA-sequencing data, miR-185-3p was found to be downregulated and associated with chemosensitivity [45, 46] and cell proliferation inhibition . In addition, miR-185-3p is involved in cancer progression as a key molecule in the ceRNA networks [48,49,50]. In CLL, the mechanisms of drug resistance and cell proliferation have not been fully elucidated and the potential role of miR-185-3p in CLL deserves further investigation. Study findings of predicted target genes regulated by miR-185-3p indicated that TCF7L1 is highly expressed in prostate cancer , skin squamous cell carcinoma  and gastric cancer  and is related to their occurrence and progression. The regulation of TCF7L1 by miR-185-3p may be one of its potential mechanisms to promote tumorigenesis and development in CLL. Circ_0002078 was screened by intersecting the predicted miR-185-3p-binding circRNAs with the differentially expressed circRNAs in our RNA-sequencing results.
Therefore, in the current investigation, we verified circ_0002078/miR-185-3p/TCF7L1 axis through qRT-PCR, explored their functions and regulatory relationship. The expression of miR-185-3p was down-regulated in CLL patients and its overexpression inhibited the proliferation of MEC1 cells. The expression and function of circ_0002078 were explored in CLL patients and cell lines. The isoforms of circ_0002078 have not been reported in the existing studies, so we did not detect the isoforms of circ_0002078. The overexpression of circ_0002078 in CLL and its role in promoting proliferation, inhibiting cell apoptosis and increasing cell cycle progression make it a promising target for the treatment of CLL. The dual-luciferase reporter gene assay results indicated that miR-185-3p targeted circ_0002078 and 3’UTR of TCF7L1 at the MREs we predicted. In short, these results indicated the potential role of circ_0002078/miR-185-3p/TCF7L1 aixs in CLL.
In summary, we identified the distinct expression profile of lncRNAs, circRNAs, miRNAs and mRNAs in CLL specimens and CLL cells. Functional enrichment analyses uncovered the biological functions of DEGs. Subsequently, we constructed a risk score model consisting of HTN3, IL3RA and NCK1 to predict OS in CLL. The distinct immune signatures in low- and high-risk groups provided new insight into the exploration of immune molecular underpinnings and immunotherapeutic targets discovery. In addition, we constructed ceRNA networks associated with prognosis. The current study is helpful in enriching deeper understanding of the molecular mechanism in CLL pathogenesis and providing directions for exploring novel prognostic biomarkers.
Availability of data and materials
Chronic lymphocytic leukemia
Competing endogenous RNAs
Differentially expressed miRNAs
Differentially expressed mRNAs
Long noncoding RNAs
MiRNA response elements
Differentially expressed genes
Differentially expressed lncRNAs
Differentially expressed circRNAs
Kyoto Encyclopedia of Genes and Genomes
Fetal bovine serum
Quantitative real-time polymerase chain reaction
- CCK-8 :
Cell Counting Kit-8
Receiver operating characteristic
Area under the curve
Time to first treatment
Kikushige Y. Pathogenesis of chronic lymphocytic leukemia and the development of novel therapeutic strategies. J Clin Exp Hematop. 2020;60(4):146–58.
Rozman C, Montserrat E. Chronic lymphocytic leukemia. N Engl J Med. 1995;333(16):1052–7.
Hallek M. Chronic lymphocytic leukemia: 2020 update on diagnosis, risk stratification and treatment. Am J Hematol. 2019;94(11):1266–87.
Hombach S, Kretz M. Non-coding RNAs: Classification, Biology and Functioning. Adv Exp Med Biol. 2016;937:3–17.
Yue Y, Lin X, Qiu X, Yang L, Wang R. The Molecular Roles and Clinical Implications of Non-Coding RNAs in Gastric Cancer. Front Cell Dev Biol. 2021;9:802745.
Liu Y, Liu X, Lin C, Jia X, Zhu H, Song J, et al. Noncoding RNAs regulate alternative splicing in Cancer. J Exp Clin Cancer Res. 2021;40(1):11.
Agnelli L, Bisognin A, Todoerti K, Manzoni M, Taiana E, Galletti S, et al. Expanding the repertoire of miRNAs and miRNA-offset RNAs expressed in multiple myeloma by small RNA deep sequencing. Blood Cancer J. 2019;9(3):21.
Beermann J, Piccoli MT, Viereck J, Thum T. Non-coding RNAs in Development and Disease: Background, Mechanisms, and Therapeutic Approaches. Physiol Rev. 2016;96(4):1297–325.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146(3):353–8.
Zhong Y, Du Y, Yang X, Mo Y, Fan C, Xiong F, et al. Circular RNAs function as ceRNAs to regulate and control human cancer progression. Mol Cancer. 2018;17(1):79.
Ala U. Competing endogenous RNAs, non-coding RNAs and diseases: an intertwined story. Cells. 2020;9(7):1574.
Hu Y, Zhang Y, Ding M, Xu R. LncRNA TMPO-AS1/miR-126–5p/BRCC3 axis accelerates gastric cancer progression and angiogenesis via activating PI3K/Akt/mTOR pathway. J Gastroenterol Hepatol. 2021;36(7):1877–88.
Wu J, Zheng C, Wang Y, Yang Z, Li C, Fang W, et al. LncRNA APCDD1L-AS1 induces icotinib resistance by inhibition of EGFR autophagic degradation via the miR-1322/miR-1972/miR-324-3p-SIRT5 axis in lung adenocarcinoma. Biomark Res. 2021;9(1):9.
Wang WT, Chen TQ, Zeng ZC, Pan Q, Huang W, Han C, et al. The lncRNA LAMP5-AS1 drives leukemia cell stemness by directly modulating DOT1L methyltransferase activity in MLL leukemia. J Hematol Oncol. 2020;13(1):78.
Bin X, Chen Y, Ma J, Tang R, Zhao Z, Wang K, et al. circ_0001588 Induces the Malignant Progression of Hepatocellular Carcinoma by Modulating miR-874/CDK4 Signaling. J Immunol Res. 2021;2021:3759879.
Hou W, Zhang Y. Circ_0025033 promotes the progression of ovarian cancer by activating the expression of LSM4 via targeting miR-184. Pathol Res Pract. 2021;217:153275.
He J, Chu Z, Lai W, Lan Q, Zeng Y, Lu D, et al. Circular RNA circHERC4 as a novel oncogenic driver to promote tumor metastasis via the miR-556-5p/CTBP2/E-cadherin axis in colorectal cancer. J Hematol Oncol. 2021;14(1):194.
Xia L, Wu L, Bao J, Li Q, Chen X, Xia H, et al. Circular RNA circ-CBFB promotes proliferation and inhibits apoptosis in chronic lymphocytic leukemia through regulating miR-607/FZD3/Wnt/β-catenin pathway. Biochem Biophys Res Commun. 2018;503(1):385–90.
Wu W, Wu Z, Xia Y, Qin S, Li Y, Wu J, et al. Downregulation of circ_0132266 in chronic lymphocytic leukemia promoted cell viability through miR-337-3p/PML axis. Aging (Albany NY). 2019;11(11):3561–73.
Wu Z, Sun H, Liu W, Zhu H, Fu J, Yang C, et al. Circ-RPL15: a plasma circular RNA as novel oncogenic driver to promote progression of chronic lymphocytic leukemia. Leukemia. 2020;34(3):919–23.
Ni J, Hong J, Li Q, Zeng Q, Xia R. Long non-coding RNA CRNDE suppressing cell proliferation is regulated by DNA methylation in chronic lymphocytic leukemia. Leuk Res. 2021;105:106564.
Hallek M, Cheson BD, Catovsky D, Caligaris-Cappio F, Dighiero G, Dohner H, et al. Guidelines for the diagnosis and treatment of chronic lymphocytic leukemia: a report from the International Workshop on Chronic Lymphocytic Leukemia updating the National Cancer Institute-Working Group 1996 guidelines. Blood. 2008;111(12):5446–56.
Zhang Y, Zhou X, Li Y, Xu Y, Lu K, Li P, et al. Inhibition of maternal embryonic leucine zipper kinase with OTSSP167 displays potent anti-leukemic effects in chronic lymphocytic leukemia. Oncogene. 2018;37(41):5520–33.
Lu K, Fang XS, Feng LL, Jiang YJ, Zhou XX, Liu X, et al. The STAT3 inhibitor WP1066 reverses the resistance of chronic lymphocytic leukemia cells to histone deacetylase inhibitors induced by interleukin-6. Cancer Lett. 2015;359(2):250–8.
Zhou X, Fang X, Jiang Y, Geng L, Li X, Li Y, et al. Klotho, an anti-aging gene, acts as a tumor suppressor and inhibitor of IGF-1R signaling in diffuse large B cell lymphoma. J Hematol Oncol. 2017;10(1):37.
Maiti GP, Sinha S, Mahmud H, Boysen J, Mendez MT, Vesely SK, et al. SIRT3 overexpression and epigenetic silencing of catalase regulate ROS accumulation in CLL cells activating AXL signaling axis. Blood Cancer J. 2021;11(5):93.
Öztürk S, Paul Y, Afzal S, Gil-Farina I, Jauch A, Bruch PM, et al. Longitudinal analyses of CLL in mice identify leukemia-related clonal changes including a Myc gain predicting poor outcome in patients. Leukemia. 2022;36(2):64–475.
Dighiero G, Hamblin TJ. Chronic lymphocytic leukaemia. The Lancet. 2008;371(9617):1017–29.
Tari K, Shamsi Z, Reza Ghafari H, Atashi A, Shahjahani M, Abroun S. The role of the genetic abnormalities, epigenetic and microRNA in the prognosis of chronic lymphocytic leukemia. Exp Oncol. 2018;40(4):261–7.
Tian Z, Liu M, Fang X, Zhou X, Li P, Li Y, et al. Distinct Age-Related Clinical Features and Risk Assessment in Chinese With Chronic Lymphocytic Leukemia. Front Oncol. 2022;12:885150.
Cech TR, Steitz JA. The noncoding RNA revolution-trashing old rules to forge new ones. Cell. 2014;157(1):77–94.
George J, Patel T. Noncoding RNA as Therapeutic Targets for Hepatocellular Carcinoma. Semin Liver Dis. 2015;35(01):063–74.
Kartha RV, Subramanian S. Competing endogenous RNAs (ceRNAs): new entrants to the intricacies of gene regulation. Front Genet. 2014;5:8.
Sen R, Ghosal S, Das S, Balti S, Chakrabarti J. Competing endogenous RNA: the key to posttranscriptional regulation. ScientificWorldJournal. 2014;2014:896206.
Qi X, Zhang DH, Wu N, Xiao JH, Wang X, Ma W. ceRNA in cancer: possible functions and clinical implications. J Med Genet. 2015;52(10):710–8.
Karreth FA, Pandolfi PP. ceRNA cross-talk in cancer: when ce-bling rivalries go awry. Cancer Discov. 2013;3(10):1113–21.
Brzozowa M, Mielanczyk L, Michalski M, Malinowski L, Kowalczyk-Ziomek G, Helewski K, et al. Role of Notch signaling pathway in gastric cancer pathogenesis. Contemp Oncol (Pozn). 2013;17(1):1–5.
Ou A, Ott M, Fang D, Heimberger AB. The Role and Therapeutic Targeting of JAK/STAT Signaling in Glioblastoma. Cancers (Basel). 2021;13(3):437.
Yun X, Zhang Y, Wang X. Recent progress of prognostic biomarkers and risk scoring systems in chronic lymphocytic leukemia. Biomark Res. 2020;8:40.
Wang Y, Sinha S, Wellik LE, Secreto CR, Rech KL, Call TG, et al. Distinct immune signatures in chronic lymphocytic leukemia and Richter syndrome. Blood Cancer J. 2021;11(5):86.
Sun C, Chen YC, Martinez AZ, Baptista MJ, Pittaluga S, Liu D, et al. The Immune Microenvironment Shapes Transcriptional and Genetic Heterogeneity in Chronic Lymphocytic Leukemia. Blood Adv. 2022;31:bloodadvances.202106941.
Flowers E, Allen IE, Kanaya AM, Aouizerat BE. Circulating MicroRNAs predict glycemic improvement and response to a behavioral intervention. Biomark Res. 2021;9(1):65.
Johnson DT, Davis AG, Zhou JH, Ball ED, Zhang DE. MicroRNA let-7b downregulates AML1-ETO oncogene expression in t(8;21) AML by targeting its 3’UTR. Exp Hematol Oncol. 2021;10(1):8.
Klümper T, Bruckmueller H, Diewock T, Kaehler M, Haenisch S, Pott C, et al. Expression differences of miR-142-5p between treatment-naïve chronic myeloid leukemia patients responding and non-responding to imatinib therapy suggest a link to oncogenic ABL2, SRI, cKIT and MCL1 signaling pathways critical for development of therapy resistance. Exp Hematol Oncol. 2020;9:26.
Li K, Zhu X, Yuan C. Inhibition of miR-185-3p Confers Erlotinib Resistance Through Upregulation of PFKL/MET in Lung Cancers. Front Cell Dev Biol. 2021;9:677860.
Zhou C, Kong W, Ju T, Xie Q, Zhai L. MiR-185-3p mimic promotes the chemosensitivity of CRC cells via AQP5. Cancer Biol Ther. 2020;21(9):790–8.
Zhang W, Han D. miR-185-3p targets Annexin-A8 to inhibit proliferation in cervical cancer cells. Cytotechnology. 2021;73(4):585–92.
Ou C, Sun Z, He X, Li X, Fan S, Zheng X, et al. Targeting YAP1/LINC00152/FSCN1 Signaling Axis Prevents the Progression of Colorectal Cancer. Adv Sci (Weinh). 2020;7(3):1901380.
Liu X, Meng X, Peng X, Yao Q, Zhu F, Ding Z, et al. Impaired AGO2/miR-185-3p/NRP1 axis promotes colorectal cancer metastasis. Cell Death Dis. 2021;12(4):390.
Li C, Li W, Cao S, Xu J, Qian Y, Pan X, et al. Circ_0058106 promotes proliferation, metastasis and EMT process by regulating Wnt2b/β-catenin/c-Myc pathway through miR-185-3p in hypopharyngeal squamous cell carcinoma. Cell Death Dis. 2021;12(11):1063.
Wen YC, Liu YN, Yeh HL, Chen WH, Jiang KC, Lin SR, et al. TCF7L1 regulates cytokine response and neuroendocrine differentiation of prostate cancer. Oncogenesis. 2021;10(11):81.
Ku AT, Shaver TM, Rao AS, Howard JM, Rodriguez CN, Miao Q, et al. TCF7L1 promotes skin tumorigenesis independently of β-catenin through induction of LCN2. eLife. 2017;6:e23242.
Zhang B, Wu J, Cai Y, Luo M, Wang B, Gu Y. TCF7L1 indicates prognosis and promotes proliferation through activation of Keap1/NRF2 in gastric cancer. Acta Biochim Biophys Sin (Shanghai). 2019;51(4):375–85.
This study was funded by National Natural Science Foundation (No.82000195, No.82070203 and No.81770210); Key Research and Development Program of Shandong Province (No.2018CXGC1213); Translational Research Grant of NCRCH (No.2021WWB02, No.2020ZKMB01); Taishan Scholars Program of Shandong Province; Shandong Provincial Natural Science Foundation (No.ZR2020QH094); Shandong Provincial Engineering Research Center of Lymphoma; Technology Development Project of Jinan City (No.202019182); Academic Promotion Programme of Shandong First Medical University (No.2019QL018, No.2020RC007); Shandong Provincial Hospital Youth Talent Plan.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Functional enrichment analyses of differentially expressed genes (DEGs) in chronic lymphocytic leukemia (CLL) cell lines. (A) GO analysis results showed that changes in MF, CC and BP of DEGs between CLL cell lines and control were mainly enriched in sequence−specific DNA binding, Golgi membrane and cation transmembrane transport. (B) KEGG enrichment analysis of DEGs in CLL cell lines and control were mainly enriched in mTOR signaling pathway, NF-kappaB signaling pathway and cell cycle.
The relationship between immune-related function and prognosis of CLL patients.
The expression of TRIM34, SLC30A10, HOXD4 and circ_0007675 in CLL patients and normal B cells.(A) The expression of TRIM34 in patient specimens were significantly increased. (B) The expression of SLC30A10 in patient specimens were significantly increased. (C) The expression of HOXD4 in patient specimens were significantly increased. (D) The expression of circ_0007675 in patient specimens were significantly increased. All results are expressed as mean ± SEM.
Efficiency of miR-185-3p transfection and functional enrichment analyses of target genes of circ_0002078. (A) Efficiency verification of miR-185-3p overexpression in MEC-1 cells by qRT-PCR. (B) GO analysis results demonstrated that the target genes of circ_0002078 were enriched in positive regulation of gene expression, lymphocyte proliferation, and regulation of B cell apoptotic process. (C) KEGG pathway analysis showed that the target genes were primarily relevant to JAK-STAT signaling pathway and cytokine-cytokine receptor interaction.
The top 10 up-regulation and down-regulation differentially expressed mRNAs (DEmRNAs) between chronic lymphocytic leukemia (CLL) patients and normal B cells. Table S2. The top 10 up-regulation and down-regulation differentially expressed microRNAs (DEmiRNAs) between CLL patients and normal B cells. Table S3. The top 10 up-regulation and down-regulation differentially expressed long noncoding RNAs (DElncRNAs) between CLL patients and normal B cells. Table S4. The top 10 up-regulation and down-regulation differentially expressed circular RNAs (DEcircRNAs) between CLL patients and normal B cells. Table S5. The top 10 up-regulation and down-regulation DEmRNAs between CLL cells and normal B cells. Table S6. The top 10 up-regulation and down-regulation DEmiRNAs between CLL cells and normal B cells. Table S7. The top 10 up-regulation and down-regulation DElncRNAs between CLL cells and normal B cells. Table S8. The top 10 up-regulation and down-regulation DEcircRNAs between CLL cells and normal B cells. Table S9. Primer sequences of the genes tested by qRT-PCR.
About this article
Cite this article
Zhang, X., Han, Y., Hu, X. et al. Competing endogenous RNA networks related to prognosis in chronic lymphocytic leukemia: comprehensive analyses and construction of a novel risk score model. Biomark Res 10, 75 (2022). https://doi.org/10.1186/s40364-022-00423-y