# Impact of temperature on dengue and chikungunya transmission by the mosquito Aedes albopictus

Aurélien Mercier, Thomas Obadia, Davide Carraretto, Enkelejda Velo, Gaelle Gabiane, Silvia Bino, Marie Vazeille, Giuliano Gasperi, Catherine Dauga, Anna Malacrida, et al.

► To cite this version:

Aurélien Mercier, Thomas Obadia, Davide Carraretto, Enkelejda Velo, Gaelle Gabiane, et al. Impact of temperature on dengue and chikungunya transmission by the mosquito Aedes albopictus. Scientific Reports, 2022, 12 (1), pp.6973. 10.1038/s41598-022-10977-4. pasteur-03655009

HAL Id: pasteur-03655009

https://pasteur.hal.science/pasteur-03655009

Submitted on 29 Apr 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Impact of temperature on dengue and chikungunya transmission by the mosquito *Aedes albopictus*

Aurélien Mercier<fnref n="1" />,⁷, Thomas Obadia²,³, Davide Carraretto D⁴, Enkelejda Velo⁵, Gaelle Gabiane⁶, Silvia Bino⁵, Marie Vazeille⁶, Giuliano Gasperi⁴, Catherine Dauga⁶, Anna R. Malacrida⁴, Paul Reiter<fnref n="1" />,⁸ & Anna-Bella Failloux⁶,⁸✉

The mosquito *Aedes albopictus* is an invasive species first detected in Europe in Albania in 1979, and now established in 28 European countries. Temperature is a limiting factor in mosquito activities and in the transmission of associated arboviruses namely chikungunya (CHIKV) and dengue (DENV). Since 2007, local transmissions of CHIKV and DENV have been reported in mainland Europe, mainly in South Europe. Thus, the critical question is how far north transmission could occur. In this context, the Albanian infestation by *Ae. albopictus* is of interest because the species is present up to 1200 m of altitude; this allows using altitude as a proxy for latitude. Here we show that *Ae. albopictus* can transmit CHIKV at 28 °C as well as 20 °C, however, the transmission of DENV is only observed at 28 °C. We conclude that if temperature is the key environmental factor limiting transmission, then transmission of CHIKV, but not DENV is feasible in much of Europe.

Numerous arboviruses affecting humans are already endemic in continental Europe where they cause sporadic cases during a limited time window fitting with vectors activity; those are West Nile virus, Usutu virus, Tick-borne encephalitis virus, Crimean-Congo hemorrhagic fever virus, and several phleboviruses including Toscana virus<fnref n="1" />. In the last decades, new arboviruses have been introduced in Europe where invasive mosquito species are well established². Among them, five *Aedes* species are described including *Aedes albopictus*. This mosquito has been first detected in Europe, Albania in 1979³ and is now present in 28 European countries from Spain to Romania through Germany⁴. Originally from Asia, its high ecological plasticity (i.e. colonizing artificial as well as natural breeding sites) and physiological characteristics (i.e. eggs able to resist to desiccation and freezing) allow this species to spread and successfully establish in both tropical and temperate regions⁴. As expected, *Ae. albopictus*-associated arboviruses were soon after detected in Europe: chikungunya in 2007 in Italy⁵, dengue in 2010 in France⁶ and lastly, Zika in 2019 in France⁷.

Chikungunya hit Italy in 2007 causing more than 200 laboratory-confirmed human cases in Ravenna province from July to September⁸. An East-Central-South-African genotype of chikungunya virus (CHIKV; *Alphavirus*, Togaviridae) was introduced by a traveler coming back from Kerala, India. The virus had the substitution A- > V at the position 226 in the E1 glycoprotein conferring an enhanced replication and transmission by *Ae. albopictus*⁹,¹⁰. Besides, dengue was not an unknown disease in Europe. The last documented outbreak happened in 1927–1928 in Greece causing ~ 1 million cases and ~ 1,000 deaths¹¹ with *Aedes aegypti* as the main vector. Dengue virus (DENV; *Flavivirus*, Flaviviridae) disappeared from continental Europe following the successful control of *Ae. aegypti* in 1953¹². However, this species has been detected again around the Mediterranean and Black Sea¹³. While it is described as a bad dengue vector, *Ae. albopictus* can sustain an outbreak in the absence of *Ae. aegypti*; it was involved in the first local cases of dengue in 2010, in France⁶ and Croatia¹⁴. Since 2010, autochthonous cases of chikungunya and dengue are repeatedly reported in mainland Europe: for dengue, 2013–2015¹⁵⁻¹⁷, 2018–2020¹⁸,¹⁹ in France, and for chikungunya, 2014²⁰, 2017²<fnref n="1" /> in France, and 2017 in Italy²². These European transmission episodes are usually initiated by imported cases returning from tropical regions during the summer season²³.

<footnote n="1">¹Institut Pasteur, Université Paris Cité, Insects and Infectious Diseases, 75015 Paris, France. ²Institut Pasteur, Université Paris Cité, Bioinformatics and Biostatistics Hub, 75015 Paris, France. ³Institut Pasteur, Université Paris Cité, G5 Infectious Disease Epidemiology and Analytics, 75015 Paris, France. ⁴Department of Biology and Biotechnology, University of Pavia, Pavia, Italy. ⁵Institute of Public Health, Tirana, Albania. ⁶Institut Pasteur, Université Paris Cité, Arboviruses and Insect Vectors, 75015 Paris, France. ⁷Present address: INSERM, Univ. Limoges, CHU Limoges, IRD, U1094 Neuroépidémiologie Tropicale, Institut d'Epidémiologie Et de Neurologie Tropicale, GEIST, Limoges, France. ⁸These authors jointly supervised this work: Paul Reiter and Anna-Bella Failloux. ✉email: anna-bella.failloux@pasteur.fr</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** PCoA plot obtained from the *F*<sub>ST</sub>*P* matrix using GenAlEx (a) and graphical representation of the co-ancestry percentages obtained with Distract 1.1. (b). In panel (a), highlighted populations are those showing a higher genetic correlation. In panel (b), each cluster is represented by a specific color and each bar represents a specific individual and its percentage of identity to the two clusters (K1 and K2).</figcaption></figure>

