BCL3 rearrangements in B-cell lymphoid neoplasms occur in two breakpoint clusters associated with different diseases.

The t(14;19)(q32;q13) often juxtaposes BCL3 with immunoglobulin heavy chain (IGH) resulting in overexpression of the gene. In contrast to other oncogenic translocations, BCL3 rearrangement (BCL3-R) has been associated with a broad spectrum of lymphoid neoplasms. Here we report an integrative whole-genome sequence, transcriptomic, and DNA methylation analysis of 13 lymphoid neoplasms with BCL3-R. The resolution of the breakpoints at single base-pair revealed that they occur in two clusters at 5' (n=9) and 3' (n=4) regions of BCL3 associated with two different biological and clinical entities. Both breakpoints were mediated by aberrant class switch recombination of the IGH locus. However, the 5' breakpoints (upstream) juxtaposed BCL3 next to an IGH enhancer leading to overexpression of the gene whereas the 3' breakpoints (downstream) positioned BCL3 outside the influence of the IGH and were not associated with its expression. Upstream BCL3-R tumors had unmutated IGHV, trisomy 12, and mutated genes frequently seen in chronic lymphocytic leukemia (CLL) but had an atypical CLL morphology, immunophenotype, DNA methylome, and expression profile that differ from conventional CLL. In contrast, downstream BCL3-R neoplasms were atypical splenic or nodal marginal zone lymphomas (MZL) with mutated IGHV, complex karyotypes and mutated genes typical of MZL. Two of the latter four tumors transformed to a large B-cell lymphoma. We designed a novel fluorescence in situ hybridization assay that recognizes the two different breakpoints and validated these findings in 17 independent tumors. Overall, upstream or downstream breakpoints of BCL3-R are mainly associated with two subtypes of lymphoid neoplasms with different (epi)genomic, expression, and clinicopathological features resembling atypical CLL and MZL, respectively.

Chronic lymphocytic leukemia patient-derived xenografts recapitulate clonal evolution to Richter transformation.

Chronic lymphocytic leukemia (CLL) is a B-cell neoplasm with a heterogeneous clinical behavior. In 5-10% of patients the disease transforms into a diffuse large-B cell lymphoma known as Richter transformation (RT), which is associated with dismal prognosis. Here, we aimed to establish patient-derived xenograft (PDX) models to study the molecular features and evolution of CLL and RT. We generated two PDXs by injecting CLL (PDX12) and RT (PDX19) cells into immunocompromised NSG mice. Both PDXs were morphologically and phenotypically similar to RT. Whole-genome sequencing analysis at different time points of the PDX evolution revealed a genomic landscape similar to RT tumors from both patients and uncovered an unprecedented RT subclonal heterogeneity and clonal evolution during PDX generation. In PDX12, the transformed cells expanded from a very small subclone already present at the CLL stage. Transcriptomic analysis of PDXs showed a high oxidative phosphorylation (OXPHOS) and low B-cell receptor (BCR) signaling similar to the RT in the patients. IACS-010759, an OXPHOS inhibitor, reduced proliferation, and circumvented resistance to venetoclax. In summary, we have generated new RT-PDX models, one of them from CLL cells that mimicked the evolution of CLL to RT uncovering intrinsic features of RT cells of therapeutical value.

Detection of early seeding of Richter transformation in chronic lymphocytic leukemia.

Richter transformation (RT) is a paradigmatic evolution of chronic lymphocytic leukemia (CLL) into a very aggressive large B cell lymphoma conferring a dismal prognosis. The mechanisms driving RT remain largely unknown. We characterized the whole genome, epigenome and transcriptome, combined with single-cell DNA/RNA-sequencing analyses and functional experiments, of 19 cases of CLL developing RT. Studying 54 longitudinal samples covering up to 19 years of disease course, we uncovered minute subclones carrying genomic, immunogenetic and transcriptomic features of RT cells already at CLL diagnosis, which were dormant for up to 19 years before transformation. We also identified new driver alterations, discovered a new mutational signature (SBS-RT), recognized an oxidative phosphorylation (OXPHOS)high-B cell receptor (BCR)low-signaling transcriptional axis in RT and showed that OXPHOS inhibition reduces the proliferation of RT cells. These findings demonstrate the early seeding of subclones driving advanced stages of cancer evolution and uncover potential therapeutic targets for RT.

Signatures of TOP1 transcription-associated mutagenesis in cancer and germline.

