Vertical Transmission of the Retrotransposable Elements R1 and R2 during the Evolution of the Drosophila Melanogaster Species Subgroup

R1 and R2 are non-long-terminal repeat retrotransposable elements that insert into specific sequences of insect 28S ribosomal RNA genes. These elements have been extensively described in Drosophila melanogaster. To determine whether these elements have been horizontally or vertically transmitted, we characterized R1 and R2 elements from the seven other members of the melanogaster species subgroup by genomic blotting and nucleotide sequencing. Each species was found to have homogeneous families of R1 and R2 elements with the exception of erecta and orena, which have no R2 elements. The DNA sequences of multiple R1 and R2 copies from each species indicated nucleotide divergence within each species averaged only 0.48% for R1 and 0.35% for R2, well below the level of divergence among the species. Most copies of R1 and R2 (40 of 47) sequenced from the seven species were potentially functional, as indicated by the absence of premature termination codons or translational frameshifts that would destroy the open reading frame of the element. The sequence relationships of both the R1 and R2 elements from the various members of the melanogaster subgroup closely followed that of the species phylogeny, suggesting that R1 and R2 have been stably maintained by vertical transmission since the origin of this species subgroup 17-20 million years ago. The remarkable stability of R1 and R2, compared to what has been suggested for transposable elements that insert at multiple locations in these same species, may be due to their unique specificity for sites in the rRNA gene locus. Under low copy number conditions, when it is essential for any mobile element to transpose, the insertion specificities of R1 and R2 ensure uniform developmentally regulated target sites that can be occupied with little or no detrimental effect on the host.     

Evolution of R1 and R2 in the rDNA units of the genus Drosophila

R1 and R2 are non-long terminal repeat (non-LTR) retrotransposable elements that specifically insert in the 28S ribosomal RNA (rRNA) genes of insects. Using the Drosophila genus, which includes some of the best characterized insect taxa, we have conducted a number of studies on the evolution of these elements. We find that R1 and R2 are subject to the same recombinational forces that give rise to the concerted evolution of the rDNA units. The turnover of R1 and R2 elements can be readily documented in different strains of D. melanogaster using 5′ truncated elements as restriction-length polymorphisms. This turnover leads to uniform populations of elements with nucleotide sequence divergence of different copies averaging only 0.23% for the R2 and 0.47% for the R1 elements. Molecular phylogenetic analysis of elements from 16 different species of Drosophila suggests that these elements have been stable components of the rDNA locus for the 50–70 million year history of the Drosophila genus. Using changes at synonymous positions within the protein-encoding regions as estimates of the baseline substitution rate, it could be shown that R1 and R2 are evolving at rates similar to that of typical protein encoding genes provided corrections are made for the low codon bias of the elements. R1 and R2 are clearly well-adapted for their existence in the rDNA units of their host. This revised version was published online in August 2006 with corrections to the Cover Date.     

Evolutionary stability of the R1 retrotransposable element in the genus Drosophila

R1 is a non-long terminal repeat (non-LTR) retrotransposable element that inserts into a specific sequence of insect 28S ribosomal RNA genes. We have previously shown that this element has been maintained through vertical transmission in the melanogaster species subgroup of Drosophila. To address whether R1 elements have been vertically transmitted for longer periods of evolutionary time, the analysis has been extended to 11 other species from four species groups of the genus Drosophila (melanogaster, obscura, testecea, and repleta). All sequenced elements appeared functional on the basis of the preservation of their open-reading frames and consistently higher rate of substitution at synonymous sites relative to replacement sites. The phylogenetic relationships of the R1 elements from all species analyzed were congruent with the species phylogenies, suggesting that the R1 elements have been vertically transmitted since the inception of the Drosophila genus, an estimated 50-70 Mya. The stable maintenance of R1 through the germ line appears to be the major mechanism for the widespread distribution of these elements in Drosophila. In two species, D. neotestecea of the testecea group and D. takahashii of the melanogaster group, a second family of R1 elements was also present that differed in sequence by 46% and 31%, respectively, from the family that was congruent with the species phylogeny. These second families may represent occasional horizontal transfers or, alternatively, they could reflect the ability of R1 elements to diverge into new families within a species and evolve independently.      

Synonymous substitution rates in Drosophila: Mitochondrial versus nuclear genes

