Skip to main content

Climate change projections of West Nile virus infections in Europe: implications for blood safety practices

Abstract

Background

West Nile virus (WNV) is transmitted by mosquitoes in both urban as well as in rural environments and can be pathogenic in birds, horses and humans. Extrinsic factors such as temperature and land use are determinants of WNV outbreaks in Europe, along with intrinsic factors of the vector and virus.

Methods

With a multivariate model for WNV transmission we computed the probability of WNV infection in 2014, with July 2014 temperature anomalies. We applied the July temperature anomalies under the balanced A1B climate change scenario (mix of all energy sources, fossil and non-fossil) for 2025 and 2050 to model and project the risk of WNV infection in the future. Since asymptomatic infections are common in humans (which can result in the contamination of the donated blood) we estimated the predictive prevalence of WNV infections in the blood donor population.

Results

External validation of the probability model with 2014 cases indicated good prediction, based on an Area Under Curve (AUC) of 0.871 (SD = 0.032), on the Receiver Operating Characteristic Curve (ROC). The climate change projections for 2025 reveal a higher probability of WNV infection particularly at the edges of the current transmission areas (for example in Eastern Croatia, Northeastern and Northwestern Turkey) and an even further expansion in 2050. The prevalence of infection in (blood donor) populations in the outbreak-affected districts is expected to expand in the future.

Conclusions

Predictive modelling of environmental and climatic drivers of WNV can be a valuable tool for public health practice. It can help delineate districts at risk for future transmission. These areas can be subjected to integrated disease and vector surveillance, outreach to the public and health care providers, implementation of personal protective measures, screening of blood donors, and vector abatement activities.

Background

West Nile virus (WNV) is responsible for large outbreaks of fatal neuroinvasive disease worldwide [1]. WNV infections occur predominantly through mosquito bites but also through blood transfusion or organ, tissue and cell-transplantations. Most human infections are asymptomatic and mild cases present only flu-like symptoms; more severe cases present with signs of encephalitis, meningo-encephalitis or meningitis. Globally, WNV is a widespread arthropod-borne virus, an enveloped, single-strand RNA virus of the genus Flavivirus in the family of Flaviviridae [2]–[4]. It circulates in Africa, Americas, Asia, Europe, and Australia, where it is thought have been introduced from the Middle East [4]. In Europe, Middle East and Africa, WNV has been responsible for sporadic outbreaks in the 1950s in Israel, in the 1960s in Russia and France, in the 1970s in Belarus, South Africa, and Ukraine. However, since 1996 it has caused more recurrent outbreaks in Europe and northern Africa [5]–[8]. In 2010, large outbreaks in humans occurred in Southeastern and Eastern Europe. The European outbreaks occurred in Russia, the Czech Republic, Hungary, Romania, Turkey, Greece, Italy, France, Spain and Portugal [9]. Since 2010, there have been annual outbreaks in Southeastern and Eastern Europe, suggesting an endemic transmission cycle and thus a resurgent public health problem [10].

The transmission cycle of WNV is maintained in nature in an enzootic cycle between susceptible birds and competent ornithophilic mosquitos. A number of intrinsic and extrinsic factors are responsible for the complex nature of WNV epidemiology. Intrinsic factors include vector and host competence, genetic traits, mosquito feeding rates and preferences, mosquito longevity and host immunity [11]. WNV has been detected in over 60 species of mosquitoes but in Europe Culex pipiens seems to play the most important role [12], [13]. These mosquitoes can act as bridge-vectors that feed both on birds as the principal host as well as on incidental hosts such as humans and horses. Local birds are responsible for WNV amplification while infected migratory birds account for the long-distance dispersion of WNV. The extrinsic factors include abiotic conditions such as the environmental settings that shape the transmission dynamic but also biotic factors such as the density and composition of hosts and vectors [9].

In Europe, there are sylvatic and synanthropic transmission scenarios, in rural and urban areas. They respectively (and independently) contribute to high concentrations of hosts with competent mosquito vectors that support intense local avian transmission. Rural areas with estuaries, wetlands or marshes attract migratory birds for breeding, nesting, and rearing their young. These bird habitats also attract bird-feeding mosquitoes where congregating bird populations can get infected with WNV. In fact, a number of human outbreaks have originated in European estuaries such as in the Danube delta in Romania, in the Volga delta in Russia, and in the Rhone delta in France [14]. Urbanized areas can also attract large bird and mosquito populations where humans can get exposed [15]. Moreover, urban infrastructure can serve as mosquito breeding sites, such as pools in ineffective drainage systems [16]. For example, in 1996, a large outbreak occurred in the city of Bucharest, Romania that affected 4 % of the population [7]. Over 800 patients were hospitalized during another large WNV outbreak that occurred in 1999 in Volgograd City on the west bank of the Great Volga River [17] and in 2010 WNV infections were documented for the first time in humans in the Greek city of Thessaloniki [18], [19]. Certain metropolitan areas can sustain a high breeding density of birds; for example, European Starlings thrive on urban lawns or parks where they can feed, or gulls proliferate near open water [20], [21]. Birds from these types of urban environments harbour viruses with higher genetic diversity than birds from residential areas, indicating that anthropogenic factors associated with urbanization play an important role in arboviral transmission and evolution [22]. In the United States, WNV infection rates in crows and humans are higher in more urbanized environments that are less forested [23], [24], while in Europe, WNV transmission can occur both in rural and urban areas.

