Impact of metric and sample size on determining malaria hotspot boundaries

The Potential of Surveillance Data for Dengue Risk Mapping: An Evaluation of Different Approaches in Cuba

To better guide dengue prevention and control efforts, the use of routinely collected data to develop risk maps is proposed. For this purpose, dengue experts identified indicators representative of entomological, epidemiological and demographic risks, hereafter called components, by using surveillance data aggregated at the level of Consejos Populares (CPs) in two municipalities of Cuba (Santiago de Cuba and Cienfuegos) in the period of 2010-2015. Two vulnerability models (one with equally weighted components and one with data-derived weights using Principal Component Analysis), and three incidence-based risk models were built to construct risk maps. The correlation between the two vulnerability models was high (tau > 0.89). The single-component and multicomponent incidence-based models were also highly correlated (tau ≥ 0.9). However, the agreement between the vulnerability- and the incidence-based risk maps was below 0.6 in the setting with a prolonged history of dengue transmission. This may suggest that an incidence-based approach does not fully reflect the complexity of vulnerability for future transmission. The small difference between single- and multicomponent incidence maps indicates that in a setting with a narrow availability of data, simpler models can be used. Nevertheless, the generalized linear mixed multicomponent model provides information of covariate-adjusted and spatially smoothed relative risks of disease transmission, which can be important for the prospective evaluation of an intervention strategy. In conclusion, caution is needed when interpreting risk maps, as the results vary depending on the importance given to the components involved in disease transmission. The multicomponent vulnerability mapping needs to be prospectively validated based on an intervention trial targeting high-risk areas.

Characterising spatial patterns of neglected tropical disease transmission using integrated sero-surveillance in Northern Ghana

Background

As prevalence decreases in pre-elimination settings, identifying the spatial distribution of remaining infections to target control measures becomes increasingly challenging. By measuring multiple antibody responses indicative of past exposure to different pathogens, integrated serological surveys enable simultaneous characterisation of residual transmission of multiple pathogens.

Methodology/Principal findings

Here, we combine integrated serological surveys with geostatistical modelling and remote sensing-derived environmental data to estimate the spatial distribution of exposure to multiple diseases in children in Northern Ghana. The study utilised the trachoma surveillance survey platform (cross-sectional two-stage cluster-sampled surveys) to collect information on additional identified diseases at different stages of elimination with minimal additional cost. Geostatistical modelling of serological data allowed identification of areas with high probabilities of recent exposure to diseases of interest, including areas previously unknown to control programmes. We additionally demonstrate how serological surveys can be used to identify areas with exposure to multiple diseases and to prioritise areas with high uncertainty for future surveys. Modelled estimates of cluster-level prevalence were strongly correlated with more operationally feasible metrics of antibody responses.

Conclusions/Significance

This study demonstrates the potential of integrated serological surveillance to characterise spatial distributions of exposure to multiple pathogens in low transmission and elimination settings when the probability of detecting infections is low.

A threshold-free approach with age-dependency for estimating malaria seroprevalence

Background

In malaria serology analysis, the standard approach to obtain seroprevalence, i.e the proportion of seropositive individuals in a population, is based on a threshold which is used to classify individuals as seropositive or seronegative. The choice of this threshold is often arbitrary and is based on methods that ignore the age-dependency of the antibody distribution.

Methods

Using cross-sectional antibody data from the Western Kenyan Highlands, this paper introduces a novel approach that has three main advantages over the current threshold-based approach: it avoids the use of thresholds; it accounts for the age dependency of malaria antibodies; and it allows us to propagate the uncertainty from the classification of individuals into seropositive and seronegative when estimating seroprevalence. The reversible catalytic model is used as an example for illustrating how to propagate this uncertainty into the parameter estimates of the model.

Results

This paper finds that accounting for age-dependency leads to a better fit to the data than the standard approach which uses a single threshold across all ages. Additionally, the paper also finds that the proposed threshold-free approach is more robust against the selection of different age-groups when estimating seroprevalence.

