REVEILLE2 thermosensitive splicing: a molecular basis for the integration of nocturnal temperature information by the Arabidopsis circadian clock.

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.

BaRTv2: a highly resolved barley reference transcriptome for accurate transcript-specific RNA-seq quantification.

Accurate characterisation of splice junctions (SJs) as well as transcription start and end sites in reference transcriptomes allows precise quantification of transcripts from RNA-seq data, and enables detailed investigations of transcriptional and post-transcriptional regulation. Using novel computational methods and a combination of PacBio Iso-seq and Illumina short-read sequences from 20 diverse tissues and conditions, we generated a comprehensive and highly resolved barley reference transcript dataset from the European 2-row spring barley cultivar Barke (BaRTv2.18). Stringent and thorough filtering was carried out to maintain the quality and accuracy of the SJs and transcript start and end sites. BaRTv2.18 shows increased transcript diversity and completeness compared with an earlier version, BaRTv1.0. The accuracy of transcript level quantification, SJs and transcript start and end sites have been validated extensively using parallel technologies and analysis, including high-resolution reverse transcriptase-polymerase chain reaction and 5'-RACE. BaRTv2.18 contains 39 434 genes and 148 260 transcripts, representing the most comprehensive and resolved reference transcriptome in barley to date. It provides an important and high-quality resource for advanced transcriptomic analyses, including both transcriptional and post-transcriptional regulation, with exceptional resolution and precision.

Nonsense-Mediated RNA Decay Factor UPF1 Is Critical for Posttranscriptional and Translational Gene Regulation in Arabidopsis.

Nonsense-mediated RNA decay (NMD) is an RNA control mechanism that has also been implicated in the broader regulation of gene expression. Nevertheless, a role for NMD in genome regulation has not yet been fully assessed, partially because NMD inactivation is lethal in many organisms. Here, we performed an in-depth comparative analysis of Arabidopsis (Arabidopsis thaliana) mutants lacking the NMD-related proteins UPF3, UPF1, and SMG7. We found different impacts of these proteins on NMD and the Arabidopsis transcriptome, with UPF1 having the biggest effect. Transcriptome assembly in UPF1-null plants revealed genome-wide changes in alternative splicing, suggesting that UPF1 functions in splicing. The inactivation of UPF1 led to translational repression, as manifested by a global shift in mRNAs from polysomes to monosomes and the downregulation of genes involved in translation and ribosome biogenesis. Despite these global changes, NMD targets and mRNAs expressed at low levels with short half-lives were enriched in the polysomes of upf1 mutants, indicating that UPF1/NMD suppresses the translation of aberrant RNAs. Particularly striking was an increase in the translation of TIR domain-containing, nucleotide binding, leucine-rich repeat (TNL) immune receptors. The regulation of TNLs via UPF1/NMD-mediated mRNA stability and translational derepression offers a dynamic mechanism for the rapid activation of TNLs in response to pathogen attack.

Rapid and Dynamic Alternative Splicing Impacts the Arabidopsis Cold Response Transcriptome.

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.

3D RNA-seq: a powerful and flexible tool for rapid and accurate differential expression and alternative splicing analysis of RNA-seq data for biologists.

RNA-sequencing (RNA-seq) analysis of gene expression and alternative splicing should be routine and robust but is often a bottleneck for biologists because of different and complex analysis programs and reliance on specialized bioinformatics skills. We have developed the '3D RNA-seq' App, an R shiny App and web-based pipeline for the comprehensive analysis of RNA-seq data from any organism. It represents an easy-to-use, flexible and powerful tool for analysis of both gene and transcript-level gene expression to identify differential gene/transcript expression, differential alternative splicing and differential transcript usage (3D) as well as isoform switching from RNA-seq data. 3D RNA-seq integrates state-of-the-art differential expression analysis tools and adopts best practice for RNA-seq analysis. The program is designed to be run by biologists with minimal bioinformatics experience (or by bioinformaticians) allowing lab scientists to analyse their RNA-seq data. It achieves this by operating through a user-friendly graphical interface which automates the data flow through the programs in the pipeline. The comprehensive analysis performed by 3D RNA-seq is extremely rapid and accurate, can handle complex experimental designs, allows user setting of statistical parameters, visualizes the results through graphics and tables, and generates publication quality figures such as heat-maps, expression profiles and GO enrichment plots. The utility of 3D RNA-seq is illustrated by analysis of data from a time-series of cold-treated Arabidopsis plants and from dexamethasone-treated male and female mouse cortex and hypothalamus data identifying dexamethasone-induced sex- and brain region-specific differential gene expression and alternative splicing.

Illuminating the dark side of the human transcriptome with long read transcript sequencing.