Ambient temperature is another important environmental determinant in the transmission of WNV as it has a direct impact on mosquito survival, developmental rates of immature stages, growth rates of vector populations and decrease in the interval between blood meals [25]–[28]. Moreover, temperature also affects the extrinsic incubation period (the number of days from ingestion to transmission) by influencing the viral replication rates and thus the transmission of WNV [29], [30]. In a modelling study elevated air temperature was the strongest predictor of increased infection in mosquito vectors [28]. As for Europe, it was found that the unprecedented upsurge in the number of human WNF cases in the summer of 2010 was preceded by extremely hot spells in the region [25]. Moreover, recent research analysed the status of infection by WNV in Europe and its neighbouring countries in relation to environmental and climatic risk parameters. The anomalies of temperature in July were identified as one of the main risk factors [10]. During recent decades, parts of Europe have warmed up more than the global average. Additionally, more frequent and more intense hot extremes have occurred. This trend is expected to continue while predictions suggest a further temperature increase (between 1.0 °C and 5.5 °C) by the end of the century [31]. Increases in ambient temperature due to climate change are therefore projected to impact WNV transmission in Europe and its neighbouring countries [32]–[36].

In order to examine these land-use and climatic variables as predictors of the probability of WNV infection [4] we tested the contribution of remotely sensed temperature, the state of vegetation and water bodies, and bird migratory routes in a statistical model. We also project the WNV risk in Europe into 2025 and 2050, with July temperature projections under a balanced IPCC climate change scenario (A1B) [37]. Among IPCC scenarios that cover a wide range of the main demographic, economic, and technological driving forces of greenhouse gases and sulphur emissions, the A1 scenario groups are distinguished by their technological emphasis. A1B represent a balance across all energy sources (intensive fossil and non-fossil energy). Balanced is defined as not relying too heavily on one particular energy source on the assumption that similar improvement rates apply to all energy supply and end use technologies [37].

Insights from these analyses can also be used to assess the current and future WNV risk to the safety of blood supply in the region [38]. Thus, we provide a quantification of the risk of WNV transmission through blood transfusion by estimating the prevalence of infection in (donor) populations of an area affected by a WNV outbreak, as well as the infection risk of a blood donor that visited an outbreak affected area. In the long-run, the environmental determinants identified in this model lend themselves for an integration of environmental monitoring in public health surveillance systems of human cases, serological surveillance of domestic and wild avifauna, and entomological surveillance [4], [39], [40].

Methods

Temperature data

The association between WNV outbreaks in Europe and a number of temperature parameters such as mean, minimum and maximum and their anomalies were examined by Paz et al. (2013) [25]. Positive temperature anomalies were found to be a major risk factor for WNV occurrence. Subsequently, in a study by Tran et al. (2014) remotely sensed mean temperature anomalies of July were identified as a risk factor for WNV transmission in Europe [10]. Based on these observations temperature anomalies were modelled for the present analysis. Monthly anomalies of July 2014 temperatures at the locations of WNF outbreaks reported in humans were computed from NOAA NCEP-NCAR (US National Centre for Atmospheric Research) database. Predicted surface temperatures and anomalies were extracted for 2015–2050 climate change A1B scenario, all for the gridded region 30°N-60°N and 10°W-55°E [41]. The predicted temperatures were based on the Community Climate System Model (CCSM3) of NCAR, generated on a Gaussian grid. Following NCAR methodology, the temperature anomalies data in the current study are relative to 20th Century Experiment 1980–1999 [36]. The A1B Scenario was chosen since it is a balanced scenario; its main characteristics include: low population growth, very high GDP growth, very high energy use, low-medium land use changes, medium resource (mainly oil and gas) availability, rapid pace and direction of technological change favouring balanced development. The chosen model output is the Ensemble Average which is the mean state of the climate among all model runs. Global climate simulations were produced at NCAR by the CCSM3 for the 4th Assessment report of the IPCC [42].

WNV epidemiology

The methods for the WNV epidemiological model in Europe have been described previously [10]. Briefly, the epidemiologic data of human West Nile cases were obtained from the West Nile fever surveillance conducted at ECDC during the transmission season in Europe [43]. The transmission season in Europe starts in May and extends into November. Population data for the European Union were based on the nomenclature of territorial units for statistics classification (NUTS) at level 3 with population estimate of 2010 [44]. Global Administrative Unit Layers (GAUL) were used for regions outside of the European Union, and project and population estimates for the year 2010 were derived from the Gridded Population of the World (GPW) dataset [45].

Environmental variables

To derive the Normalized difference vegetation index (NDVI) and Modified normalized difference water index (MNDWI) values, Moderate resolution imaging spectroradiometer (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).

Remotely sensed vegetation indices such as the NDVI have been identified as risk factors for WNV outbreaks occurrence in previous North American studies [46], [47]. Indeed, NDVI may serve as an indicator of environmental conditions suitable for vegetation growth and emergence of mosquito populations [46]. The presence of water bodies was identified as another environmental WNV risk factor [47], because large standing water resources may lead to an upsurge of mosquito populations. The Modified Normalized Difference Water Index (MNDWI) is particularly suited to the detection of free water [48]. This index was found to be more appropriate for the analysis than precipitation due to the inconsistent association between precipitation and WNV outbreaks [36]. The presence or absence of wetlands in a district was defined according to Ramsar Sites Information (https://rsis.ramsar.org/) and the Global Lakes and Wetlands Database (GLWD, level data 1 and 2) (http://www.worldwildlife.org/pages/global-lakes-and-wetlands-database). 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 (http://www.seen-net.eu/).

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 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. The mean anomalies of temperatures, NDVI and MNDWI were computed for each district and each month and MODIS 8-days period. The analysis was performed at the district level (n = 1113) categorized as ‘infected’ if WNF cases in humans were reported there that year, and as ‘non-infected’ otherwise.

Multivariate models

We used multivariate logistic regression models to test the status of a district as ‘infected’ or ‘not infected’ as the response variable. Thus, the probability distribution was assumed to be binomial, and the response variable a logit function. As explanatory variables we used the population, the presence of wetlands, the presence of birds’ migratory routes, the anomalies of temperature, NDVI and MNDWI as previously described [10]. We also tested as explanatory variable: the occurrence of a WNV outbreak the previous year (λ: weighted average of the number of infected districts amongst the neighbourhood the previous year), considering that WNV could persist locally through survival in overwintering mosquitoes or infected birds. A bootstrap procedure (1,000 replicates) was applied to estimate the 95 % confidence interval (95 % CI) of the parameter estimates, 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). The final model used in the current analysis is described in Table 1. This model was validated using 2012–2013 epidemiologic data that were not used for model constructions [10].

