Lineage analysis of human papillomavirus type 39 in cervical samples of Iranian women
Virology Journal volume 18, Article number: 152 (2021)
The data with regards to the regional variants of distinct HPV types is of great value. Accordance with this, this study aimed to investigate the sequence variations of E6 gene and long control region of HPV 39 among normal, premalignant and malignant cervical samples in order to characterize the frequent HPV 39 variants circulating in Tehran, Iran.
In total, 70 cervical samples (45 normal, 16 premalignant, and 9 malignant samples) infected with HPV 39 were analyzed by nested-PCR and sequencing.
Our results revealed that all samples belonged to A lineage. Almost all sequences (98.6%) were classified in A1 sublineage and only one sample (1.4%) was A2 sub lineage.
Our findings showed that lineages A, sublineage A1, is dominant in Tehran, Iran. However, the small sample size was the most important limitations of this study. Further studies with larger sample size from different geographical regions of Iran are necessary to estimate the pathogenicity risk of HPV 39 variants in this population.
Cervical cancer is reported as the fourth most frequent cancer among women worldwide . Human papillomavirus (HPV) is considered as the etiological agent of cervical cancer. Although more than 40 HPV types can infect the anogenital tract, 14 HPV types, designated as the high-risk HPV types including 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 68 and 73, associate with progression to cervical cancer [2,3,4,5,6]. Among high-risk HPV types, HPV 16 and 18 are leading cause and responsible for almost 72% of cervical cancer . Most studies on the prevalence of HPV types has indicated that HPV 16 and 18 are the most prevalent types in women with normal cervical cytology and invasive cervical cancer in Iran [8,9,10,11,12,13]. A meta-analysis was revealed that the five most common HPV types are HPV 16, 18, 31, 39, and 45 in Iranian women with cervical cancer . Indeed, HPV 16, 18, 31, 39, and 45 are account for 52%, 14.1%, 5%, 3.3%. and 3.3% of cervical cancer, respectively .
A distinct HPV type is considered when the DNA sequence of the L1 gene was differ from that of any other characterized type by at least 10%. Isolates of the same HPV type are designated to as lineage and sublineage when the nucleotide sequences of the L1 gene differ by less than 1–10% and 0.5–1%, respectively. Up to date, HPV 39 has two distinct lineages A and B; which lineage A include 2 different sublineages A1 and A2 .
While the co-evolution of HPV 16 and HPV 18 variants and human populations is well characterized [4, 16,17,18,19], the geographic associations for variants of other types remains to be inconclusive. To better understand of an association between the distinct variants of HPV 39 and ethnicity further studies are needed in the world.
The data with regards to the regional variants of distinct HPV types is of great value as it would provide a rational for future studies on their evolution, epidemiology, pathogenicity, and biology. Whereas the spreading of HPV types is identified in Iran, there is much less known about HPV variants. Previous studies were investigated the common HPV 16, 18, 31, and 45 variants in Iran [19,20,21]. Regard to above-mentioned data that HPV 39 is one of the five most common HPV types in Iranian women with cervical cancer , this study aimed to investigate the frequent HPV 39 variants circulating in Iran. According to this fact that E6 gene can only separate A and B lineages from each other, long control region (LCR) were analyzed to distinguish two distinctive sublineages A1 and A2.
Material and methods
To characterize lineages and sublineages of HPV 39 a cross-sectional study was designed from 2018 to 2020. One-hundred and fifty-six formalin-fixed paraffin-embedded (FFPE) samples (98 invasive cervical cancer samples and 58 premalignant cervical lesions) were obtained from Immam-Khomeini hospital in Tehran. Moreover, 55 ThinPrep Pap Test specimens (10 premalignant and 45 normal samples) that were HPV 39 positive were collected from three referral laboratories in Tehran from 2018 to 2020. These samples were previously genotyped by Cobas assay or INNO-LiPA® HPV Genotyping assay. All of study subjects sign the informed consent and all information regard to their medical records were stored on secure system.
Variant analysis of HPV 39 based on E6 gene and long control region (LCR)
DNA extraction from ThinPrep Pap Test specimens was carried out by the High Pure Viral Nucleic Acid Kit (Roche Diagnostics GmbH, Roche Applied Science, Mannheim, Germany) according to the manufacturer’s instruction. Genomic DNA from formalin-fixed paraffin-embedded (FFPE) samples was isolated using phenol–chloroform assay according to previously published procedure .
Detection of HPV genome was investigated in 156 FFPE samples using nested-PCR with MY09/MY11 and GP+5/GP+6 primer pairs, respectively, to amplify a 150 bp fragment of L1 gene. To genotype, all HPV positive samples were subjected to sequencing by BigDye® Terminator v3.1 Cycle Sequencing Kit and a 3130 Genetic Analyzer Automated Sequencer as specified by Applied Biosystems manuals (Foster City, CA). Nucleotide sequences were edited using Bioedit software and converted to FASTA format. Finally, edited sequences were blasted using the Blast server (http://www.ncbi.nlm.nih.gov/blast/) to find HPV genotypes. All of HPV positive samples that were genotyped by sequencing also screened by nested-PCR with sequence-specific primers of HPV 39 to amplify of E6 gene (Table 1).
The complete sequence of HPV 39 E6 gene (nucleotide [nt] 107–583) was investigated using nested-PCR with the primer pairs were shown in Table 1. The PCR reaction was achieved in a 50 μl reaction mixture including 100–200 ng of DNA template, 2 mM MgCl2, 50 μM of each dNTP, 10 pmol of each primer, and 2 U of Taq DNA polymerase. PCR amplification cycles were as follow for both round: 35 cycles of 95 °C for 20 s, 55 °C for 40 s and 72 °C for 50 s. A reaction mixture without template DNA, as a negative control, was included in every set of PCR run (Additional file 1: Figure S1).
The LCR region of HPV 39 genome was amplified by nested-PCR with primer pairs 39-LCR-F1 and 39-LCR-R1 for first round and 39-LCR-F2 and 39-LCR-R2 for the second round (Table 1). The PCR reaction was carried out in a 50 μl reaction mixture including 100–200 ng of DNA template, 50 μM of each dNTP, 1.5 mM MgCl2, 10 pmol of each primer, and 2 U of Taq DNA polymerase. PCR amplification cycles were as follow for both round: 35 cycles of 95 °C for 20 s, 56 °C for 50 s and 72 °C for 50 s. A reaction mixture without template DNA, as a negative control, was included in every set of PCR run (Additional file 2: Figure S2).
To analyze of HPV 39 variants, the PCR products were addressed to bidirectional direct sequencing as described above. All sequences obtained in this study are available at http://www.ncbi. nlm.nih.gov/ with accession numbers MW390245-MW390314 and MW390315-MW390384 for E6 gene and LCR, respectively. To identify lineages and sublieages of the studied samples, sequences were aligned to reference sequences with accession number M62849, KC470239, and KC470247 designated for A1, A2, and B, respectively . Phylogenetic tree was constructed using the maximum likelihood method by Mega software version 6 . The reliability was measured by calculation of bootstrap with 1000 replicates. The reference sequences of A1, A2, and B  were also retrieved from GenBank database.
The statistical analysis was performed by Mantel–Haenszel test (Epi Info 7, Statistical Analysis System Software) and the p-value less than 0.05 was considered statistically significant.
HPV 39 was detected in 15 FFPE samples (6 premalignant and 9 malignant samples). In total, 70 HPV 39-positive cervical samples of normal (n = 45), premalignant (n = 16), and malignant (n = 9) specimens were investigated in this study. The mean age of studied groups was 51.5, 37.5, and 30 years old in malignant, premalignant, and normal samples, respectively.
The complete E6 gene (nt 107-583) of 70 HPV 39-positive samples were sequenced and compared with the HPV 39 E6 prototype sequence (GeneBank accession number M62849). As indicated in Table 2, nucleotide substitutions at positions of G78T, G151T, C260G, T305C, and A453G are specific for lineage B. Phylogenetic tree analysis of HPV E6 gene was conducted in MEGA6 by using the Maximum Likelihood method based on the Kimura 2-parameter model . Phylogenetic tree analysis of 70 samples were indicated that no nucleotide substitutions were happened in the E6 gene and all samples were belonged to lineage A (Fig. 1). Consequently, no amino acid changes were occurred in our samples (Table 2).
To distinguish sublineages A1 and A2 from each other, the partial sequence of LCR was investigated. Regard to LCR analysis, five nucleotide substitutions were found in only one of our studied samples at positions of C7203A, G7398T, T7438G, A7571C, and A7702G. Phylogenetic tree analysis of HPV LCR was conducted in MEGA6 by using the Maximum Likelihood method based on the Kimura 2-parameter model . Phylogenetic tree analysis was revealed that almost all sequences (98.6%) were classified in sublineage A1 and only one sample (1.4%) was belonged to sublineage A2 (Fig. 2, Table 2).
As shown in Table 1 and Fig. 2, all of premalignant and malignant cervical samples were classified in sublineage A1. In normal group, 44 out of 45 samples (97.8%) were belonged to sublineage A1 and one samples (2.2%) was classified in sublineage A2. However, no statistically significant differences were found in this regard among three studied groups (P > 0.05).
As indicated in Table 3, 50 out of 70 samples (71.4%) were found to be infected with at least two different HPV types. In most of multiple HPV-infected samples (64%), another high-risk HPV types were detected and in 36% of samples a low-risk HPV type was found. Regard to histological type, multiple HPV infections were detected in 77.8%, 50%, and 77.8% of malignant, premalignant, and normal cases, respectively. Interestingly, infection with at least a high-risk HPV type in studied groups were found to be 85.7%, 50%, and 62.8% in malignant, premalignant, and normal samples, respectively (Table 3).
It is well recognized that the distribution of HPV 16 and 18 variants can be different geographically due to the evolution associated to the population ethnicity [16, 24,25,26,27]. However, the geographic relations for variants of other types remain inconclusive. In this regard, this study investigated the sequence variations of E6 gene and LCR of HPV 39 in order to find the frequent HPV 39 variants in normal samples as well as premalignant and malignant lesions of the cervix from Iranian women.
Almost uniform distribution of HPV 39 variants was found in normal and premalignant/malignant groups, as lineage A was dominant. Our result is consistent with few studies that conducted in some geographical regions of the world. One study in the United States of America was investigated the HPV 39 lineages in 479 samples as the most of samples belonged to lineage A (99%) and B lineage was only detected in 1% . Another study shown that although lineage A was dominant in Costa Rica (67% A1 and 33% A2 sublineages), only lineage B was found in Burkina Faso . In women from Shanghai, indicated that all three sublineages of HPV 39 are present as follow: sublineage A1 was dominant (81.2%) and the remaining of cases were belonged to sublineages A2 and B in 12.5% and 6.3%, respectively .
It is suggested that like HPV 16-related alpha-9 species variants, isolates of alpha-7 species have shown co-divergence of archaic Hominid and HPV variants. It seems that rapid expansion of the host population causing viral variant lineages and sublineages that had occurred almost 0.2–1.0 million years ago . Indeed, the distinction coevolution of HPV 16 variants with closely related ancestral human populations together with introgression of certain archaic alleles related to innate immunity and keratinocytes differentiation into the human genomes might let to the dominancy of HPV 16 A lineage in modern human ancestor population .
Viral variants of each high-risk HPV types show differences in infectivity, long-term persistence, development of precancer lesions, and invasion of HPV containing cell . Indeed, it is indicated that distinct HPV 16 lineages have different risks for development of cervical cancer. Several studies found that HPV 16 variants of non-European lineages (B, C, D) have stronger oncogenic potential than the European lineage A [18, 32,33,34]. In this study, however, no statistically significant differences were found by distinct lineages or sublineages of HPV 39 in three studied groups. Concordance with our result, finding of one study was revealed that no difference was observed in the increased risk of CIN2/3 regard to HPV 39 sublineages . Moreover, regard to this point that multiple infection with another high-risk HPV types, particularly HPV 16 and 18 was common in malignant samples; it seems that probably HPV 16 or 18 played an important role in development of cervical cancer in Iranian women rather than HPV 39.
An association between population-based oncogenicity of HPV 39 variants and host genetic variations, particularly HLA class I and II alleles, could be assumed. Although in this regard no data was found for HPV 39, it is indicated that the E-G350 variant of HPV 16 had 4–fivefold higher risk for cancer development among Swedish women with HLA-B*44, HLA-B*51, or HLA-B*57 alleles . A significant link between different E6 variants of HPV 16 and HLA class II alleles was also found for three distinct alleles (DRB1*1501, DRB1*1502 and DQB1*0602) in Japanese women, as DRB1*1501 and DQB1*0602 alleles were significantly higher among patients with the prototype variant while DRB1*1502 was increased in patients with D25E variants compared with controls .
The most important limitations of this study were found to be the small sample size and multiple infections.
Our findings showed that lineages A, sublineage A1, is dominant in Tehran, Iran. HPV variant studies have been made important data that could help better understanding the carcinogenicity of different HPV types at molecular level. However, the small sample size and multiple infections were the most important limitations of this study. Further studies with larger sample size of single infection with HPV 39 from different geographical regions of Iran are necessary to estimate the pathogenicity risk of HPV 39 variants in this population. The characterization of HLA in Iranian women with cervical cancer is highly recommended as it may shed a new light to understand a link between HPVs variants and genetic background.
Availability of data and materials
Data available within the article.
Long control region
(2020) International Agency for Research on Cancer. GLOBOCAN 2020: estimated cancer incidence, mortality and prevalence worldwide in 2020.
Bouvard V, Baan R, Straif K, Grosse Y, Secretan B, El GF, Brahim-Tallaa L, Guha N, Freeman C, Galichet L, Cogliano V. A review of human carcinogens—part B: biological agents. Lancet Oncol. 2009;10:321–2.
Doorbar J, Quint W, Banks L, Bravo IG, Stoler M, Broker TR, Stanley MA. The biology and life-cycle of human papillomaviruses. Vaccine. 2012;30:F55–70.
Arroyo SL, Basaras M, Arrese E, Hernaez S, Andia D, Esteban V, Garcia-Etxebarria K, Jugo BM, Cisterna R. Human Papillomavirus (HPV) genotype 18 variants in patients with clinical manifestations of HPV related infections in Bilbao Spain. Virol J. 2012;9:258.
Okunade KS. Human papillomavirus and cervical cancer. J Obst Gynaecol J Inst Obst Gynaecol. 2020;40:602–8.
Béziat V. Human genetic dissection of papillomavirus-driven diseases: new insight into their pathogenesis. Hum Genet. 2020;139:919–39.
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:927–35.
Hamkar R, Shoja Z, Ghavami N, Heydari N, Farahmand M, Jalilvand S. Type-specific human papillomavirus prevalence in Iranian women with normal cervical cytology: the impact of current HPV vaccines. Intervirology. 2017;60:125–30.
Heydari N, Oskouee MA, Vaezi T, Shoja Z, Esmaeili HA, Hamkar R, Shahmahmoodi S, Jalilvand S. Type-specific human papillomavirus prevalence in cervical intraepithelial neoplasia and cancer in Iran. J Med Virol. 2018;90:172–6.
Hamkar R, Azad TM, Mahmoodi M, Seyedirashti S, Severini A, Nategh R. Prevalence of human papillomavirus in Mazandaran province, Islamic Republic of Iran. 2002.
Hasanzadeh M, Amel Jamehdar S, Malakuti H, Bolandi S, Shakeri M, Harifi E, Radmanesh HODA. The prevalence of fourteen human papilloma virus types in women of Khorasan Razavi Province, Northeast Iran. Int J Gynecol Cancer. 2015;25:595.
Salehi-Vaziri M, Sadeghi F, Hashemi FS, Haeri H, Bokharaei-Salim F, Monavari SH, Keyvani H. Distribution of human papillomavirus genotypes in Iranian women according to the severity of the cervical lesion. Iran Red Crescent Med J. 2016; 18.
Khodakarami N, Clifford GM, Yavari P, Farzaneh F, Salehpour S, Broutet N, Bathija H, Heideman DA, van Kemenade FJ, Meijer CJ, Hosseini SJ, Franceschi S. Human papillomavirus infection in women with and without cervical cancer in Tehran, Iran. Int J Cancer. 2012;131:E156-161.
Jalilvand S, Shoja Z, Nourijelyani K, Tohidi HR, Hamkar R. Meta-analysis of type-specific human papillomavirus prevalence in Iranian women with normal cytology, precancerous cervical lesions and invasive cervical cancer: Implications for screening and vaccination. J Med Virol. 2015;87:287–95.
Burk RD, Harari A, Chen Z. Human papillomavirus genome variants. Virology. 2013;445:232–43.
Xi LF, Kiviat NB, Hildesheim A, Galloway DA, Wheeler CM, Ho J, Koutsky LA. Human papillomavirus type 16 and 18 variants: race-related distribution and persistence. J Natl Cancer Inst. 2006;98:1045–52.
Burk RD, Terai M, Gravitt PE, Brinton LA, Kurman RJ, Barnes WA, Greenberg MD, Hadjimichael OC, Fu L, McGowan L. Distribution of human papillomavirus types 16 and 18 variants in squamous cell carcinomas and adenocarcinomas of the cervix. Can Res. 2003;63:7215–20.
Freitas LB, Chen Z, Muqui EF, Boldrini NA, Miranda AE, Spano LC, Burk RD. Human papillomavirus 16 non-European variants are preferentially associated with high-grade cervical lesions. PLoS ONE. 2014;9:100746.
Vaezi T, Shoja Z, Hamkar R, Shahmahmoodi S, Nozarian Z, Marashi SM, Jalilvand S. Human papillomavirus type 16 lineage analysis based on E6 region in cervical samples of Iranian women. Infect Genet Evol J Mol Epidemiol Evolut Genet Infect Dis. 2017;55:26–30.
Salavatiha Z, Shoja Z, Heydari N, Marashi SM, Younesi S, Nozarian Z, Jalilvand S. Lineage analysis of human papillomavirus type 18 based on E6 region in cervical samples of Iranian women. J Med Virol. 2020.
Hosseini N, Shoja Z, Younesi S, Shafiei-Jandaghi NZ, Jalilvand S. Lineage analysis of human papillomavirus types 31 and 45 in cervical samples of Iranian women. J Med Virol. 2021;93:3857–64.
Jalilvand S, Tornesello ML, Buonaguro FM, Buonaguro L, Naraghi ZS, Shoja Z, Ziaee AA, Hamkar R, Shahmahmoodi S, Nategh R, Mokhtari-Azad T. Molecular epidemiology of human herpesvirus 8 variants in Kaposi’s sarcoma from Iranian patients. Virus Res. 2012;163:644–9.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Cornet I, Gheit T, Iannacone MR, Vignat J, Sylla BS, Del Mistro A, Franceschi S, Tommasino M, Clifford GM. HPV16 genetic variation and the development of cervical cancer worldwide. Br J Cancer. 2013;108:240.
Zehbe I, Mytilineos J, Wikstrom I, Henriksen R, Edler L, Tommasino M. Association between human papillomavirus 16 E6 variants and human leukocyte antigen class I polymorphism in cervical cancer of Swedish women. Hum Immunol. 2003;64:538–42.
Chen AA, Gheit T, Franceschi S, Tommasino M, Clifford GM, Group IHVS. Human papillomavirus 18 genetic variation and cervical cancer risk worldwide. J Virol. 2015;89:10680–7.
Perez S, Cid A, Inarrea A, Pato M, Lamas MJ, Couso B, Gil M, Alvarez MJ, Rey S, Lopez-Miragaya I, Melon S, Ona M. Prevalence of HPV 16 and HPV 18 lineages in Galicia, Spain. PLoS ONE. 2014;9:e104678.
Xi LF, Schiffman M, Koutsky LA, Hughes JP, Winer RL, Mao C, Hulbert A, Lee SK, Shen Z, Kiviat NB. Lineages of oncogenic human papillomavirus types other than type 16 and 18 and risk for cervical intraepithelial neoplasia. J Natl Cancer Inst. 2014; 106.
Chen Z, Schiffman M, Herrero R, DeSalle R, Anastos K, Segondy M, Sahasrabuddhe VV, Gravitt PE, Hsing AW, Burk RD. Evolution and taxonomic classification of alphapapillomavirus 7 complete genomes: HPV18, HPV39, HPV45, HPV59, HPV68 and HPV70. PLoS ONE. 2013;8:e72565.
Zhao J, Zhan Q, Guo J, Liu M, Ruan Y, Zhu T, Han L, Li F. Phylogeny and polymorphism in the E6 and E7 of human papillomavirus: alpha-9 (HPV16, 31, 33, 52, 58), alpha-5 (HPV51), alpha-6 (HPV53, 66), alpha-7 (HPV18, 39, 59, 68) and alpha-10 (HPV6, 44) in women from Shanghai. Infect Agent Cancer. 2019;14:38.
Pimenoff VN, de Oliveira CM, Bravo IG. Transmission between archaic and modern human ancestors during the evolution of the oncogenic human papillomavirus 16. Mol Biol Evol. 2017;34:4–19.
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–8.
Hildesheim A, Schiffman M, Bromley C, Wacholder S, Herrero R, Rodriguez A, Bratti MC, Sherman ME, Scarpidis U, Lin QQ, Terai M, Bromley RL, Buetow K, Apple RJ, Burk RD. Human papillomavirus type 16 variants and risk of cervical cancer. J Natl Cancer Inst. 2001;93:315–8.
Schiffman M, Rodriguez AC, Chen Z, Wacholder S, Herrero R, Hildesheim A, Desalle R, Befano B, Yu K, Safaeian M, Sherman ME, Morales J, Guillen D, Alfaro M, Hutchinson M, Solomon D, Castle PE, Burk RD. A population-based prospective study of carcinogenic human papillomavirus variant lineages, viral persistence, and cervical neoplasia. Cancer Res. 2010;70:3159–69.
Matsumoto K, Yasugi T, Nakagawa S, Okubo M, Hirata R, Maeda H, Yoshikawa H, Taketani Y. Human papillomavirus type 16 E6 variants and HLA class II alleles among Japanese women with cervical cancer. Int J Cancer. 2003;106:919–22.
This study has been funded and supported by Tehran University of Medical Sciences (TUMS), Grant No. (98-3-99-46124). It has also been part of a MSc thesis supported by Tehran University of Medical Sciences; Grant No. 240/984.
Ethics approval and consent to participate
Our research conducted ethically in accordance with the World Medical Association Declaration of Helsinki. We declare that informed consent was obtained from all study subjects and the study was approved by the local ethical committee of Tehran University of Medical Sciences (Grant No. 98-3-99-46124).
Consent to publication
The authors declare that they have no conflict of interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Hosseini, N., Shoja, Z., Arashkia, A. et al. Lineage analysis of human papillomavirus type 39 in cervical samples of Iranian women. Virol J 18, 152 (2021). https://doi.org/10.1186/s12985-021-01619-8
- Human papillomavirus
- Type 39
- Cervical cancer