Skip to main content

Main menu

  • Home
  • Content
    • Current Issue
    • Archive
    • Preview Papers
  • About
    • Editorial Board and Staff
    • About the Journal
    • Terms & Privacy
  • More
    • Alerts
    • Contact Us
  • Submit a Manuscript
    • Instructions for Authors
    • Submit a Manuscript
  • Other Publications
    • Plant Physiology
    • The Plant Cell
    • Plant Direct
    • The Arabidopsis Book
    • Teaching Tools in Plant Biology
    • ASPB
    • Plantae

User menu

  • My alerts
  • Log in

Search

  • Advanced search
Plant Cell
  • Other Publications
    • Plant Physiology
    • The Plant Cell
    • Plant Direct
    • The Arabidopsis Book
    • Teaching Tools in Plant Biology
    • ASPB
    • Plantae
  • My alerts
  • Log in
Plant Cell

Advanced Search

  • Home
  • Content
    • Current Issue
    • Archive
    • Preview Papers
  • About
    • Editorial Board and Staff
    • About the Journal
    • Terms & Privacy
  • More
    • Alerts
    • Contact Us
  • Submit a Manuscript
    • Instructions for Authors
    • Submit a Manuscript
  • Follow PlantCell on Twitter
  • Visit PlantCell on Facebook
  • Visit Plantae
Research ArticleResearch Article
You have accessRestricted Access

Global Dissection of Alternative Splicing in Paleopolyploid Soybean

Yanting Shen, Zhengkui Zhou, Zheng Wang, Weiyu Li, Chao Fang, Mian Wu, Yanming Ma, Tengfei Liu, Ling-An Kong, De-Liang Peng, Zhixi Tian
Yanting Shen
aState Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
bUniversity of Chinese Academy of Sciences, Beijing 100039, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhengkui Zhou
aState Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zheng Wang
aState Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Weiyu Li
aState Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chao Fang
aState Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
bUniversity of Chinese Academy of Sciences, Beijing 100039, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Mian Wu
aState Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Yanming Ma
aState Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tengfei Liu
aState Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
bUniversity of Chinese Academy of Sciences, Beijing 100039, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ling-An Kong
cState Key Laboratory for Biology of Plant Diseases and Insect Pest, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
De-Liang Peng
cState Key Laboratory for Biology of Plant Diseases and Insect Pest, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zhixi Tian
aState Key Laboratory of Plant Cell and Chromosome Engineering, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Beijing 100101, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: zxtian@genetics.ac.cn

Published March 2014. DOI: https://doi.org/10.1105/tpc.114.122739

  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading
  • © 2014 American Society of Plant Biologists. All rights reserved.

Abstract

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

INTRODUCTION

Alternative splicing (AS) is an important posttranscriptional regulatory mechanism for modulating gene expression and functional diversity in higher eukaryotes (Lareau et al., 2004), and it is commonly found in both animal and plant species. Previous EST cDNA-based analysis led to the identification of a number of AS events in different species. Along with technological advances, high-throughput sequencing provides huge data sets for deep analysis and greatly increases the ability to explore AS. For example, an in-depth analysis using RNA-seq suggested that 92 to 94% of human genes are subject to AS (Wang et al., 2008), a level 2- to 3-fold higher than the previous estimation (Mironov et al., 1999; Modrek et al., 2001). It has been reported that 61% of Arabidopsis thaliana genes and 21.2 to 33% of rice (Oryza sativa) genes are alternatively spliced (G.J. Zhang et al., 2010; Marquez et al., 2012).

AS may be associated with 15% of genetic diseases in humans (Krawczak et al., 1992) and is involved in a range of functions in plants, such as growth, development, signal transduction, and responses to biotic and abiotic stress (Dinesh-Kumar and Baker, 2000; Jordan et al., 2002; Kazan, 2003; Balasubramanian et al., 2006; Egawa et al., 2006; Iida and Go, 2006). However, not all of the AS events are functional (Modrek and Lee, 2002; Lewis et al., 2003). It was estimated that approximately one-third of AS events could produce premature termination codons (PTCs) (Lewis et al., 2003; Wang and Brendel, 2006), and these isoforms are potentially targeted for degradation by the nonsense-mediated mRNA decay system (Kalyna et al., 2012; Drechsel et al., 2013). In addition to causing PTCs, a large number of AS variations are predicted to alter the coding sequence, as AS variations occur within the protein coding regions (Modrek and Lee, 2002; Campbell et al., 2006).

Relative to the predominant transcript isoform, AS can be mainly classified as exon skipping (ES), intron retention (IR), alternative donor sites (AD), and alternative acceptor sites (AA) (Black, 2003). In animals and yeast, ES is the most common AS form and IR is the least common (Modrek and Lee, 2002; Sultan et al., 2008; Wang et al., 2008), whereas in plants, IR is the primary form and ES only accounts for a small portion of AS (Iida and Go, 2006; Wang and Brendel, 2006; Filichkin et al., 2010; Marquez et al., 2012). The differences in the frequencies of the AS types between plants and animals may be due to differences in their pre-mRNA splicing (Reddy, 2007). Usually, organisms that have small introns use an intron definition splicing mechanism, which preferentially leads to IR, whereas those organisms with large introns use an exon definition mechanism, which results in ES (Nakai and Sakamoto, 1994).

Some AS events are shared among different species (Modrek and Lee, 2003; Sorek et al., 2004; Birzele et al., 2008), indicating their universally conserved functions. Nevertheless, a number of AS events have been identified to be specific for individual species. Studies showed that predominant splice variants are more conserved than rare variants, and species-specific AS events may be important for a gene to gain new functions, although many of them have unknown or even no function (Kan et al., 2002; Modrek and Lee, 2003; Thanaraj et al., 2003; Lareau et al., 2004). Numerous studies have attempted to identify the regulators that affect AS (Modrek and Lee, 2002; Buratti and Baralle, 2004; Chasin, 2007; Luco et al., 2011). Typically, short conserved sequences 10 nucleotides in length located in either the exons or introns that stimulate (exonic/intronic splicing enhancers) or inhibit (exonic/intronic splicing silencers) splice site usage (Long and Caceres, 2009; Han et al., 2010) are thought to be the predominant regulators of AS. Other factors may also influence AS, such as the promoter used to drive transcription (Cramer et al., 1997, 1999; Pagani et al., 2003), chromatin modifications (Batsché et al., 2006; Schor et al., 2009; Luco et al., 2010), changes in transcription or the chromatin environment (Kornblihtt et al., 2004), hormone-responsive elements (Auboeuf et al., 2002), and the recruitment of different transcription factors or coactivators to the promoters (Nogues et al., 2002; Auboeuf et al., 2004a, 2004b). However, the mechanisms of AS are complex, and many questions remain. For instance, what are the underlying factors that influence the variations of AS frequency and AS type among tissues? Thus far, the mechanisms that determine the activation and inhibition of different splice forms remain poorly understood (Blencowe, 2006).

Gene duplication (GD) and subsequent divergence are the major evolutionary mechanisms that contribute to genome evolution and gain of novel gene functions (Flagel and Wendel, 2009; Innan and Kondrashov, 2010). Changes in the AS patterns are important for the evolution of GD pairs. The majority of studies have reported that GD decreases the AS frequency, as the number of AS isoforms per locus in singleton genes is higher than in duplicated genes and small gene families have a larger average number of AS isoforms per locus than the large gene families (Kopelman et al., 2005; Su et al., 2006; Talavera et al., 2007; P.G. Zhang et al., 2010). However, some studies have revealed that duplicated genes or paralogous protein family genes tend to have more alternative isoforms than singletons (Jin et al., 2008; Lin et al., 2008; Roux and Robinson-Rechavi, 2011). Another study showed an inverse correlation among recently duplicated genes but a positive correlation among anciently duplicated genes (Chen et al., 2011). Thus far, the cause of AS divergence between the duplicated genes remains unclear.

Soybean (Glycine max) is one of the most economically important legume species. Previous studies have shown that soybean has undergone at least two polyploidization or whole-genome duplication (WGD) events within the last 60 million years (Shoemaker et al., 2006; Schlueter et al., 2007). As a result, ∼75% of the genes are present in more than one copy in the soybean genome (Schmutz et al., 2010). It is hypothesized that a large number of the duplicated gene pairs are highly divergent (Lin et al., 2010; Roulin et al., 2012). Nevertheless, the landscape of AS and the changes in AS between the genes duplicated by WGD have not been explored in soybean. Here, using a genome-wide analysis of high throughput RNA-seq data from 28 developing tissues, we found that a large number of genes were alternatively spliced in soybean. Global dissection revealed that the changes in gene structure and genomic features were the main factors governing AS variation and the AS decrease in WGD-derived genes.

RESULTS

Overview of RNA-seq Data

To explore AS at a genome-wide level and to assess the factors that influence AS in soybean, we performed high-throughput RNA-seq using the Illumina Highseq-2000 sequencing technology for 28 samples, including cotyledon, root, shoot, leaf, flower, pod, and seed tissues at the different developmental stages of germination, trefoil, flowering, seed development, and plant senescence (Supplemental Figure 1). In total, ∼1811 million reads (∼181 Gb) were obtained, with an average of 64.7 million reads (∼6.47 Gb) per sample (Figure 1A; Supplemental Table 1). After the adaptor sequences were trimmed and the low-quality reads were filtered out, 1442 million high-quality reads, corresponding to 79.64% of the total reads, were aligned to the soybean genome using TopHat2 (Trapnell et al., 2009, 2012). Among the mapped reads, ∼89.37% were uniquely aligned, and the remaining 11% were mapped to multiple loci (Figure 1A; Supplemental Table 1).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Summary of the Read Alignments in This Study.

