Covasim: An agent-based model of COVID-19 dynamics and interventions

A data science pipeline applied to Australia's 2022 COVID-19 Omicron waves

The field of software engineering is advancing at astonishing speed, with packages now available to support many stages of data science pipelines. These packages can support infectious disease modelling to be more robust, efficient and transparent, which has been particularly important during the COVID-19 pandemic. We developed a package for the construction of infectious disease models, integrated it with several open-source libraries and applied this composite pipeline to multiple data sources that provided insights into Australia's 2022 COVID-19 epidemic. We aimed to identify the key processes relevant to COVID-19 transmission dynamics and thereby develop a model that could quantify relevant epidemiological parameters. The pipeline's advantages include markedly increased speed, an expressive application programming interface, the transparency of open-source development, easy access to a broad range of calibration and optimisation tools and consideration of the full workflow from input manipulation through to algorithmic generation of the publication materials. Extending the base model to include mobility effects slightly improved model fit to data, with this approach selected as the model configuration for further epidemiological inference. Under our assumption of widespread immunity against severe outcomes from recent vaccination, incorporating an additional effect of the main vaccination programs rolled out during 2022 on transmission did not further improve model fit. Our simulations suggested that one in every two to six COVID-19 episodes were detected, subsequently emerging Omicron subvariants escaped 30-60% of recently acquired natural immunity and that natural immunity lasted only one to eight months on average. We documented our analyses algorithmically and present our methods in conjunction with interactive online code notebooks and plots. We demonstrate the feasibility of integrating a flexible domain-specific syntax library with state-of-the-art packages in high performance computing, calibration, optimisation and visualisation to create an end-to-end pipeline for infectious disease modelling. We used the resulting platform to demonstrate key epidemiological characteristics of the transition from the emergency to the endemic phase of the COVID-19 pandemic.

Investigating the impact of vaccine hesitancy on an emerging infectious disease: a mathematical and numerical analysis

Throughout the last two centuries, vaccines have been helpful in mitigating numerous epidemic diseases. However, vaccine hesitancy has been identified as a substantial obstacle in healthcare management. We examined the epidemiological dynamics of an emerging infection under vaccination using an SVEIR model with differential morbidity. We mathematically analyzed the model, derived R 0 , and provided a complete analysis of the bifurcation at R 0 = 1 . Sensitivity analysis and numerical simulations were used to quantify the tradeoffs between vaccine efficacy and vaccine hesitancy on reducing the disease burden. Our results indicated that if the percentage of the population hesitant about taking the vaccine is 10%, then a vaccine with 94% efficacy is required to reduce the peak of infections by 40%. If 60% of the population is reluctant about being vaccinated, then even a perfect vaccine will not be able to reduce the peak of infections by 40%.

Learning from the COVID-19 pandemic: A systematic review of mathematical vaccine prioritization models

As the world becomes ever more connected, the chance of pandemics increases as well. The recent COVID-19 pandemic and the concurrent global mass vaccine roll-out provides an ideal setting to learn from and refine our understanding of infectious disease models for better future preparedness. In this review, we systematically analyze and categorize mathematical models that have been developed to design optimal vaccine prioritization strategies of an initially limited vaccine. As older individuals are disproportionately affected by COVID-19, the focus is on models that take age explicitly into account. The lower mobility and activity level of older individuals gives rise to non-trivial trade-offs. Secondary research questions concern the optimal time interval between vaccine doses and spatial vaccine distribution. This review showcases the effect of various modeling assumptions on model outcomes. A solid understanding of these relationships yields better infectious disease models and thus public health decisions during the next pandemic.

Admission screening testing of patients and staff N95 respirators are cost-effective in reducing COVID-19 hospital-acquired infections

BACKGROUND

Coronavirus disease 2019 (COVID-19) outbreaks in acute care settings can have severe consequences for patients due to their underlying vulnerabilities, and can be costly due to additional patient bed-days and the need to replace isolating staff. This study assessed the cost-effectiveness of clinical staff N95 respirators and admission screening testing of patients to reduce COVID-19 hospital-acquired infections.

METHODS

An agent-based model was calibrated to data on 178 outbreaks in acute care settings in Victoria, Australia between October 2021 and July 2023. Outbreaks were simulated under different combinations of staff masking (surgical, N95) and patient admission screening testing [none, rapid antigen test (RAT), polymerase chain reaction]. For each scenario, average diagnoses, COVID-19 deaths, quality-adjusted life years from discharged patients, and costs (masks, testing, patient COVID-19 bed-days, staff replacement costs while isolating) from acute COVID-19 were estimated over a 12-month period.

FINDINGS

Compared with no admission screening testing and staff surgical masks, all scenarios were cost saving with health gains. Staff N95 respirators + RAT admission screening of patients was the cheapest scenario, saving A$78.4M [95% uncertainty interval (UI) 44.4M-135.3M] and preventing 1543 (95% UI 1070-2146) deaths state-wide per annum. Both interventions were individually beneficial: staff N95 respirators saved A$54.7M and 854 deaths state-wide per annum, while RAT admission screening of patients saved A$57.6M and 1176 deaths state-wide per annum.

INTERPRETATION

In acute care settings, staff N95 respirators and admission screening testing of patients can reduce hospital-acquired COVID-19 and COVID-19 deaths, and are cost saving because of reduced patient bed-days and staff replacement needs.

Simulation to optimize the laboratory diagnosis of bacteremia

Blood cultures are central to the management of patients with sepsis and bloodstream infection. Clinical decisions depend on the timely availability of laboratory information, which, in turn, depends on the optimal laboratory processing of specimens. Discrete event simulation (DES) offers insights into where optimization efforts can be targeted. Here, we generate a detailed process map of blood culture processing within a laboratory and use it to build a simulator. Direct observation of laboratory staff processing blood cultures was used to generate a flowchart of the blood culture laboratory pathway. Retrospective routinely collected data were combined with direct observations to generate probability distributions over the time taken for each event. These data were used to inform the DES model. A sensitivity analysis explored the impact of staff availability on turnaround times. A flowchart of the blood culture pathway was constructed, spanning labeling, incubation, organism identification, and antimicrobial susceptibility testing. Thirteen processes in earlier stages of the pathway, not otherwise captured by routinely collected data, were timed using direct observations. Observations revealed that specimen processing is predominantly batched. Another eight processes were timed using retrospective data. A simulator was built using DES. Sensitivity analysis revealed that specimen progression through the simulation was especially sensitive to laboratory technician availability. Gram stain reporting time was also sensitive to laboratory scientist availability. Our laboratory simulation model has wide-ranging applications for the optimization of laboratory processes and effective implementation of the changes required for faster and more accurate results.

IMPORTANCE

Optimization of laboratory pathways and resource availability has a direct impact on the clinical management of patients with bloodstream infection. This research offers an insight into the laboratory processing of blood cultures at a system level and allows clinical microbiology laboratories to explore the impact of changes to processes and resources.

Effect of vaccine dose intervals: Considering immunity levels, vaccine efficacy, and strain variants for disease control strategy

In this study, we present an immuno-epidemic model to understand mitigation options during an epidemic break. The model incorporates comorbidity and multiple-vaccine doses through a system of coupled integro-differential equations to analyze the epidemic rate and intensity from a knowledge of the basic reproduction number and time-distributed rate functions. Our modeling results show that the interval between vaccine doses is a key control parameter that can be tuned to significantly influence disease spread. We show that multiple doses induce a hysteresis effect in immunity levels that offers a better mitigation alternative compared to frequent vaccination which is less cost-effective while being more intrusive. Optimal dosing intervals, emphasizing the cost-effectiveness of each vaccination effort, and determined by various factors such as the level of immunity and efficacy of vaccines against different strains, appear to be crucial in disease management. The model is sufficiently generic that can be extended to accommodate specific disease forms.

Mathematical assessment of the roles of age heterogeneity and vaccination on the dynamics and control of SARS-CoV-2

