Human activities and densities shape insecticide resistance distribution and dynamics in the virus-vector Culex pipiens mosquitoes from Morocco.

BACKGROUND: Mosquitoes of the Culex pipiens complex are widely distributed vectors for several arboviruses affecting humans. Consequently, their populations have long been controlled using insecticides, in response to which different resistance mechanisms have been selected. Moreover, their ecological preferences and broad adaptability allow C. pipiens mosquitoes to breed in highly polluted water bodies where they are exposed to many residuals from anthropogenic activities. It has been observed for several mosquito species that anthropization (in particular urbanization and agricultural lands) can lead to increased exposure to insecticides and thus to increased resistance. The main objective of the present study was to investigate whether and how urbanization and/or agricultural lands had a similar impact on C. pipiens resistance to insecticides in Morocco. METHODS: Breeding sites were sampled along several transects in four regions around major Moroccan cities, following gradients of decreasing anthropization. The imprint of anthropogenic activities was evaluated around each site as the percentage of areas classified in three categories: urban, agricultural and natural. We then assessed the frequencies of four known resistance alleles in these samples and followed their dynamics in five urban breeding sites over 4 years. RESULTS: The distribution of resistance alleles revealed a strong impact of anthropization, in both agricultural and urbanized lands, although different between resistance mutations and between Moroccan regions; we did not find any clear trend in the dynamics of these resistance alleles during the survey. CONCLUSIONS: Our study provides further evidence for the role of anthropic activities in the selection and maintenance of mutations selected for resistance to insecticides in mosquitoes. The consequences are worrying as this could decrease vector control capacities and thus result in epizootic and epidemic outbreaks. Consequently, concerted and integrated disease control strategies must be designed that include better management regarding the consequences of our activities.

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.

Despite structural identity, ace-1 heterogenous duplication resistance alleles are quite diverse in Anopheles mosquitoes.

Anopheles gambiae s.l. has been the target of intense insecticide treatment since the mid-20th century to try and control malaria. A substitution in the ace-1 locus has been rapidly selected for, allowing resistance to organophosphate and carbamate insecticides. Since then, two types of duplication of the ace-1 locus have been found in An. gambiae s.l. populations: homogeneous duplications that are composed of several resistance copies, or heterogeneous duplications that contain both resistance and susceptible copies. The substitution induces a trade-off between resistance in the presence of insecticides and disadvantages in their absence: the heterogeneous duplications allow the fixation of the intermediate heterozygote phenotype. So far, a single heterogeneous duplication has been described in An. gambiae s.l. populations (in contrast with the multiple duplicated alleles found in Culex pipiens mosquitoes). We used a new approach, combining long and short-read sequencing with Sanger sequencing to precisely identify and describe at least nine different heterogeneous duplications, in two populations of An. gambiae s.l. We show that these alleles share the same structure as the previously identified heterogeneous and homogeneous duplications, namely 203-kb tandem amplifications with conserved breakpoints. Our study sheds new light on the origin and maintenance of these alleles in An. gambiae s.l. populations, and their role in mosquito adaptation.

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.

Functional analysis of Wolbachia Cid effectors unravels cooperative interactions to target host chromatin during replication.

Wolbachia are common bacteria among terrestrial arthropods. These endosymbionts transmitted through the female germline manipulate their host reproduction through several mechanisms whose most prevalent form called Cytoplasmic Incompatibility -CI- is a conditional sterility syndrome eventually favoring the infected progeny. Upon fertilization, the sperm derived from an infected male is only compatible with an egg harboring a compatible Wolbachia strain, this sperm leading otherwise to embryonic death. The Wolbachia Cif factors CidA and CidB responsible for CI and its neutralization function as a Toxin-Antitoxin system in the mosquito host Culex pipiens. However, the mechanism of CidB toxicity and its neutralization by the CidA antitoxin remain unexplored. Using transfected insect cell lines to perform a structure-function analysis of these effectors, we show that both CidA and CidB are chromatin interactors and CidA anchors CidB to the chromatin in a cell-cycle dependent-manner. In absence of CidA, the CidB toxin localizes to its own chromatin microenvironment and acts by preventing S-phase completion, independently of its deubiquitylase -DUB- domain. Experiments with transgenic Drosophila show that CidB DUB domain is required together with CidA during spermatogenesis to stabilize the CidA-CidB complex. Our study defines CidB functional regions and paves the way to elucidate the mechanism of its toxicity.

