Life and Death of Selfish Genes: Comparative Genomics Reveals the Dynamic Evolution of Cytoplasmic Incompatibility

Low Wolbachia incidence in Bactrocera and Zeugodacus species from Thailand and genome analysis of Wolbachia associated with Zeugodacus apicalis

Wolbachia are bacterial endosymbionts found widely in arthropods and filarial nematodes. Infecting about half of all arthropod species, Wolbachia manipulate their hosts in various ways, including cytoplasmic incompatibility. Here, we investigated Wolbachia diversity in Bactrocera and Zeugodacus, two prevalent tephritid fruit fly genera, using molecular methods. Wolbachia was only detected in Zeugodacus apicalis (de Meijere) (Diptera: Tephritidae) and not in the other 7 studied species. This newly discovered strain, named wZap, belongs to supergroup B with a 1.3 Mb genome containing 1,248 genes. Phylogenetic analysis of its cytoplasmic incompatibility factor genes cifA and cifB revealed their placement within the Type I clade. Given the presence of cif genes in the wZap genome, further research into their roles in fruit flies could be crucial for developing pest control strategies that exploit CI mechanisms.

Cardinium disrupts Wolbachia-host dynamics in the domestic mite Tyrophagus putrescentiae: evidence from manipulative experiments

We investigated the tripartite interactions between two intracellular bacterial symbionts, Cardinium and Wolbachia in Tyrophagus putrescentiae. Cultures of Tyrophagus putrescentiae are typically single-infected by one intracellular symbiont. However, co-infection can be experimentally induced by mixing single-infected cultures, resulting in 10% of mite individuals being double-infected (Cardinium + Wolbachia) and a corresponding reduction in host fitness. Here, we assembled the genomes of Cardinium and Wolbachia and analyzed their gene expression in parental single-infected and mixed mite cultures using population-level samples (ranging from 7,500 to 10,000 mites). Wolbachia interacts more extensively with its mite host than Cardinium in single-infected cultures. However, in mixed cultures, (i) Wolbachia exhibited reduced regulation of the host compared with Cardinium; (ii) the gene expression profile of Cardinium shifted, increasing its interactions with the host, whereas the gene expression profile of Wolbachia remained unchanged; and (iii) Wolbachia genes exhibited a loss of interactions with mite gene expression, as indicated by reduced correlations (for example with host MAPK, endocytosis, and calcium signaling pathways). The experiments show that at the mite population level, symbiont infection disrupts gene expression interaction between the two symbionts and their host in different ways. Wolbachia was more influenced by Cardinium gene expression than vice versa. Cardinium can inhibit the growth of Wolbachia by disrupting its interaction with the host, leading to a loss of Wolbachia's influence on mite immune and regulatory pathways. The reasons for responses are due to co-infection or the reduced frequency of Wolbachia single-infected individuals due to the analyses of population-level samples.

IMPORTANCE

We found that Cardinium disrupts the interaction between Wolbachia and mite host. In Wolbachia single-infected cultures, strong correlations exist between symbiont and host gene expressions. Interestingly, although Cardinium can also interact with the host, this interaction appears weaker compared with Wolbachia in single-infected cultures. These results suggest that both symbionts affect mite host gene expression, particularly in immune and regulatory pathways. In mixed samples, Cardinium appears to outcompete Wolbachia by disrupting its host interaction. It indicates competition between these two intracellular symbionts in mite populations. Wolbachia belongs to a mite-specific supergroup Q, distinct from the more commonly studied Wolbachia supergroups. As these mite-specific bacteria exhibit pathogen-blocking effects, our findings may have relevance for other systems, such as ticks and tick-borne diseases. The study sheds light on intracellular symbiont interaction within a novel mite-symbiont model.

Evolution of Wolbachia male-killing mechanism within a host species

Male-killing bacterial symbionts, prevalent in arthropods, skew population sex ratios by selectively killing male progeny, profoundly impacting ecology and the evolution of their hosts. Male killing is a convergently evolved trait, with microbes evolving diverse male-killing mechanisms across host species with widely divergent sex determination pathways. A common evolutionary response to male-killing presence is the spread of suppressor mutations that restore male survival. In this study, we demonstrate the evolution of a novel male-killing mechanism that is insensitive to an existing male-killing suppressor. Hypolimnas bolina butterflies from Yogyakarta, Indonesia, showed extreme female-biased population sex ratio associated with high prevalence of a male-killing Wolbachia. This strain, wBol1Y, shared a very recent common ancestor with the previously characterized "suppressed" male-killing strain in the species, wBol1, but it retained its male-killing ability in the presence of the male-killing suppressor. The genome of wBol1Y differed from the suppressed wBol1 in carrying an additional prophage that included strong candidate genes for male killing. In vitro and in vivo data demonstrated that wBol1Y feminized splicing and expression of lepidopteran sex determination pathway genes and that the gene Hb-oscar-present on wBol1Y's unique prophage insert-was sufficient to disrupt the male sex determination pathway. Our study demonstrates that the diversity of male-killing mechanisms is a product both of interaction with varying insect sex determination systems and the evolution of male killing within a host species. Our data indicate that the male killer and host may be involved in escalating arms races, where spreading male-killing suppression drives the evolution of additional systems that reestablish male killing by the symbiont.

Molecular detection of Wolbachia sp. and Cytoplasmic incompatibility factors (CifA/B) in wild caught mosquitoes in Côte d'Ivoire

BACKGROUND

Wolbachia is an endosymbiont bacterium known to stimulate host immunity against arboviruses and protozoa. Côte d'Ivoire is in a malaria-endemic region, and has experienced several dengue epidemics in recent decades as well. In order to help reduce the transmission of pathogens by mosquito vectors, we studied the prevalence of Wolbachia and the distribution of Cytoplasmic incompatibility factors (Cif) genes in different mosquito species caught in the wild in Cote d'Ivoire.

METHODS AND RESULTS

Mosquitoes of the genera Anopheles, Aedes, Culex, Eretmapodites and Mansonia were captured in five cities. Mosquitoes were collected at larval stage in breeding sites and adults were captured using BG sentinel traps. The mosquitoes were identified morphologically and Wolbachia and Cif were screened using qPCR targeting the 16s rRNA gene and the CifA, B genes. A total of 518 mosquito samples belonging to 15 species and 4 genera were examined. 60% of the species were infected with Wolbachia. The three medically important mosquito species Aedes aegypti, Anopheles gambiae s.l. and Culex quinquefasciatus had a prevalence of 12.84%, 13.46% and 72.64% respectively. The Wolbachia strains infecting the different mosquito species of the genus Culex encoded 98.46% for the CifA gene and 77.69% for the CifB gene.

CONCLUSION

The presence of Wolbachia and CifA, B genes in mosquitoes of different species in Côte d'Ivoire offer a promising opportunity to reduce the competence of mosquito vectors. Characterization of Wolbachia strains and cytoplasmic incompatibility factors will provide a better understanding of these endosymbionts, enabling the development of vector control strategies.

A phylogenetic approach to comparative genomics

Comparative genomics, whereby the genomes of different species are compared, has the potential to address broad and fundamental questions at the intersection of genetics and evolution. However, species, genomes and genes cannot be considered as independent data points within statistical tests. Closely related species tend to be similar because they share genes by common descent, which must be accounted for in analyses. This problem of non-independence may be exacerbated when examining genomes or genes but can be addressed by applying phylogeny-based methods to comparative genomic analyses. Here, we review how controlling for phylogeny can change the conclusions of comparative genomics studies. We address common questions on how to apply these methods and illustrate how they can be used to test causal hypotheses. The combination of rapidly expanding genomic datasets and phylogenetic comparative methods is set to revolutionize the biological insights possible from comparative genomic studies.

Recombination, truncation and horizontal transfer shape the diversity of cytoplasmic incompatibility patterns

