1 Volume 21 2022 e226337 Original Article Braz J Oral Sci. 2022;21:e226337http://dx.doi.org/10.20396/bjos.v21i00.8666337 1 Department of Oral Medicine, School of Dentistry, Tehran University of Medical Sciences, Tehran, Iran. 2 Assistant Professor, Razi Vaccine and Serum Research Institute, Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran. 3 Associate Professor, Department of Oral Medicine, Imam Khomeini Hospital, Tehran University of Medical Science, Tehran, Iran. Corresponding author: Soheila Manifar, Associate Professor, Department of oral medicine, Imam Khomeini Hospital, Tehran University of Medical Science, Tehran, Iran. Tel: +98-2188351150, +98-912-033-0580 E-mail: mariakoopaie@gmail.com Address: North Kargar St,Tehran, Iran Po. Code: 14399-55991 P.O.BOX:14395-433. Editor: Dr. Altair A. Del Bel Cury Received: July 10, 2021 Accepted: April 2, 2022 Lineage and phylogenetic analysis of HPV-16, -18 in saliva of HNSCC patients Maryam Koopaie1 , Mohamad Amin Nematollahi1 , Maryam Dadar2 , Soheila Manifar3,* Aim: Head and Neck Squamous Cell Carcinoma (HNSCC) is a global health problem whose incidence varies by geographic region and race according to risk factors. Human papillomavirus (HPV) infection is a significant risk factor for HNSCC. HPV-16 and HPV-18 are two forms of HPV that are carcinogenic. HNSCCs that are HPV positive have a better prognosis rather than HPV negative. The purpose of this research was to characterize HPV-16, -18 variations in the saliva of HNSCC patients by examining the genetic diversity of HPV-16, -18 utilizing the full E6, E7, and L1 genes. Methods: The case-control research included 15 patients with HNSCC and 15 healthy volunteers. Unstimulated entire saliva samples were obtained from the case and control groups by spitting method. Genomic DNA was isolated from all saliva samples. A PCR reaction was used to determine the presence of HPV in saliva. HPV-positive samples were genotyped and data were analyzed. We conducted a variant study on the HPV-16, -18 E6, and E7 genes. Results: Three patients with HNSCC were HPV- positive for two HPV genotypes out of 30 people diagnosed with HPV-DNA. HPV-16 and -18 were the most common genotypes. The HPV-16, -18 E6, and E7 genes were sequenced and compared to the HPV-16, -18 (E6, E7) prototype sequence. In all, HPV-16 lineages A1 and HPV-18 lineages A3 were discovered. Conclusion: Regarding the variation of HPV found in Iranian HNSCC patients, the need for further studies in HPV genotyping was seen. Sequencing HPV genes in HNSCC may help answer questions about HPV genotyping in the Iranian population. HPV genotype analysis aids in the development of vaccinations against HNSCC, halting disease progression and preventing HPV-associated HNSCC. Keywords: Squamous cell carcinoma of head and neck. Saliva. Alphapapillomavirus. Phylogeny. https://orcid.org/0000-0002-9999-1443 https://orcid.org/0000-0003-0605-781X https://orcid.org/0000-0001-5831-801X https://orcid.org/0000-0002-5898-3893 2 Koopaie et al. Braz J Oral Sci. 2022;21:e226337 Introduction Squamous cell carcinoma (SCC) accounts for 90% of oral cavity cancers, and about 400000 new instances of head and neck squamous cell carcinoma (HNSCC) are diagnosed each year1. Men are 1.5 times more likely than women to develop HNSCC2. The prevalence of HNSCC varies by geographic region and race due to risk factor differences3. The incidence of HNSCC is higher in India, Pakistan, and Bangladesh in the southern part of Asia1,4. It may appear on the mucosa as a red, white, or red-white prominence. It may also be portrayed as a severe wound with noticeable edges. It is often linked to lymph node swelling5. HNSCC is a complex ill- ness with a variety of intrinsic and extrinsic risk factors that contribute to its devel- opment. Tobacco and alcohol usage, infections, and dietary variables are examples of extrinsic risk factors. Iron deficiency anemia and alterations in cell cycle control genes are among the inherent risk factors for HNSCC6,7. In contrast to other viruses, human papillomavirus (HPV) has a significant role in HNSCC pathogenesis among extrinsic risk factors8. HPV is a DNA virus that attacks the skin and mucous mem- branes. There are about 30 different types of HPV viruses that have been linked to benign and malignant epithelial lesions9. Furthermore, high-risk HPV strains such as HPV-16 and HPV-18 are thought to be etiologic factors for HNSCC10,11. The pres- ence of HPV-16 in exfoliated oral epithelial cells may raise the risk of HNSCC and oropharyngeal cancer by 4 and 14 times, respectively12. According to reports, HPV- 16 is the most common kind in nearly all geographical areas of Iran, followed by HPV-1813,14. HPV-16, HPV-18, HPV-11, and HPV-52 as the most frequently identi- fied HPV types in Iranian patients with HNSCC12,15. Although the data on HPV types in Iran is well-documented, there is little information on HPV variations in HNSCC. The current research sought to examine basic data of HPV-16 and HPV-18 variants in Iranian patients with HNSCC by examining the genetic variability of HPV-16 and HPV-18 utilizing the whole E6, E7, and L genes. The findings of our study might pro- vide justification for future research into HPV genetic epidemiology, pathogenicity, and evolution. Material and method Ethical statement This study was approved by the Tehran University of Medical Sciences Ethical Committee (Ethical code IR.TUMS.DENTISTRY.REC.1398.036). After describing the study objectives, written informed consent was obtained from all case and control group subjects. Study population This case-control study was performed in Imam Khomeini Hospital, Tehran Uni- versity of Medical Sciences (Tehran, Iran). 15 HNSCC patients and 15 healthy vol- unteers as the control group enrolled in the study. HNSCC of all patients was con- firmed with the histopathologic examination, and they had not undergone surgery, chemotherapy, and radiotherapy. The study participants were not affected by any 3 Koopaie et al. Braz J Oral Sci. 2022;21:e226337 other malignancy or systemic disease. In the intraoral examination, there was no active dental and periodontal infection. The control group was matched with the case group for age, sex, tobacco, and alcohol consumption. Unstimulated whole saliva samples were collected from the case and control group using the spitting method with a minimum of 5 ml saliva. All participants in the case and control group were asked not to eat, drink or brush their teeth, smoke at least one hour prior to the trial, and not consume any alcoholic beverages 24 hours prior to collection1. The authors confirm that all methods were performed in accordance with the rele- vant guidelines and regulations. Evaluation of the presence of HPV with MY09 / 11 primer Genomic DNA was extracted from saliva samples of the case (HNSCC patients) and control group using the High Pure Viral Nucleic Acid Kit (Roche Diagnostics GmbH, Roche Applied Science, Mannheim, Germany) as directed by the manufacturer. The DNA concentration was measured using the ND-1000 Nanodrop at 260/280 nm. Furthermore, the DNA integrity was tested using a 1.5 percent agarose gel and kept at 20 °C for future study. The presence of HPV in the saliva was determined using the MY09 / MY11 oligonucleotide primers 5’-CGT CCM ARR GGA WAC TGA TC-3’ and 5’-CGT CAG GGW CAT AAY AAT GG-3’, which amplify a region of about 450 pb in the genome of most HVP types16. A 25 L reaction mixture containing 10 pmol of each primer, 50 M of each dNTP, 200 ng of DNA template, 2 U of Taq DNA poly- merase, 20 mM Tris-HCl, and 3 mM MgCl2 was used for the PCR process. The PCR was started with a 5 minute denaturation at 95 °C, followed by 35 cycles of dena- turation for 30 seconds at 95 °C, annealing for 30 seconds at 55 °C, and extension for 30 seconds at 72 °C, with a final extension at 72 °C for 7 minutes. A negative PCR control, consisting of a reaction mixture devoid of template DNA, was added to each pair of PCR reactions. Electrophoresis on a 1 percent agarose gel was used to examine the PCR results. Variant analysis of HPV-16 E6 and E7 genes PCR was used to examine the HPV-16 E6 and E7 genes, with an amplicon size of 762 bp targeted (Table 1). A 50 L reaction mixture containing 10 pmol of each primer, 50 M of each dNTP, 200 ng of DNA template, 2 U of Taq DNA polymerase, 20 mM Tris-HCl, and 3 mM MgCl2 was used for the PCR process. An initial dena- turation of 5 minutes at 95 °C was followed by 40 cycles of denaturation at 95 °C for 30 seconds, annealing at 55 °C for 50 seconds, extension at 72 °C for 50 sec- onds, and final elongation at 72 °C for 10 minutes. A negative control, consist- ing of a reaction mixture devoid of template DNA, was added to each pair of PCR reactions. Electrophoresis on a 1 percent agarose gel was used to examine the PCR results. Variant analysis of HPV-16 L1 gene HPV-16 L1 gene was analyzed by PCR (Table 1), targeting an amplicon size of 7932 bp. 4 Koopaie et al. Braz J Oral Sci. 2022;21:e226337 Table 1. Primers for PCR analyzing of HPV-16 E6, E7 and L1 and HPV-18 E6, E7 and L1 Primer sequences Genes 5’-CGA AAC CGG TTA GTA TAA AAG CAG AC-3’ HPV-16 E6 5’-TAG ATT ATG GTT TCT GAG AAC A-3’ HPV-16 E7 5’-ATG TGC CTG TAT ACA CGG GTC-3’ HPV-16 L1, Forward 5’-TTA CTT CCT GGC ACG TAC ACG C HPV-16 L1, Reverse 5’-ATG GCG CGC TTT GAG GAT CC -3’ HPV-18 E6 5’- TTA CTG CTG GGA TGC ACA CC -3’ HPV-18 E7 5’- ATG TGC CTG TAT ACA CGG GTC -3’ HPV-18 L1, Forward 5’- TTA CTT CCT GGC ACG TAC ACG C -3’ HPV-18 L1, Reverse PCR reaction method is similar to that described for HPV-16 E6 and E7 genes, and the schematic of HPV-16 E6 primers design is depicted in Figure 1. Figure 1. Schematic primers design for PCR analyzing of HPV-16 E6. Variant analysis of HPV-18 E6 and E7 genes The HPV-18 E6 and E7 genes were amplified using polymerase chain reaction (PCR) with an amplicon size of 835 base pairs as the target (Table 2). A 50-L reaction mix- ture containing 10 pmol of each primer, 50 M of each dNTP, 200 ng of DNA template, 2 U of Taq DNA polymerase, 20 mM Tris-HCl, and 3 mM MgCl2 was used for the PCR process. The primers were desaturated for 5 minutes at 95 °C, followed by 40 cycles of denaturation at 95 °C for 30 seconds, annealing at 55 °C for 50 seconds, extension at 72 °C for 50 seconds, and final elongation at 72 °C for 10 minutes. Each set of PCR reactions contained a negative control consisting of a reaction mixture devoid of tem- plate DNA. Electrophoresis on a 1% agarose gel was used to verify the PCR results. 5 Koopaie et al. Braz J Oral Sci. 2022;21:e226337 Variant analysis of HPV-18 L1 gene The HPV-18 L1 gene was examined using a PCR assay with an amplicon size of 7857 bp (Table 1). The PCR reaction technique is identical to that described previ- ously for the HPV-16 E6 and E7 genes, and a schematic primer design is shown in Figure 2 for reference. Figure 2. Schematic primers design for PCR analyzing of HPV-18 L1, R. ABI Prism 377 sequencer Genetic Analyzer equipment (Applied Biosystems Fos- ter City, Canada) was used to sequence all of the PCR products from the E6 and E7 genes to evaluate the HPV-16 and HPV-18 variants. The sequences were run according to the Big-Dye Terminator protocol (Applied Biosystems) (GeneAlI, Seoul, South Korea). It was determined that the reference sequences were unique to each lineage and sublineage by aligning them to all previously known lineages and sub- lineages of A1-4, B1-4, C2-4, and D1-316. The maximum likelihood approach was employed to build the phylogenetic tree, which was done using Mega software version 6. In order to analyze the HPV-18 lineages and sublineages, all HPV-18 sequences were analyzed using the Mega software version, which was used to depict the phylogenetic tree according to the maximum likelihood method with the reliability of bootstrap on 1000 replicates using the maximum likelihood method with bootstrap reliability on 1000 replicates. The A1-5, C, and B1-3 reference sequences for HPV-18 were also obtained from the GenBank database, as were the lineage and sublineage-specific reference sequences for HPV-18. 6 Koopaie et al. Braz J Oral Sci. 2022;21:e226337 Results This research comprised 30 people diagnosed with HPV-DNA over a year (from March 2020 to March 2021). Ten patients were male and five were female in the patient group. The control group consisted of 9 males and 6 females. The mean age and standard deviation (SD) of patient and control group were 53.84±13.61 and 45.23±12.80, respectively. Four patients in the HNSCC group had tongue SCC, six had oropharyngeal SCC, three had soft palate SCC, and two had laryngeal SCC. Figure 3. HPV-16 and HPV-18 (E6, E7) were sequenced in genomes isolated from saliva of subjects. When the sequences recovered in this research were compared to the sequences of samples in the HPV-18 E6, E7, and L1 gene banks, it was discovered that the viral strain isolated from the patient’s saliva belonged to lineage A3 (Figure 4-D, E, F). Over- all, lineages A1 of HPV-16 (Figure A, B, C) and A3 of HPV-18 variations (Figure 4-D, E, F) were discovered (Figure 4), as shown by the red outlined data. 7 Koopaie et al. Braz J Oral Sci. 2022;21:e226337 A B C D E F Fi gu re 4 . T he p hy lo ge ne tic tr ee in c om pa ris on w ith th e se qu en ce o f s am pl es in th e H P V -1 6 E6 , E 7 an d L1 a nd H P V -1 8 E6 , E 7 an d L1 g en e ba nk s ho w ed th at th e vi ru s st ra in is ol at ed fr om th e pa tie nt ’s s al iv a be lo ng s to ty pe 1 6; li ne ag es A 1 (A , B a nd C ) a nd ty pe 1 8; li ne ag es A 3 (D , E a nd F ). 8 Koopaie et al. Braz J Oral Sci. 2022;21:e226337 Discussion HNSCC is a group of heterogenic cancers with multiple risk factors. HPV positiv- ity may be an important contributing factor in the histologic grade and prognosis of cases with HNSCC2. Risk factors of HPV transmission include multiple sex partners, tobacco consumption, and HIV infection. The roots of HPV transmission include sex- ual or nonsexual contact with the affected person, transmission through saliva con- taminated with the viruses, breastfeeding, transmission during labor3. In the elderly, the most important predisposing factors for HNSCC are alcohol and tobacco use, but in younger people, the role of alcohol, and tobacco in HNSCC is lower and the number of young people with HPV are more than older ones. Hence, the HPV virus may be one of the major risk factors for HNSCC in the young population4. The annual incidence of oropharyngeal SCC was increased, while SCC incidence in the other sites of oral cav- ity was decreased. These differences may be attributed to the role of HPV in HNSCC pathogenesis5. The role of HPV virus in anogenital cancers is confirmed. Furthermore, HPV has a role in 20 – 40 % of HNSCC6. PCR test to detect the expression of E6 and E7 onco- genes in HPV virus is the best diagnostic test for detection of HPV in HNSCC7. The carcinogenic’s mechanism of HNSCC is attributed to E6 and E7 oncogenes. HPV-E7 protein disrupts the pRb-mediated cell-cycle control and causes a reduction in cyclin D1 synthesis and the over-expression of the cyclin-dependent kinase 4/6 inhibitor p16INK4A. HPV E6 protein promotes the destruction of P53 tumor suppressor gene8. HPV positive HNSCC patients have a better prognosis than HPV negative patients, and their response to treatment is better9. Some studies suggested that HPV-18 is the dominant type of HPV in Iranian HNSCC patients, but these reports are controver- sial and this controversy attribute to the geographic, cultural, and habit difference10,11. Wood et al.13 showed that the prevalence of oral HPV varied significantly in different geographical areas. Seifi et al.11 introduced HPV-18 as the dominant type of virus in East Azerbaijan province of Iran. Delavarian et al.14 found no correlation between oral SCC and HPV infection. These results were in line with Sibers et al.12 and Sisk et al.15 but in contrast with Kreimer et al.17. A systematic review study by Jalilvand et al. stated that HPV was detected in 44.4% of patients with head and neck cancers16. The most common HPV types were HPV-16 and -18. There was a strong association between HPV-16 and oropharyngeal cancer. These studies reflect the controversial reports about HPV prevalence and its subtypes in Iran. This controversy recently con- firmed by a systematic review study18. In conclusion, the HPV vaccine is available for the prevention of cervical cancer. Gar- dasil is a type of HPV vaccine which induces immunity against HPV-6, -11, -16, -18. Clinical trials on the effectiveness of this vaccine showed over 98% efficacy in the prevention of cervical, anal, vulvar, and vaginal cancers19. Genomic analysis of HPV in the saliva of HNSCC patients in the present study is a starting point for designing a vaccine to prevent HPV-associated HNSCC. Disclosure Statement No potential conflict of interest was reported by the authors. 9 Koopaie et al. Braz J Oral Sci. 2022;21:e226337 Acknowledgments The authors would like to acknowledge the help and support by Imam Khomeini Hos- pital, Tehran University of Medical Sciences and School of Dentistry, Tehran Univer- sity of Medical Sciences for their constant assistance with the research. Data availability Datasets related to this article will be available upon request to the corresponding author. Authors Contribution MK and MD conceived the study idea and led data collection. MK, MAN, and MD created the study protocol and wrote the original draft. MAN and SM contributed in collecting the saliva sample. MD, MK and MAN contributed to data analysis / interpretation and preparation of the manuscript. MK, MAN and SM led the writing- review & editing. MK and MD interpreted the results. All authors read and approved the final manuscript. References 1. Jasim H, Carlsson A, Hedenberg-Magnusson B, Ghafouri B, Ernberg M. Saliva as a medium to detect and measure biomarkers related to pain. Sci Rep. 2018 Feb 19;8(1):3220. doi: 10.1038/s41598-018-21131-4. 2. Wang S, Zhuang X, Gao C, Qiao T. Expression of p16, p53, and TLR9 in HPV-Associated Head and Neck Squamous Cell Carcinoma: Clinicopathological Correlations and Potential Prognostic Significance. Onco Targets Ther. 2021 Feb 5;14:867-77. doi: 10.2147/OTT.S293163. 3. Wierzbicka M, Klussmann JP, San Giorgi MR, Wuerdemann N, Dikkers FG. Oral and laryngeal HPV infection: Incidence, prevalence and risk factors, with special regard to concurrent infection in head, neck and genitals. Vaccine. 2021 Apr 22;39(17):2344-2350. doi: 10.1016/j.vaccine.2021.03.047. 4. McDermott JD, Bowles DW. Epidemiology of head and neck squamous cell carcinomas: impact on staging and prevention strategies. Curr Treat Options Oncol. 2019;20(5):1-13. doi: 10.1007/s11864-019-0650-5. 5. Boscolo-Rizzo P, Zorzi M, Del Mistro A, Da Mosto MC, Tirelli G, Buzzoni C, et al. The evolution of the epidemiological landscape of head and neck cancer in Italy: Is there evidence for an increase in the incidence of potentially HPV-related carcinomas? PLoS One. 2018;13(2):e0192621. doi: 10.1371/journal.pone.0192621. 6. Kofler B, Laimer J, Bruckmoser E, Steinbichler TB, Runge A, Schartinger VH, et al. The role of HPV and non-HPV sexually transmitted infections in patients with oropharyngeal carcinoma: A case control study. Cancers (Basel). 2020;12(5):1192. doi: 10.3390/cancers12051192. 7. Augustin JG, Lepine C, Morini A, Brunet A, Veyer D, Brochard C, et al. HPV detection in head and neck squamous cell carcinomas: what is the issue? Front Oncol. 2020 Sep 15;10:1751. doi: 10.3389/fonc.2020.01751. 8. Rodriguez J, Herrero A, Li S, Rauch N, Quintanilla A, Wynne K, et al. PHD3 regulates p53 protein stability by hydroxylating proline 359. Cell Rep. 2018 Jul 31;24(5):1316-29. doi: 10.1016/j.celrep.2018.06.108. 9. Kobayashi K, Hisamatsu K, Suzui N, Hara A, Tomita H.Miyazaki T. A review of HPV-related head and neck cancer. Journal of clinical medicine. 2018;7(9):241. doi:10.3390/jcm7090241. 10 Koopaie et al. Braz J Oral Sci. 2022;21:e226337 10. Kermani IA, Seifi S, Dolatkhah R, Sakhinia E, Dastgiri S, Ebrahimi A, et al. Human papilloma virus in head and neck squamous cell cancer. Iranian journal of cancer prevention. 2012;5(1):21. doi: 10.1200/jco.2011.29.15_suppl.e16011. 11. Seifi S, Kermani IA, Dolatkhah R, Kermani AA, Sakhinia E, Asgarzadeh M, et al. Prevalence of oral human papilloma virus in healthy individuals in East azerbaijan province of iran. Iran J Public Health. 2013;42(1):79-85. 12. Siebers T, Merkx M, Slootweg P, Melchers W, Van Cleef P, de Wilde P. No high-risk HPV detected in SCC of the oral tongue in the absolute absence of tobacco and alcohol—a case study of seven patients. Oral Maxillofac Surg. 2008 Dec;12(4):185-8. doi: 10.1007/s10006-008-0131-7. 13. Wood ZC, Bain CJ, Smith DD, Whiteman DC.Antonsson A. Oral human papillomavirus infection incidence and clearance: a systematic review of the literature. J Gen Virol. 2017 Apr;98(4):519-26. doi: 10.1099/jgv.0.000727. 14. Delavarian Z, Pakfetrat A, Falaki F, Pazouki M.Pazouki N. The role of viruses in oral squamous cell carcinoma in young patients in Khorasan (Northeast of Iran). J Appl Sci. 2010;10(11):981-5. doi: 10.3923/jas.2010.981.985. 15. Sisk EA, Soltys SG, Zhu S, Fisher SG, Carey TE.Bradford CR. Human papillomavirus and p53 mutational status as prognostic factors in head and neck carcinoma. Head & Neck: Journal for the Sciences and Specialties of the Head and Neck. 2002;24(9):841-49. 16. Jalilvand S, Shoja Z.Hamkar R. Human papillomavirus burden in different cancers in Iran: a systematic assessment. Asian Pac J Cancer Prev. 2014;15(17):7029-35. 17. Kreimer AR, Johansson M, Waterboer T, Kaaks R, Chang-Claude J, Drogen D, et al. Evaluation of human papillomavirus antibodies and risk of subsequent head and neck cancer. J Clin Oncol. 2013;31(21):2708. 18. Aghaeipour F, Salehiniya H.Abbaszadeh H. Prevalence of human papillomavirus (HPV) in oral mucosal lesions in Iran: A systematic review and meta-analysis. J Med Virol. 2021;93(11):6089-99. 19. Toh ZQ, Kosasih J, Russell FM, Garland SM, Mulholland EK.Licciardi PV. Recombinant human papillomavirus nonavalent vaccine in the prevention of cancers caused by human papillomavirus. Infection and drug resistance. 2019;12:1951. doi:10.2147/IDR.S178381.