Characterization and comparison of intron structure and alternative splicing between Medicago truncatula, Populus trichocarpa, Arabidopsis and rice

Regulation of Flowering Time and Other Developmental Plasticities by 3' Splicing Factor-Mediated Alternative Splicing in Arabidopsis thaliana

Plants, as sessile organisms, show a high degree of plasticity in their growth and development and have various strategies to cope with these alterations under continuously changing environments and unfavorable stress conditions. In particular, the floral transition from the vegetative and reproductive phases in the shoot apical meristem (SAM) is one of the most important developmental changes in plants. In addition, meristem regions, such as the SAM and root apical meristem (RAM), which continually generate new lateral organs throughout the plant life cycle, are important sites for developmental plasticity. Recent findings have shown that the prevailing type of alternative splicing (AS) in plants is intron retention (IR) unlike in animals; thus, AS is an important regulatory mechanism conferring plasticity for plant growth and development under various environmental conditions. Although eukaryotes exhibit some similarities in the composition and dynamics of their splicing machinery, plants have differences in the 3' splicing characteristics governing AS. Here, we summarize recent findings on the roles of 3' splicing factors and their interacting partners in regulating the flowering time and other developmental plasticities in Arabidopsis thaliana.

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

Genome-Wide Identification and Characterization of the Trehalose-6-phosphate Synthetase (TPS) Gene Family in Watermelon (Citrullus lanatus) and Their Transcriptional Responses to Salt Stress

With the increase in watermelon cultivation area, there is an urgent need to explore enzymatic and genetic resources for the sustainable development of watermelon, especially under salt stress. Among the various compounds known, trehalose plays an important role in regulating abiotic stress tolerances in diverse organisms, including plants. Therefore, the present study comprehensively analyzed the trehalose-6-phosphate synthase (TPS) gene family in watermelon. The study analyzed the functional classification, evolutionary characteristics, and expression patterns of the watermelon TPS genes family. Seven ClTPSs were identified and classified into two distinct classes according to gene structure and phylogeny. Evolutionary analysis suggested the role of purifying selection in the evolution of the TPS family members. Further, cis-acting elements related to plant hormones and abiotic stress were identified in the promoter region of the TPS genes. The tissue-specific expression analysis showed that ClTPS genes were widely expressed in roots, stems, leaves, flowers, and fruits, while ClTPS3 was significantly induced under salt stress. The overexpression of ClTPS3 in Arabidopsis thaliana significantly improved salt tolerance. Finally, the STRING functional protein association networks suggested that the transcription factor ClMYB and ClbHLH regulate ClTPS3. Thus, the study indicates the critical role of ClTPS3 in watermelon response to salt stress.

PSDX: A Comprehensive Multi-Omics Association Database of Populus trichocarpa With a Focus on the Secondary Growth in Response to Stresses