Wolbachia are endosymbiotic bacteria inducing various reproductive manipulations of which cytoplasmic incompatibility (CI) is the most common. CI leads to reduced embryo viability in crosses between males carrying Wolbachia and uninfected females or those carrying an incompatible symbiont strain. In the mosquito Culex pipiens, the Wolbachia wPip causes highly complex crossing patterns. This complexity is linked to the amplification and diversification of the CI causal genes, cidA and cidB, with polymorphism located in the CidA-CidB interaction regions. We previously showed correlations between the identity of gene variants and CI patterns. However, these correlations were limited to specific crosses, and it is still unknown whether cid gene polymorphism in males' and females' Wolbachia can explain and predict the wide range of crossing types observed in C. pipiens. Taking advantage of a new method enabling full-gene acquisition, we sequenced complete cid repertoires from 45 wPip strains collected worldwide. We demonstrated that the extensive diversity of cid genes arises from recombination and horizontal transfers. We uncovered further cidB polymorphism outside the interface regions and strongly correlated with CI patterns. Most importantly, we showed that in every wPip genome, all but one cidB variant are truncated. Truncated cidBs located in palindromes are partially or completely deprived of their deubiquitinase domain, crucial for CI. The identity of the sole full-length cidB variant seems to dictate CI patterns, irrespective of the truncated cidBs present. Truncated CidBs exhibit reduced toxicity and stability in Drosophila cells, which potentially hinders their loading into sperm, essential for CI induction.

Wolbachia strain diversity in a complex group of sympatric cryptic parasitoid wasp species

BACKGROUND

Maternally-inherited symbionts can induce pre-mating and/or post-mating reproductive isolation between sympatric host lineages, and speciation, by modifying host reproductive phenotypes. The large parasitoid wasp genus Cotesia (Braconidae) includes a diversity of cryptic species, each specialized in parasitizing one to few related Lepidoptera host species. Here, we characterized the infection status of an assemblage of 21 Cotesia species from 15 countries by several microbial symbionts, as a first step toward investigating whether symbionts may provide a barrier to gene flow between these parasitoid host lineages.

RESULTS

The symbiotic microbes Arsenophonus, Cardinium, Microsporidium and Spiroplasma were not detected in the Cotesia wasps. However, the endosymbiotic bacterium Wolbachia was present in at least eight Cotesia species, and hence we concentrated on it upon screening additional DNA extracts and SRAs from NCBI. Some of the closely related Cotesia species carry similar Wolbachia strains, but most Wolbachia strains showed patterns of horizontal transfer between phylogenetically distant host lineages.

CONCLUSIONS

The lack of co-phylogenetic signal between Wolbachia and Cotesia suggests that the symbiont and hosts have not coevolved to an extent that would drive species divergence between the Cotesia host lineages. However, as the most common facultative symbiont of Cotesia species, Wolbachia may still function as a key-player in the biology of the parasitoid wasps. Its precise role in the evolution of this complex clade of cryptic species remains to be experimentally investigated.

Cell-based assays and comparative genomics revealed the conserved and hidden effects of Wolbachia on insect sex determination

It is advantageous for maternally transmitted endosymbionts to skew the sex ratio of their hosts toward females. Some endosymbiotic bacteria, such as Wolbachia, cause their insect hosts to exclusively produce female offspring through male killing (MK) or feminization. In some lepidopteran insects, MK is achieved by affecting the sex-determining process in males, and a unique mechanism of MK and its functional link with feminization have been implicated. However, comparative analysis of these phenotypes is often difficult because they have been analyzed in different host-symbiont systems, and transinfection of Wolbachia across different hosts is often challenging. In this study, we demonstrated the effects of nine Wolbachia strains on the splicing of sex-determining genes in Lepidoptera by fixing the host genetic background using a cell culture system. Cell transinfection assays confirmed that three MK-inducing Wolbachia strains and one feminization-inducing Wolbachia strain increased the female-type splicing products of the core sex-determining genes doublesex, masculinizer, and zinc finger protein 2. Regarding Wolbachia strains that do not induce MK/feminization, three had no effect on these sex-determining genes, whereas two strains induced female-type splicing of masculinizer and doublesex but not zinc finger protein 2. Comparative genomics confirmed that homologs of oscar, the Wolbachia gene responsible for MK in Ostrinia, were encoded by four MK/feminizing Wolbachia strains, but not by five non-MK/nonfeminizing strains. These results support the conserved effects underlying MK and feminization induced by oscar-bearing Wolbachia and suggested other potential mechanisms that Wolbachia might employ to manipulate host sex.

Genomic Underpinnings of Cytoplasmic Incompatibility: CIF Gene-Neighborhood Diversification Through Extensive Lateral Transfers and Recombination in Wolbachia

Cytoplasmic incompatibility (CI), a non-Mendelian genetic phenomenon, involves the manipulation of host reproduction by Wolbachia, a maternally transmitted alphaproteobacterium. The underlying mechanism is centered around the CI Factor (CIF) system governed by two genes, cifA and cifB, where cifB induces embryonic lethality, and cifA counteracts it. Recent investigations have unveiled intriguing facets of this system, including diverse cifB variants, prophage association in specific strains, copy number variation, and rapid component divergence, hinting at a complex evolutionary history. We utilized comparative genomics to systematically classify CIF systems, analyze their locus structure and domain architectures, and reconstruct their diversification and evolutionary trajectories. Our new classification identifies ten distinct CIF types, featuring not just versions present in Wolbachia, but also other intracellular bacteria, and eukaryotic hosts. Significantly, our analysis of CIF loci reveals remarkable variability in gene composition and organization, encompassing an array of diverse endonucleases, variable toxin domains, deubiquitinating peptidases (DUBs), prophages, and transposons. We present compelling evidence that the components within the loci have been diversifying their sequences and domain architectures through extensive, independent lateral transfers and interlocus recombination involving gene conversion. The association with diverse transposons and prophages, coupled with selective pressures from host immunity, likely underpins the emergence of CIF loci as recombination hotspots. Our investigation also posits the origin of CifB-REase domains from mobile elements akin to CR (Crinkler-RHS-type) effectors and Tribolium Medea1 factor, which is linked to another non-Mendelian genetic phenomenon. This comprehensive genomic analysis offers novel insights into the molecular evolution and genomic foundations of Wolbachia-mediated host reproductive control.

Culex-Transmitted Diseases: Mechanisms, Impact, and Future Control Strategies using Wolbachia

Mosquitoes of the Culex genus are responsible for a large burden of zoonotic virus transmission globally. Collectively, they play a significant role in the transmission of medically significant diseases such as Japanese encephalitis virus and West Nile virus. Climate change, global trade, habitat transformation and increased urbanisation are leading to the establishment of Culex mosquitoes in new geographical regions. These novel mosquito incursions are intensifying concerns about the emergence of Culex-transmitted diseases and outbreaks in previously unaffected areas. New mosquito control methods are currently being developed and deployed globally. Understanding the complex interaction between pathogens and mosquitoes is essential for developing new control strategies for Culex species mosquitoes. This article reviews the role of Culex mosquitos as vectors of zoonotic disease, discussing the transmission of viruses across different species, and the potential use of Wolbachia technologies to control disease spread. By leveraging the insights gained from recent successful field trials of Wolbachia against Aedes-borne diseases, we comprehensively discuss the feasibility of using this technique to control Culex mosquitoes and the potential for the development of next generational Wolbachia-based control methods.

Describing endosymbiont-host interactions within the parasitism-mutualism continuum

