Europe PMC

This website requires cookies, and the limited processing of your personal data in order to function. By using the site you are agreeing to this as outlined in our privacy notice and cookie policy.

From importation to autochthonous transmission: Drivers of chikungunya and dengue emergence in a temperate area.

Author information

Affiliations

  • 1. Santé publique France (French National Public Health Agency), Saint-Maurice, France.
      Authors
    • Jourdain F 1, 2
    • de Valk H 1
    • Noël H 1
    • Paty MC 1
    • Desenclos JC 1
    (5 authors)
  • 2. MIVEGEC Unit, IRD 224, CNRS 5290, Univ Montpellier, Montpellier, France.
      Authors
    • Jourdain F 1, 2
    • Roiz D 2
    • Roche B 2
    (3 authors)
  • 3. Entente interdépartementale pour la démoustication du littoral méditerranéen (EID Méditerranée), Montpellier, France.
      Authors
    • L'Ambert G 3
    (1 author)
  • 4. Santé publique France (French National Public Health Agency), Marseille, France.
      Authors
    • Franke F 4
    (1 author)
  • 5. Santé publique France (French National Public Health Agency), Toulouse, France.
      Authors
    • Guinard A 5
    (1 author)

ORCIDs linked to this article

Plos Neglected Tropical Diseases, 11 May 2020, 14(5):e0008320
https://doi.org/10.1371/journal.pntd.0008320 PMID: 32392224 PMCID: PMC7266344

This article is in the Europe PMC Open access subset. Refer to the copyright information in the article for licensing details.
Free full text in Europe PMC

Abstract 

Background

The global spread of Aedes albopictus has exposed new geographical areas to the risk of dengue and chikungunya virus transmission. Several autochthonous transmission events have occurred in recent decades in Southern Europe and many indicators suggest that it will become more frequent in this region in the future. Environmental, socioeconomic and climatic factors are generally considered to trigger the emergence of these viruses. Accordingly, a greater knowledge of the determinants of this emergence in a European context is necessary to develop adapted surveillance and control strategies, and public health interventions.

Methodology/principal findings

Using French surveillance data collected from between 2010 and 2018 in areas of Southern France where Ae. albopictus is already established, we assessed factors associated with the autochthonous transmission of dengue and chikungunya. Cases leading to autochthonous transmission were compared with those without subsequent transmission using binomial regression. We identified a long reporting delay (≥ 21 days) of imported cases to local health authorities as the main driver for autochthonous transmission of dengue and chikungunya in Southern France. The presence of wooded areas around the cases' place of residence and the accumulation of heat during the season also increased the risk of autochthonous arbovirus transmission.

Conclusions

Our findings could inform policy-makers when developing strategies to the emerging threats of dengue and chikungunya in Southern Europe and can be extrapolated in this area to other viruses such as Zika and yellow fever, which share the same vector. Furthermore, our results allow a more accurate characterization of the environments most at risk, and highlight the importance of implementing surveillance systems which ensure the timely reporting and of imported cases and swift interventions.

Free full text 

Logo of plosntdLink to Publisher's site
PLoS Negl Trop Dis. 2020 May; 14(5): e0008320.
Published online 2020 May 11. https://doi.org/10.1371/journal.pntd.0008320
PMCID: PMC7266344
PMID: 32392224

From importation to autochthonous transmission: Drivers of chikungunya and dengue emergence in a temperate area

Frédéric Jourdain, Data curation, Formal analysis, Methodology, Writing – original draft,1,2,* David Roiz, Methodology, Writing – review & editing,2 Henriette de Valk, Methodology, Supervision, Writing – review & editing,1 Harold Noël, Methodology, Writing – review & editing,1 Grégory L’Ambert, Resources, Writing – review & editing,3 Florian Franke, Data curation, Resources, Writing – review & editing,4 Marie-Claire Paty, Methodology, Writing – review & editing,1 Anne Guinard, Resources, Writing – review & editing,5 Jean-Claude Desenclos, Methodology, Supervision, Writing – review & editing,1 and Benjamin Roche, Methodology, Supervision, Writing – review & editing2
Louis Lambrechts, Editor
Author information Article notes Copyright and License information Disclaimer

Associated Data

Supplementary Materials
Data Availability Statement

Abstract

Background

The global spread of Aedes albopictus has exposed new geographical areas to the risk of dengue and chikungunya virus transmission. Several autochthonous transmission events have occurred in recent decades in Southern Europe and many indicators suggest that it will become more frequent in this region in the future. Environmental, socioeconomic and climatic factors are generally considered to trigger the emergence of these viruses. Accordingly, a greater knowledge of the determinants of this emergence in a European context is necessary to develop adapted surveillance and control strategies, and public health interventions.

Methodology/Principal findings

Using French surveillance data collected from between 2010 and 2018 in areas of Southern France where Ae. albopictus is already established, we assessed factors associated with the autochthonous transmission of dengue and chikungunya. Cases leading to autochthonous transmission were compared with those without subsequent transmission using binomial regression. We identified a long reporting delay (≥ 21 days) of imported cases to local health authorities as the main driver for autochthonous transmission of dengue and chikungunya in Southern France. The presence of wooded areas around the cases’ place of residence and the accumulation of heat during the season also increased the risk of autochthonous arbovirus transmission.

Conclusions

Our findings could inform policy-makers when developing strategies to the emerging threats of dengue and chikungunya in Southern Europe and can be extrapolated in this area to other viruses such as Zika and yellow fever, which share the same vector. Furthermore, our results allow a more accurate characterization of the environments most at risk, and highlight the importance of implementing surveillance systems which ensure the timely reporting and of imported cases and swift interventions.

Author summary

The dengue, chikungunya and Zika viruses have tremendously expanded their geographic range during recent decades and are now considered emerging threats in temperate areas. The increase in international travel and trade appear to be major factors, encouraging both a circulation of these viruses on a global scale and the dispersion of one of their main vectors: Aedes albopictus, the tiger mosquito. However, these two factors cannot fully explain the risk of viral emergence in temperate areas. Climatic factors, land cover, socioeconomic characteristics and public health systems also likely impact the risk of virus transmission. In this context, we performed a statistical analysis of the different imported cases of dengue and chikungunya in Southern France, where the tiger mosquito is already established. A comparison between contexts with local virus transmission and situations without local transmission made it possible to identify the factors associated with the risk of viral circulation. We found that a delay in reporting imported cases, the presence of wooded areas around the cases’ places of residence, and the accumulation of heat during the season, all increased the risk of transmission. The identification of these conditions is therefore of primary importance both to determine the risk of the emergence of these viruses, and to improve public health policies.

Introduction

The dengue (DENV) and chikungunya (CHIKV) viruses have greatly expanded their geographic range globally in recent decades [1] and are considered emerging public health threats throughout the world, including Europe [2]. The global number of dengue infections in 2010 was estimated at 390 (284–528) million per year, 96 (67–136) million cases being clinically manifested [3]. CHIKV has been responsible for two major epidemics in recent decades. The first spread in 2004 from Eastern Africa to the Indian Ocean and to South Asia. The second occurred in the Americas, with more than 1.2 million suspected cases reported for the 2013–2014 period [4]. DENV and CHIKV are mainly transmitted between humans through the bite of Aedes aegypti and Aedes albopictus mosquitoes in urban settings, and are introduced in non-endemic countries by infected returning travellers [5,6]. Autochthonous transmission can then occur in areas where a competent vector is established and where climatic conditions are favourable for transmission. In the Mediterranean and central Europe, only Ae. albopictus is present. Its expansion is a direct consequence of the globalization of trade [7]. The continued spread of this vector through trade and the constant growth in international travel will increase the risk of exotic viruses emerging in many other European areas. Italy, France, Croatia and Spain experienced several events of autochthonous DENV and CHIKV transmission between 2010 and 2018 [8–18]. Nevertheless, the number of imported cases remains well above the number of autochthonous transmission cases [19] and, to date, there is no evidence-based explanation as to why autochthonous transmission occurs in some circumstances in Europe but not in others.

