Prevalence of human papillomavirus infection and phylogenetic analysis of HPV-16 E6 variants among infected women from Northern Brazil
Infectious Agents and Cancer volume 9, Article number: 25 (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.
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.
Socio-demographic data was obtained for a total of 143 women. Data analysis revealed that these women had 34 years median age (SD ± 13.74 y; 95% CI 32-36), were married (48%), had basic education as the main school level (58%), had housewife occupation (35%), were not exposed to alcohol (58%) or to tobacco (86%), had no history of STD’s (Sexual Transmitted Diseases) (84%), had 16 years old (SD ± 2.9 y; 95% CI 15.5-16.5) as the median age of the first sexual intercourse and had a median of 10 sexual partners (SD ± 32.3 partners; 95% CI 4.3-15.3). The higher median observed on the variable number of sexual partners was probably influenced by data related to women that referred to be sex workers (n = 7). When these data was removed from analysis, observed media of sexual partners was 5.5 partners (SD ± 8.6 partners; 95% CI 4-7). There was no observed statistical association between the described variables and the occurrence of HPV infection in our population (p > 0.05) (Table 1).
All of the 143 samples included in this study had normal cytology results. 23 HPV negative samples (23/101; 23%) had exhibited no amplification of internal control gene and were considered inadequate for subsequent analysis. Linear Array HPV Genotyping Test® (Roche Molecular Systems, Alameda, CA) results revealed that 78 samples were negative (78/120; 65%) and 42 were positive (42/120; 35%) for HPV infection. Among the HPV positive samples, the high risk HPV-16 was the most prevalent type (14%), followed by HPV-52 (9.5%) and HPV-45 (7%). Among low risk HPV types HPV-62 (7%) was the most frequent. The others high and low risk HPV types were identified in low frequencies (<2,5% each). From those HPV-16 positive samples, six were HPV-16 single infection (14.3%) and five were multiple infections that included HPV-16 (12%) (Table 2).
There was a variation of HPV prevalence between different age groups (Figure 1) and age was associated with HPV prevalence (Z(U) = 14.62; p = <0.0001). HPV of low and high risk was more prevalent among women at younger ages (<25 years) (37.8%; 95% CI 31.6%-44%) and hrHPV types were more prevalent among 25-34 age-group women (31.3%; 95% CI 26.9%-36.4%) (Table 3).Phylogenetic analysis was performed for a total of 8 samples among the HPV-16 positive samples (73%). The phylogenetic tree segregated into four major branches that could be recognized as European (EUR), Asian (As), African (including Afr1 and Afr2) and American (including Asian-American and North-American) lineages. The HPV-16 variants were classified as European (87.5%) and American (12.5%) (Figure 2).
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.
Stanley M: Immunobiology of HPV and vaccines. Gynecol Oncol. 2009, 109: S15-S21.
De Villiers E-M: Cross-roads in the classification of papillomaviruses. Virology. 2013, 445: 2-10.
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.
Castellsagué X: Natural history and epidemiology of HPV infection and cervical cancer. Gynecol Oncol. 2008, 110: S4-S7.
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.
Frazer I: Correlating immunity with protection for HPV infection. Int J Infect Dis. 2007, 11 (suppl. 1): S10-S16.
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.
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.
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.
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.
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-
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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-
IARC Working Group on the Evaluation of Carcinogenic Risks to Human: A review of human carcinogens. Part B: Biological agents. 2009, Lyon, France
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-
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.
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.
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.
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.
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.
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.
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.
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. 1994
BioEstat 5.3, Instituto de Desenvolvimento Sustentável Mamirauá. 2011, [cited 2011 Sep 16]. Available from: http://www.mamiraua.org.br/download/
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).
The authors declare that they have no competing interests.
BPTL performed the genotyping of HPV variants, data and statistical analysis and lead drafting of the manuscript. ESCJ participated in the HPV variant analysis and performed phylogenetic analysis. WAM is member of the HPV at Instituto Evandro Chagas and is one of the researchers of the Juruti Project. FP participated in the study design and data interpretation. CGMS participated in the study design and data interpretation. RVDS was responsible for study design, for the initial HPV-DNA investigation, is a researcher and member of the HPV group at Instituto Evandro Chagas. All authors read and approved the final manuscript.
Bruna Pedroso Tamegão-Lopes and Rodrigo Vellasco Duarte Silvestre contributed equally to this work.
About this article
Cite this article
Tamegão-Lopes, B.P., Sousa-Júnior, E.C., Passetti, F. et al. Prevalence of human papillomavirus infection and phylogenetic analysis of HPV-16 E6 variants among infected women from Northern Brazil. Infect Agents Cancer 9, 25 (2014). https://doi.org/10.1186/1750-9378-9-25