Survey of Conserved Alternative Splicing Events of mRNAs Encoding SR Proteins in Land Plants

Pan-transcriptomic analysis reveals alternative splicing control of cold tolerance in rice

Plants have evolved complex mechanisms to adapt to harsh environmental conditions. Rice (Oryza sativa) is a staple food crop that is sensitive to low temperatures. However, its cold stress responses remain poorly understood, thus limiting possibilities for crop engineering to achieve greater cold tolerance. In this study, we constructed a rice pan-transcriptome and characterized its transcriptional regulatory landscape in response to cold stress. We performed Iso-Seq and RNA-Seq of 11 rice cultivars subjected to a time-course cold treatment. Our analyses revealed that alternative splicing-regulated gene expression plays a significant role in the cold stress response. Moreover, we identified CATALASE C (OsCATC) and Os03g0701200 as candidate genes for engineering enhanced cold tolerance. Importantly, we uncovered central roles for the 2 serine-arginine-rich proteins OsRS33 and OsRS2Z38 in cold tolerance. Our analysis of cold tolerance and resequencing data from a diverse collection of 165 rice cultivars suggested that OsRS2Z38 may be a key selection gene in japonica domestication for cold adaptation, associated with the adaptive evolution of rice. This study systematically investigated the distribution, dynamic changes, and regulatory mechanisms of alternative splicing in rice under cold stress. Overall, our work generates a rich resource with broad implications for understanding the genetic basis of cold response mechanisms in plants.

Transcription and splicing variations of SR genes accompany with genome-wide accumulation of long-introns in pine

Most of mRNAs in Eukaryote were matured after the removal of introns in their pre-mRNA transcripts. Serine/arginine-rich (SR) proteins are a group of splicing regulators regulating the splicing processes globally. Expressions of SR proteins themselves were extensively regulated, at both transcription and splicing levels, under different environmental conditions, specially heat stress conditions. The pine genome is characterized by super-long and easily methylated introns in a large number of genes that derived from the extensive accumulation of transposons (TEs). Here, we identified and analyzed the phylogenetic characteristics of 24 SR proteins and their encoding genes from the pine genome. Then we explored transcription and pre-mRNA splicing expression patterns of SR genes in P. massoniana seedlings under normal and heat stress temperature conditions. Our results showed that the transcription patterns of SR genes in pine exhibited significant changes compared to other plant species, and these changes were not strictly correlated with the intron length and DNA methylation intensity of the SR genes. Interestingly, none of the long introns of SR genes underwent alternative splicing (AS) in our experiment. Furthermore, the intensity of AS regulation may be related to the potential DNA methylation intensity of SR genes. Taken together, this study explores for the first time the characteristics of significant variations in the transcription and splicing patterns of SR proteins in a plant species with an over-accumulation of super-long introns.

OsSCL30 overexpression reduces the tolerance of rice seedlings to low temperature, drought and salt

Rice is one of the main food crops for the world population. Various abiotic stresses, such as low temperature, drought, and high salinity, affect rice during the entire growth period, determining its yield and quality, and even leading to plant death. In this study, by constructing overexpression vectors D-163 + 1300:OsSCL30 and D-163 + 1300-AcGFP:OsSCL30-GFP, the mechanism of action of OsSCL30 in various abiotic stresses was explored. Bioinformatics analysis showed that OsSCL30 was located on the chromosome 12 of rice Nipponbare, belonging to the plant-specific SCL subfamily of the SR protein family. The 1500 bp section upstream of the open reading frame start site contains stress-related cis-acting elements such as ABRE, MYC, and MYB. Under normal conditions, the expression of OsSCL30 was higher in leaves and leaf sheaths. The results of reverse transcription polymerase chain reaction showed that the expression of OsSCL30 decreased after low temperature, drought and salt treatment. In root cells OsSCL30 was localized in the nuclei. The results of the rice seedling tolerance and recovery tests showed that overexpression of OsSCL30 diminished the resistance to low temperature, drought and salt stresses in transgenic rice and resulted in larger accumulation of reactive oxygen species. This study is of great significance for exploring the response mechanisms of SR proteins under abiotic stresses.

Emerging Function of Ecotype-Specific Splicing in the Recruitment of Commensal Microbiome

In recent years, host-microbiome interactions in both animals and plants has emerged as a novel research area for studying the relationship between host organisms and their commensal microbial communities. The fitness advantages of this mutualistic interaction can be found in both plant hosts and their associated microbiome, however, the driving forces mediating this beneficial interaction are poorly understood. Alternative splicing (AS), a pivotal post-transcriptional mechanism, has been demonstrated to play a crucial role in plant development and stress responses among diverse plant ecotypes. This natural variation of plants also has an impact on their commensal microbiome. In this article, we review the current progress of plant natural variation on their microbiome community, and discuss knowledge gaps between AS regulation of plants in response to their intimately related microbiota. Through the impact of this article, an avenue could be established to study the biological mechanism of naturally varied splicing isoforms on plant-associated microbiome assembly.

Overexpression of cassava RSZ21b enhances drought tolerance in Arabidopsis

Drought is one of the major environmental constraints affecting crop productivity. Plants have to adjust their developmental and physiological processes to cope with drought. We previously identified 18 cassava serine/arginine-rich (SR) proteins that had a pivotal role in alternative splicing in response to environmental stress. However, functional characterization of SR proteins is rarely explored. Here, we characterized the RSZ subfamily gene MeRSZ21b in cassava. The RSZ21b belongs to the RSZ subfamily, which was widely distributed in major crops and was highly conserved. Quantitative RT-PCR assay showed that the expression of MeRSZ21b was significantly induced by drought. Moreover, overexpression of MeRSZ21b in Arabidopsis was hypersensitive to abscisic acid (ABA) in the phases of seed germination and post-germination seedling growth. Meantime, MeRSZ21b overexpression lines were resistant to sorbitol treatment, and quickly closed the stomata when compared with Col-0 under drought condition. Importantly, overexpression of MeRSZ21b resulted in improved drought tolerance through modulating ABA-dependent signaling. Therefore, our findings refine our knowledge of the SR protein-coding genes and provide novel insights for enhancing plant resistance to environmental stress.

The SR Splicing Factors: Providing Perspectives on Their Evolution, Expression, Alternative Splicing, and Function in Populus trichocarpa

