Constraints on Intron Evolution in the Gene Encoding the Myosin Alkali Light Chain in Drosophila

Interspecific comparisons of intron sequences reveal conserved blocks of invariant nucleotides and several other departures from the strictly neutral model of molecular evolution. To distinguish the past action of evolutionary forces in introns known to have regulatory information, we examined nucleotide sequence variation at 991 sites in a random sample of 16 Drosophila melanogaster alleles of the gene encoding the myosin alkali light chain (Mlc1). The Mlc1 gene of D. melanogaster encodes two MLC1 isoforms via developmentally regulated alternative pre-mRNA splicing. Analyses of these data reveal that introns 4 and 5, which flank the alternatively spliced exon 5, have reduced levels of both intraspecific polymorphism and interspecific divergence relative to intron 3. No polymorphism was observed in any of the exons examined in D. melanogaster. A genealogical analysis clearly demonstrates the occurrence of intragenic recombination in the ancestral history of Mlc1. Recombination events are estimated to be 13 times more likely than mutation events over the span of the sequenced region. Although there is little evidence for pairwise linkage disequilibrium in the Mlc1 region, higher order disequilibrium does seem to be present in the 5' half of the portion of the gene that was examined. Predictions of the folding free energy of the pre-mRNA reveal that sampled alleles have a significantly higher (less stable) free energy than do randomly permuted sequences. These results are consistent with the hypothesis that introns surrounding an alternatively spliced exon are subjected to additional constraints, perhaps due to specific aspects of secondary structure required for appropriate splicing of the pre-mRNA molecule.     

High rate of recent intron gain and loss in simultaneously duplicated Arabidopsis genes

We examined the gene structure of a set of 2563 Arabidopsis thaliana paralogous pairs that were duplicated simultaneously 20-60 MYA by tetraploidy. Out of a total of 23,164 introns in these genes, we found that 10,004 pairs have been conserved and 578 introns have been inserted or deleted in the time since the duplication event. This intron insertion/deletion rate of 2.7 x 10(-3) to 9.1 x 10(-4) per site per million years is high in comparison to previous studies. At least 56 introns were gained and 39 lost based on parsimony analysis of the phylogenetic distribution of these introns. We found weak evidence that genes undergoing intron gain and loss are biased with respect to gene ontology terms. Gene pairs that experienced at least 2 intron insertions or deletions show evidence of enrichment for membrane location and transport and transporter activity function. We do not find any relationship of intron flux to expression level or G + C content of the gene. Detection of a bias in the location of intron gains and losses within a gene depends on the method of measurement: an intragene method indicates that events (specifically intron losses) are biased toward the 3' end of the gene. Despite the relatively recent acquisition of these introns, we found only one case where we could identify the mechanism of intron origin--the TOUCH3 gene has experienced 2 tandem, partial, internal gene duplications that duplicated a preexisting intron and also created a novel, alternatively spliced intron that makes use of a duplicated pair of cryptic splice sites.     

Alternatively spliced p53 RNA in transformed and normal cells of different tissue types.

The alternatively spliced RNA species of tumor suppressor gene p53, containing an additional 96 bases derived from intron 10, is present at approximately 25 to 30% the level of regularly spliced p53 RNA in both normal epidermal and carcinoma cells. The presence of this alternatively spliced RNA in 10T1/2 fibroblast cells, mouse liver and testis suggests that this alternative splicing may be universal. The level of alternatively spliced p53 RNA was increased coordinately with that of regularly spliced p53 in 10T1/2 cells in response to epidermal growth factor. Immunoprecipitation analysis of epidermal cells using monoclonal antibodies which recognize different epitopes of p53 suggested that distinct p53 proteins may be translated from both RNA species. Considering previous observations on the potential importance of carboxyl terminal sequences in p53 function, knowledge of the ubiquitous presence of alternatively spliced p53 is important for future studies of p53 function in normal cells and in oncogenesis.     

Molecular population genetics of sequence length diversity in the Adh region of Drosophila pseudoobscura

