Extensive genomic diversity of closely related Wolbachia strains

Phage WO diversity and evolutionary forces associated with Wolbachia-infected crickets

Introduction

Phage WO represents the sole bacteriophage identified to infect Wolbachia, exerting a range of impacts on the ecological dynamics and evolutionary trajectories of its host. Given the extensive prevalence of Wolbachia across various species, phage WO is likely among the most prolific phage lineages within arthropod populations. To examine the diversity and evolutionary dynamics of phage WO, we conducted a screening for the presence of phage WO in Wolbachia-infected cricket species from China.

Methods

The presence of phage WO was detected using a PCR-based methodology. To elucidate the evolutionary forces driving phage WO diversity, analyses of intragenic recombination were conducted employing established recombination techniques, and horizontal transmission was investigated through comparative phylogenetic analysis of the phages and their hosts.

Results and discussion

Out of 19 cricket species infected with Wolbachia, 18 species were found to harbor phage WO. Notably, 13 of these 18 cricket species hosted multiple phage types, with the number of types ranging from two to 10, while the remaining five species harbored a single phage type. Twelve horizontal transmission events of phage WO were identified, wherein common phage WO types were shared among different Wolbachia strains. Notably, each phage WO horizontal transfer event was associated with distinct Wolbachia supergroups, specifically supergroups A, B, and F. Previous studies have found that four Wolbachia strains infect two to five species of crickets. However, among these cricket species, in addition to the shared phage WO types, all harbored species-specific phage WO types. This suggests that Wolbachia in crickets may acquire phage WO types through horizontal viral transfer between eukaryotes, independent of Wolbachia involvement. Furthermore, nine putative recombination events were identified across seven cricket species harboring multiple phage types. These findings suggest that horizontal transmission and intragenic recombination have played a significant role in the evolution of the phage WO genome, effectively enhancing the diversity of phage WO associated with crickets.

Prevalence of Wolbachia infection in field natural population of the apricot seed wasp Eurytoma samsonowi (Hymenoptera: Eurytomidae)

Obligate endosymbiont bacteria associated with insects are naturally providing their hosts with essential nutrients such as vitamins and amino acids and biological services including protection from pathogens. In this study, we aimed to investigate the presence of Wolbachia infection among males and females of the parasitic apricot seed wasp (ASW) Eurytoma samsonowi Vassiliev (Vassiliev Petrograd 11: 1-15, 1915) (Hymenoptera: Eurytomidae), a very harmful pest of apricot (Prunus armeniaca), in the oasis of Gafsa, Southern-West of Tunisia. The detection of Wolbachia infection was assessed based on the amplification of the Wolbachia surface protein (wsp) gene and a multilocus sequence typing (MLST) as a universal genotyping tool for Wolbachia involving the analyses of genes gatB, coxA, hcpA, fbpA, and ftsz. Confirming the screening results, Wolbachia was detected in the natural apricot wasp for the first time, with a significant difference between males (5%) and females (59%) based on wsp gene. All Wolbachia strains identified in E. samsonowi were clustered among supergroups B of Wolbachia.

Wolbachia Interactions with Diverse Insect Hosts: From Reproductive Modulations to Sustainable Pest Management Strategies

Effective in a variety of insect orders, including dipteran, lepidopteran, and hemipteran, Wolbachia-based control tactics are investigated, noting the importance of sterile and incompatible insect techniques. Encouraging approaches for controlling Aedes mosquitoes are necessary, as demonstrated by the evaluation of a new SIT/IIT combination and the incorporation of SIT into Drosophila suzukii management. For example, Wolbachia may protect plants from rice pests, demonstrating its potential for agricultural biological vector management. Maternal transmission and cytoplasmic incompatibility dynamics are explored, while Wolbachia phenotypic impacts on mosquito and rice pest management are examined. The importance of host evolutionary distance is emphasised in recent scale insect research that addresses host-shifting. Using greater information, a suggested method for comprehending Wolbachia host variations in various contexts emphasises ecological connectivity. Endosymbionts passed on maternally in nematodes and arthropods, Wolbachia are widely distributed around the world and have evolved both mutualistic and parasitic traits. Wolbachia is positioned as a paradigm for microbial symbiosis due to advancements in multiomics, gene functional assays, and its effect on human health. The challenges and opportunities facing Wolbachia research include scale issues, ecological implications, ethical conundrums, and the possibility of customising strains through genetic engineering. It is thought that cooperative efforts are required to include Wolbachia-based therapies into pest management techniques while ensuring responsible and sustainable ways.

Supergroup F Wolbachia with extremely reduced genome: transition to obligate insect symbionts

BACKGROUND

Wolbachia belong to highly abundant bacteria which are frequently found in invertebrate microbiomes and manifest by a broad spectrum of lifestyles from parasitism to mutualism. Wolbachia supergroup F is a particularly interesting clade as it gave rise to symbionts of both arthropods and nematodes, and some of its members are obligate mutualists. Investigations on evolutionary transitions among the different symbiotic stages have been hampered by a lack of the known diversity and genomic data for the supergroup F members.

RESULTS

Based on amplicon screening, short- and long-read WGS approaches, and laser confocal microscopy, we characterize five new supergroup F Wolbachia strains from four chewing lice species. These strains reached different evolutionary stages and represent two remarkably different types of symbiont genomes. Three of the genomes resemble other known members of Wolbachia F supergroup, while the other two show typical signs of ongoing gene inactivation and removal (genome size, coding density, low number of pseudogenes). Particularly, wMeur1, a symbiont fixed in microbiomes of Menacanthus eurysternus across four continents, possesses a highly reduced genome of 733,850 bp. The horizontally acquired capacity for pantothenate synthesis and localization in specialized bacteriocytes suggest its obligate nutritional role.

CONCLUSIONS

The genome of wMeur1 strain, from the M. eurysternus microbiome, represents the smallest currently known Wolbachia genome and the first example of Wolbachia which has completed genomic streamlining as known from the typical obligate symbionts. This points out that despite the large amount and great diversity of the known Wolbachia strains, evolutionary potential of these bacteria still remains underexplored. The diversity of the four chewing lice microbiomes indicates that this vast parasitic group may provide suitable models for further investigations. Video Abstract.

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.

Accumulation of endosymbiont genomes in an insect autosome followed by endosymbiont replacement

Eukaryotic genomes can acquire bacterial DNA via lateral gene transfer (LGT).¹ A prominent source of LGT is Wolbachia,² a widespread endosymbiont of arthropods and nematodes that is transmitted maternally through female germline cells.^3,4 The DNA transfer from the Wolbachia endosymbiont wAna to Drosophila ananassae is extensive^5-7 and has been localized to chromosome 4, contributing to chromosome expansion in this lineage.⁶ As has happened frequently with claims of bacteria-to-eukaryote LGT, the contribution of wAna transfers to the expanded size of D. ananassae chromosome 4 has been specifically contested⁸ owing to an assembly where Wolbachia sequences were classified as contaminants and removed.⁹ Here, long-read sequencing with DNA from a Wolbachia-cured line enabled assembly of 4.9 Mbp of nuclear Wolbachia transfers (nuwts) in D. ananassae and a 24-kbp nuclear mitochondrial transfer. The nuwts are <8,000 years old in at least two locations in chromosome 4 with at least one whole-genome integration followed by rapid extensive duplication of most of the genome with regions that have up to 10 copies. The genes in nuwts are accumulating small indels and mobile element insertions. Among the highly duplicated genes are cifA and cifB, two genes associated with Wolbachia-mediated Drosophila cytoplasmic incompatibility. The wAna strain that was the source of nuwts was subsequently replaced by a different wAna endosymbiont. Direct RNA Nanopore sequencing of Wolbachia-cured lines identified nuwt transcripts, including spliced transcripts, but functionality, if any, remains elusive.

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.