The COVID-19 pandemic, caused by SARS-CoV-2, disproportionately affected certain segments of society, particularly the elderly population (which suffered the brunt of the burden of the pandemic in terms of severity of the disease, hospitalization, and death). This study presents a generalized multigroup model, with m heterogeneous sub-populations, to assess the population-level impact of age heterogeneity and vaccination on the transmission dynamics and control of the SARS-CoV-2 pandemic in the United States. Rigorous analysis of the model for the homogeneous case (i.e., the model with m = 1) reveal that its disease-free equilibrium is globally-asymptotically stable for two special cases (with perfect vaccine efficacy or negligible disease-induced mortality) whenever the associated reproduction number is less than one. The model has a unique and globally-asymptotically stable endemic equilibrium, for special a case, when the associated reproduction threshold exceeds one. The homogeneous model was fitted using the observed cumulative mortality data for the United States during three distinct waves (Waves A (October 17, 2020 to April 5, 2021), B (July 9, 2021 to November 7, 2021) and C (January 1, 2022 to May 7, 2022)) chosen to align with time periods when the Alpha, Delta and Omicron were, respectively, the predominant variants in the United States. The calibrated model was used to derive a theoretical expression for achieving vaccine-derived herd immunity (needed to eliminate the disease in the United States). It was shown that, using the one-group homogeneous model, vaccine-derived herd immunity is not attainable during Wave C of the pandemic in the United States, regardless of the coverage level of the fully-vaccinated individuals. Global sensitivity analysis was carried out to determine the parameters of the model that have the most influence on the disease dynamics and burden. These analyses reveal that control and mitigation strategies that may be very effective during one wave may not be so very effective during the other wave or waves. However, strategies that target asymptomatic and pre-symptomatic infectious individuals are shown to be consistently effective across all waves. To study the impact of the disproportionate effect of COVID-19 on the elderly population, we considered the heterogeneous model for the case where the total population is subdivided into the sub-populations of individuals under 65 years of age and those that are 65 and older. The resulting two-group heterogeneous model, which was also fitted using the cumulative mortality data for wave C, was also rigorously analysed. Unlike for the case of the one-group model, it was shown, for the two-group model, that vaccine-derived herd immunity can indeed be achieved during Wave C of the pandemic if at least 61% of the populace is fully vaccinated. Thus, this study shows that adding age heterogeneity into a SARS-CoV-2 vaccination model with homogeneous mixing significantly reduces the level of vaccination coverage needed to achieve vaccine-derived herd immunity (specifically, for the heterogeneous model, herd-immunity can be attained during Wave C if a moderate proportion of susceptible individuals are fully vaccinated). The consequence of this result is that vaccination models for SARS-CoV-2 that do not explicitly account for age heterogeneity may be overestimating the level of vaccine-derived herd immunity threshold needed to eliminate the SARS-CoV-2 pandemic.

Effectiveness of probabilistic contact tracing in epidemic containment: The role of superspreaders and transmission path reconstruction

The recent COVID-19 pandemic underscores the significance of early stage nonpharmacological intervention strategies. The widespread use of masks and the systematic implementation of contact tracing strategies provide a potentially equally effective and socially less impactful alternative to more conventional approaches, such as large-scale mobility restrictions. However, manual contact tracing faces strong limitations in accessing the network of contacts, and the scalability of currently implemented protocols for smartphone-based digital contact tracing becomes impractical during the rapid expansion phases of the outbreaks, due to the surge in exposure notifications and associated tests. A substantial improvement in digital contact tracing can be obtained through the integration of probabilistic techniques for risk assessment that can more effectively guide the allocation of diagnostic tests. In this study, we first quantitatively analyze the diagnostic and social costs associated with these containment measures based on contact tracing, employing three state-of-the-art models of SARS-CoV-2 spreading. Our results suggest that probabilistic techniques allow for more effective mitigation at a lower cost. Secondly, our findings reveal a remarkable efficacy of probabilistic contact-tracing techniques in performing backward and multistep tracing and capturing superspreading events.

Role of heterogeneity: National scale data-driven agent-based modeling for the US COVID-19 Scenario Modeling Hub

UVA-EpiHiper is a national scale agent-based model to support the US COVID-19 Scenario Modeling Hub (SMH). UVA-EpiHiper uses a detailed representation of the underlying social contact network along with data measured during the course of the pandemic to initialize and calibrate the model. In this paper, we study the role of heterogeneity on model complexity and resulting epidemic dynamics using UVA-EpiHiper. We discuss various sources of heterogeneity that we encounter in the use of UVA-EpiHiper to support modeling and analysis of epidemic dynamics under various scenarios. We also discuss how this affects model complexity and computational complexity of the corresponding simulations. Using round 13 of the SMH as an example, we discuss how UVA-EpiHiper was initialized and calibrated. We then discuss how the detailed output produced by UVA-EpiHiper can be analyzed to obtain interesting insights. We find that despite the complexity in the model, the software, and the computation incurred to an agent-based model in scenario modeling, it is capable of capturing various heterogeneities of real-world systems, especially those in networks and behaviors, and enables analyzing heterogeneities in epidemiological outcomes between different demographic, geographic, and social cohorts. In applying UVA-EpiHiper to round 13 scenario modeling, we find that disease outcomes are different between and within states, and between demographic groups, which can be attributed to heterogeneities in population demographics, network structures, and initial immunity.

Dynamic Contact Networks in Confined Spaces: Synthesizing Micro-Level Encounter Patterns through Human Mobility Models from Real-World Data

This study advances the field of infectious disease forecasting by introducing a novel approach to micro-level contact modeling, leveraging human movement patterns to generate realistic temporal-dynamic networks. Through the incorporation of human mobility models and parameter tuning, this research presents an innovative method for simulating micro-level encounters that closely mirror infection dynamics within confined spaces. Central to our methodology is the application of Bayesian optimization for parameter selection, which refines our models to emulate both the properties of real-world infection curves and the characteristics of network properties. Typically, large-scale epidemiological simulations overlook the specifics of human mobility within confined spaces or rely on overly simplistic models. By focusing on the distinct aspects of infection propagation within specific locations, our approach strengthens the realism of such pandemic simulations. The resulting models shed light on the role of spatial encounters in disease spread and improve the capability to forecast and respond to infectious disease outbreaks. This work not only contributes to the scientific understanding of micro-level transmission patterns but also offers a new perspective on temporal network generation for epidemiological modeling.

Drivers of epidemic dynamics in real time from daily digital COVID-19 measurements

Understanding the drivers of respiratory pathogen spread is challenging, particularly in a timely manner during an ongoing epidemic. In this work, we present insights that we obtained using daily data from the National Health Service COVID-19 app for England and Wales and that we shared with health authorities in almost real time. Our indicator of the reproduction number R(t) was available days earlier than other estimates, with an innovative capability to decompose R(t) into contact rates and probabilities of infection. When Omicron arrived, the main epidemic driver switched from contacts to transmissibility. We separated contacts and transmissions by day of exposure and setting and found pronounced variability over days of the week and during Christmas holidays and events. For example, during the Euro football tournament in 2021, days with England matches showed sharp spikes in exposures and transmissibility. Digital contact-tracing technologies can help control epidemics not only by directly preventing transmissions but also by enabling rapid analysis at scale and with unprecedented resolution.

Intervention effect of targeted workplace closures may be approximated by single-layered networks in an individual-based model of COVID-19 control

Individual-based models of infectious disease dynamics commonly use network structures to represent human interactions. Network structures can vary in complexity, from single-layered with homogeneous mixing to multi-layered with clustering and layer-specific contact weights. Here we assessed policy-relevant consequences of network choice by simulating different network structures within an established individual-based model of SARS-CoV-2 dynamics. We determined the clustering coefficient of each network structure and compared this to several epidemiological outcomes, such as cumulative and peak infections. High-clustered networks estimate fewer cumulative infections and peak infections than less-clustered networks when transmission probabilities are equal. However, by altering transmission probabilities, we find that high-clustered networks can essentially recover the dynamics of low-clustered networks. We further assessed the effect of workplace closures as a layer-targeted intervention on epidemiological outcomes and found in this scenario a single-layered network provides a sufficient approximation of intervention effect relative to a multi-layered network when layer-specific contact weightings are equal. Overall, network structure choice within models should consider the knowledge of contact weights in different environments and pathogen mode of transmission to avoid over- or under-estimating disease burden and impact of interventions.

Evaluating equity-promoting interventions to prevent race-based inequities in influenza outcomes

Importance

Seasonal influenza hospitalizations pose a considerable burden in the United States, with BIPOC (Black, Indigenous, and other People of Color) communities being disproportionately affected.

Objective

To determine and quantify the effects of different types of mitigation strategies on inequities in influenza outcomes (symptomatic infections and hospitalizations).