Positive and negative selection on indel variation may explain the correlation between intron length and recombination levels in natural populations of Drosophila. A nucleotide sequence analysis of the 3.5 kilobase sequence of the alcohol dehydrogenase (Adh) region from 139 Drosophila pseudoobscura strains and one D. miranda strain was used to determine whether positive or negative selection acts on indel variation in a gene that experiences high levels of recombination. A total of 30 deletion and 36 insertion polymorphisms were segregating within D. pseudoobscura populations and no indels were fixed between D. pseudoobscura and its two sibling species D. miranda and D. persimilis. The ratio of Tajima's D to its theoretical minimum value (D(min)) was proposed as a metric to assess the heterogeneity in D among D. pseudoobscura loci when the number of segregating sites differs among loci. The magnitude of the D/D(min) ratio was found to increase as the rate of population expansion increases, allowing one to assess which loci have an excess of rare variants due to population expansion versus purifying selection. D. pseudoobscura populations appear to have had modest increases in size accounting for some of the observed excess of rare variants. The D/D(min) ratio rejected a neutral model for deletion polymorphisms. Linkage disequilibrium among pairs of indels was greater than between pairs of segregating nucleotides. These results suggest that purifying selection removes deletion variation from intron sequences, but not insertion polymorphisms. Genome rearrangement and size-dependent intron evolution are proposed as mechanisms that limit runaway intron expansion.     

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.      

Repetitive Indel Markers within the ALMT1 Gene Conditioning Aluminium Tolerance in Wheat (Triticum aestivum L.)

ALMT1 gene encoding a membrane protein that facilitates an aluminium stimulated malate efflux has been characterised and mapped in wheat (Triticum aestivum L.). Here, we have identified molecular markers targeting insertion/deletion (indel) and SSR repeats within intron 3 region of the ALMT1 gene. Both the markers: ALMT1-SSR3a and ALMT1-SSR3b based on repetitive indels, exhibited complete cosegregation with Al tolerance, malate efflux, and a CAPS marker discriminating ALMT1-1 and ALMT1-2 alleles, in a doubled haploid population derived from Diamondbird (Al-tolerant)/Janz (Al-sensitive). A parental screen of 20 diverse wheat genotypes with repetitive indel markers indicated that six allele variants exist at the ALMT1SSR3 locus. Sequence analysis confirmed that these variations were due to indels, copy number of SSR repeats, and base substitution within SSR repeats. The higher level of variation in intron 3 suggests that this genomic region has been constrained by indels, SSR and single nucleotide polymorphisms. Results have proven that repetitive indel markers cosegregating with the Al tolerance locus will be useful for marker assisted selection and population and evolution studies.     

Selective constraints on intron evolution in Drosophila

Intron sizes show an asymmetrical distribution in a number of organisms, with a large number of "short" introns clustered around a minimal intron length and a much broader distribution of longer introns. In Drosophila melanogaster, the short intron class is centered around 61 bp. The narrow length distribution suggests that natural selection may play a role in maintaining intron size. A comparison of 15 orthologous introns among species of the D. melanogaster subgroup indicates that, in general, short introns are not under greater DNA sequence or length constraints than long introns. There is a bias toward deletions in all introns (deletion/insertion ratio is 1.66), and the vast majority of indels are of short length (<10 bp). Indels occurring on the internal branches of the phylogenetic tree are significantly longer than those occurring on the terminal branches. These results are consistent with a compensatory model of intron length evolution in which slightly deleterious short deletions are frequently fixed within species by genetic drift, and relatively rare larger insertions that restore intron length are fixed by positive selection. A comparison of paralogous introns shared among duplicated genes suggests that length constraints differ between introns within the same gene. The janusA, janusB, and ocnus genes share two short introns derived from a common ancestor. The first of these introns shows significantly fewer indels than the second intron, although the two introns show a comparable number of substitutions. This indicates that intron-specific selective constraints have been maintained following gene duplication, which preceded the divergence of the D. melanogaster species subgroup.     

Bovine spongiform encephalopathy associated insertion/deletion polymorphisms of the prion protein gene in the four beef cattle breeds from North China

Two insertion/deletion (indel) polymorphisms of the prion protein gene (PRNP), a 23-bp indel in the putative promoter region and a 12-bp indel within intron I, are associated with the susceptibility to bovine spongiform encephalopathy (BSE) in cattle. In the present study, the polymorphism frequencies of the two indels in four main beef cattle breeds (Hereford, Simmental, Black Angus, and Mongolian) from North China were studied. The results showed that the frequencies of deletion genotypes and alleles of 23- and 12-bp indels were lower, whereas the frequencies of insertion genotypes and alleles of the two indels were higher in Mongolian cattle than in the other three cattle breeds. In Mongolian cattle, the 23-bp insertion / 12-bp insertion was the major haplotype, whereas in Hereford, Simmental, and Black Angus cattle, the 23-bp deletion / 12-bp deletion was the major haplotype. These results demonstrated that Mongolian cattle could be more resistant to BSE, compared with the other three cattle breeds, because of its relatively low frequencies of deletion genotypes and alleles of 23- and 12-bp indel polymorphisms. Thus, this race could be important for selective breeding to improve resistance against BSE in this area.     

