Skip to main content

Histone H3 lysine 36 methylation affects temperature-induced alternative splicing and flowering in plants

Abstract

Background

Global warming severely affects flowering time and reproductive success of plants. Alternative splicing of pre-messenger RNA (mRNA) is an important mechanism underlying ambient temperature-controlled responses in plants, yet its regulation is poorly understood. An increase in temperature promotes changes in plant morphology as well as the transition from the vegetative to the reproductive phase in Arabidopsis thaliana via changes in splicing of key regulatory genes. Here we investigate whether a particular histone modification affects ambient temperature-induced alternative splicing and flowering time.

Results

We use a genome-wide approach and perform RNA-sequencing (RNA-seq) analyses and histone H3 lysine 36 tri-methylation (H3K36me3) chromatin immunoprecipitation sequencing (ChIP-seq) in plants exposed to different ambient temperatures. Analysis and comparison of these datasets reveal that temperature-induced differentially spliced genes are enriched in H3K36me3. Moreover, we find that reduction of H3K36me3 deposition causes alteration in temperature-induced alternative splicing. We also show that plants with mutations in H3K36me3 writers, eraser, or readers have altered high ambient temperature-induced flowering.

Conclusions

Our results show a key role for the histone mark H3K36me3 in splicing regulation and plant plasticity to fluctuating ambient temperature. Our findings open new perspectives for the breeding of crops that can better cope with environmental changes due to climate change.

Background

Plants sense changes in the environment and adapt their development accordingly. Temperature is one of the environmental signals that strongly affects plant development. For example, plants are able to adapt their organ shape in relation to the temperature they experience, a phenomenon called “thermo-morphogenesis” [1]. As such, elevated temperature promotes hypocotyl elongation and plants also show differences in leaf morphology depending on the growing temperature [2,3,4,5]. The circadian clock, the internal clock that confers an organism the innate ability to measure time, is also influenced by temperature [6,7,8,9]. Another important plant trait controlled by temperature is flowering time, the switch from the vegetative to the reproductive phase. An increase in ambient temperature promotes flowering in Arabidopsis and, consequently, at high ambient temperatures these plants flower earlier than at colder conditions [10,11,12].

A change in ambient temperature affects gene functioning at transcription level and at post-transcriptional level, e.g. through protein conformation and stability. Recently, it has been shown that ambient temperature affects phytochrome B (PHYB) inactivation during the night [13]. Moreover, protein stability can be influenced by temperature, for example the flowering time regulator short vegetative phase (SVP) is subject to faster protein degradation at higher temperature [14]. Finally, ambient temperature also affects alternative splicing [15]. Alternative splicing (AS) is the process that produces different mature messenger RNAs (mRNAs) from a single pre-mRNA [16]. AS is a rapid and adjustable process that allows a fast response to external stimuli and facilitates plasticity in plant development [17]. In plants, key regulators of various processes such as circadian clock, flowering time, and reproduction are regulated via temperature-induced AS [14, 15, 18,19,20,21,22,23,24]. Ambient temperature changes also affect splicing of AS regulators and thereby plant development [25, 26]. As such, AS was suggested to act as “molecular thermosensor” in plants [15]; however, how temperature can affect splicing of a specific subset of intron-containing genes is poorly understood. In higher eukaryotes, it has been shown that pre-mRNA splicing occurs co-transcriptionally in the nucleus. Thus, while RNA polymerase is synthesizing the pre-mRNA, the splicing machinery is recruited to perform its function [27]. The co-transcriptional nature of splicing allows a contribution of the chromatin landscape in splicing regulation [28,29,30]. In mammals, two models have been proposed to explain the interplay between chromatin landscape and AS: the kinetic model and the chromatin-adaptor model [31]. Evidence for both modes of action has been reported. According to the kinetic model, the chromatin landscape can influence the transcription elongation rate and thereby affects splicing. For example, histone acetylation allows opening of the chromatin leading to a fast elongation rate that promotes splice site skipping. To the contrary, silencing histone marks lead to a more compact chromatin and slow elongation rate allowing recognition of weak splicing sites [29]. In the chromatin-adaptor model, histone modifications are bound by readers able to recruit splicing factors. An example is the H3K36me3-MRG15-PTB chromatin-adaptor system. The H3K36me3 histone mark serves as an anchor for binding of the adaptor protein MRG15, which in turn interacts with the splicing regulator polypyrimidine tract-binding (PTB) recruiting it to the nascent pre-mRNA [32].

In plants, it is not known if a similar link between chromatin and splicing exists and whether such a mechanism would be ambient temperature sensitive. H3K36me3 was found to be prevalent in transcribed regions and therefore a function in transcriptional elongation has been postulated for this histone mark [33]. Here we investigated the potential contribution of this mark in the regulation of ambient temperature-induced AS and flowering time control in the model plant Arabidopsis thaliana.

Results

Identification of temperature-induced differentially spliced events

To study the proposed linkage between the H3K36me3 mark and ambient temperature-induced AS, we first identified transcriptome-wide differentially spliced (DiS) events upon a temperature change and investigated the temporal dynamics of this process. Arabidopsis plants were grown at 16 °C in short day conditions (8 h light-16 h dark) for five weeks and then transferred to 25 °C. With RNA sequencing (RNA-seq), we analyzed the changes in the transcriptome at one, three, and five days after the temperature change. This experimental setup allows to monitor initial changes in splicing pattern upon a temperature change and to investigate the stability of these effects. For each time point, we quantified AS setting one isoform as reference and then calculating the “Percentage of Spliced In” (PSI) for each type of event (intron retention [IR], mutually exclusive exon [MXE], exon skipping [ES], alternative 3′ site [A3], alternative 5′ site [A5]; Additional file 1: Table S1). We defined isoforms that present a significant difference in PSI (ΔPSI) between two conditions as temperature-induced DiS. In total, we identified 683 DiS events, derived from 511 distinct genes, involved in a plethora of biological processes (Additional file 1: Table S1, Fig. 1a, and Additional file 2: Figure S1A).

Fig. 1
figure 1

