How intron splicing affects the deletion and insertion profile in Drosophila melanogaster

Studies of "dead-on-arrival" transposable elements in Drosophila melanogaster found that deletions outnumber insertions approximately 8:1 with a median size for deletions of approximately 10 bp. These results are consistent with the deletion and insertion profiles found in most other Drosophila pseudogenes. In contrast, a recent study of D. melanogaster introns found a deletion/insertion ratio of 1.35:1, with 84% of deletions being shorter than 10 bp. This discrepancy could be explained if deletions, especially long deletions, are more frequently strongly deleterious than insertions and are eliminated disproportionately from intron sequences. To test this possibility, we use analysis and simulations to examine how deletions and insertions of different lengths affect different components of splicing and determine the distribution of deletions and insertions that preserve the original exons. We find that, consistent with our predictions, longer deletions affect splicing at a much higher rate compared to insertions and short deletions. We also explore other potential constraints in introns and show that most of these also disproportionately affect large deletions. Altogether we demonstrate that constraints in introns may explain much of the difference in the pattern of deletions and insertions observed in Drosophila introns and pseudogenes.     

Molecular Evolution of the Metallothionein Gene Mtn in the Melanogaster Species Group: Results from Drosophila Ananassae

Three distinctly different alleles of the metallothionein gene Mtn have been identified in natural Drosophila melanogaster populations: Mtn(.3), Mtn(1), and Dp(Mtn(1)), where the latter designates a tandem duplication of Mtn(1). In Drosophila simulans, only Mtn(.3)-type alleles have been found. It has been suggested that Mtn(.3) is the ancestral allele and demonstrated that a presumed two-step transition from Mtn(.3) to Mtn(1) to Dp(Mtn(1)) is accompanied by an approximate 5-fold increase in RNA levels. We analyzed the evolutionary genetics of the Mtn locus of Drosophila ananassae, a distant relative of D. melanogaster and D. simulans within the melanogaster species group. The Mtn gene of D. ananassae is most similar to Mtn(.3). (i) it is identical with Mtn(.3) at the amino acid level, but differs from Mtn(1) in its terminal codon; (ii) its 3'' UTR contains a characteristic extra DNA segment of about 50 bp which is present in Mtn(.3), but lacking in Mtn(1); (iii) duplications of Mtn were not found in a worldwide sample of 110 wild D. ananassae chromosomes. However, the intron of the Mtn gene in D. ananassae is only 69 bp long, whereas the length of the Mtn(.3) and Mtn(1) introns is 265 bp; and it lacks a polypyrimidine stretch upstream of the 3'' splice site in contrast to the much greater pyrimidine-richness found in the Mtn(.3) and Mtn(1) introns. A short intron (67 bp) was also identified in a D. pseudoobscura Mtn allele, suggesting that the short intron is the ancestral form and that the transition from the short to the long intron occurred within the melanogaster species group. We discuss the significance of this observation with regard to the recently proposed classification of D. melanogaster introns into two groups: short introns (<90 bp) which tend to lack polypyrimidine stretches, and longer ones which have strong 3'' splice signals similar to mammalian introns. A database search revealed that this length dimorphism is an evolutionarily conserved feature of Drosophila introns; transitions from one size class to the other appear to be rare between closely related species (e.g., within the melanogaster subgroup).     

Growth and decline of introns

To gauge the processes that might direct the length of introns, I studied the balance of indels (insertions or deletions, determined using Alu and LINE1 retroposon repeats) and the density of these repeats in the introns of the human genome. The indel balance is biased in favour of deletions and correlated with the divergence of repeats. At fixed repeat divergence, the indel bias correlated with the intron size: the shorter the intron, the more deletions were favoured over insertions. This correlation with the intron size was stronger than with the gene-wide or isochore-wide parameters. The density of repeats (the number of repeats in a unit of intron length) correlated positively with the intron size. Thus, quite different mechanisms, the indel bias and the integration and/or persistence of retroposons, act in the same direction in regards to intron size, which suggests selection for the size of individual introns.     

The correlation between intron length and recombination in drosophila. Dynamic equilibrium between mutational and selective forces

