Application of mathematical modelling to inform national malaria intervention planning in Nigeria

Background: For their 2021–2025 National Malaria Strategic Plan (NMSP), Nigeria’s National Malaria Elimination Programme (NMEP), in partnership with the World Health Organization (WHO), developed a targeted approach to intervention deployment at the local government area (LGA) level as part of the High Burden to High Impact response. Mathematical models of malaria transmission were used to predict the impact of proposed intervention strategies on malaria burden. Methods: An agent-based model of Plasmodium falciparum transmission was used to simulate malaria morbidity and mortality in Nigeria’s 774 LGAs under four possible intervention strategies from 2020 to 2030. The scenarios represented the previously implemented plan (business-as-usual), the NMSP at an 80% or higher coverage level and two prioritized plans according to the resources available to Nigeria. LGAs were clustered into 22 epidemiological archetypes using monthly rainfall, temperature suitability index, vector abundance, pre-2010 parasite prevalence, and pre-2010 vector control coverage. Routine incidence data were used to parameterize seasonality in each archetype. Each LGA’s baseline malaria transmission intensity was calibrated to parasite prevalence in children under the age of five years measured in the 2010 Malaria Indicator Survey (MIS). Intervention coverage in the 2010–2019 period was obtained from the Demographic and Health Survey, MIS, the NMEP, and post-campaign surveys. Results: Pursuing a business-as-usual strategy was projected to result in a 5% and 9% increase in malaria incidence in 2025 and 2030 compared with 2020, while deaths were projected to remain unchanged by 2030. The greatest intervention impact was associated with the NMSP scenario with 80% or greater coverage of standard interventions coupled with intermittent preventive treatment in infants and extension of seasonal malaria chemoprevention (SMC) to 404 LGAs, compared to 80 LGAs in 2019. The budget-prioritized scenario with SMC expansion to 310 LGAs, high bed net coverage with new formulations, and increase in effective case management rate at the same pace as historical levels was adopted as an adequate alternative for the resources available. Conclusions: Dynamical models can be applied for relative assessment of the impact of intervention scenarios but improved subnational data collection systems are required to allow increased confidence in predictions at sub-national level.