Neuron

Article Dnmt3a Regulates Global Gene Expression in Olfactory Sensory Neurons and Enables Odorant-Induced Transcription Bradley M. Colquitt,1,2 Eirene Markenscoff-Papadimitriou,1,2 Rachel Duffie´,2 and Stavros Lomvardas2,* 1Neuroscience

Graduate Program, University of California, San Francisco, San Francisco, CA 94158, USA of Anatomy, University of California, San Francisco, San Francisco, CA 94158, USA *Correspondence: [email protected] http://dx.doi.org/10.1016/j.neuron.2014.07.013 2Department

SUMMARY

During differentiation, neurons exhibit a reorganization of DNA modification patterns across their genomes. The de novo DNA methyltransferase Dnmt3a is implicated in this process, but the effects of its absence have not been fully characterized in a purified neuronal population. To better understand how DNA modifications contribute to neuronal function, we performed a comprehensive analysis of the epigenetic and transcriptional landscapes of Dnmt3adeficient mature olfactory sensory neurons (mOSNs), the primary sensory neurons of the olfactory epithelium. Dnmt3a is required for both 5-methylcytosine and 5-hydroxymethylcytosine patterning within accessible genomic regions, including hundreds of neurodevelopmental genes and neural enhancers. Loss of Dnmt3a results in the global disruption of gene expression via activation of silent genes and reduction of mOSN-expressed transcripts. Importantly, the DNA modification state and inducibility of odorant-activated genes are markedly impaired in Dnmt3a knockouts, suggesting a crucial role for this enzyme in establishing an epigenetic landscape compatible with neuronal plasticity. INTRODUCTION The sculpting of DNA modification patterns during cellular differentiation is essential to the construction of cellular identity. Once thought to be largely static in somatic tissues after embryogenesis, the patterning of cytosine modifications is now appreciated to be dynamic in certain cellular and genomic contexts. In particular, 5-hydroxymethylcytosine (5hmC), an oxidized derivative of 5-methylcytosine (5mC), is highly enriched in neurons (Globisch et al., 2010; Kriaucionis and Heintz, 2009; Ruzov et al., 2011); increases in abundance during neurogenesis (Hahn et al., 2013; Szulwach et al., 2011); and is localized to gene bodies, regions upstream of transcription start sites (TSSs), and enhancer elements (Colquitt et al., 2013; Hahn et al., 2013; Melle´n et al., 2012; Song et al., 2011; Szulwach et al., 2011). To date, it is unclear how neuron-specific 5hmC patterning is established during

the transition from neuronal progenitor to differentiated neuron. Recent studies (Hahn et al., 2013; Szulwach et al., 2011) noted that the increase of 5hmC during neurodevelopment is not strictly accompanied by a reduction of 5mC, suggesting that de novo DNA methylation and cytosine oxidation, mediated by the Tet family of enzymes (Wu and Zhang, 2011), are coupled during this developmental period. Moreover, substantial de novo 5mC patterning was recently found in cortical forebrain neurons relative to nonneuronal cell types (Lister et al., 2013). Dnmt3a, one of two de novo DNA methyltransferases in mammals, is expressed in neural precursor cells and neurons during late embryogenesis and in postmitotic neurons in the postnatal CNS (Feng et al., 2005). Its ablation specifically within the CNS (Nguyen et al., 2007) recapitulates the early death (3–4 weeks) seen in mice with constitutive loss of the enzyme (Okano et al., 1999). In addition, its loss within postmitotic neurons is associated with deficits in long-term potentiation, learning, and memory (Feng et al., 2010). Recent work identified a role for Dnmt3a in both the repression and facilitation of gene expression in neural stem cells (Wu et al., 2010), suggesting that the enzyme’s effects extend beyond establishing repressive 5mC. Here, we tested the hypothesis that Dnmt3a-mediated de novo DNA methylation contributes to 5hmC patterning within neurons and is necessary to define neuronal regulatory and transcriptional states. To explore this model in vivo, we analyzed the relationships between Dnmt3a-dependent 5mC and 5hmC patterning, chromatin accessibility, histone modifications, and transcription within both the mouse main olfactory epithelium (MOE) and the resident sensory neurons of the tissue, mature olfactory sensory neurons (mOSNs). We find that Dnmt3a is required for a restricted fraction of 5mC and 5hmC patterning within mOSNs that is strongly associated with highly accessible regions, including enhancer elements and TSSs. We also identify a set of regions exhibiting broad but precisely defined Dnmt3a-dependent 5mC and 5hmC patterning, high accessibility, and an association with neurodevelopmental genes. Interestingly, these Dnmt3a-dependent regions display intermediate levels of 5mC and 5hmC in Dnmt3a heterozygotes, indicating that levels of available Dnmt3a strongly influence DNA modification patterning. Moreover, the absence of Dnmt3a in mOSNs results in widespread transcriptional alterations, including the upregulation of repressed genes and modest downregulation of expressed genes. Finally, in the absence of Dnmt3a, immediate early genes display significant reductions of 5mC and 5hmC levels and have compromised Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc. 823

Neuron Dnmt3a Regulates Odorant-Induced Transcription

A

B

C

D

E

F

G

H

Figure 1. Dnmt3a Is Necessary for a Defined Fraction of Neuronal 5mC and 5hmC Patterning (A) RNA-seq fragments per kilobase per million (FPKM) of Dnmt family members in HBCs (multipotent stem cells), GBCs (neuronal progenitors), and mOSNs (primary sensory neurons). Data from Colquitt et al. (2013). (B) Basal-apical quantification of Dnmt3a immunoreactivity in Dnmt3a +/+, +/
, and
/
MOEs. Lines are average intensity across three biological replicates, and error is 95% confidence interval. Positions of GBCs and mOSNs determined using Neurog1-EGFP and Omp-EGFP MOEs, respectively. (C) Left: immunofluorescence of 5mC and 5hmC in Dnmt3a +/+ and
/
MOE. Right: quantification of 5mC and 5hmC intensity along basal-apical axis. Lines are average intensity across three biological replicates, and error is 95% confidence interval. Scale bar, 20 mm. (D) Representative locus showing 5mC and 5hmC DIP-seq tracks from Dnmt3a WT, Het, and KO mOSNs. Regions with significantly depleted 5mC or 5hmC are indicated below each modification set. (legend continued on next page)

824 Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc.

Neuron Dnmt3a Regulates Odorant-Induced Transcription

odorant-induced gene expression. These data support a model in which Dnmt3a-mediated DNA modification patterning is necessary to establish both repressive and active transcriptional states within neurons and to generate the full set of gene responses associated with neuronal plasticity. RESULTS We find that Dnmt3a transcription increases along the developmental lineage of OSNs, from horizontal basal cells (HBCs)—the multipotent stem cells of the tissue (Leung et al., 2007)—through globose basal cells (GBCs)—the neuronal progenitors of mOSNs (Caggiano et al., 1994)—and finally to mOSNs (Figure 1A), in agreement with previous work (MacDonald et al., 2005; Watanabe et al., 2006). Notably, this developmental increase is similar to the enrichment of 5hmC along the mOSN developmental path we described in an earlier study (Colquitt et al., 2013). In contrast, Dnmt3b, the other mammalian de novo DNA methyltransferase, is weakly expressed in HBCs but is not expressed within GBCs or mOSNs. Thus, Dnmt3a is the primary source of de novo 5mC within mOSNs. In 3-week-old Dnmt3a wild-type (WT) mice, Dnmt3a protein is most abundant in the immature neuronal stage, between GBCs and mOSNs, but is absent from basal stem and apical nonneuronal layers (Figure 1B and Figure S1A available online). Importantly, we do not find significant alterations in Dnmt3a knockout (KO) MOEs of the expression of the other DNA methyltransferases (except for a 2-fold increase of Dnmt3b expression [p = 0.02] that still remains at 5% of the expression level of Dnmt3a in WT) or of the three Tet members (Figure S1B). In addition, mOSN differentiation is largely unaffected in Dnmt3a KO MOEs, as determined by unaltered percentages of mOSNs and GBCs within the tissue (Figure S1C), unaffected mitosis rates (Figure S1D), and only mild increases in apoptosis rates (0.3% of total cells in WT to 0.9% in KOs, p = 0.11; Figures S1D and S1E). To determine the effects of the absence of Dnmt3a on global 5mC and 5hmC levels, we immunostained Dnmt3a WT and KO MOEs with antibodies specific for the two modified bases (Figure 1C). Interestingly, we find no significant difference in 5mC levels along the basal-apical axis (p = 0.39, Student’s t test, n = 3). However, we detect a significant reduction of 5hmC in Dnmt3a KO MOEs within the mOSN layer (p = 0.03, Student’s t test, n = 3), despite initially reaching levels similar to those in WT within more basal layers (Figure 1C). These results indicate that Dnmt3a is necessary to provide some component of de novo 5hmC patterning during the transition from neuronal progenitor to neuron. To determine what genomic regions exhibit Dnmt3a-dependent 5mC and 5hmC patterning, we isolated mOSNs via fluorescence-activated cell sorting (FACS) from the MOEs of 3-week-old Olfactory marker

