The hitchhiker's guide to Europe: the infection dynamics of an ongoing Wolbachia invasion and mitochondrial selective sweep in Rhagoletis cerasi

Spatially Varying Wolbachia Frequencies Reveal the Invasion Origin of an Agricultural Pest Recently Introduced From Europe to North America

The introduction of non-native species across the world represents a major global challenge. Retracing invasion origin is an important first step in understanding the invasion process, often requiring detailed sampling within the native range. Insect species frequently host Wolbachia, a widespread endosymbiotic bacterium that manipulates host reproduction to increase infected female fitness. Here, we draw on the spatial variation in infection frequencies of an actively spreading Wolbachia strain wCer2 to investigate the invasion origin of the European cherry fruit fly, Rhagoletis cerasi. This pest of cherries was introduced from Europe to North America within the last decade. First, we screen the introduced fly population for the presence of Wolbachia. The introduced populations lack the wCer2 strain and the strongly associated mitochondrial haplotype, suggesting strain absence due to founder effects with invading individuals originating from wCer2-uninfected native population(s). To narrow down geographic regions of invasion origin, we perform spatial interpolation of the wCer2 infection frequency across the native range and predict the infection frequency in unsampled regions. For this, we use an extensive dataset of R. cerasi infection covering 238 populations across Europe over 25 years, complemented with 14 additional populations analyzed for this study. We find that R. cerasi was unlikely introduced from wCer2-infected populations in Central and Western Europe. We propose wCer2-uninfected populations from Eastern Europe and the Mediterranean region as the most likely candidates for the invasion origin. This work utilizes Wolbachia as an indirect instrument to provide insights into the invasion source of R. cerasi in North America, revealing yet another application for this multifaceted heritable endosymbiont. Given the prevalence of biological invasions, rapidly uncovering invasion origins gives fundamental insights into how invasive species adapt to new environments.

Modelling suggests Wolbachia-induced cytoplasmic incompatibility in oak gall wasps with cyclical parthenogenesis

Oak gall wasps typically exhibit a life cycle with one sexual and one asexual generation each year. These wasps can carry various endosymbionts, one of which is the maternally inherited bacterium Wolbachia that can induce several reproductive manipulations on its host. Cytoplasmic incompatibility (CI) has been described as the most prominent of these manipulations. CI leads to embryonic mortality in the hosts' offspring when infected males mate with either uninfected females or with females that harbour different Wolbachia strains. It has been hypothesized that Wolbachia can induce CI in oak gall wasps. To address this hypothesis, we derived a mathematical model to investigate the spread of a bacterial infection in naive populations and to determine the plausibility of CI occurrence. To validate our model, we used published data from Wolbachia-infected Belonocnema kinseyi populations in two approaches. Our first approach uses measurements of infection frequencies and maternal transmission in the sexual generation. For the second approach, we extended the model to compare predictions to estimates of mtDNA-haplotypes, which, like Wolbachia, are maternally inherited, and can therefore be associated with the infection. Both approaches indicate that CI is present in these populations. Our model can be generalized to investigate the occurrence of CI not only for oak gall wasps but also for other species.

Genetic Diversity of Dengue Vector Aedes albopictus Collected from South Korea, Japan, and Laos

Aedes albopictus is native to Southeast Asia and has emerged as a major vector for vector-borne diseases that are spreading rapidly worldwide. Recent studies have shown that Ae. albopictus populations have different genetic groups dependent on their thermal adaptations; however, studies on Korean populations are limited. In this study, we analyzed the genetic diversity and structure of two mitochondrial genes (COI and ND5) and sixteen microsatellites in mosquitoes inhabiting Korea, Japan, and Laos. The results indicate that the Korean population has low genetic diversity, with an independent cluster distinct from the Laos population. Mixed clusters have also been observed in the Korean population. On the basis of these findings, two hypotheses are proposed. First, certain Korean populations are native. Second, some subpopulations that descended from the metapopulation (East Asian countries) were introduced to Japan before migrating to Korea. Furthermore, we previously demonstrated that Ae. albopictus appears to have been imported to Korea. In conclusion, the dengue-virus-carrying mosquitoes could migrate to Korea from Southeast Asian epidemic regions, where they can survive during the severe winter months. The key findings can be used to establish an integrated pest management strategy based on population genetics for the Korean Ae. albopictus population.

Microbiome of pear psyllids: A tale about closely related species sharing their endosymbionts

Psyllids are phloem-feeding insects that can transmit plant pathogens such as phytoplasmas, intracellular bacteria causing numerous plant diseases worldwide. Their microbiomes are essential for insect physiology and may also influence the capacity of vectors to transmit pathogens. Using 16S rRNA gene metabarcoding, we compared the microbiomes of three sympatric psyllid species associated with pear trees in Central Europe. All three species are able to transmit 'Candidatus Phytoplasma pyri', albeit with different efficiencies. Our results revealed potential relationships between insect biology and microbiome composition that varied during psyllid ontogeny and between generations in Cacopsylla pyri and C. pyricola, as well as between localities in C. pyri. In contrast, no variations related to psyllid life cycle and geography were detected in C. pyrisuga. In addition to the primary endosymbiont Carsonella ruddii, we detected another highly abundant endosymbiont (unclassified Enterobacteriaceae). C. pyri and C. pyricola shared the same taxon of Enterobacteriaceae which is related to endosymbionts harboured by other psyllid species from various families. In contrast, C. pyrisuga carried a different Enterobacteriaceae taxon related to the genus Sodalis. Our study provides new insights into host-symbiont interactions in psyllids and highlights the importance of host biology and geography in shaping microbiome structure.

Multiple invasions, Wolbachia and human-aided transport drive the genetic variability of Aedes albopictus in the Iberian Peninsula

The Asian tiger mosquito, Aedes albopictus, is one of the most invasive species in the world. Native to the tropical forests of Southeast Asia, over the past 30 years it has rapidly spread throughout tropical and temperate regions of the world. Its dramatic expansion has resulted in public health concerns as a consequence of its vector competence for at least 16 viruses. Previous studies showed that Ae. albopictus spread has been facilitated by human-mediated transportation, but much remains unknown about how this has affected its genetic attributes. Here we examined the factors that contributed to shaping the current genetic constitution of Ae. albopictus in the Iberian Peninsula, where the species was first found in 2004, by combining population genetics and Bayesian modelling. We found that both mitochondrial and nuclear DNA markers showed a lack of genetic structure and the presence of worldwide dominant haplotypes, suggesting regular introductions from abroad. Mitochondrial DNA showed little genetic diversity compared to nuclear DNA, likely explained by infection with maternally transmitted bacteria of the genus Wolbachia. Multilevel models revealed that greater mosquito fluxes (estimated from commuting patterns and tiger mosquito population distribution) and spatial proximity between sampling sites were associated with lower nuclear genetic distance, suggesting that rapid short- and medium-distance dispersal is facilitated by humans through vehicular traffic. This study highlights the significant role of human transportation in shaping the genetic attributes of Ae. albopictus and promoting regional gene flow, and underscores the need for a territorially integrated surveillance across scales of this disease-carrying mosquito.

A continental-scale survey of Wolbachia infections in blue butterflies reveals evidence of interspecific transfer and invasion dynamics

Infections by maternally inherited bacterial endosymbionts, especially Wolbachia, are common in insects and other invertebrates but infection dynamics across species ranges are largely under studied. Specifically, we lack a broad understanding of the origin of Wolbachia infections in novel hosts, and the historical and geographical dynamics of infections that are critical for identifying the factors governing their spread. We used Genotype-by-Sequencing data from previous population genomics studies for range-wide surveys of Wolbachia presence and genetic diversity in North American butterflies of the genus Lycaeides. As few as one sequence read identified by assembly to a Wolbachia reference genome provided high accuracy in detecting infections in host butterflies as determined by confirmatory PCR tests, and maximum accuracy was achieved with a threshold of only 5 sequence reads per host individual. Using this threshold, we detected Wolbachia in all but 2 of the 107 sampling localities spanning the continent, with infection frequencies within populations ranging from 0% to 100% of individuals, but with most localities having high infection frequencies (mean = 91% infection rate). Three major lineages of Wolbachia were identified as separate strains that appear to represent 3 separate invasions of Lycaeides butterflies by Wolbachia. Overall, we found extensive evidence for acquisition of Wolbachia through interspecific transfer between host lineages. Strain wLycC was confined to a single butterfly taxon, hybrid lineages derived from it, and closely adjacent populations in other taxa. While the other 2 strains were detected throughout the rest of the continent, strain wLycB almost always co-occurred with wLycA. Our demographic modeling suggests wLycB is a recent invasion. Within strain wLycA, the 2 most frequent haplotypes are confined almost exclusively to separate butterfly taxa with haplotype A1 observed largely in Lycaeides melissa and haplotype A2 observed most often in Lycaeides idas localities, consistent with either cladogenic mode of infection acquisition from a common ancestor or by hybridization and accompanying mutation. More than 1 major Wolbachia strain was observed in 15 localities. These results demonstrate the utility of using resequencing data from hosts to quantify Wolbachia genetic variation and infection frequency and provide evidence of multiple colonizations of novel hosts through hybridization between butterfly lineages and complex dynamics between Wolbachia strains.

