Intron retention is a major phenomenon in alternative splicing in Arabidopsis

Cross-kingdom patterns of alternative splicing and splice recognition

BACKGROUND

Variations in transcript splicing can reveal how eukaryotes recognize intronic splice sites. Retained introns (RIs) commonly appear when the intron definition (ID) mechanism of splice site recognition inconsistently identifies intron-exon boundaries, and cassette exons (CEs) are often caused by variable recognition of splice junctions by the exon definition (ED) mechanism. We have performed a comprehensive survey of alternative splicing across 42 eukaryotes to gain insight into how spliceosomal introns are recognized.

RESULTS

All eukaryotes we studied exhibit RIs, which appear more frequently than previously thought. CEs are also present in all kingdoms and most of the organisms in our analysis. We observe that the ratio of CEs to RIs varies substantially among kingdoms, while the ratio of competing 3' acceptor and competing 5' donor sites remains nearly constant. In addition, we find the ratio of CEs to RIs in each organism correlates with the length of its introns. In all 14 fungi we examined, as well as in most of the 9 protists, RIs far outnumber CEs. This differs from the trend seen in 13 multicellular animals, where CEs occur much more frequently than RIs. The six plants we analyzed exhibit intermediate proportions of CEs and RIs.

CONCLUSION

Our results suggest that most extant eukaryotes are capable of recognizing splice sites via both ID and ED, although ED is most common in multicellular animals and ID predominates in fungi and most protists.

Monitoring changes in alternative precursor messenger RNA splicing in multiple gene transcripts

Alternative splicing (AS) increases the proteomic and functional capacity of genomes through the generation of alternative mRNA transcripts from the same gene. AS is now estimated to occur in a third of Arabidopsis and rice genes, and includes genes involved in the control of growth and development, responses to stress and signalling. Regulation of AS reflects the interactions between positive and negative cis sequences in the precursor messenger RNA and a range of trans-acting factors. The levels and activities of these factors differ in different cells and growth conditions. To identify changes in AS in multiple genes simultaneously, we have established a reproducible RT-PCR panel that can analyse 96 alternative splicing events and accurately measure the ratio of alternatively spliced products. This procedure detected statistically significant changes in AS in different plant organs, in plants grown under different light and day-length conditions, and in plants overexpressing splicing factors. The system provides a convenient, medium-throughput means of monitoring changes in AS in multiple genes. It can readily be applied to much larger or targeted sets of gene transcripts to generate information on the significance and regulation of AS in plant growth and development, specific processes and responses to external stimuli.

Utilizing tiling microarrays for whole-genome analysis in plants

The recent explosion in available genome sequence data has ushered in an era in which analysis of a whole genome can be performed in a single experiment. While DNA microarrays have long been the established technology for measuring gene expression levels, standard expression arrays use relatively few probes for each gene and are typically biased toward known and predicted gene structures. Recently, with the availability of complete genome sequences for many organisms, very-high-density oligonucleotide-based microarrays that span the entire genome have emerged as the preferred platform for genomic analysis. Whole-genome tiling microarrays can be employed for a myriad of purposes, including empirical annotation of the transcriptome, chromatin immunoprecipitation-chip studies, analysis of alternative RNA splicing, characterization of the methylation state of cytosine bases throughout a genome (methylome), and DNA polymorphism discovery. Here, we review several applications of whole-genome technology to obtain a variety of genomic-scale information in plants.

A first glimpse at the transcriptome of Physarum polycephalum

Background

Physarum polycephalum, an acellular plasmodial species belongs to the amoebozoa, a major branch in eukaryote evolution. Its complex life cycle and rich cell biology is reflected in more than 2500 publications on various aspects of its biochemistry, developmental biology, cytoskeleton, and cell motility. It now can be genetically manipulated, opening up the possibility of targeted functional analysis in this organism.

Methods

Here we describe a large fraction of the transcribed genes by sequencing a cDNA library from the plasmodial stage of the developmental cycle.

Results

In addition to the genes for the basic metabolism we found an unexpected large number of genes involved in sophisticated signaling networks and identified potential receptors for environmental signals such as light. In accordance with the various developmental options of the plasmodial cell we found that many P. polycephalum genes are alternatively spliced. Using 30 donor and 30 acceptor sites we determined the splicing signatures of this species.

Comparisons to various other organisms including Dictyostelium, the closest relative, revealed that roughly half of the transcribed genes have no detectable counterpart, thus potentially defining species specific adaptations. On the other hand, we found highly conserved proteins, which are maintained in the metazoan lineage, but absent in D. discoideum or plants. These genes arose possibly in the last common ancestor of Amoebozoa and Metazoa but were lost in D. discoideum.

Conclusion

This work provides an analysis of up to half of the protein coding genes of Physarum polycephalum. The definition of splice motifs together with the description of alternatively spliced genes will provide a valuable resource for the ongoing genome project.

Analysis of alternative splicing in plants with bioinformatics tools

Alternative splicing is a molecular mechanism utilized by a broad range of eukaryotes to extend the repertoire of functions encoded by single genes and to posttranscriptionally regulate gene expression. Recent analyses of expressed transcript sequences aligned to the complete genomes of Arabidopsis and rice indicate that alternative splicing in plants is prevalent and exhibits several features similar to other higher eukaryotes including mouse and human. This chapter reviews the computational strategies employed to study alternative splicing with bioinformatics tools and the recent findings from analyses performed on plants by applying such methods.

Regulation of splicing by protein phosphorylation

