EZH2 is required for germinal center formation and somatic EZH2 mutations promote lymphoid transformation (original) (raw)

. Author manuscript; available in PMC: 2014 May 13.

Published in final edited form as: Cancer Cell. 2013 May 13;23(5):677–692. doi: 10.1016/j.ccr.2013.04.011

Abstract

The EZH2 histone methyltransferase is highly expressed in germinal center (GC) B-cells and targeted by somatic mutations in B-cell lymphomas. Here we find that EZH2 deletion or pharmacologic inhibition suppresses GC formation and functions in mice. EZH2 represses proliferation checkpoint genes and helps establish bivalent chromatin domains at key regulatory loci to transiently suppress GC B-cell differentiation. Somatic mutations reinforce these physiological effects through enhanced silencing of EZH2 targets in B-cells, and in human B-cell lymphomas. Conditional expression of mutant EZH2 in mice induces GC hyperplasia and accelerated lymphomagenesis in cooperation with BCL2. GCB-type DLBCLs are mostly addicted to EZH2, regardless of mutation status, but not the more differentiated ABC-type DLBCLs, thus clarifying the therapeutic scope of EZH2 targeting.

INTRODUCTION

The majority of B-cell lymphomas arise from germinal center (GC) B-cells, which form part of the humoral immune response to T-cell dependent antigen stimulation (Ci et al., 2008; Klein and Dalla-Favera, 2008). GC B-cells originate from mature, VDJ rearranged, antigen-naïve B-cells present in lymphoid follicles that upon activation migrate to the center of lymphoid follicles and form GCs. GC B-cells are unique in their ability to replicate at an accelerated rate while undergoing somatic hypermutation, which requires attenuation of DNA damage sensing and replication checkpoints. Clonal expansion of these cells is followed by terminal differentiation of the subset of B-cells that encode high affinity antibodies; by contrast the remainder of GC B-cells undergo apoptosis. The rapid proliferation of GC cells with concomitant attenuation of DNA damage responses increases the risk of oncogenic mutations as a byproduct of somatic hypermutation and class switch recombination. Such mutations can lead to the development of diffuse large B-cell lymphoma (DLBCL) and follicular lymphoma (FL) (Alizadeh et al., 2000; Klein and Dalla-Favera, 2008). Understanding how B-cells impose and maintain the GC phenotype may thus provide important clues that could explain the pathogenesis of GC B cell-derived lymphomas and inform the design of rational therapeutic strategies for patients with DLBCL and FL.

Upon activation, GC B-cells upregulate and highly express EZH2, a SET domain containing histone methyltransferase that forms part of Polycomb Repressive Complex-2 (PRC2) (Raaphorst et al., 2000; Velichutina et al., 2010). EZH2 catalyzes methylation of the lysine 27 residue of histone 3 (H3K27), but is mostly catalytically active when complexed with the PRC2 components EED and SUZ12 (Chase and Cross, 2011; Mahmoudi and Verrijzer, 2001). H3K27 methylation is a repressive histone mark associated with gene repression (Cao et al., 2002; Czermin et al., 2002; Muller et al., 2002). EZH2 is critical for normal development and EZH2 knockout embryos fail to implant and undergo gastrulation (O’Carroll et al., 2001). EZH2 plays a major role in regulating gene expression patterning in embryonic and tissue specific stem cells, where it contributes to setting genes into a “poised” bivalent state characterized by the simultaneous appearance of the H3K27me3 repressive mark along with the H3K4me3 activation-associated mark (Bernstein et al., 2006; Shin et al., 2012). These poised bivalent genes can then be resolved into activated or stably repressed loci depending on lineage commitment. In tissue-specific stem cells, EZH2 and partner proteins are critical for the proper coordination of differentiation and proliferation. For example, conditional knockout studies in early stages of B-cell differentiation showed that EZH2 is important for normal immunoglobulin VDJ recombination in pre-B-cells (Su et al., 2003). However after the pre-B cell stage EZH2 expression declines and is not detectable in mature B-cells located within lymphoid tissues, until these cells enter the GC reaction (Su et al., 2003; van Galen et al., 2004; Velichutina et al., 2010). Once B-cells exit the GC reaction EZH2 is once again downregulated (Velichutina et al., 2010).

We previously found that EZH2 target genes in GC B-cells only partially overlap with EZH2 targets in embryonic stem cells, suggesting that it has GC specific functions (Velichutina et al., 2010). Moreover PRC2 components including EZH2 are often highly expressed in GC B-cell derived DLBCL (van Kemenade et al., 2001), suggesting that PRC2 functions may contribute to lineage-based proliferation of malignant B-cells. RNA-sequencing studies recently identified heterozygous somatic point mutations targeting the EZH2 SET domain, most commonly affecting the Y641 residue of EZH2, in ~30% of GCB-type DLBCL and 10% of FL (Morin et al., 2010). These mutations alter the enzymatic activity of EZH2 resulting in a protein that fails to recognize unmodified H3K27 and preferentially converts mono or dimethylated H3K27 to the trimethylated state (Sneeringer et al., 2010; Yap et al, 2011). DLBCL cell lines harboring EZH2 point mutations display aberrant accumulation of H3K27me3 compared to WT-EZH2 cells (McCabe et al., 2012a; Sneeringer et al., 2010). The importance of EZH2 in maintaining the growth of lymphoma cells was demonstrated in siRNA knockdown experiments (Velichutina et al., 2010) and more recently through the development of highly specific and potent EZH2 small molecule inhibitors (Knutson et al., 2012; McCabe et al., 2012b; Qi et al., 2012). However the role and requirement for EZH2 and its mutant alleles in the formation of GCs and pathogenesis of B-cell lymphomas has not been explored, nor is its mechanism of transformation known beyond the association with global cellular abundance of H3K27me3. Here we explore the function of normal and mutant EZH2 in mature B-cells and lymphomas using murine models and human cells. We also address the mechanisms of action of WT and mutant EZH2 and provide the basis for rational combinatorial targeted therapy for GCB-type DLBCL.

RESULTS

EZH2 Is Required For Formation Of Germinal Centers

The increased expression of EZH2 in GC B-cells suggested that this protein plays critical roles in mediating the phenotype of GC B-cells. However, EZH2 deletion is lethal in early embryonic development (O’Carroll et al., 2001), and its inducible knockout in early hematopoietic cells perturbs lymphoid differentiation at the pre-B cell stage (Su et al., 2003). We therefore crossed conditional _Ezh2_−/− mice with the Cγicre strain, which expresses cre recombinase in GC B-cells (Casola et al., 2006). Upon reaching immunological maturity, Ezh2fl/fl, Ezh2−/− and WT control mice were injected with T-cell dependent antigen sheep red blood cells (SRBC), to induce GC formation, and sacrificed 10 days later, at which time the GC reaction is at its peak. EZH2 loss resulted in a marked reduction in the number of splenic GC (GL7+/FAS+/B220+) B-cells (Figure 1A). Immunohistochemical analysis using peanut agglutinin (a GC B-cell marker) revealed a reduction in the number (p<0.0001) and size (p=0.001) of GCs in _Ezh2_−/− vs. Ezh2+/+mice, whereas there was no change in GCs in _Ezh2_fl/fl mice (Figure 1B-1C). There was also marked reduction in Ki67-positive cells consistent with loss of the proliferative GC B-cell compartment (Figure 1B). The histologic appearance of the spleen and primary lymphoid follicles was otherwise normal (Figure S1A) and marginal zone and follicular B-cells were unaffected (Figure S1B). The few residual GC B-cells were EZH2 positive, consistent with incomplete cre-mediated excision of Ezh2 and the notion that GC B cells require EZH2 (data not shown). A similar phenotype was observed when _Ezh2_−/− mice were injected with NP-KLH (nitrophenyl-keyhole limpet hemocyanin), which induces cognate high affinity antibodies (Figure S1C). The defect in GC formation was accompanied by impaired immunoglobulin affinity maturation, with reduction in formation of high affinity antibodies (Figure S1D).

Figure 1. EZH2 Is Required for Germinal Center Formation.

Figure 1

(A-C) EZH2fl/fl, EZH2−/− and WT C57BL6 control mice (n=5 per group) were immunized with SRBC to induce GC formation. (D-F) C57BL6 mice were immunized with SRBC and treated with GSK503 (150 mg/kg/day) or vehicle (n=5 per group). (A and D) Representative flow cytometric plot of splenic GC cells and quantification. (B and E) Splenic tissue was stained with PNA, EZH2, Ki67 and B220. (C and F) Quantitative imaging of PNA staining from (B) and (E) respectively. Values in A, C, D and F are shown as mean±SEM. t test ***p<0.001, **p<0.01, *p<0.05. See also Figure S1 and Tables S1-S4.

