Ten-eleven translocation (TET) enzymes oxidize 5-methylcytosine, facilitating DNA demethylation and generating new epigenetic marks. Here we show that concomitant loss of Tet2 and Tet3 in mice at early B cell stage blocked the pro- to pre-B cell transition in the bone marrow, decreased Irf4 expression and impaired the germline transcription and rearrangement of the Igκ locus. Tet2/3-deficient pro-B cells showed increased CpG methylation at the Igκ 3’ and distal enhancers that was mimicked by depletion of E2A or PU.1, as well as a global decrease in chromatin accessibility at enhancers. Importantly, re-expression of the Tet2 catalytic domain in Tet2/3-deficient B cells resulted in demethylation of the Igκ enhancers and restored their chromatin accessibility. Our data suggest that TET proteins and lineage-specific transcription factors cooperate to influence chromatin accessibility and Igκ enhancer function by modulating the modification status of DNA.
https://doi.org/10.7554/eLife.18290.001Cell lineage specification is typically accompanied by changes in DNA cytosine methylation, most of which occur at CpG dinucleotides in somatic cells (Klose and Bird, 2006; Ooi et al., 2009). This prototypical epigenetic modification, once regarded as static, is now known to be remarkably dynamic (Pastor et al., 2013; Wu and Zhang, 2014). DNA methyltransferases (DNMTs) convert cytosine to 5-methylcytosine (5mC)(Ziller et al., 2013); subsequently, proteins of the TET dioxygenase family oxidise 5mC to 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC)(He et al., 2011; Ito et al., 2011; Pastor et al., 2013; Tahiliani et al., 2009). These modified bases, together termed oxidized methylcytosines (oxi-mC), facilitate DNA demethylation and are also epigenetic marks in their own right (Mellén et al., 2012; Spruijt et al., 2013).
Much of the interest in TET proteins has centered around the possibility that oxidised methylcytosines are intermediates in one or more pathways of DNA demethylation (Pastor et al., 2013; Wu and Zhang, 2014). In one well-documented mechanism, the maintenance DNA methyltransferase complex DNMT1/UHRF1 complex restores symmetrical methylation to the hemi-methylated CpGs that are formed upon DNA replication (Ooi et al., 2009). However, 5hmC, 5fC and 5caC all inhibit this process, thus causing 'passive' (replication-dependent) DNA demethylation. In a second mechanism, 5fC and 5caC are substrates for excision by TDG and subsequent 'active' demethylation through base excision repair (He et al., 2011; Shen et al., 2013; Song et al., 2013a). Thus TET-mediated oxidation of 5mC to 5hmC, 5fC or 5caC facilitates replication-dependent DNA demethylation on the one hand, and replication-independent excision of 5fC and 5caC followed by their replacement with unmodified C on the other hand.
B cells are an essential component of the adaptive immune system, which are selected to expand and subsequently undergo somatic hypermutation and immunoglobulin (Ig) class switching on the basis of their ability to produce high affinity antibodies that provide immunity against pathogens (Alt et al., 2013; Bossen et al., 2012; Victora and Nussenzweig, 2012). B cell development in the bone marrow involves step-wise rearrangement of the genes encoding Ig heavy and light chains, which combine to form the B cell antigen receptor (BCR). The Ig heavy chain locus, Igμ, is rearranged first at the pro-B cell stage. Subsequently, the pre-BCR (formed by pairing the rearranged heavy chain with a surrogate light chain) reactivates Rag expression and stimulates rearrangement of an Ig light chain locus (Igκ or Igλ) at the pre-B cell stage; the Igκ locus predominates, accounting for over 95% of all productively rearranged light chains expressed on B cells in mice. Prior to rearrangement, there are two non-coding germline transcripts from the κ locus, Cκ and Igκ, which are closely correlated with the rearrangement potential of the locus.
As in many other cell types, extensive changes in DNA cytosine modification are observed at a genome-wide level during B cell development (Kulis et al., 2015). Specifically, Ig locus rearrangement requires activation of tissue-specific enhancers, a process typically accompanied by changes in DNA and histone modification. In mature B cells, only the rearranged allele of Igκ is demethylated, whereas the unrearranged allele remains methylated. It was thus postulated that methylation suppresses Igκ locus accessibility and that differential methylation may contribute to allelic exclusion; however, subsequent studies found that both κ alleles can be expressed in mature B cells, even though one allele is methylated (Levin-Klein and Bergman, 2014). Conversely, analysis of Dnmt1-deficient pre-B cells (which lose DNA methylation at both alleles of the Igκ locus) indicated that demethylation of the Igκ locus was not sufficient to induce locus rearrangement and expression (Cherry et al., 2000). Therefore, the relation between DNA methylation and Ig locus rearrangement remains to be fully elucidated at a molecular level.
Here we have examined the role of TET proteins in B cell development. Three TET enzymes, TET1, TET2 and TET3, are expressed in mammalian cells. In mice, Tet1 is highly expressed in the embryo and in embryonic stem cells (Kang et al., 2015; Koh et al., 2011), whereas Tet2 and Tet3 are abundantly expressed in somatic cells including hematopoietic cells. Using mice doubly deficient for Tet2 and Tet3 in early B cells, we show that TET function is required for developing B cells to transit from the pro-B to the pre-B stage. Tet2 and Tet3 regulate germline transcription and rearrangement of the Igκ light chain, and synergize with B lineage-specific transcription factors such as E2A and PU.1 to promote DNA demethylation and chromatin accessibility at B cell enhancers. TET catalytic activity is required for IRF4 expression, but IRF4 cannot restore Igκ germline transcription or the demethylation status of the Igκ enhancers in the absence of Tet2 and Tet3. Our results suggest that TET enzymes regulate multiple aspects of B cell development at the pro-B to pre-B cell transition, and demonstrate a causal relationship between DNA modification and chromatin accessibility.
Tet2 and Tet3 mRNAs are abundantly expressed at all stages of B cell development, whereas Tet1 mRNA is expressed at much lower levels (Ko et al., 2010)(Figure 1—figure supplement 1A, left, blue bars). However neither Tet2-/- mice (which are fully viable and fertile [Ko et al., 2011]) nor Tet3fl/fl Mb1Cre mice (which we generated to bypass the perinatal lethality of Tet3-/- mice [Gu et al., 2011]) displayed any striking B cell phenotypes (Figure 1—figure supplement 1B, C and D and not shown). To assess the effect of a profound loss of TET function in B cells, we generated Tet2-/- Tet3fl/fl Mb1Cre mice (here termed Tet2/3 DKO mice), in which a conditional Tet3 allele (Ko et al., 2015) is deleted in the context of a germline deletion of Tet2 at the transition from pre-pro B cells to pro-B cells (Hobeika et al., 2006). As judged by DNA dot blot using an anti-5hmC antibody, 5hmC levels were at least 4-fold lower in vitro-cultured pro-B cells of Tet2/3 DKO mice compared to wild type (WT) (Figure 1—figure supplement 1A, right).
Tet2/3 DKO mice showed a striking reduction in the percentages and numbers of B cells in the bone marrow compared to WT mice, with a partial block at the pro-B to pre-B transition (Figure 1). The percentage of B220+CD19+ cells in the Tet2/3 DKO bone marrow was substantially reduced (<50% of that in WT bone marrow) at 7–8 weeks and even more pronounced (<10%) at 11–12 weeks of age (Figure 1A). The percentages and numbers of pre-B cells (CD43lowB220+IgM-) and immature B cells (CD43lowB220+IgM+) in the Tet2/3 DKO bone marrow at 11–12 weeks were 7–20% of those in the WT bone marrow (Figure 1B–D); concomitantly, the percentages and numbers of re-circulating (mature) IgM+IgD+CD19+ B cells in the bone marrow were also greatly diminished in Tet2/3 DKO mice (Figure 1C,D). Because B220 and CD43 are co-expressed not only on B cells but also on plasmacytoid dendritic cells, we reanalyzed CD19+B220+ bone marrow cells based on c-kit and CD25 expression; this analysis confirmed that percentages and numbers of pre-B cell (IgM-CD19+B220+ckit–CD25+) were substantially reduced in Tet2/3 DKO mice (Figure 1E). In parallel, Tet2/3 DKO mice showed an increased percentage of pro-B cells (IgM-CD19+B220+ckit+CD25–) in the bone marrow (Figure 1E, left), but total pro-B cell numbers were unaltered because of the overall decrease in total B-lineage cells (Figure 1E, right). Consistent with these findings, there was a reduction in the percentage and number of mature B cells in the spleen (Figure 1F).
(A) Reduced bone marrow B cells (B220+ CD19+) in Tet2-/-Tet3fl/fl Mb1Cre (Tet2/3 DKO) mice. Total bone marrow cells from wild type (WT) or Tet2/3 DKO mice at eight weeks (upper) and 11 weeks (lower) …
(A) Right, floxed Tet3 alleles are efficiently deleted in Tet2-/-Tet3fl/fl pro-B cells transduced with Cre-IRES-GFP retrovirus. Expression of Tet1, Tet2 and Tet3 mRNA in Cre+ (GFP+) WT and Tet2-/-Tet…
A large fraction of CD19+B220+ B cells in the spleen of Tet2/3 DKO mice lacked cell-surface IgM and/or IgD expression (~25% and~45% IgM–IgD– cells in eight week-old and 11 week-old Tet2/3 DKO mice respectively; Figure 1G) These peripheral Ig-negative B cells expressed a significantly higher level of Terminal deoxynucleotidyl transferase (TdT) and pre-BCR (VpreB, also known as CD179α), thus displaying the expression profile of developing pro-B cells (Figure 1H); they also uniformly lacked expression of Tet3 mRNA, indicating complete deletion of the Tet3 allele (Figure 1—figure supplement 1E). In contrast, two out of 3 surface Ig-positive cell samples analyzed showed residual expression of Tet3 mRNA (Figure 1—figure supplement 1E). Together, these data suggest that the surface Ig-positive cells in Tet2/3 DKO mice were ‘escapees’ that had not completely deleted the Tet3 allele, and hence had expanded due to the proliferative advantage of B cells expressing a cell-surface B cell receptor (BCR) (Kraus et al., 2004).
Notably, there was a detectable expansion of CD11b+ myeloid-lineage cells in the bone marrow and spleen of Tet2/3 DKO mice (Figure 1—figure supplement 1F), resembling the myeloid skewing observed in mice deficient in Tet2 or Tet3 alone in the hematopoietic compartment (Ko et al., 2011, 2015). Moreover, all older Tet2/3 DKO mice developed B cell lymphomas with splenomegaly and lymphadenopathy by the time they were five months old (Figure 1—figure supplement 1G). These findings are reported here for completeness, but will be pursued in a separate study. The data recall the late B cell malignancies observed in germline Tet1-deficient mice (Cimmino et al., 2015), and reinforce the notion that TET loss-of-function is associated with hematopoietic malignancies (Huang and Rao, 2014; Ko et al., 2015).
The splenic Tet2/3 DKO B cells that lacked cell-surface IgM and/or IgD expression (Figure 1G) also showed greatly reduced expression of cell-surface and intracellular Igκ light chains (Figure 2A) and impaired Igκ rearrangement (Figure 2B). Furthermore, purified pre-B cells from Tet2/3 DKO mice showed a marked reduction (>50%) in Vκ-Jκ1 rearrangement compared to their WT counterparts (Figure 2C), demonstrating that TET function is important for Igκ locus rearrangement and κ chain expression in developing B cells and hence critical for the generation of mature B cells with functional B cell receptors.
(A) IgM-IgD- B cells in the Tet2/3 DKO mice have diminished immunoglobulin light chain expression. IgM-IgD- and IgM+IgD+ splenic B cells (CD19+B220+) were analyzed for intracellular (left panel) and …
TET proteins are known to facilitate both active and passive DNA demethylation at gene bodies, promoter-TSS regions and enhancers ([Hon et al., 2014; Pastor et al., 2011; 2013; Tsagaratou et al., 2016] and below). To ask whether loss of TET function in Tet2/3 DKO B cells affected Igκ expression through changes in DNA modification, we analyzed the DNA modification status of known cis-regulatory elements in the Igκ locus (Inlay et al., 2002; Zhou et al., 2010) (Figure 3A, Figure 3—figure supplement 1A). The 3’ enhancer (3’Eκ) and a more distal 3’ enhancer (dEκ) contain 2 and 3 CpG dinucleotides respectively; the intronic Eκ enhancer (iEκ) contains no CpGs and was not considered further (Figure 3A, Figure 3—figure supplement 1A). Using bisulfite sequencing (which does not distinguish 5hmC from 5mC, or 5fC and 5caC from unmodified C [Huang et al., 2010; Pastor et al., 2013]), we found that all 5 CpG sites in the 3’ and distal Eκ enhancers were fully 'methylated' (5mC + 5 hmC) in LSK (Lin-, Sca1+ c-kit+) cells, a population enriched in hematopoietic stem/precursor (HSPC) cells (Hodges et al., 2011); began to lose this modification at the pre-pro-B cell stage; were partially methylated in pro-B cells; and almost completely 'unmethylated' (C/5fC/5caC) in WT splenic CD19+ B cells (Figure 3—figure supplement 1B). Examination of publicly available data for human HSC, neutrophils and B cells showed that the corresponding human enhancers also undergo B cell-specific 'demethylation' (loss of 5mC/5hmC; Figure 3—figure supplement 1C).
(A) The Igκ locus. Diagram depicts gene segments and regulatory elements of the Igκ locus. For sequences of the 3’ and distal Eκ enhancers, see Figure 3—figure supplement 1A. (B) Flow-chart …
(A) In vitro-derived WT and Tet2/3 DKO pro-B cells express similar level of pro-B surface markers. Cell-surface expression of B220, CD43, IgM, and CD127 on WT, control (Tet2-/- Tet3f/f) and DKO (Tet2…
(A) DNA sequences of the 3’Eκ and distal Eκ enhancers. CpG dinucleotides, consensus E-box, IRF, PU.1, and FKH motifs are indicated. (B) The 3’Eκ and distal Eκ enhancers undergo B cell-specific DNA …
To analyze the defects in Tet2/3 DKO B cells at a molecular level, we used a well-characterized cell culture system (Vieira and Cumano, 2004) in which TET loss-of-function could be acutely induced. Tet2-/-Tet3fl/fl bone marrow cells were cultured for one week with IL-7 and OP-9 feeder cells to generate pro-B cells (also known as large pre-B or pre-B II cells), then retrovirally transduced with Cre-IRES-GFP or control empty IRES-GFP retrovirus to delete the floxed Tet3 allele (Figure 3B). The efficiency of Tet3 deletion was >90%, and there was only a slight compensatory upregulation of Tet1 mRNA (Figure 1—figure supplement 1A, left, red and green bars). Under these conditions, expression of surface markers in Cre-transduced Tet2-/-Tet3fl/fl pro-B cells (termed Tet2/3 DKO pro-B cells) was indistinguishable from those in control cells (WT cells transduced with Cre or Tet2-/-Tet3fl/fl cells transduced with empty vector), with all genotypes exhibiting the typical surface phenotype of pro-B cells (B220+, CD19+, CD43high, IgM–, CD127+; Figure 3—figure supplement 2A and not shown).
Bisulfite sequencing showed that in vitro-generated Tet2/3 DKO pro-B cells displayed a striking increase in DNA modification (5mC + 5 hmC) at the 3’ and distal Igκ enhancers and a parallel decrease (>90%) in Cκ transcripts and Igκ germline (GL) transcripts compared to control pro-B cells (Figure 3C and D, compare blue and green bars), consistent with the corresponding decreased Igκ expression and increased DNA modification (5mC + 5 hmC) observed in Tet2/3 DKO B cells compared to WT B cells analyzed ex vivo (Figures 2A and 3E). The increased modification indeed reflected increased DNA methylation (5mC, not 5hmC), since genome-wide 5hmC mapping by CMS-IP showed almost no remaining 5hmC in Tet2/3 DKO pro-B cells compared to WT pro-B cells (Figure 3—figure supplement 2B); this is a more reliable technique in terms of signal-to-noise ratio than DNA dot blot analysis with anti-5hmC (Figure 1—figure supplement 1A, right). In vitro-generated pro-B cells deficient in Tet2 alone (Tet2-/-Tet3fl/fl cells transduced with empty retrovirus, Figure 3C-3E, red bars) showed a partial increase in methylation (Figure 3C) and a partial reduction in Cκ and Igκ germline transcripts (Figure 3D) whereas loss of Tet2 in vivo had no apparent effect (Figure 3E), suggesting redundant functions of Tet2 and Tet3 in modulating the DNA modification status of these Igκ enhancers and showing that Tet3 can compensate fully for lack of Tet2 in vivo but not in vitro.
Igκ germline transcription and Igκ rearrangement is regulated by the concerted action of multiple transcription factors at the Eκ enhancers, including E2A, EBF1, PAX5, PU.1, IRF4 and IRF8 (Clark et al., 2014; Cobaleda et al., 2007; Hagman and Lukin, 2006; Murre, 2005)(Figure 3A). With the exception of the functionally redundant transcription factors IRF4 and IRF8 which showed decreased expression (see below), mRNA and protein levels for most of these known regulators were unchanged in Tet2/3 DKO pro-B cells compared to WT (Figure 3D, Figure 3—figure supplement 2C and not shown). Similarly, Tet2/3 DKO pre-B cells transformed with BCR-ABL (Chen et al., 1994; Klug et al., 1994) had diminished Igκ germline transcription, both under resting conditions and after treatment with the Abelson kinase inhibitor Gleevec (also known as imatinib or STI-571), which induces Igκ locus activation in BCR-Abl-transformed cells (Muljo and Schlissel, 2003) (Figure 3—figure supplement 2D). Gleevec-treated Tet2/3 DKO pre-B cells also had significantly reduced formation of Rag-induced DNA double-strand breaks at the Jκ1 segment compared to WT cells (Figure 3—figure supplement 2E). These results show that TET proteins demethylate and modulate the activity of Igκ enhancers and thus influence subsequent Igκ rearrangement.
Ig germline transcription requires locus accessibility (Yancopoulos and Alt, 1985). To test whether TET proteins had a role in regulating the accessibility of the Igκ locus, we identified accessible chromatin regions in cultured pro-B cells at a genome-wide level by ATAC-seq (assay for transposase-accessible chromatin using sequencing) (Buenrostro et al., 2013) (Figure 4). The majority of accessible regions (n = 44668, 79.8%; see Figure 4 legend for details) were ‘commonly accessible’, i.e. similarly accessible in WT compared to Tet2/3 DKO pro-B cells (Figure 4A, grey dots). Importantly, 1303 regions (2.2%) were potentially regulated by TET because they were more accessible in WT than in Tet2/3 DKO pro-B cells (Figure 4A, WT>DKO differentially accessible regions (WT>DKO DARs); blue dots). 1192 regions (2.1%) gained accessibility in the Tet2/3 DKO cells (Figure 4A, DKO>WT differentially accessible regions (DKO>WT DARs); red dots), and are likely affected indirectly by TET (see below).
(A) ATAC-seq identifies TET-regulated accessible regions in the genome. WT and Tet2/3 DKO pro-B cells were generated as described in Figure 3 and used to prepare ATAC-seq libraries with two …
Genome browser views of the Foxo1 (A) and CD79α (C) loci in WT and Tet2/3 DKO pro-B cells. (A, C) From top to bottom are ATAC-seq peaks and 5hmC (CMS) in WT and Tet2/3 DKO pro-B cell, followed by …
ChIP-seq signals (Heinz et al., 2010) of CTCF, EBF, Foxo1, and Oct2 (in mature B cells) are plotted for the following regions as in Figure 4F: TET-regulated WT>DKO differential accessible regions …
TET-regulated WT > DKO accessible regions (WT>DKO DARs; more accessible in WT than in Tet2/3 DKO pro-B cells) (A) and DKO > WT accessible regions (DKO>WT DARs) (B) were analyzed for enrichment of …
To correlate TET activity and chromatin accessibility, we analyzed genome-wide 5hmC distributed by CMS-IP (Figure 4B). Compared to randomly sampled regions (Figure 4B, dotted line), all accessible regions displayed higher levels of 5hmC (solid lines), suggesting an overall positive correlation between TET activity (5hmC deposition) and chromatin accessibility. DARs that were more accessible in WT than in DKO pro-B cells (TET-regulated DARs, WT>DKO DARs or simply, WT DARs) exhibited significantly higher levels of 5hmC compared to other accessible regions (Figure 4B), whereas commonly accessible regions and DARs that were more accessible in DKO than in WT (DKO>WT DARs or simply, DKO DARs) had lower, comparable levels of 5hmC (Figure 4B, compare red and grey lines), further supporting the notion that TET activity is required to maintain chromatin accessibility. The majority (>86%) of TET-regulated DARs exhibit the characteristics of enhancers: they are located distal to transcription start sites (Figure 4C) and are enriched for H3K4me1 over H3K4me3 (Creyghton et al., 2010; Shlyueva et al., 2014) (Figure 4D). As expected from this global analysis, TET deficiency greatly decreased the accessibility of the dEκ enhancer (Figure 4E).
We mined the differentially accessible regions for consensus transcription factor binding motifs, identified by motif enrichment analysis (Heinz et al., 2010) and ChIP-seq. Comparing TET-regulated regions (WT>DKO DARs) to commonly accessible regions, the E-box/ basic region-helix-loop-helix (bHLH) (CANNTG) motif was highly represented, with the second most prominent motif being the composite ETS:IRF (PU.1-IRF) motif (Figure 4—figure supplement 1A). Consistent with their enrichment for bHLH and ETS binding motifs, WT>DKO DARs showed strong enrichment for PU.1 and E2A binding sites identified experimentally by ChIP-seq (Heinz et al., 2010; Lin et al., 2010), compared to DARs more accessible in DKO than in WT, commonly accessible regions and control random regions (Figure 4F,G) In contrast, the distribution of binding sites for other relevant transcription factors in pro-B cells (Heinz et al., 2010; Lin et al., 2010) either did not show a consistent relation to accessibility (EBF, Oct2) or were similar between WT>DKO DARs and other accessible regions (Foxo1; Figure 4—figure supplement 2).
Other examples of decreased accessibility in Tet2/3 DKO compared to WT pro-B cells are shown for the CD79α and Foxo1 loci; some of these sites coincide with E2A and/or PU.1 binding sites identified by ChIP-seq (Heinz et al., 2010; Lin et al., 2010), and CpG sites with low methylation (5mC + 5 hmC) at these regions in WT cells show increased methylation in Tet2/3 DKO cells (Figure 4—figure supplement 3). Together these data point to a significant association between chromatin accessibility, TET protein activity (i.e. presence of 5hmC) and the binding of two transcription factors, PU.1 and E2A, with key roles in B cell development.
Notably, consistent with a previous study (Benner et al., 2015), a fraction of the commonly accessible regions (especially those with the lowest ATAC-seq signals) corresponded to regions with significant CTCF binding in pro-B cells (Lin et al., 2010), suggesting these regions are potent insulators (Song et al., 2011). In contrast, CTCF binding sites were not enriched in the WT>DKO and DKO>WT DARs (Figure 4—figure supplement 2). These results suggest that regulation of insulator accessibility is likely TET-independent and imply that multiple mechanisms exist to regulate accessibility in different chromatin contexts.
In addition to mapping 5hmC by CMS-IP, we characterized the methylome of WT and Tet2/3 DKO pro-B cells by whole-genome bisulfite sequencing (WGBS), which does not distinguish 5hmC from 5mC, or 5fC and 5caC from unmodified C (Huang et al., 2010; Pastor et al., 2013). Based on the striking global loss of 5hmC in DKO pro-B cells as judged by the highly selective CMS-IP method (Huang et al., 2012; Pastor et al., 2012)(Figure 3—figure supplement 2B), all unconverted C’s identified by WGBS were assigned to 5mC. Notably, the methylation landscapes were similar between WT and Tet2/3 DKO pro-B cells, with 872 regions showing increased methylation and 258 showing decreased methylation in DKO pro-B cells (Figure 5). Analysis of these differentially methylated regions (DMRs) showed that DMRs with increased methylation in DKO cells compared to WT (DKO>WT DMRs) were concurrently less accessible (Figure 5A, left; Figure 5—figure supplement 1A, compare first and second columns), and had originally been modified with 5hmC in WT pro-B cells (Figure 5A, right; Figure 5—figure supplement 1A, see third column), suggesting that TET activity is required for antagonizing DNA methylation and maintaining accessibility at these regions. Consistent with the role of TET in regulating enhancer accessibility (shown above), the majority of regions with increased methylation in DKO cells compared to WT (72.5%) show characteristics of enhancers: location distal to transcription start sites (Figure 5B) and relative enrichment for H3K4me1 over H3K4me3 (Figure 5C) and were enriched for E2A and PU.1 binding (Figure 5—figure supplement 1A; compared to randomly sampled regions in Figure 5—figure supplement 1B). Together these results emphasize the strong correlation between chromatin accessibility, TET enzymatic activity based on generation of 5hmC, and maintenance of DNA demethylation at the enhancer regions marked by 5hmC. It has yet to be established, however, whether loss of 5hmC or increased 5mC plays a more significant role in the decreased chromatin accessibility in Tet2/3 DKO pro-B cells.
DNA isolated from WT and Tet2/3 DKO pro-B cells (two replicates each) were bisulfite-treated and sequenced to analyze the genome-wide DNA cytosine methylation. Differential methylation regions were …
(A–B) Differentially methylated regions (DMRs) were analyzed as described in Materials and Methods, and 827 and 258 regions were found to have more methylation in DKO (upper panel in (A); ‘DKO > WT …
Our results so far pointed to a strong functional interaction between TET proteins and two key transcription factors expressed in pro-B cells, E2A and PU.1: briefly, we observed TET activity (deposition of 5hmC and subsequent demethylation) at enhancers marked by the binding of E2A and PU.1, but not EBF, Foxo1 or Oct2 (Figures 4 and 5). We confirmed the functional interaction by transducing BCR-Abl1-transformed pre-B cells with lentiviruses encoding two independent shRNAs against E2A and PU.1 respectively, then examining the DNA modification status of the 3’ and distal Eκ enhancers by bisulfite sequencing (Figure 6). The shRNAs reduced the amount of E2A and PU.1 proteins substantially (Figure 6—figure supplement 1A and B), and increased the average level of DNA methylation at the first and second CpGs of the distal Eκ enhancer by 2- to 3-fold (Figure 6A and B), a magnitude similar to that observed in Tet2/3 DKO pro-B cells (Figure 3D).
(A) Depletion of E2A increases DNA modification (5mC + 5 hmC) at the distal Eκ enhancer. WT BCR-Abl-transformed pre-B cells were transduced with two independent shRNAs against E2A or scrambled shRNA …
PU.1 binding patterns in cultured WT and Tet2/3 DKO pro-B cells were analyzed by ChIP-seq, with two biological replicates for each genotype. (A) Similar global PU.1 binding in WT and Tet2/3 DKO …
WT Abl-transformed pre-B cells were transduced with two individual shRNAs targeting E2A (A) and PU.1 (B) and the respective protein level was detected by immunoblotting. Actin (A) and GAPDH (B) were …
Based on these data, we tested the possibility of a physical interaction between TET proteins and E2A and/or PU.1. We found that indeed, TET2 co-immunoprecipitated with E2A and PU.1 (Figure 6C); the interaction was direct, since we used benzonase (a nuclease that degrades both DNA and RNA) and the DNA intercalator ethidium bromide to prevent false positive interactions indirectly mediated by contaminating nucleic acids in cell extracts. These findings suggested that E2A and PU.1 recruit TET proteins to enhancers, where they increase chromatin accessibility by depositing 5hmC and facilitating DNA demethylation.
We tested the recruitment hypothesis by using ChIP-qPCR to evaluate the presence of TET2 at the 3’ and distal Eκ enhancers, in pro-B cells depleted of E2A or PU.1 (Figure 6A and B). Indeed, recruitment of TET proteins to the distal Eκ enhancer was significantly diminished in cells depleted for either E2A or PU.1 (Figure 6D and E). Notably, recruitment of TET2 to the 3’Eκ enhancer was not significantly affected (Figure 6D and E); this is entirely consistent with the fact that we observed increased methylation under the same conditions at the distal Eκ enhancer but no increase, or a much smaller increase, at the 3’Eκ enhancer (Figure 6A and B). The difference could be due to a requirement for additional factors (other than E2A or PU.1 alone) at the 3’Eκ enhancer, or to a threshold effect stemming from the stronger binding of E2A and PU.1 to the 3’Eκ compared to the dEκ enhancer (Figure 4E).
PU.1 has been proposed to be a 'pioneer' transcription factor capable of association with 'closed' chromatin (Ghisletti et al., 2010; Heinz et al., 2010). Consistent with this notion, we found that the genome-wide binding of PU.1 was virtually indistinguishable between WT and Tet2/3 DKO pro-B cells (Figure 6—figure supplement 2), suggesting that PU.1 is upstream of TET and enlists TET activity to modulate the accessibility of bound enhancers. Nonetheless, these data suggest a tight functional relationship between TET proteins and the key transcription factors during B cell development.
The availability of the in vitro system allowed us to ask whether TET catalytic activity was required for Igκ expression (Figure 7). For unknown reasons, we could not efficiently express full-length Tet2 or Tet3 in pro-B cells, but we did achieve short-lived expression of isolated Tet2 catalytic domain (Tet2CD) when we used a retroviral vector containing a blasticidin resistance gene (Figure 7A; the cells expressed Tet2CD for ~1 week after drug selection was initiated, after which expression declined rapidly). Under these conditions, Tet2/3 DKO BCR-Abl pre-B cells reconstituted with catalytically-active Tet2CD showed increased Igκ germline and IgCκ transcription (Figure 7B); they also displayed decreased DNA methylation at the 3’ and distal Eκ enhancers, compared to mock-transduced cells (Figure 7C). Cells transduced with a catalytically-inactive version of Tet2CD (HxD mutant (Ko et al., 2010; Tahiliani et al., 2009) were functionally inert (Figure 7B and C), despite the fact that the mutant protein was recruited to the 3’ Eκ and distal Eκ enhancers as or almost as effectively as the catalytically active version (Figure 7D). Expression of Tet2CD, but not its HxD mutant, also restored chromatin accessibility at the 3’Eκ and distal Eκ enhancers (Figure 7E) as well as at a genome-wide level (Figure 7F) in Tet2/3 DKO pro-B cells. Importantly, re-expression of Tet2CD was able to restore the VκJ rearrangement in Tet2/3 DKO pre-B cells (Figure 7G), demonstrating the functional relevance of enhancer accessibility to biological outcome.
(A–B) The enzymatic activity of Tet2 is required to promote Irf4 expression and Igκ germline transcription. BCR-Abl-transformed Tet2/3 DKO pre-B cells were transduced with retrovirus expressing …
In addition to decreased levels of Igκ germline transcripts, Tet2/3 DKO cells showed a significant decrease in Irf4 mRNA and protein levels (Figure 3D, Figure 3—figure supplement 2C). Expression of IRF4 protein and mRNA required TET catalytic activity, since re-expression of Tet2CD restored mRNA and protein levels of IRF4 whereas the Tet2CDHxD mutant did not (Figures 7B and 8A). Thus in addition to VκJ rearrangement, TET proteins regulate the expression of IRF4, which binds to 3’Eκ and potentially dEκ and regulates Ig light chain rearrangement and control pre-B cell development (Ma et al., 2006). Nonetheless, over-expression of IRF4 marginally induced Igκ gene transcription and did not affect 3’Eκ and dEκ methylation in the absence of TET (Figure 8B and C). These results suggested that IRF4 is induced via TET catalytic activity and then functions downstream of TET proteins in regulating early B cell development.
(A) Immunoblotting for Tet2 and IRF4 in Tet2/3 DKO BCR-Abl cells transduced with Tet2CD or Tet2CD HxD mutant. (B) Reconstitution with Irf4 only marginally rescued Igκ expression in the Tet2/3 DKO …
We have shown that mice with a germline deletion of Tet2 in conjunction with Mb1Cre-driven deletion of Tet3 display a block in B cell development at the pro-B to pre-B transition, reflective of decreased Igκ rearrangement and expression. A profound loss of TET function is necessary, since mice lacking only Tet2 or Tet3 in B cells have no obvious B cell phenotypes. By acute Cre-mediated deletion of Tet3 in cultured pro-B cells and rescue with Tet2, we have demonstrated that TET proteins concurrently regulate Igκ germline transcription and the DNA methylation status of the 3’ and distal Eκ enhancers, through a mechanism that requires TET catalytic activity.
There are well-established correlations among histone/ DNA modifications, chromatin accessibility and B cell development (Benner et al., 2015; Thurman et al., 2012). DNA cytosine demethylation, most likely mediated by TET enzymes, has a role in organizing genome domains by affecting the binding of CTCF (Benner et al., 2015; Flavahan et al., 2016). Moreover, conditional Dnmt3a/b deletion results in early Vκ-Jκ rearrangement, increased expression of Igκ germline transcripts, and decreased DNA methylation at Igκ enhancers (Manoharan et al., 2015). BRG1, the catalytic (ATPase) component of the SWI/SNF chromatin remodeling complex, and the chromatin reader BRWD1, have both been implicated in early B cell development (Bossen et al., 2015; Mandal et al., 2015). Indeed, Yancopoulos and Alt postulated more than 30 years ago that Ig locus accessibility is critical for Ig germline transcription and V(D)J rearrangement (Yancopoulos and Alt, 1985). Our data suggest that Tet2 and Tet3 maintain chromatin accessibility at Igκ locus enhancers as well as at many other promoters and enhancers in pro-B cells. While many protein complexes can modulate chromatin accessibility, the fact that acute deletion of TET function almost eliminates Igκ germline transcription, in a manner that can be rescued by a TET catalytic domain, implies that the TET enzymatic activity has a primary and essential role.
The expected consequence of decreased accessibility is decreased binding of transcription factors to the poorly accessible sites. We have shown that genomic regions that are less accessible in Tet2/3 DKO compared to WT pro-B cells are enriched for E2A and PU.1 binding motifs, and conversely, a subset of validated E2A and PU.1 binding sites, including the Igκ enhancers, are less accessible in Tet2/3 DKO pro-B cells. Notably, both transcription factors are involved in maintaining the demethylated status of the distal Eκ enhancer. PU.1, a key transcription factor in both the myeloid and lymphoid lineages, functions in macrophages as a 'pioneer' transcription factor capable of enhancing chromatin accessibility by binding linker regions between nucleosomes, thereby facilitating the modification of neighboring histones and the binding of other transcription factors (Ghisletti et al., 2010; Heinz et al., 2010). Consistent with a previous finding in macrophages that both TET2 and DNMT3b co-immunoprecipitate with PU.1 (de la Rica et al., 2013), we report here that TET2 co-immunoprecipitates with both PU.1 and E2A in developing B cells. Loss of TET2 and TET3 has a limited effect on the binding of PU.1 (Figure 6—figure supplement 2), suggesting that PU.1, and potentially E2A, recruit TET proteins to Igκ and other enhancers where they render the enhancers permissive for additional and/or more stable binding of those and other transcription factors (Figure 8D).
IRF4 and IRF8 are functionally redundant for the rearrangement of Ig light chains during pro-B to pre-B transition (Lu et al., 2003), while EBF1 is imperative for the pre-pro-B to pro-B transition (Lin and Grosschedl, 1995). Our data show that in addition to regulating Igκ germline transcription, TET2 and TET3 proteins are required for the expression of IRF4 and, to a lesser extent, IRF8. Ectopic IRF4 expression induced Igκ germline transcription as well as premature Igκ rearrangement in WT pro-B cells (Bevington and Boyes, 2013), but ectopic expression of IRF4 in Tet2/3 DKO cells only marginally increased Igκ germline transcription. Moreover, expression of IRF4 in Tet2/3 DKO cells had no effect on the DNA modification status of the Igκ enhancers, suggesting that IRF4 and IRF8 function at a late stage, after transcription factors such as PU.1 and E2A have cooperated with TET enzymes to facilitate enhancer accessibility. Interestingly, CpG dinucleotides are present in at least one potential IRF4 binding site in both the 3’Eκ and dEκ enhancers (Figure 3—figure supplement 1A), hinting that TET-mediated demethylation may be required for efficient binding of IRF4 at these enhancers. Finally, although we did not examine EBF1 because of lack of enrichment of its binding sites in TET-regulated accessible regions of the genome, EBF1 has been shown to increase both chromatin accessibility and DNA demethylation in B progenitor cells (Boller et al., 2016), and to interact with TET2 in a cancer cell line (Guilhamon et al., 2013).
Our reconstitution experiments show clearly that TET-mediated changes in DNA modification status and chromatin accessibility can be rescued by re-expression of the TET2 catalytic domain in Tet2/3 DKO pro-B cells. Although for technical reasons, we were unable to reconstitute the cells with full-length Tet2, it is not surprising that we observed a substantial restoration with the catalytic domain alone: many documented interactions of TET proteins with their partners involve the catalytic domain (e.g. TET1-OGT (Balasubramani and Rao, 2013), TET2-IDAX (Ko et al., 2013)), and Tet2 lacks a covalently associated CXXC domain that would be expected to target it to genomic regions containing unmethylated CpGs (Ko et al., 2013).
Our current study focuses on the function of Tet2 and Tet3 in early B cell development, but our data indicate that TET proteins also have a role in mature B cells. Specifically, we observed an age-related decline of mature B cells in Tet2/3 DKO mice, suggesting that Tet2 and Tet3 are important for B cell homeostasis. Despite this, we observed a fully penetrant development of B cell lymphomas in all Tet2/3 DKO mice, consistent with our recent finding that acute deletion of Tet2 and Tet3 in hematopoietic stem/ precursor cells results in acute myeloid leukemia (An et al., 2015). A recent study characterizing DNA methylation dynamics during human B cell differentiation revealed distinct demethylation patterns between early B cells and mature B cells (Kulis et al., 2015); demethylation in early human B cell development mainly occurs at enhancers, consistent with our observations in mice.
Bestor (Bestor et al., 2015) and Schübeler (Schübeler, 2015) have persuasively argued that DNA methylation and demethylation are not 'instructive' for gene transcription but rather are byproducts of altered transcription. In contrast, our data suggest that at least at Igκ enhancers in pro-B cells, the increase in DNA methylation that results from TET loss-of-function can be a determining ('instructive') factor for the observed decrease in Igκ germline transcription. However, our data also show that acute depletion of E2A results in increased DNA methylation at the distal Eκ enhancer, thus connecting transcription factor binding to DNA demethylation as postulated (Bestor et al., 2015; Schübeler, 2015). The simplest scenario is a positive feedback loop in which E2A, PU.1 and other transcription factors (de la Rica et al., 2013; Guilhamon et al., 2013; Wang et al., 2015) recruit TET proteins to enhancers, thus facilitating TET-mediated DNA modification; the recruited TET proteins, via DNA demethylation, maintain chromatin accessibility at these regions, promoting increased binding of E2A and other transcription factors such as PU.1 (Figure 8D). Future studies will reveal how lineage-specific transcription factors, chromatin regulatory proteins and TET enzymes together orchestrate cell fate determination and gene expression in B cells and other cell types.
Tet2fl/fl, Tet2-/- and Tet3fl/fl mice were generated as described (Ko et al., 2011, 2010, 2015). Mb1Cre mice were obtained from Dr. Changchun Xiao and originally from Dr. Michael Reth (Hobeika et al., 2006). All mice were in the C57BL/6 background.
Total RNA was isolated with RNeasy plus kit (Qiagen, Germnay) and cDNA was synthesized using SuperScript III reverse transcriptase (Thermo Fisher, Waltham, MA). Quantitative PCR was performed using FastStart Universal SYBR Green Master mix (Roche, Germany) on a StepOnePlus real-time PCR system (Thermo Fisher). Gene expression was normalized to Actb. Primers were listed in Supplementary file 1.
Total bone marrow cells from Tet2-/-Tet3fl/fl mice were cultured in ɑ-MEM medium with 20% FBS and IL-7 (Peprotech, Rocky Hill, NJ; 10 ng/ml) on OP9 feeder cells for seven days, and transduced with retrovirus containing Cre-IRES-GFP. GFP+ cells were sorted 5–7 days after transduction and maintained as above. Cells were routinely monitored for mycoplasma contamination and were negative.
To obtain BCR-Abl transformed pre-B cells, pro-B cells were transduced with retrovirus encoding a BCR-ABL fusion protein (p210, Addgene). Three days after transduction, cells were removed from IL-7 and OP9. To induce Igκ locus rearrangement, cells were treated with Gleevec (10 μM) for indicated times, and double strand breaks at the Jκ1 RSS1 region were quantified by ligation-mediated PCR as described (Curry et al., 2005).
To knockdown E2A and PU.1 in WT BCR-ABL transformed pre-B cells, cells were transduced overnight with lentivirus generated from pLKO.1 based shRNAs, including shScramble, shE2A-1, shE2A-2, shPU.1–1, and shPU.1–2 similar to previously described (Martinez et al., 2015). Cells were selected with puromycin at 18 hr post-transduction and continuously for 72 hr. shRNA target sequences are as follow: shE2A-1, CTGCACCTCAAGTCGGATAAG; shE2A-2, TTTGACCCTAGCCGGACATAC; shPU.1–1, AACAGAGCTGAACAGTTTGGG; shPU.1–2,TTCTGATACGTCATGCGCTTG.
Nuclear extracts from WT and Tet2/3 DKO Abl-preB cells were immunoprecipitated with anti-PU.1 and anti-E2A in the presence of benzonase (Sigma, St. Louis, MO; final 500 U/mL) and ethithum bromide (Sigma; 100 μg/mL) to prevent indirect pulldown via DNA binding. Resulting proteins and 2.5% of input nuclear lysate were blotted with anti-TET2, -PU.1, and -E2A antibodies. HRP-conjugated polyclonal goat-anti-rabbit antibodies were used for detecting TET2. For PU.1 and E2A detection, blots were stripped, and probed with primary and then secondary monoclonal antibody specific for rabbit IgG conformation (Cell Signaling Technology, Danvers, MA) to minimize reactivity with IgH and IgL at around 50 and 20 kDa, respectively. Note that the overall signal was lower for E2A and PU.1 due to the secondary antibody used.
For cell surface staining, cells were stained with indicated antibodies in FACS buffer (1% FBS and 0.1% NaN3 in PBS) for 30 mins on ice, washed, and fixed in 1% paraformaldehyde in PBS. Cytofix/Cytoperm kit (BD Bioscience, Franklin Lakes, NJ) was used for intracellular staining. Samples were analyzed with LSR Fortessa (BD Biosciences) and FlowJo (FlowJo LLC, Ashland, OR). The following antibodies were used in this study and were purchased from BD, eBioscience (San Diego, CA), or Biolegend (San Diego, CA): CD19 (clone 6D5), B220 (clone RA-3-6B2), CD25 (PC61), IgM (R6-60), IgD (11–26c.2a), CD43 (S7), CD179α (R3), cKit (2B8), Igκ (187.1), Igμ (Il/41). The following antibodies were used for immunoblotting and ChIP: anti-E2A (Santa Cruz, Dallas, TX; V18), anti-PU.1 (Santa Cruz, T-21), anti-IRF4 (Santa Cruz, M-17) and anti-TET2 (Abcam, United Kingdom;ab124297).
To analyze Igκ rearrangements, a degenerate Vκ primer (GGC TGC AGS TTC AGT GGC AGT GGR TCW GGR AC) and an Igκ intron primer (AAC ACT GGA TAA AGC AGT TTA TGC CCT TTC) were used; Quantitation of Vκ-Jκ1 rearrangement was done by qPCR with these primers: degVκ: GGC TGC AGS TTC AGT GGC AGT GGR TCW GGR AC and κ-J1-R: AGC ATG GTC TGA GCA CCG AGT AAA GG.
CMS-IP was performed similar to pervious described (Huang et al., 2012; Pastor et al., 2011). Briefly, genomic DNA from in vitro cultured pro-B cells was isolated with PureLink Genomic DNA kit (Thermo Fisher) and were spiked-in with cl857 Sam7 λDNA (Promega,Madison, WI) and PCR-generated, hmC-containing puromycin-resistant gene at a ratio of 200:1 and 100,000:1, respectively. DNA was sheared with a Covaris E220 (Covaris), end-repaired, A-tailed, ligated with methylated Illumina adaptors (NEB, Ipswich, MA), and bisulfite-converted (MethylCode Bisulfite Conversion Kit, Thermo Fisher). Bisulfite-converted DNA was denatured and immunoprecipitated with anti-CMS serum. Immunoprecipitated DNA was PCR-amplified with barcoded primers (NEBNext Multiplex Oligos for Illumina, NEB) for 15 cycles with Kapa HiFi Uracil+ (Kapa Biosystems, Wilmington, MA). Resulting libraries were sequenced with a HiSeq2500 system for 50 bp paired-end reads (Illuminia, San Diego, CA). The sequence reads were mapped to mm9 with Bismap, and CMS-enriched genomic regions were identified using the ‘findPeaks' command in HOMER with the 'histone’ mode and default parameters (Heinz et al., 2010).
Genomic DNA was bisulfite-treated according to the manufacturer’s instruction (MethylCode Bisulfite Conversion Kit, Thermo Fisher), and target sequences were amplified by PCR. The PCR products were TA-cloned and sequenced as individual clones. Alternatively, PCR products were ligated with Illumina sequencing adaptors, re-amplified and sequenced using a Miseq platform (Illuminia). For Sanger sequencing of individual clones, at least 15 sequences were obtained from each sample and aligned to the target sequences using BLAST. For bisulfite sequencing using Miseq, reads were aligned to the targets with BS-Seeker2 (Guo et al., 2013), and the sequencing depth was adjusted to at least 20 times for each. Average bisulfite conversion efficiencies were above 99%, calculated based on C to T conversion frequencies at non-CpG sites.
Genomic DNA from in vitro cultured WT and DKO pro-B cells was prepared similar to above for CMS-IP without spiking-in hmC-containing DNA. The adapter-ligated fragments were amplified for four cycles, cleaned up, and sequenced for an approximately 30x coverage per base. We employed BSMAP v2.74 (Xi and Li, 2009) to align paired-end reads from bisulfite-treated samples to the mm9 mouse reference genome allowing four mismatches. Reads mapping to multiple locations in the reference genome with the same mapping score (multiple mappers) were removed as well as the 5’ ends bearing a quality lower than 20 (-v 4 w 2 r 0 -q 20). Bisulfite conversion efficiency was estimated based on cytosine methylation in non-CpG context. For all the samples the bisulfite conversion efficiency was higher than 0.9936. Duplicated reads caused by PCR amplification were removed by BSeQC v1.2.0 (Lin et al., 2013) applying a Poisson P-value cutoff of 1e-5. Consequently, a maximum of three stacked reads at the same genomic location were allowed and kept for further analysis. In addition, BSeQC was employed for removing DNA methylation artifacts introduced by end repair during adaptor ligation. Finally, overlapping segments of two mates of a pair were reduced to only one copy to avoid counting the same region twice during DNA methylation quantification.
To estimate CpG DNA methylation at both DNA strands, we executed the methratio.py script, from BSMAP v2.74 (Xi and Li, 2009) (-t 0 g 1 -x CG -i correct).
We used Regression Analysis of Differential Methylation (RADMeth; methpipe V.3.4.2) (Song et al., 2013b) for computing individual differentially methylated sites with the regression command, and posterior grouping into differentially methylated regions (DMRs) with the adjust and merge commands, (p.value 0.01).
In vitro cultured pro-B cells were fixed with 1% formaldehyde at room temperature for 10 min at 1×106 cells/mL in media and then quench with 125 mM glycine. Nuclei were isolated and sonicated with Covaris E220 for 1200 s. Resulting chromatin was pre-cleared, immunoprecipitated with the indicated antibodies overnight, washed, and digested with proteinase K (Thermo Fisher) and RNaseA (Qiagen) at 65°C to de-crosslink overnight. DNA was purified (Zymo ChIP DNA Clean and Concentration Kit; Zymo Research, Irvine, CA) and libraries were prepared with NEBNext Ultra kit (NEB) and sequenced on Illumina Hiseq 2500 for single-end 50 bp reads.
ATAC-seq libraries were prepared as described (Buenrostro et al., 2013). Briefly, 50,000 in vitro cultured pro-B cells were lysed in 25 μl ATAC-seq lysis buffer and tagged with transposase (Illumina, Nextera DNA Library Prep Kit). ATAC-seq libraries were sequenced on an Illumina Hiseq platform for paired-end 2×50 bp reads. The sequences were first mapped to the mouse genome (mm9) with bowtie (parameters: -m 1 -X 2000). The mapped reads were filtered to discard mitochondrial reads, and PCR duplicates were removed with Samtools rmdup. The accessible regions were then determined with MACS2 (-q 0.05), and the all peaks identified from two replicates from both WT and DKO were merged to generate a master peak set, which were used as regions of interest for subsequent analysis by MEDIPS (Lienhard et al., 2014). Differentially accessible regions were identified with FDR adjusted p value of 0.1. Commonly accessible regions used for analysis have adjusted p value higher than 0.1 and less than 4-fold difference between mean RPKM of two ATAC-seq replicates from WT and DKO. All mouse ENCODE blacklisted regions were removed (Consortium, 2012). The numbers for WT DARs, commonly accessible, and DKO DARs are 1303, 44668, 1192, respectively. Random genomic regions were generated with Bedtools shuffle with commonly accessible regions as input (n = 44668). HOMER annotatePeaks was used for identifying the distance from peaks to TSS. Motifs enriched in the preferential accessible regions were identified through HOMER (findMotifsGenome.pl –size given) with commonly accessible regions as background (Heinz et al., 2010).
To analyze the histone modifications and transcription factor binding associated with accessible regions, previous published ChIP-seq data were reanalyzed by mapping to reference genome (mm9) with Bowtie2 (Langmead and Salzberg, 2012). For histone marks enrichment analysis, differential and common peaks called by MEDIPS were extended for ±250 bp and the number of reads from H3K4me1 and H3K4me3 ChIP-seq datasets were retrieved and the RPKM was calculated by MEDIPS and used for scatter plot. Heatmaps and histograms were generated using Deeptools 2 (Ramírez et al., 2016) with a bin size of 10 bp (normalized using RPKM) and region of ±1500 bp from the peak center. Log2 ratio between WT and DKO merged ATAC-seq data was calculated by Deeptools 2 bamCompare.
The following published ChIP-seq data sets were used in this study: Oct2 (GSM537990), EBF (GSM546524), PU.1 (GSM537996), E2A (GSM546523), H3K4me1 (GSM546527), H3K4me3 (GSM546529), FOXO1 (GSM546525), CTCF (GSM546526), bisulfite sequencing (GSE31971).
All animal works were performed according to protocol (AP128-AR2-0516) approved by the Institutional Animal Care and Use Committee at La Jolla Institute.
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We would like to thank Dr. Cornelis Murre, Dr. Chris Benner, and Dr. Sven Heinz (UCSD) for comments and discussion; Dr. Changchun Xiao (Scripps Institute) for providing Mb1-Cre mice; Rao lab members for discussion and technical support; Dr. James Scott-Browne for advice on ATAC-seq analysis; Susan Togher for managing our mouse colony and for technical assistance; Cheryl Kim, Kurt van Gunst, Lara Nosworthy and Denise Hinz (LJI Flow Cytometry Core) for help with cell sorting; and Jeremy Day (LJI Functional Genomics Center) for assistance with next generation sequencing. This work was supported by NIH R01 grants AI44432 and CA151535 and LLS TRP grant 6187–12 (to AR). C-WJL was supported by a Cancer Research Institute Irvington Postdoctoral Fellowship. The authors declare no competing financial interests.
Animal experimentation: This study was performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All animal works were performed according to protocol (AP128-AR2-0516) approved by the Institutional Animal Care and Use Committee at La Jolla Institute.
© 2016, Lio et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
RNA interference (RNAi) is a conserved pathway that utilizes Argonaute proteins and their associated small RNAs to exert gene regulatory function on complementary transcripts. While the majority of germline-expressed RNAi proteins reside in perinuclear germ granules, it is unknown whether and how RNAi pathways are spatially organized in other cell types. Here, we find that the small RNA biogenesis machinery is spatially and temporally organized during Caenorhabditis elegans embryogenesis. Specifically, the RNAi factor, SIMR-1, forms visible concentrates during mid-embryogenesis that contain an RNA-dependent RNA polymerase, a poly-UG polymerase, and the unloaded nuclear Argonaute protein, NRDE-3. Curiously, coincident with the appearance of the SIMR granules, the small RNAs bound to NRDE-3 switch from predominantly CSR-class 22G-RNAs to ERGO-dependent 22G-RNAs. NRDE-3 binds ERGO-dependent 22G-RNAs in the somatic cells of larvae and adults to silence ERGO-target genes; here we further demonstrate that NRDE-3-bound, CSR-class 22G-RNAs repress transcription in oocytes. Thus, our study defines two separable roles for NRDE-3, targeting germline-expressed genes during oogenesis to promote global transcriptional repression, and switching during embryogenesis to repress recently duplicated genes and retrotransposons in somatic cells, highlighting the plasticity of Argonaute proteins and the need for more precise temporal characterization of Argonaute-small RNA interactions.
A new method for mapping torsion provides insights into the ways that the genome responds to the torsion generated by RNA polymerase II.