Splicing and alternative splicing in rice and humans

An alternative 3' splice site of PeuHKT1;3 improves the response to salt stress through enhancing affinity to K+ in Populus

Alternative splicing (AS) serves as a crucial post-transcriptional regulator in plants that contributes to the resistance to salt stress. However, the underlying mechanism is largely unknown. In this research, we identified an important AS transcript in Populus euphratica, PeuHKT1:3a, generated by alternative 3' splice site splicing mode that resulted in the removal of 252 bases at the 5' end of the first exon in PeuHKT1:3. Protein sequence comparison showed that the site of AS occurred in PeuHKT1:3 is located at a crucial Ser residue within the first pore-loop domain, which leads to inefficient K+ transport in HKT I-type transporters. Expressing PeuHKT1;3a in an axt3 mutant yeast strain can effectively compensate for the lack of intracellular K+, whereas the expression of PeuHKT1;3 cannot yield the effect. Furthermore, in transgenic Arabidopsis and poplar plants, it was observed that lines expressing PeuHKT1;3a exhibited greater salt tolerance compared to those expressing the PeuHKT1;3 strain. Analysis of ion content and flux demonstrated that the transgenic PeuHKT1;3a line exhibited significantly higher K+ content compared to the PeuHKT1;3 line, while there was no significant difference in Na+ content. In conclusion, our findings revealed that AS can give rise to novel variants of HKT I-type proteins in P. euphratica with modified K+ selectivity to keep a higher K+/Na+ ratio to enhanced salt tolerance.

Nucleotide-level distance metrics to quantify alternative splicing implemented in TranD

Advances in affordable transcriptome sequencing combined with better exon and gene prediction has motivated many to compare transcription across the tree of life. We develop a mathematical framework to calculate complexity and compare transcript models. Structural features, i.e. intron retention (IR), donor/acceptor site variation, alternative exon cassettes, alternative 5'/3' UTRs, are compared and the distance between transcript models is calculated with nucleotide level precision. All metrics are implemented in a PyPi package, TranD and output can be used to summarize splicing patterns for a transcriptome (1GTF) and between transcriptomes (2GTF). TranD output enables quantitative comparisons between: annotations augmented by empirical RNA-seq data and the original transcript models; transcript model prediction tools for longread RNA-seq (e.g. FLAIR versus Isoseq3); alternate annotations for a species (e.g. RefSeq vs Ensembl); and between closely related species. In C. elegans, Z. mays, D. melanogaster, D. simulans and H. sapiens, alternative exons were observed more frequently in combination with an alternative donor/acceptor than alone. Transcript models in RefSeq and Ensembl are linked and both have unique transcript models with empirical support. D. melanogaster and D. simulans, share many transcript models and long-read RNAseq data suggests that both species are under-annotated. We recommend combined references.

Advances in alternative splicing identification: deep learning and pantranscriptome

In plants, alternative splicing is a crucial mechanism for regulating gene expression at the post-transcriptional level, which leads to diverse proteins by generating multiple mature mRNA isoforms and diversify the gene regulation. Due to the complexity and variability of this process, accurate identification of splicing events is a vital step in studying alternative splicing. This article presents the application of alternative splicing algorithms with or without reference genomes in plants, as well as the integration of advanced deep learning techniques for improved detection accuracy. In addition, we also discuss alternative splicing studies in the pan-genomic background and the usefulness of integrated strategies for fully profiling alternative splicing.

Prediction of transcript structure and concentration using RNA-Seq data

Alternative splicing (AS) is a key post-transcriptional modification that helps in increasing protein diversity. Almost 90% of the protein-coding genes in humans are known to undergo AS and code for different transcripts. Some transcripts are associated with diseases such as breast cancer, lung cancer and glioblastoma. Hence, these transcripts can serve as novel therapeutic and prognostic targets for drug discovery. Herein, we have developed a pipeline, Finding Alternative Splicing Events (FASE), as the R package that includes modules to determine the structure and concentration of transcripts using differential AS. To predict the correct structure of expressed transcripts in given conditions, FASE combines the AS events with the information of exons, introns and junctions using graph theory. The estimated concentration of predicted transcripts is reported as the relative expression in terms of log2CPM. Using FASE, we were able to identify several unique transcripts of EMILIN1 and SLK genes in the TCGA-BRCA data, which were validated using RT-PCR. The experimental study demonstrated consistent results, which signify the high accuracy and precision of the developed methods. In conclusion, the developed pipeline, FASE, can efficiently predict novel transcripts that are missed in general transcript-level differential expression analysis. It can be applied selectively from a single gene to simple or complex genome even in multiple experimental conditions for the identification of differential AS-based biomarkers, prognostic targets and novel therapeutics.

