Genome-wide survey of Alternative Splicing in Sorghum Bicolor

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology, bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization, and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside a gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

Analysis of Whole-Transcriptome RNA-Seq Data Reveals the Involvement of Alternative Splicing in the Drought Response of Glycyrrhiza uralensis

Alternative splicing (AS) is a post-transcriptional regulatory mechanism that increases protein diversity. There is growing evidence that AS plays an important role in regulating plant stress responses. However, the mechanism by which AS coordinates with transcriptional regulation to regulate the drought response in Glycyrrhiza uralensis remains unclear. In this study, we performed a genome-wide analysis of AS events in G. uralensis at different time points under drought stress using a high-throughput RNA sequencing approach. We detected 2,479 and 2,764 AS events in the aerial parts (AP) and underground parts (UP), respectively, of drought-stressed G. uralensis. Of these, last exon AS and exon skipping were the main types of AS. Overall, 2,653 genes undergoing significant AS regulation were identified from the AP and UP of G. uralensis exposed to drought for 2, 6, 12, and 24 h. Gene Ontology analyses indicated that AS plays an important role in the regulation of nitrogen and protein metabolism in the drought response of G. uralensis. Notably, the spliceosomal pathway and basal transcription factor pathway were significantly enriched with differentially spliced genes under drought stress. Genes related to splicing regulators in the AP and UP of G. uralensis responded to drought stress and underwent AS under drought conditions. In summary, our data suggest that drought-responsive AS directly and indirectly regulates the drought response of G. uralensis. Further in-depth studies on the functions and mechanisms of AS during abiotic stresses will provide new strategies for improving plant stress resistance.

Breeding for Climate Change Resilience: A Case Study of Loblolly Pine (Pinus taeda L.) in North America

Earth's atmosphere is warming and the effects of climate change are becoming evident. A key observation is that both the average levels and the variability of temperature and precipitation are changing. Information and data from new technologies are developing in parallel to provide multidisciplinary opportunities to address and overcome the consequences of these changes in forest ecosystems. Changes in temperature and water availability impose multidimensional environmental constraints that trigger changes from the molecular to the forest stand level. These can represent a threat for the normal development of the tree from early seedling recruitment to adulthood both through direct mortality, and by increasing susceptibility to pathogens, insect attack, and fire damage. This review summarizes the strengths and shortcomings of previous work in the areas of genetic variation related to cold and drought stress in forest species with particular emphasis on loblolly pine (Pinus taeda L.), the most-planted tree species in North America. We describe and discuss the implementation of management and breeding strategies to increase resilience and adaptation, and discuss how new technologies in the areas of engineering and genomics are shaping the future of phenotype-genotype studies. Lessons learned from the study of species important in intensively-managed forest ecosystems may also prove to be of value in helping less-intensively managed forest ecosystems adapt to climate change, thereby increasing the sustainability and resilience of forestlands for the future.

Drought-Induced Root Pressure in Sorghum bicolor

Root pressure, also manifested as profusive sap flowing from cut stems, is a phenomenon in some species that has perplexed biologists for much of the last century. It is associated with increased crop production under drought, but its function and regulation remain largely unknown. In this study, we investigated the initiation, mechanisms, and possible adaptive function of root pressure in six genotypes of Sorghum bicolor during a drought experiment in the greenhouse. We observed that root pressure was induced in plants exposed to drought followed by re-watering but possibly inhibited by 100% re-watering in some genotypes. We found that root pressure in drought stressed and re-watered plants was associated with greater ratio of fine: coarse root length and shoot biomass production, indicating a possible role of root allocation in creating root pressure and adaptive benefit of root pressure for shoot biomass production. Using RNA-Seq, we identified gene transcripts that were up- and down-regulated in plants with root pressure expression, focusing on genes for aquaporins, membrane transporters, and ATPases that could regulate inter- and intra-cellular transport of water and ions to generate positive xylem pressure in root tissue.

