Drosophila Evolution over Space and Time (DEST): A New Population Genomics Resource.

Drosophila melanogaster is a leading model in population genetics and genomics, and a growing number of whole-genome data sets from natural populations of this species have been published over the last years. A major challenge is the integration of disparate data sets, often generated using different sequencing technologies and bioinformatic pipelines, which hampers our ability to address questions about the evolution of this species. Here we address these issues by developing a bioinformatics pipeline that maps pooled sequencing (Pool-Seq) reads from D. melanogaster to a hologenome consisting of fly and symbiont genomes and estimates allele frequencies using either a heuristic (PoolSNP) or a probabilistic variant caller (SNAPE-pooled). We use this pipeline to generate the largest data repository of genomic data available for D. melanogaster to date, encompassing 271 previously published and unpublished population samples from over 100 locations in >20 countries on four continents. Several of these locations have been sampled at different seasons across multiple years. This data set, which we call Drosophila Evolution over Space and Time (DEST), is coupled with sampling and environmental metadata. A web-based genome browser and web portal provide easy access to the SNP data set. We further provide guidelines on how to use Pool-Seq data for model-based demographic inference. Our aim is to provide this scalable platform as a community resource which can be easily extended via future efforts for an even more extensive cosmopolitan data set. Our resource will enable population geneticists to analyze spatiotemporal genetic patterns and evolutionary dynamics of D. melanogaster populations in unprecedented detail.

Fitness consequences of a non-recombining sex-ratio drive chromosome can explain its prevalence in the wild.

Understanding the pleiotropic consequences of gene drive systems on host fitness is essential to predict their spread through a host population. Here, we study sex-ratio (SR) X-chromosome drive in the fly Drosophila recens, where SR causes the death of Y-bearing sperm in male carriers. SR males only sire daughters, which all carry SR, thus giving the chromosome a transmission advantage. The prevalence of the SR chromosome appears stable, suggesting pleiotropic costs. It was previously shown that females homozygous for SR are sterile, and here, we test for additional fitness costs of SR. We found that females heterozygous for SR have reduced fecundity and that male SR carriers have reduced fertility in conditions of sperm competition. We then use our fitness estimates to parametrize theoretical models of SR drive and show that the decrease in fecundity and sperm competition performance can account for the observed prevalence of SR in natural populations. In addition, we found that the expected equilibrium frequency of the SR chromosome is particularly sensitive to the degree of multiple mating and performance in sperm competition. Together, our data suggest that the mating system of the organism should be carefully considered during the development of gene drive systems.

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.

Ongoing hybridization obscures phylogenetic relationships in the Drosophila subquinaria species complex.

Inferring evolutionary relationships among recently diverged lineages is necessary to understand how isolating barriers produce independent lineages. Here, we investigate the phylogenetic relationships between three incompletely isolated and closely related mushroom-feeding Drosophila species. These species form the Drosophila subquinaria species complex and consist of one Eurasian species (D. transversa) and two widespread North American species (D. subquinaria and D. recens) that are sympatric in central Canada. Although patterns of pre- and post-mating isolation among these species are well characterized, previous work on their phylogenetic relationships is limited and conflicting. In this study, we generated a multi-locus data set of 29 loci from across the genome sequenced in a population sample from each species, and then, we inferred species relationships and patterns of introgression. We find strong statistical support that D. subquinaria is paraphyletic, showing that samples from the geographic region sympatric with D. recens are most closely related to D. recens, whereas samples from the geographic region allopatric with D. recens are most closely related to D. transversa. We present several lines of evidence that both incomplete lineage sorting and gene flow are causing phylogenetic discordance. We suggest that ongoing gene flow primarily from D. recens into D. subquinaria in the sympatric part of their ranges causes phylogenetic uncertainty in the evolutionary history of these species. Our results highlight how population genetic data can be used to disentangle the sources of phylogenetic discordance among closely related species.

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 Ka /Ks 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.

Fine-scale geographic patterns of gene flow and reproductive character displacement in Drosophila subquinaria and Drosophila recens.

