Regulation of alternative splicing of pre-mRNAs by stresses

Splice Variants of Superoxide Dismutases in Rice and Their Expression Profiles under Abiotic Stresses

The superoxide dismutases (SODs) play vital roles in controlling cellular reactive oxygen species (ROS) that are generated both under optimal as well as stress conditions in plants. The rice genome harbors seven SOD genes (CSD1, CSD2, CSD3, CSD4, FSD1, FSD2, and MSD) that encode seven constitutive transcripts. Of these, five (CSD2, CSD3, CSD4, FSD1, and MSD) utilizes an alternative splicing (AS) strategy and generate seven additional splice variants (SVs) or mRNA variants, i.e., three for CSD3, and one each for CSD2, CSD4, FSD1, and MSD. The exon-intron organization of these SVs revealed variations in the number and length of exons and/or untranslated regions (UTRs). We determined the expression patterns of SVs along with their constitutive forms of SODs in rice seedlings exposed to salt, osmotic, cold, heavy metal (Cu⁺²) stresses, as well as copper-deprivation. The results revealed that all seven SVs were transcriptionally active in both roots and shoots. When compared to their corresponding constitutive transcripts, the profiles of five SVs were almost similar, while two specific SVs (CSD3-SV4 and MSD-SV2) differed significantly, and the differences were also apparent between shoots and roots suggesting that the specific SVs are likely to play important roles in a tissue-specific and stress-specific manner. Overall, the present study has provided a comprehensive analysis of the SVs of SODs and their responses to stress conditions in shoots and roots of rice seedlings.

Global insights into duplicated gene expression and alternative splicing in polyploid Brassica napus under heat, cold, and drought stress

Polyploidy has been a prevalent process during plant evolution and it has made a major impact on the structure and evolution of plant genomes. Many important crop plants are polyploid. There is considerable interest in expression patterns of duplicated genes in polyploids. Alternative splicing (AS) is a fundamental aspect of gene expression that produces multiple final transcript types from a single type of mRNAs. The effects of abiotic stress conditions on AS in polyploids has received little attention. We conducted a global transcriptome analysis of Brassica napus, an allotetraploid derived from B. rapa (A_T ) and B. oleracea (C_T ), by RNA-Seq of plants subjected to cold, heat, and drought stress treatments. Analyses of 27,360 pairs of duplicated genes revealed overall A_T subgenome biases in gene expression and C_T subgenome biases in the extent of alternative splicing under all three stress treatments. More genes increased in expression than decreased in response to the stresses. Negative correlations were found between expression levels and AS frequency for each type of AS. Cold stress produced the greatest changes in gene expression and AS. Cold-induced AS changes were more likely to be shared with those generated by drought than by heat stress. We used homeologs of FLC and CCA1 as case studies to show the dynamics of how duplicates in a polyploid respond to cold stress. Our results suggest that divergence in gene expression and AS patterns between duplicated genes may increase the flexibility of polyploids when responding to abiotic stressors.

Genome-Wide Characterization of Alternative Splicing Events and Their Responses to Cold Stress in Tilapia

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism for cells to generate transcript variability and proteome diversity. No systematic investigation of AS events among different tissues in response to stressors is available for tilapia currently. In this study, AS among different tissues was identified and the cold stress-related AS events were explored in a Nile tilapia (Oreochromis niloticus) line based on 42 RNA-seq datasets using a bioinformatics pipeline. 14,796 (82.76%; SD = 2,840) of the expression genes showed AS events. The two most abundant AS types were alternative transcription start site (TSS) and terminal site (TTS) in tilapia. Testis, brain and kidney possess the most abundant AS gene number, while the blood, muscle and liver possess the least number in each tissue. Furthermore, 208 differentially alternative splicing (DAS) genes in heart and 483 DAS in brain in response to cold stress. The number of AS types for alternative exon end, exon skipping and retention of single intron increased significantly under cold stress. GO enrichment and pathway overrepresentation analysis indicated that many DAS genes, e.g., genes in circadian clock pathway, may influence expression of downstream genes under cold stress. Our study revealed that AS exists extensively in tilapia and plays an important role in cold adaption.

Comprehensive profiling of alternative splicing landscape during cold acclimation in tea plant

BACKGROUND

Alternative splicing (AS) may generate multiple mRNA splicing isoforms from a single mRNA precursor using different splicing sites, leading to enhanced diversity of transcripts and proteins. AS has been implicated in cold acclimation by affecting gene expression in various ways, yet little information is known about how AS influences cold responses in tea plant (Camellia sinensis).

RESULTS