We next immunized a cohort of C57BL6 mice with SRBC followed by once daily treatment with 150 mg/kg/day GSK503, a specific EZH2 methyltransferase inhibitor or vehicle. GSK503 was designed by performing high-throughput screening for EZH2 methyltransferase inhibiting compounds followed by chemical optimization (Figure S1E). GSK503 inhibits the methyltransferase activity of WT and mutant EZH2 with similar potency (Kiapp=3−27nM) and is structurally related to GSK126 and GSK343 (McCabe et al., 2012b; Verma et al., 2012). GSK503 was > 200 fold selective over EZH1 (Kiapp=636nM) and > 4000 fold selective over other histone methyltransferases (Tables S1-S4; Figure S1F-S1G). GSK503 displayed favorable pharmacokinetics in mice (Figure S1H-S1I) and was thus selected for use in in vivo experiments. We assessed the effects of GSK503 mediated EZH2 inhibition in vivo in GC B-cells. GSK503, but not vehicle, prevented the formation of GC after SRBC or NP-KLH immunization, phenocopying the Ezh2 null phenotype. GSK503 treatment led to reduced numbers of GC B-cells by flow cytometry, reduced number and volume of GCs by immunohistochemistry, and impaired formation high affinity antibodies (Figure 1D-1F; Figure S1J-S1K). There was no effect on marginal zone or follicular B-cells (Figure S1L). Immunoblot analysis revealed that GSK503 reduced the level of H3K27me3 in splenocytes, confirming inhibition of the catalytic function of EZH2 in vivo (Figure S1M). Collectively these data demonstrate that EZH2 is required for the formation of GCs and immunoglobulin affinity maturation, and that this function is dependent on its histone methyltransferase activity.

Mutant EZH2 Induces Germinal Center Hyperplasia

Given that EZH2 was essential for development of GC B cells, we surmised that gain of function EZH2 mutants might reinforce these actions, with implications for lymphomagenesis. To address this point we generated mice conditionally expressing the Ezh2Y641N lymphoma allele from the ColA1 locus after excision of a Lox-Stop-Lox cassette (Figure S2A-S2C). These animals were crossed into the Cγicre background to express _Ezh2_Y641N in GC B-cells. Subsequent to immunization these mice displayed greatly increased numbers of GL7+/FAS+/B220+ GC B-cells (p<0.001, Figure 2A), and increased number and size of GCs (p<0.0005, Figure 2B-2C), while cre-negative _Ezh2_fl/fl animals showed no such effects. There was no effect on marginal zone or follicular B-cell numbers or morphology (Figure S2D). The phenotype was accompanied by a three-fold increase in the abundance of H3K27me3 in sorted GC B-cells from Ezh2Y641N vs. non-recombined mice (p<0.005, Figure 2D), analogous to what is observed in mutant EZH2 DLBCL cell lines (McCabe et al., 2012a; Sneeringer et al., 2010). Lymphoid GC hyperplasia was dependent on the increased histone methyltransferase activity as GSK503 treatment abrogated the Ezh2Y641N GC hyperplasia phenotype (p<0.001 Figure 2A-2C). These data demonstrate that Y641 mutation reinforces the GC phenotype-driving function of wild type EZH2, resulting in an expansion of these proliferative and mutagenic cells. However, Ezh2Y641N knock-on mice did not develop B-cell lymphomas (data not shown).

Figure 2. Mutant EZH2 Induces Germinal Center Hyperplasia.

Figure 2

EZH2fl/fl and EZH2Y641N mice were immunized with NP-KLH and treated with GSK503 (150 mg/kg/day) or vehicle (n=6 per group). (A) Representative flow cytometric plot of splenic GC cells and quantification. (B) Splenic tissue was stained with PNA, EZH2, Ki67 and B220. (C) Quantitative imaging of PNA from (B). (D) Immunoblotting from whole cell lysates from sorted splenic GC B cells (GL7+/FAS+/B220+) Values in A, C and D are shown as mean±SEM. t test ***p<0.001, **p<0.01. See also Figure S2.

Mutant EZH2 Enhances Proliferation Effects in part through Enhanced Repression of CDKN1A

To characterize in more detail the functions of EZH2 in mature B-cells we next expressed wild type EZH2 and EZH2Y641 in murine BCL1 cells. BCL1 cells originate from a spontaneous murine lymphoproliferative disease, exhibit clonal somatic hypermutation, and have been used to model GC B-cell biology (Blackman et al., 1986). BCL1 cells were transduced with GFP expressing retrovirus harboring FLAG-tagged Ezh2Y641N, Ezh2Y641F, or WT Ezh2 (Figure 3A). BCL1 cells transduced with EZH2Y641N and EZH2Y641F exhibited increased H3K27me3 levels whereas H3K27me3 levels were unchanged by WT EZH2 (Figure 3A-3B). EZH2Y641N and EZH2Y641F, but not WT EZH2 conferred a significant growth advantage to BCL1 cells compared to vector infected cells (p<0.05, Figure 3C). Moreover EZH2Y641N and EZH2Y641F induced >10-fold expansion in colony formation as compared to WT EZH2 or vector (Figure 3D). Treatment with the specific EZH2 inhibitor GSK343 (Verma et al., 2012 and Figure S1E-S1G) but not vehicle or the inactive chemical control GSK669 depleted H3K27me3 in BCL1 cells (Figure 3E) and abrogated colony formation induced by EZH2Y641N or EZH2Y641F (Figure 3F), confirming that enzymatic activity is required for their biological effects.

Figure 3. Mutant EZH2 Enhances Proliferation Effects Through Repression of CDKN1A.

Figure 3

(A) Immunoblotting from nuclear extracts from GFP+ BCL1 cells expressing GFP-FLAG-tagged EZH2Y641F, EZH2Y641N, or WT EZH2. (B) Mass spectrometry measurement of H3K27me3 mark. (C) Viability using trypan blue exclusion. (D and F) Colony counts and representative pictures. (E) Immunoblotting was performed as in (A) with cells treated as indicated for 3 days. (G) Immunoblotting from whole cell lysates from YFP+ DLBCL cells expressing 2 independent EZH2 shRNA or control (H) Viability of DLBCL YFP+ by flow cytometry using annexinV and DAPI exclusion. Green, WT EZH2; blue, mutant EZH2 cell lines. (I) RT-qPCR of CDKN1A in DLBCL cell lines treated with 2μM GSK for 7 days. Green and blue code as in (H). (J) Pfeiffer cells were co-transduced with the indicated shRNAs and viable YFP+ cells counted as in (H). (K) Pfeiffer cells expressing 2 independent shCDKN1A or control were treated with 2μM GSK for 6 days and viability was measured as in (H). Values in B, C, D, F, H, I, J and K are mean of duplicate or triplicate ±SD. t test ***p<0.001, **p<0.01, *p<0.05. EV, empty vector; WT, WT EZH2. See also Figure S3.

Given that BCL1 cells express endogenous EZH2 already (Figure 3A), and expression of exogenous EZH2 did not affect proliferation or colony formation in BCL1 cells, the effects we observed were specific to expression of mutant EZH2 in lymphoid cells. Recent studies have shown that EZH2 inhibitors and RNAi suppress proliferation of EZH2 mutant DLBCL cells (Knutson et al., 2012; McCabe et al., 2012b; Qi et al., 2012; Velichutina et al., 2010). Consistent with these findings, shRNA depletion of EZH2 protein (Figure 3G) in EZH2 mutant Pfeiffer and WSU-DLCL2 cell lines led to >50% suppression of cell growth relative to shRNA control. shRNA mediated depletion of EZH2 in WT EZH2 Farage cells led to a more modest and delayed 20-40% decrease of growth (Figure 3H). GSK343 2 μM even more efficiently depleted H3K27me3 and more rapidly suppressed the growth of DLBCL cell lines, with greater effect observed in those with mutant EZH2 (Figure S3A-S3B). We identified CDKN1A as a direct EZH2 target in GC B-cells and DLBCL (Velichutina et al., 2010)and data not shown). Treatment of WT EZH2 cells with GSK343 led to a 2-4 fold induction of CDKN1A while in mutant cell lines EZH2 inhibition led to 8-12 fold activation (Figure 3I). While EZH2 shRNA or GSK343 targeting of EZH2 led to a reduction in DLBCL cell viability, shRNA mediated depletion of CDKN1A partially rescued this phenotype (p<0.05, Figure 3J-3K; Figure S3C). We performed gene expression profiling in diagnostic samples from patients with DLBCL and observed a trend for CDKN1A to be expressed at a lower level in GCB-DLBCL patients harboring EZH2 mutations vs. GCB-DLBCL with wild type EZH2 (p=0.07, Figure S3D). These data are consistent with the notion that mutant EZH2 reinforces repression of GC B-cell EZH2 target genes, and that CDKN1A silencing may contribute along with other genes to the phenotypic effects of EZH2 in B-cells.

