Partial Malaria Vaccination Requires Vector Control for Elimination, Mathematical Model Shows

A mathematical modeling investigation published in Discover Public Health demonstrates that deploying partial malaria vaccination reduces overall disease burden and delays infection peaks but remains fundamentally insufficient to achieve malaria elimination unless paired with vector-control measures. Conducted by researchers Kunj Chotaliya, Anil Chavada, Nimisha Pathak, and Akanksha Singh, the study formulates an eight-compartment deterministic system of ordinary differential equations designed to capture the nuanced realities of imperfect vaccine protection. While traditional epidemiological frameworks often treat immunization as an all-or-nothing barrier or introduce vaccination exclusively among susceptible hosts, this framework explicitly tracks both vaccinated and unvaccinated cohorts across exposed, infectious, and recovered stages. By incorporating partial efficacy, breakthrough infections, vaccine-assisted recovery, and the waning of both vaccine-induced and naturally acquired immunity, the researchers evaluated how demographic, immunological, and entomological factors interact. Their quantitative findings show that although ramping up immunization suppresses transmission potential, high mosquito-to-human transmission rates can overpower vaccine coverage, establishing that sustainable elimination demands integrated suppression of the mosquito vector.

Malaria persists as a devastating global public health challenge, predominantly across tropical and subtropical regions where warm, humid climates sustain mosquito breeding and parasite incubation. The vector-borne disease is caused by Plasmodium parasites—most notably Plasmodium falciparum, Plasmodium vivax, Plasmodium malariae, and Plasmodium ovale—and is transmitted to humans through the bites of infected female Anopheles mosquitoes. According to the World Malaria Report 2025 cited by the authors, the World Health Organization African Region bears an overwhelming proportion of global malaria cases and fatalities, with children under five years old facing the highest vulnerability. Over recent decades, public health programmes have relied heavily on classical control methods, including insecticide-treated bed nets, indoor residual spraying, rapid diagnostics, and artemisinin-based combination therapies. However, the increasing spread of insecticide-resistant mosquito populations has introduced substantial operational hurdles. While the introduction of malaria vaccines such as RTS,S/AS01 has provided a critical preventive tool, clinical evidence indicates that vaccine-induced protection is partial and wanes significantly over time, allowing vaccinated individuals to contract breakthrough infections and remain capable of secondary transmission.

Motivated by the incomplete protection observed in clinical settings, the authors expanded upon earlier deterministic models to construct a comprehensive human population model partitioned into eight mutually exclusive epidemiological compartments. The total human population comprises unvaccinated susceptible individuals, vaccinated susceptible individuals, unvaccinated exposed individuals, vaccinated exposed individuals, unvaccinated infectious individuals, vaccinated infectious individuals, unvaccinated recovered individuals, and vaccinated recovered individuals. Unvaccinated susceptible individuals enter through demographic recruitment or the loss of natural immunity, exiting via infection, natural mortality, or vaccination. Vaccinated susceptible individuals receive partial protection but remain vulnerable to breakthrough infections; they leave through breakthrough infection, the waning of vaccine-induced immunity back into the unvaccinated susceptible pool, or natural death. Exposed compartments represent individuals undergoing the latent incubation phase before becoming infectious. In the infectious classes, vaccinated patients are modeled with an enhanced recovery rate reflecting the clinical benefits of prior immunization, while both infected groups share a baseline disease-induced mortality parameter to preserve mathematical tractability.

To confirm that the mathematical structure represents biologically meaningful dynamics, the authors performed rigorous qualitative assessments of the differential equation system. They proved that all state variables remain non-negative for all positive time intervals whenever initialized with non-negative conditions, demonstrating that the non-negative orthant is positively invariant. The researchers further established the existence of a bounded, positively invariant feasible region defined by the demographic equilibrium ratio of recruitment to natural mortality. By reformulating the system into an equivalent Volterra integral equation, the team showed that the vector field satisfies a global Lipschitz continuity condition within the invariant space. Applying the Banach Fixed-Point Theorem, they demonstrated that the system guarantees existence and uniqueness of solutions over local intervals, which naturally extend across all continuous time horizons without finite-time blow-up. These mathematical proofs confirm that the model behaves stably and realistically, providing a sound theoretical foundation for analyzing long-term epidemiological outcomes and intervention strategies.

