Widespread phages of endosymbionts: Phage WO genomics and the proposed taxonomic classification of Symbioviridae.

Wolbachia are the most common obligate, intracellular bacteria in animals. They exist worldwide in arthropod and nematode hosts in which they commonly act as reproductive parasites or mutualists, respectively. Bacteriophage WO, the largest of Wolbachia's mobile elements, includes reproductive parasitism genes, serves as a hotspot for genetic divergence and genomic rearrangement of the bacterial chromosome, and uniquely encodes a Eukaryotic Association Module with eukaryotic-like genes and an ensemble of putative host interaction genes. Despite WO's relevance to genome evolution, selfish genetics, and symbiotic applications, relatively little is known about its origin, host range, diversification, and taxonomic classification. Here we analyze the most comprehensive set of 150 Wolbachia and phage WO assemblies to provide a framework for discretely organizing and naming integrated phage WO genomes. We demonstrate that WO is principally in arthropod Wolbachia with relatives in diverse endosymbionts and metagenomes, organized into four variants related by gene synteny, often oriented opposite the putative origin of replication in the Wolbachia chromosome, and the large serine recombinase is an ideal typing tool to distinguish the four variants. We identify a novel, putative lytic cassette and WO's association with a conserved eleven gene island, termed Undecim Cluster, that is enriched with virulence-like genes. Finally, we evaluate WO-like Islands in the Wolbachia genome and discuss a new model in which Octomom, a notable WO-like Island, arose from a split with WO. Together, these findings establish the first comprehensive Linnaean taxonomic classification of endosymbiont phages, including non-Wolbachia phages from aquatic environments, that includes a new family and two new genera to capture the collective relatedness of these viruses.

Wolbachia 16S rRNA haplotypes detected in wild Anopheles stephensi in eastern Ethiopia.

BACKGROUND: About two out of three Ethiopians are at risk of malaria, a disease caused by the parasites Plasmodium falciparum and Plasmodium vivax. Anopheles stephensi, an invasive vector typically found in South Asia and the Middle East, was recently found to be distributed across eastern and central Ethiopia and is capable of transmitting both P. falciparum and P. vivax. The detection of this vector in the Horn of Africa (HOA) coupled with widespread insecticide resistance requires that new methods of vector control be investigated in order to control the spread of malaria. Wolbachia, a naturally occurring endosymbiotic bacterium of mosquitoes, has been identified as a potential vector control tool that can be explored for the control of malaria transmission. Wolbachia could be used to control the mosquito population through suppression or potentially decrease malaria transmission through population replacement. However, the presence of Wolbachia in wild An. stephensi in eastern Ethiopia is unknown. This study aimed to identify the presence and diversity of Wolbachia in An. stephensi across eastern Ethiopia. METHODS: DNA was extracted from An. stephensi collected from eastern Ethiopia in 2018 and screened for Wolbachia using a 16S targeted PCR assay, as well as multilocus strain typing (MLST) PCR assays. Haplotype and phylogenetic analysis of the sequenced 16S amplicons were conducted to compare with Wolbachia from countries across Africa and Asia. RESULTS: Twenty out of the 184 mosquitoes screened were positive for Wolbachia, with multiple haplotypes detected. In addition, phylogenetic analysis revealed two superclades, representing Wolbachia supergroups A and B (bootstrap values of 81 and 72, respectively) with no significant grouping of geographic location or species. A subclade with a bootstrap value of 89 separates the Ethiopian haplotype 2 from other sequences in that superclade. CONCLUSIONS: These findings provide the first evidence of natural Wolbachia populations in wild An. stephensi in the HOA. They also identify the need for further research to confirm the endosymbiotic relationship between Wolbachia and An. stephensi and to investigate its utility for malaria control in the HOA.

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.

Living in the endosymbiotic world of Wolbachia: A centennial review.

The most widespread intracellular bacteria in the animal kingdom are maternally inherited endosymbionts of the genus Wolbachia. Their prevalence in arthropods and nematodes worldwide and stunning arsenal of parasitic and mutualistic adaptations make these bacteria a biological archetype for basic studies of symbiosis and applied outcomes for curbing human and agricultural diseases. Here, we conduct a summative, centennial analysis of living in the Wolbachia world. We synthesize literature on Wolbachia's host range, phylogenetic diversity, genomics, cell biology, and applications to filarial, arboviral, and agricultural diseases. We also review the mobilome of Wolbachia including phage WO and its essentiality to hallmark reproductive phenotypes in arthropods. Finally, the Wolbachia system is an exemplar for discovery-based science education using biodiversity, biotechnology, and bioinformatics lessons. As we approach a century of Wolbachia research, the interdisciplinary science of this symbiosis stands as a model for consolidating and teaching the integrative rules of endosymbiotic life.

