Two-By-One model of cytoplasmic incompatibility: Synthetic recapitulation by transgenic expression of cifA and cifB in Drosophila

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.

Cytoplasmic incompatibility factor proteins from Wolbachia prophage are costly to sperm development in Drosophila melanogaster

The symbiosis between arthropods and Wolbachia bacteria is globally widespread, largely due to selfish-drive systems that favour the fitness of symbiont-transmitting females. The most common drive, cytoplasmic incompatibility (CI), is central to arboviral control efforts. In Drosophila melanogaster carrying wMel Wolbachia deployed in mosquito control, two prophage genes in Wolbachia, cifA and cifB, cause CI that results in a paternal-effect lethality of embryos in crosses between Wolbachia-bearing males and aposymbiotic females. While the CI mechanism by which Cif proteins alter sperm development has recently been elucidated in D. melanogaster and Aedes aegypti mosquitoes, the Cifs' extended impact on male reproductive fitness such as sperm morphology and quantity remains unclear. Here, using cytochemical, microscopic and transgenic assays in D. melanogaster, we demonstrate that both CifA and CifB cause a significant portion of defects in elongating spermatids, culminating in malformed mature sperm nuclei. Males expressing Cifs have reduced spermatid bundles and sperm counts, and transgenic expression of Cifs can occasionally result in no mature sperm formation. We reflect on Cifs' varied functional impacts on the Host Modification model of CI as well as host evolution, behaviour and vector control strategies.

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.

Gene Drive and Symbiont Technologies for Control of Mosquito-Borne Diseases

Mosquito-borne diseases, such as dengue and malaria, pose a significant burden to global health. Current control strategies with insecticides are only moderately effective. Scalable solutions are needed to reduce the transmission risk of these diseases. Symbionts and genome engineering-based mosquito control strategies have been proposed to address these problems. Bacterial, fungal, and viral symbionts affect mosquito reproduction, reduce mosquito lifespan, and block pathogen transmission. Field tests of endosymbiont Wolbachia-based methods have yielded promising results, but there are hurdles to overcome due to the large-scale rearing and accurate sex sorting required for Wolbachia-based suppression approaches and the ecological impediments to Wolbachia invasion in replacement approaches. Genome engineering-based methods, in which mosquitoes are genetically altered for the modification or suppression of wild populations, offer an additional approach for control of mosquito-borne diseases. In particular, the use of gene drive alleles that bias inheritance in their favor is a potentially powerful approach. Several drives are frequency dependent, potentially giving them broadly similar population dynamics to Wolbachia. However, public acceptance and the behavior of released drives in natural mosquito populations remain challenges. We summarize the latest developments and discuss the knowledge gaps in both symbiont- and gene drive-based methods.

Post-release monitoring pathway for the deployment of gene drive-modified mosquitoes for malaria control in Africa

BACKGROUND

Gene drive-modified mosquitoes (GDMMs) have been promoted as one of the innovative technologies that may control and eliminate malaria and other mosquito-borne diseases. Several products are in early stages of development, targeting either population suppression or population modification of the mosquito vector. However, there is no direct experience of conducting risk assessment for environmental releases and subsequent policies regarding conditions for post-release. This study was carried out to gain additional insights on the possible post-release concerns that may arise, as they may inform future risk assessment and planning for deployment.

METHODS

This study involved desktop reviews on post release monitoring experiences with previously released biological control products. Stakeholder consultations involving online surveys, and face to face workshop with experts from selected African countries from Eastern, Western, and Southern African regions was then carried out to establish post-release monitoring concerns for GDMMs.

RESULTS

Review of genetic biocontrol technologies showed only limited lessons from post-release monitoring regimes with a focus largely limited to efficacy. For genetically modified organisms general surveillance and case-specific monitoring is expected in some of the regions. A number of post-release monitoring concerns in relation to the protection goals of human and animal health, biodiversity, and water quality were identified.

CONCLUSION

Based on established- protection goals, several post-release monitoring concerns have been identified. Subject to a rigorous risk assessment process for future GDMMs products, the concerns may then be prioritized for post-release monitoring.

Three feminizing Wolbachia strains in a single host species: comparative genomics paves the way for identifying sex reversal factors

Introduction

Endosymbiotic bacteria in the genus Wolbachia have evolved numerous strategies for manipulating host reproduction in order to promote their own transmission. This includes the feminization of males into functional females, a well-studied phenotype in the isopod Armadillidium vulgare. Despite an early description of this phenotype in isopods and the development of an evolutionary model of host sex determination in the presence of Wolbachia, the underlying genetic mechanisms remain elusive.

