Genome-scale analysis of Arabidopsis splicing-related protein kinase families reveals roles in abiotic stress adaptation

Nearly 60 - 80 % of intron-containing plant genes undergo alternative splicing in response to either stress or plant developmental cues. RNA splicing is performed by a large ribonucleoprotein complex called the spliceosome in conjunction with associated subunits such as serine arginine (SR) proteins, all of which undergo extensive phosphorylation. In plants, there are three main protein kinase families suggested to phosphorylate core spliceosome subunits and related splicing factors based on orthology to human splicing-related kinases: the SERINE/ARGININE PROTEIN KINASES (SRPK), ARABIDOPSIS FUS3 COMPLEMENT (AFC), and Pre-mRNA PROCESSING FACTOR 4 (PRP4K) protein kinases. To better define the conservation and role(s) of these kinases in plants, we performed a genome-scale analysis of the three families across photosynthetic eukaryotes, followed by extensive transcriptomic and bioinformatic analysis of all Arabidopsis thaliana SRPK, AFC, and PRP4K protein kinases to elucidate their biological functions. Unexpectedly, this revealed the existence of SRPK and AFC phylogenetic groups with distinct promoter elements and patterns of transcriptional response to abiotic stress, while PRP4Ks possess no phylogenetic sub-divisions, suggestive of functional redundancy. We also reveal splicing-related kinase families are both diel and photoperiod regulated, implicating different orthologs as discrete time-of-day RNA splicing regulators. This foundational work establishes a number of new hypotheses regarding how reversible spliceosome phosphorylation contributes to both diel plant cell regulation and abiotic stress adaptation in plants.

Closing the protein gap in plant chronobiology

Our modern understanding of diel cell regulation in plants stems from foundational work in the late 1990s that analysed the dynamics of selected genes and mutants in Arabidopsis thaliana. The subsequent rise of transcriptomics technologies such as microarrays and RNA sequencing has substantially increased our understanding of anticipatory (circadian) and reactive (light- or dark-triggered) diel events in plants. However, it is also becoming clear that gene expression data fail to capture critical events in diel regulation that can only be explained by studying protein-level dynamics. Over the past decade, mass spectrometry technologies and quantitative proteomic workflows have significantly advanced, finally allowing scientists to characterise diel protein regulation at high throughput. Initial proteomic investigations suggest that the diel transcriptome and proteome generally lack synchrony and that the timing of daily regulatory events in plants is impacted by multiple levels of protein regulation (e.g., post-translational modifications [PTMs] and protein-protein interactions [PPIs]). Here, we highlight and summarise how the use of quantitative proteomics to elucidate diel plant cell regulation has advanced our understanding of these processes. We argue that this new understanding, coupled with the extraordinary developments in mass spectrometry technologies, demands greater focus on protein-level regulation of, and by, the circadian clock. This includes hitherto unexplored diel dynamics of protein turnover, PTMs, protein subcellular localisation and PPIs that can be masked by simple transcript- and protein-level changes. Finally, we propose new directions for how the latest advancements in quantitative proteomics can be utilised to answer outstanding questions in plant chronobiology.

Global transcriptome analysis reveals circadian control of splicing events in Arabidopsis thaliana

The circadian clock of Arabidopsis thaliana controls many physiological and molecular processes, allowing plants to anticipate daily changes in their environment. However, developing a detailed understanding of how oscillations in mRNA levels are connected to oscillations in co/post-transcriptional processes, such as splicing, has remained a challenge. Here we applied a combined approach using deep transcriptome sequencing and bioinformatics tools to identify novel circadian-regulated genes and splicing events. Using a stringent approach, we identified 300 intron retention, eight exon skipping, 79 alternative 3' splice site usage, 48 alternative 5' splice site usage, and 350 multiple (more than one event type) annotated events under circadian regulation. We also found seven and 721 novel alternative exonic and intronic events. Depletion of the circadian-regulated splicing factor AtSPF30 homologue resulted in the disruption of a subset of clock-controlled splicing events. Altogether, our global circadian RNA-seq coupled with an in silico, event-centred, splicing analysis tool offers a new approach for studying the interplay between the circadian clock and the splicing machinery at a global scale. The identification of many circadian-regulated splicing events broadens our current understanding of the level of control that the circadian clock has over this co/post-transcriptional regulatory layer.

