Human germinal center transcriptional programs are de-synchronized in B cell lymphoma.

Most adult B cell lymphomas originate from germinal center (GC) B cells, but it is unclear to what extent B cells in overt lymphoma retain the functional dynamics of GC B cells or are blocked at a particular stage of the GC reaction. Here we used integrative single-cell analysis of phenotype, gene expression and variable-region sequence of the immunoglobulin heavy-chain locus to track the characteristic human GC B cell program in follicular lymphoma B cells. By modeling the cyclic continuum of GC B cell transitional states, we identified characteristic patterns of synchronously expressed gene clusters. GC-specific gene-expression synchrony was lost in single lymphoma B cells. However, distinct follicular lymphoma-specific cell states co-existed within single patient biopsies. Our data show that lymphoma B cells are not blocked in a GC B cell state but might adopt new dynamic modes of functional diversity, which opens the possibility of novel definitions of lymphoma identity.

The alternative RelB NF-κB subunit is a novel critical player in diffuse large B-cell lymphoma.

Diffuse large B-cell lymphoma (DLBCL) is the most frequent lymphoid malignancy affecting adults. The NF-κB transcription factor family is activated by 2 main pathways, the canonical and the alternative NF-κB activation pathway, with different functions. The alternative NF-κB pathway leads to activation of the transcriptionally active RelB NF-κB subunit. Alternative NF-κB activation status and its role in DLBCL pathogenesis remain undefined. Here, we reveal a frequent activation of RelB in a large cohort of DLBCL patients and cell lines, independently of their activated B-cell-like or germinal center B-cell-like subtype. RelB activity defines a new subset of patients with DLBCL and a peculiar gene expression profile and mutational pattern. Importantly, RelB activation does not correlate with the MCD genetic subtype, enriched for activated B-cell-like tumors carrying MYD88L265P and CD79B mutations that cooperatively activate canonical NF-κB, thus indicating that current genetic tools to evaluate NF-κB activity in DLBCL do not provide information on the alternative NF-κB activation. Furthermore, the newly defined RelB-positive subgroup of patients with DLBCL exhibits a dismal outcome after immunochemotherapy. Functional studies revealed that RelB confers DLBCL cell resistance to DNA damage-induced apoptosis in response to doxorubicin, a genotoxic agent used in the front-line treatment of DLBCL. We also show that RelB positivity is associated with high expression of cellular inhibitor of apoptosis protein 2 (cIAP2). Altogether, RelB activation can be used to refine the prognostic stratification of DLBCL and may contribute to subvert the therapeutic DNA damage response in a segment of patients with DLBCL.

Oncogenetic landscape of lymphomagenesis in coeliac disease.

OBJECTIVE: Enteropathy-associated T-cell lymphoma (EATL) is a rare but severe complication of coeliac disease (CeD), often preceded by low-grade clonal intraepithelial lymphoproliferation, referred to as type II refractory CeD (RCDII). Knowledge on underlying oncogenic mechanisms remains scarce. Here, we analysed and compared the mutational landscape of RCDII and EATL in order to identify genetic drivers of CeD-associated lymphomagenesis. DESIGN: Pure populations of RCDII-cells derived from intestinal biopsies (n=9) or sorted from blood (n=2) were analysed by whole exome sequencing, comparative genomic hybridisation and RNA sequencing. Biopsies from RCDII (n=50), EATL (n=19), type I refractory CeD (n=7) and uncomplicated CeD (n=18) were analysed by targeted next-generation sequencing. Moreover, functional in vitro studies and drug testing were performed in RCDII-derived cell lines. RESULTS: 80% of RCDII and 90% of EATL displayed somatic gain-of-functions mutations in the JAK1-STAT3 pathway, including a remarkable p.G1097 hotspot mutation in the JAK1 kinase domain in approximately 50% of cases. Other recurrent somatic events were deleterious mutations in nuclear factor kappa-light-chain-enhancer of activated B-cells (NF-κB) regulators TNFAIP3 and TNIP3 and potentially oncogenic mutations in TET2, KMT2D and DDX3X. JAK1 inhibitors, and the proteasome inhibitor bortezomib could block survival and proliferation of malignant RCDII-cell lines. CONCLUSION: Mutations activating the JAK1-STAT3 pathway appear to be the main drivers of CeD-associated lymphomagenesis. In concert with mutations in negative regulators of NF-κB, they may favour the clonal emergence of malignant lymphocytes in the cytokine-rich coeliac intestine. The identified mutations are attractive therapeutic targets to treat RCDII and block progression towards EATL.