Most eukaryotic messenger RNAs are transcribed as precursors that necessitate specific and exact processing of intron boundaries. Furthermore, the choice of these boundaries appears to be fluid and adaptive to the rate of transcription and the developmental and physiological state of the cell. A central regulator of splicing reactions and choice are kinases that work through phosphorylation of specific factors like RNA polymerase II, which influences the pace of transcription and of SR splicing factors. While very different in their mechanisms both regulatory pathways will impact on splicing site choice. This chapter summarizes the biology of splicing-related phosphorylation activity, emphasizing plant-specific aspects in relation to the metazoan counterpart.

Alternative splicing in plant defense

Plant resistance proteins directly or indirectly perceive the presence of pathogen virulence factors and trigger an effective form of plant immunity that often includes programmed host cell death. Because the activation of resistance proteins has the potential to be detrimental to the plant, this process is tightly regulated on multiple levels. Several resistance genes have been shown to be alternatively spliced. Depending on the resistance gene, alternative transcripts are thought to limit the expression of R proteins or encode truncated R proteins with a positive role in defense activation. In addition, R gene alternative splicing is dynamic during the defense response. Possible mechanisms of R gene alternative splicing regulation and how alternative R gene transcripts fit into the current view of resistance protein-mediated defense responses are discussed.

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.

Bioinformatics detection of alternative splicing

In recent years, genome-wide detection of alternative splicing based on Expressed Sequence Tag (EST) sequence alignments with mRNA and genomic sequences has dramatically expanded our understanding of the role of alternative splicing in functional regulation. This chapter reviews the data, methodology, and technical challenges of these genome-wide analyses of alternative splicing, and briefly surveys some of the uses to which such alternative splicing databases have been put. For example, with proper alternative splicing database schema design, it is possible to query genome-wide for alternative splicing patterns that are specific to particular tissues, disease states (e.g., cancer), gender, or developmental stages. EST alignments can be used to estimate exon inclusion or exclusion level of alternatively spliced exons and evolutionary changes for various species can be inferred from exon inclusion level. Such databases can also help automate design of probes for RT-PCR and microarrays, enabling high throughput experimental measurement of alternative splicing.

Alternative splicing and mRNA levels of the disease resistance gene RPS4 are induced during defense responses

The Arabidopsis (Arabidopsis thaliana) disease resistance protein RESISTANCE TO PSEUDOMONAS SYRINGAE4 (RPS4) activates defenses in response to bacterial pathogens expressing avrRps4 in a gene-for-gene specific manner. The RPS4 gene produces multiple transcripts via alternative splicing of two regular introns flanking exon 3 and a cryptic intron within exon 3. We showed previously that RPS4-mediated resistance requires the combined presence of transcripts encoding both full-length and truncated open reading frames. Here, we demonstrate that alternative splicing of RPS4 undergoes dynamic changes specifically during the resistance response. Furthermore, RPS4 expression was induced by the presence of AvrRps4 in an EDS1-dependent manner. Interestingly, inducible alternative splicing was not limited to the avrRps4-RPS4 interaction, indicating that regulation of alternative splicing may be a general response to prime the plant stress response system. Intron-deficient transgenes lacking only one intron were previously shown to be nonfunctional. Here, we establish quantitatively that the absence of one intron had no effect on the splicing frequency of remaining introns. Given the lack of functionality of single intron-deficient transgenes, this suggests that the products of individual transcripts have distinct functions during RPS4-triggered resistance. Transient expression of truncated RPS4 proteins in Nicotiana benthamiana induced hypersensitive response-like cell death in the absence of AvrRps4. Interestingly, different truncated proteins had markedly differing stability. In summary, RPS4 function is regulated at multiple levels, including gene expression, alternative splicing, and protein stability, presumably to fine-tune activity and limit damage inflicted by activated RPS4 protein.

Ascorbate peroxidase gene family in tomato: its identification and characterization

The antioxidative response, where ascorbate peroxidase (APX) is a key enzyme, is an integral part of the plant tolerance response to environmental stresses. As a first step towards the study of the physiological role and the regulation of the members of the Apx gene family, the orthologs of the stress-sensitive cultivated tomato Solanum lycopersicum cv. M82 (Slm) and of the wild salt-tolerant species S. pennellii acc. Atico (Spa) were identified by utilizing the tomato EST database, and characterized. A redundant list of 16 virtual Apx transcripts and four singleton ESTs was shown to correspond to seven genuine Apx genes. The complete tomato Apx gene family is comprised of genes encoding three cytosolic, two peroxisomal, and two chloroplastic APXs. These genes attained differential regulatory patterns in various Slm organs. More detailed study of Apx1 and Apx2 genes, that are the products of a recent gene duplication event, shows that they have already attained differential regulation within and between Slm and Spa under control and stress conditions. It is also suggested that due to lineage-specific gene duplication and lose events, intricate phylogenetic relationships exist among the members of the Apx gene families.

Characterization of Arabidopsis thaliana genes encoding functional homologues of the yeast metal chaperone Cox19p, involved in cytochrome c oxidase biogenesis

The Arabidopsis thaliana genome contains two nearly identical genes which encode proteins showing similarity with the yeast metal chaperone Cox19p, involved in cytochrome c oxidase biogenesis. One of these genes (AtCOX19-1) produces two transcript forms that arise from an alternative splicing event and encode proteins with different N-terminal portions. Both AtCOX19 isoforms are imported into mitochondria in vitro and are found attached to the inner membrane facing the intermembrane space. The smaller AtCOX19-1 isoform, but not the larger one, is able to restore growth on non-fermentable carbon sources when expressed in a yeast cox19 null mutant. AtCOX19 transcript levels increase by treatment with copper or compounds that produce reactive oxygen species. Young roots and anthers are highly stained in AtCOX19-1::GUS plants. Expression in leaves is only observed when cuts are produced, suggesting an induction by wounding. Infection of plants with the pathogenic bacterium Pseudomonas syringae pv. tomato also induces AtCOX19 gene expression. The results suggest that AtCOX19 genes encode functional homologues of the yeast metal chaperone. Induction by biotic and abiotic stress factors may indicate a relevant role of this protein in the biogenesis of cytochrome c oxidase to replace damaged forms of the enzyme or a more general role in the response of plants to stress.

