Potential public health impact of RTS,S malaria candidate vaccine in sub-Saharan Africa: a modelling study
- Christophe J. Sauboin1Email authorView ORCID ID profile,
- Laure-Anne Van Bellinghen†2,
- Nicolas Van De Velde†1 and
- Ilse Van Vlaenderen†2
© Sauboin et al. 2015
Received: 24 June 2015
Accepted: 9 December 2015
Published: 23 December 2015
Adding malaria vaccination to existing interventions could help to reduce the health burden due to malaria. This study modelled the potential public health impact of the RTS,S candidate malaria vaccine in 42 malaria-endemic countries in sub-Saharan Africa.
An individual-based Markov cohort model was constructed with three categories of malaria transmission intensity and six successive malaria immunity levels. The cycle time was 5 days. Vaccination was assumed to reduce the risk of infection, with no other effects. Vaccine efficacy was assumed to wane exponentially over time. Malaria incidence and vaccine efficacy data were taken from a Phase III trial of the RTS,S vaccine with 18 months of follow-up (NCT00866619). The model was calibrated to reproduce the malaria incidence in the control arm of the trial in each transmission category and published age distribution data. Individual-level heterogeneity in malaria exposure and vaccine protection was accounted for. Parameter uncertainty and variability were captured by using stochastic model transitions. The model followed a cohort from birth to 10 years of age without malaria vaccination, or with RTS,S malaria vaccination administered at age 6, 10 and 14 weeks or at age 6, 7-and-a-half and 9 months. Median and 95 % confidence intervals were calculated for the number of clinical malaria cases, severe cases, malaria hospitalizations and malaria deaths expected to be averted by each vaccination strategy. Univariate sensitivity analysis was conducted by varying the values of key input parameters.
Vaccination assuming the coverage of diphtheria-tetanus-pertussis (DTP3) at age 6, 10 and 14 weeks is estimated to avert over five million clinical malaria cases, 119,000 severe malaria cases, 98,600 malaria hospitalizations and 31,000 malaria deaths in the 42 countries over the 10-year period. Vaccination at age 6, 7-and-a-half and 9 months with 75 % of DTP3 coverage is estimated to avert almost 12.5 million clinical malaria cases, 250,000 severe malaria cases, 208,000 malaria hospitalizations and 65,400 malaria deaths in the 42 countries. Univariate sensitivity analysis indicated that for both vaccination strategies, the parameters with the largest impact on the malaria mortality estimates were waning of vaccine efficacy and malaria case-fatality rate.
Addition of RTS,S malaria vaccination to existing malaria interventions is estimated to reduce substantially the incidence of clinical malaria, severe malaria, malaria hospitalizations and malaria deaths across 42 countries in sub-Saharan Africa.
KeywordsMalaria Vaccination RTS,S candidate vaccine Public health
The health burden attributable to malaria is immense. The World Health Organization (WHO) estimated that in 2013 there were 163 million malaria cases and 528,000 malaria deaths in the WHO Africa Region . In the updated Global Burden of Disease study, malaria was the second leading cause of disability-adjusted life-years (DALYs) in sub-Saharan Africa and the first cause in western Africa in 2013 . Most malaria deaths occur in children aged <5 years (estimated 437,000 malaria deaths in WHO Africa Region in 2013) . Among children aged <5 years (excluding infants in the first month of life), malaria was the leading cause of death in Africa in 2010, accounting for 15 % of deaths .
Current preventive interventions in malaria aim to reduce malaria transmission by insect vectors (insecticide spraying and the use of insecticide-treated bed nets), or to reduce the disease burden by prophylactic treatment of defined population groups such as pregnant women, infants or children [4, 5]. The main treatment for malaria is artemisinin-based combination therapy . Decreases in estimated malaria case incidence (34 % decrease between 2000 and 2013) and estimated malaria mortality rate in children aged <5 years (58 % decrease between 2000 and 2013) in the WHO Africa region have been reported, and probably reflect the combined effects of malaria intervention programmes, increased urbanization and overall economic development . However, current malaria interventions have some limitations. Insecticide resistance has been reported in 49/63 reporting countries worldwide since 2010, and artemisinin-resistant Plasmodium falciparum malaria has been reported in five countries in Southeast Asia . Insecticide-treated nets need to be routinely distributed and replaced when they wear out in order to maintain coverage and therefore effectiveness. A slower growth in insecticide-treated bed net coverage in 2012–2013 in sub-Saharan Africa was associated with a slower decline in malaria mortality between 2011 and 2013 .
A malaria vaccine could be a valuable addition to existing malaria control interventions. The RTS,S malaria candidate vaccine is a pre-erythrocytic vaccine, directed against a protein from the sporozoite stage of the P. falciparum malaria parasite. It aims to trigger an immune defence to prevent the parasites infecting the liver and red blood cells. The RTS,S malaria candidate vaccine is administered in three doses given at one-month intervals, with a possible booster dose given 18 months after the third dose. A double-blind randomized controlled Phase III trial in seven countries in Africa enrolled 15,460 infants aged 6–12 weeks and children aged 5–17 months between March 2009 and January 2011, and randomized them in a 1:1:1 ratio to receive the RTS,S candidate vaccine with a booster, without a booster, or a non-malaria comparator vaccine [7–9]. Results with 18 months of follow-up from the third dose of vaccine were published in 2014, combining the first two groups and comparing them with the third (control) group . In children, vaccination reduced clinical malaria cases by 46 % [95 % confidence interval (CI) 42, 50 %], severe malaria cases by 34 % (95 % CI 15, 48 %) and malaria hospitalizations by 41 % (95 % CI 30, 50 %). In infants, the reduction in clinical malaria cases was 27 % (95 % CI 20, 32 %), with no significant protection against severe malaria or malaria hospitalization . The highest number of malaria episodes averted was seen in areas with greatest malaria incidence .
The objective of the present analysis is to apply these 18-month clinical trial results to estimate the potential public health impact of the RTS,S candidate malaria vaccine on the number of clinical malaria cases, severe malaria cases, hospitalizations and deaths in malaria-endemic countries in sub-Saharan Africa using a cohort model.
All episodes of clinical and severe malaria occurring in the trial were reported;
There are six levels of acquired immunity to clinical malaria depending on the previous number of infections;
A fixed percentage of severe infections result in death;
The effect of the vaccine is only to reduce the risk of infection, with no additional effect on the risk of developing clinical or severe disease;
Vaccine efficacy wanes over time by exponential decay;
No seasonality is included. Only one site (Nanoro in Burkina Faso) of the 11 trial sites exhibits a strongly seasonal pattern  which limits the availability of data to calibrate the model.
Incidence of clinical and severe malaria in the control group of the Phase III trial
Low (PP ≤ 5 %)
Moderate (PP 5–40 %)
High (PP > 40 %)
Kilifi (Kenya), Korogwe (Tanzania), Manhiça (Mozambique)
Bagamoyo (Tanzania), Lambarene (Gabon), Agogo (Ghana), Kintampo (Ghana), Lilongwe (Malawi)
Siaya (Kenya), Kombewa (Kenya), Nanoro (Burkina Faso)
Control group incidence (95 % CI) of clinical malaria per child per year in each 6–month follow-up period for infants receiving the first vaccine dose at age 6–12 weeks
Control group incidence (95 % CI) of clinical malaria per child per year in each 6–month follow-up period for children receiving the first vaccine dose at age 5–17 months
Incidence of severe malaria per 1000 children (95 % CI) per year across the whole follow-up period
Infants receiving the first vaccine dose at age 6–12 weeks
Children receiving the first vaccine dose at age 5–17 months
The model was calibrated to reproduce simultaneously the incidence of malaria observed in the control arm of the Phase III trial for both age groups (infants and children) in each transmission category and each six-month follow-up period, and the age distribution up to the age of 5 years in moderate and high transmission areas with no marked seasonality reported by Carneiro et al. . Since different definitions of transmission intensity categories were used, high transmission intensity data from the model were compared against moderate transmission intensity data from Carneiro et al., and moderate transmission data from the model were compared against an average calculated from moderate and low transmission data from Carneiro et al. The age distribution in low transmission areas was not included because the definition of low transmission used by Carneiro et al. (PP ≤ 25 %)  was too different from the definition used in the model (PP ≤ 5 %).
Infection probability for each transmission category (q);
Probability of asymptomatic infection with each successive infection (a1–a6);
Susceptibility factor (s) representing an increase in susceptibility with age to account for higher probability of mosquito bites with larger body surface, is derived from a negative exponential equation of the following type: 1 − exp(−k*n), where ‘n’ is the age of the child expressed as the number of model cycles (5 days per cycle) and ‘k’ is derived by calibration;
Waning of immunity (w) implemented as a transition from Si+1 to Si;
Percentage of individuals gaining full immunity after six infections (rimm).
Calibration was performed by the least square method using Matlab (version R2013b) minimization tools (interior-point algorithm). To ensure that immunity always increased with successive infections, a constraint was set so that ai was always less than or equal to ai+1 for values of i from 1 to 5.
Capturing uncertainty and heterogeneity
A main source of uncertainty comes from the observation of incidence rates and vaccine efficacy estimates of the Phase III trial. Although the Phase III trial is powered for clinical malaria incidence and severe malaria incidence overall, larger variations appear when stratifying data by malaria transmission level and periods of follow-up. In order to capture this uncertainty, model parameter distributions are generated by repeating the calibration process with several incidence rates. These incidence rates are randomly drawn from a log-normal distribution generated using median values and 95 % CI for each 6-month follow-up period and for each age-group.
Point estimates, 95 % confidence intervals and distributions of fitted parameters
95 % CI
q, low transmission
1.224 e−3, 2.217 e−3
q, moderate transmission
19.222 e−3, 37.682 e−3
q, high transmission
87.972 e−3, 272.147 e−3
Probability of full immunity, r imm
1.93 %, 2.17 %
Age-related susceptibility factor, k
1.02 e−2, 3.63 e−2
Waning of acquired immunity, w
Probability of asymptomatic infection
0.29, 17.79 %
7.74, 43.62 %
7.94, 43.82 %
8.37, 43.93 %
15.65, 44.12 %
25.87, 61.82 %
Percentage of clinical cases that become severe with good access to care
1.19, 5,46 %
1.19, 3.04 %
1.19, 3.03 %
1.19, 2.92 %
1.19, 2.92 %
0.96, 2.12 %
Vaccine efficacy in infants aged 6–12 weeks
35.53, 41.11 %
Half-life in infants aged 6–12 weeks (months)
Vaccine efficacy in children aged 5–17 months
56.5, 59.6 %
Half-life in children aged 5–17 months (months)
The distribution of the number malaria episodes per child in the Phase III trial in each transmission category indicated an individual-level heterogeneity in the risk of infection, as the proportion of children with no malaria episodes was higher in moderate and high transmission sites than would be expected with a homogeneous risk of infection. This heterogeneity in exposure was captured by assuming an individual-level variation around a mean exposure level corresponding to each of the three transmission levels included in the model.
Severe malaria cases and hospitalizations
Data on the incidence of severe cases (defined as a parasite density of >5000) per 1000 person-years for each transmission category were obtained from the Phase III trial  (Table 1). A calibration step fitted the model parameters for the percentage of clinical cases becoming severe at each of the six modelled successive infections (f1–f6) to the incidence of severe cases observed in the Phase III trial. A constraint of diminishing percentages of clinical cases becoming severe with increasing number of previous infections (fi+1 < fi) is imposed to reflect the acquisition of immunity against severe malaria. This constraint results in a faster acquisition of immunity against severe disease than immunity against clinical disease.
The ratio of severe malaria cases to hospitalized malaria cases (defined as a parasite density of >5000 and excluding planned admissions and trauma) was calculated for each transmission category from the Phase III trial. However, the level of care available to patients in the Phase III trial was very high, with all malaria cases having access to high levels of treatment resulting in very low malaria mortality. This is not representative of the real-life situation in sub-Saharan Africa, where not all malaria cases will have access to prompt treatment and/or hospital care. Thus, compared with real-world situations, the Phase III trial would be expected to have a higher percentage of severe cases receiving hospital care, and a lower percentage of cases that become severe. The model thus required adjustment to reflect limited access to treatment for uncomplicated cases, and limited access to hospitalization for severe cases, as could be expected in real-life settings.
The proportion of cases in real-life settings with access to artemisinin-based combination treatment was taken from a published study for Ghana, Kenya and Tanzania , with the average of these countries (fixed value of 54 %) used for other countries. For the country-level estimates, the percentage of severe cases that were hospitalized in real-life settings was assumed to be the same as the percentage with access to artemisinin-based combination treatment .
The risk of developing a severe episode from a treated uncomplicated episode was calibrated from the trial data as described above. The relative risk of developing a severe episode from an untreated uncomplicated episode was calculated using the following formula:
Severe malaria in population = (Proportion untreated × relative risk + proportion treated) × severe malaria in model
Severe malaria in model = Estimated incidence of severe malaria in children aged <5 years predicted by model;
Severe malaria in population = Incidence of severe malaria in hospitals in children aged 0–9 years , corrected for the proportion of severe episodes in children aged <5 years using the age distribution from Carneiro et al.  and for the proportion of severe cases who are not hospitalized .
Values for fixed parameters
Vaccine efficacy after 1 dose
Vaccine efficacy after 2 doses
One quarter of the vaccine efficacy of a full course (see Table 2)
Maternal protection, % of infants protected
Waning of maternal protection, wm
3 months half-life
Similar to previous models
Probability of recovery from clinical disease, r
Half recover after 15 days
Assumption for a context of good access to care
Access to ACT treatment
Ghana: 55 %, Kenya: 45 %, Tanzania: 62 % and country surveys for 22 countries
Average of 31 % for the remaining 17 countries
Sicuri  and Demographic and Health Surveys and Malaria Indicators Surveys
% of severe cases hospitalized
Assumed same as access to ACT in public health facilities in Ghana, Kenya and Tanzania
Relative risk for an untreated uncomplicated episode becoming severe (compared with a treated uncomplicated episode)
1.84 (95 % CI 1.68, 2.01)
Calculated from the modelled number of severe cases, the % of severe cases hospitalized, and access to treatment 
% of severe malaria cases resulting in sequelae
1.7 % (range 0.85–2.54 %)
Case-fatality rate, CFR (% of severe cases)
13.6 % (95 % CI 8.4, 18.8 %)
3 × treated cases
40.8 % (95 % CI 25.2, 56.4 %)
Long-term sequelae and deaths
The proportion of severe cases with long-term sequelae was calculated from published data on the proportion of severe cases that are categorized as cerebral malaria  and the percentage of cerebral malaria cases with sequelae . The mean age of children presenting with cerebral malaria is about 3 years compared to 1.8 for those presenting with malarial anaemia . Malaria mortality was estimated as a fixed percentage of severe cases (case-fatality rate), with values for hospitalized and non-hospitalized severe cases obtained from Thwing et al.  (Table 3).
Calibration of vaccine effect
The vaccine effect was modelled by applying a reduction in the force of infection (i.e., a reduction in the probability of moving from state S to state I in each successive infection). This type of vaccine effect was considered in order to reproduce the mechanism of action of RTS,S targeting the parasite at the pre-erythrocytic stage. It was assumed that vaccine efficacy would wane over time following an exponential decay function.
Point estimates and 95 % CI for vaccine efficacy were calculated from the Phase III clinical trial data for each of the three transmission categories and each of the two vaccination age groups using the exact Fisher formula .
The point estimates for the initial reduction in the force of infection post-dose 3 (or vaccine efficacy against infection) and decay parameter in the model were fitted to the vaccine efficacy observed in the Phase III trial simultaneously for each transmission level, each six-month follow-up period and the entire 18-month follow-up period. The upper and lower bounds of the 95 % CI for vaccine efficacy against infection were obtained by fitting simultaneously the upper and lower values for each six-month follow-up period and the 18-month follow-up period. For the vaccine waning parameter, the upper (lower) bound was derived by combining the upper (lower) bound, point estimates and lower (upper) bound of vaccine efficacy from the three successive six-month periods and the 18-month point estimate for vaccine efficacy. The fitted parameters for vaccine efficacy and the half-life of the waning function for each vaccination age group are shown in Table 2. Vaccine efficacy rates after one or two doses were set at fixed values and are shown in Table 3.
Estimation of country-level impact
Estimates of the malaria burden in the absence of vaccination and the potential impact of vaccination were made for 42 countries in sub-Saharan Africa. For each country, the size of the annual birth cohort and all-cause child mortality rates were taken from United Nations data [24, 25]. The birth cohort was then allocated between the three transmission intensity categories using 2010 MAP data for that country . In some countries a fraction of the birth cohort was considered not at risk or with an unstable risk of malaria. No impact of vaccination was considered in that population. The percentage of cases with access to treatment was taken from a published survey  for Ghana, Kenya and Tanzania and Demographic and Health and Malaria Surveys in 22 additional countries representing about 80 % of the birth cohorts of the 42 countries. The average of 31 % access to anti-malarial treatment from the 25 countries in these surveys was used for remaining countries. The percentage of severe cases hospitalized and the case-fatality rate were set at fixed values for all countries (see Table 3). Running the model with these input data produced estimates of the numbers of clinical malaria cases, severe malaria cases, malaria hospitalizations and malaria deaths for each country. The proportion of the population in each malaria transmission category was assumed to remain stable over the projected period of 10 years. Data are presented in this paper for all 42 countries combined and for 8 individual countries: two countries with a large malaria burden (Democratic Republic of Congo and Nigeria accounting for 39 % of global malaria cases and 34 % of global malaria deaths ) and four of the countries in the Phase III trial (Burkina Faso, Ghana, Kenya, and Tanzania) and two additional countries (Senegal and Uganda). Details for each country are presented in Additional file 3: Table S1.
Without malaria vaccination;
With RTS,S malaria vaccination with doses given at age 6, 10 and 14 weeks (corresponding to the age group vaccinated at age 6–12 weeks in the trial, and to addition of the malaria vaccine to the Expanded Programme on Immunization (EPI) schedule), with coverage assumed to be the same as the diphtheria-tetanus-pertussis (DTP3) vaccination coverage from each country;
With RTS,S malaria vaccination with doses given at age 6, 7-and-a-half and 9 months (corresponding to the age group vaccinated at age 5–17 months in the trial).
The vaccine efficacy in the 5–17 months age groups was used for children receiving the first dose at 6 months of age. Efficacy did not change between children receiving the first dose before or after 12 months of age. The coverage was assumed to reach 75 % of DTP3 vaccination coverage from each country to account for the difficulty to reach children at a later age compared to DTP schedule. The modelled outcomes were the numbers of clinical malaria cases, severe malaria cases, malaria hospitalizations and malaria deaths expected in the cohort from birth to the age of 10 years. This process was repeated 50 times with and without vaccination and the values compared to estimate the number of events averted by vaccination for each outcome. At each year, the median value and the 2.5th and 97.5th percentile values (corresponding to the 95 % CI) were recorded for each event averted.
Univariate sensitivity analysis was conducted by running the model with different values for transmission, access to treatment, duration of vaccine protection, vaccine coverage, risk of mortality in the community, relative risk of severe malaria for untreated cases, and case-fatality rate. The other parameters are held constant with median values provided in Tables 2 and 3.
An additional scenario analysis was conducted to compare the model with estimates from the WHO World Malaria Report 2012 .
Model calibration and validation
Model validation, clinical and severe malaria in vaccine and control arms
Age group 6–12 weeks
Severe episodes (% of clinical episodes)
Age group 5–17 months
Severe episodes (% of clinical episodes)
Estimated public health impact of adding RTS,S vaccination to existing interventions
Estimated country-level public health impact of adding RTS,S vaccination
Estimated country-level impact of RTS,S vaccination
Number of vaccinees
Median (95 % CI)
Clinical malaria cases averted
Severe malaria cases averted
Malaria hospitalizations averted
Malaria deaths averted
Vaccination at 6, 10 and 14 weeks
234,951 (2090, 477,517)
4182 (−6846, 14,515)
3474 (−5688, 12,059)
1092 (−1788, 3,790)
Democratic Republic of Congo
590,733 (77,960, 1,114,625)
13, 943 (−22,087, 41,080)
11, 583 (−18,349, 34,128)
3641 (−5767, 10,727)
180,785 (17,756, 344,188)
4080 (−8121, 11,908)
3389 (−6747, 9893)
1065 (−2121, 3109)
96,415 (12,032, 189,548)
2805 (−6042, 9841)
2330 (−5020, 8176)
732 (−1578, 2570)
1,137,914 (106,434, 2,194,118)
25,223 (−41,239, 77,330)
20,955 (−34,261, 64,243)
6586 (−10,768, 20,192)
255,645 (37,085, 478,220)
6484 (−13,984, 19,374)
5387 (−11,618, 16,096)
1693 (−3,652, 5059)
36,048 (7275, 67,343)
1477 (−3337, 4803)
1227 (−2772, 3990)
386 (−871, 1254)
411,047 (11,695, 819,896)
7480 (−14,269, 24,160)
6214 (−11,854, 20,071)
1953 (−3726, 6309)
42 countries in sub-Saharan Africa
5,004,489 (474,867, 9,684,160)
118,736 (−214,795, 372,995)
98,642 (−178,445, 309,873)
31,004 (−56,087, 97,396)
Vaccination at 6, 7.5 and 9 months
562,795 (327,910, 669,793)
9985 (−54, 19,791)
8295 (−45, 16,441)
2607 (−14, 5168)
Democratic Republic of Congo
1,414,342 (923,277, 1,657,736)
28,922 (1270, 54,286)
24,028 (1,055, 45,099)
7552 (332, 14,175)
468,774 (323,528, 549,868)
8254 (82, 15,470)
6857 (69, 12,852)
2155 (22, 4040)
260,446 (189,489, 307,495)
5,641 (−636, 11,453)
4,686 (−528, 9,515)
1473 (−166, 2991)
2,749,346 (1,748,642, 3,234,411)
54,546 (1758, 103,606)
45,315 (1460, 86,073)
14,243 (459, 27,054)
693,506 (515,124, 811,055)
12,453 (−451, 23,670)
10,346 (−375, 19,665)
3252 (−118, 6181)
101,544 (80,459, 118,862)
2741 (−315, 5441)
2277 (−262, 4520)
716 (−82, 1421)
1,038,791 (657,896, 1,229,031)
16,671 (−72, 32,293)
13,850 (−60, 26,828)
4353 (−19, 8432)
42 countries in sub-Saharan Africa
12,484,309 (8,242,412, 14,696,728)
250,448 (−433, 481,819)
208,064 (−359, 400,280)
65,397 (−113, 125,813)
Estimated country-level impact of RTS,S vaccination per 100,000 vaccinees
Median (95 % CI)
Clinical malaria cases averted
Severe malaria cases averted
Malaria hospitalizations averted
Malaria deaths averted
Vaccination at 6, 10 and 14 weeks
34,861 (310, 70,852)
621 (−1016, 2154)
516 (−844, 1789)
162 (−265, 562)
Democratic Republic of Congo
27,132 (3581, 51,194)
640 (−1014, 1887)
532 (−843, 1567)
167 (−265, 493)
24,435 (2400, 46,521)
551 (−1098, 1610)
458 (−912, 1337)
144 (−287, 420)
6773 (845, 13,316)
197 (−424, 691)
164 (−353, 574)
51 (−111, 181)
31,252 (2923, 60,259)
693 (−1133, 2124)
576 (−941, 1764)
181 (−296, 555)
13,604 (1973, 25,448)
345 (−744, 1031)
287 (−618, 857)
90 (−194, 269)
7671 (1548, 14,331)
314 (−710, 1022)
261 (−590, 849)
82 (−185, 267)
28,104 (800, 56,058)
511 (−976, 1652)
425 (−811, 1372)
134 (−255, 431)
42 countries in sub-Saharan Africa
19,851 (1884, 38,413)
471 (−852, 1480)
391 (−708, 1229)
123 (−222, 386)
Vaccination at 6, 7.5 and 9 months
115,124 (67,076, 137,011)
2042 (−11, 4048)
1697 (−9, 3363)
533 (−3, 1057)
Democratic Republic of Congo
91,161 (59,510, 106,849)
1864 (82, 3499)
1549 (68, 2907)
487 (21, 914)
86,586 (59,758, 101,564)
1524 (15, 2857)
1267 (13, 2374)
398 (4, 746)
25,005 (18,192, 29,522)
542 (−61, 1100)
450 (−51, 914)
141 (−16, 287)
104,755 (66,626, 123,236)
2078 (67, 3948)
1727 (56, 3280)
543 (17, 1031)
50,152 (37,252, 58,652)
901 (−33, 1712)
748 (−27, 1422)
235 (−9, 447)
29,478 (23,357, 34,505)
796 (−92, 1580)
661 (−76, 1312)
208 (−24, 412)
96,898 (61,368, 114,643)
1555 (−7, 3012)
1292 (−6, 2502)
406 (−2, 787)
42 countries in sub-Saharan Africa
68,169 (45,007, 80,249)
1368 (−2, 2631)
1136 (−2, 2186)
357 (−1, 687)
Over the 42 African countries it is estimated that 2019 (−3652; 6341) cases of neurological sequelae would be averted with vaccination at the age of 6, 10 and 14 weeks and 4258 (−7; 8191) cases of neurological sequelae with vaccination at the age of 6, 7-and-a-half and 9 months.
This analysis applied 18-month follow-up data from a Phase III clinical trial to estimate the potential public health impact of the RTS,S candidate malaria vaccine in 42 malaria-endemic countries in sub-Saharan Africa using a cohort model. Two vaccination strategies were considered, with doses administered either at 6, 10 and 14 weeks (corresponding to adding the RTS,S vaccine to the EPI schedule) or at 6, 7-and-a-half and 9 months. Each strategy was compared with no vaccination to estimate the expected impact of RTS,S vaccination on averting clinical malaria cases, severe cases, malaria hospitalizations and malaria deaths in a birth cohort followed to the age of 10 years. Vaccination at age 6, 10 and 14 weeks would be expected to avert over five million cases of clinical malaria, 119,000 severe malaria cases, 98,600 malaria hospitalizations and 31,000 malaria deaths. Vaccination at age 6, 7-and-a-half and 9 months would be expected to avert almost 12.5 million cases of clinical malaria, 250,000 severe malaria cases, 208,000 malaria hospitalizations and 65,400 malaria deaths. The wide CIs reflect the combination of heterogeneity and uncertainty allowed for by the stochastic process in the model. The greatest uncertainty is around the impact on mortality.
The model is calibrated to the incidence of clinical malaria in the control arm of the Phase III trial based on the secondary case definition (parasite density of >0 and presence of fever). It may be possible that the secondary case definition would have captured cases of fever due to other illnesses and the malaria parasites observed were present as a co-morbidity rather than the primary cause. The primary case definition (parasite density of >5000 and presence of fever) would be more specific and lead to a reduction of estimates by 35 %.
Estimated malaria episodes from WHO World Malaria report may not sufficiently adjust for under-reporting, although a correction factor is used;
The model estimated malaria incidence based on reported proportion of children within categories of parasite prevalence from the MAP project. If this proportion varies over the year or the force of infection varies within a category, this is not captured by the model. Finally, the case-fatality rate for severe cases not treated in hospital is subject to a large uncertainty. An estimate from Thwing et al.  is applied in the model but this factor has a large impact on malaria mortality estimates (see Fig. 5a, b).
However it should be noted that clinical vaccine efficacy does not change with the primary or secondary case definition.
Another modelling estimate of 252 (95 % CI 171–353) million malaria episodes in sub-Saharan Africa for 2010 has been reported by Griffin et al. . These estimates include the total population and suggest the proportion of cases that are in under-fives varies from above 60 % at high transmission to below 20 %. Hay et al.  have reported estimates of 271 million of malaria cases due P. falciparum in Africa in 2007.
Obtaining a realistic estimate of the number of severe malaria cases and malaria hospitalizations is important as these outcomes contribute to understanding the health burden of malaria. However, available data are limited and generally relate to populations with access to treatment. For example, the Phase III trial was conducted in settings where all patients had access to a high level of care. It was found that about half of the patients hospitalized with malaria met the case definition for severe malaria (parasite density <5000), which is similar to the value reported in a study of hospitalized patients in Tanzania, where 4261 patients of the 9337 with positive blood slides had severe disease (46 %) . However, neither of these studies could provide data on the fraction of severe malaria cases that did not result in hospitalization, or on the risk of severe malaria developing in patients without access to treatment. The present analysis attempted to adjust for limited access to prompt and high-quality treatment for uncomplicated cases and limited access to hospitalization for severe cases.
Estimates of the potential impact of a pre-erythrocytic vaccine have been produced using other published malaria models developed by Imperial College London, the Swiss Tropical and Public Health Institute (Swiss TPH) and the Institute for Disease Modelling ([29–33]). These analyses have used a theoretical vaccine profile or calibrated the profile using results from a Phase II trial  and the 12-month follow-up data of the Phase III trial  of the RTS,S vaccine. Compared with the model presented in this study, the other models include the whole population (not just a birth cohort) and a dynamic malaria transmission process simulating vector population infectiousness and its interaction with the human host population. Therefore, the entomological inoculation rate is used as a measure of human host exposure, and can be affected by other interventions, such as insecticide-treated nets and indoor residual spraying. Those models include a mechanism of immunity acquisition acting on parasite density in individuals and reducing the probability of clinical disease, increasing tolerance to parasite and speeding-up parasite clearance. In the model described here, dynamic transmission is not included and no estimate for parasite density is made in individuals. Instead, parasite prevalence is used as a measure of exposure at the population level and a risk of infection is derived from trial incidence data. The immunity acquisition process against clinical disease depends on the number of previous infections in the individual. Treatment of malaria episodes is the only other intervention explicitly included in the model, because the effects of preventive vector-level interventions are assumed to be included in the estimates of parasite prevalence from MAP. Vaccine effect is introduced in a similar way to the other models, i.e., by reducing the initial risk of infection with no additional protection assumed. The definition of severe cases used in the present analysis is derived from the definition used in the Phase III trial, and differs from that used in an earlier analysis conducted with the Swiss TPH model . However, the incidence of severe cases in the present analysis was similar to that reported by the Swiss TPH model .
One of the strengths of the present model is that it is simpler than the dynamic models and can be populated with readily available data, such as parasite prevalence and country-specific birth cohort data. The cohort approach used for the model helps to directly compare the results with the trial data. The model can thus be easily adapted and populated for specific countries. It is also relatively straightforward for decision-makers who are not experts in epidemiological modelling to understand and use. The present analysis has fitted the model to data from the most recent and large multicountry, multicentre trial on malaria vaccination, using the data available from the trial for the 18 months follow-up period, and may therefore be more representative of the current malaria epidemiological situation than models based on studies conducted during a period of very high transmission with no artemisinin-based combination therapy or long-lasting insecticide-treated nets.
The model has some limitations. First, as it is not a dynamic model it cannot take account of herd immunity. However, herd immunity is likely to be limited because of the small fraction of the population targeted by vaccination, the limited duration of protection and hence the small fraction of the parasite reservoir in human hosts which could be reached. Third, the model does not distinguish between different types of severe malaria, such as anaemia or cerebral malaria, which may have different implications for outcomes and resource use, although it does take account of the long-term neurological sequelae that can result from cerebral malaria. Fourth, it is assumed that the current levels of malaria transmission would remain stable over the period of analysis, although additional usage of long-lasting insecticide-treated nets and access to treatment are likely to further decrease malaria burden in the coming years. Conversely, the growing threat of resistance to insecticide and first-line treatment could lead to increasing malaria burden over that period.
A further area of work is to include seasonality in the model to improve the calibration using data from Nanoro trial site and impact estimate in seasonal settings, reflecting factors such as changes in temperature and/or rainfall.
The uncertainty around impact estimates is larger for severe malaria events and malaria deaths due to the trial design and stochastic variability. This leads to having lower bounds below zero on impact estimates for these outcomes, although it is a marginal effect in the 6, 7-and-a-half and 9 months schedule. It should be noted that the approach applied to characterize the uncertainty around the decay of vaccine-induced protection could potentially lead to overestimation of its variability. Structural uncertainty was assessed mainly around the natural immunity acquisition mechanism. The calibration process has been performed for model structures involving less immunity stages showing poorer fit on age-distribution data and severe malaria events. To test more immunity stages a functional form of immunity acquisition should be assumed in order to limit the number of parameters. Alternatively, transitions for the waning of acquired immunity were included but did not improve the fit. Vaccination has a different mechanism of action than existing malaria interventions. Rather than targeting the vector or directly acting on parasites in the blood, vaccination reduces the risk of infection by increasing immunity to the malaria circumsporozoite protein. The use of interventions with complementary effects may allow greater malaria control than a single intervention alone.
Addition of RTS,S malaria vaccination to existing interventions with doses administered at either age 6, 10 and 14 weeks, or at age 6, 7-and-a-half and 9 months, would be expected to reduce substantially the incidence of clinical malaria, severe malaria, malaria hospitalizations and malaria deaths in children across 42 countries in sub-Saharan Africa, compared with no vaccination. There is a larger uncertainty for the estimates impact on of severe malaria cases and deaths.
Expanded Programme on Immunization
Malaria Atlas Project
World Health Organization
- Swiss TPH:
Swiss Tropical and Public Health Institute
CJS has provided substantial scientific input to the study, participating to the method selection and development, to the model elaboration and population, to the determination of inputs and the acquisition of data, to the statistical data analysis and the sensitivity analysis, to the elaboration of the study report and its critical review; and has written sections of the initial draft. LAVB and IVV have provided substantial scientific input to the study, critically reviewing the model inputs and the study report. NVDV has provided substantial scientific input to the study, participating to the determination of model settings and inputs. All authors read and approved the final manuscript.
This study was funded by GlaxoSmithKline Biologicals SA (internal study identifier: HO1414955), which was involved in all stages of the study conduct and analysis including: design and operation of the model; data inputs, analysis and interpretation; manuscript preparation, review and decision to submit for publication. GlaxoSmithKline Biologicals SA also funded all costs associated with the development and the publishing of the present manuscript. All authors had full access to the data and agreed with the submission of the publication. The authors would like to thank Lisa White (Mahidol-Oxford Tropical Medicine Research Unit, Centre for Tropical Medicine) for her participation to the method selection, the development of the model, and the determination of the inputs; Edouard Ledent (GSK Vaccines) for his input on the statistical analysis of the calibration, Carole Nadin (Fleetwith Ltd, on behalf of GSK Vaccines) for medical writing assistance and Fabien Debailleul (Business and Decision Life Sciences, on behalf of GSK Vaccines) for publication co-ordination.
All authors have completed the Unified Competing Interest form available at http://www.icmje.org/coi_disclosure.pdf (available on request from the corresponding author). CJS and NVDV are employees of the GSK group of companies and report ownership of stock options/restricted shares from GSK group of companies. IVV and LAVB report that GSK group of companies paid consulting fees to CHESS for its contribution in model development as well as consulting fees for other projects commissioned to CHESS by GSK group of companies.
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.
- WHO. World Malaria Report. Geneva: World Health Organization; 2014. http://www.who.int/malaria/publications/world_malaria_report_2014/en/. Accessed 16 December 2014.
- GBD 2013 DALYs and HALE Collaborators. Global, regional, and national disability-adjusted life years (DALYs) for 306 diseases and injuries and healthy life expectancy (HALE) for 188 countries, 1990–2013: quantifying the epidemiological transition. Lancet. 2015;386:2145–91.
- Liu L, Johnson HL, Cousens S, Perin J, Scott S, Lawn JE, et al. Global, regional, and national causes of child mortality: an updated systematic analysis for 2010 with time trends since 2000. Lancet. 2012;379:2151–61.PubMedView ArticleGoogle Scholar
- WHO. World Malaria Report. Geneva: World Health Organization; 2011. http://www.who.int/malaria/world_malaria_report_2011/en/index.html. Accessed 1 April 2012.
- WHO. WHO Policy Recommendation: seasonal malaria chemoprevention (SMC) for Plasmodium falciparum malaria control in highly seasonal transmission areas of the Sahel sub-region in Africa. Geneva: World Health Organization; 2012. http://www.who.int/malaria/publications/atoz/who_smc_policy_recommendation/en/index.html. Accessed 1 September 2012.
- WHO. Guidelines for the treatment of malaria, 2nd edn. Geneva: World Health Organization; 2010.
- Leach A, Vekemans J, Lievens M, Ofori-Anyinam O, Cahill C, Owusu-Agyei S, et al. Design of a phase III multicenter trial to evaluate the efficacy of the RTS, S/AS01 malaria vaccine in children across diverse transmission settings in Africa. Malar J. 2011;10:224.PubMedPubMed CentralView ArticleGoogle Scholar
- Lievens M, Aponte JJ, Williamson J, Mmbando B, Mohamed A, Bejon P, et al. Statistical methodology for the evaluation of vaccine efficacy in a phase III multi-centre trial of the RTS, S/AS01 malaria vaccine in African children. Malar J. 2011;10:222.PubMedPubMed CentralView ArticleGoogle Scholar
- Vekemans J, Marsh K, Greenwood B, Leach A, Kabore W, Soulanoudjingar S, et al. Assessment of severe malaria in a multicenter, phase III, RTS, S/AS01 malaria candidate vaccine trial: case definition, standardization of data collection and patient care. Malar J. 2011;10:221.PubMedPubMed CentralView ArticleGoogle Scholar
- The RTS,S Clinical Trials Partnership. Efficacy and safety of the RTS, S/AS01 malaria vaccine during 18 months after vaccination: a phase 3 randomized, controlled trial in children and young infants at 11 African sites. PLoS Med. 2014;11:e1001685.View ArticleGoogle Scholar
- Gething PW, Patil AP, Smith DL, Guerra CA, Elyazar IR, Johnston GL, et al. A new world malaria map: Plasmodium falciparum endemicity in 2010. Malar J. 2011;10:378.PubMedPubMed CentralView ArticleGoogle Scholar
- Roset Bahmanyar E. Plasmodium falciparum parasite prevalence in sub-Saharan Africa: an epidemiological assessment [abstract]. In: American Society of Tropical Medicine and Hygiene, 61st Annual Meeting. 2012.
- Ulm K. A simple method to calculate the confidence interval of a standardized mortality ratio (SMR). Am J Epidemiol. 1990;131:373–5.PubMedGoogle Scholar
- Carneiro I, Roca-Feltrer A, Griffin JT, Smith L, Tanner M, Schellenberg JA, et al. 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Sicuri E, Vieta A, Lindner L, Constenla D, Sauboin C. The economic costs of malaria in children in three sub-Saharan countries: Ghana, Tanzania and Kenya. Malar J. 2013;12:307.PubMedPubMed CentralView ArticleGoogle Scholar
- Cohen JM, Woolsey AM, Sabot OJ, Gething PW, Tatem AJ, Moonen B. Public health. Optimizing investments in malaria treatment and diagnosis. Science. 2012;338:612–4.PubMedView ArticleGoogle Scholar
- Cohen JM, Woolsey AM, Sabot OJ, Gething PW, Tatem AJ, Moonen B. Public health. Optimizing investments in malaria treatment and diagnosis (supplementary materials). 2015. http://www.sciencemag.org/content/338/6107/612.short. Accessed 4 June 2015.
- Marsh K, Snow RW. Malaria transmission and morbidity. Parassitologia. 1999;41:241–6.PubMedGoogle Scholar
- Tediosi F, Maire N, Smith T, Hutton G, Utzinger J, Ross A, et al. An approach to model the costs and effects of case management of Plasmodium falciparum malaria in sub-saharan Africa. Am J Trop Med Hyg. 2006;75:90–103.PubMedGoogle Scholar
- 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 severe Plasmodium falciparum malaria. JAMA. 2005;293:1461–70.PubMedView ArticleGoogle Scholar
- World Health Organization. Severe malaria. Trop Med Int Health. 2014;19(suppl 1):7–131.Google Scholar
- Thwing J, Eisele TP, Steketee RW. Protective efficacy of malaria case management and intermittent preventive treatment for preventing malaria mortality in children: a systematic review for the Lives Saved Tool. BMC Public Health. 2011;11(suppl 3):S14.PubMedPubMed CentralView ArticleGoogle Scholar
- Sahai H, Khurshid A. Statistics in epidemiology: methods, techniques and applications. Boca Raton: CRC Press; 1995.Google Scholar
- United Nations Inter-agency Group for Child Mortality Estimation. Levels and trends in child mortality: 2012 report. 2012.
- United Nations Department of Economic and Social Affairs. World Population Prospects, the 2012 revision. 2015. http://esa.un.org/unpd/wpp/index.html. Accessed 29 May 2015.
- WHO. World Malaria Report. Geneva: World Health Organization; 2012. http://www.who.int/malaria/publications/world_malaria_report_2012/en/index.html. Accessed 29 October 2013.
- Griffin JT, Ferguson NM, Ghani AC. Estimates of the changing age-burden of Plasmodium falciparum malaria disease in sub-Saharan Africa. Nat Commun. 2014;5:3136.PubMedPubMed CentralView ArticleGoogle Scholar
- Hay SI, Okiro EA, Gething PW, Patil AP, Tatem AJ, Guerra CA, et al. Estimating the global clinical burden of Plasmodium falciparum malaria in 2007. PLoS Med. 2010;7:e1000290.PubMedPubMed CentralView ArticleGoogle Scholar
- Maire N, Tediosi F, Ross A, Smith T. Predictions of the epidemiologic impact of introducing a pre-erythrocytic vaccine into the expanded program on immunization in sub-Saharan Africa. Am J Trop Med Hyg. 2006;75:111–8.PubMedGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffin JT, Ghani AC. The potential public health impact and cost-effectiveness of the RTS,S vaccine in Africa. Abstract P271 [abstract]. Multilateral initiative on malaria (MIM) 6th Pan-African Conference. Durban, South Africa; 2013, p. 320.
- Ross A, Maire N, Molineaux L, Smith T. An epidemiologic model of severe morbidity and mortality caused by Plasmodium falciparum. Am J Trop Med Hyg. 2006;75:63–73.PubMedGoogle Scholar
- McCarthy KA, Wenger EA, Huynh GH, Eckhoff PA. Calibration of an intrahost malaria model and parameter ensemble evaluation of a pre-erythrocytic vaccine. Malar J. 2015;14:6.PubMedPubMed CentralView ArticleGoogle Scholar