Discover the Microbes Within! The Wolbachia Project: Citizen Science and Student-Based Discoveries for 15 Years and Counting.

The Elizabeth W. Jones Award for Excellence in Education recognizes an individual who has had a significant impact on genetics education at any education level. Seth R. Bordenstein, Ph.D., Centennial Professor of Biological Sciences at Vanderbilt University and Founding Director of the Vanderbilt Microbiome Initiative, is the 2020 recipient in recognition of his cofounding, developing, and expanding Discover the Microbes Within! The Wolbachia Project.

The phage gene wmk is a candidate for male killing by a bacterial endosymbiont.

Wolbachia are the most widespread maternally-transmitted bacteria in the animal kingdom. Their global spread in arthropods and varied impacts on animal physiology, evolution, and vector control are in part due to parasitic drive systems that enhance the fitness of infected females, the transmitting sex of Wolbachia. Male killing is one common drive mechanism wherein the sons of infected females are selectively killed. Despite decades of research, the gene(s) underlying Wolbachia-induced male killing remain unknown. Here using comparative genomic, transgenic, and cytological approaches in fruit flies, we identify a candidate gene in the eukaryotic association module of Wolbachia prophage WO, termed WO-mediated killing (wmk), which transgenically causes male-specific lethality during early embryogenesis and cytological defects typical of the pathology of male killing. The discovery of wmk establishes new hypotheses for the potential role of phage genes in sex-specific lethality, including the control of arthropod pests and vectors.

Author Correction: The Wolbachia mobilome in Culex pipiens includes a putative plasmid.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The Wolbachia mobilome in Culex pipiens includes a putative plasmid.

Wolbachia is a genus of obligate intracellular bacteria found in nematodes and arthropods worldwide, including insect vectors that transmit dengue, West Nile, and Zika viruses. Wolbachia's unique ability to alter host reproductive behavior through its temperate bacteriophage WO has enabled the development of new vector control strategies. However, our understanding of Wolbachia's mobilome beyond its bacteriophages is incomplete. Here, we reconstruct near-complete Wolbachia genomes from individual ovary metagenomes of four wild Culex pipiens mosquitoes captured in France. In addition to viral genes missing from the Wolbachia reference genome, we identify a putative plasmid (pWCP), consisting of a 9.23-kbp circular element with 14 genes. We validate its presence in additional Culex pipiens mosquitoes using PCR, long-read sequencing, and screening of existing metagenomes. The discovery of this previously unrecognized extrachromosomal element opens additional possibilities for genetic manipulation of Wolbachia.

Finer-Scale Phylosymbiosis: Insights from Insect Viromes.

Phylosymbiosis was recently proposed to describe the eco-evolutionary pattern whereby the ecological relatedness (e.g., beta diversity relationships) of host-associated microbial communities parallels the phylogeny of the host species. Representing the most abundant biological entities on the planet and common members of the animal-associated microbiome, viruses can be influential members of host-associated microbial communities that may recapitulate, reinforce, or ablate phylosymbiosis. Here we sequence the metagenomes of purified viral communities from three different parasitic wasp Nasonia species, one cytonuclear introgression line of Nasonia, and the flour moth outgroup Ephestia kuehniella. Results demonstrate complete phylosymbiosis between the viral metagenome and insect phylogeny. Across all Nasonia contigs, 69% of the genes in the viral metagenomes are either new to the databases or uncharacterized, yet over 99% of the contigs have at least one gene with similarity to a known sequence. The core Nasonia virome spans 21% of the total contigs, and the majority of that core is likely derived from induced prophages residing in the genomes of common Nasonia-associated bacterial genera: Proteus, Providencia, and Morganella. We also assemble the first complete viral particle genomes from Nasonia-associated gut bacteria. Taken together, results reveal the first complete evidence for phylosymbiosis in viral metagenomes, new genome sequences of viral particles from Nasonia-associated gut bacteria, and a large set of novel or uncharacterized genes in the Nasonia virome. This work suggests that phylosymbiosis at the host-microbiome level will likely extend to the host-virome level in other systems as well. IMPORTANCE Viruses are the most abundant biological entity on the planet and interact with microbial communities with which they associate. The virome of animals is often dominated by bacterial viruses, known as bacteriophages or phages, which can (re)structure bacterial communities potentially vital to the animal host. Beta diversity relationships of animal-associated bacterial communities in laboratory and wild populations frequently parallel animal phylogenetic relationships, a pattern termed phylosymbiosis. However, little is known about whether viral communities also exhibit this eco-evolutionary pattern. Metagenomics of purified viruses from recently diverged species of Nasonia parasitoid wasps reared in the lab indicates for the first time that the community relationships of the virome can also exhibit complete phylosymbiosis. Therefore, viruses, particularly bacteriophages here, may also be influenced by animal evolutionary changes either directly or indirectly through the tripartite interactions among hosts, bacteria, and phage communities. Moreover, we report several new bacteriophage genomes from the common gut bacteria in Nasonia.

