Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2018 Jul 25:7:e35832.
doi: 10.7554/eLife.35832.

Quantification of anti-parasite and anti-disease immunity to malaria as a function of age and exposure

Affiliations

Quantification of anti-parasite and anti-disease immunity to malaria as a function of age and exposure

Isabel Rodriguez-Barraquer et al. Elife. .

Abstract

Fundamental gaps remain in our understanding of how immunity to malaria develops. We used detailed clinical and entomological data from parallel cohort studies conducted across the malaria transmission spectrum in Uganda to quantify the development of immunity against symptomatic P. falciparum as a function of age and transmission intensity. We focus on: anti-parasite immunity (i.e. ability to control parasite densities) and anti-disease immunity (i.e. ability to tolerate higher parasite densities without fever). Our findings suggest a strong effect of age on both types of immunity, not explained by cumulative-exposure. They also show an independent effect of exposure, where children living in moderate/high transmission settings develop immunity faster as transmission increases. Surprisingly, children in the lowest transmission setting appear to develop immunity more efficiently than those living in moderate transmission settings. Anti-parasite and anti-disease immunity develop in parallel, reducing the probability of experiencing symptomatic malaria upon each subsequent P. falciparum infection.

Keywords: P. falciparum; anti-disease immunity; anti-parasite immunity; epidemiology; global health; immunity; malaria.

PubMed Disclaimer

Conflict of interest statement

IR, EA, PJ, MK, PR, JR, GD, JN, SS, MK, CD, IS, DS, BG No competing interests declared

Figures

