Skip to content

Advertisement

  • Methodology
  • Open Access

Estimating malaria transmission intensity from Plasmodium falciparum serological data using antibody density models

  • 1, 2, 3Email author,
  • 1,
  • 4 and
  • 1
Malaria Journal201615:79

https://doi.org/10.1186/s12936-016-1121-0

  • Received: 27 November 2015
  • Accepted: 22 January 2016
  • Published:

Abstract

Background

Serological data are increasingly being used to monitor malaria transmission intensity and have been demonstrated to be particularly useful in areas of low transmission where traditional measures such as EIR and parasite prevalence are limited. The seroconversion rate (SCR) is usually estimated using catalytic models in which the measured antibody levels are used to categorize individuals as seropositive or seronegative. One limitation of this approach is the requirement to impose a fixed cut-off to distinguish seropositive and negative individuals. Furthermore, the continuous variation in antibody levels is ignored thereby potentially reducing the precision of the estimate.

Methods

An age-specific density model which mimics antibody acquisition and loss was developed to make full use of the information provided by serological measures of antibody levels. This was fitted to blood-stage antibody density data from 12 villages at varying transmission intensity in Northern Tanzania to estimate the exposure rate as an alternative measure of transmission intensity.

Results

The results show a high correlation between the exposure rate estimates obtained and the estimated SCR obtained from a catalytic model (r = 0.95) and with two derived measures of EIR (r = 0.74 and r = 0.81). Estimates of exposure rate obtained with the density model were also more precise than those derived from catalytic models.

Conclusion

This approach, if validated across different epidemiological settings, could be a useful alternative framework for quantifying transmission intensity, which makes more complete use of serological data.

Keywords

  • Malaria
  • Model
  • Serology
  • Measuring transmission intensity
  • Antibody titre
  • Cross-sectional data
  • Force of infection
  • Plasmodium falciparum

Background

Malaria remains a leading cause of morbidity and mortality worldwide [1], with heterogeneous levels of endemicity across the globe [2]. Measuring malaria transmission intensity is a key element of monitoring changes in transmission and assessing the impact of anti-malaria interventions. The reference standard historically used for reporting malaria transmission intensity is the entomological inoculation rate (EIR), defined as the number of infectious bites per person per year (ibppy), estimated by catching mosquitoes most commonly using light-traps and then dissecting these mosquitoes to estimate the sporozoite (infectious) rate. Despite its usefulness in providing a direct estimate of the force of infection, this method is time consuming, expensive and can lack precision, especially in low endemicity areas. The prevalence of individuals carrying the parasite—termed parasite prevalence (PrP)—is an alternative measurement of malaria transmission intensity. This can be, estimated in cross-sectional surveys using microscopy, rapid diagnostic test (RDT), or increasingly using PCR methods to detect infection in individuals. This approach is in widespread use, with population-level representative samples now undertaken in many malaria endemic areas as part of national Malaria Indicator Surveys [3]. However, although PrP can be estimated rapidly in populations, it is subject to seasonal variation (e.g. 30 % variation between the peak and trough in highly seasonal areas in West Africa [4]), is affected by anti-malarial treatment levels, requires highly skilled staff (for microscopy and PCR), and lacks precision in low endemicity settings.

Serological data, which measures antibody responses to one or more malaria specific antigens, offer an alternative means to estimate past exposure to malaria [5]. Such data were used historically during the Global Malaria Eradication Programme [6] and after [79] to monitor changes in transmission. There are several advantages to using serology in lower transmission settings where obtaining precise estimates of PrP can be prohibitively expensive. Serological assays such as ELISA are simple, quick and cheap to perform. Antibodies persist for months or years after infection, therefore, the effect of seasonality in transmission is smoothed out. Also, the longevity of antibodies means that seroprevalence remains sufficiently high in low transmission settings to obtain meaningful estimates with achievable sample sizes. Serological data are typically analysed using catalytic models to estimate the antibody seroconversion rate (SCR)—the rate at which seronegative individuals become seropositive—as a proxy for the force of infection [8, 9]. However, one limitation of this method is that it is necessary to distinguish seropositives from seronegatives using continuous measures of antibody levels. For malaria this is typically achieved using sera from European (i.e. unexposed) volunteers to define a cut-off. However, there may be underlying differences between the immune responses in these unexposed volunteers and those living in endemic countries. An alternative approach is to fit mixture models to the bi-modal distribution of antibody levels [1214]. However, this approach can be problematic in highly endemic areas where a large proportion of the population are seropositive.

Here is presented a continuous model of the acquisition and loss of antibodies which can be fitted to individual-level data on measured antibody levels from cross sectional surveys. An advantage of this approach is that it takes into consideration the full information contained in measurement of antibodies levels rather than reducing the data to seropositive/seronegative status. Estimates of malaria transmission intensity are thus derived without the use of cut-offs and with better precision than estimates obtained with binary seroprevalence data. The utility of the approach is illustrated by fitting the model to data from a cross-sectional survey in Tanzania across areas with differing endemicity and compare the results to traditional measures of malaria transmission intensity as well as to estimates obtained from catalytic models.

Methods

Data