When two species are incompletely isolated, strengthening premating isolation barriers in response to the production of low fitness hybrids may complete the speciation process. Here, we use the sister species Drosophila subquinaria and Drosophila recens to study the conditions under which this reinforcement of species boundaries occurs in natural populations. We first extend the region of known sympatry between these species, and then we conduct a fine-scale geographic survey of mate discrimination coupled with estimates of gene flow within and admixture between species. Within D. subquinaria, reinforcement is extremely effective: we find variation in mate discrimination both against D. recens males and against conspecific allopatric males on the scale of a few kilometres and in the face of gene flow both from conspecific populations and introgression from D. recens. In D. recens, we do not find evidence for increased mate discrimination in sympatry, even where D. recens is rare, consistent with substantial gene flow throughout the species' range. Finally, we find that introgression between species is asymmetric, with more from D. recens into D. subquinaria than vice versa. Within each species, admixture is highest in the geographic region where it is rare relative to the other species, suggesting that when hybrids are produced they are of low fitness. In sum, reinforcement within D. subquinaria is effective at maintaining species boundaries, but even when reinforcing selection is strong it may not always result in a pattern of strong reproductive character displacement due to variation in the frequency of hybridization and gene flow from neighbouring populations.

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.

An experimental test of the evolutionary consequences of sympatry in Drosophila subquinaria.

Prezygotic isolation is often stronger between sympatric as opposed to allopatric taxa, but the underlying cause can be difficult to infer from comparative studies alone. Experimental evolution, where evolutionary responses to treatments manipulating the presence/absence of heterospecific individuals are tracked, can provide a powerful complementary approach. We used experimental evolution to investigate a naturally occurring pattern of reproductive character displacement in the mushroom-feeding fly, Drosophila subquinaria. In nature, female D. subquinaria from populations sympatric with the closely related Drosophila recens discriminate more strongly against heterospecific males than do females from allopatric populations. Starting with 16 replicate allopatric populations of D. subquinaria, we manipulated the presence/absence of D. recens during mating (experimental sympatry vs. control) and, when present, we allowed hybrids to live or kill them each generation. Across 12 generations, heterospecific offspring production from no-choice mating trials between D. subquinaria females and D. recens males declined in both experimental sympatry treatments relative to the control, suggesting increased sexual isolation. Male cuticular hydrocarbon profiles also evolved, but only in the hybrids killed treatment. Our results strongly imply that the existing reproductive character displacement in wild D. subquinaria populations was an evolutionary response to selection arising from secondary contact with D. recens.

Association of polyandry and sex-ratio drive prevalence in natural populations of Drosophila neotestacea.

Selfish genetic elements bias their own transmission to the next generation, even at the expense of the fitness of their carrier. Sex-ratio (SR) meiotic drive occurs when an X-chromosome causes Y-bearing sperm to die during male spermatogenesis, so that it is passed on to all of the male's offspring, which are all daughters. How SR is maintained as a stable polymorphism in the absence of genetic suppressors of drive is unknown. Here, we investigate the potential for the female remating rate to affect SR dynamics in natural populations, using the fly Drosophila neotestacea. In controlled laboratory conditions, females from populations where SR is rare mate more often than females from populations where SR is common. Furthermore, only when males mate multiply does the average fertility of SR males relative to wild-type males decrease to a level that can prevent SR from spreading. Our results suggest that differences in the female mating rate among populations may contribute to SR dynamics in the wild, and thus also affect the outcome of this intragenomic conflict. In line with this, we also present evidence of a localized population crash due to SR that may have resulted from habitat fragmentation along with a reduced mating rate.

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.

Epicuticular compounds of Drosophila subquinaria and D. recens: identification, quantification, and their role in female mate choice.

