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.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

BACKGROUND

Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis.

RESULTS

We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage.

CONCLUSIONS

AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

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.

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.

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.

Alternative Splicing of Circadian Clock Genes Correlates With Temperature in Field-Grown Sugarcane

Alternative Splicing (AS) is a mechanism that generates different mature transcripts from precursor mRNAs (pre-mRNAs) of the same gene. In plants, a wide range of physiological and metabolic events are related to AS, as well as fast responses to changes in temperature. AS is present in around 60% of intron-containing genes in Arabidopsis, 46% in rice, and 38% in maize and it is widespread among the circadian clock genes. Little is known about how AS influences the circadian clock of C4 plants, like commercial sugarcane, a C4 crop with a complex hybrid genome. This work aims to test if the daily dynamics of AS forms of circadian clock genes are regulated by environmental factors, such as temperature, in the field. A systematic search for AS in five sugarcane clock genes, ScLHY, ScPRR37, ScPRR73, ScPRR95, and ScTOC1 using different organs of sugarcane sampled during winter, with 4 months old plants, and during summer, with 9 months old plants, revealed temperature- and organ-dependent expression of at least one alternatively spliced isoform in all genes. Expression of AS isoforms varied according to the season. Our results suggest that AS events in circadian clock genes are correlated with temperature.

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.

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.

Cold-Dependent Expression and Alternative Splicing of Arabidopsis Long Non-coding RNAs

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.

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.

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.

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.

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.

Evaluation and improvement of the regulatory inference for large co-expression networks with limited sample size

BACKGROUND

Co-expression has been widely used to identify novel regulatory relationships using high throughput measurements, such as microarray and RNA-seq data. Evaluation studies on co-expression network analysis methods mostly focus on networks of small or medium size of up to a few hundred nodes. For large networks, simulated expression data usually consist of hundreds or thousands of profiles with different perturbations or knock-outs, which is uncommon in real experiments due to their cost and the amount of work required. Thus, the performances of co-expression network analysis methods on large co-expression networks consisting of a few thousand nodes, with only a small number of profiles with a single perturbation, which more accurately reflect normal experimental conditions, are generally uncharacterized and unknown.

METHODS

We proposed a novel network inference methods based on Relevance Low order Partial Correlation (RLowPC). RLowPC method uses a two-step approach to select on the high-confidence edges first by reducing the search space by only picking the top ranked genes from an intial partial correlation analysis and, then computes the partial correlations in the confined search space by only removing the linear dependencies from the shared neighbours, largely ignoring the genes showing lower association.

RESULTS

We selected six co-expression-based methods with good performance in evaluation studies from the literature: Partial correlation, PCIT, ARACNE, MRNET, MRNETB and CLR. The evaluation of these methods was carried out on simulated time-series data with various network sizes ranging from 100 to 3000 nodes. Simulation results show low precision and recall for all of the above methods for large networks with a small number of expression profiles. We improved the inference significantly by refinement of the top weighted edges in the pre-inferred partial correlation networks using RLowPC. We found improved performance by partitioning large networks into smaller co-expressed modules when assessing the method performance within these modules.

CONCLUSIONS

The evaluation results show that current methods suffer from low precision and recall for large co-expression networks where only a small number of profiles are available. The proposed RLowPC method effectively reduces the indirect edges predicted as regulatory relationships and increases the precision of top ranked predictions. Partitioning large networks into smaller highly co-expressed modules also helps to improve the performance of network inference methods. The RLowPC R package for network construction, refinement and evaluation is available at GitHub: https://github.com/wyguo/RLowPC .

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.

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.

Monitoring Alternative Splicing Changes in Arabidopsis Circadian Clock Genes

Posttranscriptional control makes an important contribution to circadian regulation of gene expression. In higher plants, alternative splicing is particularly prevalent upon abiotic and biotic stress and in the circadian system. Here we describe in detail a high-resolution reverse transcription-PCR based panel (HR RT-PCR) to monitor alternative splicing events. The use of the panel allows the quantification of changes in the proportion of splice isoforms between different samples, e.g., different time points, different tissues, genotypes, ecotypes, or treatments.

