Helene Kretzmer

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.

metilene: Fast and sensitive calling of differentially methylated regions from bisulfite sequencing data

The detection of differentially methylated regions (DMRs) is a necessary prerequisite for characterizing different epigenetic states. We present a novel program, metilene, to identify DMRs within whole-genome and targeted data with unrivaled specificity and sensitivity. A binary segmentation algorithm combined with a two-dimensional statistical test allows the detection of DMRs in large methylation experiments with multiple groups of samples in minutes rather than days using off-the-shelf hardware. metilene outperforms other state-of-the-art tools for low coverage data and can estimate missing data. Hence, metilene is a versatile tool to study the effect of epigenetic modifications in differentiation/development, tumorigenesis, and systems biology on a global, genome-wide level. Whether in the framework of international consortia with dozens of samples per group, or even without biological replicates, it produces highly significant and reliable results.

MINCR is a MYC-induced lncRNA able to modulate MYC’s transcriptional network in Burkitt lymphoma cells

Significance Gains of the MYC gene are the most common imbalances in cancer and are associated with poor prognosis, particularly in B-cell lymphoma. Recent advances in DNA sequencing have revealed the existence of thousands of long noncoding RNAs (lncRNAs) with unknown functional relevance. We have here identified a MYC-regulated lncRNA that we named MYC-induced long noncoding RNA (MINCR) that has a strong correlation with MYC expression in cancer. We show that MINCR is functional and controls cell cycle progression by influencing the expression of MYC-regulated cell cycle genes. MINCR is, therefore, a novel player in MYC’s transcriptional network, with the potential to open new therapeutic windows in the fight against malignant lymphoma and, possibly, all cancers that rely on MYC expression. Despite the established role of the transcription factor MYC in cancer, little is known about the impact of a new class of transcriptional regulators, the long noncoding RNAs (lncRNAs), on MYC ability to influence the cellular transcriptome. Here, we have intersected RNA-sequencing data from two MYC-inducible cell lines and a cohort of 91 B-cell lymphomas with or without genetic variants resulting in MYC overexpression. We identified 13 lncRNAs differentially expressed in IG-MYC-positive Burkitt lymphoma and regulated in the same direction by MYC in the model cell lines. Among them, we focused on a lncRNA that we named MYC-induced long noncoding RNA (MINCR), showing a strong correlation with MYC expression in MYC-positive lymphomas. To understand its cellular role, we performed RNAi and found that MINCR knockdown is associated with an impairment in cell cycle progression. Differential gene expression analysis after RNAi showed a significant enrichment of cell cycle genes among the genes down-regulated after MINCR knockdown. Interestingly, these genes are enriched in MYC binding sites in their promoters, suggesting that MINCR acts as a modulator of the MYC transcriptional program. Accordingly, MINCR knockdown was associated with a reduction in MYC binding to the promoters of selected cell cycle genes. Finally, we show that down-regulation of Aurora kinases A and B and chromatin licensing and DNA replication factor 1 may explain the reduction in cellular proliferation observed on MINCR knockdown. We, therefore, suggest that MINCR is a newly identified player in the MYC transcriptional network able to control the expression of cell cycle genes.

Alterations of microRNA and microRNA-regulated messenger RNA expression in germinal center B-cell lymphomas determined by integrative sequencing analysis

MicroRNA are well-established players in post-transcriptional gene regulation. However, information on the effects of microRNA deregulation mainly relies on bioinformatic prediction of potential targets, whereas proof of the direct physical microRNA/target messenger RNA interaction is mostly lacking. Within the International Cancer Genome Consortium Project “Determining Molecular Mechanisms in Malignant Lymphoma by Sequencing”, we performed miRnome sequencing from 16 Burkitt lymphomas, 19 diffuse large B-cell lymphomas, and 21 follicular lymphomas. Twenty-two miRNA separated Burkitt lymphomas from diffuse large B-cell lymphomas/follicular lymphomas, of which 13 have shown regulation by MYC. Moreover, we found expression of three hitherto unreported microRNA. Additionally, we detected recurrent mutations of hsa-miR-142 in diffuse large B-cell lymphomas and follicular lymphomas, and editing of the hsa-miR-376 cluster, providing evidence for microRNA editing in lymphomagenesis. To interrogate the direct physical interactions of microRNA with messenger RNA, we performed Argonaute-2 photoactivatable ribonucleoside-enhanced cross-linking and immunoprecipitation experiments. MicroRNA directly targeted 208 messsenger RNA in the Burkitt lymphomas and 328 messenger RNA in the non-Burkitt lymphoma models. This integrative analysis discovered several regulatory pathways of relevance in lymphomagenesis including Ras, PI3K-Akt and MAPK signaling pathways, also recurrently deregulated in lymphomas by mutations. Our dataset reveals that messenger RNA deregulation through microRNA is a highly relevant mechanism in lymphomagenesis.

BAT: Bisulfite Analysis Toolkit

Here, we present BAT, a modular bisulfite analysis toolkit, that facilitates the analysis of bisulfite sequencing data. It covers the essential analysis steps of read alignment, quality control, extraction of methylation information, and calling of differentially methylated regions, as well as biologically relevant downstream analyses, such as data integration with gene expression, histone modification data, or transcription factor binding site annotation.

Disturbance of poised chromatin coincides with increased expression of developmental genes in cancer

Poised (bivalent or paused) chromatin comprises activating and repressing histone modifications at the same location (Voigt et al., 2013). This combination of epigenetic marks keeps genes expressed at low levels but poised for rapid activation (Margaritis and Holstege, 2008; Mikkelsen et al., 2007). Typically, DNA at poised promoters carries low levels of methylation in normal cells (Meissner et al., 2008; Roadmap Epigenomics Consortium et al., 2015), but frequently shows elevated methylation levels in cancer samples (Hinoue et al., 2012; Gal-Yam et al., 2008; Ohm et al., 2007; Rodriguez et al., 2008; Easwaran et al., 2012). Although higher levels of methylation are normally associated with transcriptional silencing, recently counter-intuitive positive correlations between methylation and expression levels have been reported for two cancer types (Hahn et al., 2014; Kretzmer et al.,). Here, we analyze one of the largest combined expression and methylation data-sets to date, comprising over 5,000 samples and demonstrate that the hypermethylation of poised chromatin in conjunction with up-regulation of the corresponding genes is a general phenomenon in cancer. This up-regulation affects developmental genes and transcription factors, including many genes implicated in cancer. From analysis of 7,000 methylation data sets, we built a universal classifier that can identify cancer samples solely from the hypermethylation status of originally poised chromatin. We reason that the alteration of the epigenetic status of poised chromatin is intimately linked to tumorigenesis.