The mutational landscape is shaped by many processes. Genic regions are vulnerable to mutation but are preferentially protected by transcription-coupled repair1. In microorganisms, transcription has been demonstrated to be mutagenic2,3; however, the impact of transcription-associated mutagenesis remains to be established in higher eukaryotes4. Here we show that ID4-a cancer insertion-deletion (indel) mutation signature of unknown aetiology5 characterized by short (2 to 5 base pair) deletions -is due to a transcription-associated mutagenesis process. We demonstrate that defective ribonucleotide excision repair in mammals is associated with the ID4 signature, with mutations occurring at a TNT sequence motif, implicating topoisomerase 1 (TOP1) activity at sites of genome-embedded ribonucleotides as a mechanistic basis. Such TOP1-mediated deletions occur somatically in cancer, and the ID-TOP1 signature is also found in physiological settings, contributing to genic de novo indel mutations in the germline. Thus, although topoisomerases protect against genome instability by relieving topological stress6, their activity may also be an important source of mutations in the human genome.

IGLV3-21R110 identifies an aggressive biological subtype of chronic lymphocytic leukemia with intermediate epigenetics.

B-cell receptor (BCR) signaling is crucial for chronic lymphocytic leukemia (CLL) biology. IGLV3-21-expressing B cells may acquire a single point mutation (R110) that triggers autonomous BCR signaling, conferring aggressive behavior. Epigenetic studies have defined 3 CLL subtypes based on methylation signatures reminiscent of naïve-like (n-CLL), intermediate (i-CLL), and memory-like (m-CLL) B cells with different biological features. i-CLL carries a borderline IGHV mutational load and significantly higher use of IGHV3-21/IGLV3-21. To determine the clinical and biological features of IGLV3-21R110 CLL and its relationship to these epigenetic subtypes, we characterized the immunoglobulin gene of 584 CLL cases using whole-genome/exome and RNA sequencing. IGLV3-21R110 was detected in 6.5% of cases: 30 (38%) of 79 i-CLLs, 5 (1.7%) of 291 m-CLLs, and 1 (0.5%) of 189 n-CLLs. All stereotype subset 2 cases carried IGLV3-21R110, whereas 62% of IGLV3-21R110 i-CLL cases had nonstereotyped BCR immunoglobulins. IGLV3-21R110 i-CLL had a significantly higher number of SF3B1 and ATM mutations and total number of driver alterations. However, the R110 mutation was the sole alteration in 1 i-CLL and was accompanied only by del(13q) in 3. Although IGHV mutational status varied, IGLV3-21R110 i-CLL transcriptomically resembled n-CLL/unmutated IGHV CLL with a specific signature including WNT5A/B overexpression. In contrast, i-CLL lacking IGLV3-21R110 mirrored m-CLL/mutated IGHV. Patients with IGLV3-21R110 i-CLL had a short time to first treatment and overall survival similar to those of n-CLL/unmutated IGHV patients, whereas patients with non-IGLV3-21R110 i-CLL had a good prognosis similar to that of patients with m-CLL/mutated IGHV. IGLV3-21R110 defines a CLL subgroup with specific biological features and an unfavorable prognosis independent of IGHV mutational status and epigenetic subtype.

IgCaller for reconstructing immunoglobulin gene rearrangements and oncogenic translocations from whole-genome sequencing in lymphoid neoplasms.

Immunoglobulin (Ig) gene rearrangements and oncogenic translocations are routinely assessed during the characterization of B cell neoplasms and stratification of patients with distinct clinical and biological features, with the assessment done using Sanger sequencing, targeted next-generation sequencing, or fluorescence in situ hybridization (FISH). Currently, a complete Ig characterization cannot be extracted from whole-genome sequencing (WGS) data due to the inherent complexity of the Ig loci. Here, we introduce IgCaller, an algorithm designed to fully characterize Ig gene rearrangements and oncogenic translocations from short-read WGS data. Using a cohort of 404 patients comprising different subtypes of B cell neoplasms, we demonstrate that IgCaller identifies both heavy and light chain rearrangements to provide additional information on their functionality, somatic mutational status, class switch recombination, and oncogenic Ig translocations. Our data thus support IgCaller to be a reliable alternative to Sanger sequencing and FISH for studying the genetic properties of the Ig loci.

Genomic and epigenomic insights into the origin, pathogenesis, and clinical behavior of mantle cell lymphoma subtypes.