Wolbachia in Black Spiny Whiteflies and Their New Parasitoid Wasp in Japan: Evidence of the Distinct Infection Status on Aleurocanthus camelliae Cryptic Species Complex

Wolbachia, an alphaproteobacterial reproductive parasite, can cause profound mitochondrial divergence in insects, which might eventually be a part of cryptic speciation. Aleurocanthus camelliae is a cryptic species complex consisting of several morphospecies and/or haplotypes that are genetically different but morphologically indistinctive. However, little is known about the Wolbachia infection status in these tea and Citrus pests. Thus, this study aimed to profile the diversity and phenotypic characteristics of Wolbachia natural infections in the A. camelliae cryptic species complex. A monophyletic strain of Wolbachia that infected the A. camelliae cryptic species complex (wAlec) with different patterns was discovered. Whiteflies that are morphologically identical to Aleurocanthus spiniferus (Aleurocanthus cf. A. spiniferus in Eurya japonica and A. spiniferus in Citrus) were grouped into uninfected populations, whereas the fixed infection was detected in A. camelliae B1 from Theaceae. The rapid evolution of wAlec was also found to occur through a high recombination event, which produced subgroups A and B in wAlec. It may also be associated with the non-cytoplasmic incompatibility (CI) phenotype of wAlec due to undetectable CI-related genes from phage WO (WOAlec). The current discovery of a novel cryptic species of A. camelliae led to a discussion about the oscillation hypothesis, which may provide insights on cryptic speciation, particularly on how specialization and host expansion have been recorded among these species. This study also identified a parasitoid wasp belonging to the genus Eretmocerus in A. camelliae, for the first time in Japan.

Wolbachia infection and genetic diversity of Italian populations of Philaenus spumarius, the main vector of Xylella fastidiosa in Europe

Philaenus spumarius is a cosmopolitan species that has become a major threat to European agriculture being recognized as the main vector of the introduced plant pathogen Xylella fastidiosa, the agent of the “olive quick decline syndrome”, a disease which is devastating olive orchards in southern Italy. Wolbachia are bacterial symbionts of many insects, frequently as reproductive parasites, sometime by establishing mutualistic relationships, able to spread within host populations. Philaenus spumarius harbors Wolbachia, but the role played by this symbiont is unknown and data on the infection prevalence within host populations are limited. Here, the Wolbachia infection rate was analyzed in relation to the geographic distribution and the genetic diversity of the Italian populations of P. spumarius. Analysis of the COI gene sequences revealed a geographically structured distribution of the three main mitochondrial lineages of P. spumarius. Wolbachia was detected in half of the populations sampled in northern Italy where most individuals belonged to the western-Mediterranean lineage. All populations sampled in southern and central Italy, where the individuals of the eastern-Mediterranean lineage were largely prevalent, were uninfected. Individuals of the north-eastern lineage were found only in populations from the Alps in the northernmost part of Italy, at high altitudes. In this area, Wolbachia infection reached the highest prevalence, with no difference between north-eastern and western-Mediterranean lineage. Analysis of molecular diversity of COI sequences suggested no significant effect of Wolbachia on population genetics of P. spumarius. Using the MLST approach, six new Wolbachia sequence types were identified. Using FISH, Wolbachia were observed within the host’s reproductive tissues and salivary glands. Results obtained led us to discuss the role of Wolbachia in P. spumarius, the factors influencing the geographic distribution of the infection, and the exploitation of Wolbachia for the control of the vector insect to reduce the spread of X. fastidiosa.

Wolbachia Infection in Native Populations of the Invasive Tawny Crazy Ant Nylanderia fulva

Antagonistic interactions can affect population growth and dispersal of an invasive species. Wolbachia are intracellular endosymbiont bacteria that infect arthropod and nematode hosts and are able to manipulate reproduction, which in some cases leads to cocladogenesis. Moreover, the presence of the strictly maternally transferred Wolbachia in a population can indirectly induce selective sweeps on the hosts' mitochondria. Ants have a Wolbachia infection rate of about 34%, which makes phylogenetic studies using mitochondrial markers vulnerable of being confounded by the effect of the endosymbiont. Nylanderia fulva is an invasive ant native to South America, considered a pest in the United States. Its distribution and biology are poorly known in its native range, and the taxonomic identity of this and its closely related species, Nylanderia pubens, has only recently been understood with the aid of molecular phylogenies. Aiming at estimating robust phylogenetic relationships of N. fulva in its native range, we investigated the presence and pattern of Wolbachia infection in populations of N. fulva from Argentina, part of its native range, to account for its possible effect on the host population structure. Using the ftsZ gene, 30 nests of N. fulva and four from sympatric Nylanderia species were screened for the presence of Wolbachia. We sequenced the MLST genes, the highly variable gene wsp, as well as glyQ, a novel target gene for which new primers were designed. Phylogeny of the ants was estimated using mtDNA (COI). We found supergroup A Wolbachia strains infecting 73% of N. fulva nests and two nests of Nylanderia sp. 1. Wolbachia phylogenetic tree inferred with MLST genes is partially congruent with the host phylogeny topology, with the exception of a lineage of strains shared by ants from different N. fulva clades. Furthermore, by comparing with Wolbachia sequences infecting other ants, we found that the strains infecting different N. fulva clades are not monophyletic. Our findings suggest there are three recent independent horizontally transmitted Wolbachia infections in N. fulva, and we found no evidence of influence of Wolbachia in the host mtDNA based phylogeny.

Substantial rearrangements, single nucleotide frameshift deletion and low diversity in mitogenome of Wolbachia-infected strepsipteran endoparasitoid in comparison to its tephritid hosts

Insect mitogenome organisation is highly conserved, yet, some insects, especially with parasitic life cycles, have rearranged mitogenomes. Furthermore, intraspecific mitochondrial diversity can be reduced by fitness-affecting bacterial endosymbionts like Wolbachia due to their maternal coinheritance with mitochondria. We have sequenced mitogenomes of the Wolbachia-infected endoparasitoid Dipterophagus daci (Strepsiptera: Halictophagidae) and four of its 22 known tephritid fruit fly host species using total genomic extracts of parasitised flies collected across > 700 km in Australia. This halictophagid mitogenome revealed extensive rearrangements relative to the four fly mitogenomes which exhibited the ancestral insect mitogenome pattern. Compared to the only four available other strepsipteran mitogenomes, the D. daci mitogenome had additional transpositions of one rRNA and two tRNA genes, and a single nucleotide frameshift deletion in nad5 requiring translational frameshifting or, alternatively, resulting in a large protein truncation. Dipterophagus daci displays an almost completely endoparasitic life cycle when compared to Strepsiptera that have maintained the ancestral state of free-living adults. Our results support the hypothesis that the transition to extreme endoparasitism evolved together with increased levels of mitogenome changes. Furthermore, intraspecific mitogenome diversity was substantially smaller in D. daci than the parasitised flies suggesting Wolbachia reduced mitochondrial diversity because of a role in D. daci fitness.

Unparalleled mitochondrial heteroplasmy and Wolbachia co-infection in the non-model bee, Amphylaeus morosus

