Molecular evolution of a Y chromosome to autosome gene duplication in Drosophila

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila, many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila. First, we inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of them could develop on diet with toxin. Finally, we asked how α-amanitin tolerance might have evolved across the immigrans-tripunctata radiation, and inferred that toxin tolerance was ancestral in mushroom-feeding Drosophila and subsequently lost multiple times. Our findings expand our understanding of toxin tolerance across the immigrans-tripunctata radiation and emphasize the uniqueness of toxin tolerance in this adaptive radiation and the complexity of biochemical adaptations.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila , many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila . We first inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of these species could develop on diet with toxin. Third, we examined fly development on a diet containing a natural mix of toxins extracted from the Death Cap Amanita phalloides mushroom. Both tolerant and susceptible species developed on diet with this mix, though tolerant species survived at significantly higher concentrations. Finally, we asked how cyclopeptide tolerance might have evolved across the immigrans-tripunctata radiation and inferred that toxin tolerance was ancestral and subsequently lost multiple times. Our results suggest the evolutionary history of cyclopeptide tolerance is complex, and simply describing this trait as present or absent does not fully capture the occurrence or impact on this adaptive radiation. More broadly, the evolution of novelty can be more complex than previously thought, and that accurate descriptions of such novelties are critical in studies examining their evolution.

Gene duplication to the Y chromosome in Trinidadian Guppies

Differences in allele frequencies at autosomal genes between males and females in a population can result from two scenarios. First, unresolved sexual conflict over survival can produce allelic differentiation between the sexes. However, given the substantial mortality costs required to produce allelic differences between males and females at each generation, it remains unclear how many loci within the genome experience significant sexual conflict over survival. Alternatively, recent studies have shown that similarity between autosomal and Y sequences can create perceived allelic differences between the sexes. However, Y duplications are most likely in species with large non-recombining regions, in part because they simply represent larger targets for duplications. We assessed the genomes of 120 wild-caught guppies, which experience extensive predation- and pathogen-induced mortality and have a relatively small ancestral Y chromosome. We identified seven autosomal genes that show allelic differences between male and female adults. Five of these genes show clear evidence of whole or partial gene duplication between the Y chromosome and the autosomes. The remaining two genes show evidence of partial homology to the Y. Overall, our findings suggest that the guppy genome experiences a very low level of unresolved sexual conflict over survival, and instead the Y chromosome, despite its small ancestral size and recent origin, may nonetheless accumulate genes with male-specific functions.

Complex Evolutionary History of the Y Chromosome in Flies of the Drosophila obscura Species Group

The Drosophila obscura species group shows dramatic variation in karyotype, including transitions among sex chromosomes. Members of the affinis and pseudoobscura subgroups contain a neo-X chromosome (a fusion of the X with an autosome), and ancestral Y genes have become autosomal in species harboring the neo-X. Detailed analysis of species in the pseudoobscura subgroup revealed that ancestral Y genes became autosomal through a translocation to the small dot chromosome. Here, we show that the Y-dot translocation is restricted to the pseudoobscura subgroup, and translocation of ancestral Y genes in the affinis subgroup likely followed a different route. We find that most ancestral Y genes have translocated to unique autosomal or X-linked locations in different taxa of the affinis subgroup, and we propose a dynamic model of sex chromosome formation and turnover in the obscura species group. Our results suggest that Y genes can find unique paths to escape unfavorable genomic environments that form after sex chromosome–autosome fusions.

Novel Y Chromosome Retrocopies in Canids Revealed through a Genome-Wide Association Study for Sex

The lack of an annotated reference sequence for the canine Y chromosome has limited evolutionary studies, as well as our understanding of the role of Y-linked sequences in phenotypes with a sex bias. In genome-wide association studies (GWASs), we observed spurious associations with autosomal SNPs when sex was unbalanced in case-control cohorts and hypothesized that a subset of SNPs mapped to autosomes are in fact sex-linked. Using the Illumina 230K CanineHD array in a GWAS for sex, we identified SNPs that amplify in both sexes but possess significant allele frequency differences between males and females. We found 48 SNPs mapping to 14 regions of eight autosomes and the X chromosome that are Y-linked, appearing heterozygous in males and monomorphic in females. Within these 14 regions are eight genes: three autosomal and five X-linked. We investigated the autosomal genes (MITF, PPP2CB, and WNK1) and determined that the SNPs are diverged nucleotides in retrocopies that have transposed to the Y chromosome. MITFY and WNK1Y are expressed and appeared recently in the Canidae lineage, whereas PPP2CBY represents a much older insertion with no evidence of expression in the dog. This work reveals novel canid Y chromosome sequences and provides evidence for gene transposition to the Y from autosomes and the X.

