Multiscale model for forecasting Sabin 2 vaccine virus household and community transmission

Review of Poliovirus Transmission and Economic Modeling to Support Global Polio Eradication: 2020-2024

Continued investment in the development and application of mathematical models of poliovirus transmission, economics, and risks leads to their use in support of polio endgame strategy development and risk management policies. This study complements an earlier review covering the period 2000-2019 and discusses the evolution of studies published since 2020 by modeling groups supported by the Global Polio Eradication Initiative (GPEI) partners and others. We systematically review modeling papers published in English in peer-reviewed journals from 2020-2024.25 that focus on poliovirus transmission and health economic analyses. In spite of the long-anticipated end of poliovirus transmission and the GPEI sunset, which would lead to the end of its support for modeling, we find that the number of modeling groups supported by GPEI partners doubled and the rate of their publications increased. Modeling continued to play a role in supporting GPEI and national/regional policies, but changes in polio eradication governance, decentralized management and decision-making, and increased heterogeneity in modeling approaches and findings decreased the overall impact of modeling results. Meanwhile, the failure of the 2016 globally coordinated cessation of type 2 oral poliovirus vaccine use for preventive immunization and the introduction of new poliovirus vaccines and formulation, increased the complexity and uncertainty of poliovirus transmission and economic models and policy recommendations during this time.

From vaccine to pathogen: Modeling Sabin 2 vaccine virus reversion and evolutionary epidemiology in Matlab, Bangladesh

The oral poliovirus vaccines (OPVs) are one of the most effective disease eradication tools in public health. However, the OPV strains are genetically unstable and can cause outbreaks of circulating, vaccine-derived Type 2 poliovirus (cVDPV2) that are clinically indistinguishable from wild poliovirus (WPV) outbreaks. Here, we developed a Sabin 2 reversion model that simulates the reversion of Sabin 2 to reacquire a WPV-like phenotype based on the clinical differences in shedding duration and infectiousness between individuals vaccinated with Sabin 2 and those infected with WPV. Genetic reversion is informed by a canonical reversion pathway defined by three gatekeeper mutations (A481G, U2909C, and U398C) and the accumulation of deleterious nonsynonymous mutations. Our model captures essential aspects of both phenotypic and molecular evolution and simulates transmission using a multiscale transmission model that consolidates the relationships among immunity, susceptibility, and transmission risk. Despite rapid Sabin 2 attenuation reversal, we show that the emergence of a revertant virus does not guarantee a cVDPV2 outbreak. When simulating outbreaks in Matlab, Bangladesh, we found that cVDPV2 outbreaks are most likely in areas with low population-level immunity and poor sanitation. In Matlab, our model predicted that declining immunity against Type 2 poliovirus following the cessation of routine OPV vaccination was not enough to promote cVDPV2 emergence. However, cVDPV2 emergencedepended on the average viral exposure dose per contact, which was modeled as a combination of the viral concentration per fecal gram and the average fecal-oral dose per contact. These results suggest that cVDPV2 emergence risk can be mitigated by reducing the amount of infectious fecal material individuals are exposed to. Thus, a combined strategy of assessing and improving sanitation levels in conjunction with high-coverage vaccination campaigns could limit the future cVDPV2 emergence.

Novel Network Method Major Minor Variation Clustering Enables Identification of Poliovirus Clusters with High-Resolution Linkages

The Global Polio Eradication Initiative uses an outbreak response protocol that defines type 2 Sabin or Sabin-like virus as those with 0-5 nucleotides diverging from their parental strain in the complete VP1 genomic region. Sabin or Sabin-like viruses share highly similar genome sequences, regardless of their origin. Thus, it is challenging to distinguish viruses at a higher resolution to detect polio clusters or trace sources for local transmissions of viruses at an early stage. To identify type 2 Sabin or Sabin-like sources and improve our ability to map viral sources to campaigns during the polio endgame, we investigated the feasibility of a new method for genetic sequence analysis. We named the method Major Minor Variation Clustering (MMVC), which uses a network model to simultaneously incorporate sequence similarity in major and minor variants in addition to onset dates to detect fine-scale polio clusters. Each identified cluster represents a collection of sequences that are highly similar in both major and minor variants, enabling the discovery of new links between viruses. By applying the method to a published data set collected in Nigeria during 2009-2012, we found that clusters identified using this method have several improvements over clusters derived from a phylogenetic tree approach. Integrative data analysis reveals that sequences in the same cluster have greater genomic similarities and better agreement with onset dates. As a complement to current phylogenetic tree approaches, MMVC has the potential to improve epidemiological surveillance and investigation precision to guide polio eradication.

A meta-epidemiological assessment of transparency indicators of infectious disease models

Mathematical models have become very influential, especially during the COVID-19 pandemic. Data and code sharing are indispensable for reproducing them, protocol registration may be useful sometimes, and declarations of conflicts of interest (COIs) and of funding are quintessential for transparency. Here, we evaluated these features in publications of infectious disease-related models and assessed whether there were differences before and during the COVID-19 pandemic and for COVID-19 models versus models for other diseases. We analysed all PubMed Central open access publications of infectious disease models published in 2019 and 2021 using previously validated text mining algorithms of transparency indicators. We evaluated 1338 articles: 216 from 2019 and 1122 from 2021 (of which 818 were on COVID-19); almost a six-fold increase in publications within the field. 511 (39.2%) were compartmental models, 337 (25.2%) were time series, 279 (20.9%) were spatiotemporal, 186 (13.9%) were agent-based and 25 (1.9%) contained multiple model types. 288 (21.5%) articles shared code, 332 (24.8%) shared data, 6 (0.4%) were registered, and 1197 (89.5%) and 1109 (82.9%) contained COI and funding statements, respectively. There was no major changes in transparency indicators between 2019 and 2021. COVID-19 articles were less likely to have funding statements and more likely to share code. Further validation was performed by manual assessment of 10% of the articles identified by text mining as fulfilling transparency indicators and of 10% of the articles lacking them. Correcting estimates for validation performance, 26.0% of papers shared code and 41.1% shared data. On manual assessment, 5/6 articles identified as registered had indeed been registered. Of articles containing COI and funding statements, 95.8% disclosed no conflict and 11.7% reported no funding. Transparency in infectious disease modelling is relatively low, especially for data and code sharing. This is concerning, considering the nature of this research and the heightened influence it has acquired.

Integration of Heterogeneous Biological Data in Multiscale Mechanistic Model Calibration: Application to Lung Adenocarcinoma

Mechanistic models are built using knowledge as the primary information source, with well-established biological and physical laws determining the causal relationships within the model. Once the causal structure of the model is determined, parameters must be defined in order to accurately reproduce relevant data. Determining parameters and their values is particularly challenging in the case of models of pathophysiology, for which data for calibration is sparse. Multiple data sources might be required, and data may not be in a uniform or desirable format. We describe a calibration strategy to address the challenges of scarcity and heterogeneity of calibration data. Our strategy focuses on parameters whose initial values cannot be easily derived from the literature, and our goal is to determine the values of these parameters via calibration with constraints set by relevant data. When combined with a covariance matrix adaptation evolution strategy (CMA-ES), this step-by-step approach can be applied to a wide range of biological models. We describe a stepwise, integrative and iterative approach to multiscale mechanistic model calibration, and provide an example of calibrating a pathophysiological lung adenocarcinoma model. Using the approach described here we illustrate the successful calibration of a complex knowledge-based mechanistic model using only the limited heterogeneous datasets publicly available in the literature.