An algorithm to predict 3 ' intron splice sites in Plasmodium falciparum genomic sequences

A new algorithm, PfAGSS, for predicting 3' splice sites in Plasmodium falciparum genomic sequences is described. Application of this program to the published P. falciparum chromosome 2 and 3 data suggests that existing programs result in a high error rate in assigning 3' intron boundaries. (C) 2001 Elsevier Science B.V. All rights reserved.

Comparison of splice sites in mammals and chicken

We have carried out an initial analysis of the dynamics of the recent evolution of the splice-sites sequences on a large collection of human, rodent (mouse and rat), and chicken introns. Our results indicate that the sequences of splice sites are largely homogeneous within tetrapoda. We have also found that orthologous splice signals between human and rodents and within rodents are more conserved than unrelated splice sites, but the additional conservation can be explained mostly by background intron conservation. In contrast, additional conservation over background is detectable in orthologous mammalian and chicken splice sites. Our results also indicate that the U2 and U12 intron classes seem to have evolved independently since the split of mammals and birds; we have not been able to find a convincing case of interconversion between these two classes in our collections of orthologous introns. Similarly, we have not found a single case of switching between AT-AC and GT-AG subtypes within U12 introns, suggesting that this event has been a rare occurrence in recent evolutionary times. Switching between GT-AG and the noncanonical GC-AG U2 subtypes, on the contrary, does not appear to be unusual; in particular, T to C mutations appear to be relatively well tolerated in GT-AG introns with very strong donor sites.

Sensitivity of splice sites to antisense oligonucleotides in vivo.

A series of HeLa cell lines which stably express beta-globin pre-mRNAs carrying point mutations at nt 654, 705, or 745 of intron 2 has been developed. The mutations generate aberrant 5' splice sites and activate a common 3' cryptic splice site upstream leading to aberrantly spliced beta-globin mRNA. Antisense oligonucleotides, which in vivo blocked aberrant splice sites and restored correct splicing of the pre-mRNA, revealed major differences in the sensitivity of these sites to antisense probes. Although the targeted pre-mRNAs differed only by single point mutations, the effective concentrations of the oligonucleotides required for correction of splicing varied up to 750-fold. The differences among the aberrant 5' splice sites affected sensitivity of both the 5' and 3' splice sites; in particular, sensitivity of both splice sites was severely reduced by modification of the aberrant 5' splice sites to the consensus sequence. These results suggest large differences in splicing of very similar pre-mRNAs in vivo. They also indicate that antisense oligonucleotides may provide useful tools for studying the interactions of splicing machinery with pre-mRNA.

Splicing of constitutive upstream introns is essential for the recognition of intra-exonic suboptimal splice sites in the thrombopoietin gene

The human thrombopoietin (TPO) gene, which codes for the principal cytokine involved in platelet maturation, shows a peculiar alternative splicing of its last exon, where an intra-exonic 116 nt alternative intron is spliced out in a fraction of its mRNA. To characterize the molecular mechanism underlying this alternative splicing, minigenes of TPO genomic constructs with variable exon–intron configurations or carrying exclusively the TPO cDNA were generated and transiently transfected in the Hep3B cell line. We have found that the final rate of the alternative intron splicing is determined by three elements: the presence of upstream constitutive introns, the suboptimal splice sites of the alternative intron and the length of the alternative intron itself. Our results indicate that the recognition of suboptimal intra-exonic splice junctions in the TPO gene is influenced by the assembly of the spliceosome complex on constitutive introns and by a qualitative scanning of the sequence by the transcriptional/splicing machinery complex primed by upstream splicing signals.

SpliceDB: database of canonical and non-canonical mammalian splice sites