A fast-evolving X-linked duplicate of importin-α2 is overexpressed in sex-ratio drive in Drosophila neotestacea

Selfish genetic elements that manipulate gametogenesis to achieve a transmission advantage are known as meiotic drivers. Sex-ratio X chromosomes (SR) are meiotic drivers that prevent the maturation of Y-bearing sperm in male carriers to result in the production of mainly female progeny. The spread of an SR chromosome can affect host genetic diversity and genome evolution, and can even cause host extinction if it reaches sufficiently high prevalence. Meiotic drivers have evolved independently many times, though only in a few cases is the underlying genetic mechanism known. In this study we use a combination of transcriptomics and population genetics to identify widespread expression differences between the standard (ST) and sex-ratio (SR) X chromosomes of the fly Drosophila neotestacea. We found the X chromosome is enriched for differentially expressed transcripts and that many of these X-linked differentially expressed transcripts had elevated K_a /K_s values between ST and SR, indicative of potential functional differences. We identified a set of candidate transcripts, including a testis-specific, X-linked duplicate of the nuclear transport gene importin-α2 that is overexpressed in SR. We find suggestions of positive selection in the lineage leading to the duplicate and that its molecular evolutionary patterns are consistent with relaxed purifying selection in ST. As these patterns are consistent with involvement in the mechanism of drive in this species, this duplicate is a strong candidate worthy of further functional investigation. Nuclear transport may be a common target for genetic conflict, as the mechanism of the autosomal Segregation Distorter drive system in D. melanogaster involves the same pathway.

An investigation of Y chromosome incorporations in 400 species of Drosophila and related genera

Y chromosomes are widely believed to evolve from a normal autosome through a process of massive gene loss (with preservation of some male genes), shaped by sex-antagonistic selection and complemented by occasional gains of male-related genes. The net result of these processes is a male-specialized chromosome. This might be expected to be an irreversible process, but it was found in 2005 that the Drosophila pseudoobscura Y chromosome was incorporated into an autosome. Y chromosome incorporations have important consequences: a formerly male-restricted chromosome reverts to autosomal inheritance, and the species may shift from an XY/XX to X0/XX sex-chromosome system. In order to assess the frequency and causes of this phenomenon we searched for Y chromosome incorporations in 400 species from Drosophila and related genera. We found one additional large scale event of Y chromosome incorporation, affecting the whole montium subgroup (40 species in our sample); overall 13% of the sampled species (52/400) have Y incorporations. While previous data indicated that after the Y incorporation the ancestral Y disappeared as a free chromosome, the much larger data set analyzed here indicates that a copy of the Y survived as a free chromosome both in montium and pseudoobscura species, and that the current Y of the pseudoobscura lineage results from a fusion between this free Y and the neoY. The 400 species sample also showed that the previously suggested causal connection between X-autosome fusions and Y incorporations is, at best, weak: the new case of Y incorporation (montium) does not have X-autosome fusion, whereas nine independent cases of X-autosome fusions were not followed by Y incorporations. Y incorporation is an underappreciated mechanism affecting Y chromosome evolution; our results show that at least in Drosophila it plays a relevant role and highlight the need of similar studies in other groups.

In contrast to other chromosomes (X and autosomes), which are present in males and females, Y chromosomes spend all time in males. Hence it is not surprising that along evolution they became male specialized, e.g., containing a disproportionate amount of male-fertility genes. Interestingly it was found in 2005 that in Drosophila pseudoobscura the Y chromosome reverted to "male-female existence", being incorporated into an autosome. These "Y chromosome incorporations" have important consequences on sex-chromosome evolution, and allow the study of the evolutionary forces that shaped Y chromosomes as they act backwards. As D. pseudoobscura was the second Drosophila species investigated in this respect, it is likely that other cases exist, and that perhaps it is a common phenomenon. In order to answer this question we studied 400 Drosophila species. We found one additional case of Y incorporation, which occurred in the ancestor of Drosophila montium, and currently affects a large number of species; overall 13% of the species we sampled (52/400) have Y incorporations. We also found that a previously suggested cause of Y incorporations (X-autosome fusions) is not a general explanation. Our results show that in Drosophila Y incorporations play a relevant role and highlight the need of similar studies in other groups.