High Levels of Multiple Phage WO Infections and Its Evolutionary Dynamics Associated With Wolbachia-Infected Butterflies

Wolbachia is a maternally inherited bacterium that is widely distributed among arthropods, in which it manipulates the reproduction of its hosts. Phage WO is the only bacteriophage known to infect Wolbachia, and may provide benefit to its host or arthropods. We screened for the presence of phage WO in Wolbachia-infected butterfly species for the first time, to investigate their diversity and evolutionary dynamics. All Wolbachia-infected butterfly species, including members of the families Hesperiidae, Lycaenidae, Nymphalidae, Papilionidae, and Pieridae, were found to harbor phage WO. Interestingly, 84% of 19 butterfly species, which were infected with a single Wolbachia strain harbored high levels of multiple phage types (ranging from 3 to 17 types), another three species harbored one or two phage types. For Wolbachia strains (ST-41, ST-19, ST-125 and ST-374) shared among various butterfly species, their host insects all harbored multiple phage types, while two Wolbachia strains (ST-297 and ST-wPcau) were found to infect one butterfly species, whose insect hosts harbored a single phage type, suggesting that horizontal transfer of Wolbachia between insects increased the likelihood of exposure to phages, resulting in increased phage genetic diversity. Twelve horizontal transmission events of phage WO were found, which shared common phage WO types among different Wolbachia strains associated with butterflies. Most horizontal transfer events involved different Wolbachia supergroups (A and B). Horizontal acquisition of phage WO might also occur between eukaryotes without Wolbachia transfer. Furthermore, 22 putative recombination events were identified in 13 of 16 butterfly species which harbored multiple phage types. These results showed that horizontal transfer of Wolbachia caused it to be exposed to the phage gene pool, and that horizontal transmission of phage WO, as well as intragenic recombination were important dynamics for phage WO genome evolution, which effectively promoted the high level of phage WO diversity associated with butterflies.

Temperature effects on cellular host-microbe interactions explain continent-wide endosymbiont prevalence

Endosymbioses influence host physiology, reproduction, and fitness, but these relationships require efficient microbe transmission between host generations to persist. Maternally transmitted Wolbachia are the most common known endosymbionts,¹ but their frequencies vary widely within and among host populations for unknown reasons.^2,3 Here, we integrate genomic, cellular, and phenotypic analyses with mathematical models to provide an unexpectedly simple explanation for global wMel Wolbachia prevalence in Drosophila melanogaster. Cooling temperatures decrease wMel cellular abundance at a key stage of host oogenesis, producing temperature-dependent variation in maternal transmission that plausibly explains latitudinal clines of wMel frequencies on multiple continents. wMel sampled from a temperate climate targets the germline more efficiently in the cold than a recently differentiated tropical variant (∼2,200 years ago), indicative of rapid wMel adaptation to climate. Genomic analyses identify a very narrow list of wMel alleles-most notably, a derived stop codon in the major Wolbachia surface protein WspB-that underlie thermal sensitivity of cellular Wolbachia abundance and covary with temperature globally. Decoupling temperate wMel and host genomes further reduces transmission in the cold, a pattern that is characteristic of host-microbe co-adaptation to a temperate climate. Complex interactions among Wolbachia, hosts, and the environment (GxGxE) mediate wMel cellular abundance and maternal transmission, implicating temperature as a key determinant of Wolbachia spread and equilibrium frequencies, in conjunction with Wolbachia effects on host fitness and reproduction.^4,5 Our results motivate the strategic use of locally selected wMel variants for Wolbachia-based biocontrol efforts, which protect millions of individuals from arboviruses that cause human disease.⁶.

Positive Selection and Horizontal Gene Transfer in the Genome of a Male-Killing Wolbachia

Wolbachia are a genus of widespread bacterial endosymbionts in which some strains can hijack or manipulate arthropod host reproduction. Male killing is one such manipulation in which these maternally transmitted bacteria benefit surviving daughters in part by removing competition with the sons for scarce resources. Despite previous findings of interesting genome features of microbial sex ratio distorters, the population genomics of male-killers remain largely uncharacterized. Here, we uncover several unique features of the genome and population genomics of four Arizonan populations of a male-killing Wolbachia strain, wInn, that infects mushroom-feeding Drosophila innubila. We first compared the wInn genome with other closely related Wolbachia genomes of Drosophila hosts in terms of genome content and confirm that the wInn genome is largely similar in overall gene content to the wMel strain infecting D. melanogaster. However, it also contains many unique genes and repetitive genetic elements that indicate lateral gene transfers between wInn and non-Drosophila eukaryotes. We also find that, in line with literature precedent, genes in the Wolbachia prophage and Octomom regions are under positive selection. Of all the genes under positive selection, many also show evidence of recent horizontal transfer among Wolbachia symbiont genomes. These dynamics of selection and horizontal gene transfer across the genomes of several Wolbachia strains and diverse host species may be important underlying factors in Wolbachia's success as a male-killer of divergent host species.

A Case of Intragenic Recombination Dramatically Impacting the Phage WO Genetic Diversity in Gall Wasps

The phage WO was characterized in Wolbachia, a strictly intracellular bacterium causing several reproductive alterations in its arthropod hosts. This study aimed to screen the presence of Wolbachia and phage WO in 15 gall wasp species from six provinces of southern China to investigate their diversity and prevalence patterns. A high incidence of Wolbachia infection was determined in the gall wasp species, with an infection rate of 86.7% (13/15). Moreover, seven species had double or multiple infections. All Wolbachia-infected gall wasp species were found to harbor phage WO. The gall wasp species infected with a single Wolbachia strain were found to harbor a single phage WO type. On the contrary, almost all species with double or multiple Wolbachia infections harbored a high level of phage WO diversity (ranging from three to 27 types). Six horizontal transfer events of phage WO in Wolbachia were found to be associated with gall wasps, which shared identical orf7 sequences among their respective accomplices. The transfer potentially took place through gall inducers and associated inquilines infected with or without Wolbachia. Furthermore, 10 putative recombination events were identified from Andricus hakonensis and Andricus sp2, which harbored multiple phage WO types, suggesting that intragenic recombination was the important evolutionary force, which effectively promoted the high level of phage WO diversity associated with gall wasps.

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.

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.

Comparative Genomics Reveals Factors Associated with Phenotypic Expression of Wolbachia

Wolbachia is a widespread, vertically transmitted bacterial endosymbiont known for manipulating arthropod reproduction. Its most common form of reproductive manipulation is cytoplasmic incompatibility (CI), observed when a modification in the male sperm leads to embryonic lethality unless a compatible rescue factor is present in the female egg. CI attracts scientific attention due to its implications for host speciation and in the use of Wolbachia for controlling vector-borne diseases. However, our understanding of CI is complicated by the complexity of the phenotype, whose expression depends on both symbiont and host factors. In the present study, we perform a comparative analysis of nine complete Wolbachia genomes with known CI properties in the same genetic host background, Drosophila simulans STC. We describe genetic differences between closely related strains and uncover evidence that phages and other mobile elements contribute to the rapid evolution of both genomes and phenotypes of Wolbachia. Additionally, we identify both known and novel genes associated with the modification and rescue functions of CI. We combine our observations with published phenotypic information and discuss how variability in cif genes, novel CI-associated genes, and Wolbachia titer might contribute to poorly understood aspects of CI such as strength and bidirectional incompatibility. We speculate that high titer CI strains could be better at invading new hosts already infected with a CI Wolbachia, due to a higher rescue potential, and suggest that titer might thus be a relevant parameter to consider for future strategies using CI Wolbachia in biological control.

Evolution of Wolbachia mutualism and reproductive parasitism: insight from two novel strains that co-infect cat fleas