Mitochondrial heteroplasmy is the occurrence of more than one type of mitochondrial DNA within a single individual. Although generally reported to occur in a small subset of individuals within a species, there are some instances of widespread heteroplasmy across entire populations. Amphylaeus morosus is an Australian native bee species in the diverse and cosmopolitan bee family Colletidae. This species has an extensive geographical range along the eastern Australian coast, from southern Queensland to western Victoria, covering approximately 2,000 km. Seventy individuals were collected from five localities across this geographical range and sequenced using Sanger sequencing for the mitochondrial cytochrome c oxidase subunit I (COI) gene. These data indicate that every individual had the same consistent heteroplasmic sites but no other nucleotide variation, suggesting two conserved and widespread heteroplasmic mitogenomes. Ion Torrent shotgun sequencing revealed that heteroplasmy occurred across multiple mitochondrial protein-coding genes and is unlikely explained by transposition of mitochondrial genes into the nuclear genome (NUMTs). DNA sequence data also demonstrated a consistent co-infection of Wolbachia across the A. morosus distribution with every individual infected with both bacterial strains. Our data are consistent with the presence of two mitogenomes within all individuals examined in this species and suggest a major divergence from standard patterns of mitochondrial inheritance. Because the host's mitogenome and the Wolbachia genome are genetically linked through maternal inheritance, we propose three possible hypotheses that could explain maintenance of the widespread and conserved co-occurring bacterial and mitochondrial genomes in this species.

Male Age and Wolbachia Dynamics: Investigating How Fast and Why Bacterial Densities and Cytoplasmic Incompatibility Strengths Vary

Endosymbionts can influence host reproduction and fitness to favor their maternal transmission. For example, endosymbiotic Wolbachia bacteria often cause cytoplasmic incompatibility (CI) that kills uninfected embryos fertilized by Wolbachia-modified sperm. Infected females can rescue CI, providing them a relative fitness advantage. Wolbachia-induced CI strength varies widely and tends to decrease as host males age. Since strong CI drives Wolbachia to high equilibrium frequencies, understanding how fast and why CI strength declines with male age is crucial to explaining age-dependent CI’s influence on Wolbachia prevalence. Here, we investigate if Wolbachia densities and/or CI gene (cif) expression covary with CI-strength variation and explore covariates of age-dependent Wolbachia-density variation in two classic CI systems. wRi CI strength decreases slowly with Drosophila simulans male age (6%/day), but wMel CI strength decreases very rapidly (19%/day), yielding statistically insignificant CI after only 3 days of Drosophila melanogaster adult emergence. Wolbachia densities and cif expression in testes decrease as wRi-infected males age, but both surprisingly increase as wMel-infected males age, and CI strength declines. We then tested if phage lysis, Octomom copy number (which impacts wMel density), or host immune expression covary with age-dependent wMel densities. Only host immune expression correlated with density. Together, our results identify how fast CI strength declines with male age in two model systems and reveal unique relationships between male age, Wolbachia densities, cif expression, and host immunity. We discuss new hypotheses about the basis of age-dependent CI strength and its contributions to Wolbachia prevalence.

Frequent Drivers, Occasional Passengers: Signals of Symbiont-Driven Seasonal Adaptation and Hitchhiking in the Pea Aphid, Acyrthosiphon pisum

Insects harbor a variety of maternally inherited bacterial symbionts. As such, variation in symbiont presence/absence, in the combinations of harbored symbionts, and in the genotypes of harbored symbiont species provide heritable genetic variation of potential use in the insects' adaptive repertoires. Understanding the natural importance of symbionts is challenging but studying their dynamics over time can help to elucidate the potential for such symbiont-driven insect adaptation. Toward this end, we studied the seasonal dynamics of six maternally transferred bacterial symbiont species in the multivoltine pea aphid (Acyrthosiphon pisum). Our sampling focused on six alfalfa fields in southeastern Pennsylvania, and spanned 14 timepoints within the 2012 growing season, in addition to two overwintering periods. To test and generate hypotheses on the natural relevance of these non-essential symbionts, we examined whether symbiont dynamics correlated with any of ten measured environmental variables from the 2012 growing season, including some of known importance in the lab. We found that five symbionts changed prevalence across one or both overwintering periods, and that the same five species underwent such frequency shifts across the 2012 growing season. Intriguingly, the frequencies of these dynamic symbionts showed robust correlations with a subset of our measured environmental variables. Several of these trends supported the natural relevance of lab-discovered symbiont roles, including anti-pathogen defense. For a seventh symbiont-Hamiltonella defensa-studied previously across the same study periods, we tested whether a reported correlation between prevalence and temperature stemmed not from thermally varying host-level fitness effects, but from selection on co-infecting symbionts or on aphid-encoded alleles associated with this bacterium. In general, such "hitchhiking" effects were not evident during times with strongly correlated Hamiltonella and temperature shifts. However, we did identify at least one time period in which Hamiltonella spread was likely driven by selection on a co-infecting symbiont-Rickettsiella viridis. Recognizing the broader potential for such hitchhiking, we explored selection on co-infecting symbionts as a possible driver behind the dynamics of the remaining six species. Out of twelve examined instances of symbiont dynamics unfolding across 2-week periods or overwintering spans, we found eight in which the focal symbiont underwent parallel frequency shifts under single infection and one or more co-infection contexts. This supported the idea that phenotypic variation created by the presence/absence of individual symbionts is a direct target for selection, and that symbiont effects can be robust under co-habitation with other symbionts. Contrastingly, in two cases, we found that selection may target phenotypes emerging from symbiont co-infections, with specific species combinations driving overall trends for the focal dynamic symbionts, without correlated change under single infection. Finally, in three cases-including the one described above for Hamiltonella-our data suggested that incidental co-infection with a (dis)favored symbiont could lead to large frequency shifts for "passenger" symbionts, conferring no apparent cost or benefit. Such hitchhiking has rarely been studied in heritable symbiont systems. We propose that it is more common than appreciated, given the widespread nature of maternally inherited bacteria, and the frequency of multi-species symbiotic communities across insects.

Testing the potential contribution of Wolbachia to speciation when cytoplasmic incompatibility becomes associated with host-related reproductive isolation

Endosymbiont‐induced cytoplasmic incompatibility (CI) may play an important role in arthropod speciation. However, whether CI consistently becomes associated or coupled with other host‐related forms of reproductive isolation (RI) to impede the transfer of endosymbionts between hybridizing populations and further the divergence process remains an open question. Here, we show that varying degrees of pre‐ and postmating RI exist among allopatric populations of two interbreeding cherry‐infesting tephritid fruit flies (Rhagoletis cingulata and R. indifferens) across North America. These flies display allochronic and sexual isolation among populations, as well as unidirectional reductions in egg hatch in hybrid crosses involving southwestern USA males. All populations are infected by a Wolbachia strain, wCin2, whereas a second strain, wCin3, only co‐infects flies from the southwest USA and Mexico. Strain wCin3 is associated with a unique mitochondrial DNA haplotype and unidirectional postmating RI, implicating the strain as the cause of CI. When coupled with nonendosymbiont RI barriers, we estimate the strength of CI associated with wCin3 would not prevent the strain from introgressing from infected southwestern to uninfected populations elsewhere in the USA if populations were to come into secondary contact and hybridize. In contrast, cytoplasmic–nuclear coupling may impede the transfer of wCin3 if Mexican and USA populations were to come into contact. We discuss our results in the context of the general paucity of examples demonstrating stable Wolbachia hybrid zones and whether the spread of Wolbachia among taxa can be constrained in natural hybrid zones long enough for the endosymbiont to participate in speciation.

Genome analyses of four Wolbachia strains and associated mitochondria of Rhagoletis cerasi expose cumulative modularity of cytoplasmic incompatibility factors and cytoplasmic hitchhiking across host populations

Background

The endosymbiont Wolbachia can manipulate arthropod reproduction and invade host populations by inducing cytoplasmic incompatibility (CI). Some host species are coinfected with multiple Wolbachia strains which may have sequentially invaded host populations by expressing different types of modular CI factor (cif) genes. The tephritid fruit fly Rhagoletis cerasi is a model for CI and Wolbachia population dynamics. It is associated with at least four Wolbachia strains in various combinations, with demonstrated (wCer2, wCer4), predicted (wCer1) or unknown (wCer5) CI phenotypes.

Results

We sequenced and assembled the draft genomes of the Wolbachia strains wCer1, wCer4 and wCer5, and compared these with the previously sequenced genome of wCer2 which currently invades R. cerasi populations. We found complete cif gene pairs in all strains: four pairs in wCer2 (three Type I; one Type V), two pairs in wCer1 (both Type I) and wCer4 (one Type I; one Type V), and one pair in wCer5 (Type IV). Wolbachia genome variant analyses across geographically and genetically distant host populations revealed the largest diversity of single nucleotide polymorphisms (SNPs) in wCer5, followed by wCer1 and then wCer2, indicative of their different lengths of host associations. Furthermore, mitogenome analyses of the Wolbachia genome-sequenced individuals in combination with SNP data from six European countries revealed polymorphic mitogenome sites that displayed reduced diversity in individuals infected with wCer2 compared to those without.