AtRTD - a comprehensive reference transcript dataset resource for accurate quantification of transcript-specific expression in Arabidopsis thaliana

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.

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.

Alternative splicing of mini-exons in the Arabidopsis leaf rust receptor-like kinase LRK10 genes affects subcellular localisation

KEY MESSAGE

AtLRK10L1.2 produces a variety of alternatively spliced variants in the region a mini-exon and skipping of the mini-exon alters the subcellular localization of the protein. We have examined expression and alternative splicing in the gene encoding Arabidopsis LRK10-like 1 (AtLRK10L1) which is most closely related to wheat leaf rust 10 disease-resistance locus receptor-like protein kinase (LRK10). AtLRK10L1 produces two different transcripts, LRK10L1.1 and 1.2 through the use of two different promoters. We found no evidence of alternative splicing for the AtLRK10L1.1 transcript but identified numerous alternative splicing variants of AtLRK10L1.2 by sequencing of cloned cDNAs prepared from RNA isolated from whole cell, nucleolar and nucleoplasmic fractions. Many of these transcripts contained unspliced introns and accumulated differentially in the nucleolus and the nucleoplasm consistent with intron retention transcripts being retained in the nucleus (Göhring et al., Plant Cell 26:754-764, 2014). We examined the fate of different alternatively spliced transcripts by fusing variants to YFP and expressing them by agroinfiltration in Nicotiana benthamiana. AtLRK10L1 contains a 45 nt mini-exon which encodes part of a putative transmembrane domain. Full-length cDNA of LRK10L1.2 fused to YFP targeted the fusion protein to the plasma membrane while expression of transcripts where the mini-exon had been deleted, altered the localization of the fusion protein to the endoplasmic reticulum. Similarly, expression of full-length and mini-exon deleted versions of three other members of the LRK10 receptor-like kinase (RLK) gene family also showed the switch in localization. Thus, the mini-exons in Arabidopsis LRK10 genes are required for localization to the plasma membrane.

Long noncoding RNA modulates alternative splicing regulators in Arabidopsis

Alternative splicing (AS) of pre-mRNA represents a major mechanism underlying increased transcriptome and proteome complexity. Here, we show that the nuclear speckle RNA-binding protein (NSR) and the AS competitor long noncoding RNA (or ASCO-lncRNA) constitute an AS regulatory module. AtNSR-GFP translational fusions are expressed in primary and lateral root (LR) meristems. Double Atnsr mutants and ASCO overexpressors exhibit an altered ability to form LRs after auxin treatment. Interestingly, auxin induces a major change in AS patterns of many genes, a response largely dependent on NSRs. RNA immunoprecipitation assays demonstrate that AtNSRs interact not only with their alternatively spliced mRNA targets but also with the ASCO-RNA in vivo. The ASCO-RNA displaces an AS target from an NSR-containing complex in vitro. Expression of ASCO-RNA in Arabidopsis affects the splicing patterns of several NSR-regulated mRNA targets. Hence, lncRNA can hijack nuclear AS regulators to modulate AS patterns during development.

Arabidopsis PTB1 and PTB2 proteins negatively regulate splicing of a mini-exon splicing reporter and affect alternative splicing of endogenous genes differentially

