Plant serine/arginine-rich proteins: roles in precursor messenger RNA splicing, plant development, and stress responses

Chemical Modulation of Alternative Splicing for Molecular-Target Identification by Potential Genetic Control in Agrochemical Research

Alternative splicing (AS), the process of removing introns from pre-mRNA and the rearrangement of exons to produce several types of mature transcripts, is a remarkable step preceding protein synthesis. In particular, it has now been conclusively shown that up to ∼95% of genes are alternatively spliced to generate a complex and diverse proteome in eukaryotic organisms. Consequently, AS is one of the determinants of the functional repertoire of cells. Many studies have revealed that AS in plants can be regulated by cell type, developmental stage, environmental stress, and the circadian clock. Moreover, increasing amounts of evidence reveal that chemical compounds can affect various steps during splicing to induce major effects on plant physiology. Hence, the chemical modulation of AS can serve as a good strategy for molecular-target identification in attempts to potentially control plant genetics. However, the kind of mechanisms involved in the chemical modulation of AS that can be used in agrochemical research remain largely unknown. This review introduces recent studies describing the specific roles AS plays in plant adaptation to environmental stressors and in the regulation of development. We also discuss recent advances in small molecules that induce alterations of AS and the possibility of using this strategy in agrochemical-target identification, giving a new direction for potential genetic control in agrochemical research.

Comparative analysis of transcriptome in two wheat genotypes with contrasting levels of drought tolerance

Drought tolerance is a complex trait that is governed by multiple genes. The study presents differential transcriptome analysis between drought-tolerant (Triticum aestivum Cv. C306) and drought-sensitive (Triticum aestivum Cv. WL711) genotypes, using Affymetrix GeneChip® Wheat Genome Array. Both genotypes exhibited diverse global transcriptional responses under control and drought conditions. Pathway analysis suggested significant induction or repression of genes involved in secondary metabolism, nucleic acid synthesis, protein synthesis, and transport in C306, as compared to WL711. Significant up- and downregulation of transcripts for enzymes, hormone metabolism, and stress response pathways were observed in C306 under drought. The elevated expression of plasma membrane intrinsic protein 1 and downregulation of late embryogenesis abundant in the leaf tissues could play an important role in delayed wilting in C306. The other regulatory genes such as MT, FT, AP2, SKP1, ABA2, ARF6, WRKY6, AOS, and LOX2 are involved in defense response in C306 genotype. Additionally, transcripts with unknown functions were identified as differentially expressed, which could participate in drought responses.

Identification of splice variant of OsGBF1 in Oryza sativa ssp. indica genotypes under salinity stress

G-box-binding factors are plant transcription factors, involved in a wide range of biological processes including abiotic stress responses. In this study, we analyzed the expression of OsGBF1 during salt stress in two contrasting Oryza sativa spp. indica genotypes, Rasi and Tellahamsa. Two-day-old seedlings were exposed to NaCl stress under two different conditions. One set was exposed to 100 mM NaCl before transferring to 250 mM (induction stress), while another set was transferred directly to 250 mM (shock stress). During early induction stress, OsGBF1 was up-regulated in Rasi when compared to Tellahamsa. We cloned full-length OsGBF1 from these two genotypes, and analyzed the sequences. Our analysis indicated the presence of transcript variants, which are designated as OsGBF1a and 1b. OsGBF1b variant retained introns, which lead to the generation of premature termination codon. OsGBF1b transcript levels were not significantly different at 12-h of induction stress in Tellahamsa and Rasi. At 24-h of shock stress, OsGBF1b was up-regulated in both genotypes and the transcript was more in Rasi. Since, OsGBF1a and 1b are predicted to be splice variants, we examined expression pattern of OsSKIP, a splicing factor and component of the spliceosome. In induction stress, OsSKIP was up-regulated at 12- and 24-h in Rasi when compared to Tellahamsa. On the contrary, at 24-h shock stress, OsSKIP was down-regulated in Rasi when compared to Tellahamsa. It is possible that OsSKIP expression was increased in Rasi during induction stress for accurate splicing that could aid in tolerance. This is the first report on OsGBF1 splice variant and the variant could have specific functions linked to stress tolerance in rice.

Identification and evaluation of PCR reference genes for host and pathogen in sugarcane-Sporisorium scitamineum interaction system

BACKGROUND

Sugarcane (Saccharum L. plant) is an important crop for sugar and bio-energy production around the world. Among sugarcane diseases, smut caused by Sporisorium scitamineum is one of the major fungal diseases causing severe losses to the sugarcane industry. The use of PCR reference genes is essential to the normalization of data on gene expression involving the sugarcane-S. scitamineum interaction system; however, no report that addresses criteria in selecting these reference genes has been published to date.

RESULTS

In this study, 10 sugarcane genes and eight S. scitamineum genes were selected as candidate PCR reference genes in the sugarcane-S. scitamineum interaction system. The stability and reliability of these 18 candidate genes were analyzed in smut-resistant (NCo376) and -susceptible (YC71-374) genotypes using the statistical algorithms geNorm, NormFinder, BestKeeper, and deltaCt method. Subsequently, the relative expression levels of the sugarcane chitinase I-3 gene and S. scitamineum chorismate mutase gene were determined to validate the applicability of these sugarcane and S. scitamineum PCR reference genes, respectively. We finally found that the acyl-CoA dehydrogenase gene (ACAD), serine/arginine repetitive matrix protein 1 gene (SARMp1), or their combination (ACAD + SARMp1) could be utilized as the most suitable reference genes for normalization of sugarcane gene expression in sugarcane bud tissues after S. scitamineum infection. Similarly, the inosine 5'-monophosphate dehydrogenase gene (S10), the SEC65-signal recognition particle subunit gene (S11), or their combination (S10 + S11) were suitable for normalization of S. scitamineum gene expression in sugarcane bud tissues.