Weighted gene co-expression network analysis identifies modules and functionally enriched pathways in the lactation process

The exponential growth in knowledge has resulted in a better understanding of the lactation process in a wide variety of animals. However, the underlying genetic mechanisms are not yet clearly known. In order to identify the mechanisms involved in the lactation process, various mehods, including meta-analysis, weighted gene co-express network analysis (WGCNA), hub genes identification, gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment at before peak (BP), peak (P), and after peak (AP) stages of the lactation processes have been employed. A total of 104, 85, and 26 differentially expressed genes were identified based on PB vs. P, BP vs. AP, and P vs. AP comparisons, respectively. GO and KEGG pathway enrichment analysis revealed that DEGs were significantly enriched in the "ubiquitin-dependent ERAD" and the "chaperone cofactor-dependent protein refolding" in BP vs. P and P vs. P, respectively. WGCNA identified five significant functional modules related to the lactation process. Moreover, GJA1, AP2A2, and NPAS3 were defined as hub genes in the identified modules, highlighting the importance of their regulatory impacts on the lactation process. The findings of this study provide new insights into the complex regulatory networks of the lactation process at three distinct stages, while suggesting several candidate genes that may be useful for future animal breeding programs. Furthermore, this study supports the notion that in combination with a meta-analysis, the WGCNA represents an opportunity to achieve a higher resolution analysis that can better predict the most important functional genes that might provide a more robust bio-signature for phenotypic traits, thus providing more suitable biomarker candidates for future studies.

Weighted gene co-expression network analysis of the salt-responsive transcriptomes reveals novel hub genes in green halophytic microalgae Dunaliella salina

Despite responses to salinity stress in Dunaliella salina, a unicellular halotolerant green alga, being subject to extensive study, but the underlying molecular mechanism remains unknown. Here, Empirical Bayes method was applied to identify the common differentially expressed genes (DEGs) between hypersaline and normal conditions. Then, using weighted gene co-expression network analysis (WGCNA), which takes advantage of a graph theoretical approach, highly correlated genes were clustered as a module. Subsequently, connectivity patterns of the identified modules in two conditions were surveyed to define preserved and non-preserved modules by combining the Zsummary and medianRank measures. Finally, common and specific hub genes in non-preserved modules were determined using Eigengene-based module connectivity or module membership (k_ME) measures and validation was performed by using leave-one-out cross-validation (LOOCV). In this study, the power of beta = 12 (scale-free R2 = 0.8) was selected as the soft-thresholding to ensure a scale-free network, which led to the identification of 15 co-expression modules. Results also indicate that green, blue, brown, and yellow modules are non-preserved in salinity stress conditions. Examples of enriched Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in non-preserved modules are Sulfur metabolism, Oxidative phosphorylation, Porphyrin and chlorophyll metabolism, Vitamin B6 metabolism. Moreover, the systems biology approach was applied here, proposed some salinity specific hub genes, such as radical-induced cell death1 protein (RCD1), mitogen-activated protein kinase kinase kinase 13 (MAP3K13), long-chain acyl-CoA synthetase (ACSL), acetyl-CoA carboxylase, biotin carboxylase subunit (AccC), and fructose-bisphosphate aldolase (ALDO), for the development of metabolites accumulating strains in D. salina.

Integrative analysis of gene expression and alternative splicing in microalgae grown under heterotrophic condition

Heterotrophic cultures are the most effective approach to overcome low growth rate challenge in the most commercial microalgae. However, the mechanism through which heterotrophic condition regulates algae metabolism are not completely clear. Alternative Splicing (AS) is a common posttranscriptional process by which transcriptome and proteome plasticity increases at different environmental conditions. To identify and characterize of AS events in Auxenochlorella protothecoides microalga grown in autotrophic and heterotrophic, RNA-Seq data were analysed. We found that AS increased with the transition from autotrophic to heterotrophic condition. 705 and 660 differentially expressed (DEG) and spliced (DAS) genes were identified for A.protothecoides was transferred from autotrophic to heterotrophic condition, respectively. Moreover, there was slight coverage between DEG and DAS genes. Furthermore, functional analysis showed that the DAS genes are most frequently related to ion binding and stimulus response. The results also indicated that prevalence of Intron retention is associated with down-regulation of the genes involved in carotenoid biosynthesis. This study provides valuable insights into transcriptional and posttranscriptional plasticity of microalgae during growth mode change.