A phylogenetic examination of host use evolution in the quinaria and testacea groups of Drosophila

Adaptive radiations provide an opportunity to examine complex evolutionary processes such as ecological specialization and speciation. While a well-resolved phylogenetic hypothesis is critical to completing such studies, the rapid rates of evolution in these groups can impede phylogenetic studies. Here we study the quinaria and testacea species groups of the immigrans-tripunctata radiation of Drosophila, which represent a recent adaptive radiation and are a developing model system for ecological genetics. We were especially interested in understanding host use evolution in these species. In order to infer a phylogenetic hypothesis for this group we sampled loci from both the nuclear genome and the mitochondrial DNA to develop a dataset of 43 protein-coding loci for these two groups along with their close relatives in the immigrans-tripunctata radiation. We used this dataset to examine their evolutionary relationships along with the evolution of feeding behavior. Our analysis recovers strong support for the monophyly of the testacea but not the quinaria group. Results from our ancestral state reconstruction analysis suggests that the ancestor of the testacea and quinaria groups exhibited mushroom-feeding. Within the quinaria group, we infer that transition to vegetative feeding occurred twice, and that this transition did not coincide with a genome-wide change in the rate of protein evolution.

High rate of translocation-based gene birth on the Drosophila Y chromosome

Using a powerful method that uses inexpensive short reads to detect Y-linked transfers, we show that gene traffic onto the Drosophila Y chromosome is 10 times more frequent than previously thought and includes the first Y-linked retrocopies discovered in these taxa. All 25 identified Y-linked gene transfers were relatively young (<1 million years old), although most appear to be pseudogenes because only three of these transfers show signs of purifying selection. Our method provides compelling evidence that the Drosophila Y chromosome is a highly challenging and dynamic genetic environment that is capable of rapidly diverging between species and promises to reveal fundamental insights into Y chromosome evolution across many taxa.

The Y chromosome is a unique genetic environment defined by a lack of recombination and male-limited inheritance. The Drosophila Y chromosome has been gradually acquiring genes from the rest of the genome, with only seven Y-linked genes being gained over the past 63 million years (0.12 gene gains per million years). Using a next-generation sequencing (NGS)-powered genomic scan, we show that gene transfers to the Y chromosome are much more common than previously suspected: at least 25 have arisen across three Drosophila species over the past 5.4 million years (1.67 per million years for each lineage). The gene transfer rate is significantly lower in Drosophila melanogaster than in the Drosophila simulans clade, primarily due to Y-linked retrotranspositions being significantly more common in the latter. Despite all Y-linked gene transfers being evolutionarily recent (<1 million years old), only three showed evidence for purifying selection (ω ≤ 0.14). Thus, although the resulting Y-linked functional gene acquisition rate (0.25 new genes per million years) is double the longer-term estimate, the fate of most new Y-linked genes is defined by rapid degeneration and pseudogenization. Our results show that Y-linked gene traffic, and the molecular mechanisms governing these transfers, can diverge rapidly between species, revealing the Drosophila Y chromosome to be more dynamic than previously appreciated. Our analytical method provides a powerful means to identify Y-linked gene transfers and will help illuminate the evolutionary dynamics of the Y chromosome in Drosophila and other species.

X chromosome drive in a widespread Palearctic woodland fly, Drosophila testacea

