Genome-wide detection of tissue-specific alternative splicing in the human transcriptome

We have developed an automated method for discovering tissue-specific regulation of alternative splicing through a genome-wide analysis of expressed sequence tags (ESTs). Using this approach, we have identified 667 tissue-specific alternative splice forms of human genes. We validated our muscle-specific and brain-specific splice forms for known genes. A high fraction (8/10) were reported to have a matching tissue specificity by independent studies in the published literature. The number of tissue-specific alternative splice forms is highest in brain, while eye-retina, muscle, skin, testis and lymph have the greatest enrichment of tissue-specific splicing. Overall, 10-30% of human alternatively spliced genes in our data show evidence of tissue-specific splice forms. Seventy-eight percent of our tissue-specific alternative splices appear to be novel discoveries. We present bioinformatics analysis of several tissue-specific splice forms, including automated protein isoform sequence and domain prediction, showing how our data can provide valuable insights into gene function in different tissues. For example, we have discovered a novel kidney-specific alternative splice form of the WNK1 gene, which appears to specifically disrupt its N-terminal kinase domain and may play a role in PHAII hypertension. Our database greatly expands knowledge of tissue-specific alternative splicing and provides a comprehensive dataset for investigating its functional roles and regulation in different human tissues.     

Conserved Alternative Splicing Patterns and Splicing Signals in the Drosophila Sodium Channel Gene Para

We cloned genomic DNA corresponding to the Drosophila virilis homologue of para, a gene encoding a sodium channel α-subunit, and obtained many partial cDNA clones from embryos and adults. Para protein has been well conserved, and the optional elements at six different sites of alternative splicing in D. melanogaster are present in D. virilis, in addition to one new optional exon. Among 31 different splice-types observed in D. virilis, the stage-specific pattern of alternative splicing seen in D. melanogaster is also conserved. Comparison of genomic DNA sequence revealed three aspects that vary between alternatively and constitutively used exon sequences. Sixteen short blocks (10-75 bp), the only recognizably conserved intron sequence, were disproportionately associated with alternatively used splice sites. Silent site substitutions were found much less frequently in alternative than constitutive exon elements, and the degree of match to the Drosophila splice site consensus tended to be lower at less frequently selected alternative splice junctions. This study shows that the developmentally regulated variability of para products is highly conserved and therefore likely to be of functional significance and suggests that a variety of different sequence-dependent mechanisms may regulate this pattern of alternative splicing.     

Splicing Reporter Mice Revealed the Evolutionally Conserved Switching Mechanism of Tissue-Specific Alternative Exon Selection

Since alternative splicing of pre-mRNAs is essential for generating tissue-specific diversity in proteome, elucidating its regulatory mechanism is indispensable to understand developmental process or tissue-specific functions. We have been focusing on tissue-specific regulation of mutually exclusive selection of alternative exons because this implies the typical molecular mechanism of alternative splicing regulation and also can be good examples to elicit general rule of “splice code”. So far, mutually exclusive splicing regulation has been explained by the outcome from the balance of multiple regulators that enhance or repress either of alternative exons discretely. However, this “balance” model is open to questions of how to ensure the selection of only one appropriate exon out of several candidates and how to switch them. To answer these questions, we generated an original bichromatic fluorescent splicing reporter system for mammals using fibroblast growth factor-receptor 2 (FGFR2) gene as model. By using this splicing reporter, we demonstrated that FGFR2 gene is regulated by the “switch-like” mechanism, in which key regulators modify the ordered splice-site recognition of two mutually exclusive exons, eventually ensure single exon selection and their distinct switching. Also this finding elucidated the evolutionally conserved “splice code,” in which combination of tissue-specific and broadly expressed RNA binding proteins regulate alternative splicing of specific gene in a tissue-specific manner. These findings provide the significant cue to understand how a number of spliced genes are regulated in various tissue-specific manners by a limited number of regulators, eventually to understand developmental process or tissue-specific functions.     

ProSplicer: a database of putative alternative splicing information derived from protein,mRNA and expressed sequence tag sequence data

ProSplicer is a database of putative alternative splicing information derived from the alignment of proteins, mRNA sequences and expressed sequence tags (ESTs) against human genomic DNA sequences. Proteins, mRNA and ESTs provide valuable evidence that can reveal splice variants of genes. The alternative splicing information in the database can help users investigate the alternative splicing and tissue-specific expression of genes.     

Genomic splice-site analysis reveals frequent alternative splicing close to the dominant splice site

Alternative pre-mRNA splicing may be the most efficient and widespread mechanism to generate multiple protein isoforms from single genes. Here, we describe the genomic analysis of one of the most frequent types of alternative pre-mRNA splicing, alternative 5'- and 3'-splice-site selection. Using an EST-based alternative splicing database recording >47,000 alternative splicing events, we determined the frequency and location of alternative 5'- and 3'-splice sites within the human genome. The most common alternative splice sites used in the human genome are located within 6 nucleotides (nt) of the dominant splice site. We show that the EST database overrepresents alternative splicing events that maintain the reading frame, thus supporting the concept that RNA quality-control steps ensure that mRNAs that encode for potentially harmful protein products are destroyed and do not serve as templates for translation. The most frequent location for alternative 5'-splice sites is 4 nt upstream or downstream from the dominant splice site. Sequence analysis suggests that this preference is a consequence of the U1 snRNP binding sequence at the 5'-splice site, which frequently contains a GU dinucleotide 4 nt downstream from the dominant splice site. Surprisingly, approximately 50% of duplicated 3'-YAG splice junctions are subject to alternative splicing. This high probability of alternative 3'-splice-site activation in close proximity of the dominant 3'-splice site suggests that the second step of the splicing may be prone to violate splicing fidelity.     