This paper examines the function of Arabidopsis thaliana AtPTB1 and AtPTB2 as plant splicing factors. The effect on splicing of overexpression of AtPTB1 and AtPTB2 was analysed in an in vivo protoplast transient expression system with a novel mini-exon splicing reporter. A range of mutations in pyrimidine-rich sequences were compared with and without AtPTB and NpU2AF65 overexpression. Splicing analyses of constructs in protoplasts and RNA from overexpression lines used high-resolution reverse transcription polymerase chain reaction (RT-PCR). AtPTB1 and AtPTB2 reduced inclusion/splicing of the potato invertase mini-exon splicing reporter, indicating that these proteins can repress plant intron splicing. Mutation of the polypyrimidine tract and closely associated Cytosine and Uracil-rich (CU-rich) sequences, upstream of the mini-exon, altered repression by AtPTB1 and AtPTB2. Coexpression of a plant orthologue of U2AF65 alleviated the splicing repression of AtPTB1. Mutation of a second CU-rich upstream of the mini-exon 3' splice site led to a decline in mini-exon splicing, indicating the presence of a splicing enhancer sequence. Finally, RT-PCR of AtPTB overexpression lines with c. 90 known alternative splicing (AS) events showed that AtPTBs significantly altered AS of over half the events. AtPTB1 and AtPTB2 are splicing factors that influence alternative splicing. This occurs in the potato invertase mini-exon via the polypyrimidine tract and associated pyrimidine-rich sequence.

A chloroplast retrograde signal regulates nuclear alternative splicing

Light is a source of energy and also a regulator of plant physiological adaptations. We show here that light/dark conditions affect alternative splicing of a subset of Arabidopsis genes preferentially encoding proteins involved in RNA processing. The effect requires functional chloroplasts and is also observed in roots when the communication with the photosynthetic tissues is not interrupted, suggesting that a signaling molecule travels through the plant. Using photosynthetic electron transfer inhibitors with different mechanisms of action, we deduce that the reduced pool of plastoquinones initiates a chloroplast retrograde signaling that regulates nuclear alternative splicing and is necessary for proper plant responses to varying light 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.

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.

Plant SILAC: stable-isotope labelling with amino acids of arabidopsis seedlings for quantitative proteomics

Stable Isotope Labelling by Amino acids in Cell culture (SILAC) is a powerful technique for comparative quantitative proteomics, which has recently been applied to a number of different eukaryotic organisms. Inefficient incorporation of labelled amino acids in cell cultures of Arabidopsis thaliana has led to very limited use of SILAC in plant systems. We present a method allowing, for the first time, efficient labelling with stable isotope-containing arginine and lysine of whole Arabidopsis seedlings. To illustrate the utility of this method, we have combined the high labelling efficiency (>95%) with quantitative proteomics analyses of seedlings exposed to increased salt concentration. In plants treated for 7 days with 80 mM NaCl, a relatively mild salt stress, 215 proteins were identified whose expression levels changed significantly compared to untreated seedling controls. The 92 up-regulated proteins included proteins involved in abiotic stress responses and photosynthesis, while the 123 down-regulated proteins were enriched in proteins involved in reduction of oxidative stress and other stress responses, respectively. Efficient labelling of whole Arabidopsis seedlings by this modified SILAC method opens new opportunities to exploit the genetic resources of Arabidopsis and analyse the impact of mutations on quantitative protein dynamics in vivo.

An hnRNP-like RNA-binding protein affects alternative splicing by in vivo interaction with transcripts in Arabidopsis thaliana

Alternative splicing (AS) of pre-mRNAs is an important regulatory mechanism shaping the transcriptome. In plants, only few RNA-binding proteins are known to affect AS. Here, we show that the glycine-rich RNA-binding protein AtGRP7 influences AS in Arabidopsis thaliana. Using a high-resolution RT–PCR-based AS panel, we found significant changes in the ratios of AS isoforms for 59 of 288 analyzed AS events upon ectopic AtGRP7 expression. In particular, AtGRP7 affected the choice of alternative 5′ splice sites preferentially. About half of the events are also influenced by the paralog AtGRP8, indicating that AtGRP7 and AtGRP8 share a network of downstream targets. For 10 events, the AS patterns were altered in opposite directions in plants with elevated AtGRP7 level or lacking AtGRP7. Importantly, RNA immunoprecipitation from plant extracts showed that several transcripts are bound by AtGRP7 in vivo and indeed represent direct targets. Furthermore, the effect of AtGRP7 on these AS events was abrogated by mutation of a single arginine that is required for its RNA-binding activity. This indicates that AtGRP7 impacts AS of these transcripts via direct interaction. As several of the AS events are also controlled by other splicing regulators, our data begin to provide insights into an AS network in Arabidopsis.