Methods

Here we present the first complete genomes of the three feminizing Wolbachia (wVulC, wVulP, and wVulM) known to date in A. vulgare. These genomes, belonging to Wolbachia B supergroup, contain a large number of mobile elements such as WO prophages with eukaryotic association modules. Taking advantage of these data and those of another Wolbachia-derived feminizing factor integrated into the host genome (f element), we used a comparative genomics approach to identify putative feminizing factors.

Results

This strategy has enabled us to identify three prophage-associated genes secreted by the Type IV Secretion System: one ankyrin repeat domain-containing protein, one helix-turn-helix transcriptional regulator and one hypothetical protein. In addition, a latrotoxin-related protein, associated with phage relic genes, was shared by all three genomes and the f element.

Conclusion

These putative feminization-inducing proteins shared canonical interaction features with eukaryotic proteins. These results pave the way for further research into the underlying functional interactions.

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.

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.

PSRs: Selfish chromosomes that manipulate reproductive development

B chromosomes are intriguing "selfish" genetic elements, many of which exhibit higher-than-Mendelian transmission. This perspective highlights a group of B chromosomes known as Paternal Sex Ratio chromosomes (PSRs), which are found in several insects with haplo-diploid reproduction. PSRs harshly alter the organism's reproduction to facilitate their own inheritance. A manifestation of this effect is the conversion of female destined individuals into males. Key to this conversion is the mysterious ability of PSRs to cause elimination of the sperm-inherited half of the genome during zygote formation. Here we discuss how PSRs were discovered, what is known about how they alter paternal chromatin dynamics to cause sex conversion, and how PSR-induced genome elimination is different from other forms of programmed genome elimination in different insects. PSRs also stand out because their DNA sequence compositions differ in remarkable ways from their insect's essential chromosomes, a characteristic suggestive of interspecies origins. Broadly, we also highlight poorly understood aspects of PSR dynamics that need to be investigated.

Prophage proteins alter long noncoding RNA and DNA of developing sperm to induce a paternal-effect lethality

The extent to which prophage proteins interact with eukaryotic macromolecules is largely unknown. In this work, we show that cytoplasmic incompatibility factor A (CifA) and B (CifB) proteins, encoded by prophage WO of the endosymbiont Wolbachia, alter long noncoding RNA (lncRNA) and DNA during Drosophila sperm development to establish a paternal-effect embryonic lethality known as cytoplasmic incompatibility (CI). CifA is a ribonuclease (RNase) that depletes a spermatocyte lncRNA important for the histone-to-protamine transition of spermiogenesis. Both CifA and CifB are deoxyribonucleases (DNases) that elevate DNA damage in late spermiogenesis. lncRNA knockdown enhances CI, and mutagenesis links lncRNA depletion and subsequent sperm chromatin integrity changes to embryonic DNA damage and CI. Hence, prophage proteins interact with eukaryotic macromolecules during gametogenesis to create a symbiosis that is fundamental to insect evolution and vector control.

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.

The mechanism of cytoplasmic incompatibility is conserved in Wolbachia-infected Aedes aegypti mosquitoes deployed for arbovirus control

The rising interest and success in deploying inherited microorganisms and cytoplasmic incompatibility (CI) for vector control strategies necessitate an explanation of the CI mechanism. Wolbachia-induced CI manifests in the form of embryonic lethality when sperm from Wolbachia-bearing testes fertilize eggs from uninfected females. Embryos from infected females however survive to sustain the maternally inherited symbiont. Previously in Drosophila melanogaster flies, we demonstrated that CI modifies chromatin integrity in developing sperm to bestow the embryonic lethality. Here, we validate these findings using wMel-transinfected Aedes aegypti mosquitoes released to control vector-borne diseases. Once again, the prophage WO CI proteins, CifA and CifB, target male gametic nuclei to modify chromatin integrity via an aberrant histone-to-protamine transition. Cifs are not detected in the embryo, and thus elicit CI via the nucleoprotein modifications established pre-fertilization. The rescue protein CifA in oogenesis localizes to stem cell, nurse cell, and oocyte nuclei, as well as embryonic DNA during embryogenesis. Discovery of the nuclear targeting Cifs and altered histone-to-protamine transition in both Aedes aegypti mosquitoes and D. melanogaster flies affirm the Host Modification Model of CI is conserved across these host species. The study also newly uncovers the cell biology of Cif proteins in the ovaries, CifA localization in the embryos, and an impaired histone-to-protamine transition during spermiogenesis of any mosquito species. Overall, these sperm modification findings may enable future optimization of CI efficacy in vectors or pests that are refractory to Wolbachia transinfections.

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.

