Genome-wide survey of alternative splicing in the grass Brachypodium distachyon: a emerging model biosystem for plant functional genomics

Investigation for a multi-silique trait in Brassica napus by alternative splicing analysis

Background

Flower and fruit development are vital stages of the angiosperm lifecycle. We previously investigated the multi-silique trait in the rapeseed (Brassica napus) line zws-ms on a genomic and transcriptomic level, leading to the identification of two genomic regions and several candidate genes associated with this trait. However, some events on the transcriptome level, like alternative splicing, were poorly understood.

Methods

Plants from zws-ms and its near-isogenic line (NIL) zws-217 were both grown in Xindu with normal conditions and a colder area Ma'erkang. Buds from the two lines were sampled and RNA was isolated to perform the transcriptomic sequencing. The numbers and types of alternative splicing (AS) events from the two lines were counted and classified. Genes with AS events and expressed differentially between the two lines, as well as genes with AS events which occurred in only one line were emphasized. Their annotations were further studied.

Results

From the plants in Xindu District, an average of 205,496 AS events, which could be sorted into five AS types, were identified. zws-ms and zws-217 shared highly similar ratios of each AS type: The alternative 5' and 3' splice site types were the most common, while the exon skipping type was observed least often. Eleven differentially expressed AS genes were identified, of which four were upregulated and seven were downregulated in zws-ms. Their annotations implied that five of these genes were directly associated with the multi-silique trait. While samples from colder area Ma'erkang generated generally reduced number of each type of AS events except for Intron Retention; but the number of differentially expressed AS genes increased significantly. Further analysis found that among the 11 differentially expressed AS genes from Xindu, three of them maintained the same expression models, while the other eight genes did not show significant difference between the two lines in expression level. Additionally, the 205 line-specific expressed AS genes were analyzed, of which 187 could be annotated, and two were considered to be important.

Discussion

This study provides new insights into the molecular mechanism of the agronomically important multi-silique trait in rapeseed on the transcriptome level and screens out some environment-responding candidate genes.

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

Transcriptome analysis of two contrasting rice cultivars during alkaline stress

Soil alkalinity greatly affects plant growth and crop productivity. Although RNA-Seq analyses have been conducted to investigate genome-wide gene expression in response to alkaline stress in many plants, the expressions of alkali-responsive genes in rice have not been previously investigated. In this study, the transcriptomic data between an alkaline-tolerant (WD20342) and an alkaline-sensitive (Caidao) rice cultivar were compared under alkaline stress conditions. A total of 962 important alkali-responsive (IAR) genes from highly differentially expressed genes (DEGs) were identified, including 28 alkaline-resistant cultivar-related genes, 771 alkaline-sensitive cultivar-related genes and 163 cultivar-non-specific genes. Gene ontology (GO) analysis indicated the enrichment of IAR genes involved in various stimulus or stress responses. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, the IAR genes were related primarily to plant hormone signal transduction and biosynthesis of secondary metabolites. Additionally, among these 962 IAR genes, 74 were transcription factors and 15 occurred with differential alternative splicing between the different samples after alkaline treatment. Our results provide a valuable resource on alkali-responsive genes and should benefit the improvement of alkaline stress tolerance in rice.

Advances in Transcriptomics of Plants

The current global population of 7.3 billion is estimated to reach 9.7 billion in the year 2050. Rapid population growth is driving up global food demand. Additionally, global climate change, environmental degradation, drought, emerging diseases, and salty soils are the current threats to global food security. In order to mitigate the adverse effects of these diverse agricultural productivity constraints and enhance crop yield and stress-tolerance in plants, we need to go beyond traditional and molecular plant breeding. The powerful new tools for genome editing, Transcription Activator-Like Effector Nucleases (TALENs) and Clustered Regulatory Interspaced Short Palindromic Repeats (CRISPR)/Cas systems (CRISPR-Cas9), have been hailed as a quantum leap forward in the development of stress-resistant plants. Plant breeding techniques, however, have several drawbacks. Hence, identification of transcriptional regulatory elements and deciphering mechanisms underlying transcriptional regulation are crucial to avoiding unintended consequences in modified crop plants, which could ultimately have negative impacts on human health. RNA splicing as an essential regulated post-transcriptional process, alternative polyadenylation as an RNA-processing mechanism, along with non-coding RNAs (microRNAs, small interfering RNAs and long non-coding RNAs) have been identified as major players in gene regulation. In this chapter, we highlight new findings on the essential roles of alternative splicing and alternative polyadenylation in plant development and response to biotic and abiotic stresses. We also discuss biogenesis and the functions of microRNAs (miRNAs) and small interfering RNAs (siRNAs) in plants and recent advances in our knowledge of the roles of miRNAs and siRNAs in plant stress response. Graphical Abstract.