BACKGROUND: The human transcriptome annotation is regarded as one of the most complete of any eukaryotic species. However, limitations in sequencing technologies have biased the annotation toward multi-exonic protein coding genes. Accurate high-throughput long read transcript sequencing can now provide additional evidence for rare transcripts and genes such as mono-exonic and non-coding genes that were previously either undetectable or impossible to differentiate from sequencing noise. RESULTS: We developed the Transcriptome Annotation by Modular Algorithms (TAMA) software to leverage the power of long read transcript sequencing and address the issues with current data processing pipelines. TAMA achieved high sensitivity and precision for gene and transcript model predictions in both reference guided and unguided approaches in our benchmark tests using simulated Pacific Biosciences (PacBio) and Nanopore sequencing data and real PacBio datasets. By analyzing PacBio Sequel II Iso-Seq sequencing data of the Universal Human Reference RNA (UHRR) using TAMA and other commonly used tools, we found that the convention of using alignment identity to measure error correction performance does not reflect actual gain in accuracy of predicted transcript models. In addition, inter-read error correction can cause major changes to read mapping, resulting in potentially over 6 K erroneous gene model predictions in the Iso-Seq based human genome annotation. Using TAMA's genome assembly based error correction and gene feature evidence, we predicted 2566 putative novel non-coding genes and 1557 putative novel protein coding gene models. CONCLUSIONS: Long read transcript sequencing data has the power to identify novel genes within the highly annotated human genome. The use of parameter tuning and extensive output information of the TAMA software package allows for in depth exploration of eukaryotic transcriptomes. We have found long read data based evidence for thousands of unannotated genes within the human genome. More development in sequencing library preparation and data processing are required for differentiating sequencing noise from real genes in long read RNA sequencing data.

BaRTv1.0: an improved barley reference transcript dataset to determine accurate changes in the barley transcriptome using RNA-seq.

BACKGROUND: The time required to analyse RNA-seq data varies considerably, due to discrete steps for computational assembly, quantification of gene expression and splicing analysis. Recent fast non-alignment tools such as Kallisto and Salmon overcome these problems, but these tools require a high quality, comprehensive reference transcripts dataset (RTD), which are rarely available in plants. RESULTS: A high-quality, non-redundant barley gene RTD and database (Barley Reference Transcripts - BaRTv1.0) has been generated. BaRTv1.0, was constructed from a range of tissues, cultivars and abiotic treatments and transcripts assembled and aligned to the barley cv. Morex reference genome (Mascher et al. Nature; 544: 427-433, 2017). Full-length cDNAs from the barley variety Haruna nijo (Matsumoto et al. Plant Physiol; 156: 20-28, 2011) determined transcript coverage, and high-resolution RT-PCR validated alternatively spliced (AS) transcripts of 86 genes in five different organs and tissue. These methods were used as benchmarks to select an optimal barley RTD. BaRTv1.0-Quantification of Alternatively Spliced Isoforms (QUASI) was also made to overcome inaccurate quantification due to variation in 5' and 3' UTR ends of transcripts. BaRTv1.0-QUASI was used for accurate transcript quantification of RNA-seq data of five barley organs/tissues. This analysis identified 20,972 significant differentially expressed genes, 2791 differentially alternatively spliced genes and 2768 transcripts with differential transcript usage. CONCLUSION: A high confidence barley reference transcript dataset consisting of 60,444 genes with 177,240 transcripts has been generated. Compared to current barley transcripts, BaRTv1.0 transcripts are generally longer, have less fragmentation and improved gene models that are well supported by splice junction reads. Precise transcript quantification using BaRTv1.0 allows routine analysis of gene expression and AS.

The Arabidopsis SR45 Splicing Factor, a Negative Regulator of Sugar Signaling, Modulates SNF1-Related Protein Kinase 1 Stability.

The ability to sense and respond to sugar signals allows plants to cope with environmental and metabolic changes by adjusting growth and development accordingly. We previously reported that the SR45 splicing factor negatively regulates glucose signaling during early seedling development in Arabidopsis thaliana Here, we show that under glucose-fed conditions, the Arabidopsis sr45-1 loss-of-function mutant contains higher amounts of the energy-sensing SNF1-Related Protein Kinase 1 (SnRK1) despite unaffected SnRK1 transcript levels. In agreement, marker genes for SnRK1 activity are upregulated in sr45-1 plants, and the glucose hypersensitivity of sr45-1 is attenuated by disruption of the SnRK1 gene. Using a high-resolution RT-PCR panel, we found that the sr45-1 mutation broadly targets alternative splicing in vivo, including that of the SR45 pre-mRNA itself. Importantly, the enhanced SnRK1 levels in sr45-1 are suppressed by a proteasome inhibitor, indicating that SR45 promotes targeting of the SnRK1 protein for proteasomal destruction. Finally, we demonstrate that SR45 regulates alternative splicing of the Arabidopsis 5PTase13 gene, which encodes an inositol polyphosphate 5-phosphatase previously shown to interact with and regulate the stability of SnRK1 in vitro, thus providing a mechanistic link between SR45 function and the modulation of degradation of the SnRK1 energy sensor in response to sugars.