CONCLUSIONS

The PCR reference genes ACAD, SARMp1, S10, and S11 may be employed in gene transcriptional studies involving the sugarcane-S. scitamineum interaction system.

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

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

Thermopriming triggers splicing memory in Arabidopsis

Abiotic and biotic stresses limit crop productivity. Exposure to a non-lethal stress, referred to as priming, can allow plants to survive subsequent and otherwise lethal conditions; the priming effect persists even after a prolonged stress-free period. However, the molecular mechanisms underlying priming are not fully understood. Here, we investigated the molecular basis of heat-shock memory and the role of priming in Arabidopsis thaliana. Comprehensive analysis of transcriptome-wide changes in gene expression and alternative splicing in primed and non-primed plants revealed that alternative splicing functions as a novel component of heat-shock memory. We show that priming of plants with a non-lethal heat stress results in de-repression of splicing after a second exposure to heat stress. By contrast, non-primed plants showed significant repression of splicing. These observations link 'splicing memory' to the ability of plants to survive subsequent and otherwise lethal heat stress. This newly discovered priming-induced splicing memory may represent a general feature of heat-stress responses in plants and other organisms as many of the key components are conserved among eukaryotes. Furthermore, this finding could facilitate the development of novel approaches to improve plant survival under extreme heat stress.

Global profiling of alternative splicing landscape responsive to drought, heat and their combination in wheat (Triticum aestivum L.)

Plant can acquire tolerance to environmental stresses via transcriptome reprogramming at transcriptional and alternative splicing (AS) levels. However, how AS coordinates with transcriptional regulation to contribute to abiotic stresses responses is still ambiguous. In this study, we performed genome-wide analyses of AS responses to drought stress (DS), heat stress (HS) and their combination (HD) in wheat seedlings, and further compared them with transcriptional responses. In total, we found 200, 3576 and 4056 genes exhibiting significant AS pattern changes in response to DS, HS and HD, respectively, and combined drought and heat stress can induce specific AS compared with individual one. In addition, wheat homeologous genes exhibited differential AS responses under stress conditions that more AS events occurred on B subgenome than on A and D genomes. Comparison of genes regulated at AS and transcriptional levels showed that only 12% of DS-induced AS genes were subjected to transcriptional regulation, whereas the proportion increased to ~40% under HS and HD. Functional enrichment analysis revealed that abiotic stress-responsive pathways tended to be highly overrepresented among these overlapped genes under HS and HD. Thus, we proposed that transcriptional regulation may play a major role in response to DS, which coordinates with AS regulation to contribute to HS and HD tolerance in wheat.

Genome-wide association studies of doubled haploid exotic introgression lines for root system architecture traits in maize (Zea mays L.)

Root system architecture (RSA) is becoming recognized as important for water and nutrient acquisition in plants. This study focuses on finding single nucleotide polymorphisms (SNPs) associated with seedling RSA traits from 300 doubled haploid (DH) lines derived from crosses between Germplasm Enhancement of Maize (GEM) accessions and inbred lines PHB47 and PHZ51. These DH lines were genotyped using 62,077 SNP markers, while root and shoot phenotype data were collected from 14-day old seedlings. Genome-wide association studies (GWAS) were conducted using three models to offset false positives/negatives. Multiple SNPs associated with seedling root traits were detected, some of which were within or linked to gene models that showed expression in seedling roots. Significant trait associations involving the SNP S5_152926936 on Chromosome 5 were detected in all three models, particularly the trait network area. The SNP is within the gene model GRMZM2G021110, which is expressed in roots at seedling stage. SNPs that were significantly associated with seedling root traits, and closely linked to gene models that encode proteins associated with root development were also detected. This study shows that the GEM-DH panel may be a source of allelic diversity for genes controlling seedling root development.

RNA-seq analysis reveals alternative splicing under salt stress in cotton, Gossypium davidsonii

BACKGROUND

Numerous studies have focused on the regulation of gene expression in response to salt stress at the transcriptional level; however, little is known about this process at the post-transcriptional level.

RESULTS

Using a diploid D genome wild salinity-tolerant cotton species, Gossypium davidsonii, we analyzed alternative splicing (AS) of genes related to salt stress by comparing high-throughput transcriptomes from salt-treated and well-watered roots and leaves. A total of 14,172 AS events were identified involving 6798 genes, of which intron retention (35.73%) was the most frequent, being detected in 3492 genes. Under salt stress, 1287 and 1228 differential alternative splicing (DAS) events were identified in roots and leaves, respectively. These DAS genes were associated with specific functional pathways, such as "responses to stress", "metabolic process" and "RNA splicing", implying that AS represents an important pathway of gene regulation in response to salt stress. Several salt response genes, such as pyrroline-5-carboxylate synthase (P5CS), K⁺ channel outward (KCO1), plasma membrane intrinsic protein (PIP) and WRKY33 which were involved in osmotic balance, ion homeostasis, water transportation and transcriptional regulation, respectively, were identified with differential alternative splicing under salt stress. Moreover, we revealed that 13 genes encoding Ser/Arg-rich (SR) proteins related to AS regulation were differentially alternatively spliced under salt stress.

CONCLUSION

This study first provide a comprehensive view of AS in G. davidsonii, and highlight novel insights into the potential roles of AS in plant responses to salt stress.

Combining Single Nucleotide Polymorphism Genotyping Array with Bulked Segregant Analysis to Map a Gene Controlling Adult Plant Resistance to Stripe Rust in Wheat Line 03031-1-5 H62

