Extent and diversity of human alternative splicing established by complementary database annotation and microarray analysis

Alternative splicing generates functional diversity in higher organisms through alternative first and last exons, skipped and included exons, intron retentions and alternative donor, and acceptor sites. In large-scale microarray studies in humans and the mouse, emphasis so far has been placed on exon-skip events, leaving the prevalence and importance of other splice types largely unexplored. Using a new human splice variant database and a genome-wide microarray to probes thousands of splice events of each type, we measured differential expression of splice types across six pair of diverse cell lines and validated the database annotation process. Results suggest that splicing in humans is more complex than simple exon-skip events, which account for a minority of splicing differences. The relative frequency of differential expression of the splice types correlates with what is found by our annotation efforts. In conclusion, alternative splicing in human cells is considerably more complex than the canonical example of the exon skip. The complementary approaches of genome-wide annotation of alternative splicing in human and design of genome-wide splicing microarrays to measure differential splicing in biological samples provide a powerful high-throughput tool to study the role of alternative splicing in human biology.     

StressDB: a locally installable web-based relational microarray database designed for small user communities

We have built a microarray database, StressDB, for management of microarray data from our studies on stress-modulated genes in Arabidopsis. StressDB provides small user groups with a locally installable web-based relational microarray database. It has a simple and intuitive architecture and has been designed for cDNA microarray technology users. StressDB uses Windows(trade mark) 2000 as the centralized database server with Oracle(trade mark) 8i as the relational database management system. It allows users to manage microarray data and data-related biological information over the Internet using a web browser. The source-code is currently available on request from the authors and will soon be made freely available for downloading from our website athttp://arastressdb.cac.psu.edu.     

AVATAR: a database for genome-wide alternative splicing event detection using large scale ESTs and mRNAs

In the past years, identification of alternative splicing (AS) variants has been gaining momentum. We developed AVATAR, a database for documenting AS using 5,469,433 human EST sequences and 26,159 human mRNA sequences. AVATAR contains 12000 alternative splicing sites identified by mapping ESTs and mRNAs with the whole human genome sequence. AVATAR also contains AS information for 6 eukaryotes. We mapped EST alignment information into a graph model where exons and introns are represented with vertices and edges, respectively. AVATAR can be queried using, (1) gene names, (2) number of identified AS events in a gene, (3) minimal number of ESTs supporting a splicing site, etc. as search parameters. The system provides visualized AS information for queried genes.

Availability     

Blood coagulation and alternative pre-mRNA splicing: an overview

Throughout the 20th century, great advances were made in understanding of how blood coagulation occurs, what physiological and biochemical mechanisms are responsible for its regulation, and what genes and their protein products comprise the essential components of the hemostatic network. Recently, complete sequencing of the human genome revealed that the structural diversity of higher eukaryotes cannot be solely attributed to the number of protein-encoding genes, whereas tools of molecular biology helped establish that pre-mRNAs produced by most protein-encoding genes undergo alternative splicing, a mechanism that enables production of multiple protein isoforms by a single gene. Research in the field of thrombosis and hemostasis revealed that the genes encoding several critical proteins at various junctures of the coagulation cascade produce alternatively spliced protein isoforms with distinct structural and biochemical characteristics, revealing a principally novel dimension in the regulation of blood clotting and, possibly, a few novel therapeutic approaches to treatment of abnormal hemostasis. This review summarizes recently published data pertaining to biosynthesis of the alternatively spliced isoforms of tissue factor (TF, or coagulation factor III), tissue factor pathway inhibitor (TFPI), and coagulation factor XI (FXI), and discusses future directions of this continuously evolving area of biomedical research, with an emphasis on molecular mechanics responsible for regulation of constitutive as well as alternative pre-mRNA splicing.     

HMM sampling and applications to gene finding and alternative splicing

The standard method of applying hidden Markov models to biological problems is to find a Viterbi (maximal weight) path through the HMM graph. The Viterbi algorithm reduces the problem of finding the most likely hidden state sequence that explains given observations, to a dynamic programming problem for corresponding directed acyclic graphs. For example, in the gene finding application, the HMM is used to find the most likely underlying gene structure given a DNA sequence. In this note we discuss the applications of sampling methods for HMMs. The standard sampling algorithm for HMMs is a variant of the common forward-backward and backtrack algorithms, and has already been applied in the context of Gibbs sampling methods. Nevetheless, the practice of sampling state paths from HMMs does not seem to have been widely adopted, and important applications have been overlooked. We show how sampling can be used for finding alternative splicings for genes, including alternative splicings that are conserved between genes from related organisms. We also show how sampling from the posterior distribution is a natural way to compute probabilities for predicted exons and gene structures being correct under the assumed model. Finally, we describe a new memory efficient sampling algorithm for certain classes of HMMs which provides a practical sampling alternative to the Hirschberg algorithm for optimal alignment. The ideas presented have applications not only to gene finding and HMMs but more generally to stochastic context free grammars and RNA structure prediction.     

Inferring global levels of alternative splicing isoforms using a generative model of microarray data

MOTIVATION: Alternative splicing (AS) is a frequent step in metozoan gene expression whereby the exons of genes are spliced in different combinations to generate multiple isoforms of mature mRNA. AS functions to enrich an organism's proteomic complexity and regulates gene expression. Despite its importance, the mechanisms underlying AS and its regulation are not well understood, especially in the context of global gene expression patterns. We present here an algorithm referred to as the Generative model for the Alternative Splicing Array Platform (GenASAP) that can predict the levels of AS for thousands of exon skipping events using data generated from custom microarrays. GenASAP uses Bayesian learning in an unsupervised probability model to accurately predict AS levels from the microarray data. GenASAP is capable of learning the hybridization profiles of microarray data, while modeling noise processes and missing or aberrant data. GenASAP has been successfully applied to the global discovery and analysis of AS in mammalian cells and tissues. RESULTS: GenASAP was applied to data obtained from a custom microarray designed for the monitoring of 3126 AS events in mouse cells and tissues. The microarray design included probes specific for exon body and junction sequences formed by the splicing of exons. Our results show that GenASAP provides accurate predictions for over one-third of the total events, as verified by independent RT-PCR assays. SUPPLEMENTARY INFORMATION: http://www.psi.toronto.edu/GenASAP.     

An alternative protein splicing mechanism for inteins lacking an N-terminal nucleophile

Variations in the intein-mediated protein splicing mechanism are becoming more apparent as polymorphisms in conserved catalytic residues are identified. The conserved Ser or Cys at the intein N-terminus and the conserved intein penultimate His are absent in the KlbA family of inteins. These inteins were predicted to be inactive, since an N-terminal Ala cannot perform the initial reaction of the standard protein splicing pathway to yield the requisite N-terminal splice junction (thio)ester. Despite the presence of an N-terminal Ala and a penultimate Ser, the KlbA inteins splice efficiently using an alternative protein splicing mechanism. In this non-canonical pathway, the C-extein nucleophile attacks a peptide bond at the N-terminal splice junction rather than a (thio)ester bond, alleviating the need to form the initial (thio)ester at the N-terminal splice junction. The remainder of the two pathways is the same: branch resolution by Asn cyclization is followed by an acyl rearrangement to form a native peptide bond between the ligated exteins.