- Open Access
HELQ and EGR3 expression correlate with IGHV mutation status and prognosis in chronic lymphocytic leukemia
Journal of Translational Medicine volume 19, Article number: 42 (2021)
IGHV mutation status is a crucial prognostic biomarker for CLL. In the present study, we investigated the transcriptomic signatures associating with IGHV mutation status and CLL prognosis.
The co-expression modules and hub genes correlating with IGHV status, were identified using the GSE28654, by ‘WGCNA’ package and R software (version 4.0.2). The over-representation analysis was performed to reveal enriched cell pathways for genes of correlating modules. Then 9 external cohorts were used to validate the correlation of hub genes expression with IGHV status or clinical features (treatment response, transformation to Richter syndrome, etc.). Moreover, to elucidate the significance of hub genes on disease course and prognosis of CLL patients, the Kaplan–Meier analysis for the OS and TTFT of were performed between subgroups dichotomized by the median expression value of individual hub genes.
2 co-expression modules and 9 hub genes ((FCRL1/FCRL2/HELQ/EGR3LPL/LDOC1/ZNF667/SOWAHC/SEPTIN10) correlating with IGHV status were identified by WGCNA, and validated by external datasets. The modules were found to be enriched in NF-kappaB, HIF-1 and other important pathways, involving cell proliferation and apoptosis. The expression of hub genes was revealed to be significantly different, not only between CLL and normal B cell, but also between various types of lymphoid neoplasms. HELQ expression was found to be related with response of immunochemotherapy treatment significantly (p = 0.0413), while HELQ and ZNF667 were expressed differently between stable CLL and Richter syndrome patients (p < 0.0001 and p = 0.0278, respectively). By survival analysis of subgroups, EGR3 expression was indicated to be significantly associated with TTFT by 2 independent cohorts (GSE39671, p = 0.0311; GSE22762, p = 0.0135). While the expression of HELQ and EGR3 was found to be associated with OS (p = 0.0291 and 0.0114 respectively).The Kras, Hedgehog and IL6-JAK-STAT3 pathways were found to be associating with the expression of hub genes, resulting from GSEA.
The expression of HELQ and EGR3 were correlated with IGHV mutation status in CLL patients. Additionally, the expression of HELQ/EGR3 were prognostic markers for CLL associating with targetable cell signaling pathways.
CLL (chronic lymphocytic leukemia) is characterized by uncontrolled proliferation of monoclonal B cells, and resistance to cell apoptosis. CLL is the most prevalent adult leukemia in Europe and America. The age-adjusted incidence in United States is 4.1 per 100,000 inhabitants, with 4500 estimated deaths . The diagnosis requires monoclonal B cells count more than 5 \(\times \) 103/L in peripheral blood, with characteristic morphology and immunophenotype (typically positive for CD5, CD23, CD19, CD20 and CD200). The disease course of CLL is heterogenous, and the treatment was initiated only in patients with advancing or symptomatic disease. The immunochemotherapy, including anti-CD20 monoclonal antibody and cytotoxic agents (fludarabine and cyclophosphamide, etc.) was the traditional choice. While novel agents, including BTK (Bruton’s Tyrosine Kinase) inhibitor (ibrutinib, zanubrutinib) and BCL2 (B cell lymphoma 2) inhibitor (venetoclax), have greatly improved the survival of CLL patients.
Due to the high heterogeneity of CLL, several risk score systems have bene established, among which 2 most widely accepted systems are Rai and Binet risk stratification developed about 40 years ago . These clinical staging systems seemed insufficient for clinical practice, due to the rapid progress of treatment. More molecular and cytogenetic markers were included in the current scoring system, such as CLL International Prognosis Index (CLL-IPI) , which included TP53 gene mutation/deletion, IGHV mutation status, serum beta2 microglobulin, clinical stage and age. The IGHV (immunoglobulin heavy-chain variable region) gene mutation status is demonstrated to be a pivotal prognostic marker for CLL. The U-CLL (unmutated CLL) patients have shorter lymphocyte doubling time and higher expression of CD38, and correlated more aggressive disease course with shorter TTFT (time to first treatment) in comparison with that of M-CLL (mutated CLL) patients [3,4,5]. Moreover, the unmutated IGHV status predicts unfavorable OS (overall survival) for CLL patients receiving immunochemotherapy . However, the expression signature associating with IGHV mutation has been rarely been investigated.
The transcriptomic analysis based on microarray and RNAseq methods, have preliminarily revealed the expression profiles and altered signaling pathways for CLL, which have provided potential biomarkers and therapeutic targets [7,8,9,10]. By WGCNA (weighted gene co-expression network analysis), the hub genes and co-expression modules associating with IGHV status, were identified and validated. The expression signature of hub genes were also associated with clinical features (response to immunochemotherapy and Richter transformation) and clinical outcomes (OS and TTFT). The flowchart for the overall design of this work was shown in Fig. 1. Our work revealed the transcriptomic signature characterized by co-expression modules, and provided insights and rationales to utilize HELQ/EGR3 expression as prognostic markers for CLL.
The expression matrix and clinical/genetic data was downloaded from GEO database repository (https://www.ncbi.nlm.nih.gov/gds/). The details of GEO datasets used in this study were summarized in Table 1 [7,8,9, 11,12,13,14,15,16]. Since the WGCNA required representative samples and expression information of full-scale genome. Among the 6 independent datasets including IGHV mutation information (Table 1), GSE38611/40570/51529 included only early stage patients (Binet A). And GSE69034 used total lymphocytes instead of purified B cells, which may lead to bias to transcriptomic analysis. GSE9992 was performed in GPL96 platform, which were only able to detect 12,402 genes and resulted in missing value for many genes. Therefore, GSE28654 was selected for WGCNA analysis, which avoided abovementioned problems. The last access to the GEO database is on 2020.10.15.
Weighted gene co-expression network analysis
The weighted gene co-expression network was based on the expression data of the whole genome, by ‘WGCNA’ package  and R software (version 4.0.2). The outliers among samples were detected by hierarchical clustering by average link. The soft threshold power was defined as the minimal beta value which set the scale-free R2 > 0.85. Then inter-gene correlating coefficients were calculated by Pearson’s method, which constructed the matrix of gene adjacency and turned into TOM (topological overlap matrix) sequentially. The minimal size of co-expression modules was limited to 30 genes. Genes in the whole genome were classified into co-expression modules according to TOM-based dissimilarity by average linkage hierarchical clustering method. The first principal component of expression matrix is set as module eigengenes. The module membership of individual genes was defined as the correlation coefficients between gene expression and eigengene of the module. The gene significance of individual genes was defined as the correlation coefficients between gene expression and IGHV status. The modules were identified as targets, with the highest Pearson’s coefficient with IGHV status. Within target modules, the hub genes were defined as gene module membership ≥ 0.8, weighted q value < 0.01, and gene significance ≥ 0.2.
Additionally, PPI (protein–protein interaction) network was established for target modules by STRING database (https://string-db.org/) based on the previous evidence and experiments. The genes of target modules were mapped into STRING, and the criteria of extracting PPI pairs was that confidence ≥ 0.4.
Over-representation analysis for the correlated modules
The ORA (over-representation analysis) was performed for genes in the ‘black’ and ‘purple’ modules using hypergeometric distribution method. The analysis was performed by CPDB online tools  (http://cpdb.molgen.mpg.de/) based on GO (gene ontology) database [19, 20] and KEGG (Kyoto Encyclopedia of Genes and Genomes) database [21, 22] respectively. P value < 0.05 was set as the criteria of enriched pathways.
Validation of hub genes expression signature
To validate the association of hub genes and IGHV mutation status, the expression data was normalized by ‘limma’ package in R software (version 4.0.2), gene expression was compared between M-CLL and U-CLL groups using 6 independent datasets (GSE38611, GSE40570, GSE51529, GSE69034 and GSE9992). Due to the EGR3 expression data was missing in GSE9992, only HELQ/ZNF667/SOWAHC were analyzed for this dataset. Then the comparison of hub genes expression was performed across different lymphoid neoplasms (CLL/FL/DLBCL/MALT-MZL/NMZL/MCL) and normal lymphoid tissues, using GSE50006 and GSE32018 .
The association of hub genes and clinical features of CLL patients
The expression level of hub genes was extracted from the whole transcriptomic datasets. GSE58211  was used to compare the hub genes expression between different Binet stages. GSE10138 dataset was utilized for comparison of hub genes expression between response groups to immunochemotherapy. Moreover, the hub genes transcription was compared between Richter transformed and non-transformed CLL patients of GSE103265.
Survival analysis for hub genes expression in CLL
The GSE39761 provided individual TTFT data for 130 untreated CLL patients. Meanwhile, the GSE22762 consisted of individual TTFT data of 70 CLL patients, and OS data for 107 CLL patients. The cohorts were dichotomized into low and high expression groups by the median value of hub genes expression.
Genome-wide gene expression profile associating with hub genes
Due to the prognostic significance of HELQ and EGR3 genes, the expression correlation analysis was performed to uncover the associating genes. We calculated Pearson’s coefficients for individual genes of the whole genome. R software (version 4.0.2) and ‘stats’ package was utilized for the calculation. The criteria of associating genes were defined as |R2|> 0.45 and p value < 0.05. Then, The GSEA was performed to access the enrichment of HELQ/EGR3 associating genes on signaling pathways, according to MSigDB database [25,26,27] (http://software.broadinstitute.org/gsea/msigdb). The significantly enriched pathways were selected with |NES (normalized enrichment score)|> 1 and q value < 0.05.
The expression data was normalized by ‘normalizeBetweenArrays’ function in ‘limma’ package from R software (version 4.0.2). The unpaired t test was utilized to compare the normalized continuous variables between 2 subgroups. In the situation of comparing between more than 2 groups, the ordinary one-way ANOVA test was used. The logrank test was utilized to test the survival difference between subgroups, using p < 0.05 as cut-off value.
The results of WGCNA
The GSE28654 cohort included 61 M-CLL, 28 U-CLL and 23 CLL with undetermined IGHV status (excluded from WGCNA) . The median age at diagnosis was 61 years old, and the majority of the cohort were in early stages (104 Binet stage A and 8 Binet stage B patients) and untreated (78 untreated vs 34 treated patients) . No outliers were detected by hierarchically clustering with average distance (Additional file 1: Fig. S1). We selected 3 as the soft threshold power according to Fig. 2a. The maximal dissimilarity was set as 15% for merging similar modules, resulting in a total of 14 co-expression modules. The topological overlap degree of individual modules was shown in Fig. 2c, which was generated in the form of topological overlap heatmap after grouping 4 hundred randomly selected genes into modules. Additionally, the eigengene adjacency heatmap demonstrated the relationship between modules (Fig. 2d). The Pearson’s coefficients and p values, generated from the correlation analysis between individual module eigengenes and IGHV status, were shown in Fig. 2e. The correlation between the ‘black’ module and M-CLL was the most significant (R2 = 0.59, p = 9e−12). Whereas the ‘purple’ module was the most correlated with U-CLL (R2 = 0.85, p = 1e−32). 4 genes (FCRL1/FCRL2/HELQ/EGR3) of the ‘black’ module, and 5 genes (LPL/LDOC1/ZNF667/SOWAHC/SEPTIN10) of the ‘purple’ module were identified as hub genes.
To our best knowledge, HELQ/EGR3/ZNF667/SOWAHC were not described to be relevant with CLL previously. As an important regulation factor of DNA repair pathway, HELQ (Helicase POLQ-like protein) expression is demonstrated as an indicator of resistance to platinum based chemotherapy in epithelial ovarian cancer . Hui Cheng et al. demonstrate that EGR3 (Early growth response protein 3) is a strong limiting factor for potential of hematopoietic stem cell proliferation . Recently, the under-expression of EGR3 is demonstrated to be an independent risk factor for metastatic prostate cancer . ZNF667, encoding Zinc finger protein 667, may involve in transcription regulation, the aberrantly hypermethylation of which promoted progression of laryngeal and esophageal squamous cell carcinoma [31, 32]. SOWAHC, encoding Ankyrin repeat domain-containing protein SOWAHC, is demonstrated to be prognostic in bladder cancer  and lung squamous cell carcinoma . The 4 genes were selected as target genes in the following analysis, according to prognostic value in abovementioned studies and correlation with IGHV status in WGCNA. Moreover, the genes of ‘black’ and ‘purple’ modules were mapped into STRING tools, to establish the PPI network (Additional file 2: Fig. S2).
The results of ORA for co-expression modules
For the ‘black’ module, the genes were mainly enriched in biological processes involving with B cells, like B cell activation, B cell proliferation, B cell apoptotic process, negative regulation of B cell receptor signaling pathway, etc. (Fig. 3a). And the molecular functions of the genes were enriched in binding, catalytic activity, etc. whereas the products of genes were predominantly located in cell, membrane, etc. Based on KEGG database, the genes were predominantly enriched in NF-kappa B signaling, JAK-STAT signaling, p53 signaling pathway, etc. (Fig. 3b).
For the ‘purple’ module, the genes were mainly enriched in biological processes, like regulation of activated CD8-positive alpha–beta T cell apoptotic process, signal transduction by p53 class mediator resulting in transcription of p21 class mediator, etc. (Fig. 3c). And the molecular functions of the genes were enriched in transport activity, catalytic activity, etc. whereas the products of genes were predominantly located in cell, membrane, etc. Based on KEGG database, the genes were predominantly enriched in mTOR signaling, Wnt signaling, p53 signaling pathway, etc. (Fig. 3d).
Validation of hub genes by additional independent CLL cohorts
The comparison of hub genes expression level between M-CLL and U-CLL patients based on GSE38211 was shown in Fig. 4, in which the HELQ and EGR3 expression was significantly higher in M-CLL group than that in U-CLL group. While ZNF667 and SOWAHC expression was significantly lower in M-CLL than that in U-CLL group. Similar results were calculated and obtained from other independent datasets (Table 2). The expression of HELQ, and ZNF667 was also significantly different between CLL leukemic cells and normal B cells (Fig. 5, Table 2). The expression was significantly differential across various types of lymphoid neoplasms for HELQ/EGR3/ZNF667 (Fig. 6, Table 2).
Expression of hub genes was associated with treatment response and Richter transformation of CLL patients
No significant association between the Binet clinical stage and expression of all hub genes, was uncovered in our analysis based on GSE58211 (data not shown). The higher expression of HELQ predicted stable disease instead of progressive disease for CLL in GSE10138 (Fig. 7a, p = 0.0413), while other hub genes were not significantly expressed differently (Fig. 7b–d). In Richter transformed CLL patients, the expression of HELQ significantly decreased (p < 0.0001), while the expression of ZNF667 significantly increased (p = 0.0278) in comparison with non-transformed cases (Fig. 8).
HELQ and EGR3 predicted TTFT and OS for CLL patients
The analysis on impact of hub genes expression on TTFT was shown in Fig. 9. EGR3-low group had significantly shorter TTFT in comparison with EGR3-high group based on GSE39671 (Fig. 9b) and GSE22762 (Fig. 9f). Although it’s not significant, the TTFT of HELQ-low group tended to be inferior to HELQ-high group (Fig. 9a, e). No significantly results of TTFT analysis were found for ZNF667 (Fig. 9c, g) and SOWAHC (Fig. 9d, h), either. However, intriguing trends were revealed for ZNF667 (Fig. 9c) and SOWAHC (Fig. 9h), which demonstrated that the survival of subgroups was not significantly different (p > 0.05), but the Kaplan–Meier curves were separated (Fig. 9c, h).
The survival analysis using GSE22762 dataset was shown in Fig. 10. The HELQ-low group showed significantly unfavorable OS in comparison with the HELQ-high group (Fig. 10a). The OS of the EGR3-low group was also inferior over the EGR3-high group (Fig. 10b).
The results of GSEA for HELQ and EGR3
The results of GSEA were shown in Table 3. The activated pathways associating with HELQ expression, included MYC (Myc proto-oncogene protein) targets, E2F (Transcription factor E2F1 targets and DNA repair pathways, etc. while the suppressed pathways associating with HELQ expression, included Hedgehog signaling, Kras signaling pathway etc. (Fig. 11a). The activated pathways correlating with EGR3 expression, included MYC targets, E2F targets, PI3K-Akt-mTOR signaling pathways, etc. while the suppressed pathways correlating with EGR3 expression, included IL6-JAK-STAT3 signaling, Kras signaling pathway etc. (Fig. 11b).
The associated pathways of HELQ and EGR3 were overlapped, among which E2F signaling, MYC signaling and DNA repair pathways were activated in both sets, whereas Kras signaling/inflammatory response were both suppressed.
The expression signature of patients harboring different status of IGHV mutation, has been explored previously [35,36,37,38]. But no researchers investigated the co-expression modules correlating with IGHV status using WGCNA. WGCNA is a clustering method to investigate the scale-free property of gene expression network, in which TOM was used to evaluate the connectivity of individual genes . Then the co-expression modules were identified, genes of which is highly co-expressed and strongly connected in network.
In the results of WGCNA, the ‘black’ module was significantly correlated with M-CLL, and the ‘purple’ module was significantly correlated with U-CLL. Notably, the 2 modules were connected in the protein–protein interaction analysis based on pre-existing experimental evidence in STRING database (Additional file 2: Fig. S2). Moreover, ORA revealed several overlapped pathways enriched by both modules, such as mTOR, Hippo and p53 signaling (Fig. 3). In accordance with our results, mTOR signaling was activated in the majority of CLL patients, but in 2 different manners . In most of U-CLL, mTOR activation occurs downstream of BCR signaling. While a subset of M-CLL is driven by mTOR signaling in a non-BCR-dependent manner.
The ORA based on GO database demonstrated that genes in the ‘black’ module were mainly enriched in biological processes related with B cell activation/proliferation/apoptosis. Of note, negative regulation of B cell receptor signaling were found to be enriched by gene set of the ‘black’ module (p = 0.00185), which involved PLCL2 and FCRL3. PLCL2 expresses in hematopoietic cells, encoding phospholipase C-L2 protein. The B cells in PLCL2-knockout mice were highly proliferative to cross-linking of BCR signaling , suggesting a negative regulation role in BCR signaling. FCRL3, encoding Fc receptor-like protein 3, is reported to be highly expressed in M-CLL instead of U-CLL , which demonstrates the inhibitory potential on BCR signaling . And deletion or downregulation of FCRL3 predicts poor prognosis for CLL patients .
For the ‘black’ module, ORA according to KEGG database, demonstrated the NF-kappaB, HIF-1 and AMPK signaling pathways were enriched for the ‘black’ module. The activity of NF-kappaB signaling is demonstrated to be variable but overall increased in CLL leukemic cells [46,47,48], which is also crucial for survival of leukemic cells  and potentially targetable for CLL . HIF-1a expression and HIF-1 signaling is demonstrated to promote the interaction of CLL leukemic cells and microenvironment , which facilitates the survival and propagation of CLL. The transcription of HIF-1 is increased in TP53-disrupted CLL patients, while the HIF-1 induced interaction, between leukemic cells and stromal cells, is independent of TP53 status . The inhibition of HIF-1 signaling is potentially therapeutic, especially in M-CLL based on our analysis. The previous report indicates AMPK signaling is in the control of apoptosis for CLL cells [53, 54], which can be activated by acadesine in a p53 independent way. Additionally, AMPK signaling can be activated by ATP deleting agents, such as metformin and 8-chloro-adenosine, which will switch the energy-generating pathways, leading to autophagy both in vitro and in vivo [55, 56]. Therefore, the AMPK signaling and autophagy are potential crucial pathways in M-CLL.
In ORA for the ‘purple’ module, Wnt signaling was enriched for U-CLL patients. The Wnt signaling is required for survival of CLL leukemic cells in the functional study , which is activated by somatic mutations  or Notch2 activity from BMSC . Combining with our ORA analysis for the ‘purple’ module, inhibition of Wnt signaling or Notch 2 was potential therapeutic for U-CLL patients.
A total of 9 hub genes were identified by WGCNA, including FCRL1, FCRL2, HELQ, EGR3, LPL, LDOC1, ZNF667, SOWAHC and SEPTIN10. The impact of Fc receptor like molecules (including FCRL1 and FCRL2) have been elucidated in CLL, which predicts the IGHV mutation status and clinical progression , and can act as potential immunotherapeutic targets [61, 62]. LPL, encoding lipoprotein lipase, is expressed in CLL patients with aggressive clinical properties, which promotes activating ligands for PPARα (Peroxisome proliferator-activated receptor alpha) and switch energy source to fatty acid . Additionally, the LPL expression is reported to be a strong prognostic indicator . Hatice Duzkale et al. demonstrate that mRNA expression of LDOC1 is correlated with prognostic markers (cytogenetic markers, IGHV mutation status, and ZAP-70 expression), also a predictor of OS for CLL patients . The expression of SEPTIN10 is reported to an independent prognostic factor for survival  and TTFT  of CLL patients. 4 other genes were not reported in previous studies, and selected as target genes in our study, due to the prognostic value in cancers and correlation with IGHV status.
The correlation of hub genes expression with IGHV status in external cohorts was consistent with the result obtained from WGCNA (Fig. 4, Table 2), in which the HELQ and EGR3 were overexpressed in M-CLL than that of U-CLL significantly, while ZNF667 and SOWAHC were under-expressed.. The differential expression was found to be significant for HELQ and ZNF667 between CLL cells and normal B cells (Fig. 5). Moreover, HELQ, EGR3 and ZNF667 were expressed differentially among diverse types of lymphoid neoplasm (Fig. 6). The above results suggested that the expression signature of hub genes is CLL-specific.
The analysis on GSE10138 indicated higher expression of HELQ correlated with better response to immunochemotherapy (Fig. 7). Additionally, the results obtained from GSE103265 suggested that HELQ and ZNF may serve as potential indicators for Richter transformation (Fig. 8), which may help to predict high risk CLL patients. The association of hub genes expression with clinical features suggested that these genes may involve in pathogenesis of CLL.
As the individual course of early stage CLL is heterogenous, the anticipation for urgency and probability of more aggressive intervention is still not solved. Our results demonstrated that EGR3-high group had significant longer TTFT and more indolent disease course, which were consistent with International Prognostic Score for Early-stage CLL (IPS-E) system, in which U-CLL predicts unfavorable TTFT. Furthermore, our analysis on OS confirmed the survival advantage of HELQ-high as well as EGR3-high group (GSE22762). Although novel target agents have greatly improved the prognosis of CLL patients than traditional immunochemotherapy, the prognostic value of HELQ and EGR3 is still potentially crucial, the re-evaluation of which is worthy to conduct in the novel drug era.
Due to the relevance of HELQ and EGR3 with disease course and prognosis of CLL, we performed GSEA to reveal the significantly activated and/or suppressed signaling pathways correlating with HELQ/EGR3 overexpression (Fig. 11 and Table 3). The GSEA indicated that Kras signaling and Hedgehog signaling was negatively correlated with HELQ expression. The somatic mutations in RAS signaling pathway, including KRAS, occurs in a subset of CLL cases, who more frequently have unmutated IGHV gene and worse TTFT . While the functional activation, instead of activating mutations, of Kras signaling has not been reported yet. The inhibitors of Kras signaling, such as ulixertinib , may provide a new treatment option for HELQ-low or U-CLL patients. The Hedgehog signaling is implicated in the initiation, maintenance , and survival  of CLL cells. CLL patients with activation of Hedgehog pathway are associated with a shorter median TTFT , which may attribute to the inferior clinical outcomes in HELQ-low group. Although vismodegib, a Hedgehog inhibitor, significantly suppressed the Hedgehog signaling in CLL patients, no patients response to vismodegib treatment in a phase II clinical trials . Limited number of included patients and ‘ligand-independent’ activation of Hedgehog pathway may contribute to the failure.
The IL6-JAK-STAT3 signaling was also demonstrated to be negatively correlated with EGR3 expression (Table 3). In CLL cells, extracellular IL6 or BCR signaling induces tyrosine phosphorylation of STAT3 [73, 74], leading to upregulated of anti-apoptosis genes and a survival advantage. The BMSC are reported to interact with CLL cells by modulating JAK2/STAT3 signaling , protecting from CLL cells from ibrutinib. This effect can be reversed by combination of ibrutinib and JAK2 inhibitor (AG490), which triggers apoptosis of CLL cells even in the presence of BMSC. Therefore, additional inhibition of IL6-JAK-STAT pathway may be a potential option for EGR3-low or U-CLL patients, which facilitate the clearance of residual CLL cells in protective bone marrow microenvironment after ibrutinib treatment. The PI3K-Akt-mTOR signaling was positively correlated with EGR3 expression. According to the DNA perturbation based stratification study, signaling mediated via mTOR plays a greater role than canonical BCR signaling for survival/proliferation of a M-CLL subset, as the effect of inhibiting mTOR is greater than BTK . So, for M-CLL patients, high expression of EGR3 may be an indicator for usage of additional mTOR inhibitor.
The MYC target signaling was activated in both HELQ-high and EGR3-high group (Table 3). MYC and downstream targets play a role in antigen induced CLL proliferation [76, 77]. The GSEA indicated the possible pathways related with HELQ/EGR3 expression and provided insights to investigation of personalized therapy.
We identified the co-expression modules and hub genes correlating with IGHV mutation status in CLL patients. The differential expression of hub genes was validated by external cohorts, and associated with clinical features like treatment response and Richter transformation. HELQ and EGR3 expression were prognostic markers and predicted TTFT and OS.
Availability of data and materials
The data that support the findings of this study are available from GEO database (https://www.ncbi.nlm.nih.gov/gds/), which are all publicly available.
Chronic lymphocytic leukemia
Immunoglobulin heavy-chain variable region
Weighed gene co-expression network analysis
Time to first treatment
Gene set enrichment analysis
IGHV gene mutated CLL
IGHV gene unmutated
B cell signaling
Bruton's Tyrosine Kinase
B cell lymphoma 2
Fc receptor-like protein 1
Fc receptor-like protein 2
Early growth response protein 3
Zinc finger protein 667
Ankyrin repeat domain-containing protein SOWAHC
Gene expression omnibus
Topological overlap matrix
Diffuse large B cell lymphoma
Mantle cell lymphomas
Mucosa associating lymphoid tissue
Marginal zone lymphoma
Nodal marginal zone lymphoma
Search Tool for the Retrieval of Interacting Genes/Proteins
Kyoto Encyclopedia of Genes and Genomes
Normalized enrichment score
False discovery rate
Peroxisome proliferator activated receptor
Myc proto-oncogene protein
Hallek M. Chronic lymphocytic leukemia: 2020 update on diagnosis, risk stratification and treatment. Am J Hematol. 2019;94(11):1266–87.
International CLLIPIwg: An international prognostic index for patients with chronic lymphocytic leukaemia (CLL-IPI): a meta-analysis of individual patient data. Lancet Oncol. 2016; 17(6):779–790.
Hamblin TJ, Davis Z, Gardiner A, Oscier DG, Stevenson FK. Unmutated Ig V(H) genes are associated with a more aggressive form of chronic lymphocytic leukemia. Blood. 1999;94(6):1848–54.
Krober A, Seiler T, Benner A, Bullinger L, Bruckle E, Lichter P, Dohner H, Stilgenbauer S. V(H) mutation status, CD38 expression level, genomic aberrations, and survival in chronic lymphocytic leukemia. Blood. 2002;100(4):1410–6.
Hu B, Patel KP, Chen HC, Wang X, Luthra R, Routbort MJ, Kanagal-Shamanna R, Medeiros LJ, Yin CC, Zuo Z, et al. Association of gene mutations with time-to-first treatment in 384 treatment-naive chronic lymphocytic leukaemia patients. Br J Haematol. 2019;187(3):307–18.
Rossi D, Terzi-di-Bergamo L, De Paoli L, Cerri M, Ghilardi G, Chiarenza A, Bulian P, Visco C, Mauro FR, Morabito F, et al. Molecular prediction of durable remission after first-line fludarabine-cyclophosphamide-rituximab in chronic lymphocytic leukemia. Blood. 2015;126(16):1921–4.
Trojani A, Di Camillo B, Tedeschi A, Lodola M, Montesano S, Ricci F, Vismara E, Greco A, Veronese S, Orlacchio A, et al. Gene expression profiling identifies ARSD as a new marker of disease progression and the sphingolipid metabolism as a potential novel metabolism in chronic lymphocytic leukemia. Cancer Biomark. 2011;11(1):15–28.
Fabris S, Mosca L, Cutrona G, Lionetti M, Agnelli L, Ciceri G, Barbieri M, Maura F, Matis S, Colombo M, et al. Chromosome 2p gain in monoclonal B-cell lymphocytosis and in early stage chronic lymphocytic leukemia. Am J Hematol. 2013;88(1):24–31.
Maura F, Cutrona G, Mosca L, Matis S, Lionetti M, Fabris S, Agnelli L, Colombo M, Massucco C, Ferracin M, et al. Association between gene and miRNA expression profiles and stereotyped subset #4 B-cell receptor in chronic lymphocytic leukemia. Leuk Lymphoma. 2015;56(11):3150–8.
Maura F, Mosca L, Fabris S, Cutrona G, Matis S, Lionetti M, Agnelli L, Barbieri M, D’Anca M, Manzoni M, et al. Insulin growth factor 1 receptor expression is associated with NOTCH1 mutation, trisomy 12 and aggressive clinical course in chronic lymphocytic leukaemia. PLoS ONE. 2015;10(3):e0118801.
Morabito F, Mosca L, Cutrona G, Agnelli L, Tuana G, Ferracin M, Zagatti B, Lionetti M, Fabris S, Maura F, et al. Clinical monoclonal B lymphocytosis versus Rai 0 chronic lymphocytic leukemia: A comparison of cellular, cytogenetic, molecular, and clinical features. Clin Cancer Res. 2013;19(21):5890–900.
Fabris S, Mosca L, Todoerti K, Cutrona G, Lionetti M, Intini D, Matis S, Colombo M, Agnelli L, Gentile M, et al. Molecular and transcriptional characterization of 17p loss in B-cell chronic lymphocytic leukemia. Genes Chromosomes Cancer. 2008;47(9):781–93.
Weisser M, Yeh RF, Duchateau-Nguyen G, Palermo G, Nguyen TQ, Shi X, Stinson SY, Yu N, Dufour A, Robak T, et al. PTK2 expression and immunochemotherapy outcome in chronic lymphocytic leukemia. Blood. 2014;124(3):420–5.
Friedman DR, Weinberg JB, Barry WT, Goodman BK, Volkheimer AD, Bond KM, Chen Y, Jiang N, Moore JO, Gockerman JP, et al. A genomic approach to improve prognosis and predict therapeutic response in chronic lymphocytic leukemia. Clin Cancer Res. 2009;15(22):6947–55.
Chuang HY, Rassenti L, Salcedo M, Licon K, Kohlmann A, Haferlach T, Foa R, Ideker T, Kipps TJ. Subnetwork-based analysis of chronic lymphocytic leukemia identifies pathways that associate with disease progression. Blood. 2012;120(13):2639–49.
Herold T, Mulaw MA, Jurinovic V, Seiler T, Metzeler KH, Dufour A, Schneider S, Kakadia PM, Spiekermann K, Mansmann U, et al. High expression of MZB1 predicts adverse prognosis in chronic lymphocytic leukemia, follicular lymphoma and diffuse large B-cell lymphoma and is associated with a unique gene expression signature. Leuk Lymphoma. 2013;54(8):1652–7.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinf. 2008;9:559.
Pentchev K, Ono K, Herwig R, Ideker T, Kamburov A. Evidence mining and novelty assessment of protein-protein interactions with the ConsensusPathDB plugin for Cytoscape. Bioinformatics. 2010;26(21):2796–7.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.
The Gene Ontology C. The Gene Ontology Resource: 20 years and still GOing strong. Nucleic Acids Res. 2019;47(1):D330–D338.
Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 1999;27(1):29–34.
Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016;44(D1):D457-462.
Gomez-Abad C, Pisonero H, Blanco-Aparicio C, Roncador G, Gonzalez-Menchen A, Martinez-Climent JA, Mata E, Rodriguez ME, Munoz-Gonzalez G, Sanchez-Beato M, et al. PIM2 inhibition as a rational therapeutic approach in B-cell lymphoma. Blood. 2011;118(20):5517–27.
Garcia-Barchino MJ, Sarasquete ME, Panizo C, Morscio J, Martinez A, Alcoceba M, Fresquet V, Gonzalez-Farre B, Paiva B, Young KH, et al. Richter transformation driven by Epstein–Barr virus reactivation during therapy-related immunosuppression in chronic lymphocytic leukaemia. J Pathol. 2018;245(1):61–73.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003;34(3):267–73.
Dweep H, Sticht C, Gretz N. In-silico algorithms for the screening of possible microRNA binding sites and their interactions. Curr Genomics. 2013;14(2):127–36.
Long J, Zhu JY, Liu YB, Fu K, Tian Y, Li PY, Yang WQ, Yang SY, Yin JY, Yin G, et al. Helicase POLQ-like (HELQ) as a novel indicator of platinum-based chemoresistance for epithelial ovarian cancer. Gynecol Oncol. 2018;149(2):341–9.
Cheng H, Hao S, Liu Y, Pang Y, Ma S, Dong F, Xu J, Zheng G, Li S, Yuan W, et al. Leukemic marrow infiltration reveals a novel role for Egr3 as a potent inhibitor of normal hematopoietic stem cell proliferation. Blood. 2015;126(11):1302–13.
Shin SH, Kim I, Lee JE, Lee M, Park JW. Loss of EGR3 is an independent risk factor for metastatic progression in prostate cancer. Oncogene. 2020;39(36):5839–54.
Dong Z, Li S, Wu X, Niu Y, Liang X, Yang L, Guo Y, Shen S, Liang J, Guo W. Aberrant hypermethylation-mediated downregulation of antisense lncRNA ZNF667-AS1 and its sense gene ZNF667 correlate with progression and prognosis of esophageal squamous cell carcinoma. Cell Death Dis. 2019;10(12):930.
Meng W, Cui W, Zhao L, Chi W, Cao H, Wang B. Aberrant methylation and downregulation of ZNF667-AS1 and ZNF667 promote the malignant progression of laryngeal squamous cell carcinoma. J Biomed Sci. 2019;26(1):13.
Yang Z, Liu A, Xiong Q, Xue Y, Liu F, Zeng S, Zhang Z, Li Y, Sun Y, Xu C. Prognostic value of differentially methylated gene profiles in bladder cancer. J Cell Physiol. 2019;234(10):18763–72.
Zhu Q, Wang J, Zhang Q, Wang F, Fang L, Song B, Xie C, Liu J. Methylationdriven genes PMPCAP1, SOWAHC and ZNF454 as potential prognostic biomarkers in lung squamous cell carcinoma. Mol Med Rep. 2020;21(3):1285–95.
Tsanousa A, Ntoufa S, Papakonstantinou N, Stamatopoulos K, Angelis L. Study of gene expressions’ correlation structures in subgroups of Chronic Lymphocytic Leukemia Patients. J Biomed Inform. 2019;95:103211.
Prieto D, Seija N, Uriepero A, Souto-Padron T, Oliver C, Irigoin V, Guillermo C, Navarrete MA, Ines Landoni A, Dighiero G, et al. LPL protein in chronic lymphocytic leukaemia have different origins in mutated and unmutated patients. Advances for a new prognostic marker in CLL. Br J Haematol. 2018;182(4):521–5.
Morabito F, Cutrona G, Mosca L, D’Anca M, Matis S, Gentile M, Vigna E, Colombo M, Recchia AG, Bossio S, et al. Surrogate molecular markers for IGHV mutational status in chronic lymphocytic leukemia for predicting time to first treatment. Leuk Res. 2015;39(8):840–5.
Blachly JS, Ruppert AS, Zhao W, Long S, Flynn J, Flinn I, Jones J, Maddocks K, Andritsos L, Ghia EM, et al. Immunoglobulin transcript sequence and somatic hypermutation computation from unselected RNA-seq reads in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A. 2015;112(14):4322–7.
Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabasi AL. Hierarchical organization of modularity in metabolic networks. Science. 2002;297(5586):1551–5.
Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4:17.
Dietrich S, Oles M, Lu J, Sellner L, Anders S, Velten B, Wu B, Hullein J, da Silva LM, Walther T, et al. Drug-perturbation-based stratification of blood cancer. J Clin Invest. 2018;128(1):427–45.
Takenaka K, Fukami K, Otsuki M, Nakamura Y, Kataoka Y, Wada M, Tsuji K, Nishikawa S, Yoshida N, Takenawa T. Role of phospholipase C-L2, a novel phospholipase C-like protein that lacks lipase activity. B-cell receptor signaling. Mol Cell Biol. 2003;23(20):7329–38.
Li FJ, Ding S, Pan J, Shakhmatov MA, Kashentseva E, Wu J, Li Y, Soong SJ, Chiorazzi N, Davis RS. FCRL2 expression predicts IGHV mutation status and clinical progression in chronic lymphocytic leukemia. Blood. 2008;112(1):179–87.
Kochi Y, Myouzen K, Yamada R, Suzuki A, Kurosaki T, Nakamura Y, Yamamoto K. FCRL3, an autoimmune susceptibility gene, has inhibitory potential on B-cell receptor-mediated signaling. J Immunol. 2009;183(9):5502–10.
Daniele G, Abbate A, Turchiano A, Palumbo O, Carella M, Lo Cunsolo C, Iuzzolino P, Lonoce A, Hernandez-Sanchez M, Minoia C, et al. 1q23.1 homozygous deletion and downregulation of Fc receptor-like family genes confer poor prognosis in chronic lymphocytic leukemia. Clin Exp Med. 2019;19(2):261–7.
Furman RR, Asgary Z, Mascarenhas JO, Liou HC, Schattner EJ. Modulation of NF-kappa B activity and apoptosis in chronic lymphocytic leukemia B cells. J Immunol. 2000;164(4):2200–6.
Cuni S, Perez-Aciego P, Perez-Chacon G, Vargas JA, Sanchez A, Martin-Saavedra FM, Ballester S, Garcia-Marco J, Jorda J, Durantez A. A sustained activation of PI3K/NF-kappaB pathway is critical for the survival of chronic lymphocytic leukemia B cells. Leukemia. 2004;18(8):1391–400.
Liu Z, Hazan-Halevy I, Harris DM, Li P, Ferrajoli A, Faderl S, Keating MJ, Estrov Z. STAT-3 activates NF-kappaB in chronic lymphocytic leukemia cells. Mol Cancer Res. 2011;9(4):507–15.
Hewamana S, Alghazal S, Lin TT, Clement M, Jenkins C, Guzman ML, Jordan CT, Neelakantan S, Crooks PA, Burnett AK, et al. The NF-kappaB subunit Rel A is associated with in vitro survival and clinical disease progression in chronic lymphocytic leukemia and represents a promising therapeutic target. Blood. 2008;111(9):4681–9.
Chen Y, Chen L, Yu J, Ghia EM, Choi MY, Zhang L, Zhang S, Sanchez-Lopez E, Widhopf GF 2nd, Messer K, et al. Cirmtuzumab blocks Wnt5a/ROR1 stimulation of NF-kappaB to repress autocrine STAT3 activation in chronic lymphocytic leukemia. Blood. 2019;134(13):1084–94.
Valsecchi R, Coltella N, Belloni D, Ponente M, Ten Hacken E, Scielzo C, Scarfo L, Bertilaccio MT, Brambilla P, Lenti E, et al. HIF-1alpha regulates the interaction of chronic lymphocytic leukemia cells with the tumor microenvironment. Blood. 2016;127(16):1987–97.
Griggio V, Vitale C, Todaro M, Riganti C, Kopecka J, Salvetti C, Bomben R, Bo MD, Magliulo D, Rossi D, et al. HIF-1alpha is over-expressed in leukemic cells from TP53-disrupted patients and is a promising therapeutic target in chronic lymphocytic leukemia. Haematologica. 2020;105(4):1042–54.
Campas C, Lopez JM, Santidrian AF, Barragan M, Bellosillo B, Colomer D, Gil J. Acadesine activates AMPK and induces apoptosis in B-cell chronic lymphocytic leukemia cells but not in T lymphocytes. Blood. 2003;101(9):3674–80.
Van Den Neste E, Van den Berghe G, Bontemps F. AICA-riboside (acadesine), an activator of AMP-activated protein kinase with potential for application in hematologic malignancies. Expert Opin Investig Drugs. 2010;19(4):571–8.
Hardie DG. AMPK and autophagy get connected. EMBO J. 2011;30(4):634–5.
Stellrecht CM, Chen LS, Ayres ML, Dennison JB, Shentu S, Chen Y, Keating MJ, Wierda WG, Gandhi V. Chlorinated adenosine analogue induces AMPK and autophagy in chronic lymphocytic leukaemia cells during therapy. Br J Haematol. 2017;179(2):266–71.
Bologna C, Buonincontri R, Serra S, Vaisitti T, Audrito V, Brusa D, Pagnani A, Coscia M, D’Arena G, Mereu E, et al. SLAMF1 regulation of chemotaxis and autophagy determines CLL patient response. J Clin Invest. 2016;126(1):181–94.
Gandhirajan RK, Staib PA, Minke K, Gehrke I, Plickert G, Schlosser A, Schmitt EK, Hallek M, Kreuzer KA. Small molecule inhibitors of Wnt/beta-catenin/lef-1 signaling induces apoptosis in chronic lymphocytic leukemia cells in vitro and in vivo. Neoplasia. 2010;12(4):326–35.
Wang L, Shalek AK, Lawrence M, Ding R, Gaublomme JT, Pochet N, Stojanov P, Sougnez C, Shukla SA, Stevenson KE, et al. Somatic mutation as a mechanism of Wnt/beta-catenin pathway activation in CLL. Blood. 2014;124(7):1089–98.
Mangolini M, Gotte F, Moore A, Ammon T, Oelsner M, Lutzny-Geier G, Klein-Hitpass L, Williamson JC, Lehner PJ, Durig J, et al. Notch2 controls non-autonomous Wnt-signalling in chronic lymphocytic leukaemia. Nat Commun. 2018;9(1):3839.
Du X, Nagata S, Ise T, Stetler-Stevenson M, Pastan I. FCRL1 on chronic lymphocytic leukemia, hairy cell leukemia, and B-cell non-Hodgkin lymphoma as a target of immunotoxins. Blood. 2008;111(1):338–43.
Kazemi T, Asgarian-Omran H, Hojjat-Farsangi M, Shabani M, Memarian A, Sharifian RA, Razavi SM, Jeddi-Tehrani M, Rabbani H, Shokri F. Fc receptor-like 1–5 molecules are similarly expressed in progressive and indolent clinical subtypes of B-cell chronic lymphocytic leukemia. Int J Cancer. 2008;123(9):2113–9.
Spaner DE, Lee E, Shi Y, Wen F, Li Y, Tung S, McCaw L, Wong K, Gary-Gouy H, Dalloul A, et al. PPAR-alpha is a therapeutic target for chronic lymphocytic leukemia. Leukemia. 2013;27(5):1090–9.
Oppezzo P, Vasconcelos Y, Settegrana C, Jeannel D, Vuillier F, Legarff-Tavernier M, Kimura EY, Bechet S, Dumas G, Brissard M, et al. The LPL/ADAM29 expression ratio is a novel prognosis indicator in chronic lymphocytic leukemia. Blood. 2005;106(2):650–7.
Duzkale H, Schweighofer CD, Coombes KR, Barron LL, Ferrajoli A, O’Brien S, Wierda WG, Pfeifer J, Majewski T, Czerniak BA, et al. LDOC1 mRNA is differentially expressed in chronic lymphocytic leukemia and predicts overall survival in untreated patients. Blood. 2011;117(15):4076–84.
Veer MB, Brooijmans AM, Langerak AW, Verhaaf B, Goudswaard CS, Graveland WJ, Lom K, Valk PJ. The predictive value of lipoprotein lipase for survival in chronic lymphocytic leukemia. Haematologica. 2006;91(1):56–63.
Condoluci A, Terzidi Bergamo L, Langerbeins P, Hoechstetter MA, Herling CD, De Paoli L, Delgado J, Rabe KG, Gentile M, Doubek M, et al. International prognostic score for asymptomatic early-stage chronic lymphocytic leukemia. Blood. 2020;135(21):1859–69.
Gimenez N, Martinez-Trillos A, Montraveta A, Lopez-Guerra M, Rosich L, Nadeu F, Valero JG, Aymerich M, Magnano L, Rozman M, et al. Mutations in the RAS-BRAF-MAPK-ERK pathway define a specific subgroup of patients with adverse clinical features and provide new therapeutic options in chronic lymphocytic leukemia. Haematologica. 2019;104(3):576–86.
Kern D, Regl G, Hofbauer SW, Altenhofer P, Achatz G, Dlugosz A, Schnidar H, Greil R, Hartmann TN, Aberger F. Hedgehog/GLI and PI3K signaling in the initiation and maintenance of chronic lymphocytic leukemia. Oncogene. 2015;34(42):5341–51.
Hegde GV, Peterson KJ, Emanuel K, Mittal AK, Joshi AD, Dickinson JD, Kollessery GJ, Bociek RG, Bierman P, Vose JM, et al. Hedgehog-induced survival of B-cell chronic lymphocytic leukemia cells in a stromal cell microenvironment: a potential new therapeutic target. Mol Cancer Res. 2008;6(12):1928–36.
Ghia EM, Rassenti LZ, Neuberg DS, Blanco A, Yousif F, Smith EN, McPherson JD, Hudson TJ, Harismendy O, et al. Activation of hedgehog signaling associates with early disease progression in chronic lymphocytic leukemia. Blood. 2019;133(25):2651–63.
Houot R, Soussain C, Tilly H, Haioun C, Thieblemont C, Casasnovas O, Bouabdallah K, Morschhauser F, Le Gouill S, Salles GA, et al. Inhibition of Hedgehog signaling for the treatment of lymphoma and CLL: a phase II study from the LYSA. Ann Oncol. 2016;27(7):1349–50.
Rozovski U, Wu JY, Harris DM, Liu Z, Li P, Hazan-Halevy I, Ferrajoli A, Burger JA, O’Brien S, Jain N, et al. Stimulation of the B-cell receptor activates the JAK2/STAT3 signaling pathway in chronic lymphocytic leukemia cells. Blood. 2014;123(24):3797–802.
Rozovski U, Harris DM, Li P, Liu Z, Jain P, Veletic I, Ferrajoli A, Burger J, Thompson P, Jain N, et al. Activation of the B-cell receptor successively activates NF-kappaB and STAT3 in chronic lymphocytic leukemia cells. Int J Cancer. 2017;141(10):2076–81.
Severin F, Frezzato F, Visentin A, Martini V, Trimarco V, Carraro S, Tibaldi E, Brunati AM, Piazza F, Semenzato G, et al. In chronic lymphocytic leukemia the JAK2/STAT3 pathway is constitutively activated and its inhibition leads to CLL cell death unaffected by the protective bone marrow microenvironment. Cancers. 2019;11:12.
Krysov S, Dias S, Paterson A, Mockridge CI, Potter KN, Smith KA, Ashton-Key M, Stevenson FK, Packham G. Surface IgM stimulation induces MEK1/2-dependent MYC expression in chronic lymphocytic leukemia cells. Blood. 2012;119(1):170–9.
Yeomans A, Thirdborough SM, Valle-Argos B, Linley A, Krysov S, Hidalgo MS, Leonard E, Ishfaq M, Wagner SD, Willis AE, et al. Engagement of the B-cell receptor of chronic lymphocytic leukemia cells drives global and MYC-specific mRNA translation. Blood. 2016;127(4):449–57.
This work was supported by the China-Japan Friendship Hospital Fund (2019-1-QN-52).
Ethics approval and consent to participate
This study does not involve animal and/or human tissue/individual data/ participants, there is no ethics related issues. No permissions were required to use any repository data involved in the present study.
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.
The plot of sample clustering to detect outliers for GSE22654.
The protein-protein interaction network for the ‘black’ and ‘purple’ modules. The color of node symbol is black for genes from the ‘black’ module, as well as purple symbol for genes from the ‘purple’ module. The color of node indicated the connectivity degrees (red for high degrees, blue for low degrees).
About this article
Cite this article
Guo, C., Gao, Yy., Ju, Qq. et al. HELQ and EGR3 expression correlate with IGHV mutation status and prognosis in chronic lymphocytic leukemia. J Transl Med 19, 42 (2021). https://doi.org/10.1186/s12967-021-02708-6
- Chronic lymphocytic leukemia
- Survival analysis
- Gene expression profile
- IGHV mutation