Length variation and secondary structure of introns in the Mlc1 gene in six species of Drosophila

A nearly universal feature of intron sequences is that even closely related species exhibit a large number of insertion/deletion differences. The goal of the analysis described here is to test whether the observed pattern of insertion/deletion events in the genealogy of the myosin alkali light chain (Mlc1) gene is consistent with neutrality, and if not, to determine the underlying forces of evolutionary change. Mlc1 pre-mRNA is alternatively spliced, and one constraint is that signals necessary for tissue-specificity of directed splicing must be conserved. If the total length of an intron is functionally constrained, then the distribution of indels on branches of the gene genealogy should reflect a departure from randomness. Here we perform a phylogenetic analysis, inferring ancestral states wherever possible on a phylogeny of 29 alleles of Mlc1 from six species of Drosophila. Observed patterns of indels on the genealogy were compared to those from simulated data, with the result that we cannot reject the null hypothesis of neutrality. A clear departure from a neutral prediction was seen in the excess folding free energy predicted for the introns flanking the alternatively spliced exon. Relative rate tests also suggest a retardation in the rate of Mlc1 sequence evolution in the simulans clade.      

Conservation of alternative splicing and genomic organization of the myosin alkali light-chain (Mlc1) gene among Drosophila species

The Mlc1 gene of Drosophila melanogaster encodes two MLC1 isoforms via developmentally regulated alternative pre-mRNA splicing. In larval muscle and tubular and abdominal muscles of adults, all of the six exons are included in the spliced mRNA, whereas, in the fibrillar indirect flight muscle of adult, exon 5 is excluded from the mRNA. We show that this tissue-specific pattern of alternative splicing of the Mlc1 pre-mRNA is conserved in D. simulans, D. pseudoobscura, and D. virilis. Isolation and sequencing of the Mlc1 genes from these three other Drosophila species have revealed that the overall organization of the genes is identical and that the genes have maintained a very high level of sequence identity within the coding region. Pairwise amino acid identities are 94%-99%, and there are no charge changes among the proteins. Total nucleotide divergence within the coding region of the four genes supports the accepted genealogy of these species, but the data indicate a significantly higher rate of amino acid replacement in the branch leading to D. pseudoobscura. A comparison of nucleotide substitutions in the coding portions of exon 5 and exon 6, which encode the alternative carboxyl termini of the two MLC1 isoforms, suggests that exon 5 is subject to greater evolutionary constraints than is exon 6. In addition to the coding sequences, there is significant sequence conservation within the 5' and 3' noncoding DNA and two of the introns, including one that flanks exon 5. These regions are candidates for cis- regulatory elements. Our results suggest that evolutionary constraints are acting on both the coding and noncoding sequences of the Mlc1 gene to maintain proper expression and function of the two MLC1 polypeptides.      

The mechanism of sex-specific splicing at the doublesex gene is different between Drosophila melanogaster and Bombyx mori

We have previously reported that Bmdsx, a homologue of the sex-determining gene, doublesex (dsx), was found to be sex-specifically expressed in various tissues at larval, pupal, and adult stages in the silkworm, Bombyx mori, and was alternatively spliced to yield male- and female-specific mRNAs. To reveal sex-specific differences in splicing patterns of Bmdsx pre-mRNA, the genomic sequence was determined and compared with male- and female-specific Bmdsx cDNA sequences. The open reading frame (ORF) consisted of five exons. Exons 3 and 4 were specifically incorporated into the female type of Bmdsx mRNA. On the other hand, exon 2 was spliced to exon 5 to produce the male type mRNA of Bmdsx. As in the case of Drosophila dsx, the OD2 domain was separated by a female-specific intron into sex-independent and sex-dependent regions. Sex-specific splicing occurred in equivalent positions in the Drosophila dsx gene. However, unlike Drosophila dsx, the female-specific introns showed no weak 3′ splice sites, and the TRA/TRA-2 binding site related sequences were not found in the female-specific exon, nor even in any other regions of the Bmdsx gene. Moreover, an in vitro splicing reaction consisting of HeLa cell nuclear extracts showed that the female-type of Bmdsx mRNA represented the default splicing. These findings suggest that the structural features of the sex-specific splicing patterns of Bmdsx pre-mRNA are similar to those of Drosophila dsx but the regulation of sex-specific alternative splicing of Bmdsx pre-mRNA is different.     

