Male Circumcision and Human Papillomavirus ... - Semantic Scholar

6 downloads 241 Views 751KB Size Report
of human papillomavirus infection among fishermen along Lake. Victoria Shore in the ... Svare EI, Kjaer SK, Worm AM, Osterlind A, Meijer CJ, van den Brule. AJ.
MAJOR ARTICLE

Male Circumcision and Human Papillomavirus Infection in Men: A Systematic Review and MetaAnalysis Natasha Larke,1,2 Sara L. Thomas,2 Isabel dos Santos Silva,3 and Helen A. Weiss1,2 1Medical Research Council Tropical Epidemiology Group, 2Department of Infectious Disease Epidemiology, and 3Faculty of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, United Kingdom

Background. We systematically reviewed the evidence for an association between male circumcision and Human Papillomavirus (HPV) infection and genital warts in men. Methods. PubMed and Embase were searched to 15 September 2010. The measure of effect was the adjusted odds ratio (OR) or rate ratio (RR) when present and the crude estimate otherwise. Random effects meta-analyses were used to calculate summary measures of effect. Results. We identified 23 papers about the association between circumcision and HPV DNA. Circumcised men were less likely to have prevalent genital HPV infection than uncircumcised men (summary OR, 0.57, 95% confidence interval [CI], 0.45–0.71) with between-study heterogeneity (P-heterogeneity 5 0.006; I2 5 50.5%; 19 studies). Similar summary associations were seen in clinical and methodological subgroups. The effect of circumcision was stronger at the glans/corona (OR, 0.47; 95% CI, 0.37–0.60) and urethra (OR, 0.35; 95% CI, 0.12– 1.05) compared with sites more distal to the foreskin. There was weak evidence that circumcision was associated with decreased HPV incidence (summary RR, 0.75, 95% CI, 0.57–0.99; 3 studies) and increased HPV clearance (summary RR, 1.33; 95% CI, 0.89–1.98; 3 studies) but no evidence of an association with prevalent genital warts (OR, 0.93, 95% CI, 0.65–1.33; 15 studies). Conclusions. Several countries are expanding access to voluntary medical male circumcision to reduce HIV prevalence. This could provide additional benefit in reducing HPV prevalence.

Three randomized controlled trials (RCTs) have shown that circumcision reduces the risk of HIV acquisition in heterosexual men by about 60%, and WHO/UNAIDS recommend that circumcision be added to current HIV prevention strategies [1]. Consequently, several countries with high prevalence of HIV are expanding access to safe male circumcision for HIV prevention.

Received 8 January 2011; accepted 26 July 2011. Presented in part: International Union Against Sexually Transmitted Infections, Cape Town, South Africa, 10 November 2009. Correspondence: Natasha Larke, DPhil, MRC Tropical Epidemiology Group, Department of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, Keppel St, London, WC1E 7HT, United Kingdom (natasha.larke@ lshtm.ac.uk). The Journal of Infectious Diseases 2011;204:1375–90 Ó The Author 2011. Published by Oxford University Press on behalf of the Infectious Diseases Society of America. All rights reserved. For Permissions, please e-mail: [email protected] 0022-1899 (print)/1537-6613 (online)/2011/2049-0011$14.00 DOI: 10.1093/infdis/jir523

Circumcision also reduces the risk of other sexually transmitted infections (STIs) including trichomonas vaginalis, herpes simplex virus type 2, and genital ulcer disease [2–5]. Proposed biological mechanisms include reduction of pathogen entry through abrasions in the thinly keratinized inner mucosal surface of the foreskin [6] and removal of the moist environment under the foreskin, which may favor pathogen survival and replication [7]. Human papillomavirus (HPV) infection is the major cause of genital warts and a necessary cause for cervical cancer. In men, HPV infection is strongly associated with anal cancer and some penile cancers [8]. There is inconclusive evidence to date that circumcision protects against genital HPV infection. A systematic review of studies published to March 2006 found no evidence of an association between circumcision and genital HPV (OR, 1.20; 95% CI, 0.80–1.79) [9]. However, this review was criticized methodologically, and a reanalysis of the same studies found a strongly protective effect (OR, 0.56; 95% CI, 0.39–0.82) [10]. An updated Male Circumcision and HPV Infection

d

JID 2011:204 (1 November)

d

1375

meta-analysis using the same search strategy and including studies published to September 2007 found a similar effect (OR, 0.52; 95% CI, 0.33–0.82) [11]. The discrepancy between these findings is partly due to methodological challenges including the choice of anatomical sites for HPV sampling. For example, 1 study among HPVinfected men detected HPV more frequently at the penile shaft of circumcised than uncircumcised men [12]; it has been suggested that failure to sample the penile shaft results in underdetection of HPV in circumcised men compared with uncircumcised men [9]. However, other data suggest that observed differences in genital HPV prevalence between circumcised and uncircumcised men are not due to detection biases [10, 13]. Other methodological difficulties include the definition and methods used to ascertain the outcome. Previous reviews included circumcision as a keyword in the search strategy. This may preferentially identify studies that found an association with circumcision status. In this paper, we are reporting a systematic review and metaanalyses of the effect of circumcision on HPV and genital warts in men. This adds to previous reviews in several ways. We included several recent studies, including 3 RCTs; we also considered HPV incidence and clearance as outcomes; and we examined for the first time the association between circumcision and genital warts in detail. In addition, we used a comprehensive search strategy that did not include the keyword circumcision to avoid selection bias, and we undertook a priori defined subgroup analyses and meta-regression to allow in-depth analysis and to address methodological issues. METHODS Search Strategy

PubMed and Embase were searched up to 15 September 2010 (see Supplementary material for search terms). There were no language or publication date restrictions. All abstracts were reviewed independently by 2 authors (N. L. and H. W.), and papers likely to contain information on HPV risk factors were deemed potentially relevant. Full-text copies of these papers were obtained and reviewed. We also searched reference lists of relevant papers. We e-mailed authors when information was missing from published papers. Inclusion Criteria

Studies of HPV infection were restricted to those detecting HPV DNA using Hybrid Capture 2 or PCR because of the low sensitivity and specificity of other DNA detection methods and the low sensitivity of serological HPV tests [14]. Eligible samples included HPV DNA sampled through exfoliated cytology of the male anogenital sites as well as semen, urine, and swabs of lesions or biopsy material. Studies of genital warts were eligible if they detected genital lesions through either self-report or clinical 1376

d

JID 2011:204 (1 November)

d

Larke et al

exam. Studies among men who had sex with men were excluded because circumcision is only likely to have a beneficial effect on those who practice insertive anal intercourse. Data Extraction

A standardized prepiloted data extraction form in Microsoft ACCESS was used. Data were extracted independently by 2 authors (N. L. and H. W.), and any inconsistencies were discussed to reach consensus. Study populations described in more than 1 paper were included only once; we used data from the paper with most detailed information. When more than 1 outcome (HPV prevalence, HPV incidence, HPV clearance, or genital warts) was evaluated in 1 study population, these were classified as separate studies. Data Analysis