(A) The table shows the statistics of the aligned reads to the soybean genome.

(B) The aligned reads number and the gene length across each chromosome. The columns represent the read number, and the line represents the gene length.

(C) The distribution of reads based on soybean genome annotation.

The aligned read number for a particular chromosome was highly associated with the total gene length in this chromosome (Figure 1B). An overview of the distribution of the aligned reads in 100-kb contiguous subregions along each of the 20 chromosomes showed that most of the reads were located in the chromosome arm regions and were closely associated with the gene distribution and genetic recombination (GR) rates (Supplemental Figure 2). Further investigations showed that the majority of the mapped reads were in exonic regions, and the coverage was ∼1600× per sample. Additionally, approximately one-fifth of the aligned reads were mapped to intronic regions, and a small number of the reads were mapped to nongene regions (Figure 1C). The aligned reads that were located in exonic and intronic regions were used to detect gene expression and to identify the AS events for each gene.

Large Numbers of AS Genes Are Identified in the Soybean Genome

Cufflinks (Trapnell et al., 2010) was used to assemble the transcripts and to estimate the transcriptional expression. The AS events were identified using ASTALAVISTA (Foissac and Sammeth, 2007). In total, 45,083 annotated genes in the soybean genome were expressed in the 28 samples. In the soybean genome annotation, 8355 genes, which accounted for 18.53% of the expressed genes, were predicted to be subject to AS. In our analysis, a total of 154,469 AS events were identified. The AS events were distributed in 23,764 genes, which accounted for 52.71% of the total expressed genes (Figure 2A). Only 7130 genes were identified as AS genes by both our analysis and the genome annotation. A large number of AS genes (16,634 genes) were newly identified in this study. Moreover, among the 45,083 expressed genes, only 37,677 had two or more exons (referred to as multiexonic genes). Therefore, in this aspect, ∼63% of the multiexonic genes were alternatively spliced in soybean (Figure 2A).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Numbers of AS Genes and GO Analysis of Continuously Alternatively Spliced Genes.

(A) Statistics of AS genes according to genome annotation and detected in this study.

(B) Enriched GO terms among AS genes that were alternatively spliced in all tissues examined (details in Supplemental Figure 6): biological process (GO:0008150); metabolic process (GO:0008152); macromolecule metabolic process (GO:0043170); nitrogen compound metabolic process (GO:0006807); gene expression (GO:0010467); cellular macromolecule metabolic process (GO:0044260); nucleobase, nucleoside, nucleotide, and nucleic acid metabolic process (GO:0006139); RNA metabolic process (GO:0016070); RNA processing (GO:0006396); regulation of macromolecule metabolic process (GO:0060255); RNA splicing (GO:0008380); mRNA metabolic process (GO:0016071); regulation of RNA metabolic process (GO:0051252); mRNA processing (GO:0006397); RNA splicing, via transesterification reactions (GO:0000375); regulation of gene expression (GO:0010468); RNA splicing, via transesterification reactions with bulged adenosine as nucleophile (GO:0000377); regulation of RNA splicing (GO:0043484); regulation of mRNA processing (GO:0050684); nuclear mRNA splicing via spliceosome (GO:0000398); and regulation of nuclear RNA splicing via spliceosome (GO:0048024).

To check the accuracy of our identification, 109 AS genes distributed across the 20 chromosomes (Supplemental Figure 3A) were used for validation. These 109 AS genes contained 116 events, including 43 IR, 31 AA, 17 AD, and 25 ES (Supplemental Figure 3B). Fragments containing the AS events were amplified using RT-PCR (primer information is listed in Supplemental Table 2) and were submitted for sequencing using the Illumina Miseq analyzer. After quality control checks, the raw sequencing reads were aligned to the sequences of the validation genes to check the corresponding AS events, such as IR in Glyma10g31240 (Supplemental Figure 3C), ES in Glyma03g00970 (Supplemental Figure 3D), AD in Glyma03g36910 (Supplemental Figure 3E), and AA in Glyma03g28320 (Supplemental Figure 3F). Of the 116 AS events, 109 (94%) were supported by our validation experiments. Our validation results further confirmed that a large number of AS events occur in soybean.

The number of AS genes varied widely among the individual tissues. There were more AS genes in the meristem or tissues from younger developmental stages, and the number of AS genes decreased with plant development in the same organ (Supplemental Figure 4A). We presumed that this unequal distribution might be associated with the total number of expressed genes in the individual tissues. A correlation analysis using data from 28 individual samples showed that the AS gene number was highly correlated with the expressed gene number (Supplemental Figure 4B). To eliminate the bias effect of the expressed gene number, we then calculated the AS gene ratio (AS gene number/total expressed gene number) for each tissue. Interestingly, we found that the AS gene ratios were also higher in the meristem and younger tissues (Supplemental Figure 4C). Further detection suggested that the AS frequency might be associated with the level of gene expression (Supplemental Figure 4D).

Previous studies had indicated that many AS events would change the coding frame and potentially result in nonsense-mediated mRNA decay (Modrek and Lee, 2002; Lewis et al., 2003; Wang and Brendel, 2006; Kalyna et al., 2012; Drechsel et al., 2013). Our analysis showed that ∼10% of the AS events shifted the transcript frame, creating what we term frameshift transcripts. More intriguingly, the frameshift transcript number was also positively correlated with the level of gene expression (Supplemental Figure 4E). We also characterized the domains in each protein isoform of individual genes by RPS-BLAST using the CDD–45746 PSSMs as database (Marchler-Bauer et al., 2011). Consistent with the study in rice and Arabidopsis (Chen et al., 2011), our results demonstrated that the AS genes had more domains than the non-AS genes (Supplemental Figure 5A). Comparison of the genic features between the AS and non-AS genes revealed that AS genes and their exons were longer (Supplemental Figures 5B and 5D) and that AS genes had more exons (Supplemental Figure 5C), which indicated that the variation in domain number between AS genes and non-AS genes may be related to the gene structure. Further investigation showed that within the AS genes, the domain number was also positively correlated with gene length, exon number, and exon length (Supplemental Table 3).

Among the 23,764 AS genes, 1440 were found to be alternatively spliced in all of the samples in this study. To detect the biological processes in which these continuously alternatively spliced genes might be involved, the Gene Ontology (GO) term enrichment was investigated with agriGO (Du et al., 2010) using all soybean genes as the background. The results (Supplemental Figure 6) showed that these AS genes were highly enriched in processes related to RNA splicing, such as mRNA processing, RNA processing, RNA splicing, and nuclear mRNA splicing (Figure 2B).

Characteristics of Different AS Types

Of the 154,469 AS events that we identified, IR represented 26.47% of the total and was the most abundant type (Figure 3A), followed by AA (24.89%), AD (11.28%), and ES (8.89%). The results were in accord with the observations in Arabidopsis (Wang and Brendel, 2006; Filichkin et al., 2010; Marquez et al., 2012). Of all the AS events, 99,847 were identified through comparisons of transcript isoforms within individual tissues (referred to as within-tissue). The other 54,622 events were identified through comparisons of the transcript isoforms among different tissues (referred to as among-tissue) (Supplemental Figure 7A). These numbers indicate that approximately one-third of the AS events were tissue differential AS events. The ratios of the four AS types were also different when compared within tissues and among tissues. Of the within-tissue AS events, AA was the most abundant, whereas of the among-tissue AS events, IR was the most common. Almost the same numbers of ES and AD events were identified in comparisons among tissues, but more AD events were identified in comparisons within tissues (Supplemental Figures 7B and 7C). The four AS types changed remarkably in the individual tissues. Overall, they had relatively lower numbers in the tissues from older developmental stages than those from younger stages (Supplemental Figure 8).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Statistics of the Different AS Events and AS Distribution in the Soybean Genome.

(A) Statistics of the different AS types.

(B) AS, transcriptome profiling, and genomic feature distribution across the soybean genome. From outside to inside, each circle represents the chromosome structure, gene location, gene number, gene length, GC content, and GR rates, aligned reads number and AS event number.

The usage of splice sites was another important feature for AS. Among the AS events detected in this study, a total of 121,048 GT-AG splice sites (97.17%), 2848 GC-AG splice sites (2.29%), 288 AT-AC splice sites (0.23%), and 390 splice sites of other types (0.31%) were identified. The GT-AG site represented the highest proportion of all of the four AS types, whereas GC-AG was the second most common splice site (Supplemental Figure 9), which is consistent with findings from other species (Thanaraj and Clark, 2001; Churbanov et al., 2008; Filichkin et al., 2010). However, the ratios of splice site usage showed divergence among the different AS types. The GC-AG splice site in AD was more common than the other three types (Supplemental Figure 9). In both ES and AA, the proportion of GT-AG was ∼98%, which was higher than that in IR and AD (∼95%). Another interesting phenomenon was that in the soybean genome, ∼4.13% of the splice sites were AT-AC (Supplemental Figure 9E), and this splice site was used by few AS events (Supplemental Figure 9), indicating that AT-AC might be conserved to inhibit AS.