Synonymous substitution rates in mitochondrial and nuclear genes of Drosophila were compared. To make accurate comparisons, we considered the following: (1) relative synonymous rates, which do not require divergence time estimates, should be used; (2) methods estimating divergence should take into account base composition; (3) only very closely related species should be used to avoid effects of saturation; (4) the heterogeneity of rates should be examined. We modified the methods estimating synonymous substitution numbers to account for base composition bias. By using these methods, we found that mitochondrial genes have 1.7–3.4 times higher synonymous substitution rates than the fastest nuclear genes or 4.5–9.0 times higher rates than the average nuclear genes. The average rate of synonymous transversions was 2.7 (estimated from the melanogaster species subgroup) or 2.9 (estimated from the obscura group) times higher in mitochondrial genes than in nuclear genes. Synonymous transversions in mitochondrial genes occurred at an approximately equivalent rate to those in the fastest nuclear genes. This last result is not consistent with the hypothesis that the difference in turnover rates between mitochondrial and nuclear genomes is the major factor determining higher synonymous substitution rates in mtDNA. We conclude that the difference in synonymous substitution rates is due to a combination of two factors: a higher transitional mutation rate in mtDNA and constraints on nuclear genes due to selection for codon usage. Received: 27 November 1996 / Accepted: 8 May 1997     

Nucleotide variation of seven genes in Drosophila kikkawai

We examined levels and patterns of nucleotide variation in 21 strains of Drosophila kikkawai from Miyako island, Japan for the partial regions of the following seven nuclear genes: Adh, Ddc, esc, ksr, Pgi, su(f), and Tpi. The nucleotide variation at total sites (pi(t)) ranged from 0.0013 in the ksr, to 0.0173 in the Adh. The nucleotide divergence at total sites (K(t)) between D. kikkawai and D. lini ranged from 0.0286 in the Tpi to 0.0687 in the su(f). The levels of nucleotide polymorphism and divergence were heterogeneous among the investigated gene regions. The HKA test, which tests imbalance between the intra and interspecific nucleotide variation, showed that the intraspecific nucleotide variation in the Pgi region was much lower than the interspecific variation, while intraspecific variation in the Tpi region was only slightly lower than interspecific variation. The MK test showed an excess of low frequency replacement polymorphic changes in the Adh region, suggesting that most replacement mutations are deleterious. Fay and Wu's test detected an excess of newly arisen variants in the Ddc region. In total, four of the seven gene regions showed significant deviation from the neutrality.     

Variable Rates of Evolution among Drosophila Opsin Genes

DNA sequences and chromosomal locations of four Drosophila pseudoobscura opsin genes were compared with those from Drosophila melanogaster, to determine factors that influence the evolution of multigene families. Although the opsin proteins perform the same primary functions, the comparisons reveal a wide range of evolutionary rates. Amino acid identities for the opsins range from 90% for Rh2 to more than 95% for Rh1 and Rh4. Variation in the rate of synonymous site substitution is especially striking: the major opsin, encoded by the Rh1 locus, differs at only 26.1% of synonymous sites between D. pseudoobscura and D. melanogaster, while the other opsin loci differ by as much as 39.2% at synonymous sites. Rh3 and Rh4 have similar levels of synonymous nucleotide substitution but significantly different amounts of amino acid replacement. This decoupling of nucleotide substitution and amino acid replacement suggests that different selective pressures are acting on these similar genes. There is significant heterogeneity in base composition and codon usage bias among the opsin genes in both species, but there are no consistent relationships between these factors and the rate of evolution of the opsins. In addition to exhibiting variation in evolutionary rates, the opsin loci in these species reveal rearrangements of chromosome elements.     

A single lineage of r2 retrotransposable elements is an active, evolutionarily stable component of the Drosophila rDNA locus

R2 elements are non-long-terminal-repeat (non-LTR) retrotransposons that insert specifically in the 28S rRNA genes of many insects. Previous reports concerning this element in the genus Drosophila have suggested that R2 elements are absent from many species of this genus, particularly those species from the subgenus Drosophila. In this report, we present an extensive study of the distribution and evolution of R2 elements in Drosophila. A PCR survey of 59 species from 23 species groups of the two major Drosophila subgenera found that R2 elements are present in all but two species of the melanogaster species subgroup. Phylogenetic analysis based on partial nucleotide sequences of R2 elements from 23 species demonstrates that the relationships of R2 elements are congruent with those of the Drosophila species phylogeny, suggesting that these elements have been vertically inherited since the divergence of this genus some 60 MYA. Sequence variation between different copies of R2 elements within each species was less than 0.16%, indicating that these elements are undergoing concerted evolution similar to that of the 28S genes. Several properties of the R2 sequences suggest that these elements depend on retrotransposition in addition to simple recombination to remain within the rDNA locus: the rates of synonymous substitutions averaged 4.8 times the rate of replacement substitutions, 82 of 83 R2 copies partially sequenced contained intact open reading frames, and, finally, length variation associated with the poly(A) 3' tails indicated that many R2 copies are the direct result of retrotransposition.      

