-
Loading metrics
Open Access
Peer-reviewed
Research Article
Performance characteristics and potential public health impact of improved pre-erythrocytic malaria vaccines targeting childhood burden
Performance characteristics and potential public health impact of improved pre-erythrocytic malaria vaccines targeting childhood burden
- Josephine Malinga,
- Lydia Braunack-Mayer,
- Thiery Masserey,
- Aurélien Cavelan,
- Daniel Chandramohan,
- Alassane Dicko,
- Jean-Bosco Ouédraogo,
- Sherrie L. Kelly,
- Epke A. Le Rutte,
- Narimane Nekkab
- Published: August 4, 2025
- https://doi.org/10.1371/journal.pgph.0004549
Figures
Abstract
New malaria vaccine development builds on groundbreaking recommendations and roll-out of two approved pre-erythrocytic vaccines (PEVs); RTS,S/AS01 and R21/Matrix-M. Whilst these vaccines are effective in reducing childhood malaria within yearly routine immunization programs or seasonal vaccination, there is little evidence on how different PEV efficacies, durations of protection, and spacing between doses influence the potential to avert uncomplicated and severe childhood malaria. Mainly, lacking understanding of the required vaccine properties and delivery strategies that lead to an effective childhood vaccine with multi-year protection. We used an individual-based model of malaria transmission informed by trial data to quantify trade-offs between PEV performance properties and impact across different endemicities, deployment schedules, and coverage levels. We found that deploying a vaccine with 90% initial efficacy against infection, with a six to 12-month half-life duration of protection, co-administered with a blood-stage drug, followed by yearly boosters, results in 60–80% yearly incidence reduction, consistent with seasonal RTS,S and R21 trials. Halting vaccination after five years, leads to sustained protection of at least a 35% incidence reduction in children less than six years in the 12 months following cessation in settings where PfPR2–10 < 30%. Increasing the half-life duration to 12–18 months or reaching more children provides the same health impact with lower vaccine efficacy. Without a booster (fourth dose), high efficacy against infection (>90%) and longer half-life duration (>12 months) are required to sustain impact beyond primary vaccination, averting up to half the preceding year’s burden. The contribution of each property to the overall impact varies by setting and clinical endpoint, indicating that public health goals should dictate key vaccine performance criteria. Overall, our findings support the need for well-defined target product profiles for long duration vaccines linking priority use cases of where, how, and to whom to deploy new malaria vaccines, to maximize public health impact.
Citation: Malinga J, Braunack-Mayer L, Masserey T, Cavelan A, Chandramohan D, Dicko A, et al. (2025) Performance characteristics and potential public health impact of improved pre-erythrocytic malaria vaccines targeting childhood burden. PLOS Glob Public Health 5(8): e0004549. https://doi.org/10.1371/journal.pgph.0004549
Editor: Paulo Jorge Gonçalves de Bettencourt, Universidade Catolica Portuguesa, PORTUGAL
Received: January 8, 2025; Accepted: July 11, 2025; Published: August 4, 2025
Copyright: © 2025 Malinga et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The datasets, code, and plot scripts used and analyzed in our study are all open access and available via the following link: https://doi.org/10.5281/zenodo.14012674.
Funding: This study was funded by the Bill and Melinda Gates Foundation (INV-002562 to MAP). MAP and TM acknowledge support from the Swiss National Science Foundation (SNF Professorship PP00P3_170702 and PP00P3_203450 to MAP). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: I have read the journal’s policy and the authors of this manuscript have the following competing interests: AD is currently on the editorial board of this journal. All other authors declare no competing interests. This does not alter our adherence to PLOS Global Public Health policies on sharing data and materials.
Introduction
As of 2024, the World Health Organization (WHO) has recommended two pre-erythrocytic vaccines (PEVs), RTS,S/AS01 and R21/Matrix M, for global use against Plasmodium falciparum (Pf) malaria in pediatric populations [1]. It is anticipated that vaccine delivery will be aligned with existing routine childhood immunization platforms in perennial settings, or as seasonal mass vaccination before peak transmission among children in areas with moderate to high malaria transmission [1,2]. RTS,S and R21 both act at the pre-erythrocytic stage by targeting the sporozoite surface antigen of the Plasmodium falciparum parasite to prevent infection [3,4]. The groundbreaking WHO recommendation for RTS,S followed three extensive studies confirming the safety and efficacy of the vaccine. These include the RTS,S Phase 3 trials [5–7], an implementation program for vaccination of children aged five to nine months via the expanded program of immunization (EPI) with a fourth dose at 27 months [8], and an implementation study of seasonal use case of RTS,S deployed in combination with, or as an alternative to, seasonal malaria chemoprevention (SMC) with sulfadoxine-pyrimethamine and amodiaquine (SPAQ) [9,10]. Most recently, promising evidence from Phase 2 and 3 clinical trials for R21 for both seasonal and perennial use cases reports protective efficacy of more than 75% over 12 months (comparable to RTS,S protective efficacy six months after administration)[11], with further analysis of the follow-up results pending [12–14]. In the short-term, these are likely to be the only malaria vaccines in use, with the choice of implementation strategy predominantly depending on their supply and operational system factors.
Besides RTS,S and R21, other vaccine candidates are in pre-clinical or clinical trial stages, such as the whole sporozoite vaccine PfSPZ [15] and the blood-stage protein vaccine RH5 VLP [16]. There is also a renewed interest and investment in developing novel malaria vaccines, including mRNA vaccines [17]. Development of these new vaccine candidates comes at a crucial time when global progress in the malaria response has stalled [2,18]. Due to drug-resistant parasites, insecticide-resistant mosquitoes, funding needs, climate change and other factors, many African countries are off-track to meet the 2016 − 2030 Global Technical Strategy (GTS) targets to reduce global malaria incidence and mortality rates by at least 90% by 2030 over 2015 levels [19]. Consequently, the WHO and partners have called for revitalized efforts and the use of new tools to maintain the substantial gains witnessed in previous years [2]. Incorporating a vaccine into the existing and diverse malaria toolbox of interventions is a major milestone that could aid in achieving these targets and increase the proportion of children covered by any intervention [19]. Therefore, there is a need to optimize current vaccine implementation using existing delivery strategies and to understand the preferred vaccine properties, such as efficacy, duration of protection and dosage intervals of new and improved vaccines on their own, as well as alongside other novel interventions for malaria control [2]. Assessing how such vaccine properties are linked to public health benefits and understanding vaccine performance early in clinical development, including the vaccine's mode of action and immunogenicity, is essential to support new vaccines to achieve more significant impact. This will enable stakeholders to make informed investment decisions and streamline candidate selection in the Research and Development (R&D) phase of vaccine development.
In 2022, the WHO issued preferred product characteristics (PPCs) for malaria vaccines, providing updated advice on requirements for new vaccine candidates [20]. Informed by multiple stakeholders and public consultation, three strategic goals were identified, the first of which is to develop malaria vaccines that reduce morbidity and mortality in individuals at risk. While the document did not specify strict modes of action, it is outlined in the strategic goals that vaccines are envisaged to provide immunological protection against clinical and severe malaria targeting pre-erythrocytic or blood-stage antigens. Strict efficacy and duration requirements for burden reduction only vaccines were not explicitly defined, though preferred targets against clinical malaria over 12 months were identified. Of note, the PPC highlights that lower clinical efficacy thresholds can be justified in parallel with longer duration of protection, as well as other key drivers of public health impact, including vaccination coverage [20].
The PPCs document also identified the role of mathematical transmission modelling to support and guide discussions around vaccine impact and performance characteristics [20]. To date, mathematical modelling groups have provided a range of quantitative analyses to support thinking and policy decisions on malaria vaccines. Several studies have used results from RTS,S clinical trials to inform detailed models of malaria transmission and intervention dynamics, predicting the likely population-level health impact and cost-effectiveness of such vaccines. These studies examined vaccines deployed alone or in combination with other malaria interventions [7], as part of EPI for infants and children [5], mass vaccination [21], or seasonal use [7]. Other modelling studies have explored target efficacy profiles and decay properties of vaccines for mass vaccination with expanded age groups to support significant prevalence reduction [[22]], vaccines for transmission-blocking [23], or for improved childhood vaccines for EPI use [24]. However, there is still limited evidence from modelling studies of how improving pre-erythrocytic vaccine performance properties and optimizing vaccine deployment could increase public health impact. This includes understanding the impact of vaccine delivery, where boosters are given every other year rather than yearly, and how implementation factors drive vaccine impact.
In this study, we link the full range of vaccine properties, deployment schedules and vaccination coverage to different health outcomes using a detailed simulation model of malaria transmission and vaccines. Firstly, we identify the impact of improving the initial efficacy against infection and duration of vaccine protection on different clinical outcomes, such as averting uncomplicated and severe childhood malaria. Secondly, we investigate how vaccine impact differs by delivery strategy or how impact is driven by system factors such as coverage. More specifically, we focus on understanding the public health impact of implementing improved PEVs, with duration of protection longer than existing vaccines (such as RTS,S), delivered via routine immunization or mass vaccination campaigns followed by annual boosters for five years. We assess vaccine impact by predicting the reduction in infection prevalence, and incidence of clinical and severe cases achieved over the 12 months following the final annual booster in the fifth year. To explore the potential for multi-year vaccine impact, we evaluate the extended protection in the year following primary vaccination for children who do not receive any booster. Through simulation and sensitivity analysis, we provide a quantitative understanding of the trade-offs between vaccine performance and implementation impact-drivers across different endemicities and delivery schedules.
The novel insights from our modelling around the impact of improved PEVs can support the optimization of new malaria vaccine development. Moreover, our findings support improved understanding of current PEVs and their population impact, as well as our understanding of the potential public health benefits of deploying improved PEVs.
Materials and methods
Malaria transmission model
Model simulations were performed using a validated stochastic, individual-based model of malaria transmission in humans, linked to a deterministic model of malaria in mosquitoes, known as OpenMalaria [25–29]. The model is fully open access, has previously been described in [29], and the details regarding our current application with this model are summarized in S1 Text. OpenMalaria facilitates impact predictions for a wide range of interventions that target different stages in the parasite life cycle, including bed nets, chemoprevention, and vaccines. Intervention impact can be assessed for various health outcomes, including infection prevalence, uncomplicated and severe malaria cases and deaths. Malaria vaccine strategies have previously been simulated using this model, informed by and calibrated against estimates for RTS,S vaccine's protective efficacy from clinical trial data [5,30] and a range of other vaccine trial and implementation data across varying transmission settings [5,21,22]. As outlined below, we simulated a range of vaccine properties, namely, potential efficacies and durations of protection for improved PEVs for different endemicities, seasonal profiles and deployment strategies.
Simulated model scenarios, settings, and intervention dynamics
The scenarios modelled in this study include a range of vaccine properties of probable values for initial efficacy against infection and half-life duration for PEVs (Table 1). The range of settings represent different archetypal transmission profiles (short season, long season and constant transmission), prevalence levels and intervention coverage levels (reflecting access and uptake). The simulated vaccine deployment schedules encompass vaccine dosage either through EPI or through yearly mass vaccination, co-administration with or without curative malaria treatment and probabilities for accessing antimalarials. These scenarios were developed and refined based on estimates from modelling studies and stakeholder engagement. We further informed the scenarios given our insilico vaccine dose efficacy validation exercise, which used data from a clinical trial on seasonal vaccine use [9] (S1 Text). By simulating these scenarios with wide-ranging parameter values, we captured an extensive spectrum of epidemiological malaria dynamics.
Vaccine properties.
The initial vaccine efficacy against the pre-erythrocytic stage is assumed to reach its maximum level following the primary series. For RTS,S and R21, the primary series includes three vaccine doses, but for new vaccines the primary series may involve fewer doses. In our current study, the vaccine-induced efficacy against infection is assumed to be negligible before administering the final (third) dose in a primary series. Booster doses are assumed to restore the waning vaccine efficacy against infection to levels between 50% and 90% lower than the initial maximum efficacy against infection reached, based on reported values from previous studies [30] and our validation results (see next section and S1 Text). The decay in protective efficacy over time is assumed to follow a Weibull function with a biphasic shape parameter k = 0.69, with a rapid decay in the initial months, which is identical both after the initial dose as well as the booster doses as estimated previously [30]. Additionally, in a field trial in seasonal settings, the RTS,S protective efficacy against clinical malaria was shown to decline over a three-year trial period, with more rapid decline in the initial six months [9,11]. A two year extension of the same trial also showed sustained protection following additional annual booster doses given before peak transmission until children reached five years of age [10].
Vaccine booster efficacy against infection validation.
Using a Bayesian optimization approach (S1 Text), we used results from a recently completed Phase three 3 clinical trial for seasonal vaccination with RTS,S in Mali and Burkina Faso [9] to validate the properties of seasonal vaccine boosters. More explicitly, this meant better understanding how the vaccine efficacy against infection following the boosters contrasted against data from the pivotal RTS,S clinical trial, conducted in seven African countries [6]. This comparison informed model assumptions around efficacy against infection estimates of the first booster dose given 12 months after the three-dose primary series, rather than 18 months after, as implemented in the original RTS,S trial [6,30]. The initial vaccine induced efficacy against infection following the primary series, 91.1% [95% CI 74.5–99.7%] [30], was used as a model input and this efficacy against infection was separately reproduced in silico using the OpenMalaria model to match the trial results (S1 Text). The booster doses, given 12 and 24 months following the primary series, were assumed to partially restore waning vaccine efficacy, although it remains unclear by how much and for how long protection is extended. We used our Bayesian optimization approach, which utilizes a Gaussian process regression model as the objective function [31], to find the optimal parameter values for booster efficacies against infection for the different trial arms, that minimize the residual sum of squares between the observed data and modelled outputs.
SMC with SPAQ half-life validation.
Similarly, using the Bayesian optimization approach described and the results from the same Phase 3 clinical trial in Mali and Burkina Faso [9], we sought to calibrate the preventive half-life duration of seasonal chemoprevention using SPAQ. In the control arm, SPAQ was deployed to a cohort of children alone then in another arm in combination with seasonal vaccination. We incorporated past estimates of initial efficacy against infection following dosing with SPAQ as inputs for our model [32]. We assumed that SPAQ acts by first clearing all blood stage infections, followed by preventive action represented by a Weibull decay function with shape parameter k = 5.40 [32]. We also validated the shape parameter, k.
Vaccine deployment.
As described above, we defined primary vaccination as receiving the first series of doses (for example, up to three doses) and full vaccination as receiving the primary series of doses and annual boosters up to age five. In this study, all vaccinations are delivered through two approaches to allow comparison between deployment schedules. In the first approach, which we refer to as hybrid vaccination, we deploy the three-dose primary series as part of the age-based immunization schedule. The initial vaccine doses are given continuously during the intervention period to children aged five, seven and a half, and nine months. In the second approach, which we refer to as mass vaccination, children aged between five and 17 months receive the three-dose primary series through a mass campaign, timed so that the third dose is given one month before the transmission season's peak. For both approaches, additional annual boosters are deployed over four years to the same children up to five years of age, one month prior to peak malaria transmission in seasonal settings (Fig 1). In settings where the transmission is constant throughout the year, we follow the same delivery schedules for the primary series with the additional booster doses given through yearly mass campaigns (Fig A in S1 Text). The dropout rate is assumed to be 20% for all booster doses compared to the primary series. Vaccine doses are deployed singly or co-administered with a highly efficacious antimalarial treatment, modelled as blood-stage parasite clearance over five days to represent treatment with artemether-lumefantrine.
Illustration shows the five-year vaccination program with a primary series dosage, the timing of the annual boosters and at the bottom of the figure the age groups at each dose for both vaccine delivery schedules, hybrid vaccination (grey boxes) and mass vaccination (yellow boxes). For the mass vaccination schedule the primary series and annual booster doses are deployed before the peak transmission season. For the hybrid vaccination schedule doses for the primary series are deployed continuously as part of an age-based immunization schedule while the annual boosters are also deployed before the peak transmission season. The zero to 18 months period (green shaded areas) shows when the multi-seasonal or multi-year vaccine impact is evaluated by comparing a cohort of children who received the primary series doses and the first booster (dose four) and those who only received the primary series doses. The 12-month period (grey shaded area) follows the final annual booster dose (dose seven) in the fifth year where the public health vaccine impact is evaluated for children who received the primary series followed by annual boosters.
Endpoints to assess vaccine impact
Three public health outcomes were evaluated, including the relative reduction in infection prevalence, incidence of clinical cases, and incidence of severe cases, all compared to a no-intervention counterfactual. All three health outcomes were evaluated for two target age groups and follow-up periods. First, the vaccine impact was evaluated in children aged six years and below, 12 months following the final annual booster dose (dose seven) in the fifth year. We estimated vaccine impact across all the simulated scenarios for the two deployment schedules (Fig 1). Second, we evaluated the multi-year vaccine impact in the two years following the three-dose primary series by comparing a cohort of children who received only the primary series with those who received the primary series plus one booster dose (dose four). For the multi-year impact, the 24-month follow-up period was divided into six-month intervals. This multi-year impact is intended to assess a vaccine's extended protection in the first and second years if children do not receive booster doses, particularly in seasonal settings.
Statistical and global sensitivity analyses to evaluate vaccine impact
Using our individual-based stochastic malaria transmission model, we simulated experiments matching the scenarios described above. For the vaccine properties (initial efficacy against infection and half-life duration) and coverage, we generated a Latin hypercube of 1000 samples, and for each combination simulated outcomes for five replicates. We used heteroskedastic Gaussian process regression (hetgp package in R [33]) for each scenario to fit a model emulator to our database of model simulations. This emulator could capture the relationship between key vaccine performance properties and other factors, such as coverage or access to treatment, as regression inputs, and health outcomes, as outputs [34]. Exploring the entire parameter space for the different combinations of vaccine properties, deployment schedules, endemicities, and seasonality profiles requires a large number of simulations, which is computationally intensive. At low computational cost, emulators captured the relationships between vaccine properties and system factors, as well as the predicted health outcomes [22]. These emulators were then used to predict the vaccine's impact over the different follow-up periods and target age groups. We evaluated emulator performance by testing 10% of the total simulations against 90% used in the training set.
To identify the most important drivers of vaccine impact for different settings, health outcomes, and follow-up periods, we conducted a global sensitivity analysis of our Gaussian process regression model results using the Sobol method [35] and reported total effect indices. To calculate the relative contribution of each property, the total-order effect indices were normalised. The sensitivity analysis measured the extent to which a small change in an intervention's key performance property corresponded to a change in its impact. For example, an increase in the initial maximum vaccine efficacy against infection from 50% to 90% may lead to a larger change in the achievable clinical incidence reduction than a six-month increase in the vaccine's duration of protection. All analyses were conducted in R-software (version 4.1.0) [36].
Results
Validation of vaccine booster efficacy against infection using clinical trial data
Using the earlier described Bayesian optimization approach [31], we could determine booster dose efficacy against infection estimates from the clinical trial arm where only the RTS,S vaccine was deployed and in the arm where RTS,S was combined with SMC. By incorporating past estimates of the vaccine's initial efficacy against infection following the three-dose primary series of 91.1% and the half-life duration of 7.32 months as model inputs [30], we were able to approximate observed data from the trial to match our simulated modelling outcomes. Assuming a biphasic decay in vaccine efficacy and our Bayesian optimization approach (S1 Text), we estimated a best fit to incidence data for the vaccine's maximum efficacy against infection after the first booster (dose 4 given 12 months after dose 3) of 78.96% [95%CI 70.74 – 87.18%] in the vaccine only arm and 77.45% [95%CI 68.26 – 86.62%] when the vaccine was combined with SMC. We used trial data aggregated from both countries, Burkina Faso (high prevalence, seasonal transmission) and Mali (moderate prevalence, highly seasonal transmission). Following the second booster (dose 5, given 12 months after dose 4), the maximum boosting efficacy against infection was estimated to be 57.31% [95%CI 44.05 – 70.56%] in the vaccine only arm and 48.76% [95%CI 34.57 – 62.95%] in the combination arm. Fig 2 shows the monthly clinical incidence throughout the study period, reflecting the known estimates for initial efficacy against infection achieved following primary vaccination and the estimated assumptions for the efficacy against infection achieved following the two booster doses (doses four and five) administered prior to the peak malaria season in Burkina Faso and Mali that best fit the clinical trial data.
A) Best fit parameters for Burkina Faso (high prevalence and seasonal transmission), Mali (moderate prevalence and highly seasonal transmission) and both countries aggregated shown for the three trial arms, B) Hazard ratios between the trial arms for both countries separately and combined. The black dots shown with 95% confidence intervals represent the trial field data and the blue lines/triangles illustrate the modelled output from the simulations with the shaded region showing the confidence intervals averaged over 100 seeds. In this figure, the parameters were optimized for the arms where seasonal chemoprevention or vaccination were given alone and used to simulate the model results for the arm where the vaccine and chemoprevention were combined.
While there were slight variations in the vaccine boosting efficacy against infection estimates between the trial arms, there were overlaps within the error margins of uncertainty showing good agreement. The estimates ranged from 5% to 30% lower than the initial efficacy against infection reached after primary vaccination (Fig 2 and Fig C in S1 Text). Differences in vaccine impact between the two countries could also be observed as seen in the trials which could be attributed to reported variations in malaria epidemiology, as measured in children aged six to 12 years old at the end of the peak transmission seasons during the trials, but also likely due to differential profiles of exposure and acquisition of immunity in the two countries [9]. Thus, understanding of country- or endemicity-specific evaluation of vaccine performance in addition to global or archetypical estimates of vaccine efficacy and duration of protection is important. Hazard ratios calculated from both countries were also shown to fit the clinical trial data for both countries across all the trial years (Fig 2, Fig D in S1 Text).
Validation of SPAQ preventive half-life using clinical trial data
We replicated the optimization approach for the chemoprevention only arm where SPAQ was deployed as part of seasonal chemoprevention for four monthly cycles each year of the trial and in the arm where chemoprevention was combined with seasonal vaccination. We estimated the optimal range for the half-life duration of protection to be 24.6 days [95%CI 21.3 – 27.9 days] in the SMC only arm and 25.3 days [95%CI 22.9 – 27.7 days] in the arm where SMC was combined with RTS,S. Similarly, we used trial data aggregated for both countries, Burkina Faso and Mali. These estimates were lower than what has been reported previously from clinical studies [37,38] and through modelling [32]. We were also able to approximate the decay shape parameter, k, which matched our simulated model results in both countries of 4.05 [95%CI 3.33 – 4.76] in the SMC only arm and 3.82 [95%CI 2.92 – 4.71] when SMC was combined with seasonal vaccination. Fig 2 shows the monthly clinical incidence rate and hazard ratio estimates throughout the study period, reflecting the best fit assumptions for the half-life following administration of SPAQ administered prior to the peak malaria season in Burkina Faso and Mali. Results for the trial arm where RTS,S was deployed in combination with SMC are also shown in Fig 2, with additional results on best fit in Fig C in S1 Text.
Public health impact of improved PEVs on Plasmodium falciparum malaria burden
Modelling results indicate that implementing an improved PEV, targeting children aged over five months at recruitment, is expected to yield substantial impact on reducing infection prevalence, followed by decreases in clinical incidence, and then severe disease, when assessed among children under six years of age (Fig 3). As described above, vaccine impact is evaluated in the 12 months following the final annual booster dose in the fifth year. Across all modelled scenarios, our findings consistently confirm that co-administering a PEV with a blood-stage parasite clearance drug leads to a substantially greater reduction in disease burden compared to vaccination alone (Fig E in S1 Text). Subsequent sections present results for scenarios where each vaccine dose was co-administered with a blood-stage clearance drug.
A) Target reduction (%) in infection prevalence (top row), clinical incidence (middle row), and severe disease (bottom row) illustrating trade-offs between initial efficacy against infection and half-life duration of protection in settings where baseline PfPR2-10 ranged between 10% and 20%. The initial efficacy against infection ranged from 50% to 100%, half-life duration from six to 18 months, and assuming a primary series vaccination coverage of 90%. Burden reduction improved with increasing efficacy against infection and longer half-life duration of protection. B) Median (interquartile range (IQR)) relative reduction in infection prevalence (top row), clinical incidence (middle row), and severe disease (bottom row), considering varying levels of baseline PfPR2-10, coverage, and initial efficacy against infection for a long duration vaccine with a half-life duration between 12 and 18 months. Burden reduction was highest with high vaccination coverage and high initial efficacy against infection. Results are shown for a PEV co-administered with a blood stage clearance drug, for a four-month short seasonality profile, for both deployment schedules, in settings with a 30% probability of accessing curative treatment within 14 days of symptom onset.
The most notable relative reduction in disease burden occurred in areas with low and low-moderate transmission (PfPR2–10 < 30%) and decreased with increasing transmission (Fig 3B). Vaccine impact defined as the relative reduction in burden 12 months after the final annual booster, increased with improved vaccine performance and higher vaccination coverage (Fig 3A, 3B). Improving the initial efficacy against infection by, for instance, increasing it from 50% to 90% or by increasing vaccination coverage from 50% to 90%, resulted in almost a twofold increase in impact. Vaccines with extended half-life durations could offer protection for multiple years, while booster doses strengthen this protection, albeit contingent on dropout rates. From our trade-off analysis, our results show that we still need high initial efficacy against infection if the half-life duration is less than 12 months. Increasing this half-life to between 12 and 18 months or reaching more children can allow us to achieve the same impact with lower levels of vaccine efficacy (Fig 3, Fig H in S1 Text).
Across all contexts, burden reduction varied only slightly between the two delivery schedules (Fig 3). In highly seasonal settings, the hybrid vaccination schedule demonstrated a marginally better potential when compared to settings with longer seasons where the mass vaccination schedule showed a higher impact (Fig F in S1 Text). This preference stems from the fact that when vaccines are administered through mass campaigns, they simultaneously provide maximum protection to more children and this, in turn, leads to a greater impact on the transmission dynamics and, as such, to less infections. Conversely, administering the primary series continuously through a hybrid schedule could provide more protection, especially during short peak seasonal transmission. We observed similar findings in settings where transmission is assumed to be constant throughout the year. Burden reduction was slightly more pronounced when all vaccine doses were deployed through yearly mass campaigns, as opposed to hybrid vaccination (Fig F in S1 Text). We also found that, in all settings, relative burden reduction was substantially lower if health systems were stronger, particularly where higher levels of access to curative malaria treatment for a clinical case were available (Fig G in S1 Text).
Factors influencing vaccine impact following the five-year vaccination program
Reaching more children with a vaccine leads to greater impact. From our sensitivity analysis, the most important driver of impact is the proportion of children reached with the primary series doses, or vaccination coverage, followed by the vaccine's initial efficacy against infection and half-life duration. This is especially the case for infection-related endpoints and on a smaller scale, when vaccine administration is through hybrid vaccination compared to the mass vaccination schedule (Fig 4).
Our key model findings also underscore the crucial role played by both the initial efficacy against infection and the half-life duration of protection. We found that the relative contribution of each vaccine property to the overall impact depends on the clinical endpoint of interest, the seasonality patterns, the timing and length of the evaluation period, and whether the vaccine was co-administered with a blood-stage clearance drug (Fig 4 and Fig I, J in S1 Text). For instance, when endpoints are evaluated closer to the biological time of action, such as the 12 months following the final annual booster, burden reduction is primarily driven by the initial efficacy against infection, as shown in Fig 4. However, by enhancing the initial efficacy against infection to >90% and achieving high coverage of >80% with primary doses, the half-life duration becomes the predominant driver for burden reduction following the final annual booster, which in some instances accounts for over 60% of the total impact. Our findings regarding the significance of initial efficacy against infection and half-life duration hold consistent across different delivery schedules. However, our results suggest higher requirements for vaccine performance when existing infections are not pre-cleared during vaccination. The vaccine properties drive most of the impact in these cases, particularly for severe clinical outcomes (Fig I in S1 Text).
Multi-year vaccine impact in the two years following primary vaccination
To estimate the extended protection in the second year following primary vaccination, we evaluated impact in the 24-month period by following children who did not receive booster doses (Fig 5). Multi-year impact was evaluated by comparing children who received only the primary series doses in the first year to those who received the primary series doses plus one booster (fourth dose) in the second year (Fig 5, Fig M in S1 Text). Without the booster, a modeled PEV with a half-life between six and 12 months still provided extended protection beyond the first year following primary vaccination, covering part of the subsequent year’s second season (Fig 5). Notably, when the half-life duration of protection was longer than 12 months, a substantially higher burden reduction was predicted in the second year. This was due to the vaccine's protection extending to encompass the entirety of the second season, although at reduced efficacy. This phenomenon was particularly apparent in settings with pronounced seasonality, characterized by shorter periods of higher transmission, suggesting that the vaccine could be classed as multi-seasonal or multi-year (Fig 5). With yearly boosters given through mass vaccination in seasonal settings, our results show a 60–80% incidence reduction each year (Fig M in S1 Text), consistent with seasonal RTS,S and R21 trials in settings where PfPR2–10 < 30%. For both cohorts, our findings show a higher impact during the second year when the primary doses were deployed through the mass vaccination schedule before the first year’s peak season, in comparison to hybrid vaccination. This could be attributed to the mass deployment of primary series doses to more children simultaneously in the first year, preventing a higher proportion of infections even when booster doses were not given (Fig 5, Fig L in S1 Text).
A) Trade-offs between initial efficacy against infection and half-life duration of protection for clinical incidence reduction in settings where baseline PfPR2-10 ranged between 10% and 20%, initial efficacy against infection ranged from 50% to 100%, and half-life duration from six to 18 months, assuming a primary series vaccination coverage of 100%. B) Median (interquartile range (IQR)) relative reduction in clinical incidence for different levels of PfPR2-10 and initial efficacy against infection, for a long duration vaccine with a half-life between 12 and 18 months and primary series vaccination coverage of 100%. C) The bars represent the total Sobol effect indices which quantify the variance in predictions of relative reduction in clinical incidence. These indices can be interpreted as the proportion of variation in the outcome attributed to changes in each variable. Results are shown over different parameter ranges for initial efficacy against infection (75% to 100%) and half-life duration (six to 18 months), where vaccination coverage was fixed at 100%. Results are shown for a PEV co-administered with a blood stage clearance drug, for the zero- to six-, six- to 12-, 12- to 18- and 18- to 24-month periods following primary vaccination, for a four-month short seasonality profile, for the mass vaccination deployment schedule, in settings with a 30% likelihood of accessing curative treatment within 14 days of symptom onset.
Compared to the factors driving vaccine impact after five years of vaccination with boosters in previous sections, our analysis of multi-year impact confirms the vital role played by both the initial efficacy against infection and the half-life duration of protection (Fig 5c). Vaccinating all children in the cohort with the primary doses, we found that the importance of each vaccine property in driving impact varied by the delivery schedule and follow-up period. The impact of a multi-year PEV with an initial efficacy against infection >75% (comparable to existing vaccines like RTS,S) during the 12–24 months following primary vaccination is determined by the duration of antibody protection (Fig 5c and Fig L in S1 Text). However, while the initial efficacy against infection is a less influential driver for extended protection in the second year, a multi-year vaccine still needs to have a relatively high initial efficacy against infection to achieve adequate burden reduction (Fig 5). Such a vaccine would require a half-life duration longer than existing vaccines (at least 12 months) to achieve a burden reduction in the second season that is at least half of that estimated in the first season (Fig 5A, 5B and Fig L in S1 Text). For instance, in both the deployment schedules and without boosters, to achieve a > 30% reduction in clinical incidence in the 12- to 18-month period following a > 60% reduction during the zero- to six-month period after primary vaccination, a half-life duration of at least 12 months and an initial vaccine efficacy against infection of >70% are both required, if the vaccination coverage is 100%. The initial efficacy against infection drove most of the impact in the six-month period following the mass vaccination schedule (Fig 5C). However, this was not the case for the hybrid vaccination approach, where the main driver of impact was the half-life duration across all the follow-up periods (Fig L in S1 Text).
Discussion
Our modelling results show that PEVs with improved duration of protection and vaccine-induced protective efficacy have the potential for increased benefit in reducing childhood malaria. Moving beyond current CSP-targeting vaccines (i.e., RTS,S), we provide additional evidence on the relationships and trade-offs between vaccine performance and drivers of impact for these improved vaccines. Longer duration PEVs could provide extended protection in the years following primary vaccination, and while vaccine efficacy wanes over time, boosting can restore protection. Our results suggest that deploying a PEV with a half-life duration that could be extended by three to five months compared to that of RTS,S may result in sustained impact into the second and third year following primary vaccination. If the initial efficacy against infection for a new vaccine candidate is more than 90%, improving the half-life duration to more than 12 months can achieve a burden reduction of more than half relative to the preceding year’s level, even when a booster is not given in the second year. Moreover, if the duration of protection can be increased further, trade-offs could be made with lower vaccine efficacy for a similar impact. This means there is a crucial need to reliably measure the duration of protection of vaccines. While it is currently challenging to adequately measure duration before large scale clinical trials, early evidence could be measured in a controlled human malaria infection (CHMI) trial alongside reliable immune correlates of protection [39].
Recent evidence has shown promising results for the RTS,S vaccine delivered through EPI or seasonal mass campaigns [8,9,10], with additional encouraging results reported for R21 [13,14]. While these might be the only vaccines available for a few years, our findings indicate that use cases for improved PEVs should adapt and capitalize on the benefits of multi-year protection. If a new vaccine can provide multi-year protection, a program can reduce the number of booster doses to be delivered to at least every two years instead of annually. In addition, when implementing multi-year PEVs with a half-life duration of protection greater than 12 months, children missed during annual seasonal vaccination will remain partially protected. If more children are reached with the vaccine, individuals could also benefit from the indirect effects of reduced malaria transmission across the entire population. Developing such a vaccine and combining it with other malaria interventions, such as vector control or chemoprevention, can potentially reduce severe disease burden, save lives, and accelerate elimination efforts.
We parameterized boosting efficacies against infection of RTS,S and SMC properties in our model using trial results of a seasonal vaccination and chemoprevention trial [9]. Accounting for trial site characteristics we were able to match our modelling results to trial data from both Mali and Burkina Faso. We reproduced trial findings that showed increased impact following seasonal boosting, similar to our modelling results, when the booster dose was given 12 and 24 months following primary dosing. We also provided new estimates for the half-life duration and the protective decay of SPAQ adding onto the evidence base. These findings emphasize the need to understand how vaccines work in settings with varying transmission dynamics or how vaccines interact with other interventions. In particular, we recognize the importance of clearly defining temporal transmission profiles for seasonal vaccines or chemoprevention to ensure deployment before or during the maximum risk period in seasonal settings for increased impact. We carefully explored the parameterization of the RTS,S vaccine given alone or in combination with SMC in these seasonal settings, thus also validating intervention properties. Similar to clinical trial findings [9], we found higher impact when combining the PEV with antimalarials. From our modelled scenarios, co-administering a PEV with a blood-stage clearance drug to clear existing infections leads to higher impact than when the vaccine is deployed alone. It has previously been shown that combining vaccines with mass drug administration in the final stages of an elimination program, may enhance the success of interrupting transmission as opposed to single deployments [21]. Where the goal has been predominantly to reduce burden, higher impact was seen when RTS,S was deployed in combination with antimalarials, as shown in a Phase 3 trial for seasonal vaccination plus SMC using SPAQ [9,10] and in a Phase 2 trial combining RTS,S with dihydroartemisinin, piperaquine, or primaquine [40]. Interactions with combination drugs should be investigated and resolved to understand their relative and joined efficacies, as well as to cultivate community uptake of such treatment strategies.
Overall, we show that the impact of PEVs on reducing burden is consistent across transmission settings and higher for infection endpoints than clinical and severe outcomes. This indicates that, considering alternative clinical endpoints may allow better evaluation of vaccine efficacy in Phase 2 and 3 clinical trials conducted in diverse settings. For PEVs, an infection endpoint is closer to the mode-of-action of the vaccine and it can be translated across settings and, furthermore, offers a systematic and unbiased way to compare different malaria vaccines. At present, efficacy endpoints for current vaccines are generally reported against uncomplicated and severe disease. Endpoints from existing trials are difficult to interpret and compare, as they are measured for varying transmission settings, follow-up periods, and underlying age patterns of disease [6,14,15]. Future clinical trials could evaluate all infection endpoints with sensitive diagnostics or appropriate serological monitoring, particularly if vaccines are to be considered as part of the toolkit for achieving malaria elimination. Tracking the incidence of malaria infections in a smaller group during a clinical trial will help confirm the underlying efficacy of a vaccine against infection. Additionally, identifying and validating immune correlates including CHMIs, can help inform the early phases of clinical trials. For instance, this could support the evaluation of dosing regimens prior to conducting larger scale, later-stage clinical trials. This approach will be essential for reducing the time required to evaluate longer duration vaccines.
Our modelling results show that, for all health outcomes examined, impact increased with improvements in vaccine performance and the proportion of children vaccinated. Notably, our findings show more pronounced changes to vaccine impact for duration of protection improvements. While it may take several years to develop and approve a next-generation malaria vaccine (such as multi-stage vaccines or non PEV), improving the efficacy against infection, duration, or the delivery of PEVs could improve their potential health impact. Our findings also show that, modelling new candidates with optimized deployment strategies, such as timed seasonal deployment to increase coverage, will aid in the assessment of the required performance characteristics for their suitability as novel vaccines. Additional considerations for a prospective vaccine candidate should be evaluated in parallel with determining their efficacy profile, duration of protection, and the scheduling of booster doses. For instance, properties which could influence vaccine coverage include its formulation and the number of doses required, as they could impact acceptability and adherence, cold storage requirements and a potential low frequency of adverse effects [41]. Based on the results from our trade-offs analysis, we find that a vaccine candidate with a low initial efficacy against infection of 50%, less than for RTS,S, would likely have substantially lower impact, an outcome which is contingent on vaccine durability as well as factors such as health system capacity in the affected community. However, an implementation strategy could prioritize coverage by having fewer vaccine doses per child, at the cost of reducing vaccine efficacy as less doses are required during the primary series. Such a strategy may still translate to an equivalent public health benefit to that of a lower coverage, higher efficacy vaccine. This implies that continued efforts to develop vaccine candidates that are less efficacious than RTS,S remain worthwhile, even in the context of the current focus for vaccine development strategies to maintain protectiveness over time.
From our results, reaching children with a vaccine is the most important determinant of vaccine impact, particularly for low efficacy or shorter duration PEVs. This underscores the importance of evaluating and understanding both vaccine properties and operational factors that influence intervention access, delivery, and uptake. Particularly, in addition to reducing primary series doses, if vaccines are delivered through routine health systems. Improving coverage, for instance, by matching current booster dose timing to routine immunization schedules, fractional dosing, reducing supply and demand gaps, or sub-national targeted vaccination, could improve the expected impact [24,42]. Several challenges, including health system constraints, gaps in communication and engagement between stakeholders, and inadequate training and community sensitization, were identified in the piloting of RTS,S in Ghana, Kenya, and Malawi [43–45]. These are postulated to have contributed to inefficiencies in delivery, impeding high vaccine uptake [19,44,45]. The timing of the RTS,S fourth dose was also not aligned with existing childhood immunization schedules (the booster was given at 27 months of age whilst the more common measles and meningococcal vaccines are given at 18 months of age). Additionally, the age eligibility criteria did not match the burden in some countries [19,44,45]. We did not explicitly explore the influence of the efficacy of all doses in the primary series, or of reducing the number of primary doses in our study. However, a competitive edge could be provided by a PEV with fewer doses at a lower cost, with an advantage in supply, or with safety and efficacy demonstrated in school-going children and adults.
Our study has several limitations to consider when interpreting the results. First, we based our model parameterizations for an improved PEV on data and estimates from two clinical trials: a Phase 3 trial of RTS,S that included the likely protective efficacy decay profile after a three-dose primary series [30], and a trial of seasonal vaccination with RTS,S in two countries [9] to validate model assumptions of boosting efficacy. Estimates surrounding boosting efficacy against infection have yet to be exhaustively validated in the field and are hence uncertain. While our validation exercise captured some of these uncertainties, future modelling studies must describe these ambiguities fully. In particular, it will be critical to onboard new clinical evidence or knowledge around a new vaccine candidate's decay in protection. Better informed efficacy decay parameter estimates will directly improve model predictions of the potential PEV impact. Second, our findings are influenced by our assumptions on vaccine performance properties, including the selected parameter ranges for initial efficacy against infection and half-life duration of protection. We also evaluated limited values for boosting efficacy against infection and we did not examine vaccines with a half-life greater than 18 months. Third, our modelling scenarios are composed of archetypal seasonal transmission profiles and health system characteristics that broadly indicate the range of results for a particular prevalence setting. Moreover, we did not account for the heterogeneities in transmission or care-seeking likely to occur in endemic malaria settings [46]. Previous comparisons between model estimates for geographic specific locations and setting archetypes do, however, show that estimates are similar for childhood vaccination with limited indirect benefits [45]. Fourth, while our results suggest that combining PEVs with blood-stage clearance drugs yield higher impact, our modelling approach does not account for the emergence of drug resistance. Future modelling studies should consider the synergy and acceptability of such combinations. Lastly, vaccination coverage does not account for the nuances of access and uptake and only represents a broad metric of real-life implementation.
Conclusions
The development of a highly efficacious, durable vaccine remains a priority for the malaria vaccine research and development community. Thus, obtaining an early understanding of a vaccine's duration of protection and its efficacy decay profile is crucial. It is critical to incorporate an understanding of the duration of protection with appropriate clinical trial endpoints for burden reduction or infection prevention, alongside correlates of immunity. This will allow developers and stakeholders to assess and prioritize use cases aiming for greater public health benefits. Our modelling suggests that PEVs with high initial efficacy against infection of more than 90% and a half-life duration of protective efficacy greater than 12 months offer opportunities for protection over multiple years, suggesting that such a vaccine could reduce the need for yearly booster doses, particularly in lower transmission settings. Our modelling results also provide a better understanding of trade-offs between vaccine performance properties, health system, and programmatic factors, and could support decision making for both clinical investment in and recommendations for new or next-generation malaria vaccines.
Supporting information
S1 Text. Supporting Information for ‘Performance characteristics and potential public health impact of improved pre-erythrocytic malaria vaccines targeting childhood malaria burden’.
https://doi.org/10.1371/journal.pgph.0004549.s001
(DOCX)
Acknowledgments
We would like to thank our colleagues in the Global Disease Modelling (GDM) team at The Kids Research Institute Australia, Nedlands, Western Australia, and the Disease Modelling research unit at the Swiss Tropical and Public Health Institute, Allschwil and University of Basel, Basel, Switzerland for the support and discussions. We acknowledge and appreciate all the participants and investigators of the seasonal vaccination with RTS,S/AS01 and seasonal malaria chemoprevention with SP + AQ clinical trial in Mali and Burkina Faso. We would also like to acknowledge all participants of the 2021 ‘Malaria Prevention: Shaping Next-Gen Medical Interventions’ convening for their discussions and feedback. We also thank Jean-Luc Bodmer (Gates Foundation), Scott Miller (Gates Medical Research Institute), and Mary Hamel and Lindsay Wu (World Health Organization) for their support and discussions. Calculations were performed at sciCORE (http://scicore.unibas.ch/) scientific computing center at the University of Basel.
References
- 1. WHO. WHO recommends groundbreaking malaria vaccine for children at risk. 2021 [cited 14 Apr 2022]. Available from: https://www.who.int/news/item/06-10-2021-who-recommends-groundbreaking-malaria-vaccine-for-children-at-risk
- 2. WHO. World malaria report 2021. Geneva: World Health Organization; 2022. Available from: https://apps.who.int/iris/handle/10665/350147
- 3. Gordon DM, McGovern TW, Krzych U, Cohen JC, Schneider I, LaChance R, et al. Safety, immunogenicity, and efficacy of a recombinantly produced Plasmodium falciparum circumsporozoite protein-hepatitis B surface antigen subunit vaccine. J Infect Dis. 1995;171(6):1576–85. pmid:7769295
- 4. Collins KA, Snaith R, Cottingham MG, Gilbert SC, Hill AVS. Enhancing protective immunity to malaria with a highly immunogenic virus-like particle vaccine. Sci Rep. 2017;7:46621. pmid:28422178
- 5. Penny MA, Verity R, Bever CA, Sauboin C, Galactionova K, Flasche S, et al. Public health impact and cost-effectiveness of the RTS,S/AS01 malaria vaccine: a systematic comparison of predictions from four mathematical models. Lancet. 2016;387(10016):367–75. pmid:26549466
- 6. RTS,S Clinical Trials Partnership. Efficacy and safety of RTS,S/AS01 malaria vaccine with or without a booster dose in infants and children in Africa: final results of a phase 3, individually randomised, controlled trial. Lancet. 2015;386(9988):31–45. pmid:25913272
- 7. Thompson HA, Hogan AB, Walker PGT, Winskill P, Zongo I, Sagara I, et al. Seasonal use case for the RTS,S/AS01 malaria vaccine: a mathematical modelling study. Lancet Glob Health. 2022;10(12):e1782–92. pmid:36400084
- 8. WHO. Malaria: The malaria vaccine implementation programme (MVIP). 2021 [cited 18 May 2021]. Available from: https://www.who.int/news-room/q-a-detail/malaria-vaccine-implementation-programme
- 9. Chandramohan D, Zongo I, Sagara I, Cairns M, Yerbanga R-S, Diarra M, et al. Seasonal Malaria Vaccination with or without Seasonal Malaria Chemoprevention. N Engl J Med. 2021;385(11):1005–17. pmid:34432975
- 10. Dicko A, Ouedraogo J-B, Zongo I, Sagara I, Cairns M, Yerbanga RS, et al. Seasonal vaccination with RTS,S/AS01E vaccine with or without seasonal malaria chemoprevention in children up to the age of 5 years in Burkina Faso and Mali: a double-blind, randomised, controlled, phase 3 trial. Lancet Infect Dis. 2024;24(1):75–86. pmid:37625434
- 11. Cairns M, Barry A, Zongo I, Sagara I, Yerbanga SR, Diarra M, et al. The duration of protection against clinical malaria provided by the combination of seasonal RTS,S/AS01E vaccination and seasonal malaria chemoprevention versus either intervention given alone. BMC Med. 2022;20(1):352. pmid:36203149
- 12. Datoo MS, Natama MH, Somé A, Traoré O, Rouamba T, Bellamy D, et al. Efficacy of a low-dose candidate malaria vaccine, R21 in adjuvant Matrix-M, with seasonal administration to children in Burkina Faso: a randomised controlled trial. Lancet. 2021;397(10287):1809–18. pmid:33964223
- 13. Datoo MS, Natama HM, Somé A, Bellamy D, Traoré O, Rouamba T, et al. Efficacy and immunogenicity of R21/Matrix-M vaccine against clinical malaria after 2 years’ follow-up in children in Burkina Faso: a phase 1/2b randomised controlled trial. Lancet Infect Dis. 2022;22(12):1728–36. pmid:36087586
- 14. Datoo MS, Dicko A, Tinto H, Ouédraogo J-B, Hamaluba M, Olotu A, et al. Safety and efficacy of malaria vaccine candidate R21/Matrix-M in African children: a multicentre, double-blind, randomised, phase 3 trial. Lancet. 2024;403(10426):533–44. pmid:38310910
- 15. Sissoko MS, Healy SA, Katile A, Zaidi I, Hu Z, Kamate B, et al. Safety and efficacy of a three-dose regimen of Plasmodium falciparum sporozoite vaccine in adults during an intense malaria transmission season in Mali: a randomised, controlled phase 1 trial. Lancet Infect Dis. 2022;22(3):377–89. pmid:34801112
- 16. Payne RO, Silk SE, Elias SC, Miura K, Diouf A, Galaway F, et al. Human vaccination against RH5 induces neutralizing antimalarial antibodies that inhibit RH5 invasion complex interactions. JCI Insight. 2017;2(21):e96381. pmid:29093263
- 17. BioNTech. BioNTech Initiates Phase 1 Clinical Trial for Malaria Vaccine Program BNT165 | BioNTech. 2022 [cited 26 Jan 2023]. Available from: https://investors.biontech.de/news-releases/news-release-details/biontech-initiates-phase-1-clinical-trial-malaria-vaccine/
- 18. Noor AM, Alonso PL. The message on malaria is clear: progress has stalled. Lancet. 2022;399(10337):1777. pmid:35461617
- 19. WHO. World malaria report 2024. 2024 [cited 2 Apr 2025]. Available from: https://www.who.int/teams/global-malaria-programme/reports/world-malaria-report-2024
- 20. WHO. Malaria vaccines: preferred product characteristics and clinical development considerations. 2022 [cited 26 Jan 2023]. Available from: https://www.who.int/publications-detail-redirect/9789240057463
- 21. Camponovo F, Ockenhouse CF, Lee C, Penny MA. Mass campaigns combining antimalarial drugs and anti-infective vaccines as seasonal interventions for malaria control, elimination and prevention of resurgence: a modelling study. BMC Infect Dis. 2019;19(1):920. pmid:31664924
- 22. Golumbeanu M, Yang G-J, Camponovo F, Stuckey EM, Hamon N, Mondy M, et al. Leveraging mathematical models of disease dynamics and machine learning to improve development of novel malaria interventions. Infect Dis Poverty. 2022;11(1):61. pmid:35659301
- 23. Sherrard-Smith E, Sala KA, Betancourt M, Upton LM, Angrisano F, Morin MJ, et al. Synergy in anti-malarial pre-erythrocytic and transmission-blocking antibodies is achieved by reducing parasite density. Elife. 2018;7:e35213. pmid:29914622
- 24. Hogan AB, Winskill P, Ghani AC. Estimated impact of RTS,S/AS01 malaria vaccine allocation strategies in sub-Saharan Africa: A modelling study. PLoS Med. 2020;17(11):e1003377. pmid:33253211
- 25. Smith T, Ross A, Maire N, Chitnis N, Studer A, Hardy D, et al. Ensemble modeling of the likely public health impact of a pre-erythrocytic malaria vaccine. PLoS Med. 2012;9(1):e1001157. pmid:22272189
- 26. Smith T, Killeen GF, Maire N, Ross A, Molineaux L, Tediosi F, et al. Mathematical modeling of the impact of malaria vaccines on the clinical epidemiology and natural history of Plasmodium falciparum malaria: Overview. Am J Trop Med Hyg. 2006;75(2 Suppl):1–10. pmid:16931810
- 27. Smith T, Maire N, Ross A, Penny M, Chitnis N, Schapira A, et al. Towards a comprehensive simulation model of malaria epidemiology and control. Parasitology. 2008;135(13):1507–16. pmid:18694530
- 28. Chitnis N, Hardy D, Smith T. A periodically-forced mathematical model for the seasonal dynamics of malaria in mosquitoes. Bull Math Biol. 2012;74(5):1098–124. pmid:22218880
- 29. Reiker T, Golumbeanu M, Shattock A, Burgert L, Smith TA, Filippi S, et al. Emulator-based Bayesian optimization for efficient multi-objective calibration of an individual-based model of malaria. Nat Commun. 2021;12(1):7212. pmid:34893600
- 30. Penny MA, Pemberton-Ross P, Smith TA. The time-course of protection of the RTS,S vaccine against malaria infections and clinical disease. Malar J. 2015;14:437. pmid:26537608
- 31. Balandat M, Karrer B, Jiang D, Daulton S, Letham B, Wilson AG, et al. BoTorch: A Framework for Efficient Monte-Carlo Bayesian Optimization. Advances in Neural Information Processing Systems. Curran Associates, Inc.; 2020. p. 21524–21538. Available from: https://proceedings.neurips.cc/paper/2020/hash/f5b1b89d98b7286673128a5fb112cb9a-Abstract.html
- 32. Burgert L, Reiker T, Golumbeanu M, Möhrle JJ, Penny MA. Model-informed target product profiles of long-acting-injectables for use as seasonal malaria prevention. PLOS Glob Public Health. 2022;2(3):e0000211. pmid:36962305
- 33. Binois M, Gramacy RB. hetGP: Heteroskedastic Gaussian Process Modeling and Sequential Design in R. J Stat Soft. 2021;98(13):1–44.
- 34. Rasmussen CE. Gaussian Processes in Machine Learning. In: Bousquet O, von Luxburg U, Rätsch G, editors. Advanced Lectures on Machine Learning: ML Summer Schools 2003, Canberra, Australia, February 2 - 14, 2003, Tübingen, Germany, August 4 - 16, 2003, Revised Lectures. Berlin, Heidelberg: Springer; 2004. p. 63–71.
- 35. Sobol′ IM. Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates. Math Comput Simul. 2001;55(1–3):271–80.
- 36. R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2021. Available from: https://www.R-project.org
- 37. Zongo I, Milligan P, Compaore YD, Some AF, Greenwood B, Tarning J, et al. Randomized Noninferiority Trial of Dihydroartemisinin-Piperaquine Compared with Sulfadoxine-Pyrimethamine plus Amodiaquine for Seasonal Malaria Chemoprevention in Burkina Faso. Antimicrob Agents Chemother. 2015;59(8):4387–96. pmid:25918149
- 38. Mousa A, Cuomo-Dannenburg G, Thompson HA, Bell DJ, D’Alessandro U, Gosling R, et al. Estimating the Impact of DHPS Mutations on Sulfadoxine-Pyrimethamine Protective Efficacy: A Pooled Analysis of Individual Patient Data and Implications for Malaria Chemoprevention in Sub-Saharan Africa. Rochester, NY; 2024.
- 39. Thompson HA, Hogan AB, Walker PGT, White MT, Cunnington AJ, Ockenhouse CF, et al. Modelling the roles of antibody titre and avidity in protection from Plasmodium falciparum malaria infection following RTS,S/AS01 vaccination. Vaccine. 2020;38(47):7498–507. pmid:33041104
- 40. von Seidlein L, Hanboonkunupakarn B, Jittamala P, Pongsuwan P, Chotivanich K, Tarning J, et al. Combining antimalarial drugs and vaccine for malaria elimination campaigns: a randomized safety and immunogenicity trial of RTS,S/AS01 administered with dihydroartemisinin, piperaquine, and primaquine in healthy Thai adult volunteers. Hum Vaccin Immunother. 2020;16(1):33–41. pmid:31306084
- 41. Praet N, Asante KP, Bozonnat M-C, Akité EJ, Ansah PO, Baril L, et al. Assessing the safety, impact and effectiveness of RTS,S/AS01E malaria vaccine following its introduction in three sub-Saharan African countries: methodological approaches and study set-up. Malar J. 2022;21(1):132. pmid:35468801
- 42. Regules JA, Cicatelli SB, Bennett JW, Paolino KM, Twomey PS, Moon JE, et al. Fractional Third and Fourth Dose of RTS,S/AS01 Malaria Candidate Vaccine: A Phase 2a Controlled Human Malaria Parasite Infection and Immunogenicity Study. J Infect Dis. 2016;214(5):762–71. pmid:27296848
- 43. Mwangoka G, Ogutu B, Msambichaka B, Mzee T, Salim N, Kafuruki S, et al. Experience and challenges from clinical trials with malaria vaccines in Africa. Malar J. 2013;12:86. pmid:23496910
- 44. Guignard A, Praet N, Jusot V, Bakker M, Baril L. Introducing new vaccines in low- and middle-income countries: challenges and approaches. Expert Rev Vaccines. 2019;18(2):119–31. pmid:30689472
- 45. Grant J, Gyan T, Agbokey F, Webster J, Greenwood B, Asante KP. Challenges and lessons learned during the planning and early implementation of the RTS,S/AS01E malaria vaccine in three regions of Ghana: a qualitative study. Malar J. 2022;21(1):147. pmid:35550113
- 46. Galactionova K, Tediosi F, de Savigny D, Smith T, Tanner M. Effective coverage and systems effectiveness for malaria case management in sub-Saharan African countries. PLoS One. 2015;10(5):e0127818. pmid:26000856