Reported associations of male circumcision (MC) with human papillomavirus (HPV) infection in men have been inconsistent.
4,033 healthy men were examined every six months for a median of 17.5 months. In each study visit, exfoliated cell specimens from the coronal sulcus/glans penis, penile shaft, and scrotum were collected and combined into one sample per person for HPV DNA detection. Samples were tested for 37 HPV types. Cox proportional hazards models were used to evaluate the association between MC and the incidence and clearance of HPV infections and specific genotypes.
The overall incidence of new HPV infections did not differ by MC status (for any HPV, adjusted hazard ratio (aHR) 1.08, 95% confidence interval (CI) 0.91-1.27). However, incidence was significantly lower among circumcised versus uncircumcised men for HPV types 58 (p = 0.01), 68 (p < 0.001), 42 (p = 0.01), 61 (p < 0.001), 71 (p < 0.001), 81 (p = 0.04), and IS39 (p = 0.01), and higher for HPV types 39 (p = 0.01) and 51 (p = 0.02). Despite the lack of an overall association in the risk of HPV clearance by MC (for any HPV, aHR 0.95, 95% CI 0.88-1.02), median times to clearance were significantly shorter among circumcised than uncircumcised men for HPV types 33 (p = 0.02) and 64 (p = 0.04), and longer for HPV types 6 (p < 0.001), 16 (p < 0.001), and 51 (p = 0.02).
MC is not associated with the incidence and clearance of genital HPV detection, except for certain HPV types. The use of a single combined sample from the penis and scrotum for HPV DNA detection likely limited our ability to identify a true effect of MC at the distal penis.
Keywords:Male circumcision; Genital; HPV; Incidence; Clearance
The majority of sexually active men and women will acquire genital human papillomavirus (HPV) infection at some point during their lifetime. In 2008, an estimated 610,000 cancers in men and women were attributed to HPV infection worldwide . The natural history of cervical HPV infection is well characterized; however, little is known about genital HPV infection in men. A recent randomized controlled trial (RCT) of male circumcision (MC) in Uganda (Africa) showed that MC reduced the incidence of high-risk HPV (HR-HPV) infections and increased clearance of HR-HPV infections at the coronal sulcus . However, findings regarding the role of MC in the incidence and clearance of genital HPV infections have not been consistent across studies [3-5].
Our previously published data regarding the prevalence of genital HPV among all men enrolled in the HPV Infection in Men (HIM) Study showed no overall association between MC and genital HPV infections, except for certain HPV types . The purpose of this study was to determine whether MC affects the incidence and clearance of genital HPV infections in a large, multinational cohort study of healthy men in Brazil, Mexico, and the United States (USA).
From June 2005 through September 2009, healthy men were enrolled in an ongoing longitudinal study of HPV infection in men, the HIM cohort study. Men were recruited from the general population, universities, and organized health-care systems. Details of the study have been previously described . In brief, men were included if they 1) were aged 18–70 years, 2) were residents of southern Florida, USA, Cuernavaca, Mexico, or São Paulo, Brazil, 3) had no prior diagnosis of penile or anal cancer, 4) had no prior diagnosis of genital or anal warts, 5) had no current diagnosis or symptoms of sexually transmitted infections (STIs), 6) had not participated in an HPV vaccine study, 7) reported no previous diagnosis of HIV, 8) had not been imprisoned, homeless, or received drug treatment during the past six months, 9) were willing to commit to 10 scheduled visits every six months, and 10) had no plans to relocate within the next four years. Men who were eligible to participate reviewed and signed a written informed consent form. Before study initiation, the Human Subjects Committees of the University of South Florida, the Centro de Referencia e Tratamento de Doencas Sexualmente Transmissiveis e AIDS, Brazil, and The National Institute of Public Health of Mexico approved the research protocol.
Consenting participants completed a pre-enrollment visit, an enrollment visit scheduled approximately two weeks later, and eight additional visits after enrollment that occurred every six months over a period of four years. At each study visit, participants completed an 88-item computer-assisted self-interview (CASI) which collected information regarding sociodemographic characteristics, tobacco consumption, and sexual behavior. The primary languages spoken by participants from Brazil, Mexico, and the United States were Portuguese, Spanish, and English, respectively, and the CASI was available in these languages. After the interview, men underwent a clinical examination at which time circumcision status was assessed. Participants with full or partial circumcision were considered circumcised.
Penile and scrotal sampling
Sampling techniques have been described in detail previously . Briefly, exfoliated epithelial cells from the coronal sulcus/glans penis, penile shaft, and scrotum were collected using three different saline pre-wetted Dacron swabs and combined into one sample before DNA extraction. Among uncircumcised men, the foreskin was sampled at the time of collecting the coronal sulcus/glans penis specimen. All HPV samples were stored at -80°C prior to polymerase chain reaction (PCR) analysis and genotyping. These sampling procedures were standardized across the three countries, with one clinical site per country. Staff at each site collecting specimens and conducting physical exams were either medical doctors or nurse practitioners.
The detailed protocol for HPV analysis has been previously described . Briefly, DNA was extracted using the Media Kit (QIAGen, Valencia, CA, USA), according to the manufacturer’s instructions. HPV testing was undertaken by use of PCR for amplification of a fragment of the HPV L1 gene. Specimens were tested for the presence of HPV with the PGMY09/11 L1 consensus primer system . HPV genotyping was conducted on all specimens, regardless of the HPV PCR result, using the Linear Array method (Roche Molecular Diagnostics, Alameda, CA, USA)  to detect 37 HPV types. The 13 HPV types classified as oncogenic included: 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59 and 68. Non-oncogenic HPV types: 6, 11, 26, 40, 42, 53, 54, 55, 61, 62, 64, 66, 67, 69, 70, 71, 72, 73, 81, 82, 83, 84, IS39, and 89 (CP6108) . Samples were considered valid if they were positive for β-globin and/or any HPV DNA genotype.
Differences in the distribution of sociodemographic characteristics, tobacco consumption, and sexual behavior characteristics by circumcision status were compared using the Pearson’s chi-square test. Differences in the median time to follow-up by MC status were compared using a nonparametric equality of medians test.
The classification of “any HPV” was defined as a positive test result for at least one of the 37 HPV genotypes included in the Linear Array test. HPV detection of single or multiple oncogenic HPV types was classified as “oncogenic HPV”. Similarly, HPV detection of single or multiple non-oncogenic HPV types was classified as “non-oncogenic HPV”.
HPV incidence by circumcision status was estimated according to different classifications of HPV types: any HPV, oncogenic HPV, non-oncogenic HPV, and for each specific HPV type. For each of the above analyses, only participants free of the relevant HPV type at enrollment were included. Men with concomitant incident HPV infections could contribute to several HPV classifications. Time to newly acquired HPV was estimated using the time from the date of study entry to the date of the first detection of HPV DNA, assuming a new infection arose at the date of detection. The exact 95% confidence intervals (CIs) for incidence estimates were based on the number of events modeled as a Poisson variable for the total person-months. The unit of observation was the individual participant, and each man with an incident HPV infection was counted only once during follow-up using only data from the first occurrence. Cumulative incidence of any HPV, HPV 16, oncogenic HPV, and non-oncogenic HPV types was estimated for circumcised and uncircumcised men using the Kaplan-Meier method.
Cox proportional hazards models were used to assess the association between HPV incidence and circumcision status . Hazards ratios (HRs) and 95% CIs were used as measures of association for comparison of circumcised versus uncircumcised men. The proportional hazards assumption for the Cox models was tested , and no violations were found, except in the non-oncogenic HPV model. The non-oncogenic HPV incidence model was stratified according to age to reduce violations of the proportional hazards assumption.
HPV clearance was defined as a participant testing HPV negative at two subsequent consecutive visits following a positive HPV test result, excluding those testing HPV positive for the first time at a participant’s final visit. Time to HPV clearance was estimated by using the time at which the participant first tested positive to the date of the first negative test. Median time to clearance was estimated among incident infections using the Kaplan-Meier method for any HPV, oncogenic HPV, non-oncogenic HPV, and for each specific HPV type, according to circumcision status. Analyses were performed for each individual HPV type. Men whose HPV infections did not clear were censored in the analysis.
Cox proportional hazards models with the robust covariance matrix estimator to account for within-subject correlation were used to assess the association between HPV clearance and circumcision status . Men with HPV infections, regardless of baseline HPV status, were included in the models. HRs and 95% CIs were used as measures of association for the comparison of circumcised versus uncircumcised men. The proportional hazards assumption for the Cox models was tested , and HPV clearance models were stratified according to age and country to reduce violations of the proportional hazards assumption.
An HPV infection was considered persistent if a man was HPV DNA-positive at two or more consecutive visits with the same specific HPV type, and an HPV infection was considered transient if a man was positive only once. When a participant missed a study visit, the results from the next visit were used.
The same variable selection procedure was used to evaluate factors associated with incidence and clearance of HPV infections. Factors that had a p-value < 0.10 were considered covariates. Backwards selection methods, with a significance threshold of 0.05, were used to identify covariates for inclusion in the final multivariable model. Candidate variables included education, marital status, smoking status, lifetime number of female sexual partners, lifetime number of male anal sex partners, number of female sexual partners in the past 3–6 months, number of male anal sex partners in the past 3 months, and six-month visit compliance status (i.e. whether the elapsed time between follow-up visits was longer than 6.5 months). Country (USA, Brazil, and Mexico) and age (categorical) were included in all models as study design factors. In addition, HPV status at baseline was a candidate variable in the HPV clearance models.
Statistical analyses were conducted using R 2.13.0 (R Development Core Team) and SAS software, version 9.2 (SAS Institute Inc., Cary, NC). Statistical tests were two-sided, with a significance threshold of 0.05.
Of the 4,074 initial HIM Study participants, 4,033 contributed valid samples for HPV DNA detection and were included in the analyses. Thus, 19 circumcised and 22 uncircumcised men (chi-square p-value = 0.03) contributed inadequate samples, as determined by both lack of β-globin detection and absence of any HPV genotype. The median duration of follow-up was 17.5 months (interquartile range [IQR], 6.9 – 31.0 months) for this analysis. There were no significant differences in the median duration of follow-up for circumcised (17.9 months) versus uncircumcised men (17.1 months) (p = 0.5 for nonparametric equality of medians test). More than 60% of men contributed at least 4 visits (IQR, 3 – 6 visits). The majority of men were uncircumcised (63.6%) (Table 1). Sociodemographic and behavioral characteristics of men varied by circumcision status. Circumcised men were generally younger (mean age: 31 vs. 33 years) and were residents of the USA (72.7%). Statistically significant differences by circumcision status were observed for marital status, sexual orientation, lifetime number of female sexual partners, number of female sexual partners in the past 3–6 months, lifetime number of male sexual partners, number of male anal sexual partners in past 3 months, self-reported diagnosis of STIs, and six-month visit compliance status (Table 1).
Table 1. Baseline characteristics of participants in the study cohort
HPV incidence by MC status
Kaplan-Meier analyses showed that the overall incidence of any HPV did not differ by circumcision status (P = 0.287, log-rank test). Similarly, incidence of HPV 16, oncogenic HPV, and non-oncogenic HPV did not differ by circumcision status (Figure 1). Estimates for incidence per 100 person-years of any HPV, oncogenic HPV, and non-oncogenic HPV were 50.5, 28.4, and 39.8, respectively, among uncircumcised men, and 45.6, 28.7, and 34.0, respectively, among circumcised men (Table 2). However, the incidence of oncogenic HPV types 39, 51, 58, and 68, and of non-oncogenic HPV types 42, 61, 71, 81 and IS39, did differ between circumcised and uncircumcised participants. Incidence was significantly higher among circumcised men than uncircumcised men for HPV types 39 (HR 1.47, 95% CI 1.11-1.94) and 51 (HR 1.30, 95% CI 1.04-1.61). In contrast, incidence was significantly lower among circumcised men than uncircumcised men for HPV types 58 (HR 0.61, 95% CI 0.43-0.88), 68 (HR 0.56, 95% CI 0.40-0.78), 42 (HR 0.59, 95% CI 0.38-0.90), 61 (HR 0.59, 95% CI 0.44-0.78), 71 (HR 0.37, 95% CI 0.21-0.65), 81 (HR 0.70, 95% CI 0.50-0.98), and IS39 (HR 0.41, 95% CI 0.20-0.86). The 12-month cumulative incidence of any HPV, oncogenic HPV, and non-oncogenic HPV among uncircumcised men was 38.7%, 23.7%, and 32.5%, respectively, compared with 35.2%, 23.3%, and 28.9%, respectively, among circumcised men.
Figure 1. Kaplan-Meier estimates of the cumulative incidence of human papillomavirus (HPV) infections. Kaplan-Meier estimates of the cumulative incidence of human papillomavirus (HPV) infections for any HPV, oncogenic HPV, non-oncogenic HPV, and HPV 16 by male circumcision status: A) Incidence of any HPV; B) Incidence of oncogenic HPV; C) Incidence of non-oncogenic HPV; D) Incidence of HPV 16.
Table 2. Incidence of HPV by male circumcision status
In univariate and multivariate analyses, no differences in the risk of HPV incidence by circumcision status were observed (Table 3). A total of 359 circumcised men (48.2%) had an incident infection, compared with 555 uncircumcised men (47.3%). For any HPV, no differences in incidence were found by circumcision status in univariate and multivariate analyses (adjusted hazard ratio [aHR] 1.08, 95% CI 0.91-1.27). Similarly, no differences in incidence were found by circumcision status for oncogenic HPV (aHR 1.11, 95% CI 0.94-1.31) and non-oncogenic HPV (aHR 1.11, 95% CI 0.94-1.30). The magnitude of the associations remained the same after adjustment for potential confounders, including sociodemographic characteristics and sexual behavior, or when changing the unit of analysis from men to infection (data not shown). Results were consistent when analyses were stratified by country of residence (data not shown).
Table 3. Univariate and multivariate hazard ratios for the association between male circumcision and incidence of genital HPV infection
HPV clearance by MC status
Kaplan-Meier analyses showed that the median time to clearance of any HPV infection was significantly longer among circumcised men than uncircumcised men (P < 0.0001, log-rank test) (Figure 2). Similarly, median time to clearance of oncogenic HPV and non-oncogenic HPV types was significantly shorter among uncircumcised men than circumcised men. Median time to clearance was also significantly shorter among uncircumcised men than circumcised men for HPV types 6 (HR 0.60, 95% CI 0.44-0.81), 16 (HR 0.56, 95% CI 0.42-0.75), and 51 (HR 0.72, 95% CI 0.54-0.95) (Table 4). Correspondingly, the probability of clearing an HPV infection was significantly lower among circumcised men compared to uncircumcised men. However, median time to clearance was significantly shorter among circumcised men compared to uncircumcised men for HPV 33 (HR 2.54, 95% CI 1.12-5.73). Although the median time to clearance for HPV 64 was significantly shorter among circumcised versus uncircumcised men, the corresponding HR did not reach statistical significance (HR 6.06, 95% CI 0.84-43.74).
Figure 2. Kaplan-Meier estimates of time to clearance of human papillomavirus (HPV) incident infections. Kaplan-Meier estimates of time to clearance of human papillomavirus (HPV) incident infections for any HPV, oncogenic HPV, non-oncogenic HPV, and HPV 16 by male circumcision status: A) Clearance of any HPV; B) Clearance of oncogenic HPV; C) Clearance of non-oncogenic HPV; D) Clearance of HPV 16.
Table 4. Clearance of HPV by male circumcision status
Univariate and multivariate analyses for clearance are presented in Table 5. HPV clearance was significantly decreased for any HPV infection among circumcised men compared to uncircumcised men in the univariate analysis (HR 0.92, 95% CI 0.86-0.99). However, after adjustment, this difference was no longer significant (aHR 0.95, 95% CI 0.88-1.02). Similarly, there were no differences in the risk of clearance for oncogenic HPV (aHR 0.90, 95% CI 0.81-1.00) or non-oncogenic HPV (aHR 0.98, 95% CI 0.89-1.07) by circumcision status. The magnitude of the associations with MC did not change substantially for transient (Table 6) versus persistent (Table 7) HPV infections. Moreover, similar results were observed when analyses were stratified by country of residence (data not shown).
Table 5. Univariate and multivariate hazard ratios for the association between male circumcision and clearance of genital HPV infection
Table 6. Univariate and multivariate hazard ratios for the association between male circumcision and clearance of transient genital HPV infection
Table 7. Univariate and multivariate hazard ratios for the association between male circumcision and clearance of persistent genital HPV infection
MC was not associated with an overall reduction in the incidence of genital HPV detection in men. Other longitudinal studies also found no differences in the risk of HPV incidence by MC [3-5]. Lu et al.  and Hernandez et al.  used combined samples from the coronal sulcus, glans penis, shaft of the penis, and scrotum for HPV detection and found no differences in HPV acquisition by MC. Although VanBuskirk et al.  observed site-specific differences in HPV incidence by MC, they found no differences in overall HPV incidence by MC at any site (penile shaft, glans/corona, and scrotum). Our findings are not consistent with those from a RCT conducted in Uganda, which found that MC reduced the incidence of HR-HPV infections (rate ratio [RR] 0.67, 95% CI 0.50-0.90) . However, this RCT was not comparable to our study, as they reported associations with specimens collected from only the coronal sulcus .
Despite the lack of evidence for an overall association, we did find higher HPV incidences rates in uncircumcised men than in circumcised men for most genotypes. These differences were statistically significant for HPV types 58, 68, 42, 61, 71, 81, and IS39. In contrast, the incidence of oncogenic HPV types 39 and 51 were significantly higher among circumcised than uncircumcised men. Consistent with our study, VanBuskirk et al.  reported higher HPV incidence among uncircumcised than circumcised men for >50% of the 21 specific HPV types assessed. Moreover, the RCT in Uganda found that HPV incidence was higher among uncircumcised men than among circumcised men for all HR-HPV genotypes, being statistically significant for HR-HPV types 18 and 33 .
No differences in overall HPV clearance by MC were observed in our study. This finding is in agreement with findings from other longitudinal studies [4,5]. VanBuskirk et al.  found that MC had no effect on the likelihood of detecting a persistent versus transient HPV infection over an 8-month period. Hernandez et al.  found that the duration of HPV infection did not vary by MC for the penile shaft, scrotum, or for all genital sites combined. However, the median duration of HPV infection at the glans/coronal sulcus was significantly longer in uncircumcised men that in circumcised men . Similarly, MC increased clearance rates of HR-HPV infections at the coronal sulcus in the Uganda RCT . In the study by Lu et al. , circumcised men were more likely to clear infection with any HPV and oncogenic HPV types using a combined HPV sample from the coronal sulcus, glans penis, shaft, and scrotum for HPV detection.
Differences in HPV clearance by MC for certain HPV types were observed in the present study. Median times to clearance for HPV types 33 and 64 were significantly shorter among circumcised men compared to uncircumcised men. The RCT in Uganda reported higher rates of clearance among circumcised men than uncircumcised men at the coronal sulcus for most HR-HPV types, but the differences were only statistically significant for HR-HPV types 39, 51, and 58 . In contrast, we found that the median times to clearance for HPV types 16, 51, and 6 were significantly longer among circumcised men than uncircumcised men.
Even though our study was not designed to assess the effects of MC and HPV by anatomic genital subsite, there is evidence from the literature that the association between MC and HPV DNA detection varies according to the subsite sampled. Four out of five studies show that the effect of MC in reducing HPV infection is stronger at the glans/corona than at other, more distant, subsites [14-18]. Weaver et al. , found that HPV prevalence was significantly lower among circumcised men (17%) compared to uncircumcised men (32%) on the glans and foreskin. Interestingly, and consistent with our findings, when combining results from all subsites, overall HPV prevalence was similar between circumcised (31%) and uncircumcised (29%) men. Consistent with these results, a study by Nielson et al.  reported that HPV prevalence at the glans penis/coronal sulcus was significantly lower among circumcised men (29.8%) compared to uncircumcised men (35.2%). However, HPV prevalence combining the glans, shaft, and scrotum was similar in circumcised (48.3%) and uncircumcised men (44.8%). In another US study, Hernandez et al.  reported a significantly lower HPV prevalence at the glans penis/coronal sulcus among circumcised men (29%) compared to uncircumcised men (46%). However, HPV prevalence at the external genitalia was similar among circumcised (63%) versus uncircumcised (71%) men. At year 1 in the RCT in Uganda , HR-HPV types were more frequently detected on the coronal sulcus than on the penile shaft among both circumcised and uncircumcised men. In addition, HR-HPV prevalence at year 1 on the penile coronal sulcus was significantly lower among circumcised men (21.5%) as compared to uncircumcised men (36.3%). Finally, a study by Partridge et al. , conducted among male university students in the USA, found that circumcision status was not associated with incident HPV infection for any of the individual subsites or for all sites combined.
Further supporting the evidence for an effect of MC on HPV infection on the proximal region of the penis, three additional studies reported specific inverse associations between MC and HPV infection in the glans, coronal sulcus, and the distal urethra , the glans and coronal sulcus , or the urethra .
In addition to the subsites included for HPV DNA detection, several other factors may explain the contrasting effects observed across studies, such as differences in participants’ ages, the proportion of circumcised men, the number of HPV genotypes included in testing, the time intervals between visits, and the age at which circumcision was performed.
To our knowledge, this is the first international study reporting HPV incidence and clearance by specific HPV types in circumcised and uncircumcised men. With a sample size of 4,033 men, it is also the largest study exploring these associations. A key strength of the present study is its longitudinal design, which allowed repeated measures of genital HPV DNA status for each participant, with a median follow-up of 17.5 months and scheduled visits every six months.
Our study, despite being large, multinational, and prospective, is not free of some limitations, including: (1) the potential for misclassification of HPV infection, as we combined samples from different anatomic subsites; (2) the potential for misclassification of prevalent infections as incident infections, which may inflate the estimates of cumulative incidence; (3) the possibility of not distinguishing HPV clearance from failure to detect a true HPV infection, even though two HPV negative results following one positive result were required to define HPV clearance; (4) potentially limited generalizability, as only men willing to comply with multiple clinical visits over four years were included in the study; and (5) the potential for confounding, as there may have been relevant factors related to religious or cultural practices that may also be associated with both MC status and HPV incidence or clearance that were not taken into account.
Overall, our study found that circumcision status was not associated with incidence and clearance of genital HPV detection, except for among specific HPV types. However, the use in this study of a single combined sample from the penis and scrotum likely limited our ability to identify a true effect at the distal penis. Further research is needed to determine how male circumcision impacts HPV incidence and clearance by specific anatomic subsite and HPV types.
HPV: Genital human papillomavirus; MC: Male circumcision; RCT: Randomized controlled trials; HR-HPV: High-risk human papillomavirus; HIM Study: HPV Infection in men study; STIs: Sexually transmitted infections; CASI: Computer-assisted self-interview; PCR: Polymerase chain reaction; CIs: Confidence intervals; HRs: Hazard ratios; IQR: Interquartile range; aHR: Adjusted hazard ratio; RR: Rate ratio.
GA: Travel grants to conferences/meetings are occasionally granted by either GlaxoSmithKline, Merck, Sanofi Pasteur MSD, Roche or Qiagen.
GA, XC, and FXB: My research unit is involved in HPV vaccine trials organized by GlaxoSmithKline, Merck and Sanofi Pasteur MSD.
XC: Travel grants to scientific meetings and honorarium for consultancy are occasionally granted by either GlaxoSmithKline, Merck, Sanofi Pasteur MSD.
ARG: has received research grant support from Merck & Co., Inc. and GSK.
ARG and LLV: are on the speakers’ bureau for Merck and are members of its advisory board.
AGN: has received research funding from Merck & Co.
FXB: Travel grants to conferences/symposia/meetings and honorarium are occasionally granted by either GlaxoSmithKline, Merck, Sanofi Pasteur MSD, Roche or Qiagen.
All other authors declare that they have no competing interests.
ARG, LLV, and EL-P were the principal investigators, conceived the study, wrote the protocols, assured funding, identified clinical investigators and study personnel, and made critical revision of the manuscript. GA was responsible for conduct analysis, interpretation of data, and drafting and revision of the manuscript. WF, H-YL, and MRP made the statistical analyses, and made substantial comments to the manuscript. MA trained and supervised study clinical staff, and made substantial contributions to the manuscript. JS and MQ coordinated the field work in terms of study implementation and data collection, and made substantial comments to the manuscript. AGN, BL, XC, and FXB participated in the interpretation of data, and critical revision of the manuscript. All authors read and approved the final manuscript.
The authors would like to thank the HIM Study Teams in the USA (Moffitt Cancer Center, Tampa, FL: AR Giuliano, JL Messina, C Gage, K Eyring, K Kennedy, K Isaacs, A Bobanic, BA Sirak, MT O’Keefe, DJ Ingles, MB Schabath, CM Pierce Campbell; Johns Hopkins, Baltimore, MD: RP Viscidi), Brazil (Centro de Referência e Treinamento em DST/AIDS, Fundação Faculdade de Medicina Instituto do Câncer do Estado de São Paulo, and Ludwig Institute for Cancer Research, São Paulo: LL Villa, ML Baggio, L Galan, RJ Carvalho da Silva, L Sichero, E Gomes, E Brito, F Cernicchiaro, R Cintra, R Cunha, R Matsuo, V Souza, B Fietzek, R Hessel, V Relvas, F Silva, J Antunes, G Ribeiro, R Bocalon, R Otero, R Terreri, S Araujo, M Ishibashi, CRT-DST/AIDS nursing team), and Mexico (Instituto Mexicano del Seguro Social and Instituto Nacional de Salud Pública, Cuernavaca: E Lazcano-Ponce, A Cruz Valdez, R Alvear Vásquez, O Rojas Juárez, R González Sosa, A Salgado Morales, A Rodríguez Galván, P Román Rodríguez, A Landa Vélez, M Zepeda Mendoza, G Díaz García, V Chávez Abarca, J Ruiz Sotelo, A Gutiérrez Luna, M Hernández Nevárez, G Sánchez Martínez, A Ortiz Rojas, C Barrera Flores). The authors would also like to thank the men who provided personal information and biological samples for the study. Finally, we would like to thank Donna J. Ingles for her work copyediting the manuscript.
Gray RH, Serwadda D, Kong X, Makumbi F, Kigozi G, Gravitt PE, Watya S, Nalugoda F, Ssempijja V, Tobian AA, Kiwanuka N, Moulton LH, Sewankambo NK, Reynolds SJ, Quinn TC, Iga B, Laeyendecker O, Oliver AE, Wawer MJ: Male circumcision decreases acquisition and increases clearance of high-risk human papillomavirus in HIV-negative men: a randomized trial in Rakai, Uganda.
Lu B, Wu Y, Nielson CM, Flores R, Abrahamsen M, Papenfuss M, Harris RB, Giuliano AR: Factors associated with acquisition and clearance of human papillomavirus infection in a cohort of US men: a prospective study.
Albero G, Villa LL, Lazcano-Ponce E, Fulp W, Papenfuss MR, Nyitray AG, Lu B, Castellsague X, Abrahamsen M, Smith D, Bosch FX, Salmeron J, Quiterio M, Giuliano AR: Male circumcision and prevalence of genital human papillomavirus infection in men: a multinational study.
Giuliano AR, Lee JH, Fulp W, Villa LL, Lazcano E, Papenfuss MR, Abrahamsen M, Salmeron J, Anic GM, Rollison DE, Smith D: Incidence and clearance of genital human papillomavirus infection in men (HIM): a cohort study.
Giuliano AR, Lazcano-Ponce E, Villa LL, Flores R, Salmeron J, Lee JH, Papenfuss MR, Abrahamsen M, Jolles E, Nielson CM, Baggio ML, Silva R, Quiterio M: The human papillomavirus infection in men study: human papillomavirus prevalence and type distribution among men residing in Brazil, Mexico, and the United States.
J Am Stat Assoc 1989, 84:1074-1078. Publisher Full Text
Biometrika 1993, 80:557-572. Publisher Full Text
Nielson CM, Schiaffino MK, Dunne EF, Salemi JL, Giuliano AR: Associations between male anogenital human papillomavirus infection and circumcision by anatomic site sampled and lifetime number of female sex partners.
Tobian AA, Kong X, Gravitt PE, Eaton KP, Kigozi G, Serwadda D, Oliver AE, Nalugoda F, Makumbi F, Chen MZ, Wawer MJ, Quinn TC, Gray RH: Male circumcision and anatomic sites of penile high-risk human papillomavirus in Rakai, Uganda.
Partridge JM, Hughes JP, Feng Q, Winer RL, Weaver BA, Xi LF, Stern ME, Lee SK, O'Reilly SF, Hawes SE, Kiviat NB, Koutsky LA: Genital human papillomavirus infection in men: incidence and risk factors in a cohort of university students.
Castellsague X, Bosch FX, Munoz N, Meijer CJ, Shah KV, de Sanjose S, Eluf-Neto J, Ngelangel CA, Chichareon S, Smith JS, Herrero R, Moreno V, Franceschi S: Male circumcision, penile human papillomavirus infection, and cervical cancer in female partners.
Baldwin SB, Wallace DR, Papenfuss MR, Abrahamsen M, Vaught LC, Giuliano AR: Condom use and other factors affecting penile human papillomavirus detection in men attending a sexually transmitted disease clinic.
Auvert B, Sobngwi-Tambekou J, Cutler E, Nieuwoudt M, Lissouba P, Puren A, Taljaard D: Effect of male circumcision on the prevalence of high-risk human papillomavirus in young men: results of a randomized controlled trial conducted in Orange Farm, South Africa.
The pre-publication history for this paper can be accessed here: