Potential Risk of Dengue and Chikungunya Outbreaks in Northern Italy Based on a Population Model of Aedes albopictus (Diptera: Culicidae)

Giorgio Guzzetta; Fabrizio Montarsi; Frédéric Alexandre Baldacchino; Markus Metz; Gioia Capelli; Annapaola Rizzoli; Andrea Pugliese; Roberto Rosà; Piero Poletti; Stefano Merler

doi:10.1371/journal.pntd.0004762

Abstract

The rapid invasion and spread of Aedes albopictus (Skuse, 1894) within new continents and climatic ranges has created favorable conditions for the emergence of tropical arboviral diseases in the invaded areas. We used mosquito abundance data from 2014 collected across ten sites in northern Italy to calibrate a population model for Aedes albopictus and estimate the potential of imported human cases of chikungunya or dengue to generate the condition for their autochthonous transmission in the absence of control interventions. The model captured intra-year seasonality and heterogeneity across sites in mosquito abundance, based on local temperature patterns and the estimated site-specific mosquito habitat suitability. A robust negative correlation was found between the latter and local late spring precipitations, indicating a possible washout effect on larval breeding sites. The model predicts a significant risk of chikungunya outbreaks in most sites if a case is imported between the beginning of summer and up to mid-November, with an average outbreak probability between 4.9% and 25%, depending on the site. A lower risk is predicted for dengue, with an average probability between 4.2% and 10.8% for cases imported between mid-July and mid-September. This study shows the importance of an integrated entomological and medical surveillance for the evaluation of arboviral disease risk, which is a precondition for designing cost-effective vector control programs.

Author Summary

The tropical mosquito Aedes albopictus is a vector for several viruses that can be transmitted to human hosts through infectious bites. In the last 20 years, Aedes albopictus has rapidly expanded its habitat to temperate climates and it is now widespread in a large part of southern Europe and the United States. Here, it may cause outbreaks of local transmission of exotic diseases, such as chikungunya and dengue, after the arrival of infected patients through international travels. In this study, we estimate Ae. albopictus abundance in a subalpine region of Italy, by means of a mathematical model calibrated to data from adult mosquito captures throughout 2014; we then proceed to estimate the probability that an outbreak of local transmission of chikungunya and dengue will occur, following the importation from abroad of a single infectious case at different times in the year. These estimates are useful for planning preventive actions, such as programs of vector population control, as well as for increasing the preparedness of public health authorities in responding to an outbreak.

Citation: Guzzetta G, Montarsi F, Baldacchino FA, Metz M, Capelli G, Rizzoli A, et al. (2016) Potential Risk of Dengue and Chikungunya Outbreaks in Northern Italy Based on a Population Model of Aedes albopictus (Diptera: Culicidae). PLoS Negl Trop Dis 10(6): e0004762. https://doi.org/10.1371/journal.pntd.0004762

Editor: Samuel V. Scarpino, Santa Fe Institute, UNITED STATES

Received: December 2, 2015; Accepted: May 14, 2016; Published: June 15, 2016

Copyright: © 2016 Guzzetta et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information files.

Funding: This study was funded by the Autonomous Province of Trento, Research fund “Grandi Progetti 2012”, project LExEM (Laboratory of excellence for epidemiology and modelling) awarded to AP. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

The invasive mosquito species Aedes (Stegomyia) albopictus (Skuse, 1894) (syn. Stegomyia albopicta) has spread widely in Europe over the last twenty years, especially in the Mediterranean region [1, 2]. Aedes albopictus is a highly competent vector for many arboviruses including chikungunya [3, 4], dengue [3, 4] and Zika virus [5]. Therefore, its broad diffusion provides favorable conditions for the potential spread of diseases that have been until now confined to tropical regions [6], and which are continuously imported into temperate regions by infected international travelers. A major outbreak of chikungunya occurred in Italy in 2007, with over 200 cases [7], and sporadic local transmission has been recorded elsewhere in Europe [8, 9]. A small outbreak of dengue (15 cases) took place in Croatia in 2010 [10], followed by a much larger one, mainly transmitted by Aedes aegypti L., in the Portuguese island of Madeira in 2012 [11], with over 2000 cases; autochthonous dengue transmission has been repeatedly reported in France as well [12, 13]. In Italy, chikungunya and dengue are included in the mandatory medical surveillance system [14] to apply timely interventions aimed at limiting the impact of imported tropical infections [9]. In addition to medical surveillance, entomological surveillance is crucial to evaluate potential risks and identify optimal preventive and control strategies targeting the vector population. Aedes albopictus populations are highly variable across different sites, depending among other things on local temperature patterns, which strongly influence their life cycle [15], and on the amount and quality of breeding sites available for egg hatching and development of immature stages. The amount of breeding sites also depends directly on many complex factors, including precipitation patterns, air humidity, availability of natural or artificial containers, presence of shaded areas, and indirectly on land use (vegetation index, urbanization level). In particular, precipitations fill natural and artificial containers, creating potential habitats for the aquatic stages; on the other hand, excessive rainfalls may flush larval habitats reducing the overall population of adult mosquitoes [16]. Temperature and air humidity influence the evaporation of water, which is detrimental to the productivity of the breeding site [17] and may lead to complete drying out of the container. The suitability of a breeding site for oviposition may depend, among other factors, on: water surface and depth [18] which vary with evaporation and meteoric water input; nutrient concentration [19], depending on the presence, density and type of vegetation surrounding the container; and shading [20] which depends on land use and vegetation. The combination of all these dynamic variables determines a high spatial and temporal heterogeneity of mosquito population densities.

The purpose of this study is to estimate the potential risk of outbreaks of chikungunya and dengue in northern Italy related to the abundance of Ae. albopictus mosquitoes. To this aim, we estimated vector abundance by fitting a mathematical model for mosquito population dynamics to novel data obtained from entomological surveillance.

Methods

Study area and mosquito data

The study was carried out in the provinces of Belluno (46°08’27”N, 12°12’56”E) and Trento (46°04’00”N, 11°07’00”E), northern Italy. This mountainous area covers a large part of the Dolomites and the Southern Alps. More than 70% of the territory lies over 1,000 m a.s.l. and about 55% is covered by coniferous and deciduous forests [21]. The climate is temperate-oceanic with four main areas: sub- Mediterranean (close to Lake Garda with mild winters), subcontinental (the main valleys with more severe winters), continental (the alpine valleys) and alpine (the areas above the tree line) [21]. The human population as of 2014 is around 208,000 in the Belluno Province and 537,000 in the Trento Province with the majority of the population concentrated in the valley floors.

An entomological surveillance has been initiated in northeastern Italy since the introduction and establishment of Ae. albopictus in 1991 [22]. Furthermore, Aedes koreicus (Edwards, 1917), another invasive mosquito species, has been detected in the Belluno province in 2011 and in the Trento province in 2013 [23, 24]. Therefore, as part of the LExEM project (Laboratory of Excellence for Epidemiology and Modeling, www.lexem.eu), a mosquito monitoring focused on these two invasive Aedes species has been carried out in several localities of the provinces of Belluno and Trento ranging from 74m a.s.l. to 650m a.s.l. (see S1 Text). Seventy trapping locations in ten municipalities (Feltre, Povo, Riva del Garda, Santa Giustina, Strigno, Tenno, Tezze, Trento, Belluno and Rovereto) were selected by skilled entomologists in different places such as houses, schools, cemeteries, garden centers, public buildings, farms and tyre storage (Fig 1).

Download:

Fig 1. Geographic distribution of trap locations (red circles) and municipalities in the study area.

https://doi.org/10.1371/journal.pntd.0004762.g001

Mosquitoes were collected using Biogents Sentinel traps (Biogents AG, Regensburg, Germany) baited with BG lures and CO₂ from dry ice. Traps were placed in shaded positions sheltered from wind and rainfall as recommended by the manufacturer. All the collections were conducted in 2014. Each trap was set fortnightly for 24h from mid-April to the beginning of November except the trap of Rovereto which was set monthly from the end of June to the end of September. After each trapping session, mosquitoes were killed by freezing at -20°C and identified using taxonomic keys [25, 26] and confirmed by PCR if found in a location for the first time [26].

Land surface temperature

The MODIS sensors onboard the sun-synchronous polar-orbiting NASA satellites AQUA and TERRA are designed for global environmental monitoring, providing four records per day with a spatial resolution of 250 m to 1000 m. Most MODIS land products are accompanied by a detailed pixel-wise quality assessment, which facilitates further automated processing such as filtering and enhancement. Land surface temperature (LST) data for 2014 were obtained from the MODIS version 5 LST products MOD11A1 and MYD11A1 with four records per day and a spatial resolution of 1000 m. We applied outlier filtering and gap filling to these data and enhanced the spatial resolution to 250 m [27]. Here we used daily averages calculated from the four records per day of the reconstructed dataset.

Population model and calibration

We developed a mathematical model representing the dynamics of Ae. albopictus populations throughout the whole developmental cycle (eggs, larvae, pupae, female adults) within a single mosquito season (April to November). We adopted a compartmental structure similar to other models commonly used for Ae. albopictus [28–30] and other related species [31]. We used previous estimates of temperature-dependent rates of mortality, development from one stage to the next, and gonotrophic cycle [28] based on experimental data specific for Ae. albopictus [15]. The model is initialized at April 1^st 2014 with a fixed number of 10,000 mosquito eggs. The model robustness with respect to this assumption was verified by means of a sensitivity analysis. Abundance variability across sites was captured by local temperatures as well as by the estimation of a site-specific parameter summarizing the quantity and quality of breeding sites and representing the carrying capacity of the larval population. In other words, larval carrying capacity represents a measure of habitat suitability in our analysis. The process of capture of adults was modeled with a common capture rate across sites that was also estimated during the calibration process. We do not explicitly consider mosquito dispersal, because the average flight range of Ae. albopictus is much smaller (a few hundreds of meters [32]) than the scale of our risk predictions, which are given at the municipality level. However, the mosquito flight range was implicitly considered in the modeling of the capture process (see S1 Text). The number of captures estimated by the model was calibrated to experimental data by a Monte Carlo Markov Chain approach based on a Poisson likelihood combined over the capture sessions and ten study sites. All details are reported in the S1 Text.

