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

Transcriptome sequencing and comparative analysis of cucumber flowers with different sex types

Background

Cucumber, Cucumis sativus L., is an economically and nutritionally important crop of the Cucurbitaceae family and has long served as a primary model system for sex determination studies. Recently, the sequencing of its whole genome has been completed. However, transcriptome information of this species is still scarce, with a total of around 8,000 Expressed Sequence Tag (EST) and mRNA sequences currently available in GenBank. In order to gain more insights into molecular mechanisms of plant sex determination and provide the community a functional genomics resource that will facilitate cucurbit research and breeding, we performed transcriptome sequencing of cucumber flower buds of two near-isogenic lines, WI1983G, a gynoecious plant which bears only pistillate flowers, and WI1983H, a hermaphroditic plant which bears only bisexual flowers.

Result

Using Roche-454 massive parallel pyrosequencing technology, we generated a total of 353,941 high quality EST sequences with an average length of 175bp, among which 188,255 were from gynoecious flowers and 165,686 from hermaphroditic flowers. These EST sequences, together with ~5,600 high quality cucumber EST and mRNA sequences available in GenBank, were clustered and assembled into 81,401 unigenes, of which 28,452 were contigs and 52,949 were singletons. The unigenes and ESTs were further mapped to the cucumber genome and more than 500 alternative splicing events were identified in 443 cucumber genes. The unigenes were further functionally annotated by comparing their sequences to different protein and functional domain databases and assigned with Gene Ontology (GO) terms. A biochemical pathway database containing 343 predicted pathways was also created based on the annotations of the unigenes. Digital expression analysis identified ~200 differentially expressed genes between flowers of WI1983G and WI1983H and provided novel insights into molecular mechanisms of plant sex determination process. Furthermore, a set of SSR motifs and high confidence SNPs between WI1983G and WI1983H were identified from the ESTs, which provided the material basis for future genetic linkage and QTL analysis.

Conclusion

A large set of EST sequences were generated from cucumber flower buds of two different sex types. Differentially expressed genes between these two different sex-type flowers, as well as putative SSR and SNP markers, were identified. These EST sequences provide valuable information to further understand molecular mechanisms of plant sex determination process and forms a rich resource for future functional genomics analysis, marker development and cucumber breeding.

Extensive divergence in alternative splicing patterns after gene and genome duplication during the evolutionary history of Arabidopsis

Gene duplication at various scales, from single gene duplication to whole-genome (WG) duplication, has occurred throughout eukaryotic evolution and contributed greatly to the large number of duplicated genes in the genomes of many eukaryotes. Previous studies have shown divergence in expression patterns of many duplicated genes at various evolutionary time scales and cases of gain of a new function or expression pattern by one duplicate or partitioning of functions or expression patterns between duplicates. Alternative splicing (AS) is a fundamental aspect of the expression of many genes that can increase gene product diversity and affect gene regulation. However, the evolution of AS patterns of genes duplicated by polyploidy, as well as in a sizable number of duplicated gene pairs in plants, has not been examined. Here, we have characterized conservation and divergence in AS patterns in genes duplicated by a polyploidy event during the evolutionary history of Arabidopsis thaliana. We used reverse transcription-polymerase chain reaction to assay 104 WG duplicates in six organ types and in plants grown under three abiotic stress treatments to detect organ- and stress-specific patterns of AS. Differences in splicing patterns in one or more organs, or under stress conditions, were found between the genes in a large majority of the duplicated pairs. In a few cases, AS patterns were the same between duplicates only under one or more abiotic stress treatments and not under normal growing conditions or vice versa. We also examined AS in 42 tandem duplicates and we found patterns of AS roughly comparable with the genes duplicated by polyploidy. The alternatively spliced forms in some of the genes created premature stop codons that would result in missing or partial functional domains if the transcripts are translated, which could affect gene function and cause functional divergence between duplicates. Our results indicate that AS patterns have diverged considerably after gene and genome duplication during the evolutionary history of the Arabidopsis lineage, sometimes in an organ- or stress-specific manner. AS divergence between duplicated genes may have contributed to gene functional evolution and led to preservation of some duplicated genes.

Genome-wide analysis of alternative splicing in Chlamydomonas reinhardtii

Background

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

Results

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

Conclusions

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

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.

Genome-wide cloning and sequence analysis of leucine-rich repeat receptor-like protein kinase genes in Arabidopsis thaliana

BACKGROUND