Gene discovery and transcript analyses in the corn smut pathogen Ustilago maydis: expressed sequence tag and genome sequence comparison

Background

Ustilago maydis is the basidiomycete fungus responsible for common smut of corn and is a model organism for the study of fungal phytopathogenesis. To aid in the annotation of the genome sequence of this organism, several expressed sequence tag (EST) libraries were generated from a variety of U. maydis cell types. In addition to utility in the context of gene identification and structure annotation, the ESTs were analyzed to identify differentially abundant transcripts and to detect evidence of alternative splicing and anti-sense transcription.

Results

Four cDNA libraries were constructed using RNA isolated from U. maydis diploid teliospores (U. maydis strains 518 × 521) and haploid cells of strain 521 grown under nutrient rich, carbon starved, and nitrogen starved conditions. Using the genome sequence as a scaffold, the 15,901 ESTs were assembled into 6,101 contiguous expressed sequences (contigs); among these, 5,482 corresponded to predicted genes in the MUMDB (MIPS Ustilago maydis database), while 619 aligned to regions of the genome not yet designated as genes in MUMDB. A comparison of EST abundance identified numerous genes that may be regulated in a cell type or starvation-specific manner. The transcriptional response to nitrogen starvation was assessed using RT-qPCR. The results of this suggest that there may be cross-talk between the nitrogen and carbon signalling pathways in U. maydis. Bioinformatic analysis identified numerous examples of alternative splicing and anti-sense transcription. While intron retention was the predominant form of alternative splicing in U. maydis, other varieties were also evident (e.g. exon skipping). Selected instances of both alternative splicing and anti-sense transcription were independently confirmed using RT-PCR.

Conclusion

Through this work: 1) substantial sequence information has been provided for U. maydis genome annotation; 2) new genes were identified through the discovery of 619 contigs that had previously escaped annotation; 3) evidence is provided that suggests the regulation of nitrogen metabolism in U. maydis differs from that of other model fungi, and 4) Alternative splicing and anti-sense transcription were identified in U. maydis and, amid similar observations in other basidiomycetes, this suggests these phenomena may be widespread in this group of fungi. These advances emphasize the importance of EST analysis in genome annotation.

Comparative cross-species alternative splicing in plants

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

Alternative splicing of pre-mRNA in plants

Alternative splicing is an important cellular mechanism that increases the diversity of gene products. The study of alternatively spliced genes reported so far in plants is far less documented than that in mammals, but considerable results have been reported, showing the role of these genes in regulating mechanisms, influencing factors, and specificities and function of alternative splicing in plants. This review summarizes briefly the major progress made on alternative splicing in plants.

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

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

Plant Gene and Alternatively Spliced Variant Annotator. A plant genome annotation pipeline for rice gene and alternatively spliced variant identification with cross-species expressed sequence tag conservation from seven plant species

The completion of the rice (Oryza sativa) genome draft has brought unprecedented opportunities for genomic studies of the world's most important food crop. Previous rice gene annotations have relied mainly on ab initio methods, which usually yield a high rate of false-positive predictions and give only limited information regarding alternative splicing in rice genes. Comparative approaches based on expressed sequence tags (ESTs) can compensate for the drawbacks of ab initio methods because they can simultaneously identify experimental data-supported genes and alternatively spliced transcripts. Furthermore, cross-species EST information can be used to not only offset the insufficiency of same-species ESTs but also derive evolutionary implications. In this study, we used ESTs from seven plant species, rice, wheat (Triticum aestivum), maize (Zea mays), barley (Hordeum vulgare), sorghum (Sorghum bicolor), soybean (Glycine max), and Arabidopsis (Arabidopsis thaliana), to annotate the rice genome. We developed a plant genome annotation pipeline, Plant Gene and Alternatively Spliced Variant Annotator (PGAA). Using this approach, we identified 852 genes (931 isoforms) not annotated in other widely used databases (i.e. the Institute for Genomic Research, National Center for Biotechnology Information, and Rice Annotation Project) and found 87% of them supported by both rice and nonrice EST evidence. PGAA also identified more than 44,000 alternatively spliced events, of which approximately 20% are not observed in the other three annotations. These novel annotations represent rich opportunities for rice genome research, because the functions of most of our annotated genes are currently unknown. Also, in the PGAA annotation, the isoforms with non-rice-EST-supported exons are significantly enriched in transporter activity but significantly underrepresented in transcription regulator activity. We have also identified potential lineage-specific and conserved isoforms, which are important markers in evolutionary studies. The data and the Web-based interface, RiceViewer, are available for public access at http://RiceViewer.genomics.sinica.edu.tw/.

Retained introns increase putative microRNA targets within 3' UTRs of human mRNA