In this study, the AS transcriptional landscape was characterized in the tea plant genome using high-throughput RNA-seq during cold acclimation. We found that more than 41% (14,103) of genes underwent AS events. We summarize the possible existence of 11 types of AS events, including the four common types of intron retention (IR), exon skipping (ES), alternative 5' splice site (A5SS), and alternative 3' splice site (A3SS); of these, IR was the major type in all samples. The number of AS events increased rapidly during cold treatment, but decreased significantly following de-acclimation (DA). It is notable that the number of differential AS genes gradually increased during cold acclimation, and these genes were enriched in pathways relating to oxidoreductase activity and sugar metabolism during acclimation and de-acclimation. Remarkably, the AS isoforms of bHLH transcription factors showed higher expression levels than their full-length ones during cold acclimation. Interestingly, the expression pattern of some AS transcripts of raffinose and sucrose synthase genes were significantly correlated with sugar contents.

CONCLUSION

Our findings demonstrated that changes in AS numbers and transcript expression may contribute to rapid changes in gene expression and metabolite profile during cold acclimation, suggesting that AS events play an important regulatory role in response to cold acclimation in tea plant.

The Effect of Drought on Transcriptome and Hormonal Profiles in Barley Genotypes With Contrasting Drought Tolerance

Like many cereal crops, barley is also negatively affected by drought stress. However, due to its simple genome as well as enhanced stress resilient nature compared to rice and wheat, barley has been considered as a model to decipher drought tolerance in cereals. In the present study, transcriptomic and hormonal profiles along with several biochemical features were compared between drought-tolerant (Otis) and drought-sensitive (Baronesse) barley genotypes subjected to drought to identify molecular and biochemical differences between the genotypes. The drought-induced decrease in the leaf relative water content, net photosynthesis, and biomass accumulation was relatively low in Otis compared to Baronesse. The hormonal profiles did not reveal significant differences for majority of the compounds other than the GA20 and the cis-zeatin-o-glucoside (c-ZOG), whose levels were greatly increased in Otis compared to Baronesse under drought. The major differences that emerged from the transcriptome analysis are; (1), the overall number of differentially expressed genes was relatively low in drought-tolerant Otis compared to drought-sensitive Baronesse; (2), a wax biosynthesis gene (CER1), and NAC transcription factors were specifically induced in Otis but not in Baronesse; (3), the degree of upregulation of betaine aldehyde dehydrogenase and a homeobox transcription factor (genes with proven roles in imparting drought tolerance), was greater in Otis compared to Baronesse; (4) the extent of downregulation of gene expression profiles for proteins of the reaction center photosystem II (PSII) (D1 and D2) was low in Otis compared to Baronesse; and, (5), alternative splicing (AS) was also found to differ between the genotypes under drought. Taken together, the overall transcriptional responses were low in drought-tolerant Otis but the genes that could confer drought tolerance were either specifically induced or greatly upregulated in the tolerant genotype and these differences could be important for drought tolerance in barley.

Identification and characterization of jasmonic acid- and linolenic acid-mediated transcriptional regulation of secondary laticifer differentiation in Hevea brasiliensis

Hevea brasiliensis remains the primary crop commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. Here, we described the transcriptional events related to jasmonic acid (JA)- and linolenic acid (LA)-induced secondary laticifer differentiation (SLD) in H. brasiliensis clone RRIM 600 based on RNA-seq approach. Histochemical approach proved that JA- and LA-treated samples resulted in SLD in H. brasiliensis when compared to ethephon and untreated control. RNA-seq data resulted in 86,614 unigenes, of which 2,664 genes were differentially expressed in JA and LA-induced secondary laticifer harvested from H. brasiliensis bark samples. Among these, 450 genes were unique to JA and LA as they were not differentially expressed in ethephon-treated samples compared with the untreated samples. Most transcription factors from the JA- and LA-specific dataset were classified under MYB, APETALA2/ethylene response factor (AP2/ERF), and basic-helix-loop-helix (bHLH) gene families that were involved in tissue developmental pathways, and we proposed that Bel5-GA2 oxidase 1-KNOTTED-like homeobox complex are likely involved in JA- and LA-induced SLD in H. brasiliensis. We also discovered alternative spliced transcripts, putative novel transcripts, and cis-natural antisense transcript pairs related to SLD event. This study has advanced understanding on the transcriptional regulatory network of SLD in H. brasiliensis.

Genome-wide characterization, molecular evolution and expression profiling of the metacaspases in potato (Solanum tuberosum L.)