Transcriptome analysis of Brachypodium during fungal pathogen infection reveals both shared and distinct defense responses with wheat

Fusarium crown rot (FCR) of wheat and barley, predominantly caused by the fungal pathogen Fusarium pseudograminearum, is a disease of economic significance. The quantitative nature of FCR resistance within cultivated wheat germplasm has significantly limited breeding efforts to enhanced FCR resistance in wheat. In this study, we characterized the molecular responses of Brachypodium distachyon (Brachypodium hereafter) to F. pseudograminearum infection using RNA-seq to determine whether Brachypodium can be exploited as a model system towards better understanding of F. pseudograminearum-wheat interaction. The transcriptional response to infection in Brachypodium was strikingly similar to that previously reported in wheat, both in shared expression patterns of wheat homologs of Brachypodium genes and functional overlap revealed through comparative gene ontology analysis in both species. Metabolites produced by various biosynthetic pathways induced in both wheat and Brachypodium were quantified, revealing a high degree of overlap between these two species in metabolic response to infection but also showed Brachypodium does not produce certain defence-related metabolites found in wheat. Functional analyses of candidate genes identified in this study will improve our understanding of resistance mechanisms and may lead to the development of new strategies to protect cereal crops from pathogen infection.

An event of alternative splicing affects the expression of the NTRC gene, encoding NADPH-thioredoxin reductase C, in seed plants

The NTRC gene encodes a NADPH-dependent thioredoxin reductase with a joint thioredoxin domain, exclusive of photosynthetic organisms. An updated search shows that although most species harbor a single copy of the NTRC gene, two copies were identified in different species of the genus Solanum, Glycine max and the moss Physcomitrella patens. The phylogenetic analysis of NTRCs from different sources produced a tree with the major groups of photosynthetic organisms: cyanobacteria, algae and land plants, indicating the evolutionary success of the NTRC gene among photosynthetic eukaryotes. An event of alternative splicing affecting the expression of the NTRC gene was identified, which is conserved in seed plants but not in algae, bryophytes and lycophytes. The alternative splicing event results in a transcript with premature stop codon, which would produce a truncated form of the enzyme. The standard splicing/alternative splicing (SS/AS) transcripts ratio was higher in photosynthetic tissues from Arabidopsis, Brachypodium and tomato, in line with the higher content of the NTRC polypeptide in these tissues. Moreover, environmental stresses such as cold or high salt affected the SS/AS ratio of the NTRC gene transcripts in Brachypodium seedlings. These results suggest that the alternative splicing of the NTRC gene might be an additional mechanism for modulating the content of NTRC in photosynthetic and non-photosynthetic tissues of seed plants.

Heat Stress Regulates the Expression of Genes at Transcriptional and Post-Transcriptional Levels, Revealed by RNA-seq in Brachypodium distachyon