Stripe rust, caused by Puccinia striiformis f. sp. tritici, is one of the most devastating diseases of wheat worldwide. Growing resistant cultivars is considered the best approach to manage this disease. In order to identify the resistance gene(s) in wheat line 03031-1-5 H62, which displayed high resistance to stripe rust at adult plant stage, a cross was made between 03031-1-5 H62 and susceptible cultivar Avocet S. The mapping population was tested with Chinese P. striiformis f. sp. tritici race CYR32 through artificial inoculation in a field in Yangling, Shaanxi Province and under natural infection in Tianshui, Gansu Province. The segregation ratios indicated that the resistance was conferred by a single dominant gene, temporarily designated as YrH62. A combination of bulked segregant analysis (BSA) with wheat 90K single nucleotide polymorphism (SNP) array was used to identify molecular markers linked to YrH62. A total of 376 polymorphic SNP loci identified from the BSA analysis were located on chromosome 1B, from which 35 kompetitive allele-specific PCR (KASP) markers selected together with 84 simple sequence repeat (SSR) markers on 1B were used to screen polymorphism and a chromosome region associated with rust resistance was identified. To saturate the chromosomal region covering the YrH62 locus, a 660K SNP array was used to identify more SNP markers. To develop tightly linked markers for marker-assisted selection of YrH62 in wheat breeding, 18 SNPs were converted into KASP markers. A final linkage map consisting of 15 KASP and 3 SSR markers was constructed with KASP markers AX-109352427 and AX-109862469 flanking the YrH62 locus in a 1.0 cM interval. YrH62 explained 63.8 and 69.3% of the phenotypic variation for disease severity and infection type, respectively. YrH62 was located near the centromeric region of chromosome 1BS based on the positions of the SSR markers in 1B deletion bins. Based on the origin, responses to P. striiformis f. sp. tritici races, and marker distances, YrH62 is likely different from the other reported stripe rust resistance genes/quantitative trait loci on 1B. The gene and tightly linked KASP markers will be useful for breeding wheat cultivars with resistance to stripe rust.

Abiotic Stresses Modulate Landscape of Poplar Transcriptome via Alternative Splicing, Differential Intron Retention, and Isoform Ratio Switching

Abiotic stresses affect plant physiology, development, growth, and alter pre-mRNA splicing. Western poplar is a model woody tree and a potential bioenergy feedstock. To investigate the extent of stress-regulated alternative splicing (AS), we conducted an in-depth survey of leaf, root, and stem xylem transcriptomes under drought, salt, or temperature stress. Analysis of approximately one billion of genome-aligned RNA-Seq reads from tissue- or stress-specific libraries revealed over fifteen millions of novel splice junctions. Transcript models supported by both RNA-Seq and single molecule isoform sequencing (Iso-Seq) data revealed a broad array of novel stress- and/or tissue-specific isoforms. Analysis of Iso-Seq data also resulted in the discovery of 15,087 novel transcribed regions of which 164 show AS. Our findings demonstrate that abiotic stresses profoundly perturb transcript isoform profiles and trigger widespread intron retention (IR) events. Stress treatments often increased or decreased retention of specific introns - a phenomenon described here as differential intron retention (DIR). Many differentially retained introns were regulated in a stress- and/or tissue-specific manner. A subset of transcripts harboring super stress-responsive DIR events showed persisting fluctuations in the degree of IR across all treatments and tissue types. To investigate coordinated dynamics of intron-containing transcripts in the study we quantified absolute copy number of isoforms of two conserved transcription factors (TFs) using Droplet Digital PCR. This case study suggests that stress treatments can be associated with coordinated switches in relative ratios between fully spliced and intron-retaining isoforms and may play a role in adjusting transcriptome to abiotic stresses.

Herboxidiene triggers splicing repression and abiotic stress responses in plants

Background

Constitutive and alternative splicing of pre-mRNAs from multiexonic genes controls the diversity of the proteome; these precisely regulated processes also fine-tune responses to cues related to growth, development, and stresses. Small-molecule inhibitors that perturb splicing provide invaluable tools for use as chemical probes to uncover the molecular underpinnings of splicing regulation and as potential anticancer compounds.

Results

Here, we show that herboxidiene (GEX1A) inhibits both constitutive and alternative splicing. Moreover, GEX1A activates genome-wide transcriptional patterns involved in abiotic stress responses in plants. GEX1A treatment -activated ABA-inducible promoters, and led to stomatal closure. Interestingly, GEX1A and pladienolide B (PB) elicited similar cellular changes, including alterations in the patterns of transcription and splicing, suggesting that these compounds might target the same spliceosome complex in plant cells.

Conclusions

Our study establishes GEX1A as a potent splicing inhibitor in plants that can be used to probe the assembly, dynamics, and molecular functions of the spliceosome and to study the interplay between splicing stress and abiotic stresses, as well as having potential biotechnological applications.

Electronic supplementary material

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

Depletion of Arabidopsis SC35 and SC35-like serine/arginine-rich proteins affects the transcription and splicing of a subset of genes

Serine/arginine-rich (SR) proteins are important splicing factors which play significant roles in spliceosome assembly and splicing regulation. However, little is known regarding their biological functions in plants. Here, we analyzed the phenotypes of mutants upon depleting different subfamilies of Arabidopsis SR proteins. We found that loss of the functions of SC35 and SC35-like (SCL) proteins cause pleiotropic changes in plant morphology and development, including serrated leaves, late flowering, shorter roots and abnormal silique phyllotaxy. Using RNA-seq, we found that SC35 and SCL proteins play roles in the pre-mRNA splicing. Motif analysis revealed that SC35 and SCL proteins preferentially bind to a specific RNA sequence containing the AGAAGA motif. In addition, the transcriptions of a subset of genes are affected by the deletion of SC35 and SCL proteins which interact with NRPB4, a specific subunit of RNA polymerase II. The splicing of FLOWERING LOCUS C (FLC) intron1 and transcription of FLC were significantly regulated by SC35 and SCL proteins to control Arabidopsis flowering. Therefore, our findings provide mechanistic insight into the functions of plant SC35 and SCL proteins in the regulation of splicing and transcription in a direct or indirect manner to maintain the proper expression of genes and development.

