Competing endogenous RNA networks related to prognosis in chronic lymphocytic leukemia: comprehensive analyses and construction of a novel risk score model

Abstract

Abstract Background 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. Methods 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. Results 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. Conclusions 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.