Abstract

The objective of our study was to describe oral and genital human papillomavirus (HPV) infection prevalence and concordance by sexual orientation among US men using a nationally representative sample. We conducted a retrospective cross-sectional analysis of the 2013-2016 National Health and Nutrition Examination Survey. The survey conducts a physical examination and collects oral rinse and genital swab specimens; demographic and health behaviors are self-reported. We used descriptive statistics and multivariate regression models to estimate HPV infection prevalence and the likelihood of HPV infection, respectively. All analyses were adjusted for National Health and Nutrition Examination Survey design and weights, and statistical significance was tested at a 2-sided P value of less than .05. Men who have sex with men had a statistically significantly higher prevalence of oral HPV (high-risk, 9-valent, 4-valent, and HPV 16 and 18), genital HPV (9-valent, 4-valent, and HPV 16 and 18), and concordant oral and genital HPV (high-risk and 9-valent) infections compared with heterosexual men. Improved HPV prevention among men is needed.

The incidence of human papillomavirus (HPV)–associated cancers, particularly oropharyngeal and anal cancers, is rising rapidly among US men (1-4). Gay, bisexual, and other men who have sex with men, collectively referred to as MSM, represent a large (nearly 7 million) and growing population at elevated risk of developing anogenital and oral HPV-associated cancers (5-7). The study of HPV infection occurring at oral and genital sites among MSM and the extent to which HPV infection risk is elevated compared with heterosexual men is crucial. Furthermore, given the likelihood of HPV transmission that may occur between oral and genital sites, the study of the risk of genital and oral HPV infection concordance (ie, presence of the same HPV genotype at oral and genital sites) among MSM compared with heterosexual men is also important, although it is yet to be quantified. This study aims to describe oral and genital HPV infection prevalence and concordance by sexual orientation among US men.

We conducted a retrospective cross-sectional analysis of the 2013-2016 National Health and Nutrition Examination Survey (NHANES), a nationally representative survey of noninstitutionalized individuals in the United States. A clustered multistage probability sample of participants is identified, and information on sociodemographic and health and sexual behaviors is collected by trained interviewers during the home interview or in mobile examination centers. The survey also conducts a physical examination and collects oral rinse and genital swab specimens in the mobile examination centers. The oral rinse sample (10 mL sample of mouthwash or saline) and self-collected external genital swabs are refrigerated at 4°C and shipped by the NHANES personnel to a laboratory to perform polymerase chain reaction assay for genotyping. Details regarding the survey are available on the NHANES website (8).

In the current study, men aged 18-59 years of age with nonmissing data on sexual orientation and oral and/or genital HPV tests were identified. Sexual orientation was self-reported (male sex partners contact and sexual identity). The HPV genotypes were identified from assay results and classified (as per prior studies) as type 16, type 18, 4-valent vaccine types (HPV 6, 11, 16, 18); 9-valent vaccine types (HPV 6, 11, 16, 18, 31, 33, 45, 52, 58); and any high-risk types (HPV 16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73, 82) (9,10). Concordant infection was defined as the detection of the same HPV genotype at oral and genital sites. The difference in infection prevalence was examined using the Rao-Scott χ2 test. Multivariable logistic regression models, adjusted for age, race, and ethnicity (Black, Hispanic, Other, and White); cigarette use; lifetime number of sex partners; history of sexually transmitted infections; age at first sex; and circumcision, estimated the odds ratio for HPV infection and concordant infection. Statistical significance was tested at a 2-sided P value less than .05. All analyses were performed with SAS software version 9.4 (SAS Institute, Cary, NC, USA) using SAS PROC SURVEY procedures to incorporate sampling weights and to adjust for the complex survey design. The study was deemed exempt from review by the institutional review board committee as it utilizes deidentified publicly available dataset.

