Alternative splice variants encoding unstable protein domains exist in the human brain

Alternative splicing has been recognized as a major mechanism by which protein diversity is increased without significantly increasing genome size in animals and has crucial medical implications, as many alternative splice variants are known to cause diseases. Despite the importance of knowing what structural changes alternative splicing introduces to the encoded proteins for the consideration of its significance, the problem has not been adequately explored. Therefore, we systematically examined the structures of the proteins encoded by the alternative splice variants in the HUGE protein database derived from long (>4 kb) human brain cDNAs. Limiting our analyses to reliable alternative splice junctions, we found alternative splice junctions to have a slight tendency to avoid the interior of SCOP domains and a strong statistically significant tendency to coincide with SCOP domain boundaries. These findings reflect the occurrence of some alternative splicing events that utilize protein structural units as a cassette. However, 50 cases were identified in which SCOP domains are disrupted in the middle by alternative splicing. In six of the cases, insertions are introduced at the molecular surface, presumably affecting protein functions, while in 11 of the cases alternatively spliced variants were found to encode pairs of stable and unstable proteins. The mRNAs encoding such unstable proteins are much less abundant than those encoding stable proteins and tend not to have corresponding mRNAs in non-primate species. We propose that most unstable proteins encoded by alternative splice variants lack normal functions and are an evolutionary dead-end.     

Structural implication of splicing stochastics

Even though nearly every human gene has at least one alternative splice form, very little is so far known about the structure and function of resulting protein products. It is becoming increasingly clear that a significant fraction of all isoforms are products of noisy selection of splice sites and thus contribute little to actual functional diversity, and may potentially be deleterious. In this study, we examine the impact of alternative splicing on protein sequence and structure in three datasets: alternative splicing events conserved across multiple species, alternative splicing events in genes that are strongly linked to disease and all observed alternative splicing events. We find that the vast majority of all alternative isoforms result in unstable protein conformations. In contrast to that, the small subset of isoforms conserved across species tends to maintain protein structural integrity to a greater extent. Alternative splicing in disease-associated genes produces unstable structures just as frequently as all other genes, indicating that selection to reduce the effects of alternative splicing on this set is not especially pronounced. Overall, the properties of alternative spliced proteins are consistent with the outcome of noisy selection of splice sites by splicing machinery.     

Verification of alternative splicing variants based on domain integrity, truncation length and intrinsic protein disorder

According to current estimations ～95% of multi-exonic human protein-coding genes undergo alternative splicing (AS). However, for 4000 human proteins in PDB, only 14 human proteins have structures of at least two alternative isoforms. Surveying these structural isoforms revealed that the maximum insertion accommodated by an isoform of a fully ordered protein domain was 5 amino acids, other instances of domain changes involved intrinsic structural disorder. After collecting 505 minor isoforms of human proteins with evidence for their existence we analyzed their length, protein disorder and exposed hydrophobic surface. We found that strict rules govern the selection of alternative splice variants aimed to preserve the integrity of globular domains: alternative splice sites (i) tend to avoid globular domains or (ii) affect them only marginally or (iii) tend to coincide with a location where the exposed hydrophobic surface is minimal or (iv) the protein is disordered. We also observed an inverse correlation between the domain fraction lost and the full length of the minor isoform containing the domain, possibly indicating a buffering effect for the isoform protein counteracting the domain truncation effect. These observations provide the basis for a prediction method (currently under development) to predict the viability of splice variants.     

Altered Motor Activity of Alternative Splice Variants of the Mammalian Kinesin-3 Protein KIF1B

Several mammalian kinesin motor proteins exist as multiple isoforms that arise from alternative splicing of a single gene. However, the roles of many motor protein splice variants remain unclear. The kinesin-3 motor protein KIF1B has alternatively spliced isoforms distinguished by the presence or absence of insertion sequences in the conserved amino-terminal region of the protein. The insertions are located in the loop region containing the lysine-rich cluster, also known as the K-loop, and in the hinge region adjacent to the motor domain. To clarify the functions of these alternative splice variants of KIF1B, we examined the biochemical properties of recombinant KIF1B with and without insertion sequences. In a microtubule-dependent ATPase assay, KIF1B variants that contained both insertions had higher activity and affinity for microtubules than KIF1B variants that contained no insertions. Mutational analysis of the K-loop insertion revealed that variants with a longer insertion sequence at this site had higher activity. However, the velocity of movement in motility assays was similar between KIF1B with and without insertion sequences. Our results indicate that splicing isoforms of KIF1B that vary in their insertion sequences have different motor activities.     