The event length of the different AS types also showed divergence. For IR, the retained intron length ranged from 51 to 7080 bp, with the largest number of events at 81 bp (Supplemental Figure 10A), which was shorter than that in rice (G.J. Zhang et al., 2010). Our analysis also showed that the most frequent AA length was 3 bp (Supplemental Figure 10C), and the most frequent AD length was 4 bp (Supplemental Figure 10D), which is in agreement with previous findings from other species (Zavolan et al., 2003; Akerman and Mandel-Gutfreund, 2006; Campbell et al., 2006). The skipping length in ES peaked at ∼73 bp, but this peak was not as sharp as the other types of AS events. At the same time, unlike the smooth curves of the other types, the skipping length varied widely (Supplemental Figure 10B).

AS Frequency Is Correlated with Intron Length, Exon Number, Gene Transcriptional Level, GC Content, and GR Rates

To obtain an overview of the landscape of AS and its relationships with other genomic features, the distributions of the AS events and other gene features were pooled in 1-Mb contiguous subregions across each of the 20 chromosomes, as illustrated in Supplemental Figure 2. Overall, AS was mainly located in the arm regions of the chromosome and was associated with gene distribution (Figure 3B).

To identify the factors that influence AS, correlations between the AS distribution and intron length, exon number, transcriptional level (aligned reads), GC content, and GR rate were computed in contiguous 1-Mb windows across the genome (see Methods). The AS distribution was significantly positively correlated with the intron length, exon number, aligned reads, and GR rate and was significantly negatively correlated with GC content (Table 1). Because a large part of the soybean genome is occupied by the pericentromere, where fewer genes are embedded (Schmutz et al., 2010), we also analyzed data sets in which the pericentromeric regions were excluded to eliminate the potential influence of the “centromere effect” (Tian et al., 2009). Similar correlation patterns were observed. Most of the correlation coefficients were slightly decreased, but the correlation between AS event and the transcription level increased (Table 1).

View this table:
  • View inline
  • View popup
Table 1. Correlations of AS Distribution with Genome Features in 1-Mb Contiguous Windows

The above analysis using contiguous windows suggested that gene features, such as intron length, transcriptional level, GC content, exon number, and GR rate, regulate AS in soybeans. To further explore these relationships, we performed correlation analysis using data from individual genes, and similar results were obtained (Supplemental Table 4). As the AS event number per AS gene ranged from 1 to 339, we hypothesized that the AS frequency in the different genes may also be related to these gene features. To further confirm this correlation, we classified all expression genes into three groups: high AS gene (more than five AS events/gene), low AS gene (one to four AS events per gene), and non-AS gene (no AS in the gene). We subsequently compared the characteristics among these three groups. Significant differences were detected between the AS and non-AS genes. Consistent with the above analysis, the AS genes (both the high AS gene group and low AS gene group) had significantly longer introns (Figure 4A), more exons (Figure 4B), higher expression levels (Figure 4C), lower GC content (Figure 4D), and higher GR rates (Figure 4E) than the non-AS genes. Although no significant differences were detected for the expression level (Figure 4C) and GR rates (Figure 4E) between the high AS genes and the low AS genes, significant differences were observed in their intron lengths, exon numbers, and GC content. The high AS genes had longer introns (Figure 4A), more exons (Figure 4B), and lower GC content (Figure 4D) than the low AS genes.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Comparisons of Gene Features with Different AS Events.

Intron length (A), exon number (B), gene expression level (C), GC content (D), and GR rates (E). The letter above each column represents its significance level. The different letters indicate significant differences at the level of P < 0.001. The bars indicate se.

Therefore, the above analysis suggested that the AS frequency was positively correlated with the intron length, exon number, gene transcriptional level, and GR rates and was negatively correlated with the GC content. To determine whether these correlations are common in different species, we performed the correlation analyses using the recently identified AS events in Arabidopsis (Marquez et al., 2012). The same correlation patterns were observed (Supplemental Table 5).

Intron Length, Exon Number, Gene Transcriptional Level, GC Content, and GR Rates Have Different Influences on the Different Types of AS Events

We also determined the correlations between intron length, exon number, transcriptional level, GR rate, and GC content and the numbers of the different types of AS events. The results showed that these four individual AS types showed similar patterns. All of these events were significantly positively correlated with the gene expression level, intron length, and exon number and were significantly negatively correlated with the GC content (Supplemental Table 6), which was similar to what was found for the total AS events (Table 1).

Because the AS events in one AS gene can consist of different AS types, many AS genes were repeatedly used in the above correlation analyses for the individual AS types. To clarify the feature differences of the different AS types, the genes that exhibited only one of the four AS types were chosen and classified into four groups: IR, AA, AD, and ES. The features of these groups were then compared. We found that the IR, AA, and AD types had relatively higher transcriptional levels than the ES, and the IR and AD had a higher GC content than the ES and AA. Additionally, the IR types had higher GR rates than the ES, whereas the ES had longer introns than the AA, AD, and IR; finally, the IR types had more exons than the ES, AA, and AD genes (Table 2). These results not only suggested that different AS types have different characteristics, but also indicated that these features may differentially influence the individual AS types. For instance, an increased number of exons may accelerate the IR frequency. To confirm this hypothesis, we classified the AS genes that contained only one AS type into three groups based on exon number: high exon genes (more than 10), medium exon genes (three to 10), and low exon genes (less than three). Then, the frequencies of each AS type within these groups (AS number of each type divided by the total number for this category) were calculated. The results showed that, along with the increase in the exon number, the proportion of IR also increased, while the proportions of AA, AD, and ES events decreased (Supplemental Figure 11A). Similar computations were conducted for the intron length (Supplemental Figure 11B) and gene expression (Supplemental Figure 11C). Overall, with an increase in the intron length, the proportion of ES increased, and the proportion of IR decreased (Supplemental Figure 11B); with an increase in the gene expression, the proportion of IR increased, and the proportion of ES decreased (Supplemental Figure 11C).

View this table:
  • View inline
  • View popup
Table 2. Comparison of Gene Features among Genes with Different AS Type

Decrease of AS Frequency in WGD-Derived Duplicated Genes Is Related to the Changes in Intron Length, Exon Number, and Transcriptional Level

GD is important for gene evolution to gain novel functions. We classified the gene families into three groups: small gene families (containing two to three members), medium gene families (containing four to nine members), and large gene families (containing more than 10 members). Consistent with previous findings (Kopelman et al., 2005; Su et al., 2006; Talavera et al., 2007; P.G. Zhang et al., 2010), we found that the AS event number per gene was inversely correlated with the size of a gene’s family (Supplemental Figure 12A). Comparisons showed that the smaller gene families had longer introns, more exons, higher transcriptional levels, lower GC content, and higher GR rates than the larger gene families (Supplemental Table 7).

To elucidate the factors that drive the decrease of AS events in duplicated genes, the AS variation between the WGD paired genes and its relationship to the changes in gene features were investigated. Among the WGD gene pairs, 17,078 pairs included genes that both were expressed, and the remaining pairs included genes that both were not expressed or only one was expressed (Supplemental Figure 12B). Among the 17,078 expressed pairs, there were 3936 pairs in which only one gene underwent AS and the other did not and 6397 pairs in which the paired genes had different numbers of AS events; among the remainder, there were 5517 pairs in which both genes did not exhibit AS and 1228 pairs in which both genes had the same number of AS events (Supplemental Figure 12B). To identify the factors that lead to AS decrease in the duplicated genes, the transcriptional level and gene features of the AS genes and the non-AS genes were compared using those 3936 WGD gene pairs. A paired t test showed that the AS genes had significantly longer introns, more exons, higher transcriptional levels, and lower GC content than the non-AS genes (Table 3). To further confirm this hypothesis, comparisons were performed between the high AS genes and low AS genes using the 6397 WGD pairs (Table 3). The t test results also showed that the high AS genes had significantly longer introns, more exons, and higher transcriptional levels than the low AS genes (Table 3). These comparisons indicated that the decreases in intron length, exon number, and transcriptional level may reduce the AS frequency of the WGD genes.

View this table:
  • View inline
  • View popup
Table 3. Feature Comparisons between AS Genes and Non-AS Genes, High-AS Genes, and Low-AS Genes within WGD Pair Genes

The 1228 pairs of duplicated genes in which both genes had the same number of AS events (Supplemental Figure 12B) can be further classified into two groups: 441 pairs with the same number of AS events and same AS types, and 787 pairs with the same number of AS events but different AS types. To determine whether the AS type variations in duplicated genes were associated with the genomic features, we classified 787 pairs with the same number of AS events but different AS types into low-IR and high-IR, low-ES and high-ES, low-AA and high-AA, and low-AD and high-AD. The differences in genomic features between these pairs of genes were compared. t test analysis showed that there was no significant association between the changes in AS type and genomic features (Supplemental Table 8). The insignificant difference may indicate that the AS type variation in duplicated genes is complicated or be caused by the small number of duplicated pairing genes used for this computation.

DISCUSSION

AS Is Common in the Soybean Genome and Shows Variation in Different Tissues

