Email updates

Keep up to date with the latest news and content from BMC Public Health and BioMed Central.

Open Access Highly Accessed Research article

Climate change and the emergence of vector-borne diseases in Europe: case study of dengue fever

Maha Bouzid1*, Felipe J Colón-González23, Tobias Lung45, Iain R Lake2 and Paul R Hunter1

Author Affiliations

1 Norwich Medical School, University of East Anglia, Norwich, UK

2 School of Environmental Sciences, University of East Anglia, Norwich, UK

3 The Abdus Salam International Centre for Theoretical Physics, Earth System Physics Section, Trieste, Italy

4 Joint Research Centre, European Commission, Institute for Environment and Sustainability, Ispra, Italy

5 Present address: European Environment Agency, Copenhagen, Denmark

For all author emails, please log on.

BMC Public Health 2014, 14:781  doi:10.1186/1471-2458-14-781

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2458/14/781


Received:2 April 2014
Accepted:24 July 2014
Published:22 August 2014

© 2014 Bouzid et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.

Abstract

Background

Dengue fever is the most prevalent mosquito-borne viral disease worldwide. Dengue transmission is critically dependent on climatic factors and there is much concern as to whether climate change would spread the disease to areas currently unaffected. The occurrence of autochthonous infections in Croatia and France in 2010 has raised concerns about a potential re-emergence of dengue in Europe. The objective of this study is to estimate dengue risk in Europe under climate change scenarios.

Methods

We used a Generalized Additive Model (GAM) to estimate dengue fever risk as a function of climatic variables (maximum temperature, minimum temperature, precipitation, humidity) and socioeconomic factors (population density, urbanisation, GDP per capita and population size), under contemporary conditions (1985–2007) in Mexico. We then used our model estimates to project dengue incidence under baseline conditions (1961–1990) and three climate change scenarios: short-term 2011–2040, medium-term 2041–2070 and long-term 2071–2100 across Europe. The model was used to calculate average number of yearly dengue cases at a spatial resolution of 10 × 10 km grid covering all land surface of the currently 27 EU member states. To our knowledge, this is the first attempt to model dengue fever risk in Europe in terms of disease occurrence rather than mosquito presence.

Results

The results were presented using Geographical Information System (GIS) and allowed identification of areas at high risk. Dengue fever hot spots were clustered around the coastal areas of the Mediterranean and Adriatic seas and the Po Valley in northern Italy.

Conclusions

This risk assessment study is likely to be a valuable tool assisting effective and targeted adaptation responses to reduce the likely increased burden of dengue fever in a warmer world.

Keywords:
Dengue fever; Climate change; Europe; Disease risk

Background

Several vector-borne diseases are spread in Europe and the effect of climate change on disease distribution has been extensively discussed [1-5]. Most authors consider that climate change is likely to have greatest impact on dengue fever, West Nile fever, chikungunya fever, malaria, leishmaniasis, tick-borne encephalitis, Lyme borreliosis, Crimean-Congo haemorrhagic fever, spotted fever rickettsioses, Yellow fever and Rift Valley fever. One disease that has received much interest in recent years is dengue fever. Dengue is a mosquito-borne disease caused by an RNA virus of the genus Flavivirus. Uncomplicated dengue can present with fever, headache and muscle and joint pains. A proportion of infections can develop into severe forms namely dengue haemorrhagic fever and dengue shock syndrome, which are associated with higher mortality rates. Dengue fever is endemic in over 100 countries in Africa, the Americas, the Eastern Mediterranean, South-east Asia and the Western Pacific, with the last two regions being the most seriously affected [6]. It is estimated that over 50 million new dengue fever infections and approximately 12,000 deaths, mainly among children, occur worldwide every year [7].

There has been a significant global increase in dengue incidence and it is currently considered the most important human arboviral disease worldwide. The successful spread of dengue has been attributed to various factors including population growth, urbanization, global travel, and environmental conditions. In Europe, dengue fever is rare but cases are imported every year by tourists returning from endemic areas. Recently, autochthonous dengue cases have been reported in Croatia and France, highlighting the suitability of these regions for dengue transmission [8,9]. These cases have raised concerns about the potential for the emergence of dengue fever in Europe especially with predicted climate change.