Conclusions

Coinfections with Wolbachia are common in arthropods and affect options for Wolbachia-based management strategies of pest and vector species already infected by Wolbachia. Our analyses of Wolbachia genomes of a host naturally coinfected by several strains unravelled signatures of the evolutionary dynamics in both Wolbachia and host mitochondrial genomes as a consequence of repeated invasions. Invasion of already infected populations by new Wolbachia strains requires new sets of functionally different cif genes and thereby may select for a cumulative modularity of cif gene diversity in invading strains. Furthermore, we demonstrated at the mitogenomic scale that repeated CI-driven Wolbachia invasions of hosts result in reduced mitochondrial diversity and hitchhiking effects. Already resident Wolbachia strains may experience similar cytoplasmic hitchhiking effects caused by the invading Wolbachia strain.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07906-6.

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.

Parallel Sequencing of Wolbachia wCer2 from Donor and Novel Hosts Reveals Multiple Incompatibility Factors and Genome Stability after Host Transfers

The application of Wolbachia in insect pest and vector control requires the establishment of genotypically stable host associations. The cytoplasmic incompatibility (CI) inducing Wolbachia strain wCer2 naturally occurs in the cherry fruit fly Rhagoletis cerasi as co-infection with other strains and was transferred to other fruit fly species by embryonic microinjections. We obtained wCer2 genome data from its native and three novel hosts, Drosophila simulans, Drosophila melanogaster, and Ceratitis capitata and assessed its genome stability, characteristics, and CI factor (cif) genes. De novo assembly was successful from Wolbachia cell-enriched singly infected D. simulans embryos, with minimal host and other bacterial genome traces. The low yield of Wolbachia sequence reads from total genomic extracts of one multiply infected R. cerasi pupa and one singly infected C. capitata adult limited de novo assemblies but was sufficient for comparative analyses. Across hosts wCer2 was stable in genome synteny and content. Polymorphic nucleotide sites were found in wCer2 of each host; however, only one nucleotide was different between R. cerasi and C. capitata, and none between replicated D. simulans lines. The wCer2 genome is highly similar to wAu (D. simulans), wMel (D. melanogaster), and wRec (Drosophila recens). In contrast to wMel and wRec (each with one cif gene pair) and wAu (without any cif genes), wCer2 has three pairs of Type I cif genes, and one Type V cifB gene without a cifA complement. This may explain previously reported CI patterns of wCer2, including incomplete rescue of its own CI modification in three novel host species.

Long-term spatiotemporal genetic structure of an accidental parasitoid introduction, and local changes in prevalence of its associated Wolbachia symbiont

Population bottlenecks associated with founder events strongly impact the establishment and genetic makeup of populations. In addition to their genotype, founding individuals also bring along parasites, as well as symbionts that can manipulate the phenotype of their host, affecting the host population establishment, dynamics and evolution. Thus, to understand introduction, invasion, and spread, we should identify the roles played by accompanying symbionts. In 1991, the parasitoid wasp, Hyposoter horticola, and its associated hyperparasitoid were accidentally introduced from the main Åland islands, Finland, to an isolated island in the archipelago, along with their host, the Glanville fritillary butterfly. Though the receiving island was unoccupied, the butterfly was present on some of the small islands in the vicinity. The three introduced species have persisted locally ever since. A strain of the endosymbiotic bacterium Wolbachia has an intermediate prevalence in the parasitoid H. horticola across the main Åland population. The infection increases its susceptibility of to hyperparasitism. We investigated the establishment and spread of the parasitoid, along with patterns of prevalence of its symbiont using 323 specimens collected between 1992 and 2013, from five localities across Åland, including the source and introduced populations. Using 14 microsatellites and one mitochondrial marker, we suggest that the relatively diverse founding population and occasional migration between islands might have facilitated the persistence of all isolated populations, despite multiple local population crashes. We also show that where the hyperparasitoid is absent, and thus selection against infected wasp genotypes is relaxed, there is near-fixation of Wolbachia.

Comparative genome sequencing reveals insights into the dynamics of Wolbachia in native and invasive cherry fruit flies

Wolbachia is a maternally inherited obligate endosymbiont that can induce a wide spectrum of effects in its host, ranging from mutualism to reproductive parasitism. At the genomic level, recombination within and between strains, transposable elements, and horizontal transfer of strains between host species make Wolbachia an evolutionarily dynamic bacterial system. The invasive cherry fruit fly Rhagoletis cingulata arrived in Europe from North America ~40 years ago, where it now co‐occurs with the native cherry pest R. cerasi. This shared distribution has been proposed to have led to the horizontal transfer of different Wolbachia strains between the two species. To better understand transmission dynamics, we performed a comparative genome study of the strain wCin2 in its native United States and invasive European populations of R. cingulata with wCer2 in European R. cerasi. Previous multilocus sequence genotyping (MLST) of six genes implied that the source of wCer2 in R. cerasi was wCin2 from R. cingulata. However, we report genomic evidence discounting the recent horizontal transfer hypothesis for the origin of wCer2. Despite near identical sequences for the MLST markers, substantial sequence differences for other loci were found between wCer2 and wCin2, as well as structural rearrangements, and differences in prophage, repetitive element, gene content, and cytoplasmic incompatibility inducing genes. Our study highlights the need for whole‐genome sequencing rather than relying on MLST markers for resolving Wolbachia strains and assessing their evolutionary dynamics.

Host plants and Wolbachia shape the population genetics of sympatric herbivore populations

Changing climate and land-use practices have the potential to bring previously isolated populations of pest insects into new sympatry. This heightens the need to better understand how differing patterns of host-plant association, and unique endosymbionts, serve to promote genetic isolation or integration. We addressed these factors in populations of potato psyllid, Bactericera cockerelli (Šulc), a generalist herbivore that vectors a bacterial pathogen (Candidatus Liberibacter solanacearum, causal pathogen of zebra chip disease) of potato (Solanum tuberosum L.). Genome-wide SNP data revealed two major genetic clusters-psyllids collected from potato crops were genetically similar to psyllids found on a common weed, Lycium spp., but dissimilar from those found on another common non-crop host, Solanum dulcamara L. Most psyllids found on Lycium spp. and potato represented a single mitochondrial cytochrome oxidase I (COI) haplotype that has been suggested to not be native to the region, and whose arrival may have been concurrent with zebra chip disease first emerging. The putatively introduced COI haplotype usually co-occurred with endosymbiotic Wolbachia, while the putatively resident COI haplotype generally did not. Genetic intermediates between the two genetic populations of insects were rare, consistent with recent sympatry or reproductive isolation, although admixture patterns of apparent hybrids were consistent with introgression of genes from introduced into resident populations. Our results suggest that both host-plant associations and endosymbionts are shaping the population genetic structure of sympatric psyllid populations associated with different non-crop hosts. It is of future interest to explicitly examine vectorial capacity of the two populations and their potential hybrids, as population structure and hybridization might alter regional vector capacity and disease outbreaks.

Evidence for spatial clines and mixed geographic modes of speciation for North American cherry-infesting Rhagoletis (Diptera: Tephritidae) flies

An important criterion for understanding speciation is the geographic context of population divergence. Three major modes of allopatric, parapatric, and sympatric speciation define the extent of spatial overlap and gene flow between diverging populations. However, mixed modes of speciation are also possible, whereby populations experience periods of allopatry, parapatry, and/or sympatry at different times as they diverge. Here, we report clinal patterns of variation for 21 nuclear-encoded microsatellites and a wing spot phenotype for cherry-infesting Rhagoletis (Diptera: Tephritidae) across North America consistent with these flies having initially diverged in parapatry followed by a period of allopatric differentiation in the early Holocene. However, mitochondrial DNA (mtDNA) displays a different pattern; cherry flies at the ends of the clines in the eastern USA and Pacific Northwest share identical haplotypes, while centrally located populations in the southwestern USA and Mexico possess a different haplotype. We hypothesize that the mitochondrial difference could be due to lineage sorting but more likely reflects a selective sweep of a favorable mtDNA variant or the spread of an endosymbiont. The estimated divergence time for mtDNA suggests possible past allopatry, secondary contact, and subsequent isolation between USA and Mexican fly populations initiated before the Wisconsin glaciation. Thus, the current genetics of cherry flies may involve different mixed modes of divergence occurring in different portions of the fly's range. We discuss the need for additional DNA sequencing and quantification of prezygotic and postzygotic reproductive isolation to verify the multiple mixed-mode hypothesis for cherry flies and draw parallels from other systems to assess the generality that speciation may commonly involve complex biogeographies of varying combinations of allopatric, parapatric, and sympatric divergence.