Alternative splicing of beta-tropomyosin pre-mRNA: multiple cis-elements can contribute to the use of the 5'- and 3'-splice sites of the nonmuscle/smooth muscle exon 6.

We previously found that the splicing of exon 5 to exon 6 in the rat beta-TM gene required that exon 6 first be joined to the downstream common exon 8 (Helfman et al., Genes and Dev. 2, 1627-1638, 1988). Pre-mRNAs containing exon 5, intron 5 and exon 6 are not normally spliced in vitro. We have carried out a mutational analysis to determine which sequences in the pre-mRNA contribute to the inability of this precursor to be spliced in vitro. We found that mutations in two regions of the pre-mRNA led to activation of the 3'-splice site of exon 6, without first joining exon 6 to exon 8. First, introduction of a nine nucleotide poly U tract upstream of the 3'-splice site of exon 6 results in the splicing of exon 5 to exon 6 with as little as 35 nucleotides of exon 6. Second, introduction of a consensus 5'-splice site in exon 6 led to splicing of exon 5 to exon 6. Thus, three distinct elements can act independently to activate the use of the 3'-splice site of exon 6: (1) the sequences contained within exon 8 when joined to exon 6, (2) a poly U tract in intron 5, and (3) a consensus 5'-splice site in exon 6. Using biochemical assays, we have determined that these sequence elements interact with distinct cellular factors for 3'-splice site utilization. Although HeLa cell nuclear extracts were able to splice all three types of pre-mRNAs mentioned above, a cytoplasmic S100 fraction supplemented with SR proteins was unable to efficiently splice exon 5 to exon 6 using precursors in which exon 6 was joined to exon 8. We also studied how these elements contribute to alternative splice site selection using precursors containing the mutually exclusive, alternatively spliced cassette comprised of exons 5 through 8. Introduction of the poly U tract upstream of exon 6, and changing the 5'-splice site of exon 6 to a consensus sequence, either alone or in combination, facilitated the use of exon 6 in vitro, such that exon 6 was spliced more efficiently to exon 8. These data show that intron sequences upstream of an exon can contribute to the use of the downstream 5'-splice, and that sequences surrounding exon 6 can contribute to tissue-specific alternative splice site selection.     

The mechanism of sex-specific splicing at the doublesex gene is different between Drosophila melanogaster and Bombyx mori

We have previously reported that Bmdsx, a homologue of the sex-determining gene, doublesex (dsx), was found to be sex-specifically expressed in various tissues at larval, pupal, and adult stages in the silkworm, Bombyx mori, and was alternatively spliced to yield male- and female-specific mRNAs. To reveal sex-specific differences in splicing patterns of Bmdsx pre-mRNA, the genomic sequence was determined and compared with male- and female-specific Bmdsx cDNA sequences. The open reading frame (ORF) consisted of five exons. Exons 3 and 4 were specifically incorporated into the female type of Bmdsx mRNA. On the other hand, exon 2 was spliced to exon 5 to produce the male type mRNA of Bmdsx. As in the case of Drosophila dsx, the OD2 domain was separated by a female-specific intron into sex-independent and sex-dependent regions. Sex-specific splicing occurred in equivalent positions in the Drosophila dsx gene. However, unlike Drosophila dsx, the female-specific introns showed no weak 3′ splice sites, and the TRA/TRA-2 binding site related sequences were not found in the female-specific exon, nor even in any other regions of the Bmdsx gene. Moreover, an in vitro splicing reaction consisting of HeLa cell nuclear extracts showed that the female-type of Bmdsx mRNA represented the default splicing. These findings suggest that the structural features of the sex-specific splicing patterns of Bmdsx pre-mRNA are similar to those of Drosophila dsx but the regulation of sex-specific alternative splicing of Bmdsx pre-mRNA is different.     