Serine/arginine-rich (SR) proteins are important splicing factors in plant development and abiotic/hormone-related stresses. However, evidence that SR proteins contribute to the process in woody plants has been lacking. Using phylogenetics, gene synteny, transgenic experiments, and RNA-seq analysis, we identified 24 PtSR genes and explored their evolution, expression, and function in Popolus trichocarpa. The PtSR genes were divided into six subfamilies, generated by at least two events of genome triplication and duplication. Notably, they were constitutively expressed in roots, stems, and leaves, demonstrating their fundamental role in P. trichocarpa. Additionally, most PtSR genes (~83%) responded to at least one stress (cold, drought, salt, SA, MeJA, or ABA), and, especially, cold stress induced a dramatic perturbation in the expression and/or alternative splicing (AS) of 18 PtSR genes (~75%). Evidentially, the overexpression of PtSCL30 in Arabidopsis decreased freezing tolerance, which probably resulted from AS changes of the genes (e.g., ICE2 and COR15A) critical for cold tolerance. Moreover, the transgenic plants were salt-hypersensitive at the germination stage. These indicate that PtSCL30 may act as a negative regulator under cold and salt stress. Altogether, this study sheds light on the evolution, expression, and AS of PtSR genes, and the functional mechanisms of PtSCL30 in woody plants.

Comparative genomics of the coconut crab and other decapod crustaceans: exploring the molecular basis of terrestrial adaptation

Background

The complex life cycle of the coconut crab, Birgus latro, begins when an obligate terrestrial adult female visits the intertidal to hatch zoea larvae into the surf. After drifting for several weeks in the ocean, the post-larval glaucothoes settle in the shallow subtidal zone, undergo metamorphosis, and the early juveniles then subsequently make their way to land where they undergo further physiological changes that prevent them from ever entering the sea again. Here, we sequenced, assembled and analyzed the coconut crab genome to shed light on its adaptation to terrestrial life. For comparison, we also assembled the genomes of the long-tailed marine-living ornate spiny lobster, Panulirus ornatus, and the short-tailed marine-living red king crab, Paralithodes camtschaticus. Our selection of the latter two organisms furthermore allowed us to explore parallel evolution of the crab-like form in anomurans.

Results

All three assembled genomes are large, repeat-rich and AT-rich. Functional analysis reveals that the coconut crab has undergone proliferation of genes involved in the visual, respiratory, olfactory and cytoskeletal systems. Given that the coconut crab has atypical mitochondrial DNA compared to other anomurans, we argue that an abundance of kif22 and other significantly proliferated genes annotated with mitochondrial and microtubule functions, point to unique mechanisms involved in providing cellular energy via nuclear protein-coding genes supplementing mitochondrial and microtubule function. We furthermore detected in the coconut crab a significantly proliferated HOX gene, caudal, that has been associated with posterior development in Drosophila, but we could not definitively associate this gene with carcinization in the Anomura since it is also significantly proliferated in the ornate spiny lobster. However, a cuticle-associated coatomer gene, gammacop, that is significantly proliferated in the coconut crab, may play a role in hardening of the adult coconut crab abdomen in order to mitigate desiccation in terrestrial environments.

Conclusion

The abundance of genomic features in the three assembled genomes serve as a source of hypotheses for future studies of anomuran environmental adaptations such as shell-utilization, perception of visual and olfactory cues in terrestrial environments, and cuticle sclerotization. We hypothesize that the coconut crab exhibits gene proliferation in lieu of alternative splicing as a terrestrial adaptation mechanism and propose life-stage transcriptomic assays to test this hypothesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07636-9.

Differential Transcription and Alternative Splicing in Cotton Underly Specialized Defense Responses Against Pests

The green mirid bug (Apolygus lucorum) and the cotton bollworm (Helicoverpa armigera) are both preferred to live on cotton but cause different symptoms, suggesting specialized responses of cotton to the two insects. In this study, we investigated differential molecular mechanisms underlying cotton plant defenses against A. lucorum and H. armigera via transcriptomic analyses. At the transcription level, jasmonate (JA) signaling was dominated in defense against H. armigera whereas salicylic acid (SA) signaling was more significant in defense against A. lucorum. A set of pathogenesis-related (PR) genes and protease inhibitor genes were differentially induced by the two insects. Insect infestations also had an impact on alternative splicing (AS), which was altered more significantly by the H. armigera than A. lucorum. Interestingly, most differential AS (DAS) genes had no obvious change at the transcription level. GO analysis revealed that biological process termed "RNA splicing" and "cellular response to abiotic stimulus" were enriched only in DAS genes from the H. armigera infested samples. Furthermore, insect infestations induced the retained intron of GhJAZs transcripts, which produced a truncated protein lacking the intact Jas motif. Taken together, our data demonstrate that the specialized cotton response to different insects is regulated by gene transcription and AS as well.

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.

Genomic survey of RNA recognition motif (RRM) containing RNA binding proteins from barley (Hordeum vulgare ssp. vulgare)

One of the major mechanisms of post-transcriptional gene regulation is achieved by proteins bearing well-defined sequence motifs involved in 'RNA binding'. In eukaryotes, RNA binding proteins (RBPs) are key players of RNA metabolism that includes synthesis, processing, editing, modifying, transport, storage and stability of RNA. In plants, the family of RBPs is vastly expanded compared to other eukaryotes including humans. In this study we identified 363 RBPs in the barley genome. Gene ontology enrichment analysis of barley RBPs indicated these proteins were in all the major cellular compartments and associated with key biological processes including translation, splicing, seed development and stress signaling. Members with the classical RNA binding motifs such as the RNA recognition motif (RRM), KH domain, Helicase, CRM, dsRNA and Pumilio were identified in the repertoire of barley RBPs. Similar to Arabidopsis, the RRM containing RBPs were the most abundant in barley genome. In-depth analysis of the RRM containing proteins - polyA binding proteins, Ser/Arg rich proteins and Glycine-rich RBPs were undertaken. Reanalysis of the proteome dataset of various stages during barley malting identified 38 RBPs suggesting an important role for these proteins during the malting process. This survey provides a systematic analysis of barley RBPs and serves as the basis for the further functional characterization of this important family of proteins.

Genome-Wide Identification of Splicing Quantitative Trait Loci (sQTLs) in Diverse Ecotypes of Arabidopsis thaliana

Alternative splicing (AS) of pre-mRNAs contributes to transcriptome diversity and enables plants to generate different protein isoforms from a single gene and/or fine-tune gene expression during different development stages and environmental changes. Although AS is pervasive, the genetic basis for differential isoform usage in plants is still emerging. In this study, we performed genome-wide analysis in 666 geographically distributed diverse ecotypes of Arabidopsis thaliana to identify genomic regions [splicing quantitative trait loci (sQTLs)] that may regulate differential AS. These ecotypes belong to different microclimatic conditions and are part of the relict and non-relict populations. Although sQTLs were spread across the genome, we observed enrichment for trans-sQTL (trans-sQTLs hotspots) on chromosome one. Furthermore, we identified several sQTL (911) that co-localized with trait-linked single nucleotide polymorphisms (SNP) identified in the Arabidopsis genome-wide association studies (AraGWAS). Many sQTLs were enriched among circadian clock, flowering, and stress-responsive genes, suggesting a role for differential isoform usage in regulating these important processes in diverse ecotypes of Arabidopsis. In conclusion, the current study provides a deep insight into SNPs affecting isoform ratios/genes and facilitates a better mechanistic understanding of trait-associated SNPs in GWAS studies. To the best of our knowledge, this is the first report of sQTL analysis in a large set of Arabidopsis ecotypes and can be used as a reference to perform sQTL analysis in the Brassicaceae family. Since whole genome and transcriptome datasets are available for these diverse ecotypes, it could serve as a powerful resource for the biological interpretation of trait-associated loci, splice isoform ratios, and their phenotypic consequences to help produce more resilient and high yield crop varieties.

