DNA methylome analysis in Burkitt and follicular lymphomas identifies differentially methylated regions linked to somatic mutation and transcriptional control

Although Burkitt lymphomas and follicular lymphomas both have features of germinal center B cells, they are biologically and clinically quite distinct. Here we performed whole-genome bisulfite, genome and transcriptome sequencing in 13 IG-MYC translocation-positive Burkitt lymphoma, nine BCL2 translocation-positive follicular lymphoma and four normal germinal center B cell samples. Comparison of Burkitt and follicular lymphoma samples showed differential methylation of intragenic regions that strongly correlated with expression of associated genes, for example, genes active in germinal center dark-zone and light-zone B cells. Integrative pathway analyses of regions differentially methylated in Burkitt and follicular lymphomas implicated DNA methylation as cooperating with somatic mutation of sphingosine phosphate signaling, as well as the TCF3-ID3 and SWI/SNF complexes, in a large fraction of Burkitt lymphomas. Taken together, our results demonstrate a tight connection between somatic mutation, DNA methylation and transcriptional control in key B cell pathways deregulated differentially in Burkitt lymphoma and other germinal center B cell lymphomas

Systematic Evaluation of Somatic <i>Cis</i>-Regulatory Mutations in Follicular Lymphoma

Background: Follicular lymphoma (FL) is an incurable indolent B cell malignancy characterized in the majority of cases by the t(14;18) translocation. While the mutational landscape of the coding genome is nearing completion, less is known about the characteristics of its noncoding genome. Our expectation is that the distribution of noncoding mutations will be non-random, determined by epigenomics features, such as chromatin modification and accessibility. Our strategy is therefore to integrate whole genome epigenomic marks and mutations in order to enrich for variants with regulatory potential and identify unique mutational processes specific to thesecis-regulatory elements (CREs).Allele-specific expression (ASE) patterns allow further refinement to resolve functional CREs and bona fide mutations associated with changes in gene expression. Methods: The H3K27Ac consensus CREs were determined using the ChIP-Seq data of 9 FL patient primary cells (Koues et al., Immunity2015). DNase I-hypersensitive sites (DHSs) and Hi-C of the B-lymphocyte cell line GM12878 was downloaded from ENCODE. Whole genome sequencing and RNA-Seq of 70 FL patients with relative high tumour cellularity (≥30% for DNA tumour purity and ≥25% for RNA B-cell content) were obtained from the International Cancer Genome Consortium project ICGC MMML-Seq. The mutation rate within CREs and DHSs in comparison to their flanking sequences was estimated using a previous pipeline (Sabarinathan et al., Nature2016). We further developed cis-ASE, an integrated analytic pipeline for the identification of recurrent ASE genes, and significantly associated CREs and mutations (Fig1.A). Results: In total 1.04 million noncoding mutations, corresponding to 14.8K mutations per sample with a median of 9,991 noncoding mutations were identified in our series of 70 FL samples. 62K (6.0%) mutations were located within H3K27Ac bound CREs and there was an elevated mutation rate in H3K27Ac CREs compared to the corresponding flanking regions of 1kb up and downstream (χ² test p = 1.95e-19, Fig1.B). For DHSs, we observed the opposite pattern, with a lower mutation rate in DHSs than in the flanking regions (p = 0.04, Fig1.C), most likely reflecting the higher accessibility to global genome repair machinery in DHSs in relation to flanking sequences. Dividing mutations into high (≥0.3) and low (&amp;lt;0.3) adjusted variant allele frequency (VAF) groups (accounting for tumour purity), we observed significantly higher mutation rate in CREs than in the flanking sequences, that was specific for mutations with higher VAFs (p= 6.95e−38), as the difference was much weaker for low VAF mutations (p= 0.02). Mutation signatures 6 and 20, linked to defective DNA mismatch repair, were highly enriched for mutations within H3K27Ac CREs (44%) versus other regions outside (27%) (p&amp;lt; 0.001). ASE was assessed using 45.6k informative SNPs (1.7%) per sample wherein matched genotype data was available from WGS and RNA-Seq profiles. Cis-ASE identified on average 480 ASE genes per sample (binomial test, adjusted p&amp;lt;0.05), corresponding to 1,943 recurrent ASE genes with a minimum threshold of ≥5 samples. These ASE events were not significantly influenced by local copy number changes or promoter methylation. KEGG pathway analysis of recurrent ASE genes identified adherens junction, B cell receptor signaling pathway and Fc gamma R-mediated phagocytosis as the most overrepresented pathways. Cis-ASE further identified 18 ASE-CRE interactions where CRE mutations were significantly correlated with an imbalance in elevated alternative allele ratios. These ASE genes included recognized lymphoma related genes including BCL2, STAT6, MAF, and additional novel targets. The pattern of aberrant somatic hypermutation (aSHM) was assessed for these 18 CREs, and we narrowed down to 15 significant ASE-CRE interactions not strongly affected by aSHM, consisting of 92 bona fide mutation candidates present in 37 FL samples of our cohort (52.9%). Conclusion: Our study identified unique mutation processes operating in H3K27Ac CREs. Using an integrated genomic approach of whole genome mutations, chromatin marks and RNA-Seq, we explored ASE-CRE interactions, and identified 15 H3K27Ac bound CREs enriched for cis-regulatory mutations significantly associated with ASE and total expression of targeted genes, deserving for further exploration to enrich our understanding of FL noncoding genome. Disclosures Okosun: BeiGene: Membership on an entity's Board of Directors or advisory committees, Research Funding; Gilead Sciences: Research Funding. </jats:sec

