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.

Light regulates alternative splicing outcomes via the TOR kinase pathway.

For plants, light is the source of energy and the most relevant regulator of growth and adaptations to the environment by inducing changes in gene expression at various levels, including alternative splicing. Light-triggered chloroplast retrograde signals control alternative splicing in Arabidopsis thaliana. Here, we provide evidence that light regulates the expression of a core set of splicing-related factors in roots. Alternative splicing responses in roots are not directly caused by light but are instead most likely triggered by photosynthesized sugars. The target of rapamycin (TOR) kinase plays a key role in this shoot-to-root signaling pathway. Knocking down TOR expression or pharmacologically inhibiting TOR activity disrupts the alternative splicing responses to light and exogenous sugars in roots. Consistently, splicing decisions are modulated by mitochondrial activity in roots. In conclusion, by activating the TOR pathway, sugars act as mobile signals to coordinate alternative splicing responses to light throughout the whole plant.

Targeting alternative splicing by RNAi: from the differential impact on splice variants to triggering artificial pre-mRNA splicing.

Alternative splicing generates multiple transcript and protein isoforms from a single gene and controls transcript intracellular localization and stability by coupling to mRNA export and nonsense-mediated mRNA decay (NMD). RNA interference (RNAi) is a potent mechanism to modulate gene expression. However, its interactions with alternative splicing are poorly understood. We used artificial microRNAs (amiRNAs, also termed shRNAmiR) to knockdown all splice variants of selected target genes in Arabidopsis thaliana. We found that splice variants, which vary by their protein-coding capacity, subcellular localization and sensitivity to NMD, are affected differentially by an amiRNA, although all of them contain the target site. Particular transcript isoforms escape amiRNA-mediated degradation due to their nuclear localization. The nuclear and NMD-sensitive isoforms mask RNAi action in alternatively spliced genes. Interestingly, Arabidopsis SPL genes, which undergo alternative splicing and are targets of miR156, are regulated in the same manner. Moreover, similar results were obtained in mammalian cells using siRNAs, indicating cross-kingdom conservation of these interactions among RNAi and splicing isoforms. Furthermore, we report that amiRNA can trigger artificial alternative splicing, thus expanding the RNAi functional repertoire. Our findings unveil novel interactions between different post-transcriptional processes in defining transcript fates and regulating gene expression.

The cyclin-dependent kinase G group defines a thermo-sensitive alternative splicing circuit modulating the expression of Arabidopsis ATU2AF65A.

The ability to adapt growth and development to temperature variations is crucial to generate plant varieties resilient to predicted temperature changes. However, the mechanisms underlying plant response to progressive increases in temperature have just started to be elucidated. Here, we report that the cyclin-dependent kinase G1 (CDKG1) is a central element in a thermo-sensitive mRNA splicing cascade that transduces changes in ambient temperature into differential expression of the fundamental spliceosome component, ATU2AF65A. CDKG1 is alternatively spliced in a temperature-dependent manner. We found that this process is partly dependent on both the cyclin-dependent kinase G2 (CDKG2) and the interacting co-factor CYCLIN L1 (CYCL1), resulting in two distinct messenger RNAs. The relative abundance of both CDKG1 transcripts correlates with ambient temperature and possibly with different expression levels of the associated protein isoforms. Both CDKG1 alternative transcripts are necessary to fully complement the expression of ATU2AF65A across the temperature range. Our data support a previously unidentified temperature-dependent mechanism based on the alternative splicing (AS) of CDKG1 and regulated by CDKG2 and CYCL1. We propose that changes in ambient temperature affect the relative abundance of CDKG1 transcripts, and this in turn translates into differential CDKG1 protein expression coordinating the AS of ATU2AF65A.

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.

The spliceosome-associated protein Nrl1 suppresses homologous recombination-dependent R-loop formation in fission yeast.

The formation of RNA-DNA hybrids, referred to as R-loops, can promote genome instability and cancer development. Yet the mechanisms by which R-loops compromise genome instability are poorly understood. Here, we establish roles for the evolutionarily conserved Nrl1 protein in pre-mRNA splicing regulation, R-loop suppression and in maintaining genome stability. nrl1Δ mutants exhibit endogenous DNA damage, are sensitive to exogenous DNA damage, and have defects in homologous recombination (HR) repair. Concomitantly, nrl1Δ cells display significant changes in gene expression, similar to those induced by DNA damage in wild-type cells. Further, we find that nrl1Δ cells accumulate high levels of R-loops, which co-localize with HR repair factors and require Rad51 and Rad52 for their formation. Together, our findings support a model in which R-loop accumulation and subsequent DNA damage sequesters HR factors, thereby compromising HR repair at endogenously or exogenously induced DNA damage sites, leading to genome instability.

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.

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.