Skip to main content

Modelling dynamic change of malaria transmission in holoendemic setting (Dielmo, Senegal) using longitudinal measures of antibody prevalence to Plasmodium falciparum crude schizonts extract

Abstract

Background

Evaluation of local Plasmodium falciparum malaria transmission has been investigated previously using the reversible catalytic model based on prevalence of antibody responses to single antigen to estimate seroconversion rates. High correlations were observed between seroconversion rates and entomological inoculation rates (EIR). However, in this model, the effects of malaria control interventions and clinical episodes on serological measurements were not assessed. This study monitors the use of antibody responses to P. falciparum crude extracts for assessing malaria transmission, compares seroconversion rates estimated from longitudinal data to those derived from cross-sectional surveys and investigates the effects of malaria control interventions on these measures in an area of declining malaria transmission. In addition, the validity of this model was evaluated by comparison with the alternative model.

Methods

Five cross-sectional surveys were carried out at the end of the wet season in Dielmo, a malaria-endemic Senegalese rural area in 2000, 2002, 2008, 2010 and 2012. Antibodies against schizonts crude extract of a local P. falciparum strain adapted to culture (Pf 07/03) were measured by ELISA. Age-specific seroprevalence model was used both for cross-sectional surveys and longitudinal data (combined data of all surveys).

Results

A total of 1504 plasma samples obtained through several years follow-up of 350 subjects was used in this study. Seroconversion rates based on P. falciparum schizonts crude extract were estimated for each cross-sectional survey and were found strongly correlated with EIR. High variability between SCRs from cross-sectional and longitudinal surveys was observed. In longitudinal studies, the alternative catalytic reversible model adjusted better with serological data than the catalytic model. Clinical malaria attacks and malaria control interventions were found to have significant effect on seroconversion.

Discussion

The results of the study suggested that crude extract was a good serological tool that could be used to assess the level of malaria exposure in areas where malaria transmission is declining. However, additional parameters such as clinical malaria and malaria control interventions must be taken into account for determining serological measurements for more accuracy in transmission assessment.

Background

Falciparum malaria remains a public health priority and a major cause of morbidity and mortality in tropical areas [1]. The observed pathophysiology of Plasmodium falciparum malaria infection is strongly dependent upon endemicity, age and level of immunity [2]. In malaria endemic areas, where transmission is perennial and stable, parasitic tolerance has been described and explained by a partial immunity acquired over years [3, 4]. Age is, therefore, a major indicator of the duration of exposure to malaria parasite. However, acquired immunity is not completely protective but can be effective against clinical symptoms and severe form of the disease [5].

Several interventions have been implemented during the two last decades in the fight against malaria. These interventions include treatment by artemisinin-based combination therapy (ACT), use of long-lasting insecticide-impregnated bed nets for exposed populations and use of rapid diagnostic tests for malaria diagnosis. These interventions led to considerable reduction in the number of clinical episodes and deaths due to malaria [1]. Despite the significant progress achieved, the disease remains a major problem. The World Health Organization (WHO) reported 212 million clinical malaria cases and 429,000 deaths in 2015. Children under 5 years old are particularly susceptible to malaria illness, infection and death. In 2015, malaria killed an estimated 303,000 under-fives globally, including 292,000 in the African Region [1].

Monitoring malaria transmission intensity and strengthening malaria control measures are necessary since incidence of severe disease and mortality increase with transmission intensity [6,7,8]. The ongoing changes in malaria endemicity become a key determinant of the progress achieved. Moreover these changes determine the time required to reach the step where elimination could be foreseen [9]. Standard measurements of malaria transmission based on entomological inoculation rate (EIR) and parasite prevalence are expensive, time-consuming and lack of precision because of micro-heterogeneity of malaria transmission [10,11,12]. In addition, both EIR and parasite prevalence are affected by seasonality [13, 14]. Furthermore, assessing malaria transmission intensity and evaluating the impact of interventions are complicated in areas where transmission has been substantially reduced. Therefore, alternative approaches are required to assess malaria transmission and evaluate intervention programmes.

The use of mathematical models and prevalence of anti-malarial antibodies constitute alternative approaches to evaluate malaria endemicity [15, 16]. Mathematical models allow the determination of seroconversion rate (SCR) and seroreversion rate (SRR) which are, for a given time interval, the rates with which a seronegative subject becomes seropositive and a seropositive subject turns back to seronegative, respectively. Predictive seroprevalence (prevalence of antibody responses) can be calculated by using maximum likelihood methods. Seroprevalence reflects cumulative exposure and thus is less affected by seasonality because antibody responses can persist for years after infection. Serological markers have been previously used to assess malaria transmission intensity [15,16,17], to detect recent changes in malaria endemicity [17, 18], to evaluate effectiveness of malaria eradication efforts [19, 20], and SCR has been shown to correlate with EIR [15]. However, these single antigen-based serological approaches have significant limitations related to the sensitivity of the assay and the antigenic polymorphism that may affect detectable antibody responses [21,22,23].

An alternative method could be the use of a combination of antigens rather than a single one since recent studies have shown that responses to crude parasite extract or multiple antigens can overcome [17, 24, 25] these limitations.

Almost all reported mathematical approaches used the catalytic model without considering the effect of several factors (covariates). Taking these covariates into account could improve the effectiveness of serological models in assessing the intensity of malaria transmission [15, 16, 18, 25,26,27,28,29].

Incorporation of covariates such as the number of clinical attack or the use of long-lasting insecticide-impregnated bed nets (LLIN) in the model, allows to take into account the heterogeneity of the population and the effect of interventions affecting malaria transmission.