Gene alterations in epigenetic modifiers and JAK-STAT signaling are frequent in breast implant-associated ALCL.

The oncogenic events involved in breast implant-associated anaplastic large cell lymphoma (BI-ALCL) remain elusive. To clarify this point, we have characterized the genomic landscape of 34 BI-ALCLs (15 tumor and 19 in situ subtypes) collected from 54 BI-ALCL patients diagnosed through the French Lymphopath network. Whole-exome sequencing (n = 22, with paired tumor/germline DNA) and/or targeted deep sequencing (n = 24) showed recurrent mutations of epigenetic modifiers in 74% of cases, involving notably KMT2C (26%), KMT2D (9%), CHD2 (15%), and CREBBP (15%). KMT2D and KMT2C mutations correlated with a loss of H3K4 mono- and trimethylation by immunohistochemistry. Twenty cases (59%) showed mutations in ≥1 member of the JAK/STAT pathway, including STAT3 (38%), JAK1 (18%), and STAT5B (3%), and in negative regulators, including SOCS3 (6%), SOCS1 (3%), and PTPN1 (3%). These mutations were more frequent in tumor-type samples than in situ samples (P = .038). All BI-ALCLs expressed pSTAT3, regardless of the mutational status of genes in the JAK/STAT pathway. Mutations in the EOMES gene (12%) involved in lymphocyte development, PI3K-AKT/mTOR (6%), and loss-of-function mutations in TP53 (12%) were also identified. Copy-number aberration (CNA) analysis identified recurrent alterations, including gains on chromosomes 2, 9p, 12p, and 21 and losses on 4q, 8p, 15, 16, and 20. Regions of CNA encompassed genes involved in the JAK/STAT pathway and epigenetic regulators. Our results show that the BI-ALCL genomic landscape is characterized by not only JAK/STAT activating mutations but also loss-of-function alterations of epigenetic modifiers.

Adverse outcome in follicular lymphoma is associated with MYC rearrangements but not MYC extra copies.

Follicular lymphomas (FLs) with MYC rearrangements (MYC-R) and extra copies of MYC (MYC-EC) are rare and the prognosis impact is uncertain. We conducted a retrospective study including 321 FL patients, among whom 259 (81%) had no 8q24 alterations and 62 (19%) were assigned to 8qAlt. Forty-five cases were classified as MYC-EC and six as MYC-R. MYC-R patients were significantly older (P = 0·008), had higher follicular lymphoma international prognostic index (FLIPI) stage (P = 0·05) and ß2-microglobulin (ß2m; P = 0·05). Among patients treated with immuno-chemotherapy, four presented a MYC-R and 25 a MYC-EC. Univariate analysis showed the absence of significant difference between MYC-EC and normal MYC (MYC-NL) regarding progression-free survival (PFS; HR1·3; 95% CI [0·4-1·6]) and specific overall survival (SOS; HR 1·6; 95% CI [0·4-5·7]). Those results were compared to data from the PRIMA trial. This confirmed that MYC-EC had no impact on PFS (P = 0·86) or SOS (P = 0·9). Conversely, MYC-R was associated with a trend to inferior outcome regarding PFS (HR : 6·1; 95% CI [2·2-17·1]; P = 0·00026), lymphoma-related death (SOS; HR 13·6; 95% CI [2·9-65]; P = 0·00014) and risk of transformation (transformation-free survival (TFS); HR 82·7; 95% CI [14·8-463·4]; P < 0·0001). In conclusion, MYC-EC has no prognostic impact in FL but MYC-R FL tended to be associated with an increased risk of transformation and poorer outcome.