protein (Omp)-ires-EGFP and Dnmt3a WT, heterozygote (Het), and KO mice and performed DNA immunoprecipitation against 5mC and 5hmC coupled to deep sequencing (DNA immunoprecipitation [DIP-seq]; Table S1). The specificity of these antibodies has been extensively demonstrated by us and other laboratories (Colquitt et al., 2013; Ito et al., 2010; Williams et al., 2011). We used two methods to identify genomic regions that contain Dnmt3a-dependent 5mC and 5hmC—peak analysis using a small peak size (500 bp) and differential enrichment analysis in broader windows (5–50 kb)—and found that the loss of Dnmt3a affects the modification status of a restricted set of genomic regions. For the first method, biological variation was controlled for by intersecting differential peak sets from two independent biological replicates. For the second method, variation across biological replicates is directly modeled during differential enrichment identification (Anders and Huber, 2010). Both sets of Dnmt3a-dependent 5mC (N = 955) and 5hmC (N = 5,767) 500 bp regions constitute only a small fraction of the total number of 5mC (N = 91,859) and 5hmC-enriched (N = 166,572) regions in mOSNs, 1% and 3%, respectively. A reciprocal analysis using identical peak cutoffs indicates that only two regions significantly gain 5mC and one region gains 5hmC in Dnmt3a KO relative to WT mOSNs. This specificity is readily seen in a representative 1.5 Mb locus in which the majority of 5mC and 5hmC patterning is unaffected by the loss of Dnmt3a except for a few sparsely distributed regions (Figure 1D). Genomewide, 5mC and 5hmC levels are highly correlated between Dnmt3a genotypes; however, 5mC and 5hmC levels in Het and KO mOSNs are reduced in regions that are highly modified in WT mOSNs (linear model slopes < 1; Figure 1E). We find little overlap between Dnmt3a-dependent 5mC and 5hmC peak sets (Figure 1F), reflecting their complementary distribution. However, regions containing Dnmt3a-dependent 5mC display mildly reduced 5hmC levels and vice versa (Figure 1G). Relative to total 5mC and 5hmC peak sets, Dnmt3a-mediated 5mC peaks are localized to CpG islands (CGIs), whereas Dnmt3a-mediated 5hmC peaks are associated not only with CGIs, but also with regions upstream of TSSs, and conserved intergenic regions (Figure 1H). In addition, both Dnmt3a-dependent 5mC and 5hmC are enriched at intergenic MOE enhancers, as described in more detail in Figure 5. Using the second differential enrichment method, we find large domains (average size 9 kb) of Dnmt3a-dependent 5mC (N = 1,162) and 5hmC (N = 1,805) patterning that are strongly associated with gene regions (2,025 of 2,549, p = 4 3 10
58, Fisher’s exact test). These domains, which we refer to as ‘‘Dnmt3a-dependent blocks,’’ are significantly associated with gene regions that are implicated in neuronal projection development and function (e.g., Ephrin family members,

(E) Scatter plots of Dnmt3a 5mC and 5hmC levels in Dnmt3a WT, Het, and KO mOSNs. Each point is the average RPKM within nonoverlapping genome-wide 10 kb windows. Shown is a random 10% subset of the data. Range of linear slopes across biological replicates is given for each comparison. (F) Overlap between Dnmt3a-dependent 5mC and 5hmC peak sets. (G) Levels of 5mC and 5hmC in either Dnmt3a-dependent 5mC or 5hmC peak sets. RPKM values have been square-root transformed for visual clarity. (H) Genomic feature-peak intersections represented as the ratio of the observed number of peaks that intersect with a given feature set over the expected number of peaks (modeled as peaks randomly positioned across the genome). Shown are intersections for total and Dnmt3a-dependent 5mC or 5hmC peak sets. Red labels indicate features with statistically greater proportional number of Dnmt3a-dependent peaks than total peaks (p < 0.01, Fisher’s exact test). See also Figure S1.

Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc. 825

Neuron Dnmt3a Regulates Odorant-Induced Transcription

A

B

C

E

D

F

Figure 2. Dnmt3a Is Required for 5mC and 5hmC Patterning within Genomic Blocks Associated with Key Neurodevelopmental Genes (A) Log2(KO/WT RPKM) of average 5mC and 5hmC RPKM within the union of the 5–50 kb regions containing Dnmt3a-dependent 5mC or 5hmC patterning, ‘‘Dnmt3a-dependent blocks’’ (N = 2,549). Genes of interest that intersect with these regions are indicated. (B) Densities of average vertebrate phyloP conservation across Dnmt3a-dependent blocks, 20 kb flanking both sides of each block, or an equal number of random genomic regions. (C) Densities of average DNase I hypersensitivity (DHS) RPKM in Dnmt3a WT MOEs across the same regions in (B). (legend continued on next page)

826 Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc.

Neuron Dnmt3a Regulates Odorant-Induced Transcription

Robo1 and Robo2, and amyloid beta precursor protein), cell adhesion (Pcdha cluster and Dscaml1), and neurodevelopmental transcriptional regulation (Foxg1, Otx2, Olig2, and Klf5; Figure 2A; Table S2). Of these, approximately 200 are transcription factors. Interestingly, 94 of the 1,944 associated genes are entirely contained within the region of Dnmt3a dependency. Dnmt3a-dependent blocks are highly conserved, containing significantly higher vertebrate phyloP scores (median ± SD, 0.11 ± 0.17) than flanking regions (0.09 ± 0.08, p = 0.02, twosided Wilcoxon) or an equal number of random genomic regions (0.05 ± 0.11, p = 8 3 10
6; Figure 2B). Similarly, these blocks are significantly more accessible to DNase I cleavage, a hallmark of active genomic sites, in Dnmt3a WT MOEs (median reads per kilobase per million [RPKM] ± SD, 0.67 ± 0.31) relative to flanking (0.55 ± 0.18, p < 2 3 10
16) and random regions (0.35 ± 0.36, p < 2 3 10
16; Figure 2C). To integrate these Dnmt3a-dependent blocks into the broader epigenetic landscape of the MOE, we performed chromatin immunoprecipitation sequencing (ChIP-seq) against H3K4me1 (a mark associated with genes and enhancer elements), H3K27ac (TSSs and active enhancer elements), and H3K27me3 (repressed loci) using native chromatin preparations isolated from the MOEs of 3-week-old Dnmt3a WT mice. Although the MOE is composed of multiple cell types, 45% of the tissue consists of mOSNs, making it the dominant cell type in this tissue (Figure S1C). We hierarchically clustered Dnmt3a-dependent blocks by their average H3K27ac and H3K27me3 levels, producing three clusters (enriched for H3K27ac only, H3K27me3 only, or both modifications; Figure 2D). This grouping by histone modifications resulted in a partial segregation of 5mC and 5hmC levels (Figure 2D): regions enriched for H3K27ac (cluster 1, example in Figure 2E, N = 1,219) are enriched for 5hmC, whereas those enriched for H3K27me3 (cluster 2, example in Figure 2E, N = 317) are correspondingly enriched for 5mC. Interestingly, a substantial fraction of these regions (N = 913) is enriched for both histone modifications and possesses high levels of both DNA modifications (cluster 3, examples in Figure 2F; Figures S2A–S2D). Notably, the set of genes associated with this dually modified cluster is enriched for genes implicated in neurodevelopment, including several transcription factors—Neurog1, Dlx5, Lhx2, and Emx2—known to be important for OSN maturation or olfactory receptor expression (Table S2) (Hirota and Mombaerts, 2004; Nicolay et al., 2006). In addition, DNase I accessibility within Dnmt3a WT MOEs is elevated precisely within the area of H3K27ac, H3K27me3, and Dnmt3a-dependent 5mC and 5hmC enrichment (Figures 2E and 2F; Figures S2A–S2D). This dual presence of histone modifications associated with activation and repression within single loci may be due to developmental heterogeneity within the tissue or may reflect differential monoallelic epigenetic states.