Intron length is negatively correlated with recombination in both Drosophila melanogaster and humans. This correlation is not likely to be the result of mutational processes alone: evolutionary analysis of intron length polymorphism in D. melanogaster reveals equivalent ratios of deletion to insertion in regions of high and low recombination. The polymorphism data do reveal, however, an excess of deletions relative to insertions (i.e., a deletion bias), with an overall deletion-to-insertion events ratio of 1.35. We propose two types of selection favoring longer intron lengths. First, the natural mutational bias toward deletion must be opposed by strong selection in very short introns to maintain the minimum intron length needed for the intron splicing reaction. Second, selection will favor insertions in introns that increase recombination between mutations under the influence of selection in adjacent exons. Mutations that increase recombination, even slightly, will be selectively favored because they reduce interference among selected mutations. Interference selection acting on intron length mutations must be very weak, as indicated by frequency spectrum analysis of Drosophila intron length polymorphism, making the equilibrium for intron length sensitive to changes in the recombinational environment and population size. One consequence of this sensitivity is that the advantage of longer introns is expected to decrease inversely with the rate of recombination, thus leading to a negative correlation between intron length and recombination rate. Also in accord with this model, intron length differs between closely related Drosophila species, with the longest variant present more often in D. melanogaster than in D. simulans. We suggest that the study of the proposed dynamic model, taking into account interference among selected sites, might shed light on many aspects of the comparative biology of genome sizes including the C value paradox.     

Polymorphism and structural variation of rps16 group-II intron in the Solanum species

Nucleotide sequences of the self-splicing group-II intron of rps16 have first been determined in nine species of the Solanum genus. It was found that the observed variations in the intron length (855–864 bp) was associated with indels of 1 to 9 bp. Altogether, five indels and 50 nucleotide substitutions were detected, which were used to identify six Solanum haplotypes. Although the intron sequence was in general fairly well conserved, the distribution of the described mutations among its structural elements corresponding to six pre-RNA domains was qualitatively and quantitatively nonuniform. The highest polymorphism levels were observed in domains I, II, and IV. The sequence of domain V was absolutely invariable, which is in agreement with its functional significance. The chloroplast rpS16 intron sequences have been characterized in nine Solanum species. The intron length ranged from 855 bp to 864 bp, which is associated with 1–9-nucleotide indels. In total five indels and 50 nucleotide substitutions have been detected and six Solanum haplotypes have been revealed. Solanum rpS16 introns has been characterized by mutation rate heterogeneity between structure regions of all six domains its pre-RNA. Intron domains I, II, IV are shown to be more variable. Sequences of the domain V are invariant, that agrees with its functional significance.     

Heterogeneous yet similar introns reside in identical positions of the rRNA genes in natural isolates of the archaeon Aeropyrum pernix

Some archaeal ribosomal DNA (rDNA) introns carry homing endonuclease-like genes and are therefore assumed to propagate by "intron homing". A previous study demonstrated that three introns are located within the rRNA operon (arnSL) of Aeropyrum pernix strain K1, two of which, Ialpha and Igamma, harbor open reading frames (ORFs) encoding putative LAGLIDADG-type endonucleases. In an effort to understand further the rDNA intron distribution in natural A. pernix populations, 11 A. pernix strains were isolated from marine hydrothermal biotopes, and comparative nucleotide sequence analysis of the arnSL alleles was performed. Of the 11 isolates, eight contained multiple introns, and three patterns of intron insertion were found. Three novel introns, Idelta (62 bp in length), Ivarepsilon (122 bp) and Izeta (57 bp) were identified. They were all ORF-less, but their predicted RNA secondary structure at the exon-intron junctions was consistent with the bulge-helix-bulge motif. The insertion positions and the terminal inverted repeat sequences of Idelta and Izeta were in agreement with those of Ialpha and Igamma, respectively. This suggests that these intron variants were generated by large indels (insertions/deletions) during their evolution.     

Complex selection on intron size in Cryptococcus neoformans

We conducted a genome-wide analysis of the roles of mutation and selection in sculpting intron size in the fungal pathogen Cryptococcus neoformans. We find that deletion rate is positively associated with intron length and that insertion rate exhibits a weak negative association with intron length. These patterns suggest that long introns as well as extremely short introns in this unusually intron-rich fungal genome are in mutation-selection disequilibrium and that the proportion of constrained functional sequence in introns does not scale linearly with size. We find that untranslated region introns are longer than coding-region introns and that first introns are substantially longer than subsequent introns, suggesting heterogeneous distribution of constrained functional sequence and/or selective pressures on intron size within genes. In contrast to Drosophila, we find a positive correlation between d(N) and first intron or last intron length and a negative correlation between d(N) and internal intron length. This contrasting pattern may indicate that terminal introns and internal introns are differentially subject to hypothesized selection pressures modulating intron size and provides further evidence of widespread selective constraints on noncoding sequences.     

