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.

Developing an appropriate evolutionary baseline model for the study of SARS-CoV-2 patient samples.

Over the past 3 years, Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spread through human populations in several waves, resulting in a global health crisis. In response, genomic surveillance efforts have proliferated in the hopes of tracking and anticipating the evolution of this virus, resulting in millions of patient isolates now being available in public databases. Yet, while there is a tremendous focus on identifying newly emerging adaptive viral variants, this quantification is far from trivial. Specifically, multiple co-occurring and interacting evolutionary processes are constantly in operation and must be jointly considered and modeled in order to perform accurate inference. We here outline critical individual components of such an evolutionary baseline model-mutation rates, recombination rates, the distribution of fitness effects, infection dynamics, and compartmentalization-and describe the current state of knowledge pertaining to the related parameters of each in SARS-CoV-2. We close with a series of recommendations for future clinical sampling, model construction, and statistical analysis.

A single mutation weakens symbiont-induced reproductive manipulation through reductions in deubiquitylation efficiency.

Animals interact with microbes that affect their performance and fitness, including endosymbionts that reside inside their cells. Maternally transmitted Wolbachia bacteria are the most common known endosymbionts, in large part because of their manipulation of host reproduction. For example, many Wolbachia cause cytoplasmic incompatibility (CI) that reduces host embryonic viability when Wolbachia-modified sperm fertilize uninfected eggs. Operons termed cifs control CI, and a single factor (cifA) rescues it, providing Wolbachia-infected females a fitness advantage. Despite CI's prevalence in nature, theory indicates that natural selection does not act to maintain CI, which varies widely in strength. Here, we investigate the genetic and functional basis of CI-strength variation observed among sister Wolbachia that infect Drosophila melanogaster subgroup hosts. We cloned, Sanger sequenced, and expressed cif repertoires from weak CI-causing wYak in Drosophila yakuba, revealing mutations suspected to weaken CI relative to model wMel in D. melanogaster A single valine-to-leucine mutation within the deubiquitylating (DUB) domain of the wYak cifB homolog (cidB) ablates a CI-like phenotype in yeast. The same mutation reduces both DUB efficiency in vitro and transgenic CI strength in the fly, each by about twofold. Our results map hypomorphic transgenic CI to reduced DUB activity and indicate that deubiquitylation is central to CI induction in cid systems. We also characterize effects of other genetic variation distinguishing wMel-like cifs Importantly, CI strength determines Wolbachia prevalence in natural systems and directly influences the efficacy of Wolbachia biocontrol strategies in transinfected mosquito systems. These approaches rely on strong CI to reduce human disease.

Patterns of Population Structure and Introgression Among Recently Differentiated Drosophila melanogaster Populations.

Despite a century of genetic analysis, the evolutionary processes that have generated the patterns of exceptional genetic and phenotypic variation in the model organism Drosophila melanogaster remains poorly understood. In particular, how genetic variation is partitioned within its putative ancestral range in Southern Africa remains unresolved. Here, we study patterns of population genetic structure, admixture, and the spatial structuring of candidate incompatibility alleles across a global sample, including 223 new accessions, predominantly from remote regions in Southern Africa. We identify nine major ancestries, six that primarily occur in Africa and one that has not been previously described. We find evidence for both contemporary and historical admixture between ancestries, with admixture rates varying both within and between continents. For example, while previous work has highlighted an admixture zone between broadly defined African and European ancestries in the Caribbean and southeastern USA, we identify West African ancestry as the most likely African contributor. Moreover, loci showing the strongest signal of introgression between West Africa and the Caribbean/southeastern USA include several genes relating to neurological development and male courtship behavior, in line with previous work showing shared mating behaviors between these regions. Finally, while we hypothesized that potential incompatibility loci may contribute to population genetic structure across the range of D. melanogaster; these loci are, on average, not highly differentiated between ancestries. This work contributes to our understanding of the evolutionary history of a key model system, and provides insight into the partitioning of diversity across its range.

Recurrent Collection of Drosophila melanogaster from Wild African Environments and Genomic Insights into Species History.

