Understanding the factors underlying the spatio-temporal distribution of infectious diseases provides useful information regarding their prevention and control. Dengue fever spatio-temporal patterns result from complex interactions between the virus, the host, and the vector. These interactions can be influenced by environmental conditions. Our objectives were to analyse dengue fever spatial distribution over New Caledonia during epidemic years, to identify some of the main underlying factors, and to predict the spatial evolution of dengue fever under changing climatic conditions, at the 2100 horizon. We used principal component analysis and support vector machines to analyse and model the influence of climate and socio-economic variables on the mean spatial distribution of 24,272 dengue cases reported from 1995 to 2012 in thirty-three communes of New Caledonia. We then modelled and estimated the future evolution of dengue incidence rates using a regional downscaling of future climate projections. The spatial distribution of dengue fever cases is highly heterogeneous. The variables most associated with this observed heterogeneity are the mean temperature, the mean number of people per premise, and the mean percentage of unemployed people, a variable highly correlated with people' s way of life. Rainfall does not seem to play an important role in the spatial distribution of dengue cases during epidemics. By the end of the 21st century, if temperature increases by approximately 3°C, mean incidence rates during epidemics could double. In New Caledonia, a subtropical insular environment, both temperature and socio-economic conditions are influencing the spatial spread of dengue fever. Extension of this study to other countries worldwide should improve the knowledge about climate influence on dengue burden and about the complex interplay between different factors. This study presents a methodology that can be used as a step by step guide to model dengue spatial heterogeneity in other countries.
Dengue fever is the most important mosquito-borne viral disease, with an estimated 50 million people being infected each year and 2. 5 billion people living in areas at risk of dengue worldwide [1]. The true burden of clinically apparent dengue could be twice as high, and the total burden of dengue fever infections could reach 390 million people when including asymptomatic cases [2]. Whereas only nine countries were affected by dengue epidemics in the 1970' s, more than a hundred countries are now reporting dengue outbreaks on a regular basis, making dengue fever the most rapidly spreading mosquito-borne viral disease in the world [1,2]. This rapid global spatial spread over the past 40 years probably results from recent socio-economic changes such as global population growth and uncontrolled urbanisation. Lack of effective mosquito control in endemic areas, increased international air traffic or decay in public health infra-structure in developing countries are also important factors that could explain the rapid regional spread of the disease [3–6]. However, in a given country where there are sufficient numbers of susceptible hosts, these factors need to be associated with suitable climate conditions before dengue fever can establish, since it is transmitted by mosquito species whose life cycle is influenced by temperature, humidity and rainfall [7–9]. Indeed, several studies have pointed out that the current geographic distribution of dengue fever or its vector worldwide could be predicted accurately based on climate variables using either statistical models [10] or deterministic models [11–13]. Other studies have pointed out that climate change could have profound consequences on the epidemiology of dengue fever, because increased temperature and rainfall could facilitate viral transmission and could lead to the geographic expansion of the mosquito species responsible for its transmission [11,14–17]. The complex interplay and relative importance of each factor in the occurrence and spread of dengue fever epidemics might differ from one country to another, depending on the specific climate conditions, cultural and socio-economic environment the virus circulates in [18]. Identifying the factors limiting dengue fever spatial spread at a national level could help understanding the worldwide pattern of dengue disease, could help predicting its future spatial distribution, and could provide national decisions-makers with useful information regarding the appropriate control measures to be implemented. Most studies trying to identify dengue risk factors spatially were performed at a city scale or a local scale (< 80 km) [19–34]. Among these studies, some have identified risk factors for the presence of Aedes mosquito species, such as socio-economic factors [29,30], proximity to specific plantations [28,32], proximity of potential breeding sites [22,28,29,32], or human behaviour [30,32]. Some studies highlighted the importance of human movement [23,26,31,33] or population immunity [34] in shaping the spatial transmission of dengue fever at small spatial scales. High dengue incidence rates have also been found in neighbourhoods with low social income [20,26], difficult access to piped water [21,26], or no implementation of mosquito protection measures [21,27]. Spatial analyses at a country or territorial scale (> 200 km) are scarce. Some of these studies focused on the spatio-temporal dynamics of the disease only [35–37] and proposed hypotheses about the underlying processes, but did not include analyses of potential explicative factors. To our knowledge, there are only five studies to date identifying and quantifying spatial risk factors for dengue at a “national” scale > ~200 km and < ~1000 km. Four studies identified temperature as having a major influence on the spatial distribution of locally acquired dengue cases [38–41], the last one did not assess the role of climate factors [42]. The role of other factors in the spatial distribution of dengue cases varied from one place to another. For example, in Australia (Queensland) [38] and Taiwan [40], rainfall seemed to play a minor role whereas in Brazil, rainfall played a major role [39]. In Taiwan [40] and Argentina [41], urbanization level was a key factor in dengue fever spatial distribution, and in one province of Thailand [42], the main factor identified was the proximity to major urban centres. No association was found with socio-economic covariates in Argentina [41] or Australia [38]. New Caledonia, where the present study takes place, has a unique situation: it is a developed insular territory located in the inter-tropical area of the South Pacific where the access to high quality data and the lack of terrestrial borders with other countries make it a natural laboratory to study dengue dynamics. A gross average of ten imported cases is detected each year by the Public Health authorities. However, large dengue epidemics develop only every three to five years, sometimes causing the circulation of the same serotype during two consecutive years [43]. A recent study analysing the temporal relationship between dengue epidemic occurrence and climate variables at an inter-annual scale showed that the development of an epidemic in New Caledonia needs precise climate conditions relying on both temperature and relative humidity [43]. The objectives of the present study were i) to characterise the spatial distribution of dengue cases in New Caledonia once an epidemic spreads over the territory; ii) to determine which of possible covariates are shaping the observed distribution; iii) to explore the potential spatial distribution of dengue cases under future climate projections. We present a complete methodology, from data collection, data transformation, variable selection, and application to future climate projections. We address a number of methodological issues such as spatial autocorrelation, correlation between explanatory variables, or potential non-linearity between epidemiological data and explanatory covariates.
New Caledonia is a French territory, located in the Pacific Ocean 1,500 km East of Australia. It is divided into 33 communes covering 18,576 km2. Out of the 245,344 inhabitants (2009), around 58% (147,365 people) live in Noumea, the main city, and its surroundings. The rest of the population is scattered in small towns of about 2,000 people, or live in rural areas, including traditional Melanesian settlements locally called “tribes” (Fig 1). Although the average population density outside Noumea is very low (5. 3 inhabitants per km2), local densities can be high as people gather in small settlements. New Caledonia is located at the limit of the tropical zone between latitudes 19° and 23° South. The East coast and the West coast are separated by a mountain range culminating at 1629 m. Climate is heterogeneous: the East coast and the southern tip of the main island get more rain than the West coast, as mountains provide a vertical lift to the warm and humid air brought by trade winds. Average rainfall range from 800 mm/year in some western weather stations to 3200 mm/year in the East. Temperature can drop below 10°C during the cool season on clear nights and sometimes rise above 35°C due to the influence of tropical air masses [45]. From an oversimplified point of view, there are three population groups, having different cultural and social habits: Melanesian people, people of French descent who migrated two hundred years ago, and people from various origins who migrated recently. Although the three groups are spatially partially mixed, Melanesian people live mostly on the East coast, whereas the second group live mostly on the West coast and the third group live mainly in Noumea. In New Caledonia, dengue represents a major public health problem with large epidemics affecting the territory every three to five years and involving a succession of all four serotypes [43,46–48]. Co-circulation of different dengue virus serotypes (DENV1-4) during major epidemics is rare, and has been observed only once (2009). Before 2003, vector control measures consisted in systematic chemical control of adult mosquitoes covering large areas during the warm season, independently from the occurrence of dengue cases. Since 2003, systematic spreading of adulticide has been stopped, and vector control measures include continuous large communication and prevention campaigns fostering source reduction aimed at all citizens, as well as focal chemical control of adult mosquitoes 100 m around declared cases within 24 h of notification. Public Health infrastructure is reliable, and the surveillance system for dengue fever has been efficient for many years. All people have access to medical care, even though people living in remote areas might have more difficulties to reach local health centres. As described below in the results section, temperature is a key factor determining dengue spatial variability over New Caledonia. We thus decided to explore the evolution of dengue average annual incidence rates during epidemics under changing climate conditions (considering all others variables as remaining constant), by applying the best explicative multivariable model with inputs from maps of temperature for the future (see methods/data/climate covariates: assessing the trends of future mean temperature in New Caledonia). Because the use of kernels in non–linear SVM models impairs predictions outside the observed range of explanatory variables, we built a linear approximation of the best SVM model on present observed data. The linear approximation consists of a simple linear model linking the two best explanatory variables to observed dengue age-standardised average (across epidemic years) annual incidence rates as the response variable. Normality and homoscedasticity of residuals were confirmed by the Shapiro-Wilks' test and the Bartlett' s test respectively [67]. To evaluate the error in incidence rates predictions due to the inter-GCM variability of mean temperature increase projections, we calculated, for each time-period and each scenario, the average annual incidence rates during epidemics as predicted by each GCM, and then calculated a standard deviation of predicted annual incidence rates across the different models.
Fig 3 shows that once an epidemic spreads over the territory, dengue cases are distributed heterogeneously. Mean annual age-standardised incidence rates across epidemic years range from 22 to 375 cases per 10,000 people per year, with a mean across communes of 168 cases per 10,000 people per year and a standard deviation across communes of 83 cases per 10,000 people per year. On average the East coast is more affected than the West coast. We can also see that the North-eastern corner of New Caledonia is heavily affected, with dengue incidence rates two to three times higher than in the rest of the territory. By definition, the average across epidemic years of age-standardised annual incidence rates reflects mainly the spatial pattern of severe epidemics, i. e. epidemics of years 1995,1998,2003 and 2009. During years 1995,1998 and 2003, the North-eastern corner was the most affected. During the 2009 epidemic, the most affected communes were Voh and Koné, on the West coast (see Fig 3 for the location of these communes), but the North Eastern corner was still severely affected [72]. The semi-variograms of dengue incidence rates did not reveal any significant spatial autocorrelation, whether they were calculated for each epidemic year separately or on the average incidence rates across years. This suggests that the local spread of dengue viruses around a case imported in a commune do not exceed the mean radius of a commune in New Caledonia, e. g. approximately 13 kilometres. Hence, we did not incorporate any spatial structure into the subsequent models. Table 1 shows Pearson' s correlation coefficients (rho) between each explanatory variable and dengue age-standardised annual incidence rates averaged across epidemic years. Dengue is spatially positively correlated with variables related to temperature and precipitation, but is negatively correlated with variables reflecting mean thermal range or extreme thermal conditions (see" Isothermality" ," Temp range" or the number of days when maximum temperature exceeds 32°C in January, February and March in Table 1). This suggests that, in a given commune, marked temporal variations of temperature is a factor limiting viral circulation. Based on the linear dependence measure of correlation, dengue is also more strongly associated with temperature than with precipitation. Socio-economic variables are highly spatially correlated to dengue average (across epidemic years) annual incidence rates. Variables reflecting people' s way of life (e. g. place of birth), local human density (e. g. mean number of people per household, percentage of premises under 40 m2), or human movement are more correlated with dengue average (across epidemic years) annual incidence rates than variables related to the housing type (e. g. premises with inside toilets) (absolute value of rho up to 0. 75 for the former and 0. 58 for the latter). In particular, the fact that the place where people were born is spatially significantly associated with dengue fever incidence rates (correlation coefficient around 0. 5 for people born in New Caledonia and– 0. 5 for people born elsewhere) whereas the type of premise is not (absolute correlation coefficient lower than 0. 3 for variables describing access to water or electricity) suggests that individual behaviours have a stronger influence on incidence rates than local housing conditions. Fig 4 shows the PCA results. For clarity reasons, we only show the results of PCA performed on the variables most spatially correlated with dengue average (across epidemic years) annual incidence rates, with an absolute Pearson correlation coefficient over 0. 6 for socio-economic variables, and over 0. 4 for climate variables (these thresholds were selected after verifying that they did not modify the variable pre-selection results). PCA of climate variables (Fig 4A) shows that in New Caledonia, temperature is the factor accounting for most of the spatial climatic variability among communes. Temperature is highly correlated with the first PCA axis which represents 68% of the total climatic variance. Temperature and rainfall are not spatially correlated at the commune level. In each group of temperature or rainfall variables, the variables most spatially correlated with dengue average (across epidemic years) annual incidence rates were the average mean temperature (Mean temp) and the mean daily rainfall during the wettest quarter of the year (Wettest quarter) (see Table 1). In addition to these two variables, we decided to keep a third variable, the average daily rainfall, for further statistical modelling, as this variable is more easily available in other countries or climate model simulations. Fig 4B shows that the spatial variability of socio-economic factors mainly reflects the spatial distribution of people with different cultural habits. Communes where a high proportion of inhabitants live in a tribal way, in small premises, with few means of transportation and a high percentage of unemployment are opposed to communes where many people live a western way of life, in permanent buildings, using air conditioning and getting around using cars. Even though the number of people per premise seems to be correlated with the proportion of people living in tribes, we kept this variable as it stands out of the cluster of variables representing the way of life. We thus decided to keep the percentage of unemployed people and the mean number of inhabitants per housing as representative of socio-economic factors for further statistical modelling. Table 2 shows the RMSE of the optimised models built on all possible combinations of one, two or three of the five selected explanatory variables (Mean temperature, daily rainfall averaged over the wettest quarter, average daily rainfall, number of people per household and fraction of unemployed people). When looking at univariable non-linear SVM models, the best variable explaining the spatial heterogeneity of dengue average (across epidemic years) annual incidence rates is the percentage of unemployed people per commune. The second most important explanatory variable is the mean temperature. Rainfall is the least explanatory variable of those selected for multivariable regression modelling. Moreover, the RMSE of models based on observed rainfall almost equal the initial standard deviation (across the territory) of dengue average annual incidence rates, which means that rainfall are poor predictors of dengue average annual incidence rates during epidemic years. The relationship between dengue average annual incidence rates and each of the explanatory variables is linear, except for the fraction of unemployed people (S4 Fig). When looking at the spatial structure of dengue average annual incidence rates predicted by SVM models based only on one of the selected variables (S5 Fig), we see that temperature captures mainly the South to North gradient of increasing incidence rates (S5B Fig) whereas socio-economic variables captures the spatial heterogeneity between the West coast and the East coast (S5E and S5F Fig). Temperature seems to have no influence in communes located below 21°S (S5B Fig). All models based on two explanatory variables and including at least one variable related to rainfall (best RMSE of ~58 cases per 10,000 people per year) performed worse than the best univariable model (RMSE of ~53 cases per 10,000 people per year). This suggests that in New Caledonia, rainfall has little influence on the spatial variability of dengue viral circulation at the commune level. Models combining two explanatory variables (excluding rainfall) performed better than models based on only one variable. The addition of a third explanatory variable did not improve significantly model performances. Hence we focused our attention on models combining two explanatory variables. The best explicative model is a model predicting increasing average annual incidence rates during epidemics in communes where the mean temperature and the mean number of people per premise increase (see Fig 5). The influence of these two variables on the spatial structure of dengue incidence rates is close to linear as shown by almost parallel contour lines on Fig 5A. This model accurately predicts the sharp mean increase in incidence in the three communes of the North East of New Caledonia (Hienghène, Ouégoa and Pouébo). The maximal error of the model is observed for Farino (West coast), which is the only commune where all inhabitants live at an altitude higher than 200 m above sea level. Fig 6 shows the spatial structure of observed average (across epidemic years) annual incidence rates (Fig 6A) and of average annual incidence rates as predicted by the best SVM model based on two explanatory variables (Fig 6B). This model captures the observed spatial heterogeneity in average annual incidence rates between the East coast and the West coast, as well as the sharp increase in the three communes of the North East. Table 3 shows, for both climate change scenarios, the average increase of mean temperature for the two selected 20-year periods compared to the 1980–1999 historical simulations. All models predict that the mean temperature will increase over time, with projections being more pessimistic for RCP 8. 5 simulations. The CMIP5-AR4 inter-model variability in temperature increase is presented in Table 3 and S3 Fig. According to the RCP 8. 5 scenario, temperature could increase by more than 3°C by the end of the next century, with a standard deviation across models of only 0. 6°C, showing the strong coherency in different model projections. Fig 6 shows a comparison of the average (across epidemic years) annual dengue incidence rates predicted by the SVM model (panel 6B) or the linear model (panel 6C). The SVM model performs slightly better than the linear one: the correlation coefficient between observed and predicted incidence rates are 0. 89 (SVM) and 0. 85 (linear), and the RMSE are 42 and 43 cases/10,000 people/year respectively for the SVM and the linear model. The low RMSE of the linear model (~43 cases per 10,000 people per year) shows that the linear model based on the two best explanatory variables is suitable. The Shapiro-Wilks and the Bartlett' s test confirmed the normality and homoscedasticity of residuals. Fig 6D and 6E show the potential future spatial distribution of dengue incidence rates during epidemics according to the RCP 4. 5 and RCP 8. 5 emission scenarios. By the end of the century, dengue incidence rates during epidemic years could reach a maximum of 378 cases per 10,000 people per year in the most affected commune under the RCP 4. 5 scenario (Fig 6D), and 454 cases per 10,000 people per year in the most affected commune under the RCP 8. 5 scenario (Fig 6E). Under the RCP 8. 5 scenario, communes at low risk now might experience a sharp increase in dengue incidence rates during epidemic years from 64 to more than 200 cases per 10,000 people per year. According to RCP 8. 5 climate projections, the average (across communes) dengue mean annual incidence rates during epidemic years could raise by 29 cases per 10,000 people per year for the 2010–2029 period, and by 149 cases per 10,000 people per year for the 2080–2099 period, almost doubling dengue burden in New Caledonia by the end of the century (Table 3).
The spatial association found between temperature and dengue incidence rates during epidemics in New Caledonia can be explained by the influence of temperature on the life cycle of the mosquito transmitting the virus in New Caledonia, Aedes aegypti. High temperatures increase the productivity of the breeding sites through an acceleration of the metabolism of the mosquito, and a faster development of the micro-organisms the larvae feed on, resulting in a higher vector density even with the same number of breeding sites [7,73–75]. High temperatures also speed up the extrinsic incubation period [7,76,77], with the effect that an increased proportion of females Ae. aegypti can reach the infectious stage before dying. Finally, warmer temperatures accelerate the mosquito gonotrophic cycle, and make females Ae. aegypti more aggressive [7,74,78–80], increasing the biting rate and the frequency of potential transmission of viral particles to susceptible hosts. Regarding the effect of increasing temperatures on the mortality of Ae. aegypti adult mosquitoes, a review of 50 field mark-release-recapture studies has shown that in the field, unless temperatures become extreme (over 35°C or less than 5°C), temperature has little effect on daily mortality rate [81], highlighting the central importance of the length of the extrinsic incubation period in the ability of adult mosquitoes to transmit dengue viruses. In Noumea, the main city, precise climate variables and important thresholds values have been identified as necessary conditions to trigger an epidemic (e. g. number of days when maximal temperature exceeds 32°C in January/February/March, and number of days when maximal relative humidity exceeds 95% during January [43]). At the scale of the entire territory, we found that the spatial distribution of dengue cases during epidemic years is strongly influenced by the average mean temperature. These results suggest that temperature has a major role in dengue dynamics in an insular territory characterised by climate seasonality. However, we did not find a strong association between the spatial distribution of dengue cases during epidemics and average rainfall or with the number of days when maximal temperature exceeds 32°C. The variables influencing either the triggering of an epidemic [43] or its spatial distribution are not the same. Our findings highlight the complexity of studying and understanding dengue dynamics, the importance of well separating the two epidemiological processes of epidemic triggering in a susceptible population, and its intensity once it has started by clearly defining the modelling target (incidence rates for epidemic intensity, or dummy variables for epidemic triggering), and the importance of well defining the scale of study (temporal evolution, or spatial distribution). The positive association found between the mean temperature and dengue incidence rates is consistent with the one found in previous studies having analysed the spatial distribution of dengue cases at spatial scale > 200 km [38–42]. In these studies as well as in ours, all regions were located between 10° and 25° of latitude, at the fringe of the tropical area, except Argentina, where the region studied extends to 35° South. In the 10° ˗ 25° latitudinal band, annual mean temperature lies in a range of temperature where the life cycle of the mosquito is very sensitive to temperature changes [7]. Some Aedes species, including Ae. aegypti, are able to breed in very small amounts of water, e. g. snails’ shells. Rainfall can play a role in dengue transmission cycle by filling up potential breeding sites [7], thus influencing the vector density. Rainfall also increases the relative humidity, which extends the mosquitoes’ lifespan and therefore the likelihood of those who had an infectious blood meal to reach the infectious stage. However, our study suggests that in New Caledonia, there is no strong association between rainfall and the spatial distribution of cases during epidemics. A plausible explanation can be the multi-factorial nature of dengue fever, and the relative influence each factor plays on dengue dynamics: despite suitable rainfall conditions, dengue might not circulate well if other factors are limiting dengue viral circulation, such as some human behaviour influencing the contact between vector and host. This aspect has been highlighted very clearly in the United States [27]. Some studies have found that the effect of rainfall on vector density can be modulated by human activities such as water storage practices [82]. However, in New Caledonia, we are not aware of specific practices to store water that could explain the lack of association between rainfall and virus circulation intensity. Another potential explanation could be that in dry areas, breeding sites are filled up by other non-climatic mechanisms, such as automatic irrigation or plant watering. Worldwide, the spatial association between rainfall and the spatial distribution of dengue cases at a “national” scale (> 200 km) is not as clear as the one for temperature: one spatial study did not find any association between dengue incidence rates and rainfall [40], whereas two others did [38,39]. Other factors influencing dengue transmission (e. g. anthropogenic factors influencing the availability of filled breeding sites) and not included in the different studies might blur the rainfall signal. It would be interesting to perform the same kind of multi-factorial spatial analysis in areas of epidemic or endemic transmission located closer to the equator, where the mean temperature is higher, to see what climatic factors impact the spatial distribution of cases. This kind of study could help understand better the complex interplay between the different factors (climate, socio-economic, immunologic, viral, entomologic…) associated with dengue fever transmission. Regarding the link between socio-economic variables and dengue incidence rates during epidemics, a limitation of this study is the absence of historical time series of socio-economic variables. We then had to assume that the data retrieved from the 2009 census is representative of the mean socio-economic spatial pattern over the epidemic years of the 1995–2012 period. As there has been no major historical event leading to population migration in New Caledonia during this time period and as socio-economic variables represent mainly people’s way of life, we think this assumption is realistic. Our results are consistent with previous studies that have pointed out the importance of socio-economic factors on the spatial distribution of dengue cases, whatever the spatial scale studied: national (>200 km) [39–42] or local (<10 km) [20,26,29,83]. The spatial association between the percentage of unemployed people and dengue in New Caledonia cannot be interpreted in terms of lack of economic activity only, as shown by the PCA on socio-economic factors. This variable we selected as input for the models is highly correlated with other variables reflecting the way of life, socio-economic and cultural differences existing in New Caledonia, which are in turn highly correlated to housing type. Therefore, at this spatial scale in New Caledonia, it is difficult to statistically differentiate the role played by human behaviour, human activity or housing type in dengue fever transmission. However, those three factors influence the contact rate between viraemic patients or susceptible hosts on one hand, and mosquitoes on the other hand. This highlights the importance of limiting the contact between humans and vectors and should lead local authorities to strengthen communication campaigns about personal protection measures towards populations at risk. Regarding the spatial association found between the fraction of unemployed people (i. e. people’s way of life) and dengue incidence rates during epidemics in New Caledonia, it is interesting to point out that on the East coast, a larger fraction of inhabitants are Melanesian people living in tribes, whereas on the West coast, the majority of people are people from French descent having a western way of life. It would be interesting to perform sociologic studies to precisely identify which human behaviour leads to an increased risk of catching dengue fever. Such information would be useful to define communication messages towards at risk populations. The spatial association between the number of people per household and dengue incidence rates can be explained by the short flight range of Ae. aegypti mosquitoes. These mosquitoes are often captured in the very house where they emerged or in the neighbouring houses, flying an average of 40 to 80 m during their life [84–87]. Hence, dengue outbreaks involving Ae. aegypti as the main vector are known to be highly spatially focal, with dengue cases usually clustering within 200 m to 800 m of each other [23,33,34,88–94]. Our results suggest that, in New Caledonia, dengue cases probably cluster within houses. Sick people should protect themselves until they are no longer vireamic to avoid human to mosquito transmission, and people living around a case should protect themselves to avoid getting infected while infectious mosquitoes are still active in the neighbourhood. Taking such individual actions could reduce the intensity of dengue transmission and reduce dengue burden over the territory. This message could be strengthened in the recommendations given by the authorities. The results about climate change must be interpreted keeping in mind that they represent a climate risk only, and that the spatial association between dengue incidence rates during epidemics and temperature might change over time depending on socio-demographic changes, or changes in dengue control strategy. Assuming all other factors remain constant in time, our results suggest that Public Health authorities can expect the dengue burden to raise significantly during the next century over the territory, and can expect the dengue spatial range to increase. As the GCM projections are spatially homogeneous over the territory, and as the model used to predict dengue incidence rates in the future is linear and is based on only one climate variable, the predicted absolute increase in dengue incidence rates is currently the same for all communes. This highlights the need for spatially downscaling GCM projections to gain a better understanding of the impact of climate change in the future. Communes that are already severely affected by dengue epidemics will have to prepare to face higher burden of dengue fever. For communes that are at low risk now, we can see that in the future they might be affected as severely as communes at high risk now. These communes might not be prepared now to face severe epidemics of dengue fever, and they will probably need support for adaptation. As said earlier, the positive association found between temperature and dengue incidence rates during epidemics for mean temperatures ranging from 22°C to 25°C can be explained by the effect of temperature on the mosquito life cycle and duration of extrinsic incubation period. Here, by applying a statistical model built using current observed temperature to future projections, we make the assumption that the biological effect of temperature on the mosquito life cycle and on the extrinsic incubation will remain the same under the range of temperature that might be observed in the future. For most parameters influencing transmission, this assumption is reasonable. For example, we know that in Thaïland, for DENV-2, the extrinsic incubation period is reduced from 15 days at 30°C to 7 days at 32–35°C [77], which is in support of increasing temperatures inducing an increase in dengue incidence rates under future climate. However, because we used a statistical model, we were not able to incorporate the known negative effect that an increase in temperature might have on dengue transmission when temperature reaches extremes. For example, a review of fifty mark-release-recapture studies has shown that the survival and longevity of Ae. Aegypti mosquitoes is highly reduced when temperatures exceed a threshold, which might be around 35°C [81]. It would be interesting to develop models that are able to integrate these negative effects in the future in order to gain a better understanding of the effects of climate change on dengue transmission. Here we used data collected routinely by the Direction of Sanitary and Social Affairs. As any surveillance system, it is highly probable that not all dengue cases have been recorded. However, in New Caledonia, the data is of high quality, and the spatial standardisation of the surveillance system (i. e. all the actions taken to be able to compare data collected by different people, at different places [95]) is good, which means that the proportion of cases that are not recorded by the surveillance system are probably comparable from one commune to another. Hence, maps of incidence rates calculated from routinely collected data can be used to study the spatial variability of true incidence rates. To calculate mean incidence rates, we have used the consultation date of cases. Consultation can occur 1 to 5 days after the onset of symptoms, and the incubation period lasts 4 to7 days on average [3], which means that the consultation date can differ from one to two weeks from the date of infection. This loss of temporal precision is not important here to calculate maps of incidence rates, as for each commune, we have averaged incidence rates across many years. As in any epidemiological spatial study using routinely collected surveillance data, it is possible that some spatial bias has been introduced due to the fact that the spatial data recorded is the commune where people live, which can sometimes differ from the commune where they got infected. Another limitation of using routinely collected data is that only clinically apparent cases are recorded, dismissing clinically inapparent cases, whose proportion can vary in time and space [96,97]. A seroprevalence survey is currently undertaken by the Public health authorities. It will be interesting to compare the spatial distribution of seroprevalence to the spatial distribution of average incidence. The analysis of the spatial pattern of infectious diseases, in relation with environmental or socio-economic factors raises a number of methodological issues, such as the presence of spatial auto-correlation, the spatial scale of aggregation of the data, the existence of possible non-linear links between the response and the explanatory variables, or the presence of multi-collinearity between the response variables. Most issues have already been addressed in the past, and solutions already exist to handle them. For example, the reviews by Dormann et al. deal with multi-collinearity [66] or spatial-autocorrelation [98]. The main issue we have been confronted with in our study was the spatial upscaling of meteorological data observed at precise locations to the same spatial level of aggregation as the epidemiological data. In existing spatial studies of dengue fever at a national scale, authors have geo-spatially interpolated climate variables on regular grids using kriging methods, and have averaged gridded values over a given administrative division [38–41]. This approach has two drawbacks in New Caledonia. Simple kriging models do not take into account the potential elevation between two given points, leading to biased estimates of temperature in mountainous regions. Moreover, the traditional approach used in climatology, which consist in aggregating temperature over grid points taken uniformly over the whole aggregative area makes the implicit assumption that people at risk are distributed homogeneously over the aggregative area. This is particularly problematic in New Caledonia where large areas are not inhabited. In our approach, as epidemiological data are collected at the individual level, we tried to estimate the climate conditions for each individual (and therefore for the mosquitoes surrounding each individual). However, the algorithm used introduces some noise, due to the fact that the weather stations are sometimes kilometres away from some towns or tribes. High spatial resolution climate data obtained from high resolution modelling of atmospheric conditions could be used, but some noise will be introduced by the modelling error compared to the observed data. This issue needs further attention in the future to increase the quality of spatial epidemiological and environmental studies. Some factors that could influence the spatial distribution of dengue cases during epidemics have not been taken into account in this study: the location of the first cases introduced each year, the spatial variability in population immunity, viral factors such as the serotype circulating, or factors associated to the mosquito such as the spatial variability in vector competence, or dengue vector control measures. We decided not to include the serotype, as we performed the analysis on averages over several years, and as, except in 2009, there was no co-circulation of different serotypes over the territory. Therefore, spatial differences in the level of viral circulation cannot be associated with genetic differences between serotypes. A territorial seroprevalence survey to assess population immunity has been implemented recently in New Caledonia, but data are not available yet. It could be interesting to include environmental variables derived from GIS data or remote-sensing in this kind of study. For example, GIS data about built areas could be used to create indicators of the proximity of houses to reflect the fragmentation of the Ae. aegypti habitat per commune, given that fragmentation of this habitat could potentially slow down viral circulation. The amount of vector control effort implemented in each commune is heterogeneous on the territory, as this activity falls within the commune’s authority, and each commune is free to implement or not the territorial guidelines. The vector control effort in each commune is thus difficult to quantify, and data were not available yet at the time of analysis. As soon as these data will be collected by local authorities, they could be incorporated in the modelling process to assess the efficiency of vector control measures. The study we present here is about the spatial heterogeneity of dengue incidence rates across epidemic years, independently of the inter-annual variability of dengue incidence rates from one epidemic year to another. The mean spatial pattern studied is very robust to changes in the definition of an epidemic year, as severe epidemics will always be considered in the calculation of the mean, whatever the threshold used to distinguish epidemic from non-epidemic years. It would be interesting to know whether the spatial association found here between the severity of dengue epidemics, temperature, local people’s density and people’s way of life is consistent through time or not, and to identify the factors associated to the temporal variability of spatial patterns. Two other viruses transmitted by Ae. aegypti like dengue virus caused outbreaks recently in New Caledonia. Chikungunya virus has been introduced on four occasions since 2011 but in each case, the outbreaks were limited to a few cases in Noumea and surroundings. Conversely, Zika virus caused large epidemics over the territory in 2014 and 2015, with more than 1,500 confirmed cases and more than 11,000 estimated cases. Although these viruses are transmitted by the same mosquito as dengue fever, no sufficient data are available to know if the socio-economic and climatic factors driving epidemics are the same. It is likely that local vector competence and population immunity represent major limiting factors. Although dengue has caused major outbreaks in NC in 2013, chikungunya viruses have only caused a limited number of cases for reasons that remain unexplained today and despite the competence of local Ae. aegypti for chikungunya virus transmission [99]. It is likely that climatic factors and interactions between viruses circulating together between human-hosts and mosquito-vectors influence the epidemiology of arboviruses in New Caledonia. A comparative analysis of the spatio-temporal distribution of these three arboviruses in an insular territory accommodating only Ae. aegypti represents an important issue to understand and predict outbreaks.
Dengue fever is the most important viral arthropod-borne disease worldwide and its geographical expansion during the past decades has been of growing concern for scientists and public health authorities because of its heavy sanitary burden and economic impacts. In the absence of an effective vaccine, control is currently limited to vector-control measures. In this context, understanding the sociologic, entomologic and environmental factors underlying dengue dynamics is essential and can provide public health authorities with sound information about control measures to be implemented. In this study, we analyse socio-economic, climatic and epidemiological data to understand the impact of the studied factors on the spatial distribution of dengue cases during epidemic years in New Caledonia, a French island located in the South Pacific. We identify at risk areas, and find that temperature and people’s way of life are key factors determining the level of viral circulation in New Caledonia. Hence, communication campaigns fostering individual protection measures against mosquito bites could help reduce dengue burden in New Caledonia. Using projections of temperature under different scenarios of climate change, we find that dengue incidence rates during epidemics could double by the end of the century, with areas at low risk of dengue fever being highly affected in the future.
Socio-economic and Climate Factors Associated with Dengue Fever Spatial Heterogeneity: A Worked Example in New Caledonia