Skip to main content

Environmental predictors of West Nile fever risk in Europe

Abstract

Background

West Nile virus (WNV) is a mosquito-borne pathogen of global public health importance. Transmission of WNV is determined by abiotic and biotic factors. The objective of this study was to examine environmental variables as predictors of WNV risk in Europe and neighboring countries, considering the anomalies of remotely sensed water and vegetation indices and of temperature at the locations of West Nile fever (WNF) outbreaks reported in humans between 2002 and 2013.

Methods

The status of infection by WNV in relationship to environmental and climatic risk factors was analyzed at the district level using logistic regression models. Temperature, remotely sensed Normalized Difference Vegetation Index (NDVI) and Modified Normalized Difference Water Index (MNDWI) anomalies, as well as population, birds’ migratory routes, and presence of wetlands were considered as explanatory variables.

Results

The anomalies of temperature in July, of MNDWI in early June, the presence of wetlands, the location under migratory routes, and the occurrence of a WNF outbreak the previous year were identified as risk factors. The best statistical model according to the Akaike Information Criterion was used to map WNF risk areas in 2012 and 2013. Model validations showed a good level of prediction: area under Receiver Operator Characteristic curve = 0.854 (95% Confidence Interval 0.850-0.856) for internal validation and 0.819 (95% Confidence Interval 0.814-0.823) (2012) and 0.853 (95% Confidence Interval 0.850-0.855) (2013) for external validations, respectively.

Conclusions

WNF incidence is increasing in Europe and WNV is expanding into new areas where it had never been observed before. Our model can be used to direct surveillance activities and public health interventions for the upcoming WNF season.

Background

West Nile virus (WNV), a member of the Flavivirus genus (family Flaviviridae), is responsible for West Nile disease (WND), which causes considerable morbidity and mortality worldwide. WND, including West Nile fever (WNF) and West Nile neuro invasive disease (WNND), is an emerging arbovirosis originating from the Old World which was introduced to the Americas in 1999[1]. Today WNV is the most widespread flavivirus in the world, with a geographic distribution covering Africa, Eurasia, and the Americas. The public health burden with the associated human and economic costs present a strong case for antecedent estimates of WND risk. Wild birds are the principal hosts and critical to support enzootic cycles of WNV and virus dispersal[2]. WNV is transmitted from birds to birds by ornithophilic mosquitoes, predominantly of the Culex genus[3]. WNV is pathogenic for horses and humans who are accidental hosts[4].

The transmission, epidemiology and geographic distribution of WND is a complex function of biotic (i.e. host abundance and diversity, bird migration, mosquito distribution, etc.) and abiotic (i.e. physical features of the environment and climate) factors. Both biological factors (i.e. vectors diversity, survival of pathogens within the vectors) and environmental determinants (i.e. climatic conditions, availability of water bodies) have impacts on the transmission cycles[2, 5].

WN virus circulation has been confirmed at erratic intervals in the last decades in several countries of the Mediterranean basin (e.g.[6–9]) and in Eastern and Central Europe (e.g.[10]). However, in the summer of 2010, the number of WND cases in previously uninfected areas in Europe and its neighboring countries was the highest ever reported[11]. In a recent study, positive temperature anomalies were identified as major risk factors in WND outbreak occurrence in 2010 in Europe[12], but environmental variables derived from Earth Observation data had not been tested yet in the European context over large areas[13].

Thus, in this study we examined meteorological and environmental variables related to the temperature, the state of vegetation and water bodies, and birds’ migratory routes as predictors of WND risk.

Results and discussion

WND outbreaks, 2002–2013

Between 2002 and 2009 WND outbreaks in humans were confined to a restricted number of districts in Europe and neighboring countries (Table 1). Sustained transmission within one district over a number of years was rarely observed. However, in the summer of 2010, Europe and its neighboring countries experienced an unprecedented upsurge in the number of WND cases: in Eastern Europe, particularly in Russia, Greece, Romania, and Hungary, and in two Mediterranean countries, Israel and Turkey (Figure 1). In Western Europe, two countries, Spain and Italy, reported a limited number of cases.The outbreaks in Europe and Eurasia during the summers of 2011 and 2012 followed most of the disease locations of 2010 (except in Spain), and reached new areas such as Ukraine, the Western Balkans, Tunisia, and Algeria (Figure 1). As in 2010, most of the cases were reported between the end of July and the end of September. In 2013, 97 districts were affected notably in Italy and in central Europe (Hungary, Romania and Serbia) with persistent circulation of WNV in Russia. Tunisia was also affected.

Table 1 West Nile Fever outbreaks, Europe and neighboring countries, 2002–2013
Figure 1
figure 1

Map of West Nile fever outbreaks by year, Europe and neighboring countries between 2002 and 2011. Note: The map represents the number of WNF presence notification by year over the study period (2002–2011).

Temperature anomalies

