Genome-wide analysis of alternative splicing of pre-mRNA under salt stress in Arabidopsis

Background

Alternative splicing (AS) of precursor mRNA (pre-mRNA) is an important gene regulation process that potentially regulates many physiological processes in plants, including the response to abiotic stresses such as salt stress.

Results

To analyze global changes in AS under salt stress, we obtained high-coverage (~200 times) RNA sequencing data from Arabidopsis thaliana seedlings that were treated with different concentrations of NaCl. We detected that ~49% of all intron-containing genes were alternatively spliced under salt stress, 10% of which experienced significant differential alternative splicing (DAS). Furthermore, AS increased significantly under salt stress compared with under unstressed conditions. We demonstrated that most DAS genes were not differentially regulated by salt stress, suggesting that AS may represent an independent layer of gene regulation in response to stress. Our analysis of functional categories suggested that DAS genes were associated with specific functional pathways, such as the pathways for the responses to stresses and RNA splicing. We revealed that serine/arginine-rich (SR) splicing factors were frequently and specifically regulated in AS under salt stresses, suggesting a complex loop in AS regulation for stress adaptation. We also showed that alternative splicing site selection (SS) occurred most frequently at 4 nucleotides upstream or downstream of the dominant sites and that exon skipping tended to link with alternative SS.

Conclusions

Our study provided a comprehensive view of AS under salt stress and revealed novel insights into the potential roles of AS in plant response to salt stress.

Electronic supplementary material

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

Exploring Splicing Variants and Novel Genes in Sacred Lotus Based on RNA-seq Data

Sacred lotus (Nelumbo nucifera) is a typical aquatic plant, belonging to basal eudicot plant, which is ideal for genome and genetic evolutionary study. Understanding lotus gene diversity is important for the study of molecular genetics and breeding. In this research, public RNA-seq data and the annotated reference genome were used to identify the genes in lotus. A total of 26,819 consensus and 1,081 novel genes were identified. Meanwhile, a comprehensive analysis of gene alternative splicing events was conducted, and a total of 19,983 “internal” alternative splicing (AS) events and 14,070 “complete” AS events were detected in 5,878 and 5,881 multi-exon expression genes, respectively. Observations made from the AS events show the predominance of intron retention (IR) subtype of AS events representing 33%. IR is followed by alternative acceptor (AltA), alternative donor (AltD) and exon skipping (ES), highlighting the universality of the intron definition model in plants. In addition, functional annotations of the gene with AS indicated its relationship to a number of biological processes such as cellular process and metabolic process, showing the key role for alternative splicing in influencing the growth and development of lotus. The results contribute to a better understanding of the current gene diversity in lotus, and provide an abundant resource for future functional genome analysis in lotus.     

Identification,characterization and molecular phylogeny of U12-dependent introns in the Arabidopsis thaliana genome

U12-dependent introns are spliced by the minor U12-type spliceosome and occur in a variety of eukaryotic organisms, including Arabidopsis. In this study, a set of putative U12-dependent introns was compiled from a large collection of cDNA/EST- confirmed introns in the Arabidopsis thaliana genome by means of high-throughput bioinformatic analysis combined with manual scrutiny. A total of 165 U12-type introns were identified based upon stringent criteria. This number of sequences well exceeds the total number of U12-type introns previously reported for plants and allows a more thorough statistical analysis of U12-type signals. Of particular note is the discovery that the distance between the branch site adenosine and the acceptor site ranges from 10 to 39 nt, significantly longer than the previously postulated limit of 21 bp. Further analysis indicates that, in addition to the spacing constraint, the sequence context of the potential acceptor site may have an important role in 3′ splice site selection. Several alternative splicing events involving U12-type introns were also captured in this study, providing evidence that U12-dependent acceptor sites can also be recognized by the U2-type spliceosome. Furthermore, phylogenetic analysis suggests that both U12-type AT-AC and U12-type GT-AG introns occurred in Na⁺/H⁺ antiporters in a progenitor of animals and plants.     

Comprehensive splice-site analysis using comparative genomics