The epicuticle of various Drosophila species consists of long-chain cuticular hydrocarbons (CHCs) and their derivatives that play a role in waterproofing and a dynamic means of chemical communication. Here, via gas chromatography and mass spectrometry, we identified and quantified the epicuticular composition of D. recens and D. subquinaria, two closely related species that show a pattern of reproductive character displacement in nature. Twenty-four compounds were identified with the most abundant, 11-cis-Vaccenyl acetate, present only in males of each species. Also exclusive to males were five tri-acylglycerides. The 18 remaining compounds were CHCs, all shared between the sexes and species. These CHCs were composed of odd carbon numbers (C29, C31, C33, and C35), with an increase in structural isomers in the C33 and C35 groups. Saturated hydrocarbons comprise only methyl-branched alkanes and were found only in the C29 and C31 groups. Alkenes were the least prevalent, with alkadienes dominating the chromatographic landscape in the longer chain lengths. Sexual dimorphism was extensive with 6/8 of the logcontrast CHCs differing significantly in relative concentration between males and females in D. recens and D. subquinaria, respectively. Males of the two species also differed significantly in relative concentration of six CHCs, while females differed in none. Female-choice mating trials revealed directional sexual selection on male CHCs in a population of each species, consistent with female mate preferences for these traits. The sexual selection vectors differed significantly in multivariate trait space, suggesting that different pheromone blends determine male attractiveness in each species.

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.

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

In contrast to the rest of the genome, the Y chromosome is restricted to males and lacks recombination. As a result, Y chromosomes are unable to respond efficiently to selection, and newly formed Y chromosomes degenerate until few genes remain. The rapid loss of genes from newly formed Y chromosomes has been well studied, but gene loss from highly degenerate Y chromosomes has only recently received attention. Here, we identify and characterize a Y to autosome duplication of the male fertility gene kl-5 that occurred during the evolution of the testacea group species of Drosophila. The duplication was likely DNA based, as other Y-linked genes remain on the Y chromosome, the locations of introns are conserved, and expression analyses suggest that regulatory elements remain linked. Genetic mapping reveals that the autosomal copy of kl-5 resides on the dot chromosome, a tiny autosome with strongly suppressed recombination. Molecular evolutionary analyses show that autosomal copies of kl-5 have reduced polymorphism and little recombination. Importantly, the rate of protein evolution of kl-5 has increased significantly in lineages where it is on the dot versus Y linked. Further analyses suggest this pattern is a consequence of relaxed purifying selection, rather than adaptive evolution. Thus, although the initial fixation of the kl-5 duplication may have been advantageous, slightly deleterious mutations have accumulated in the dot-linked copies of kl-5 faster than in the Y-linked copies. Because the dot chromosome contains seven times more genes than the Y and is exposed to selection in both males and females, these results suggest that the dot suffers the deleterious effects of genetic linkage to more selective targets compared with the Y chromosome. Thus, a highly degenerate Y chromosome may not be the worst environment in the genome, as is generally thought, but may in fact be protected from the accumulation of deleterious mutations relative to other nonrecombining regions that contain more genes.

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.

Asymmetrical reinforcement and Wolbachia infection in Drosophila.

Reinforcement refers to the evolution of increased mating discrimination against heterospecific individuals in zones of geographic overlap and can be considered a final stage in the speciation process. One the factors that may affect reinforcement is the degree to which hybrid matings result in the permanent loss of genes from a species' gene pool. Matings between females of Drosophila subquinaria and males of D. recens result in high levels of offspring mortality, due to interspecific cytoplasmic incompatibility caused by Wolbachia infection of D. recens. Such hybrid inviability is not manifested in matings between D. recens females and D. subquinaria males. Here we ask whether the asymmetrical hybrid inviability is associated with a corresponding asymmetry in the level of reinforcement. The geographic ranges of D. recens and D. subquinaria were found to overlap across a broad belt of boreal forest in central Canada. Females of D. subquinaria from the zone of sympatry exhibit much stronger levels of discrimination against males of D. recens than do females from allopatric populations. In contrast, such reproductive character displacement is not evident in D. recens, consistent with the expected effects of unidirectional cytoplasmic incompatibility. Furthermore, there is substantial behavioral isolation within D. subquinaria, because females from populations sympatric with D. recens discriminate against allopatric conspecific males, whereas females from populations allopatric with D. recens show no discrimination against any conspecific males. Patterns of general genetic differentiation among populations are not consistent with patterns of behavioral discrimination, which suggests that the behavioral isolation within D. subquinaria results from selection against mating with Wolbachia-infected D. recens. Interspecific cytoplasmic incompatibility may contribute not only to post-mating isolation, an effect already widely recognized, but also to reinforcement, particularly in the uninfected species. The resulting reproductive character displacement not only increases behavioral isolation from the Wolbachia-infected species, but may also lead to behavioral isolation between populations of the uninfected species. Given the widespread occurrence of Wolbachia among insects, it thus appears that there are multiple ways by which these endosymbionts may directly and indirectly contribute to reproductive isolation and speciation.

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.

