Distribution and Evolution of the Bacteriophage WO and Its Antagonism With Wolbachia

Both symbionts and environmental factors contribute to shape the microbiota in a pest insect, Sogatella furcifera

Introduction

Bacterial symbionts are prevalent in arthropods globally and play a vital role in the fitness and resistance of hosts. While several symbiont infections have been identified in the white-backed planthopper Sogatella furcifera, the impact of environmental factors on the microbiota within S. furcifera remains elusive.

Methods

In this study, a total of 142 S. furcifera individuals from 18 populations were collected from 14 locations across six countries (China, Thailand, Myanmar, Cambodia, Vietnam, and Laos) analyzed with 2bRAD-M sequencing, to examine the effects of symbionts on the microbiota in the S. furcifera population, as well as the vital effects of environmental factors on the bacterial communities.

Results and discussion

Based on the results, in S. furcifera, the presence of symbionts Wolbachia and Cardinium negatively influenced the abundance of other bacteria, including Enterobacter, Acinetobacter, and Lysinibacillus, while Wolbachia infection significantly decreased the diversity of the microbial community. Moreover, several environmental factors, including longitude, latitude, temperature, and precipitation, affected the abundance of symbionts and microbiota diversity in S. furcifera. These results collectively highlight the vital role of Wolbachia in S. furcifera microbiota, as well as the intricate effects of environmental factors on the bacterial communities of S. furcifera.

Evolution of Wolbachia reproductive and nutritional mutualism: insights from the genomes of two novel strains that double infect the pollinator of dioecious Ficus hirta

Wolbachia is a genus of maternally inherited endosymbionts that can affect reproduction of their hosts and influence metabolic processes. The pollinator, Valisia javana, is common in the male syconium of the dioecious fig Ficus hirta. Based on a high-quality chromosome-level V. javana genome with PacBio long-read and Illumina short-read sequencing, we discovered a sizeable proportion of Wolbachia sequences and used these to assemble two novel Wolbachia strains belonging to supergroup A. We explored its phylogenetic relationship with described Wolbachia strains based on MLST sequences and the possibility of induction of CI (cytoplasmic incompatibility) in this strain by examining the presence of cif genes known to be responsible for CI in other insects. We also identified mobile genetic elements including prophages and insertion sequences, genes related to biotin synthesis and metabolism. A total of two prophages and 256 insertion sequences were found. The prophage WOjav1 is cryptic (structure incomplete) and WOjav2 is relatively intact. IS5 is the dominant transposon family. At least three pairs of type I cif genes with three copies were found which may cause strong CI although this needs experimental verification; we also considered possible nutritional effects of the Wolbachia by identifying genes related to biotin production, absorption and metabolism. This study provides a resource for further studies of Wolbachia-pollinator-host plant interactions.

Phylogenomic analysis of Wolbachia genomes from the Darwin Tree of Life biodiversity genomics project

The Darwin Tree of Life (DToL) project aims to sequence all described terrestrial and aquatic eukaryotic species found in Britain and Ireland. Reference genome sequences are generated from single individuals for each target species. In addition to the target genome, sequenced samples often contain genetic material from microbiomes, endosymbionts, parasites, and other cobionts. Wolbachia endosymbiotic bacteria are found in a diversity of terrestrial arthropods and nematodes, with supergroups A and B the most common in insects. We identified and assembled 110 complete Wolbachia genomes from 93 host species spanning 92 families by filtering data from 368 insect species generated by the DToL project. From 15 infected species, we assembled more than one Wolbachia genome, including cases where individuals carried simultaneous supergroup A and B infections. Different insect orders had distinct patterns of infection, with Lepidopteran hosts mostly infected with supergroup B, while infections in Diptera and Hymenoptera were dominated by A-type Wolbachia. Other than these large-scale order-level associations, host and Wolbachia phylogenies revealed no (or very limited) cophylogeny. This points to the occurrence of frequent host switching events, including between insect orders, in the evolutionary history of the Wolbachia pandemic. While supergroup A and B genomes had distinct GC% and GC skew, and B genomes had a larger core gene set and tended to be longer, it was the abundance of copies of bacteriophage WO who was a strong determinant of Wolbachia genome size. Mining raw genome data generated for reference genome assemblies is a robust way of identifying and analysing cobiont genomes and giving greater ecological context for their hosts.