Relevance of ID3-TCF3-CCND3 pathway mutations in pediatric aggressive B-cell lymphoma treated according to the non-Hodgkin Lymphoma Berlin-Frankfurt-Münster protocols

Mature B-cell non-Hodgkin lymphoma is the most common subtype of non-Hodgkin lymphoma in childhood and adolescence. B-cell non-Hodgkin lymphomas are further classified into histological subtypes, with Burkitt lymphoma and Diffuse large B-cell lymphoma being the most common subgroups in pediatric patients. Translocations involving the MYC oncogene are known as relevant but not sufficient for Burkitt lymphoma pathogenesis. Recently published large-scale next-generation sequencing studies unveiled sets of additional recurrently mutated genes in samples of pediatric and adult B-cell non-Hodgkin lymphoma patients. ID3, TCF3 and CCND3 are potential drivers of Burkitt lymphomagenesis. In the study herein, frequency and clinical relevance of mutations in ID3, TCF3 and CCND3 were analyzed within a well-defined cohort of 84 uniformly diagnosed and treated pediatric B-cell non-Hodgkin lymphoma patients of the Berlin-Frankfurt-Münster group. Mutation frequency was 78% (ID3), 13% (TCF3) and 36% (CCND3) in Burkitt lymphoma (including Burkitt leukemia). ID3 and CCND3 mutations were associated with more advanced stages of the disease in MYC rearrangement positive Burkitt lymphoma. In conclusion, ID3-TCF3-CCND3 pathway genes are mutated in more than 88% of MYC-rearranged pediatric B-cell non-Hodgkin lymphoma and the pathway may represent a highly relevant second hit of Burkitt lymphoma pathogenesis, especially in children and adolescents

DNA methylome analysis in Burkitt and follicular lymphomas identifies differentially methylated regions linked to somatic mutation and transcriptional control

Although Burkitt lymphomas and follicular lymphomas both have features of germinal center B cells, they are biologically and clinically quite distinct. Here we performed whole-genome bisulfite, genome and transcriptome sequencing in 13 IG-MYC translocation-positive Burkitt lymphoma, nine BCL2 translocation-positive follicular lymphoma and four normal germinal center B cell samples. Comparison of Burkitt and follicular lymphoma samples showed differential methylation of intragenic regions that strongly correlated with expression of associated genes, for example, genes active in germinal center dark-zone and light-zone B cells. Integrative pathway analyses of regions differentially methylated in Burkitt and follicular lymphomas implicated DNA methylation as cooperating with somatic mutation of sphingosine phosphate signaling, as well as the TCF3-ID3 and SWI/SNF complexes, in a large fraction of Burkitt lymphomas. Taken together, our results demonstrate a tight connection between somatic mutation, DNA methylation and transcriptional control in key B cell pathways deregulated differentially in Burkitt lymphoma and other germinal center B cell lymphomas