The final sample included 3232 (156 MSM and 3076 heterosexual) men with information on oral HPV infection, 2954 (113 MSM and 2841 heterosexual) men with data on genital HPV infection, and 2883 (109 MSM and 2774 heterosexual) men with information on both infections (Supplementary Table 1, available online). Prevalence of oral HPV infection was statistically significantly higher among MSM compared with heterosexual men—any high-risk (12.0% vs 6.2%; P =.009), 9-valent types (9.5% vs 2.9%; P < .001), 4-valent (8.2% vs 1.6%; P < .001), and type 16 and 18 (6.9% vs 1.3%; P < .001) (Table 1). Prevalence of genital 9-valent (24.4% vs 13.3%; P = .006), 4-valent (19.0% vs 7.2%; P =.001), and type 16 and 18 (14.3% vs 5.0%; P =.01) HPV types was also statistically significantly higher among MSM compared with heterosexual men. Concordant oral and genital infection prevalence of any high-risk (5.5% vs. 1.4%; P < .001) and 9-valent (2.5% vs 0.4%; P =.07) HPV types was more common among MSM compared with heterosexual males. The likelihood of oral and genital any high-risk, 9-valent, 4-valent, and types 16 and 18 HPV infection was statistically significantly higher for MSM compared with heterosexual men (Supplementary Table 2, available online). The likelihood of high-risk oral and genital HPV concordance was statistically significantly greater for MSM compared with heterosexual men.

Table 1.

Prevalence and concordance of oral and genital HPV infection among US men by sexual orientation, NHANES 2013-2016

HPV typesOral HPV
Genital HPV
Oral- and genital HPV concordancea
MSMHeterosexualPbMSMHeterosexualPbMSMHeterosexualPb
(n = 156)(n = 3076)(n = 113)(n = 2841)(n = 109)(n = 2774)
HPV types 16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73, 82 (any high risk)
 Unweighted, with infection, No.19189.00935795.63644<.001
 Prevalence (95% CI), %12 (5.6 to 18.3)6.2 (5.0 to 7.4)30.1 (19.1 to 41.1)27.6 (24.9 to 30.3)5.5 (0 to 11.4)1.4 (0.9 to 1.9)
 Weighted, with infection/total population, No.c491 395/4 103 0154 401 882/71 090 232951 862/3 160 06618 111 189/65 614 189167 187/3 020 778657 335/64 089 450
HPV types 6, 11, 16, 18, 31, 33, 45, 52, 58 (9 valent)
 Unweighted, with infection, No.1186<.00127405.006416.007
 Prevalence (95% CI), %9.5 (3.7 to 15.3)2.9 (2.2 to 3.6)24.4 (14.5 to 34.3)13.3 (11.6 to 15.0)2.5 (0 to 5.7)0.4 (0.2 to 0.7)
 Weighted, with infection/total population, No.c389 535/4 103 0152 078 201/71 090 232770 065/3 160 0668 732 472/65 614 18974 921/3 020 778278 506/64 089 450
HPV types 6, 11, 16, 18 (4 valent)
 Unweighted, with infection, No.955<.00119214.001213.22
 Prevalence (95% CI), %8.2 (2.6 to 13.8)1.6 (1.2 to 2.0)19.0 (8.7 to 29.2)7.2 (6.1 to 8.4)0.7 (0 to 1.9)0.3 (0.1 to 0.4)
 Weighted, with infection/total population, No.c336 024/4 103 0151 126 245/71 090 232599 461/3 160 0664 753 939/65 614 18921 410/3 020 778169 349/64 089 450
HPV types 16, 18
 Unweighted, with infection, No.842<.00114139.01210.12
 Prevalence (95% CI), %6.9 (2.0 to 11.7)1.3 (0.9 to 1.7)14.3 (3.4 to 25.2)5.0 (4.0 to 6.0)0.7 (0 to 1.9)0.2 (0.1 to 0.3)
 Weighted, with infection/total population, No.c281 970/4 103 015928 793/71 090 232451 593/3 160 0663 295 664/65 614 18921 410/3 020 778134 225/64 089 450
HPV typesOral HPV
Genital HPV
Oral- and genital HPV concordancea
MSMHeterosexualPbMSMHeterosexualPbMSMHeterosexualPb
(n = 156)(n = 3076)(n = 113)(n = 2841)(n = 109)(n = 2774)
HPV types 16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73, 82 (any high risk)
 Unweighted, with infection, No.19189.00935795.63644<.001
 Prevalence (95% CI), %12 (5.6 to 18.3)6.2 (5.0 to 7.4)30.1 (19.1 to 41.1)27.6 (24.9 to 30.3)5.5 (0 to 11.4)1.4 (0.9 to 1.9)
 Weighted, with infection/total population, No.c491 395/4 103 0154 401 882/71 090 232951 862/3 160 06618 111 189/65 614 189167 187/3 020 778657 335/64 089 450
