Molecular phylogeny and divergence times of drosophilid species

The phylogenetic relationships and divergence times of 39 drosophilid species were studied by using the coding region of the Adh gene. Four genera--Scaptodrosophila, Zaprionus, Drosophila, and Scaptomyza (from Hawaii)--and three Drosophila subgenera--Drosophila, Engiscaptomyza, and Sophophora--were included. After conducting statistical analyses of the nucleotide sequences of the Adh, Adhr (Adh-related gene), and nuclear rRNA genes and a 905-bp segment of mitochondrial DNA, we used Scaptodrosophila as the outgroup. The phylogenetic tree obtained showed that the first major division of drosophilid species occurs between subgenus Sophophora (genus Drosophila) and the group including subgenera Drosophila and Engiscaptomyza plus the genera Zaprionus and Scaptomyza. Subgenus Sophophora is then divided into D. willistoni and the clade of D. obscura and D. melanogaster species groups. In the other major drosophilid group, Zaprionus first separates from the other species, and then D. immigrans leaves the remaining group of species. This remaining group then splits into the D. repleta group and the Hawaiian drosophilid cluster (Hawaiian Drosophila, Engiscaptomyza, and Scaptomyza). Engiscaptomyza and Scaptomyza are tightly clustered. Each of the D. repleta, D. obscura, and D. melanogaster groups is monophyletic. The splitting of subgenera Drosophila and Sophophora apparently occurred about 40 Mya, whereas the D. repleta group and the Hawaiian drosophilid cluster separated about 32 Mya. By contrast, the splitting of Engiscaptomyza and Scaptomyza occurred only about 11 Mya, suggesting that Scaptomyza experienced a rapid morphological evolution. The D. obscura and D. melanogaster groups apparently diverged about 25 Mya. Many of the D. repleta group species studied here have two functional Adh genes (Adh-1 and Adh-2), and these duplicated genes can be explained by two duplication events.      

Evolution of the Adh locus in the Drosophila willistoni group: the loss of an intron, and shift in codon usage

We report here the DNA sequence of the alcohol dehydrogenase gene (Adh) cloned from Drosophila willistoni. The three major findings are as follows: (1) Relative to all other Adh genes known from Drosophila, D. willistoni Adh has the last intron precisely deleted; PCR directly from total genomic DNA indicates that the deletion exists in all members of the willistoni group but not in any other group, including the closely related saltans group. Otherwise the structure and predicted protein are very similar to those of other species. (2) There is a significant shift in codon usage, especially compared with that in D. melanogaster Adh. The most striking shift is from C to U in the wobble position (both third and first position). Unlike the codon-usage-bias pattern typical of highly biased genes in D. melanogaster, including Adh, D. willistoni has nearly 50% G + C in the third position. (3) The phylogenetic information provided by this new sequence is in agreement with almost all other molecular and morphological data, in placing the obscura group closer to the melanogaster group, with the willistoni group farther distant but still clearly within the subgenus Sophophora.      

Phylogeny of Drosophila and related genera: conflict between molecular and anatomical analyses

Drosophila species are extensively used in biological research; yet, important phylogenetic relationships within the genus and with related genera remain unresolved. The combined data for three genes (Adh, Sod, and Gpdh) statistically resolves outstanding issues. We define the genus Drosophila inclusively so as to include Scaptomyza and Zaprionus (considered distinct genera in the taxonomy of Wheeler, 1981) but excluding Scaptodrosophila. The genus Drosophila so defined is monophyletic. The subgenus Sophophora (including the melanogaster, obscura, and willistoni groups) is monophyletic and the sister clade to all other Drosophila subgenera. The Hawaiian Drosophila (including Scaptomyza) is a monophyletic group, but the subgenus Drosophila is not monophyletic, because the immigrans group is more closely related to the subgenus Hirtodrosophila than to other species of the subgenus Drosophila, such as the virilis and repleta groups.     

Synonymous Substitutions in the Xdh Gene of Drosophila: Heterogeneous Distribution along the Coding Region

