With limited resources and spatio-temporal heterogeneity of malaria in developing countries, it is still difficult to assess the real impact of socioeconomic and environmental factors in order to set up targeted campaigns against malaria at an accurate scale. Our goal was to detect malaria hotspots in rural area and assess the extent to which household socioeconomic status and meteorological recordings may explain the occurrence and evolution of these hotspots. Data on malaria cases from 2010 to 2014 and on socioeconomic and meteorological factors were acquired from four health facilities within the Nanoro demographic surveillance area. Statistical cross correlation was used to quantify the temporal association between weekly malaria incidence and meteorological factors. Local spatial autocorrelation analysis was performed and restricted to each transmission period using Kulldorff’s elliptic spatial scan statistic. Univariate and multivariable analysis were used to assess the principal socioeconomic and meteorological determinants of malaria hotspots using a Generalized Estimating Equation (GEE) approach. Rainfall and temperature were positively and significantly associated with malaria incidence, with a lag time of 9 and 14 weeks, respectively. Spatial analysis showed a spatial autocorrelation of malaria incidence and significant hotspots which was relatively stable throughout the study period. Furthermore, low socioeconomic status households were strongly associated with malaria hotspots (aOR = 1.21, 95% confidence interval: 1.03–1.40). These fine-scale findings highlight a relatively stable spatio-temporal pattern of malaria risk and indicate that social and environmental factors play an important role in malaria incidence. Integrating data on these factors into existing malaria struggle tools would help in the development of sustainable bottleneck strategies adapted to the local context for malaria control.

Background: According to the World Health Organization, there were more than 228 million cases of malaria globally in 2018, with 93% of cases occurring in Africa; in Mali, a 13% increase in the number of cases was observed between 2015 and 2018; this study aimed to evaluate the impact of meteorological and environmental factors on the geo-epidemiology of malaria in the health district of Dire, Mali. Methods: Meteorological and environmental variables were synthesized using principal component analysis and multiple correspondence analysis, the relationship between malaria incidence and synthetic indicators was determined using a multivariate general additive model; hotspots were detected by SaTScan. Results: Malaria incidence showed high inter and intra-annual variability; the period of high transmission lasted from September to February; health areas characterized by proximity to the river, propensity for flooding and high agricultural yield were the most at risk, with an incidence rate ratio of 2.21 with confidence intervals (95% CI: 1.85–2.58); malaria incidence in Dire declined from 120 to 20 cases per 10,000 person-weeks between 2013 and 2017. Conclusion: The identification of areas and periods of high transmission can help improve malaria control strategies.

Background: In malaria endemic areas, identifying spatio-temporal hotspots is becoming an important element of innovative control strategies targeting transmission bottlenecks. The aim of this work was to describe the spatio-temporal variation of malaria hotspots in central Senegal and to identify the meteorological, environmental, and preventive factors that influence this variation. Methods: This study analysed the weekly incidence of malaria cases recorded from 2008 to 2012 in 575 villages of central Senegal (total population approximately 500,000) as part of a trial of seasonal malaria chemoprevention (SMC). Data on weekly rainfall and annual vegetation types were obtained for each village through remote sensing. The time series of weekly malaria incidence for the entire study area was divided into periods of high and low transmission using change-point analysis. Malaria hotspots were detected during each transmission period with the SaTScan method. The effects of rainfall, vegetation type, and SMC intervention on the spatio-temporal variation of malaria hotspots were assessed using a General Additive Mixed Model. Results: The malaria incidence for the entire area varied between 0 and 115.34 cases/100,000 person weeks during the study period. During high transmission periods, the cumulative malaria incidence rate varied between 7.53 and 38.1 cases/100,000 person-weeks, and the number of hotspot villages varied between 62 and 147. During low transmission periods, the cumulative malaria incidence rate varied between 0.83 and 2.73 cases/100,000 person-weeks, and the number of hotspot villages varied between 10 and 43. Villages with SMC were less likely to be hotspots (OR = 0.48, IC95%: 0.33–0.68). The association between rainfall and hotspot status was non-linear and depended on both vegetation type and amount of rainfall. The association between village location in the study area and hotspot status was also shown. Conclusion: In our study, malaria hotspots varied over space and time according to a combination of meteorological, environmental, and preventive factors. By taking into consideration the environmental and meteorological characteristics common to all hotspots, monitoring of these factors could lead targeted public health interventions at the local level. Moreover, spatial hotspots and foci of malaria persisting during LTPs need to be further addressed. Trial registration: The data used in this work were obtained from a clinical trial registered on July 10, 2008 at www.clinicaltrials.gov under NCT00712374.