EZH2 Mediates Differentiation Blockade in DLBCL Cells

GCB-DLBCLs exhibit differentiation blockade locking them into the GC phenotype (Rui et al., 2011). Given that gain-of-function mutants induce an expansion in the population of GC B-cells we postulated that EZH2 might play a critical role in suppressing differentiation. We first compared RNA-seq gene expression profiles in BCL1 cells transduced with WT or mutant EZH2. This analysis revealed that in BCL1 cells EZH2Y641N or EZH2Y641F increased repression of genes that are normally expressed upon terminal differentiation into plasma and memory cells when compared to EZH2-WT expressing cells (Figure 4A). We induced differentiation in BCL1 cells using IL2 and IL5. These cytokines typically drive BCL1 cells towards a partial plasma cell phenotype accompanied by proliferation arrest (Blackman et al., 1986). Expression of EZH2Y641N or EZH2Y641F but not WT-EZH2 significantly impaired cytokine-induced proliferation arrest (p<0.05, Figure 4B) and more profoundly suppressed activation of Prdm1, a key plasma cell differentiation gene (p<0.01, Figure 4C).

Figure 4. EZH2 Mediates Differentiation Blockade.

Figure 4

(A) GSEA in WT vs. mutant EZH2 BCL1 transduced cells. (B) Cells were treated with 20ng/mL IL2 and IL5 for 5 days and viability measured using trypan blue exclusion. (C) RT-qPCR of Prdm1. (D) Cells were treated with 2μM GSK343 or control GSK669 for 7 days and Ig expression levels were examined by flow cytometry. (E) Representative pictures of cells treated as in (D). (F) RT-qPCR of indicated mRNAs in cells treated with 2μM GSK for 7 days. Green, WT EZH2; blue, mutant EZH2 cell lines. (G) Heat map of over-represented gene categories among genes up-regulated by EZH2 shRNA or 2μM GSK343 for 7 days in Ly7, Ly1, SUDHL5, Farage, WSU-DLCL2 and Pfeiffer cells. Enrichment measured using hypergeometric p-values. (H) GSEA in GCB-DLBCL patient samples with WT vs. mutant EZH2. Values in B, C, D and F are mean of duplicate or triplicate ±SD. t test ***p<0.001, **p<0.01, *p<0.05. See also Figure S4 and Tables S5-S6.

Based on these data we next examined whether EZH2 inhibition by GSK343 induces differentiation in human DLBCL cells. In GCB-type DLBCL cell lines (two mutant and one WT), EZH2 inhibitor treatment induced upregulation of the plasma cell proteins Syndecan/CD138 and p63, upregulation of the memory cell marker CD27, and downregulation of the B-cell surface marker CD20 (Figure S4A). Treatment with GSK343, but not GSK669, induced functional evidence of differentiation manifested by increased light chain and immunoglobulin production (Figure 4D), and morphologic changes consistent with plasma cell differentiation (Figure 4E). Gene expression analysis by QPCR revealed upregulation of differentiation related genes including PRDM1, IRF4 and syndecan CD138 (Figure 4F) with greater effect generally noted in EZH2 mutant DLBCL cells (p<0.05). Analysis of RNA-seq expression profiles induced by GSK343 in three DLBCL cell lines (2 mutant and 1 WT EZH2) showed enrichment for transcriptional programs involved in plasma and memory B cell differentiation (FDR q<0.001, Figure S4B). We also found strong enrichment in gene sets involved in exit from the GC reaction including genes induced by CD40, IRF4, IL10, and NFkB, and notably the ABC-DLBCL signature (which consists of activation of these post-GC genes), as well as genes involved in apoptosis and immune responses (Figure 4G). Memory and plasma cell differentiation genes were among genes preferentially repressed in GCB-type DLBCL patients with EZH2 mutations as compared to GCB-DLBCL patients with WT EZH2 (FDR q<0.001, Figure 4H; Tables S5-S6). IRF4 was expressed at a significantly lower level in EZH2 mutant vs. WT GCB-DLBCL patients (p=0.01, Figure S4C). Collectively these data suggest that EZH2 maintains the GC phenotype by suppressing transcriptional programs required for exiting the GC reaction and terminal differentiation. Mutant EZH2 augments these same functions by reinforcing repression of B-cell differentiation genes.

Mutant EZH2 Aberrantly Represses Genes Through Increased Promoter H3K27 Trimethylation

In order to gain insight into how mutant EZH2 might reinforce or alter normal EZH2 dependent transcriptional programs we performed H3K27me3 ChIP-seq in BCL1 cells after transduction with EZH2Y641N, EZH2Y641F, WT EZH2 or vector and determined how gene expression in these cells responded to EZH2 inhibitors . Cells transduced with either EZH2Y641N or EZH2Y641F showed a marked increase in the abundance of H3K27me3 at gene promoters compared to EZH2 WT transduced cells (p<2×10−16 for both mutants, Figure 5A) and genes with abundant H3K27me3 were generally repressed (p<2×10−16, data not shown). In addition, loci displaying increased H3K27me3 were more potently induced upon exposure of the BCL1- EZH2Y641N/EZH2Y641F cells to GSK343 than EZH2 WT cells (p<5x10−11 for EZH2Y641N and EZH2Y641F Figure 5B). Remarkably, those genes with increased H3K27me3 due to mutant EZH2 and more strongly induced by GSK343 corresponded to the set of human EZH2 target genes that acquire de novo H3K27me3 in human GC B-cells (FDR q<0.001, Figure 5C). Hence mutant EZH2 appears to function at least in part by exaggerating the epigenetic silencing of normal GC B-cell targets of wild type EZH2.

Figure 5. Mutant EZH2 Aberrantly Represses Genes by Increasing H3K27me3 at Promoters and EZH2 Generates GC B-cell Specific Bivalent Genes Involved in Differentiation.

Figure 5

(A) ChIP-seq density plot. (B) Heat map of change in expression of genes with increased H3K27me3 within promoters in BCL1 cells treated with 0.5μM GSK for 3 days. (C and G) GSEA in WT vs. mutant EZH2 BCL1 cells. (D) State of GCB cell bivalent promoters within naïve B cells (NBC). (E) Density strip representations of normalized EZH2 ChIP-seq reads within DLBCL cell line promoters. (F) Gene expression level of bivalent, H3K4me3 monovalent, H3K27me3 monovalent, and promoters not exhibiting either of these marks in GCB cells. (H) GSEA in GCB-DLBCL patient samples with WT vs. mutant EZH2. (I) Heat map of relative expression of GCB cell bivalent genes repressed in WT and mutant EZH2 GCB-DLBCL patient samples. See also Figure S5.

EZH2 is Linked to Formation of GC B-Cell Specific Bivalent Genes Involved in Differentiation

In stem cells a subset of genes that are H3K27 trimethylated by EZH2 are also marked by the activating chromatin modification H3K4me3. These so-called “bivalent (H3K4me3/H3K27me3) domains” are believed to represent a mechanism for “poising” key lineage transcription factors so that they can be either activated or repressed during subsequent differentiation (Bernstein et al., 2006). We examined H3K27me3 and H3K4me3 distribution in purified primary human naïve B (NB)-cells and GC B-cells by ChIP-seq. Notably, even though GC B-cells are mature, committed cells (i.e. in contrast to stem cells), we find that they gain 1026 new bivalent domains at promoters that are not found in NB cells (Figure 5D). Of these 910 (88%) bivalent promoters originate from H3K4me3 only promoters in NB cells (Figure 5D), indicating that most bivalent loci in GC B-cells occur due to acquisition of H3K27me3, concordant with upregulation of EZH2 in these cells. ChIP-seq profiles from seven human GCB-DLBCL cell lines showed that these bivalent promoters were preferentially bound by EZH2 (p<1×10−16, Figure 5E). In primary human GC B-cells bivalent genes were expressed at lower levels than loci with H3K4me3 alone, but at higher levels than genes marked only by H3K27me3 (p<1×10−300, Figure 5F). This is consistent with the proposed “poised” nature of bivalent genes in stem cells. GC B-cell bivalent genes were highly enriched in gene sets associated with termination of the GC reaction such as IRF4 induced, ABC-DLBCL, memory cell upregulated genes, as well as other GC B-cell relevant genes such as negative regulation of cell cycle (Figure S5A). The key regulatory transcription factors IRF4 and PRDM1 were among the GC B-cell specific bivalent genes, consistent with the notion that these genes are marked for dynamic activation upon exit from the GC reaction (Figure S5B) and that bivalent marks at this stage of B-cell maturation function to transiently suppress terminal differentiation.

