Environmental effects on cytoplasmic incompatibility and bacterial load in Wolbachia-infected Drosophila simulans

The effects of high temperatures, antibiotics, nutrition and larval density on cytoplasmic incompatibility caused by a Wolbachia infection were investigated in Drosophila simulans. Exposure of larvae from an infected stock to moderate doses of tetracycline led to complete incompatibility when treated females were crossed to infected males; the same doses only caused a partial restoration of compatibility when treated males were crossed to uninfected females. In crosses with treated females, there was a strong correlation between dose effects on hatch rates and infection levels in embryos produced by these females. Ageing and rearing males at a high temperature led to increased compatibility. However, exposing infected females to a high temperature did not influence their compatibility with infected males. Male temperature effects depended on conditions experienced at the larval stage but not the pupal stage. Exposure to 25 °C reduced the density of Wolbachia in embryos compared with a 19 °C treatment. Low levels of nutrition led to increased compatibility, but no effect of larval crowding was detected. These findings show the ways environmental factors can influence the expression of cytoplasmic incompatibility and suggest that environmental effects may be mediated by bacterial density.     

Contrasting patterns of geographic variation in the cosmopolitan sibling species Drosophila melanogaster and Drosophila simulans

An electrophoretic study was carried out to compare the geographic pattern of genetic variation in Drosophila simulans with that of its sibling species, Drosophila melanogaster. An identical set of 32 gene-protein loci was studied in four geographically distant populations of D. simulans and two populations of D. melanogaster, all originating from Europe and Africa. The comparison yielded the following results: (1) tropical populations of D. simulans were, in terms of the number of unique alleles, average heterozygosity per locus, and percentage of loci polymorphic, more variable than conspecific-temperate populations; (2) some loci in both species showed interpopulation differences in allele frequencies that suggest latitudinal clines; and (3) temperate-tropical genetic differentiation between populations was much less in D. simulans than in D. melanogaster. Similar differences between these two species have previously been shown for chromosomal, quantitative, physiological, and middle-repetitive DNA variation. Estimates of N _m (number of migrants per generation) from the spatial distribution of rare alleles suggest that both species have similar levels of interpopulation gene flow. These observations lead us to propose two competing hypotheses: the low level of geographic differentiation in D. simulans is due to its evolutionarily recent worldwide colonization and, alternatively, D. simulans has a narrower niche than D. melanogaster. Geographic variation data on different genetic elements (e.g., mitochondrial DNA, two-dimensional proteins, etc.) are required before these hypotheses can be adequately tested.We thank the Natural Science and Engineering Research Council of Canada for financial support (Grant A0235 to R.S.S.).     

Mitochondrial genotype affects fitness in Drosophila simulans

Drosophila simulans is known to harbor three distinct mitochondrial DNA (mtDNA) haplotype groups (siI, -II, and -III) with nearly 3.0% interhaplotypic divergence but <0.06% intrahaplotypic diversity. With the large amount of genetic variation in this system, the potential power to detect intraspecific fitness differences in fly lines that carry distinct haplotypes is great. We test three life-history traits on fly lines with known sequence differences in the mtDNA genome after controlling the nuclear genome by backcrossing. We find that flies with the siI haplotype are fastest developing and have the lowest probability of surviving to three experimental periods (2-6, 12-17, and 34-39 days of age). Wild-type males with siIII mtDNA were more active while disruption of specific coadapted nucleo-mitochondrial complexes caused a significant decrease in activity. These results are discussed in the context of the geographic distribution of each haplotype.     

A re-examination of Wolbachia-induced cytoplasmic incompatibility in California Drosophila simulans

Background

In California Drosophila simulans, the maternally inherited Riverside strain Wolbachia infection (wRi) provides a paradigm for rapid spread of Wolbachia in nature and rapid evolutionary change. wRi induces cytoplasmic incompatibility (CI), where crosses between infected males and uninfected females produce reduced egg-hatch. The three parameters governing wRi infection-frequency dynamics quantify: the fidelity of maternal transmission, the level of cytoplasmic incompatibility, and the relative fecundity of infected females. We last estimated these parameters in nature in 1993. Here we provide new estimates, under both field and laboratory conditions. Five years ago, we found that wRi had apparently evolved over 15 years to enhance the fecundity of infected females; here we examine whether CI intensity has also evolved.

Methodology/Principal Findings

New estimates using wild-caught flies indicate that the three key parameters have remained relatively stable since the early 1990s. As predicted by our three-parameter model using field-estimated parameter values, population infection frequencies remain about 93%. Despite this relative stability, laboratory data based on reciprocal crosses and introgression suggest that wRi may have evolved to produce less intense CI (i.e., higher egg hatch from incompatible crosses). In contrast, we find no evidence that D. simulans has evolved to lower the susceptibility of uninfected females to CI.

