Effects of wolbachia on genetic divergence between populations: mainland-island model

The Transmission Patterns of the Endosymbiont Wolbachia within the Hawaiian Drosophilidae Adaptive Radiation

The evolution of endosymbionts and their hosts can lead to highly dynamic interactions with varying fitness effects for both the endosymbiont and host species. Wolbachia, a ubiquitous endosymbiont of arthropods and nematodes, can have both beneficial and detrimental effects on host fitness. We documented the occurrence and patterns of transmission of Wolbachia within the Hawaiian Drosophilidae and examined the potential contributions of Wolbachia to the rapid diversification of their hosts. Screens for Wolbachia infections across a minimum of 140 species of Hawaiian Drosophila and Scaptomyza revealed species-level infections of 20.0%, and across all 399 samples, a general infection rate of 10.3%. Among the 44 Wolbachia strains we identified using a modified Wolbachia multi-locus strain typing scheme, 30 (68.18%) belonged to supergroup B, five (11.36%) belonged to supergroup A, and nine (20.45%) had alleles with conflicting supergroup assignments. Co-phylogenetic reconciliation analysis indicated that Wolbachia strain diversity within their endemic Hawaiian Drosophilidae hosts can be explained by vertical (e.g., co-speciation) and horizontal (e.g., host switch) modes of transmission. Results from stochastic character trait mapping suggest that horizontal transmission is associated with the preferred oviposition substrate of the host, but not the host's plant family or island of occurrence. For Hawaiian Drosophilid species of conservation concern, with 13 species listed as endangered and 1 listed as threatened, knowledge of Wolbachia strain types, infection status, and potential for superinfection could assist with conservation breeding programs designed to bolster population sizes, especially when wild populations are supplemented with laboratory-reared, translocated individuals. Future research aimed at improving the understanding of the mechanisms of Wolbachia transmission in nature, their impact on the host, and their role in host species formation may shed light on the influence of Wolbachia as an evolutionary driver, especially in Hawaiian ecosystems.

Cytoplasmic incompatibility in hybrid zones: infection dynamics and resistance evolution

Cytoplasmic incompatibility is an endosymbiont-induced mating incompatibility common in arthropods. Unidirectional cytoplasmic incompatibility impairs crosses between infected males and uninfected females, whereas bidirectional cytoplasmic incompatibility occurs when two host lineages are infected with reciprocally in compatible endosymbionts. Bidirectional cytoplasmic incompatibility is unstable in unstructured populations, but may be stable in hybrid zones. Stable coexistence of incompatible host lineages should generate frequent incompatible crosses. Therefore, hosts are expected to be under selection to resist their endosymbionts. Here, we for mulate a mathematical model of hybrid zones where two bidirectionally incompatible host lineages meet. We expand this model to consider the invasion of a hypothetical resistance allele. To corroborate our mathematical predictions, we test each prediction with stochastic, individual-based simulations. Our models suggest that hybrid zones may sustain stable coinfections of bidirectionally incompatible endosymbiont strains. Over a range of conditions, host are under selection for resistance against cytoplasmic incompatibility. Under asymetric migration, a resistance allele can facilitate infection turnover and subsequently either persist or become lost. The predictions we present may inform our understanding of the cophylogenetic relationship between the endosym biont Wolbachia and its hosts.

Genetic stability of Aedes aegypti populations following invasion by wMel Wolbachia

Background

Wolbachia wMel is the most commonly used strain in rear and release strategies for Aedes aegypti mosquitoes that aim to inhibit the transmission of arboviruses such as dengue, Zika, Chikungunya and yellow fever. However, the long-term establishment of wMel in natural Ae. aegypti populations raises concerns that interactions between Wolbachia wMel and Ae. aegypti may lead to changes in the host genome, which could affect useful attributes of Wolbachia that allow it to invade and suppress disease transmission.

Results

We applied an evolve-and-resequence approach to study genome-wide genetic changes in Ae. aegypti from the Cairns region, Australia, where Wolbachia wMel was first introduced more than 10 years ago. Mosquito samples were collected at three different time points in Gordonvale, Australia, covering the phase before (2010) and after (2013 and 2018) Wolbachia releases. An additional three locations where Wolbachia replacement happened at different times across the last decade were also sampled in 2018. We found that the genomes of mosquito populations mostly remained stable after Wolbachia release, with population differences tending to reflect the geographic location of the populations rather than Wolbachia infection status. However, outlier analysis suggests that Wolbachia may have had an influence on some genes related to immune response, development, recognition and behavior.

Conclusions