A central milestone in the analysis is the derivation of the basic reproduction number, denoted as R0, which quantifies the expected number of secondary cases generated by a single infectious human introduced into a fully disease-free community. Utilizing the next-generation matrix operator framework, the researchers decomposed the transmission dynamics into matrices describing newly generated infections and compartmental transition flows. The resulting reproduction number resolves into the sum of two distinct transmission pathways: one driven by unvaccinated infectious individuals and the other driven by vaccinated individuals suffering breakthrough infections, weighted by their relative infectiousness parameter. Under the theoretical condition where vaccination coverage is zero, the formula collapses cleanly into the standard classical reproduction number for unvaccinated malaria transmission. When evaluated with baseline parameters taken from published literature—including a baseline recruitment rate of 1.2 individuals per day, a natural mortality rate of 0.03 per day, and an initial transmission contact rate of 0.001—the baseline reproduction number was calculated to be approximately 0.1737, placing the baseline system safely below the epidemic threshold.

The researchers subsequently established the equilibrium points and their stability criteria, demonstrating how the basic reproduction number governs the ultimate fate of malaria within a host community. When the basic reproduction number is strictly less than one, the disease-free equilibrium is proven to be both locally asymptotically stable and globally asymptotically stable throughout the feasible domain. Under this condition, the number of exposed and infectious hosts inevitably decays toward zero over time, confirming that the infection is eliminated. Conversely, when the reproduction number exceeds unity, the disease-free state becomes unstable, and the system admits a unique, positive endemic equilibrium where malaria persists indefinitely. The authors verified the global asymptotic stability of this endemic state using a nonlinear Volterra-type Lyapunov function combined with LaSalle’s invariance principle. Under verified Lyapunov coefficient constraints, all trajectories originating within the interior of the feasible region converge inexorably toward the positive endemic steady state, demonstrating that disease persistence is structurally robust whenever transmission exceeds the threshold.

To explore the precise nature of the transition between disease extinction and disease persistence, the team conducted a center-manifold bifurcation analysis around the critical threshold where R0 equals one. Treating the effective mosquito-to-human transmission rate as the bifurcation parameter, the authors calculated the critical transmission threshold at approximately 0.0058 under baseline conditions. By determining the second-order partial derivatives of the system and evaluating them along the right and left eigenvectors of the Jacobian matrix, they calculated the center-manifold coefficients, demonstrating that the primary coefficient is strictly negative while the secondary coefficient is strictly positive. According to bifurcation theory, this configuration rules out subcritical backward bifurcations and confirms the presence of a forward transcritical bifurcation. From an epidemiological standpoint, this finding indicates that the endemic equilibrium branches continuously and stably away from the disease-free equilibrium as the reproduction number surpasses unity, ensuring that reducing the basic reproduction number below one is sufficient to eliminate the parasite without risking sudden bistable epidemic resurgences.

Understanding which biological and intervention parameters most heavily dictate transmission led the researchers to execute extensive local and normalized sensitivity analyses. The normalized forward sensitivity index revealed that the effective infection rate exhibits an index of exactly positive one, meaning that any percentage adjustment in transmission intensity produces an identical percentage change in the basic reproduction number. In contrast, the disease progression rate, which governs the speed at which exposed individuals become infectious, displayed a normalized sensitivity index of approximately 0.231 at baseline, showing a diminishing marginal impact as the latent period shortens. Bivariate sensitivity analysis evaluating the infection rate simultaneously against the vaccination rate illustrated that increasing vaccination coverage significantly widens the parameter space where the reproduction number remains below unity. However, the threshold boundary becomes nearly vertical as the transmission rate approaches roughly 0.0081, demonstrating that even universal vaccination coverage cannot suppress the basic reproduction number below one if baseline mosquito contact rates remain excessively elevated.