A long-standing enigma concerns the geographic and ecological origins of the intensively studied vinegar fly, Drosophila melanogaster. This globally distributed human commensal is thought to originate from sub-Saharan Africa, yet until recently, it had never been reported from undisturbed wilderness environments that could reflect its precommensal niche. Here, we document the collection of 288 D. melanogaster individuals from multiple African wilderness areas in Zambia, Zimbabwe, and Namibia. The presence of D. melanogaster in these remote woodland environments is consistent with an ancestral range in southern-central Africa, as opposed to equatorial regions. After sequencing the genomes of 17 wilderness-collected flies collected from Kafue National Park in Zambia, we found reduced genetic diversity relative to town populations, elevated chromosomal inversion frequencies, and strong differences at specific genes including known insecticide targets. Combining these genomes with existing data, we probed the history of this species' geographic expansion. Demographic estimates indicated that expansion from southern-central Africa began ∼10,000 years ago, with a Saharan crossing soon after, but expansion from the Middle East into Europe did not begin until roughly 1,400 years ago. This improved model of demographic history will provide an important resource for future evolutionary and genomic studies of this key model organism. Our findings add context to the history of D. melanogaster, while opening the door for future studies on the biological basis of adaptation to human environments.

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.

Pervasive effects of Wolbachia on host activity.

Heritable symbionts have diverse effects on the physiology, reproduction and fitness of their hosts. Maternally transmitted Wolbachia are one of the most common endosymbionts in nature, infecting about half of all insect species. We test the hypothesis that Wolbachia alter host behaviour by assessing the effects of 14 different Wolbachia strains on the locomotor activity of nine Drosophila host species. We find that Wolbachia alter the activity of six different host genotypes, including all hosts in our assay infected with wRi-like Wolbachia strains (wRi, wSuz and wAur), which have rapidly spread among Drosophila species in about the last 14 000 years. While Wolbachia effects on host activity were common, the direction of these effects varied unpredictably and sometimes depended on host sex. We hypothesize that the prominent effects of wRi-like Wolbachia may be explained by patterns of Wolbachia titre and localization within host somatic tissues, particularly in the central nervous system. Our findings support the view that Wolbachia have wide-ranging effects on host behaviour. The fitness consequences of these behavioural modifications are important for understanding the evolution of host-symbiont interactions, including how Wolbachia spread within host populations.

Revisiting classic clines in Drosophila melanogaster in the age of genomics.

Adaptation to spatially varying environments has been studied for decades, but advances in sequencing technology are now enabling researchers to investigate the landscape of genetic variation underlying this adaptation genome wide. In this review we highlight some of the decades-long research on local adaptation in Drosophila melanogaster from well-studied clines in North America and Australia. We explore the evidence for parallel adaptation and identify commonalities in the genes responding to clinal selection across continents as well as discussing instances where patterns differ among clines. We also investigate recent studies utilizing whole-genome data to identify clines in D. melanogaster and several other systems. Although connecting segregating genomic variation to variation in phenotypes and fitness remains challenging, clinal genomics is poised to increase our understanding of local adaptation and the selective pressures that drive the extensive phenotypic diversity observed in nature.

Aedes albopictus is present in the lowlands of southern Zambia.

Identifying the current geographic range of disease vectors is a critical first step towards determining effective mechanisms for controlling and potentially eradicating them. This is particularly true given that historical vector ranges may expand due to changing climates and human activity. The Aedes subgenus Stegomyia contains over 100 species, and among them, Ae. aegypti and Ae. albopictus mosquitoes represent the largest concern for public health, spreading dengue, chikungunya, and zika viruses. While Ae. aegypti has been observed in the country of Zambia for decades, Ae. albopictus has not. In 2015 we sampled four urban and three rural areas in Zambia for Aedes species. Using DNA barcoding, we confirmed the presence of immature and adult Ae. albopictus at two sites: Siavonga and Livingstone. These genotypes seem most closely related to specimens previously collected in Mozambique based on mtDNA barcoding. We resampled Siavonga and Livingstone sites in 2019, again observing immature and adult Ae. albopictus at both sites. Relative Ae. albopictus frequencies were similar between sites, with the exception of immature life stages, which were higher in Siavonga than in Livingstone in 2019. While Ae. albopictus frequencies did not vary through time in Livingstone, both immature and adult frequencies increased through time in Siavonga. This report serves to document the presence of Ae. albopictus in Zambia, which will contribute to understanding the potential public health implications of this disease vector in southern Africa.