Heat stress greatly affects plant growth/development and influences the output of crops. With the increased occurrence of extreme high temperature, the negative influence on cereal products from heat stress becomes severer and severer. It is urgent to reveal the molecular mechanism in response to heat stress in plants. In this research, we used RNA-seq technology to identify differentially expressed genes (DEGs) in leaves of seedlings, leaves and inflorescences at heading stage of Brachypodium distachyon, one model plant of grasses. Results showed many genes in responding to heat stress. Of them, the expression level of 656 DEGs were altered in three groups of samples treated with high temperature. Gene ontology (GO) analysis showed that the highly enriched DEGs were responsible for heat stress and protein folding. According to KEGG pathway analysis, the DEGs were related mainly to photosynthesis-antenna proteins, the endoplasmic reticulum, and the spliceosome. Additionally, the expression level of 454 transcription factors belonging to 49 gene families was altered, as well as 1,973 splicing events occurred after treatment with high temperature. This research lays a foundation for characterizing the molecular mechanism of heat stress response and identifying key genes for those responses in plants. These findings also clearly show that heat stress regulates the expression of genes not only at transcriptional level, but also at post-transcriptional level.

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.

nagnag: Identification and quantification of NAGNAG alternative splicing using RNA-Seq data

Regulation of proteome diversity by alternative splicing has been widely demonstrated in plants and animals. NAGNAG splicing, which was recently defined as a tissue specific event, results in the production of two distinct isoforms that are distinguished by three nucleotides (NAG) as a consequence of the intron proximal or distal to the splice site. Since the NAGNAG mechanism is not well characterized, tools for the identification and quantification of NAGNAG splicing events remain under-developed. Here we report nagnag, an R-based NAGNAG splicing detection tool, which accurately identifies and quantifies NAGNAG splicing events using RNA-Seq. Overall, nagnag produces user-friendly visualization reports and highlights differences between the DNA/RNA/protein across the identified isoforms of the reported gene. The package is available on https://sourceforge.net/projects/nagnag/files/; or http://genome.sdau.edu.cn/research/software/nagnag.html.

SplicingTypesAnno: annotating and quantifying alternative splicing events for RNA-Seq data

Alternative splicing plays a key role in the regulation of the central dogma. Four major types of alternative splicing have been classified as intron retention, exon skipping, alternative 5 splice sites or alternative donor sites, and alternative 3 splice sites or alternative acceptor sites. A few algorithms have been developed to detect splice junctions from RNA-Seq reads. However, there are few tools targeting at the major alternative splicing types at the exon/intron level. This type of analysis may reveal subtle, yet important events of alternative splicing, and thus help gain deeper understanding of the mechanism of alternative splicing. This paper describes a user-friendly R package, extracting, annotating and analyzing alternative splicing types for sequence alignment files from RNA-Seq. SplicingTypesAnno can: (1) provide annotation for major alternative splicing at exon/intron level. By comparing the annotation from GTF/GFF file, it identifies the novel alternative splicing sites; (2) offer a convenient two-level analysis: genome-scale annotation for users with high performance computing environment, and gene-scale annotation for users with personal computers; (3) generate a user-friendly web report and additional BED files for IGV visualization. SplicingTypesAnno is a user-friendly R package for extracting, annotating and analyzing alternative splicing types at exon/intron level for sequence alignment files from RNA-Seq. It is publically available at https://sourceforge.net/projects/splicingtypes/files/ or http://genome.sdau.edu.cn/research/software/SplicingTypesAnno.html.

BARE retrotransposons are translated and replicated via distinct RNA pools

The replication of Long Terminal Repeat (LTR) retrotransposons, which can constitute over 80% of higher plant genomes, resembles that of retroviruses. A major question for retrotransposons and retroviruses is how the two conflicting roles of their transcripts, in translation and reverse transcription, are balanced. Here, we show that the BARE retrotransposon, despite its organization into just one open reading frame, produces three distinct classes of transcripts. One is capped, polyadenylated, and translated, but cannot be copied into cDNA. The second is not capped or polyadenylated, but is destined for packaging and ultimate reverse transcription. The third class is capped, polyadenylated, and spliced to favor production of a subgenomic RNA encoding only Gag, the protein forming virus-like particles. Moreover, the BARE2 subfamily, which cannot synthesize Gag and is parasitic on BARE1, does not produce the spliced sub-genomic RNA for translation but does make the replication competent transcripts, which are packaged into BARE1 particles. To our knowledge, this is first demonstration of distinct RNA pools for translation and transcription for any retrotransposon.

Genome-wide landscape of alternative splicing events in Brachypodium distachyon