Use of Drosophila Transgenics to Identify Functions for Symbiont Effectors

Wolbachia, one of the most successful and studied insect symbionts, and Drosophila, one of the most understood model insects, can be exploited as complementary tools to unravel mechanisms of insect symbiosis. Although Wolbachia itself cannot be grown axenically as clonal isolates or genetically manipulated by standard methods, its reproductive phenotypes, including cytoplasmic incompatibility (CI), have been elucidated using well-developed molecular tools and precise transgenic manipulations available for Drosophila melanogaster. Current research only scratches the surface of how Drosophila can provide a tool for understanding Wolbachia's evolutionary success and the molecular roles of its genetic elements. Here, we briefly outline basic methodologies inherent to transgenic Drosophila systems that have already contributed significant advances in understanding CI, but may be unfamiliar to those who lack experience in Drosophila genetics. In the future, these approaches will continue providing significant insights into Wolbachia that undoubtedly will be extended to other insect symbionts and their biological capabilities.

Using Baker's Yeast to Determine Functions of Novel Wolbachia (and Other Prokaryotic) Effectors

Yeasts are single-celled eukaryotic organisms classified as fungi, mostly in the phylum Ascomycota. Of about 1500 named species, Saccharomyces cerevisiae, also known as baker's yeast, domesticated by humans in the context of cooking and brewing, is a profound genetic tool for exploring functions of novel effector proteins from Wolbachia and prokaryotes in general. Wolbachia is a Gram-negative alpha-proteobacterium that infects up to ~75% of all insects as an obligate intracellular microbe (Jeyaprakash A, Hoy MA, Insect Mol Biol 9:393-405, 2000). Wolbachia's lifestyle presents unique challenges for researchers. Wolbachia cannot be axenically cultured and has never been genetically manipulated. Furthermore, many Wolbachia genes have no known function or well-annotated orthologs in other genomes. Yet given the effects of Wolbachia on host phenotypes, which have considerable practical applications for pest control, they undoubtedly involve secreted effector proteins that interact with host gene products. Studying these effectors is challenging with Wolbachia's current genetic limitations. However, some of the constraints to working with Wolbachia can be overcome by expressing candidate proteins in S. cerevisiae. This approach capitalizes on yeast's small genome (~6500 genes), typical eukaryotic cellular organization, and the sophisticated suite of genetic tools available for its manipulation in culture. Thus, yeast can serve as a powerful mock eukaryotic host background to study Wolbachia effector function. Specifically, yeast is used for recombinant protein expression, drug discovery, protein localization studies, protein interaction mapping (yeast two-hybrid system), modeling chromosomal evolution, and examining interactions between proteins responsible for complex phenotypes in less tractable prokaryotic systems. As an example, the paired genes responsible for Wolbachia-mediated cytoplasmic incompatibility (CI) encode novel proteins with limited homology to other known proteins, and no obvious function. This article details how S. cerevisiae was used as an initial staging ground to explore the molecular basis of one of Wolbachia's trademark phenotypes (CI).

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.

Making waves: Comparative analysis of gene drive spread characteristics in a continuous space model

With their ability to rapidly increase in frequency, gene drives can be used to modify or suppress target populations after an initial release of drive individuals. Recent advances have revealed many possibilities for different types of drives, and several of these have been realized in experiments. These drives have advantages and disadvantages related to their ease of construction, confinement and capacity to be used for modification or suppression. Though characteristics of these drives have been explored in modelling studies, assessment in continuous space environments has been limited, often focusing on outcomes rather than fundamental properties. Here, we conduct a comparative analysis of many different gene drive types that have the capacity to form a wave of advance in continuous space using individual-based simulations in continuous space. We evaluate the drive wave speed as a function of drive performance and ecological parameters, which reveals substantial differences between drive performance in panmictic versus spatial environments. In particular, we find that suppression drive waves are uniquely vulnerable to fitness costs and undesired CRISPR cleavage activity in embryos by maternal deposition. Some drives, however, retain robust performance even with widely varying efficiency parameters. To gain a better understanding of drive waves, we compare their panmictic performance and find that the rate of wild-type allele removal is correlated with drive wave speed, though this is also affected by other factors. Overall, our results provide a useful resource for understanding the performance of drives in spatially continuous environments, which may be most representative of potential drive deployment in many relevant scenarios.

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.