Conclusion

The novel threshold-free approach presented in this paper provides a statistically principled and more objective approach to estimating malaria seroprevalence. The introduced statistical framework also provides a means to compare results across studies which may use different age ranges for the estimation of seroprevalence.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-04022-4.

Spatial cluster analysis of Plasmodium vivax and P. malariae exposure using serological data among Haitian school children sampled between 2014 and 2016

BACKGROUND

Estimation of malaria prevalence in very low transmission settings is difficult by even the most advanced diagnostic tests. Antibodies against malaria antigens provide an indicator of active or past exposure to these parasites. The prominent malaria species within Haiti is Plasmodium falciparum, but P. vivax and P. malariae infections are also known to be endemic.

METHODOLOGY/PRINCIPAL FINDINGS

From 2014-2016, 28,681 Haitian children were enrolled in school-based serosurveys and were asked to provide a blood sample for detection of antibodies against multiple infectious diseases. IgG against the P. falciparum, P. vivax, and P. malariae merozoite surface protein 19kD subunit (MSP119) antigens was detected by a multiplex bead assay (MBA). A subset of samples was also tested for Plasmodium DNA by PCR assays, and for Plasmodium antigens by a multiplex antigen detection assay. Geospatial clustering of high seroprevalence areas for P. vivax and P. malariae antigens was assessed by both Ripley's K-function and Kulldorff's spatial scan statistic. Of 21,719 children enrolled in 680 schools in Haiti who provided samples to assay for IgG against PmMSP119, 278 (1.27%) were seropositive. Of 24,559 children enrolled in 788 schools providing samples for PvMSP119 serology, 113 (0.46%) were seropositive. Two significant clusters of seropositivity were identified throughout the country for P. malariae exposure, and two identified for P. vivax. No samples were found to be positive for Plasmodium DNA or antigens.

CONCLUSIONS/SIGNIFICANCE

From school-based surveys conducted from 2014 to 2016, very few Haitian children had evidence of exposure to P. vivax or P. malariae, with no children testing positive for active infection. Spatial scan statistics identified non-overlapping areas of the country with higher seroprevalence for these two malarias. Serological data provides useful information of exposure to very low endemic malaria species in a population that is unlikely to present to clinics with symptomatic infections.

Identifying Plasmodium falciparum transmission patterns through parasite prevalence and entomological inoculation rate

Background

Monitoring malaria transmission is a critical component of efforts to achieve targets for elimination and eradication. Two commonly monitored metrics of transmission intensity are parasite prevalence (PR) and the entomological inoculation rate (EIR). Comparing the spatial and temporal variations in the PR and EIR of a given geographical region and modelling the relationship between the two metrics may provide a fuller picture of the malaria epidemiology of the region to inform control activities.

Methods

Using geostatistical methods, we compare the spatial and temporal patterns of Plasmodium falciparum EIR and PR using data collected over 38 months in a rural area of Malawi. We then quantify the relationship between EIR and PR by using empirical and mechanistic statistical models.

Results

Hotspots identified through the EIR and PR partly overlapped during high transmission seasons but not during low transmission seasons. The estimated relationship showed a 1-month delayed effect of EIR on PR such that at lower levels of EIR, increases in EIR are associated with rapid rise in PR, whereas at higher levels of EIR, changes in EIR do not translate into notable changes in PR.

Conclusions

Our study emphasises the need for integrated malaria control strategies that combine vector and human host managements monitored by both entomological and parasitaemia indices.

Funding

This work was supported by Stichting Dioraphte grant number 13050800.

A unified and flexible modelling framework for the analysis of malaria serology data

