Variation in Wolbachia cidB gene, but not cidA, is associated with cytoplasmic incompatibility mod phenotype diversity in Culex pipiens

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.

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.

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).

Molecular Biology of Cytoplasmic Incompatibility Caused by Wolbachia Endosymbionts

Among endosymbiotic bacteria living within eukaryotic cells, Wolbachia is exceptionally widespread, particularly in arthropods. Inherited through the female germline, it has evolved ways to increase the fraction of bacterially infected offspring by inducing parthenogenesis, feminization, male killing, or, most commonly, cytoplasmic incompatibility (CI). In CI, Wolbachia infection of males causes embryonic lethality unless they mate with similarly infected females, creating a relative reproductive advantage for infected females. A set of related Wolbachia bicistronic operons encodes the CI-inducing factors. The downstream gene encodes a deubiquitylase or nuclease and is responsible for CI induction by males, while the upstream product when expressed in females binds its sperm-introduced cognate partner and rescues viability. Both toxin-antidote and host-modification mechanisms have been proposed to explain CI. Interestingly, male killing by either Spiroplasma or Wolbachia endosymbionts involves deubiquitylases as well. Interference with the host ubiquitin system may therefore be a common theme among endosymbiont-mediated reproductive alterations.

Nanopore sequencing of PCR products enables multicopy gene family reconstruction

The importance of gene amplifications in evolution is more and more recognized. Yet, tools to study multi-copy gene families are still scarce, and many such families are overlooked using common sequencing methods. Haplotype reconstruction is even harder for polymorphic multi-copy gene families. Here, we show that all variants (or haplotypes) of a multi-copy gene family present in a single genome, can be obtained using Oxford Nanopore Technologies sequencing of PCR products, followed by steps of mapping, SNP calling and haplotyping. As a proof of concept, we acquired the sequences of highly similar variants of the cidA and cidB genes present in the genome of the Wolbachia wPip, a bacterium infecting Culex pipiens mosquitoes. Our method relies on a wide database of cid genes, previously acquired by cloning and Sanger sequencing. We addressed problems commonly faced when using mapping approaches for multi-copy gene families with highly similar variants. In addition, we confirmed that PCR amplification causes frequent chimeras which have to be carefully considered when working on families of recombinant genes. We tested the robustness of the method using a combination of bioinformatics (read simulations) and molecular biology approaches (sequence acquisitions through cloning and Sanger sequencing, specific PCRs and digital droplet PCR). When different haplotypes present within a single genome cannot be reconstructed from short reads sequencing, this pipeline confers a high throughput acquisition, gives reliable results as well as insights of the relative copy numbers of the different variants.

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.

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.

No association between habitat, autogeny and genetics in Moroccan Culex pipiens populations

Background

Mosquitoes of the Culex pipiens complex are found across the globe and are the focus of many research studies. Among the temperate species C. pipiens sensu stricto (s.s.), two forms are usually described: molestus and pipiens. These two forms are indistinguishable in terms of morphology but show behavioral and physiological differences that may have consequences for their associated epidemiology. The two forms are well defined in the northern part of the species distribution, where autogeny is strictly associated with the molestus form. However, whether the two remain distinct and show the characteristic differences in behavior is less clear in North Africa, at the southern edge of their range.

Methods

The association between autogeny, as determined by ovarian dissection, and molecular forms, based on the CQ11 microsatellite marker, was studied in six Moroccan populations of C. pipiens.

Results

An overall low prevalence of autogeny was found at three of the Moroccan regions studied, although it reached 17.5% in the Agadir population. The prevalence of form-specific CQ11 alleles was quite similar across all populations, with the molestus allele being rarer (approx. 15%), except in the Agadir population where it reached 43.3%. We found significant deficits in heterozygotes at the diagnostic CQ11 locus in three populations, but the three other populations showed no significant departure from panmixia, which is in line with the results of a retrospective analysis of the published data. More importantly, we found no association between the autogeny status and CQ11 genotypes, despite the many females analyzed.

Conclusions

There was limited evidence for two discrete forms in Morocco, where individuals carrying pipiens and molestus alleles breed and mate in the same sites and are equally likely to be capable of autogeny. These observations are discussed in the epidemiological context of Morocco, where C. pipiens is the main vector of several arboviruses.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05469-3.

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.

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.

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.