Meta-Analysis The association between circumcision and genital HPV prevalence at any anatomical site was estimated using ORs. When ORs were not presented, crude ORs and 95% CI were calculated from the raw data. The effects of circumcision on HPV incidence and clearance were estimated using ORs or rate ratios (RR). A best estimate of the effect for each study was selected to be the adjusted estimate when present and the crude estimate otherwise. Summary measures of effect were obtained from random effects meta-analysis. Meta-analyses of HPV incidence and clearance were limited to studies that presented the RR, because HPV infection was relatively common in all studies and the OR was likely to overestimate the RR (ie, lie further from unity). The effects of circumcision on prevalence of high- and low-risk HPV genotypes for cervical cancer were also evaluated separately. The effect of circumcision on HPV included studies in which the majority of men were HIV- uninfected. The one study among HIV-positive men is presented separately, because the effect of circumcision on HPV may differ between HIV-positive and -negative men. The effect of circumcision on HPV infection and genital warts may be biased if noncases are symptomatic for STIs because circumcision protects against some other STIs [2–5]. Therefore, study populations were categorized as high- or low-risk for other STIs (see Supplementary methods). Publication bias was assessed via funnel plots and the Egger funnel plot asymmetry test [15]. Data were analyzed using Stata 11.1 (StataCorp, Texas). Heterogeneity Between-study heterogeneity was evaluated using I2 and the P value for heterogeneity (Cochrane’s Q statistic) [16]. To evaluate potential sources of heterogeneity, we compared I2 values between clinical and methodological subgroups [16]. Summary estimates from subgroups were formally compared using metaregression. We compared subgroups according to (1) adjustment for covariates, (2) ascertainment of circumcision status

(clinical exam vs self-report), and (3) risk of STIs in the study population. For HPV studies, we also assessed (4) the range of HPV genotypes tested. For genital warts studies, we also assessed (5) current warts (vs history of warts) and (6) ascertainment of genital warts (clinical exam or cytology vs self-report) Anatomical Sampling of HPV Because anatomical distribution of HPV may vary by circumcision status [12, 17, 18], a sensitivity analysis was restricted to studies that measured HPV at a minimum at the following sites: glans/corona, penile shaft, and scrotum (sites with the highest prevalence). The effect of circumcision on HPV may also vary by anatomical site, representing a further potential source of heterogeneity. Therefore, we analyzed the effect of circumcision on prevalent HPV at different anatomical sites.

Results Results of Search Strategy

We identified 3366 papers from the database searches and obtained the full text for 306 based on abstract reviews. Of the 306 papers, 269 were excluded based on the full text (Figure 1). An additional 6 papers were identified from references and review

papers. Of the 43 papers for which data were extracted, 23 contained information on the association between circumcision and HPV DNA, and 16 contained information on genital warts. (One paper contained information on both outcomes.) A further 5 papers contained information only on penile cancer, the results of which are reported elsewhere. Description of Eligible HPV DNA Studies

The 23 papers (30 studies) evaluating the effect of circumcision on HPV DNA prevalence, incidence, and/or clearance included subgroups from 3 RCTs and 4 cohort studies (Tables 1 and 2). Twenty-seven studies (22 papers) were among predominantly HIV-negative populations, and 3 studies (1 paper) were among exclusively HIV-positive men [33]. Twenty studies provided information on the association of circumcision and HPV prevalence [4, 12, 17–32] (Table 1), 5 studies on HPV incidence [25, 34, 36, 37] and 5 studies on HPV clearance [25, 34–36] (Table 2). Study populations were from North America (10 studies), Central and South America (7 studies), Europe (5 studies), Africa (9 studies), and Asia (2 studies). Circumcision prevalence varied from 6% (the Netherlands) to more than 80% (the United States and South Korea). Circumcision status was

Figure 1. Flow chart of study selection for inclusion in the systematic review. Numbers in brackets represent number of studies. * Data extracted for a separate review (reported elsewhere).   Studies where majority of men are HIV negative. Male Circumcision and HPV Infection

d

JID 2011:204 (1 November)

d

1377

1378

Table 1. Summary of Studies of the Association of Male Circumcision and HPV Prevalence

d

JID 2011:204 (1 November)

First author

Study design

Country and year of study

Study population (high/low riska)

Study size (% HPV circumcised) prevalence

Sites sampled for HPV DNAb

Assessment of circumcision

Crude OR (95% CI)

Adjusted OR (95% CI)

HIV-negative populationsc Randomized controlled trials

d

Larke et al

Auvert [19]

RCT

d

South Africa, 2002–2004

1264 (51%)

19%e

Urethra

Clinical exam

0.54 (0.40–0.73)f

Tobian [4]

RCT

g

Uganda, 2003–2006 Men ages 15–49 y HIV negative at baseline (low risk)

520 (45%)

44%h

Glans/corona

Clinical exam

0.52 (0.36–0.76)f

Men ages 18–24 y HIV negative at baseline, HIV prevalence at follow-up 7.3% (low risk)

Observational studies Aynaud [20]

Cross-sectional

France, published in 2002

Partners of women with HPV-associated genital lesions (low risk)

111 (20%)

23%i

Semen

Clinical exam

0.95 (0.25–3.13)

Baldwin [21]

Cross-sectional

United States, 2000–2001

STI clinic attendees (high risk)

344 (67%)

27%j

Glans/corona

Clinical exam

0.35 (0.2120.57) 0.34 (0.20–0.57)k

Bleeker [22]

Cross-sectional

Netherlands, 2002

Dermatology clinic attendees (low risk)

83 (16%)

25%h

Foreskin and glans/corona

Clinical exam

0.87 (0.14–3.90)

Bleeker [22]

Cross-sectional

Netherlands, 1995–2002

Partners of women with CIN (low risk)

170 (6%)

59%h

Foreskin and glans/corona

Clinical exam

2.88 (0.55–28.5)

1139 (26%)

16%l

Urethra and glans/corona

Self-reported

0.24 (0.13–0.41) 0.37 (0.16–0.85)m

988 (40%)

59%h

Foreskin, glans/ corona, shaft and scrotum

Clinical exam

1.03 (0.72–1.47) 0.70 (0.52–0.94)n

University population, HIV negative (low risk)

254 (77%)

83%h

Glans/corona, shaft, scrotum, urine, and semeno

Clinical exam

0.74 (0.31–1.65) 0.49 (0.19–1.28)p

Castellsague [23] Cross-sectional

Spain, Colombia, Husbands of women Brazil, Thailand, from case control and Philippinnes, study of cervical 1985–1993 cancer (low risk) Brazil, Mexico, General, HIV negative and US, (low risk) 2005–2006

Giuliano [24]

Cross-sectional

Hernandez [17]

Cross-sectional

Lajous [25]

Baseline data Mexico, from cohort study 2000–2003

Military (high risk)

925 (10%)

42%j

Urethra, glans/ corona, shaft, and scrotum

Self- reported

0.47 (0.29–0.75)q 0.48 (0.30–0.77)r

Mu¨ller [26]

Cross-sectional

South Africa, 2006–2008

Sexual health clinic attendees, HIV prevalence 49.5% (high risk)

208 (26%)

77%h

Glans/corona, and shaft

Clinical exam

0.43 (0.21–0.85) 0.51 (0.21–1.25)s

Ng’ayo [27]

Cross-sectional

Kenya, published 2008

Men working in the fishing industry, HIV prevalence 25.6% (high risk)

250 (7%)

58%h

Glans/corona, shaft, scrotum, and perianal area

Probably clinician 0.56 (0.19–1.66)

US, 2004–2006

Table 1. Continued

First author Nielson [18]

Study design Cross-sectional

Country and year of study

Study population (high/low riska)

Study size (% HPV circumcised) prevalence

Sites sampled for HPV DNAb

Assessment of circumcision

Crude OR (95% CI)

United States, 2002–2006

General, HIV negative (low risk)

463 (84%)

51%h

Foreskin, urethra, glans/corona, shaft, scrotum, perianal area, anus, and semen

Clinical exam

0.99 (.58–1.68)