A database (SpliceDB) of known mammalian splice site sequences has been developed. We extracted 43 337 splice pairs from mammalian divisions of the gene-centered Infogene database, including sites from incomplete or alternatively spliced genes. Known EST sequences supported 22 815 of them. After discarding sequences with putative errors and ambiguous location of splice junctions the verified dataset includes 22 489 entries. Of these, 98.71% contain canonical GT–AG junctions (22 199 entries) and 0.56% have non-canonical GC–AG splice site pairs. The remainder (0.73%) occurs in a lot of small groups (with a maximum size of 0.05%). We especially studied non-canonical splice sites, which comprise 3.73% of GenBank annotated splice pairs. EST alignments allowed us to verify only the exonic part of splice sites. To check the conservative dinucleotides we compared sequences of human non-canonical splice sites with sequences from the high throughput genome sequencing project (HTG). Out of 171 human non-canonical and EST-supported splice pairs, 156 (91.23%) had a clear match in the human HTG. They can be classified after sequence analysis as: 79 GC–AG pairs (of which one was an error that corrected to GC–AG), 61 errors corrected to GT–AG canonical pairs, six AT–AC pairs (of which two were errors corrected to AT–AC), one case was produced from a non-existent intron, seven cases were found in HTG that were deposited to GenBank and finally there were only two other cases left of supported non-canonical splice pairs. The information about verified splice site sequences for canonical and non-canonical sites is presented in SpliceDB with the supporting evidence. We also built weight matrices for the major splice groups, which can be incorporated into gene prediction programs. SpliceDB is available at the computational genomic Web server of the Sanger Centre: http://genomic.sanger.ac.uk/spldb/SpliceDB.html and at http://www.softberry.com/spldb/SpliceDB.html.

Enhancer-dependent interaction between 5' and 3' splice sites in trans.

Splice-site selection and alternative splicing of nuclear pre-mRNAs can be controlled by splicing enhancers that act by promoting the activity of upstream splice sites. Here we show that RNA molecules containing a 3' splice site and enhancer sequence are efficiently spliced in trans to RNA molecules containing normally cis-spliced 5' splice sites or to normally trans-spliced spliced leader RNAs from lower eukaryotes. In addition, we show that this reaction is stimulated by (Ser + Arg)-rich splicing factors that are known to promote protein-protein interactions in the cis-splicing reaction. Thus, splicing enhancers facilitate the assembly of protein complexes on RNAs containing a 3' splice site, and this complex is sufficiently stable to functionally interact with 5' splice sites located on separate RNAs. This trans-splicing is mediated by interactions between (Ser + Arg)-rich splicing factors bound to the enhancer and general splicing factors bound to the 5' and 3' splice sites. These same interactions are likely to play a crucial role in alternative splicing and splice-site selection in cis.

The choice of alternative 5' splice sites in influenza virus M1 mRNA is regulated by the viral polymerase complex.

The influenza virus M1 mRNA has two alternative 5' splice sites: a distal 5' splice site producing mRNA3 that has the coding potential for 9 amino acids and a proximal 5' splice site producing M2 mRNA encoding the essential M2 ion-channel protein. Only mRNA3 was made in uninfected cells transfected with DNA expressing M1 mRNA. Similarly, using nuclear extracts from uninfected cells, in vitro splicing of M1 mRNA yielded only mRNA3. Only when the mRNA3 5' splice site was inactivated by mutation was M2 mRNA made in uninfected cells and in uninfected cell extracts. In influenza virus-infected cells, M2 mRNA was made, but only after a delay, suggesting that newly synthesized viral gene product(s) were needed to activate the M2 5' splice site. We present strong evidence that these gene products are the complex of the three polymerase proteins, the same complex that functions in the transcription and replication of the viral genome. Gel shift experiments showed that the viral polymerase complex bound to the 5' end of the viral M1 mRNA in a sequence-specific and cap-dependent manner. During in vitro splicing catalyzed by uninfected cell extracts, the binding of the viral polymerase complex blocked the mRNA3 5' splice site, resulting in the switch to the M2 mRNA 5' splice site and the production of M2 mRNA.

