Research | Open | Published:
Geographical patterns of malaria transmission based on serological markers for falciparum and vivax malaria in Ratanakiri, Cambodia
Malaria Journalvolume 15, Article number: 510 (2016)
Malaria transmission is highly heterogeneous, especially in low endemic countries, such as Cambodia. This results in geographical clusters of residual transmission in the dry, low transmission season, which can fuel the transmission to wider areas or populations during the wet season. A better understanding of spatial clustering of malaria can lead to a more efficient, targeted strategy to reduce malaria transmission. This study aims to evaluate the potential of the use of serological markers to define spatial patterns in malaria exposure.
Blood samples collected in a community-based randomized trial performed in 98 high endemic communities in Ratanakiri province, north-eastern Cambodia, were screened with a multiplex serological assay for five serological markers (three Plasmodium falciparum and two Plasmodium vivax). The antibody half-lives range from approximately six months until more than two years. Geographical heterogeneity in malaria transmission was examined using a spatial scan statistic on serology, PCR prevalence and malaria incidence rate data. Furthermore, to identify behavioural patterns or intrinsic factors associated with malaria exposure (antibody levels), risk factor analyses were performed by using multivariable random effect logistic regression models. The serological outcomes were then compared to PCR prevalence and malaria incidence data.
A total of 6502 samples from two surveys were screened in an area where the average parasite prevalence estimated by PCR among the selected villages is 3.4 %. High-risk malaria pockets were observed adjacent to the ‘Tonle San River’ and neighbouring Vietnam for all three sets of data (serology, PCR prevalence and malaria incidence rates). The main risk factors for all P. falciparum antigens and P. vivax MSP1.19 are age, ethnicity and staying overnight at the plot hut.
It is possible to identify similar malaria pockets of higher malaria transmission together with the potential risk factors by using serology instead of PCR prevalence or malaria incidence data. In north-eastern Cambodia, the serological markers show that malaria transmission occurs mainly in adults staying overnight in plot huts in the field. Pf.GLURP.R2 showed a shrinking pocket of malaria transmission over time, and Pf.MSP1.19, CSP, PvAMA1 were also informative for current infection to a lesser extent. Therefore, serology could contribute in future research. However, further in-depth research in selecting the best combination of antigens is required.
The annual malaria incidence and mortality have steadily declined in the Greater Mekong Subregion over the past 15 years . In particular, in the Kingdom of Cambodia malaria cases were reduced by more than 75 %, due to the improvements in malaria control by the National Malaria Control Programme, such as free distribution of long-lasting insecticidal nets (LLINs), performing better case management (i.e. diagnosis by RDT and treatment with ACT) [2, 3]. Although malaria transmission in Cambodia nowadays is low [4, 5], still 21 out of 25 provinces are endemic, of which the northeast region (Ratanakiri) accounts for more than 70 % of the malaria burden [3, 6]. Moreover, a shift towards more heterogeneous malaria transmission has been observed. This results in areas that support malaria transmission, which are referred to as foci [7, 8]. Within these foci, elevations in malaria transmission in small areas (sometimes <1 km2 ) or populations can be identified, which are respectively called hotspots or hotpops, presenting a higher risk of infection as compared to the rest of the focus [7, 9]. Some studies have shown that stable hotspots with a permanent transmission of the parasite over consecutive dry seasons mainly consist of asymptomatic carriers . Consequently, in low endemic areas as Cambodia, in the wet high transmission season, when the malaria vector population expands, these remaining reservoirs tend to fuel malaria transmission to surrounding areas or populations [7, 11]. Despite considerable malaria control efforts in Cambodia, this persisting transmission makes it impossible to eliminate malaria by current control strategies. Therefore, it is necessary to move forward from these strategies by creating a novel strategic plan based on the understanding of the hotspots/hotpops [3, 12]. It is assumed that there will be no persistence of malaria transmission once these recurrent sources of infection are eliminated by targeted interventions to these hotspots and hotpops [3, 8].
So far, hotspot identification has been carried out by various approaches, including microscopy, rapid diagnostic tests (RDTs), PCR detection and entomological surveys. Microscopy and RDTs cannot detect low-density infections . In addition, among the PCR-based assays, real-time PCR is a highly sensitive technique capable of detecting higher numbers of infected individuals, including asymptomatic infections. However, this technique has a couple of drawbacks when used in the field e.g. high costs and complexity of its applicability . Serological markers of malaria exposure, specifically antibodies (Abs) against Plasmodium antigens (Ags), are appropriate to use when detecting stable hotspots of malaria transmission in low endemic areas . These Ab-responses increase by cumulative exposure and the longevity of the Abs depends on the Ag . Therefore, this method can provide an indication of past and recent malaria exposure that can be used to pick up temporal and spatial trends in malaria transmission [7, 14]. Moreover, previous elimination programs have already observed that the absence of Ab-titers in the youngest age-groups could be used as proof of the cessation of malaria transmission . However, in Southeast Asia, primo infection by malaria parasites may be delayed to adolescence, due to behavioural and occupational activities . Serology has already been used to detect spatial trends by previous studies in high endemic settings [8, 11, 16–20]. On the contrary, the serological value for detecting spatial clustering of malaria exposure in low endemic areas has not yet been completely confirmed [7, 8].
The proposed study aims to further validate five Plasmodium markers for their potential to detect recent infection  by defining spatial patterns in malaria exposure over two different surveys in comparison with PCR prevalence and malaria incidence data. As in this study, community (= cluster) based data were used, the outcomes were defined as ‘malaria pockets’  referring to an area in between a hotspot (<1 km2 ) and a foci (>1 province) where the malaria exposure is higher than the surrounding areas.
Ratanakiri province (13°44′N, 107°0′E) bordering Lao PDR and Vietnam, is located 520 km from Phnom Penh in the northeast of Cambodia. The area has a monsoonal climate, resulting in perennial malaria transmission with a peak during the rainy season (April until October) [5, 15]. Ratanakiri has a population of 149,997 individuals spread over 240 villages , of which approximately 70 % is living in the highlands and 30 % in the urbanized towns. This study area is largely inhabited by the ethnic minority (e.g. Jarai, Kreung, Tumpuon) as opposed to the Khmer in the rest of the Country. The ethnic minority generate revenue by subsistence slash-and-burn farming. Therefore, they own plot huts located near or inside the forests, where they stay during most of the rainy season [23, 24].
Samples used were derived from a community (=cluster) based randomized trial (MalaResT project—NCT01663831) that aims to evaluate the effectiveness of topical repellents, in addition to long-lasting insecticidal nets, on malaria prevalence and incidence [5, 21, 25]). For the purpose of this study two cross-sectional surveys carried out in November 2012 and 2013 were included. No differences were observed between the control and the intervention arm for PCR prevalence, serological indicators and malaria incidence. Blood samples were collected by a finger prick on filter paper, and immediately screened by real-time PCR to determine the presence or absence of parasites . In the MalaResT project 65 people per community were randomly selected, and in case of a low attendance rate an additional set of 15 randomly selected people was added to reach at least 50 participants per community. Throughout the entire sampling process an over 70 % success rate was reached, in which all age groups were proportionally covered . From a total population of 48,838 individuals residing in all 113 villages grouped in 98 communities (i.e. clusters; 88 single villages +25 neighbouring villages having a distance of <1 km were grouped into 10 communities), a total of 6640 and 6715 were randomly recruited by community in November 2012 and 2013 of which respectively 4996 and 5431 were sampled . Serology was performed on 3264 (2012) and 3238 (2013) on randomly selected samples. The PCR Plasmodium prevalence was 4.9 and 3.4 % in November 2012 and 2013 respectively (Table 1) .
Initially, each Ag was coupled to paramagnetic beads (MagPlex microspheres, Luminex Corp, Austin, TX, USA) as described earlier . All beads with different Ags were put together to prepare a microsphere working mixture at a concentration of 1000 beads/Ag/well. Bloodspot samples were then analysed in 96-well plates in duplicate. Positive controls (pool of four P. falciparum and two P. vivax infected individuals), negative control serum and blanco (PBS-CR) were added in duplicate to each 96-well plate . The MAGPIX®-system was set for reading a minimum of 400 beads per spectral address and results were expressed as median fluorescent intensity (MFI) .
The selection of the five Ags used in this study was based on the half-lives estimated on 20 different Ags previously . For P. falciparum three out of the six serological markers that were most likely to be reflective for recent exposure were chosen. These were Pf.GLURP.R2, Pf.MSP1.19 and CSP, showing a half-life of respectively six months, ~8 months and ~1 year. For P. vivax, the shortest Ab half-life found was more than 1.5 years (PvAMA1 and Pv.MSP1.19).
Passive case detection (confirmed malaria cases) was reinforced for the purpose of the study and rely on the National Health System (Village Malaria Workers, health centres and hospitals) . For each case the living place was recorded in the database. Incidence was estimated by community.
The serology outcome data (median fluorescent intensity—MFI) were processed and analysed in R-version 3.1.0 . Spatial clustering of serological data was detected for each survey by using spatial scan statistic (SaTScan) [28, 29], which is the preferred software in low risk settings according to Aamodt et al. [25, 30]. The spatial scan statistic was used on the natural logarithm transformed median fluorescent intensity values (ln(MFI)—normal probability model) to define spatial patterns in malaria exposure and to compare these outcomes to malaria pockets obtained by the PCR prevalence (binominal, Bernoulli model) and malaria incidence data recorded by community in all ages (counts, Poisson model).
The R package ‘rsatscan’  was used to prepare the data, and R package ‘PlotKML’  was used to be able to plot the 98 communities on the Google Earth images created with SaTScan. By using the spatial scan statistic, SaTScan version 9.4.2, 64-bit [28, 29], the study area (province of Ratanakiri) was systematically screened for circular windows of higher MFI values. Ratanakiri is a large area (10,782 km2). To avoid detection of too large pockets that cover almost the entire area and hide small homogeneous malaria pockets within the larger pockets, the maximum of the total population at risk within a community was set at 20 %, in accordance with the observations in the study of Mosha et al. . This allows that both small and large malaria pockets can be detected . Observed and expected means of ln (MFI) values inside and outside each window at each location were calculated [8, 11, 29]. The areas with the maximum likelihood were defined as the malaria pockets . These pockets were then examined based on 999 Monte Carlo simulations. The malaria pockets were considered statistically significant with a p value <0.05. Circular windows were visualized separately for November 2012 and 2013. The median age inside and outside the malaria pockets was estimated for the serological data.
The outcomes perceived with the serological markers were compared to PCR prevalence (performed on the same surveys) and malaria incidence rate (recorded the same year) . When these outcomes correspond to the results observed in this study, ongoing malaria transmission at the localities could be confirmed. Therefore, these malaria pockets were detected by using a similar approach as detailed above, but with the PCR prevalence data (Bernoulli model) from November 2012 and 2013, whereby the maximum allowed population size was set to 20 % as well. Thereby, pockets were significant with a p-value below 0.05. The same was done for the malaria incidence data of symptomatic cases based on the health information system (Poisson model) from 2012 and 2013. Pockets were considered significant when presenting a p-value below 0.05. A sensitivity and specificity analysis was done to assess the prediction of serological based malaria pockets in comparison to PCR or Incidence malaria pockets.
Moreover, a risk factor analysis was performed to detect a pattern in behavioral or intrinsic factors associated with high endemic pockets. The risk factor analysis was carried out in various steps as described previously , with minor adjustments. Firstly, a univariate analysis was performed on all explanatory variables: age, gender, ethnicity (khmer vs ethnic minority), axillary temperature, plot hut (a human behaviour associated with indigenous farming), sleeping in the forest (indoor/outdoor) and repellent use [control (only bed net) vs intervention (bed net and repellent use)] . The risk factors were analysed by fitting linear mixed effect models (lmer function in the ‘lme4’ package applied in R version 3.1.0 ) with ln(MFI) values as outcome variable taking into account community (=cluster) within survey (year) as random effect . Incidence rate ratios (IRR) and 95 % confidence intervals (95 % CI) were estimated by exponentiation of the model coefficients per variable and per Ag. Statistical significance was evaluated based on the p-values below 0.10. Next, residual plots were made to evaluate how well the models fitted the data and how the data meet the assumptions of the model. To check whether variables needed to be omitted due to co-linearity for the multivariable model, the Variance inflations factors (VIF) of each linear model were evaluated (R packages ‘MASS’ and ‘car’ [36, 37]). Lastly, multivariable random effect logistic regression models were fitted (lmer function in R package ‘lme4’ ). Model selection through stepwise backward deletion, starting from the full model, was built based on the outcomes of the univariate analysis to define independent relations between the variables and the strength of the Ab-response [38, 39].
For this study a random selection of 6502 samples were screened with the Luminex technology after being tested for parasite infection using real-time PCR. An overview of the descriptive statistics is given in Table 1.
Geographical clustering based on serology (Plasmodium falciparum and Plasmodium vivax antigens)
A first pocket was detected around the most northerly point of the ‘Tonle San River’ for all five Ags (P. falciparum and P. vivax) in November 2012 (Fig. 1; Table 2). For four out of five Ags (CSP, Pf.MSP1.19, Pf.GLURP.R2 and Pv.MSP1.19), the mean geographical radius was estimated 21.33 km, consisting of approximately 18 communities (p = 0.001). Moreover, for PvAMA1 one pocket was located at the same place, but showed a smaller radius of 16.45 km, with 8 communities (p = 0.001). In the same year, another pocket (Fig. 1; Table 2) was perceived around the ‘Tonle San River’ nearby the border with Vietnam for Pf.MSP1.19, Pf.GLURP.R2 and PvAMA1, with a radius of respectively 5.97 km (5 communities), 2.90 km (2 communities) and 10.67 km (7 communities) (p = 0.001). For the remaining two Ags (CSP and PvMSP1.19), a malaria pocket was detected adjacent to Vietnam as well, but was distant from the ‘Tonle San River’. These malaria pockets had a radius of 0 km, as these consist of only one community (p = 0.002). PvAMA1 showed a fourth pocket more inland (3 communities, radius 7 km), and for Pf.MSP1.19 a pocket close to the capital of Ratanakiri, ‘Ban Lung’, was observed (1 community, radius 0 km).
In November 2014, a malaria pocket was observed again for both Ags Pf.GLURP.R2 and Pv.AMA1, located at the northerly point of the ‘Tonle San River’ (Fig. 1; Table 2). For Pf.GLURP.R2, the pocket size was reduced by approximately 35 % (from 18 to 6 communities in the malaria pockets) as compared to November 2012, whereas for PvAMA1 the radius was reduced by 8 % (from 8 to 6 communities in the malaria pockets), however, for PvAMA1 the communities neighbouring Vietnam increased with 16.8 % (from 7 to 12 communities in the malaria pockets) and showed a more northern spread compared to November 2012. A fourth pocket for Pf.GLURP.R2 was perceived closer to ‘Ban Lung’ (5 communities, radius 7.54 km).
A different pattern was perceived for CSP, showing an eastward shift (37.35 km) in the malaria transmission locality between November 2012 and 2013 (Fig. 1; Table 2). This resulted in only one large malaria pocket (>1 community) in November 2013, neighbouring Vietnam and having a radius of 12.46 km consisting of 12 communities (p = 0.001).
Primary malaria pockets obtained for Pf.MSP1.19 and Pv.MSP1.19 were similar in November 2012 and 2013 (Fig. 1; Table 2) having a mean radius of 23.18 km, and 19 communities (p < 0.001) with a shift of the main pocket towards the Vietnamese border in November 2013 (mean radius of 27.49 km, 19 communities (p < 0.001). Additionally, in November 2013 a smaller malaria pocket appeared for both Ags with a mean radius of 6.44 km, including 6 communities (p < 0.001) more to the northwest side of the ‘Tonle San River’.
For all Ags the age distribution inside and outside the pockets was shown to be similar (Table 2).
Confirmation of the ongoing malaria transmission
Plasmodium prevalence and malaria incidence were determined by PCR (prevalence per survey, Fig. 2; Table 3)  and passive case detection (yearly incidence rate, Fig. 2; Table 4). The same pockets as the ones detected by serological markers were observed with PCR in the most northerly area of the ‘Tonle San River’ in November 2012 for P. falciparum, whereas no PCR prevalence pockets were detected in November 2013. Looking at the P. falciparum malaria incidence data, malaria pockets did not move between 2012 and 2013, but slightly shifted to the western part of the river as compared to the main pocket acquired with PCR prevalence and serological markers.
For P. vivax sero-reactive pockets were located nearby the ‘Tonle San River’ and the border of Vietnam in November 2012. Malaria pockets detected by PCR prevalence data show a shift between the surveys from the river-side towards the south (inland) of Ratanakiri in November 2013, the same pocket was observed for PvAMA1 in November 2012. This shift to the south is also seen with the malaria incidence data but on a different locality namely at the border of Vietnam. However, in 2013 (malaria incidence data) the main pocket persisted at the northern site of the ‘Tonle San River’, which was not the case for PCR prevalence data.
Overall, malaria pockets found by analysing malaria incidence data and ln(MFI) data are similar. The largest serologically measured pocket of each Ag consistently overlaps with the largest incidence based pockets. The malaria pockets adjacent to Vietnam overlap for all three approaches. Moreover, serologically obtained pockets from November 2012 show overlap with both P. falciparum and P. vivax pockets found by PCR. In November 2013, no significant PCR prevalence based pockets were observed for P. falciparum and the largest PCR prevalence based pocket for P. vivax was located more inland.
Sensitivity and specificity analysis (Additional file 1: Table S1)
Serology with the three P. falciparum Ags provides a relative good specificity (>74 %) and sensibility (>72 %) to identify positive PCR falciparum pockets. This was particularly true in November 2012 with the Pf.GLURP.R2 Ag having 95 % for both specificity and sensitivity. As no PCR falciparum pockets could be identified for November 2014, sensitivity could not be calculated. To predict P. falciparum incidence pockets, serological data provide a much lower sensitivity (between 13 and 60 %) although specificity was around 80 %. Serological data with the two P. vivax Ags provides a low sensitivity (between 0 and 38 %) but a relatively high specificity (69 and 82 %) to predict PCR vivax pockets. Results were slightly better in predicting P. vivax incidence pockets (sensitivity 10 and 57 %; specificity around 80 %). For comparison P. falciparum and P. vivax incidence data to predict PCR falciparum and vivax pockets sensitivity was around 50 % in 2012, but only 10 % in 2013 for P. vivax. Specificity lies between 70 and 80 %.
Risk factor analysis
Significant variables (age, gender, ethnicity, plot hut and sleeping in the forest) observed in the univariate analysis (Additional file 2: Table S2) were further explored in the multivariable regression model (Table 5). In the full model used for the multivariable analysis age, gender, ethnicity, plot hut and sleeping in the forest were included as dependent variables, and community within survey as random effect. Axillary temperature and repellent were not significant for all five Ags in the univariate exploration and were, therefore, not included.
Differences in the presence of Abs were perceived between age categories for all Ags. Ethnic minority showed a higher presence of Ab-levels for all Ags (mean IRR 2.051 95 % CI [1.788–2.354]). Significant differences were observed for ‘overnight stay in a plot hut’ (p < 0.10) with a mean IRR and 95 % CI of 1.126 [1.067–1.188]. Furthermore, gender appeared to be significantly different for CSP, Pf.MSP1.19 and Pv.MSP1.19, and sleeping in the forest for CSP and Pv.MSP1.19. However, the IRR of these last two risk factors are so close to one, that this is probably negligible.
Methods that can identify stable areas of transmission over time are suggested to be most effective for assessing geographical variations in malaria exposure. Therefore, Ab-responses acquired with cumulative malaria exposure, measured over several seasons, were recommended for implementation in geographical clustering analyses . During the ‘90s geographical cluster analyses mainly relied on symptomatic cases with accurate details about the place of infection and/or residence [20, 40], and were merely based on passive case detection (PCD) . It was not until the 21th century that due to a lack of information about the parasite reservoir in asymptomatic cases , new studies arose focussing on the spatial distribution estimated with PCR-prevalence data of species-specific geographical areas of infections based on asymptomatic carriers . This approach was especially important in countries with a low endemicity where the majority of infected people are asymptomatic carriers . Another innovative approach is the application of serological markers in defining these geographical areas. Serology already proved its ability to improve predictions of low transmission risk [26, 41].
Where most studies only focused on PfAMA1 and PfMSP1.19 [8, 11, 16–19], the advantage of this study is the amount of additional Ags from both falciparum and vivax malaria investigated compared to most other studies. Only one previous geostatistical study has used several Ab markers (namely PvAMA1, PvMSP1.19, PfAMA1 and Pf.GLURP.R2). However, in contrast with the current study, the researchers considered an individual positive when it responded for any of the two Ags for each species, not taking into account the differences in biological activity (e.g. longevity) among these Ags .
The previous study performed by Kerkhof et al.  has led to identification of serological markers with a relatively short half-life that were most likely to be reflective for recent exposure, such as P. falciparum Ags Pf.GLURP.R2, Pf.MSP1.19 and CSP. These serological markers could map the transmission risks with more precision and accuracy, as they provide the ability of distinguishing recent from past exposure [26, 41]. The current study, presented here, explored whether or not the use of serological markers is comparable to the use of PCR prevalence (asymptomatic cases) and malaria incidence data (symptomatic cases) to investigate spatial patterns in malaria transmission.
In Ratanakiri, significant malaria pockets were observed for both P. falciparum and P. vivax Ags. The largest pockets were located around the most northerly site of the ‘Tonle San River’ for all Ags. In comparison with the PCR prevalence data, P. falciparum exhibited similar pockets, whereas for P. vivax differences were seen. The similar pockets found between the PCR prevalence rates and sero-reactivity are in line with a study performed by Bousema et al. , that observed tight correlations as well.
When comparing the serologically based pockets with the incidence based pockets, the pockets neighbouring Vietnam were comparable, while the most northerly pockets at the ‘Tonle San River’ were slightly shifted to the West. There are malaria incidence based pockets found for P. falciparum situated around the capital ‘Ban Lung’ of the Ratanakiri province. Different studies [24, 42, 43] investigated the movement of individuals between villages, districts and countries. This might explain the malaria pockets seen around ‘Ban Lung’ raising the possibility that these individuals travel occasionally towards communities nearby the river or to remote areas. Overall, overlap was seen in the serological based pockets compared to the malaria incidence and parasite prevalence data.
That most pockets were perceived around the river confirms findings from other studies that also found more malaria pockets along open water bodies [19, 20, 44, 45]. The same pattern was observed by Sluydts et al.  who suggests that this is perhaps associated with increased movements of infected individuals and mosquito populations along the ‘Tonle San River’, and with the more remote location of these villages .
The specificity of serological markers for P. falciparum and P. vivax was acceptable (between 72 and 95 %) to predict PCR and malaria incidence pockets. However, sensitivity was in general much lower, except in predicting P. falciparum PCR pockets (between 74 and 95 %). In comparison, sensitivity of P. falciparum and P. vivax incidence in predicting PCR pockets lies between 10 and 50 %.
When looking at the different serological markers, variable patterns were observed, going from malaria pockets that move between the east and west in November 2012 and 2013 (CSP, Pf.MSP1.19 and Pv.MSP1.19) to lasting pockets that became smaller (Pf.GLURP.R2) or remained similar in size (PvAMA1). These varying patterns require further investigation related to the differences in immunogenicity and persistence of the Ab-responses [11, 46]. The only Ag that follows an expected altering pattern over time was Pf.GLURP.R2, which seems to correlate best with the PCR-prevalence and malaria incidence data. The latter is probably explained by the fact that this is a blood stage Ag with a short estimated Ab half-life . This might reflect recent exposure with observing pockets that decline over time, suggesting that this serological marker might have potential in evaluating targeted malaria control efforts.
Risk factors related to sero-reactivity were identified by univariate and multivariable analyses. Significant elevated risks for P. falciparum malaria were seen for age, ethnicity and overnight stay at the plot hut. There were also differences observed in gender and sleeping in the forest, however, this is most probably negligible, as the IRR was very close to 1. Significant elevated risks for vivax markers were seen for age, whereas staying in plot huts showed to be a risk factor for Pv.MSP1.19 only. The P. falciparum outcomes are in line with a previous study performed in the same area by Sluydts et al. . In this PCR prevalence based study that was performed on the baseline survey during the dry season, the most important risk factor detected was the overnight stay in the plot hut, based on both univariate and multivariable analyses. However, in the current serologically based study, it seems that age, concerning the older age groups, was the most important factor determining Ab-levels, compatible with cumulative exposure [7, 17, 26]. When immunity is acquired these Abs can persist for several years. This is caused by the presence of long-lived plasma cells that with every new exposure rapidly produce Abs against these parasites . Therefore, when defining current exposure it is important to observe the Ab-levels in especially the younger age groups [17, 47]. Differences between P. falciparum and P. vivax could be explained by that fact that P. vivax shows relapse patterns that influence the serological outcomes, and that longer half-lives were observed for the P. vivax Ags in a previous study .
Overall, these outcomes confirm the findings of Sluydts et al. , and are also in line with the findings of Incardona et al. . These researchers mentioned that entire families go together to the field and sleep in plot huts resulting in an increased exposure risk [25, 48]. Although this is not related to the age differences, as the age composition was similar inside and outside the pockets. However, this can be explained by the immunological maturity-status where children that acquire a malaria infection have the ability to boost their IgG titers, followed by a rapid decay . These outcomes explain the population characteristics in the Greater Mekong Subregion, where ethnic minority groups, forest workers (of all ages) and migrants are known as the most widely recognized groups at risk [25, 50].
This study contributes in the validation of serological markers to distinguish very recent from past exposure, as suggested by Sturrock et al. . This is especially the case for Pf.GLURP.R2, but also for Pf.MSP1.19, CSP and PvAMA1. By this means, combining more different Ags, covering the entire Plasmodium life cycle and having a longevity ranging from very short (~1 months) to long (year round), might lead to other promising results. However, methods to acquire the exact Ab-persistence are still in its infancy . While PvAMA1 showed stable malaria pockets, Pf.GLURP.R2 suggest a decline in the remaining malaria pockets. The stability of pockets was also observed by Mosha et al.  on P. falciparum Ag AMA1 in the high endemic setting of Tanzania. Further development in quantifying exposure over different timescales, as well as the measurement of very recent exposure, serological approaches will provide a major contribution in estimating spatio-temporal patterns of risk . The use of serology could benefit future malaria control programmes, since the use of serological markers can more precisely identify variation in transmission in low endemic areas. It should be noted that more serological markers that are competent to estimate exposure over different time-scales are required, as at present Pf.GLURP.R2 is most informative , as well as Pf.MSP1.19, CSP and PvAMA1 to a lesser extent. However, Pv.MSP1.19 should certainly not be ruled out, as it probably reflects transmission in the former past for which no PCR-prevalence data may be available.
Identification of pockets with higher malaria transmission would be essential when adopting a malaria elimination strategy. Present study shows that PCR parasite prevalence, malaria incidence rates or serology, show equivalent results in identifying malaria pockets. The attempt to validate serological markers that are most likely to reflect current and past exposure led to Pf.GLURP.R2 showing a shrinking malaria pocket over time. Moreover, Pf.MSP1.19, CSP and PvAMA1 are also reflective for recent malaria transmission to a lesser extent. This means that serology can provide promising information for future research, especially in evaluating short-term interventions for malaria elimination. However, there is still more research required in selecting a promising combination of Ags, in particular for P. vivax.
WHO: Strategy for malaria elimination in the Greater Mekong Subregion (2015–2030). Geneva: World Health Organization. 2015.
WHO: World Malaria Report 2015. Geneva: World Health Organization. 2015.
Ministry of Health. Cambodia malaria elimination action framework 2016–2020. Phnom Penh, 2016.
Imwong M, Nguyen TN, Tripura R, Peto TJ, Lee SJ, Lwin KM, et al. The epidemiology of subclinical malaria infections in South-East Asia: findings from cross-sectional surveys in Thailand-Myanmar border areas, Cambodia, and Vietnam. Malar J. 2015;14:381.
Canier L, Khim N, Kim S, Sluydts V, Heng S, Dourng D, et al. An innovative tool for moving malaria PCR detection of parasite reservoir into the field. Malar J. 2013;12:405.
Maude RJ, Nguon C, Ly P, Bunkea T, Ngor P, Canavati de la Torre SE, et al. Spatial and temporal epidemiology of clinical malaria in Cambodia 2004–2013. Malar J. 2014;13:385.
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.
Bousema T, Stevenson J, Baidjoe A, Stresman G, Griffin JT, Kleinschmidt I, et al. The impact of hotspot-targeted interventions on malaria transmission: study protocol for a cluster-randomized controlled trial. Trials. 2013;14:36.
Mosha JF, Sturrock HJW, Greenhouse B, Greenwood B, Sutherland CJ, Gadalla N, et al. Epidemiology of subpatent Plasmodium falciparum infection: implications for detection of hotspots with imperfect diagnostics. Malar J. 2013;12:221.
Bejon P, Williams TN, Liljander A, Noor AM, Wambua J, Ogada E, et al. Stable and unstable malaria hotspots in longitudinal cohort studies in Kenya. PLoS Med. 2010;7:e1000304.
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.
Cui L, Yan G, Sattabongkot J, Chen B, Cao Y, Fan Q, et al. Challenges and prospects for malaria elimination in the Greater Mekong Subregion. Acta Trop. 2012;121:240–5.
Kumar S, Kumari R, Pandey R. New insight-guided approaches to detect, cure, prevent and eliminate malaria. Protoplasma. 2015;252:717–53.
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.
Gryseels C, Durnez L, Gerrets R, Uk S, Suon S, Set S, et al. Re-imagining malaria: heterogeneity of human and mosquito behaviour in relation to residual malaria transmission in Cambodia. Malar J. 2015;14:165.
Bousema T, Drakeley C, Gesase S, Hashim R, Magesa S, Mosha F, et al. Identification of hot spots of malaria transmission for targeted malaria control. J Inf Dis. 2010;201:1764.
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:392.
Kangoye DT, Noor A, Midega J, Mwongeli J, Mkabili D, Mogeni P, et al. Malaria hotspots defined by clinical malaria, asymptomatic carriage, PCR and vector numbers in a low transmission area on the Kenyan Coast. Malar J. 2016;15:213.
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;6:1251.
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.
Kerkhof K, Canier L, Kim S, Heng S, Sochantha T, Sovannaroth S, et al. Implementation and application of a multiplex assay to detect malaria-specific antibodies: a promising tool for assessing malaria transmission in Southeast Asian pre-elimination areas. Malar J. 2015;14:338.
Gething PW, Elyazar IRF, Moyes CL, Smith DL, Battle KE, Guerra CA, et al. A long neglected world malaria map: Plasmodium vivax endemicity in 2010. PLoS Negl Trop Dis. 2012;6:e1814.
Sluydts V, Durnez L, Somony H, Gryseels C, Canier L, Kim S, et al. Efficacy of topical mosquito repellent (Picaridin) plus long-lasting insecticidal nets versus long-lasting insecticidal nets alone for control of malaria: a cluster randomised controlled trial. Lancet Infect Dis. 2016;16:1169–77.
Gryseels C, Grietens KP, Dierickx S, Bannister-Tyrrell M, Trienekens S, Xuan XN, et al. High mobility and low use of malaria preventive measures among the jarai male youth along the Cambodia-Vietnam border. Am J Trop Med Hyg. 2015;93:810.
Sluydts V, Heng S, Coosemans M, Van Roey K, Gryseels C, Canier L, et al. Spatial clustering and risk factors of malaria infections in Ratanakiri Province, Cambodia. Malar J. 2014;13:387.
Kerkhof K, Sluydts V, Willen L, Saorin K, Canier L, Heng S. Serological markers to measure recent changes in malaria transmission at population level in Cambodia. Malar J. 2016. Revised.
R Core Team. R: a language and environment for statistical computing. R Foundation for statistical computing. Vienna, Austria. 2015. https://www.r-project.org/. Accessed 1999.
Kulldorf M. SaTScan™. Software for the spatial, temporal, and space-time scan statistics. Boston: Harvard Medical School and Harvard Pilgrim Health Care; 2015. http://www.satscan.org/. Accessed 2005.
Kulldorf M. A spatial scan statistics. Commun Stat Theor M. 1997;26:1481.
Aamodt G, Samuelsen SO, Skrondal A. A simulation study of three methods for detecting disease clusters. Int J Health Geogr. 2006;5:15.
Kleinman K. rsatscan—tools, classes, and methods for interfacing with SaTScan stand-alone software. 2016. http://www.satscan.org/. Accessed 17 Mar 2015.
Heng T, Roudier P, Beaudette D, Pebesma E. Scientific visualization of spatio-temporal data. J Stat Softw. 2015;5:1.
Chen J, Roth RE, Naito AT, Lengerich EJ, Maceachren AM. Geovisual analytics to enhance spatial scan statistic interpretation: an analysis of US cervical cancer mortality. Int J Health Geogr. 2008;7:57.
Bates D, Martin M. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1.
Khaireh BA, Briolant S, Pascual A, Mokrane M, Machault V, Travaillé C, et al. Plasmodium vivax and Plasmodium falciparum infections in the Republic of Djibouti: evaluation of their prevalence and potential determinants. Malar J. 2012;11:395.
Fox J, Weisberg S. An R companion to applied regression. 2nd ed. Thousand Oaks: SAGE Publications, Inc; 2011.
Venabled W, Ripley B. Modern applied statistics with S. 4th ed. New York: Springer; 2002.
Kuznetsova A, Brockhoff B, Bojesen Christensen H. lmerTest. Tests in linear mixed effect models. 2016. https://cran.r-project.org/. Accessed 23 Jun 2016.
Souza-Silva FA, da Silva-Nunes M, Sanchez BAM, Ceravolo IP, Malafronte RS, Brito CFA, et al. Naturally acquired antibodies to Plasmodium vivax Duffy binding protein (DBP) in Brazilian Amazon. Am J Trop Med Hyg. 2010;185:193.
Gamage-Mendis AC, Carter R, Mendis C, De Zoysa APK, Herath PRJ, Mendis KN. Clustering of malaria infections within an endemic population: risk of malaria associated with the type of housing construction. Am J Trop Med Hyg. 2010;82:185.
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:1.
Guyant P, Canavati SE, Chea N, Ly P, Whittaker MA, Roca-Feltrer A, et al. Malaria and the mobile and migrant population in Cambodia: a population movement framework to inform strategies for malaria control and elimination. Malar J. 2015;14:252.
Parker DM, Matthews SA, Yan G, Zhou G, Lee MC, Sirichaisinthop J, et al. Microgeography and molecular epidemiology of malaria at the Thailand-Myanmar border in the malaria pre-elimination phase. Malar J. 2015;14:198.
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.
Brooker S, Clarke S, Njagi JK, Polack S, Mugo B, Estambale B, et al. Spatial clustering of malaria and associated risk factors during an epidemic in a highland area of Western Kenya. Trop Med Int Health. 2004;9:757.
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.
Hviid L, Barfod L, Fowkes FJI. Trying to remember: immunological B cell memory to malaria. Trends Parasitol. 2015;31:89.
Incardona S, Vong S, Chiv L, Lim P, Nhem S, Sem R, et al. Large-scale malaria survey in Cambodia: novel insights on species distribution and risk factors. Malar J. 2007;6:37.
Cook J, Speybroeck N, Sochanta T, Somony H, Sokny M, Claes F, et al. Sero-epidemiological evaluation of changes in Plasmodium falciparum and Plasmodium vivax transmission patterns over the rainy season in Cambodia. Malar J. 2012;11:86.
WHO: Malaria in the Greater Mekong Subregion: regional and country profiles. Geneva: World Health Organization. 2010.
Sample collection of the survey was performed by the CNM team (National Centre for Parasitology, Entomology and Malaria Control), the ITM Antwerp team and the Institut Pasteur du Cambodge (IPC) team in Ratanakiri Province in Cambodia. Screening of the blood spot samples was performed at the IPC by KK and LW and LC. The study design and literature research has been performed at the ITM by KK, MP, LD, VS and MC. Data entry and quality control was performed at the ITM by KK, LD and VS. Statistical analyses were performed by KK, LD and VS. The first draft of the manuscript was written by KK. All authors read and approved the final manuscript.
We would like to thank all colleagues at the ITM and IPC for all their support and contributions to fulfil this manuscript. We would like to thank Chris Drakeley for providing recombinant proteins PvMSP1-19 and PfMSP1-19 and Takafumi Tsuboi for providing recombinant protein PvAMA1.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets supporting the conclusions of this paper are included within the paper. Raw data may be obtained from the corresponding author on request.
Ethics approval and consent to participate
The study protocol was reviewed and approved by the Cambodian National Ethics Committee on Health Research (Approval 265 NECHR), the Institutional Review Board of the Institute of Tropical Medicine Antwerp (Approval IRB/AB/ac/154) and the Ethics Committee of the University of Antwerp (Approval B300201112714). Gatekeepers provided informed written consent for the participation of their village. The survey participant or his/her parents or guardian provided informed written consent for individual participation.
This work is part of a larger project “MalaResT” funded by the Bill & Melinda Gates Foundation under the Global Health Grant number OPP1032354. This project aims to evaluate the effectiveness of topical repellents, in addition to long-lasting insecticidal nets, on malaria prevalence and incidence. Next, this work fits within the large-scale research project for the SOFI research program “MalaSpot” entitled ‘Eco-epidemiology of asymptomatic malaria hotspots, a prerequisite for malaria elimination?’. Furthermore, KK was supported by the association of Les Amis des Instituts Pasteur à Bruxelles.