Mathematical analysis of a two-sex Human Papillomavirus (HPV) model

A two-sex deterministic model for Human Papillomavirus (HPV) that assesses the impact of treatment and vaccination on its transmission dynamics is designed and rigorously analyzed. The model is shown to exhibit the phenomenon of backward bifurcation, caused by the imperfect vaccine as well as the re-infection of individuals who recover from a previous infection, when the associated reproduction number is less than unity. Analysis of the reproduction number reveals that the impact of treatment on effective control of the disease is conditional, and depends on the sign of a certain threshold unlike when preventive measures are implemented (i.e. condom use and vaccination of both males and females). Numerical simulations of the model showed that, based on the parameter values used therein, a vaccine (with 75% efficacy) for male population with about 40% condom compliance by females will result in a significant reduction in the disease burden in the population. Also, the numerical simulations of the model reveal that with 70% condom compliance by the male population, administering female vaccine (with 45% efficacy) is sufficient for effective control of the disease.