Evolutionary Genetics of Cytoplasmic Incompatibility Genes cifA and cifB in Prophage WO of Wolbachia.

The bacterial endosymbiont Wolbachia manipulates arthropod reproduction to facilitate its maternal spread through host populations. The most common manipulation is cytoplasmic incompatibility (CI): Wolbachia-infected males produce modified sperm that cause embryonic mortality, unless rescued by embryos harboring the same Wolbachia. The genes underlying CI, cifA and cifB, were recently identified in the eukaryotic association module of Wolbachia's prophage WO. Here, we use transcriptomic and genomic approaches to address three important evolutionary facets of the cif genes. First, we assess whether or not cifA and cifB comprise a classic toxin-antitoxin operon in wMel and show that the two genes exhibit striking, transcriptional differences across host development. They can produce a bicistronic message despite a predicted hairpin termination element in their intergenic region. Second, cifA and cifB strongly coevolve across the diversity of phage WO. Third, we provide new domain and functional predictions across homologs within Wolbachia, and show that amino acid sequences vary substantially across the genus. Finally, we investigate conservation of cifA and cifB and find frequent degradation and loss of the genes in strains that no longer induce CI. Taken together, we demonstrate that cifA and cifB exhibit complex transcriptional regulation in wMel, provide functional annotations that broaden the potential mechanisms of CI induction, and report recurrent erosion of cifA and cifB in non-CI strains, thus expanding our understanding of the most widespread form of reproductive parasitism.

Prophage WO genes recapitulate and enhance Wolbachia-induced cytoplasmic incompatibility.

The genus Wolbachia is an archetype of maternally inherited intracellular bacteria that infect the germline of numerous invertebrate species worldwide. They can selfishly alter arthropod sex ratios and reproductive strategies to increase the proportion of the infected matriline in the population. The most common reproductive manipulation is cytoplasmic incompatibility, which results in embryonic lethality in crosses between infected males and uninfected females. Females infected with the same Wolbachia strain rescue this lethality. Despite more than 40 years of research and relevance to symbiont-induced speciation, as well as control of arbovirus vectors and agricultural pests, the bacterial genes underlying cytoplasmic incompatibility remain unknown. Here we use comparative and transgenic approaches to demonstrate that two differentially transcribed, co-diverging genes in the eukaryotic association module of prophage WO from Wolbachia strain wMel recapitulate and enhance cytoplasmic incompatibility. Dual expression in transgenic, uninfected males of Drosophila melanogaster crossed to uninfected females causes embryonic lethality. Each gene additively augments embryonic lethality in crosses between infected males and uninfected females. Lethality associates with embryonic defects that parallel those of wild-type cytoplasmic incompatibility and is notably rescued by wMel-infected embryos in all cases. The discovery of cytoplasmic incompatibility factor genes cifA and cifB pioneers genetic studies of prophage WO-induced reproductive manipulations and informs the continuing use of Wolbachia to control dengue and Zika virus transmission to humans.

Eukaryotic association module in phage WO genomes from Wolbachia.

Viruses are trifurcated into eukaryotic, archaeal and bacterial categories. This domain-specific ecology underscores why eukaryotic viruses typically co-opt eukaryotic genes and bacteriophages commonly harbour bacterial genes. However, the presence of bacteriophages in obligate intracellular bacteria of eukaryotes may promote DNA transfers between eukaryotes and bacteriophages. Here we report a metagenomic analysis of purified bacteriophage WO particles of Wolbachia and uncover a eukaryotic association module in the complete WO genome. It harbours predicted domains, such as the black widow latrotoxin C-terminal domain, that are uninterrupted in bacteriophage genomes, enriched with eukaryotic protease cleavage sites and combined with additional domains to forge one of the largest bacteriophage genes to date (14,256 bp). To the best of our knowledge, these eukaryotic-like domains have never before been reported in packaged bacteriophages and their phylogeny, distribution and sequence diversity imply lateral transfers between bacteriophage/prophage and animal genomes. Finally, the WO genome sequences and identification of attachment sites will potentially advance genetic manipulation of Wolbachia.

Recent genome reduction of Wolbachia in Drosophila recens targets phage WO and narrows candidates for reproductive parasitism.