The Xdh (rosy) region of Drosophila subobscura has been sequenced and compared to the homologous region of D. pseudoobscura and D. melanogaster. Estimates of the numbers of synonymous substitutions per site (Ks) confirm that Xdh has a high synonymous substitution rate. The distributions of both nonsynonymous and synonymous substitutions along the coding region were found to be heterogeneous. Also, no relationship has been detected between Ks estimates and codon usage bias along the gene, in contrast with the generally observed relationship among genes. This heterogeneous distribution of synonymous substitutions along the Xdh gene, which is expression-level independent, could be explained by a differential selection pressure on synonymous sites along the coding region acting on mRNA secondary structure. The synonymous rate in the Xdh coding region is lower in the D. subobscura than in the D. pseudoobscura lineage, whereas the reverse is true for the Adh gene.     

Differentiation of Muller''s Chromosomal Elements D and E in the Obscura Group of Drosophila

Twenty-two markers located on Muller's elements D or E have been mapped by in situ hybridization in six species of the obscura group of Drosophila and in D. melanogaster. The obscura species can be grouped into a Palearctic cluster (D. subobscura, D. madeirensis and D. guanche) and a Nearctic one (D. pseudoobscura, D. persimilis and D. miranda). Eleven of the probes contain known genes: E74, Acp70A, Est5, hsp28/23, hsp83, emc, hsp70, Xdh, Acph-1, Cec and rp49. The remaining probes are recombinant phages isolated from a D. subobscura genomic library. All these markers hybridize to the putative homologous chromosome or chromosomal arm of elements D and E. Thus, these elements have conserved their genic content during species divergence. Chromosomal homologies proposed previously for each element among the species of the same cluster have been compared with the present results. The distribution of markers within each element has changed considerably as inferred from pairwise comparisons of obscura species included in the two different clusters. Only chromosomal segments defined by closely linked markers have been conserved: one such segment has been detected in element D and three in element E between D. subobscura and D. pseudoobscura.     

Phylogeny of the subgenus sophophora (Diptera: drosophilidae) based on combined analysis of nuclear and mitochondrial sequences

Sequences from the nuclear (nu) alcohol dehydrogenase gene, the nu 28S ribosomal RNA locus, and the mitochondrial cytochrome oxidase II gene were used both individually and in combined analyses to infer the phylogeny of the subgenus Sophophora (Diptera: Drosophilidae). We used several optimality criteria, including maximum likelihood, maximum parsimony, and minimum evolution, to analyze these partitions to test the monophyly of the subgenus Sophophora and its four largest species groups, melanogaster, obscura, saltans, and willistoni. Our results suggest that the melanogaster and obscura species groups are each monophyletic and form a closely related clade. The Neotropical clade, containing the saltans and willistoni species groups, is also recovered, as previous studies have suggested. While the saltans species group is strongly supported as monophyletic, the results of several analyses indicate that the willistoni species group may be paraphyletic with respect to the saltans species group.     

The Evolution of Duplicate Glyceraldehyde-3-Phosphate Dehydrogenase Genes in Drosophila

In Drosophila melanogaster there are two genes which encode the enzyme glyceraldehyde-3-phosphate dehydrogenase (GAPDH), Gapdh-43E and Gapdh-13F. We have shown that Gapdh-43E codes for the GAPDH subunit with an apparently larger molecular weight while Gapdh-13F encodes the GAPDH subunit having an apparently smaller molecular weight. Immunoblots of sodium dodecyl sulfate gels were used to survey species from throughout the genus and results indicated that two classes of GAPDH subunits are present only in Drosophila species of the melanogaster and takahashi subgroups of the melanogaster group. Only the smaller subunit is found in species of the obscura group while all other species have only a large subunit. Drosophila hydei was analyzed at the DNA level as a representative species of the subgenus Drosophila. The genome of this species has a single Gapdh gene which is localized at a cytogenetic position likely to be homologous to Gapdh-43 E of D. melanogaster. Comparison of its sequence with the sequence of the D. melanogaster Gapdh genes indicates that the two genes of D. melanogaster are more similar to one another than either is to the gene from D. hydei. The Gapdh gene from D. hydei contains an intron following codon 29. Neither Gapdh gene of D. melanogaster has an intron within the coding region. Southern blots of genomic DNA were used to determine which species have duplicate Gapdh genomic sequences. Gene amplification was used to determine which species have a Gapdh gene that is interrupted by an intron. Species of the subgenus Drosophila have a single Gapdh gene with an intron. Species of the willistoni and saltans groups have a single Gapdh gene that does not contain an intron.(ABSTRACT TRUNCATED AT 250 WORDS)     