Since Walter Gilbert proposed that it might be possible for a gene to produce different mRNA isoforms through the splicing recombination of different exons in 1978 (Gilbert, 1978), increasing evidence has shown that AS commonly occurs in higher eukaryotes. In Drosophila melanogaster, the Dscam gene has the potential to produce more than 38,000 alternatively spliced variants (Graveley, 2005). In human, it was reported that the ratio of AS genes is as high as 92 to 94% (Wang et al., 2008). Studies have predicted that ∼60% of intron-containing genes are alternatively spliced in plants (Reddy et al., 2013). Here, using high-throughput RNA sequencing and comprehensive analyses, we demonstrated that ∼63% of multiexonic genes (intron containing genes) are alternatively spliced in soybean, and, on average, each AS gene contain six to seven AS events (154,469 AS events/23,764 AS genes). Moreover, the samples used in this study were collected mainly from different developmental stages and did not include samples from extreme circumstances, such as stress conditions. It has been suggested that numerous AS events are induced only by abiotic and biotic stresses (Staiger and Brown, 2013); therefore, AS in soybean might be more prevalent than we detected.

Our analysis demonstrated that some genes were constitutively alternatively spliced in all of the samples that we examined. GO analysis indicated that these constitutive AS genes were highly enriched in pathways related to RNA splicing. This result may be partially supported by the observation that in addition to playing important roles in regulating the AS of other genes, the RNA splicing-related genes are also regulated by AS by themselves or other splicing factors (Saltzman et al., 2011; Thomas et al., 2012).

However, numerous observations have suggested that there are also many AS events in plants that are regulated by tissue-specific cues (Yoshimura et al., 2002). Our study demonstrated that approximately one-third of these events showed differential tissue expression, and the number of AS events and the frequency of AS types varied dramatically in the different tissues. The AS frequency was higher in the rapidly developing tissues, such as shoot meristem or young seeds, than it was in other tissues. Similar phenomena were found in functionally complex tissues with diverse cell types that provide diverse functions in metazoans, such as the brain, testis, and immune system (Modrek et al., 2001; Johnson et al., 2003; Yeo et al., 2004; Watson et al., 2005). It was thought that in addition to introducing proteome diversity, AS is an important factor in regulating gene expression (Reddy et al., 2013). These tissue-specific AS events may function in a coordinated manner in specific pathways or interaction networks as groups of genes that are coregulated at the transcriptional level (Blencowe, 2006).

Changes in Gene Structure Have Significant Effects on AS

The gene structures of plants and animals are quite different in terms of the average size of introns. The average intron length in animals is relatively large; for example, introns can be as long as 5 kb in humans (Sakharkar et al., 2004). In plants, the average intron length is much smaller (Iwata and Gotoh, 2011; Marquez et al., 2012); for instance, the average intron lengths are 539 and 470 bp in soybeans and in rice, respectively. The difference in intron size between plants and animals may be related to their different principles of intron processing (Brown et al., 1986; Hartmuth and Barta, 1986), which in turn results in different proportions of the AS types between plants and animals. In animals, ES is more prevalent than other AS types (Modrek and Lee, 2002; Sultan et al., 2008; Wang et al., 2008), whereas in plants, IR is the most common type (Iida and Go, 2006; Wang and Brendel, 2006; Filichkin et al., 2010; Marquez et al., 2012). Our analysis showed that, similar to other plants, IR was the most abundant form of AS in soybean. Furthermore, our results also suggested that the ratios of the different AS types in plants may be different between dicots and monocots. Although it is common for IR to account for the majority of AS events in plant genomes, the ratios of the different AS types show variations in different plant species. In Arabidopsis, the proportions of AA and AD were higher than that of ES (Wang and Brendel, 2006; Filichkin et al., 2010; Marquez et al., 2012), whereas ES was more common than AA and AD in rice (Wang and Brendel, 2006; G.J. Zhang et al., 2010). In our analysis, the proportions of AA and AD were higher than that of ES, which is consistent with results from Arabidopsis (Wang and Brendel, 2006; Filichkin et al., 2010; Marquez et al., 2012). A study on alternative first exon clusters also indicated that alternative first exons may have evolved independently between rice and Arabidopsis (Chen et al., 2007). The difference between the proportions of each AS type in Arabidopsis, rice, and soybean indicate that the principles of intron processing or AS may have diverged between dicots and monocots.

The acquisition of splice variants is thought to occur concomitantly with an increase in the number of exons (Kopelman et al., 2005). Our analysis suggested that the gene structure could greatly affect the AS frequency. Analyses using different data sets, including contiguous windows across whole chromosomes, contiguous windows only from chromosome arms, and individual genes, all showed that the AS frequency was positively correlated with the intron length and exon number. A further analysis demonstrated that the genes with longer introns and more exons exhibited more AS events. Several studies have found that the AS frequency decreases after GD (Kopelman et al., 2005; Su et al., 2006; Talavera et al., 2007; P.G. Zhang et al., 2010), but other studies have reported inconsistent results (Jin et al., 2008; Lin et al., 2008; Roux and Robinson-Rechavi, 2011). Chen et al. (2011) found that the AS variation among duplicated genes was age dependent: There was an inverse correlation among recently duplicated genes, but a positive correlation among anciently duplicated genes. Our results revealed that the AS event number per gene was inversely correlated with the size of a gene’s family in soybean. Statistical analysis found that AS genes and high AS genes had longer introns and more exons than the non-AS genes and the low AS genes, respectively, indicating that the change in intron length and exon number might be the main factors associated with AS variation.

These results suggested that changes in gene structure, especially intron length and exon number, have significant effects on AS frequency.

The Effect of Gene Transcriptional Level on AS Might Be Related to the Efficiency of Splicing

Our analysis suggested that the gene transcriptional level also affects the occurrence of AS, but not as strongly as changes in gene structure. The AS frequency was positively correlated with the transcriptional level, and AS preferentially occurred in the tissues with higher gene expression levels. Additionally, our results suggested that the frameshift transcripts caused by AS were increased in the tissues with higher gene expression. Although some AS isoforms containing PTCs may be functional, most of these transcripts are thought to be nonsense transcripts (Lewis et al., 2003; Filichkin et al., 2010).

We proposed that the high proportion of AS and frameshift transcripts occurring in the tissues with higher gene expression levels might be due to the low or insufficient amount of spliceosome complexes in these tissues. Functional analyses of pol II and its effects on splicing may partially support our hypothesis. It was presumed that similar to other macromolecular machines, the spliceosome is not perfect and will occasionally make mistakes (Graveley, 2001). Transcription and pre-mRNA splicing are coupled through the interaction of RNA Pol II and the SFs, and transcription has been shown to affect AS (reviewed in Luco et al., 2011). The pol II C4 mutant exhibits a lower elongation rate (Chen et al., 1996), and this mutant can inhibit the skipping of the fibronectin EDI exon and stimulate the resplicing of the Hox gene Ultrabi (de la Mata et al., 2003). Additional proof was collected using an artificial cassette exon that was created by mutating the branch point upstream of the constitutive internal exon of the DYN2 gene in yeast. The skipping of this exon can be partially prevented when this gene is expressed in a yeast mutant carrying a slow pol II or when expressed in the presence of elongation inhibitors (Howe et al., 2003). Moreover, if the SFs are mutated, then additional AS events can occur (reviewed in Staiger and Brown, 2013). These results indicate that the balance between the relative rates of spliceosome formation and pol II processivity are important for AS.

No single mechanism can fully explain the precise regulation of AS (Fox-Walsh and Hertel, 2009). Many factors contribute to the control of AS (Luco et al., 2011; Witten and Ule, 2011), and splicing decisions result from the integration of multiple variables (Shin and Manley, 2004). Here, through the comprehensive investigation of the relationships between AS variation and changes in genomic features at a genome-wide level, our results elucidate how AS is affected by the changes in gene structure and expression level and provide clues for future functional analyses.

METHODS

Plant Growth Conditions and Material Collection

Soybean plants (Glycine max cv Williams82) were grown during the growing season at the experimental station of the Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, in 2011. Twenty-eight tissues, including cotyledon, root, shoot, leaf, flower, pod, and seed at the different developmental stages of germination, trefoil, flowering, seed development, and plant senescence, were collected. Each sample was collected from at least five individual plants, and the mixed bundle was quickly frozen in liquid nitrogen for RNA isolation.

RNA Isolation, cDNA Library Construction, and High-Throughput Sequencing

The total RNA was isolated using the TRIzol reagent (Invitrogen) according to the manufacturer’s instructions. To remove the residual DNA, the extracted RNA was treated with RNase-free DNase I (New England Biolabs) for 30 min at 37°C. cDNA libraries were constructed as previously described (Severin et al., 2010) and sequenced using a HiSeq 2000 analyzer.

Read Alignment to the Reference Soybean Genome and Gene Expression Estimation

Soybean genome sequences and annotation (version 109) were downloaded from http://www.phytozome.net/soybean. After removing the low-quality reads and clipping the adapter sequences, the raw RNA sequence data from each sample were mapped to the reference soybean genome using TopHat2 (Trapnell et al., 2009). As the tissues used in this study were obtained from the reference cultivar, Williams 82, the parameter of no mismatch was used during the mapping. The maximum alignment number for each read was set to 2, and the inner distance between the mating pairs was set to 230 bp. The rest of the parameters were set as default. The gene expression was estimated using Cufflinks (Trapnell et al., 2010) and was normalized using the numbers of reads per kilobase of exon sequence in a gene per million mapped reads (Mortazavi et al., 2008). GO analysis was performed using agriGO (Du et al., 2010).