Comparative analysis of Wolbachia maternal transmission and localization in host ovaries.

Many insects and other animals carry microbial endosymbionts that influence their reproduction and fitness. These relationships only persist if endosymbionts are reliably transmitted from one host generation to the next. Wolbachia are maternally transmitted endosymbionts found in most insect species, but transmission rates can vary across environments. Maternal transmission of wMel Wolbachia depends on temperature in natural Drosophila melanogaster hosts and in transinfected Aedes aegypti, where wMel is used to block pathogens that cause human disease. In D. melanogaster, wMel transmission declines in the cold as Wolbachia become less abundant in host ovaries and at the posterior pole plasm (the site of germline formation) in mature oocytes. Here, we assess how temperature affects maternal transmission and underlying patterns of Wolbachia localization across 10 Wolbachia strains diverged up to 50 million years-including strains closely related to wMel-and their natural Drosophila hosts. Many Wolbachia maintain high transmission rates across temperatures, despite highly variable (and sometimes low) levels of Wolbachia in the ovaries and at the developing germline in late-stage oocytes. Identifying strains like closely related wMel-like Wolbachia with stable transmission across variable environmental conditions may improve the efficacy of Wolbachia-based biocontrol efforts as they expand into globally diverse environments.

Comparative analysis of Wolbachia maternal transmission and localization in host ovaries.

Many insects and other animals carry microbial endosymbionts that influence their reproduction and fitness. These relationships only persist if endosymbionts are reliably transmitted from one host generation to the next. Wolbachia are maternally transmitted endosymbionts found in most insect species, but transmission rates can vary across environments. Maternal transmission of wMel Wolbachia depends on temperature in natural Drosophila melanogaster hosts and in transinfected Aedes aegypti, where wMel is used to block pathogens that cause human disease. In D. melanogaster, wMel transmission declines in the cold as Wolbachia become less abundant in host ovaries and at the posterior pole plasm (the site of germline formation) in mature oocytes. Here, we assess how temperature affects maternal transmission and underlying patterns of Wolbachia localization across 10 Wolbachia strains diverged up to 50 million years-including strains closely related to wMel-and their natural Drosophila hosts. Many Wolbachia maintain high transmission rates across temperatures, despite highly variable (and sometimes low) levels of Wolbachia in the ovaries and at the developing germline in late-stage oocytes. Identifying strains like closely related wMel-like Wolbachia with stable transmission across variable environmental conditions may improve the efficacy of Wolbachia-based biocontrol efforts as they expand into globally diverse environments.

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.

Flies developed small bodies and small cells in warm and in thermally fluctuating environments.

Although plasma membranes benefit cells by regulating the flux of materials to and from the environment, these membranes cost energy to maintain. Because smaller cells provide relatively more membrane area for transport, ectotherms that develop in warm environments should consist of small cells despite the energetic cost. Effects of constant temperatures on cell size qualitatively match this prediction, but effects of thermal fluctuations on cell size are unknown. Thermal fluctuations could favour either small or large cells; small cells facilitate transport during peaks in metabolic demand whereas large cells minimize the resources needed for homeoviscous adaptation. To explore this problem, we examined effects of thermal fluctuations during development on the size of epidermal cells in the wings of Drosophila melanogaster. Flies derived from a temperate population were raised at two mean temperatures (18 and 25°C), with either no variation or a daily variation of ±4°C. Flies developed faster at a mean temperature of 25°C. Thermal fluctuations sped development, but only at 18°C. An increase in the mean and variance of temperature caused flies to develop smaller cells and wings. Thermal fluctuations reduced the size of males at 18°C and the size of females at 25°C. The thorax, the wings and the cells decreased with an increase in the mean and in the variance of temperature, but the response of cells was the strongest. Based on this pattern, we hypothesize that development of the greater area of membranes under thermal fluctuations provides a metabolic advantage that outweighs the greater energetic cost of remodelling membranes.