Table 1 Multivariate logistic regression model parameter of the risk of WNV infection at district level, EU and neighbouring countries [10]

The average, standard deviation, and anomalies of the 2014 temperature was computed and applied to the model described above [10]. The results were compared with the actual occurrence of WNV in 2014. Validity of the model was assessed based on the ability of the model to distinguish between districts with and without WNV, using sensitivity, specificity, and Area Under Curve (AUC) of the Receiver Operator Characteristic Curve (ROC).

Climate change projections

The model was used to predict the probability of WNV infection by applying the projected July temperatures for 2025 and 2050. We extracted gridded temperature projections for the A1B Scenario, which is a balanced climate change scenario [42]. July temperature anomalies were entered for each year between 2015 and 2050 in order to compute the probability of WNV infections per district and per year which is one of the parameters in the model (infection the previous year) [10].

Projections of the prevalence of infection

To assess the safety of the blood supply from a WNV infection we used the European Up-Front Risk Assessment Tool (EUFRAT) developed by ECDC [49]. The risk to transfusion recipients in outbreak affected regions can be computed with this tool, quantifying the risk of transmission of emerging pathogens through blood transfusion. The risk can be calculated for any infectious disease for which a number of characterizing parameters are known or, as in this study, estimated (i.e. population size, prevalence, duration of epidemic, proportion of undetected cases, etc.). We calculated the prevalence of infection in (donor) population in affected areas and the probability of the blood donor to become infected after visiting such an area. For the calculation, the absolute number of donors infected with WNV in the outbreak-affected area was used from the notification/surveillance systems reported to ECDC. This number refers to individuals in the population of the outbreak-affected areas that are reported/ notified as cases when they seek health care during the epidemic. We used 60 % as a conservative estimate of the proportion of undetected cases which includes WNV infected individuals that may not develop symptoms, did not seek for health care, or were misdiagnosed with other diseases and were thus not reported as cases. The population size was entered as the population number in the outbreak-affected areas where the notified cases were recorded. The proportion of chronic infection was set to zero. The duration of the epidemic was the length of period since the first case was reported up until the last day of reported cases and was set to 300 days. The prevalence of infection in the (donor) population was calculated based on these parameters for 2014.

Results

Temperature anomalies in Europe for July 2014 are shown in Fig. 1. The locations of 210 probable and confirmed cases of WNV infections in Europe and neighbouring countries for 2014 are illustrated in Fig. 2. An affected area was defined as an area with one or more autochthonous human WNV cases which were recorded in a number of countries: Austria (1); Greece (15); Hungary (11); Italy (24); Romania (23); Bosnia and Herzegovina (13); Israel (17); Palestine (1); Russian Federation (29); and Serbia (76). We used the multivariate logistic regression model developed by Tran et al. (2014) to compute the probability of WNV occurrence based on July temperatures anomalies, the anomaly of the Modified Normalized Difference Water Index (MNDWI) in early June, an outbreak of the previous year, the size of the human population, wetlands and the type of avian flyways [10]. We applied a standard threshold of 0.5 to classify infected (with probability of infection greater than 0.5) and non-infected (with probability of infection lower than 0.5) districts. The probability map with the July 2014 temperature deviation is presented in Fig. 3. Geographic areas of predicted probability of high WNV infection were North-eastern Greece, central Hungary, North-eastern Italy, Eastern Romania, central Serbia, and large areas of Southern Russia. This probability map was compared with the actual occurrence of WNV cases in 2014 (Fig. 2) and a concordance was observed. Validity of the model was assessed based on the ability of the model to distinguish between districts with and without WNV. External validation of the model with the 2014 data indicated good prediction based on the Receiver Operating Characteristic (ROC) Curve (Fig. 4). The Area Under Curve (AUC) of the model reached 0.871 (SD = 0.032). Sensitivity was calculated as the proportion of WNV positive districts correctly identified by the model and yielded 0.236 (SD = 0.074). Specificity was calculated as the proportion of WNV negative districts correctly identified by the model and yielded 0.982 (SD = 0.006).

Fig. 1
figure 1

Temperature anomalies for July 2014

Fig. 2
figure 2

Districts with probable and confirmed cases of West Nile Virus infections, as of 20/11/2014. Note: An affected area is defined as an area with one or more autochthonous human WNV cases (neuro-invasive and non neuro-invasive), meeting laboratory criteria as per EU case definition’ (Directive 2008/426/EC). Probable and confirmed: A probable case is any person meeting the clinical criteria AND with at least one of the following two: − an epidemiological link; − a laboratory test for a probable case. WNV cases by country: Austria (1); Greece (15); Hungary (11); Italy (24); Romania (23); Bosnia and Herzegovina (13); Israel (17); Palestine (1); Russian Federation (29); and Serbia (76)