Serology data are an increasingly important tool in malaria surveillance, especially in low transmission settings where the estimation of parasite-based indicators is often problematic. Existing methods rely on the use of thresholds to identify seropositive individuals and estimate transmission intensity, while making assumptions about the temporal dynamics of malaria transmission that are rarely questioned. Here, we present a novel threshold-free approach for the analysis of malaria serology data which avoids dichotomization of continuous antibody measurements and allows us to model changes in the antibody distribution across age in a more flexible way. The proposed unified mechanistic model combines the properties of reversible catalytic and antibody acquisition models, and allows for temporally varying boosting and seroconversion rates. Additionally, as an alternative to the unified mechanistic model, we also propose an empirical approach to analysis where modelling of the age-dependency is informed by the data rather than biological assumptions. Using serology data from Western Kenya, we demonstrate both the usefulness and limitations of the novel modelling framework.

A Model for Assessing the Quantitative Effects of Heterogeneous Affinity in Malaria Transmission along with Ivermectin Mass Administration

Using an agent-based model of malaria, we present numerical evidence that in communities of individuals having an affinity varying within a broad range of values, disease transmission may increase up to 300%. Moreover, our findings provide new insight into how to combine different strategies for the prevention of malaria transmission. In particular, we uncover a relationship between the level of heterogeneity and the level of conventional and unconventional anti-malarial drug administration (ivermectin and gametocidal agents), which, when taken together, will define a control parameter, tuning between disease persistence and elimination. Finally, we also provide evidence that the entomological inoculation rate, as well as the product between parasite and sporozoite rates are both good indicators of malaria incidence in the presence of heterogeneity in disease transmission and may configure a possible improvement in that setting, upon classical standard measures such as the basic reproductive number.

Discovering, Defining, and Summarizing Persistent Hotspots in SCORE Studies

The Schistosomiasis Consortium for Operational Research and Evaluation (SCORE) conducted large field studies on schistosomiasis control and elimination in Africa. All of these studies, carried out in low-, moderate-, and high-prevalence areas, resulted in a reduction in prevalence and intensity of Schistosoma infection after repeated mass drug administration (MDA). However, in all studies, there were locations that experienced minimal or no decline or even increased in prevalence and/or intensity. These areas are termed persistent hotspots (PHS). In SCORE studies in medium- to high-prevalence areas, at least 30% of study villages were PHS. There was no consistent relationship between PHS and the type or frequency of intervention, adequacy of reported MDA coverage, and prevalence or intensity of infection at baseline. In a series of small studies, factors that differed between PHS and villages that responded to repeated MDA as expected included sources of water for personal use, sanitation, and hygiene. SCORE studies comparing PHS with villages that responded to MDA suggest the potential for PHS to be identified after a few years of MDA. However, additional studies in different social-ecological settings are needed to develop generalizable approaches that program managers can use to identify and address PHS. This is essential if goals for schistosomiasis control and elimination are to be achieved.

Using health facility-based serological surveillance to predict receptive areas at risk of malaria outbreaks in elimination areas

BACKGROUND

In order to improve malaria burden estimates in low transmission settings, more sensitive tools and efficient sampling strategies are required. This study evaluated the use of serological measures from repeated health facility-based cross-sectional surveys to investigate Plasmodium falciparum and Plasmodium vivax transmission dynamics in an area nearing elimination in Indonesia.

METHODS

Quarterly surveys were conducted in eight public health facilities in Kulon Progo District, Indonesia, from May 2017 to April 2018. Demographic data were collected from all clinic patients and their companions, with household coordinates collected using participatory mapping methods. In addition to standard microscopy tests, bead-based serological assays were performed on finger-prick bloodspot samples from 9453 people. Seroconversion rates (SCR, i.e. the proportion of people in the population who are expected to seroconvert per year) were estimated by fitting a simple reversible catalytic model to seroprevalence data. Mixed effects logistic regression was used to examine factors associated with malaria exposure, and spatial analysis was performed to identify areas with clustering of high antibody responses.

RESULTS