Aedes albopictus is present in the lowlands of southern Zambia.

Identifying the current geographic range of disease vectors is a critical first step towards determining effective mechanisms for controlling and potentially eradicating them. This is particularly true given that historical vector ranges may expand due to changing climates and human activity. The Aedes subgenus Stegomyia contains over 100 species, and among them, Ae. aegypti and Ae. albopictus mosquitoes represent the largest concern for public health, spreading dengue, chikungunya, and Zika viruses. While Ae. aegypti has been observed in the country of Zambia for decades, Ae. albopictus has not. In 2015 we sampled four urban and two rural areas in Zambia for Aedes species. Using DNA barcoding, we confirmed the presence of immature and adult Ae. albopictus at two rural sites: Siavonga and Livingstone. These genotypes seem most closely related to specimens previously collected in Mozambique based on CO1 sequence from mtDNA. We resampled Siavonga and Livingstone sites in 2019, again observing immature and adult Ae. albopictus at both sites. Relative Ae. albopictus frequencies were similar between sites, with the exception of immature life stages, which were higher in Siavonga than in Livingstone in 2019. While Ae. albopictus frequencies did not vary through time in Livingstone, both immature and adult frequencies increased through time in Siavonga. This report serves to document the presence of Ae. albopictus in Zambia, which will contribute to the process of determining the potential public health implications of this disease vector in Central Africa.

A simulation framework for modeling the within-patient evolutionary dynamics of SARS-CoV-2.

The global impact of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has led to considerable interest in detecting novel beneficial mutations and other genomic changes that may signal the development of variants of concern (VOCs). The ability to accurately detect these changes within individual patient samples is important in enabling early detection of VOCs. Such genomic scans for positive selection are best performed via comparison of empirical data to simulated data wherein evolutionary factors, including mutation and recombination rates, reproductive and infection dynamics, and purifying and background selection, can be carefully accounted for and parameterized. While there has been work to quantify these factors in SARS-CoV-2, they have yet to be integrated into a baseline model describing intra-host evolutionary dynamics. To construct such a baseline model, we develop a simulation framework that enables one to establish expectations for underlying levels and patterns of patient-level variation. By varying eight key parameters, we evaluated 12,096 different model-parameter combinations and compared them to existing empirical data. Of these, 592 models (~5%) were plausible based on the resulting mean expected number of segregating variants. These plausible models shared several commonalities shedding light on intra-host SARS-CoV-2 evolutionary dynamics: severe infection bottlenecks, low levels of reproductive skew, and a distribution of fitness effects skewed towards strongly deleterious mutations. We also describe important areas of model uncertainty and highlight additional sequence data that may help to further refine a baseline model. This study lays the groundwork for the improved analysis of existing and future SARS-CoV-2 within-patient data.

A Simulation Framework for Modeling the Within-Patient Evolutionary Dynamics of SARS-CoV-2.

The global impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to considerable interest in detecting novel beneficial mutations and other genomic changes that may signal the development of variants of concern (VOCs). The ability to accurately detect these changes within individual patient samples is important in enabling early detection of VOCs. Such genomic scans for rarely acting positive selection are best performed via comparison of empirical data with simulated data wherein commonly acting evolutionary factors, including mutation and recombination, reproductive and infection dynamics, and purifying and background selection, can be carefully accounted for and parameterized. Although there has been work to quantify these factors in SARS-CoV-2, they have yet to be integrated into a baseline model describing intrahost evolutionary dynamics. To construct such a baseline model, we develop a simulation framework that enables one to establish expectations for underlying levels and patterns of patient-level variation. By varying eight key parameters, we evaluated 12,096 different model-parameter combinations and compared them with existing empirical data. Of these, 592 models (∼5%) were plausible based on the resulting mean expected number of segregating variants. These plausible models shared several commonalities shedding light on intrahost SARS-CoV-2 evolutionary dynamics: severe infection bottlenecks, low levels of reproductive skew, and a distribution of fitness effects skewed toward strongly deleterious mutations. We also describe important areas of model uncertainty and highlight additional sequence data that may help to further refine a baseline model. This study lays the groundwork for the improved analysis of existing and future SARS-CoV-2 within-patient data.

