Research | Open | Published:
Influence of climate and river level on the incidence of malaria in Cacao, French Guiana
Malaria Journalvolume 10, Article number: 26 (2011)
The epidemiological profiles of vector-borne diseases, such as malaria, are strongly associated with environmental conditions. An understanding of the effect of the climate on the occurrence of malaria may provide indirect insight into the anopheles mosquito vectors endemic to a particular region. The association between meteorological and hydrographical factors and the occurrence of malaria was studied in a village in French Guiana during an epidemic caused essentially by Plasmodium vivax.
A cohort of confirmed cases of P. vivax malaria occurring between 2002 and 2007 was studied to search for an association between the number of new infection episodes occurring each month, mean, maximum and minimum monthly temperatures, cumulative rainfall for the month and the mean monthly height of the river bordering the village, with the aid of time series. Cross-correlation analysis revealed that these meteorological factors had large effects on the number of episodes, over a study period of 12 months.
Climatic factors supporting the continuance of the epidemic were identified in the short-term (low minimum temperatures during the month), medium-term (low maximum temperatures two months before) and long-term (low maximum temperatures nine months before and high lowest level of the river 12 months before). Cross-correlation analysis showed that the effects of these factors were greatest at the beginning of the short rainy season.
The association between the river level and the number of malaria attacks provides clues to better understand the environment of malaria transmission and the ecological characteristics of the vectors in the region.
Malaria may be seen as the combination of a Plasmodium, an Anopheles mosquito and a human host. Factors modifying the presence, survival or abundance of any one of these actors are likely to favour the maintenance or breaking of the parasite cycle. Environmental conditions are among the factors likely to affect this parasite cycle. The temperature of the water collecting in hollows, the presence of vegetation or predators and altitude may all favour Anopheles mosquitoes or limit their development within a given geographic zone [1, 2].
The ecological characteristics of the vectors of malaria in French Guiana have not been completely elucidated [3, 4]. In the absence of entomological data, the goal of this study was to investigate the environmental factors conditioning the encounter between humans and parasites, focusing on the incidence of malaria episode in the village of Cacao in French Guiana. Over a period of 14 years, there were no cases of malaria in this village. Then, between 2002 and 2007, an epidemic caused essentially by Plasmodium vivax occurred, affecting almost 43% of the inhabitants of the village.
To better understand the determinants of this epidemic, this manuscript presents results studying the association between meteorological and hydrographical factors and the occurrence of malaria among the inhabitants of Cacao.
A cohort of patients from the Cacao health centre presenting between January 1st 2002 and December 31st 2007 with an episode of vivax malaria confirmed by biological tests was studied. The only health facility within a radius of 80 kms was located in the village. A malaria case corresponding to an attack of malaria was defined as an episode of fever (fever recorded or history of fever in the previous 48 hours, associated to the finding of asexual forms of Plasmodium in the blood using Giemsa-stained thick and thin blood smear. in the absence of any other clear aetiology. The parasitological results were validated by a biologist from Cayenne University Hospital Parasitology and Mycology Laboratory (LHUPM). Only patients residing permanently in the village during the study period were included. To ensure data quality, the population that was living temporarily in Cacao (seasonal workers, visitors, tourists...) was not included in the analysis. Therefore, the sample considered in the analysis, remained stable. Patients were included in the cohort from the day the positive sample for Plasmodium vivax was obtained. The primary endpoint was the occurrence of a new episode of P. vivax malaria (detected by the health centre in Cacao and subsequently biologically confirmed, excluding cases of treatment failure). Uncomplicated malaria cases were treated with chloroquine. None of the patients received 8-aminoquinoline. Relapses were excluded by applying the rule of three months as shown recently in French Guiana for local strains of P. vivax. Indeed, in a study carried out in Camopi, a village located in an endemic area, the secondary P. vivax attack was considered as a relapse when it occurred during the first three months after the first attack .
The association between climatic factors and the monthly incidence of P. vivax malaria was investigated by a time-series model of the autoregregressive integrated moving average (ARIMA) type. The series was non-stationary and the variance grew exponentially. First, logarithmic transformation ensured the normality and homogeneity of the variance of the residuals. Then, to make the series stationary, a first-order differentiation was performed . The parameters of the ARIMA and the seasonal parameters were determined with the aid of Akaike's information criterion (AIC), the autocorrelation function (ACF) and the partial autocorrelation function (PACF). The model chosen was ARIMA(1,1,1)(0,1,1)12. The ARIMA residuals were not (auto-) correlated and were normally distributed according to the residuals plot and the autocorrelogram of residuals. The Ljung-Box Q test was also applied to confirm that the residual series were white noise.
The meteorological factors studied (measured at the Météo France weather station at the village) were selected on the basis of the specific climatic features in the study zone. This zone has a humid tropical environment and is covered by a dense forest with large amounts of aerial biomass. Humidity is about 90% all year round. The village of Cacao is situated about 80 km east of Cayenne, the main city of French Guiana (Figure 1). The mountain tops are low, at an altitude of only 30 to 400 m. The village of Cacao is bordered to the North West by a river, the Comté. The region has four seasons: i) a short rainy season (November to January); ii) a "short summer" (February and March); iii) a long rainy season (mid-March to mid-July), and iv) the dry season (August to November).
The following factors were studied: i) mean minimum temperature for each month; ii) mean maximum temperature for each month; iii) cumulative rainfall for each month; iv) the lowest level of the Comté for each month and v) the highest level of the Comté for each month. These climatic data were measured daily and introduced into the model in the form of the mean monthly value for a period of up to 12 months. These factors were tested in the univariate analysis, with all factors identified as significant at the 10% level being introduced into the multivariate model. A stepwise descending approach was used to obtain the final model. The residuals of the models were visualised graphically (by autocorrelation and non systematic tendency). Spearman's rank correlation analysis was carried out to assess the association between the occurrence of malaria attacks and meteorological factors.
Statistical significance was set at p < 0.05. The data was analyzed using STATA 10.0 (STATA Corporation, College Station, Texas).
The village of Cacao is inhabited by 839 people, 359 of whom presented at least one attack of P. vivax malaria between 2002 and 2007. The Cacao health centre detected 817 attacks of P. vivax malaria and 40 attacks of Plasmodium falciparum malaria among the inhabitants of this village during this epidemic. The incidence of attacks due to new infections (616 episodes) followed a seasonal pattern over the 12 months of the year, with a peak in April (Figure 2).
Between 2001 and 2007, the months of March, April, May and June had the highest levels of rainfall during the year (cumulative rainfall levels greater than 400 mm per month) and, in these months, the height of the river exceeded 4 m. Minimum and maximum temperatures remained fairly stable, with a slight increase in maximum temperatures observed between August and November.
The mean minimum temperature fell below 21°C between June and September 2003, reaching a trough of 19.4°C in July 2003. The mean maximum temperature exceeded 35°C during the following periods: September and October 2002, from August to November in 2003 and September and October in 2004 and 2005, with a peak at 36.3°C during August 2003 (Figure 3b). Precipitation levels followed a strongly seasonal pattern (Figure 3a) during the study period. Rainfall levels were highest in April 2002, May 2003 and April 2005 (cumulative precipitation of 893.2 mm for this month).
The introduction of meterological factors into the univariate ARIMA model (1,1,1)(0,1,1)12 identified nine relevant factors: i) mean minimum temperatures at time t and t-12 months (factor 1 and 2); ii) mean maximum temperature at t-1, t-2 and t-9 months (factor 3,4,5); iii) the minimum level of the Comté at t-11 and t-12 months (factor 6 and 7); and iv) the maximum level of the Comté at t-8 and t-12 months (factor 8 and 9).
The multivariate ARIMA model obtained with the inclusion of these nine variables, applying a stepwise descending procedure, is summarized in Table 1. The incidence of a malaria episode at time t was positively correlated with the minimum level of the Comté at t-12 months and inversely correlated with the minimum temperature at time t and with maximum temperature at times t-2 and t-9 months.
The relationship between the incidence of malaria and the minimum temperature (Figure 4a) was strongest in December (rho = -0.7143; p = 0.1108). Correlations between the number of attacks and maximum temperatures two months before (Figure 4b) and nine months before (Figure 4c) were strongest in January (rho = -0.71; p = 0.117 and rho = - 0.62; p = 0.1911, respectively). The same was true for the minimum level of the Comté river (Figure 4d; rho = 0.53; p = 0.27)
This study was carried out on a small cohort, limiting the power of statistical analyses, including cross-correlations in particular. However, it should be noted that the medical data obtained are exhaustive for the village, homogeneous (malaria caused by P. vivax only) and high quality (identical diagnostic and therapeutic management over a period of five years).
Very little entomological knowledge is currently available for the vector of malaria at Cacao . Vector control programmes are, therefore, based on data from the other countries of South America, for Anopheles darlingi[7–11]. The lack of local specificity of the knowledge regarding the vector results in intervention programmes needing to be massive, with poor cost/efficacy ratios.
This study could indirectly help to improve this ratio, by providing clues to better understand the transmission of the parasite between vector and host in this region. Indeed, the results of this study showed certain climatic factors could participate in the maintenance of the epidemic in the short-, medium- or long-term.
First, the positive association between the river level and the incidence of malaria could be explained by the location of potential larval habitats. The banks of the Comté river are encumbered with voluminous tree roots and are, therefore, not straight. During periods of high water, these features favour the creation of pockets of water with weak currents, bordered by vegetation and ideal for the development of anopheles mosquito larvae [11–13]. The persistence of a minimum level of water in the river ensures the maintenance of these larval habitats and adequate development of the surrounding vegetation. Without entomological evidence, this hypothesis could not be confirmed. However, others authors in America have also proposed this. In the north of Peru, a number of climatic factors (temperature and rainfall figures) have been studied, but only the level of the river Nanay was found to be significantly associated with the incidence of malaria due to P. vivax. Similarly, a study in Suriname identified associations between mosquito density and the level of the river . In Brazil, flooding of the fields bordering the river Branco has also been suggested as a reason for increases in the number of An. darlingi bites per person during the rainy season, close to the river . In Venezuela, the abundance of An. darlingi was positively correlated with the maximum level of the Orinoco river because of the creation of lagoons during the river overflow .
The time lag between the occurrence of malaria and the associated meteorological factors, such as maximum temperature nine months previously and the level of the Comté 12 months previously, would suggest that meteorological conditions in a given year may affect malaria in the following year. With persistent larval habitats in all seasons due to sufficient river levels and low temperatures, the transmission of malaria is maintained into the next year. These environmental factors could create the appropriate conditions to preserve mosquito larvae with sufficient humidity in the dry season. Another hypothesis would be the need for a number of parasitic amplification cycles to make malaria transmission patent. The long-term impact of climate on malaria has already been reported in the literature. Indeed, a seven-month delayed effect of fog was associated with malaria incidence in a Chinese study. The authors explained this delay with the same biological hypothesis . Similarly, significant increases of malaria mortality and morbidity have been showed in years following recognized El Niño events .
The intertropical convergence zone (ITCZ) or meteorological equator determines the seasons in French Guiana. It descends the northern hemisphere, crossing French Guiana in mid-December and bringing with it cumulonimbus clouds, which generate rainfall. Cross-correlation analysis showed that the time at which temperatures and the level of the Comté river seemed to influence the occurrence of cases coincided with the period during which the ITCZ passes through French Guiana for the first time. The resulting increase in precipitation, with high humidity, favours an increase in the number of mosquito breeding sites and the development of the aquatic stages of the mosquito .
Finally, an increase in temperature (both minimum and maximum) has often been identified as a factor accelerating the gonotrophic cycle of anopheles mosquitoes, decreasing the interval between egg-laying episodes and increasing vector abundance [20–23]. The development cycle of the parasite is also accelerated by high temperatures. However, if the temperature is too high, it may also jeopardise the survival of adult mosquitoes and dry up breeding sites more rapidly, particularly during the dry season. Indeed, in this analysis, the number of malaria cases in 2003 and 2004 fell sharply when temperatures were particularly high.
A better understanding of the ecology of the anopheles mosquitoes endemic to the environment should make it possible to concentrate human and material resources more effectively in vector control programmes. In the absence of entomological and behavioural data, analyses of the incidence of malaria over time has provided information about the transmission of P. vivax in the region, and this information may be useful for the prevention and control of future epidemics.
Tangad MC, Ngundu WI, Judith N, Mbuh Tendongfor N, Simard F, Wanji S: Climate change and altitudinal structuring of malaria vectors in south-western Cameroon: their relation to malaria transmission. Trans R Soc Trop Med Hyg. 2010, 104: 453-460. 10.1016/j.trstmh.2010.02.006.
Sinka ME, Rubio-Palis Y, Manguin S, Patil AP, Temperley WH, Gething PW, Van Boeckel T, Kabaria CW, Harbach RE, Hay SI: The dominant Anopheles vectors of human malaria in the Americas: occurrence data, distribution maps and bionomic précis. Parasite Vectors. 2010, 3: 72-10.1186/1756-3305-3-72.
Girod R, Gaborit P, Issaly J, Carinci R, Basurko C: Anophèles et transmission du paludisme dans le village et les environs de Cacao (Guyane). Bull Soc Pathol Exot. 2009, 102: 199-205.
Dusfour I, Carinci R, Gaborit P, Issaly J, Girod R: Evaluation of four methods for collecting malaria vectors in French Guiana. J Econ Entomol. 2010, 103: 973-976. 10.1603/EC09328.
Hanf M, Stéphani A, Basurko C, Nacher M, Carme B: Determination of the Plasmodium vivax relapse pattern in Camopi French Guiana. Malar J. 2009, 8: 278-10.1186/1475-2875-8-278.
Tian L, Bi Y, Ho S, Liu W, Liang S, Gobbins WB, Chan EY, Zhoo S, Sung JJ: One-year delayed effect of fog on malaria transmission: a time-series analysis in the rain forest area of Mengla County, south-west China. Malar J. 2008, 7: 110-10.1186/1475-2875-7-110.
Roper M, Torres R, Goicochea C, Andersen E, Guarda J, Calampa C, Hightower A, Magill A: The epidemiology of malaria in an epidemic area of the Peruvian Amazon. Am J Trop Med Hyg. 2000, 62: 247-
Zeilhofer P, Santos E, Ribeiro A, Miyazaki R, Santos M: Habitat suitability mapping of Anopheles darlingi in the surroundings of the Manso hydropower plant reservoir, Mato Grosso, Central Brazil. Int J Health Geogr. 2007, 6: 7-10.1186/1476-072X-6-7.
Charlwood JD: Biological variation in Anopheles darlingi Root. Memórias do Instituto Oswaldo Cruz. 1996, 91: 391-398.
Giglioli G: Nation-wide malaria eradication projects in the Americas. III Eradication of Anopheles darlingi from the inhabited areas of British Guiana by DDT residual spraying. J Nat Mal Soc. 1951, 10: 142-161.
Vittor AY, Pan W, Gilman RH, Tielsch J, Glass G, Shields T, Sánchez-Lozano W, Pinedo VV, Salas-Cobos E, Flores S, Patz JA: Linking deforestation to malaria in the Amazon: characterization of the breeding habitat of the principal malaria vector, Anopheles darlingi. Am J Trop Med Hyg. 2009, 81: 5-12.
Grieco JP, Rejmánková E, Achee NL, Klein CN, Andre R, Roberts D: Habitat suitability for three species of Anopheles mosquitoes: larval growth and survival in reciprocal placement experiments. J Vector Ecol. 2007, 32: 176-87. 10.3376/1081-1710(2007)32[176:HSFTSO]2.0.CO;2.
Manguin S, Roberts DR, Andre RG, Rejmankova E, Hakre S: Characterization of Anopheles darlingi (Diptera: Culicidae)larval habitats in Belize, Central America. J Med Entomol. 1996, 33: 205-211.
Bautista C, Chan A, Ryan J, Calampa C, Roper M, Hightower A, Magill A: Epidemiology and spatial analysis of malaria in the northern Peruvian Amazon. Am J Trop Hyg. 2006, 75: 1216-1222.
Rozendaal J: Relations between Anopheles darlingi breeding habitats, rainfall, river level and malaria transmission rates in the rain forest of Suriname. Med Vet Entomol. 1992, 6: 16-22. 10.1111/j.1365-2915.1992.tb00029.x.
Barros F, Honório N: Man biting rate seasonal variation of malaria vectors in Roraima, Brazil. Memórias do Instituto Oswaldo Cruz. 2007, 102: 299-302.
Magris M, Rubio-Palis Y, Menares C, Villegas L: Vector bionomics and malaria transmission in the Upper Orinoco River, Southern Venezuela. Mem Inst Oswaldo Cruz. 2007, 102: 303-311. 10.1590/S0074-02762007005000049.
Bouma MJ, Dye C: Cycles of malaria associated with El Niño in Venezuela. JAMA. 1997, 278: 1772-1774. 10.1001/jama.278.21.1772.
Gomez-Elipe A, Otero A, Van Herp M, Aguirre-Jaime A: Forecasting malaria incidence based on monthly case reports and environmental factors in Karuzi, Burundi, 1997 - 2003. Malar J. 2007, 6: 129-10.1186/1475-2875-6-129.
Haghdoost A, Alexander N, Cox J: Modelling of malaria temporal variations in Iran. Trop Med Int Health. 2008, 13: 1501-1508. 10.1111/j.1365-3156.2008.02166.x.
Mabaso M, Vounatsou P, Midzi S, Da Silva J, Smith T: Spatio-temporal analysis of the role of climate in inter-annual variation of malaria incidence in Zimbabwe. Int J Health Geogr. 2006, 5: 20-10.1186/1476-072X-5-20.
Clements A, Barnett A, Cheng Z, Snow R, Zhou H: Space-time variation of malaria incidence in Yunnan province, China. Malar J. 2009, 8: 180-10.1186/1475-2875-8-180.
Kristan M, Abeku TA, Beard J, Okia M, Rapuoda B, Sang J, Cox J: Variations in entomological indices in relation to weather patterns and malaria incidence in East African highlands: implications for epidemic prevention and control. Malar J. 2008, 7: 231-10.1186/1475-2875-7-231.
We thank the staff of the Health Centre of Cacao and the staff of CIC-EC for their assistance.
The authors declare that they have no competing interests.
CB participated in the research design, performed data analysis and interpretation, and prepared the manuscript. MH participated in the statistical analysis. RHS, FG, CG and MJ participated in interpretation of data and manuscript revision as well. SG participated in research method. MN provided guidance on data analyses and was involved in the interpretation of data and manuscript revision. BC was responsible for data collection, initiated the study, and was involved in the interpretation of data and manuscript revision. All authors read and approved the final manuscript.