A two-sex model of human papillomavirus infection: Vaccination strategies and a case study

Vaccination is effective in preventing human papillomavirus (HPV) infection. It still remains debatable whether males should be included in a vaccination program and unclear how to allocate the vaccine in genders to achieve the maximum benefits. In this paper, we use a two-sex model to assess HPV vaccination strategies and use the data from Guangxi Province in China as a case study. Both mathematical analysis and numerical simulations show that the basic reproduction number, an important indicator of the transmission potential of the infection, achieves its minimum when the priority of vaccination is given to the gender with a smaller recruit rate. Given a fixed amount of vaccine, splitting the vaccine evenly usually leads to a larger basic reproduction number and a higher prevalence of infection. Vaccination becomes less effective in reducing the infection once the vaccine amount exceeds the smaller recruit rate of the two genders. In the case study, we estimate the basic reproduction number is 1.0333 for HPV 16/18 in people aged 15-55. The minimal bivalent HPV vaccine needed for the disease prevalence to be below 0.05% is 24050 per year, which should be given to females. However, with this vaccination strategy it would require a very long time and a large amount of vaccine to achieve the goal. In contrast with allocating the same vaccine amount every year, we find that a variable vaccination strategy with more vaccine given in the beginning followed by less vaccine in later years can save time and total vaccine amount. The variable vaccination strategy illustrated in this study can help to better distribute the vaccine to reduce the HPV prevalence. Although this work is for HPV infection and the case study is for a province in China, the model, analysis and conclusions may be applicable to other sexually transmitted diseases in other regions or countries.