Populus trichocarpa (P. trichocarpa) is a model tree for the investigation of wood formation. In recent years, researchers have generated a large number of high-throughput sequencing data in P. trichocarpa. However, no comprehensive database that provides multi-omics associations for the investigation of secondary growth in response to diverse stresses has been reported. Therefore, we developed a public repository that presents comprehensive measurements of gene expression and post-transcriptional regulation by integrating 144 RNA-Seq, 33 ChIP-seq, and six single-molecule real-time (SMRT) isoform sequencing (Iso-seq) libraries prepared from tissues subjected to different stresses. All the samples from different studies were analyzed to obtain gene expression, co-expression network, and differentially expressed genes (DEG) using unified parameters, which allowed comparison of results from different studies and treatments. In addition to gene expression, we also identified and deposited pre-processed data about alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI). The post-transcriptional regulation, differential expression, and co-expression network datasets were integrated into a new P. trichocarpa Stem Differentiating Xylem (PSDX) database (http://forestry.fafu.edu.cn/db/SDX), which further highlights gene families of RNA-binding proteins and stress-related genes. The PSDX also provides tools for data query, visualization, a genome browser, and the BLAST option for sequence-based query. Much of the data is also available for bulk download. The availability of PSDX contributes to the research related to the secondary growth in response to stresses in P. trichocarpa, which will provide new insights that can be useful for the improvement of stress tolerance in woody plants.

Genome-wide profiling of alternative splicing genes in hybrid poplar (P.alba×P.glandulosa cv.84K) leaves

Alternative splicing (AS) is a post-transcriptional process common in plants and essential for regulation of environmental fitness of plants. In the present study, we focus on the AS events in poplar leaves to understand their effects on plant growth and development. The hybrid poplar (P.alba×P.glandulosa cv.84K) leaves were collected for RNA extraction. The extracted RNA was sequenced using on an Illumina HiSeq™ 2000 platform. Using the Populus trichocarpa genome as the reference, a total of 3810 AS genes were identified (9225 AS events), which accounted for 13.51% of all the expressed genes. Intron retention was the most common AS event, accounting for 43.86% of all the AS events, followed by alternative 3′ splice sites (23.75%), alternative 5′ splice sites (23.71%), and exon skipping (8.68%). Chromosomes 10 had the most condensed AS events (33.67 events/Mb) and chromosome 19 had the least (12.42 events/Mb). Association analysis showed that AS in the poplar leaves was positively correlated with intron length, exon number, exon length, and gene expression level, and was negatively correlated with GC content. AS genes in the poplar leaves were associated mainly with inositol phosphate metabolism and phosphatidylinositol signaling system pathways that would be significant on wooden plant production.

The interplay between microRNA and alternative splicing of linear and circular RNAs in eleven plant species

MOTIVATION

MicroRNA (miRNA) and alternative splicing (AS)-mediated post-transcriptional regulation has been extensively studied in most eukaryotes. However, the interplay between AS and miRNAs has not been explored in plants. To our knowledge, the overall profile of miRNA target sites in circular RNAs (circRNA) generated by alternative back splicing has never been reported previously. To address the challenge, we identified miRNA target sites located in alternatively spliced regions of the linear and circular splice isoforms using the up-to-date single-molecule real-time (SMRT) isoform sequencing (Iso-Seq) and Illumina sequencing data in eleven plant species.

RESULTS

In total, we identified 399 401 and 114 574 AS events from linear and circular RNAs, respectively. Among them, there were 64 781 and 41 146 miRNA target sites located in linear and circular AS region, respectively. In addition, we found 38 913 circRNAs to be overlapping with 45 648 AS events of its own parent isoforms, suggesting circRNA regulation of AS of linear RNAs by forming R-loop with the genomic locus. Here, we present a comprehensive database of miRNA targets in alternatively spliced linear and circRNAs (ASmiR) and a web server for deposition and identification of miRNA target sites located in the alternatively spliced region of linear and circular RNAs. This database is accompanied by an easy-to-use web query interface for meaningful downstream analysis. Plant research community can submit user-defined datasets to the web service to search AS regions harboring small RNA target sites. In conclusion, this study provides an unprecedented resource to understand regulatory relationships between miRNAs and AS in both gymnosperms and angiosperms.

AVAILABILITY AND IMPLEMENTATION

The readily accessible database and web-based tools are available at http://forestry.fafu.edu.cn/bioinfor/db/ASmiR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

Genome-wide identification and comparative analysis of alternative splicing across four legume species

MAIN CONCLUSION

Alternative splicing EVENTS were genome-wide identified for four legume species, and nitrogen fixation-related gene families and evolutionary analysis was also performed. Alternative splicing (AS) is a key regulatory mechanism that contributes to transcriptome and proteome diversity. Investigation of the genome-wide conserved AS events across different species will help with the understanding of the evolution of the functional diversity in legumes, allowing for genetic improvement. Genome-wide identification and characterization of AS were performed using the publically available mRNA, EST, and RNA-Seq data for four important legume species. A total of 15,165 AS genes in Glycine max, 6077 in Cicer arietinum, 7240 in Medicago truncatula, and 7358 in Lotus japonicus were identified. Intron retention (IntronR) was the dominant AS type among the identified events, with IntronR occurring from 53.76% in M. truncatula to 43.91% in C. arietinum. We identified 1159 AS genes that were conserved among four species. Furthermore, nine nitrogen fixation-related gene families with 237 genes were identified, and 80 of them were AS, accounting for the 43.48% in G. max and 27.78% in C. arietinum. An evolutionary analysis showed that these AS genes tended to be located adjacent to each other in the evolutionary tree and are unbalanced in the distribution in the sub-family. This study provides a foundation for future studies on transcription complexity, evolution, and the role of AS on plant functional regulation.

Computational analysis of alternative splicing in plant genomes

Computational analyses play crucial roles in characterizing splicing isoforms in plant genomes. In this review, we provide a survey of computational tools used in recently published, genome-scale splicing analyses in plants. We summarize the commonly used software and pipelines for read mapping, isoform reconstruction, isoform quantification, and differential expression analysis. We also discuss methods for analyzing long reads and the strategies to combine long and short reads in identifying splicing isoforms. We review several tools for characterizing local splicing events, splicing graphs, coding potential, and visualizing splicing isoforms. We further discuss the procedures for identifying conserved splicing isoforms across plant species. Finally, we discuss the outlook of integrating other genomic data with splicing analyses to identify regulatory mechanisms of AS on genome-wide scale.

Identification and characterization of evolutionarily conserved alternative splicing events in a mangrove genus Sonneratia

Alternative splicing (AS), which produces multiple mRNA transcripts from a single gene, plays crucial roles in plant growth, development and environmental stress responses. Functional significances of conserved AS events among congeneric species have not been well characterized. In this study, we performed transcriptome sequencing to characterize AS events in four common species of Sonneratia, a mangrove genus excellently adaptive to intertidal zones. 7,248 to 12,623 AS events were identified in approximately 25% to 35% expressed genes in the roots of the four species. The frequency of AS events in Sonneratia was associated with genomic features, including gene expression level and intron/exon number and length. Among the four species, 1,355 evolutionarily conserved AS (ECAS) events were identified from 1,170 genes. Compared with non-ECAS events, ECAS events are of shorter length and less possibility to introduce premature stop codons (PTCs) and frameshifts. Functional annotations of the genes containing ECAS events showed that four of the 26 enriched Gene Ontology (GO) terms are involved in proton transport, signal transduction and carbon metabolism, and 60 genes from another three GO terms are implicated in responses to osmotic, oxidative and heat stresses, which may contribute to the adaptation of Sonneratia species to harsh intertidal environments.

Wild Carrot Differentiation in Europe and Selection at DcAOX1 Gene?

By definition, the domestication process leads to an overall reduction of crop genetic diversity. This lead to the current search of genomic regions in wild crop relatives (CWR), an important task for modern carrot breeding. Nowadays massive sequencing possibilities can allow for discovery of novel genetic resources in wild populations, but this quest could be aided by the use of a surrogate gene (to first identify and prioritize novel wild populations for increased sequencing effort). Alternative oxidase (AOX) gene family seems to be linked to all kinds of abiotic and biotic stress reactions in various organisms and thus have the potential to be used in the identification of CWR hotspots of environment-adapted diversity. High variability of DcAOX1 was found in populations of wild carrot sampled across a West-European environmental gradient. Even though no direct relation was found with the analyzed climatic conditions or with physical distance, population differentiation exists and results mainly from the polymorphisms associated with DcAOX1 exon 1 and intron 1. The relatively high number of amino acid changes and the identification of several unusually variable positions (through a likelihood ratio test), suggests that DcAOX1 gene might be under positive selection. However, if positive selection is considered, it only acts on some specific populations (i.e. is in the form of adaptive differences in different population locations) given the observed high genetic diversity. We were able to identify two populations with higher levels of differentiation which are promising as hot spots of specific functional diversity.

The origin of exon 3 skipping of paternal GLOBOSA pre-mRNA in some Nicotiana tabacum lines correlates with a point mutation of the very last nucleotide of the exon

In plants, genome duplication followed by genome diversification and selection is recognized as a major evolutionary process. Rapid epigenetic and genetic changes that affect the transcription of parental genes are frequently observed after polyploidization. The pattern of alternative splicing is also frequently altered, yet the related molecular processes remain largely unresolved. Here, we study the inheritance and expression of parental variants of three floral organ identity genes in allotetraploid tobacco. DEFICIENS and GLOBOSA are B-class genes, and AGAMOUS is a C-class gene. Parental variants of these genes were found to be maintained in the tobacco genome, and the respective mRNAs were present in flower buds in comparable amounts. However, among five tobacco cultivars, we identified two in which the majority of paternal GLOBOSA pre-mRNA transcripts undergo exon 3 skipping, producing an mRNA with a premature termination codon. At the DNA level, we identified a G-A transition at the very last position of exon 3 in both cultivars. Although alternative splicing resulted in a dramatic decrease in full-length paternal GLOBOSA mRNA, no phenotypic effect was observed. Our finding likely serves as an example of the initiation of homoeolog diversification in a relatively young polyploid genome.

Genome-wide cataloging and analysis of alternatively spliced genes in cereal crops

Background

Protein functional diversity at the post-transcriptional level is regulated through spliceosome mediated pre-mRNA alternative splicing (AS) events and that has been widely demonstrated to be a key player in regulating the functional diversity in plants. Identification and analysis of AS genes in cereal crop plants are critical for crop improvement and understanding regulatory mechanisms.

Results

We carried out the comparative analyses of the functional landscapes of the AS using the consensus assembly of expressed sequence tags and available mRNA sequences in four cereal plants. We identified a total of 8,734 in Oryza sativa subspecies (ssp) japonica, 2,657 in O. sativa ssp indica, 3,971 in Sorghum bicolor, and 10,687 in Zea mays AS genes. Among the identified AS events, intron retention remains to be the dominant type accounting for 23.5 % in S. bicolor, and up to 55.8 % in O. sativa ssp indica. We identified a total of 887 AS genes that were conserved among Z. mays, S. bicolor, and O. sativa ssp japonica; and 248 AS genes were found to be conserved among all four studied species or ssp. Furthermore, we identified 53 AS genes conserved with Brachypodium distachyon. Gene Ontology classification of AS genes revealed functional assignment of these genes in many biological processes with diverse molecular functions.

Conclusions

AS is common in cereal plants. The AS genes identified in four cereal crops in this work provide the foundation for further studying the roles of AS in regulation of cereal plant growth and development. The data can be accessed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1914-5) contains supplementary material, which is available to authorized users.