The Ecology and Evolutionary Dynamics of Meiotic Drive.

Meiotic drivers are genetic variants that selfishly manipulate the production of gametes to increase their own rate of transmission, often to the detriment of the rest of the genome and the individual that carries them. This genomic conflict potentially occurs whenever a diploid organism produces a haploid stage, and can have profound evolutionary impacts on gametogenesis, fertility, individual behaviour, mating system, population survival, and reproductive isolation. Multiple research teams are developing artificial drive systems for pest control, utilising the transmission advantage of drive to alter or exterminate target species. Here, we review current knowledge of how natural drive systems function, how drivers spread through natural populations, and the factors that limit their invasion.

Evolutionary dynamics of a spatially structured host-parasite association: Drosophila innubila and male-killing Wolbachia.

The mode and tempo of host-parasite evolution depend on population structure and history and the strength of selection that the species exert on each other. Here we genetically and epidemiologically characterize populations of the mycophagous fly Drosophila innubila and its male-killing Wolbachia endosymbiont, with the aim of integrating the local through global nature of this association. Drosophila innubila inhabit the forested "sky island" regions of the of the southwestern United States and northern Mexico, where its distribution is highly fragmented. We examine geographically isolated sky island populations of D. innubila, surveying the frequency and expression of Wolbachia infection as well as the distribution of genetic variation within and among populations of the host and parasite. In all populations, Wolbachia infection is associated with virtually complete male-killing, thus providing no evidence for the evolution of population-specific interaction phenotypes or local resistance. Although Wolbachia infection occurs in each of the main populations, there is variation among populations in the prevalence of infection and the resulting population-level sex ratio of D. innubila. Among these populations, the nuclear genes of D. innubila show moderate, though significant, differentiation. In contrast, the host mitochondrial DNA (mtDNA), which shares transmission with Wolbachia, exhibits substantially greater geographic differentiation, even after accounting for differences in transmission between nuclear and mitochondrial genes. We suggest that this pattern is caused by local Wolbachia--but not D. innubila--fluctuations in prevalence that increase the severity of drift experienced only by the mtDNA. Overall, our data suggest that the association between D. innubila and male-killing Wolbachia is ecologically dynamic within local populations, but evolutionarily coherent across the species as a whole.

Expression and modulation of embryonic male-killing in Drosophila innubila: opportunities for multilevel selection.

Organisms and the symbionts they harbor may experience opposing forces of selection. In particular, the contrasting inheritance patterns of maternally transmitted symbionts and their host's nuclear genes can engender conflict among organizational levels over the optimal host offspring sex ratio. This study uses a male-killing Wolbachia endosymbiont and its host Drosophila innubila to experimentally address the potential for multilevel selection in a host-symbiont system. We show that bacterial density can vary among infected females, and that females with a higher density have a more female-biased offspring sex ratio. Furthermore, bacterial density is an epigenetic and heritable trait: females with a low bacterial load have daughters with a lower-than-average bacterial density, whose offspring then experience less severe male-killing. For infected sons, the probability of embryonic mortality increases with the bacterial density in their mothers. The frequency distribution of Wolbachia density among individual D. innubila females, and therefore the dynamics of infection within populations of these flies, results both from processes affecting the growth and regulation of bacterial populations within cytoplasmic lineages and from selection among cytoplasmic lineages that vary in bacterial density. Estimates of effective population size of Wolbachia within cytoplasmic lineages and of D. innubila at the host population level suggest that selection among cytoplasmic lineages is likely to overwhelm the results of selection within lineages.