Phylogenomic Analysis of Wolbachia Strains Reveals Patterns of Genome Evolution and Recombination

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

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

Mixed Wolbachia infections resolve rapidly during in vitro evolution

The intracellular symbiont Wolbachia pipientis evolved after the divergence of arthropods and nematodes, but it reached high prevalence in many of these taxa through its abilities to infect new hosts and their germlines. Some strains exhibit long-term patterns of co-evolution with their hosts, while other strains are capable of switching hosts. This makes strain selection an important factor in symbiont-based biological control. However, little is known about the ecological and evolutionary interactions that occur when a promiscuous strain colonizes an infected host. Here, we study what occurs when two strains come into contact in host cells following horizontal transmission and infection. We focus on the faithful wMel strain from Drosophila melanogaster and the promiscuous wRi strain from Drosophila simulans using an in vitro cell culture system with multiple host cell types and combinatorial infection states. Mixing D. melanogaster cell lines stably infected with wMel and wRi revealed that wMel outcompetes wRi quickly and reproducibly. Furthermore, wMel was able to competitively exclude wRi even from minuscule starting quantities, indicating that this is a nearly deterministic outcome, independent of the starting infection frequency. This competitive advantage was not exclusive to wMel's native D. melanogaster cell background, as wMel also outgrew wRi in D. simulans cells. Overall, wRi is less adept at in vitro growth and survival than wMel and its in vivo state, revealing differences between cellular and humoral regulation. These attributes may underlie the observed low rate of mixed infections in nature and the relatively rare rate of host-switching in most strains. Our in vitro experimental framework for estimating cellular growth dynamics of Wolbachia strains in different host species, tissues, and cell types provides the first strategy for parameterizing endosymbiont and host cell biology at high resolution. This toolset will be crucial to our application of these bacteria as biological control agents in novel hosts and ecosystems.

Wolbachia supergroup A in Enoplognatha latimana (Araneae: Theridiidae) in Poland as an example of possible horizontal transfer of bacteria

Wolbachia (phylum Pseudomonadota, class Alfaproteobacteria, order Rickettsiales, family Ehrlichiaceae) is a maternally inherited bacterial symbiont infecting more than half of arthropod species worldwide and constituting an important force in the evolution, biology, and ecology of invertebrate hosts. Our study contributes to the limited knowledge regarding the presence of intracellular symbiotic bacteria in spiders. Specifically, we investigated the occurrence of Wolbachia infection in the spider species Enoplognatha latimana Hippa and Oksala, 1982 (Araneae: Theridiidae) using a sample collected in north-western Poland. To the best of our knowledge, this is the first report of Wolbachia infection in E. latimana. A phylogeny based on the sequence analysis of multiple genes, including 16S rRNA, coxA, fbpA, ftsZ, gatB, gltA, groEL, hcpA, and wsp revealed that Wolbachia from the spider represented supergroup A and was related to bacterial endosymbionts discovered in other spider hosts, as well as insects of the orders Diptera and Hymenoptera. A sequence unique for Wolbachia supergroup A was detected for the ftsZ gene. The sequences of Wolbachia housekeeping genes have been deposited in publicly available databases and are an important source of molecular data for comparative studies. The etiology of Wolbachia infection in E. latimana is discussed.

Detection and quantification of natural Wolbachia in Aedes aegypti in Metropolitan Manila, Philippines using locally designed primers

Background

The Philippines bears health and economic burden caused by high dengue cases annually. Presently, the Philippines still lack an effective and sustainable vector management. The use of Wolbachia, a maternally transmitted bacterium, that mitigate arbovirus transmission has been recommended. Cytoplasmic incompatibility and viral blocking, two characteristics that make Wolbachia suitable for vector control, depend on infection prevalence and density. There are no current Wolbachia release programs in the Philippines, and studies regarding the safety of this intervention. Here, we screened for Wolbachia in Aedes aegypti collected from Metropolitan Manila, Philippines. We designed location-specific primers for qPCR to test whether this improved Wolbachia detection in Ae. aegypti. We explored if host sex and Wolbachia strain could be potential factors affecting Wolbachia density.

Methods

Ae. aegypti mosquitoes (n=429) were screened for natural Wolbachia by taqman qPCR using location-specific Wolbachia surface protein primers (wspAAML) and known 16S rRNA primers. Samples positive for wspAAML (n=267) were processed for Sanger sequencing. We constructed a phylogenetic tree using IQ-TREE 2 to further characterize Wolbachia present in the Philippine Ae. aegypti. We then compared Wolbachia densities between Wolbachia groups and host sex. Statistical analyses were done using GraphPad Prism 9.0.

