  The crisis caused by the COVID-19 outbreak around the globe raised an
increasing concern about the ongoing emergence of variants of SARS-CoV-2 that
may evade the immune response provided by vaccines. New variants appear due to
mutation, and as the cases accumulate, the probability of the emergence of a
variant of concern increases. In this article, we propose a modified SIR model
with waning immunity that captures the competition of two strain classes of an
infectious disease under the effect of vaccination with a highly contagious and
deadly strain class emerging from a prior strain due to mutation. When these
strains compete for a limited supply of susceptible individuals, changes in the
efficiency of vaccines may affect the behaviour of the disease in a non-trivial
way, resulting in complex outcomes. We characterise the parameter space
including intrinsic parameters of the disease, and using the vaccine
efficiencies as control variables. We find different types of transcritical
bifurcations between endemic fixed points and a disease-free equilibrium and
identify a region of strain competition where the two strain classes coexist
during a transient period. We show that a strain can be extinguished either due
to strain competition or vaccination, and we obtain the critical values of the
efficiency of vaccines to eradicate the disease. Numerical studies using
parameters estimated from publicly reported data agree with our theoretical
results. Our mathematical model could be a tool to assess quantitatively the
vaccination policies of competing and emerging strains using the dynamics in
epidemics of infectious diseases.