Design

In this simulation study, we fit a race-stratified agent-based model of influenza transmission to demographic and hospitalization data of the United States.

Participants

We consider five racial-ethnic groups: non-Hispanic White persons, non- Hispanic Black persons, non-Hispanic Asian persons, non-Hispanic American Indian or Alaska Native persons, and Hispanic or Latino persons.

Setting

We tested five idealized equity-promoting interventions to determine their effectiveness in reducing inequity in influenza outcomes. The interventions assumed (i) equalized vaccination rates, (ii) equalized comorbidities, (iii) work-risk distribution proportional to the distribution of the population, (iv) reduced work contacts for all, or (v) a combination of equalizing vaccination rates and comorbidities and reducing work contacts.

Main Outcomes and Measures

Reduction in symptomatic or hospitalization risk ratios, defined as the ratio of the number of symptomatic infections (hospitalizations respectively) in each age- and racial-ethnic group and their corresponding white counterpart. We also evaluated the reduction in the absolute mean number of symptomatic infections or hospitalizations in each age- and racial-ethnic group compared to the fitted scenario (baseline).

Results

Our analysis suggests that symptomatic infections were equalized and reduced (by up to 17% in BIPOC adults aged 18-49) by strategies reducing work contacts or equalizing vaccination rates. Reducing comorbidities resulted in significant decreases in hospitalizations, with a reduction of over 40% in BIPOC groups. All tested interventions reduced the inequity in influenza hospitalizations in all racial-ethnic groups, but interventions reducing comorbidities in marginalized populations were the most effective. Notably, these interventions resulted in better outcomes across all racial-ethnic groups, not only those prioritized by the interventions.

Conclusions and Relevance

In this simulation modeling study, equalizing vaccination rates and reducing number of work contacts (which are relatively simple strategies to implement) reduced the both the inequity in hospitalizations and the absolute number of symptomatic infections and hospitalizations in all age and racial-ethnic groups.

Analysis of the impact of COVID-19 variants and vaccination on the time-varying reproduction number: statistical methods

Introduction

The COVID-19 pandemic has profoundly impacted global health systems, requiring the monitoring of infection waves and strategies to control transmission. Estimating the time-varying reproduction number is crucial for understanding the epidemic and guiding interventions.

Methods

Probability distributions of serial interval are estimated for Pre-Delta and Delta periods. We conducted a comparative analysis of time-varying reproduction numbers, taking into account population immunity and variant differences. We incorporated the regional heterogeneity and age distribution of the population, as well as the evolving variants and vaccination rates over time. COVID-19 transmission dynamics were analyzed with variants and vaccination.

Results

The reproduction number is computed with and without considering variant-based immunity. In addition, values of reproduction number significantly differed by variants, emphasizing immunity's importance. Enhanced vaccination efforts and stringent control measures were effective in reducing the transmission of the Delta variant. Conversely, Pre-Delta variant appeared less influenced by immunity levels, due to lower vaccination rates. Furthermore, during the Pre-Delta period, there was a significant difference between the region-specific and the non-region-specific reproduction numbers, with particularly distinct pattern differences observed in Gangwon, Gyeongbuk, and Jeju in Korea.

Discussion

This research elucidates the dynamics of COVID-19 transmission concerning the dominance of the Delta variant, the efficacy of vaccinations, and the influence of immunity levels. It highlights the necessity for targeted interventions and extensive vaccination coverage. This study makes a significant contribution to the understanding of disease transmission mechanisms and informs public health strategies.

Model-driven decision support: A community-based meta-implementation strategy to predict population impact

PURPOSE

Standard tools for public health decision making such as data dashboards, trial repositories, and intervention briefs may be necessary but insufficient for guiding community leaders in optimizing local public health strategy. Predictive modeling decision support tools may be the missing link that allows community level decision makers to confidently direct funding and other resources to interventions and implementation strategies that will improve upon the status quo.

METHODS

We describe a community-based model-driven decision support (MDDS) approach that requires community engagement, local data, and predictive modeling tools (agent-based modeling in our case studies) to improve decision-making on implementing strategies to address complex public health problems such as overdose deaths. We refer to our approach as a meta-implementation strategy as it provides guidance to a community on what intervention combinations and their required implementation strategies are needed to achieve desired outcomes. We use standard implementation measures including the Stages of Implementation Completion to assess adoption of this meta-implementation approach.

RESULTS

Using two case studies, we illustrate how MDDS can be used to support decision making related to HIV prevention and reductions in overdose deaths at the city and county level. Even when community acceptance seems high, data acquisition and diffuse responsibility for implementing specific strategies recommended by modeling are barriers to adoption.

CONCLUSIONS

MDDS has the capacity to improve community decision makers use of scientific knowledge by providing projections of the impact of intervention strategies under various scenarios. Further research is necessary to assess its effectiveness and the best strategies to implement it.

Leveraging artificial intelligence and software engineering methods in epidemiology for the co-creation of decision-support tools based on mechanistic models

Epidemiological modeling is a key lever for infectious disease control and prevention on farms. It makes it possible to understand the spread of pathogens, but also to compare intervention scenarios even in counterfactual situations. However, the actual capability of decision makers to use mechanistic models to support timely interventions is limited. This study demonstrates how artificial intelligence (AI) techniques can make mechanistic epidemiological models more accessible to farmers and veterinarians, and how to transform such models into user-friendly decision-support tools (DST). By leveraging knowledge representation methods, such as the textual formalization of model components through a domain-specific language (DSL), the co-design of mechanistic models and DST becomes more efficient and collaborative. This facilitates the integration of explicit expert knowledge and practical insights into the modeling process. Furthermore, the utilization of AI and software engineering enables the automation of web application generation based on existing mechanistic models. This automation simplifies the development of DST, as tool designers can focus on identifying users' needs and specifying expected features and meaningful presentations of outcomes, instead of wasting time in writing code to wrap models into web apps. To illustrate the practical application of this approach, we consider the example of Bovine Respiratory Disease (BRD), a tough challenge in fattening farms where young beef bulls often develop BRD shortly after being allocated into pens. BRD is a multi-factorial, multi-pathogen disease that is difficult to anticipate and control, often resulting in the massive use of antimicrobials to mitigate its impact on animal health, welfare, and economic losses. The DST developed from an existing mechanistic BRD model empowers users, including farmers and veterinarians, to customize scenarios based on their specific farm conditions. It enables them to anticipate the effects of various pathogens, compare the epidemiological and economic outcomes associated with different farming practices, and decide how to balance the reduction of disease impact and the reduction of antimicrobial usage (AMU). The generic method presented in this article illustrates the potential of artificial intelligence (AI) and software engineering methods to enhance the co-creation of DST based on mechanistic models in veterinary epidemiology. The corresponding pipeline is distributed as an open-source software. By leveraging these advancements, this research aims to bridge the gap between theoretical models and the practical usage of their outcomes on the field.

HPVsim: An agent-based model of HPV transmission and cervical disease

In 2020, the WHO launched its first global strategy to accelerate the elimination of cervical cancer, outlining an ambitious set of targets for countries to achieve over the next decade. At the same time, new tools, technologies, and strategies are in the pipeline that may improve screening performance, expand the reach of prophylactic vaccines, and prevent the acquisition, persistence and progression of oncogenic HPV. Detailed mechanistic modelling can help identify the combinations of current and future strategies to combat cervical cancer. Open-source modelling tools are needed to shift the capacity for such evaluations in-country. Here, we introduce the Human papillomavirus simulator (HPVsim), a new open-source software package for creating flexible agent-based models parameterised with country-specific vital dynamics, structured sexual networks, and co-transmitting HPV genotypes. HPVsim includes a novel methodology for modelling cervical disease progression, designed to be readily adaptable to new forms of screening. The software itself is implemented in Python, has built-in tools for simulating commonly-used interventions, includes a comprehensive set of tests and documentation, and runs quickly (seconds to minutes) on a laptop. Performance is greatly enhanced by HPVsim's multiscale modelling functionality. HPVsim is open source under the MIT License and available via both the Python Package Index (via pip install) and GitHub (hpvsim.org).

Agent-based modeling to estimate the impact of lockdown scenarios and events on a pandemic exemplified on SARS-CoV-2