MicroRNAs (miRNAs) are a class of non-coding RNA that post-transcriptionally regulates the expression of target genes by binding to mRNAs. As one form of alternative splicing, intron retention has influence upon mRNA modification and protein encoding. The effect of miRNA on mRNA containing retained intron within 3' UTR, however, has not been systematically elucidated. Here, we examined a total of 2864 human genes which contain at least one retained intron from the MAASE and ASD databases and found 387 genes having contained retained introns within 3' UTR. The effect of retained introns upon miRNA targets was explored with three web-based programs for miRNA prediction including miRanda, TargetScanS and PicTar. The results showed that retained introns can increase putative miRNA targets in human mRNA. Retained introns have higher chances than other regions of 3' UTR in involving the site of miRNAs targets of most genes which contain putative miRNA targets within it. Furthermore, some transcripts contain miRNA targets solely because of the retained introns in 3' UTR. In addition, we examined those 'Ignored' retained introns by miRanda software and the results indicated that miRNAs may contain many more putative targets.

Signalling through kinase-defective domains: the prevalence of atypical receptor-like kinases in plants

The structure of plant receptor-like kinases (RLKs) is similar to that of animal receptor tyrosine kinases (RTKs), and consists of an extracellular domain, a transmembrane span, and a cytoplasmic domain containing the conserved kinase domain. The mechanism by which animal RTKs, and probably plant RLKs, signal includes the dimerization of the receptor, their intermolecular phosphorylation, and the phosphorylation of downstream signalling proteins. However, atypical RTKs with a kinase-dead domain that signal through phosphorylation-independent mechanisms have also been described in animals. In the last few years, some atypical RLKs have also been reported in plants. Here these examples and their possible signalling mechanisms are reviewed. Plant genomes contain a much larger number of genes coding for receptor kinases than other organisms. The prevalence of atypical RLKs in plants is analysed here. A sequence analysis of the Arabidopsis kinome revealed that 13% of the kinase genes do not retain some of the residues that are considered as invariant within kinase catalytic domains, and are thus putatively kinase-defective. This percentage rises to close to 20% when analysing RLKs, suggesting that phosphorylation-independent mechanisms mediated by atypical RLKs are particularly important for signal transduction in plants.

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.

Regulation of the beta-hydroxyacyl ACP dehydratase gene of Picea mariana by alternative splicing

The gene for beta-hydroxyacyl ACP dehydratase, a de novo fatty acid biosynthetic enzyme, was cloned from Picea mariana (black spruce) and consists of five exons and four introns. The first intron of the beta-hydroxyacyl ACP dehydratase mRNA is alternatively spliced. Retention of intron 1 in splice variants results in truncation of the beta-hydroxyacyl ACP dehydratase ORF at a premature termination codon. In addition, splicing of intron 1 was found to be associated with cold temperature. mRNAs retaining intron 1 increase with seed imbibition at 22 degrees C but not 4 degrees C, whereas, splicing of intron 1 increases in winter weeks with temperatures below freezing. These results provide evidence that alternative splicing may also contribute to regulation of lipid biosynthesis in Picea mariana.

Comprehensive analysis of alternative splicing in rice and comparative analyses with Arabidopsis

Background

Recently, genomic sequencing efforts were finished for Oryza sativa (cultivated rice) and Arabidopsis thaliana (Arabidopsis). Additionally, these two plant species have extensive cDNA and expressed sequence tag (EST) libraries. We employed the Program to Assemble Spliced Alignments (PASA) to identify and analyze alternatively spliced isoforms in both species.

Results

A comprehensive analysis of alternative splicing was performed in rice that started with >1.1 million publicly available spliced ESTs and over 30,000 full length cDNAs in conjunction with the newly enhanced PASA software. A parallel analysis was performed with Arabidopsis to compare and ascertain potential differences between monocots and dicots. Alternative splicing is a widespread phenomenon (observed in greater than 30% of the loci with transcript support) and we have described nine alternative splicing variations. While alternative splicing has the potential to create many RNA isoforms from a single locus, the majority of loci generate only two or three isoforms and transcript support indicates that these isoforms are generally not rare events. For the alternate donor (AD) and acceptor (AA) classes, the distance between the splice sites for the majority of events was found to be less than 50 basepairs (bp). In both species, the most frequent distance between AA is 3 bp, consistent with reports in mammalian systems. Conversely, the most frequent distance between AD is 4 bp in both plant species, as previously observed in mouse. Most alternative splicing variations are localized to the protein coding sequence and are predicted to significantly alter the coding sequence.

Conclusion

Alternative splicing is widespread in both rice and Arabidopsis and these species share many common features. Interestingly, alternative splicing may play a role beyond creating novel combinations of transcripts that expand the proteome. Many isoforms will presumably have negative consequences for protein structure and function, suggesting that their biological role involves post-transcriptional regulation of gene expression.

Different levels of alternative splicing among eukaryotes

Alternative splicing increases transcriptome and proteome diversification. Previous analyses aiming at comparing the rate of alternative splicing between different organisms provided contradicting results. These contradicting results were attributed to the fact that both analyses were dependent on the expressed sequence tag (EST) coverage, which varies greatly between the tested organisms. In this study we compare the level of alternative splicing among eight different organisms. By employing an EST independent approach we reveal that the percentage of genes and exons undergoing alternative splicing is higher in vertebrates compared with invertebrates. We also find that alternative exons of the skipping type are flanked by longer introns compared to constitutive ones, whereas alternative 5′ and 3′ splice sites events are generally not. In addition, although the regulation of alternative splicing and sizes of introns and exons have changed during metazoan evolution, intron retention remained the rarest type of alternative splicing, whereas exon skipping is more prevalent and exhibits a slight increase, from invertebrates to vertebrates. The difference in the level of alternative splicing suggests that alternative splicing may contribute greatly to the mammal higher level of phenotypic complexity, and that accumulation of introns confers an evolutionary advantage as it allows increasing the number of alternative splicing forms.