Association between precipitation and habitat suitability

We investigated the relationship between site-specific habitat suitability (encoded by values of the larval carrying capacity estimated by the population model) and two precipitation variables: the average daily rainfall (r_A) and the number of rainy days (d_r), calculated for each study site. We sampled 10,000 sets of values for the carrying capacities of the ten sites from their posterior distributions, and computed the precipitations occurring at the corresponding study sites within a given time window. We performed a correlation test with Spearman rank method between each set of carrying capacity values and the two precipitation variables. If the correlation was statistically significant (p<0.05), the value of Spearman’s rank correlation coefficient (ρ) was stored for that test; otherwise we considered the correlation coefficient to be zero. We then computed the 95% credible interval (CI) of the distribution of 10,000 correlation coefficients and considered the correlation to be robust if both boundaries resulted to have the same sign (i.e. either strictly positive or strictly negative). This procedure was iterated over 340 different temporal windows, starting on 34 different dates (between January 1^st and August 20^th 2014, with interval of one week) and with ten different window lengths (between 4 and 13 weeks). For instance, the first temporal window of 4 weeks starts on January 1^st and ends on January 28^th while the last temporal window of ten weeks starts on August 20^th and ends on November 10^th.

Risk of arbovirus outbreaks

We assessed the potential risk of outbreaks of chikungunya and dengue caused by the introduction in the study sites of a single infectious individual, in the absence of control interventions. To this aim, we used mathematical formulae previously derived by the analysis of host-vector transmission dynamic models [33, 34]. These models assume the same SEI-SEIR epidemiological structure for both chikungunya and dengue: Ae. albopictus female adults may become infected after biting an infectious human host and start transmitting to other hosts after an incubation period and throughout the rest of their life (SEI). Human hosts may be infected by a bite of an infectious mosquito, develop symptoms after a latent period and may transmit to other mosquitoes for a given time, called infectious period (SEIR). The average number of mosquitoes infected by a single infectious human host in a population of fully susceptible mosquitoes and hosts is given by: where k is the biting rate of mosquitoes, β_V is the probability of transmission to mosquitoes per bite, τ is the infectious period of a human host, V is the density of female adult vectors, H is the density of human hosts, ω_V is the inverse of the mosquito incubation period and μ_V is the mosquito mortality rate. Similarly, the average number of hosts infected by a single infectious mosquito introduced in a population of fully susceptible mosquitoes and hosts is given by: where β_H is the probability of transmission to humans per bite.

The basic reproduction number (i.e. the average number of secondary human infections caused by the introduction of a single infectious host in a population of fully susceptible mosquitoes and hosts) can be simply computed as the product of these two quantities: When a single infectious host is introduced in a population, the transmission process is stochastic and the host can heal before transmitting the infection even for values of R₀ much above the epidemic threshold, purely because of chance. In this case, transmission dies out early without causing an epidemic outbreak. The probability that an outbreak will actually take place can thus be defined, and for a SEI-SEIR model is given by [35]:

We compute the reproduction numbers and probabilities of outbreak for both chikungunya and dengue in the ten study sites at different times of introduction of the index case. In particular, the values for the transmission probabilities, host infectious period and mosquito incubation period were taken from [28] for chikungunya, and from [36] for dengue. The mosquito incubation period and the transmission probabilities for dengue (but not for chikungunya) have been shown to depend on temperature [36], according to formulae initially estimated by Focks and colleagues on Ae. aegypti [37]. We used these formulae to compute those parameter values over time based on the daily temperature recorded at each trap locations. For both chikungunya and dengue, we assumed the same biting rate of 0.09 bites per mosquito per day estimated for the chikungunya outbreak occurred in 2007 in Italy [28]. The density of human hosts was estimated by dividing the census population of municipalities considered in the study by the corresponding urban surface. Finally, the density of female adult mosquitoes was obtained by dividing model estimates of the total population by the surface of the modeled capture area (about 56 ha, see S1 Text). Differences across sites in human density, local temperatures influencing epidemiological parameters, and the estimated vector abundance translate into a geographical heterogeneity in the risk of outbreaks.

Results

Model calibration

In Fig 2, we show the density of female adults over time as estimated by the model for the different study sites. The curves are substantially synchronized across sites due to the relative homogeneity of temperature patterns over time, although quantitative differences in average temperatures (see S1 Text) contribute to the geographical heterogeneity in mosquito abundance. The average peak density ranges from 17.0 (95% CI: 10.3–25.5) adult female mosquitoes per hectare in Strigno to 365 (95% CI: 264–562) in Riva del Garda. These estimates were quite robust with respect to the model assumption on the initial number of eggs (see S1 Text).

Download:

Fig 2. Model estimates for densities of adult females by study site; grey bars represent 95% CI calculated over 200 samples of the posterior distribution of parameters and 50 stochastic iterations.

https://doi.org/10.1371/journal.pntd.0004762.g002

Overall, the general seasonal pattern is well reproduced both qualitatively and quantitatively by the model (see S1 Text). Fig 3 shows boxplots of the distribution of errors between the observed number of captured females and corresponding model estimates (with 95% CI) for the ten study sites, normalized by the corresponding maximum number of captured females during the season to allow comparability across sites. In Table 1, we reported two metrics to assess the goodness of fit: the value of R² between predicted and observed values and the root mean squared error, normalized by the maximum number of captured females (nRMSE). Both metrics are obtained by averaging over the 10,000 stochastic simulations. Two sites, Strigno and Belluno, are characterized by a lower R² and higher nRMSE. The poorer fit at these two sites can be explained with the low estimated local mosquito abundance (see Fig 2), resulting in a stronger variability driven by stochastic effects in the capture process. Posterior means and 95% CI of estimates for the capture rate (α₀) and the site-specific larval carrying capacities (a_s), as estimated by the MCMC procedures, are reported in Table 2.

Download:

Fig 3. Boxplots of the errors between observed and simulated captures at each capture session in the ten study sites.

The errors are normalized by the maximum number of captured females throughout the season.

https://doi.org/10.1371/journal.pntd.0004762.g003

Download:

Table 1. Metrics for the goodness of fit between predicted and observed mosquito abundances.

https://doi.org/10.1371/journal.pntd.0004762.t001

Download:

Table 2. Estimates of model parameters from the MCMC calibration procedure.

https://doi.org/10.1371/journal.pntd.0004762.t002

Association between local factors, habitat suitability and transmissibility

A robust association between the average daily amount of precipitations and the estimated carrying capacities was found only in one time window, covering the period May 1^st—June 5^th, with a strong negative correlation (average Spearman’s ρ = -0.76; 95% CI: -0.69 –-0.82). When considering the number of rainy days, a negative correlation was also found for 17 time windows, covering periods between April 17^th and July 17^th (average ρ ranging between -0.70 and -0.84, depending on the time window, see S1 Text). Thus, late spring/early summer precipitations contribute negatively to local mosquito abundance and therefore to transmissibility of arboviruses. Figs 4 and 5 show the predicted peak value of R₀ in the ten study sites for chikungunya and dengue respectively, against four site-specific variables: local altitude (computed as an average between the altitudes of active traps), average daily rainfall between May 1^st and June 5^th, average temperature between April 1^st and October 31^st, and urban density. These figures summarize the effects of different local environmental factors on the geographic heterogeneity in the transmissibility of arboviral diseases. A Spearman rank correlation test between the peak value of R₀ for each disease and each variable was computed, and the results are summarized in Table 3. Lower altitudes and higher temperatures significantly increase R₀ by creating favorable conditions for the vector population. Conversely, an excess of precipitations in late spring is associated with a reduction in the larval carrying capacity and, therefore, in the amount of mosquitoes available for transmission. Finally, urban population density does not seem to have a strong effect on the transmissibility of disease, despite the fact that this variable appears directly in the denominator of R₀. This is explained by considering that, in a mountainous region such as the one under study, densely populated areas are concentrated at the bottom of the valleys, characterized by lower altitudes and warmer temperatures. Therefore, the negative effect of higher urban density on transmission is offset by the presence in the same sites of higher vector densities.

Download:

Fig 4. Scatterplots of the predicted peak of R₀ for chikungunya at different sites and corresponding local variables; dashed lines represent regression lines.

https://doi.org/10.1371/journal.pntd.0004762.g004

Download:

Fig 5. Scatterplots of the predicted peak of R₀ for dengue at different sites and corresponding local variables; dashed lines represent regression lines.

https://doi.org/10.1371/journal.pntd.0004762.g005

Download:

Table 3. Spearman’s correlation coefficients (ρ) between environmental variables and predicted peak of R₀ for chikungunya and dengue.

https://doi.org/10.1371/journal.pntd.0004762.t003

Risk of arbovirus outbreaks

We estimated the basic reproduction number of a potential chikungunya outbreak originated by a single imported case. The peak value of R₀ is highly heterogeneous across sites, ranging from 0.25 (95%CI: 0.15–0.38) in Strigno to 2.9 (95% CI: 2.0–4.2) in Feltre and occurring, depending on the site, between mid-August and mid-September. The latter compares well with the R₀ estimated for the chikungunya outbreak occurred in the Italian region of Emilia-Romagna in 2007 [28] and suggests the existence of further areas where similar outbreaks are likely to occur. Fig 6 reports the corresponding probability of observing an outbreak for different times of importation of the index case in the ten study sites. The three sites characterized by lowest estimated mosquito abundance (Strigno, Tenno and Belluno) are, according to the model, virtually safe from possible outbreaks. For all other sites, we can define a “chikungunya season” as the period over which the probability of observing an outbreak is higher than zero. The probability of an outbreak, averaged within each site’s chikungunya season, ranges from 4.9% in Tezze over a period of just below two months (end July to mid-September), to 25% in Feltre for a duration of 4.5 months covering the whole summer and up to mid-November. The peak probability within these seasons ranges from 8.6% (95% CI: 0.0–23.2%) in Tezze to 38.3% in Feltre (95% CI: 25.0–51.2%).