Adjusted OR (95% CI) 0.53 (.28–.99)t

Male Circumcision and HPV Infection

Ogilvie [28]

Cross-sectional

Canada, published 2009

STI clinic attendees (high risk)

262 (50%)

70%h

Foreskin, glans/ corona, shaft, and scrotum

Clinical exam

1.14 (.65–2.00)

Rombaldi [29]

Cross-sectional

Brazil, 2003–2004

Partners of women with CIN (low risk)

99 (10%)

55%u

Not reported

2.09 (.44–13.20)

Shin [30]

Cross-sectional

South Korea, 2002

University and college students (low risk)

368 (88%)

9%v

Foreskin, urethra, glans/corona, and shaft Foreskin, urethra, glans/corona, shaft, and scrotum

Self-reported

1.31 (.38–7.00)

1.80 (.40–8.20)w

Svare [31]

Cross-sectional

Denmark, 1993

STI clinic attendees (high risk)

198 (12%)

45%x

Glans/corona, shaft, Not reported scrotum, and perianal area

0.21 (.05–.68)

0.20 (.06–.60)y

Vaccarella [32]

Cross-sectional

Mexico, 2003–2004

Vasectomy clinics (low risk)

779 (32%)

9%h

Glans/corona, shaft, Clinical exam and scrotum

0.19 (.07–.44)

0.20 (.10–.40)w

Weaver [12]

Cross-sectional

United States, 2001–2002

University students (low risk)

317 (81%)

31%j,*

Foreskin, glans/ corona, shaft, scrotum, and urine

Clinical exam

1.05 (.53–2.15)

Uganda, 2003–2007

Men ages 15–49 y, HIV positive at baseline with CD4 cell counts 1350 cells/mL (low risk)

191 (45%)

64% 

Glans/corona

Clinical exam

0.49 (.26–.93)à

HIV-positive populationsc Serwadda [33]

RCT

d

JID 2011:204 (1 November)

Abbreviations: CD4 cells, T cells with CD4 receptor; CI, confidence interval; CIN, cervical intraepithelial neoplasia; OR, odds ratio; RCT, randomized controlled trial; STI, sexually transmitted infection. a

Study populations categorized as high or low risk for STIs (see Supplementary methods for details).

b

Exfoliated cells collected by study clinician/nurse unless otherwise stated. It was assumed that HPV was sampled by clinician when the person sampling was not explicitly stated in the article but a clinical exam was undertaken. All semen samples were self collected. c

HIV negative populations 5 Populations in which the majority of participants are HIV negative; HIV positive populations 5 Populations that are all HIV positive.

d

A subgroup (56%) of the randomized population was tested for HPV infection at 21-month visit.

e

Thirteen HPV genotypes high-risk for cervical cancer were evaluated using reverse line blot (Roche Diagnostics) (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68).

f

An adjusted prevalence risk ratio was presented in the paper; adjustment did not affect the magnitude or precision of the risk ratio compared to the crude estimate. In the South African trial, adjustment factors included ethnicity, age, education, lifetime number of sex partners, marital status, number of nonspousal partners in past 12 months, condom use in past 12 months, number of sex acts in past 12 months, HIV status, and in the Ugandan trial enrollment characteristics, rates of sex practices, and symptoms of STIs.

d

1379

g

A subgroup (50%) of the randomized population was tested for HPV infection at the 24-month visit.

h

Thirty-seven HPV genotypes evaluated using reverse line blot (Roche Diagnostics).

JID 2011:204 (1 November)

Larke et al

The effect of circumcision on the prevalence of low-risk HPV was also evaluated: crude OR, 0.29; 95% CI, .15–.57. à

Amongst those with all sites sufficient for detection of HPV.

Fourteen high-risk HPV genotypes evaluated using reverse line blot (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, and 68) (Roche Diagnostics).

Adjusted for age, lifetime number of sex partners, sex partners in last year, and genital warts.

To assess whether failure to sample sufficient sites may underestimate HPV prevalence, we compared prevalence at specific anatomic sites by circumcision status among the 3 studies with this information [12, 17, 18]. The anatomical distribution of HPV was roughly similar between circumcised and uncircumcised men in 2 studies [17, 18] (Figure 2). The shaft was the most common site for infection in both circumcised (40% and 50%) and uncircumcised men (41% and 60%) followed by the glans/corona, scrotum, or foreskin (when sampled). In the third study [12], the shaft was the most common site among circumcised men (27%), followed by the scrotum and glans/corona. HPV prevalence was considerably higher at the foreskin (30%) than other sites in uncircumcised men.

 

*

y

Adjusted for age and lifetime number of sex partners.

Six HPV genotypes detected using HPV genotype specific PCR (6, 11, 16, 18, 31, and 33). x

w

Twenty-five HPV genotypes detected using reverse line blot.

Six HPV genotypes evaluated using RFLP of PCR product (6, 11, 16, 40, 61, 84). d

v

u

Adjusted for patient group and HIV status.

Adjusted for date of analysis, smoking status, lifetime number of sex partners, and condom use in past 3 months. t

s

Age adjusted.

Adjusted for age, socioeconomic status, and lifetime number of sex partners. r

Adjusted for age, birthplace, ethnicity, education, lifetime number of sex partners, history of sex with men, age at first sexual encounter, condom use, history of genital warts, and history of cigarette smoking.

HPV was sampled at the foreskin, but this sample was not included in the analyses of HPV at ‘‘any site.’’

Anatomical Distribution of HPV

q

p

o

Adjusted for age, study location, education, age at first sex, lifetime number of sex partners, and frequency of genital washing after sex.

Adjusted for lifetime number of sex partners, sex in past 3 months, and number of sex partners in the past 3 months. n

m

At least 6 HPV genotypes (6, 11, 16, 18, 31, and 33) detected using oligonucleotide probes.

Twenty-seven HPV genotypes detected using reverse line blot (Roche Diagnostics).

Six HPV genotypes evaluated (6, 11, 16, 18, 33, and 42). Method not given.

Adjusted frequency of sex, genital warts, condom use in past 3 months, and having a regular partner. d

l

k

j

i

1380

evaluated by a clinician in 22 studies, self-reported in 5 studies, and not stated in 3 studies. All studies measured HPV DNA using PCR, and 23 measured more than 24 high- and low-risk HPV genotypes for cervical cancer.

Circumcision and HPV Prevalence

HPV prevalence was lower in circumcised than in uncircumcised men in 14 of the 19 studies among predominately HIV-negative men (Table 1). ORs ranged from 0.20 to 2.88. Meta-analysis indicated a strong protective effect overall (summary OR, 0.57; 95% CI, 0.45–0.71) but with substantial between-study heterogeneity (P-heterogeneity 5 .006; I2 5 50.5%; Figure 3A). There was no evidence of publication bias from the funnel plot (available from the author on request) and Egger test (P 5 0.56). Estimates for high-risk HPV genotypes (summary OR, 0.59; 95% CI, 0.49–0.70) and low-risk (summary OR, 0.56; 95% CI, 0.45–0.70) were similar to those for any HPV genotype (see Supplementary material). Heterogeneity

The strong protective effect of circumcision on HPV prevalence was seen consistently within the a priori defined subgroups (Figure 4). The adjusted estimates demonstrated a stronger effect (OR, 0.45; 95% CI, 0.33–0.61) compared with crude estimates (OR, 0.75; 95% CI, 0.55–1.02; P value for meta-regression 5 .03), although heterogeneity remained among adjusted estimates (P-heterogeneity 5 0.02; I2 5 54%). Similar summary estimates were observed among studies that ascertained circumcision status clinically and from self-report, and heterogeneity was observed within these 2 subgroups (I2 5 56% and 39%, respectively) (Figure 4). There was no evidence that the effect of circumcision varied according to the risk of STIs in the population (P 5 0.49) or the range of HPV genotypes examined (P 5 0.72); considerable heterogeneity remained within these subgroups.