Parasite prevalence by microscopy was extremely low (0.06% (95% confidence interval 0.03-0.14, n = 6) and 0 for P. vivax and P. falciparum, respectively). However, spatial analysis of P. vivax antibody responses identified high-risk areas that were subsequently the site of a P. vivax outbreak in August 2017 (62 cases detected through passive and reactive detection systems). These areas overlapped with P. falciparum high-risk areas and were detected in each survey. General low transmission was confirmed by the SCR estimated from a pool of the four surveys in people aged 15 years old and under (0.020 (95% confidence interval 0.017-0.024) and 0.005 (95% confidence interval 0.003-0.008) for P. vivax and P. falciparum, respectively). The SCR estimates in those over 15 years old were 0.066 (95% confidence interval 0.041-0.105) and 0.032 (95% confidence interval 0.015-0.069) for P. vivax and P. falciparum, respectively.

CONCLUSIONS

These findings demonstrate the potential use of health facility-based serological surveillance to better identify and target areas still receptive to malaria in an elimination setting. Further implementation research is needed to enable integration of these methods with existing surveillance systems.

Malaria Hotspots: Is There Epidemiological Evidence for Fine-Scale Spatial Targeting of Interventions?

As data at progressively granular spatial scales become available, the temptation is to target interventions to areas with higher malaria transmission - so-called hotspots - with the aim of reducing transmission in the wider community. This paper reviews literature to determine if hotspots are an intrinsic feature of malaria epidemiology and whether current evidence supports hotspot-targeted interventions. Hotspots are a consistent feature of malaria transmission at all endemicities. The smallest spatial unit capable of supporting transmission is the household, where peri-domestic transmission occurs. Whilst the value of focusing interventions to high-burden areas is evident, there is currently limited evidence that local-scale hotspots fuel transmission. As boundaries are often uncertain, there is no conclusive evidence that hotspot-targeted interventions accelerate malaria elimination.

Heterogeneous exposure and hotspots for malaria vectors at three study sites in Uganda

Background: Heterogeneity in malaria transmission has household, temporal, and spatial components. These factors are relevant for improving the efficiency of malaria control by targeting heterogeneity. To quantify variation, we analyzed mosquito counts from entomological surveillance conducted at three study sites in Uganda that varied in malaria transmission intensity. Mosquito biting or exposure is a risk factor for malaria transmission.

Methods: Using a Bayesian zero-inflated negative binomial model, validated via a comprehensive simulation study, we quantified household differences in malaria vector density and examined its spatial distribution. We introduced a novel approach for identifying changes in vector abundance hotspots over time by computing the Getis-Ord statistic on ratios of household biting propensities for different scenarios. We also explored the association of household biting propensities with housing and environmental covariates.

Results: In each site, there was evidence for hot and cold spots of vector abundance, and spatial patterns associated with urbanicity, elevation, or other environmental covariates. We found some differences in the hotspots in rainy vs. dry seasons or before vs. after the application of control interventions. Housing quality explained a portion of the variation among households in mosquito counts.

Conclusion: This work provided an improved understanding of heterogeneity in malaria vector density at the three study sites in Uganda and offered a valuable opportunity for assessing whether interventions could be spatially targeted to be aimed at abundance hotspots which may increase malaria risk. Indoor residual spraying was shown to be a successful measure of vector control interventions in Tororo, Uganda.  Cement walls, brick floors, closed eaves, screened airbricks, and tiled roofs were features of a house that had shown reduction of household biting propensity. Improvements in house quality should be recommended as a supplementary measure for malaria control reducing risk of infection.

Do hotspots fuel malaria transmission: a village-scale spatio-temporal analysis of a 2-year cohort study in The Gambia

BACKGROUND

Despite the biological plausibility of hotspots fueling malaria transmission, the evidence to support this concept has been mixed. If transmission spreads from high burden to low burden households in a consistent manner, then this could have important implications for control and elimination program development.

METHODS

Data from a longitudinal cohort in The Gambia was analyzed. All consenting individuals residing in 12 villages across the country were sampled monthly from June (dry season) to December 2013 (wet season), in April 2014 (mid dry season), and monthly from June to December 2014. A study nurse stationed within each village recorded passively detected malaria episodes between visits. Plasmodium falciparum infections were determined by polymerase chain reaction and analyzed using a geostatistical model.