Download:

Fig 6. Model predictions for the probability of occurrence (average and 95% CI) of an outbreak of sustained chikungunya transmission caused by a single importation of an infected case occurred at different times of the year in the ten study sites.

https://doi.org/10.1371/journal.pntd.0004762.g006

Equivalent model predictions for dengue suggest generally lower values of R₀ than for chikungunya. This is due to lower average transmissibility (confirmed experimentally in Ae. albopictus from northern Italy [38] and France [3]) and longer incubation periods [37, 38], resulting in longer generation times, so that a higher proportion of mosquitoes die while in the exposed (non-infectious) phase. Peak values of R₀ range from 0.15 (95% CI: 0.09–0.22) in Strigno to 2.3 (95% CI: 1.7–3.3) in Riva del Garda and occur generally at the end of July, i.e. earlier than chikungunya. The anticipated peak of dengue transmissibility depends on the steady decrease of temperatures registered in 2014 after the end of July across all study sites (see S1 Text), which has reduced the estimated transmission probability and increased the incubation period [37], thereby offsetting the continued growth of vector abundance throughout August. Fig 7 shows the predicted probability of observing a dengue outbreak at different times of importation. Predictions are much more erratic over time than for chikungunya due to the strong temperature dependence assumed for dengue epidemiological parameters. Tezze adds to the list of virtually outbreak-free study sites, and the probability of outbreaks averaged over dengue seasons ranges from 4.2% in Trento, with a season length of one month across the end of July, to 10.8% in Riva del Garda over a two-month season that extends up to mid-September. The peak probability ranges from 12.1% (95% CI: 1.2–24.8%) in Trento to 27.8% (95% CI: 16.5–40.9%) in Riva del Garda.

Download:

Fig 7. Model predictions for the probability of occurrence (average and 95% CI) of an outbreak of sustained dengue transmission caused by a single importation of an infected case occurred at different times of the year in the ten study sites.

https://doi.org/10.1371/journal.pntd.0004762.g007

Spatializing the risk of outbreaks

A spatial interpolation of R₀ is not possible, given the small number of data points (ten, i.e. one for each study site). For the same reason, a multivariate analysis to assess the contribution of environmental factors to local mosquito abundance or transmissibility would be statistically underpowered. Therefore, to provide spatial estimates of possible ranges of the outbreak risk, we applied our model to predict vector densities within the provinces of Trento and Belluno, based on local temperature records from satellite data and considering spatially constant values of the larval carrying capacity (see S1 Text). Afterwards, we computed the corresponding outbreak risks using population density maps derived from high-resolution census data. In Fig 8, we show the peak values of the outbreak risks for both infections in the worst-case scenario, i.e. using the maximum value of the larval carrying capacity estimated among our ten study sites. Fig 8 shows that areas with highest risk are rural villages placed at the bottom of the main valleys, where vector populations are predicted to be most abundant while human density is relatively lower. Minor valleys may be considered potentially free from the risk of outbreaks because of their lower mosquito abundance, which is in turn due to lower average temperatures. In the worst-case scenario, the fraction of the population living in areas with peak outbreak risk higher than zero is up to 74% of the total for chikungunya and 50% for dengue. However, the importation of a case in the study area at the exact time of the peak has a limited probability of generating an outbreak: on average 33% for chikungunya and 16% for dengue. The fraction of the population living in areas with a risk of outbreak higher than 50% after importation of a case at the time of highest risk is 25% for chikungunya and 6% for dengue. These results are in agreement with the specific findings obtained for the 10 study sites where entomological surveillance was carried out.

Download:

Fig 8. Peak values of the predicted outbreak risk over the study region under the maximum value of the larval carrying capacity estimated by our model in the 10 capture sites (a = 95).

A) chikungunya; B) dengue.

https://doi.org/10.1371/journal.pntd.0004762.g008

Discussion

The diffusion of Ae. albopictus in temperate climates has created the conditions for local transmission of typically tropical pathogens, such as dengue and chikungunya viruses. These pathogens are continually imported into Europe through returning infected travelers from endemic areas, and may occasionally cause outbreaks of autochthonous transmission [7, 8, 10, 11]. In this work, we explore the probability of occurrence of dengue and chikungunya outbreaks in ten different sites from northern Italy, using a mosquito population model calibrated to capture data collected throughout 2014. The model was able to reproduce seasonal patterns in mosquito abundance, as well as variability across study sites. This geographic heterogeneity was represented in the model by local temperatures driving mosquito life cycles and by site-specific habitat suitability, summarized by a parameter encoding the larval carrying capacity. Increased habitat suitability was strongly associated with reduced precipitations in late spring and early summer. This negative correlation can be interpreted as a washout effect of rainfalls on breeding sites, previously observed in field experiments [39] and entomological observations [40], although other studies did not find any association between precipitations and Ae. albopictus abundance [41]. Given the strong inter-annual variability of precipitation patterns, our result will need to be validated by further studies.

Our analysis also provides estimates on the transmission potential and on the probability of outbreaks of chikungunya and dengue in the absence of control interventions, for different dates of importation of the index case. Based on model estimates of site-specific mosquito abundance and previous estimates for parameters driving host-vector transmission dynamics [28, 36], results suggest that six sites out of ten (Feltre, Povo, Riva del Garda, Santa Giustina, Trento and Rovereto) are affected by a high and sustained risk of chikungunya outbreak, as well as by a comparatively lower risk of dengue that extends over a shorter period. One study site (Tezze) had a negligible risk of a dengue outbreak but a significant risk of a chikungunya outbreak; the remaining three sites (Belluno, Strigno, and Tenno) were characterized by lowest observed abundances and had a negligible risk of outbreak for both diseases independently of the time of importation. Based on risk maps of potential outbreak in the study region, results suggest that the density of mosquitoes is expected to be higher in the main valleys. Within these areas, medium-sized towns (5,000–20,000 inhabitants) may be the local hotspots of an outbreak risk due to their lower urban densities compared to larger towns and cities. The limitation of these risk maps is the use of spatially-constant values for the parameter encoding the local availability and productivity of breeding sites. A punctual estimation of the actual outbreak risk in a given location would require additional entomological surveillance data similar to those collected from our 10 study sites. However, these maps are a useful tool for providing scenario predictions of the potential risk for Chikungunya and Dengue in our study region.

For what concerns the temporal variability of the risk of outbreak, it is interesting to compare our predictions with the observed distribution of the dates of importation in Italy in 2006–2014 ([42] and C. Rizzo, personal communication), which shows a clear seasonal pattern. For chikungunya, imported cases occur with equal probability throughout the year, except between June and September, which account for over 60% of imported cases overall. The peak of importations overlaps significantly with the peak of the outbreak risk for chikungunya (August and September, see Fig 6). Similarly, the rate of imported cases for dengue is about constant throughout the year, except for the months of August and September, representing alone about one sixth and one fourth of all imported cases respectively. In this case the peak of importations occurs slightly late with respect to the peak in outbreak risk (mid-July to mid-August, Fig 7), meaning that only a small fraction of imported cases occur when the risk of dengue is sufficiently high for an outbreak. This may explain why, despite the much higher number of imported cases of dengue (502 between 2006 and 2014) relatively to chikungunya (95 in the same period), Italy has already experienced a large chikungunya outbreak but only sporadic cases of autochthonous dengue transmission. However, these conclusions need to be taken with caution: the timing of importation of cases depends on international travels and on the timing of epidemics in countries from which the infection is imported. In addition, model predictions are strongly dependent on the recorded temperature patterns in 2014 in terms of the estimated mosquito abundance and, for dengue, also in terms of epidemiological parameters. Yearly variations in temperature records over time may provide different results for both the estimated risk and its timing. In the S1 Text, we propose a sensitivity analysis where a constant temperature variation between -2°C and +2°C is applied to the observed daily temperature. The analysis shows that model predictions are robust in terms of sites considered at risk of an outbreak, but the probability of outbreak and season length may vary considerably for sustained daily temperature variations higher than ±0.5°C.

Another key source of uncertainty for model estimates is the value of the biting rate. We assumed it to be equal to that estimated during the 2007 chikungunya outbreak [28] and performed a sensitivity analysis using boundaries of the estimated credible interval in the same study [28]. Resulting predictions range from complete absence of outbreak risk for both chikungunya and dengue in all sites, to a doubling of estimates of the peak probability of outbreak for both diseases with respect to the main analysis. However, even in the high biting rate scenario, sites with lowest abundances maintained a very low or negligible risk of chikungunya or dengue outbreaks. Full results of the sensitivity analysis are reported in the S1 Text.

In temperate climates, the population dynamics of Ae. albopictus during winters (overwintering) remains poorly understood [43]. Shortening photoperiods (duration of daylight) induce diapause in female adults, leading to generation of eggs that are particularly resistant to desiccation and cold [44]. It is thought that the hatching of diapausing eggs at low temperatures is slowed down, until warmer temperatures in spring allow survival of the larvae. However, other mechanisms have been shown to contribute to overwintering of Ae. albopictus, such as the availability of refuges with warmer microclimates for both adults and eggs [45], persisting oviposition at very low temperatures [43] and mosquito dispersal [46]. Lack of data on the abundance of vectors during winter in our study area prevents us from trying to reproduce the mechanism of persistence; therefore, the dynamics of mosquito abundance is simulated only within a single season. Although the population model for Ae. albopictus was able to reproduce well the mosquito capture data from the ten sites, further specifications on some of the biological modifications of the mosquito life cycle at low temperatures and short photoperiods may improve model accuracy. For example, the model has a general tendency to underestimate the sharp drop in the observed number of captured females between the end of October and the start of November. One reason for this may be that the model neglects the reduced developmental rates [47] and increased durations of gonotrophic cycles [48] caused by diapause. In addition, we may be underestimating adult mortality rates at temperatures below 10°C [49], since our mortality curves are based on data collected at a minimum temperature of 15°C [15].