Table 2. Summary of Studies on the Association of Male Circumcision and HPV Incidence and Clearance HPV incidence

First author

Study design, country, and year of study

Study population

Cumulative Assessment Frequency HPV Crude effect of of HPV Loss to Study size incidence estimate circumcision sampling follow-up (% circumcised) (12m) (95% CI)

HPV clearance Adjusted effect estimate (95% CI)

Study size Clearance of (% circumcised) infections

Crude effect estimate (95% CI)

Adjusted effect estimate (95% CI)

HIV-negative populationsa Clinical exam Baseline, 6 m, General, 12 m, and HIV-negative 24 mb,c married men

n/ad

RCT, Uganda, 2003–2006

Hernandez [35]

Clinician Cohort, United University students, States, HIV negative 2004–2006

Every 2 m (mean follow-up 431d)j

20%k

Lajous [25]

Cohort, Mexico, Military, Self report 2000–2003 HIV negative

Baseline and 12 m follow-upn

67%o

210 (17%)

21%p

OR, 1.02 (0.42–2.50)

Lu [36]

Cohort, United General, Clinical exam Every 6 m for States, HIV negative 18 mt 2003–2005

24%

285 (8%)

29%u

RR, 1.10 (0.50–2.30)

Male Circumcision and HPV Infection

Partridge [37] Cohort, United University Clinical exam Every 4 m for Not given States, students ages 3 y periodw 2003–2006 18–20 y, HIV negative

840 (53%)

24%e, f

Gray [34]

RR, 0.67 (0.50–0.90)g

645 (41%)

RR, 1.36 18.0 per (1.13–1.63) 10 pyr (over 24 m)

RR, 1.39 (1.17–1.64)i

357 (81%)

Not given

Not given

RR, 0.96 (0.71–1.32)m

OR, 1.12 (0.45–2.80)q

105 (11%)

64.8% (over 1y)

OR, 6.67 (0.1–N)r

OR, 8.33 (1.15–N)s

RR, 0.80 (0.40–1.90)v

285 (88%)

9.35 per 10 pyr RR, 2.7 cleared (1.3–5.7) infections (median follow-up 15.5 m)

RR, 3.1 (1.2–8.2)v

RR, 0.67 (0.50–0.91)h

Data not givenl

240 (77%)

35%e,

174 (47%)

50%e

x

RR, 1.10 (0.60–2.00)

n/a

n/a

n/a

n/a

HIV-negative populationsa Serwadda [33] RCT, Uganda, 2003–2007

Men, ages 15–49 y, HIV positive at baseline

Clinical exam Baseline and 24 m

n/ay

RR, 0.74 (0.57–1.01)

RR, 0.68 (0.44–1.04)z

174 (47%)

73.6% Not given (over 24 m)

RR, 1.09 (0.94–1.27)*

Abbreviations. CI, confidence interval; n/a, not applicable; OR, odds ratio; pyr 5 person years; RCT, randomized controlled trail; RR, rate ratio. a

HIV-negative populations 5 Populations where majority of participants are HIV negative; HIV-positive populations 5 populations that are all HIV positive.

b

All participants were assayed at baseline and 24 months. 330/441 (75%) of participants randomly selected from the intervention arm were assayed at 6 and 12 months, and 321/399 (80%) and 314/399 (77%) of participants randomly selected from the control arm were assayed at 6 and 12 months respectively.

d

c

HPV sampled from glans/corona.

JID 2011:204 (1 November)

d

Men with both baseline and 24-month HPV samples were randomly selected: 441/835 (52.8%) of HIV-negative married men with samples at both time points in the intervention arm; 399/803 (49.7%) of HIV-negative married men with participants at both time points in the control arm. e

Fourteen high-risk HPV genotypes evaluated using reverse line blot (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, and 68) (Roche Diagnostics).

f

Twelve-month cumulative incidence risk of high-risk HPV isotpyes measured among 644 randomly selected samples.

g

Results presented for high-risk HPV. Marginal evidence of an effect of circumcision on the incidence of low-risk HPV was also observed (crude RR, 0.84; 95% CI , 0.66–1.10).

h

Adjusted for age, education, condom use, alcohol consumption with sex, and number of sex partners.

i

Adjusted for age, education, number of sex partners, and condom use.

j

HPV sampled from glans/coronal sulcus, shaft, and scrotum.

k

d

l

1381

m

Proportion of men who enrolled but did not complete more than 1 visit. Authors report no difference in HPV acquisition between circumcised and uncircumcised men, but data were not shown. Adjusted for age, ethnicity, birthplace, education, lifetime number of partners, history of sex with men, condom use, and history of genital warts.

Men with both baseline and 24-month HPV samples were randomly selected: 21.7% of men in the intervention arm and 23.9% of men in the control arm.

To address possible underdetection of HPV, we first restricted analysis to the 10 studies that measured HPV at all of the glans/ corona, penile shaft and scrotum. The effect was the same as for all studies (OR, 0.58, 95% CI, 0.40–.82). However there as was still considerable between-study heterogeneity (P 5 0.003; I2 5 63.8%). The magnitude of the effect estimate varied by anatomical site of HPV sampling (Figure 4 and Supplementary material). There was some evidence of effect at the urethra (OR, 0.35; 95% CI, 0.12–1.05) and strong evidence of effect at the glans/corona (OR, 0.47, 95% CI, .37–.60). There was less evidence of an association at sites more distal to the foreskin (shaft, scrotum, and anus) or in semen. Notably, there was very little heterogeneity between estimates from the glans/corona, shaft, and semen.

JID 2011:204 (1 November)

d

Adjusted for correlation between multiple clearance observed within the same individual. *

Adjusted for marital status and number of sexual partners.

Circumcision and HPV Incidence

z

y

HPV sampled from foreskin, glans/corona, shaft, and scrotum.

Approximate 12-month cumulative incidence (taken from Figure 1 in article). The 24-month cumulative incidence was 62.5%. x

w

Thirty-seven HPV genotypes evaluated using reverse line blot (6, 11, 16, 18, 26, 31, 33, 35, 39, 40, 42, 45, 51–59, 61, 62, 64, 66–73, 81–84, CP6108, and IS39) (Roche Diagnostics).

Adjusted for cigarette smoking and lifetime number of sex partners. v

u

Adjusted for age, religion, intercourse with prostitute, anal intercourse with women, and baseline infections.

HPV sampled from glans/corona, shaft, and scrotum. t

s

Adjusted for age, socioeconomic status, lifetime number of sex partners, and anal intercourse with men.

Age adjusted. r

Twenty-seven HPV genotypes detected using reverse line blot (Roche Diagnostics) (6, 11, 16, 18, 26, 31, 33, 35, 39, 40, 42, 45, 51–59, 66, 68, 73, 82–84). d

q

694/1030 participants were members of a regiment that did not remain in the military area for the full year and so were not included in the 1-year follow-up cohort.

1382

p

o

n

HPV sampled from urethra, glans/corona, shaft, and scrotum.

Anatomical Site of HPV Sampling

Larke et al