Previous studies have estimated SCR with reversible catalytic model using cross-sectional surveys and high correlation was shown between serological measurements of transmission and EIR [15, 16, 18]. However, given the fact that longitudinal surveys are often considered in sero-epidemiological studies and the increasing use of cross-sectional serological data to investigate malaria transmission, it would be important to assess consistency of cross-sectional results with SCR estimated from longitudinal survey in order to validate this approach. In this context, Arnold et al. [30] investigated catalytic model in longitudinal cohort in order to validate this approach. The authors compared SCR using antibody responses against the merozoite surface protein1–19 antigen measured in a prospective longitudinal cohort of children aged 0–11 years old to SCR calculated in a cross-sectional survey in a different group ranging from 0 to 90 years old. The comparison between the two approaches was done exclusively in children up to 11 years.

One of the objectives of the present study was to explore if the incorporation of some variations such as allowing age-varying to SCR or adding covariates in the reversible catalytic model could improve the precision or accuracy of exposure estimates using serological measurement of antibody responses to whole parasite extract. The availability of longitudinal data on a cohort whose age ranged from 0 to 90 years allowed the analysis of SCR and SRR values in different age group as studies have shown that antibody responses conversion and reversion varied with age and antigen [31,32,33,34].

A second objective of this study was to investigate the utility of P. falciparum crude schizonts extract which is a mixture of antigens, to monitor malaria transmission intensity and to detect temporal changes in malaria epidemiology in Dielmo, a Senegalese rural holoendemic area, using the age-specific reversible catalytic model described in previous studies [15, 16]. In Dielmo village, various malaria control interventions have been implemented since 2000, leading to significant reduction in malaria transmission and incidence [35]. Thirdly, the effects of potential covariates such as malaria control interventions and clinical malaria attacks on the SCR and SRR were investigated. Lastly, the use of cross-sectional data versus longitudinal data was also analysed and a comparison was performed between reversible catalytic model and alternative catalytic model in longitudinal data for assessing transmission evolution through SCR or SRR.

Methods

Study area and population

In this study, sera collected during five cross-sectional surveys in 2000, 2002, 2008, 2010 and 2012 from individuals living in Dielmo were used. Villagers are followed up longitudinally since 1990 as part of the Dielmo project [36, 37].

Different interventions were implemented based on recommendations from the National Malaria Control Programme. Treatments included the use of chloroquine in 2000, replaced by the association amodiaquine plus sulfadoxine–pyrimethamine from November 2003 to May 2006, followed by the introduction of artemisinin-based combination therapy (ACT) since June 2006. LLINs have been deployed since 2008. As a result, the village has experienced a dramatic decrease of the different malaria indicators monitored [35]. The incidence of clinical malaria dropped from 771 in 2000 to only 17 cases in 2012 [35]. Malaria transmission decreased considerably with EIR dropping from 500 infective bites/person/year in 2000 to 7.6 infective bites/person/year in 2012 [35]. For the longitudinal study, cross-sectional data were combined and only people who have at least two sera observations during the survey periods were considered in the longitudinal analysis. A total of 1147 sera obtained from 359 subjects aged 0–91 years old were used in the longitudinal analysis.

Laboratory methods

Plasmodium falciparum crude schizonts extract obtained from a culture-adapted P. falciparum strain (Pf 07/03) isolated from an inhabitant of Dielmo, was used as antigens for serological investigations. Procedures of P. falciparum culture and crude extract schizonts (Sch 07/03) preparations have been described previously [38]. Antibody responses were assayed by ELISA, the cut-off value for positive antibody response was defined as previously described [38]. Briefly, antibody responses were considered positive if the optical density ratio (OD ratio = mean OD value of sample divided by mean OD value of naïve sera) was above two.

Statistical analysis

Antibody levels comparison

Comparisons of OD ratio of antibody responses were assessed using Kruskal–Wallis rank test and this has been presented in a previous study [39]. Three age groups were studied: under 5, 5–15 years old and above 15 years old. Odds ratio for prevalence was determined using a simple logistic regression. Seroprevalence was estimated by the proportion of seropositive participants. Kruskal–Wallis rank test was used to compare ODs between surveys. Chi Square and Fischer exact tests were used to compare seroprevalences between age groups. All statistics were done using R software [40].

Reversible catalytic prevalence models for cross-sectional data

Seroconversion and seroreversion rates were calculated using the five observed cross-sectional data separately. At first instance, a previously described age specific reversible catalytic conversion model [15, 16, 18] was adjusted to the binary sero-immunological data (antibody responses to crude P. falciparum schizonts extract) to estimate seroconversion (λ) and seroreversion (ρ) rate, respectively. For each cross-sectional survey, seroprevalence has been calculated as described by Corran et al. [15] according to the following formula:

$$P\left( a \right) = \frac{\lambda }{\lambda + \rho }\left( {1 - \exp \left( { - \left( {\lambda + \rho } \right)a} \right)} \right)$$

λ: the annual seroconversion rate, ρ: the annual seroreversion rate, P(a): seroprevalence for a given age a.

Seroconversion and seroreversion rates represent respectively the force of exposure to malaria parasites over time relative to the intensity of infection in malaria endemic areas and the persistence of antibody response [16]. Secondly, an age specific variation in seroconversion rate (λ) and seroreversion rate (ρ) was investigated by adjusting an alternative catalytic model in which λ and ρ were allowed to change at time-point. Comparison of these models was performed using likelihood ratio test.

Estimation of seroprevalence and seroconversion rates with longitudinal data using reversible catalytic and alternative catalytic prevalence models

