SciELO - Scientific Electronic Library Online

vol.112 issue7Detection and molecular characterisation of a diagnosis escape variant associated with occult hepatitis B virus in BrazilParticipation of the arcRACME protein in self-activation of the arc operon located in the arginine catabolism mobile element in pandemic clone USA300 author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand




Related links


Memórias do Instituto Oswaldo Cruz

Print version ISSN 0074-0276On-line version ISSN 1678-8060

Mem. Inst. Oswaldo Cruz vol.112 no.7 Rio de Janeiro July 2017 


Human papillomavirus type distribution and HPV16 intratype diversity in southern Brazil in women with and without cervical lesions

Gisele R de Oliveira1  2  + 

Valdimara C Vieira2 

Emiliana C Ávila1 

Fabiana Finger-Jardim1 

Thaís DM Caldeira1 

Fabiane AA Gatti1 

Carla V Gonçalves3 

Sandro G Oliveira4 

Vanusa P da Hora1 

Marcelo A Soares2 

Ana MB de Martinez1 

1Universidade Federal do Rio Grande, Escola de Medicina, Laboratório de Biologia Molecular, Rio Grande, RS, Brasil

2Instituto Nacional de Câncer, Programa de Oncovirologia, Rio de Janeiro, RJ, Brasil

3Universidade Federal do Rio Grande, Escola de Medicina, Centro de Obstetrícia e Ginecologia, Rio Grande, RS, Brasil

4Policlínica de Assistência Médica, Rio Grande, RS, Brasil



Increasing evidence suggests that human papillomavirus (HPV) intratype variants (specific lineages and sublineages) are associated with pathogenesis and progression from HPV infection to persistence and the development of cervical cancer.


This study aimed to verify the prevalence of HPV infection and distribution of HPV types and HPV16 variants in southern Brazil in women with normal cytology or intraepithelial lesions.


HPV typing was determined by L1 gene sequencing. To identify HPV16 variants, the LCR and E6 regions were sequenced, and characteristic single nucleotide variants were identified.


A total of 445 samples were studied, with 355 from cervical scrapes and 90 from cervical biopsies. HPV was detected in 24% and 91% of these samples, respectively. The most prevalent HPV types observed were 16 (cervical, 24%; biopsies, 57%) and 58 (cervical, 12%; biopsies, 12%). Seventy-five percent of the HPV16-positive samples were classified into lineages, with 88% defined as lineage A, 10% as lineage D, and 2% as lineage B.


This study identified a high frequency of European and North American HPV16 lineages, consistent with the genetic background of the human population in southern Brazil.

Key words: HPV; lineage; HPV16; intratype variants; Brazil

Human papillomavirus (HPV) infection is one of the main risk factors for the development of cervical cancer (zur Hausen 1976). Cervical cancer is the fourth most common type of malignancy worldwide, with approximately 90% of related deaths occurring in developing countries (Globocan 2012). According to the Brazilian National Cancer Institute (INCA 2016), this is the third most common type of cancer among women and fourth most common cause of cancer-related deaths in Brazil.

Over 200 HPV types have been identified to date, and among these, approximately 60 infect the human genital tract and can be classified according their oncogenic potential as low- or high-risk HPV types. Low-risk types are more frequently associated with genital warts and low-grade squamous intraepithelial lesions (LSIL), whereas high-risk types are more frequently associated with high-grade squamous intraepithelial lesions (HSIL) and cervical cancer (Muñoz et al. 2003, Bernard et al. 2010). HPV6 and 11 are the most prevalent low-risk types in genital warts, and HPV16 and 18 are the most frequent high-risk types found in malignant lesions and are associated with 70% of all cervical cancer cases worldwide. Nonetheless, HPV types 31, 33, 35, 45, 52, and 58 are the six most common types in the world after HPV16 and 18, and are responsible for 30% of cervical cancer cases (de Sanjosé et al. 2010).

Phylogenetic classification of HPV types is based on differences in the L1 gene sequence, which encodes the main capsid protein, with different HPV types showing at least 10% divergence from one another (de Villiers et al. 2004). HPV types can be further classified into lineages, with genome sequences that diverge by 1-10%, or sub-lineages, with genomes that are 0.5-1.0% divergent (Burk et al. 2013).