Figure 1.
Figure 1.. Incidence and prevalence of malaria as a function of age.
(a) and prevalence of asymptomatic parasitemia (b) in the three study sites as a function of age, modeled using generalized additive models (GAMS). Shaded areas represent 95% confidence bounds.
Figure 2.
Figure 2.. aEIR as a metric of individual exposure.
(a) Distribution of the average annual entomological inoculation rate (aEIR) experienced by the study households in the three study sites. (b) Correlation between the measured aEIRs and the estimated individual hazards of infection.
Figure 3.
Figure 3.. Changes in parasite densities and objective temprature as a function of age and exposure.
Trends in parasite densities. (a) recorded during symptomatic (passive surveillance) infections and routine (active surveillance) visits as a function of age (left) and aEIR (right); and trends in the objective temperature (b) recorded during visits in which participants were found to have a parasite density between 50,000 parasites/μL and 200,000 parasites/μL, as a function of age (left) and aEIR (right). Each point represents a measurement obtained during a study visit. The median and interquartile range are shown in black.
Figure 4.
Figure 4.. Results of best fitting models quantifying immunity.
(a), anti-disease immunity (b) and overall immunity against symptomatic malaria (c). Each plot shows, for specific ages and aEIRs, the expected parasite density (/μL) (a), objective temperature given a density of 40,000 parasites/μL (b) and the probability of developing symptomatic malaria upon infection (c), estimated using the best fitting model. 95% confidence intervals of the estimates are also shown.
Figure 5.
Figure 5.. Anti-parasite and Anti-disease immunity as a function of age and exposure.
(a) and anti-disease immunity (b). These results are similar to those presented in Figure 4, but for the full range of ages and aEIRs included in the data. Panel (a) shows expected parasite densities (parasites/μL, log 10) upon infection for different ages and levels of exposure (aEIR). Panel (b) shows the ‘fever threshold’ or ‘pyrogenic density’, the minimum parasite densities (parasites/μL, log 10) associated with fever (temperature 38°C or greater), again as a function of age and exposure.
Figure 5—figure supplement 1.
Figure 5—figure supplement 1.. Confidence bounds of best fitting model quantifying anti-parasite immunity.
Left panel (a) is equivalent to Figure 5a, and shows the expected parasite densities (log 10, parasites/μL) after infection at different ages and levels of exposure (aEIR). Panel (b) shows the lower and upper 95% confidence bounds. Panel (c) shows the width of the confidence intervals (2*1.96*standard error) for different levels of age and exposure.
Figure 5—figure supplement 2.
Figure 5—figure supplement 2.. Adjusting for cumulative exposure.
To further explore the independent contribution of age on the development of anti-parasite and anti-disease immunity, we fit models where, instead of adjusting for aEIR, we explicitly adjusted for the cumulative aEIR (cumEIR) as a metric of cumulative exposure. Cumulative aEIR was calculated as the product of age and aEIR. We fit the following models. Anti-parasite immunity.  Log10(Parasite density)ijk = f(ageijk, Log2cumEIRijk)+ui+γj. Anti-disease immunity. Temperatureijk=f(ageijk, Log2cumEIRijk, Log10Parasite densityijk)+ui+γj. Figure shows the results of these models. These results are similar to those presented in Figure 5, but adjusted for cumulative aEIR rather than for EIR. Left panel shows expected parasite densities (log 10, parasites/μL) after infection at different ages and levels of cumulative exposure. Right panel shows the expected fever thresholds (parasite densities required to develop a temperature 38°C or greater).
Figure 5—figure supplement 3.
Figure 5—figure supplement 3.. Adjusting for the number of infections in the past 3 months: Results of models.
quantifying anti-parasite (left) and anti-disease immunity (right). These results are similar to those presented in Figure 5, but adjusted for the number of P. falciparum positive visits in the last 3 months. Left panel shows expected parasite densities (log 10, parasites/μL) after infection at different ages and levels of exposure (aEIR). Right panel shows the expected fever thresholds (parasite densities required to develop a temperature 38°C or greater). Appendix 2 shows the estimated coefficients for these and additional models exploring the impact of recent exposure.
Figure 5—figure supplement 4.
Figure 5—figure supplement 4.. Adjusting for the probability of non-malaria fevers.
Since it is known that the probability of developing fever (not due to malaria) decreases with age, we performed sensitivity analyses where we adjusted for the probability of non-malaria fevers. Using data from instances when participants consulted for fever and/or were objectively febrile but were found to be smear negative, we fit a model to estimate the monthly probability of non-malaria fever as a function of age. This model was specified as: LogOddsnonmalariafeverij=f(ageij) where i is an index for individuals and j is an index for month. The model also included a random effect to account for clustering within individuals. Figure shows the distribution of estimated probabilities of non-malaria for the different age-groups in the dataset. We then used this model to predict the probability of non-malaria fever for each individual at each time point, and re-estimated our anti-parasite and anti-disease models adjusting for this probability (on the logit scale). See Figure 5—figure supplement 5.
Figure 5—figure supplement 5.
Figure 5—figure supplement 5.. Adjusting for the probability of non-malaria fevers.
Results of models quantifying anti-parasite (left) and anti-disease immunity (right). Left panel shows expected parasite densities (log 10, parasites/μL) after infection at different ages and levels of exposure (aEIR). Right panel shows the expected fever thresholds (parasite densities required to develop a temperature 38°C or greater). These results are similar to those presented in Figure 5, but adjusting for the probability of non-malaria fever experienced by the different age-groups. See Figure 5—figure supplement 4.
Figure 5—figure supplement 6.
Figure 5—figure supplement 6.. Adjusting for the probability of observation.
Since the study design only included 'routine' samples collections every 3 months, it is likely that we missed several events of asymptomatic infection, particularly in the moderate and high transmission settings. Moreover, since infections were detected using microscopy, we lack knowledge about the complexity of each infection. To account for the lower probability of observing asymptomatic, as compared to symptomatic infections, we conducted sensitivity analyses where we up-weighted asymptomatic observations to account for potentially missed infections. Weights were calculated as the expected number of infections experienced by the participant during the period of time comprised between the last visit (where the participant was tested for parasites) and the current visit. The expected number of infections (M) were derived using estimates of the average daily individual hazard of infection (λ) as. M=D(1-exp-λ) where D is the number of days between the last visit when the participant was tested and (1-exp-λ) is the daily probability of infection. Figure shows results of models quantifying anti-parasite (left) and anti-disease immunity (right). Left panel shows expected parasite densities (log 10, parasites/μL) after infection at different ages and levels of exposure (aEIR). Right panel shows the expected fever thresholds (parasite densities required to develop a temperature 38°C or greater). These results are similar to those presented in Figure 5, but weighted for the probability of observation of asymptomatic infections.
Figure 5—figure supplement 7.
Figure 5—figure supplement 7.. Limiting the analysis to ‘incident’ infections.
Since we don’t have data on the genotypes of parasites, it is also possible that consecutive asymptomatic infections represent persistent, rather than new infections. To assess the potential impact of including these persistent infections in the analyses, we conducted sensitivity analyses where we limited the dataset to ‘incident’ infections. We considered the following samples as belonging to ‘Incident’ infections. All instances of symptomatic malaria infection not preceded by an asymptomatic infection within 10 days. All instances of asymptomatic infection when. The previous blood smear (from routine or active surveillance) was negative. The previous blood smear (from routine surveillance) was positive, but the participant received malaria treatment between that visit and the current one. Limiting the analyses to incident infections implied excluding 700/1431 instances of asymptomatic parasitemia. Figure shows results of models quantifying anti-parasite (left) and anti-disease immunity (right). Left panel shows expected parasite densities (log 10, parasites/μL) after infection at different ages and levels of exposure (aEIR). Right panel shows the expected fever thresholds (parasite densities required to develop a temperature 38°C or greater). These results are similar to those presented in Figure 5, but limiting the dataset to ‘incident’ infections.
Figure 5—figure supplement 8.
Figure 5—figure supplement 8.. Limiting the analysis to individuals without the sickle hemoglobin mutation (β globin E6V).
This analysis only included data from 618/773 individuals without the sickle hemoglobin mutation (β globin E6V). Figure shows results of models quantifying anti-parasite (left) and anti-disease immunity (right). Left panel shows expected parasite densities (log 10, parasites/μL) after infection at different ages and levels of exposure (aEIR). Right panel shows the expected fever thresholds (parasite densities required to develop a temperature 38°C or greater). These results are similar to those presented in Figure 5, but limited to children without the sickle hemoglobin mutation (β globin E6V).
Figure 5—figure supplement 9.
Figure 5—figure supplement 9.. Limiting the analysis to individuals from Tororo and Kanungu.
This analysis only included data from 634/773 individuals living in Tororo (Nagongera) and Kanungu (Kihihi). Figure shows results of models quantifying anti-parasite (left) and anti-disease immunity (right). Left panel shows expected parasite densities (log 10, parasites/μL) after infection at different ages and levels of exposure (aEIR). Right panel shows the expected fever thresholds (parasite densities required to develop a temperature 38°C or greater). These results are similar to those presented in Figure 5, but limited to children living in Tororo (Nagongera) and Kanungu (Kihihi).
Figure 5—figure supplement 10.
Figure 5—figure supplement 10.. Limiting the analysis to individuals living in settings with aEIR ≥5.
This analysis only included data from 554/773 individuals living in households with aEIR >5. Figure shows results of models quantifying anti-parasite (left) and anti-disease immunity (right). Left panel shows expected parasite densities (log 10, parasites/μL) after infection at different ages and levels of exposure (aEIR). Right panel shows the expected fever thresholds (parasite densities required to develop a temperature 38°C or greater). These results are similar to those presented in Figure 5, but limited to children living in households with aEIR ≥5.
Figure 5—figure supplement 11.
Figure 5—figure supplement 11.. Using a metric of aEIR that only includes prior observations.
In our main analysis, the metric of exposure used was the mean aEIR over the whole study period, as a proxy of the mean exposure that each participant has experienced over their life-time. We performed a sensitivity analysis where instead, aEIR is calculated as the mean aEIR up to each point in time. This analysis only excluded data from 344/5640 observations where we did not have prior exposure information. Figure shows results of models quantifying anti-parasite (left) and anti-disease immunity (right). Left panel shows expected parasite densities (log 10, parasites/μL) after infection at different ages and levels of exposure (aEIR). Right panel shows the expected fever thresholds (parasite densities required to develop a temperature 38°C or greater). These results are similar to those presented in Figure 5, but with an alternative metric of aEIR as described above.
Figure 5—figure supplement 12.
Figure 5—figure supplement 12.. Predicted individual trajectories of anti-disease and anti-parasite immunity.
Figure showing predicted individual trajectories in the development of anti-parasite (a) and anti-disease (b) for a sample of study participants from each study site. (a) shows expected parasite densities (parasites/μL, log 10) after infection at different ages and levels of exposure (aEIR). (b) shows the expected fever thresholds (parasite densities required to develop a temperature 38°C or greater). Each line represents the predicted trajectory, but the solid portion represents the period of time when the individual contributed to the dataset.
Figure 5—figure supplement 13.
Figure 5—figure supplement 13.. Model checks.
Residual plots and response vs. fitted plots for the best fitting models (Models AP4 and AD4).
Figure 6.
Figure 6.. Anti-disease immunity as a function of age and exposure.
Each panel shows how the expected objective temperature (°C) varies as a function of age and parasite density, for different transmission settings. (a) aEIR = 2; (b) aEIR = 10; (c) aEIR = 50; (d) aEIR = 200. Contours indicating the fever threshold (38°C) are also shown. Confidence bounds for these plots are presented in Figure 6—figure supplement 1.
Figure 6—figure supplement 1.
Figure 6—figure supplement 1.. Association between parasite density (parasites/μL) and objective temperature (in °C).
Figure 6—figure supplement 2.
Figure 6—figure supplement 2.. Confidence bounds of best fitting model quantifying anti-disease immunity.
Left panel (a) is equivalent to Figure 6c, and shows how the expected objective temperature (°C) varies as a function of age and parasite density, in a setting with aEIR of 50. Panel (b) shows the lower and upper 95% confidence bounds. Panel (c) shows the width of the confidence intervals (2*1.96*standard error) for different levels of age and exposure.
Figure 7.
Figure 7.. Results of the best model quantifying overall immunity against symptomatic malaria (upon microscopically detectable infection), as a function of age and exposure (aEIR).
Colors represent the expected probability of developing symptomatic malaria upon infection. Confidence bounds for these plots are presented in Figure 7—figure supplement 1.
Figure 7—figure supplement 1.
Figure 7—figure supplement 1.. Confidence bounds of best fitting model quantifying overall immunity against symptomatic malaria.
Left panel (a) is equivalent to Figure 7. Colors represent the expected probability of developing symptomatic malaria upon infection, as a function of age and exposure. Panel (b) shows the lower and upper 95% confidence bounds. Panel (c) shows the width of the confidence intervals (2*1.96*standard error), on the logit scale, for different levels of age and exposure.
Author response image 1.
Author response image 1.. Anti-parasite immunity.
Results of models quantifying anti-parasite immunity. Left panel shows result when models are fit to data from routine visits only (Active) while right panel shows result when models are fit to non-routine (Passive) visits.
Author response image 2.
Author response image 2.. Anti-disease immunity.
Results of models quantifying disease immunity when fit to data from non-routine (Passive) visits only. It is not possible to fit this model using only data from routine visits as the majority (~95%) of participants are asymptomatic during these visits and therefore it is not possible to estimate a fever threshold.
Author response image 3.
Author response image 3.. Trends in parasite densities recorded during symptomatic (passive surveillance) infections and routine (active surveillance) visits as a function of age; and trends in the objective temperature recorded during visits in which participants were found to have a parasite density between 50,000 parasites/μL and 200,000 parasites/μL, as a function of age Each point represents a measurement obtained during a study visit.
Dark squares represent the mean of the observations from each site.
Author response image 4.
Author response image 4.. Figure summarizing evidence for anti-parasite immunity (left) and anti-disease immunity (right) in the raw data.
Each plot shows, for specific ages and aEIRs, the mean parasite density (/μL) (left), and the mean objective temperature given a parasite density between 10,000 to 80,000 parasites/μL (right). Means and 95% confidence intervals are shown. This plot corresponds to Figure 4 in the manuscript but was created using raw data (rather than model output).
Author response image 5.
Author response image 5.. Results of models quantifying anti-parasite and anti-disease immunity.
This figure is similar to figure 5 in the manuscript, but with the analysis limited to 2238/5640 observations from Walukuba and Kihihi (thus the limited aEIR range).