Dynamics of temperature-induced DiS events. a Heat map showing the percentage of spliced in (PSI) values obtained at each condition for significant DiS events. b Partial gene model depicting the detected DiS events of FLOWERING LOCUS M (FLM), MADS AFFECTING FLOWERING 2 (MAF2), and ARGININE/SERINE-RICH SPLICING FACTOR 40 (RS40), and corresponding plot area showing RNA-seq results. Coverage plot reads within the AS regions are shown in black, other reads are in gray. We detected a significant difference in IR events for FLM [35] and MAF2 [34] and MXE events for FLM [21] and RS40 upon the temperature change. c Venn diagram shows the overlap between DiS and totally differentially expressed genes (DEGs) one day upon the temperature change. DEGs are defined as genes with a change in expression of log2 Fold Change |1| and adjusted p value according to the BH method for controlling false discovery < 0.05

We found DiS events in genes previously reported to be subject to temperature-induced AS, such as in the flowering time regulators FLOWERING LOCUS M (FLM) [14, 21] and MADS AFFECTING FLOWERING 2 (MAF2) [22, 34] (Additional file 2: Figure S1B) and the clock components PSEUDO-RESPONSE REGULATOR 3 (PRR3) and PRR7 [8, 18]. For FLM, we observed an increase in the retention of the second intron and a MXE event [14, 21, 35]. The retention of the second intron leads to a transcript targeted for degradation by nonsense-mediated mRNA decay (NMD) and as such decreased amount of functional protein at higher temperature [35], while the MXE event is supposed to result in the production of two proteins with opposite function in flowering time regulation [21]. The IR at MAF2 locus generates an isoform with a premature stop codon that is most likely targeted for NMD and the increased abundance of this isoform at higher temperature leads to a decrease of the functional flowering time repressor [34]. PRR3 and PRR7 were shown previously to undergo AS upon a change to low temperature from 20 °C to 4 °C [18]. In our experiment, we see a change in retention of the third intron of PRR7 as previously reported for a change to low temperature, as well as new events for PRR7 and PRR3 that appear to be specific to the increase in ambient temperature in our experimental setup (Additional file 1: Table S1).

Moreover, we found splicing regulators DiS upon the temperature change, such as RS40 (Fig. 1b), SC35-LIKE SPLICING FACTOR 33 (SCL33), SERINE/ARGININE-RICH PROTEIN SPLICING FACTOR 34 (SR34), and the U2 snRNP auxilliary factors U2AF65A and U2AF65B.

Additionally, we also identify novel temperature-induced splicing events in master regulators of various biological processes, such as the flowering time regulator SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 2 (SPL2) [36] and MODIFIER OF SNC 1 (MOS1) [37], the regulators of auxin signaling AUXIN RESPONSE FACTOR 2 (ARF2) [38] and PIN-FORMED 7 (PIN7) [39], and the chromatin remodelers VIN3-LIKE 2 (VIL2) [40] and SPLAYED (SYD) [41,42,43] (Fig. 1b and Additional file 1: Table S1).

The majority of the genes show DiS within one day upon the temperature change and after that their PSI remained stable over time (Fig. 1a, b), revealing a response to a change in ambient temperature within 24 h. Surprisingly, we found only a small overlap (13 genes) between differentially expressed genes (DEGs) (Additional file 3: Table S2) and DiS genes (Fig. 1c), indicating a different basis for the effect of ambient temperature on these two molecular processes.

In conclusion, a change in ambient temperature triggers fast and stable changes in splicing of many key regulatory genes, including various flowering time regulators.

DiS genes are enriched in H3K36me3 modification

In mammalian studies, the H3K36me3 mark was found to be enriched in genes undergoing AS and therefore was suggested to play a predominant role in the AS regulation. We hypothesize that the H3K36me3 deposition contributes to the regulation of ambient temperature-induced AS in plants. To investigate the putative role of this mark, we first identified H3K36me3 modified genomic regions in plants growing at 16 °C and in plants that were transferred from 16 °C to 25 °C. Since we observed that most of the DiS events already occur one day after the temperature shift and then remain present (Fig. 1), we selected this time point to investigate the H3K36me3 profile. H3K36me3 chromatin immunoprecipitation sequencing (ChIP-seq) experiments were performed on four biological replicates for each temperature and H3K36me3-marked regions were identified. We found a total of 59,736 regions to be enriched in H3K36me3 (Additional file 4: Table S3). Analyzing the distribution of H3K36me3 along the whole gene body, we observed that H3K36me3 containing histones are prevalent at the beginning of the gene body, with a peak approximately 0.5 kb after the transcription start site (TSS) and predominantly in expressed genes (Fig. 2a). Remarkably, we also observed that DiS genes are enriched in the H3K36me3 mark, while this is not the case for DEGs (Fig. 2a). Around 96% of the DiS genes have an H3K36me3 enriched region in the gene body, while this is only 65% for the DEGs (Additional file 2: Figure S2A). Among the H3K36me3-marked DiS genes, we found key regulators of ambient temperature-sensitive processes, such as the flowering time regulators FLM, MAF2, and FCA, the circadian clock regulators PRR3 and PRR7, and several splicing regulators, such as RS40 and SC35-LIKE SPLICING FACTOR 33 (SCL33) (Fig. 2b).

Fig. 2
figure 2

Temperature-induced DiS genes are enriched in H3K36me3. a H3K36me3 profile around the TSS and the transcription termination site (TTS), for different classes of genes. b H3K36me3 profile for selected genes with temperature-induced DiS events

The quantitative comparison between the profile obtained before and after the temperature switch (Additional file 4: Table S3) revealed a higher level of H3K36me3 located around 0.5 kb after the TSS in plants exposed at 25 °C (Additional file 2: Figure S2B). For example, we observed higher levels of H3K36me3 at the TSS of the DiS genes PIN7 and VIL2 at 25 °C compared with 16 °C (Additional file 4: Table S3).

H3K36me3 regions tend to be more broad at 25 °C compared with 16 °C (Additional file 2: Figure S2B and S3C). For example, we observed an increased in the peak width at FLM, SPL2, SCL33, and at RUBISCO ACTIVASE (RCA) loci (Additional file 4: Table S3).

A few genes showed a quantitative change in H3K36me3 at the location of a DiS event under these conditions, for example, for the circadian clock genes PRR3 and PRR7 (Additional file 2: Figure S3); however, for the majority of cases there was no perfect co-localization.

Thus, overrepresentation of the H3K36me3 histone mark was observed in DiS genes and changes in the mark deposition were observed after increasing temperature, suggesting a role for the mark in the regulation of ambient temperature-induced AS.

H3K36me3 contributes to regulation of temperature-induced DiS events