Endosymbionts are widespread in arthropods, living in host cells with effects that extend from parasitic to mutualistic. Newly acquired endosymbionts tend to be parasitic, but vertical transmission favors coevolution toward mutualism, with hosts sometimes developing dependency. Endosymbionts negatively affecting host fitness may still spread by impacting host reproductive traits, referred to as reproductive "manipulation," although costs for hosts are often assumed rather than demonstrated. For cytoplasmic incompatibility (CI) that involves endosymbiont-mediated embryo death, theory predicts directional shifts away from "manipulation" toward reduced CI strength; moreover, CI-causing endosymbionts need to increase host fitness to initially spread. In nature, endosymbiont-host interactions and dynamics are complex, often depending on environmental conditions and evolutionary history. We advocate for capturing this complexity through appropriate datasets, rather than relying on terms like "manipulation." Such imprecision can lead to the misclassification of endosymbionts along the parasitism-mutualism continuum.

Rickettsia induces strong cytoplasmic incompatibility in a predatory insect

Rickettsia, a group of intracellular bacteria found in eukaryotes, exhibits diverse lifestyles, with some acting as vertebrate pathogens transmitted by arthropod vectors and others serving as maternally transmitted arthropod endosymbionts, some of which manipulate host reproduction for their own benefit. Two phenotypes, namely male-killing and parthenogenesis induction are known as Rickettsia-induced host reproductive manipulations, but it remains unknown whether Rickettsia can induce other types of host manipulation. In this study, we discovered that Rickettsia induced strong cytoplasmic incompatibility (CI), in which uninfected females produce no offspring when mated with infected males, in the predatory insect Nesidiocoris tenuis (Hemiptera: Miridae). Molecular phylogenetic analysis revealed that the Rickettsia strain was related to Rickettsia bellii, a common insect endosymbiont. Notably, this strain carried plasmid-encoded homologues of the CI-inducing factors (namely cifA-like and cifB-like genes), typically found in Wolbachia, which are well-known CI-inducing endosymbionts. Protein domain prediction revealed that the cifB-like gene encodes PD-(D/E)XK nuclease and deubiquitinase domains, which are responsible for Wolbachia-induced CI, as well as ovarian tumour-like (OTU-like) cysteine protease and ankyrin repeat domains. These findings suggest that Rickettsia and Wolbachia endosymbionts share underlying mechanisms of CI and that CI-inducing ability was acquired by microbes through horizontal plasmid transfer.

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

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

Comparative analysis of Wolbachia maternal transmission and localization in host ovaries

Many insects and other animals carry microbial endosymbionts that influence their reproduction and fitness. These relationships only persist if endosymbionts are reliably transmitted from one host generation to the next. Wolbachia are maternally transmitted endosymbionts found in most insect species, but transmission rates can vary across environments. Maternal transmission of wMel Wolbachia depends on temperature in natural Drosophila melanogaster hosts and in transinfected Aedes aegypti, where wMel is used to block pathogens that cause human disease. In D. melanogaster, wMel transmission declines in the cold as Wolbachia become less abundant in host ovaries and at the posterior pole plasm (the site of germline formation) in mature oocytes. Here, we assess how temperature affects maternal transmission and underlying patterns of Wolbachia localization across 10 Wolbachia strains diverged up to 50 million years-including strains closely related to wMel-and their natural Drosophila hosts. Many Wolbachia maintain high transmission rates across temperatures, despite highly variable (and sometimes low) levels of Wolbachia in the ovaries and at the developing germline in late-stage oocytes. Identifying strains like closely related wMel-like Wolbachia with stable transmission across variable environmental conditions may improve the efficacy of Wolbachia-based biocontrol efforts as they expand into globally diverse environments.

Wolbachia symbionts control sex in a parasitoid wasp using a horizontally acquired gene

Host reproduction can be manipulated by bacterial symbionts in various ways. Parthenogenesis induction is the most effective type of reproduction manipulation by symbionts for their transmission. Insect sex is determined by regulation of doublesex (dsx) splicing through transformer2 (tra2) and transformer (tra) interaction. Although parthenogenesis induction by symbionts has been studied since the 1970s, its underlying molecular mechanism is unknown. Here we identify a Wolbachia parthenogenesis-induction feminization factor gene (piff) that targets sex-determining genes and causes female-producing parthenogenesis in the haplodiploid parasitoid Encarsia formosa. We found that Wolbachia elimination repressed expression of female-specific dsx and enhanced expression of male-specific dsx, which led to the production of wasp haploid male offspring. Furthermore, we found that E. formosa tra is truncated and non-functional, and Wolbachia has a functional tra homolog, termed piff, with an insect origin. Wolbachia PIFF can colocalize and interact with wasp TRA2. Moreover, Wolbachia piff has coordinated expression with tra2 and dsx of E. formosa. Our results demonstrate the bacterial symbiont Wolbachia has acquired an insect gene to manipulate the host sex determination cascade and induce parthenogenesis in wasps. This study reveals insect-to-bacteria horizontal gene transfer drives the evolution of animal sex determination systems, elucidating a striking mechanism of insect-microbe symbiosis.

Wolbachia strain wAlbB shows favourable characteristics for dengue control use in Aedes aegypti from Burkina Faso

Dengue represents an increasing public health burden worldwide. In Africa, underreporting and misdiagnosis often mask its true epidemiology, and dengue is likely to be both more widespread than reported data suggest and increasing in incidence and distribution. Wolbachia-based dengue control is underway in Asia and the Americas but has not to date been deployed in Africa. Due to the genetic heterogeneity of African Aedes aegypti populations and the complexity of the host-symbiont interactions, characterization of key parameters of Wolbachia-carrying mosquitoes is paramount for determining the potential of the system as a control tool for dengue in Africa. The wAlbB Wolbachia strain was stably introduced into an African Ae. aegypti population by introgression, and showed high intracellular density in whole bodies and different mosquito tissues; high intracellular density was also maintained following larval rearing at high temperatures. No effect on the adult lifespan induced by Wolbachia presence was detected. Moreover, the ability of this strain to strongly inhibit DENV-2 dissemination and transmission in the host was also demonstrated in the African background. Our findings suggest the potential of harnessing Wolbachia for dengue control for African populations of Ae. aegypti.

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

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

Intra-lineage microevolution of Wolbachia leads to the emergence of new cytoplasmic incompatibility patterns

Mosquitoes of the Culex pipiens complex are worldwide vectors of arbovirus, filarial nematodes, and avian malaria agents. In these hosts, the endosymbiotic bacteria Wolbachia induce cytoplasmic incompatibility (CI), i.e., reduced embryo viability in so-called incompatible crosses. Wolbachia infecting Culex pipiens (wPip) cause CI patterns of unparalleled complexity, associated with the amplification and diversification of cidA and cidB genes, with up to 6 different gene copies described in a single wPip genome. In wPip, CI is thought to function as a toxin-antidote (TA) system where compatibility relies on having the right antidotes (CidA) in the female to bind and neutralize the male's toxins (CidB). By repeating crosses between Culex isofemale lines over a 17 years period, we documented the emergence of a new compatibility type in real time and linked it to a change in cid genes genotype. We showed that loss of specific cidA gene copies in some wPip genomes results in a loss of compatibility. More precisely, we found that this lost antidote had an original sequence at its binding interface, corresponding to the original sequence at the toxin's binding interface. We showed that these original cid variants are recombinant, supporting a role for recombination rather than point mutations in rapid CI evolution. These results strongly support the TA model in natura, adding to all previous data acquired with transgenes expression.

Transgenic expression of cif genes from Wolbachia strain wAlbB recapitulates cytoplasmic incompatibility in Aedes aegypti

The endosymbiotic bacteria Wolbachia can invade insect populations by modifying host reproduction through cytoplasmic incompatibility (CI), an effect that results in embryonic lethality when Wolbachia-carrying males mate with Wolbachia-free females. Here we describe a transgenic system for recreating CI in the major arbovirus vector Aedes aegypti using CI factor (cif) genes from wAlbB, a Wolbachia strain currently being deployed to reduce dengue transmission. CI-like sterility is induced when cifA and cifB are co-expressed in testes; this sterility is rescued by maternal cifA expression, thereby reproducing the pattern of Wolbachia-induced CI. Expression of cifB alone is associated with extensive DNA damage and disrupted spermatogenesis. The strength of rescue by maternal cifA expression is dependent on the comparative levels of cifA/cifB expression in males. These findings are consistent with CifB acting as a toxin and CifA as an antitoxin, with CifA attenuating CifB toxicity in both the male germline and in developing embryos. These findings provide important insights into the interactions between cif genes and their mechanism of activity and provide a foundation for the building of a cif gene-based drive system in Ae. aegypti.