In actual pandemic situations like COVID-19, it is important to understand the influence of single mitigation measures as well as combinations to create most dynamic impact for lockdown scenarios. Therefore we created an agent-based model (ABM) to simulate the spread of SARS-CoV-2 in an abstract city model with several types of places and agents. In comparison to infection numbers in Germany our ABM could be shown to behave similarly during the first wave. In our model, we implemented the possibility to test the effectiveness of mitigation measures and lockdown scenarios on the course of the pandemic. In this context, we focused on parameters of local events as possible mitigation measures and ran simulations, including varying size, duration, frequency and the proportion of events. The majority of changes to single event parameters, with the exception of frequency, showed only a small influence on the overall course of the pandemic. By applying different lockdown scenarios in our simulations, we could observe drastic changes in the number of infections per day. Depending on the lockdown strategy, we even observed a delayed peak in infection numbers of the second wave. As an advantage of the developed ABM, it is possible to analyze the individual risk of single agents during the pandemic. In contrast to standard or adjusted ODEs, we observed a 21% (with masks) / 48% (without masks) increased risk for single reappearing participants on local events, with a linearly increasing risk based on the length of the events.

Effectiveness of interventions to reduce COVID-19 transmission in schools

School reopenings in 2021 and 2022 coincided with the rapid emergence of new SARS-CoV-2 variants in the United States. In-school mitigation efforts varied, depending on local COVID-19 mandates and resources. Using a stochastic age-stratified agent-based model of SARS-CoV-2 transmission, we estimate the impacts of multiple in-school strategies on both infection rates and absenteeism, relative to a baseline scenario in which only symptomatic cases are tested and positive tests trigger a 10-day isolation of the case and 10-day quarantine of their household and classroom. We find that monthly asymptomatic screening coupled with the 10-day isolation and quarantine period is expected to avert 55.4% of infections while increasing absenteeism by 104.3%. Replacing quarantine with test-to-stay would reduce absenteeism by 66.3% (while hardly impacting infection rates), but would require roughly 10-fold more testing resources. Alternatively, vaccination or mask wearing by 50% of the student body is expected to avert 54.1% or 43.1% of infections while decreasing absenteeism by 34.1% or 27.4%, respectively. Separating students into classrooms based on mask usage is expected to reduce infection risks among those who wear masks (by 23.1%), exacerbate risks among those who do not (by 27.8%), but have little impact on overall risk. A combined strategy of monthly screening, household and classroom quarantine, a 50% vaccination rate, and a 50% masking rate (in mixed classrooms) is expected to avert 81.7% of infections while increasing absenteeism by 90.6%. During future public health emergencies, such analyses can inform the rapid design of resource-constrained strategies that mitigate both public health and educational risks.

Projecting the future impact of emerging SARS-CoV-2 variants under uncertainty: Modeling the initial Omicron outbreak

Over the past several years, the emergence of novel SARS-CoV-2 variants has led to multiple waves of increased COVID-19 incidence. When the Omicron variant emerged, there was considerable concern about its potential impact in the winter of 2021-2022 due to its increased fitness. However, there was also considerable uncertainty regarding its likely impact due to questions about its relative transmissibility, severity, and degree of immune escape. We sought to evaluate the ability of an agent-based model to forecast incidence in the context of this emerging pathogen variant. To project COVID-19 cases and deaths in Indiana, we calibrated our model to COVID-19 hospitalizations, deaths, and test-positivity rates through November 2021, and then projected COVID-19 incidence through April 2022 under four different scenarios that covered the plausible ranges of Omicron's severity, transmissibility, and degree of immune escape. Our initial projections from December 2021 through March 2022 indicated that under a pessimistic scenario with high disease severity, the peak in weekly COVID-19 deaths in Indiana would be larger than the previous peak in December 2020. However, retrospective analyses indicate that Omicron's severity was closer to the optimistic scenario, and even though cases and hospitalizations reached a new peak, fewer deaths occurred than during the previous peak. According to our results, Omicron's rapid spread was consistent with a combination of higher transmissibility and immune escape relative to earlier variants. Our updated projections starting in January 2022 accurately predicted that cases would peak in mid-January and decline rapidly over the next several months. The performance of our projections shows that following the emergence of a new pathogen variant, models can help quantify the potential range of outbreak magnitudes and trajectories. Agent-based models are particularly useful in these scenarios because they can efficiently track individual vaccination and infection histories with multiple variants with varying degrees of cross-protection.

flepiMoP: The evolution of a flexible infectious disease modeling pipeline during the COVID-19 pandemic

The COVID-19 pandemic led to an unprecedented demand for projections of disease burden and healthcare utilization under scenarios ranging from unmitigated spread to strict social distancing policies. In response, members of the Johns Hopkins Infectious Disease Dynamics Group developed flepiMoP (formerly called the COVID Scenario Modeling Pipeline), a comprehensive open-source software pipeline designed for creating and simulating compartmental models of infectious disease transmission and inferring parameters through these models. The framework has been used extensively to produce short-term forecasts and longer-term scenario projections of COVID-19 at the state and county level in the US, for COVID-19 in other countries at various geographic scales, and more recently for seasonal influenza. In this paper, we highlight how the flepiMoP has evolved throughout the COVID-19 pandemic to address changing epidemiological dynamics, new interventions, and shifts in policy-relevant model outputs. As the framework has reached a mature state, we provide a detailed overview of flepiMoP's key features and remaining limitations, thereby distributing flepiMoP and its documentation as a flexible and powerful tool for researchers and public health professionals to rapidly build and deploy large-scale complex infectious disease models for any pathogen and demographic setup.

COVID-19 outbreaks in residential aged care facilities: an agent-based modeling study

Introduction

A disproportionate number of COVID-19 deaths occur in Residential Aged Care Facilities (RACFs), where better evidence is needed to target COVID-19 interventions to prevent mortality. This study used an agent-based model to assess the role of community prevalence, vaccination strategies, and non-pharmaceutical interventions (NPIs) on COVID-19 outcomes in RACFs in Victoria, Australia.

Methods

The model simulated outbreaks in RACFs over time, and was calibrated to distributions for outbreak size, outbreak duration, and case fatality rate in Victorian RACFs over 2022. The number of incursions to RACFs per day were estimated to fit total deaths and diagnoses over time and community prevalence.Total infections, diagnoses, and deaths in RACFs were estimated over July 2023-June 2024 under scenarios of different: community epidemic wave assumptions (magnitude and frequency); RACF vaccination strategies (6-monthly, 12-monthly, no further vaccines); additional non-pharmaceutical interventions (10, 25, 50% efficacy); and reduction in incursions (30% or 60%).

Results

Total RACF outcomes were proportional to cumulative community infections and incursion rates, suggesting potential for strategic visitation/staff policies or community-based interventions to reduce deaths. Recency of vaccination when epidemic waves occurred was critical; compared with 6-monthly boosters, 12-monthly boosters had approximately 1.2 times more deaths and no further boosters had approximately 1.6 times more deaths over July 2023-June 2024. Additional NPIs, even with only 10-25% efficacy, could lead to a 13-31% reduction in deaths in RACFs.

Conclusion

Future community epidemic wave patterns are unknown but will be major drivers of outcomes in RACFs. Maintaining high coverage of recent vaccination, minimizing incursions, and increasing NPIs can have a major impact on cumulative infections and deaths.

Using an epidemiological model to explore the interplay between sharing and advertising in viral videos

How to exploit social networks to make internet content spread rapidly and consistently is an interesting question in marketing management. Although epidemic models have been employed to comprehend the spread dynamics of internet content, such as viral videos, the effects of advertising and individual sharing on information dissemination are difficult to distinguish. This gap forbids us to evaluate the efficiency of marketing strategies. In this paper, we modify a classic mean-field SIR (susceptible-infected-recovered) model, incorporating the influences of sharing and advertising in viral videos. We mathematically analyze the global stability of the system and propose an agent-based modeling approach to evaluate the efficiency of sharing and advertising. We further provide a case study of music videos on YouTube to show the validity of our model.

Task-oriented machine learning surrogates for tipping points of agent-based models