RESULTS

Household-level observed monthly incidence ranged from 0 to 0.50 infection per person (interquartile range = 0.02-0.10) across the sampling months, and high burden households exist across all study villages. There was limited evidence of a spatio-temporal pattern at the monthly timescale irrespective of transmission intensity. Within-household transmission was the most plausible hypothesis examined to explain the observed heterogeneity in infections.

CONCLUSIONS

Within-village malaria transmission patterns are concentrated in a small proportion of high burden households, but patterns are stochastic regardless of endemicity. Our findings support the notion of transmission occurring at the household and village scales but not the use of a targeted approach to interrupt spreading of infections from high to low burden areas within villages in this setting.

Micro-epidemiology of mixed-species malaria infections in a rural population living in the Colombian Amazon region

Malaria outbreaks have been reported in recent years in the Colombian Amazon region, malaria has been re-emerging in areas where it was previously controlled. Information from malaria transmission networks and knowledge about the population characteristics influencing the dispersal of parasite species is limited. This study aimed to determine the distribution patterns of Plasmodium vivax, P. malariae and P. falciparum single and mixed infections, as well as the significant socio-spatial groupings relating to the appearance of such infections. An active search in 57 localities resulted in 2,106 symptomatic patients being enrolled. Parasitaemia levels were assessed by optical microscopy, and parasites were detected by PCR. The association between mixed infections (in 43.2% of the population) and socio-spatial factors was modelled using logistic regression and multiple correspondence analyses. P. vivax occurred most frequently (71.0%), followed by P. malariae (43.2%), in all localities. The results suggest that a parasite density-dependent regulation model (with fever playing a central role) was appropriate for modelling the frequency of mixed species infections in this population. This study highlights the under-reporting of Plasmodium spp. mixed infections in the malaria-endemic area of the Colombian Amazon region and the association between causative and environmental factors in such areas.

Using non-exceedance probabilities of policy-relevant malaria prevalence thresholds to identify areas of low transmission in Somalia

Background

Countries planning malaria elimination must adapt from sustaining universal control to targeted intervention and surveillance. Decisions to make this transition require interpretable information, including malaria parasite survey data. As transmission declines, observed parasite prevalence becomes highly heterogeneous with most communities reporting estimates close to zero. Absolute estimates of prevalence become hard to interpret as a measure of transmission intensity and suitable statistical methods are required to handle uncertainty of area-wide predictions that are programmatically relevant.

Methods

A spatio-temporal geostatistical binomial model for Plasmodium falciparum prevalence (PfPR) was developed using data from cross-sectional surveys conducted in Somalia in 2005, 2007–2011 and 2014. The fitted model was then used to generate maps of non-exceedance probabilities, i.e. the predictive probability that the region-wide population-weighted average PfPR for children between 2 and 10 years (PfPR_2–10) lies below 1 and 5%. A comparison was carried out with the decision-making outcomes from those of standard approaches that ignore uncertainty in prevalence estimates.

Results

By 2010, most regions in Somalia were at least 70% likely to be below 5% PfPR_2–10 and, by 2014, 17 regions were below 5% PfPR_2–10 with a probability greater than 90%. Larger uncertainty is observed using a threshold of 1%. By 2011, only two regions were more than 90% likely of being < 1% PfPR_2–10 and, by 2014, only three regions showed such low level of uncertainty. The use of non-exceedance probabilities indicated that there was weak evidence to classify 10 out of the 18 regions as < 1% in 2014, when a greater than 90% non-exceedance probability was required.

Conclusion

Unlike standard approaches, non-exceedance probabilities of spatially modelled PfPR_2–10 allow to quantify uncertainty of prevalence estimates in relation to policy relevant intervention thresholds, providing programmatically relevant metrics to make decisions on transitioning from sustained malaria control to strategies that encompass methods of malaria elimination.

Electronic supplementary material

The online version of this article (10.1186/s12936-018-2238-0) contains supplementary material, which is available to authorized users.