HPV types 6, 11, 16, 18, 31, 33, 45, 52, 58 (9 valent)
 Unweighted, with infection, No.1186<.00127405.006416.007
 Prevalence (95% CI), %9.5 (3.7 to 15.3)2.9 (2.2 to 3.6)24.4 (14.5 to 34.3)13.3 (11.6 to 15.0)2.5 (0 to 5.7)0.4 (0.2 to 0.7)
 Weighted, with infection/total population, No.c389 535/4 103 0152 078 201/71 090 232770 065/3 160 0668 732 472/65 614 18974 921/3 020 778278 506/64 089 450
HPV types 6, 11, 16, 18 (4 valent)
 Unweighted, with infection, No.955<.00119214.001213.22
 Prevalence (95% CI), %8.2 (2.6 to 13.8)1.6 (1.2 to 2.0)19.0 (8.7 to 29.2)7.2 (6.1 to 8.4)0.7 (0 to 1.9)0.3 (0.1 to 0.4)
 Weighted, with infection/total population, No.c336 024/4 103 0151 126 245/71 090 232599 461/3 160 0664 753 939/65 614 18921 410/3 020 778169 349/64 089 450
HPV types 16, 18
 Unweighted, with infection, No.842<.00114139.01210.12
 Prevalence (95% CI), %6.9 (2.0 to 11.7)1.3 (0.9 to 1.7)14.3 (3.4 to 25.2)5.0 (4.0 to 6.0)0.7 (0 to 1.9)0.2 (0.1 to 0.3)
 Weighted, with infection/total population, No.c281 970/4 103 015928 793/71 090 232451 593/3 160 0663 295 664/65 614 18921 410/3 020 778134 225/64 089 450
a

Concordant infection was defined as the detection of the same genotype of the HPV present at oral and genital sites. CI = confidence interval; HPV = human papillomavirus; MSM = men who have sex with men; NHANES = National Health and Nutrition Examination Survey.

b

P value for Rao-Scott χ2 test accounting for survey design and adjusted for weights.

c

Weighted No. and total population represent the number of men with infection and the total population in the United States, respectively, estimated using the NHANES sampling weights.

Table 1.

Prevalence and concordance of oral and genital HPV infection among US men by sexual orientation, NHANES 2013-2016

HPV typesOral HPV
Genital HPV
Oral- and genital HPV concordancea
MSMHeterosexualPbMSMHeterosexualPbMSMHeterosexualPb
(n = 156)(n = 3076)(n = 113)(n = 2841)(n = 109)(n = 2774)
HPV types 16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73, 82 (any high risk)
 Unweighted, with infection, No.19189.00935795.63644<.001
 Prevalence (95% CI), %12 (5.6 to 18.3)6.2 (5.0 to 7.4)30.1 (19.1 to 41.1)27.6 (24.9 to 30.3)5.5 (0 to 11.4)1.4 (0.9 to 1.9)
 Weighted, with infection/total population, No.c491 395/4 103 0154 401 882/71 090 232951 862/3 160 06618 111 189/65 614 189167 187/3 020 778657 335/64 089 450
HPV types 6, 11, 16, 18, 31, 33, 45, 52, 58 (9 valent)
 Unweighted, with infection, No.1186<.00127405.006416.007
 Prevalence (95% CI), %9.5 (3.7 to 15.3)2.9 (2.2 to 3.6)24.4 (14.5 to 34.3)13.3 (11.6 to 15.0)2.5 (0 to 5.7)0.4 (0.2 to 0.7)
 Weighted, with infection/total population, No.c389 535/4 103 0152 078 201/71 090 232770 065/3 160 0668 732 472/65 614 18974 921/3 020 778278 506/64 089 450
