- Research article
- Open Access
- Open Peer Review
High genital prevalence of cutaneous human papillomavirus DNA on male genital skin: the HPV Infection in Men Study
BMC Infectious Diseasesvolume 14, Article number: 677 (2014)
The genital skin of males hosts a diversity of HPV genotypes and uncharacterized HPV genotypes. Previously we demonstrated that a specific viral genotype was not identified in 14% of all genital specimens (i.e., HPV unclassified specimens) using the Roche Linear Array method. Our goal was to identify and assess the prevalence of individual HPV types among genital HPV unclassified specimens collected in the HIM Study population, at enrollment, and examine associations with socio-demographic and behavioral characteristics.
Genital skin specimens of men that were considered unclassified (HPV PCR positive, no genotype specified) at enrollment were typed by sequencing amplified PGMY09/11 products or cloning of PGMY/GP+ nested amplicons followed by sequencing. PGMY/GP+ negative specimens were further analyzed using FAP primers. HPV type classification was conducted through comparisons with sequences in the GenBank database.
Readable nucleotide sequences were generated for the majority of previously unclassified specimens (66%), including both characterized (77%) and yet uncharacterized (23%) HPV types. Of the characterized HPV types, most (73%) were Beta [β]-HPVs, primarily from β-1 and β-2 species, followed by Alpha [α]-HPVs (20%). Smokers (current and former) were significantly more likely to have an α-HPV infection, compared with any other genus; no other factors were associated with specific HPV genera or specific β-HPV species.
Male genital skin harbor a large number of β-HPV types. Knowledge concerning the prevalence of the diverse HPV types in the men genital is important to better understand the transmission of these viruses.
Human papillomaviruses (HPVs) are a diverse group of circular double-strand DNA viruses of approximately 8.000pb. Persistent infection with high-risk HPV is not only the major etiological factor for cervical cancer development but also for a high proportion of tumors of the vagina, vulva, penis, anal canal and oropharynx .
A novel papillomavirus is defined whenever the complete L1 sequence of a viral genome differs by at least 10% from that of all characterized types . To date, more than 180 different HPV types have been fully sequenced and characterized. The majority of HPVs comprise three genera: Alpha [α]-, Beta [β]- and Gamma [γ]-papillomavirus. α-HPV types have been predominantly isolated from mucosal and genital lesions, and include the 13 viral types classified as high-risk oncogenic (carcinogen type I) by the World Health Organization due to their increased prevalence in cervical cancer. β- and γ-HPV types were mainly isolated from the skin and for this reason have been grouped together as cutaneous HPV. β-HPVs are additionally divided into 6 different species including 45 HPV types, of which β-1 and β-2 are the most diverse and prevalent. Similarly, γ-HPVs are separated into 20 species, which together include 50 different HPV types. Partial DNA sequence information originating from several studies points towards the existence of hundreds of putative novel HPV types of the β- and γ-HPV genera ,. Recent data describe high prevalence of cutaneous HPV types at body sites that are different from those in which they were originally isolated -.
During the last decade, there has been growing interest in understanding HPV natural history and related diseases among men. The HPV Infection in Men (HIM) Study is an ongoing, prospective genital HPV natural history study of over 4.000 men aged 18–70 years residing in Brazil (São Paulo), Mexico (Cuernavaca), and the United States (USA; Tampa). In this study, over 66% of men tested positive for HPV at their first study visit; however, the viral type could not be identified in 14% of the genital specimens using the Linear Array method (i.e., HPV PCR positive, no genotype specified) . Our goal was to examine the prevalence of HPV types among genital HPV unclassified specimens collected among HIM Study participants at enrollment, and examine associations with socio-demographic and behavioral characteristics.
Men were enrolled between 2005 and 2009 and reported no prior diagnosis of genital warts or anogenital cancers, and had no recent symptoms of or treatment for a sexually transmitted infection, including HIV/AIDS. Men completed a pre-enrollment (baseline) visit, were enrolled on completion of their second (enrollment) visit two weeks post-baseline, and subsequently followed approximately every six months for up to four years. Details of the HIM Study cohort are described elsewhere ,. The present cross-sectional prevalence study was conducted among the first 3,105 men who completed their enrollment visits between July 2005 and December 2008. The IRD of the University of South Florida, Tampa, USA, the Instituto Mexicano del Seguro Social, and the Instituto de Salud Publica de Mexico, Mexico, and the Centro de Referência e Treinamento DST//Aids, Brazil, approved all study procedures, and all participants provided written informed consent.
Genital specimen collection, DNA extraction, and HPV detection
At each study visit, participants completed a risk factor questionnaire and underwent a clinical examination of the external genital skin. Exfoliated cells were obtained from the coronal sulcus/glans penis, penile shaft, and scrotum using Dacron swabs (Digene, Gaithersburg, MD, USA), and these samples were further combined and stored at −80°C. HPV DNA was extracted using the QIAamp Media MDx Kit (Qiagen, Valencia, CA, USA). Samples were PCR amplified using PGMY09/11 generic primers and genotyped using the Roche Linear Array (Roche Molecular Diagnostics, Alameda, CA, USA). This technique allows testing for the presence of 37 α-HPV types commonly detected at the cervix (high-risk: 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 68; low-risk: 6, 11, 26, 40, 42, 53, 54, 55, 61, 62, 64, 66, 67, 69, 70, 71, 72, 73, 81, 82, 82 subtype IS39, 83, 84, 89 [CP6108]) . Samples that tested PCR-positive and Linear Array-negative were considered unclassified and underwent additional typing.
Typing of unclassified samples
Initially, purified DNAs were HPV genotyped by direct sequencing of PGMY09/11 primers PCR amplimers or cloning of these amplicons followed by sequencing. Next, 1 μl of PGMY09/11 negative product was used in a nested PCR using GP5+/6+ primers  and positive samples were cloned and sequenced. Finally, nested PGMY09/11-GP5+/6+ PCR negative samples were submitted to a novel amplification reaction employing FAP59/64 primers , and amplimers were analyzed solely by direct sequencing. All PCRs were performed using AmpliTaq Gold polymerase (Perkin-Elmer, Foster City, CA, USA). Before direct sequencing, PCR products were purified using the EXO SAP-IT (GE Healthcare, Buckinghamshire, UK). Sequencing reactions were performed in an ABI 3130XL Genetic Analyzer (AB Applied Biosystems, CA, USA) using the BigDye Terminator v3.1 Cycle Sequencing kit (AB Applied Biosystems, CA, USA). Sequence identity was determined by comparison with the BlastN database; sequences with identity scores higher than 90% within at least 200 bp were conclusively typed.
The prevalence and genotype-distribution of unclassified genital HPV infections were estimated among all men enrolled in the HIM Study before December 2008. Infections were considered HPV-negative, HPV-positive for a characterized type, HPV-positive for an uncharacterized type, or untyped (i.e., inconclusive). Characterized HPV types were further classified according to genus and species. Socio-demographic and behavioral risk factors thought to be associated with presence of 1) a specific genus (α-, β-, γ-, or other HPV types), and 2) a specific β-HPV species were evaluated using exact Pearson chi-square tests and Monte Carlo methods. All statistical tests were two-sided and attained statistical significance at α = 0.05. Analyses were performed using SAS version 9.3 (SAS Institute, Cary, NC, USA).
Of the 3105 men enrolled in the HIM Study through December 2008, 1126 (36.3%) were HPV negative at their enrollment visit, 1572 (50.6%) had at least one of the 37 α-HPV types detected in the Linear Array assay, and 407 (13.1%) had an unclassified HPV type detected. Of the 407 samples with unclassified HPV at enrollment, 404 were available for additional characterization. These include 165, 134, and 105 men from the USA, Brazil and Mexico, respectively. Table 1 presents socio-demographic characteristics of participants included in this analysis, by country. Men with unclassified HPV were more likely to be younger (18–30 years, 52.5%), white (54.7%), non-Hispanic (63.0%), single (52.2%), uncircumcised (56.7%), never smokers (66.2%), moderate drinkers (41.6 and men who have sex with women (89.2%).
HPV type distribution of previously unclassified HPV infections is shown in Table 2. HPV could not be genotyped among 121 (30.0%) FAP59/64 PCR-positive specimens; direct sequencing was inconclusive due to overlapping peak patterns. Among 265 (65.6 specimens with successful sequences, 60 (22.6%) were considered uncharacterized; sequence comparison analysis showed high identity to partial HPV nucleotide sequences in GenBank without type designation at this time. Among the other 205 (77.4%) specimens with already characterized HPV types, 41 (20.0%) belonged to the α-HPV genus, 150 (73.2%) to the β-HPV genus, and 14 (6.8%) to the γ-HPV genus. The distribution of HPV types for each genus was similar across countries, though a higher proportion of γ-HPV types were detected in Brazil and Mexico as compared to the USA.
The genotype-specific distribution of characterized HPV types is found in Figure 1. A total of 27 α-HPV types, 32 β-HPV types, and 14 γ-HPV types were detected. The most common α-HPV types included HPV-2 and HPV-6 (2.0 4 each), the most common β-HPVs included HPV-107 (11.7%), HPV-38 (7.8%), HPV-120 (6.8%), and HPV-17 (4.8%), and the most common γ-HPVs included HPV-4, HPV-134 and HPV-147 (1.0%, 2 each).
Of the risk factors, only cigarette smoking was significantly associated with detection of α-, β-, or γ-HPV types (Table 3); smokers (current and former) were significantly more likely to have an α-HPV infection, compared with any other genus (p = 0.039). No other factors were significantly associated with the detection of a specific HPV genus. Nevertheless, the prevalence of α-HPV types was non-significantly highest among participants living in the USA, and those aged 18–30 years, white, non-Hispanic, single, circumcised, and reported ≥8 lifetime sex partners. In contrast, the prevalence of β-HPV types was equally as high among men who were single and married/cohabitating, and highest among men who were uncircumcised and reported <8 lifetime sex partners. γ-HPV types were more common in Brazil and Mexico than USA, and among men who were married/cohabitating, uncircumcised, and reported <8 lifetime sex partners. Further, none of the risk factors examined were significantly associated with particular β-HPV species detection.
Recently we reported our initial observations concerning the broad distribution of β- and γ-HPV types in genital specimens of men participating in the HIM Study at enrollment through three years of follow-up . We then analyzed the β-HPV type distribution in external genital lesion (EGL) specimens and preceding normal genital skin specimens among HIM participants who developed external genital lesions . In the present study, we expand our previous work by presenting viral sequence characterization of all unclassified HPV genital specimens at the time of participant enrollment and by evaluating the association of HPV detection with participant characteristics and sexual behavior. Some studies have reported the presence of cutaneous HPV types in penile carcinomas  and in cervical and penile condylomas ,. However, to our knowledge the HIM study is unique in searching for cutaneous HPV in normal male genital skin specimens.
Among the 404 samples analyzed, 4.5% were truly HPV negative, i.e. PCR negative using any of the three primers sets, suggesting that a small proportion of unclassified HPV detection presumably represents spurious amplification. In 30.0% of the specimens, unreadable sequences were generated after direct sequencing of FAP54/69 PCR products; among these the presence of more than one HPV type is suspected . This hypothesis has been confirmed through sequencing of a number of clones from each sample whenever overlapping peak patterns were observed and also by the Luminex technology ,. DNA from multiple cutaneous HPV types has been commonly detected not only in the normal skin and cutaneous tumors ,, but also in oral cavity specimens , and in the anogenital tract . Together these data suggest that cutaneous HPV types may be a commensal component of the microbiological flora of the human skin ,-. We also observed that nucleotide sequences of 60 specimens mostly consisted of partial FAP59/64 sequences previously deposited in the GenBank database ,, and were therefore considered uncharacterized HPV types. Currently, more than 200 different putative HPV types are pending full genome characterization, with the majority clustering within the γ-HPV genus . We detected several α-HPVs that should have been identified with the Linear-array genotyping test. These false-negative results are likely due to low viral copy numbers which were detected on a different PCR reaction followed by sequencing, a more sensitive technical approach.
Overall, unclassified HPV detection in HIM participants was significantly associated with younger age, white race, and non-Hispanic ethnicity. β-HPV types were detected among individuals residing in three countries and with diverse habits; however, slight differences in participant characteristics were demonstrated in men with α-HPVs versus β-HPVs, though none of these differences achieved statistical significance. In contrast to α-HPV prevalence in the genital region of men from the HIM study , β-HPV DNA detection was not significantly associated with sexual behaviors, indicating that viral types from this genus may not be transmitted sexually. The risk factor questionnaire completed by HIM participants focuses mostly on penetrative sexual behavior and lacks information on masturbation practices and non-penetrative sexual activities, which may be helpful in understanding the mechanism of transmission for this group of viruses. It is plausible that direct hand-to-penis contact may favor viral transmission of cutaneous HPV to the genital area. In fact, evidence of hand transmission to the anogenital region has been reported for α-HPV types . Furthermore, α-HPVs have been detected on fingers of women who are infected at the cervix . Nevertheless, the possibility remains that the detection of β-HPVs in male genitals may also represent deposition of virions shed from other body sites with productive infections.
HPV characterization of skin specimens from healthy individuals and squamous cell carcinomas of the skin (SCC) indicated that HPV types of β-2 species predominates in tumor samples . In the present study, the analysis of participant characteristics and sexual behavior stratified by β-HPV species found no significant associations. As with all studies, there are limitations that could influence the interpretation of the results obtained. This study was conceived to identify HPV types that were grouped together as unclassified HPV types at enrollment in the HIM Study. For this reason, this analysis excluded specimens for which an HPV type was already assigned by Linear Array, although these specimens may have also harbored β- and γ-HPV types. Co-detection with α- and β-HPV types was shown to be common in the oral cavity  and in penile condylomas .
The role of HPV in the development of non-melanoma skin cancer is not well defined. Furthermore, it is unknown whether HPV DNA detected in healthy skin represents viral particles or superficial cells containing episomal HPV genomes . It has been suggested that skin infected with β-HPV types harbor increased susceptibility to ultraviolet radiation induced DNA damage, ultimately resulting in skin cancer . At present, it is consensus that HPV-5 and HPV-8 are associated with benign and malignant lesions of the cutaneous disease epidermodysplasia verruciformis . Recent data show the presence of β-HPV types in the oral cavity ,, cervical and penile condylomas , and also in cervicovaginal cells ,. Nevertheless, at least for the oral cavity, only viral transcripts of α-HPV were detected in non-malignant lesions and cancer co-harboring β- or γ-HPV DNA . Most studies examining the role of HPV in the development of male EGL have evaluated mucosal HPV types of the α-HPV genus. We observed that β-HPV DNA is prevalent in EGLs (61.1%) and that most viral DNAs are detected on the normal genital skin before EGL development . At present it is unclear how β-HPV types may influence lesion development in the genital region of men. Although detection of cutaneous HPV types in the genital region may simply reflect virions released from other parts of the body once β- and γ-HPV are ubiquitous throughout the body, a possible role for these viruses as co-factors in the oncogenesis of penile cancer remains to be clarified. Studies conducted prospectively in a larger number of patients will contribute to this knowledge.
Detection of HPV types not previously assigned by Roche´s Linear-Array in the genitals of participants of a natural history study of HPV in men was significantly associated with younger age, white race, and non-Hispanic ethnicity. Among these HPVs, we detected a broad range of α- and β-HPVs. Further research is needed to determine how these cutaneous HPV are transmitted.
Ethical approval was given by the University of South Florida (IRB# 102660), Ludwig Institute for Cancer Research, Centro de Referência e Treinamento em Doenças Sexualmente Transmissíveis e AIDS, and Instituto Nacional de Salud Pública de México.
- HIM Study:
HPV Infection in Men Study
Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, El Ghissassi F, Benbrahim-Tallaa L, Guha N, Freeman C, Galichet L, Cogliano V: WHO International Agency for Research on Cancer Monograph Working Group. A review of human carcinogens–Part B: biological agents. Lancet Oncol. 2009, 10: 321-322. 10.1016/S1470-2045(09)70096-8.
Bernard HU, Burk RD, Chen Z, van Doorslaer K, zur Hausen H, de Villiers EM: Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments. Virology. 2010, 401: 70-79. 10.1016/j.virol.2010.02.002.
Chouhy D, Bolatti EM, Pérez GR, Giri AA: Analysis of the genetic diversity and phylogenetic relationships of putative human papillomavirus types. J Gen Virol. 2013, 94: 2480-2488. 10.1099/vir.0.055137-0.
Forslund O: Genetic diversity of cutaneous human papillomaviruses. J Gen Virol. 2007, 88: 2662-2669. 10.1099/vir.0.82911-0.
Bottalico D, Chen Z, Dunne A, Ostoloza J, McKinney S, Sun C, Schlecht NF, Fatahzadeh M, Herrero R, Schiffman M, Burk RD: The oral cavity contains abundant known and novel human papillomaviruses from the Betapapillomavirus and Gammapapillomavirus genera. J Infect Dis. 2011, 204: 787-792. 10.1093/infdis/jir383.
Ekström J, Forslund O, Dillner J: Three novel papillomaviruses (HPV109, HPV112 and HPV114) and their presence in cutaneous and mucosal samples. Virology. 2010, 397: 331-336. 10.1016/j.virol.2009.11.027.
Pierce Campbell CM, Messina JL, Stoler MH, Jukic DM, Tommasino M, Gheit T, Rollison DE, Sichero L, Sirak BA, Ingles DJ, Abrahamsen M, Lu B, Villa LL, Lazcano-Ponce E, Giuliano AR: Cutaneous human papillomavirus types detected on the surface of male external genital lesions: a case series within the HPV Infection in Men Study. J Clin Virol. 2013, 58: 652-659. 10.1016/j.jcv.2013.10.011.
Sichero L, Pierce Campbell CM, Ferreira S, Sobrinho JS, Luiza Baggio M, Galan L, Silva RC, Lazcano-Ponce E, Giuliano AR, Villa LL: Broad HPV distribution in the genital region of men from the HPV infection in men (HIM) study. Virology. 2013, 443: 214-217. 10.1016/j.virol.2013.04.024.
Akogbe GO, Ajidahun A, Sirak B, Anic GM, Papenfuss MR, Fulp WJ, Lin HY, Abrahamsen M, Villa LL, Lazcano-Ponce E, Quiterio M, Smith D, Schabath MB, Salmeron J, Giuliano AR: Race and prevalence of human papillomavirus infection among men residing in Brazil, Mexico and the United States. Int J Cancer. 2012, 131: E282-E291. 10.1002/ijc.27397.
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. Cancer Epidemiol Biomarkers Prev. 2008, 17: 2036-2043. 10.1158/1055-9965.EPI-08-0151.
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. Lancet. 2011, 377: 932-940. 10.1016/S0140-6736(10)62342-2.
Gravitt PE, Peyton CL, Apple RJ, Wheeler CM: Genotyping of 27 human papillomavirus types by using L1 consensus PCR products by a single-hybridization, reverse line blot detection method. J Clin Microbiol. 1998, 36: 3020-3027.
de Roda Husman AM, Walboomers JM, van den Brule AJ, Meijer CJ, Snijders PJ: The use of general primers GP5 and GP6 elongated at their 3' ends with adjacent highly conserved sequences improves human papillomavirus detection by PCR. J Gen Virol. 1995, 76: 1057-1062. 10.1099/0022-1317-76-4-1057.
Forslund O, Antonsson A, Nordin P, Stenquist B, Hansson BG: A broad range of human papillomavirus types detected with a general PCR method suitable for analysis of cutaneous tumours and normal skin. J Gen Virol. 1999, 80: 2437-2443.
del Pino M, Bleeker MC, Quint WG, Snijders PJ, Meijer CJ, Steenbergen RD: Comprehensive analysis of human papillomavirus prevalence and the potential role of low-risk types in verrucous carcinoma. Mod Pathol. 2012, 25: 1354-1363. 10.1038/modpathol.2012.91.
Ekström J, Bzhalava D, Svenback D, Forslund O, Dillner J: High throughput sequencing reveals diversity of Human Papillomaviruses in cutaneous lesions. Int J Cancer. 2011, 129: 2643-2650. 10.1002/ijc.26204.
Sturegard E, Johansson H, Ekström J, Hansson BG, Johnsson A, Gustafsson E, Dillner J, Forslund O: Human papillomavirus typing in reporting of condyloma. Sex Transm Dis. 2013, 40: 123-129.
Berkhout RJ, Tieben LM, Smits HL, Bavinck JN, Vermeer BJ, ter Schegget J: Nested PCR approach for detection and typing of epidermodysplasia verruciformis-associated human papillomavirus types in cutaneous cancers from renal transplant recipients. J Clin Microbiol. 1995, 33: 690-695.
Forslund O, Lindelöf B, Hradil E, Nordin P, Stenquist B, Kirnbauer R, Slupetzky K, Dillner J: High prevalence of cutaneous human papillomavirus DNA on the top of skin tumors but not in "Stripped" biopsies from the same tumors. J Invest Dermatol. 2004, 123: 388-394. 10.1111/j.0022-202X.2004.23205.x.
Antonsson A, Forslund O, Ekberg H, Sterner G, Hansson BG: The ubiquity and impressive genomic diversity of human skin papillomaviruses suggest a commensalic nature of these viruses. J Virol. 2000, 74: 11636-11641. 10.1128/JVI.74.24.11636-11641.2000.
Astori G, Lavergne D, Benton C, Höckmayr B, Egawa K, Garbe C, de Villiers EM: Human papillomaviruses are commonly found in normal skin of immunocompetent hosts. J Invest Dermatol. 1998, 110: 752-755. 10.1046/j.1523-1747.1998.00191.x.
Boxman IL, Berkhout RJ, Mulder LH, Wolkers MC, Bouwes Bavinck JN, Vermeer BJ, ter Schegget J: Detection of human papillomavirus DNA in plucked hairs from renal transplant recipients and healthy volunteers. J Invest Dermatol. 1997, 108: 712-715. 10.1111/1523-1747.ep12292090.
de Villiers EM, Lavergne D, McLaren K, Benton EC: Prevailing papillomavirus types in non-melanoma carcinomas of the skin in renal allograft recipients. Int J Cancer. 1997, 73: 356-361. 10.1002/(SICI)1097-0215(19971104)73:3<356::AID-IJC9>3.0.CO;2-Z.
Forslund O, Iftner T, Andersson K, Lindelof B, Hradil E, Nordin P, Stenquist B, Kirnbauer R, Dillner J, de Villiers EM: Cutaneous human papillomaviruses found in sun-exposed skin: Beta-papillomavirus species 2 predominates in squamous cell carcinoma. J Infect Dis. 2007, 196: 876-883. 10.1086/521031.
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. J Infect Dis. 2009, 199: 362-371. 10.1086/596050.
Sonnex C, Strauss S, Gray JJ: Detection of human papillomavirus DNA on the fingers of patients with genital warts. Sex Transm Infect. 1999, 75: 317-319. 10.1136/sti.75.5.317.
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. J Infect Dis. 2007, 196: 1128-1136. 10.1086/521192.
Heideman DA, Waterboer T, Pawlita M, Delis-van Diemen P, Nindl I, Leijte JA, Bonfrer JM, Horenblas S, Meijer CJ, Snijders PJ: Human papillomavirus-16 is the predominant type etiologically involved in penile squamous cell carcinoma. J Clin Oncol. 2007, 25: 4550-4556. 10.1200/JCO.2007.12.3182.
Wang J, Aldabagh B, Yu J, Arron ST: Role of human papillomavirus in cutaneous squamous cell carcinoma: A meta-analysis. J Am Acad Dermatol. 2014, 70: 621-629. 10.1016/j.jaad.2014.01.857.
Viarisio D, Decker KM, Aengeneyndt B, Flechtenmacher C, Gissmann L, Tommasino M: Human papillomavirus type 38 E6 and E7 act as tumour promoters during chemically induced skin carcinogenesis. J Gen Virol. 2013, 94: 749-752. 10.1099/vir.0.048991-0.
Orth G: Genetics and susceptibility to human papillomaviruses: epidermodysplasia verruciformis, a disease model. Bull Acad Natl Med. 2010, 194: 923-940.
Lindel K, Helmke B, Simon C, Weber KJ, Debus J, de Villiers EM: Cutaneous human papillomavirus in head and neck squamous cell carcinomas. Cancer Invest. 2009, 27: 781-787. 10.1080/07357900802653456.
Johansson H, Bzhalava D, Ekström J, Hultin E, Dillner J, Forslund O: Metagenomic sequencing of "HPV-negative" condylomas detects novel putative HPV types. Virology. 2013, 440: 1-7. 10.1016/j.virol.2013.01.023.
Chen Z, Schiffman M, Herrero R, Desalle R, Burk RD: Human papillomavirus (HPV) types 101 and 103 isolated from cervicovaginal cells lack an E6 open reading frame (ORF) and are related to gamma-papillomaviruses. Virology. 2007, 360: 447-453. 10.1016/j.virol.2006.10.022.
Paolini F, Rizzo C, Sperduti I, Pichi B, Mafera B, Rahimi SS, Vigili MG, Venuti A: Both mucosal and cutaneous papillomaviruses are in the oral cavity but only alpha genus seems to be associated with cancer. J Clin Virol. 2013, 56: 72-76. 10.1016/j.jcv.2012.09.016.
The authors would like to thank the HIM Study Teams in the USA (Moffitt Cancer Center, Tampa, FL: HY Lin, C Gage, K Eyring, K Kennedy, K Isaacs, A Bobanic, DJ Ingles, MT O’Keefe, MR Papenfuss, MB Schabath, AG Nyitray), 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: 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: J Salmerón, M Quiterio Trenado, 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).
Financial Support: Ludwig Institute for Cancer Research; Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) [Grants numbers 10/15282-3 to LS and 08/57889-1 to LLV]; Conselho Nacional de Desenvolvimento Científico e Tencnológico (CNPq) [Grant number 573799/2008-3 to LLV]. The infrastructure of the HIM Study cohort was supported through a grant from the National Cancer Institute (NCI), National Institutes of Health (NIH), to ARG (CA R01CA098803). CMPC was supported by a Postdoctoral Fellowship (PF-13-222-01 - MPC) from the American Cancer Society.
ARG receives research funding from Merck Sharp & Dohme. LLV and ARG are consultants of Merck Sharp & Dohme for HPV vaccines. None of the other authors have conflicts of interest to report.
ARG, LLV, and EL-P were the principal and co-principal investigators, conceived the study, assured funding, identified clinical investigators and study personnel, and made critical revision of the manuscript. LS and CMPC contributed to the conception and design of this study were responsible for analysis and interpretation of data, drafting and revision of the manuscript. WP did the statistical analysis and made comments to the manuscript. SF and JSS wrote the protocols, conducted all molecular analysis and participated in the interpretation of the data and revision of the manuscript. MLB, LG and RCS coordinated the field work, study implementation and data collection, and made critical revision of the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.