scTAM-seq enables targeted high-confidence analysis of DNA methylation in single cells.

Single-cell DNA methylation profiling currently suffers from excessive noise and/or limited cellular throughput. We developed scTAM-seq, a targeted bisulfite-free method for profiling up to 650 CpGs in up to 10,000 cells per experiment, with a dropout rate as low as 7%. We demonstrate that scTAM-seq can resolve DNA methylation dynamics across B-cell differentiation in blood and bone marrow, identifying intermediate differentiation states that were previously masked. scTAM-seq additionally queries surface-protein expression, thus enabling integration of single-cell DNA methylation information with cell atlas data. In summary, scTAM-seq is a high-throughput, high-confidence method for analyzing DNA methylation at single-CpG resolution across thousands of single cells.

Mutational mechanisms shaping the coding and noncoding genome of germinal center derived B-cell lymphomas.

B cells have the unique property to somatically alter their immunoglobulin (IG) genes by V(D)J recombination, somatic hypermutation (SHM) and class-switch recombination (CSR). Aberrant targeting of these mechanisms is implicated in lymphomagenesis, but the mutational processes are poorly understood. By performing whole genome and transcriptome sequencing of 181 germinal center derived B-cell lymphomas (gcBCL) we identified distinct mutational signatures linked to SHM and CSR. We show that not only SHM, but presumably also CSR causes off-target mutations in non-IG genes. Kataegis clusters with high mutational density mainly affected early replicating regions and were enriched for SHM- and CSR-mediated off-target mutations. Moreover, they often co-occurred in loci physically interacting in the nucleus, suggesting that mutation hotspots promote increased mutation targeting of spatially co-localized loci (termed hypermutation by proxy). Only around 1% of somatic small variants were in protein coding sequences, but in about half of the driver genes, a contribution of B-cell specific mutational processes to their mutations was found. The B-cell-specific mutational processes contribute to both lymphoma initiation and intratumoral heterogeneity. Overall, we demonstrate that mutational processes involved in the development of gcBCL are more complex than previously appreciated, and that B cell-specific mutational processes contribute via diverse mechanisms to lymphomagenesis.

Dynamics of genome architecture and chromatin function during human B cell differentiation and neoplastic transformation.

To investigate the three-dimensional (3D) genome architecture across normal B cell differentiation and in neoplastic cells from different subtypes of chronic lymphocytic leukemia and mantle cell lymphoma patients, here we integrate in situ Hi-C and nine additional omics layers. Beyond conventional active (A) and inactive (B) compartments, we uncover a highly-dynamic intermediate compartment enriched in poised and polycomb-repressed chromatin. During B cell development, 28% of the compartments change, mostly involving a widespread chromatin activation from naive to germinal center B cells and a reversal to the naive state upon further maturation into memory B cells. B cell neoplasms are characterized by both entity and subtype-specific alterations in 3D genome organization, including large chromatin blocks spanning key disease-specific genes. This study indicates that 3D genome interactions are extensively modulated during normal B cell differentiation and that the genome of B cell neoplasias acquires a tumor-specific 3D genome architecture.

Chromatin activation as a unifying principle underlying pathogenic mechanisms in multiple myeloma.

Multiple myeloma (MM) is a plasma cell neoplasm associated with a broad variety of genetic lesions. In spite of this genetic heterogeneity, MMs share a characteristic malignant phenotype whose underlying molecular basis remains poorly characterized. In the present study, we examined plasma cells from MM using a multi-epigenomics approach and demonstrated that, when compared to normal B cells, malignant plasma cells showed an extensive activation of regulatory elements, in part affecting coregulated adjacent genes. Among target genes up-regulated by this process, we found members of the NOTCH, NF-kB, MTOR signaling, and TP53 signaling pathways. Other activated genes included sets involved in osteoblast differentiation and response to oxidative stress, all of which have been shown to be associated with the MM phenotype and clinical behavior. We functionally characterized MM-specific active distant enhancers controlling the expression of thioredoxin (TXN), a major regulator of cellular redox status and, in addition, identified PRDM5 as a novel essential gene for MM. Collectively, our data indicate that aberrant chromatin activation is a unifying feature underlying the malignant plasma cell phenotype.

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.

Insight into genetic predisposition to chronic lymphocytic leukemia from integrative epigenomics.

Genome-wide association studies have provided evidence for inherited genetic predisposition to chronic lymphocytic leukemia (CLL). To gain insight into the mechanisms underlying CLL risk we analyze chromatin accessibility, active regulatory elements marked by H3K27ac, and DNA methylation at 42 risk loci in up to 486 primary CLLs. We identify that risk loci are significantly enriched for active chromatin in CLL with evidence of being CLL-specific or differentially regulated in normal B-cell development. We then use in situ promoter capture Hi-C, in conjunction with gene expression data to reveal likely target genes of the risk loci. Candidate target genes are enriched for pathways related to B-cell development such as MYC and BCL2 signalling. At 14 loci the analysis highlights 63 variants as the probable functional basis of CLL risk. By integrating genetic and epigenetic information our analysis reveals novel insights into the relationship between inherited predisposition and the regulatory chromatin landscape of CLL.

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.

SOX11, a key oncogenic factor in mantle cell lymphoma.