In Europe, the epidemiological landscape of mosquito-borne diseases has drastically changed these last decades mainly as a consequence of the establishment of *Ae. albopictus* in part of Europe in addition to the growing number of people travelling between arbovirus-endemic countries and Europe, and suitable climate conditions conducive to *Ae. albopictus* expansion and proliferation¹. The first record of *Ae. albopictus* in Europe was reported in Albania in 1979 probably introduced from China³. In Albania, the species is commonly found everywhere and even succeeded in reaching high altitudes up to 1200 m²⁴. Altitude can be a proxy for temperature which is a critical driver of mosquito activities requiring genetic selection and physiological adaptation of mosquitoes to the cold. Temperature determines mosquito biology, ecology, behavior, and more notably, transmission of arboviruses²⁵. In particular, temperature influences significantly adult size²⁶ and mosquito females can be larger at high altitudes²⁷. Larger females absorb more blood than smaller females²⁸, thus potentially increasing the quantity of viral particles ingested by mosquitoes and consequently, the vector competence²⁹,³⁰. Alterations in mosquito gene expression and physiology could be observed between small and large mosquitoes³¹. Thus, as mosquitoes may extend their distribution northward, arboviral diseases are expected to emerge in much more of Europe³². The critical question is how far north transmission could occur?

In this context, the Albanian *Ae. albopictus* populations are of interest as they colonize high altitudes. These altitude-adapted populations could be a good indicator of latitude-adapted mosquito populations as this invasive species moves northwards into northern Europe. Here we study the effect of temperature on transmission of CHIKV and DENV by *Ae. albopictus* collected at different altitudes in Dajti mountain in Central Albania. We showed that *Ae. albopictus* can transmit CHIKV with similar viral loads in mosquito saliva independently of the incubation temperature tested. In contrast, the transmission of DENV is only observed at 28 °C. Altogether, our study shows that CHIKV transmission can occur at 20 °C (tested as a proxy of high altitudes) increasing the risk of epidemics beyond its current range.

## Results

**Aedes albopictus** samples collected at high altitudes are genetically different from populations at low altitudes. The 11 SSR used to genotype the four populations of mosquitoes originating from different altitudes of Tirana and the Dajti mountain did not display signs of linkage disequilibrium between any pairs, suggesting that populations were independent, and their variability might well reflect genome-wide patterns. Across these loci, mosquito samples displayed similar levels of variability (Table S1) suggesting no recent bottleneck effect. However, a certain degree of inter-population differentiation was detected which increased along with altitude. Indeed, high F<sub>ST</sub>P values (0.107, 0.122 and 0.134) were estimated between the population of Tirana at 149 m and the others (542 m, 762 m and 1209 m, respectively; Table S2). The two populations at the highest altitudes (762 m and 1209 m) were more genetically similar to each other (F<sub>ST</sub>P value of 0.022) than with others. This differentiation pattern was confirmed from the PCoA (Fig. 1a) which separated the population 149 m from the two populations 762 m and 1209 m along the first axis (representative for 58.4% of variance). The population 542 m was at an intermediate position on the second axis (depictive for 34.3% of variance).
The Bayesian clustering analysis implemented in STRUCTURE, based on Evanno's method, indicated that two lineages (K1 and K2, Fig. 1b) represented the most parsimonious partitioning of the ancestry among individuals from the four populations (Table S3). The population 149 m was mainly structured in K1 (0.972%); starting from this altitude, we observed a gradual transitioning in genetic composition towards the second lineage or K2. Indeed, the population 542 m displayed an ancestry profile structured for 0.799% in K1 and 0.201% in K2. At higher altitudes, the populations 762 m and 1209 m presented the highest membership in K2 (0.895% and 0.910% respectively, Fig. 1b). When comparing with other neighbouring populations in Europe, the four Albanian populations were genetically different from the other populations (Fig. S1). The population from Tirana was genetically related to mosquitoes from China while mosquitoes from Brescia in the region of Lombardy in Northern Italy were genetically close to mosquitoes from La Reunion. Surprisingly, the population from Cesena in the Emilia-Romagna region in Italy was a mix of mosquitoes from Tirana and La Reunion. Mosquitoes from Greece were genetically distinct from the other populations (Fig. S1).

**Ae. albopictus** mosquitoes from high altitudes disseminate and transmit **CHIKV**. To examine whether viral dissemination and transmission varied according to the incubation temperature and the day post-infection (dpi), dissemination efficiency (DE) and transmission efficiency (TE) were estimated for each population infected with **CHIKV** (Fig. 2a,b; Table S4).