The Arabidopsis LDL1/2-HDA6 histone modification complex is functionally associated with CCA1/LHY in regulation of circadian clock genes

In Arabidopsis, the circadian clock central oscillator genes are important cellular components to generate and maintain circadian rhythms. There is a negative feedback loop between the morning expressed CCA1 (CIRCADIAN CLOCK ASSOCIATED 1)/LHY (LATE ELONGATED HYPOCOTYL) and evening expressed TOC1 (TIMING OF CAB EXPRESSION 1). CCA1 and LHY negatively regulate the expression of TOC1, while TOC1 also binds to the promoters of CCA1 and LHY to repress their expression. Recent studies indicate that histone modifications play an important role in the regulation of the central oscillators. However, the regulatory relationship between histone modifications and the circadian clock genes remains largely unclear. In this study, we found that the Lysine-Specific Demethylase 1 (LSD1)-like histone demethylases, LDL1 and LDL2, can interact with CCA1/LHY to repress the expression of TOC1. ChIP-Seq analysis indicated that LDL1 targets a subset of genes involved in the circadian rhythm regulated by CCA1. Furthermore, LDL1 and LDL2 interact with the histone deacetylase HDA6 and co-regulate TOC1 by histone demetylation and deacetylaion. These results provide new insight into the molecular mechanism of how the circadian clock central oscillator genes are regulated through histone modifications.

Spliceosome factors target timeless (tim) mRNA to control clock protein accumulation and circadian behavior in Drosophila

Transcription-translation feedback loops that comprise eukaryotic circadian clocks rely upon temporal delays that separate the phase of active transcription of clock genes, such as Drosophila period (per) and timeless (tim), from negative feedback by the two proteins. However, our understanding of the mechanisms involved is incomplete. Through an RNA interference screen, we found that pre-mRNA processing 4 (PRP4) kinase, a component of the U4/U5.U6 triple small nuclear ribonucleoprotein (tri-snRNP) spliceosome, and other tri-snRNP components regulate cycling of the molecular clock as well as rest:activity rhythms. Unbiased RNA-Sequencing uncovered an alternatively spliced intron in tim whose increased retention upon prp4 downregulation leads to decreased TIM levels. We demonstrate that the splicing of tim is rhythmic with a phase that parallels delayed accumulation of the protein in a 24 hr cycle. We propose that alternative splicing constitutes an important clock mechanism for delaying the daily accumulation of clock proteins, and thereby negative feedback by them.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

From network to phenotype: the dynamic wiring of an Arabidopsis transcriptional network induced by osmotic stress

Plants have established different mechanisms to cope with environmental fluctuations and accordingly fine-tune their growth and development through the regulation of complex molecular networks. It is largely unknown how the network architectures change and what the key regulators in stress responses and plant growth are. Here, we investigated a complex, highly interconnected network of 20 Arabidopsis transcription factors (TFs) at the basis of leaf growth inhibition upon mild osmotic stress. We tracked the dynamic behavior of the stress-responsive TFs over time, showing the rapid induction following stress treatment, specifically in growing leaves. The connections between the TFs were uncovered using inducible overexpression lines and were validated with transient expression assays. This study resulted in the identification of a core network, composed of ERF6, ERF8, ERF9, ERF59, and ERF98, which is responsible for most transcriptional connections. The analyses highlight the biological function of this core network in environmental adaptation and its redundancy. Finally, a phenotypic analysis of loss-of-function and gain-of-function lines of the transcription factors established multiple connections between the stress-responsive network and leaf growth.