A highly stable duplex structure sequesters the 5' splice site region of hnRNP A1 alternative exon 7B.

Exon 7B in the hnRNP A1 pre-mRNA is alternatively spliced to yield A1 and A1(B), two proteins that differ in their ability to modulate 5' splice site selection. Sequencing the murine intron downstream of exon 7B revealed the existence of several regions of similarity to the corresponding human intron. In vitro splicing assays indicate that an 84-nt region (CE6IO) decreases splicing to the proximal 5' splice site in a pre-mRNA carrying the 5' splice sites of exon 7 and 7B. In vivo, the CE6IO element promotes exon 7B skipping in pre-mRNAs expressed from a mini-gene containing the hnRNP A1 alternative splicing unit. Using oligonucleotide-targeted RNase H cleavage assays, we provide support for the existence of highly stable base pairing interactions between CE6IO and the 5' splice site region of exon 7B. Duplex formation occurs in naked pre-mRNA, resists incubation in splicing extracts, and is associated with a reduction in the assembly of U1 snRNP-dependent complexes to the 5' splice site of exon 7B. Our results demonstrate that pre-mRNA secondary structure plays an important role in promoting exon 7B skipping in the A1 pre-mRNA.     

Intronic alternative splicing regulators identified by comparative genomics in nematodes

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high- scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.     

Cis requirements for alternative splicing of the cardiac troponin T pre-mRNA.

The cardiac troponin T (cTNT) pre-mRNA splices 17 exons contiguously but alternatively splices (includes or excludes) the fifth exon. Because both alternative splice products are processed from the same pre-mRNA species, the cTNT pre-mRNA must contain cis-acting sequences which specify exon 5 as an alternative exon. A cTNT minigene (SM-1) transfected into cultured cells produces mRNAs both including and excluding exon 5. The junctions of exons 4-5-6 and 4-6 in the cTNT minigene mRNAs are identical to those of endogenous cTNT mRNAs and no other exons are alternatively spliced. Thus, the SM-1 pre-mRNA is correctly alternatively spliced in transfected cells. To circumscribe the pre-mRNA regions which are required for the alternative nature of exon 5, we have constructed a systematic series of deletion mutants of SM-1. Transfection of this series demonstrates that a 1200 nt pre-mRNA region containing exons 4, 5, and 6 is sufficient to direct alternative splicing of exon 5. Within this region are two relatively large inverted repeats which potentially sequester the alternative exon via intramolecular base-pairing. Such sequestration of an alternative exon is consistent with models which propose pre-mRNA conformation as being determinative for alternative splicing of some pre-mRNAs. However, deletion mutants which remove the majority of each of the inverted repeats retain the ability to alternatively splice exon 5 demonstrating that neither is required for cTNT alternative splice site selection. Taken together, deletion analysis has limited cis elements required for alternative splicing to three small regions of the pre-mRNA containing exons 4, 5, and 6. In addition, the cTNT minigene pre-mRNA expresses both alternative splice products in a wide variety of cultured non-muscle cells as well as in cultured striated muscle cells, although expression of the cTNT pre-mRNA is normally restricted to striated muscle. This indicates that cis elements involved in defining the cTNT exon 5 as an alternative exon do not require muscle-specific factors in trans to function.     

Rapid changes in the exon/intron structure of a mammalian thrombin inhibitor gene

The genomic organization of the heparin cofactor II (HCII) gene from rat and mouse was investigated and compared with their human counterpart. The genes share a common core structure consisting of five exons interrupted by four introns, but the mouse and rat gene reveal individual additional features. A unique differentially spliced exon is present in the 5'-untranslated region of the rat gene, which most probably has arisen de novo by point mutations in intronic sequences of the ancestor gene. In the mouse HCII gene, a novel intron/exon boundary has been created due to the presence of an additional DNA segment, which simultaneously provides a 3'-splice site and a polypyrimidine stretch leading to an alternatively used exon of increased size. Our data suggest that, in contrast to most other mammalian genes, the exon/intron pattern of the gene coding for HCII is in dynamic evolution.     