Alternative Splicing of TaHsfA2-7 Is Involved in the Improvement of Thermotolerance in Wheat

High temperature has severely affected plant growth and development, resulting in reduced production of crops worldwide, especially wheat. Alternative splicing (AS), a crucial post-transcriptional regulatory mechanism, is involved in the growth and development of eukaryotes and the adaptation to environmental changes. Previous transcriptome data suggested that heat shock transcription factor (Hsf) TaHsfA2-7 may form different transcripts by AS. However, it remains unclear whether this post-transcriptional regulatory mechanism of TaHsfA2-7 is related to thermotolerance in wheat (Triticum aestivum). Here, we identified a novel splice variant, TaHsfA2-7-AS, which was induced by high temperature and played a positive role in thermotolerance regulation in wheat. Moreover, TaHsfA2-7-AS is predicted to encode a small truncated TaHsfA2-7 isoform, retaining only part of the DNA-binding domain (DBD). TaHsfA2-7-AS is constitutively expressed in various tissues of wheat. Notably, the expression level of TaHsfA2-7-AS is significantly up-regulated by heat shock (HS) during flowering and grain-filling stages in wheat. Further studies showed that TaHsfA2-7-AS was localized in the nucleus but lacked transcriptional activation activity. Ectopic expression of TaHsfA2-7-AS in yeast exhibited improved thermotolerance. Compared to non-transgenic plants, overexpression of TaHsfA2-7-AS in Arabidopsis results in enhanced tolerance to heat stress. Simultaneously, we also found that TaHsfA1 is directly involved in the transcriptional regulation of TaHsfA2-7 and TaHsfA2-7-AS. In summary, our findings demonstrate the function of TaHsfA2-7-AS splicing variant in response to heat stress and establish a link between regulatory mechanisms of AS and the improvement of thermotolerance in wheat.

Alternative Splicing and Its Roles in Plant Metabolism

Plant metabolism, including primary metabolism such as tricarboxylic acid cycle, glycolysis, shikimate and amino acid pathways as well as specialized metabolism such as biosynthesis of phenolics, alkaloids and saponins, contributes to plant survival, growth, development and interactions with the environment. To this end, these metabolic processes are tightly and finely regulated transcriptionally, post-transcriptionally, translationally and post-translationally in response to different growth and developmental stages as well as the constantly changing environment. In this review, we summarize and describe the current knowledge of the regulation of plant metabolism by alternative splicing, a post-transcriptional regulatory mechanism that generates multiple protein isoforms from a single gene by using alternative splice sites during splicing. Numerous genes in plant metabolism have been shown to be alternatively spliced under different developmental stages and stress conditions. In particular, alternative splicing serves as a regulatory mechanism to fine-tune plant metabolism by altering biochemical activities, interaction and subcellular localization of proteins encoded by splice isoforms of various genes.

AS-Quant: Detection and Visualization of Alternative Splicing Events with RNA-seq Data

(1) Background: A simplistic understanding of the central dogma falls short in correlating the number of genes in the genome to the number of proteins in the proteome. Post-transcriptional alternative splicing contributes to the complexity of the proteome and is critical in understanding gene expression. mRNA-sequencing (RNA-seq) has been widely used to study the transcriptome and provides opportunity to detect alternative splicing events among different biological conditions. Despite the popularity of studying transcriptome variants with RNA-seq, few efficient and user-friendly bioinformatics tools have been developed for the genome-wide detection and visualization of alternative splicing events. (2) Results: We propose AS-Quant, (Alternative Splicing Quantitation), a robust program to identify alternative splicing events from RNA-seq data. We then extended AS-Quant to visualize the splicing events with short-read coverage plots along with complete gene annotation. The tool works in three major steps: (i) calculate the read coverage of the potential spliced exons and the corresponding gene; (ii) categorize the events into five different categories according to the annotation, and assess the significance of the events between two biological conditions; (iii) generate the short reads coverage plot for user specified splicing events. Our extensive experiments on simulated and real datasets demonstrate that AS-Quant outperforms the other three widely used baselines, SUPPA2, rMATS, and diffSplice for detecting alternative splicing events. Moreover, the significant alternative splicing events identified by AS-Quant between two biological contexts were validated by RT-PCR experiment. (3) Availability: AS-Quant is implemented in Python 3.0. Source code and a comprehensive user's manual are freely available online.

