Background

In Australia, high uptake of the quadrivalent human papillomavirus (4vHPV) vaccine has led to reductions in the prevalence of human papillomavirus (HPV) genotypes 6, 11, 16, and 18 in women and girls aged ≤25 years. We evaluated the impact of the program impact on HPV prevalence in unvaccinated male subjects.

Methods

Sexually active heterosexual male subjects aged 16–35 years were recruited in 2014–2016. Participants provided a self-collected penile swab sample for HPV genotyping (Roche Linear Array) and completed a demographic and risk factor questionnaire.

Results

The prevalence of 4vHPV genotypes among 511 unvaccinated male subjects was significantly lower in those aged ≤25 than in those aged >25 years: 3.1% (95% confidence interval, 1.5%–5.7%) versus 13.7% (8.9%–20.1%), respectively (P < .001); adjusted prevalence ratio, 0.22 (.09–.51; P < .001). By contrast, the prevalence of high-risk HPV genotypes other than 16 and 18 remained the same across age groups: 16.8% (95% confidence interval, 12.6%–21.9%) in men aged ≤25 years and 17.9% (12.4%–25.0%) in those aged >25 years (P = .76); adjusted prevalence ratio, 0.98, (.57–1.37; P = .58).

Conclusions

A 78% lower prevalence of 4vHPV genotypes was observed among younger male subjects. These data suggest that unvaccinated men may have benefited from herd protection as much as women from a female-only HPV vaccination program with high coverage.

Human papillomavirus (HPV) vaccination programs have been implemented among adolescent girls and young women in an increasing number of countries using the highly effective prophylactic quadrivalent and/or bivalent HPV vaccines [1]. Both vaccines protect against infection with HPV genotypes 16 and 18, which cause 70%–80% of cervical cancers globally [2, 3], as well as 70% of vaginal cancers, 43% of vulvar cancers, 50% of penile cancers, 88% of anal cancers, and 13%–56% of cancers of the base of tongue and oropharynx [4]. The quadrivalent HPV (4vHPV) vaccine also protects against infection with 2 additional HPV genotypes—6 and 11—that cause 85%–95% of anogenital warts [5].

In Australia, between 2007 and 2009, all girls and women aged 12–26 years were offered free vaccination using a 3-dose course of the 4vHPV vaccine. The program was delivered through schools for girls aged 12–18 years, and through community providers for women up to the age of 26 years. From 2010 the program continued for 12–13-year-old girls, and in 2013, it was extended to include 12–13-year-old boys with a 2-year catch-up for those aged 14–15 years [6, 7]. A National HPV Vaccination Program Register (NHVPR) was established to monitor vaccine coverage and to support the management and evaluation of the program [8]. In 2013, data from the NHVPR showed that 86% of adolescent girls and 78% of boys aged 12–13 years received ≥1 vaccine dose, with 77% and 67%, respectively, receiving all 3 doses [9, 10].

Surveillance data from several countries have provided growing evidence for the population-level benefits of HPV vaccination in young women, both direct, and through herd protection of those who remain unvaccinated. These benefits include rapid and substantial reductions in the prevalence of 4vHPV genotypes (HPV-6/11/16/18) and diagnoses of genital warts [7, 1013]. The demonstrated benefit for young unvaccinated women is likely to be attributable to reductions in the risk of HPV acquisition and hence prevalence among their male sexual partners. Indeed, significant decreases in genital wart diagnoses have been observed in young Australian heterosexual men [9, 14].

A recent study from a single clinic of stored specimens (mostly [97%] urine specimens) collected between 2004 and 2015, from Australian-born men with chlamydia aged <25 years, reported a significant fall in 4vHPV genotypes after the introduction of female vaccination, suggesting herd protection [15]. To further investigate the evidence for male herd protection through the female vaccination program, we estimated the age-specific prevalence of penile HPV genotypes among 16–35-year-old unvaccinated heterosexual male subjects recruited through a network of clinic and community based settings across Australia.

METHODS

Study Population

Participants in this study were years recruited from 2 national cross-sectional HPV prevalence studies conducted between January 2014 and May 2016: the National HPV Monitoring Program (IMPACT) and the Impact of HPV Vaccination Research Study (IMPRESS). IMPACT is an ongoing HPV surveillance system that has established a national network of sentinel sites to monitor circulating HPV genotypes in the Australian population and is aimed at evaluating the effectiveness of the National HPV Vaccination Program. Heterosexual male subjects aged 16–35 years were recruited for HPV testing at 10 participating sites, including 7 sexual health services in Melbourne, Adelaide, Hobart, Cairns, Fremantle, and Sydney and 4 general practice clinics in Melbourne and regional Victoria. In addition, participants were recruited outside the clinic setting through Facebook, using strategies developed for the Vaccine Against Cervical Cancer Impact and Effectiveness (VACCINE) Study, as described elsewhere [16, 17].

IMPRESS is a study of 17–19-year-old heterosexual male subjects recruited from the 7 aforementioned sexual health clinics and through advertisements placed on notice boards at a university and family planning clinic in Melbourne. This study aimed to determine the impact of male HPV vaccination on the prevalence of penile HPV in teenage male subjects.

Subjects were eligible for either study if they reported sex with ≥1 woman during their lifetime. To obtain samples that reflected local prevalence and not that seen in international travelers or recent migrants, inclusion was limited to Australian citizens or men resident in Australia from age 12 years or younger. Participants were excluded if they were unable to provide informed consent or if they reported ever having sex with a man. Ethical approval for this study was granted from 8 human research and ethics committees governing recruitment sites and Facebook, and all participants provided written informed consent.

Procedures

Subjects completed a brief questionnaire on their sociodemographic characteristics, smoking status, age at first sexual intercourse, number of female sexual partners (lifetime and in the previous 12 months), and whether or not they had received HPV vaccination. They were instructed by the research assistant to self-collect a sample for HPV using a moistened flocked swab that was firmly rubbed over the entire surface of the penile shaft, coronal sulcus, glans penis, and if present, the retracted prepuce. Men were also provided with illustrated step-by-step instructions. Sampling of the penile shaft, the glans, and coronal sulcus has been reported to provide the highest detection rates for assessing HPV status because of direct contact with the female genital tract [18]. Furthermore, genital self-sampling methods have been shown to be comparable to sampling by clinicians with respect to specimen adequacy and HPV detection [19, 20].