Interestingly, 5mC and 5hmC levels within Dnmt3a-dependent blocks exhibit an intermediate reduction in Dnmt3a Het mOSNs (Figures 3A and 3B). This observation could result from an admixture of highly modified and unmodified sequences, as might be generated by tissue or allelic heterogeneity, or by sequences with partially reduced methylation across individual alleles, which would indicate that the DNA modification patterning in these affected regions is sensitive to Dnmt3a dosage. To distinguish between these two models, we performed bisulfite sequencing using genomic DNA (gDNA) from Dnmt3a WT, Het, and KO MOEs on four gene regions that contain Dnmt3a-dependent 5mC or 5hmC (Egr1, Nr4a2, Fos, and the Dlx5 and Dlx6 intergenic region), one that contains Dnmt3a-independent 5mC and 5hmC (Hspa1a), and three that have enrichment for neither modification (TSSs of Egr1, Nr4a2, and Hspa1a; Figure 3C; Table S3). As bisulfite sequencing cannot distinguish between 5mC and 5hmC (Huang et al., 2010), these data provide a conflated measure of 5mC and 5hmC status. In agreement with the DIPseq data, Dnmt3a-dependent 5(h)mC exhibits intermediate levels in Dnmt3a Het MOE relative to WT and KO MOEs (Figure 3D). Moreover, individual sequencing clones from Dnmt3a Het MOEs display a partial loss of modification and are not in either fully modified or unmodified states (Figures 3D–3F; Figures S3A). In particular, this partial modification loss does not appear to be stereotyped across clones, consistent with a model in which Dnmt3a protein is limiting and acts by a nonprocessive mechanism (Gowher and Jeltsch, 2001). Although we detect little 5(h)mC in the CH context, known to be generated by Dnmt3a (Gowher and Jeltsch, 2001; Meissner et al., 2005; Ramsahoye et al., 2000), we do observe Dnmt3a-dependent 5(h)mC within the CA context at the Dlx5 and Dlx6 intergenic region (roughly corresponding to a previously characterized enhancer [Zerucha et al., 2000]). Our finding that Dnmt3a-dependent patterning is strongly associated with gene elements and with accessible transcription factor domains suggests that the enzyme generally modifies accessible regions of the neuronal genome. A recent study described a positive relationship between chromatin accessibility and neuronal DNA methylation in the CH dinucleotide context (Lister et al., 2013), which was hypothesized to be mediated by Dnmt3a. To examine this association, we compared 5mC and 5hmC levels in Dnmt3a WT and KO mOSNs with DNase I accessibility in Dnmt3a WT MOEs. In the absence of Dnmt3a, regions flanking DNase I hypersensitive sites in the MOE exhibit significant depletions of 5mC (Figure 4A; mean log2(KO/WT) ± SD, RPKM =
0.09 ± 0.52, p = 1 3 10
21, twosided Wilcoxon) and 5hmC (
0.15 ± 0.39, p = 2 3 10
48), and Dnmt3a-dependent 5mC and 5hmC peak sets significantly overlap with the set of accessible regions (Figure 4B; p = 1 3 10
9 and 3 3 10
259, respectively, Fisher’s exact test). Moreover,

(D) Hierarchical clustering of the union of Dnmt3a-dependent 5mC or 5hmC regions by average H3K27ac and H3K27me3 levels. Displayed are average 5mC and 5hmC RPKM in Dnmt3a WT and KO mOSNs and average DHS, H3K27ac, and H3K27me3 RPKM in Dnmt3a WT MOE. Each track is input subtracted and scaled to a common range. Numbered boxes indicate positions of three broad clusters from hierarchical clustering. Representative genes associated with each cluster are indicated at far right. (E) Tracks of representative genes found in cluster 1 (Rarg) and cluster 2 (Emx1). Shaded boxes denote regions containing Dnmt3a-dependent 5mC and/or 5hmC. y axis in RPKM. See also Figure S2. (F) As (E) with representative genes found in cluster 3 (Dlx5 and Dlx6). See also Figure S2.

Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc. 827

Neuron Dnmt3a Regulates Odorant-Induced Transcription

5mC

4

B

5hmC

log2 RPKM relative to +/+

A

Dnmt3a +/+

Density

3

2

1

0

C

0.4

0.8

1.2

1.6 0.5 RPKM

1.0

1.5

5mC

1

5hmC

0

+/+

+/-

-/-

+/+

+/-

-/-

Dnmt3a

+/+ 5hmC

+/-/+/+

5mC

+/-

-/Bisulfite region

D

Egr1 TSS

% (h)mCG

100

Egr1

Nr4a2

Nr4a2 TSS

Fos

Dlx5/6 intergenic

Hspa1a

Hspa1a TSS

Dnmt3a

75

+/+

50 25 0

E

CG

CA CC CG CT

(h)mCG 100 80 60 40

+/+

20 0

+/-

% (h)mC

Fos

Dnmt3a

100 80 60 40 20 0

100 80 60

-/-

40 20 0

chr12 86,816,841

86,817,584

F 100 80 60

+/+

40 20

+/-

% (h)mC

Dlx5/6

Dnmt3a

0 100 80 60 40 20 0 100

-/-

80 60 40 20 0

chr6

6,819,266

6,819,990

(legend on next page)

828 Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc.

Neuron Dnmt3a Regulates Odorant-Induced Transcription

global chromatin accessibility positively correlates with Dnmt3adependent 5mC (10 kb windows, Pearson R = 0.22, p < 1 3 10
4, permutation test) and 5hmC (R = 0.36, p < 1 3 10
4) patterning (Figures 4C and 4D). In a previous study (Colquitt et al., 2013), we overexpressed Tet3 in mOSNs (Tet3-tg) and found that regions with high levels of 5hmC in WT mOSNs exhibited depletions of the mark in Tet3-tg mOSNs, consistent with the role of Tet3 in progressive oxidation of modified cytosines (Wu and Zhang, 2011). Incorporating these data with the Dnmt3a data here, we find a positive correlation between Tet3-tg/WT differential 5hmC and Dnmt3a-dependent 5hmC (Pearson R = 0.46, p < 1 3 10
4; Figures 4C and 4D), suggesting that Dnmt3a operates on genomic loci that are also sensitive to neuronal Tet-mediated oxidation. As noted above, Dnmt3a-dependent 5mC and 5hmC are significantly associated with potential MOE enhancers (Figure 1H), defined as intergenic H3K4me1 peaks. We find positive correlations at these regions between flanking 5hmC levels and both H3K27ac levels (R = 0.33, p < 1 3 10
4, permutation test) and DNase hypersensitivity (R = 0.22, p < 1 3 10
4; Figure 5A) but only weak correlations between these features and flanking 5mC levels (R =
0.04 and 0.07, respectively, p < 1 3 10
4 for both). To further investigate the relationship between enhancer activity and 5hmC levels, we compared enhancer 5hmC levels and the expression of nearby genes along the mOSN developmental lineage. We clustered the 5,000 potential enhancers with the highest 5hmC variance across HBCs, GBCs, and mOSNs by their 5hmC levels (data from Colquitt et al., 2013), producing eight clusters each representing a developmental 5hmC path (Figure 5B) and computed the median expression levels of nearby genes. In support of an active role for 5hmC in enhancer function, we find a striking correspondence between developmental enhancer 5hmC patterning and nearby gene expression. Enhancers have been previously classified as ‘‘poised’’ (competent for activation, enriched for H3K4me1) or ‘‘active’’ (engaged with its target gene, enriched for both H3K4me1 and H3K27ac) (Creyghton et al., 2010). To generate a similar classification for MOE enhancers, H3K27ac peaks were intersected with intergenic H3K4me1 peaks, yielding H3K27ac-negative H3K4me1 peaks (poised) and H3K27ac-positive peaks (active). We find that 5mC levels in Dnmt3a KO mOSNs are reduced relative to WT flanking enhancers regardless of classification (Figure 5C). Similarly, 5hmC is depleted from the flanking regions of both poised and active enhancers in Dnmt3a KOs, although the depletion is strongest at regions flanking active enhancers. As with DNase hypersensitive sites in general, 5hmC in Dnmt3a Het mOSNs is present at intermediate levels flanking both poised and active enhancers.

The reduction of 5mC and 5hmC levels at active and poised enhancers in the absence of Dnmt3a may indicate that the enzyme modifies all enhancers or only those that are active within neural tissues. To resolve this, we computed Dnmt3a KO versus WT 5mC and 5hmC ratios at VISTA enhancers, functionally validated regulatory elements whose tissue activities have been determined by mouse reporter transgenics (Visel et al., 2007), and we found that those with Dnmt3a-dependent 5mC or 5hmC patterning are strongly associated with activity in the nose and other neural or face tissues (forebrain, midbrain, hindbrain, and neural tube; Figure 5D). Conditioning these enhancers by their tissue activity indicates that, although most enhancer types exhibit reduced 5mC and 5hmC in Dnmt3a KO mOSNs, those with neural activity display most substantial reductions (Figure 5E). Moreover, VISTA enhancers with at least 2-fold reduction of 5hmC are significantly associated with homeodomain transcription factors implicated in neurodevelopment (Table S4). To determine if Dnmt3a-dependent 5mC and 5hmC patterning is critical for the proper regulation of neuronal gene expression programs, we prepared RNA-seq libraries from FACS-isolated Dnmt3a WT, Het, and KO mOSNs with two biological replicates for each condition (Tables S1 and S5). Overall, RPKM values from each genotype are highly correlated (Figure S4A). To assess global gene expression, we averaged RPKM values across genotypes and replicates to obtain a transcriptional index. Plotting expression fold-change against this index, we find a striking bidirectional change of expression in Het and KO relative to WT mOSNs that is dependent on transcriptional output (Figures 6A and 6B). On average, lowly expressed or silent genes (as determined by Gaussian-based clustering of RPKMs) are transcriptionally upregulated in KO mOSNs (mean log2(KO/WT RPKM) ± SD = 0.9 ± 1.1), whereas highly expressed genes are weakly downregulated (
0.3 ± 0.4). Moreover, genes with an absolute KO-WT fold-change of at least 2-fold exhibit an intermediate disruption of expression levels in Het mOSNs, in concordance with Dnmt3a dosage response of 5mC and 5hmC levels (Figure 6C). To assess how Dnmt3a-sensitive genes are regulated during mOSN development, we clustered these gene sets by their expression in Dnmt3a WT mOSNs, GBCs, and HBCs, using previously published data sets (Figures S4B and S4C) (Colquitt et al., 2013; Magklara et al., 2011) and find that downregulated genes (N = 387) are expressed in mOSNs (Figure S4B), whereas upregulated genes (N = 1,845) are typically expressed in HBCs, GBCs, or not at all in the MOE (Figure S4C). To independently verify these effects, and to gauge if the general effects seen in the RNA-seq data are due to normalization artifacts, we performed RT-quantitative (q)PCR using RNA isolated from FACS-sorted Dnmt3a WT, Het, and KO mOSNs

