Catsnap: a user-friendly algorithm for determining the conservation of protein variants reveals extensive parallelisms in the evolution of alternative splicing.

Understanding the evolutionary conservation of complex eukaryotic transcriptomes significantly illuminates the physiological relevance of alternative splicing (AS). Examining the evolutionary depth of a given AS event with ordinary homology searches is generally challenging and time-consuming. Here, we present Catsnap, an algorithmic pipeline for assessing the conservation of putative protein isoforms generated by AS. It employs a machine learning approach following a database search with the provided pair of protein sequences. We used the Catsnap algorithm for analyzing the conservation of emerging experimentally characterized alternative proteins from plants and animals. Indeed, most of them are conserved among other species. Catsnap can detect the conserved functional protein isoforms regardless of the AS type by which they are generated. Notably, we found that while the primary amino acid sequence is maintained, the type of AS determining the inclusion or exclusion of protein regions varies throughout plant phylogenetic lineages in these proteins. We also document that this phenomenon is less seen among animals. In sum, our algorithm highlights the presence of unexpectedly frequent hotspots where protein isoforms recurrently arise to carry physiologically relevant functions. The user web interface is available at https://catsnap.cesnet.cz/.

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.

Xyloglucan Remodeling Defines Auxin-Dependent Differential Tissue Expansion in Plants.

Size control is a fundamental question in biology, showing incremental complexity in plants, whose cells possess a rigid cell wall. The phytohormone auxin is a vital growth regulator with central importance for differential growth control. Our results indicate that auxin-reliant growth programs affect the molecular complexity of xyloglucans, the major type of cell wall hemicellulose in eudicots. Auxin-dependent induction and repression of growth coincide with reduced and enhanced molecular complexity of xyloglucans, respectively. In agreement with a proposed function in growth control, genetic interference with xyloglucan side decorations distinctly modulates auxin-dependent differential growth rates. Our work proposes that auxin-dependent growth programs have a spatially defined effect on xyloglucan's molecular structure, which in turn affects cell wall mechanics and specifies differential, gravitropic hypocotyl growth.

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.

Differential nucleosome occupancy modulates alternative splicing in Arabidopsis thaliana.

Alternative splicing (AS) is a major gene regulatory mechanism in plants. Recent evidence supports co-transcriptional splicing in plants, hence the chromatin state can impact AS. However, how dynamic changes in the chromatin state such as nucleosome occupancy influence the cold-induced AS remains poorly understood. Here, we generated transcriptome (RNA-Seq) and nucleosome positioning (MNase-Seq) data for Arabidopsis thaliana to understand how nucleosome positioning modulates cold-induced AS. Our results show that characteristic nucleosome occupancy levels are strongly associated with the type and abundance of various AS events under normal and cold temperature conditions in Arabidopsis. Intriguingly, exitrons, alternatively spliced internal regions of protein-coding exons, exhibit distinctive nucleosome positioning pattern compared to other alternatively spliced regions. Likewise, nucleosome patterns differ between exitrons and retained introns, pointing to their distinct regulation. Collectively, our data show that characteristic changes in nucleosome positioning modulate AS in plants in response to cold.

FRUITFULL Is a Repressor of Apical Hook Opening in Arabidopsis thaliana.

Plants adjust their architecture to a constantly changing environment, requiring adaptation of differential growth. Despite their importance, molecular switches, which define growth transitions, are largely unknown. Apical hook development in dark grown Arabidopsis thaliana (A. thaliana) seedlings serves as a suitable model for differential growth transition in plants. Here, we show that the phytohormone auxin counteracts the light-induced growth transition during apical hook opening. We, subsequently, identified genes which are inversely regulated by light and auxin. We used in silico analysis of the regulatory elements in this set of genes and subsequently used natural variation in gene expression to uncover correlations between underlying transcription factors and the in silico predicted target genes. This approach uncovered that MADS box transcription factor AGAMOUS-LIKE 8 (AGL8)/FRUITFULL (FUL) modulates apical hook opening. Our data shows that transient FUL expression represses the expression of growth stimulating genes during early phases of apical hook development and therewith guards the transition to growth promotion for apical hook opening. Here, we propose a role for FUL in setting tissue identity, thereby regulating differential growth during apical hook development.

A Collection of Pre-mRNA Splicing Mutants in Arabidopsis thaliana.