Four studies of predominately HIV-negative men evaluated the effect of circumcision on HPV incidence (Table 2), including 1 RCT from Uganda and cohort studies from Mexico and the United States [25, 34, 36, 37]. HPV incidence was relatively high in all 4 populations (21%–35% cumulative risk over 12 months). Male circumcision was protective against incident infection in the RCT (adjusted RR, 0.67; 95% CI, 0.50–0.91) but not in the cohort studies [25, 36, 37]. The meta-analysis was restricted to the 3 studies, which reported a RR (1 study presented an OR) and showed some evidence of a protective effect of circumcision on HPV incidence (summary RR, 0.75; 95% CI, 0.57–0.99; P-heterogeneity 5 0.35; I2 5 5.8%). Circumcision and HPV Clearance

Four studies looked at clearance of HPV infections where the majority of men were HIV negative [25, 34–36] (Table 2). Cleared infections were defined as follows: positive HPV status at baseline and negative HPV status for the same genotype 12 months later [25]; negative status following an initial positive status for the same HPV genotype over any 6month period during 18 months of follow-up [36]; acquisition of an incident infection during follow-up followed by negative status for the same HPV genotype at $2 subsequent visits (every 2 months) [35]; and negative status following an initial positive status for the same HPV genotype between any 2 sequential visits (0, 6, 12, and 24 months) [34]. Three studies [25, 34, 36] showed strong evidence that circumcision was associated with increased clearance (Table 2). No effect of circumcision was observed (RR, 0.96; 95% CI, .71– 1.32) in the remaining study [35]. The summary estimate for those studies that reported a RR (3 out of 4 studies) was 1.33 (95% CI, .89–1.98) with substantial between-study heterogeneity (P-heterogenity 5 0.03; I2 5 73.0%).

Figure 2. Prevalence of human papillomavirus (HPV) by external anogenital anatomical site. HPV sampled at the urethra, glans/corona, shaft, scrotum, perianal area, anus, and foreskin in Nielson (2009) study [18] and the glans/corona, shaft, scrotum, and foreskin in the Hernandez (2008) [17] and Weaver (2004) [12] studies.

Association of Circumcision and HPV Among HIV-Infected Men

One RCT of male circumcision was conducted among HIVinfected men [33] (Tables 1 and 2). This showed an effect of circumcision on prevalence of high-risk HPV genotypes 24 months after the procedure (OR, 0.49; 95% CI, 0.26–0.93). This study also showed a high incidence of HPV (50% over 24 months) and borderline evidence of a protective effect of circumcision on HPV incidence (RR, 0.68; 95% CI, 0.44–1.04) (Table 2). Cleared HPV infections were defined as participants who cleared an infection with the same HPV genotype between enrollment and 24-month follow-up. There was no evidence that circumcision was protective against clearance (RR, 1.09; 95% CI, 0.94–1.27) (Table 2). Description of Genital Warts Studies

Sixteen studies from 15 papers (11 cross-sectional and 5 case control) evaluated the association of circumcision with prevalence of genital warts [22, 38–50] along with 1 cohort study on incidence [44] (Table 3). Circumcision prevalence varied from 21% (United Kingdom) to 87% (Kenya). Circumcision status was evaluated by a clinician in 11 studies, self-reported in 4, and the method not reported in 2. Warts were diagnosed clinically or cytologically in 14 studies and self-reported in 3 studies. Prevalence of warts ranged from

3.6%–62%. Only 2 studies adjusted for covariates, which included age and at least 1 measure of sexual behavior [40, 48]. One study presented insufficient information to calculate a 95% CI [51]. Association of Circumcision and Genital Warts

Overall, there was no evidence of an association of circumcision and prevalent genital warts, although there was marked between-study heterogeneity (summary OR, 0.93; 95% CI, 0.65–1.33; P-heterogeneity , .001; I2 5, 78.1%; Figure 3B). Similarly, no association was seen among the subgroups of studies evaluated (Figure 4). Considerable heterogeneity was observed between the estimates from many subgroups, although there was less heterogeneity among the studies that presented adjusted estimates (I2 5 10.9%) and among populations at high risk for STIs (I2 5 0%). There was no evidence of publication bias (Egger test P 5 0.13). DISCUSSION Circumcised men are at substantially lower risk for prevalent genital HPV infection than non-circumcised men. Few studies assessed the effect of circumcision on HPV incidence or clearance, but data suggest strong evidence of an effect of

Male Circumcision and HPV Infection

d

JID 2011:204 (1 November)

d

1383

Figure 3. Association between male circumcision and prevalence of any penile HPV genotypes (A) or genital warts (B ) among predominantly HIVnegative populations. A subgroup of the randomized populations from the RCTs was tested for HPV infection (56% at 21-month visit in South Africa [19] and 50% at 24-month visit in Uganda [4, 34]).

1384

d

JID 2011:204 (1 November)

d

Larke et al