SR proteins were identified to be important splicing factors. This work generated mutants of different subfamilies of the classic Arabidopsis SR proteins. Genetic analysis revealed that loss of the function of SC35/SCL proteins influences the plant development. This study revealed SC35/SCL proteins regulate alternative splicing, preferentially bind a specific RNA motif, interact with NRPB4, and affect the transcription of a subset of genes. This study further revealed that SC35/SCL proteins control flowering by regulating the splicing and transcription of FLC. These results shed light on the functions of SR proteins in plants.

Light Regulation of Alternative Pre-mRNA Splicing in Plants

Alternative splicing (AS) is a major post-transcriptional mechanism to enhance the diversity of proteome in response to environmental signals. Among the numerous external signals perceived by plants, light is the most crucial one. Plants utilize complex photoreceptor signaling networks to sense different light conditions and adjust their growth and development accordingly. Although light-mediated gene expression has been widely investigated, little is known regarding the mechanism of light affecting AS to modulate mRNA at the post-transcriptional level. In this minireview, we summarize current progresses on how light affects AS, and how sensory photoreceptors and retrograde signaling pathways may coordinately regulate AS of pre-mRNAs. In addition, we also discuss the possibility that AS of the mRNAs encoding photoreceptors may be involved in feedback control of AS. We hypothesize that light regulation of the expression and activity of splicing factors would be a major mechanism of light-mediated AS. The combination of genetic study and high-throughput analyses of AS and splicing complexes in response to light is likely to further advance our understanding of the molecular mechanisms underlying light control of AS and plant development.

Functions of long intergenic non-coding (linc) RNAs in plants

Whole transcriptome analyses in many organisms have revealed that most transcribed RNAs do not encode proteins. These non-coding RNAs likely contribute to the regulation of gene expression during the development of multicellular organisms. In eukaryotes, the roles of small RNAs, one class of non-coding RNAs, in transcriptional and post-transcriptional regulation have been well characterized. However, the functions of a second class of non-coding RNAs, long intergenic noncoding (linc) RNAs, are relatively unknown, especially in plants. Recent advances in RNA-seq and tiling microarray technologies have revealed the presence of many lincRNAs across plant species. This review focuses on the functions of lincRNAs that have been recently reported in plants. One of the most well characterized functions of lincRNAs is to epigenetically regulate gene expression by recruiting proteins for chromosome modification to specific loci. Second, lincRNAs are known to inhibit the physical interaction between microRNAs (miRNAs) and their target mRNAs thus controling protein levels of the target mRNAs. Lastly, lincRNAs control alternative splicing by binding and sequestering the proteins required for alternative splicing.

Pre-mRNA splicing repression triggers abiotic stress signaling in plants

Alternative splicing (AS) of precursor RNAs enhances transcriptome plasticity and proteome diversity in response to diverse growth and stress cues. Recent work has shown that AS is pervasive across plant species, with more than 60% of intron-containing genes producing different isoforms. Mammalian cell-based assays have discovered various inhibitors of AS. Here, we show that the macrolide pladienolide B (PB) inhibits constitutive splicing and AS in plants. Also, our RNA sequencing (RNA-seq) data revealed that PB mimics abiotic stress signals including salt, drought and abscisic acid (ABA). PB activates the abiotic stress- and ABA-responsive reporters RD29A::LUC and MAPKKK18::uidA in Arabidopsis thaliana and mimics the effects of ABA on stomatal aperture. Genome-wide analysis of AS by RNA-seq revealed that PB perturbs the splicing machinery and leads to a striking increase in intron retention and a reduction in other forms of AS. Interestingly, PB treatment activates the ABA signaling pathway by inhibiting the splicing of clade A PP2C phosphatases while still maintaining to some extent the splicing of ABA-activated SnRK2 kinases. Taken together, our data establish PB as an inhibitor and modulator of splicing and a mimic of abiotic stress signals in plants. Thus, PB reveals the molecular underpinnings of the interplay between stress responses, ABA signaling and post-transcriptional regulation in plants.

High-Resolution Expression Map of the Arabidopsis Root Reveals Alternative Splicing and lincRNA Regulation

The extent to which alternative splicing and long intergenic noncoding RNAs (lincRNAs) contribute to the specialized functions of cells within an organ is poorly understood. We generated a comprehensive dataset of gene expression from individual cell types of the Arabidopsis root. Comparisons across cell types revealed that alternative splicing tends to remove parts of coding regions from a longer, major isoform, providing evidence for a progressive mechanism of splicing. Cell-type-specific intron retention suggested a possible origin for this common form of alternative splicing. Coordinated alternative splicing across developmental stages pointed to a role in regulating differentiation. Consistent with this hypothesis, distinct isoforms of a transcription factor were shown to control developmental transitions. lincRNAs were generally lowly expressed at the level of individual cell types, but co-expression clusters provided clues as to their function. Our results highlight insights gained from analysis of expression at the level of individual cell types.

Alternative splicing of basic chitinase gene PR3b in the low-nicotine mutants of Nicotiana tabacum L. cv. Burley 21

