1
|
Chen K, Wei F, Zhang X, Jin H, Zhou R, Zuo Y, Fan K. Dynamics of an SVEIR transmission model with protection awareness and two strains. Infect Dis Model 2025; 10:207-228. [PMID: 39469221 PMCID: PMC11513685 DOI: 10.1016/j.idm.2024.10.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2024] [Revised: 09/30/2024] [Accepted: 10/01/2024] [Indexed: 10/30/2024] Open
Abstract
As of May 2024, the main strains of COVID-19 caused hundreds of millions of infection cases and millions of deaths worldwide. In this study, we consider the COVID-19 epidemics with the main strains in the Chinese mainland. We study complex interactions among hosts, non-pharmaceutical interventions, and vaccinations for the main strains by a differential equation model called SVEIR. The disease transmission model incorporates two strains and protection awareness of the susceptible population. Results of this study show that the protection awareness plays a crucial role against infection of the population, and that the vaccines are effective against the circulation of the earlier strains, but ineffective for emerging strains. By using the next generation matrix method, the basic reproduction number of the SVEIR model is firstly obtained. Our analysis by Hurwitz criterion and LaSalle's invariance principle shows that the disease free-equilibrium point is locally and globally asymptotically stable when the threshold value is below one. The existences of endemic equilibrium points are also established, and the global asymptotic stabilities are analyzed using the Lyapunov function method. Further, the SVEIR model is confirmed to satisfy the principle of competitive exclusion, of which the strain with the larger value of the basic reproduction number is dominant. Numerically, the surveillance data with the Omicron strain and the XBB strain are split by the cubic spline interpolation method. The fitting curves against the surveillance data are plotted using the least-squares method from MATLAB. The results indicate that the XBB strain dominates in this study. Moreover, a global sensitivity analysis of the key parameters is performed by using of PRCC. The numerical simulations imply that combination control strategy positively impacts on the infection scale than what separate control strategy does, and that the earlier time producing protection awareness for the public creates less infection scale, further that the increment of protection awareness also reduces the infection scale. Therefore, the policymakers of the local government are suggested to concern the changes of protection awareness of the public.
Collapse
Affiliation(s)
- Kaijing Chen
- School of Mathematics and Statistics, Fuzhou University, Fuzhou, 350116, Fujian, China
| | - Fengying Wei
- School of Mathematics and Statistics, Fuzhou University, Fuzhou, 350116, Fujian, China
- Key Laboratory of Operations Research and Control of Universities in Fujian, Fuzhou University, Fuzhou, 350116, Fujian, China
- Center for Applied Mathematics of Fujian Province, Fuzhou University, Fuzhou, 350116, Fujian, China
| | - Xinyan Zhang
- Jinzhou Center for Disease Control and Prevention, Jinzhou, 121000, Liaoning, China
| | - Hao Jin
- Jinzhou Center for Disease Control and Prevention, Jinzhou, 121000, Liaoning, China
| | - Ruiyang Zhou
- School of Mathematics and Statistics, Fuzhou University, Fuzhou, 350116, Fujian, China
| | - Yue Zuo
- Jinzhou Center for Disease Control and Prevention, Jinzhou, 121000, Liaoning, China
| | - Kai Fan
- Jinzhou Center for Disease Control and Prevention, Jinzhou, 121000, Liaoning, China
| |
Collapse
|
2
|
Saldaña F, Stollenwerk N, Aguiar M. Modelling COVID-19 mutant dynamics: understanding the interplay between viral evolution and disease transmission dynamics. ROYAL SOCIETY OPEN SCIENCE 2024; 11:240919. [PMID: 39493297 PMCID: PMC11529628 DOI: 10.1098/rsos.240919] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/03/2024] [Revised: 08/15/2024] [Accepted: 09/13/2024] [Indexed: 11/05/2024]
Abstract
Understanding virus mutations is critical for shaping public health interventions. These mutations lead to complex multi-strain dynamics often under-represented in models. Aiming to understand the factors influencing variants' fitness and evolution, we explore several scenarios of virus spreading to gain qualitative insight into the factors dictating which variants ultimately predominate at the population level. To this end, we propose a two-strain stochastic model that accounts for asymptomatic transmission, mutations and the possibility of disease import. We find that variants with milder symptoms are likely to spread faster than those with severe symptoms. This is because severe variants can prompt affected individuals to seek medical help earlier, potentially leading to quicker identification and isolation of cases. However, milder or asymptomatic cases may spread more widely, making it harder to control the spread. Therefore, increased transmissibility of milder variants can still result in higher hospitalizations and fatalities due to widespread infection. The proposed model highlights the interplay between viral evolution and transmission dynamics. Offering a nuanced view of factors influencing variant spread, the model provides a foundation for further investigation into mitigating strategies and public health interventions.
Collapse
Affiliation(s)
| | | | - Maíra Aguiar
- Basque Center for Applied Mathematics (BCAM), Bilbao, Spain
- Ikerbasque, Basque Foundation for Science, Bilbao, Spain
| |
Collapse
|
3
|
Bugalia S, Tripathi JP, Wang H. Mutations make pandemics worse or better: modeling SARS-CoV-2 variants and imperfect vaccination. J Math Biol 2024; 88:45. [PMID: 38507066 DOI: 10.1007/s00285-024-02068-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Revised: 07/04/2023] [Accepted: 02/18/2024] [Indexed: 03/22/2024]
Abstract
COVID-19 is a respiratory disease triggered by an RNA virus inclined to mutations. Since December 2020, variants of COVID-19 (especially Delta and Omicron) continuously appeared with different characteristics that influenced death and transmissibility emerged around the world. To address the novel dynamics of the disease, we propose and analyze a dynamical model of two strains, namely native and mutant, transmission dynamics with mutation and imperfect vaccination. It is also assumed that the recuperated individuals from the native strain can be infected with mutant strain through the direct contact with individual or contaminated surfaces or aerosols. We compute the basic reproduction number, R 0 , which is the maximum of the basic reproduction numbers of native and mutant strains. We prove the nonexistence of backward bifurcation using the center manifold theory, and global stability of disease-free equilibrium whenR 0 < 1 , that is, vaccine is effective enough to eliminate the native and mutant strains even if it cannot provide full protection. Hopf bifurcation appears when the endemic equilibrium loses its stability. An intermediate mutation rate ν 1 leads to oscillations. When ν 1 increases over a threshold, the system regains its stability and exhibits an interesting dynamics called endemic bubble. An analytical expression for vaccine-induced herd immunity is derived. The epidemiological implication of the herd immunity threshold is that the disease may effectively be eradicated if the minimum herd immunity threshold is attained in the community. Furthermore, the model is parameterized using the Indian data of the cumulative number of confirmed cases and deaths of COVID-19 from March 1 to September 27 in 2021, using MCMC method. The cumulative cases and deaths can be reduced by increasing the vaccine efficacies to both native and mutant strains. We observe that by considering the vaccine efficacy against native strain as 90%, both cumulative cases and deaths would be reduced by 0.40%. It is concluded that increasing immunity against mutant strain is more influential than the vaccine efficacy against it in controlling the total cases. Our study demonstrates that the COVID-19 pandemic may be worse due to the occurrence of oscillations for certain mutation rates (i.e., outbreaks will occur repeatedly) but better due to stability at a lower infection level with a larger mutation rate. We perform sensitivity analysis using the Latin Hypercube Sampling methodology and partial rank correlation coefficients to illustrate the impact of parameters on the basic reproduction number, the number of cumulative cases and deaths, which ultimately sheds light on disease mitigation.
Collapse
Affiliation(s)
- Sarita Bugalia
- Department of Mathematics, Central University of Rajasthan, Bandar Sindri, Kishangarh, Ajmer, Rajasthan, 305817, India
| | - Jai Prakash Tripathi
- Department of Mathematics, Central University of Rajasthan, Bandar Sindri, Kishangarh, Ajmer, Rajasthan, 305817, India.
| | - Hao Wang
- Department of Mathematical and Statistical Sciences, University of Alberta, Edmonton, AB, T6G 2G1, Canada
| |
Collapse
|
4
|
Rahman A, Kuddus MA, Paul AK, Hasan MZ. The impact of triple doses vaccination and other interventions for controlling the outbreak of COVID-19 cases and mortality in Australia: A modelling study. Heliyon 2024; 10:e25945. [PMID: 38384567 PMCID: PMC10878934 DOI: 10.1016/j.heliyon.2024.e25945] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2023] [Revised: 01/29/2024] [Accepted: 02/05/2024] [Indexed: 02/23/2024] Open
Abstract
COVID-19 is a significant public health problem around the globe, including in Australia. Despite this, Australia's Ministry of Health has expanded COVID-19 control measures widely, logistical trials exist, and the disease burden still needs more clarity. One of the best methods to comprehend the dynamics of disease transmission is by mathematical modeling of COVID-19, which also makes it possible to quantify factors in many places, including Australia. In order to understand the dynamics of COVID-19 in Australia, we examine a mathematical modeling framework for the virus in this study. Australian COVID-19 actual incidence data from January to December 2021 was used to calibrate the model. We also performed a sensitivity analysis of the model parameters and found that the COVID-19 transmission rate was the primary factor in determining the basic reproduction number (R0). Gradually influential intervention policies were established, with accurate effect and coverage regulated with the help of COVID-19 experts in Australia. We simulated data for the period from April 2022 to August 2023. To ascertain which of these outcomes is most effective in lowering the COVID-19 burden, we here assessed the COVID-19 burden (as shown by the number of incident cases and mortality) under a range of intervention scenarios. Regarding the policy of single intervention, the fastest and most efficient way to lower the incidence of COVID-19 is via increasing the first-dose immunization rate, while an improved treatment rate for the afflicted population is also helps to lower mortality in Australia. Furthermore, our results imply that integrating more therapies at the same time increases their efficacy, particularly for mortality, which significantly reduced with a moderate effort, while lowering the number of COVID-19 instances necessitates a major and ongoing commitment.
Collapse
Affiliation(s)
- Azizur Rahman
- School of Computing, Mathematics and Engineering, Charles Sturt University, NSW 2678, Australia
| | - Md Abdul Kuddus
- School of Computing, Mathematics and Engineering, Charles Sturt University, NSW 2678, Australia
- Australian Institute of Tropical Health and Medicine, James Cook University, Townsville, QLD 4810, Australia
- Department of Mathematics, University of Rajshahi, Rajshahi-6205, Bangladesh
| | - Anip Kumar Paul
- Department of Mathematics, University of Rajshahi, Rajshahi-6205, Bangladesh
| | - Md Zobaer Hasan
- School of Computing, Mathematics and Engineering, Charles Sturt University, NSW 2678, Australia
- School of Science, Monash University Malaysia, Jalan Lagoon Selatan, 47500 Bandar Sunway, Selangor D. E., Malaysia
| |
Collapse
|
5
|
Pell B, Brozak S, Phan T, Wu F, Kuang Y. The emergence of a virus variant: dynamics of a competition model with cross-immunity time-delay validated by wastewater surveillance data for COVID-19. J Math Biol 2023; 86:63. [PMID: 36988621 PMCID: PMC10054223 DOI: 10.1007/s00285-023-01900-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2022] [Revised: 12/28/2022] [Accepted: 03/12/2023] [Indexed: 03/30/2023]
Abstract
We consider the dynamics of a virus spreading through a population that produces a mutant strain with the ability to infect individuals that were infected with the established strain. Temporary cross-immunity is included using a time delay, but is found to be a harmless delay. We provide some sufficient conditions that guarantee local and global asymptotic stability of the disease-free equilibrium and the two boundary equilibria when the two strains outcompete one another. It is shown that, due to the immune evasion of the emerging strain, the reproduction number of the emerging strain must be significantly lower than that of the established strain for the local stability of the established-strain-only boundary equilibrium. To analyze the unique coexistence equilibrium we apply a quasi steady-state argument to reduce the full model to a two-dimensional one that exhibits a global asymptotically stable established-strain-only equilibrium or global asymptotically stable coexistence equilibrium. Our results indicate that the basic reproduction numbers of both strains govern the overall dynamics, but in nontrivial ways due to the inclusion of cross-immunity. The model is applied to study the emergence of the SARS-CoV-2 Delta variant in the presence of the Alpha variant using wastewater surveillance data from the Deer Island Treatment Plant in Massachusetts, USA.
Collapse
Affiliation(s)
- Bruce Pell
- Mathematics and Computer Science Department, Lawrence Technological University, 21000 W. 10 Mile Rd, Southfield, MI, 48075, USA.
| | - Samantha Brozak
- School of Mathematical and Statistical Sciences, Arizona State University, 901 S. Palm Walk, Tempe, AZ, 85287-1804, USA
| | - Tin Phan
- Theoretical Biology and Biophysics Group, Houston, Los Alamos, NM, 87545, USA
| | - Fuqing Wu
- Texas Epidemic Public Health Institute, Houston, TX, 77030, USA
- Department of Epidemiology, Human Genetics and Environmental Sciences, School of Public Health, The University of Texas Health Science Center at Houston, Houston, Texas, 77030, USA
| | - Yang Kuang
- School of Mathematical and Statistical Sciences, Arizona State University, 901 S. Palm Walk, Tempe, AZ, 85287-1804, USA
| |
Collapse
|
6
|
Reis RF, Congdon P. Editorial: Epidemiological considerations in COVID-19 forecasting. FRONTIERS IN EPIDEMIOLOGY 2023; 2:1119559. [PMID: 38455284 PMCID: PMC10910939 DOI: 10.3389/fepid.2022.1119559] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/08/2022] [Accepted: 12/20/2022] [Indexed: 03/09/2024]
Affiliation(s)
- Ruy Freitas Reis
- Instituto de Ciências Exatas, Departamento de Ciência da Computação, Universidade Federal de Juiz de Fora, Juiz de Fora, Brazil
- Pós-Graduação em Modelagem Computacional, Universidade Federal de Juiz de Fora, Juiz de Fora, Brazil
| | - Peter Congdon
- Queen Mary University of London, London, United Kingdom
| |
Collapse
|
7
|
Phan T, Brozak S, Pell B, Gitter A, Mena KD, Kuang Y, Wu F. A simple SEIR-V model to estimate COVID-19 prevalence and predict SARS-CoV-2 transmission using wastewater-based surveillance data. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2022:2022.07.17.22277721. [PMID: 35898336 PMCID: PMC9327624 DOI: 10.1101/2022.07.17.22277721] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Abstract
Wastewater-based surveillance (WBS) has been widely used as a public health tool to monitor SARS-CoV-2 transmission. However, epidemiological inference from WBS data remains understudied and limits its application. In this study, we have established a quantitative framework to estimate COVID-19 prevalence and predict SARS-CoV-2 transmission through integrating WBS data into an SEIR-V model. We conceptually divide the individual-level viral shedding course into exposed, infectious, and recovery phases as an analogy to the compartments in population-level SEIR model. We demonstrated that the temperature effect on viral losses in the sewer can be straightforwardly incorporated in our framework. Using WBS data from the second wave of the pandemic (Oct 02, 2020 â€" Jan 25, 2021) in the Great Boston area, we showed that the SEIR-V model successfully recapitulates the temporal dynamics of viral load in wastewater and predicts the true number of cases peaked earlier and higher than the number of reported cases by 16 days and 8.6 folds ( R = 0.93), respectively. This work showcases a simple, yet effective method to bridge WBS and quantitative epidemiological modeling to estimate the prevalence and transmission of SARS-CoV-2 in the sewershed, which could facilitate the application of wastewater surveillance of infectious diseases for epidemiological inference and inform public health actions.
Collapse
Affiliation(s)
- Tin Phan
- Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, New Mexico, USA
| | - Samantha Brozak
- School of Mathematical and Statistical Sciences, Arizona State University, Arizona, USA
| | - Bruce Pell
- Department of Mathematics and Computer Science, Lawrence Technological University, MI, USA
| | - Anna Gitter
- The University of Texas Health Science Center at Houston, School of Public Health, Houston, Texas, USA 77030
| | - Kristina D. Mena
- The University of Texas Health Science Center at Houston, School of Public Health, Houston, Texas, USA 77030
| | - Yang Kuang
- School of Mathematical and Statistical Sciences, Arizona State University, Arizona, USA
| | - Fuqing Wu
- The University of Texas Health Science Center at Houston, School of Public Health, Houston, Texas, USA 77030
| |
Collapse
|
8
|
Mathematical Modeling: Global Stability Analysis of Super Spreading Transmission of Respiratory Syncytial Virus (RSV) Disease. COMPUTATION 2022. [DOI: 10.3390/computation10070120] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
In this paper, a model for the transmission of respiratory syncytial virus (RSV) in a constant human population in which there exist super spreading infected individuals (who infect many people during a single encounter) is considered. It has been observed in the epidemiological data for the diseases caused by this virus that there are cases where some individuals are super-spreaders of the virus. We formulate a simply SEIrIsR (susceptible–exposed–regular infected–super-spreading infected–recovered) mathematical model to describe the dynamics of the transmission of this disease. The proposed model is analyzed using the standard stability method by using Routh-Hurwitz criteria. We obtain the basic reproductive number (R0) using the next generation method. We establish that when R0<1, the disease-free state is locally asymptotically stable and the disease endemic state is unstable. The reverse is true when R0>1, the disease endemic state becomes the locally asymptotically stable state and the disease-free state becomes unstable. It is also established that the two equilibrium states are globally asymptotically stable. The numerical simulations show how the dynamics of the disease change as values of the parameters in the SEIrIsR are varied.
Collapse
|
9
|
Massard M, Eftimie R, Perasso A, Saussereau B. A multi-strain epidemic model for COVID-19 with infected and asymptomatic cases: application to French data. J Theor Biol 2022; 545:111117. [PMID: 35513167 PMCID: PMC9059428 DOI: 10.1016/j.jtbi.2022.111117] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Revised: 03/22/2022] [Accepted: 03/29/2022] [Indexed: 10/27/2022]
Abstract
Many SARS-CoV-2 variants have appeared over the last months, and many more will continue to appear. Understanding the competition between these different variants could help make future predictions on the evolution of epidemics. In this study we use a mathematical model to investigate the impact of three different SARS-CoV-2 variants on the spread of COVID-19 across France, between January-May 2021 (before vaccination was extended to the full population). To this end, we use the data from Geodes (produced by Public Health France) and a particle swarm optimisation algorithm, to estimate the model parameters and further calculate a value for the basic reproduction number R0. Sensitivity and uncertainty analysis is then used to better understand the impact of estimated parameter values on the number of infections leading to both symptomatic and asymptomatic individuals. The results confirmed that, as expected, the alpha, beta and gamma variants are more transmissible than the original viral strain. In addition, the sensitivity results showed that the beta/gamma variants could have lead to a larger number of infections in France (of both symptomatic and asymptomatic people).
Collapse
Affiliation(s)
- Mathilde Massard
- Laboratoire de Mathématiques de Besançon, University of Bourgogne Franche-Comté, UFR ST 16 route de GRAY, 25030 Besançon, France.
| | - Raluca Eftimie
- Laboratoire de Mathématiques de Besançon, University of Bourgogne Franche-Comté, UFR ST 16 route de GRAY, 25030 Besançon, France.
| | - Antoine Perasso
- Laboratoire Chrono-environnement, University of Bourgogne Franche-Comté, UFR ST 16 route de GRAY, 25030 Besançon, France.
| | - Bruno Saussereau
- Laboratoire de Mathématiques de Besançon, University of Bourgogne Franche-Comté, UFR ST 16 route de GRAY, 25030 Besançon, France.
| |
Collapse
|
10
|
Estimating the Number of COVID-19 Cases and Impact of New COVID-19 Variants and Vaccination on the Population in Kerman, Iran: A Mathematical Modeling Study. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2022; 2022:6624471. [PMID: 35495892 PMCID: PMC9039779 DOI: 10.1155/2022/6624471] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/01/2021] [Revised: 11/24/2021] [Accepted: 02/09/2022] [Indexed: 12/23/2022]
Abstract
COVID-19 is spreading all over Iran, and Kerman is one of the most affected cities. We conducted this study to predict COVID-19-related deaths, hospitalization, and infected cases under different scenarios (scenarios A, B, and C) by 31 December 2021 in Kerman. We also aimed to assess the impact of new COVID-19 variants and vaccination on the total number of COVID-19 cases, deaths, and hospitalizations (scenarios D, E, and F) using the modified susceptible-exposed-infected-removed (SEIR) model. We calibrated the model using deaths reported from the start of the epidemic to August 30, 2021. A Monte Carlo Markov Chain (MCMC) uncertainty analysis was used to estimate 95% uncertainty intervals (UI). We also calculated the time-varying reproductive number (Rt) following time-dependent methods. Under the worst-case scenario (scenario A; contact rate = 10, self‐isolation rate = 30%, and average vaccination shots per day = 5,000), the total number of infections by December 31, 2021, would be 1,625,000 (95% UI: 1,112,000–1,898,000) with 6,700 deaths (95% UI: 5,200–8,700). With the presence of alpha and delta variants without vaccine (scenario D), the total number of infected cases and the death toll were estimated to be 957,000 (95% UI: 208,000–1,463,000) and 4,500 (95% UI: 1,500–7,000), respectively. If 70% of the population were vaccinated when the alpha variant was dominant (scenario E), the total number of infected cases and deaths would be 608,000 (95% UI: 122,000–743,000) and 2,700 (95% UI: 700–4,000), respectively. The Rt was ≥1 almost every day during the epidemic. Our results suggest that policymakers should concentrate on improving vaccination and interventions, such as reducing social contacts, stricter limitations for gathering, public education to promote social distancing, incensing case finding and contact tracing, effective isolation, and quarantine to prevent more COVID-19 cases, hospitalizations, and deaths in Kerman.
Collapse
|
11
|
de León UAP, Avila-Vales E, Huang KL. Modeling COVID-19 dynamic using a two-strain model with vaccination. CHAOS, SOLITONS, AND FRACTALS 2022; 157:111927. [PMID: 35185299 PMCID: PMC8847090 DOI: 10.1016/j.chaos.2022.111927] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/27/2021] [Revised: 02/05/2022] [Accepted: 02/15/2022] [Indexed: 05/21/2023]
Abstract
Multiple strains of the SARS-CoV-2 have arisen and jointly influence the trajectory of the coronavirus disease (COVID-19) pandemic. However, current models rarely account for this multi-strain dynamics and their different transmission rate and response to vaccines. We propose a new mathematical model that accounts for two virus variants and the deployment of a vaccination program. To demonstrate utility, we applied the model to determine the control reproduction number ( R c ) and the per day infection, death and recovery rates of each strain in the US pandemic. The model dynamics predicted the rise of the alpha variant and shed light on potential impact of the delta variant in 2021. We obtained the minimum percentage of fully vaccinated individuals to reduce the spread of the variants in combination with other intervention strategies to deaccelerate the rise of a multi-strain pandemic.
Collapse
Affiliation(s)
| | - Eric Avila-Vales
- Facultad de Matemáticas, Universidad Autónoma de Yucatán. Anillo Periférico Norte. Tablaje Catastral 13615, C.P. 97119. Mérida, Yucatán
| | - Kuan-Lin Huang
- Department of Genetics and Genomic Sciences. Center for Transformative Disease Modeling. Tisch Cancer Institute. Icahn Institute for Data Science and Genomic Technology. Icahn School of Medicine at Mount Sinai. New York. NY 10029. USA
| |
Collapse
|
12
|
Mathematical Modeling to Study Optimal Allocation of Vaccines against COVID-19 Using an Age-Structured Population. AXIOMS 2022. [DOI: 10.3390/axioms11030109] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Vaccination against the coronavirus disease 2019 (COVID-19) started in early December of 2020 in the USA. The efficacy of the vaccines vary depending on the SARS-CoV-2 variant. Some countries have been able to deploy strong vaccination programs, and large proportions of their populations have been fully vaccinated. In other countries, low proportions of their populations have been vaccinated, due to different factors. For instance, countries such as Afghanistan, Cameroon, Ghana, Haiti and Syria have less than 10% of their populations fully vaccinated at this time. Implementing an optimal vaccination program is a very complex process due to a variety of variables that affect the programs. Besides, science, policy and ethics are all involved in the determination of the main objectives of the vaccination program. We present two nonlinear mathematical models that allow us to gain insight into the optimal vaccination strategy under different situations, taking into account the case fatality rate and age-structure of the population. We study scenarios with different availabilities and efficacies of the vaccines. The results of this study show that for most scenarios, the optimal allocation of vaccines is to first give the doses to people in the 55+ age group. However, in some situations the optimal strategy is to first allocate vaccines to the 15–54 age group. This situation occurs whenever the SARS-CoV-2 transmission rate is relatively high and the people in the 55+ age group have a transmission rate 50% or less that of those in the 15–54 age group. This study and similar ones can provide scientific recommendations for countries where the proportion of vaccinated individuals is relatively small or for future pandemics.
Collapse
|
13
|
Impact of Infective Immigrants on COVID-19 Dynamics. MATHEMATICAL AND COMPUTATIONAL APPLICATIONS 2022; 27. [DOI: 10.3390/mca27010011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The COVID-19 epidemic is an unprecedented and major social and economic challenge worldwide due to the various restrictions. Inflow of infective immigrants have not been given prominence in several mathematical and epidemiological models. To investigate the impact of imported infection on the number of deaths, cumulative infected and cumulative asymptomatic, we formulate a mathematical model with infective immigrants and considering vaccination. The basic reproduction number of the special case of the model without immigration of infective people is derived. We varied two key factors that affect the transmission of COVID-19, namely the immigration and vaccination rates. In addition, we considered two different SARS-CoV-2 transmissibilities in order to account for new more contagious variants such as Omicron. Numerical simulations using initial conditions approximating the situation in the US when the vaccination program was starting show that increasing the vaccination rate significantly improves the outcomes regarding the number of deaths, cumulative infected and cumulative asymptomatic. Other factors are the natural recovery rates of infected and asymptomatic individuals, the waning rate of the vaccine and the vaccination rate. When the immigration rate is increased significantly, the number of deaths, cumulative infected and cumulative asymptomatic increase. Consequently, accounting for the level of inflow of infective immigrants may help health policy/decision-makers to formulate policies for public health prevention programs, especially with respect to the implementation of the stringent preventive lock down measure.
Collapse
|
14
|
Modeling the Transmission of the SARS-CoV-2 Delta Variant in a Partially Vaccinated Population. Viruses 2022; 14:v14010158. [PMID: 35062363 PMCID: PMC8781299 DOI: 10.3390/v14010158] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Revised: 12/16/2021] [Accepted: 01/12/2022] [Indexed: 12/19/2022] Open
Abstract
In a population with ongoing vaccination, the trajectory of a pandemic is determined by how the virus spreads in unvaccinated and vaccinated individuals that exhibit distinct transmission dynamics based on different levels of natural and vaccine-induced immunity. We developed a mathematical model that considers both subpopulations and immunity parameters, including vaccination rates, vaccine effectiveness, and a gradual loss of protection. The model forecasted the spread of the SARS-CoV-2 delta variant in the US under varied transmission and vaccination rates. We further obtained the control reproduction number and conducted sensitivity analyses to determine how each parameter may affect virus transmission. Although our model has several limitations, the number of infected individuals was shown to be a magnitude greater (~10×) in the unvaccinated subpopulation compared to the vaccinated subpopulation. Our results show that a combination of strengthening vaccine-induced immunity and preventative behavioral measures like face mask-wearing and contact tracing will likely be required to deaccelerate the spread of infectious SARS-CoV-2 variants.
Collapse
|
15
|
Abstract
INTRODUCTION Appearances of SARS-CoV-2 variants have created havoc and additional challenges for the ongoing vaccination drive against pandemic COVID-19. Interestingly, several vaccine platforms are showing great potential to produce successful vaccines against SARS-CoV-2 and its variants. Billions of COVID-19 vaccine doses have been administered worldwide. Mix-and-Match COVID-19 vaccines involving the mixing of the same platform vaccines and also two different vaccine platforms may provide greater protection against SARS-CoV-2 and its variants. COVID-19 vaccines have become one of the most important tools to mitigate the ongoing pandemic COVID-19. AREAS COVERED We describe SARS-Cov-2 variants, their impact on the population, COVID-19 vaccines, diverse vaccine platforms, doses of vaccines, the efficacy of vaccines against SARS-CoV-2 and its variants, mitigation of the COVID-19 transmission- alternatives to vaccines. EXPERT OPINION Diverse vaccine platforms may safeguard against ongoing, deadly SARS-CoV-2 and its infectious variants. The efficacies of COVID-19 vaccines are significantly high after the administration of the second dose. Further, it protects individuals including vulnerable patients with co-morbidities from SARS-CoV-2 and its variants. The hospitalizations and deaths of the individuals may be prevented by COVID-19 vaccines.
Collapse
Affiliation(s)
- Bhaswati Chatterjee
- Chemical Science, National Institute of Pharmaceutical Education and Research, Hyderabad, India,CONTACT Bhaswati Chatterjee Chemical Science National Institute of Pharmaceutical Education and Research, India
| | - Suman S. Thakur
- Proteomics and Cell Signaling, Centre for Cellular and Molecular Biology, Uppal Road, Hyderabad, India,Suman S. Thakur Principal Scientist, Proteomics and Cell Signaling, Lab W110 Centre for Cellular and Molecular Biology Uppal Road, Hyderabad-500007, India
| |
Collapse
|
16
|
Hohan R, Milu P, Paraschiv S, Casangiu C, Tudor A, Vlaicu O, Banica L, Surleac M, Florea D, Otelea D. The Predictive Value of Mutation Screening for Anticipating COVID-19 Waves. Pathogens 2021; 10:pathogens10111464. [PMID: 34832619 PMCID: PMC8622441 DOI: 10.3390/pathogens10111464] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2021] [Revised: 11/05/2021] [Accepted: 11/09/2021] [Indexed: 11/18/2022] Open
Abstract
Emerging SARS-CoV-2 strains continue to generate difficulties for authorities and health care professionals worldwide due to enhanced transmissibility and/or immune response evasion. The appearance of the Alpha and Delta strains has been associated with substantial increases in the number of COVID-19 cases and associated deaths. Whole Genome Sequencing (WGS) continues to be the gold standard for molecular surveillance of the pandemics but other assays such as mutation genotyping can be used to reduce costs and allocated time. This study investigates the efficiency of mutation screening tests compared to WGS and their predictive value to anticipate future waves. A very high degree of fidelity for this type of assay was found, regardless of the method used. The positive predictive value (PPV) of 4/5 markers was over 95% for the detection of Alpha and Delta variants. By estimating the prevalence of the Alpha and Delta strains using genotyping assays and fitting the data to a mathematical model, a five week period between the point of exponential growth of variant prevalence and a drastic increase in case numbers was found. For that reason, raising awareness about the efficacy of mutation screening could help authorities adopt better measures in the future.
Collapse
Affiliation(s)
- Robert Hohan
- “Prof. Dr. Matei Bals” National Institute for Infectious Diseases, 021105 Bucharest, Romania; (R.H.); (P.M.); (S.P.); (C.C.); (O.V.); (L.B.); (M.S.); (D.F.)
| | - Petre Milu
- “Prof. Dr. Matei Bals” National Institute for Infectious Diseases, 021105 Bucharest, Romania; (R.H.); (P.M.); (S.P.); (C.C.); (O.V.); (L.B.); (M.S.); (D.F.)
- Faculty of Medicine, “Carol Davila” University of Medicine and Pharmacy, 050474 Bucharest, Romania
| | - Simona Paraschiv
- “Prof. Dr. Matei Bals” National Institute for Infectious Diseases, 021105 Bucharest, Romania; (R.H.); (P.M.); (S.P.); (C.C.); (O.V.); (L.B.); (M.S.); (D.F.)
- Faculty of Medicine, “Carol Davila” University of Medicine and Pharmacy, 050474 Bucharest, Romania
| | - Corina Casangiu
- “Prof. Dr. Matei Bals” National Institute for Infectious Diseases, 021105 Bucharest, Romania; (R.H.); (P.M.); (S.P.); (C.C.); (O.V.); (L.B.); (M.S.); (D.F.)
- Faculty of Biology, University of Bucharest, 077120 Bucharest, Romania
| | - Andreea Tudor
- “Prof. Dr. C.C. Iliescu” Emergency Institute for Cardiovascular Diseases, 021105 Bucharest, Romania;
| | - Ovidiu Vlaicu
- “Prof. Dr. Matei Bals” National Institute for Infectious Diseases, 021105 Bucharest, Romania; (R.H.); (P.M.); (S.P.); (C.C.); (O.V.); (L.B.); (M.S.); (D.F.)
| | - Leontina Banica
- “Prof. Dr. Matei Bals” National Institute for Infectious Diseases, 021105 Bucharest, Romania; (R.H.); (P.M.); (S.P.); (C.C.); (O.V.); (L.B.); (M.S.); (D.F.)
| | - Marius Surleac
- “Prof. Dr. Matei Bals” National Institute for Infectious Diseases, 021105 Bucharest, Romania; (R.H.); (P.M.); (S.P.); (C.C.); (O.V.); (L.B.); (M.S.); (D.F.)
- Research Institute of the University of Bucharest (ICUB), University of Bucharest, 060031 Bucharest, Romania
| | - Dragos Florea
- “Prof. Dr. Matei Bals” National Institute for Infectious Diseases, 021105 Bucharest, Romania; (R.H.); (P.M.); (S.P.); (C.C.); (O.V.); (L.B.); (M.S.); (D.F.)
- Faculty of Medicine, “Carol Davila” University of Medicine and Pharmacy, 050474 Bucharest, Romania
| | - Dan Otelea
- “Prof. Dr. Matei Bals” National Institute for Infectious Diseases, 021105 Bucharest, Romania; (R.H.); (P.M.); (S.P.); (C.C.); (O.V.); (L.B.); (M.S.); (D.F.)
- Correspondence:
| |
Collapse
|
17
|
Mancuso M, Eikenberry SE, Gumel AB. Will vaccine-derived protective immunity curtail COVID-19 variants in the US? Infect Dis Model 2021; 6:1110-1134. [PMID: 34518808 PMCID: PMC8426325 DOI: 10.1016/j.idm.2021.08.008] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2021] [Revised: 08/12/2021] [Accepted: 08/31/2021] [Indexed: 12/14/2022] Open
Abstract
Multiple effective vaccines are currently being deployed to combat the COVID-19 pandemic, and are viewed as the major factor in marked reductions of disease burden in regions with moderate to high vaccination coverage. The effectiveness of COVID-19 vaccination programs is, however, significantly threatened by the emergence of new SARS-COV-2 variants that, in addition to being more transmissible than the wild-type (original) strain, may at least partially evade existing vaccines. A two-strain (one wild-type, one variant) and two-group (vaccinated or otherwise) mechanistic mathematical model is designed and used to assess the impact of the vaccine-induced cross-protective efficacy on the spread the COVID-19 pandemic in the United States. Rigorous analysis of the model shows that, in the absence of any co-circulating SARS-CoV-2 variant, the vaccine-derived herd immunity threshold needed to eliminate the wild-type strain can be achieved if 59% of the US population is fully-vaccinated with either the Pfizer or Moderna vaccine. This threshold increases to 76% if the wild-type strain is co-circulating with the Alpha variant (a SARS-CoV-2 variant that is 56% more transmissible than the wild-type strain). If the wild-type strain is co-circulating with the Delta variant (which is estimated to be 100% more transmissible than the wild-type strain), up to 82% of the US population needs to be vaccinated with either of the aforementioned vaccines to achieve the vaccine-derived herd immunity. Global sensitivity analysis of the model reveal the following four parameters as the most influential in driving the value of the reproduction number of the variant strain (hence, COVID-19 dynamics) in the US: (a) the infectiousness of the co-circulating SARS-CoV-2 variant, (b) the proportion of individuals fully vaccinated (using Pfizer or Moderna vaccine) against the wild-type strain, (c) the cross-protective efficacy the vaccines offer against the variant strain and (d) the modification parameter accounting for the reduced infectiousness of fully-vaccinated individuals experiencing breakthrough infection. Specifically, numerical simulations of the model show that future waves or surges of the COVID-19 pandemic can be prevented in the US if the two vaccines offer moderate level of cross-protection against the variant (at least 67%). This study further suggests that a new SARS-CoV-2 variant can cause a significant disease surge in the US if (i) the vaccine coverage against the wild-type strain is low (roughly <66%) (ii) the variant is much more transmissible (e.g., 100% more transmissible), than the wild-type strain, or (iii) the level of cross-protection offered by the vaccine is relatively low (e.g., less than 50%). A new SARS-CoV-2 variant will not cause such surge in the US if it is only moderately more transmissible (e.g., the Alpha variant, which is 56% more transmissible) than the wild-type strain, at least 66% of the population of the US is fully vaccinated, and the three vaccines being deployed in the US (Pfizer, Moderna, and Johnson & Johnson) offer a moderate level of cross-protection against the variant.
Collapse
Affiliation(s)
- Marina Mancuso
- Arizona State University, School of Mathematical and Statistical Sciences, Tempe, AZ, 85287, USA
| | - Steffen E. Eikenberry
- Arizona State University, School of Mathematical and Statistical Sciences, Tempe, AZ, 85287, USA
| | - Abba B. Gumel
- Arizona State University, School of Mathematical and Statistical Sciences, Tempe, AZ, 85287, USA
- Department of Mathematics and Applied Mathematics, University of Pretoria, Pretoria, 0002, South Africa
| |
Collapse
|