Large scale genome reconstructions illuminate Wolbachia evolution

Wolbachia is an iconic example of a successful intracellular bacterium. Despite its importance as a manipulator of invertebrate biology, its evolutionary dynamics have been poorly studied from a genomic viewpoint. To expand the number of Wolbachia genomes, we screen over 30,000 publicly available shotgun DNA sequencing samples from 500 hosts. By assembling over 1000 Wolbachia genomes, we provide a substantial increase in host representation. Our phylogenies based on both core-genome and gene content provide a robust reference for future studies, support new strains in model organisms, and reveal recent horizontal transfers amongst distantly related hosts. We find various instances of gene function gains and losses in different super-groups and in cytoplasmic incompatibility inducing strains. Our Wolbachia-host co-phylogenies indicate that horizontal transmission is widespread at the host intraspecific level and that there is no support for a general Wolbachia-mitochondrial synchronous divergence.

By greatly expanding the number of assembled genomes for Wolbachia (a group of intracellular bacteria) and constructing robust phylogenies, this study finds strong rate heterogeneity among Wolbachiapopulations and no support for synchronous divergence between Wolbachia and host mitochondria.

Ongoing Coevolution of Wolbachia and a Widespread Invasive Ant, Anoplolepis gracilipes

While Wolbachia are commonly found among arthropods, intraspecific infection rates can vary substantially across the geographic populations. Here we report nearly 100% prevalence of Wolbachia in the global populations of the yellow crazy ant, Anoplolepis gracilipes. To understand coevolutionary history between Wolbachia and A. gracilipes, we identified single nucleotide polymorphisms (SNPs) in Wolbachia from the ant across 12 geographical regions and compared the phylogeny of SNP-based Wolbachia to patterns of the ant's mitochondrial DNA (mtDNA) variation. Our results revealed a strong concordance between phylogenies of Wolbachia and host mtDNA, providing immediate evidence of co-divergence. Among eight identified SNP loci separating the genetic clusters of Wolbachia, seven loci are located in potential protein-coding genes, three of which being non-synonymous SNPs that may influence gene functions. We found a Wolbachia hypothetical protein gene with signature of positive selection. These findings jointly allow us to characterize Wolbachia-ant coevolution and also raise a question about mechanism(s) underlying maintenance of high prevalence of Wolbachia during the colonization of this invasive ant.

Wolbachia-Mitochondrial DNA Associations in Transitional Populations of Rhagoletis cerasi

The endosymbiont Wolbachia can manipulate arthropod host reproduction by inducing cytoplasmic incompatibility (CI), which results in embryonic mortality when infected males mate with uninfected females. A CI-driven invasion of Wolbachia can result in a selective sweep of associated mitochondrial haplotype. The co-inheritance of Wolbachia and host mitochondrial DNA can therefore provide significant information on the dynamics of an ongoing Wolbachia invasion. Therefore, transition zones (i.e., regions where a Wolbachia strain is currently spreading from infected to uninfected populations) represent an ideal area to investigate the relationship between Wolbachia and host mitochondrial haplotype. Here, we studied Wolbachia-mitochondrial haplotype associations in the European cherry fruit fly, Rhagoletis cerasi, in two transition zones in the Czech Republic and Hungary, where the CI-inducing strain wCer2 is currently spreading. The wCer2-infection status of 881 individuals was compared with the two known R. cerasi mitochondrial haplotypes, HT1 and HT2. In accordance with previous studies, wCer2-uninfected individuals were associated with HT1, and wCer2-infected individuals were mainly associated with HT2. We found misassociations only within the transition zones, where HT2 flies were wCer2-uninfected, suggesting the occurrence of imperfect maternal transmission. We did not find any HT1 flies that were wCer2-infected, suggesting that Wolbachia was not acquired horizontally. Our study provides new insights into the dynamics of the early phase of a Wolbachia invasion.

Wolbachia have made it twice: Hybrid introgression between two sister species of Eurema butterflies

Wolbachia, cytoplasmically inherited endosymbionts of arthropods, are known to hijack their host reproduction in various ways to increase their own vertical transmission. This may lead to the selective sweep of associated mitochondria, which can have a large impact on the evolution of mitochondrial lineages. In Japan, two different Wolbacahia strains (wCI and wFem) are found in two sister species of pierid butterflies, Eurema mandarina and Eurema hecabe. In both species, females infected with wCI (C females) produce offspring with a nearly 1:1 sex ratio, while females infected with both wCI and wFem (CF females) produce all-female offspring. Previous studies have suggested the historical occurrence of hybrid introgression in C individuals between the two species. Furthermore, hybrid introgression in CF individuals is suggested by the distinct mitochondrial lineages between C females and CF females of E. mandarina. In this study, we performed phylogenetic analyses based on nuclear DNA and mitochondrial DNA markers of E. hecabe with previously published data on E. mandarina. We found that the nuclear DNA of this species significantly diverged from that of E. mandarina. By contrast, mitochondrial DNA haplotypes comprised two clades, mostly reflecting Wolbachia infection status rather than the individual species. Collectively, our results support the previously suggested occurrence of two independent historical events wherein the cytoplasms of CF females and C females moved between E. hecabe and E. mandarina through hybrid introgression.

Wolbachia pipientis Associated With Tephritid Fruit Fly Pests: From Basic Research to Applications

Members of the true fruit flies (family Tephritidae) are among the most serious agricultural pests worldwide, whose control and management demands large and costly international efforts. The need for cost-effective and environmentally friendly integrated pest management (IPM) has led to the development and implementation of autocidal control strategies. These approaches include the widely used sterile insect technique and the incompatible insect technique (IIT). IIT relies on maternally transmitted bacteria (namely Wolbachia) to cause a conditional sterility in crosses between released mass-reared Wolbachia-infected males and wild females, which are either uninfected or infected with a different Wolbachia strain (i.e., cytoplasmic incompatibility; CI). Herein, we review the current state of knowledge on Wolbachia-tephritid interactions including infection prevalence in wild populations, phenotypic consequences, and their impact on life history traits. Numerous pest tephritid species are reported to harbor Wolbachia infections, with a subset exhibiting high prevalence. The phenotypic effects of Wolbachia have been assessed in very few tephritid species, due in part to the difficulty of manipulating Wolbachia infection (removal or transinfection). Based on recent methodological advances (high-throughput DNA sequencing) and breakthroughs concerning the mechanistic basis of CI, we suggest research avenues that could accelerate generation of necessary knowledge for the potential use of Wolbachia-based IIT in area-wide integrated pest management (AW-IPM) strategies for the population control of tephritid pests.

Multiplex PCR for identification of two butterfly sister species: Eurema mandarina and Eurema hecabe

Objective

In insects, closely related species are often difficult or impossible to distinguish solely by morphological traits. Mitochondrial DNA (mtDNA) markers are often useful and reliable for distinguishing closely related species. However, useful mtDNA markers can be unavailable, particularly when such species pairs experienced hybrid introgression in the past. Although polymorphic nuclear DNA markers would be necessary to distinguish such species pairs, recombination, multiple copies, and slower mutation rates of the nuclear DNA compared with those of mtDNA often make it challenging. The objective of this study was to develop a multiplex polymerase chain reaction that can reliably amplify and distinguish the Tpi sequences of Eurema mandarina and Eurema hecabe.

Results

We successfully analyzed the nucleotide sequences of the Z chromosome-linked triose phosphate isomerase (Tpi) gene to develop a multiplex polymerase chain reaction (PCR) that amplified ca. 120-bp products for E. mandarina and ca. 375-bp products for E. hecabe. We suggest that multiplex PCR using Tpi with appropriately designed primers can be used to accurately and reliably distinguish between other closely related Lepidoptera species.

Wolbachia-infected ant colonies have increased reproductive investment and an accelerated life cycle