Evolutionary trade-offs associated with copy number variations in resistance alleles in Culex pipiens mosquitoes.

Organophosphate and carbamate insecticides have largely been used worldwide to control mosquito populations. As a response, the same amino acid substitution in the ace-1 gene (G119S), conferring resistance to both insecticides, has been selected independently in many mosquito species. In Anopheles gambiae, it has recently been shown that the G119S mutation is actually part of homogeneous duplications that associate multiple resistance copies of the ace-1 gene. In this study, we showed that duplications of resistance copies of the ace-1 gene also exist in the Culex pipiens species complex. The number of copies is variable, and different numbers of copies are associated with different phenotypic trade-offs: we used a combination of bioassays and competition in population cages to show that having more resistance copies conferred higher resistance levels, but was also associated with higher selective disadvantage (or cost) in the absence of insecticide. These results further show the versatility of the genetic architecture of resistance to organophosphate and carbamate insecticides around the ace-1 locus and its role in fine-tuned adaptation to insecticide treatment variations.

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.

Spatio-Temporal Dynamics of a Dieldrin Resistance Gene in Aedes albopictus and Culex quinquefasciatus Populations From Reunion Island.

The control of mosquito populations using insecticides is increasingly threatened by the spread of resistance mechanisms. Dieldrin resistance, conferred by point mutations in the Rdl gene encoding the γ-aminobutyric acid receptor, has been reported at high prevalence in mosquito populations in response to selective pressures. In this study, we monitored spatio-temporal dynamics of the resistance-conferring RdlR allele in Aedes (Stegomyia) albopictus (Skuse, 1895) and Culex (Culex) quinquefasciatus (Say, 1823) populations from Reunion Island. Specimens of both mosquito species were sampled over a 12-month period in three cities and in sites located at lower (<61 m) and higher (between 503 and 564 m) altitudes. Mosquitoes were genotyped using a molecular test detecting the alanine to serine substitution (A302S) in the Rdl gene. Overall, the RdlR frequencies were higher in Cx. quinquefasciatus than Ae. albopictus. For both mosquito species, the RdlR frequencies were significantly influenced by location and altitude with higher RdlR frequencies in the most urbanized areas and at lower altitudes. This study highlights environmental factors that influence the dynamics of insecticide resistance genes, which is critical for the management of insecticide resistance and the implementation of alternative and efficient vector control strategies.

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.

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.

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

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.IMPORTANCECulex 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. Recent studies revealed that (i) two genes, cidA and cidB, are central in Wolbachia-CI mechanisms, and (ii) compatibility versus incompatibility between mosquito lines depends on the wPip phylogenetic groups at play. Here, we studied CI variations in relation to wPip groups and cid genes divergence. We showed, as expected, that the crosses involving wPip from the same group were compatible. In contrast, 78% of the crosses involving different wPip groups were partially or fully incompatible. In such crosses, we reported defects during the first zygotic division, a hallmark of CI. We showed that CI was more severe and frequent in crosses involving wPip-IV strains exhibiting cid variants, which markedly diverge from those of other wPip groups.

Distribution of insecticide resistance and molecular mechanisms involved in the West Nile vector Culex pipiens in Morocco.