An intron with a constitutive transport element is retained in a Tap messenger RNA

Alternative splicing is a key factor contributing to genetic diversity and evolution. Intron retention, one form of alternative splicing, is common in plants but rare in higher eukaryotes, because messenger RNAs with retained introns are subject to cellular restriction at the level of cytoplasmic export and expression. Often, retention of internal introns restricts the export of these mRNAs and makes them the targets for degradation by the cellular nonsense-mediated decay machinery if they contain premature stop codons. In fact, many of the database entries for complementary DNAs with retained introns represent them as artefacts that would not affect the proteome. Retroviruses are important model systems in studies of regulation of RNAs with retained introns, because their genomic and mRNAs contain one or more unspliced introns. For example, Mason-Pfizer monkey virus overcomes cellular restrictions by using a cis-acting RNA element known as the constitutive transport element (CTE). The CTE interacts directly with the Tap protein (also known as nuclear RNA export factor 1, encoded by NXF1), which is thought to be a principal export receptor for cellular mRNA, leading to the hypothesis that cellular mRNAs with retained introns use cellular CTE equivalents to overcome restrictions to their expression. Here we show that the Tap gene contains a functional CTE in its alternatively spliced intron 10. Tap mRNA containing this intron is exported to the cytoplasm and is present in polyribosomes. A small Tap protein is encoded by this mRNA and can be detected in human and monkey cells. Our results indicate that Tap regulates expression of its own intron-containing RNA through a CTE-mediated mechanism. Thus, CTEs are likely to be important elements that facilitate efficient expression of mammalian mRNAs with retained introns.

Stress- and development-induced expression of spliced and unspliced transcripts from two highly similar dehydrin 1 genes in V. riparia and V. vinifera

Dehydrins are proteins that accumulate in vegetative tissues subjected to various dehydrating stress conditions such as cold, drought, and salinity and in seeds at later stages of embryogenesis. Here, we report on two highly identical dehydrin genes, DHN1a and DHN1b, in wild and cultivated grapes, Vitis riparia and Vitis vinifera, and their expression in different tissues and under different environmental conditions. The two genes and their transcripts can easily be distinguished by RT-PCR because DHN1b has an 18 bp deletion compared to DHN1a. V. riparia expressed only DHN1a; V. vinifera expressed both DHN1a and DHN1b. Spliced transcripts, DHN1-S, encoding a putative YSK(2)-type dehydrin were present in low amounts in control leaves, but in high amounts in buds and seeds. Unspliced transcripts, DHN1-U, accumulated to high levels in buds and seeds. Cold, drought, and ABA treatment increased accumulation of both DHN1-S and DHN1-U in leaves, whereas short-day treatment increased only DHN1-S. The possible relation of these results with the difference in freezing stress tolerance between V. riparia and V. vinifera is discussed.

BIGPETALp, a bHLH transcription factor is involved in the control of Arabidopsis petal size

In Arabidopsis, APETALA1, PISTILLATA, APETALA3 and SEPALLATA interact to form multimeric protein complexes required to specify petal identity. However, the downstream events that lead to petal specific shape and size remain largely unknown. Organ final size can be influenced by cell number or cell expansion or both. To date, no gene that specifically limits petal size by controlling postmitotic cell expansion has been identified. Here we have identified a novel petal-expressed, basic helix-loop-helix encoding gene (BIGPETAL, BPE) that is involved in the control of petal size. BPE is expressed via two mRNAs derived from an alternative splicing event. The BPEub transcript is expressed ubiquitously, whereas the BPEp transcript is preferentially expressed in petals. We demonstrate that BPEp is positively regulated downstream of APETALA3, PISTILLATA, APETALA1 and PISTILLATA3 and is negatively regulated downstream of AGAMOUS. Plants that lack the petal-expressed variant BPEp have larger petals as a result of increased cell size, showing that BPEp interferes with postmitotic cell expansion. BPEp is therefore a part of the network that links the patterning genes to final morphogenesis.

Identification and characterization of syntaxin 1 antisense variants in Limulus polyphemus

1.A Limulus SMART(TM) cDNA library screening resulted in the cloning of four syntaxin 1 homologs (referred to as Limulus syntaxin [Lim-syn] 1A, 1B, 1C, and 1D) (Wang, Y., Cao, Z., Xu, W., Kemp, M. D., McAdory, B. S., Newkirk, R. F., Ivy, M. T., and Townsel, J. G. (2004). Gene 326:189-199) and two novel intron-retaining syntaxin 1-like variants, designated Limulus syntaxin variant [Lim-synV] 1A/1C and Lim-synV 1B/1D. 2.The variants exhibited high amino acid sequence identity with the four syntaxin 1 homologs. Specifically, Lim-synV 1A/1C and Lim-synV 1B/1D were homologous to Lim-syn 1A/1C and Lim-syn 1B/1D, respectively. Surprisingly, both Lim-synV 1A/1C and 1B/1D are unusual in that each has a poly A+ tail, an intron, and the common splice motif "GT-AG" at the intron-exon boundary. Exons one and two on the complementary transcript of Lim-synV 1B/1D are separated by a 150 bp intron beginning at #95/96 of the predicted sequences for Lim-syn 1B and 1D, respectively. 3. In contrast, examination of the approximately 3.17 kb Lim-synV 1A/1C clone indicated the inclusion of an insert of 1120 base pairs (bp) beginning at codon #37/38 of the predicted Lim-syn 1A and 1C cDNAs' open reading frames (ORFs). Further, the intron sequence of Lim-synV 1A/C contained multiple stop codons and showed no significant homology to other known sequences as determined by a search of the GenBank database. Thus, the focus of this paper will be Lim-synV 1B/D exclusively. 4. To substantiate that an intron is retained in the full-length mRNA, two types of syntaxin cDNA fragments for Lim-syn 1B/D were generated by RT-PCR and analyzed on Northern blots. The products generated were a mixture of intron-retaining, as well as intron-spliced products. The syntaxin-like variants that retained the intron presumably are derived from a mRNA molecule that has not undergone splicing.5. Although the significance of such intron-containing mRNAs in Limulus has not yet been elucidated, future studies of such variants may serve to broaden our knowledge concerning established splicing mechanisms as well as to focus attention on nonconventional concepts about gene product regulation.

