Computational comparative analyses of alternative splicing regulation using full-length cDNA of various eukaryotes

We previously reported a computational approach to infer alternative splicing patterns from Mus musculus full-length cDNA clones and microarray data. Although we predicted a large number of unreported splice variants, the general mechanisms regulating alternative splicing were yet unknown. In the present study, we compared alternative exons and constitutive exons in terms of splice-site strength and frequency of potential regulatory sequences. These regulatory features were further compared among five different species: Homo sapiens, M. musculus, Arabidopsis thaliana, Oryza sativa, and Drosophila melanogaster. Solid statistical validations of our comparative analyses indicated that alternative exons have (1) weaker splice sites and (2) more potential regulatory sequences than constitutive exons. Based on our observations, we propose a combinatorial model of alternative splicing mechanisms, which suggests that alternative exons contain weak splice sites regulated alternatively by potential regulatory sequences on the exons.     

A developmental model for generating frequency maps in the reptilian and avian cochleas.

Hair cells in the turtle cochlea are frequency-tuned by a mechanism involving the combined activation of voltage-sensitive Ca2+ channels and Ca(2+)-activated K+ (KCa) channels. The main determinants of a hair cell's characteristic frequency (Fo) are the KCa channels' density and kinetics, both of which change systematically with location in the cochlea in conjunction with the observed frequency map. We have developed a model based on the differential expression of two KCa channel subunits, which when accompanied by concurrent changes in other properties (e.g., density of Ca2+ channels and inwardly rectifying K+ channels), will generate sharp tuning at frequencies from 40 to 600 Hz. The kinetic properties of the two subunits were derived from previous single-channel analysis, and it was assumed that the subunits (A and B) combine to form five species of tetrameric channel (A4, A3B, A2B2, AB3, and B4) with intermediate kinetics and overlapping distribution. Expression of KCa and other channels was assumed to be regulated by diffusional gradients in either one or two chemicals. The results are consistent with both current- and voltage-clamp data on turtle hair cells, and they show that five channel species are sufficient to produce smooth changes in both Fo and kinetics of the macroscopic KCa current. Other schemes for varying KCa channel kinetics are examined, including one that allows extension of the model to the chick cochlea to produce hair cells with Fo's from 130 to 4000 Hz. A necessary assumption in all models is a gradient in the values of the parameters identified with the cell's cytoplasmic Ca2+ buffer.     

Identification of common genetic variation that modulates alternative splicing

Alternative splicing of genes is an efficient means of generating variation in protein function. Several disease states have been associated with rare genetic variants that affect splicing patterns. Conversely, splicing efficiency of some genes is known to vary between individuals without apparent ill effects. What is not clear is whether commonly observed phenotypic variation in splicing patterns, and hence potential variation in protein function, is to a significant extent determined by naturally occurring DNA sequence variation and in particular by single nucleotide polymorphisms (SNPs). In this study, we surveyed the splicing patterns of 250 exons in 22 individuals who had been previously genotyped by the International HapMap Project. We identified 70 simple cassette exon alternative splicing events in our experimental system; for six of these, we detected consistent differences in splicing pattern between individuals, with a highly significant association between splice phenotype and neighbouring SNPs. Remarkably, for five out of six of these events, the strongest correlation was found with the SNP closest to the intron–exon boundary, although the distance between these SNPs and the intron–exon boundary ranged from 2 bp to greater than 1,000 bp. Two of these SNPs were further investigated using a minigene splicing system, and in each case the SNPs were found to exert cis-acting effects on exon splicing efficiency in vitro. The functional consequences of these SNPs could not be predicted using bioinformatic algorithms. Our findings suggest that phenotypic variation in splicing patterns is determined by the presence of SNPs within flanking introns or exons. Effects on splicing may represent an important mechanism by which SNPs influence gene function.     

The human polycystic kidney disease 2-like (PKDL) gene: exon/intron structure and evidence for a novel splicing mechanism