Age-specific antibody seroprevalence and seroconversion in longitudinal study were estimated by combining measurements obtained from all five cross sectional surveys. In order to select the best model for the serological longitudinal data, two approaches were used. First, a catalytic model and an alternative catalytic model were applied and λ (seroconversion rate) and ρ (seroreversion rate) were allowed to change at time-point without covariates. Second, these two models with covariates were applied to investigate factors that could influence both seroconversion (λ) and seroreversion (ρ) rates. In this study, covariates which were incorporated in the models were clinical malaria episode and use of bed nets. In these models seroconversion (λ) and seroreversion (ρ) rates were calculated respectively by:

$$\uplambda( {\text{t}}) = \uplambda_{0} *{ \exp }( \beta_{1} *{\text{clinical episode }} + \beta_{2} *{\text{use of mosquito bed nets}})$$
$$\uprho ({\text{t}}) = \uprho_{0} * \exp ( \beta^{\prime }_{1} *{\text{clinical}}\;{\text{episode }} + \beta^{\prime }_{2} *{\text{use}}\;{\text{of}}\;{\text{mosquito}}\;{\text{bed}}\;{\text{nets}})$$

where λ0(t) and ρ0(t) were annual conversion rate from seronegative to seropositive (without covariate effect) and annual reversion rate from seropositive to seronegative (without covariate effect), respectively. This model is similar to the one described by Sepulveda et al. [41], though the authors did not describe a model that allowed the SCR to vary by covariates. The regression coefficients (β 1, β 2, β′ 1, β′ 2)can be interpreted in terms of relative risk (RR) similarly to regression coefficient in the proportional hazards regression model of Cox [42]. Likelihood ratio test (LRT) and Wald test were used to assess the effect of covariates and to select the parsimonious model. Significant cut-off was fixed at 0.05 (p < 0.05). The goal here was to assess the effect of clinical malaria attacks and LLIN-based malaria control interventions on SCR. The reversible catalytic and alternative reversible catalytic models are detailed in Additional file 1.

Comparison between observed and predicted seroprevalence

The observed seroprevalence and predicted seroprevalence were compared using Log likelihood ratio test. Predicted seroprevalence curves and 95% confidence intervals were plotted and the corresponding λ was calculated for each cross-sectional period and both cross-sectional surveys combined. Seroprevalence using observed positive responses data were also represented. Wald test and likelihood ratio test (LRT) were used to test effects of covariates. All analyses were conducted using R software version 3.1.2 [40].

Results

Characteristics of the surveyed population

Table 1 presents epidemiological, demographical and clinical characteristics of the surveyed population. The age characteristics of the surveyed population were described in a previous study [39] Both total and mean number of clinical malaria episodes/person decreased over time. Prior to the introduction of LLINs in Dielmo in 2008, the maximum number of clinical malaria attacks varied between 14 and 17 episodes/person/year). Importantly, the maximum number of clinical malaria episodes/person decreased dramatically after the introduction of LLINs in the study site in 2008 [35] falling to 2 and 1 episode in 2010 and 2012, respectively. Table 2 describes the studied cohort; 35.1% of subjects were observed three times and 15.9% were observed five times during the study. Age medians were relatively similar for all screened times groups (Table 2), except the group that was observed five times during the survey period.

Table 1 Characteristics of survey participants
Table 2 Characteristics of the longitudinal cohort

Prevalence of antibody responses to Sch 07/03 amongst age groups

Prevalence of antibodies responses against P. falciparum crude schizonts extract decreased significantly over surveys. Figure 1 shows the seroprevalence for each survey and by age group. Seroprevalence increased significantly with age per survey. Seroprevalence decreased over time during surveys in both children and adults. The reduction in seroprevalence was greater after the introduction of bet nets in the village in late 2008 for children and older children (5–15 years). For adults, reduction in seroprevalence was observed later after 2010. The proportion of seropositive subjects varied from 97% in 2000 to 72% in 2012. However, during the high transmission period in 2000, antibody responses were not different between older children (5–15 years) and adults whereas significant differences in seroprevalence were observed between these age groups in 2002, 2008, 2010 and 2012. Prevalence of antibodies was lowest among younger children (0–5 years). In this age group (0–5 years), the highest percentage was observed in 2002 (4.57%), significantly different to the two other age groups. Of note, in the youngest age group, the already low seroprevalence continued to decrease after the introduction of mosquito nets in 2008 despite a small increase in 2012 (Fig. 1).

Fig. 1
figure 1

Seroprevalence to P. falciparum malaria schizonts sch07/03 with age groups by survey

In older children (5–15 years), seroprevalence decreased from 2000 to 2010 with a high rate of decay between 2008 and 2010 (79.27% of responders in 2008 and 50.98% in 2010). The rebound observed in 2012 in the young children (0–5 years) was also found in the older children group (5–15 years) (Fig. 1). For the adult group, seroprevalence decreased slightly between 2000 and 2012. Despite the decrease observed before and after 2008, the difference was not significant (p = 0.12).

Analysis of seroconversion rate and seroprevalence curves by cross-sectional survey

Seroprevalence was estimated using an age-dependent model, a reversible catalytic conversion model described previously [16]. However, the SCR was not fixed to one value, but could be modified [16, 25]. Figure 2 shows an increase in proportion of seropositivity with participant’s age. Comparison of the different predicted curves at different time points showed a difference in the rate of acquisition of positive responses with a decreased rate from 2000 to 2012. However, visual assessment of the plot suggested a poor fitting of the model for all surveyed time points. For instance, the observed seroprevalence was lower than the predicted one for the younger age group.

Fig. 2
figure 2