The mosquito microbiome includes habitat-specific but rare symbionts

Microbial communities are known to influence mosquito lifestyles by modifying essential metabolic and behavioral processes that affect reproduction, development, immunity, digestion, egg survival, and the ability to transmit pathogens. Many studies have used 16S rRNA gene amplicons to characterize mosquito microbiota and investigate factors that influence host-microbiota dynamics. However, a relatively low taxonomic resolution due to clustering methods based on arbitrary threshold and the overall dominance of Wolbachia or Asaia symbionts obscured the investigation of rare members of mosquito microbiota in previous studies. Here, we used high resolution Shannon entropy-based oligotyping approach to analyze the microbiota of Culex pipiens, Culex quinquefasciatus and Aedes individuals from continental Southern France and overseas Guadeloupe as well as from laboratories with or without antibiotics treatment. Our experimental design that resulted in a series of mosquito samples with a gradient of Wolbachia density and relative abundance along with high-resolution analyses of amplicon sequences enabled the recovery of a robust signal from typically less accessible bacterial taxa. Our data confirm species-specific mosquito-bacteria associations with geography as a primary factor that influences bacterial community structure. But interestingly, they also reveal co-occurring symbiotic bacterial variants within single individuals for both Elizabethkingia and Erwinia genera, distinct and specific Asaia and Chryseobacterium in continental and overseas territories, and a putative rare Wolbachia variant. Overall, our study reveals the presence of previously overlooked microdiversity and multiple closely related symbiotic strains within mosquito individuals with a remarkable habitat-specificity.

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.

Parallel Sequencing of Wolbachia wCer2 from Donor and Novel Hosts Reveals Multiple Incompatibility Factors and Genome Stability after Host Transfers

The application of Wolbachia in insect pest and vector control requires the establishment of genotypically stable host associations. The cytoplasmic incompatibility (CI) inducing Wolbachia strain wCer2 naturally occurs in the cherry fruit fly Rhagoletis cerasi as co-infection with other strains and was transferred to other fruit fly species by embryonic microinjections. We obtained wCer2 genome data from its native and three novel hosts, Drosophila simulans, Drosophila melanogaster, and Ceratitis capitata and assessed its genome stability, characteristics, and CI factor (cif) genes. De novo assembly was successful from Wolbachia cell-enriched singly infected D. simulans embryos, with minimal host and other bacterial genome traces. The low yield of Wolbachia sequence reads from total genomic extracts of one multiply infected R. cerasi pupa and one singly infected C. capitata adult limited de novo assemblies but was sufficient for comparative analyses. Across hosts wCer2 was stable in genome synteny and content. Polymorphic nucleotide sites were found in wCer2 of each host; however, only one nucleotide was different between R. cerasi and C. capitata, and none between replicated D. simulans lines. The wCer2 genome is highly similar to wAu (D. simulans), wMel (D. melanogaster), and wRec (Drosophila recens). In contrast to wMel and wRec (each with one cif gene pair) and wAu (without any cif genes), wCer2 has three pairs of Type I cif genes, and one Type V cifB gene without a cifA complement. This may explain previously reported CI patterns of wCer2, including incomplete rescue of its own CI modification in three novel host species.

Cytoplasmic Incompatibility Variations in Relation with Wolbachia cid Genes Divergence in Culex pipiens

Culex pipiens mosquitoes are infected with wPip. These endosymbionts induce a conditional sterility called CI resulting from embryonic deaths, which constitutes a cornerstone for Wolbachia antivectorial methods.

