Cold stress is one of the major environmental factors that limit growth and yield of plants. However, it is still not fully understood how plants account for daily temperature fluctuations, nor how these temperature changes are integrated with other regulatory systems such as the circadian clock. We demonstrate that REVEILLE2 undergoes alternative splicing after chilling that increases accumulation of a transcript isoform encoding a MYB-like transcription factor. We explore the biological function of REVEILLE2 in Arabidopsis thaliana using a combination of molecular genetics, transcriptomics, and physiology. Disruption of REVEILLE2 alternative splicing alters regulatory gene expression, impairs circadian timing, and improves photosynthetic capacity. Changes in nuclear gene expression are particularly apparent in the initial hours following chilling, with chloroplast gene expression subsequently upregulated. The response of REVEILLE2 to chilling extends our understanding of plants immediate response to cooling. We propose that the circadian component REVEILLE2 restricts plants responses to nocturnal reductions in temperature, thereby enabling appropriate responses to daily environmental changes.

RNA-sequencing (RNA-seq) is currently the method of choice for analysis of differential gene expression. To fully exploit the wealth of data generated from genome-wide transcriptomic approaches, the initial design of the experiment is of paramount importance. Biological rhythms in nature are pervasive and are driven by endogenous gene networks collectively known as circadian clocks. Measuring circadian gene expression requires time-course experiments which take into account time-of-day factors influencing variability in expression levels. We describe here an approach for characterizing diurnal changes in expression and alternative splicing for plants undergoing cooling. The method uses inexpensive everyday laboratory equipment and utilizes an RNA-seq application (3D RNA-seq) that can handle complex experimental designs and requires little or no prior bioinformatics expertise.

The core of the plant circadian clock involves multiple interlocking gene expression loops and post-translational controls along with inputs from light and metabolism. The complexity of the interactions is such that few specific functions can be ascribed to single components. In previous work, we reported differences in the operation of the clocks in Arabidopsis shoots and roots, including the effects of mutations of key clock components. Here, we have used luciferase imaging to study prr7 mutants expressing CCA1::LUC and GI::LUC markers. In mature shoots expressing CCA1::LUC, loss of PRR7 radically altered behaviour in light:dark cycles and caused loss of rhythmicity in constant light but had little effect on roots. In contrast, in mature plants expressing GI::LUC, loss of PRR7 had little effect in light:dark cycles but in constant light increased the circadian period in shoots and reduced it in roots. We conclude that most or all of the circadian input to the CCA1 promoter in shoots is mediated by PRR7 and that loss of PRR7 has organ-specific effects. The results emphasise the differences in operation of the shoot and root clocks, and the importance of studying clock mutants in both light:dark cycles and constant light.

The circadian clock regulates the timing of many aspects of plant physiology, and this requires entrainment of the clock to the prevailing day:night cycle. Different plant cells and tissues can oscillate with different free-running periods, so coordination of timing across the plant is crucial. Previous work showed that a major difference between the clock in mature shoots and roots involves light inputs. The objective of this work was to define, in Arabidopsis thaliana, the operation of the root clock in more detail, and in particular how it responds to light quality. Luciferase imaging was used to study the shoot and root clocks in several null mutants of clock components and in lines with aberrant expression of phytochromes. Mutations in each of the components of the evening complex (EARLY FLOWERING 3 and 4, and LUX ARRHYTHMO) were found to have specific effects on roots, by affecting either rhythmicity or period and its response to light quality. The data suggest that the evening complex is a key part of the light input mechanism that differs between shoots and roots and show that roots sense red light via phytochrome B.

Plants re-program their gene expression when responding to changing environmental conditions. Besides differential gene expression, extensive alternative splicing (AS) of pre-mRNAs and changes in expression of long non-coding RNAs (lncRNAs) are associated with stress responses. RNA-sequencing of a diel time-series of the initial response of Arabidopsis thaliana rosettes to low temperature showed massive and rapid waves of both transcriptional and AS activity in protein-coding genes. We exploited the high diversity of transcript isoforms in AtRTD2 to examine regulation and post-transcriptional regulation of lncRNA gene expression in response to cold stress. We identified 135 lncRNA genes with cold-dependent differential expression (DE) and/or differential alternative splicing (DAS) of lncRNAs including natural antisense RNAs, sORF lncRNAs, and precursors of microRNAs (miRNAs) and trans-acting small-interfering RNAs (tasiRNAs). The high resolution (HR) of the time-series allowed the dynamics of changes in transcription and AS to be determined and identified early and adaptive transcriptional and AS changes in the cold response. Some lncRNA genes were regulated only at the level of AS and using plants grown at different temperatures and a HR time-course of the first 3 h of temperature reduction, we demonstrated that the AS of some lncRNAs is highly sensitive to small temperature changes suggesting tight regulation of expression. In particular, a splicing event in TAS1a which removed an intron that contained the miR173 processing and phased siRNAs generation sites was differentially alternatively spliced in response to cold. The cold-induced reduction of the spliced form of TAS1a and of the tasiRNAs suggests that splicing may enhance production of the siRNAs. Our results identify candidate lncRNAs that may contribute to the regulation of expression that determines the physiological processes essential for acclimation and freezing tolerance.