HPV types 6, 11, 16, 18 (4 valent)
 Unweighted, with infection, No.955<.00119214.001213.22
 Prevalence (95% CI), %8.2 (2.6 to 13.8)1.6 (1.2 to 2.0)19.0 (8.7 to 29.2)7.2 (6.1 to 8.4)0.7 (0 to 1.9)0.3 (0.1 to 0.4)
 Weighted, with infection/total population, No.c336 024/4 103 0151 126 245/71 090 232599 461/3 160 0664 753 939/65 614 18921 410/3 020 778169 349/64 089 450
HPV types 16, 18
 Unweighted, with infection, No.842<.00114139.01210.12
 Prevalence (95% CI), %6.9 (2.0 to 11.7)1.3 (0.9 to 1.7)14.3 (3.4 to 25.2)5.0 (4.0 to 6.0)0.7 (0 to 1.9)0.2 (0.1 to 0.3)
 Weighted, with infection/total population, No.c281 970/4 103 015928 793/71 090 232451 593/3 160 0663 295 664/65 614 18921 410/3 020 778134 225/64 089 450
HPV typesOral HPV
Genital HPV
Oral- and genital HPV concordancea
MSMHeterosexualPbMSMHeterosexualPbMSMHeterosexualPb
(n = 156)(n = 3076)(n = 113)(n = 2841)(n = 109)(n = 2774)
HPV types 16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73, 82 (any high risk)
 Unweighted, with infection, No.19189.00935795.63644<.001
 Prevalence (95% CI), %12 (5.6 to 18.3)6.2 (5.0 to 7.4)30.1 (19.1 to 41.1)27.6 (24.9 to 30.3)5.5 (0 to 11.4)1.4 (0.9 to 1.9)
 Weighted, with infection/total population, No.c491 395/4 103 0154 401 882/71 090 232951 862/3 160 06618 111 189/65 614 189167 187/3 020 778657 335/64 089 450
HPV types 6, 11, 16, 18, 31, 33, 45, 52, 58 (9 valent)
 Unweighted, with infection, No.1186<.00127405.006416.007
 Prevalence (95% CI), %9.5 (3.7 to 15.3)2.9 (2.2 to 3.6)24.4 (14.5 to 34.3)13.3 (11.6 to 15.0)2.5 (0 to 5.7)0.4 (0.2 to 0.7)
 Weighted, with infection/total population, No.c389 535/4 103 0152 078 201/71 090 232770 065/3 160 0668 732 472/65 614 18974 921/3 020 778278 506/64 089 450
HPV types 6, 11, 16, 18 (4 valent)
 Unweighted, with infection, No.955<.00119214.001213.22
 Prevalence (95% CI), %8.2 (2.6 to 13.8)1.6 (1.2 to 2.0)19.0 (8.7 to 29.2)7.2 (6.1 to 8.4)0.7 (0 to 1.9)0.3 (0.1 to 0.4)
 Weighted, with infection/total population, No.c336 024/4 103 0151 126 245/71 090 232599 461/3 160 0664 753 939/65 614 18921 410/3 020 778169 349/64 089 450
HPV types 16, 18
 Unweighted, with infection, No.842<.00114139.01210.12
 Prevalence (95% CI), %6.9 (2.0 to 11.7)1.3 (0.9 to 1.7)14.3 (3.4 to 25.2)5.0 (4.0 to 6.0)0.7 (0 to 1.9)0.2 (0.1 to 0.3)
 Weighted, with infection/total population, No.c281 970/4 103 015928 793/71 090 232451 593/3 160 0663 295 664/65 614 18921 410/3 020 778134 225/64 089 450
a

Concordant infection was defined as the detection of the same genotype of the HPV present at oral and genital sites. CI = confidence interval; HPV = human papillomavirus; MSM = men who have sex with men; NHANES = National Health and Nutrition Examination Survey.

b

P value for Rao-Scott χ2 test accounting for survey design and adjusted for weights.

c

Weighted No. and total population represent the number of men with infection and the total population in the United States, respectively, estimated using the NHANES sampling weights.

