Coupled, multi-strain epidemic models of mutating pathogens

Stochastic Analysis for the Dual Virus Parallel Transmission Model with Immunity Delay

In this article, the qualitative properties of a stochastic dual virus parallel transmission model with immunity delay are analyzed. First, we use Lyapunov theory to study the existence and uniqueness of the global positive solution of the proposed model. Second, the threshold values of the persistence and extinction of two viruses were obtained. Finally, the numerical simulation verifies the theoretical results. The results show that the immunity delay and the intensity of noise have important effects on the two diseases spreading in parallel.

Characterization of differential susceptibility and differential infectivity epidemic models

Heterogeneity in susceptibility and infectivity is a central issue in epidemiology. Although the latter has received some attention recently, the former is often neglected in modeling of epidemic systems. Moreover, very few studies consider both of these heterogeneities. This paper is concerned with the characterization of epidemic models with differential susceptibility and differential infectivity under a general setup. Specifically, we investigate the global asymptotic behavior of equilibria of these systems when the network configuration of the Susceptible-Infectious interactions is strongly connected. These results prove two conjectures by Bonzi et al. (J Math Biol 62:39-64, 2011) and Hyman and Li (Math Biosci Eng 3:89-100, 2006). Moreover, we consider the scenario in which the strong connectivity hypothesis is dropped. In this case, the model exhibits a wider range of dynamical behavior, including the rise of boundary and interior equilibria, all based on the topology of network connectivity. Finally, a model with multidirectional transitions between infectious classes is presented and completely analyzed.

Vaccine breakthrough and rebound infections modeling: Analysis for the United States and the ten U.S. HHS regions

A vaccine breakthrough infection and a rebound infection cases of COVID-19 are studied and analyzed for the ten U.S. Department of Health and Human Services (HHS) regions and the United States as a nation in this work. An innovative multi-strain susceptible-vaccinated-exposed-asymptomatic-symptomatic-recovered (SVEAIR) epidemic model is developed for this purpose for a population assumed to be susceptible to n-different variants of the disease, and those who are vaccinated and recovered from a specific strain k(k ≤ n) of the disease are immune to present strain and its predecessors j = 1, 2, …, k, but can still be infected by newer emerging strains j = k + 1, k + 2, …, n. The model is used to estimate epidemiological parameters, namely, the latent and infectious periods, the transmission rates, vaccination rates, recovery rates for each of the Delta B.1.617.2, Omicron B.1.1.529, and lineages BA.2, BA.2.12.1, BA.4, BA.5, BA.1.1, BA.4.6, and BA.5.2.6 for the United States and for each of the ten HHS regions. The transmission rate is estimated for both the asymptomatic and symptomatic cases. The effect of vaccines on each strain is analyzed. Condition that guarantees existence of an endemic with certain number of strains is derived and used to describe the endemic state of the population.

Stochastic disease spreading and containment policies under state-dependent probabilities

We analyze the role of disease containment policy in the form of treatment in a stochastic economic-epidemiological framework in which the probability of the occurrence of random shocks is state-dependent, namely it is related to the level of disease prevalence. Random shocks are associated with the diffusion of a new strain of the disease which affects both the number of infectives and the growth rate of infection, and the probability of such shocks realization may be either increasing or decreasing in the number of infectives. We determine the optimal policy and the steady state of such a stochastic framework, which is characterized by an invariant measure supported on strictly positive prevalence levels, suggesting that complete eradication is never a possible long run outcome where instead endemicity will prevail. Our results show that: (i) independently of the features of the state-dependent probabilities, treatment allows to shift leftward the support of the invariant measure; and (ii) the features of the state-dependent probabilities affect the shape and spread of the distribution of disease prevalence over its support, allowing for a steady state outcome characterized by a distribution alternatively highly concentrated over low prevalence levels or more spread out over a larger range of prevalence (possibly higher) levels.

Analysis of multi-strain infection of vaccinated and recovered population through epidemic model: Application to COVID-19

In this work, an innovative multi-strain SV EAIR epidemic model is developed for the study of the spread of a multi-strain infectious disease in a population infected by mutations of the disease. The population is assumed to be completely susceptible to n different variants of the disease, and those who are vaccinated and recovered from a specific strain k (k ≤ n) are immune to previous and present strains j = 1, 2, ⋯, k, but can still be infected by newer emerging strains j = k + 1, k + 2, ⋯, n. The model is designed to simulate the emergence and dissemination of viral strains. All the equilibrium points of the system are calculated and the conditions for existence and global stability of these points are investigated and used to answer the question as to whether it is possible for the population to have an endemic with more than one strain. An interesting result that shows that a strain with a reproduction number greater than one can still die out on the long run if a newer emerging strain has a greater reproduction number is verified numerically. The effect of vaccines on the population is also analyzed and a bound for the herd immunity threshold is calculated. The validity of the work done is verified through numerical simulations by applying the proposed model and strategy to analyze the multi-strains of the COVID-19 virus, in particular, the Delta and the Omicron variants, in the United State.

The Global dynamics of a SIR model considering competitions among multiple strains in patchy environments