Figure 2a presents the monthly anomalies of temperatures for selected month (July) for the years 2010–2013. The maps show that while spring 2010 was very hot in North Africa, the summer of 2010 was severely warm in Eurasia, with extreme anomalies in July (above 6°C) and August. Spring 2011 was warm above normal in Western Europe and the summer was warm in June and August in the central and southern parts of the continent. Again, Eurasia was very warm in July. In 2012, the temperatures of the hot season were above the perennial averages mainly in Eurasia and Southern Europe.

Figure 2
figure 2

Anomalies of temperatures and MNDWI values, 2010–2013. a) Monthly anomalies from the perennial mean monthly temperature (July). b) MNDWI anomalies from the 2002–2011 average (21st MODIS period).

Water index anomalies

Figure 2b presents the anomalies of Modified Normalized Difference Water Index (MNDWI) for one selected Moderate Resolution Imaging Spectroradiometer (MODIS) composite 8-days period (21st MODIS 8-day period: 9th – 16th June) for the years 2010–2013. Maps highlight a high heterogeneity in the anomalies of the water index derived from MODIS imagery throughout Europe. Higher anomalies were observed in Southern and Eastern Europe in June 2010.

Environmental risk factors of WNV infection in Europe

After the univariate analysis, 12 covariates were kept for multi-variate analysis (see Additional file1: Table S1). According to the Akaike Information Criterion (AIC) values, the probability of infection by WNV at the district level is better explained as a function of the anomalies of temperature in July, of MNDWI in the 21st MODIS-8 day period (i.e. June 9th-16th), the weighted average of the number of infected districts amongst the neighborhood the previous year (λ), the presence of wetlands, the type of passerine migratory routes, the population, all positively and highly significantly (p < 0.01) correlated with the probability of infection (Table 2). Normalized Difference Vegetation Index (NDVI) variables were not retained as risk factor for WNV infection. Inspection of the residuals from the final model showed no significant spatial autocorrelation (Durbin-Watson test = 1.99, p-value = 0.23).Area Under Curve (AUC) of the model was 0.854 (95% Confidence Interval (CI) 0.850-0.856) for internal validation, 0.819 (95% CI 0.814-0.823) (2012 data) and 0.853 (95% CI 0.850-0.855) (2013 data) for external validation, indicating good predictions (Receiver operator characteristic curve (ROC) curves are presented in Figure 3).

Table 2 Multivariate logistic regression model parameter of the risk of WNV infection at district level, EU and neighbouring countries
Figure 3
figure 3

ROC curves of the best model of the probability of WNV infection.

Applying the model to 2012 and 2013 data with the best threshold identified through the analysis of the ROC curve (threshold = 0.014, determined by the observed prevalence method), specificity was 0.639 (95% CI 0.633-0.646) and sensitivity was 0.872 (95% CI 0.841-0.875) for 2012 data and specificity was 0.739 (95% CI 0.729-0.749) and sensitivity was 0.788 (95% CI 0.782-0.805) for 2013 data.

Figure 4 shows the probability of WNV infection per district in 2012 and 2013 as predicted by the model (see Additional file2: Figure S1, for 2002–2011 risk maps). The maps highlight spatial heterogeneity of the risk of WNV occurrence throughout Europe. In 2012 higher risk values were predicted for countries of Central and Eastern Europe, Turkey, Israel, and Tunisia. WND cases were reported in all of the predicted high risk areas, excepted in Ukraine, and Turkey. In 2013, Tunisia, Northern Italy, Northern Greece, Central Europe and South Russia presented the highest predicted values in agreement with main areas of transmission.

Figure 4
figure 4

Map of predicted probability of WNV infection based on environmental predictors, Europe and neighboring countries, 2012 and 2013.

WNV is a vector-borne pathogen of global importance since it is the most widely distributed of the encephalitic flaviviruses[14]. Thus, in this study we tested environmental variables as risk factors for WNV transmission in Europe and neighboring countries. Our results identified the anomaly of summer temperatures (July), the anomaly of the MNDWI in early June (June 9th-16th), the occurrence of a WND outbreak the year before, the human population, the presence of wetlands and the type of passerine flyways as risk factors for WNV transmission in Europe.

Thus, our study based on 2002–2011 WND data confirmed the conclusions of Paz et al. (2013) on the linkage between the 2010 heat and the WND eruptions in Europe and its neighboring countries. Similar associations between high temperatures and WND eruption were found in other locations, especially in studies from North America. In a comprehensive study of 17 states in the United States, positive associations with increasing temperature was detected over each of the four weeks prior to symptom onset[15]. Positive relations with heat were found in additional studies for Georgia[16], California[17], Illinois[18], Connecticut[19], and even northwards in the Canadian Prairies[20].