Muramidase, nuclease, or hypothetical protein genes intervene between paired genes encoding DNA packaging terminase and portal proteins in Wolbachia phages and prophages

Genomes of the obligate intracellular alpha proteobacterium Wolbachia pipientis often encode prophage-like regions, and in a few cases, purified particles have been recovered. Because the structure of a conserved WO phage genome has been difficult to establish, we examined paired terminase and portal genes in Wolbachia phages and prophages, relative to those encoded by the gene transfer agent RcGTA from the free-living alpha proteobacterium Rhodobacter capsulatus. Terminase and portal proteins from Wolbachia have higher similarity to orthologs encoded by RcGTA than to orthologs encoded by bacteriophage lambda. In lambdoid phages, these proteins play key roles in assembly of mature phage particles, while in less well-studied gene transfer agents, terminase and portal proteins package random fragments of bacterial DNA, which could confound elucidation of WO phage genomes. In WO phages and prophages, terminase genes followed by a short gpW gene may be separated from the downstream portal gene by open-reading frames encoding a GH_25 hydrolase/muramidase, a PD-(D/E)XK nuclease, a hypothetical protein and/or a RelE/ParE toxin-antitoxin module. These aspects of gene organization, coupled with evidence for a low, non-inducible yield of WO phages, and the small size of WO phage particles described in the literature raise the possibility that Wolbachia prophage regions participate in processes that extend beyond conventional bacteriophage lysogeny and lytic replication. These intervening genes, and their possible relation to functions associated with GTAs, may contribute to variability among WO phage genomes recovered from physical particles and impact the ability of WO phages to act as transducing agents.

Wolbachia endosymbionts in two Anopheles species indicates independent acquisitions and lack of prophage elements

Wolbachia is a genus of obligate bacterial endosymbionts that infect a diverse range of arthropod species as well as filarial nematodes, with its single described species, Wolbachia pipientis, divided into several ‘supergroups’ based on multilocus sequence typing. Wolbachia strains in mosquitoes have been shown to inhibit the transmission of human pathogens, including Plasmodium malaria parasites and arboviruses. Despite their large host range, Wolbachia strains within the major malaria vectors of the Anopheles gambiae and Anopheles funestus complexes appear at low density, established solely on PCR-based methods. Questions have been raised as to whether this represents a true endosymbiotic relationship. However, recent definitive evidence for two distinct, high-density strains of supergroup B Wolbachia within Anopheles demeilloni and Anopheles moucheti has opened exciting possibilities to explore naturally occurring Wolbachia endosymbionts in Anopheles for biocontrol strategies to block Plasmodium transmission. Here, we utilize genomic analyses to demonstrate that both Wolbachia strains have retained all key metabolic and transport pathways despite their smaller genome size, with this reduction potentially attributable to degenerated prophage regions. Even with this reduction, we confirmed the presence of cytoplasmic incompatibility (CI) factor genes within both strains, with wAnD maintaining intact copies of these genes while the cifB gene was interrupted in wAnM, so functional analysis is required to determine whether wAnM can induce CI. Additionally, phylogenetic analysis indicates that these Wolbachia strains may have been introduced into these two Anopheles species via horizontal transmission events, rather than by ancestral acquisition and subsequent loss events in the Anopheles gambiae species complex. These are the first Wolbachia genomes, to our knowledge, that enable us to study the relationship between natural strain Plasmodium malaria parasites and their anopheline hosts.

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.

The impact of artificial selection for Wolbachia-mediated dengue virus blocking on phage WO