Polycystin-L is a member of the expanding family of polycystins. Mutations in polycystin-1 or -2 cause human autosomal dominant polycystic kidney disease (ADPKD). The mouse ortholog of PKDL, Pkdl, is deleted in a mouse line with renal and retinal defects. We recently have shown that polycystin-L has calcium channel properties. In the current study, we determined the exon/intron organization of the PKDL gene and its alternative splicing. We show that PKDL has 16 exons. All splice acceptor/donor sites for these exons conform to the GT-AG rule. The positions of introns and the sizes of exons in the PKDL gene are very similar to those of PKD2, except for the last two 3′ end exons. RT-PCR demonstrates the existence of at least three polycystin-L splice variants: PKDL(Δ5), PKDL(Δ456), and PKDL(Δ15) that are expressed in a tissue-specific manner. In addition, we have localized polymorphic marker D10S603 to intron 4 and exon 5 of PKDL. Elucidation of the gene structure, exact location, and alternative splicing patterns of PKDL will facilitate its evaluation as a candidate gene in cystic or other genetic disorders. Received: 26 July 1999 / Accepted: 16 September 1999     

Evidence that RNA editing modulates splice site selection in the 5-HT2C receptor gene

Adenosine to inosine editing of mRNA from the human 5-HT2C receptor gene (HTR2C) occurs at five exonic positions (A–E) in a stable stem–loop that includes the normal 5′ splice site of intron 5 and is flanked by two alternative splice sites. Using in vitro editing, we identified a novel editing site (F) located in the intronic part of the stem–loop and demonstrated editing at this site in human brain. We have shown that in cell culture, base substitutions to mimic editing at different combinations of the six sites profoundly affect relative splicing at the normal and the upstream alternative splice site, but splicing at the downstream alternative splice site was consistently rare. Editing combinations in different splice variants from human brain were determined and are consistent with the effects of editing on splicing observed in cell culture. As RNA editing usually occurs close to exon/intron boundaries, this is likely to be a general phenomenon and suggests an important novel role for RNA editing.     

A splice variant of Arabidopsis mitogen-activated protein kinase and its regulatory function in the MKK6–MPK13 pathway

Alternative splicing (AS) generates multiple forms of proteins. A role for AS in the plant mitogen-activated protein kinase (MAPK) cascade has not been clarified. In this study, we analyzed expression of 20 Arabidopsis MAPK genes by RT-PCR and found five that generate splice variants. The MPK13 gene, with six exons and five introns, generates at least three splice variants, one in which complete splicing of five introns occurs (MPK13 Full), and ones in which the 4th and 5th introns are retained (MPK13 I4 and I5). Translation products of the splice variants MPK13 Full and I4, were found in Arabidopsis tissues by Western blot. Yeast two-hybrid analysis and protein kinase assays of recombinant proteins showed that neither I4 nor I5 interacted with upstream MAPKKs, and neither had protein kinase activity. However, MPK13 I4 protein enhanced the MKK6-dependent activation of MPK13 Full, indicating the possibility of an additional mechanism to regulate the MAPK cascade by AS.     

Tissue-specific alternative splicing of rat brain fructose 6-phosphate 2-kinase/fructose 2,6-bisphosphatase

We have reported the occurrence of eight splice variants of rat brain fructose 6-phosphate 2-kinase/fructose 2,6-bisphosphatase (RB2K). In the present study, we quantified these splice variants in various tissues using a RNAse protection assay and found a tissue-specific pattern of alternative splicing of the RB2K gene. Splice variants containing exon F were specifically expressed in brain. Moreover, exons D and E were spliced in brain, skeletal muscle and heart. Consequently, eight, six, four and two splice variants were expressed in brain, skeletal muscle, heart and liver plus testis, respectively. These results suggest that distinct RB2K isoforms could be involved in regulation of glycolysis in a tissue-specific manner.