A gene-expression profiling score for prediction of outcome in patients with follicular lymphoma: a retrospective training and validation analysis in three international cohorts.

BACKGROUND: Patients with follicular lymphoma have heterogeneous outcomes. Predictor models to distinguish, at diagnosis, between patients at high and low risk of progression are needed. The objective of this study was to use gene-expression profiling data to build and validate a predictive model of outcome for patients treated in the rituximab era. METHODS: A training set of fresh-frozen tumour biopsies was prospectively obtained from 160 untreated patients with high-tumour-burden follicular lymphoma enrolled in the phase 3 randomised PRIMA trial, in which rituximab maintenance was evaluated after rituximab plus chemotherapy induction (median follow-up 6·6 years [IQR 6·0-7·0]). RNA of sufficient quality was obtained for 149 of 160 cases, and Affymetrix U133 Plus 2.0 microarrays were used for gene-expression profiling. We did a multivariate Cox regression analysis to identify genes with expression levels associated with progression-free survival independently of maintenance treatment in a subgroup of 134 randomised patients. Expression levels from 95 curated genes were then determined by digital expression profiling (NanoString technology) in 53 formalin-fixed paraffin-embedded samples of the training set to compare the technical reproducibility of expression levels for each gene between technologies. Genes with high correlation (>0·75) were included in an L2-penalised Cox model adjusted on rituximab maintenance to build a predictive score for progression-free survival. The model was validated using NanoString technology to digitally quantify gene expression in 488 formalin-fixed, paraffin-embedded samples from three independent international patient cohorts from the PRIMA trial (n=178; distinct from the training cohort), the University of Iowa/Mayo Clinic Lymphoma SPORE project (n=201), and the Barcelona Hospital Clinic (n=109). All tissue samples consisted of pretreatment diagnostic biopsies and were confirmed as follicular lymphoma grade 1-3a. The patients were all treated with regimens containing rituximab and chemotherapy, possibly followed by either rituximab maintenance or ibritumomab-tiuxetan consolidation. We determined an optimum threshold on the score to predict patients at low risk and high risk of progression. The model, including the multigene score and the threshold, was initially evaluated in the three validation cohorts separately. The sensitivity and specificity of the score for the prediction of the risk of lymphoma progression at 2 years were assessed on the combined validation cohorts. FINDINGS: In the training cohort, the expression levels of 395 genes were associated with a risk of progression. 23 genes reflecting both B-cell biology and tumour microenvironment with correlation coefficients greater than 0·75 between the two technologies and sample types were retained to build a predictive model that identified a population at an increased risk of progression (p<0·0001). In a multivariate Cox model for progression-free survival adjusted on rituximab maintenance treatment and Follicular Lymphoma International Prognostic Index 1 (FLIPI-1) score, this predictor independently predicted progression (adjusted hazard ratio [aHR] of the high-risk group compared with the low-risk group 3·68, 95% CI 2·19-6·17 [p<0·0001]). The 5-year progression-free survival was 26% (95% CI 16-43) in the high-risk group and 73% (64-83) in the low-risk group. The predictor performances were confirmed in each of the individual validation cohorts (aHR comparing high-risk to low-risk groups 2·57 [95% CI 1·65-4·01] in cohort 1; 2·12 [1·32-3·39] in cohort 2; and 2·11 [1·01-4·41] in cohort 3). In the combined validation cohort, the median progression-free survival was 3·1 years (95% CI 2·4-4·8) in the high-risk group and 10·8 years (10·1-not reached) in the low-risk group (p<0·0001). The risk of lymphoma progression at 2 years was 38% (95% CI 29-46) in the high-risk group and 19% (15-24) in the low-risk group. In a multivariate analysis, the score predicted progression-free survival independently of anti-CD20 maintenance treatment and of the FLIPI score (aHR for the combined cohort 2·30, 95% CI 1·72-3·07). INTERPRETATION: We developed and validated a robust 23-gene expression-based predictor of progression-free survival that is applicable to routinely available formalin-fixed, paraffin-embedded tumour biopsies from patients with follicular lymphoma at time of diagnosis. Applying this score could allow individualised therapy for patients according to their risk category. FUNDING: Roche, SIRIC Lyric, LYSARC, National Institutes of Health, the Henry J Predolin Foundation, and the Spanish Plan Nacional de Investigacion.