Transmembrane receptor kinases play critical roles in both animal and plant signaling pathways regulating growth, development, differentiation, cell death, and pathogenic defense responses. In Arabidopsis thaliana, there are at least 223 Leucine-rich repeat receptor-like kinases (LRR-RLKs), representing one of the largest protein families. Although functional roles for a handful of LRR-RLKs have been revealed, the functions of the majority of members in this protein family have not been elucidated.

RESULTS

As a resource for the in-depth analysis of this important protein family, the complementary DNA sequences (cDNAs) of 194 LRR-RLKs were cloned into the Gateway donor vector pDONR/Zeo and analyzed by DNA sequencing. Among them, 157 clones showed sequences identical to the predictions in the Arabidopsis sequence resource, TAIR8. The other 37 cDNAs showed gene structures distinct from the predictions of TAIR8, which was mainly caused by alternative splicing of pre-mRNA. Most of the genes have been further cloned into Gateway destination vectors with GFP or FLAG epitope tags and have been transformed into Arabidopsis for in planta functional analysis. All clones from this study have been submitted to the Arabidopsis Biological Resource Center (ABRC) at Ohio State University for full accessibility by the Arabidopsis research community.

CONCLUSIONS

Most of the Arabidopsis LRR-RLK genes have been isolated and the sequence analysis showed a number of alternatively spliced variants. The generated resources, including cDNA entry clones, expression constructs and transgenic plants, will facilitate further functional analysis of the members of this important gene family.

Genome-wide cloning and sequence analysis of leucine-rich repeat receptor-like protein kinase genes in Arabidopsis thaliana

Background

Transmembrane receptor kinases play critical roles in both animal and plant signaling pathways regulating growth, development, differentiation, cell death, and pathogenic defense responses. In Arabidopsis thaliana, there are at least 223 Leucine-rich repeat receptor-like kinases (LRR-RLKs), representing one of the largest protein families. Although functional roles for a handful of LRR-RLKs have been revealed, the functions of the majority of members in this protein family have not been elucidated.

Results

As a resource for the in-depth analysis of this important protein family, the complementary DNA sequences (cDNAs) of 194 LRR-RLKs were cloned into the Gateway^R donor vector pDONR/Zeo^R and analyzed by DNA sequencing. Among them, 157 clones showed sequences identical to the predictions in the Arabidopsis sequence resource, TAIR8. The other 37 cDNAs showed gene structures distinct from the predictions of TAIR8, which was mainly caused by alternative splicing of pre-mRNA. Most of the genes have been further cloned into Gateway^R destination vectors with GFP or FLAG epitope tags and have been transformed into Arabidopsis for in planta functional analysis. All clones from this study have been submitted to the Arabidopsis Biological Resource Center (ABRC) at Ohio State University for full accessibility by the Arabidopsis research community.

Conclusions

Most of the Arabidopsis LRR-RLK genes have been isolated and the sequence analysis showed a number of alternatively spliced variants. The generated resources, including cDNA entry clones, expression constructs and transgenic plants, will facilitate further functional analysis of the members of this important gene family.

Gene regulation during cold stress acclimation in plants

Cold stress adversely affects plant growth and development and thus limits crop productivity. Diverse plant species tolerate cold stress to a varying degree, which depends on reprogramming gene expression to modify their physiology, metabolism, and growth. Cold signal in plants is transmitted to activate CBF-dependent (C-repeat/drought-responsive element binding factor-dependent) and CBF-independent transcriptional pathway, of which CBF-dependent pathway activates CBF regulon. CBF transcription factor genes are induced by the constitutively expressed ICE1 (inducer of CBF expression 1) by binding to the CBF promoter. ICE1-CBF cold response pathway is conserved in diverse plant species. Transgenic analysis in different plant species revealed that cold tolerance can be significantly enhanced by genetic engineering CBF pathway. Posttranscriptional regulation at pre-mRNA processing and export from nucleus plays a role in cold acclimation. Small noncoding RNAs, namely micro-RNAs (miRNAs) and small interfering RNAs (siRNAs), are emerging as key players of posttranscriptional gene silencing. Cold stress-regulated miRNAs have been identified in Arabidopsis and rice. In this chapter, recent advances on cold stress signaling and tolerance are highlighted.

Extensive coupling of alternative splicing of pre-mRNAs of serine/arginine (SR) genes with nonsense-mediated decay

