Skip to main content

Analysis of serological data to investigate heterogeneity of malaria transmission: a community-based cross-sectional study in an area conducting elimination in Indonesia

Abstract

Background

Analysis of anti-malarial antibody responses has the potential to improve characterization of the variation in exposure to infection in low transmission settings, where conventional measures, such as entomological estimates and parasitaemia point prevalence become less sensitive and expensive to measure. This study evaluates the use of sero-epidemiological analysis to investigate heterogeneity of transmission in area conducting elimination in Indonesia.

Methods

Filter paper bloodspots and epidemiological data were collected through a community-based cross-sectional study conducted in two sub-districts in Sabang municipality, Aceh province, Indonesia in 2013. Antibody responses to merozoite surface protein 1 (MSP-119) and apical membrane antigen 1 (AMA-1) for Plasmodium falciparum and Plasmodium vivax were measured using indirect enzyme-linked immunosorbent assay (ELISA). Seroconversion rates (SCR) were estimated by fitting a simple reversible catalytic model to seroprevalence data for each antibody. Spatial analysis was performed using a Normal model (SaTScan v.9.4.2) to identify the clustering of higher values of household antibody responses. Multiple logistic regression was used to investigate factors associated with exposure.

Results

1624 samples were collected from 605 households. Seroprevalence to any P. falciparum antigen was higher than to any P. vivax antigen, 6.9% (95% CI 5.8–8.2) vs 2.0% (95% CI 1.4–2.8). SCR estimates suggest that there was a significant change in P. falciparum transmission with no exposure seen in children under 5 years old. Plasmodium falciparum SCR in over 5 years old was 0.008 (95% CI 0.003–0.017) and 0.012 (95% CI 0.005–0.030) in Sukakarya and Sukajaya sub-districts, respectively. Clusters of exposure were detected for both P. falciparum and P. vivax, most of them in Sukajaya sub-district. Higher age, P. vivax seropositivity and use of long-lasting insecticide-treated bed net (LLIN) were associated with higher P. falciparum exposure.

Conclusion

Analysis of community-based serological data helps describe the level of transmission, heterogeneity and factors associated with malaria transmission in Sabang. This approach could be an important additional tool for malaria monitoring and surveillance in low transmission settings in Indonesia.

Background

In recent years, there has been a decline in malaria transmission in many regions, leading to optimism that malaria elimination might be achieved in numerous countries [1,2,3,4,5,6,7,8]. As transmission declines, monitoring changes in malaria transmission intensity and disease prevalence through surveillance systems becomes increasingly important to allow the evaluation of health services and control programs [9, 10]. The latest World Health Organization (WHO) malaria surveillance manual confirms that improved surveillance is a major component of the WHO strategy [11]. However, conventional measures such as entomological estimates and parasitaemia point prevalence become less sensitive and relatively more expensive as transmission declines [12, 13]. Disease surveillance is further compounded by difficult access to remote and isolated communities, increased risks in forest workers and other highly mobile populations and the difficulties of tracking cross-border movements [14,15,16,17,18,19,20].

An additional approach to measure malaria transmission is to detect anti-malarial antibodies, which provide a marker for exposure to malaria [9]. Malaria infections generate antibodies which can be detected for several months and years after the infection has been resolved. Although serology is unlikely to be useful for diagnosing actively infected individuals because antibodies take days to develop and then persist after infection [9, 13], detection of these antibodies indicates previous exposure and offers an additional, more sensitive measure of infection and transmission, particularly in low endemic settings where the sensitivity of parasitological tools is inadequate [21,22,23,24,25] and gold standard tests like the parasite rate and the entomological inoculation rate (EIR), may have insufficient statistical power unless the sampling is intensively done [26,27,28]. This approach has been utilized in several countries and reported as a more sensitive tool to assess population-level malaria exposure in low-transmission settings [9, 13].

Seroconversion rate (the proportion of people in the population who are expected to seroconvert each year) is a serological parameter used to understand malaria transmission dynamics. Previous studies found that seroconversion rate (SCR) provides a proxy measure for estimating the transmission intensity in a community as it was strongly correlated with the EIR and annual parasite incidence collected by the malaria surveillance programme [10, 14]. Serological estimates of transmission have been utilized in many low endemic settings, including Indonesia [29, 30], and have additionally been used to identify populations at higher risk of malaria exposure [9, 31], foci of transmission [32, 33] and to describe historical changes in disease burden [25]. While there is great promise in this approach, it needs further refinement.

Recent studies have reported the potential use of recombinant Merozoite Surface Protein 1 (PfMSP-119) and Apical Membrane Antigen 1 (PfAMA-1) as serological parameters to assess malaria transmission intensity in Indonesia. First, a population-based cross-sectional study conducted in three different endemicity areas showed the potential application of these methods for detecting changes in transmission exposure, particularly in lower transmission settings and with less immunogenic antigens (such as PfMSP-119) [30]. Second, a cohort study of Indonesian schoolchildren found that it is possible to assess the interruption of transmission by measuring seroconversion rates from individual-level longitudinal data on antibody titres [29]. These studies suggested serological analysis has the potential to assess malaria burden and heterogeneity of infections in the Indonesian population. As antibodies to AMA-1 and MSP-119 antigens have been reported to persist for several years after infection and in the absence of reinfection, any antibodies detected in younger children would reflect more recent infection in low transmission settings [10]. Therefore, as Indonesia aims to eliminate malaria by 2030, further implementation and evaluation of sero-epidemiological analysis in areas moving towards elimination would garner valuable information for malaria control programmes. This study explores the use of sero-epidemiological analysis for assessing the intensity and heterogeneity of malaria transmission as well as factors associated with malaria exposure in an area conducting elimination in Indonesia.

