Rapid Global Spread of wRi-like Wolbachia across Multiple Drosophila

Maternally transmitted Wolbachia, Spiroplasma, and Cardinium bacteria are common in insects [1], but their interspecific spread is poorly understood. Endosymbionts can spread rapidly within host species by manipulating host reproduction, as typified by the global spread of wRi Wolbachia observed in Drosophila simulans [2, 3]. However, because Wolbachia cannot survive outside host cells, spread between distantly related host species requires horizontal transfers that are presumably rare [4-7]. Here, we document spread of wRi-like Wolbachia among eight highly diverged Drosophila hosts (10-50 million years) over only about 14,000 years (5,000-27,000). Comparing 110 wRi-like genomes, we find ≤0.02% divergence from the wRi variant that spread rapidly through California populations of D. simulans. The hosts include both globally invasive species (D. simulans, D. suzukii, and D. ananassae) and narrowly distributed Australian endemics (D. anomalata and D. pandora) [8]. Phylogenetic analyses that include mtDNA genomes indicate introgressive transfer of wRi-like Wolbachia between closely related species D. ananassae, D. anomalata, and D. pandora but no horizontal transmission within species. Our analyses suggest D. ananassae as the Wolbachia source for the recent wRi invasion of D. simulans and D. suzukii as the source of Wolbachia in its sister species D. subpulchrella. Although six of these wRi-like variants cause strong cytoplasmic incompatibility, two cause no detectable reproductive effects, indicating that pervasive mutualistic effects [9, 10] complement the reproductive manipulations for which Wolbachia are best known. "Super spreader" variants like wRi may be particularly useful for controlling insect pests and vector-borne diseases with Wolbachia transinfections [11].

Persistence of a Wolbachia infection frequency cline in Drosophila melanogaster and the possible role of reproductive dormancy

Field populations of arthropods are often polymorphic for Wolbachia but the factors maintaining intermediate Wolbachia frequencies are generally not understood. In Drosophila melanogaster, Wolbachia frequencies are highly variable across the globe. We document the persistence of a Wolbachia infection frequency cline in D. melanogaster populations from eastern Australia across at least 20 years, with frequencies generally high in the tropics but lower in cool temperate regions. The results are interpreted using a model of frequency dynamics incorporating cytoplasmic incompatibility (CI), imperfect maternal transmission and Wolbachia effects on fitness. Clinal variation is less pronounced in eastern North America which may reflect annual recolonization at higher latitudes. Limited samples from Africa from latitudes matching our tropical and subtropical samples from Australia and North America show comparably high infection frequencies, but some equatorial samples show lower frequencies. Adult dormancy across cold periods may contribute to the Australian Wolbachia cline. Infected flies exposed to cold conditions for an extended period had reduced fecundity and viability, an effect not evident in unexposed controls. These fitness costs may contribute to the relatively low Wolbachia frequencies in Australian temperate areas; whereas different processes, including CI induced by young males, may contribute to higher frequencies in tropical locations.

Loss of cytoplasmic incompatibility and minimal fecundity effects explain relatively low Wolbachia frequencies in Drosophila mauritiana

Maternally transmitted Wolbachia bacteria infect about half of all insect species. Many Wolbachia cause cytoplasmic incompatibility (CI) and reduced egg hatch when uninfected females mate with infected males. Although CI produces a frequency-dependent fitness advantage that leads to high equilibrium Wolbachia frequencies, it does not aid Wolbachia spread from low frequencies. Indeed, the fitness advantages that produce initial Wolbachia spread and maintain non-CI Wolbachia remain elusive. wMau Wolbachia infecting Drosophila mauritiana do not cause CI, despite being very similar to CI-causing wNo from Drosophila simulans (0.068% sequence divergence over 682,494 bp), suggesting recent CI loss. Using draft wMau genomes, we identify a deletion in a CI-associated gene, consistent with theory predicting that selection within host lineages does not act to increase or maintain CI. In the laboratory, wMau shows near-perfect maternal transmission; but we find no significant effect on host fecundity, in contrast to published data. Intermediate wMau frequencies on the island of Mauritius are consistent with a balance between unidentified small, positive fitness effects and imperfect maternal transmission. Our phylogenomic analyses suggest that group-B Wolbachia, including wMau and wPip, diverged from group-A Wolbachia, such as wMel and wRi, 6-46 million years ago, more recently than previously estimated.

Wolbachia Acquisition by Drosophila yakuba-Clade Hosts and Transfer of Incompatibility Loci Between Distantly Related Wolbachia