Figure 3. 5mC and 5hmC Patterning in Dnmt3a-Dependent Regions Is Sensitive to Dnmt3a Heterozygosity (A) Densities of Dnmt3a WT, Het, and KO 5mC and 5hmC square-root RPKM (transformed for better visualization) within Dnmt3a-dependent blocks. Ticks indicate distribution medians. (B) Average 5mC and 5hmC RPKM normalized to WT within Dnmt3a-dependent blocks. Lines intersect with genotype medians. (C) 5hmC and 5mC DIP-seq tracks from Dnmt3a WT, Het, and KO mOSNs. Regions analyzed by bisulfite sequencing are indicated below each track. (D) Summary of percent modified CGs for each tested locus. Each point is the average CG modification percentage across individual sequencing clones. (E) Left: diagram depicting 5(h)mC status of individual CGs (columns) across sequencing clones (rows) for the Fos locus. Right: percent modification of cytosines across sequencing clones within all four dinucleotide contexts. See also Figure S3. (F) As (E) for the Dlx5 and Dlx6 intergenic region. See also Figure S3.

Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc. 829

Neuron Dnmt3a Regulates Odorant-Induced Transcription

A

C

B

D

Figure 4. Dnmt3a-Dependent 5mC and 5hmC Patterning Occurs in Accessible Chromatin (A) Average 5hmC and 5mC levels spanning MOE Dnmt3a WT DNase hypersensitive sites (N = 12,419). Error is bootstrapped 95% confidence intervals. (B) Overlap of Dnmt3a-dependent 5mC and 5hmC peak sets with DHS regions. (C) Chromatin accessibility and Dnmt3a-dependent DNA modification patterning across chromosome 19. Outside row: MOE Dnmt3a WT DNase accessibility. Second row from outside: difference between 5hmC levels in Tet3-tg mOSNs (Omp-tTA 3 tetO-Tet3-GFP) and control mOSNs (Omp-GFP). Data from Colquitt et al. (2013). Dnmt3a KO minus WT 5hmC (third row from outside) and 5mC (fourth row from outside) levels in mOSNs are shown. 5mC and 5hmC values are averages across biological replicates. y axes in RPKM. (D) Dnmt3a KO
WT 5mC (left) and 5hmC (center) RPKM versus Dnmt3a WT DHS. Right: Tet3-tg
control mOSN 5hmC RPKM (Colquitt et al., 2013) versus Dnmt3a KO
WT 5hmC RPKM.

830 Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc.

Neuron Dnmt3a Regulates Odorant-Induced Transcription

A

B

Enhancer

C

Gene

6

5

5

4

4

3

3

0.4 HBC GBC mOSN

0.8

3.5

2 1

1.5 HBC GBC mOSN

D

log2(FPKM+1)

1 0

5hmC 5mC -5 -5 kb b H3K4me1 H3K4me 3 3K4me1 1 +5 5 kb b -5 kb b H3K4me1 H3K4me 3K4me1 +5 +5 kb b k peak peakk

5hmC RPKM

1.6

2

0.5

RPK R RPKM PK PK

RPKM

6

1.0 0.1

7

H3K4me1 peak

E r r

MB

r

hear

hear

NEG

e NEG

HB

HB

r

e DRG

e DRG f

h

5hmC

CN MB

r e CN

cr

5mC

H3K4me1 peak

FB nose

r FB nose

f

Active

RPKM

7

+/+

0.2

2 0

DHS

0.4

H3K27ac

RPKM PKM P KM

0

8

RPKM RPK R PK PK

1

RPKM PKM P KM

8

Dnmt3a Poised 0.6

Closest

FM

FM

0

1

0

1

Figure 5. Dnmt3a Is Required for 5mC and 5hmC Patterning at Neuronal Enhancers (A) MOE H3K27ac, MOE DHS, mOSN 5hmC, and mOSN 5mC levels over MOE intergenic H3K4me1 peaks ordered by H3K27ac levels. (B) Analysis of the relationship between enhancer 5hmC levels (left) and the expression of closest genes (within 50 kb, right) within HBCs, GBCs, and mOSNs. (C) Average 5mC and 5hmC levels flanking poised and active enhancers (defined in text) in mOSNs from Dnmt3a +/+, +/
, and
/
MOEs. Error bars are bootstrapped 95% confidence intervals. (D) Association of enhancer activity with Dnmt3a-dependent 5mC and 5hmC. Shown are matrices of VISTA enhancer activities (rows) ordered by 5mC (left) or 5hmC (right) log2(Dnmt3a KO/WT RPKM) (columns). Each row was then loess smoothed and normalized, and the matrix was hierarchically clustered. (E) Distributions of average 5mC and 5hmC log2(KO/WT RPKM) across VISTA enhancers conditioned by tissue activity.

(n = 3 for each genotype). We tested the expression of 11 genes predicted to be downregulated, four genes predicted to be upregulated, and four predicted to remain the same according to KO/WT fold-changes from the RNA-seq data sets (Figure 6D). Of these genes, the expression of seven downregulated genes, one upregulated gene, and no unchanged genes displayed significant dependencies on Dnmt3a genotype (p < 0.05, one-way ANOVA). The remaining three upregulated genes exhibit a trend toward increased expression in Dnmt3a-deficient mOSNs by RT-qPCR. Furthermore, expression estimates from the RNAseq and RT-qPCR data sets are strongly correlated (Figure 6E; R = 0.76, p = 1 3 10
4, Pearson correlation test). Finally, in situ hybridizations (ISHs) using probes against genes that are expressed in mOSNs yet exhibit reduced expression in Dnmt3a KOs by RNA-seq (Omp and Lrrc23) display reduced expression (Figures 6F and S4D). As described above, Dnmt3a-dependent 5mC and 5hmC patterning is strongly associated with gene elements. To determine if modification changes correlate with transcriptional alterations, we examined the relationships between genic 5mC and

5hmC levels and expression. Grouping genes by Dnmt3a WT mOSN expression quartiles indicates that 5mC and 5hmC levels are largely maintained across genes (Figure S4E). However, TSS-flanking 5mC is decreased in Dnmt3a KO mOSNs at moderately expressed genes. In contrast, 5hmC levels are steadily reduced across the gene body (with a bias toward the 50 end) with increasing levels of expression (Figure S4E). To determine how the alteration of 5mC and 5hmC levels across genes correlates with changes in expression, we computed correlation values between 5mC or 5hmC KO/WT fold-change across genes and the corresponding expression fold-change (Figure 6G). After adjusting for multiple comparisons, we find that Dnmt3a-dependent 5mC levels surrounding the TSS are negatively correlated with expression changes, consistent with its repressive role. In contrast, Dnmt3a-dependent 5hmC levels across the gene body are positively correlated with expression changes in Dnmt3a KO mOSNs, whereas those present at the TSS are negatively correlated, consistent with a dual regulatory function (Colquitt et al., 2013; Wu et al., 2011). Ordering average 5mC levels at the TSS and 5hmC levels at the gene body by Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc. 831

Neuron Dnmt3a Regulates Odorant-Induced Transcription

A

B

C

D

G

E

H

I

F

J