Cloning of a differentially expressed tropomyosin isoform from cultured rabbit aortic smooth muscle cells

The four known tropomyosin genes have highly conserved DNA and amino acid sequences, and at least 18 isoforms are generated by alternative RNA splicing in muscle and non-muscle cells. No rabbit tropomyosin nucleotide sequences are known, although protein sequences for alpha- and beta-tropomyosin expressed by rabbit skeletal muscle have been described. Subtractive hybridisation was used to select for genes differentially expressed in rabbit aortic smooth muscle cells (SMC), during the change in cell phenotype in primary culture that is characterised by a loss of cytoskeletal filaments and contractile proteins. This led to the cloning of a tropomyosin gene predominantly expressed in rabbit SMC during this change. The full-length cDNA clone, designated "rabbit TM-beta", contains an open reading frame of 284 amino acids, 5' untranslated region (UTR) of 117 base pairs and 3' UTR of 79 base pairs. It is closely related to the beta-gene isoforms in other species, with the highest homology in DNA and protein sequences to the human fibroblast isoform TM-1 (91.7% identity in 1035 bp and 93.3% identity in the entire 284 amino acid sequence of the protein). It differs from rabbit skeletal muscle beta-tropomyosin (81.7% homology at the protein level) mainly in two regions at amino acids 189-213 and 258-283 suggesting alternative splicing of exons 6a for 6b and 9d for 9a. Since this TM-beta gene was the only gene strongly enough expressed in SMC changing phenotype to be observed by the subtractive hybridisation screen, it likely plays a significant role in this process.     

Alternative splicing and protein structure evolution

Alternative splicing is thought to be one of the major sources for functional diversity in higher eukaryotes. Interestingly, when mapping splicing events onto protein structures, about half of the events affect structured and even highly conserved regions i.e. are non-trivial on the structure level. This has led to the controversial hypothesis that such splice variants result in nonsense-mediated mRNA decay or non-functional, unstructured proteins, which do not contribute to the functional diversity of an organism. Here we show in a comprehensive study on alternative splicing that proteins appear to be much more tolerant to structural deletions, insertions and replacements than previously thought. We find literature evidence that such non-trivial splicing isoforms exhibit different functional properties compared to their native counterparts and allow for interesting regulatory patterns on the protein network level. We provide examples that splicing events may represent transitions between different folds in the protein sequence–structure space and explain these links by a common genetic mechanism. Taken together, those findings hint to a more prominent role of splicing in protein structure evolution and to a different view of phenotypic plasticity of protein structures.     

Identification and characterization of novel splice variants of the human EPM2A gene mutated in Lafora progressive myoclonus epilepsy

The EPM2A gene, defective in the fatal neurodegenerative disorder Lafora disease (LD), is known to encode two distinct proteins by differential splicing; a phosphatase active cytoplasmic isoform and a phosphatase inactive nuclear isoform. We report here the identification of three novel EPM2A splice variants with potential to code for five distinct proteins in alternate reading frames. These novel isoforms, when ectopically expressed in cell lines, show distinct subcellular localization, interact with and serve as substrates of malin ubiquitin ligase-the second protein defective in LD. Two phosphatase active isoforms interact to form a heterodimeric complex that is inactive as a phosphatase in vitro, suggesting an antagonistic function for laforin isoforms if expressed endogenously in significant amounts in human tissues. Thus alternative splicing could possibly be one of the mechanisms by which EPM2A may regulate the cellular functions of the proteins it codes for.     

Structure and expression of two kininogen genes in mice