The role of alternative Polyadenylation in regulation of rhythmic gene expression

Background

Alternative transcription is common in eukaryotic cells and plays important role in regulation of cellular processes. Alternative polyadenylation results from ambiguous PolyA signals in 3′ untranslated region (UTR) of a gene. Such alternative transcripts share the same coding part, but differ by a stretch of UTR that may contain important functional sites.

Methods

The methodoogy of this study is based on mathematical modeling, analytical solution, and subsequent validation by datamining in multiple independent experimental data from previously published studies.

Results

In this study we propose a mathematical model that describes the population dynamics of alternatively polyadenylated transcripts in conjunction with rhythmic expression such as transcription oscillation driven by circadian or metabolic oscillators. Analysis of the model shows that alternative transcripts with different turnover rates acquire a phase shift if the transcript decay rate is different. Difference in decay rate is one of the consequences of alternative polyadenylation. Phase shift can reach values equal to half the period of oscillation, which makes alternative transcripts oscillate in abundance in counter-phase to each other. Since counter-phased transcripts share the coding part, the rate of translation becomes constant. We have analyzed a few data sets collected in circadian timeline for the occurrence of transcript behavior that fits the mathematical model.

Conclusion

Alternative transcripts with different turnover rate create the effect of rectifier. This “molecular diode” moderates or completely eliminates oscillation of individual transcripts and stabilizes overall protein production rate. In our observation this phenomenon is very common in different tissues in plants, mice, and humans. The occurrence of counter-phased alternative transcripts is also tissue-specific and affects functions of multiple biological pathways. Accounting for this mechanism is important for understanding the natural and engineering the synthetic cellular circuits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3958-1) contains supplementary material, which is available to authorized users.

The alternative splicing of EAM8 contributes to early flowering and short-season adaptation in a landrace barley from the Qinghai-Tibetan Plateau

The early flowering of Lalu was determined to be due to a novel spontaneous eam8 mutation, which resulted in intron retention and the formation of a putative truncated protein. Barley is a staple crop grown over an extensive area in the Qinghai-Tibetan Plateau. Understanding the genetic mechanism for its success in a high altitude is important for crop improvement in marginal environments. Early flowering is a critical adaptive trait that strongly influences reproductive fitness in a short growing season. Loss-of-function mutations at the circadian clock gene EARLY MATURITY 8 (EAM8) promote rapid flowering. In this study, we identified a novel, spontaneous recessive eam8 mutant with an early flowering phenotype in a Tibetan barley landrace Lalu, which is natively grown at a high altitude of approximately 4000 m asl. The co-segregation analysis in a F₂ population derived from the cross Lalu (early flowering) × Diqing 1 (late flowering) confirmed that early flowering of Lalu was determined to be due to an allele at EAM8. The eam8 allele from Lalu carries an A/G alternative splicing mutation at position 3257 in intron 3, designated eam8.l; this alternative splicing event leads to intron retention and a putative truncated protein. Of the 134 sequenced barley accessions, which are primarily native to the Qinghai-Tibet Plateau, three accessions carried this mutation. The eam8.l mutation was likely to have originated in wild barley due to the presence of the Lalu haplotype in H. spontaneum from Tibet. Overall, alternative splicing has contributed to the evolution of the barley circadian clock and in the short-season adaptation of local barley germplasm. The study has also identified a novel donor of early-flowering barley which will be useful for barley improvement.

Optimizing Rubisco and its regulation for greater resource use efficiency