Genome-wide identification, evolution, expression, and alternative splicing profiles of peroxiredoxin genes in cotton

Peroxiredoxin (PRX) is a ubiquitous thioredoxin-dependent peroxidase that can eliminate excessive free radicals produced by stress and protect cells from oxidative damage. PRXs are also involved in reactive oxygen species (ROS)- and redox-dependent signaling by performing redox interactions with other proteins and modify their redox status. At present, PRX family identification, evolution and regulation research has been conducted in some plants; however, systematic research about this family is lacking in cotton. In this study, a total of 44 PRXs were identified in the cotton genome. Phylogenetic and conserved active site analyses showed that the PRXs were divided into six subfamilies according to the conserved site (PxxxTxxC…S…W/F) and conserved cysteinyl residues positions. Segmental duplication and polyploid events were the main methods for PRX family expansion, and the PRXs of diploid G. arboreum were the donors of PRXs in the D subgenomes of allotetraploid G. hirsutum and G. barbadense during the evolution of the PRX family. qRT-PCR analysis confirmed that cis-acting elements play important roles in regulating the expression of PRXs. Alternative splicing events occurred in GhPRX14-D that can increased the complexity of transcripts in G. hirsutum. Subcellular localization showed that most PRX members were located in chloroplasts, the cytoplasmic membrane and the nucleus. Our results provide systematic support for a better understanding of PRXs in cotton and a starting point for further studies of the specific functions of PRXs in cotton.

A genome-wide survey of copy number variations reveals an asymmetric evolution of duplicated genes in rice

Background

Copy number variations (CNVs) are an important type of structural variations in the genome that usually affect gene expression levels by gene dosage effect. Understanding CNVs as part of genome evolution may provide insights into the genetic basis of important agricultural traits and contribute to the crop breeding in the future. While available methods to detect CNVs utilizing next-generation sequencing technology have helped shed light on prevalence and effects of CNVs, the complexity of crop genomes poses a major challenge and requires development of additional tools.

Results

Here, we generated genomic and transcriptomic data of 93 rice (Oryza sativa L.) accessions and developed a comprehensive pipeline to call CNVs in this large-scale dataset. We analyzed the correlation between CNVs and gene expression levels and found that approximately 13% of the identified genes showed a significant correlation between their expression levels and copy numbers. Further analysis showed that about 36% of duplicate pairs were involved in pseudogenetic events while only 5% of them showed functional differentiation. Moreover, the offspring copy mainly contributed to the expression levels and seemed more likely to become a pseudogene, whereas the parent copy tended to maintain the function of ancestral gene.

Conclusion

We provide a high-accuracy CNV dataset that will contribute to functional genomics studies and molecular breeding in rice. We also showed that gene dosage effect of CNVs in rice is not exponential or linear. Our work demonstrates that the evolution of duplicated genes is asymmetric in both expression levels and gene fates, shedding a new insight into the evolution of duplicated genes.

A single mutation in the ACTR8 gene associated with lineage-specific expression in primates

Background

Alternative splicing (AS) generates various transcripts from a single gene and thus plays a significant role in transcriptomic diversity and proteomic complexity. Alu elements are primate-specific transposable elements (TEs) and can provide a donor or acceptor site for AS. In a study on TE-mediated AS, we recently identified a novel AluSz6-exonized ACTR8 transcript of the crab-eating monkey (Macaca fascicularis). In the present study, we sought to determine the molecular mechanism of AluSz6 exonization of the ACTR8 gene and investigate its evolutionary and functional consequences in the crab-eating monkey.

