- Open Access
Wolbachia diversity and cytoplasmic incompatibility patterns in Culex pipiens populations in Turkey
Parasites & Vectorsvolume 11, Article number: 198 (2018)
Wolbachia are maternally transmitted bacteria that can manipulate their hosts’ reproduction causing cytoplasmic incompatibility (CI). CI is a sperm-egg incompatibility resulting in embryonic death. Due to this sterilising effect on mosquitoes, Wolbachia are considered for vector control strategies. Important vectors for arboviruses, filarial nematodes and avian malaria, mosquitoes of Culex pipiens complex are suitable for Wolbachia-based vector control. They are infected with Wolbachia wPip strains belonging to five genetically distinct groups (wPip-I to V) within the Wolbachia B supergroup. CI properties of wPip strongly correlate with this genetic diversity: mosquitoes infected with wPip strains from a different wPip group are more likely to be incompatible with each other. Turkey is a critical spot for vector-borne diseases due to its unique geographical position as a natural bridge between Asia, Europe and Africa. However, general wPip diversity, distribution and CI patterns in natural Cx. pipiens (s.l.) populations in the region are unknown. In this study, we first identified wPip diversity in Turkish Cx. pipiens (s.l.) populations, by assigning them to one of the five groups within wPip (wPip-Ito V). We further investigated CI properties between different wPip strains from this region.
We showed a wPip fixation in Cx. pipiens (s.l.) populations in Turkey by analysing 753 samples from 59 sampling sites. Three wPip groups were detected in the region: wPip-I, wPip-II and wPip-IV. The most dominant group was wPip-II. While wPip-IV was restricted to only two locations, wPip-I and wPip-II had wider distributions. Individuals infected with wPip-II were found co-existing with individuals infected with wPip-I or wPip-IV in some sampling sites. Two mosquito isofemale lines harbouring either a wPip-I or a wPip-II strain were established from a population in northwestern Turkey. Reciprocal crosses between these lines showed that they were fully compatible with each other but bidirectionally incompatible with wPip-IV Istanbul infected line.
Our findings reveal a high diversity of wPip and CI properties in Cx. pipiens (s.l.) populations in Turkey. Knowledge on naturally occurring CI patterns caused by wPip diversity in Turkey might be useful for Cx. pipiens (s.l.) control in the region.
First discovered in Culex pipiens (s.l.) mosquitoes , the α-proteobacterium Wolbachia pipientis is one of the most common vertically transmitted cytoplasmic symbionts. Indeed, meta-analysis predicts Wolbachia infection in up to 50% of the arthropod species . The success of their vertical transmission mainly relies on their host reproduction manipulation strategies; parthenogenesis induction, feminization, male killing and cytoplasmic incompatibility (CI) . CI, the most common of these strategies, is modelled by a modification-rescue (mod-resc) system where Wolbachia modifies sperm of infected males (mod function), and only a compatible Wolbachia strain in the eggs can rescue (resc function) this modification . Consequently, Wolbachia causes conditional sterility in crosses either between uninfected females and infected males  or like in the case of Cx. pipiens (s.l.), between females and males infected with incompatible Wolbachia strains [6, 7].
Consisting of several species, including important disease vectors with worldwide distribution (e.g. Cx. pipiens and Cx. quinquefasciatus) , mosquitoes of Cx. pipiens complex have a unique relationship with their endosymbiont Wolbachia (wPip). wPip is fixed in natural Cx. pipiens (s.l.) populations where they induce the most complex CI relationships yet described among arthropods, including uni and bidirectional incompatibility [7, 9,10,11,12]. So far, reciprocal crosses between many isofemale lines and four reference lines showed eight different mod and four different resc functions in Cx. pipiens (s.l.), resulting in the definition of 14 different cytotypes throughout the world .
In contrast to this observed diversity of CI patterns, Cx. pipiens wPip strains are closely related, and all belong to a clade within Wolbachia B supergroup [12,13,14]. However, recent studies of fast evolving markers showed the presence of many genetically distinct wPip strains in Cx. pipiens (s.l.) mosquitoes [10, 12, 15] distributed in five distinct phylogenetic groups (wPip-I to V) . Using a PCR/RFLP assay based on pk1 gene, encoding proteins with ankyrin motifs, a wPip strain can be assigned to one of these five groups [12, 16]. A study of the wPip worldwide distribution showed an important spatial structure of wPip groups . For instance, only wPip-I was found in sub-Saharan Africa, South America and Southeast Asia, while wPip-III was mainly observed in North America. Strains belonging to the wPip-II group were mostly found in western Europe and wPip-V in Asia. wPip-IV group strains exhibit a patchy distribution in Europe, North Africa and Asia . Also, Wolbachia genetic diversity and their CI patterns strongly correlate; most wPip strains from the same group render their host compatible with each other (except few unidirectional incompatibilities) whereas those from different groups often lead to unidirectional or bidirectional incompatibilities . Recently this huge diversity of CI patterns observed in Cx. pipiens has been explained by the amplification and the diversity of an operon in wPip strains’ genomes  composed of cidA and cidB genes involved in Wolbachia induced CI [18, 19]. No effect of host genetic background on the CI patterns  and no multiple infections by several strains have ever been shown [12, 16, 21].
Being a natural bridge between Africa, Asia and Europe, Turkey is a critical spot for many emerging and re-emerging vector-borne diseases [22, 23] and for the diversity of the vectors that transmit these diseases . For instance, high diversity and abundance of Cx. pipiens (s.l.) species have been recorded in the area including Cx. quinquefasciatus, Cx. pipiens and its physiological variant Cx. pipiens f. molestus . Arboviruses such as West Nile virus, mainly transmitted by these mosquitoes, have also been shown to circulate in Turkey [25,26,27,28]. Therefore, understanding wPip diversity and their CI properties of wPip to control Cx. pipiens (s.l.) populations in Turkey is a cornerstone for vector control in the region and prevention of putative epidemics extending through Europe, Asia and northern Africa. This knowledge can contribute to the biological vector control techniques using CI properties such as incompatible insect technique (IIT). IIT, the mass release of males harbouring incompatible Wolbachia into focal populations, has been shown to successfully decrease the female reproduction by sterilisation and reduce the pest/vector populations [29,30,31,32,33,34]. Nevertheless, our knowledge of the wPip genetic diversity and the CI patterns in Turkey is yet limited to only one line established with samples collected in Istanbul in 2003 .
Here, we collected and analysed 753 Cx. pipiens (s.l.) individuals (larvae and adults) from natural populations across Turkey. We studied (i) the wPip diversity in this geographically critical region in the crossroads of three continents, (ii) the CI relationships between Cx. pipiens lines from Turkey and (iii) the CI relationships between Turkish lines and reference lines to compare their CI properties to previously characterized mod (male crossing type) and resc (female crossing type) functions. Taken together, these results might be used in integrated vector control programs against Cx. pipiens (s.l.) in Turkey.
Sample collection and identification
A total of 753 samples from 59 different sampling sites in Turkey were tested for Wolbachia diversity. Most of these samples (n = 677) were collected during the larval stage, between July to September 2016 (Table 1). The rest of the samples has been collected as adults, using adult light traps, from May to September (2012–2015) (Table 1, Fig. 1). All of the sampling sites were situated outdoors with the only exception of sample site 16 (Table 1). Collected larvae and adults were morphologically identified as Culex pipiens (s.l.) / Cx. torrentium . As the samples were only morphologically identified, we used Cx. pipiens (s.l.) to refer to Cx pipiens assemblage that includes both Cx. pipiens and Cx. quinquefasciatus, their hybrids and physiological forms [21, 35]. Samples were stored in 70% ethanol until DNA extraction before testing them for the wPip presence and genetic characterisation.
For analysing the CI patterns induced by the wPip strains belonging to different groups found in Turkey, egg rafts and larvae were collected from a population in Thrace region of Turkey, in Tekirdag (Table 1 sampling site 52, Fig. 1). Collected larvae were reared to adults in insectary conditions (at 25 ± 2 °C and 75 ± 2% relative humidity and a 12:12 h photoperiod) and fed with a mixture of shrimp powder and rabbit pellets. To establish isofemale lines, females were fed with turkey blood using a Hemotek membrane feeding system (Discovery Workshops, Blackburn, United Kingdom) and were allowed to lay eggs five days later. Each egg raft (100–300 eggs) was isolated for hatching, and the isofemale line was established using resulting sibling larvae. A pool of first-instar larvae (L1) was tested to identify the wPip group. Two isofemale lines each harbouring either a wPip-I or wPip-II strain were reared for further crossing experiments in 65 dm3 cages in insectary conditions and were fed with a honey solution and a weekly blood meal.
Crossing experiments between Turkish Culex pipiens lines
Tek wPip-I and Tek wPip-II isofemale mosquito lines were reared for at least four generations in insectary conditions to allow their acclimatisation before crossing experiments. Mosquitoes were isolated during pupal stage, and emerging adults were sexed. Then, 2–5 days old virgin males and females (n = 25–50) were used to carry out reciprocal crosses between them and with Istanbul wPip-IV line. Females were fed with turkey blood using a Hemotek membrane feeding system (Discovery Workshops, United Kingdom) on the sixth day following caging and were allowed to lay eggs five days after the blood meal. Egg rafts were then isolated individually in 24 well plates filled with tap water until hatching. Embryonic development of all the unhatched egg rafts was verified to differentiate between non-fertilized egg rafts and CI induced embryonic death as previously described .
The crossing relationships were identified as following : compatible (C) when > 90% of the rafts hatched in the two reciprocal crosses; and incompatible (IC), with two CI patterns: (i) unidirectionally incompatible crosses: when between 0–10% of the rafts hatched in one of the reciprocal crosses and > 90% in the other; and (ii) bidirectionally incompatible crosses: when less than 10% of the rafts hatched in both reciprocal crosses.
Crossing experiments to infer mod and resc functions
The mod (male crossing type) and resc (female crossing type) functions caused by many wPip strains, which belong to different wPip groups (I-V), have been identified by reciprocal crosses with 4 reference lines: LaVar (wPip-II), MaClo, Slab (wPip-III) and Istanbul (wPip-IV) . Here, we used same four reference lines to define the mod-resc functions of Tek wPip-I and Tek wPip-II isofemale lines and to compare them to previously defined ones .
Identification of Wolbachia diversity
DNA was extracted from the samples using CTAB method . PCR assays were conducted using pk1 primers (PK1 Forward: 5'-CCA CTA CAT TGC GCT ATA GA-3' and PK1 Reverse: 5'-ACA GTA GAA CTA CAC TCC TCC A-3'-AM397079 ), which amplify a 1.3-kilobase (kb) fragment from ankyrin domain coding gene of Wolbachia. PCR amplifications were made in following conditions: initial denaturation for 5 min at 94 °C, followed by 35 cycles of denaturation, annealing and elongation respectively at 94 °C for 30 s, 52 °C for 30s, and 72 °C for 90 s, and a final elongation at 72 °C for 5mn. Resulting PCR products then used in RFLP assays first with TaqαI enzyme to discriminate specific wPip alleles “a” or “e” (wPip- I or wPip-V; 991, 251, 107 bp), “b” (wPip-III; 669, 665 bp), “c” (wPip-II; 851, 498 bp) and “d” (wPip-IV; 497, 251, 107 bp) [7, 16]. Secondly, since TaqαI digestion of “a” and “e” alleles show the same digestion pattern, pk1 PCR products of the samples showing this pattern were digested with PstI enzyme to further discriminate “a” (wPip- I; 903, 303, 141 bp) and “e” (wPip-V; 903, 430 bp) alleles [7, 16]. Digested amplified fragments were separated by agarose gel electrophoresis (2%), stained with ethidium bromide (1 μg/ml) and visualized with UV light. Samples from laboratory mosquito lines with different wPip groups and tetracycline-cured Wolbachia negative lines were included in every reaction as positive and negative controls, respectively, and always gave the expected result.
The occurrence of different wPip groups was compared by a Chi-square test using R software (version 3.3.1).
Diversity and distribution of wPip groups
Wolbachia wPip was present in all of the 753 Cx. pipiens (s.l.) individuals tested and they were further identifiable to one of the five previously described groups (wPip-I to V). Co-infection of one individual by different wPip groups was never observed. Out of five wPip groups identified so far in the world, three of them (i.e. wPip-I-II and IV) were represented in the studied area (Fig. 1). The abundance of these groups was significantly different from each other (χ2 = 474.99, df = 2, P < 0.0001). While wPip-II was the most dominant (n = 500, 66% of the samples, Table 1, Fig. 1) and widespread group (39 sampling sites out of 59 total) in Turkey; wPip-IV was found only in two locations, both in Thrace Region (in Edirne, sampling site 21 and Istanbul sampling site 34, Table 1; Fig. 1), and was the least abundant group (n = 12, 1% of the samples, Table 1; Fig. 1). wPip-I was found in 31 locations and a total of 241 individuals.
Co-existence of wPip strains in different individuals from the same sampling sites
In 20 % of the sampling sites wPip-I and wPip-II co-existed (Table 1, Fig. 1). wPip-IV was only found co-existing with wPip-II in one sampling site but never found in the same sampling site with wPip-I, even though they were sampled from nearby sites (~8 km) in north western Turkey (in Edirne, sampling site 21 and 23, Table 1; Fig. 1).
Naturally occurring CI patterns in Turkey
Two isofemale lines (Tek wPip-I and Tek wPip-II), harbouring two different wPip strains from two different groups, were established from north western Turkey (Tekirdag Province, Table 1 sampling site 52, Fig. 1) to identify CI patterns caused by different wPip groups in the region. Reciprocal crosses between these lines showed that Tek wPip-I and Tek wPip-II were fully compatible with each other (Table 2). Both lines were bidirectionally incompatible with the line harbouring Istanbul strain (wPip-IV, Table 2).
Mod and resc properties of Turkish wPip strains
To compare mod and resc functions of Turkish wPip strains with mod and resc functions of worldwide collected wPip strains, we performed reciprocal crosses of Tek wPip-I and Tek wPip-II lines with the four reference lines [LaVar (wPip-II), MaClo (wPip-III), Slab (wPip-III) and Istanbul (wPip-IV)]. Tek wPip-I males were compatible with LaVar (wPip-II) and MaClo (wPip-III) females while incompatible with Slab (wPip-III) and Istanbul (wPip-IV) females (Table 2). This type of mod property, inferred from similar crosses, has already been shown for the wPip-I group from Tunisia; numbered “vi” . Contrarily, Tek wPip-II males demonstrated a new mod property, as they were incompatible with LaVar and Istanbul, and compatible with MaClo and Slab females (Table 2). We numbered this new mod as “ix” to continue the previously published numeration . Both Tek wPip-I and Tek wPip-II lines showed the same resc type, which was characterised by the compatible crosses of females of these lines with all the males from the reference lines except Istanbul (Table 2). This resc type (resc “2”) is the most common resc type found worldwide for wPip-I and wPip-II groups .
In Turkey, all tested Cx. pipiens were infected with Wolbachia wPip. Such fixation of wPip has been demonstrated worldwide, including in the neighbouring country Iran , in Cx. pipiens and Cx. quinquefasciatus populations [35, 40,41,42,43,44]. A previous study in Turkey, however, showed a lower prevalence of wPip . This might be caused by the misidentification of a recently described cryptic species within Cx. pipiens complex that has been shown to lack wPip infection and to be reproductively isolated from the other members of the complex [46, 47]. Similarly, Cx. torrentium, which is difficult to differentiate morphologically from Cx. pipiens (s.l.) mosquitoes is not infected with Wolbachia [41, 44, 48]. Therefore 100% wPip infection rate of our samples confirmed that we only analysed Cx. pipiens complex members (excluding both previously mentioned cryptic species and Cx. torrentium) in the present study.
The identification of the Cx. pipiens taxa was left out of the scope of this study for several reasons. Previous studies on the diversification of wPip in Cx. pipiens (s.l.) have proved that their diversity is not directly related to the nuclear genetic background of the mosquitoes, meaning that no wPip group was specific for a Cx. pipiens sibling species [10, 21, 49]. It rather follows the same distribution as mitochondrial diversity (mtDNA) of mosquitoes, as wPip are maternally transmitted to the next generation through the egg cytoplasm along with mitochondria [16, 46]. Moreover, CI properties are independent of the genetic background of Cx. pipiens (s.l.) and directly dictated by their Wolbachia . Recent studies on Cx. pipiens (s.l.) in Turkey had shown that both Cx. quinquefasciatus, Cx. pipiens and its form Cx. pipiens f. molestus, were present in Turkey . The co-existence of these sibling species in same sampling sites [24, 50] and the existence of hybrids [50,51,52] suggest that they can exchange wPip strains easily in natural populations.
We have identified three different wPip groups, i.e. wPip-I-II and IV in Cx. pipiens (s.l.) mosquito populations in Turkey. The only previous sample from Turkey, which has been assigned to wPip groups, was a wPip-IV group strain collected in Istanbul in 2003 [16, 35]. Other than this single case, the wPip diversity in Turkish Cx. pipiens populations was to date completely unknown. Although wPip diversity was investigated in regions around Turkey . Dumas et al.  have found wPip-I strains in Middle East (Lebanon, Israel, Jordan) and in northern Africa (Tunisia), wPip-II strains widely distributed in eastern Europe and Cyprus, and wPip-IV strains in a patchy distribution in Europe -in areas dominated by other wPip groups. We demonstrated that wPip strains belonging to wPip-I, wPip-II and wPip-IV, previously identified near Turkey, are all present in the Cx. pipiens (s.l.) populations within this country, suggesting that Turkey is a crossroads for wPip strains from eastern Europe, Africa and Middle East as for their vector hosts.
The most widespread groups in Turkey, wPip-I and wPip-II, induce reciprocal compatibility between their hosts and co-exist in many populations. Indeed, different wPip strains can co-exist in a single natural Cx. pipiens (s.l.) population [12, 35, 43, 53] and these coexisting strains are usually compatible with each other . Mathematical models confirm that only compatible strains can stably coexist in unstructured and panmictic host populations when the fitness costs related to infection by different wPip strains are the same [54, 55]. When bidirectional incompatibility inducing wPip strains co-exist in one population the most prevalent strain is expected to eventually invade the population [54, 55] and when unidirectional incompatibility inducing wPip strains co-exist, CI-inducing strain is expected to invade the population once above a frequency threshold [49, 55]. Therefore, a stable co-existence of incompatible strains is predicted to be rare. However, an example of the co-existence of unidirectional incompatibility inducing strains, belonging to a wPip-I group and wPip-IV group, has been shown in Tunisia while mathematical models predicted that wPip-I should have invaded this area in only 4 generations [5, 49]. Atyame et al.  hypothesized that low dispersal and extinction-recolonization events could explain this stable co-existence. In one site in Turkey, we observed the co-existence of bidirectional incompatibility inducing strains (wPip-IV and wPip-II). The low prevalence of wPip-IV and its incompatibilities suggest that it should disappear from the population. However, we have evidence that wPip-IV strains were already present at least 13 years ago at Istanbul since it has been sampled in 2003 . This persistence of wPip-IV, at low frequencies, could be explained by higher fitness costs associated with wPip-I and wPip-II infections or by extinction-recolonization events of wPip-IV-infected individuals as it has been suspected in Tunisia [43, 49, 55]. Fitness difference could be, for instance, linked to differences in fecundity [56, 57] or to a possible ability of the different wPip strains to protect their hosts against other microbial infections [58,59,60,61,62]. Further studies on the differences between wPip strains in terms of infection costs and pathogen protection might help to understand stable co-existence of bidirectionally incompatible wPip-IV strains observed in Turkey.
To study the phenotypical diversity of crossing types in Turkey, we crossed Turkish isofemale lines harbouring wPip-I and wPip-II strains with four reference lines defined by Atyame et al. . We inferred both their mod and resc functions and compared them to the eight mod and four resc functions already described worldwide. The Tek wPip-I line showed the most common resc functions for a wPip-I infected line (i.e. resc 2) but a rare mod function previously defined in few lines harbouring wPip-I or wPip-II strains (i.e. mod “vi”). The Tek wPip-II line showed the same resc function as wPip-I (i.e. resc 2) but a totally new mod function (i.e. mod “ix”). Our findings are consistent with theoretical predictions and empirical data suggesting new mod functions can more easily evolve and spread in the population than new resc functions [7, 63].
Natural CI properties induced by Wolbachia can be used to control the vector populations: the mass release of males harbouring incompatible Wolbachia into the natural populations can decrease the female reproduction and eradicate the pest/ vector populations (IIT) [29, 30]. Indeed, wPip induced CI has been used against Cx. quinquefasciatus (formerly named Cx. pipiens fatigans) for the first time in 1967 to control filariasis in Southeast Asia . More recently, natural CI properties caused by wPip infection have been found promising to control Cx. pallens (no longer considered as a valid species) in China  and Cx. pipiens populations in La Réunion Island [33, 34]. In the latter study, a wPip-IV strain from Istanbul has been successfully used to sterilise wPip-I females in semi-field conditions. We demonstrated that this Istanbul strain also induces bidirectional incompatibility with mosquitoes harbouring wPip-I or wPip-II in Turkey. This means that most Cx. pipiens females in Turkey, except in few sites in the Thrace region, can be sterilised by the release of males infected with Istanbul strain. Although further studies on intrapopulation CI variability, mating choice, hatching rate and population dynamics in semi field populations are needed for wPip-IV Istanbul to be used in future vector control programs in Turkey, a critical region for vector-borne diseases, our results suggest that it could constitute a good candidate.
We identified wPip diversity in natural Cx. pipiens (s.l.) populations in Turkey. The previously described wPip-IV group was in fact restricted to only two populations while wPip-I and wPip-II group are widely distributed and coexist in many populations all over the country. The wPip-IV strain Istanbul was found bidirectionally incompatible with individuals harbouring wPip-I or wPip-II from Turkey. This highlights the potential of wPip-IV harbouring males as a vector control to sterilise local Cx. pipiens populations, particularly where only wPip-I or wPip-II harbouring females were found.
- Mod :
- Resc :
Hertig M. The Rickettsia, Wolbachia pipientis (gen. et sp. n.) and associated inclusions of the mosquito, Culex pipiens. Parasitology. 1936;28:453–86.
Weinert LA, Araujo-Jnr EV, Ahmed MZ, Welch JJ. The incidence of bacterial endosymbionts in terrestrial arthropods. P Roy Soc B Biol Sci. 2015;282:20150249.
Werren JH, Baldo L, Clark ME. Wolbachia: master manipulators of invertebrate biology. Nat Rev Microbiol. 2008;6:741–51.
Werren JH. Biology of Wolbachia. Annu Rev Entomol. 1997;124:587–609.
Hoffmann a. a., Turelli M, Harshman LG. Factors affecting the distribution of cytoplasmic incompatibility in Drosophila simulans. Genetics. 1990;126:933–48.
Yen JH, Barr AR. New hypothesis of the cause of cytoplasmic incompatibility in Culex pipiens L. Nature. 1971;232:657–8.
Atyame CM, Labbé P, Dumas E, Milesi P, Charlat S, Fort P, et al. Wolbachia divergence and the evolution of cytoplasmic incompatibility in Culex pipiens. PLoS One. 2014;9:e87336.
Farajollahi A, Fonseca DM, Kramer LD, Marm Kilpatrick A. “Bird biting” mosquitoes and human disease: a review of the role of Culex pipiens complex mosquitoes in epidemiology. Infect Genet Evol. 2011;11:1577–85.
Laven H. Speciation and evolution in Culex pipiens. In: Wright JW, Pal R, editors. Genetics of Insect Vectors of Disease. Amsterdam: Elsevier; 1967. p. 251–75.
Duron O, Bernard C, Unal S, Berthomieu A, Berticat C, Weill M. Tracking factors modulating cytoplasmic incompatibilities in the mosquito Culex pipiens. Mol Ecol. 2006;15:3061–71.
Wolfgang A, Markus R, Dimitrios NA, Christian S. Evidence for low-titre infections in insect symbiosis: Wolbachia in the bark beetle Pityogenes chalcographus (Coleoptera, Scolytinae). Environ Microbiol. 2009;11:1923–33.
Atyame CM, Delsuc F, Pasteur N, Weill M, Duron O. Diversification of Wolbachia endosymbiont in the Culex pipiens mosquito. Mol Biol Evol. 2011;28:2761–72.
Guillemaud T, Pasteur N, Rousset F. Contrasting levels of variability between cytoplasmic genomes and incompatibility types in the mosquito Culex pipiens. P Roy Soc B Biol Sci. 1997;264:245–51.
Baldo L, Hotopp JCDD, Jolley KA, Bordenstein SR, Biber SA, Choudhury RR, et al. Multilocus sequence typing system for the endosymbiont Wolbachia pipientis. Appl Environ Microbiol. 2006;72:7098–110.
Duron O, Boureux A, Echaubard P, Berthomieu A, Berticat C, Fort P, et al. Variability and expression of ankyrin domain genes in Wolbachia variants infecting the mosquito Culex pipiens. J Bacteriol. 2007;189:4442–8.
Dumas E, Atyame CM, Milesi P, Fonseca DM, Shaikevich EV, Unal S, et al. Population structure of Wolbachia and cytoplasmic introgression in a complex of mosquito species. BMC Evol Biol. 2013;13:181.
Bonneau M, Atyame C, Beji M, Justy F, Cohen-Gonsaud M, Sicard M, et al. Culex pipiens crossing type diversity is governed by an amplified and polymorphic operon in Wolbachia genome. Nat Commun. 2018;9:319.
LePage DP, Metcalf JA, Bordenstein SR, On J, Perlmutter JI, Shropshire JD, et al. Prophage WO genes recapitulate and enhance Wolbachia-induced cytoplasmic incompatibility. Nature. 2017;543:243.
Beckmann JF, Ronau JA, Hochstrasser M, Tillett D, Ginalski K. A Wolbachia deubiquitylating enzyme induces cytoplasmic incompatibility. Nat Microbiol. 2017;2:17007.
Duron O, Bernard J, Atyame CM, Dumas E, Weill M. Rapid evolution of Wolbachia incompatibility types. P Roy Soc B Biol Sci. 2012;279:4473–80.
Atyame CM, Duron O, Tortosa P, Pasteur N, Fort P, Weill M. Multiple Wolbachia determinants control the evolution of cytoplasmic incompatibilities in Culex pipiens mosquito populations. Mol Ecol. 2011;20:286–98.
Inci A, Yildirim A, Duzlu O, Doganay M, Aksoy S. Tick-borne diseases in Turkey: A review based on one health perspective. PLoS Neglect Trop D. 2016;10:e0005021.
Failloux A-B, Bouattour A, Faraj C, Gunay F, Haddad N, Harrat Z, et al. Surveillance of arthropod-borne viruses and their vectors in the Mediterranean and Black Sea regions within the MediLabSecure Network. Curr Trop Med Rep. 2017;4:27–39.
Gunay F, Alten B, Simsek F, Aldemir A, Linton Y-M. Barcoding Turkish Culex mosquitoes to facilitate arbovirus vector incrimination studies reveals hidden diversity and new potential vectors. Acta Trop. 2015;143:112–20.
Ergunay K, Gunay F, Oter K, Kasap OE, Orsten S, Akkutay AZ, et al. Arboviral surveillance of field-collected mosquitoes reveals circulation of West Nile virus lineage 1 strains in Eastern Thrace, Turkey. Vector-Borne Zoonot. 2013;13:744–52.
Ergunay K, Gunay F, Erisoz Kasap O, Oter K, Gargari S, Karaoglu T, et al. Serological, molecular and entomological surveillance demonstrates widespread circulation of West Nile virus in Turkey. PLoS Neglect Trop D. 2014;8:1–10.
Ergunay K, Bakonyi T, Nowotny N, Ozkul A. Close relationship between West Nile virus from Turkey and lineage 1 strain from Central African Republic. Emerg Infect Dis. 2015;21:352–5.
Öncü C, Brinkmann A, Günay F, Kar S, Öter K, Sarıkaya Y, et al. West Nile virus, Anopheles flavivirus, a novel flavivirus as well as Merida-like rhabdovirus Turkey in field-collected mosquitoes from Thrace and Anatolia. Infect Genet Evol. 2018;57:36–45.
Knipling EF, Laven H, Craig GB, Pal R, Kitzmiller JB, Smith CN, et al. Genetic control of insects of public health importance. Bull World Health Organ. 1968;38:421–38.
Zabalou S, Riegler M, Theodorakopoulou M, Stauffer C, Savakis C, Bourtzis K. Wolbachia-induced cytoplasmic incompatibility as a means for insect pest population control. Proc Natl Acad Sci USA. 2004;101:15042–5.
Laven H. Eradication of Culex pipiens fatigans through cytoplasmic incompatibility. Nature. 1967;216:383–4.
Chen L, Zhu C, Zhang D. Naturally occurring incompatibilities between different Culex pipiens pallens populations as the basis of potential mosquito control measures. PLoS Neglect Trop D. 2013;7:e2030.
Atyame CM, Pasteur N, Dumas E, Tortosa P, Tantely ML, Pocquet N, et al. Cytoplasmic incompatibility as a means of controlling Culex pipiens quinquefasciatus mosquito in the islands of the South-Western Indian Ocean. PLoS Neglect Trop D. 2011;5:e1440.
Atyame CM, Cattel J, Lebon C, Flores O, Dehecq J-S, Weill M, et al. Wolbachia-based population control strategy targeting Culex quinquefasciatus mosquitoes proves efficient under semi-field conditions. PLoS One. 2015;10:e0119288.
Duron O, Lagnel J, Raymond M, Bourtzis K, Fort P, Weill M. Transposable element polymorphism of Wolbachia in the mosquito Culex pipiens: evidence of genetic diversity, superinfection and recombination. Mol Ecol. 2005;14:1561–73.
Schaffner E, Angel G, Geoffroy B, Hervy JP, Rhaiem A, Brunhes J. The mosquitoes of Europe. Montpellier. France: Institute de Recherche pour le Développement (IRD); 2001.
Duron O, Weill M. Wolbachia infection influences the development of Culex pipiens embryo in incompatible crosses. Heredity. 2006;96:493–500.
Rogers SO, Bendich AJ. Extraction of DNA from plant tissues. In: Gelvin SB, Schilperoort RA, Verma DPS, editors. Plant Molecular Biology Manuel. Dordrecht: Springer Netherlands; 1989. p. 73–83.
Karami M, Moosa-kazemi SH, Oshaghi MA, Vatandoost H, Sedaghat MM, Rajabnia R. Wolbachia endobacteria in natural populations of Culex pipiens of Iran and its phylogenetic congruence. J Arthropod-Borne Di. 2016;10:349–65.
Rasgon JL, Scott TW. Cytoplasmic incompatibility in the California Culex pipiens mosquito species complex: parameter estimates and infection dynamics in natural populations. Genetics. 2003;165:2029–38.
Fedorova MV, Shaikevich EV. Morphological and molecular-genetic distinctions between adult mosquitoes Culex torrentium Martini and C. pipiens Linnaeus (Diptera, Culicidae) from Moscow Province. Entomol Rev. 2007;87:127–35.
Khrabrova NV, Bukhanskaya ED, Sibataev AK, Volkova TV. The distribution of strains of endosymbiotic bacteria Wolbachia pipientis in natural populations of Culex pipiens mosquitoes (Diptera: Culicidae). Eur Mosq Bull. 2009;27:18–22.
Duron O, Raymond M, Weill M. Many compatible Wolbachia strains coexist within natural populations of Culex pipiens mosquito. Heredity. 2011;106:986–93.
Raharimalala FN, Boukraa S, Bawin T, Boyer S, Francis F. Molecular detection of six (endo-) symbiotic bacteria in Belgian mosquitoes: first step towards the selection of appropriate paratransgenesis candidates. Parasitol Res. 2016;115:1391–9.
Yildirim A, Inci A, Duzlu O, Onder Z, Ciloglu A. Detection and molecular characterization of the Wolbachia endobacteria in the Culex pipiens (Diptera: Culicidae) specimens collected from Kayseri province of Turkey. Ankara Univ Vet Fak. 2013;60:189–94.
Rasgon JL, Cornel AJ, Scott TW. Evolutionary history of a mosquito endosymbiont revealed through mitochondrial hitchhiking. P Roy Soc B Biol Sci. 2006;273:1603–11.
Dumas E, Atyame CM, Malcolm CA, Le Goff G, Unal S, Makoundou P, et al. Molecular data reveal a cryptic species within the Culex pipiens mosquito complex. Insect Mol Biol. 2016;25:800–9.
Leggewie M, Krumkamp R, Badusche M, Heitmann A, Jansen S, Schmidt-Chanasit J, et al. Culex torrentium mosquitoes from Germany are negative for Wolbachia. Med Vet Entomol. 2017:1–6.
Atyame CM, Labbé P, Rousset F, Beji M, Makoundou P, Duron O, et al. Stable coexistence of incompatible Wolbachia along a narrow contact zone in mosquito field populations. Mol Ecol. 2015;24:508–21.
Shaikevich EV, Vinogradova EB, Bouattour A, Gouveia de Almeida AP. Genetic diversity of Culex pipiens mosquitoes in distinct populations from Europe: contribution of Cx. quinquefasciatus in Mediterranean populations. Parasit Vectors. 2016;9:47.
Shaikevich EV, Vinogradova EB. The discovery of a hybrid population of mosquitoes of the Culex pipiens L. complex (Diptera, Culicidae) on the Kos Island (Greece) by means of molecular markers. Entomol Rev. 2014;94:35–9.
Smith JL, Fonseca DM. Rapid assays for identification of members of the Culex (Culex) pipiens complex, their hybrids, and other sibling species (Diptera: Culicidae). Am J Trop Med Hyg. 2004;70:339–45.
Duron O, Fort P, Weill M. Hypervariable prophage WO sequences describe an unexpected high number of Wolbachia variants in the mosquito Culex pipiens. P Roy Soc B Biol Sci. 2006;273:495–502.
Rousset F, Raymond M, Kjellberg F. Cytoplasmic incompatibility in the mosquito Culex pipiens: how to explain a cytotype polymorphism? J Evol Biol. 1991;4:69–81.
Engelstädter J, Telschow A. Cytoplasmic incompatibility and host population structure. Heredity. 2009;103:196–207.
Poinsot D, Merçot H. Wolbachia infection in Drosophila simulans: does the female host bear a physiological cost? Evolution. 1997;51:180–6.
Bordenstein SR, Werren JH. Do Wolbachia influence fecundity in Nasonia vitripennis? Heredity. 2000;84:54–62.
Hedges LM, Brownlie JC, O’Neill SL, Johnson KN. Wolbachia and virus protection in insects. Science. 2008;322:702.
Brownlie JC, Johnson KN. Symbiont-mediated protection in insect hosts. Trends Microbiol. 2009;17:348–54.
Hughes GL, Koga R, Xue P, Fukatsu T, Rasgon JL. Wolbachia infections are virulent and inhibit the human malaria parasite Plasmodium falciparum in Anopheles gambiae. PLoS Pathog. 2011;7:3–10.
Braquart-Varnier C, Altinli M, Pigeault R, Chevalier FD, Grève P, Bouchon D, et al. The mutualistic side of Wolbachia-isopod interactions: Wolbachia mediated protection against pathogenic intracellular bacteria. Front Microbiol. 2015;6:1388.
Teixeira L, Ferreira Á, Ashburner M. The bacterial symbiont Wolbachia induces resistance to RNA viral infections in Drosophila melanogaster. PLoS Biol. 2008;6:2753–63.
Charlat S, Calme C, Merçot H. On the mod resc model and the evolution of Wolbachia compatibility types. Genetics. 2001;159:1415–22.
We are thankful to Dr Hilal Bedir (Kafkas University, Kars, Turkey), Dr Berna Demirci (Kafkas University, Kars, Turkey) and Dr Mustafa Akiner (RTE University, Rize, Turkey) for their guidance during fieldwork, to Dr Sırrı Kar (Namık Kemal University, Tekirdağ, Turkey) for the collection of Cx. pipiens from Tekirdağ (Turkey), to Patrick Makoundou, Sandra Unal and Marco Perriat-Sanguinet from team EVAS (ISEM) for technical support.
This work was funded by the French ANR program (project ANR-16-CE02-0006-01 "CIAWOL", 2017-2020). Mine Altinli received support from Ministère des affaires étrangères (Campus France) via the French Embassy in Turkey and Project "ISEM Sud", 2016 for her PhD.
Availability of data and materials
All data generated or analysed during this study are included in this published article.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.