Anke K. Bergmann

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.

A B-cell epigenetic signature defines three biologic subgroups of chronic lymphocytic leukemia with clinical impact.

Prospective identification of patients with chronic lymphocytic leukemia (CLL) destined to progress would greatly facilitate their clinical management. Recently, whole-genome DNA methylation analyses identified three clinicobiologic CLL subgroups with an epigenetic signature related to different normal B-cell counterparts. Here, we developed a clinically applicable method to identify these subgroups and to study their clinical relevance. Using a support vector machine approach, we built a prediction model using five epigenetic biomarkers that was able to classify CLL patients accurately into the three subgroups, namely naive B-cell-like, intermediate and memory B-cell-like CLL. DNA methylation was quantified by highly reproducible bisulfite pyrosequencing assays in two independent CLL series. In the initial series (n=211), the three subgroups showed differential levels of IGHV (immunoglobulin heavy-chain locus) mutation (P<0.001) and VH usage (P<0.03), as well as different clinical features and outcome in terms of time to first treatment (TTT) and overall survival (P<0.001). A multivariate Cox model showed that epigenetic classification was the strongest predictor of TTT (P<0.001) along with Binet stage (P<0.001). These findings were corroborated in a validation series (n=97). In this study, we developed a simple and robust method using epigenetic biomarkers to categorize CLLs into three subgroups with different clinicobiologic features and outcome.

Recurrent mutation of JAK3 in T-cell prolymphocytic leukemia.

T‐cell prolymphocytic leukemia (T‐PLL) is an aggressive post‐thymic T‐cell malignancy characterized by the recurrent inv(14)(q11q32)/t(14;14)(q11;q32) or t(X;14)(q28;q11) leading to activation of either the TCL1 or MTCP1 gene, respectively. However, these primary genetic events are insufficient to drive leukemogenesis. Recently, activating mutations in JAK3 have been identified in other T‐cell malignancies. Since JAK3 is essential for T‐cell maturation, we analyzed a cohort of 32 T‐PLL patients for mutational hot spots in the JAK3 gene using a step‐wise screening approach. We identified 14 mutations in 11 of 32 patients (34%). The most frequently detected mutation in our cohort was M511I (seen in 57% of cases) previously described as an activating change in other T‐cell malignancies. Three patients carried two mutations in JAK3. In two patients M511I and R657Q were simultaneously detected and in another patient V674F and V678L. In the latter case we could demonstrate that the mutations were on the same allele in cis. Protein modeling and homology analyses of mutations present in other members of the JAK family suggested that these mutations likely activate JAK3, possibly by disrupting the activation loop and the interface between N and C lobes, increasing the accessibility of the catalytic loop. In addition, four of the 21 patients lacking a JAK3 point mutation presented an aberrant karyotype involving the chromosomal band 19p13 harboring the JAK3 locus. The finding of recurrent activating JAK3 mutations in patients with T‐PLL could enable the use of JAK3 inhibitors to treat patients with this unfavorable malignancy who otherwise have a very poor prognosis.

eFORGE: A Tool for Identifying Cell Type-Specific Signal in Epigenomic Data