Maternally transmitted Wolbachia infect about half of insect species, yet the predominant mode(s) of Wolbachia acquisition remains uncertain. Species-specific associations could be old, with Wolbachia and hosts codiversifying (i.e., cladogenic acquisition), or relatively young and acquired by horizontal transfer or introgression. The three Drosophila yakuba-clade hosts [(D. santomea, D. yakuba) D. teissieri] diverged ∼3 MYA and currently hybridize on the West African islands Bioko and São Tomé. Each species is polymorphic for nearly identical Wolbachia that cause weak cytoplasmic incompatibility (CI)-reduced egg hatch when uninfected females mate with infected males. D. yakuba-clade Wolbachia are closely related to wMel, globally polymorphic in D. melanogaster We use draft Wolbachia and mitochondrial genomes to demonstrate that D. yakuba-clade phylogenies for Wolbachia and mitochondria tend to follow host nuclear phylogenies. However, roughly half of D. santomea individuals, sampled both inside and outside of the São Tomé hybrid zone, have introgressed D. yakuba mitochondria. Both mitochondria and Wolbachia possess far more recent common ancestors than the bulk of the host nuclear genomes, precluding cladogenic Wolbachia acquisition. General concordance of Wolbachia and mitochondrial phylogenies suggests that horizontal transmission is rare, but varying relative rates of molecular divergence complicate chronogram-based statistical tests. Loci that cause CI in wMel are disrupted in D. yakuba-clade Wolbachia; but a second set of loci predicted to cause CI are located in the same WO prophage region. These alternative CI loci seem to have been acquired horizontally from distantly related Wolbachia, with transfer mediated by flanking Wolbachia-specific ISWpi1 transposons.

Genome comparisons indicate recent transfer of wRi-like Wolbachia between sister species Drosophila suzukii and D. subpulchrella

Wolbachia endosymbionts may be acquired by horizontal transfer, by introgression through hybridization between closely related species, or by cladogenic retention during speciation. All three modes of acquisition have been demonstrated, but their relative frequency is largely unknown. Drosophila suzukii and its sister species D. subpulchrella harbor Wolbachia, denoted wSuz and wSpc, very closely related to wRi, identified in California populations of D. simulans. However, these variants differ in their induced phenotypes: wRi causes significant cytoplasmic incompatibility (CI) in D. simulans, but CI has not been detected in D. suzukii or D. subpulchrella. Our draft genomes of wSuz and wSpc contain full-length copies of 703 of the 734 single-copy genes found in wRi. Over these coding sequences, wSuz and wSpc differ by only 0.004% (i.e., 28 of 704,883 bp); they are sisters relative to wRi, from which each differs by 0.014%-0.015%. Using published data from D. melanogaster, Nasonia wasps and Nomada bees to calibrate relative rates of Wolbachia versus host nuclear divergence, we conclude that wSuz and wSpc are too similar-by at least a factor of 100-to be plausible candidates for cladogenic transmission. These three wRi-like Wolbachia, which differ in CI phenotype in their native hosts, have different numbers of orthologs of genes postulated to contribute to CI; and the CI loci differ at several nucleotides that may account for the CI difference. We discuss the general problem of distinguishing alternative modes of Wolbachia acquisition, focusing on the difficulties posed by limited knowledge of variation in absolute and relative rates of molecular evolution for host nuclear genomes, mitochondria, and Wolbachia.

Genome-wide divergence among invasive populations of Aedes aegypti in California

BACKGROUND

In the summer of 2013, Aedes aegypti Linnaeus was first detected in three cities in central California (Clovis, Madera and Menlo Park). It has now been detected in multiple locations in central and southern CA as far south as San Diego and Imperial Counties. A number of published reports suggest that CA populations have been established from multiple independent introductions.

RESULTS

Here we report the first population genomics analyses of Ae. aegypti based on individual, field collected whole genome sequences. We analyzed 46 Ae. aegypti genomes to establish genetic relationships among populations from sites in California, Florida and South Africa. Based on 4.65 million high quality biallelic SNPs, we identified 3 major genetic clusters within California; one that includes all sample sites in the southern part of the state (South of Tehachapi mountain range) plus the town of Exeter in central California and two additional clusters in central California.

CONCLUSIONS

