Skip to main content

Geospatial analysis of Plasmodium falciparum serological indicators: school versus community sampling in a low-transmission malaria setting

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.

Peer Review reports

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:

$${Y}_{i}|P\left({x}_{i}\right)\sim Binomial\left({N}_{i}, P\left({x}_{i}\right)\right)$$

The seroprevalence at each location was then linked to its linear predictors through the logit:

$$logit\left(P\left({x}_{i}\right)={\beta }_{0}+d{\left({x}_{i}\right)}^{\prime}\beta +{S(x}_{i}\right),$$

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).

Fig. 1
figure 1

Study sites across Haiti. Tracking Results Continuously (TRaC; n = 324) sites are represented by yellow triangles and transmission assessment survey (TAS; n = 1163) sites by blue dots. Names of Haitian departments are provided; outlines represent department boundaries

Table 1 Raw data summary of P. falciparum serology samples collected from two study populations in Haiti

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.

Fig. 2
figure 2

Observed seroprevalence by study site. Maps are grouped based on seroprevalence to A liver stage antigen-1 (LSA-1), B apical membrane antigen 1 (AMA1), and C merozoite surface protein 1 (MSP1). Tracking Results Continuously (TRaC) sites are on the left-hand side and transmission assessment survey (TAS) on the right. Circle sizes are proportional to the number of participants surveyed at each site. Darker red shading indicates higher seroprevalence; white fill indicates zero participants were seropositive at a given site. Department boundaries are outlined

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).

Fig. 3
figure 3

Mean estimated coefficients of covariates with 95% credible intervals for seroprevalence to liver stage antigen-1 (LSA-1), apical membrane antigen 1 (AMA1), and merozoite surface protein 1 (MSP1). All covariates have been centered and standardized. Accessibility to cities is represented by mean travel time to urban centers, such that an increase in this variable indicates a more rural location. Dots represent mean estimated beta coefficients for covariates used in integrated nested Laplace approximation (INLA) models with the stochastic partial difference equation (SPDE) approach; bars are the coefficients’ 95% credible intervals. A vertical line at 0.0 represents a null effect

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).

Fig. 4
figure 4

Predicted seroprevalence across Haiti to A liver stage antigen-1 (LSA-1), B apical membrane antigen 1 (AMA1), and C merozoite surface protein 1 (MSP1). Tracking Results Continuously (TRaC) predictions are on the left panels and transmission assessment survey (TAS) on the right. Darker red to purple shading indicates higher seroprevalence, while paler shades of red to yellow indicate lower seroprevalence. Scale of shading is consistent for the same antigen target between two survey types. Department boundaries are outlined, and scale bar and north arrow illustrated. Pixel size is approximately 0.9 km2. Predicted seroprevalence was modeled by linking the observed seroprevalence at TAS or TRaC survey locations to environmental covariate data spanning all coordinate points across Haiti. Predicted seroprevalence was extracted as a continuous surface for mapping

Table 2 Spearman’s rank correlation coefficient, ρ, and 95% confidence intervals for TAS versus TRaC surveys

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

  1. World Health Organization. World malaria report 2022. Geneva: World Health Organization; 2022.

    Google Scholar 

  2. 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.

  3. World Health Organization. Zeroing in on malaria elimination. Geneva: World Health Organization; 2021.

    Google Scholar 

  4. Hay SI, Smith DL, Snow RW. Measuring malaria endemicity from intense to interrupted transmission. Lancet Infect Dis. 2008;8(6):369–78.

    Article  PubMed  PubMed Central  Google Scholar 

  5. World Health Organization. Disease surveillance for malaria elimination: an operational manual. Geneva: World Health Organization; 2012.

    Google Scholar 

  6. World Health Organization. Malaria surveillance, monitoring & evaluation: a reference manual. Geneva: World Health Organization; 2018.

    Google Scholar 

  7. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  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.

    Article  PubMed  PubMed Central  Google Scholar 

  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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  11. 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.

  12. Corran P, Coleman P, Riley E, Drakeley C. Serology: a robust indicator of malaria transmission intensity? Trends Parasitol. 2007;23(12):575–82.

    Article  PubMed  Google Scholar 

  13. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  14. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  16. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  17. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 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.

  21. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  24. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  25. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. 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.

  27. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  28. WHO. Lymphatic filariasis: monitoring and epidemiological assessment of mass drug administration. https://apps.who.int/iris/handle/10665/44580. 2011.

  29. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  31. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  32. 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.

  33. 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.

    Article  Google Scholar 

  34. 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.

    Article  PubMed  Google Scholar 

  35. 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.

    Article  PubMed  Google Scholar 

  36. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  38. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  39. 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.

    Article  PubMed  Google Scholar 

  40. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Tatem AJ, Gething PW, Smith DL, Hay SI. Urbanization and the global malaria recession. Malar J. 2013;12(1):133.

    Article  PubMed  PubMed Central  Google Scholar 

  44. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 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.

    Article  Google Scholar 

  47. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  48. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  49. 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.

  50. 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.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

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

Authors

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

Correspondence to Eric Rogier.

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-023-03145-6

Keywords