Regulation of the kinase activity of the MIK GCK-like MAP4K by alternative splicing

Alternative splicing of introns is essential to ensure the complexity of mammalian genome functions. In particular, the generation of a high number of different isoforms by alternative splicing is an important characteristic of genes coding for signalling proteins such as mitogen activated protein kinases (MAPKs). This is thought to allow these proteins to transduce multiple stimuli in a highly regulated manner. Plant genes are also subjected to alternative splicing. Nevertheless, clear examples of the functional consequences of this phenomenon are still scarce in plants. MIK is a maize gene coding for a GCK-like MAP4K that can be activated by interaction with maize atypical receptor kinase (MARK), an atypical receptor kinase. Here we show that MIK is subjected to alternative splicing. Expression of MIK leads to, at least, 4 different mature mRNAs that accumulate with particular expression profiles during maize development. Our results show that the polypeptides encoded by the different MIK mRNAs display different kinase activity and are differentially activated by interaction with the MARK receptor. Two MIK isoforms display constitutive kinase activity, one isoform is inactive but can be activated by MARK, and the fourth MIK isoform is inactive and cannot be activated by MARK. Our results constitute a clear example of the biochemical consequences of alternative splicing in plants. The selective conservation during evolution of the intron-exon structure of the region coding for the regulator domain of MIK, as well as the maintenance in maize, rice and Arabidopsis of the alternative splicing of some of these introns, are strong indications of its functional importance.

The Wilms' tumor 1 (WT1) gene (+KTS isoform) functions with a CTE to enhance translation from an unspliced RNA with a retained intron

The Wilms' tumor 1 (WT1) gene plays an important role in mammalian urogenital development, and dysregulation of this gene is observed in many human cancers. Alternative splicing of WT1 RNA leads to the expression of two major protein isoforms, WT1(+KTS) and WT1(-KTS). Whereas WT1(-KTS) acts as a transcriptional regulator, no clear function has been ascribed to WT1(+KTS), despite the fact that this protein is crucial for normal development. Here we show that WT1(+KTS) functions to enhance expression from RNA possessing a retained intron and containing either a cellular or viral constitutive transport element (CTE). WT1(+KTS) expression increases the levels of unspliced RNA containing a CTE and specifically promotes the association of this RNA with polyribosomes. These studies provide further support for links between different steps in RNA metabolism and for the existence of post-transcriptional operons.

Genomewide comparative analysis of alternative splicing in plants

Alternative splicing (AS) has been extensively studied in mammalian systems but much less in plants. Here we report AS events deduced from EST/cDNA analysis in two model plants: Arabidopsis and rice. In Arabidopsis, 4,707 (21.8%) of the genes with EST/cDNA evidence show 8,264 AS events. Approximately 56% of these events are intron retention (IntronR), and only 8% are exon skipping. In rice, 6,568 (21.2%) of the expressed genes display 14,542 AS events, of which 53.5% are IntronR and 13.8% are exon skipping. The consistent high frequency of IntronR suggests prevalence of splice site recognition by intron definition in plants. Different AS events within a given gene occur, for the most part, independently. In total, 36-43% of the AS events produce transcripts that would be targets of the non-sense-mediated decay pathway, if that pathway were to operate in plants as in humans. Forty percent of Arabidopsis AS genes are alternatively spliced also in rice, with some examples strongly suggesting a role of the AS event as an evolutionary conserved mechanism of posttranscriptional regulation. We created a comprehensive web-interfaced database to compile and visualize the evidence for alternative splicing in plants (Alternative Splicing in Plants, available at www.plantgdb.org/ASIP).

Intron retention as an alternative splice variant of the rat urocortin 1 gene

Urocortin 1, highly conserved metazoan gene of the corticotropin-releasing hormone family, is a simple gene structured in two exons and the corresponding intron. The urocortin 1 prepropeptide is entirely coded in the second exon. Preliminary non-isotopic in situ hybridization experiments with an oligonucleotide complementary to an intron sequence of the urocortin 1 gene showed a significant cytoplasmic-like staining, suggesting the occurrence of an intron-retained urocortin 1 transcript. This observation prompted us to study whether the urocortin 1 gene presents alternative splicing by intron retention event. Confocal fluorescent in situ hybridization for urocortin 1 RNA and the use of the specific DNA dye TOPRO-3 allowed us to show significant expression of the intron-retained urocortin 1 transcript that did not colocalize with TOPRO-3 staining indicating a cytoplasmic localization for the intron-retained urocortin 1 transcript. The natural occurrence of a polyadenylated intron-retained urocortin 1 RNA was further documented by reverse transcriptase polymerase chain reaction (PCR), primed with oligo(dT), of total RNA extracted from three brain regions, a midbrain region containing the Edinger-Westphal nucleus, cerebellum and prefrontal cortex. In the three brain regions studied, it was possible to amplify both intron-less as well as intron-retained urocortin 1 transcripts. The use of PCR primers that simultaneously amplify both urocortin 1 transcripts allowed us to show that the expression of both urocortin 1 transcripts differs among the brain regions analyzed, suggesting a tissue specific regulation of this alternative splicing. In silico analysis of the five known mammalian urocortin 1 genomic sequences showed high conservation of the urocortin 1 intron sequence. Further studies should investigate the regulation of this intron retention event and its consequence for the functionality of the urocortin 1 gene.