In Arabidopsis, pre-mRNAs encoding serine/arginine (SR) proteins, key regulators of constitutive and alternative splicing, are extensively alternatively spliced. In seedlings, 13 SR genes are alternatively spliced to generate 75 transcripts, of which 53 contain a premature termination codon (PTC). However, it is not known if any of the PTC-containing splice variants are the targets of nonsense-mediated decay (NMD) and if there is any link between NMD and the abundance of functional transcripts. Here, we analyzed the abundances of all splice variants for each alternatively spliced gene in an Arabidopsis mutant that lacks UPF3, one of the core components of NMD machinery, to determine if the PTC-containing transcripts are degraded by NMD. Our results show that about half of the 53 splice variants with a PTC are the targets of degradation by NMD. The accumulation of PTC-containing transcripts resulted in concomitant reduction in the amount of functional transcript. These results show widespread coupling of alternative splicing with NMD in the SR gene family, suggesting a strong link between unproductive splicing and the abundance of functional transcripts.

Jasmonate controls late development stages of petal growth in Arabidopsis thaliana

In Arabidopsis, four homeotic gene classes, A, B, C and E, are required for the patterning of floral organs. However, very little is known about how the activity of these master genes is translated into regulatory processes leading to specific growth patterns and the formation of organs with specific shapes and sizes. Previously we showed that the transcript variant BPEp encodes a bHLH transcription factor that is involved in limiting petal size by controlling post-mitotic cell expansion. Here we show that the phytohormone jasmonate is required for control of BPEp expression. Expression of BPEp was negatively regulated in opr3 mutant flowers that are deficient in jasmonate synthesis. Moreover, the expression of BPEp was restored in opr3 flowers following exogenous jasmonate treatments. Expression of the second transcript variant BPEub, which originates from the same gene as BPEp via an alternative splicing event, was not affected, indicating that BPEp accumulation triggered by jasmonate occurs at the post-transcriptional level. Consistent with these data, opr3 exhibited an increase in petal size as a result of increased cell size, as well as a modified vein pattern, phenotypes that are similar to those of the bpe-1 mutant. Furthermore, exogenous treatments with jasmonate rescued petal phenotypes associated with loss of function of OPR3. Our data demonstrate that jasmonate signaling downstream of OPR3 is involved in the control of cell expansion and in limiting petal size, and that BPEp is a downstream target that functions as a component mediating jasmonate signaling during petal growth.

Genome-wide mapping of alternative splicing in Arabidopsis thaliana

Alternative splicing can enhance transcriptome plasticity and proteome diversity. In plants, alternative splicing can be manifested at different developmental stages, and is frequently associated with specific tissue types or environmental conditions such as abiotic stress. We mapped the Arabidopsis transcriptome at single-base resolution using the Illumina platform for ultrahigh-throughput RNA sequencing (RNA-seq). Deep transcriptome sequencing confirmed a majority of annotated introns and identified thousands of novel alternatively spliced mRNA isoforms. Our analysis suggests that at least approximately 42% of intron-containing genes in Arabidopsis are alternatively spliced; this is significantly higher than previous estimates based on cDNA/expressed sequence tag sequencing. Random validation confirmed that novel splice isoforms empirically predicted by RNA-seq can be detected in vivo. Novel introns detected by RNA-seq were substantially enriched in nonconsensus terminal dinucleotide splice signals. Alternative isoforms with premature termination codons (PTCs) comprised the majority of alternatively spliced transcripts. Using an example of an essential circadian clock gene, we show that intron retention can generate relatively abundant PTC(+) isoforms and that this specific event is highly conserved among diverse plant species. Alternatively spliced PTC(+) isoforms can be potentially targeted for degradation by the nonsense mediated mRNA decay (NMD) surveillance machinery or regulate the level of functional transcripts by the mechanism of regulated unproductive splicing and translation (RUST). We demonstrate that the relative ratios of the PTC(+) and reference isoforms for several key regulatory genes can be considerably shifted under abiotic stress treatments. Taken together, our results suggest that like in animals, NMD and RUST may be widespread in plants and may play important roles in regulating gene expression.

Zn stress-induced inhibition of bean PvSR2-GUS fusion gene splicing is gene-specific in transgenic tobacco

The stress-related gene no. 2 of Phaseolus vulgaris (PvSR2) is metal inducible and contains a single intron. Here, we report that Zn stress inhibited the splicing of the PvSR2-beta-glucuronidase (GUS) fusion gene in a concentration- and time-dependent manner in tobacco seedlings. The inhibition appears to be specific for the PvSR2-GUS transgene: splicing of four endogenous tobacco genes was unaffected by Zn stress. Our results provide in vivo evidence that Zn stress-dependent intron retention is transgene specific in plants.

Post-transcriptional regulation of gene expression in plants during abiotic stress