BCL2 mutations do not confer adverse prognosis in follicular lymphoma patients treated with rituximab.

BCL2 mutations have been suggested to confer an adverse prognosis to follicular lymphoma (FL) patients, but their prognostic value has not been assessed in patients treated with a rituximab-containing regimen. Here we evaluated the prognostic value of BCL2 mutations in a large prospective cohort of 252 patients with FL treated with immunochemotherapy in the PRIMA randomized trial. Using a DNA-targeted sequencing approach, we detected amino acid altering mutations in 135 patients (54%) and showed that these mutations were probably mediated by the over-activation of AICDA (activation-induced cytidine deaminase) in the context of the t(14;18) translocation. The BCL2 variants identified in PRIMA patients affected the BH1, BH2, and BH3 functional motifs at a lower frequency than the N-terminus and flexible loop domain, with mostly conservative aminoacid changes. With a median follow-up of 6.7 years, we did not observe any impact of BCL2 mutations either on overall survival or progression-free survival.

Whole exome sequencing of relapsed/refractory patients expands the repertoire of somatic mutations in diffuse large B-cell lymphoma.

Despite the many efforts already spent to enumerate somatic mutations in diffuse large B-cell lymphoma (DLBCL), previous whole-genome and whole-exome studies conducted on patients of mixed outcomes failed at characterizing the 30% of patients who will relapse or resist current immunochemotherapies. To address this issue, we performed whole-exome sequencing of normal/tumoral DNA pairs in 14 relapsed/refractory (R/R) patients subclassified by full-transcriptome arrays (six activated B-cell like, three germinal center B-cell like, and five primary mediastinal B-cell lymphomas), from the LNH-03 LYSA clinical trial program. Aside from well-known DLBCL features, gene and pathway level recurrence analyses proposed several interesting leads including TBL1XR1 and activating mutations in IRF4 or in the insulin regulation pathway. Sequencing-based copy number analysis defined 23 short recurrently altered regions involving genes such as REL, CDKN2A, HYAL2, and TP53. Moreover, it highlighted mutations in genes such as GNA13, CARD11, MFHAS1, and PCLO as associated with secondary variant allele amplification events. The five primary mediastinal B-cell lymphomas (PMBL), while unexpected in a R/R cohort, showed a significantly higher mutation rate (P = 0.003) and provided many insights on this classical Hodgkin lymphoma related subtype. Novel genes such as XPO1, MFHAS1, and ITPKB were found particularly mutated, along with various cytokine-based signaling pathways. Among these analyses, somatic events in the NF-κB pathway were found preponderant in the three DLBCL subtypes, confirming its major implication in DLBCL aggressiveness and pinpointing several new candidate genes.

Detection of miRNA regulatory effect on triple negative breast cancer transcriptome.