Rubisco catalyses the carboxylation of ribulose-1,5-bisphosphate (RuBP), enabling net CO2 assimilation in photosynthesis. The properties and regulation of Rubisco are not optimal for biomass production in current and projected future environments. Rubisco is relatively inefficient, and large amounts of the enzyme are needed to support photosynthesis, requiring large investments in nitrogen. The competing oxygenation of RuBP by Rubisco decreases photosynthetic efficiency. Additionally, Rubisco is inhibited by some sugar phosphates and depends upon interaction with Rubisco activase (Rca) to be reactivated. Rca activity is modulated by the chloroplast redox status and ADP/ATP ratios, thereby mediating Rubisco activation and photosynthetic induction in response to irradiance. The extreme thermal sensitivity of Rca compromises net CO2 assimilation at moderately high temperatures. Given its central role in carbon assimilation, the improvement of Rubisco function and regulation is tightly linked with irradiance, nitrogen and water use efficiencies. Although past attempts have had limited success, novel technologies and an expanding knowledge base make the challenge of improving Rubisco activity in crops an achievable goal. Strategies to optimize Rubisco and its regulation are addressed in relation to their potential to improve crop resource use efficiency and climate resilience of photosynthesis.

Environmental stresses modulate abundance and timing of alternatively spliced circadian transcripts in Arabidopsis

Environmental stresses profoundly altered accumulation of nonsense mRNAs including intron-retaining (IR) transcripts in Arabidopsis. Temporal patterns of stress-induced IR mRNAs were dissected using both oscillating and non-oscillating transcripts. Broad-range thermal cycles triggered a sharp increase in the long IR CCA1 isoforms and altered their phasing to different times of day. Both abiotic and biotic stresses such as drought or Pseudomonas syringae infection induced a similar increase. Thermal stress induced a time delay in accumulation of CCA1 I4Rb transcripts, whereas functional mRNA showed steady oscillations. Our data favor a hypothesis that stress-induced instabilities of the central oscillator can be in part compensated through fluctuations in abundance and out-of-phase oscillations of CCA1 IR transcripts. Taken together, our results support a concept that mRNA abundance can be modulated through altering ratios between functional and nonsense/IR transcripts. SR45 protein specifically bound to the retained CCA1 intron in vitro, suggesting that this splicing factor could be involved in regulation of intron retention. Transcriptomes of nonsense-mediated mRNA decay (NMD)-impaired and heat-stressed plants shared a set of retained introns associated with stress- and defense-inducible transcripts. Constitutive activation of certain stress response networks in an NMD mutant could be linked to disequilibrium between functional and nonsense mRNAs.

Circadian rhythms and post-transcriptional regulation in higher plants

The circadian clock of plants allows them to cope with daily changes in their environment. This is accomplished by the rhythmic regulation of gene expression, in a process that involves many regulatory steps. One of the key steps involved at the RNA level is post-transcriptional regulation, which ensures a correct control on the different amounts and types of mRNA that will ultimately define the current physiological state of the plant cell. Recent advances in the study of the processes of regulation of pre-mRNA processing, RNA turn-over and surveillance, regulation of translation, function of lncRNAs, biogenesis and function of small RNAs, and the development of bioinformatics tools have helped to vastly expand our understanding of how this regulatory step performs its role. In this work we review the current progress in circadian regulation at the post-transcriptional level research in plants. It is the continuous interaction of all the information flow control post-transcriptional processes that allow a plant to precisely time and predict daily environmental changes.

RNA around the clock - regulation at the RNA level in biological timing

The circadian timing system in plants synchronizes their physiological functions with the environment. This is achieved by a global control of gene expression programs with a considerable part of the transcriptome undergoing 24-h oscillations in steady-state abundance. These circadian oscillations are driven by a set of core clock proteins that generate their own 24-h rhythm through periodic feedback on their own transcription. Additionally, post-transcriptional events are instrumental for oscillations of core clock genes and genes in clock output. Here we provide an update on molecular events at the RNA level that contribute to the 24-h rhythm of the core clock proteins and shape the circadian transcriptome. We focus on the circadian system of the model plant Arabidopsis thaliana but also discuss selected regulatory principles in other organisms.

Circadian expression profiles of chromatin remodeling factor genes in Arabidopsis