Wolbachia is a widespread genus of maternally transmitted endosymbiotic bacteria that often manipulates the reproductive strategy and life history of its hosts to favor its own transmission. Wolbachia-mediated phenotypic effects are well characterized in solitary hosts, but effects in social hosts are unclear. The invasive pharaoh ant, Monomorium pharaonis, shows natural variation in Wolbachia infection between colonies and can be readily bred under laboratory conditions. We previously showed that Wolbachia-infected pharaoh ant colonies had more queen-biased sex ratios than uninfected colonies, which is expected to favor the spread of maternally transmitted Wolbachia Here, we further characterize the effects of Wolbachia on the short- and longer-term reproductive and life history traits of pharaoh ant colonies. First, we characterized the reproductive differences between naturally infected and uninfected colonies at three discrete time points and found that infected colonies had higher reproductive investment (i.e. infected colonies produced more new queens), particularly when existing colony queens were 3 months old. Next, we compared the long-term growth and reproduction dynamics of infected and uninfected colonies across their whole life cycle. Infected colonies had increased colony-level growth and early colony reproduction, resulting in a shorter colony life cycle, when compared with uninfected colonies.

New Data on the Range Expansion of the Thaumetopoea pityocampa (Lepidoptera: Notodontidae) 'ENA clade' in Greece: The Role of Bacterial Endosymbionts

The pine processionary moth, Thaumetopoea pityocampa (Denis and Schiffermüller), is an important insect in the Mediterranean region, as it defoliates pines and its urticating hairs can cause allergic reactions in humans and animals. Moreover, this species exhibits an interesting genetic structure as recently a distinct East-North African mtDNA lineage ('ENA clade') has been described. This clade has been recently detected in Greek populations where it has currently expanded its range by replacing the 'endemic' T. pityocampa lineages. Here, we report new data on the rapid spread of 'ENA clade' in the Greek island Evoia in only a few years. As the underlying mechanisms of the 'ENA clade' range expansion has not been studied so far, we screened T. pityocampa for an infection with the heritable bacterial endosymbionts Wolbachia (Bacteria: Anaplasmataceae), Cardinium (Bacteria: Bacteroidaceae), Rickettsia (Bacteria: Rickettsiaceae) and Spiroplasma (Bacteria: Spiroplasmataceae). These bacteria can manipulate the reproduction of infected hosts, something that could potentially explain the rapid spread of 'ENA clade' lineage. Therefore, we screened 28 individuals that exhibited T. pityocampa 'ENA clade' and 'endemic' T. pityocampa haplotypes from nine populations scattered all over Greece. None of them was infected with any of the four endosymbionts, suggesting that these bacteria do not cause reproductive manipulations in T. pityocampa lineages and, thus, other factors should be explored in future research efforts.

Host plant-related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi

Elucidating the mechanisms and conditions facilitating the formation of biodiversity are central topics in evolutionary biology. A growing number of studies imply that divergent ecological selection may often play a critical role in speciation by counteracting the homogenising effects of gene flow. Several examples involve phytophagous insects, where divergent selection pressures associated with host plant shifts may generate reproductive isolation, promoting speciation. Here, we use ddRADseq to assess the population structure and to test for host‐related genomic differentiation in the European cherry fruit fly, Rhagoletis cerasi (L., 1758) (Diptera: Tephritidae). This tephritid is distributed throughout Europe and western Asia, and has adapted to two different genera of host plants, Prunus spp. (cherries) and Lonicera spp. (honeysuckle). Our data imply that geographic distance and geomorphic barriers serve as the primary factors shaping genetic population structure across the species range. Locally, however, flies genetically cluster according to host plant, with consistent allele frequency differences displayed by a subset of loci between Prunus and Lonicera flies across four sites surveyed in Germany and Norway. These 17 loci display significantly higher F_ST values between host plants than others. They also showed high levels of linkage disequilibrium within and between Prunus and Lonicera flies, supporting host‐related selection and reduced gene flow. Our findings support the existence of sympatric host races in R. cerasi embedded within broader patterns of geographic variation in the fly, similar to the related apple maggot, Rhagoletis pomonella, in North America.

Genetic Diversity and Wolbachia Infection Patterns in a Globally Distributed Invasive Ant

Understanding the phylogeographic history of an invasive species may facilitate reconstructing the history and routes of its invasion. The longhorn crazy ant, Paratrechina longicornis, is a ubiquitous agricultural and household pest throughout much of the tropics and subtropics, but little is known about the history of its spread. Here, we examine worldwide genetic variation in P. longicornis and its associated Wolbachia bacterial symbionts. Analyses of mtDNA sequences of 248 P. longicornis workers (one per colony) from 13 geographic regions reveal two highly diverged mtDNA clades that co-occur in most of the geographic regions. These two mtDNA clades are associated with different Wolbachia infection patterns, but are not congruent with patterns of nDNA (microsatellite) variation. Multilocus sequence typing reveals two distinct Wolbachia strains in P. longicornis, namely, wLonA and wLonF. The evolutionary histories of these two strains differ; wLonA appears to be primarily transmitted maternally, and patterns of mtDNA and nDNA variation and wLonA infection status are consistent with a relatively recent Wolbachia-induced selective sweep. In contrast, the observed patterns of mtDNA variation and wLonF infections suggest frequent horizontal transfer and losses of wLonF infections. The lack of nDNA structure among sampled geographic regions coupled with the finding that numerous mtDNA haplotypes are shared among regions implies that inadvertent long-distance movement through human commerce is common in P. longicornis and has shaped the genetic structure of this invasive ant worldwide.

Geographic and Ecological Dimensions of Host Plant-Associated Genetic Differentiation and Speciation in the Rhagoletis cingulata (Diptera: Tephritidae) Sibling Species Group

Ascertaining the causes of adaptive radiation is central to understanding how new species arise and come to vary with their resources. The ecological theory posits adaptive radiation via divergent natural selection associated with novel resource use; an alternative suggests character displacement following speciation in allopatry and then secondary contact of reproductively isolated but ecologically similar species. Discriminating between hypotheses, therefore, requires the establishment of a key role for ecological diversification in initiating speciation versus a secondary role in facilitating co-existence. Here, we characterize patterns of genetic variation and postzygotic reproductive isolation for tephritid fruit flies in the Rhagoletis cingulata sibling species group to assess the significance of ecology, geography, and non-adaptive processes for their divergence. Our results support the ecological theory: no evidence for intrinsic postzygotic reproductive isolation was found between two populations of allopatric species, while nuclear-encoded microsatellites implied strong ecologically based reproductive isolation among sympatric species infesting different host plants. Analysis of mitochondrial DNA suggested, however, that cytoplasmic-related reproductive isolation may also exist between two geographically isolated populations within R cingulata. Thus, ecology associated with sympatric host shifts and cytoplasmic effects possibly associated with an endosymbiont may be the key initial drivers of the radiation of the R. cingulata group.

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.

Cuticular hydrocarbons as potential mediators of cryptic species divergence in a mutualistic ant association

Upon advances in sequencing techniques, more and more morphologically identical organisms are identified as cryptic species. Often, mutualistic interactions are proposed as drivers of diversification. Species of the neotropical parabiotic ant association between Crematogaster levior and Camponotus femoratus are known for highly diverse cuticular hydrocarbon (CHC) profiles, which in insects serve as desiccation barrier but also as communication cues. In the present study, we investigated the association of the ants' CHC profiles with genotypes and morphological traits, and discovered cryptic species pairs in both genera. To assess putative niche differentiation between the cryptic species, we conducted an environmental association study that included various climate variables, canopy cover, and mutualistic plant species. Although mostly sympatric, the two Camponotus species seem to prefer different climate niches. However in the two Crematogaster species, we could not detect any differences in niche preference. The strong differentiation in the CHC profiles may thus suggest a possible role during speciation itself either by inducing assortative mating or by reinforcing sexual selection after the speciation event. We did not detect any further niche differences in the environmental parameters tested. Thus, it remains open how the cryptic species avoid competitive exclusion, with scope for further investigations.

The association between mitochondrial genetic variation and reduced colony fitness in an invasive wasp