Ability of a selfish B chromosome to evade genome elimination in the jewel wasp, Nasonia vitripennis

B chromosomes are non-essential, extra chromosomes that can exhibit transmission-enhancing behaviors, including meiotic drive, mitotic drive, and induction of genome elimination, in plants and animals. A fundamental but poorly understood question is what characteristics allow B chromosomes to exhibit these extraordinary behaviors. The jewel wasp, Nasonia vitripennis, harbors a heterochromatic, paternally transmitted B chromosome known as paternal sex ratio (PSR), which causes complete elimination of the sperm-contributed half of the genome during the first mitotic division of fertilized embryos. This genome elimination event may result from specific, previously observed alterations of the paternal chromatin. Due to the haplo-diploid reproduction of the wasp, genome elimination by PSR causes female-destined embryos to develop as haploid males that transmit PSR. PSR does not undergo self-elimination despite its presence with the paternal chromatin until the elimination event. Here we performed fluorescence microscopic analyses aimed at understanding this unexplained property. Our results show that PSR, like the rest of the genome, participates in the histone-to-protamine transition, arguing that PSR does not avoid this transition to escape self-elimination. In addition, PSR partially escapes the chromatin-altering activity of the intracellular bacterium, Wolbachia, demonstrating that this ability to evade chromatin alteration is not limited to PSR's own activity. Finally, we observed that the rDNA locus and other unidentified heterochromatic regions of the wasp's genome also seem to evade chromatin disruption by PSR, suggesting that PSR's genome-eliminating activity does not affect heterochromatin. Thus, PSR may target an aspect of euchromatin to cause genome elimination.

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.

Convergent Aedes and Drosophila CidB interactomes suggest cytoplasmic incompatibility targets are conserved

Wolbachia-mediated cytoplasmic incompatibility (CI) is a conditional embryonic lethality induced when Wolbachia-modified sperm fertilizes an uninfected egg. The Wolbachia proteins, CidA and CidB control CI. CidA is a rescue factor that reverses lethality. CidA binds to CidB. CidB contains a deubiquitinating enzyme and induces CI. Precisely how CidB induces CI and what it targets are unknown. Likewise, how CidA prevents sterilization by CidB is not clear. To identify CidB substrates in mosquitos we conducted pull-down assays using recombinant CidA and CidB mixed with Aedes aegypti lysates to identify the protein interactomes of CidB and the CidB/CidA protein complex. Our data allow us to cross compare CidB interactomes across taxa for Aedes and Drosophila. Our data replicate several convergent interactions, suggesting that CI targets conserved substrates across insects. Our data support a hypothesis that CidA rescues CI by tethering CidB away from its substrates. Specifically, we identify ten convergent candidate substrates including P32 (protamine-histone exchange factor), karyopherin alpha, ubiquitin-conjugating enzyme, and bicoid stabilizing factor. Future analysis on how these candidates contribute to CI will clarify mechanisms.

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.

Harnessing Wolbachia cytoplasmic incompatibility alleles for confined gene drive: A modeling study

Wolbachia are maternally-inherited bacteria, which can spread rapidly in populations by manipulating reproduction. cifA and cifB are genes found in Wolbachia phage that are responsible for cytoplasmic incompatibility, the most common type of Wolbachia reproductive interference. In this phenomenon, no viable offspring are produced when a male with both cifA and cifB (or just cifB in some systems) mates with a female lacking cifA. Utilizing this feature, we propose new types of toxin-antidote gene drives that can be constructed with only these two genes in an insect genome, instead of the whole Wolbachia bacteria. By using both mathematical and simulation models, we found that a drive containing cifA and cifB together creates a confined drive with a moderate to high introduction threshold. When introduced separately, they act as a self-limiting drive. We observed that the performance of these drives is substantially influenced by various ecological parameters and drive characteristics. Extending our models to continuous space, we found that the drive individual release distribution has a critical impact on drive persistence. Our results suggest that these new types of drives based on Wolbachia transgenes are safe and flexible candidates for genetic modification of populations.

Single-cell transcriptome sequencing reveals Wolbachia-mediated modification in early stages of Drosophila spermatogenesis