Recently, Brachypodium distachyon has emerged as a model plant for studying monocot grasses and cereal crops. Using assembled expressed transcript sequences and subsequent mapping to the corresponding genome, we identified 1219 alternative splicing (AS) events spanning across 2021 putatively assembled transcripts generated from 941 genes. Approximately, 6.3% of expressed genes are alternatively spliced in B. distachyon. We observed that a majority of the identified AS events were related to retained introns (55.5%), followed by alternative acceptor sites (16.7%). We also observed a low percentage of exon skipping (5.0%) and alternative donor site events (8.8%). The 'complex event' that consists of a combination of two or more basic splicing events accounted for ∼14.0%. Comparative AS transcript analysis revealed 163 and 39 homologous pairs between B. distachyon and Oryza sativa and between B. distachyon and Arabidopsis thaliana, respectively. In all, we found 16 AS transcripts to be conserved in all 3 species. AS events and related putative assembled transcripts annotation can be systematically browsed at Plant Alternative Splicing Database (http://proteomics.ysu.edu/altsplice/plant/).

Alternative splicing: enhancing ability to cope with stress via transcriptome plasticity

Alternative splicing is a mechanism for the regulation of gene expression that is widespread in higher eukaryotes. Genome-wide approaches, based on comparison of expressed and genomic sequences, on tiling arrays, and on next-generation sequencing, have provided growing knowledge of the extent, distribution and association of alternative splicing with stress-related genes in plants. The functional meaning of alternative splicing in response to stress has been defined for many genes, and in particular for those involved in the regulation of the stress responses, such as protein kinases, transcription factors, splicing regulators and pathogen-resistance genes. The production of proteins with diverse domain rearrangements from the same gene is the main alternative splicing mechanism for pathogen-resistance genes. The plant response to abiotic stress is also characterized by a second mechanism, which consists of the expression of alternative transcripts that are targeted to nonsense-mediated decay. These quantitatively regulate stress-related gene expression. Many alternative splicing events are well conserved among plant species, and also across kingdoms, especially those observed in response to stress, for genes encoding splicing regulators, and other classes of RNA-binding proteins. Nevertheless, non-conserved events indicate that alternative splicing represents an evolutionary strategy that rapidly increases genome plasticity and develops new gene functions, along with other mechanisms such as gene duplication. Finally, the study of the naturally occurring variability of alternative splicing and the identification of genomic regions involved in the regulation of alternative splicing in crops are proposed as strategies for selecting genotypes with superior performance under adverse environmental conditions.

Deciphering the plant splicing code: experimental and computational approaches for predicting alternative splicing and splicing regulatory elements

Extensive alternative splicing (AS) of precursor mRNAs (pre-mRNAs) in multicellular eukaryotes increases the protein-coding capacity of a genome and allows novel ways to regulate gene expression. In flowering plants, up to 48% of intron-containing genes exhibit AS. However, the full extent of AS in plants is not yet known, as only a few high-throughput RNA-Seq studies have been performed. As the cost of obtaining RNA-Seq reads continues to fall, it is anticipated that huge amounts of plant sequence data will accumulate and help in obtaining a more complete picture of AS in plants. Although it is not an onerous task to obtain hundreds of millions of reads using high-throughput sequencing technologies, computational tools to accurately predict and visualize AS are still being developed and refined. This review will discuss the tools to predict and visualize transcriptome-wide AS in plants using short-reads and highlight their limitations. Comparative studies of AS events between plants and animals have revealed that there are major differences in the most prevalent types of AS events, suggesting that plants and animals differ in the way they recognize exons and introns. Extensive studies have been performed in animals to identify cis-elements involved in regulating AS, especially in exon skipping. However, few such studies have been carried out in plants. Here, we review the current state of research on splicing regulatory elements (SREs) and briefly discuss emerging experimental and computational tools to identify cis-elements involved in regulation of AS in plants. The availability of curated alternative splice forms in plants makes it possible to use computational tools to predict SREs involved in AS regulation, which can then be verified experimentally. Such studies will permit identification of plant-specific features involved in AS regulation and contribute to deciphering the splicing code in plants.