Reisen et al. (2006) indicated that WNV dispersed into new areas during years with above-normal temperatures and that implication during the following year occurred in summers with above- or normal temperatures. This was the case in Europe, as a comparison between the WND distributions in 2010, 2011 and 2012 (Figure 1, Table 1) and the July temperature anomalies (Figure 2a) suggests that the limited appearance of WND in 2011 was related to less warm summer conditions. When the temperature increased in summer 2012, WND number of cases raised (in 2011, 127 and 212 WND cases were reported in the EU and neighboring countries, respectively; increasing in 2012 to 242 and 693 respectively)[11].

These observations are consistent with our results identifying λ, the weighted average of the number of infected districts in a close neighborhood the year before as risk factor for WNV transmission. They support the hypothesis that WNV has the capacity to persist locally in Europe after a first introduction, probably through survival in overwintering mosquitoes[21, 22] or in wild bird populations[23, 24]. This epidemiological pattern was only recently observed in Europe, the Italian 2008 outbreak marking a change in the epidemiology of the disease[25].

Our study stressed the importance of water bodies in the risk of WNV transmission, our results showing that areas including wetlands with positive anomalies of the MNDWI in June are more at risk. The role of water bodies in WNV transmission is related to the availability of breeding sites for the mosquito vector populations and may be complex: on one hand, large areas of surfacewater may favor vector abundance and WNV transmission as demonstrated in other studies[26], but in drought conditions, the reduced number of water pools may lead to higher interactions between wild birds and mosquitoes[27]. Our findings support the first hypothesis, suggesting that in Europe, above average surface of water in June, favoring mosquito proliferation before the summer months, is a significant risk factor for WNV transmission. Yet, entomological field observations are needed at the European scale to assess i) the distribution and dynamics of WNV vectors in Europe and ii) the relationship between mosquitoes and MNDWI dynamics, assessing whether an increase of MNDWI, which is related to an increase of water bodies, is a reliable proxy for an increase of available breeding sites of WNV vectors. On the other hand, our study highlights that although water surfaces are generally closely linked to rainfall events, water indices such as the MNDWI might be more relevant as proxies of WND risk, as the association between WND outbreaks occurrence with precipitation was not consistent[28]. The inconsistent linkages between WND eruptions and precipitation amounts appear in the literature for other regions. For instance, Frost et al. (2012) noted that the eruption of acute equine encephalitis (virulent strain of WNVKUN) in 2011 in Australia was followed by extensive floods across eastern Australia that created ideal breeding conditions for freshwater Culex annulirostris[29]. In studies from the USA, Soverow et al. (2009) found positive associations with rainfall in the months preceding the disease outbreaks; in contrast Shaman et al. (2002) demonstrated that heavy rainfall might lead to a negative association by flushing the ditches and drainage channels used by Culex larvae[30]. On the other side, drought can facilitate population outbreaks of some species of mosquitoes. It was shown by DeGroote et al. that prior drought contributed to the initial United States WNV outbreak[31].The predictions of the best model applied to 2012 and 2013 data performed well according to the AUC, sensitivity and specificity indices (Figure 3). Thus, an important part of the variability in WNV transmission in Europe at the district level can be captured by environmental drivers. Yet, other factors related to the presence and abundance of hosts (wild birds) and vectors (mosquito species) of WNV should be studied to better understand the relation between WNV transmission and environmental conditions. To achieve this, entomological and ornithological field studies are required to identify the potential hosts and vectors of WNV in different European locations. Serological surveys on wild avifauna in European and African countries are also needed to confirm and better understand the importance of the type of passerine flyways (Eastern or Western pathway) to predict the risk of WND occurrence in Europe. Expansion of this work should incorporate entomological and ornithological information to improve our understanding of WNV transmission cycles, and model predictions.

Environmental indices derived from satellite remote sensing have been widely used to study environmental risk factors of vector-borne diseases in the last decades[32, 33]. Our study also demonstrates that remote sensing medium spatial resolution imagery as MODIS can provide relevant environmental indicators for studies on WNV risk, as it was also recently shown in Northern America[34]. This analysis should include high resolution geolocation of WND cases (i.e. sub-district level) to better address relationship between outbreaks and monitoring of vegetation and water resources.

As in most epidemiological studies, our analysis is dependent on the quality of the original dataset, which suffers from biases related to possible under-detection and under-reporting of WND outbreaks. Indeed, presence prediction can be wrongly interpreted as false positives due to under reporting by national surveillance systems. To limit these biases of possible heterogeneity in space and time of WND surveillance systems, we considered only presence/absence data per district instead of incidence data, and a bootstrap approach was applied to construct confidence intervals. In the future, harmonization of WND surveillance systems across European member states will contribute to improve the epidemiological datasets for analysis.

Conclusion

This study contributes to a better understanding of the environmental and climatic drivers of WNV transmission in Europe and provides the basis for further integration of environmental information together with WNV surveillance programs in Europe, including surveillance of human cases, serological surveillance of domestic and wild avifauna, and entomological surveillance, as recommended by[13, 35]. Our results suggest that risk maps for WNV transmission could be built based on updated anomalies of temperatures and MNDWI. Use of those two environmental risk factors would allow for a better integration of environmental monitoring into public health surveillance systems.