Cauliflower mosaic virus Transcriptome Reveals a Complex Alternative Splicing Pattern

The plant pararetrovirus Cauliflower mosaic virus (CaMV) uses alternative splicing to generate several isoforms from its polycistronic pregenomic 35S RNA. This pro-cess has been shown to be essential for infectivity. Previous works have identified four splice donor sites and a single splice acceptor site in the 35S RNA 5' region and suggested that the main role of CaMV splicing is to downregulate expression of open reading frames (ORFs) I and II. In this study, we show that alternative splicing is a conserved process among CaMV isolates. In Cabb B-JI and Cabb-S isolates, splicing frequently leads to different fusion between ORFs, particularly between ORF I and II. The corresponding P1P2 fusion proteins expressed in E. coli interact with viral proteins P2 and P3 in vitro. However, they are detected neither during infection nor upon transient expression in planta, which suggests rapid degradation after synthesis and no important biological role in the CaMV infectious cycle. To gain a better understanding of the functional relevance of 35S RNA alternative splicing in CaMV infectivity, we inactivated the previously described splice sites. All the splicing mutants were as pathogenic as the corresponding wild-type isolate. Through RT-PCR-based analysis we demonstrate that CaMV 35S RNA exhibits a complex splicing pattern, as we identify new splice donor and acceptor sites whose selection leads to more than thirteen 35S RNA isoforms in infected turnip plants. Inactivating splice donor or acceptor sites is not lethal for the virus, since disrupted sites are systematically rescued by the activation of cryptic and/or seldom used splice sites. Taken together, our data depict a conserved, complex and flexible process, involving multiple sites, that ensures splicing of 35S RNA.

