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.     

The phylogenetic relationships of flies in the family Drosophilidae deduced from mtDNA sequences.

The phylogenetic relationships of several taxa from representative genera, subgenera, groups, and subgroups in the Drosophilidae were examined using sequences from a 905-bp mtDNA fragment. Conventional cloning and sequencing techniques were used to obtain nucleotide sequences. In addition, polymerase chain reaction primers were designed for the rapid amplification and sequencing of this region for the species examined in the Drosophilidae. Phylogenetic analysis was done by cladistic techniques. Because of the coding nature of the 905-bp mtDNA fragment, several separate analyses of these sequences were performed. The genera Scaptomyza and Hirtodrosophila occupy ancestral branching positions in the molecular phylogeny. The genera Chymomyza and Zaprionus have intermediate branching positions, while the subgenera Drosophila and Sophophora are in the most derived position in the molecular phylogeny. Within the subgenus Sophophora, there is little resolution using these sequences, while within the subgenus Drosophila, D. melanica, D. funebris, and D. pinicola form a clade in a derived part of the phylogeny, with D. robusta and D. immigrans branching in an intermediate position in the phylogeny. D. mercatorum, a member of the repleta species group, occupies an ancestral position in the molecular phylogeny.     

A phylogeny of Drosophilidae using the Amyrel gene: questioning the Drosophila melanogaster species group boundaries

In this study, the phylogenetic relationships of 164 species of the family Drosophilidae are discussed, using the Amyrel gene, a member of the α -amylase multigene family. This study focuses on numerous species groups in the subgenera Sophophora and Drosophila of the genus Drosophila but also includes other closely related genera. Nucleotide data were analysed by several methods: maximum parsimony, neighbour joining, maximum likelihood and Bayesian inference. Heterogeneity of base composition (mainly low GC contents in the species groups willistoni and saltans ) has been addressed. In all analyses, the genus Drosophila appeared paraphyletic. The subgenus Sophophora clearly appeared to be a monophyletic group, showing well-resolved clades, with the Neotropical groups arising in a basal position. Here, it is proposed to raise the species subgroups ananassae and montium to the rank of species group, and to restrict the melanogaster species group to the melanogaster subgroup plus the 'Oriental' subgroups, among which the suzukii subgroup is polyphyletic. Some related genera such as Zaprionus , Liodrosophila , Scaptomyza and Hirtodrosophila are clustered with, or inside the subgenus Drosophila , which is therefore paraphyletic and should be reviewed.     

Switch in codon bias and increased rates of amino acid substitution in the Drosophila saltans species group.

We investigated the nucleotide composition of five genes, Xdh, Adh, Sod, Per, and 28SrRNA, in nine species of Drosophila (subgenus Sophophora) and one of Scaptodrosophila. The six species of the Drosophila saltans group markedly differ from the others in GC content and codon use bias. The GC content in the third codon position, and to a lesser extent in the first position and the introns, is higher in the D. melanogaster and D. obscura groups than in the D. saltans group (in Scaptodrosophila it is intermediate but closer to the melanogaster and obscura species). Differences are greater for Xdh than for Adh, Sod, Per, and 28SrRNA, which are functionally more constrained. We infer that rapid evolution of GC content in the saltans lineage is largely due to a shift in mutation pressure, which may have been associated with diminished natural selection due to smaller effective population numbers rather than reduced recombination rates. The rate of GC content evolution impacts the rate of protein evolution and may distort phylogenetic inferences. Previous observations suggesting that GC content evolution is very limited in Drosophila may have been distorted due to the restricted number of genes and species (mostly D. melanogaster) investigated.     

Molecular evolution of Drosophila odorant receptor genes

A total of 752 odorant receptor (Or) genes, including pseudogenes, were identified in 11 Drosophila species and named after their orthologs in Drosophila melanogaster. The 813 Or genes, including 61 from D. melanogaster, were classified into 59 orthologous groups that are well supported by gene phylogeny. By reconciling with the gene family phylogeny, we estimated the number of gene duplication/loss events and intron gain/loss events in the species phylogeny. We found that these events are particularly frequent in Drosophila grimshawi, Drosophila willistoni, and obscura group. More than half of the duplicated genes stay as tandem arrays, whose size range from 2 to 8. These genes vary in sequence and some likely underwent positive selection, indicating that the gene duplication was important for flies to acquire new olfactory functions. We hypothesize that Or genes conferred the basic olfactory repertoire to ancestral flies before the speciation of the Drosophila and Sophophora subgenera about 40 Mya. This repertoire has been largely maintained in the current species, whereas lineage-specific gene duplication seems to have led to additional specialization in some species in response to specific ecological conditions.     

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 the Drosophila immigrans species group (Diptera: Drosophilidae) based on Adh and Gpdh sequences

The immigrans species group in the Drosophilinae is one of the representative species groups of Drosophila in East Asia. Although this group constitutes a significant part of the drosophilid fauna in the Old World, only a few species have been analyzed in previous molecular phylogenetic studies. To study the phylogeny of the immigrans group, we analyzed the nucleotide sequences of two nuclear genes, alcohol dehydrogenase (Adh) and glycerol-3-phosphate dehydrogenase (Gpdh), for 36 drosophilid species, including 12 species of the immigrans group. In the resultant phylogenetic trees, 10 species of the immigrans group (D. immigrans, D. formosana, D. ruberrima, D. albomicans, D. nasuta, D. neonasuta, D. pallidifrons, D. hypocausta, D. neohypocausta, D. siamana) consistently formed a clade (the immigrans group proper), although the phylogeny within this clade did not exactly correspond to the classification of species subgroups. However, D. annulipes and D. quadrilineata, both of which belong to the quadrilineata subgroup of the immigrans group, were not included in the immigrans group proper. Furthermore, we obtained the unexpected result that D. annulipes was included in a clade comprising Scaptomyza and Hawaiian Drosophila, together with D. maculinotata of the funebris group, although the phylogenetic relationships within this clade remain uncertain and need to be substantiated with further studies. Thus, according to the present study, the immigrans group is polyphyletic.     