Fig. 3
figure 3

Predicted probability of districts with West Nile Virus infections for 2014

Fig. 4
figure 4

Receiver operating characteristic (ROC) curve of the probability of West Nile Virus infections in 2014

Probability maps were generated with the A1B temperature projections (Fig. 5). The results reveal a progressive expansion of areas with an elevated probability for WNV infections, particularly at the edges of the transmission areas. For example, the 2025 map reveals a higher probability of WNV infection in Northern Serbia, Central Hungary, North-Eastern Greece, Eastern and Western parts of Romania and Northwestern Turkey. These are also the areas that recorded new districts affected by WNV infections in 2025 compared to 2014 (Fig. 6), which indicates the expansion of the epidemic. According to this analysis, total of 81 and 268 new districts recorded WNV infections for the first time in 2025 and 2050, respectively compared to 2014 that were situated on the perimeter of the transmission areas. In 2050, the area with a higher probability will have expanded even more with a total of 147 and 405 districts being affected in 2025 and 2050, respectively. WNV infections are a significant concern to the safety of the blood supply, because the blood donated by asymptomatic carriers might inadvertently contaminate the blood supply [38]. Therefore, we calculated the prevalence of infection in (donor) populations in the outbreak-affected areas. We mapped the donor population infectivity (Fig. 7) and found an extended area of elevated WNV infection hazard for the safety of the blood supply in 2025 compared to 2014.

Fig. 5
figure 5

Predicted probability of districts with West Nile Virus infections based on July temperatures for A1B scenario projections for 2025 (a) and 2050 (b). Note: Among IPCC scenarios, the A1 scenario groups are distinguished by their technological emphasis. A1B represent a balance across all energy sources (intensive fossil and non-fossil energy)

Fig. 6
figure 6

New districts affected by West Nile Virus infections in 2025 compared to 2014. Note: Panel a: Confirmed: A confirmed case is any person meeting laboratory criteria for case confirmation. Panel b: Total (confirmed and probable): A probable case is any person meeting the clinical criteria AND with at least one of the following two: − an epidemiological link; − a laboratory test for a probable case

Fig. 7
figure 7

Estimated prevalence of West Nile Virus infections in the blood donor population (per 100,000) by districts for 2014 (a) and for 2025 (b). Note: The prevalence of infection in the (donor) population was calculated based on the European Up-Front Risk Assessment Tool (EUFRAT) developed by ECDC [49]

Discussion

Projections of WNV risk

The multivariate model for WNV outbreaks developed for the period 2002 to 2011 and validated with 2012 and 2013 data [10] was applied here to July 2014 temperature anomalies. The probability map revealed a good agreement with the actual WNV outbreaks in 2014, based on the ROC curve. Under the A1B climate change scenario for 2025, this model projects a higher probability of WNV infection in Northern Serbia, Central Hungary, North-Eastern Greece, Eastern and Western parts of Romania, and Northwestern Turkey. These are also the areas where new districts will be affected by WNV infections.

The ability of the model to detect a district with WNV infections was moderate (sensitivity = 0.236) while the ability of the model to correctly detect a negative district was very good (specificity = 0.982). The probability cut off was chose based on public health grounds. For example, the threshold value obtained using the observed prevalence method (p = 0.07) would generate a sensitivity of 0.706 and a specificity of 0.915. Instead we opted for a higher specificity, which implies (given the compromise Se/Sp) a low sensitivity. In that sense, our predictions are rather optimistic. By increasing the sensitivity the number of WNV positive districts would have increased, potentially “overestimating” the number of WNV positive districts. The distribution of WNV positive districts does not follow a Gaussian distribution, but is skewed to the right. Consequently, a small change in the threshold would translate into a significant increase in the number of WNV positive districts even with potentially few WNV cases per district (It is important to note that the model predicts the probability for a district to be positive for WNV, but not the number of cases). In light of limited resources in public health we chose a more conservative approach trying to identify those districts that are predicted to be truly negative for WNV (high specificity).

These findings indicate that July temperature anomalies can be used in our prediction model as an early warning of the coming WNV season and help delineate areas of imminent transmission. Positive temperature anomalies may lead to an upsurge in the growth rates of vector populations, a decrease in the interval between blood meals, a shortening of the incubation time in mosquitoes, virus evolution rate accelerating and an increase of the viral transmission efficiency to birds [9], [29], [36]. These areas can then be targeted for enhanced epidemiologic surveillance of neuro-invasive illness (which may be suggestive of WNV infection), awareness rising among healthcare workers for the clinical presentation of WNV infection, and reinforced laboratory diagnostic capacity. Moreover, these areas can be subjected to entomological surveillance to characterize mosquito breeding sites, and vector abatement measures to diminish mosquito densities. Passive surveillance in domestic birds and equine populations to monitor the dispersion of WNV might also be considered.

This model was also applied to July temperature projections for 2025 and 2050 in order to quantify the future burden of WNV. These results indicate that 81 and 268 districts will be impacted by WNV outbreaks by 2025 and 2050, respectively. Districts adjacent to districts with current transmission are at elevated risk and should therefore also be considered for the public health interventions listed above.

Blood supply safety