Wolbachia are the most widely distributed intracellular bacteria, and their most common effect on host phenotype is cytoplasmic incompatibility (CI). A variety of models have been proposed to decipher the molecular mechanism of CI, among which the host modification (HM) model predicts that Wolbachia effectors play an important role in sperm modification. However, owing to the complexity of spermatogenesis and testicular cell-type heterogeneity, whether Wolbachia have different effects on cells at different stages of spermatogenesis or whether these effects are linked with CI remains unknown. Therefore, we used single-cell RNA sequencing to analyse gene expression profiles in adult male Drosophila testes that were infected or uninfected by Wolbachia. We found that Wolbachia significantly affected the proportion of different types of germ cells and affected multiple metabolic pathways in germ cells. Most importantly, Wolbachia had the greatest impact on germline stem cells, resulting in dysregulated expression of genes related to DNA compaction, and Wolbachia infection also influenced the histone-to-protamine transition in the late stage of sperm development. These results support the HM model and suggest that future studies on Wolbachia-induced CI should focus on cells in the early stages of spermatogenesis.

Wolbachia causes cytoplasmic incompatibility but not male-killing in a grain pest beetle

The endosymbiotic Wolbachia is one of the most common intracellular bacteria known in arthropods and nematodes. Its ability for reproductive manipulation can cause unequal inheritance to male and female offspring, allowing the manipulator to spread, but potentially also impact the evolutionary dynamics of infected hosts. Estimated to be present in up to 66% of insect species, little is known about the phenotypic impact of Wolbachia within the order Coleoptera. Here, we describe the reproductive manipulation by the Wolbachia strain wSur harboured by the sawtoothed grain beetle Oryzaephilus surinamensis (Coleoptera, Silvanidae), through a combination of genomics approaches and bioassays. The Wolbachia strain wSur belongs to supergroup B that contains well-described reproductive manipulators of insects and encodes a pair of cytoplasmic incompatibility factor (cif) genes, as well as multiple homologues of the WO-mediated killing (wmk) gene. A phylogenetic comparison with wmk homologues of wMel of Drosophila melanogaster identified 18 wmk copies in wSur, including one that is closely related to the wMel male-killing homologue. However, further analysis of this particular wmk gene revealed an eight-nucleotide deletion leading to a stop-codon and subsequent reading frame shift midsequence, probably rendering it nonfunctional. Concordantly, utilizing a Wolbachia-deprived O. surinamensis population and controlled mating pairs of wSur-infected and noninfected partners, we found no experimental evidence for male-killing. However, a significant ~50% reduction of hatching rates in hybrid crosses of uninfected females with infected males indicates that wSur is causing cytoplasmic incompatibility. Thus, Wolbachia also represents an important determinant of host fitness in Coleoptera.

Transgenic cytoplasmic incompatibility persists across age and temperature variation in Drosophila melanogaster

Environmental stressors can impact the basic biology and applications of host-microbe symbioses. For example, Wolbachia symbiont densities and cytoplasmic incompatibility (CI) levels can decline in response to extreme temperatures and host aging. To investigate whether transgenic expression of CI-causing cif genes overcomes the environmental sensitivity of CI, we exposed transgenic male flies to low and high temperatures as well as aging treatments. Our results indicate that transgenic cif expression induces nearly complete CI regardless of temperature and aging, despite severe weakening of Wolbachia-based wild-type CI. Strong CI levels correlate with higher levels of cif transgene expression in young males. Altogether, our results highlight that transgenic CI persists against common environmental pressures and may be relevant for future control applications involving the cifA and cifB transgenes.

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.

Penetrance during Wolbachia-mediated parthenogenesis of Trichogramma wasps is reduced by continuous oviposition, associated with exhaustion of Wolbachia titers in ovary and offspring eggs

BACKGROUND

Thelytokous Wolbachia-infected Trichogramma wasps are superior to bisexual uninfected wasps regarding biological control programs. However, continuous oviposition weakens the parthenogenesis-inducing (PI) strength of Wolbachia. Whether this reduced PI strength relates to decreases in the titer of Wolbachia in the ovary and offspring eggs of Trichogramma remains unclear. Here, using fluorescence in situ hybridization (FISH) and reverse transcription quantitative polymerase chain reaction (RT-qPCR) methods, we investigated how the penetrance of Wolbachia-mediated parthenogenesis, Wolbachia density, and distributions of two Wolbachia-infected Trichogramma species, T. pretiosum (TP) and T. dendrolimi (TD), were influenced by different host access treatments [newly-emerged virgin females (NE), 7-day-old females without access to host eggs (NAH), and 7-day-old virgin females with access to host eggs (AH)].

RESULTS