Distributions of exons and introns in the human genome

The human genome is revisited using exon and intron distribution profiles. The 26,564 annotated genes in the human genome (build October, 2003) contain 233,785 exons and 207,344 introns. On average, there are 8.8 exons and 7.8 introns per gene. About 80% of the exons on each chromosome are < 200 bp in length. < 0.01% of the introns are < 20 bp in length and < 10% of introns are more than 11,000 bp in length. These results suggest constraints on the splicing machinery to splice out very long or very short introns and provide insight to optimal intron length selection. Interestingly, the total length in introns and intergenic DNA on each chromosome is significantly correlated to the determined chromosome size with a coefficient of correlation r = 0.95 and r = 0.97, respectively. These results suggest their implication in genome design.     

The genomic landscape of short insertion and deletion polymorphisms in the chicken (Gallus gallus) Genome: a high frequency of deletions in tandem duplicates

It is increasingly recognized that insertions and deletions (indels) are an important source of genetic as well as phenotypic divergence and diversity. We analyzed length polymorphisms identified through partial (0.25x) shotgun sequencing of three breeds of domestic chicken made by the International Chicken Polymorphism Map Consortium. A data set of 140,484 short indel polymorphisms in unique DNA was identified after filtering for microsatellite structures. There was a significant excess of tandem duplicates at indel sites, with deletions of a duplicate motif outnumbering the generation of duplicates through insertion. Indel density was lower in microchromosomes than in macrochromosomes, in the Z chromosome than in autosomes, and in 100 bp of upstream sequence, 5'-UTR, and first introns than in intergenic DNA and in other introns. Indel density was highly correlated with single nucleotide polymorphism (SNP) density. The mean density of indels in pairwise sequence comparisons was 1.9 x 10(-4) indel events/bp, approximately 5% the density of SNPs segregating in the chicken genome. The great majority of indels involved a limited number of nucleotides (median 1 bp), with A-rich motifs being overrepresented at indel sites. The overrepresentation of deletions at tandem duplicates indicates that replication slippage in duplicate sequences is a common mechanism behind indel mutation. The correlation between indel and SNP density indicates common effects of mutation and/or selection on the occurrence of indels and point mutations.     

Indel evolution of mammalian introns and the utility of non-coding nuclear markers in eutherian phylogenetics

Nuclear DNA intron sequences are increasingly used to investigate evolutionary relationships among closely related organisms. The phylogenetic usefulness of intron sequences at higher taxonomic levels has, however, not been firmly established and very few studies have used these markers to address evolutionary questions above the family level. In addition, the mechanisms driving intron evolution are not well understood. We compared DNA sequence data derived from three presumably independently segregating introns (THY, PRKC I and MGF) across 158 mammalian species. All currently recognized extant eutherian mammalian orders were included with the exception of Cingulata, Dermoptera and Scandentia. The total aligned length of the data was 6366 base pairs (bp); after the exclusion of autapomorphic insertions, 1511 bp were analyzed. In many instances the Bayesian and parsimony analyses were complementary and gave significant posterior probability and bootstrap support (>80) for the monophyly of Afrotheria, Euarchontoglires, Laurasiatheria and Boreoeutheria. Apart from finding congruent support when using these methods, the intron data also provided several indels longer than 3 bp that support, among others, the monophyly of Afrotheria, Paenungulata, Ferae and Boreoeutheria. A quantitative analysis of insertions and deletions suggested that there was a 75% bias towards deletions. The average insertion size in the mammalian data set was 16.49 bp +/- 57.70 while the average deletion was much smaller (4.47 bp +/- 14.17). The tendency towards large insertions and small deletions is highlighted by the observation that out of a total of 17 indels larger than 100 bp, 15 were insertions. The majority of indels (>60% of all events) were 1 or 2 bp changes. Although the average overall indel substitution rate of 0.00559 per site is comparable to that previously reported for rodents and primates, individual analyses among different evolutionary lineages provide evidence for differences in the formation rate of indels among the different mammalian groups.     

果蝇基因组中内含子数目随其长度的分布研究

为了研究内含子可能储存的有关生命进化信息,本研究以果蝇基因组中的内含子为样本,统计了内含子在不同长度区间内的相对频数,得到了内含子随其长度的分布,发现这种分布呈现出一定的规律性,即长度在1~80 bp内的内含子数目随其长度增加而增加,长度大于80 bp的内含子数目随其长度的增加而减少;本研究推论这种分布规律应该与某种分布模型一致,经过各种分布的拟合,最后发现这种分布与Γ(α,β)分布相符合;另外,将果蝇内含子分布规律和相应的外显分布进行了比较分析,发现虽然两类序列的分布规律都符合Γ(α,β)分布,但也存在明显的区别,并就此讨论了两类序列在生物进化中的意义。     

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.     

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.      