The molecular evolution of the alcohol dehydrogenase and alcohol dehydrogenase-related genes in the Drosophila melanogaster species subgroup

The DNA sequences of the Adh genes of three members of the Drosophila melanogaster species subgroup have been determined. This completes the Adh sequences of the eight species of this subgroup. Two species, D. yakuba and D. teissieri, possess processed Adh pseudogenes. In all of the species of the subgroup, a gene of unknown function, Adhr, is located about 300 bp 3' to Adh. Although this gene is experiencing a higher rate of synonymous substitution than Adh, it is more constrained at the amino acid level. Phylogenetic relationships between all eight members of the melanogaster subgroup have been analyzed using a variety of methods. All analyses suggested that the D. yakuba and D. teissieri pseudogenes have a single common ancestor, rather than evolving independently in each species, and that D. melanogaster is the sister species to D. simulans, D. sechellia, and D. mauritiana. The evolutionary relationships of the latter three species remain equivocal.      

Low rates of silent substitution in nuclear genes of two distantly related Scrophulariaceae (Antirrhinum and Verbascum)

Low levels of genetic diversity and divergence at nuclear loci have previously been observed for cycloidea and fil1-like genes within and between several Antirrhinum species, and divergence at these loci is also low between species in genera at different levels of relatedness in the former family Scrophulariaceae (Digitalis and Verbascum). The low divergence values are surprising, because (based on the sequences of chloroplast loci) the Scrophulariaceae are thought to be polyphyletic, with two anciently diverged clades, and the species we compared belonged to the two different clades. Here, we extend our studies of sequence divergence to more nuclear genes: fil2, far, globosa, and ADH: Detailed studies revealed that in Antirrhinum these genes belong to gene families. Low levels of divergence between Antirrhinum and Verbascum were observed for four of the loci studied, fil2-1, fil2-2, far-L, and globosa, similar to our previous observations. We discuss hypotheses to explain these low synonymous divergence values. For Adh, no cases of very similar sequences were found, but, rather, our sequences from the three different genera (Antirrhinum, Digitalis, and Verbascum) were all very diverged. Repeated gene duplication and loss of elements in the Adh gene family is likely in these lineages, making it impossible to determine orthology of the Adh genes.     

Nucleotide Substitution Rate of Mammalian Mitochondrial Genomes

We present here for the first time a comprehensive study based on the analysis of closely related organisms to provide an accurate determination of the nucleotide substitution rate in mammalian mitochondrial genomes. This study examines the evolutionary pattern of the different functional mtDNA regions as accurately as possible on the grounds of available data, revealing some important ``genomic laws.' The main conclusions can be summarized as follows. (1) High intragenomic variability in the evolutionary dynamic of mtDNA was found. The substitution rate is strongly dependent on the region considered, and slow- and fast-evolving regions can be identified. Nonsynonymous sites, the D-loop central domain, and tRNA and rRNA genes evolve much more slowly than synonymous sites and the two peripheral D-loop region domains. The synonymous rate is fairly uniform over the genome, whereas the rate of nonsynonymous sites depends on functional constraints and therefore differs considerably between genes. (2) The commonly accepted statement that mtDNA evolves more rapidly than nuclear DNA is valid only for some regions, thus it should be referred to specific mitochondrial components. In particular, nonsynonymous sites show comparable rates in mitochondrial and nuclear genes; synonymous sites and small rRNA evolve about 20 times more rapidly and tRNAs about 100 times more rapidly in mitochondria than in their nuclear counterpart. (3) A species-specific evolution is particularly evident in the D-loop region. As the divergence times of the organism pairs under consideration are known with sufficient accuracy, absolute nucleotide substitution rates are also provided. Received: 11 May 1998 / Accepted: 2 September 1998     

Rates of DNA change and phylogeny from the DNA sequences of the alcohol dehydrogenase gene for five closely related species of Hawaiian Drosophila

