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.     

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 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.      

Distinct activation patterns of EGF receptor signaling in the homoplastic evolution of eggshell morphology in genus Drosophila

Homoplasy is a phenomenon in which organisms in different phylogenetic groups independently acquire similar traits. However, it is largely unknown how developmental mechanisms are altered to give rise to homoplasy. In the genus Drosophila, all species of the subgenus Sophophora, including Drosophila (D.) melanogaster, have eggshells with two dorsal appendages (DAs); most species in the subgenus Drosophila, including D. virilis, and in the subgenus Dorsilopha, have four-DAs. D. melanica belongs to the Drosophila subgenus, but has two-DAs, and phylogenetic analyses suggest that it acquired this characteristic independently. The patterning of the DAs is tightly regulated by epidermal growth factor receptor (EGFR) signaling in D. melanogaster. Previous studies suggested that a change in the EGFR signal activation pattern could have led to the divergence in DA number between D. melanogaster and D. virilis. Here, we compared the patterns of EGFR signal activation across the Drosophila subgenera by immunostaining for anti-activated MAP kinase (MAPK). Our analysis revealed distinct patterns of EGFR signal activation in each subgenus that was consistent with their phylogenetic relationship. In addition, the number of DAs always corresponded to the number of EGFR signaling activation domains in two, three, and four-DA species. Despite their common two-DA characteristic, the EGFR signaling activation pattern in D. melanica diverged significantly from that of species in the subgenus Sophophora. Our results suggest that acquisition of the homoplastic two-DA characteristic could be explained by modifications of the EGFR signaling system in the genus Drosophila that occurred independently and at least twice during evolution.     

Interspecific analysis of the glycosidases of the sperm plasma membrane in Drosophila

We have studied the presence of four sperm glycosidases, alpha-mannosidase, alpha-L- fucosidase and two beta-hexosaminidase isoforms, in 11 species of the genus Drosophila spanning approximately an evolutionary 60 MY period, and in Scaptodrosophila lebanonensis, belonging to the ancestor genus Scaptodrosophila. These enzymes had been previously identified in Drosophila melanogaster as putative receptors for glycoconjugates of the egg surface. Alpha-mannosidase and beta-hexosaminidases are intrinsic proteins of the sperm plasma membrane in species closely related to D. melanogaster as well as in the divergent species D. willistoni, D. hydei, D. virilis, and S. lebanonensis. Alpha-L-fucosidase is restricted to the species of the genus Drosophila. Alpha-mannosidase and beta-hexosaminidases have been purified and characterized in all species. Their molecular masses and optimal pHs are similar in all the species, whereas interspecific differences in enzyme activities were detected. Cross-species comparison of kinetic parameters indicated a relationship between enzyme efficiency and phylogenetic relatedness. Beta-hexosaminidases were the most efficient enzymes. Lectin cytochemistry suggested the presence of carbohydrate residues complementary to the glycosidases on the eggshell at the site of sperm entry in all species. Bioinformatic analysis of the coding sequences of beta-hexosaminases and alpha-L-fucosidase and of their predicted products showed no evidence of positive selection of the genes coding for these enzymes and a high degree of sequence identity of the predicted polypeptides among the species of the genus Drosophila. Collectively, our findings indicate that the Drosophila sperm glycosidases are structurally and functionally conserved and strengthen the hypothesis of their involvement in the interactions with the egg surface.     

Convergent intron loss of MRP1 in Drosophila and mosquito species

Previous study revealed that the MRP1 gene ortholog DMRP1/CG6214 of Drosophila melanogaster contains 12 exons in the coding region. In the current study, the genes of DMRP1/CG6214 from D. melanogaster and Drosophila virilis were compared, and the result indicated that D. virilis had an extra intron located in exon 2, implying that intron loss or gain might have occurred at this locus. To track the evolution of the extra intron (Intron Z), orthologous nucleotide sequences of 37 arthropod species were cloned or annotated. Based on phylogenetic analysis, we found that Intron Z should present in the common ancestor of arthropod species, more than 420 Ma. In addition, we found that Sophophora subgenus species and mosquito (Culex pipiens) lost Intron Z independently, showing evolutionary convergence.     

Ultrastructural features of spermatozoa and their phylogenetic application in Zaprionus (Diptera,Drosophilidae)