A high quality Arabidopsis transcriptome for accurate transcript-level analysis of alternative splicing.

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.

TSIS: an R package to infer alternative splicing isoform switches for time-series data.

SUMMARY: An alternative splicing isoform switch is where a pair of transcript isoforms reverse their relative expression abundances in response to external or internal stimuli. Although computational methods are available to study differential alternative splicing, few tools for detection of isoform switches exist and these are based on pairwise comparisons. Here, we provide the TSIS R package, which is the first tool for detecting significant transcript isoform switches in time-series data. The main steps of TSIS are to search for the isoform switch points in the time-series, characterize the switches and filter the results with user input parameters. All the functions are integrated into a Shiny App for ease of implementation of the analysis. AVAILABILITY AND IMPLEMENTATION: The TSIS package is available on GitHub: https://github.com/wyguo/TSIS. CONTACT: runxuan.zhang@hutton.ac.uk.

How does temperature affect splicing events? Isoform switching of splicing factors regulates splicing of LATE ELONGATED HYPOCOTYL (LHY).

One 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.

Lost in Translation: Pitfalls in Deciphering Plant Alternative Splicing Transcripts.

Transcript annotation in plant databases is incomplete and often inaccurate, leading to misinterpretation. As more and more RNA-seq data are generated, plant scientists need to be aware of potential pitfalls and understand the nature and impact of specific alternative splicing transcripts on protein production. A primary area of concern and the topic of this article is the (mis)annotation of open reading frames and premature termination codons. The basic message is that to adequately address expression and functions of transcript isoforms, it is necessary to be able to predict their fate in terms of whether protein isoforms are generated or specific transcripts are unproductive or degraded.

A physical, genetic and functional sequence assembly of the barley genome.

Barley (Hordeum vulgare L.) is among the world's earliest domesticated and most important crop plants. It is diploid with a large haploid genome of 5.1 gigabases (Gb). Here we present an integrated and ordered physical, genetic and functional sequence resource that describes the barley gene-space in a structured whole-genome context. We developed a physical map of 4.98 Gb, with more than 3.90 Gb anchored to a high-resolution genetic map. Projecting a deep whole-genome shotgun assembly, complementary DNA and deep RNA sequence data onto this framework supports 79,379 transcript clusters, including 26,159 'high-confidence' genes with homology support from other plant genomes. Abundant alternative splicing, premature termination codons and novel transcriptionally active regions suggest that post-transcriptional processing forms an important regulatory layer. Survey sequences from diverse accessions reveal a landscape of extensive single-nucleotide variation. Our data provide a platform for both genome-assisted research and enabling contemporary crop improvement.

The spliceosome assembly factor GEMIN2 attenuates the effects of temperature on alternative splicing and circadian rhythms.

The mechanisms by which poikilothermic organisms ensure that biological processes are robust to temperature changes are largely unknown. Temperature compensation, the ability of circadian rhythms to maintain a relatively constant period over the broad range of temperatures resulting from seasonal fluctuations in environmental conditions, is a defining property of circadian networks. Temperature affects the alternative splicing (AS) of several clock genes in fungi, plants, and flies, but the splicing factors that modulate these effects to ensure clock accuracy throughout the year remain to be identified. Here we show that GEMIN2, a spliceosomal small nuclear ribonucleoprotein assembly factor conserved from yeast to humans, modulates low temperature effects on a large subset of pre-mRNA splicing events. In particular, GEMIN2 controls the AS of several clock genes and attenuates the effects of temperature on the circadian period in Arabidopsis thaliana. We conclude that GEMIN2 is a key component of a posttranscriptional regulatory mechanism that ensures the appropriate acclimation of plants to daily and seasonal changes in temperature conditions.

High-quality reference transcript datasets hold the key to transcript-specific RNA-sequencing analysis in plants.

525 I. 525 II. 526 III. 527 IV. 527 V. 529 VI. 529 529 References 529 SUMMARY: Re-programming of the transcriptome involves both transcription and alternative splicing (AS). Some genes are regulated only at the AS level with no change in expression at the gene level. AS data must be incorporated as an essential aspect of the regulation of gene expression. RNA-sequencing (RNA-seq) can deliver both transcriptional and AS information, but accurate methods to analyse the added complexity in RNA-seq data are needed. The construction of a comprehensive reference transcript dataset (RTD) for a specific plant species, variety or accession, from all available sequence data, will immediately allow more robust analysis of RNA-seq data. RTDs will continually evolve and improve, a process that will be more efficient if resources across a community are shared and pooled.