Methods

Study site

The study was conducted in Sabang municipality, Aceh province, Indonesia (Fig. 1). The municipality is located at the north-westernmost part of Indonesia and is part of Aceh province. The municipality has an area of 153 km2 covering five islands but only the largest island, Weh, is permanently inhabited. The population on Weh island is approximately 30,000 and it is divided administratively into two sub-districts (Sukakarya and Sukajaya) with 18 villages. Sabang has a very low-level annual parasite incidence, 0.13 per thousand population in 2011. Based on its geographic position at the western end of the archipelago, its diverse mosquito fauna, the presence of both major malaria parasites, and its strong local government, Sabang municipality was considered as an appropriate place to pilot malaria elimination in Indonesia [34].

Fig. 1
figure 1

Map showing study sites in two sub-districts in Sabang municipality, Aceh province, Indonesia (a). Inset maps showing geographical location of Sabang municipality within Aceh province (b), and location of Aceh province within Indonesian archipelago (c)

Study design and data collection

Community-based cross-sectional surveys were performed during the malaria transmission season between October and December 2013. Households list were obtained from local authorities and were arbitrarily assigned numbers according to their geographic location. Households were randomly selected and invited to participate in the study. Households with no adult present were excluded from the survey and were replaced by the neighbouring households. Individual signed consent was obtained from all adults or guardians of household member under 18 years of age. Samples were collected from all household members present aged over 6 months. The minimum sample size of 439 individuals per sub-district was met to ensure the antibody SCR of 0.0036 could be estimated with a precision level of ± 0.0013 [35]. Standard microscopy blood smears were collected as per routine national diagnostic standards. Filter paper bloodspots were collected on Whatman 3 M paper (Whatman, UK) as described by Corran et al. [36] and stored at − 20 °C until transferred to the Parasitology Laboratory at Department of Parasitology, UGM, Yogyakarta. Data on age, gender, education, occupation, long-lasting insecticide-treated bed net (LLIN) use, indoor residual spraying (IRS) in last 12 months and auxiliary temperature were recorded using a short questionnaire form, and household GPS coordinates were collected using handheld GPS.

Laboratory methods

