Evolutionary Character of Alternative Splicing in Plants

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Temporal regulation of alternative splicing events in rice memory under drought stress

Plant adaptation to drought stress is essential for plant survival and crop yield. Recently, harnessing drought memory, which is induced by repeated stress and recovery cycles, was suggested as a means to improve drought resistance at the transcriptional level. However, the genetic mechanism underlying drought memory is unclear. Here, we carried out a quantitative analysis of alternative splicing (AS) events in rice memory under drought stress, generating 12 transcriptome datasets. Notably, we identified exon skipping (ES) as the predominant AS type (>80%) in differential alternative splicing (DAS) in response to drought stress. Applying our analysis pipeline to investigate DAS events following drought stress in six other plant species revealed variable ES frequencies ranging from 9.94% to 60.70% depending on the species, suggesting that the relative frequency of DAS types in plants is likely to be species-specific. The dinucleotide sequence at AS splice sites in rice following drought stress was preferentially GC-AG and AT-AC. Since U12-type splicing uses the AT-AC site, this suggests that drought stress may increase U12-type splicing, and thus increase ES frequency. We hypothesize that multiple isoforms derived from exon skipping may be induced by drought stress in rice. We also identified 20 transcription factors and three highly connected hub genes with potential roles in drought memory that may be good targets for plant breeding.

Applications of Multi-Omics Technologies for Crop Improvement

Multiple "omics" approaches have emerged as successful technologies for plant systems over the last few decades. Advances in next-generation sequencing (NGS) have paved a way for a new generation of different omics, such as genomics, transcriptomics, and proteomics. However, metabolomics, ionomics, and phenomics have also been well-documented in crop science. Multi-omics approaches with high throughput techniques have played an important role in elucidating growth, senescence, yield, and the responses to biotic and abiotic stress in numerous crops. These omics approaches have been implemented in some important crops including wheat (Triticum aestivum L.), soybean (Glycine max), tomato (Solanum lycopersicum), barley (Hordeum vulgare L.), maize (Zea mays L.), millet (Setaria italica L.), cotton (Gossypium hirsutum L.), Medicago truncatula, and rice (Oryza sativa L.). The integration of functional genomics with other omics highlights the relationships between crop genomes and phenotypes under specific physiological and environmental conditions. The purpose of this review is to dissect the role and integration of multi-omics technologies for crop breeding science. We highlight the applications of various omics approaches, such as genomics, transcriptomics, proteomics, metabolomics, phenomics, and ionomics, and the implementation of robust methods to improve crop genetics and breeding science. Potential challenges that confront the integration of multi-omics with regard to the functional analysis of genes and their networks as well as the development of potential traits for crop improvement are discussed. The panomics platform allows for the integration of complex omics to construct models that can be used to predict complex traits. Systems biology integration with multi-omics datasets can enhance our understanding of molecular regulator networks for crop improvement. In this context, we suggest the integration of entire omics by employing the "phenotype to genotype" and "genotype to phenotype" concept. Hence, top-down (phenotype to genotype) and bottom-up (genotype to phenotype) model through integration of multi-omics with systems biology may be beneficial for crop breeding improvement under conditions of environmental stresses.

Changes in Alternative Splicing in Response to Domestication and Polyploidization in Wheat

Alternative splicing (AS) occurs extensively in eukaryotes as an important mechanism for regulating transcriptome complexity and proteome diversity, but variation in the AS landscape in response to domestication and polyploidization in crops is unclear. Hexaploid wheat (AABBDD, Triticum aestivum) has undergone two separate allopolyploidization events, providing an ideal model for studying AS changes during domestication and polyploidization events. In this study, we performed high-throughput transcriptome sequencing of roots and leaves from wheat species with varied ploidies, including wild diploids (A^bA^b, Triticum boeoticum) and tetraploids (AABB, Triticum dicoccoides), domesticated diploids (A^mA^m, Triticum monococcum) and tetraploids (AABB, Triticum dicoccum), hexaploid wheat (AABBDD, T aestivum), as well as newly synthesized hexaploids together with their parents. Approximately 22.1% of genes exhibited AS, with the major AS type being intron retention. The number of AS events decreased after domestication in both diploids and tetraploids. Moreover, the frequency of AS occurrence tended to decrease after polyploidization, consistent with the functional sharing model that proposes AS and duplicated genes are complementary in regulating transcriptome plasticity in polyploid crops. In addition, the subgenomes exhibited biased AS responses to polyploidization, and ∼87.1% of homeologs showed AS partitioning in hexaploid wheat. Interestingly, substitution of the D-subgenome modified 42.8% of AS patterns of the A- and B-subgenomes, indicating subgenome interplay reprograms AS profiles at a genome-wide level, although the causal-consequence relationship requires further study. Conclusively, our study shows that AS variation occurs extensively after polyploidization and domestication in wheat species.