Land plants are anchored in one place for most of their life cycle and therefore must constantly adapt their growth and metabolism to abiotic stresses such as light intensity, temperature and the availability of water and essential minerals. Thus, plants’ subsistence depends on their ability to regulate rapidly gene expression in order to adapt their physiology to their environment. Recent studies indicate that post-transcriptional regulations of gene expression play an important role in how plants respond to abiotic stresses. We will review the different mechanisms of post-transcriptional regulation of nuclear genes expression including messenger RNA (mRNA) processing, stability, localization and protein translation, and discuss their relative importance for plant adaptation to abiotic stress.

Two alternatively spliced isoforms of the Arabidopsis SR45 protein have distinct roles during normal plant development

The serine-arginine-rich (SR) proteins constitute a conserved family of pre-mRNA splicing factors. In Arabidopsis (Arabidopsis thaliana), they are encoded by 19 genes, most of which are themselves alternatively spliced. In the case of SR45, the use of alternative 3' splice sites 21 nucleotides apart generates two alternatively spliced isoforms. Isoform 1 (SR45.1) has an insertion relative to isoform 2 (SR45.2) that replaces a single arginine with eight amino acids (TSPQRKTG). The biological implications of SR45 alternative splicing have been unclear. A previously described loss-of-function mutant affecting both isoforms, sr45-1, shows several developmental defects, including defects in petal development and root growth. We found that the SR45 promoter is highly active in regions with actively growing and dividing cells. We also tested the ability of each SR45 isoform to complement the sr45-1 mutant by overexpression of isoform-specific green fluorescent protein (GFP) fusion proteins. As expected, transgenic plants overexpressing either isoform displayed both nuclear speckles and GFP fluorescence throughout the nucleoplasm. We found that SR45.1-GFP complements the flower petal phenotype, but not the root growth phenotype. Conversely, SR45.2-GFP complements root growth but not floral morphology. Mutation of a predicted phosphorylation site within the alternatively spliced segment, SR45.1-S219A-GFP, does not affect complementation. However, a double mutation affecting both serine-219 and the adjacent threonine-218 (SR45.1-T218A + S219A-GFP) behaves like isoform 2, complementing the root but not the floral phenotype. In conclusion, our study provides evidence that the two alternatively spliced isoforms of SR45 have distinct biological functions.

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

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

Genetic analysis of functional redundancy of BRM ATPase and ATSWI3C subunits of Arabidopsis SWI/SNF chromatin remodelling complexes

In yeast and mammals, ATP-dependent chromatin remodelling complexes of the SWI/SNF family play critical roles in the regulation of transcription, cell proliferation, differentiation and development. Homologues of conserved subunits of SWI/SNF-type complexes, including Snf2-type ATPases and SWI3-type proteins, participate in analogous processes in Arabidopsis. Recent studies indicate a remarkable similarity between phenotypic effects of mutations in the SWI3 homologue ATSWI3C and bromodomain-ATPase BRM genes. To verify the extent of functional similarity between BRM and ATSWI3C, we have constructed atswi3c brm double mutants and compared their phenotypic traits to those of simultaneously grown single atswi3c and brm mutants. In addition to inheritance of characteristic developmental abnormalities shared by atswi3c and brm mutants, some additive brm-specific traits were also observed in the atswi3c brm double mutants. Unlike atswi3c, the brm mutation results in the enhancement of abnormal carpel development and pollen abortion leading to complete male sterility. Despite the overall similarity of brm and atswi3c phenotypes, a critical requirement for BRM in the differentiation of reproductive organs suggests that its regulatory functions do not entirely overlap those of ATSWI3C. The detection of two different transcript isoforms indicates that BRM is regulated by alternative splicing that creates an in-frame premature translation stop codon in its SNF2-like ATPase coding domain. The analysis of Arabidopsis mutants in nonsense-mediated decay suggests an involvement of this pathway in the control of alternative BRM transcript level.

Flooding stress-induced glycine-rich RNA-binding protein from Nicotiana tabacum

A cDNA clone for a transcript preferentially expressed during an early phase of flooding was isolated from Nicotiana tabacum. Nucleotide sequencing of the cDNA clone identified an open reading frame that has high homology to the previously reported glycine-rich RNA-binding proteins. The open reading frame consists of 157 amino acids with an N-terminal RNA-recognition motif and a C-terminal glycine-rich domain, and thus the cDNA clone was designated as Nicotiana tabaccum glycine-rich RNA-binding protein-1 (NtGRP1). Expression of NtGRP1 was upregulated under flooding stress and also increased, but at much lower levels, under conditions of cold, drought, heat, high salt content, and abscisic acid treatment. RNA homopolymer-binding assay showed that NtGRP1 binds to all the RNA homopolymers tested with a higher affinity to poly r(G) and poly r(A) than to poly r(U) and poly r(C). Nucleic acid-binding assays showed that NtGRP1 binds to ssDNA, dsDNA, and mRNA. NtGRP1 suppressed expression of the fire luciferase gene in vitro, and the suppression of luciferase gene expression could be rescued by addition of oligonucleotides. Collectively, the data suggest NtGRP1 as a negative modulator of gene expression by binding to DNA or RNA in bulk that could be advantageous for plants in a stress condition like flooding.