Metacaspases are distant relatives of animal caspases found in plants, protozoa and fungi. Some recent studies have demonstrated that metacaspases are involved in regulating the developmental and environmentally induced programmed cell death in plants. In this study, we identified metacaspase gene family in potato (Solanum tuberosum L.) and analyzed their expression pattern in various developmental tissues and stress responses of plants. There were eight metacaspase genes identified in the Peptidase (Cysteine protease) C14 family and based upon sequence alignment and phylogenetic analysis, a systematic nomenclature of potato metacaspases (SotubMCs) has been proposed. Three of the eight candidate genes showing homology with Arabidopsis thaliana type I metacaspase, AtMC1 were given name SotubMC1, SotubMC2 and SotubMC3 as per the degree of relatedness. Similarly, the next three being homologous to A. thaliana type I metacaspase, AtMC3 were named SotubMC4, SotubMC5, and SotubMC6. The remaining two were named SotubMC7 and SotubMC8, showing significant similarity with type II metacaspases of A. thaliana, AtMC4 and AtMC9, respectively. Evolutionary divergence analysis of SotubMCs from its orthologs in seven other members of Solanaceae family as well as with A. thaliana, Vitis vinifera and Oryza sativa was also carried out. The d_N/d_S ratios of the orthologous pairs suggested the SotubMCs were under purifying (negative) selection in course of plant evolution. Splicing patterns of potato metacaspases were also analyzed. Amongst all SotubMCs, SotubMC2, SotubMC4, SotubMC6 and SotubMC7 genes appeared to produce multiple alternative spliced variants of different lengths. Furthermore using protein modeling tools, we have predicted the protein structure of identified metacaspases. The cis-regulatory elements analysis was also performed exhibiting the presence of development, stress and hormones related cis-elements in the promoter regions of the SotubMCs. This indicates that potato metacaspases might be playing important roles in the development, stress and hormone responsive pathways. Moreover, relative expression analysis of identified genes was carried out using qRT-PCR in various developmental tissues that also include stolons and tubers. The eight metacaspases showed differential expression in different tissues. Some of the tissues such as leaf undergoing senescence among different leaf developmental stages (immature, mature and senescent) displayed higher relative expression of some of the metacaspases, implying their involvement in leaf senescence. The expression pattern of SotubMCs under various abiotic, biotic and hormonal stresses was also analysed. The results showed that many members of the potato metacaspase gene family displayed differential expression patterns under various stress conditions. Taken together, the study could provide crucial resources for further investigations to understand the functional roles of the identified metacaspases in potato.

Transcript-level expression control of plant NLR genes

Plant NLR genes encode sensitive immune receptors that can mediate the specific recognition of pathogen avirulence effectors and activate a strong defence response, termed effector-triggered immunity. The expression of NLRs requires strict regulation, as their ability to trigger immunity is dependent on their dose, and overexpression of NLRs results in autoimmunity and massive fitness costs. An elaborate interplay of different mechanisms controlling NLR transcript levels allows plants to maximize their defence capacity, whilst limiting negative impact on their fitness. Global suppression of NLR transcripts may be a prerequisite for the fast evolution of new NLR variants and the expansion of this gene family. Here, we summarize recent progress made towards a comprehensive understanding of NLR transcript-level expression control. Multiple mechanistic steps, including transcription as well as co-/post-transcriptional processing and transcript turn-over, contribute to balanced base levels of NLR transcripts and allow for dynamic adjustments to defence situations.

A genetic screen implicates a CWC16/Yju2/CCDC130 protein and SMU1 in alternative splicing in Arabidopsis thaliana

To identify regulators of pre-mRNA splicing in plants, we developed a forward genetic screen based on an alternatively spliced GFP reporter gene in Arabidopsis thaliana In wild-type plants, three major splice variants issue from the GFP gene but only one represents a translatable GFP mRNA. Compared to wild-type seedlings, which exhibit an intermediate level of GFP expression, mutants identified in the screen feature either a "GFP-weak" or "Hyper-GFP" phenotype depending on the ratio of the three splice variants. GFP-weak mutants, including previously identified prp8 and rtf2, contain a higher proportion of unspliced transcript or canonically spliced transcript, neither of which is translatable into GFP protein. In contrast, the coilin-deficient hyper-gfp1 (hgf1) mutant displays a higher proportion of translatable GFP mRNA, which arises from enhanced splicing of a U2-type intron with noncanonical AT-AC splice sites. Here we report three new hgf mutants that are defective, respectively, in spliceosome-associated proteins SMU1, SmF, and CWC16, an Yju2/CCDC130-related protein that has not yet been described in plants. The smu1 and cwc16 mutants have substantially increased levels of translatable GFP transcript owing to preferential splicing of the U2-type AT-AC intron, suggesting that SMU1 and CWC16 influence splice site selection in GFP pre-mRNA. Genome-wide analyses of splicing in smu1 and cwc16 mutants revealed a number of introns that were variably spliced from endogenous pre-mRNAs. These results indicate that SMU1 and CWC16, which are predicted to act directly prior to and during the first catalytic step of splicing, respectively, function more generally to modulate splicing patterns in plants.