Two age-stratified cross-sectional surveys were conducted in Tanzania. Individuals aged 0–46 years of age were sampled from 12 villages across three altitudes transects in North Pare, South Pare and West Usambara. In each transect, villages were selected at high, medium and low altitude, reflecting a range in malaria transmission intensity [10]. As infants may present with maternal antibodies, only individuals between one and 46 years were included in the analysis. Previous publications have already described the study design and epidemiology of malaria in this area [15]. Ethical approval was obtained from the institutional review boards of the National Institute of Medical Research of Tanzania, Kilimanjaro Christian Medical Centre, and the London School of Hygiene and Tropical Medicine. Individuals’ sera were collected and antibodies to the asexual stage merozoite antigens, Merozoite Surface Protein (MSP-119) and Apical Membrane Protein (AMA-1), determined using ELISA [16]. Only anti-MSP-119 antibodies were considered for these analyses, as these are assumed to be more immunogenic than anti-AMA-1 antibodies [1719]. Measurements were recorded as optical densities which were log-transformed prior to analysis with measurements below the limit of detection (LoD) assigned an approximate LoD of 0.01.

Density model

Model specification

A mathematical model was developed to describe the dynamics of acquisition and loss of antibodies in the population. The model assumes that, following exposure to an infectious bite which occurs at rate λ, an individual’s antibody level is boosted by δ(x t ) where x t is the base-10 logarithm of antibody density. In the absence of exposure, antibodies are assumed to decay exponentially at a constant rate ρ. Let y(x, t) denote the proportion of the population with antibody level xat time t and K(x *, x) denote the probability that individuals with (log10) antibody level x at time t are boosted to level x *(x * > x) on exposure to an infectious bite between t and δt with K(x *, x) = 0 if x * ≤ x (∫ x K(x *x) dx * = 1). Then the distribution of antibody levels in the population is given by:
$$\begin{aligned} \frac{\partial y(x,t)}{\partial t}& = \int_{0}^{\infty } {\lambda \left[ {K(x,x^{*} ),y(x^{*} ,t) - K(x^{*} ,x)y(x,t)} \right]} dx^{*}\\ &\quad + \rho \frac{\partial y(x,t)}{\partial x} \end{aligned}$$

With ∫λK(x, x *) y(x *) dx * corresponding to all the individuals with lower antibody levels who get boosted to level xupon exposure. Similarly ∫λK(x *, x) y (x)dx * corresponds to the individuals with antibody levels xwho get boosted to a higher level upon exposure and \(\rho \frac{\partial y(x)}{\partial x}\) corresponds to the individuals losing their antibodies.

The model was numerically approximated by a version in which the log10 antibody density variable, x, was discretized by dividing the range of the variable into N compartments each of width D, with x i denoting the value of (log10) antibody density at the mid-point of antibody class i. The first class represents measurements below the LoD, x min. Here N = 51, with D = 0.052 and x min  = −2. The resulting discrete model describes the dynamics of the proportion of the population in each antibody density category i, denotedy i , and is defined by the following set of ordinary differential equations:

$$\frac{{dy_{i} }}{dt} = \lambda \sum\limits_{\begin{subarray}{l} j < i \\ i \ne 1 \end{subarray} } {k_{ij} y_{j} + \frac{\rho }{\varDelta }y_{j + 1} - \lambda \sum\limits_{\begin{subarray}{l} h > i \\ i \ne N \end{subarray} } {k_{hi} y_{i} - } } \frac{\rho }{\varDelta }y_{i} \begin{array}{*{20}c} {} & {} \\ \end{array} 1 \le i \le N$$
where h, i, j index the N antibody level classes. The rates of exposure and decay of antibodies, λ and ρ, are assumed to be independent of antibody density and age. The probability that following exposure, antibody levels are boosted to class i from class j, k ij , is distributed according to a discretized lognormal distribution:
$$k_{ij} = \left\{ \begin{array}{ll} 0 &\quad{{if }\quad}i < j\\ F(x_{i} + \varDelta /2 - x_{j} ;\delta (x_{j} ),S) - F(x_{i} - \varDelta /2 - x_{j} ;\delta (x_{j} ),S)&\quad{{if }\quad}j \le i < N \\ 1 - F\left( {x_{N - 1} + \varDelta /2 - x_{j} ;\delta (x_{j} ),S} \right) &\quad{{if }\quad}i = N \end{array} \right.$$
where F(z, δ(x), S)is the cumulative density function at point z of the lognormal distribution with mean δ(x) and standard deviation S. δ(x) is the mean boost size, a function of the current log10 antibody level, x, assumed to be given by:
$$\delta (x) = \left\{ \begin{array}{ll} ae^{-bx} &\quad {{if }\quad}x > x_{{\rm min} } \\ \eta &\quad {{otherwise}} \end{array} \right.$$
where a, b and η are parameters. This model assumes that exposure increases the log of antibody density by a decreasing amount as current density increases.

The model is run at equilibrium and constant malaria exposure over the years is assumed. As a result, age of individuals is considered as a proxy for time.

Parameter estimation

A Bayesian approach was used to estimate the model parameters, summarized in Table 1, by fitting the model to the data from the 12 villages simultaneously, allowing only the exposure, λv, to vary by village. The rate of decay of antibodies was fixed to 0.7 years−1. Using θ to denote the estimated parameter vector and D the data, the multinomial log-likelihood is given by:
Table 1

Model parameters and their maximum values

Parameter

Description

Maximum value

Estimates median (95 % CrI)

λ v

Rate of exposure for village v

100

See Table 2

a

Maximum antibody boost size on exposure

10

0.242 (0.237–0.246)

b

Slope of dependence of antibody boost on current log10 optical density

10

0.090 (0.035–0.121)

s

Standard deviation for boost size distribution

10

0.018 (0.006–0.032)

η

Mean boost for individuals with no current antibody

100

0.026 (0.025–0.029)

$$l = \log (P(\,D|\theta )) = \sum\limits_{v} {\sum\limits_{t \ne 0} {\sum\limits_{i} {n_{i,t,v} \log (} } } y_{i,t,v} )$$

Here n i,t,v and y i,t,v are, respectively, the observed number and predicted proportion of individuals in antibody category i in village v at age t. The differential equations were numerically integrated in C using the Runge–Kutta method [20]. MCMC methods (using a Metropolis–Hastings algorithm [21]) were used to calculate the posterior distribution of the parameters. As all parameters were strictly positive we used a log-normal random walk proposal density and assumed uniform priors on [0, max], where max was the maximum permitted value of each parameter as listed in Table 1. Two runs of 500,000 iterations were performed for each run of the MCMC algorithm with a burn-in period of 50,000 steps. Chain convergence was checked visually. The output was then recorded every 200 iterations to generate a sample from the posterior distribution. The standard deviation of the proposal distribution was tuned in order to achieve appropriate mixing of the chains and an acceptance rate close to 20 % [22].

Catalytic model

A comparison of the estimates with those obtained using a previously described catalytic model [23] was performed. In this simple model the proportion of individuals who are seropositive at age t is given by:
$$y_{c} (t) = \frac{{\lambda_{c} }}{{\lambda_{c} + \rho_{c} }}\left[ {1 - e^{{ - (\lambda_{c} + \rho_{c} )t}} } \right]$$
where λ C is the mean annual rate of conversion from seronegative to seropositive and ρ C the mean annual rate of reversion from seropositive to seronegative. It should be noted that (for similar decay rates, r) the SCR estimated with the catalytic model is expected to be smaller in absolute magnitude than the exposure rate calculated from the density model, as exposure in the density model may not always lead to an antibody boost sufficient to caused seroconversion by the criterion used in the catalytic model. A Bayesian MCMC approach as described above was used for parameter estimation. The model was fitted to all villages simultaneously, again allowing λ C to vary by village but with the constraint of estimating a single value forρ C across all villages. Two methods were considered to define individual’s seropositivity. In the first a fixed cut-off value of antibody density of 0.5 was used based on data from non-exposed European sera [10]. In the second a mixture model was fitted to the antibody level distribution for all the village’s data combined across all age groups. The mixture model assumes that the population is composed of a subpopulation of seropositive individuals making up a proportion p of the whole population, and a seronegative subpopulation containing the rest of the population. Antibody levels of individuals in each sub-population are normally distributed with parameters (μ 1, σ 1) for the seronegative group and (μ 2, σ 2) for the seropositive group. The parameters of the mixture model were estimated using MCMC methods with the following likelihood:
$$L(\mu_{1} ,\sigma_{1} ,\mu_{2} ,\sigma_{2} ,p;X) = \mathop \varPi \nolimits_{i = 1}^{n} \{ (1 - p)\phi (X_{i} ,\mu_{1} ,\sigma_{1} ) + p\phi (X_{i} ,\mu_{2} ,\sigma_{2} )\}$$
where \(X_{i}\) is the value of antibody levels X i {X 1, …, X n } and φ(X, μ, σ) represents the probability density function for an observation X from a Normal distribution \(N\) with mean μ and standard deviation σ. So that the overall distribution of antibody levels is (1 − p) N(μ 1, σ 1 2 ) + pN(μ 2, σ 2 2 ). For comparison with studies that use fixed cutoff, analyses were also conducted in which individuals with an antibody level greater than μ 1 + 3σ 1 are defined as seropositive.

Results

Data from 5227 individuals in the 12 villages were included in the analysis. The antibody levels of individuals in each village stratified by age are shown in Fig. 1. The overall trend shows an increase in mean antibody level in adults in the village with decreasing altitude and hence increasing transmission intensity in this setting [2426]. As previously noted [10], two villages—Ngulu and Funta—have higher than expected antibody densities, suggesting higher transmission despite being at medium altitude. As expected, the trend is for antibody density to increase with age in each village, representing cumulative exposure to infection.
Fig. 1
Fig. 1

Antibody levels associated with age and village altitude for MSP-119 antigens. Black dots indicate median estimates of antibody density for the actual data for each village, while median and 95 % credible intervals are represented in red for the model fit (line for median and shaded area for the credible interval). In each transect (each column), villages are presented with decreasing altitude/increasing transmission intensity from top to bottom

The fitted density model is able to capture antibody density patterns across most of the villages (Fig. 1). In each transect, North Pare, South Pare and West Usambara, our estimate of the exposure rate increased with increasing transmission intensity (as indicated by decreasing altitude), with the exception, as expected, of Ngulu and Funta villages (Fig. 2).
Fig. 2
Fig. 2

Estimated exposure rate (median ± 95 % CrI) by village. In each transect (presented by boxes) North Pare, (a), South Pare (b) and West Usambara, (c), the altitude for each village decreases from left to right and hence transmission intensity increases from left to right

The maximum boost size, a, was estimated to be 0.242 [95 % credible interval (CrI) 0.237–0.246] for individuals with antibody level above the detection threshold, and the slope of decline of boost with increasing current antibody level to be b = 0.09 (95 % CrI 0.035–0.121). As illustrated in Fig. 3, the low value of b means that the estimated antibody mean boost size declines approximately linearly with the current log10 antibody level. For the previously unexposed population, the mean boost size was estimated to be η = 0.026 (95 % CrI 0.025–0.029). The lack of overlapping CrIs for the estimates of a and η indicate that data allow us to distinguish the antibody responses of individuals who have never been exposed from those previously exposed who have low levels of antibodies.
Fig. 3
Fig. 3

The estimated dependence of the mean antibody boost size (δ) on the exposed individual’s current log10 antibody level. Its 95 % Credible Interval is represented by the shaded area. η denotes the antibody boost size for individual lacking current circulating antibodies (Median and 95 % CrI)

A high correlation was observed between the estimates of exposure rates obtained using the density model, and that estimated by fitting a catalytic model to the data, using either European controls or a mixture model to derive the cut-off (see Fig. 4a). As anticipated, villages with high altitude (Bwambo, Kilomeni, Mpinji and Kwadoe) had lower estimates of exposure while villages with lower altitude (Kadando and Mgila) had higher estimates of exposure, and this was consistent for estimates obtained with both the density and catalytic models.
Fig. 4
Fig. 4

Comparison with seroconversion rate a Association of median estimates of exposure rates for different villages estimated using the antibody density model (x-axis) and the catalytic model (y-axis) using both European control (squares) and a mixture model to define the cut-off (circles). Note that the two measures are not equivalent but are highly correlated. b Plot of coefficient of variation (standard deviation of posterior/mean of posterior) of the exposure rate estimated from both the density model (black circles) and the catalytic model (x) using European controls. Data are presented for each village categorized by transect

One of the advantages of estimating the exposure rate using the density model rather than SCR is that it makes fuller use of the continuous nature of the data, thus potentially increasing inferential power. The coefficient of variation of the exposure rate estimates, a measure of the precision of those estimates, is consistently smaller for the density model estimates than for the catalytic model estimates (Fig. 4b), demonstrating an increase in the precision of the estimate of transmission intensity obtained by fitting a density model rather than a catalytic model. This result is more marked for villages with lower transmission rates (higher altitude).

Figure 5 shows that the estimates of exposure rates were also highly correlated with the two different estimates of the EIR available for the study villages (and listed in Table 2).
Fig. 5
Fig. 5

Comparison of metrics across villages. a EIR1 (calculated from altitude) with EIR2 (calculated from parasite prevalence); b exposure rate estimated with the antibody density model with EIR1; c exposure rate estimated with the antibody density model with EIR2. Speaman’s correlation coefficientis denoted by r

Table 2

Summary of study data

Transect

Village

Altitude

N

PrP

EIR1

EIR2

Estimated λv median (95 % CrI)

North Pare

Kilomeni (Ki)

1556

411

1

0.047

0.08

3.087 (2.995–3.178)

Lambo (La)

1188

355

10

0.38

1.11

3.287 (3.195–3.379)

Ngulu (Ng)

832

486

8

2.9

0.84

4.109 (4.004–4.217)

Kambi ya Simba (KyS)

746

494

10

4.7

1.11

4.071 (3.972–4.173)

South Pare

Bwanbo (Bw)

1598

485

3

0.037

0.26

3.085 (3.001–3.174)

Mpinji (Mp)

1445

461

2

0.088

0.17

3.071 (2.983–3.158)

Goha (Go)

1163

453

13

0.44

1.57

3.522 (3.437–3.607)

Kadando (Ka)

528

382

25

16

4.45

4.339 (4.212–4.477)

West Usambara

Kwadoe (Kw)

1564

357

4

0.045

0.37

3.138 (3.037–3.236)

Funta (Fu)

1240

429

17

0.28

2.38

4.364 (4.243–4.494)

Tamota (Ta)

1055

449

19

0.81

2.74

3.857 (3.765–3.95)

Mgila (Mg)

375

465

34

39

8.31

4.326 (4.198–4.463)

PrP denotes parasite prevalence measured with microscopy in children 0–4 years of age

EIR1 is estimated from altitude [EIR = 331.5 exp (−0.0057 × altitude) [10]]

EIR2 was estimated from the parasite prevalence data using a previously published relationship [39]

Discussion

The utility of serological data to measure malaria transmission intensity has gained recognition in recent years [10, 11, 27] and is increasingly being incorporated in cross-sectional and longitudinal studies to monitor changes in transmission [2832], identify “hotspots” of transmission [3335] and to identify high risk groups [36]. One of the key advantages of the methods is their ease of use in the field, with new laboratory techniques enabling serological responses to multiple antigens to be made from dried blood spots that can be stored and transported without the need for refrigeration [37]. Classically the approach to analysing such data has been to distinguish seropositives from seronegatives using a cut-off value informed by unexposed European control populations [10, 11]. This has recognized limitations as the European control population may differ genetically in their immunological response to infection from the populations being analysed [38]. To avoid the need to incorporate a cut-off independent of the antibody background level, a density model was developed and fitted to serological data from a malaria endemic setting in Northern Tanzania.

The results demonstrate that estimates of the exposure rate obtained from fitting such models correlate highly both with previous estimates of the SCR obtained from the catalytic model as well as traditional measures of transmission intensity (the EIR) derived from altitude or from PrP data [10, 39]. Overall, the model fits the field data very well though less for older individuals in some areas, perhaps due to the small sample size in this age group. Simulation studies (not presented here) also showed a good fit to the data. The model, therefore, provides an additional method to estimate transmission intensity from serological data that avoids the need to determine a cut-off between seropositivity and seronegativity disregarding the background antibody level. One alternative approach to using European controls to define a cut-off has been to using a mixture model [12, 13]. Since this method can take into consideration the potential for misclassification of seropositive and seronegative individuals and in addition does not require an external dataset for standardization it provides an alternative appealing method for analysing serological data. However, one limitation of the mixture method approach is that it would not be appropriate in high transmission settings where the antibody distribution of seropositive and seronegative individuals becomes very difficult to distinguish. In contrast, the density model presented here performed equally well across all transmission settings.

An additional advantage of the density model over the catalytic model, as demonstrated here, is the improvement in the precision of the estimate of transmission intensity obtained by fitting a density model compared with that obtained from fitting a catalytic model. In particular, this was notably better in the lower transmission settings. This is not surprising, as by using a density model we use a greater degree of information in the data. However, it is also relevant from a practical perspective, as serological measures are likely to be of greatest use in areas of low transmission intensity where other commonly used measures (in particular PrP) lack precision. Thus by utilizing the full data set, this increased precision is likely to improve the ability to distinguish temporal and spatial trends in settings in which malaria has recently fallen to low levels [40, 41] as well as to monitor low-level transmission in countries working towards local elimination of the parasite [42, 43]. The ability to determine significant changes in transmission such as drops in antibody levels would be an important development for the evaluation of interventions [44, 45]. However, a counter argument to this is that the difference in precision may not be attributable solely to more effective use of continuous data but reflects inherent uncertainty about how exposure translates into prevalence, incorporated in the SCR but not in the exposure rate derived with the density model.

The point of comparison for the newly developed density model was chosen to be what is commonly used in practice, i.e. methods deriving thresholds to define seropositivity [10, 27, 4648]. However, mixture models have long been used to also avoid arbitrary misclassification and directly derive seroprevalences and SCRs [12, 4951]. Surprisingly these methods, despite being increasingly studied in the field of malaria [52] have not yet been adopted as common practices.

Additionally, note that the rate of seroconversion from the catalytic model and exposure rate from the density model measure different quantities. Indeed, the exposure rate measures the incidence of blood stage infections while the rate of seroconversion measures the rate at which individuals become seropositive, i.e. their antibody level increase above the cut-off value. One can, therefore, argue on the limited relevance in comparing these estimates. However it is their correlation with exposure measurements, which is of particular interest for comparing transmission intensity from different settings. Applying this approach to other serological datasets where additional measures of exposure are available is the obvious next step.

This model was developed using anti-MSP-1 antibodies as a proof-of-concept, as theses were assumed to be more immunogenic than other antibodies such as anti-AMA-1. However, such analysis could similarly be performed independently to other antibodies and the model could also be further developed to account for the dynamics of simultaneous antibodies.

Whilst the density model clearly has many advantages over the classical catalytic model, there are some limitations are worth noting. Firstly the density model outlined here, whilst biologically grounded, is a simplistic representation of the true process of antibody-acquisition and loss: it does not take into consideration more complex immune responses (such as interaction between antibody responses to different blood-stage antigens and/or between antibody- and cell-mediated immune responses [53]). Nevertheless, it would be interesting to check what effect, if any, incorporating these factors would have on the estimates of transmission as such data become available. Also, while providing a simplistic representation of the complex process, the model does not consider any age related changes in the affinity of the response [54, 55] nor in the age dependency of biting rates [56]. Whilst such aspects are clearly important, they cannot be estimated from cross-sectional data.

The process of loss of individual’s antibodies assumes a constant rate of decay of antibodies fixed to 0.7 years−1 corresponding to a half-life around 360 days. This value, which originated from a longitudinal study [57], made the assumption that antibodies detected for serological studies produced by long lived plasma cells. However, short-lived antibody responses to merozoite antigens are mostly observed [55, 58, 59] and by ignoring these, the model might over-estimate antibody levels. Sensitivity analyses on the rate of decay (not presented here) have shown a high correlation with exposure rate. As a result, the absolute value for the exposure might not be accurate depending on the assumption made for the rate of decay of antibodies but its usefulness in ranking settings according to their transmission intensity still remain.

Additionally, fitting a density model to serological data is computationally intensive and this may limit its wider utility. However, this can to some extent be overcome given the wider availability of statistical packages that can be used to perform such analyses and through sharing of code. To this end the R code used for this analysis is available from the authors on request. For wider use, the practicalities of the assay generating the serological data would need to include a strong component of standardization for the data between laboratories or conducted at different times to be compared confidently.

Conclusion

In summary, the density model presented here provides a new method for analysing serological data that complements existing widely utilized tools for measuring malaria transmission intensity. No gold standard method has yet been developed for estimating exposure from serological data and this newly develop model represents an additional approach to do this. However, future development of this method is needed to test it against a wider set of transmission settings, incorporate methods for assessing spatial and temporal variations in exposure and to assess its utility in capturing changes in transmission following scaling up of malaria interventions.

Abbreviations

SCR: 

seroconversion rate

PrP: 

parasite prevalence

EIR: 

entomological inoculation rate

CrI: 

credible interval

Declarations

Authors’ contributions

This manuscript has not been submitted or accepted for publication elsewhere. EP performed the analyses, implemented the R code and wrote the manuscript. ACG and NMF designed the study, supervised the analyses and helped preparing the manuscript. CJD provided the datasets, helped coordinate the study and revised the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank all members of Joint Malaria Program community study teams and the laboratory staff at Kilimanjaro Christian Medical College and the National Institute for Medical Research who contributed to the data collection. We are also grateful to Simon Cauchemez, Anne Cori and Jeff Eaton for helpful technical discussions. EP is funded by the Medical Research Council UK. CD is funded by the Wellcome Trust (091924). The funders had no role in the design, execution, or publication of the study.

Competing interests

The authors declare that they have no competing interests.

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

Authors’ Affiliations

(1)
Department of Infectious Disease Epidemiology, MRC Centre for Outbreak Analysis and Modelling, Imperial College London, London, UK
(2)
Department of Epidemiology and Public Health, Swiss Tropical and Public Health Institute, Basel, Switzerland
(3)
University of Basel, Basel, Switzerland
(4)
Department of Immunology, London School of Hygiene and Tropical Medicine, London, UK

References

  1. WHO. World Malaria Report 2014. Geneva: World Health Organization; 2014.Google Scholar
  2. Gething PW, Patil AP, Smith DL, Guerra CA, Elyazar IRF, Johnston GL, et al. A new world malaria map: Plasmodium falciparum endemicity in 2010. Malar J. 2011;10:378.PubMed CentralView ArticlePubMedGoogle Scholar
  3. Malaria Indicators Surveys [http://www.malariasurveys.org].
  4. Oduro AR, Conway DJ, Schellenberg D, Satoguina J, Greenwood BM, Bojang KA. Seroepidemiological and parasitological evaluation of the heterogeneity of malaria infection in the Gambia. Malar J. 2013;12:222.PubMed CentralView ArticlePubMedGoogle Scholar
  5. Voller A, Draper CC. Immunodiagnosis and sero-epidemiology of malaria. Br Med Bull. 1982;38:173–7.PubMedGoogle Scholar
  6. Voller A, Bruce-Chwatt LJ. Serological malaria surveys in Nigeria. Bull World Health Organ. 1968;39:883–97.PubMed CentralPubMedGoogle Scholar
  7. Bruce-Chwatt LJ, Draper CC, Konfortion P, Drapper CC. Sero-epidemiological evidence of eradication of malaria from Mauritius. Lancet. 1973;2:547–51.View ArticlePubMedGoogle Scholar
  8. Ambroise-Thomas P, Wernsdorfer WH, Grab B, Cullen J, Bertagna P. Etude sero-seropidemiologique longitudinale sur le paludisme en Tunisie. Bull World Health Organ. 1976;54:355–67.PubMed CentralPubMedGoogle Scholar
  9. Voller A, Storey J, Molineaux L. A longitudinal study of Plasmodium falciparum malaria in the West African savanna using the ELISA technique. World Health. 1980;58:429–38.Google Scholar
  10. Drakeley CJ, Corran PH, Coleman PG, Tongren JE, McDonald SL, Carneiro I, et al. Estimating medium- and long-term trends in malaria transmission by using serological markers of malaria exposure. Proc Natl Acad Sci U S A. 2005;102:5108–13.PubMed CentralView ArticlePubMedGoogle Scholar
  11. Corran PH, Coleman PG, Riley E, Drakeley CJ. Serology: a robust indicator of malaria transmission intensity? Trends Parasitol. 2007;23:575–82.View ArticlePubMedGoogle Scholar
  12. Vyse AJ, Gay NJ, Hesketh LM, Pebody R, Morgan-Capner P, Miller E. Interpreting serological surveys using mixture models: the seroepidemiology of measles, mumps and rubella in England and Wales at the beginning of the 21st century. Epidemiol Infect. 2006;134:1303–12.PubMed CentralView ArticlePubMedGoogle Scholar
  13. Irion A, Beck HP, Smith TA. Assessment of positivity in immuno-assays with variability in background measurements: a new approach applied to the antibody response to Plasmodium falciparum MSP2. J Immunol Methods. 2002;259:111–8.View ArticlePubMedGoogle Scholar
  14. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  15. Drakeley CJ, Carneiro I, Reyburn H, Malima R, Lusingu JPA, Cox J, et al. Altitude-dependent and -independent variations in Plasmodium falciparum prevalence in northeastern Tanzania. J Infect Dis. 2005;191:1589–98.View ArticlePubMedGoogle Scholar
  16. Bruce-Chwatt LJ, Dodge JS, Draper CC, Topley E, Voller A. Sero-epidemiological studies on population groups previously exposed to malaria. Lancet. 1972;1:512–5.View ArticlePubMedGoogle Scholar
  17. Clark EH, Silva CJ, Weiss GE, Li S, Padilla C, Crompton PD, et al. Plasmodium falciparum malaria in the Peruvian Amazon, a region of low transmission, is associated with immunologic memory. Infect Immun. 2012;80:1583–92.PubMed CentralView ArticlePubMedGoogle Scholar
  18. Arnot DE, Cavanagh DR, Remarque EJ, Creasey AM, Sowa MPK, Morgan WD, et al. Comparative testing of six antigen-based malaria vaccine candidates directed toward merozoite-stage Plasmodium falciparum. Clin Vaccine Immunol. 2008;15:1345–55.PubMed CentralView ArticlePubMedGoogle Scholar
  19. Greenhouse B, Ho B, Hubbard A, Njama-Meya D, Narum DL, Lanar DE, et al. Antibodies to Plasmodium falciparum antigens predict a higher risk of malaria but protection from symptoms once parasitemic. J Infect Dis. 2011;204:19–26.PubMed CentralView ArticlePubMedGoogle Scholar
  20. Press WH, Teukolsky S, Vetterling W, Flannery B. Numerical recipes in C: the art of Scientific Computing. 1992.Google Scholar
  21. Gilks WR, Richardson S, Spiegelhalter DJ. Markov Chain Monte Carlo in Practice. Chapman & Hall; 1996:486.Google Scholar
  22. Roberts GO, Rosenthal JS. Optimal scaling for various Metropolis–Hastings algorithms. Stat Sci. 2001;16:351–67.View ArticleGoogle Scholar
  23. Muench H. Catalytic Models in Epidemiology. Cambridge: Harvard University Press; 1959.View ArticleGoogle Scholar
  24. Bødker AR, Akida J, Shayo D, Kisinza W, Msangeni HA, Lindsay SW. Relationship between altitude and intensity of malaria transmission in the Usambara Mountains. Tanzania. Am J Trop Med Hyg. 2003;40:706–17.Google Scholar
  25. Maxwell CA, Chambo W, Mwaimu M, Magogo F, Carneiro IA, Curtis CF. Variation of malaria transmission and morbidity with altitude in Tanzania and with introduction of alphacypermethrin treated nets. Malar J. 2003;2:28.PubMed CentralView ArticlePubMedGoogle Scholar
  26. Kulkarni AMA, Kweka E, Nyale E, Lyatuu E, Chandramohan D, Rau ME, et al. Entomological evaluation of malaria vectors at different altitudes in Hai District. Northeastern Tanzania. J Med Entomol. 2006;43:580–8.PubMedGoogle Scholar
  27. Stewart L, Gosling R, Griffin J, Gesase S, Campo J, Hashim R, et al. Rapid assessment of malaria transmission using age-specific sero-conversion rates. PLoS One. 2009;4:e6083.PubMed CentralView ArticlePubMedGoogle Scholar
  28. Hay SI, Smith DL, Snow RW. Measuring malaria endemicity from intense to interrupted transmission. Lancet Infect Dis. 2008;8:369–78.PubMed CentralView ArticlePubMedGoogle Scholar
  29. Greenwood BM. Control to elimination: implications for malaria research. Trends Parasitol. 2008;24:449–54.View ArticlePubMedGoogle Scholar
  30. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  31. O’Meara WP, Bejon P, Mwangi TW, Okiro EA, Peshu N, Snow RW, et al. Effect of a fall in malaria transmission on morbidity and mortality in Kilifi, Kenya. Lancet. 2008;372:1555–62.PubMed CentralView ArticlePubMedGoogle Scholar
  32. Snow RW, Amratia P, Kabaria CW, Noor AM, Marsh K. The changing limits and incidence of malaria in Africa: 1939–2009. Adv Parasitol. 2012;78:169–262.PubMed CentralView ArticlePubMedGoogle Scholar
  33. 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 Infect Dis. 2010;201:1764–74.View ArticlePubMedGoogle Scholar
  34. Bousema JT, 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.PubMed CentralView ArticlePubMedGoogle Scholar
  35. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  36. The malERA Consultative group on Monitoring. Evaluation and Surveillance: a research agenda for malaria eradication: monitoring, evaluation, and surveillance. PLoS Med. 2011;8:e1000400.View ArticleGoogle Scholar
  37. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  38. Kwiatkowski DP. How malaria has affected the human genome and what human genetics can teach us about malaria. Am J Hum Genet. 2005;77:171–92.PubMed CentralView ArticlePubMedGoogle Scholar
  39. Griffin JT, Hollingsworth TD, Okell LC, Churcher TS, White M, Hinsley W, et al. Reducing Plasmodium falciparum malaria transmission in Africa: a model-based evaluation of intervention strategies. PLoS Med. 2010;7:e1000324.PubMed CentralView ArticlePubMedGoogle Scholar
  40. Noor AM, Clements AC, Gething PW, Moloney G, Borle M, Shewchuk T, et al. Spatial prediction of Plasmodium falciparum prevalence in Somalia. Malar J. 2008;7:159.PubMed CentralView ArticlePubMedGoogle Scholar
  41. Oesterholt MJAM, Bousema JT, Mwerinde OK, Harris C, Lushino P, Masokoto A, et al. Spatial and temporal variation in malaria transmission in a low endemicity area in northern Tanzania. Malar J. 2006;5:98.PubMed CentralView ArticlePubMedGoogle Scholar
  42. Steenkeste N, Rogers WO, Okell L, Jeanne I, Incardona S, Duval L, et al. Sub-microscopic malaria cases and mixed malaria infection in a remote area of high malaria endemicity in Rattanakiri province, Cambodia: implication for malaria elimination. Malar J. 2010;9:108.PubMed CentralView ArticlePubMedGoogle Scholar
  43. Harris I, Sharrock WW, Bain LM, Gray K-A, Bobogare A, Boaz L, et al. A large proportion of asymptomatic Plasmodium infections with low and sub-microscopic parasite densities in the low transmission setting of Temotu Province, Solomon Islands: challenges for malaria diagnostics in an elimination setting. Malar J. 2010;9:254.PubMed CentralView ArticlePubMedGoogle Scholar
  44. Cook J, Kleinschmidt I, Schwabe C, Nseng G, Bousema T, Corran PH, et al. Serological markers suggest heterogeneity of effectiveness of malaria control interventions on Bioko Island, equatorial Guinea. PLoS One. 2011;6:e25137.PubMed CentralView ArticlePubMedGoogle Scholar
  45. Molineaux L, Gramiccia G. The Garki Project: Research on the Epidemiology and Control of Malaria in the Sudan Savanna of West Africa. Chap 10.Google Scholar
  46. 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:1–15.Google Scholar
  47. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  48. 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–9.PubMed CentralView ArticlePubMedGoogle Scholar
  49. Gay NJ. Analysis of serological surveys using mixture models: application to a survey of parvovirus B19. Stat Med. 1996;15:1567–73.View ArticlePubMedGoogle Scholar
  50. Bollaerts K, Aerts M, Shkedy Z, Faes C, Van der Stede Y, Beutels P, Hens N. Estimating the population prevalence and force of infection directly from antibody titres. Stat Modelling. 2012;12:441–62.View ArticleGoogle Scholar
  51. Gay NJ, Vyse AJ, Enquselassie F, Nigatu W, Nokes DJ. Improving sensitivity of oral fluid testing in IgG prevalence studies: application of mixture models to a rubella antibody survey. Epidemiol Infect. 2003;130:285–91.PubMed CentralView ArticlePubMedGoogle Scholar
  52. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  53. Taylor RR, Egan A, McGuinness D, Jepson A, Adair R, Drakeley C, et al. Selective recognition of malaria antigens by human serum antibodies is not genetically determined but demonstrates some features of clonal imprinting. Int Immunol. 1996;8:905–15.View ArticlePubMedGoogle Scholar
  54. Khosravi A, Hommel M, Sayemiri K. Age-dependent antibody response to Plasmodium falciparum merozoite surface protein 2 (MSP-2). Parasite Immunol. 2011;33:145–57.View ArticlePubMedGoogle Scholar
  55. Akpogheneta OJ, Duah NO, Tetteh KA, Dunyo S, Lanar DE, Pinder M, et al. Duration of naturally acquired antibody responses to blood-stage Plasmodium falciparum is age dependent and antigen specific. Infect Immun. 2008;76:1748–55.PubMed CentralView ArticlePubMedGoogle Scholar
  56. Carnevale P, Frezil JL, Bosseno MF, Le Pont F, Lancien J. Study of agressivity of Anopheles gambiae A in relation to age and sex of human subjects. Bull World Health Organ. 1978;56:147–54.PubMed CentralPubMedGoogle Scholar
  57. White M. Models for measuring and predicting malaria vaccine efficacy. Thesis submitted for PhD examination Department of Infectious Disease Epidemiology. 2012.Google Scholar
  58. Kinyanjui SM, Conway DJ, Lanar DE, Marsh K. IgG antibody responses to Plasmodium falciparum merozoite antigens in Kenyan children have a short half-life. Malar J. 2007;6:82.PubMed CentralView ArticlePubMedGoogle Scholar
  59. Cavanagh DR, Elhassan IM, Roper C, Robinson VJ, Giha H, Holder AA, et al. A longitudinal study of type-specific antibody responses to Plasmodium falciparum merozoite surface protein-1 in an area of unstable malaria in Sudan. J Immunol. 1998;161:347–59.PubMedGoogle Scholar

Copyright

© Pothin et al. 2016

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.

Advertisement