Fitting seroprevalence by reversible catalytic model against Pf Sh07/03 for each cross-sectional survey. Black dots represent the observed values, red lines indicate estimated seroconversion with the model and black broken lines indicate the 95% confidence interval (CI) with the likelihood ratio test. The sixth plot present the correlation between the seroconversion rate and the observed EIR given by entomologists

Significant difference of SCRs were observed between 2000 and 2008 (p = 0.024) and between 2002 and 2008 (0.011). SCRs were highest before the introduction of LLINs in 2008 with SCR values of 0.22 [95% CI (0.173–0.280)] and 0.126 [95% CI (0.099–0.159)] in 2000 and 2002 respectively compared to 0.094 [95% CI (0.079–0.112)] in 2008. In addition, a decrease of SCRs was observed after 2008 (Table 3). When an alternative catalytic model that included age specific variation for seroconversion rate was used, the best adjustment to the real data (LRT, p < 0.001) was observed as compared to the catalytic model with constant SCR and SRR.

Table 3 Seroconversion rates from catalytic models in cross-sectional surveys, their corresponding EIR [15] calculated from the serological measures and observed EIR given by entomologists

Model applicability: correlation between SCR and EIR

A positive correlation was observed between EIR and SCR (λ) of antibodies against Pfsch07/03 (r-square = 0.905, p = 0.024). This trend was also shown in Fig. 3a, which presents the evolution of observed EIRs and corresponding EIRs calculated from serological measures using the method described in previous studies, respectively [15, 16].

Fig. 3
figure 3

Trends of observed and estimated entomological inoculate rate (EIR), seroconversion rate (SCR) and seroreversion rate (SRR): a dynamic of observed entomological inoculate rate, EIR (black) calculated from parasitological measurements, seroconversion rate of anti-Pfsch07/03 and their corresponding entomological inoculate obtained by the catalytic model described previously [15], b comparison between observed and corresponding EIR obtained by the log–log calibration [15, 16], c evolution of SCR and SRR by years

The concomitant decrease of both EIR and SCR (λ) during survey period suggested that serological data could be used to assess the level of malaria transmission intensity. Figure 3b shows EIRs estimated from SCRs of antibody responses to Pfsch07/03—compared to observed EIRs measured with entomological methods in the study site (Dielmo) over the 5 years of surveys [35]. Estimated EIR showed a greater variability across survey periods. Despite slight difference between observed EIR and estimated EIRs from all period except for 2002, the corresponding EIRs closely followed up trend of observed EIRs over the survey periods (Fig. 3a).

Figure 3c describes dynamics of both SCR and SRR over the period. Dramatic reduction was observed in estimates seroconversion rates (SCR). In seroreversion rates (SRR), a small reduction was observed over time until 2010 where a slight increase was found.

Analysis of seroconversion rate and seroprevalence curves in longitudinal survey

Seroprevalence results from catalytic model applied to longitudinal data are shown in Fig. 4a. A poor fitting of age-seroprevalence of anti-Pfsch0703 was observed. Seroprevalence in adults group was under-estimated compared to younger age groups, which seemed to be over-estimated. Seroconversion and seroreversion rates estimated from the catalytic model was 0.284 and 0.066, respectively (Table 4).

Fig. 4
figure 4

Fitting seroprevalence by reversible catalytic model (a) and alternative catalytic model (b) for anti-Pf Sh07/03. For all these figures, the red curves were the seroprevalence obtained from models without covariates and their 95% likelihood ratio test confidence intervals and the green ones these from models according to the clinical attacks and the LLINs usage. The triangle points were the seroprevalences calculated among the observed antibodies response

Table 4 Estimation of seroconversion and seroreversion rates from longitudinal and cross-sectional studies

Slight variation of SCRs and SRRs were observed when covariates such as clinical malaria attacks and LLINs usage were taken into account in the catalytic model in which SCR and SRR were estimated respectively at 0.284 (0.230–0.351) and 0.041 (0.030–0.055). The LRT showed a significant difference between these models in favor of the model with covariates (p < 0.001).

Figure 4b presents results of the alternative catalytic model in which age-varying SCR and SRR are allowed. The LRT showed that the alternative model was significantly better than the catalytic model in which SCR and SRR were constant (X 2 = 25.67, p < 0.001). Current SCR (2012) using all sample surveys was calculated and found to be equal to 0.056 [95% CI (0.04–0.11)] from the alternative catalytic model and the corresponding SRR was 0.01 [95% CI (0.0007–0.08)].

Comparison of SCR and SRR from longitudinal study versus cross-sectional surveys