To elucidate the contribution of H3K36me3 in temperature-induced splicing regulation, we investigated the effect of decreased levels of H3K36me3 on temperature-induced AS. In Arabidopsis, the deposition of H3K36me3 is regulated by the histone methyltransferases SDG8 and SDG26 and, accordingly, levels of H3K36me3 are low in sdg8-2 and sdg26-1 mutant plants [44, 45]. We performed an RNA-seq experiment on wild-type (WT) and mutant plants, continuously grown at 16 °C, or one day after a temperature shift to 25 °C. For each genotype, we retrieved DiS events upon the temperature change (Additional file 5: Table S4). Significant differences in temperature-induced DiS events were observed between WT and the H3K36me3 mutants upon the temperature change (Fig. 3). In the sdg8-2 and sdg26-1 mutant different types of splicing events are affected by the temperature change (Additional file 2: Figure S4A), indicating that the histone mark is not involved in the regulation of a specific type of event.

Fig. 3
figure 3

H3K36me3 contributes to regulation of temperature-induced splicing. a The graph shows distribution of differences in PSI between 16 °C and 25 °C (ΔPSI) of significant DiS events identified in WT and corresponding ΔPSI distributions in sdg8-2 and sdg26-1 (Additional file 2: Figure S2C). b Graph shows plot of empiric cumulative distribution function (ecdf) for absolute ΔPSI of WT, sdg8-2, and sdg26-1 for significant DiS events in the comparison of WT 16 °C and 25 °C. Distribution of ΔPSI are significantly different between WT and sdg8-2 or sdg26-1 according to the Kolmogorov–Smirnov test c Partial gene model of two examples depicting the DiS event and corresponding plot area showing RNA-seq results. Reads within the AS regions are shown in black, other reads are in gray. At3g08840 is only significantly DiS in the WT while At3g55080 is significantly DiS in WT as well as in sdg26-1. d Gene ontology enrichment for genes of which the splicing is affected by ambient temperature conditions in WT, sdg8-2, and sdg26-1. The histogram includes all overrepresented categories with at least five genes and p value < 0.05 for at least one genotype

Consistent with the contribution of H3K36me3 in temperature-induce splicing, only a minority of DiS genes in WT were also DiS in sdg8-2 and sdg26-1 (Additional file 2: Figure S4B), while many events are affected by the decreased levels of H3K36me3 (Additional file 5: Table S4). For example, reduction of H3K36me3 in the mutants disturb temperature-induced AS of flowering time genes, such as FLM and MAF3, the splicing regulators SCL33 and U2AF65B, and the circadian clock gene PRR7 (Additional file 5: Table S4).

Because we noticed various cases in which DiS and H3K36me3 deposition were not perfectly co-localized, we investigated if there is a correlation between the location of temperature-dependent H3K36me3 deposition and temperature-dependent AS. We categorized temperature-induced AS events as H3K36me3-dependent if the events are affected by loss of H3K36me3 (significantly DiS in the WT but not in the sdg8 and sdg26 mutants) and H3K36me3-independent if the events are not affected by the loss of H3K36me3 (significantly DiS in all the genotypes). H3K36me3-dependent events tend to be located closer to an H3K36me3 region compared with H3K36me3-independent events (Additional file 2: Figure S4C).

We further characterized to which biological processes the genes are associated that undergo DiS in the sdg8-2 and sdg26-1 mutants (Fig. 3d). We observed GO categories such as “regulation of flower development” or “gene silencing” to be no longer overrepresented in sdg8-2 and sdg26-1, while GO categories related to metabolic processes appeared higher overrepresented in the mutants.

Finally, we investigated if sdg8-2 and sdg26-1 mutants are only affected in temperature-induced AS or generally in the regulation of splicing. We calculated ΔPSI between WT and sdg8-2 or sdg26-1 growing at 16 °C. Only a small set of genes show changes in splicing pattern in the mutants (Additional file 6: Table S5). These results highlight a specific function for H3K36me3 in the regulation of high temperature-induced AS.

In conclusion, the higher levels of H3K36me3 in high ambient temperature-induced DiS genes (Fig. 2) and the altered temperature-induced AS for several genes in H3K36me3 mutants (Fig. 3) demonstrate a role of H3K36me3 in regulating temperature-induced AS in Arabidopsis.

Mutations in H3K36me3 writers, eraser, or readers affected high temperature-induced flowering

In order to investigate the biological relevance of the observed molecular events, we investigated whether a lack of H3K36me3 deposition or lack of proteins recognizing this histone modification affects ambient temperature-controlled flowering. Arabidopsis WT and the H3K36me3 methyltransferase mutants, sdg8-2 and sdg26-1, were grown at 16 °C for five weeks and then maintained at 16 °C or transferred to 25 °C. We quantified temperature response as the ratio in days to bolting (DTB, Fig. 4) or rosette leaf number (RLN, Additional file 2: Figure S5) at low ambient temperatures over the DTB or RLN at high ambient temperatures. Thus, for a genotype that accelerates flowering in higher ambient temperatures, this ratio is higher than 1, as we observed in WT, while for a plant not responding to temperature, this ratio would be 1 [46]. As expected, the increase in temperature accelerated flowering in WT plants (Fig. 4 and Additional file 2: Figure S5). In contrast, sdg8-2 and sdg26-1 plants flowered at similar time points (ratio close to 1) at both temperatures (Fig. 4b and Additional file 2: Figure S5A). The significant difference in temperature response between WT plants and the H3K36me3 methyltransferase mutants (Fig. 4b and Additional file 2: Figure S5A) indicates that H3K36me3 is required for temperature-induced flowering.

Fig. 4
figure 4

Temperature induced-flowering is compromised in mutant of H3K36me3 writers, eraser, and readers. a Cartoon illustrates regulation of H3K36me3 and splicing according to the adaptor-complex model in Arabidopsis. H3K36me3 modified histones are bound by an adapter, such as MRG1 and MRG2, that interacts with splicing regulators (SRs), such as PTB. The splicing regulators act in promoting or preventing AS. b Flowering time of WT, sdg8-2, and sdg26-1 plants growing at 16 °C constantly or moved to 25 °C after five weeks. In the top histogram, flowering time is expressed as numbers of DTB. In the bottom histogram, the ratio of DTB between the two conditions is shown. In plants with low levels of H3K36me3, flowering is not promoted by higher temperature. ** indicates significant differences at p value < 0.001 while * indicates p value < 0.05 according to the Student’s t-test. c Flowering time of WT, jmj30-1, and JMJ30-OE plants is shown as in (b). Loss of JMJ30 increases the response to temperature-induced flowering, while an increase in JMJ30 activity causes an opposite effect. d Flowering time of WT, mrg1-1 mrg2-3 plants are shown as in (b). mrg1-1 mrg2-3 plants are less responsive to temperature