Finally, a general limitation of population dynamics models applied to arbovirus vectors is the lack of biological data for the endemic mosquito strain under local, possibly natural, environmental conditions [50]. This limitation applies also to this work, where parameters for the mosquito life cycle were estimated from laboratory experiments on Ae. albopictus from Reunion [15]. Collection of analogous data from strains adapted to temperate areas is therefore needed to improve the accuracy of model estimates on vector abundance. As for epidemiological parameters, we used estimates from the 2007 chikungunya outbreak in Italy [28] and from the 2014 dengue outbreak in Madeira [36]. Although the latter epidemics was spread by Ae. aegypti, there are no available estimates for actual outbreaks caused by Ae. albopictus. We preferred to avoid using data from laboratory experiments on Ae. albopictus [3], given the uncertainties in extending these findings to natural conditions. The competence of Ae. albopictus for transmitting dengue is generally considered lower than that of Ae. aegypti, the dominant vector of these diseases in endemic areas [51]. However, recent experimental results on mosquitoes from southern France [3] found a transmission efficiency for dengue of up to 16–28% in Ae. albopictus and up to 15–18% in Ae. aegypti; however the shorter extrinsic incubation period in the latter species may compensate the lower transmission efficiency in terms of epidemic potential [3]. The high vectorial competence of Ae. albopictus from temperate climates suggests the possibility that a rapid evolutionary adaptation towards a higher transmissibility is ongoing [3,4]; in this case the risk of arboviral diseases would be higher than estimated. Furthermore, the Madeira outbreak demonstrated that Ae. aegypti poses an additional serious threat to public health [36] in Europe. Aedes aegypti has recently re-established also on the European mainland (namely in territories surrounding the Black Sea) after an absence of several decades [52, 53]. International travel of persons and goods increases episodes of introduction of Ae. aegypti in other areas, as documented in the Netherlands [52]. Entomological surveillance and mosquito control are therefore critical to prevent its re-establishment, and the ensuing costs to health systems [44], in other European countries with a suitable climate [54].

Medical surveillance against arboviral diseases has been introduced in Italy since the chikungunya outbreak in the Emilia-Romagna region in 2007 [7]. When a case is diagnosed during the mosquito season, a number of preventive measures are applied, including quarantine of the infected individuals and vector control interventions (e.g. adulticide spraying in combination with larvicide treatment of public catch basins) in the neighborhood of the imported case’s residence. Furthermore, many local authorities in areas heavily infested by Ae. albopictus have introduced routine vector control to reduce the transmission risk even in the absence of imported cases. Thanks to preventive and reactive interventions, autochthonous transmission of arboviral diseases has not been observed in Italy since the 2007 outbreak, despite the high risk estimated by our model and the high importation rates (48 chikungunya and 78 dengue cases between May and October 2014 alone). This public health success highlights once again the importance of surveillance and preventive actions against arboviral diseases in temperate climates.

There are other arbovirus infections for which Ae. albopictus is competent and which may be of great public health interest. Aedes albopictus can enhance the spread of West Nile Virus [55], which is mainly transmitted by endemic mosquito species such as Culex pipiens L. [56]. Importation of Zika virus infection [57] raises concerns as the pathogen, which generally causes mild symptoms, is suspected of being associated with Guillain-Barre’ syndrome, a severe autoimmune neurologic disease [58], and with congenital microcephaly, a disabling and potentially lethal condition in newborns [59]. We did not include in our analysis the risk of outbreak for West Nile Virus because of the significant complications in the transmission dynamics [55] due to the presence of other animal hosts such as avian species and of other vector species such as Cx. pipiens. Furthermore, we did not consider the potential risk of Zika virus outbreaks, given that estimates of epidemiological parameters from actual outbreaks [59–61] are not yet available for this infection.

The continued expansion and environmental adaptation of Ae. albopictus is making chikungunya, dengue and other tropical diseases a growing threat for public health authorities in temperate climates. Expected trends in climate and land use and more frequent pathogen introductions due to worldwide intensification of travels will exacerbate this problem in the years to come. For example, the “explosive epidemic” of Zika virus ongoing in Brazil [59] poses specific concerns related to the massive inflow of international tourists expected during the upcoming Olympics in Rio de Janeiro in the second half of August. The return of infected travelers to countries in the northern hemisphere, where the mosquito season will be at its peak, may facilitate the geographical spread of Zika virus. Enhanced entomological and medical surveillance will be especially important to prevent such a scenario; population and transmission dynamic models can assist surveillance policies by quantifying the risks of arboviral infections, designing effective intervention strategies and optimizing resource allocation for the control of Ae. albopictus populations.

Supporting Information

S1 Dataset. Capture data from the ten municipalities.

For each date and municipality, the number of active traps and the total number of female Ae. albopictus captured is reported.

https://doi.org/10.1371/journal.pntd.0004762.s001

(CSV)

S1 Text. Supplementary Material.

Further details on model implementation, validation and sensitivity analyses.

https://doi.org/10.1371/journal.pntd.0004762.s002

(PDF)

Acknowledgments

We thank for mosquito collection: Daniele Arnoldi, Francesca Bussola, Charlotte Lapère, Sara Carlin, Graziana Da Rold, Marco Dal Pont, Nicola Ferro Milone, Luca Tripepi.

Author Contributions

Conceived and designed the experiments: GG GC AR RR PP SM. Performed the experiments: GG MM FM FAB. Analyzed the data: GG MM AP RR PP SM. Wrote the paper: GG FM FAB MM GC AR AP RR PP SM.

References

