A co-infection model for human papillomavirus (HPV) and syphilis with cost-effectiveness optimal control analysis is developed and presented. The full co-infection model is shown to undergo the phenomenon of backward bifurcation when a certain condition is satisfied. The global asymptotic stability of the disease-free equilibrium of the full model is shown not to exist when the associated reproduction number is less than unity. The existence of endemic equilibrium of the syphilis-only sub-model is shown to exist and the global asymptotic stability of the disease-free and endemic equilibria of the syphilis-only sub-model was established, for a special case. Sensitivity analysis is also carried out on the parameters of the model. Using the syphilis associated reproduction number, [Formula: see text], as the response function, it is observed that the five-ranked parameters that drive the dynamics of the co-infection model are the demographic parameter [Formula: see text], the effective contact rate for syphilis transmission, [Formula: see text], the progression rate to late stage of syphilis [Formula: see text], and syphilis treatment rates: [Formula: see text] and [Formula: see text] for co-infected individuals in compartments [Formula: see text] and [Formula: see text], respectively. Moreover, when the HPV associated reproduction number, [Formula: see text], is used as the response function, the five most dominant parameters that drive the dynamics of the model are the demographic parameter [Formula: see text], the effective contact rate for HPV transmission, [Formula: see text], the fraction of HPV infected who develop persistent HPV [Formula: see text], the fraction of individuals vaccinated against incident HPV infection [Formula: see text] and the HPV vaccine efficacy [Formula: see text]. Numerical simulations of the optimal control model showed that the optimal control strategy which implements syphilis treatment controls for singly infected individuals is the most cost-effective of all the control strategies in reducing the burden of HPV and syphilis co-infections.

In this work, we develop and present a co-infection model for human papillomavirus (HPV) and syphilis with cost-effectiveness optimal control analysis. The full co-infection model is shown to undergo the phenomenon of backward bifurcation when a certain condition is satisfied. The global asymptotic stability of the disease-free equilibrium of the full model is shown \textbf{not to exist}, when the associated reproduction number is less than unity. The existence of endemic equilibrium of the syphilis-only sub-model is shown to exist and the global asymptotic stability of the disease-free and endemic equilibria of both the syphilis-only sub-model and HPV-only sub-model were established. The global asymptotic stability of disease-free equilibrium of the HPV-only sub-model is also proven. Numerical simulations of the optimal control model showed that the optimal control strategy which implements syphilis treatment controls for singly infected individuals is the most cost-effective of all the control strategies in reducing the burden of HPV and syphilis co-infections.