Conclusions/Significance

Evolution of wRi that reduces CI is consistent with counterintuitive theoretical predictions that within-population selection on CI-causing Wolbachia does not act to increase CI. Within taxa, CI is likely to evolve mainly via pleiotropic effects associated with the primary targets of selection on Wolbachia, i.e., host fecundity and transmission fidelity. Despite continuous, strong selection, D. simulans has not evolved appreciably to suppress CI. Our data demonstrate a lack of standing genetic variation for CI resistance in the host.     

Inbreeding alters intersexual fitness correlations in Drosophila simulans.

Intralocus sexual conflict results from sexually antagonistic selection on traits shared by the sexes. This can displace males and females from their respective fitness optima, and negative intersexual correlations (r_mf) for fitness are the unequivocal indicator of this evolutionary conflict. It has recently been suggested that intersexual fitness correlations can vary depending on the segregating genetic variation present in a population, and one way to alter genetic variation and test this idea is via inbreeding. Here, we test whether intersexual correlations for fitness vary with inbreeding in Drosophila simulans isolines reared under homogenous conditions. We measured male and female fitness at different times following the establishment of isofemale lines and found that the sign of the association between the two measures varied with time after initial inbreeding. Our results are consistent with suggestions that the type of genetic variation segregating within a population can determine the extent of intralocus sexual conflict and also support the idea that sexually antagonistic alleles segregate for longer in populations than alleles with sexually concordant effects.     

Evolution of Wolbachia-induced cytoplasmic incompatibility in Drosophila simulans and D. sechellia

Abstract.— The intracellular bacterium Wolbachia invades arthropod host populations through various mechanisms, the most common of which being cytoplasmic incompatibility (CI). CI involves elevated embryo mortality when infected males mate with uninfected females or females infected with different, incompatible Wolbachia strains. The present study focuses on this phenomenon in two Drosophila species: D. simulans and D. sechellia . Drosophila simulans populations are infected by several Wolbachia strains, including w Ha and w No. Drosophila sechellia is infected by only two Wolbachia : w Sh and w Sn. In both Drosophila species, double infections with Wolbachia are found. As indicated by several molecular markers, w Ha is closely related to w Sh, and w No to w Sn. Furthermore, the double infections in the two host species are associated with closely related mitochondrial haplotypes, namely si I (associated with w Ha and w No in D. simulans ) and se (associated with w Sh and w Sn in D. sechellia ). To test the theoretical prediction that Wolbachia compatibility types can diverge rapidly, we injected w Sh and w Sn into D. simulans , to compare their CI properties to those of their sister strains w Ha and w No, respectively, in the same host genetic background. We found that within each pair of sister strains CI levels were similar and that sister strains were fully compatible. We conclude that the short period for which the Wolbachia sister strains have been evolving separated from each other was not sufficient for their CI properties to diverge significantly.     

Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.     

A Wolbachia-associated fitness benefit depends on genetic background in Drosophila simulans

The alpha-proteobacteria Wolbachia infect a number of insect species and influence host reproduction to favour the spread of infected females through a population. The fitness effect of this infection is important in understanding the spread and maintenance of Wolbachia within and among host populations. However, a full elucidation of fitness effect requires careful control of host genetic background. Here, I transferred a single clone of Wolbachia (the wHa strain) into three genetically distinct isofemale lines of the fly Drosophila simulans using microinjection methodology. These lines carried one of the three described mitochondrial haplogroups (siI, siII or siIII) and differ in nuclear genome as well. Population cage assays showed that wHa-infected siIII flies enjoyed a dramatic fitness benefit compared to uninfected siIII. In contrast, wHa did not affect the fitness of siI or siII flies. This study points to the importance of host-by-symbiont interaction terms that may play an important role in organismal-fitness.     

Nucleotide variation at the runt locus in Drosophila melanogaster and Drosophila simulans.

Intra- and interspecific nucleotide variation for the major developmental gene runt in Drosophila was studied in D. melanogaster and D. simulans. The 1.5-kb protein-coding region and the 0.4-kb intron of the runt gene were sequenced for 11 alleles in each species. The D. melanogaster alleles originated from east Africa. Estimated parameters of intraspecific variation in D. melanogaster (exons: theta = 0.018, pi = 0.018; intron: theta = 0.014, pi = 0.014) and D. simulans (exons: theta = 0.007, pi = 0.005; intron: theta = 0.008, pi = 0.005) were below average for other X-linked genes, while divergence between species (exons: D = 0.094; intron: D = 0.069) fell within the normal range for both silent and replacement changes. This estimate for runt, along with published values for three other genes in regions of normal recombination, show east African D. melanogaster to be roughly twice as polymorphic as D. simulans. The majority of nucleotide variation, silent and replacement, in both species was found to be selectively neutral using various statistical tests (HKA, McDonald-Kreitman, Tajima, and Fu and Li tests). Monte Carlo simulations of the coalescent process significantly rejected a Wright-Fisher model with respect to an amino acid polymorphism and the distribution of polymorphic sites among the D. simulans lines. This indicated an old lineage and may reflect ancestral population substructuring in D. simulans.     