1. Knudsen AB, Romi R, Majori G. Occurrence and spread in Italy of Aedes albopictus, with implications for its introduction into other parts of Europe. J Am Mosq Control Assoc. 1996; 12(2 Pt 1):177–83. pmid:8827590
- View Article
- PubMed/NCBI
- Google Scholar
2. Kraemer MU, Sinka ME, Duda KA, Mylne A, Shearer FM, Brady OJ, Messina JP, Barker CM, Moore CG, Carvalho RG, Coelho GE, Van Bortel W, Hendrickx G, Schaffner F, Wint GR, Elyazar IR, Teng HJ, Hay SI. The global compendium of Aedes aegypti and Ae. albopictus occurrence. Sci Data. 2015; 2:150035. pmid:26175912
- View Article
- PubMed/NCBI
- Google Scholar
3. Vega-Rúa A, Zouache K, Caro V, Diancourt L, Delaunay P, Grandadam M, Failloux AB. High efficiency of temperate Aedes albopictus to transmit chikungunya and dengue viruses in the Southeast of France. PLoS One. 2013; 8(3):e59716. pmid:23527259
- View Article
- PubMed/NCBI
- Google Scholar
4. Vega-Rúa A, Zouache K, Girod R, Failloux AB, Lourenço-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. 2014; 88(11):6294–306. pmid:24672026
- View Article
- PubMed/NCBI
- Google Scholar
5. Wong PS, Li MZ, Chong CS, Ng LC, Tan CH. Aedes (Stegomyia) albopictus (Skuse): a potential vector of Zika virus in Singapore. PLoS Negl Trop Dis. 2013; 7(8):e2348. pmid:23936579
- View Article
- PubMed/NCBI
- Google Scholar
6. Schaffner F, Medlock JM, Van Bortel W. Public health significance of invasive mosquitoes in Europe. Clin Microbiol Infect. 2013; 19(8):685–92. pmid:23574618
- View Article
- PubMed/NCBI
- Google Scholar
7. Rezza G, Nicoletti L, Angelini R, Romi R, Finarelli AC, Panning M, Cordioli P, Fortuna C, Boros S, Magurano F, Silvi G, Angelini P, Dottori M, Ciufolini MG, Majori GC, Cassone A; CHIKV study group. Infection with chikungunya virus in Italy: an outbreak in a temperate region. Lancet. 2007; 370(9602):1840–6. pmid:18061059
- View Article
- PubMed/NCBI
- Google Scholar
8. Delisle E, Rousseau C, Broche B, Leparc-Goffart I, L'Ambert G, Cochet A, Prat C, Foulongne V, Ferre JB, Catelinois O, Flusin O, Tchernonog E, Moussion IE, Wiegandt A, Septfons A, Mendy A, Moyano MB, Laporte L, Maurel J, Jourdain F, Reynes J, Paty MC, Golliot F. chikungunya outbreak in Montpellier, France, September to October 2014. Euro Surveill. 2015; 20(17). pii: 21108. pmid:25955774
- View Article
- PubMed/NCBI
- Google Scholar
9. Tomasello D, Schlagenhauf P. chikungunya and dengue autochthonous cases in Europe, 2007–2012. Travel Med Infect Dis. 2013; 11(5):274–84. pmid:23962447
- View Article
- PubMed/NCBI
- Google Scholar
10. Gjenero-Margan I, Aleraj B, Krajcar D, Lesnikar V, Klobučar A, Pem-Novosel I, Kurečić-Filipović S, Komparak S, Martić R, Duričić S, Betica-Radić L, Okmadžić J, Vilibić-Čavlek T, Babić-Erceg A, Turković B, Avsić-Županc T, Radić I, Ljubić M, Sarac K, Benić N, Mlinarić-Galinović G. Autochthonous dengue fever in Croatia, August-September 2010. Euro Surveill. 2011; 16(9). pii: 19805. pmid:21392489
- View Article
- PubMed/NCBI
- Google Scholar
11. Wilder-Smith A, Quam M, Sessions O, Rocklov J, Liu-Helmersson J, Franco L, Khan K. The 2012 dengue outbreak in Madeira: exploring the origins. Euro Surveill. 2014; 19(8):20718. pmid:24602277
- View Article
- PubMed/NCBI
- Google Scholar
12. Gould EA, Gallian P, De Lamballerie X, Charrel RN. First cases of autochthonous dengue fever and chikungunya fever in France: from bad dream to reality! Clin Microbiol Infect. 2010; 16(12):1702–4. pmid:21040155
- View Article
- PubMed/NCBI
- Google Scholar
13. Marchand E, Prat C, Jeannin C, Lafont E, Bergmann T, Flusin O, Rizzi J, Roux N, Busso V, Deniau J, Noel H, Vaillant V, Leparc-Goffart I, Six C, Paty MC. Autochthonous case of dengue in France, October 2013. Euro Surveill. 2013; 18(50):20661. pmid:24342514
- View Article
- PubMed/NCBI
- Google Scholar
14. Ministero della Salute. Sorveglianza dei casi umani di chikungunya, Dengue, West Nile Disease ed altre arbovirosi e valutazione del rischio di trasmissione in Italia– 2015.
15. Delatte H, Gimonneau G, Triboire A, Fontenille D. Influence of temperature on immature development, survival, longevity, fecundity, and gonotrophic cycles of Aedes albopictus, vector of chikungunya and dengue in the Indian Ocean. J Med Entomol. 2009; 46(1):33–41. pmid:19198515
- View Article
- PubMed/NCBI
- Google Scholar
16. Morin CW, Comrie AC, Ernst K. Climate and dengue transmission: evidence and implications. Environ Health Perspect. 2013;121(11–12):1264–72. pmid:24058050
- View Article
- PubMed/NCBI
- Google Scholar
17. Alto BW, Juliano SA. Precipitation and Temperature Effects on Populations of Aedes albopictus (Diptera: Culicidae): Implications for Range Expansion. J Med Entomol. 2001;38(5):646–56. pmid:11580037
- View Article
- PubMed/NCBI
- Google Scholar
18. Reiskind MH, Zarrabi AA. Water Surface Area and Depth Determine Oviposition Choice in Aedes albopictus (Diptera: Culicidae). J Med Entomol. 2012;49(1):71–6. pmid:22308773
- View Article
- PubMed/NCBI
- Google Scholar
19. Trexler JD, Apperson CS, Schal C. Laboratory and Field Evaluations of Oviposition Responses of Aedes albopictus and Aedes triseriatus (Diptera: Culicidae) to Oak Leaf Infusions. J Med Entomol. 1998;35(6):967–76. pmid:9835688
- View Article
- PubMed/NCBI
- Google Scholar
20. Delatte H, Dehecq JS, Thiria J, Domerg C, Paupy C, Fontenille D. Delatte 2008 Geographic Distribution and Developmental Sites of Aedes albopictus (Diptera: Culicidae) During a Chikungunya Epidemic Event. Vector Borne Zoonotic Dis. 2008;8(1):25–34.
- View Article
- Google Scholar
21. Neteler M, David R, Rocchini D, Castellani C & Rizzoli A. Terra and Acqua satellites track tiger mosquito invasion: modelling the potential distribution of Aedes albopictus in northeastern Italy. International Journal of Health Geographics 2011; 10, 49. pmid:21812983
- View Article
- PubMed/NCBI
- Google Scholar
22. Capelli G, Drago A, Martini S, Montarsi F, Soppelsa M, Delai N, Ravagnan S, Mazzon L, Schaffner F, Mathis A, Di Luca M, Romi R, Russo F. First report in Italy of the exotic mosquito species Aedes (Finlaya) koreicus, a potential vector of arboviruses and filariae. Parasite & Vectors 2011, 4, 188.
- View Article
- Google Scholar
23. Capelli G, Drago A, Martini S, Montarsi F, Soppelsa M, Delai N, Ravagnan S, Mazzon L, Schaffner F, Mathis A, Di Luca M, Romi R, Russo F. First report in Italy of the exotic mosquito species Aedes (Finlaya) koreicus, a potential vector of arboviruses and filariae. Parasit Vectors. 2011; 4:188. pmid:21951867
- View Article
- PubMed/NCBI
- Google Scholar
24. Montarsi F, Drago A, Martini S, Calzolari M, De Filippo F, Bianchi A, Mazzucato M, Ciocchetta S, Arnoldi D, Baldacchino F, Rizzoli A, Capelli G “Current distribution of the invasive mosquito species, Aedes koreicus [Hulecoeteomyia koreica] in northern Italy” Parasites and Vectors. 2015; 8:614. pmid:26626019
- View Article
- PubMed/NCBI
- Google Scholar
25. Becker N, Petric D, Zgomba M, Boase C, Madon M, Dahl C, Kaiser A. (2010) Mosquitoes and Their Control. 2nd edition. Heidelberg, Springer.
26. Montarsi F, Martini S, Dal Pont M, Delai N, Milone NF, Mazzucato M, Soppelsa F, Cazzola L, Cazzin S, Ravagnan S, Ciocchetta S, Russo F, Capelli G. Distribution and habitat characterization of the recently introduced invasive mosquito Aedes koreicus (Hulecoeteomyia koreica), a new potential vector and pest in north-eastern Italy. Parasites & Vectors. 2013; 6, 292.
- View Article
- Google Scholar
27. Metz M, Rocchini D, Neteler M. Surface temperatures at continental scale: tracking changes with remote sensing at unprecedented detail. Remote Sensing 2014, 6, 3822–3840.
- View Article
- Google Scholar
28. Poletti P, Messeri G, Ajelli M, Vallorani R, Rizzo C, Merler S. Transmission potential of chikungunya virus and control measures: the case of Italy. PLoS One. 2011 May 3;6(5):e18860. pmid:21559329
- View Article
- PubMed/NCBI
- Google Scholar
29. Erickson RA, Presley SM, Allen LJ, Long KR, Cox SB. A stage-structured, Aedes albopictus population model. Ecological Modelling. 2010 May 10;221(9):1273–82.
- View Article
- Google Scholar
30. Tran A, L'Ambert G, Lacour G, Benoît R, Demarchi M, Cros M, Cailly P, Aubry-Kientz M, Balenghien T, Ezanno P. A rainfall-and temperature-driven abundance model for Aedes albopictus populations. International journal of environmental research and public health. 2013 Apr 26;10(5):1698–719. pmid:23624579
- View Article
- PubMed/NCBI
- Google Scholar
31. Otero M, Solari HG, Schweigmann N. A stochastic population dynamics model for Aedes aegypti: formulation and application to a city with temperate climate. Bulletin of mathematical biology. 2006 Nov 1;68(8):1945–74. pmid:16832731
- View Article
- PubMed/NCBI
- Google Scholar
32. Marini F, Caputo B, Pombi M, Tarsitani G, della Torre A. Study of Aedes albopictus dispersal in Rome, Italy, using sticky traps in mark-release-recapture experiments. Med Vet Entomol 2010; 24(4):361–8. pmid:20666995
- View Article
- PubMed/NCBI
- Google Scholar
33. Smith DL, Battle KE, Hay SI, Barker CM, Scott TW, McKenzie FE. Ross, Macdonald, and a Theory for the Dynamics and Control of Mosquito-Transmitted Pathogens. PLoS Pathog. 2012;8(4):e1002588. pmid:22496640
- View Article
- PubMed/NCBI
- Google Scholar
34. Reiner RC, Perkins TA, Barker CM, Niu T, Chaves LF, Ellis AM, et al. A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970–2010. J. R. Soc. Interface 2013; 10.
- View Article
- Google Scholar
35. Lloyd AL, Zhang J, Root AM. Stochasticity and heterogeneity in host–vector models. J R Soc Interface 2007; 4: 851–863 pmid:17580290
- View Article
- PubMed/NCBI
- Google Scholar
36. Lourenço J, Recker M. The 2012 Madeira dengue outbreak: epidemiological determinants and future epidemic potential. PLoS Negl Trop Dis. 2014; 8(8):e3083. pmid:25144749
- View Article
- PubMed/NCBI
- Google Scholar
37. Focks D, Daniels E, Haile DG, Keesling JE (1995) A simulation model of the epidemiology of urban dengue fever: literature analysis, model development, preliminary validation, and samples of simulation results. The American journal of tropical medicine and hygiene 53: 489–506. pmid:7485707
- View Article
- PubMed/NCBI
- Google Scholar
38. Talbalaghi A, Moutailler S, Vazeille M, Failloux AB. Are Aedes albopictus or other mosquito species from northern Italy competent to sustain new arboviral outbreaks? Med Vet Entomol. 2010; 24(1):83–7. pmid:20377735
- View Article
- PubMed/NCBI
- Google Scholar
39. Dieng H, Rahman GM, Abu Hassan A, Che Salmah MR, Satho T, Miake F, Boots M, Sazaly A. The effects of simulated rainfall on immature population dynamics of Aedes albopictus and female oviposition. Int J Biometeorol 2012; 56(1):113–20. pmid:21267602
- View Article
- PubMed/NCBI
- Google Scholar
40. Roiz D, Rosa R, Arnoldi D, Rizzoli A. Effects of temperature and rainfall on the activity and dynamics of host-seeking Aedes albopictus females in northern Italy. Vector Borne Zoonotic Dis. 2010; 10: 811–816. pmid:20059318
- View Article
- PubMed/NCBI
- Google Scholar
41. Toma L, Severini F, Di Luca M, Bella A, Romi R. Seasonal patterns of oviposition and egg hatching rate of Aedes albopictus in Rome. J Am Mosq Control Assoc. 2003; 19: 19–22. pmid:12674530
- View Article
- PubMed/NCBI
- Google Scholar
42. Napoli C, Salcuni P, Pompa MG, Declich S, Rizzo C. Estimated imported infections of chikungunya and Dengue in Italy, 2008 to 2011. J Travel Med 2012; 19(5):294–7. pmid:22943269
- View Article
- PubMed/NCBI
- Google Scholar
43. Romi R, Severini F, Toma L. Cold acclimation and overwintering of female Aedes albopictus in Roma. Journal of the American Mosquito Control Association. 2006 Mar;22(1):149–51. pmid:16646341
- View Article
- PubMed/NCBI
- Google Scholar
44. Hanson SM, Craig GB Jr. 1994. Cold acclimation, diapause, and geographic origin affect cold hardiness in eggs of Aedes albopictus (Diptera:Culicidae). J Med Entomol 31:192–201. pmid:8189409
- View Article
- PubMed/NCBI
- Google Scholar
45. Tsunoda T, Cuong TC, Dong TD, Yen NT, Le NH, Phong TV, Minakawa N. Winter refuge for Aedes aegypti and Ae. albopictus mosquitoes in Hanoi during Winter. PloS one. 2014 Apr 21;9(4):e95606. pmid:24752230
- View Article
- PubMed/NCBI
- Google Scholar
46. Otero M, Schweigmann N, Solari HG. A stochastic spatial dynamical model for Aedes aegypti. Bulletin of mathematical biology. 2008 Jul 1;70(5):1297–325. pmid:18293043
- View Article
- PubMed/NCBI
- Google Scholar
47. Yee DA, Juliano SA, Vamosi SM. Seasonal photoperiods alter developmental time and mass of an invasive mosquito, Aedes albopictus (Diptera: Culicidae), across its north-south range in the United States. J Med Entomol. 2012; 49(4):825–32. pmid:22897042
- View Article
- PubMed/NCBI
- Google Scholar
48. Lacour G, Vernichon F, Cadilhac N, Boyer S, Lagneau C, Hance T. When mothers anticipate: effects of the prediapause stage on embryo development time and of maternal photoperiod on eggs of a temperate and a tropical strain of Aedes albopictus (Diptera: Culicidae). J Insect Physiol. 2014; 71:87–96. pmid:25450563
- View Article
- PubMed/NCBI
- Google Scholar
49. Brady OJ, Johansson MA, Guerra CA, Bhatt S, Golding N, Pigott DM, Delatte H, Grech MG, Leisnham PT, Maciel-de-Freitas R, Styer LM, Smith DL, Scott TW, Gething PW, Hay SI. Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings. Parasit Vectors. 2013; 6:351. pmid:24330720
- View Article
- PubMed/NCBI
- Google Scholar
50. Aznar VR, Otero M, De Majo MS, Fischer S, Solari HG. Modeling the complex hatching and development of Aedes aegypti in temperate climates. Ecological modelling. 2013 Mar 24;253:44–55.
- View Article
- Google Scholar
51. Lambrechts L, Scott TW, Gubler DJ. Consequences of the Expanding Global Distribution of Aedes albopictus for Dengue Virus Transmission. PLoS Negl Trop Dis. 2010;4(5):e646. pmid:20520794
- View Article
- PubMed/NCBI
- Google Scholar
52. Medlock JM, Hansford KM, Schaffner F, Versteirt V, Hendrickx G, Zeller H, Van Bortel W. A review of the invasive mosquitoes in Europe: ecology, public health risks, and control options. Vector Borne Zoonotic Dis. 2012 12(6):435–47. pmid:22448724
- View Article
- PubMed/NCBI
- Google Scholar
53. Schaffner F, Mathis A. Dengue and dengue vectors in the WHO European region: past, present, and scenarios for the future. Lancet Infect Dis. 2014;14(12):1271–80. pmid:25172160
- View Article
- PubMed/NCBI
- Google Scholar
54. Kraemer MU, Sinka ME, Duda KA, Mylne AQ, Shearer FM, Barker CM, Moore CG, Carvalho RG, Coelho GE, Van Bortel W, Hendrickx G, Schaffner F, Elyazar IR, Teng HJ, Brady OJ, Messina JP, Pigott DM, Scott TW, Smith DL, Wint GR, Golding N, Hay SI. The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. Elife. 2015;4:e08347. pmid:26126267
- View Article
- PubMed/NCBI
- Google Scholar
55. Fortuna C, Remoli ME, Severini F, DI Luca M, Toma L, Fois F, Bucci P, Boccolini D, Romi R,Ciufolini MG. Evaluation of vector competence for West Nile virus in Italian Stegomyia albopicta (= Aedes albopictus) mosquitoes. Med Vet Entomol, 2015. [Epub ahead of print]
- View Article
- Google Scholar
56. Angelini P, Tamba M, Finarelli AC, Bellini R, Albieri A, Bonilauri P, Cavrini F, Dottori M, Gaibani P, Martini E, Mattivi A, Pierro AM, Rugna G, Sambri V, Squintani G, Macini P. West Nile virus circulation in Emilia-Romagna, Italy: the integrated surveillance system 2009. Euro Surveill. 2010; 15(16).
- View Article
- Google Scholar
57. Tappe D, Rissland J, Gabriel M, Emmerich P, Gunther S, Held G, Smola S, Schmidt-Chanasit J. First case of laboratory-confirmed Zika virus infection imported into Europe, November 2013. Euro Surveill. 2014; 19(4). pii: 20685. pmid:24507467
- View Article
- PubMed/NCBI
- Google Scholar
58. Oehler E, Watrin L, Larre P, Leparc-Goffart I, Lastere S, Valour F, Baudouin L, Mallet H, Musso D, Ghawche F. Zika virus infection complicated by Guillain-Barre syndrome—case report, French Polynesia, December 2013. Euro Surveill. 2014; 19(9). pii: 20720. pmid:24626205
- View Article
- PubMed/NCBI
- Google Scholar
59. Fauci AS, Morens DM. Zika Virus in the Americas—Yet Another Arbovirus Threat. N Engl J Med. 2016;
- View Article
- Google Scholar
60. Duffy MR, Chen TH, Hancock WT, Powers AM, Kool JL, Lanciotti RS, Pretrick M, Marfel M, Holzbauer S, Dubray C, Guillaumot L, Griggs A, Bel M, Lambert AJ, Laven J, Kosoy O, Panella A, Biggerstaff BJ, Fischer M, Hayes EB. Zika virus outbreak on Yap Island, Federated States of Micronesia. N Engl J Med. 2009; 360(24):2536–43. pmid:19516034
- View Article
- PubMed/NCBI
- Google Scholar
61. Ioos S, Mallet HP, Leparc Goffart I, Gauthier V, Cardoso T, Herida M. Current Zika virus epidemiology and recent epidemics. Med Mal Infect. 2014; 44(7):302–7. pmid:25001879
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Knudsen AB, Romi R, Majori G. Occurrence and spread in Italy of Aedes albopictus, with implications for its introduction into other parts of Europe. J Am Mosq Control Assoc. 1996; 12(2 Pt 1):177–83. pmid:8827590
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Kraemer MU, Sinka ME, Duda KA, Mylne A, Shearer FM, Brady OJ, Messina JP, Barker CM, Moore CG, Carvalho RG, Coelho GE, Van Bortel W, Hendrickx G, Schaffner F, Wint GR, Elyazar IR, Teng HJ, Hay SI. The global compendium of Aedes aegypti and Ae. albopictus occurrence. Sci Data. 2015; 2:150035. pmid:26175912
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Vega-Rúa A, Zouache K, Caro V, Diancourt L, Delaunay P, Grandadam M, Failloux AB. High efficiency of temperate Aedes albopictus to transmit chikungunya and dengue viruses in the Southeast of France. PLoS One. 2013; 8(3):e59716. pmid:23527259
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Vega-Rúa A, Zouache K, Girod R, Failloux AB, Lourenço-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. 2014; 88(11):6294–306. pmid:24672026
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Wong PS, Li MZ, Chong CS, Ng LC, Tan CH. Aedes (Stegomyia) albopictus (Skuse): a potential vector of Zika virus in Singapore. PLoS Negl Trop Dis. 2013; 7(8):e2348. pmid:23936579
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Schaffner F, Medlock JM, Van Bortel W. Public health significance of invasive mosquitoes in Europe. Clin Microbiol Infect. 2013; 19(8):685–92. pmid:23574618
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Rezza G, Nicoletti L, Angelini R, Romi R, Finarelli AC, Panning M, Cordioli P, Fortuna C, Boros S, Magurano F, Silvi G, Angelini P, Dottori M, Ciufolini MG, Majori GC, Cassone A; CHIKV study group. Infection with chikungunya virus in Italy: an outbreak in a temperate region. Lancet. 2007; 370(9602):1840–6. pmid:18061059
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Delisle E, Rousseau C, Broche B, Leparc-Goffart I, L'Ambert G, Cochet A, Prat C, Foulongne V, Ferre JB, Catelinois O, Flusin O, Tchernonog E, Moussion IE, Wiegandt A, Septfons A, Mendy A, Moyano MB, Laporte L, Maurel J, Jourdain F, Reynes J, Paty MC, Golliot F. chikungunya outbreak in Montpellier, France, September to October 2014. Euro Surveill. 2015; 20(17). pii: 21108. pmid:25955774
View Article
PubMed/NCBI
Google Scholar