A combinatorial role for exon, intron and splice site sequences in splicing in maize

Plant introns are typically AU-rich or U-rich, and this feature has been shown to be important for splicing. In maize, however, about 20% of the introns exceed 50% GC, and most of them are efficiently spliced. A series of constructs has been designed to analyze the cis requirements for splicing of the GC-rich Bz2 maize intron and two other GC-rich intron derivatives. By manipulating exon, intron and splice site sequences it is shown that exons can play an important role in intron definition: changes in exon sequences can increase splicing efficiency of a GC-rich intron from 17% to 86%. The relative difference, or base compositional contrast, in GC and U content between exon and intron sequences in the vicinity of splice sites, rather than the absolute base-content of the intron or exons, correlates with splicing efficiency. It is also shown that GC-rich intron constructs that are poorly spliced can be partially rescued by an improved 3' splice site.     

Mutation of putative branchpoint consensus sequences in plant introns reduces splicing efficiency

Intron lariat formation between the 5' end of an intron and a branchpoint adenosine is a fundamental aspect of the first step in animal and yeast nuclear pre-mRNA splicing. Despite similarities in intron sequence requirements and the components of splicing, differences exist between the splicing of plant and vertebrate introns. The identification of AU-rich sequences as major functional elements in plant introns and the demonstration that a branchpoint consensus sequence was not required for splicing have led to the suggestion that the transition from AU-rich intron to GC-rich exon is a major potential signal by which plant pre-mRNA splice sites are recognized. The role of putative branchpoint sequences as an internal signal in plant intron recognition/definition has been re-examined. Single nucleotide mutations in putative branchpoint adenosines contained within CUNAN sequences in four different plant introns all significantly reduced splicing efficiency. These results provide the most direct evidence to date for preferred branchpoint sequences being required for the efficient splicing of at least some plant introns in addition to the important role played by AU sequences in dicot intron recognition. The observed patterns of 3' splice site selection in the introns studied are consistent with the scanning model described for animal intron 3' splice site selection. It is suggested that, despite the clear importance of AU sequences for plant intron splicing, the fundamental processes of splice site selection and splicing in plants are similar to those in animals.     

Aberrant pre-mRNA maturation is caused by LINE insertions into introns of the white gene of Drosophila melanogaster.

Insertional mutagenesis screens have provided thousands of mutant alleles for analysing genes of varied functions in Drosophila melanogaster. We here document mechanisms of insertional mutagenesis by a LINE element, the I factor, by determining the molecular structure of RNAs produced from two alleles of the white gene of D.melanogaster, wIR1 and wIR6. These alleles result from insertion of the I factor into introns of the gene. We show that sequences present within the element direct aberrant splicing and termination events. When the I factor is inserted within the white first intron it may lead to the use of a cryptic 3' splice site which does not contain the dinucleotide AG. This splicing gives rise to a chimeric messenger RNA whose synthesis is controlled differently in tissues where the mutated gene is expressed. When the I factor is inserted within the white last intron it induces synthesis of truncated mRNAs. These results provide, for the first time, mechanisms for I factor insertional mutagenesis. They are discussed in the more general context of RNA processing in Drosophila and the evolution of eukaryotic gene introns.     

Molecular variation of the shuttle craft and Lim3 genes, controlling the development of the nervous system, in a natural Drosophila melanogaster population

Comparative polymorphism of the first exon and first intron of the shuttle craft (stc) and Lim3 genes and their putative regulatory 5'-flanking sequences was analyzed using 20 sequenced natural alleles. A comparison of the stc and Lim3 genes showed that the extent of polymorphism was similar in their introns and corresponded to the variation level characteristic of Drosophila melanogaster, while the putative regulatory region and first intron of the stc gene proved to be more variable than the corresponding regions of the Lim3 gene. Since the genes under study occurred on the same chromosomes isolated from one population and were close together in a region having a high recombination rate, the difference in the extent of polymorphism between the regulatory and coding regions was explained by individual characteristics of each gene. The results made it possible to assume that the extent of polymorphism of the coding gene regions is maintained by balancing selection.     