A lack of concordance between mitochondrial and nuclear genealogies suggests that the three founding populations were polymorphic for two main mitochondrial haplotypes prior to being introduced to California. One of these has been lost in the Clovis populations, possibly by a founder effect. Genome-wide comparisons indicate extensive differentiation between genetic clusters. Our observations support recent introductions of Ae. aegypti into California from multiple, genetically diverged source populations. Our data reveal signs of hybridization among diverged populations within CA. Genetic markers identified in this study will be of great value in pursuing classical population genetic studies which require larger sample sizes.

A phylogeny for the Drosophila montium species group: A model clade for comparative analyses

The Drosophila montium species group is a clade of 94 named species, closely related to the model species D. melanogaster. The montium species group is distributed over a broad geographic range throughout Asia, Africa, and Australasia. Species of this group possess a wide range of morphologies, mating behaviors, and endosymbiont associations, making this clade useful for comparative analyses. We use genomic data from 42 available species to estimate the phylogeny and relative divergence times within the montium species group, and its relative divergence time from D. melanogaster. To assess the robustness of our phylogenetic inferences, we use 3 non-overlapping sets of 20 single-copy coding sequences and analyze all 60 genes with both Bayesian and maximum likelihood methods. Our analyses support monophyly of the group. Apart from the uncertain placement of a single species, D. baimaii, our analyses also support the monophyly of all seven subgroups proposed within the montium group. Our phylograms and relative chronograms provide a highly resolved species tree, with discordance restricted to estimates of relatively short branches deep in the tree. In contrast, age estimates for the montium crown group, relative to its divergence from D. melanogaster, depend critically on prior assumptions concerning variation in rates of molecular evolution across branches, and hence have not been reliably determined. We discuss methodological issues that limit phylogenetic resolution - even when complete genome sequences are available - as well as the utility of the current phylogeny for understanding the evolutionary and biogeographic history of this clade.

Population genomics of Drosophila suzukii reveal longitudinal population structure and signals of migrations in and out of the continental United States

Drosophila suzukii, or spotted-wing drosophila, is now an established pest in many parts of the world, causing significant damage to numerous fruit crop industries. Native to East Asia, D. suzukii infestations started in the United States a decade ago, occupying a wide range of climates. To better understand invasion ecology of this pest, knowledge of past migration events, population structure, and genetic diversity is needed. In this study, we sequenced whole genomes of 237 individual flies collected across the continental United States, as well as several sites in Europe, Brazil, and Asia, to identify and analyze hundreds of thousands of genetic markers. We observed strong population structure between Western and Eastern US populations, but no evidence of any population structure between different latitudes within the continental United States, suggesting that there are no broad-scale adaptations occurring in response to differences in winter climates. We detect admixture from Hawaii to the Western United States and from the Eastern United States to Europe, in agreement with previously identified introduction routes inferred from microsatellite analysis. We also detect potential signals of admixture from the Western United States back to Asia, which could have important implications for shipping and quarantine policies for exported agriculture. We anticipate this large genomic dataset will spur future research into the genomic adaptations underlying D. suzukii pest activity and development of novel control methods for this agricultural pest.

A male-killing Wolbachia endosymbiont is concealed by another endosymbiont and a nuclear suppressor

Bacteria that live inside the cells of insect hosts (endosymbionts) can alter the reproduction of their hosts, including the killing of male offspring (male killing, MK). MK has only been described in a few insects, but this may reflect challenges in detecting MK rather than its rarity. Here, we identify MK Wolbachia at a low frequency (around 4%) in natural populations of Drosophila pseudotakahashii. MK Wolbachia had a stable density and maternal transmission during laboratory culture, but the MK phenotype which manifested mainly at the larval stage was lost rapidly. MK Wolbachia occurred alongside a second Wolbachia strain expressing a different reproductive manipulation, cytoplasmic incompatibility (CI). A genomic analysis highlighted Wolbachia regions diverged between the 2 strains involving 17 genes, and homologs of the wmk and cif genes implicated in MK and CI were identified in the Wolbachia assembly. Doubly infected males induced CI with uninfected females but not females singly infected with CI-causing Wolbachia. A rapidly spreading dominant nuclear suppressor genetic element affecting MK was identified through backcrossing and subsequent analysis with ddRAD SNPs of the D. pseudotakahashii genome. These findings highlight the complexity of nuclear and microbial components affecting MK endosymbiont detection and dynamics in populations and the challenges of making connections between endosymbionts and the host phenotypes affected by them.

Wolbachia in the spittlebug Prosapia ignipectus: Variable infection frequencies, but no apparent effect on host reproductive isolation

