- Open Access
Serology reveals heterogeneity of Plasmodium falciparum transmission in northeastern South Africa: implications for malaria elimination
Malaria Journalvolume 16, Article number: 48 (2017)
It is widely acknowledged that modifications to existing control interventions are required if South Africa is to achieve malaria elimination. Targeting indoor residual spraying (IRS) to areas where cases have been detected is one strategy currently under investigation in northeastern South Africa. This seroprevalence baseline study, nested within a targeted IRS trial, was undertaken to provide insights into malaria transmission dynamics in South Africa and evaluate whether sero-epidemiological practices have the potential to be routinely incorporated into elimination programmes.
Filter-paper blood spots, demographic and household survey data were collected from 2710 randomly selected households in 56 study wards located in the municipalities of Ba-Phalaborwa and Bushbuckridge. Blood spots were assayed for Plasmodium falciparum apical membrane antigen-1 and merozoite surface protein-119 blood-stage antigens using an enzyme linked immunosorbent assay. Seroprevalence data were analysed using a reverse catalytic model to determine malaria seroconversion rates (SCR). Geospatial cluster analysis was used to investigate transmission heterogeneity while random effects logistic regression identified risk factors associated with malaria exposure.
The overall SCR across the entire study site was 0.012 (95% CI 0.008–0.017) per year. Contrasting SCRs, corresponding to distinct geographical regions across the study site, ranging from <0.001 (95% CI <0.001–0.005) to 0.022 (95% CI 0.008–0.062) per annum revealed prominent transmission heterogeneity. Geospatial cluster analysis of household seroprevalence and age-adjusted antibody responses detected statistically significant (p < 0.05) spatial clusters of P. falciparum exposure. Formal secondary education was associated with lower malaria exposure in the sampled population (AOR 0.72, 95% CI 0.56–0.95, p = 0.018).
Although overall transmission intensity and exposure to malaria was low across both study sites, malaria transmission intensity was highly heterogeneous and associated with low socio-economic status in the region. Findings suggest focal targeting of interventions has the potential to be an appropriate strategy to deploy in South Africa. Furthermore, routinely incorporating sero-epidemiological practices into elimination programmes may prove useful in monitoring malaria transmission intensity in South Africa, and other countries striving for malaria elimination.
Uninterrupted implementation of effective control measures since 2001 has substantially reduced malaria transmission intensity (MTI) in South Africa . In 2012, with a national incidence below the World Health Organization (WHO) malaria elimination threshold of <one case per 1000 population at risk , the Malaria Directorate reoriented programme focus from malaria control to elimination. In support of this shift in focus, a strategic plan aiming to halt local malaria transmission within South African borders by 2018 is now in place . To date, this elimination programme has made steady progress but still faces challenges; one acknowledged obstacle is the identification and targeting of malaria hotspots of transmission as a means of combating the disease .
The ability to stratify malaria risk and target areas with persistent, higher transmission intensities is integral to any effective elimination strategy [5, 6]. Several methods to determine MTI have been developed, although each is associated with its own set of shortcomings and challenges [7–9]. Health facility data, entomological inoculation rates (EIRs) and parasite prevalence are three of the most commonly used surrogate indicators of MTI. However, a reliance on the detection of parasitic material in either humans or mosquitoes reduces the sensitivity of these methods in unstable and low transmission areas, where parasite carriage is generally low [10–12].
Since malaria antibodies persist for extended periods , often irrespective of parasite presence, and can be detected using a high throughput quantitative enzyme-linked immunosorbent assay (ELISA), determining the prevalence of malaria antibodies is a cost-effective, sensitive option for assessing MTI in low-transmission settings [14–17]. Blood stage anti-malaria antibodies have been shown to be reliable markers of malaria exposure  and have been used in numerous sero-epidemiological studies investigating the level of, and changes in, MTI along with geographical patterns in malaria exposure [15, 19–23]. The main, serologically derived, measure of MTI in these studies is the seroconversion rate (SCR). The SCR is the average annual rate at which a population converts from being seronegative to seropositive to malaria antigen(s), and has been used in various geographical settings to accurately determine both current and historical changes in MTI [7, 10, 18, 19]. Across a range of geographical settings, the SCR has been shown to correlate well with both malaria incidence and EIRs [7, 24, 25], although this association has not been validated among very low transmission settings where few serological, entomological and health metric surveys have been conducted simultaneously.
The study reported here forms part of a larger cluster, randomized, clinical trial investigating the efficacy of a targeted method of vector control delivery in response to passively detected cases . This baseline, cross-sectional, seroprevalence survey was undertaken to determine the level and spatial characteristics of malaria transmission, as well as provide insights into the risk factors associated with malaria exposure in northeastern South Africa.
Study site location
The study was conducted in the municipalities of Ba-Phalaborwa and Bushbuckridge, situated in the malaria-endemic provinces of Limpopo and Mpumalanga, South Africa, respectively (Fig. 1). Malaria transmission occurs predominately during the summer months from September to May , with the dominant parasite and vector being Plasmodium falciparum and Anopheles arabiensis, respectively . The primary vector control intervention is annual generalized indoor residual spraying (IRS) of households in all areas in which malaria cases occur, using pyrethroids or dichlorodiphenyltrichloroethane (DDT) .
For the purposes of the targeted IRS, cluster randomized trial, the study area was divided into 56 study wards (13 in Ba-Phalaborwa, 43 in Bushbuckridge) (Fig. 1). Health facilities which operated as passive reporting systems in these wards, accumulated the number of reported malaria cases between 2010 and 2015, which enabled historical incidence to be calculated in the study area. Following collection over 5 years, historical malaria incidence was then averaged in each study ward. Between April and June 2015, a cross-sectional knowledge, attitudes and practices (KAP) and serological survey was undertaken among 60–80 randomly preselected households in each study site ward. Study personnel visited each of the selected households, obtained consent, then administered questionnaires concerning IRS practices, travel history and recent outdoor activity to the head of each household or their representative. Households with no adults present were excluded from the survey. Answers to survey questions were referred to the other consenting individuals residing in the household. Finger-prick dried blood spots (3MM Whatmann Paper) were collected from consenting adults (aged >18 years) and assenting children (aged >2–18 years) along with basic demographic information (age, gender education, employment status). Children under 2 years old were excluded from the survey. Parasite prevalence was not investigated in this study due to the very low reported mean incidence in the region (approximately one per 1000 per annum), which makes prevalence surveys impractical.
Assay of P. falciparum antibodies
Bloodspots were stored at 4 °C as described previously . Sera, eluted from the filter paper blood spots, were subjected to a previously described quantitative indirect ELISA to detect IgG antibodies to blood stage P. falciparum apical membrane antigen-1 (PfAMA-1) and merozoite surface protein-119 (PfMSP-119) antigens . In short, antibodies in 3-mm circular cuttings from filter paper blood spots were eluted in a 1/200 dilution of reconstitution buffer (phosphate buffer saline + 0.5% Tween-20 + 0.05% sodium azide). Sera were then assayed in duplicate for antibodies against both PfAMA-1 and PfMSP-119 in 1/2000 and 1/1000 dilutions, respectively, to obtain two optical density (OD) values per sample. A serial dilution of pooled sera from a malaria hyperendemic African region was used as a positive control to standardize OD values.
Duplicate sample OD values were averaged and normalized against a positive control curve generated from hyperimmune sera. Sample OD values that differed more than 50% were dropped, and when possible repeated. Samples were then dichotomized as either seropositive or seronegative using a two component finite mixture model with five standard deviations as described in [7, 25]. Seropositive was defined as being positive to either PfAMA-1 and/or PfMSP-119. SCRs were obtained from age-specific seroprevalence curves generated from reverse catalytic models, while age-adjusted antibody responses for PfAMA-1 and PfMSP-119 were derived from log-transformed, normalized OD values as previously described in [15, 29]. Age-adjusted antibody responses were then averaged per household and categorized, based equally around the median, as ‘lower than average’, ‘average’, ‘slightly higher than average’, ‘higher than average’, and ‘much higher than average’ to generate an antibody response heat map. Study ward percentage seroprevalence corresponds to the percentage of sampled seropositive households (households that contained at least one member seropositive to PfAMA-1 and/or PfMSP-119) in each ward.
Statistical analysis was performed using Stata 13.0 (Stata Corp, College Station, TX, USA) to identify potential risk factors associated with P. falciparum exposure among sampled participants. Odds ratios (ORs) associated with being seropositive to either antigen were derived from random effects logistic regression models which adjusted for correlation at the study ward level. Correlation was adjusted for at the study ward level, opposed to the household level, as any correlation at household level is reflected in the model. Correlation at ward level therefore provides a more conservative assessment of significance. Adjusted odds ratios (AORs) were derived using a multivariable model, including the following covariates: age, gender, education status, employment status, IRS practices, travel history, night-time outdoor activity, elevation, and study site (Ba-Phalaborwa and Bushbuckridge). Each statistically significant multivariate model was developed using the forward stepwise approach which compared multivariate models to the most significant univariate model using p-values derived from likelihood ratio tests.
Household elevation was determined using point sampling analyses in ArcGIS (v10.3.1). Sampled households with corresponding GPS coordinates were provided with elevation values derived from a 90-m resolution raster elevation data file (DIVA-GIS). Using the same raster elevation data file, average study ward elevation was calculated using zonal statistical methods in ArcGIS (v10.3.1). For geospatial analysis of seroprevalence, households containing only seronegative individuals, at least one seropositive individual as well as household averaged, age-adjusted antibody responses to P. falciparum PfAMA-1 and PfMSP-119 were geographically plotted using ArcGIS software (v10.3.1). The spatial software SaTScan (v.9.4.2) was used to detect clusters of seropositive households and higher than average age-adjusted antibody responses to PfAMA-1 and PfMSP-119. In order to detect clusters of seropositive households (households with at least one seropositive individual), seropositive households were used as cases, while seronegative households were used as controls using the Bernoulli model. Higher than average, age-adjusted PfAMA-1 antibody responses were detected using the Normal model. The scan statistic involves a scanning window, which enumerates both the observed and expected cases (seropositive households or higher than average individual age adjusted responses) across the study area to calculate non-overlapping, statistically significant (p < 0.05) clusters with a maximum set radius of 3 km.
Of the 3522 visited households, 76.9% (2710) agreed to participate in the survey. Within these households, 56.7% (4948/8728) of the eligible participants were present and consented to filter-paper blood spot and demographic data collection. Among the non-consenting participants, 42.5% (425/999) refused without providing a reason. Serological analysis was successfully conducted on 97% (4783/4948) of the blood samples that were correctly labelled. Once assayed, 94.1% (4499/4783) of the serological data correctly matched with demographic data to enable serological analysis. This merged data then matched with 99.5% (4477/4499) of the survey data and were subjected to further statistical and geospatial analysis.
Adult females were over-represented in this study with 75% (2651/3521) of those aged over 18 years being female and only 20% (956/4477) of all participants aged between two and 17 years. Unemployment was high among the working-age population with 85% (2601/3057) of those aged between 18 and 65 years reporting as unemployed. Formal education was limited among the adult population with 40% (1401/3521) of the study participants over 18 years having received no formal secondary education. Night-time outdoor activity was uncommon with only 21% (945/4477) of the surveyed population reported undertaking night-time outdoor activities the evening before the survey. Only 1.5% (65/4477) of the surveyed individuals reported travelling outside South Africa in the past 6 months. Lastly IRS activity was not widespread among the randomly selected households, only 41% (1109/2698) of households reported receiving IRS during the previous malaria season (Table 1).
Plasmodium falciparum antibody responses and associated risk factors
Seropositive cut-offs corresponding to the OD values of 0.067 and 0.103 for PfAMA-1 and PfMSP-119, respectively, dichotomized the sample population as seropositive or seronegative to P. falciparum. Overall antibody responses to both PfAMA-1 and PfMSP-119 were low among the 4499 study participants as no OD value exceeded a value of 0.5 and many seropositive OD values, particularly among those under 10 years, were just above the seropositive thresholds (Additional file 1). Approximately 13% (587/4499) of all sampled individuals were seropositive to PfAMA-1 and/or PfMSP-119, with 9% (403/4499) and 6% (279/4499) seropositive to just PfAMA-1 and PfMSP-119 antigens, respectively, (Additional file 2).
Seroprevalence increased with age across both study sites. Statistically, individuals aged between 2 and 5 years were less likely to be P. falciparum seropositive compared to individuals aged ≥41 years (AOR 0.22, 95% CI 0.13–0.38, p < 0.001) (Table 1). P. falciparum seroprevalence was lower among individuals who attained formal secondary education or above than those who achieved up to formal primary education (AOR 0.72, 95% CI 0.56–0.95, p = 0.018) (Table 1). There was weak evidence for an association between P. falciparum seroprevalence and being female, but no evidence for an association with employment status, recent IRS, recent travel outside of South Africa, night-time outdoor activity or household elevation (Table 1). Furthermore, no statistical association was observed between these risk factors and seroprevalence to just AMA-1 or MSP-119 (Additional file 3).
Plasmodium falciparum transmission intensity and spatial distribution
MTI, as expressed by SCR, was estimated at 0.012 per annum across the entire study site (95% CI 0.008–0.017) (Fig. 2a). SCRs for separate antigens: PfAMA-1 and PfMSP-119 equated to 0.011 (95% CI 0.007–0.017) and 0.003 (95% CI 0.002–0.005), respectively, per year (Additional file 2). MTI however, was not uniform across the study site, as SCRs corresponding to geographical distinct regions (grouped study site wards) varied considerably, ranging from <0.001 (95% CI <0.001–0.005) to 0.022 (95% CI 0.008–0.062) per annum (Fig. 2b, c). In the Ba-Phalaborwa study site alone, two regional SCR values were significantly different, equalling <0.001 (95% CI <0.001–0.005) and 0.009 (95% CI 0.006–0.016) per year (Fig. 2b).
Across both study sites, five statistically significant (p < 0.05) clusters of seropositive households (households containing at least one seropositive individual) were detected (Fig. 3). Clusters of higher than average, age-adjusted antibody responses to both PfAMA-1 and PfMSP-119 were revealed in similar locations to clusters of seropositive households (Fig. 4). In total, 4/5 and 3/4 clusters of higher than average, age-adjusted antibody responses to PfAMA-1 and PfMSP-119, respectively, spatially overlapped with clusters of seropositive households. This spatial overlap was less apparent between age-adjusted antibody responses to either PfAMA-1 or PfMSP-119. In Ba-Phalaborwa, a cluster of higher than average, age-adjusted antibody responses to PfAMA-1 was identified but not to PfMSP-119. Unlike seroprevalence among all ages, malaria seroprevalence in children under 5 years showed no statistical evidence (p > 0.05) of clustering. Nonetheless most children under 5 years were surveyed in houses situated in clusters of seropositive households containing all ages (Fig. 3a, b).
Ward-level malaria incidence, elevation and Plasmodium falciparum seroprevalence
Very low malaria incidence between 2010 and 2015, averaging 0.95 cases per 1000 population, was reported across the whole study region, with values ranging between 0.1 and 3.8 per 1000 population among individual study site wards. Moreover, historical malaria incidence was heterogenous across the study region as wards, that reported <one case per 1000 population, were often situated adjacent to each other (Fig. 1). Variation in elevation was also observed across the study site, with household elevation ranging from 363 to 815 m. Average study ward historical incidence and elevation was inversely associated (slope: −0.006, r: −0.58), highlighting lower incidence at higher ground elevation (Fig. 5). No statistically significant linear relationship was observed between historical malaria incidence and ward-level P. falciparum seroprevalence (Fig. 6).
This cross-sectional, baseline survey, nested within an ongoing targeted IRS trial , aimed to investigate the transmission dynamics, spatial distribution and risk factors associated with P. falciparum exposure within the study area using serological markers, the results of which are to be used to both inform elimination efforts in South Africa and assess whether sero-epidemiological methods have the potential to be routinely incorporated into elimination programmes.
The overall low malaria risk in the northeast border regions of South Africa has been well documented , however until now no detailed population-based assessment of transmission intensity had been conducted. This first sero-epidemiological, cross-sectional study investigating exposure to malaria, estimated a P. falciparum SCR of 0.012 per year, which translates to only 1.2/100 sampled participants becoming seropositive to malaria every year. This result suggests a low force of malaria infection in this region of South Africa, which is expected given the very low reported incidence in the area of approximately one case per 1000 per year and is similar to serological findings in nearby Swaziland . Although a SCR revealed overall MTI to be low across the entire study area, variations in malaria incidence among study site wards and contrasting SCRs between geographically distinct geographical regions highlight prominent heterogeneity of MTI in this area of South Africa. Similar spatial patterns in malaria transmission have been observed in other low transmission settings [5, 18, 21, 24, 32], although few have been conducted in very low endemic settings. It is therefore important to note that while measures of MTI, including SCRs, EIRs and incidence data, recorded over vast geographical areas provide a useful overall picture of malaria transmission, these measures do not necessarily reflect transmission intensity at the micro-epidemiological level. In the Ba-Phalaborwa study site for instance, an area that spans approximately 10 × 20 km, most serological evidence of elevated transmission intensity is focused in the northern half of this region. This suggests a MTI value across this area would not be representative of true transmission at ground level.
Spatial clusters of P. falciparum seropositive households and elevated age-adjusted P. falciparum antibody responses suggest malaria transmission is not only heterogeneous across this study region, but concentrated in malaria hotspots. Hotspots are defined as areas which experience a significantly higher burden of malaria compared to the surrounding vicinity and may act as geographical reservoirs of the disease [6, 15, 29, 33]. To further characterize identified hotspots in the study region, malaria seroprevalence in study participants aged 5 years and under was used as a proxy for detecting recent malaria exposure. Households with seropositive individuals aged ≤5 years were predominately located within hotspots of seropositive households containing all ages, implying these malaria hotspots are likely persistent contributors to malaria transmission and not representative of historical exposure. It should be noted however that many OD values among seropositive children are only just above the seropositive cut-off value. Moreover, elevated antibody responses in children could be attributed to hypergammaglobulinaemia , a process in which antibody production is accelerated due to immunological exposure to other antigens.
Previous studies have identified numerous favourable abiotic and biotic conditions that are believed to facilitate the existence of malaria hotspots, ranging from mosquito breeding site proximity to poor health care practices [5, 35–37]. In this study, poor formal education, a likely proxy for low socio-economic status (SES), was associated with increased malaria exposure. This association is consistent with the rationale that poor quality housing better enables mosquitoes to infect inhabitants [36, 38]. Across the study region, increasing elevation was inversely associated with increasing historical malaria incidence, which is expected given lower temperatures at higher attitudes are believed to impede parasite sporogony in the mosquito gut leading to reduced transmission . The association between elevation and malaria seroprevalence in this study however was more ambiguous, as no declining trend in malaria seroprevalence with increasing elevation was observed. This observation may suggest other factors, aside from elevation, including malaria importation, and successful control efforts have greater impacts on transmission dynamics in lower transmission settings, although this requires further validation. Despite the ambiguity, it should be noted that both historical malaria incidence and malaria seroprevalence was reduced, although not significantly, above altitudes of 650 m.
In this survey, historical malaria incidence represents symptomatic cases that reported to health facilities between 2010 and 2015 while seroprevalence reflects those who have been exposed to the P. falciparum parasite, potentially decades prior. Thus, the observed non-linear relationship between increasing average ward-level incidence and increasing ward-level percentage seroprevalence could be explained by this study detecting previously exposed and/or asymptomatic individuals who may not have reported to health facilities between 2010 and 2015. If this discordance between exposure and incidence is confirmed in more extensive surveys, it raises a question for targeted vector control strategies: should these be directed at foci of higher case incidence or should they be directed at foci of higher exposure to malaria?
Following antigen stratification, seroprevalence and overall seroconversion was higher for PfAMA-1 compared to PfMSP-119, a result shown elsewhere and thought to be the result of contrasting immunogenicities between the two antigens [7, 32]. This disparity may account for the observed spatial discordance in age-adjusted antibody responses towards each antigen in certain study site areas. In the Ba-Phalaborwa study site for instance, a cluster of elevated PfAMA-1, but not PfMSP-119, antibody responses was detected. This discordance would explain why seroprevalence was lower in the Ba-Phalaborwa study site and highlights the importance of using both these serological markers to more effectively characterize malaria exposure, particularly in low endemic settings where antibody responses are low. Nevertheless, given antigenic discordance has been observed elsewhere , this puts into question whether these two serological markers, alone, truly characterize malaria exposure correctly. Fortunately, identifying novel serological markers that accurately represent recent exposure to malaria is an area of continued investigation .
Oversampling of unemployed adult females in this study may have resulted in results that are not fully representative of the study population. Most notable is recent cross-border travel not being significantly associated with malaria exposure, despite a higher proportion of reported travellers being seropositive to malaria. This is attributed to few reports of travel, likely due to the lack of adult males surveyed who are thought to undertake more cross-border travel than adult females . The low sampling rate of males may also have led to the finding of weak evidence that seroprevalence is higher in females than in males. Oversampling is likely a result of survey timing, as most surveys were conducted during weekdays when children and employed adult males were often absent.
This study verifies that serology is effective in characterizing malaria transmission dynamics, even in a very hypo-endemic setting where antibody responses to malaria are low. Therefore, by routinely incorporating sero-epidemiological practices into elimination programmes, this offers countries embarking on malaria elimination the opportunity to both characterize and monitor malaria transmission. This survey, conducted in the northeastern region of South Africa, revealed MTI is both low and spatially heterogenous. This reinforces the notion that malaria elimination in South Africa is achievable, provided amendments to existing methods of control are adopted. Such prominent transmission heterogeneity suggests widespread control interventions are likely to have varying degrees of success while targeting interventions has the potential to be a more appropriate, cost-effective and sustainable strategy for eliminating malaria. However, whether this strategy can cope with malaria epidemics remains unknown, and should be taken into consideration when South Africa pilots a new enhanced surveillance and response strategy in selected malaria hotspots.
adjusted odds ratio
entomological inoculation rate
enzyme-linked immunosorbent assay
indoor residual spraying
knowledge, attitudes and practice
malaria transmission intensity
Plasmodium falciparum apical membrane antigen-1
- PfMSP-119 :
Plasmodium falciparum merozoite surface protein-119
World Health Organization
Maharaj R, Morris N, Seocharan I, Kruger P, Moonasar D, Mabuza A, et al. The feasibility of malaria elimination in South Africa. Malar J. 2012;11:423.
WHO. Malaria elimination: a field manual for low and moderate endemic countries. Geneva: World Health Organization; 2007. http://apps.who.int/iris/bitstream/10665/43796/1/9789241596084_eng.pdf. Accessed 20 July 2016.
South Africa National Department of Health. Malaria elimination strategy for South Africa 2012–2018, Pretoria: NDoH; 2012 http://www.ccoat.uct.ac.za/sites/default/files/image_tool/images/229/SA%20Malaria%20Elimination%20Strategic%20Plan%20%202012.pdf. Accessed 17 July 2016.
Raman J, Morris N, Frean J, Brooke B, Blumberg L, Kruger P, et al. Reviewing South Africa’s malaria elimination strategy (2012–2018): progress, challenges and priorities. Malar J. 2016;15:438.
Rosas-Aguirre A, Speybroeck N, Llanos-Cuentas A, Rosanas-Urgell A, Carrasco-Escobar G, Rodriguez H, et al. Hotspots of malaria transmission in the Peruvian amazon: rapid assessment through a parasitological and serological survey. PLoS ONE. 2015;10:e0137458.
Bousema T, Griffin JT, Sauerwein RW, Smith DL, Churcher TS, Takken W, et al. Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Med. 2012;9:e1001165.
Stewart L, Gosling R, Griffin J, Gesase S, Campo J, Hashim R, et al. Rapid assessment of malaria transmission using age-specific sero-conversion rates. PLoS ONE. 2009;4:e6083.
Wu L, van den Hoogen LL, Slater H, Walker PGT, Ghani AC, Drakeley CJ, et al. Comparison of diagnostics for the detection of asymptomatic Plasmodium falciparum infections to inform control and elimination strategies. Nature. 2015;528:86–93.
Kelly-Hope LA, McKenzie FE. The multiplicity of malaria transmission: a review of entomological inoculation rate measurements and methods across sub-Saharan Africa. Malar J. 2009;8:19.
Corran P, Coleman P, Riley E, Drakeley C. Serology: a robust indicator of malaria transmission intensity? Trends Parasitol. 2007;23:575–82.
Snow RW, Omumbo JA, Lowe B, Molyneux CS, Obiero JO, Palmer A, et al. Relation between severe malaria morbidity in children and level of Plasmodium falciparum transmission in Africa. Lancet. 1997;349:1650–4.
Drakeley C, Schellenberg D, Kihonda J, Sousa CA, Arez AP, Lopes D, et al. An estimation of the entomological inoculation rate for Ifakara: a semi-urban area in a region of intense malaria transmission in Tanzania. Trop Med Int Health. 2003;8:767–74.
Actman A, Bull P, Stephens R, Langhorne J. Longevity of the immune response and memory to blood-stage malaria infection. Curr Top Micobiol Immunol. 2005;297:71–205.
van den Hoogen LL, Griffin JT, Cook J, Sepúlveda N, Corran P, Conway DJ, et al. Serology describes a profile of declining malaria transmission in Farafenni, The Gambia. Malar J. 2015;14:416.
Cook J, Kleinschmidt I, Schwabe C, Nseng G, Bousema T, Corran PH, et al. Serological markers suggest heterogeneity of effectiveness of malaria control interventions on Bioko Island, Equatorial Guinea. PLoS ONE. 2011;6:e25137.
Teirlinck AC, McCall MBB, Roestenberg M, Scholzen A, Woestenenk R, de Mast Q, et al. Longevity and composition of cellular immune responses following experimental Plasmodium falciparum malaria infection in humans. PLoS Pathog. 2011;7:e1002389.
Yeka A, Nankabirwa J, Mpimbaza A, Kigozi R, Arinaitwe E, Drakeley C, et al. Factors associated with malaria parasitemia, anemia and serological responses in a spectrum of epidemiological settings in Uganda. PLoS ONE. 2015;10:e0118901.
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;13:392–9.
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.
Arnold BF, Priest JW, Hamlin KL, Moss DM, Colford JM, Lammie PJ. Serological measures of malaria transmission in haiti: comparison of longitudinal and cross-sectional methods. PLoS ONE. 2014;9:e93684.
Lynch C, Cook J, Nanyunja S, Bruce J, Bhasin A, Drakeley C, et al. Application of serological tools and spatial analysis to investigate malaria transmission dynamics in highland areas of Southwest Uganda. Am J Trop Med Hyg. 2016;94:1251–8.
Rosas-Aguirre A, Llanos-Cuentas A, Speybroeck N, Cook J, Contreras-Mancilla J, Soto V, et al. Assessing malaria transmission in a low endemicity area of north-western Peru. Malar J. 2013;12:339.
von Fricken ME, Weppelmann TA, Lam B, Eaton WT, Schick L, Masse R, et al. Age-specific malaria seroprevalence rates: a cross-sectional analysis of malaria transmission in the Ouest and Sud-Est departments of Haiti. Malar J. 2014;13:361.
Oduro AR, Conway DJ, Schellenberg D, Satoguina J, Greenwood BM, Bojang KA. Seroepidemiological and parasitological evaluation of the heterogeneity of malaria infection in the Gambia. Malar J. 2013;12:222.
Corran PH, Cook J, Lynch C, Leendertse H, Manjurano A, Griffin J, et al. Dried blood spots as a source of anti-malarial antibodies for epidemiological studies. Malar J. 2008;7:195.
Kleinschmidt I, Coetzee M. From malaria control to sustainable elimination: cluster randomised trial comparing targeted versus generalised vector control in South Africa [Internet]. 2015. ClinicalTrials.gov:NCT02556242.
Khosa E, Kuonza LR, Kruger P, Maimela E. Towards the elimination of malaria in South Africa: a review of surveillance data in Mutale Municipality, Limpopo Province, 2005 to 2010. Malar J. 2013;12:7.
Maharaj R, Raman J, Morris N, Moonasar D, Durrheim DN, Seocharan I, et al. Epidemiology of malaria in South Africa: from control to elimination. S Afr Med J. 2013;103:779–83.
Wilson S, Booth M, Jones FM, Mwatha JK, Kimani G, Kariuki HC, et al. Age-adjusted Plasmodium falciparum antibody levels in school-aged children are a stable marker of microgeographical variations in exposure to Plasmodium infection. BMC Infect Dis. 2007;7:67.
Morris N, Frean J, Baker L, Ukpe IS, Barnes KI, Kruger P, et al. Re-defining the extent of malaria transmission in South Africa: implications for chemoprophylaxis. S Afr Med J. 2013;103:861–4.
Hsiang MS, Hwang J, Kunene S, Drakeley C, Kandula D, Novotny J, et al. Surveillance for malaria elimination in Swaziland: a national cross-sectional study using pooled PCR and serology. PLoS ONE. 2012;7:e29550.
Wong J, Hamel MJ, Drakeley CJ, Kariuki S, Shi YP, Lal AA, et al. Serological markers for monitoring historical changes in malaria transmission intensity in a highly endemic region of Western Kenya, 1994–2009. Malar J. 2014;13:451.
Sturrock HJW, Hsiang MS, Cohen JM, Smith DL, Greenhouse B, Bousema T, et al. Targeting asymptomatic malaria infections: active surveillance in control and elimination. PLoS Med. 2013;10:e1001467.
Langhorne J, Ndungu FM, Sponaas AM, Marsh K. Immunity to malaria: more questions than answers. Nat Immunol. 2008;9:725–32.
Ndiath M, Faye B, Cisse B, Ndiaye JL, Gomis JF, Dia AT, et al. Identifying malaria hotspots in Keur Soce health and demographic surveillance site in context of low transmission. Malar J. 2014;13:453.
Kumar DS, Andimuthu R, Rajan R, Venkatesan MS. Spatial trend, environmental and socioeconomic factors associated with malaria prevalence in Chennai. Malar J. 2014;13:14.
Ernst KC, Adoka SO, Kowuor DO, Wilson ML, John CC. Malaria hotspot areas in a highland Kenya site are consistent in epidemic and non-epidemic years and are associated with ecological factors. Malar J. 2006;5:78.
Tusting LS, Willey B, Lucas H, Thompson J, Kafy HT, Smith R, et al. Socioeconomic development as an intervention against malaria: a systematic review and meta-analysis. Lancet. 2013;382:963–72.
Bødker R, Akida J, Shayo D, Kisinza W, Msangeni HA, Pedersen EM, et al. Relationship between altitude and intensity of malaria transmission in the Usambara Mountains, Tanzania. J Med Entomol. 2003;40:706–17.
Mosha JF, Sturrock HJW, Greenwood B, Sutherland CJ, Gadalla NB, Atwal S, et al. Hotspot or not: a comparison of spatial statistical methods to predict prospective malaria infections. Malar J. 2014;13:53.
Helb DA, Tetteh KKA, 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 USA. 2015;112:e4438–47.
Cotter C, Sturrock HJW, Hsiang MS, Liu J, Philips AA, Hwang J, et al. The changing epidemiology of malaria elimination : new strategies for new challenges. Lancet. 2013;6736:1–12.
Morris N. Health GIS Centre, South Africa Medical Research Council, Durban. 2015.
JR, IK and CD designed the study. JR, KH, NM, IS, PK, AM, AZ, EM, and MC assisted with study implementation. JB, EA and JR processed serological data while NM produced incidence calculations. Data analyses were carried out by JB with assistance from JC, IK and CD. JB and JR drafted the manuscript with help from JC, IK, CD, MC, NM, and KH. All authors read and approved the final manuscript.
The authors wish to thank the communities of Ba-Phalaborwa and Bushbuckridge for their support during this study, the KAP survey field team for their diligent work during the KAP survey, Mr Itumeleng Moletsane, Mr Matimba Zitha and the malaria control staff from Limpopo and Mpumalanga Provinces for logistic support, Tom Hall for his serological advice and the South African Medical Research Council for malaria case data.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets during and/or analysed during the current study are available from IK on reasonable request.
Ethical approval and consent to participate
Written informed consent was obtained from study participants older than 18 years prior to blood sample collection. Although a parent or guardian gave informed written consent for blood collection from individuals younger than 18 years, blood was only collected if assent was obtained from the minor. Approval for the randomized trial (including this baseline survey) was obtained from the Provincial Departments of Health of Mpumalanga and Limpopo Provinces (4/2/2), as well as ethics review boards of the London School of Hygiene and Tropical Medicine (7396-1) and the University of the Witwatersrand (M140762).
This study was financially supported by a UK-MRC/DIFD/Wellcome Trust grant awarded to IK and MC. The funders did not participate in study design, data collection, data analysis and interpretation or in the drafting of the manuscript. Funding was provided by Medical Research Council (Grant No. MR/L00433X/1.)
Joseph Biggs and Jaishree Raman are co-primary authors