Figure 6. Loss of Dnmt3a Results in the Bidirectional Disruption of Transcription in mOSNs (A and B) Fold-change expression in Dnmt3a Het (A) and KO (B) mOSNs relative to Dnmt3a WT mOSNs. x axis: average log2(RPKM + 1) across Dnmt3a WT and Het (A) or WT and KO (B) for each gene. y axis: log2 fold-change for each gene. (C) RNA log2(RPKM + 1) values in Dnmt3a WT, Het, and KO mOSNs normalized to WT. Displayed are values from genes that have an absolute log2 foldchange > 1. (D) RT-qPCR using RNA from Dnmt3a WT, Het, and KO mOSNs. Values are normalized to Tubb5 expression then to WT for each gene. *p < 0.05, **p < 0.01, ***p < 0.001, one-way ANOVA, n = 3 each genotype. (E) Comparison of fold-change expression estimates by RT-qPCR (y axis) and RNA-seq (x axis). Error bars are SEM for RT-qPCR and the range of two replicates for RNA-seq. (F) ISH of Dnmt3a WT and KO coronal MOE sections using probe against Omp. (G) Spearman rho correlations between Dnmt3a KO/WT fold-change 5mC or 5hmC at each gene position and fold-change gene expression. Dashed blue lines indicate rho at which the FDR-adjusted p value < 0.01. Shaded regions indicate positioned used for analysis in (H). TES, transcription end site. (H) Relationships between average TSS 5mC and gene body 5hmC levels and KO/WT fold-change expression. Lines are loess-smoothed averages, and error is 95% confidence interval. (I and J) Dnmt3a WT, Het, and KO 5mC and 5hmC levels across an upregulated gene (Ugt2a1, I) and a downregulated gene (Kirrel, J). Shades regions correspond to those in (G). See also Figure S4.

832 Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc.

Neuron Dnmt3a Regulates Odorant-Induced Transcription

KO/WT expression fold-change demonstrates striking differences between the two modifications (Figure 6H): 5mC levels are mildly reduced in Dnmt3a KO mOSNs at the TSSs of upregulated genes, whereas 5hmC levels display Dnmt3a dosagesensitive depletions within downregulated genes. Two representative upregulated (Utg2a1) and downregulated (Kirrel) genes (RT-qPCR expression levels found in Figure 6D) provide examples of this relationship (Figures 6I and 6J): 5mC levels are reduced at the TSS of Utg2a1 and 5hmC levels are reduced within the gene body of Kirrel. Hierarchical clustering of Dnmt3a-dependent block gene expression demonstrates a similar bidirectional disruption of transcription (Figure S4F). The two main fold-change expression clusters (upregulated, cluster 1; downregulated, cluster 2) are not significantly enriched for particular gene ontology terms with respect to each other, but, as with the full set of genes, lowly expressed genes tend to exhibit increased expression, whereas active genes tend to display reduced expression in Dnmt3a KO mOSNs (Figure S4G). Moreover, in agreement with the association of 5mC with repressed and 5hmC with active transcriptional states, 5mC is enriched in regions associated with cluster 1 (repressed and upregulated), whereas 5hmC is enriched in regions associated with cluster 2 (active and downregulated; Figure S4G). A previous study observed increased levels of the repressive histone modification H3K27me3 at the TSSs of genes downregulated in neural stem cells in the absence of Dnmt3a (Wu et al., 2010). This increase was argued to result from the inappropriate spread of Polycomb complex activity, which generates H3K27me3, onto TSSs that are normally protected by Dnmt3amediated 5mC. To determine if a similar increase occurs within the olfactory epithelium, we performed H3K27me3 ChIP-seq using Dnmt3a KO MOEs. In contrast to the previous report, we do not observe an increase of H3K27me3 levels at genes with reduced expression in KO mOSNs and instead detect a modest depletion of H3K27me3 at the TSSs of genes with increased expression (Figure S4H; median RPKM ± SD: WT = 0.37 ± 1.33, KO = 0.31 ± 1.09, p = 1 3 10
4, Wilcoxon two-sided). These expression defects led us to investigate transcriptional responses to odorant-induced neuronal activation in Dnmt3adeficient animals. In particular, we hypothesized that the mild transcriptional effects observed in Dnmt3a KO mOSNs would be exacerbated upon gene induction. Initially, we exposed Dnmt3a WT mice to a concentrated mixture of previously characterized odorants (1:1:1 amyl acetate:acetophenone:octanal [am:ac:oc]) (Mombaerts, 2004; Rubin and Katz, 1999) and assayed transcriptional levels in their MOEs relative to those of mice exposed to water (Figure 7A). We detected maximal induction of the immediate early genes Fos, Egr1, and Nr4a1 after 1 hr exposure (Figures 7B and S5A). The induction of Egr1 occurs broadly in the mOSN layer (Figure 7C). To globally assay odorant-induced transcription, we prepared two replicate RNA-seq libraries from MOEs of Dnmt3a WT, Het, and KO mice that had been exposed to 1 hr of am:ac:oc or water (Tables S1 and S6). Across genotypes, we identified 23 genes that were significantly upregulated and one gene that was downregulated (false discovery rate [FDR] < 0.05) in response to odorant exposure (Table S7). Among these genes was the transcription factor

Nr4a2 that strikingly exhibits both depletions of 5mC and 5hmC levels in its gene body and the reduction of odorant-induced expression in Dnmt3a KOs (Figure 7D). In contrast, Pcdh10, which is induced to WT levels in the Dnmt3a KO (Figures 7D and S5B), has WT levels of DNA modifications in the KO OSNs (Figure 7D). On average, the fold transcriptional response to odorant of the set of 23 genes is significantly reduced in Dnmt3a KOs (median log2(Het/WT) RPKM ± SD:
0.93 ± 1.30, p = 1 3 10
3, two-sided Wilcoxon; Figure S5C), and their gene bodies are significantly depleted of both 5mC (p = 4 3 10
5) and 5hmC (p = 2 3 10
4) (Figure S5D). Interestingly, although odorant-induced transcription levels are reduced in Dnmt3a Hets (Figure 7E; median log2(Het/WT) RPKM ± SD:
0.22 ± 0.29, p = 3 3 10
3, paired two-sided Wilcoxon) and KOs (
0.94 ± 1.23, p = 7 3 10
4), basal levels in both genotypes remain unchanged (Het:
0.08 ± 0.46, p = 0.16; KO: 0.03 ± 0.39, p = 0.84). This association between reduced 5mC, 5hmC, and reduced gene activation generalizes across odorant-induced genes (Figures 7F and 7G). Dnmt3a KO/WT levels of 5mC and (to a lesser extent) 5hmC in the bodies of odorant-induced genes are positively correlated with KO/WT fold-change odorant-induced RNA levels (5mC: Pearson’s R = 0.45, p = 0.02, permutation test; 5hmC: R = 0.16, p = 0.21). Strikingly, unlike the mild global reduction of expressed genes in Dnmt3a KO mOSNs (log2(KO/ WT RPKM) of
0.3), we find a median
0.9 log2 reduction of odorant-induced expression in Dnmt3a KO relative to WT (p = 2 3 10
3, Wilcoxon two-sided), indicating that activity-dependent responses are particularly sensitive to the loss of the enzyme (Figure 7H). Several lines of evidence argue against these effects being a result of a general disruption of odorant detection. The loss of Dnmt3a does not affect the protein levels of several markers of mOSN identity, including the major mOSN signaling components (Adcy3 and Gnal), the mOSN axon guidance molecule Cntn4, and the glutamate vesicular transporter Slc17a6 (Figure S5E). Similarly, glomerulus formation of the most highly expressed olfactory receptor, Olfr1507, is unaffected by the absence of Dnmt3a (Figure S5F). Finally, phosphorylation of the S6 ribosomal subunit, a marker for neuronal activity (Knight et al., 2012), is induced to similar levels in both Dnmt3a WT and KO animals exposed to am:ac:oc for 1 hr (Figure S5G). DISCUSSION The patterning of 5mC and 5hmC during neurodevelopment is dynamic, and it is likely that these changes inform the development and function of neurons. Our comprehensive analysis of how the lack of Dnmt3a, a de novo DNA methyltransferase strongly implicated in neurodevelopment, influences the epigenetic and transcriptional states of mOSNs provides insight into how DNA modifications contribute to defining the neuronal state. Although de novo DNA methyltransferase activity is completely absent in Dnmt3a KO mOSNs, the vast majority of 5mC and 5hmC patterning is unaffected by the loss of the enzyme. Instead, Dnmt3a is required to pattern defined components of the 5mC and 5hmC landscapes of mOSNs that are strongly associated with accessible regions of the genome. Recent work described a global buildup of methylated (m)CH within Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc. 833

Neuron Dnmt3a Regulates Odorant-Induced Transcription

A

B

C

D

E

F

G

H

Figure 7. Loss of Dnmt3a Disrupts DNA Modification States and Inducibilities of Odorant-Sensitive Genes (A) Schematic of odorant exposure experiments. (B) RT-qPCR of immediate early gene Egr1 expression in MOEs from 3-week-old WT mice that had been exposed to either water or 1:1:1 am:ac:oc for 15, 30, 60, or 120 min or exposed for 60 min then allowed to recover in a new cage for 60 min (60+/
). Expression levels are expressed as the fraction of Actb expression ± SEM. The p values are from Student’s t test. (C) ISHs against Egr1 on coronal MOE sections from Dnmt3a WT mice exposed to either water or am:ac:oc for 1 hr. (D) 5hmC, 5mC, and RNA tracks of two representative odorant-induced genes, Nr4a2 and Pcdh10. One-hour exposure to am:ac:oc (odor +) or water (odor
). (E) WT-normalized RPKM values of 23 odorant-induced genes in either water or am:ac:oc conditions. Blue bars intersect with the median of each genotype. Values are the averages across two biological replicates. (legend continued on next page)