The arrival and dispersal of tropical pathogens to Europe and its neighbouring countries commonly associated with warmer temperatures pose a threat to the supply of safe blood products, particularly if they are unknown or without diagnostic tests [38], [50]. Thus, emerging infectious diseases will continue to pose a threat to transfusion safety on European and international levels [51]–[53]. Specifically, a progressive expansion of areas with an elevated probability for WNV infections will increase the threat to the safety of blood transfusion. In geographically larger areas affected by WNV, a higher number of blood donors will be exposed to infection for a longer time period (if the duration of the annual mosquito activity season will be prolonged). Our climate change projections of the predicted probability of WNV infection in Europe have far reaching implications for public health in the future, because the findings in this paper can contribute to WNV preparedness activities. Besides the cases of primary WNV infections, secondary infections through contaminated blood products are of increasing concern to threaten the safety of the blood supply [38]. The asymptomatic blood-borne phase of a WNV infection increases the potential for transmission by transfusion, even if it is relatively short, compared to Hepatitis B virus or HIV [54]. Moreover, WNV has the ability to survive and persist in collected blood and stored blood components and subsequently cause an infection through the intravenous application. Thus, we calculated the prevalence of infection in the donor populations in the outbreak-affected areas (Fig. 7). However, this instantaneous estimate may underestimate the true prevalence of infection if the timing of the WNV epidemic is at the peak of the epidemic curve. Nevertheless, the map reveals considerable vulnerabilities in Southeastern Europe when it comes to the safety of the blood supply. These insights can help transfusion service and clinical staff identify, manage and plan for transfusion-transmitted infections. The overall management of blood safety should be addressed at the institutional level, specifically at regulatory agencies or professional organizations [54]. To preserve the number of eligible donors and an adequate blood supply, authorities will inevitably reassess risk reduction interventions [38]. These include deferral strategies [55], [56], screening strategies and triggers [57]–[59], and also pathogen reduction technologies [38], [60]–[62]. Moreover, it might be necessary to distribute blood components to outbreak areas from unaffected areas in order to prevent intermittent shortfalls in the blood supply [63]. Such coordination requires supra-national inventories of blood products that can be dispatched upon demand. Therefore, the projected temperature change with elevated probability for WNV infections and possible increased prevalence of WNV infected blood donors should be taken into account in developing preparedness plans for the WNV safety of the blood supply. This could substantially increase the costs of blood transfusion therapy. The projections presented here may therefore offer an insight into future developments in the risk of transfusion from WNV infections; it also provides an opportunity to timely define the optimal strategy of blood safety in the face of limited available resources.

Limitations

The projections presented in this paper are computed on the assumption that the other variables remained stable. For example, the population size was maintained constant in these projections; however the contribution of the population size to the model output is very minimal, based on a low parameter value (1.66 10−7). Moreover, the Eastern bird migration path was retained in the model, as well as the presence of wetlands, and the MNDWI21 under the assumption that they will not change significantly. Yet, other WNV determinants that were not part of our model could change over the given time period. For example, mosquito control measures could differ at the area and country level and change over time, with implications for WNV incidence. Similarly, the quality of surveillance clinical case detection might differ too and change over time. WNV might spread at different rates in local, amplifying bird populations that might eventually develop herd immunity. Moreover, other environmental determinates might vary over time as well. We used a statistical approach which involved modelling the 2002–2013 WNV distribution to assess the eco-climatic constraints influencing WNV transmission. The advantage of this statistical approach is that the most significant temperature and environmental conditions for determining WNV presence and absence can be identified. In contrast, biological, mechanistic, or deterministic models compute the basic reproduction rate (R0) with biological thresholds relevant to WNV transmission. However, parameterization of these biological processes has proven to be difficult; moreover, biological models do not explicitly model the socioeconomic and public health contexts in which WNV transmission occurs. Although there are a number of limitations to our model, its validation performed rather well based on the sensitivity and specificity calculations for WNV infections in 2012, 2013 [10] and 2014.

Conclusion

In several countries of Southeastern Europe, WNV transmission is now established. Our predictive model suggests further WNV dispersal in the coming years to adjacent districts. Monitoring and modelling climatic and environmental conditions permissive for the interaction of migratory birds, resident birds, competent mosquito vectors and humans can help delineate districts at risk of transmission. In light of the rapid pace of urbanization internationally [64], areas at risk of current and future WNV transmission can be targeted for integrated surveillance, vector control measures, outreach to the public and health care sector, strengthened laboratory capacity for reliable WNV diagnosis, and systematic screening of blood donors [40]. These activities call for inter-sectorial collaboration to tackle the challenges of WNV transmission.

Additional file