Selfish genes that bias their own transmission during meiosis can spread rapidly in populations, even if they contribute negatively to the fitness of their host. Driving X chromosomes provide a clear example of this type of selfish propagation. These chromosomes have important evolutionary and ecological consequences, and can be found in a broad range of taxa including plants, mammals and insects. Here, we report a new case of X chromosome drive (X drive) in a widespread woodland fly, Drosophila testacea. We show that males carrying the driving X (SR males) sire 80-100% female offspring and possess a diagnostic X chromosome haplotype that is perfectly associated with the sex ratio distortion phenotype. We find that the majority of sons produced by SR males are sterile and appear to lack a Y chromosome, suggesting that meiotic defects involving the Y chromosome may underlie X drive in this species. Abnormalities in sperm cysts of SR males reflect that some spermatids are failing to develop properly, confirming that drive is acting during gametogenesis. By screening wild-caught flies using progeny sex ratios and a diagnostic marker, we demonstrate that the driving X is present in wild populations at a frequency of ~ 10% and that suppressors of drive are segregating in the same population. The testacea species group appears to be a hot spot for X drive, and D. testacea is a promising model to compare driving X chromosomes in closely related species, some of which may even be younger than the chromosomes themselves.

Genomic changes following the reversal of a Y chromosome to an autosome in Drosophila pseudoobscura

Robertsonian translocations resulting in fusions between sex chromosomes and autosomes shape karyotype evolution by creating new sex chromosomes from autosomes. These translocations can also reverse sex chromosomes back into autosomes, which is especially intriguing given the dramatic differences between autosomes and sex chromosomes. To study the genomic events following a Y chromosome reversal, we investigated an autosome‐Y translocation in Drosophila pseudoobscura. The ancestral Y chromosome fused to a small autosome (the dot chromosome) approximately 10–15 Mya. We used single molecule real‐time sequencing reads to assemble the D. pseudoobscura dot chromosome, including this Y‐to‐dot translocation. We find that the intervening sequence between the ancestral Y and the rest of the dot chromosome is only ∼78 Kb and is not repeat‐dense, suggesting that the centromere now falls outside, rather than between, the fused chromosomes. The Y‐to‐dot region is 100 times smaller than the D. melanogaster Y chromosome, owing to changes in repeat landscape. However, we do not find a consistent reduction in intron sizes across the Y‐to‐dot region. Instead, deletions in intergenic regions and possibly a small ancestral Y chromosome size may explain the compact size of the Y‐to‐dot translocation.

Occasional recombination of a selfish X-chromosome may permit its persistence at high frequencies in the wild

The sex-ratio X-chromosome (SR) is a selfish chromosome that promotes its own transmission to the next generation by destroying Y-bearing sperm in the testes of carrier males. In some natural populations of the fly Drosophila neotestacea, up to 30% of the X-chromosomes are SR chromosomes. To investigate the molecular evolutionary history and consequences of SR, we sequenced SR and standard (ST) males at 11 X-linked loci that span the ST X-chromosome and at seven arbitrarily chosen autosomal loci from a sample of D. neotestacea males from throughout the species range. We found that the evolutionary relationship between ST and SR varies among individual markers, but genetic differentiation between SR and ST is chromosome-wide and likely due to large chromosomal inversions that suppress recombination. However, SR does not consist of a single multilocus haplotype: we find evidence for gene flow between ST and SR at every locus assayed. Furthermore, we do not find long-distance linkage disequilibrium within SR chromosomes, suggesting that recombination occurs in females homozygous for SR. Finally, polymorphism on SR is reduced compared to that on ST, and loci displaying signatures of selection on ST do not show similar patterns on SR. Thus, even if selection is less effective on SR, our results suggest that gene flow with ST and recombination between SR chromosomes may prevent the accumulation of deleterious mutations and allow its long-term persistence at relatively high frequencies.

Patterns of reproductive isolation in the Drosophila subquinaria complex: can reinforced premating isolation cascade to other species?

The reinforcement of premating barriers due to reduced hybrid fitness in sympatry may cause secondary sexual isolation within a species as a by-product. Consistent with this, in the fly Drosophila subquinaria, females that are sympatric with D. recens mate at very low rates not only with D. recens, but also with conspecific D. subquinaria males from allopatry. Here, we ask if these effects of reinforcement cascade more broadly to affect sexual isolation with other closely related species. We assay reproductive isolation of these species with D. transversa and find that choosy D. subquinaria females from the region sympatric with D. recens discriminate strongly against male D. transversa, whereas D. subquinaria from the allopatric region do not. This increased sexual isolation cannot be explained by natural selection to avoid mating with this species, as they are allopatric in geographic range and we do not identify any intrinsic postzygotic isolation between D. subquinaria and D. transversa. Variation in epicuticular hydrocarbons, which are used as mating signals in D. subquinaria, follow patterns of premating isolation: D. transversa and allopatric D. subquinaria are most similar to each other and differ from sympatric D. subquinaria, and those of D. recens are distinct from the other two species. We suggest that the secondary effects of reinforcement may cascade to strengthen reproductive isolation with other species that were not a target of selection. These effects may enhance the divergence that occurs in allopatry to help explain why some species are already sexually isolated upon secondary contact.