Evolution of Alternative Splicing in Eudicots

Alternative pre-mRNA splicing (AS) is prevalent in plants and is involved in many interactions between plants and environmental stresses. However, the patterns and underlying mechanisms of AS evolution in plants remain unclear. By analyzing the transcriptomes of four eudicot species, we revealed that the divergence of AS is largely due to the gains and losses of AS events among orthologous genes. Furthermore, based on a subset of AS, in which AS can be directly associated with specific transcripts, we found that AS that generates transcripts containing premature termination codons (PTC), are likely more conserved than those that generate non-PTC containing transcripts. This suggests that AS coupled with nonsense-mediated decay (NMD) might play an important role in affecting mRNA levels post-transcriptionally. To understand the mechanisms underlying the divergence of AS, we analyzed the key determinants of AS using a machine learning approach. We found that the presence/absence of alternative splice site (SS) within the junction, the distance between the authentic SS and the nearest alternative SS, the size of exon-exon junctions were the major determinants for both alternative 5' donor site and 3' acceptor site among the studied species, suggesting a relatively conserved AS mechanism. The comparative analysis further demonstrated that variations of the identified AS determinants significantly contributed to the AS divergence among closely related species in both Solanaceae and Brassicaceae taxa. Together, these results provide detailed insights into the evolution of AS in plants.

Serine/Arginine-rich protein family of splicing regulators: New approaches to study splice isoform functions

Serine/arginine-rich (SR) proteins are conserved RNA-binding proteins that play major roles in RNA metabolism. They function as molecular adaptors, facilitate spliceosome assembly and modulate constitutive and alternative splicing of pre-mRNAs. Pre-mRNAs encoding SR proteins and many other proteins involved in stress responses are extensively alternatively spliced in response to diverse stresses. Hence, it is proposed that stress-induced changes in splice isoforms contribute to the adaptation of plants to stress responses. However, functions of most SR genes and their splice isoforms in stress responses are not known. Lack of easy and robust tools hindered the progress in this area. Emerging technologies such as CRISPR/Cas9 will facilitate studies of SR function by enabling the generation of single and multiple knock-out mutants of SR subfamily members. Moreover, CRISPR/Cas13 allows targeted manipulation of splice isoforms from SR and other genes in a constitutive or tissue-specific manner to evaluate functions of individual splice variants. Identification of the in vivo targets of SR proteins and their splice variants using the recently developed TRIBE (Targets of RNA-binding proteins Identified By Editing) and other methods will help unravel their mode of action and splicing regulatory elements under various conditions. These new approaches are expected to provide significant new insights into the roles of SRs and splice isoforms in plants adaptation to diverse stresses.

Transcriptome response of roots to salt stress in a salinity-tolerant bread wheat cultivar

Salt stress is one of the major adverse environmental factors limiting crop productivity. Considering Iran as one of the bread wheat origins, we sequenced root transcriptome of an Iranian salt tolerant cultivar, Arg, under salt stress to extend our knowledge of the molecular basis of salinity tolerance in Triticum aestivum. RNA sequencing resulted in more than 113 million reads and about 104013 genes were obtained, among which 26171 novel transcripts were identified. A comparison of abundances showed that 5128 genes were differentially expressed due to salt stress. The differentially expressed genes (DEGs) were annotated with Gene Ontology terms, and the key pathways were identified using Kyoto Encyclopedia of Gene and Genomes (KEGG) pathway mapping. The DEGs could be classified into 227 KEGG pathways among which transporters, phenylpropanoid biosynthesis, transcription factors, glycosyltransferases, glutathione metabolism and plant hormone signal transduction represented the most significant pathways. Furthermore, the expression pattern of nine genes involved in salt stress response was compared between the salt tolerant (Arg) and susceptible (Moghan3) cultivars. A panel of novel genes and transcripts is found in this research to be differentially expressed under salinity in Arg cultivar and a model is proposed for salt stress response in this salt tolerant cultivar of wheat employing the DEGs. The achieved results can be beneficial for better understanding and improvement of salt tolerance in wheat.

The chromatin remodeler ZmCHB101 impacts alternative splicing contexts in response to osmotic stress

Maize SWI3-type chromatin remodeler impacts alternative splicing contexts in response to osmotic stress by altering nucleosome density and affecting transcriptional elongation rate. Alternative splicing (AS) is commonly found in higher eukaryotes and is an important posttranscriptional regulatory mechanism to generate transcript diversity. AS has been widely accepted as playing essential roles in different biological processes including growth, development, signal transduction and responses to biotic and abiotic stresses in plants. However, whether and how chromatin remodeling complex functions in AS in plant under osmotic stress remains unknown. Here, we show that a maize SWI3D protein, ZmCHB101, impacts AS contexts in response to osmotic stress. Genome-wide analysis of mRNA contexts in response to osmotic stress using ZmCHB101-RNAi lines reveals that ZmCHB101 impacts alternative splicing contexts of a subset of osmotic stress-responsive genes. Intriguingly, ZmCHB101-mediated regulation of gene expression and AS is largely uncoupled, pointing to diverse molecular functions of ZmCHB101 in transcriptional and posttranscriptional regulation. We further found ZmCHB101 impacts the alternative splicing contexts by influencing alteration of chromatin and histone modification status as well as transcriptional elongation rates mediated by RNA polymerase II. Taken together, our findings suggest a novel insight of how plant chromatin remodeling complex impacts AS under osmotic stress .

Global dissection of alternative splicing uncovers transcriptional diversity in tissues and associates with the flavonoid pathway in tea plant (Camellia sinensis)

BACKGROUND

Alternative splicing (AS) regulates mRNA at the post-transcriptional level to change gene function in organisms. However, little is known about the AS and its roles in tea plant (Camellia sinensis), widely cultivated for making a popular beverage tea.

RESULTS

In our study, the AS landscape and dynamics were characterized in eight tissues (bud, young leaf, summer mature leaf, winter old leaf, stem, root, flower, fruit) of tea plant by Illumina RNA-Seq and confirmed by Iso-Seq. The most abundant AS (~ 20%) was intron retention and involved in RNA processes. The some alternative splicings were found to be tissue specific in stem and root etc. Thirteen co-expressed modules of AS transcripts were identified, which revealed a similar pattern between the bud and young leaves as well as a distinct pattern between seasons. AS events of structural genes including anthocyanidin reductase and MYB transcription factors were involved in biosynthesis of flavonoid, especially in vegetative tissues. The AS isoforms rather than the full-length ones were the major transcripts involved in flavonoid synthesis pathway, and is positively correlated with the catechins content conferring the tea taste. We propose that the AS is an important functional mechanism in regulating flavonoid metabolites.