References

  1. Petersen LR, Brault AC, Nasci RS: West Nile virus: review of the literature. JAMA. 2013, 310: 308-315. 10.1001/jama.2013.8042.

    Article  CAS  Google Scholar 

  2. Ciota AT, Kramer LD: Vector-virus interactions and transmission dynamics of West Nile virus. Viruses. 2013, 5: 3021-3047. 10.3390/v5123021.

    Article  Google Scholar 

  3. 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: 2964-2974. 10.1128/JVI.01963-10.

    Article  CAS  Google Scholar 

  4. 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: 5399-5432. 10.3390/ijerph10115399.

    Article  Google Scholar 

  5. Hayes C: West Nile Fever. The Arboviruses: Epidemiology and Ecology. Edited by: Monath TP. 1989, CRC Press, Boca Raton

    Google Scholar 

  6. Karabatsos N: International Catalogue of Arbovirus Including Certain Other Viruses of Vertebrate. 1985, American Society of Tropical Medicine and Hygiene, San Antonio

    Google Scholar 

  7. Tsai TF, Popovici F, Cernescu C, Campbell GL, Nedelcu NI: West Nile encephalitis epidemic in southeastern Romania. Lancet. 1998, 352: 767-771. 10.1016/S0140-6736(98)03538-7.

    Article  CAS  Google Scholar 

  8. Lvov DK, Butenko AM, Gromashevsky VL, Larichev VP, Gaidamovich SY, Vyshemirsky OI, et al: Isolation of two strains of West Nile virus during an outbreak in southern Russia, 1999. Emerg Infect Dis. 2000, 6: 373-376. 10.3201/eid0604.000408.

    Article  CAS  Google Scholar 

  9. 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: 3543-3562. 10.3390/ijerph10083543.

    Article  Google Scholar 

  10. Tran A, Sudre B, Paz S, Rossi M, Desbrosse A, Chevalier V, et al: Environmental predictors of West Nile fever risk in Europe. Int J Health Geogr. 2014, 13: 26-10.1186/1476-072X-13-26.

    Article  Google Scholar 

  11. Hardy JL, Houk EJ, Kramer LD, Reeves WC: Intrinsic factors affecting vector competence of mosquitoes for arboviruses. Annu Rev Entomol. 1983, 28: 229-262. 10.1146/annurev.en.28.010183.001305.

    Article  CAS  Google Scholar 

  12. Gomes B, Sousa CA, Vicente JL, Pinho L, Calderon I, Arez E, et al: Feeding patterns of molestus and pipiens forms of Culex pipiens (Diptera: Culicidae) in a region of high hybridization. Parasit Vectors. 2013, 6: 93-10.1186/1756-3305-6-93.

    Article  Google Scholar 

  13. Molaei G, Andreadis TG, Armstrong PM, Anderson JF, Vossbrinck CR: Host feeding patterns of Culex mosquitoes and West Nile virus transmission, northeastern United States. Emerg Infect Dis. 2006, 12: 468-474. 10.3201/eid1203.051004.

    Article  Google Scholar 

  14. Hubalek Z, Halouzka J: West Nile fever--a reemerging mosquito-borne viral disease in Europe. Emerg Infect Dis. 1999, 5: 643-650. 10.3201/eid0505.990505.

    Article  CAS  Google Scholar 

  15. Deichmeister JM, Telang A: Abundance of West Nile virus mosquito vectors in relation to climate and landscape variables. J Vector Ecol. 2011, 36: 75-85. 10.1111/j.1948-7134.2011.00143.x.

    Article  Google Scholar 

  16. Epstein PR: West Nile virus and the climate. J Urban Health. 2001, 78: 367-371. 10.1093/jurban/78.2.367.

    Article  CAS  Google Scholar 

  17. Platonov AE, Shipulin GA, Shipulina OY, Tyutyunnik EN, Frolochkina TI, Lanciotti RS, et al: Outbreak of West Nile virus infection, Volgograd Region, Russia, 1999. Emerg Infect Dis. 2001, 7: 128-132. 10.3201/eid0701.010118.

    Article  CAS  Google Scholar 

  18. Danis K, Papa A, Papanikolaou E, Dougas G, Terzaki I, Baka A, et al. Ongoing outbreak of West Nile virus infection in humans, Greece, July to August 2011. Euro Surveill. 2011;16(34).

    Google Scholar 

  19. Danis K, Papa A, Theocharopoulos G, Dougas G, Athanasiou M, Detsis M, et al: Outbreak of West Nile virus infection in Greece, 2010. Emerg Infect Dis. 2011, 17: 1868-1872. 10.3201/eid1710.110525.

    Article  Google Scholar 

  20. Savard JL, Clergeau P, Mennechez G: Biodiversity concepts and urban ecosystems. Landsc Urban Plan. 2000, 48: 131-142. 10.1016/S0169-2046(00)00037-2.

    Article  Google Scholar 

  21. Hodgson JC, Spielman A, Komar N, Krahforst CF, Wallace GT, Pollack RJ: Interrupted blood-feeding by Culiseta melanura (Diptera: Culicidae) on European starlings. J Med Entomol. 2001, 38: 59-66. 10.1603/0022-2585-38.1.59.

    Article  CAS  Google Scholar 

  22. Bertolotti L, Kitron UD, Walker ED, Ruiz MO, Brawn JD, Loss SR, et al: Fine-scale genetic variation and evolution of West Nile Virus in a transmission “hot spot” in suburban Chicago, USA. Virology. 2008, 374: 381-389. 10.1016/j.virol.2007.12.040.

    Article  CAS  Google Scholar 

  23. LaDeau SL, Calder CA, Doran PJ, Marra PP: West Nile Virus impacts in american crow populations are associated with human land use and climate. Ecol Res. 2011, 26: 909-916. 10.1007/s11284-010-0725-z.

    Article  Google Scholar 

  24. Brown H, Duik-Wasser M, Andreadis T, Fish D: Remotely-sensed vegetation indices identify mosquito clusters of West Nile virus vectors in an urban landscape in the northeastern United States. Vector Borne Zoonotic Dis. 2008, 8: 197-206. 10.1089/vbz.2007.0154.

    Article  Google Scholar 

  25. Paz S, Malkinson D, Green MS, Tsioni G, Papa A, Danis K, et al: Permissive summer temperatures of the 2010 European West Nile fever upsurge. PLoS One. 2013, 8: e56398-10.1371/journal.pone.0056398.

    Article  CAS  Google Scholar 

  26. Paz S, Albersheim I: Influence of warming tendency on Culex pipiens population abundance and on the probability of West Nile fever outbreaks (Israeli Case Study: 2001–2005). Ecohealth. 2008, 5: 40-48. 10.1007/s10393-007-0150-0.

    Article  Google Scholar 

  27. Meyer RP, Hardy JL, Reisen WK: Diel changes in adult mosquito microhabitat temperatures and their relationship to the extrinsic incubation of arboviruses in mosquitoes in Kern County, California. J Med Entomol. 1990, 27: 607-614. 10.1093/jmedent/27.4.607.

    Article  CAS  Google Scholar 

  28. Ruiz MO, Chaves LF, Hamer GL, Sun T, Brown WM, Walker ED, et al: Local impact of temperature and precipitation on West Nile virus infection in Culex species mosquitoes in northeast Illinois, USA. Parasit Vectors. 2010, 3: 19-10.1186/1756-3305-3-19.

    Article  Google Scholar 

  29. 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: 309-317. 10.1603/0022-2585(2006)043[0309:EOTOTT]2.0.CO;2.

    Article  Google Scholar 

  30. Andrade CC, Maharaj PD, Reisen WK, Brault AC: North American West Nile virus genotype isolates demonstrate differential replicative capacities in response to temperature. J Gen Virol. 2011, 92: 2523-2533. 10.1099/vir.0.032318-0.

    Article  CAS  Google Scholar 

  31. EEA. Climate change, impacts and vulnerability in Europe. 2012, http://www.ecdc.europa.eu/en/healthtopics/west_nile_fever/west-nile-fever-maps/pages/index.aspx, . An indicator-based report. No 12/2012. 2012.

  32. Greer A, Ng V, Fisman D: Climate change and infectious diseases in North America: the road ahead. CMAJ. 2008, 178: 715-722.

    Google Scholar 

  33. Platonov AE, Fedorova MV, Karan LS, Shopenskaya TA, Platonova OV, Zhuravlev VI: Epidemiology of West Nile infection in Volgograd, Russia, in relation to climate change and mosquito (Diptera: Culicidae) bionomics. Parasitol Res. 2008, 103 (Suppl 1): S45-53. 10.1007/s00436-008-1050-0.

    Article  Google Scholar 

  34. Morin CW, Comrie AC: Modeled response of the West Nile virus vector Culex quinquefasciatus to changing climate using the dynamic mosquito simulation model. Int J Biometeorol. 2010, 54: 517-529. 10.1007/s00484-010-0349-6.

    Article  Google Scholar 

  35. Morin CW, Comrie AC: Regional and seasonal response of a West Nile virus vector to climate change. Proc Natl Acad Sci U S A. 2013, 110: 15620-15625. 10.1073/pnas.1307135110.

    Article  CAS  Google Scholar 

  36. Paz S. Climate change impacts on West Nile virus transmission in a global context. Philos Trans R Soc Lond B Biol Sci. 2015;370(1665). doi: 10.1098/rstb.2013.0561.

  37. IPCC. Special Report – Emission Scenarios. 2000, https://www.ipcc.ch/pdf/special-reports/spm/sres-en.pdf, . Summary for policymakers:

  38. Semenza JC, Domanović D: Blood supply under threat. Nature Climate Change. 2013, 3: 432-435. 10.1038/nclimate1867.

    Article  Google Scholar 

  39. 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: 1298-1306. 10.3201/eid1808.111558.

    Article  Google Scholar 

  40. Semenza JC, Zeller H. Integrated surveillance for prevention and control of emerging vector-borne diseases in Europe. Euro Surveill. 2014;19(13).

  41. NCAR Climate change scenarios:. https://gisclimatechange.ucar.edu/

  42. Pachauri RKaR, A. (Eds.). Contribution of Working Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. 2007, http://www.ipcc.ch/publications_and_data/publications_ipcc_fourth_assessment_report_synthesis_report.htm, .

  43. European Center for Disese Prevention and Control (ECDC). West Nile fever maps. 2014, http://www.ecdc.europa.eu/en/healthtopics/west_nile_fever/west-nile-fever-maps/pages/index.aspx, .

  44. EUROSTAT. Nomenclature of territorial units for statistics. 2014, http://www.ecdc.europa.eu/en/healthtopics/west_nile_fever/west-nile-fever-maps/pages/index.aspx, .

  45. NASA. Socioeconomic Data and Application Centre (SEDAC): Gridded population of the world (gpw), v3. 2014, http://www.ecdc.europa.eu/en/healthtopics/west_nile_fever/west-nile-fever-maps/pages/index.aspx, .

  46. 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: e46882-10.1371/journal.pone.0046882.

    Article  CAS  Google Scholar 

  47. 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: 20-31. 10.1016/j.prevetmed.2006.11.008.

    Article  CAS  Google Scholar 

  48. Ward MP: Equine West Nile virus disease occurrence and the Normalized Difference Vegetation Index. Prev Vet Med. 2009, 88: 205-212. 10.1016/j.prevetmed.2008.10.003.

    Article  Google Scholar 

  49. ECDC European Up-Front Risk Assessment Tool (EUFRAT). http://eufrattool.ecdc.europa.eu/

  50. Semenza JC, Sudre B, Miniota J, Rossi M, Hu W, Kossowsky D, et al: International dispersal of dengue through air travel: importation risk for Europe. PLoS Negl Trop Dis. 2014, 8: e3278-10.1371/journal.pntd.0003278.

    Article  Google Scholar 

  51. Dodd RY, Leiby DA: Emerging infectious threats to the blood supply. Annu Rev Med. 2004, 55: 191-207. 10.1146/annurev.med.55.091902.103831.

    Article  CAS  Google Scholar 

  52. Alter HJ, Stramer SL, Dodd RY: Emerging infectious diseases that threaten the blood supply. Semin Hematol. 2007, 44: 32-41. 10.1053/j.seminhematol.2006.09.016.

    Article  Google Scholar 

  53. Stramer SL: Current perspectives in transfusion-transmitted infectious diseases: emerging and re-emerging infections. ISBT Sci Ser. 2014, 9: 30-36. 10.1111/voxs.12070.

    Article  Google Scholar 

  54. Stramer SL, Hollinger FB, Katz LM, Kleinman S, Metzel PS, Gregory KR, et al: Emerging infectious disease agents and their potential threat to transfusion safety. Transfusion. 2009, 49 (Suppl 2): 1S-29S. 10.1111/j.1537-2995.2009.02279.x.

    Article  Google Scholar 

  55. Orton SL, Stramer SL, Dodd RY: Self-reported symptoms associated with West Nile virus infection in RNA-positive blood donors. Transfusion. 2006, 46: 272-277. 10.1111/j.1537-2995.2006.00710.x.

    Article  CAS  Google Scholar 

  56. Lieshout-Krikke RW, Zaaijer HL, Prinsze FJ: The yield of temporary exclusion of blood donors, exposed to emerging infections abroad. Vox Sang. 2013, 104: 12-18. 10.1111/j.1423-0410.2012.01631.x.

    Article  CAS  Google Scholar 

  57. Busch MP, Caglioti S, Robertson EF, McAuley JD, Tobler LH, Kamel H, et al: Screening the blood supply for West Nile virus RNA by nucleic acid amplification testing. N Engl J Med. 2005, 353: 460-467. 10.1056/NEJMoa044029.

    Article  CAS  Google Scholar 

  58. Kleinman SH, Williams JD, Robertson G, Caglioti S, Williams RC, Spizman R, et al: West Nile virus testing experience in 2007: evaluation of different criteria for triggering individual-donation nucleic acid testing. Transfusion. 2009, 49: 1160-1170. 10.1111/j.1537-2995.2009.02127.x.

    Article  Google Scholar 

  59. Custer B, Tomasulo PA, Murphy EL, Caglioti S, Harpool D, McEvoy P, et al: Triggers for switching from minipool testing by nucleic acid technology to individual-donation nucleic acid testing for West Nile virus: analysis of 2003 data to inform 2004 decision making. Transfusion. 2004, 44: 1547-1554. 10.1111/j.0041-1132.2004.04227.x.

    Article  Google Scholar 

  60. Gallian P, Vignoli C, Dombey AM, Mayaudon V, Lin L, Galichet V, et al: Inactivation of a European strain of West Nile virus in single- donor platelet concentrate using the INTERCEPT blood system. Vox Sang. 2006, 91: 345-347. 10.1111/j.1423-0410.2006.00844.x.

    Article  CAS  Google Scholar 

  61. Vanlandingham DL, Keil SD, Horne KM, Pyles R, Goodrich RP, Higgs S: Photochemical inactivation of chikungunya virus in plasma and platelets using the Mirasol pathogen reduction technology system. Transfusion. 2013, 53: 284-290. 10.1111/j.1537-2995.2012.03717.x.

    Article  CAS  Google Scholar 

  62. Burnouf T, Chou ML, Cheng LH, Li ZR, Wu YW, El-Ekiaby M, et al: Dengue virus inactivation by minipool TnBP/Triton X-45 treatment of plasma and cryoprecipitate. Vox Sang. 2013, 104: 1-6. 10.1111/j.1423-0410.2012.01621.x.

    Article  CAS  Google Scholar 

  63. Bambrick HJ, Woodruff RE, Hanigan IC. Climate change could threaten blood supply by altering the distribution of vector-borne disease: an Australian case-study. Glob Health Action. 2009;2. doi: 10.3402/gha.v2i0.2059.

  64. Xinhu L, Jinchao S, Tao L, Dixon J, Guoqin Z, Hong Y. Urbanization and health in China, thinking at the national, local and individual levels. Environ Health. 2016;15(Suppl 1):xx.

    Google Scholar 