Asymmetrical sexual isolation but no postmating isolation between the closely related species Drosophila suboccidentalis and Drosophila occidentalis

BACKGROUND

During the speciation process several types of isolating barriers can arise that limit gene flow between diverging populations. Studying recently isolated species can inform our understanding of how and when these barriers arise, and which barriers may be most important to limiting gene flow. Here we focus on Drosophila suboccidentalis and D. occidentalis, which are closely related mushroom-feeding species that inhabit western North America and are not known to overlap in geographic range. We investigate patterns of reproductive isolation between these species, including premating, postmating prezygotic, and postzygotic barriers to gene flow.

RESULTS

Using flies that originate from a single population of each species, we find that the strength of premating sexual isolation between these species is asymmetric: while D. occidentalis females mate with D. suboccidentalis males at a reduced but moderate rate, D. suboccidentalis females discriminate strongly against mating with D. occidentalis males. Female hybrids will mate at high rates with males of either species, indicating that this discrimination has a recessive genetic basis. Hybrid males are accepted by females of both species. We do not find evidence for postmating prezygotic or postzygotic isolating barriers, as females use the sperm of heterospecific males and both male and female hybrids are fully fertile.

CONCLUSIONS

Premating isolation is substantial but incomplete, and appears to be the primary form of reproductive isolation between these species. If these species do hybridize, the lack of postzygotic barriers may allow for gene flow between them. Given that these species are recently diverged and are not known to be sympatric, the level of premating isolation is relatively strong given the lack of intrinsic postzygotic isolation. Further work is necessary to characterize the geographic and genetic variation in reproductive isolating barriers, as well as to determine the factors that drive reproductive isolation and the consequences that isolating barriers as well as geographic isolation have had on patterns of gene flow between these species.

The relations between recombination rate and patterns of molecular variation and evolution in Drosophila

Genetic recombination affects levels of variability and the efficacy of selection because natural selection acting at one site affects evolutionary processes at linked sites. The variation in local recombination rates across the Drosophila genome provides excellent material for testing hypotheses concerning the evolutionary consequences of recombination. The current state of knowledge from studies of Drosophila genomics and population genetics is reviewed here. Selection at linked sites has influenced the relations between recombination rates and patterns of molecular variation and evolution, such that higher rates of recombination are associated with both higher levels of variability and a greater efficacy of selection. It seems likely that background selection against deleterious mutations is a major factor contributing to these patterns in genome regions in which crossing over is rare or absent, whereas selective sweeps of positively selected mutations probably play an important role in regions with crossing over.

Positive and purifying selection on the Drosophila Y chromosome

Y chromosomes, with their reduced effective population size, lack of recombination, and male-limited transmission, present a unique collection of constraints for the operation of natural selection. Male-limited transmission may greatly increase the efficacy of selection for male-beneficial mutations, but the reduced effective size also inflates the role of random genetic drift. Together, these defining features of the Y chromosome are expected to influence rates and patterns of molecular evolution on the Y as compared with X-linked or autosomal loci. Here, we use sequence data from 11 genes in 9 Drosophila species to gain insight into the efficacy of natural selection on the Drosophila Y relative to the rest of the genome. Drosophila is an ideal system for assessing the consequences of Y-linkage for molecular evolution in part because the gene content of Drosophila Y chromosomes is highly dynamic, with orthologous genes being Y-linked in some species whereas autosomal in others. Our results confirm the expectation that the efficacy of natural selection at weakly selected sites is reduced on the Y chromosome. In contrast, purifying selection on the Y chromosome for strongly deleterious mutations does not appear to be compromised. Finally, we find evidence of recurrent positive selection for 4 of the 11 genes studied here. Our results thus highlight the variable nature of the mode and impact of natural selection on the Drosophila Y chromosome.