Despite the mitochondrion's long-recognized role in energy production, mitochondrial DNA (mtDNA) variation commonly found in natural populations was assumed to be effectively neutral. However, variation in mtDNA has now been increasingly linked to phenotypic variation in life history traits and fitness. We examined whether the relative fitness in native and invasive common wasp (Vespula vulgaris) populations in Belgium and New Zealand (NZ), respectively, can be linked to mtDNA variation. Social wasp colonies in NZ were smaller with comparatively fewer queen cells, indicating a reduced relative fitness in the invaded range. Interestingly, queen cells in this population were significantly larger leading to larger queen offspring. By sequencing 1,872 bp of the mitochondrial genome, we determined mitochondrial haplotypes and detected reduced genetic diversity in NZ. Three common haplotypes in NZ frequently produced many queens, whereas the four rare haplotypes produced significantly fewer or no queens. The entire mitochondrial genome for each of these haplotypes was sequenced to identify polymorphisms associated with fitness reduction. We found 16 variable sites; however, no nonsynonymous mutation that was clearly causing impaired mitochondrial function was detected. We discuss how detected variants may alter secondary structures, gene expression or mito-nuclear interactions, or could be associated with nuclear-encoded variation. Whatever the ultimate mechanism, we show reduced fitness and mtDNA variation in an invasive wasp population as well as specific mtDNA variants associated with fitness variation within this population. Ours is one of only a few studies that confirm fitness impacts of mtDNA variation in wild nonmodel populations.

Dynamics of an Ongoing Wolbachia Spread in the European Cherry Fruit Fly, Rhagoletis cerasi (Diptera: Tephritidae)

Numerous terrestrial arthropods are infected with the alphaproteobacterium Wolbachia. This endosymbiont is usually transmitted vertically from infected females to their offspring and can alter the reproduction of hosts through various manipulations, like cytoplasmic incompatibility (CI), enhancing its spread in new host populations. Studies on the spatial and temporal dynamics of Wolbachia under natural conditions are scarce. Here, we analyzed Wolbachia infection frequencies in populations of the European cherry fruit fly, Rhagoletis cerasi (L.), in central Germany—an area of an ongoing spread of the CI-inducing strain wCer2. In total, 295 individuals from 19 populations were PCR-screened for the presence of wCer2 and their mitochondrial haplotype. Results were compared with historic data to understand the infection dynamics of the ongoing wCer2 invasion. An overall wCer2 infection frequency of about 30% was found, ranging from 0% to 100% per population. In contrast to an expected smooth transition from wCer2-infected to completely wCer2-uninfected populations, a relatively scattered infection pattern across geography was observed. Moreover, a strong Wolbachia-haplotype association was detected, with only a few rare misassociations. Our results show a complex dynamic of an ongoing Wolbachia spread in natural field populations of R. cerasi.

Spatial spread of Wolbachia in Rhagoletis cerasi populations

The bacterial endosymbiont Wolbachia has been used to control insect pests owing to its ability to manipulate their life history and suppress infectious diseases. Therefore, knowledge on Wolbachia dynamics in natural populations is fundamental. The European cherry fruit fly, Rhagoletis cerasi, is infected with the Wolbachia strain wCer2, mainly present in southern and central European populations, and is currently spreading into wCer2-uninfected populations driven by high unidirectional cytoplasmic incompatibility. Here, we describe the distribution of wCer2 along two transition zones where the infection is spreading into wCer2-uninfected R. cerasi populations. Fine-scale sampling of 19 populations in the Czech Republic showed a smooth decrease of wCer2 frequency from south to north within a distance of less than 20 km. Sampling of 12 Hungarian populations, however, showed a sharp decline of wCer2 infection frequency within a few kilometres. We fitted a standard wave equation to our empirical data and estimated a Wolbachia wave speed of 1.9 km yr⁻¹ in the Czech Republic and 1.0 km yr⁻¹ in Hungary. Considering the univoltine life cycle and limited dispersal ability of R. cerasi, our study highlights a rapid Wolbachia spread in natural host populations.

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.

Experimental evidence that thermal selection shapes mitochondrial genome evolution

Mitochondria are essential organelles, found within eukaryotic cells, which contain their own DNA. Mitochondrial DNA (mtDNA) has traditionally been used in population genetic and biogeographic studies as a maternally-inherited and evolutionary-neutral genetic marker. However, it is now clear that polymorphisms within the mtDNA sequence are routinely non-neutral, and furthermore several studies have suggested that such mtDNA polymorphisms are also sensitive to thermal selection. These observations led to the formulation of the “mitochondrial climatic adaptation” hypothesis, for which all published evidence to date is correlational. Here, we use laboratory-based experimental evolution in the fruit fly, Drosophila melanogaster, to test whether thermal selection can shift population frequencies of two mtDNA haplogroups whose natural frequencies exhibit clinal associations with latitude along the Australian east-coast. We present experimental evidence that the thermal regime in which the laboratory populations were maintained drove changes in haplogroup frequencies across generations. Our results strengthen the emerging view that intra-specific mtDNA variants are sensitive to selection, and suggest spatial distributions of mtDNA variants in natural populations of metazoans might reflect adaptation to climatic environments rather than within-population coalescence and diffusion of selectively-neutral haplotypes across populations.

Co-speciation in bedbug Wolbachia parallel the pattern in nematode hosts

Wolbachia bacteria, vertically transmitted intracellular endosymbionts, are associated with two major host taxa in which they show strikingly different symbiotic modes. In some taxa of filarial nematodes, where Wolbachia are strictly obligately beneficial to the host, they show complete within- and among-species prevalence as well as co-phylogeny with their hosts. In arthropods, Wolbachia usually are parasitic; if beneficial effects occurs, they can be facultative or obligate, related to host reproduction. In arthropods, the prevalence of Wolbachia varies within and among taxa, and no co-speciation events are known. However, one arthropod species, the common bedbug Cimex lectularius was recently found to be dependent on the provision of biotin and riboflavin by Wolbachia, representing a unique case of Wolbachia providing nutritional and obligate benefits to an arthropod host, perhaps even in a mutualistic manner. Using the presence of presumably functional biotin gene copies, our study demonstrates that the obligate relationship is maintained at least in 10 out of 15 species of the genera Cimex and Paracimex. The remaining five species harboured Wolbachia as well, demonstrating the first known case of 100% prevalence of Wolbachia among higher arthropod taxa. Moreover, we show the predicted co-cladogenesis between Wolbachia and their bedbug hosts, also as the first described case of Wolbachia co-speciation in arthropods.

Reproductive Manipulators in the Bark Beetle Pityogenes chalcographus (Coleoptera: Curculionidae)-The Role of Cardinium, Rickettsia, Spiroplasma, and Wolbachia

Heritable bacterial endosymbionts can alter the biology of numerous arthropods. They can influence the reproductive outcome of infected hosts, thus affecting the ecology and evolution of various arthropod species. The spruce bark beetle Pityogenes chalcographus (L.) (Coleoptera: Curculionidae: Scolytinae) was reported to express partial, unidirectional crossing incompatibilities among certain European populations. Knowledge on the background of these findings is lacking; however, bacterial endosymbionts have been assumed to manipulate the reproduction of this beetle. Previous work reported low-density and low-frequency Wolbachia infections of P. chalcographus but found it unlikely that this infection results in reproductive alterations. The aim of this study was to test the hypothesis of an endosymbiont-driven incompatibility, other than Wolbachia, reflected by an infection pattern on a wide geographic scale. We performed a polymerase chain reaction (PCR) screening of 226 individuals from 18 European populations for the presence of the endosymbionts Cardinium, Rickettsia, and Spiroplasma, and additionally screened these individuals for Wolbachia. Positive PCR products were sequenced to characterize these bacteria. Our study shows a low prevalence of these four endosymbionts in P. chalcographus. We detected a yet undescribed Spiroplasma strain in a single individual from Greece. This is the first time that this endosymbiont has been found in a bark beetle. Further, Wolbachia was detected in three beetles from two Scandinavian populations and two new Wolbachia strains were described. None of the individuals analyzed were infected with Cardinium and Rickettsia. The low prevalence of bacteria found here does not support the hypothesis of an endosymbiont-driven reproductive incompatibility in P. chalcographus.

Diversity and distribution of Wolbachia in relation to geography, host plant affiliation and life cycle of a heterogonic gall wasp

BACKGROUND

The maternally inherited endosymbiont Wolbachia is widespread in arthropods and nematodes and can play an important role in the ecology and evolution of its host through reproductive manipulation. Here, we survey Wolbachia in Belonocnema treatae, a widely distributed North American cynipid gall forming wasp that exhibits regional host specialization on three species of oaks and alternation of sexually and asexually reproducing generations. We investigated whether patterns of Wolbachia infection and diversity in B. treatae are associated with the insect's geographic distribution, host plant association, life cycle, and mitochondrial evolutionary history.

RESULTS

