Several variants of the SARS-CoV-2 virus have been detected during the COVID-19 pandemic. Some of these new variants have been of health public concern due to their higher infectiousness. We propose a theoretical mathematical model based on differential equations to study the effect of introducing a new, more transmissible SARS-CoV-2 variant in a population. The mathematical model is formulated in such a way that it takes into account the higher transmission rate of the new SARS-CoV-2 strain and the subpopulation of asymptomatic carriers. We find the basic reproduction number using the method of the next generation matrix. This threshold parameter is crucial since it indicates what parameters play an important role in the outcome of the COVID-19 pandemic. We study the local stability of the infection-free and endemic equilibrium states, which are potential outcomes of a pandemic. Moreover, by using a suitable Lyapunov functional and the LaSalle invariant principle, it is proved that if the basic reproduction number is less than unity, the infection-free equilibrium is globally asymptotically stable. Our study shows that the new more transmissible SARS-CoV-2 variant will prevail and the prevalence of the preexistent variant would decrease and eventually disappear. We perform numerical simulations to support the analytic results and to show some effects of a new more transmissible SARS-CoV-2 variant in a population.
Keywords: Lyapunov functions; SARS-CoV-2 virus; basic reproduction number; global stability analysis; variants.