To investigate factors influencing pre-mRNA splicing in plants, we conducted a forward genetic screen using an alternatively-spliced GFP reporter gene in Arabidopsis thaliana This effort generated a collection of sixteen mutants impaired in various splicing-related proteins, many of which had not been recovered in any prior genetic screen or implicated in splicing in plants. The factors are predicted to act at different steps of the spliceosomal cycle, snRNP biogenesis pathway, transcription, and mRNA transport. We have described eleven of the mutants in recent publications. Here we present the final five mutants, which are defective, respectively, in RNA-BINDING PROTEIN 45D (rbp45d), DIGEORGE SYNDROME CRITICAL REGION 14 (dgcr14), CYCLIN-DEPENDENT KINASE G2 (cdkg2), INTERACTS WITH SPT6 (iws1) and CAP BINDING PROTEIN 80 (cbp80). We provide RNA-sequencing data and analyses of differential gene expression and alternative splicing patterns for the cbp80 mutant and for several previously published mutants, including smfa and new alleles of cwc16a, for which such information was not yet available. Sequencing of small RNAs from the cbp80 mutant highlighted the necessity of wild-type CBP80 for processing of microRNA (miRNA) precursors into mature miRNAs. Redundancy tests of paralogs encoding several of the splicing factors revealed their functional non-equivalence in the GFP reporter gene system. We discuss the cumulative findings and their implications for the regulation of pre-mRNA splicing efficiency and alternative splicing in plants. The mutant collection provides a unique resource for further studies on a coherent set of splicing factors and their roles in gene expression, alternative splicing and plant development.

Current Challenges in Studying Alternative Splicing in Plants: The Case of Physcomitrella patens SR Proteins.

To colonize different terrestrial habitats, early land plants had to overcome the challenge of coping with harsh new environments. Alternative splicing - an RNA processing mechanism through which splice sites are differentially recognized, originating multiple transcripts and potentially different proteins from the same gene - can be key for plant stress tolerance. Serine/arginine-rich (SR) proteins constitute an evolutionarily conserved family of major alternative splicing regulators that in plants subdivides into six subfamilies. Despite being well studied in animals and a few plant species, such as the model angiosperm Arabidopsis thaliana and the crop Oryza sativa, little is known of these splicing factors in early land plants. Establishing the whole complement of SR proteins in different species is essential to understand the functional and evolutionary significance of alternative splicing. An in silico search for SR proteins in the extant moss Physcomitrella patens revealed inconsistencies both in the published data and available databases, likely arising from automatic annotation lacking adequate manual curation. These misannotations interfere with the description not only of the number and subfamily classification of Physcomitrella SR proteins but also of their domain architecture, potentially hindering the elucidation of their molecular functions. We therefore advise caution when looking into P. patens genomic resources. Our systematic survey nonetheless confidently identified 16 P. patens SR proteins that fall into the six described subfamilies and represent counterparts of well-established members in Arabidopsis and rice. Intensified research efforts should disclose whether SR proteins were already determining alternative splicing modulation and stress tolerance in early land plants.

Alternative Splicing and DNA Damage Response in Plants.

Plants are exposed to a variety of abiotic and biotic stresses that may result in DNA damage. Endogenous processes - such as DNA replication, DNA recombination, respiration, or photosynthesis - are also a threat to DNA integrity. It is therefore essential to understand the strategies plants have developed for DNA damage detection, signaling, and repair. Alternative splicing (AS) is a key post-transcriptional process with a role in regulation of gene expression. Recent studies demonstrate that the majority of intron-containing genes in plants are alternatively spliced, highlighting the importance of AS in plant development and stress response. Not only does AS ensure a versatile proteome and influence the abundance and availability of proteins greatly, it has also emerged as an important player in the DNA damage response (DDR) in animals. Despite extensive studies of DDR carried out in plants, its regulation at the level of AS has not been comprehensively addressed. Here, we provide some insights into the interplay between AS and DDR in plants.

Does co-transcriptional regulation of alternative splicing mediate plant stress responses?

Plants display exquisite control over gene expression to elicit appropriate responses under normal and stress conditions. Alternative splicing (AS) of pre-mRNAs, a process that generates two or more transcripts from multi-exon genes, adds another layer of regulation to fine-tune condition-specific gene expression in animals and plants. However, exactly how plants control splice isoform ratios and the timing of this regulation in response to environmental signals remains elusive. In mammals, recent evidence indicate that epigenetic and epitranscriptome changes, such as DNA methylation, chromatin modifications and RNA methylation, regulate RNA polymerase II processivity, co-transcriptional splicing, and stability and translation efficiency of splice isoforms. In plants, the role of epigenetic modifications in regulating transcription rate and mRNA abundance under stress is beginning to emerge. However, the mechanisms by which epigenetic and epitranscriptomic modifications regulate AS and translation efficiency require further research. Dynamic changes in the chromatin landscape in response to stress may provide a scaffold around which gene expression, AS and translation are orchestrated. Finally, we discuss CRISPR/Cas-based strategies for engineering chromatin architecture to manipulate AS patterns (or splice isoforms levels) to obtain insight into the epigenetic regulation of AS.