We have collected over half a million splice sites from five species-Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans and Arabidopsis thaliana-and classified them into four subtypes: U2-type GT-AG and GC-AG and U12-type GT-AG and AT-AC. We have also found new examples of rare splice-site categories, such as U12-type introns without canonical borders, and U2-dependent AT-AC introns. The splice-site sequences and several tools to explore them are available on a public website (SpliceRack). For the U12-type introns, we find several features conserved across species, as well as a clustering of these introns on genes. Using the information content of the splice-site motifs, and the phylogenetic distance between them, we identify: (i) a higher degree of conservation in the exonic portion of the U2-type splice sites in more complex organisms; (ii) conservation of exonic nucleotides for U12-type splice sites; (iii) divergent evolution of C.elegans 3' splice sites (3'ss) and (iv) distinct evolutionary histories of 5' and 3'ss. Our study proves that the identification of broad patterns in naturally-occurring splice sites, through the analysis of genomic datasets, provides mechanistic and evolutionary insights into pre-mRNA splicing.     

Alternative splicing and bioinformatic analysis of human U12-type introns

U12-type introns exist, albeit rarely, in a variety of multicellular organisms. Splicing of U12 intron-containing precursor mRNAs takes place in the U12-type spliceosome that is distinct from the major U2-type spliceosome. Due to incompatibility of these two spliceosomes, alternative splicing involving a U12-type intron may give rise to a relatively complicated impact on gene expression. We studied alternative U12-type intron splicing in an attempt to gain more mechanistic insights. First, we characterized mutually exclusive exon selection of the human JNK2 gene, which involves an unusual intron possessing the U12-type 5′ splice site and the U2-type 3′ splice site. We demonstrated that the long and evolutionary conserved polypyrimidine tract of this hybrid intron provides important signals for inclusion of its downstream alternative exon. In addition, we examined the effects of single nucleotide polymorphisms in the human WDFY1 U12-type intron on pre-mRNA splicing. These results provide mechanistic implications on splice-site selection of U12-type intron splicing. We finally discuss the potential effects of splicing of a U12-type intron with genetic defects or within a set of genes encoding RNA processing factors on global gene expression.     

水稻NBS-LRR基因选择性剪接的全基因组检测及分析

选择性剪接是促进基因组复杂性和蛋白质组多样性的一种主要机制,但是对水稻NBS-LRR序列选择性剪接的全基因组分析却未见报道。通过隐马尔柯夫模型搜索,从TIGR数据库里得到了855条编码NBS-LRR基序的序列。利用这些序列在KOME、TIGR基因索引及UniProt三个数据库中进行同源搜索,获得同源的完整cDNA序列、假设一致性序列和蛋白质序列。再利用Spidey和SIM4程序把完整cDNA序列和假设一致性序列联配到相应的BAC序列上来预测选择性剪接。蛋白质序列和基因组序列之间的联配使用tBLASTn。在这875个NBS-LRR基因中,119个基因具有选择性剪接现象,其中包括71内含子保留,20个外显子跳跃,25个选择性起始,16个选择性终止,12个5′端的选择性剪接和16个3′端选择性剪接。大多数选择性剪接都为两个和多个转录本所支持。可以通过访问http://www.bioinfor.org查询这些数据。进而通过生物信息学分析剪接边界发现外显子跳跃和内含子保留的‘GT…AG’的规则不如组成型的保守。这暗示了它们是通过不同的调控机制来指导剪接变构体的形成。通过分析内含子保留对蛋白质的影响,发现选择性剪接的蛋白更倾向于改变其C端氨基酸序列。最后对选择性剪接的组织分布和蛋白质定位进行分析,结果表明选择性剪接的最大类的组织分布是根和愈伤组织。超过1/3剪接变构体的蛋白质定位是质膜和细胞质。这些选择性剪接蛋白可能在抗病信号转导中起到重要作用。     

基于RNA-Seq数据识别果蝇剪接位点和可变剪接事件

完整基因结构的预测是当前生命科学研究的一个重要基础课题,其中一个关键环节是剪接位点和各种可变剪接事件的精确识别.基于转录组测序（RNA-seq）数据,识别剪接位点和可变剪接事件是近几年随着新一代测序技术发展起来的新技术策略和方法.本工作基于黑腹果蝇睾丸RNA-seq数据,使用TopHat软件成功识别出39718个果蝇剪接位点,其中有10584个新剪接位点.同时,基于剪接位点的不同组合,针对各类型可变剪接特征开发出计算识别算法,成功识别了8477个可变剪接事件（其中新识别的可变剪接事件3922个）,包括可变供体位点、可变受体位点、内含子保留和外显子缺失4种类型.RT-PCR实验验证了2个果蝇基因上新识别的可变剪接事件,发现了全新的剪接异构体.进一步表明,RNA-seq数据可有效应用于识别剪接位点和可变剪接事件,为深入揭示剪接机制及可变剪接生物学功能提供新思路和新手段.     

A computational scan for U12-dependent introns in the human genome sequence