PURPOSE OF REVIEW: SOX11 has emerged as a key transcription factor in the pathogenesis of mantle cell lymphoma (MCL) whereas it is not expressed in normal B cells or virtually in any other mature B-cell neoplasm. This review will examine the role of SOX11 as a biomarker in MCL, the new information on its transcriptional targets, and the mechanisms regulating its expression in MCL. RECENT FINDINGS: SOX11 is highly expressed in conventional MCL, including cyclin D1-negative cases, but it is not expressed in the indolent leukemic nonnodal MCL subtype. These two MCL subtypes also differ in their cell-of-origin, IGHV mutational status and genomic instability. SOX11 promotes tumor growth of MCL cells in vivo and regulates a broad transcriptional program that includes B-cell differentiation pathways and tumor-microenvironment interactions, among others. The mechanisms upregulating SOX11 in MCL are not well understood but are mediated in part by the three-dimensional reconfiguration of the DNA, bringing together a distant enhancer region and the SOX11 promoter. SUMMARY: SOX11 is a relevant element in the pathogenesis of MCL and has been instrumental to identify two distinct clinicobiological subtypes of this tumor. Further studies should clarify the mechanisms mediating its oncogenic potential and leading to its intriguing expression in these tumors.

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.

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.

Whole-genome fingerprint of the DNA methylome during human B cell differentiation.

We analyzed the DNA methylome of ten subpopulations spanning the entire B cell differentiation program by whole-genome bisulfite sequencing and high-density microarrays. We observed that non-CpG methylation disappeared upon B cell commitment, whereas CpG methylation changed extensively during B cell maturation, showing an accumulative pattern and affecting around 30% of all measured CpG sites. Early differentiation stages mainly displayed enhancer demethylation, which was associated with upregulation of key B cell transcription factors and affected multiple genes involved in B cell biology. Late differentiation stages, in contrast, showed extensive demethylation of heterochromatin and methylation gain at Polycomb-repressed areas, and genes with apparent functional impact in B cells were not affected. This signature, which has previously been linked to aging and cancer, was particularly widespread in mature cells with an extended lifespan. Comparing B cell neoplasms with their normal counterparts, we determined that they frequently acquire methylation changes in regions already undergoing dynamic methylation during normal B cell differentiation.

Whole-epigenome analysis in multiple myeloma reveals DNA hypermethylation of B cell-specific enhancers.

While analyzing the DNA methylome of multiple myeloma (MM), a plasma cell neoplasm, by whole-genome bisulfite sequencing and high-density arrays, we observed a highly heterogeneous pattern globally characterized by regional DNA hypermethylation embedded in extensive hypomethylation. In contrast to the widely reported DNA hypermethylation of promoter-associated CpG islands (CGIs) in cancer, hypermethylated sites in MM, as opposed to normal plasma cells, were located outside CpG islands and were unexpectedly associated with intronic enhancer regions defined in normal B cells and plasma cells. Both RNA-seq and in vitro reporter assays indicated that enhancer hypermethylation is globally associated with down-regulation of its host genes. ChIP-seq and DNase-seq further revealed that DNA hypermethylation in these regions is related to enhancer decommissioning. Hypermethylated enhancer regions overlapped with binding sites of B cell-specific transcription factors (TFs) and the degree of enhancer methylation inversely correlated with expression levels of these TFs in MM. Furthermore, hypermethylated regions in MM were methylated in stem cells and gradually became demethylated during normal B-cell differentiation, suggesting that MM cells either reacquire epigenetic features of undifferentiated cells or maintain an epigenetic signature of a putative myeloma stem cell progenitor. Overall, we have identified DNA hypermethylation of developmentally regulated enhancers as a new type of epigenetic modification associated with the pathogenesis of MM.

Cooperativity of RUNX1 and CSF3R mutations in severe congenital neutropenia: a unique pathway in myeloid leukemogenesis.

Severe congenital neutropenia (CN) is a preleukemic bone marrow failure syndrome with a 20% risk of evolving into leukemia or myelodysplastic syndrome (MDS). Patterns of acquisition of leukemia-associated mutations were investigated using next-generation deep-sequencing in 31 CN patients who developed leukemia or MDS. Twenty (64.5%) of the 31 patients had mutations in RUNX1. A majority of patients with RUNX1 mutations (80.5%) also had acquired CSF3R mutations. In contrast to their high frequency in CN patients who developed leukemia or MDS, RUNX1 mutations were found in only 9 of 307 (2.9%) patients with de novo pediatric acute myeloid leukemia. A sequential analysis at stages prior to overt leukemia revealed RUNX1 mutations to be late events in leukemic transformation. Single-cell analyses in 2 patients showed that RUNX1 and CSF3R mutations were present in the same malignant clone. Functional studies demonstrated elevated granulocyte colony-stimulating factor (G-CSF)-induced proliferation with diminished myeloid differentiation of hematopoietic CD34(+) cells coexpressing mutated forms of RUNX1 and CSF3R. The high frequency of cooperating RUNX1 and CSF3R mutations in CN patients suggests a novel molecular pathway of leukemogenesis: mutations in the hematopoietic cytokine receptor (G-CSFR) in combination with the second mutations in the downstream hematopoietic transcription fator (RUNX1). The detection of both RUNX1 and CSF3R mutations could be used as a marker for identifying CN patients with a high risk of progressing to leukemia or MDS.