Detection and Analysis of Wolbachia in Plant-Parasitic Nematodes and Insights into Wolbachia Evolution

Since the discovery of Wolbachia in plant-parasitic nematodes (PPNs), there has been increased interest in this earliest branching clade that may hold important clues to early transitions in Wolbachia function in the Ecdysozoa. However, due to the specialized skills and equipment of nematology and the difficulty in culturing most PPNs, these PPN-type Wolbachia remain undersampled and poorly understood. To date, there are few established laboratory methods for working with PPN-type Wolbachia strains, and most research has relied on chance discovery and comparative genomics. Here, we address this challenge by providing detailed methods to assist researchers with more efficiently collecting PPNs and screen these communities, populations, or single nematodes with a newly developed PPN-type Wolbachia-specific PCR assay. We provide an overview of the typical yields and outcomes of these methods, to facilitate further targeted cultivation or experimental methods, and finally we provide a short introduction to some of the specific challenges and solutions in following through with comparative or population genomics on PPN-type Wolbachia strains.

NDUFA8 potentially rescues Wolbachia-induced cytoplasmic incompatibility in Laodelphax striatellus

The endosymbiont Wolbachia manipulates host reproduction by several strategies, one of the most important of which is cytoplasmic incompatibility (CI). CI can be rescued when Wolbachia-infected males mate with females infected with the same Wolbachia strain. However, the potential rescue mechanism of CI in the small brown planthopper Laodelphax striatellus is unclear. In this study, comparative transcriptome analysis was applied to explore the effect of Wolbachia on L. striatellus eggs. A total of 1387 differentially expressed genes were identified. RNA interference of 7 Wolbachia-upregulated key planthopper genes reduced egg reproduction, suggesting that Wolbachia might improve fecundity in L. striatellus by affecting these 7 genes. Suppressing the expression of another upregulated gene, NDUFA8 (encoding NADH dehydrogenase [ubiquinone] 1 α subcomplex subunit 8-like) by RNA interference significantly increased the mortality of early embryos without affecting the number of deposited eggs. Wolbachia infection upregulated the mRNA level of NDUFA8, and dsNDUFA8 treatment of Wolbachia-infected females recreated CI-like symptoms, suggesting that NDUFA8 is associated with the rescue phenotype. Because all L. striatellus populations worldwide are infected with Wolbachia, NDUFA8 is a potential pest control target.

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

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

Cardinium symbionts are pervasive in Iranian populations of the spider mite Panonychus ulmi despite inducing an infection cost and no demonstrable reproductive phenotypes when Wolbachia is a symbiotic partner

Maternally transmitted symbionts such as Cardinium and Wolbachia are widespread in arthropods. Both Cardinium and Wolbachia can cause cytoplasmic incompatibility, a reproductive phenotype that interferes with the development of uninfected eggs that are fertilized by infected sperm. In haplodiploid hosts, these symbionts can also distort sex allocation to facilitate their spread through host populations. Without other fitness effects, symbionts that induce strong reproductive phenotypes tend to spread to high and stable infection frequencies, whereas variants that induce weak reproductive phenotypes are typically associated with intermediate and variable frequencies. To study the spread of Cardinium in a haplodiploid host, we sampled Iranian populations of the economically important spider mite Panonychus ulmi in apple orchards. Within several field populations, we also studied the Wolbachia infection frequencies. All P. ulmi field populations carried a Cardinium infection and exhibited high infection frequencies. In contrast, Wolbachia frequency ranged between ca. 10% and ca. 70% and was only found in co-infected mites. To test whether Cardinium induce reproductive phenotypes in P. ulmi, a Cardinium-cured derived line was generated by antibiotic treatment from a co-infected field population. Genetic crosses indicated that Cardinium do not induce demonstrable levels of cytoplasmic incompatibility and sex allocation distortion in co-infected P. ulmi. However, Cardinium infection was associated with a longer developmental time and reduced total fecundity for co-infected females. We hypothesize that Cardinium spread through P. ulmi populations via uncharacterized fitness effects and that co-infection with Wolbachia might impact these drive mechanisms.

Transinfected Wolbachia strains induce a complex of cytoplasmic incompatibility phenotypes: Roles of CI factor genes

Wolbachia can modulate the reproductive development of their hosts in multiple modes, and cytoplasmic incompatibility (CI) is the most well-studied phenotype. The whitefly Bemisia tabaci is highly receptive to different Wolbachia strains: wCcep strain from the rice moth Corcyra cephalonica and wMel strain from the fruit fly Drosophila melanogaster could successfully establish and induce CI in transinfected whiteflies. Nevertheless, it is unknown what will happen when these two exogenous Wolbachia strains are co-transinfected into a new host. Here, we artificially transinferred wCcep and wMel into the whitefly and established double- and singly-transinfected B. tabaci isofemale lines. Reciprocal crossing experiments showed that wCcep and wMel induced a complex of CI phenotypes in the recipient host, including unidirectional and bidirectional CI. We next sequenced the whole genome of wCcep and performed a comparative analysis of the CI factor genes between wCcep and wMel, indicating that their cif genes were phylogenetically and structurally divergent, which can explain the crossing results. The amino acid sequence identity and structural features of Cif proteins may be useful parameters for predicting their function. Structural comparisons between CifA and CifB provide valuable clues for explaining the induction or rescue of CI observed in crossing experiments between transinfected hosts.

Complete De Novo Assembly of Wolbachia Endosymbiont of Frankliniella intonsa

As an endosymbiont, Wolbachia exerts significant effects on the host, including on reproduction, immunity, and metabolism. However, the study of Wolbachia in Thysanopteran insects, such as flower thrips Frankliniella intonsa, remains limited. Here, we assembled a gap-free looped genome assembly of Wolbachia strain wFI in a length of 1,463,884 bp (GC content 33.80%), using Nanopore long reads and Illumina short reads. The annotation of wFI identified a total of 1838 protein-coding genes (including 85 pseudogenes), 3 ribosomal RNAs (rRNAs), 35 transfer RNAs (tRNAs), and 1 transfer-messenger RNA (tmRNA). Beyond this basic description, we identified mobile genetic elements, such as prophage and insertion sequences (ISs), which make up 17% of the entire wFI genome, as well as genes involved in riboflavin and biotin synthesis and metabolism. This research lays the foundation for understanding the nutritional mutualism between Wolbachia and flower thrips. It also serves as a valuable resource for future studies delving into the intricate interactions between Wolbachia and its host.

Metagenomics of culture isolates and insect tissue illuminate the evolution of Wolbachia, Rickettsia and Bartonella symbionts in Ctenocephalides spp. fleas