U12-dependent introns are found in small numbers in most eukaryotic genomes, but their scarcity makes accurate characterisation of their properties challenging. A computational search for U12-dependent introns was performed using the draft version of the human genome sequence. Human expressed sequences confirmed 404 U12-dependent introns within the human genome, a 6-fold increase over the total number of non-redundant U12-dependent introns previously identified in all genomes. Although most of these introns had AT-AC or GT-AG terminal dinucleotides, small numbers of introns with a surprising diversity of termini were found, suggesting that many of the non-canonical introns found in the human genome may be variants of U12-dependent introns and, thus, spliced by the minor spliceosome. Comparisons with U2-dependent introns revealed that the U12-dependent intron set lacks the ‘short intron’ peak characteristic of U2-dependent introns. Analysis of this U12-dependent intron set confirmed reports of a biased distribution of U12-dependent introns in the genome and allowed the identification of several alternative splicing events as well as a surprising number of apparent splicing errors. This new larger reference set of U12-dependent introns will serve as a resource for future studies of both the properties and evolution of the U12 spliceosome.     

An LKB1 AT-AC intron mutation causes Peutz-Jeghers syndrome via splicing at noncanonical cryptic splice sites

Peutz-Jeghers syndrome (PJS) is an autosomal dominant disorder associated with gastrointestinal polyposis and an increased cancer risk. PJS is caused by germline mutations in the tumor suppressor gene LKB1. One such mutation, IVS2+1A>G, alters the second intron 5' splice site, which has sequence features of a U12-type AT-AC intron. We report that in patients, LKB1 RNA splicing occurs from the mutated 5' splice site to several cryptic, noncanonical 3' splice sites immediately adjacent to the normal 3' splice site. In vitro splicing analysis demonstrates that this aberrant splicing is mediated by the U12-dependent spliceosome. The results indicate that the minor spliceosome can use a variety of 3' splice site sequences to pair to a given 5' splice site, albeit with tight constraints for maintaining the 3' splice site position. The unusual splicing defect associated with this PJS-causing mutation uncovers differences in splice-site recognition between the major and minor pre-mRNA splicing pathways.     

An intronic enhancer regulates splicing of the twintron of Drosophila melanogaster prospero pre-mRNA by two different spliceosomes

We have examined the alternative splicing of the Drosophila melanogaster prospero twintron, which contains splice sites for both the U2- and U12-type spliceosome and generates two forms of mRNA, pros-L (U2-type product) and pros-S (U12-type product). We find that twintron splicing is developmentally regulated: pros-L is abundant in early embryogenesis while pros-S displays the opposite pattern. We have established a Kc cell in vitro splicing system that accurately splices a minimal pros substrate containing the twintron and have examined the sequence requirements for pros twintron splicing. Systematic deletion and mutation analysis of intron sequences established that twintron splicing requires a 46-nucleotide purine-rich element located 32 nucleotides downstream of the U2-type 5' splice site. While this element regulates both splicing pathways, its alteration showed the severest effects on the U2-type splicing pathway. Addition of an RNA competitor containing the wild-type purine-rich element to the Kc extract abolished U2-type splicing and slightly repressed U12-type splicing, suggesting that a trans-acting factor(s) binds the enhancer element to stimulate twintron splicing. Thus, we have identified an intron region critical for prospero twintron splicing as a first step towards elucidating the molecular mechanism of splicing regulation involving competition between the two kinds of spliceosomes.     

Splicing and alternative splicing in rice and humans

Rice is a monocot gramineous crop, and one of the most important staple foods. Rice is considered a model species for most gramineous crops. Extensive research on rice has provided critical guidance for other crops, such as maize and wheat. In recent years, climate change and exacerbated soil degradation have resulted in a variety of abiotic stresses, such as greenhouse effects, lower temperatures, drought, floods, soil salinization and heavy metal pollution. As such, there is an extremely high demand for additional research, in order to address these negative factors. Studies have shown that the alternative splicing of many genes in rice is affected by stress conditions, suggesting that manipulation of the alternative splicing of specific genes may be an effective approach for rice to adapt to abiotic stress. With the advancement of microarrays, and more recently, next generation sequencing technology, several studies have shown that more than half of the genes in the rice genome undergo alternative splicing. This mini-review summarizes the latest progress in the research of splicing and alternative splicing in rice, compared to splicing in humans. Furthermore, we discuss how additional studies may change the landscape of investigation of rice functional genomics and genetically improved rice. [BMB Reports 2013; 46(9): 439-447]