Results

Wolbachia prevalence for 16S rRNA (40%) and wspAAML (62%) markers were high. Wolbachia relative densities for 16S rRNA ranged from -3.84 to 2.71 and wspAAML from -4.02 to 1.81. Densities were higher in male than female mosquitoes. Wolbachia strains detected in Ae. aegypti clustered into supergroup B. Some 54% (123/226) of these sequences clustered under a group referred to here as "wAegML," that belongs to the supergroup B, which had a significantly lower density than wAegB/wAlbB, and wAlbA strains.

Conclusion

Location-specific primers improved detection of natural Wolbachia in Ae. aegypti and allowed for relative quantification. Wolbachia density is relatively low, and differed between host sexes and Wolbachia strains. An economical way of confirming sporadic or transient Wolbachia in Ae. aegypti is necessary while considering host sex and bacterial strain.

Bacterial community and genome analysis of cytoplasmic incompatibility-inducing Wolbachia in American serpentine leafminer, Liriomyza trifolii

Liriomyza trifolii, an agricultural pest, is occasionally infected by Wolbachia. A Wolbachia strain present in Liriomyza trifolii is associated with cytoplasmic incompatibility (CI) effects, leading to the death of embryos resulting from incompatible crosses between antibiotic-treated or naturally Wolbachia-free strain females and Wolbachia-infected males. In this study, high-throughput sequencing of hypervariable rRNA genes was employed to characterize the bacterial community in Wolbachia-infected L. trifolii without antibiotic treatment. The analysis revealed that Wolbachia dominates the bacterial community in L. trifolii, with minor presence of Acinetobacter, Pseudomonas, and Limnobacter. To elucidate the genetic basis of the CI phenotype, metagenomic sequencing was also conducted to assemble the genome of the Wolbachia strain. The draft-genome of the Wolbachia strain wLtri was 1.35 Mbp with 34% GC content and contained 1,487 predicted genes. Notably, within the wLtri genome, there are three distinct types of cytoplasmic incompatibility factor (cif) genes: Type I, Type III, and Type V cifA;B. These genes are likely responsible for inducing the strong cytoplasmic incompatibility observed in L. trifolii.

Population Genomics of the Mostly Thelytokous Diplolepis rosae (Linnaeus, 1758) (Hymenoptera: Cynipidae) Reveals Population-specific Selection for Sex