While fleas are often perceived simply as a biting nuisance and a cause of allergic dermatitis, they represent important disease vectors worldwide, especially for bacterial zoonoses such as plague (transmitted by rodent fleas) and some of the rickettsioses and bartonelloses. The cosmopolitan cat (Ctenocephalides felis ) and dog (Ctenocephalides canis ) fleas, as well as Ctenocephalides orientis (restricted to tropical and subtropical Asia), breed in human dwellings and are vectors of cat-scratch fever (caused by Bartonella spp.) and Rickettsia spp., including Rickettsia felis (agent of flea-borne spotted fever) and Rickettsia asembonensis , a suspected pathogen. These Rickettsia spp. are members of a phylogenetic clade known as the ‘transitional group’, which includes both human pathogens and arthropod-specific endosymbionts. The relatively depauperate flea microbiome can also contain other endosymbionts, including a diverse range of Wolbachia strains. Here, we present circularized genome assemblies for two C. orientis -derived pathogens (Bartonella clarridgeiae and R. asembonensis ) from Malaysia, a novel Wolbachia strain (w Cori), and the C. orientis mitochondrion; all were obtained by direct metagenomic sequencing of flea tissues. Moreover, we isolated two Wolbachia strains from Malaysian C. felis into tick cell culture and recovered circularized genome assemblies for both, one of which (w CfeF) is newly sequenced. We demonstrate that the three Wolbachia strains are representatives of different major clades (‘supergroups’), two of which appear to be flea-specific. These Wolbachia genomes exhibit unique combinations of features associated with reproductive parasitism or mutualism, including prophage WO, cytoplasmic incompatibility factors and the biotin operon of obligate intracellular microbes. The first circularized assembly for R. asembonensis includes a plasmid with a markedly different structure and gene content compared to the published plasmid; moreover, this novel plasmid was also detected in cat flea metagenomes from the USA. Analysis of loci under positive selection in the transitional group revealed genes involved in host–pathogen interactions that may facilitate host switching. Finally, the first B. clarridgeiae genome from Asia exhibited large-scale genome stability compared to isolates from other continents, except for SNPs in regions predicted to mediate interactions with the vertebrate host. These findings highlight the paucity of data on the genomic diversity of Ctenocephalides -associated bacteria and raise questions regarding how interactions between members of the flea microbiome might influence vector competence.

Modeling emergence of Wolbachia toxin-antidote protein functions with an evolutionary algorithm

Evolutionary algorithms (EAs) simulate Darwinian evolution and adeptly mimic natural evolution. Most EA applications in biology encode high levels of abstraction in top-down population ecology models. In contrast, our research merges protein alignment algorithms from bioinformatics into codon based EAs that simulate molecular protein string evolution from the bottom up. We apply our EA to reconcile a problem in the field of Wolbachia induced cytoplasmic incompatibility (CI). Wolbachia is a microbial endosymbiont that lives inside insect cells. CI is conditional insect sterility that operates as a toxin antidote (TA) system. Although, CI exhibits complex phenotypes not fully explained under a single discrete model. We instantiate in-silico genes that control CI, CI factors (cifs), as strings within the EA chromosome. We monitor the evolution of their enzymatic activity, binding, and cellular localization by applying selective pressure on their primary amino acid strings. Our model helps rationalize why two distinct mechanisms of CI induction might coexist in nature. We find that nuclear localization signals (NLS) and Type IV secretion system signals (T4SS) are of low complexity and evolve fast, whereas binding interactions have intermediate complexity, and enzymatic activity is the most complex. Our model predicts that as ancestral TA systems evolve into eukaryotic CI systems, the placement of NLS or T4SS signals can stochastically vary, imparting effects that might impact CI induction mechanics. Our model highlights how preconditions and sequence length can bias evolution of cifs toward one mechanism or another.

The toxins of vertically transmitted Spiroplasma

Vertically transmitted (VT) microbial symbionts play a vital role in the evolution of their insect hosts. A longstanding question in symbiont research is what genes help promote long-term stability of vertically transmitted lifestyles. Symbiont success in insect hosts is due in part to expression of beneficial or manipulative phenotypes that favor symbiont persistence in host populations. In Spiroplasma, these phenotypes have been linked to toxin and virulence domains among a few related strains. However, these domains also appear frequently in phylogenetically distant Spiroplasma, and little is known about their distribution across the Spiroplasma genus. In this study, we present the complete genome sequence of the Spiroplasma symbiont of Drosophila atripex, a non-manipulating member of the Ixodetis clade of Spiroplasma, for which genomic data are still limited. We perform a genus-wide comparative analysis of toxin domains implicated in defensive and reproductive phenotypes. From 12 VT and 31 non-VT Spiroplasma genomes, ribosome-inactivating proteins (RIPs), OTU-like cysteine proteases (OTUs), ankyrins, and ETX/MTX2 domains show high propensity for VT Spiroplasma compared to non-VT Spiroplasma. Specifically, OTU and ankyrin domains can be found only in VT-Spiroplasma, and RIP domains are found in all VT Spiroplasma and three non-VT Spiroplasma. These domains are frequently associated with Spiroplasma plasmids, suggesting a possible mechanism for dispersal and maintenance among heritable strains. Searching insect genome assemblies available on public databases uncovered uncharacterized Spiroplasma genomes from which we identified several spaid-like genes encoding RIP, OTU, and ankyrin domains, suggesting functional interactions among those domain types. Our results suggest a conserved core of symbiont domains play an important role in the evolution and persistence of VT Spiroplasma in insects.

Variation in sex determination mechanisms may constrain parthenogenesis-induction by endosymbionts in haplodiploid systems

Endosymbionts are maternally transmitted, and therefore benefit from maximizing female offspring numbers. Parthenogenesis-induction (PI) is the most effective type of manipulation for transmission, but has solely been detected in haplodiploid species, whereas cytoplasmic incompatibility (CI) is detected frequently across the arthropod phylum, including haplodiploids. This puzzling observation led us to hypothesize that the molecular sex-determination mechanism of the haplodiploid host may be a constraining factor in the ability of endosymbionts to induce parthenogenesis. Recent insights indicate that PI-endosymbionts may be able to directly manipulate sex-determination genes to induce the necessary steps required for PI in haplodiploids. However, sex-determination cascades vary extensively, so PI-induction would require a specialized and host-dependent tool set. Contrastingly, CI-related genes target conserved cell-cycle mechanisms, are located on mobile elements, and spread easily. Finally, endosymbiont-manipulations may have a strong impact on the effectiveness of haplodiploid biocontrol agents, but can also be used to enhance their efficacy.

Modelling Emergence of Wolbachia Toxin-Antidote Protein Functions with an Evolutionary Algorithm

Evolutionary algorithms (EAs) simulate Darwinian evolution and adeptly mimic natural evolution. Most EA applications in biology encode high levels of abstraction in top-down ecological population models. In contrast, our research merges protein alignment algorithms from bioinformatics into codon based EAs that simulate molecular protein string evolution from the bottom up. We apply our EA to reconcile a problem in the field of Wolbachia induced cytoplasmic incompatibility (CI). Wolbachia is a microbial endosymbiont that lives inside insect cells. CI is conditional insect sterility that operates as a toxin antidote (TA) system. Although, CI exhibits complex phenotypes not fully explained under a single discrete model. We instantiate in-silico genes that control CI, CI factors ( cifs ), as strings within the EA chromosome. We monitor the evolution of their enzymatic activity, binding, and cellular localization by applying selective pressure on their primary amino acid strings. Our model helps rationalize why two distinct mechanisms of CI induction might coexist in nature. We find that nuclear localization signals (NLS) and Type IV secretion system signals (T4SS) are of low complexity and evolve fast, whereas binding interactions have intermediate complexity, and enzymatic activity is the most complex. Our model predicts that as ancestral TA systems evolve into eukaryotic CI systems, the placement of NLS or T4SS signals can stochastically vary, imparting effects that might impact CI induction mechanics. Our model highlights how preconditions, genetic diversity, and sequence length can bias evolution of cifs towards one mechanism or another.

A male-killing Wolbachia endosymbiont is concealed by another endosymbiont and a nuclear suppressor