Sequence divergence at the putative flowering time locus COL1 in Brassicaceae

An insertion/deletion polymorphism (Ind2) in the Brassica nigra CONSTANS LIKE 1 (Bni COL1) gene was previously found to be associated with variation in flowering time. In the present study we examine the inter-specific divergence of COL1 in the family Brassicaceae. Analysis of codon substitution models did not reveal evidence of positive Darwinian selection, but comparisons of the COL1 gene in different species revealed a surprising number of indels. A total of 24 indels were found in the 650 bp of the middle variable region of the gene. This high number of indels could reflect a lack of constraint on length of this region of the protein, or the effect of positive selection. The number of indels was close to that expected in non-coding DNA, but the indels were longer in COL1 than those observed in non-coding regions. Reconstruction of indel evolution indicated that most indels resulted from deletions rather than insertions. The Ind2 indel that has shown association with flowering time in Brassica nigra exhibited a remarkable distribution in the Brassicaceae family, indicating that the polymorphism may have persisted more than ten million years. Considering presumed historic populations sizes of Brassicaceae species, such a long persistence time seems unlikely for a neutral polymorphism.     

Evolution of the mitochondrial rps3 intron in perennial and annual angiosperms and homology to nad5 intron 1.

The plant mitochondrial rps3 intron was analyzed for substitution and indel rate variation among 15 monocot and dicot angiosperms from 10 genera, including perennial and annual taxa. Overall, the intron sequence was very conserved among angiosperms. Based on length polymorphism, 10 different alleles were identified among the 10 genera. These allelic differences were mainly attributable to large indels. An insertion of 133 nucleotides, observed in the Alnus intron was partially or completely absent in the other lineages of the family Betulaceae. This insertion was located within domain IV of the secondary-structure model of this group IIA intron. A mobile element of 47 nucleotides that showed homology to sequences located in rice rps3 intron and in intergenic plant mitochondrial genomes was found within this insertion. Both substitution and indel rates were low among the Betulaceae sequences, but substitution rates were increasingly larger than indel rates in comparisons involving more distantly related taxa. From a secondary-structure model, regions involved in helical structures were shown to be well preserved from indels as compared to substitutions, but compensatory changes were not observed among the angiosperm sequences analyzed. Using approximate divergence times based on the fossil record, substitution and indel rate heterogeneity was observed between different pairs of annual and perennial taxa. In particular, the annual petunia and primrose evolved more than 15 and 10 times faster, for substitution and indel rates respectively, than the perennial birch and alder. This is the first demonstration of an evolutionary rate difference between perennial and annual forms in noncoding DNA, lending support to neutral causes such as the generation time, population size, and speciation rate effects to explain such rate heterogeneity. Surprisingly, the sequence from the rps3 intron had a high identity with the sequence of intron 1 from the angiosperm mitochondrial nad5 gene, suggesting a common origin of these two group IIA introns.     

Phylogenetic relationships among <Emphasis Type="Italic">Fagopyrum</Emphasis> species based on analysis of the <Emphasis Type="Italic">nad1</Emphasis> Gene b/c intron

For the first time, the b/c intron of the nad1 gene was characterized in 13 species of Fagopyrum and the applicability of the intron was demonstrated for identification of buckwheat species. The length of the b/c intron of the nad1 gene varied from 1217 bp in F. tataricum and F. cymosum to 1239 bp in F. capillatum. Species-specific SNPs and indels were identified. The genus phylogeny that was determined with the molecular approach correlated with phylogenetic data that are based on morphological characteristics. The NJ and MP dendrograms revealed distinct division of the analyzed species into two main groups, cymosum and urophyllum. The species of the urophyllum group were characterized by the presence of several synapomorphic substitutions and indels, including the insertion of 11 nucleotides, which is flanked by direct repeats.