Kininogens serve dual functions by forming a scaffold for the assembly of the protein complex initiating the surface-activated blood coagulation cascade and as precursors for the kinin hormones. While rats have three kininogen genes, for mice, cattle, and humans only one gene has been described. Here, we present sequence and expression data of a second mouse kininogen gene. The two genes, kininogen-I and kininogen-II, are located in close proximity on chromosome 16 in a head-to-head arrangement. In liver and kidney, both genes are expressed and for each gene three alternative splice variants are synthesized. Two of them are the expected high and low molecular weight isoforms known from all mammalian kininogens. However, for both genes also a third, hitherto unknown splice variant was detected which lacks part of the high molecular weight mRNA due to splicing from the low molecular weight donor site to alternative splice acceptor sites in exon 10. The physiological functions of the six kininogen isoforms predicted by these findings need to be determined.     

水稻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剪接变构体的蛋白质定位是质膜和细胞质。这些选择性剪接蛋白可能在抗病信号转导中起到重要作用。     

Theoretical analysis of alternative splice forms using computational methods

Nowadays understanding alternative splicing is one of the greatest challenges in biology, because it is a genetic process much more important than thought at the time of its discovery. In this paper, we explain the approach of using the different available databases and software tools to start a large scale investigation of alternative splice forms. To collect information about alternative splicing we investigated known data in the databases using different computational methods. The investigations proceeded from the genomic sequence data to structural protein data. Then, we interpreted those data to find the relationship between alternative splice forms and protein function and structure. We found some interesting features of alternative splicing which are presented here. We discuss the results of one chosen example. They concern the coverage quality of the protein sequence of a known structure, an EST analysis, the validation of splice variants, the determination of the alternative splice type, and finally the link between alternative splicing and disease.     

Tumor microenvironment–associated modifications of alternative splicing

Pre-mRNA alternative splicing is modified in cancer, but the origin and specificity of these changes remain unclear. Here, we probed ovarian tumors to identify cancer-associated splicing isoforms and define the mechanism by which splicing is modified in cancer cells. Using high-throughput quantitative PCR, we monitored the expression of splice variants in laser-dissected tissues from ovarian tumors. Surprisingly, changes in alternative splicing were not limited to the tumor tissues but were also found in the tumor microenvironment. Changes in the tumor-associated splicing events were found to be regulated by splicing factors that are differentially expressed in cancer tissues. Overall, ∼20% of the alternative splicing events affected by the down-regulation of the splicing factors QKI and RBFOX2 were altered in the microenvironment of ovarian tumors. Together, our results indicate that the tumor microenvironment undergoes specific changes in alternative splicing orchestrated by a limited number of splicing factors.     

Identification of potential proteins translated from circular RNA splice variants

Circular RNAs (circRNAs) are covalently closed RNA molecules generated from precursor RNAs by the head-to-tail backsplicing of exons. Hundreds of studies demonstrated that circRNAs are ubiquitously expressed and regulate cellular events by modulating microRNA (miRNA) and RNA-binding protein (RBP) activities. A few circRNAs are also known to translate into functional polypeptides regulating cellular physiology. All these functions primarily depend on the full-length sequence of the circRNAs. CircRNA backsplice junction sequence is the key to identifying circRNAs and their full-length mature sequence. However, some multi-exonic circRNAs exist in different isoforms sharing identical backsplice junction sequences and are termed circRNA splice variants. Here, we analyzed the previously published HeLa cell RNA-seq datasets to identify circRNA splice variants using the de novo module of the CIRCexplorer2 circRNA annotation pipeline. A subset of circRNAs with splice variants was validated by the circRNA-rolling circle amplification (circRNA-RCA) method. Interestingly, several validated circRNAs were predicted to translate into proteins by the riboCIRC database. Furthermore, polyribosome fractionation followed by quantitative PCR confirmed the association of a subset of circRNAs with polyribosome supporting their protein-coding potential. Finally, bioinformatics analysis of proteins derived from splice variants of circCORO1C and circASPH suggested altered protein sequences and structures that could affect their physiological functions. Together, our study identified novel circRNA splice variants and their potential translation into protein isoforms which may regulate various physiological processes.