CONCLUSION

Our study provides the insight into the AS events underlying tea plant's uniquely different developmental process and highlights the important contribution and efficacy of alternative splicing regulatory function to biosynthesis of flavonoids.

SINC-seq: correlation of transient gene expressions between nucleus and cytoplasm reflects single-cell physiology

We report a microfluidic system that physically separates nuclear RNA (nucRNA) and cytoplasmic RNA (cytRNA) from a single cell and enables single-cell integrated nucRNA and cytRNA-sequencing (SINC-seq). SINC-seq constructs two individual RNA-seq libraries, nucRNA and cytRNA, per cell, quantifies gene expression in the subcellular compartments, and combines them to create novel single-cell RNA-seq data. Leveraging SINC-seq, we discover distinct natures of correlation among cytRNA and nucRNA that reflect the transient physiological state of single cells. These data provide unique insights into the regulatory network of messenger RNA from the nucleus toward the cytoplasm at the single-cell level.

Alternative Splicing of Barley Clock Genes in Response to Low Temperature

Alternative splicing (AS) is a regulated mechanism that generates multiple transcripts from individual genes. It is widespread in eukaryotic genomes and provides an effective way to control gene expression. At low temperatures, AS regulates Arabidopsis clock genes through dynamic changes in the levels of productive mRNAs. We examined AS in barley clock genes to assess whether temperature-dependent AS responses also occur in a monocotyledonous crop species. We identify changes in AS of various barley core clock genes including the barley orthologues of Arabidopsis AtLHY and AtPRR7 which showed the most pronounced AS changes in response to low temperature. The AS events modulate the levels of functional and translatable mRNAs, and potentially protein levels, upon transition to cold. There is some conservation of AS events and/or splicing behaviour of clock genes between Arabidopsis and barley. In addition, novel temperature-dependent AS of the core clock gene HvPPD-H1 (a major determinant of photoperiod response and AtPRR7 orthologue) is conserved in monocots. HvPPD-H1 showed a rapid, temperature-sensitive isoform switch which resulted in changes in abundance of AS variants encoding different protein isoforms. This novel layer of low temperature control of clock gene expression, observed in two very different species, will help our understanding of plant adaptation to different environments and ultimately offer a new range of targets for plant improvement.

Functional analysis of alternative splicing of the FLOWERING LOCUS T orthologous gene in Chrysanthemum morifolium

As the junction of floral development pathways, the FLOWERING LOCUS T (FT) protein called 'florigen' plays an important role in the process of plant flowering through signal integration. We isolated four transcripts encoding different isoforms of a FT orthologous gene CmFTL1, from Chrysanthemum morifolium cultivar 'Jimba'. Sequence alignments suggested that the four transcripts are related to the intron 1. Expression analysis showed that four alternative splicing (AS) forms of CmFTL1 varied depending on the developmental stage of the flower. The functional complement experiment using an Arabidopsis mutant ft-10 revealed that the archetypal and AS forms of CmFTL1 had the function of complementing late flower phenotype in different levels. In addition, transgenic confirmation at transcript level showed CmFTL1 and CmFTL1ast coexist in the same tissue type at the same developmental stage, indicating a post-transcriptional modification of CmFTL1 in Arabidopsis. Moreover, ectopic expression of different AS forms in chrysanthemum resulted in the development of multiple altered phenotypes, varying degrees of early flowering. We found that an alternative splicing form (CmFTL1-astE134) without the exon 2 lacked the ability causing the earlier flower phenotype. The evidence in this study indicates that complex alternative processing of CmFTL1 transcripts in C. morifolium may be associated with flowering regulation and hold some potential for biotechnical engineering to create early-flowering phenotypes in ornamental cultivars.

Pre-mRNA Splicing in Plants: In Vivo Functions of RNA-Binding Proteins Implicated in the Splicing Process

Alternative pre-messenger RNA splicing in higher plants emerges as an important layer of regulation upon exposure to exogenous and endogenous cues. Accordingly, mutants defective in RNA-binding proteins predicted to function in the splicing process show severe phenotypic alterations. Among those are developmental defects, impaired responses to pathogen threat or abiotic stress factors, and misregulation of the circadian timing system. A suite of splicing factors has been identified in the model plant Arabidopsis thaliana. Here we summarize recent insights on how defects in these splicing factors impair plant performance.

nagnag: Identification and quantification of NAGNAG alternative splicing using RNA-Seq data

Regulation of proteome diversity by alternative splicing has been widely demonstrated in plants and animals. NAGNAG splicing, which was recently defined as a tissue specific event, results in the production of two distinct isoforms that are distinguished by three nucleotides (NAG) as a consequence of the intron proximal or distal to the splice site. Since the NAGNAG mechanism is not well characterized, tools for the identification and quantification of NAGNAG splicing events remain under-developed. Here we report nagnag, an R-based NAGNAG splicing detection tool, which accurately identifies and quantifies NAGNAG splicing events using RNA-Seq. Overall, nagnag produces user-friendly visualization reports and highlights differences between the DNA/RNA/protein across the identified isoforms of the reported gene. The package is available on https://sourceforge.net/projects/nagnag/files/; or http://genome.sdau.edu.cn/research/software/nagnag.html.

Polyribosomal RNA-Seq reveals the decreased complexity and diversity of the Arabidopsis translatome

Recent RNA-seq studies reveal that the transcriptomes in animals and plants are more complex than previously thought, leading to the inclusion of many more splice isoforms in annotated genomes. However, it is possible that a significant proportion of the transcripts are spurious isoforms that do not contribute to functional proteins. One of the current hypotheses is that commonly used mRNA extraction methods isolate both pre-mature (nuclear) mRNA and mature (cytoplasmic) mRNA, and these incompletely spliced pre-mature mRNAs may contribute to a large proportion of these spurious transcripts. To investigate this, we compared a traditional RNA-seq dataset (total RNA-seq) and a ribosome-bound RNA-seq dataset (polyribosomal RNA-seq) from Arabidopsis thaliana. An integrative framework that combined de novo assembly and genome-guided assembly was applied to reconstruct transcriptomes for the two datasets. Up to 44.8% of the de novo assembled transcripts in total RNA-seq sample were of low abundance, whereas only 0.09% in polyribosomal RNA-seq de novo assembly were of low abundance. The final round of assembly using PASA (Program to Assemble Spliced Alignments) resulted in more transcript assemblies in the total RNA-seq than those in polyribosomal sample. Comparison of alternative splicing (AS) patterns between total and polyribosomal RNA-seq showed a significant difference (G-test, p-value<0.01) in intron retention events: 46.4% of AS events in the total sample were intron retention, whereas only 23.5% showed evidence of intron retention in the polyribosomal sample. It is likely that a large proportion of retained introns in total RNA-seq result from incompletely spliced pre-mature mRNA. Overall, this study demonstrated that polyribosomal RNA-seq technology decreased the complexity and diversity of the coding transcriptome by eliminating pre-mature mRNAs, especially those of low abundance.

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.