Reciprocal regulation of glycine-rich RNA-binding proteins via an interlocked feedback loop coupling alternative splicing to nonsense-mediated decay in Arabidopsis

The Arabidopsis RNA-binding protein AtGRP8 undergoes negative autoregulation at the post-transcriptional level. An elevated AtGRP8 protein level promotes the use of a cryptic 5' splice site to generate an alternatively spliced transcript, as_AtGRP8, retaining the 5' half of the intron with a premature termination codon. In mutants defective in nonsense-mediated decay (NMD) abundance of as_AtGRP8 but not its pre-mRNA is elevated, indicating that as_AtGRP8 is a direct NMD target, thus limiting the production of functional AtGRP8 protein. In addition to its own pre-mRNA, AtGRP8 negatively regulates the AtGRP7 transcript through promoting the formation of the equivalent alternatively spliced as_AtGRP7 transcript, leading to a decrease in AtGRP7 abundance. Recombinant AtGRP8 binds to its own and the AtGRP7 pre-mRNA, suggesting that this interaction is relevant for the splicing decision in vivo. AtGRP7 itself is part of a negative autoregulatory circuit that influences circadian oscillations of its own and the AtGRP8 transcript through alternative splicing linked to NMD. Thus, we identify an interlocked feedback loop through which two RNA-binding proteins autoregulate and reciprocally crossregulate by coupling unproductive splicing to NMD. A high degree of evolutionary sequence conservation in the introns retained in as_AtGRP8 or as_AtGRP7 points to an important function of these sequences.

Differential SAGE analysis in Arabidopsis uncovers increased transcriptome complexity in response to low temperature

Background

Abiotic stress, including low temperature, limits the productivity and geographical distribution of plants, which has led to significant interest in understanding the complex processes that allow plants to adapt to such stresses. The wide range of physiological, biochemical and molecular changes that occur in plants exposed to low temperature require a robust global approach to studying the response. We have employed Serial Analysis of Gene Expression (SAGE) to uncover changes in the transcriptome of Arabidopsis thaliana over a time course of low temperature stress.

Results

Five SAGE libraries were generated from A. thaliana leaf tissue collected at time points ranging from 30 minutes to one week of low temperature treatment (4°C). Over 240,000 high quality SAGE tags, corresponding to 16,629 annotated genes, provided a comprehensive survey of changes in the transcriptome in response to low temperature, from perception of the stress to acquisition of freezing tolerance. Interpretation of these data was facilitated by representing the SAGE data by gene identifier, allowing more robust statistical analysis, cross-platform comparisons and the identification of genes sharing common expression profiles. Simultaneous statistical calculations across all five libraries identified 920 low temperature responsive genes, only 24% of which overlapped with previous global expression analysis performed using microarrays, although similar functional categories were affected. Clustering of the differentially regulated genes facilitated the identification of novel loci correlated with the development of freezing tolerance. Analysis of their promoter sequences revealed subsets of genes that were independent of CBF and ABA regulation and could provide a mechanism for elucidating complementary signalling pathways. The SAGE data emphasised the complexity of the plant response, with alternate pre-mRNA processing events increasing at low temperatures and antisense transcription being repressed.

Conclusion

Alternate transcript processing appears to play an important role in enhancing the plasticity of the stress induced transcriptome. Novel genes and cis-acting sequences have been identified as compelling targets to allow manipulation of the plant's ability to protect against low temperature stress. The analyses performed provide a contextual framework for the interpretation of quantitative sequence tag based transcriptome analysis which will prevail with the application of next generation sequencing technology.

Heat stress-responsive transcriptome analysis in heat susceptible and tolerant wheat (Triticum aestivum L.) by using Wheat Genome Array

Background

Wheat is a major crop in the world, and the high temperature stress can reduce the yield of wheat by as much as 15%. The molecular changes in response to heat stress are poorly understood. Using GeneChip^® Wheat Genome Array, we analyzed genome-wide gene expression profiles in the leaves of two wheat genotypes, namely, heat susceptible 'Chinese Spring' (CS) and heat tolerant 'TAM107' (TAM).

Results