Bacteria that live inside the cells of insect hosts (endosymbionts) can alter the reproduction of their hosts, including the killing of male offspring (male killing, MK). MK has only been described in a few insects, but this may reflect challenges in detecting MK rather than its rarity. Here, we identify MK Wolbachia at a low frequency (around 4%) in natural populations of Drosophila pseudotakahashii. MK Wolbachia had a stable density and maternal transmission during laboratory culture, but the MK phenotype which manifested mainly at the larval stage was lost rapidly. MK Wolbachia occurred alongside a second Wolbachia strain expressing a different reproductive manipulation, cytoplasmic incompatibility (CI). A genomic analysis highlighted Wolbachia regions diverged between the 2 strains involving 17 genes, and homologs of the wmk and cif genes implicated in MK and CI were identified in the Wolbachia assembly. Doubly infected males induced CI with uninfected females but not females singly infected with CI-causing Wolbachia. A rapidly spreading dominant nuclear suppressor genetic element affecting MK was identified through backcrossing and subsequent analysis with ddRAD SNPs of the D. pseudotakahashii genome. These findings highlight the complexity of nuclear and microbial components affecting MK endosymbiont detection and dynamics in populations and the challenges of making connections between endosymbionts and the host phenotypes affected by them.

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

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

Genomic and Phenotypic Comparisons Reveal Distinct Variants of Wolbachia Strain wAlbB

The intracellular bacterium Wolbachia inhibits virus replication and is being harnessed around the world to fight mosquito-borne diseases through releases of mosquitoes carrying the symbiont. Wolbachia strains vary in their ability to invade mosquito populations and suppress viruses in part due to differences in their density within the insect and associated fitness costs. Using whole-genome sequencing, we demonstrate the existence of two variants in wAlbB, a Wolbachia strain being released in natural populations of Aedes aegypti mosquitoes. The two variants display striking differences in genome architecture and gene content. Differences in the presence/absence of 52 genes between variants include genes located in prophage regions and others potentially involved in controlling the symbiont's density. Importantly, we show that these genetic differences correlate with variation in wAlbB density and its tolerance to heat stress, suggesting that different wAlbB variants may be better suited for field deployment depending on local environmental conditions. Finally, we found that the wAlbB genome remained stable following its introduction in a Malaysian mosquito population. Our results highlight the need for further genomic and phenotypic characterization of Wolbachia strains in order to inform ongoing Wolbachia-based programs and improve the selection of optimal strains in future field interventions. IMPORTANCE Dengue is a viral disease transmitted by Aedes mosquitoes that threatens around half of the world population. Recent advances in dengue control involve the introduction of Wolbachia bacterial symbionts with antiviral properties into mosquito populations, which can lead to dramatic decreases in the incidence of the disease. In light of these promising results, there is a crucial need to better understand the factors affecting the success of such strategies, in particular the choice of Wolbachia strain for field releases and the potential for evolutionary changes. Here, we characterized two variants of a Wolbachia strain used for dengue control that differ at the genomic level and in their ability to replicate within the mosquito. We also found no evidence for the evolution of the symbiont within the 2 years following its deployment in Malaysia. Our results have implications for current and future Wolbachia-based health interventions.

Wolbachia wAlbB remains stable in Aedes aegypti over 15 years but exhibits genetic background-dependent variation in virus blocking

The ability of the maternally transmitted endosymbiotic bacterium Wolbachia to induce cytoplasmic incompatibility (CI) and virus blocking makes it a promising weapon for combatting mosquito-borne diseases through either suppression or replacement of wild-type populations. Recent field trials show that both approaches significantly reduce the incidence of dengue fever in humans. However, new questions emerge about how Wolbachia-mosquito associations will co-evolve over time and whether Wolbachia-mediated virus blocking will be affected by the genetic diversity of mosquitoes and arboviruses in the real world. Here, we have compared the Wolbachia density and CI expression of two wAlbB-infected Aedes aegypti lines transinfected 15 years apart. We have also assessed wAlbB-mediated virus blocking against dengue (DENV), Zika (ZIKV), and Chikungunya (CHIKV) viruses and examined whether host genetic backgrounds modulate viral blocking effects by comparing ZIKV infection in mosquitoes with a Mexican genetic background to those with a Singaporean background. Our results show that over 15 years, wAlbB maintained the capacity to form a stable association with Ae. aegypti in terms of both density and CI expression. There were variations in wAlbB-induced virus blocking against CHIKV, DENV, and ZIKV, and higher inhibitory effects on ZIKV in mosquitoes on the Singaporean genetic background than on the Mexican background. These results provide important information concerning the robustness and long-term stability of Wolbachia as a biocontrol agent for arbovirus disease control.

Genome sequencing and comparative analysis of Wolbachia strain wAlbA reveals Wolbachia-associated plasmids are common

Wolbachia are widespread maternally-transmitted bacteria of arthropods that often spread by manipulating their host's reproduction through cytoplasmic incompatibility (CI). Their invasive potential is currently being harnessed in field trials aiming to control mosquito-borne diseases. Wolbachia genomes commonly harbour prophage regions encoding the cif genes which confer their ability to induce CI. Recently, a plasmid-like element was discovered in wPip, a Wolbachia strain infecting Culex mosquitoes; however, it is unclear how common such extra-chromosomal elements are in Wolbachia. Here we sequenced the complete genome of wAlbA, a strain of the symbiont found in Aedes albopictus, after eliminating the co-infecting and higher density wAlbB strain that previously made sequencing of wAlbA challenging. We show that wAlbA is associated with two new plasmids and identified additional Wolbachia plasmids and related chromosomal islands in over 20% of publicly available Wolbachia genome datasets. These plasmids encode a variety of accessory genes, including several phage-like DNA packaging genes as well as genes potentially contributing to host-symbiont interactions. In particular, we recovered divergent homologues of the cif genes in both Wolbachia- and Rickettsia-associated plasmids. Our results indicate that plasmids are common in Wolbachia and raise fundamental questions around their role in symbiosis. In addition, our comparative analysis provides useful information for the future development of genetic tools to manipulate and study Wolbachia symbionts.

Toward an accurate mechanistic understanding of Wolbachia-induced cytoplasmic incompatibility

Wolbachia are the most successful intracellular bacteria in arthropods. They can manipulate host reproduction to favour infected females, which transmit Wolbachia to their progeny and increase the presence of Wolbachia in the population. The reproductive alterations caused by Wolbachia include feminization, parthenogenesis, male killing and cytoplasmic incompatibility (CI), among which CI is the most common. CI leads to embryonic lethality when Wolbachia-infected males mate with uninfected females or those infected with an incompatible strain. This lethality can be rescued if females are infected with a compatible strain. Although CI was described in the 1960s and its connection to Wolbachia was made in the 1970s, the genes responsible for CI, called CI factors, were not identified until recently. Since then, significant progress has been made in understanding the molecular mechanism of CI using a combination of genetic, phylogenetic, biochemical and structural approaches. The detailed molecular mechanisms behind this fascinating endosymbiotic bacteria-induced phenotype have begun to emerge. Here, we summarize recent progress in understanding the molecular mechanism of CI, especially focusing on the recently solved CI factor structures and discussing what these new structures brought in terms of CI mechanism.

cifB-transcript levels largely explain cytoplasmic incompatibility variation across divergent Wolbachia

Divergent hosts often associate with intracellular microbes that influence their fitness. Maternally transmitted Wolbachia bacteria are the most common of these endosymbionts, due largely to cytoplasmic incompatibility (CI) that kills uninfected embryos fertilized by Wolbachia-infected males. Closely related infections in females rescue CI, providing a relative fitness advantage that drives Wolbachia to high frequencies. One prophage-associated gene (cifA) governs rescue, and two contribute to CI (cifA and cifB), but CI strength ranges from very strong to very weak for unknown reasons. Here, we investigate CI-strength variation and its mechanistic underpinnings in a phylogenetic context across 20 million years (MY) of Wolbachia evolution in Drosophila hosts diverged up to 50 MY. These Wolbachia encode diverse Cif proteins (100% to 7.4% pairwise similarity), and AlphaFold structural analyses suggest that CifB sequence similarities do not predict structural similarities. We demonstrate that cifB-transcript levels in testes explain CI strength across all but two focal systems. Despite phylogenetic discordance among cifs and the bulk of the Wolbachia genome, closely related Wolbachia tend to cause similar CI strengths and transcribe cifB at similar levels. This indicates that other non-cif regions of the Wolbachia genome modulate cif-transcript levels. CI strength also increases with the length of the host's larval life stage, presumably due to prolonged cif action. Our findings reveal that cifB-transcript levels largely explain CI strength, while highlighting other covariates. Elucidating CI's mechanism contributes to our understanding of Wolbachia spread in natural systems and to improving the efficacy of CI-based biocontrol of arboviruses and agricultural pests globally.

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.