Optimal intron analyses in the trimeresurus radiation of Asian pitvipers

Nuclear introns are commonly used as phylogenetic markers, but a number of issues related to alignment strategies, indel treatments, and the incorporation of length-variant heterozygotes (LVHs) are not routinely addressed when generating phylogenetic hypotheses. Topological congruence in relation to an extensive mitochondrial DNA multigene phylogeny (derived from 2,423 bp of 12S, 16S, ND4, and CYTB genes) of the Asian pitviper Trimeresurus radiation was used to compare combinations of "by eye" and edited and unedited ClustalX 1.8 alignments of two nuclear introns. Indels were treated as missing data, fifth character states, and assigned simple and multistate codes. Upon recovery of the optimal alignment and indel treatment strategy, a total evidence approach was used to investigate the phylogenetic utility of the indels and test new generic arrangements within Trimeresurus. Approximately one third of the intron data partitions exhibited LVHs, suggesting that they are common in introns. Furthermore, a simple concatenation approach can facilitate the incorporation of LVHs into phylogenetic analyses to make use of all available data and investigate mechanisms of molecular evolution. Analyses of ClustalX 1.8-assisted alignments were generally more congruent than the "by eye" alignment and the analysis of a simple coded, edited ClustalX 1.8 (gap opening cost 5, gap extension cost 1) alignment revealed the most congruent tree. The total evidence approach supported the new arrangements within Trimeresurus, suggesting that the phylogeny should be considered as a working benchmark in Asian pitviper systematics. Finally, a critical appraisal of the diverse array of indels (56 to 57 per intron, ranging from 1 to 151 bp in length) suggested that they are a combination of Hennigian and homoplasious events unrelated to indel size or location within the intron. [Alignment; indels; intron analysis; length-variant heterozygotes; Trimeresurus.].     

Distribution and evolution of introns in drosophila amylase genes

While the two amylase genes of Drosophila melanogaster are intronless, the three genes of D. pseudoobscura harbor a short intron. This raises the question of the common structure of the Amy gene in Drosophila species. We have investigated the presence or absence of an intron in the amylase genes of 150 species of Drosophilids. Using polymerase chain reaction (PCR), we have amplified a region that surrounds the intron site reported in D. pseudoobscura and a few other species. The results revealed that most species contain an intron, with a variable size ranging from 50 to 750 bp, although the very majoritary size was around 60–80 bp. Several species belonging to different lineages were found to lack an intron. This loss of intervening sequence was likely due to evolutionarily independent and rather frequent events. Some other species had both types of genes: In the obscura group, and to a lesser extent in the ananassae subgroup, intronless copies had much diverged from intron-containing genes. Base composition of short introns was found to be variable and correlated with that of the surrounding exons, whereas long introns were all A-T rich. We have extended our study to non-Drosophilid insects. In species from other orders of Holometaboles, Lepidoptera and Hymenoptera, an intron was found at an identical position in the Amy gene, suggesting that the intron was ancestral. Received: 23 October 1995 / Accepted: 5 March 1996     

Nucleotide diversity and molecular evolution of the WAG-2 gene in common wheat (Triticum aestivum L) and its relatives

In this work, we examined the genetic diversity and evolution of the WAG-2 gene based on new WAG-2 alleles isolated from wheat and its relatives. Only single nucleotide polymorphisms (SNP) and no insertions and deletions (indels) were found in exon sequences of WAG-2 from different species. More SNPs and indels occurred in introns than in exons. For exons, exons+introns and introns, the nucleotide polymorphism π decreased from diploid and tetraploid genotypes to hexaploid genotypes. This finding indicated that the diversity of WAG-2 in diploids was greater than in hexaploids because of the strong selection pressure on the latter. All dn/ds ratios were < 1.0, indicating that WAG-2 belongs to a conserved gene affected by negative selection. Thirty-nine of the 57 particular SNPs and eight of the 10 indels were detected in diploid species. The degree of divergence in intron length among WAG-2 clones and phylogenetic tree topology suggested the existence of three homoeologs in the A, B or D genome of common wheat. Wheat AG-like genes were divided into WAG-1 and WAG-2 clades. The latter clade contained WAG-2, OsMADS3 and ZMM2 genes, indicating functional homoeology among them.