Results

We performed RT-PCR and genomic PCR to analyze AluSz6 exonization in the ACTR8 gene and the expression of the AluSz6-exonized transcript in nine primate samples, including prosimians, New world monkeys, Old world monkeys, and hominoids. AluSz6 integration was estimated to have occurred before the divergence of simians and prosimians. The Alu-exonized transcript obtained by AS was lineage-specific and expressed only in Old world monkeys and apes, and humans. This lineage-specific expression was caused by a single G duplication in AluSz6, which provides a new canonical 5′ splicing site. We further identified other alternative transcripts that were unaffected by the AluSz6 insertion. Finally, we observed that the alternative transcripts were transcribed into new isoforms with C-terminus deletion, and in silico analysis showed that these isoforms do not have a destructive function.

Conclusions

The single G duplication in the TE sequence is the source of TE exonization and AS, and this mutation may suffer a different fate of ACTR8 gene expression during primate evolution.

A large-scale comparative study of isoform expressions measured on four platforms

BACKGROUND

Most eukaryotic genes produce different transcripts of multiple isoforms by inclusion or exclusion of particular exons. The isoforms of a gene often play diverse functional roles, and thus it is necessary to accurately measure isoform expressions as well as gene expressions. While previous studies have demonstrated the strong agreement between mRNA sequencing (RNA-seq) and array-based gene and/or isoform quantification platforms (Microarray gene expression and Exon-array), the more recently developed NanoString platform has not been systematically evaluated and compared, especially in large-scale studies across different cancer domains.

RESULTS

In this paper, we present a large-scale comparative study among RNA-seq, NanoString, array-based, and RT-qPCR platforms using 46 cancer cell lines across different cancer types. The goal is to understand and evaluate the calibers of the platforms for measuring gene and isoform expressions in cancer studies. We first performed NanoString experiments on 59 cancer cell lines with 404 custom-designed probes for measuring the expressions of 478 isoforms in 155 genes, and additional RT-qPCR experiments for a subset of the measured isoforms in 13 cell lines. We then combined the data with the matched RNA-seq, Exon-array, and Microarray data of 46 of the 59 cell lines for the comparative analysis.

CONCLUSION

In the comparisons of the platforms for measuring the expressions at both isoform and gene levels, we found that (1) the agreement on isoform expressions is lower than the agreement on gene expressions across the four platforms; (2) NanoString and Exon-array are not consistent on isoform quantification even though both techniques are based on hybridization reactions; (3) RT-qPCR experiments are more consistent with RNA-seq and Exon-array than NanoString in isoform quantification; (4) different RNA-seq isoform quantification methods show varying estimation results, and among the methods, Net-RSTQ and eXpress are more consistent across the platforms; and (5) RNA-seq has the best overall consistency with the other platforms on gene expression quantification.

Insights into the transcriptional and post-transcriptional regulation of the rice SUMOylation machinery and into the role of two rice SUMO proteases

BACKGROUND

SUMOylation is an essential eukaryotic post-translation modification that, in plants, regulates numerous cellular processes, ranging from seed development to stress response. Using rice as a model crop plant, we searched for potential regulatory points that may influence the activity of the rice SUMOylation machinery genes.

RESULTS

We analyzed the presence of putative cis-acting regulatory elements (CREs) within the promoter regions of the rice SUMOylation machinery genes and found CREs related to different cellular processes, including hormone signaling. We confirmed that the transcript levels of genes involved in target-SUMOylation, containing ABA- and GA-related CREs, are responsive to treatments with these hormones. Transcriptional analysis in Nipponbare (spp. japonica) and LC-93-4 (spp. indica), showed that the transcript levels of all studied genes are maintained in the two subspecies, under normal growth. OsSUMO3 is an exceptional case since it is expressed at low levels or is not detectable at all in LC-93-4 roots and shoots, respectively. We revealed post-transcriptional regulation by alternative splicing (AS) for all genes studied, except for SUMO coding genes, OsSIZ2, OsOTS3, and OsELS2. Some AS forms have the potential to alter protein domains and catalytic centers. We also performed the molecular and phenotypic characterization of T-DNA insertion lines of some of the genes under study. Knockouts of OsFUG1 and OsELS1 showed increased SUMOylation levels and non-overlapping phenotypes. The fug1 line showed a dwarf phenotype, and significant defects in fertility, seed weight, and panicle architecture, while the els1 line showed early flowering and decreased plant height. We suggest that OsELS1 is an ortholog of AtEsd4, which was also supported by our phylogenetic analysis.