Two unlinked semi-dominant loci, A (NIC1) and B (NIC2), control nicotine and related alkaloid biosynthesis in Burley tobaccos. Mutations in either or both loci (nic1 and nic2) lead to low nicotine phenotypes with altered environmental stress responses. Here we show that the transcripts derived from the pathogenesis-related (PR) protein gene PR3b are alternatively spliced to a greater extent in the nic1 and nic2 mutants of Burley 21 tobacco and the nic1nic2 double mutant. The alternative splicing results in a deletion of 65 nucleotides and introduces a premature stop codon into the coding region of PR3b that leads to a significant reduction of PR3b specific chitinase activity. Assays of PR3b splicing in F2 individuals derived from crosses between nic1 and nic2 mutants and wild-type plants showed that the splicing phenotype is controlled by the NIC1 and NIC2 loci, even though NIC1 and NIC2 are unlinked loci. Moreover, the transcriptional analyses showed that the splicing patterns of PR3b in the low-nicotine mutants were differentially regulated by jasmonate (JA) and ethylene (ET). These data suggest that the NIC1 and NIC2 loci display differential roles in regulating the alternative splicing of PR3b in Burley 21. The findings in this study have provided valuable information for extending our understanding of the broader effects of the low-nicotine mutants of Burley 21 and the mechanism by which JA and ET signalling pathways post-transcriptionally regulate the activity of PR3b protein.

NMD mechanism and the functions of Upf proteins in plant

Nonsense-mediated decay (NMD) mechanism, also called mRNA surveillance, is a universal mRNA degradation pathway in eukaryotes. Hundreds of genes can be regulated by NMD whether in single-celled or higher organisms. There have been many studies on NMD and NMD factors (Upf proteins) with regard to their crucial roles in mRNA decay, especially in mammals and yeast. However, research focusing on NMD in plant is still lacking compared to the research that has been dedicated to NMD in mammals and yeast. Even so, recent study has shown that NMD factors in Arabidopsis can provide resistance against biotic and abiotic stresses. This discovery and its associated developments have given plant NMD mechanism a new outlook and since then, more and more research has focused on this area. In this review, we focused mainly on the distinctive NMD micromechanism and functions of Upf proteins in plant with references to the role of mRNA surveillance in mammals and yeast. We also highlighted recent insights into the roles of premature termination codon location, trans-elements and functions of other NMD factors to emphasize the particularity of plant NMD. Furthermore, we also discussed conventional approaches and neoteric methods used in plant NMD researches.

Insect herbivory elicits genome-wide alternative splicing responses in Nicotiana attenuata

Changes in gene expression and alternative splicing (AS) are involved in many responses to abiotic and biotic stresses in eukaryotic organisms. In response to attack and oviposition by insect herbivores, plants elicit rapid changes in gene expression which are essential for the activation of plant defenses; however, the herbivory-induced changes in AS remain unstudied. Using mRNA sequencing, we performed a genome-wide analysis on tobacco hornworm (Manduca sexta) feeding-induced AS in both leaves and roots of Nicotiana attenuata. Feeding by M. sexta for 5 h reduced total AS events by 7.3% in leaves but increased them in roots by 8.0% and significantly changed AS patterns in leaves and roots of existing AS genes. Feeding by M. sexta also resulted in increased (in roots) and decreased (in leaves) transcript levels of the serine/arginine-rich (SR) proteins that are involved in the AS machinery of plants and induced changes in SR gene expression that were jasmonic acid (JA)-independent in leaves but JA-dependent in roots. Changes in AS and gene expression elicited by M. sexta feeding were regulated independently in both tissues. This study provides genome-wide evidence that insect herbivory induces changes not only in the levels of gene expression but also in their splicing, which might contribute to defense against and/or tolerance of herbivory.

Low temperature conditioning of garlic (Allium sativum L.) "seed" cloves induces alterations in sprouts proteome

Low-temperature conditioning of garlic "seed" cloves substitutes the initial climatic requirements of the crop and accelerates the cycle. We have reported that "seed" bulbs from "Coreano" variety conditioned at 5°C for 5 weeks reduces growth and plant weight as well as the crop yields and increases the synthesis of phenolic compounds and anthocyanins. Therefore, this treatment suggests a cold stress. Plant acclimation to stress is associated with deep changes in proteome composition. Since proteins are directly involved in plant stress response, proteomics studies can significantly contribute to unravel the possible relationships between protein abundance and plant stress acclimation. The aim of this work was to study the changes in the protein profiles of garlic "seed" cloves subjected to conditioning at low-temperature using proteomics approach. Two sets of garlic bulbs were used, one set was stored at room temperature (23°C), and the other was conditioned at low temperature (5°C) for 5 weeks. Total soluble proteins were extracted from sprouts of cloves and separated by two-dimensional gel electrophoresis. Protein spots showing statistically significant changes in abundance were analyzed by LC-ESI-MS/MS and identified by database search analysis using the Mascot search engine. The results revealed that low-temperature conditioning of garlic "seed" cloves causes alterations in the accumulation of proteins involved in different physiological processes such as cellular growth, antioxidative/oxidative state, macromolecules transport, protein folding and transcription regulation process. The metabolic pathways affected include protein biosynthesis and quality control system, photosynthesis, photorespiration, energy production, and carbohydrate and nucleotide metabolism. These processes can work cooperatively to establish a new cellular homeostasis that might be related with the physiological and biochemical changes observed in previous studies.

MAF2 Is Regulated by Temperature-Dependent Splicing and Represses Flowering at Low Temperatures in Parallel with FLM