In Hymenoptera, arrhenotokous parthenogenesis (arrhenotoky) is a common reproductive mode. Thelytokous parthenogenesis (thelytoky), when virgin females produce only females, is less common and is found in several taxa. In our study, we assessed the efficacy of recombination and the effect of thelytoky on the genome structure of Diplolepis rosae, a gall wasp-producing bedeguars in dog roses. We assembled a high-quality reference genome using Oxford Nanopore long-read technology and sequenced 17 samples collected in France with high-coverage Illumina reads. We found two D. rosae peripatric lineages that differed in the level of recombination and homozygosity. One of the D. rosae lineages showed a recombination rate that was 13.2 times higher and per-individual heterozygosity that was 1.6 times higher. In the more recombining lineage, the genes enriched in functions related to male traits ('sperm competition", "insemination", and "copulation" gene ontology terms) showed signals of purifying selection, whereas in the less recombining lineage, the same genes showed traces pointing towards balancing or relaxed selection. Thus, although D. rosae reproduces mainly by thelytoky, selection may act to maintain sexual reproduction.

Recombination in Bacterial Genomes: Evolutionary Trends

Bacterial organisms have undergone homologous recombination (HR) and horizontal gene transfer (HGT) multiple times during their history. These processes could increase fitness to new environments, cause specialization, the emergence of new species, and changes in virulence. Therefore, comprehensive knowledge of the impact and intensity of genetic exchanges and the location of recombination hotspots on the genome is necessary for understanding the dynamics of adaptation to various conditions. To this end, we aimed to characterize the functional impact and genomic context of computationally detected recombination events by analyzing genomic studies of any bacterial species, for which events have been detected in the last 30 years. Genomic loci where the transfer of DNA was detected pertained to mobile genetic elements (MGEs) housing genes that code for proteins engaged in distinct cellular processes, such as secretion systems, toxins, infection effectors, biosynthesis enzymes, etc. We found that all inferences fall into three main lifestyle categories, namely, ecological diversification, pathogenesis, and symbiosis. The latter primarily exhibits ancestral events, thus, possibly indicating that adaptation appears to be governed by similar recombination-dependent mechanisms.

The prevalence of Wolbachia in multiple cockroach species and its implication for urban insect management

Cockroach management relies heavily on the use of conventional insecticides in urban settings, which no longer provide the anticipated level of control. Knowledge of cockroach endosymbionts, like Wolbachia, might provide novel avenues for control. Therefore, we screened 16 cockroach species belonging to 3 families (Ectobiidae, Blattidae, and Blaberidae) for the presence of Wolbachia. We mapped the evolution of Wolbachia-cockroach relationships based on maximum likelihood phylogeny and phylogenetic species clustering on a multi-loci sequence dataset (i.e., coxA, virD4, hcpA, and gatB) of Wolbachia genes. We confirmed the previous report of Wolbachia in 1 Ectobiid species; Supella longipalpa (Fab.), and detected the presence of Wolbachia in 2 Ectobiid species; Balta notulata (Stål) and Pseudomops septentrionalis Hebard, and 1 Blaberid species; Gromphadorhina portentosa (Schaum). All cockroach-associated Wolbachia herein detected were clustered with the ancestor of F clade Wolbachia of Cimex lectularius L. (bed bugs). Since Wolbachia provision C. lectularius with biotin vitamins that confer reproductive fitness, we screened the cockroach-associated Wolbachia for the presence of biotin genes. In toto, our results reveal 2 important findings: (i) Wolbachia is relatively uncommon among cockroach species infecting about 25% of species investigated, and (ii) cockroach-associated Wolbachia have biotin genes that likely provide nutritional benefits to their hosts. Thus, we discuss the potential of exploring Wolbachia as a tool for urban insect management.

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.

Host phylogeny and ecological associations best explain Wolbachia host shifts in scale insects

Wolbachia are among the most prevalent and widespread endosymbiotic bacteria on Earth. Wolbachia's success in infecting an enormous number of arthropod species is attributed to two features: the range of phenotypes they induce in their hosts, and their ability to switch between host species. Whilst much progress has been made in elucidating their induced phenotypes, our understanding of Wolbachia host-shifting is still very limited: we lack answers to even fundamental questions concerning Wolbachia's routes of transfer and the importance of factors influencing host shifts. Here, we investigate the diversity and host-shift patterns of Wolbachia in scale insects, a group of arthropods with intimate associations with other insects that make them well suited to studying host shifts. Using Illumina multitarget amplicon sequencing of Wolbachia-infected scale insects and their direct associates we determined the identity of all Wolbachia strains. We then fitted a generalized additive mixed model to our data to estimate the influence of host phylogeny and the geographical distribution on Wolbachia strain sharing among scale insect species. The model predicts no significant contribution of host geography but strong effects of host phylogeny, with high rates of Wolbachia sharing among closely related species and a sudden drop-off in sharing with increasing phylogenetic distance. We also detected the same Wolbachia strain in scale insects and several intimately associated species (ants, wasps and flies). This indicates putative host shifts and potential routes of transfers via these associates and highlights the importance of ecological connectivity in Wolbachia host-shifting.

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.

Narrow Genetic Diversity of Wolbachia Symbionts in Acrididae Grasshopper Hosts (Insecta, Orthoptera)

Bacteria of the Wolbachia genus are maternally inherited symbionts of Nematoda and numerous Arthropoda hosts. There are approximately 20 lineages of Wolbachia, which are called supergroups, and they are designated alphabetically. Wolbachia strains of the supergroups A and B are predominant in arthropods, especially in insects, and supergroup F seems to rank third. Host taxa have been studied very unevenly for Wolbachia symbionts, and here, we turn to one of largely unexplored insect families: Acrididae. On the basis of five genes subject to multilocus sequence typing, we investigated the incidence and genetic diversity of Wolbachia in 41 species belonging three subfamilies (Gomphocerinae, Oedipodinae, and Podisminae) collected in Turkey, Kazakhstan, Tajikistan, Russia, and Japan, making 501 specimens in total. Our results revealed a high incidence and very narrow genetic diversity of Wolbachia. Although only the strains belonging to supergroups A and B are commonly present in present, the Acrididae hosts here proved to be infected with supergroups B and F without A-supergroup variants. The only trace of an A-supergroup lineage was noted in one case of an inter-supergroup recombinant haplotype, where the ftsZ gene came from supergroup A, and the others from supergroup B. Variation in the Wolbachia haplotypes in Acrididae hosts within supergroups B and F was extremely low. A comprehensive genetic analysis of Wolbachia diversity confirmed specific features of the Wolbachia allelic set in Acrididae hosts. This result can help to elucidate the crucial issue of Wolbachia biology: the route(s) and mechanism(s) of Wolbachia horizontal transmission.

Complete de novo assembly of Wolbachia endosymbiont of Diaphorina citri Kuwayama (Hemiptera: Liviidae) using long-read genome sequencing

Wolbachia, a gram-negative \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{\alpha }$$\end{document}α-proteobacterium, is an endosymbiont found in some arthropods and nematodes. Diaphorina citri Kuwayama, the vector of ‘Candidatus Liberibacter asiaticus’ (CLas), are naturally infected with a strain of Wolbachia (wDi), which has been shown to colocalize with the bacteria pathogens CLas, the pathogen associated with huanglongbing (HLB) disease of citrus. The relationship between wDi and CLas is poorly understood in part because the complete genome of wDi has not been available. Using high-quality long-read PacBio circular consensus sequences, we present the largest complete circular wDi genome among supergroup-B members. The assembled circular chromosome is 1.52 megabases with 95.7% genome completeness with contamination of 1.45%, as assessed by checkM. We identified Insertion Sequences (ISs) and prophage genes scattered throughout the genomes. The proteins were annotated using Pfam, eggNOG, and COG that assigned unique domains and functions. The wDi genome was compared with previously sequenced Wolbachia genomes using pangenome and phylogenetic analyses. The availability of a complete circular chromosome of wDi will facilitate understanding of its role within the insect vector, which may assist in developing tools for disease management. This information also provides a baseline for understanding phylogenetic relationships among Wolbachia of other insect vectors.

Are ecological communities the seat of endosymbiont horizontal transfer and diversification? A case study with soil arthropod community

Maternally inherited endosymbionts of arthropods are one of the most abundant and diverse group of bacteria. These bacterial endosymbionts also show extensive horizontal transfer to taxonomically unrelated hosts and widespread recombination in their genomes. Such horizontal transfers can be enhanced when different arthropod hosts come in contact like in an ecological community. Higher rates of horizontal transfer can also increase the probability of recombination between endosymbionts, as they now share the same host cytoplasm. However, reports of community-wide endosymbiont data are rare as most studies choose few host taxa and specific ecological interactions among the hosts. To better understand endosymbiont spread within host populations, we investigated the incidence, diversity, extent of horizontal transfer, and recombination of three endosymbionts (Wolbachia, Cardinium, and Arsenophonus) in a specific soil arthropod community. Wolbachia strains were characterized with MLST genes whereas 16S rRNA gene was used for Cardinium and Arsenophonus. Among 3,509 individual host arthropods, belonging to 390 morphospecies, 12.05% were infected with Wolbachia, 2.82% with Cardinium and 2.05% with Arsenophonus. Phylogenetic incongruence between host and endosymbiont indicated extensive horizontal transfer of endosymbionts within this community. Three cases of recombination between Wolbachia supergroups and eight incidences of within-supergroup recombination were also found. Statistical tests of similarity indicated supergroup A Wolbachia and Cardinium show a pattern consistent with extensive horizontal transfer within the community but not for supergroup B Wolbachia and Arsenophonus. We highlight the importance of extensive community-wide studies for a better understanding of the spread of endosymbionts across global arthropod communities.

Comparative analysis reveals the expansion of mitochondrial DNA control region containing unusually high G-C tandem repeat arrays in Nasonia vitripennis

Insect mitochondrial DNA (mtDNA) ranges from 14 to 19 kbp, and the size difference is attributed to the AT-rich control region. Jewel wasps have a parasitoid lifestyle, which may affect mitochondria function and evolution. We sequenced, assembled, and annotated mitochondrial genomes in Nasonia and outgroup species. Gene composition and order are conserved within Nasonia, but they differ from other parasitoids by two large inversion events that were not reported before. We observed a much higher substitution rate relative to the nuclear genome and mitochondrial introgression between N. giraulti and N. oneida, which is consistent with previous studies. Most strikingly, N. vitripennis mtDNA has an extremely long control region (7665 bp), containing twenty-nine 217 bp tandem repeats and can fold into a super-cruciform structure. In contrast to tandem repeats commonly found in other mitochondria, these high-copy repeats are highly conserved (98.7% sequence identity), much longer in length (approximately 8 Kb), extremely GC-rich (50.7%), and CpG-rich (percent CpG 19.4% vs. 1.1% in coding region), resulting in a 23 kbp mtDNA beyond the typical size range in insects. These N. vitripennis-specific mitochondrial repeats are not related to any known sequences in insect mitochondria. Their evolutionary origin and functional consequences warrant further investigations.