Genomic-based root plasticity to enhance abiotic stress adaptation and edible yield in grain crops

Phenotypic plasticity refers to changes expressed by a genotype across different environments and is one of the major means by which plants cope with environmental variability. Multi-fold differences in phenotypic plasticity have been noted across crops, with wild ancestors and landraces being more plastic than crops when under stress. Plasticity in response to abiotic stress adaptation, plant architecture, physio-reproductive and quality traits are multi-genic (QTL). Plasticity QTL (pQTL) were either collocated with main effect QTL and QEI (QTL × environment interaction) or located independently from the main effect QTL. For example, variations in root plasticity have been successfully introgressed to enhance abiotic stress adaptation in rice. The independence of genetic control of a trait and of its plasticity suggests that breeders may select for high or low plasticity in combination with high or low performance of economically important traits. Trait plasticity in stressful environments may be harnessed through breeding stress-tolerant crops. There exists a genetic cost associated with plasticity, so a better understanding of the trade-offs between plasticity and productivity is warranted prior to undertaking breeding for plasticity traits together with productivity in stress environments.

Integration of Cross Species RNA-seq Meta-Analysis and Machine-Learning Models Identifies the Most Important Salt Stress-Responsive Pathways in Microalga Dunaliella

Photosynthetic microalgae are potentially yielding sources of different high-value secondary metabolites. Salinity is a complex stress that influences various metabolite-related pathways in microalgae. To obtain a clear view of the underlying metabolic pathways and resolve contradictory information concerning the transcriptional regulation of Dunaliella species in salt stress conditions, RNA-seq meta-analysis along with systems levels analysis was conducted. A p-value combination technique with Fisher method was used for cross species meta-analysis on the transcriptomes of two Dunaliella salina and Dunaliellatertiolecta species. The potential functional impacts of core meta-genes were surveyed based on gene ontology and network analysis. In the current study, the integration of supervised machine-learning algorithms with RNA-seq meta-analysis was performed. The analysis shows that the lipid and nitrogen metabolism, structural proteins of photosynthesis apparatus, chaperone-mediated autophagy, and ROS-related genes are the keys and core elements of the Dunaliella salt stress response system. Cross-talk between Ca²⁺ signal transduction, lipid accumulation, and ROS signaling network in salt stress conditions are also proposed. Our novel approach opens new avenues for better understanding of microalgae stress response mechanisms and for selection of candidate gene targets for metabolite production in microalgae.

Genome-wide identification and co-expression network analysis of nuclear factor-Y in barley revealed potential functions in salt stress

Nuclear factor-Ys (NF-Ys) were previously shown to have important regulatory impacts in different developmental and physiological process. However, in barley the function of the NF-Y genes at system levels is not well known. To identify barley NF-Ys, Arabidopsis and wheat NF-Y protein sequences were retrieved and the BLAST program along with the hidden Markov model were used. Multiple sequence alignments of identified NF-Ys were constructed using ClustalW. Expression patterns of the NF-Ys at different physiological and developmental conditions were also surveyed based on microarray datasets in public databases and subsequently co-expression network were constructed. Validation of in silico expression analysis was performed by real-time qPCR under salt stress condition. In total, 23 barley NF-Ys (8 NF-YA, 11 NF-YB and 4 NF-YC) were identified. Based on the sequence homology, the subunits of the NF-Y complex were divided into three to five groups. Structural analysis highlighted the conserved domains of HvNF-YA, HvNF-YB and HvNF-YC. Co-expression network analysis indicated the potential functions of HvNF-Ys in photosynthesis, starch biosynthesis and osmotic stress tolerance. The results of qRT-PCR also confirmed the HvNF-Ys roles in adaptation responses of barley to salt stress. We identified some potential candidate genes which could be used for improvements of cereals tolerance to salinity stress.

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).

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).