We present a machine learning framework bridging manifold learning, neural networks, Gaussian processes, and Equation-Free multiscale approach, for the construction of different types of effective reduced order models from detailed agent-based simulators and the systematic multiscale numerical analysis of their emergent dynamics. The specific tasks of interest here include the detection of tipping points, and the uncertainty quantification of rare events near them. Our illustrative examples are an event-driven, stochastic financial market model describing the mimetic behavior of traders, and a compartmental stochastic epidemic model on an Erdös-Rényi network. We contrast the pros and cons of the different types of surrogate models and the effort involved in learning them. Importantly, the proposed framework reveals that, around the tipping points, the emergent dynamics of both benchmark examples can be effectively described by a one-dimensional stochastic differential equation, thus revealing the intrinsic dimensionality of the normal form of the specific type of the tipping point. This allows a significant reduction in the computational cost of the tasks of interest.

Autonomous and policy-induced behavior change during the COVID-19 pandemic: Towards understanding and modeling the interplay of behavioral adaptation

Changes in human behaviors, such as reductions of physical contacts and the adoption of preventive measures, impact the transmission of infectious diseases considerably. Behavioral adaptations may be the result of individuals aiming to protect themselves or mere responses to public containment measures, or a combination of both. What drives autonomous and policy-induced adaptation, how they are related and change over time is insufficiently understood. Here, we develop a framework for more precise analysis of behavioral adaptation, focusing on confluence, interactions and time variance of autonomous and policy-induced adaptation. We carry out an empirical analysis of Germany during the fall of 2020 and beyond. Subsequently, we discuss how behavioral adaptation processes can be better represented in behavioral-epidemiological models. We find that our framework is useful to understand the interplay of autonomous and policy-induced adaptation as a "moving target". Our empirical analysis suggests that mobility patterns in Germany changed significantly due to both autonomous and policy-induced adaption, with potentially weaker effects over time due to decreasing risk signals, diminishing risk perceptions and an erosion of trust in the government. We find that while a number of simulation and prediction models have made great efforts to represent behavioral adaptation, the interplay of autonomous and policy-induced adaption needs to be better understood to construct convincing counterfactual scenarios for policy analysis. The insights presented here are of interest to modelers and policy makers aiming to understand and account for behaviors during a pandemic response more accurately.

The impact of quality-adjusted life years on evaluating COVID-19 mitigation strategies: lessons from age-specific vaccination roll-out and variants of concern in Belgium (2020-2022)

BACKGROUND

When formulating and evaluating COVID-19 vaccination strategies, an emphasis has been placed on preventing severe disease that overburdens healthcare systems and leads to mortality. However, more conventional outcomes such as quality-adjusted life years (QALYs) and inequality indicators are warranted as additional information for policymakers.

METHODS

We adopted a mathematical transmission model to describe the infectious disease dynamics of SARS-COV-2, including disease mortality and morbidity, and to evaluate (non)pharmaceutical interventions. Therefore, we considered temporal immunity levels, together with the distinct transmissibility of variants of concern (VOCs) and their corresponding vaccine effectiveness. We included both general and age-specific characteristics related to SARS-CoV-2 vaccination. Our scenario study is informed by data from Belgium, focusing on the period from August 2021 until February 2022, when vaccination for children aged 5-11 years was initially not yet licensed and first booster doses were administered to adults. More specifically, we investigated the potential impact of an earlier vaccination programme for children and increased or reduced historical adult booster dose uptake.

RESULTS

Through simulations, we demonstrate that increasing vaccine uptake in children aged 5-11 years in August-September 2021 could have led to reduced disease incidence and ICU occupancy, which was an essential indicator for implementing non-pharmaceutical interventions and maintaining healthcare system functionality. However, an enhanced booster dose regimen for adults from November 2021 onward could have resulted in more substantial cumulative QALY gains, particularly through the prevention of elevated levels of infection and disease incidence associated with the emergence of Omicron VOC. In both scenarios, the need for non-pharmaceutical interventions could have decreased, potentially boosting economic activity and mental well-being.

CONCLUSIONS

When calculating the impact of measures to mitigate disease spread in terms of life years lost due to COVID-19 mortality, we highlight the impact of COVID-19 on the health-related quality of life of survivors. Our study underscores that disease-related morbidity could constitute a significant part of the overall health burden. Our quantitative findings depend on the specific setup of the interventions under review, which is open to debate or should be contextualised within future situations.

A Framework for Assessing the Impact of Outbreak Response Immunization Programs

The impact of outbreak response immunization (ORI) can be estimated by comparing observed outcomes to modelled counterfactual scenarios without ORI, but the most appropriate metrics depend on stakeholder needs and data availability. This study developed a framework for using mathematical models to assess the impact of ORI for vaccine-preventable diseases. Framework development involved (1) the assessment of impact metrics based on stakeholder interviews and literature reviews determining data availability and capacity to capture as model outcomes; (2) mapping investment in ORI elements to model parameters to define scenarios; (3) developing a system for engaging stakeholders and formulating model questions, performing analyses, and interpreting results; and (4) example applications for different settings and pathogens. The metrics identified as most useful were health impacts, economic impacts, and the risk of severe outbreaks. Scenario categories included investment in the response scale, response speed, and vaccine targeting. The framework defines four phases: (1) problem framing and data sourcing (identification of stakeholder needs, metrics, and scenarios); (2) model choice; (3) model implementation; and (4) interpretation and communication. The use of the framework is demonstrated by application to two outbreaks, measles in Papua New Guinea and Ebola in the Democratic Republic of the Congo. The framework is a systematic way to engage with stakeholders and ensure that an analysis is fit for purpose, makes the best use of available data, and uses suitable modelling methodology.

AI-powered COVID-19 forecasting: a comprehensive comparison of advanced deep learning methods

BACKGROUND

The coronavirus disease 2019 (COVID-19) pandemic continues to pose significant challenges to the public health sector, including that of the United Arab Emirates (UAE). The objective of this study was to assess the efficiency and accuracy of various deep-learning models in forecasting COVID-19 cases within the UAE, thereby aiding the nation's public health authorities in informed decision-making.

METHODS

This study utilized a comprehensive dataset encompassing confirmed COVID-19 cases, demographic statistics, and socioeconomic indicators. Several advanced deep learning models, including long short-term memory (LSTM), bidirectional LSTM, convolutional neural network (CNN), CNN-LSTM, multilayer perceptron, and recurrent neural network (RNN) models, were trained and evaluated. Bayesian optimization was also implemented to fine-tune these models.

RESULTS

The evaluation framework revealed that each model exhibited different levels of predictive accuracy and precision. Specifically, the RNN model outperformed the other architectures even without optimization. Comprehensive predictive and perspective analytics were conducted to scrutinize the COVID-19 dataset.

CONCLUSION

This study transcends academic boundaries by offering critical insights that enable public health authorities in the UAE to deploy targeted data-driven interventions. The RNN model, which was identified as the most reliable and accurate for this specific context, can significantly influence public health decisions. Moreover, the broader implications of this research validate the capability of deep learning techniques in handling complex datasets, thus offering the transformative potential for predictive accuracy in the public health and healthcare sectors.

Learning from the COVID-19 pandemic: a systematic review of mathematical vaccine prioritization models

As the world becomes ever more connected, the chance of pandemics increases as well. The recent COVID-19 pandemic and the concurrent global mass vaccine roll-out provides an ideal setting to learn from and refine our understanding of infectious disease models for better future preparedness. In this review, we systematically analyze and categorize mathematical models that have been developed to design optimal vaccine prioritization strategies of an initially limited vaccine. As older individuals are disproportionately affected by COVID-19, the focus is on models that take age explicitly into account. The lower mobility and activity level of older individuals gives rise to non-trivial trade-offs. Secondary research questions concern the optimal time interval between vaccine doses and spatial vaccine distribution. This review showcases the effect of various modeling assumptions on model outcomes. A solid understanding of these relationships yields better infectious disease models and thus public health decisions during the next pandemic.

Incorporating temporal information during feature engineering bolsters emulation of spatio-temporal emergence

MOTIVATION

Emergent biological dynamics derive from the evolution of lower-level spatial and temporal processes. A long-standing challenge for scientists and engineers is identifying simple low-level rules that give rise to complex higher-level dynamics. High-resolution biological data acquisition enables this identification and has evolved at a rapid pace for both experimental and computational approaches. Simultaneously harnessing the resolution and managing the expense of emerging technologies-e.g. live cell imaging, scRNAseq, agent-based models-requires a deeper understanding of how spatial and temporal axes impact biological systems. Effective emulation is a promising solution to manage the expense of increasingly complex high-resolution computational models. In this research, we focus on the emulation of a tumor microenvironment agent-based model to examine the relationship between spatial and temporal environment features, and emergent tumor properties.