Wolbachiae are obligate intracellular bacteria that infect arthropods and certain nematodes. Usually maternally inherited, they may provision nutrients to (mutualism) or alter sexual biology of (reproductive parasitism) their invertebrate hosts. We report the assembly of closed genomes for two novel wolbachiae, wCfeT and wCfeJ, found co-infecting cat fleas (Ctenocephalides felis) of the Elward Laboratory colony (Soquel, CA, USA). wCfeT is basal to nearly all described Wolbachia supergroups, while wCfeJ is related to supergroups C, D and F. Both genomes contain laterally transferred genes that inform on the evolution of Wolbachia host associations. wCfeT carries the Biotin synthesis Operon of Obligate intracellular Microbes (BOOM); our analyses reveal five independent acquisitions of BOOM across the Wolbachia tree, indicating parallel evolution towards mutualism. Alternately, wCfeJ harbors a toxin-antidote operon analogous to the wPip cinAB operon recently characterized as an inducer of cytoplasmic incompatibility (CI) in flies. wCfeJ cinB and three adjacent genes are collectively similar to large modular toxins encoded in CI-like operons of certain Wolbachia strains and Rickettsia species, signifying that CI toxins streamline by fission of large modular toxins. Remarkably, the C. felis genome itself contains two CI-like antidote genes, divergent from wCfeJ cinA, revealing episodic reproductive parasitism in cat fleas and evidencing mobility of CI loci independent of WO-phage. Additional screening revealed predominant co-infection (wCfeT/wCfeJ) amongst C. felis colonies, though fleas in wild populations mostly harbor wCfeT alone. Collectively, genomes of wCfeT, wCfeJ, and their cat flea host supply instances of lateral gene transfers that could drive transitions between parasitism and mutualism.

Endosymbiotic Rickettsiella causes cytoplasmic incompatibility in a spider host

Many arthropod hosts are infected with bacterial endosymbionts that manipulate host reproduction, but few bacterial taxa have been shown to cause such manipulations. Here, we show that a bacterial strain in the genus Rickettsiella causes cytoplasmic incompatibility (CI) between infected and uninfected hosts. We first surveyed the bacterial community of the agricultural spider Mermessus fradeorum (Linyphiidae) using high throughput sequencing and found that individual spiders can be infected with up to five different strains of maternally inherited symbiont from the genera Wolbachia, Rickettsia, and Rickettsiella. The Rickettsiella strain was pervasive, found in all 23 tested spider matrilines. We used antibiotic curing to generate uninfected matrilines that we reciprocally crossed with individuals infected only with Rickettsiella. We found that only 13% of eggs hatched when uninfected females were mated with Rickettsiella-infected males; in contrast, at least 83% of eggs hatched in the other cross types. This is the first documentation of Rickettsiella, or any Gammaproteobacteria, causing CI. We speculate that induction of CI may be much more widespread among maternally inherited bacteria than previously appreciated. Further, our results reinforce the importance of thoroughly characterizing and assessing the inherited microbiome before attributing observed host phenotypes to well-characterized symbionts such as Wolbachia.

Endosymbiont diversity in natural populations of Tetranychus mites is rapidly lost under laboratory conditions

Although the diversity of bacterial endosymbionts in arthropods is well documented, whether and how such diversity is maintained remains an open question. We investigated the temporal changes occurring in the prevalence and composition of endosymbionts after transferring natural populations of Tetranychus spider mites from the field to the laboratory. These populations, belonging to three different Tetranychus species (T. urticae, T. ludeni and T. evansi) carried variable infection frequencies of Wolbachia, Cardinium, and Rickettsia. We report a rapid change of the infection status of these populations after only 6 months of laboratory rearing, with an apparent loss of Rickettsia and Cardinium, while Wolbachia apparently either reached fixation or was lost. We show that Wolbachia had variable effects on host longevity and fecundity, and induced variable levels of cytoplasmic incompatibility (CI) in each fully infected population, despite no sequence divergence in the markers used and full CI rescue between all populations. This suggests that such effects are largely dependent upon the host genotype. Subsequently, we used these data to parameterize a theoretical model for the invasion of CI-inducing symbionts in haplodiploids, which shows that symbiont effects are sufficient to explain their dynamics in the laboratory. This further suggests that symbiont diversity and prevalence in the field are likely maintained by environmental heterogeneity, which is reduced in the laboratory. Overall, this study highlights the lability of endosymbiont infections and draws attention to the limitations of laboratory studies to understand host-symbiont interactions in natural populations.

Wolbachia Endosymbiont of the Horn Fly (Haematobia irritans irritans): a Supergroup A Strain with Multiple Horizontally Acquired Cytoplasmic Incompatibility Genes

The horn fly, Haematobia irritans irritans, is a hematophagous parasite of livestock distributed throughout Europe, Africa, Asia, and the Americas. Welfare losses on livestock due to horn fly infestation are estimated to cost between $1 billion and $2.5 billion (U.S. dollars) annually in North America and Brazil. The endosymbiotic bacterium Wolbachia pipientis is a maternally inherited manipulator of reproductive biology in arthropods and naturally infects laboratory colonies of horn flies from Kerrville, TX, and Alberta, Canada, but it has also been identified in wild-caught samples from Canada, the United States, Mexico, and Hungary. Reassembly of PacBio long-read and Illumina genomic DNA libraries from the Kerrville H. i. irritans genome project allowed for a complete and circularized 1.3-Mb Wolbachia genome (wIrr). Annotation of wIrr yielded 1,249 coding genes, 34 tRNAs, 3 rRNAs, and 5 prophage regions. Comparative genomics and whole-genome Bayesian evolutionary analysis of wIrr compared to published Wolbachia genomes suggested that wIrr is most closely related to and diverged from Wolbachia supergroup A strains known to infect Drosophila spp. Whole-genome synteny analyses between wIrr and closely related genomes indicated that wIrr has undergone significant genome rearrangements while maintaining high nucleotide identity. Comparative analysis of the cytoplasmic incompatibility (CI) genes of wIrr suggested two phylogenetically distinct CI loci and acquisition of another cifB homolog from phylogenetically distant supergroup A Wolbachia strains, suggesting horizontal acquisition of these loci. The wIrr genome provides a resource for future examination of the impact Wolbachia may have in both biocontrol and potential insecticide resistance of horn flies.IMPORTANCE Horn flies, Haematobia irritans irritans, are obligate hematophagous parasites of cattle having significant effects on production and animal welfare. Control of horn flies mainly relies on the use of insecticides, but issues with resistance have increased interest in development of alternative means of control. Wolbachia pipientis is an endosymbiont bacterium known to have a range of effects on host reproduction, such as induction of cytoplasmic incompatibility, feminization, male killing, and also impacts vector transmission. These characteristics of Wolbachia have been exploited in biological control approaches for a range of insect pests. Here we report the assembly and annotation of the circular genome of the Wolbachia strain of the Kerrville, TX, horn fly (wIrr). Annotation of wIrr suggests its unique features, including the horizontal acquisition of additional transcriptionally active cytoplasmic incompatibility loci. This study provides the foundation for future studies of Wolbachia-induced biological effects for control of horn flies.

Transgenic Testing Does Not Support a Role for Additional Candidate Genes in Wolbachia Male Killing or Cytoplasmic Incompatibility

Wolbachia are widespread bacterial endosymbionts that manipulate the reproduction of diverse arthropods to spread through a population and can substantially shape host evolution. Recently, reports identified three prophage WO genes (wmk, cifA, and cifB) that transgenically recapitulate many aspects of reproductive manipulation in Drosophila melanogaster. Here, we transgenically tested 10 additional gene candidates for CI and/or male killing in flies. The results yield no evidence for the involvement of these gene candidates in reproductive parasitism, bolstering the evidence for identification of the cif and wmk genes as the major factors involved in their phenotypes. In addition, evidence supports new hypotheses for prediction of male-killing phenotypes or lack thereof based on wmk transcript length and copy number. These experiments inform efforts to understand the full basis of reproductive parasitism for basic and applied purposes and lay the foundation for future work on the function of an interesting group of Wolbachia and phage WO genes.