In arthropods, Wolbachia endosymbionts induce conditional sterility, called cytoplasmic incompatibility (CI), resulting from embryonic lethality. CI penetrance (i.e., embryonic death rate) varies depending on host species and Wolbachia strains involved. All Culex pipiens mosquitoes are infected by the endosymbiotic alphaproteobacteria Wolbachia wPip. CI in Culex, characterized as a binary “compatible/incompatible” phenomenon, revealed an unparalleled diversity of patterns linked to the amplification-diversification of cidA and cidB genes. Here, we accurately studied CI penetrance variations in the light of cid genes divergence by generating a C. pipiens compatibility matrix between 11 lines hosting different phylogenetic wPip groups and exhibiting distinct cid gene repertoires. We showed, as expected, that crosses involving wPip from the same group were mostly compatible. In contrast, only 22% of the crosses involving different wPip groups were compatible, while 54% were fully incompatible. For the remaining 24% of the crosses, “intermediate” compatibilities were reported, and a cytological observation of the first zygotic division confirmed the occurrence of “canonical” CI phenotypes in a fraction of the eggs. Backcross experiments demonstrated that intermediate compatibilities were not linked to host genetic background but to the Wolbachia strains involved. This previously unstudied intermediate penetrance CI was more severe and frequent in crosses involving wPip-IV strains exhibiting cid variants markedly divergent from other wPip groups. Our data demonstrate that CI is not always a binary compatible/incompatible phenomenon in C. pipiens but that intermediate compatibilities putatively resulting from partial mismatch due to Cid proteins divergence exist in this species complex.

A Role for Maternal Factors in Suppressing Cytoplasmic Incompatibility

Wolbachia are maternally transmitted bacterial endosymbionts, carried by approximately half of all insect species. Wolbachia prevalence in nature stems from manipulation of host reproduction to favor the success of infected females. The best known reproductive modification induced by Wolbachia is referred to as sperm-egg Cytoplasmic Incompatibility (CI). In CI, the sperm of Wolbachia-infected males cause embryonic lethality, attributed to paternal chromatin segregation defects during early mitotic divisions. Remarkably, the embryos of Wolbachia-infected females “rescue” CI lethality, yielding egg hatch rates equivalent to uninfected female crosses. Several models have been discussed as the basis for Rescue, and functional evidence indicates a major contribution by Wolbachia CI factors. A role for host contributions to Rescue remains largely untested. In this study, we used a chemical feeding approach to test for CI suppression capabilities by Drosophila simulans. We found that uninfected females exhibited significantly higher CI egg hatch rates in response to seven chemical treatments that affect DNA integrity, cell cycle control, and protein turnover. Three of these treatments suppressed CI induced by endogenous wRi Wolbachia, as well as an ectopic wMel Wolbachia infection. The results implicate DNA integrity as a focal aspect of CI suppression for different Wolbachia strains. The framework presented here, applied to diverse CI models, will further enrich our understanding of host reproductive manipulation by insect endosymbionts.

Symbiont-mediated cytoplasmic incompatibility: what have we learned in 50 years?

Cytoplasmic incompatibility (CI) is the most common symbiont-induced reproductive manipulation. Specifically, symbiont-induced sperm modifications cause catastrophic mitotic defects in the fertilized embryo and ensuing lethality in crosses between symbiotic males and either aposymbiotic females or females harboring a different symbiont strain. However, if the female carries the same symbiont strain, then embryos develop properly, thereby imparting a relative fitness benefit to symbiont-transmitting mothers. Thus, CI drives maternally-transmitted bacteria to high frequencies in arthropods worldwide. In the past two decades, CI experienced a boom in interest due to its (i) deployment in worldwide efforts to curb mosquito-borne diseases, (ii) causation by bacteriophage genes, cifA and cifB, that modify sexual reproduction, and (iii) important impacts on arthropod speciation. This review serves as a gateway to experimental, conceptual, and quantitative themes of CI and outlines significant gaps in understanding CI’s mechanism that are ripe for investigation from diverse subdisciplines in the life sciences.

The Biochemistry of Cytoplasmic Incompatibility Caused by Endosymbiotic Bacteria