Alternative splicing complexity contributes to genetic improvement of drought resistance in the rice maintainer HuHan2B

Water-saving and drought-resistantce rice (WDR) breeding practices have greatly increased grain yield and drought resistance. To study the genetic basis of adaptation to drought, transcriptome sequences from the WDR maintainer line HuHan2B and the recurrent parent HanFengB were analyzed for alternative splicing (AS) complexity. Intron retention, the dominant AS type, accounted for 42% of the observed AS events. Differential expression analysis revealed transcripts were preferentially expressed in different varieties and conditions. Based on gene ontology predictions, the biological functions of drought-induced transcripts were significantly enriched in genes involved in transcription regulation, chloroplast components and response to abiotic stimulus in HuHan2B, whereas developmental processes for reproduction were primarily enriched in HanFengB. The regulatory network of transcription factors was driven by cohorts of transcript splicing targets, resulting in more diversified regulatory relationships due to AS complexity than in our previous findings. Moreover, several genes were validated to accumulate novel splicing transcripts in a drought-induced manner. Together, these results suggest that HuHan2B and HanFengB share similar AS features but that a subset of genes with increased levels of AS involved in transcription regulatory networks may contribute an additional level of control for genetic improvement of drought resistance in rice maintainer HuHan2B through breeding.

Domestication reduces alternative splicing expression variations in sorghum

Domestication is known to strongly reduce genomic diversity through population bottlenecks. The resulting loss of polymorphism has been thoroughly documented in numerous cultivated species. Here we investigate the impact of domestication on the diversity of alternative transcript expressions using RNAseq data obtained on cultivated and wild sorghum accessions (ten accessions for each pool). In that aim, we focus on genes expressing two isoforms in sorghum and estimate the ratio between expression levels of those isoforms in each accession. Noticeably, for a given gene, one isoform can either be overexpressed or underexpressed in some wild accessions, whereas in the cultivated accessions, the balance between the two isoforms of the same gene appears to be much more homogenous. Indeed, we observe in sorghum significantly more variation in isoform expression balance among wild accessions than among domesticated accessions. The possibility exists that the loss of nucleotide diversity due to domestication could affect regulatory elements, controlling transcription or degradation of these isoforms. Impact on the isoform expression balance is discussed. As far as we know, this is the first time that the impact of domestication on transcript isoform balance has been studied at the genomic scale. This could pave the way towards the identification of key domestication genes with finely tuned isoform expressions in domesticated accessions while being highly variable in their wild relatives.

Evolutionarily Conserved Alternative Splicing Across Monocots

One difficulty when identifying alternative splicing (AS) events in plants is distinguishing functional AS from splicing noise. One way to add confidence to the validity of a splice isoform is to observe that it is conserved across evolutionarily related species. We use a high throughput method to identify junction-based conserved AS events from RNA-Seq data across nine plant species, including five grass monocots (maize, sorghum, rice, Brachpodium, and foxtail millet), plus two nongrass monocots (banana and African oil palm), the eudicot Arabidopsis, and the basal angiosperm Amborella In total, 9804 AS events were found to be conserved between two or more species studied. In grasses containing large regions of conserved synteny, the frequency of conserved AS events is twice that observed for genes outside of conserved synteny blocks. In plant-specific RS and RS2Z subfamilies of the serine/arginine (SR) splice-factor proteins, we observe both conservation and divergence of AS events after the whole genome duplication in maize. In addition, plant-specific RS and RS2Z splice-factor subfamilies are highly connected with R2R3-MYB in STRING functional protein association networks built using genes exhibiting conserved AS. Furthermore, we discovered that functional protein association networks constructed around genes harboring conserved AS events are enriched for phosphatases, kinases, and ubiquitylation genes, which suggests that AS may participate in regulating signaling pathways. These data lay the foundation for identifying and studying conserved AS events in the monocots, particularly across grass species, and this conserved AS resource identifies an additional layer between genotype to phenotype that may impact future crop improvement efforts.

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.

