Multivariate transcriptome analysis identifies networks and key drivers of chronic lymphocytic leukemia relapse risk and patient survival.

BACKGROUND: Chronic lymphocytic leukemia (CLL) is an indolent heme malignancy characterized by the accumulation of CD5+ CD19+ B cells and episodes of relapse. The biological signaling that influence episodes of relapse in CLL are not fully described. Here, we identify gene networks associated with CLL relapse and survival risk. METHODS: Networks were investigated by using a novel weighted gene network co-expression analysis method and examining overrepresentation of upstream regulators and signaling pathways within co-expressed transcriptome modules across clinically annotated transcriptomes from CLL patients (N = 203). Gene Ontology analysis was used to identify biological functions overrepresented in each module. Differential Expression of modules and individual genes was assessed using an ANOVA (Binet Stage A and B relapsed patients) or T-test (SF3B1 mutations). The clinical relevance of biomarker candidates was evaluated using log-rank Kaplan Meier (survival and relapse interval) and ROC tests. RESULTS: Eight distinct modules (M2, M3, M4, M7, M9, M10, M11, M13) were significantly correlated with relapse and differentially expressed between relapsed and non-relapsed Binet Stage A CLL patients. The biological functions of modules positively correlated with relapse were carbohydrate and mRNA metabolism, whereas negatively correlated modules to relapse were protein translation associated. Additionally, M1, M3, M7, and M13 modules negatively correlated with overall survival. CLL biomarkers BTK, BCL2, and TP53 were co-expressed, while unmutated IGHV biomarker ZAP70 and cell survival-associated NOTCH1 were co-expressed in modules positively correlated with relapse and negatively correlated with survival days. CONCLUSIONS: This study provides novel insights into CLL relapse biology and pathways associated with known and novel biomarkers for relapse and overall survival. The modules associated with relapse and overall survival represented both known and novel pathways associated with CLL pathogenesis and can be a resource for the CLL research community. The hub genes of these modules, e.g., ARHGAP27P2, C1S, CASC2, CLEC3B, CRY1, CXCR5, FUT5, MID1IP1, and URAHP, can be studied further as new therapeutic targets or clinical markers to predict CLL patient outcomes.

Transcriptome Network Analysis Identifies CXCL13-CXCR5 Signaling Modules in the Prostate Tumor Immune Microenvironment.

The tumor immune microenvironment (TIME) consists of multiple cell types that contribute to the heterogeneity and complexity of prostate cancer (PCa). In this study, we sought to understand the gene-expression signature of patients with primary prostate tumors by investigating the co-expression profiles of patient samples and their corresponding clinical outcomes, in particular "disease-free months" and "disease reoccurrence". We tested the hypothesis that the CXCL13-CXCR5 axis is co-expressed with factors supporting TIME and PCa progression. Gene expression counts, with clinical attributes from PCa patients, were acquired from TCGA. Profiles of PCa patients were used to identify key drivers that influence or regulate CXCL13-CXCR5 signaling. Weighted gene co-expression network analysis (WGCNA) was applied to identify co-expression patterns among CXCL13-CXCR5, associated genes, and key genetic drivers within the CXCL13-CXCR5 signaling pathway. The processing of downloaded data files began with quality checks using NOISeq, followed by WGCNA. Our results confirmed the quality of the TCGA transcriptome data, identified 12 co-expression networks, and demonstrated that CXCL13, CXCR5 and associated genes are members of signaling networks (modules) associated with G protein coupled receptor (GPCR) responsiveness, invasion/migration, immune checkpoint, and innate immunity. We also identified top canonical pathways and upstream regulators associated with CXCL13-CXCR5 expression and function.