Understanding the mechanisms of dormancy in an invasive alien Sycamore lace bug, Corythucha ciliata through transcript and metabolite profiling

The sycamore lace bug, Corythucha ciliata, is a pest of sycamore trees. In China, it is found in the most northern border where it has been known to become dormant during harsh winters. But the molecular and metabolic basis for dormancy in this insect is still unknown. In this study, we analyzed the transcript and metabolite profiles of this bug to identify key genes and metabolites that are significantly regulated during dormancy in adult females and males. In total, 149 differentially expressed genes (DEGs) were significantly up-regulated and 337 DEGs were significantly down-regulated in dormant adults (both females and males). We found major differences in heat shock protein (HSPs), immunity-responsive genes, NAD-dependent deacetylase sirtuin-1 (SIRT1) and genes involved in the spliceosome pathway that is known to regulate stress. Among the 62 metabolites identified by GC-MS, 12 metabolites including glycerol, trehalose, and alanine were significantly increased during C. ciliata dormancy. By integrating the transcriptome and metabolite datasets, we found that the metabolites in glycolysis/gluconeogenesis and citrate cycle (TCA) were significantly reduced. This study is the first to report both transcript and metabolite profiles of the overwintering responses of C. ciliata to cold stress at the molecular level.

Alternative splicing of the receptor-like kinase Nt-Sd-RLK in tobacco cells responding to lipopolysaccharides: suggestive of a role in pathogen surveillance and perception?

Nt-Sd-RLK encodes an S-domain lectin receptor-like kinase that is induced in response to microbe-associated molecular pattern molecules (MAMPs) such as lipopolysaccharide (LPS). In this study, we investigated the alternative splicing of Nt-Sd-RLK in response to LPS stimulation. Our data indicate that in nonstimulated cells, a shorter transcript of Nt-Sd-RLK is generated and that in response to LPS, alternative splicing produces the full-length transcript. We propose that the extracellular domain of Nt-Sd-RLK encoded by the shorter transcript functions in pathogen surveillance. Once this domain binds LPS, alternative splicing generates the kinase domain-containing Nt-Sd-RLK that activates downstream signalling leading to a defence response. Thus, our findings suggest that plant defence signalling may be regulated through the alternative splicing of receptor-like kinases involved in pathogen recognition.

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.

Computational challenges, tools, and resources for analyzing co- and post-transcriptional events in high throughput

Co- and post-transcriptional regulation of gene expression is complex and multifaceted, spanning the complete RNA lifecycle from genesis to decay. High-throughput profiling of the constituent events and processes is achieved through a range of technologies that continue to expand and evolve. Fully leveraging the resulting data is nontrivial, and requires the use of computational methods and tools carefully crafted for specific data sources and often intended to probe particular biological processes. Drawing upon databases of information pre-compiled by other researchers can further elevate analyses. Within this review, we describe the major co- and post-transcriptional events in the RNA lifecycle that are amenable to high-throughput profiling. We place specific emphasis on the analysis of the resulting data, in particular the computational tools and resources available, as well as looking toward future challenges that remain to be addressed.

Stress-induced endogenous siRNAs targeting regulatory intron sequences in Brachypodium

Exposure to abiotic stresses triggers global changes in the expression of thousands of eukaryotic genes at the transcriptional and post-transcriptional levels. Small RNA (smRNA) pathways and splicing both function as crucial mechanisms regulating stress-responsive gene expression. However, examples of smRNAs regulating gene expression remain largely limited to effects on mRNA stability, translation, and epigenetic regulation. Also, our understanding of the networks controlling plant gene expression in response to environmental changes, and examples of these regulatory pathways intersecting, remains limited. Here, to investigate the role of smRNAs in stress responses we examined smRNA transcriptomes of Brachypodium distachyon plants subjected to various abiotic stresses. We found that exposure to different abiotic stresses specifically induced a group of novel, endogenous small interfering RNAs (stress-induced, UTR-derived siRNAs, or sutr-siRNAs) that originate from the 3' UTRs of a subset of coding genes. Our bioinformatics analyses predicted that sutr-siRNAs have potential regulatory functions and that over 90% of sutr-siRNAs target intronic regions of many mRNAs in trans. Importantly, a subgroup of these sutr-siRNAs target the important intron regulatory regions, such as branch point sequences, that could affect splicing. Our study indicates that in Brachypodium, sutr-siRNAs may affect splicing by masking or changing accessibility of specific cis-elements through base-pairing interactions to mediate gene expression in response to stresses. We hypothesize that this mode of regulation of gene expression may also serve as a general mechanism for regulation of gene expression in plants and potentially in other eukaryotes.