Global dissection of alternative splicing in paleopolyploid soybean

Alternative splicing (AS) is common in higher eukaryotes and plays an important role in gene posttranscriptional regulation. It has been suggested that AS varies dramatically among species, tissues, and duplicated gene families of different sizes. However, the genomic forces that govern AS variation remain poorly understood. Here, through genome-wide identification of AS events in the soybean (Glycine max) genome using high-throughput RNA sequencing of 28 samples from different developmental stages, we found that more than 63% of multiexonic genes underwent AS. More AS events occurred in the younger developmental stages than in the older developmental stages for the same type of tissue, and the four main AS types, exon skipping, intron retention, alternative donor sites, and alternative acceptor sites, exhibited different characteristics. Global computational analysis demonstrated that the variations of AS frequency and AS types were significantly correlated with the changes of gene features and gene transcriptional level. Further investigation suggested that the decrease of AS within the genome-wide duplicated genes were due to the diminution of intron length, exon number, and transcriptional level. Altogether, our study revealed that a large number of genes were alternatively spliced in the soybean genome and that variations in gene structure and transcriptional level may play important roles in regulating AS.

Discovery and expression analysis of alternative splicing events conserved among plant SR proteins

The high frequency of alternative splicing among the serine/arginine-rich (SR) family of proteins in plants has been linked to important roles in gene regulation during development and in response to environmental stress. In this article, we have searched and manually annotated all the SR proteins in the genomes of maize and sorghum. The experimental validation of gene structure by reverse transcription-polymerase chain reaction (RT-PCR) analysis revealed, with few exceptions, that SR genes produced multiple isoforms of transcripts by alternative splicing. Despite sharing high structural similarity and conserved positions of the introns, the profile of alternative splicing diverged significantly between maize and sorghum for the vast majority of SR genes. These include many transcript isoforms discovered by RT-PCR and not represented in extant expressed sequence tag (EST) collection. However, we report the occurrence of various maize and sorghum SR mRNA isoforms that display evolutionary conservation of splicing events with their homologous SR genes in Arabidopsis and moss. Our data also indicate an important role of both 5' and 3' untranslated regions in the regulation of SR gene expression. These observations have potentially important implications for the processes of evolution and adaptation of plants to land.

Complexity of the alternative splicing landscape in plants

Alternative splicing (AS) of precursor mRNAs (pre-mRNAs) from multiexon genes allows organisms to increase their coding potential and regulate gene expression through multiple mechanisms. Recent transcriptome-wide analysis of AS using RNA sequencing has revealed that AS is highly pervasive in plants. Pre-mRNAs from over 60% of intron-containing genes undergo AS to produce a vast repertoire of mRNA isoforms. The functions of most splice variants are unknown. However, emerging evidence indicates that splice variants increase the functional diversity of proteins. Furthermore, AS is coupled to transcript stability and translation through nonsense-mediated decay and microRNA-mediated gene regulation. Widespread changes in AS in response to developmental cues and stresses suggest a role for regulated splicing in plant development and stress responses. Here, we review recent progress in uncovering the extent and complexity of the AS landscape in plants, its regulation, and the roles of AS in gene regulation. The prevalence of AS in plants has raised many new questions that require additional studies. New tools based on recent technological advances are allowing genome-wide analysis of RNA elements in transcripts and of chromatin modifications that regulate AS. Application of these tools in plants will provide significant new insights into AS regulation and crosstalk between AS and other layers of gene regulation.

Lariat sequencing in a unicellular yeast identifies regulated alternative splicing of exons that are evolutionarily conserved with humans

Alternative splicing is a potent regulator of gene expression that vastly increases proteomic diversity in multicellular eukaryotes and is associated with organismal complexity. Although alternative splicing is widespread in vertebrates, little is known about the evolutionary origins of this process, in part because of the absence of phylogenetically conserved events that cross major eukaryotic clades. Here we describe a lariat-sequencing approach, which offers high sensitivity for detecting splicing events, and its application to the unicellular fungus, Schizosaccharomyces pombe, an organism that shares many of the hallmarks of alternative splicing in mammalian systems but for which no previous examples of exon-skipping had been demonstrated. Over 200 previously unannotated splicing events were identified, including examples of regulated alternative splicing. Remarkably, an evolutionary analysis of four of the exons identified here as subject to skipping in S. pombe reveals high sequence conservation and perfect length conservation with their homologs in scores of plants, animals, and fungi. Moreover, alternative splicing of two of these exons have been documented in multiple vertebrate organisms, making these the first demonstrations of identical alternative-splicing patterns in species that are separated by over 1 billion y of evolution.

Features of evolutionarily conserved alternative splicing events between Brassica and Arabidopsis

Alternative splicing (AS) generates multiple types of mRNA from a single type of pre-mRNA by differential intron splicing. It can result in new protein isoforms or down-regulation of gene expression by transcript decay. The evolutionary conservation of AS events in plants is largely unexplored and only a small number of AS events have been identified as conserved between divergent species. We performed a large-scale analysis of cDNA data from Brassica and Arabidopsis to identify and further characterize conserved AS events. We identified 537 conserved AS events in 485 genes. Alternative donor and acceptor events are significantly overrepresented among conserved events, whereas intron retention and exon skipping events are underrepresented. Conserved AS events are significantly shorter, less likely to be in the 3'UTR, and they are enriched for genes whose products function in the chloroplast. AS modified a functional domain for about half of the genes with conserved events. We further characterized three genes with conserved AS events. This study identifies many AS events that are conserved between Brassica and Arabidopsis, revealing features of conserved AS events. Many of the conserved AS events may have important, but uncharacterized, functions.

Genome-wide landscape of alternative splicing events in Brachypodium distachyon

Recently, Brachypodium distachyon has emerged as a model plant for studying monocot grasses and cereal crops. Using assembled expressed transcript sequences and subsequent mapping to the corresponding genome, we identified 1219 alternative splicing (AS) events spanning across 2021 putatively assembled transcripts generated from 941 genes. Approximately, 6.3% of expressed genes are alternatively spliced in B. distachyon. We observed that a majority of the identified AS events were related to retained introns (55.5%), followed by alternative acceptor sites (16.7%). We also observed a low percentage of exon skipping (5.0%) and alternative donor site events (8.8%). The 'complex event' that consists of a combination of two or more basic splicing events accounted for ∼14.0%. Comparative AS transcript analysis revealed 163 and 39 homologous pairs between B. distachyon and Oryza sativa and between B. distachyon and Arabidopsis thaliana, respectively. In all, we found 16 AS transcripts to be conserved in all 3 species. AS events and related putative assembled transcripts annotation can be systematically browsed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/plant/).

Alternative splicing in plants--coming of age