We also characterized the flowering time response upon a temperature change for the H3K36me3 demethylase mutant jmj30-1 and for a previously published JMJ30 overexpression line (JMJ30-OE) [47] (Fig. 4c and Additional file 2: Figure S5B). In agreement with the previous experiment we found that jmj30-1 is more sensitive to the temperature while JMJ30-OE is less responsive in temperature-induced flowering (Fig. 4c and Additional file 2: Figure S5B).

In mammals, H3K36me3 has been shown to regulate AS via the adapter-complex model: the protein complex H3K36me3–MRG15–PTB forms a splicing–chromatin co-transcriptional architecture that can bring the splicing machinery in the vicinity of chromatin (Fig. 4a) [32]. In Arabidopsis, the MRG15 homologs MORF RELATED GENE 1 (MRG1) and MRG2 are able to bind to H3K36m3 modified histones, as is the case for the human MRG15 [48]. Moreover, the characterization of PTB homologs in Arabidopsis (PTB1, PTB2, and PTB3), reveals their role in the regulation of splicing of flowering time genes, including FLM [49]. In the proposed adaptor-complex model [31], the H3K36me3 readers MRG1 and MRG2 proteins play a pivotal role. In agreement with this model, we found that mrg1-1 mrg2-3 mutant plants are less sensitive to temperature-induced flowering than the WT (Fig. 4d and Additional file 2: Figure S5C).

Remarkably, none of the investigated mutant genotypes show strong pleiotropic defects throughout plant development (Additional file 2: Figure S6), in contrast to what was previously observed for core-splicing factor mutants [25, 50,51,52,53,54]. This is in line with our finding that in the H3K36me3 mutants splicing as such is not much affected, but that the defects are specifically for high ambient temperature-induced DiS.

In conclusion, mutations in H3K36me3 writers, an eraser, and readers affect plant response to ambient temperature changes, suggesting a role of this histone mark in plasticity to this environmental cue.

Discussion

In plants, AS of pre-mRNA is a key mechanism to generate protein diversity and also allows rapid changes in transcript abundance via NMD [55]. In temperature sensing, AS represent a fast and flexible way to confer plasticity in plant development. Our time-series experiments highlight a rapid and rather stable response of changes in splicing upon a temperature change. We found genes belonging to many biological processed to be regulated by temperature-induced AS, such as circadian rhythm and flowering time control, as well as many genes coding for splicing regulators.

The co-transcriptional nature of splicing creates opportunities for cross-regulation between splicing and chromatin landscape. Recent finding indicated that introns are spliced immediately after their transcription [56], indicating the need of a fast recruitment of splicing regulators to the nascent pre-mRNA. In mammals, histone modifications contribute to the regulation of splicing, but whether epigenetic regulation of splicing occurs in plants is unknown. Our results indicate a direct role for histone modifications in AS regulation in plants, highlighting conservation between mammals and plants. However, while in mammals H3K36me3 prevalently regulates ES [32], we found that H3K36me3 is equally involved in the regulation of all type of AS events in plants. Furthermore, in Arabidopsis H3K36me3 specifically affects ambient temperature-induced splicing regulation.

While our findings highlight a key role for H3K36me3 in the regulation of ambient temperature-induced splicing, the exact molecular mode of action still remains to be elucidated. Our results do not exclude the kinetic model, neither the adapter-complex model, which have been proposed to explain the linkage between chromatin modifications and AS [28, 31]. Instead, our data suggest that both mechanisms can play a role in H3K36me3 mediated temperature-induced AS. A reduced response to temperature in flowering induction observed in the mrg1-1 mrg2-3 double mutant suggests that H3K36me3 presence affects splicing outcome at least partially by influencing the recruitment of splicing regulators via a chromatin-binding protein (adapter-complex model). Nevertheless, the preferential location of the H3K36me3 mark in the proximity of H3K36me3 dependent AS events could also point to an effect on transcription speed in support of the kinetic model. In agreement, transcription rate was found to be affected by temperature in Arabidopsis [57]. Additionally, deposition of the Histone 2 variant H2A.Z, which was associated to temperature-mediated flowering time control in Arabidopsis [58], was recently found to affect transcription elongation kinetics and AS in yeast, and especially under destabilizing conditions, such as high temperature [59, 60]. These observations and our results are in favor of the kinetic model. Nevertheless, further research will be essential to characterize the molecular mode(s) of action of H3K36me3 temperature-induced AS regulation in depth.

We show that H3K36me3 plays a key role in the regulation of temperature-induced flowering. Flowering time is a key developmental process, because a proper timing in the switch from the vegetative to the reproductive phase ensures plant fitness. We found that key regulators in the ambient temperature-controlled flowering time pathway are subject to H3K36me3 level dependent AS. However, H3K36me3 may also affect temperature-induced flowering indirectly by affecting splicing of upstream regulators of flowering time genes. For example, the temperature DiS gene VIL2 regulates the ambient temperature-associated flowering time gene MAF5 [40]. Alternatively, the effects on ambient-temperature dependent flowering time can be the result of differential gene expression in the various studied H3K36me3-linked mutants.

Conclusions

AS of pre-mRNA is a key mechanism to generate protein diversity or to regulate transcript isoform abundance and it allows rapid changes in the transcriptome. In temperature sensing, AS represents a fast and flexible way to confer plasticity in plant development. The co-transcriptional nature of splicing creates opportunities for cross-regulation between splicing and chromatin landscape. In mammals, histone modifications contribute to the regulation of splicing. Here we show a direct role for histone modifications in AS regulation in plants. We reveal a key role for H3K36me3 in the control of temperature-induced AS and in temperature-dependent flowering time control in Arabidopsis.

Besides that AS of key regulatory genes confers plasticity, because it is a fast and reversible response to changing environmental conditions, our findings on the role of H3K36me3 in this process open interesting prospective for a memory of temperature mediated AS. It is tempting to speculate that temperature variation is memorized in the chromatin landscape via H3K36me3 deposition, leading to a specific splicing pattern, which allows the memory of environmental conditions during the life cycle and a possible adaptation to a rapidly changing climate.

Methods

Plant material and growth conditions