Ae. aegypti populations remained geographically distinct after Wolbachia wMel releases in North Australia despite their Wolbachia infection status. At some specific genomic loci, we found signs of selection associated with Wolbachia, suggesting potential evolutionary impacts can happen in the future and further monitoring is warranted.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08200-1.

Wolbachia and host intrinsic reproductive barriers contribute additively to postmating isolation in spider mites

Wolbachia are maternally-inherited bacteria that induce cytoplasmic incompatibility in many arthropod species. However, the ubiquity of this isolation mechanism for host speciation processes remains elusive, as only few studies have examined Wolbachia-induced incompatibilities when host populations are not genetically compatible. Here, we used three populations of two genetically differentiated colour forms of the haplodiploid spider mite Tetranychus urticae to dissect the interaction between Wolbachia-induced and host-associated incompatibilities, and their relative contribution to postmating isolation. We found that these two sources of incompatibility act through different mechanisms in an additive fashion. Host-associated incompatibility contributes 1.5 times more than Wolbachia-induced incompatibility in reducing hybrid production, the former through an overproduction of haploid sons at the expense of diploid daughters (ca. 75% decrease) and the latter by increasing the embryonic mortality of daughters (by ca. 49%). Furthermore, regardless of cross direction, we observed near-complete F1 hybrid sterility and complete F2 hybrid breakdown between populations of the two forms, but Wolbachia did not contribute to this outcome. We thus show mechanistic independence and an additive nature of host-intrinsic and Wolbachia-induced sources of isolation. Wolbachia may contribute to reproductive isolation in this system, thereby potentially affecting host differentiation and distribution in the field.

Geographic and Temporal Variation of Distinct Intracellular Endosymbiont Strains of Wolbachia sp. in the Grasshopper Chorthippus parallelus: a Frequency-Dependent Mechanism?

Wolbachia is an intracellular endosymbiont that can produce a range of effects on host fitness, but the temporal dynamics of Wolbachia strains have rarely been experimentally evaluated. We compare interannual strain frequencies along a geographical region for understanding the forces that shape Wolbachia strain frequency in natural populations of its host, Chorthippus parallelus (Orthoptera, Acrididae). General linear models show that strain frequency changes significantly across geographical and temporal scales. Computer simulation allows to reject the compatibility of the observed patterns with either genetic drift or sampling errors. We use consecutive years to estimate total Wolbachia strain fitness. Our estimation of Wolbachia fitness is significant in most cases, within locality and between consecutive years, following a negatively frequency-dependent trend. Wolbachia spp. B and F strains show a temporal pattern of variation that is compatible with a negative frequency-dependent natural selection mechanism. Our results suggest that such a mechanism should be at least considered in future experimental and theoretical research strategies that attempt to understand Wolbachia biodiversity.

Evidence for horizontal transmission of secondary endosymbionts in the Bemisia tabaci cryptic species complex

Bemisia tabaci (Hemiptera: Aleyrodidae) is a globally distributed pest composed of at least 34 morphologically indistinguishable cryptic species. At least seven species of endosymbiont have been found infecting some or all members of the complex. The origin(s) of the associations between specific endosymbionts and their whitefly hosts is unknown. Infection is normally vertical, but horizontal transmission does occur and is one way for new infections to be introduced into individuals. The relationships between the different members of the cryptic species complex and the endosymbionts have not been well explored. In this study, the phylogenies of different cryptic species of the host with those of their endosymbionts were compared. Of particular interest was whether there was evidence for both coevolution and horizontal transmission. Congruence was observed for the primary endosymbiont, Portiera aleyrodidarum, and partial incongruence in the case of two secondary endosymbionts, Arsenophonus and Cardinium and incongruence for a third, Wolbachia. The patterns observed for the primary endosymbiont supported cospeciation with the host while the patterns for the secondary endosymbionts, and especially Wolbachia showed evidence of host shifts and extinctions through horizontal transmission rather than cospeciation. Of particular note is the observation of several very recent host shift events in China between exotic invader and indigenous members of the complex. These shifts were from indigenous members of the complex to the invader as well as from the invader to indigenous relatives.

Maintenance of adaptive differentiation by Wolbachia induced bidirectional cytoplasmic incompatibility: the importance of sib-mating and genetic systems

BACKGROUND

Bacteria of the genus Wolbachia are reproductive parasites widespread among arthropods. The most common effect arising from the presence of Wolbachia in a population is Cytoplasmic Incompatibility (CI), whereby postmating reproductive isolation occurs in crosses between an infected male and an uninfected female, or when a male is infected with a different strain of Wolbachia to that of the female (bidirectional CI). Previous theoretical models have demonstrated that bidirectional CI can contribute to the genetic divergence of populations in haploid and diploid organisms. However, haplodiploid organisms were not considered in these models even though they include Nasonia parasitoid wasps - the best example of the implication of Wolbachia in ongoing speciation. Moreover, previous work did not investigate inbreeding mating systems, which are frequently observed in arthropod species.