How plants perceive and respond to temperature remains an important question in the plant sciences. Temperature perception and signal transduction may occur through temperature-sensitive intramolecular folding of primary mRNA transcripts. Recent studies suggested a role for retention of the first intron in the 5'UTR of the clock component LATE ELONGATED HYPOCOTYL (LHY) in response to changes in temperature. Here, we identified a set of haplotypes in the LHY 5'UTR, examined their global spatial distribution, and obtained evidence that haplotype can affect temperature-dependent splicing of LHY transcripts. Correlations of haplotype spatial distributions with global bioclimatic variables and altitude point to associations with annual mean temperature and temperature fluctuation. Relatively rare relict type accessions correlate with lower mean temperature and greater temperature fluctuation and the spatial distribution of other haplotypes may be informative of evolutionary processes driving colonization of ecosystems. We propose that haplotypes may possess distinct 5'UTR pre-mRNA folding thermodynamics and/or specific biological stabilities based around the binding of trans-acting RNA splicing factors, a consequence of which is scalable splicing sensitivity of a central clock component that is likely tuned to specific temperature environments.

AbstractOne of the ways in which plants can respond to temperature is via alternative splicing (AS). Previous work showed that temperature changes affected the splicing of several circadian clock gene transcripts. Here, we investigated the role of RNA‐binding splicing factors (SFs) in temperature‐sensitive AS of the clock gene LATE ELONGATED HYPOCOTYL (LHY). We characterized, in wild type plants, temperature‐associated isoform switching and expression patterns for SF transcripts from a high‐resolution temperature and time series RNA‐seq experiment. In addition, we employed quantitative RT‐PCR of SF mutant plants to explore the role of the SFs in cooling‐associated AS of LHY. We show that the splicing and expression of several SFs responds sufficiently, rapidly, and sensitively to temperature changes to contribute to the splicing of the 5′UTR of LHY. Moreover, the choice of splice site in LHY was altered in some SF mutants. The splicing of the 5′UTR region of LHY has characteristics of a molecular thermostat, where the ratio of transcript isoforms is sensitive to temperature changes as modest as 2 °C and is scalable over a wide dynamic range of temperature. Our work provides novel insight into SF‐mediated coupling of the perception of temperature to post‐transcriptional regulation of the clock.

Plants have adapted to tolerate and survive constantly changing environmental conditions by reprogramming gene expression The dynamics of the contribution of alternative splicing (AS) to stress responses are unknown. RNA-sequencing of a time-series of Arabidopsis thaliana plants exposed to cold determines the timing of significant AS changes. This shows a massive and rapid AS response with coincident waves of transcriptional and AS activity occurring in the first few hours of temperature reduction and further AS throughout the cold. In particular, hundreds of genes showed changes in expression due to rapidly occurring AS in response to cold ("early AS" genes); these included numerous novel cold-responsive transcription factors and splicing factors/RNA binding proteins regulated only by AS. The speed and sensitivity to small temperature changes of AS of some of these genes suggest that fine-tuning expression via AS pathways contributes to the thermo-plasticity of expression. Four early AS splicing regulatory genes have been shown previously to be required for freezing tolerance and acclimation; we provide evidence of a fifth gene, U2B"-LIKE Such factors likely drive cascades of AS of downstream genes that, alongside transcription, modulate transcriptome reprogramming that together govern the physiological and survival responses of plants to low temperature.

The full-text may be used and/or reproduced, and given to third parties in any format or medium, without prior permission or charge, for personal research or study, educational, or not-for-pro t purposes provided that: • a full bibliographic reference is made to the original source • a link is made to the metadata record in DRO • the full-text is not changed in any way The full-text must not be sold in any format or medium without the formal permission of the copyright holders. Please consult the full DRO policy for further details.