Genome-wide identification of evolutionarily conserved alternative splicing events in flowering plants

Alternative splicing (AS) plays important roles in many plant functions, but its conservation across the plant kingdom is not known. We describe a methodology to identify AS events and identify conserved AS events across large phylogenetic distances using RNA-Seq datasets. We applied this methodology to transcriptome data from nine angiosperms including Amborella, the single sister species to all other extant flowering plants. AS events within 40–70% of the expressed multi-exonic genes per species were found, 27,120 of which are conserved among two or more of the taxa studied. While many events are species specific, many others are shared across long evolutionary distances suggesting they have functional significance. Conservation of AS event data provides an estimate of the number of ancestral AS events present at each node of the tree representing the nine species studied. Furthermore, the presence or absence of AS isoforms between species with different whole genome duplication (WGD) histories provides the opportunity to examine the impact of WDG on AS potential. Examining AS in gene families identifies those with high rates of AS, and conservation can distinguish ancient events vs. recent or species specific adaptations. The MADS-box and SR protein families are found to represent families with low and high occurrences of AS, respectively, yet their AS events were likely present in the MRCA of angiosperms.

Genome-wide analysis of alternative splicing in Volvox carteri

Background

Alternative splicing is an essential mechanism for increasing transcriptome and proteome diversity in eukaryotes. Particularly in multicellular eukaryotes, this mechanism is involved in the regulation of developmental and physiological processes like growth, differentiation and signal transduction.

Results

Here we report the genome-wide analysis of alternative splicing in the multicellular green alga Volvox carteri. The bioinformatic analysis of 132,038 expressed sequence tags (ESTs) identified 580 alternative splicing events in a total of 426 genes. The predominant type of alternative splicing in Volvox is intron retention (46.5%) followed by alternative 5′ (17.9%) and 3′ (21.9%) splice sites and exon skipping (9.5%). Our analysis shows that in Volvox at least ~2.9% of the intron-containing genes are subject to alternative splicing. Considering the total number of sequenced ESTs, the Volvox genome seems to provide more favorable conditions (e.g., regarding length and GC content of introns) for the occurrence of alternative splicing than the genome of its close unicellular relative Chlamydomonas. Moreover, many randomly chosen alternatively spliced genes of Volvox do not show alternative splicing in Chlamydomonas. Since the Volvox genome contains about the same number of protein-coding genes as the Chlamydomonas genome (~14,500 protein-coding genes), we assumed that alternative splicing may play a key role in generation of genomic diversity, which is required to evolve from a simple one-cell ancestor to a multicellular organism with differentiated cell types (Mol Biol Evol 31:1402-1413, 2014). To confirm the alternative splicing events identified by bioinformatic analysis, several genes with different types of alternatively splicing have been selected followed by experimental verification of the predicted splice variants by RT-PCR.

Conclusions

The results show that our approach for prediction of alternative splicing events in Volvox was accurate and reliable. Moreover, quantitative real-time RT-PCR appears to be useful in Volvox for analyses of relationships between the appearance of specific alternative splicing variants and different kinds of physiological, metabolic and developmental processes as well as responses to environmental changes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1117) contains supplementary material, which is available to authorized users.

Genome-wide survey of Alternative Splicing in Sorghum Bicolor

Sorghum bicolor is a member of grass family which is an attractive model plant for genome study due to interesting genome features like low genome size. In this research, we performed comprehensive investigation of Alternative Splicing and ontology aspects of genes those have undergone these events in sorghum bicolor. We used homology based alignments between gene rich transcripts, represented by tentative consensus (TC) transcript sequences, and genomic scaffolds to deduce the structure of genes and identify alternatively spliced transcripts in sorghum. Using homology mapping of assembled expressed sequence tags with genomics data, we identified 2,137 Alternative Splicing events in S. bicolor. Our study showed that complex events and intron retention are the main types of Alternative Splicing events in S. bicolor and highlights the prevalence of splicing site recognition for definition of introns in this plant. Annotations of the alternatively spliced genes revealed that they represent diverse biological process and molecular functions, suggesting a fundamental role for Alternative Splicing in affecting the development and physiology of S. bicolor.

The extent of Ds1 transposon to enrich transcriptomes and proteomes by exonization

BACKGROUND

Exonization is an event which an intronic transposed element (TE) provides splice sites and leads to alternatively spliced cassette exons. Without disrupting of the inserted gene's function, TEs can expand the proteome diversity by adding the splice variant that encodes a different, yet functional protein. Previously, we found that the main contribution of Ds exonization for gene divergence is not providing genetic messages but incorporating the intron sequences with different reading frame patterns to enrich the plant proteome. Ds1, another member of Ac/Ds transposon system, differs from Ds by providing 3 splice donor sites and 2 acceptor sites for alternative splicing, which may greatly increase the extent for proteome expansion.

RESULTS

In this study, we performed a genome-wide survey of Ds1 exonization events to assess its extent to enrich proteomes in plants. Each Ds1 insertion yielded 11 transcript isoforms by integrating the splice donor and/or acceptor sites, which composed a bulk of all exonized transcript orthologs from the dicot Arabidopsis thaliana and the monocot Oryza sativa (rice). The exonized transcripts were analyzed by the locations of the termination codon (PTC) and the putative targets for the nonsense-mediated decay (NMD) pathway were then excluded. Compared with the Ds element, Ds1 harbors more contents of non-NMD transcripts for protein isoforms.