834 Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc.

Neuron Dnmt3a Regulates Odorant-Induced Transcription

neurons in the forebrain (Lister et al., 2013) and dentate gyrus (Guo et al., 2014) that parallels Dnmt3a expression. Moreover, in the dentate gyrus methylation in this dinucleotide context but not the CG context is partially dependent on Dnmt3a. The apparent stability of much of mOSN 5mC patterning in the absence of Dnmt3a could be due to relatively low levels of mCH within this peripheral neuronal population or to the poor sensitivity of the antibody-based approach used here in detecting low levels of modification. Moreover, we cannot exclude that the antibody we used in these experiments has significantly lower affinity for mCH versus methylated CG. Genome-wide bisulfite sequencing will be needed to resolve this issue. In addition, it will be important to assay DNA modification patterning in GBCs (OSN progenitors) in the absence of Dnmt3a to specifically identify sites of Dnmt3a activity during differentiation. Furthermore, because the Dnmt3a mutant allele is constitutive, we cannot rule out possible compensatory effects of Dnmt3b, the other mammalian de novo DNA methyltransferase, during earlier developmental time points, including in HBCs, the multipotent stem cell population of MOE. However, both Dnmt1 and Dnmt3b expression is largely unaltered in the MOEs of Dnmt3a KO mice, making it unlikely that the persistence of 5mC or 5hmC in this tissue is due to atypical activities of the other methyltransferases. Nonetheless, the constitutive mutation makes it difficult to formally rule out the influence of other factors, such as the alteration of epigenetic patterning during early development, that may have indirectly produced the observations presented here. Interestingly, Dnmt3a is required to fully generate the elevated 5hmC levels established during mOSN development. This dependency suggests that de novo 5mC in neurons is targeted by the Tet oxidation machinery. In support of this model, we also find that Dnmt3a-dependent 5mC and 5hmC patterning is positively correlated with Tet3-sensitive 5hmC patterning. In neurons, 5hmC occurs nearly exclusively on cytosines in CG dinucleotides, whereas de novo 5mC occurs broadly in multiple contexts (Lister et al., 2013). This suggests that Dnmt3a, which has been shown to provide both CG and CH methylation in nonneural systems (Meissner et al., 2005; Ramsahoye et al., 2000), generates a 5mC platform on which neuronal Tet activity can provide elevated levels of 5hmC. Alternatively, 5mC oxidation may occur broadly across cytosine contexts but with less processivity in CG contexts, resulting in higher levels of 5hmC, than on those in CH contexts. Recent structural data and in vitro enzymatic assays of human TET2 indicate that the Tet enzymatic domain specifically oxidizes 5mC within the CG context, in support of the former model (Hu et al., 2013). However, these assays have not been extended to mouse Tet3, which is the dominant Tet family member in mOSNs (Colquitt et al., 2013), or to Tet1. We describe a set of broad domains that contain Dnmt3adependent 5mC and 5hmC patterning, which we term Dnmt3a-dependent blocks. These regions are strongly associ-

ated with key genomic loci, including genes and regulatory elements associated with neuronal development. In particular, approximately 200 transcription factors are embedded within or near these Dnmt3a-dependent blocks. The spatial correspondence between accessibility and Dnmt3a-mediated patterning suggests either that Dnmt3a is only able to modify open genomic regions or that Dnmt3a itself establishes zones of chromatin accessibility. This patterning is reminiscent of recently characterized ‘‘super-enhancers,’’ which are defined by the extensive occupancy of transcription factors and structural proteins implicated in mediating enhancer-promoter contacts as well as the enrichment of active histone modifications (Hnisz et al., 2013; Whyte et al., 2013). Although the function of these regions is unclear, superenhancers are strongly associated with transcription factors involved in the specification of cell identity. Moreover, much like Dnmt3a-dependent blocks, superenhancers often span transcription factor loci, thus establishing large, permissive regulatory domains across these regions. These likenesses suggest that Dnmt3a-mediated DNA modification patterning is an additional feature of superenhancers; however, what role this patterning plays in defining the function of these enhancers remains to be determined. The finding that 5mC and 5hmC levels depend on Dnmt3a dosage, evidenced by intermediate 5mC and 5hmC levels in heterozygotes, argues that either Dnmt3a protein level is limiting or that Dnmt3a-modified sites are engaged in a methylation-demethylation cycle, rendering them vulnerable to fluctuations in Dnmt3a activity. Interestingly, we find that these DNA modification effects are propagated to intermediate levels of gene expression in Dnmt3a Het mOSNs. This is surprising considering standard models of DNA methylation function in which the methylation of key CGs alters the affinity of protein trans factors. Instead, these data support a model in which the density of 5mC and its derivative 5hmC—and not the modification of key positions—is the functionally relevant output of Dnmt3a-mediated patterning. Indeed, the mosaicism of 5mC or 5hmC in Dnmt3a heterozygotes, the distributive activity of the enzyme (Gowher and Jeltsch, 2001), as well as its ability to methylate cytosines in multiple contexts suggest that Dnmt3a stochastically methylates targeted sites. This behavior may permit the fine-tuning of gene expression levels through intermediate DNA modification patterning. In keeping with such fine regulation, the absence of Dnmt3a in mOSNs causes bidirectional alterations to gene expression profiles, through the upregulation of repressed genes and downregulation of expressed genes. Traditionally associated with transcriptional repression, DNA methylation is increasingly appreciated as a bifunctional modification depending on its genomic context and interpretation. In particular, gene body methylation in Arabidopsis thaliana (Lister et al., 2008; Zilberman et al., 2008), the mammalian X chromosome (Hellman and Chess, 2007), and embryonic stem cells (Lister et al., 2009) is

(F and G) Gene-body 5mC (F) and 5hmC (G) loss is correlated with reduced odorant-induced expression. x axis: average log2(KO/WT RPKM) 5mC or 5hmC across gene-bodies. y axis: log2(KO/WT) RPKM of induced RNA expression of odorant-induced genes. (H) Comparison of KO versus WT fold-change for total mOSN-expressed genes (as defined in text) and odorant-induced genes. The p value is from Wilcoxon two-sided. See also Figure S5.

Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc. 835

Neuron Dnmt3a Regulates Odorant-Induced Transcription

positively correlated with expression. Similarly, although 5hmC at the TSS is negatively correlated with transcription in neurons, 5hmC levels flanking TSSs and within gene bodies generally positively correlate, and to a greater degree than 5mC, with transcriptional output (Colquitt et al., 2013; Hahn et al., 2013; Melle´n et al., 2012; Song et al., 2011; Szulwach et al., 2011). We observed in previous work that overexpressing Tet3 in mOSNs resulted in either the down- or upregulation of transcription depending on whether 5hmC was decreased or increased, respectively, in gene bodies (Colquitt et al., 2013). Accordingly, reduced gene body 5hmC in Dnmt3a KO mOSNs correlates with mildly reduced transcriptional output, consistent with 5hmC playing a facilitating role in transcription. Here, despite the reduction of 5hmC levels in the absence of de novo DNA methyltransferase activity, the bulk of 5hmC patterning in mOSNs is unaltered, indicating that most of this patterning is stable in this neuronal population. This strongly suggests that neuronal 5hmC is not exclusively an intermediate in Tet-mediated demethylation but is itself stably maintained in neurons. Recent characterizations of proteins that bind 5hmC, including MBD3 (Yildirim et al., 2011), methyl-CpG binding protein 2 (MeCP2) (Melle´n et al., 2012), and a diverse array of other factors (Spruijt et al., 2013), suggest that the modified base is actively interpreted and may directly influence gene regulation. Finally, our finding that odorant-induced transcriptional responses require Dnmt3a, combined with the observation that substantial de novo 5mC occurs in neurons around the time of synaptogenesis (Lister et al., 2013), argues that de novo neuronal DNA modification patterns are critical to regulate activity-sensitive transcriptional programs. This connection is supported by experiments indicating that MeCP2, which binds both 5mC and 5hmC in neurons (Melle´n et al., 2012), is phosphorylated by neuronal activity (Cohen et al., 2011; Zhou et al., 2006), indicating that a neuron’s activity state is communicated to one of the primary interpreters of the epigenome. Furthermore, DNA modification patterns associated with activity-dependent genes are sensitive to neuronal activity in the hippocampus (Guo et al., 2011a, 2011b). Interestingly, despite significant reductions of 5mC and 5hmC levels on activity-regulated genes in the absence of odorant stimulation, basal transcription levels are unaffected by Dnmt3a deletion. Transcriptional perturbations of these genes only become apparent with gene induction. This relationship suggests that a Dnmt3a-mediated DNA modification signature on these activity-regulated genes is required for enhanced transcription upon neuronal stimulation. Taking into account recent observations suggesting unexpected plasticity in the peripheral olfactory system (Abraham et al., 2014; Kass et al., 2013), our data may provide insight into the molecular underpinnings of activity-induced transcription in olfactory neurons. In effect, Dnmt3a may establish an epigenetic landscape that prepares the neuronal genome for future transcriptional responses, supporting a crucial role for DNA modifications in neuronal plasticity. EXPERIMENTAL PROCEDURES Animal Care and Use Mice were treated in compliance with the rules and regulations of Institutional Animal Care and Use Committee under protocol approval AN099395-01.

