Dynamics of a two-strain vaccination model for polio

A new deterministic model for the transmission dynamics of two strains of polio, the vaccine-derived polio virus (VDPV) and the wild polio virus (WPV), in a population is designed and rigorously analysed. It is shown that Oral Polio Vaccine (OPV) reversion (leading to increased incidences of WPV and VDPV strains), together with the combined effect of vaccinating a fraction of the unvaccinated susceptible and missed susceptible children, could induce the phenomenon of backward bifurcation when the associated reproduction number of the model is less than unity. Furthermore, the model undergoes competitive exclusion, where the strain with the higher reproduction number (greater than unity) drives the other (with reproduction number less than unity) to extinction. In the absence of OPV reversions (leading to the co-existence of both strains in the population), it is shown that the disease-free equilibrium of the model is globally-asymptotically stable whenever the associated reproduction number is less than unity. Numerical simulations of the model suggest that the model undergoes the phenomenon of competitive exclusion, where the strain with the higher reproduction number (greater than unity) drives the other to extinction. Furthermore, co-existence of the two strains is feasible if their respective reproduction number are equal or approximately equal (and greater than unity).