Two genomic regions together cause dark abdominal pigmentation in Drosophila tenebrosa

Pigmentation is a rapidly evolving trait that is under both natural and sexual selection in many organisms. In the quinaria group of Drosophila, nearly all of the 30 species have an abdomen that is light in color with distinct markings; D. tenebrosa is the exception in that it has a completely melanic abdomen with no visible markings. In this study, we use a combination of quantitative genetic and candidate gene approaches to investigate the genetic basis of abdominal pigmentation in D. tenebrosa. We find that abdominal pigmentation is invariant across wild-caught lines of D. tenebrosa and is not sexually dimorphic. Quantitative genetic mapping utilizing crosses between D. tenebrosa and the light-colored D. suboccidentalis indicates that two genomic regions together underlie abdominal pigmentation, including the X-chromosome and an autosome (Muller Element C/E). Further support for their central importance in pigmentation is that experimental introgression of one phenotype into the other species, in either direction, results in introgression of these two genomic regions. Finally, the expression of the X-linked gene yellow in the pupae exactly foreshadows the adult melanization pattern in the abdomen of both species, suggesting that changes in the regulation of yellow are important for the phenotypic divergence of D. tenebrosa from the rest of the quinaria group. These results contribute to a body of work that demonstrates how changes in expression of highly conserved genes can cause substantial phenotypic differences even between closely related species.

Genomic conflict drives patterns of X-linked population structure in Drosophila neotestacea

Intragenomic conflict has the potential to cause widespread changes in patterns of genetic diversity and genome evolution. In this study, we investigate the consequences of sex-ratio (SR) drive on the population genetic patterns of the X-chromosome in Drosophila neotestacea. An SR X-chromosome prevents the maturation of Y-bearing sperm during male spermatogenesis and thus is transmitted to ~100% of the offspring, nearly all of which are daughters. Selection on the rest of the genome to suppress SR can be strong, and the resulting conflict over the offspring sex ratio can result in the accumulation of multiple loci on the X-chromosome that are necessary for the expression of drive. We surveyed variation at 12 random X-linked microsatellites across 16 populations of D. neotestacea that range in SR frequency from 0% to 30%. First, every locus was differentiated between SR and wild-type chromosomes, and this drives genetic structure at the X-chromosome. Once the association with SR is accounted for, the patterns of differentiation among populations are similar to the autosomes. Second, within wild-type chromosomes, the relative heterozygosity is reduced in populations with an increased prevalence of drive, and the heterozygosity of SR chromosomes is higher than expected based on its prevalence. The combination of the relatively high prevalence of SR drive and the structuring of polymorphism between the SR and wild-type chromosomes suggests that genetic conflict because of SR drive has had significant consequences on the patterns of X-linked polymorphism and thus also probably affects the tempo of X-chromosome evolution in D. neotestacea.

Wolbachia-mediated persistence of mtDNA from a potentially extinct species

Drosophila quinaria is polymorphic for infection with Wolbachia, a maternally transmitted endosymbiont. Wolbachia-infected individuals carry mtDNA that is only distantly related to the mtDNA of uninfected individuals, and the clade encompassing all mtDNA haplotypes within D. quinaria also includes the mtDNA of several other species of Drosophila. Nuclear gene variation reveals no difference between the Wolbachia-infected and uninfected individuals of D. quinaria, indicating that they all belong to the same interbreeding biological species. We suggest that the Wolbachia and the mtDNA with which it is associated were derived via interspecific hybridization and introgression. The sequences in the Wolbachia and the associated mtDNA are ≥6% divergent from those of any known Drosophila species. Thus, in spite of nearly complete species sampling, the sequences from which these mitochondria were derived remain unknown, raising the possibility that the donor species is extinct. The association between Wolbachia infection and mtDNA type within D. quinaria suggests that Wolbachia may be required for the continued persistence of the mtDNA from an otherwise extinct Drosophila species. We hypothesize that pathogen-protective effects conferred by Wolbachia operate in a negative frequency-dependent manner, thus bringing about a stable polymorphism for Wolbachia infection.