Genome-wide analysis of alternative splicing in Volvox carteri

Background

Alternative splicing is an essential mechanism for increasing transcriptome and proteome diversity in eukaryotes. Particularly in multicellular eukaryotes, this mechanism is involved in the regulation of developmental and physiological processes like growth, differentiation and signal transduction.

Results

Here we report the genome-wide analysis of alternative splicing in the multicellular green alga Volvox carteri. The bioinformatic analysis of 132,038 expressed sequence tags (ESTs) identified 580 alternative splicing events in a total of 426 genes. The predominant type of alternative splicing in Volvox is intron retention (46.5%) followed by alternative 5′ (17.9%) and 3′ (21.9%) splice sites and exon skipping (9.5%). Our analysis shows that in Volvox at least ~2.9% of the intron-containing genes are subject to alternative splicing. Considering the total number of sequenced ESTs, the Volvox genome seems to provide more favorable conditions (e.g., regarding length and GC content of introns) for the occurrence of alternative splicing than the genome of its close unicellular relative Chlamydomonas. Moreover, many randomly chosen alternatively spliced genes of Volvox do not show alternative splicing in Chlamydomonas. Since the Volvox genome contains about the same number of protein-coding genes as the Chlamydomonas genome (~14,500 protein-coding genes), we assumed that alternative splicing may play a key role in generation of genomic diversity, which is required to evolve from a simple one-cell ancestor to a multicellular organism with differentiated cell types (Mol Biol Evol 31:1402-1413, 2014). To confirm the alternative splicing events identified by bioinformatic analysis, several genes with different types of alternatively splicing have been selected followed by experimental verification of the predicted splice variants by RT-PCR.

Conclusions

The results show that our approach for prediction of alternative splicing events in Volvox was accurate and reliable. Moreover, quantitative real-time RT-PCR appears to be useful in Volvox for analyses of relationships between the appearance of specific alternative splicing variants and different kinds of physiological, metabolic and developmental processes as well as responses to environmental changes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1117) contains supplementary material, which is available to authorized users.

Spatial analysis of gene regulation reveals new insights into the molecular basis of upper thermal limits

The cellular stress response has long been the primary model for studying the molecular basis of thermal adaptation, yet the link between gene expression, RNA metabolism and physiological responses to thermal stress remains largely unexplored. We address this by comparing the transcriptional and physiological responses of three geographically distinct populations of Drosophila melanogaster from eastern Australia in response to, and recovery from, a severe heat stress with and without a prestress hardening treatment. We focus on starvin (stv), recently identified as an important thermally responsive gene. Intriguingly, stv encodes seven transcripts from alternative transcription sites and alternative splicing, yet appears to be rapidly heat inducible. First, we show genetic differences in upper thermal limits of the populations tested. We then demonstrate that the stv locus does not ubiquitously respond to thermal stress but is expressed as three distinct thermal and temporal RNA phenotypes (isoforms). The shorter transcript isoforms are rapidly upregulated under stress in all populations and show similar molecular signatures to heat-shock proteins. Multiple stress exposures seem to generate a reserve of pre-mRNAs, effectively 'priming' the cells for subsequent stress. Remarkably, we demonstrate a bypass in the splicing blockade in these isoforms, suggesting an essential role for these transcripts under heat stress. Temporal profiles for the weakly heat responsive stv isoform subset show opposing patterns in the two most divergent populations. Innate and induced transcriptome responses to hyperthermia are complex, and warrant moving beyond gene-level analyses.

Abscisic acid (ABA) regulation of Arabidopsis SR protein gene expression