836 Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc.

Dnmt3a constitutive null mice (Okano et al., 1999) were obtained from MGI (2182412). FACS Mature OSNs were isolated from the MOEs of Omp-ires-GFP knockin animals (Shykind et al., 2004). Single-cell suspensions from 1 to 10 MOEs were made using the papain dissociation system (Worthington Biochemical): MOEs were dissected into Earl’s buffered saline solution, minced in papain dissociation solution (Worthington Biochem), and incubated at 37 C for 45 min. Cells were washed once in 1:10 inhibitor solution and once in PBS. For a given assay (DIP-seq or RNA-seq), 2 3 105 to 1 3 106 cells were used. Population purity was assayed by resorting a small fraction of the sorted cells. Purities of 90%–95% were routinely achieved. DIP-Seq Genomic DNA was isolated using the DNeasy Blood and Tissue Kit (QIAGEN) from the cell populations described in the main text, prepared for paired-end Illumina sequencing using 75–200 ng of gDNA and standard protocols, and immunoprecipitated as described in Weber et al. (2005), with minor modifications. See Supplemental Experimental Procedures for details. ChIP-Seq Nuclei were isolated from the MOEs of 3-week-old Dnmt3a WT and KO mice, and native ChIP was performed as described previously (Magklara et al., 2011) using anti-H3K27ac (Millipore, cma309), anti-H3K4me1 (Abcam, ab8895), and anti-H3K27me3 (Millipore, 07-449) antibodies. For library preparation, ChIP DNA was sonicated for 120–180 s on a Covaris S220 focusedultrasonicator and prepared for sequencing using the Ovation Ultralow Library Kit (NuGEN). Odorant Exposure Mice were divided into two clean cages containing no bedding and placed into separate dark cabinets. After 10 min in these cages, a 4 cm 3 4 cm piece of Whatman paper with either 80 ml of water or a 1:1:1 mixture of pure am:ac:oc (all purchased from Sigma-Aldrich) was placed at one end of the cage. Details for the timing of odorant exposure and associated RNA analysis are given in the Supplemental Experimental Procedures. ACCESSION NUMBERS The Gene Expression Omnibus accession number for all data sets reported in this paper is GSE52464. SUPPLEMENTAL INFORMATION Supplemental Information includes Supplemental Experimental Procedures, Figures S1–S5, and Tables S1–S7 and can be found with this article online at http://dx.doi.org/10.1016/j.neuron.2014.07.013. AUTHOR CONTRIBUTIONS B.M.C. designed the research project, performed most experiments, analyzed the data, and wrote the manuscript. E.M.-P. performed the DNase-seq and several ChIP-seq assays. R.D. performed the bisulfite sequencing and analysis. S.L. designed the research project and wrote the manuscript. ACKNOWLEDGMENTS Sequencing data are deposited in GEO under accession GSE52464. This work was supported by the NIH (DP2 OD006667, R01MH091661), NSF (Graduate Research Fellowship), and the Rett Syndrome Trust. Accepted: July 8, 2014 Published: August 7, 2014

Neuron Dnmt3a Regulates Odorant-Induced Transcription

REFERENCES Abraham, N.M., Vincis, R., Lagier, S., Rodriguez, I., and Carleton, A. (2014). Long term functional plasticity of sensory inputs mediated by olfactory learning. eLife 3, e02109. Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11, R106. Caggiano, M., Kauer, J.S., and Hunter, D.D. (1994). Globose basal cells are neuronal progenitors in the olfactory epithelium: a lineage analysis using a replication-incompetent retrovirus. Neuron 13, 339–352. Cohen, S., Gabel, H.W., Hemberg, M., Hutchinson, A.N., Sadacca, L.A., Ebert, D.H., Harmin, D.A., Greenberg, R.S., Verdine, V.K., Zhou, Z., et al. (2011). Genome-wide activity-dependent MeCP2 phosphorylation regulates nervous system development and function. Neuron 72, 72–85. Colquitt, B.M., Allen, W.E., Barnea, G., and Lomvardas, S. (2013). Alteration of genic 5-hydroxymethylcytosine patterning in olfactory neurons correlates with changes in gene expression and cell identity. Proc. Natl. Acad. Sci. USA 110, 14682–14687. Creyghton, M.P., Cheng, A.W., Welstead, G.G., Kooistra, T., Carey, B.W., Steine, E.J., Hanna, J., Lodato, M.A., Frampton, G.M., Sharp, P.A., et al. (2010). Histone H3K27ac separates active from poised enhancers and predicts developmental state. Proc. Natl. Acad. Sci. USA 107, 21931–21936. Feng, J., Chang, H., Li, E., and Fan, G. (2005). Dynamic expression of de novo DNA methyltransferases Dnmt3a and Dnmt3b in the central nervous system. J. Neurosci. Res. 79, 734–746. Feng, J., Zhou, Y., Campbell, S.L., Le, T., Li, E., Sweatt, J.D., Silva, A.J., and Fan, G. (2010). Dnmt1 and Dnmt3a maintain DNA methylation and regulate synaptic function in adult forebrain neurons. Nat. Neurosci. 13, 423–430. Globisch, D., Mu¨nzel, M., Mu¨ller, M., Michalakis, S., Wagner, M., Koch, S., Bru¨ckl, T., Biel, M., and Carell, T. (2010). Tissue distribution of 5-hydroxymethylcytosine and search for active demethylation intermediates. PLoS ONE 5, e15367–e15367. Gowher, H., and Jeltsch, A. (2001). Enzymatic properties of recombinant Dnmt3a DNA methyltransferase from mouse: the enzyme modifies DNA in a non-processive manner and also methylates non-CpG [correction of nonCpA] sites. J. Mol. Biol. 309, 1201–1208. Guo, J.U., Ma, D.K., Mo, H., Ball, M.P., Jang, M.-H., Bonaguidi, M.A., Balazer, J.A., Eaves, H.L., Xie, B., Ford, E., et al. (2011a). Neuronal activity modifies the DNA methylation landscape in the adult brain. Nat. Neurosci. 14, 1345–1351. Guo, J.U., Su, Y., Zhong, C., Ming, G.L., and Song, H. (2011b). Hydroxylation of 5-methylcytosine by TET1 promotes active DNA demethylation in the adult brain. Cell 145, 423–434. Guo, J.U., Su, Y., Shin, J.H., Shin, J., Li, H., Xie, B., Zhong, C., Hu, S., Le, T., Fan, G., et al. (2014). Distribution, recognition and regulation of non-CpG methylation in the adult mammalian brain. Nat. Neurosci. 17, 215–222. Hahn, M.A., Qiu, R., Wu, X., Li, A.X., Zhang, H., Wang, J., Jui, J., Jin, S.-G., Jiang, Y., Pfeifer, G.P., and Lu, Q. (2013). Dynamics of 5-hydroxymethylcytosine and chromatin marks in Mammalian neurogenesis. Cell Rep. 3, 291–300. Hellman, A., and Chess, A. (2007). Gene body-specific methylation on the active X chromosome. Science 315, 1141–1143. Hirota, J., and Mombaerts, P. (2004). The LIM-homeodomain protein Lhx2 is required for complete development of mouse olfactory sensory neurons. Proc. Natl. Acad. Sci. USA 101, 8751–8755. Hnisz, D., Abraham, B.J., Lee, T.I., Lau, A., Saint-Andre´, V., Sigova, A.A., Hoke, H.A., and Young, R.A. (2013). Super-enhancers in the control of cell identity and disease. Cell 155, 934–947. Hu, L., Li, Z., Cheng, J., Rao, Q., Gong, W., Liu, M., Shi, Y.G., Zhu, J., Wang, P., and Xu, Y. (2013). Crystal structure of TET2-DNA complex: insight into TETmediated 5mC oxidation. Cell 155, 1545–1555. Huang, Y., Pastor, W.A., Shen, Y., Tahiliani, M., Liu, D.R., and Rao, A. (2010). The behaviour of 5-hydroxymethylcytosine in bisulfite sequencing. PLoS ONE 5, e8888.