CONCLUSIONS

Overall, we provide a comprehensive analysis of the rice SUMOylation machinery and discuss possible effects of the regulation of these genes at the transcriptional and post-transcriptional level. We also contribute to the characterization of two rice SUMO proteases, OsELS1 and OsFUG1.

Insertion of 275-bp SINE into first intron of PDIA4 gene is associated with litter size in Xiang pigs

The aim of the study was to investigate the SINE polymorphism in Xiang, Kele, Qianbei black, Jiangkouluobo, Large White, and Duroc pig breeds. The PCR based detection method was conducted to identify the short interspersed nuclear element (SINE) polymorphism in the PDIA4 gene. There were greater frequencies of the SINE-/- genotypes in Xiang pigs (55.9%) as compared with other pig breed groups. There was an association between this 275 bp SINE polymorphism and litter size (P = 0.003). The homozygous SINE+/+ genotype of the PDIA4 gene had a 1.45-piglets larger litter sizes compared to those with the homozygous SINE-/- genotype. Furthermore, there were assessments of mRNA and protein abundances as a result of PDIA4 gene expression in the ovaries of Xiang pigs for the three different SINE genotypes, and the results indicated that relative abundances of PDIA4 mRNA and protein was greater for the SINE-/- and SINE-/+ genotypes compared with the SINE+/+ genotype (P < 0.05). These findings suggested that the 275 bp SINE polymorphism might change the expression of the PDIA4 gene and could be used as a candidate DNA marker for the selection for litter size in Xiang pigs.

Improving Salt Tolerance of Chickpea Using Modern Genomics Tools and Molecular Breeding

INTRODUCTION

The high protein value, essential minerals, dietary fibre and notable ability to fix atmospheric nitrogen make chickpea a highly remunerative crop, particularly in low-input food production systems. Of the variety of constraints challenging chickpea productivity worldwide, salinity remains of prime concern owing to the intrinsic sensitivity of the crop. In view of the projected expansion of chickpea into arable and salt-stressed land by 2050, increasing attention is being placed on improving the salt tolerance of this crop. Considerable effort is currently underway to address salinity stress and substantial breeding progress is being made despite the seemingly highly-complex and environment-dependent nature of the tolerance trait.

CONCLUSION

This review aims to provide a holistic view of recent advances in breeding chickpea for salt tolerance. Initially, we focus on the identification of novel genetic resources for salt tolerance via extensive germplasm screening. We then expand on the use of genome-wide and cost-effective techniques to gain new insights into the genetic control of salt tolerance, including the responsive genes/QTL(s), gene(s) networks/cross talk and intricate signalling cascades.

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.

The Evolutionary Relationship between Alternative Splicing and Gene Duplication

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

Analysis of Ribosome-Associated mRNAs in Rice Reveals the Importance of Transcript Size and GC Content in Translation

Gene expression is controlled at transcriptional and post-transcriptional levels including decoding of messenger RNA (mRNA) into polypeptides via ribosome-mediated translation. Translational regulation has been intensively studied in the model dicot plant Arabidopsis thaliana, and in this study, we assessed the translational status [proportion of steady-state mRNA associated with ribosomes] of mRNAs by Translating Ribosome Affinity Purification followed by mRNA-sequencing (TRAP-seq) in rice (Oryza sativa), a model monocot plant and the most important food crop. A survey of three tissues found that most transcribed rice genes are translated whereas few transposable elements are associated with ribosomes. Genes with short and GC-rich coding regions are overrepresented in ribosome-associated mRNAs, suggesting that the GC-richness characteristic of coding sequences in grasses may be an adaptation that favors efficient translation. Transcripts with retained introns and extended 5′ untranslated regions are underrepresented on ribosomes, and rice genes belonging to different evolutionary lineages exhibited differential enrichment on the ribosomes that was associated with GC content. Genes involved in photosynthesis and stress responses are preferentially associated with ribosomes, whereas genes in epigenetic regulation pathways are the least enriched on ribosomes. Such variation is more dramatic in rice than that in Arabidopsis and is correlated with the wide variation of GC content of transcripts in rice. Taken together, variation in the translation status of individual transcripts reflects important mechanisms of gene regulation, which may have a role in evolution and diversification.