Correct operation of the plant circadian clock is crucial for optimal growth and development. Recent evidence has shown that the plant clock is tissue specific and potentially hierarchical, implying that there are signalling mechanisms that can synchronise the clock in different tissues. Here, I have addressed the mechanism that allows the shoot and root clocks to be synchronised in light:dark cycles but not in continuous light. Luciferase imaging data from 2 different Arabidopsis accessions with 2 different markers show that the period of the root clock is much less sensitive to blue light than to red light. Decapitated roots were imaged either in darkness or with the top section of root tissue exposed to light. Exposure to red light reduced the period of the root tissue maintained in darkness, whereas exposure to blue light did not. The data indicate that light can be piped through root tissue to affect the circadian period of tissue in darkness. I propose that the synchronisation of shoots and roots in light:dark cycles is achieved by light piping from shoots to roots.

Alternative splicing generates multiple transcript and protein isoforms from the same gene and thus is important in gene expression regulation. To date, RNA-sequencing (RNA-seq) is the standard method for quantifying changes in alternative splicing on a genome-wide scale. Understanding the current limitations of RNA-seq is crucial for reliable analysis and the lack of high quality, comprehensive transcriptomes for most species, including model organisms such as Arabidopsis, is a major constraint in accurate quantification of transcript isoforms. To address this, we designed a novel pipeline with stringent filters and assembled a comprehensive Reference Transcript Dataset for Arabidopsis (AtRTD2) containing 82,190 non-redundant transcripts from 34 212 genes. Extensive experimental validation showed that AtRTD2 and its modified version, AtRTD2-QUASI, for use in Quantification of Alternatively Spliced Isoforms, outperform other available transcriptomes in RNA-seq analysis. This strategy can be implemented in other species to build a pipeline for transcript-level expression and alternative splicing analyses.

Summary Circadian clocks allow the temporal compartmentalization of biological processes. In Arabidopsis, circadian rhythms display organ specificity but the underlying molecular causes have not been identified. We investigated the mechanisms responsible for the similarities and differences between the clocks of mature shoots and roots in constant conditions and in light : dark cycles. We developed an imaging system to monitor clock gene expression in shoots and light‐ or dark‐grown roots, modified a recent mathematical model of the Arabidopsis clock and used this to simulate our new data. We showed that the shoot and root circadian clocks have different rhythmic properties (period and amplitude) and respond differently to light quality. The root clock was entrained by direct exposure to low‐intensity light, even in antiphase to the illumination of shoots. Differences between the clocks were more pronounced in conditions where light was present than in constant darkness, and persisted in the presence of sucrose. We simulated the data successfully by modifying those parameters of a clock model that are related to light inputs. We conclude that differences and similarities between the shoot and root clocks can largely be explained by organ‐specific light inputs. This provides mechanistic insight into the developing field of organ‐specific clocks.

Background Alternative splicing is the major post-transcriptional mechanism by which gene expression is regulated and affects a wide range of processes and responses in most eukaryotic organisms. RNA-sequencing (RNA-seq) can generate genome-wide quantification of individual transcript isoforms to identify changes in expression and alternative splicing. RNA-seq is an essential modern tool but its ability to accurately quantify transcript isoforms depends on the diversity, completeness and quality of the transcript information. Results We have developed a new Reference Transcript Dataset for Arabidopsis (AtRTD2) for RNA-seq analysis containing over 82k non-redundant transcripts, whereby 74,194 transcripts originate from 27,667 protein-coding genes. A total of 13,524 protein-coding genes have at least one alternatively spliced transcript in AtRTD2 such that about 60% of the 22,453 protein-coding, intron-containing genes in Arabidopsis undergo alternative splicing. More than 600 putative U12 introns were identified in more than 2,000 transcripts. AtRTD2 was generated from transcript assemblies of ca. 8.5 billion pairs of reads from 285 RNA-seq data sets obtained from 129 RNA-seq libraries and merged along with the previous version, AtRTD, and Araport11 transcript assemblies. AtRTD2 increases the diversity of transcripts and through application of stringent filters represents the most extensive and accurate transcript collection for Arabidopsis to date. We have demonstrated a generally good correlation of alternative splicing ratios from RNA-seq data analysed by Salmon and experimental data from high resolution RT-PCR. However, we have observed inaccurate quantification of transcript isoforms for genes with multiple transcripts which have variation in the lengths of their UTRs. This variation is not effectively corrected in RNA-seq analysis programmes and will therefore impact RNA-seq analyses generally. To address this, we have tested different genome-wide modifications of AtRTD2 to improve transcript quantification and alternative splicing analysis. As a result, we release AtRTD2-QUASI specifically for use in Quantification of Alternatively Spliced Isoforms and demonstrate that it out-performs other available transcriptomes for RNA-seq analysis. Conclusions We have generated a new transcriptome resource for RNA-seq analyses in Arabidopsis (AtRTD2) designed to address quantification of different isoforms and alternative splicing in gene expression studies. Experimental validation of alternative splicing changes identified inaccuracies in transcript quantification due to UTR length variation. To solve this problem, we also release a modified reference transcriptome, AtRTD2-QUASI for quantification of transcript isoforms, which shows high correlation with experimental data.