Serine/arginine-rich (SR) proteins are major modulators of alternative splicing, a key generator of proteomic diversity and flexible means of regulating gene expression likely to be crucial in plant environmental responses. Indeed, mounting evidence implicates splicing factors in signal transduction of the abscisic acid (ABA) phytohormone, which plays pivotal roles in the response to various abiotic stresses. Using real-time RT-qPCR, we analyzed total steady-state transcript levels of the 18 SR and two SR-like genes from Arabidopsis thaliana in seedlings treated with ABA and in genetic backgrounds with altered expression of the ABA-biosynthesis ABA2 and the ABA-signaling ABI1 and ABI4 genes. We also searched for ABA-responsive cis elements in the upstream regions of the 20 genes. We found that members of the plant-specific SC35-Like (SCL) Arabidopsis SR protein subfamily are distinctively responsive to exogenous ABA, while the expression of seven SR and SR-related genes is affected by alterations in key components of the ABA pathway. Finally, despite pervasiveness of established ABA-responsive promoter elements in Arabidopsis SR and SR-like genes, their expression is likely governed by additional, yet unidentified cis-acting elements. Overall, this study pinpoints SR34, SR34b, SCL30a, SCL28, SCL33, RS40, SR45 and SR45a as promising candidates for involvement in ABA-mediated stress responses.

Alternative splicing in plant immunity

Alternative splicing (AS) occurs widely in plants and can provide the main source of transcriptome and proteome diversity in an organism. AS functions in a range of physiological processes, including plant disease resistance, but its biological roles and functional mechanisms remain poorly understood. Many plant disease resistance (R) genes undergo AS, and several R genes require alternatively spliced transcripts to produce R proteins that can specifically recognize pathogen invasion. In the finely-tuned process of R protein activation, the truncated isoforms generated by AS may participate in plant disease resistance either by suppressing the negative regulation of initiation of immunity, or by directly engaging in effector-triggered signaling. Although emerging research has shown the functional significance of AS in plant biotic stress responses, many aspects of this topic remain to be understood. Several interesting issues surrounding the AS of R genes, especially regarding its functional roles and regulation, will require innovative techniques and additional research to unravel.

Transcriptomic characterization of cold acclimation in larval zebrafish

Background

Temperature is one of key environmental parameters that affect the whole life of fishes and an increasing number of studies have been directed towards understanding the mechanisms of cold acclimation in fish. However, the adaptation of larvae to cold stress and the cold-specific transcriptional alterations in fish larvae remain largely unknown. In this study, we characterized the development of cold-tolerance in zebrafish larvae and investigated the transcriptional profiles under cold stress using RNA-seq.

Results

Pre-exposure of 96 hpf zebrafish larvae to cold stress (16°C) for 24 h significantly increased their survival rates under severe cold stress (12°C). RNA-seq generated 272 million raw reads from six sequencing libraries and about 92% of the processed reads were mapped to the reference genome of zebrafish. Differential expression analysis identified 1,431 up- and 399 down-regulated genes. Gene ontology enrichment analysis of cold-induced genes revealed that RNA splicing, ribosome biogenesis and protein catabolic process were the most highly overrepresented biological processes. Spliceosome, proteasome, eukaryotic ribosome biogenesis and RNA transport were the most highly enriched pathways for genes up-regulated by cold stress. Moreover, alternative splicing of 197 genes and promoter switching of 64 genes were found to be regulated by cold stress. A shorter isoform of stk16 that lacks 67 amino acids at the N-terminus was specifically generated by skipping the second exon in cold-treated larvae. Alternative promoter usage was detected for per3 gene under cold stress, which leading to a highly up-regulated transcript encoding a truncated protein lacking the C-terminal domains.

Conclusions

These findings indicate that zebrafish larvae possess the ability to build cold-tolerance under mild low temperature and transcriptional and post-transcriptional regulations are extensively involved in this acclimation process.

New levels of transcriptome complexity at upper thermal limits in wild Drosophila revealed by exon expression analysis

While the cellular heat-shock response has been a paradigm for studying the impact of thermal stress on RNA metabolism and gene expression, the genome-wide response to thermal stress and its connection to physiological stress resistance remain largely unexplored. Here, we address this issue using an array-based exon expression analysis to interrogate the transcriptome in recently established Drosophila melanogaster stocks during severe thermal stress and recovery. We first demonstrated the efficacy of exon-level analyses to reveal a level of thermally induced transcriptome complexity extending well beyond gene-level analyses. Next, we showed that the upper range of both the cellular and physiological thermal stress response profoundly affected message expression and processing in D. melanogaster, limiting expression to a small subset of transcripts, many that share features of known rapidly responding stress genes. As predicted from cellular heat-shock research, constitutive splicing was blocked in a set of novel genes; we did not detect changes to alternative splicing during heat stress, but rather induction of intronless isoforms of known heat-responsive genes. We observed transcriptome plasticity in the form of differential isoform expression during recovery from heat shock, mediated by multiple mechanisms including alternative transcription and alternative splicing. This affected genes involved in DNA regulation, immune response, and thermotolerance. These patterns highlight the complex nature of innate transcriptome responses under stress and potential for adaptive shifts through plasticity and evolved genetic responses at different hierarchical levels.

