A mathematical model for COVID-19 considering waning immunity, vaccination and control measures

Kinetics of Viral Shedding for Outbreak Surveillance of Emerging Infectious Diseases: Modeling Approach to SARS-CoV-2 Alpha and Omicron Infection

BACKGROUND

Previous studies have highlighted the importance of viral shedding using cycle threshold (Ct) values obtained via reverse transcription polymerase chain reaction to understand the epidemic trajectories of SARS-CoV-2 infections. However, it is rare to elucidate the transition kinetics of Ct values from the asymptomatic or presymptomatic phase to the symptomatic phase before recovery using individual repeated Ct values.

OBJECTIVE

This study proposes a novel Ct-enshrined compartment model to provide a series of quantitative measures for delineating the full trajectories of the dynamics of viral load from infection until recovery.

METHODS

This Ct-enshrined compartment model was constructed by leveraging Ct-classified states within and between presymptomatic and symptomatic compartments before recovery or death among people with infections. A series of recovery indices were developed to assess the net kinetic movement of Ct-up toward and Ct-down off recovery. The model was applied to (1) a small-scale community-acquired Alpha variant outbreak under the "zero-COVID-19" policy without vaccines in May 2021 and (2) a large-scale community-acquired Omicron variant outbreak with high booster vaccination rates following the lifting of the "zero-COVID-19" policy in April 2022 in Taiwan. The model used Bayesian Markov chain Monte Carlo methods with the Metropolis-Hastings algorithm for parameter estimation. Sensitivity analyses were conducted by varying Ct cutoff values to assess the robustness of the model.

RESULTS

The kinetic indicators revealed a marked difference in viral shedding dynamics between the Alpha and Omicron variants. The Alpha variant exhibited slower viral shedding and lower recovery rates, but the Omicron variant demonstrated swifter viral shedding and higher recovery rates. Specifically, the Alpha variant showed gradual Ct-up transitions and moderate recovery rates, yielding a presymptomatic recovery index slightly higher than 1 (1.10), whereas the Omicron variant had remarkable Ct-up transitions and significantly higher asymptomatic recovery rates, resulting in a presymptomatic recovery index much higher than 1 (152.5). Sensitivity analysis confirmed the robustness of the chosen Ct values of 18 and 25 across different recovery phases. Regarding the impact of vaccination, individuals without booster vaccination had a 19% higher presymptomatic incidence rate compared to those with booster vaccination. Breakthrough infections in boosted individuals initially showed similar Ct-up transition rates but higher rates in later stages compared to nonboosted individuals. Overall, booster vaccination improved recovery rates, particularly during the symptomatic phase, although recovery rates for persistent asymptomatic infection were similar regardless of vaccination status once the Ct level exceeded 25.

CONCLUSIONS

The study provides new insights into dynamic Ct transitions, with the notable finding that Ct-up transitions toward recovery outpaced Ct-down and symptom-surfacing transitions during the presymptomatic phase. The Ct-up against Ct-down transition varies with variants and vaccination status. The proposed Ct-enshrined compartment model is useful for the surveillance of emerging infectious diseases in the future to prevent community-acquired outbreaks.

A stochastic approach for co-evolution process of virus and human immune system

Infectious diseases have long been a shaping force in human history, necessitating a comprehensive understanding of their dynamics. This study introduces a co-evolution model that integrates both epidemiological and evolutionary dynamics. Utilizing a system of differential equations, the model represents the interactions among susceptible, infected, and recovered populations for both ancestral and evolved viral strains. Methodologically rigorous, the model's existence and uniqueness have been verified, and it accommodates both deterministic and stochastic cases. A myriad of graphical techniques have been employed to elucidate the model's dynamics. Beyond its theoretical contributions, this model serves as a critical instrument for public health strategy, particularly predicting future outbreaks in scenarios where viral mutations compromise existing interventions.

An Epidemic Model with Infection Age and Vaccination Age Structure

A model of epidemic dynamics is developed that incorporates continuous variables for infection age and vaccination age. The model analyzes pre-symptomatic and symptomatic periods of an infected individual in terms of infection age. This property is shown to be of major importance in the severity of the epidemic, when the infectious period of an infected individual precedes the symptomatic period. The model also analyzes the efficacy of vaccination in terms of vaccination age. The immunity to infection of vaccinated individuals varies with vaccination age and is also of major significance in the severity of the epidemic. Application of the model to the 2003 SARS epidemic in Taiwan and the COVID-19 epidemic in New York provides insights into the dynamics of these diseases. It is shown that the SARS outbreak was effectively contained due to the complete overlap of infectious and symptomatic periods, allowing for the timely isolation of affected individuals. In contrast, the pre-symptomatic spread of COVID-19 in New York led to a rapid, uncontrolled epidemic. These findings underscore the critical importance of the pre-symptomatic infectious period and the vaccination strategies in influencing the dynamics of an epidemic.

Mathematical Modeling Evaluates How Vaccinations Affected the Course of COVID-19 Disease Progression

The regulation policies implemented, the characteristics of vaccines, and the evolution of the virus continue to play a significant role in the progression of the SARS-CoV-2 pandemic. Numerous research articles have proposed using mathematical models to predict the outcomes of different scenarios, with the aim of improving awareness and informing policy-making. In this work, we propose an expansion to the classical SEIR epidemiological model that is designed to fit the complex epidemiological data of COVID-19. The model includes compartments for vaccinated, asymptomatic, hospitalized, and deceased individuals, splitting the population into two branches based on the severity of progression. In order to investigate the impact of the vaccination program on the spread of COVID-19 in Greece, this study takes into account the realistic vaccination program implemented in Greece, which includes various vaccination rates, different dosages, and the administration of booster shots. It also examines for the first time policy scenarios at crucial time-intervention points for Greece. In particular, we explore how alterations in the vaccination rate, immunity loss, and relaxation of measures regarding the vaccinated individuals affect the dynamics of COVID-19 spread. The modeling parameters revealed an alarming increase in the death rate during the dominance of the delta variant and before the initiation of the booster shot program in Greece. The existing probability of vaccinated people becoming infected and transmitting the virus sets them as catalytic players in COVID-19 progression. Overall, the modeling observations showcase how the criticism of different intervention measures, the vaccination program, and the virus evolution has been present throughout the various stages of the pandemic. As long as immunity declines, new variants emerge, and vaccine protection in reducing transmission remains incompetent; monitoring the complex vaccine and virus evolution is critical to respond proactively in the future.