More than 60% of intron-containing genes undergo alternative splicing (AS) in plants. This number will increase when AS in different tissues, developmental stages, and environmental conditions are explored. Although the functional impact of AS on protein complexity is still understudied in plants, recent examples demonstrate its importance in regulating plant processes. AS also regulates transcript levels and the link with nonsense-mediated decay and generation of unproductive mRNAs illustrate the need for both transcriptional and AS data in gene expression analyses. AS has influenced the evolution of the complex networks of regulation of gene expression and variation in AS contributed to adaptation of plants to their environment and therefore will impact strategies for improving plant and crop phenotypes.

Rice LGD1 containing RNA binding activity affects growth and development through alternative promoters

Tiller initiation and panicle development are important agronomical traits for grain production in Oryza sativa L. (rice), but their regulatory mechanisms are not yet fully understood. In this study, T-DNA mutant and RNAi transgenic approaches were used to functionally characterize a unique rice gene, LAGGING GROWTH AND DEVELOPMENT 1 (LGD1). The lgd1 mutant showed slow growth, reduced tiller number and plant height, altered panicle architecture and reduced grain yield. The fewer unelongated internodes and cells in lgd1 led to respective reductions in tiller number and to semi-dwarfism. Several independent LGD1-RNAi lines exhibited defective phenotypes similar to those observed in lgd1. Interestingly, LGD1 encodes multiple transcripts with different transcription start sites (TSSs), which were validated by RNA ligase-mediated rapid amplification of 5' and 3' cDNA ends (RLM-RACE). Additionally, GUS assays and a luciferase promoter assay confirmed the promoter activities of LGD1.1 and LGD1.5. LGD1 encoding a von Willebrand factor type A (vWA) domain containing protein is a single gene in rice that is seemingly specific to grasses. GFP-tagged LGD1 isoforms were predominantly detected in the nucleus, and weakly in the cytoplasm. In vitro northwestern analysis showed the RNA-binding activity of the recombinant C-terminal LGD1 protein. Our results demonstrated that LGD1 pleiotropically regulated rice vegetative growth and development through both the distinct spatiotemporal expression patterns of its multiple transcripts and RNA binding activity. Hence, the study of LGD1 will strengthen our understanding of the molecular basis of the multiple transcripts, and their corresponding polypeptides with RNA binding activity, that regulate pleiotropic effects in rice.

Comparative analysis of serine/arginine-rich proteins across 27 eukaryotes: insights into sub-family classification and extent of alternative splicing