Animals serve as hosts for complex communities of microorganisms, including endosymbionts that live inside their cells. Wolbachia bacteria are perhaps the most common endosymbionts, manipulating host reproduction to propagate. Many Wolbachia cause cytoplasmic incompatibility (CI), which results in reduced egg hatch when uninfected females mate with infected males. Wolbachia that cause intense CI spread to high and relatively stable frequencies, while strains that cause weak or no CI tend to persist at intermediate, often variable, frequencies. Wolbachia could also contribute to host reproductive isolation (RI), although current support for such contributions is limited to a few systems. To test for Wolbachia frequency variation and effects on host RI, we sampled several local Prosapia ignipectus (Fitch) (Hemiptera: Cercopidae) spittlebug populations in the northeastern United States over two years, including closely juxtaposed Maine populations with different monomorphic color forms, "black" and "lined." We discovered a group-B Wolbachia (wPig) infecting P. ignipectus that diverged from group-A Wolbachia-like model wMel and wRi strains in Drosophila-6 to 46 MYA. Populations of the sister species Prosapia bicincta (Say) from Hawaii and Florida are uninfected, suggesting that P. ignipectus acquired wPig after their initial divergence. wPig frequencies were generally high and variable among sites and between years. While phenotyping wPig effects on host reproduction is not currently feasible, the wPig genome contains three divergent sets of CI loci, consistent with high wPig frequencies. Finally, Maine monomorphic black and monomorphic lined populations of P. ignipectus share both wPig and mtDNA haplotypes, implying no apparent effect of wPig on the maintenance of this morphological contact zone. We hypothesize P. ignipectus acquired wPig horizontally as observed for many Drosophila species, and that significant CI and variable transmission produce high but variable wPig frequencies.

cifB-transcript levels largely explain cytoplasmic incompatibility variation across divergent Wolbachia

Divergent hosts often associate with intracellular microbes that influence their fitness. Maternally transmitted Wolbachia bacteria are the most common of these endosymbionts, due largely to cytoplasmic incompatibility (CI) that kills uninfected embryos fertilized by Wolbachia-infected males. Closely related infections in females rescue CI, providing a relative fitness advantage that drives Wolbachia to high frequencies. One prophage-associated gene (cifA) governs rescue, and two contribute to CI (cifA and cifB), but CI strength ranges from very strong to very weak for unknown reasons. Here, we investigate CI-strength variation and its mechanistic underpinnings in a phylogenetic context across 20 million years (MY) of Wolbachia evolution in Drosophila hosts diverged up to 50 MY. These Wolbachia encode diverse Cif proteins (100% to 7.4% pairwise similarity), and AlphaFold structural analyses suggest that CifB sequence similarities do not predict structural similarities. We demonstrate that cifB-transcript levels in testes explain CI strength across all but two focal systems. Despite phylogenetic discordance among cifs and the bulk of the Wolbachia genome, closely related Wolbachia tend to cause similar CI strengths and transcribe cifB at similar levels. This indicates that other non-cif regions of the Wolbachia genome modulate cif-transcript levels. CI strength also increases with the length of the host's larval life stage, presumably due to prolonged cif action. Our findings reveal that cifB-transcript levels largely explain CI strength, while highlighting other covariates. Elucidating CI's mechanism contributes to our understanding of Wolbachia spread in natural systems and to improving the efficacy of CI-based biocontrol of arboviruses and agricultural pests globally.

Ancient origin of an urban underground mosquito

Understanding how life is adapting to urban environments represents an important challenge in evolutionary biology. Here we investigate a widely cited example of urban adaptation, Culex pipiens form molestus, also known as the London Underground Mosquito. Population genomic analysis of ~350 contemporary and historical samples counter the popular hypothesis that molestus originated belowground in London less than 200 years ago. Instead, we show that molestus first adapted to human environments aboveground in the Middle East over the course of >1000 years, likely in concert with the rise of agricultural civilizations. Our results highlight the role of early human society in priming taxa for contemporary urban evolution and have important implications for understanding arbovirus transmission.

Rapid host switching of Wolbachia and even more rapid turnover of their phages and incompatibility-causing loci