RESULTS

We developed a stochastic two-island model which simulated three genetic scenarios, diploidy, haploidy, and haplodiploidy, with two CI phenotypes being considered for the latter: (1) male development of female progeny; and (2) mortality of fertilized eggs. We also investigated the effect of varying the proportion of sib mating. In the model each allopatric population was initially fixed for a single allele at a nuclear locus under positive selection and infected with one strain of Wolbachia. Each simulation presupposed that the two populations were fixed for a different allele and a different strain of Wolbachia. The degree of genetic differentiation observed in the locus under selection due to bidirectional CI was much lower for the two haplodiploid phenotypes than for either diploids or haploids. Furthermore, we demonstrated that sib-mating may compensate for the lower efficiency of bidirectional CI in haplodiploids by maintaining genetic divergence.

CONCLUSION

Our model suggests that maintenance of genetic differentiation facilitated by Wolbachia is more likely to occur in diploids and haploids than in haplodiploids. However, increasing the level of sib-mating may compensate for the weak effect of bidirectional CI in haplodiploids. Our work therefore gives a potential explanation for why the haplodiploid Nasonia species, which are infected with bidirectionally incompatible Wolbachia strains and undergo sib-mating, have differentiated genetically and maintained this differentiation without premating isolation.

The neutral effective migration rate in a mainland-island context

Genetic influx into a population often does not correspond to the real migration rate (m) of individuals, due to class structure within the population. The effective migration rate (m(e)) is a concept to measure gene flow in such a situation. The ratio of the effective migration rate to the real migration rate (m(e)/m) is called the gene flow factor, and represents the degree of gene flow modification. Prior authors proposed different definitions of the effective migration rate. These may be categorized into two groups: the neutral effective migration rate and the selective effective migration rate. In this article, we construct a general model of a class-structured population with a mainland-island structure. Using the model, we prove that the gene flow factor of the neutral effective migration rate converges to the mean reproductive value of immigrants if the limit is taken with the real migration rate converging to zero. This limit theorem provides a novel interpretation of gene flow and can be used to derive approximation formulae of the neutral effective migration rate. We illustrate this method analyzing two examples, sex ratio distortion due to extrinsic factors and hybrid zones with underdominance.

Wolbachia-induced unidirectional cytoplasmic incompatibility and speciation: mainland-island model

Bacteria of the genus Wolbachia are among the most common endosymbionts in the world. In many insect species these bacteria induce a sperm-egg incompatibility between the gametes of infected males and uninfected females, commonly called unidirectional cytoplasmic incompatibility (CI). It is generally believed that unidirectional CI cannot promote speciation in hosts because infection differences between populations will be unstable and subsequent gene flow will eliminate genetic differences between diverging populations. In the present study we investigate this question theoretically in a mainland-island model with migration from mainland to island. Our analysis shows that (a) the infection polymorphism is stable below a critical migration rate, (b) an (initially) uninfected “island” can better maintain divergence at a selected locus (e.g. can adapt locally) in the presence of CI, and (c) unidirectional CI selects for premating isolation in (initially) uninfected island populations if they receive migration from a Wolbachia-infected mainland. Interestingly, premating isolation is most likely to evolve if levels of incompatibility are intermediate and if either the infection causes fecundity reductions or Wolbachia transmission is incomplete. This is because under these circumstances an infection pattern with an infected mainland and a mostly uninfected island can persist in the face of comparably high migration. We present analytical results for all three findings: (a) a lower estimation of the critical migration rate in the presence of local adaptation, (b) an analytical approximation for the gene flow reduction caused by unidirectional CI, and (c) a heuristic formula describing the invasion success of mutants at a mate preference locus. These findings generally suggest that Wolbachia-induced unidirectional CI can be a factor in divergence and speciation of hosts.

Wolbachia-induced unidirectional cytoplasmic incompatibility and the stability of infection polymorphism in parapatric host populations

Wolbachia are intracellular, maternally inherited bacteria that are widespread among arthropods and commonly induce a reproductive incompatibility between infected male and uninfected female hosts known as unidirectional cytoplasmic incompatibility (CI). If infected and uninfected populations occur parapatrically, CI acts as a post-zygotic isolation barrier. We investigate the stability of such infection polymorphisms in a mathematical model with two populations linked by migration. We determine critical migration rates below which infected and uninfected populations can coexist. Analytical solutions of the critical migration rate are presented for mainland-island models. These serve as lower estimations for a more general model with two-way migration. The critical migration rate is positive if either Wolbachia causes a fecundity reduction in infected female hosts or its transmission is incomplete, and is highest for intermediate levels of CI. We discuss our results with respect to local adaptations of the Wolbachia host, speciation, and pest control.

