- Research article
- Open access
- Published:
Geospatial analysis of Plasmodium falciparum serological indicators: school versus community sampling in a low-transmission malaria setting
BMC Medicine volume 22, Article number: 31 (2024)
Abstract
Background
Due to low numbers of active infections and persons presenting to health facilities for malaria treatment, case-based surveillance is inefficient for understanding the remaining disease burden in low malaria transmission settings. Serological data through the detection of IgG antibodies from previous malaria parasite exposure can fill this gap by providing a nuanced picture of where sustained transmission remains. Study enrollment at sites of gathering provides a potential approach to spatially estimate malaria exposure and could preclude the need for more intensive community-based sampling.
Methods
This study compared spatial estimates of malaria exposure from cross-sectional school- and community-based sampling in Haiti. A total of 52,405 blood samples were collected from 2012 to 2017. Multiplex bead assays (MBAs) tested IgG against P. falciparum liver stage antigen-1 (LSA-1), apical membrane antigen 1 (AMA1), and merozoite surface protein 1 (MSP1). Predictive geospatial models of seropositivity adjusted for environmental covariates, and results were compared using correlations by coordinate points and communes across Haiti.
Results
Consistent directional associations were observed between seroprevalence and environmental covariates for elevation (negative), air temperature (negative), and travel time to urban centers (positive). Spearman’s rank correlation for predicted seroprevalence at coordinate points was lowest for LSA-1 (ρ = 0.10, 95% CI: 0.09–0.11), but improved for AMA1 (ρ = 0.36, 95% CI: 0.35–0.37) and MSP1 (ρ = 0.48, 95% CI: 0.47–0.49).
Conclusions
In settings approaching P. falciparum elimination, case-based prevalence data does not provide a resolution of ongoing malaria transmission in the population. Immunogenic antigen targets (e.g., AMA1, MSP1) that give higher population rates of seropositivity provide moderate correlation to gold standard community sampling designs and are a feasible approach to discern foci of residual P. falciparum transmission in an area.
Background
Malaria is a disease caused by infection with Plasmodium parasites, and global cases and deaths have declined substantially since 2000 [1], with 25 countries achieving zero locally acquired cases and many more approaching this significant milestone [2, 3]. For countries approaching elimination, lower overall parasite burden in the population and symptomatic disease prevalence in the population increases difficulty in estimating the remaining areas of Plasmodium transmission [4]. P. falciparum incidence is often estimated through passive case detection from patients presenting to health facilities, but also is accomplished through active case detection in community settings that can capture the infected, but non-treatment–seeking population [5, 6]. Active case detection is typically conducted through diagnostic testing of high-risk individuals with microscopy or rapid diagnostic tests (RDTs), but low-density parasite infections below diagnostic test limit of detection, or history of past infections, cannot be captured through these methods [7]. Similarly, parasite prevalence estimates lack sensitivity and statistical or spatial precision when numbers of parasite-positive samples are low [8, 9]. Thus, in low-transmission malaria settings, more sensitive tools are required to understand a population’s infection burden and where transmission foci remain.
Assays detecting antibodies elicited in response to P. falciparum infection are a sensitive indicator of personal malaria infection history, and multiple antigen targets have been identified with varying durability of seropositivity in humans [10,11,12,13]. Antibodies produced during sporozoite and liver stages of the parasite life cycle (such as circumsporozoite protein (CSP) and liver stage antigen-1 (LSA-1), respectively) can be used to estimate an individual’s P. falciparum exposure within the past few months [14, 15]. Conversely, IgG generated against blood stage antigens (including apical membrane antigen 1 (AMA1) and merozoite surface protein 1 (MSP1)) can remain for years, or even decades, post-exposure [14, 15]. Serological data are more commonly analyzed at the population level in order to garner historical estimates of malaria exposure in a setting (or a sub-population) [8, 12, 16]. Recent advancements in serological assay technology have allowed measurement of IgG against multiple antigens simultaneously in a serum sample, and the multiplex bead assay (MBA) has become a popular tool for this objective [17, 18].
Serological data have been collected from various survey designs to estimate malaria exposure, including population-representative surveys [19, 20] as well as those conducted at venues where persons gather [21, 22]. The convenience and accessibility of venue-based sampling presents an attractive methodology for enrollment of individuals for malaria estimates [23], but an important question for investigation is how estimates from venue-based sampling would compare to those from surveys designed to representatively sample the population.
Haiti is one of two Caribbean countries endemic for malaria, and recent cross-sectional surveys have found low levels of nationwide transmission in the country with approximately 1% or less of persons with active P. falciparum infection [24,25,26]. Investments into the nation’s health system after the 2010 earthquake and subsequent cholera epidemic resulted in intensified malaria elimination efforts [26], and surveys at sites of gathering have complemented population-representative surveys in efforts to estimate remaining burden [27]. This study compares MBA IgG data collected from population-representative community surveys in Haiti to those from school-based surveys conducted from 2012 to 2017. Geostatistical predictive modeling yielded estimates of school-based seroprevalence for comparison with those from community-based surveys, which were considered the gold standard for representative population sampling.
Methods
TRaC and TAS study designs and ethics
Tracking Results Continuously (TRaC) were community-based surveys conducted in 2012, 2015, and 2017 as part of The Global Fund round 8 malaria surveillance activities in Haiti. Surveys were implemented by Population Services International (PSI)-Haiti to characterize the burden of malaria across Haiti. Signed informed consent forms were obtained from all participants; for children under 15 years, verbal assent was given with documented parental consent. Enumeration areas throughout the country (sections d'énumération, SDE) were chosen on a proportional sampling by predicted malaria risk strata as determined by predictive modeling. A target of 20 households was randomly selected by field teams within each SDE, and all members of the household were offered the opportunity to participate. Blood was collected by fingerprick on Whatman 903 Protein Saver cards (GE Healthcare, Chicago, IL), dried overnight, and individually stored in plastic bags with desiccant at − 20 °C until blood elution and IgG assay. RDTs (CareStart HRP2; AccessBio, Somerset, NJ) were provided to all participants, and individuals testing positive were treated according to national guidelines. Blood samples were assigned unique identification numbers that were not traceable to the individual. The study protocols were approved by the Haitian Ministry of Health and approved as a non-research activity by the Center for Global Health Human Research Protection Office (HRPO), US Centers for Disease Control and Prevention (CDC; Center for Global Health determination #2015–04).
Lymphatic filariasis transmission assessment surveys (TASs) were conducted in Haiti from 2014 to 2017 with integration of malaria RDTs (First Response Malaria Histidine-Rich Protein II (HRP2), II3FRC30, Premier Medical Corporation, New Jersey) into surveys in the years 2016 and 2017 [27]. As described previously [25, 27], TASs were conducted in evaluation units (EUs) that had met World Health Organization (WHO) criteria to conduct lymphatic filariasis school-based TASs (where the net primary-school enrollment rate is ≥ 75%), and target sample size was determined by WHO-provided tables [28]. For enrollment into a TAS, school headmasters were contacted in advance regarding the survey’s purpose and asked to notify parents; children’s verbal assent was provided before sample collection. Approximately 60 µL of fingerprick blood was collected on filter papers (TropBio filter wheels, Cellabs, Sydney, Australia), dried to create a dried blood spot (DBS), and packaged individually with desiccant for later laboratory analysis at the Centers for Disease Control and Prevention in Atlanta, GA. DBS were kept at − 20 °C until laboratory analysis. The study protocol was approved by the National Bioethics Committee of Haiti, and this activity was considered a program evaluation activity by the CDC Human subjects office (#2014–256). Persons consented to future laboratory testing of DBS, and CDC laboratory staff did not have access to any personal identifiers.
Plasmodium falciparum recombinant antigens
Purified recombinant antigens were covalently linked to xMAP microspheres (Luminex Corp, Austin, TX), as described previously [29]. Three antigens were utilized for all surveys: P. falciparum liver stage antigen 1 Pl1043 epitope peptide (LSA-1, coupled at pH 5 at 60 µg/mL) [30]; recombinant merozoite surface protein 1 19kD fragment (MSP1, coupled at pH 5 at 20 µg/mL) [31]; and recombinant apical membrane antigen 1 N-terminal region (AMA1, coupled at pH 5 at 20 µg/mL) [31]. These three P. falciparum antigens have been widely used for previous malaria serological studies by multiple groups, and represent the evaluation of long-term IgG responses (MSP1 19kD region, AMA1 N-terminal region) [8, 10, 14,15,16] as well as host IgG responses with shorter duration (LSA-1 Pl1043 epitope) [14, 15].
Sample preparation and multiplex IgG detection assay
A single six-mm DBS punch of the Whatman 903 card for TRaC surveys or a single DBS tab of the TropBio filter wheel for TAS surveys was rehydrated in blocking buffer: PBS pH 7.2; 0.5% polyvinyl alcohol (Sigma-Aldrich, St. Louis, MO); 0.5% polyvinylpyrrolidine (Sigma-Aldrich); 0.1% casein (ThermoFisher, Waltham, MA); 0.5% BSA (Sigma-Aldrich); 0.3% Tween-20; 0.05% sodium azide; and 3 µg/mL E. coli extract to prevent non-specific binding. For the multiplex IgG detection assay, samples were diluted to a final concentration of 1:200 whole blood, which is approximately a 1:400 serum dilution with the assumption of 50% hematocrit in whole blood.
IgG detection through MBA was performed as described previously [32]. Briefly, samples were incubated with a mixture of all three bead regions in reagent buffer (PBS, 0.05% Tween-20, 0.5% BSA, 0.02% NaN3) for 90 min at room temperature under gentle shaking protected from light in MultiScreen-BV filter plates (MilliporeSigma, Burlington, MA) or Bio-Plex Pro plates (Bio-Rad, Hercules, CA). After three washes (wash buffer: PBS, 0.05% Tween-20), beads were incubated with 50 µL biotinylated detection antibody (a mixture of 1:500 anti-hIgG and 1:625 anti-hIgG4, both produced by Southern Biotech, Birmingham, AL) for 45 min with the same incubation conditions as above. After three washes, 50 µL of a 1:250 × dilution of streptavidin–phycoerythrin (Invitrogen, Waltham, MA) was added to all wells for a 30-min incubation. After three washes, sample beads were incubated with 50 µL reagent buffer for 30 min to remove loosely-bound IgG, washed once, and resuspended in 100 µL PBS. Assay plates were briefly shaken and read on a Bio-Plex 200 machine (Bio-Rad) at a high RP1 target or on a MAGPIX machine (Luminex Corp) by generating the median fluorescence intensity (MFI) for 50 beads. The final measure, denoted as MFI-bg, was reported by subtracting MFI values from beads on each plate only exposed to sample diluent during the sample incubation step. Malaria IgG positive and negative samples were included in duplicate on each assay plate to ensure appropriate plate reading. The MFI-bg threshold for true positive IgG assay signal for each survey was ascertained if the sample MFI-bg was higher than the mean + 3 SD of the MFI-bg signal of a panel of known IgG negative DBS samples.
Statistical analysis
Data analysis was performed in SAS (version 9.4; SAS Institute Inc., Cary, USA). For comparison purposes, covariates of interest included in the analysis were only those that had been collected for both TAS and TRaC studies. For the TRaC study, age was classified into approximately evenly distributed groups. For TASs, only children ages 6 and 7 were included. To assess the assumption that TRaC seroprevalence was consistent across the 2012–2017 study period, a reversible catalytic model was fit to seropositivity by age data for AMA1 and MSP1 using 15 age categories. Estimates for serological conversion rate (SCR) and serological reversion rate (SRR) for each survey year were directly derived from the likelihood model [16].
Geostatistical analysis of seroprevalence to each antigen was conducted using the R-INLA package [33] in R Statistical Software (version 4.2.0; R Foundation for Statistical Computing, Vienna, Austria). Covariates have known or putative effects on malaria exposure (Additional file 1: Tab S1). Data values for each covariate were sampled at each TAS/TRaC survey location’s GPS coordinates using QGIS (version 3.20.3-Odense). For temporal covariates, annual averages from the latest year of data available were used. All covariates were standardized and centered to facilitate model convergence, and none were omitted from the models to allow for direct comparison between studies and to evaluate all effects on malaria exposure.
Geostatistical models of seroprevalence to LSA-1, AMA1, and MSP1 for each study.
were fit using a Bayesian framework, whereby P(xi) was seroprevalence at survey locations xi, i = 1, …, n, and the number of seropositive participants, Yi, out of the number of Ni participants surveyed at a location was assumed to follow a binomial distribution:
The seroprevalence at each location was then linked to its linear predictors through the logit:
where, for each model, β0 was the intercept, \(d{\left({x}_{i}\right)}^{\prime}\beta\) represented the vector of location-specific covariate effects, and \({S(x}_{i})\) was the spatial random effect, modeled using a Matérn covariance function with a stochastic partial difference equation (SPDE) (Additional file 2: Fig. S1) in integrated nested Laplace approximation (INLA). Models were validated splitting datasets into 75% training values and 25% test values. Results were extracted as mean estimated beta coefficients and as continuous surfaces of seroprevalence with a pixel size of 0.9 km2. The latter were imported into QGIS for sampling of predicted values for all coordinates in Haiti (n = 32,735) and estimating mean seroprevalence within communes (n = 134). Spearman’s rank correlation coefficients compared predictions at points and communes.
Results
From the 2014–2017 Haiti TASs, 36,094 children ages 6–7 years from 1163 schools provided DBS for IgG detection (Fig. 1, Table 1). Through Haiti TRaC surveys in 2012, 2015, and 2017, 16,311 participants of all ages enrolled from 324 communities (SDEs) (Fig. 1, Table 1). Of 16,063 children with P. falciparum RDT results from TASs, 18 (0.1%, 95% CI: 0.0–0.1%) tested positive; of 15,754 persons from TRaC with RDT results, 99 (0.6%, 0.5–0.8%) tested positive. Due to these low numbers, no clear spatial pattern was observed in RDT positives from both surveys (Additional file 2: Fig. S2). For the years when tests were performed, RDT prevalence was low (≤ 1% positivity) across time for both survey types (Additional file 1: Tab S2).
IgG seropositivity to the liver-stage antigen LSA-1 was the lowest among all antigens included in the multiplex assay at 1.5% (95% CI: 1.3–1.7%) for TRaC and 0.4% (0.3–0.4%) for TAS. Seropositivity to AMA1 was 15.0% (14.5–15.6%) for TRaC, and 3.8% (3.6–4.0%) for TAS; similarly, seropositivity to MSP1 was 17.8% (17.2–18.4%) for TRaC and 4.5% (4.3–4.7%) for TAS. Seroprevalence curves and seropositivity estimates were consistent among TRaC survey years (Additional file 2: Fig. S3, Additional file 1: Tab S3). Most TAS schools (91%) and TRaC SDEs (60%) had no participants seropositive to LSA-1, whereas at least 44% of schools and SDEs had persons seropositive to AMA1 or MSP1 (Fig. 2). TRaC sites had higher median population density, sparser vegetation, and lower travel time to urban centers (Additional files 1: Tab S4 and S5). Elevation, temperature, rainfall, and distance to the nearest river, lake, or stream were generally comparable between studies. Individual-level correlations in the IgG level among LSA-1, MSP1, and AMA1 antigens are shown in Additional file 2: Fig. S4 and showed no marked correlation among IgG levels to these three P. falciparum targets.
From spatial predictive modeling, many covariates had similar relationships with seroprevalence outcomes between the two survey types (Fig. 3, Additional file 1: Tab S6). Increased air temperature was consistently associated with decreased seroprevalence to P. falciparum targets, while greater travel time to urban centers (i.e., a more rural setting) was positively associated with seroprevalence. There was a negative association between elevation and seroprevalence across both survey types. Increases in rainfall and distance to the nearest body of water were generally associated with increases in seroprevalence. Based on cross-validation results, the R-INLA model of seroprevalence to LSA-1 in the TAS venue-based survey had the lowest prediction ability (Additional file 2: Fig. S5).
Due to the enrollment of older ages in TRaC, observed and predicted seroprevalence was higher and more geographically widespread for TRaC (Fig. 4). Correlations were weakest for LSA-1 (ρpoints = 0.10, 95% CI: 0.09–0.11; ρcommunes = -0.03, 95% CI: − 0.20–0.14) and strongest for MSP1 (ρpoints = 0.48, 0.47–0.49; ρcommunes = 0.49, 0.34–0.61) when examining strength of correlations between predicted seroprevalence values (Additional file 2: Fig S6). The MSP1 target had the highest observed and predicted seroprevalence across both studies. The AMA1 target, with lower values relative to MSP1, had lower correlations (ρpoints = 0.36, 0.35–0.37; ρcommunes = 0.40, 0.25–0.53) (Table 2). While higher population seropositivity was observed if persons were considered “seropositive” when IgG positive to either AMA1 or MSP1 targets (23.5% for TRaC and 6.8% for TAS), there was no improvement in correlations of predicted estimates (ρpoints = 0.43, 0.42–0.43; ρcommunes = 0.46, 0.32–0.58).
Discussion
This study utilized geospatial modeling to predict seroprevalence to different P. falciparum antigens to compare two sampling methodologies: population-representative household-based enrollment (TRaC surveys), and school-based enrollment (TASs). Regions approaching P. falciparum elimination are faced with epidemiological challenges to accurately identify the remaining foci of transmission [2, 4, 7]. In low-transmission settings where persons presenting to healthcare settings with malaria infection is rare, seropositivity to P. falciparum antigens extends the window of a “positive” malaria exposure and can indicate where (and within whom) infections are still occurring [8, 34, 35]. As malaria transmission decreases in a population, resources become diverted away from malaria programs in favor of matters considered more important to public health. For this and other reasons, generating accurate serological estimates from venue-based sampling or survey designs implemented by non-malaria groups could increase efficiency and augment malaria epidemiology in these settings.
Among the years of sample collection, tests for active P. falciparum infection by RDT confirmed a low prevalence of P. falciparum in Haiti, and each cross-sectional survey found RDT positivity of 1.0% or less. The same paucity of infections through cross-sectional surveys with RDTs has been observed in other low-transmission settings [37, 38]. These studies have shown reductions in RDT diagnostic performance among non-treatment–seeking persons [36, 37], and demonstrated that RDTs cannot reliably detect asymptomatic infections accounting for a significant proportion of P. falciparum burden in these populations [38, 39]. These RDT results diverge with serological findings from this study, which identified persons with anti-P. falciparum IgG throughout Haiti and provided augmented numbers of positives to allow for statistical modeling. Much like RDT positivity, seropositivity to LSA-1, reflecting P. falciparum exposure in the past few months [15], was low across both survey types, reaching a maximum seropositivity of 1.5% for persons enrolled in TRaC. Cumulative exposure to malaria over time increases the proportion of persons seropositive to P. falciparum antigens in a population, and this is especially profound for IgG targets such as AMA1 and MSP1 that are known to provide “long-lived” B cell responses and seropositivity for years after an infection [8, 16, 40]. This was noted by the AMA1 and MSP1 seroprevalence curves for the TRaC surveys, which showed no change in exposure estimates from 2012 to 2017. Because TASs only included children 6 and 7 years of age, seroprevalence to all P. falciparum antigens was overall lower when compared to persons enrolled in TRaC surveys, which were inclusive of all ages.
Results from INLA-SPDE models showed that certain environmental variables yielded consistent relationships between survey types. Elevation had a negative association with all serological outcomes. Multiple studies have also observed this association with increased elevation and decreased parasite prevalence [34], malaria incidence [41], and seroprevalence to P. falciparum antigens [16, 19], as lower elevations allow for more suitable temperature and moisture conditions for the Anopheles vector to breed [42]. Increased travel time to urban centers also had a consistently positive relationship with seropositivity. Greater mean travel time to urban centers would indicate more rural survey locations, and rural residence has been consistently associated with greater Anopheline habitat and malaria transmission when compared to urban settings, as well as generally lower socioeconomic indicators for rural residents [43, 44]. In line with improved conditions for mosquito habitat, rainfall had a positive relationship with P. falciparum seropositivity for all models except for one. This positive relationship has also been noted in multiple malaria studies [45, 46]. Air temperature was found to have a consistently negative association with seropositivity. Previous studies have found positive and negative associations between air temperature and malaria incidence; given that temperature interacts with other factors influencing malaria prevalence, these conflicting findings are not unexpected [47, 48].
Spatial patterns predicted from this analysis are reflected in maps produced by Cameron et al., which show higher clinical incidence in central and southwestern Haiti compared to other areas [49]. A higher risk of transmission has also previously been reported in the southern coastal regions of Haiti [50]. Between the TRaC and TAS surveys, the AMA1 and MSP1 IgG targets had both the highest observed seropositivity and highest correlations when comparing predicted seropositivity. However, if defining “seropositivity” as IgG positive to either long-term IgG target (in the same manner as previous studies [21, 35]), while the proportion of those who were seropositive increased, a correlation between predictions was not improved. This may reflect a tradeoff between sensitivity and specificity of a true “seropositive” designation, whereas a higher number of positives may be compounded by more false positives through a multi-antigen approach. While the magnitude of the correlation between TRaC and TAS predicted seroprevalence was similar whether comparing points or communes, correlation coefficients for coordinate points had narrower confidence intervals than those for communes. This emphasizes the utility of high-resolution data points in geostatistical analysis over data aggregated to administrative boundaries. Moderate agreement in predictive spatial outputs between the two survey types shows some utility in using venue-based sampling methods as a proxy for sampling designed to be representative of an entire population.
Limitations to this study include limited sociological and individual covariates that were common between TRaC and TAS surveys that could be used for predictive modeling. While TRaC surveys were designed to be nationally representative with an exclusive focus on malaria, TAS surveys were powered based on estimated lymphatic filariasis antigenemia within enumeration units throughout Haiti. More comprehensive TAS sampling with consideration for malaria may have improved correlations between predictive outputs, in addition to geographical representation of the same study areas and inclusion of a wider age range for more extensive serological data analysis. A wider age range would also be ideal when using venue-based sampling methods as a proxy for population-level transmission or exposure estimates. This study used a limited panel of P. falciparum IgG targets, and dozens of potential targets are now available [41]. An expanded panel may find other antigen targets (or combinations thereof) with better correlations among survey types.
In settings of low P. falciparum transmission, serological data can augment epidemiological information beyond what tests for active malaria infection can provide. While various methodologies can use serological data to produce malaria exposure estimates, sampling at places of congregation could provide a pragmatic approach as countries seek elimination strategies.
Conclusions
Along with carefully considered epidemiological designs, potential assays that include more immunogenic antigen targets, which yield higher population rates of seropositivity, are a feasible approach to discern the foci of residual P. falciparum transmission in an area. Venue-based sampling methods to generate serological data may also have utility for non-Plasmodium infectious diseases, and this strategy should be tested further in other global settings.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Abbreviations
- AMA1:
-
Apical membrane antigen 1
- CDC:
-
Centers for disease control and prevention
- CSP:
-
Circumsporozoite protein
- DBS:
-
Dried blood spot
- IgG:
-
Immunoglobulin G
- INLA:
-
Integrated nested Laplace approximation
- LSA-1:
-
Liver-stage antigen-1
- MBA:
-
Multiplex bead assay
- MFI:
-
Median fluorescence intensity
- MFI-bg:
-
Median fluorescence intensity minus background
- MSP1:
-
Merozoite surface protein 1
- PSI:
-
Population Services International
- RDT:
-
Rapid diagnostic test
- SCR:
-
Serological conversion rate
- SDE:
-
Sections d'énumération
- SPDE:
-
Stochastic partial difference equation
- SRR:
-
Serological reversion rate
- TAS:
-
Transmission assessment survey
- TRaC:
-
Tracking Results Continuously
- WHO:
-
World Health Organization
References
World Health Organization. World malaria report 2022. Geneva: World Health Organization; 2022.
Shretta R, Liu J, Cotter C, Cohen J, Dolenz C, Makomva K, et al. Malaria elimination and eradication. In: Holmes KK, Bertozzi S, Bloom BR, Jha P, editors. Major Infectious Diseases. 3rd ed. Washington (DC)2017.
World Health Organization. Zeroing in on malaria elimination. Geneva: World Health Organization; 2021.
Hay SI, Smith DL, Snow RW. Measuring malaria endemicity from intense to interrupted transmission. Lancet Infect Dis. 2008;8(6):369–78.
World Health Organization. Disease surveillance for malaria elimination: an operational manual. Geneva: World Health Organization; 2012.
World Health Organization. Malaria surveillance, monitoring & evaluation: a reference manual. Geneva: World Health Organization; 2018.
Cotter C, Sturrock HJW, Hsiang MS, Liu J, Phillips AA, Hwang J, et al. The changing epidemiology of malaria elimination: new strategies for new challenges. Lancet. 2013;382(9895):900–11.
Bousema T, Youssef RM, Cook J, Cox J, Alegana VA, Amran J, et al. Serologic markers for detecting malaria in areas of low endemicity, Somalia, 2008. Emerg Infect Dis. 2010;16(3):392–9.
James S, Takken W, Collins FH, Gottlieb M. Needs for monitoring mosquito transmission of malaria in a pre-elimination world. Am J Trop Med Hyg. 2014;90(1):6–10.
Akpogheneta OJ, Duah NO, Tetteh KK, Dunyo S, Lanar DE, Pinder M, et al. Duration of naturally acquired antibody responses to blood-stage Plasmodium falciparum is age dependent and antigen specific. Infect Immun. 2008;76(4):1748–55.
França CT, Li Wai Suen CSN, Carmagnac A, Lin E, Kiniboro B, Siba P, et al. IgG antibodies to synthetic GPI are biomarkers of immune-status to both Plasmodium falciparum and Plasmodium vivax malaria in young children. Malar J. 2017;16(1):386.
Corran P, Coleman P, Riley E, Drakeley C. Serology: a robust indicator of malaria transmission intensity? Trends Parasitol. 2007;23(12):575–82.
Greenhouse B, Smith DL, Rodríguez-Barraquer I, Mueller I, Drakeley CJ. Taking sharper pictures of Malaria with CAMERAs: Combined antibodies to measure exposure recency assays. Am J Trop Med Hyg. 2018;99(5):1120–7.
Ondigo BN, Hodges JS, Ireland KF, Magak NG, Lanar DE, Dutta S, et al. Estimation of recent and long-term malaria transmission in a population by antibody testing to multiple Plasmodium falciparum antigens. J Infect Dis. 2014;210(7):1123–32.
Rogier E, Nace D, Dimbu PR, Wakeman B, Pohl J, Beeson JG, et al. Framework for characterizing longitudinal antibody response in children after plasmodium falciparum infection. Front Immunol. 2021;12: 617951.
Drakeley CJ, Corran PH, Coleman PG, Tongren JE, McDonald SL, Carneiro I, et al. Estimating medium- and long-term trends in malaria transmission by using serological markers of malaria exposure. Proc Natl Acad Sci U S A. 2005;102(14):5108–13.
Fujii Y, Kaneko S, Nzou SM, Mwau M, Njenga SM, Tanigawa C, et al. Serological surveillance development for tropical infectious diseases using simultaneous microsphere-based multiplex assays and finite mixture models. PLoS Negl Trop Dis. 2014;8(7): e3040.
Sarr JB, Orlandi-Pradines E, Fortin S, Sow C, Cornelie S, Rogerie F, et al. Assessment of exposure to Plasmodium falciparum transmission in a low endemicity area by using multiplex fluorescent microsphere-based serological assays. Parasit Vectors. 2011;4:212.
Assefa A, Ali Ahmed A, Deressa W, Sime H, Mohammed H, Kebede A, et al. Multiplex serology demonstrate cumulative prevalence and spatial distribution of malaria in Ethiopia. Malar J. 2019;18(1):246.
Edwards HM, Dixon R, Zegers de Beyl C, Celhay O, Rahman M, Myint Oo M, et al. Prevalence and seroprevalence of Plasmodium infection in Myanmar reveals highly heterogeneous transmission and a large hidden reservoir of infection. PLoS One. 2021;16(6):e0252957.
Stevenson JC, Stresman GH, Gitonga CW, Gillig J, Owaga C, Marube E, et al. Reliability of school surveys in estimating geographic variation in malaria transmission in the western Kenyan highlands. PLoS ONE. 2013;8(10): e77641.
Takem EN, Affara M, Amambua-Ngwa A, Okebe J, Ceesay SJ, Jawara M, et al. Detecting foci of Malaria transmission with school surveys: A pilot study in the Gambia. PLoS ONE. 2013;8(6): e67108.
Brooker S, Kolaczinski JH, Gitonga CW, Noor AM, Snow RW. The use of schools for malaria surveillance and programme evaluation in Africa. Malar J. 2009;8(1):231.
Lucchi NW, Karell MA, Journel I, Rogier E, Goldman I, Ljolje D, et al. PET-PCR method for the molecular detection of malaria parasites in a national malaria surveillance study in Haiti, 2011. Malar J. 2014;13:462.
Oviedo A, Knipes A, Worrell C, Fox LM, Desir L, Fayette C, et al. Combination of Serological, Antigen Detection, and DNA data for plasmodium falciparum provides robust geospatial estimates for Malaria transmission in Haiti. Sci Rep. 2020;10(1):8443.
Lemoine JF, Boncy J, Filler S, Kachur SP, Fitter D, Chang MA. Haiti's Commitment to Malaria Elimination: Progress in the Face of Challenges, 2010–2016. Am J Trop Med Hyg. 2017;97(4_Suppl):43–8.
Knipes AK, Lemoine JF, Monestime F, Fayette CR, Direny AN, Desir L, et al. Partnering for impact: Integrated transmission assessment surveys for lymphatic filariasis, soil transmitted helminths and malaria in Haiti. PLoS Negl Trop Dis. 2017;11(2): e0005387.
WHO. Lymphatic filariasis: monitoring and epidemiological assessment of mass drug administration. https://apps.who.int/iris/handle/10665/44580. 2011.
Priest JW, Jenks MH, Moss DM, Mao B, Buth S, Wannemuehler K, et al. Integration of multiplex bead assays for parasitic diseases into a national, population-based serosurvey of women 15–39 years of age in Cambodia. PLoS Negl Trop Dis. 2016;10(5): e0004699.
Plucinski MM, Candrinho B, Chambe G, Muchanga J, Muguande O, Matsinhe G, et al. Multiplex serology for impact evaluation of bed net distribution on burden of lymphatic filariasis and four species of human malaria in northern Mozambique. PLoS Negl Trop Dis. 2018;12(2): e0006278.
Rogier E, Wiegand R, Moss D, Priest J, Angov E, Dutta S, et al. Multiple comparisons analysis of serological data from an area of low Plasmodium falciparum transmission. Malar J. 2015;14:436.
Oviedo A, Herman C, Knipes A, Worrell CM, Fox LM, Desir L, et al. Spatial cluster analysis of Plasmodium vivax and P. malariae exposure using serological data among Haitian school children sampled between 2014 and 2016. PLoS Negl Trop Dis. 2022;16(1):e0010049.
Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 2009;71(2):319–92.
Drakeley CJ, Carneiro I, Reyburn H, Malima R, Lusingu JP, Cox J, et al. Altitude-dependent and -independent variations in Plasmodium falciparum prevalence in northeastern Tanzania. J Infect Dis. 2005;191(10):1589–98.
Byrne I, Cramer E, Nelli L, Rerolle F, Wu L, Patterson C, et al. Characterizing the spatial distribution of multiple malaria diagnostic endpoints in a low-transmission setting in Lao PDR. Front Med (Lausanne). 2022;9: 929366.
Mwenda MC, Fola AA, Ciubotariu II, Mulube C, Mambwe B, Kasaro R, et al. Performance evaluation of RDT, light microscopy, and PET-PCR for detecting Plasmodium falciparum malaria infections in the 2018 Zambia National Malaria Indicator Survey. Malar J. 2021;20(1):386.
Watson OJ, Sumner KM, Janko M, Goel V, Winskill P, Slater HC, et al. False-negative malaria rapid diagnostic test results and their impact on community-based malaria surveys in sub-Saharan Africa. BMJ Glob Health. 2019;4(4): e001582.
Tadesse FG, Pett H, Baidjoe A, Lanke K, Grignard L, Sutherland C, et al. Submicroscopic carriage of Plasmodium falciparum and Plasmodium vivax in a low endemic area in Ethiopia where no parasitaemia was detected by microscopy or rapid diagnostic test. Malar J. 2015;14(1):303.
Okell LC, Bousema T, Griffin JT, Ouédraogo AL, Ghani AC, Drakeley CJ. Factors determining the occurrence of submicroscopic malaria infections and their relevance for control. Nat Commun. 2012;3(1):1237.
Cook J, Reid H, Iavro J, Kuwahata M, Taleo G, Clements A, et al. Using serological measures to monitor changes in malaria transmission in Vanuatu. Malar J. 2010;9:169.
Helb DA, Tetteh KK, Felgner PL, Skinner J, Hubbard A, Arinaitwe E, et al. Novel serologic biomarkers provide accurate estimates of recent Plasmodium falciparum exposure for individuals and communities. Proc Natl Acad Sci U S A. 2015;112(32):E4438–47.
Mattah PA, Futagbi G, Amekudzi LK, Mattah MM, de Souza DK, Kartey-Attipoe WD, et al. Diversity in breeding sites and distribution of Anopheles mosquitoes in selected urban areas of southern Ghana. Parasit Vectors. 2017;10(1):25.
Tatem AJ, Gething PW, Smith DL, Hay SI. Urbanization and the global malaria recession. Malar J. 2013;12(1):133.
Wilson ML, Krogstad DJ, Arinaitwe E, Arevalo-Herrera M, Chery L, Ferreira MU, et al. Urban Malaria: understanding its epidemiology, ecology, and transmission across seven diverse ICEMR Network Sites. Am J Trop Med Hyg. 2015;93(3 Suppl):110–23.
Abiodun GJ, Maharaj R, Witbooi P, Okosun KO. Modelling the influence of temperature and rainfall on the population dynamics of Anopheles arabiensis. Malar J. 2016;15(1):364.
Oheneba-Dornyo TV, Amuzu S, Maccagnan A, Taylor T. Estimating the impact of temperature and rainfall on malaria incidence in Ghana from 2012 to 2017. Environ Model Assess. 2022;27(3):473–89.
Arab A, Jackson MC, Kongoli C. Modelling the effects of weather and climate on malaria distributions in West Africa. Malar J. 2014;13(1):126.
Dabaro D, Birhanu Z, Negash A, Hawaria D, Yewhalaw D. Effects of rainfall, temperature and topography on malaria incidence in elimination targeted district of Ethiopia. Malar J. 2021;20(1):104.
Cameron E, Young AJ, Twohig KA, Pothin E, Bhavnani D, Dismer A, et al. Mapping the endemicity and seasonality of clinical malaria for intervention targeting in Haiti using routine case data. eLife. 2021;10:e62122.
Jules JR, Alencar J, Suárez-Mutis MC, Baptiste EJ, Albuquerque H, Rosa-Freitas MG, et al. Malaria in Haiti: A descriptive study on spatial and temporal profile from 2009 to 2018. Rev Soc Bras Med Trop. 2022;55: e0355.
Acknowledgements
We would like to thank the Haitian field teams for their efforts during the surveys. Patrick Lammie (The Task Force for Global Health) deserves special thanks for his contributions to project coordination.
Disclaimer
The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention (CDC).
Funding
This study was supported by The Global Fund for HIV, Malaria, and Tuberculosis in collaboration with Population Services International and the US CDC.
Author information
Authors and Affiliations
Contributions
AJU, formal analysis and manuscript preparation; SEJ, project coordination; DN, data collection; CH, data collection; AK, project coordination; CMW, project coordination, LMF, project coordination; LD, project coordination; CF, project coordination; AJ, project coordination; FM, project coordination; KEM, project coordination; MAC, project coordination; JFL, project coordination; KYW, project coordination; VU, project coordination; ESE, data collection; KR, data collection; ER, project coordination and manuscript preparation. All authors read and approved the final manuscript. ER was employed by the US Centers for Disease Control and Prevention at the time of this work.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The study protocol was approved by the National Bioethics Committee of Haiti, and this activity was considered a program evaluation activity by the CDC Human subjects office (#2014–256 and #2015–04). Signed informed consent forms were obtained from all Tracking Results Continuously participants; verbal assent was given by children under 15 years with documented parental consent. All Transmission Assessment Survey participants provided verbal assent; school headmasters were contacted in advance regarding the survey’s purpose and asked to notify parents.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1: Table S1.
Remote sensing data: resolutions, units, and sources. Table S2. Prevalence of positive rapid diagnostic test results with 95% confidence intervals by survey year for Tracking Results Continuously versus Transmission Assessment Surveys. Table S3. Observed seropositivity to P. falciparum antigens by survey year with 95% confidence intervals for community members surveyed for Tracking Results Continuously, 2012-2017. Table S4. Descriptive statistics of survey-site temporal covariates. Table S5. Descriptive statistics of survey-site static covariates. Table S6. Mean estimated coefficients of covariates with 95% credible intervals.
Additional file 2: Fig S1.
Triangulated meshes used to build the stochastic partial difference equation (SPDE) models. Figure S2. Number of participants with a positive rapid diagnostic test (RDT) for malaria by study site. Figure S3. IgG seropositivity by age for Tracking Results Continuously (TRaC) community surveys in Haiti, 2012-2017. Figure S4. Correlation in the observed IgG levels among the LSA-1, MSP1, and AMA1 antigens utilized in this study for both survey types. Figure S5. Correlations between observed and predicted seroprevalence to (A) LSA-1, (B) AMA1, and (C) MSP1 based on model validation. Figure S6. Predicted seroprevalence across Haiti to (A) LSA-1, (B) AMA1, and (C) MSP1, aggregated for each commune.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Jaramillo-Underwood, A., Herman, C., Jean, S.E. et al. Geospatial analysis of Plasmodium falciparum serological indicators: school versus community sampling in a low-transmission malaria setting. BMC Med 22, 31 (2024). https://doi.org/10.1186/s12916-023-03145-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12916-023-03145-6