A total of 6560 (~10.7%) probe sets displayed 2-fold or more changes in expression in at least one heat treatment (false discovery rate, FDR, α = 0.001). Except for heat shock protein (HSP) and heat shock factor (HSF) genes, these putative heat responsive genes encode transcription factors and proteins involved in phytohormone biosynthesis/signaling, calcium and sugar signal pathways, RNA metabolism, ribosomal proteins, primary and secondary metabolisms, as well as proteins related to other stresses. A total of 313 probe sets were differentially expressed between the two genotypes, which could be responsible for the difference in heat tolerance of the two genotypes. Moreover, 1314 were differentially expressed between the heat treatments with and without pre-acclimation, and 4533 were differentially expressed between short and prolonged heat treatments.

Conclusion

The differences in heat tolerance in different wheat genotypes may be associated with multiple processes and mechanisms involving HSPs, transcription factors, and other stress related genes. Heat acclimation has little effects on gene expression under prolonged treatments but affects gene expression in wheat under short-term heat stress. The heat stress responsive genes identified in this study will facilitate our understanding of molecular basis for heat tolerance in different wheat genotypes and future improvement of heat tolerance in wheat and other cereals.

Spliceosomal proteins in plants

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

Spatiotemporal organization of pre-mRNA splicing proteins in plants

The general organization ofeukaryotic nuclei, including plant nuclei, into functional domains is now widely recognized. Conventional immunocytochemistry and visualization of proteins fused to fluorescent proteins (FP) have revealed that in plants, RNA and protein components of pre-mRNA splicing are spatially organized depending on the stage of cell cycle, development, and the cell's physiological state. Application of some of the latest microscopy techniques, which reveal biophysical properties such as diffusion and interaction properties of proteins, has begun to provide important insights into the functional organization of spliceosomal proteins in plants. Although some progress has been made in understanding the spatial and temporal organization of splicing machinery in plants, the mechanisms that regulate this organization and its functional consequences remain unresolved.

Regulation of alternative splicing of pre-mRNAs by stresses

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

Genome-wide analyses of alternative splicing in plants: opportunities and challenges

Alternative splicing (AS) creates multiple mRNA transcripts from a single gene. While AS is known to contribute to gene regulation and proteome diversity in animals, the study of its importance in plants is in its early stages. However, recently available plant genome and transcript sequence data sets are enabling a global analysis of AS in many plant species. Results of genome analysis have revealed differences between animals and plants in the frequency of alternative splicing. The proportion of plant genes that have one or more alternative transcript isoforms is approximately 20%, indicating that AS in plants is not rare, although this rate is approximately one-third of that observed in human. The majority of plant AS events have not been functionally characterized, but evidence suggests that AS participates in important plant functions, including stress response, and may impact domestication and trait selection. The increasing availability of plant genome sequence data will enable larger comparative analyses that will identify functionally important plant AS events based on their evolutionary conservation, determine the influence of genome duplication on the evolution of AS, and discover plant-specific cis-elements that regulate AS. This review summarizes recent analyses of AS in plants, discusses the importance of further analysis, and suggests directions for future efforts.

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

Alignment of transcripts and genome sequences yielded a set of alternatively spliced transcripts in four angiosperm genomes: three dicotyledon species Medicago truncatula (Medicago), Populus trichocarpa (poplar) and Arabidopsis thaliana (Arabidopsis), and the monocotyledon Oryzae sativa (rice). Intron retention was the predominant mode of alternative splicing (AS) in each species, consistent with previous reports for Arabidopsis and rice. We analyzed the structure of 5'-splice junctions and observed commonalities between species. There was dependency of base composition between sites flanking the 5'-splice junction, with the potential to create a subset of splice sites that interact more weakly or strongly than average with U1 snRNA. Such altered nucleotide composition was correlated with splicing fidelity in all four species. For Medicago, poplar and Arabidopsis, but not in rice, alternative splicing was most prevalent for introns with decreased UA content, consistent with lower UA content for monocot introns and potentially reflecting evolved differences in splicing mechanisms. Similarly, the occurrence of AS between transcript Gene Ontology categories was positively correlated between Arabidopsis and Medicago, with no correlation between dicots and rice. Analysis of within-species paralogs and between-species reciprocal best-hit homologs yielded rare cases of potentially conserved AS events. Reverse transcriptase PCR and amplicon sequencing were used to confirm a subset of the in silico-predicted AS events within Medicago, as well as to characterize conserved AS events between Medicago and Arabidopsis.

Alternative splicing at NAGNAG acceptors in Arabidopsis thaliana SR and SR-related protein-coding genes

Background