Retention of spliceosomal components along ligated exons ensures efficient removal of multiple introns

The majority of mammalian pre-mRNAs contains multiple introns that are excised prior to export and translation. After intron excision, ligated exon intermediates participate in subsequent intron excisions. However, exon ligation generates an exon of increased size, a feature of pre-mRNA splicing that can interfere with downstream splicing events. These considerations raise the question of whether unique mechanisms exist that permit efficient removal of introns neighboring ligated exons. Kinetic analyses of multiple intron-containing pre-mRNAs revealed that splicing is more efficient following an initial intron removal event, suggesting that either the recruitment of the exon junction complex (EJC) to ligated exons increases the efficiency of multiple intron excisions or that the initial definition of splice sites is sufficient to permit efficient splicing of introns neighboring ligated exons. Knockdown experiments show that the deposition of the EJC does not affect subsequent splicing kinetics. Instead, spliceosomal components that are not involved in the initial splicing event remain associated with the pre-mRNA to ensure efficient removal of neighboring introns. Thus, ligated exons do not require redefinition, providing an additional kinetic advantage for exon defined splice sites.     

Ovomucoid intervening sequences specify functional domains and generate protein polymorphism

Two thirds of the natural chicken ovomucoid gene has been sequenced, including all exons and the intron sequences surrounding all fourteen intron/ exon junctions. The junction sequences surrounding four of the introns are redundant; however, the sequences surrounding the other three introns contain no redundancies and thus the splicing sites at either end of these three introns are unambiguous. The splicing in all cases conforms to the GT-AG rule. The ovomucoid gene sequence around intron F can be used to predict the cause of an internal deletion polymorphism in the ovomucoid protein, which is an apparent error in the processing of the ovomucoid pre-mRNA. We also compare the structural organization of the ovomucoid gene with the ovomucoid protein sequence to examine theories of the evolution of ovomucoids as well as the origin of intervening sequences. This analysis suggests that the present ovomucoid gene evolved from a primordial ovomucoid gene by two separate intragenic duplications. Furthermore, sequence analyses suggest that introns were present in the primordial ovomucoid gene before birds and mammals diverged, about 300 million years ago. Finally, the positions of the introns within the ovomucoid gene support the theory that introns separate gene segments that code for functional domains of proteins and provide insight on the manner by which eucaryotic genes were constructed during the process of evolution.     

Combinatorial splicing of exon pairs by two-site binding of U1 small nuclear ribonucleoprotein particle.

A two-site model for the binding of U1 small nuclear ribonucleoprotein particle (U1 snRNP) was tested in order to understand how exon partners are selected in complex pre-mRNAs containing alternative exons. In this model, it is proposed that two U1 snRNPs define a functional unit of splicing by base pairing to the 3' boundary of the downstream exon as well as the 5' boundary of the intron to be spliced. Three-exon substrates contained the alternatively spliced exon 4 (E4) region of the preprotachykinin gene. Combined 5' splice site mutations at neighboring exons demonstrate that weakened binding of U1 snRNP at the downstream site and improved U1 snRNP binding at the upstream site result in the failure to rescue splicing of the intron between the mutations. These results indicate the stringency of the requirement for binding a second U1 snRNP to the downstream 5' splice site for these substrates as opposed to an alternative model in which a certain threshold level of U1 snRNP can be provided at either site. Further support for the two-site model is provided by single-site mutations in the 5' splice site of the third exon, E5, that weaken base complementarity to U1 RNA. These mutations block E5 branchpoint formation and, surprisingly, generate novel branchpoints that are specified chiefly by their proximity to a cryptic 5' splice site located at the 3' terminus of the pre-mRNA. The experiments shown here demonstrate a true stimulation of 3' splice site activity by the downstream binding of U1 snRNP and suggest a possible mechanism by which combinatorial patterns of exon selection are achieved for alternatively spliced pre-mRNAs.