Variable modes of inheritance of morphometrical traits in hybrids between Drosophila melanogaster and Drosophila simulans.

We investigated body-size inheritance in interspecific sterile hybrids by crossing a Drosophila simulans strain with 13 strains of Drosophila melanogaster, which were of various origins and chosen for their broad range of genetic variation. A highly significant parent-offspring correlation was observed, showing that the D. melanogaster genes for size are still expressed in a hybrid background. Superimposed on to this additive inheritance, the size of hybrids was always less than the mid-parent value. This phenomenon, which at first sight might be described as dominance or overdominance, is more precisely interpreted as a consequence of a hybrid breakdown, that is, a dysfunction of the parental genes for size when put to work together. This interpretation is enforced by the fact that phenotypic variability was much more prevalent in hybrids than in parents. We also analysed body pigmentation inheritance in the same crosses and got a very different picture. There was no increase in the phenotypic variance of F(1) hybrids and only a low parent-offspring correlation. Apparent overdominance could be observed but in opposite directions, with no evidence of hybrid breakdown. Our data point to the possibility of analysing a diversity of quantitative traits in interspecific hybrids, and indicate that breakdown might be restricted to some traits only.     

Contrasting patterns of X-linked and autosomal nucleotide variation in Drosophila melanogaster and Drosophila simulans

Surveys of molecular variation in Drosophila melanogaster and Drosophila simulans have suggested that diversity outside of Africa is a subset of that within Africa. It has been argued that reduced levels of diversity in non-African populations reflect a population bottleneck, adaptation to temperate climates, or both. Here, I summarize the available single-nucleotide polymorphism data for both species. A simple "out of Africa" bottleneck scenario is consistent with geographic patterns for loci on the X chromosome but not with loci on the autosomes. Interestingly, there is a trend toward lower nucleotide diversity on the X chromosome relative to autosomes in non-African populations of D. melanogaster, but the opposite trend is seen in African populations. In African populations, autosomal inversion polymorphisms in D. melanogaster may contribute to reduced autosome diversity relative to the X chromosome. To elucidate the role that selection might play in shaping patterns of variability, I present a summary of within- and between-species patterns of synonymous and replacement variation in both species. Overall, D. melanogaster autosomes harbor an excess of amino acid replacement polymorphisms relative to D. simulans. Interestingly, range expansion from Africa appears to have had little effect on synonymous-to-replacement polymorphism ratios.     

The frequency distribution of nucleotide variation in Drosophila simulans.

Patterns of codon bias in Drosophila suggest that silent mutations can be classified into two types: unpreferred (slightly deleterious) and preferred (slightly beneficial). Results of previous analyses of polymorphism and divergence in Drosophila simulans were interpreted as supporting a mutation-selection-drift model in which slightly deleterious, silent mutants make significantly greater contributions to polymorphism than to divergence. Frequencies of unpreferred polymorphisms were inferred to be lower than frequencies of other silent polymorphisms. Here, I analyzed additional D. simulans data to reevaluate the support for these ideas. I found that D. simulans has fixed more unpreferred than preferred mutations, suggesting that this lineage has not been at mutation-selection-drift equilibrium at silent sites. Frequencies of polarized unpreferred polymorphisms are not skewed toward rare alleles. However, frequencies of unpolarized unpreferred codons are lower in high-bias genes than in low-bias genes. This supports the idea that unpreferred codons are borderline deleterious mutations. Purifying selection on silent sites appears to be stronger at twofold-degenerate codons than at fourfold-degenerate codons. Finally, I found that X-linked polymorphisms occur at a higher average frequency than polymorphisms on chromosome arm 3R, even though an average X-linked site is significantly less likely to be polymorphic than an average site on 3R. This result supports a previous analysis of D. simulans indicating different population genetics of X-linked versus autosomal mutations.     

Nucleotide variation and recombination along the fourth chromosome in Drosophila simulans