While the presence of an established vector population and virus introduction by infected travellers are necessary conditions for the emergence of these infections, they may not be sufficient for arbovirus transmission. Indeed, effective transmission is multifactorial and results from complex interactions between mosquito vectors, the human population, viral agents, their environment and climate. Genetics play an important role in fostering the transmission of some viral genotypes by locally established vector populations [21,22]. Socioeconomic and environmental factors influence the epidemiology of the disease by affecting the introduction of the virus, the contact between vectors and hosts, vector-pathogen interactions, as well as vector population distribution and dynamics [3,23–25]. Finally, public health interventions are likely to alter the dynamics of infection transmission [26].

Aedes albopictus became established in France in 2004 and has since spread throughout a large part of the country [20]. The French population can be considered fully susceptible to DENV and CHIKV infection. There is no specific antiviral drug treatment or recommended vaccine in France for DENV or CHIKV infection. Therefore, prevention and control of these infections is based on i) larval control to reduce the vector population as a preventive measure, ii) the surveillance of human imported cases, iii) the early detection of any local transmission, and iv) the implementation of proportionate vector control measures to prevent and contain autochthonous transmission [27]. A national preparedness and response plan has been implemented since 2006 [28].

Few autochthonous transmission events have occurred in France in recent years [8,13–18] but the situation continues to evolve. Aedes albopictus is still spreading across the country, leading to a greater proportion of the population exposed to DENV and CHIKV transmission risk. Other emerging viruses may also prove to be a challenge for the country’s preparedness and response systems, something already observed with the Zika virus when cases were reported in France in 2016 following the epidemic in the Americas [29]. Although substantial human and logistical resources are already mobilized every year, the challenges raised by emerging viruses will only lead to further costs and this raises the question of the future sustainability of France’s surveillance and control system [19]. This imminent situation underlines the need for a better understanding of the factors that favour DENV and CHIKV autochthonous transmission in France.

Several statistical and mathematical models have been used to identify the determinants of the distribution and abundance of Ae. albopictus, as well as the associated transmission risk of DENV and CHIKV [30–34]. Whereas contexts where transmission is high have been widely investigated (especially for DENV), the number of studies in places with sporadic and limited transmission remain scarce [35–37].

The present study aimed to identify and quantify the relative importance of the factors associated with DENV and CHIKV autochthonous transmission events in mainland France following the introduction of a viremic traveller in areas where Ae. albopictus is already established. We based our analysis on the enhanced surveillance of imported cases and performed an in-depth study of the different transmission events identified in France during recent years by the national arbovirosis surveillance system. This work was carried out within the context of improving the country’s preparedness and response system for arboviral risks in temperate areas.

Materials and methods

Binomial regression was used to compare cases leading to autochthonous transmission with those which had no subsequent transmission.

Case definitions

The different case definitions adopted in the study are indicated in Table 1. An autochthonous transmission event was defined as the occurrence of at least one autochthonous case in the study area and during the study period. The duration of viremia was fixed at 10 days: corresponding to a period between two days before the date of symptom onset and 7 days after this date [43,44]. Only cases present in an area colonized by Ae. albopictus during their viremia were included.

Table 1

Case definitions for dengue and chikungunya virus infection.
DefinitionDengueChikungunya
Suspected caseFever of 38.5°C or higher and at least one of the following symptoms not explained by other medical conditions: headache, back pain, retro-orbital pain, myalgia, and arthralgia.
Probable caseSuspected case with positive IgM antibodies in a single sample
Confirmed caseSuspected case with at least one of the following biological results:

- Positive RT-PCR


- seroconversion


- positive NS1 test


- 4-fold increase in IgG antibodies

Suspected case with at least one of the following biological results:

- Positive RT-PCR


- seroconversion


- 4-fold increase in IgG antibodies

Imported caseCase with travel history in the 15 days before the onset of symptoms in an area known for DENV or CHIKV circulation
Autochthonous caseCase without travel history in the 15 days before the onset of symptoms in an area known for DENV or CHIKV circulation

Study area and period

The study area comprises the following five French administrative districts (départements in French) (with European NUTS 3 statistical classification [38]) along the Mediterranean coast colonized by Ae. albopictus: Alpes-Maritimes, Var, Bouches-du-Rhône, Gard and Hérault. It covers a total area of 27 436 km2 with a population of 6.1 million people at the end of 2018 [39].

This area has the oldest recorded Ae. albopictus establishment in France, with progressive colonization from East to West since 2004 (Fig 1). In 2019, the invasive process is still ongoing in the western and northern parts of the study area. The area is characterized by a typically Mediterranean climate in the coastal region, with dry summers, mild winters and irregular rainfall concentrated mainly in autumn, and potential spring downpours [40].

An external file that holds a picture, illustration, etc. Object name is pntd.0008320.g001.jpg
Spatial and temporal distribution of Aedes albopictus in the study area.

Source of data: national surveillance of invasive mosquitoes, EID-Méditerranée.

The study period included the season from 1 May through 30 November—representing the seasonal activity of Ae. albopictus in the South of France [41]–for the years 2010, 2013, 2014, 2015, 2017 and 2018. The three years not included were due to the lack of epidemiological (2012) and entomological data (2011, 2016). Our study period included all the years when an autochthonous event was observed (Table 2).

Table 2

Dengue and chikungunya autochthonous transmission events in France between 2010 and 2018.
YearLocalityDépartementVirusNumber of autochthonous casesIdentification of the imported source caseRef.
2010NiceAlpes-MaritimesDengue 12Yes[18]
2010FréjusVarChikungunya2Yes[17]
2011Absence of autochthonous transmission event
2012Absence of autochthonous transmission event
2013VenellesBouches-du-RhôneDengue 21Yes[16]
2014AubagneBouches-du-RhôneDengue 22Yes[42]
2014ToulonVarDengue 11No[42]
2014ToulonVarDengue 21No[42]
2014MontpellierHéraultChikungunya12Yes[15]
2015NîmesGardDengue 17Yes[14]
2016Absence of autochthonous transmission event
2017Le Cannet-des-maures (*)VarChikungunya11Yes[13]
2017Taradeau (*)VarChikungunya6Yes[13]
2018Saint-Laurent-du-VarAlpes-MaritimesDengue 25No[8]
2018NîmesGardDengue 11Yes[8]
2018ClapiersHéraultDengue 12No[8]

(*) These transmission foci were geographically separate but an epidemiological link was established.

Data sources

Data sources which had both with nationwide coverage and regular updates were preferred. The different data sources are summarized in Table 3 and detailed in the following sections.

Table 3

Data sources.
Data typeSource
EpidemiologicalData from the national arbovirus surveillance, Santé publique France (French Public Health Agency) https://www.santepubliquefrance.fr/
Vector control interventionsEntente Interdépartementale pour la démoustication du littoral méditerranéen—Interdepartmental Agency for mosquito control on the Mediterranean coast, (EID Méditerranée) http://www.eid-med.org/
Rainfall and temperatureFrench meteorological agency. Météo France http://www.meteofrance.com/accueil
NDVILandsat 7 and Landsat 8 data processed by the THEIA Land Data Centre (surface reflectance corrected for atmospheric effects and cloud cover level: Level 2A) https://theia.cnes.fr/ based on images acquired by the United States Geological Survey (USGS)
Land coverLand cover map produced by the Center for the Study of the Biosphere from Space (CESBIO), THEIA Land Data Centre http://www.cesbio.ups-tlse.fr/index_us.htm
Housing and vegetationBD TOPO, Institut national de l’information géographique et forestière http://www.ign.fr/
DemographicFrench National Institute of Statistics and Economic Studies (INSEE) https://www.insee.fr/en/accueil
SocioeconomicFrench National Institute of Statistics and Economic Studies (INSEE) https://www.insee.fr/en/accueil
Social Deprivation indexFrench National Institute for Health and Medical Research (Inserm) https://geo.data.gouv.fr/fr/datasets/9c6009a2bb10c4d69a15d399def4770b038be18a