RESULTS

Despite significant feature engineering, we find limited predictive capacity of tumor properties from initial system representations. However, incorporating temporal information derived from intermediate simulation states dramatically improves the predictive performance of machine learning models. We train a deep-learning emulator on intermediate simulation states and observe promising enhancements over emulators trained solely on initial conditions. Our results underscore the importance of incorporating temporal information in the evaluation of spatio-temporal emergent behavior. Nevertheless, the emulators exhibit inconsistent performance, suggesting that the underlying model characterizes unique cell populations dynamics that are not easily replaced.

AVAILABILITY AND IMPLEMENTATION

All source codes for the agent-based model, emulation, and analyses are publicly available at the corresponding DOIs: 10.5281/zenodo.10622155, 10.5281/zenodo.10611675, 10.5281/zenodo.10621244, respectively.

The multi-dimensional challenges of controlling respiratory virus transmission in indoor spaces: Insights from the linkage of a microscopic pedestrian simulation and SARS-CoV-2 transmission model

SARS-CoV-2 transmission in indoor spaces, where most infection events occur, depends on the types and duration of human interactions, among others. Understanding how these human behaviours interface with virus characteristics to drive pathogen transmission and dictate the outcomes of non-pharmaceutical interventions is important for the informed and safe use of indoor spaces. To better understand these complex interactions, we developed the Pedestrian Dynamics-Virus Spread model (PeDViS), an individual-based model that combines pedestrian behaviour models with virus spread models incorporating direct and indirect transmission routes. We explored the relationships between virus exposure and the duration, distance, respiratory behaviour, and environment in which interactions between infected and uninfected individuals took place and compared this to benchmark 'at risk' interactions (1.5 metres for 15 minutes). When considering aerosol transmission, individuals adhering to distancing measures may be at risk due to the buildup of airborne virus in the environment when infected individuals spend prolonged time indoors. In our restaurant case, guests seated at tables near infected individuals were at limited risk of infection but could, particularly in poorly ventilated places, experience risks that surpass that of benchmark interactions. Combining interventions that target different transmission routes can aid in accumulating impact, for instance by combining ventilation with face masks. The impact of such combined interventions depends on the relative importance of transmission routes, which is hard to disentangle and highly context dependent. This uncertainty should be considered when assessing transmission risks upon different types of human interactions in indoor spaces. We illustrated the multi-dimensionality of indoor SARS-CoV-2 transmission that emerges from the interplay of human behaviour and the spread of respiratory viruses. A modelling strategy that incorporates this in risk assessments can help inform policy makers and citizens on the safe use of indoor spaces with varying inter-human interactions.

COVSIM: A stochastic agent-based COVID-19 SIMulation model for North Carolina

We document the evolution and use of the stochastic agent-based COVID-19 simulation model (COVSIM) to study the impact of population behaviors and public health policy on disease spread within age, race/ethnicity, and urbanicity subpopulations in North Carolina. We detail the methodologies used to model the complexities of COVID-19, including multiple agent attributes (i.e., age, race/ethnicity, high-risk medical status), census tract-level interaction network, disease state network, agent behavior (i.e., masking, pharmaceutical intervention (PI) uptake, quarantine, mobility), and variants. We describe its uses outside of the COVID-19 Scenario Modeling Hub (CSMH), which has focused on the interplay of nonpharmaceutical and pharmaceutical interventions, equitability of vaccine distribution, and supporting local county decision-makers in North Carolina. This work has led to multiple publications and meetings with a variety of local stakeholders. When COVSIM joined the CSMH in January 2022, we found it was a sustainable way to support new COVID-19 challenges and allowed the group to focus on broader scientific questions. The CSMH has informed adaptions to our modeling approach, including redesigning our high-performance computing implementation.

Seasonal variations in social contact patterns in a rural population in north India: Implications for pandemic control

Social contact mixing patterns are critical to model the transmission of communicable diseases, and have been employed to model disease outbreaks including COVID-19. Nonetheless, there is a paucity of studies on contact mixing in low and middle-income countries such as India. Furthermore, mathematical models of disease outbreaks do not account for the temporal nature of social contacts. We conducted a longitudinal study of social contacts in rural north India across three seasons and analysed the temporal differences in contact patterns. A contact diary survey was performed across three seasons from October 2015-16, in which participants were queried on the number, duration, and characteristics of contacts that occurred on the previous day. A total of 8,421 responses from 3,052 respondents (49% females) recorded characteristics of 180,073 contacts. Respondents reported a significantly higher number and duration of contacts in the winter, followed by the summer and the monsoon season (Nemenyi post-hoc, p<0.001). Participants aged 0-9 years and 10-19 years of age reported the highest median number of contacts (16 (IQR 12-21), 17 (IQR 13-24) respectively) and were found to have the highest node centrality in the social network of the region (pageranks = 0.20, 0.17). A large proportion (>80%) of contacts that were reported in schools or on public transport involved physical contact. To the best of our knowledge, our study is the first from India to show that contact mixing patterns vary by the time of the year and provides useful implications for pandemic control. We compared the differences in the number, duration and location of contacts by age-group and gender, and studied the impact of the season, age-group, employment and day of the week on the number and duration of contacts using multivariate negative binomial regression. We created a social network to further understand the age and gender-specific contact patterns, and used the contact matrices in each season to parameterise a nine-compartment agent-based model for simulating a COVID-19 epidemic in each season. Our results can be used to parameterize more accurate mathematical models for prediction of epidemiological trends of infections in rural India.

Stochastic transmission in epidemiological models

Recent empirical evidence suggests that the transmission coefficient in susceptible-exposed-infected-removed-like (SEIR-like) models evolves with time, presenting random patterns, and some stylized facts, such as mean-reversion and jumps. To address such observations we propose the use of jump-diffusion stochastic processes to parameterize the transmission coefficient in an SEIR-like model that accounts for death and time-dependent parameters. We provide a detailed theoretical analysis of the proposed model proving the existence and uniqueness of solutions as well as studying its asymptotic behavior. We also compare the proposed model with some variations possibly including jumps. The forecast performance of the considered models, using reported COVID-19 infections from New York City, is then tested in different scenarios. Despite the simplicity of the epidemiological model, by considering stochastic transmission, the forecasted scenarios were fairly accurate.

Generating synthetic population for simulating the spatiotemporal dynamics of epidemics

Agent-based models have gained traction in exploring the intricate processes governing the spread of infectious diseases, particularly due to their proficiency in capturing nonlinear interaction dynamics. The fidelity of agent-based models in replicating real-world epidemic scenarios hinges on the accurate portrayal of both population-wide and individual-level interactions. In situations where comprehensive population data are lacking, synthetic populations serve as a vital input to agent-based models, approximating real-world demographic structures. While some current population synthesizers consider the structural relationships among agents from the same household, there remains room for refinement in this domain, which could potentially introduce biases in subsequent disease transmission simulations. In response, this study unveils a novel methodology for generating synthetic populations tailored for infectious disease transmission simulations. By integrating insights from microsample-derived household structures, we employ a heuristic combinatorial optimizer to recalibrate these structures, subsequently yielding synthetic populations that faithfully represent agent structural relationships. Implementing this technique, we successfully generated a spatially-explicit synthetic population encompassing over 17 million agents for Shenzhen, China. The findings affirm the method's efficacy in delineating the inherent statistical structural relationship patterns, aligning well with demographic benchmarks at both city and subzone tiers. Moreover, when assessed against a stochastic agent-based Susceptible-Exposed-Infectious-Recovered model, our results pinpointed that variations in population synthesizers can notably alter epidemic projections, influencing both the peak incidence rate and its onset.

SARS-CoV-2 Delta and Omicron community transmission networks as added value to contact tracing

OBJECTIVES

Calculations of SARS-CoV-2 transmission networks at a population level have been limited. Networks that estimate infections between individuals and whether this results in a mutation, can be a way to evaluate fitness of a mutational clone by how much it expands in number as well as determining the likelihood a transmission results in a new variant.

METHODS

Australian Delta and Omicron SARS-CoV-2 sequences were downloaded from GISAID. Transmission networks of infection between individuals were estimated using a novel mathematical method.

RESULTS