The circadian clock is a biological time keeper mechanism that regulates biological rhythms to a period of approximately 24 h. The circadian clock enables organisms to anticipate environmental cycles and coordinates internal cellular physiology with external environmental cues. In plants, correct matching of the clock with the environment confers fitness advantages to plant survival and reproduction. Therefore, circadian clock components are regulated at multiple layers to fine-tune the circadian oscillation. Epigenetic regulation provides an additional layer of circadian control. However, little is known about which chromatin remodeling factors are responsible for circadian control. In this work, we analyzed circadian expression of 109 chromatin remodeling factor genes and identified 17 genes that display circadian oscillation. In addition, we also found that a candidate interacts with a core clock component, supporting that clock activity is regulated in part by chromatin modification. As an initial attempt to elucidate the relationship between chromatin modification and circadian oscillation, we identified novel regulatory candidates that provide a platform for future investigations of chromatin regulation of the circadian clock.

A protocol for visual analysis of alternative splicing in RNA-Seq data using integrated genome browser

Ultrahigh-throughput sequencing of cDNA (RNA-Seq) is an invaluable resource for investigating alternative splicing in an organism. Alternative splicing is a form of posttranscriptional regulation in which primary RNA transcripts from a single gene can be spliced in multiple ways leading to different RNA and protein products. In plants and other species, it has been shown that many genes involved in circadian regulation are alternatively spliced. As new RNA-Seq data sets become available, these data will lead to new insights into links between regulation RNA splicing and the circadian system. Analyzing RNA-Seq data sets requires software tools that can display RNA-Seq read alignments alongside gene models, enabling assessment of how treatments or developmental stages affect splicing patterns and production of novel variants. The Integrated Genome Browser (IGB) software program is a free and flexible desktop tool that enables discovery and quantification of alternative splicing. In this protocol, we use IGB and a cold-stress RNA-Seq data set to examine alternative splicing of Arabidopsis thaliana LHY, a circadian clock regulator. IGB is freely available from http://www.bioviz.org .

Role for LSM genes in the regulation of circadian rhythms

Growing evidence suggests that core spliceosomal components differentially affect RNA processing of specific genes; however, whether changes in the levels or activities of these factors control specific signaling pathways is largely unknown. Here we show that some SM-like (LSM) genes, which encode core components of the spliceosomal U6 small nuclear ribonucleoprotein complex, regulate circadian rhythms in plants and mammals. We found that the circadian clock regulates the expression of LSM5 in Arabidopsis plants and several LSM genes in mouse suprachiasmatic nucleus. Further, mutations in LSM5 or LSM4 in Arabidopsis, or down-regulation of LSM3, LSM5, or LSM7 expression in human cells, lengthens the circadian period. Although we identified changes in the expression and alternative splicing of some core clock genes in Arabidopsis lsm5 mutants, the precise molecular mechanism causing period lengthening remains to be identified. Genome-wide expression analysis of either a weak lsm5 or a strong lsm4 mutant allele in Arabidopsis revealed larger effects on alternative splicing than on constitutive splicing. Remarkably, large splicing defects were not observed in most of the introns evaluated using RNA-seq in the strong lsm4 mutant allele used in this study. These findings support the idea that some LSM genes play both regulatory and constitutive roles in RNA processing, contributing to the fine-tuning of specific signaling pathways.

Circadian regulation of gene expression: at the crossroads of transcriptional and post-transcriptional regulatory networks

Gene expression programs activated in response to, or in anticipation of, environmental changes involve sequential steps, from transcription and RNA processing to nuclear export and translation. Here we review recent advances in our understanding of the multiple regulatory layers that control the oscillations in gene expression associated with daily rhythms in metabolism and physiology across eukaryotic organisms. Whereas many genes show coordinated oscillations in transcription, RNA processing and translation, others show significant temporal disconnections between these processes. Thus, circadian oscillations constitute an ideal system for examining how multiple transcriptional and post-transcriptional regulatory steps are integrated to maximize organismal adjustments to environmental conditions.