Tissue-specific expression, developmentally and spatially regulated alternative splicing, and protein subcellular localization of OsLpa rice

The OsLpa1 gene (LOC_Os57400) was identified to be involved in phytic acid (PA) metabolism because its knockout and missense mutants reduce PA content in rice grain. However, little is known about the molecular characteristics of OsLpa rice and of its homologues in other plants. In the present study, the spatial pattern of OsLpa1 expression was revealed using OsLpa1 promoter::GUS transgenic plants (GUS: β-glucuronidase); GUS histochemical assay showed that OsLpa1 was strongly expressed in stem, leaf, and root tissues, but in floral organ it is expressed mainly and strongly in filaments. In seeds, GUS staining was concentrated in the aleurone layers; a few blue spots were observed in the outer layers of embryo, but no staining was observed in the endosperm. Three OsLpa1 transcripts (OsLpa1.1, OsLpa1.2, OsLpa1.3) are produced due to alternative splicing; quantitative reverse-transcriptase polymerase chain reaction (RT-PCR) analysis revealed that the abundance of OsLpa1.3 was negligible compared with OsLpa1.1 and OsLpa all tissues. OsLpa1.2 is predominant in germinating seeds (about 5 times that of OsLpa1.1), but its abundance decreases quickly with the development of seedlings and plants, whereas the abundance of OsLpa1.1 rises and falls, reaching its highest level in 45-d-old plants, with abundance greater than that of OsLpa both leaves and roots. In seeds, the abundance of OsLpa1 continuously increases with seed growth, being 27.5 and 15 times greater in 28-DAF (day after flowering) seeds than in 7-DAF seeds for OsLpa1.1 and OsLpa1.2, respectively. Transient expression of chimeric genes with green fluorescence protein (GFP) in rice protoplasts demonstrated that all proteins encoded by the three OsLpa1 transcripts are localized to the chloroplast.

CoSREM: a graph mining algorithm for the discovery of combinatorial splicing regulatory elements

BACKGROUND

Alternative splicing (AS) is a post-transcriptional regulatory mechanism for gene expression regulation. Splicing decisions are affected by the combinatorial behavior of different splicing factors that bind to multiple binding sites in exons and introns. These binding sites are called splicing regulatory elements (SREs). Here we develop CoSREM (Combinatorial SRE Miner), a graph mining algorithm to discover combinatorial SREs in human exons. Our model does not assume a fixed length of SREs and incorporates experimental evidence as well to increase accuracy. CoSREM is able to identify sets of SREs and is not limited to SRE pairs as are current approaches.

RESULTS

We identified 37 SRE sets that include both enhancer and silencer elements. We show that our results intersect with previous results, including some that are experimental. We also show that the SRE set GGGAGG and GAGGAC identified by CoSREM may play a role in exon skipping events in several tumor samples. We applied CoSREM to RNA-Seq data for multiple tissues to identify combinatorial SREs which may be responsible for exon inclusion or exclusion across tissues.

CONCLUSION

The new algorithm can identify different combinations of splicing enhancers and silencers without assuming a predefined size or limiting the algorithm to find only pairs of SREs. Our approach opens new directions to study SREs and the roles that AS may play in diseases and tissue specificity.

Identifying splicing regulatory elements with de Bruijn graphs

Splicing regulatory elements (SREs) are short, degenerate sequences on pre-mRNA molecules that enhance or inhibit the splicing process via the binding of splicing factors, proteins that regulate the functioning of the spliceosome. Existing methods for identifying SREs in a genome are either experimental or computational. Here, we propose a formalism based on de Bruijn graphs that combines genomic structure, word count enrichment analysis, and experimental evidence to identify SREs found in exons. In our approach, SREs are not restricted to a fixed length (i.e., k-mers, for a fixed k). As a result, we identify 2001 putative exonic enhancers and 3080 putative exonic silencers for human genes, with lengths varying from 6 to 15 nucleotides. Many of the predicted SREs overlap with experimentally verified binding sites. Our model provides a novel method to predict variable length putative regulatory elements computationally for further experimental investigation.