On the physiological significance of alternative splicing events in higher plants

Alternative splicing, which generates multiple transcripts from the same gene and potentially different protein isoforms, is a key posttranscriptional regulatory mechanism for expanding proteomic diversity and functional complexity in higher eukaryotes. The most recent estimates, based on whole transcriptome sequencing, indicate that about 95 % of human and 60 % of Arabidopsis multi-exon genes undergo alternative splicing, suggesting important roles for this mechanism in biological processes. However, while the misregulation of alternative splicing has been associated with many human diseases, its biological relevance in plant systems is just beginning to unfold. We review here the few plant genes for which the production of multiple splice isoforms has been reported to have a clear in vivo functional impact. These case studies implicate alternative splicing in the control of a wide range of physiological and developmental processes, including photosynthetic and starch metabolism, hormone signaling, seed germination, root growth and flowering, as well as in biotic and abiotic stress responses. Future functional characterization of alternative splicing events and identification of the transcripts targeted by major regulators of this versatile means of modulating gene expression should uncover the breadth of its physiological significance in higher plants.

The role of canonical and noncanonical pre-mRNA splicing in plant stress responses

Plants are sessile organisms capable of adapting to various environmental constraints, such as high or low temperatures, drought, soil salinity, or pathogen attack. To survive the unfavorable conditions, plants actively employ pre-mRNA splicing as a mechanism to regulate expression of stress-responsive genes and reprogram intracellular regulatory networks. There is a growing evidence that various stresses strongly affect the frequency and diversity of alternative splicing events in the stress-responsive genes and lead to an increased accumulation of mRNAs containing premature stop codons, which in turn have an impact on plant stress response. A number of studies revealed that some mRNAs involved in plant stress response are spliced counter to the traditional conception of alternative splicing. Such noncanonical mRNA splicing events include trans-splicing, intraexonic deletions, or variations affecting multiple exons and often require short direct repeats to occur. The noncanonical alternative splicing, along with common splicing events, targets the spliced transcripts to degradation through nonsense-mediated mRNA decay or leads to translation of truncated proteins. Investigation of the diversity, biological consequences, and mechanisms of the canonical and noncanonical alternative splicing events will help one to identify those transcripts which are promising for using in genetic engineering and selection of stress-tolerant plants.

Transgene regulation in plants by alternative splicing of a suicide exon

Compared to transcriptional activation, other mechanisms of gene regulation have not been widely exploited for the control of transgenes. One barrier to the general use and application of alternative splicing is that splicing-regulated transgenes have not been shown to be reliably and simply designed. Here, we demonstrate that a cassette bearing a suicide exon can be inserted into a variety of open reading frames (ORFs), generating transgenes whose expression is activated by exon skipping in response to a specific protein inducer. The surprisingly minimal sequence requirements for the maintenance of splicing fidelity and regulation indicate that this splicing cassette can be used to regulate any ORF containing one of the amino acids Glu, Gln or Lys. Furthermore, a single copy of the splicing cassette was optimized by rational design to confer robust gene activation with no background expression in plants. Thus, conditional splicing has the potential to be generally useful for transgene regulation.

Transcriptome survey reveals increased complexity of the alternative splicing landscape in Arabidopsis

Alternative splicing (AS) is a key regulatory mechanism that contributes to transcriptome and proteome diversity. As very few genome-wide studies analyzing AS in plants are available, we have performed high-throughput sequencing of a normalized cDNA library which resulted in a high coverage transcriptome map of Arabidopsis. We detect ∼150,000 splice junctions derived mostly from typical plant introns, including an eightfold increase in the number of U12 introns (2069). Around 61% of multiexonic genes are alternatively spliced under normal growth conditions. Moreover, we provide experimental validation of 540 AS transcripts (from 256 genes coding for important regulatory factors) using high-resolution RT-PCR and Sanger sequencing. Intron retention (IR) is the most frequent AS event (∼40%), but many IRs have relatively low read coverage and are less well-represented in assembled transcripts. Additionally, ∼51% of Arabidopsis genes produce AS transcripts which do not involve IR. Therefore, the significance of IR in generating transcript diversity was generally overestimated in previous assessments. IR analysis allowed the identification of a large set of cryptic introns inside annotated coding exons. Importantly, a significant fraction of these cryptic introns are spliced out in frame, indicating a role in protein diversity. Furthermore, we show extensive AS coupled to nonsense-mediated decay in AFC2, encoding a highly conserved LAMMER kinase which phosphorylates splicing factors, thus establishing a complex loop in AS regulation. We provide the most comprehensive analysis of AS to date which will serve as a valuable resource for the plant community to study transcriptome complexity and gene regulation.