Alternative splicing at the intersection of biological timing, development, and stress responses.

High-throughput sequencing for transcript profiling in plants has revealed that alternative splicing (AS) affects a much higher proportion of the transcriptome than was previously assumed. AS is involved in most plant processes and is particularly prevalent in plants exposed to environmental stress. The identification of mutations in predicted splicing factors and spliceosomal proteins that affect cell fate, the circadian clock, plant defense, and tolerance/sensitivity to abiotic stress all point to a fundamental role of splicing/AS in plant growth, development, and responses to external cues. Splicing factors affect the AS of multiple downstream target genes, thereby transferring signals to alter gene expression via splicing factor/AS networks. The last two to three years have seen an ever-increasing number of examples of functional AS. At a time when the identification of AS in individual genes and at a global level is exploding, this review aims to bring together such examples to illustrate the extent and importance of AS, which are not always obvious from individual publications. It also aims to ensure that plant scientists are aware that AS is likely to occur in the genes that they study and that dynamic changes in AS and its consequences need to be considered routinely.

A methyl transferase links the circadian clock to the regulation of alternative splicing.

Circadian rhythms allow organisms to time biological processes to the most appropriate phases of the day-night cycle. Post-transcriptional regulation is emerging as an important component of circadian networks, but the molecular mechanisms linking the circadian clock to the control of RNA processing are largely unknown. Here we show that PROTEIN ARGININE METHYL TRANSFERASE 5 (PRMT5), which transfers methyl groups to arginine residues present in histones and Sm spliceosomal proteins, links the circadian clock to the control of alternative splicing in plants. Mutations in PRMT5 impair several circadian rhythms in Arabidopsis thaliana and this phenotype is caused, at least in part, by a strong alteration in alternative splicing of the core-clock gene PSEUDO RESPONSE REGULATOR 9 (PRR9). Furthermore, genome-wide studies show that PRMT5 contributes to the regulation of many pre-messenger-RNA splicing events, probably by modulating 5'-splice-site recognition. PRMT5 expression shows daily and circadian oscillations, and this contributes to the mediation of the circadian regulation of expression and alternative splicing of a subset of genes. Circadian rhythms in locomotor activity are also disrupted in dart5-1, a mutant affected in the Drosophila melanogaster PRMT5 homologue, and this is associated with alterations in splicing of the core-clock gene period and several clock-associated genes. Our results demonstrate a key role for PRMT5 in the regulation of alternative splicing and indicate that the interplay between the circadian clock and the regulation of alternative splicing by PRMT5 constitutes a common mechanism that helps organisms to synchronize physiological processes with daily changes in environmental conditions.

The SERRATE protein is involved in alternative splicing in Arabidopsis thaliana.

How alternative splicing (AS) is regulated in plants has not yet been elucidated. Previously, we have shown that the nuclear cap-binding protein complex (AtCBC) is involved in AS in Arabidopsis thaliana. Here we show that both subunits of AtCBC (AtCBP20 and AtCBP80) interact with SERRATE (AtSE), a protein involved in the microRNA biogenesis pathway. Moreover, using a high-resolution reverse transcriptase-polymerase chain reaction AS system we have found that AtSE influences AS in a similar way to the cap-binding complex (CBC), preferentially affecting selection of 5' splice site of first introns. The AtSE protein acts in cooperation with AtCBC: many changes observed in the mutant lacking the correct SERRATE activity were common to those observed in the cbp mutants. Interestingly, significant changes in AS of some genes were also observed in other mutants of plant microRNA biogenesis pathway, hyl1-2 and dcl1-7, but a majority of them did not correspond to the changes observed in the se-1 mutant. Thus, the role of SERRATE in AS regulation is distinct from that of HYL1 and DCL1, and is similar to the regulation of AS in which CBC is involved.

Evolutionary relationships among barley and Arabidopsis core circadian clock and clock-associated genes.

The circadian clock regulates a multitude of plant developmental and metabolic processes. In crop species, it contributes significantly to plant performance and productivity and to the adaptation and geographical range over which crops can be grown. To understand the clock in barley and how it relates to the components in the Arabidopsis thaliana clock, we have performed a systematic analysis of core circadian clock and clock-associated genes in barley, Arabidopsis and another eight species including tomato, potato, a range of monocotyledonous species and the moss, Physcomitrella patens. We have identified orthologues and paralogues of Arabidopsis genes which are conserved in all species, monocot/dicot differences, species-specific differences and variation in gene copy number (e.g. gene duplications among the various species). We propose that the common ancestor of barley and Arabidopsis had two-thirds of the key clock components identified in Arabidopsis prior to the separation of the monocot/dicot groups. After this separation, multiple independent gene duplication events took place in both monocot and dicot ancestors.