A major challenge in human papillomavirus (HPV) vaccine programs is the universal gender-neutral recommendation, as well as estimation of its long-term effect. The objective of this study is to predict the added benefit of male vaccination, especially in men who have sex with men (MSM), and to analyze the impact of the program on society. We propose a mathematical model of the HPV infection based on a network paradigm. Data from Spain allowed constructing the sexual network. HPV force of infection was taken from literature. Different scenarios using variable vaccine coverage in both males and females were studied. Strong herd immunity is shown in the heterosexual population, with an important decrease of HPV 6/11 infections both in men and in unvaccinated women with an only-women vaccination at 14 years of age. No impact of this program occurred in the infection incidence in MSM. This group would only benefit from a vaccination program that includes males. However, the impact at short term would be lower than in heterosexual men. The protection of MSM can only be achieved by direct vaccination of males. This may have important consequences for public health.
Keywords: epidemiologic transition; human papillomavirus; men who have sex with men; network model; sexual behavior; sexual partners; vaccination strategies.