Highly specific gene silencing by artificial microRNAs in Arabidopsis

Plant microRNAs (miRNAs) affect only a small number of targets with high sequence complementarity, while animal miRNAs usually have hundreds of targets with limited complementarity. We used artificial miRNAs (amiRNAs) to determine whether the narrow action spectrum of natural plant miRNAs reflects only intrinsic properties of the plant miRNA machinery or whether it is also due to past selection against natural miRNAs with broader specificity. amiRNAs were designed to target individual genes or groups of endogenous genes. Like natural miRNAs, they had varying numbers of target mismatches. Previously determined parameters of target selection for natural miRNAs could accurately predict direct targets of amiRNAs. The specificity of amiRNAs, as deduced from genome-wide expression profiling, was as high as that of natural plant miRNAs, supporting the notion that extensive base pairing with targets is required for plant miRNA function. amiRNAs make an effective tool for specific gene silencing in plants, especially when several related, but not identical, target genes need to be downregulated. We demonstrate that amiRNAs are also active when expressed under tissue-specific or inducible promoters, with limited nonautonomous effects. The design principles for amiRNAs have been generalized and integrated into a Web-based tool (http://wmd.weigelworld.org).

Whole-Genome Microarray in Arabidopsis Facilitates Global Analysis of Retained Introns

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism that can increase protein diversity and affect mRNA stability. Different types of AS have been observed; these include exon skipping, alternative donor or acceptor site and intron retention. In humans, exon skipping is the most common type while intron retention is rare. In contrast, in Arabidopsis, intron retention is the most prevalent AS type (approximately 40%). Here we show that direct transcript expression analysis using high-density oligonucleotide-based whole-genome microarrays (WGAs) is particularly amenable for assessing global intron retention in Arabidopsis. By applying a novel algorithm retained introns are detected in 8% of the transcripts examined. A sampling of 14 transcripts showed that 86% can be confirmed by RT-PCR. This rate of detection predicts an overall total AS rate of 20% for Arabidopsis compared with 10-22% based on EST/cDNA-based analysis. These findings will facilitate monitoring constitutive and dynamic whole-genome splicing on the next generation WGA slides.

Features of Arabidopsis genes and genome discovered using full-length cDNAs

Arabidopsis is currently the reference genome for higher plants. A new, more detailed statistical analysis of Arabidopsis gene structure is presented including intron and exon lengths, intergenic distances, features of promoters, and variant 5'-ends of mRNAs transcribed from the same transcription unit. We also provide a statistical characterization of Arabidopsis transcripts in terms of their size, UTR lengths, 3'-end cleavage sites, splicing variants, and coding potential. These analyses were facilitated by scrutiny of our collection of sequenced full-length cDNAs and much larger collection of 5'-ESTs, together with another set of full-length cDNAs from Salk/Stanford/Plant Gene Expression Center/RIKEN. Examples of alternative splicing are observed for transcripts from 7% of the genes and many of these genes display multiple spliced isoforms. Most splicing variants lie in non-coding regions of the transcripts. Non-canonical splice sites constitute less than 1% of all splice sites. Genes with fewer than four introns display reduced average mRNA levels. Putative alternative transcription start sites were observed in 30% of highly expressed genes and in more than 50% of the genes with low expression. Transcription start sites correlate remarkably well with a CG skew peak in the DNA sequences. The intergenic distances vary considerably, those where genes are transcribed towards one another being significantly shorter. New transcripts, missing in the current TIGR genome annotation and ESTs that are non-coding, including those antisense to known genes, are derived and cataloged in the Supplementary Material. They identify 148 new loci in the Arabidopsis genome. The conclusions drawn provide a better understanding of the Arabidopsis genome and how the gene transcripts are processed. The results also allow better predictions to be made for, as yet, poorly defined genes and provide a reference for comparisons with other plant genomes whose complete sequences are currently being determined. Some comparisons with rice are included in this paper.

Differential mRNA processing in hematopoietic stem cells

Hematopoietic stem cells (HSCs) maintain tissue homeostasis by rapidly responding to environmental changes. Although this function is well understood, the molecular mechanisms governing this characteristic are largely unknown. We used a sequenced-based strategy to explore the role of both transcriptional and post-transcriptional regulation in HSC biology. We characterized the gene expression differences between HSCs, both quiescent and proliferating, and their differentiated progeny. This analysis revealed a large fraction of sequence tags aligned to intronic sequences, which we showed were derived from unspliced transcripts. A comparison of the biological properties of the observed spliced versus unspliced transcripts in HSCs showed that the unspliced transcripts were enriched in genes involved in DNA binding and RNA processing. In addition, levels of unspliced message decreased in a transcript-specific fashion after HSC activation in vivo. This change in unspliced transcript level coordinated with increases in gene expression of splicing machinery components. Combined, these results suggest that post-transcriptional regulation is important in HSC activation in vivo.

Analysis of the cDNAs of hypothetical genes on Arabidopsis chromosome 2 reveals numerous transcript variants

In the fully sequenced Arabidopsis (Arabidopsis thaliana) genome, many gene models are annotated as "hypothetical protein," whose gene structures are predicted solely by computer algorithms with no support from either expressed sequence matches from Arabidopsis, or nucleic acid or protein homologs from other species. In order to confirm their existence and predicted gene structures, a high-throughput method of rapid amplification of cDNA ends (RACE) was used to obtain their cDNA sequences from 11 cDNA populations. Primers from all of the 797 hypothetical genes on chromosome 2 were designed, and, through 5' and 3' RACE, clones from 506 genes were sequenced and cDNA sequences from 399 target genes were recovered. The cDNA sequences were obtained by assembling their 5' and 3' RACE polymerase chain reaction products. These sequences revealed that (1) the structures of 151 hypothetical genes were different from their predictions; (2) 116 hypothetical genes had alternatively spliced transcripts and 187 genes displayed polyadenylation sites; and (3) there were transcripts arising from both strands, from the strand opposite to that of the prediction and possible dicistronic transcripts. Promoters from five randomly chosen hypothetical genes (At2g02540, At2g31270, At2g33640, At2g35550, and At2g36340) were cloned into report constructs, and their expressions are tissue or development stage specific. Our results indicate at least 50% of hypothetical genes on chromosome 2 are expressed in the cDNA populations with about 38% of the gene structures differing from their predictions. Thus, by using this targeted approach, high-throughput RACE, we revealed numerous transcripts including many uncharacterized variants from these hypothetical genes.

Species-specific variation of alternative splicing and transcriptional initiation in six eukaryotes

The genome-wide detection of alternative splicing and transcriptional initiation (ASTI) was conducted in six eukaryotes (human, mouse, fruit fly, nematode, cress and rice) whose genome sequencing has been completed or nearly completed. Transcriptional isoforms were collected by mapping a batch of full-length cDNA sequences onto the respective cognate genomic sequences. Isoforms mapped on the same gene locus were compared pair-wise, ASTI patterns were segmented into minimal spans, and then the minimal patterns (ASTI units) were classified into unique types, such as the cassette type or the alternative donor site. All these procedures were performed automatically under the same conditions so that the results obtained from different species could be compared directly. The fraction of loci that underwent ASTI of the total mapped loci was the largest for mammals and fruit fly, and the smallest for plants. Exactly the same trend was observed for the number of unique ASTI types found in each species. The observed fractional representations of the ASTI types were similar between evolutionarily close species, such as human and mouse or cress and rice. On the other hand, the relative orders of abundance in individual ASTI type were considerably different between evolutionarily distant species, such as between mammals and plants. In human and mouse, alternative splicing other than the retained introns tended to occur within the protein coding sequence (CDS) regions rather than within the untranslated regions (UTRs), whereas this tendency was obscure in the other four species. In all the species examined, the difference in alternative exon lengths was most likely in multiples of three, and this tendency was most prominent when the alternative exons were embedded within the CDSs. These observations are generally consistent with the idea that higher organisms utilize the ASTI mechanisms more extensively and in a more complicated manner than lower organisms, and that ASTI actively participates in the enhancement of the functional and structural diversity of products generated from a limited number of genes on a genome.

Alternative splicing of the Vupur3 transcript in cowpea produces multiple mRNA species with a single protein product that is present in both plastids and mitochondria

A heterogeneous population of cDNAs (designated Vupur3) encoding phosphoribosylglycinamide formyltransferase (GART; EC 2.1.2.2) was isolated from a cowpea (Vigna unguiculata L. Walp.) nodule library. Three classes of cDNA with the same ORF, but differing in their 3'-UTRs, were identified. Southern analysis and sequencing of genomic DNA confirmed that these differences result from alternative splicing of the primary transcript of a single Vupur3 gene. Alternative splicing does not appear to play a role in the production of soybean (Glycine max Merrill.) pur3 transcripts. The presence of the protein product of the Vupur3 gene, GART, in plastids and mitochondria was confirmed by immunoblotting with antibodies raised against the recombinant protein. The antibodies recognised two proteins with apparent molecular masses of 27 and 27.5 kDa in both mitochondria and plastids. All Vupur3 transcripts have two in-frame start codons that are active in wheatgerm in vitro transcription / translation experiments suggesting a mechanism by which the gene product could be targeted to two organelles. Like other genes encoding enzymes for purine synthesis, Vupur3 is expressed in nodules before nitrogen fixation begins but in contrast to these genes its expression does not increase markedly after nitrogen fixation begins.

Ultraconserved elements in insect genomes: a highly conserved intronic sequence implicated in the control of homothorax mRNA splicing

Recently, we identified a large number of ultraconserved (uc) sequences in noncoding regions of human, mouse, and rat genomes that appear to be essential for vertebrate and amniote ontogeny. Here, we used similar methods to identify ultraconserved genomic regions between the insect species Drosophila melanogaster and Drosophila pseudoobscura, as well as the more distantly related Anopheles gambiae. As with vertebrates, ultraconserved sequences in insects appear to occur primarily in intergenic and intronic sequences, and at intron-exon junctions. The sequences are significantly associated with genes encoding developmental regulators and transcription factors, but are less frequent and are smaller in size than in vertebrates. The longest identical, nongapped orthologous match between the three genomes was found within the homothorax (hth) gene. This sequence spans an internal exon-intron junction, with the majority located within the intron, and is predicted to form a highly stable stem-loop RNA structure. Real-time quantitative PCR analysis of different hth splice isoforms and Northern blotting showed that the conserved element is associated with a high incidence of intron retention in hth pre-mRNA, suggesting that the conserved intronic element is critically important in the post-transcriptional regulation of hth expression in Diptera.