Thermoplasticity in the plant circadian clock: how plants tell the time-perature

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.

Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis

Alternative splicing (AS) is a key regulatory mechanism that contributes to transcriptome and proteome diversity. As very few genome-wide studies analyzing AS in plants are available, we have performed high-throughput sequencing of a normalized cDNA library which resulted in a high coverage transcriptome map of Arabidopsis. We detect ∼150,000 splice junctions derived mostly from typical plant introns, including an eightfold increase in the number of U12 introns (2069). Around 61% of multiexonic genes are alternatively spliced under normal growth conditions. Moreover, we provide experimental validation of 540 AS transcripts (from 256 genes coding for important regulatory factors) using high-resolution RT-PCR and Sanger sequencing. Intron retention (IR) is the most frequent AS event (∼40%), but many IRs have relatively low read coverage and are less well-represented in assembled transcripts. Additionally, ∼51% of Arabidopsis genes produce AS transcripts which do not involve IR. Therefore, the significance of IR in generating transcript diversity was generally overestimated in previous assessments. IR analysis allowed the identification of a large set of cryptic introns inside annotated coding exons. Importantly, a significant fraction of these cryptic introns are spliced out in frame, indicating a role in protein diversity. Furthermore, we show extensive AS coupled to nonsense-mediated decay in AFC2, encoding a highly conserved LAMMER kinase which phosphorylates splicing factors, thus establishing a complex loop in AS regulation. We provide the most comprehensive analysis of AS to date which will serve as a valuable resource for the plant community to study transcriptome complexity and gene regulation.

Alternative splicing and nonsense-mediated decay modulate expression of important regulatory genes in Arabidopsis

Alternative splicing (AS) coupled to nonsense-mediated decay (NMD) is a post-transcriptional mechanism for regulating gene expression. We have used a high-resolution AS RT-PCR panel to identify endogenous AS isoforms which increase in abundance when NMD is impaired in the Arabidopsis NMD factor mutants, upf1-5 and upf3-1. Of 270 AS genes (950 transcripts) on the panel, 102 transcripts from 97 genes (32%) were identified as NMD targets. Extrapolating from these data around 13% of intron-containing genes in the Arabidopsis genome are potentially regulated by AS/NMD. This cohort of naturally occurring NMD-sensitive AS transcripts also allowed the analysis of the signals for NMD in plants. We show the importance of AS in introns in 5' or 3'UTRs in modulating NMD-sensitivity of mRNA transcripts. In particular, we identified upstream open reading frames overlapping the main start codon as a new trigger for NMD in plants and determined that NMD is induced if 3'-UTRs were >350 nt. Unexpectedly, although many intron retention transcripts possess NMD features, they are not sensitive to NMD. Finally, we have shown that AS/NMD regulates the abundance of transcripts of many genes important for plant development and adaptation including transcription factors, RNA processing factors and stress response genes.

Involvement of the plant nucleolus in virus and viroid infections: parallels with animal pathosystems

The nucleolus is a dynamic subnuclear body with roles in ribosome subunit biogenesis, mediation of cell-stress responses, and regulation of cell growth. An increasing number of reports reveal that similar to the proteins of animal viruses, many plant virus proteins localize in the nucleolus to divert host nucleolar proteins from their natural functions in order to exert novel role(s) in the virus infection cycle. This chapter will highlight studies showing how plant viruses recruit nucleolar functions to facilitate virus translation and replication, virus movement and assembly of virus-specific ribonucleoprotein (RNP) particles, and to counteract plant host defense responses. Plant viruses also provide a valuable tool to gain new insights into novel nucleolar functions and processes. Investigating the interactions between plant viruses and the nucleolus will facilitate the design of novel strategies to control plant virus infections.

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.

Plant U13 orthologues and orphan snoRNAs identified by RNomics of RNA from Arabidopsis nucleoli