Many of the sequences were identical, with clone sizes following power law distributions driven by negative binomial probability distributions for both the number of infections per individual and the number of mutations per transmission (median 0.74 nucleotide changes for Delta and 0.71 for Omicron). Using these distributions, an agent-based model was able to replicate the observed clonal network structure, providing a basis for more detailed COVID-19 modelling. Possible recombination events, tracked by insertion/deletion (indel) patterns, were identified for each variant in these outbreaks.

CONCLUSIONS

This modelling approach reveals key transmission characteristics of SARS-CoV-2 and may complement traditional contact tracing. This methodology can also be applied to other diseases as genetic sequencing of viruses becomes more commonplace.

Predictive models for health outcomes due to SARS-CoV-2, including the effect of vaccination: a systematic review

BACKGROUND

The interaction between modelers and policymakers is becoming more common due to the increase in computing speed seen in recent decades. The recent pandemic caused by the SARS-CoV-2 virus was no exception. Thus, this study aims to identify and assess epidemiological mathematical models of SARS-CoV-2 applied to real-world data, including immunization for coronavirus 2019 (COVID-19).

METHODOLOGY

PubMed, JSTOR, medRxiv, LILACS, EconLit, and other databases were searched for studies employing epidemiological mathematical models of SARS-CoV-2 applied to real-world data. We summarized the information qualitatively, and each article included was assessed for bias risk using the Joanna Briggs Institute (JBI) and PROBAST checklist tool. The PROSPERO registration number is CRD42022344542.

FINDINGS

In total, 5646 articles were retrieved, of which 411 were included. Most of the information was published in 2021. The countries with the highest number of studies were the United States, Canada, China, and the United Kingdom; no studies were found in low-income countries. The SEIR model (susceptible, exposed, infectious, and recovered) was the most frequently used approach, followed by agent-based modeling. Moreover, the most commonly used software were R, Matlab, and Python, with the most recurring health outcomes being death and recovery. According to the JBI assessment, 61.4% of articles were considered to have a low risk of bias.

INTERPRETATION

The utilization of mathematical models increased following the onset of the SARS-CoV-2 pandemic. Stakeholders have begun to incorporate these analytical tools more extensively into public policy, enabling the construction of various scenarios for public health. This contribution adds value to informed decision-making. Therefore, understanding their advancements, strengths, and limitations is essential.

Priority populations' experiences of isolation, quarantine and distancing for COVID-19: protocol for a longitudinal cohort study (Optimise Study)

INTRODUCTION

Longitudinal studies can provide timely and accurate information to evaluate and inform COVID-19 control and mitigation strategies and future pandemic preparedness. The Optimise Study is a multidisciplinary research platform established in the Australian state of Victoria in September 2020 to collect epidemiological, social, psychological and behavioural data from priority populations. It aims to understand changing public attitudes, behaviours and experiences of COVID-19 and inform epidemic modelling and support responsive government policy.

METHODS AND ANALYSIS

This protocol paper describes the data collection procedures for the Optimise Study, an ongoing longitudinal cohort of ~1000 Victorian adults and their social networks. Participants are recruited using snowball sampling with a set of seeds and two waves of snowball recruitment. Seeds are purposively selected from priority groups, including recent COVID-19 cases and close contacts and people at heightened risk of infection and/or adverse outcomes of COVID-19 infection and/or public health measures. Participants complete a schedule of monthly quantitative surveys and daily diaries for up to 24 months, plus additional surveys annually for up to 48 months. Cohort participants are recruited for qualitative interviews at key time points to enable in-depth exploration of people's lived experiences. Separately, community representatives are invited to participate in community engagement groups, which review and interpret research findings to inform policy and practice recommendations.

ETHICS AND DISSEMINATION