A stress responsive alternative splicing mechanism in Citrus clementina leaves

Chitinases are often considered pathogenesis-related proteins since their activity can be induced by viral infections, fungal and bacterial cell wall components, and also by more general sources of stress such as wounding, salicylic acid, ethylene, auxins and cytokinins. In the present study, comparative proteomic analysis showed the defense-related acidic chitinase II to be specifically induced in Citrus clementina leaves infested by the two-spotted spider mite Tetranychus urticae or treated with MeJA. In parallel, changes in the mRNA profiles of two partially homologous chitinase forms were shown by RT-PCR. In particular, the appearance of an additional cDNA chitinase fragment in T. urticae-infested and MeJA-treated leaves was observed. This finding may indicate a specific regulatory mechanism of chitinase expression. We report evidence for alternative splicing in T. urticae-infested C. clementina, where a premature stop codon after the first 135 amino acids was introduced. We observed inducible chitinase activity after MeJA treatment, indicative of a rapid plant response to infestation. This work provides the first evidence of chitinase alternative splicing in C. clementina. In addition, the presence of the dual-band pattern for chitinase cDNA by RT-PCR may represent a suitable predictive marker for early diagnosis of plant biotic stress.

Genome-wide analysis of alternative splicing in Chlamydomonas reinhardtii

Background

Genome-wide computational analysis of alternative splicing (AS) in several flowering plants has revealed that pre-mRNAs from about 30% of genes undergo AS. Chlamydomonas, a simple unicellular green alga, is part of the lineage that includes land plants. However, it diverged from land plants about one billion years ago. Hence, it serves as a good model system to study alternative splicing in early photosynthetic eukaryotes, to obtain insights into the evolution of this process in plants, and to compare splicing in simple unicellular photosynthetic and non-photosynthetic eukaryotes. We performed a global analysis of alternative splicing in Chlamydomonas reinhardtii using its recently completed genome sequence and all available ESTs and cDNAs.

Results

Our analysis of AS using BLAT and a modified version of the Sircah tool revealed AS of 498 transcriptional units with 611 events, representing about 3% of the total number of genes. As in land plants, intron retention is the most prevalent form of AS. Retained introns and skipped exons tend to be shorter than their counterparts in constitutively spliced genes. The splice site signals in all types of AS events are weaker than those in constitutively spliced genes. Furthermore, in alternatively spliced genes, the prevalent splice form has a stronger splice site signal than the non-prevalent form. Analysis of constitutively spliced introns revealed an over-abundance of motifs with simple repetitive elements in comparison to introns involved in intron retention. In almost all cases, AS results in a truncated ORF, leading to a coding sequence that is around 50% shorter than the prevalent splice form. Using RT-PCR we verified AS of two genes and show that they produce more isoforms than indicated by EST data. All cDNA/EST alignments and splice graphs are provided in a website at http://combi.cs.colostate.edu/as/chlamy.

Conclusions

The extent of AS in Chlamydomonas that we observed is much smaller than observed in land plants, but is much higher than in simple unicellular heterotrophic eukaryotes. The percentage of different alternative splicing events is similar to flowering plants. Prevalence of constitutive and alternative splicing in Chlamydomonas, together with its simplicity, many available public resources, and well developed genetic and molecular tools for this organism make it an excellent model system to elucidate the mechanisms involved in regulated splicing in photosynthetic eukaryotes.

Characterization of the Sesbania rostrata phytochelatin synthase gene: alternative splicing and function of four isoforms

Phytochelatins (PCs) play an important role in detoxification of heavy metals in plants. PCs are synthesized from glutathione by phytochelatin synthase (PCS), a dipeptidyltransferase. Sesbania rostrata is a tropical legume plant that can tolerate high concentrations of Cd and Zn. In this study, the S. rostrata PCS gene (SrPCS) and cDNAs were isolated and characterized. Southern blot and sequence analysis revealed that a single copy of the SrPCS gene occurs in the S. rostrata genome, and produces four different SrPCS mRNAs and proteins, SrPCS1–SrPCS4, by alternative splicing of the SrPCS pre-mRNA. The SrPCS1 and SrPCS3 proteins conferred Cd tolerance when expressed in yeast cells, whereas the SrPCS2 and SrPCS4 proteins, which lack the catalytic triad and the N-terminal domains, did not. These results suggested that SrPCS1 and SrPCS3 have potential applications in genetic engineering of plants for enhancing heavy metal tolerance and phytoremediation of contaminated soils.