We wondered whether somatic mutation of EZH2 might lock these bivalent promoters into a more repressed configuration, perhaps helping to explain the irreversible differentiation blockade induced by mutant EZH2. Indeed we observed that GC B-cell bivalent genes were even more highly repressed in murine BCL1 cells transduced with EZH2Y641N and EZH2Y641F than in WT EZH2 transduced cells (FDR q<0.001, Figure 5G). GC B-cell bivalent domain genes were also significantly more repressed in human GCB-DLBCL patients with EZH2 somatic mutations compared to GCB-DLBCL with WT EZH2 (FDR q<0.001, Wilcoxon p<0.0006 Figure 5H-5I). GC B-cell bivalent genes involved in GC B-cell differentiation such as IRF4 induced genes, CD40 induced genes and plasma/memory cell genes were especially enriched among genes differentially expressed in mutant EZH2 DLBCL patients vs. GCB-DLBCL specimens with wild type EZH2 (Figure S5C-D). GC B-cell bivalent genes were also significantly enriched among genes induced by GSK343 treatment of individual human DLBCL cell lines (GSEA FDR q<0.001 Figure S5E, and hypergeometric distribution test, p=3.4×10−48). ChIP-seq showed enrichment of H3K27me3 at GCB bivalent genes in DLCBL cell lines (p<2×10−16, data not shown). To confirm that bivalent marks are indeed occurring at the same chromatin regions within GC derived DLBCL cells, we performed ChIP re-ChIP assays. The key regulatory transcription factors necessary for GC exit, IRF4 and PRDM1, as well as the proliferation checkpoint gene promoter CDKN1B, all shown to be putative GC B-cell bivalent genes by our ChIP-seq experiments, were confirmed to be significantly co-occupied by H3K27me3 and H3K4me3 marks (Figure S5F). In contrast HOXA7 (which is silenced in mature B-cells) was exclusively H3K27 trimethylated and BCL6 (which is highly expressed in GC B-cells) was almost exclusively H3K4 trimethylated. Collectively, these data suggest that in normal GC B-cell development EZH2 reversibly suppresses terminal differentiation by forming or maintaining bivalent domains at specific loci that are also marked by H3K4me3, in addition to epigenetically silencing other genes through pure H3K27me3. In B-cell lymphomas mutant EZH2 reinforces silencing of these genes, perhaps by increasing the frequency with which these genes are H3K27methylated among tumor cells, or by more subtle stoichiometric effects tipping the balance of H3K27 towards trimethylation and away from demethylation.

EZH2 Cooperates With BCL2 to Generate Germinal Center Derived Lymphomas

We next investigated whether mutant EZH2 might cooperate with other GC B-cell lymphoma oncoproteins, such as BCL2, which is frequently translocated in patients with EZH2Y641 mutations (Morin et al., 2011). We transduced bone marrow of vav-Bcl2 mice with retrovirus expressing GFP and encoding WT EZH2, EZH2Y641F or GFP alone, and transplanted them into lethally irradiated cohorts of ten recipients each (Figure 6A). Animals were immunized with SRBC every four weeks to ensure continuous formation of GCs. Macroscopic examination of spleens showed marked splenomegaly in EZH2Y641F vs. EZH2 WT or empty vector (Figure 6B). Immunoblot analysis of splenic extracts showed similar levels of expression of EZH2 in the EZH2 WT and EZH2Y641F mice (Figure 6C). By contrast H3K27me3 abundance was only elevated in splenocytes from EZH2Y641F mice, whereas H3K27me3 levels in EZH2 WT mice were similar to controls (Figure 6C). The livers of EZH2Y641F but not EZH2 WT mice were also significantly enlarged vs. vector control (Figure 6D).

Figure 6. EZH2 Cooperates With BCL2 To Generate Germinal Center Derived Lymphomas.

Figure 6

(A) Bone marrow transplantation was performed using VavP-Bcl2 transgenic donor mice. (B) Representative pictures of spleens from mice sacrificed 111 days after transplantation and quantification of the spleen weight (n=10 per group). (C) Immunoblotting from whole cell lysates from splenocytes of transplanted mice. (D) Quantification of the liver weight of transplanted mice (n=10 per group). (E) Splenic tissue from transplanted mice was stained with H&E, B220 and EZH2. (F) Tumor clonality analysis by RT-PCR performed in B220+ sorted splenocytes. Sample lanes separated by thin white lines were run on the same gel but were noncontiguous. (G) Dose reduction plot for Obatoclax and ABT737 at GI90 after exposure of cells to increasing concentrations of GSK343 for 6 days. Data represent mean of triplicate experiments. (H) Area under the curve (AUC) of the tumor growth curves for 20 days in SUDHL4 and SUDHL6 xenografted mice treated with vehicle (n=7), Obatoclax (2 mg/kg/d, n=8 in SUDHL4 and n=7 in SUDHL6), GSK503 (150 mg/kg/d, n=8 in SUDHL4 and n=7 in SUDHL6), or the combination of Obatoclax and GSK503 (n=8 in SUDHL4 and n=7 in SUDHL6). P values were calculated by t test. Values in B, C, D and H are mean±SEM. t test ***p<0.001, **p<0.01, *p<0.05. BM, bone marrow. See also Figure S6.

Detailed histopathologic examination and determination of clonality by immunoglobulin gene rearrangements indicated that whereas 70% of EZH2Y641F mice developed a B-cell lymphoma at 111 days, by contrast only 20% of the WT EZH2 and none of the EV mice showed evidence of lymphoma at this timepoint. All EV and WT EZH2 mice displayed evidence of follicular hyperplasia as expected in Bcl2 transgenic mice (Figure 6E). However examination of spleens in EZH2Y641F/Bcl2 mice revealed disruption of splenic architecture by neoplastic-appearing B220+ B-cells, most of which were EZH2 positive (Figure 6E; Figure S6A). Neoplastic B cells were large (3-4x the size of reactive lymphocytes), with pleomorphic nuclear morphology, fine chromatin and several nucleoli, ~60% Ki-67 expression, and were similar morphologically to human DLBCLs with centroblastic morphology (Figure S6A-S6B).

EZH2Y641F/Bcl2 mice displayed extensive infiltration of multiple tissues by neoplastic B cells, including lung, liver, lymph nodes, thymus, peri-cardiac soft tissue, pancreas, kidneys, salivary glands, urinary bladder, thyroid, ovarian bursa, bone marrow, skeletal muscles, small intestine, stomach (Figure S6C-S6D and data not shown). Histologic examination of lungs revealed marked peri-vascular and peri-bronchiolar and mild interstitial infiltration of neoplastic B-cells (Figure S6C). In the liver infiltration was also extensive with large periportal and to lesser extent perivenular and canalicular infiltrates (Figure S6D). Lymph nodes in EZH2Y641F/Bcl2 mice were also infiltrated by neoplastic cells with effacement of normal architecture and a mixture of nodular and diffuse patterns (data not shown). PCR amplification of IgVH regions from purified B220+ splenocytes revealed a dominant band in the background of few very weak bands suggestive of a mostly monoclonal population of B-cells (Figure 6F). Overall the histopathological analysis was consistent with an aggressive DLBCL-like disease.

The two EZH2 WT animals that developed lymphoid malignancy included one with T-cell and one with B-cell disease, and manifested a more indolent phenotype including better preservation of the lymphoid architecture, minimal tissue infiltration in lung and liver, and were composed of plasma cells and reactive-appearing small lymphocytes (data not shown). A second cohort of _Bcl2_-transduced mice (EZH2Y641F n=12, WT EZH2 n=8 and GFP alone n=9) were observed for survival. EZH2Y641F/Bcl2 manifested an accelerated lethal phenotype, with deaths due to progressive lymphoma beginning at day 90, whereas WT and EV mice began to die at day 140 (Figure S6E). The EZH2Y641F mice displayed a trend towards reduced median survival (median survival EZH2Y641F: 203 days, WT EZH2: 252, EV: 243).

Finally, given that BCL2 and EZH2 can cooperate in lymphomagenesis we hypothesized that BH3 mimetic drugs such as Obatoclax and ABT737 that block BCL2 function might enhance the activity of EZH2 inhibitors. We exposed a panel of GCB DLBCL cells to increasing concentrations of GSK343 or GSK503, in combination with ABT737 or Obatoclax. In almost every cell line tested the concentration of GSK343 or GSK503 required to yield 50% growth inhibition was reduced when cells were concomitantly treated with BH3 mimetics (Figure 6G; Figure S6F). In order to determine the impact of EZH2 inhibitor combinatorial therapy in a preclinical model, we evaluated the action of GSK503 and Obatoclax alone or in combination in mice bearing human DLBCL cell line (SUDHL4 and SUDHL6) xenografts. Although both GSK503 and Obatoclax inhibited tumor growth alone, the combination of these inhibitors again more potently and significantly suppressed tumor xenograft growth (Figure 6H; Figure S6G). EZH2 mutants therefore enable, accelerate and maintain malignant transformation of GC B-cells in cooperation with BCL2.