[30] View Article

[31] PubMed/NCBI

[32] Google Scholar

[ref9] 9. Tomasello D, Schlagenhauf P. chikungunya and dengue autochthonous cases in Europe, 2007–2012. Travel Med Infect Dis. 2013; 11(5):274–84. pmid:23962447
View Article
PubMed/NCBI
Google Scholar

[34] View Article

[35] PubMed/NCBI

[36] Google Scholar

[ref10] 10. Gjenero-Margan I, Aleraj B, Krajcar D, Lesnikar V, Klobučar A, Pem-Novosel I, Kurečić-Filipović S, Komparak S, Martić R, Duričić S, Betica-Radić L, Okmadžić J, Vilibić-Čavlek T, Babić-Erceg A, Turković B, Avsić-Županc T, Radić I, Ljubić M, Sarac K, Benić N, Mlinarić-Galinović G. Autochthonous dengue fever in Croatia, August-September 2010. Euro Surveill. 2011; 16(9). pii: 19805. pmid:21392489
View Article
PubMed/NCBI
Google Scholar

[38] View Article

[39] PubMed/NCBI

[40] Google Scholar

[ref11] 11. Wilder-Smith A, Quam M, Sessions O, Rocklov J, Liu-Helmersson J, Franco L, Khan K. The 2012 dengue outbreak in Madeira: exploring the origins. Euro Surveill. 2014; 19(8):20718. pmid:24602277
View Article
PubMed/NCBI
Google Scholar

[42] View Article

[43] PubMed/NCBI

[44] Google Scholar

[ref12] 12. Gould EA, Gallian P, De Lamballerie X, Charrel RN. First cases of autochthonous dengue fever and chikungunya fever in France: from bad dream to reality! Clin Microbiol Infect. 2010; 16(12):1702–4. pmid:21040155
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref13] 13. Marchand E, Prat C, Jeannin C, Lafont E, Bergmann T, Flusin O, Rizzi J, Roux N, Busso V, Deniau J, Noel H, Vaillant V, Leparc-Goffart I, Six C, Paty MC. Autochthonous case of dengue in France, October 2013. Euro Surveill. 2013; 18(50):20661. pmid:24342514
View Article
PubMed/NCBI
Google Scholar