Ho et al. (1993) classified HPV16 variants from samples from around the world using long control region (LCR) sequences. Four major branches, associated with geographic distribution, were identified: European-Asian (EAS) composed of European (EUR) and Asian (As) sublineages, African 1 (AFR1), African 2 (AFR2), and Asian-American (AA) lineages. A few years later, an additional lineage (North-American) was described (Yamada et al. 1997). A 2013 study by Burk and colleagues updated HPV16 lineages and sublineages using HPV16 complete genomic sequences and an alphanumeric nomenclature: lineage A for lineages previously identified as European and Asian (EAS); lineage B for African 1 (AF1), lineage C for African 2 (AF2), and lineage D for North-American and Asian-American lineages.

Interestingly, the genetic differences between HPV16 sublineages may be associated with their distinct oncogenic potentials. HPV variants may differ biologically and etiologically, leading to differences in the development and behaviour of tumours (Bernard 2002). For example, non-European variants have a general tendency to cause persistent infections and are associated with cervical lesions (Villa et al. 2000). In addition, infections with European variants that carry the T350G polymorphism, which leads to a valine replacing a leucine in the E6 protein, was suggested as an additional risk factor for persistent infection and cervical lesions (Zhang 2015). However, few reports on the distribution of HPV16 intratype variants circulating in Brazil, and especially in southern Brazil, exist. The aim of this study was to analyse the prevalence of HPV infection and the distribution of HPV types and HPV16 genetic variants in HIV-negative and HIV-positive women with or without cervical intraepithelial lesions and carcinoma and who attended centres and hospitals in southern Brazil.


Study population - A descriptive cross-sectional study was carried out in women attending the Gynecological Service of the University Hospital of Federal University of Rio Grande (HU-FURG), a regional reference centre for the treatment of HIV patients and where the Papanicolaou test was performed. Cervical scrapes were collected using a cytobrush (Vagispec®). Women diagnosed with HSIL in the public health system of Rio Grande city also participated in this study, with cervical biopsy samples collected at the Medical Assistance Polyclinic. Biopsy samples were collected directly from lesions. Cervical scrapes and biopsy samples were collected between March 2013 and March 2016.

Informed consent for participation in the study was obtained from all participants in strict compliance with the ethical guidelines for studies involving human subjects in Brazil as required by Resolution #466/2012 of the National Health Council. The study was approved by the Institutional Review Boards (IRB) of Federal University of Rio Grande and Prefeitura Municipal do Rio Grande (protocols # 057/2013 and 92A/2013, respectively), Brazil.

Molecular analysis - After collection, cervical scrapes were stored in cryogenic tubes containing 1 mL of TE buffer (10 mM Tris-HCl pH 8.5; 1 mM EDTA). Cervical biopsies were stored in tubes containing 1 mL of RNAlater®. All samples were stored at -20ºC until DNA extraction and polymerase chain reaction (PCR) amplification of HPV genomic fragments.

Cervical scrapes were centrifuged at 14,000 rpm for 15 min, supernatants were discarded, and cell pellets were used for DNA extraction. Extractions were performed using the Purelink Genomic DNA Kit (Life Technologies, Carlsbad, CA) according to manufacturer’s protocol.

For detection of HPV DNA, nested PCR was performed. The first-round PCR was performed with the external primers MY09/MY11, described by Manos et al. (1989), which amplify a 450-bp fragment of the viral L1 region. The second round was performed with the internal primers GP5/6 (Van den Brule et al. 1990), which amplify a 150-bp internal fragment of L1 (Table I). The reaction mixtures contained 5 µL of DNA or 3 µL of first-round product, respectively, 1× PCR buffer, 2 mM MgCl2, 0.5 mM dNTPs, 1 U of recombinant Taq DNA polymerase (Life Technologies, São Paulo, Brazil), 0.5 mM primers MY09/MY11 or GP5/6, and water to a final volume of 50 µL. The amplified PCR products were subjected to agarose gel electrophoresis, visualised under an ultraviolet light transilluminator, and images were captured. A positive control from a sample confirmed to be HPV-positive by sequencing, a negative control (negative by molecular and cytological assays), and a blank reaction (with no DNA) were used.

TABLE I Primers and cycling conditions used for polymerase chain reaction (PCR) amplification of the L1 region to detect human papillomavirus and the E6 and LCR regions of HPV16 