Small nucleolar RNAs (snoRNAs) and small Cajal body-specific RNAs (scaRNAs) are non-coding RNAs whose main function in eukaryotes is to guide the modification of nucleotides in ribosomal and spliceosomal small nuclear RNAs, respectively. Full-length sequences of Arabidopsis snoRNAs and scaRNAs have been obtained from cDNA libraries of capped and uncapped small RNAs using RNA from isolated nucleoli from Arabidopsis cell cultures. We have identified 31 novel snoRNA genes (9 box C/D and 22 box H/ACA) and 15 new variants of previously described snoRNAs. Three related capped snoRNAs with a distinct gene organization and structure were identified as orthologues of animal U13snoRNAs. In addition, eight of the novel genes had no complementarity to rRNAs or snRNAs and are therefore putative orphan snoRNAs potentially reflecting wider functions for these RNAs. The nucleolar localization of a number of the snoRNAs and the localization to nuclear bodies of two putative scaRNAs was confirmed by in situ hybridization. The majority of the novel snoRNA genes were found in new gene clusters or as part of previously described clusters. These results expand the repertoire of Arabidopsis snoRNAs to 188 snoRNA genes with 294 gene variants.

Localization of eIF4A-III in the nucleolus and splicing speckles is an indicator of plant stress

The mechanisms of long-term adaptation to low oxygen environment are quite well studied, but little is known about the sensing of oxygen shortage, the signal transduction and the shortterm effects of hypoxia in plant cells. We have found that an RNA helicase eIF4A-III, a putative component of the Exon Junction Complex, rapidly changes its pattern of localisation in the plant nucleus under hypoxic conditions. In normal cell growth conditions GFPeIF4A-III was mainly nucleoplasmic, but in hypoxia stress conditions it moved to the nucleolus and splicing speckles. This transition occurred within 15-20 min in Arabidopsis culture cells and seedling root cells, but took more than 2 h in tobacco BY-2 culture cells. Inhibition of respiration, transcription or phosphorylation in cells and ethanol treatment had similar effects to hypoxia. The most likely consequence is that a certain mRNA population will remain bound to the eIF4A-III and other mRNA processing proteins, rather than being transported from the nucleus to the cytoplasm, and thus its translation will be suspended.

Filtering of deep sequencing data reveals the existence of abundant Dicer-dependent small RNAs derived from tRNAs

Deep sequencing technologies such as Illumina, SOLiD, and 454 platforms have become very powerful tools in discovering and quantifying small RNAs in diverse organisms. Sequencing small RNA fractions always identifies RNAs derived from abundant RNA species such as rRNAs, tRNAs, snRNA, and snoRNA, and they are widely considered to be random degradation products. We carried out bioinformatic analysis of deep sequenced HeLa RNA and after quality filtering, identified highly abundant small RNA fragments, derived from mature tRNAs that are likely produced by specific processing rather than from random degradation. Moreover, we showed that the processing of small RNAs derived from tRNA(Gln) is dependent on Dicer in vivo and that Dicer cleaves the tRNA in vitro.

Involvement of the nuclear cap-binding protein complex in alternative splicing in Arabidopsis thaliana

The nuclear cap-binding protein complex (CBC) participates in 5′ splice site selection of introns that are proximal to the mRNA cap. However, it is not known whether CBC has a role in alternative splicing. Using an RT–PCR alternative splicing panel, we analysed 435 alternative splicing events in Arabidopsis thaliana genes, encoding mainly transcription factors, splicing factors and stress-related proteins. Splicing profiles were determined in wild type plants, the cbp20 and cbp80(abh1) single mutants and the cbp20/80 double mutant. The alternative splicing events included alternative 5′ and 3′ splice site selection, exon skipping and intron retention. Significant changes in the ratios of alternative splicing isoforms were found in 101 genes. Of these, 41% were common to all three CBC mutants and 15% were observed only in the double mutant. The cbp80(abh1) and cbp20/80 mutants had many more changes in alternative splicing in common than did cbp20 and cbp20/80 suggesting that CBP80 plays a more significant role in alternative splicing than CBP20, probably being a platform for interactions with other splicing factors. Cap-binding proteins and the CBC are therefore directly involved in alternative splicing of some Arabidopsis genes and in most cases influenced alternative splicing of the first intron, particularly at the 5′ splice site.