Epidemiological data

Dengue and chikungunya case data were obtained through the national surveillance system for arbovirosis, which comprises active human surveillance based on the reporting of suspected cases of dengue and/or chikungunya to public health authorities, followed by timely biological diagnosis. The human surveillance system was previously described by Paty and coll. [19]. Only confirmed or probable imported cases were included in the present analysis. We used the earliest date of presence of the case during viremia (EDP) in the study area to calculate the different delays (i.e., sampling, reporting, etc.). EDP could be either the date of symptom onset or the return date from travel (in the case of symptom onset outside the study area). The sampling delay was defined as the period of time between EDP and the first date of blood sampling for biological testing. The reporting delay was defined as the period of time between EDP and the moment health authorities received report of the case. The exposure duration was defined as the time interval between EDP and the end of viremia of the imported case. Imported cases arriving in the study area 7 days or more after the date of symptom onset were excluded since they were considered to be no longer viremic.

Vector control interventions data

Each observation was completed with data derived from details regarding entomological operations (entomological survey, vector control). Data were obtained from the French mosquito control agency, EID-Méditerranée. The intervention delay was defined as the period of time between EDP and the first intervention focusing on vectors, whether an entomological survey or vector control. Entomological surveys aim to identify the presence of the vector–pre-imaginal stages or adults—while control measures are taken when the vector’s presence is confirmed. Vector control is implemented according to national guidelines and depends on the local entomological, epidemiological and environmental context. The intervention may consist of one or more of the following components: source reduction of breeding sites, application of larvicides, hand-held thermal fogging and vehicle-mounted ultra-low volume fogging. Bacillus thuringiensis israelensis and diflubenzuron are used as larvicides while deltamethrin is the main sprayed adulticide.

Social, environmental and meteorological data

Observations were georeferenced within the French national address database https://adresse.data.gouv.fr/) depending on the available information. When a specific address was not available, the localization of an observation was characterized at the municipality level within the specific département. Social, environmental and meteorological explanatory variables were selected according to results from other studies in the literature [30,45–47].

Socioeconomic and demographic data were locally obtained at the IRIS Census unit level for georeferenced observations. IRIS are aggregated units used in France for statistical purposes. They constitute the smallest geographical unit for which population census data with housing and socioeconomic details are available in France. They are defined with a target size of 2000 inhabitants per basic unit and data for each IRIS are supplied by the French National Institute of Statistics and Economic Studies (INSEE, www.insee.fr). They provide information on the following variables: the characteristics of each household and of the persons who compose that household (marital status, activity/inactivity, socio-professional category), the characteristics of each family (number of children), the use of residences (principal residences, secondary residences, vacant) as well as the proportions of the type of dwelling (houses versus apartments). We used the French social deprivation index (FDep09), a socioeconomic indicator which is available at the IRIS scale [48], to assess the influence of social inequity.

Normalized difference vegetation indexes (NDVI) were derived from Landsat 7 and Landsat 8 satellite imagery data (level 2A) processed by the THEIA Land Data Centre. For each georeferenced observation, an image was selected for NDVI calculation, with the image shooting date as close as possible to the date of symptom onset of the case and with cloud cover less than 10%. Mean NDVI were calculated within buffers of 300m around each geo-located observation.

Daily meteorological measurements (precipitation, minimal and maximal temperature, diurnal temperature range) were obtained from the French meteorological agency (Météo France) for the period from 1 March to 30 November of each study year. For each observation, the closest meteorological station was selected from among the 147 stations located within or in the direct vicinity of the study area. The average distance between case location and the associated meteorological station was 6.4 km (sd = 3.6 km) for temperature data and 6.2 km (sd = 3.6) for rainfall data. Mean (Tmean), maximal (Tmax) and minimal (Tmin) temperatures were computed for the 7 and 10 days before and after EDP, respectively. A bounded accumulated Growing Degree Days (GDD) index was also calculated with a baseline temperature of 11°C and a maximum threshold of 1350°C as proposed by Roiz et al. [36] at EDP (GDD0) and 10 days after EDP (GDD10). Different time windows were considered for calculation of weekly accumulated rainfall in the 1 to 4 weeks before EDP. Weekly diurnal temperature ranges (DTR) were also computed after EDP. A period of 7 to 10 days was chosen for the construction of different temporal variables, as this period corresponds to the duration proposed in the literature for the extrinsic incubation of DENV and CHIKV in Ae. albopictus [49].

Land cover data were extracted from a map produced from optical imaging by the THEIA Land Data Centre [50]. This map is available at the national scale with 17 land cover classes and 10 m spatial resolution. Housing and vegetation were also characterized with the BD TOPO produced by the French National Geographic Institute [51]. Land cover information was extracted in buffers of different sizes, i.e. 100 m, 200 m and 300 m radius zones. This range of buffers was chosen according to existing knowledge on Ae. albopictus dispersal [52].

Statistical analysis

To identify factors associated with the emergence of autochthonous transmission events of DENV and CHIKV in mainland France, we based our analysis on all imported cases identified by the national surveillance system within the study area during the study period. For every imported case, we considered each place where the case was present for at least an hour during viremia as an observation. This implies that a single notified case could result in more than one observation, depending upon the case’s movements and the number of places where the case stayed during the infectious period. Only observations in localities known to be colonized by Ae. albopictus during stays of notified cases were included.

Finally, we compared observations that led to an autochthonous transmission event (n = 13) with observations which did not (n = 844) in order to assess the influence of the above possible explanatory variables on the occurrence of dengue and chikungunya autochthonous events. Comparison was made using complementary log-log regression. The binary response variable of interest was therefore the occurrence of an autochthonous event, coded as “1” in the case of autochthonous transmission and coded as “0” in the absence of autochthonous transmission. We used a generalized linear model with a ‘complementary log-log’ link function, as the probability of the event occurring appeared small.

Univariate binomial regression models were first used to identify candidate variables for the multivariate analysis, while avoiding collinearity. Variables with a p-value <0.25 (arbitrarily chosen) in the univariate analyses were selected as candidates for the multivariate analysis. Spearman’s correlation coefficient was calculated for pairwise variables. Based on the correlation analysis, different sets of variables were defined to ensure that collinearity was reduced. Multivariate models were then built based on univariate analyses using the various sets of variables that did not have statistically significant pairwise correlations. Variance Inflation Factors (VIF) were used to assess multicollinearity between selected variables [53]. Prior to performing a global multivariate analysis, different multivariate models were built for each type of factor (surveillance, climatic, socioeconomic and environmental factors). These analyses (hereinafter referred to as “multivariate sectoral analyses”) were performed to assess the variance explained by each category of variable and to take into account the difference in the number of observations for each type of factor. Variables were selected for multivariate analyses using forward and backward selection. Best-fit models were selected on the basis of Akaike Information Criteria. Statistical analyses was performed using R software [54] with MASS [55] and MuMIn packages [56] for model selection. A general diagram of the analysis strategy is shown in Fig 2.

An external file that holds a picture, illustration, etc. Object name is pntd.0008320.g002.jpg
General outline of the analysis strategy.

Management of missing data

Social, surveillance, environmental and meteorological data were missing for some observations as no imported cases were identified for four events of autochthonous dengue transmission (Table 2). Two of these four events led to more than one autochthonous case. For these clusters, a likely place of transmission was identified and the associated focal point of transmission consequently set at the centre of the transmission area. This hypothesis seems reasonable since all clusters occurred in an area with a radius of less than 300 meters. For events with only one autochthonous case, two scenarios were formulated, the first (Sc1) at the residence of the autochthonous case, while the second (Sc2) considered the absence of georeferenced information. For meteorological data, the time of introduction was set at 30 days prior to the date of symptoms onset of the first autochthonous case. This period was chosen as it is the mean time interval between the occurrence of the first autochthonous case for the five dengue transmission events for which an imported index case was identified and the date of return of imported primary dengue cases (n = 5, Table 2). We only relied on dengue cases as this situation (i.e., a single autochthonous case in the absence of any identified imported case) only occurred for dengue.