Seeds of knock-out lines, sdg8-2 Salk_026442 [45], sdg26-1 Salk_013895 [44], jmj30-1 (N67787) [47], and the overexpression line JMJ30 (N67797) [47], were obtained from the Nottingham Arabidopsis Stock Center (NASC). Seeds of knock-out line mrg1-1 mrg2-3 [48] were a gift from Toshito Ito (Temasek Life Sciences Laboratory, Singapore). All plants were grown under short day conditions (8 h light, 16 h dark) on rock-wool at 16 °C for five weeks. Plants were left at the same temperature or moved at 25 °C (6 h after lights on). WT Col-0 plants for the time series RNA-seq experiments were grown in growth cabinets with LED lamps with light at intensity of 200 μmol m-2 s-1 and humidity of 75%. WT and mutant lines for the RNA-seq and flowering time assay were grown in walk-in climate growth chambers with fluorescent lamps with light at intensity of 185 μmol m-2 s-1 and humidity of 65%.

RNA-seq experiments

Three biological samples were generated for each condition. For each sample, tissue from approximately ten plants was collected 4 h after light on for the time-series experiment and between 4 h and 7 h after light on for the mutant experiment. Using jeweler’s forceps, leaves were removed to obtain SAM enriched tissue. Total RNA was extracted using the InviTrap Spin Plant RNA Mini Kit (REF: 1064100300) from Invitek according to the manufacturer’s protocol. DNAse treatment was performed to remove genomic DNA. DNase I digestion was performed on total RNA using Turbo DNase from Ambion according to the manufacturer’s protocol. RNA integrity was checked on 1% (w/v) agarose gel after DNase I treatment. Samples were prepared for Illumina sequencing using the TruSeq Stranded mRNA Sample Prep kit (REF: 15032613) from Illumina according to the manufacturer’s protocol. Libraries were analyzed on the Bioanalizer and quantified with the qBit before pooling for sequencing on HiSeq2500. Two lanes on a 125pb PE run were used.

RNA-seq data analysis

RNA-seq reads were mapped against the Arabidopsis thaliana genome version TAIR10 (www.arabidopsis.org) using TopHat2 [61]. To improve TAIR10 transcript annotation, reference-based full-length transcript isoform reconstruction was performed for each experiment separately using Cufflinks [62]. Cuffmerge, which is part of the Cufflinks package, was finally used for merging the individual cufflinks results into an overall set of full-length transcripts. Custom python scripts were used for detecting the following type of alternative splicing events: IR, ES, A5, A3, and MXE. MISO software [63] was used to quantify AS-events (PSI values) in each individual sample and in pooled samples that were generated by merging the replicates of each condition. MISO, which also used for the differential splicing analyses, does not have built-in methods for analyzing experimental replicates. As suggested by the authors of MISO, we used the pooled samples for the actual differential splicing test and the individual replicates for filtering AS events based on the following two rules: first, we only considered those AS events that were supported by at least 20 isoform-specific reads in all at least two of the replicates of the conditions under comparison; second, the within-condition PSI differences were required to be smaller than the between-condition PSI differences. Finally, AS events that met the criteria and for which the compare-sample module of MISO returned a Bayes factor of at least 5 were considered significant.

To retrieve differentially expressed genes upon the temperature changes, the number of fragments mapping to TAIR10 annotated genes was determined using HTseq count [64]. Differentially expressed genes were detected using DESeq2 [65].

GO analysis

GO analysis was performed using the Cytoscape plugin BiNGO [66] using the most recent go-basic ontology file and Arabidopsis gene association file from The Gene Ontology Consortium [67]. Overrepresentation was tested using a hypergeometric statistical test and the Benjamini–Hochberg false discovery rate correction was used. GO category was considered significantly overrepresented if corrected p value < 0.05 and at least five genes were present.

ChIP-seq experiments

Four independent biological replicates for each temperature were generated. For each sample, 0.5 g of plant material was used per biological replicate. Material was collected from plants growing at 16 °C, as well as from plants one day after the plants were moved to 25 °C (5–6 h after lights on). Using jeweler’s forceps, leaves were removed to obtain SAM-enriched tissue. ChIP experiments were performed following a previously published protocol [68]. After the preclearing step, the sample was split, one aliquot was incubated with anti-H3K36me3 antibody (ab9050, abcam) [69, 70], and the other with anti-H3 antibody (ab1791, abcam). Samples were prepared for Illumina sequencing using the ThruPLEX DNA-seq Kit (Cat. No: R400429) from Rubicon Genomics according to the manufacturer’s protocol. Libraries were analyzed on the Bioanlizer and quantified with the qBit before pooling for sequencing on HiSeq2500. In total three lanes on the chip were used and sequenced 50 bp single run (SR).

ChIP-seq data analysis