Several recent studies indicate that alternative splicing in Arabidopsis and other plants is a common mechanism for post-transcriptional modulation of gene expression. However, few analyses have been done so far to elucidate the functional relevance of alternative splicing in higher plants. Representing a frequent and universal subtle alternative splicing event among eukaryotes, alternative splicing at NAGNAG acceptors contributes to transcriptome diversity and therefore, proteome plasticity. Alternatively spliced NAGNAG acceptors are overrepresented in genes coding for proteins with RNA-recognition motifs (RRMs). As SR proteins, a family of RRM-containing important splicing factors, are known to be extensively alternatively spliced in Arabidopsis, we analyzed alternative splicing at NAGNAG acceptors in SR and SR-related genes.

Results

In a comprehensive analysis of the Arabidopsis thaliana genome, we identified 6,772 introns that exhibit a NAGNAG acceptor motif. Alternative splicing at these acceptors was assessed using available EST data, complemented by a sequence-based prediction method. Of the 36 identified introns within 30 SR and SR-related protein-coding genes that have a NAGNAG acceptor, we selected 15 candidates for an experimental analysis of alternative splicing under several conditions. We provide experimental evidence for 8 of these candidates being alternatively spliced. Quantifying the ratio of NAGNAG-derived splice variants under several conditions, we found organ-specific splicing ratios in adult plants and changes in seedlings of different ages. Splicing ratio changes were observed in response to heat shock and most strikingly, cold shock. Interestingly, the patterns of differential splicing ratios are similar for all analyzed genes.

Conclusion

NAGNAG acceptors frequently occur in the Arabidopsis genome and are particularly prevalent in SR and SR-related protein-coding genes. A lack of extensive EST coverage can be compensated by using the proposed sequence-based method to predict alternative splicing at these acceptors. Our findings indicate that the differential effects on NAGNAG alternative splicing in SR and SR-related genes are organ- and condition-specific rather than gene-specific.

Redox and heavy metal effects on the biochemical activities of an Arabidopsis polyadenylation factor subunit

The Arabidopsis CPSF30 ortholog (AtCPSF30) is an RNA-binding endonuclease that is part of the plant polyadenylation complex. Previous work (B. Addepalli, A.G. Hunt, Nucleic Acids Res. 35 (2007) 4453-4463) demonstrated that different zinc finger motifs in the protein were responsible for RNA-binding and nuclease activity, respectively. In this study, a more detailed functional map of AtCPSF30 is presented, a map that includes descriptions of novel biochemical activities. Elevated temperatures, the specific zinc chelator 1,10-phenanthroline, and the sulfhydryl reagent dithiothreitol all had differential inhibitory effects on the RNA-binding and nuclease activities. The endonuclease activity of AtCPSF30 was inhibited by relatively high (>100muM) concentrations of zinc, and this inhibition required a plant-specific N-terminal domain apart from the zinc finger core of the protein. ATP stimulated the nuclease activity in the presence of zinc, and this stimulation required a plant-specific C-terminal domain, again apart from the zinc finger core. These studies reveal a subtle and unexpected complexity to AtCPSF30, and raise the possibility that multiple avenues of regulation may impinge on this protein through different functional domains.

Alternative splicing regulation during C. elegans development: splicing factors as regulated targets

Alternative splicing generates protein diversity and allows for post-transcriptional gene regulation. Estimates suggest that 10% of the genes in Caenorhabditis elegans undergo alternative splicing. We constructed a splicing-sensitive microarray to detect alternative splicing for 352 cassette exons and tested for changes in alternative splicing of these genes during development. We found that the microarray data predicted that 62/352 (∼18%) of the alternative splicing events studied show a strong change in the relative levels of the spliced isoforms (>4-fold) during development. Confirmation of the microarray data by RT-PCR was obtained for 70% of randomly selected genes tested. Among the genes with the most developmentally regulated alternatively splicing was the hnRNP F/H splicing factor homolog, W02D3.11 – now named hrpf-1. For the cassette exon of hrpf-1, the inclusion isoform comprises 65% of hrpf-1 steady state messages in embryos but only 0.1% in the first larval stage. This dramatic change in the alternative splicing of an alternative splicing factor suggests a complex cascade of splicing regulation during development. We analyzed splicing in embryos from a strain with a mutation in the splicing factor sym-2, another hnRNP F/H homolog. We found that approximately half of the genes with large alternative splicing changes between the embryo and L1 stages are regulated by sym-2 in embryos. An analysis of the role of nonsense-mediated decay in regulating steady-state alternative mRNA isoforms was performed. We found that 8% of the 352 events studied have alternative isoforms whose relative steady-state levels in embryos change more than 4-fold in a nonsense-mediated decay mutant, including hrpf-1. Strikingly, 53% of these alternative splicing events that are affected by NMD in our experiment are not obvious substrates for NMD based on the presence of premature termination codons. This suggests that the targeting of splicing factors by NMD may have downstream effects on alternative splicing regulation.