BACKGROUND: Mosquitoes of the Culex pipiens complex are the vectors of several arboviruses and are thus subjected to insecticide control worldwide. However, overuse of insecticides selects for resistance. While assessing the resistance status of the vectors is required for effective and sustainable disease control, resistance has so far only been sparsely studied in Morocco. In this study, we establish a first countrywide assessment of the levels of resistance to various insecticides and the potential responsible mechanisms involved. Cx. pipiens larvae were collected from natural populations of five regions of Morocco, and their taxonomic status was determined (molecular forms). The level of their susceptibility to insecticides was assessed by single-diagnostic-dose bioassays. Molecular identification of known resistance alleles was investigated to determine the frequency of target-site mutations. RESULTS: This study confirms that Moroccan populations are an interbreeding mix of pipiens and molestus forms, with large gene flow for the resistance alleles. We also found that Cx. pipiens mosquitoes are resistant to all insecticide families, all over Morocco: resistance is high for insecticides used in mosquito control, but also present for other pesticides. Resistance alleles are similarly more frequent for mosquito control insecticides. However, their distribution is heterogeneous in the five regions, with significant genetic differentiation between populations, revealing the crucial role of local insecticide treatment practices. CONCLUSION: This study provides reference countrywide data that highlight the need for further research to refine the distribution of resistance in Morocco and to understand the role of agriculture/urban residuals in its spread. © 2020 Society of Chemical Industry.

Wolbachia modulates prevalence and viral load of Culex pipiens densoviruses in natural populations.

The inadequacy of standard mosquito control strategies calls for ecologically safe novel approaches, for example the use of biological agents such as the endosymbiotic α-proteobacteria Wolbachia or insect-specific viruses (ISVs). Understanding the ecological interactions between these "biocontrol endosymbionts" is thus a fundamental step. Wolbachia are transmitted vertically from mother to offspring and modify their hosts' phenotypes, including reproduction (e.g., cytoplasmic incompatibility) and survival (e.g., viral interference). In nature, Culex pipiens (sensu lato) mosquitoes are always found infected with genetically diverse Wolbachia called wPip that belong to five phylogenetic groups. In recent years, ISVs have also been discovered in these mosquito species, although their interactions with Wolbachia in nature are unknown. Here, we studied the interactions between a widely prevalent ISV, the Culex pipiens densovirus (CpDV, Densovirinae), and Wolbachia in northern Tunisian C. pipiens populations. We showed an influence of different Wolbachia groups on CpDV prevalence and a general positive correlation between Wolbachia and CpDV loads. By investigating the putative relationship between CpDV diversification and wPip groups in the different sites, we detected a signal linked to wPip groups in CpDV phylogeny in sites where all larvae were infected by the same wPip group. However, no such signal was detected where the wPip groups coexisted, suggesting CpDV horizontal transfer between hosts. Overall, our results provide good evidence for an ecological influence of Wolbachia on an ISV, CpDV, in natural populations and highlight the importance of integrating Wolbachia in our understanding of ISV ecology in nature.

Identification and detection of a novel point mutation in the Chitin Synthase gene of Culex pipiens associated with diflubenzuron resistance.

BACKGROUND: Diflubenzuron (DFB) is one of the most used insecticides in mosquito larval control including that of Culex pipiens, the proven vector of the recent West Nile Virus epidemics in Europe. Two mutations (I1043L and I1043M) in the chitin synthase (CHS) putative binding site of DFB have been previously reported in Cx. pipiens from Italy and associated with high levels of resistance against this larvicide. METHODOLOGY/PRINCIPAL FINDINGS: Here we report the identification of a third mutation at the same I1043 position of the CHS gene resulting in the substitution of Isoleucine to Phenylalanine (I1043F). This mutation has also been found in agricultural pests and has been functionally validated with genome editing in Drosophila, showing to confer striking levels (>15,000 fold) of DFB resistance. The frequency of the I1043F mutation was found to be substantially higher in Cx. pipiens mosquitoes surviving DFB doses largely exceeding the recommended field dose, raising concerns about the future efficient use of this insecticide. We monitored the presence and frequency of DFB mutations in Cx. pipiens mosquitoes from several Mediterranean countries, including Italy, France, Greece, Portugal and Israel. Among the Cx. pipiens populations collected in Northern Italy all but one had at least one of the three DFB mutations at allele frequencies reaching 93.3% for the I1043M, 64.8% for the I1043L and 10% for the I1043F. The newly reported I1043F mutation was also identified in two heterozygote individuals from France (4.2% allelic frequency). In contrast to Italy and France, no DFB resistant mutations were identified in the Cx. pipiens mosquitoes sampled from Greece, Portugal and Israel. CONCLUSIONS/SIGNIFICANCE: The findings of our study are of major concern for mosquito control programs in Europe, that rely on the use of a limited number of available larvicides, and highlight the necessity for the development of appropriate Insecticide Resistance Management (IRM) programs, to ensure the sustainable use of DFB.