Distinct functions of SR proteins in recruitment of U1 small nuclear ribonucleoprotein to alternative 5' splice sites.

Alternative splicing of precursor messenger RNAs (pre-mRNAs) is an important mechanism for the regulation of gene expression. The members of the SR protein family of pre-mRNA splicing factors have distinct functions in promoting alternative splice site usage. Here we show that SR proteins are required for the first step of spliceosome assembly, interaction of the U1 small nuclear ribonucleoprotein complex (U1 snRNP) with the 5' splice site of the pre-mRNA. Further, we find that individual SR proteins have distinct abilities to promote interaction of U1 snRNP with alternative 5' splice junctions. These results suggest that SR proteins direct 5' splice site selection by regulation of U1 snRNP assembly onto the pre-mRNA.

Evolutionary origin and functions of retrogene introns.

Retroposed genes (retrogenes) originate via the reverse transcription of mature messenger RNAs from parental source genes and are therefore usually devoid of introns. Here, we characterize a particular set of mammalian retrogenes that acquired introns upon their emergence and thus represent rare cases of intron gain in mammals. We find that although a few retrogenes evolved introns in their coding or 3' untranslated regions (untranslated region, UTR), most introns originated together with untranslated exons in the 5' flanking regions of the retrogene insertion site. They emerged either de novo or through fusions with 5' UTR exons of host genes into which the retrogenes inserted. Generally, retrogenes with introns display high transcription levels and show broader spatial expression patterns than other retrogenes. Our experimental expression analyses of individual intron-containing retrogenes show that 5' UTR introns may indeed promote higher expression levels, at least in part through encoded regulatory elements. By contrast, 3' UTR introns may lead to downregulation of expression levels via nonsense-mediated decay mechanisms. Notably, the majority of retrogenes with introns in their 5' flanks depend on distant, sometimes bidirectional CpG dinucleotide-enriched promoters for their expression that may be recruited from other genes in the genomic vicinity. We thus propose a scenario where the acquisition of new 5' exon-intron structures was directly linked to the recruitment of distant promoters by these retrogenes, a process potentially facilitated by the presence of proto-splice sites in the genomic vicinity of retrogene insertion sites. Thus, the primary role and selective benefit of new 5' introns (and UTR exons) was probably initially to span the often substantial distances to potent CpG promoters driving retrogene transcription. Later in evolution, these introns then obtained additional regulatory roles in fine tuning retrogene expression levels. Our study provides novel insights regarding mechanisms underlying the origin of new introns, the evolutionary relevance of intron gain, and the origin of new gene promoters.

Oscheius tipulae as an example of eEF1A gene diversity in nematodes

We characterized four eEF1A genes in the alternative rhabditid nematode model organism Oscheius tipulae. This is twice the copy number of eEF1A genes in C. elegans, C. briggsae, and, probably, many other free-living and parasitic nematodes. The introns show features remarkably different from those of other metazoan eEF1A genes. Most of the introns in the eEF1A genes are specific to O. tipulae and are not shared with any of the other genes described in metazoans. Most of the introns are phase 0 (inserted between two codons), and few are inserted in protosplice sites (introns inserted between the nucleotide sequence A/CAG and G/A). Two of these phase 0 introns are conserved in sequence in two or more of the four eEF1A gene copies, and are inserted in the same position in the genes. Neither of these characteristics has been detected in any of the nematode eEF1A genes characterized to date. The coding sequences were also compared with other eEF1A cDNAs from 11 different nematodes to determine the variability of these genes within the phylum Nematoda. Parsimony and distance trees yielded similar topologies, which were similar to those created using other molecular markers. The presence of more than one copy of the eEF1A gene with nearly identical coding regions makes it difficult to define the orthologous cDNAs. As shown by our data on O. tipulae, careful and extensive examination of intron positions in the eEF1A gene across the phylum is necessary to define their potential for use as valid phylogenetic markers.