The fourth chromosome of Drosophila melanogaster and its sister species are believed to be nonrecombining and have been a model system for testing predictions of the effects of selection on linked, neutral variation. We recently examined nucleotide variation along the chromosome of D. melanogaster and revealed that a low average level of recombination could be associated with considerably high levels of nucleotide variation. In this report, we further investigate the variation along the fourth chromosome of D. simulans. We sequenced 12 gene regions evenly distributed along the fourth chromosome for a worldwide collection of 11 isofemale lines and 5 gene regions in a local population of 10 isofemale lines from South America. In contrast to predictions for regions of very low recombination, these data reveal that the variation levels in many gene regions, including an intron region of the ci gene, vary considerably along the fourth chromosome. Nucleotide diversity ranged from 0.0010 to 0.0074 in 9 gene regions interspersed with several regions of greatly reduced variation. Tests of recombination indicate that the recombination level is not as low as previously thought, likely an order of magnitude higher than that in D. melanogaster. Finally, estimates of the recombination parameters are shown to support a crossover-plus-conversion model.     

Evolution and inheritance of early embryonic patterning in Drosophila simulans and D. sechellia

Pattern formation in Drosophila is a widely studied example of a robust developmental system. Such robust systems pose a challenge to adaptive evolution, as they mask variation that selection may otherwise act upon. Yet we find variation in the localization of expression domains (henceforth "stripe allometry") in the pattern formation pathway. Specifically, we characterize differences in the gap genes giant and Kruppel, and the pair-rule gene even-skipped, which differ between the sibling species Drosophila simulans and D. sechellia. In a double-backcross experiment, stripe allometry is consistent with maternal inheritance of stripe positioning and multiple genetic factors, with a distinct genetic basis from embryo length. Embryos produced by F1 and F2 backcross mothers exhibit novel spatial patterns of gene expression relative to the parental species, with no measurable increase in positional variance among individuals. Buffering of novel spatial patterns in the backcross genotypes suggests that robustness need not be disrupted in order for the trait to evolve, and perhaps the system is incapable of evolving to prevent the expression of all genetic variation. This limitation, and the ability of natural selection to act on minute genetic differences that are within the "margin of error" for the buffering mechanism, indicates that developmentally buffered traits can evolve without disruption of robustness.     

Genetic mapping of adaptive wing size variation in Drosophila simulans

Many ecologically important traits exhibit latitudinal variation. Body size clines have been described repeatedly in insects across multiple continents, suggesting that similar selective forces are shaping these geographical gradients. It is unknown whether these parallel clinal patterns are controlled by the same or different genetic mechanism(s). We present here, quantitative trait loci (QTL) analysis of wing size variation in Drosophila simulans. Our results show that much of the wing size variation is controlled by a QTL on Chr 3L with relatively minor contribution from other chromosome arms. Comparative analysis of the genomic positions of the QTL indicates that the major QTL on Chr 3 are distinct in D. simulans and D. melanogaster, whereas the QTL on Chr 2R might overlap between species. Our results suggest that parallel evolution of wing size clines could be driven by non-identical genetic mechanisms but in both cases involve a major QTL as well as smaller effects of other genomic regions.     

Genomic effects of nucleotide substitutions in Drosophila simulans

Selective fixation of beneficial mutations reduces levels of linked, neutral variation. The magnitude of this "hitchhiking effect" is determined by the strength of selection and the recombination rate between selected and neutral sites. Thus, depending on the values of these parameters and the frequency with which directional selection occurs, the genomic scale over which directional selection reduces levels of linked variation may vary widely. Here we present a permutation-based analysis of nucleotide polymorphisms and fixations in Drosophila simulans. We show evidence of pervasive small-scale hitchhiking effects in this lineage. Furthermore, our results reveal that different types of fixations are associated with different levels of linked variation.     

Characterization of non-cytoplasmic incompatibility inducing Wolbachia in two continental African populations of Drosophila simulans

Wolbachia is an endocellular bacterium infecting arthropods and nematodes. In arthropods, it invades host populations through various mechanisms, affecting host reproduction, the most common of which being cytoplasmic incompatibility (CI). CI is an embryonic mortality occurring when infected males mate with uninfected females or females infected by a different Wolbachia strain. This phenomenon is observed in Drosophila simulans, an intensively studied Wolbachia host, harbouring at least five distinct bacterial strains. In this study, we investigate various aspects of the Wolbachia infections occurring in two continental African populations of D. simulans: CI phenotype, phylogenetic position based on the wsp gene and associated mitochondrial haplotype. From the East African population (Tanzania), we show that (i) the siIII mitochondrial haplotype occurs in continental populations, which was unexpected based on the current views of D. simulans biogeography, (ii) the wKi strain (that rescues from CI while being unable to induce it) is very closely related to the CI-inducing strain wNo, (iii) wKi and wNo might not derive from a unique infection event, and (iv) wKi is likely to represent the same entity as the previously described wMa variant. In the West African population (Cameroon), the Wolbachia infection was found identical to the previously described wAu, which does not induce CI. This finding supports the view that wAu might be an ancient infection in D. simulans.