The structure of the Adh locus of Drosophila mettleri: an intermediate in the evolution of the Adh locus in the repleta group of Drosophila.

Members of species of the mulleri and hydei subgroups of the repleta group of Drosophila have duplicate Adh genes. The Adh regions of D. mojavensis, D. mulleri, and D. hydei contain three genes--a pseudogene, Adh-2, and Adh-1--arranged 5' to 3'. To understand the evolution of the triplicate Adh structure, we have cloned and sequenced the Adh locus of D. mettleri. This region consists of a 5' pseudogene and a 3' functional Adh gene. On the basis of the structure and nucleotide sequence comparisons of Adh genes of D. mettleri and other species, we propose that an initial duplication of the ancestral Adh gene generated two Adh genes arranged in tandem. The more 5' Adh gene became a pseudogene, while the more 3' gene remained functional through all the developmental stages. A second duplication of this 3' gene resulted in Adh regions with three genes--a pseudogene, Adh-2, and Adh-1.     

Distribution and conservation of sequences homologous to the 1731 retrotransposon in Drosophila

The distribution of 1731 retrotransposon-hybridizing sequences in the family Drosophilidae has been studied using a 1731 probe from Drosophila melanogaster. Squash blot and Southern blot analyses of 42 species reveal that the 1731 sequences are widespread within both the Sophophora and Drosophila subgenera and are also present in the genera Scaptomyza and Zaprionus. Hence the 1731 retrotransposon family appears to have a long evolutionary history in the Drosophilidae genome. Differences of hybridization signal intensity suggested that the 1731 sequence is well conserved only in the three species most closely related to D. melanogaster (D. simulans, D. mauritiana, and D. sechellia). A survey of insertion sites in numerous different populations of the previous four species by in situ hybridization to polytene chromosomes has shown in all cases both chromocentric hybridizations and a low number of sites (0-5) on the chromosomal arms. This number of sites is among the lowest observed in D. melanogaster and D. simulans when 1731 is compared with other retrotransposon families. In addition, we have observed species-specific patterns of the chromocentric hybridization signal, suggesting rapid modifications of the beta-heterochromatin components since the radiation of the melanogaster subgroup.      

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.     

Discovery of three new species of Drosophila obscura species group (Diptera: Drosophilidae) from Mount Kinabalu in Borneo

We describe three new species of the Drosophila obscura species group, D. hypercephala , D. hideakii and D. quadrangula Gao & Toda, spp. nov., all discovered from high altitudes (>1500 m) on Mt Kinabalu, Sabah, Malaysia. These are the southernmost distribution records of the D. obscura group in the Oriental Region. The three new species are each morphologically similar to different species of the Old World D. obscura species subgroup. This suggests that they are descendants of different founder species that independently colonized the tropics from the subtropical or temperate zones during the glacial ages. The new species D. hypercephala is hypercephalic not only in males but also, although less distinct, in females, unlike other hypercephalic drosophilid species in which only males are hypercephalic. This new finding of a hypercephalic species from the subgenus Sophophora , of which species are relatively easy to culture in the laboratory and include the model species such as D. melanogaster and D. pseudoobscura , will promote genetic analyses on the development and the evolution of hypercephaly.     

Regulatory elements involved in Drosophila Adh gene expression are conserved in divergent species and separate elements mediate expression in different tissues.

Alcohol dehydrogenase (ADH) is expressed in a complex temporal and spatial pattern from tandem promoters (proximal and distal) in Drosophila melanogaster, and from two closely linked genes (Adh-1 and Adh-2) in D. mulleri. The expression patterns of Adh-1 and the proximal promoter, and Adh-2 and the distal promoter are similar, but not identical. We show that the mulleri Adh genes are appropriately expressed when introduced into the melanogaster genome, indicating that the cis- and trans-acting elements which regulate the corresponding promoters are functionally equivalent in the two species. By analyzing the expression of in vitro generated mutants of the mulleri Adh locus, we define at least three regulatory regions of the mulleri Adh genes and show that different control elements mediate the expression of Adh in different tissues.     

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.      

Structure and Evolution of the Adh Genes of Drosophila Mojavensis

The nucleotide sequence of the Adh region of Drosophila mojavensis has been completed and the region found to contain a pseudogene, Adh-2 and Adh-1 arranged in that order. Comparison of the sequence divergence of these genes to one another and to the Adh region of Drosophila mulleri and other species has allowed the development of a model for the evolution of the duplication of the Adh genes. There have been two major events. An initial duplication of an Adh gene whose dual promoter structure was similar to Drosophila melanogaster, resulted in a species with two Adh genes, one of which may have had only a proximal promoter. A second duplication of this gene generated an Adh region containing three genes. It is proposed that one of these is the ancestral gene having dual promoters, while the other two possess only proximal promoters. Subsequent events have resulted in both a change in the regulation of Adh-2 such that it is expressed as if it had a "distal" type promoter and the mutational inactivation of the most upstream gene resulting in the creation of a pseudogene. The sequence of the D. mojavensis Adh region has also revealed the presence of an element which is composed of juxtaposed inverted imperfectly repeated elements. There is a surprising and not fully explainable strong similarity of the nucleotide sequence of the 5' flanking region of the pseudogene in D. mojavensis and D. mulleri.