EZH2 Targeted Therapy Preferentially Affects GCB But Not ABC DLBCL Cells

EZH2 gain-of-function somatic mutations are restricted to GCB type DLBCLs (Morin et al., 2010). Moreover whereas EZH2 is a critical mediator of the GC B-cell phenotype which is reflected by the phenotype of GCB-DLBCL, it also represses genes and pathways that drive the phenotype of ABC-DLBCLs (e.g. IRF4, NFkB pathway genes). We questioned whether these GCB specific functions would translate to a specific role for EZH2 in the pathogenesis and therapeutic targeting of GCB vs. ABC subtypes of DLBCL. Along these lines we noted that EZH2 target genes were significantly more repressed in GCB than in ABC DLBCLs, although even more significantly repressed when comparing ABC vs. mutant EZH2 DLBCLs (Figure 7A). Gene set enrichment analysis revealed that EZH2 target genes (Figure 7B-7C), including GCB bivalent genes (Figure S7A) were overrepresented among genes upregulated in ABC vs. GCB DLBCLs (FDR q<0.001), with enrichment scores even higher when comparing ABC vs. mutant EZH2 DLBCLs (Figure S7A). Hence EZH2 target genes, including those with bivalent marks at the GC stage of development are expressed at relatively higher levels in ABC DLBCL cells, suggesting that EZH2 does not play a key role in their regulation in this form of lymphoma associated with the transformation of a post-germinal center B cell.

Figure 7. EZH2 Targeted Therapy Preferentially Affects GCB But Not ABC DLBCL Cells.

Figure 7

(A) Heat map of relative expression of GCB-specific H3K27me3 target genes repressed in mutant EZH2 GCB-DLBCL patient samples. (B-C) GSEA in ABC-DLBCL vs. WT EZH2 GCB-DLBCL (B) and mutant EZH2 GCB-DLBCL (C). (D-E) 4 ABC-DLBCL (HBL-1, Ly3, U2932, TMD8) and 9 GCB-DLBCL cell lines (5 WT EZH2: Ly7, Ly19, Farage, Ly18, SUDHL5, and 4 mutant EZH2: Pfeiffer, WSU-DLCL2, SUDHL4, SUDHL6) were exposed to increasing concentrations of GSK343 and GSK669 for 6 days. (D) Concentration of GSK343 required to inhibit 50% of growth (GI50) relative to GSK669. (E) Cell viability at 10μM GSK. Data are mean with 95% confidence interval for duplicate. (F) When naïve B-cells are activated EZH2 expression is highly induced. EZH2 is required for germinal center (GC) formation and Ig affinity maturation. EZH2 levels decrease as B cells exit the GC reaction, enabling expression of genes that mediate terminal differentiation (e.g. IRF4, PRDM1, NFkB). B-cells in this “exiting” phase such as plasmablasts are believed to give rise to ABC-DLBCLs. However, the occurrence of EZH2 somatic mutations aberrantly sustains repression of proliferation checkpoint and differentiation genes, resulting in GC hyperplasia, and the presence of other oncogenic hits, such as BCL2, enables transformation to GCB-type DLBCL or follicular lymphoma (FL). A possible alternative route leading to GCB-DLBCL could involve overexpression or aberrant maintenance of WT EZH2 expression. GCB-DLBCLs and FLs but not ABC-DLBCLs require EZH2 to maintain their proliferation and survival. Thus, EZH2 methyltransferase inhibitors suppress GCB-DLBCLs but not ABC-DLBCLs. See also Figure S7.

If EZH2 is not critical for the suppression of bivalent genes that drive proliferation once a cell exits the GC, then the biological effects of EZH2 inhibition would be predicted to be significantly different for ABC vs. GCB lymphoma cells and ABC cells might be predicted to be relatively insensitive to EZH2 inhibitors. Therefore, we treated a panel of ABC and GCB DLBCL cell lines with increasing concentrations of GSK343 or GSK503. Strikingly, the IC50 for GCB lymphoma cell lines were in the 0.5-20uM range, while for ABC cells no significant inhibition of cell growth was observed (p=0.0004, Figure 7D; Figure S7B-S7D). Exposure of EZH2 WT GCB cells to a fixed dose of 10μM GSK343 or GSK503 led to 30-75% reduction in viable cell number while EZH2 mutant GCB cells were inhibited 50-99% suggesting a trend toward increased sensitivity. By contrast this dose of drug led to no killing of ABC lymphoma cells (p=0.0001 vs GCB, Figure 7E; Figure S7E). EZH2 inhibitors completely demethylated H3K27me3 in ABC-DLBCL cells indicating that resistance is biological and not due to failure of the drug to inhibit its target (Figure S7F). Collectively these data suggest that EZH2 inhibitors may be useful for GCB DLBCL either with or without EZH2 mutations but are likely to be ineffective for ABC-type DLBCL.

DISCUSSION

In this manuscript we show that EZH2 is a master regulator of the GC B-cell phenotype, a function that is aberrantly reinforced by mutant EZH2 lymphoma disease alleles (Figure 7F). We find that EZH2 mediates its effects in GC B-cells by repressing target genes involved in proliferation checkpoints (e.g. CDKN1A) and exit from the GC and terminal differentiation (e.g. IRF4 and PRDM1). For immunoglobulin affinity maturation to occur, GC B-cells must maintain their phenotype long enough to transit repeated rounds of division and somatic hypermutation. Hence in the absence of EZH2 function, which is needed to support the GC B-cell phenotype, mice display defective immunoglobulin affinity maturation. Normally, EZH2 levels decrease as B cells exit the GC reaction, enabling expression of genes that mediate terminal differentiation (Velichutina et al., 2010). However in the presence of somatically mutated EZH2, suppression of GC exit genes and checkpoints persists resulting in hyperplasia, and the presence of other oncogenic hits enables transformation to GCB-type DLBCL (Figure 7F). An alternative route leading to GCB-DLBCL could involve overexpression or aberrant maintenance of WT EZH2. Indeed the highest quartile WT EZH2 expressing GCB-DLBCLs display a trend towards increased repression of EZH2 targets and tend to cluster together with EZH2 mutant patients (Figure S7G). The role of EZH2 in lymphomas is in large part to maintain or exaggerate (i.e. to hijack) the same EZH2 transcriptional program that is also required for development of normal GCs and immunoglobulin affinity maturation. Hence GC B-cells exhibit both oncogene and non-oncogene addiction to EZH2.

The actions of EZH2 appear linked in part to de novo formation of bivalent chromatin domains, whereby genes marked by H3K4me3 in naïve follicular B-cells acquire H3K27me3 in GC B-cells concordant with upregulation of EZH2. In stem cells bivalent (H3K4me3/H3K27me3) domains are hypothesized to maintain genes in a repressed but poised conformation, which can be subsequently dynamically activated or repressed according to lineage specific differentiation programs (Bernstein et al., 2006; Shin et al., 2012). GC B-cells represent a potentially unique situation where bivalent chromatin marks at specific genes are gained during differentiation; while in ES and tissue stem cell differentiation they are lost. We however cannot exclude that other cell types and tumors that may over-express EZH2 might also form bivalent domains at promoters. GC bivalent promoters are enriched in key gene sets involved in exit from the GC such as IRF4 and CD40 induced genes, and genes upregulated in memory and plasma cells. As long as EZH2 maintains H3K27me3 at these loci, these genes are expressed at low levels, maintaining the GC phenotype. Repression of these genes is terminated as B-cells exit the GC reaction and are selected for terminal differentiation. By contrast EZH2 mutants disrupt the equilibrium of bivalent domains, enabling aberrant, persistent epigenetic silencing of genes, which in turn allows the GC B-cell phenotype to persist, facilitating lymphoid transformation. Consistent with this hypothesis, we found that bivalent domain genes are the most aberrantly repressed in loci in EZH2 mutant DLBCL patient specimens. These data, along with the recent identification of frequent loss-of-function mutations in the histone methyltransferase protein MLL2 in B-cell lymphomas (Lohr et al., 2012; Morin et al., 2011), suggest that the balance of H3K4me3 and H3K27me3 is disrupted and represents a therapeutic target in lymphoid malignancies. Mutations in MLL2 and EZH2 in lymphoma are not mutually exclusive (Morin et al., 2011), suggesting these mutations may cooperate to deregulate bivalent domains in GC B-cells, or that they also have independent roles in lymphomagenesis.