Methods

Epidemiological data

Data series of confirmed cases of WND in humans reported in Europe and its neighboring countries between 2002 and 2013 were compiled at the district level to define presence of human case (Table 1, Figure 5). The WND data was assembled through a systematic review of the scientific literature using the MEDLINE database, Embase, and Scopus. Abstracts were screened to identify epidemiological information and the full text articles were retrieved to obtain the precise location of WND cases. In addition, a systematic screening was conducted using the Global Infectious Diseases and Epidemiology Online Network, Promed for central Europe and meeting reports and newsletters of scientific projects on WND of the Framework Programme for Research and Technological Development funded by European union[36].

Figure 5
figure 5

Study area: location of the West Nile fever outbreaks, Europe and neighbouring countries, 2002–2011, and extent of environmental datasets. Note: the brown box represents the extent of temperature dataset; the green boxes represent the extent of the MODIS tiles.

The screening described above was merged with West Nile human cases on West Nile fever surveillance conducted by ECDC during the transmission season in Europe[11].

Population data

Within European Union, the nomenclature of territorial units for statistics classification (NUTS) was used for administrative unit at level 3 with population estimate of 2010[37]. For outside of the European Union, administrative units data were provided by the Global Administrative Unit Layers (GAUL) project and population estimates for the year 2010 were derived from the Gridded Population of the World (GPW) dataset[38].

Ecological and climatic data birds’ migration routes

Passerine fly ways were digitized in order to categorize administrative units in two categories of migration flyway (western and eastern) according to the migration flyways of Western Palearctic Passerines South Eastern European bird migration network[39].

Wetlands

The presence or absence of wetlands in a district was defined according to Ramsar Sites Information Service form the Convention on Wetlands of International Importance[40]. Quality control was made by image interpretation using Google Earth and the Global Lakes and Wetlands Database (GLWD, level data 1 and 2) maintained by the World Wide Fund for Nature and the Center for Environmental Systems Research, University of Kassel, Germany[41].

Time-dependent variables

We examined three time-dependent environmental variables as predictors of WNV risk: the temperature, an index related to the state of the vegetation, the NDVI, and an index related to the state of the water bodies, the MNDWI. The two latter were derived from MODIS imagery. These three indices vary seasonally and inter-annually, and are hypothesized to influence mosquito population dynamics, and WNV transmission[31, 34, 42].

Temperature has been found to be associated with WND outbreaks in the Old[12] and New World[34]. Indeed, elevated temperatures favor vector competence, accelerate mosquito development and reproduction rates, thus influencing mosquito-to-vertebrate transmission rates[43, 44]. The gridded data of monthly mean of air temperature for the region between 30°N-60°N and 10°W-55°E was obtained from the NOAA NCEP-NCAR database[45, 46] for each month from January 1981 to December 2013, (Table 3, Figure 5).

Table 3 Characteristics of the environmental datasets

On the other hand, vegetation indices such as the NDVI have been identified as risk factors for WND outbreaks occurrence in previous North American studies[34, 47]. NDVI may serve as an indicator of environmental conditions suitable for vegetation growth and emergence of mosquito populations[34].

The presence of water bodies was identified as another environmental WND risk factor[26], because large standing water resources may lead to an upsurge of mosquito populations. Several methods are used in radar and optical remote sensing to delineate water bodies and map flooded areas. The Modified Normalized Difference Water Index (MNDWI) is particularly suited to the detection of free water[48].

To derive NDVI and MNDWI values, MODIS data products were acquired from Land Process Distributed Active Archive Center (LP DAAC). MODIS Terra 8-day composite images of surface reflectance estimates at 500 m spatial resolution (product MOD09A1) were acquired for all WNV infected countries for a twelve years period (2002–2013) (Table 3, Figure 5). Preprocessing steps consisted in masking the low quality pixels and the pixels covered by clouds using the surface reflectance quality file, and performing a linear temporal interpolation of the masked pixels using the Time Series Generator (TiSeG) freeware[49].

NDVI and MNDWI temporal series were computed from the reflectance values of the cleaned images according to Equations 1 and 2:

NDVI = NIR − R / NIR + R
(1)
MNDWI = G − SWIR / G + SWIR
(2)

with NIR: reflectance in near infrared range, R: reflectance in red wavelength; G: reflectance in green wavelength; SWIR: reflectance in short wave infrared range.

Image processing was performed using ENVI IDL software 4.8 (Exelis, Boulder, C0, USA).

Quantifying anomalies of temperature, vegetation and water indices relative to the long term average