The sequence of 1.6 kb of DNA surrounding the alcohol dehydrogenase (Adh) gene from five species of the Planitibia subgroup of the Hawaiian picture-winged Drosophila, with estimated divergence times of 0.4-5.1 Myr, has been determined. The gene trees which were found by using the sequence divergence from different regions of the sequences are generally in accord with the phylogeny proposed for these species when chromosomal inversions and island of origin are used. One of the species (D. picticornis) appears to be more distant from the other species in this group than they are from a member of the Grimshawi group (D. affinidisjuncta) which is chromosomally more distant. Two of the species (D. differens and D. plantibia) show heterogeneity in the nucleotide changes in the Adh coding region, heterogeneity which is interpreted to be due to a gene conversion or recombination after hybridization between the two species. The minimal rate of nucleotide substitution of synonymous nucleotides and of nontranscribed nucleotides downstream from the coding region is estimated as 1.5 x 10(-8) and 1.1 x 10(-8) substitutions/nucleotide/year, respectively. This rate is two to three times the maximal rate estimated for mammalian synonymous substitutions.     

Rates of Synonymous Substitution and Base Composition of Nuclear Genes in Drosophila

We compared the rates of synonymous (silent) substitution among various genes in a number of species of Drosophila. First, we found that even for a particular gene, the rate of synonymous substitution varied considerably with Drosophila lineages. Second, we showed a large variation in synonymous substitution rates among nuclear genes in Drosophila. These rates of synonymous substitution were correlated negatively with C content and positively with A content at the third codon positions. Nucleotide sequences were also compared between pseudogenes and their functional homologs. The C content of the pseudogenes was lower than that of the functional genes and the A content of the former was higher than that of the latter. Because the synonymous substitution for functional genes and the nucleotide substitution for pseudogenes are exempted from any selective constraint at the protein level, these observations could be explained by a biased pattern of mutation in the Drosophila nuclear genome. Such a bias in the mutation pattern may affect the molecular clock (local clock) of each nuclear gene of each species. Finally, we obtained the average rates of synonymous substitution for three gene groups in Drosophila; 11.0 x 10(-9), 17.5 x 10(-9) and 27.1 x 10(-9)/site/year.     

Codon Usage Bias and Base Composition of Nuclear Genes in Drosophila

The nuclear genes of Drosophila evolve at various rates. This variation seems to correlate with codon-usage bias. In order to elucidate the determining factors of the various evolutionary rates and codon-usage bias in the Drosophila nuclear genome, we compared patterns of codon-usage bias with base compositions of exons and introns. Our results clearly show the existence of selective constraints at the translational level for synonymous (silent) sites and, on the other hand, the neutrality or near neutrality of long stretches of nucleotide sequence within noncoding regions. These features were found for comparisons among nuclear genes in a particular species (Drosophila melanogaster, Drosophila pseudoobscura and Drosophila virilis) as well as in a particular gene (alcohol dehydrogenase) among different species in the genus Drosophila. The patterns of evolution of synonymous sites in Drosophila are more similar to those in the prokaryotes than they are to those in mammals. If a difference in the level of expression of each gene is a main reason for the difference in the degree of selective constraint, the evolution of synonymous sites of Drosophila genes would be sensitive to the level of expression among genes and would change as the level of expression becomes altered in different species. Our analysis verifies these predictions and also identifies additional selective constraints at the translational level in Drosophila.     

Rates of R1 and R2 retrotransposition and elimination from the rDNA locus of Drosophila melanogaster

R1 and R2 elements are non-LTR retrotransposons that insert specifically into the 28S rRNA genes of arthropods. The process of concerted evolution of the rDNA locus should give rise to rapid turnover of these mobile elements compared to elements that insert at sites throughout a genome. To estimate the rate of R1 and R2 turnover we have examined the insertion of new elements and elimination of old elements in the Harwich mutation accumulation lines of Drosophila melanogaster, a set of inbred lines maintained for >350 generations. Nearly 300 new insertion and elimination events were observed in the 19 Harwich lines. The retrotransposition rate for R1 was 18 times higher than the retrotransposition rate for R2. Both rates were within the range previously found for retrotransposons that insert outside the rDNA loci in D. melanogaster. The elimination rates of R1 and R2 from the rDNA locus were similar to each other but over two orders of magnitude higher than that found for other retrotransposons. The high rates of R1 and R2 elimination from the rDNA locus confirm that these elements must maintain relatively high rates of retrotransposition to ensure their continued presence in this locus.     

Mitochondrial genome sequence evolution in Chlamydomonas