[50] View Article

[51] PubMed/NCBI

[52] Google Scholar

[ref14] 14. Ministero della Salute. Sorveglianza dei casi umani di chikungunya, Dengue, West Nile Disease ed altre arbovirosi e valutazione del rischio di trasmissione in Italia– 2015.

[ref15] 15. Delatte H, Gimonneau G, Triboire A, Fontenille D. Influence of temperature on immature development, survival, longevity, fecundity, and gonotrophic cycles of Aedes albopictus, vector of chikungunya and dengue in the Indian Ocean. J Med Entomol. 2009; 46(1):33–41. pmid:19198515
View Article
PubMed/NCBI
Google Scholar

[55] View Article

[56] PubMed/NCBI

[57] Google Scholar

[ref16] 16. Morin CW, Comrie AC, Ernst K. Climate and dengue transmission: evidence and implications. Environ Health Perspect. 2013;121(11–12):1264–72. pmid:24058050
View Article
PubMed/NCBI
Google Scholar

[59] View Article

[60] PubMed/NCBI

[61] Google Scholar

[ref17] 17. Alto BW, Juliano SA. Precipitation and Temperature Effects on Populations of Aedes albopictus (Diptera: Culicidae): Implications for Range Expansion. J Med Entomol. 2001;38(5):646–56. pmid:11580037
View Article
PubMed/NCBI
Google Scholar

[63] View Article

[64] PubMed/NCBI

[65] Google Scholar

[ref18] 18. Reiskind MH, Zarrabi AA. Water Surface Area and Depth Determine Oviposition Choice in Aedes albopictus (Diptera: Culicidae). J Med Entomol. 2012;49(1):71–6. pmid:22308773
View Article
PubMed/NCBI
Google Scholar

[67] View Article

[68] PubMed/NCBI

[69] Google Scholar

[ref19] 19. Trexler JD, Apperson CS, Schal C. Laboratory and Field Evaluations of Oviposition Responses of Aedes albopictus and Aedes triseriatus (Diptera: Culicidae) to Oak Leaf Infusions. J Med Entomol. 1998;35(6):967–76. pmid:9835688
View Article
PubMed/NCBI
Google Scholar

[71] View Article

[72] PubMed/NCBI

[73] Google Scholar

[ref20] 20. Delatte H, Dehecq JS, Thiria J, Domerg C, Paupy C, Fontenille D. Delatte 2008 Geographic Distribution and Developmental Sites of Aedes albopictus (Diptera: Culicidae) During a Chikungunya Epidemic Event. Vector Borne Zoonotic Dis. 2008;8(1):25–34.
View Article
Google Scholar

[75] View Article

[76] Google Scholar

[ref21] 21. Neteler M, David R, Rocchini D, Castellani C & Rizzoli A. Terra and Acqua satellites track tiger mosquito invasion: modelling the potential distribution of Aedes albopictus in northeastern Italy. International Journal of Health Geographics 2011; 10, 49. pmid:21812983
View Article
PubMed/NCBI
Google Scholar

[78] View Article

[79] PubMed/NCBI

[80] Google Scholar

[ref22] 22. Capelli G, Drago A, Martini S, Montarsi F, Soppelsa M, Delai N, Ravagnan S, Mazzon L, Schaffner F, Mathis A, Di Luca M, Romi R, Russo F. First report in Italy of the exotic mosquito species Aedes (Finlaya) koreicus, a potential vector of arboviruses and filariae. Parasite & Vectors 2011, 4, 188.
View Article
Google Scholar

[82] View Article

[83] Google Scholar

[ref23] 23. Capelli G, Drago A, Martini S, Montarsi F, Soppelsa M, Delai N, Ravagnan S, Mazzon L, Schaffner F, Mathis A, Di Luca M, Romi R, Russo F. First report in Italy of the exotic mosquito species Aedes (Finlaya) koreicus, a potential vector of arboviruses and filariae. Parasit Vectors. 2011; 4:188. pmid:21951867
View Article
PubMed/NCBI
Google Scholar

[85] View Article

[86] PubMed/NCBI

[87] Google Scholar

[ref24] 24. Montarsi F, Drago A, Martini S, Calzolari M, De Filippo F, Bianchi A, Mazzucato M, Ciocchetta S, Arnoldi D, Baldacchino F, Rizzoli A, Capelli G “Current distribution of the invasive mosquito species, Aedes koreicus [Hulecoeteomyia koreica] in northern Italy” Parasites and Vectors. 2015; 8:614. pmid:26626019
View Article
PubMed/NCBI
Google Scholar

[89] View Article

[90] PubMed/NCBI

[91] Google Scholar

[ref25] 25. Becker N, Petric D, Zgomba M, Boase C, Madon M, Dahl C, Kaiser A. (2010) Mosquitoes and Their Control. 2nd edition. Heidelberg, Springer.

[ref26] 26. Montarsi F, Martini S, Dal Pont M, Delai N, Milone NF, Mazzucato M, Soppelsa F, Cazzola L, Cazzin S, Ravagnan S, Ciocchetta S, Russo F, Capelli G. Distribution and habitat characterization of the recently introduced invasive mosquito Aedes koreicus (Hulecoeteomyia koreica), a new potential vector and pest in north-eastern Italy. Parasites & Vectors. 2013; 6, 292.
View Article
Google Scholar

[94] View Article

[95] Google Scholar

[ref27] 27. Metz M, Rocchini D, Neteler M. Surface temperatures at continental scale: tracking changes with remote sensing at unprecedented detail. Remote Sensing 2014, 6, 3822–3840.
View Article
Google Scholar

[97] View Article

[98] Google Scholar

[ref28] 28. Poletti P, Messeri G, Ajelli M, Vallorani R, Rizzo C, Merler S. Transmission potential of chikungunya virus and control measures: the case of Italy. PLoS One. 2011 May 3;6(5):e18860. pmid:21559329
View Article
PubMed/NCBI
Google Scholar

[100] View Article

[101] PubMed/NCBI

[102] Google Scholar

[ref29] 29. Erickson RA, Presley SM, Allen LJ, Long KR, Cox SB. A stage-structured, Aedes albopictus population model. Ecological Modelling. 2010 May 10;221(9):1273–82.
View Article
Google Scholar

[104] View Article

[105] Google Scholar

[ref30] 30. Tran A, L'Ambert G, Lacour G, Benoît R, Demarchi M, Cros M, Cailly P, Aubry-Kientz M, Balenghien T, Ezanno P. A rainfall-and temperature-driven abundance model for Aedes albopictus populations. International journal of environmental research and public health. 2013 Apr 26;10(5):1698–719. pmid:23624579
View Article
PubMed/NCBI
Google Scholar

[107] View Article

[108] PubMed/NCBI

[109] Google Scholar

[ref31] 31. Otero M, Solari HG, Schweigmann N. A stochastic population dynamics model for Aedes aegypti: formulation and application to a city with temperate climate. Bulletin of mathematical biology. 2006 Nov 1;68(8):1945–74. pmid:16832731
View Article
PubMed/NCBI
Google Scholar

[111] View Article

[112] PubMed/NCBI

[113] Google Scholar

[ref32] 32. Marini F, Caputo B, Pombi M, Tarsitani G, della Torre A. Study of Aedes albopictus dispersal in Rome, Italy, using sticky traps in mark-release-recapture experiments. Med Vet Entomol 2010; 24(4):361–8. pmid:20666995
View Article
PubMed/NCBI
Google Scholar

[115] View Article

[116] PubMed/NCBI

[117] Google Scholar

[ref33] 33. Smith DL, Battle KE, Hay SI, Barker CM, Scott TW, McKenzie FE. Ross, Macdonald, and a Theory for the Dynamics and Control of Mosquito-Transmitted Pathogens. PLoS Pathog. 2012;8(4):e1002588. pmid:22496640
View Article
PubMed/NCBI
Google Scholar

[119] View Article

[120] PubMed/NCBI

[121] Google Scholar

[ref34] 34. Reiner RC, Perkins TA, Barker CM, Niu T, Chaves LF, Ellis AM, et al. A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970–2010. J. R. Soc. Interface 2013; 10.
View Article
Google Scholar

[123] View Article

[124] Google Scholar

[ref35] 35. Lloyd AL, Zhang J, Root AM. Stochasticity and heterogeneity in host–vector models. J R Soc Interface 2007; 4: 851–863 pmid:17580290
View Article
PubMed/NCBI
Google Scholar

[126] View Article

[127] PubMed/NCBI

[128] Google Scholar

[ref36] 36. Lourenço J, Recker M. The 2012 Madeira dengue outbreak: epidemiological determinants and future epidemic potential. PLoS Negl Trop Dis. 2014; 8(8):e3083. pmid:25144749
View Article
PubMed/NCBI
Google Scholar

[130] View Article

[131] PubMed/NCBI

[132] Google Scholar

[ref37] 37. Focks D, Daniels E, Haile DG, Keesling JE (1995) A simulation model of the epidemiology of urban dengue fever: literature analysis, model development, preliminary validation, and samples of simulation results. The American journal of tropical medicine and hygiene 53: 489–506. pmid:7485707
View Article
PubMed/NCBI
Google Scholar

[134] View Article

[135] PubMed/NCBI

[136] Google Scholar

[ref38] 38. Talbalaghi A, Moutailler S, Vazeille M, Failloux AB. Are Aedes albopictus or other mosquito species from northern Italy competent to sustain new arboviral outbreaks? Med Vet Entomol. 2010; 24(1):83–7. pmid:20377735
View Article
PubMed/NCBI
Google Scholar

[138] View Article

[139] PubMed/NCBI

[140] Google Scholar