Mantle cell lymphoma (MCL) is a mature B-cell neoplasm initially driven by CCND1 rearrangement with 2 molecular subtypes, conventional MCL (cMCL) and leukemic non-nodal MCL (nnMCL), that differ in their clinicobiological behavior. To identify the genetic and epigenetic alterations determining this diversity, we used whole-genome (n = 61) and exome (n = 21) sequencing (74% cMCL, 26% nnMCL) combined with transcriptome and DNA methylation profiles in the context of 5 MCL reference epigenomes. We identified that open and active chromatin at the major translocation cluster locus might facilitate the t(11;14)(q13;32), which modifies the 3-dimensional structure of the involved regions. This translocation is mainly acquired in precursor B cells mediated by recombination-activating genes in both MCL subtypes, whereas in 8% of cases the translocation occurs in mature B cells mediated by activation-induced cytidine deaminase. We identified novel recurrent MCL drivers, including CDKN1B, SAMHD1, BCOR, SYNE1, HNRNPH1, SMARCB1, and DAZAP1. Complex structural alterations emerge as a relevant early oncogenic mechanism in MCL, targeting key driver genes. Breakage-fusion-bridge cycles and translocations activated oncogenes (BMI1, MIR17HG, TERT, MYC, and MYCN), generating gene amplifications and remodeling regulatory regions. cMCL carried significant higher numbers of structural variants, copy number alterations, and driver changes than nnMCL, with exclusive alterations of ATM in cMCL, whereas TP53 and TERT alterations were slightly enriched in nnMCL. Several drivers had prognostic impact, but only TP53 and MYC aberrations added value independently of genomic complexity. An increasing genomic complexity, together with the presence of breakage-fusion-bridge cycles and high DNA methylation changes related to the proliferative cell history, defines patients with different clinical evolution.

The proliferative history shapes the DNA methylome of B-cell tumors and predicts clinical outcome.

We report a systematic analysis of the DNA methylation variability in 1,595 samples of normal cell subpopulations and 14 tumor subtypes spanning the entire human B-cell lineage. Differential methylation among tumor entities relates to differences in cellular origin and to de novo epigenetic alterations, which allowed us to build an accurate machine learning-based diagnostic algorithm. We identify extensive patient-specific methylation variability in silenced chromatin associated with the proliferative history of normal and neoplastic B cells. Mitotic activity generally leaves both hyper- and hypomethylation imprints, but some B-cell neoplasms preferentially gain or lose DNA methylation. Subsequently, we construct a DNA methylation-based mitotic clock called epiCMIT, whose lapse magnitude represents a strong independent prognostic variable in B-cell tumors and is associated with particular driver genetic alterations. Our findings reveal DNA methylation as a holistic tracer of B-cell tumor developmental history, with implications in the differential diagnosis and prediction of clinical outcome.

A practical guide for mutational signature analysis in hematological malignancies.

Analysis of mutational signatures is becoming routine in cancer genomics, with implications for pathogenesis, classification, prognosis, and even treatment decisions. However, the field lacks a consensus on analysis and result interpretation. Using whole-genome sequencing of multiple myeloma (MM), chronic lymphocytic leukemia (CLL) and acute myeloid leukemia, we compare the performance of public signature analysis tools. We describe caveats and pitfalls of de novo signature extraction and fitting approaches, reporting on common inaccuracies: erroneous signature assignment, identification of localized hyper-mutational processes, overcalling of signatures. We provide reproducible solutions to solve these issues and use orthogonal approaches to validate our results. We show how a comprehensive mutational signature analysis may provide relevant biological insights, reporting evidence of c-AID activity among unmutated CLL cases or the absence of BRCA1/BRCA2-mediated homologous recombination deficiency in a MM cohort. Finally, we propose a general analysis framework to ensure production of accurate and reproducible mutational signature data.

Author Correction: A practical guide for mutational signature analysis in hematological malignancies.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The reference epigenome and regulatory chromatin landscape of chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is a frequent hematological neoplasm in which underlying epigenetic alterations are only partially understood. Here, we analyze the reference epigenome of seven primary CLLs and the regulatory chromatin landscape of 107 primary cases in the context of normal B cell differentiation. We identify that the CLL chromatin landscape is largely influenced by distinct dynamics during normal B cell maturation. Beyond this, we define extensive catalogues of regulatory elements de novo reprogrammed in CLL as a whole and in its major clinico-biological subtypes classified by IGHV somatic hypermutation levels. We uncover that IGHV-unmutated CLLs harbor more active and open chromatin than IGHV-mutated cases. Furthermore, we show that de novo active regions in CLL are enriched for NFAT, FOX and TCF/LEF transcription factor family binding sites. Although most genetic alterations are not associated with consistent epigenetic profiles, CLLs with MYD88 mutations and trisomy 12 show distinct chromatin configurations. Furthermore, we observe that non-coding mutations in IGHV-mutated CLLs are enriched in H3K27ac-associated regulatory elements outside accessible chromatin. Overall, this study provides an integrative portrait of the CLL epigenome, identifies extensive networks of altered regulatory elements and sheds light on the relationship between the genetic and epigenetic architecture of the disease.

