Mutations make pandemics worse or better: modeling SARS-CoV-2 variants and imperfect vaccination

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, $$ , 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 $$ , 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 $$ leads to oscillations. When $$ 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.