Endosymbiotic bacteria in the genus Wolbachia remarkably infect nearly half of all arthropod species. They spread in part because of manipulations of host sexual reproduction that enhance the maternal transmission of the bacteria, including male killing (death of infected males) and unidirectional cytoplasmic incompatibility (CI; death of offspring from infected fathers and uninfected mothers). Recent discoveries identified several genes in prophage WO of Wolbachia (wmk, cifA, and cifB) that fully or partially recapitulate male killing or CI when transgenically expressed in Drosophila melanogaster. However, it is not yet fully resolved if other gene candidates contribute to these phenotypes. Here, we transgenically tested 10 additional gene candidates for their involvement in male killing and/or CI. The results show that despite sequence and protein architecture similarities or comparative associations with reproductive parasitism, transgenic expression of the candidates does not recapitulate male killing or CI. Sequence analysis across Wmk and its closest relatives reveals amino acids that may be important to its function. In addition, evidence is presented to propose new hypotheses regarding the relationship between wmk transcript length and its ability to kill a given host, as well as copy number of wmk homologs within a bacterial strain, which may be predictive of host resistance. Together, these analyses continue to build the evidence for identification of wmk, cifA, and cifB as the major genes that have thus far been shown to cause reproductive parasitism in Wolbachia, and the transgenic resources provide a basis for further functional study of phage WO genes.

IMPORTANCEWolbachia are widespread bacterial endosymbionts that manipulate the reproduction of diverse arthropods to spread through a population and can substantially shape host evolution. Recently, reports identified three prophage WO genes (wmk, cifA, and cifB) that transgenically recapitulate many aspects of reproductive manipulation in Drosophila melanogaster. Here, we transgenically tested 10 additional gene candidates for CI and/or male killing in flies. The results yield no evidence for the involvement of these gene candidates in reproductive parasitism, bolstering the evidence for identification of the cif and wmk genes as the major factors involved in their phenotypes. In addition, evidence supports new hypotheses for prediction of male-killing phenotypes or lack thereof based on wmk transcript length and copy number. These experiments inform efforts to understand the full basis of reproductive parasitism for basic and applied purposes and lay the foundation for future work on the function of an interesting group of Wolbachia and phage WO genes.

Wolbachia Horizontal Transmission Events in Ants: What Do We Know and What Can We Learn?

While strict vertical transmission insures the durability of intracellular symbioses, phylogenetic incongruences between hosts and endosymbionts suggest horizontal transmission must also occur. These horizontal acquisitions can have important implications for the biology of the host. Wolbachia is one of the most ecologically successful prokaryotes in arthropods, infecting an estimated 50–70% of all insect species. Much of this success is likely due to the fact that, in arthropods, Wolbachia is notorious for manipulating host reproduction to favor transmission through the female germline. However, its natural potential for horizontal transmission remains poorly understood. Here we evaluate the fundamental prerequisites for successful horizontal transfer, including necessary environmental conditions, genetic potential of bacterial strains, and means of mediating transfers. Furthermore, we revisit the relatedness of Wolbachia strains infecting the Panamanian leaf-cutting ant, Acromyrmex echinatior, and its inquiline social parasite, Acromyrmex insinuator, and compare our results to a study published more than 15 years ago by Van Borm et al. (2003). The results of this pilot study prompt us to reevaluate previous notions that obligate social parasitism reliably facilitates horizontal transfer and suggest that not all Wolbachia strains associated with ants have the same genetic potential for horizontal transmission.

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.

Contrasting Patterns of Virus Protection and Functional Incompatibility Genes in Two Conspecific Wolbachia Strains from Drosophila pandora

Wolbachia infections can present different phenotypes in hosts, including different forms of reproductive manipulation and antiviral protection, which may influence infection dynamics within host populations. In populations of Drosophila pandora two distinct Wolbachia strains coexist, each manipulating host reproduction: strain wPanCI causes cytoplasmic incompatibility (CI), whereas strain wPanMK causes male killing (MK). CI occurs when a Wolbachia-infected male mates with a female not infected with a compatible type of Wolbachia, leading to nonviable offspring. wPanMK can rescue wPanCI-induced CI but is unable to induce CI. The antiviral protection phenotypes provided by the wPanCI and wPanMK infections were characterized; the strains showed differential protection phenotypes, whereby cricket paralysis virus (CrPV)-induced mortality was delayed in flies infected with wPanMK but enhanced in flies infected with wPanCI compared to their respective Wolbachia-cured counterparts. Homologs of the cifA and cifB genes involved in CI identified in wPanMK and wPanCI showed a high degree of conservation; however, the CifB protein in wPanMK is truncated and is likely nonfunctional. The presence of a likely functional CifA in wPanMK and wPanMK's ability to rescue wPanCI-induced CI are consistent with the recent confirmation of CifA's involvement in CI rescue, and the absence of a functional CifB protein further supports its involvement as a CI modification factor. Taken together, these findings indicate that wPanCI and wPanMK have different relationships with their hosts in terms of their protective and CI phenotypes. It is therefore likely that different factors influence the prevalence and dynamics of these coinfections in natural Drosophila pandora hosts.IMPORTANCEWolbachia strains are common endosymbionts in insects, with multiple strains often coexisting in the same species. The coexistence of multiple strains is poorly understood but may rely on Wolbachia organisms having diverse phenotypic effects on their hosts. As Wolbachia is increasingly being developed as a tool to control disease transmission and suppress pest populations, it is important to understand the ways in which multiple Wolbachia strains persist in natural populations and how these might then be manipulated. We have therefore investigated viral protection and the molecular basis of cytoplasmic incompatibility in two coexisting Wolbachia strains with contrasting effects on host reproduction.

Genome organisation and comparative genomics of four novel Wolbachia genome assemblies from Indian Drosophila host

Wolbachia has long been known to share an endosymbiotic relationship with its host as an obligate intracellular organism. Wolbachia diversity as different supergroups is found to be host-specific in most cases except a few, where the host species is seen to accommodate multiple strains. Besides, the Wolbachia genome must have undergone several changes in response to the evolving host genome in order to adapt and establish a strong association with its host, thus making a distinctive Wolbachia-host alliance. The present study focusses on four novel genome assembly and genome-wide sequence variations of Indian Wolbachia strains, i.e. wMel and wRi isolated from two different Drosophila hosts. The genome assembly has an average size of ~ 1.1 Mb and contains ~ 1100 genes, which is comparable with the previously sequenced Wolbachia genomes. The comparative genomics analysis of these genomes and sequence-wide comparison of some functionally significant genes, i.e. ankyrin repeats, Wsp and T4SS, highlight their sequence similarities and dissimilarities, further supporting the strain-specific association of Wolbachia to its host. Interestingly, some of the sequence variations are also found to be restricted to only Indian Wolbachia strains. Further analysis of prophage and their flanking regions in the Wolbachia genome reveals the presence of several functional genes which may assist the phage to reside inside the bacterial host, thus providing a trade-off for the endosymbiont-host association. Understanding this endosymbiont genome in different eco-geographical conditions has become imperative for the recent use of Wolbachia in medical entomology as a vector-control agent.

Biological Control of Mosquito-Borne Diseases: The Potential of Wolbachia-Based Interventions in an IVM Framework

People living in the tropical and subtropical regions of the world face an enormous health burden due to mosquito-borne diseases such as malaria, dengue fever, and filariasis. Historically and today, targeting mosquito vectors with, primarily, insecticide-based control strategies have been a key control strategy against major mosquito-borne diseases. However, the success to date of such approaches is under threat from multiple insecticide resistance mechanisms while vector control (VC) options are still limited. The situation therefore requires the development of innovative control measures against major mosquito-borne diseases. Transinfecting mosquitos with symbiotic bacteria that can compete with targeted pathogens or manipulate host biology to reduce their vectorial capacity are a promising and innovative biological control approach. In this review, we discuss the current state of knowledge about the association between mosquitoes and Wolbachia, emphasizing the limitations of different mosquito control strategies and the use of mosquitoes' commensal microbiota as innovative approaches to control mosquito-borne diseases.

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.

