What this is
- Circadian rhythms regulate physiological processes in organisms through transcriptional feedback loops.
- This research investigates the relationship between circadian clock genes and histone modifications in Arabidopsis.
- It identifies how clock genes influence the expression of histone modifiers, impacting gene expression patterns.
Essence
- Circadian clock genes CCA1 and LHY regulate the expression of histone modifiers SDG2 and JMJ14, which in turn influence rhythmic histone modifications and gene expression in Arabidopsis.
Key takeaways
- CCA1 and LHY directly regulate the expression of the histone demethylase JMJ14 and indirectly regulate the methyltransferase SDG2. This regulation leads to rhythmic changes in histone modification levels, specifically , affecting target gene expression.
- Diurnal rhythms of and are associated with 20-30% of expressed genes, indicating a significant role of histone modifications in regulating gene expression throughout the day.
- Only 13% of morning-phased genes and 22% of evening-phased genes exhibit both and peaks, suggesting distinct functions of these histone modifications in diurnal gene expression.
Caveats
- The study does not clarify the exact mechanisms by which histone modifiers are recruited to their target genes, leaving gaps in understanding their functional dynamics.
- The overlap between CCA1 target genes and those with reduced levels is small, which may indicate that other factors influence histone modification patterns.
Definitions
- H3K4me3: Trimethylation of the fourth lysine residue on histone H3, associated with active transcription.
- H3K9ac: Acetylation of the ninth lysine residue on histone H3, often linked to gene activation.
AI simplified
Background
Circadian clocks anticipate diurnal cycles to coordinate physiological and developmental processes in all organisms from cyanobacteria, flies, to humans and plants [1–3], which are achieved through feedback transcriptional-translational regulatory modules [4–7]. In the transcriptional module, chromatin modifications regulate central circadian oscillators in mammals [8–10] and plants [7, 11–14]. In mouse liver, there are time-dependent patterns of transcription factor binding, RNA polymerase II recruitment, RNA abundance, and chromatin states, which account for 22% of mRNA cycling genes, suggesting a transcriptional mechanism for regulating the mammal circadian network [8].
In Arabidopsis, the central transcriptional loop involves two morning-phased Myb-like transcription factors, CIRCADIAN CLOCK ASSOCIATED 1 (CCA1) and LATE ELONGATED HYPOCOTYL (LHY) [15, 16], which repress the evening-phased regulator, timing of CAB expression 1 (TOC1), by directly binding to its promoter [17, 18]. In turn, TOC1 indirectly promotes the expression of CCA1 and LHY by inhibiting the repressive PSEUDO-RESPONSE REGULATOR (PRR) genes or directly represses expression of CCA1 and LHY by binding their promoters [19, 20]. Integrating other clock regulators including GI and ELF4–ELF3–LUX complex with the central CCA1/LHY/TOC1 loop provides a classic feedback regulatory mechanism for the circadian clock in Arabidopsis [21, 22], which includes a number of core molecular oscillators that constitute interlocking transcriptional-translational negative feedback loops (see reviews in [2, 4, 7, 23, 24]).
Using candidate gene approaches, several studies have indicated links between histone modifications and expression of circadian clock genes. For example, induction of the morning-phased TOC1 expression wave form is related to histone acetylation, which is repressed by CCA1 and coincident with histone deacetylase activities and repressive chromatin status at dawn [12]. Rhythmic expression of CCA1 and LHY and their reciprocal regulator TOC1 is associated with H3K4me3 and H3K9/14 ac levels in the translational start sites of these genes [11, 13], but correlated negatively with H3K36me2 levels [11]. Moreover, changes in diurnal transcript levels correlate with H3K9ac, H3K27ac, and H3S28p levels between end-of-day and end-of-night [25]. Inhibiting acetylation and H3K4me3 abolishes rhythmic expression of circadian oscillators, while blocking H3K4me3 leads to increased levels of clock-repressor binding, suggesting a transitional role for the H3K4me3 mark in modulating clock gene expression from activation to repression [14]. JMJD5, a histone demethylase in humans and plants, is co-regulated with evening-phased clock regulators and promotes expression of clock genes at dawn [26]. However, neither the basis for histone mark rhythms nor the relationship between histone modifications and circadian clock regulators is clearly defined. Here we report dynamic interactions between chromatin modification of circadian clock gene expression and diurnal regulation of histone modifications, which regulate diurnal gene expression networks in Arabidopsis.
Results
Diurnal regulation of histone methyltransferase and demethylase by clock genes
In contrast to SDG2, JMJ14 expression peaked at night (ZT40) and reached the valley after dawn (ZT28) in the WT and was significantly upregulated in the cca1 lhy mutant and downregulated in the CCA1-OX line under diurnal condition (Fig. 1b). The data suggest that CCA1/LHY exerts temporal effects on positive and negative correlation, respectively, with SDG2 and JMJ14 expression levels. Although both SDG2 and JMJ14 have one CCA1-binding site (CBS) in their promoter sequences, SDG2 and JMJ14 were not reported as target genes of CCA1 according to the published dataset of chromatin immunoprecipitation sequencing (ChIP-seq), which was performed using CCA1p::CCA1-GFP transgenic plants at ZT2 and ZT14 [39]. To determine whether CCA1 can bind promoters of SDG2 and JMJ14 at different time points, we performed ChIP-qPCR using antibodies against CCA1 at ZT0 when CCA expression peaks. ChIP-qPCR showed significant enrichment of CCA1 in the promoters of JMJ14 and TOC1 (normalized to UBQ10), but not in the promoter of SDG2 (Fig. 1c). Although we could not exclude a possibility of non-specific binding of CCA1-antibodies, this result indicates that CCA1/LHY may directly bind to the JMJ14 promoter and regulates its expression, as other members (JMJ30 and JMJD5) of the gene family that are co-regulated with evening-phased genes [26, 40], but may indirectly regulate SDG2 expression through other factors such as CCA1-mediated genes or in the CCA1 complex. Alternatively, regulation of SDG2 by CCA1 could occur at a different time of the day. In addition, other H3K4 methyltransferase ATX1 and H3K4 demethylase JMJ15 were also down- and upregulated, respectively, in the cca1 lhy mutant (Additional file 1: Figure S1c, d). The data suggest that CCA1/LHY can upregulate expression of H3K4me3 writers (SDG2 and ATX1) and downregulate H3K4me3 erasers (JMJ14 and JMJ15), which would result in the overall decreased distributions of H3K4me3 in the cca1 lhy mutant.
Using ChIP-seq analysis with antibodies against H3K4me3, we analyzed H3K4me3 distribution patterns in the rosette leaves of wild type (Ws) and cca1 lhy mutant plants at dawn (ZT0) under diurnal conditions, when CCA1/LHY expression level peaked [41]. Compared to the wild type, H3K4me3 levels near the transcription start sites of genic regions were damped in the cca1 lhy mutant (Fig. 1d). A total of 880 genes showed decreased levels of H3K4me3 peaks (one-way ANOVA, P < 0.05) (Additional file 2: Table S1), including two examples shown in Fig. 1e. A previous study has identified CCA1 target genes by ChIP-seq using 2-week-old seedlings grown under 12 h light/12 h dark for 12 days and then transferred to constant light for 2 days [39]. The overlapping fraction is small between the CCA1 target genes (5%, 78/1433) and those of reduced H3K4me3 levels in the cca1 lhy mutant (P > 0.5, hypergeometric test). The results indicate that CCA1 affects overall H3K4me3 accumulation levels but does not directly participate in the establishment of H3K4me3 marks for CCA1 target genes. However, different growth conditions and tissue stages in these two studies may also contribute to the small overlapping fraction between the CCA1 target genes and genes with reduced H3K4me3 levels in the cca1 lhy mutant.
The clock regulates expression of histone methyltransferase (SDG2) and demethylase (JMJ14) genes.,Relative expression levels (R.E.L.) of() and() in the wild type (WT, Ws),mutant (Ws) (upper panel), and WT (Col-0) and CCA1-OX (Col-0) (lower panel) plants under a diurnal cycle (16 h light/8 h dark). Black and white bars indicate dark (ZT16, 20, and 24) and light (ZT28, 32, 36, and 40), respectively.ChIP-qPCR showed the fold enrichment (-axis) of CCA1-binding fractions in the promoters ofandusing rosette tissues from 3-week-oldat ZT0. Thepromoter containing evening elements was used as the positive control.was used as a nonbinding control for normalization. Error bars in–indicate standard deviation of biological replicates (= 3) with an asterisk showing the statistical significance level (< 0.05, Studenttest).Binding profiles of H3K4me3 at the TSS ± 1.5 kb in the wild type andmutant with two biological replicates.Examples showing decreased binding of H3K4me3 to the genes(upper panel) and(lower panel) inmutant compared with the WT (Ws) a b a b c a c d e SDG2 JMJ14 cca1 lhy Y SDG2 JMJ14 A. thaliana TOC1 UBQ10 n P t cca1 lhy AT1G21270 AT3G05880 cca1 lhy
Disrupted histone modifications lead to altered expressions of clock genes
Disrupting histone modifications affects expression of key circadian clock regulators.ChIP-seq assays showing decreased binding of H3K4me3 to the genic regions of(top) and(bottom) in themutant compared with the wild type (WT). Green arrows next to the gene diagrams indicate transcriptional direction.Relative expression levels (R.E.L.) of(top) and(bottom) in themutant compared with the WT at ZT0, ZT12, and ZT24 under a diurnal cycle (16 h light/8 h dark).-axis was broken as a result from large variation of expression levels.ChIP-seq assays showing increased binding of H3K4me3 to the genic regions of(top) and(bottom) in themutant compared with the WT.R.E.L of(upper) and(lower) in themutant compared with the WT under a diurnal cycle (16 h light/8 h dark). Error bars inandindicate standard deviation of biological replicates (= 3) with an asterisk showing the statistical significance level (< 0.05, Studenttest) a b c d b d CCA1 LHY sdg2 CCA1 LHY sdg2 Y CCA1 LHY jmj14 CCA1 LHY jmj14 n P t
Dynamic diurnal rhythms of histone modifications in Arabidopsis
The phase distribution analysis indicated genome-wide binding rhythms of H3K4me3 with a major peak at ZT12 and a minor peak at ZT6, whereas the major peak of H3K9ac binding rhythms was delayed towards ZT15 and the minor peak was moved forward to ZT3 (Additional file 1: Figure S4a). This is consistent with total levels of H3K4me3 and H3K9ac as determined by the Western blot analysis (Additional file 1: Figure S4b). The ChIP-seq data further showed a diurnal rhythm of the total H3K4me3 levels with the peak at ZT12, whereas the total H3K9ac levels did not show obvious peaks during the day and night (Additional file 1: Figure S4b). The inconsistency between diurnal H3K9ac distribution and global H3K9ac content may be related to inconspicuous H3K9ac distribution variation in the time points examined in total H3K9ac levels (ZT0, 6, 12, 18, 24). Alternatively, non-rhythmic H3K9ac distribution may account for the majority of overall H3K9ac variation because of its role in other aspects of plant growth and development [47, 48].
The ChIP-seq data were further classified into the morning-phased (between ZT0 and ZT6) and evening-phased (between ZT9 and ZT18) peaks using k-medoids clustering [49] (Fig. 3a and Additional file 4: Table S3). A majority (86%, 14,029/16,391) of H3K4me3- or H3K9ac-enriched genes had both H3K4me3 and H3K9ac peaks, and only a small proportion of genes was associated with only one histone mark (Fig. 3b). This suggests a coordinated mechanism for H3K4me3 and H3K9ac to regulate gene expression [45]. However, among morning-phased genes, only 13% (313/2335) were associated with both morning-phased H3K4me3 and H3K9ac peaks. Among evening-phased genes, 696 genes (22%) were related to both evening-phased H3K4me3 and H3K9ac peaks (Fig. 3b). The proportion of genes associated with both histone marks in morning-phased (13%) or evening-phased (22%) peaks was significantly lower than that enriched with either H3K4me3 or H3K9ac peaks (P < 2e−100, hypergeometric test).
Gene ontology (GO) analysis suggested that the morning- and evening-phased H3K4me3 and H3K9ac marks were enriched in genes with different gene ontology (GO) groups (Fig. 3c). The morning-phased genes are enriched with photosynthetic activities and responses to a variety of stresses, while the evening-phased genes are enriched with DNA and RNA metabolism and protein ubiquitination and folding (Fig. 3c). These data implied a specificity for diurnal H3K4me3 and H3K9ac rhythms in the expression of genes involved in different biological networks. A total of 22% (958/4367) and 28% (1208/4367) genes with rhythmic H3K4me3 peaks in the sdg2 and jmj14 mutants, respectively, showed altered H3K4me3 levels (Fig. 3d), which was 1.5-fold higher than those genes (12%, 1252/10,503 and 17%, 1810/10,503) with non-rhythmic H3K4me3 peaks in the mutants, respectively, displaying altered H3K4me3 levels (P < 1e−10, Fisher’s exact test). The data suggest that the genes with rhythmic histone modifications are more sensitive to the loss of histone methyltransferase and demethylase activities than other genes.
To examine the correlation between oscillating histone modifications and gene expression, we analyzed expression patterns of genes using published data of expression profiling in Arabidopsis under diurnal (16 h light/8 h dark) and circadian (constant light) conditions [41]. For the genes with both diurnal and circadian expression patterns, 14% (612/4377) had both rhythmic H3K4me3 and H3K9ac marks (Additional file 4: Table S3). However, the fraction of both rhythmic H3K4me3 and H3K9ac marks was much lower for the genes with circadian but not diurnal expression patterns (4%, 150/4134) or for the genes with diurnal but not circadian expression patterns (4%, 139/3121). These results suggest different impact of histone-modifying enzymes in diurnal and/or circadian expression patterns. JMJ proteins are circadian-mediated probably through coregulation with the clock components [26, 40], while SGD2 expression is indirectly regulated by CCA1 (Fig. 1c and Additional file 1: Figure S1). At the genome-wide level, ~ 48% (2104/4367) of the genes with oscillating H3K4me3 and ~ 47% (1552/3276) of the genes with oscillating H3K9ac displayed rhythmic expression with distinct morning-phased and evening-phased patterns [41] (Fig. 3e), indicating a strong correlation between diurnal histone modifications and rhythmic gene expression.
Dynamics of diurnal histone modifications in.Heat-map view of H3K4me3 and H3K9ac binding signals every 3 h in a 24-h period (16 h light/8 h dark). The peaks are partitioned into the morning and evening clusters. The number of peaks in each cluster is indicated at the left side of the panel.Overlap of genes with H3K4me3 and H3K9ac peaks.Gene Ontology (GO) analysis of the genes with diurnal H3K4me3 or H3K9ac binding peaks.Overlap between rhythmic H3K4me3 peaks and peaks with altered H3K4me3 enrichment in theormutant compared with the WT. The asterisk indicates the statistical significance level for overlap with< 0.05.Phase distributions of the genes with rhythmic H3K4me3 or H3K9ac peaks in the morning phase (blue line) and evening phase (black line) Arabidopsis sdg2 jmj14 P a b c d e
Roles of histone modifications in core clock genes
A possibility for expression rhythms and chromatin modifications is that the core circadian clock genes may interact with chromatin-modifying enzymes to alter histone modifications and regulate target gene expression. To test this, we analyzed binding targets of CCA1 and TOC1 from the published ChIP-seq datasets [20, 39], and examined enrichment levels of H3K4me3 and H3K9ac in the target genes. As CCA1 and LHY often repress expression of evening-phased genes by binding their promoters [39], the majority of CCA1 targets showed rhythmic expression patterns with an evening phase. The results showed that approximately 33% and 20% of the CCA1 targets had evening-phased binding patterns for H3K4me3 and H3K9ac (P < 1e−50, hypergeometric test), respectively, while fewer than 7% of CCA1 targets exhibited morning-phased profiles for H3K4me3 and H3K9ac (Fig. 4e). This finding is consistent with the expression patterns of those CCA1 target genes in the evening phase (Additional file 1: Figure S7), although phased distributions of morning and evening genes could change under long- or short-day conditions [41, 50]. In contrast, 23% and 13% of TOC1 target genes correlated with H3K4me3 and H3K9ac marks in the morning phase (Fig. 4f) (P < 1e−5, hypergeometric test), while a smaller proportion of target genes correlated with those marks in the evening. This is consistent with TOC1 being a regulator of the morning-phased genes [19, 20] (Additional file 1: Figure S7). To further test the roles of SDG2 and JMJ14 in circadian-mediated gene expression, we randomly selected three CCA1 target genes and examined enrichment of H3K4me3 in genic region of these genes in wild type and sdg2 and jmj14 mutants at ZT0, 12, and 24. Enrichment of H3K4me3 in all three CCA1 target genes was significantly decreased in the sdg2 mutant at similar levels in the time points tested but increased in the jmj14 mutant in most time points (Additional file 1: Figure S8). The effect in a few time points was moderate in the jmj14, probably due to genetic redundancy, as observed for the clock genes (Fig. 2c and Additional file 1: Figure S3). Notably, the H3K4me3 levels in the sdg2 mutant are reduced at all time points tested (Additional file 1: Figure S8a), consistent with that SDG2 is not directly regulated by CCA1 (Fig. 1c). JMJ14, however, showed increased H3K4me3 levels with relative rhythms (Additional file 1: Figure S8b), consistent with a role for JMJ proteins in co-regulating with the circadian clock and its target genes [26, 32]. These data indicate that SDG2 and JMJ14 are involved in regulation of rhythmic deposition of H3K4me3 in CCA1 target genes.
Histone modifications mediate regulation of the clock genes.,Binding signals of H3K4me3 and H3K9ac at the genic regions of() and() at different time points. Each track represents the normalized ChIP-seq read coverage at a single time point. The-axis scales are 0–15 for H3K4me3 and 0–10 for H3K9ac. ZT time is indicated on the left.,Quantifications of the ChIP-seq signals at diurnal peaks located in() and(), respectively. Relative expression levels (R.E.L.) ofandare shown (blue line).was used as an internal control for gene expression.,Enrichment of() and() binding targets in the genes with rhythmic H3K4me3 or H3K9ac peaks. Error bars indicate the standard deviation of biological replicates (= 3) with an asterisk showing the significant enrichment (< 1e−5, hypergeometric test) a c a c b d b d e f e f CCA1 TOC1 Y CCA1 TOC1 CCA1 TOC1 ACT7 CCA1 TOC1 n P
Discussion and conclusions
A large proportion of histone modification sites exhibit rhythmic patterns, while different histone modifications (H3K4me3 and H3K9ac) have distinct phases for their target genes, which is consistent with the report in mammals [8]. The data collectively indicate a conserved correlation between histone modifications and circadian clock regulation in plants and mammals, as well as a specific distribution of different histone modifications in morning-phased and evening-phased genes. This could be related to distinct peaks of H3K4me3 and H3K9ac in the morning and evening phases, respectively, which correspond to expression of morning- and evening-phased genes (Fig. 3). We also noticed both H3K4me3 and H3K9ac showed lower levels at the night (ZT18-24) than during the day (Fig. 3a). One possibility is that total RNA transcription is lower at night than during the day, resulting in lower levels of H3K4me3 and H3K9ac modifications at night. Alternatively, activation of those genes at the night may not require both H3K4me3 and H3K9ac modifications. Other histone modifications could also be involved in evening-phased genes. Further examination of additional histone marks will provide a comprehensive view for diurnal rhythms of histone modifications and gene expression.
For H3K4me3, the rhythm for regulating circadian-mediated gene expression is associated with expression of the H3K4me3 methyltransferase gene SDG2, which is expressed at significantly higher levels during the day than at night; on the other hand, the expression level of the H3K4me3 demethylase gene (JMJ14) is low during the day but high at night. This provides a basis for transcriptional regulation of rhythmic histone modifications. Indeed, both global H3K4me3 and morning-phased H3K4me3 levels peak during the day (ZT12) (Additional file 1: Figure S4), suggesting that H3K4me3 is largely established during the day. However, the global H3K4me3 peak is delayed 8 h, compared to the expression peak of H3K4me3 methyltransferase SDG2 (ZT4). This could be related to the protein and mRNA that accumulate at different levels in a diurnal cycle. For example, CCA1 binds to the JMJ14 promoter at ZT0 in our ChIP-qPCR study but JMJ14 is not a target gene at ZT2 or ZT14 in the ChIP-seq analysis [39]. CCA1 protein peaks 3 h later than mRNA in maize [56]. The decrease of JMJ14 transcripts during the day period may result in the presence of H3K4me3 for a longer time. In addition, the dynamic histone modifications depend on maximal activity of histone-modifying enzymes that may change over time. Notably, other members of H3K4me3 writers (e.g., ATX1) and erasers (e.g., JMJ15) may also contribute to expression rhythms of circadian-mediated regulatory networks.
How histone methyltransferases and histone demethylases are recruited to their target genes remains unknown. One possibility is that histone methyltransferases such as ATX1 can be recruited by the transcriptional machinery through interaction with the Ser5-phosphorylated CTD of Pol II for transcription initiation and elongation [46]. Alternatively, histone modifications could be established through the interaction of sequence-specific transcription factors with histone methyltransferases/demethylases to regulate target gene expression [57, 58]. The targeting and binding time of transcription regulators for specific sequences could affect the time of H3K4me3 deposition. A recent study reported that lysine-specific demethylase 1 (LSD1)-like histone demethylases (LDL1 and LDL2) can interact with CCA1/LHY to repress TOC1 expression [59]. Moreover, LDL1 and LDL2 interact with HDA6 to co-regulate TOC1 expression by histone demethylation and deacetylation. It is probable that this type of histone complex can include JMJ14 and other factors to regulate expression of clock genes. It is equally plausible that histone methyltransferases and acetyltransferases could be in the same complex to carry out the opposite direction of circadian-mediated gene expression. It will be interesting to further test how specific transcription factors including circadian clock regulators recruit other histone methyltransferases and demethylases to establish rhythmic histone modifications that exert growth and developmental regulation in plants and animals.
Histone modifications regulate expression of key circadian clock regulators. A model showing feedback regulation between histone H3K4 trimethylation (H3K4me3) and circadian clock regulation through interactions with central clock oscillators (CCA1/LHY) and H3K4me3 writers such as SDG2 and erasers such as JMJ14. Solid lines demonstrated direct effects of positive (arrow) and negative (stop) interactions between CCA1 and JMJ14 (red) and SDG2 (green); dashed lines proposed indirect effects of positive (arrow) and negative (stop) interactions between CCA1 and JMJ14 (red) and SDG2 (green). CCA1/LHY can form homodimers for autoregulation (gray lines) in positive (arrow) and negative (stop) directions
Methods
Plant materials and growth conditions
A. thaliana ecotypes Col-0 and Ws were used in this study. cca1 lhy mutant in Ws background and CCA1-OX transgenic lines in Col-0 background were generated as previously described [27]. The mutants sdg2 (SALK-021008) and jmj14 (SALK_135712) were obtained from Arabidopsis Biological Resource Center (ABRC). The primer sequences for genotyping are listed in Additional file 5: Table S4. For diurnal conditions, plants were grown under the light/dark (L/D) cycle of 16 h/L and 8 h/D at 22 °C.
RNA extraction and qRT-PCR
Total RNA was isolated from aerial rosette leaves of ~ 3-week-old plants using Plant RNA Reagent (Thermo Fisher Scientific, Waltham, Massachusetts). After digestion by RNase-Free DNase (Promega, Madison, Wisconsin), total RNA (1 μg) was used to produce first-strand cDNA with the Omniscript RT Kit (Qiagen, Valencia, California). The cDNA was used as the template for qRT-PCR using FastStart Universal SYBR Green Master in a LightCycler®96 System (Roche, Indianapolis, Indiana). The relative expression level was quantified using the internal control ACT7 (AT5G09810). Three biological replicates were performed for each experiment, and three technical replicates were used for each biological replicate in qRT-PCR analysis. Student’s t test and calculation of error bars were performed to determine the significance level in each comparison using three biological replicates. The primer sequences are listed in Additional file 5: Table S4.
Histone extraction and Western blot analysis
Histone extraction was performed as described previously [60]. An aliquot (1 g) of aerial rosette leaves from ~ 3-week-old plants were ground to powder with liquid nitrogen and then suspended in NIB buffer (250 mM sucrose, 60 mM KCl, 15 mM NaCl, 5 mM MgCl2, 1 mM CaCl2, 15 mM PIPES, pH 6.8, 0.8% Triton X-100, 1 mM PMSF, and 1x cocktail). The homogenized solution was filtered through four layers of cheesecloth; the supernatant was centrifuged at 3000g for 20 min at 4 °C. Pellets were resuspended in 800 μl 0.4 N H2SO4 and incubated with a rotator overnight at 4 °C. After centrifugation at 16,000g at 4 °C for 10 min, the supernatant was transferred into a new tube by adding 264 μl trichloroacetic acid. Samples were incubated overnight at 4 °C and centrifuged again at 16,000g at 4 °C for 10 min to precipitate histone proteins. After washing by acetone, proteins were dissolved in water for Western blot analysis. The antibodies were anti-histone H3 (Abcam; Ab1791, 1:5000 dilution), anti-H3K4me3 (Abcam; Ab8580, 1:5000 dilution), and anti-H3K9Ac (Abcam; Ab10812, 1:5000 dilution).
Chromatin immunoprecipitation (ChIP) and ChIP-seq
Chromatin immunoprecipitation was performed following the published protocol [61]. Briefly, 2 g of aerial rosette tissues from ~ 3-week-old A. thaliana (Col-0) were crosslinked in 1% formaldehyde, and the chromatin fractions were isolated and purified. Antibodies against histone H3 (Abcam; ab1791), H3K4me3 (Abcam; Ab8580), H3K9ac (Abcam; Ab10812), or CCA1 (Abiocode, R1234-3) were added to the chromatin extracts (with a 1:1000 dilution) and incubated overnight at 4 °C with gentle rotation. Immunoprecipitated DNA was extracted using Protein A/G Magnetic Beads (Thermo Fisher Scientific, Waltham, Massachusetts). Immunoprecipitated DNA (IP DNA) and input DNA (without IP but with incubation of beads) were subjected to reverse cross-linking and purified using QIAquick PCR Purification Kit (Qiagen, Valencia, California). For ChIP-seq, libraries with three biological replicates were constructed from IP DNA and non-IP DNA using NEBNext® Ultra™ II DNA Library Prep Kit for Illumina (NEB, Ipswich, Massachusetts). For ChIP-qPCR, immunoprecipitated DNA was used as the template for qPCR using FastStart Universal SYBR Green Master in a LightCycler®96 System (Roche, Indianapolis, Indiana). The relative enrichment was quantified using the internal control AT2G26560. Three biological replicates were performed for each experiment with three technical replicates for each biological replicate. Student’s t test was used to determine the significance level in each comparison.
ChIP-seq data analysis
ChIP-seq dataset generated in this study and the published ChIP-seq datasets in sdg2 and jmj14 mutants were analyzed as follows. After removing adapters, low-quality reads, and PCR duplicates, the filtered pair-end reads were aligned to A. thaliana reference genome (TAIR10) using Bowtie2 (Version 2.3.0) with parameters “bowtie2 --no-mixed --no-discordant” [62]. Only uniquely and concordantly mapped reads were used for further analysis. In order to adjust different sequencing depths among different samples, the uniquely mapped reads were “down-sampled” to the lowest read number of samples, as previously described [56]. Enriched peaks were identified using Model-based Analysis of ChIP-Seq (MACS, version 2.1.0) with the option “macs2 callpeak -g dm -q 0.01” [63]. Peaks were used for further analysis only if they were present in two biological replicates. To make a master-peak list for all time points, the peaks obtained from each time point were merged. When different samples were compared, the read coverage was normalized to 1× sequencing depth to eliminate bias of different insertion lengths in pair-end reads. The genes overlapping with peaks in the gene body or ± 1-kb regions were defined as peak-overlapped genes. The peaks were compared between the wild type and cca1 lhy mutant using one-way ANOVA, and their significance was determined with P < 0.05. As no biological replicate was used from the ChIP-seq datasets in sdg2 and jmj14 mutants, we used 1.5-fold changes to determine differential peaks between the wild type and mutants. The non-parametric algorithm JTK_CYCLE was used to detect rhythmic peaks and identify the phase of each peak with Padj < 0.01 [64]. Gene Ontology (GO) analysis was performed using DAVID [65].
Additional files
Additional file 1:Supplementary Figure S1–S8. (PS 3131 kb)Additional file 2:Table S1. List of peaks with different H3K4me3 enrichment between WT and cca1 lhy mutant. (XLSX 101 kb)Additional file 3:Table S2. List of all identified H3K4me3 and H3K9ac peaks with coverage levels at all time points. (XLSX 11750 kb)Additional file 4:Table S3. The list of diurnal H3K4me3 and H3K9ac peaks with coverage levels at all time points. (XLSX 2978 kb)Additional file 5:Table S4. Primers used in this study. (DOCX 14 kb)