Written consent was obtained to validate self-reported HPV vaccination status through the NHVPR, which is required to collect individual records on all HPV vaccination doses delivered through schools and primary care in Australia. Participants were classified as being unvaccinated if they reported not receiving the vaccine and if the NHVPR had no record of any HPV vaccine dose being administered. If participants reported receiving the vaccine but had no NHVPR record, registry staff contacted the healthcare provider identified by the participant to ascertain the number of doses and date of vaccination. Participants were excluded if the register was unable to confirm receipt of the vaccine with the provider.

HPV Testing

All samples underwent HPV DNA detection and genotyping at the HPV Regional Labnet, located at the Molecular Microbiology Department, Royal Women's Hospital, Melbourne, Australia. Swab samples were first vigorously agitated in 500 µL of phosphate-buffered saline to release the cellular material. Initially, as an in-house modification, a 100-µL aliquot was added to 1 mL of PreservCyt solution (Hologic) and screened for the presence of 14 high-risk HPV genotypes using the Cobas HPV test (Roche Diagnostics). Concurrently, DNA was extracted from a 200-µL aliquot of the original cell suspension (MagNA Pure 96 DNA and Viral Nucleic Acid Small Volume Kit, Pathogen Universal 200 protocol; Roche Molecular Diagnostics) and eluted in 100 µL.

Extracted DNA was assessed for adequacy by quantitative polymerase chain reaction (PCR) amplification of a 260–base pair segment of the human β-globin gene [21], and a 110–base pair segment to improve detection from keratinized cells, which may contain degraded DNA [22]. DNA extracts from samples negative with the Cobas test were screened for the presence of mucosal HPV DNA by PCR amplification using L1 consensus primer set PGMY09-PGMY11, followed by detection with enzyme-linked immunosorbent assay using a set of biotin-labeled probes [23, 24]. Samples positive for HPV by either screening test were subsequently genotyped using the Linear Array HPV genotyping test (Roche Molecular Diagnostics). Owing to possible cross-reactivity of the HPV-52 probe with types 33, 35, and 58 amplicons, samples positive for ≥1 of these 3 probes were further tested for HPV-52 using a type-specific PCR assay [25].

Statistical Analyses

For men who consented to NHVPR verification and had no evidence of vaccination, HPV prevalence estimates and 95% confidence intervals (CIs) in each of 4 age groups (16–20, 21–25, 26–30, and 31–35 years) were calculated according to the following: detection of any of the 37 HPV genotypes tested; 4vHPV genotypes (HPV-6/11/16/18); any high-risk HPV genotype (HPV-16/18/31/33/35/39/45/51/52/56/58/59/68); any high-risk HPV genotype other than vaccine-targeted genotypes 16 and 18; and the 5 additional HPV types covered by the HPV vaccine (9vHPV; HPV-31/33/45/52/58) [26]. To assess the potential herd effect of female vaccination on penile HPV detection, we compared HPV prevalence between male subjects aged 16–25 and those aged 26–35 years, because (1) population surveys have shown Australian men and their female sexual partners to be of similar age [27, 28] and (2) the low prevalence of vaccine-targeted HPV genotypes in Australian women aged ≤25 years reported in several studies [12, 14, 17], consistent with the higher vaccine coverage achieved in this age group compared with older women [29].