Ito, S., D’Alessio, A.C., Taranova, O.V., Hong, K., Sowers, L.C., and Zhang, Y. (2010). Role of Tet proteins in 5mC to 5hmC conversion, ES-cell self-renewal and inner cell mass specification. Nature 466, 1129–1133. Kass, M.D., Rosenthal, M.C., Pottackal, J., and McGann, J.P. (2013). Fear learning enhances neural responses to threat-predictive sensory stimuli. Science 342, 1389–1392. Knight, Z.A., Tan, K., Birsoy, K., Schmidt, S., Garrison, J.L., Wysocki, R.W., Emiliano, A., Ekstrand, M.I., and Friedman, J.M. (2012). Molecular profiling of activated neurons by phosphorylated ribosome capture. Cell 151, 1126– 1137. Kriaucionis, S., and Heintz, N. (2009). The nuclear DNA base 5-hydroxymethylcytosine is present in Purkinje neurons and the brain. Science 324, 929–930. Leung, C.T., Coulombe, P.A., and Reed, R.R. (2007). Contribution of olfactory neural stem cells to tissue maintenance and regeneration. Nat. Neurosci. 10, 720–726. Lister, R., O’Malley, R.C., Tonti-Filippini, J., Gregory, B.D., Berry, C.C., Millar, A.H., and Ecker, J.R. (2008). Highly integrated single-base resolution maps of the epigenome in Arabidopsis. Cell 133, 523–536. Lister, R., Pelizzola, M., Dowen, R.H., Hawkins, R.D., Hon, G., Tonti-Filippini, J., Nery, J.R., Lee, L., Ye, Z., Ngo, Q.-M., et al. (2009). Human DNA methylomes at base resolution show widespread epigenomic differences. Nature 462, 315–322. Lister, R., Mukamel, E.A., Nery, J.R., Urich, M., Puddifoot, C.A., Johnson, N.D., Lucero, J., Huang, Y., Dwork, A.J., Schultz, M.D., et al. (2013). Global epigenomic reconfiguration during mammalian brain development. Science 341, 1237905. MacDonald, J.L., Gin, C.S.Y., and Roskams, A.J. (2005). Stage-specific induction of DNA methyltransferases in olfactory receptor neuron development. Dev. Biol. 288, 461–473. Magklara, A., Yen, A., Colquitt, B.M., Clowney, E.J., Allen, W., MarkenscoffPapadimitriou, E., Evans, Z.A., Kheradpour, P., Mountoufaris, G., Carey, C., et al. (2011). An epigenetic signature for monoallelic olfactory receptor expression. Cell 145, 555–570. Meissner, A., Gnirke, A., Bell, G.W., Ramsahoye, B., Lander, E.S., and Jaenisch, R. (2005). Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis. Nucleic Acids Res. 33, 5868–5877. Melle´n, M., Ayata, P., Dewell, S., Kriaucionis, S., and Heintz, N. (2012). MeCP2 binds to 5hmC enriched within active genes and accessible chromatin in the nervous system. Cell 151, 1417–1430. Mombaerts, P. (2004). Genes and ligands for odorant, vomeronasal and taste receptors. Nat. Rev. Neurosci. 5, 263–278. Nguyen, S., Meletis, K., Fu, D., Jhaveri, S., and Jaenisch, R. (2007). Ablation of de novo DNA methyltransferase Dnmt3a in the nervous system leads to neuromuscular defects and shortened lifespan. Dev. Dyn. 236, 1663–1676. Nicolay, D.J., Doucette, J.R., and Nazarali, A.J. (2006). Transcriptional regulation of neurogenesis in the olfactory epithelium. Cell. Mol. Neurobiol. 26, 803–821. Okano, M., Bell, D.W., Haber, D.A., and Li, E. (1999). DNA methyltransferases Dnmt3a and Dnmt3b are essential for de novo methylation and mammalian development. Cell 99, 247–257. Ramsahoye, B.H., Biniszkiewicz, D., Lyko, F., Clark, V., Bird, A.P., and Jaenisch, R. (2000). Non-CpG methylation is prevalent in embryonic stem cells and may be mediated by DNA methyltransferase 3a. Proc. Natl. Acad. Sci. USA 97, 5237–5242. Rubin, B.D., and Katz, L.C. (1999). Optical imaging of odorant representations in the mammalian olfactory bulb. Neuron 23, 499–511. Ruzov, A., Tsenkina, Y., Serio, A., Dudnakova, T., Fletcher, J., Bai, Y., Chebotareva, T., Pells, S., Hannoun, Z., Sullivan, G., et al. (2011). Lineage-specific distribution of high levels of genomic 5-hydroxymethylcytosine in mammalian development. Cell Res. 21, 1332–1342.

Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc. 837

Neuron Dnmt3a Regulates Odorant-Induced Transcription

Shykind, B.M., Rohani, S.C., O’Donnell, S., Nemes, A., Mendelsohn, M., Sun, Y., Axel, R., and Barnea, G. (2004). Gene switching and the stability of odorant receptor gene choice. Cell 117, 801–815.

Williams, K., Christensen, J., Pedersen, M.T., Johansen, J.V., Cloos, P.A.C., Rappsilber, J., and Helin, K. (2011). TET1 and hydroxymethylcytosine in transcription and DNA methylation fidelity. Nature 473, 343–348.

Song, C.-X., Szulwach, K.E., Fu, Y., Dai, Q., Yi, C., Li, X., Li, Y., Chen, C.-H., Zhang, W., Jian, X., et al. (2011). Selective chemical labeling reveals the genome-wide distribution of 5-hydroxymethylcytosine. Nat. Biotechnol. 29, 68–72.

Wu, H., and Zhang, Y. (2011). Mechanisms and functions of Tet protein-mediated 5-methylcytosine oxidation. Genes Dev. 25, 2436–2452.

Spruijt, C.G., Gnerlich, F., Smits, A.H., Pfaffeneder, T., Jansen, P.W.T.C., Bauer, C., Mu¨nzel, M., Wagner, M., Mu¨ller, M., Khan, F., et al. (2013). Dynamic readers for 5-(hydroxy)methylcytosine and its oxidized derivatives. Cell 152, 1146–1159. Szulwach, K.E., Li, X., Li, Y., Song, C.-X., Wu, H., Dai, Q., Irier, H., Upadhyay, A.K., Gearing, M., Levey, A.I., et al. (2011). 5-hmC-mediated epigenetic dynamics during postnatal neurodevelopment and aging. Nat. Neurosci. 14, 1607–1616. Visel, A., Minovitsky, S., Dubchak, I., and Pennacchio, L.A. (2007). VISTA Enhancer Browser—a database of tissue-specific human enhancers. Nucleic Acids Res. 35, D88–D92. Watanabe, D., Uchiyama, K., and Hanaoka, K. (2006). Transition of mouse de novo methyltransferases expression from Dnmt3b to Dnmt3a during neural progenitor cell development. Neuroscience 142, 727–737. Weber, M., Davies, J.J., Wittig, D., Oakeley, E.J., Haase, M., Lam, W.L., and Schu¨beler, D. (2005). Chromosome-wide and promoter-specific analyses identify sites of differential DNA methylation in normal and transformed human cells. Nat. Genet. 37, 853–862. Whyte, W.A., Orlando, D.A., Hnisz, D., Abraham, B.J., Lin, C.Y., Kagey, M.H., Rahl, P.B., Lee, T.I., and Young, R.A. (2013). Master transcription factors and mediator establish super-enhancers at key cell identity genes. Cell 153, 307–319.

838 Neuron 83, 823–838, August 20, 2014 ª2014 Elsevier Inc.

Wu, H., Coskun, V., Tao, J., Xie, W., Ge, W., Yoshikawa, K., Li, E., Zhang, Y., and Sun, Y.E. (2010). Dnmt3a-dependent nonpromoter DNA methylation facilitates transcription of neurogenic genes. Science 329, 444–448. Wu, H., D’Alessio, A.C., Ito, S., Wang, Z., Cui, K., Zhao, K., Sun, Y.E., and Zhang, Y. (2011). Genome-wide analysis of 5-hydroxymethylcytosine distribution reveals its dual function in transcriptional regulation in mouse embryonic stem cells. Genes Dev. 25, 679–684. Yildirim, O., Li, R., Hung, J.-H., Chen, P.B., Dong, X., Ee, L.-S., Weng, Z., Rando, O.J., and Fazzio, T.G. (2011). Mbd3/NURD complex regulates expression of 5-hydroxymethylcytosine marked genes in embryonic stem cells. Cell 147, 1498–1510. Zerucha, T., Stu¨hmer, T., Hatch, G., Park, B.K., Long, Q., Yu, G., Gambarotta, A., Schultz, J.R., Rubenstein, J.L.R., and Ekker, M. (2000). A highly conserved enhancer in the Dlx5/Dlx6 intergenic region is the site of cross-regulatory interactions between Dlx genes in the embryonic forebrain. J. Neurosci. 20, 709–721. Zhou, Z., Hong, E.J., Cohen, S., Zhao, W.-N., Ho, H.-Y.H., Schmidt, L., Chen, W.G., Lin, Y., Savner, E., Griffith, E.C., et al. (2006). Brain-specific phosphorylation of MeCP2 regulates activity-dependent Bdnf transcription, dendritic growth, and spine maturation. Neuron 52, 255–269. Zilberman, D., Coleman-Derr, D., Ballinger, T., and Henikoff, S. (2008). Histone H2A.Z and DNA methylation are mutually antagonistic chromatin marks. Nature 456, 125–129.