Identifying key microRNAs (miRNAs) contributing to the genesis and development of a particular disease is a focus of many recent studies. We introduce here a rank-based algorithm to detect miRNA regulatory activity in cancer-derived tissue samples which combines measurements of gene and miRNA expression levels and sequence-based target predictions. The method is designed to detect modest but coordinated changes in the expression of sequence-based predicted target genes. We applied our algorithm to a cohort of 129 tumour and healthy breast tissues and showed its effectiveness in identifying functional miRNAs possibly involved in the disease. These observations have been validated using an independent publicly available breast cancer dataset from The Cancer Genome Atlas. We focused on the triple negative breast cancer subtype to highlight potentially relevant miRNAs in this tumour subtype. For those miRNAs identified as potential regulators, we characterize the function of affected target genes by enrichment analysis. In the two independent datasets, the affected targets are not necessarily the same, but display similar enriched categories, including breast cancer related processes like cell substrate adherens junction, regulation of cell migration, nuclear pore complex and integrin pathway. The R script implementing our method together with the datasets used in the study can be downloaded here (http://bioinfo-out.curie.fr/projects/targetrunningsum).

Comprehensive genetic profiling reveals frequent alterations of driver genes on the X chromosome in extranodal NK/T-cell lymphoma.

Extranodal NK/T-cell lymphoma (ENKTCL) is an Epstein-Barr virus (EBV)-related neoplasm with male dominance and a poor prognosis. A better understanding of the genetic alterations and their functional roles in ENKTCL could help improve patient stratification and treatments. Here, we performed comprehensive genetic analysis of 177 ENKTCL cases to delineate the landscape of mutations, copy number alterations (CNAs), and structural variations, identifying 34 driver genes including six previously unappreciated ones, namely HLA-B, HLA-C, ROBO1, CD58, POT1, and MAP2K1. Among them, CD274 (24%) was the most frequently altered, followed by TP53 (20%), CDKN2A (19%), ARID1A (15%), HLA-A (15%), BCOR (14%), and MSN (14%). Chromosome X (chrX) losses were the most common arm-level CNAs in females (~40%), and alterations of four X-linked driver genes (MSN, BCOR, DDX3X, and KDM6A) were more frequent in males and females harboring chrX losses. Among X-linked drivers, MSN was the most recurrently altered, and its expression was lost in approximately one-third of cases using immunohistochemical analysis. Functional studies of human cell lines demonstrated that MSN disruption promoted cell proliferation and NF-κB activation. Moreover, MSN inactivation increased sensitivity to NF-κB inhibition in vitro and in vivo. In addition, recurrent deletions were observed at the origin of replication in the EBV genome (6%). Finally, by integrating the 34 drivers and 19 significant arm-level CNAs, non-negative matrix factorization and consensus clustering identified two molecular groups with different genetic features and prognosis irrespective of clinical prognostic factors. Together, these findings could help improve diagnostic and therapeutic strategies in ENKTCL.

Critical reasoning on causal inference in genome-wide linkage and association studies.

Genome-wide linkage and association studies of tens of thousands of clinical and molecular traits are currently underway, offering rich data for inferring causality between traits and genetic variation. However, the inference process is based on discovering subtle patterns in the correlation between traits and is therefore challenging and could create a flood of untrustworthy causal inferences. Here we introduce the concerns and show that they are already valid in simple scenarios of two traits linked to or associated with the same genomic region. We argue that more comprehensive analysis and Bayesian reasoning are needed and that these can overcome some of the pitfalls, although not in every conceivable case. We conclude that causal inference methods can still be of use in the iterative process of mathematical modeling and biological validation.

Expression quantitative trait loci are highly sensitive to cellular differentiation state.

Genetical genomics is a strategy for mapping gene expression variation to expression quantitative trait loci (eQTLs). We performed a genetical genomics experiment in four functionally distinct but developmentally closely related hematopoietic cell populations isolated from the BXD panel of recombinant inbred mouse strains. This analysis allowed us to analyze eQTL robustness/sensitivity across different cellular differentiation states. Although we identified a large number (365) of "static" eQTLs that were consistently active in all four cell types, we found a much larger number (1,283) of "dynamic" eQTLs showing cell-type-dependence. Of these, 140, 45, 531, and 295 were preferentially active in stem, progenitor, erythroid, and myeloid cells, respectively. A detailed investigation of those dynamic eQTLs showed that in many cases the eQTL specificity was associated with expression changes in the target gene. We found no evidence for target genes that were regulated by distinct eQTLs in different cell types, suggesting that large-scale changes within functional regulatory networks are uncommon. Our results demonstrate that heritable differences in gene expression are highly sensitive to the developmental stage of the cell population under study. Therefore, future genetical genomics studies should aim at studying multiple well-defined and highly purified cell types in order to construct as comprehensive a picture of the changing functional regulatory relationships as possible.

Undifferentiated embryonic cell transcription factor 1 regulates ESC chromatin organization and gene expression.

Previous reports showed that embryonic stem (ES) cells contain hyperdynamic and globally transcribed chromatin-properties that are important for ES cell pluripotency and differentiation. Here, we demonstrate a role for undifferentiated embryonic cell transcription factor 1 (UTF1) in regulating ES cell chromatin structure. Using chromatin immunoprecipitation-on-chip analysis, we identified >1,700 UTF1 target genes that significantly overlap with previously identified Nanog, Oct4, Klf-4, c-Myc, and Rex1 targets. Gene expression profiling showed that UTF1 knock down results in increased expression of a large set of genes, including a significant number of UTF1 targets. UTF1 knock down (KD) ES cells are, irrespective of the increased expression of several self-renewal genes, Leukemia inhibitory factor (LIF) dependent. However, UTF1 KD ES cells are perturbed in their differentiation in response to dimethyl sulfoxide (DMSO) or after LIF withdrawal and display increased colony formation. UTF1 KD ES cells display extensive chromatin decondensation, reflected by a dramatic increase in nucleosome release on micrococcal nuclease (MNase) treatment and enhanced MNase sensitivity of UTF1 target genes in UTF1 KD ES cells. Summarizing, our data show that UTF1 is a key chromatin component in ES cells, preventing ES cell chromatin decondensation, and aberrant gene expression; both essential for proper initiation of lineage-specific differentiation of ES cells.

High-grade Follicular Lymphomas Exhibit Clinicopathologic, Cytogenetic, and Molecular Diversity Extending Beyond Grades 3A and 3B.

Although follicular lymphoma (FL) is usually graded as FL1-2, FL3A, and FL3B, some borderline cases can be observed and led us to investigate the clinicopathologic diversity of grade 3 FL (FL3). Among 2449 FL patients enrolled in Lymphoma Study Association (LYSA) trials, 1921 cases with sufficient material underwent a central pathologic review. The resulting diagnoses comprised 89.6% FL1-2 (n=1723), 7.2% FL3A (n=138), and 0.5% purely follicular FL3B (n=9). The remaining 51 unclassifiable cases (2.7%) exhibited high-grade features but did not meet WHO criteria for either FL3A or FL3B; and were considered as "unconventional" high-grade FL (FL3U). FL3U morphological pattern consisted of nodular proliferation of large cleaved cells or small-sized to medium-sized blast cells. Compared with FL3A, FL3U exhibited higher MUM1 and Ki67 expression, less BCL2 breaks and more BCL6 rearrangements, together with a higher number of cases without any BCL2, BCL6 or MYC rearrangement. FL3U harbored less frequent mutations in BCL2, KMT2D, KMT2B, and CREBBP than FL3A. MYC and BCL2 were less frequently mutated in FL3U than FL3B. Rituximab cyclophosphamide, doxorubicin, vincristine, and prednisone treated FL3U patients had a worse survival than FL1-2 patients with similar follicular lymphoma international prognostic index and treatment. These results suggest that high-grade FLs encompass a heterogeneous spectrum of tumors with variable morphology and genomic alterations, including FL3U cases that do not strictly fit WHO criteria for either FL3A or FL3B, and display a worse outcome than FL1-2. The distinction of FL3U may be useful to allow a better comprehension of high-grade FLs and to design clinical trials.

Lack of intrafollicular memory CD4 + T cells is predictive of early clinical failure in newly diagnosed follicular lymphoma.

Despite a characteristic indolent course, a substantial subset of follicular lymphoma (FL) patients has an early relapse with a poor outcome. Cells in the microenvironment may be a key contributor to treatment failure. We used a discovery and validation study design to identify microenvironmental determinants of early failure and then integrated these results into the FLIPI. In total, 496 newly diagnosed FL grade 1-3 A patients who were prospectively enrolled into the MER cohort from 2002 to 2012 were evaluated. Tissue microarrays were stained for CD4, CD8, FOXP3, CD32b, CD14, CD68, CD70, SIRP-α, TIM3, PD-1, and PD-L1. Early failure was defined as failing to achieve event-free survival at 24 months (EFS24) in immunochemotherapy-treated patients and EFS12 in all others. CyTOF and CODEX analysis were performed to characterize intratumoral immunophenotypes. Lack of intrafollicular CD4 expression was the only predictor of early failure that replicated with a pooled OR 2.37 (95%CI 1.48-3.79). We next developed a bio-clinical risk model (BioFLIPI), where lack of CD4 intrafollicular expression moved patients up one FLIPI risk group, adding a new fourth high-risk group. Compared with BioFLIPI score of 1, patients with a score of 2 (OR 2.17; 95% CI 1.08-4.69), 3 (OR 3.53; 95% CI 1.78-7.54), and 4 (OR 8.92; 95% CI 4.00-21.1) had increasing risk of early failure. The favorable intrafollicular CD4 T cells were identified as activated central memory T cells, whose prognostic value was independent from genetic features. In conclusion, lack of intrafollicular CD4 expression predicts early failure in FL and combined with FLIPI improves identification of high-risk patients; however, independent validation is warranted.

DiffCoEx: a simple and sensitive method to find differentially coexpressed gene modules.

BACKGROUND: Large microarray datasets have enabled gene regulation to be studied through coexpression analysis. While numerous methods have been developed for identifying differentially expressed genes between two conditions, the field of differential coexpression analysis is still relatively new. More specifically, there is so far no sensitive and untargeted method to identify gene modules (also known as gene sets or clusters) that are differentially coexpressed between two conditions. Here, sensitive and untargeted means that the method should be able to construct de novo modules by grouping genes based on shared, but subtle, differential correlation patterns. RESULTS: We present DiffCoEx, a novel method for identifying correlation pattern changes, which builds on the commonly used Weighted Gene Coexpression Network Analysis (WGCNA) framework for coexpression analysis. We demonstrate its usefulness by identifying biologically relevant, differentially coexpressed modules in a rat cancer dataset. CONCLUSIONS: DiffCoEx is a simple and sensitive method to identify gene coexpression differences between multiple conditions.

Refining diffuse large B-cell lymphoma subgroups using integrated analysis of molecular profiles.

BACKGROUND: Gene expression profiling (GEP), next-generation sequencing (NGS) and copy number variation (CNV) analysis have led to an increasingly detailed characterization of the genomic profiles of DLBCL. The aim of this study was to perform a fully integrated analysis of mutational, genomic, and expression profiles to refine DLBCL subtypes. A comparison of our model with two recently published integrative DLBCL classifiers was carried out, in order to best reflect the current state of genomic subtypes. METHODS: 223 patients with de novo DLBCL from the prospective, multicenter and randomized LNH-03B LYSA clinical trials were included. GEP data was obtained using Affymetrix GeneChip arrays, mutational profiles were established by Lymphopanel NGS targeting 34 key genes, CNV analysis was obtained by array CGH, and FISH and IHC were performed. Unsupervised independent component analysis (ICA) was applied to GEP data and integrated analysis of multi-level molecular data associated with each component (gene signature) was performed. FINDINGS: ICA identified 38 components reflecting transcriptomic variability across our DLBCL cohort. Many of the components were closely related to well-known DLBCL features such as cell-of-origin, stromal and MYC signatures. A component linked to gain of 19q13 locus, among other genomic alterations, was significantly correlated with poor OS and PFS. Through this integrated analysis, a high degree of heterogeneity was highlighted among previously described DLBCL subtypes. INTERPRETATION: The results of this integrated analysis enable a global and multi-level view of DLBCL, as well as improve our understanding of DLBCL subgroups.