Alternative splicing (AS) of pre-mRNA is a fundamental molecular process that generates diversity in the transcriptome and proteome of eukaryotic organisms. SR proteins, a family of splicing regulators with one or two RNA recognition motifs (RRMs) at the N-terminus and an arg/ser-rich domain at the C-terminus, function in both constitutive and alternative splicing. We identified SR proteins in 27 eukaryotic species, which include plants, animals, fungi and "basal" eukaryotes that lie outside of these lineages. Using RNA recognition motifs (RRMs) as a phylogenetic marker, we classified 272 SR genes into robust sub-families. The SR gene family can be split into five major groupings, which can be further separated into 11 distinct sub-families. Most flowering plants have double or nearly double the number of SR genes found in vertebrates. The majority of plant SR genes are under purifying selection. Moreover, in all paralogous SR genes in Arabidopsis, rice, soybean and maize, one of the two paralogs is preferentially expressed throughout plant development. We also assessed the extent of AS in SR genes based on a splice graph approach (http://combi.cs.colostate.edu/as/gmap_SRgenes). AS of SR genes is a widespread phenomenon throughout multiple lineages, with alternative 3' or 5' splicing events being the most prominent type of event. However, plant-enriched sub-families have 57%-88% of their SR genes experiencing some type of AS compared to the 40%-54% seen in other sub-families. The SR gene family is pervasive throughout multiple eukaryotic lineages, conserved in sequence and domain organization, but differs in gene number across lineages with an abundance of SR genes in flowering plants. The higher number of alternatively spliced SR genes in plants emphasizes the importance of AS in generating splice variants in these organisms.

ARTADE2DB: improved statistical inferences for Arabidopsis gene functions and structure predictions by dynamic structure-based dynamic expression (DSDE) analyses

Recent advances in technologies for observing high-resolution genomic activities, such as whole-genome tiling arrays and high-throughput sequencers, provide detailed information for understanding genome functions. However, the functions of 50% of known Arabidopsis thaliana genes remain unknown or are annotated only on the basis of static analyses such as protein motifs or similarities. In this paper, we describe dynamic structure-based dynamic expression (DSDE) analysis, which sequentially predicts both structural and functional features of transcripts. We show that DSDE analysis inferred gene functions 12% more precisely than static structure-based dynamic expression (SSDE) analysis or conventional co-expression analysis based on previously determined gene structures of A. thaliana. This result suggests that more precise structural information than the fixed conventional annotated structures is crucial for co-expression analysis in systems biology of transcriptional regulation and dynamics. Our DSDE method, ARabidopsis Tiling-Array-based Detection of Exons version 2 and over-representation analysis (ARTADE2-ORA), precisely predicts each gene structure by combining two statistical analyses: a probe-wise co-expression analysis of multiple transcriptome measurements and a Markov model analysis of genome sequences. ARTADE2-ORA successfully identified the true functions of about 90% of functionally annotated genes, inferred the functions of 98% of functionally unknown genes and predicted 1,489 new gene structures and functions. We developed a database ARTADE2DB that integrates not only the information predicted by ARTADE2-ORA but also annotations and other functional information, such as phenotypes and literature citations, and is expected to contribute to the study of the functional genomics of A. thaliana. URL: http://artade.org.

A role for SR proteins in plant stress responses

Members of the SR (serine/arginine-rich) protein gene family are key players in the regulation of alternative splicing, an important means of generating proteome diversity and regulating gene expression. In plants, marked changes in alternative splicing are induced by a wide variety of abiotic stresses, suggesting a role for this highly versatile gene regulation mechanism in the response to environmental cues. In support of this notion, the expression of plant SR proteins is stress-regulated at multiple levels, with environmental signals controlling their own alternative splicing patterns, phosphorylation status and subcellular distribution. Most importantly, functional links between these RNA-binding proteins and plant stress tolerance are beginning to emerge, including a role in the regulation of abscisic acid (ABA) signaling. Future identification of the physiological mRNA targets of plant SR proteins holds much promise for the elucidation of the molecular mechanisms underlying their role in the response to abiotic stress.

Genome-wide analysis of alternative splicing in Chlamydomonas reinhardtii

Background

Genome-wide computational analysis of alternative splicing (AS) in several flowering plants has revealed that pre-mRNAs from about 30% of genes undergo AS. Chlamydomonas, a simple unicellular green alga, is part of the lineage that includes land plants. However, it diverged from land plants about one billion years ago. Hence, it serves as a good model system to study alternative splicing in early photosynthetic eukaryotes, to obtain insights into the evolution of this process in plants, and to compare splicing in simple unicellular photosynthetic and non-photosynthetic eukaryotes. We performed a global analysis of alternative splicing in Chlamydomonas reinhardtii using its recently completed genome sequence and all available ESTs and cDNAs.

Results

Our analysis of AS using BLAT and a modified version of the Sircah tool revealed AS of 498 transcriptional units with 611 events, representing about 3% of the total number of genes. As in land plants, intron retention is the most prevalent form of AS. Retained introns and skipped exons tend to be shorter than their counterparts in constitutively spliced genes. The splice site signals in all types of AS events are weaker than those in constitutively spliced genes. Furthermore, in alternatively spliced genes, the prevalent splice form has a stronger splice site signal than the non-prevalent form. Analysis of constitutively spliced introns revealed an over-abundance of motifs with simple repetitive elements in comparison to introns involved in intron retention. In almost all cases, AS results in a truncated ORF, leading to a coding sequence that is around 50% shorter than the prevalent splice form. Using RT-PCR we verified AS of two genes and show that they produce more isoforms than indicated by EST data. All cDNA/EST alignments and splice graphs are provided in a website at http://combi.cs.colostate.edu/as/chlamy.

Conclusions

The extent of AS in Chlamydomonas that we observed is much smaller than observed in land plants, but is much higher than in simple unicellular heterotrophic eukaryotes. The percentage of different alternative splicing events is similar to flowering plants. Prevalence of constitutive and alternative splicing in Chlamydomonas, together with its simplicity, many available public resources, and well developed genetic and molecular tools for this organism make it an excellent model system to elucidate the mechanisms involved in regulated splicing in photosynthetic eukaryotes.

Analysis of multiple occurrences of alternative splicing events in Arabidopsis thaliana using novel sequenced full-length cDNAs

Alternative splicing (AS) is a mechanism by which multiple types of mature mRNAs are generated from a single pre-mature mRNA. In this study, we completely sequenced 1800 full-length cDNAs from Arabidopsis thaliana, which had 5′ and/or 3′ sequences that were previously found to have AS events or alternative transcription start sites. Unexpectedly, these sequences gave us further evidence of AS, as 601 out of 1800 transcripts showed novel AS events. We focused on the combination patterns of multiple AS events within individual genes. Interestingly, some specific AS event combination patterns tended to appear more frequently than expected. The two most common patterns were: (i) alternative donor–0∼12 times of exon skips–alternative acceptor and (ii) several times (∼8) of retained introns. We also found that multiple AS events in a transcript tend to have the same effects concerning the length of the mature mRNA. Our current results are consistent with our previous observations, which showed changes in AS profiles under different conditions, and suggest the involvement of hypothetical cis- and trans-acting factors in the regulation of AS events.

Comparative analysis indicates that alternative splicing in plants has a limited role in functional expansion of the proteome

Background

Alternative splicing (AS) is a widespread phenomenon in higher eukaryotes but the extent to which it leads to functional protein isoforms and to proteome expansion at large is still a matter of debate. In contrast to animal species, for which AS has been studied extensively at the protein and functional level, protein-centered studies of AS in plant species are scarce. Here we investigate the functional impact of AS in dicot and monocot plant species using a comparative approach.

Results

Detailed comparison of AS events in alternative spliced orthologs from the dicot Arabidopsis thaliana and the monocot Oryza sativa (rice) revealed that the vast majority of AS events in both species do not result from functional conservation. Transcript isoforms that are putative targets for the nonsense-mediated decay (NMD) pathway are as likely to contain conserved AS events as isoforms that are translated into proteins. Similar results were obtained when the same comparison was performed between the two more closely related monocot species rice and Zea mays (maize).

Genome-wide computational analysis of functional protein domains encoded in alternatively and constitutively spliced genes revealed that only the RNA recognition motif (RRM) is overrepresented in alternatively spliced genes in all species analyzed. In contrast, three domain types were overrepresented in constitutively spliced genes. AS events were found to be less frequent within than outside predicted protein domains and no domain type was found to be enriched with AS introns. Analysis of AS events that result in the removal of complete protein domains revealed that only a small number of domain types is spliced-out in all species analyzed. Finally, in a substantial fraction of cases where a domain is completely removed, this domain appeared to be a unit of a tandem repeat.

Conclusion

The results from the ortholog comparisons suggest that the ability of a gene to produce more than one functional protein through AS does not persist during evolution. Cross-species comparison of the results of the protein-domain oriented analyses indicates little correspondence between the analyzed species. Based on the premise that functional genetic features are most likely to be conserved during evolution, we conclude that AS has only a limited role in functional expansion of the proteome in plants.

Highly diversified molecular evolution of downstream transcription start sites in rice and Arabidopsis

Alternative usage of transcription start sites (TSSs) is one of the key mechanisms to generate gene variation in eukaryotes. Here, we show diversified molecular evolution of TSSs in remotely related flowering plants, rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), by comprehensive analyses of large collections of full-length cDNAs and genome sequences. We determined 45,917 representative TSSs within 23,445 loci of rice and 35,313 TSSs within 16,964 loci of Arabidopsis, about two TSSs per locus in either species. The nucleotide features around TSSs displayed distinct patterns when the most upstream TSSs were compared with downstream TSSs. We found that CG-skew and AT-skew were clearly different between upstream and downstream TSSs, and that this difference was commonly observed in rice and Arabidopsis. Relative entropy analysis revealed that the most upstream TSSs had retained canonical cis elements, whereas downstream TSSs showed atypical nucleotide features. Expression patterns were distinguishable between upstream and downstream TSSs. These results indicate that plant TSSs were generally diversified in downstream regions, resulting in the development of new gene expression patterns. Furthermore, our comparative analysis of TSS variation between the species showed a positive correlation between TSS number and gene conservation. Rice and Arabidopsis might have evolved novel TSSs in an independent manner, which led to diversification of these two species.

Spliceosomal proteins in plants

The spliceosome is a large nuclear structure consisting of dynamically interacting RNAs and proteins. This chapter briefly reviews some of the known components and their interactions. Large-scale proteomics and gene expression studies may be required to unravel the many intricate mechanisms involved in splice site recognition and selection.

Regulation of alternative splicing of pre-mRNAs by stresses

A substantial fraction (approximately 30%) of plant genes is alternatively spliced, but how alternative splicing is regulated remains unknown. Many plant genes undergo alternative splicing in response to a variety of stresses. Large-scale computational analyses and experimental approaches focused on select genes are beginning to reveal that alternative splicing constitutes an integral part of gene regulation in stress responses. Based on the studies discussed in this chapter, it appears that alternative splicing generates transcriptome/proteome complexity that is likely to be important for stress adaptation. However, the signaling pathways that relay stress conditions to splicing machinery and if and how the alternative spliced products confer adaptive advantages to plants are poorly understood.

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.

Cross-species EST alignments reveal novel and conserved alternative splicing events in legumes