First, the long-term average and standard deviation of each of the environmental indices were computed on monthly bases for the temperature data, and on 8-days interval bases for the MODIS NDVI and MNDWI data. The long-term statistics were computed on a 20 years period (1981–2010) and on a 10 years period (2002–2011) for the temperatures and for the MODIS-derived indices, respectively.

Second, the anomaly (z) of temperatures, NDVI and MNDWI was calculated for each date i (month or 8-days period) as a function of the annual indices x i and their long-term average and standard deviation values (Equation 3). Study period was selected to begin in March (before WND outbreaks) and end in August.

z i = x i − x ¯ i σ i
(3)

Finally, the mean anomalies of temperatures, NDVI and MNDWI were computed for each district and each month and MODIS 8-days period.

Analyses were performed using ArcGIS 10.1 and Spatial analyst extension (ESRI, Redlands, CA).

Statistical analysis

The status of infection by WNV in relationship to ecological and climatic variables was analyzed at the district level (n = 1113 spatial units). For each year, each district was categorized as ‘infected’ if WND human cases were reported there that year, and as ‘non infected’ otherwise. We used the data from 2002 to 2011 to fit the models, and the data for 2012 and 2013 for external validation.

The probability of a district to be infected by WNV was assessed using logistic regression models, with the status of infection as the response variable, and as explanatory variables the population, the presence of wetlands, the presence of birds’ migratory routes, the anomalies of temperature, NDVI and MNDWI. We also tested as explanatory variable the occurrence of a WND outbreak the previous year, considering that WNV could persist locally through survival in overwintering mosquitoes[25]: for each year i and each district j, a synthetic index (λ) of WND outbreak occurrence in the neighborhood of district j was defined as the weighted average of the number of infected districts amongst a set of n j neighbours of district j according to Equation 4:

λ i , j = ∑ k = 1 n j y i − 1 , k . w jk n j + 1
(4)

with yi,j: infected status of district j, year i (0: ‘non infected’; 1: ‘infected’); n j : number of neighbouring districts of district j; the weight w jk given to district k is 1 if k ≠ j and 2 otherwise.

Univariate analysis. Explanatory variables were tested one at a time to test associations with district WNV status. Significant variables in this preliminary univariate screening analysis at 0.05 p-value were kept for analysis of co-linearity.

Multivariate analysis. Multivariate logistic models were built to examine the role of explanatory variables, having adjusted for other variables. All possible models including variables significant in the univariate analysis were fitted. In case of co-linearity, the variable corresponding to the earlier date was kept. We used backward model selection based on AIC to select the best model based on both model fit and model parsimony[50]. A bootstrap procedure (1,000 replicates) was applied to estimate the 95% confidence interval (95% CI) of the logistic regression models’ coefficients, selecting randomly each time from the original set of 1113 districts 90% of infected districts between 2002 and 2011 (n = 98) and 90% of non-infected districts (n = 903).

Validation. The predictive accuracy of the final model was assessed using the ROC (Receiver Operating Characteristic) curve[51]. The area under curve (AUC) of the ROC curve, sensitivity and specificity indices were computed using data from 2002 to 2011 (internal validation) and data from 2012 and 2013 (external validation). The greater the AUC, the closer the predictions are to the observed data. 95% CI were estimated using the same bootstrap procedure described above (1,000 replicates).

Generation of predictive maps. The final model was used to predict the probability of WNV infection for the entire Europe and neighboring countries by applying the model to observed data from all districts of the studied area from 2002 to 2013.

Analyses were performed using the ‘R’ freeware and additional ‘spdep’, ‘lme4’, and ‘PresenceAbsence’ packages[52].

Abbreviations

AIC:

Akaike information criterion

AUC:

Area under curve

CI:

Confidence interval

GAUL:

Global administrative unit layers

GPW:

Gridded population of the world

GLWD:

Google earth and the global lakes and wetlands database

NDVI:

Normalized difference vegetation index

NUTS:

Nomenclature of territorial units for statistics classification

MODIS:

Moderate resolution imaging spectroradiometer

MNDWI:

Modified normalized difference water index

ROC curve:

Receiver operator characteristic curve

WND:

West nile disease

WNND:

West nile neuro invasive disease

WNF:

West nile fever

WNV:

West nile virus.