EZH2 mutants may also mediate their actions through additional mechanisms. We and others have noted that H3K27 ChIP-seq profiles as well as gene expression profiles induced by EZH2 inhibitors are variable between cell lines (McCabe et al., 2012b), suggesting that hyperactive EZH2 molecule might lead to epigenetic instability and stochastic aberrant epigenetic silencing of different gene sets. On the other hand variability may be due to the fact that cell lines are genetically quite diverse and under the stress of continuous passage in vitro tend to drift apart epigenetically. This has been demonstrated in the case of cytosine methylation profiles, which differ considerably between lymphoma cell lines and primary lymphoma specimens (De et al., 2013). By contrast primary human DLBCL specimens with mutant EZH2 showed a robust signature consisting of greater repression of EZH2 target genes including bivalent domains and with similarity to the EZH2 mutant signature induced in the isogenic BCL1 experimental model. Moreover analysis of individual DLBCL cell lines after exposure to GSK343 showed upregulation of GC B-cell EZH2 targets including bivalent genes as well as phenotypic effects consistent with derepression of GC B-cell EZH2 target genes. Notably these studies underline that EZH2 mediated epigenetic effects are reversible in lymphomas, which is consistent with data in the prostate cancer field also showing that suppression of EZH2 can result in reactivation of genes with tumor-suppressing activity (Taniguchi et al., 2012; Yu et al., 2007).

Expression of mutant EZH2 alone in GC B-cells was insufficient to induce development of DLBCL. In this way EZH2 mutation appears analogous to many of the somatic mutations in acute myeloid leukemia (AML), which induce a myeloproliferative phenotype when expressed in murine hematopoietic stem cells but when expressed together cooperate to form AML (Shih et al., 2012). Mutant EZH2 induces a lymphoproliferative phenotype with expansion of the proliferative GC B-cell compartment. Constitutive expression of the GC B-cell oncoproteins BCL6 and BCL2 also manifest GC hyperplasia and a partially penetrant GC B-cell lymphoma phenotype (Cattoretti et al., 2005; Egle et al., 2004). Hence even though normal GC B-cells exhibit features of partially transformed cells such as suppression of proliferative checkpoints and attenuated DNA damage response, multiple oncogenic hits are still required for overt lymphomagenesis. This concept is supported by mutational profiling studies revealing multiple concurrent somatic mutations in DLBCL specimens (Lohr et al., 2012; Morin et al., 2011) as well as our transplantation studies demonstrating cooperation between mutant EZH2 and BCL2 in accelerating lymphomagenesis in mice. Knowledge of the genetic composition of lymphomas and how these cooperate to transform B-cells affords the opportunity to rationally design combinatorial therapies. The enhanced anti-lymphoma activities of GSK343 and GSK503 in combination with anti-BCL2 therapies support this notion and point towards design of clinical trials geared towards the underlying biology of DLBCL and a reduced reliance on relatively non-specific cytotoxic chemotherapy.

While previous reports have focused on the role of EZH2 inhibitors against mutant-EZH2 DLBCL (Knutson et al., 2012; McCabe et al., 2012b; Qi et al., 2012), our data indicate that EZH2 is a relevant target beyond those tumors. We find that GCB type DLBCLs are dependent on wild type EZH2 for their proliferation and survival regardless of somatic mutation, although response to EZH2 inhibitors is slightly delayed as compared to mutant EZH2 DLBCL cells. This finding is in agreement with the absolute requirement of normal GC B-cells for EZH2 and indicates that EZH2 is a lineage factor to which GC type DLBCLs are addicted (Figure 7F). By contrast ABC-DLBCLs do not require EZH2 to maintain their proliferation and survival. Indeed the target genes repressed by EZH2 consist of many of same genes that define the ABC-DLBCL subtype. These data provide a mechanism by which increased EZH2 mediated repression of target genes can impair B-cell differentiation, and demonstrate that therapeutic targeting of EZH2 in GCB DLBCLs can induce differentiation and abrogate proliferation of GCB-DLBCLs with mutant or wild-type EZH2. We thus provide the basis for the expanded clinical translation of EZH2 inhibitors for the treatment of GCB type DLBCLs. Clinical studies with pharmacologic EZH2 inhibitors will determine if this approach can improve outcomes for lymphoma patients. Although many GCB type DLBCLs can be cured with combination chemo-immunotherapy regimens, these involve the use of toxic drugs that carry a life long risk of developing second malignancies. The current standard therapy for B-cell lymphomas (RCHOP) wipes out the entire bone marrow of patients, and the B-cell lineage is completely eradicated for months by rituximab (anti-CD20 antibody). In comparison the transient suppression of high affinity antibody formation by EZH2 inhibitors seems less likely to be clinically significant. By targeting the oncoproteins that drive and define the GCB DLBCL phenotype it may be possible to reduce our reliance on cytotoxic drugs to eradicate this disease.

EXPERIMENTAL PROCEDURES

A more detailed description of the experimental procedures and reagents used in this study can be found in the Supplemental Experimental Procedures

Generation of conditional Ezh2Y641N knock-on mice

Ezh2Y641N mice were generated by Frt-mediated targeting of Ezh2Y641N into the ColA1 locus. The strategy for FLPe-mediated recombination (shown in Figure S2) was done as described previously (Beard et al., 2006). Briefly, homologous recombination was used to place an frt-hygro-pA cassette downstream of the ColA1 locus. Ezh2Y641N cDNA was then targeted to the modified locus by coelectroporation of a pgkATGfrt vector carrying Ezh2Y641N and an FLPe transient expression vector. This pgkATGfrt contains the PGK promoter, an ATG initiation codon, an frt site, and a CAG promoter driving the expression of Ezh2Y641N. Interspersed between the CAG promoter and Ezh2Y641N cDNA is a lox-stop-lox (LSL) cassette ensuring that Ezh2Y641N is only expressed in the presence of cre-mediated excision of the LSL cassette. Interstrand intrachromosomal recombination at the frt sites results in loss of the PGKneopA cassette and insertion of Ezh2Y641N and pgkATG cassette to restore and confer hygromycin resistance. The integration site is located ~500bp downstream of the 3′ untranslated region of the ColA1 locus. This left the coding portion and 3′ UTR transcriptional control region of the ColA1 gene intact. V6.5 ES cells were electroporated with this homologous targeting construct and stable G418 clones were derived.

Conditional transgenic murine models

The Research Animal Resource Center of the Weill Cornell Medical College of Medicine approved all mouse procedures. Conditional Ezh2 knockout mice (_loxP_-flanked EZH2 allele, EZH2fl/fl) were a generous gift of Dr. Alexander Tarakhovsky, The Rockefeller University (Su et al., 2003). By crossing _Ezh2_fl/fl with the transgenic Cγicre strain (The Jackson Laboratory, 010611) we generated heterozygous mice, which were crossed to yield _Ezh2_−/− mice. As control group, we used _Ezh2_fl/fl Cγicre negative littermates. _Ezh2_Y641N transgenic mice were crossed with Cγicre strain. The Cγicre negative littermates were used as control group. Conditional Ezh2 knockout and _Ezh2_Y641N knockin mice were used for assessment of the germinal center formation, which were induced with SRBC or NP-KLH.

DLBCL patient samples

Patient-deidentified left-over tissues were obtained at diagnosis from 69 patients with de novo DLBCL in Vancouver at the British Columbia Cancer Agency. Cases were selected on the basis of the presence of at least 80% of the neoplastic cells within the tumor section. The use of human tissue was approved by the research ethics board of the Vancouver Cancer Center/University of British Columbia and the Weill Cornell Medical Center. For additional information see Table S6 and Supplemental Experimental Procedures.

ChIP-seq and mRNA-seq library preparation and Illumina sequencing processing

ChIP-seq and RNA-seq libraries were prepared using the Illumina ChIP-Seq and TruSeq RNA sample kits, respectively, according to the manufacturer. Libraries were validated using the Agilent Technologies 2100 Bioanalyzer and Quant-iT™ dsDNA HS Assay (Life Technologies), and 8-10 pM sequenced on HiSeq2000 sequencer as follows: ChIP-seq, 1 × 50; mRNA-seq, 2 × 50. RNA-seq data was aligned to transcripts using TopHat (Trapnell et al., 2009), gene expression quantified in RPKM (Mortazavi et al., 2008) and differentially expressed genes identified. ChIP-seq experiments from human B cells, human cell lines and murine cells were aligned to the hg18, hg19 and mm9 genome, respectively using ELAND. H3K4me3 ChIP-seq reads were called into peaks using the ChIPseeqer framework (p<10−15 and fold-change threshold 2) (Giannopoulou and Elemento, 2011) and H3K27me3 and EZH2 ChIP-seq reads were quantified in 1kb bins genome-wide, identifying regions of enrichment as consecutive bins with read counts greater than one standard deviation of the genome-wide mean. Data analysis is described in Supplemental Experimental Procedures.

Supplementary Material

01

02