Many species of arthropods carry maternally inherited bacterial endosymbionts that can influence host sexual reproduction to benefit the bacterium. The most well-known of such reproductive parasites is Wolbachia pipientis. Wolbachia are obligate intracellular α-proteobacteria found in nearly half of all arthropod species. This success has been attributed in part to their ability to manipulate host reproduction to favor infected females. Cytoplasmic incompatibility (CI), a phenomenon wherein Wolbachia infection renders males sterile when they mate with uninfected females, but not infected females (the rescue mating), appears to be the most common. CI provides a reproductive advantage to infected females in the presence of a threshold level of infected males. The molecular mechanisms of CI and other reproductive manipulations, such as male killing, parthenogenesis, and feminization, have remained mysterious for many decades. It had been proposed by Werren more than two decades ago that CI is caused by a Wolbachia-mediated sperm modification and that rescue is achieved by a Wolbachia-encoded rescue factor in the infected egg. In the past few years, new research has highlighted a set of syntenic Wolbachia gene pairs encoding CI-inducing factors (Cifs) as the key players for the induction of CI and its rescue. Within each Cif pair, the protein encoded by the upstream gene is denoted A and the downstream gene B. To date, two types of Cifs have been characterized based on the enzymatic activity identified in the B protein of each protein pair; one type encodes a deubiquitylase (thus named CI-inducing deubiquitylase or cid), and a second type encodes a nuclease (named CI-inducing nuclease or cin). The CidA and CinA proteins bind tightly and specifically to their respective CidB and CinB partners. In transgenic Drosophila melanogaster, the expression of either the Cid or Cin protein pair in the male germline induces CI and the expression of the cognate A protein in females is sufficient for rescue. With the identity of the Wolbachia CI induction and rescue factors now known, research in the field has turned to directed studies on the molecular mechanisms of CI, which we review here.

Wolbachia pipientis Associated With Tephritid Fruit Fly Pests: From Basic Research to Applications

Members of the true fruit flies (family Tephritidae) are among the most serious agricultural pests worldwide, whose control and management demands large and costly international efforts. The need for cost-effective and environmentally friendly integrated pest management (IPM) has led to the development and implementation of autocidal control strategies. These approaches include the widely used sterile insect technique and the incompatible insect technique (IIT). IIT relies on maternally transmitted bacteria (namely Wolbachia) to cause a conditional sterility in crosses between released mass-reared Wolbachia-infected males and wild females, which are either uninfected or infected with a different Wolbachia strain (i.e., cytoplasmic incompatibility; CI). Herein, we review the current state of knowledge on Wolbachia-tephritid interactions including infection prevalence in wild populations, phenotypic consequences, and their impact on life history traits. Numerous pest tephritid species are reported to harbor Wolbachia infections, with a subset exhibiting high prevalence. The phenotypic effects of Wolbachia have been assessed in very few tephritid species, due in part to the difficulty of manipulating Wolbachia infection (removal or transinfection). Based on recent methodological advances (high-throughput DNA sequencing) and breakthroughs concerning the mechanistic basis of CI, we suggest research avenues that could accelerate generation of necessary knowledge for the potential use of Wolbachia-based IIT in area-wide integrated pest management (AW-IPM) strategies for the population control of tephritid pests.

Variation in Wolbachia cidB gene, but not cidA, is associated with cytoplasmic incompatibility mod phenotype diversity in Culex pipiens

Endosymbiotic Wolbachia bacteria are, to date, considered the most widespread symbionts in arthropods and are the cornerstone of major biological control strategies. Such a high prevalence is based on the ability of Wolbachia to manipulate their hosts' reproduction. One manipulation called cytoplasmic incompatibility (CI) is based on the death of the embryos generated by crosses between infected males and uninfected females or between individuals infected with incompatible Wolbachia strains. CI can be seen as a modification‐rescue system (or mod‐resc) in which paternal Wolbachia produce mod factors, inducing embryonic defects, unless the maternal Wolbachia produce compatible resc factors. Transgenic experiments in Drosophila melanogaster and Saccharomyces cerevisiae converged towards a model where the cidB Wolbachia gene is involved in the mod function while cidA is involved in the resc function. However, as cidA expression in Drosophila males was required to observe CI, it has been proposed that cidA could be involved in both resc and mod functions. A recent correlative study in natural Culex pipiens mosquito populations has revealed an association between specific cidA and cidB variations and changes in mod phenotype, also suggesting a role for both these genes in mod diversity. Here, by studying cidA and cidB genomic repertoires of individuals from newly sampled natural C. pipiens populations harbouring wPipIV strains from North Italy, we reinforce the link between cidB variation and mod phenotype variation fostering the involvement of cidB in the mod phenotype diversity. However, no association between any cidA variants or combination of cidA variants and mod phenotype variation was observed. Taken together our results in natural C. pipiens populations do not support the involvement of cidA in mod phenotype variation.