[ref39] 39. Dieng H, Rahman GM, Abu Hassan A, Che Salmah MR, Satho T, Miake F, Boots M, Sazaly A. The effects of simulated rainfall on immature population dynamics of Aedes albopictus and female oviposition. Int J Biometeorol 2012; 56(1):113–20. pmid:21267602
View Article
PubMed/NCBI
Google Scholar

[142] View Article

[143] PubMed/NCBI

[144] Google Scholar

[ref40] 40. Roiz D, Rosa R, Arnoldi D, Rizzoli A. Effects of temperature and rainfall on the activity and dynamics of host-seeking Aedes albopictus females in northern Italy. Vector Borne Zoonotic Dis. 2010; 10: 811–816. pmid:20059318
View Article
PubMed/NCBI
Google Scholar

[146] View Article

[147] PubMed/NCBI

[148] Google Scholar

[ref41] 41. Toma L, Severini F, Di Luca M, Bella A, Romi R. Seasonal patterns of oviposition and egg hatching rate of Aedes albopictus in Rome. J Am Mosq Control Assoc. 2003; 19: 19–22. pmid:12674530
View Article
PubMed/NCBI
Google Scholar

[150] View Article

[151] PubMed/NCBI

[152] Google Scholar

[ref42] 42. Napoli C, Salcuni P, Pompa MG, Declich S, Rizzo C. Estimated imported infections of chikungunya and Dengue in Italy, 2008 to 2011. J Travel Med 2012; 19(5):294–7. pmid:22943269
View Article
PubMed/NCBI
Google Scholar

[154] View Article

[155] PubMed/NCBI

[156] Google Scholar

[ref43] 43. Romi R, Severini F, Toma L. Cold acclimation and overwintering of female Aedes albopictus in Roma. Journal of the American Mosquito Control Association. 2006 Mar;22(1):149–51. pmid:16646341
View Article
PubMed/NCBI
Google Scholar

[158] View Article

[159] PubMed/NCBI

[160] Google Scholar

[ref44] 44. Hanson SM, Craig GB Jr. 1994. Cold acclimation, diapause, and geographic origin affect cold hardiness in eggs of Aedes albopictus (Diptera:Culicidae). J Med Entomol 31:192–201. pmid:8189409
View Article
PubMed/NCBI
Google Scholar

[162] View Article

[163] PubMed/NCBI

[164] Google Scholar

[ref45] 45. Tsunoda T, Cuong TC, Dong TD, Yen NT, Le NH, Phong TV, Minakawa N. Winter refuge for Aedes aegypti and Ae. albopictus mosquitoes in Hanoi during Winter. PloS one. 2014 Apr 21;9(4):e95606. pmid:24752230
View Article
PubMed/NCBI
Google Scholar

[166] View Article

[167] PubMed/NCBI

[168] Google Scholar

[ref46] 46. Otero M, Schweigmann N, Solari HG. A stochastic spatial dynamical model for Aedes aegypti. Bulletin of mathematical biology. 2008 Jul 1;70(5):1297–325. pmid:18293043
View Article
PubMed/NCBI
Google Scholar

[170] View Article

[171] PubMed/NCBI

[172] Google Scholar

[ref47] 47. Yee DA, Juliano SA, Vamosi SM. Seasonal photoperiods alter developmental time and mass of an invasive mosquito, Aedes albopictus (Diptera: Culicidae), across its north-south range in the United States. J Med Entomol. 2012; 49(4):825–32. pmid:22897042
View Article
PubMed/NCBI
Google Scholar

[174] View Article

[175] PubMed/NCBI

[176] Google Scholar

[ref48] 48. Lacour G, Vernichon F, Cadilhac N, Boyer S, Lagneau C, Hance T. When mothers anticipate: effects of the prediapause stage on embryo development time and of maternal photoperiod on eggs of a temperate and a tropical strain of Aedes albopictus (Diptera: Culicidae). J Insect Physiol. 2014; 71:87–96. pmid:25450563
View Article
PubMed/NCBI
Google Scholar

[178] View Article

[179] PubMed/NCBI

[180] Google Scholar

[ref49] 49. Brady OJ, Johansson MA, Guerra CA, Bhatt S, Golding N, Pigott DM, Delatte H, Grech MG, Leisnham PT, Maciel-de-Freitas R, Styer LM, Smith DL, Scott TW, Gething PW, Hay SI. Modelling adult Aedes aegypti and Aedes albopictus survival at different temperatures in laboratory and field settings. Parasit Vectors. 2013; 6:351. pmid:24330720
View Article
PubMed/NCBI
Google Scholar

[182] View Article

[183] PubMed/NCBI

[184] Google Scholar

[ref50] 50. Aznar VR, Otero M, De Majo MS, Fischer S, Solari HG. Modeling the complex hatching and development of Aedes aegypti in temperate climates. Ecological modelling. 2013 Mar 24;253:44–55.
View Article
Google Scholar

[186] View Article

[187] Google Scholar

[ref51] 51. Lambrechts L, Scott TW, Gubler DJ. Consequences of the Expanding Global Distribution of Aedes albopictus for Dengue Virus Transmission. PLoS Negl Trop Dis. 2010;4(5):e646. pmid:20520794
View Article
PubMed/NCBI
Google Scholar

[189] View Article

[190] PubMed/NCBI

[191] Google Scholar

[ref52] 52. Medlock JM, Hansford KM, Schaffner F, Versteirt V, Hendrickx G, Zeller H, Van Bortel W. A review of the invasive mosquitoes in Europe: ecology, public health risks, and control options. Vector Borne Zoonotic Dis. 2012 12(6):435–47. pmid:22448724
View Article
PubMed/NCBI
Google Scholar

[193] View Article

[194] PubMed/NCBI

[195] Google Scholar

[ref53] 53. Schaffner F, Mathis A. Dengue and dengue vectors in the WHO European region: past, present, and scenarios for the future. Lancet Infect Dis. 2014;14(12):1271–80. pmid:25172160
View Article
PubMed/NCBI
Google Scholar

[197] View Article

[198] PubMed/NCBI

[199] Google Scholar

[ref54] 54. Kraemer MU, Sinka ME, Duda KA, Mylne AQ, Shearer FM, Barker CM, Moore CG, Carvalho RG, Coelho GE, Van Bortel W, Hendrickx G, Schaffner F, Elyazar IR, Teng HJ, Brady OJ, Messina JP, Pigott DM, Scott TW, Smith DL, Wint GR, Golding N, Hay SI. The global distribution of the arbovirus vectors Aedes aegypti and Ae. albopictus. Elife. 2015;4:e08347. pmid:26126267
View Article
PubMed/NCBI
Google Scholar

[201] View Article

[202] PubMed/NCBI

[203] Google Scholar

[ref55] 55. Fortuna C, Remoli ME, Severini F, DI Luca M, Toma L, Fois F, Bucci P, Boccolini D, Romi R,Ciufolini MG. Evaluation of vector competence for West Nile virus in Italian Stegomyia albopicta (= Aedes albopictus) mosquitoes. Med Vet Entomol, 2015. [Epub ahead of print]
View Article
Google Scholar

[205] View Article

[206] Google Scholar

[ref56] 56. Angelini P, Tamba M, Finarelli AC, Bellini R, Albieri A, Bonilauri P, Cavrini F, Dottori M, Gaibani P, Martini E, Mattivi A, Pierro AM, Rugna G, Sambri V, Squintani G, Macini P. West Nile virus circulation in Emilia-Romagna, Italy: the integrated surveillance system 2009. Euro Surveill. 2010; 15(16).
View Article
Google Scholar

[208] View Article

[209] Google Scholar

[ref57] 57. Tappe D, Rissland J, Gabriel M, Emmerich P, Gunther S, Held G, Smola S, Schmidt-Chanasit J. First case of laboratory-confirmed Zika virus infection imported into Europe, November 2013. Euro Surveill. 2014; 19(4). pii: 20685. pmid:24507467
View Article
PubMed/NCBI
Google Scholar

[211] View Article

[212] PubMed/NCBI

[213] Google Scholar

[ref58] 58. Oehler E, Watrin L, Larre P, Leparc-Goffart I, Lastere S, Valour F, Baudouin L, Mallet H, Musso D, Ghawche F. Zika virus infection complicated by Guillain-Barre syndrome—case report, French Polynesia, December 2013. Euro Surveill. 2014; 19(9). pii: 20720. pmid:24626205
View Article
PubMed/NCBI
Google Scholar

[215] View Article

[216] PubMed/NCBI

[217] Google Scholar

[ref59] 59. Fauci AS, Morens DM. Zika Virus in the Americas—Yet Another Arbovirus Threat. N Engl J Med. 2016;
View Article
Google Scholar

[219] View Article

[220] Google Scholar

[ref60] 60. Duffy MR, Chen TH, Hancock WT, Powers AM, Kool JL, Lanciotti RS, Pretrick M, Marfel M, Holzbauer S, Dubray C, Guillaumot L, Griggs A, Bel M, Lambert AJ, Laven J, Kosoy O, Panella A, Biggerstaff BJ, Fischer M, Hayes EB. Zika virus outbreak on Yap Island, Federated States of Micronesia. N Engl J Med. 2009; 360(24):2536–43. pmid:19516034
View Article
PubMed/NCBI
Google Scholar

[222] View Article

[223] PubMed/NCBI

[224] Google Scholar

[ref61] 61. Ioos S, Mallet HP, Leparc Goffart I, Gauthier V, Cardoso T, Herida M. Current Zika virus epidemiology and recent epidemics. Med Mal Infect. 2014; 44(7):302–7. pmid:25001879
View Article
PubMed/NCBI
Google Scholar

[226] View Article

[227] PubMed/NCBI

[228] Google Scholar

Figures

Abstract

Author Summary

Introduction

Methods

Study area and mosquito data

Land surface temperature

Population model and calibration

Association between precipitation and habitat suitability

Risk of arbovirus outbreaks

Results

Model calibration

Association between local factors, habitat suitability and transmissibility

Risk of arbovirus outbreaks

Spatializing the risk of outbreaks

Discussion

Supporting Information

S1 Dataset. Capture data from the ten municipalities.

S1 Text. Supplementary Material.

Acknowledgments

Author Contributions

References