Family level variation in Wolbachia-mediated dengue virus blocking in Aedes aegypti

BACKGROUND

The mosquito vector Aedes aegypti is responsible for transmitting a range of arboviruses including dengue (DENV) and Zika (ZIKV). The global reach of these viruses is increasing due to an expansion of the mosquito's geographic range and increasing urbanization and human travel. Vector control remains the primary means for limiting these diseases. Wolbachia pipientis is an endosymbiotic bacterium of insects that has the ability to block the replication of pathogens, including flaviviruses such as DENV or ZIKV, inside the body of the vector. A strain of Wolbachia called wMel is currently being released into wild mosquito populations to test its potential to limit virus transmission to humans. The mechanism that underpins the virus blocking effect, however, remains elusive.

METHODS

We used a modified full-sib breeding design in conjunction with vector competence assays in wildtype and wMel-infected Aedes aegypti collected from the field. All individuals were injected with DENV-2 intrathoracically at 5-6 days of age. Tissues were dissected 7 days post-infection to allow quantification of DENV and Wolbachia loads.

RESULTS

We show the first evidence of family level variation in Wolbachia-mediated blocking in mosquitoes. This variation may stem from either genetic contributions from the mosquito and Wolbachia genomes or environmental influences on Wolbachia. In these families, we also tested for correlations between strength of blocking and expression level for several insect immunity genes with possible roles in blocking, identifying two genes of interest (AGO2 and SCP-2).

CONCLUSIONS

In this study we show variation in Wolbachia-mediated DENV blocking in Aedes aegypti that may arise from genetic contributions and environmental influences on the mosquito-Wolbachia association. This suggests that Wolbachia-mediated blocking may have the ability to evolve through time or be expressed differentially across environments. The long-term efficacy of Wolbachia in the field will be dependent on the stability of blocking. Understanding the mechanism of blocking will be necessary for successful development of strategies that counter the emergence of evolved resistance or variation in its expression under diverse field conditions.

Large-Scale Identification of Wolbachia pipientis Effectors

Wolbachia pipientis is an intracellular symbiont of arthropods well known for the reproductive manipulations induced in the host and, more recently, for the ability of Wolbachia to block virus replication in insect vectors. Since Wolbachia cannot yet be genetically manipulated, and due to the constraints imposed when working with an intracellular symbiont, little is known about mechanisms used by Wolbachia for host interaction. Here we employed a bioinformatics pipeline and identified 163 candidate effectors, potentially secreted by Wolbachia into the host cell. A total of 84 of these candidates were then subjected to a screen of growth defects induced in yeast upon heterologous expression which identified 14 top candidates likely secreted by Wolbachia. These predicted secreted effectors may function in concert as we find that their native expression is correlated and is highly upregulated at specific time points during Drosophila development. In addition, the evolutionary histories of some of these predicted effectors are also correlated, suggesting they may function together, or in the same pathway, during host infection. Similarly, most of these predicted effectors are limited to one or two Wolbachia strains—perhaps reflecting shared evolutionary history and strain specific functions in host manipulation. Identification of these Wolbachia candidate effectors is the first step in dissecting the mechanisms of symbiont–host interaction in this important system.

Insertion sequence polymorphism and genomic rearrangements uncover hidden Wolbachia diversity in Drosophila suzukii and D. subpulchrella

Ability to distinguish between closely related Wolbachia strains is crucial for understanding the evolution of Wolbachia-host interactions and the diversity of Wolbachia-induced phenotypes. A useful model to tackle these issues is the Drosophila suzukii - Wolbachia association. D. suzukii, a destructive insect pest, harbor a non-CI inducing Wolbachia 'wSuz' closely related to the strong CI-inducing wRi strain. Multi locus sequence typing (MLST) suggests presence of genetic homogeneity across wSuz strains infecting European and American D. suzukii populations, although different Wolbachia infection frequencies and host fecundity levels have been observed in both populations. Currently, it is not clear if these differences are due to cryptic wSuz polymorphism, host background, geographical factors or a combination of all of them. Here, we have identified geographical diversity in wSuz in D. suzukii populations from different continents using a highly diagnostic set of markers based on insertion sequence (IS) site polymorphism and genomic rearrangements (GR). We further identified inter-strain diversity between Wolbachia infecting D. suzukii and its sister species D. subpulchrella (wSpc). Based on our results, we speculate that discernible wSuz variants may associate with different observed host phenotypes, a hypothesis that demands future investigation. More generally, our results demonstrate the utility of IS and GRs in discriminating closely related Wolbachia strains.

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 (CI), 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^2,3, as well as control of arbovirus vectors^4,5,6 and agricultural pests⁷, the bacterial genes underlying CI remain unknown. Here, we use comparative and transgenic approaches to demonstrate that two differentially transcribed, codiverging genes in the eukaryotic association module of prophage WO⁸ from Wolbachia strain wMel recapitulate and enhance CI. Dual expression in transgenic, uninfected males of Drosophila melanogaster crossed to uninfected females causes embryonic lethality. Each gene additively augments embryonic lethality in infected males crossed to uninfected females. Lethality associates with embryonic defects that parallel those of wild type CI 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 Wolbachia’s ongoing utility to control dengue and Zika transmission to humans.

Molecular diversity of Wolbachia in Lepidoptera: Prevalent allelic content and high recombination of MLST genes

Wolbachia are common endosymbiotic bacteria of Arthropoda and Nematoda that are ordinarily transmitted vertically in host lineages through the egg cytoplasm. Despite the great interest in the Wolbachia symbiont, many issues of its biology remain unclear, including its evolutionary history, routes of transfer among species, and the molecular mechanisms underlying the symbiont's effect on its host. In this report, we present data relating to Wolbachia infection in 120 species of 13 Lepidoptera families, mostly butterflies, from West Siberian localities based on Multilocus sequence typing (MLST) and the wsp locus and perform a comprehensive survey of the distribution of Wolbachia and its genetic diversity in Lepidoptera worldwide. We observed a high infection incidence in the studied region; this finding is probably also true for other temperate latitude regions because many studied species have broad Palearctic and even Holarctic distribution. Although 40 new MLST alleles and 31 new STs were described, there was no noticeable difference in the MLST allele content in butterflies and probably also in moths worldwide. A genetic analysis of Wolbachia strains revealed the MLST allele core in lepidopteran hosts worldwide, viz. the ST-41 allele content. The key finding of our study was the detection of rampant recombination among MLST haplotypes. High rates of homologous recombination between Wolbachia strains indicate a substantial contribution of genetic exchanges to the generation of new STs. This finding should be considered when discussing issues related to the reconstruction of Wolbachia evolution, divergence time, and the routes of Wolbachia transmission across arthropod hosts.

Bacteriophage WO Can Mediate Horizontal Gene Transfer in Endosymbiotic Wolbachia Genomes

Phage-mediated horizontal gene transfer (HGT) is common in free-living bacteria, and many transferred genes can play a significant role in their new bacterial hosts. However, there are few reports concerning phage-mediated HGT in endosymbionts (obligate intracellular bacteria within animal or plant hosts), such as Wolbachia. The Wolbachia-infecting temperate phage WO can actively shift among Wolbachia genomes and has the potential to mediate HGT between Wolbachia strains. In the present study, we extend previous findings by validating that the phage WO can mediate transfer of non-phage genes. To do so, we utilized bioinformatic, phylogenetic, and molecular analyses based on all sequenced Wolbachia and phage WO genomes. Our results show that the phage WO can mediate HGT between Wolbachia strains, regardless of whether the transferred genes originate from Wolbachia or other unrelated bacteria.