Dynamic of resistance alleles of two major insecticide targets in Anopheles gambiae (s.l.) populations from Benin, West Africa.

BACKGROUND: Insecticide resistance is a growing concern for malaria control and vector control effectiveness relies on assessing it distribution and understanding its evolution. METHODS: We assessed resistance levels and the frequencies of two major target-site mutations, L1014F-VGSC and G119S-ace-1, conferring resistance to pyrethroids (PYRs) and carbamates/organophosphates (CXs/OPs) insecticides. These data were compared to those acquired between 2006 and 2010 to follow resistance evolutionary trends over ten years. RESULTS: We report the results of a 3-year survey (2013-2015) of insecticide resistance in 13 localities across the whole country of Benin. Permethrin (PYR) resistance was found in all populations tested, L1014F-VGSC being almost fixed everywhere, while bendiocarb resistance was limited to a few localities, G119S-ace-1 remaining rare, with very limited variations during surveyed period. Interestingly, we found no effect of the type of insecticide pressure on the dynamics of these mutations. CONCLUSIONS: These results confirm both the high prevalence of PYR resistance and the potential of CXs/OPs as short- to medium-term alternatives in Benin. They also underline the need for regular resistance monitoring and informed management in their usage, as the G119S-ace-1 mutation is already present in Benin and surrounding countries. Their unwise usage would rapidly lead to its spread, which would jeopardize PYR-resistant Anopheles control.

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

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

Evolution and phylogeography of Culex pipiens densovirus.

Viruses of the Parvoviridae family infect a wide range of animals including vertebrates and invertebrates. So far, our understanding of parvovirus diversity is biased towards medically or economically important viruses mainly infecting vertebrate hosts, while invertebrate infecting parvoviruses-namely densoviruses-have been largely neglected. Here, we investigated the prevalence and the evolution of the only mosquito-infecting ambidensovirus, Culex pipiens densovirus (CpDV), from laboratory mosquito lines and natural populations collected worldwide. CpDV diversity generally grouped in two clades, here named CpDV-1 and -2. The incongruence of the different gene trees for some samples suggested the possibility of recombination events between strains from different clades. We further investigated the role of selection on the evolution of CpDV genome and detected many individual sites under purifying selection both in non-structural and structural genes. However, some sites in structural genes were under diversifying selection, especially during the divergence of CpDV-1 and -2 clades. These substitutions between CpDV-1 and -2 clades were mostly located in the capsid protein encoding region and might cause changes in host specificity or pathogenicity of CpDV strains from the two clades. However, additional functional and experimental studies are necessary to fully understand the protein conformations and the resulting phenotype of these substitutions between clades of CpDV.

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.

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

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

Wolbachia prevalence, diversity, and ability to induce cytoplasmic incompatibility in mosquitoes.

To protect humans and domestic animals from mosquito borne diseases, alternative methods to chemical insecticides have to be found. Pilot studies using the vertically transmitted bacterial endosymbiont Wolbachia were already launched in different parts of the world. Wolbachia can be used either in Incompatible Insect Technique (IIT), to decrease mosquito population, or to decrease the ability of mosquitoes to transmit pathogens. Not all mosquito species are naturally infected with Wolbachia: while in Culex pipiens and Aedes albopictus almost all individuals harbor Wolbachia, putative infections have to be further investigated in Anopheles species and in Aedes aegypti. All Wolbachia-based control methods rely on the ability of Wolbachia to induce cytoplasmic incompatibility (CI) resulting in embryonic death in incompatible crossings. Knowledge on CI diversity in mosquito is required to find the better Wolbachia-mosquito associations to optimize the success of both 'sterile insect' and 'pathogen blocking' Wolbachia-based methods.