This study documents a higher prevalence of oral and genital HPV infection among MSM compared with heterosexual men. Our study further reports for the first time that concordant oral and genital HPV (high-risk and 9-valent) infection prevalence is also greater among MSM compared with heterosexual men. A previous study reported higher genital HPV prevalence among men having sex with men and women than among MSM or heterosexual men (11). A prior NHANES study also reported greater high-risk oral HPV infection prevalence among MSM compared with heterosexual men (12). Consistent with these studies, we found a higher prevalence of 9-valent, 4-valent, and type 16 and 18 oral and genital HPV infection and higher genotypes-specific concordance among MSM. Greater concordance of oral and genital HPV suggests that the transmission of genital and oral HPV may be more likely among MSM, possibly from autoinoculation by individuals through fingers or greater transmission risk through partners. Further research is needed to understand the factors that increase the susceptibility of HPV infection among MSM and the contribution of bidirectional transmission to cancer risk among MSM and heterosexual males. Careful consideration of the self-reported nature of the data and the limited sample for concordance analysis should be given when interpreting the findings of this study. The cross-sectional nature of our analysis also precludes causal inferences.

Our study has important cancer prevention and policy implications. The higher prevalence of oral and genital HPV among MSM underscores the importance of intensive efforts for HPV prevention. Unlike heterosexual males, the MSM population may not benefit from herd protection through female HPV vaccination, which increases the importance of achieving high vaccination coverage among males. Studies have documented that HPV vaccination is effective for protection against anogenital and oral HPV among MSM, particularly if initiated at younger ages (13,14). Unfortunately, HPV vaccination coverage among adolescent and young adult males is lower compared with females (15,16). Furthermore, only 32.8% of MSM aged 18-26 years were reportedly vaccinated in 2017 (17). Continued efforts are needed to improve vaccination coverage among adolescent and young adult boys and the currently unvaccinated MSM population to reduce their risk of developing HPV-associated cancers in future decades.

Funding

Research reported in this publication was supported by the National Institute on Minority Health and Health Disparities (Award Numbers K01MD016440) and the National Cancer Institute of the National Institutes of Health (Award Numbers R01CA256660 and R01CA232888).

Notes

Role of the funder: The funder had no role in the design of the study; collection, analysis, or interpretation of the data; the preparation of the manuscript or decision to submit it for publication.

Disclosures: Dr Sonawane has received consulting fees from Value Analytics Labs on unrelated projects. Dr Deshmukh has received consulting fees from Merck on unrelated projects. The remaining authors have no disclosures.

Author contributions: KS (Conceptualization, Formal Analysis, Funding acquisition, Resources, Supervision, Writing-original draft, Writing- reviewing & editing); SSS (Data Curation, Formal Analysis, Visualization); HD (Validation, Writing-original draft, Writing- reviewing & editing); RL (Formal Analysis, Visualization, Validation); AGN (Conceptualization, Writing-original draft, Writing- reviewing & editing); AAD (Conceptualization, Resources, Funding acquisition, Writing-original draft, Writing- reviewing & editing).

Acknowledgements: We acknowledge the National Center for Health Statistics (NCHS) of the Centers for Disease Control and Prevention (CDC) for making the National Health and Nutrition Examination Survey (NHANES) publicly available and men who participated in the survey.

Disclaimer: The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Data availability

Data is available publicly on the NHANES website.

References

1

Deshmukh
AA
,
Suk
R
,
Shiels
MS
, et al.
Recent trends in squamous cell carcinoma of the anus incidence and mortality in the United States, 2001-2015
.
J Natl Cancer Inst
.
2020
;
112
(
8
):
829
-
838
. doi:.

2

Damgacioglu
H
,
Sonawane
K
,
Zhu
Y
, et al.
Oropharyngeal cancer incidence and mortality trends in all 50 states in the US, 2001-2017
.
JAMA Otolaryngol Head Neck Surg
.
2022
;
148
(
2
):
155
-
165
. doi:.

3

Damgacioglu
H
,
Sonawane
K
,
Chhatwal
J
, et al.
Long-term impact of HPV vaccination and COVID-19 pandemic on oropharyngeal cancer incidence and burden among men in the USA: a modeling study
.
Lancet Reg Health Am
.
2022
;
8
:
100143
. doi:.

