Introduction
The Asian tiger mosquito,
Aedes albopictus, is a highly invasive mosquito endemic to Southeast Asia, and its adaptation to urban environments has led to its global spread [
1]. The first detection of this mosquito in Greece occurred in Corfu and Thesprotia in 2003/4, and since then, it has been found in nearly every prefecture [
2,
3]. Due to its ability to transmit many important human pathogens, such as Zika, Chinkungunya (CHIKV), Dengue (DENV- 1, 2, 3, and 4), and
West Nile (WNV) viruses, the Asian tiger mosquito has attracted substantial research interest [
4,
5]. With traditional use of chemical insecticides raising serious concern about toxic effects on organisms and ecosystems, as well as, about the increasing mosquito resistance, alternative environment-friendly methods for the control of mosquito populations have arisen [
6]. The Sterile Insect Technique (SIT), a genetic control method that has been used against agricultural pests since decades, is reconsidered as a promising approach against disease transmitting mosquito species, including
Ae. albopictus, exhibiting considerable progress towards this goal in recent years [
7]. Furthermore, symbiont-based approaches are gaining more interest as promising tools for vector control [
8].
Wolbachia pipientis is a Gram-negative, intracellular, endosymbiotic bacterium found in several arthropod species and filarial nematodes. It is classified into 17 supergroups A-R (except G) based on genomic diversity [
9]. Supergroups A and B are the most prevalent found in arthropod species [
10]. Initial reports proposed that 20% of arthropod species are infected by
Wolbachia [
11,
12], but a later meta-analysis suggested that the number is closer to 60% [
13]. The empirical evidence regarding the intraspecific prevalence of
Wolbachia infection seems to follow an “all or nothing” pattern, with either very low (<10%) or very high (>90) infection rates in each species [
13]. The spread of
Wolbachia in arthropod species is closely linked to the manipulations and phenotypes that the bacterium induces in its hosts. The most frequently observed reproductive phenotype is cytoplasmic incompatibility (CI). The maternally inherited
Wolbachia strain induces a modification during spermatogenesis in males, which can be reversed if a compatible strain from the mated female is present in the embryo. If the strains are not compatible, the eggs will not hatch. It is important to note that a non-infected male is compatible with any infected female. This phenotype explains the ability of
Wolbachia to establish infections in natural populations with a high prevalence, close to 100% [
10,
14,
15]. In mosquitoes, another important Wolbachia-associated observed phenotype is pathogen interference (PI). Certain
Wolbachia strains can reduce the replication and transmission of specific viruses and other pathogens. This phenotype has been observed in natural infections and has been most successfully demonstrated through transfections [
10,
16,
17]. [
10,
16,
17] In recent years,
Wolbachia has emerged as a promising tool for reducing mosquito-transmitted diseases, through the introduction of appropriately infected mosquitoes that can exert CI or PI effects in natural populations [
8]. It is widely accepted that careful consideration and thorough research on the effects of a particular
Wolbachia strain in a specific mosquito species must precede the release of mosquitoes in nature to prevent unintended consequences. Furthermore, extensive screening of natural populations for
Wolbachia prevalence and the strain composition is an equally important step before any
Wolbachia-based control effort can be considered [
8].
In the natural populations of
Ae. albopictus tested, mainly in Asia and south America,
Wolbachia infection is established with a high prevalence (>95%) [
18,
19,
20,
21]. Two commonly found
Wolbachia strains in
Ae. albopictus are
wAlbA and
wAlbB, belonging to supergroups A and B, respectively [
22]. These strains often co-infect the same individual, although single infections with
wAlbB, more commonly observed in males, have also been reported [
23]. Enviromental factors have been shown to affect the density of the
Wolbachia strains in
A. albopictus, with higher temperature and lower nutrition resulting to lower
Wolbachia density [
24]. Genetic diversity of the
wAlbA and
wAlbB strains has been studied by the MLST (MultiLocus Sequence Typing) system [
25], a typing system based on the polymorphism of five conserved genes (
gatB, coxA, hcpA, ftsZ and fbpA), in natural populations [
25]. Both strains were found stable with
wAlbA presenting a unique profile and
wAlbB two very similar ones [
18,
25,
26,
27]. In addition to the naturally occurring
Wolbachia strains in
Ae. albopictus,
wAlbA and
wAlbB, several
Wolbachia transfections have been explored in the species, including: (a) stable transfection with the wMel strain from
Drosophila melanogaster, which showed positive results in blocking DENV and CHIKV transmission [
28,
29], (b) the
wAu strain from
Drosophila simulans, which displayed resistance to Zika and DENV with a moderate fitness cost to the mosquitoes [
30], and (c) the
wPip strain from
Culex pipiens, which led to a reduction in field populations when triple-infected (
wAlbA,
wAlbB and
wPip) (line HC) or single infected (
wPip)(line ARwP) male mosquitoes were released in large field trials [
31,
32,
33,
34].
To date, the prevalence and genetic diversity of
Wolbachia strains in
Ae. albopictus populations in Europe remain poorly understood. To our knowledge only one study reports a prevalence of 79% in field population of the Iberian Peninsula, with no information provided about the type of stains (
wAlbA or
wAlbB) or genotypic diversity [
35]. Apart from a field trial of the IIT method with the ARwP line in Italy [
34], no
Wolbachia-based management practices are currently underway in Europe. Notably, most cases of viruses transmitted by
Ae. albopictus are not autochthonous [
36,
37]. However, sporadic autochthonous cases of DENV-1 in Italy, Croatia, and southern France [
38,
39,
40,
41], as well as CHIKV outbreaks in Italy and France [
40,
41,
42,
43,
44], and WNV cases in Italy, Germany, the Netherlands, France [
45,
46,
47,
48,
49,
50,
51], have been reported. The frequency of
Ae. albopictus-transmitted diseases is expected to increase in Europe in the coming years, driven by climate change [
52,
53,
54], and thus increased efforts should be made to study and control the mosquito populations, including the use of
Wolbachia-based methods.
In this study, our aim was to elucidate the unknown prevalence of
Wolbachia in
Ae. albopictus mosquitoes in Greece. As mentioned,
Ae. albopictus was first reported in Greece relatively recently (around 2003), and in the subsequent years, its populations have expanded throughout the country [
3]. We also aimed to determine the detected strains,
wAlbA or
wAlbB, and the occurrence of double and single infections. Additionally, we sought to shed light on the genetic diversity of
wAlbA or
wAlbB strains in northeastern Europe by comparing our MLST genotypes with those in other studies. Furthermore, we recorded the sex of the mosquitoes and the temperature environment of the studied populations to provide empirical evidence regarding the association between these factors and the propensity for double or single infections. Our experimental design involved a field sampling campaign across eight different prefectures in Greece between July and September 2021, during which we collected a total of 114 adult
Ae. albopictus mosquitoes. The sex of each mosquito was recorded, and we employed PCR and Sanger sequencing techniques to detect and genotype the
Wolbachia strains. To assess the impact of environmental temperature, we collected temperature recordings four times a day for 30 days prior to the sampling date. We then analyzed this data in relation to the observed patterns of occurrence of double/single infections. Overall, our findings provide valuable insights into the distribution and genetic diversity of
Wolbachia in Greek
Ae. albopictus populations, Furthermore, they highlight the importance of understanding the organisms’ short-term plastic and adaptive responses to environmental stressors and rapid climate change. We anticipate these results to contribute to the development of effective
Wolbachia-based approaches for controlling mosquito-borne diseases in
Ae. albopictus in Greece but also in other regions.
Discussion
High Wolbachia infection rate in Ae. albopictus population in Greece. The high infection rate (95%) in the natural population of
Ae. albopictus in Greece supports previous studies from other regions and confirms that
Wolbachia infection has been fixed in
Ae. albopictus worldwide [
13,
21,
65]. However, we observed slight differences in the infection rates among different prefectures. All the examined prefectures of Northern Greece exhibited an infection rate of 100% unlike the prefectures in central and southern Greece (88% mean rate). Variance in infection rates due to different geographical regions has been previously reported [
18]. Understanding the infection dynamics of
Wolbachia in
Ae. albopictus holds the promise of establishing biologically based management practices, as the density of the endosymbiont
Wolbachia plays an important role in crossing sterility (e.g. Werren et al. 1995 [
12]).
Low MLST profile diversity for Wolbachia strains in Greek populations and worldwide | The MLST profiles for the
wAlbA and
wAlbB strains were monomorphic in the tested Greek populations, which is consistent with the low diversity observed in the MLST genes in the
Ae. albopictus Wolbachia strains worldwide [
18,
25,
26,
27]. Specifically, the sole identified profile for the
wAlbB strain (242, 229, 166, 210, and 27 as per the numbering of the MLST database) was also the most common profile identified in almost all populations studied from China, Russia, and Argentina [
18,
25,
26,
27]. There is only one other profile reported for the
wAlbB strain in China, which differs only in the
gatB allele [
18]. The sole identified
wAlbA profile (3, 229, 2, 10, and 3) shares identical alleles at four out of the five loci with samples from China and Thailand [
18,
66]. It only differs in the
coxA locus, where we detected the same allele (229) as in the
wAlbB strain. Since it has been proposed that the 229 allele has moved from supergroup-A to supergroup-B through recombination [
26], it may be that the supergroup-A-specific primers mis-amplified the 229 allele of the
wAlbB strain as it is expected to be more abundant in the double infected individuals [
67]. Notably, as mentioned in Badieritakis et al. (2018) [
3], the Greek populations of
Ae. albopictus have been found to have some genetic affinity with those of Thailand [
68]. As such, our results seem to align with those earlier findings. The Greek populations of
Ae. albopictus have also been established relatively recently. For example, the first reports in Serres, Aetolia-Acarnania, and Thessaloniki were made between 2007 and 2009 [
3]. This recent occurrence of populations may have a role in the observed low MLST diversity due to founder effects and limited time for divergence.
Differential strain distribution of Wolbachia by sex | By examining the double and single infection patterns in relation to the sex of the individual mosquitoes, we noted that almost all female mosquitoes were double infected (92%), whereas only about half of males carried both strains (52%). We observed that in males, the single infection with the
wAlbB strain is more common in marginally non-significant level (
Figure 2B), which agrees with the scientific literature [
19,
20,
21,
23,
69]. There are different hypotheses for the detection of both strains in females and only the
wAlbB strain in males. It is believed that male mosquitoes lose their
wAlbA infection as they get older as an evolutionary selection of nuclear countermeasures to
Wolbachia manipulation [
23]. Female
Ae. albopictus carry a fixed infection with both strains, with rare instances of single infected females [
19,
20,
21,
23,
69]. The mechanism of CI suggests that females need both strains, to reverse the sperm modification each
Wolbachia strain induces in males to maximize mating compatibility [
10,
70]. However, males with a single or no infection remain compatible with double infected females [
10]. Male embryos carry both infections since their mothers are almost always double infected, but they tend to lose the
wAlbA infection as they age. The same pattern is not observed in the
wAlbB strain. The first detected
Wolbachia infected
Ae. albopictus (before 1970) carried only the
wAlbA strain, indicating that the infection with the
wAlbB strain is more recent [
65].
Differential strain distribution of Wolbachia by environmental temperature | The prefecture-specific infection rates highlighted that non-infected individuals were only found in the three southernmost prefectures examined (ACH, AET, KAR). Furthermore, single-infected individuals were also mainly found in these three prefectures (with one single-infected individual found in Imathia as well) (
Table 1). This observation aligns with previous results from China showing that
Wolbachia prevalence in
Ae albopictus natural populations is affected by geographical location and climate type [
18]. It also agrees with hypotheses that higher temperatures lead to decreased
Wolbachia density in laboratory populations, affecting both infections but particularly
wAlbA [
24,
71]. Our statistical analysis strongly supported this hypothesis, as the effect of the mean temperature from 30 days prior to the sampling date had a highly significant impact on the ratio of single infections in the mosquito populations (P = 6.50E-07;
Figure 2A). The mechanistic underpinnings of this observation are not fully understood. Some insight could be provided by studies concluding that maternal transmission under high temperatures and larval density is more often unsuccessful for the
wAlbA than for the
wAlbB strain [
67,
71]. Even so, when scoring single- to double-infected individuals it should not be disregarded that the density of the
wAlbA strain is smaller than that of
wAlbB in double-infected individuals [
18,
67] and, thus, it is possible that what appears like absence of the
wAlbA infection could in fact reflect a reduction beyond PCR detection limits [
15,
71]. Nevertheless, the above findings raise important points regarding the molecular ecology of the
Wolbachia infections, which may also be of applied interest.
Conclusions: Altogether: this study provides important insights into the distribution and diversity of Wolbachia strains in Ae. albopictus populations in Greece and the broader northeastern European region. We report the first MLST genotype information of Wolbachia strains in Ae. albopictus in Europe and identify a high infection rate in the natural population in Greece, confirming that Wolbachia infection has become fixed in populations worldwide. We found a positive correlation between higher environmental temperatures and single Wolbachia infections in Ae. albopictus mosquitoes, demonstrating that ecological and evolutionary factors play a role in Wolbachia dynamics. Overall, this study contributes to our understanding of Wolbachia dynamics in natural populations of Ae. albopictus. We anticipate that future research with even denser sampling will evaluate our findings, and it will have profound implications for the potential use of Wolbachia-based approaches in controlling mosquito-borne diseases in the area.