Nucleotide sequence of the Xdh region in Drosophila pseudoobscura and an analysis of the evolution of synonymous codons

The nucleotide sequence of the Xdh region of Drosophila pseudoobscura is presented. The Xdh gene structure and organization are compared with the homologous region in D. melanogaster. This locus is shown to have similar organization in the two species, although an additional intron and three insertion/deletion events are described for the D. pseudoobscura coding region. The encoded proteins are predicted to have very similar charges and hydrophobic/hydrophilic domains even though 11% of the amino acids are different. A gene 5' to Xdh, putative l(3)s12, is suggested from sequence similarity between the species. Synonymous differences at the Xdh locus between the two species are analyzed using a new method described in the preceding paper by Lewontin. This analysis shows that synonymous positions within the Xdh locus are evolving at very different rates, being dependent on level of codon redundancy. A comparison of synonymous divergence between D. melanogaster and D. pseudoobscura in five additional genes reveals variation in the level of synonymous substitution.      

On the detection of nonrandom associations between DNA polymorphisms in natural populations of Drosophila

The capacity to detect nonrandom associations between restriction-map variants was examined in eight gene regions of Drosophila melanogaster (yellow-achaetescute, white, Zw, Adh, Est6, and rosy) and D. pseudoobscura (Adh and Xdh), on the basis of published population data. The statistical power from individual pairwise tests was both heterogeneous and generally low across gene regions. Sample sizes larger than those currently being used are needed to ensure any power to detect disequilibrium by individual tests. It is found that the heterogeneity in power is mostly explained by large differences in the intensity of sample disequilibrium among regions. The yellow-achaete- scute, Zw, and Adh loci of D. melanogaster displayed both the highest mean power (approximately 0.4) and a very great disequilibrium (mean absolute values of D' were 0.8-1). By contrast, all the other gene regions exhibited lower mean power (approximately 0.2) and moderate levels of disequilibrium (0.4-0.6). Although the proportion of significant pairwise associations, especially for white, Est6, and rosy in D. melanogaster and for Adh and Xdh in D. pseudoobscura, is more or less close to the type I error, simultaneous-inference significance tests show that gametic disequilibrium is occurring at the eight DNA regions examined.      

Changes in the recombinational environment affect divergence in the yellow gene of Drosophila

The complete coding region of the yellow (y) gene was sequenced in different Drosophila species. In the species of the melanogaster subgroup (D. melanogaster, D. simulans, D. mauritiana, D. yakuba, and D. erecta), this gene is located at the tip of the X chromosome in a region with a strong reduction in recombination rate. In contrast, in D. ananassae (included in the ananassae subgroup of the melanogaster group) and in the obscura group species (D. subobscura, D. madeirensis, D. guanche, and D. pseudoobscura), the y gene is located in regions with normal recombination rates. As predicted by the hitchhiking and background selection models, this change in the recombinational environment affected synonymous divergence in the y-gene-coding region. Estimates of the number of synonymous substitutions per site were much lower between the obscura group species and D. ananassae than between the species of the obscura group and the melanogaster subgroup. In fact, a highly significant increase in the rate of synonymous substitution was detected in all lineages leading to the species of the melanogaster subgroup relative to the D. ananassae lineage. This increase can be explained by a higher fixation rate of mutations from preferred to unpreferred codons (slightly deleterious mutations). The lower codon bias detected in all species of the melanogaster subgroup relative to D. ananassae (or to the obscura group species) would be consistent with this proposal. Therefore, at least in Drosophila, changes in the recombination rate in different lineages might cause deviations of the molecular-clock hypothesis and contribute to the overdispersion of the rate of synonymous substitution. In contrast, the change in the recombinational environment of the y gene has no detectable effect on the rate of amino acid replacement in the Yellow protein.     

Phylogenetic relationships between Sophophora and Lordiphosa, with proposition of a hypothesis on the vicariant divergences of tropical lineages between the Old and New Worlds in the family Drosophilidae

Despite many studies on the phylogeny of the subgenus Sophophora, its monophyly has not been established, especially in relation to its putative relative, the genus Lordiphosa. We analyzed their phylogenetic relationships using DNA sequence data of two mitochondrial genes (ND2 and COII) and two nuclear genes (Adh and 28SrRNA). In constructing phylogenetic trees, we accounted for the problem of among-taxa nucleotide compositional heterogeneity, and took a sequence-partitioning approach to allow multiple substitution models for nucleotide sequences that have evolved under different evolutionary processes, particularly developing a novel, sequence-partitioning procedure for Neighbor Joining (NJ) tree construction. Trees constructed by different methods showed an almost identical and strongly supported topology in which Sophophora was paraphyletic: Lordiphosa was placed as the sister to the Neotropical Sophophora consisting of the saltans and willistoni groups, and Sophophora was divided into the clade of Lordiphosa+Neotropical Sophophora and the clade of the obscura+melanogaster groups. Based on the estimated time, 45.9 Mya, of divergence between the Old World Lordiphosa and the Neotropical Sophophora and evidence from paleontology, paleo-geography and -climatology, we propose a hypothesis that this vicariant divergence should have occurred when the North Atlantic Land Bridge between Europe and North America broke in the middle Eocene Epoch.     

Local recombination and mutation effects on molecular evolution in Drosophila.

I studied the cause of the significant difference in the synonymous-substitution pattern found in the achaete-scute complex genes in two Drosophila lineages, higher codon bias in Drosophila yakuba, and lower bias in D. melanogaster. Besides these genes, the functionally unrelated yellow gene showed the same substitution pattern, suggesting a region-dependent phenomenon in the X-chromosome telomere. Because the numbers of A/T --> G/C substitutions were not significantly different from those of G/C --> A/T in the yellow noncoding regions of these species, a AT/GC mutational bias could not completely account for the synonymous-substitution biases. In contrast, we did find an approximately 14-fold difference in recombination rates in the X-chromosome telomere regions between the two species, suggesting that the reduction of recombination rates in this region resulted in the reduction of the efficacy of selection in D. melanogaster. In addition, the D. orena yellow showed a 5% increase in the G + C content at silent sites in the coding and noncoding regions since the divergence from D. erecta. This pattern was significantly different from those at the orena Adh and Amy loci. These results suggest that local changes in recombination rates and mutational pressures are contributing to the irregular synonymous-substitution patterns in Drosophila.     

Evolution of gypsy endogenous retrovirus in the Drosophila obscura species group

The Ty3/gypsy family of retroelements is closely related to retroviruses, and some of their members have an open reading frame resembling the retroviral gene env. Sequences homologous to the gypsy element from Drosophila melanogaster are widely distributed among Drosophila species. In this work, we report a phylogenetic study based mainly on the analysis of the 5' region of the env gene from several species of the obscura group, and also from sequences already reported of D. melanogaster, Drosophila virilis, and Drosophila hydei. Our results indicate that the gypsy elements from species of the obscura group constitute a monophyletic group which has strongly diverged from the prototypic D. melanogaster gypsy element. Phylogenetic relationships between gypsy sequences from the obscura group are consistent with those of their hosts, indicating vertical transmission. However, D. hydei and D. virilis gypsy sequences are closely related to those of the affinis subgroup, which could be indicative of horizontal transmission.