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.

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.

The evolutionary fate of heterogeneous gene duplications: A precarious overdominant equilibrium between environment, sublethality and complementation.

Gene duplications occur at a high rate. Although most appear detrimental, some homogeneous duplications (identical gene copies) can be selected for beneficial increase in produced proteins. Heterogeneous duplications, which combine divergent alleles of a single locus, are seldom studied due to the paucity of empirical data. We investigated their role in an ongoing adaptive process at the ace-1 locus in Culex pipiens mosquitoes. We assessed the worldwide diversity of the ace-1 alleles (single-copy, susceptible S and insecticide-resistant R, and duplicated D that pair one S and one R copy), analysed their phylogeography and measured their fitness to understand their early dynamics using population genetics models. It provides a coherent and comprehensive evolutionary scenario. We show that D alleles are present in most resistant populations and display a higher diversity than R alleles (27 vs. 4). Most appear to result from independent unequal crossing-overs between local single-copy alleles, suggesting a recurrent process. Most duplicated alleles have a limited geographic distribution, probably resulting from their homozygous sublethality (HS phenotype). In addition, heterozygotes carrying different HS D alleles showed complementation, indicating different recessive lethal mutations. Due to mosaic insecticide control practices, balancing selection (overdominance) plays a key role in the early dynamics of heterogeneous duplicated alleles; it also favours a high local polymorphism of HS D alleles in natural populations (overdominance reinforced by complementation). Overall, our study shows that the evolutionary fate of heterogeneous duplications (and their long-term role) depends on finely balanced selective pressures due to the environment and to their genomic structure.

Population structure of Wolbachia and cytoplasmic introgression in a complex of mosquito species.

BACKGROUND: The maternally inherited bacterium Wolbachia often acts as a subtle parasite that manipulates insect reproduction, resulting potentially in reproductive isolation between host populations. Whilst distinct Wolbachia strains are documented in a group of evolutionarily closely related mosquitoes known as the Culex pipiens complex, their impact on mosquito population genetics remains unclear. To this aim, we developed a PCR-RFLP test that discriminates the five known Wolbachia groups found in this host complex. We further examined the Wolbachia genetic diversity, the variability in the coinherited host mitochondria and their partitioning among members of the Cx. pipiens complex, in order to assess the impact of Wolbachia on host population structure. RESULTS: There was a strong association between Wolbachia and mitochondrial haplotypes indicating a stable co-transmission in mosquito populations. Despite evidence that members of the Cx. pipiens complex are genetically distinct on the basis of nuclear DNA, the association of Wolbachia and mtDNA with members of the Cx. pipiens complex were limited. The Wolbachia wPip-I group, by far the most common, was associated with divergent Cx. pipiens members, including Cx. quinquefasciatus, Cx. pipiens pipiens form pipiens and Cx. pipiens pipiens form molestus. Four other wPip groups were also found in mosquito populations and all were shared between diverse Cx. pipiens members. CONCLUSION: This data overall supports the hypothesis that wPip infections, and their allied mitochondria, are associated with regular transfers between Cx. pipiens members rather than specific host associations. Overall, this is suggestive of a recent and likely ongoing cytoplasmic introgression through hybridization events across the Cx. pipiens complex.

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.

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.

Forty years of erratic insecticide resistance evolution in the mosquito Culex pipiens.

One view of adaptation is that it proceeds by the slow and steady accumulation of beneficial mutations with small effects. It is difficult to test this model, since in most cases the genetic basis of adaptation can only be studied a posteriori with traits that have evolved for a long period of time through an unknown sequence of steps. In this paper, we show how ace-1, a gene involved in resistance to organophosphorous insecticide in the mosquito Culex pipiens, has evolved during 40 years of an insecticide control program. Initially, a major resistance allele with strong deleterious side effects spread through the population. Later, a duplication combining a susceptible and a resistance ace-1 allele began to spread but did not replace the original resistance allele, as it is sublethal when homozygous. Last, a second duplication, (also sublethal when homozygous) began to spread because heterozygotes for the two duplications do not exhibit deleterious pleiotropic effects. Double overdominance now maintains these four alleles across treated and nontreated areas. Thus, ace-1 evolution does not proceed via the steady accumulation of beneficial mutations. Instead, resistance evolution has been an erratic combination of mutation, positive selection, and the rearrangement of existing variation leading to complex genetic architecture.

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.

Multiple insecticide resistances in the disease vector Culex p. quinquefasciatus from Western Indian Ocean.

Several mosquito-borne diseases affect the Western Indian Ocean islands. Culex pipiens quinquefasciatus is one of these vectors and transmits filariasis, Rift Valley and West Nile viruses and the Japanese encephalitis. To limit the impact of these diseases on public health, considerable vector control efforts have been implemented since the 50s, mainly through the use of neurotoxic insecticides belonging to Organochlorines (OC), Organophosphates (OP) and pyrethroids (PYR) families. However, mosquito control failures have been reported on site, and they were probably due to the selection of resistant individuals in response to insecticide exposure. In this study, we used different approaches to establish a first regional assessment of the levels and mechanisms of resistance to various insecticides. Bioassays were used to evaluate resistance to various insecticides, enzyme activity was measured to assess the presence of metabolic resistances through elevated detoxification, and molecular identification of known resistance alleles was investigated to determine the frequency of target-site mutations. These complementary approaches showed that resistance to the most used insecticides families (OC, OP and PYR) is widespread at a regional scale. However, the distribution of the different resistance genes is quite heterogeneous among the islands, some being found at high frequencies everywhere, others being frequent in some islands and absent in others. Moreover, two resistance alleles displayed clinal distributions in Mayotte and La Réunion, probably as a result of a heterogeneous selection due to local treatment practices. These widespread and diverse resistance mechanisms reduce the capacity of resistance management through classical strategies (e.g. insecticide rotation). In case of a disease outbreak, it could undermine the efforts of the vector control services, as only few compounds could be used. It thus becomes urgent to find alternatives to control populations of Cx. p. quinquefasciatus in the Indian Ocean.

Tracking factors modulating cytoplasmic incompatibilities in the mosquito Culex pipiens.

Wolbachia are maternally inherited endosymbiotic bacteria that infect many arthropod species and may induce cytoplasmic incompatibility (CI), resulting in abortive embryonic development. One Wolbachia host, Culex pipiens complex mosquitoes, displays high levels of variability in both CI crossing types (cytotypes) and DNA markers. We report here an analysis of 14 mosquito strains, containing 13 Wolbachia variants, and with 13 different cytotypes. Cytotypes were Wolbachia-dependent, as antibiotic treatment rendered all strains tested compatible. Cytotype distributions were independent of geographical distance between sampling sites and host subspecies, suggesting that Wolbachia does not promote a reproductive isolation depending on these parameters. Backcross analysis demonstrated a mild restoring effect of the nuclear genome, indicating that CI is mostly cytoplasmically determined for some crosses. No correlation was found between the phenotypic and genotypic variability of 16 WO prophage and transposon markers, except for the WO prophage Gp15 gene, which encodes a protein similar to a bacterial virulence factor. However, Gp15 is partially correlated with CI expression, suggesting that it could be just linked to a CI gene.