Table 4 shows comparison of SCR and SRR estimated from longitudinal surveys to those from the five cross-sectional surveys by both catalytic models with and without covariates. SCRs from longitudinal approach [SCR (without covariates) = 0.284 (0.230–0.351) and SCR (with covariates) = 0.215 (0.171, 0.269)] were almost equal to the SCR (0.220 (0.173, 0.280) from the cross-sectional survey in 2000. For other cross-sectional surveys (2002, 2008, 2010 and 2012), there were a high variability between SCRs from cross-sectional and longitudinal surveys. In contrast, for the alternative catalytic model with covariates current SCR estimated to 0.046 (0.003–0.052) was comparable to SCR value of 0.042 (0.036–0.049) from the cross-sectional survey in 2012 (Fig. 4b).

Effects of clinical malaria attack and malaria control interventions on seroconversion

To investigate respectively the impact of LLINs usage and clinical malaria episodes on SCR and SRR, catalytic age reversion model with covariates has been applied on data combined from all surveys. Table 5 summarizes the estimation of these parameters. The mosquito nets intervention reduced the seroconversion rate, nevertheless it increased the seroreversion rate. A significant association was found between clinical episodes and SCR with a p value < 0.001.

Table 5 Estimation of parameters [respectively the relative risk (RR)] influencing seroconversion and seroreversion rates from longitudinal studies

Discussion

In malaria pre-elimination context with reduced clinical episodes and transmission intensity due to implemented strategies against the disease, standard assessment tools lack of sensitivity and more sensitive tools are necessary to accurately evaluate both the level of exposure and the impact of malaria control interventions. These tools can also help to detect resurgence in transmission intensity. In recent years, several studies described the estimation of malaria transmission using model based on the seroprevalence of antibody responses against recombinant blood stage antigens, such as MSP1, MSP2, AMA-1, MSP-1_19, GLURP [15,16,17,18, 27,28,29,30, 43, 44] and sporozoite stage antigens of Plasmodium [45]. These serological markers could be potentially sensitive, especially in areas of low transmission [29] as antibody responses, in particular to blood stage antigens, have been shown to persist for several years after transmission has ceased [46,47,48,49]. Thus, antibody responses against blood stage antigens could be detected even if transmission is lower. However, the use of specific single antigen might induce a lack of sensitivity and lead to a significant underestimation of the immune response to P. falciparum infection. The variation in individual immune reactivity as well as the different duration of antibody responses against different antigens [34, 50] and antigenic polymorphism may affect serological responses, thus the results obtained with specific recombinant antigens [21,22,23, 51]. As an example, Ondigo et al. showed that antibody half-life may be short, intermediate or long depending on age and antigen and concluded for a need to associate antigen with similar half-life to improve precision and accuracy of malaria exposure estimates [52]. Recently, Helb et al., in a study aimed at improving tools to reliably measure P. falciparum exposure, observed that P. falciparum-specific antibody responses differ by antigen and, in contrast with Ondigo et al., concluded to the need of selecting and combining different antigens with different kinetics for improving estimates of P. falciparum exposure [53].

These combinations have been demonstrated to improve the accuracy and the precision of malaria infection exposure. However the technologies used (multiplex, microarrays or cytometry) are quite expensive and the required equipment is not always available (or only in very few, high-level laboratories) in the countries where monitoring of malaria infection is needed. In a country with limited resources, there is probably a space for the use of less sophisticated tools, such as ELISA assays, and based on the results of this study, the use of crude extract could be a good choice. In fact, crude extract containing a plurality of antigens can overcome the constraints associated to the choice of specific antigens and give a better sensitivity, allowing detection of P. falciparum exposure in areas of low-level transmission [22, 24].

In this current study, a crude schizont extract of PfSch07/03 was used in ELISA to examine the dynamic of antibody responses in individuals living in Dielmo (Senegal), an area where malaria transmission has substantially decreased [35]. An estimation transmission model based on prevalence of antibody responses against PfSch07/03 was applied to assess the applicability of this model in the study area.

Results showed a significant increase both in level and prevalence of antibody responses with age within surveys. These observations confirmed the general trend that anti-malaria immune responses varied with age in endemic areas [38]. For all age groups, comparisons of both OD magnitude and seropositivity prevalence showed that antibody level and seroprevalence were generally significantly higher in 2000 and decreased over time to 2012. The decreased of antibody responses correlated with the decrease of transmission level induced by control interventions deployed in Dielmo [35]. These results confirmed the high sensitivity of serological responses to detect changes in exposure to malaria infection [15, 16].

Maximum likelihood method from age-specific seroprevalence reversible catalytic model was used to estimate seroconversion rate (SCR). SCR estimates were shown to correlate closely with EIR, the gold standard measurement of malaria transmission. The result was consistent with previous studies that have shown good correlation between EIR and SCR in low transmission areas [15]. Importantly, heterogeneity of SCR by time-point survey corresponded to that of transmission in the study site. Indeed, the study locality has for long been a holoendemic malaria area in which transmission has varied from 142.5 (in 1990) to 482 (in 2000) infective bites/person/year [35]. A recent study in this locality has shown that both transmission level and malaria attacks dropped since the introduction of LLINs in 2008 combined with ACT [35]. The data showed the recent change about the level of transmission in Dielmo. Some cautiousness should be considered on the interpretation of the relationship between SCR and EIR. The variability of these measures should be taken into account in order to avoid misinterpretation. However, estimating the level of malaria transmission with the serological tool is relatively simple and could have an interest in the assessment of control interventions all the more so the antibodies response remain detectable for years.

SCRs obtained by combining all sectional data were statistically different from those obtained with time-point survey. In cross-sectional studies, SCRs gave precise information in the level of transmission and detected previous changes in malaria exposure. This difference could be explained by the heterogeneity in the transmission in Dielmo during these periods [35] and by antibody responses variability [38]. In contrast, the study of Arnold et al. [30] found a similarity between the prospective estimate of the seroconversion rate from the longitudinal data and that estimated from the cross-sectional survey using data that covered the same period (cross-sectional sample in 2012). These contrasting results could be explained partly by the difference of samples in these two studies. Indeed the samples used by Arnold et al. [30] included only children aged 0–11 years old. Another explanation could be the dramatic reduction in malaria transmission over the 12 years study period, while in the study by Arnold et al. [30]. the transmission was constant over the study period.

The observation that SCR averaged over the entire period in the longitudinal analysis differed from the weighted average of the stratum-specific SCRs could be the consequence of the large change in SCR values over the study period. In view of the results in the present study. It seemed that the use for longitudinal data of catalytic models with constant SCR and SRR is not suitable in areas with rapid decline in malaria transmission. In these areas models allowed for time-varying SCR and SRR could be more adapted for longitudinal data in order to take into account the changes in SCRs following the reduction of malaria infection exposure.

Investigation of factors which might affect the SCR and SRR showed that both clinical episodes and the use of mosquito nets had significant impact on these parameters. In agreement with previous studies [18] showing that the incidence of the disease has a significant influence on the estimation of SCR, clinical episodes were positively associated with the SCR of antibody responses to PfSch07/03. Similarly, the use of mosquito nets reduced the risk of antibodies jumping from seronegative to seropositive state. This is consistent with previous studies that have shown significant decrease in the level of immunity in the locality after the introduction of LLINs [38].

Results of LRT and Wald test were in favor of the alternative model suggesting that the previously described catalytic model [15, 16] was restrictive both in cross-sectional and in longitudinal studies in Dielmo. However, this model gave some interesting results, among which, the correlation between corresponding EIRs obtained by serological data and observed EIR (entomological measurement). These results suggested the applicability of the catalytic model to serological data. However, as described by Yman et al. [54] the use of catalytic model has some weakness such as the conversion of serological data to binary outcomes through a threshold model and the choice of this threshold. Furthermore, dichotomizing serological measures (seronegative vs seropositive) may lead to a loss of information due to the high range of serological values of the seropositive group. Several authors proposed alternative models using continuous serological data, as the density model [41, 42] or the antibody acquisition model [54]. Interestingly Pothin et al. have shown that while seroconversion rate from the catalytic model and exposure rate from the density model measure different quantities, a high correlation has been observed between these two measures [43]. This consolidates the use of the catalytic model in this study. However, in the antibody acquisition model, Yman et al. found that in settings of moderate and intense transmission, this model based on continuous data gave better precision and accuracy compared to catalytic models [54], even if the authors agreed that this approach has some disadvantages in particular if the data are not log-normally distributed [54]. Arnold et al. proposed the use of quantitative antibody levels as a valuable method to measure changes in transmission or difference in exposure for pathogens that elicit a transient antibody response or for monitoring populations with very high- or very low transmission [55]. Indeed, in previous studies on the population of the present study, in parallel to seroprevalence estimations, quantitative antibody levels were measured in different age groups and compared between two very different periods of malaria transmission [38]. The results showed that the levels of antibody responses added complementary and valuable information to seroprevalence values.

Comparisons of the different methods with the same data could be interesting to envisage.

Conclusions

The present study underlines the capacity of serological responses using crude extracts of P. falciparum to detect changes in malaria transmission and suggests to take into account malaria control interventions in the estimation of the serological measurements of malaria transmission. Further investigations comparing SCRs obtained from recombinant antigens and especially vaccine candidates antigens to those obtained from crude extracts would provide the basis for the choice and validation of a reference serological marker to assess the level of malaria transmission in relation to changing epidemiology. Additionally, estimations of seroconversion rates should take into account the heterogeneity of the transmission within a study period and the methodology used for these estimations should be adapted to the epidemiology context.

References

  1. WHO. World malaria report 2016. Geneva: World Health Organization; 2016.

    Google Scholar 

  2. Snow RW, Omumbo JA, Lowe B, Molyneux CS, Obiero JO, Palmer A, et al. Relation between severe malaria morbidity in children and level of Plasmodium falciparum transmission in Africa. Lancet. 1997;349:1650–4.

    Article  CAS  PubMed  Google Scholar 

  3. Rogier C, Commenges D, Trape JF. Evidence for an age-dependent pyrogenic threshold of Plasmodium falciparum parasitemia in highly endemic populations. Am J Trop Med Hyg. 1996;54:613–9.

    Article  CAS  PubMed  Google Scholar 

  4. Ladeia-Andrade S, Ferreira MU, de Carvalho ME, Curado I, Coura JR. Age-dependent acquisition of protective immunity to malaria in riverine populations of the Amazon Basin of Brazil. Am J Trop Med Hyg. 2009;80:452–9.

    CAS  PubMed  Google Scholar 

  5. Trape JF, Pison G, Spiegel A, Enel C, Rogier C. Combating malaria in Africa. Trends Parasitol. 2002;18:224–30.

    Article  PubMed  Google Scholar 

  6. Chandramohan D, Greenwood B, Cox J, Mswia R, Setel P. Relationship between malaria endemicity and acute febrile illness mortality in children. Bull World Health Organ. 2001;79:375–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. Reyburn H, Mbatia R, Drakeley C, Bruce J, Carneiro I, Olomi R, et al. Association of transmission intensity and age with clinical manifestations and case fatality of several Plasmodium falciparum malaria. J Am Med Assoc. 2005;293:1461–70.

    Article  CAS  Google Scholar 

  8. Snow RW, Marsh K. The consequences of reducing transmission of Plasmodium falciparum in Africa. Adv Parasitol. 2002;52:235–64.

    Article  PubMed  Google Scholar 

  9. Smith DL, Hay SI, Noor AM, Snow RW. Predicting changing malaria risk after expanded insecticide-treated net coverage in Africa. Trends Parasitol. 2009;25:511–6.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Gu W, Mbogo CM, Githure JI, Regens JL, Killeen GF, Swalm CM, et al. Low recovery rates stabilize malaria endemicity in areas of low transmission in coastal Kenya. Acta Trop. 2003;86:71–81.

    Article  PubMed  Google Scholar 

  11. Shililu J, Ghebremeskel T, Mengistu S, Fekadu H, Zerom M, Mbogo C, et al. High seasonal variation in entomologic inoculation rates in Eritrea, a semi-arid region of unstable malaria in Africa. Am J Trop Med Hyg. 2003;69:607–13.

    PubMed  Google Scholar 

  12. Mbogo CM, Mwangangi JM, Nzovu J, Gu W, Yan G, Gunter JT, et al. Spatial and temporal heterogeneity of Anopheles mosquitoes and Plasmodium falciparum transmission along the Kenyan coast. Am J Trop Med Hyg. 2003;68:734–42.

    PubMed  Google Scholar 

  13. Kelly-Hope LA, McKenzie FE. The multiplicity of malaria transmission: a review of entomological inoculation rate measurements and methods across sub-Saharan Africa. Malar J. 2009;8:19.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Robert V, Le Goff G, Andrianaivolambo L, Randimby FM, Domarle O, Randrianarivelojosia M, et al. Moderate transmission but high prevalence of malaria in Madagascar. Int J Parasitol. 2006;36:1273–81.

    Article  PubMed  Google Scholar 

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

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

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Cook J, Reid H, Lavro 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 

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

    Article  PubMed  PubMed Central  Google Scholar 

  19. Bruce-Chwatt LJ. Eradication of malaria in Africa. Ann Parasitol Hum Comp. 1973;48:221–9.

    CAS  PubMed  Google Scholar 

  20. Ambroise-Thomas P, Wernsdorfer WH, Grab B, Cullen J, Bertagna P. Longitudinal sero-epidemiologic study of malaria in Tunisia. Bull World Health Organ. 1976;54:355–67.

    CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  22. Takala SL, Plowe CV. Genetic diversity and malaria vaccine design, testing and efficacy: preventing and overcoming ‘vaccine resistant malaria’. Parasite Immunol. 2009;31:560–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Takala SL, Coulibaly D, Thera MA, Batchelor AH, Cummings MP, Escalante AA, et al. Extreme polymorphism in a vaccine antigen and risk of clinical malaria: implications for vaccine development. Sci Transl Med. 2009;1:2ra5.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Ambrosino E, Chloé D, Orlandi-Pradines E, Remoue F, Toure-Baldé A, Tall A, et al. A multiplex assay for the simultaneous detection of antibodies against 15 Plasmodium falciparum and Anopheles gambiae saliva antigens. Malar J. 2010;9:317.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Wong J, Hamel MJ, Drakeley CJ, Kariuki S, Shi YP, Lal AA, et al. Serological markers for monitoring historical changes in malaria transmission intensity in a highly endemic region of Western Kenya, 1994–2009. Malar J. 2014;13:451.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

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

  28. Rosas-Aguirre A, Llanos-Cuentas A, Speybroeck N, Cook J, Contreras-Mancilla J, Soto V, et al. Assessing malaria transmission in a low endemicity area of north-western Peru. Malar J. 2013;12:339.

    Article  PubMed  PubMed Central  Google Scholar 

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

    Article  PubMed  PubMed Central  Google Scholar 

  30. Arnold BF, Priest JW, Hamlin KL, Moss DM, Colford JM Jr, Lammie PJ. Serological measures of malaria transmission in Haiti: comparison of longitudinal and cross-sectional methods. PLoS ONE. 2014;9:e93684.

    Article  PubMed  PubMed Central  Google Scholar 

  31. McCallum FJ, Persson KEM, Fowkes FJI, Reiling L, Mugyenyi CK, Richards JS, et al. Differing rates of antibody acquisition to merozoite antigens in malaria: implications for immunity and surveillance. J Leukoc Biol. 2017;14:913–25.

    Article  Google Scholar 

  32. Perraut R, Varela ML, Loucoubar C, Niass O, Sidibe A, Tall A, et al. Serological signatures of declining exposure following intensification of integrated malaria control in two rural Senegalese communities. PLoS ONE. 2017;12:e01791.

    Article  Google Scholar 

  33. White MT, Griffin JT, Akpogheneta O, Conway DJ, Koram KA, Riley EM, et al. Dynamics of the antibody response to Plasmodium falciparum infection in African children. J Infect Dis. 2014;11:15–22.

    Google Scholar 

  34. Akpogheneta O, Duah NO, Tetteh KKA, Dunyo S, Lanar D, 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Trape JF, Tall A, Sokhna C, Ly AB, Diagne N, Ndiath O, et al. The rise and fall of malaria in a West African rural community, Dielmo, Senegal, from 1990 to 2012: a 22 year longitudinal study. Lancet Infect Dis. 2014;14:476–88.

    Article  PubMed  Google Scholar 

  36. Trape JF, Rogier C, Konate L, Diagne N, Bouganali H, Canque B, et al. The Dielmo project: a longitudinal study of natural malaria infection and the mechanisms of protective immunity in a community living in a holoendemic area of Senegal. Am J Trop Med Hyg. 1994;51:123–37.

    Article  CAS  PubMed  Google Scholar 

  37. Rogier C, Ly AB, Tall A, Cisse B, Trape JF. Plasmodium falciparum clinical malaria in Dielmo, a holoendemic area in Senegal: no influence of acquired immunity on initial symptomatology and severity of malaria attacks. Am J Trop Med Hyg. 1999;60:410–20.

    Article  CAS  PubMed  Google Scholar 

  38. Diop F, Richard V, Diouf B, Sokhna C, Diagne N, Trape JF, et al. Dramatic declines in seropositivity as determined with crude extracts of Plasmodium falciparum schizonts between 2000 and 2010 in Dielmo and Ndiop, Senegal. Malar J. 2014;13:83.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Niang M, Niass O, Diop F, Diouf B, Faye MM, Diene-Sarr F, et al. Temporal analysis of IgG antibody responses to Plasmodium falciparum antigens in relation to changing malaria epidemiology in a West African setting. Malar J. 2017;16:283.

    Article  PubMed  PubMed Central  Google Scholar 

  40. http://www.R-project.org. Accessed 29 Sept 2017.

  41. Sepulveda N, Stresman G, White MT, Drakeley CJ. Current mathematical models for analyzing anti-malaria antibody data with an eye to malaria elimination and eradication. J Immun Res. 2015;2015:e738030.

    Google Scholar 

  42. Cox DR. Regression models and life tables (with discussion). J R Stat Soc B. 1972;34:187–220.

    Google Scholar 

  43. Pothin E, Ferguson MN, Drakeley CJ, Ghani AC. Estimating malaria transmission intensity from Plasmodium falciparum serological data using antibody density models. Malar J. 2016;15:79.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Dewasurendra RL, Dias JN, Sepulveda N, Gunawardena GSA, Chandrasekharan N, Drakeley C, et al. Effectiveness of a serological tool to predict malaria transmission intensity in an elimination setting. BMC Infect Dis. 2017;17:49.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Kusi KA, Bosomprah S, Dodoo D, Kyei-Baafour E, Dickson EK, Mensah D, et al. Anti-sporozoite antibodies as alternative markers for malaria transmission intensity estimation. Malar J. 2014;13:103.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Druilhe P, Pradier O, Marc JP, Miltgen F, Mazier D, Parent G. Levels of antibodies to Plasmodium falciparum sporozoite surface antigens reflect malaria transmission rates and are persistent in the absence of reinfection. Infect Immun. 1986;53:393–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Wipasa J, Suphavilai C, Okell LC, Cook J, Corran PH, Thaikla K, et al. Long-lived antibody and B Cell memory responses to the human malaria parasites, Plasmodium falciparum and Plasmodium vivax. PLoS Pathog. 2010;6:e1000770.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Morais CG, Soares IS, Carvalho LH, Fontes CJ, Krettli AU, Braga EM. Antibodies to Plasmodium vivax apical membrane antigen 1: persistence and correlation with malaria transmission intensity. Am J Trop Med Hyg. 2006;75:582–7.

    CAS  PubMed  Google Scholar 

  49. Braga EM, Fontes CJ, Krettli AU. Persistence of humoral response against sporozoite and blood-stage malaria antigens 7 years after a brief exposure to Plasmodium vivax. J Infect Dis. 1998;177:1132–5.

    Article  CAS  PubMed  Google Scholar 

  50. Kinyanjui SM, Bejon P, Osier FH, Bull PC, Marsh K. What you see is not what you get: implications of the brevity of antibody responses to malaria antigens and transmission heterogeneity in longitudinal studies of malaria immunity. Malar J. 2009;8:242.

    Article  PubMed  PubMed Central  Google Scholar 

  51. Gray JC, Corran PH, Mangia E, Gaunt MW, Li Q, Tetteh KKA, et al. Profiling the antibody immune response against blood stage malaria vaccine candidates. Clin Chem. 2007;53:1244–53.

    Article  CAS  PubMed  Google Scholar 

  52. Ondigo BN, Hodges JS, Ireland KF, Magak NG, Lanar DE, Dutta SJ, et al. Estimation of recent and long-term malaria transmission in a population by antibody testing to multiple Plasmodium falciparum antigens. J Infect Dis. 2014;210:1123–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Helb DA, Tetteh KKA, Felgner PL, Skinner J, Hubbard A, Arinaitwe E, et al. Novel serological 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  PubMed Central  Google Scholar 

  54. Yman V, White MT, Rono J, Arcà B, Osier FH, Troye-Blomberg M, et al. Antibody acquisition models: a new tool for serological surveillance of malaria transmission intensity. Sci Rep. 2016;6:19472.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Arnold BF, Van der Laan MJ, Hubbard AE, Steel C, Kubofcik J, Hamlin KL, et al. Measuring changes in transmission of neglected tropical diseases, malaria, and enteric pathogens from quantitative antibody levels. PLoS Negl Trop Dis. 2017;11:e0005616.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

AKD and ATB designed the study. FD, BD, MMF performed the laboratory experiments. JF, ND, CS, JFT, RP, AT and ATB coordinated the fieldwork, sampling, archiving and withdrawal. ON, PSP, AKD and ATB performed the modelisation studies. ON, PSP, AKD, MN and ATB wrote the manuscript. All authors contributed to the revision of the manuscript, and have seen and approved the final version. All authors read and approved the final manuscript

Acknowledgements

The authors express their gratitude to the population of Dielmo and Ndiop for their participation in the study. The healthcare workers in Dielmo and Ndiop are dully acknowledged for their cooperation in conducting this study. Moreover, the authors thank the CEA-MITIC and the SCAC embassy of French for their grant.

Competing interests

The authors declare that they have no competing interests.

Funding

The funder has no role in the design of the study and data collection, analysis, and interpretation of data and in writing the manuscript.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Aïssatou Toure Balde.

Additional file

12936_2017_2052_MOESM1_ESM.docx

Additional file 1. This supporting information contains statistical details. It includes information about models used in the analysis and the estimates methodologies.

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

Niass, O., Saint-Pierre, P., Niang, M. et al. Modelling dynamic change of malaria transmission in holoendemic setting (Dielmo, Senegal) using longitudinal measures of antibody prevalence to Plasmodium falciparum crude schizonts extract. Malar J 16, 409 (2017). https://doi.org/10.1186/s12936-017-2052-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-017-2052-0

Keywords