PRP4KA, a Putative Spliceosomal Protein Kinase, Is Important for Alternative Splicing and Development in Arabidopsis thaliana.

Splicing of precursor messenger RNAs (pre-mRNAs) is an essential step in the expression of most eukaryotic genes. Both constitutive splicing and alternative splicing, which produces multiple messenger RNA (mRNA) isoforms from a single primary transcript, are modulated by reversible protein phosphorylation. Although the plant splicing machinery is known to be a target for phosphorylation, the protein kinases involved remain to be fully defined. We report here the identification of pre-mRNA processing 4 (PRP4) KINASE A (PRP4KA) in a forward genetic screen based on an alternatively spliced GFP reporter gene in Arabidopsis thaliana (Arabidopsis). Prp4 kinase is the first spliceosome-associated kinase shown to regulate splicing in fungi and mammals but it has not yet been studied in plants. In the same screen we identified mutants defective in SAC3A, a putative mRNA export factor that is highly coexpressed with PRP4KA in Arabidopsis Whereas the sac3a mutants appear normal, the prp4ka mutants display a pleiotropic phenotype featuring atypical rosettes, late flowering, tall final stature, reduced branching, and lowered seed set. Analysis of RNA-sequencing data from prp4ka and sac3a mutants identified widespread and partially overlapping perturbations in alternative splicing in the two mutants. Quantitative phosphoproteomic profiling of a prp4ka mutant detected phosphorylation changes in several serine/arginine-rich proteins, which regulate constitutive and alternative splicing, and other splicing-related factors. Tests of PRP4KB, the paralog of PRP4KA, indicated that the two genes are not functionally redundant. The results demonstrate the importance of PRP4KA for alternative splicing and plant phenotype, and suggest that PRP4KA may influence alternative splicing patterns by phosphorylating a subset of splicing regulators.

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.

Generating Targeted Gene Knockout Lines in Physcomitrella patens to Study Evolution of Stress-Responsive Mechanisms.

The moss Physcomitrella patens possesses highly efficient homologous recombination allowing targeted gene manipulations and displays many features of the early land plants including high tolerance to abiotic stresses. It is therefore an invaluable model organism for studies of gene functions and comparative studies of evolution of stress responses in plants. Here, we describe a method for generating targeted gene knockout lines in P. patens using a polyethylene glycol-mediated transformation of protoplasts including basic in vitro growth, propagation, and maintenance techniques.

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.

Unmasking alternative splicing inside protein-coding exons defines exitrons and their role in proteome plasticity.

Alternative splicing (AS) diversifies transcriptomes and proteomes and is widely recognized as a key mechanism for regulating gene expression. Previously, in an analysis of intron retention events in Arabidopsis, we found unusual AS events inside annotated protein-coding exons. Here, we also identify such AS events in human and use these two sets to analyse their features, regulation, functional impact, and evolutionary origin. As these events involve introns with features of both introns and protein-coding exons, we name them exitrons (exonic introns). Though exitrons were detected as a subset of retained introns, they are clearly distinguishable, and their splicing results in transcripts with different fates. About half of the 1002 Arabidopsis and 923 human exitrons have sizes of multiples of 3 nucleotides (nt). Splicing of these exitrons results in internally deleted proteins and affects protein domains, disordered regions, and various post-translational modification sites, thus broadly impacting protein function. Exitron splicing is regulated across tissues, in response to stress and in carcinogenesis. Intriguingly, annotated intronless genes can be also alternatively spliced via exitron usage. We demonstrate that at least some exitrons originate from ancestral coding exons. Based on our findings, we propose a "splicing memory" hypothesis whereby upon intron loss imprints of former exon borders defined by vestigial splicing regulatory elements could drive the evolution of exitron splicing. Altogether, our studies show that exitron splicing is a conserved strategy for increasing proteome plasticity in plants and animals, complementing the repertoire of AS events.

Shedding light on the chloroplast as a remote control of nuclear gene expression.

Plants rely on a sophisticated light sensing and signaling system that allows them to respond to environmental changes. Photosensory protein systems -phytochromes, cryptochromes, phototropins, and ultraviolet (UV)-B photoreceptors- have evolved to let plants monitor light conditions and regulate different levels of gene expression and developmental processes. However, even though photoreceptor proteins are best characterized and deeply studied, it is also known that chloroplasts are able to sense light conditions and communicate the variations to the nucleus that adjust its transcriptome to the changing environment. The redox state of components of the photosynthetic electron transport chain works as a sensor of photosynthetic activity and can affect nuclear gene expression by a retrograde signaling pathway. Recently, our groups showed that a retrograde signaling pathway can modulate the alternative splicing process, revealing a novel layer of gene expression control by chloroplast retrograde signaling.

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.