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 when \(R_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 \(\nu _1\) leads to oscillations. When \(\nu _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.
Similar content being viewed by others
References
Ahmed SE (2008). Markov chain monte carlo: Stochastic simulation for bayesian inference. https://doi.org/10.1198/tech.2008.s542
Alexander ME, Bowman C, Moghadas SM, Summers R, Gumel AB, Sahai BM (2004) A vaccination model for transmission dynamics of influenza. SIAM J Appl Dyn Syst 3(4):503–524. https://doi.org/10.1137/030600370
Anderson RM (1992) The concept of herd immunity and the design of community-based immunization programmes. Vaccine 10:928e935. https://doi.org/10.1016/0264-410X(92)90327-G
Anderson RM, May RM (1985) Vaccination and herd immunity to infectious diseases. Nature 318:323e329. https://doi.org/10.1038/318323a0
Arino J, McCluskey CC, van den Driessche P (2003) Global results for an epidemic model with vaccination that exhibits backward bifurcation. SIAM J Appl Math 64(1):260–276. https://doi.org/10.1137/S0036139902413829
Arruda EF, Pastore DH, Dias CM, Das SS (2021) Modelling and optimal control of multi strain epidemics, with application to COVID-19. https://doi.org/10.1371/journal.pone.0257512; arXiv preprint arXiv:2101:08137
Bajiya VP, Bugalia S, Tripathi JP (2020) Mathematical modeling of COVID-19: impact of non-pharmaceutical interventions in India. Chaos: Interdiscip J Nonlinear Sci 30(11):113143. https://doi.org/10.1063/5.0021353
Bonhoeffer S, Nowak MA (1994) Mutation and the evolution of virulence. Proc R Soc Lond Ser B: Biol Sci 258(1352):133–140. https://doi.org/10.1098/rspb.1994.0153
Bugalia S, Bajiya VP, Tripathi JP, Li MT, Sun GQ (2020) Mathematical modeling of COVID-19 transmission: the roles of intervention strategies and lockdown. Math Biosci Eng 17(5):5961–5986. https://doi.org/10.3934/mbe.2020318
Bugalia S, Tripathi JP, Wang H (2021) Mathematical modeling of intervention and low medical resource availability with delays: Applications to COVID-19 outbreaks in Spain and Italy. Math Biosci Eng 18(5):5865–5920. https://doi.org/10.3934/mbe.2021295
Bugalia S, Tripathi JP, Wang H (2023) Estimating the time-dependent effective reproduction number and vaccination rate for COVID-19 in the USA and India. Math Biosci Eng 20(3):4673–4689. https://doi.org/10.3934/mbe.2023216
Bussiness Insider (2021) https://tinyurl.com/6q3g47b8. Accessed 2 Nov 2021
Butler G, Freedman HI, Waltman P (1986) Uniformly persistent systems. Proc Am Math Soc 96(3):425–430. https://doi.org/10.2307/2046588
Cable News Network (2020) Fauci says covid-19 vaccine may not get us to herd immunity if too many people refuse to get it, CNN. https://edition.cnn.com/2020/06/28/health/fauci-coronavirus-vaccine-contact-tracing-aspen/index.html. Accessed 10 Sept 2021
Cai L, Xiang J, Li X, Lashari AA (2012) A two-strain epidemic model with mutant strain and vaccination. J Appl Math Comput 40(1):125–142. https://doi.org/10.1007/s12190-012-0580-x
Cai LM, Li Z, Song X (2018) Global analysis of an epidemic model with vaccination. J Appl Math Comput 57(1):605–628. https://doi.org/10.1007/s12190-017-1124-1
Castillo-Chavez C, Blower S, Van den Driessche P, Kirschner D, Yakubu AA, (eds) (2002) Mathematical approaches for emerging and reemerging infectious diseases: models, methods, and theory. Springer, Berlin. https://doi.org/10.1007/978-1-4613-0065-6
Castillo-Chavez C, Feng Z, Huang W (2002) On the computation of \(R_0\) and its role in global stability. IMA Vol Math Appl 125:229–250
Castillo-Chavez C, Song B (2004) Dynamical models of tuberculosis and their applications. Math Biosci Eng 1(2):361. https://doi.org/10.3934/mbe.2004.1.361
Centers for Disease Control and Prevention (CDC) (2021a). https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-classifications.html. Accessed 13 Oct 2021
Centers for Disease Control and Prevention (CDC) (2021b). https://www.cdc.gov/coronavirus/2019-ncov/your-health/reinfection.html. Accessed 13 Oct 2021
Chowdhury MM, Islam MR, Hossain MS, Tabassum N, Peace A (2022) Incorporating the mutational landscape of SARS-COV-2 variants and case-dependent vaccination rates into epidemic models. Infect Disease Model 7(2):75–82. https://doi.org/10.1016/j.idm.2022.02.003
Deng X, Garcia-Knight MA, Khalid MM, Servellita V, Wang C, Morris MK, Sotomayor-González A, Glasner DR, Reyes KR, Gliwa AS, Reddy NP (2021) Transmission, infectivity, and neutralization of a spike L452R SARS-CoV-2 variant. Cell. https://doi.org/10.1016/j.cell.2021.04.025
Diekmann O, Heesterbeek JA, Roberts MG (2010) The construction of next-generation matrices for compartmental epidemic models. J R Soc Interface 7(47):873–885. https://doi.org/10.1098/rsif.2009.0386
Duong BV, Larpruenrudee P, Fang T, Hossain SI, Saha SC, Gu Y, Islam MS (2022) Is the SARS CoV-2 omicron variant deadlier and more transmissible than delta variant? Int J Environ Res Public Health 19(8):4586. https://doi.org/10.3390/ijerph19084586
Elbasha EH, Gumel AB (2021) Vaccination and herd immunity thresholds in heterogeneous populations. J Math Biol. https://doi.org/10.1007/s00285-021-01686-z
Eletreby R, Zhuang Y, Carley KM, Yağan O, Poor HV (2020) The effects of evolutionary adaptations on spreading processes in complex networks. Proc Natl Acad Sci 117(11):5664–5670. https://doi.org/10.1073/pnas.1918529117
Eron JJ, Vernazza PL, Johnston DM, Seillier-Moiseiwitsch F, Alcorn TM, Fiscus SA, Cohen MS (1998) Resistance of HIV-1 to antiretroviral agents in blood and seminal plasma: implications for transmission. AIDS 12(15):F181–F189. https://doi.org/10.1097/00002030-199815000-00003
Freedman HI, Waltman P (1985) Persistence in a model of three competitive populations. Math Biosci 73(1):89–101. https://doi.org/10.1016/0025-5564(85)90078-1
Friedman TL (2020) Is Sweden doing it right? New York Times. https://www.nytimes.com/2020/04/28/opinion/coronavirus-sweden.html. Accessed 3 Oct 2021
Fudolig M, Howard R (2020) The local stability of a modified multi-strain SIR model for emerging viral strains. PLoS ONE 15(12):e0243408. https://doi.org/10.1371/journal.pone.0243408
Gonzalez-Parra G, Martínez-Rodríguez D, Villanueva-Micó RJ (2021) Impact of a new SARS-CoV-2 variant on the population: A mathematical modeling approach. Math Comput Appl 26(2):25. https://doi.org/10.3390/mca26020025
Gukenheimer J, Holmes P (1983) Nonlinear oscillations, Dynamical Systems, and Bifurcation of Vector Fields. Springer, New York. https://doi.org/10.1007/978-1-4612-1140-2
Gumel AB, McCluskey CC, Watmough J (2006) An SVEIR model for assessing potential impact of an imperfect anti-SARS vaccine. Math Biosci Eng 3(3):485. https://doi.org/10.3934/mbe.2006.3.485
Gupta S, Ferguson NM, Anderson RM (1997) Vaccination and the population structure of antigenically diverse pathogens that exchange genetic material. Proc R Soc Lond B 264(1387):1435–1443. https://doi.org/10.1098/rspb.1997.0200
Haario H, Laine M, Mira A, Saksman E (2006) DRAM: efficient adaptive MCMC. Stat Comput 16(4):339–354. https://doi.org/10.1007/s11222-006-9438-0
Haario H, Saksman E, Tamminen J (2001) An adaptive Metropolis algorithm. Bernoulli 7(2):223–242. https://doi.org/10.2307/3318737
Health, The Sciences (2021) https://science.thewire.in/the-sciences/covid-19-reinfection-hong-kong-man-immunity/. Accessed 10 Nov 2021
Hethcote HW (2000) The mathematics of infectious diseases. SIAM Rev Soc Ind Appl Math 42(4):599–653. https://doi.org/10.1137/S0036144500371907
How Dangerous Is the Delta Variant (B.1.617.2)? (2023) American society of microbiology. https://asm.org/Articles/2021/July/How-Dangerous-is-the-Delta-Variant-B-1-617-2. Accessed 20 June 2023
Iboi EA, Ngonghala CN, Gumel AB (2020) Will an imperfect vaccine curtail the COVID-19 pandemic in the US? Infect Disease Model 5:510–524. https://doi.org/10.1016/j.idm.2020.07.006
Iboi EA, Sharomi OO, Ngonghala CN, Gumel AB (2020) Mathematical modeling and analysis of COVID-19 pandemic in Nigeria. Math Biosci Eng 17(6):7192–7220. https://doi.org/10.3934/mbe.2020369
Iwami S, Takeuchi Y, Liu X (2007) Avian-human influenza epidemic model. Math Biosci 207(1):1–25. https://doi.org/10.1016/j.mbs.2006.08.001
Johns Hopkins Medicine (2021) https://www.hopkinsmedicine.org/health/conditions-and-diseases/coronavirus/a-new-strain-of-coronavirus-what-youshould-know. Accessed 17 Sept 2021
Kermack WO, McKendrick AG (1927) A contribution to the mathematical theory of epidemics. Proc R Soc Lond Ser A Contain Pap Math Phys Charact 115(772):700–721. https://doi.org/10.1098/rspa.1927.0118
Khyar O, Allali K (2020) Global dynamics of a multi-strain SEIR epidemic model with general incidence rates: application to COVID-19 pandemic. Nonlinear Dyn 102(1):489–509. https://doi.org/10.1007/s11071-020-05929-4
Korber B, Fischer WM, Gnanakaran S, Yoon H, Theiler J, Abfalterer W, Hengartner N, Giorgi EE, Bhattacharya T, Foley B, Hastie KM (2020) Tracking changes in SARS-CoV-2 Spike: evidence that D614G increases infectivity of the COVID-19 virus. Cell 182(4):812–827. https://doi.org/10.1016/j.cell.2020.06.043
LaSalle JP (1976) The stability of dynamical systems. Soc Ind Appl Math. https://doi.org/10.1137/1.9781611970432
Lauer SA, Grantz KH, Bi Q, Jones FK, Zheng Q, Meredith HR, Azman AS, Reich NG, Lessler J (2020) The incubation period of coronavirus disease 2019 (COVID-19) from publicly reported confirmed cases: estimation and application. Ann Intern Med 172(9):577–582. https://doi.org/10.7326/M20-0504
Lemieux JE, Li JZ (2021) Uncovering ways that emerging SARS-CoV-2 lineages may increase transmissibility. J Infect Dis 223(10):1663–1665. https://doi.org/10.1093/infdis/jiab083
Li J, Zhou Y, Ma Z, Hyman JM (2004) Epidemiological models for mutating pathogens. SIAM J Appl Math 65(1):1–23. https://doi.org/10.1137/S0036139903430185
Liu L, Ren X, Liu X (2018) Dynamical behaviors of an influenza epidemic model with virus mutation. J Biol Syst 26(03):455–472. https://doi.org/10.1142/S0218339018500201
Liu M, Liz E, Rost G (2015) Endemic bubbles generated by delayed behavioral response: global stability and bifurcation switches in an SIS model. SIAM J Appl Math 75(1):75–91. https://doi.org/10.1137/140972652
Liu WM (1994) Criterion of Hopf bifurcations without using eigenvalues. J Math Anal Appl 182(1):250–256. https://doi.org/10.1006/jmaa.1994.1079
Liu X, Takeuchi Y, Iwami S (2008) SVIR epidemic models with vaccination strategies. J Theor Biol 253(1):1–11. https://doi.org/10.1016/j.jtbi.2007.10.014
MacIntyre CR, Costantino V, Trent MJ (2021) Modelling of COVID-19 vaccination strategies and herd immunity, in scenarios of limited and full vaccine supply in NSW. Vaccine, Australia. https://doi.org/10.1016/j.vaccine.2021.04.042
Marino S, Hogue IB, Ray CJ, Kirschner DE (2008) A methodology for performing global uncertainty and sensitivity analysis in systems biology. J Theor Biol 254(1):178–196. https://doi.org/10.1016/j.jtbi.2008.04.011
Martcheva M (2015) An introduction to mathematical epidemiology. Springer, New York. https://doi.org/10.1007/978-1-4899-7612-3
Martcheva M, Iannelli M, Li XZ (2007) Subthreshold coexistence of strains: the impact of vaccination and mutation. Math Biosci Eng 4(2):287. https://doi.org/10.3934/mbe.2007.4.287
May RM, Nowak MA (1995) Coinfection and the evolution of parasite virulence. Proc R Soc Lond B 261(1361):209–215. https://doi.org/10.1098/rspb.1995.0138
May RM, Nowak MA (1994) Superinfection, metapopulation dynamics, and the evolution of diversity. J Theor Biol 170(1):95–114. https://doi.org/10.1006/jtbi.1994.1171
McLean AR (1995) Vaccination, evolution and changes in the efficacy of vaccines: a theoretical framework. Proc R Soc Lond B 261(1362):389–393. https://doi.org/10.1098/rspb.1995.0164
Ministry of Health and Family Welfare, Government of India (2021). https://www.mohfw.gov.in/covid_vaccination/vaccination/faqs.html#what-to-expect-before-vaccination-2. Accessed 16 Oct 2021
Nature news (2021) https://www.nature.com/articles/d41586-021-01059-y. Accessed 18 Oct 2021
Olliaro P, Torreele E, Vaillant M (2021) COVID-19 vaccine efficacy and effectiveness-the elephant (not) in the room. Lancet Microbe. https://doi.org/10.1016/S2666-5247(21)00069-0
Palese P, Young JF (1982) Variation of influenza A, B, and C viruses. Science 215(4539):1468–1474. https://doi.org/10.1126/science.7038875
Parton R, Hall E, Wardlaw AC (1994) Responses to Bordetella pertussis mutant strains and to vaccination in the coughing rat model of pertussis. J Med Microbiol 40(5):307–312. https://doi.org/10.1099/00222615-40-5-307
Paul S, Mahata A, Mukherjee S, Roy B, Salimi M, Ahmadian A (2022) Study of fractional order SEIR epidemic model and effect of vaccination on the spread of COVID-19. Int J Appl Comput Math 8(5):237. https://doi.org/10.1007/s40819-022-01411-4
Porco TC, Blower SM (1998) Designing HIV vaccination policies: subtypes and cross-immunity. Interfaces 28(3):167–190. https://doi.org/10.1287/inte.28.3.167
Porco TC, Blower SM (2000) HIV vaccines: the effect of the mode of action on the coexistence of HIV subtypes. Math Popul Stud 8(2):205–229. https://doi.org/10.1080/08898480009525481
Public Health England (2021) investigation-of-novel-sars-cov-2-variant-variant-of-concern-20201201. https://www.gov.uk/government/publications/investigation-of-novel-sars-cov-2-variant-variant-of-concern-20201201. Accessed 27 Sept 2021
Rahimi F, Abadi AT (2021) Implications of the emergence of a new variant of SARS-CoV-2, VUI-202012/01. Arch Med Res 52(5):569–571. https://doi.org/10.1016/j.arcmed.2021.01.001
Sansonetti PJ, Arondel J (1989) Construction and evaluation of a double mutant of Shigella flexneri as a candidate for oral vaccination against shigellosis. Vaccine 7(5):443–450. https://doi.org/10.1016/0264-410X(89)90160-6
Sato S, Suzuki K, Akahane Y, Akamatsu K, Akiyama K, Yunomura K, Tsuda F, Tanaka T, Okamoto H, Miyakawa Y, Mayumi M (1995) Hepatitis B virus strains with mutations in the core promoter in patients with fulminant hepatitis. Ann Intern Med 122(4):241–248. https://doi.org/10.7326/0003-4819-122-4-199502150-00001
Scherer A, McLean A (2002) Mathematical models of vaccination. Br Med Bull 62(1):187–199. https://doi.org/10.1093/bmb/62.1.187
Shastri J, Parikh S, Aggarwal V, Agrawal S, Chatterjee N, Shah R, Devi P, Mehta P, Pandey R (2021) Severe SARS-CoV-2 breakthrough reinfection with Delta variant after recovery from breakthrough infection by Alpha variant in a fully vaccinated health worker. Front Med 1:1379. https://doi.org/10.3389/fmed.2021.737007
Van den Driessche P, Watmough J (2002) Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci 180(1–2):29–48. https://doi.org/10.1016/S0025-5564(02)00108-6
Martcheva M (2014) Avian flu: modeling and implications for control. J Biol Syst 22(01):151–75. https://doi.org/10.1142/S0218339014500090
World Health Organization (2021a) https://covid19.who.int/region/searo/country/in. Accessed 29 Oct 2021
World Health Organization (2021b) https://covid19.who.int/. Accessed 29 Oct 2021
World Health Organization (2021c) WHO advisory committee on variola virus research: report of the thirteenth meeting report (World Health Organization, 2011). https://apps.who.int/iris/handle/10665/70778. Accessed 29 Oct 2021
World Health Organization (2021d) Tracking SARS-CoV-2 Variants. https://www.who.int/en/activities/tracking-SARS-CoV-2-variants/. Accessed 29 Oct 2021
Xue L, Jing SL, Wang H (2022) Evaluating the impacts of non-pharmaceutical interventions on the transmission dynamics of COVID-19 in Canada based on mobile network. PLoS ONE 16(12):e0261424. https://doi.org/10.1371/journal.pone.0261424
Yagan O, Sridhar A, Eletreby R, Levin S, Plotkin JB, Poor HV (2021) Modeling and analysis of the spread of COVID-19 under a multiple-strain model with mutations. Harvard Data Sci Rev 4:1. https://doi.org/10.1162/99608f92.a11bf693
Acknowledgements
The research work of Sarita Bugalia is supported by the Council of Scientific & Industrial Research (CSIR), India [File No. 09/1131(0025)/2018-EMR-I]. The research work of Jai Prakash Tripathi is supported by the Science and Engineering Research Board (SERB), India [File No. ECR/2017/002786] and [File No. MTR/2022/001028]. We are highly thankful to all the anonymous reviewers and editor for their insightful comments and suggestions, which helped us to improve the manuscript considerably.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix A
Appendix A
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Bugalia, S., Tripathi, J.P. & Wang, H. Mutations make pandemics worse or better: modeling SARS-CoV-2 variants and imperfect vaccination. J. Math. Biol. 88, 45 (2024). https://doi.org/10.1007/s00285-024-02068-x
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s00285-024-02068-x
Keywords
- COVID-19
- Endemic bubble
- Mutation
- Imperfect vaccination
- MCMC
- Two strain dynamics
- Hopf bifurcation
- Transcritical bifurcation
- Sensitivity analysis