Putative Transcript Assembly and AS Event Identification

After the reads were aligned, empirical transcript construction for each sample was performed using Cufflinks (Trapnell et al., 2010). The minimum isoform fraction was set to 0.05, the small anchor fraction for spliced reads was set to 0.01, the precursor mRNA fraction was set to 0.15, and the minimum fragments supporting assembled transfrag was set to 10. The minimum and maximum intron lengths were set to 30 and 15,000 bp, respectively. The rest of the parameters were then set as default.

To identify AS events, first, the different assembled transcript isoforms were mapped to the corresponding gene model using Cuffcompare, which is included in the Cufflinks program, as previously described (Trapnell et al., 2010). Then, AS identification was performed using ASTALAVISTA (Foissac and Sammeth, 2007). Four main types of AS events, including IR, ES, AA, and AD, were analyzed as previously described (Sammeth et al., 2008).

Validation of AS Events

A total of 116 AS events distributed across the 20 chromosomes were chosen to validate the accuracy of the AS prediction in this study. The 116 transcript fragments, each of which contained a predicted AS event, were amplified using cDNA as template. The fragment length was between 100 and 500 bp. The primer information is listed in Supplemental Table 2. After amplification, all of the PCR products were mixed together and used for Miseq library construction according to the protocol provided by Illumina. The library was sequenced using a Miseq Benchtop Sequencer. After removing the low quality reads and clipping the adapter sequences, the raw Miseq reads were mapped to the sequences of the corresponding transcript fragment. To determine if the predicted AS events were supported by the Miseq reads, Tablet (Milne et al., 2010) was used to visualize the mapped results.

Statistical Analysis

The local GR rate estimation was performed using MareyMap (Rezvoy et al., 2007) as previously described (Tian et al., 2012). For the statistical analysis at the genome-wide level, each chromosome was split into contiguous 1-Mb windows from the beginning to the end, and the genomic features within each window were calculated as previously described (Tian et al., 2012). The GR rates were obtained for each window and were plotted on the basis of their midpoints. For genes that have more than two gene models in the genome annotation, the longest one was used to calculate the intron length, exon number, and GC content.

The correlations between the investigated features were assessed using Pearson’s correlation with 10,000 bootstrap resamplings as previously described (Zhang and Gaut, 2003). Comparisons of genomic features and expression levels between the non-AS and AS genes and between the high AS and low AS genes of the WGD gene pairs were conducted using an ANOVA t test. Comparisons of different features among genes with different AS events and among different types of AS events were conducted using multivariate analysis with Fisher’s LSD test.

Accession Numbers

RNA-seq raw sequence data for the 28 samples from this article can be found in the NCBI Short Read Archive database under the following accession numbers: SRX474427, SRX474441, SRX474445, SRX474430, SRX474431, SRX474433, SRX474432, SRX474439, SRX474442, SRX474419, SRX474428, SRX474440, SRX474443, SRX474424, SRX474423, SRX474422, SRX474434, SRX474436, SRX474437, SRX474416, SRX474435, SRX474438, SRX474421, SRX474420, SRX474446, SRX474444, SRX474426, and SRX474429.

Supplemental Data

The following materials are available in the online version of this article.

  • Supplemental Figure 1. Tissues Collected from the Different Developmental Stages for RNA-seq.

  • Supplemental Figure 2. Distributions of the Aligned Reads, Genomic Features, and Estimated Local GR Rates along the 20 Soybean Chromosomes.

  • Supplemental Figure 3. Validation of AS.

  • Supplemental Figure 4. AS Gene Distributions and Their Relationships with the Expressed Gene Number and Expression Level in Different Tissues.

  • Supplemental Figure 5. Comparisons between AS and Non-AS Genes.

  • Supplemental Figure 6. Enriched Biological Process GO Terms in Genes That Were Alternatively Spliced in All Tissues.

  • Supplemental Figure 7. AS Events within Tissues and among Different Tissues.

  • Supplemental Figure 8. Numbers of Individual AS Types in Different Tissues.

  • Supplemental Figure 9. Splice Site Usage of the Different Types of AS Events.

  • Supplemental Figure 10. Sequence Length Distribution of the Different Types of AS Events.

  • Supplemental Figure 11. Frequencies of the Individual AS Types in the Different Groups Characterized by Different Features.

  • Supplemental Figure 12. AS in Singleton and Duplicated Genes.

  • Supplemental Table 1. Statistics of Aligned Reads in Each Sample.

  • Supplemental Table 2. Primer Information for AS Validation.

  • Supplemental Table 3. Correlations of Domain Number with Gene Features.

  • Supplemental Table 4. Correlations of AS Number with Other Features for Individual Genes.

  • Supplemental Table 5. Correlations of AS Number with Gene Features in Arabidopsis.

  • Supplemental Table 6. Correlations of AS Number with Gene Features for Individual AS Types.

  • Supplemental Table 7. Comparisons of Different Features between Gene Families with Different Sizes.

  • Supplemental Table 8. Feature Comparisons between Duplicated Genes with the Same Number of AS Events but Different AS Types.

Acknowledgments

We thank the two anonymous reviewers for constructive comments. This work was supported by the National Natural Science Foundation of China (Grants 91131005 and 31222042) and “One-hundred talents” Startup Funds from the Chinese Academy of Sciences. The authors declare no conflict of interest.

AUTHOR CONTRIBUTIONS

Y.S., Z.Z., and Z.T. conceived this project and designed all of the experiments. Y.S., Z.W., C.F., M.W., Y.M., L.A.K., and D.L.P. performed the experiments. Y.S., Z.Z., T.L., W.L., and Z.T. analyzed the data. Y.S., Z.Z., and Z.T. wrote the article.

Footnotes

  • www.plantcell.org/cgi/doi/10.1105/tpc.114.122739

  • The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantcell.org) is: Zhixi Tian (zxtian{at}genetics.ac.cn).

  • ↵1 These authors contributed equally to this work.

  • ↵[W] Online version contains Web-only data.

Glossary

AS
alternative splicing
PTC
premature termination codon
ES
exon skipping
IR
intron retention
AD
alternative donor sites
AA
alternative acceptor sites
GD
gene duplication
WGD
whole-genome duplication
GR
genomic recombination
GO
Gene Ontology
  • Received January 8, 2014.
  • Revised March 5, 2014.
  • Accepted March 12, 2014.
  • Published March 28, 2014.