When considering DE, most populations started viral dissemination at 3 dpi though with overall low success rates except when mosquitoes were incubated at 28 °C (Fig. 2a). DEs were high whatever the population, reaching 83.33% for 149 m (7 dpi, 28 °C), 90% for 542 m (21 dpi, 20 °C or 28 °C), 100% for 762 m (21 dpi, 20 °C variable), and 95% for 1209 m (21 dpi, 20 °C or 28 °C) (Table S4). When comparing DEs obtained at the three incubation temperatures for each dpi and each population, only one comparison (population 762 m, 3 dpi; DE: 0% at 20 °C, 5.26% at 20 °C Var, and 57.89% at 28 °C) among 20 was statistically significant (after Bonferroni correction) (Fisher's exact test: p < 0.05).

For TE, most populations started viral transmission from 3 dpi when incubated at 28 °C, and 7 dpi at 20 °C variable and 20 °C (except population 542 m at 3 dpi) (Fig. 2b). TEs were lower reaching 66.66% for 149 m (7 dpi, 28 °C), 61.90% for 542 m (14 dpi, 20 °C variable), 72.22% for 762 m (7 dpi, 28 °C), and 86.67% for 1209 m (10 dpi, 20 °C) (Table S4). When comparing TEs at the three incubation temperatures for each dpi and each population, two comparisons among 20 were significant (after Bonferroni correction) (Fisher's exact test: p < 0.05): populations 149 m at 7 dpi and 762 m at 3 dpi. The efficiency of the salivary glands as barrier to viral transmission can be measured by estimating the transmission rate (TR). For each population, TRs were comparable whatever the dpi and the incubation temperature (Fisher's exact test: p > 0.05) with TR reaching 80.0% for 149 m (7 dpi, 28 °C), 66.7% for 542 m (21 dpi, 20 °C), 81.2% for 762 m (7 dpi, 28 °C), and 92.9% for 1209 m (10 dpi, 20 °C) (Table S5). When examining the extrinsic incubation period (EIP), all four populations were able to transmit from 3 dpi when incubated at 28 °C and only one population (542 m) from 3 dpi at 20 °C and 20 °C Variable.

In Fig. 2a, a sharp difference could be observed between the magnitudes of DE at 28 °C compared to that at lower temperatures, especially in the early dpi. This trend can be captured statistically by including interaction terms to the regression models, allowing the effect of incubation temperatures to vary according to mosquito populations. Figure 3 shows the predicted value of the outcome (Fig. 3a: DE; Fig. 3b: TE) on the logit scale according to dpi for the four populations. DE was favoured for populations at higher altitudes and an altitude cut-off was observed between populations 542 m and 762 m while no significant difference was identified between the two populations at highest altitudes and the two populations at lowest altitudes. The rate at which DE increased according to dpi was lower at 28 °C, reaching high values earlier after infection. For TE however, this cut-off was not found.

**No significantly different CHIKV loads in *Ae. albopictus* saliva.** To test whether the number of viral particles delivered by mosquitoes varied according to the incubation temperature and the dpi, saliva collected from each mosquito was titrated. When examining the four populations at a given dpi and incubation temperature, no significant difference of saliva viral titers was detected (Kruskal–Wallis test: p > 0.05; Fig. 4a-o). Likewise, no significant difference was detected according to the dpi when considering each population at a given incubation temperature (Table S6). The highest mean number of viral particles was detected in two mosquito saliva of the population 762 m incubated for 3 days at 28 °C: $10^{3.18\pm3.32}$ (1504 ± 2116) infectious particles.

**Ae. albopictus** only disseminates and transmits **DENV at 28 °C.** To measure the vector competence of *Ae. albopictus* to DENV, viral dissemination and transmission were estimated according to the incubation temperature and the dpi (Fig. 5a; Table S4). Viral dissemination was mainly detected at 28 °C, with DE starting from 7 dpi and increasing significantly along with dpi (Fisher's exact test: p < 0.05) (Fig. 5a). DE reached the highest values at 21 dpi (85% for 149 m, 80% for 542 m, 94.73% for 762 m, and 65% for 1209 m). When comparing DEs obtained at the three incubation temperatures for each dpi and each population, 10 comparisons among 16 were statistically significant (after Bonferroni correction) (Fisher's exact test: p < 0.05) with a viral dissemination mainly observed at 28 °C (Table S4). Consequently, viral transmission was only detected at 28 °C reaching highest TEs at 21 dpi (33.33% for 149 m, 15% for 542 m, 42.10% for 762 m) (Fig. 5b). No transmission was observed for population 1209 m at 21 dpi whatever the incubation temperature (Table S4). When comparing TE at the three incubation temperatures for each dpi and each population, only one among 8 comparisons was significantly different (after Bonferroni correction) (Fisher's exact test: p < 0.05): population 762 m at 21 dpi (Table S4). When examining TR, values were low with a maximum of 42.10% for population 762 m (21 dpi, 28 °C), meaning that more than half of mosquitoes with disseminated infection were not able to transmit DENV (Table S5). When examining the EIP, all four populations only transmitted at 28 °C: 149 m from 14 dpi,
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Dissemination and transmission efficiencies of *Ae. albopictus* populations infected with CHIKV examined at different days (3, 7, 10, 14, and 21) after incubation at different temperatures (20 °C, 20 °C variable, 28 °C). After infection, mosquito heads and saliva were titrated. (a) Dissemination efficiency corresponds to the proportion of mosquitoes with infected heads. (b) Transmission efficiency refers to the proportion of mosquitoes with infectious saliva. Error bars show the exact Binomial confidence interval (95%).</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Predicted value of Dissemination Efficiency (DE, panel a) and Transmission Efficiency (TE, panel b) according to days-post infection, modeled using logistic regression. In both panels, graphs show the expected marginal means (dots) with corresponding 95% confidence intervals (bars) on the logit scale, highlighting the interacting effect of temperature on days post-infection.</figcaption></figure>

542 m from 10 dpi, 762 m from 7 dpi, and 1209 m at 10 dpi (Fig. S2). The highest mean number of viral particles detected in mosquito saliva was $10^{2.04}$ (110) infectious particles for one mosquito of the population 1209 m (10 dpi, 28 °C) (Table S6).

Discussion

In this study, we showed that the four *Ae. albopictus* populations collected at 149 m, 542 m, 762 m, and 1209 m were able to transmit CHIKV at 20 °C as at 28 °C. Contrariwise, DENV transmission only occurs at 28 °C; we observed that transmission is less likely for the population at 1209 m altitude. These results indicate that CHIKV transmission may occur at high altitudes. So as the distribution of *Ae. albopictus* is expanding, the risk of CHIKV outbreaks can extend beyond South Europe and reach northern European countries.

Two waves of introduction of *Ae. albopictus* in Europe are mainly described: from China to Tirana in Albania in 1979³³, and from Japan via the United States and La Reunion Island in Italy³⁴. We showed that the two *Ae. albopictus* samples collected at high altitudes (762 m and 1209 m) in Dajti mountain in Central Albania were genetically different from the two populations at low altitudes (149 m and 542 m) and all four populations from Dajti mountain are genetically distinct from the population of Tirana³. It is believed that *Ae. albopictus* from Tirana experienced a different colonization dynamic than those from Italy, partly related to strong bottlenecks and/or differences of populations in the ability to adapt to local environmental conditions³⁵; while the species
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Viral titers in saliva of *Ae. albopictus* infected with CHIKV and examined at different days (3 (a,f,k), 7 (b,g,l), 10 (c,h,m), 14 (d,l,n), and 21(e,j,o)) after incubation at different temperatures (20 °C, 20 °C variable, 28 °C). Saliva was collected from individual females using the forced salivation technique and titrated on C6/36 *Ae. albopictus* cells. Bars indicate the mean.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Dissemination and transmission efficiencies of *Ae. albopictus* populations infected with DENV examined at different days (3, 7, 10, 14, and 21) after incubation at different temperatures (20 °C, 20 °C variable, 28 °C). After infection, mosquito heads and saliva were titrated. (a) Dissemination efficiency corresponds to the proportion of mosquitoes with infected heads. (b) Transmission efficiency refers to the proportion of mosquitoes with infectious saliva. Error bars show the exact Binomial confidence interval (95%).</figcaption></figure>
spread rapidly and widely in Italy after its first detection in 1990, *Ae. albopictus* remains confined to Albania after its introduction in 1979.

The genetic variability of founder populations is expected to be extremely low as usually, only few individuals are introduced and among them, few survive³³; this hypothesis could be raised to explain the significant genetic divergence between populations 149 m and 542 m with populations 762 m and 1209 m. The ability of *Ae. albopictus* to colonize high altitudes is presumably related to the cold hardiness of eggs. This phenotypic trait results in a modification of the egg endochorion³⁶; a physical change in the chorion structure (more compact wax layers or additional inter-membrane space between different layers) protects the embryo against freezing. Epigenetic control of the expression of these genes may allow rapid adaptation to cold; small RNAs and histone modifications are involved in diapause mechanisms³⁶⁻³⁸. To note, we used the generation F6 for genotyping which could introduce some distortions in the genetic analysis.

Vector competence data is a reliable predictor for emergence of arboviral diseases as it has been demonstrated with *Ae. albopictus* populations from southern Europe³⁹. As each pairing vector and virus genotypes lead to a specific outcome of infection which is also modulated by the temperature, all described as genotype-by-genotype-by-environment (G x G x E) interactions⁴⁰, our study will feed a risk map of CHIKV outbreaks in Europe associated with *Ae. albopictus*. We show that all four *Ae. albopictus* populations transmit CHIKV at 20 °C and 28 °C and DENV only at 28 °C; these temperatures could serve as a proxy for the transmission of CHIKV and DENV at higher latitudes in Europe. To obtain adult mosquitoes of same age and size, immature stages and adults were exposed to a controlled temperature of 24 ± 1 °C before experimental infections. Lower temperature may generate small mosquitoes ingesting a lower amount of blood and then lower number of viral particles with subsequent effects on mosquito susceptibility to virus infection⁴¹. The different mosquito genetic backgrounds may have impacted competence for arboviruses³³. Mosquitoes may regulate gene expression to adapt their phenotypes and maintain fitness in response to stressors such as temperature changes⁴²,⁴³. Temperature may alter mosquito gene expression including genes related to mosquito antiviral responses such as RNA interference⁴⁴. RNA interference can be disrupted when mosquitoes are exposed to cooler temperatures⁴⁵. Impairment of immune barriers may affect susceptibility to arboviruses including CHIKV⁴⁵,⁴⁶. Moreover, temperature may also affect mosquito microbial communities (protozoans, fungi, bacteria and viruses) which intervene in protection against pathogens²⁵. Mosquito microbiota plays an important role in defining differences in vector competence⁴⁷. At 20 °C, some commensal bacteria may promote arboviral infection by facilitating virus entry in the mosquito gut epithelium⁴⁸.

As the main vector is *Ae. albopictus*, control strategies in Europe should be more feasible than in endemic regions where multiple genotypes and vectors are involved. Likewise, innovative mosquito control strategies can also be explored such as the use of trans-infected endosymbiotic *Wolbachia* bacteria⁴⁹ and possibly lead to a successful control of *Ae. albopictus* as it has happened in the past for *Ae. aegypti*. Nevertheless, as temperature on the field may alter microbiota composition and mosquito immune responses during viral infections, the relevance and robustness of such method should be evaluating considering the different environmental temperatures⁵⁰. Thus, the Albanian experience teaches us that chikungunya is no longer a tropical disease and its expansion will follow that of its vector. Active cooperation between European countries is critical in succeeding to control *Ae. albopictus* which does not recognize borders.

## Materials and methods

**Ethic statements.** Animals (mice as blood source for mosquito rearing and rabbit for blood in mosquito experimental infections) were housed in the Institut Pasteur animal facilities (Paris) accredited by the French Ministry of Agriculture for performing experiments on live rodents. Work on animals was performed in compliance with French and European regulations on care and protection of laboratory animals (EC Directive 2010/63, French Law 2013-118, February 6th, 2013). All experiments were approved by the Ethics Committee #89 and registered under the reference APAFIS (Autorisation de Projet utilisant des Animaux à des FINS Scientifiques)#6573-2016061412077987 v2. The study was carried out in compliance with the ARRIVE guidelines.

**Mosquito rearing.** *Aedes albopictus* eggs were collected in Central Albania at Dajti mountain located a few kilometers from Tirana, the capital city of Albania. Ovitraps were placed at four altitudes (149 m, 542 m, 762 m, and 1209 m) (Fig. 6) and eggs were weekly collected in 2015 during the period of optimal activities of *Ae. albopictus*²⁷, from June to mid-August. Batches of eggs were sent to the Institut Pasteur and rearing was carried out in controlled conditions (24 ± 1 °C, 70% relative humidity, a 12:12 h (Light:Dark) photoperiod). Larvae were distributed in pans (200 larvae/pan) containing 1.5 L of dechlorinated tap water supplemented with yeast tablets. Obtained adults were placed in cages and daily provided with 10% sucrose solution. Generations F4/F5 mosquitoes were used for experimental infections and F6 for genotyping.

**Microsatellite genotyping of *Ae. albopictus* populations sampled at different altitudes.** Genomic DNA was individually extracted from 120 mosquitoes (15 males and 15 females for each population) and genotyped at 11 polymorphic simple sequence repeats (SSR) loci: Aal albmic1, 2, 3, 5, 6, 9, 11, 14, 15, 16 and 17 as described in⁵¹. PCR amplifications and fragment identifications were performed as described in⁵². Each PCR product was then diluted 1:10 in ddH₂O water and 2 µL of this dilution was added to 10 µL of a mixture of deionized formamide and GeneScan-500 ROX size standard (Applied Biosystems, CA, USA). Genotyping was processed in an ABI3730XL sequence analyser (Applied Biosystems) and data analysed using GeneScan and Genemapper software. The first step of the genetic analysis was to bin the genotyped raw data using TANDEM V. 1.09⁵³; this program helps to overcome problems related with genotyping errors. When microsatellite amplification was not successful or scoring was uncertain, re-extraction of DNA was performed.
<figure><img src="image_9.png" /></figure>

<table><thead><tr><td>Population name</td><td>Geographic coordinates (Latitude, Longitude)</td><td>Locality name</td><td>Habitat#</td></tr></thead><tbody><tr><td>1209m</td><td>19° 54' 46.197" E 41° 22' 1,110" N</td><td>Gov Residence</td><td>Forest and semi natural areas</td></tr><tr><td>762m</td><td>19°, 56' 2.357" E 41° 19' 33,738" N</td><td>Forest checkpoint (Trau)</td><td>Forest and semi natural areas</td></tr><tr><td>542m</td><td>19° 55' 56.638" E 41° 18' 39.343" N</td><td>Prisk village down</td><td>Rural area, agricultural area, arable land, permanently irrigated land</td></tr><tr><td>149m</td><td>19°49' 49.383" E 41° 18' 55.483" N</td><td>Urban Park (Dreri)</td><td>Urban, forest and semi natural areas</td></tr></tbody></table>

#, extracted from CORINE Land Cover; https://land.copernicus.eu/pan-european/corine-land-cover

<caption>**Figure 6.** *Aedes albopictus* populations collected in Albania. Mosquitoes were processed for genotyping with microsatellites and vector competence studies to DENV and CHIKV. The figure was prepared by Eng. Migel Ali, GIS expert at the Company Geo Consulting Albania (https://gc-al.com/en/staff_mali.php).</caption>

**Variability of mosquito populations and genetic structure analyses.** The genetic variability of the four sampled populations was estimated in terms of mean number of alleles ($N_m$), mean number of effective alleles ($N_e$), observed and expected heterozygosity ($H_o$ and $H_e$), unbiased expected heterozygosity ($uH_e$), inbreeding coefficient (F) and pairwise $F_{ST}$ ($F_{ST}$P) using GenAlEx 6.51⁵⁴. The statistical significance of each $F_{ST}$P value was assessed by comparison of the observed value with the values obtained in 10,000 matrix permutations. Linkage disequilibrium between pairs of loci in each sample and deviations from Hardy-Weinberg equilibrium (HWE) at each locus/sample combination were examined with GENEPOP V.4.7.5⁵⁵,⁵⁶ and the statistical significance was assessed following Bonferroni corrections. The relationships among populations were appraised using Principal Coordinate Analysis (PCoA) in GenAlEx⁵⁴.
Genetic population structure was assessed using Bayesian clustering method proposed in the software STRU CTURE v.2.3.4⁵⁷, using the admixture model and assuming independent allele frequencies. The burn-in was set to 500,000 steps and was followed by 1,000,000 Markov Chain Monte Carlo replications. All runs were repeated 20 times for each number of possible clusters (K), set between 1 to 8 (i.e. twice the number of populations). The proper number of genetic clusters was determined by plotting the log probability (L(K)) and ΔK across multiple runs as implemented in STRUCTURE HARVESTER (http://taylor0.biology.ucla.edu/structureHarvester/). The Greedy algorithm in CLUMPP V.1.1.2⁵⁸ was used to merge the 20 independent runs and the graphical representa- tion of the co-ancestry percentages obtained was plotted using DISTRUCT 1.1⁵⁹. The four populations studied were compared to populations from neighbouring countries in Europe (Italy and Greece) and from China and La Reunion, both locations being the source of colonisation of Albania and Italy by *Ae. albopictus*, respectively³³.

**Viral strains and mosquito experimental infections.** Mosquito experimental infections used CHIKV 06.21 and DENV-2 strains. CHIKV 06.21 (accession number AM258992) isolated in 2005 from a patient on La Reunion belongs to the East-Central-South African (ECSA) lineage and contains the E1-A226V mutation⁶⁰. DENV-2 strain (accession number: MK268692) was isolated in 1974 from a patient in Bangkok (Thailand)⁶¹. Both strains are our references for each virus. Viral stocks were produced after passages on C6/36 cells.

Batches of 6/10-day-old females were fed for 15 min through a pig intestine membrane covering the base of a feeder (Hemotek® membrane feeding system, UK) containing 1.4 mL of washed rabbit erythrocytes, 0.7 mL of viral suspension, and 10 mM of ATP. The titer of the blood-meal was 10⁷ FFU/mL for both viruses. Engorged females were transferred in cardboard containers and maintained with 10% sucrose in climatic chambers (KB 53, Binder, Tuttlingen, Germany) under three different incubation temperatures: (i) a constant temperature of 28 °C ± 0.1 °C, (ii) a constant temperature of 20 °C ± 0.1 °C or (iii) at temperatures displaying daily fluctuations between 17 °C ± 0.1 °C and 23 °C ± 0.1 °C (average: 20 °C ± 0.1 °C). The temperature of 28 °C represents a mean temperature in tropical regions and 20 °C, a mean temperature in France during the 2010 CHIKV outbreak⁴⁰. These temperatures correspond to those recorded during the period of the year when mosquitoes laid the most eggs in ovitraps (Fig. S3). To limit confounding factors, mosquitoes (immature stages and adults) were reared at 24 ± 1 °C before experimental infections; adult females reared at 32 °C are larger and better transmit CHIKV than females at 18 °C²⁹.

**Analysis of mosquito susceptibility.** Batches of 20 females were analysed at 3, 7, 10, 14, 21 day post-infection (dpi). After cold anaesthesia, wings and legs of each mosquito were removed and the proboscis was inserted into 20 µL tip filled with 5 µL of Fetal Bovine Serum (FBS) for saliva collection⁶². After 30–45 min, medium containing saliva was expelled into 1.5 mL tube containing 45 µL of Leibovitz L15 medium (Invitrogen, CA, USA). Head and body were separated from each mosquito and ground individually in 300 µL of L15 medium supplemented with 3% FBS.

For saliva and head samples, infectious particles were detected using focus fluorescent assay on C6/36 *Ae. albopictus* cells. Samples were inoculated onto monolayers of C6/36 *Ae. albopictus* cell culture in 96-well plates. After incubation at 28 °C for 3 and 5 days for CHIKV and DENV respectively, plates were fixed with 10% formaldehyde, washed and stained using hyper-immune ascetic fluid as primary antibody and Alexa Fluor 488 goat anti-mouse IgG as the second antibody (Life technologies, CA, USA).

Two parameters were used to describe the viral dissemination and transmission. Dissemination efficiency (DE;⁶³) gives the proportion of mosquitoes with infected head (i.e. mosquitoes able to disseminate the virus from the midgut into the mosquito general cavity). Transmission efficiency (TE;⁶³) refers to the proportion of mosquitoes with infectious saliva (i.e. mosquitoes with virus having succeeded in replicating in salivary glands and released with saliva during mosquito blood feeding). Viral loads in mosquito saliva were estimated. To measure whether the salivary glands can act as a barrier to the release of viral particles in saliva, the transmission rate (TR) was also calculated, which corresponds to the proportion of mosquitoes with virus detected in saliva among mosquitoes with disseminated virus in head. Moreover, the extrinsic incubation period (EIP) was calculated, which refers to the time from ingestion of an infectious blood meal to transmission of virus.

**Statistical analysis.** Logistic regressions were used to model DE and TE according to mosquito populations, incubation temperatures, and dpi. All three covariates (mosquito population, temperature, and dpi) were coded as categorical and for each, the lowest value served as reference level (i.e. 20 °C constant, 149 m and 3 dpi). Akaike Information Criterion was used as a guide for model selection, investigating combinations of fixed effects and interactions between the three covariates of interest. The low dimensionality of the data also allowed for conveniently visualizing the possible interactions of covariates over two-dimensions graphics (outcome according to dpi) stratified over incubation temperatures and mosquito populations.

Statistical tests were conducted using the STATA software (StataCorp LP, Texas, USA) and R 4.0.3. P-values above 0.05 were considered non-significant. If necessary, the significance level of each test was adjusted based on the number of tests run, according to the sequential method of Bonferroni⁶⁴ or using Tukey's range test for multiple comparisons in generalized linear models (R package emmeans 1.5.2–1;⁶⁵).
References

1. Barzon, L. Ongoing and emerging arbovirus threats in Europe. *J. Clin. Virol.* **107**, 38-47. https://doi.org/10.1016/j.jcv.2018.08.007 (2018).

2. Papa, A. Emerging arboviruses of medical importance in the Mediterranean region. *J. Clin. Virol.* **115**, 5-10. https://doi.org/10.1016/j.jcv.2019.03.007 (2019).

3. Adhami, J. & Reiter, P. Introduction and establishment of Aedes (Stegomyia) albopictus skuse (Diptera: Culicidae) in Albania. *J. Am. Mosq. Control Assoc.* **14**, 340-343 (1998).

4. Bonizzoni, M., Gasperi, G., Chen, X. & James, A. A. The invasive mosquito species Aedes albopictus: Current knowledge and future perspectives. *Trends Parasitol.* **29**, 460-468. https://doi.org/10.1016/j.pt.2013.07.003 (2013).

5. Rezza, G. et al. Infection with chikungunya virus in Italy: An outbreak in a temperate region. *Lancet* **370**, 1840-1846. https://doi.org/10.1016/S0140-6736(07)61779-6 (2007).

6. La Ruche, G. et al. First two autochthonous dengue virus infections in metropolitan France, September 2010. *Euro Surveill.* **15**, 19676 (2010).

7. Giron, S. et al. Vector-borne transmission of Zika virus in Europe, southern France, August 2019. *Euro Surveill.* https://doi.org/10.2807/1560-7917.ES.2019.24.45.1900655 (2019).

8. Angelini, R. et al. An outbreak of chikungunya fever in the province of Ravenna, Italy. *Euro Surveill.* **12**, E070906 (2007).

9. Tsetsarkin, K. A., Vanlandingham, D. L., McGee, C. E. & Higgs, S. A single mutation in chikungunya virus affects vector specificity and epidemic potential. *PLoS Pathog.* **3**, e201. https://doi.org/10.1371/journal.ppat.0030201 (2007).

10. Vazeille, M. et al. Two Chikungunya isolates from the outbreak of La Reunion (Indian Ocean) exhibit different patterns of infection in the mosquito Aedes albopictus. *PLoS ONE* **2**, e1168. https://doi.org/10.1371/journal.pone.0001168 (2007).

11. Anon,. The dengue epidemic in Greece. *League Nations Monthly Epidemiol. Rep.* **7**, 334 (1928).

12. Curtin, T. J. Status of Aedes aegypti in the Eastern Mediterranean. *J. Med. Entomol.* **4**, 48-50 (1967).

13. Kotsakiozi, P., Gloria-Soria, A., Schaffner, F., Robert, V. & Powell, J. R. Aedes aegypti in the Black Sea: Recent introduction or ancient remnant?. *Parasit. Vectors* **11**, 396. https://doi.org/10.1186/s13071-018-2933-2 (2018).

14. Gjenero-Margan, I. et al. Autochthonous dengue fever in Croatia, August-September 2010. *Euro Surveill.* **16**, 19805 (2011).

15. Marchand, E. et al. Autochthonous case of dengue in France, October 2013. *Euro Surveill.* **18**, 20661 (2013).

16. Giron, S. et al. Nouvelles apparitions de cas autochtones de dengue en région Provence-Alpes-Côte d'Azur, France, août-septembre 2014. *Bull. Epidémiol. Hebd.* **13-14**, 217-225 (2015).

17. Succo, T. et al. Autochthonous dengue outbreak in Nimes, South of France, July to September 2015. *Euro Surveill.* **21**, 5-11. https://doi.org/10.2807/1560-7917.ES.2016.21.30240 (2016).

18. Terrien, E. et al. Surveillance du chikungunya, de la dengue et du virus Zika en France métropolitaine, 2018. *Bull. Epidémiol. Hebd.* **19-20**, 363-373 (2019).

19. Santé Publique France. Chikungunya, dengue et zika - Données de la surveillance renforcée en France métropolitaine en 2020. (2020). <https://www.santepubliquefrance.fr/maladies-et-traumatismes/maladies-a-transmission-vectorielle/chikungunya/articles/donnees-en-france-metropolitaine/chikungunya-dengue-et-zika-donnees-de-la-surveillance-renforcee-en-france-metropolitaine-en-2019>.

20. Delisle, E. et al. Chikungunya outbreak in Montpellier, France, September to October 2014. *Euro Surveill.* **20**, 2 (2015).

21. Calba, C. et al. Preliminary report of an autochthonous chikungunya outbreak in France, July to September 2017. *Euro Surveill.* https://doi.org/10.2807/1560-7917.ES.2017.22.39.17-00647 (2017).

22. Venturi, G. et al. Detection of a chikungunya outbreak in Central Italy, August to September 2017. *Euro Surveill.* https://doi.org/10.2807/1560-7917.ES.2017.22.39.17-00646 (2017).

23. Jourdain, F. et al. From importation to autochthonous transmission: Drivers of chikungunya and dengue emergence in a temperate area. *PLoS Negl. Trop. Dis.* **14**, e0008320. https://doi.org/10.1371/journal.pntd.0008320 (2020).

24. Tisseuil, C. et al. Forecasting the spatial and seasonal dynamic of Aedes albopictus oviposition activity in Albania and Balkan countries. *PLoS Negl. Trop. Dis.* **12**, e0006236. https://doi.org/10.1371/journal.pntd.0006236 (2018).

25. Bellone, R. & Failloux, A. B. The role of temperature in shaping mosquito-borne viruses transmission. *Front. Microbiol.* **11**, 584846. https://doi.org/10.3389/fmicb.2020.584846 (2020).

26. Mohammed, A. & Chadee, D. D. Effects of different temperature regimes on the development of Aedes aegypti (L.) (Diptera: Culicidae) mosquitoes. *Acta Trop.* **119**, 38-43. https://doi.org/10.1016/j.actatropica.2011.04.004 (2011).

27. Prudhomme, J. et al. The native European Aedes geniculatus mosquito species can transmit chikungunya virus. *Emerg. Microbes Infect.* **8**, 962-972. https://doi.org/10.1080/22221751.2019.1634489 (2019).

28. Briegel, H. Fecundity, metabolism, and body size in Anopheles (Diptera: Culicidae), vectors of malaria. *J. Med. Entomol.* **27**, 839-850. https://doi.org/10.1093/jmedent/27.5.839 (1990).

29. Westbrook, C. J., Reiskind, M. H., Pesko, K. N., Greene, K. E. & Lounibos, L. P. Larval environmental temperature and the susceptibility of Aedes albopictus Skuse (Diptera: Culicidae) to Chikungunya virus. *Vector Borne Zoonotic Dis.* **10**, 241-247. https://doi.org/10.1089/vbz.2009.0035 (2010).

30. Alto, B. W., Reiskind, M. H. & Lounibos, L. P. Size alters susceptibility of vectors to dengue virus infection and dissemination. *Am. J. Trop. Med. Hyg.* **79**, 688-695 (2008).

31. Price, D. P., Schilkey, F. D., Ulanov, A. & Hansen, I. A. Small mosquitoes, large implications: crowding and starvation affects gene expression and nutrient accumulation in Aedes aegypti. *Parasit. Vectors* **8**, 252. https://doi.org/10.1186/s13071-015-0863-9 (2015).

32. Liu-Helmersson, J. et al. Climate change and aedes vectors: 21st century projections for dengue transmission in Europe. *EBio-Medicine* **7**, 267-277. https://doi.org/10.1016/j.ebiom.2016.03.046 (2016).

33. Vega-Rua, A. et al. Vector competence of Aedes albopictus populations for chikungunya virus is shaped by their demographic history. *Commun. Biol.* **3**, 326. https://doi.org/10.1038/s42003-020-1046-6 (2020).

34. Manni, M. et al. Genetic evidence for a worldwide chaotic dispersion pattern of the arbovirus vector *Aedes albopictus*. *PLoS Negl Trop Dis* **11**, e0005332. https://doi.org/10.1371/journal.pntd.0005332 (2017).

35. Pichler, V. et al. Complex interplay of evolutionary forces shaping population genomic structure of invasive Aedes albopictus in southern Europe. *PLoS Negl. Trop. Dis.* **13**, e0007554. https://doi.org/10.1371/journal.pntd.0007554 (2019).

36. Kress, A., Kuch, U., Oehlmann, J. & Muller, R. Effects of diapause and cold acclimation on egg ultrastructure: new insights into the cold hardiness mechanisms of the Asian tiger mosquito Aedes (Stegomyia) albopictus. *J. Vector Ecol.* **41**, 142-150. https://doi.org/10.1111/jvec.12206 (2016).

37. Reynolds, J. A., Bautista-Jimenez, R. & Denlinger, D. L. Changes in histone acetylation as potential mediators of pupal diapause in the flesh fly Sarcophaga bullata. *Insect Biochem. Mol. Biol.* **76**, 29-37. https://doi.org/10.1016/j.ibmb.2016.06.012 (2016).

38. Poupardin, R. et al. Early transcriptional events linked to induction of diapause revealed by RNAseq in larvae of drosophilid fly Chemyomyza costata. *BMC Genom.* **16**, 720. https://doi.org/10.1186/s12864-015-1907-4 (2015).

39. Mariconti, M. et al. Estimating the risk of arbovirus transmission in Southern Europe using vector competence data. *Sci. Rep.* **9**, 17852. https://doi.org/10.1038/s41598-019-54395-5 (2019).

40. Zouache, K. et al. Three-way interactions between mosquito population, viral strain and temperature underlying chikungunya virus transmission potential. *Proc. Biol. Sci. R. Soc.* https://doi.org/10.1098/rspb.2014.1078 (2014).
41. Alto, B. W. & Bettinardi, D. Temperature and dengue virus infection in mosquitoes: Independent effects on the immature and adult stages. *Am. J. Trop. Med. Hyg.* **88**, 497–505. https://doi.org/10.4269/ajtmh.12-0421 (2013).

42. Piglucci, M. How organisms respond to environmental changes: From phenotypes to molecules (and vice versa). *Trends Ecol. Evol.* **11**, 168–173. https://doi.org/10.1016/0169-5347(96)10008-2 (1996).

43. Sorensen, J. G. et al. Pronounced plastic and evolutionary responses to unpredictable thermal fluctuations in drosophila simulans. *Front. Genet.* **11**, 555843. https://doi.org/10.3389/fgene.2020.555843 (2020).

44. Blair, C. D. & Olson, K. E. The role of RNA interference (RNAi) in arbovirus-vector interactions. *Viruses* **7**, 820–843. https://doi.org/10.3390/v7020820 (2015).

45. Adelman, Z. N. et al. Cooler temperatures destabilize RNA interference and increase susceptibility of disease vector mosquitoes to viral infection. *PLoS Negl. Trop. Dis.* **7**, e2239. https://doi.org/10.1371/journal.pntd.0002239 (2013).

46. Kay, B. H. & Jennings, C. D. Enhancement or modulation of the vector competence of Ochlerotatus vigilax (Diptera: Culicidae) for ross river virus by temperature. *J. Med. Entomol.* **39**, 99–105. https://doi.org/10.1603/0022-2585-39.1.99 (2002).

47. Cirimotich, C. M., Ramirez, J. L. & Dimopoulos, G. Native microbiota shape insect vector competence for human pathogens. *Cell Host Microbe* **10**, 307–310. https://doi.org/10.1016/j.chom.2011.09.006 (2011).

48. Wu, P. et al. A gut commensal bacterium promotes mosquito permissiveness to arboviruses. *Cell Host Microbe* **25**, 101–112. https://doi.org/10.1016/j.chom.2018.11.004 (2019).

49. Moretti, R. et al. Combining Wolbachia-induced sterility and virus protection to fight Aedes albopictus-borne viruses. *PLoS Negl. Trop. Dis.* **12**, e0006626. https://doi.org/10.1371/journal.pntd.0006626 (2018).

50. Zouache, K., Michelland, R. J., Failloux, A. B., Grundmann, G. L. & Mavingui, P. Chikungunya virus impacts the diversity of symbiotic bacteria in mosquito vector. *Mol. Ecol.* **21**, 2297–2309. https://doi.org/10.1111/j.1365-294X.2012.05526.x (2012).

51. Houe, V., Bonizzoni, M. & Failloux, A. B. Endogenous non-retroviral elements in genomes of Aedes mosquitoes and vector competence. *Emerg. Microbes Infect.* **8**, 542–555. https://doi.org/10.1080/22221751.2019.1599302 (2019).

52. Manni, M. et al. Molecular markers for analyses of intraspecific genetic diversity in the Asian Tiger mosquito Aedes albopictus. *Parasit. Vectors* **8**, 188. https://doi.org/10.1186/s13071-015-0794-5 (2015).

53. Matschiner, M. & Salzburger, W. TANDEM: Integrating automated allele binning into genetics and genomics workflows. *Bioinformatics* **25**, 1982–1983. https://doi.org/10.1093/bioinformatics/btp303 (2009).

54. Peakall, R. & Smouse, P. E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—an update. *Bioinformatics* **28**, 2537–2539. https://doi.org/10.1093/bioinformatics/bts460 (2012).

55. Raymond, M. & Rousset, F. GENEPOP (Version 1.2): Population genetics software for exact tests and ecumenicism. *J. Hered.* **86**, 248–249. https://doi.org/10.1093/oxfordjournals.jhered.a111573 (1995).

56. Rousset, F. *genepop'007*: A complete re-implementation of the genepop software for Windows and Linux. *Mol. Ecol. Resour.* **8**, 103–106. https://doi.org/10.1111/j.1471-8286.2007.01931.x (2008).

57. Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. *Genetics* **155**, 945–959 (2000).

58. Jakobsson, M. & Rosenberg, N. A. CLUMPP: A cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. *Bioinformatics* **23**, 1801–1806. https://doi.org/10.1093/bioinformatics/btm233 (2007).

59. Rosenberg, N. DISTRUCT: A program for the graphical display of population structure. *Mol. Ecol. Notes* **4**, 137–138. https://doi.org/10.1046/j.1471-8286.2003.00566.x (2004).

60. Schuffenecker, I. et al. Genome microevolution of chikungunya viruses causing the Indian Ocean outbreak. *PLoS Med.* **3**, e263. https://doi.org/10.1371/journal.pmed.0030263 (2006).

61. Vazeille-Falcoz, M., Mousson, L., Rodhain, F., Chungue, E. & Failloux, A.-B. Variation in oral susceptibility to dengue type 2 virus of populations of Aedes aegypti from the islands of Tahiti and Moorea, French Polynesia. *Am. J. Trop. Med. Hyg.* **60**, 292–299 (1999).

62. Dubrulle, M., Mousson, L., Moutailler, S., Vazeille, M. & Failloux, A.-B. Chikungunya virus and Aedes mosquitoes: Saliva is infectious as soon as two days after oral infection. *PLoS ONE* **4**, e5895 (2009).

63. Vega-Rua, A., Zouache, K., Girod, R., Failloux, A. B. & Lourenco-de-Oliveira, R. High level of vector competence of *Aedes aegypti* and *Aedes albopictus* from ten American countries as a crucial factor in the spread of Chikungunya virus. *J Virol* **88**, 6294–6306. https://doi.org/10.1128/JVI.00370-14 (2014).

64. Holm, S. A simple sequentially rejective multiple test procedure. *Scand. J. Stat.* **6**, 65–70 (1979).

65. emmeans: Estimated Marginal Means, aka Least-Squares Means. R package version 1.5.2-1. (2020).

## Acknowledgements

We thank Charlotte Balière and Valérie Caro for technical assistance in mosquito genotyping. We also thank Parim Kadriaj and Erjon Muhaxhiri for their contribution in collecting mosquito eggs. This work was supported by the French Government's Investissement d'Avenir program, Laboratoire d'Excellence Integrative Biology of Emerging Infectious Diseases (grant ANR-10-LABX-62-IBEID and the European Commission under the 7th Framework under DengueTools Grant Agreement no. 282589.

## Author contributions

A.M. performed mosquito infection experiments, analyzed the data, and participated in the interpretation of the results. T.O. did statistical analyses. D.C., G.G.A. and A.R.M. did population genetic analyses. E.V. and S.B. collected mosquito samples. C.D. did phylogenetic analyses. G.G. performed mosquito genotyping. M.V. provided technical help. PR conceptualized the study, provided resources and edited the manuscript. ABF analyzed, visualized the data, wrote the first draft and edited the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-10977-4.

**Correspondence** and requests for materials should be addressed to A.-B.F.

**Reprints and permissions information** is available at www.nature.com/reprints.
**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022