Primer Sequence Size (bp) Cycling conditions
MY09/MY11 5’CGTCCMAARGGAWACTGATC3’ 5’GCMCAGGGWCATAAYAATGG3’ 450 35 cycles of denaturation at 94ºC/40 s, annealing at 55ºC/30 s, and extension at 72ºC/1 min, and a final step at 72ºC/10 min.
GP5/6 5’TTTGTTACTCTGGTAGATAC3’ 5’GAAAAATAAACTGTAAATCA3’ 150 40 cycles, with denaturation at 95ºC/40 s, annealing at 40ºC/40 s, and extension at 72ºC/40 s.
HPV16 E6 F HPV16 E6R 5’ CGAAACCGGTTAGTATAA 5’ GTATCTCCATGCATGATT 524 35 cycles with denaturation at 94ºC/30 s, annealing at 49ºC/30 s, and extension at 72ºC/45 s. Final extension at 72ºC/5 min.
HPV16 LCR out F HPV16 LCR out R 5’TGCTTGCCAACCATTCCATTGTTTTTT 5’TTGCTTGCAGTACACACATTCTAATATT 531 35 cycles with denaturation at 56ºC/30 s, annealing at 49ºC/30 s, and extension at 72ºC/40 s. Final extension at 72ºC/10 min.
HPV16 LCR in F HPV16 LCR in R 5’GTACATTGTGTCATATAAAATAAATCA 5’TGTGGGTCCTGAAACATTGCAGTTCTCTTTT 382 30 cycles with denaturation at 94ºC/30 s, annealing at 55ºC/30 s, and extension at 72ºC/45 s. Final extension at 72ºC/10 min.

