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

  1. Isabel Rodriguez-Barraquer  Is a corresponding author
  2. Emmanuel Arinaitwe
  3. Prasanna Jagannathan
  4. Moses R Kamya
  5. Phillip J Rosenthal
  6. John Rek
  7. Grant Dorsey
  8. Joaniter Nankabirwa
  9. Sarah G Staedke
  10. Maxwell Kilama
  11. Chris Drakeley
  12. Isaac Ssewanyana
  13. David L Smith
  14. Bryan Greenhouse
  1. San Francisco General Hospital, University of California, San Francisco, United States
  2. Infectious Diseases Research Collaboration, Uganda
  3. Stanford University, United States
  4. Makerere University College of Health Sciences, Uganda
  5. London School of Hygiene and Tropical Medicine, United Kingdom
  6. University of Washington, United States
12 figures, 7 tables and 1 additional file

Figures

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.

https://doi.org/10.7554/eLife.35832.004
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 …

https://doi.org/10.7554/eLife.35832.005
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 …

https://doi.org/10.7554/eLife.35832.006
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 …

https://doi.org/10.7554/eLife.35832.008
Figure 5 with 13 supplements
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 …

https://doi.org/10.7554/eLife.35832.009
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 …

https://doi.org/10.7554/eLife.35832.010
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 …

https://doi.org/10.7554/eLife.35832.011
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 …

https://doi.org/10.7554/eLife.35832.012
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. …

https://doi.org/10.7554/eLife.35832.013
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 …

https://doi.org/10.7554/eLife.35832.014
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 …

https://doi.org/10.7554/eLife.35832.015
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 …

https://doi.org/10.7554/eLife.35832.016
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 …

https://doi.org/10.7554/eLife.35832.017
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 …

https://doi.org/10.7554/eLife.35832.018
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 …

https://doi.org/10.7554/eLife.35832.019
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 …

https://doi.org/10.7554/eLife.35832.020
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 …

https://doi.org/10.7554/eLife.35832.021
Figure 5—figure supplement 13
Model checks.

Residual plots and response vs. fitted plots for the best fitting models (Models AP4 and AD4).

https://doi.org/10.7554/eLife.35832.022
Figure 6 with 2 supplements
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) …

https://doi.org/10.7554/eLife.35832.023
Figure 6—figure supplement 1
Association between parasite density (parasites/μL) and objective temperature (in °C).
https://doi.org/10.7554/eLife.35832.024
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 …

https://doi.org/10.7554/eLife.35832.025
Figure 7 with 1 supplement
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.

https://doi.org/10.7554/eLife.35832.026
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 …

https://doi.org/10.7554/eLife.35832.027
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 …

https://doi.org/10.7554/eLife.35832.037
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 …

https://doi.org/10.7554/eLife.35832.038
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.

https://doi.org/10.7554/eLife.35832.039
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 …

https://doi.org/10.7554/eLife.35832.040
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).

https://doi.org/10.7554/eLife.35832.041

Tables

Table 1
Characteristics of the study participants.
https://doi.org/10.7554/eLife.35832.003
CharacteristicNagongeraKihihiWalukuba
Number of households10610076
Number of children329305139
Female, n (%)151(46)150 (49)66 (47)
Mean age at enrollment, years (sd)4.4 (2.7)4.6 (2.6)4.3 (2.6)
Mean follow up time, months (range)23.5 (0, 38.8)24.4 (0.8, 38.8)22.1 (2.3, 3.9)
Symptomatic malaria
Symptomatic Malaria episodes, n24471555207
Median number of symptomatic malaria episodes/child, n (range)6 (0, 29)4 (0, 30)1 (0. 12)
Median incidence of symptomatic malaria episodes ppy (range)2.6 (0, 10)1.6 (0, 15.2)0.6 (0, 5.1)
Asymptomatic parasitemia
Asymptomatic parasitemia episodes, n955331145
Median number of asymptomatic parasitemia episode/child, n (range)2 (0, 12)0 (0, 11)1 (0, 10)
Median prevalence of asymptomatic parasitemia (range)0.12 (0.07–0.17)0.05 (0.02–0.10)0.07 (0.03–0.11)
Household malaria exposure
Household aEIR, median (range)51 (10–582)7.7 (3.6–47)2.1 (1.5–8.1)
Table 2
Results of linear models quantifying the association between age, aEIR and immunity outcomes.
https://doi.org/10.7554/eLife.35832.007
All dataaEIR ≥ 5 (n = 5047)aEIR < 5 (n = 593)
Anti-parasite immunityFold change in parasite density (95% CI)
Age (years)0.78 (0.77, 0.79)0.76 (0.75, 0.77)0.87 (0.83, 0.90)
Log2 aEIR0.82 (0.79, 0.84)0.73 (0.69, 0.77)1.92 (1.69, 2.15)
Anti-disease immunity*Change in objective temperature C (95% CI)
Age (years)−0.07 (–0.06, –0.08)−0.08 (–0.07, –0.1)−0.04 (–0.07, –0.01)
Log2 aEIR−0.02 (–0.04, 0.0)−0.07 (–0.05, –0.1)0.27 (0.11, 0.44)
Overall immunity against symptomatic malariaOdds ratio of symptomatic disease (95% CI)
Age (years)0.78 (0.75, 0.82)0.77 (0.74, 0.80)0.90 (0.83, 0.99)
Log2 aEIR0.91 (0.74, 1.13)0.62 (0.48, 0.80)3.83 (1.39, 10.6)
  1. *Model adjusted as well for Log parasite density.

Appendix 1—table 1
Anti-parasite immunity Model fit.
https://doi.org/10.7554/eLife.35832.030
ModelDF% Deviance ExplainedAIC
AP1418.733.214052
AP2414.233.114048
AP34083314046
AP44043314032
Appendix 1—table 2
Anti-disease immunity Model fit.
https://doi.org/10.7554/eLife.35832.031
ModelDF% Deviance ExplainedAIC
AD131537.714946
AD228539.014774
AD331438.314894
AD429639.414758
Appendix 1—table 3
Overall immunity against symptomatic malaria Model fit.
https://doi.org/10.7554/eLife.35832.032
ModelDF% Deviance ExplainedAIC
SM1399.8285382
SM2385.3285373
SM3385.527.85380
SM4362.627.35369
Appendix 2—table 1
Anti-parasite immunity.
https://doi.org/10.7554/eLife.35832.034
Time windowFunctional formCoefficient (95% CI)p-value
3 monthsLinear−0.003 (−0.03–0.03)0.83
6 monthsLinear0.003 (−0.02–0.02)0.75
6 monthsSmooth-0.77
Appendix 2—table 2
Anti-disease immunity.
https://doi.org/10.7554/eLife.35832.035
Time windowFunctional formCoefficient (95% CI)p-value
3 monthsLinear−0.01 (−0.04–0.02)0.38
6 monthsLinear−0.01 (−0.03–0.01)0.15
6 monthsSmooth-0.15

Additional files

Download links