References

  1. Lanciotti RS, Roehrig JT, Deubel V, Smith J, Parker M, Steele K, Crise B, Volpe KE, Crabtree MB, Scherret JH, Hall RA, MacKenzie JS, Cropp CB, Panigrahy B, Ostlund E, Schmitt B, Malkinson M, Banet C, Weissman J, Komar N, Savage HM, Stone W, McNamara T, Gubler DJ: Origin of the West Nile virus responsible for an outbreak of encephalitis in the northeastern United States. Science. 1999, 286 (5448): 2333-2337.

    Article  CAS  PubMed  Google Scholar 

  2. Owen J, Moore F, Panella N, Edwards E, Bru R, Hughes M, Komar N: Migrating birds as dispersal vehicles for West Nile virus. Ecohealth. 2006, 3 (2): 79-85.

    Article  Google Scholar 

  3. Campbell GL, Marfin AA, Lanciotti RS, Gubler DJ: West Nile virus. Lancet Infect Dis. 2002, 2 (9): 519-529.

    Article  PubMed  Google Scholar 

  4. Dauphin G, Zientara S, Zeller H, Murgue B: West Nile: worldwide current situation in animals and humans. Comp Immunol Microbiol Infect Dis. 2004, 27 (5): 343-355.

    Article  PubMed  Google Scholar 

  5. Paz S, Semenza JC: Environmental drivers of West Nile fever epidemiology in Europe and Western Asia–a review. Int J Environ Res Public Health. 2013, 10 (8): 3543-3562.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Tber Abdelhaq A: West Nile fever in horses in Morocco. Bulletin De l'Office International des Epizooties. 1996, 108: 867-869.

    Google Scholar 

  7. Autorino GL, Battisti A, Deubel V, Ferrari G, Forletta R, Giovannini A, Lelli R, Murri S, Scicluna MT: West Nile virus epidemic in horses, Tuscany region Italy. Emerg Infect Dis. 2002, 8 (12): 1372-1378.

    Article  PubMed Central  PubMed  Google Scholar 

  8. Savini G, Monaco F, Calistri P, Lelli R: Phylogenetic analysis of West Nile virus isolated in Italy in 2008. Euro surveill. 2009, 13 (48): 1-2.

    Google Scholar 

  9. Murgue B, Murri S, Zientara S, Durand B, Durand JP, Zeller H: West Nile outbreak in horses in Southern France, 2000: the return after 35 years. Emerg Infect Dis. 2001, 7 (4): 692-696.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Hubalek Z, Halouzka J, Juricova Z: West Nile fever in Czechland. Emerg Infect Dis. 1999, 5 (4): 594-595.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. European Center for Disese Prevention and Control (ECDC): West Nile fever maps. [http://www.ecdc.europa.eu/en/healthtopics/west_nile_fever/west-nile-fever-maps/pages/index.aspx]. Updated: November 6, 2013. Updated: November 6, 2013

  12. Paz S, Malkinson D, Green MS, Tsioni G, Papa A, Danis K, Sirbu A, Ceianu C, Katalin K, Ferenczi E, Zeller H, Semenza JC: Permissive summer temperatures of the 2010 European West Nile fever upsurge. PLoS One. 2013, 8 (2): e56398-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Ozdenerol E, Taff GN, Akkus C: Exploring the spatio-temporal dynamics of reservoir hosts, vectors, and human hosts of West Nile virus: a review of the recent literature. Int J Environ Res Public Health. 2013, 10 (11): 5399-5432.

    Article  PubMed Central  PubMed  Google Scholar 

  14. May FJ, Davis CT, Tesh RB, Barrett AD: Phylogeography of West Nile virus: from the cradle of evolution in Africa to Eurasia, Australia, and the Americas. J Virol. 2011, 85 (6): 2964-2974.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Soverow JE, Wellenius GA, Fisman DN, Mittleman MA: Infectious disease in a warming world: how weather influenced West Nile virus in the United States (2001–2005). Environ Health Perspect. 2009, 117 (7): 1049-1052.

    Article  PubMed Central  PubMed  Google Scholar 

  16. Gibbs SEJ, Wimberly MC, Madden M, Masour J, Yabsley MJ, Stallknecht DE: Factors affecting the geographic distribution of West Nile virus in Georgia, USA: 2002–2004. Vector Borne Zoonot. 2006, 6 (1): 73-82.

    Article  Google Scholar 

  17. Reisen WK, Fang Y, Martinez VM: Effects of temperature on the transmission of West Nile virus by Culex tarsalis (Diptera : Culicidae). J Med Entomol. 2006, 43 (2): 309-317.

    Article  PubMed  Google Scholar 

  18. Kunkel KE, Novak RJ, Lampman RL, Gu W: Modeling the impact of variable climatic factors on the crossover of culex restauns and Culex pipiens (Diptera: culicidae), vectors of West Nile virus in Illinois. Am J Trop Med Hyg. 2006, 74 (1): 168-173.

    PubMed  Google Scholar 

  19. Liu A, Lee V, Galusha D, Slade MD, Diuk-Wasser M, Andreadis T, Scotch M, Rabinowitz PM: Risk factors for human infection with West Nile virus in connecticut: a multi-year analysis. Int J Health Geogr. 2009, 8: 67-

    Article  PubMed Central  PubMed  Google Scholar 

  20. Chen CC, Jenkins E, Epp T, Waldner C, Curry PS, Soos C: Climate change and West Nile virus in a highly endemic region of North America. Int J Environ Res Public Health. 2013, 10 (7): 3052-3071.

    Article  PubMed Central  PubMed  Google Scholar 

  21. Reisen WK, Fang Y, Lothrop HD, Martinez VM, Wilson J, Oconnor P, Carney R, Cahoon-Young B, Shafii M, Brault AC: Overwintering of West Nile virus in Southern California. J Med Entomol. 2006, 43 (2): 344-355.

    Article  PubMed  Google Scholar 

  22. Shaman J, Day JF, Stieglitz M: Drought-induced amplification and epidemic transmission of West Nile virus in southern Florida. J Med Entomol. 2005, 42 (2): 134-141.

    Article  PubMed  Google Scholar 

  23. Shelite TR, Rogers CM, Litzner BR, Johnson RR, Schneegurt MA: West Nile virus antibodies in permanent resident and overwintering migrant birds in south-central Kansas. Vector Borne Zoonotic Dis. 2008, 8 (3): 321-329.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Tesh RB, Parsons R, Siirin M, Randle Y, Sargent C, Guzman H, Wuithiranyagool T, Higgs S, Vanlandingham DL, Bala AA, Haas K, Brian Zerinque B: Year-round West Nile virus activity, Gulf Coast region Texas and Louisiana. Emerg Infect Dis. 2004, 10 (9): 1649-1652.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Monaco F, Savini G, Calistri P, Polci A, Pinoni C, Bruno R, Lelli R: 2009 West Nile disease epidemic in Italy: first evidence of overwintering in Western Europe?. Res Vet Sci. 2011, 91 (2): 321-326.

    Article  CAS  PubMed  Google Scholar 

  26. Leblond A, Sandoz A, Lefebvre G, Zeller H, Bicout DJ: Remote sensing based identification of environmental risk factors associated with West Nile disease in horses in Camargue France. Prev Vet Med. 2007, 79 (1): 20-31.

    Article  CAS  PubMed  Google Scholar 

  27. Chase JM, Knight TM: Drought-induced mosquito outbreaks in wetlands. Ecol Lett. 2003, 6 (11): 1017-1024.

    Article  Google Scholar 

  28. Paz S, Malkinson D, Green MS, Tsioni G, Papa A, Danis K, Sirbu A, Ceianu C, Krisztalovics K, CsohĂ¡n A, EmÅ‘ke F, Zeller H, Semenza J: West Nile Fever Upsurge in Europe and its Neighbouring Countries and its Possible Association with Elevated Temperatures. Proceedings of the GRF One Health Summit: 2012. 2012, Switzerland: Davos, 170-173.

    Google Scholar 

  29. Frost MJ, Zhang J, Edmonds JH, Prow NA, Gu X, Davis R, Hornitzky C, Arzey KE, Finlaison D, Hick P, Read A, Hobson-Peters J, May FJ, Doggett SL, Haniotis J, Russell RC, Hall RA, Khromykh AA, Kirkland PD: Characterization of virulent West Nile virus Kunjin strain, Australia, 2011. Emerg Infect Dis. 2012, 18 (5): 792-800.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Shaman J, Stieglitz M, Stark C, Le Blancq S, Cane M: Using a dynamic hydrology model to predict mosquito abundances in flood and swamp water. Emerg Infect Dis. 2002, 8 (1): 6-13.

    Article  PubMed  Google Scholar 

  31. DeGroote JP, Sugumaran R, Brend SM, Tucker BJ, Bartholomay LC: Landscape, demographic, entomological, and climatic associations with human disease incidence of West Nile virus in the state of Iowa USA. Int J Health Geogr. 2008, 7: 19-

    Article  PubMed Central  PubMed  Google Scholar 

  32. Kalluri S, Gilruth P, Rogers D, Szczur M: Surveillance of arthropod vector-borne infectious diseases using remote sensing techniques: A review. Plos Pathogens. 2007, 3 (10): 1361-1371.

    Article  CAS  PubMed  Google Scholar 

  33. Lambin EF, Tran A, Vanwambeke SO, Linard C, Soti V: Pathogenic landscapes: interactions between land, people, disease vectors, and their animal hosts. Int J Health Geogr. 2010, 9: 54-

    Article  PubMed Central  PubMed  Google Scholar 

  34. Chuang TW, Wimberly MC: Remote Sensing of Climatic Anomalies and West Nile Virus Incidence in the Northern Great Plains of the United States. Plos One. 2012, 7 (10): e46882-

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Kwan JL, Park BK, Carpenter TE, Ngo V, Civen R, Reisen WK: Comparison of enzootic risk measures for predicting West Nile disease, Los Angeles, California, USA, 2004–2010. Emerg Infect Dis. 2012, 18 (8): 1298-1306.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Global infections diseases and epidemiology online network (GIDEON). [http://www.gideononline.com/]. Updated: May 6, 2014. Updated: May 6, 2014

  37. EUROSTAT: Nomenclature of territorial units for statistics. 2012, [http://epp.eurostat.ec.europa.eu/portal/page/portal/nuts_nomenclature/introduction]. Updated: November 9, 2012. Updated: November 9, 2012.

    Google Scholar 

  38. NASA, Socioeconomic Data and Application Centre (SEDAC): Gridded population of the world (gpw), v3. [http://sedac.ciesin.columbia.edu/data/set/gpw-v3-population-density]

  39. SE European Bird Migration Network: Migration flyways of western paleartic passerines. [http://www.seen-net.eu/]

  40. The Ramsar Convention of Wetlands: Ramsar information information service. [http://www.ramsar.org]

  41. Global Lakes and Wetlands Database. [http://www.ramsar.org]

  42. Calistri P, Ippoliti C, Candeloro L, Benjelloun A, El Harrak M, Bouchra B, Danzetta ML, Di Sabatino D, Conte A: Analysis of climatic and environmental variables associated with the occurrence of West Nile virus in Morocco. Prev Vet Med. 2013, 110 (3–4): 549-553.

    Article  PubMed  Google Scholar 

  43. Epstein PR: West Nile virus and the climate. J Urban Health. 2001, 78 (2): 367-371.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Kilpatrick AM, Meola MA, Moudy RM, Kramer LD: Temperature, viral genetics, and the transmission of West Nile virus by Culex pipiens mosquitoes. PLoS Pathog. 2008, 4 (6): e1000092-

    Article  PubMed Central  PubMed  Google Scholar 

  45. Kalnay E, Kanamitsu M, Kistler R, Collins W, Deaven D, Gandin L: The ncep/ncar 40-year reanalysis project. Bull Am Meteorol Soc. 1996, 77: 437-471.

    Article  Google Scholar 

  46. NOAA: NOAANCEP-NCAR CDAS-1 Monthly Diagnostic Surface temp: Temperature Data. [http://iridl.ldeo.columbia.edu/SOURCES/.NOAA/.NCEP-NCAR/.CDAS-1/.MONTHLY/.Diagnostic/.surface/.temp/]. Updated: May 2, 2014. Updated: May 2, 2014

  47. Ward MP: Equine West Nile virus disease occurrence and the Normalized Difference Vegetation Index. Prev Vet Med. 2009, 88 (3): 205-212.

    Article  PubMed  Google Scholar 

  48. Xu HQ: Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery. International Journal of Remote Sensing. 2006, 27 (14): 3025-3033.

    Article  Google Scholar 

  49. Colditz RR, Conrad C, Wehrmann T, Schmidt M, Dech S: TiSeG: A Flexible Software Tool for Time-Series Generation of MODIS Data Utilizing the Quality Assessment Science Data Set. Ieee T Geosci Remote. 2008, 46 (10): 3296-3308.

    Article  Google Scholar 

  50. Burnham K, Anderson D: Multi model inference understanding AIC and BIC in model selection. Sociological Methods & Res. 2004, 33 (2): 261-304.

    Article  Google Scholar 

  51. Park SH, Goo JM, Jo CH: Receiver operating characteristic (ROC) curve: practical review for radiologists. Korean J Radiol. 2004, 5 (1): 11-18.

    Article  PubMed Central  PubMed  Google Scholar 

  52. The R Development Core Team: R: A language and environment for statistical computing. 2008, Vienna: Austria: R: Foundation for Statistical Computing

    Google Scholar 

Download references

Acknowledgements

This study was financed by ECDC (OJ/2012/02/16 – PROC/2012/015). We are grateful to the support from G. Balança, N. Gaidet, M. Gely, T. Baldet, and T. Balenghien from CIRAD, J. Figuerola from EstaciĂ³n BiolĂ³gica Doñana, CSIC, L. Marrama, and A. Donachie from ECDC, and Pshem Busse from the Eastern European bird migration network for migration flyways information.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jan C Semenza.

Additional information

Competing interests

Competing financial interests declaration: none declared.

Authors’ contributions

JCS initiated the project entitled: Mosquito-Borne Diseases Determinants (OJ/2012/02/16 – PROC/2012/015) and developed the strategy for a NWF analysis. BS compiled the epidemiologic data and AT led the analysis with contributions from BS, SP, MR, AD, VC and JCS. All authors read and approved the final manuscript.

Electronic supplementary material

12942_2014_598_MOESM1_ESM.doc

Additional file 1: Table S1: Results of univariate analysis: significant variables in univariate screening analysis at 0.05 p-value. (DOC 35 KB)

12942_2014_598_MOESM2_ESM.zip

Additional file 2: Figure S1: Maps of predicted probability of WNV infection based on environmental predictors, and West Nile fever outbreaks, Europe and neighbouring countries, 2002–2013. (ZIP 14 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tran, A., Sudre, B., Paz, S. et al. Environmental predictors of West Nile fever risk in Europe. Int J Health Geogr 13, 26 (2014). https://doi.org/10.1186/1476-072X-13-26

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1476-072X-13-26

Keywords