Distinct Wolbachia localization patterns in oocytes of diverse host species reveal multiple strategies of maternal transmission.

A broad array of endosymbionts radiate through host populations via vertical transmission, yet much remains unknown concerning the cellular basis, diversity, and routes underlying this transmission strategy. Here, we address these issues, by examining the cellular distributions of Wolbachia strains that diverged up to 50 million years ago in the oocytes of 18 divergent Drosophila species. This analysis revealed 3 Wolbachia distribution patterns: (1) a tight clustering at the posterior pole plasm (the site of germline formation); (2) a concentration at the posterior pole plasm, but with a significant bacteria population distributed throughout the oocyte; and (3) a distribution throughout the oocyte, with none or very few located at the posterior pole plasm. Examination of this latter class indicates Wolbachia accesses the posterior pole plasm during the interval between late oogenesis and the blastoderm formation. We also find that 1 Wolbachia strain in this class concentrates in the posterior somatic follicle cells that encompass the pole plasm of the developing oocyte. In contrast, strains in which Wolbachia concentrate at the posterior pole plasm generally exhibit no or few Wolbachia in the follicle cells associated with the pole plasm. Taken together, these studies suggest that for some Drosophila species, Wolbachia invade the germline from neighboring somatic follicle cells. Phylogenomic analysis indicates that closely related Wolbachia strains tend to exhibit similar patterns of posterior localization, suggesting that specific localization strategies are a function of Wolbachia-associated factors. Previous studies revealed that endosymbionts rely on 1 of 2 distinct routes of vertical transmission: continuous maintenance in the germline (germline-to-germline) or a more circuitous route via the soma (germline-to-soma-to-germline). Here, we provide compelling evidence that Wolbachia strains infecting Drosophila species maintain the diverse arrays of cellular mechanisms necessary for both of these distinct transmission routes. This characteristic may account for its ability to infect and spread globally through a vast range of host insect species.

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.

Temperature effects on cellular host-microbe interactions explain continent-wide endosymbiont prevalence.

Endosymbioses influence host physiology, reproduction, and fitness, but these relationships require efficient microbe transmission between host generations to persist. Maternally transmitted Wolbachia are the most common known endosymbionts,1 but their frequencies vary widely within and among host populations for unknown reasons.2,3 Here, we integrate genomic, cellular, and phenotypic analyses with mathematical models to provide an unexpectedly simple explanation for global wMel Wolbachia prevalence in Drosophila melanogaster. Cooling temperatures decrease wMel cellular abundance at a key stage of host oogenesis, producing temperature-dependent variation in maternal transmission that plausibly explains latitudinal clines of wMel frequencies on multiple continents. wMel sampled from a temperate climate targets the germline more efficiently in the cold than a recently differentiated tropical variant (∼2,200 years ago), indicative of rapid wMel adaptation to climate. Genomic analyses identify a very narrow list of wMel alleles-most notably, a derived stop codon in the major Wolbachia surface protein WspB-that underlie thermal sensitivity of cellular Wolbachia abundance and covary with temperature globally. Decoupling temperate wMel and host genomes further reduces transmission in the cold, a pattern that is characteristic of host-microbe co-adaptation to a temperate climate. Complex interactions among Wolbachia, hosts, and the environment (GxGxE) mediate wMel cellular abundance and maternal transmission, implicating temperature as a key determinant of Wolbachia spread and equilibrium frequencies, in conjunction with Wolbachia effects on host fitness and reproduction.4,5 Our results motivate the strategic use of locally selected wMel variants for Wolbachia-based biocontrol efforts, which protect millions of individuals from arboviruses that cause human disease.6.