SIGNIFICANCE.

We show that EZH2 is required for the humoral immune response. EZH2 function is linked to de novo formation of bivalent chromatin domains at genes required for exit of B-cells from the GC. Mutant EZH2 alleles tip the balance of these domains towards increased H3K27 methylation leading to more potent and permanent silencing. This manifests as differentiation blockade, GC hyperplasia and accelerated lymphomagenesis in cooperation with BCL2. Wild type and mutant EZH2 functions are essential for GCB type DLBCLs but dispensable for ABC-DLBCLs. These data reveal the function of EZH2 and its mutants in mature B-cells and lymphomas and provide the basis for specific rational combinatorial targeted therapy for GCB type DLBCLs.

HIGHLIGHTS.

EZH2 is required for germinal center formation and immunoglobulin affinity maturation.

Conditional expression of an EZH2 mutant lymphoma allele drives GC hyperplasia in vivo.

Mutant EZH2 functions in part through aberrant repression of GC specific bivalent promoters.

EZH2 cooperates with BCL2 to initiate and maintain diffuse large B-cell lymphomas.

ACKNOWLEDGEMENTS

We thank Dr. Alexander Tarakhovsky (Rockefeller University) for sharing the EZH2 conditional knockout mouse strain. AM is supported by the Burroughs Wellcome Foundation and Chemotherapy Foundation. AM, RL and JDL are supported by the Samuel Waxman Cancer Research Foundation and a collaborative trans-network grant from the National Cancer Institute Physical Sciences in Oncology Center program (U54 CA143869 and CA143879). Supported by a Leukemia and Lymphoma Society Specialized Center of Research Excellence (JDL, AM, OE) and the T & C Schwartz Family Foundation (AM and JDL). OAW is an American Society of Hematology Basic Research Fellow and is supported by a grant from the NIH K08 Clinical Investigator Award (1K08CA160647-01). LCC is a Raymond and Beverly Sackler Scholar and Scholar of the American Society of Hematology. DWS is supported by a CIHR postdoctoral fellowship award and RDG is supported by New Frontiers in Cancer Terry Fox Program project grant (#019001). OE and AM are supported by NCI R01 CA104348 . OE is supported by the NSF CAREER grant and grants from the Starr Cancer Consortium. This work was enabled by the Beverly and Raymond Sackler Center for Physical and Biomedical Sciences as well as the Weill Cornell Epigenomics Core Facility. SKV, MTMC, HMO, GSVA, RGK, YL, CFMH, and CLC are employees of GlaxoSmithKline.

Footnotes

Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

ACCESSION NUMBER The Gene Expression Omnibus accession numbers for the ChIP-seq, RNA-seq and the expression microarray data reported in this paper are GSE45982 and GSE23501.

REFERENCES

  1. Alizadeh AA, Eisen MB, Davis RE, Ma C, Lossos IS, Rosenwald A, Boldrick JC, Sabet H, Tran T, Yu X, et al. Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature. 2000;403:503–511. doi: 10.1038/35000501. [DOI] [PubMed] [Google Scholar]
  2. Beard C, Hochedlinger K, Plath K, Wutz A, Jaenisch R. Efficient method to generate single-copy transgenic mice by site-specific integration in embryonic stem cells. Genesis. 2006;44:23–28. doi: 10.1002/gene.20180. [DOI] [PubMed] [Google Scholar]
  3. Bernstein BE, Mikkelsen TS, Xie X, Kamal M, Huebert DJ, Cuff J, Fry B, Meissner A, Wernig M, Plath K, et al. A bivalent chromatin structure marks key developmental genes in embryonic stem cells. Cell. 2006;125:315–326. doi: 10.1016/j.cell.2006.02.041. [DOI] [PubMed] [Google Scholar]
  4. Blackman MA, Tigges MA, Minie ME, Koshland ME. A model system for peptide hormone action in differentiation: interleukin 2 induces a B lymphoma to transcribe the J chain gene. Cell. 1986;47:609–617. doi: 10.1016/0092-8674(86)90625-2. [DOI] [PubMed] [Google Scholar]
  5. Cao R, Wang L, Wang H, Xia L, Erdjument-Bromage H, Tempst P, Jones RS, Zhang Y. Role of histone H3 lysine 27 methylation in Polycomb-group silencing. Science. 2002;298:1039–1043. doi: 10.1126/science.1076997. [DOI] [PubMed] [Google Scholar]
  6. Casola S, Cattoretti G, Uyttersprot N, Koralov SB, Seagal J, Hao Z, Waisman A, Egert A, Ghitza D, Rajewsky K. Tracking germinal center B cells expressing germ-line immunoglobulin gamma1 transcripts by conditional gene targeting. Proceedings of the National Academy of Sciences of the United States of America. 2006;103:7396–7401. doi: 10.1073/pnas.0602353103. [DOI] [PMC free article] [PubMed] [Google Scholar]
  7. Cattoretti G, Pasqualucci L, Ballon G, Tam W, Nandula SV, Shen Q, Mo T, Murty VV, Dalla-Favera R. Deregulated BCL6 expression recapitulates the pathogenesis of human diffuse large B cell lymphomas in mice. Cancer Cell. 2005;7:445–455. doi: 10.1016/j.ccr.2005.03.037. [DOI] [PubMed] [Google Scholar]
  8. Chase A, Cross NC. Aberrations of EZH2 in cancer. Clinical cancer research : an official journal of the American Association for Cancer Research. 2011;17:2613–2618. doi: 10.1158/1078-0432.CCR-10-2156. [DOI] [PubMed] [Google Scholar]
  9. Ci W, Polo JM, Cerchietti L, Shaknovich R, Wang L, Yang SN, Ye K, Farinha P, Horsman DE, Gascoyne RD, et al. The BCL6 transcriptional program features repression of multiple oncogenes in primary B cells and is deregulated in DLBCL. Blood. 2009;113:5536–5548. doi: 10.1182/blood-2008-12-193037. [DOI] [PMC free article] [PubMed] [Google Scholar]
  10. Ci W, Polo JM, Melnick A. B-cell lymphoma 6 and the molecular pathogenesis of diffuse large B-cell lymphoma. Curr Opin Hematol. 2008;15:381–390. doi: 10.1097/MOH.0b013e328302c7df. [DOI] [PMC free article] [PubMed] [Google Scholar]
  11. Czermin B, Melfi R, McCabe D, Seitz V, Imhof A, Pirrotta V. Drosophila enhancer of Zeste/ESC complexes have a histone H3 methyltransferase activity that marks chromosomal Polycomb sites. Cell. 2002;111:185–196. doi: 10.1016/s0092-8674(02)00975-3. [DOI] [PubMed] [Google Scholar]
  12. De S, Shaknovich R, Riester M, Elemento O, Geng H, Kormaksson M, Jiang Y, Woolcock B, Johnson N, Polo JM, et al. Aberration in DNA Methylation in B-Cell Lymphomas Has a Complex Origin and Increases with Disease Severity. PLoS Genetics. 2013;9:e1003137. doi: 10.1371/journal.pgen.1003137. [DOI] [PMC free article] [PubMed] [Google Scholar]
  13. Egle A, Harris AW, Bath ML, O’Reilly L, Cory S. VavP-Bcl2 transgenic mice develop follicular lymphoma preceded by germinal center hyperplasia. Blood. 2004;103:2276–2283. doi: 10.1182/blood-2003-07-2469. [DOI] [PubMed] [Google Scholar]
  14. Giannopoulou EG, Elemento O. An integrated ChIP-seq analysis platform with customizable workflows. BMC bioinformatics. 2011;12:277. doi: 10.1186/1471-2105-12-277. [DOI] [PMC free article] [PubMed] [Google Scholar]
  15. Klein U, Dalla-Favera R. Germinal centres: role in B-cell physiology and malignancy. Nat Rev Immunol. 2008;8:22–33. doi: 10.1038/nri2217. [DOI] [PubMed] [Google Scholar]
  16. Knutson SK, Wigle TJ, Warholic NM, Sneeringer CJ, Allain CJ, Klaus CR, Sacks JD, Raimondi A, Majer CR, Song J, et al. A selective inhibitor of EZH2 blocks H3K27 methylation and kills mutant lymphoma cells. Nature chemical biology. 2012;8:890–896. doi: 10.1038/nchembio.1084. [DOI] [PubMed] [Google Scholar]
  17. Lohr JG, Stojanov P, Lawrence MS, Auclair D, Chapuy B, Sougnez C, Cruz-Gordillo P, Knoechel B, Asmann YW, Slager SL, et al. Discovery and prioritization of somatic mutations in diffuse large B-cell lymphoma (DLBCL) by whole-exome sequencing. Proceedings of the National Academy of Sciences of the United States of America. 2012;109:3879–3884. doi: 10.1073/pnas.1121343109. [DOI] [PMC free article] [PubMed] [Google Scholar]
  18. Mahmoudi T, Verrijzer CP. Chromatin silencing and activation by Polycomb and trithorax group proteins. Oncogene. 2001;20:3055–3066. doi: 10.1038/sj.onc.1204330. [DOI] [PubMed] [Google Scholar]
  19. McCabe MT, Graves AP, Ganji G, Diaz E, Halsey WS, Jiang Y, Smitheman KN, Ott HM, Pappalardi MB, Allen KE, et al. Mutation of A677 in histone methyltransferase EZH2 in human B-cell lymphoma promotes hypertrimethylation of histone H3 on lysine 27 (H3K27) Proceedings of the National Academy of Sciences of the United States of America. 2012a;109:2989–2994. doi: 10.1073/pnas.1116418109. [DOI] [PMC free article] [PubMed] [Google Scholar]
  20. McCabe MT, Ott HM, Ganji G, Korenchuk S, Thompson C, Van Aller GS, Liu Y, Graves AP, Della Pietra A, 3rd, Diaz E, et al. EZH2 inhibition as a therapeutic strategy for lymphoma with EZH2-activating mutations. Nature. 2012b;492:108–112. doi: 10.1038/nature11606. [DOI] [PubMed] [Google Scholar]
  21. Morin RD, Johnson NA, Severson TM, Mungall AJ, An J, Goya R, Paul JE, Boyle M, Woolcock BW, Kuchenbauer F, et al. Somatic mutations altering EZH2 (Tyr641) in follicular and diffuse large B-cell lymphomas of germinal-center origin. Nature genetics. 2010;42:181–185. doi: 10.1038/ng.518. [DOI] [PMC free article] [PubMed] [Google Scholar]
  22. Morin RD, Mendez-Lago M, Mungall AJ, Goya R, Mungall KL, Corbett RD, Johnson NA, Severson TM, Chiu R, Field M, et al. Frequent mutation of histone-modifying genes in non-Hodgkin lymphoma. Nature. 2011;476:298–303. doi: 10.1038/nature10351. [DOI] [PMC free article] [PubMed] [Google Scholar]
  23. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature methods. 2008;5:621–628. doi: 10.1038/nmeth.1226. [DOI] [PubMed] [Google Scholar]
  24. Muller J, Hart CM, Francis NJ, Vargas ML, Sengupta A, Wild B, Miller EL, O’Connor MB, Kingston RE, Simon JA. Histone methyltransferase activity of a Drosophila Polycomb group repressor complex. Cell. 2002;111:197–208. doi: 10.1016/s0092-8674(02)00976-5. [DOI] [PubMed] [Google Scholar]
  25. O’Carroll D, Erhardt S, Pagani M, Barton SC, Surani MA, Jenuwein T. The polycomb-group gene Ezh2 is required for early mouse development. Molecular and cellular biology. 2001;21:4330–4336. doi: 10.1128/MCB.21.13.4330-4336.2001. [DOI] [PMC free article] [PubMed] [Google Scholar]
  26. Qi W, Chan H, Teng L, Li L, Chuai S, Zhang R, Zeng J, Li M, Fan H, Lin Y, et al. Selective inhibition of Ezh2 by a small molecule inhibitor blocks tumor cells proliferation. Proceedings of the National Academy of Sciences of the United States of America. 2012;109:21360–21365. doi: 10.1073/pnas.1210371110. [DOI] [PMC free article] [PubMed] [Google Scholar]
  27. Raaphorst FM, van Kemenade FJ, Fieret E, Hamer KM, Satijn DP, Otte AP, Meijer CJ. Cutting edge: polycomb gene expression patterns reflect distinct B cell differentiation stages in human germinal centers. J Immunol. 2000;164:1–4. doi: 10.4049/jimmunol.164.1.1. [DOI] [PubMed] [Google Scholar]
  28. Rui L, Schmitz R, Ceribelli M, Staudt LM. Malignant pirates of the immune system. Nature immunology. 2011;12:933–940. doi: 10.1038/ni.2094. [DOI] [PubMed] [Google Scholar]
  29. Shih AH, Abdel-Wahab O, Patel JP, Levine RL. The role of mutations in epigenetic regulators in myeloid malignancies. Nature reviews Cancer. 2012;12:599–612. doi: 10.1038/nrc3343. [DOI] [PubMed] [Google Scholar]
  30. Shin DM, Liu R, Wu W, Waigel SJ, Zacharias W, Ratajczak MZ, Kucia M. Global gene expression analysis of very small embryonic-like stem cells reveals that the Ezh2-dependent bivalent domain mechanism contributes to their pluripotent state. Stem cells and development. 2012;21:1639–1652. doi: 10.1089/scd.2011.0389. [DOI] [PMC free article] [PubMed] [Google Scholar]
  31. Sneeringer CJ, Scott MP, Kuntz KW, Knutson SK, Pollock RM, Richon VM, Copeland RA. Coordinated activities of wild-type plus mutant EZH2 drive tumor-associated hypertrimethylation of lysine 27 on histone H3 (H3K27) in human B-cell lymphomas. Proceedings of the National Academy of Sciences of the United States of America. 2010;107:20980–20985. doi: 10.1073/pnas.1012525107. [DOI] [PMC free article] [PubMed] [Google Scholar]
  32. Su IH, Basavaraj A, Krutchinsky AN, Hobert O, Ullrich A, Chait BT, Tarakhovsky A. Ezh2 controls B cell development through histone H3 methylation and Igh rearrangement. Nature immunology. 2003;4:124–131. doi: 10.1038/ni876. [DOI] [PubMed] [Google Scholar]
  33. Taniguchi H, Jacinto FV, Villanueva A, Fernandez AF, Yamamoto H, Carmona FJ, Puertas S, Marquez VE, Shinomura Y, Imai K, et al. Silencing of Kruppel-like factor 2 by the histone methyltransferase EZH2 in human cancer. Oncogene. 2012;31:1988–1994. doi: 10.1038/onc.2011.387. [DOI] [PMC free article] [PubMed] [Google Scholar]
  34. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–1111. doi: 10.1093/bioinformatics/btp120. [DOI] [PMC free article] [PubMed] [Google Scholar]
  35. van Galen JC, Dukers DF, Giroth C, Sewalt RG, Otte AP, Meijer CJ, Raaphorst FM. Distinct expression patterns of polycomb oncoproteins and their binding partners during the germinal center reaction. European journal of immunology. 2004;34:1870–1881. doi: 10.1002/eji.200424985. [DOI] [PubMed] [Google Scholar]
  36. van Kemenade FJ, Raaphorst FM, Blokzijl T, Fieret E, Hamer KM, Satijn DP, Otte AP, Meijer CJ. Coexpression of BMI-1 and EZH2 polycomb-group proteins is associated with cycling cells and degree of malignancy in B-cell non-Hodgkin lymphoma. Blood. 2001;97:3896–3901. doi: 10.1182/blood.v97.12.3896. [DOI] [PubMed] [Google Scholar]
  37. Velichutina I, Shaknovich R, Geng H, Johnson NA, Gascoyne RD, Melnick AM, Elemento O. EZH2-mediated epigenetic silencing in germinal center B cells contributes to proliferation and lymphomagenesis. Blood. 2010;116:5247–5255. doi: 10.1182/blood-2010-04-280149. [DOI] [PMC free article] [PubMed] [Google Scholar]
  38. Verma SK, Tian X, LaFrance LV, Duquenne C, Suarez DP, Newlander KA, Romeril SP, Burgess JL, Grant SW, Brackley JA, et al. Identification of Potent, Selective, Cell-Active Inhibitors of the Histone Lysine Methyltransferase EZH2. Acs Med Chem Lett. 2012;3:1091–1096. doi: 10.1021/ml3003346. [DOI] [PMC free article] [PubMed] [Google Scholar]
  39. Yap DB, Chu J, Berg T, Schapira M, Cheng SW, Moradian A, Morin RD, Mungall AJ, Meissner B, Boyle M, et al. Somatic mutations at EZH2 Y641 act dominantly through a mechanism of selectively altered PRC2 catalytic activity, to increase H3K27 trimethylation. Blood. 2011;117:2451–2459. doi: 10.1182/blood-2010-11-321208. [DOI] [PMC free article] [PubMed] [Google Scholar]
  40. Yu J, Cao Q, Mehra R, Laxman B, Yu J, Tomlins SA, Creighton CJ, Dhanasekaran SM, Shen R, Chen G, et al. Integrative genomics analysis reveals silencing of beta-adrenergic signaling by polycomb in prostate cancer. Cancer Cell. 2007;12:419–431. doi: 10.1016/j.ccr.2007.10.016. [DOI] [PubMed] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

01

02