The genus Zaprionus consists of approximately 60 species of drosophilids that are native to the Afrotropical region. The phylogenetic position of Zaprionus within the Drosophilidae family is still unresolved. In the present study, ultrastructural features of spermatozoa of 6 species of Zaprionus as well as the species Drosophila willistoni and Scaptodrosophila latifasciaeformis were analyzed. The ultrastructure revealed that the species have the same flagellar ultrastructure. Two mitochondrial derivatives, one larger than the other, close to the axoneme were present, primarily in D. willistoni (subgenus Sophophora). Except for Z. davidi and Z. tuberculatus, the analyzed species had paracrystalline material in both mitochondrial derivatives. Moreover, the testes showed 64 spermatozoa per bundle in all of the species. In the cluster analysis, 6 Zaprionus species were grouped closely, but there were some incongruent positions in the cladogram. The results indicated that sperm ultrastructure is an important tool for elucidating the phylogeny and taxonomy of insects.     

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.     

X chromosome DNA variation in Drosophila virilis.

Comparisons of polymorphism patterns between distantly related species are essential in order to determine their generality. However, most work on the genus Drosophila has been done only with species of the subgenus Sophophora. In the present work, we have sequenced one intron and surrounding coding sequences of 6 X-linked genes (chorion protein s36, elav, fused, runt, suppressor of sable and zeste) from 21 strains of wild-type Drosophila virilis (subgenus Drosophila). From these data, we have estimated the average level of DNA polymorphism, inferred the effective population size and population structure of this species, and compared the results with those obtained for other Drosophila species. There is no reduction in variation at two loci close to the centromeric heterochromatin, in contrast to Drosophila melanogaster.     

Two independent duplications forming the Cyp307a genes in Drosophila

The conserved relationship between orthologs of many cytochrome P450 genes involved in ecdysone synthesis is not reflected in the evolution of the Drosophila Cyp307a genes. In Drosophila melanogaster Cyp307a1 (spook) and Cyp307a2 (spookier) both play essential roles in ecdysone synthesis and may possess biochemically redundant functions. Using phylogenetic analyses we show that the Drosophila Cyp307a genes were formed from two independent duplication events depicting a complicated evolutionary scenario. An initial duplication, from a Cyp307a2 ancestral gene produced the Cyp307a1 gene that has been maintained only in the Sophophoran subgenus. A second duplication in the Drosophila subgenus formed an additional paralog, Cyp307a3. Microsynteny is conserved for Cyp307a2 throughout the Drosophila species, but is not conserved between Cyp307a1 and Cyp307a3. These are located in different genomic positions in the Sophophora and Drosophila subgenera, respectively. Cyp307a3 appears to encode a functional gene product and is expressed in a different spatial and temporal manner to Cyp307a1. This suggests some level of functional divergence between the Cyp307a paralogs in different Drosophila species.     

Chromosomal evolution of elements B and C in the Sophophora subgenus of Drosophila: evolutionary rate and polymorphism

The locations of 77 markers along the chromosomal elements B (41 markers) and C (36 markers) of Drosophila subobscura, D. pseudoobscura, and D. melanogaster were obtained by in situ hybridization on polytene chromosomes. In comparisons between D. subobscura and D. pseudoobscura, 10 conserved segments (accounting for 32% of the chromosomal length) were detected on element B and eight (17% of the chromosomal length) on element C. The fixation rate of paracentric inversions inferred by a maximum likelihood approach differs significantly between elements. Muller's element C (0.17 breakpoints/Mb/million years) is evolving two times faster than element B (0.08 breakpoints/Mb/million years). This difference in the evolutionary rate is paralleled by differences in the extent of chromosomal polymorphism in the corresponding lineages. Element C is highly polymorphic in D. subobscura, D. pseudoobscura, and in other obscura group species such as D. obscura and D. athabasca. In contrast, the level of polymorphism in element B is much lower in these species. The fixation rates of paracentric inversions estimated in the present study between species of the Sophophora subgenus are the highest estimates so far reported in the genus for the autosomes. At the subgenus level, there is also a parallelism between the high fixation rate and the classical observation that the species of the Sophophora subgenus tend to be more polymorphic than the species of the Drosophila subgenus. Therefore, the detected relationship between level of polymorphism and evolutionary rate might be a general characteristic of chromosomal evolution in the genus Drosophila.     

Drosophila mediopunctata P elements: a new example of horizontal transfer.

