Hypervariable prophage WO sequences describe an unexpected high number of Wolbachia variants 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. Among all the described host species, mosquitoes of the Culex pipiens complex display the highest variability of CI crossing types. Paradoxically, searches for polymorphism in Wolbachia infecting strains and field populations hitherto failed or produced very few markers. Here, we show that an abundant source of the long-sought polymorphism lies in WO prophage sequences present in multiple copies dispersed in the genome of Wolbachia infecting C. pipiens (wPip). We identified up to 66 different Wolbachia variants in C. pipiens strains and field populations and no occurrence of superinfection was observed. At least 49 different Wolbachia occurred in Southern Europe C. pipiens populations, and up to 10 different Wolbachia were even detected in a single population. This is in sharp contrast with North African and Cretan samples, which exhibited only six variants. The WO polymorphism appeared stable over time, and was exclusively transferred maternally. Interestingly, we found that the CI pattern previously described correlates with the variability of Gp15, a prophage protein similar to a bacterial virulence protein. WO prophage sequences thus represent variable markers that now open routes for approaching the molecular basis of CI, the host effects, the structure and dynamics of Wolbachia populations.