Prevalence of human papillomavirus infection and phylogenetic analysis of HPV-16 E6 variants among infected women from Northern Brazil
- Bruna Pedroso Tamegão-Lopes†1,
- Edivaldo Costa Sousa-Júnior1,
- Fabio Passetti2,
- Carlos Gil Ferreira2,
- Wyller Alencar de Mello1 and
- Rodrigo Vellasco Duarte Silvestre†1Email author
© Tamegão-Lopes et al.; licensee BioMed Central Ltd. 2014
Received: 14 March 2014
Accepted: 8 July 2014
Published: 5 August 2014
The main cause of cervical cancer in the world is high risks human papillomavirus infection (mainly represented by HPV-16 and HPV-18), that are associated to the development of malign transformation of the epithelium. HPV prevalence exhibits a wide geographical variability and HPV-16 variants have been related to an increased risk of developing cervical intraepithelial lesion. The aim of this study was to describe DNA-HPV prevalence and HPV-16 variants among a women population from Northern Brazil.
One hundred and forty three women, during routine cervical cancer screening, at Juruti Project, fulfilled an epidemiological inquiry and were screened through a molecular HPV test. HPV-16 variants were determined by sequencing the HPV-16 E6 open reading frame.
Forty two samples were considered HPV positive (29.4%). None of those had abnormal cytology results. HPV prevalence varied between different age groups (Z(U) = 14.62; p = <0.0001) and high-risk HPVs were more frequent among younger ages. The most prevalent type was HPV-16 (14%) and it variants were classified, predominantly, as European (87.5%).
HPV prevalence in our population was higher than described by others and the most prevalent HPV types were high-risk HPVs. The European HPV-16 variant was the most prevalent among HPV-16 positive samples. Our study reinforces the fact that women with normal cytology and a positive molecular test for high-risk HPVs should be submitted to continuous follow up, in order to verify persistence of infection, promoting an early diagnosis of cervical cancer and/or its precursors.
KeywordsPrevalence HPV HPV-16 Variants
Human papillomavirus (HPV) is a large group of epitheliotropic viruses of more than 160 different types. HPVs are very successful infectious agents – inducing chronic infections that have no systemic sequelae, rarely killing the host and shedding large amounts of infectious particles for transmission. However, a dozen of the 40 HPV types that infect the human body are associated with probable or definite oncogenic risk.
High-risk HPVs (hrHPV), which includes HPV type 16 (HPV-16), are etiological agents of cervical cancer[3–5], being the persistent of hrHPV infection, known as the integrated HPV DNA into host genome, in the presence of other environmental and host factors, a necessary step for the development of high-grade cervical intraepithelial neoplasia and invasive cervical cancer[6–9]. Cervical cancer is the third most common cancer in women and more than 85% of the global burden occurred in developing countries, causing 275,000 deaths in 2008, with 31,400 of those deaths occurring in Latin America.
HPV molecular variants are defined taxonomically based on L1 coding region of the DNA sequence. The comparative nucleotide sequence analysis, that includes ORF E6, have shown that HPV-16 variants, which differ in nucleotide sequence by no more than 2% in the coding regions and until a 5% in the non coding region, evolved into six major phylogenetic branches: European (E), Asian (As), Asian-American (AA), African-1 (Af1), African-2 (Af2) and North-American (NA1)[5, 12, 13]. Instead only a small portion of HPV-16 infections persist, several studies from Europe and the Americas, based principally on the sequencing of E6 and/or the LCR, had suggested that HPV-16 variants can influence viral persistence and the development of cervical cancer[14, 15].
The aim of this study was to describe HPV infection prevalence and related epidemiological variables, as well as, the HPV-16 variants circulating in a population of women from Northern Brazil.
Distribution of socio-demographic characteristics of the 143 women included in this study
N = 143 n (%)
Single or live with partner
Exposition to tobacco?
History of STD?
Condyloma: 3 (27.5)
Syphilis: 3 (27,5)
Vaginal discharge: 1 (9)
Gonorrhea: 1 (9)
Trichomoniasis: 1 (9)
Not informed: 2 (18)
Use of condom in sexual intercourse?
Distribution of HPV types according to Array HPV Genotyping Test® results
HPV types (N = 42)
16 + 18
16 + 31
16 + 59
16 + 31 + 52
16 + 51 + 56
52 + 58
56 + 39
61 + 81
35 + 52 + 59
45 + 51 + 62
45 + 55 + 61
45 + 58 + 61
11 + 53 + 84 + GTCP6108
Distribution of age-group according to HPV infection
Age at sample collection
Nº of women
HPV + (%) (CI 95)
hrHPV + (%) (CI 95)
17 (37.8) (31.6%-44%)
11 (24.4) (20.4%-28.4%)
12 (37.5) (31.4%-43.6%)
10 (31.3) (26.9%-36.4%)
7 (18.4) (15.4%-21.4%)
5 (13.2) (11%-15.4%)
3 (17.6) (14.7%-20.5%)
2 (11.8) (9.9%-13.7%)
3 (27.3) (22.8%-31.8%)
1 (9.1) (7.6%-10.6%)
In the present study we found a 29.4% HPV prevalence between women without abnormal cytology screened by Pap Smear Test, at Juruti, Northern Brazil. This prevalence was higher than that described by SanJosé et al, in meta-analysis of a systematic literature review of HPV prevalence among women with normal cervical cytology worldwide, where global and South America prevalence’s were 10.4% and 12.3%, respectively. As well as, that described by Bruni et al, in meta-analysis of one million women with normal cytological findings, where global and South America prevalence’s were 11.7% and 15.3%, respectively. A broader range of HPV prevalence (from 10.4% to 24.5%) has been described in Brazilian studies including women with normal cytology. Our result was similar to that described by the Latin America Screening Study (24.5%) and to data related to Amazonas State (29.2%). However, it was lower than described in others Brazilian studies including women populations from Rio Grande do Sul (34%), São Paulo (58.4%), Distrito Federal (62%) and Rio de Janeiro (63.6%)[19–21].
A marked variation of HPV prevalence among age groups was observed in our population. Women with less than 25 years and comprehend between 25-34 years were more predisposal to have a positive HPV-DNA test, including for hrHPV types. Some studies in Central and South America have demonstrated a first peak prevalence of HPV infection in women below age 25 and a second peak after age 55[22, 23] which corroborates our data. The first peak suggest that women acquire infection in the first few years after become sexually active and the second peak suggest that middle-age women harbor long-term HPV persistent infections or new incident infections, which might also reflect changes in hormonal or immunological status[24, 25].
HPV-16 was the most prevalent type observed in our study (14.3%), followed by HPV-52 and HPV-45. The most prevalent HPV types described in women with normal cytology are HPV-16, -58, -31, -18, -45 and -52, however, it is well documented that this prevalence is lower than that observed among women with HSIL (high-grade lesions) and SCC (squamous-cell carcinoma). HPV-16 prevalence in our population was higher than that described in Natal (Northeast Brazil) (12.7%) and lower than that observed in Rio Grande do Sul (19.6%), Rio de Janeiro (19%), São Paulo (18.2%) and Amazonas (31.8%)[18–21].
Phylogenetic analysis based on E6 sequences of eight HPV-16 isolates allowed the identification of seven European and only one non-European variants. A previous study realized in Belém (Northern Brazil) exhibited a higher prevalence of Asian-American (43.2%) and European (42%) variants among women with cervical cancer. Previous data suggests that women harboring non-European HPV-16 variants have an 2- to 9- fold increased risk of cervical cancer[29, 30], however, it has been reported that HPV-16 European variants persists longer in white women[31–33].
Genital HPV infection is one of the most common sexually transmitted infections worldwide. It has been estimated that approximately 10% of women worldwide with normal cytological findings carry a detectable cervical HPV infection, depending on the HPV testing technology, study size, and the age groups and geographical region studied. It is necessary to emphasize that women with one positive DNA test for HPV-16, as well as, for others hrHPV, need to be followed-up to confirm persistent infection, to monitor and to avoid late diagnosis of possible cervical cancer and its precursors, since it is well documented that HPV-16 and HPV-18 account for 70% of all cervical cancer cases worldwide. This will be only achieved through consistent cervical screening associated with alternative screening tools as molecular techniques.
A major limitation of the present study is represented by the modest sample size and by the absence of a follow-up period to confirm persistence of HPV infection or clearance through timeline, as well as, the lack of ethnicity data, which was expected to be associated with HPV-16 variants. Another limitation was related with no amplification of internal control (β-globin gene) between 23 cervical samples on the Linear Array HPV Genotyping Test®. Those results were probably caused by inadequate sample cellularity (yield DNA specimens), as the DNA extraction technique was held as a constant, with DNA isolated using the same commercial kit.
HPV prevalence (29.4%) in a population of women with normal cytology from Juruti (Northern Brazil) is higher than described by others with the most prevalent HPV types being hrHPVs (HPV-16, -52 and -45). The European HPV-16 variant (87.5%) is the most prevalent among HPV-16 positive samples. HPV-16 is defined as a human carcinogenic agent and is associated with cervical cancer among carriers. Despite being descriptive, the present study reinforces the necessity to follow up women with normal cytology, meanwhile, with a positive molecular test for hrHPVs, to determine persistence of infection and to promote early diagnosis of cervical cancer and its precursors.
Samples were collected from patients included on the Juruti Project, a research project designed to investigate the impact of migratory flow on the incidence and prevalence of sexual transmissible diseases. Juruti is a county located in the West Region of the State of Pará that had an estimated population of 34.338 individuals in the year of 2007. The main economical activities of Juruti are agriculture and mining and this county is marked by migratory flow of people from the States of Amazonas (Northern Brazil), Maranhão and Ceará (Northeastern Brazil).
During fifteen days of March (2007) women participating on the Juruti Project were offered systematic cervical cancer screening (Pap smear test) combining conventional cytology through Cytobrush Plus Cell Collector® (Coopersurgical, Trumbull, Connecticut, US) for optimal collection of exfoliated cells and a commercialized version of the line blot assay – Linear Array HPV Genotyping Test® (Roche Molecular Systems, Alameda, CA) – when they had their routine cervical smear performed in a public hospital. The Juruti Project has been approved by the Instituto Evandro Chagas, statement nº30/05, CAAE 0013.0.072.000-06DC-2008-374. All women were informed and gave their written consent to participate in the study.
HPV cytological and molecular analysis
After samples collection cytobrush was immersed into 15 ml falcon tubes with 3ml of phosphate buffer saline (PBS: 137 mM NaCl, 10 mM Phosphate buffer, 2,7 mM KCl, pH 7.4) and kept at 4ºC according to the manufacturer’s instructions. At the HPV laboratory the falcon tubes were submitted to manual agitation, followed by cytobrush discard and centrifugation of the liquid per fifteen minutes (1200 rpm). Supernatant was discarded and cell pellet was resuspended in 3 mL of PBS solution.
The conventional cytology results were classified according to the recommendations of Brazilian Ministry of Health and Brazilian Society of Cytology which is based on Bethesda’s definition and were performed at Juruti by the local routine professionals. Cervical specimen DNA extraction was performed according to the Linear Array HPV Genotyping Test® kit protocol (Roche Molecular Systems, Alameda, CA), which can detect until thirty seven distinct HPV types in a unique reaction, which includes high and low risk and multiple infections. Only the strips that had an internal positive control test (β-globin gene) were considered proper to genotyping analysis.
The PCR and sequencing of the HPV 16 E6 gene (nt 46-256) was performed as described by Chopjitt et al. The amplification of E6 fragment was confirmed as described beforey. The obtained nucleotide sequences of the E6 region were edited and assembled using UGENE software. Sequences were globally aligned and compared with HPV variants sequences available at GenBank database (GenBank accession number: K02718, NC_001526, AF125673, AF536179, AF486310, AF472508, AF536180, AF472509, AF486324, AF404705, AF534061, AF486299, AF402678, AF486325) using MAFFT version 7 software. Phylogenetic tree was constructed and edited using RaxML version 7.0.4 (1.000 bootstraps) and Dendroscope version 3 software’s, respectively.
Data were analyzed with the software BioEstat 5.3. Fisher test was used to compare categorical variables between two groups (HPV-infected and non-infected women). Wilcoxon-Mann-Whitney test and G test were used to verify a possible association between age and HPV prevalence. After plotting HPV prevalence and hrHPV prevalence among HPV-infected women, in each age-group, a polynomial fit was assessed as the best model for logistic regression.
At the Juruti Project part of the equipments and consumption were sponsored by Alcoa Alumínio S.A., together with Brazilian Ministry of Health. BPTL, FP and CGF acknowledge the support of Fundação do Câncer. FP acknowledges the support of Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).
- Stanley M: Immunobiology of HPV and vaccines. Gynecol Oncol. 2009, 109: S15-S21.View ArticleGoogle Scholar
- De Villiers E-M: Cross-roads in the classification of papillomaviruses. Virology. 2013, 445: 2-10.PubMedView ArticleGoogle Scholar
- Gross G, Pfister H: Role of human papillomaviruses in penile cancer, penile intraepithelial squamous cell neoplasias and in genital warts. Med Microbiol Immunol. 2004, 193: 35-44.PubMedView ArticleGoogle Scholar
- Castellsagué X: Natural history and epidemiology of HPV infection and cervical cancer. Gynecol Oncol. 2008, 110: S4-S7.PubMedView ArticleGoogle Scholar
- Chopjitt P, Ekalaksananan T, Pietong C, Kongyingyoes B, Kleebkaow P, Charoensri N: Prevalence of human papillomavirus type 16 and its variants in abnormal squamous cervical cells in Northeast Thailand. Int J Infect Dis. 2009, 13: 212-219.PubMedView ArticleGoogle Scholar
- Frazer I: Correlating immunity with protection for HPV infection. Int J Infect Dis. 2007, 11 (suppl. 1): S10-S16.PubMedView ArticleGoogle Scholar
- Tornesello ML, Duraturo ML, Giorgi-Rossi P, Sansone M, Piccoli R, Buonaguro L, Buonaguro FM: Human papillomavirus (HPV) genotypes and HPV16 variants in human immunodeficiency virus-positive Italian women. J Gen Virol. 2008, 89: 1380-1389.PubMedView ArticleGoogle Scholar
- Fernandes JV, Meissner RV, Carvalho MGF, Fernandes TAAM, Azevedo PRM, Villa LL: Prevalence of HPV infection by cervical cytologic status in Brazil. Int J Gynaecol Obstet. 2009, 105: 21-24.PubMedView ArticleGoogle Scholar
- Zuna RE, Tuller E, Wentzensen N, Mathews C, Allen RA, Shanesmith R, Dunn ST, Gold MA, Wang SS, Walker J, Schiffman M: HPV16 variant lineage, clinical stage, and survival in women with invasive cervical cancer. Infect Agent Cancer. 2011, 6: 1-9.View ArticleGoogle Scholar
- Ferlay J, Shin H-R, Bray F, Forman D, Mathers C, Parkin DM: Estimates of worldwide burden of cancer in, GLOBOCAN 2008. Int J Cancer. 2008, 2010 (127): 2893-2917.Google Scholar
- Chen Z, Schiffman M, Herrero R, DeSalle R, Anastos K, Segongy M, Sahasrabuddhe VV, Gravitt PE, Hsing AW, Burk RD: Evolution and Taxonomic Classification of Human Papillomavirus 16 (HPV-16)-Related Variant Genomes: HPV31, HPV33, HPV35, HPV52, HPV58 and HPV67. PLoS One. 2011, 6 (5): e20183-PubMedPubMed CentralView ArticleGoogle Scholar
- Alencar TR, Cerqueira DM, Cruz MR, Wyant PS, Ramalho ED, Martins CRF: New HPV-16 European and non-European variants in Central Brazil. Virus Genes. 2007, 35: 1-4.PubMedView ArticleGoogle Scholar
- Huertas-Salgado A, Martín-Gámez DC, Moreno P, Murillo R, Bravo MM, Villa L, Molano M: E6 molecular variants of human paillomavirus (HPV) type 16: An updated and unified criterion for clustering and nomenclatura. Virol. 2011, 410: 201-215.View ArticleGoogle Scholar
- Trottier H, Mahmud S, Prado JCM, Sobrinho JS, Costa MC, Rohan TE, Villa L, Franco EL: Type-Specific Duration of Human Papillomavirus Infection: Implications for Human Papillomavirus Screening and Vaccination. J Infect Dis. 2008, 197: 1436-1447.PubMedView ArticleGoogle Scholar
- Cornet I, Gheit T, Franceschi S, Vignat J, Burk RD, Sylla BS, Tommasino M, Clifford GM, IARC HPV Variant Study Group: Human Papillomavirus Type 16 Genetic Variants: Phylogeny and Classification Based on E6 and LCR. J Virol. 2012, 86 (12): 6855-6861.PubMedPubMed CentralView ArticleGoogle Scholar
- de SanJosé S, Diaz M, Castellsagué X, Clifford G, Bruni L, Muñoz N, Bosch FX: Worldwide prevalence and genotype distribution of cervical human papillomavirus DNA in women with normal cytology: a meta-analysis. Lancet Infect Dis. 2007, 7: 453-459.PubMedView ArticleGoogle Scholar
- Bruni L, Diaz M, Castellsagué X, Ferrer E, Bosch FX, de Sanjosé S: Cervical Human Papillomavirus Prevalence in 5 Continents: Meta-Analysis of 1 Million Women with Normal Cytological Findings. J Infect Dis. 2010, 202 (12): 1789-1799.PubMedView ArticleGoogle Scholar
- Rocha DAP, Filho RAAB, Queiroz AF, Santos CMB: High Prevalence and Genotypic Diversity of the Human Papillomavirus in Amazonian Women, Brazil. Infect Dis Obstet Gynecol. 2013, ID514859: 1-5.Google Scholar
- Coser J, Rocha Boeira T, Simon D, Kazantzi Fonseca AS, Ikuta N, Lunge VR: Prevalence and genotypic diversity of cervical human papillomavirus infection among women from an urban Center in Brazil. Genet Mol Res. 2013, 12 (4): 4276-4285.PubMedView ArticleGoogle Scholar
- Lippman SA, Sucupira MCA, Jones HE, Luppi CG, Palefsky J, van de Wijgert JHHM, Oliveira RLS, Dias RS: Prevalence, distribution and correlates of endocervical human papillomavirus types in Brazilian women. Int J STD AIDS. 2010, 21 (2): 105-109.PubMedPubMed CentralView ArticleGoogle Scholar
- Oliveira-Silva M, Lordello CX, Zardo LMG, Bonvicino CR, Moreira MAM: Human Papillomavirus in Brazilian women with and without cervical lesions. Virol J. 2011, 8: 1-4.View ArticleGoogle Scholar
- Smith JS, Melendy A, Rana RK, Pimenta JM: Age-Specific Prevalence of Infection with Human Papillomavirus in Females: A Global Review. J Adolesc Health. 2008, 43: S5-S25.PubMedView ArticleGoogle Scholar
- Franceschi S, Herrero R, Clifford GM, Snijders PJF, Arslan A, Anh PTH, Bosch FX, Ferreccio C, Hieu NT, Lazcano-Ponce E, Matos E, Molano M, Qiao Y-L, Rajkumar R, Ronco G, de Sanjosé S, Shin H-R, Sukvirach S, Thomas JO, Meijer CJLM, Muñoz N, the IARC HPV Prevalence Surveys Study Group: Variations in the age-specific curves of human papillomavirus prevalence in women worldwide. Int J Cancer. 2006, 119: 2677-2684.PubMedView ArticleGoogle Scholar
- Brismar-Wendel S, Froberg M, Hjerpe A, Anderson S, Johansson B: Age-specific prevalence of HPV genotypes in cervical cytology samples with equivocal or low-grade lesions. Br J Cancer. 2009, 101: 511-517.PubMedPubMed CentralView ArticleGoogle Scholar
- Baussano I, Franceschi S, Gillio-Tos A, Carozzi F, Confortini M, Palma PD, de Lillo M, del Mistro A, de Marco L, Naldoni C, Pierotti P, Schincaglia P, Segnan N, Zorzi M, Giorgi-Rossi P, Ronco G: Difference in overall and age-specific prevalence of high-risk human papillomavirus infection in Italy: evidence from NTCC Trial. BMC Infect Dis. 2013, 13: 238-PubMedPubMed CentralView ArticleGoogle Scholar
- IARC Working Group on the Evaluation of Carcinogenic Risks to Human: A review of human carcinogens. Part B: Biological agents. 2009, Lyon, FranceGoogle Scholar
- Fernandes JV, Meissner RV, Carvalho MGF, Fernandes TAAM, Azevedo PRM, Sobrinho JS, Prado JCM, Villa LL: Prevalence of human papillomavirus in archival samples obtained from patients with cervical pre-malignant and malignant lesions from Northeast Brazil. BMC Res Notes. 2010, 3: 96-PubMedPubMed CentralView ArticleGoogle Scholar
- Junes-Gill K, Sichero L, Maciag PC, Mello W, Noronha V, Villa LL: Human Papillomavirus Type 16 Variants in Cervical Cancer from an Admixtured Population in Brazil. J Med Virol. 2008, 80: 1639-1645.PubMedView ArticleGoogle Scholar
- Villa L, Sichero L, Rahal P, Caballero O, Ferenczy A, Rohan T, Franco EL: Molecular variants of human papillomavirus types 16 and 18 preferentially associated with cervical neoplasia. J Gen Virol. 2000, 81: 2959-2968.PubMedView ArticleGoogle Scholar
- Zehbe I, Lichtig H, Westerback A, Lambert PF, Tommasino M, Sherman L: Rare human papillomavirus 16 E6 variants reveal significant oncogenic potential. Mol Cancer. 2011, 10: 1-13.View ArticleGoogle Scholar
- Muñoz N, Franco EL, Herrero R, Andrus JK, de Quadros C, Goldie SJ, Bosch FX: Recommendations for Cervical Cancer Prevention in Latin America and the Caribbean. Vaccine. 2008, 26S: L96-L107.View ArticleGoogle Scholar
- Bernard H-U, Calleja-Macias IE, Dunn ST: Genome variation of human papillomavirus types: Phylogenetic and medical implications. Int J Cancer. 2006, 118: 1071-1076.PubMedView ArticleGoogle Scholar
- Sichero L, Ferreira S, Trottier H, Duarte-Franco E, Ferenczy A, Franco EL, Villa LL: High grade cervical lesions are caused preferentially by non-European variants of HPVs 16 and 18. Int J Cancer. 2007, 120: 1763-1768.PubMedView ArticleGoogle Scholar
- Clifford G, Franceschi S, Diaz M, Muñoz N, Villa LL: Chapter 3: HPV type-distribution in women with and without cervical neoplastic diseases. Vaccine. 2006, 24S3: S3/26-S3/34.Google Scholar
- Kurman RJ, Solomon D: The Bethesda System for Reporting Cervical/Vaginal Cytologic Diagnoses Definitions, Criteria, and Explanatory Notes for Terminology and Specimen Adequacy. ISBN: 978-0-387-94077-9 (Print) 978-1-4684-0201-8 (Online). Springer- Verlag, New York Inc. 1994Google Scholar
- BioEstat 5.3, Instituto de Desenvolvimento Sustentável Mamirauá. 2011, [cited 2011 Sep 16]. Available from: http://www.mamiraua.org.br/download/
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.