CONCLUSIONS

The contribution of Ds1 exonization for gene divergence is incorporating the intron sequences with different reading frame patterns to enrich the plant proteome. All these simulation results direct new experimental analysis at the molecular level.

The SERRATE protein is involved in alternative splicing in Arabidopsis thaliana

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

Comprehensive analysis of RNA-seq data reveals the complexity of the transcriptome in Brassica rapa

Background

The species Brassica rapa (2n=20, AA) is an important vegetable and oilseed crop, and serves as an excellent model for genomic and evolutionary research in Brassica species. With the availability of whole genome sequence of B. rapa, it is essential to further determine the activity of all functional elements of the B. rapa genome and explore the transcriptome on a genome-wide scale. Here, RNA-seq data was employed to provide a genome-wide transcriptional landscape and characterization of the annotated and novel transcripts and alternative splicing events across tissues.

Results

RNA-seq reads were generated using the Illumina platform from six different tissues (root, stem, leaf, flower, silique and callus) of the B. rapa accession Chiifu-401-42, the same line used for whole genome sequencing. First, these data detected the widespread transcription of the B. rapa genome, leading to the identification of numerous novel transcripts and definition of 5'/3' UTRs of known genes. Second, 78.8% of the total annotated genes were detected as expressed and 45.8% were constitutively expressed across all tissues. We further defined several groups of genes: housekeeping genes, tissue-specific expressed genes and co-expressed genes across tissues, which will serve as a valuable repository for future crop functional genomics research. Third, alternative splicing (AS) is estimated to occur in more than 29.4% of intron-containing B. rapa genes, and 65% of them were commonly detected in more than two tissues. Interestingly, genes with high rate of AS were over-represented in GO categories relating to transcriptional regulation and signal transduction, suggesting potential importance of AS for playing regulatory role in these genes. Further, we observed that intron retention (IR) is predominant in the AS events and seems to preferentially occurred in genes with short introns.

Conclusions

The high-resolution RNA-seq analysis provides a global transcriptional landscape as a complement to the B. rapa genome sequence, which will advance our understanding of the dynamics and complexity of the B. rapa transcriptome. The atlas of gene expression in different tissues will be useful for accelerating research on functional genomics and genome evolution in Brassica species.

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.

BARE retrotransposons are translated and replicated via distinct RNA pools

The replication of Long Terminal Repeat (LTR) retrotransposons, which can constitute over 80% of higher plant genomes, resembles that of retroviruses. A major question for retrotransposons and retroviruses is how the two conflicting roles of their transcripts, in translation and reverse transcription, are balanced. Here, we show that the BARE retrotransposon, despite its organization into just one open reading frame, produces three distinct classes of transcripts. One is capped, polyadenylated, and translated, but cannot be copied into cDNA. The second is not capped or polyadenylated, but is destined for packaging and ultimate reverse transcription. The third class is capped, polyadenylated, and spliced to favor production of a subgenomic RNA encoding only Gag, the protein forming virus-like particles. Moreover, the BARE2 subfamily, which cannot synthesize Gag and is parasitic on BARE1, does not produce the spliced sub-genomic RNA for translation but does make the replication competent transcripts, which are packaged into BARE1 particles. To our knowledge, this is first demonstration of distinct RNA pools for translation and transcription for any retrotransposon.

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.

The developing xylem transcriptome and genome-wide analysis of alternative splicing in Populus trichocarpa (black cottonwood) populations

Background

Alternative splicing (AS) of genes is an efficient means of generating variation in protein structure and function. AS variation has been observed between tissues, cell types, and different treatments in non-woody plants such as Arabidopsis thaliana (Arabidopsis) and rice. However, little is known about AS patterns in wood-forming tissues and how much AS variation exists within plant populations.

Results

Here we used high-throughput RNA sequencing to analyze the Populus trichocarpa (P. trichocarpa) xylem transcriptome in 20 individuals from different populations across much of its range in western North America. Deep transcriptome sequencing and mapping of reads to the P. trichocarpa reference genome identified a suite of xylem-expressed genes common to all accessions. Our analysis suggests that at least 36% of the xylem-expressed genes in P. trichocarpa are alternatively spliced. Extensive AS was observed in cell-wall biosynthesis related genes such as glycosyl transferases and C2H2 transcription factors. 27902 AS events were documented and most of these events were not conserved across individuals. Differences in isoform-specific read densities indicated that 7% and 13% of AS events showed significant differences between individuals within geographically separated southern and northern populations, a level that is in general agreement with AS variation in human populations.

Conclusions

This genome-wide analysis of alternative splicing reveals high levels of AS in P. trichocarpa and extensive inter-individual AS variation. We provide the most comprehensive analysis of AS in P. trichocarpa to date, which will serve as a valuable resource for the plant community to study transcriptome complexity and AS regulation during wood formation.

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/).

Genome-wide survey of ds exonization to enrich transcriptomes and proteomes in plants