References

  1. ↵
    1. Akerman M.,
    2. Mandel-Gutfreund Y.
    (2006). Alternative splicing regulation at tandem 3′ splice sites. Nucleic Acids Res. 34: 23–31.
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Auboeuf D.,
    2. Dowhan D.H.,
    3. Kang Y.K.,
    4. Larkin K.,
    5. Lee J.W.,
    6. Berget S.M.,
    7. O’Malley B.W.
    (2004b). Differential recruitment of nuclear receptor coactivators may determine alternative RNA splice site choice in target genes. Proc. Natl. Acad. Sci. USA 101: 2270–2274.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Auboeuf D.,
    2. Dowhan D.H.,
    3. Li X.T.,
    4. Larkin K.,
    5. Ko L.,
    6. Berget S.M.,
    7. O’Malley B.W.
    (2004a). CoAA, a nuclear receptor coactivator protein at the interface of transcriptional coactivation and RNA splicing. Mol. Cell. Biol. 24: 442–453.
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Auboeuf D.,
    2. Hönig A.,
    3. Berget S.M.,
    4. O’Malley B.W.
    (2002). Coordinate regulation of transcription and splicing by steroid receptor coregulators. Science 298: 416–419.
    OpenUrlAbstract/FREE Full Text
  5. ↵
    1. Balasubramanian S.,
    2. Sureshkumar S.,
    3. Lempe J.,
    4. Weigel D.
    (2006). Potent induction of Arabidopsis thaliana flowering by elevated growth temperature. PLoS Genet. 2: e106.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Batsché E.,
    2. Yaniv M.,
    3. Muchardt C.
    (2006). The human SWI/SNF subunit Brm is a regulator of alternative splicing. Nat. Struct. Mol. Biol. 13: 22–29.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Birzele F.,
    2. Csaba G.,
    3. Zimmer R.
    (2008). Alternative splicing and protein structure evolution. Nucleic Acids Res. 36: 550–558.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Black D.L.
    (2003). Mechanisms of alternative pre-messenger RNA splicing. Annu. Rev. Biochem. 72: 291–336.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Blencowe B.J.
    (2006). Alternative splicing: New insights from global analyses. Cell 126: 37–47.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Brown J.W.S.,
    2. Feix G.,
    3. Frendewey D.
    (1986). Accurate in vitro splicing of two pre-mRNA plant introns in a HeLa cell nuclear extract. EMBO J. 5: 2749–2758.
    OpenUrlPubMed
  11. ↵
    1. Buratti E.,
    2. Baralle F.E.
    (2004). Influence of RNA secondary structure on the pre-mRNA splicing process. Mol. Cell. Biol. 24: 10505–10514.
    OpenUrlFREE Full Text
  12. ↵
    1. Campbell M.A.,
    2. Haas B.J.,
    3. Hamilton J.P.,
    4. Mount S.M.,
    5. Buell C.R.
    (2006). Comprehensive analysis of alternative splicing in rice and comparative analyses with Arabidopsis. BMC Genomics 7: 327.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Chasin L.A.
    (2007). Searching for splicing motifs. Adv. Exp. Med. Biol. 623: 85–106.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Chen T.W.,
    2. Wu T.H.,
    3. Ng W.V.,
    4. Lin W.C.
    (2011). Interrogation of alternative splicing events in duplicated genes during evolution. BMC Genomics 12 (suppl. 3): S16.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Chen W.H.,
    2. Lv G.,
    3. Lv C.,
    4. Zeng C.,
    5. Hu S.
    (2007). Systematic analysis of alternative first exons in plant genomes. BMC Plant Biol. 7: 55.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Chen Y.,
    2. Chafin D.,
    3. Price D.H.,
    4. Greenleaf A.L.
    (1996). Drosophila RNA polymerase II mutants that affect transcription elongation. J. Biol. Chem. 271: 5993–5999.
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Churbanov A.,
    2. Winters-Hilt S.,
    3. Koonin E.V.,
    4. Rogozin I.B.
    (2008). Accumulation of GC donor splice signals in mammals. Biol. Direct 3: 30.
    OpenUrlCrossRefPubMed
  18. ↵
    1. Cramer P.,
    2. Cáceres J.F.,
    3. Cazalla D.,
    4. Kadener S.,
    5. Muro A.F.,
    6. Baralle F.E.,
    7. Kornblihtt A.R.
    (1999). Coupling of transcription with alternative splicing: RNA pol II promoters modulate SF2/ASF and 9G8 effects on an exonic splicing enhancer. Mol. Cell 4: 251–258.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Cramer P.,
    2. Pesce C.G.,
    3. Baralle F.E.,
    4. Kornblihtt A.R.
    (1997). Functional association between promoter structure and transcript alternative splicing. Proc. Natl. Acad. Sci. USA 94: 11456–11460.
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. de la Mata M.,
    2. Alonso C.R.,
    3. Kadener S.,
    4. Fededa J.P.,
    5. Blaustein M.,
    6. Pelisch F.,
    7. Cramer P.,
    8. Bentley D.,
    9. Kornblihtt A.R.
    (2003). A slow RNA polymerase II affects alternative splicing in vivo. Mol. Cell 12: 525–532.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Dinesh-Kumar S.P.,
    2. Baker B.J.
    (2000). Alternatively spliced N resistance gene transcripts: Their possible role in tobacco mosaic virus resistance. Proc. Natl. Acad. Sci. USA 97: 1908–1913.
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Du Z.,
    2. Zhou X.,
    3. Ling Y.,
    4. Zhang Z.,
    5. Su Z.
    (2010). agriGO: A GO analysis toolkit for the agricultural community. Nucleic Acids Res. 38: W64–W70.
    OpenUrl
  23. ↵
    1. Drechsel G.,
    2. Kahles A.,
    3. Kesarwani A.K.,
    4. Stauffer E.,
    5. Behr J.,
    6. Drewe P.,
    7. Rätsch G.,
    8. Wachter A.
    (2013). Nonsense-mediated decay of alternative precursor mRNA splicing variants is a major determinant of the Arabidopsis steady state transcriptome. Plant Cell 25: 3726–3742.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Egawa C.,
    2. Kobayashi F.,
    3. Ishibashi M.,
    4. Nakamura T.,
    5. Nakamura C.,
    6. Takumi S.
    (2006). Differential regulation of transcript accumulation and alternative splicing of a DREB2 homolog under abiotic stress conditions in common wheat. Genes Genet. Syst. 81: 77–91.
    OpenUrlCrossRefPubMed
  25. ↵
    1. Filichkin S.A.,
    2. Priest H.D.,
    3. Givan S.A.,
    4. Shen R.,
    5. Bryant D.W.,
    6. Fox S.E.,
    7. Wong W.K.,
    8. Mockler T.C.
    (2010). Genome-wide mapping of alternative splicing in Arabidopsis thaliana. Genome Res. 20: 45–58.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Flagel L.E.,
    2. Wendel J.F.
    (2009). Gene duplication and evolutionary novelty in plants. New Phytol. 183: 557–564.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Foissac S.,
    2. Sammeth M.
    (2007). ASTALAVISTA: Dynamic and flexible analysis of alternative splicing events in custom gene datasets. Nucleic Acids Res. 35: W297–W299.
    OpenUrl
  28. ↵
    1. Fox-Walsh K.L.,
    2. Hertel K.J.
    (2009). Splice-site pairing is an intrinsically high fidelity process. Proc. Natl. Acad. Sci. USA 106: 1766–1771.
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Gilbert W.
    (1978). Why genes in pieces? Nature 271: 501.
    OpenUrlCrossRefPubMed
  30. ↵
    1. Graveley B.R.
    (2001). Alternative splicing: Increasing diversity in the proteomic world. Trends Genet. 17: 100–107.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Graveley B.R.
    (2005). Mutually exclusive splicing of the insect Dscam pre-mRNA directed by competing intronic RNA secondary structures. Cell 123: 65–73.
    OpenUrlCrossRefPubMed
  32. ↵
    1. Han S.P.,
    2. Tang Y.H.,
    3. Smith R.
    (2010). Functional diversity of the hnRNPs: Past, present and perspectives. Biochem. J. 430: 379–392.
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Hartmuth K.,
    2. Barta A.
    (1986). In vitro processing of a plant pre-mRNA in a HeLa cell nuclear extract. Nucleic Acids Res. 14: 7513–7528.
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Howe K.J.,
    2. Kane C.M.,
    3. Ares M. Jr.
    . (2003). Perturbation of transcription elongation influences the fidelity of internal exon inclusion in Saccharomyces cerevisiae. RNA 9: 993–1006.
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Iida K.,
    2. Go M.
    (2006). Survey of conserved alternative splicing events of mRNAs encoding SR proteins in land plants. Mol. Biol. Evol. 23: 1085–1094.
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Innan H.,
    2. Kondrashov F.
    (2010). The evolution of gene duplications: Classifying and distinguishing between models. Nat. Rev. Genet. 11: 97–108.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Iwata H.,
    2. Gotoh O.
    (2011). Comparative analysis of information contents relevant to recognition of introns in many species. BMC Genomics 12: 45.
    OpenUrlCrossRefPubMed
  38. ↵
    1. Jin L.,
    2. Kryukov K.,
    3. Clemente J.C.,
    4. Komiyama T.,
    5. Suzuki Y.,
    6. Imanishi T.,
    7. Ikeo K.,
    8. Gojobori T.
    (2008). The evolutionary relationship between gene duplication and alternative splicing. Gene 427: 19–31.
    OpenUrlCrossRefPubMed
  39. ↵
    1. Johnson J.M.,
    2. Castle J.,
    3. Garrett-Engele P.,
    4. Kan Z.,
    5. Loerch P.M.,
    6. Armour C.D.,
    7. Santos R.,
    8. Schadt E.E.,
    9. Stoughton R.,
    10. Shoemaker D.D.
    (2003). Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays. Science 302: 2141–2144.
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Jordan T.,
    2. Schornack S.,
    3. Lahaye T.
    (2002). Alternative splicing of transcripts encoding Toll-like plant resistance proteins - what’s the functional relevance to innate immunity? Trends Plant Sci. 7: 392–398.
    OpenUrlCrossRefPubMed
  41. ↵
    1. Kalyna M.,
    2. et al
    . (2012). Alternative splicing and nonsense-mediated decay modulate expression of important regulatory genes in Arabidopsis. Nucleic Acids Res. 40: 2454–2469.
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Kan Z.,
    2. States D.,
    3. Gish W.
    (2002). Selecting for functional alternative splices in ESTs. Genome Res. 12: 1837–1845.
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Kazan K.
    (2003). Alternative splicing and proteome diversity in plants: The tip of the iceberg has just emerged. Trends Plant Sci. 8: 468–471.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Kopelman N.M.,
    2. Lancet D.,
    3. Yanai I.
    (2005). Alternative splicing and gene duplication are inversely correlated evolutionary mechanisms. Nat. Genet. 37: 588–589.
    OpenUrlCrossRefPubMed
  45. ↵
    1. Kornblihtt A.R.,
    2. de la Mata M.,
    3. Fededa J.P.,
    4. Munoz M.J.,
    5. Nogues G.
    (2004). Multiple links between transcription and splicing. RNA 10: 1489–1498.
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. Krawczak M.,
    2. Reiss J.,
    3. Cooper D.N.
    (1992). The mutational spectrum of single base-pair substitutions in mRNA splice junctions of human genes: causes and consequences. Hum. Genet. 90: 41–54.
    OpenUrlPubMed
  47. ↵
    1. Lareau L.F.,
    2. Green R.E.,
    3. Bhatnagar R.S.,
    4. Brenner S.E.
    (2004). The evolving roles of alternative splicing. Curr. Opin. Struct. Biol. 14: 273–282.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Lewis B.P.,
    2. Green R.E.,
    3. Brenner S.E.
    (2003). Evidence for the widespread coupling of alternative splicing and nonsense-mediated mRNA decay in humans. Proc. Natl. Acad. Sci. USA 100: 189–192.
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Lin H.,
    2. Ouyang S.,
    3. Egan A.,
    4. Nobuta K.,
    5. Haas B.J.,
    6. Zhu W.,
    7. Gu X.,
    8. Silva J.C.,
    9. Meyers B.C.,
    10. Buell C.R.
    (2008). Characterization of paralogous protein families in rice. BMC Plant Biol. 8: 18.
    OpenUrlCrossRefPubMed
  50. ↵
    1. Lin J.Y.,
    2. Stupar R.M.,
    3. Hans C.,
    4. Hyten D.L.,
    5. Jackson S.A.
    (2010). Structural and functional divergence of a 1-Mb duplicated region in the soybean (Glycine max) genome and comparison to an orthologous region from Phaseolus vulgaris. Plant Cell 22: 2545–2561.
    OpenUrlAbstract/FREE Full Text
  51. ↵
    1. Long J.C.,
    2. Caceres J.F.
    (2009). The SR protein family of splicing factors: Master regulators of gene expression. Biochem. J. 417: 15–27.
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Luco R.F.,
    2. Allo M.,
    3. Schor I.E.,
    4. Kornblihtt A.R.,
    5. Misteli T.
    (2011). Epigenetics in alternative pre-mRNA splicing. Cell 144: 16–26.
    OpenUrlCrossRefPubMed
  53. ↵
    1. Luco R.F.,
    2. Pan Q.,
    3. Tominaga K.,
    4. Blencowe B.J.,
    5. Pereira-Smith O.M.,
    6. Misteli T.
    (2010). Regulation of alternative splicing by histone modifications. Science 327: 996–1000.
    OpenUrlAbstract/FREE Full Text
  54. ↵
    1. Marchler-Bauer A.,
    2. et al
    . (2011). CDD: A Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 39: D225–D229.
    OpenUrlAbstract/FREE Full Text
  55. ↵
    1. Marquez Y.,
    2. Brown J.W.,
    3. Simpson C.,
    4. Barta A.,
    5. Kalyna M.
    (2012). Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis. Genome Res. 22: 1184–1195.
    OpenUrlAbstract/FREE Full Text
  56. ↵
    1. Milne I.,
    2. Bayer M.,
    3. Cardle L.,
    4. Shaw P.,
    5. Stephen G.,
    6. Wright F.,
    7. Marshall D.
    (2010). Tablet—Next generation sequence assembly visualization. Bioinformatics 26: 401–402.
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Mironov A.A.,
    2. Fickett J.W.,
    3. Gelfand M.S.
    (1999). Frequent alternative splicing of human genes. Genome Res. 9: 1288–1293.
    OpenUrlAbstract/FREE Full Text
  58. ↵
    1. Modrek B.,
    2. Lee C.
    (2002). A genomic view of alternative splicing. Nat. Genet. 30: 13–19.
    OpenUrlCrossRefPubMed
  59. ↵
    1. Modrek B.,
    2. Lee C.J.
    (2003). Alternative splicing in the human, mouse and rat genomes is associated with an increased frequency of exon creation and/or loss. Nat. Genet. 34: 177–180.
    OpenUrlCrossRefPubMed
  60. ↵
    1. Modrek B.,
    2. Resch A.,
    3. Grasso C.,
    4. Lee C.
    (2001). Genome-wide detection of alternative splicing in expressed sequences of human genes. Nucleic Acids Res. 29: 2850–2859.
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. Mortazavi A.,
    2. Williams B.A.,
    3. McCue K.,
    4. Schaeffer L.,
    5. Wold B.
    (2008). Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 5: 621–628.
    OpenUrlCrossRefPubMed
  62. ↵
    1. Nakai K.,
    2. Sakamoto H.
    (1994). Construction of a novel database containing aberrant splicing mutations of mammalian genes. Gene 141: 171–177.
    OpenUrlCrossRefPubMed
  63. ↵
    1. Nogues G.,
    2. Kadener S.,
    3. Cramer P.,
    4. Bentley D.,
    5. Kornblihtt A.R.
    (2002). Transcriptional activators differ in their abilities to control alternative splicing. J. Biol. Chem. 277: 43110–43114.
    OpenUrlAbstract/FREE Full Text
  64. ↵
    1. Pagani F.,
    2. Stuani C.,
    3. Zuccato E.,
    4. Kornblihtt A.R.,
    5. Baralle F.E.
    (2003). Promoter architecture modulates CFTR exon 9 skipping. J. Biol. Chem. 278: 1511–1517.
    OpenUrlAbstract/FREE Full Text
  65. ↵
    1. Reddy A.S.
    (2007). Alternative splicing of pre-messenger RNAs in plants in the genomic era. Annu. Rev. Plant Biol. 58: 267–294.
    OpenUrlCrossRefPubMed
  66. ↵
    1. Reddy A.S.,
    2. Marquez Y.,
    3. Kalyna M.,
    4. Barta A.
    (2013). Complexity of the alternative splicing landscape in plants. Plant Cell 25: 3657–3683.
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Rezvoy C.,
    2. Charif D.,
    3. Guéguen L.,
    4. Marais G.A.
    (2007). MareyMap: An R-based tool with graphical interface for estimating recombination rates. Bioinformatics 23: 2188–2189.
    OpenUrlAbstract/FREE Full Text
  68. ↵
    1. Roulin A.,
    2. Auer P.L.,
    3. Libault M.,
    4. Schlueter J.,
    5. Farmer A.,
    6. May G.,
    7. Stacey G.,
    8. Doerge R.W.,
    9. Jackson S.A.
    (2012). The fate of duplicated genes in a polyploid plant genome. Plant J. 73: 143–153.
    OpenUrl
  69. ↵
    1. Roux J.,
    2. Robinson-Rechavi M.
    (2011). Age-dependent gain of alternative splice forms and biased duplication explain the relation between splicing and duplication. Genome Res. 21: 357–363.
    OpenUrlAbstract/FREE Full Text
  70. ↵
    1. Sakharkar M.K.,
    2. Chow V.T.,
    3. Kangueane P.
    (2004). Distributions of exons and introns in the human genome. In Silico Biol. (Gedrukt) 4: 387–393.
    OpenUrlPubMed
  71. ↵
    1. Saltzman A.L.,
    2. Pan Q.,
    3. Blencowe B.J.
    (2011). Regulation of alternative splicing by the core spliceosomal machinery. Genes Dev. 25: 373–384.
    OpenUrlAbstract/FREE Full Text
  72. ↵
    1. Sammeth M.,
    2. Foissac S.,
    3. Guigó R.
    (2008). A general definition and nomenclature for alternative splicing events. PLOS Comput. Biol. 4: e1000147.
    OpenUrlCrossRefPubMed
  73. ↵
    1. Schlueter J.A.,
    2. Lin J.Y.,
    3. Schlueter S.D.,
    4. Vasylenko-Sanders I.F.,
    5. Deshpande S.,
    6. Yi J.,
    7. O’Bleness M.,
    8. Roe B.A.,
    9. Nelson R.T.,
    10. Scheffler B.E.,
    11. Jackson S.A.,
    12. Shoemaker R.C.
    (2007). Gene duplication and paleopolyploidy in soybean and the implications for whole genome sequencing. BMC Genomics 8: 330.
    OpenUrlCrossRefPubMed
  74. ↵
    1. Schmutz J.,
    2. et al
    . (2010). Genome sequence of the palaeopolyploid soybean. Nature 463: 178–183.
    OpenUrlCrossRefPubMed
  75. ↵
    1. Schor I.E.,
    2. Rascovan N.,
    3. Pelisch F.,
    4. Alló M.,
    5. Kornblihtt A.R.
    (2009). Neuronal cell depolarization induces intragenic chromatin modifications affecting NCAM alternative splicing. Proc. Natl. Acad. Sci. USA 106: 4325–4330.
    OpenUrlAbstract/FREE Full Text
  76. ↵
    1. Severin A.J.,
    2. et al
    . (2010). RNA-Seq Atlas of Glycine max: A guide to the soybean transcriptome. BMC Plant Biol. 10: 160.
    OpenUrlCrossRefPubMed
  77. ↵
    1. Shin C.,
    2. Manley J.L.
    (2004). Cell signalling and the control of pre-mRNA splicing. Nat. Rev. Mol. Cell Biol. 5: 727–738.
    OpenUrlCrossRefPubMed
  78. ↵
    1. Shoemaker R.C.,
    2. Schlueter J.,
    3. Doyle J.J.
    (2006). Paleopolyploidy and gene duplication in soybean and other legumes. Curr. Opin. Plant Biol. 9: 104–109.
    OpenUrlCrossRefPubMed
  79. ↵
    1. Sorek R.,
    2. Shamir R.,
    3. Ast G.
    (2004). How prevalent is functional alternative splicing in the human genome? Trends Genet. 20: 68–71.
    OpenUrlCrossRefPubMed
  80. ↵
    1. Staiger D.,
    2. Brown J.W.S.
    (2013). Alternative splicing at the intersection of biological timing, development, and stress responses. Plant Cell 25: 3640–3656.
    OpenUrlAbstract/FREE Full Text
  81. ↵
    1. Su Z.,
    2. Wang J.,
    3. Yu J.,
    4. Huang X.,
    5. Gu X.
    (2006). Evolution of alternative splicing after gene duplication. Genome Res. 16: 182–189.
    OpenUrlAbstract/FREE Full Text
  82. ↵
    1. Sultan M.,
    2. et al
    . (2008). A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science 321: 956–960.
    OpenUrlAbstract/FREE Full Text
  83. ↵
    1. Talavera D.,
    2. Vogel C.,
    3. Orozco M.,
    4. Teichmann S.A.,
    5. de la Cruz X.
    (2007). The (in)dependence of alternative splicing and gene duplication. PLOS Comput. Biol. 3: e33.
    OpenUrlCrossRefPubMed
  84. ↵
    1. Thanaraj T.A.,
    2. Clark F.
    (2001). Human GC-AG alternative intron isoforms with weak donor sites show enhanced consensus at acceptor exon positions. Nucleic Acids Res. 29: 2581–2593.
    OpenUrlAbstract/FREE Full Text
  85. ↵
    1. Thanaraj T.A.,
    2. Clark F.,
    3. Muilu J.
    (2003). Conservation of human alternative splice events in mouse. Nucleic Acids Res. 31: 2544–2552.
    OpenUrlAbstract/FREE Full Text
  86. ↵
    1. Thomas J.,
    2. Palusa S.G.,
    3. Prasad K.V.,
    4. Ali G.S.,
    5. Surabhi G.K.,
    6. Ben-Hur A.,
    7. Abdel-Ghany S.E.,
    8. Reddy A.S.
    (2012). Identification of an intronic splicing regulatory element involved in auto-regulation of alternative splicing of SCL33 pre-mRNA. Plant J.
  87. ↵
    1. Tian Z.,
    2. Rizzon C.,
    3. Du J.,
    4. Zhu L.,
    5. Bennetzen J.L.,
    6. Jackson S.A.,
    7. Gaut B.S.,
    8. Ma J.
    (2009). Do genetic recombination and gene density shape the pattern of DNA elimination in rice long terminal repeat retrotransposons? Genome Res. 19: 2221–2230.
    OpenUrlAbstract/FREE Full Text
  88. ↵
    1. Tian Z.,
    2. Zhao M.,
    3. She M.,
    4. Du J.,
    5. Cannon S.B.,
    6. Liu X.,
    7. Xu X.,
    8. Qi X.,
    9. Li M.W.,
    10. Lam H.M.,
    11. Ma J.
    (2012). Genome-wide characterization of nonreference transposons reveals evolutionary propensities of transposons in soybean. Plant Cell 24: 4422–4436.
    OpenUrlAbstract/FREE Full Text
  89. ↵
    1. Trapnell C.,
    2. Pachter L.,
    3. Salzberg S.L.
    (2009). TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics 25: 1105–1111.
    OpenUrlAbstract/FREE Full Text
  90. ↵
    1. Trapnell C.,
    2. Roberts A.,
    3. Goff L.,
    4. Pertea G.,
    5. Kim D.,
    6. Kelley D.R.,
    7. Pimentel H.,
    8. Salzberg S.L.,
    9. Rinn J.L.,
    10. Pachter L.
    (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7: 562–578.
    OpenUrlCrossRefPubMed
  91. ↵
    1. Trapnell C.,
    2. Williams B.A.,
    3. Pertea G.,
    4. Mortazavi A.,
    5. Kwan G.,
    6. van Baren M.J.,
    7. Salzberg S.L.,
    8. Wold B.J.,
    9. Pachter L.
    (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28: 511–515.
    OpenUrlCrossRefPubMed
  92. ↵
    1. Wang B.B.,
    2. Brendel V.
    (2006). Genomewide comparative analysis of alternative splicing in plants. Proc. Natl. Acad. Sci. USA 103: 7175–7180.
    OpenUrlAbstract/FREE Full Text
  93. ↵
    1. Wang E.T.,
    2. Sandberg R.,
    3. Luo S.,
    4. Khrebtukova I.,
    5. Zhang L.,
    6. Mayr C.,
    7. Kingsmore S.F.,
    8. Schroth G.P.,
    9. Burge C.B.
    (2008). Alternative isoform regulation in human tissue transcriptomes. Nature 456: 470–476.
    OpenUrlCrossRefPubMed
  94. ↵
    1. Watson F.L.,
    2. Püttmann-Holgado R.,
    3. Thomas F.,
    4. Lamar D.L.,
    5. Hughes M.,
    6. Kondo M.,
    7. Rebel V.I.,
    8. Schmucker D.
    (2005). Extensive diversity of Ig-superfamily proteins in the immune system of insects. Science 309: 1874–1878.
    OpenUrlAbstract/FREE Full Text
  95. ↵
    1. Witten J.T.,
    2. Ule J.
    (2011). Understanding splicing regulation through RNA splicing maps. Trends Genet. 27: 89–97.
    OpenUrlCrossRefPubMed
  96. ↵
    1. Yeo G.,
    2. Holste D.,
    3. Kreiman G.,
    4. Burge C.B.
    (2004). Variation in alternative splicing across human tissues. Genome Biol. 5: R74.
    OpenUrlCrossRefPubMed
  97. ↵
    1. Yoshimura K.,
    2. Yabuta Y.,
    3. Ishikawa T.,
    4. Shigeoka S.
    (2002). Identification of a cis element for tissue-specific alternative splicing of chloroplast ascorbate peroxidase pre-mRNA in higher plants. J. Biol. Chem. 277: 40623–40632.
    OpenUrlAbstract/FREE Full Text
  98. ↵
    1. Zavolan M.,
    2. Kondo S.,
    3. Schonbach C.,
    4. Adachi J.,
    5. Hume D.A.,
    6. Hayashizaki Y.,
    7. Gaasterland T.,
    8. Group R.G.,
    9. RIKEN GER Group,
    10. GSL Members
    (2003). Impact of alternative initiation, splicing, and termination on the diversity of the mRNA transcripts encoded by the mouse transcriptome. Genome Res. 13: 1290–1300.
    OpenUrlAbstract/FREE Full Text
  99. ↵
    1. Zhang G.J.,
    2. et al
    . (2010). Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome. Genome Res. 20: 646–654.
    OpenUrlAbstract/FREE Full Text
  100. ↵
    1. Zhang L.,
    2. Gaut B.S.
    (2003). Does recombination shape the distribution and evolution of tandemly arrayed genes (TAGs) in the Arabidopsis thaliana genome? Genome Res. 13: 2533–2540.
    OpenUrlAbstract/FREE Full Text
    1. Zhang P.G.,
    2. Huang S.Z.,
    3. Pin A.L.,
    4. Adams K.L.
    (2010). Extensive divergence in alternative splicing patterns after gene and genome duplication during the evolutionary history of Arabidopsis. Mol. Biol. Evol. 27: 1686–1697.
    OpenUrlAbstract/FREE Full Text
View Abstract
PreviousNext
Back to top

Table of Contents

Print
Download PDF
Email Article

Thank you for your interest in spreading the word on Plant Cell.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Global Dissection of Alternative Splicing in Paleopolyploid Soybean
(Your Name) has sent you a message from Plant Cell
(Your Name) thought you would like to see the Plant Cell web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Global Dissection of Alternative Splicing in Paleopolyploid Soybean
Yanting Shen, Zhengkui Zhou, Zheng Wang, Weiyu Li, Chao Fang, Mian Wu, Yanming Ma, Tengfei Liu, Ling-An Kong, De-Liang Peng, Zhixi Tian
The Plant Cell Mar 2014, 26 (3) 996-1008; DOI: 10.1105/tpc.114.122739

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Global Dissection of Alternative Splicing in Paleopolyploid Soybean
Yanting Shen, Zhengkui Zhou, Zheng Wang, Weiyu Li, Chao Fang, Mian Wu, Yanming Ma, Tengfei Liu, Ling-An Kong, De-Liang Peng, Zhixi Tian
The Plant Cell Mar 2014, 26 (3) 996-1008; DOI: 10.1105/tpc.114.122739
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • METHODS
    • Acknowledgments
    • AUTHOR CONTRIBUTIONS
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF

In this issue

The Plant Cell Online: 26 (3)
The Plant Cell
Vol. 26, Issue 3
Mar 2014
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Front Matter (PDF)
View this article with LENS

More in this TOC Section

  • SAUR17 and SAUR50 Differentially Regulate PP2C-D1 during Apical Hook Development and Cotyledon Opening in Arabidopsis
  • AUTOPHAGY-RELATED14 and Its Associated Phosphatidylinositol 3-Kinase Complex Promote Autophagy in Arabidopsis
  • Grass-Specific EPAD1 Is Essential for Pollen Exine Patterning in Rice
Show more RESEARCH ARTICLES

Similar Articles

Our Content

  • Home
  • Current Issue
  • Plant Cell Preview
  • Archive
  • Teaching Tools in Plant Biology
  • Plant Physiology
  • Plant Direct
  • Plantae
  • ASPB

For Authors

  • Instructions
  • Submit a Manuscript
  • Editorial Board and Staff
  • Policies
  • Recognizing our Authors

For Reviewers

  • Instructions
  • Peer Review Reports
  • Journal Miles
  • Transfer of reviews to Plant Direct
  • Policies

Other Services

  • Permissions
  • Librarian resources
  • Advertise in our journals
  • Alerts
  • RSS Feeds
  • Contact Us

Copyright © 2021 by The American Society of Plant Biologists

Powered by HighWire