Modifying mosquitoes to suppress disease transmission: Is the long wait over?

For more than 50 years it has been a dream of medical entomologists and public health workers to control diseases like malaria and dengue fever by modifying, through genetics and other methods, the arthropods that transmit them to humans. A brief synopsis of the history of these efforts as applied to mosquitoes is presented; none proved to be effective in reducing disease prevalence. Only in the last few years have novel approaches been developed or proposed that indicate the long wait may be over. Three recent developments are particularly promising: CRISPR-Cas9 driven genetic modification, shifting naturally occurring allele frequencies, and microbe-based modifications. The last is the furthest along in implementation. Dengue fever incidence has been reduced between 40% and 96% in 4 different regions of the world where Wolbachia-infected Aedes aegypti have been established in the field. It is not yet clear how sustainable such control programs will prove to be, but there is good reason for optimism. In light of this, the time is ripe for reinvigorated research on vectors, especially genetics. Vector-borne diseases primarily affect under-developed countries and thus have not received the attention they deserve from wealthier countries with well-developed and funded biomedical research establishments.

Restriction of Wolbachia Bacteria in Early Embryogenesis of Neotropical Drosophila Species via Endoplasmic Reticulum-Mediated Autophagy

Wolbachia are maternally transmitted intracellular bacteria that are not only restricted to the reproductive organs but also found in various somatic tissues of their native hosts. The abundance of the endosymbiont in the soma, usually a dead end for vertically transmitted bacteria, causes a multitude of effects on life history traits of their hosts, which are still not well understood. Thus, deciphering the host-symbiont interactions on a cellular level throughout a host's life cycle is of great importance to understand their homeostatic nature, persistence, and spreading success. Using fluorescent and transmission electron microscopy, we conducted a comprehensive analysis of Wolbachia tropism in soma and germ line of six Drosophila species at the intracellular level during host development. Our data uncovered diagnostic patterns of infections to embryonic primordial germ cells and to particular cells of the soma in three different neotropical Drosophila species that have apparently evolved independently. We further found that restricted patterns of Wolbachia tropism are determined in early embryogenesis via selective autophagy, and their spatially restricted infection patterns are preserved in adult flies. We observed tight interactions of Wolbachia with membranes of the endoplasmic reticulum, which might play a scaffolding role for autophagosome formation and subsequent elimination of the endosymbiont. Finally, by analyzing D. simulans lines transinfected with nonnative Wolbachia, we uncovered that the host genetic background regulates tissue tropism of infection. Our data demonstrate a novel and peculiar mechanism to limit and spatially restrict bacterial infection in the soma during a very early stage of host development. IMPORTANCE All organisms are living in close and intimate interactions with microbes that cause conflicts but also cooperation between both unequal genetic partners due to their different innate interests of primarily enhancing their own fitness. However, stable symbioses often result in homeostatic interaction, named mutualism, by balancing costs and benefits, where both partners profit. Mechanisms that have evolved to balance and stably maintain homeostasis in mutualistic relationships are still quite understudied; one strategy is to "domesticate" potentially beneficial symbionts by actively controlling their replication rate below a critical and, hence, costly threshold, and/or to spatially and temporally restrict their localization in the host organism, which, in the latter case, in its most extreme form, is the formation of a specialized housing organ for the microbe (bacteriome). However, questions remain: how do these mutualistic associations become established in their first place, and what are the mechanisms for symbiont control and restriction in their early stages? Here, we have uncovered an unprecedented symbiont control mechanism in neotropical Drosophila species during early embryogenesis. The fruit fly evolved selective autophagy to restrict and control the proliferation of its intracellular endosymbiont Wolbachia in a defined subset of the stem cells as soon as the host's zygotic genome is activated.

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.

Predicting distributions of Wolbachia strains through host ecological contact—Who's manipulating whom?

Reproductive isolation in response to divergent selection is often mediated via third‐party interactions. Under these conditions, speciation is inextricably linked to ecological context. We present a novel framework for understanding arthropod speciation as mediated by Wolbachia, a microbial endosymbiont capable of causing host cytoplasmic incompatibility (CI). We predict that sympatric host sister‐species harbor paraphyletic Wolbachia strains that provide CI, while well‐defined congeners in ecological contact and recently diverged noninteracting congeners are uninfected due to Wolbachia redundancy. We argue that Wolbachia provides an adaptive advantage when coupled with reduced hybrid fitness, facilitating assortative mating between co‐occurring divergent phenotypes—the contact contingency hypothesis. To test this, we applied a predictive algorithm to empirical pollinating fig wasp data, achieving up to 91.60% accuracy. We further postulate that observed temporal decay of Wolbachia incidence results from adaptive host purging—adaptive decay hypothesis—but implementation failed to predict systematic patterns. We then account for post‐zygotic offspring mortality during CI mating, modeling fitness clines across developmental resources—the fecundity trade‐off hypothesis. This model regularly favored CI despite fecundity losses. We demonstrate that a rules‐based algorithm accurately predicts Wolbachia infection status. This has implications among other systems where closely related sympatric species encounter adaptive disadvantage through hybridization.

Crystal Structures of Wolbachia CidA and CidB Reveal Determinants of Bacteria-induced Cytoplasmic Incompatibility and Rescue

Cytoplasmic incompatibility (CI) results when Wolbachia bacteria-infected male insects mate with uninfected females, leading to embryonic lethality. “Rescue” of viability occurs if the female harbors the same Wolbachia strain. CI is caused by linked pairs of Wolbachia genes called CI factors (CifA and CifB). The co-evolution of CifA-CifB pairs may account in part for the incompatibility patterns documented in insects infected with different Wolbachia strains, but the molecular mechanisms remain elusive. Here, we use X-ray crystallography and AlphaFold to analyze the CI factors from Wolbachia strain wMel called CidA^wMel and CidB^wMel. Substituting CidA^wMel interface residues with those from CidA^wPip (from strain wPip) enables the mutant protein to bind CidB^wPip and rescue CidB^wPip-induced yeast growth defects, supporting the importance of CifA-CifB interaction in CI rescue. Sequence divergence in CidA^wPip and CidB^wPip proteins affects their pairwise interactions, which may help explain the complex incompatibility patterns of mosquitoes infected with different wPip strains.

Wolbachia induced cytoplasmic incompatibility (CI) is caused by linked pairs of genes named cifA and cifB. Here, authors show that the residues at interfaces of the CidA-CidB complex is crucial for their binding and contribute to the diversity of CI.

The CinB Nuclease from wNo Wolbachia Is Sufficient for Induction of Cytoplasmic Incompatibility in Drosophila