The mitochondrial genomes of the Chlorophyta exhibit significant diversity with respect to gene content and genome compactness; however, quantitative data on the rates of nucleotide substitution in mitochondrial DNA, which might help explain the origin of this diversity, are lacking. To gain insight into the evolutionary forces responsible for mitochondrial genome diversification, we sequenced to near completion the mitochondrial genome of the chlorophyte Chlamydomonas incerta, estimated the evolutionary divergence between Chlamydomonas reinhardtii and C. incerta mitochondrial protein-coding genes and rRNA-coding regions, and compared the relative evolutionary rates in mitochondrial and nuclear genes. Synonymous and nonsynonymous substitution rates do not differ significantly between the mitochondrial and nuclear protein-coding genes. The mitochondrial rRNA-coding regions, however, are evolving much faster than their nuclear counterparts, and this difference might be explained by relaxed functional constraints on the mitochondrial translational apparatus due to the small number of proteins synthesized in Chlamydomonas mitochondria. Substitution rates at synonymous sites in a nonstandard mitochondrial gene (rtl) and at intronic and synonymous sites in nuclear genes expressed at low levels suggest that the mutation rate is similar in these two genetic compartments. Potential evolutionary forces shaping mitochondrial genome evolution in Chlamydomonas are discussed.     

Nucleotide polymorphism in the Adh2 region of the wild rice Oryza rufipogon

DNA variation in the alcohol dehydrogenase (Adh2) region of the wild rice Oryza rufipogon and its related species was analyzed to clarify maintenance mechanisms of the DNA variation in these species. A dimorphic pattern was detected in the Adh2 region of O. rufipogon. The silent nucleotide diversity (π) in the Adh2 region in O. rufipogon was 0.011, which was higher than that of the Adh1 region in O. rufipogon. Especially, a high nucleotide diversity was detected at synonymous sites of the catalytic domain 1. Average nucleotide diversity at silent sites within each of the dimorphic sequence types of the Adh2 region was similar to that in the Adh1 region, indicating that the high level of silent polymorphism in the Adh2 region was caused by the difference between the dimorphic sequence types. On the other hand, the level of replacement polymorphism in the Adh2 region was as low as that in the Adh1 region. The neutrality test of Fu and Li indicated significantly negative deviation from the neutral mutation model for the replacement sites of the Adh2 region. This result suggests purifying selection on the replacement sites of the Adh2 region, as detected for the Adh1 region. Significant linkage disequilibria (16.4% of the tests) were detected between the Adh1 and Adh2 regions. Even when nonrandom association was tested for the strains belonging to one of the divergent sequence types of the Adh2 region, significant interlocus linkage disequilibria were detected. The close physical distance and/or epistasis between the two Adh regions could be invoked to explain these nonrandom associations.     

Contrasting selection modes at the Adh1 locus in outcrossing Miscanthus sinensis vs. inbreeding Miscanthus condensatus (Poaceae)

We estimated DNA sequence variation of the Adh1 locus in the outcrossing Miscanthus sinensis (Poaceae) and its close selfing relative, M. condensatus. Tajima's test of selection is significantly negative for both overall exons and replacement sites in M. sinensis. Among its entire sample, nucleotide diversity of nonsynonymous sites is higher than that of synonymous sites. A McDonald and Kreitman test of neutrality indicates an excess of intraspecific replacement polymorphisms, suggesting possible directional selection toward advantageous mutants. However, frequent intragenic recombination suggests both purifying and positive selection is unlikely. Recent demographic expansions coupled with relaxation of purifying selection may have resulted in elevated genetic diversity at the Adh1 locus as well as the trnL-trnF intergenic spacer of cpDNA in this outcrossing species. In contrast, low levels of genetic diversity were detected at both the Adh1 locus and the cpDNA spacer in M. condensatus, consistent with bottlenecks associated with selfing in all populations. While Tajima's D and Fu and Li's F statistics did not reveal deviation from neutrality at the Adh1 locus in M. condensatus, 12 replacements vs. 10 synonymous changes were detected. Based on pairwise comparisons of the d(N)/d(S) ratio, lineages of closely related populations of the species distributed along saline habitats appeared to be under directional selection.     

Molecular Population Genetics of an Electrophoretically Monomorphic Protein in the Alcohol Dehydrogenase Region of Drosophila Pseudoobscura

Nucleotide sequence data from the alcohol dehydrogenase (Adh) region of 18 isochromosomal strains of Drosophila pseudoobscura were used to determine whether the lack of amino acid polymorphism in ADH results from a low neutral mutation rate or a recent directional selection event. We estimated the neutral mutation parameter, 4Nmu, in synonymous sites for 17 subregions of Adh. The nucleotide diversity data were tested for departures from an equilibrium neutral model with two statistical tests. The Tajima test and the Hudson, Kreitman and Aguade test each failed to reject a neutral model. These results suggest that the ADH enzyme of D. pseudoobscura lacks amino acid polymorphisms because the neutral mutation rate of nonsynonymous sites is low. The neutral mutation parameter for synonymous sites is heterogeneous between domains of the Adh region. These data indicate that selective constrains on synonymous sites can vary between functional domains.