Insertion of transposable elements (TEs) into introns can lead to their activation as alternatively spliced cassette exons, an event called exonization which can enrich the complexity of transcriptomes and proteomes. Previously, we performed the first experimental assessment of TE exonization by inserting a Ds element into each intron of the rice epsps gene. Exonization of Ds in plants was biased toward providing splice donor sites from the beginning of the inserted Ds sequence. Additionally, Ds inserted in the reverse direction resulted in a continuous splice donor consensus region by offering 4 donor sites in the same intron. The current study involved genome-wide computational analysis of Ds exonization events in the dicot Arabidopsis thaliana and the monocot Oryza sativa (rice). Up to 71% of the exonized transcripts were putative targets for the nonsense-mediated decay (NMD) pathway. The insertion patterns of Ds and the polymorphic splice donor sites increased the transcripts and subsequent protein isoforms. Protein isoforms contain protein sequence due to unspliced intron-TE region and/or a shift of the reading frame. The number of interior protein isoforms would be twice that of C-terminal isoforms, on average. TE exonization provides a promising way for functional expansion of the plant proteome.

Alternatively spliced transcripts of group 3 late embryogenesis abundant protein from Pogonatherum paniceum confer different abiotic stress tolerance in Escherichia coli

The gene encoding a group 3 late embryogenesis abundant (LEA) protein was cloned from callus of the drought-tolerant grass Pogonatherum paniceum. Three alternatively spliced transcripts of this gene were amplified by RT-PCR. According to the bioinformatics analysis, the gene contained three exons and two introns. The PpLEA3.1 transcript which was the most abundant one in P. paniceum contained all three exons, and the PpLEA3.2 transcript lacked fragments of the first two exons which encoded a 41-amino acid-long region of the PpLEA3 protein. The PpLEA3.3 transcript retained the second intron. The three splicing patterns resulted in changes in the number of repeats of an 11-amino acid motif, hydropathy and the predicted 3-dimensional structure. When expressed in Escherichia coli, the three proteins differentially affected growth responses to salt, cold and heat stress. These results confirmed that the complete motif repeat structures and an appropriate hydrophilic/hydrophobic balance are important for the LEA protein in providing protection against various forms of stresses.

Ds transposon is biased towards providing splice donor sites for exonization in transgenic tobacco

Insertion of transposed elements into introns can lead to their activation as alternatively spliced cassette exons, an event called exonization, which can enrich the complexity of transcriptomes and proteomes. In this study, the first exonization event was detected when the modified rice EPSPS marker gene was inserted with the Ac transposon 5' end, which provided a splice donor site to yield abundant novel transcripts. To assess the contribution of splice donor and acceptor sites of transposon sequences, we inserted a Ds element into each intron of the EPSPS marker gene. This process yielded 14 constructs, with the Ds transposon inserted in the forward and reverse direction in each of the 7 introns of the EPSPS marker gene. The constructs were transformed into tobacco plants, and novel transcripts were identified by RT-PCR with specific primers. Exonization of Ds in EPSPS was biased towards providing splice donor sites of the inserted Ds sequence. Additionally, when the Ds inserted in reverse direction, a continuous splice donor consensus region was determined by offering 4 donor sites in the same intron. Information on these exonization events may help enhance gene divergence and functional genomic studies.

Alternative splicing: enhancing ability to cope with stress via transcriptome plasticity

Alternative splicing is a mechanism for the regulation of gene expression that is widespread in higher eukaryotes. Genome-wide approaches, based on comparison of expressed and genomic sequences, on tiling arrays, and on next-generation sequencing, have provided growing knowledge of the extent, distribution and association of alternative splicing with stress-related genes in plants. The functional meaning of alternative splicing in response to stress has been defined for many genes, and in particular for those involved in the regulation of the stress responses, such as protein kinases, transcription factors, splicing regulators and pathogen-resistance genes. The production of proteins with diverse domain rearrangements from the same gene is the main alternative splicing mechanism for pathogen-resistance genes. The plant response to abiotic stress is also characterized by a second mechanism, which consists of the expression of alternative transcripts that are targeted to nonsense-mediated decay. These quantitatively regulate stress-related gene expression. Many alternative splicing events are well conserved among plant species, and also across kingdoms, especially those observed in response to stress, for genes encoding splicing regulators, and other classes of RNA-binding proteins. Nevertheless, non-conserved events indicate that alternative splicing represents an evolutionary strategy that rapidly increases genome plasticity and develops new gene functions, along with other mechanisms such as gene duplication. Finally, the study of the naturally occurring variability of alternative splicing and the identification of genomic regions involved in the regulation of alternative splicing in crops are proposed as strategies for selecting genotypes with superior performance under adverse environmental conditions.

Deciphering the plant splicing code: experimental and computational approaches for predicting alternative splicing and splicing regulatory elements

Extensive alternative splicing (AS) of precursor mRNAs (pre-mRNAs) in multicellular eukaryotes increases the protein-coding capacity of a genome and allows novel ways to regulate gene expression. In flowering plants, up to 48% of intron-containing genes exhibit AS. However, the full extent of AS in plants is not yet known, as only a few high-throughput RNA-Seq studies have been performed. As the cost of obtaining RNA-Seq reads continues to fall, it is anticipated that huge amounts of plant sequence data will accumulate and help in obtaining a more complete picture of AS in plants. Although it is not an onerous task to obtain hundreds of millions of reads using high-throughput sequencing technologies, computational tools to accurately predict and visualize AS are still being developed and refined. This review will discuss the tools to predict and visualize transcriptome-wide AS in plants using short-reads and highlight their limitations. Comparative studies of AS events between plants and animals have revealed that there are major differences in the most prevalent types of AS events, suggesting that plants and animals differ in the way they recognize exons and introns. Extensive studies have been performed in animals to identify cis-elements involved in regulating AS, especially in exon skipping. However, few such studies have been carried out in plants. Here, we review the current state of research on splicing regulatory elements (SREs) and briefly discuss emerging experimental and computational tools to identify cis-elements involved in regulation of AS in plants. The availability of curated alternative splice forms in plants makes it possible to use computational tools to predict SREs involved in AS regulation, which can then be verified experimentally. Such studies will permit identification of plant-specific features involved in AS regulation and contribute to deciphering the splicing code in plants.