Screening of 463 individuals from 23 populations including sexual and asexual generations from all three host plants across the southern U.S. showed an average infection rate of 56% with three common Wolbachia strains: wTre1-3 and an additional rare variant wTre4. Phylogenetic analysis based on wsp showed that these strains are unrelated and likely independently inherited. We found no difference in Wolbachia infection frequency among host plant associated populations or between the asexual and sexual generations, or between males and females of the sexual generation. Partially incomplete Wolbachia transmission rates might explain the occurrence of uninfected individuals. A parallel analysis of the mitochondrial cytochrome oxidase I gene in B. treatae showed high mtDNA haplotype diversity in both infected and uninfected populations suggesting an ancestral infection by Wolbachia as well as a clear split between eastern and western B. treatae mtDNA clades with a sequence divergence of > 6%. The strain wTre1 was present almost exclusively in the western clade while wTre2 and wTre3 occur almost exclusively in eastern populations. In contrast, the same strains co-occur as double-infections in Georgia and triple-infections in two populations in central Florida.

CONCLUSIONS

The diversity of Wolbachia across geographically and genetically distinct populations of B. treatae and the co-occurrence of the same strains within three populations highlights the complex infection dynamics in this system. Moreover, the association of distinct Wolbachia strains with mitochondrial haplotypes of its host in populations infected by different Wolbachia strains suggests a potential role of the endosymbiont in reproductive isolation in B. treatae.

DNA barcoding mosquitoes: advice for potential prospectors

Mosquitoes’ importance as vectors of pathogens that drive disease underscores the importance of precise and comparable methods of taxa identification among their species. While several molecular targets have been used to study mosquitoes since the initiation of PCR in the 1980s, its application to mosquito identification took off in the early 1990s. This review follows the research's recent journey into the use of mitochondrial DNA (mtDNA) cytochrome oxidase 1 (COI or COX1) as a DNA barcode target for mosquito species identification – a target whose utility for discriminating mosquitoes is now escalating. The pros and cons of using a mitochondrial genome target are discussed with a broad sweep of the mosquito literature suggesting that nuclear introgressions of mtDNA sequences appear to be uncommon and that the COI works well for distantly related taxa and shows encouraging utility in discriminating more closely related species such as cryptic/sibling species groups. However, the utility of COI in discriminating some closely related groups can be problematic and investigators are advised to proceed with caution as problems with incomplete lineage sorting and introgression events can result in indistinguishable COI sequences appearing in reproductively independent populations. In these – if not all – cases, it is advisable to run a nuclear marker alongside the mtDNA and thus the utility of the ribosomal DNA – and in particular the internal transcribed spacer 2 – is also briefly discussed as a useful counterpoint to the COI.

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].

Wolbachia-induced meiotic drive and feminization is associated with an independent occurrence of selective mitochondrial sweep in a butterfly

Maternally inherited Wolbachia endosymbionts manipulate arthropod reproduction in various ways. In the butterfly Eurema mandarina, a cytoplasmic incompatibility-inducing Wolbachia strain wCI and the associated mtDNA haplotypes are known to originate from the sister species Eurema hecabe, which offered a good case study for microbe-mediated hybrid introgression. Besides wCI, some females with the Z0 karyotype harbour a distinct Wolbachia strain wFem, which causes all-female production by meiotic drive and feminization. We report that a considerable proportion of E. mandarina females (65.7%) were infected with both wCI and wFem (CF) on Tanegashima Island. While females singly infected with wCI (C) produced offspring at a 1 : 1 sex ratio, CF females produced only females. Although Z-linked sequence polymorphism showed no signs of divergence between C and CF females, mtDNA split into two discrete clades; one consisted of C females and the other CF females, both of which formed a clade with E. hecabe but not with uninfected E. mandarina This suggests that CF matrilines also, but independently, experienced a selective sweep after hybrid introgression from E. hecabe Distinct evolutionary forces were suggested to have caused C and CF matrilines to diverge, which would be irreversible because of the particular phenotype of wFem.

Relative Abundance and Strain Diversity in the Bacterial Endosymbiont Community of a Sap-Feeding Insect Across Its Native and Introduced Geographic Range

Most insects are associated with bacterial symbionts. The bacterial diversity and community composition within hosts may play an important role in shaping insect population biology, ecology and evolution. We focussed on the bacterial microbiome of the Australian fig homotomid Mycopsylla fici (Hemiptera: Psylloidea), which can cause defoliation of its only host tree, Ficus macrophylla. This sap-feeding insect is native to mainland Australia and Lord Howe Island (LHI) but also occurs where its host has been planted, notably in New Zealand. By using a high-throughput 16S rDNA amplicon sequencing approach, we compared the bacterial diversity and community composition in individual adult males of four host populations, Sydney, Brisbane, LHI and Auckland. We also compared males, females and nymphs of the Sydney population. The microbiome of M. fici was simple and consisted mostly of the following three maternally inherited endosymbiont species: the primary endosymbiont Carsonella, a secondary (S-) endosymbiont and Wolbachia. However, the relative abundance of their sequence reads varied between host populations, except for similarities between Sydney and Auckland. In addition, insects from Sydney and Auckland had identical bacterial strains supporting the hypothesis that Sydney is the source population for Auckland. In contrast, mainland and LHI populations harboured the same S-endosymbiont, co-diverged Carsonella but different Wolbachia strains. Besides detecting endosymbiont-specific patterns of either co-evolution or horizontal acquisition, our study highlights that relative abundance of maternally inherited endosymbionts should also be taken into account when studying bacterial communities across host populations, as variations in bacterial density may impact host biology and ecology.

Wolbachia in Parasitoids Attacking Native European and Introduced Eastern Cherry Fruit Flies in Europe

The eastern cherry fruit fly, Rhagoletis cingulata Loew (Diptera: Tephritidae), is an economically important pest of cherries in North America. In 1983 it was first reported in Europe where it shares its ecological niche with the native European cherry fruit fly, Rhagoletis cerasi L. (Diptera: Tephritidae). Their coexistence in Europe led to the recent horizontal transmission of the Wolbachia strain wCer1 from R. cerasi to R. cingulata Horizontal Wolbachia transmission is mediated by either sharing of ecological niches or by interacting species such as parasitoids. Here we describe for the first time that two braconid wasps, Psyttalia rhagoleticola Sachtleben (Hymenoptera: Braconidae) and Utetes magnus Fischer (Hymenoptera: Braconidae), naturally parasitizing R. cerasi, use the invasive R. cingulata in Europe as a new host. In contrast, no parasitoids that parasitize R. cingulata in its native American range were detected in the introduced European range. Diagnostic Wolbachia PCR screening and sequence analyses demonstrated that all P. rhagoleticola individuals were infected with the newly described Wolbachia strain wRha while all U. magnus individuals were uninfected. wRha is different from wCer1 but had an Wolbachia surface protein (wsp) gene sequence that was identical to wCer2 of R. cerasi and wCin2 of R. cingulata. However, multi locus sequence typing revealed differences in all loci between wRha and the tephritid's strains. The horizontal transmission of wCer1 between the two tephritid species did not result in fixed heritable infections in the parasitoids. However, the parasitoids may have acted as a transient wCer1 vector.

Phylogeographic analyses of bacterial endosymbionts in fig homotomids (Hemiptera: Psylloidea) reveal codiversification of both primary and secondary endosymbionts

While obligate primary (P-) endosymbionts usually cospeciate with their insect hosts, less is known about codiversification of secondary (S-) endosymbionts that are generally considered facultative. Typically, insects of the superfamily Psylloidea harbour one P- (Carsonella) and at least one S-endosymbiont, thought to compensate for Carsonella genome reduction. Most codiversification studies have used phylogenies of psyllids and their endosymbionts across and within host families or genera, but few have explored patterns within species. We focussed on P- and S-endosymbionts of three Mycopsylla (Homotomidae) species to explore whether they have congruent phylogenies and within-species geographic structures. The P-endosymbiont Carsonella, a S-endosymbiont and Wolbachia all had 100% prevalence, while Arsenophonus was only found in one species at low prevalence. Congruent phylogenies of Mycopsylla and P-endosymbionts across populations and species support strict cospeciation. S-endosymbiont phylogenies were also congruent across host species but low genetic variation in the S-endosymbiont was not correlated with host phylogeography, possibly due to a shorter evolutionary association. Between species, Wolbachia and Mycopsylla phylogenies were incongruent, probably due to horizontal transmission events. Our study is the first to explore endosymbionts of Mycopsylla and further supports the codivergence of Psylloidea hosts and P-endosymbionts, with obligate host interactions for both P- and S-endosymbionts.