E2A is an essential regulator of early B cell development. Here, we have demonstrated that E2A together with E2-2 controlled germinal center (GC) B cell and plasma cell development. As shown by the identification of regulated E2A,E2-2 target genes in activated B cells, these E-proteins directly activated genes with important functions in GC B cells and plasma cells by inducing and maintaining DNase I hypersensitive sites. Through binding to multiple enhancers in the Igh 3′ regulatory region and Aicda locus, E-proteins regulated class switch recombination by inducing both Igh germline transcription and AID expression. By regulating 3′ Igk and Igh enhancers and a distal element at the Prdm1 (Blimp1) locus, E-proteins contributed to Igk, Igh, and Prdm1 activation in plasmablasts. Together, these data identified E2A and E2-2 as central regulators of B cell immunity.
B cell immunity provides acute and long-term protection of the host against infections through the generation and secretion of high-affinity antibodies that recognize a shear unlimited number of pathogens. This enormous adaptive potential of B cells is brought about by V(D)J recombination of the immunoglobulin heavy chain (Igh) and light chain (Igk and Igl) genes in early B cell development, and by subsequent affinity maturation of the Ig heavy chain in late B cell differentiation (Victora and Nussenzweig, 2012). Somatic hypermutation alters the antigen-binding VH sequences of the Ig heavy-chain, whereas class switch recombination (CSR) exchanges the CH exons to generate Ig isotypes with distinct effector functions (Chaudhuri and Alt, 2004). Whereas the activation-induced deaminase (AID) is an essential regulator of both processes (Muramatsu et al., 2000), somatic hypermutation only takes place in germinal centers (GCs), which are formed upon antigen exposure by the interplay of T follicular helper (Tfh) cells and follicular (FO) B cells in secondary lymphoid organs (Victora and Nussenzweig, 2012). Affinity-based selection in this specialized compartment leads to clonal expansion of B cells expressing high-affinity B cell receptors, which subsequently differentiate to proliferating, antibody-secreting plasmablasts (Victora and Nussenzweig, 2012). Upon migration to specialized bone marrow niches, plasmablasts differentiate into long-lived quiescent plasma cells secreting high amounts of antibodies (Nutt et al., 2015). While many transcription factors are involved in coordinating these B cell responses, we have here studied the role of E-proteins in the regulation of these processes.
Basic helix–loop–helix (bHLH) transcription factors can be subdivided into different classes based on biochemical and functional properties (Murre, 2005). Class I bHLH proteins, also known as E-proteins, consist of the three members, E2A (Tcf3), E2-2 (Tcf4), and HEB (Tcf12; Murre, 2005), which bind the E-box (CANNTG) motif with similar sequence specificity (Fig. S1 A). E-proteins are broadly expressed and heterodimerize with class II bHLH proteins in nonlymphoid cell types. Within the lymphoid system, E-proteins function as homodimers or heterodimers with a different E-protein (Bain et al., 1993; Shen and Kadesch, 1995). E-proteins are thought to mainly function as transcriptional activators, as they interact with the co-activators p300 and CBP (Bradney et al., 2003; Bayly et al., 2004), as well as the promoter recognition factor TFIID (Chen et al., 2013). The activity of E-proteins is controlled by the inhibitor of DNA binding proteins, which are HLH proteins lacking the basic DNA-binding domain, and are thus capable of sequestering E-proteins into DNA-binding–incompetent heterodimers (Kee, 2009).
E-proteins control different aspects of B cell development (Murre, 2005). E2A is required for the commitment of lymphoid progenitors to the B cell lineage (Bain et al., 1994; Zhuang et al., 1994). E2A is also essential for Vκ-Jκ recombination (Inlay et al., 2004) and early B cell development (Kwon et al., 2008). Inactivation of all E-proteins in activated B cells by overexpression of the antagonist ID3 revealed an important role for these transcription factors in promoting CSR to different IgG isotypes (Quong et al., 1999) and activating the Aicda (AID) gene (Sayegh et al., 2003). As shown by conditional Tcf3 inactivation, E2A is largely dispensable for the formation and function of different mature B cell types and plasma cells, except for GC B cell differentiation, which is reduced but not lost in the absence of E2A (Kwon et al., 2008). It is, however, possible that the activity of another E-protein may compensate for the loss of E2A in late B cell differentiation in analogy to the cooperative function of E2A and HEB in T cell development (Jones-Mason et al., 2012).
Here, we have used conditional mutagenesis to demonstrate a cooperative role of E2A and E2-2 in controlling GC B cell and plasma cell development. Using genome-wide approaches, we comprehensively analyzed the molecular role of E2A and E2-2 in late B cell development, which revealed that these E-proteins directly control many essential functions of GC B cells and plasma cells. Hence, these experiments identified E2A and E2-2 as central regulators of B cell immunity.
Efficient generation of mature B cells upon combined loss of E2A and E2-2
As shown by RNA-seq, Tcf3 (E2A) was highly expressed in FO and GC B cells compared with Tcf4 (E2-2) and Tcf12 (HEB; Fig. 1 A). Tcf4 was, however, similarly expressed like Tcf3 in bone marrow plasma cells in contrast to Tcf12. We therefore hypothesized that Tcf4 likely compensates for the loss of E2A in late B cell development. To test this hypothesis, we used the Cd23-Cre line, which initiates Cre-mediated deletion in immature B cells of the spleen (Kwon et al., 2008), the floxed Tcf3 allele (Tcf3fl), which expresses a nonfunctional E2A-GFP fusion protein upon Cre-mediated elimination of the bHLH domain-encoding exons (Kwon et al., 2008), and the Tcf4fl allele (Bergqvist et al., 2000). We thus generated Tcf3fl/fl Tcf4fl/fl mice (referred to as Tcf3,4fl/fl or ‘WT’ mice) and Cd23-Cre Tcf3fl/fl Tcf4fl/fl mice (referred to as Cd23-Cre Tcf3,4fl/fl for the mice and DKO for the respective B cells). As shown by flow cytometric analysis, mature B cells (B220+CD19+IgMloIgDhi), FO B cells (B220+CD19+CD21intCD23hi), and marginal zone (MZ) B cells (B220+CD19+CD21hiCD23lo/–) were present at similar or slightly reduced numbers in the spleen of Cd23-Cre Tcf3,4fl/fl mice compared with Tcf3,4fl/fl littermates (Fig. 1 B). GFP expression furthermore suggested complete Tcf3 deletion in FO and MZ B cells of Cd23-Cre Tcf3,4fl/fl mice (Fig. 1 B), which was confirmed by PCR genotyping and immunoblot analysis with an E2A antibody (Fig. 1, C and D). In contrast, Tcf3 and Tcf4 were only partially deleted in splenic and peritoneal B-1 cells (B220loCD19+) of the Cd23-Cre Tcf3,4fl/fl genotype (Fig. 1 C and not depicted). Hence, FO and MZ B cells were efficiently generated in the absence of E2A and E2-2.
Loss of GC B cell differentiation in the absence of E2A and E2-2
To study the role of E2A and E2-2 in GC B cell development, we immunized mice with 4-hydroxy-3-nitrophenylacetyl-conjugated keyhole limpet hemocyanin (NP-KLH). 7 d after immunization, GC B cells could be detected in the spleen of Cd23-Cre Tcf4fl/fl and control Tcf3,4fl/fl mice as Fas+GL7+CD19+B220+ cells by flow cytometry (Fig. 2 A) and as GL7+ cells on histological sections (Fig. 2 B). As previously shown (Kwon et al., 2008), the GC B cell number and GC size were strongly reduced in Cd23-Cre Tcf3fl/fl mice (Fig. 2, A and B), consistent with a prominent role of E2A in GC B cell development. GC B cells were, however, completely absent in Cd23-Cre Tcf3,4fl/fl mice (Fig. 2, A and B). These observations were confirmed by analyzing splenic GC B cells at day 14 and lymph node GC B cells at day 7 after NP-KLH immunization (Fig. 2, C and D), as well as by investigating splenic GC B cells at day 14 after immunization with sheep RBCs (SRBCs; Fig. 2 E). Hence, the development of GC B cells critically depends on both E2A and E2-2.
Gene regulation by E2A and E2-2 in anti-CD40 and IL-4–activated B cells
We next determined the E2A,E2-2–dependent gene expression program in activated B cells that were stimulated with anti-CD40 antibodies and IL-4, which mimics the T cell help required for GC B cell differentiation. First, we determined the genome-wide pattern of E-protein binding by chromatin immunoprecipitation (ChIP)-seq of unstimulated FO B cells and activated B cells after stimulation with anti-CD40 and IL-4 for 2 d. ChIP was performed with an E2A antibody after cross-linking of the chromatin with formaldehyde alone (single cross-linking [s]) or with a combination of disuccinimidyl glutarate and formaldehyde (double cross-linking [d]). E2A binding at the Pou2af1 and Mef2b loci (Fig. 3 A) and genome-wide E2A peak calling (Fig. S1, B and C) revealed a significant overlap of the peaks identified with both cross-linking methods. We therefore called an E2A peak only if it was detected by both single and double cross-linking in the same cell type. Based on this stringent peak calling, we determined 5,314 and 7,571 E2A peaks in FO (day 0) and activated B cells (day 2), respectively, which contained the consensus E2A-binding motif (Figs. 3 B and S1, D and E). Peak-to-gene assignment defined 3,893 and 4,946 E2A-bound target genes in FO and activated B cells, respectively (Fig. 3 B).
Second, we determined the transcriptome of FO B cells of the Tcf3,4fl/fl (WT) and Cd23-Cre Tcf3,4fl/fl (DKO) genotypes before or after anti-CD40 and IL-4 stimulation for 1, 2, or 3 d by RNA-seq. Gene expression comparison of activated WT and DKO B cells at day 3 identified 157 activated and 108 repressed genes, which were selected for an expression difference of greater than threefold, an adjusted p-value of <0.05, and an RPKM value of >3 in stimulated WT (activated) or DKO (repressed) B cells, respectively (Fig. S1 F). By determining the overlap between the E2A-bound genes (Fig. 3 B) and E2A,E2-2-regulated genes (Fig. S1 F), we identified 120 potentially directly activated and 32 potentially directly repressed E2A,E2-2 target genes (Fig. 3 C and Table S1). Notably, the vast majority (76.4%) of all activated genes (157) were also bound by E2A, in contrast to only 29.6% of all repressed genes (108; Fig. 3 C). Hence, E2A and E2-2 primarily activate gene transcription, whereas they repress genes mainly in an indirect manner in activated B cells.
More than half of all activated E-protein target genes code for cell surface receptors (20), signal transducers (26), and transcriptional regulators (22), which suggests a role for E2A and E2-2 in B cell signaling (Fig. 3 D). Importantly, four activated target genes (Icosl, Mef2b, Pou2af1, and Neil1) are known to play important roles in GC B cell differentiation (Fig. 3 E). E2A and E2-2 also activated three transcription factor genes (Prdm1 [Blimp1], Xbp1, and Eaf2) with important functions in plasma cell differentiation (Fig. 3 F), as well as four genes (Mzb1, Edem1, Fcrla, and Wfs1) that contribute to the homeostatic control of the ER (Fig. S2 A). The E2A,E2-2–dependent regulation of additional transcription factor genes with known functions in B cells (Bhlhe41, Id3, Bhlha15, Klf2, Arid3a, and Bcl3) and T cells (Aire, Ahr, Bcl6b, and Hivep3) are shown in Fig. S2 B. Another functional class of eight activated target genes codes for potentially inhibitory (Sit1, Lax1, Rasal1, Ptpn3, and Dusp6) or stimulatory (Lat, Src, and Ralgds) signal transducers of BCR signaling (Figs. 3 G and S2 C). Finally, we identified 23 activated target genes coding for cell surface proteins, signal transducers, and cytoskeletal proteins involved in cell adhesion and migration (Figs. 3 H and S2 D), indicating that E2A and E2-2 may control the migratory or sessile behavior of activated B cells.
E-proteins induce and maintain open chromatin at activated target genes
To investigate a possible role of E-proteins in controlling regulatory elements at activated target genes, we mapped open chromatin regions (known as DNase I hypersensitive [DHS] sites) by assay for transposase-accessible chromatin (ATAC)-seq (Buenrostro et al., 2013) in activated WT and DKO B cells at day 2 of anti-CD40 and IL-4 stimulation. The chromatin accessibility, which was measured as read density of each DHS site, was increased at E2A peaks of activated genes in WT B cells compared with DKO B cells (Fig. 4 A, left). In contrast, the densities of DHS sites lacking E2A binding at activated genes were similar to those of all DHS sites in both cell types (Fig. 4 A, right). This suggests that E-proteins are essential for inducing or maintaining DHS sites at activated target genes. As exemplified for the Lmo7 gene, the induction of a DHS site correlated with increased E2A binding upon stimulation of WT B cells, but not of DKO B cells (Fig. 4 B). Conversely, an E2A-bound DHS site, which was detected at the Selplg gene in FO and activated B cells, was specifically lost in activated DKO B cells (Fig. 4 B). We conclude therefore that E-proteins are required for the induction and maintenance of DHS sites at activated target genes.
E2A and E2-2 control CSR by regulating
We next investigated CSR to IgG1 in GC B cells (Fas+GL7+), which were still formed in Peyer’s patches in the absence of E2A and E2-2 (Fig. 5 A). IgG1+Fas+ GC B cells were absent in Peyer’s patches of Cd23-Cre Tcf3,4fl/fl mice (Fig. 5 A). Moreover, IgG1+Fas+ B cells were more strongly reduced in Cd23-Cre Tcf3fl/fl mice compared with Cd23-Cre Tcf4fl/fl mice, indicating that E2A is the dominant E-protein controlling IgG1 CSR in Peyer’s patches (Fig. 5 A). A similar situation was observed upon in vitro stimulation of FO B cells from WT and mutant lymph nodes for 4 d with anti-CD40 and IL-4 (Fig. 5 B). CSR to IgG1 was lost in DKO B cells, although the proliferation of these cells was only minimally reduced compared with WT B cells, as shown by dilution of the CellTrace Violet reagent (Fig. 5 C). Notably, treatment of FO B cells with LPS and IL-4 for 4 d revealed an equally important role of E2A and E2-2 in controlling CSR to IgG1 (Fig. 5 D) in contrast to the dominant role of E2A observed upon anti-CD40 plus IL-4 stimulation (Fig. 5 B). Hence, both E2A and E2-2 are essential for IgG1 CSR.
The switch regions of Igh constant genes are made accessible for CSR by germline transcription from an upstream I promoter that is activated upon signaling by specific cytokines (Chaudhuri and Alt, 2004). Stimulation of WT B cells with anti-CD40 and IL-4 for 3 d strongly induced the expression of Iγ1 and Iε germline transcripts (GLTs) in addition to the constitutively expressed Iμ transcript (Fig. 5 E). Importantly, the Iμ and Iε GLTs were 2.3- and 7-fold reduced in activated DKO B cells (Fig. 5 E). Unexpectedly, the Iγ1 GLT was only minimally decreased upon loss of E2A and E2-2, indicating that the absence of IgG1 CSR in DKO B cells is caused by another defect (Fig. 5 E). E2A binding was detected at the Iγ1 promoter and a downstream enhancer (Fig. 5 F), as well as at the Iε promoter (Fig. 5 G) in activated B cells. At both genes, DHS sites remained unaffected, and active (H3K27ac) chromatin was minimally reduced in activated DKO B cells, suggesting that the decrease of Iε GLTs is caused by a defect at a distant regulatory element. Indeed, the abundance of H3K27ac was strongly reduced, and the DHS sites were lost at three enhancers (HS3A, HS1,2, and HS3B) of the Igh 3′ regulatory region (3′RR) in activated DKO B cells (Fig. 5 H). Notably, E2A binding was observed at all three enhancers (Fig. 5 H). As the 3′RR is important for CSR to all Igh isotypes (Vincent-Fabert et al., 2010; Saintamand et al., 2015), we conclude that E2A and E2-2 regulate Igh germline transcription, and thus CSR, by controlling the activity of the 3′RR.
E-proteins regulate the
Aicda locus by binding to multiple enhancers
Retroviral overexpression of the antagonist ID3 in activated B cells demonstrated a role for E-proteins in controlling Aicda expression (Sayegh et al., 2003). We confirmed this finding by providing genetic evidence that the loss of E2A and E2-2 prevented Aicda activation in response to anti-CD40 and IL-4 stimulation (Fig. 6 A). The Aicda locus contains five enhancer regions (E1-E5), which interact with each other and the promoter (P1) to activate Aicda expression (Kieffer-Kwon et al., 2013). We detected E2A binding not only at the previously identified site in the intronic enhancer E4 (Sayegh et al., 2003), but also at the upstream enhancers E1 and E2 and downstream enhancer E5 in activated B cells at day 2 of anti-CD40 and IL-4 treatment (Fig. 6 B). Active chromatin (H3K27ac) was strongly reduced at all E2A-bound enhancers (E1, E2, E4, and E5) in activated DKO B cells relative to WT B cells, but the abundance of H3K27ac was similar in both cell types at the enhancer E3 (Fig. 6 B) and in a genomic control region (Fig. S1 G). The density of the DHS sites at all four E2A-bound enhancers was also decreased in activated DKO B cells (Fig. 6 B). Collectively, these data indicate that E-proteins activate Aicda transcription by controlling multiple enhancers.
As E-proteins affect CSR by regulating Aicda and Igh GLTs, we discriminated between these two functions by performing retroviral rescue experiments. FO B cells from Cd23-Cre Tcf3,4fl/fl or Aicda−/− mice were infected with an empty retrovirus (MiCD2) expressing the human CD2 indicator protein or with retroviruses additionally expressing E2A (MiCD2-E2A), E2-2 (MiCD2-E2-2), or AID (MiCD2-AID). Flow cytometric analysis of the infected hCD2+ B cells at day 4 after LPS and IL-4 stimulation revealed that CSR to IgG1 and IgE was rescued by expression of E2A and E2-2 in DKO B cells, but not in Aicda−/− B cells, as expected (Fig. 6 C). Interestingly, AID expression restored CSR to IgG1 in activated DKO B cells, indicating that the loss of AID expression is the main reason for the absence of IgG1 CSR in these mutant B cells (Fig. 6 C). In contrast, AID expression could not rescue CSR to IgE in DKO cells (Fig. 6 C), as this switching process additionally depends on the E-protein–dependent expression of Iε GLTs (Fig. 5 E). Hence, these experiments further demonstrated that E2A and E2-2 control CSR by activating both Aicda expression and Igh germline transcription.
To investigate whether E-proteins are equally important for Aicda activation by different stimuli, we treated WT and DKO B cells for 3 d with different stimulation conditions before RT-qPCR analysis of Aicda mRNA. Aicda induction in response to anti-CD40, IL-4, and IL-5 or anti-CD40 and IL-4 stimulation required the presence of E2A and E2-2 (Fig. 6 D), consistent with the RNA-seq data of Fig. 6 A. Unexpectedly, Aicda expression was equally well induced in DKO and WT B cells upon stimulation with anti-CD40 and IL-21 (Fig. 6 D). We next cultured DKO B cells on stromal 40LB cells (Nojima et al., 2011), expressing BAFF and the membrane-bound CD40 ligand, in the presence of IL-4 for 4 d, which failed to activate CSR to IgG1 (Fig. 6 E), consistent with the data of Fig. 5 B. In contrast, DKO B cells underwent efficient IgG1 CSR upon further stimulation with IL-21 (replacing IL-4) for another 4 d (Fig. 6 E). We conclude therefore that Aicda activation differs in the requirement for E2A and E2-2, depending on the stimulation conditions.
Loss of plasma cells in the absence of E2A and E2-2
Long-lived memory plasma cells (CD138+CD28+B220intLin–; Delogu et al., 2006) were present at similar numbers in the bone marrow of nonimmunized Cd23-Cre Tcf3fl/fl (E2A KO), Cd23-Cre Tcf3,4fl/fl and control Tcf3,4fl/fl mice (Fig. 7 A). However, plasma cells were threefold reduced in the bone marrow of Cd23-Cre Tcf4fl/fl (E2-2 KO) mice (Fig. 7 A), indicating that E2-2 is the dominant E-protein consistent with its high expression in plasma cells (Fig. 1 A). Only half of all plasma cells in DKO mice expressed GFP (Fig. 7 A). Moreover, the second Tcf3 allele was incompletely deleted and the intact floxed Tcf4 alleles were retained in GFP+ DKO plasma cells, as shown by PCR genotyping (Fig. 7 B). Hence, the strong counterselection against Tcf3 and Tcf4 deletion indicates an essential role for E2A and E2-2 in plasma cell development. Immunization with the T cell–independent antigen TNP-Ficoll resulted at day 14 in a significant reduction of plasma cells in the spleen of both E2A KO and E2-2 KO mice. Although the DKO plasma cells were present in higher numbers, they largely failed to express GFP, further indicating efficient counterselection against the combined loss of E2A and E2-2 in plasma cells (Fig. 7 C). A similar picture was observed at day 14 after immunization with the T cell-dependent antigen NP-KLH (Fig. 7 D), which also demonstrated incomplete Tcf3 deletion and retention of the intact floxed Tcf4 alleles in GFP+ DKO plasma cells (Fig. 7 E). Hence, E2A and E2-2 are essential regulators of plasma cell development.
Loss of E2A and E2-2 arrests plasmablast differentiation at an activated B cell stage
We next investigated the role of E2A and E2-2 in an in vitro plasmablast differentiation system. For this, FO B cells from WT mice were stimulated with LPS for 4 d, which generated activated B cells (CD22+CD138–), preplasmablasts (CD22–CD138–), and plasmablasts (CD22–CD138+; Minnich et al., 2016; Fig. 8 A). In contrast, LPS-induced differentiation of DKO B cells resulted in a major CD22+CD138– cell population with high GFP expression and in a small population (0.6%) of DKO plasmablasts, which was largely GFP– and thus failed to delete Tcf3 (Fig. 8 A). Notably, the DKO B cells exhibited only a minor decrease in proliferation (Fig. 8 B), but contained 10-fold less anti-IgM–secreting cells compared with WT cells (Fig. 8 C). Collectively, these results demonstrate that the loss of E2A and E2-2 stringently arrests plasmablast differentiation at a CD22+CD138– cell stage.
We next performed RNA-seq with activated WT and DKO B cells at day 2 of LPS stimulation (Fig. S3 A), and with activated WT B cells, preplasmablasts, plasmablasts, and the arrested DKO cells at day 3 of LPS stimulation (Fig. S3 B). Principal component analysis of these RNA-seq data demonstrated that the DKO cells were most closely related to the activated WT B cells (Fig. 8 D). Moreover, the B cell–specific genes Pax5, Spib, Bcl11a, Cd22, and Cd40, which are normally repressed during plasmablast differentiation, were still highly expressed in the DKO cells (Fig. S4 A). Together, these findings indicate that the loss of E2A and E2-2 blocks LPS-induced differentiation at an activated B cells stage before the onset of plasmablast formation.
Analysis of the RNA-seq data of activated WT and DKO cells resulted in 41 activated and 45 repressed genes in activated B cells at day 2, as well as in 127 activated and 106 repressed genes at day 3 of LPS stimulation (Fig. S3, C and D). ChIP-seq analysis further revealed 4,092 E2A-bound target genes in LPS-activated B cells at day 2, as well as 5,488 and 6,985 E2A target genes in LPS-differentiated preplasmablast and plasmablasts at day 4, respectively (Fig. S3, E–I). By determining the overlap between the E2A,E2-2-regulated genes (Fig. S3 D) and E2A-bound genes (Fig. S3 I), we identified 110 potentially directly activated and 42 potentially directly repressed E2A,E2-2 target genes in activated B cells at day 3 of LPS stimulation (Fig. 8 E and Table S2). Whereas both anti-CD40 plus IL-4 and LPS stimulation conditions defined a common set of 33 activated E2A,E2-2 target genes, only three repressed target genes were common to both treatments (Fig. 8 F). Moreover, 46% of all activated target genes code for 14 surface receptors, 17 signal transducers, 16 transcriptional regulators, and 4 proteins involved in ER function (Fig. 8 G). We further divided the activated target genes of these four classes according to their expression during plasmablast differentiation. Activated target genes, which were not at all or only weakly up-regulated during plasmablast differentiation, code for 12 cell surface receptors (Tbxa2r, Slamf7, Plxnd1, Il9r, Cd2, Sdc4, Il6ra, Crim1, Cxcr4, Sell [CD62L], Gpr183 [EBI2], and Cd9), 10 signal transducers (Sit1, Rap1gap2, Pik3r5, Sik1, Rasgrp2, Pim1, Ticam2, Spred2, Map3k8, and Ralgds), and 8 transcriptional regulators (Id3, Cbfa2t3, Trp73, Id2, Mef2b, Bhlhe41, Hivep3, and Sox4; Fig. S4, B–D). The second class of activated target genes was more than threefold up-regulated during plasmablast differentiation and codes for two cell surface receptors (Ccr9 and Itgb3), seven signal transducers (Gnaz, Irs2, Plcd3, Lax1, Dusp5, Dusp14, and Myzap), eight transcriptional regulators (Blimp1 [Prdm1], Xbp1, Eaf2, Bhlha15, Hes1, Cebpb, Creb3l2, and Cbx4), and four molecules involved in protein secretion and homeostatic control of the ER (Tram2, Edem3, Edem1, Wfs1; Fig. 8, H–K). It is important to note that the strong induction of these genes in activated B cells at day 3 of LPS stimulation cannot be explained by contaminating preplasmablasts, as nonregulated genes with high expression in preplasmablasts were similarly expressed in activated WT and DKO B cells (Fig. S4 E). Hence, the identified activated target genes likely contribute to the stringent arrest of plasmablast differentiation in the absence of E2A and E2-2.
E-protein–dependent control of 3′ enhancers at the
Igh and Igk loci
We next studied the E-protein–dependent regulation of Igh and Igk gene expression in view of the fact that increased expression and secretion of immunoglobulins is a hallmark of plasma cells (Nutt et al., 2015). Whereas the expression of Igh and Igk transcripts was strongly increased during LPS-induced differentiation of activated B cells to plasmablasts, this increase was not observed in activated DKO cells (Fig. S5, A and B). Moreover, the mRNAs encoding the secreted Igμs, Igγ2bs, and Igγ3s proteins were strongly increased in WT preplasmablasts and plasmablasts in contrast to activated DKO B cells, indicating that the posttranscriptional switch to the Igμs, Igγ2bs, and Igγ3s transcripts did not take place in the absence of E2A and E2-2 (Fig. S5, B and C).
High expression of Ig heavy-chain proteins in plasma cells depends on the four enhancers (HS1/2, HS3A, HS3B, and HS4) in the 3′RR of the Igh locus (Vincent-Fabert et al., 2010). As shown by DHS site analysis, all four enhancers were already present in open chromatin in LPS-activated WT B cells (Fig. S5 D). In contrast, the accessibility at the HS3A, HS1,2, and HS3B enhancers was strongly reduced in activated DKO B cells (Fig. S5 D). Notably, E2A bound to all four enhancers in the 3′RR (Fig. S5 D). Hence, these data suggest that E2A and E2-2 are directly responsible for inducing open chromatin at the 3′RR in LPS-activated B cells, similar to the situation observed with anti-CD40 and IL-4–stimulated B cells (Fig. 5 H). Likewise, the loss of E2A and E2-2 in activated DKO B cells led to reduced chromatin accessibility at three E2A-bound enhancers (iEκ, 3′Eκ, and Ed) in the 3′ region of the Igk locus (Fig. S5 E). Hence, E-proteins regulate the activity of 3′ Igh and Igk enhancers in activated B cells.
Regulation of the
Prdm1 gene by a distant E-protein–dependent enhancer
Both LPS and anti-CD40 plus IL-4 stimulation identified Xbp1 and Prdm1 (Blimp1) as activated E2A,E2-2 target genes, which have important functions in plasma cells (Figs. 3 F and 8 I). The Xbp1 locus, which codes for an essential regulator of immunoglobulin secretion (Reimold et al., 2001), was bound by E2A at several putative upstream enhancers, some of which lost their open chromatin in activated DKO cells (Fig. S5 F). As Blimp1 is an essential regulator of plasma cell development (Martins and Calame, 2008), we next investigated whether retroviral Blimp1 expression in DKO B cells could restore plasmablast formation. Lymph node B cells from Tcf3,4fl/fl, Cd23-Cre Tcf3,4fl/fl, and Cd23-Cre Prdm1Gfp/fl mice were infected with the control retrovirus MiCD2 or the retroviruses MiCD2-E2-2 or MiCD2-Blimp1, followed by flow cytometric analysis of the infected hCD2+ B cells at day 4 after LPS stimulation (Fig. 9 A). As expected, E2-2 expression rescued the differentiation of DKO B cells to plasmablasts (Fig. 9 A), which efficiently secreted IgM antibodies (Fig. 9 B). Although retroviral Blimp1 expression could restore plasmablast differentiation of Blimp1 KO B cells, it was unable to restore the differentiation of E2A,E2-2–deficient DKO B cells to IgM-secreting plasmablasts (Fig. 9, A and B). Hence, Blimp1 expression is not sufficient to overcome the developmental arrest of activated DKO B cells.
Mapping of DHS sites in LPS-activated B cells and plasmablasts identified several potential enhancers in the upstream region of the Prdm1 gene (Fig. 10 A). Whereas E2A bound weakly to some of these elements, the most prominent E2A-binding site was detected at a putative enhancer located 272-kb upstream of the Prdm1 gene (Fig. 10 A). Importantly, the DHS site at this putative enhancer (referred to as site H) was lost in activated DKO B cells (Fig. 10 A). To determine a potential interaction between the Prdm1 promoter and DHS site H, we analyzed the three-dimensional (3D) architecture of the Prdm1 locus by performing chromosome conformation capture (3C) experiments with LPS-stimulated WT plasmablasts and in vitro–cultured control pro–B cells. The 3C-qPCR analysis was performed with primers located in the reference HindIII fragment at the Prdm1 promoter and in HindIII fragments containing different upstream DHS sites (Fig. 10 B). Relative cross-linking frequencies that were high in plasmablasts and low in pro–B cells identified long-range interactions between the Prdm1 promoter and the DHS sites B (−145 kb), C (−170 kb), D (−230 kb), F (−250 kb), and H (−272 kb) in plasmablasts (Fig. 10 B). To determine a possible role of the E2A-binding site at DHS site H in long-range looping, we generated a mouse carrying a 3-bp mutation (Mut) in the consensus E2A-binding sequence of DHS site H by CRISPR/Cas9-mediated genome engineering (Fig. 10 C). E2A binding to the mutated DHS site H was reduced by sixfold but not abolished in Prdm1Mut/Mut plasmablasts (Fig. 10 D). Notably, the partial loss of E2A binding specifically reduced the long-range interaction of DHS site H with the Prdm1 promoter in Prdm1Mut/Mut plasmablasts (Fig. 10 B) and resulted in a significant decrease of nascent Prdm1 transcripts in Prdm1Mut/Mut preplasmablast and plasmablasts (Fig. 10 E). We therefore conclude that E-proteins contribute to the 3D architecture of the Prdm1 locus and regulate Prdm1 transcription by activating the distal enhancer H.
E2A is an essential regulator of early B cell development (Murre, 2005; Kwon et al., 2008). Here, we have demonstrated that E2A cooperates with E2-2 in regulating late B cell development and B cell immunity, as both transcription factors are strictly required for the development of GC B cells and plasma cells in response to immunization. E2A proved to be the dominant E-protein for GC B cell differentiation, and E2-2 was more important for plasma cell development, in agreement with their relative expression in the two cell types.
The main function of E2A and E2-2 in activated B cells appears to be transcriptional activation rather than repression, as indicated by the fact that E2A bound to most activated genes in contrast to only one third of the repressed genes. Moreover, few repressed target genes were highly regulated by the two E-proteins in contrast to a large proportion of the activated target genes. In support of transcriptional activation, E-proteins contain 3 activation domains (AD1, AD2, and AD3), which interact with the coactivators p300 and CBP (Bradney et al., 2003; Bayly et al., 2004) and the promoter recognition factor TFIID (Chen et al., 2013). At the chromatin level, we identified a novel role for E-proteins in shaping the enhancer landscape at its activated target genes, as E2A and E2-2 are responsible for the induction and maintenance of DHS sites containing E2A-binding sites. This novel function, which is best exemplified by the E-protein–dependent enhancers located at the Aicda locus and in the 3′ regions of the Igh and Igk loci, may be mediated by recruitment of the histone acetyltransferases p300 and CBP to E2A-binding sites at activated target genes.
Conditional E2A loss in mature B cells is known to impair the development of GC B cells (Kwon et al., 2008). Here, we demonstrate that the combined loss of E2A and E2-2 entirely prevents GC B cell development. The loss of GC B cells is likely explained by the reduced expression of the activated E2A,E2-2 target genes Icosl, Pou2af1, and Mef2b. GC formation strictly depends on the interplay between Tfh and B cells, which is mediated by interaction of the ICOS receptor on Tfh cells with the ICOS ligand (Icosl) on B cells (Nurieva et al., 2008). GC B cells also fail to develop in the absence of the transcriptional co-activator OBF-1/OCA-B (Pou2af1; Schubart et al., 1996; Qin et al., 1998). Finally, the transcription factor MEF2B (Mef2b) was shown to directly activate the Bcl6 gene, which itself codes for an essential regulator of GC B cell development (Ying et al., 2013).
Retroviral overexpression of ID3 in activated B cells implicated E-proteins in the control of CSR (Quong et al., 1999) and AID expression (Sayegh et al., 2003). Here, we have significantly extended these findings by demonstrating that E-proteins control CSR by regulating Igh germline transcription through the 3′RR domain and by controlling Aicda expression through multiple enhancers. E2A controls three of the four enhancers at the Igh 3′RR in activated B cells, as shown by the following evidence. First, E2A bound to all four enhancers. Second, the DHS sites were strongly reduced or lost at three enhancers (HS3A, HS1,2, and HS3B) in activated DKO B cells. Third, active chromatin (H3K27ac) was strongly decreased in the entire 3′RR domain, possibly as a result of inefficient recruitment of the histone acetyltransferases p300 and CBP in the absence of E2A and E2-2. The loss of 3′RR enhancer activity likely explains the strong decrease of Iε germline transcription in DKO B cells stimulated with anti-CD40 and IL-4. Surprisingly however, Iγ1 germline transcription was minimally influenced by the loss of both E-proteins, consistent with the observation that deletion of the entire 3′RR leads to a relatively small decrease of Iγ1 GLTs and a partial loss of IgG1 CSR (Vincent-Fabert et al., 2010; Saintamand et al., 2015). In agreement with this finding, retroviral restoration of AID expression in DKO B cells was sufficient to rescue CSR to IgG1, but not to IgE, demonstrating that IgE CSR requires the E-protein–dependent activation of both the 3′RR enhancers and Aicda gene.
The Aicda locus contains five distinct enhancer regions (E1-E5), which interact with each other and the Aicda promoter to form a local promoter-enhancer interactome that functions as one cooperative regulatory unit to induce Aicda expression (Kieffer-Kwon et al., 2013). In support of this idea, each enhancer is essential for efficient Aicda induction (Crouch et al., 2007; Huong et al., 2013; Kieffer-Kwon et al., 2013). Here we have demonstrated that E-proteins bind to four Aicda enhancers (E1, E2, E4, and E5) and contribute to Aicda regulation by establishing active chromatin and DHS sites at these enhancers.
The loss of plasma cells in Cd23-Cre Tcf3,4fl/fl mice identified a novel role for E-proteins in terminal B cell differentiation. The strict dependency of plasma cell development on E2A and E2-2 is best documented by the strong counterselection against Tcf3 and Tcf4 deletion in the residual plasma cells of Cd23-Cre Tcf3,4fl/fl mice. Moreover, the loss of E2A and E2-2 stringently arrests LPS-induced plasmablast differentiation at an activated B cell stage. Whereas plasmablasts and plasma cells are characterized by a high rate of immunoglobulin secretion (Nutt et al., 2015), the arrested DKO B cells fail to up-regulate Igh and Igk transcription and to undergo the posttranscriptional expression switch from the membrane-bound to secreted immunoglobulin heavy-chain. The inability to further activate Igh gene transcription is likely caused by the observed loss of open chromatin at the three Igh 3′ enhancers HS3A, HS1,2, and HS3B in activated DKO B cells. E2A was originally discovered as a DNA-binding protein interacting with the iEκ enhancer (Murre et al., 1989). Later, E2A was shown to regulate Vκ-Jκ rearrangement in pre–B cells (Inlay et al., 2004). Here, we have demonstrated an important role for E2A and E2-2 in the control of Igk gene transcription by inducing DHS sites at the iEκ, 3′Eκ, and Ed enhancers in LPS-activated B cells. This evidence strongly suggests that E-proteins contribute to the strong up-regulation of immunoglobulin expression in plasmablasts and plasma cells by establishing functional enhancers at the 3′ end of the Igh and Igk loci in activated B cells, consistent with the requirement of an intact 3′RR for promoting high secretion of Ig heavy-chain proteins in plasma cells (Vincent-Fabert et al., 2010).
Five transcription factors, IRF4, Aiolos, Ikaros, Blimp1, and XBP1, are currently known to play essential roles in plasma cell development, in addition to the E-proteins described here (Nutt et al., 2015). Of these transcription factor genes, only Xbp1 and Prdm1 (Blimp1) are directly activated by E2A and E2-2. As XBP1 is required for antibody secretion, but not for plasmablast differentiation (Taubenheim et al., 2012), its loss cannot account for the E-protein–dependent block of plasma cell development. In contrast, loss of Blimp1 stringently arrests plasmablast differentiation at an early stage (Shapiro-Shelef et al., 2003; Kallies et al., 2007). Although we identified E2A binding at multiple DHS sites in the 5′ region of the Prdm1 locus, the most prominent E2A-binding region was detected at DHS site H, located 272-kb upstream of the Prdm1 promoter. This distal DHS site is not only lost in activated DKO B cells, but also interacts with the Prdm1 promoter. Importantly, mutation of a consensus E2A-binding sequence in DHS site H reduced not only the long-range interaction with the promoter, but also decreased Prdm1 transcription, indicating that DHS site H functions as a distal enhancer to activate the Prdm1 gene in plasmablasts. Retroviral restoration of Blimp1 expression in DKO B cells was, however, not sufficient to rescue plasmablast differentiation in the absence of E2A and E2-2, indicating that the reduced expression of some of the other 109 activated E2A,E2-2 target genes in DKO B cells also contributes to the differentiation block. Hence, E-proteins regulate plasma cell development in a pleiotropic manner, possibly as a result of their important function in shaping the enhancer landscape at their target genes during terminal B cell differentiation.
MATERIALS AND METHODS
The Tcf3fl/fl mice (Kwon et al., 2008), Tcf4fl/fl mice (Bergqvist et al., 2000), Prdm1fl/fl mice (Ohinata et al., 2005), Prdm1Gfp/+ mice (Kallies et al., 2004), Aicda−/− mice (Muramatsu et al., 2000), and transgenic Cd23-Cre mice (Kwon et al., 2008) were maintained on the C57BL/6 genetic background. All animal experiments were performed according to valid project licenses, which were approved and regularly controlled by the Austrian Veterinary Authorities.
The E2A-binding site in DHS site H was mutated in the endogenous Prdm1 locus by CRISPR-Cas9–mediated genome editing (Yang et al., 2013). For this, Cas9 mRNA was co-injected with a specific sgRNA (linked to the scaffold tracrRNA) and a single-stranded repair oligonucleotide (171 nucleotides) into mouse zygotes (C57BL/6 x CBA), as previously described (Yang et al., 2013). The sgRNA, repair oligonucleotide, and PCR genotyping primers are shown in Table S3. The 1,056-bp PCR product amplified from the Prdm1Mut allele was cleaved with XbaI, yielding 586-bp and 470-bp fragments, in contrast to the PCR fragment of the WT allele.
The following monoclonal antibodies were used for flow cytometric analysis of mouse spleen, lymph node, and bone marrow: B220/CD45R (RA3-6B2), CD4 (GK1.5), CD8a (53–6.7), CD11b/Mac1 (M1/70), CD19 (1D3), CD21/CD35 (7G6), CD22 (Cy34.1), CD23 (B3B4), CD28 (37.51), CD38 (90), CD49b (DX5), CD95/Fas (Jo2), CD138 (281–2), F4/80 (CI:A3-1), GL7 (GL-7), IgD (11.26c), IgE (R35-72), IgG1 (A85-1), IgG2b (R12-3), IgG3 (R40-82), IgM (II/41), and human CD2 (RPA-2.10) antibodies.
The rabbit polyconal anti-H3K27ac antibody (ab4729; Abcam) was used for ChIP experiments, and a rabbit polyconal anti-E2A antibody (Kwon et al., 2008) was used for ChIP and immunoblot analysis. The anti-E2A antibody (produced in-house) was directed against the N-terminal peptide DRPSSGSWGSSDQNSSSFDP of the mouse E2A protein, which is absent in the HEB and E2-2 proteins.
Definition and flow cytometric sorting of B cells, plasmablasts, and plasma cells
Mature B cells from lymph nodes, long-lived plasma cells from the bone marrow and in vitro–differentiated activated B cells and plasmablasts were sorted with a FACSAria machine (BD) as follows: immature B (B220+CD19+IgMhiIgDloCD21−), mature B (B220+CD19+IgMloIgDhi), FO B (B220+CD19+CD21intCD23hi), MZ B (B220+CD19+CD21hiCD23lo/−), B-1 (B220loCD19+), GC B (B220+CD19+GL7+Fas+), plasma cells (Lin−B220intCD138hiCD28+), in vitro–activated B cells (CD22+CD138−), preplasmablasts (CD22−CD138−), and plasmablasts (CD22−CD138+). The Lin marker antibodies contained anti-CD4, anti-CD8a, anti-CD11b, anti-CD21, and anti-DX5 antibodies for the analysis of plasma cells.
In vitro B cell stimulation experiments
FO B cells from the lymph nodes (RNA- and ATAC-seq) or spleen and lymph nodes (E2A-ChIP and 3C analyses) were isolated as CD43− B cells by immunomagnetic sorting. For LPS stimulation, FO B cells were seeded at a density of 5–50 × 104 cells/ml in IMDM medium containing 10% fetal calf serum (A15-101; GE Healthcare), 1 mM glutamine, 50 µM β-mercaptoethanol, and 25 µg/ml LPS (L4130; Sigma-Aldrich). For anti-CD40 plus IL-4 stimulation experiments, FO B cells were plated at 5 × 105 cells/ml in IMDM medium supplemented with 10% fetal calf serum, 1 mM glutamine, and 50 µM β-mercaptoethanol containing anti-CD40 antibody (2 µg/ml; HM40-3; eBioscience) and IL-4 (20 ng/ml) for up to 3 d. For cell proliferation analysis, the purified B cells were first stained with the CellTrace Violet reagent (5 µM; Invitrogen) before stimulation as described above. For LPS plus IL-4 stimulation, lymph node B cells were treated with LPS (25 µg/ml) and IL-4 (20 ng/ml) for 4 d. For the stimulation experiments shown in Fig. 6 D, the reagents were used at the following concentrations: anti-CD40 antibody (2 µg/ml), IL-4 (20 ng/ml), IL-5 (10 ng/ml), and IL-21 (10 ng/ml).
For experiments using the iGB system (Nojima et al., 2011), 40LB cells were cultured in DMEM medium supplemented with 10% fetal calf serum. Splenic B cells (105) were plated on irradiated 40LB feeder cells in one well of a 6-well plate in RPMI medium supplemented with 10% fetal calf serum, 1 mM glutamine, 50 µM β-mercaptoethanol, 10 mM Hepes, and 1 mM sodium pyruvate, followed by stimulation with IL-4 (20 ng/ml) for 4 d. At day 4, activated B cells (8 × 104) were transferred to fresh 40LB cells in one well of a 6-well plate and stimulated for another 4 d with IL-21 (10 ng/ml).
Sheep RBCs (SRBCs) were washed in PBS and resuspended at 109 cells/ml, followed by intraperitoneal injection of 100 µl into an adult mouse. To study the immune responses to a T cell–independent antigen, mice were intraperitoneally injected with 10 µg of TNP (2,4,6-trinitrophenyl)-Ficoll (Biosearch Technology) in PBS. The immune response to a T cell–dependent antigen was investigated by intraperitoneal injection of 100 µg NP-KLH in alum.
The frequencies of IgM antibody-secreting cells (ASCs) were determined by enzyme-linked immunospot (ELISPOT) assay as described (Smith et al., 1997). Goat anti-mouse IgM antibody-coated plates were used for capturing IgM antibodies secreted by individual cells, respectively. Spots were visualized with goat anti-mouse IgM antibodies conjugated to alkaline phosphatase (Southern Biotech), and color was developed by the addition of BCIP/NBT Plus solution (Southern Biotech). After extensive washing, the spots were counted with an AID ELIspot reader system (Autoimmun Diagnostika).
Cryosections of the spleen from immunized mice were stained with a biotinylated anti-IgD antibody (1.19; produced in-house) and eFluor660-labeled GL7 antibody (GL-7; eBioscience). The biotinylated anti-IgD antibody was visualized by incubation with Cy3-streptavidin (Jackson ImmunoResearch Laboratories).
Cloning of retroviral expression constructs
For retroviral infection experiments, we cloned full-length mouse cDNA into the retroviral vector MiCD2 (Heavey et al., 2003) upstream of the IRES-hCD2 indicator gene to generate MiCD2-Blimp1, MiCD2-AID, MiCD2-V5-E2-2, and MiCD2-V5-E2A (containing full-length E47 cDNA). V5 refers to an N-terminal insertion of the V5 epitope tag.
Retroviruses were produced by transfecting 20 µg of the retroviral expression vector together with 10 µg of the retroviral helper vector pCMV-Gag-Pol into Plat-E packaging cells using standard calcium phosphate transfection in the presence of 25 µM chloroquine. 24 h after transfection, the high-titer viral supernatant was collected in B cell medium (IMDM supplemented with 10% fetal calf serum, 1 mM glutamine, and 50 µM β-mercaptoethanol). The infection was performed in a 6-well plate by spinning for 45 min at 2,400 rpm. Each well contained 106 lymph node B cells in 1 ml of freshly collected viral supernatant in B cell medium supplemented with LPS (25 µg/ml) or LPS (25 µg/ml) and IL-4 (20 ng/ml). The infection was repeated four times at intervals of 3–4 h in the presence of 4 µg/ml of polybrene except for the last infection. 4–12 h after the last infection, 5 ml of B cell medium containing LPS (25 µg/ml) and IL-4 (20 ng/ml) was added. The infected cells were analyzed 4 d after the start of culture.
RT-qPCR analysis of mRNA and nascent transcripts
Total RNA was prepared from sorted activated B cells, preplasmablasts, and plasmablasts by using the RNeasy Mini kit (QIAGEN). Genomic DNA was eliminated by using a gDNA eliminator spin column (QIAGEN). The cDNA was synthesized using the Random Primer Mix (New England Biolabs) and SuperScript III Reverse transcription (Life Technologies). Aicda transcripts were measured by qPCR using primers shown in Table S3, and were normalized against the Tbp transcripts. Nascent Prdm1 transcripts were measured by qPCR using primers that are located in intron 2 (Table S3), and were normalized against nascent Tbp transcripts (intron 1).
The 3C templates of activated B cells and plasmablasts or pro–B cells were prepared by using HindIII as the restriction enzyme as previously described (Medvedovic et al., 2013). DNA (200–400 ng) of individual 3C templates was subjected to quantitative TaqMan PCR analysis (qPCR; Hagège et al., 2007). The SensiFAST Probe No-ROX kit (Bioline Reagents Ltd.) was used for qPCR amplification, which was performed with the Bio-Rad CFX Connect Real-Time PCR Detection System using region-specific primers and a corresponding LNA Double-Dye probes with FAM and BHQ1 at the 5′ and 3′ ends, respectively (Table S3). As internal control for the quality of the 3C template, the ubiquitously expressed Ercc3 (XPB) locus was analyzed by qPCR (Splinter et al., 2006). The cross-linking frequencies at the Prdm1 and control Ercc3 loci were calculated by using HindIII-digested and randomly ligated BAC DNA of these loci as a standard for PCR amplification. The relative cross-linking frequency was determined as the ratio of the cross-linking frequency at the Prdm1 locus relative to the cross-linking frequency at the Ercc3 gene.
Mapping of open chromatin regions
Open chromatin regions (referred to as DHS sites) were mapped in FO B cells, activated B cells and preplasmablasts by the ATAC-seq method as described (Buenrostro et al., 2013) with the following modification. The nuclei were prepared by incubating cells with nuclear preparation buffer (0.30 M sucrose, 10 mM Tris pH 7.5, 60 mM KCl, 15 mM NaCl, 5 mM MgCl2, 0.1 mM EGTA, 0.1% NP-40, 0.15 mM spermine, 0.5 mM spermidine, and 2 mM 6AA) before the Tn5 treatment (4 µl of Nextera Tn5 transposase per 20,000 cells).
ChIP-seq analysis of histone modifications
For ChIP-seq analysis, activated B cells after 2 d of anti-CD40 plus IL-4 stimulation were used for ChIP with an anti-H3K27ac antibody (see Antibodies), as previously described (Schebesta et al., 2007). The ChIP-precipitated DNA (5 ng) was used for library preparation.
ChIP-seq analysis of E2A binding
In vitro–stimulated cells were subjected to cross-linking at room temperature either for 10 min with 1% formaldehyde (single cross-linking) or for 45 min with 2 mM disuccinimidyl glutarate (Sigma-Aldrich), followed by 10 min with 1% formaldehyde (double cross-linking). The chromatin was prepared as previously described (Kohwi-Shigematsu et al., 2012) with the following minor modifications. In brief, the nuclei were prepared from fixed cells before lysis with the 4% SDS solution. The lysed nuclei were subjected to 8M-urea gradient centrifugation (Kohwi-Shigematsu et al., 2012). The pelleted genomic DNA, cross-linked with proteins, were sheared with a Bioruptor sonicator (Diagenode) followed by immunoprecipitation using the anti-E2A antibody described in Antibodies. The ChIP-precipitated DNA (1–2 ng) was used for library preparation.
cDNA preparation for RNA sequencing
cDNA was prepared from RNA of in vitro–differentiated and ex vivo–sorted cells as previously described (Minnich et al., 2016).
Library preparation and Illumina deep sequencing
Sequence reads that passed the Illumina quality filtering were considered for alignment. For ChIP-seq, the reads were aligned to the mouse genome assembly version of July 2007 (NCBI37/mm9), using the Bowtie program version 12.5. For alignment of ATAC-seq reads, the Bowtie version 2.1.0 was used with the additional parameter −sensitive −X 5000. For RNA-seq, the reads corresponding to mouse ribosomal RNAs (NCBI GenBank and RefSeq accession nos. BK000964.1 and NR_046144.1, respectively) were removed. The remaining reads were cut down to a read length of 44-nt and aligned to the mouse transcriptome (genome assembly version of July 2007 NCBI37/mm9) using TopHat version 1.4.1 (Trapnell et al., 2009).
Database of RefSeq-annotated genes
Peak-to-gene assignment and calculation of RNA expression values were all based on the RefSeq database, which was downloaded from UCSC on January 10th, 2014. The annotation of immunoglobulin and T cell receptor genes were incorporated from Ensembl release 67 (Cunningham et al., 2015). Genes with overlapping exons were flagged and double entries (i.e., exactly the same gene at two different genomic locations) were renamed. Identical genes with more than one assigned gene symbol were flagged. Genes with several transcripts were merged to consensus genes consisting of a union of all underlying exons using the fuge software, which resulted in 25,726 gene models.
Peak calling of E2A ChIP-seq data and target gene assignment
E2A peaks were determined in the following manner. First, peaks were called in the double-cross-linked E2A ChIP-seq sample with a p-value of <10−10 and in the corresponding single-cross-linked E2A ChIP-seq sample with a p-value of <10−5 by using the MACS program version 188.8.131.52 (Zhang et al., 2008) with default parameters, a genome size of 2,654,911,517 bp (mm9), and a mature B cell input sample (14,951). Second, the overlap of the called double- and single-cross-linked peaks was determined with the Multovl program (Aszódi, 2012) by using a minimal overlap length of one and allowing for all possible overlaps, with the results being parsed and converted to tables with custom-made bash, perl, and R scripts. The E2A peaks identified by this overlap analysis were then assigned to target genes as described (Revilla-i-Domingo et al., 2012). For comparisons of E2A ChIP-seq data between different cell types, we down-sampled the reads of all ChIP-seq experiments to the ChIP-seq experiment with the lowest number of aligned reads before peak calling.
Analysis of ATAC-seq data
The ATAC-seq data were analyzed as follows. The DHS sites (open chromatin) were identified by peak calling using the MACS 2.0.10 program. The peaks from WT and DKO cells were compiled. The peak regions were defined as peak summit ± 500 bp. Peaks were assigned to genes as described for the ChIP-seq analysis. The read density was calculated as RPKM value for each peak.
Read density heat maps
Read densities were calculated as previously described (Minnich et al., 2016).
For motif discovery, we used the MEME-ChIP suite (version 4.9.1; Machanick and Bailey, 2011) to predict the most significant motifs present in the 300 bp centered at the peak summit of the top 300 sequences, as sorted by the fold enrichment score of the MACS program.
Analysis of RNA-seq data
For analysis of differential gene expression, the number of reads per gene was counted using HTseq version 0.5.3 (Anders et al., 2015) with the overlap resolution mode set to union. The datasets were grouped according to the performed stimulation experiments and analyzed using the R package DESeq2 version 1.4.1 (Love et al., 2014). Sample normalizations and dispersion estimations were conducted using the default DESeq2 settings. In detail, the following DESeq2 analyses were performed: FO B cell and all anti-CD40 and IL-4 stimulation samples were analyzed together considering the genotype, time and genotype over time effects (model design formula: “~genotype + time + genotype:time”). The genotype over time effect for each day was then tested based on log ratio tests (with a reduced formula of “~genotype + time"). FO B cell and 2-d LPS stimulation data were also grouped and analyzed in the same manner. 3-d LPS-stimulated activated B cell (WT and DKO), preplasmablast (WT), and plasmablast (WT) samples were considered for sample normalization, dispersion estimation, rlog normalization of the gene counts, and the comparison of 3-d LPS-stimulated DKO versus 3-d LPS-stimulated WT activated B cells (Wald test). The 500 most varying rlog normalized counts (option blind) of the 3-d LPS stimulation dataset were used for PCA in Fig. 8 D. In general, genes with an adjusted p-value <0.05 and a fold change >3, as well as a mean RPKM (averaged within conditions) >3 were considered to be significantly expressed. Immunoglobulin and T cell receptor genes were filtered from the list of significantly expressed genes and were also disregarded in the RPKM calculations.
For detailed expression analysis of the Igh and Igk loci (Figs. 5 E and S5, A and C), IghM, IghD, IghG, IghE, and IghA were manually curated to provide detailed exon annotations and were split into constant, membrane, secreted, hinge, and I exon regions. Ratios between secreted and membrane Igh transcripts were estimated by using the RPKM values of the secreted and membrane regions.
The RNA-seq, ChIP-seq, and ATAC-seq data (Table S4) are available at the Gene Expression Omnibus (GEO) repository under the accession no. GSE77744.
Online supplemental material
Figs. S1 and S3 describe the E2A-binding analysis of B cells stimulated with anti-CD40 plus IL-4 or LPS, respectively. Figs. S2 and S4 show the expression patterns of interesting E2A,E2-2 target genes identified in B cells stimulated with anti-CD40 plus IL-4 or LPS, respectively. Fig. S5 deals with the E-protein–dependent regulation of Igh and Igk during LPS-induced plasmablast differentiation. Tables S1 and S2 contain the RNA-seq data of all regulated E2A,E2-2 target genes identified in B cells stimulated with anti-CD40 plus IL-4 or LPS, respectively. Table S3 provides the oligonucleotide primer information. Table S4 describes all Illumina sequencing experiments generated for this study.
We thank D. Holmberg (Umeå University) for providing Tcf3fl/fl mice, S.L. Nutt (WEHI Melbourne) for Prdm1Gfp/+ mice, A. Tarakhovsky (Rockefeller University) for Prdm1fl/fl mice, D. Kitamura (Tokyo University) for 40LB cells, C. Theussl for blastocyst injection, K. Aumayr and her team for flow cytometric sorting and A. Sommer and his team (Vienna Biocenter Core Facility) for Illumina sequencing.
This research was supported by Boehringer Ingelheim, a European Research Council Advanced grant (291740-LymphoControl) from the European Union Seventh Framework Program, the Austrian Industrial Research Promotion Agency, and a German Research Foundation fellowship (WO 1972/1-1; M. Wöhner).
The authors declare no competing financial interests.
Author contributions: M. Wöhner did most experiments; H. Tagoh performed all E2A ChIP-seq, ATAC-seq, 3C-qPCR, and nascent transcript analyses; I. Bilic performed the initial phenotypic analyses; D. Kostanova Poliakova generated retroviral constructs and the Prdm1Mut/Mut mouse; M. Fischer and M. Jaritz performed the bioinformatic analysis of all RNA-seq and ChIP-seq data, respectively; M. Wöhner, H. Tagoh, and M. Busslinger planned the project, designed the experiments, and wrote the manuscript.
chromosome conformation capture
assay for transposase-accessible chromatin
class switch recombination
DNase I hypersensitive
4-hydroxy-3-nitrophenylacetyl-conjugated keyhole limpet hemocyanin
reads per kilobase of exon per million mapped sequence reads
reads per gene per million mapped sequence reads
T follicular helper
H. Tagoh and I. Bilic contributed equally to this paper.
I. Bilic’s present address is Baxalta Innovations GmbH, A-1221 Vienna, Austria.