Aberrant mRNA transcripts and the nonsense-mediated decay proteins UPF2 and UPF3 are enriched in the Arabidopsis nucleolus

The eukaryotic nucleolus is multifunctional and involved in the metabolism and assembly of many different RNAs and ribonucleoprotein particles as well as in cellular functions, such as cell division and transcriptional silencing in plants. We previously showed that Arabidopsis thaliana exon junction complex proteins associate with the nucleolus, suggesting a role for the nucleolus in mRNA production. Here, we report that the plant nucleolus contains mRNAs, including fully spliced, aberrantly spliced, and single exon gene transcripts. Aberrant mRNAs are much more abundant in nucleolar fractions, while fully spliced products are more abundant in nucleoplasmic fractions. The majority of the aberrant transcripts contain premature termination codons and have characteristics of nonsense-mediated decay (NMD) substrates. A direct link between NMD and the nucleolus is shown by increased levels of the same aberrant transcripts in both the nucleolus and in Up-frameshift (upf) mutants impaired in NMD. In addition, the NMD factors UPF3 and UPF2 localize to the nucleolus, suggesting that the Arabidopsis nucleolus is therefore involved in identifying aberrant mRNAs and NMD.

Dynamic behavior of Arabidopsis eIF4A-III, putative core protein of exon junction complex: fast relocation to nucleolus and splicing speckles under hypoxia

Here, we identify the Arabidopsis thaliana ortholog of the mammalian DEAD box helicase, eIF4A-III, the putative anchor protein of exon junction complex (EJC) on mRNA. Arabidopsis eIF4A-III interacts with an ortholog of the core EJC component, ALY/Ref, and colocalizes with other EJC components, such as Mago, Y14, and RNPS1, suggesting a similar function in EJC assembly to animal eIF4A-III. A green fluorescent protein (GFP)-eIF4A-III fusion protein showed localization to several subnuclear domains: to the nucleoplasm during normal growth and to the nucleolus and splicing speckles in response to hypoxia. Treatment with the respiratory inhibitor sodium azide produced an identical response to the hypoxia stress. Treatment with the proteasome inhibitor MG132 led to accumulation of GFP-eIF4A-III mainly in the nucleolus, suggesting that transition of eIF4A-III between subnuclear domains and/or accumulation in nuclear speckles is controlled by proteolysis-labile factors. As revealed by fluorescence recovery after photobleaching analysis, the nucleoplasmic fraction was highly mobile, while the speckles were the least mobile fractions, and the nucleolar fraction had an intermediate mobility. Sequestration of eIF4A-III into nuclear pools with different mobility is likely to reflect the transcriptional and mRNA processing state of the cell.

Intronic noncoding RNAs and splicing

The gene organization of small nucleolar RNAs (snoRNAs) and microRNAs (miRNAs) varies within and among different organisms. This diversity is reflected in the maturation pathways of these small noncoding RNAs (ncRNAs). The presence of noncoding RNAs in introns has implications for the biogenesis of both mature small RNAs and host mRNA. The balance of the interactions between the processing or ribonucleoprotein assembly of intronic noncoding RNAs and the splicing process can regulate the levels of ncRNA and host mRNA. The processing of snoRNAs - both intronic and non-intronic - is well characterised in yeast, plants and animals and provides a basis for examining how intronic plant miRNAs are processed.

Monitoring changes in alternative precursor messenger RNA splicing in multiple gene transcripts