Disentangling a Holobiont - Recent Advances and Perspectives in Nasonia Wasps

The parasitoid wasp genus Nasonia (Hymenoptera: Chalcidoidea) is a well-established model organism for insect development, evolutionary genetics, speciation, and symbiosis. The host-microbiota assemblage which constitutes the Nasonia holobiont (a host together with all of its associated microbes) consists of viruses, two heritable bacterial symbionts and a bacterial community dominated in abundance by a few taxa in the gut. In the wild, all four Nasonia species are systematically infected with the obligate intracellular bacterium Wolbachia and can additionally be co-infected with Arsenophonus nasoniae. These two reproductive parasites have different transmission modes and host manipulations (cytoplasmic incompatibility vs. male-killing, respectively). Pioneering studies on Wolbachia in Nasonia demonstrated that closely related Nasonia species harbor multiple and mutually incompatible Wolbachia strains, resulting in strong symbiont-mediated reproductive barriers that evolved early in the speciation process. Moreover, research on host-symbiont interactions and speciation has recently broadened from its historical focus on heritable symbionts to the entire microbial community. In this context, each Nasonia species hosts a distinguishable community of gut bacteria that experiences a temporal succession during host development and members of this bacterial community cause strong hybrid lethality during larval development. In this review, we present the Nasonia species complex as a model system to experimentally investigate questions regarding: (i) the impact of different microbes, including (but not limited to) heritable endosymbionts, on the extended phenotype of the holobiont, (ii) the establishment and regulation of a species-specific microbiota, (iii) the role of the microbiota in speciation, and (iv) the resilience and adaptability of the microbiota in wild populations subjected to different environmental pressures. We discuss the potential for easy microbiota manipulations in Nasonia as a promising experimental approach to address these fundamental aspects.

Wolbachia from Drosophila incompta: just a hitchhiker shared by Drosophila in the New and Old World?

Wolbachia are intracellular endosymbionts that infect arthropods and filarial nematodes, occasionally causing a wide variety of modifications in host biology, such as male-killing and cytoplasmic incompatibility (CI), amongst others. This study assembled draft genomes for Wolbachia infecting Drosophila incompta, a species that uses flowers as exclusive breeding and feeding sites, in two distinct Brazilian populations. The absence of four genes involved in CI from this genome, together with literature reports of low frequencies of infected flies in wild populations that contain high mitogenome polymorphism, suggests that this bacterium does not induce CI in D. incompta. Phylogenomic analysis placed Wolbachia infecting D. incompta as closely related to the wMel strain which received such name since it was originally detected in Drosophila melanogaster. In addition, phylogenetic analysis using the Wolbachia surface protein gene and five genes used for multilocus sequence typing of Wolbachia found infecting Drosophila and other arthropod species of Old and New World displayed a complex evolutionary scenario involving recent horizontal transfer bursts in all major clades of Wolbachia pipens belonging to the supergroup A in both geographical regions.

Comparative Genomics of a Parthenogenesis-Inducing Wolbachia Symbiont

Wolbachia is an intracellular symbiont of invertebrates responsible for inducing a wide variety of phenotypes in its host. These host-Wolbachia relationships span the continuum from reproductive parasitism to obligate mutualism, and provide a unique system to study genomic changes associated with the evolution of symbiosis. We present the genome sequence from a parthenogenesis-inducing Wolbachia strain (wTpre) infecting the minute parasitoid wasp Trichogramma pretiosum. The wTpre genome is the most complete parthenogenesis-inducing Wolbachia genome available to date. We used comparative genomics across 16 Wolbachia strains, representing five supergroups, to identify a core Wolbachia genome of 496 sets of orthologous genes. Only 14 of these sets are unique to Wolbachia when compared to other bacteria from the Rickettsiales. We show that the B supergroup of Wolbachia, of which wTpre is a member, contains a significantly higher number of ankyrin repeat-containing genes than other supergroups. In the wTpre genome, there is evidence for truncation of the protein coding sequences in 20% of ORFs, mostly as a result of frameshift mutations. The wTpre strain represents a conversion from cytoplasmic incompatibility to a parthenogenesis-inducing lifestyle, and is required for reproduction in the Trichogramma host it infects. We hypothesize that the large number of coding frame truncations has accompanied the change in reproductive mode of the wTpre strain.

Identification and Characterization of a Candidate Wolbachia pipientis Type IV Effector That Interacts with the Actin Cytoskeleton

Many bacteria live as intracellular symbionts, causing persistent infections within insects. One extraordinarily common infection is that of Wolbachia pipientis, which infects 40% of insect species and induces reproductive effects. The bacteria are passed from generation to generation both vertically (through the oocyte) and horizontally (by environmental transmission). Maintenance of the infection within Drosophila melanogaster is sensitive to the regulation of actin, as Wolbachia inefficiently colonizes strains hemizygous for the profilin or villin genes. Therefore, we hypothesized that Wolbachia must depend on the host actin cytoskeleton. In this study, we identify and characterize a Wolbachia protein (WD0830) that is predicted to be secreted by the bacterial parasite. Expression of WD0830 in a model eukaryote (the yeast Saccharomyces cerevisiae) induces a growth defect associated with the appearance of aberrant, filamentous structures which colocalize with rhodamine-phalloidin-stained actin. Purified WD0830 bundles actin in vitro and cosediments with actin filaments, suggesting a direct interaction of the two proteins. We characterized the expression of WD0830 throughout Drosophila development and found it to be upregulated in third-instar larvae, peaking in early pupation, during the critical formation of adult tissues, including the reproductive system. In transgenic flies, heterologously expressed WD0830 localizes to the developing oocyte. Additionally, overexpression of WD0830 results in increased Wolbachia titers in whole flies, in stage 9 and 10 oocytes, and in embryos, compared to controls, suggesting that the protein may facilitate Wolbachia's replication or transmission. Therefore, this candidate secreted effector may play a role in Wolbachia's infection of and persistence within host niches.

IMPORTANCE

The obligate intracellular Wolbachia pipientis is a ubiquitous alphaproteobacterial symbiont of arthropods and nematodes and is related to the rickettsial pathogens Ehrlichia spp. and Anaplasma spp. Studies of Wolbachia cell biology suggest that this bacterium relies on host actin for efficient proliferation and transmission between generations. Here, we identified and characterized a Wolbachia protein that localizes to and manipulates the eukaryotic actin cytoskeleton, is expressed by Wolbachia during host development, and alters Wolbachia titers and localization in transgenic fruit flies. We hypothesize that WD0830 may be utilized by the bacterium to facilitate replication in or invasion of different niches during host development.

Comparative Genomics of Two Closely Related Wolbachia with Different Reproductive Effects on Hosts

Wolbachia pipientis are obligate intracellular bacteria commonly found in many arthropods. They can induce various reproductive alterations in hosts, including cytoplasmic incompatibility, male-killing, feminization, and parthenogenetic development, and can provide host protection against some viruses and other pathogens. Wolbachia differ from many other primary endosymbionts in arthropods because they undergo frequent horizontal transmission between hosts and are well known for an abundance of mobile elements and relatively high recombination rates. Here, we compare the genomes of two closely related Wolbachia (with 0.57% genome-wide synonymous divergence) that differ in their reproductive effects on hosts. wVitA induces a sperm-egg incompatibility (also known as cytoplasmic incompatibility) in the parasitoid insect Nasonia vitripennis, whereas wUni causes parthenogenetic development in a different parasitoid, Muscidifurax uniraptor Although these bacteria are closely related, the genomic comparison reveals rampant rearrangements, protein truncations (particularly in proteins predicted to be secreted), and elevated substitution rates. These changes occur predominantly in the wUni lineage, and may be due in part to adaptations by wUni to a new host environment, or its phenotypic shift to parthenogenesis induction. However, we conclude that the approximately 8-fold elevated synonymous substitution rate in wUni is due to a either an elevated mutation rate or a greater number of generations per year in wUni, which occurs in semitropical host species. We identify a set of genes whose loss or pseudogenization in the wUni lineage implicates them in the phenotypic shift from cytoplasmic incompatibility to parthenogenesis induction. Finally, comparison of these closely related strains allows us to determine the fine-scale mutation patterns in Wolbachia Although Wolbachia are AT rich, mutation probabilities estimated from 4-fold degenerate sites are not AT biased, and predict an equilibrium AT content much less biased than observed (57-50% AT predicted vs. 76% current content at degenerate sites genome wide). The contrast suggests selection for increased AT content within Wolbachia genomes.