Figure 4. Subgroup analyses: Male circumcision and HPV prevalence (A) or genital warts. (B) * Method not stated for 3 papers.   The study-specific estimates that contributed to these meta-analyses differed from those presented in Table 1 and are summarized in the supplementary material. à Among men who reported sex with women only. § Method not reported for one paper. { Clinical exam with or without acetic acid and penoscopy or cytology of exfoliated cells

circumcision on reducing HPV incidence and weak evidence of an effect on increasing HPV clearance. There was no evidence of an association of circumcision with genital warts. Heterogeneity

Most studies showed a protective effect of circumcision, but we observed considerable heterogeneity in the magnitude of the effect.

Anatomical Site of HPV Sampling

A key methodological issue for assessing the association of circumcision and genital HPV is the anatomical site of sampling [9, 10, 13]. A strong effect of circumcision amongst studies that measured HPV at all of the glans/corona, penile shaft and scrotum was observed (OR, 0.58; 95% CI, 0.40–0.82; n 5 10), which indicated that there is still a strong overall effect even if there is a change in the distribution of HPV following circumcision.

Male Circumcision and HPV Infection

d

JID 2011:204 (1 November)

d

1385

1386

Table 3.

Summary of Studies on the Association of Male Circumcision and Genital Warts

d

JID 2011:204 (1 November)

Country and year of study

Study population Study size (high/low riska) (% circumcisedb)

First Author

Study design

Aynaud [39]

CrossFrance, Male partners sectional published 1999 of women with genital HPV infection (low risk)

Aynaud [38]

CrossFrance, sectional 1991–1992

Male partners of women with genital condyloma or intraepithelial neoplasia (low risk)

Bleeker [22]

CrossNetherlands, sectional 2002

Dermatology clinic attendees (non-STI clinic) (low risk)

Bleeker [22]

CrossNetherlands, sectional 1995–2002

Partners of women with CIN (low risk)

224 (4%)

Cook [40]

CrossUnited States, sectional 1990

STI clinic attendees (low risk)f,g

Dave [41]

Genital warts

Current/past Assessment of Assessment history of warts genital warts of circumcision

Crude OR (95% CI)

Current warts

Clinical exam Clinical exam with acetic acid and penoscopyc,d

0.53 (0.28–1.02)

1000 (26%) 50%

Current warts

Clinical exam with acetic acid and penoscopyd

Clinical exam

0.76 (0.57–1.02)

112 (17%) 14.2%e

Current warts

Clinical exam with acetic acid and penoscopyd

Clinical exam

0.00 (0.00–1.0)

62.0%e

Current warts

Clinical exam with acetic acid and penoscopyd

Clinical exam

1.20 (0.30–0.51)

1448 (83%) 16.7%

Current warts

Clinical examh

Clinical exam

2.03 (1.45–2.89)

CrossUnited Kingdom, General (low risk) sectional 2000

4777 (21%) 3.6%

History of warts Self- reportedj

Self-reported

1.04 (0.67–1.63)

Dinh [42]

CrossUnited States, sectional 1999–2004

General (low risk)

4110 (72%) 3.9%

History of warts Self-reportedj,k

Self-reported

1.94 (1.28–3.05)

Donovan [43]

CrossAustralia, sectional 1990–1991

STI clinic attendees (high risk)

History of warts Clinical examh

Clinical exam

0.92 (0.47–1.79)

Hart[51]

CrossAustralia, sectional 1988–1991

STI clinic attendees 12 170 (29%) 15% (high risk)

Current warts

Clinical exam

Not stated

1.25 (P , 0.001)m

Lavreys [44]

Cohort

Trucking company employees (high risk)

Current warts

Clinical examh

Clinical exam

RR, 0.77o (0.23–2.5)

Mallon [50]

Case control United Kingdom, Dermatology 1994–1997 clinic attendees (low risk) CrossUnited Kingdom, STI clinic sectional published 1991 attendees (high risk)

d

210 (23%) 54%

Larke et al

Mandal [45]

Kenya, 1993–1997

300 (62%) 16.7%l

746 (87%) 1.4 per 100 pyrn

343 (48%) Not available History of warts Clinical examh,p Clinical exam

0.34 (0.14–0.77)

105 (16%) 27%q

1.60 (0.44–5.50)

Current warts

Cytology of Not reported exfoliated cells

Adjusted OR (95% CI)

1.43 (1.11–2.00)i

Table 3. Continued Country and year of study

Study population Study size (high/low riska) (% circumcisedb)

Genital warts

Current/past Assessment of Assessment history of warts genital warts of circumcision

Crude OR (95% CI)

First Author

Study design

Oriel [52]

Case control United Kingdom, STI Clinic 1967–1970 attendees (high risk)r CrossAustralia, 1981 STI clinic sectional attendees (low risk)f

263 (28%) Not available

Current warts

Clinical exam

Clinical exam

0.66 (0.36–1.22)

568 (61%) 7.3%

Current warts

Clinical examh

Clinical exam

1.54 (0.98–2.42)s

Case controlt

General (low risk)

100 (43%) 6.0%

History of wartsu

Self reportedj

Self-reported

0.96 (0.12–7.53)

General (low risk)v

237 (83%) Not available

Current wartsw Clinical examh

Self-reported

1.02 (0.47–2.27)

Clinical examh

Clinical exam

0.00 (0.00–0.23)

Parker [46]

Tseng [47]

United States, 1975–1985

Van Den Eeden [48] Case control United States, 1987–1991 Wilson [49]

Case control Canada, paper Military presented 1945 (low risk)

1018 (48%) Not available

Current warts

Abbreviations: CI, confidence interval; CIN, cervical intraepithelial neoplasia; OR, odds ratio; pyr, person years; RR, rate ratio; STI, sexually tranmitted infection. a

Study populations categorized as high or low risk for STIs (see Supplementary methods for details).

b

Proportion of entire study population circumcised for cross-section and cohort studies; proportion of controls circumcised for case-control studies.

c

Warts were diagnosed by clinician, but the patients reported if this was the first time they sought medical attention for warts.

d

Both microscopic and macroscopic warts evaluated.

d

q

All cases of warts were microscopic and had no macroscopic evidence of warts; exfoliated cells showed warty atypia with or without dyskaryosis through cytology.

r

Comparison group comprised patients with gonorrhea but no genital warts.

JID 2011:204 (1 November)

Male Circumcision and HPV Infection

e

Penile flat lesions.

f

Controls were groups of patients without any STDs.

g

Participants with STIs other than genital warts were not included in the analysis.

h

Macroscopic warts evaluated.

i

Adjusted for age, ethnicity, number of sex partners in past month, place of residence, and other sexually transmitted infections.

j

Type of warts (micro-or macroscopic) not reported.

k l m

Patients reported that a doctor or healthcare professional had diagnosed them with warts. Proportion with history of genital warts. 95% confidence interval not reported.

n

Incidence of genital warts (First appearance of warts during follow-up), where median follow-up was 21 and 20 months for uncircumcised and circumcised men, respectively.

o

Incidence rate ratio.

p

Warts diagnosed clinically either on history or examination.

s

Age adjusted.

t

Patients identified through cancer registry; neighborhood controls. Data presented from controls only.

u

Warts that occurred more than 2 years before time of interview.

v

Patients included men with incident genital warts; controls were age-matched men randomly selected from the same health cooperative.

w x

Incident cases during 1987–1991, in which patients reported this as the first time medical attention was sought for warts. Adjusted for age, clinic, marital status, reference year, and number of sex partners.

Adjusted OR (95% CI)

0.90 (0.40–2.00)x

d

1387

Circumcised men were consistently less likely than uncircumcised men to have HPV infection at the glans/corona and the urethra. These site-specific effects possibly occur because the foreskin provides a suitable environment around the glans for HPV infection [13], and HPV type-specific concordance has been shown between the glans/corona and foreskin in uncircumcised men that possibly reflects simultaneous infection or autoinoculation [17]. The effect observed at the urethra is unlikely due to detection bias in uncircumcised men. In a nested substudy among men in the control group (who were offered circumcision at the end of the trial), Auvert et al (2009) showed that HPV prevalence did not differ in samples taken from the urethra before and after circumcision. There was less evidence of an effect at sites distal to the foreskin and semen. This may in part be due to lower statistical power because HPV prevalence is lower at the anus and in semen, or it may be a true site-specific differential effect of circumcision. Less heterogeneity was observed in some of the site-specific estimates compared with the estimates of HPV at any anatomical site. Assuming a true difference in the effect of circumcision by anatomical site, this would explain some of the heterogeneity in the effect estimates. The residual heterogeneity seen in subgroup analyses is likely a reflection of the diversity in the methods used and sites sampled to assess HPV infection. The method used to sample HPV may be a source of heterogeneity. However, it would be difficult to interpret any observed variation between sampling method because this could be due to variation in (1) the person using the method, (2) the effectiveness of methods at different anatomical sites, or (3) the sampling method itself. Included Studies

Our review included 14 additional papers not included in the most recent systematic review of circumcision and HPV infection [11], including 4 additional papers on incidence/clearance [33–36] and 1 among HIV-positive men [33]. We included 2 papers that were not identified by previous reviews because these papers did not include circumcision as a keyword [20, 29]. These 2 papers found no evidence of an effect on HPV prevalence, suggesting possible ascertainment bias when circumcision was used as a search term in a systematic search. While the RCTs provided important additional information, the trials were designed to evaluate the association between circumcision and HIV incidence, not HPV. Only subgroups of men were tested for HPV infection (56% at a 21-month visit in South Africa [19], 50% (HIV negative) and 22% (HIV positive) at a 24-month visit in Uganda [33, 34]). HPV was not measured at baseline in the South African trial, so the sequence of circumcision and outcome cannot be determined. In women, the median time to clearance of prevalent infections is 1 year [53], although 2 studies suggest it is ,6 months in men [35, 54]. Nonetheless, the prevalence of HPV in the Ugandan 1388

d

JID 2011:204 (1 November)

d

Larke et al

trial declined considerably between baseline and follow up (62%–36%) among circumcised men but remained relatively high in uncircumcised men (63%–51%). In addition to the published results from the Ugandan and South African RCTs, data from the Kenyan trial also showed a strong effect of male circumcision on HPV incidence (OR, 0.4; 95% CI, 0.3-0.5) and clearance (OR, 0.3; 95% CI; 0.3-0.5) [55]. Genital Warts

No effect of circumcision was observed on genital warts. The method for diagnosis of genital warts and histology of warts varied considerably between studies. However, we found a strong effect of circumcision on low-risk HPV genotypes, which included genotypes 6 and 11 that are commonly associated with warts, with little heterogeneity. There was insufficient information to look at individual HPV genotypes. The lack of effect on genital warts may in part be due to detection bias if genital warts are more commonly reported and/or detected in circumcised men. Implications for Women

HPV prevalence is estimated to be between 15% and 71% among penile cancer cases [8], but it is not thought to be a necessary cause for penile cancer. However, HPV infection in men has consequences for female partners. The association between male circumcision and cervical cancer and HPV infection in female partners has been evaluated in a number of observational studies. Some studies showed a beneficial effect [23, 56], while others showed no evidence of any effect [57–60]. However, secondary analyses from a RCT of male circumcision on HIV infection among men and their female partners showed that female partners of circumcised HIV-negative men were at substantially reduced risk of becoming infected with high-risk HPV genotypes as compared with partners of uncircumcised men (incidence RR, 0.77; 95% CI, .63–.93) [61]. There was some evidence that circumcised men were at lower risk of incident infection with HPV genotypes 16 and 18, which are implicated in cervical neoplasia. Circumcised men were also more likely to clear these infections, although these results were not statistically significant [34]. Women would also benefit indirectly from a lower prevalence of HPV in their male partners. CONCLUSIONS Circumcision services are being expanded as an HIV prevention strategy, especially in countries in sub-Saharan Africa with high HIV prevalence. Sub-Saharan Africa has a high prevalence of HPV infection and a high incidence of invasive cervical cancer [62], but HPV vaccines are not yet available in many settings. Expansion of circumcision services in this region represents an opportunity to reduce HPV infection as well as HIV in men, with resulting benefit to women.

Supplementary Data Supplementary data are available at http://www.jid.oxfordjournals.org/ online. Supplementary materials consist of data provided by the author that are published to benefit the reader. The posted materials are not copyedited. The contents of all supplementary data are the sole responsibility of the authors. Questions or messages regarding errors should be addressed to the author.

Notes Acknowledgments. N. L., H. W., and S. T. designed the search strategy. N. L. and H. W. reviewed abstracts and papers and extracted the data. N. L. performed the meta-analysis and wrote the first draft of the paper. H. W., S. T., and I. S. commented in detail on the drafts and approved the final version. Financial support. This work was supported by the World Health Organization, which provided funding for N.L. and S.T. H.W. was funded by the Medical Research Council (United Kingdom). Potential conflicts of interest. None reported. 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. WHO/UNAIDS. New data on male circumcision and HIV prevention: Policy and programme implications: conclusions and recommendations. Montreux: WHO/UNAIDS, 2007. 2. Gray RH, Serwadda D, Tobian AA, et al. Effects of genital ulcer disease and herpes simplex virus type 2 on the efficacy of male circumcision for HIV prevention: analyses from the Rakai trials. PLoS Med 2009; 6:e1000187. 3. Sobngwi-Tambekou J, Taljaard D, Nieuwoudt M, Lissouba P, Puren A, Auvert B. Male circumcision and Neisseria gonorrhoeae, Chlamydia trachomatis and Trichomonas vaginalis: observations after a randomised controlled trial for HIV prevention. Sex Transm Infect 2009; 85:116–20. 4. Tobian AA, Serwadda D, Quinn TC, et al. Male circumcision for the prevention of HSV-2 and HPV infections and syphilis. N Engl J Med 2009; 360:1298–309. 5. Weiss HA, Thomas SL, Munabi SK, Hayes RJ. Male circumcision and risk of syphilis, chancroid, and genital herpes: a systematic review and meta-analysis. Sex Transm Infect 2006; 82:101–9; discussion 110. 6. Szabo R, Short RV. How does male circumcision protect against HIV infection? BMJ 2000; 320:1592–4. 7. Price L, Johnson K, Rattray R, et al. Circumcision is associated with significant changes in the penis bacterial microbiota [abstract 1062]. In: Program and abstracts of the 16th Conference on Retroviruses and Opportunistic infections (Montreal, Canada). Alexandria, VA: CROI: 2009. 8. Rubin MA, Kleter B, Zhou M, et al. Detection and typing of human papillomavirus DNA in penile carcinoma: evidence for multiple independent pathways of penile carcinogenesis. Am J Pathol 2001; 159:1211–8. 9. Van Howe RS. Human papillomavirus and circumcision: a metaanalysis. J Infect 2007; 54:490–6. 10. Castellsague X, Albero G, Cleries R, Bosch FX. HPV and circumcision: a biased, inaccurate and misleading meta-analysis. J Infect 2007; 55: 91–3; author reply 93–6. 11. Bosch FX, Albero G, Castellsague X. Male circumcision, human papillomavirus and cervical cancer: from evidence to intervention. J Fam Plann Reprod Health Care 2009; 35:5–7. 12. Weaver BA, Feng Q, Holmes KK, et al. Evaluation of genital sites and sampling techniques for detection of human papillomavirus DNA in men. J Infect Dis 2004; 189:677–85. 13. Auvert B, Lissouba P, Sobngwi-Tambekou J. Reply to Van Howe. J Infect Dis 2009; 200:833.

14. Molijn A, Kleter B, Quint W, van Doorn LJ. Molecular diagnosis of human papillomavirus (HPV) infections. J Clin Virol 2005; 32(Suppl 1):S43–51. 15. Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ 1997; 315:629–34. 16. Higgins JP, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in meta-analyses. BMJ 2003; 327:557–60. 17. Hernandez BY, Wilkens LR, Zhu X, et al. Circumcision and human papillomavirus infection in men: a site-specific comparison. J Infect Dis 2008; 197:787–94. 18. 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. J Infect Dis 2009; 199:7–13. 19. Auvert B, Sobngwi-Tambekou J, Cutler E, et al. 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. J Infect Dis 2009; 199:14–9. 20. Aynaud O, Poveda JD, Huynh B, Guillemotonia A, Barrasso R. Frequency of herpes simplex virus, cytomegalovirus and human papillomavirus DNA in semen. Int J STD AIDS 2002; 13:547–50. 21. 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. Sex Transm Dis 2004; 31:601–7. 22. Bleeker MC, Hogewoning CJ, Voorhorst FJ, et al. HPV-associated flat penile lesions in men of a non-STD hospital population: less frequent and smaller in size than in male sexual partners of women with CIN. Int J Cancer 2005; 113:36–41. 23. Castellsague X, Bosch FX, Munoz N, et al. Male circumcision, penile human papillomavirus infection, and cervical cancer in female partners. N Engl J Med 2002; 346:1105–2. 24. Giuliano AR, Lazcano E, Villa LL, et al. Circumcision and sexual behavior: factors independently associated with human papillomavirus detection among men in the HIM study. Int J Cancer 2009; 124:1251–7. 25. Lajous M, Mueller N, Cruz-Valdez A, et al. Determinants of prevalence, acquisition, and persistence of human papillomavirus in healthy Mexican military men. Cancer Epidemiol Biomarkers Prev 2005; 14:1710–6. 26. Muller EE, Chirwa TF, Lewis DA. Human papillomavirus (HPV) infection in heterosexual South African men attending sexual health services: associations between HPV and HIV serostatus. Sex Transm Infect 2010; 86:175–80. 27. Ng’ayo MO, Bukusi E, Rowhani-Rahbar A, et al. Epidemiology of human papillomavirus infection among fishermen along Lake Victoria Shore in the Kisumu District, Kenya. Sex Transm Infect 2008; 84:62–6. 28. Ogilvie GS, Taylor DL, Achen M, Cook D, Krajden M. Self-collection of genital human papillomavirus specimens in heterosexual men. Sex Transm Infect 2009; 85:221–5. 29. Rombaldi RL, Serafini EP, Villa LL, et al. Infection with human papillomaviruses of sexual partners of women having cervical intraepithelial neoplasia. Braz J Med Biol Res 2006; 39:177–87. 30. Shin HR, Franceschi S, Vaccarella S, et al. Prevalence and determinants of genital infection with papillomavirus, in female and male university students in Busan, South Korea. J Infect Dis 2004; 190:468–76. 31. Svare EI, Kjaer SK, Worm AM, Osterlind A, Meijer CJ, van den Brule AJ. Risk factors for genital HPV DNA in men resemble those found in women: a study of male attendees at a Danish STD clinic. Sex Transm Infect 2002; 78:215–8. 32. Vaccarella S, Lazcano-Ponce E, Castro-Garduno JA, et al. Prevalence and determinants of human papillomavirus infection in men attending vasectomy clinics in Mexico. Int J Cancer 2006; 119:1934–9. 33. Serwadda D, Wawer M, Makumbi F, et al. Circumcision of HIVinfected men: effects on high risk human papillomavirus infections in a randomized trial in Rakai, Uganda. J Infect Dis 2010; 201:1463–9. 34. Gray RH, Serwadda D, Kong X, et al. Male circumcision decreases acquisition and increases clearance of high risk human papillomavirus

Male Circumcision and HPV Infection

d

JID 2011:204 (1 November)

d

1389

35.

36.

37.

38.

39.

40. 41.

42.

43.

44.

45.

46. 47.

48.

in HIV-negative men: a randomized trial in Rakai, Uganda. J Infect Dis 2010; 201:1455–62. Hernandez BY, Shvetsov YB, Goodman MT, et al. Reduced clearance of penile human papillomavirus infection in uncircumcised men. J Infect Dis 2010; 201:1340–3. Lu B, Wu Y, Nielson CM, et al. Factors associated with acquisition and clearance of human papillomavirus infection in a cohort of US men: a prospective study. J Infect Dis 2009; 199:362–71. Partridge JM, Hughes JP, Feng Q, et al. Genital human papillomavirus infection in men: incidence and risk factors in a cohort of university students. J Infect Dis 2007; 196:1128–36. Aynaud O, Ionesco M, Barrasso R. Penile intraepithelial neoplasia. Specific clinical features correlate with histologic and virologic findings. Cancer 1994; 74:1762–7. Aynaud O, Piron D, Bijaoui G, Casanova JM. Developmental factors of urethral human papillomavirus lesions: correlation with circumcision. BJU Int 1999; 84:57–60. Cook LS, Koutsky LA, Holmes KK. Circumcision and sexually transmitted diseases. Am J Public Health 1994; 84:197–201. Dave SS, Fenton KA, Mercer CH, Erens B, Wellings K, Johnson AM. Male circumcision in Britain: findings from a national probability sample survey. Sex Transm Infect 2003; 79:499–500. Dinh TH, Sternberg M, Dunne EF, Markowitz LE. Genital warts among 18- to 59-year-olds in the United States, national health and nutrition examination survey, 1999–2004. Sex Transm Dis 2008; 35:357–60. Donovan B, Bassett I, Bodsworth NJ. Male circumcision and common sexually transmissible diseases in a developed nation setting. Genitourin Med 1994; 70:317–20. Lavreys L, Rakwar JP, Thompson ML, et al. Effect of circumcision on incidence of human immunodeficiency virus type 1 and other sexually transmitted diseases: a prospective cohort study of trucking company employees in Kenya. J Infect Dis 1999; 180:330–6. Mandal D, Haye KR, Ray TK, Goorney BP, Stanbridge CM, Corbitt G. Prevalence of occult human papillomavirus infection, determined by cytology and DNA hybridization, in heterosexual men attending a genitourinary medicine clinic. Int J STD AIDS 1991; 2:351–5. Parker SW, Stewart AJ, Wren MN, Gollow MM, Straton JA. Circumcision and sexually transmissible disease. Med J Aust 1983; 2:288–90. Tseng HF, Morgenstern H, Mack T, Peters RK. Risk factors for penile cancer: results of a population-based case-control study in Los Angeles County (United States). Cancer Causes Control 2001; 12:267–77. Van Den Eeden SK, Habel LA, Sherman KJ, McKnight B, Stergachis A, Daling JR. Risk factors for incident and recurrent condylomata

1390

d

JID 2011:204 (1 November)

d

Larke et al

49. 50. 51. 52. 53.

54.

55.

56.

57.

58.

59.

60.

61.

62.

acuminata among men. A population-based study. Sex Transm Dis 1998; 25:278–84. Wilson R. Circumcision and venereal disease. Can Med Assoc J 1947; 56:54–6. Mallon E, Hawkins D, Dinneen M, et al. Circumcision and genital dermatoses. Arch Dermatol 2000; 136:350–4. Hart G. Risk profiles and epidemiologic interrelationships of sexually transmitted diseases. Sex Transm Dis 1993; 20:126–36. Oriel JD. Natural history of genital warts. Br J Vener Dis 1971; 47:1–13. Mueller NE, Birmann B, Parsonnet J, Schiffman M, Stuver S. Chapter 26. Infectious agents. In Schottenfeld D, Fraumeni JF eds: Cancer epidemiology and prevention, 3rd ed. Oxford: Oxford University Press, 2006: xviii, 1392 p. Giuliano AR, Lu B, Nielson CM, et al. Age-specific prevalence, incidence, and duration of human papillomavirus infections in a cohort of 290 US men. J Infect Dis 2008; 198:827–35. Smith JS, Backes DM, Bailey RC, et al. The effect of male circumcision on incident and persistent penile HPV infections in men from Kenya. In: 26th International Papillomavirus Conference & Clinical and Public Health Workshops; Montreal, QC. Jul 3-8 2011. Agarwal SS, Sehgal A, Sardana S, Kumar A, Luthra UK. Role of male behavior in cervical carcinogenesis among women with one lifetime sexual partner. Cancer 1993; 72:1666–9. Asiimwe S, Whalen CC, Tisch DJ, Tumwesigye E, Sethi AK. Prevalence and predictors of high-risk human papillomavirus infection in a population-based sample of women in rural Uganda. Int J STD AIDS 2008; 19:605–10. Hernandez BY, Wilkens LR, Zhu X, et al. Transmission of human papillomavirus in heterosexual couples. Emerg Infect Dis 2008; 14: 888–94. Kjaer SK, de Villiers EM, Dahl C, et al. Case-control study of risk factors for cervical neoplasia in Denmark. I: Role of the ‘‘male factor’’ in women with one lifetime sexual partner. Int J Cancer 1991; 48:39–44. Winer RL, Hughes JP, Feng Q, et al. Condom use and the risk of genital human papillomavirus infection in young women. N Engl J Med 2006; 354:2645–54. Wawer MJ, Tobian AAR, Kigozi G, et al. Circumcision of HIV-infected men and transmission of human papillomavirus to female partners: analyses of data from a randomised trial in Rakai, Uganda. Lancet 2011; 11:604–12. Louie KS, de Sanjose S, Mayaud P. Epidemiology and prevention of human papillomavirus and cervical cancer in sub-Saharan Africa: a comprehensive review. Trop Med Int Health 2009; 14:1287–302.