4

Damgacioglu
H
,
Lin
Y-Y
,
Ortiz
AP
, et al.
State variation in squamous cell carcinoma of the anus incidence and mortality, and association with HIV/AIDS and smoking in the United States
.
J Clin Oncol
.
2022
:JCO2201390. doi:.

5

Heck
JE
,
Berthiller
J
,
Vaccarella
S
, et al.
Sexual behaviours and the risk of head and neck cancers: a pooled analysis in the International Head and Neck Cancer Epidemiology (INHANCE) consortium
.
Int J Epidemiol
.
2010
;
39
(
1
):
166
-
181
. doi:.

6

Clifford
GM
,
Georges
D
,
Shiels
MS
, et al.
A meta-analysis of anal cancer incidence by risk group: toward a unified anal cancer risk scale
.
Int J Cancer
.
2021
;
148
(
1
):
38
-
47
. doi:.

7

Mosher
WD
,
Chandra
A
,
Jones
J.
Sexual behavior and selected health measures: men and women 15-44 years of age, United States, 2002
.
Adv Data
.
2005
;(
362
):
1
-
55
.

8

Centers for Disease Control and Prevention (CDC)
; National Center for Health Statistics (NCHS).
National Health and Nutrition Examination Survey Data
.
Hyattsville, MD
:
U.S. Department of Health and Human Services, Centers for Disease Control and Prevention
. https://wwwn.cdc.gov/nchs/nhanes/Default.aspx. Accessed March 2022.

9

Munoz
N
,
Bosch
FX
,
de Sanjose
S
, et al. ; for the
International Agency for Research on Cancer Multicenter Cervical Cancer Study Group
.
Epidemiologic classification of human papillomavirus types associated with cervical cancer
.
N Engl J Med
.
2003
;
348
(
6
):
518
-
527
. doi:.

10

Gillison
ML
,
Broutian
T
,
Pickard
RK
, et al.
Prevalence of oral HPV infection in the United States, 2009-2010
.
JAMA
.
2012
;
307
(
7
):
693
-
703
. doi:.

11

Nyitray
AG
,
da Silva
RJ
,
Baggio
ML
, et al.
The prevalence of genital HPV and factors associated with oncogenic HPV among men having sex with men and men having sex with women and men: the HIM study
.
Sex Transm Dis
.
2011
;
38
(
10
):
932
-
940
. doi:.

12

Sonawane
K
,
Suk
R
,
Chiao
EY
, et al.
Oral human papillomavirus infection: differences in prevalence between sexes and concordance with genital human papillomavirus infection, NHANES 2011 to 2014
.
Ann Intern Med
.
2017
;
167
(
10
):
714
-
724
. doi:.

13

Meites
E
,
Winer
RL
,
Newcomb
ME
, et al.
Vaccine effectiveness against prevalent anal and oral human papillomavirus infection among men who have sex with men–United States, 2016-2018
.
J Infect Dis
.
2020
;
222
(
12
):
2052
-
2060
. doi:.

14

Giuliano
AR
,
Palefsky
JM
,
Goldstone
S
, et al.
Efficacy of quadrivalent HPV vaccine against HPV infection and disease in males
.
N Engl J Med
.
2011
;
364
(
5
):
401
-
411
. doi:.

15

Sonawane
K
,
Zhu
Y
,
Montealegre
JR
, et al.
Parental intent to initiate and complete the human papillomavirus vaccine series in the USA: a nationwide, cross-sectional survey
.
Lancet Public Health
.
2020
;
5
(
9
):
e484
-
e492
. doi:.

16

Chen
MM
,
Mott
N
,
Clark
SJ
, et al.
HPV vaccination among young adults in the US
.
JAMA
.
2021
;
325
(
16
):
1673
-
1674
. doi:.

17

Loretan
C
,
Chamberlain
AT
,
Sanchez
T
,
Zlotorzynska
M
,
Jones
J.
Trends and characteristics associated with human papillomavirus vaccination uptake among men who have sex with men in the United States, 2014-2017
.
Sex Transm Dis
.
2019
;
46
(
7
):
465
-
473
. doi:.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.

Supplementary data