Predicting the impact of alternative splicing on plant MADS domain protein function

Several genome-wide studies demonstrated that alternative splicing (AS) significantly increases the transcriptome complexity in plants. However, the impact of AS on the functional diversity of proteins is difficult to assess using genome-wide approaches. The availability of detailed sequence annotations for specific genes and gene families allows for a more detailed assessment of the potential effect of AS on their function. One example is the plant MADS-domain transcription factor family, members of which interact to form protein complexes that function in transcription regulation. Here, we perform an in silico analysis of the potential impact of AS on the protein-protein interaction capabilities of MIKC-type MADS-domain proteins. We first confirmed the expression of transcript isoforms resulting from predicted AS events. Expressed transcript isoforms were considered functional if they were likely to be translated and if their corresponding AS events either had an effect on predicted dimerisation motifs or occurred in regions known to be involved in multimeric complex formation, or otherwise, if their effect was conserved in different species. Nine out of twelve MIKC MADS-box genes predicted to produce multiple protein isoforms harbored putative functional AS events according to those criteria. AS events with conserved effects were only found at the borders of or within the K-box domain. We illustrate how AS can contribute to the evolution of interaction networks through an example of selective inclusion of a recently evolved interaction motif in the MADS AFFECTING FLOWERING1-3 (MAF1-3) subclade. Furthermore, we demonstrate the potential effect of an AS event in SHORT VEGETATIVE PHASE (SVP), resulting in the deletion of a short sequence stretch including a predicted interaction motif, by overexpression of the fully spliced and the alternatively spliced SVP transcripts. For most of the AS events we were able to formulate hypotheses about the potential impact on the interaction capabilities of the encoded MIKC proteins.

Genome-wide survey of alternative splicing in the grass Brachypodium distachyon: a emerging model biosystem for plant functional genomics

A draft sequence of the genome of Brachypodium distachyon, the emerging grass model, was recently released. This represents a unique opportunity to determine its functional diversity compared to the genomes of other model species. Using homology mapping of assembled expressed sequence tags with chromosome scale pseudomolecules, we identified 128 alternative splicing events in B. distachyon. Our study identified that retention of introns is the major type of alternative splicing events (53%) in this plant and highlights the prevalence of splicing site recognition for definition of introns in plants. We have analyzed the compositional profiles of exon-intron junctions by base-pairing nucleotides with U1 snRNA which serves as a model for describing the possibility of sequence conservation. The alternative splicing isoforms identified in this study are novel and represent one of the potentially biologically significant means by which B. distachyon controls the function of its genes. Our observations serve as a basis to understand alternative splicing events of cereal crops with more complex genomes, like wheat or barley.

Evolution and expression of class III peroxidases

Class III peroxidases are members of a large multigenic family, only detected in the plant kingdom and absent from green algae sensu stricto (chlorophyte algae or Chlorophyta). Their evolution is thought to be related to the emergence of the land plants. However class III peroxidases are present in a lower copy number in some basal Streptophytes (Charapyceae), which predate land colonization. Gene structures are variable among organisms and within species with respect to the number of introns, but their positions are highly conserved. Their high copy number, as well as their conservation could be related to plant complexity and adaptation to increasing stresses. No specific function has been assigned to respective isoforms, but in large multigenic families, particular structure-function relations can be expected. Plant peroxidase sequences contain highly conserved residues and motifs, variable domains surrounded by conserved residues and present a low identity level among their promoter regions, further suggesting the existence of sub-functionalization of the different isoforms.

An interspecific plant hybrid shows novel changes in parental splice forms of genes for splicing factors

Interspecific hybridization plays an important role in plant adaptive evolution and speciation, and the process often results in phenotypic novelty. Hybrids can show changes in genome structure and gene expression compared with their parents including chromosomal rearrangments, changes in cytosine methylation, up- and downregulation of gene expression, and gene silencing. Alternative splicing (AS) is a fundamental aspect of the expression of many genes. However alternative splicing patterns have not been examined in multiple genes in an interspecific plant hybrid compared with its parents. Here we studied alternative splicing patterns in an interspecific Populus hybrid and its parents by assaying 40 genes using reverse transcription PCR. Most of the genes showed identical alternative splicing patterns between the parents and the hybrid. We found new alternative splicing variants present in the hybrid in two SR genes involved in the regulation of splicing and alternative splicing. The novel alternative splicing patterns included changes in donor and acceptor sites to create a new exon in one allele of PtRSZ22 in the hybrid and retention of an intron in both alleles of PtSR34a.1 in the hybrid, with effects on the function of the corresponding truncated proteins, if present. Our results suggest that novel alternative splicing patterns are present in a small percentage of genes in hybrids, but they could make a considerable impact on the expression of some genes. Changes in alternative splicing are likely to be an important component of the genetic changes that occur upon interspecific hybridization.