Expression of a cyclophilin OsCyp2-P isolated from a salt-tolerant landrace of rice in tobacco alleviates stress via ion homeostasis and limiting ROS accumulation

Cyclophilins are a set of ubiquitous proteins present in all subcellular compartments, involved in a wide variety of cellular processes. Comparative bioinformatics analysis of the rice and Arabidopsis genomes led us to identify novel putative cyclophilin gene family members in both the genomes not reported previously. We grouped cyclophilin members with similar molecular weight and subtypes together in the phylogenetic tree which indicated their co-evolution in rice and Arabidopsis. We also characterized a rice cyclophilin gene, OsCyp2-P (Os02g0121300), isolated from a salinity-tolerant landrace, Pokkali. Publicly available massively parallel signature sequencing (MPSS) and microarray data, besides our quantitative real time PCR (qRT-PCR) data suggest that transcript abundance of OsCyp2-P is regulated under different stress conditions in a developmental and organ specific manner. Ectopic expression of OsCyp2-P imparted multiple abiotic stress tolerance to transgenic tobacco plants as evidenced by higher root length, shoot length, chlorophyll content, and K(+)/Na(+) ratio under stress conditions. Transgenic plants also showed reduced lipid peroxidase content, electrolyte leakage, and superoxide content under stress conditions suggesting better ion homeostasis than WT plants. Localization studies confirmed that OsCyp2-P is localized in both cytosol and nucleus, indicating its possible interaction with several other proteins. The overall results suggest the explicit role of OsCyp2-P in bestowing multiple abiotic stress tolerance at the whole plant level. OsCyp2-P operates via reactive oxygen species (ROS) scavenging and ion homeostasis and thus is a promising candidate gene for enhancing multiple abiotic stress tolerance in crop plants.

Comparative genomics of grass EST libraries reveals previously uncharacterized splicing events in crop plants

BACKGROUND

Crop plants such as rice, maize and sorghum play economically-important roles as main sources of food, fuel, and animal feed. However, current genome annotations of crop plants still suffer false-positive predictions; a more comprehensive registry of alternative splicing (AS) events is also in demand. Comparative genomics of crop plants is largely unexplored.

RESULTS

We performed a large-scale comparative analysis (ExonFinder) of the expressed sequence tag (EST) library from nine grass plants against three crop genomes (rice, maize, and sorghum) and identified 2,879 previously-unannotated exons (i.e., novel exons) in the three crops. We validated 81% of the tested exons by RT-PCR-sequencing, supporting the effectiveness of our in silico strategy. Evolutionary analysis reveals that the novel exons, comparing with their flanking annotated ones, are generally under weaker selection pressure at the protein level, but under stronger pressure at the RNA level, suggesting that most of the novel exons also represent novel alternatively spliced variants (ASVs). However, we also observed the consistency of evolutionary rates between certain novel exons and their flanking exons, which provided further evidence of their co-occurrence in the transcripts, suggesting that previously-annotated isoforms might be subject to erroneous predictions. Our validation showed that 54% of the tested genes expressed the newly-identified isoforms that contained the novel exons, rather than the previously-annotated isoforms that excluded them. The consistent results were steadily observed across cultivated (Oryza sativa and O. glaberrima) and wild (O. rufipogon and O. nivara) rice species, asserting the necessity of our curation of the crop genome annotations. Our comparative analyses also inferred the common ancestral transcriptome of grass plants and gain- and loss-of-ASV events.

CONCLUSIONS

We have reannotated the rice, maize, and sorghum genomes, and showed that evolutionary rates might serve as an indicator for determining whether the identified exons were alternatively spliced. This study not only presents an effective in silico strategy for the improvement of plant annotations, but also provides further insights into the role of AS events in the evolution and domestication of crop plants. ExonFinder and the novel exons/ASVs identified are publicly accessible at http://exonfinder.sourceforge.net/ .