The Optimise longitudinal cohort and qualitative interviews are approved by the Alfred Hospital Human Research Ethics Committee (# 333/20). The Optimise Study CEG is approved by the La Trobe University Human Ethics Committee (# HEC20532). All participants provide informed verbal consent to enter the cohort, with additional consent provided prior to any of the sub studies. Study findings will be disseminated through public website (https://optimisecovid.com.au/study-findings/) and through peer-reviewed publications.

TRIAL REGISTRATION NUMBER

NCT05323799.

Differentiating behavioral impact with or without vaccination certification under mass vaccination and non-pharmaceutical interventions on mitigating COVID-19

As COVID-19 vaccines became widely available worldwide, many countries implemented vaccination certification, also known as a "green pass", to promote and expedite vaccination on containing virus spread from the latter half of 2021. This policy allowed those vaccinated to have more freedom in public activities compared to more constraints on the unvaccinated in addition to existing non-pharmaceutical interventions (NPIs). Accordingly, the vaccination certification also induced heterogeneous behaviors of unvaccinated and vaccinated groups. This makes it essential yet challenging to model the behavioral impact of vaccination certification on the two groups and the transmission dynamics of COVID-19 within and between the groups. Very limited quantitative work is available for addressing these purposes. Here we propose an extended epidemiological model SEIQRD[Formula: see text] to effectively distinguish the behavioral impact of vaccination certification on unvaccinated and vaccinated groups through incorporating two contrastive transmission chains. SEIQRD[Formula: see text] also quantifies the impact of the green pass policy. With the resurgence of COVID-19 in Greece, Austria, and Israel in 2021, our simulation results indicate that their implementation of vaccination certification brought about more than a 14-fold decrease in the total number of infections and deaths as compared to a scenario with no such a policy. Additionally, a green pass policy may offer a reasonable practical solution to strike the balance between public health and individual's freedom during the pandemic.

Modeling the Impact of COVID-19 Mitigation Strategies in Pennsylvania, USA

Purpose. To estimate the impact on mortality of nonpharmaceutical interventions (NPIs) implemented early in the COVID-19 pandemic. Methods. We implemented an agent-based modified SEIR model of COVID-19, calibrated to match death numbers reported in Pennsylvania from January 2020 to April 2021 and including representations of NPIs implemented in Pennsylvania. To investigate the impact of these strategies, we ran the calibrated model with no interventions and with varying combinations, timings, and levels of interventions. Results. The model closely replicated death outcomes data for Pennsylvania. Without NPIs, deaths in the early months of the pandemic were estimated to be much higher (67,718 deaths compared to actual 6,969). Voluntary interventions alone were relatively ineffective at decreasing mortality. Delaying implementation of interventions led to higher deaths (∼9,000 more deaths with just a 1-week delay). School closure was insufficient as a single intervention but was an important part of a combined intervention strategy. Conclusions. NPIs were effective at reducing deaths early in the COVID-19 pandemic. Agent-based models can incorporate substantial detail on infectious disease spread and the impact of mitigations. Policy Implications. The model supports the importance and effectiveness of NPIs to decrease morbidity from respiratory pathogens. This is particularly important for emerging pathogens for which no vaccines or treatments exist, but such strategies are applicable to a variety of respiratory pathogens.

Highlights

Nonpharmaceutical interventions were used extensively during the early period of the COVID-19 pandemic, but their use has remained controversial.Agent-based modeling of the impact of these mitigation strategies early in the COVID-19 pandemic supports the effectiveness of nonpharmaceutical interventions at decreasing mortality.Since such interventions are not specific to a particular pathogen, they can be used to protect against any respiratory pathogen, known or emerging. They can be applied rapidly when conditions warrant.

Accuracy of Inferences About the Reproductive Number and Superspreading Potential of SARS-CoV-2 with Incomplete Contact Tracing Data

The basic reproductive number (R0) and superspreading potential ( k ) are key epidemiological parameters that inform our understanding of a disease's transmission. Often these values are estimated using the data obtained from contact tracing studies. Here we performed a simulation study to understand how incomplete data due to preferential contact tracing impacted the accuracy and inferences about the transmission of SARS-CoV-2. Our results indicate that as the number of positive contacts traced decreases, our estimates of R0 tend to decrease and our estimates of ktend to increase. Notably, when there are large amounts of positive contacts missed in the tracing process, we can conclude that there is no indication of superspreading even if we know there is. The results of this study highlight the need for a unified public health response to transmissible diseases.

Susceptible-Infected-Susceptible type COVID-19 spread with collective effects

Many models developed to forecast and attempt to understand the COVID-19 pandemic are highly complex, and few take collective behavior into account. As the pandemic progressed individual recurrent infection was observed and simpler susceptible-infected type models were introduced. However, these do not include mechanisms to model collective behavior. Here, we introduce an extension of the SIS model that accounts for collective behavior and show that it has four equilibria. Two of the equilibria are the standard SIS model equilibria, a third is always unstable, and a fourth where collective behavior and infection prevalence interact to produce either node-like or oscillatory dynamics. We then parameterized the model using estimates of the transmission and recovery rates for COVID-19 and present phase diagrams for fixed recovery rate and free transmission rate, and both rates fixed. We observe that regions of oscillatory dynamics exist in both cases and that the collective behavior parameter regulates their extent. Finally, we show that the system exhibits hysteresis when the collective behavior parameter varies over time. This model provides a minimal framework for explaining oscillatory phenomena such as recurring waves of infection and hysteresis effects observed in COVID-19, and other SIS-type epidemics, in terms of collective behavior.

Predicting COVID-19 pandemic waves including vaccination data with deep learning

Introduction

During the recent COVID-19 pandemics, many models were developed to predict the number of new infections. After almost a year, models had also the challenge to include information about the waning effect of vaccines and by infection, and also how this effect start to disappear.

Methods

We present a deep learning-based approach to predict the number of daily COVID-19 cases in 30 countries, considering the non-pharmaceutical interventions (NPIs) applied in those countries and including vaccination data of the most used vaccines.

Results

We empirically validate the proposed approach for 4 months between January and April 2021, once vaccination was available and applied to the population and the COVID-19 variants were closer to the one considered for developing the vaccines. With the predictions of new cases, we can prescribe NPIs plans that present the best trade-off between the expected number of COVID-19 cases and the social and economic cost of applying such interventions.

Discussion

Whereas, mathematical models which include the effect of vaccines in the spread of the SARS-COV-2 pandemic are available, to the best of our knowledge we are the first to propose a data driven method based on recurrent neural networks that considers the waning effect of the immunization acquired either by vaccine administration or by recovering from the illness. This work contributes with an accurate, scalable, data-driven approach to modeling the pandemic curves of cases when vaccination data is available.

Estimating the potential impact and diagnostic requirements for SARS-CoV-2 test-and-treat programs

Oral antivirals have the potential to reduce the public health burden of COVID-19. However, now that we have exited the emergency-phase of the COVID-19 pandemic, declining SARS-CoV-2 clinical testing rates (average testing rates = [Formula: see text]10 tests/100,000 people/day in low-and-middle income countries; <100 tests/100,000 people/day in high-income countries; September 2023) make the development of effective test-and-treat programs challenging. We used an agent-based model to investigate how testing rates and strategies affect the use and effectiveness of oral antiviral test-to-treat programs in four country archetypes of different income levels and demographies. We find that in the post-emergency-phase of the pandemic, in countries where low testing rates are driven by limited testing capacity, significant population-level impact of test-and-treat programs can only be achieved by both increasing testing rates and prioritizing individuals with greater risk of severe disease. However, for all countries, significant reductions in severe cases with antivirals are only possible if testing rates were substantially increased with high willingness of people to seek testing. Comparing the potential population-level reductions in severe disease outcomes of test-to-treat programs and vaccination shows that test-and-treat strategies are likely substantially more resource intensive requiring very high levels of testing (≫100 tests/100,000 people/day) and antiviral use suggesting that vaccination should be a higher priority.

The application of simulation methods during the COVID-19 pandemic: A scoping review

With the outbreak of COVID-19 pandemic, simulation modelling approaches have become effective tools to simulate the potential effects of different intervention measures and predict the dynamic COVID-19 trends. In this scoping review, Studies published between February 2020 and May 2022 that investigated the spread of COVID-19 using four common simulation modeling methods were systematically reported and summarized. Publication trend, characteristics, software, and code availability of included articles were analyzed. Among the included 340 studies, most articles used agent-based model (ABM; n = 258; 75.9 %), followed by the models of system dynamics (n = 42; 12.4 %), discrete event simulation (n = 25; 7.4 %), and hybrid simulation (n = 15; 4.4 %). Furthermore, our review emphasized the purposes and sample time period of included articles. We classified the purpose of the 340 included studies into five categories, most studies mainly analyzed the spread of COVID-19 under policy interventions. For the sample time period analysis, most included studies analyzed the COVID-19 spread in the second wave. Our findings play a crucial role for policymakers to make evidence-based decisions in preventing the spread of COVID-19 pandemic and help in providing scientific decision-makings resilient to similar events and infectious diseases in the future.

The Epidemiological and Economic Impact of COVID-19 in Kazakhstan: An Agent-Based Modeling

BACKGROUND

Our study aimed to assess how effective the preventative measures taken by the state authorities during the pandemic were in terms of public health protection and the rational use of material and human resources.

MATERIALS AND METHODS

We utilized a stochastic agent-based model for COVID-19's spread combined with the WHO-recommended COVID-ESFT version 2.0 tool for material and labor cost estimation.

RESULTS

Our long-term forecasts (up to 50 days) showed satisfactory results with a steady trend in the total cases. However, the short-term forecasts (up to 10 days) were more accurate during periods of relative stability interrupted by sudden outbreaks. The simulations indicated that the infection's spread was highest within families, with most COVID-19 cases occurring in the 26-59 age group. Government interventions resulted in 3.2 times fewer cases in Karaganda than predicted under a "no intervention" scenario, yielding an estimated economic benefit of 40%.

CONCLUSION

The combined tool we propose can accurately forecast the progression of the infection, enabling health organizations to allocate specialists and material resources in a timely manner.

Evaluation of COVID-19 intervention policies in South Korea using the stochastic individual-based model

The COVID-19 pandemic has swept the globe, and countries have responded with various intervention policies to prevent its spread. In this study, we aim to analyze the effectiveness of intervention policies implemented in South Korea. We use a stochastic individual-based model (IBM) with a synthetic population to simulate the spread of COVID-19. Using statistical data, we make the synthetic population and assign sociodemographic attributes to each individual. Individuals go about their daily lives based on their assigned characteristics, and encountering infectors in their daily lives stochastically determines whether they are infected. We reproduce the transmission of COVID-19 using the IBM simulation from November 2020 to February 2021 when three phases of increasingly stringent intervention policies were implemented, and then assess their effectiveness. Additionally, we predict how the spread of infection would have been different if these policies had been implemented in January 2022. This study offers valuable insights into the effectiveness of intervention policies in South Korea, which can assist policymakers and public health officials in their decision-making process.

Individual variation in vaccine immune response can produce bimodal distributions of protection

The ability for vaccines to protect against infectious diseases varies among individuals, but computational models employed to inform policy typically do not account for this variation. Here we examine this issue: we implement a model of vaccine efficacy developed in the context of SARS-CoV-2 in order to evaluate the general implications of modelling correlates of protection on the individual level. Due to high levels of variation in immune response, the distributions of individual-level protection emerging from this model tend to be highly dispersed, and are often bimodal. We describe the specification of the model, provide an intuitive parameterisation, and comment on its general robustness. We show that the model can be viewed as an intermediate between the typical approaches that consider the mode of vaccine action to be either "all-or-nothing" or "leaky". Our view based on this analysis is that individual variation in correlates of protection is an important consideration that may be crucial to designing and implementing models for estimating population-level impacts of vaccination programs.

Framework for converting mechanistic network models to probabilistic models

There are two prominent paradigms for the modelling of networks: in the first, referred to as the mechanistic approach, one specifies a set of domain-specific mechanistic rules that are used to grow or evolve the network over time; in the second, referred to as the probabilistic approach, one describes a model that specifies the likelihood of observing a given network. Mechanistic models (models developed based on the mechanistic approach) are appealing because they capture scientific processes that are believed to be responsible for network generation; however, they do not easily lend themselves to the use of inferential techniques when compared with probabilistic models. We introduce a general framework for converting a mechanistic network model (MNM) to a probabilistic network model (PNM). The proposed framework makes it possible to identify the essential network properties and their joint probability distribution for some MNMs; doing so makes it possible to address questions such as whether two different mechanistic models generate networks with identical distributions of properties, or whether a network property, such as clustering, is over- or under-represented in the networks generated by the model of interest compared with a reference model. The proposed framework is intended to bridge some of the gap that currently exists between the formulation and representation of mechanistic and PNMs. We also highlight limitations of PNMs that need to be addressed in order to close this gap.