A survey of the sorghum transcriptome using single-molecule long reads

Alternative splicing and alternative polyadenylation (APA) of pre-mRNAs greatly contribute to transcriptome diversity, coding capacity of a genome and gene regulatory mechanisms in eukaryotes. Second-generation sequencing technologies have been extensively used to analyse transcriptomes. However, a major limitation of short-read data is that it is difficult to accurately predict full-length splice isoforms. Here we sequenced the sorghum transcriptome using Pacific Biosciences single-molecule real-time long-read isoform sequencing and developed a pipeline called TAPIS (Transcriptome Analysis Pipeline for Isoform Sequencing) to identify full-length splice isoforms and APA sites. Our analysis reveals transcriptome-wide full-length isoforms at an unprecedented scale with over 11,000 novel splice isoforms. Additionally, we uncover APA of ∼11,000 expressed genes and more than 2,100 novel genes. These results greatly enhance sorghum gene annotations and aid in studying gene regulation in this important bioenergy crop. The TAPIS pipeline will serve as a useful tool to analyse Iso-Seq data from any organism.

Alternative splicing and alternative polyadenylation (APA) contribute to mRNA diversity but are difficult to assess using short read RNA-seq data. Here, the authors use single molecule long-read isoform sequencing and develop a computational pipeline to identify full-length splice isoforms and APA sites in sorghum.

Evolutionary Character of Alternative Splicing in Plants

Alternative splicing (AS) is one of the most important ways to enhance the functional diversity of genes. Huge amounts of data have been produced by microarray, expressed sequence tag, and RNA-seq, and plenty of methods have been developed specifically for this task. The most frequently asked questions in previous research were as follows. What is the content rate of AS genes among the whole gene set? How many AS types are presented in the genome, and which type is dominant? How about the conservation ability of AS among different species? Which kinds of isoforms from some genes have the environmental response to help individual adaptation? Based on this background, we collected analysis results from 17 species to try to map out the landscape of AS studies in plants. We have noted the shortages of previous results, and we appeal to all scientists working in the AS field to make a standard protocol so that analyses between different projects are comparable.

Genome-wide cataloging and analysis of alternatively spliced genes in cereal crops

Background

Protein functional diversity at the post-transcriptional level is regulated through spliceosome mediated pre-mRNA alternative splicing (AS) events and that has been widely demonstrated to be a key player in regulating the functional diversity in plants. Identification and analysis of AS genes in cereal crop plants are critical for crop improvement and understanding regulatory mechanisms.

Results

We carried out the comparative analyses of the functional landscapes of the AS using the consensus assembly of expressed sequence tags and available mRNA sequences in four cereal plants. We identified a total of 8,734 in Oryza sativa subspecies (ssp) japonica, 2,657 in O. sativa ssp indica, 3,971 in Sorghum bicolor, and 10,687 in Zea mays AS genes. Among the identified AS events, intron retention remains to be the dominant type accounting for 23.5 % in S. bicolor, and up to 55.8 % in O. sativa ssp indica. We identified a total of 887 AS genes that were conserved among Z. mays, S. bicolor, and O. sativa ssp japonica; and 248 AS genes were found to be conserved among all four studied species or ssp. Furthermore, we identified 53 AS genes conserved with Brachypodium distachyon. Gene Ontology classification of AS genes revealed functional assignment of these genes in many biological processes with diverse molecular functions.

Conclusions

AS is common in cereal plants. The AS genes identified in four cereal crops in this work provide the foundation for further studying the roles of AS in regulation of cereal plant growth and development. The data can be accessed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/).

Electronic supplementary material

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