Wolbachia is currently at the forefront of global efforts to control arbovirus transmission from the vector Aedes aegypti. The use of Wolbachia relies on two phenotypes—cytoplasmic incompatibility (CI), conferred by cifA and cifB genes in prophage WO, and Wolbachia-mediated pathogen blocking (WMPB). These traits allow for local, self-sustaining reductions in transmission of dengue (DENV) following release of Wolbachia-infected A. aegypti. Here, aided by previous artificial selection experiment that generated Low and High pathogen blocking lines, we examined the potential link between WMPB and phage WO. We found no evidence that Wolbachia or phage WO relative densities predict DENV blocking strength across selected lines. However, selection resulted in reduced phage WO relative density for the Low WMPB line. The Low blocking line was previously shown to have reduced fitness as a result of selection. Through subsequent genomic analyses, we demonstrate that SNP variation underpinning selection for low blocking led to elevated frequency of potential deleterious SNPs on chromosome 1. The key region on chromosome 1 contains genes relating to cell cycle regulation, oxidative stress, transcriptional pausing, among others, that may have cascading effects on Wolbachia intracellular environment. We hypothesize that reduction in phage WO may be driven by changes in the loci directly under selection for blocking, or by the accumulation of predicted deleterious alleles in linkage disequilibrium with blocking loci resulting from hitchhiking. For the Low line with fewer phage WO, we also detected reduced expression of cifA and cifB CI genes, with patterns of expression varying between somatic and reproductive tissues. In conclusion, we propose that artificial selection for WMPB trait had corresponding impacts on phage WO densities, and also the transcription of CI-causing genes. Future studies may include a more detailed analysis of the regions the A. aegypti chromosome 1’s ability to affect WMPB and other Wolbachia-associated intrinsic factors such as phage WO.

Wolbachia are widespread endosymbiotic bacteria of insects that cause Wolbachia-mediated pathogen blocking (WMPB) and cytoplasmic incompatibility (CI). The latter mediated by cif genes localized in the prophage WO region. Because of that, Wolbachia-infected mosquitoes are currently being used in field to fight the transmission of vector-borne viruses such as Dengue (DENV) to human populations. Aided by a previous artificial selection experiment that generated lines with variable (High and Low) DENV blocking strength, we tested for a potential link between WMPB and phage WO. There was no evidence that Wolbachia nor phage WO densities predict DENV blocking strength. However, we found that the Low blocking line had reduced phage WO density, and lower expression of the cif genes in a tissue-specific manner. We demonstrate that in addition to previous report of reduced fitness, the Low blocking line also exhibited increased frequency of potential deleterious SNPs on chromosome 1. Our hypotheses are that reduction in phage WO may have resulted from changes in the loci directly under selection for blocking, or by linkage disequilibrium events linked to the accumulation of mosquito predicted deleterious alleles. Our results highlight the importance of chromosome 1 for WMPB and its potential impact for other Wolbachia-associated factors like phage WO.

Wolbachia as translational science: controlling mosquito-borne pathogens

In this review we examine how exploiting the Wolbachia-mosquito relationship has become an increasingly popular strategy for controlling arbovirus transmission. Field deployments of Wolbachia-infected mosquitoes have led to significant decreases in dengue virus incidence via high levels of mosquito population suppression and replacement, emphasizing the success of Wolbachia approaches. Here, we examine how improved knowledge of Wolbachia-host interactions has provided key insight into the mechanisms of the essential phenotypes of pathogen blocking and cytoplasmic incompatibility. And we discuss recent studies demonstrating that extrinsic factors, such as ambient temperature, can modulate Wolbachia density and maternal transmission. Finally, we assess the prospects of using Wolbachia to control other vectors and agricultural pest species.

Bacteriophage-Host Association in the Phytoplasma Insect Vector Euscelidius variegatus

Insect vectors transmit viruses and bacteria that can cause severe diseases in plants and economic losses due to a decrease in crop production. Insect vectors, like all other organisms, are colonized by a community of various microorganisms, which can influence their physiology, ecology, evolution, and also their competence as vectors. The important ecological meaning of bacteriophages in various ecosystems and their role in microbial communities has emerged in the past decade. However, only a few phages have been described so far in insect microbiomes. The leafhopper Euscelidius variegatus is a laboratory vector of the phytoplasma causing Flavescence dorée, a severe grapevine disease that threatens viticulture in Europe. Here, the presence of a temperate bacteriophage in E. variegatus (named Euscelidius variegatus phage 1, EVP-1) was revealed through both insect transcriptome analyses and electron microscopic observations. The bacterial host was isolated in axenic culture and identified as the bacterial endosymbiont of E. variegatus (BEV), recently assigned to the genus Candidatus Symbiopectobacterium. BEV harbors multiple prophages that become active in culture, suggesting that different environments can trigger different mechanisms, finely regulating the interactions among phages. Understanding the complex relationships within insect vector microbiomes may help in revealing possible microbe influences on pathogen transmission, and it is a crucial step toward innovative sustainable strategies for disease management in agriculture.

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.