Download references

Acknowledgements

The authors would like to thank Mikel Santamaria Guerro for his assistance with the maps and Drs. Yvan Hutin and Mike Catchpole at ECDC for critical feedback on the manuscript. We would also like to thank the editors and anonymous reviewers for their constructive feedback on our manuscript.

Declarations

The study was funded by ECDC including the publication charges and several co-authors (JCS, LE, BS, DD) are employed by ECDC. The authors declare no financial competing interests. The opinions in this paper reflect that of the authors and not necessarily that of ECDC.

This article has been published as part of Environmental Health Volume 15 Suppl 1, 2016: Healthy-Polis: Challenges and Opportunities for Urban Environmental Health and Sustainability. The full contents of the supplement can be found at http://www.ehjournal.net/supplements/15/S1.

Peer review

Peer review reports for this article are attached as Additional file 1.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Jan C. Semenza.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JCS conceived, designed the study and wrote the paper. AT developed and ran the multivariate model. LE, DD and JCS conducted the assessment of the safety of the blood supply. BS contributed to the epidemiological data collection and model development. SP computed the temperature abnormalities and climate change projections. All authors reviewed the draft manuscript and contributed to the writing and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Semenza, J.C., Tran, A., Espinosa, L. et al. Climate change projections of West Nile virus infections in Europe: implications for blood safety practices. Environ Health 15 (Suppl 1), S28 (2016). https://0-doi-org.brum.beds.ac.uk/10.1186/s12940-016-0105-4

Download citation

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12940-016-0105-4

Keywords