Dynamics of Wolbachia pipientis Gene Expression Across the Drosophila melanogaster Life Cycle

Symbiotic interactions between microbes and their multicellular hosts have manifold biological consequences. To better understand how bacteria maintain symbiotic associations with animal hosts, we analyzed genome-wide gene expression for the endosymbiotic α-proteobacteria Wolbachia pipientis across the entire life cycle of Drosophila melanogaster. We found that the majority of Wolbachia genes are expressed stably across the D. melanogaster life cycle, but that 7.8% of Wolbachia genes exhibit robust stage- or sex-specific expression differences when studied in the whole-organism context. Differentially-expressed Wolbachia genes are typically up-regulated after Drosophila embryogenesis and include many bacterial membrane, secretion system, and ankyrin repeat-containing proteins. Sex-biased genes are often organized as small operons of uncharacterized genes and are mainly up-regulated in adult Drosophila males in an age-dependent manner. We also systematically investigated expression levels of previously-reported candidate genes thought to be involved in host-microbe interaction, including those in the WO-A and WO-B prophages and in the Octomom region, which has been implicated in regulating bacterial titer and pathogenicity. Our work provides comprehensive insight into the developmental dynamics of gene expression for a widespread endosymbiont in its natural host context, and shows that public gene expression data harbor rich resources to probe the functional basis of the Wolbachia-Drosophila symbiosis and annotate the transcriptional outputs of the Wolbachia genome.

Microbial ecology-based methods to characterize the bacterial communities of non-model insects

Among the animals of the Kingdom Animalia, insects are unparalleled for their widespread diffusion, diversity and number of occupied ecological niches. In recent years they have raised researcher interest not only because of their importance as human and agricultural pests, disease vectors and as useful breeding species (e.g. honeybee and silkworm), but also because of their suitability as animal models. It is now fully recognized that microorganisms form symbiotic relationships with insects, influencing their survival, fitness, development, mating habits and the immune system and other aspects of the biology and ecology of the insect host. Thus, any research aimed at deepening the knowledge of any given insect species (perhaps species of applied interest or species emerging as novel pests or vectors) must consider the characterization of the associated microbiome. The present review critically examines the microbiology and molecular ecology techniques that can be applied to the taxonomical and functional analysis of the microbiome of non-model insects. Our goal is to provide an overview of current approaches and methods addressing the ecology and functions of microorganisms and microbiomes associated with insects. Our focus is on operational details, aiming to provide a concise guide to currently available advanced techniques, in an effort to extend insect microbiome research beyond simple descriptions of microbial communities.

The Wolbachia WO bacteriophage proteome in the Aedes albopictus C/wStr1 cell line: evidence for lytic activity?

Wolbachia pipientis (Rickettsiales), an obligate intracellular alphaproteobacterium in insects, manipulates host reproduction to maximize invasion of uninfected insect populations. Modification of host population structure has potential applications for control of pest species, particularly if Wolbachia can be maintained, manipulated, and genetically engineered in vitro. Although Wolbachia maintains an obligate mutualism with genome stability in nematodes, arthropods can be co-infected with distinct Wolbachia strains, and horizontal gene transfer between strains is potentially mediated by WO phages encoded within Wolbachia genomes. Proteomic analysis of a robust, persistent infection of a mosquito cell line with wStr from the planthopper, Laodelphax striatellus, revealed expression of a full array of WO phage genes, as well as nine of ten non-phage genes that occur between two distinct clusters of WOMelB genes in the genome of wMel, which infects Drosophila melanogaster. These non-phage genes encode potential host-adaptive proteins and are expressed in wStr at higher levels than phage structural proteins. A subset of seven of the non-phage genes is flanked by highly conserved non-coding sequences, including a putative promoter element, that are not present in a syntenically arranged array of homologs in plasmids from three tick-associated Rickettsia spp. These studies expand our understanding of wStr in a host cell line derived from the mosquito, Aedes albopictus, and provide a basis for investigating conditions that favor the lytic phase of the WO phage life cycle and recovery of infectious phage particles.

Passage of Wolbachia pipientis through mutant drosophila melanogaster induces phenotypic and genomic changes

Wolbachia pipientis is a nearly ubiquitous, maternally transmitted bacterium that infects the germ line of insect hosts. Estimates are that Wolbachia infects 40 to 60% of insect species on the planet, making it one of the most prevalent infections on Earth. However, we know surprisingly little about the molecular mechanisms used by Wolbachia to infect its hosts. We passaged Wolbachia through normally restrictive Drosophila melanogaster hosts, bottlenecking Wolbachia through stochastic segregation while simultaneously selecting for mutants that could recolonize these previously restrictive hosts. Here, we show that Wolbachia alters its behavior when passaged through heterozygous mutant flies. After only three generations, Wolbachia was able to colonize the previously restrictive hosts at control titers. Additionally, the Wolbachia organisms passaged through heterozygous mutant D. melanogaster alter their pattern of tissue-specific Wsp protein production, suggesting a behavioral response to the host genotype. Using whole-genome resequencing, we identified the mutations accumulated by these lineages of Wolbachia and confirmed the existence and persistence of the mutations through clone library Sanger sequencing. Our results suggest that Wolbachia can quickly adapt to new host contexts, with genomic mutants arising after only two generations.

Comparative genomics of first available bovine Anaplasma phagocytophilum genome obtained with targeted sequence capture

BACKGROUND

Anaplasma phagocytophilum is a zoonotic and obligate intracellular bacterium transmitted by ticks. In domestic ruminants, it is the causative agent of tick-borne fever, which causes significant economic losses in Europe. As A. phagocytophilum is difficult to isolate and cultivate, only nine genome sequences have been published to date, none of which originate from a bovine strain.Our goals were to; 1/ develop a sequencing methodology which efficiently circumvents the difficulties associated with A. phagocytophilum isolation and culture; 2/ describe the first genome of a bovine strain; and 3/ compare it with available genomes, in order to both explore key genomic features at the species level, and to identify candidate genes that could be specific to bovine strains.

RESULTS

DNA was extracted from a bovine blood sample infected by A. phagocytophilum. Following a whole genome capture approach, A. phagocytophilum DNA was enriched 197-fold in the sample and then sequenced using Illumina technology. In total, 58.9% of obtained reads corresponded to the A. phagocytophilum genome, covering 85.3% of the HZ genome. Then by performing comparisons with nine previously-sequenced A. phagocytophilum genomes, we determined the core genome of these ten strains. Following analysis, 1281 coding DNA sequences, including 1001 complete sequences, were detected in the A. phagocytophilum bovine genome, of which four appeared to be unique to the bovine isolate. These four coding DNA sequences coded for "hypothetical proteins of unknown function" and require further analysis. We also identified nine proteins common to both European domestic ruminants tested.

CONCLUSION

Using a whole genome capture approach, we have sequenced the first A. phagocytophilum genome isolated from a cow. To the best of our knowledge, this is the first time that this method has been used to selectively enrich pathogenic bacterial DNA from samples also containing host DNA. The four proteins unique to the A. phagocytophilum bovine genome could be involved in host tropism, therefore their functions need to be explored.

Proteomic profiling of a robust Wolbachia infection in an Aedes albopictus mosquito cell line