Decoding the DNA Methylome of Mantle Cell Lymphoma in the Light of the Entire B Cell Lineage.

We analyzed the in silico purified DNA methylation signatures of 82 mantle cell lymphomas (MCL) in comparison with cell subpopulations spanning the entire B cell lineage. We identified two MCL subgroups, respectively carrying epigenetic imprints of germinal-center-inexperienced and germinal-center-experienced B cells, and we found that DNA methylation profiles during lymphomagenesis are largely influenced by the methylation dynamics in normal B cells. An integrative epigenomic approach revealed 10,504 differentially methylated regions in regulatory elements marked by H3K27ac in MCL primary cases, including a distant enhancer showing de novo looping to the MCL oncogene SOX11. Finally, we observed that the magnitude of DNA methylation changes per case is highly variable and serves as an independent prognostic factor for MCL outcome.

The BLUEPRINT Data Analysis Portal.

The impact of large and complex epigenomic datasets on biological insights or clinical applications is limited by the lack of accessibility by easy, intuitive, and fast tools. Here, we describe an epigenomics comparative cyber-infrastructure (EPICO), an open-access reference set of libraries to develop comparative epigenomic data portals. Using EPICO, large epigenome projects can make available their rich datasets to the community without requiring specific technical skills. As a first instance of EPICO, we implemented the BLUEPRINT Data Analysis Portal (BDAP). BDAP provides a desktop for the comparative analysis of epigenomes of hematopoietic cell types based on results, such as the position of epigenetic features, from basic analysis pipelines. The BDAP interface facilitates interactive exploration of genomic regions, genes, and pathways in the context of differentiation of hematopoietic lineages. This work represents initial steps toward broadly accessible integrative analysis of epigenomic data across international consortia. EPICO can be accessed at https://github.com/inab, and BDAP can be accessed at http://blueprint-data.bsc.es.

Epigenomic analysis detects aberrant super-enhancer DNA methylation in human cancer.

BACKGROUND: One of the hallmarks of cancer is the disruption of gene expression patterns. Many molecular lesions contribute to this phenotype, and the importance of aberrant DNA methylation profiles is increasingly recognized. Much of the research effort in this area has examined proximal promoter regions and epigenetic alterations at other loci are not well characterized. RESULTS: Using whole genome bisulfite sequencing to examine uncharted regions of the epigenome, we identify a type of far-reaching DNA methylation alteration in cancer cells of the distal regulatory sequences described as super-enhancers. Human tumors undergo a shift in super-enhancer DNA methylation profiles that is associated with the transcriptional silencing or the overactivation of the corresponding target genes. Intriguingly, we observe locally active fractions of super-enhancers detectable through hypomethylated regions that suggest spatial variability within the large enhancer clusters. Functionally, the DNA methylomes obtained suggest that transcription factors contribute to this local activity of super-enhancers and that trans-acting factors modulate DNA methylation profiles with impact on transforming processes during carcinogenesis. CONCLUSIONS: We develop an extensive catalogue of human DNA methylomes at base resolution to better understand the regulatory functions of DNA methylation beyond those of proximal promoter gene regions. CpG methylation status in normal cells points to locally active regulatory sites at super-enhancers, which are targeted by specific aberrant DNA methylation events in cancer, with putative effects on the expression of downstream genes.

Non-coding recurrent mutations in chronic lymphocytic leukaemia.

Chronic lymphocytic leukaemia (CLL) is a frequent disease in which the genetic alterations determining the clinicobiological behaviour are not fully understood. Here we describe a comprehensive evaluation of the genomic landscape of 452 CLL cases and 54 patients with monoclonal B-lymphocytosis, a precursor disorder. We extend the number of CLL driver alterations, including changes in ZNF292, ZMYM3, ARID1A and PTPN11. We also identify novel recurrent mutations in non-coding regions, including the 3' region of NOTCH1, which cause aberrant splicing events, increase NOTCH1 activity and result in a more aggressive disease. In addition, mutations in an enhancer located on chromosome 9p13 result in reduced expression of the B-cell-specific transcription factor PAX5. The accumulative number of driver alterations (0 to ≥4) discriminated between patients with differences in clinical behaviour. This study provides an integrated portrait of the CLL genomic landscape, identifies new recurrent driver mutations of the disease, and suggests clinical interventions that may improve the management of this neoplasia.