Alternative splicing (AS) increases the proteomic and functional capacity of genomes through the generation of alternative mRNA transcripts from the same gene. AS is now estimated to occur in a third of Arabidopsis and rice genes, and includes genes involved in the control of growth and development, responses to stress and signalling. Regulation of AS reflects the interactions between positive and negative cis sequences in the precursor messenger RNA and a range of trans-acting factors. The levels and activities of these factors differ in different cells and growth conditions. To identify changes in AS in multiple genes simultaneously, we have established a reproducible RT-PCR panel that can analyse 96 alternative splicing events and accurately measure the ratio of alternatively spliced products. This procedure detected statistically significant changes in AS in different plant organs, in plants grown under different light and day-length conditions, and in plants overexpressing splicing factors. The system provides a convenient, medium-throughput means of monitoring changes in AS in multiple genes. It can readily be applied to much larger or targeted sets of gene transcripts to generate information on the significance and regulation of AS in plant growth and development, specific processes and responses to external stimuli.

A plant virus movement protein forms ringlike complexes with the major nucleolar protein, fibrillarin, in vitro

Fibrillarin, one of the major proteins of the nucleolus, has methyltransferase activity directing 2′-O-ribose methylation of rRNA and snRNAs and is required for rRNA processing. The ability of the plant umbravirus, groundnut rosette virus, to move long distances through the phloem, the specialized plant vascular system, has been shown to strictly depend on the interaction of one of its proteins, the ORF3 protein (protein encoded by open reading frame 3), with fibrillarin. This interaction is essential for several stages in the groundnut rosette virus life cycle such as nucleolar import of the ORF3 protein via Cajal bodies, relocalization of some fibrillarin from the nucleolus to cytoplasm, and assembly of cytoplasmic umbraviral ribonucleoprotein particles that are themselves required for the long-distance spread of the virus and systemic infection. Here, using atomic force microscopy, we determine the architecture of these complexes as single-layered ringlike structures with a diameter of 18–22 nm and a height of 2.0 ± 0.4 nm, which consist of several (n = 6–8) distinct protein granules. We also estimate the molar ratio of fibrillarin to ORF3 protein in the complexes as approximately 1:1. Based on these data, we propose a model of the structural organization of fibrillarin–ORF3 protein complexes and discuss potential mechanistic and functional implications that may also apply to other viruses.

Interaction of a plant virus-encoded protein with the major nucleolar protein fibrillarin is required for systemic virus infection

The nucleolus and specific nucleolar proteins are involved in the life cycles of some plant and animal viruses, but the functions of these proteins and of nucleolar trafficking in virus infections are largely unknown. The ORF3 protein of the plant virus, groundnut rosette virus (an umbravirus), has been shown to cycle through the nucleus, passing through Cajal bodies to the nucleolus and then exiting back into the cytoplasm. This journey is absolutely required for the formation of viral ribonucleoprotein particles (RNPs) that, themselves, are essential for the spread of the virus to noninoculated leaves of the shoot tip. Here, we show that these processes rely on the interaction of the ORF3 protein with fibrillarin, a major nucleolar protein. Silencing of the fibrillarin gene prevents long-distance movement of groundnut rosette virus but does not affect viral replication or cell-to-cell movement. Repressing fibrillarin production also localizes the ORF3 protein to multiple Cajal body-like aggregates that fail to fuse with the nucleolus. Umbraviral ORF3 protein and fibrillarin interact in vitro and, when mixed with umbravirus RNA, form an RNP complex. This complex has a filamentous structure with some regular helical features, resembling the RNP complex formed in vivo during umbravirus infection. The filaments formed in vitro are infectious when inoculated to plants, and their infectivity is resistant to RNase. These results demonstrate previously undescribed functions for fibrillarin as an essential component of translocatable viral RNPs and may have implications for other plant and animal viruses that interact with the nucleolus.

Cajal bodies and the nucleolus are required for a plant virus systemic infection

The nucleolus and Cajal bodies (CBs) are prominent interacting subnuclear domains involved in a number of crucial aspects of cell function. Certain viruses interact with these compartments but the functions of such interactions are largely uncharacterized. Here, we show that the ability of the groundnut rosette virus open reading frame (ORF) 3 protein to move viral RNA long distances through the phloem strictly depends on its interaction with CBs and the nucleolus. The ORF3 protein targets and reorganizes CBs into multiple CB-like structures and then enters the nucleolus by causing fusion of these structures with the nucleolus. The nucleolar localization of the ORF3 protein is essential for subsequent formation of viral ribonucleoprotein (RNP) particles capable of virus long-distance movement and systemic infection. We provide a model whereby the ORF3 protein utilizes trafficking pathways involving CBs to enter the nucleolus and, along with fibrillarin, exit the nucleus to form viral 'transport-competent' RNP particles in the cytoplasm.