Comment in

  • Age, exposure and immunity.
    White M, Watson J. White M, et al. Elife. 2018 Aug 21;7:e40150. doi: 10.7554/eLife.40150. Elife. 2018. PMID: 30129437 Free PMC article.

References

    1. Baird JK. Age-dependent characteristics of protection v. susceptibility to Plasmodium falciparum. Annals of Tropical Medicine and Parasitology. 1998;92:367–390. doi: 10.1080/00034989859366. - DOI - PubMed
    1. Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, Battle K, Moyes CL, Henry A, Eckhoff PA, Wenger EA, Briët O, Penny MA, Smith TA, Bennett A, Yukich J, Eisele TP, Griffin JT, Fergus CA, Lynch M, Lindgren F, Cohen JM, Murray CLJ, Smith DL, Hay SI, Cibulskis RE, Gething PW. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–211. doi: 10.1038/nature15535. - DOI - PMC - PubMed
    1. Bull PC, Lowe BS, Kortok M, Molyneux CS, Newbold CI, Marsh K. Parasite antigens on the infected red cell surface are targets for naturally acquired immunity to malaria. Nature Medicine. 1998;4:358–360. doi: 10.1038/nm0398-358. - DOI - PMC - PubMed
    1. Carneiro I, Roca-Feltrer A, Griffin JT, Smith L, Tanner M, Schellenberg JA, Greenwood B, Schellenberg D. Age-patterns of malaria vary with severity, transmission intensity and seasonality in sub-Saharan Africa: a systematic review and pooled analysis. PLoS One. 2010;5:e8988. doi: 10.1371/journal.pone.0008988. - DOI - PMC - PubMed
    1. Filipe JA, Riley EM, Drakeley CJ, Sutherland CJ, Ghani AC. Determination of the processes driving the acquisition of immunity to malaria using a mathematical transmission model. PLoS Computational Biology. 2007;3:e255. doi: 10.1371/journal.pcbi.0030255. - DOI - PMC - PubMed

Publication types