Deep mtDNA divergences indicate cryptic species in a fig-pollinating wasp

BACKGROUND

Figs and fig-pollinating wasps are obligate mutualists that have coevolved for ca. 90 million years. They have radiated together, but do not show strict cospeciation. In particular, it is now clear that many fig species host two wasp species, so there is more wasp speciation than fig speciation. However, little is known about how fig wasps speciate.

RESULTS

We studied variation in 71 fig-pollinating wasps from across the large geographic range of Ficus rubiginosa in Australia. All wasps sampled belong to one morphological species (Pleistodontes imperialis), but we found four deep mtDNA clades that differed from each other by 9-17% nucleotides. As these genetic distances exceed those normally found within species and overlap those (10-26%) found between morphologically distinct Pleistodontes species, they strongly suggest cryptic fig wasp species. mtDNA clade diversity declines from all four present in Northern Queensland to just one in Sydney, near the southern range limit. However, at most sites multiple clades coexist and can be found in the same tree or even the same fig fruit and there is no evidence for parallel sub-division of the host fig species. Both mtDNA data and sequences from two nuclear genes support the monophyly of the "P. imperialis complex" relative to other Pleistodontes species, suggesting that fig wasp divergence has occurred without any host plant shift. Wasps in clade 3 were infected by a single strain (W1) of Wolbachia bacteria, while those in other clades carried a double infection (W2+W3) of two other strains.

CONCLUSION

Our study indicates that cryptic fig-pollinating wasp species have developed on a single host plant species, without the involvement of host plant shifts, or parallel host plant divergence. Despite extensive evidence for coevolution between figs and fig wasps, wasp speciation may not always be linked strongly with fig speciation.

Asymmetric gene flow and constraints on adaptation caused by sex ratio distorters

Asymmetric gene flow is generally believed to oppose natural selection and potentially impede adaptation. Whilst the cause of asymmetric gene flow has been seen largely in terms of variation in population density over space, asymmetric gene flow can also result from varying sex ratios across subpopulations with similar population sizes. We model the process of adaptation in a scenario in which two adjacent subpopulations have different sex ratios, associated with different levels of infection with maternally inherited endosymbionts that selectively kill male hosts. Two models are analyzed in detail. First, we consider one host locus with two alleles, each of which possesses a selective advantage in one of the subpopulations. We found that local adaptation can strongly be impeded in the subpopulation with the more female biased population sex ratio. Second, we analyze host alleles that provide resistance against the male-killing (MK) endosymbionts and show that asymmetric gene flow can prevent the spread of such alleles under certain conditions. These results might have important implications for the coevolution of MK bacteria and their hosts.

Bidirectional cytoplasmic incompatibility and the stable coexistence of two Wolbachia strains in parapatric host populations

Wolbachiaare intracellular bacteria which are very widely distributed among arthropods. In many insect species Wolbachiaare known to induce cytoplasmic mating incompatibility (CI). It has been suggested that Wolbachiacould promote speciation in their hosts if parapatric host populations are infected with two different Wolbachiastrains causing bidirectional mating incompatibilities. A necessary condition for this speciation scenario to work is that the two Wolbachiastrains can stably coexist. The following study investigates this problem analysing a mathematical model with two host populations and migration between them. We show that the stability of bidirectional CI can be fully described in terms of a critical migration rate which is defined as the highest migration below which a stable coexistence of two Wolbachiastrains is possible. For some special cases we could derive analytical solutions for the critical migration rate; for the general case estimations of the critical migration rate are given. Our main finding is that bidirectional CI can stably persist in the face of high migration and can be as high as over 15% per generation for CI levels observed in nature. These results have implications for the potential of Wolbachiato promote genetic divergence and speciation in their hosts.

The invasion and coexistence of competing Wolbachia strains

Cytoplasmic incompatibility between arthropods infected with different strains of Wolbachia has been proposed as an important mechanism for speciation. However, a basic requirement for this mechanism is the coexistence of different strains in neighbouring populations. Here we test whether this required coexistence is possible in a spatial context. Continuous-time models for the behaviour of one and two strains of Wolbachia within a single well-mixed population demonstrate the Allee effect and founder control, such that one strain is always driven extinct. In contrast, discretised spatial models show patchy persistence of the two strains although coexistence within the same habitat is rare. A simplified model of such founder control suggests that it is fragmentation of (or barriers within) the habitat rather than space itself that leads to persistence.