Continuous oviposition decreased Wolbachia PI strength and titers in TP and TD. Continuous oviposition in AH decreased Wolbachia titers in abdomen and offspring eggs of TP and TD females, compared with NAH and NE; NAH had a lower thorax Wolbachia titer than NE. The numbers of parasitized host eggs and offspring wasps, and emergence rates of offspring deposited by AH were lower than those of NE and NAH, for either species.

CONCLUSION

Weakened PI strength, driven by continuous oviposition in Trichogramma wasps, is associated with Wolbachia titer exhaustion in ovary and offspring eggs. Wolbachia density is dependent on PI strength in Trichogramma wasps, highlighting the side effects of continuous oviposition regarding thelytokous Wolbachia-infected Trichogramma in biological control programs. © 2022 Society of Chemical Industry.

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.

Infection Dynamics of Cotransmitted Reproductive Symbionts Are Mediated by Sex, Tissue, and Development

One of the most prevalent intracellular infections on earth is with Wolbachia, a bacterium in the Rickettsiales that infects a range of insects, crustaceans, chelicerates, and nematodes. Wolbachia is maternally transmitted to offspring and has profound effects on the reproduction and physiology of its hosts, which can result in reproductive isolation, altered vectorial capacity, mitochondrial sweeps, and even host speciation. Some populations stably harbor multiple Wolbachia strains, which can further contribute to reproductive isolation and altered host physiology. However, almost nothing is known about the requirements for multiple intracellular microbes to be stably maintained across generations while they likely compete for space and resources. Here, we use a coinfection of two Wolbachia strains (“wHa” and “wNo”) in Drosophila simulans to define the infection and transmission dynamics of an evolutionarily stable double infection. We find that a combination of sex, tissue, and host development contributes to the infection dynamics of the two microbes and that these infections exhibit a degree of niche partitioning across host tissues. wHa is present at a significantly higher titer than wNo in most tissues and developmental stages, but wNo is uniquely dominant in ovaries. Unexpectedly, the ratio of wHa to wNo in embryos does not reflect those observed in the ovaries, indicative of strain-specific transmission dynamics. Understanding how Wolbachia strains interact to establish and maintain stable infections has important implications for the development and effective implementation of Wolbachia-based vector biocontrol strategies, as well as more broadly defining how cooperation and conflict shape intracellular communities.

IMPORTANCEWolbachia is a maternally transmitted intracellular bacterium that manipulates the reproduction and physiology of arthropods, resulting in drastic effects on the fitness, evolution, and even speciation of its hosts. Some hosts naturally harbor multiple strains of Wolbachia that are stably transmitted across generations, but almost nothing is known about the factors that limit or promote these coinfections, which can have profound effects on the host’s biology and evolution and are under consideration as an insect-management tool. Here, we define the infection dynamics of a known stably transmitted double infection in Drosophila simulans with an eye toward understanding the patterns of infection that might facilitate compatibility between the two microbes. We find that a combination of sex, tissue, and development all contributes to infection dynamics of the coinfection.

Perplexing dynamics of Wolbachia proteins for cytoplasmic incompatibility

The mechanism of symbiont-induced cytoplasmic incompatibility has been a long-lasting mystery. This Primer explores a new study on Wolbachia’s Cif proteins in PLOS Biology that provides supportive evidence for the “Host-Modification Model,” although the alternative “Toxin-Antidote Model” is still in the running.

The Cif proteins from Wolbachia prophage WO modify sperm genome integrity to establish cytoplasmic incompatibility

Inherited microorganisms can selfishly manipulate host reproduction to drive through populations. In Drosophila melanogaster, germline expression of the native Wolbachia prophage WO proteins CifA and CifB cause cytoplasmic incompatibility (CI) in which embryos from infected males and uninfected females suffer catastrophic mitotic defects and lethality; however, in infected females, CifA expression rescues the embryonic lethality and thus imparts a fitness advantage to the maternally transmitted Wolbachia. Despite widespread relevance to sex determination, evolution, and vector control, the mechanisms underlying when and how CI impairs male reproduction remain unknown and a topic of debate. Here, we use cytochemical, microscopic, and transgenic assays in D. melanogaster to demonstrate that CifA and CifB proteins of wMel localize to nuclear DNA throughout the process of spermatogenesis. Cif proteins cause abnormal histone retention in elongating spermatids and protamine deficiency in mature sperms that travel to the female reproductive tract with Cif proteins. Notably, protamine gene knockouts enhance wild-type CI. In ovaries, CifA localizes to germ cell nuclei and cytoplasm of early-stage egg chambers; however, Cifs are absent in late-stage oocytes and subsequently in fertilized embryos. Finally, CI and rescue are contingent upon a newly annotated CifA bipartite nuclear localization sequence. Together, our results strongly support the Host modification model of CI in which Cifs initially modify the paternal and maternal gametes to bestow CI-defining embryonic lethality and rescue.

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.

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.