Wolbachia pipientis, a widespread vertically transmitted intracellular bacterium, provides a tool for insect control through manipulation of host-microbe interactions. We report proteomic characterization of wStr, a Wolbachia strain associated with a strong cytoplasmic incompatibility phenotype in its native host, Laodelphax striatellus. In the Aedes albopictus C/wStr1 mosquito cell line, wStr maintains a robust, persistent infection. MS/MS analyses of gel bands revealed a protein 'footprint' dominated by Wolbachia-encoded chaperones, stress response and cell membrane proteins, including the surface antigen WspA, a peptidoglycan-associated lipoprotein and a 73 kDa outer membrane protein. Functional classifications and estimated abundance levels of 790 identified proteins suggested that expression, stabilization and secretion of proteins predominate over bacterial genome replication and cell division. High relative abundances of cysteine desulphurase, serine/glycine hydroxymethyl transferase, and components of the α-ketoglutarate dehydrogenase complex in conjunction with above average abundances of glutamate dehydrogenase and proline utilization protein A support Wolbachia genome-based predictions for amino acid metabolism as a primary energy source. wStr expresses 15 Vir proteins of a Type IV secretion system and its transcriptional regulator. Proteomic characterization of a robust insect-associated Wolbachia strain provides baseline information that will inform further development of in vitro protocols for Wolbachia manipulation.

Wolbachia variants induce differential protection to viruses in Drosophila melanogaster: a phenotypic and phylogenomic analysis

Wolbachia are intracellular bacterial symbionts that are able to protect various insect hosts from viral infections. This tripartite interaction was initially described in Drosophila melanogaster carrying wMel, its natural Wolbachia strain. wMel has been shown to be genetically polymorphic and there has been a recent change in variant frequencies in natural populations. We have compared the antiviral protection conferred by different wMel variants, their titres and influence on host longevity, in a genetically identical D. melanogaster host. The phenotypes cluster the variants into two groups--wMelCS-like and wMel-like. wMelCS-like variants give stronger protection against Drosophila C virus and Flock House virus, reach higher titres and often shorten the host lifespan. We have sequenced and assembled the genomes of these Wolbachia, and shown that the two phenotypic groups are two monophyletic groups. We have also analysed a virulent and over-replicating variant, wMelPop, which protects D. melanogaster even better than the closely related wMelCS. We have found that a ~21 kb region of the genome, encoding eight genes, is amplified seven times in wMelPop and may be the cause of its phenotypes. Our results indicate that the more protective wMelCS-like variants, which sometimes have a cost, were replaced by the less protective but more benign wMel-like variants. This has resulted in a recent reduction in virus resistance in D. melanogaster in natural populations worldwide. Our work helps to understand the natural variation in wMel and its evolutionary dynamics, and inform the use of Wolbachia in arthropod-borne disease control.

Detection of the Wolbachia protein WPIP0282 in mosquito spermathecae: implications for cytoplasmic incompatibility

Cytoplasmic incompatibility (CI) is a conditional sterility induced by the bacterium Wolbachia pipientis that infects reproductive tissues in many arthropods. Although CI provides a potential tool to control insect vectors of arthropod-borne diseases, the molecular basis for CI induction is unknown. We hypothesized that a Wolbachia-encoded, CI-inducing factor would be enriched in sperm recovered from spermathecae of female mosquitoes. Using SDS-PAGE and mass spectrometry, we detected peptides from the 56 kDa hypothetical protein, encoded by wPip_0282, associated with sperm transferred to females by Wolbachia infected males. We also detected peptides from the same protein in Wolbachia infected ovaries. Homologs of wPip_0282 and the co-transcribed downstream gene, wPip_0283, occur as multiple divergent copies in genomes of CI-inducing strains of Wolbachia. The operon is located in a genomic context that includes mobile genetic elements. The absence of wPip_0282 and wPip_0283 homologs from genomes of Wolbachia in filarial nematodes, as well as other members of the Rickettsiales, suggests a role as a candidate CI effector.

Gene frequency distributions reject a neutral model of genome evolution

Evolution of prokaryotes involves extensive loss and gain of genes, which lead to substantial differences in the gene repertoires even among closely related organisms. Through a wide range of phylogenetic depths, gene frequency distributions in prokaryotic pangenomes bear a characteristic, asymmetrical U-shape, with a core of (nearly) universal genes, a “shell” of moderately common genes, and a “cloud” of rare genes. We employ mathematical modeling to investigate evolutionary processes that might underlie this universal pattern. Gene frequency distributions for almost 400 groups of 10 bacterial or archaeal species each over a broad range of evolutionary distances were fit to steady-state, infinite allele models based on the distribution of gene replacement rates and the phylogenetic tree relating the species in each group. The fits of the theoretical frequency distributions to the empirical ones yield model parameters and estimates of the goodness of fit. Using the Akaike Information Criterion, we show that the neutral model of genome evolution, with the same replacement rate for all genes, can be confidently rejected. Of the three tested models with purifying selection, the one in which the distribution of replacement rates is derived from a stochastic population model with additive per-gene fitness yields the best fits to the data. The selection strength estimated from the fits declines with evolutionary divergence while staying well outside the neutral regime. These findings indicate that, unlike some other universal distributions of genomic variables, for example, the distribution of paralogous gene family membership, the gene frequency distribution is substantially affected by selection.

Lateral transfers of insertion sequences between Wolbachia, Cardinium and Rickettsia bacterial endosymbionts

Various bacteria live exclusively within arthropod cells and collectively act as an important driver of arthropod evolutionary ecology. Whereas rampant intra-generic DNA transfers were recently shown to have a pivotal role in the evolution of the most common of these endosymbionts, Wolbachia, the present study show that inter-generic DNA transfers also commonly take place, constituting a potent source of rapid genomic change. Bioinformatic, molecular and phylogenetic data provide evidence that a selfish genetic element, the insertion sequence ISRpe1, is widespread in the Wolbachia, Cardinium and Rickettsia endosymbionts and experiences recent (and likely ongoing) transfers over long evolutionary distances. Although many ISRpe1 copies were clearly expanding and leading to rapid endosymbiont diversification, degraded copies are also frequently found, constituting an unusual genomic fossil record suggestive of ancient ISRpe1 expansions. Overall, the present data highlight how ecological connections within the arthropod intracellular environment facilitate lateral DNA transfers between distantly related bacterial lineages.

Draft genome sequence of the male-killing Wolbachia strain wBol1 reveals recent horizontal gene transfers from diverse sources

Background

The endosymbiont Wolbachia pipientis causes diverse and sometimes dramatic phenotypes in its invertebrate hosts. Four Wolbachia strains sequenced to date indicate that the constitution of the genome is dynamic, but these strains are quite divergent and do not allow resolution of genome diversification over shorter time periods. We have sequenced the genome of the strain wBol1-b, found in the butterfly Hypolimnas bolina, which kills the male offspring of infected hosts during embyronic development and is closely related to the non-male-killing strain wPip from Culex pipiens.

Results

The genomes of wBol1-b and wPip are similar in genomic organisation, sequence and gene content, but show substantial differences at some rapidly evolving regions of the genome, primarily associated with prophage and repetitive elements. We identified 44 genes in wBol1-b that do not have homologs in any previously sequenced strains, indicating that Wolbachia’s non-core genome diversifies rapidly. These wBol1-b specific genes include a number that have been recently horizontally transferred from phylogenetically distant bacterial taxa. We further report a second possible case of horizontal gene transfer from a eukaryote into Wolbachia.

Conclusions

Our analyses support the developing view that many endosymbiotic genomes are highly dynamic, and are exposed and receptive to exogenous genetic material from a wide range of sources. These data also suggest either that this bacterial species is particularly permissive for eukaryote-to-prokaryote gene transfers, or that these transfers may be more common than previously believed. The wBol1-b-specific genes we have identified provide candidates for further investigations of the genomic bases of phenotypic differences between closely-related Wolbachia strains.