Summary Epigenome-wide association studies (EWAS) provide an alternative approach for studying human disease through consideration of non-genetic variants such as altered DNA methylation. To advance the complex interpretation of EWAS, we developed eFORGE (http://eforge.cs.ucl.ac.uk/), a new standalone and web-based tool for the analysis and interpretation of EWAS data. eFORGE determines the cell type-specific regulatory component of a set of EWAS-identified differentially methylated positions. This is achieved by detecting enrichment of overlap with DNase I hypersensitive sites across 454 samples (tissues, primary cell types, and cell lines) from the ENCODE, Roadmap Epigenomics, and BLUEPRINT projects. Application of eFORGE to 20 publicly available EWAS datasets identified disease-relevant cell types for several common diseases, a stem cell-like signature in cancer, and demonstrated the ability to detect cell-composition effects for EWAS performed on heterogeneous tissues. Our approach bridges the gap between large-scale epigenomics data and EWAS-derived target selection to yield insight into disease etiology.

Genes encoding members of the JAK-STAT pathway or epigenetic regulators are recurrently mutated in T-cell prolymphocytic leukaemia.

T‐cell prolymphocytic leukaemia (T‐PLL) is an aggressive leukaemia. The primary genetic alteration in T‐PLL are the inv(14)(q11q32)/t(14;14)(q11;q32) leading to TRD/TRA‐TCL1A fusion, or the t(X;14)(q28;q11) associated with TRD/TRA‐MTCP1 fusion. However, additional cooperating abnormalities are necessary for emergence of the full neoplastic phenotype. Though the pattern of secondary chromosomal aberrations is remarkably conserved, targets of the changes are largely unknown. We analysed a cohort of 43 well‐characterized T‐PLL for hotspot mutations in the genes JAK3, STAT5B and RHOA. Additionally, we selected a subset of 23 T‐PLL cases for mutational screening of 54 genes known to be recurrently mutated in T‐cell and other haematological neoplasms. Activating mutations in the investigated regions of the JAK3 and STAT5B genes were detected in 30% (13/43) and 21% (8/39) of the cases, respectively, and were mutually exclusive. Further, we identified mutations in the genes encoding the epigenetic regulators EZH2 in 13% (3/23), TET2 in 17% (4/23) and BCOR in 9% (2/23) of the cases. We confirmed that the JAK‐STAT pathway is a major mutational target, and identified epigenetic regulators recurrently mutated in T‐PLL. These findings complement the mutational spectrum of secondary aberrations in T‐PLL and underscore the potential therapeutical relevance of epigenetic regulators in T‐PLL.

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.

DNA methylation profiling of pediatric B-cell lymphoblastic leukemia with KMT2A rearrangement identifies hypomethylation at enhancer sites.

Deregulation of the epigenome is an important pathogenetic mechanism in acute lymphoblastic leukemia (ALL) with lysine (K)‐specific methyltransferase 2A rearrangement (KMT2Ar). We performed array‐based DNA methylation profiling of KMT2Ar ALL cells from 26 children in comparison to normal B‐cell precursors. Significant changes in DNA methylation in KMT2Ar ALL were identified in 2,545 CpG loci, influenced by age and the translocation partners AFF1 and MLLT1. In KMT2Ar ALL, DNA methylation loss was enriched at enhancers and for certain transcription factor binding sites such as BCL11A, EBF, and MEF2A. In summary, DNA methylation changes in KMT2Ar ALL target enhancers, genes involved in leukemogenesis and normal hematopoiesis, as well as transcription factor networks.

Recurrent alterations of TNFAIP3 (A20) in T-cell large granular lymphocytic leukemia

The pathogenesis of T‐cell large granular lymphocytic leukemia (T‐LGL) is poorly understood, as STAT3 mutations are the only known frequent genetic lesions. Here, we identified non‐synonymous alterations in the TNFAIP3 tumor suppressor gene in 3 of 39 T‐LGL. In two cases these were somatic mutations, in one case the somatic origin was likely. A further case harbored a SNP that is a known risk allele for autoimmune diseases and B cell lymphomas. Thus, TNFAIP3 mutations represent recurrent genetic lesions in T‐LGL that affect about 8% of cases, likely contributing to deregulated NF‐κB activity in this leukemia.

SAMHD1 IS FREQUENTLY INVOLVED IN T‐CELL PROLYMPHOCYTIC LEUKEMIA (T‐PLL) PATHOGENESIS

growth to Ba/F3 cells, increased not only p‐STAT3 but also PD‐L1 expression. Consistent with these findings, PD‐L1 was over expressed in cell lines harboring STAT3 mutations and a positive correlation between PD‐L1 and p‐STAT3 expression was observed in NKTL tumor tissue (R = 0.51, P = .02). Conclusions: We characterized a novel‐activating STAT3 mutant and demonstrated its ability to drive PD‐L1 expression, which may promote tumor evasion from the antitumor immune response. The combination of PD‐1/PD‐L1 antibodies and STAT3 inhibitors might be a promising and novel therapeutic approach for NKTL in the future.