Heteroduplex formation and S1 digestion for mapping alternative splicing sites

The identification of alternatively spliced transcripts has contributed to a better comprehension of developmental mechanisms, tissue-specific physiological processes and human diseases. Polymerase chain reaction amplification of alternatively spliced variants commonly leads to the formation of heteroduplexes as a result of base pairing involving exons common between the two variants. S1 nuclease cleaves single-stranded loops of heteroduplexes and also nicks the opposite DNA strand. In order to establish a strategy for mapping alternative splice-prone sites in the whole transcriptome, we developed a method combining the formation of heteroduplexes between 2 distinct splicing variants and S1 nuclease digestion. For 20 consensuses identified here using this methodology, 5 revealed a conserved splice site after inspection of the cDNA alignment against the human genome (exact splice sites). For 8 other consensuses, conserved splice sites were mapped at 2 to 30 bp from the border, called proximal splice sites; for the other 7 consensuses, conserved splice sites were mapped at 40 to 800 bp, called distal splice sites. These latter cases showed a nonspecific activity of S1 nuclease in digesting double-strand DNA. From the 20 consensuses identified here, 5 were selected for reverse transcription-polymerase chain reaction validation, confirming the splice sites. These data showed the potential of the strategy in mapping splice sites. However, the lack of specificity of the S1 nuclease enzyme is a significant obstacle that impedes the use of this strategy in large-scale studies.

Identification of a novel splice-site mutation in the CYP1A2 gene.

AIMS: To identify the molecular basis for a low CYP1A2 metabolic status, as determined by a caffeine phenotyping test, in a 71-year-old, nonsmoking, Caucasian woman who presented with very high clozapine concentrations despite being administered a standard dose of the drug. METHODS: The nucleotide sequence of the 7 exons, exon-intron boundaries and 5'-flanking region of the CYP1A2 gene was analysed by direct sequencing. RESULTS: Only one heterozygous point mutation was identified in the donor splice site of intron 6 (3534G > A) of CYP1A2. This mutation could cause abnormal RNA splicing and therefore lead to a truncated nonfunctional enzyme. No other carrier of this mutation was identified in a population of 100 unrelated healthy Caucasians. CONCLUSIONS: This is the first report of a splice-site mutation affecting the CYP1A2 gene. This polymorphism is a likely explanation for the low CYP1A2 activity associated with high clozapine concentrations in this patient.

Évolution à fine échelle des sites d'épissage des introns dans les gènes des oomycètes

Les introns sont des portions de gènes transcrites dans l’ARN messager, mais retirées pendant l’épissage avant la synthèse des produits du gène. Chez les eucaryotes, on rencontre les introns splicéosomaux, qui sont retirés de l’ARN messager par des splicéosomes. Les introns permettent plusieurs processus importants, tels que l'épissage alternatif, la dégradation des ARNs messagers non-sens, et l'encodage d'ARNs fonctionnels. Leurs rôles nous interrogent sur l'influence de la sélection naturelle sur leur évolution. Nous nous intéressons aux mutations qui peuvent modifier les produits d'un gène en changeant les sites d'épissage des introns. Ces mutations peuvent influencer le fonctionnement d'un organisme, et constituent donc un sujet d'étude intéressant, mais il n'existe actuellement pas de logiciels permettant de les étudier convenablement. Le but de notre projet était donc de concevoir une méthode pour détecter et analyser les changements des sites d'épissage des introns splicéosomaux. Nous avons finalement développé une méthode qui repère les évènements évolutifs qui affectent les introns splicéosomaux dans un jeu d'espèces données. La méthode a été exécutée sur un ensemble d'espèces d'oomycètes. Plusieurs évènements détectés ont changé les sites d’épissage et les protéines, mais de nombreux évènements trouvés ont modifié les introns sans affecter les produits des gènes. Il manque à notre méthode une étape finale d'analyse approfondie des données récoltées. Cependant, la méthode actuelle est facilement reproductible et automatise l'analyse des génomes pour la détection des évènements. Les fichiers produits peuvent ensuite être analysés dans chaque étude pour répondre à des questions spécifiques.