Sequences homologous to the P element of Drosophila melanogaster were previously identified in Drosophila mediopunctata, a member of the tripunctata group, subgenus Drosophila. We report here that the P element is present in about three to five copies in the D. mediopunctata genome. While one of the insertion sites appears to be fixed, others may be polymorphic, indicating relatively recent P element activity. Phylogenetic analysis revealed that the D. mediopunctata element belongs to the canonical subfamily of P elements and that divergence of the D. mediopunctata element from other members of this subfamily ranges from 2% to 5% at the nucleotide level. This is the first report of a canonical P element outside the subgenus Sophophora. Based primarily on the striking incongruence between P element and host species phylogenies, the presence of a canonical P element in D. mediopunctata is most likely explained by horizontal transfer between species.     

The rate and pattern of nucleotide substitution in Drosophila mitochondrial DNA.

The nucleotide sequences of a segment of mitochondrial DNA (mtDNA) have been determined for nine species or subspecies of the subgenus Drosophila of the genus Drosophila. This segment contains two complete protein-coding genes (i.e., NADH dehydrogenase subunit 1 and cytochrome b) and a transfer RNA gene (tRNA(ser)). The G+C content at third-codon positions for the two protein-coding genes was 1.5 times higher than that in the D. melanogaster species group, which belongs to the subgenus Sophophora. However, there was a substantial difference between the nucleotide frequencies of G and C. The number of nucleotide substitutions per silent site was more than three times higher than that for nuclear DNA, although it was only 60% of that for mammalian mtDNA. Both parametric and nonparametric analyses revealed a strong transition-transversion bias in nucleotide substitution, as was observed in mammalian mtDNA. Moreover, the rate of substitution of A and T for G and C is higher than that for the opposite direction. This bias seems to be responsible for the extremely A+T-rich base composition of Drosophila mtDNA. It is also noted that the rate of transitional change between A and G is higher than that between T and C.     

Duplicative and conservative transpositions of larval serum protein 1 genes in the genus Drosophila

Interspecific comparative molecular analyses of transposed genes and their flanking regions can help to elucidate the time, direction, and mechanism of gene transposition. In the Drosophila melanogaster genome, three Larval serum protein 1 (Lsp1) genes (alpha, beta and gamma) are present and each of them is located on a different chromosome, suggesting multiple transposition events. We have characterized the molecular organization of Lsp1 genes in D. buzzatii, a species of the Drosophila subgenus and in D. pseudoobscura, a species of the Sophophora subgenus. Our results show that only two Lsp1 genes (beta and gamma) exist in these two species. The same chromosomal localization and genomic organization, different from that of D. melanogaster, is found in both species for the Lsp1beta and Lsp1gamma genes. Overall, at least two duplicative and two conservative transpositions are necessary to explain the present chromosomal distribution of Lsp1 genes in the three Drosophila species. Clear evidence for implication of snRNA genes in the transposition of Lsp1beta in Drosophila has been found. We suggest that an ectopic exchange between highly similar snRNA sequences was responsible for the transposition of this gene. We have also identified the putative cis-acting regulatory regions of these genes, which seemingly transposed along with the coding sequences.     

An alpha-L-fucosidase potentially involved in fertilization is present on Drosophila spermatozoa surface

Drosophila is emerging as a model organism to investigate egg fertilization in insects and the possible conservation of molecular mechanisms of gamete interactions demonstrated in higher organisms. This study shows that the spermatozoa of several species of Drosophila belonging to the melanogaster group have a plasma membrane associated alpha-L-fucosidase with features in common with alpha-L-fucosidases from sperm of other animals, including mammals. The enzyme has been purified and completely characterized in D. ananassae, because of its stability in this species. The sperm alpha-L-fucosidase is an integral protein terminally mannosylated, with the catalytic site oriented toward the extracellular space. It has a M(r) of 256 kDa and a multimeric structure made up by subunits of 48 and 55 kDa. Enzyme characterization included kinetic properties, pI, optimal pH, and thermal stability. A soluble form of the enzyme similar to the sperm associated alpha-L-fucosidase is secreted by the seminal vesicles. Synthetic peptides designed from the deduced product of the D. melanogaster gene encoding an alpha-L-fucosidase were used to raise a specific polyclonal antibody. Immunofluorescence labeling of spermatozoa showed that the enzyme is present in the sperm plasma membrane overlying the acrosome and the tail. Lectin cytochemistry analysis of the egg surface indicated that alpha-L-fucose terminal residues are present on the chorion with a strongly polarized localization on the micropyle. The alpha-L-fucosidase of Drosophila sperm plasma membrane appears to be potentially involved in gamete recognition by interacting with its glycoside ligands present on the egg surface at the site of sperm entry.     

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.