FASTQ files were mapped to the Arabidopsis thaliana genome TAIR10, using Bowtie [61], using the default parameters. Reproducibility between biological replicates was assessed using the Pearson correlation for the genome-wide reads distribution at each pair of replicates using DeepTool [71]. We used “multiBamSummary” function with default parameters except for “bin size” that was set to 1 kb and “plotCorrelation” function of deepTools2 in Galaxy (http://deeptools.ie-freiburg.mpg.de/). Regions with very large counts were removed to avoid bias in the Pearson correlation (Additional file 2: Figure S7). To identify H3K36me3 modified nucleosome positions, we used DANPOS2 [72]. We used the “Dpeak” function in DANPOS2 with default parameters, except for the parameter – l (read extension length) that was set to 150 bp, the size of mononucleosomal DNA. Genomic regions were associated with genes if located within the start and the end of the gene using the “intersect” function of BEDTools [73]. The “Profile” function from DANPOS2 was used to retrieve H3K36me3 occupancy profiles around the TSS and the transcription termination site (TTS).

Flowering time assay

For flowering time analyses, plants were randomly arranged in trays and grown under white light (185 μmol m-2 s-1) in short day conditions (8 h light, 16 h dark) in walk-in climate chambers with humidity of 65%. Trays were rearranged every two days and nutrient was supplied by sub-irrigation. After five weeks of growth at 16 °C, half of the plants were moved to the 25 °C growth chamber. Flowering time was quantified by determining the time until the macroscopic appearance of the first flower bud (DTB), screening was done every two days and by counting RLN after bolting of all plants in the tray. The experiment was performed in three biological replicates with 15 plants in each replicate for each genotype/condition. Data are shown as average of ratios between the replicates. Student’s t-tests were calculated with GraphPad QuickCalcs.

References

  1. Quint M, Delker C, Franklin KA, Wigge PA, Halliday KJ, van Zanten M. Molecular and genetic control of plant thermomorphogenesis. Nat Plants. 2016;2:15190.

    Article  CAS  PubMed  Google Scholar 

  2. Crawford AJ, McLachlan DH, Hetherington AM, Franklin KA. High temperature exposure increases plant cooling capacity. Curr Biol. 2012;22:R396–7.

    Article  CAS  PubMed  Google Scholar 

  3. Gray WM, Ostin A, Sandberg G, Romano CP, Estelle M. High temperature promotes auxin-mediated hypocotyl elongation in Arabidopsis. Proc Natl Acad Sci U S A. 1998;95:7197–202.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Koini MA, Alvey L, Allen T, Tilley CA, Harberd NP, Whitelam GC, et al. High temperature-mediated adaptations in plant architecture require the bHLH transcription factor PIF4. Curr Biol. 2009;19:408–13.

    Article  CAS  PubMed  Google Scholar 

  5. Sun JQ, Qi LL, Li YN, Chu JF, Li CY. PIF4-mediated activation of YUCCA8 expression integrates temperature into the auxin pathway in regulating Arabidopsis hypocotyl growth. PLoS Genet. 2012;8:e1002594.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Box MS, Huang BE, Domijan M, Jaeger KE, Khattak AK, Yoo SJ, et al. ELF3 Controls thermoresponsive growth in Arabidopsis. Curr Biol. 2015;25:194–9.

    Article  CAS  PubMed  Google Scholar 

  7. Gould PD, Locke JCW, Larue C, Southern MM, Davis SJ, Hanano S, et al. The molecular basis of temperature compensation in the Arabidopsis circadian clock. Plant Cell. 2006;18:1177–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Mizuno T, Nomoto Y, Oka H, Kitayama M, Takeuchi A, Tsubouchi M, et al. Ambient temperature signal feeds into the circadian clock transcriptional circuitry through the EC night-time repressor in Arabidopsis thaliana. Plant Cell Physiol. 2014;55:958–76.

    Article  CAS  PubMed  Google Scholar 

  9. Salome PA, Weigel D, McClung CR. The role of the Arabidopsis morning loop components CCA1, LHY, PRR7, and PRR9 in temperature compensation. Plant Cell. 2010;22:3650–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Balasubramanian S, Sureshkumar S, Lempe J, Weigel D. Potent induction of Arabidopsis thaliana flowering by elevated growth temperature. PLoS Genet. 2006;2:980–9.

    Article  CAS  Google Scholar 

  11. Capovilla G, Schmid M, Pose D. Control of flowering by ambient temperature. J Exp Bot. 2015;66:59–69.

    Article  CAS  PubMed  Google Scholar 

  12. Verhage L, Angenent GC, Immink RG. Research on floral timing by ambient temperature comes into blossom. Trends Plant Sci. 2014;19:583–91.

    Article  CAS  PubMed  Google Scholar 

  13. Jung JH, Domijan M, Klose C, Biswas S, Ezer D, Gao MJ, et al. Phytochromes function as thermosensors in Arabidopsis. Science. 2016;354:886–9.

    Article  CAS  PubMed  Google Scholar 

  14. Lee JH, Ryu HS, Chung KS, Pose D, Kim S, Schmid M, et al. Regulation of temperature-responsive flowering by MADS-Box transcription factor repressors. Science. 2013;342:628–32.

    Article  CAS  PubMed  Google Scholar 

  15. Capovilla G, Pajoro A, Immink RG, Schmid M. Role of alternative pre-mRNA splicing in temperature signaling. Curr Opin Plant Biol. 2015;27:97–103.

    Article  CAS  PubMed  Google Scholar 

  16. Modrek B, Lee C. A genomic view of alternative splicing. Nat Genet. 2002;30:13–9.

    Article  CAS  PubMed  Google Scholar 

  17. Syed NH, Kalyna M, Marquez Y, Barta A, Brown JWS. Alternative splicing in plants - coming of age. Trends Plant Sci. 2012;17:616–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. James AB, Syed NH, Bordage S, Marshall J, Nimmo GA, Jenkins GI, et al. Alternative splicing mediates responses of the Arabidopsis circadian clock to temperature changes. Plant Cell. 2012;24:961–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Macknight R, Duroux M, Laurie R, Dijkwel P, Simpson G, Dean C. Functional significance of the alternative transcript processing of the Arabidopsis floral promoter FCA. Plant Cell. 2002;14:877–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Marquardt S, Raitskin O, Wu Z, Liu FQ, Sun QW, Dean C. Functional consequences of splicing of the antisense transcript COOLAIR on FLC transcription. Mol Cell. 2014;54:156–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Pose D, Verhage L, Ott F, Yant L, Mathieu J, Angenent GC, Immink RGH, Schmid M. Temperature-dependent regulation of flowering by antagonistic FLM variants. Nature. 2013;503:414.

    Article  CAS  PubMed  Google Scholar 

  22. Rosloski SM, Singh A, Jali SS, Balasubramanian S, Weigel D, Grbic V. Functional analysis of splice variant expression of MADS AFFECTING FLOWERING 2 of Arabidopsis thaliana. Plant Mol Biol. 2013;81:57–69.

    Article  CAS  PubMed  Google Scholar 

  23. Seo PJ, Park MJ, Lim MH, Kim SG, Lee M, Baldwin IT, et al. A self-regulatory circuit of CIRCADIAN CLOCK-ASSOCIATED1 underlies the circadian clock regulation of temperature responses in Arabidopsis. Plant Cell. 2012;24:2427–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Staiger D, Brown JWS. Alternative splicing at the intersection of biological timing, development, and stress responses. Plant Cell. 2013;25:3640–56.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Pei YX, Niu LF, Lu FL, Liu CY, Zhai JX, Kong XF, et al. Mutations in the type II protein arginine methyltransferase AtPRMT5 result in pleiotropic developmental defects in Arabidopsis(1[C][OA]). Plant Physiol. 2007;144:1913–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Verhage L, Severing EI, Bucher J, Lammers M, Busscher-Lange J, Bonnema G, et al. Splicing-related genes are alternatively spliced upon changes in ambient temperatures in plants. PLoS One. 2017;12:e0172950.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Naftelberg S, Schor IE, Ast G, Kornblihtt AR. Regulation of alternative splicing through coupling with transcription and chromatin structure. Annu Rev Biochem. 2015;84:165–98.

    Article  CAS  PubMed  Google Scholar 

  28. Chen M, Manley JL. Mechanisms of alternative splicing regulation: insights from molecular and genomics approaches. Nat Rev Mol Cell Biol. 2009;10:741–54.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Zhou HL, Luo G, Wise JA, Lou H. Regulation of alternative splicing by local histone modifications: potential roles for RNA-guided mechanisms. Nucleic Acids Res. 2014;42:701–13.

    Article  CAS  PubMed  Google Scholar 

  30. Brown SJ, Stoilov P, Xing Y. Chromatin and epigenetic regulation of pre-mRNA processing. Hum Mol Genet. 2012;21:R90–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Luco RF, Allo M, Schor IE, Kornblihtt AR, Misteli T. Epigenetics in alternative pre-mRNA splicing. Cell. 2011;144:16–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Luco RF, Pan Q, Tominaga K, Blencowe BJ, Pereira-Smith OM, Misteli T. Regulation of alternative splicing by histone modifications. Science. 2010;327:996–1000.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Roudier F, Ahmed I, Berard C, Sarazin A, Mary-Huard T, Cortijo S, et al. Integrative epigenomic mapping defines four main chromatin states in Arabidopsis. Embo J. 2011;30:1928–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Airoldi CA, Mckay M, Davies B. MAF2 is regulated by temperature-dependent splicing and represses flowering at low temperatures in parallel with FLM. PLoS One. 2015;10:e0126516.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Sureshkumar S, Dent C, Seleznev A, Tasset C, Balasubramanian S. Nonsense-mediated mRNA decay modulates FLM-dependent thermosensory flowering response in Arabidopsis. Nature Plants. 2016;2:16055.

    Article  CAS  PubMed  Google Scholar 

  36. Xu ML, Hu TQ, Zhao JF, Park MY, Earley KW, Wu G, et al. Developmental functions of miR156-regulated SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes in Arabidopsis thaliana. PLoS Genet. 2016;12:e1006263.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Bao ZL, Zhang N, Hua J. Endopolyploidization and flowering time are antagonistically regulated by checkpoint component MAD1 and immunity modulator MOS1. Nat Commun. 2014;5:5628.

    Article  CAS  PubMed  Google Scholar 

  38. Lim PO, Lee IC, Kim J, Kim HJ, Ryu JS, Woo HR, et al. Auxin response factor 2 (ARF2) plays a major role in regulating auxin-mediated leaf longevity. J Exp Bot. 2010;61:1419–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Friml J, Vieten A, Sauer M, Weijers D, Schwarz H, Hamann T, et al. Efflux-dependent auxin gradients establish the apical-basal axis of Arabidopsis. Nature. 2003;426:147–53.

    Article  CAS  PubMed  Google Scholar 

  40. Kim DH, Sung S. The Plant Homeo Domain finger protein, VIN3-LIKE 2, is necessary for photoperiod-mediated epigenetic regulation of the floral repressor, MAF5. Proc Natl Acad Sci U S A. 2010;107:17029–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Bezhani S, Winter C, Hershman S, Wagner JD, Kennedy JF, Kwon CS, et al. Unique, shared, and redundant roles for the Arabidopsis SWI/SNF chromatin remodeling ATPases BRAHMA and SPLAYED. Plant Cell. 2007;19:403–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Kwon CS, Chen CB, Wagner D. WUSCHEL is a primary target for transcriptional regulation by SPLAYED in dynamic control of stem cell fate in Arabidopsis. Genes Dev. 2005;19:992–1003.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Vercruyssen L, Verkest A, Gonzalez N, Heyndrickx KS, Eeckhout D, Han SK, et al. ANGUSTIFOLIA3 binds to SWI/SNF chromatin remodeling complexes to regulate transcription during arabidopsis leaf development. Plant Cell. 2014;26:210–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Berr A, Shafiq S, Pinon V, Dong AW, Shen WH. The trxG family histone methyltransferase SET DOMAIN GROUP 26 promotes flowering via a distinctive genetic pathway. Plant J. 2015;81:316–28.

    Article  CAS  PubMed  Google Scholar 

  45. Zhao Z, Yu Y, Meyer D, Wu CJ, Shen WH. Prevention of early flowering by expression of FLOWERING LOCUS C requires methylation of histone H3K36. Nat Cell Biol. 2005;7:1256–60.

    Article  PubMed  Google Scholar 

  46. Pajoro A, Verhage L, Immink RGH. Plasticity versus adaptation of ambient-temperature flowering response. Trends Plant Sci. 2016;21:6–8.

    Article  CAS  PubMed  Google Scholar 

  47. Lu SX, Knowles SM, Webb CJ, Celaya RB, Cha CA, Siu JP, et al. The Jumonji C domain-containing protein JMJ30 regulates period length in the Arabidopsis circadian clock. Plant Physiol. 2011;155:906–15.

    Article  CAS  PubMed  Google Scholar 

  48. Xu YF, Gan ES, Zhou J, Wee WY, Zhang XY, Ito T. Arabidopsis MRG domain proteins bridge two histone modifications to elevate expression of flowering genes. Nucleic Acids Res. 2014;42:10960–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Ruhl C, Stauffer E, Kahles A, Wagner G, Drechsel G, Ratsch G, et al. Polypyrimidine tract binding protein homologs from Arabidopsis are key regulators of alternative splicing with implications in fundamental developmental processes. Plant Cell. 2012;24:4360–75.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Ali GS, Palusa SG, Golovkin M, Prasad J, Manley JL, Reddy ASN. Regulation of plant developmental processes by a novel splicing factor. PLoS One. 2007;2:e471.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Arciga-Reyes L, Wootton L, Kieffer M, Davies B. UPF1 is required for nonsense-mediated mRNA decay (NMD) and RNAi in Arabidopsis. Plant J. 2006;47:480–9.

    Article  CAS  PubMed  Google Scholar 

  52. Kalyna M, Lopato S, Barta A. Ectopic expression of atRSZ33 reveals its function in splicing and causes pleiotropic changes in development. Mol Biol Cell. 2003;14:3565–77.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Lopato S, Kalyna M, Dorner S, Kobayashi R, Krainer AR, Barta A. atSRp30, one of two SF2/ASF-like proteins from Arabidopsis thaliana, regulates splicing of specific plant genes. Genes Dev. 1999;13:987–1001.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Swaraz AM, Park YD, Hur Y. Knock-out mutations of Arabidopsis SmD3-b induce pleotropic phenotypes through altered transcript splicing. Plant Sci. 2011;180:661–71.

    Article  CAS  PubMed  Google Scholar 

  55. Shaul O. Unique aspects of plant nonsense-mediated mRNA decay. Trends Plant Sci. 2015;20:767–79.

    Article  CAS  PubMed  Google Scholar 

  56. Oesterreich FC, Herzel L, Straube K, Hujer K, Howard J, Neugebauer KM. Splicing of nascent RNA coincides with intron exit from RNA polymerase II. Cell. 2016;165:372–81.

    Article  PubMed Central  Google Scholar 

  57. Sidaway-Lee K, Costa MJ, Rand DA, Finkenstadt B, Penfield S. Direct measurement of transcription rates reveals multiple mechanisms for configuration of the Arabidopsis ambient temperature response. Genome Biol. 2014;15:R45.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Kumar SV, Wigge PA. H2A.Z-containing nucleosomes mediate the thermosensory response in Arabidopsis. Cell. 2010;140:136–47.

    Article  CAS  PubMed  Google Scholar 

  59. Neves LT, Douglass S, Spreafico R, Venkataramanan S, Kress TL, Johnson TL. The histone variant H2A.Z promotes efficient cotranscriptional splicing in S. cerevisiae. Genes Dev. 2017;31:702–17.

    Article  CAS  PubMed  Google Scholar 

  60. Nissen KE, Homer CM, Ryan CJ, Shales M, Krogan NJ, Patrick KL, et al. The histone variant H2A.Z promotes splicing of weak introns. Genes Dev. 2017;31:688–701.

    Article  CAS  PubMed  Google Scholar 

  61. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–U174.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Katz Y, Wang ET, Airoldi EM, Burge CB. Analysis and design of RNA sequencing experiments for identifying isoform regulation. Nat Methods. 2010;7:1009–U1101.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Anders S, Pyl PT, Huber W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.

    Article  CAS  PubMed  Google Scholar 

  65. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Maere S, Heymans K, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of Gene Ontology categories in Biological Networks. Bioinformatics. 2005;21:3448–9.

    Article  CAS  PubMed  Google Scholar 

  67. Blake JA, Christie KR, Dolan ME, Drabkin HJ, Hill DP, Ni L, et al. Gene Ontology Consortium: going forward. Nucleic Acids Res. 2015;43:D1049–56.

    Article  Google Scholar 

  68. Kaufmann K, Muino JM, Osteras M, Farinelli L, Krajewski P, Angenent GC. Chromatin immunoprecipitation (ChIP) of plant transcription factors followed by sequencing (ChIP-SEQ) or hybridization to whole genome arrays (ChIP-CHIP). Nat Protoc. 2010;5:457–72.

    Article  CAS  PubMed  Google Scholar 

  69. Luo CY, Sidote DJ, Zhang Y, Kerstetter RA, Michael TP, Lam E. Integrative analysis of chromatin states in Arabidopsis identified potential regulatory mechanisms for natural antisense transcript production. Plant J. 2013;73:77–90.

    Article  CAS  PubMed  Google Scholar 

  70. Yang HC, Howard M, Dean C. Antagonistic roles for H3K36me3 and H3K27me3 in the cold-induced epigenetic switch at Arabidopsis FLC. Curr Biol. 2014;24:1793–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Ramirez F, Dundar F, Diehl S, Gruning BA, Manke T. deepTools: a flexible platform for exploring deep-sequencing data. Nucleic Acids Res. 2014;42:W187–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Chen K, Xi Y, Pan X, Li Z, Kaestner K, Tyler J, et al. DANPOS: dynamic analysis of nucleosome position and occupancy by sequencing. Genome Res. 2013;23:341–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

We are grateful to B. Davies and B. Thomma for critical comments on the manuscript. We thank L. Faino for sharing R scripts and T. Ito for the mrg1-1 mrg2-3 seeds.

Funding

This work was supported by the NWO grant (849.13.005) for the project ERACAPS13.012, FLOWPLAST. The RNA-seq experiments were supported by a ZonMw Enabling Technology Hotel project (40-43500-98-054).

Availability of data and materials

Sequence data are available in the Gene Expression Omnibus database (GEO accession GSE85282).

Authors’ contributions

AP and RGHI conceived and designed the study with help from GA. AP performed the experiments. AP and ES analyzed and interpreted data. AP, GA, and RGHI wrote the paper. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to R. G. H. Immink.

Additional files

Additional file 1: Table S1.

List of significant (see Methods) DiS events in at least one 25 °C time point (day 1, 3, or 5) compared to the 16 °C. (XLSX 75 kb)

Additional file 2: Supplementary Figure S1-S7.

Temperature induced AS events. Figure S2. H3K36me3 dynamics upon the temperature change. Figure S3. Quantitative change in H3K36me3 (ChIP-seq output) at the same location of a DiS event (RNA-seq output) for the circadian clock genes PRR3 and PRR7. Figure S4. Temperature-induced AS in sdg8-2 and sdg26-1. Figure S5. Changes in H3K36me3 levels and mrg1 mrg2 affect temperature induction of flowering. Figure S6. Changes in H3K36me3 levels do not cause pleotropic effects on plant development. Figure S7. Pearson correlation of ChIP-seq experiments. (PDF 1186 kb)

Additional file 3: Table S2.

DEGs between 16 °C and 25 °C day 1 (sheet1), 16 °C and 25 °C day 3 (sheet2), 16 °C and 25 °C day 5 (sheet3). DEGs are defined as genes with a change in expression of log2 Fold Change |1| and adjusted p value according to the BH method for controlling false discovery < 0.05. (XLSX 518 kb)

Additional file 4: Table S3.

H3K36me3 modified regions at 16 °C or 25 °C day 1 and quantitative comparison obtained with DANPOS2. (XLSX 15623 kb)

Additional file 5: Table S4.

List of significant (see “Methods”) DiS events between 16 °C and 25 °C day 1 in at least one mutant genotype. (XLSX 5555 kb)

Additional file 6: Table S5.

List of significant (see “Methods”) DiS events between WT and sdg8-2 and events between WT and sdg26-1 at 16 °C. (XLSX 17 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pajoro, A., Severing, E., Angenent, G.C. et al. Histone H3 lysine 36 methylation affects temperature-induced alternative splicing and flowering in plants. Genome Biol 18, 102 (2017). https://doi.org/10.1186/s13059-017-1235-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13059-017-1235-x

Keywords