We introduce an approach based on functional data analysis to identify patterns of malaria incidence to guide effective targeting of malaria control in a seasonal transmission area. Using functional data method, a smooth function (functional data or curve) was fitted from the time series of observed malaria incidence for each of 575 villages in west-central Senegal from 2008 to 2012. These 575 smooth functions were classified using hierarchical clustering (Ward’s method), and several different dissimilarity measures. Validity indices were used to determine the number of distinct temporal patterns of malaria incidence. Epidemiological indicators characterizing the resulting malaria incidence patterns were determined from the velocity and acceleration of their incidences over time. We identified three distinct patterns of malaria incidence: high-, intermediate-, and low-incidence patterns in respectively 2% (12/575), 17% (97/575), and 81% (466/575) of villages. Epidemiological indicators characterizing the fluctuations in malaria incidence showed that seasonal outbreaks started later, and ended earlier, in the low-incidence pattern. Functional data analysis can be used to identify patterns of malaria incidence, by considering their temporal dynamics. Epidemiological indicators derived from their velocities and accelerations, may guide to target control measures according to patterns.

In the context of environmentally influenced communicable diseases, proximity to environmental sources results in spatial heterogeneity of risk, which is sometimes difficult to measure in the field. Most prevention trials use randomization to achieve comparability between groups, thus failing to account for heterogeneity. This study aimed to determine under what conditions spatial heterogeneity biases the results of randomized prevention trials, and to compare different approaches to modeling this heterogeneity. Using the example of a malaria prevention trial, simulations were performed to quantify the impact of spatial heterogeneity and to compare different models. Simulated scenarios combined variation in baseline risk, a continuous protective factor (age), a non-related factor (sex), and a binary protective factor (preventive treatment). Simulated spatial heterogeneity scenarios combined variation in breeding site density and effect, location, and population density. The performances of the following five statistical models were assessed: a non-spatial Cox Proportional Hazard (Cox-PH) model and four models accounting for spatial heterogeneity—i.e., a Data-Generating Model, a Generalized Additive Model (GAM), and two Stochastic Partial Differential Equation (SPDE) models, one modeling survival time and the other the number of events. Using a Bayesian approach, we estimated the SPDE models with an Integrated Nested Laplace Approximation algorithm. For each factor (age, sex, treatment), model performances were assessed by quantifying parameter estimation biases, mean square errors, confidence interval coverage rates (CRs), and significance rates. The four models were applied to data from a malaria transmission blocking vaccine candidate. The level of baseline risk did not affect our estimates. However, with a high breeding site density and a strong breeding site effect, the Cox-PH and GAM models underestimated the age and treatment effects (but not the sex effect) with a low CR. When population density was low, the Cox-SPDE model slightly overestimated the effect of related factors (age, treatment). The two SPDE models corrected the impact of spatial heterogeneity, thus providing the best estimates. Our results show that when spatial heterogeneity is important but not measured, randomization alone cannot achieve comparability between groups. In such cases, prevention trials should model spatial heterogeneity with an adapted method. The dataset used for the application example was extracted from Vaccine Trial #NCT02334462 ( ClinicalTrials.gov registry).