BACKGROUND

Although originally thought to be less frequent in plants than in animals, alternative splicing (AS) is now known to be widespread in plants. Here we report the characteristics of AS in legumes, one of the largest and most important plant families, based on EST alignments to the genome sequences of Medicago truncatula (Mt) and Lotus japonicus (Lj).

RESULTS

Based on cognate EST alignments alone, the observed frequency of alternatively spliced genes is lower in Mt (approximately 10%, 1,107 genes) and Lj (approximately 3%, 92 genes) than in Arabidopsis and rice (both around 20%). However, AS frequencies are comparable in all four species if EST levels are normalized. Intron retention is the most common form of AS in all four plant species (~50%), with slightly lower frequency in legumes compared to Arabidopsis and rice. This differs notably from vertebrates, where exon skipping is most common. To uncover additional AS events, we aligned ESTs from other legume species against the Mt genome sequence. In this way, 248 additional Mt genes were predicted to be alternatively spliced. We also identified 22 AS events completely conserved in two or more plant species.

CONCLUSION

This study extends the range of plant taxa shown to have high levels of AS, confirms the importance of intron retention in plants, and demonstrates the utility of using ESTs from related species in order to identify novel and conserved AS events. The results also indicate that the frequency of AS in plants is comparable to that observed in mammals. Finally, our results highlight the importance of normalizing EST levels when estimating the frequency of alternative splicing.

Plant SR proteins and their functions

SR proteins are a family of splicing factors important for splice site recognition and spliceosome assembly. Their ability to bind to RNA and to interact with proteins as well identifies them as important players in splice site choice and alternative splicing. Plants possess twice as many SR proteins as animals, and some of the subfamilies are plant specific. Arabidopsis SR proteins are involved in different aspects of plant growth and development as well as in responses to environmental cues. The plant-specific subfamilies have been shown to be regulated by alternative splicing events, which are highly conserved in evolution. The tight regulation of splicing factors by alternative splicing might allow coordinated responses of their target genes.

Comparative cross-species alternative splicing in plants

Alternative splicing (AS) can add significantly to genome complexity. Plants are thought to exhibit less AS than animals. An algorithm, based on expressed sequence tag (EST) pairs gapped alignment, was developed that takes advantage of the relatively small intron and exon size in plants and directly compares pairs of ESTs to search for AS. EST pairs gapped alignment was first evaluated in Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and tomato (Solanum lycopersicum) for which annotated genome sequence is available and was shown to accurately predict splicing events. The method was then applied to 11 plant species that include 17 cultivars for which enough ESTs are available. The results show a large, 3.7-fold difference in AS rates between plant species with Arabidopsis and rice in the lower range and lettuce (Lactuca sativa) and sorghum (Sorghum bicolor) in the upper range. Hence, compared to higher animals, plants show a much greater degree of variety in their AS rates and in some plant species the rates of animal and plant AS are comparable although the distribution of AS types may differ. In eudicots but not monocots, a correlation between genome size and AS rates was detected, implying that in eudicots the mechanisms that lead to larger genomes are a driving force for the evolution of AS.

Regulation of plant developmental processes by a novel splicing factor

Serine/arginine-rich (SR) proteins play important roles in constitutive and alternative splicing and other aspects of mRNA metabolism. We have previously isolated a unique plant SR protein (SR45) with atypical domain organization. However, the biological and molecular functions of this novel SR protein are not known. Here, we report biological and molecular functions of this protein. Using an in vitro splicing complementation assay, we showed that SR45 functions as an essential splicing factor. Furthermore, the alternative splicing pattern of transcripts of several other SR genes was altered in a mutant, sr45-1, suggesting that the observed phenotypic abnormalities in sr45-1 are likely due to altered levels of SR protein isoforms, which in turn modulate splicing of other pre-mRNAs. sr45-1 exhibited developmental abnormalities, including delayed flowering, narrow leaves and altered number of petals and stamens. The late flowering phenotype was observed under both long days and short days and was rescued by vernalization. FLC, a key flowering repressor, is up-regulated in sr45-1 demonstrating that SR45 influences the autonomous flowering pathway. Changes in the alternative splicing of SR genes and the phenotypic defects in the mutant were rescued by SR45 cDNA, further confirming that the observed defects in the mutant are due to the lack of SR45. These results indicate that SR45 is a novel plant-specific splicing factor that plays a crucial role in regulating developmental processes.

Unproductive splicing of SR genes associated with highly conserved and ultraconserved DNA elements

The human and mouse genomes share a number of long, perfectly conserved nucleotide sequences, termed ultraconserved elements. Whereas these regions can act as transcriptional enhancers when upstream of genes, those within genes are less well understood. In particular, the function of ultraconserved elements that overlap alternatively spliced exons of genes encoding RNA-binding proteins is unknown. Here we report that in every member of the human SR family of splicing regulators, highly or ultraconserved elements are alternatively spliced, either as alternative 'poison cassette exons' containing early in-frame stop codons, or as alternative introns in the 3' untranslated region. These alternative splicing events target the resulting messenger RNAs for degradation by means of an RNA surveillance pathway called nonsense-mediated mRNA decay. Mouse orthologues of the human SR proteins exhibit the same unproductive splicing patterns. Three SR proteins have been previously shown to direct splicing of their own transcripts, and one of these is known to autoregulate its expression by coupling alternative splicing with decay; our results suggest that unproductive splicing is important for regulation of the entire SR family. We find that unproductive splicing associated with conserved regions has arisen independently in different SR genes, suggesting that splicing factors may readily acquire this form of regulation.

Alternative splicing of pre-mRNAs of Arabidopsis serine/arginine-rich proteins: regulation by hormones and stresses

Precursor mRNAs with introns can undergo alternative splicing (AS) to produce structurally and functionally different proteins from the same gene. Here, we show that the pre-mRNAs of Arabidopsis genes that encode serine/arginine-rich (SR) proteins, a conserved family of splicing regulators in eukaryotes, are extensively alternatively spliced. Remarkably about 95 transcripts are produced from only 15 genes, thereby increasing the complexity of the SR gene family transcriptome by six-fold. The AS of some SR genes is controlled in a developmental and tissue-specific manner. Interestingly, among the various hormones and abiotic stresses tested, temperature stress (cold and heat) dramatically altered the AS of pre-mRNAs of several SR genes, whereas hormones altered the splicing of only three SR genes. These results indicate that abiotic stresses regulate the AS of the pre-mRNAs of SR genes to produce different isoforms of SR proteins that are likely to have altered function(s) in pre-mRNA splicing. Sequence analysis of splice variants revealed that predicted proteins from a majority of these variants either lack one or more modular domains or contain truncated domains. Because of the modular nature of the various domains in SR proteins, the proteins produced from splice variants are likely to have distinct functions. Together our results indicate that Arabidopsis SR genes generate surprisingly large transcriptome complexity, which is altered by stresses and hormones.