About half of all insect species carry maternally inherited Wolbachia alphaproteobacteria, making Wolbachia the most common endosymbionts known in nature. Often Wolbachia spread to high frequencies within populations due to cytoplasmic incompatibility (CI), a Wolbachia-induced sperm modification caused by prophage-associated genes (cifs) that kill embryos without Wolbachia. Several Wolbachia variants also block viruses, including wMel from Drosophila melanogaster when transinfected into the mosquito Aedes aegypti. CI enables the establishment and stable maintenance of pathogen-blocking wMel in natural Ae. aegypti populations. These transinfections are reducing dengue disease incidence on multiple continents. While it has long been known that closely related Wolbachia occupy distantly related hosts, the timing of Wolbachia host switching and molecular evolution has not been widely quantified. We provide a new, conservative calibration for Wolbachia chronograms based on examples of co-divergence of Wolbachia and their insect hosts. Synthesizing publicly available and new genomic data, we use our calibration to demonstrate that wMel-like variants separated by only about 370,000 years have naturally colonized holometabolous dipteran and hymenopteran insects that diverged approximately 350 million years ago. Data from Wolbachia variants closely related to those currently dominant in D. melanogaster and D. simulans illustrate that cifs are rapidly acquired and lost among Wolbachia genomes, on a time scale of 104-105 years. This turnover occurs with and without the Wovirus prophages that contain them, with closely related cifs found in distantly related phages and distantly related cifs found in closely related phages. We present evidence for purifying selection on CI rescue function and on particular Cif protein domains. Our results quantify the tempo and mode of rapid host switching and horizontal gene transfer that underlie the spread and diversity of Wolbachia sampled from diverse host species. The wMel variants we highlight from hosts in different climates may offer new options for broadening Wolbachia-based biocontrol of diseases and pests.

Counting rare Wolbachia endosymbionts using digital droplet PCR

Wolbachia is the most widespread animal-associated intracellular microbe, living within the cells of over half of insect species. Since they can suppress pathogen replication and spread rapidly through insect populations, Wolbachia is at the vanguard of public health initiatives to control mosquito-borne diseases. Wolbachia's abilities to block pathogens and spread quickly are closely linked to their abundance in host tissues. The most common method for counting Wolbachia is quantitative polymerase chain reaction (qPCR), yet qPCR can be insufficient to count rare Wolbachia, necessitating tissue pooling and consequently compromising individual-level resolution of Wolbachia dynamics. Digital droplet PCR (ddPCR) offers superior sensitivity, enabling the detection of rare targets and eliminating the need for sample pooling. Here, we report three ddPCR assays to measure total Wolbachia abundance, Wolbachia abundance adjusted for DNA extraction efficiency, and Wolbachia density relative to host genome copies. Using Drosophila melanogaster with wMel Wolbachia as a model, we show these ddPCR assays can reliably detect as few as 7 to 12 Wolbachia gene copies in a 20 μL reaction. The designed oligos are homologous to sequences from at least 106 Wolbachia strains across Supergroup A and 53 host species from the Drosophila, Scaptomyza, and Zaprionus genera, suggesting broad utility. These highly sensitive ddPCR assays are expected to significantly advance Wolbachia-host interactions research by enabling the collection of molecular data from individual insect tissues. Their ability to detect rare Wolbachia will be especially valuable in applied and natural field settings where pooling samples could obscure important variation.

Fexinidazole and Corallopyronin A target Wolbachia-infected sheath cells present in filarial nematodes

The discovery of the endosymbiotic bacteria Wolbachia as an obligate symbiont of filarial nematodes has led to antibiotic-based treatments for filarial diseases. While lab and clinical studies have yielded promising results, recent animal studies reveal that Wolbachia levels may rebound following treatment with suboptimal doses of the antibiotic rifampicin. Previous work showed that a likely source of the bacterial rebound in females were dense clusters of Wolbachia in ovarian tissue. The number, size, and density of these Wolbachia clusters were not diminished despite antibiotic treatment. Here we define the cellular characteristics of the Wolbachia clusters in Brugia pahangi (wBp) and identify drugs that also target them. We have evidence that the Wolbachia clusters originate from newly formed sheath cells adjacent to the ovarian Distal Tip Cells. The dramatically enlarged volume of an infected sheath cell is strikingly similar to endosymbiont-induced bacteriocytes found in many insect species. Ultrastructural analysis reveals that the clustered Wolbachia present within the sheath cells exhibit a distinct morphology and form direct connections with the oocyte membrane and possibly the cytoplasm. This includes membrane-based channels providing a connection between Wolbachia-infected sheath cells and oocytes. We also determined that the Wolbachia within the sheath cells are either quiescent or replicating at a very low rate. Screens of known antibiotics and other drugs revealed that two drugs, Fexinidazole and Corallopyronin A, significantly reduced the number of clustered Wolbachia located within the sheath cells.