Plants enter their reproductive phase when the environmental conditions are favourable for the successful production of progeny. The transition from vegetative to reproductive phase is influenced by several environmental factors including ambient temperature. In the model plant Arabidopsis thaliana, SHORT VEGETATIVE PHASE (SVP) is critical for this pathway; svp mutants cannot modify their flowering time in response to ambient temperature. SVP encodes a MADS-box transcription factor that directly represses genes that promote flowering. SVP binds DNA in complexes with other MADS-box transcription factors, including FLOWERING LOCUS M (FLM), which acts with SVP to repress the floral transition at low temperatures. Small temperature changes post-transcriptionally regulate FLM through temperature-dependent alternative splicing (TD-AS). As ambient temperature increases, the predominant FLM splice isoform shifts to encode a protein incapable of exerting a repressive effect on flowering. Here we characterize a closely related MADS-box transcription factor, MADS AFFECTING FLOWERING2 (MAF2), which has independently evolved TD-AS. At low temperatures the most abundant MAF2 splice variant encodes a protein that interacts with SVP to repress flowering. At increased temperature the relative abundance of splice isoforms shifts in favour of an intron-retaining variant that introduces a premature termination codon. We show that this isoform encodes a protein that cannot interact with SVP or repress flowering. At lower temperatures MAF2 and SVP repress flowering in parallel with FLM and SVP, providing an additional input to sense ambient temperature for the control of flowering.

A serine-arginine-rich (SR) splicing factor modulates alternative splicing of over a thousand genes in Toxoplasma gondii

Single genes are often subject to alternative splicing, which generates alternative mature mRNAs. This phenomenon is widespread in animals, and observed in over 90% of human genes. Recent data suggest it may also be common in Apicomplexa. These parasites have small genomes, and economy of DNA is evolutionarily favoured in this phylum. We investigated the mechanism of alternative splicing in Toxoplasma gondii, and have identified and localized TgSR3, a homologue of ASF/SF2 (alternative-splicing factor/splicing factor 2, a serine-arginine–rich, or SR protein) to a subnuclear compartment. In addition, we conditionally overexpressed this protein, which was deleterious to growth. qRT-PCR was used to confirm perturbation of splicing in a known alternatively-spliced gene. We performed high-throughput RNA-seq to determine the extent of splicing modulated by this protein. Current RNA-seq algorithms are poorly suited to compact parasite genomes, and hence we complemented existing tools by writing a new program, GeneGuillotine, that addresses this deficiency by segregating overlapping reads into distinct genes. In order to identify the extent of alternative splicing, we released another program, JunctionJuror, that detects changes in intron junctions. Using this program, we identified about 2000 genes that were constitutively alternatively spliced in T. gondii. Overexpressing the splice regulator TgSR3 perturbed alternative splicing in over 1000 genes.

Alternative splicing in plants: directing traffic at the crossroads of adaptation and environmental stress

In recent years, high-throughput sequencing-based analysis of plant transcriptomes has suggested that up to ∼60% of plant gene loci encode alternatively spliced mature transcripts. These studies have also revealed that alternative splicing in plants can be regulated by cell type, developmental stage, the environment, and the circadian clock. Alternative splicing is coupled to RNA surveillance and processing mechanisms, including nonsense mediated decay. Recently, non-protein-coding transcripts have also been shown to undergo alternative splicing. These discoveries collectively describe a robust system of post-transcriptional regulatory feedback loops which influence RNA abundance. In this review, we summarize recent studies describing the specific roles alternative splicing and RNA surveillance play in plant adaptation to environmental stresses and the regulation of the circadian clock.

Differential recruitment of splice variants from SR pre-mRNAs to polysomes during development and in response to stresses

We have previously shown that precursor mRNAs (pre-mRNAs) of serine/arginine-rich (SR) proteins are extensively alternatively spliced to generate approximately 100 distinct splice variants from 14 SR genes and that the splicing pattern of SR pre-mRNAs changes in different organs and in response to abiotic stresses. About half of the splice variants are potential targets of nonsense-mediated decay (NMD) and 25 splice forms were confirmed to be real NMD targets. However, it is not known whether (i) all splice variants are recruited to polysomes for translation; (ii) there is a preferential recruitment of specific splice isoforms to polysomes; and (iii) there is a differential recruitment of splice variants during development and in response to stresses. To address these questions, we analyzed the association of SR splice variants with polysomes from seedlings, different organs and seedlings exposed to heat and cold stress. In seedlings, about one-third of the splice variants (22 out of 72) are not recruited to polysomes. Among those associated with polysomes, the functional isoforms that code for full-length proteins and some candidate putative and confirmed NMD targets were identified. There was preferential recruitment of some splice forms over others. Predominant recruitment of functional isoforms along with a few NMD candidates was found in different organs. Furthermore, we observed differential recruitment of isoforms in different organs. Heat and cold stress enhanced or reduced recruitment of specific splice variants. Our studies reveal differential recruitment of SR splice variants to polysomes under normal conditions, during development and in response to stresses.

Environmental stresses modulate abundance and timing of alternatively spliced circadian transcripts in Arabidopsis

Environmental stresses profoundly altered accumulation of nonsense mRNAs including intron-retaining (IR) transcripts in Arabidopsis. Temporal patterns of stress-induced IR mRNAs were dissected using both oscillating and non-oscillating transcripts. Broad-range thermal cycles triggered a sharp increase in the long IR CCA1 isoforms and altered their phasing to different times of day. Both abiotic and biotic stresses such as drought or Pseudomonas syringae infection induced a similar increase. Thermal stress induced a time delay in accumulation of CCA1 I4Rb transcripts, whereas functional mRNA showed steady oscillations. Our data favor a hypothesis that stress-induced instabilities of the central oscillator can be in part compensated through fluctuations in abundance and out-of-phase oscillations of CCA1 IR transcripts. Taken together, our results support a concept that mRNA abundance can be modulated through altering ratios between functional and nonsense/IR transcripts. SR45 protein specifically bound to the retained CCA1 intron in vitro, suggesting that this splicing factor could be involved in regulation of intron retention. Transcriptomes of nonsense-mediated mRNA decay (NMD)-impaired and heat-stressed plants shared a set of retained introns associated with stress- and defense-inducible transcripts. Constitutive activation of certain stress response networks in an NMD mutant could be linked to disequilibrium between functional and nonsense mRNAs.