Carrot alternative oxidase gene AOX2a demonstrates allelic and genotypic polymorphisms in intron 3

Single nucleotide polymorphisms (SNPs) and insertion-deletions (InDels) are becoming important genetic markers for major crop species. In this study, we focus on variations at genomic level of the Daucus carota L. AOX2a gene. The use of gene-specific primers designed in exon regions on the boundaries of introns permitted to recognize intron length polymorphism (ILP) in intron 3 AOX2a by simple polymerase chain reaction (PCR) assays. The length of intron 3 can vary in individual carrot plants. Thus, allelic variation can be used as a tool to discriminate between single plant genotypes. Using this approach, individual plants from cv. Rotin and from diverse breeding lines and cultivars were identified that showed genetic variability by AOX2a ILPs. Repetitive patterns of intron length variation have been observed which allows grouping of genotypes. Polymorphic and identical PCR fragments revealed underlying high levels of sequence polymorphism. Variability was due to InDel events and intron single nucleotide polymorphisms (ISNPs), with a repetitive deletion in intron 3 affecting a putative pre-miRNA site. The results suggest that high AOX2a gene diversity in D. carota can be explored for the development of functional markers related to agronomic traits.

Intron polymorphism pattern in AOX1b of wild St John's wort (Hypericum perforatum) allows discrimination between individual plants

The present paper deals with the analysis of natural polymorphism in a selected alternative oxidase (AOX) gene of the medicinal plant, St John's wort. Four partial AOX gene sequences were isolated from the genomic DNA of a wild plant of Hypericum perforatum L. Three genes belong to the subfamily AOX1 (HpAOX1a, b and c) and one to the subfamily AOX2 (HpAOX2). The partial sequence of HpAOX1b showed polymerase chain reaction (PCR) fragment size variation as a result of variable lengths in two introns. PCR performed by Exon Primed Intron Crossing (EPIC)-PCR displayed the same two-band pattern in six plants from a collection. Both fragments showed identical sequences for all exons. However, each of the two introns showed an insertion/deletion (InDel) in identical positions for all plants that counted for the difference in the two fragment sizes. The InDel in intron 1 influenced the predictability of a pre-microRNA site. The almost identical PCR fragment pattern was characterized by a high variability in the sequences. The InDels in both introns were linked to repetitive intron single nucleotide polymorphisms (ISNP)s. The polymorphic pattern obtained by InDels and ISNPs from both fragments together was appropriate to discriminate between all individual plants. We suggest that AOX sequence polymorphism in H. perforatum can be used for studies on gene diversity and biodiversity. Further, we conclude that AOX sequence polymorphism of individual plants should be considered in biological studies on AOX activity to exclude the influence of genetic diversity. The identified polymorphic fragments are available to be explored in future experiments as a potential source for functional marker development related to the characterization of origins/accessions and agronomic traits such as plant growth, development and yield stability.

Cloning, functional expression and phylogenetic analysis of plant sterol 24C-methyltransferases involved in sitosterol biosynthesis

Sterol 24C-methyltransferases (SMTs) constitute a group of sequence-related proteins that catalyze the distinct patterns of 24-alkyl sterols that occur throughout nature. Two SMT cDNAs (SMT2-1 and SMT2-2) were cloned by homology based PCR methods from young leaves of Glycine max (soybean) and the corresponding enzymes were expressed functionally in Escherichia coli. The full-length cDNA for SMT2-1 and SMT2-2 have open reading frames of 1086 bp and 1092 bp, respectively, and encode proteins of 361 and 363 residues with a calculated molecular mass of 40.3 and 40.4 kDa, respectively. The substrate preference of the two isoforms was similar yet they differed from SMT1; kinetically SMT2-1 and SMT2-2 generated k(cat) values for the optimal substrate 24(28)methylene lophenol of 0.8 min(-1) and 1.34 min(-1), respectively, compared to the activity of SMT1 that generated a k(cat) for the optimal substrate cycloartenol of 0.6 min(-1). SMT2-2 was purified to homogeneity and the subunit organization shown to be tetrameric in similar fashion to other cloned SMTs. Analysis of the accumulated products catalyzed by the recombinant enzymes demonstrated that soybean SMT2-1 and SMT2-2 operate transalkylation activities analogous to the soybean plant SMT1. Metabolite analyses correlated with transcript profiling of the three SMT isoforms during soybean maturation clearly demonstrated that SMT isoform expression determines specific C24-methyl to C24-ethyl ratios to flowering whereas with seed development there is a disconnection such that the SMT transcript levels decrease against an increase in sterol content; generally SMT2-2 is expressed more than SMT2-1 or SMT1. These observations suggest that the genes that encode SMT1 and SMT2 in sitosterol biosynthesis may have undergone divergent evolution. In support of this proposition, the genomic organization for SMT1 of fungi and protozoa align very closely with one another and to those of the plant SMT2; both sets of SMTs lack introns. Unexpectedly, the SMT1 from Glycine max and other embryophytes of diverse origin possess disparate intron-exon characteristics that can be shown relates back to the algae. Our results suggest that the order of SMT1 appearing before SMT2 in phytosterol synthesis arose recently in plant evolution in response to duplication of a more primitive SMT gene likely to have been bifunctional and catalytically promiscuous.

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

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

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.