Crude and adjusted prevalence ratios and 95% CIs for each HPV genotype group were calculated to compare prevalence rates between men aged ≤25 years and those aged >25 years. Prevalence ratios were adjusted for source of recruitment; country of birth; education; socioeconomic status (upper or lower 50th centile, based on the Australian Bureau of Statistics Index of Relative Socioeconomic Disadvantage for each individual's residential postal code [12]), residential area (major city or other, based on the Accessibility/Remoteness Index of Australia classification [12]), smoking status, age at first vaginal sex, number of lifetime female partners, and female partners in the previous 12 months. We used χ2 tests to examine associations between age group and sociodemographic and behavioral characteristics. All estimates and 95% CIs were calculated using Poisson regression with robust variance [30]. Statistical significance was tested at the P < .05 level. Analyses were conducted using Stata software, version 14.1 (StataCorp).

RESULTS

From January 2014 to May 2016, a total of 653 men were recruited, including 449 (68.7%) recruited through the IMPACT study and 204 (31.3%) through IMPRESS. Data on 142 men were excluded from further analysis because of nonconsent to access the NHVPR (n = 33); documented receipt of ≥1 HPV vaccine dose (n = 16), or penile sample negative for human β-globin and therefore inadequate for HPV testing (n = 93).

The demographic and sexual behavioral characteristics among the remaining 511 unvaccinated heterosexual men are presented in Table 1. The median age of men included in this study was 23 years (interquartile range, 20–28 years). Overall, 68.7% of men were recruited from sexual health clinics, 42.1% were aged ≤16 years at first vaginal sex, and 66.9% reported ≥2 female sexual partners in the previous 12 months (Table 1).

Table 1.

Demographic and Behavioral Characteristics of 511 Register-Confirmed Unvaccinated Australian Heterosexual Male Subjects

VariableSubjects, No. (%)a
Age group
 16–20 y185 (36.2)
 21–25 y136 (26.6)
 26–30 y108 (21.1)
 31–35 y82 (16.1)
Source of recruitment
 Sexual health clinics351 (68.7)
 Otherb160 (31.3)
Country of birth
 Australia320 (62.6)
 Other183 (35.8)
Educational level
 Secondary or less189 (37.0)
 Tertiary310 (60.8)
Area of residence
 Major cities458 (89.6)
 Other areas42 (8.2)
Socioeconomic status
 More disadvantaged86 (16.8)
 Less disadvantaged414 (81.0)
Current smoker
 No362 (70.8)
 Yes139 (27.2)
Age at first vaginal sex
 ≤16 y215 (42.1)
 >16 y285 (55.8)
Lifetime No. of female partners
 <5147 (28.8)
 5–10122 (23.9)
 11–2094 (18.4)
 >20135 (26.4)
No. of female partners in past 12 mo
 013 (2.5)
 1149 (29.2)
 ≥2342 (66.9)
VariableSubjects, No. (%)a
Age group
 16–20 y185 (36.2)
 21–25 y136 (26.6)
 26–30 y108 (21.1)
 31–35 y82 (16.1)
Source of recruitment
 Sexual health clinics351 (68.7)
 Otherb160 (31.3)
Country of birth
 Australia320 (62.6)
 Other183 (35.8)
Educational level
 Secondary or less189 (37.0)
 Tertiary310 (60.8)
Area of residence
 Major cities458 (89.6)
 Other areas42 (8.2)
Socioeconomic status
 More disadvantaged86 (16.8)
 Less disadvantaged414 (81.0)
Current smoker
 No362 (70.8)
 Yes139 (27.2)
Age at first vaginal sex
 ≤16 y215 (42.1)
 >16 y285 (55.8)
Lifetime No. of female partners
 <5147 (28.8)
 5–10122 (23.9)
 11–2094 (18.4)
 >20135 (26.4)
No. of female partners in past 12 mo
 013 (2.5)
 1149 (29.2)
 ≥2342 (66.9)

a Totals do not always equal 511 because of a small amount of missing data. All men were confirmed to be unvaccinated with the National HPV Program Register.

b Other sources of recruitment included Facebook (19.0%), notice boards at a university (6.5%), and general practice (4.7%), and family planning (1.4%) clinics.

Table 1.

Demographic and Behavioral Characteristics of 511 Register-Confirmed Unvaccinated Australian Heterosexual Male Subjects

VariableSubjects, No. (%)a
Age group
 16–20 y185 (36.2)
 21–25 y136 (26.6)
 26–30 y108 (21.1)
 31–35 y82 (16.1)
Source of recruitment
 Sexual health clinics351 (68.7)
 Otherb160 (31.3)
Country of birth
 Australia320 (62.6)
 Other183 (35.8)
Educational level
 Secondary or less189 (37.0)
 Tertiary310 (60.8)
Area of residence
 Major cities458 (89.6)
 Other areas42 (8.2)
Socioeconomic status
 More disadvantaged86 (16.8)
 Less disadvantaged414 (81.0)
Current smoker
 No362 (70.8)
 Yes139 (27.2)
Age at first vaginal sex
 ≤16 y215 (42.1)
 >16 y285 (55.8)
Lifetime No. of female partners
 <5147 (28.8)
 5–10122 (23.9)
 11–2094 (18.4)
 >20135 (26.4)
No. of female partners in past 12 mo
 013 (2.5)
 1149 (29.2)
 ≥2342 (66.9)
VariableSubjects, No. (%)a
Age group
 16–20 y185 (36.2)
 21–25 y136 (26.6)
 26–30 y108 (21.1)
 31–35 y82 (16.1)
Source of recruitment
 Sexual health clinics351 (68.7)
 Otherb160 (31.3)
Country of birth
 Australia320 (62.6)
 Other183 (35.8)
Educational level
 Secondary or less189 (37.0)
 Tertiary310 (60.8)
Area of residence
 Major cities458 (89.6)
 Other areas42 (8.2)
Socioeconomic status
 More disadvantaged86 (16.8)
 Less disadvantaged414 (81.0)
Current smoker
 No362 (70.8)
 Yes139 (27.2)
Age at first vaginal sex
 ≤16 y215 (42.1)
 >16 y285 (55.8)
Lifetime No. of female partners
 <5147 (28.8)
 5–10122 (23.9)
 11–2094 (18.4)
 >20135 (26.4)
No. of female partners in past 12 mo
 013 (2.5)
 1149 (29.2)
 ≥2342 (66.9)

a Totals do not always equal 511 because of a small amount of missing data. All men were confirmed to be unvaccinated with the National HPV Program Register.

b Other sources of recruitment included Facebook (19.0%), notice boards at a university (6.5%), and general practice (4.7%), and family planning (1.4%) clinics.

The prevalence of penile HPV genotypes stratified by age group (16–20, 21–25, 26–30, or 31–35 years) are shown in Table 2 and illustrated in Figure 1. The combined prevalence of 4vHPV-targeted genotypes was low, at 2.7% and 3.7%, respectively, among men in the 16–20- and 21–25-year age groups, increasing significantly to 16.7% and 9.8% among those aged 25–29 or 30–35 years (Ptrend = .001). The prevalence of low-risk HPV genotypes 6 and 11 was 1.1% and 1.5% among men aged 16–20 or 21–25 years, increasing significantly to 9.3% and 6.1% among those aged 26–30 or 31–35 years, respectively (Ptrend = .003). The prevalence of high-risk HPV genotypes 16 and 18 was 1.6% and 2.2% among men aged 16–20 and 21–25 years, increasing significantly to 9.3% and 4.9% among those aged 26–30 and 31–35 years, respectively (Ptrend = .02). By contrast, the prevalence of high-risk HPV genotypes other than HPV-16/18 did not change with increasing age group (Ptrend = .77). There was no significant trend across age groups in the prevalence of any HPV genotype (Ptrend = .11) or in the prevalence of additional HPV genotypes that a 9vHPV vaccine would provide protection against (Ptrend= .50).
Table 2.

Prevalence of Penile HPV Genotypes by Age Group Among 511 Register-Confirmed Unvaccinated Australian Heterosexual Male Subjects

HPV GenotypeAge 16–20 y (n = 185)
Age 21–25 y (n = 136)
Age 26–30 y (n = 108)
Age 31–35 y (n = 82)
Ptrend
No.% (95% CI)No.% (95% CI)No.% (95% CI)No.% (95% CI)
Any3921.1 (15.0–28.8)4230.1 (22.3–41.7)4037.0 (26.5–50.4)2226.8 (16.8–40.6).11
Vaccine targeted (6, 11, 16, or 18)52.7 (.9–6.3)53.7 (1.2–8.6)1816.7 (9.9–26.3)89.8 (4.2–19.2).001
Any high risk3016.2 (10.9–23.1)2619.1 (12.5–28.0)2825.9 (17.2–37.5)1518.3 (10.2–30.2).30
High risk, other than 16 and 182815.1 (10.1–21.9)2619.1 (12.5–28.0)2220.4 (12.8–30.8)1214.6 (7.6–25.6).77
6 or 1121.1 (.1–3.9)21.5 (.2–5.3)109.3 (4.4–17.0)56.1 (2.0–14.2).003
16 or 1831.6 (.3–4.7)32.2 (.5–6.4)109.3 (4.4–17.0)44.9 (1.3–12.5).02
31, 33, 45, 52, or 58126.5 (3.4–11.3)1611.8 (6.7–19.1)87.4 (3.2–14.6)33.7 (.8–10.7).50
HPV GenotypeAge 16–20 y (n = 185)
Age 21–25 y (n = 136)
Age 26–30 y (n = 108)
Age 31–35 y (n = 82)
Ptrend
No.% (95% CI)No.% (95% CI)No.% (95% CI)No.% (95% CI)
Any3921.1 (15.0–28.8)4230.1 (22.3–41.7)4037.0 (26.5–50.4)2226.8 (16.8–40.6).11
Vaccine targeted (6, 11, 16, or 18)52.7 (.9–6.3)53.7 (1.2–8.6)1816.7 (9.9–26.3)89.8 (4.2–19.2).001
Any high risk3016.2 (10.9–23.1)2619.1 (12.5–28.0)2825.9 (17.2–37.5)1518.3 (10.2–30.2).30
High risk, other than 16 and 182815.1 (10.1–21.9)2619.1 (12.5–28.0)2220.4 (12.8–30.8)1214.6 (7.6–25.6).77
6 or 1121.1 (.1–3.9)21.5 (.2–5.3)109.3 (4.4–17.0)56.1 (2.0–14.2).003
16 or 1831.6 (.3–4.7)32.2 (.5–6.4)109.3 (4.4–17.0)44.9 (1.3–12.5).02
31, 33, 45, 52, or 58126.5 (3.4–11.3)1611.8 (6.7–19.1)87.4 (3.2–14.6)33.7 (.8–10.7).50

Abbreviations: CI, confidence interval; HPV, human papillomavirus.

Table 2.

Prevalence of Penile HPV Genotypes by Age Group Among 511 Register-Confirmed Unvaccinated Australian Heterosexual Male Subjects

HPV GenotypeAge 16–20 y (n = 185)
Age 21–25 y (n = 136)
Age 26–30 y (n = 108)
Age 31–35 y (n = 82)
Ptrend
No.% (95% CI)No.% (95% CI)No.% (95% CI)No.% (95% CI)
Any3921.1 (15.0–28.8)4230.1 (22.3–41.7)4037.0 (26.5–50.4)2226.8 (16.8–40.6).11
Vaccine targeted (6, 11, 16, or 18)52.7 (.9–6.3)53.7 (1.2–8.6)1816.7 (9.9–26.3)89.8 (4.2–19.2).001
Any high risk3016.2 (10.9–23.1)2619.1 (12.5–28.0)2825.9 (17.2–37.5)1518.3 (10.2–30.2).30
High risk, other than 16 and 182815.1 (10.1–21.9)2619.1 (12.5–28.0)2220.4 (12.8–30.8)1214.6 (7.6–25.6).77
6 or 1121.1 (.1–3.9)21.5 (.2–5.3)109.3 (4.4–17.0)56.1 (2.0–14.2).003
16 or 1831.6 (.3–4.7)32.2 (.5–6.4)109.3 (4.4–17.0)44.9 (1.3–12.5).02
31, 33, 45, 52, or 58126.5 (3.4–11.3)1611.8 (6.7–19.1)87.4 (3.2–14.6)33.7 (.8–10.7).50
HPV GenotypeAge 16–20 y (n = 185)
Age 21–25 y (n = 136)
Age 26–30 y (n = 108)
Age 31–35 y (n = 82)
Ptrend
No.% (95% CI)No.% (95% CI)No.% (95% CI)No.% (95% CI)
Any3921.1 (15.0–28.8)4230.1 (22.3–41.7)4037.0 (26.5–50.4)2226.8 (16.8–40.6).11
Vaccine targeted (6, 11, 16, or 18)52.7 (.9–6.3)53.7 (1.2–8.6)1816.7 (9.9–26.3)89.8 (4.2–19.2).001
Any high risk3016.2 (10.9–23.1)2619.1 (12.5–28.0)2825.9 (17.2–37.5)1518.3 (10.2–30.2).30
High risk, other than 16 and 182815.1 (10.1–21.9)2619.1 (12.5–28.0)2220.4 (12.8–30.8)1214.6 (7.6–25.6).77
6 or 1121.1 (.1–3.9)21.5 (.2–5.3)109.3 (4.4–17.0)56.1 (2.0–14.2).003
16 or 1831.6 (.3–4.7)32.2 (.5–6.4)109.3 (4.4–17.0)44.9 (1.3–12.5).02
31, 33, 45, 52, or 58126.5 (3.4–11.3)1611.8 (6.7–19.1)87.4 (3.2–14.6)33.7 (.8–10.7).50

Abbreviations: CI, confidence interval; HPV, human papillomavirus.

Prevalence of penile human papillomavirus (HPV) genotypes by age group among 511 register-confirmed unvaccinated Australian heterosexual male subjects.
Figure 1.

Prevalence of penile human papillomavirus (HPV) genotypes by age group among 511 register-confirmed unvaccinated Australian heterosexual male subjects.

Compared with men aged >25 years, those aged ≤25 years were less likely to be Australian born (P < .001), less likely to be tertiary educated (P < .001), more likely to report younger age at first vaginal sex (P = .048), and reported fewer female sexual partners in their lifetime (P < .001). Men in each age group were similar with respect to source of recruitment, area of residency, socioeconomic status, smoking, and number of female partners in the previous 12 months.

Crude and adjusted prevalence ratios for penile HPV genotypes, adjusted for source of recruitment, country of birth, education, lifetime female partners, and age at first vaginal sex, are presented in Table 3. The combined prevalence of 4vHPV-targeted genotypes was significantly lower in men aged ≤25 years than in those aged >25 years (3.1% vs 13.7%, respectively; P < .001), with an adjusted prevalence ratio of 0.22 (95% CI, .09–.51; P < .001). The adjusted prevalence ratios were 0.14 (95% CI, .04–.62; P = .009) for low-risk HPV genotypes 6 and 11, and 0.29 (95% CI, .11–.74; P = .009) for high-risk HPV genotypes 16 and 18 when these categories were examined separately. By contrast, the prevalence of high risk-HPV genotypes other than HPV-16/18 did not change with increasing age (P = .76), with an adjusted prevalence ratio of 0.98 (95% CI, .57–1.37; P = .58). Finally, the adjusted prevalence ratios for any HPV and for additional HPV genotypes that the 9vHPV vaccine would provide protection against were 0.79 (95% CI, .57–1.10; P = .16) and 2.09 (95% CI, .99–4.45; P = .054), respectively.

Table 3.

Crude PRs and aPRs for Penile HPV Genotypes Among 511 Register-Confirmed Unvaccinated Australian Heterosexual Male Subjects

HPV GenotypeAge ≤25 y (n = 321)
Age >25 y (n = 190)
PR (95% CI)P ValueaPRa (95% CI)P Value
No.% (95% CI)No.% (95% CI)
Any8125.2 (20.0–31.4)6232.6 (25.0–41.8)0.77 (.59–1.02).070.79 (.57–1.10).16
Vaccine targeted (6, 11, 16, and 18)103.1 (1.5–5.7)2613.7 (8.9–20.1)0.23 (.11–.46)<.0010.22 (.09–.51)<.001
Any high risk5617.5 (13.2–22.7)4322.6 (16.4–30.5)0.77 (.54–1.10).150.77 (.51–1.16).21
High risk, other than 16 and 185416.8 (12.6–21.9)3417.9 (12.4–25.0)0.94 (.64–1.39).760.98 (.57–1.37).58
6 or 1141.2 (0.3–3.2)157.9 (4.4–13.0)0.16 (.05–.47).0010.14 (.04–.62).009
16 or 1861.9 (0.06–4.1)147.4 (4.0–12.4)0.25 (.10–.65).0040.29 (.11–.74).009
31, 33, 45, 52, or 58288.7 (5.8–12.6)115.8 (2.9–10.4)1.51 (.77–2.96).232.09 (.99–4.45).054
HPV GenotypeAge ≤25 y (n = 321)
Age >25 y (n = 190)
PR (95% CI)P ValueaPRa (95% CI)P Value
No.% (95% CI)No.% (95% CI)
Any8125.2 (20.0–31.4)6232.6 (25.0–41.8)0.77 (.59–1.02).070.79 (.57–1.10).16
Vaccine targeted (6, 11, 16, and 18)103.1 (1.5–5.7)2613.7 (8.9–20.1)0.23 (.11–.46)<.0010.22 (.09–.51)<.001
Any high risk5617.5 (13.2–22.7)4322.6 (16.4–30.5)0.77 (.54–1.10).150.77 (.51–1.16).21
High risk, other than 16 and 185416.8 (12.6–21.9)3417.9 (12.4–25.0)0.94 (.64–1.39).760.98 (.57–1.37).58
6 or 1141.2 (0.3–3.2)157.9 (4.4–13.0)0.16 (.05–.47).0010.14 (.04–.62).009
16 or 1861.9 (0.06–4.1)147.4 (4.0–12.4)0.25 (.10–.65).0040.29 (.11–.74).009
31, 33, 45, 52, or 58288.7 (5.8–12.6)115.8 (2.9–10.4)1.51 (.77–2.96).232.09 (.99–4.45).054

Abbreviations: aPR, adjusted prevalence ratio; CI, confidence intervals; HPV, human papillomavirus; PR, prevalence ratio.

a Adjusted for source of recruitment, country of birth, highest level of education, lifetime number of female partners, and age at first vaginal sex. The older group was coded as the reference category.

Table 3.

Crude PRs and aPRs for Penile HPV Genotypes Among 511 Register-Confirmed Unvaccinated Australian Heterosexual Male Subjects

HPV GenotypeAge ≤25 y (n = 321)
Age >25 y (n = 190)
PR (95% CI)P ValueaPRa (95% CI)P Value
No.% (95% CI)No.% (95% CI)
Any8125.2 (20.0–31.4)6232.6 (25.0–41.8)0.77 (.59–1.02).070.79 (.57–1.10).16
Vaccine targeted (6, 11, 16, and 18)103.1 (1.5–5.7)2613.7 (8.9–20.1)0.23 (.11–.46)<.0010.22 (.09–.51)<.001
Any high risk5617.5 (13.2–22.7)4322.6 (16.4–30.5)0.77 (.54–1.10).150.77 (.51–1.16).21
High risk, other than 16 and 185416.8 (12.6–21.9)3417.9 (12.4–25.0)0.94 (.64–1.39).760.98 (.57–1.37).58
6 or 1141.2 (0.3–3.2)157.9 (4.4–13.0)0.16 (.05–.47).0010.14 (.04–.62).009
16 or 1861.9 (0.06–4.1)147.4 (4.0–12.4)0.25 (.10–.65).0040.29 (.11–.74).009
31, 33, 45, 52, or 58288.7 (5.8–12.6)115.8 (2.9–10.4)1.51 (.77–2.96).232.09 (.99–4.45).054
HPV GenotypeAge ≤25 y (n = 321)
Age >25 y (n = 190)
PR (95% CI)P ValueaPRa (95% CI)P Value
No.% (95% CI)No.% (95% CI)
Any8125.2 (20.0–31.4)6232.6 (25.0–41.8)0.77 (.59–1.02).070.79 (.57–1.10).16
Vaccine targeted (6, 11, 16, and 18)103.1 (1.5–5.7)2613.7 (8.9–20.1)0.23 (.11–.46)<.0010.22 (.09–.51)<.001
Any high risk5617.5 (13.2–22.7)4322.6 (16.4–30.5)0.77 (.54–1.10).150.77 (.51–1.16).21
High risk, other than 16 and 185416.8 (12.6–21.9)3417.9 (12.4–25.0)0.94 (.64–1.39).760.98 (.57–1.37).58
6 or 1141.2 (0.3–3.2)157.9 (4.4–13.0)0.16 (.05–.47).0010.14 (.04–.62).009
16 or 1861.9 (0.06–4.1)147.4 (4.0–12.4)0.25 (.10–.65).0040.29 (.11–.74).009
31, 33, 45, 52, or 58288.7 (5.8–12.6)115.8 (2.9–10.4)1.51 (.77–2.96).232.09 (.99–4.45).054

Abbreviations: aPR, adjusted prevalence ratio; CI, confidence intervals; HPV, human papillomavirus; PR, prevalence ratio.

a Adjusted for source of recruitment, country of birth, highest level of education, lifetime number of female partners, and age at first vaginal sex. The older group was coded as the reference category.

Finally, we examined the distribution of genotypes among HPV-positive men only, under the assumption that even if the absolute prevalence changed by age group owing to changes in sexual activity, the distribution of types among HPV-positive men should not. Similar patterns to that seen in the overall study population were observed (Supplementary Figure 1 and Table 1).

DISCUSSION

This is the first study to our knowledge to investigate herd protection by comparing prevalence of HPV in younger unvaccinated men, who are likely to have received the greatest benefit from vaccination of women in their age cohort, with that of older men, whose female partners are less likely to have been vaccinated. We found the prevalence of penile HPV genotypes targeted by the 4vHPV vaccine (HPV-6/11/16/18) was very low (3%) among male subjects aged ≤25 years, 78% lower than that among those aged >25 years (14%). By contrast, there was no difference in the prevalence of other high-risk HPV genotypes not targeted by the vaccine.

These differences in HPV prevalence could not be explained by differences in other characteristics that might influence rates of HPV infection. Given the absence of a difference in the prevalence for nonvaccine genotypes, the most likely explanation is the herd effect of HPV vaccination of girls and women. Our study is the first to provide evidence from multiple sites, suggesting that herd protection of men by female HPV vaccination has occurred across Australia. With the introduction of male vaccination, and with cohorts of vaccinated women now covering an increased age range, our study also provides baseline data for ongoing monitoring of the impact of the National HPV Vaccination Program on HPV infection in men.

Australia has relatively high HPV vaccine coverage of women, with 74% coverage for 1 dose and 55% for 3 doses among 18–26-year-old women in 2011 [31]. Previous Australian sentinel surveillance data have shown that this vaccine coverage was sufficient to reduce the combined prevalence of 4vHPV-targeted genotypes by 78% in the overall study population and 95% among fully vaccinated women [12]. In the same study, reductions in prevalence of 75% and 35% were found in partially vaccinated and unvaccinated women, respectively [12]. The 78% lower combined prevalence of the vaccine-targeted HPV genotypes among younger men in this study suggests that reductions in HPV prevalence in unvaccinated men may be of a similar magnitude to the reductions achieved in women.

No comparison data are available on the age-specific prevalence of HPV among men in Australia before the implementation of the female vaccination program. One hypothetical explanation for the higher prevalence of vaccine-targeted HPV genotypes in older men observed in our study is that HPV clears more slowly with increasing age. However, data from a large prospective study of genital HPV natural history suggest this is unlikely to be the case [32]. Moreover, the relatively flat age-specific prevalence curves observed in our study for HPV genotypes not targeted by the vaccine are consistent with patterns of HPV by age, reported in multiple studies of the epidemiology of genital HPV infection in heterosexual men [33, 34].

A number of factors should be considered when interpreting the results from this study. First, we excluded 15% of participants with penile swab samples that were inadequate for HPV testing [19, 20]. This rate of unassessable samples was higher than previously reported [19] and probably represents inadequate self-sampling by some participants. Nevertheless, these men did not differ from those included in this study with respect to demographic or HPV risk factors. Furthermore, the proportions of younger versus older men with penile samples inadequate for HPV testing did not differ (15% vs 16%, respectively; P = .78).

Second, our estimate of the extent of herd protection, which is based on a cross-sectional comparison of HPV prevalence across age groups, is likely to underestimate the true magnitude of benefit to men from herd effects [35]. This is because men in the older group are also likely to have already benefited from herd protection. An estimated 55% of 18–26-year-old women, corresponding to those aged 26–34 years in 2015, received ≥1 dose of the HPV vaccine in Australia's 2007–2009 catch up program, with an estimated 32% receiving all 3 doses [29]. Third, whereas men were recruited from diverse sources, most were sourced from sexual health clinics, which may have biased sampling toward higher-risk men and therefore a higher prevalence of HPV infection. Nevertheless, even if absolute prevalence rates are not generalizable to all Australian heterosexual men, the comparison between age-specific prevalence may well be. Finally, the results may not be generalizable to countries where female vaccine uptake is lower or among subgroups such as men who have sex with men [36, 37].

There is ongoing global debate as to whether boys should be included in national HPV vaccination programs. At the time of this study only a handful of countries (Australia, Austria, Israel, the United States, and 6 provinces of Canada) had implemented male HPV vaccination programs in addition to existing female vaccination programs. Cost-effectiveness studies have been conducted to help determine whether any further health benefits would be obtained from the addition of male HPV vaccination [38]. However, these have been based on mathematical modeling without actual male prevalence data such as those presented in this study. The World Health Organization recommends the establishment of sentinel surveillance to monitor the impact of HPV vaccination programs on HPV prevalence [39]. With the introduction of male vaccination in Australia, ongoing surveillance of HPV prevalence with linkage to the national vaccine register will allow monitoring of the impact of HPV vaccination over time, informing whether routine vaccination of boys will result in any benefits over and above those that have already been achieved through the existing female vaccination program.

Notes

Acknowledgments. The National HPV Monitoring Program (IMPACT) is led by the World Health Organization Reference Laboratory for HPV detection and genotyping, located at the Royal Women's Hospital in collaboration with a consortium of member organizations. The member organizations include the Kirby Institute at the University of New South Wales, the University of Melbourne, Melbourne Sexual Health Centre Alfred Health, Family Planning Victoria, Family Planning New South Wales, and the University of Sydney, Department of Child Health. The IMPACT study team includes Deborah Bateson, Catriona Bradshaw, M. Y. C., A. M. C., Basil Donovan, C. K. F., S. M. G., Andrew Grulich, Richard Hillman, J. S. H., Fengyi Jin, J. M. K., Marlene Kong, Matthew Law, D. A. M., Skye McGregor, Kathleen McNamee, Robert Monaghan, Samuel Phillips, Marin Poljak, Mary Poynten, David Regan, S. Rachel Skinner, Mary Stewart, H. W., R. W., S. N. T. (principal investigator for IMPACT), and the external advisory group: Linda Selvey, Vicky Sheppeard, David Smith, Kaushi Kogar, and Marilyn Jodie Clarke. The IMPRESS Study is led by the Melbourne Sexual Health Centre, Alfred Health. The study team includes R. W., Catriona Bradshaw, Julia Brotherton, C. K. F., S. M. G., J. M. K., J. S. H., David Regan, S. N. T., and M. Y. C. (principal investigator for IMPRESS). We would like to thank all the site investigators, C. O., L. M., C. B., A. M., Louise Owen, Darren Russell, and the clinic staff: Melissa Power, Jane Gilbert, Carol Pavitt, Emma Clements, El Thompson, Barb Lennox, Fiona Anderson, Collette Cashman, Michelle Andrews, Alison Beverley, and Mandy Johnson. We also thank Edmund Molesworth, Belinda Hengel, Kathy Rowed, Samantha Sukkel, and Jessica Nia for their assistance in recruitment; and Jennifer Brosi, Lisette Bicknell, and Chantal Kim at the National HPV Vaccination Program Register for assistance with validating vaccination status. Finally, we thank all the IMPACT and IMPRESS participants.

Author contribution. All authors contributed to the study design. D. A. M., E. P. F. C., M. Y. C., and S. N. T. contributed to data analysis. D. A. M., E. P. F. C., J. M. K., M. Y. C., and S. N. T. contributed to data interpretation, and writing of the report. S. M. G., R. W., A. M. C., C. K. F., J. S. H., H. W., A. M., C. B., L. M., and C. O. contributed to site and participant enrollment and reviewed the manuscript for important intellectual content. All authors approved the final version of the manuscript.

Disclaimer. The views expressed in this publication are those of the authors and do not necessarily represent the position of the Australian Government nor those of Merck Sharp & Dohme.

Financial support.The IMPRESS Study is supported by a Merck Investigator Initiated Studies Program (protocol 50939). The National HPV Monitoring Program is supported by the Australian Government Department of Health (HPV Surveillance Fund H1314G010). E. P. F. C. is supported by the Australian National Health and Medical Research Council (Early Career Fellowship 1091226).

Potential conflicts of interest. D. A. M. and E. P. F. C. have each received educational grants from Seqirus (formerly bioCSL) to assist with education, training, and academic purposes in the area of HPV. S. M. G. reports grants from Merck, GlaxoSmithKline, CSL, and the Commonwealth Department of Health outside the submitted work; has received nonfinancial support from Merck; and has delivered lectures and received speaking fees from Merck Sharp & Dohme and Sanofi Pasteur MSD for work performed in her personal time. C. K. F. has received research funding from Merck outside the submitted work, and owns shares in CSL Biotherapies, the manufacturer of Gardasil. M. Y. C. has received grants from Merck Sharp & Dohme, during the conduct of the study. S. N. T. has received grants from bioCSL, outside the submitted work. All other authors report no potential conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed.

References

1

Markowitz
LE
,
Tsu
V
,
Deeks
SL
et al. .
Human papillomavirus vaccine introduction—the first five years
.
Vaccine
2012
;
30
(
suppl 5
):
F139
48
.

2

de Sanjose
S
,
Alemany
L
,
Ordi
J
et al. .
Worldwide human papillomavirus genotype attribution in over 2000 cases of intraepithelial and invasive lesions of the vulva
.
Eur J Cancer
2013
;
49
:
3450
61
.

3

de Sanjose
S
,
Quint
WG
,
Alemany
L
et al. .
Human papillomavirus genotype attribution in invasive cervical cancer: a retrospective cross-sectional worldwide study
.
Lancet Oncol
2010
;
11
:
1048
56
.

4

Forman
D
,
de Martel
C
,
Lacey
CJ
et al. .
Global burden of human papillomavirus and related diseases
.
Vaccine
2012
;
30
(
suppl 5
):
F12
23
.

5

Garland
SM
,
Steben
M
,
Sings
HL
et al. .
Natural history of genital warts: analysis of the placebo arm of 2 randomized phase III trials of a quadrivalent human papillomavirus (types 6, 11, 16, and 18) vaccine
.
J Infect Dis
2009
;
199
:
805
14
.

6

Brotherton
JM
.
Human papillomavirus vaccination: where are we now?
J Paediatr Child Health
2014
;
50
:
959
65
.

7

Markowitz
LE
,
Liu
G
,
Hariri
S
,
Steinau
M
,
Dunne
EF
,
Unger
ER
.
Prevalence of HPV after introduction of the vaccination program in the United States
.
Pediatrics
2016
;
137
:
1
9
.

8

Gertig
DM
,
Brotherton
JM
,
Saville
M
.
Measuring human papillomavirus (HPV) vaccination coverage and the role of the National HPV Vaccination Program Register, Australia
.
Sex Health
2011
;
8
:
171
8
.

9

Ali
H
,
Donovan
B
,
Wand
H
et al. .
Genital warts in young Australians five years into national human papillomavirus vaccination programme: national surveillance data
.
BMJ
2013
;
346
:
f2032
.

10

Smith
MA
,
Liu
B
,
McIntyre
P
,
Menzies
R
,
Dey
A
,
Canfell
K
.
Fall in genital warts diagnoses in the general and indigenous Australian population following implementation of a national human papillomavirus vaccination program: analysis of routinely collected national hospital data
.
J Infect Dis
2015
;
211
:
91
9
.

11

Drolet
M
,
Benard
E
,
Boily
MC
et al. .
Population-level impact and herd effects following human papillomavirus vaccination programmes: a systematic review and meta-analysis
.
Lancet Infect Dis
2015
;
15
:
565
80
.

12

Tabrizi
SN
,
Brotherton
JM
,
Kaldor
JM
et al. .
Assessment of herd immunity and cross-protection after a human papillomavirus vaccination programme in Australia: a repeat cross-sectional study
.
Lancet Infect Dis
2014
;
14
:
958
66
.

13

Hariri
S
,
Bennett
NM
,
Niccolai
LM
et al. .
Reduction in HPV 16/18-associated high grade cervical lesions following HPV vaccine introduction in the United States - 2008–2012
.
Vaccine
2015
;
33
:
1608
13
.

14

Chow
EP
,
Read
TR
,
Wigan
R
et al. .
Ongoing decline in genital warts among young heterosexuals 7 years after the Australian human papillomavirus (HPV) vaccination programme
.
Sex Transm Infect
2015
;
91
:
214
9
.

15

Chow EP, Machalek DA, Tabrizi SN, et al. Quadrivalent vaccine targeted human 44 papillomavirus genotypes in heterosexual men following the Australian female HPV 45 vaccination programme: an observational study from 2004 to 2015. Lancet Infect Dis 2016; doi:10.1016/S1473-3099(16)30116-5
.

16

Fenner
Y
,
Garland
SM
,
Moore
EE
et al. .
Web-based recruiting for health research using a social networking site: an exploratory study
.
J Med Internet Res
2012
;
14
:
e20
.

17

Osborne
SL
,
Tabrizi
SN
,
Brotherton
JM
et al. .
Assessing genital human papillomavirus genoprevalence in young Australian women following the introduction of a national vaccination program
.
Vaccine
2015
;
33
:
201
8
.

18

Giuliano
AR
,
Nielson
CM
,
Flores
R
et al. .
The optimal anatomic sites for sampling heterosexual men for human papillomavirus (HPV) detection: the HPV detection in men study
.
J Infect Dis
2007
;
196
:
1146
52
.

19

Hernandez
BY
,
McDuffie
K
,
Goodman
MT
et al. .
Comparison of physician- and self-collected genital specimens for detection of human papillomavirus in men
.
J Clin Microbiol
2006
;
44
:
513
7
.

20

Hernandez
BY
,
Wilkens
LR
,
Unger
ER
et al. .
Evaluation of genital self-sampling methods for HPV detection in males
.
J Clin Virol
2013
;
58
:
168
75
.

21

Resnick
RM
,
Cornelissen
MT
,
Wright
DK
et al. .
Detection and typing of human papillomavirus in archival cervical cancer specimens by DNA amplification with consensus primers
.
J Natl Cancer Inst
1990
;
82
:
1477
84
.

22

Cornall
AM
,
Roberts
JM
,
Garland
SM
,
Hillman
RJ
,
Grulich
AE
,
Tabrizi
SN
.
Anal and perianal squamous carcinomas and high-grade intraepithelial lesions exclusively associated with “low-risk” HPV genotypes 6 and 11
.
Int J Cancer
2013
;
133
:
2253
8
.

23

Gravitt
PE
,
Peyton
CL
,
Alessi
TQ
et al. .
Improved amplification of genital human papillomaviruses
.
J Clin Microbiol
2000
;
38
:
357
61
.

24

Layton-Henry
J
,
Scurry
J
,
Planner
R
et al. .
Cervical adenoid basal carcinoma, five cases and literature review
.
Int J Gynecol Cancer
1996
;
6
:
193
9
.

25

Stevens
MP
,
Garland
SM
,
Tabrizi
SN
.
Development and validation of a real-time PCR assay specifically detecting human papillomavirus 52 using the Roche LightCycler 480 system
.
J Virol Methods
2008
;
147
:
290
6
.

26

Joura
EA
,
Giuliano
AR
,
Iversen
OE
et al. .
A 9-valent HPV vaccine against infection and intraepithelial neoplasia in women
.
N Engl J Med
2015
;
372
:
711
23
.

27

Rissel
C
,
Badcock
PB
,
Smith
AM
et al. .
Heterosexual experience and recent heterosexual encounters among Australian adults: the Second Australian Study of Health and Relationships
.
Sex Health
2014
;
11
:
416
26
.

28

Chow
EP
,
Read
TR
,
Law
MG
,
Chen
MY
,
Bradshaw
CS
,
Fairley
CK
.
Assortative sexual mixing patterns in male-female and male-male partnerships in Melbourne, Australia: implications for HIV and sexually transmissible infection transmission
.
Sex Health
2016
; .

29

Brotherton
JM
,
Liu
B
,
Donovan
B
,
Kaldor
JM
,
Saville
M
.
Human papillomavirus (HPV) vaccination coverage in young Australian women is higher than previously estimated: independent estimates from a nationally representative mobile phone survey
.
Vaccine
2014
;
32
:
592
7
.

30

Coutinho
LM
,
Scazufca
M
,
Menezes
PR
.
Methods for estimating prevalence ratios in cross-sectional studies
.
Rev Saude Publica
2008
;
42
:
992
8
.

31

Brotherton
J
,
Gertig
D
,
Chappell
G
,
Rowlands
L
,
Saville
M
.
Catching up with the catch-up: HPV vaccination coverage data for Australian women aged 18–26 years from the National HPV Vaccination Program Register
.
Commun Dis Intell Q Rep
2011
;
35
:
197
201
.

32

Giuliano
AR
,
Lee
JH
,
Fulp
W
et al. .
Incidence and clearance of genital human papillomavirus infection in men (HIM): a cohort study
.
Lancet
2011
;
377
:
932
40
.

33

Smith
JS
,
Gilbert
PA
,
Melendy
A
,
Rana
RK
,
Pimenta
JM
.
Age-specific prevalence of human papillomavirus infection in males: a global review
.
J Adolesc Health
2011
;
48
:
540
52
.

34

Giuliano
AR
,
Lazcano-Ponce
E
,
Villa
LL
et al. .
The human papillomavirus infection in men study: human papillomavirus prevalence and type distribution among men residing in Brazil, Mexico, and the United States
.
Cancer Epidemiol Biomarkers Prev
2008
;
17
:
2036
43
.

35

Garland
SM
,
Molesworth
EG
,
Machalek
DA
,
Cornall
AM
,
Tabrizi
SN
.
How to best measure the effectiveness of male human papillomavirus vaccine programmes?
Clin Microbiol Infect
2015
;
21
:
834
41
.

36

Machalek
DA
,
Poynten
M
,
Jin
F
et al. .
Anal human papillomavirus infection and associated neoplastic lesions in men who have sex with men: a systematic review and meta-analysis
.
Lancet Oncol
2012
;
13
:
487
500
.

37

Zou
H
,
Tabrizi
SN
,
Grulich
AE
et al. .
Early acquisition of anogenital human papillomavirus among teenage men who have sex with men
.
J Infect Dis
2014
;
209
:
642
51
.

38

Georgousakis
M
,
Jayasinghe
S
,
Brotherton
J
,
Gilroy
N
,
Chiu
C
,
Macartney
K
.
Population-wide vaccination against human papillomavirus in adolescent boys: Australia as a case study
.
Lancet Infect Dis
2012
;
12
:
627
34
.

39

World Health Organization, Department of Immunization, Vaccines and Biologicals
.
Report of the meeting on HPV Vaccine Coverage and Impact Monitoring
. .
Accessed 5 March 2016
.

Author notes

Presented in part: European Research Organisation on Genital Infection and Neoplasia (EUROGIN) Conference, Salzburg, Austria, 15–18 July 2016.

a

D. A. M. and E. P. F. C. contributed equally to this work.

b

M. Y. C. and S. N. T. contributed equally to this work.

Supplementary data