Giemsa-stained thick and thin malaria films reading was performed by trained laboratory technicians to identify active infections. For serological assays, the recombinant proteins Plasmodium falciparum MSP-119, P. falciparum AMA-1, Plasmodium vivax MSP-119 and P. vivax AMA-1 were used as antigens in indirect enzyme-linked immunosorbent assay (ELISA) as described in [9]. Briefly, antigens were coated on 96 well plates at the concentration of 0.5 µg/mL in coating buffer and incubated at 4 °C overnight. The plates were washed in phosphate buffered saline with tween (PBST) and blocked with 1% (w/v) skimmed milk solution for 3 h. After washing, samples were added in duplicate at a final dilution of 1:1000 to each plate using a pool of hyperimmune serum as a positive control and the plates were incubated overnight at 4 °C. The plates were washed and 50 µL of HRP-conjugated rabbit anti Human IgG (DAKO, #P0214) were added into each well and incubated for 3 h. After a further series of washes substrate solution (OPD, Sigma #P8287, in PBS) was added and the reaction was allowed to develop for 15–20 min before addition of stopping solution (2 M H2SO4). The optical density was read using ELISA reader at 450 nm. All serology was performed by trained laboratory technicians at the Department of Parasitology, UGM, Yogyakarta.

Statistical analysis

All statistical analyses were conducted in Stata IC 15 (Stata Corp, College Station, TX, USA). Infants under 1 year of age were excluded from each dataset to remove any influence of maternally derived antibodies [10]. Raw OD measurements were averaged and normalized against the positive control curve on each plate. A cut-off for seropositivity was determined for each antigen by calculating the mean plus 3 standard deviation values of OD values from serum samples of 40 Javanese individuals who had no history of travel to malaria endemic areas in Indonesia. Cut-offs were generated separately for each antigen [13]. Individuals were categorized as seropositive for P. falciparum if their antibody responses were above the cut-off for PfAMA-1 and/or PfMSP-119 and seropositive for P. vivax if their antibody responses were above the cut-off for PvAMA-1 and/or PvMSP-119. Seroconversion rates were estimated by fitting a simple reversible catalytic model to seroprevalence data for each antibody [10]. Models with two SCRs allowing detection of changes in SCR were fitted and a likelihood test ratio was performed to decide the most appropriate model. Bivariate and multivariable analysis were performed to identify potential factors associated with P. falciparum (and P. vivax) exposure among study participants. Logistic regression models were performed to estimate odds ratios (ORs) of factors associated with being seropositive to P. falciparum or P. vivax, respectively. Adjusted odds ratios (aORs) were obtained using a multivariable model, including the following covariates: age, gender, seropositivity to P. vivax, education status, employment status, LLIN use, IRS in last 12 months, fever status, and altitude. Samples from participants aged under 18 years old were excluded from analysis of education and employment status. Statistically significant variables (p < 0.05) detected in bivariate analysis were included in a multivariable model. The final model was developed using the forward stepwise approach which compared multivariable models to the most significant bivariate model using p-values calculated from likelihood ratio tests. Scatter plots matrix and coefficient correlation analysis were done to assess potential cross-reactivity between P. falciparum and P. vivax antigens.

Spatial analysis

The spatial software SaTScan (v.9.4.2) was used to detect clusters of higher than average age-adjusted antibody responses to PfAMA-1, PfMSP-119, PvAMA-1 and PvMSP-119. The Normal model was used to detect clusters of households with higher than average age-adjusted antibody responses to PfAMA-1, PfMSP-119, PvAMA-1, PvMSP-119 antibody responses. This method has been previously utilized in several studies investigating malaria transmission heterogeneity in low endemic setting [25, 37]. Antibody responses data were first log10 transformed and then adjusted for age. The residuals from linear regression (log titre regressed against age in years, performed in Stata IC 15) were used to determine whether antibody responses were higher or lower than expected for any given age assuming a homogeneous distribution of risk. Residuals less than zero represent individuals whose responses were lower than or average for their age group whilst residuals above zero represent individuals whose responses were higher than average. These data 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. The scan statistic was set to calculate non-overlapping, statistically significant (p < 0.05) clusters with a maximum set radius of 3 km and with minimum 2 observations detected in a cluster. Data generated from SatScan were then plotted using ArcGIS software (v10.5).

Results

Study population

General characteristics of the sampled population is presented in Table 1. There were 1624 samples collected in the surveys from 605 households. The average number of people sampled per household was 3 (SD: 1.64). The proportion of females sampled (61%) was slightly higher than males. The majority of the samples came from Sukajaya sub-district (63%), and the median age of participants was 22 years (IQR: 9–38). Educational attainment was high, with only 0.34% of adults ≥ 18 years old who had not completed primary education. More than half (57%) of the working-age population (≥ 18 years old) were unemployed. The population LLIN coverage was 60%, with 68% of those who owned nets reporting to have slept under it the night before. Only 15% of study households had received IRS in the previous 12 months. 9% of the population had fever with body temperature reading > 37.5 °C. Examination of microscopy slides found no malaria infections.

Table 1 Demographic characteristics and factors associated with P. falciparum transmission in Sabang, Indonesia, 2013

Seroprevalence and associated factors

Seropositivity to P. falciparum antigens was higher than seropositivity to P. vivax antigens, with seroprevalence 6.89% (95% CI 5.76–8.24) and 1.97% (95% CI 1.39–2.77), respectively. Seroprevalence ranged from 1.2 to 11.4% for P. falciparum and 0.5 to 2.8% for P. vivax across age groups. Notably, there were no seropositive individuals aged under 5 years old identified for either P. falciparum or P. vivax (0/210). Seroprevalence to each antigen can be found in Additional file 1.

Multivariable analysis in Table 1 shows that age, seropositivity to P. vivax and use of LLINs were significantly associated with P. falciparum seropositivity, after controlling for other covariates. As would be expected, seroprevalence increased with age. Adults were more likely to be seropositive compared to children under 15 years old, with adjusted OR 5.69 (95% CI 2.43–13.37), 12.05 (95% CI 5.59–25.94) and 10.27 (95% CI 4.74–22.27) for age group 16–24, 25–40 and over 40 years old, respectively. Seropositivity to P. falciparum was also significantly associated with higher proportion of LLIN use, with adjusted OR 1.80 (95%: 1.20–2.72). In addition, people who were seropositive to any P. vivax antigen were 3 times more likely to be seropositive for P. falciparum, with adjusted OR 3.47, (95% CI 1.48–8.12). Other factors such as gender, residence, education, employment, IRS, fever and altitude were not significantly associated with P. falciparum seropositivity. Multivariable logistic regression revealed that there were no factors significantly associated with P. vivax seropositivity (Additional file 2).

Plasmodium falciparum and Plasmodium vivax transmission intensity

Figure 2 describes the SCR estimates for P. falciparum and P. vivax in Sukakarya and Sukajaya sub-districts, Sabang municipality, Indonesia in 2013. The SCR estimates suggested that there was a significant change in P. falciparum transmission in both Sukakarya and Sukajaya sub-districts, with no exposure seen in children under 5 years old. The data suggested that the P. falciparum transmission intensity in people aged over 5 years old was SCR 0.008 (95% CI 0.003–0.017) and SCR 0.012 (95% CI 0.005–0.030) in Sukakarya and Sukajaya, respectively. The SCR estimates for P. vivax (Fig. 2c, d) also suggested a very low level of transmission, SCR 0.001 (95% CI 0.000–0.005) and 0.002 (95% CI 0.001–0.006), respectively. There was no evidence for a difference in SCR between people aged under 5 and over 5 years old in either Sukakarya or Sukajaya. Overall, these model SCRs estimates suggested that the magnitude of transmission in population level was likely to be similarly very low for P. falciparum and P vivax.

Fig. 2
figure 2

Age-seroprevalence plots for P. falciparum in Sukakarya (a), Sukajaya (b), for P. vivax in Sukakarya (c) and in Sukajaya (d), 2013. Solid lines represent the fitted probability for being seropositive to either MSP-119 or AMA-1 antigen, dashed lines represent the 95% confidence interval of these fits and red triangles represent the observed proportion of seropositive per age decile. SCR value represent the average annual rate at which the population become seropositive to any of P. falciparum or P. vivax antigen

Heterogeneity of P. falciparum and P. vivax

Spatial analysis of higher than average age-adjusted antibody responses identified 5 significant clusters for PfAMA-1 and 3 clusters for PfMSP-119. All 5 of the PfAMA-1 clusters were seen in Sukajaya (Fig. 3a), whilst 2 of 3 PfMSP-119 clusters seen in Sukajaya and spatially overlapped with the PfAMA-1 clusters (Fig. 3b). One additional PfMSP-119 cluster was identified in Sukakarya. The analysis of age adjusted antibody responses to P. vivax antigens identified 2 clusters for PvAMA-1 and 3 clusters for PvMSP-119 in Sukajaya (Fig. 4). The clusters identified for PvAMA-1 spatially overlapped the PvMSP-119 clusters. All of these P. vivax clusters were seen in Sukajaya, whilst no clusters were identified in Sukakarya. Overall, the clusters identified for P. falciparum and P. vivax were seen in the same areas.

Fig. 3
figure 3

Spatial distribution of household-averaged, age-adjusted antibody responses to a PfAMA-1 and to b PfMSP-119 in Sukakarya and Sukajaya sub-districts, Sabang, Indonesia. The resultants residual values were categorized as: ‘lower than average’ (− 4.326 to − 0.499), ‘average’ (− 0.500 to 0.500), ‘slightly higher than average’ (0.501 to 1.000), ‘higher than average’ (1.001 to 1.500) and much higher than average (1.501 to 2.842). Black circle indicates a cluster of significantly higher than expected antibody responses detected using SaTScan (p value < 0.05)

Fig. 4
figure 4

Spatial distribution of household-averaged, age-adjusted antibody responses to a PvAMA-1 and to b PvMSP-119 in Sukakarya and Sukajaya sub-districts, Sabang, Indonesia. The resultants residual values were categorized as: ‘lower than average’ (− 4.933 to − 0.499), ‘average’ (− 0.500 to 0.500), ‘slightly higher than average’ (0.501 to 1.000), ‘higher than average’ (1.001 to 1.500) and much higher than average (1.501 to 2.117). Black circle indicates a cluster of significantly higher than expected antibody responses detected using SaTScan (p value < 0.05)

Discussion

This study describes the analysis of community-based serological data to investigate malaria transmission dynamics in a low transmission setting, Sabang, Indonesia. The seroprevalence and SCR data represent exposure to infection and demonstrate that the population level of transmission intensity were similarly very low for both P. falciparum and P. vivax. The seroprevalence in children under 15 years old was negligible, 1.2% and 0.5% for P. falciparum and P. vivax, respectively. The spatial analysis of household-level data on antibody responses to any of the antigens tested describe the heterogeneity of both P. falciparum and P. vivax exposure in the study area. These results supported previous utilization of sero-epidemiological analysis in assessing population–level transmission intensity and differentiating between areas of different endemicity in Indonesia [30]. Moreover, multivariable analysis utilizing serological and epidemiological data collected through community-based survey identified that age, P. vivax seropositivity and LLIN use were significantly associated with P. falciparum seropositivity. These associations are likely related to historical exposure as P. falciparum seroprevalence was estimated to be low and parasite screening found no active infections detected by microscopy. Although sub-microscopic infections might present in the community, a previous study suggested that the proportion of sub-microscopic infections detected via PCR (polymerase chain reaction) was very low 0.07% (11/16,229) in the region [34].

The P. falciparum SCR estimates suggest that there was no exposure seen in children under 5 years old in both sub-districts in Sabang municipality. These results could represent a step change in P. falciparum transmission due to the successful impact of malaria control programme implemented in the study area, evidenced by lower antibody prevalence in children born after the intervention scale-up. This finding was supported by a previous study documenting a significant drop in malaria cases after the launch of the control program in 2004. Malaria cases in Sabang declined from 88 cases per 1000 population in 2004 to 1 per 1000 by 2010. The decline in malaria transmission in Sabang is likely related to an extensive IRS programme immediately following the tsunami in 2004, large scale LLIN distribution, and a change in malaria treatment policy to artemisinin-based combined therapy as first-line treatment for uncomplicated malaria [34]. Sabang was certified as a malaria-free region by the Indonesian Government as a result of successfully maintaining zero cases since the last locally transmitted case reported in 2011. Since then, the surveillance system detected 12 imported cases consisting of 6 P. vivax, 4 P. falciparum and 2 mixed P. vivax and P. falciparum infections from 2011 to 2013, with no local transmission. However, the surveillance system detected 15 PCR confirmed Plasmodium knowlesi infections that classified as an outbreak in 2014 [38].

Consistent with the higher P. falciparum SCR estimates in people over 5 years old, multivariable analysis revealed that adults were more likely to be seropositive compared to children under 15 years old. This is likely the result of higher exposure by staying overnight in high-risk areas. A recent study revealed that the clusters of malaria (P. knowlesi) infections in Sabang was associated with people who had a history of staying overnight in the forest, without protection from mosquitoes, in an area where macaques are common [38]. Unfortunately, data on travel behaviour and occupation in these surveys were not recorded to enable testing of these hypotheses. Future research would need to include more detailed questions regarding travel behaviour, occupation and other essential risk factor data such as travel history to high-risk areas, night outdoor activities, sleeping in plantation or forest, housing, personal protection, etc. Several programme initiatives, for example a multi-country study on vector control tools to address outdoor transmission and project management quality improvement for national malaria program workforce carried out under the Asia Pacific Malaria Elimination Network would be beneficial for the malaria elimination effort in the region. In addition, the use of LLIN was almost two times higher in area where P. falciparum seroprevalence was higher. Consistent with previous report suggesting high coverage of LLINs (over 75%) in six malaria focal villages in Sabang, this finding suggests that people living in higher risk of exposure were aware of the importance of LLIN to prevent malaria transmission in those areas [34].

The estimated age-seroprevalence curves and SCR value suggested that age was not associated with P. vivax transmission in either sub-district in Sabang. Plasmodium vivax seroprevalence was very low (2.0%) and, therefore, the absence of any associations is likely due to the statistical limitations of the low number of seropositive samples. The other possible explanation is that P. vivax infections may induce lower antibody responses or shorter-lived responses which the current assay may miss. Work is ongoing to identify P. vivax antigens that elicit short-term responses for easy identification of very recent exposure [39, 40]. The need for testing more potential P. vivax antigens is supported by a previous study showing that the number of P. vivax cases tend to be higher than the number of P. falciparum cases in Sabang [34].

The spatial analysis of age-adjusted antibody responses to either antigen (AMA-1 or MSP-119) identified significant clusters of higher exposure (hotspots) for both P. falciparum and P. vivax exposure across the study areas. Although multivariable risk factors analysis found there was no significant association between residence and higher seroprevalence to P. falciparum and P. vivax, the spatial analysis suggested that the risk of malaria transmission in the study setting is heterogeneous with people experiencing higher exposure in Sukajaya sub-district. The spatial analysis also suggest that the clusters identified for P. falciparum and P. vivax were seen in the same areas. Being able to characterize the micro-epidemiology of malaria exposure could assist malaria control programme to better allocate resources and target the intervention to achieve their goal of elimination. Targeting hotspots could be a highly efficient way to reduce malaria transmission at all levels of transmission intensity [41]. Although this study identified potential high-risk areas using historical data collected in 2013, being able to identify areas which had the most recent exposure is useful for malaria surveillance. A recent study suggested that one of two clusters of P. knowlesi infections in Sabang were identified in similar high-risk areas identified in this study [38]. As suggested in the latest WHO malaria surveillance manual [11], maintaining surveillance activities in the most receptive areas could be useful to prevent potential reintroduction or resurgence of the disease in the future. Therefore, utilizing antibody responses data to identify recent or historical hotspots of transmission could be a powerful alternative approach where gaining direct evidence of an increased exposure to infectious mosquito bites is no longer ideal in low transmission settings.

Finally, people who were seropositive to any P. vixax antigen were 3 times more likely to be P. falciparum seropositive, after controlling for age, gender, residence, employment, education, IRS, fever status, and altitude. In addition, clusters of high antibody responses suggest that P. falciparum and P. vivax receptive areas were seen in the same areas. As there was no cross-reactivity evident from the serological data (Additional file 3), these findings could suggest that people were historically exposed to both infections, potentially due to the presence of efficient vectors in those identified areas.

Findings in this study are based on community-based samples and data collected during the malaria transmission season. Although this study describes the potential use of serological data analysis in estimating malaria transmission intensity, heterogeneity and factors associated to disease exposure, the results generated would need to be carefully interpreted. Previous studies suggested that malaria transmission in other areas of Indonesia was affected by seasonality [30, 34, 42,43,44] and behavioural factors such as farm or forest-related night outdoor activity (e.g. sleeping in forest gardens) [45, 46] and domestic travel to higher endemic areas [47]. However, due to limited data collected, our study could not examine the effect of behavioural factors such as forest-related activities or recent travel history to high-risk areas outside Sabang. Therefore, future studies measuring population level antibody responses coupled with collecting more data that could describe behavioural factors associated to higher risk of exposure would be more epidemiologically informative to assist malaria surveillance and control programme to achieve elimination in the region.

Conclusion

In conclusion, these data add to the body of evidence that sero-epidemiological analysis of community-based surveys are an important additional tool to investigate malaria transmission dynamics in area aiming for elimination in Indonesia. Recent identification of alternative antigens associated with short-lived antibody responses suggests a potentially key indicator of very recent exposure which would be a very important information for public health surveillance [48]. The addition of a novel panel of P. knowlesi antigens [49] would enhance understanding of malaria transmission dynamics as recent studies reported that although laboratory identification of P. knowlesi in Indonesia is challenging [50], surprisingly, there were two clusters of P. knowlesi cases detected in Sabang after the municipality successfully eliminated P. falciparum and P. vivax cases [38]. Moreover, another recent study also reported there was a considerable proportion of P. knowlesi infection in another western part of Indonesia, in North Sumatera province [51]. Exploratory work employing techniques such as multiplex fluorescent magnetic bead-based serological assay to investigate and validate a panel of potential antigens for these applications is underway [40, 52]. The development and validation of a standardized serological sample and data collection methods utilizing existing public health surveillance system, for example as described in [53] will also facilitate the optimization of serological surveillance in understanding transmission dynamics to support malaria control programme in achieving elimination.

Availability of data and materials

The datasets used and analysed during this study are not publicly available due to the inclusion of identifying information on individuals but are available from the corresponding author on reasonable request and approval from relevant ethics committees.

Abbreviations

CI:

confidence interval

EIR:

entomological inoculation rate

IRS:

indoor residual spraying

LLIN:

long-lasting insecticide-treated bed nets

PfAMA-119 :

P. falciparum apical membrane antigen 1

PfMSP-1:

P. falciparum merozoite surface protein 1

PvAMA-119 :

P. vivax apical membrane antigen 1

PvMSP-1:

P. vivax merozoite surface protein 1

SCR:

seroconversion rates

WHO:

World Health Organization

References

  1. WHO. World malaria report 2018. Geneva: World Health Organization; 2018. p. 238.

    Google Scholar 

  2. Feachem RGA, Phillips AA, Hwang J, Cotter C, Wielgosz B, Greenwood BM, et al. Shrinking the malaria map: progress and prospects. Lancet. 2010;376:1566–78.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Kelly GC, Hale E, Donald W, Batarii W, Bugoro H, Nausien J, et al. A high-resolution geospatial surveillance-response system for malaria elimination in Solomon Islands and Vanuatu. Malar J. 2013;12:108.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Greenwood BM. Control to elimination: implications for malaria research. Trends Parasitol. 2008;24:449–54.

    Article  PubMed  Google Scholar 

  5. O’Meara WP, Mangeni JN, Steketee R, Greenwood B. Changes in the burden of malaria in sub-Saharan Africa. Lancet Infect Dis. 2010;10:545–55.

    Article  PubMed  Google Scholar 

  6. Tatem AJ, Smith DL, Gething PW, Kabaria CW, Snow RW, Hay SI. Ranking of elimination feasibility between malaria-endemic countries. Lancet. 2010;376:1579–91.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Ceesay SJ, Casals-Pascual C, Erskine J, Anya SE, Duah NO, Fulford AJ, et al. Changes in malaria indices between 1999 and 2007 in The Gambia: a retrospective analysis. Lancet. 2008;372:1545–54.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Okiro EA, Alegana VA, Noor AM, Snow RW. Changing malaria intervention coverage, transmission and hospitalization in Kenya. Malar J. 2010;9:285.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  Google Scholar 

  10. Drakeley CJ, Corran PH, Coleman PG, Tongren JE, McDonald SLR, Carneiro I, et al. Estimating medium- and long-term trends in malaria transmission by using serological markers of malaria exposure. Proc Natl Acad Sci USA. 2005;102:5108–13.

    Article  CAS  PubMed  Google Scholar 

  11. WHO. Malaria surveillance, monitoring and evaluation: a reference manual. Geneva: World Health Organization; 2018.

    Google Scholar 

  12. 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:1237.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Drakeley C, Cook J. Potential contribution of sero-epidemiological analysis for monitoring malaria control and elimination: historical and current perspectives. Adv Parasitol. 2009;69:299–352.

    Article  PubMed  Google Scholar 

  14. Cunha MG, Silva ES, Sepulveda N, Costa SPT, Saboia TC, Guerreiro JF, et al. Serologically defined variations in malaria endemicity in Para State, Brazil. PLoS One. 2014;9:e113357.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Zhang Q, Lai S, Zheng C, Zhang H, Zhou S, Hu W, et al. The epidemiology of Plasmodium vivax and Plasmodium falciparum malaria in China, 2004–2012: from intensified control to elimination. Malar J. 2014;13:419.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Pullan RL, Sturrock HJW, Soares Magalhães RJ, Clements ACA, Brooker SJ. Spatial parasite ecology and epidemiology: a review of methods and applications. Parasitology. 2012;139:1870–87.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  18. Wangdi K, Gatton ML, Kelly GC, Clements ACA. Cross-border malaria: a major obstacle for malaria elimination. Adv Parasitol. 2015;89:79–107.

    Article  PubMed  Google Scholar 

  19. Pindolia DK, Garcia AJ, Huang Z, Fik T, Smith DL, Tatem AJ. Quantifying cross-border movements and migrations for guiding the strategic planning of malaria control and elimination. Malar J. 2014;13:169.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Wesolowski A, Eagle N, Tatem AJ, Smith DL, Noor AM, Snow RW, et al. Quantifying the impact of human mobility on malaria. Science. 2013;338:267–70.

    Article  Google Scholar 

  21. Shekalaghe S, Alifrangis M, Mwanziva C, Enevold A, Mwakalinga S, Mkali H, et al. Low density parasitaemia, red blood cell polymorphisms and Plasmodium falciparum specific immune responses in a low endemic area in northern Tanzania. BMC Infect Dis. 2009;9:69.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Rolfes MA, McCarra M, Magak NG, Ernst KC, Dent AE, Lindblade KA, et al. Development of clinical immunity to malaria in highland areas of low and unstable transmission. Am J Trop Med Hyg. 2012;87:806–12.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Proietti C, Verra F, Bretscher MT, Stone W, Kanoi BN, Balikagala B, et al. Influence of infection on malaria-specific antibody dynamics in a cohort exposed to intense malaria transmission in northern Uganda. Parasite Immunol. 2013;35:164–73.

    Article  CAS  PubMed  Google Scholar 

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

  25. 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:1–9.

    Google Scholar 

  26. Tusting L, Bousema T, Smith D, Drakeley C. Measuring changes in Plasmodium falciparum transmission: precision, accuracy and costs of metrics. Adv Parasitol. 2014;84:151–208.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Bousema T, Okell L, Felger I, Drakeley C. Asymptomatic malaria infections: detectability, transmissibility and public health relevance. Nat Rev Microbiol. 2014;12:833–40.

    Article  CAS  PubMed  Google Scholar 

  28. Kilama M, Smith DL, Hutchinson R, Kigozi R, Yeka A, Lavoy G, et al. Estimating the annual entomological inoculation rate for Plasmodium falciparum transmitted by Anopheles gambiae s.l. using three sampling methods in three sites in Uganda. Malar J. 2014;13:111.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Bretscher MT, Supargiyono S, Wijayanti MA, Nugraheni D, Widyastuti AN, Lobo NF, et al. Measurement of Plasmodium falciparum transmission intensity using serological cohort data from Indonesian schoolchildren. Malar J. 2013;12:21.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Supargiyono S, Bretscher MT, Wijayanti MA, Sutanto I, Nugraheni D, Rozqie R, et al. Seasonal changes in the antibody responses against Plasmodium falciparum merozoite surface antigens in areas of differing malaria endemicity in Indonesia. Malar J. 2013;12:444.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  32. Bejon P, Turner L, Lavstsen T, Cham G, Olotu A, Drakeley CJ, et al. Serological evidence of discrete spatial clusters of Plasmodium falciparum parasites. PLoS One. 2011;6:e21711.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Mosha JF, Sturrock HJW, Greenwood B, Sutherland CJ, Gadalla NB, Atwal S, et al. Hot spot or not: a comparison of spatial statistical methods to predict prospective malaria infections. Malar J. 2014;13:53.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Herdiana H, Fuad A, Asih PB, Zubaedah S, Arisanti RR, Syafruddin D, et al. Progress towards malaria elimination in Sabang Municipality, Aceh, Indonesia. Malar J. 2013;12:42.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Sepúlveda N, Drakeley C. Sample size determination for estimating antibody seroconversion rate under stable malaria transmission intensity. Malar J. 2015;14:141.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  37. Biggs J, Raman J, Cook J, Hlongwana K, Drakeley C, Morris N, et al. Serology reveals heterogeneity of Plasmodium falciparum transmission in northeastern South Africa: implications for malaria elimination. Malar J. 2017;16:1–13.

    Article  Google Scholar 

  38. Herdiana H, Irnawati I, Coutrier FN, Munthe A, Mardiati M, Yuniarti T, et al. Two clusters of Plasmodium knowlesi cases in a malaria elimination area, Sabang Municipality, Aceh, Indonesia. Malar J. 2018;17:1–10.

    Article  Google Scholar 

  39. França CT, He W-Q, Gruszczyk J, Lim NTY, Lin E, Kiniboro B, et al. Plasmodium vivax reticulocyte binding proteins are key targets of naturally acquired immunity in Young Papua New Guinean Children. PLoS Negl Trop Dis. 2016;1:1–17.

    Google Scholar 

  40. França CT, Hostetler JB, Sharma S, White MT, Lin E, Kiniboro B, et al. An antibody screen of a Plasmodium vivax antigen library identifies novel merozoite proteins associated with clinical protection. PLoS Negl Trop Dis. 2016;10:e0004639.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  42. Syafruddin D, Krisin, Asih P, Sekartuti, Dewi RM, Coutrier F, et al. Seasonal prevalence of malaria in West Sumba district, Indonesia. Malar J. 2009;8:8.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Sutanto I, Freisleben H, Pribadi W, Atmosoedjono S, Bandi R, Purnomo. Efficacy of premethrin-impregnated bed nets on malaria control in hyperendemic area in Irian Jaya, Indonesia: influence of seasonal rainfall fluctuations. Southeast Asian J Trop Med Public Health. 1999;30:432–9.

    CAS  PubMed  Google Scholar 

  44. Ndoen EML. Environmental factors and an eco-epidemiological model of malaria in Indonesia. Griffith: Griffith University; 2009.

    Google Scholar 

  45. Dale P, Sipe N, Anto S, Hutajulu B, Ndoen E, Papayungan M, et al. Malaria in Indonesia: a summary of recent research into its environmental relationships. Southeast Asian J Trop Med Public Health. 2005;36:1–13.

    PubMed  Google Scholar 

  46. Herdiana H, Cotter C, Coutrier FN, Zarlinda I, Zelman BW, Tirta YK, et al. Malaria risk factor assessment using active and passive surveillance data from Aceh Besar, Indonesia, a low endemic, malaria elimination setting with Plasmodium knowlesi, Plasmodium vivax, and Plasmodium falciparum. Malar J. 2016;15:468.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Baird JK, Basri H, Weina P, MaGuire JD, Barcus MJ, Picarema H, et al. Adult Javanese migrants to Indonesian Papua at high risk of severe disease caused by malaria. Epidemiol Infect. 2003;131:791–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  49. Herman LS, Fornace K, Phelan J, Grigg MJ, Anstey NM, William T, et al. Identification and validation of a novel panel of Plasmodium knowlesi biomarkers of serological exposure. PLoS Negl Trop Dis. 2018;12:1–21.

    Article  Google Scholar 

  50. Coutrier FN, Tirta YK, Cotter C, Zarlinda I, Gonzalez IJ, Schwartz A, et al. Laboratory challenges of Plasmodium species identification in Aceh Province, Indonesia, a malaria elimination setting with newly discovered P. knowlesi. PLoS Negl Trop Dis. 2018;12:e0006924.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Lubis IND, Wijaya H, Lubis M, Lubis CP, Divis PCS, Beshir KB, et al. Contribution of Plasmodium knowlesi to multispecies human malaria infections in North Sumatera, Indonesia. J Infect Dis. 2017;215:1148–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Koffi D, Touré AO, Varela M-L, Vigan-Womas I, Béourou S, Brou S, et al. Analysis of antibody profiles in symptomatic malaria in three sentinel sites of Ivory Coast by using multiplex, fluorescent, magnetic, bead-based serological assay (MAGPIX™). Malar J. 2015;14:509.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Fornace KM, Surendra H, Abidin TR, Reyes R, Macalinao MLM, Stresman G, et al. Use of mobile technology—based participatory mapping approaches to geolocate health facility attendees for disease surveillance in low resource settings. Int J Health Geogr. 2018;17:21.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank the staff of the Sabang Municipal Health Office and all Sabang’s primary health centres for supporting data collection activities. We would also like to thank the project staff of and project participants for contributing to this work.

Funding

The authors acknowledge the Bill and Melinda Gate Foundation for funding received for the Indonesian Malaria Transmissium Consortium (MTC), the Indonesia Endowment Fund for Education (Lembaga Pengelola Dana Pendidikan, Indonesia) for funding received for OPSIN: Optimising Serological Surveillance for Malaria in Indonesia. The funders had no role in the design of the study, collection, analysis and interpretation of data or writing the manuscript. Henry Surendra is funded by the Indonesia Endowment Fund for Education (Lembaga Pengelola Dana Pendidikan Indonesia) (20151022084537).

Author information

Authors and Affiliations

Authors

Contributions

HS collated and analysed the data and wrote the manuscript; MAW contributed to study design, data collection and reviewed the manuscript; EHM contributed to blood samples collection and supervised the serological assay, MES, NFL and WAH contributed to the study design and field work supervision; I, TY, M and H carried out field work; JC and CD provided the antigens, standardized serological assay protocols, interpreted findings and reviewed the manuscript, S designed the study, coordinated field data collection, supervised serological assay and reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Henry Surendra.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board of the Medical Health Research Ethics Committee, Gadjah Mada University, Indonesia (KE/FK/927/EC) and the Research Ethics Committee of the London School of Hygiene and Tropical Medicine (11944). Written informed consent was obtained from all participants in this study.

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.

Additional files

Additional file 1.

Seroprevalence for each antigen studied.

Additional file 2.

Demographic characteristics and factors associated with P. vivax transmission in Sabang, Indonesia, 2013.

Additional file 3.

Scatter plots matrix of antibody responses (optical density) to P. falciparum and to P. vivax antigens tested in the study describing the absence of cross-reactivity between the P. falciparum and P. vivax antigens.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Surendra, H., Wijayanti, M.A., Murhandarwati, E.H. et al. Analysis of serological data to investigate heterogeneity of malaria transmission: a community-based cross-sectional study in an area conducting elimination in Indonesia. Malar J 18, 227 (2019). https://doi.org/10.1186/s12936-019-2866-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-019-2866-z

Keywords