Wolbachia is an obligate intracellular bacterium that can alter reproduction of its arthropod hosts, often through a mechanism called cytoplasmic incompatibility (CI). In CI, uninfected females fertilized by infected males yield few offspring, but if both are similarly infected, normal embryo viability results (called "rescue"). CI factors (Cifs) responsible for CI are pairs of proteins encoded by linked genes. The downstream gene in each pair encodes either a deubiquitylase (CidB) or a nuclease (CinB). The upstream gene products, CidA and CinA, bind their cognate enzymes with high specificity. Expression of CidB or CinB in yeast inhibits growth, but growth is rescued by expression of the cognate CifA protein. By contrast, transgenic Drosophila male germ line expression of both cifA and cifB was reported to be necessary to induce CI-like embryonic arrest; cifA expression alone in females is sufficient for rescue. This pattern, seen with genes from several Wolbachia strains, has been called the "2-by-1" model. Here, we show that male germ line expression of the cinB gene alone, from a distinct clade of cif genes from wNo Wolbachia, is sufficient to induce nearly complete loss of embryo viability. This male sterility is fully rescued by cognate cinA^wNo expression in the female germ line. The proteins behave similarly in yeast. CinB^wNo toxicity depends on its nuclease active site. These results demonstrate that highly divergent CinB nucleases can induce CI, that rescue by cognate CifA factors is a general feature of Wolbachia CI systems, and that CifA is not strictly required in males for CI induction. IMPORTANCE Wolbachia bacteria live within the cells of many insects. Like mitochondria, they are only inherited from females. Wolbachia often increases the number of infected females to promote spread of infection using a type of male sterility called cytoplasmic incompatibility (CI): when uninfected females mate with infected males, most embryos die; if both are similarly infected, embryos develop normally, giving infected females an advantage in producing offspring. CI is being used against disease-carrying mosquitoes and agricultural pests. Wolbachia proteins called CifA and CifB, which bind one another, cause CI, but how they work has been unclear. Here, we show that a CifB protein singly produced in fruit fly males causes sterility in crosses to normal females, but this is rescued if the females produce the CifA partner. These findings clarify a broad range of observations on CI and will allow more rational approaches to using it for insect control.

From Wolbachia genomics to phenotype: molecular models of cytoplasmic incompatibility must account for the multiplicity of compatibility types

Wolbachia endosymbionts commonly induce cytoplasmic incompatibility, making infected males' sperm lethal to the embryos unless these are rescued by the same bacterium, inherited from their mother. Causal genes were recently identified but two families of mechanistic models are still opposed. In the toxin-antidote model, interaction between the toxin and the antidote is required for rescuing the embryos. In host modification models, a host factor is misregulated in sperm and rescue occurs through compensation or withdrawal of this modification. While these models have been thoroughly discussed, the multiplicity of compatibility types, that is, the existence of many mutually incompatible strains, as seen in Culex mosquitoes, has not received sufficient attention. To explain such a fact, host modification models must posit that the same embryonic defects can be induced and rescued through a large variety of host targets. Conversely, the toxin-antidote model simply accommodates this pattern in a lock-key fashion, through variations in the toxin-antidote interaction sites.

Molecular detection and characterization of the endosymbiont Wolbachia in the European hedgehog flea, Archaeopsylla erinacei

Wolbachia, the endosymbiont of arthropods and onchocercid nematodes is present in many medically important insect species, being also considered for the indirect control of parasitic ones. Archaeopsylla erinacei is a flea species infesting hedgehogs acting as vector of Rickettsia felis, Bartonella henselae, and Rickettsia helvetica, thus having public health relevance. The Wolbachia surface protein (wsp) and 16S rRNA genes were used to determine the presence, prevalence and molecular typing of Wolbachia in this flea species collected in two regions of southern Italy. Of the 45 fleas tested (n = 16 males, 35.6%; n = 29 females, 64.4%), 43 (95.6%; 95% CI: 84.8-99.2) scored positive for Wolbachia, of which 15 (33.3%) and 28 (62.2%) were males and females, respectively. The sex-wise prevalence of this endosymbiont was almost equal in both sexes (males 93.8%; 95% CI: 69.5-99.7; females 96.7%; 95% CI: 83.1-99.8). Single locus sequence analysis (SLST) of Wolbachia revealed two sequence types for 16S rRNA gene, named as wAr_15227 and wAr_15234, which came from two different areas, equally distributed in male and female fleas, whilst only one sequence type was identified for wsp gene. The phylogenetic analysis placed the two 16S rRNA sequence types in paraphyletic clades belonging to the supergroup A and B, respectively. Whilst, the tree of wsp gene clustered the corresponding sequence in the same clade including those of Wolbachia supergroup A. In MLST analyses, both Wolbachia sequence types clustered in a monophyletic clade with Drosophila nikananu (wNik) and Drosophila sturtevanti (wStv) from supergroup A. ClonalFrame analysis revealed a recombination event in the wAr_15234 strain which came from Apulia region. Scientific knowledge of the presence/prevalence of Wolbachia among medically important fleas, may contribute to develop an alternative biological method for the vector control.

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

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

Structural and mechanistic insights into the complexes formed by Wolbachia cytoplasmic incompatibility factors

Wolbachia bacteria, inherited through the female germ line, infect a large fraction of arthropod species. Many Wolbachia strains manipulate host reproduction, most commonly through cytoplasmic incompatibility (CI). CI, a conditional male sterility, results when Wolbachia-infected male insects mate with uninfected females; viability is restored if the female is similarly infected (called "rescue"). CI is used to help control mosquito-borne viruses such as dengue and Zika, but its mechanisms remain unknown. The coexpressed CI factors CifA and CifB form stable complexes in vitro, but the timing and function of this interaction in the insect are unresolved. CifA expression in the female germ line is sufficient for rescue. We report high-resolution structures of a CI-factor complex, CinA-CinB, which utilizes a unique binding mode between the CinA rescue factor and the CinB nuclease; the structures were validated by biochemical and yeast growth analyses. Importantly, transgenic expression in Drosophila of a nonbinding CinA mutant, designed based on the CinA-CinB structure, suggests CinA expressed in females must bind CinB imported by sperm in order to rescue embryonic viability. Binding between cognate factors is conserved in an enzymatically distinct CI system, CidA-CidB, suggesting universal features in Wolbachia CI induction and rescue.

A single mutation weakens symbiont-induced reproductive manipulation through reductions in deubiquitylation efficiency

Animals interact with microbes that affect their performance and fitness, including endosymbionts that reside inside their cells. Maternally transmitted Wolbachia bacteria are the most common known endosymbionts, in large part because of their manipulation of host reproduction. For example, many Wolbachia cause cytoplasmic incompatibility (CI) that reduces host embryonic viability when Wolbachia-modified sperm fertilize uninfected eggs. Operons termed cifs control CI, and a single factor (cifA) rescues it, providing Wolbachia-infected females a fitness advantage. Despite CI's prevalence in nature, theory indicates that natural selection does not act to maintain CI, which varies widely in strength. Here, we investigate the genetic and functional basis of CI-strength variation observed among sister Wolbachia that infect Drosophila melanogaster subgroup hosts. We cloned, Sanger sequenced, and expressed cif repertoires from weak CI-causing wYak in Drosophila yakuba, revealing mutations suspected to weaken CI relative to model wMel in D. melanogaster A single valine-to-leucine mutation within the deubiquitylating (DUB) domain of the wYak cifB homolog (cidB) ablates a CI-like phenotype in yeast. The same mutation reduces both DUB efficiency in vitro and transgenic CI strength in the fly, each by about twofold. Our results map hypomorphic transgenic CI to reduced DUB activity and indicate that deubiquitylation is central to CI induction in cid systems. We also characterize effects of other genetic variation distinguishing wMel-like cifs Importantly, CI strength determines Wolbachia prevalence in natural systems and directly influences the efficacy of Wolbachia biocontrol strategies in transinfected mosquito systems. These approaches rely on strong CI to reduce human disease.

Genetic innovations in animal-microbe symbioses

Animal hosts have initiated myriad symbiotic associations with microorganisms and often have maintained these symbioses for millions of years, spanning drastic changes in ecological conditions and lifestyles. The establishment and persistence of these relationships require genetic innovations on the parts of both symbionts and hosts. The nature of symbiont innovations depends on their genetic population structure, categorized here as open, closed or mixed. These categories reflect modes of inter-host transmission that result in distinct genomic features, or genomic syndromes, in symbionts. Although less studied, hosts also innovate in order to preserve and control symbiotic partnerships. New capabilities to sequence host-associated microbial communities and to experimentally manipulate both hosts and symbionts are providing unprecedented insights into how genetic innovations arise under different symbiont population structures and how these innovations function to support symbiotic relationships.