We then considered different options (detailed below) to assess the model's response to variations in the “reporting delay” variable. A threshold in the reporting delay was defined arbitrarily at 21 days. This duration corresponds to the average cycle of virus transmission from the mosquito infective blood meal to the end of viremia in the first autochthonous human case. We considered the following options:

  • Option 1: only values recorded prior to the identification of the autochthonous focus retained;

  • Option 2: values recorded after identification of the outbreak introduced;

  • Option 3: all missing values and delays exceeding 21 days set at 21 days;

  • Option 4: this scenario is similar to that of scenario 3, but the variable is converted to a categorical variable. Delay during [0;21 [set as “short to medium” and delay ≥ 21 days as “long”;

  • Option 5: all missing data (including reports of imported cases after autochthonous circulation) classified in a dedicated group, “missing” group.

These different options were considered using univariate analysis to select the one that would be retained in subsequent multivariate analyses.

Ethics statements

This study was part of French national public health surveillance program for vector-borne diseases at Santé publique France (the French Agency for Public Health), a governmental agency reporting to the French Ministry of Health. All data were anonymized. Data collection through the epidemiological surveillance system was approved by the French Commission on Information Technology and Liberties (CNIL), with authorization n°911185.

Results

A total of 857 observations were included in the different analyses. This number corresponds to imported cases who met all four of the following criteria: (i) stayed for more than an hour in areas where the vector was present during the time they were there, (ii) were still viremic, (iii) were classified either as a confirmed or probable case, and (iv) had at least one spatial indication of stay at the municipality level.

Univariate analyses

Twenty-seven variables with a p value < 0.25 (Table 4) in the univariate analyses were selected for multivariate analysis (S1 Table).

Table 4

Results of the univariate binomial regression analyses sorted according to their p-values.
VariableType(1)D2p-valueOR (95% CI)
Tmean 7 days before EDP (Tmean7)M8.8%0.00131.395 (1.148–1.723)
Tmin 7 days before EDP (Tmin7)M7.4%0.00221.346 (1.115–1.642)
Reporting delay (RD)SC7.6%0.00241.042 (1.008–1.065)
Intervention delaySC9.6%0.00261.068 (1.014–1.111)
Tmax 7 days before EDP (Tmax7)M6.9%0.00321.264 (1.086–1.483)
Tmax 10 days after EDP (Tmax10)M6.5%0.00471.273 (1.082–1.507)
Tmean 10 days after EDP (Tmean10)M6.2%0.00731.345 (1.094–1.683)
Discontinuous urban fabrics within a radius of 100 m (DUF100)LC6.9%0.00771.011 (1.004–1.020)
Percentage of “house”- type residences at the IRIS scaleSE5.7%0.010311.42 (1.96–88.26)
Percentage of vegetation within a radius of 300 m (Vegetation)LC3.4%0.01751.033 (1.003–1.058)
Bounded growing degree days 10 days after EDP (GDD10)M6.4%0.01811.003 (1.001–1.005)
Sampling delaySC4.7%0.01851.037 (0.993–1.061)
Bounded growing degree days until EDP (GDD0)M5.8%0.02031.002 (1.001–1.005)
Discontinuous urban fabrics within a radius of 300 m (DUF300)LC3.8%0.03321.001 (1.000–1.002)
Tmin 10 days after EDP (Tmin10)M3.3%0.04151.214 (1.012–1.470)
Percentage of families in householdsSE3.2%0.058268.25 (1.21–8364.41)
NDVI within a radius of 300 mM/LC2.5%0.007493.15 (0.70–14501.45)
Diurnal temperature range (DTR)M2.1%0.08941.185 (0.973–1.444)
Day of the yearM2.1%0.09691.010 (0.998–1.023)
Weekly rainfall 3 weeks before EDP (Rain3w)M1.1%0.11951.009 (0.992–1.017)
Percentage of main residences (Main Res.)SE2.7%0.1212210.5 (0.8–597711.4)
Weekly rainfall 2 weeks before EDP (Rain2w)M2.9%0.15150.938 (0.833–1.000)
Percentage of vacant residencesSE1.8%0.15950.000 (0.000–14.634)
Number of buildings within a radius of 300 m (Buildings)LC2.0%0.16100.999 (0.997–1.000)
Length of viremia (days) in the study areaSC2.4%0.18691.301 (0.942–2.131)
Continuous urban fabric within a radius of 300 mLC2.6%0.20000.998 (0.993–1.000)
Weekly rainfall 1 week before EDP (Rain1w)M1.9%0.22870.949 (0.844–1.008)

Only variables below the selected threshold are presented. D2: explained deviance.

DUF: discontinuous urban fabrics (between 30 to 80% of the total surface is impermeable, covered by buildings, roads and artificially surfaced areas); EDP: earliest date of presence of an imported case during viremia in the study area; GDD: bounded growing degree days; NDVI: normalized difference vegetation index.

(1) The different variables are classified as follows. LC: variable related to land cover. S: variable related to surveillance and control. SE: socioeconomic variable. M: meteorological variable. RD: reporting delay (in days), as the period between the earliest date of presence of an imported case and the date of case reporting.

Multivariate analyses

Collinearity was mainly present due to the construction of the variables to be studied (e.g., temperatures and GDD with different time steps or land cover characteristics with different buffers). The different variables linked to surveillance and response activities (i.e., “Reporting delay”, “Sampling delay” and “intervention delay”) were strongly correlated. Considering variables linked to surveillance and response, using reporting delay instead of sampling delay or intervention delay resulted in a lower Aikake Information Criterion (AIC) and more explained deviance. In subsequent analyses we therefore chose to focus on the reporting delay (RD). Moreover, the RD reflects the time required for a case to be registered in the public health system and constitutes an operational reality.

Option 5 was selected for the transformation of the RD variable. The RD was therefore converted into a categorical variable (RDC) for the subsequent analyses as (1) it exhibited the lowest AIC score among comparable models, (2) had high explained deviance, and (4) allowed us to keep the entire dataset without imputation of missing data for transmission events. All the results of the binomial regression univariate model for different options of the RD on the risk of autochthonous arbovirus transmission are shown in S1 Appendix.

Different models were therefore built by combination of uncorrelated variables. For all the different models, the VIF for each predictor was less than 3, which is the value proposed as a threshold below which any effect of multicollinearity can be considered negligible [57]. Global multivariate models were built as a final step by combining the different type of variables with the reporting delay converted as a categorical variable. Identical global multivariate models were obtained after stepwise selection for scenarios Sc1 and Sc2. However, global multivariate models for scenario Sc1 (scenario for which the transmission location is fixed at the residence of the autochthonous case) exhibits the highest explained deviance. The different models obtained for scenario Sc1 are presented in Table 5. Models obtained for scenario Sc2 are shown in S2 Table.

Table 5

Results of binomial regression of emergence of autochthonous arbovirus infections for scenario Sc1 (place of transmission at the residence of the autochthonous case).
ModelVariablesVariable categorydflogLikAICcdeltaweightD2
Global multivariate models     
Sc1.1RD, Vegetation, GDD0All5-24.94559.970.000.6862.2%
Sc1.2RD, Vegetation, GDD10All5-25.74561.571.600.3161.0%
Multivariate sectoral models   
Sc1.3RDS3-31.20468.448.470.01052.8%
Sc1.4Vegetation, Main Res., DUF100LC/SE4-57.17122.3962.432x10-1413.4%
Sc1.5GDD0, Rain3w, Tmean10M4-57.291122.6362.672x10-1413.3%
Sc1.6Vegetation, Main Res., Houses, DUF300LC/SE5-56.543123.1663.1910−1414.4%
Sc1.7GDD10, Rain3w, Tmean10M4-57.560123.1763.2010−1412.9%
Sc1.8GDD0, Rain3w, DTR, Tmin10M5-56.903123.8863.929x10-1513.9%
Sc1.9GDD0, Rain3w, Tmax10M4-57.937123.9363.969x10-1512.3%
Sc1.10GDD0, Rain3w, Tmax7M4-57.937123.9363.969x10-1512.3%
Sc1.11Vegetation, Main Res., DUF300LC/SE4-57.952123.9563.999x10-1512.3%
Sc1.12GDD10, Rain3w, Tmax10M4-58.101124.2564.288x10-1512.0%
Sc1.13GDD10, Rain3w, Tmax7M4-58.101124.2564.288x10-1512.0%
Sc1.14GDD10, Rain3w, DTR, Tmin10M5-57.176124.4364.467x10-1513.4%
Sc1.15GDD0, Rain3w, Tmean7M4-58.305124.6664.696x10-1511.7%
Sc1.16GDD10, Rain3w, Tmean7M4-58.424124.9064.935x10-1511.5%
Sc1.17GDD0, Rain3w, DTR Tmin7M5-57.905125.8965.923x10-1512.3%
Sc1.18GDD10, Rain3w, DTR Tmin7M5-58.076126.2366.263x10-1512.1%

AICc: Aikake Information Criterion with a correction for small sample sizes; All: all categories of variables are included in the model (surveillance, meteorological, land cover and socioeconomic data); D2: explained deviance

LC/SE: only land cover and socioeconomic data are included as explanatory variables

M: only the meteorological data are included as explanatory variables; S: only the surveillance data are included as explanatory variables; Other variables used are provided in Table 4 and are detailed in S1.

We retained the global model Sc1.1 (Table 6) as the final model as it had the lowest AIC value among all models. The details of the other models, for scenarios Sc1 and Sc2, are provided in S3 Table.

Table 6

Results of binomial regression of autochthonous arboviral case emergence for model Sc1.1.
VariableCoeff.SEz-valueOR (95% CI)p-value
Intercept-10.4002.443-4.2573x10-5 (5x10-8–0.001)<0.001
RD “long” #2.9640.8203.61519.4 (3.4–112.4)<0.001
RD “missing” #9.570121.7840.079-NS
GDD00.0040.0022.0211.004 (1.001–1.009)<0.05
Vegetation0.0470.0172.7651.05 (1.008–1.08)<0.01

(#) “Short” RD (<21 days) as reference level. NS: non-significant

Explained deviance of the model: 62.2%

The RD of imported cases was the factor with the greatest explained deviance: between 48.9 and 52.8% according to the chosen scenario (Table 5). The proportion of variance explained was of the same order of magnitude for meteorological and land cover variables with values between 10 and 15% according to the specific model and scenario. Results of Sc1 and Sc2 were very similar. However, the sectoral analysis showed that explanatory variables related to land use were more sensitive to a lack of observations. For meteorological variables, no major difference was observed between the two constructions of GDD (i.e., at EDP and 10 days after EDP).

Discussion

DENV and CHIKV are emerging threats in Europe. A better understanding of the respective contribution of the main determinants and drivers of emergence of these viruses is needed to identify most-at-risk situations, to prioritize interventions and, ultimately, to adopt a proactive surveillance scheme and implement an adequate public health response. To the best of our knowledge, this is the first comprehensive work exploring different determinants of DENV and CHIKV emergence in temperate settings using an epidemiological dataset.

Reporting-based failures constitute a major factor in the occurrence of autochthonous transmission. A delay in case identification has already been indicated as a contributing factor in contexts of extended viral circulation [58], and remains the most important factor for the occurrence of foci of limited transmission. Our current findings are consistent with our previous work based on a mechanistic approach [59]. Reporting delay is the combination of two elements: (1) how promptly the case-patient seeks medical care and (2) how responsive surveillance partners (e.g., medical analysis laboratories, hospitals, general practitioners, etc.) are in reporting cases. Any action which can positively impact one or the other of these two elements can contribute to the improvement of the arbovirosis surveillance system. Examples of such actions include the continuing effort to raise awareness of this surveillance system in healthcare professionals, the ongoing endeavour to consolidate a robust network of reporting diagnostic laboratories [8,19,29] and raising patient awareness of the importance of seeking medical consultation for non-specific febrile syndromes. The latter action however is more challenging as it is covers both the promotion of travel health consultation—which is not sought by the majority of tourists visiting at risk areas [60,61]—and the information provided to travellers in key locations such as airports. One avenue for research in this area would be to study whether socioeconomic determinants can explain delay or absence in seeking consultation [62], with a view to implementing targeted awareness-raising actions. Moreover, the surveillance system does not allow for the exhaustive detection of all imported viremic infections because of the existence of subclinical infections which are able to infect mosquitoes [63] and which may constitute the source of autochthonous transmission. This may partly explain why, for several autochthonous transmission events studied in France to date, a primary case could not be identified despite reinforced surveillance which carried out door-to-door surveys to identify cases [42]. Preparing for such a contingency (i.e., lack of primary case identification because of subclinical infections) requires maintaining a high level of vigilance so that any autochthonous transmission event can be promptly detected. Moreover, reducing the risk of transmission implies keeping vector populations as low as possible through encouraging the general public to reduce mosquito breeding sites, as larval control remains the main sustainable measure of prevention and control [27]. To a lesser extent, meteorological and environmental factors are associated with autochthonous transmission. The influence of temperature and consequently, the value of using Growing Degree Days (GDD) to describe vector population dynamics has already been reported in different European settings [36,64–66]. However, to the best of our knowledge, the present work constitutes the first time that this indicator has been used together with epidemiological data. GDD can be considered a proxy of the vector population density and reflects the importance of this parameter for estimating transmission risk. The presence of wooded areas is a landscape factor that could explain a favourable environment for Ae. albopictus. It is important to highlight that such areas were associated with anthropization in the present work, since they were located in the vicinity of the places where cases stayed. These green areas, or vegetated areas, are located either inside or are in direct proximity to peri-urban and residential areas which are known to provide suitable conditions for Ae. albopictus proliferation. Numerous rainwater collection containers dedicated to gardening are present and represent the most productive breeding sites for Ae. albopictus in temperate settings. Moreover, vegetation cover maintains relative humidity and provides shelter for adult mosquitoes to rest. Plant residues constitute a major food resource for mosquito larvae [67] and plant sugars are needed as an energy source for both male and female mosquitoes [68]. Finally, the presence of gardens, terraces and other green areas promotes outdoor living and therefore greater human exposure to the exophilic species that is Ae. albopictus [69].

Our study has limitations. We focused on the identification of drivers of dengue and chikungunya emergence. All events of autochthonous transmission were therefore considered, irrespective of their size. However, one may suppose that the factors explaining the size of the foci of transmission may be different from factors which only explain the occurrence of transmission. Furthermore, DENV and CHIKV were considered to be similar, despite differences in transmission efficacy. However, these vector-virus interactions are already highly heterogeneous for each of these pathogens due to the existence of different viral genotypes [70,71] and these two arboviruses exhibit a similar ecology calling for a similar public health response. These different considerations justify, in our view, treating both similarly. Furthermore, if we had had more events, we could have had studied the importance of virus genetic factors in more detail. Host and viral genetic factors may have a major impact on arbovirus transmission, as illustrated by the increased efficacy of transmission observed after the adaptation of chikungunya virus to Ae. albopictus [72]. The results of this work could be extended to the risk of transmission of other arboviral diseases, especially Zika, as the first case of autochthonous vector-borne transmission was recently reported in France [73]. Only events detected by the national surveillance system are reported here, and we cannot exclude that autochthonous transmission events—particularly of modest size—may have gone unnoticed, due to the absence of symptoms, medical consultation, proper diagnosis and reporting. The subclinical infection rates of dengue and chikungunya infections is also a possible limitation of the study, as the real incidence rates of imported cases were probably higher than those observed. Accordingly, the different risk factors identified were potentially overestimated. However, we assume that the impact of under-detection and under-diagnosis was limited given the efficiency of the surveillance system dedicated to dengue and chikungunya infections. Aedes albopictus is still expanding its range throughout the world, including Europe. However, it is difficult to extrapolate our results to other bioclimatic zones. It is therefore essential to update this type of approach in the light of new transmission events and in different environmental and climatic settings. Moreover, the presence of nested data (as a unique case introduction could lead to several observations) would have justified the use of mixed methods, but the limited number of transmission events did not contain any sufficient information for parameters estimation in such a framework. This limited number of transmission events is explained by the epidemiological situation and requires to consider the results with caution. In this global perspective, both routine documentation of viremic cases (from a meteorological, environmental, socio-economic point of view) and routine integration of data into surveillance system databases is paramount. The development of these anticipatory tools will be useful to reduce the risk of multiple autochthonous arboviral transmission events.

Supporting information

S1 Table

Description of the different explanatory variables.

(DOCX)

S2 Table

Results of binomial regression of autochthonous arboviral case emergence for scenario Sc2.

(DOCX)

S3 Table

Results of binomial regression of autochthonous arboviral case emergence for the different global multivariate models.

(DOCX)

S1 Appendix

Management of missing data of the reporting delay.

(DOCX)

Acknowledgments

Our thanks to Grégoire Rey, from the French Institute of Health and Medical Research (INSERM,) for the provision of the most recent French Deprivation index (FDep09) data. We are also grateful to Joël Deniau and Sandra Giron (Santé publique France, Provence-Alpes-Côte-d’Azur-Corse), Damien Mouly and Amandine Cochet (Santé publique France, Occitanie) for epidemiological data collection and management. We take this opportunity to thank all the health professionals involved in France’s surveillance and response system.

Funding Statement

FJ is funded by Santé publique France (www.santepubliquefrance.fr) through a doctoral fellowship. Meteorological data were provided by Météo-France, the national meteorological agency, within the framework of the agreement linking the agency with Santé Publique France. Access to satellite images used in this work was supported by public funds received under the framework of GEOSUD, a project (ANR-10-EQPX-20) with the "Investissements d'Avenir" program, managed by the French National Research Agency. Funders do not play any role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data cannot be shared publicly because of restrictions by the French National Commission on Informatics and Liberty (CNIL) concerning the access to personnal and health data. Data can be requested by interested parties at Bruno Coignard, [email protected], Director, Infectious Diseases at Santé publique France.

References

1. Mayer S V, Tesh RB, Vasilakis N. The emergence of arthropod-borne viral diseases: A global prospective on dengue, chikungunya and Zika fevers. Acta Trop. 2017;166: 155–163. 10.1016/j.actatropica.2016.11.020 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
2. Sigfrid L, Reusken C, Eckerle I, Nussenblatt V, Lipworth S, Messina J, et al. Preparing clinicians for (re-)emerging arbovirus infectious diseases in Europe. Clin Microbiol Infect. 2018;24: 229–239. 10.1016/j.cmi.2017.05.029 [Abstract] [CrossRef] [Google Scholar]
3. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496: 504–507. 10.1038/nature12060 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
4. Zeller H, Van Bortel W, Sudre B. Chikungunya: Its History in Africa and Asia and Its Spread to New Regions in 2013–2014. J Infect Dis. 2016;214: S436–S440. 10.1093/infdis/jiw391 [Abstract] [CrossRef] [Google Scholar]
5. Wilder-Smith A, Gubler DJ. Geographic expansion of dengue: the impact of international travel. Med Clin North Am. 2008;92: 1377–90, x. 10.1016/j.mcna.2008.07.002 [Abstract] [CrossRef] [Google Scholar]
6. Van Bortel W, Dorleans F, Rosine J, Blateau A, Rousseau D, Matheus S, et al. Chikungunya outbreak in the Caribbean region, December 2013 to March 2014, and the significance for Europe. Eurosurveillance. European Centre for Disease Prevention and Control; 2014;19: 20759 10.2807/1560-7917.ES2014.19.13.20759 [Abstract] [CrossRef] [Google Scholar]
7. Medlock JM, Hansford KM, Schaffner F, Versteirt V, Hendrickx G, Zeller H, et al. A Review of the Invasive Mosquitoes in Europe: Ecology, Public Health Risks, and Control Options. Vector-Borne Zoonotic Dis. 2012;12: 435–447. 10.1089/vbz.2011.0814 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
8. Franke F, Giron S, Cochet A, Jeannin C, Leparc-Goffart I, de Valk H, et al. [Autochthonous chikungunya and dengue fever outbreak in Mainland France, 2010–2018]. Bull Epidémiol Hebd. 2019; 374–82. [Google Scholar]
9. Angelini R, Finarelli AC, Angelini P, Po C, Petropulacos K, Macini P, et al. An outbreak of chikungunya fever in the province of Ravenna, Italy. Wkly releases. 2007;12: E070906.1 10.2807/esw.12.36.03260-en [Abstract] [CrossRef] [Google Scholar]
10. European Centre for Disease Prevention and Control (ECDC). Local transmission of dengue fever in France and Spain– 2018–22 October 2018 [Internet]. 2018 [cited 20 Nov 2019]. Available: https://www.ecdc.europa.eu/sites/portal/files/documents/08-10-2018-RRA-Dengue-France.pdf.
11. Venturi G, Di Luca M, Fortuna C, Remoli ME, Riccardo F, Severini F, et al. Detection of a chikungunya outbreak in Central Italy, August to September 2017. Eurosurveillance. 2017;22 10.2807/1560-7917.ES.2017.22.39.17-00646 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
12. Gjenero-Margan I, Aleraj B, Krajcar D, Lesnikar V, Klobučar A, Pem-Novosel I, et al. Autochthonous dengue fever in Croatia, August-September 2010. Euro Surveill. 2011;16. [Abstract] [Google Scholar]
13. Calba C, Guerbois-Galla M, Franke F, Jeannin C, Auzet-Caillaud M, Grard G, et al. Preliminary report of an autochthonous chikungunya outbreak in France, July to September 2017. Eurosurveillance. 2017;22 10.2807/1560-7917.ES.2017.22.39.17-00647 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
14. Succo T, Leparc-Goffart I, Ferré J, Roiz D, Broche B, Maquart M, et al. Autochthonous dengue outbreak in Nîmes, South of France, July to September 2015. Eurosurveillance. 2016;21: 1–7. 10.2807/1560-7917.ES.2016.21.21.30240 [Abstract] [CrossRef] [Google Scholar]
15. Delisle E, Rousseau C, Broche B, Leparc-Goffart I, L’Ambert G, Cochet A, et al. Chikungunya outbreak in Montpellier, France, September to October 2014. Euro Surveill. 2015;20. [Abstract] [Google Scholar]
16. Marchand E, Prat C, Jeannin C, Lafont E, Bergmann T, Flusin O, et al. Autochthonous case of dengue in France, October 2013. Euro Surveill. 2013;18: 20661 10.2807/1560-7917.es2013.18.50.20661 [Abstract] [CrossRef] [Google Scholar]
17. Grandadam M, Caro V, Plumet S, Thiberge J-M, Souarès Y, Failloux A-B, et al. Chikungunya Virus, Southeastern France. Emerg Infect Dis. 2011;17: 910–913. 10.3201/eid1705.101873 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
18. La Ruche G, Souarès Y, Armengaud A, Peloux-Petiot F, Delaunay P, Desprès P, et al. First two autochthonous dengue virus infections in metropolitan France, September 2010. Euro Surveill. 2010;15: 19676 [Abstract] [Google Scholar]
19. Paty MC, Six C, Charlet F, Heuzé G, Cochet A, Wiegandt A, et al. Large number of imported chikungunya cases in mainland France, 2014: a challenge for surveillance and response. Euro Surveill. 2014;19: 20856 10.2807/1560-7917.es2014.19.28.20856 [Abstract] [CrossRef] [Google Scholar]
20. European Centre for Disease Prevention and Control (ECDC), European Food Safety Authority (EFSA). Mosquito maps. Stockholm. [Internet]. 2018 [cited 1 Aug 2019]. Available: https://ecdc.europa.eu/en/disease-vectors/surveillance-and-disease-data/mosquito-maps
21. Zouache K, Fontaine A, Vega-Rua A, Mousson L, Thiberge JM, Lourenco-De-Oliveira R, et al. Three-way interactions between mosquito population, viral strain and temperature underlying chikungunya virus transmission potential. Proc Biol Sci. 2014;281 10.1098/rspb.2014.1078 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
22. Lambrechts L, Chevillon C, Albright RG, Thaisomboonsuk B, Richardson JH, Jarman RG, et al. Genetic specificity and potential for local adaptation between dengue viruses and mosquito vectors. BMC Evol Biol. 2009;9 10.1186/1471-2148-9-160 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
23. Jourdain F, Roiz D, Perrin Y, Grucker K, Simard F, Paupy C. [Entomological factors of arboviruses emergences]. Transfus Clin Biol. 2015/07/05. 2015;22: 101–106. 10.1016/j.tracli.2015.06.001 [Abstract] [CrossRef] [Google Scholar]
24. Semenza JC, Lindgren E, Balkanyi L, Espinosa L, Almqvist MS, Penttinen P, et al. Determinants and Drivers of Infectious Disease Threat Events in Europe. Emerg Infect Dis. 2016;22: 581–9. 10.3201/eid2204 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
25. Semenza JC, Menne B. Climate change and infectious diseases in Europe. Lancet Infect Dis. 2009;9: 365–75. 10.1016/S1473-3099(09)70104-5 [Abstract] [CrossRef] [Google Scholar]
26. Liu T, Zhu G, He J, Song T, Zhang M, Lin H, et al. Early rigorous control interventions can largely reduce dengue outbreak magnitude: experience from Chaozhou, China. BMC Public Health. 2018;18: 90 10.1186/s12889-017-4616-x [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
27. Roiz D, Wilson AL, Scott TW, Fonseca DM, Jourdain F, Müller P, et al. Integrated Aedes management for the control of Aedes-borne diseases. Horstick O, editor. PLoS Negl Trop Dis. 2018;12: e0006845 10.1371/journal.pntd.0006845 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
28. Terrien E, Fournet N, Giron S, Franke F, Cochet A, Calba C, et al. [Chikungunya, dengue and Zika virus surveillance in metropolitan France, 2018]. Bull Epidémiol Hebd. 2019; 363–73. [Google Scholar]
29. Septfons A, Leparc-Goffart I, Couturier E, Franke F, Deniau J, Balestier A, et al. Travel-associated and autochthonous Zika virus infection in mainland France, 1 January to 15 July 2016. Euro Surveill. European Centre for Disease Prevention and Control; 2016;21 10.2807/1560-7917.ES.2016.21.32.30315 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
30. Sallam MF, Fizer C, Pilant AN, Whung P-Y. Systematic Review: Land Cover, Meteorological, and Socioeconomic Determinants of Aedes Mosquito Habitat for Risk Mapping. Int J Environ Res Public Health. Multidisciplinary Digital Publishing Institute (MDPI); 2017;14: 1230 10.3390/ijerph14101230 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
31. Fan J, Wei W, Bai Z, Fan C, Li S, Liu Q, et al. A Systematic Review and Meta-Analysis of Dengue Risk with Temperature Change. Int J Environ Res Public Health. 2014;12: 1–15. 10.3390/ijerph120100001 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
32. Banu S, Hu W, Hurst C, Tong S. Dengue transmission in the Asia-Pacific region: impact of climate change and socio-environmental factors. Trop Med Int Heal. 2011;16: 598–607. 10.1111/j.1365-3156.2011.02734.x [Abstract] [CrossRef] [Google Scholar]
33. Ren H, Zheng L, Li Q, Yuan W, Lu L. Exploring Determinants of Spatial Variations in the Dengue Fever Epidemic Using Geographically Weighted Regression Model: A Case Study in the Joint Guangzhou-Foshan Area, China, 2014. Int J Environ Res Public Health. 2017;14: 1518 10.3390/ijerph14121518 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
34. Waldock J, Chandra NL, Lelieveld J, Proestos Y, Michael E, Christophides G, et al. The role of environmental variables on Aedes albopictus biology and chikungunya epidemiology. Pathog Glob Health. 2013;107: 224–241. 10.1179/2047773213Y.0000000100 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
35. Little E, Bajwa W, Shaman J. Local environmental and meteorological conditions influencing the invasive mosquito Ae. albopictus and arbovirus transmission risk in New York City. Lenhart A, editor. PLoS Negl Trop Dis. 2017;11: e0005828 10.1371/journal.pntd.0005828 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
36. Roiz D, Boussès P, Simard F, Paupy C, Fontenille D. Autochthonous Chikungunya transmission and extreme climate events in Southern France. Carvalho MS, editor. PLoS Negl Trop Dis. 2015;9: e0003854 10.1371/journal.pntd.0003854 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
37. Roche B, Léger L, L’Ambert G, Lacour G, Foussadier R, Besnard G, et al. The spread of Aedes albopictus in Metropolitan France: Contribution of environmental drivers and human activities and predictions for a near future. PLoS One. 2015;10 10.1371/journal.pone.0125600 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
38. European Commission. NUTS—Nomenclature of territorial units for statistics [Internet]. [cited 6 Jun 2019]. Available: https://ec.europa.eu/eurostat/web/nuts/background
39. French National Institute of Statistics and Economic Studies (INSEE). Estimation de la population au 1er janvier 2019 [Internet]. [cited 6 Jun 2019]. Available: https://www.insee.fr/fr/statistiques
40. Joly D, Brossard T, Cardot H, Cavailhes J, Hilal M, Wavresky P. Les types de climats en France, une construction spatiale. Cybergeo Eur J Geogr. 2010; 10.4000/cybergeo.23155 [CrossRef] [Google Scholar]
41. Lacour G, Chanaud L, L’Ambert G, Hance T. Seasonal Synchronization of Diapause Phases in Aedes albopictus (Diptera: Culicidae). Benoit JB, editor. PLoS One. 2015;10: e0145311 10.1371/journal.pone.0145311 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
42. Giron S, Rizzi S, Leparc-Goffart I, Septfons A, Tine R, Cadiou B, et al. Nouvelles apparitions de cas autochtones de dengue en région Provence-Alpes-Côte d’Azur, France, août-septembre 2014. Bull épidémiologique Hebd. 2015;13–14: 217–223. [Google Scholar]
43. Thiberville S-D, Moyen N, Dupuis-Maguiraga L, Nougairede A, Gould EA, Roques P, et al. Chikungunya fever: Epidemiology, clinical syndrome, pathogenesis and therapy. Antiviral Res. 2013;99: 345–370. 10.1016/j.antiviral.2013.06.009 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
44. Chan M, Johansson MA. The incubation periods of Dengue viruses. PLoS One. 2012;7: e50972—e50972. 10.1371/journal.pone.0050972 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
45. Aswi A, Cramb SM, Moraga P, Mengersen K. Bayesian spatial and spatio-temporal approaches to modelling dengue fever: a systematic review. Epidemiol Infect. 2019;147: e33 10.1017/S0950268818002807 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
46. Carbajo AE, Cardo M V., Guimarey PC, Lizuain AA, Buyayisqui MP, Varela T, et al. Is autumn the key for dengue epidemics in non endemic regions? The case of Argentina. PeerJ. 2018;6: e5196 10.7717/peerj.5196 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
47. Qu Y, Shi X, Wang Y, Li R, Lu L, Liu Q. Effects of socio-economic and environmental factors on the spatial heterogeneity of dengue fever investigated at a fine scale. Geospat Health. 2018;13 10.4081/gh.2018.682 [Abstract] [CrossRef] [Google Scholar]
48. Rey G, Jougla E, Fouillet A, Hémon D. Ecological association between a deprivation index and mortality in France over the period 1997–2001: variations with spatial scale, degree of urbanicity, age, gender and cause of death. BMC Public Health. 2009;9: 33 10.1186/1471-2458-9-33 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
49. Vega-Rua A, Zouache K, Caro V, Diancourt L, Delaunay P, Grandadam M, et al. High Efficiency of Temperate Aedes albopictus to Transmit Chikungunya and Dengue Viruses in the Southeast of France. Moreira LA, editor. PLoS One. 2013;8: e59716 10.1371/journal.pone.0059716 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
50. Inglada J, Vincent A, Arias M, Tardy B, Morin D, Rodes I. Operational High Resolution Land Cover Map Production at the Country Scale Using Satellite Image Time Series. Remote Sens. 2017;9: 95 10.3390/rs9010095 [CrossRef] [Google Scholar]
51. Maugeais E, Lecordix F, Halbecq X, Braun A. Dérivation cartographique multi échelles de la BDTopo de l’IGN France: mise en oeuvre du processus de production de la Nouvelle Carte de Base. Proc 25th Int Cartogr Conf Paris. 2011; 3–8. [Google Scholar]
52. Marini F, Caputo B, Pombi M, Travaglio M, Montarsi F, Drago A, et al. Estimating Spatio-Temporal Dynamics of Aedes Albopictus Dispersal to Guide Control Interventions in Case of Exotic Arboviruses in Temperate Regions. Sci Rep. Nature Publishing Group; 2019;9: 10281 10.1038/s41598-019-46466-4 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
53. Zuur AF, Ieno EN, Elphick CS. A protocol for data exploration to avoid common statistical problems. Methods Ecol Evol. 2010;1: 3–14. 10.1111/j.2041-210X.2009.00001.x [CrossRef] [Google Scholar]
54. R Core Development Team. R: A language and environment for statistical computing R Foundation for Statistical Computing, Vienna, Austria: Version 3.5.3. [Internet]. 2019 [cited 13 Jun 2019]. Available: https://www.r-project.org/ [Google Scholar]
55. Ripley B, Venables B, Bates DM, Hornik K, Gebhardt A, Firth D. Package MASS. Support Functions and Datasets for Venables and Ripley’s MASS. R Packag version 73–511. 2018; [Google Scholar]
56. Barton K. Package MuMIn. Model selection and model averaging based on information criteria. In: R package version 1.42.1. 2018.
57. Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM. Mixed effects models and extensions in ecology with R. Springer; New York; 2009. [Google Scholar]
58. Liu K-K, Wang T, Huang X-D, WangANG G-L, XiaIA Y, Zhang Y-T, et al. Risk assessment of dengue fever in Zhongshan, China: a time-series regression tree analysis. Epidemiol Infect. 2017;145: 451–461. 10.1017/S095026881600265X [Abstract] [CrossRef] [Google Scholar]
59. Sochacki T, Jourdain F, Perrin Y, Noel H, Paty M, de Valk H, et al. Imported chikungunya cases in an area newly colonized by Aedes albopictus: Mathematical assessment of the best public health strategy. Eurosurveillance. 2016;21 10.2807/1560-7917.ES.2016.21.18.30221 [Abstract] [CrossRef] [Google Scholar]
60. Adou AA, Napolitano F, Vastola A, Angelillo IF. Travelers’ knowledge, attitudes, and behavior related to infectious diseases in Italy. Nisha C, editor. PLoS One. 2019;14: e0215252 10.1371/journal.pone.0215252 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
61. Van Herck K, Castelli F, Zuckerman J, Nothdurft H, Van Damme P, Dahlgren A-L, et al. Knowledge, Attitudes and Practices in Travel-related Infectious Diseases: The European Airport Survey. J Travel Med. 2004;11: 3–8. 10.2310/7060.2004.13609 [Abstract] [CrossRef] [Google Scholar]
62. Yoo H-S, Park O, Park H-K, Lee E-G, Jeong E-K, Lee J-K, et al. Timeliness of national notifiable diseases surveillance system in Korea: a cross-sectional study. BMC Public Health. 2009;9: 93 10.1186/1471-2458-9-93 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
63. Duong V, Lambrechts L, Paul RE, Ly S, Lay RS, Long KC, et al. Asymptomatic humans transmit dengue virus to mosquitoes. Proc Natl Acad Sci U S A. 2015;112: 14688–93. 10.1073/pnas.1508114112 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
64. Neteler M, Roiz D, Rocchini D, Castellani C, Rizzoli A. Terra and Aqua satellites track tiger mosquito invasion: modelling the potential distribution of Aedes albopictus in north-eastern Italy. Int J Heal Geogr. 2011;10: 49. [Europe PMC free article] [Abstract] [Google Scholar]
65. Neteler M, Metz M, Rocchini D, Rizzoli A, Flacio E, Engeler L, et al. Is Switzerland suitable for the invasion of Aedes albopictus? PLoS One. 2013/12/19. 2013;8: e82090 10.1371/journal.pone.0082090 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
66. Manica M, Filipponi F, D’Alessandro A, Screti A, Neteler M, Rosà R, et al. Spatial and Temporal Hot Spots of Aedes albopictus Abundance inside and outside a South European Metropolitan Area. McCall PJ, editor. PLoS Negl Trop Dis. 2016;10: e0004758 10.1371/journal.pntd.0004758 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
67. Yee DA, Allgood D, Kneitel JM, Kuehn KA. Constitutive differences between natural and artificial container mosquito habitats: vector communities, resources, microorganisms, and habitat parameters. J Med Entomol. 2012;49: 482–91. 10.1603/me11227 [Abstract] [CrossRef] [Google Scholar]
68. Samson DM, Qualls WA, Roque D, Naranjo DP, Alimi T, Arheart KL, et al. Resting and energy reserves of Aedes albopictus collected in common landscaping vegetation in St. Augustine, Florida. J Am Mosq Control Assoc. 2013;29: 231–236. 10.2987/13-6347R.1 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
69. Halasa YA, Shepard DS, Fonseca D, Farajollahi A, Healy S, Gaugler R, et al. Quantifying impact of mosquitoes on quality of life. Am J Trop Med Hyg. 2013;89: 175. [Google Scholar]
70. Vega-Rua A, Zouache K, Girod R, Failloux AB, 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. 2014/03/29. 2014;88: 6294–6306. 10.1128/JVI.00370-14 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
71. Niu C, Huang Y, Wang M, Huang D, Li J, Huang S, et al. Differences in the Transmission of Dengue Fever by Different Serotypes of Dengue Virus. Vector-Borne Zoonotic Dis. 2019; vbz.2019.2477. 10.1089/vbz.2019.2477 [Abstract] [CrossRef] [Google Scholar]
72. Tsetsarkin KA, Chen R, Yun R, Rossi SL, Plante KS, Guerbois M, et al. Multi-peaked adaptive landscape for chikungunya virus evolution predicts continued fitness optimization in Aedes albopictus mosquitoes. Nat Commun. 2014;5: 4084 10.1038/ncomms5084 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
73. Giron S, Franke F, Decoppet A, Cadiou B, Travaglini T, Thirion L, et al. Vector-borne transmission of Zika virus in Europe, southern France, August 2019. Eurosurveillance. 2019;24 10.2807/1560-7917.ES.2019.24.45.1900655 [Europe PMC free article] [Abstract] [CrossRef] [Google Scholar]
Articles from PLOS Neglected Tropical Diseases are provided here courtesy of PLOS

Full text links 

Read article at publisher's site: https://doi.org/10.1371/journal.pntd.0008320

Read article for free, from open access legal sources, via Unpaywall: https://journals.plos.org/plosntds/article/file?id=10.1371/journal.pntd.0008320&type=printableUnpaywall

Citations & impact 

Impact metrics

27
Citations
Jump to Citations

Citations of article over time

Alternative metrics

Altmetric item for https://www.altmetric.com/details/81848257
Altmetric
Discover the attention surrounding your research
https://www.altmetric.com/details/81848257

Article citations

Go to all (27) article citations

Data 

Data behind the article

This data has been text mined from the article, or deposited into data resources.

BioStudies: supplemental material and supporting data

Similar Articles 

To arrive at the top five similar articles we use a word-weighted algorithm to compare words from the Title and Abstract of each citation.