One of the reasons for these concerns is that dengue vectors are already present within Europe. Aedes (Stegomyia) aegypti (Linneaus) is the major urban vector of dengue worldwide [10]. A. aegypti is closely associated with humans and human habitations. Female mosquitoes lay their eggs on or near water surface in natural or artificial containers [10,11]. Aedes albopictus is the secondary vector of dengue fever and is adapted to the peridomestic environment [12]. According to the “European Network for arthropod vector surveillance for human public health” (VBORNET) (http://www.vbornet.eu/ webcite, last accessed June 2014), A. albopictus is present in many European countries: Spain, France, Switzerland, Italy, Slovenia, Croatia, Bosnia and Herzegovina, Serbia, Montenegro, Albania, Greece, Monaco, San Marino, Bulgaria, the Netherlands and Russia. By contrast, A. aegypti has only been reported from Madeira, the Netherlands, Georgia and southern Russia.

There is much debate about how future climate change will affect dengue risk, especially in countries where the disease is not currently endemic [10,13,14]. Recent studies have modelled the future dengue distribution under predicted climate change either on a global scale [15-19] or in endemic countries [20,21]. These models have suggested a latitudinal and altitudinal expansion of the geographical range of dengue.

In Europe, there have been too few dengue cases to conduct a rigorous analysis. Consequently, estimation of dengue risk has so far relied on past, current and projected future distribution of A. albopictus[22]. Although presence of the vector is necessary for dengue to become endemic, vector presence is not sufficient in itself to determine disease occurrence [23]. The objective of this study is to model dengue risk based on clinical data. We have used one of the largest and more spatially diverse dengue dataset yet assembled to compute significant relationships between dengue and weather parameters [24]. Subsequently, the model outputs were used to project dengue risk across Europe under climate change scenarios.

Methods

Mexican data

The dengue dataset was primarily developed for a study of the effects of weather on dengue incidence across Mexico [24]. Dengue data comprised state-specific monthly reports of laboratory confirmed dengue cases, retrieved from the Mexican Health Secretariat (http://www.epidemiologia.salud.gob.mx/anuario/html/anuarios.html webcite, last accessed June 2014) for the period January 1985 to December 2007. Monthly average minimum and maximum temperatures and monthly precipitation for each state were provided by the Mexican National Meteorological Service. Monthly mean humidity was retrieved from the National Centers for Environmental Prediction and National Center for Atmospheric Research (NCEP/NCAR) “Reanalysis 1” (http://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.pressure.html webcite, last accessed June 2014). Yearly Gross domestic product (GDP) per capita (PPP in constant 2005 international dollars) was obtained from the World Bank at the national level (http://data.worldbank.org/country/mexico webcite, last accessed June 2014). State-specific GDP estimates were computed as previously described [24]. The proportion of people living in urban areas was retrieved from the Mexican Chamber of Deputies (http://www.cefp.gob.mx/intr/bancosdeinformacion/estatales/indicadores_socioeconomicos/is003.xls webcite, last accessed June 2014). Population density was calculated by normalising population to state area size. Table 1 presents the summary statistics for these variables.

Table 1. Summary of statistical characteristics of the climatic and socioeconomic variables used for this study in Mexico and Europe

Model calibration

Generalized Additive Models (GAMs) are semi-parametric extensions of the generalized linear model (GLM), where the linear predictor ΣβjXj is replaced by a sum (hence the name additive) of smooth functions of covariates Σsj(Xj) [25]. Like in GLMs, GAMs allow the exploration of nonlinear data structures in the context of exponential family distributions (e.g. Poisson and Binomial), and use link functions to establish relationships between the mean of the outcome variable and the predictors [26,27]. Unlike GLMs, GAMs automatically identify and estimate the optimal degree of nonlinearity of the model directly from the data [28]. In our study, the expected number of dengue cases E(yti) ≡ μit for State i at time t was assumed to follow an overdispersed Poisson distribution described by:

<a onClick="popup('http://www.biomedcentral.com/1471-2458/14/781/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2458/14/781/mathml/M1">View MathML</a>

where g(.) is a log link function of the expectation μit ≡ E(yti) with yti denoting the time series of dengue counts. The logarithm of the population (ξ) at time t and state i is included as an exposure variable to standardise the dengue data by population. Weather has a delayed effect on dengue incidence. Therefore, we specified our j-th meteorological variables Xjti within biologically and physically plausible time lags based on literature reports in Mexico [29-31]. Weather variables comprised average monthly minimum (Tmin) and maximum (Tmax) temperatures, monthly precipitation (Precip) and average monthly relative humidity (Humid). All weather parameters were lagged 1 and 2 months (Tmin1:2, Tmax1:2, Precip1:2, Humid1:2). The term sj(.) corresponds to univariate smooth functions defined by penalized cubic regression splines. We adjusted our model for the effects of socioeconomic variables Zkit represented by GDP per capita, the proportion of the population living in urban settlements and population density. Socioeconomic variables entered the model linearly. Analyses were conducted in R version 2.15.0 [32].

Many epidemiological datasets are likely to be dominated by long-term and seasonal trends. Therefore, adjusting the regression models for these patterns is necessary to separate them from the effects of weather parameters on the health variable [33]. Our model does not account for seasonal trends as seasonality for Europe is unlikely to be similar to Mexico given the wider range of temperatures between summer and winter. Although mosquito presence is a key factor in the epidemiology and occurrence of the disease, to our knowledge there are no state-specific long-term time series of mosquito presence across Mexico. Consequently, data on mosquito presence could not be incorporated into our model. The GAM-estimated relationships between dengue, weather and socioeconomic development in Mexico were then used to project dengue fever risk across Europe.

European data and dengue fever risk modelling

European climate data were retrieved from the regional climate model COSMO-CLM (CCLM), forced with output from the coupled atmosphere–ocean global climate model (GCM) ECHAM5/MPIOM [34]. These regional simulations represent aerosol and GHG forcing according to the A1B scenario of the Special Report on Emissions Scenarios (SRES) of the IPCC [35]. A1B corresponds to a projected increase in global surface temperature of 2.8°C in 2090–2099 (relative to 1980–1999) and a likely range of up to 4.4°C [36]. It assumes rapid economic growth, rapid introduction of efficient technologies, convergence among regions and a balance across energy sources. The regional climate data correspond to the period 1961–2100, with a domain covering the entire European continent at a resolution of about 18 × 18 km. Data were re-scaled to a grid cell size of 10 × 10 km for the purpose of this study. The same four monthly climatic variables (Tmin, Tmax, Precip, Humid) lagged 1 and 2 months, as used for model calibration with the Mexican data, were calculated over four time periods, (a) baseline 1961–1990, (b) short-term scenario 2011–2040, (c) medium-term scenario 2041–2070, and (d) long-term scenario 2071–2100.

GDP per capita data were retrieved from EUROSTAT (in Euros) and converted into constant 2005 international dollars to be concordant with the Mexican data used for model calibration. Country level data from the World Development Indicators dataset (http://databank.worldbank.org/data/Databases.aspx webcite, last accessed June 2014) were disaggregated to NUTS-3 level (Nomenclature of territorial units for statistics) by using the NUTS-3 level shares for each country as calculated from EUROSTAT. Then, an areal weighting approach was employed to convert the NUTS-3 data into the 10 × 10 km grid (see Additional file 1). Areal weighting is commonly used to transform administrative boundary data to raster format, whereby each grid cell is assigned a value according to the percentage of its area covered by the overlying administrative region [37].

Additional file 1. Maps of Gross domestic product (GDP) of 2006 at NUTS-3 level calculated from Word Bank and EUROSTAT data (Left), and degree of urbanisation as derived from the EFGS GEOSTAT 2006 population grid dataset (right).

Format: PNG Size: 2.4MB Download fileOpen Data

The proportion of population living in urban areas and total population data were retrieved from the GEOSTAT 2006 population grid dataset of the European Forum for Geostatistics (EFGS) (http://www.efgs.info/ webcite, last accessed June 2014) at a spatial resolution of 1 × 1 km. Urban clusters were defined by two criteria first each grid cell of 1 × 1 km must have a minimum population density of 300 people per km2 and second clusters of adjoining grid cells must accommodate at least 5000 people, in line with the definitions used by the European Commission [38]. The total number of urban population for each 10 × 10 km grid cell was extracted and divided by total population to obtain proportion of population in urban area (see Additional file 1). Due to the lack of projections both in terms of SRES scenario and spatial detail, the socioeconomic variables were held constant at their mean value for baseline conditions in order to isolate the effects of climate.

Mapping of dengue fever risk

The model was used to project monthly dengue cases, which were aggregated to calculate the average number of cases per year for each time period. These were used to generate dengue risk maps using ArcGIS 10.1 (http://www.esri.com/software/arcgis webcite, last accessed June 2014). In total, four maps were produced corresponding to the time periods of study. Identical class sizes were applied across all four time periods in order to ensure that value changes could be observed over time. The map series employ a bipolar hue progression [39] ranging from green (no/low risk areas) to bright red (areas with the highest dengue risk). Moreover, a second map series was generated, that normalises dengue cases by total population to derive dengue incidence. We used colours ranging from blue for no/low risk to cherry brown for high risk areas. In addition, standard error for each grid cell was calculated. Standard error values were subjected to the same aggregation and averaging procedure as for dengue number of cases and dengue incidence and were used to produce maps of uncertainty.

Results

The dengue-weather relationships estimated by our Poisson GAM based on the Mexican data are presented in Table 2. The model explained 44% of the deviance within the Mexican dataset. Figure 1 shows the estimated effects of weather variables on number of dengue cases. All climate parameters were statistically significant in a highly non-linear way. The greatest effect was associated with monthly average of minimum temperature followed by monthly relative humidity (both variables lagged one and two months). For socioeconomic variables, population density, degree of urbanization and log population were all significantly associated with dengue incidence (Table 2).The GAM estimated relationships were used to project dengue fever risk in Europe under climate change conditions expressed as dengue cases (cases/year/10×10 km grid). For the baseline period (1960–1990), number of dengue cases are between 0 and 0.6 for most European areas, corresponding to an incidence of less than 2 per 100 000 inhabitants (Figures 2 and 3). Over time, an increase in dengue risk is projected, with highest incidence rates found for the long-term scenario 2070–2100. Indeed, for the baseline period hardly any grid cell had incidence rate exceeding 10 per 100 000 inhabitants, while a substantial amount of grid cells are within this category when considering the long-term-scenario, mostly localised in southern Europe. For each estimated dengue incidence, standard error was calculated and values presented as maps of uncertainty (Figure 4). The general trend was that standard errors tend to correlate with incidence rates, as would be expected from a Poisson model. The maps also highlight that the standard errors are not consistent across the continent.

Table 2. Model estimates of the effects of weather and socioeconomic variables on dengue

thumbnailFigure 1. GAM-estimated relationships between average monthly dengue cases and average monthly Tmin (A), Tmax (B), Humidity (C), and precipitation (D), all lagged 1 and 2 months. The x axis represents increasing variations in the meteorological covariates. The y axis indicates the contribution of the smoother to the fitted values. The y axis is labelled s(cov, edf), where cov indicates the name of the covariate, and edf represents the estimated degrees of freedom of the smooth function used to represent its relationship with number of dengue cases. The red lines indicate the maximum likelihood estimates, and the grey shaded areas represent the 95% confidence intervals. The rug at the bottom of the figures indicate observed values of the covariates.

thumbnailFigure 2. Average expected number of dengue cases in Europe modelled using GAM model for baseline conditions and climate change scenarios for early, medium and late century. Number of cases was calculated for each 10 × 10 km grid.

It can be seen from Figures 2 and 3 that the risk is not equally distributed across Europe. Generally, southern Europe appears at higher risk, with most of the coastal areas being particularly affected. In contrast, northern Europe, the British Isles and most of central Europe show virtually no risk in the baseline period. Over time, the risk in southern Europe and in particular along the Mediterranean coast is projected to increase considerably, with highest incidence revealed along the Italian coast, the Po Valley region, the Spanish Mediterranean coast and southern Spain in general (Figures 2 and 3). In other parts of Europe a change from virtually no risk for the baseline period to incidence rate of up to 10 cases per 100 000 inhabitants are projected, such as in large parts of France, south-western Germany, Hungary, and the Balkan region. By contrast, in northern Europe, for most of the British Isles and the Baltic states, the risk is projected to remain virtually zero even for the long-term scenario 2070–2100.

thumbnailFigure 3. Dengue fever incidence rate expressed as number of cases per 100 000 inhabitants per year for baseline conditions and climate change scenarios for early, medium and late century.

thumbnailFigure 4. Maps of uncertainty showing standard errors for projected average number of dengue cases for baseline conditions and climate change scenarios for early, medium and late century.

Discussion

We have presented the first ever projections of future dengue fever risk in Europe under climate change based on empirical modelling of laboratory confirmed dengue incidence as a function of climate and socioeconomic variables. Our study has shown that the risk of dengue fever is likely to increase in Europe under climate change, but that almost all of the excess risk will fall on the coastal areas of the Mediterranean and Adriatic seas and the North Eastern part of Italy, particularly the Po Valley. Although we have only modelled dengue fever, our findings may have implications for other mosquito-borne diseases such as Chikungunya, which share the same vector species and may have similar transmission patterns to dengue fever.

Previous work in the area has primarily focused on the expected future distribution of the Aedes vector [40-44]. Whilst such studies have proven helpful to determine the potential presence of the vector in a given area, the presence of the vector does not necessary translate into disease occurrence. Because our model is based on disease occurrence, the GAM-estimated relationships are likely to be more useful for estimating dengue risk across Europe than models based only on the mere presence of the mosquito. Nevertheless, our results are generally in agreement with the conclusions of projections based on vector distribution. One such example of concordance is Italy, which was identified as a potential dengue hot spot using both approaches (i.e. vector presence and dengue risk). The most noticeable discrepancies for dengue fever risk, however, were associated with Spain and France. In our study, southern and eastern Spain is associated with increased dengue risk, however, A. albopictus maps show that this area is likely to be unsuitable by mid/end century [40-43], probably because of hotter and drier weather conditions. On the contrary, France is considered here at medium risk (excluding Mediterranean areas), while it is considered an area of high suitability for A. albopictus[40,41,43]. Possible explanations for these discrepancies include different climate change scenarios and dissimilar climatic variables incorporated in modelling approaches. A recent study by Rogers and colleagues established a dengue risk map for Europe based on a global risk map taking into account vector presence, disease occurrence and various environmental factors [45]. They showed that while most Europe is at low risk, most major cities combining warmer temperatures and higher population density are highly suitable for dengue transmission. However, as the authors acknowledged, several other factors can influence disease occurrence and transmission. In our study, most cities are not at an increased risk until mid to end of the century. Understandably, it is difficult to compare two risk maps generated using different methodologies and data sets. Nevertheless, European dengue hot spots identified in this and other studies should be made aware of the projected risk.

Herein, we identified a large European geographical area permissive for dengue fever transmission. Whether dengue will become endemic in a particular area depends on various climatic and non-climatic factors, in addition to disease risk in neighbouring areas, which makes estimation of actual incidence problematic. Nevertheless, dengue is unlikely to become endemic in areas of moderate risk, especially if nearby areas have low risk. Consequently we would hypothesise that, should dengue fever become endemic in Europe, it is likely to be primarily in the Mediterranean and Adriatic coastal plains and the Po Valley area of Italy. This does not mean that localised outbreaks of dengue would not occur elsewhere, but that if they did, they would be less likely to be self-sustaining.

For the purpose of this study, we have used Mexican dengue fever surveillance data to project dengue incidence in Europe. This could be considered a limitation because of applicability to European settings and transferability issues. It would clearly have been preferable to use European dengue data but at present there have been too few cases in Europe for any meaningful analysis. In the absence of worthwhile European data, the Mexican dengue surveillance dataset is without doubt the best alternate source of empirical data. The Mexican dataset comprises the largest such set yet assembled with monthly data and sub-national resolution. Furthermore, Mexico is a large country comprising multiple climate zones and thus providing the opportunity of modelling climate impacts on dengue over a wide range of climatic conditions. Mexico is a middle income country, whose socioeconomic status could be comparable to some European countries http://data.worldbank.org/data-catalog/world-development-indicators webcite (last accessed June 2014). Nevertheless, there are several mismatches. One particular issue is seasonal climatic variation, as seasons in Mexico may not match those in Europe. For example, Mexico is highly affected by the El Niño–Southern Oscillation [46], which does not have such an impact in Europe. Similarly, winter climatic conditions can be very different between Mexico and Europe, therefore influencing overwintering and survival of Aedes populations. In order to test for this variability, we fitted the model with and without the seasonal term and found that the shape of the estimated relationships graph was similar with a slight difference in order of magnitude. Additional file 2 shows the range of the two weather variables most significantly associated with dengue risk (Tmin and humidity). There is a major area of overlap between the two sufficient to give validity to the use of the Mexico data. The major European areas that do not overlap with Mexico are unlikely to be in areas at risk from dengue. Another important point to consider is the different socioeconomic conditions and cultural habits between Mexico and Europe. Whilst we were not able to model many of these differences due to the absence of adequate data, some of the variables may have little impact on overall dengue risk. For instance, GDP was not found to be significantly associated with dengue fever risk. However, this could be due to data constrains because GDP were available as yearly data that arose from linear interpolations based on 5 years interval. Taking into account these constraints, we have been careful not to be too specific about how many extra cases of dengue we are likely to see in Europe, rather on identifying high risk areas. Clearly even in highly conducive areas, dengue fever will not become endemic if it is not introduced at some point and so the development and spread of dengue endemicity is likely to be a stochastic process, nevertheless it is relying on vector presence, virus introduction and host susceptibility.

Additional file 2. Range of Tmin and humidity, the two weather variables most significantly associated with dengue risk, in Mexico and Europe.

Format: PNG Size: 77KB Download fileOpen Data

An important issue is that the model is produced for a country where dengue is endemic. If dengue was introduced into Europe, then it could spread rapidly in the early years of its establishment and become endemic. This is because almost all Europeans would be immunologically naïve and therefore actual cases could outstrip our projections. A further source of uncertainty would be the adaptation of the European health authorities to the emergence of dengue, including health practitioners’ awareness and effective diagnostic and treatment measures. However, even well-staffed health services with adequate infrastructure could struggle to manage dengue fever [23]. Further response measures that could lower transmission rate, while awaiting the development of an effective vaccine, should include integrated vector management. However, the effectiveness of vector control strategies is not always supported by adequate evidence based evaluations [47]. In addition, unless vector control is performed in a sustainable manner, it is most likely to be inefficient.

Another limitation of the study is related to the mosquito vector. Although, A. albopictus is present in Mexico, dengue is mainly transmitted by A. aegypti. This primary dengue vector is responsible for major dengue epidemics and the severe life- threatening form of the disease [12]. In Europe, A. aegypti is only present in Madeira, where it caused sporadic cases and a sustained dengue outbreak [48]. The main Aedes species in Europe is A. albopictus. This species is associated with sporadic dengue cases due to its limited competence related to its feeding behaviour and its relative recent adaptation to flaviviruses (including dengue fever virus) [49]. Consequently, our results are likely to over-estimate dengue risk in Europe. Additionally, some dengue virus serotypes were shown to cause more severe symptoms and spread more easily [50], therefore the impact of dengue introduction in Europe and subsequent transmission is influenced by virus-vector interaction and the associated risk and severity could either increase or decrease accordingly. Nevertheless, living standards in Europe are likely to limit dengue spread as has been reported in Texas [51], consequently, the actual dengue incidence could be much lower than projected using our model.

One finding that may cause some concern is that under baseline conditions some areas are identified as being at increased risk of dengue fever, when dengue in Europe is effectively non-existent. This is legitimate because the model is projecting areas where, given the provided meteorological and socioeconomic conditions, dengue fever may occur, independently of other confounding factors including vector presence, virus circulation and control measures. Undoubtedly, some European areas are permissive for dengue fever transmission as supported by endogenous cases in France and Croatia in 2010 [8,9]. The French cases were recorded in Nice, which was indeed highlighted as high risk area in our model. Croatian cases were on the Peljesac peninsula and the island of Korcula (outside the EU27 and therefore not modelled here). Furthermore, in 2007, there was an outbreak of Chikungunya (another viral disease spread by A. albopictus mosquitoes) that affected north eastern Italy on the Adriatric coast [52,53]. Although Chikungunya virus shows adaptive mutation for A. albopictus[49], not observed for dengue virus, this area is clearly permissive for the mosquito vector and is associated with the highest projected dengue fever risk in our model. These sporadic cases and outbreaks confirm the general spatial pattern of dengue risk as estimated by our model.

An issue valid for both Mexico and Europe is that Mexican dengue fever data is based on laboratory confirmed reported cases of infection. While on one hand, and in particular for Europe, our approach of modelling dengue risk based on reported cases is considered novel and unique, on the other hand such reports are known to substantially underestimate the actual number of cases because a significant proportion of infections are not diagnosed and reported. A limitation of dengue surveillance in Mexico is its reliance on a passive surveillance system based on unspecific symptoms, coupled with low awareness of health practitioners and limited access to reliable diagnostic tests. Estimation of the sensitivity of dengue surveillance systems varies inter and intra countries but it has been reported that for every recorded case, there may be somewhere between 10–27 cases that go unreported [54]. This could mean that the number of dengue fever cases could be substantially higher than estimated here. It is not known what the sensitivity of dengue surveillance systems in Europe would be, however, asymptomatic and mild cases could go undiagnosed.

In order to assess the impact of climate change on dengue risk in Europe, we used predictions based on the A1B scenario because it is considered more realistic in light of the current global emissions. This is particularly relevant when compared to other more extreme scenarios (such as the low emission SRES B1, or the newer RCP2.6 scenarios) that assume drastic CO2 reductions globally in the coming decades, which is unlikely to happen. Running the current model using additional climate change scenarios would add value to the predictions and allow comparative analysis and could be done as a subsequent study. In order to assess the effect of climate change on dengue risk, all non-climatic variables were assumed to remain at their baseline levels, while this could be considered a limitation, some socioeconomic variables were not significantly associated with dengue risk. In addition, population and urbanisation projections for Europe show that minor changes are expected (with some local variation), especially when compared to other parts of the world (where significant increase in population size and urbanisation are expected until 2100) (http://esa.un.org/wpp/ webcite and http://esa.un.org/unpd/wup/ webcite last accessed August 2014).

A recent systematic review of quantitative models assessing the impact of climate change on dengue transmission by Naish and colleagues [55] found that despite using different methodologies, most models consider that temperature is the most important climatic factor driving dengue transmission but that precipitation and humidity are also important, which is in accordance with our model. Despite some methodological issues, most models report increased climatic suitability and expansion of geographical range under various climate change scenarios and in different regions of the world [55]. Improved climate change scenarios and better understanding of vector-borne diseases biology and transmission are likely to contribute to more accurate disease risk models in the future.

Conclusions

This study allowed modelling of dengue fever risk in Europe based on actual clinical data. The model calibrated under Mexican conditions resulted in reliable and geographically meaningful patterns of projected dengue fever risk in Europe. The risk maps indicate that climate change is likely to contribute to increased dengue risk (and possibly other mosquito-borne diseases) in many parts of Europe, especially towards the end of the century. The areas of greatest increased risk are projected to be clustered around the Mediterranean and Adriatic coasts and in northern Italy. The exact incidence is dependent on several other factors, some of which we were unable to model at this stage (such as vaccine development). Nevertheless public health agencies in high risk areas need to plan, implement and evaluate active entomological reporting and sentinel clinical surveillance and should aim to improve awareness of the increased risk amongst health practitioners and the general public.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MB and PRH designed the study and performed data analysis and interpretation. FJCC provided Mexican data and undertook the modelling. TL provided European data and performed geographical mapping of the model results. IRL was involved with modelling and results interpretation. All authors read and approved the final manuscript.

Acknowledgments

This work was partially supported by the European Union under the RESPONSES project (Grant Agreement number 244092). Professor Hunter and Dr. Lake are part funded by The National Institute for Health Research Health Protection Research Unit in Emergency Preparedness and Response at King’s College London. FJCG was supported by the Mexican National Council of Science and Technology (CONACYT).

References

  1. DH: Health effects of climate change in the UK 2008: an update of the department of health report 2001/2002. http://webarchive.nationalarchives.gov.uk/20130107105354/http://www.dh.gov.uk/en/Publicationsandstatistics/Publications/PublicationsPolicyAndGuidance/DH_080702 webcite

  2. Hunter PR: Climate change and waterborne and vector-borne disease.

    J Appl Microbiol 2003, 94(Suppl):37S-46S. PubMed Abstract | Publisher Full Text OpenURL

  3. McMichael AJ, Woodruff RE, Hales S: Climate change and human health: present and future risks.

    Lancet 2006, 367:859-869. PubMed Abstract | Publisher Full Text OpenURL

  4. Semenza JC, Menne B: Climate change and infectious diseases in Europe.

    Lancet Infect Dis 2009, 9:365-375. PubMed Abstract | Publisher Full Text OpenURL

  5. WHO: Protecting health from climate change: global research priorities. http://www.who.int/globalchange/publications/9789241598187/en/ webcite

  6. WHO: Dengue media centre: dengue and haemorrhagic fever. http://www.who.int/mediacentre/factsheets/fs117/en/ webcite 2009

  7. Lawrence J: DengueNet--WHO’s internet-based system for the global surveillance of dengue fever and dengue haemorrhagic fever (dengue/DHF).

    Euro Surveill Wkly Release 2002, 6:1883. OpenURL

  8. Gjenero-Margan I, Aleraj B, Krajcar D, Lesnikar V, Klobucar A, Pem-Novosel I, Kurecic-Filipovic S, Komparak S, Martic R, Duricic S, Betica-Radic L, Okmadzic J, Vilibic-Cavlek T, Babic-Erceg A, Turkovic B, Avsic-Zupanc T, Radic I, Ljubic M, Sarac K, Benic N, Mlinaric-Galinovic G: Autochthonous dengue fever in Croatia, August-September 2010.

    Euro Surveill 2011, 16:pii=19805. OpenURL

  9. La Ruche G, Souares Y, Armengaud A, Peloux-Petiot F, Delaunay P, Despres P, Lenglet A, Jourdain F, Leparc-Goffart I, Charlet F, Ollier L, Mantey K, Mollet T, Fournier JP, Torrents R, Leitmeyer K, Hilairet P, Zeller H, Van Bortel W, Dejour-Salamanca D, Grandadam M, Gastellu-Etchegorry M: First two autochthonous dengue virus infections in metropolitan France, September 2010.

    Euro Surveill 2010, 15:19676. PubMed Abstract | Publisher Full Text OpenURL

  10. Jansen CC, Beebe NW: The dengue vector Aedes aegypti: what comes next.

    Microbes Infect 2010, 12:272-279. PubMed Abstract | Publisher Full Text OpenURL

  11. Rigau-Perez JG, Clark GG, Gubler DJ, Reiter P, Sanders RJ, Vorndam AV: Dengue and dengue haemorrhagic fever.

    Lancet 1998, 352:971-977. PubMed Abstract | Publisher Full Text OpenURL

  12. Lambrechts L, Scott TW, Gubler DJ: Consequences of the expanding global distribution of aedes albopictus for dengue virus transmission.

    PLoS Negl Trop Dis 2010, 4:e646. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Hunter PR: Vector Borne Disease and Climate Change. In Encyclopedia of Environmental Health. Edited by Nriagu JNJ. Surrey, UK: Elsevier; 2011. OpenURL

  14. Rogers DJ, Randolph SE: Climate change and vector-borne diseases.

    Adv Parasitol 2006, 62:345-381. PubMed Abstract | Publisher Full Text OpenURL

  15. Degallier N, Favier C, Menkes C, Lengaigne M, Ramalho WM, Souza R, Servain J, Boulanger JP: Toward an early warning system for dengue prevention: modeling climate impact on dengue transmission.

    Clim Chang 2010, 98:581-592. Publisher Full Text OpenURL

  16. Hales S, de Wet N, Maindonald J, Woodward A: Potential effect of population and climate changes on global distribution of dengue fever: an empirical model.

    Lancet 2002, 360:830-834. PubMed Abstract | Publisher Full Text OpenURL

  17. Hopp MJ, Foley JA: Worldwide fluctuations in dengue fever cases related to climate variability.

    Clim Res 2003, 25:85-94. OpenURL

  18. Rogers DJ, Wilson AJ, Hay SI, Graham AJ: The global distribution of yellow fever and dengue.

    Adv Parasitol 2006, 62:181-220. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Thai KTD, Anders KL: The role of climate variability and change in the transmission dynamics and geographic distribution of dengue.

    Exp Biol Med 2011, 236:944-954. Publisher Full Text OpenURL

  20. Schmidt WP, Suzuki M, Thiem VD, White RG, Tsuzuki A, Yoshida LM, Yanai H, Haque U, Tho LH, Anh DD, Ariyoshi K: Population density, water supply, and the risk of dengue fever in Vietnam: cohort study and spatial analysis.

    Plos Med 2011, 8:e1001082. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Sriprom M, Chalvet-Monfray K, Chaimane T, Vongsawat K, Bicout DJ: Monthly district level risk of dengue occurrences in Sakon Nakhon Province, Thailand.

    Sci Total Environ 2010, 408:5521-5528. PubMed Abstract | Publisher Full Text OpenURL

  22. Caminade C, Medlock JM, Ducheyne E, McIntyre KM, Leach S, Baylis M, Morse AP: Suitability of European climate for the Asian tiger mosquito Aedes albopictus: recent trends and future scenarios.

    J R Soc Interface 2012, 9:2708-2717. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Reiter P: Climate change and mosquito-borne disease.

    Environ Health Perspect 2001, 109(Suppl 1):141-161. PubMed Abstract | PubMed Central Full Text OpenURL

  24. Colon-Gonzalez FJ, Fezzi C, Lake IR, Hunter PR: The effects of weather and climate change on dengue.

    PLoS Negl Trop Dis 2013, 7:e2503. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Hastie T, Tibshirani R: Generalized additive models.

    Stat Sci 1986, 1:297-318. Publisher Full Text OpenURL

  26. Guisan A, Edwards TC, Hastie T: Generalized linear and generalized additive models in studies of species distributions: setting the scene.

    Ecol Model 2002, 157:89-100. Publisher Full Text OpenURL

  27. Schimek MG: Semiparametric penalized generalized additive models for environmental research and epidemiology.

    Environmetrics 2009, 20:699-717. Publisher Full Text OpenURL

  28. Wood S: Generalized Additive Models: An Introduction With R. Boca Raton, FL, USA: CRC Press; 2006. OpenURL

  29. Brunkard JM, Cifuentes E, Rothenberg SJ: Assessing the roles of temperature, precipitation, and ENSO in dengue re-emergence on the Texas-Mexico border region.

    Salud Publica Mex 2008, 50:227-234. PubMed Abstract | Publisher Full Text OpenURL

  30. Colon-Gonzalez FJ, Lake IR, Bentham G: Climate variability and dengue fever in warm and humid Mexico.

    Am J Trop Med Hyg 2011, 84:757-763. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Hurtado-Diaz M, Riojas-Rodriguez H, Rothenberg SJ, Gomez-Dantes H, Cifuentes E: Short communication: impact of climate variability on the incidence of dengue in Mexico.

    Trop Med Int Health 2007, 12:1327-1337. PubMed Abstract | Publisher Full Text OpenURL

  32. The R Development Core Team: R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2010. OpenURL

  33. Bhaskaran K, Gasparrini A, Hajat S, Smeeth L, Armstrong B: Time series regression studies in environmental epidemiology.

    Int J Epidemiol 2013, 42:1187-1195. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Rockel B, Will A, Hense A: The Regional Climate Model COSMO-CLM(CCLM).

    Meteorol Z 2008, 17:347-348. Publisher Full Text OpenURL

  35. IPCC: Special Report on Emissions Scenarios: A Special Report of Working Group III of the Intergovernmental Panel on Climate Change. Cambridge; New York: Cambridge University Press; 2000. OpenURL

  36. IPCC: Climate Change 2007: The Physical Science Basis. Cambridge: Cambridge University Press; 2007. OpenURL

  37. Mennis J: Generating surface models of population using dasymetric mapping.

    Prof Geographer 2003, 55:31-42. OpenURL

  38. Dijkstra L, Poelman H: Regional Typologies: A Compilation. Regional Focus – A Series of Short Papers on Regional Research and Indicators Produced by the Directorate-General for Regional Policy. Brussels, Belgium: European Commission Directorate-General for Regional Policy; 2011. OpenURL

  39. Robinson AH: Elements of Cartography. 6th edition. New York: Wiley; 1995. OpenURL

  40. European Centre for Disease Prevention and Control: Development of Aedes Albopictus Risk Maps. Technical Report. Stockholm: European Centre for Disease Prevention and Control; 2009. OpenURL

  41. Benedict MQ, Levine RS, Hawley WA, Lounibos LP: Spread of the tiger: global risk of invasion by the mosquito Aedes albopictus.

    Vector Borne Zoonotic Dis 2007, 7:76-85. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Eritja R, Escosa R, Lucientes J, Marques E, Molina R, Roiz D, Ruiz S: Worldwide invasion of vector mosquitoes: present European distribution and challenges for Spain.

    Biol Invasions 2005, 7:87-97. Publisher Full Text OpenURL

  43. Fischer D, Thomas SM, Niemitz F, Reineking B, Beierkuhnlein C: Projection of climatic suitability for Aedes albopictus Skuse (Culicidae) in Europe under climate change conditions.

    Glob Planet Chang 2011, 78:54-64. Publisher Full Text OpenURL

  44. Knudsen AB, Romi R, Majori G: Occurrence and spread in Italy of Aedes albopictus, with implications for its introduction into other parts of Europe.

    J Am Mosq Control Assoc 1996, 12:177-183. PubMed Abstract OpenURL

  45. Rogers DJ, Suk JE, Semenza JC: Using global maps to predict the risk of dengue in Europe.

    Acta Trop 2014, 129:1-14. PubMed Abstract | Publisher Full Text OpenURL

  46. Hales S, Weinstein P, Woodward A: Dengue fever epidemics in the South Pacific: driven by El Nino Southern Oscillation?

    Lancet 1996, 348:1664-1665. OpenURL

  47. Bouzid M, Hooper L, Hunter PR: The effectiveness of public health interventions to reduce the health impact of climate change: a systematic review of systematic reviews.

    PLoS One 2013, 8:e62041. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  48. Sousa CA, Clairouin M, Seixas G, Viveiros B, Novo MT, Silva AC, Escoval MT, Economopoulou A: Ongoing outbreak of dengue type 1 in the Autonomous Region of Madeira, Portugal: preliminary report.

    Euro Surveill 2012, 17:pii=20333. OpenURL

  49. Randolph SE, Rogers DJ: The arrival, establishment and spread of exotic diseases: patterns and predictions.

    Nat Rev Microbiol 2010, 8:361-371. PubMed Abstract | Publisher Full Text OpenURL

  50. Murrell S, Wu SC, Butler M: Review of dengue virus and the development of a vaccine.

    Biotechnol Adv 2011, 29:239-247. PubMed Abstract | Publisher Full Text OpenURL

  51. Reiter P, Lathrop S, Bunning M, Biggerstaff B, Singer D, Tiwari T, Baber L, Amador M, Thirion J, Hayes J, Seca C, Mendez J, Ramirez B, Robinson J, Rawlings J, Vorndam V, Waterman S, Gubler D, Clark G, Hayes E: Texas lifestyle limits transmission of dengue virus.

    Emerg Infect Dis 2003, 9:86-89. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Angelini R, Finarelli AC, Angelini P, Po C, Petropulacos K, Macini P, Fiorentini C, Fortuna C, Venturi G, Romi R, Majori G, Nicoletti L, Rezza G, Cassone A: An outbreak of chikungunya fever in the province of Ravenna, Italy.

    Euro Surveill 2007, 12:E070906 070901. OpenURL

  53. Rezza G, Nicoletti L, Angelini R, Romi R, Finarelli AC, Panning M, Cordioli P, Fortuna C, Boros S, Magurano F, Silvi G, Angelini P, Dottori M, Ciufolini MG, Majori GC, Cassone A: Infection with chikungunya virus in Italy: an outbreak in a temperate region.

    Lancet 2007, 370:1840-1846. PubMed Abstract | Publisher Full Text OpenURL

  54. Clark DV, Mammen MP Jr, Nisalak A, Puthimethee V, Endy TP: Economic impact of dengue fever/dengue hemorrhagic fever in Thailand at the family and population levels.

    Am J Trop Med Hyg 2005, 72:786-791. PubMed Abstract | Publisher Full Text OpenURL

  55. Naish S, Dale P, Mackenzie JS, McBride J, Mengersen K, Tong SL: Climate change and dengue: a critical and systematic review of quantitative modelling approaches.

    BMC Infect Dis 2014, 14:167. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

Pre-publication history

The pre-publication history for this paper can be accessed here:

http://www.biomedcentral.com/1471-2458/14/781/prepub