Summary RNA‐sequencing (RNA‐seq) allows global gene expression analysis at the individual transcript level. Accurate quantification of transcript variants generated by alternative splicing (AS) remains a challenge. We have developed a comprehensive, nonredundant Arabidopsis reference transcript dataset (AtRTD) containing over 74 000 transcripts for use with algorithms to quantify AS transcript isoforms in RNA‐seq. The AtRTD was formed by merging transcripts from TAIR10 and novel transcripts identified in an AS discovery project. We have estimated transcript abundance in RNA‐seq data using the transcriptome‐based alignment‐free programmes Sailfish and Salmon and have validated quantification of splicing ratios from RNA‐seq by high resolution reverse transcription polymerase chain reaction (HR RT‐PCR). Good correlations between splicing ratios from RNA‐seq and HR RT‐PCR were obtained demonstrating the accuracy of abundances calculated for individual transcripts in RNA‐seq. The AtRTD is a resource that will have immediate utility in analysing Arabidopsis RNA‐seq data to quantify differential transcript abundance and expression.

In the March 2012 issue of The Plant Cell we describe extensive alternative splicing (AS) of Arabidopsis circadian clock genes. Notably these distinct post-transcriptional events associate with different steady-state temperatures and also with plants undergoing temperature transitions leading us to propose that temperature-associated AS is an additional mechanism involved in the operation and control of the plant circadian clock. Here we show that temperature associated AS also extends to REVEILLE 8 (RVE8), demonstrating a hitherto unrecognized link between the expression of this clock associated gene and temperature. Finally we discuss our observations of the plastic nature of clock gene expression at the post-transcriptional level in the context of the ongoing fascination of how plants respond to temperature.

Alternative splicing plays crucial roles by influencing the diversity of the transcriptome and proteome and regulating protein structure/function and gene expression. It is widespread in plants, and alteration of the levels of splicing factors leads to a wide variety of growth and developmental phenotypes. The circadian clock is a complex piece of cellular machinery that can regulate physiology and behavior to anticipate predictable environmental changes on a revolving planet. We have performed a system-wide analysis of alternative splicing in clock components in Arabidopsis thaliana plants acclimated to different steady state temperatures or undergoing temperature transitions. This revealed extensive alternative splicing in clock genes and dynamic changes in alternatively spliced transcripts. Several of these changes, notably those affecting the circadian clock genes late elongated hypocotyl (LHY) and pseudo response regulator7, are temperature-dependent and contribute markedly to functionally important changes in clock gene expression in temperature transitions by producing nonfunctional transcripts and/or inducing nonsense-mediated decay. Temperature effects on alternative splicing contribute to a decline in LHY transcript abundance on cooling, but LHY promoter strength is not affected. We propose that temperature-associated alternative splicing is an additional mechanism involved in the operation and regulation of the plant circadian clock.

The circadian oscillator in eukaryotes consists of several interlocking feedback loops through which the expression of clock genes is controlled. It is generally assumed that all plant cells contain essentially identical and cell- autonomous multiloop clocks. Here, we show that the circadian clock in the roots of mature Arabidopsis plants differs markedly from that in the shoots and that the root clock is synchronized by a photosynthesis- related signal from the shoot. Two of the feedback loops of the plant circadian clock are disengaged in roots, because two key clock components, the transcription factors CCA1 and LHY, are able to inhibit gene expression in shoots but not in roots. Thus, the plant clock is organ- specific but not organ- autonomous.