Isolated autosomal dominant growth hormone deficiency: stimulating mutant GH-1 gene expression drives GH-1 splice-site selection, cell proliferation, and apoptosis

The majority of mutations that cause isolated GH deficiency type II (IGHD II) affect splicing of GH-1 transcripts and produce a dominant-negative GH isoform lacking exon 3 resulting in a 17.5-kDa isoform, which further leads to disruption of the GH secretory pathway. A clinical variability in the severity of the IGHD II phenotype depending on the GH-1 gene alteration has been reported, and in vitro and transgenic animal data suggest that the onset and severity of the phenotype relates to the proportion of 17.5-kDa produced. The removal of GH in IGHD creates a positive feedback loop driving more GH expression, which may itself increase 17.5-kDa isoform productions from alternate splice sites in the mutated GH-1 allele. In this study, we aimed to test this idea by comparing the impact of stimulated expression by glucocorticoids on the production of different GH isoforms from wild-type (wt) and mutant GH-1 genes, relying on the glucocorticoid regulatory element within intron 1 in the GH-1 gene. AtT-20 cells were transfected with wt-GH or mutated GH-1 variants (5'IVS-3 + 2-bp T->C; 5'IVS-3 + 6 bp T->C; ISEm1: IVS-3 + 28 G->A) known to cause clinical IGHD II of varying severity. Cells were stimulated with 1 and 10 mum dexamethasone (DEX) for 24 h, after which the relative amounts of GH-1 splice variants were determined by semiquantitative and quantitative (TaqMan) RT-PCR. In the absence of DEX, only around 1% wt-GH-1 transcripts were the 17.5-kDa isoform, whereas the three mutant GH-1 variants produced 29, 39, and 78% of the 17.5-kDa isoform. DEX stimulated total GH-1 gene transcription from all constructs. Notably, however, DEX increased the amount of 17.5-kDa GH isoform relative to the 22- and 20-kDa isoforms produced from the mutated GH-1 variants, but not from wt-GH-1. This DEX-induced enhancement of 17.5-kDa GH isoform production, up to 100% in the most severe case, was completely blocked by the addition of RU486. In other studies, we measured cell proliferation rates, annexin V staining, and DNA fragmentation in cells transfected with the same GH-1 constructs. The results showed that that the 5'IVS-3 + 2-bp GH-1 gene mutation had a more severe impact on those measures than the splice site mutations within 5'IVS-3 + 6 bp or ISE +28, in line with the clinical severity observed with these mutations. Our findings that the proportion of 17.5-kDa produced from mutant GH-1 alleles increases with increased drive for gene expression may help to explain the variable onset progression, and severity observed in IGHD II.

Role of adenine functional groups in the recognition of the 3′-splice-site AG during the second step of pre-mRNA splicing

The AG dinucleotide at the 3′ splice sites of metazoan nuclear pre-mRNAs plays a critical role in catalytic step II of the splicing reaction. Previous studies have shown that replacement of the guanine by adenine in the AG (AG → GG) inhibits this step. We find that the second step was even more severely inhibited by cytosine (AG → CG) or uracil (AG → UG) substitutions at this position. By contrast, a relatively moderate inhibition was observed with a hypoxanthine substitution (AG → HG). When adenine was replaced by a purine base (AG → PG) or by 7-deazaadenine (AG → c7AG), little effect on the second step was observed, suggesting that the 6-NH2 and N7 groups do not play a critical role in adenine recognition. Finally, replacement of adenine by 2-aminopurine (AG → 2-APG) had no effect on the second step. Taken together, our results suggest that the N1 group of adenine functions as an essential determinant in adenine recognition during the second step of pre-mRNA splicing.