PCR products were purified using an Illustra GFX PCR DNA and Gel Band Purification Kit (GE Healthcare Life Sciences, Piscataway, USA) and sequenced using an ABI Prism BigDye Terminator Cycle Sequencing Ready Reaction Kit (Applied Biosystems, Foster City, USA) with an automated ABI 3130XL sequencer (Thermo Scientific, Waltham, USA). The presence of multiple overlapping electropherogram peaks in the sequencing data were suggestive of multiple HPV infections and were not further typed, whereas clear sequences were assembled and typed. For alignment and editing of sequences and contig assemblies, SeqMan and EditSeq software of the DNA Star package (Laser Gene Inc., Madison, USA) and BioEdit Sequence Alignment Editor were used. The generated sequences were subjected to a search using the BLASTn algorithm (Available from: to determine HPV types based on the closest sequences in the public database.

After genotyping, HPV16-positive samples were used for the assessment of genetic variability. Nested PCR was conducted to amplify the E6 and LCR regions of HPV16 with specific primers and cycling conditions shown in Table I. Reactions were performed in a final volume of 50 µL. PCR products were subjected to electrophoresis in agarose gels, as described above. Fragments were purified and sequenced using the same primers with the kit and platform described above. Sequences were assembled as described and were aligned to HPV16 reference sequences as described by Burk et al. (2013) using MEGA v.4.1. The reference sequences used can be identified by accession numbers K02718 (lineage A), AF536179 (lineage A), HQ644236 (lineage A), AF534061 (lineage A), AF536180 (lineage B), HQ644298 (lineage B), AF472509 (lineage C), HQ644257 (lineage D), AY686579 (lineage D), and AF402678 (lineage D).

To identify HPV16 variants, only samples yielding high quality E6 and LCR sequences were used to identify specific single nucleotide polymorphisms (SNPs) as established by Cornet et al. (2012), which allows the identification of the lineages and sublineages proposed by Burk et al. (2013).

All DNA sequences generated in the present study were submitted to GenBank under accession numbers KY306724-KY306804.


We analysed 445 samples, 355 from cervical scrapes and 90 from cervical biopsies. The average age of participants was 31 years old (SD, ± 12.8 years). Women who were subjected to cervical scrapes were 29.4 years of age, on average (SD ± 10.2 years), whereas those who had cervical biopsy samples collected were an average age of 37.9 years (SD ± 9.5 years). Of the women included in this study, 61.3% were self-described as white, 64.7% as living with a steady partner, 24.4% as having finished high school, and 37.9% as having an income between one and two minimum wages (one Brazilian minimum wage corresponds to US$ 250.00).

Among women who attended the Gynecological Service of HU-FURG, 23.6% (84/355) were HIV-1-positive, whereas all of the women who attended a public clinic were HIV-1-negative. Of the cervical scrape samples analysed, 24% (86/355) were positive for HPV DNA. The HPV prevalence among HIV-1-positive and negative women was 28.5% (24/84) and 23.2% (62/271), respectively.

Among the women who were positive for HPV infection and negative for HIV-1, LSIL were present in six samples (9.6%) and HSIL were found in one sample (1.6%). From the HIV-1-positive women, one sample (4.1%) showed HSIL.

All 86 smear samples found positive for HPV DNA were sequenced for HPV genotyping. In the HIV-1-negative samples with LSIL (n = 6), the genotypes detected were HPV6 (1/6), HPV31 (1/6), HPV33 (3/6), and HPV82 (1/6); in the sample with HSIL (n = 1), HPV16 was identified. From the HIV-1-positive women, HPV16 DNA was detected in the sample with HSIL. The prevalence of HPV genotypes in samples with negative cytological results (n = 78) from HIV-1-positive and -negative women is shown in Table II.

TABLE II Human papillomavirus (HPV) genotyping of samples negative for HPV cytology from HIV-1-positive and -negative women 

HPV Types HIV-1-negative N = 55 HIV-1-positive N = 23 Total N = 78
HPV6 7 (13%) 2 (9%) 9 (14%)
HPV16 12 (21%) 7 (31%) 19 (24%)
HPV18 3 (5%) - 3 (4%)
HPV31 2 (4%) 2 (9%) 4 (5%)
HPV33 3 (5%) 2 (9%) 5 (7%)
HPV35 1 (2%) - 1 (1%)
HPV44 2 (4%) 1 (4%) 3 (4%)
HPV45 7 (13%) - 7 (9%)
HPV53 2 (4%) - 2 (2%)
HPV58 8 (14%) 3 (14%) 11 (14%)
HPV61 1 (2%) 1 (4%) 2 (2%)
HPV67 1 (2%) 1 (4%) 2 (2%)
HPV68 - 1 (4%) 1 (1%)
HPV70 1 (2%) 1 (4%) 2 (2%)
HPV82 4 (7%) 1 (4%) 5 (7%)
HPV83 - 1 (4%) 1 (1%)
HPV85 1 (2%) - 1 (1%)

Of the 90 cervical biopsy samples analysed histopathologically, 33% (30/90) had HSIL, 33% (30/90) had LSIL, and 33% (30/90) were negative for lesions. HPV DNA was detected in 91% (82/90) of the biopsies, with all samples showing LSIL and HSIL determined HPV-positive; among women negative for histopathological lesions, the HPV prevalence was 73% (22/30), with 96% of these women infected with high-risk oncogenic HPV types.

Seventy-seven (94%) of the 82 HPV-positive biopsy samples were genotyped. In the majority of these samples (72/77), only one genotype was found; however, in five samples, multiple HPV infections were indicated by overlapping electropherogram peaks observed in the DNA sequence data. As shown in Table III, HPV16 was the most prevalent genotype detected in women both with and without squamous intraepithelial lesions (SIL) (57%); HPV58 was the second most frequent (12%).

TABLE III Human papillomavirus (HPV) genotyping of biopsy samples (n = 77) 

Genotype No lesions (N = 17) LSIL (N = 30) HSIL (N = 30) Total (N = 77)
HPV11 1 (6%) - - 1 (1%)
HPV16 10 (58%) 14 (47%) 20 (66%) 44 (57%)
HPV18 1 (6%) 6 (20%) 2 (7%) 9 (12%)
HPV31 - 2 (7%) 2 (7%) 4 (5%)
HPV33 1 (6%) - - 1 (1%)
HPV53 1 (6%) 1 (3%) - 2 (3%)
HPV58 2 (12%) 4 (13%) 4 (13%) 10 (13%)
HPV70 1 (6%) - - 1 (1%)
Multiple infections - 3 (10%) 2 (7%) 5 (7%)

HSIL: high-grade squamous intraepithelial lesions; LSIL: low-grade squamous intraepithelial lesions.

Of the 65 HPV16-positive samples, 49 (75%) were classified into lineages according to signatures proposed by Cornet et al. (2012). Forty-three samples (87%) were found to belong to lineage A (sublineages A1, A2 or A3). Only one (2%) analysed sequence was identified as lineage B, and lineage C was not identified in this study. Five variants (10%) were classified as lineage D (Table IV). For this analysis, we considered all SNPs in the E6 region that were proposed as criteria to classify lineages. However, in the LCR region, only two SNPs (7764C and 7786C) were used to classify HPV sequences into lineage A, and three SNPs (A7233C, A7485C, G7489A) were used to classify HPV sequences into lineage D because the amplified LCR fragments did not cover all lineage-defining SNPs in this genomic region. We failed to PCR amplify and sequence the E6 and LCR regions for three and four samples, respectively.

TABLE IV HPV16 lineage distribution in cervical samples according to histopathology 

HPV16 lineage Normal histology (N = 21) LSIL (N = 12) HSIL (N = 16) Total (N = 49)
A 19 (90%) 10 (84%) 14 (88%) 43 (88%)
B - 1 (8%) - 1 (2%)
D 2 (10%) 1 (8%) 2 (12%) 5 (10%)

HSIL: high-grade squamous intraepithelial lesions; LSIL: low-grade squamous intraepithelial lesions.

Regarding demographic and clinical data for women included in the study, HPV16 lineage A was the most prevalent in women both with and without SIL. Among those infected with lineage A HPV strains, 74% of women were white and 26% were black. Among those infected with lineage D, the second most prevalent lineage in groups both with and without SIL, 60% were black and 40% were white. HPV16 lineage B was only found in one white patient with LSIL (Table IV).

Among the HPV16 lineage A-infected participants, four women carried the T350G (L83V) polymorphism in E6, and all lesions were classified as HSIL.


In the present study, we determined the prevalence and genotypes of HPV, as well as HPV16 genetic variants, in HIV-negative and HIV-positive women in a major city of southern Brazil. This is the first study to describe the intratypic HPV distribution in women with and without intraepithelial lesions in this geographic area of the country.

HPV DNA was detected in 24% of the cervical smear samples analysed. A similar prevalence was found in other studies conducted in Brazil: in the north, 25.3% (Vieira et al. 2015); in the northeast, 24.5% (Fernandes et al. 2010); and in the south, 24.6% (Rosa et al. 2008). Cervical biopsies revealed the presence of HPV in all women with LSIL or HSIL. This relationship between HPV infection, cervical SIL, and cervical cancer has been well documented worldwide, since the presence of HPV is a necessary condition for the development of this neoplasm (Walboomers et al. 1999).

The biopsies analysed were from women preliminarily diagnosed with HSIL by cytological examination, and a high prevalence of high-risk HPV types was determined by PCR. However, one-third of these women were SIL-negative by histopathological analysis, showing a large difference in the specificity of these two techniques as conducted in our setting. Thus, the diagnostic triad (cytology, colposcopy, and histopathology), in association with molecular biology techniques, can be used for early identification of patients at high risk of developing cervical cancer. In fact, the American College of Obstetrics and Gynecology (ACOG) recommends the use of combined HPV tests in women over 30 years of age (ACOG 2012).

In this study, HPV16 was the most prevalent type, followed by HPV58, in women both with and without SIL. HPV16 causes approximately 70% of cervical cancers around the world and is the most prevalent type detected in samples from women with normal cytology and either low- or high-grade SIL (de Sanjosé et al. 2010, Bruni et al. 2015). The prevalence of other HPV types varies widely by geographic origin of the study population. The six HPV types most prevalent in South America are HPV16, 18, 31, 45, 33, and 58. HPV58 is the third most prevalent type in Brazil among women with either low- or high-grade SIL (ICO Information Centre on HPV and Cancer) (Li et al. 2011). In a study in Natal (northeastern Brazil), Fernandes et al. (2008) showed that HPV16 was the most prevalent genotype and that HPV58 was the second most prevalent type in women with normal cytology and LSIL; however, HPV58 showed a prevalence equal to that of HPV18 in patients with HSIL.

In addition to HPV16 and 58, types 31, 33, and 45 were also prevalent in this study in women either with or without SIL. Li et al. (2011) and de Sanjosé et al. (2010) showed through two meta-analyses that HPVs 31, 33, 45, and 58 are among the six typed most frequently reported to be associated with cervical cancer worldwide, after HPV16 and 18, but these types are not yet included in the bivalent or quadrivalent vaccines currently available in the developing world. Further studies are needed to clarify whether the eradication of oncogenic HPV16 and 18 with current vaccines can generate selective pressure, thereby increasing the prevalence of minority HPV types associated with carcinogenesis (Haug 2008). This study found a lower prevalence of HPV18 than that found in another study carried out in southern Brazil, which identified this type as the second most prevalent, followed by HPV16 (Entiauspe et al. 2014).

The association between HPV and HIV infections is well documented (Houlihan et al. 2012), with the prevalence of HPV infection higher in HIV-positive women (Corrêa et al. 2011, Menon et al. 2016). In this study, the frequency of HPV infection in HIV-positive women was also higher than that in HIV-negative women. A meta-analysis of data from various countries in North America, Africa, Asia, Europe, and South-Central America determined the prevalence and genotypes of HPV in 5,578 HIV-positive women. HPV16 was the most prevalent type in this population (Clifford at al. 2006). In Brazil, HPV16 is also the most common type detected in HIV-positive women (Corrêa et al. 2011), as we found in the current study.

In this study, the most prevalent variants of HPV16 were classified as lineage A, followed by lineages D and B. Ho et al. (1993) showed that HPV16 lineages A, B, and D are the most prevalent in populations of European, African, and Amerindian ancestry, respectively. In southern Brazil (Rio Grande do Sul), the site of this study, the population is largely of European descent, with less African and Amerindian ancestry (Manta et al. 2013). A higher prevalence of lineage A was expected and observed in our study.

In a recent meta-analysis of the worldwide distribution of HPV16 variants, lineage A was found to represent the majority of isolates and was common in all regions, except for sub-Saharan Africa and East Asia. Lineage D was also detected in all regions except for sub-Saharan Africa, whereas lineages B/C prevailed in North and Sub-Saharan Africa (Cornet et al. 2013). In Brazil, the high frequency of lineage A was also observed in women without cervical SIL in northern and southeastern Brazil (Tamegão-Lopes et al. 2014). The results by Ortiz-Ortiz et al. (2015) in southern Mexico in women with normal cytology, LSIL, HSIL, and cancer were similar to the results presented here, with an 82.1% prevalence for lineage A, 17.5% for lineage D, and 0.3% for lineage B.

HPV16 lineage D is reported to be the second most prevalent lineage after A in women with cervical SIL and invasive cancer (Vidal et al. 2016). Burk et al. (2003) reported that non-European variants were associated with a higher risk of developing cervical cancer, although much of this association was specifically related to an Asian-American variant (lineage D). In a recent study conducted by Mirabello et al. (2016), the D2/D3 sublineage was associated with increased risk of cervical intraepithelial neoplasia grade 3 and cervical cancer. In particular, D2 was associated with increased risk of glandular lesions, adenocarcinomas, and adenocarcinomas in situ.

The T350G variant of HPV16 lineage A was detected in the present study in women with LSIL and in women with HSIL. An earlier study suggests that this polymorphism may affect viral persistence and risk for progression to cervical lesions (Gheit et al. 2011). The E350G variant isolated in South/Central America has been associated with cervical cancer (Cornet et al. 2013), and one study reported that it has a higher capacity to transform human keratinocytes in vitro than the E350T variant (Sichero et al. 2012).

An important limitation of the present study was the inability to identify HPV types in co-infections. In five samples with indications of multiple infections, it was not possible to identify more than one HPV type because of the low sensitivity of the sequencing technique used (Sanger), as described by da Fonseca et al. (2016).

In conclusion, this study identified highly diverse HPV oncogenic types (HPV 16, 18, 31, 33, 45, and 82) and various HPV16 lineages in southern Brazil, which agrees with expected types based on the genetic background of this population. Studies on the prevalence of different HPV types in various regions of the world are essential to estimating the potential impact of vaccines to prevent cervical cancer currently available in the developing world and on the design of screening programs. For vaccine development, it is necessary to consider the distribution of HPV genotypes, as well as HPV intratype variants, in Brazil and worldwide, because genetic differences can determine the severity and outcome of HPV infections.


To the Service of Gynecology and Obstetrics of the University Hospital, Federal University of Rio Grande (HU-FURG), and to the Policlínica de Assistência Médica of Rio Grande city, for collecting and allowing an analysis of the samples used in this study. We also want to thank the Genetics Program of the Brazilian National Cancer Institute (INCA), for helping us with sequencing analyses.


ACOG - The American College of Obstetricians and Gynecologists. ACOG practice bulletin. Clinical management guidelines for obstetrician-gynecologists: screening for cervical cancer. 2012. Available from: [ Links ]

Bernard HU, Burk RD, Chen Z, Vandoorslaer K, Hausen H, Devilliers EM. Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments. Virology. 2010; 401(1): 70-9. [ Links ]

Bernard HU. Gene expression of genital human papillomavirus and considerations on potential antiviral approaches. Antivir Ther. 2002; 7(4): 219-37. [ Links ]

Bruni L, Barrionuevo-Rosas L, Albero G, Aldea M, Serrano B, Valencia S, et al. Human papillomavirus and related diseases in the world. Summary Report 2015. ICO Information Centre on HPV and Cancer (HPV Information Centre). [ Links ]

Burk RD, Harari A, Chen Z. Human papillomavirus genome variants. Virology. 2013; 445(1): 232-43. [ Links ]

Burk RD, Teraii M, Gravitt PE, Brinton LA, Kurman RJ, Barnes WA, et al. Distribution of human papillomavirus types 16 and 18 variants in squamous cell carcinomas and adenocarcinomas of the cervix. Cancer Res. 2003; 63(21): 7215-20. [ Links ]

Clifford GM, Gonçalves MA, Franceschi S. Human papillomavirus types among women infected with HIV: a meta-analysis. AIDS. 2006; 20(18): 2337-44. [ Links ]

Cornet I, Gheit T, Franceschi S, Vignat J, Burk RD, Sylla, BS, et al. Human papillomavirus type 16 genetic variants: phylogeny and classification based on E6 and LCR. J Virol. 2012; 86(12): 6855-61. [ Links ]

Cornet I, Gheit T, Liannacone MR, Vignat J, Sylla BS, del Mistro A, et al. HPV16 genetic variation and the development of cervical cancer worldwide. Br J Cancer. 2013; 108(1): 240-4. [ Links ]

Corrêa CM, Teixeira NCP, de Araújo ACL, Carvalho NO, Del Castillo DM, Campos RR, et al. Prevalence and multiplicity of HPV in HIV women in Minas Gerais, Brazil. Rev da Assoc Med Bras. 2011; 57(4): 425-30. [ Links ]

da Fonseca AJ, Galvão RS, Miranda AE, Ferreira LCDL, Chen Z. Comparison of three human papillomavirus DNA detection methods: next generation sequencing, multiplex-PCR and nested-PCR followed by sanger based sequencing. J Med Virol. 2016; 88(5): 888-94. [ Links ]

de Sanjosé S, Quint WG, Alemany L, Geraets DT, Klaustermeier JE, Lloveras B, et al. Retrospective international surveyand HPV time trends study group human papillomavirus genotype attribution in invasive cervical cancer: a retrospective cross-sectional worldwide study. Lancet Oncol. 2010; 11(11): 1048-56. [ Links ]

de Villiers EM, Fauquet C, Broker TR, Bernard HU, zur Hausen H. Classification of papillomaviruses. Virology. 2004; 324(1): 17-27. [ Links ]

Entiauspe LG, Silveira M, Nunes EM, Basgalupp SP, Stauffert D, Dellagostini OA, et al. High incidence of oncogenic HPV genotypes found in women from southern Brazil. Braz J Microbiol. 2014; 45(2): 689-94. [ Links ]

Fernandes JV, Meissner RV, Carvalho MG, Fernandes TA, Azevedo PR, Sobrinho JS, et al. 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(1): 96. [ Links ]

Fernandes TAADM, Meissner RDV, Bezerra LF, Azevedo PRMD, Fernandes JV. Human papillomavirus infection in women attended at a cervical cancer screening service in Natal, Brazil. Braz J Microbiol. 2008; 39(3): 573-8. [ Links ]

Gheit T, Cornet I, Clifford GM, Iftner T, Munk C, Tommasino M, et al. Risks for persistence and progression by human papillomavirus type 16 variant lineages among a population-based sample of Danish women. Cancer Epidemiol Biomarkers Prev. 2011; 20(7): 1315-21. [ Links ]

Globocan. Estimated cancer incidence, mortality and prevalence worldwide in 2012. Updated 2017. Available from: [ Links ]

Haug CJ. Human papillomavirus vaccination - reasons for caution. N Engl J Med. 2008; 359(8): 861-2. [ Links ]

Ho L, Chan SY, Burk RD, Das BC, Fujinaga K, Icenogle JP, et al. The genetic drift of human papillomavirus type 16 is a means of reconstructing prehistoric viral spread and the movement of ancient human populations. J Virol. 1993; 67(11): 6413-23. [ Links ]

Houlihan CF, Larke NL, Watson-Jones D, Smith-McCune KK, Shiboski S, Gravitt PE, et al. Human papillomavirus infection and increased risk of HIV acquisition: a systematic review and meta-analysis. AIDS. 2012; 26(17): 2211-22. [ Links ]

INCA - Instituto Nacional de Câncer. Estimativa 2016: incidência de câncer no Brasil. 2016. Available from: [ Links ]

Li N, Franceschi S, Howell-Jones R, Snijders PJ, Clifford GM. Human papillomavirus type distribution in 30,848 invasive cervical cancers worldwide: variation by geographical region, histological type and year of publication. Int J Cancer. 2011; 128(4): 927-35. [ Links ]

Manos MM, Ting Y, Wright DK, Lewis AJ, Broker TR, Wolinsky SM. Use of polymerase chain reaction amplification for the detection of genital human papillomaviruses. Cancer Cells. 1989; 7(17): 209-14. [ Links ]

Manta FS, Pereira R, Vianna R, de Araújo ARB, Gitaí DLG, da Silva DA, et al. Revisiting the genetic nacestry of Brazilians using autosomal AIM-Indels. PLoS ONE. 2013; 8(9): 75145. [ Links ]

Menon S, Wusiman A, Boily MC, Kariisa M, Mabeya H, Luchters S, et al. Epidemiology of HPV genotypes among HIV positive women in Kenya: a systematic review and meta-analysis. PLoS ONE. 2016; 11(10): 0163965. [ Links ]

Mirabello L, Yeager M, Cullen M, Boland JF, Chen Z, Wentzensen N, et al. HPV16 sublineage associations with histology-specific cancer risk using HPV whole-genome sequences in 3200 women. J Natl Cancer Inst. 2016; 108(9): 1-9. [ Links ]

Muñoz N, Bosch FX, de Sanjosé S, Herrero R, Castellsagué X, Shah KV, et al. Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003; 348(6): 518-27. [ Links ]

Ortiz-Ortiz J, Alarcón-Romero LC, Jiménez-López MA, Garzón-Barrientos VH, Calleja-Macías I, Barrera-Saldaña HA, et al. Association of human papillomavirus 16 E6 variants with cervical carcinoma and precursor lesions in women from southern Mexico. J Virol. 2015; 12(1): 29. [ Links ]

Rosa MI, Fachel JM, Rosa DD, Medeiros LR, Igansi CN, Bozzetti MC. Persistence and clearance of human papillomavirus infection: prospective cohort study. Am J Obstet Gynecol. 2008; 199(6): 1-7. [ Links ]

Sichero L, Sobrinho JSS, Villa Ll. Oncogenic potential diverge among human papillomavirus type 16 natural variants. Virology. 2012; 432(1): 127-32. [ Links ]

Tamegão-Lopes BP, Sousa-Júnior EC, Passetti F, Ferreira CG, de Mello WA, Silvestre RVD. Prevalence of human papillomavirus infection and phylogenetic analysis of HPV-16 E6 variants among infected women from northern Brazil. Infect Agent Cancer. 2014; 9(1): 25. [ Links ]

Van den Brule AJ, Snijders PJ, Gordijn RL, Bleker OP, Meijer CJ, Walboomers JM. General primer-mediated polymerase chain reaction permits the detection of sequenced and still unsequenced human papillomavirus genotypes in cervical scrapes and carcinomas. Int J Cancer. 1990; 45(4): 644-9. [ Links ]

Vidal JPCB, Felix SP, Chaves CBP, Patury P, Franco VF, de Morais EA, et al. Genetic diversity of HPV16 and HPV18 in Brazilian patients with invasive cervical cancer. J Med Virol. 2016; 88(7): 1279-87. [ Links ]

Vieira RC, Monteiro JDSV, Manso EP, dos Santos MRM, Tsutsumi MY, Ishikawa EAY, et al. Prevalence of type-specific HPV among female university students from northern Brazil. Infect Agent Cancer. 2015; 10(1): 21. [ Links ]

Villa LL, Sichero L, Rahal P, Caballero O, Ferenczy A, Rohan T, et al. Molecular variants of human papillomavirus types 16 and 18 preferentially associated with cervical neoplasia. J Gen Virol. 2000; 81(12): 2959-68. [ Links ]

Walboomers MM, Acos MV, Manos MM, Bosch FX, Kummer JA. Human papillomavirus is a necessary cause of invasive cervical cancer. J Pathol. 1999; 189(1): 12-9. [ Links ]

Yamada T, Manos MM, Peto J, Greer CE, Muñoz N, Bosch FX, et al. Human papillomavirus type 16 sequence variation in cervical cancers: a worldwide perspective. J Virol. 1997; 71(3): 2463-72. [ Links ]

Zhang L, Liao H, Yang B, Geffre CP, Zhang A, Zhou A, et al. Variants of human papillomavirus type 16 predispose toward persistent infection. Int J Clin Exp Pathol. 2015; 8(7): 8453-9. [ Links ]

zur Hausen H. Condylomata acuminata and human genital cancer. Cancer Res. 1976; 36(2): 794. [ Links ]

Financial support: CAPES (552361/2011-9), CNPq (305765/2015-9), FAPERJ (201.331/2014).

Received: December 12, 2016; Accepted: March 4, 2017

+ Corresponding author:


GRO, ECA and SGO - Organised and collected samples; GRO, ECA, FF-J, TDMC and FAAG - performed all molecular biological assays; GRO and VCV - performed the data analyses; CVG - supervised all sample collections and clinical evaluation of the study subjects; VPH, MAS and AMBM - conceived,supervised and provided infrastructure for the entire study; GRO, VCV, FF-J, MAS and AMBM - wrote the manuscript. All authors read and agreed with the contents and submission of this manuscript.

Creative Commons License  This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.