Genome-wide analysis of alternative splicing of pre-mRNA under salt stress in Arabidopsis

Background

Alternative splicing (AS) of precursor mRNA (pre-mRNA) is an important gene regulation process that potentially regulates many physiological processes in plants, including the response to abiotic stresses such as salt stress.

Results

To analyze global changes in AS under salt stress, we obtained high-coverage (~200 times) RNA sequencing data from Arabidopsis thaliana seedlings that were treated with different concentrations of NaCl. We detected that ~49% of all intron-containing genes were alternatively spliced under salt stress, 10% of which experienced significant differential alternative splicing (DAS). Furthermore, AS increased significantly under salt stress compared with under unstressed conditions. We demonstrated that most DAS genes were not differentially regulated by salt stress, suggesting that AS may represent an independent layer of gene regulation in response to stress. Our analysis of functional categories suggested that DAS genes were associated with specific functional pathways, such as the pathways for the responses to stresses and RNA splicing. We revealed that serine/arginine-rich (SR) splicing factors were frequently and specifically regulated in AS under salt stresses, suggesting a complex loop in AS regulation for stress adaptation. We also showed that alternative splicing site selection (SS) occurred most frequently at 4 nucleotides upstream or downstream of the dominant sites and that exon skipping tended to link with alternative SS.

Conclusions

Our study provided a comprehensive view of AS under salt stress and revealed novel insights into the potential roles of AS in plant response to salt stress.

Electronic supplementary material

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

Genome-wide analysis of light-regulated alternative splicing mediated by photoreceptors in Physcomitrella patens

Background

Light is one of the most important factors regulating plant growth and development. Light-sensing photoreceptors tightly regulate gene expression to control photomorphogenic responses. Although many levels of gene expression are modulated by photoreceptors, regulation at the mRNA splicing step remains unclear.

Results

We performed high-throughput mRNA sequencing to analyze light-responsive changes in alternative splicing in the moss Physcomitrella patens, and found that a large number of alternative splicing events were induced by light in the moss protonema. Light-responsive intron retention preferentially occurred in transcripts involved in photosynthesis and translation. Many of the alternatively spliced transcripts were expressed from genes with a function relating to splicing or light signaling, suggesting a potential impact on pre-mRNA splicing and photomorphogenic gene regulation in response to light. Moreover, most light-regulated intron retention was induced immediately upon light exposure, while motif analysis identified a repetitive GAA motif that may function as an exonic regulatory cis element in light-mediated alternative splicing. Further analysis in gene-disrupted mutants was consistent with a function for multiple red-light photoreceptors in the upstream regulation of light-responsive alternative splicing.

Conclusions

Our results indicate that intensive alternative splicing occurs in non-vascular plants and that, during photomorphogenesis, light regulates alternative splicing with transcript selectivity. We further suggest that alternative splicing is rapidly fine-tuned by light to modulate gene expression and reorganize metabolic processes, and that pre-mRNA cis elements are involved in photoreceptor-mediated splicing regulation.

Diversity and evolution of spliceosomal systems

The intron-exon structures of eukaryotic nuclear genomes exhibit tremendous diversity across different species. The availability of many genomes from diverse eukaryotic species now allows for the reconstruction of the evolutionary history of this diversity. Consideration of spliceosomal systems in comparative context reveals a surprising and very complex portrait: in contrast to many expectations, gene structures in early eukaryotic ancestors were highly complex and "animal or plant-like" in many of their spliceosomal structures has occurred; pronounced simplification of gene structures, splicing signals, and spliceosomal machinery occurring independently in many lineages. In addition, next-generation sequencing of transcripts has revealed that alternative splicing is more common across eukaryotes than previously thought. However, much alternative splicing in diverse eukaryotes appears to play a regulatory role: alternative splicing fulfilling the most famous role for alternative splicing-production of multiple different proteins from a single gene-appears to be much more common in animal species than in nearly any other lineage.

Coexpression network analysis associated with call of rice seedlings for encountering heat stress

Coexpression network analysis is useful tool for identification of functional association of coexpressed genes. We developed a coexpression network of rice from heat stress transcriptome data. Global transcriptome of rice leaf tissues was performed by microarray at three time points--post 10 and 60 min heat stress at 42 °C and 30 min recovery at 26 °C following 60 min 42 °C heat stress to investigate specifically the early events in the heat stress and recovery response. The transcriptome profile was significantly modulated within 10 min of heat stress. Strikingly, the number of up-regulated genes was higher than the number of down-regulated genes in 10 min of heat stress. The enrichment of GO terms protein kinase activity/protein serine threonine kinase activity, response to heat and reactive oxygen species in up-regulated genes after 10 min signifies the role of signal transduction events and reactive oxygen species during early heat stress. The enrichment of transcription factor (TF) binding sites for heat shock factors, bZIPs and DREBs coupled with up-regulation of TFs of different families suggests that the heat stress response in rice involves integration of various regulatory networks. The interpretation of microarray data in the context of coexpression network analysis identified several functionally correlated genes consisting of previously documented heat upregulated genes as well as new genes that can be implicated in heat stress. Based on the findings on parallel analysis of growth of seedlings, associated changes in transcripts of selected Hsps, genome-wide microarray profiling and the coexpression network analysis, this study is a step forward in understanding heat response of rice, the world's most important food crop.

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.

The regulatory role of Pcf11-similar-4 (PCFS4) in Arabidopsis development by genome-wide physical interactions with target loci

Background