Diversity and function of arthropod endosymbiont toxins

Bacterial endosymbionts induce dramatic phenotypes in their arthropod hosts, including cytoplasmic incompatibility, feminization, parthenogenesis, male killing, parasitoid defense, and pathogen blocking. The molecular mechanisms underlying these effects remain largely unknown but recent evidence suggests that protein toxins secreted by the endosymbionts play a role. Here, we describe the diversity and function of endosymbiont proteins with homology to known bacterial toxins. We focus on maternally transmitted endosymbionts belonging to the Wolbachia, Rickettsia, Arsenophonus, Hamiltonella, Spiroplasma, and Cardinium genera because of their ability to induce the above phenotypes. We identify at least 16 distinct toxin families with diverse enzymatic activities, including AMPylases, nucleases, proteases, and glycosyltransferases. Notably, several annotated toxins contain domains with homology to eukaryotic proteins, suggesting that arthropod endosymbionts mimic host biochemistry to manipulate host physiology, similar to bacterial pathogens.

Environmental Temperature, but Not Male Age, Affects Wolbachia and Prophage WO Thereby Modulating Cytoplasmic Incompatibility in the Parasitoid Wasp, Habrobracon Hebetor

Wolbachia is an endosymbiotic bacterium found in many species of arthropods and manipulates its host reproduction. Cytoplasmic incompatibility (CI) is one of the most common manipulations that is induced when an uninfected female mates with a Wolbachia-infected male. The CI factors (cifA and cifB genes) are encoded by phage WO that naturally infects Wolbachia. Here, we questioned whether an environmental factor (temperature) or host factor (male age) affected the strength of the CI phenotype in the ectoparasitoid wasp, Habrobracon hebetor. We found that temperature, but not male age, results in reduced CI penetrance. Consistent with these results, we also found that the expression of the cif CI factors decreased in temperature-exposed males but was consistent across aging male wasps. Similar to studies of other insect systems, cifA showed a higher expression level than cifB, and male hosts showed increased cif expression relative to females. Our results suggest that prophage WO is present in the Wolbachia-infected wasps and expression of cif genes contributes to the induction of CI in this insect. It seems that male aging has no effect on the intensity of CI; however, temperature affects Wolbachia and prophage WO titers as well as expression levels of cif genes, which modulate the CI level.

Cytoplasmic incompatibility in hybrid zones: infection dynamics and resistance evolution

Cytoplasmic incompatibility is an endosymbiont-induced mating incompatibility common in arthropods. Unidirectional cytoplasmic incompatibility impairs crosses between infected males and uninfected females, whereas bidirectional cytoplasmic incompatibility occurs when two host lineages are infected with reciprocally in compatible endosymbionts. Bidirectional cytoplasmic incompatibility is unstable in unstructured populations, but may be stable in hybrid zones. Stable coexistence of incompatible host lineages should generate frequent incompatible crosses. Therefore, hosts are expected to be under selection to resist their endosymbionts. Here, we for mulate a mathematical model of hybrid zones where two bidirectionally incompatible host lineages meet. We expand this model to consider the invasion of a hypothetical resistance allele. To corroborate our mathematical predictions, we test each prediction with stochastic, individual-based simulations. Our models suggest that hybrid zones may sustain stable coinfections of bidirectionally incompatible endosymbiont strains. Over a range of conditions, host are under selection for resistance against cytoplasmic incompatibility. Under asymetric migration, a resistance allele can facilitate infection turnover and subsequently either persist or become lost. The predictions we present may inform our understanding of the cophylogenetic relationship between the endosym biont Wolbachia and its hosts.

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.

wMel Wolbachia genome remains stable after 7 years in Australian Aedes aegypti field populations