Distinct Evolutionary Origins of Intron Retention Splicing Events in NHX1 Antiporter Transcripts Relate to Sequence Specific Distinctions in Oryza Species

The genome of Asian cultivated rice (Oryza sativa L.) shows the presence of six organelle-specific and one plasma membrane (OsNHX1-7) NHX-type cation proton antiporters. Of these, vacuolar-localized OsNHX1 is extensively characterized. The genus Oryza consists of 27 species and 11 genome-types, with cultivated rice, diploid O. sativa, having an AA-type genome. Oryza NHX1 orthologous regions (gene organization, 5' upstream cis elements, amino acid residues/motifs) from closely related Oryza AA genomes cluster distinctly from NHX1 regions from more ancestral Oryza BB, FF and KKLL genomes. These sequence-specific distinctions also extend to two separate intron retention (IR) events involving Oryza NHX1 transcripts that occur at the 5' and 3' ends of the NHX1 transcripts. We demonstrate that the IR event involving the 5' UTR is present only in more recently evolved Oryza AA genomes while the IR event governing retention of the 13th intron of Oryza NHX1 (terminal intron) is more ancient in origin, also occurring in halophytic wild rice, Oryza coarctata (KKLL). We also report presence of a retro-copy of the OcNHX1 cDNA in the genome of O. coarctata (rOcNHX1). Preferential species and tissue specific up- or down-regulation of the correctly spliced NHX1 transcript/5' UTR/13th intron-retaining splice variants under salinity was observed. The implications of IR on NHX1 mRNA stability and ORF diversity in Oryza spp. is discussed.

Comparative analysis of the accelerated aged seed transcriptome profiles of two maize chromosome segment substitution lines

Seed longevity is one of the most essential characteristics of seed quality. Two chromosome segment substitution lines, I178 and X178, which show significant differences in seed longevity, were subjected to transcriptome sequencing before and after five days of accelerated aging (AA) treatments. Compared to the non-aging treatment, 286 and 220 differentially expressed genes (DEGs) were identified after 5 days of aging treatment in I178 and X178, respectively. Of these DEGs, 98 were detected in both I178 and X178, which were enriched in Gene Ontology (GO) terms of the cellular component of the nuclear part, intracellular part, organelle and membrane. Only 86 commonly downregulated genes were enriched in GO terms of the carbohydrate derivative catabolic process. Additionally, transcriptome analysis of alternative splicing (AS) events in I178 and X178 showed that 63.6% of transcript isoforms occurred AS in all samples, and only 1.6% of transcript isoforms contained 169 genes that exhibited aging-specific AS arising after aging treatment. Combined with the reported QTL mapping result, 7 DEGs exhibited AS after aging treatment, and 13 DEGs in mapping interval were potential candidates that were directly or indirectly related to seed longevity.

Hybrid Sequencing of Full-Length cDNA Transcripts of the Medicinal Plant Scutellaria baicalensis

Scutellaria baicalensis is a well-known medicinal plant that produces biologically active flavonoids, such as baicalin, baicalein, and wogonin. Pharmacological studies have shown that these compounds have anti-inflammatory, anti-bacterial, and anti-cancer activities. Therefore, it is of great significance to investigate the genetic information of S. baicalensis, particularly the genes related to the biosynthetic pathways of these compounds. Here, we constructed the full-length transcriptome of S. baicalensis using a hybrid sequencing strategy and acquired 338,136 full-length sequences, accounting for 93.3% of the total reads. After the removal of redundancy and correction with Illumina short reads, 75,785 nonredundant transcripts were generated, among which approximately 98% were annotated with significant hits in the protein databases, and 11,135 sequences were classified as lncRNAs. Differentially expressed gene (DEG) analysis showed that most of the genes related to flavonoid biosynthesis were highly expressed in the roots, consistent with previous reports that the flavonoids were mainly synthesized and accumulated in the roots of S. baicalensis. By constructing unique transcription models, a total of 44,071 alternative splicing (AS) events were identified, with intron retention (IR) accounting for the highest proportion (44.5%). A total of 94 AS events were present in five key genes related to flavonoid biosynthesis, suggesting that AS may play important roles in the regulation of flavonoid biosynthesis in S. baicalensis. This study provided a large number of highly accurate full-length transcripts, which represents a valuable genetic resource for further research of the molecular biology of S. baicalensis, such as the development, breeding, and biosynthesis of active ingredients.

RNA-seq reveals the involvement of key genes for aerobic adaptation in rice

Adaptation of rice to the aerobic condition is needed to cope with the water scarcity as well as to ensure sustainable yield in future. To understand the molecular mechanisms responsible for aerobic adaptation in rice, we performed RNA-seq analysis of root and shoot i.e. developing panicle tissues at panicle initiation stage in two cultivars adapted to aerobic (CR Dhan 202) and traditional transplanted anaerobic (BPT 5204) conditions. The RNA-seq data emanated from 1.65 billion clean reads with approximately 37 million reads per sample. The number of differentially expressed transcripts was higher in the root than that in the shoot under both aerobic and anaerobic conditions. The transcription factors viz. MADS4, MADS5, MADS6, MADS7, MADS15 and transporters involved in sugar (SWEET3A) and nutrient uptake (PHT1;6, MDR-like ABC and vacuolar iron transporter homolog 2) were highly and uniquely expressed in the aerobic adapted cultivar (AAC) CR Dhan 202 under aerobic condition indicating their role in adaptation. The hormones such as ethylene and abscisic acid might be significantly involved in imparting aerobic adaptation. The higher expression of root related genes in the AAC under aerobic conditions suggests the involvement and sensitivity of roots to the water limiting condition. The metabolic activities are also more pronounced in the roots which impart rigorous plant establishment under the aerobic condition. The presence of alternative splice variants in the transcripts viz. Tetratrico peptide repeat (TPR) domain containing protein and GOLDEN2-LIKE1 (GLK1) additionally confirms that post transcriptional regulation is also crucial for aerobic adaptation. The QTLs related to root traits and stress tolerance harboring the uniquely expressed genes, which were identified in the present study can be deployed in molecular breeding programs to develop elite, high yielding aerobic rice cultivars.

Iso-Seq Allows Genome-Independent Transcriptome Profiling of Grape Berry Development

Transcriptomics has been widely applied to study grape berry development. With few exceptions, transcriptomic studies in grape are performed using the available genome sequence, PN40024, as reference. However, differences in gene content among grape accessions, which contribute to phenotypic differences among cultivars, suggest that a single reference genome does not represent the species' entire gene space. Though whole genome assembly and annotation can reveal the relatively unique or "private" gene space of any particular cultivar, transcriptome reconstruction is a more rapid, less costly, and less computationally intensive strategy to accomplish the same goal. In this study, we used single molecule-real time sequencing (SMRT) to sequence full-length cDNA (Iso-Seq) and reconstruct the transcriptome of Cabernet Sauvignon berries during berry ripening. In addition, short reads from ripening berries were used to error-correct low-expression isoforms and to profile isoform expression. By comparing the annotated gene space of Cabernet Sauvignon to other grape cultivars, we demonstrate that the transcriptome reference built with Iso-Seq data represents most of the expressed genes in the grape berries and includes 1,501 cultivar-specific genes. Iso-Seq produced transcriptome profiles similar to those obtained after mapping on a complete genome reference. Together, these results justify the application of Iso-Seq to identify cultivar-specific genes and build a comprehensive reference for transcriptional profiling that circumvents the necessity of a genome reference with its associated costs and computational weight.

Expanding Alternative Splicing Identification by Integrating Multiple Sources of Transcription Data in Tomato

Tomato (Solanum lycopersicum) is an important vegetable and fruit crop. Its genome was completely sequenced and there are also a large amount of available expressed sequence tags (ESTs) and short reads generated by RNA sequencing (RNA-seq) technologies. Mapping transcripts including mRNA sequences, ESTs, and RNA-seq reads to the genome allows identifying pre-mRNA alternative splicing (AS), a post-transcriptional process generating two or more RNA isoforms from one pre-mRNA transcript. We comprehensively analyzed the AS landscape in tomato by integrating genome mapping information of all available mRNA and ESTs with mapping information of RNA-seq reads which were collected from 27 published projects. A total of 369,911 AS events were identified from 34,419 genomic loci involving 161,913 transcripts. Within the basic AS events, intron retention is the prevalent type (18.9%), followed by alternative acceptor site (12.9%) and alternative donor site (7.3%), with exon skipping as the least type (6.0%). Complex AS types having two or more basic event accounted for 54.9% of total AS events. Within 35,768 annotated protein-coding gene models, 23,233 gene models were found having pre-mRNAs generating AS isoform transcripts. Thus the estimated AS rate was 65.0% in tomato. The list of identified AS genes with their corresponding transcript isoforms serves as a catalog for further detailed examination of gene functions in tomato biology. The post-transcriptional information is also expected to be useful in improving the predicted gene models in tomato. The sequence and annotation information can be accessed at plant alternative splicing database (http://proteomics.ysu.edu/altsplice).

Origin of exon skipping-rich transcriptomes in animals driven by evolution of gene architecture

Background

Alternative splicing, particularly through intron retention and exon skipping, is a major layer of pre-translational regulation in eukaryotes. While intron retention is believed to be the most prevalent mode across non-animal eukaryotes, animals have unusually high rates of exon skipping. However, when and how this high prevalence of exon skipping evolved is unknown. Since exon skipping can greatly expand proteomes, answering these questions sheds light on the evolution of higher organismal complexity in metazoans.

Results

We used RNA-seq data to quantify exon skipping and intron retention frequencies across 65 eukaryotic species, with particular focus on early branching animals and unicellular holozoans. We found that only bilaterians have significantly increased their exon skipping frequencies compared to all other eukaryotic groups. Unlike in other eukaryotes, however, exon skipping in nearly all animals, including non-bilaterians, is strongly enriched for frame-preserving sequences, suggesting that exon skipping involvement in proteome expansion predated the increase in frequency. We also identified architectural features consistently associated with higher exon skipping rates within all studied eukaryotic genomes. Remarkably, these architectures became more prevalent during animal evolution, indicating co-evolution between genome architectures and exon skipping frequencies.

Conclusion

We suggest that the increase of exon skipping rates in animals followed a two-step process. First, exon skipping in early animals became enriched for frame-preserving events. Second, bilaterian ancestors dramatically increased their exon skipping frequencies, likely driven by the interplay between a shift in their genome architectures towards more exon definition and recruitment of frame-preserving exon skipping events to functionally diversify their cell-specific proteomes.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1499-9) contains supplementary material, which is available to authorized users.

The comparison of alternative splicing among the multiple tissues in cucumber

BACKGROUND

Alternative splicing (AS) is an important post-transcriptional process. It has been suggested that most AS events are subject to tissue-specific regulation. However, the global dynamics of AS in different tissues are poorly explored.

RESULTS

To analyse global changes in AS in multiple tissues, we identified the AS events and constructed a comprehensive catalogue of AS events within each tissue based on the genome-wide RNA-seq reads from ten tissues in cucumber. First, we found that 58% of the multi-exon genes underwent AS. We further obtained 565 genes with significantly more AS events compared with random genes. These genes were found significant enrichment in biological processes related to the regulation of actin filament length. Second, significantly different AS event profiles among ten tissues were found. The tissues with the same origin of development are more likely to have a relatively similar AS profile. Moreover, 7370 genes showed tissue-specific AS events and were highly enriched in biological processes related to the positive regulation of cellular component organization. Root-specificity AS genes were related to the cellular response to DNA damage stimulus. Third, the genes with different intron retention (IR) patterns among the ten tissues showed significant difference in GC percentages of the retained intron, and the number of exons and FPKM of the major transcripts.

CONCLUSIONS

Our study provided a comprehensive view of AS in multiple tissues. We revealed novel insights into the patterns of AS in multiple tissues and the tissue-specific AS in cucumber.

Alternative splicing and the evolution of phenotypic novelty

Alternative splicing, a mechanism of post-transcriptional RNA processing whereby a single gene can encode multiple distinct transcripts, has been proposed to underlie morphological innovations in multicellular organisms. Genes with developmental functions are enriched for alternative splicing events, suggestive of a contribution of alternative splicing to developmental programmes. The role of alternative splicing as a source of transcript diversification has previously been compared to that of gene duplication, with the relationship between the two extensively explored. Alternative splicing is reduced following gene duplication with the retention of duplicate copies higher for genes which were alternatively spliced prior to duplication. Furthermore, and unlike the case for overall gene number, the proportion of alternatively spliced genes has also increased in line with the evolutionary diversification of cell types, suggesting alternative splicing may contribute to the complexity of developmental programmes. Together these observations suggest a prominent role for alternative splicing as a source of functional innovation. However, it is unknown whether the proliferation of alternative splicing events indeed reflects a functional expansion of the transcriptome or instead results from weaker selection acting on larger species, which tend to have a higher number of cell types and lower population sizes.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Isoform Evolution in Primates through Independent Combination of Alternative RNA Processing Events

Recent RNA-seq technology revealed thousands of splicing events that are under rapid evolution in primates, whereas the reliability of these events, as well as their combination on the isoform level, have not been adequately addressed due to its limited sequencing length. Here, we performed comparative transcriptome analyses in human and rhesus macaque cerebellum using single molecule long-read sequencing (Iso-seq) and matched RNA-seq. Besides 359 million RNA-seq reads, 4,165,527 Iso-seq reads were generated with a mean length of 14,875 bp, covering 11,466 human genes, and 10,159 macaque genes. With Iso-seq data, we substantially expanded the repertoire of alternative RNA processing events in primates, and found that intron retention and alternative polyadenylation are surprisingly more prevalent in primates than previously estimated. We then investigated the combinatorial mode of these alternative events at the whole-transcript level, and found that the combination of these events is largely independent along the transcript, leading to thousands of novel isoforms missed by current annotations. Notably, these novel isoforms are selectively constrained in general, and 1,119 isoforms have even higher expression than the previously annotated major isoforms in human, indicating that the complexity of the human transcriptome is still significantly underestimated. Comparative transcriptome analysis further revealed 502 genes encoding selectively constrained, lineage-specific isoforms in human but not in rhesus macaque, linking them to some lineage-specific functions. Overall, we propose that the independent combination of alternative RNA processing events has contributed to complex isoform evolution in primates, which provides a new foundation for the study of phenotypic difference among primates.

50/50 Expressional Odds of Retention Signifies the Distinction between Retained Introns and Constitutively Spliced Introns in Arabidopsis thaliana

Intron retention, one of the most prevalent alternative splicing events in plants, can lead to introns retained in mature mRNAs. However, in comparison with constitutively spliced introns (CSIs), the relevantly distinguishable features for retained introns (RIs) are still poorly understood. This work proposes a computational pipeline to discover novel RIs from multiple next-generation RNA sequencing (RNA-Seq) datasets of Arabidopsis thaliana. Using this pipeline, we detected 3,472 novel RIs from 18 RNA-Seq datasets and re-confirmed 1,384 RIs which are currently annotated in the TAIR10 database. We also use the expression of intron-containing isoforms as a new feature in addition to the conventional features. Based on these features, RIs are highly distinguishable from CSIs by machine learning methods, especially when the expressional odds of retention (i.e., the expression ratio of the RI-containing isoforms relative to the isoforms without RIs for the same gene) reaches to or larger than 50/50. In this case, the RIs and CSIs can be clearly separated by the Random Forest with an outstanding performance of 0.95 on AUC (the area under a receiver operating characteristics curve). The closely related characteristics to the RIs include the low strength of splice sites, high similarity with the flanking exon sequences, low occurrence percentage of YTRAY near the acceptor site, existence of putative intronic splicing silencers (ISSs, i.e., AG/GA-rich motifs) and intronic splicing enhancers (ISEs, i.e., TTTT-containing motifs), and enrichment of Serine/Arginine-Rich (SR) proteins and heterogeneous nuclear ribonucleoparticle proteins (hnRNPs).

Identification and analysis of alternative splicing events in Phaseolus vulgaris and Glycine max

BACKGROUND

The vast diversification of proteins in eukaryotic cells has been related with multiple transcript isoforms from a single gene that result in alternative splicing (AS) of primary transcripts. Analysis of RNA sequencing data from expressed sequence tags and next generation RNA sequencing has been crucial for AS identification and genome-wide AS studies. For the identification of AS events from the related legume species Phaseolus vulgaris and Glycine max, 157 and 88 publicly available RNA-seq libraries, respectively, were analyzed.

RESULTS

We identified 85,570 AS events from P. vulgaris in 72% of expressed genes and 134,316 AS events in 70% of expressed genes from G. max. These were categorized in seven AS event types with intron retention being the most abundant followed by alternative acceptor and alternative donor, representing ~75% of all AS events in both plants. Conservation of AS events in homologous genes between the two species was analyzed where an overrepresentation of AS affecting 5'UTR regions was observed for certain types of AS events. The conservation of AS events was experimentally validated for 8 selected genes, through RT-PCR analysis. The different types of AS events also varied by relative position in the genes. The results were consistent in both species.

CONCLUSIONS

The identification and analysis of AS events are first steps to understand their biological relevance. The results presented here from two related legume species reveal high conservation, over ~15-20 MY of divergence, and may point to the biological relevance of AS.

An improved assembly and annotation of the allohexaploid wheat genome identifies complete families of agronomic genes and provides genomic evidence for chromosomal translocations

Advances in genome sequencing and assembly technologies are generating many high-quality genome sequences, but assemblies of large, repeat-rich polyploid genomes, such as that of bread wheat, remain fragmented and incomplete. We have generated a new wheat whole-genome shotgun sequence assembly using a combination of optimized data types and an assembly algorithm designed to deal with large and complex genomes. The new assembly represents >78% of the genome with a scaffold N50 of 88.8 kb that has a high fidelity to the input data. Our new annotation combines strand-specific Illumina RNA-seq and Pacific Biosciences (PacBio) full-length cDNAs to identify 104,091 high-confidence protein-coding genes and 10,156 noncoding RNA genes. We confirmed three known and identified one novel genome rearrangements. Our approach enables the rapid and scalable assembly of wheat genomes, the identification of structural variants, and the definition of complete gene models, all powerful resources for trait analysis and breeding of this key global crop.

New insights into the operative network of FaEO, an enone oxidoreductase from Fragaria x ananassa Duch

The 2-methylene-furan-3-one reductase or Fragaria x ananassa Enone Oxidoreductase (FaEO) catalyses the last reductive step in the biosynthesis of 4-hydroxy-2,5-dimethyl-3(2H)-furanone, a major component in the characteristic flavour of strawberries. In the present work, we describe the association between FaEO and the vacuolar membrane of strawberry fruits. Even if FaEO lacks epitopes for stable or transient membrane-interactions, it contains a calmodulin-binding region, suggesting that in vivo FaEO may be associated with the membrane via a peripheral protein complex with calmodulin. Moreover, we also found that FaEO occurs in dimeric form in vivo and, as frequently observed for calmodulin-regulated proteins, it may be expressed in different isoforms by alternative gene splicing. Further mass spectrometry analysis confirmed that the isolated FaEO consists in the already known isoform and that it is the most characteristic during ripening. Finally, a characterization by absorption spectroscopy showed that FaEO has specific flavoprotein features. The relevance of these findings and their possible physiological implications are discussed.

The Evolutionary Relationship between Alternative Splicing and Gene Duplication

The protein diversity that exists today has resulted from various evolutionary processes. It is well known that gene duplication (GD) along with the accumulation of mutations are responsible, among other factors, for an increase in the number of different proteins. The gene structure in eukaryotes requires the removal of non-coding sequences, introns, to produce mature mRNAs. This process, known as cis-splicing, referred to here as splicing, is regulated by several factors which can lead to numerous splicing arrangements, commonly designated as alternative splicing (AS). AS, producing several transcripts isoforms form a single gene, also increases the protein diversity. However, the evolution and manner for increasing protein variation differs between AS and GD. An important question is how are patterns of AS affected after a GD event. Here, we review the current knowledge of AS and GD, focusing on their evolutionary relationship. These two processes are now considered the main contributors to the increasing protein diversity and therefore their relationship is a relevant, yet understudied, area of evolutionary study.

Plant Proteins Are Smaller Because They Are Encoded by Fewer Exons than Animal Proteins

Protein size is an important biochemical feature since longer proteins can harbor more domains and therefore can display more biological functionalities than shorter proteins. We found remarkable differences in protein length, exon structure, and domain count among different phylogenetic lineages. While eukaryotic proteins have an average size of 472 amino acid residues (aa), average protein sizes in plant genomes are smaller than those of animals and fungi. Proteins unique to plants are ∼81aa shorter than plant proteins conserved among other eukaryotic lineages. The smaller average size of plant proteins could neither be explained by endosymbiosis nor subcellular compartmentation nor exon size, but rather due to exon number. Metazoan proteins are encoded on average by ∼10 exons of small size [∼176 nucleotides (nt)]. Streptophyta have on average only ∼5.7 exons of medium size (∼230nt). Multicellular species code for large proteins by increasing the exon number, while most unicellular organisms employ rather larger exons (>400nt). Among subcellular compartments, membrane proteins are the largest (∼520aa), whereas the smallest proteins correspond to the gene ontology group of ribosome (∼240aa). Plant genes are encoded by half the number of exons and also contain fewer domains than animal proteins on average. Interestingly, endosymbiotic proteins that migrated to the plant nucleus became larger than their cyanobacterial orthologs. We thus conclude that plants have proteins larger than bacteria but smaller than animals or fungi. Compared to the average of eukaryotic species, plants have ∼34% more but ∼20% smaller proteins. This suggests that photosynthetic organisms are unique and deserve therefore special attention with regard to the evolutionary forces acting on their genomes and proteomes.

Variation in Key Flavonoid Biosynthetic Enzymes and Phytochemicals in 'Rio Red' Grapefruit (Citrus paradisi Macf.) during Fruit Development

In the current study, the phytochemical contents and expression of genes involved in flavonoid biosynthesis in Rio Red grapefruit were studied at different developmental and maturity stages for the first time. Grapefruit were harvested in June, August, November, January, and April and analyzed for the levels of carotenoids, vitamin C, limonoids, flavonoids, and furocoumarins by HPLC. In addition, genes encoding for phenylalanine ammonia-lyase (PAL), chalcone synthase (CHS), chalcone isomerase (CHI), and 1,2-rhamnosyltransferase (2RT) were isolated, and their expression in grapefruit juice vesicles was studied. Fruit maturity had significant influence on the expression of the genes, with PAL, CHS, and CHI having higher expression in immature fruits (June), whereas 2RT expression was higher in mature fruits (November and January). The levels of flavonoids (except naringin and poncirin), vitamin C, and furocoumarins gradually decreased from June to April. Furthermore, limonin levels sharply decreased in January. Lycopene decreased whereas β-carotene gradually increased with fruit maturity. Naringin did not exactly follow the pattern of 2RT or of PAL, CHS, and CHI expression, indicating that the four genes may have complementary effects on the level of naringin. Nevertheless, of the marketable fruit stages, early-season grapefruits harvested in November contained more beneficial phytochemicals as compared to mid- and late-season fruits harvested in January and April, respectively.

Factors Affecting the Radiosensitivity of Hexaploid Wheat to γ-Irradiation: Radiosensitivity of Hexaploid Wheat (Triticum aestivum L.) [Corrected]

Understanding the radiosensitivity of plants, an important factor in crop mutation breeding programs, requires a thorough investigation of the factors that contribute to this trait. In this study, we used the highly radiosensitive wheat (Triticum aestivum L.) variety HY1 and J411, a γ-irradiation-insensitive control, which were screened from a natural population, to examine the factors affecting radiosensitivity, including free radical content and total antioxidant capacity, as well as the expression of TaKu70 and TaKu80 (DNA repair-related genes) as measured by real-time PCR. We also investigated the alternative splicing of this gene in the wild-type wheat ecotype by sequence analysis. Free radical contents and total antioxidant capacity significantly increased upon exposure of HY1 wheat to γ-irradiation in a dose-dependent manner. By contrast, in J411, the free radical contents exhibited a similar trend, but the total antioxidant capacity exhibited a downward trend upon increasing γ-irradiation. Additionally, we detected dose-dependent increases in TaKu70 and TaKu80 expression levels in γ-irradiated HY1, while in J411, TaKu70 expression levels increased, followed by a decline. We also detected alternative splicing of TaKu70 mRNA, namely, intron retention, in HY1 but not in J411. Our findings indicate that γ-irradiation induces oxidative stress and DNA damage in hexaploid wheat, resulting in growth retardation of seedlings, and they suggest that TaKu70 may play a causal role in radiosensitivity in HY1. Further studies are required to exploit these factors to improve radiosensitivity in other wheat varieties.