Pandemics, such as Covid-19 and AIDS, tend to be highly contagious and have the characteristics of global spread and existence of multiple virus strains. To analyze the competition among different strains, a high dimensional SIR model studying multiple strains' competition in patchy environments is introduced in this work. By introducing the basic reproductive number of different strains, we found global stability conditions of disease-free equilibrium and persistence conditions of the model. The competition exclusion conditions of that model are also given. This work gives some insights into the properties of the multiple strain patchy model and all of the analysis methods used in this work could be used in other related high dimension systems.

The effect of human vaccination behaviour on strain competition in an infectious disease: An imitation dynamic approach

Strain competition plays an important role in shaping the dynamics of multiple pathogen outbreaks in a population. Competition may lead to exclusion of some pathogens, while it may influence the invasion of an emerging mutant in the population. However, little emphasis has been given to understand the influence of human vaccination choice on pathogen competition or strain invasion for vaccine-preventable infectious diseases. Coupling game dynamic framework of vaccination choice and compartmental disease transmission model of two strains, we explore invasion and persistence of a mutant in the population despite having a lower reproduction rate than the resident one. We illustrate that higher perceived strain severity and lower perceived vaccine efficacy are necessary conditions for the persistence of a mutant strain. The numerical simulation also extends these invasion and persistence analyses under asymmetric cross-protective immunity of these strains. We show that the dynamics of this cross-immunity model under human vaccination choices is determined by the interplay of parameters defining the cross-immune response function, perceived risk of infection, and vaccine efficacy, and it can exhibit invasion and persistence of mutant strain, even complete exclusion of resident strain in the regime of sufficiently high perceived risk. We conclude by discussing public health implications of the results, that proper risk communication in public about the severity of the disease is an important task to reduce the chance of mutant invasion. Thus, understanding pathogen competitions under social interactions and choices may be an important component for policymakers for strategic decision-making.

Threshold dynamics and optimal control on an age-structured SIRS epidemic model with vaccination

We consider a vaccination control into a age-structured susceptible-infective-recovered-susceptible (SIRS) model and study the global stability of the endemic equilibrium by the iterative method. The basic reproduction number $ R_0 $ is obtained. It is shown that if $ R_0 < 1 $, then the disease-free equilibrium is globally asymptotically stable, if $ R_0 > 1 $, then the disease-free and endemic equilibrium coexist simultaneously, and the global asymptotic stability of endemic equilibrium is also shown. Additionally, the Hamilton-Jacobi-Bellman (HJB) equation is given by employing the Bellman's principle of optimality. Through proving the existence of viscosity solution for HJB equation, we obtain the optimal vaccination control strategy. Finally, numerical simulations are performed to illustrate the corresponding analytical results.

An epidemic model for an evolving pathogen with strain-dependent immunity

Between pandemics, the influenza virus exhibits periods of incremental evolution via a process known as antigenic drift. This process gives rise to a sequence of strains of the pathogen that are continuously replaced by newer strains, preventing a build up of immunity in the host population. In this paper, a parsimonious epidemic model is defined that attempts to capture the dynamics of evolving strains within a host population. The 'evolving strains' epidemic model has many properties that lie in-between the Susceptible-Infected-Susceptible and the Susceptible-Infected-Removed epidemic models, due to the fact that individuals can only be infected by each strain once, but remain susceptible to reinfection by newly emerged strains. Coupling results are used to identify key properties, such as the time to extinction. A range of reproduction numbers are explored to characterise the model, including a novel quasi-stationary reproduction number that can be used to describe the re-emergence of the pathogen into a population with 'average' levels of strain immunity, analogous to the beginning of the winter peak in influenza. Finally the quasi-stationary distribution of the evolving strains model is explored via simulation.

On the probability of strain invasion in endemic settings: Accounting for individual heterogeneity and control in multi-strain dynamics

•

Endemic infection can insulate host populations from invasion by mutant variants.

•

The timing of control implementation strongly influences its efficacy.

•

Controls that exacerbate host heterogeneity outperform those that curtail it.

•

Differential control can facilitate strain invasion and eventual replacement.

Pathogen evolution is an imminent threat to global health that has warranted, and duly received, considerable attention within the medical, microbiological and modelling communities. Outbreaks of new pathogens are often ignited by the emergence and transmission of mutant variants descended from wild-type strains circulating in the community. In this work we investigate the stochastic dynamics of the emergence of a novel disease strain, introduced into a population in which it must compete with an existing endemic strain. In analogy with past work on single-strain epidemic outbreaks, we apply a branching process approximation to calculate the probability that the new strain becomes established. As expected, a critical determinant of the survival prospects of any invading strain is the magnitude of its reproduction number relative to that of the background endemic strain. Whilst in most circumstances this ratio must exceed unity in order for invasion to be viable, we show that differential control scenarios can lead to less-fit novel strains invading populations hosting a fitter endemic one. This analysis and the accompanying findings will inform our understanding of the mechanisms that have led to past instances of successful strain invasion, and provide valuable lessons for thwarting future drug-resistant strain incursions.