Infection of wMel Wolbachia in Aedes aegypti imparts two signature features that enable its application for biocontrol of dengue. First, the susceptibility of mosquitoes to viruses such as dengue and Zika is reduced. Second, a reproductive manipulation is caused that enables wMel introgression into wild-type mosquito populations. The long-term success of this method relies, in part, on evolution of the wMel genome not compromising the critical features that make it an attractive biocontrol tool. This study compared the wMel Wolbachia genome at the time of initial releases and 1-7 years post-release in Cairns, Australia. Our results show the wMel genome remains highly conserved up to 7 years post-release in gene sequence, content, synteny and structure. This work suggests the wMel genome is stable in its new mosquito host and, therefore, provides reassurance on the potential for wMel to deliver long-term public-health impacts.

Genome analyses of four Wolbachia strains and associated mitochondria of Rhagoletis cerasi expose cumulative modularity of cytoplasmic incompatibility factors and cytoplasmic hitchhiking across host populations

Background

The endosymbiont Wolbachia can manipulate arthropod reproduction and invade host populations by inducing cytoplasmic incompatibility (CI). Some host species are coinfected with multiple Wolbachia strains which may have sequentially invaded host populations by expressing different types of modular CI factor (cif) genes. The tephritid fruit fly Rhagoletis cerasi is a model for CI and Wolbachia population dynamics. It is associated with at least four Wolbachia strains in various combinations, with demonstrated (wCer2, wCer4), predicted (wCer1) or unknown (wCer5) CI phenotypes.

Results

We sequenced and assembled the draft genomes of the Wolbachia strains wCer1, wCer4 and wCer5, and compared these with the previously sequenced genome of wCer2 which currently invades R. cerasi populations. We found complete cif gene pairs in all strains: four pairs in wCer2 (three Type I; one Type V), two pairs in wCer1 (both Type I) and wCer4 (one Type I; one Type V), and one pair in wCer5 (Type IV). Wolbachia genome variant analyses across geographically and genetically distant host populations revealed the largest diversity of single nucleotide polymorphisms (SNPs) in wCer5, followed by wCer1 and then wCer2, indicative of their different lengths of host associations. Furthermore, mitogenome analyses of the Wolbachia genome-sequenced individuals in combination with SNP data from six European countries revealed polymorphic mitogenome sites that displayed reduced diversity in individuals infected with wCer2 compared to those without.

Conclusions

Coinfections with Wolbachia are common in arthropods and affect options for Wolbachia-based management strategies of pest and vector species already infected by Wolbachia. Our analyses of Wolbachia genomes of a host naturally coinfected by several strains unravelled signatures of the evolutionary dynamics in both Wolbachia and host mitochondrial genomes as a consequence of repeated invasions. Invasion of already infected populations by new Wolbachia strains requires new sets of functionally different cif genes and thereby may select for a cumulative modularity of cif gene diversity in invading strains. Furthermore, we demonstrated at the mitogenomic scale that repeated CI-driven Wolbachia invasions of hosts result in reduced mitochondrial diversity and hitchhiking effects. Already resident Wolbachia strains may experience similar cytoplasmic hitchhiking effects caused by the invading Wolbachia strain.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07906-6.

The impacts of cytoplasmic incompatibility factor (cifA and cifB) genetic variation on phenotypes

Wolbachia are maternally transmitted, intracellular bacteria that can often selfishly spread through arthropod populations via cytoplasmic incompatibility (CI). CI manifests as embryonic death when males expressing prophage WO genes cifA and cifB mate with uninfected females or females harboring an incompatible Wolbachia strain. Females with a compatible cifA-expressing strain rescue CI. Thus, cif-mediated CI confers a relative fitness advantage to females transmitting Wolbachia. However, whether cif sequence variation underpins incompatibilities between Wolbachia strains and variation in CI penetrance remains unknown. Here, we engineer Drosophila melanogaster to transgenically express cognate and non-cognate cif homologs and assess their CI and rescue capability. Cognate expression revealed that cifA;B native to D. melanogaster causes strong CI, and cognate cifA;B homologs from two other Drosophila-associated Wolbachia cause weak transgenic CI, including the first demonstration of phylogenetic type 2 cifA;B CI. Intriguingly, non-cognate expression of cifA and cifB alleles from different strains revealed that cifA homologs generally contribute to strong transgenic CI and interchangeable rescue despite their evolutionary divergence, and cifB genetic divergence contributes to weak or no transgenic CI. Finally, we find that a type 1 cifA can rescue CI caused by a genetically divergent type 2 cifA;B in a manner consistent with unidirectional incompatibility. By genetically dissecting individual CI functions for type 1 and 2 cifA and cifB, this work illuminates new relationships between cif genotype and CI phenotype. We discuss the relevance of these findings to CI’s genetic basis, phenotypic variation patterns, and mechanism.