Arabidopsis nucleolar protein database (AtNoPDB)

The Arabidopsis Nucleolar Protein Database (http://bioinf.scri.sari.ac.uk/cgi-bin/atnopdb/home) provides information on 217 proteins identified in a proteomic analysis of nucleoli isolated from Arabidopsis cell culture. The database is organized on the basis of the Arabidopsis gene identifier number. The information provided includes protein description, protein class, whether or not the plant protein has a homologue in the most recent human nucleolar proteome and the results of reciprocal BLAST analysis of the human proteome. In addition, for one-third of the 217 Arabidopsis nucleolar proteins, localization images are available from analysis of full-length cDNA–green fluorescent protein (GFP) fusions and the strength of signal in different parts of the cell—nucleolus, nucleolus-associated structures, nucleoplasm, nuclear bodies and extra-nuclear—is provided. For each protein, the most likely human and yeast orthologues, where identifiable through BLASTX analysis, are given with links to relevant information sources.

Proteomic analysis of the Arabidopsis nucleolus suggests novel nucleolar functions

The eukaryotic nucleolus is involved in ribosome biogenesis and a wide range of other RNA metabolism and cellular functions. An important step in the functional analysis of the nucleolus is to determine the complement of proteins of this nuclear compartment. Here, we describe the first proteomic analysis of plant (Arabidopsis thaliana) nucleoli, in which we have identified 217 proteins. This allows a direct comparison of the proteomes of an important nuclear structure between two widely divergent species: human and Arabidopsis. The comparison identified many common proteins, plant-specific proteins, proteins of unknown function found in both proteomes, and proteins that were nucleolar in plants but nonnucleolar in human. Seventy-two proteins were expressed as GFP fusions and 87% showed nucleolar or nucleolar-associated localization. In a striking and unexpected finding, we have identified six components of the postsplicing exon-junction complex (EJC) involved in mRNA export and nonsense-mediated decay (NMD)/mRNA surveillance. This association was confirmed by GFP-fusion protein localization. These results raise the possibility that in plants, nucleoli may have additional functions in mRNA export or surveillance.

Plant nuclear bodies

Knowledge of the organization of transcription, RNA processing and transport, and the assembly of complexes such as the ribosome, spliceosome and other RNPs is essential to understanding gene expression. Over several years, the nucleolus and Cajal bodies have been examined in plants, and recently, various other sub-nuclear domains that are involved in RNA metabolism and hormonal responses have been discovered. These novel domains illustrate the complexity and subtlety of expression control and herald a new era of research on the molecular and cell biology of plant nuclei.

Involvement of the nucleolus in plant virus systemic infection

The nucleolus is a prominent subnuclear domain and is classically regarded as the site of transcription of rRNA, processing of the precursor rRNAs and biogenesis of pre-ribosomal particles. In addition to these traditionally recognized activities, the nucleolus also participates in many other aspects of cell function. The umbravirus-encoded ORF3 protein is a multifunctional RNA-binding protein involved in long-distance RNA movement, and protection of viral RNA from RNase attack, including possibly small interfering RNA-guided RNA silencing. In addition to its presence in cytoplasmic ribonucleoprotein particles containing viral RNA, the umbraviral ORF3 protein accumulates in nuclei, preferentially targeting nucleoli. The ORF3 protein domains involved in the localization of the protein to the nucleolus were identified. Functional analysis of the mutants revealed the correlation between the ORF3 protein nucleolar localization and its ability to form the cytoplasmic ribonucleoprotein particles and transport viral RNA long distances via the phloem. Possible mechanisms of the nucleolar involvement in systemic virus infection are discussed.