Alternative splicing is a mechanism for generating more than one messenger RNA from a given gene. The alternative transcripts can encode different proteins that share some regions in common but have modified functions, thus increasing the number of proteins encoded by the genome. Alternative splicing can also lead to the production of mRNA isoforms that are then subject to degradation by the nonsense-mediated decay pathway, thus providing a mechanism to down-regulate gene expression without decreasing transcription. Examples of cell type-specific, hormone-responsive, and developmentally-regulated alternative splicing have been described. We decided to measure the extent of developmentally regulated alternative splicing in the nematode model organism Caenorhabditis elegans. We developed a DNA microarray that can measure the alternative splicing of 352 cassette exons simultaneously and used it to probe alternative splicing in RNA extracted from embryos, the four larval stages, and adults. We show that 18% of the alternatively spliced genes tested show >4-fold changes in alternative splicing during development. In addition, we show that one of the most regulated genes is itself a splicing factor, providing support for a model in which a cascade of alternative splicing regulation occurs during development.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Plant SR proteins and their functions

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

Cold stress regulation of gene expression in plants

Cold stress adversely affects plant growth and development. Most temperate plants acquire freezing tolerance by a process called cold acclimation. Here, we focus on recent progress in transcriptional, post-transcriptional and post-translational regulation of gene expression that is critical for cold acclimation. Transcriptional regulation is mediated by the inducer of C-repeat binding factor (CBF) expression 1 (ICE1), the CBF transcriptional cascade and CBF-independent regulons during cold acclimation. ICE1 is negatively regulated by ubiquitination-mediated proteolysis and positively regulated by SUMO (small ubiquitin-related modifier) E3 ligase-catalyzed sumoylation. Post-transcriptional regulatory mechanisms, such as pre-mRNA splicing, mRNA export and small RNA-directed mRNA degradation, also play important roles in cold stress responses.

Differential expression of alternatively spliced mRNAs of Arabidopsis SR protein homologs, atSR30 and atSR45a, in response to environmental stress

Serine/arginine-rich (SR) proteins are associated with either the regulation or the execution of both constitutive splicing and the selection of alternative splice sites in animals and plants. We demonstrated the molecular characterization of a homolog of SR protein, atSR45a, in Arabidopsis plants. Six types of mRNA variants (atSR45a-1a-e and atSR45a-2) were generated by the alternative selection of transcriptional initiation sites and the alternative splicing of introns in atSR45a pre-mRNA. The atSR45a-1a and -2 proteins, presumed mature forms, were located in the nucleus and interacted with U1-70K, suggesting that these proteins function as a splicing factor in Arabidopsis. The levels of the transcripts atSR45a and atSR30, SF2/ASF-like SR proteins, were increased by various types of stress, such as high-light irradiation and salinity. Furthermore, the splicing patterns of atSR45a and atSR30 pre-mRNA themselves were altered under these stressful conditions. In particular, the expression of atSR45a-1a, atSR45a-2, atSR30 mRNA1 and atSR30 mRNA3 was greatly increased by high-light irradiation. These results indicate that the regulation of transcription and alternative splicing of atSR45a and atSR30 is responsive to various stressful conditions.

Regulation of plant developmental processes by a novel splicing factor

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

Alternative splicing of pre-messenger RNAs in plants in the genomic era

Primary transcripts (precursor-mRNAs) with introns can undergo alternative splicing to produce multiple transcripts from a single gene by differential use of splice sites, thereby increasing the transcriptome and proteome complexity within and between cells and tissues. Alternative splicing in plants is largely an unexplored area of gene expression, as this phenomenon used to be considered rare. However, recent genome-wide computational analyses have revealed that alternative splicing in flowering plants is far more prevalent than previously thought. Interestingly, pre-mRNAs of many spliceosomal proteins, especially serine/arginine-rich (SR) proteins, are extensively alternatively spliced. Furthermore, stresses have a dramatic effect on alternative splicing of pre-mRNAs including those that encode many spliceosomal proteins. Although the mechanisms that regulate alternative splicing in plants are largely unknown, several reports strongly suggest a key role for SR proteins in spliceosome assembly and regulated splicing. Recent studies suggest that alternative splicing in plants is an important posttranscriptional regulatory mechanism in modulating gene expression and eventually plant form and function.