To evaluate parameter interactions without fixing baseline assumptions, the researchers performed a comprehensive global sensitivity analysis employing Latin Hypercube Sampling paired with Partial Rank Correlation Coefficients across 1,000 multidimensional parameter iterations varied within twenty-five percent of baseline values. Among parameters driving increased transmission, the human recruitment rate and the effective infection rate yielded the highest positive correlations with the reproduction number, recording coefficients of positive 0.929 and positive 0.926, respectively. The relative infectiousness of vaccinated individuals experiencing breakthrough infections ranked next with a positive correlation of 0.726, followed by the disease progression rate at positive 0.504. Among protective intervention levers, the vaccination rate produced a notable negative correlation of negative 0.395, while enhanced recovery rates among vaccinated and unvaccinated individuals exhibited minor negative correlations of negative 0.094 and negative 0.062. The waning rate of vaccine immunity had a modest negative coefficient of negative 0.013. These correlations underline that controlling mosquito biting rates and breakthrough infectiousness yields far greater epidemiological leverage than altering recovery dynamics.

The analytical theorems and sensitivity rankings were validated through comprehensive numerical simulations executed via Python’s scientific integration libraries over a 300-day observational window. Under the baseline setting where the infection rate was 0.001 and the reproduction number was approximately 0.1737, populations initialized with active malaria cases saw exposed and infectious compartments rapidly diminish toward zero while susceptible groups stabilized at their disease-free equilibrium values. When researchers simulated high-transmission endemic environments by elevating the effective infection rate to 0.015, yielding a reproduction number of roughly 2.605, the disease-related compartments surged before settling into persistent positive endemic plateaus. Incremental adjustments across five distinct infection rates confirmed that higher transmission rates produce substantially larger and earlier infection peaks, whereas varying progression rates primarily modulated the speed of the transition from latency to active infectiousness. Three-dimensional phase portraits tracking susceptible, exposed, and infected populations across diverse initial distributions showed all trajectories converging to the identical endemic coordinate, corroborating the global stability proof.

While the mathematical framework provides rigorous insights into malaria dynamics, the authors acknowledged several methodological limitations inherent to deterministic compartmental modeling. The model presumes homogeneous mixing between human hosts and mosquito vectors, omitting spatial heterogeneity, human migration patterns, and localized variations in mosquito biting behavior. Furthermore, the ordinary differential equations utilize integer-order derivatives that cannot capture memory-dependent immune responses or historical lag effects, unlike emerging fractional-order formulations. The system also omits climatic variables and seasonal fluctuations in mosquito abundance, and it assumes identical disease-induced mortality rates between vaccinated and unvaccinated patients to preserve tractability. Crucially, the authors emphasized that the baseline parameter values were compiled from existing literature and plausible biological ranges rather than calibrated directly against regional surveillance datasets, meaning numerical outputs illustrate theoretical dynamics rather than direct epidemiological forecasts. Consequently, future research will aim to incorporate field surveillance calibration, age-structured demographics, and vector population dynamics to better guide public health policy.

Subject of Research: Mathematical modeling of partial malaria vaccination and vector control

Article Title: Modeling the public health impact of partial vaccination and vector control strategies on malaria transmission

Article References: Chotaliya, K., Chavada, A., Pathak, N., & Singh, A. (2026). Modeling the public health impact of partial vaccination and vector control strategies on malaria transmission. Discover Public Health, 23(1), Article 1735. https://doi.org/10.1186/s12982-026-02937-4

Image Credits: AI Generated

DOI: 10.1186/s12982-026-02937-4

Keywords: malaria transmission, mathematical modeling, partial vaccination, vector control, basic reproduction number, ordinary differential equations, breakthrough infections, waning immunity, forward bifurcation, sensitivity analysis, Latin hypercube sampling, partial rank correlation coefficients, Plasmodium falciparum, Anopheles mosquitoes, disease-free equilibrium, endemic equilibrium, insecticide-treated nets, indoor residual spraying