The yeast and human Pcf11 functions in both constitutive and regulated transcription and pre-mRNA processing. The constitutive roles of PCF11 are largely mediated by its direct interaction with RNA Polymerase II C-terminal domain and a polyadenylation factor, Clp1. However, little is known about the mechanism of the regulatory roles of Pcf11. Though similar to Pcf11 in multiple aspects, Arabidopsis Pcf11-similar-4 protein (PCFS4) plays only a regulatory role in Arabidopsis gene expression. Towards understanding how PCFS4 regulates the expression of its direct target genes in a genome level, ChIP-Seq approach was employed in this study to identify PCFS4 enrichment sites (ES) and the ES-linked genes within the Arabidopsis genome.

Results

A total of 892 PCFS4 ES sites linked to 839 genes were identified. Distribution analysis of the ES sites along the gene bodies suggested that PCFS4 is preferentially located on the coding sequences of the genes, consistent with its regulatory role in transcription and pre-mRNA processing. Gene ontology (GO) analysis revealed that the ES-linked genes were specifically enriched in a few GO terms, including those categories of known PCFS4 functions in Arabidopsis development. More interestingly, GO analysis suggested novel roles of PCFS4. An example is its role in circadian rhythm, which was experimentally verified herein. ES site sequences analysis identified some over-represented sequence motifs shared by subsets of ES sites. The motifs may explain the specificity of PCFS4 on its target genes and the PCFS4's functions in multiple aspects of Arabidopsis development and behavior.

Conclusions

Arabidopsis PCFS4 has been shown to specifically target on, and physically interact with, the subsets of genes. Its targeting specificity is likely mediated by cis-elements shared by the genes of each subset. The potential regulation on both transcription and mRNA processing levels of each subset of the genes may explain the functions of PCFS4 in multiple aspects of Arabidopsis development and behavior.

The role of canonical and noncanonical pre-mRNA splicing in plant stress responses

Plants are sessile organisms capable of adapting to various environmental constraints, such as high or low temperatures, drought, soil salinity, or pathogen attack. To survive the unfavorable conditions, plants actively employ pre-mRNA splicing as a mechanism to regulate expression of stress-responsive genes and reprogram intracellular regulatory networks. There is a growing evidence that various stresses strongly affect the frequency and diversity of alternative splicing events in the stress-responsive genes and lead to an increased accumulation of mRNAs containing premature stop codons, which in turn have an impact on plant stress response. A number of studies revealed that some mRNAs involved in plant stress response are spliced counter to the traditional conception of alternative splicing. Such noncanonical mRNA splicing events include trans-splicing, intraexonic deletions, or variations affecting multiple exons and often require short direct repeats to occur. The noncanonical alternative splicing, along with common splicing events, targets the spliced transcripts to degradation through nonsense-mediated mRNA decay or leads to translation of truncated proteins. Investigation of the diversity, biological consequences, and mechanisms of the canonical and noncanonical alternative splicing events will help one to identify those transcripts which are promising for using in genetic engineering and selection of stress-tolerant plants.

Identification of an intronic splicing regulatory element involved in auto-regulation of alternative splicing of SCL33 pre-mRNA

In Arabidopsis, pre-mRNAs of serine/arginine-rich (SR) proteins undergo extensive alternative splicing (AS). However, little is known about the cis-elements and trans-acting proteins involved in regulating AS. Using a splicing reporter (GFP-intron-GFP), consisting of the GFP coding sequence interrupted by an alternatively spliced intron of SCL33, we investigated whether cis-elements within this intron are sufficient for AS, and which SR proteins are necessary for regulated AS. Expression of the splicing reporter in protoplasts faithfully produced all splice variants from the intron, suggesting that cis-elements required for AS reside within the intron. To determine which SR proteins are responsible for AS, the splicing pattern of the GFP-intron-GFP reporter was investigated in protoplasts of three single and three double mutants of SR genes. These analyses revealed that SCL33 and a closely related paralog, SCL30a, are functionally redundant in generating specific splice variants from this intron. Furthermore, SCL33 protein bound to a conserved sequence in this intron, indicating auto-regulation of AS. Mutations in four GAAG repeats within the conserved region impaired generation of the same splice variants that are affected in the scl33 scl30a double mutant. In conclusion, we have identified the first intronic cis-element involved in AS of a plant SR gene, and elucidated a mechanism for auto-regulation of AS of this intron.

Interactions of SR45, an SR-like protein, with spliceosomal proteins and an intronic sequence: insights into regulated splicing

SR45 is a serine/arginine-rich (SR)-like protein with two arginine/serine-rich (RS) domains. We have previously shown that SR45 regulates alternative splicing (AS) by differential selection of 5' and 3' splice sites. However, it is unknown how SR45 regulates AS. To gain mechanistic insights into the roles of SR45 in splicing, we screened a yeast two-hybrid library with SR45. This screening resulted in the isolation of two spliceosomal proteins, U1-70K and U2AF(35) b that are known to function in 5' and 3' splice site selection, respectively. This screen not only confirmed our prior observation that U1-70K and SR45 interact, but also helped to identify an additional interacting partner (U2AF(35) ). In vitro and in vivo analyses revealed an interaction of SR45 with both paralogs of U2AF(35) . Furthermore, we show that the RS1 and RS2 domains of SR45, and not the RNA recognition motif (RRM) domain, associate independently with both U2AF(35) proteins. Interaction studies among U2AF(35) paralogs and between U2AF(35) and U1-70K revealed that U2AF(35) can form homo- or heterodimers and that U2AF(35) proteins can associate with U1-70K. Using RNA probes from SR30 intron 10, whose splicing is altered in the sr45 mutant, we show that SR45 and U2AF(35) b bind to different parts of the intron, with a binding site for SR45 in the 5' region and two binding regions, each ending with a known 3' splice site, for U2AF(35) b. These results suggest that SR45 recruits U1snRNP and U2AF to 5' and 3' splice sites, respectively, by interacting with pre-mRNA, U1-70K and U2AF(35) and modulates AS.

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.