Human and chick alpha 2(I) collagen mRNA: comparison of the 5' end in osteoblasts and fibroblasts

Type I collagen, a heterotrimer of two alpha 1(I) chains and one alpha 2(I) chain, is the major structural protein of bone, skin, and tendon. The collagen of patients with bone diseases has been studied in skin fibroblasts instead of osteoblasts because the genes for type I collagen are single-copy genes. While these studies should detect structural changes in the gene, they might fail to detect defects in processes which are dependent on tissue-specific expression. The studies reported here sought to determine whether the expression of type I collagen in skin and bone was characterized by the use of alternate promoters or alternative splicing in the N-propeptide region. Primer extension and nuclease S1 protection experiments were used to analyze the structure of the alpha 2(I) mRNA from the 5' end of the gene through the N-telopeptide coding region (exons 1-6) in human and chick osteoblasts and fibroblasts. The protection and primer extension experiments using human and chick mRNA demonstrated identical routes of splicing in skin and bone at the first five splice junctions. These studies provide reassurance that information obtained from the study of type I collagen in fibroblasts may be extrapolated to bone.     

Divergent Functions Through Alternative Splicing: The Drosophila CRMP Gene in Pyrimidine Metabolism, Brain, and Behavior

DHP and CRMP proteins comprise a family of structurally similar proteins that perform divergent functions, DHP in pyrimidine catabolism in most organisms and CRMP in neuronal dynamics in animals. In vertebrates, one DHP and five CRMP proteins are products of six genes; however, Drosophila melanogaster has a single CRMP gene that encodes one DHP and one CRMP protein through tissue-specific, alternative splicing of a pair of paralogous exons. The proteins derived from the fly gene are identical over 90% of their lengths, suggesting that unique, novel functions of these proteins derive from the segment corresponding to the paralogous exons. Functional homologies of the Drosophila and mammalian CRMP proteins are revealed by several types of evidence. Loss-of-function CRMP mutation modifies both Ras and Rac misexpression phenotypes during fly eye development in a manner that is consistent with the roles of CRMP in Ras and Rac signaling pathways in mammalian neurons. In both mice and flies, CRMP mutation impairs learning and memory. CRMP mutant flies are defective in circadian activity rhythm. Thus, DHP and CRMP proteins are derived by different processes in flies (tissue-specific, alternative splicing of paralogous exons of a single gene) and vertebrates (tissue-specific expression of different genes), indicating that diverse genetic mechanisms have mediated the evolution of this protein family in animals.     

Noisy splicing drives mRNA isoform diversity in human cells

While the majority of multiexonic human genes show some evidence of alternative splicing, it is unclear what fraction of observed splice forms is functionally relevant. In this study, we examine the extent of alternative splicing in human cells using deep RNA sequencing and de novo identification of splice junctions. We demonstrate the existence of a large class of low abundance isoforms, encompassing approximately 150,000 previously unannotated splice junctions in our data. Newly-identified splice sites show little evidence of evolutionary conservation, suggesting that the majority are due to erroneous splice site choice. We show that sequence motifs involved in the recognition of exons are enriched in the vicinity of unconserved splice sites. We estimate that the average intron has a splicing error rate of approximately 0.7% and show that introns in highly expressed genes are spliced more accurately, likely due to their shorter length. These results implicate noisy splicing as an important property of genome evolution.     

ASAP: the Alternative Splicing Annotation Project

Recently, genomics analyses have demonstrated that alternative splicing is widespread in mammalian genomes (30-60% of genes reported to have multiple isoforms), and may be one of their most important mechanisms of functional regulation. However, by comparison with other genomics data such as genome annotation, SNPs, or gene expression, there exists relatively little database infrastructure for the study of alternative splicing. We have constructed an online database ASAP (the Alternative Splicing Annotation Project) for biologists to access and mine the enormous wealth of alternative splicing information coming from genomics and proteomics. ASAP is based on genome-wide analyses of alternative splicing in human (30 793 alternative splice relationships found) from detailed alignment of expressed sequences onto the genomic sequence. ASAP provides precise gene exon-intron structure, alternative splicing, tissue specificity of alternative splice forms, and protein isoform sequences resulting from alternative splicing. Moreover, it can help biologists design probe sequences for distinguishing specific mRNA isoforms. ASAP is intended to be a community resource for collaborative annotation of alternative splice forms, their regulation, and biological functions. The URL for ASAP is http://www.bioinformatics.ucla.edu/ASAP.     

Splice junctions are constrained by protein disorder

We have discovered that positions of splice junctions in genes are constrained by the tolerance for disorder-promoting amino acids in the translated protein region. It is known that efficient splicing requires nucleotide bias at the splice junction; the preferred usage produces a distribution of amino acids that is disorder-promoting. We observe that efficiency of splicing, as seen in the amino-acid distribution, is not compromised to accommodate globular structure. Thus we infer that it is the positions of splice junctions in the gene that must be under constraint by the local protein environment. Examining exonic splicing enhancers found near the splice junction in the gene, reveals that these (short DNA motifs) are more prevalent in exons that encode disordered protein regions than exons encoding structured regions. Thus we also conclude that local protein features constrain efficient splicing more in structure than in disorder.