Wolbachia are maternally transmitted endosymbionts that often alter their arthropod hosts' biology to favor the success of infected females, and they may also serve as a speciation microbe driving reproductive isolation. Two of these host manipulations include killing males outright and reducing offspring survival when infected males mate with uninfected females, a phenomenon known as cytoplasmic incompatibility. Little is known about the mechanisms behind these phenotypes, but interestingly either effect can be caused by the same Wolbachia strain when infecting different hosts. For instance, wRec causes cytoplasmic incompatibility in its native host Drosophila recens and male killing in D. subquinaria. The discovery of prophage WO elements in most arthropod Wolbachia has generated the hypothesis that WO may encode genes involved in these reproductive manipulations. However, PCR screens for the WO minor capsid gene indicated that wRec lacks phage WO. Thus, wRec seemed to provide an example where phage WO is not needed for Wolbachia-induced reproductive manipulation. To enable investigation of the mechanism of phenotype switching in different host backgrounds, and to examine the unexpected absence of phage WO, we sequenced the genome of wRec. Analyses reveal that wRec diverged from wMel approximately 350,000 years ago, mainly by genome reduction in the phage regions. While it lost the minor capsid gene used in standard PCR screens for phage WO, it retained two regions encompassing 33 genes, several of which have previously been associated with reproductive parasitism. Thus, WO gene involvement in reproductive manipulation cannot be excluded and reliance on single gene PCR should not be used to rule out the presence of phage WO in Wolbachia. Additionally, the genome sequence for wRec will enable transcriptomic and proteomic studies that may help elucidate the Wolbachia mechanisms of altered reproductive manipulations associated with host switching, perhaps among the 33 remaining phage genes.

Temperature affects the tripartite interactions between bacteriophage WO, Wolbachia, and cytoplasmic incompatibility.

Wolbachia infections are a model for understanding intracellular, bacterial symbioses. While the symbiosis is often studied from a binary perspective of host and bacteria, it is increasingly apparent that additional trophic levels can influence the symbiosis. For example, Wolbachia in arthropods harbor a widespread temperate bacteriophage, termed WO, that forms virions and rampantly transfers between coinfections. Here we test the hypothesis that temperatures at the extreme edges of an insect's habitable range alter bacteriophage WO inducibility and in turn, Wolbachia densities and the penetrance of cytoplasmic incompatibility. We report four key findings using the model wasp, Nasonia vitripennis: First, both cold treatment at 18 C and heat treatment at 30 C reduce Wolbachia densities by as much as 74% relative to wasps reared at 25 C. Second, in all cases where Wolbachia densities decline due to temperature changes, phage WO densities increase and inversely associate with Wolbachia densities. Heat has a marked effect on phage WO, yielding phage densities that are 552% higher than the room temperature control. Third, there is a significant affect of insect family on phage WO and endoysmbiont densities. Fourth, at extreme temperatures, there was a temperature-mediated adjustment to the density threshold at which Wolbachia cause complete cytoplasmic incompatibility. Taken together, these results demonstrate that temperature simultaneously affects phage WO densities, endosymbiont densities, and the penetrance of cytoplasmic incompatibility. While temperature shock enhances bacteriophage inducibility and the ensuing bacterial mortality in a wide range of medically and industrially-important bacteria, this is the first investigation of the associations in an obligate intracellular bacteria. Implications to a SOS global sensing feedback mechanism in Wolbachia are discussed.

Decoupling of host-symbiont-phage coadaptations following transfer between insect species.

Transferring endosymbiotic bacteria between different host species can perturb the coordinated regulation of the host and bacterial genomes. Here we use the most common maternally transmitted bacteria, Wolbachia pipientis, to test the consequences of host genetic background on infection densities and the processes underlying those changes in the parasitoid wasp genus Nasonia. Introgressing the genome of Nasonia giraulti into the infected cytoplasm of N. vitripennis causes a two-order-of-magnitude increase in bacterial loads in adults and a proliferation of the infection to somatic tissues. The host effect on W. pipientis distribution and densities is associated with a twofold decrease in densities of the temperate phage WO-B. Returning the bacteria from the new host species back to the resident host species restores the bacteria and phage to their native densities. To our knowledge, this is the first study to report a host-microbe genetic interaction that affects the densities of both W. pipientis and bacteriophage WO-B. The consequences of the increased bacterial density include a reduction in fecundity, an increase in levels of cytoplasmic incompatibility (CI), and unexpectedly, male-to-female transfer of the bacteria to uninfected females and an increased acceptance of densely infected females to interspecific mates. While paternal inheritance of the W. pipientis was not observed, the high incidence of male-to-female transfer in the introgressed background raises the possibility that paternal transmission could be more likely in hybrids where paternal leakage of other cytoplasmic elements is also known to occur. Taken together, these results establish a major change in W. pipientis densities and tissue tropism between closely related species and support a model in which phage WO, Wolbachia, and arthropods form a tripartite symbiotic association in which all three are integral to understanding the biology of this widespread endosymbiosis.