Association of viral load and physical status of HPV-16 with survival of patients with head and neck cancer
Dayahindara Veitía1, Juan Liuzzi2, Maira Ávila1, Idamelys Rodriguez3, Felix Toro3 and María Correnti1
1Molecular Genetics Laboratory, Oncology and Haematology Institute, Ministry of Popular Power for Health (MPPS), 1053, Venezuela
2Head and Neck Service, “Padre Machado” Oncology Hospital, Ministry of Popular Power for Health (MPPS), 1053, Venezuela
3Molecular Biology Laboratory, Immunology Institute, University of Central Venezuela (UCV), 1053, Venezuela
Abstract
Introduction: Head and neck cancers (NHCs) are of multifaceted origins, and tobacco and alcohol are the primary risk factors. Currently, other factors associated with the genesis of these tumours are being considered, among these viral infections, especially human papillomavirus (HPV) infection.
Objective: The objective was to evaluate HPV infection, HPV-16 E6 load and its physical status in patients with squamous cell carcinoma in the head and neck and evaluate its effects in the survival of these patients.
Methodology: A total of 80 fresh biopsies of HNC were evaluated. The genetic material was extracted using the commercial kit QIAGEN. The detection and classification of HPV were carried out using INNO-LiPA, whereas the quantification and analysis of integration of the viral genome into the host cell were carried out using real-time PCR.
Results: The average age of the patients included was 60.34 ± 14.48 years, with a predominance of the male gender. The most frequent HPV infection was genotype 16 (52.8%), with an average of 10 copies of the HPV-16 E6/β-globin gene. Furthermore, an integration of the viral genome in the host cell was observed in 86% of cases with a statistically significant relationship between the location of the tumour and the viral load (p < 0.05).
Conclusions: HPV-16 is the most common infection, and its physical status in the host cell is the determining factor in establishing response to treatment. However, more studies are needed to demonstrate the role of HPV infection in carcinogenesis.
Keywords: head and neck, HPV, viral load, integration, survival
Correspondence to: Dayahindara Veitía
Email: dayah_20@hotmail.com
Published: 30/07/2020
Received: 27/04/2020
Publication costs for this article were supported by ecancer (UK Charity number 1176307).
Copyright: © the authors; licensee ecancermedicalscience. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Introduction
Head and neck tumours of any oncological pathology are of multifactorial origin, and tobacco and alcohol are the primary risk factors. It is estimated that between 85% and 90% of cases of head and neck squamous cell carcinomas (HNSCC) are explained by exposure to these factors; the risk is being proportional to the intensity of the exposure [1]. In recent decades, viral infections have been suggested as the risk factors associated with the development of this malignity. It has also been seen that the presence of oncogenic viruses and their interaction with the genetic material of the host provoke can trigger the alteration in the genes, which regulate the cellular cycle of the infected cells [2].
Currently, there is epidemiological and molecular evidence that human papillomavirus (HPV) is related to squamous cell carcinomas, such as in the case of cervical cancer [3] and oropharyngeal cancer [4]. This infection, probably acquired from sexual contact, could alter the epidemiology and demographics of cancer of the aerodigestive tract. In consequence, a diagnosis of squamous cell cancer positive for HPV could be clinically relevant for prognostic purposes, as well as having future diagnostic and therapeutic implications, including its prevention [5]. The numerous molecular studies suggest that the high-risk oncogenic HPV types, especially genotypes 16 and 18, are aetiologically related to a subgroup of head and neck cancers, especially oropharyngeal squamous cell carcinomas. This indicates that HPV infection in the oropharynx could lead to a predisposition to tumourogenesis [4, 6, 7, 8]. It has been possible to determine the physical status of HPV in the host cells using enlargement techniques such as real-time PCR, considering that it is integrated when the E2 gene breaks off. It is found in its episomal form when there are a number of copies equivalent to the E2 and E6 genes. It is also described as a mixed form, where there are fewer copies of the E2 gene than of E6 [9]. This is of great importance since it has been considered that the integration of the viral genome into the cellular genome is an important step to consider in the malignant transformation; the breaking off of gene E2, which controls the expression of genes E6 and E7, results in an increase in the expression of said oncogenes. E6 is related to the degradation of P53, whereas E7 with the degradation of the retinoblastoma protein (pRB); both are important in regulating the cell cycle [6].
Driven by everything set out previously, the objective of this investigation is to evaluate HPV infection, HPV-16 E6 load and its physical status in patients with squamous cell carcinoma in the head and neck and evaluate its effects in the survival of these patients.
Materials and methods
Biological material
A total of 80 fresh biopsies from HNSCC tumours were evaluated. These were gathered by Dr Juan F. Liuzzi of the Head and Neck Service in the Oncology Hospital ‘Padre Machado’ in Caracas.
Each patient included in this study had to give their approval by signing an informed consent as well as responding to a survey, which can be found as part of the FONACIT projects G-2005000408, 2013000413 and Misión Ciencia (Science Mission) LPL-20074001088.
Extraction of the genetic material
The extraction of the nucleic acids was carried out using the commercial DNA Mini Kit QIAGen (250) following the firm’s specifications.
From each biopsy, a piece of tissue approximately 25 mg was cut, and this, in turn, was cut into small pieces and incubated in ATL buffer (180 μL) and 20 μL of proteinase K overnight with agitation at 56ºC. Afterwards, 200 μL of buffer AL was added, vortex mixed for 15 sec and incubated at 70ºC for 10 min.
200 μl of absolute ethanol was added, and then, it was transferred to a tube and centrifuged at 6,000 xg for 1 min; Nearly 500 μL of the buffer of washed AW1 was added and centrifuged at 6,000 xg for 1 min; 500 μL of the buffer of AW2 was added and centrifuged at 20,000 xg for 3 min. Then, 200 μL of AE buffer was added, and it was incubated at room temperature for 1 min and centrifuged at 6,000 xg for 1 min. The genetic material gathered was stored at −80ºC until its subsequent use.
Detection and genotyping of the human papillomavirus
HPV detection was performed using the INNO-LiPA HPV Genotyping Extra Kit (Innogenetics), following the commercial specifications. The kit allows specific detection of 25 HPV genotypes (HPV types 6, 11, 16, 18, 31, 33, 35, 39, 40, 42, 43, 44, 45, 51, 52, 53, 54, 56, 58, 59, 66, 68, 70, 73, 74), and it is based on reverse hybridisation. The test included a PCR amplification of a 65 bp fragment within the L1 region of the HPV genome using the broad-spectrum SPF10 biotinylated primers. Biotinylated amplicons are subsequently hybridised with HPV type-specific oligonucleotide probes which are immobilised as parallel lines on membrane strips. After hybridisation and stringent washing, streptavidin-conjugated alkaline phosphatase is added and bound to any biotinylated hybrid formed. Incubation with BCIP/NBT chromogen yields a purple precipitate, and the results can be visually interpreted.
Quantification and evaluation of the physical status of the HPV-16 genome
The quantification of HPV-16 E6 as well as evaluation of the genome’s integration was carried out using real-time PCR according to the specifications described by Koskinen et al. [9].
Real-time PCR was performed with the 7500 Real-Time PCR System (Applied Biosystem) in the Molecular Biology Laboratory at the Immunology Institute of the Central University of Venezuela. The sequences for the initiators are shown in Table 1. The amplification and quantification of the E6 and E2 genes were carried out simultaneously in separate reaction tubes. The reaction mixture consisted of 50 μL made up of 1X SSoFast Evagreen (BIORAD), 900 nM of each of the initiators and approximately 50 ng of the DNA. The amplification process consisted of 1 cycle at 95°C for 10 min, 44 denaturing cycles at 95°C for 15 sec and the binding and extension phase at 60°C for 1 min, followed by the dissociation curve.
All the samples that tested positive for HPV genotype 16 were evaluated to check whether gene E2 had broken off, as well as quantifying genes E2 and E6 of HPV genotype 16.
The calibration curve was performed using DNA from HeLa cells, the initiators GP5 /GP6 were used, which allow for the detection of the HPV genome, and the initiators PC04 and GH20 were used to amplify the β-globin human gene (shown in Table 1). The extracted DNA was quantified (Beckman Instruments DV series Spectrophotometer. 260–280 nm). Afterwards, the serial dilutions of this material were made which were used for the respective qRT-PCR. For HPV, 1.5 μm of each initiator was used, 1X SsoFast Eva Green in a final volume of 20 μL. For the amplification and quantification of the β-globin gene, 1 μM of the initiators PC04 y GH20, 2X SsoFast Evagreen were used in a final volume of 50 μl. The amplification process used was 95°C followed by 40 cycles at 95°C, hybridisation at 46°C for 1 min and extension at 72°C over 30 sec. Once the amplification processes had been performed both for the standard curve and the viral load of the E6 and E2 genes, electrophoresis was carried out in a 2% agarose gel dyed with ethidium bromide (0.2 μL of a solution of 1%) and exposed to UV light on an imaging camera ChemiDoc™ XRS (Bio-Rad) as a photographic record, corroborating that the process had been carried out successfully.
Table 1. Sequences of the initiators used in the real-time PCR to quantify the viral genome.
Statistical analysis
The relationship between the viral infection and the development of tumours was evaluated using the Chi-square test (X2), and the statistical significance was considered for values of p < 0.05.
For the analysis and estimation of the patient survival curves, the Kaplan–Meier method was used. The statistical significance was considered for values of p < 0.05.
In both the cases, the statistical program IBM SPSS Statistics 25 was used.
Results
The average age of the patients was 60.34 ± 14.48 years (range: 19–89). The most frequent gender was males with 68.75% (p = 0.001). Regarding the associated risk factors, 73.5% (p = 0.000) of the patients were smokers and 70% consumed alcohol (p = 0.000). The most frequent location of the tumour was the oral cavity (33.8%), followed by the larynx (30%) and the oropharynx (22.5%). In addition, the paranasal sinuses (5%), nasopharynx (3.8%), hypopharynx (2.5%) and nasal cavities (2.5%) are also presented as the locations of the primary tumour (Table 2).
Regarding the histopathological analysis of the tumours, 40% corresponded to extensive tumours (T4) and moderately differentiated (43.75%), in the advanced states of the illness (EIII and EIV).
Sixty-six percent of the samples evaluated were positive for HPV infection (p = 0.004), and HPV-16 was the most frequent form of single infection (17%). It could be noted that the infection for this viral genotype was present in the form of a coinfection or mixed infection in a total of 52.8% of those testing positive for HPV-16 (Table 3).
Regarding the presence of the infection according to the anatomical location of the primary tumour, it was observed that the oral cavity, larynx and oropharynx were the most frequent with 34%, 28% and 23%, respectively; the majority of the HPV-positive tumours were found in advanced states EIII (26%) and EIV (55%) (Table 3). Of these patients, 60% were smokers and 58% consumed alcohol (data not shown).
The average viral load was 10 copies of DNA of HPV for each copy of β-globin (range between 2 and 20 copies of DNA of HPV), and the oropharynx was the anatomical location with the largest number of copies. At this site, with an elevated copies number of E6-HPV-16, the viral genome was found in an episomal or mixed form. The oral cavity and the larynx, on the contrary, presented low copies numbers of E6-HPV-16, and in these, the viral DNA was found in an integrated form (p = 0.00) (Table 4). The average frequency of integration was 86%; a statistically significant relationship was observed between the low viral load and the anatomical location with integration (p = 0.016).
After following the patients for 120 months, it was observed that there was a tendency towards the reduction in the survival of the male patients, there was some engagement in other risk behaviours (alcohol and/or tobacco) and they were diagnosed in advanced stages. In addition, it could be noted that the physical status of the viral genome in the host cell played an important role; in those cases, when the viral DNA presented in an episomal form, they had a more constant survival than those patients whose viral DNA presented in an integrated form although without a statistically significant relationship (p > 0.05) (Figures 1 and 2).
Discussion
HNSCC presents primarily in males, and the age of appearance is being an important factor [10, 11]. In this investigation, the average age was maintained within the range described of 60.34 ± 14.48 years (range 19–89), and in agreement with the previous reports, this type of malignity presents most frequently between the fifth and seventh decades of life, with a mean average of 62 ± 12 years [12].
Table 2. Characteristics of the samples in the study.
Sixty-six percent of the samples evaluated were positive for the presence of the HPV genome, this value is higher than that reported by the other works which report the viral detection averages of 39%–40% [13, 14], 60% according to Yang et al. (2019) [15] and even 61% according to Koskinen et al. [9]; however, it remains within the averages reported by the other authors [16, 17]. Regarding genotype, HPV-16 was the most frequent in single infection form (17%), noting that the infection in this genotype be it in the form of co-infection or mixed infection was present in a considerable number of the analysed cases. This is within the ranges described by Koskinen et al [9] and Tao et al [18], who reported 27% and 85% positivity for HPV, respectively. Another important aspect which could be noted was that, in initial stages of the illness, the frequency of HPV infection is low (4%), and this increases as the illness stage advances (55% in advanced stages). This supports the claims carried out in other anatomical areas such as the cervix, indicating that a persistent infection with high-risk oncogenic genotypes, amongst them HPV-16, is an important risk factor for the genesis and progression of a tumour [19]. In addition, 60% of patients with HPV engage in other risk factors such as alcohol and tobacco, suggesting that these could also act as the risk factors for HPV infection and that combined (alcohol, tobacco, HPV and other infections) they could all act in a synergistic way in the process of development and progression of the malignity.
The oral cavity, larynx and oropharynx were, in this order, the anatomical locations most frequently infected with HPV. This is contrary to publications which describe the oropharynx as the area of the head and neck with the greatest infection rate [13]. However, the results of various investigations agree that these three areas are those most frequently infected by HPV. On the other hand, it suggests that the oral cavity is ‘the door of entry of the virus’, and this could coincide with an increased frequency of viral infection. It also suggests that HPV infections are a biological factor, which plays an important role in the carcinogenesis of these anatomical locations, based on epidemiological and clinicopathological evidence [20]. HPV-16 is being identified between 29% and 90% of these evidence [21–24]. It also proposed that the presence of this viral genotype in exfoliated cells of the oral cavity increases the risk of oropharyngeal cancer by 14 times [25].
On the contrary, the result of the average viral load investigation agrees with the previous works which explain that it can vary in these types of tumours and can be anywhere from 3.4 × 10-4 up to values as high as 739 copies/reference gene [9, 26]. If the results are compared with those of other anatomical locations, such as the oesophagus, for example, it could be said that the number of copies of HPV-16 E6 is similar, approximately 2.55 ± 3.19 copies per cell [27]. In this vein, Si et al [28] reported that, in China, the average viral load varies between ˂1 and 157 copies per cell in this anatomical location; as it can be seen, these are significantly low values if they are compared with those described for the cervix, vagina and vulva. whose reports reach values of more than 1 million copies of the HPV-16 E6 gene [27, 29].
Table 3. HPV infection in patients with HNSCC.
Table 4. HPV load and physical status in patients with HNSCC.
Figure 1. Survival curves of patients with squamous cell carcinoma in the head and neck. A) Patients with HNSCC in relation to tobacco consumption (n = 80; p = 0.187). B) Patients with HNSCC in relation to tumour stage (n = 80; p = 0.087). C) Patients with HNSCC in relation to gender (n = 80. p = 0.171).
As can be seen, viral load is not constant and presents notable differences between anatomical locations. It has been of great use evaluating the integration of the viral genome to the host cell genome since the previous studies reveal that the state of the viral genome is an important prognostic factor that, together with epigenetic factors, they are related with the cancer progression [30]. In this sense, it has been reported that those cases with episomal HPV show a favourable prognosis, suggesting that integration could, in some way, modify the biology of the tumours [31]. Eighty-six percent of the samples included in this study showed the integration of the viral genome with the host cell genome. This is consistent with other investigations, such as carried out by Olthof et al [32]; Otlthoft et al [33] and Mellin et al [31], who all indicate a high percentage of integration. Equally, they show that viral load is variable in HNSCC and that a high viral load is associated with the preferably episomal status, which is related with a significantly better prognosis, whereas low loads are related with integration and a poor prognosis. This was observed, in this study, in the cases of oral cavity and larynx tumours when the viral load was low, where the largest percentage of integration was detected. On the contrary, in the oropharyngeal tumours, the viral load was high in relation to the other anatomical locations evaluated, and the viral DNA was found in the preferable episomal form, which was related with a better prognosis (p = 0.016).
Figure 2. Analysis of survival of patients with squamous cell carcinoma in the head and neck in relation to: A) HPV genome presence (without distinguishing genotype) (n = 53; p > 0.05). B) Integration of the viral genome into the host cell genome (n = 28. p > 0.05).
The investigations attribute these differences to the fact that the patients with high viral load possibly induce a better immune response against cancer although this hypothesis has not been confirmed and requires more study to support it. It is important to mention that although there are notable variations in the viral load between different anatomical locations, the average percentage of integration is similar to that described in the investigation carried out in the cervix by Groves and Coleman [34]. It is also important to mention that this process of integration requires the breaking of both the viral genome and the host cell’s genome, and it is believed that the rate of integration is proportional to the level of damage suffered by the DNA. These damages have been described as endogenous and exogenous, including the effects of alcohol and tobacco as well as HPV infection or other infectious agents, such as the Epstein–Barr virus (EBV), for example. All are the risk factors for the development of HNSCC, and this triggers an activation in the DNA repair mechanisms, as well as the process of chromosome alteration both of them could contribute to the process of integration [35].
After following the survival of the patients related to each one of the variables, for a period of 120 months, it could be noted that there was a downward trend in the survival of the male patients, something that was also noted in those patients who engage in risk activities such as consumption of alcohol or tobacco as well as those who presented with a diagnosis in an advanced stage (EIII and EIV). Equally, it was not possible to demonstrate that the presence of HPV infection would have a direct effect on survival. However, it would appear that an integration of the viral genome could be the determining factor since in those cases, where the HPV-16 genome was integrated, the survival was lower than in those cases, where it was in an episomal form. Although there was a statistically significant evidence to support this, this could be attributed to the low number of samples included (Figures 1 and 2).
Conclusions
The oral cavity is the most common anatomical location for tumour growth, which could be due to it being in direct contact with aetiological agents. On the other hand, the presence of HPV-16 is established as the most common viral infection, considering the physical status of the viral genome as the determining factor for being able to establish the response to the employed therapies. It is suggested that HPV-positive patients, especially with oropharyngeal cancer, have a better prognosis than those who develop this type of cancer in other anatomical locations in the upper aerodigestive tract; these differences in prognosis could be due to the physical status of the viral genome in the host cell. However, the need to carry out more studies should be highlighted, to be able to evidence the role of HPV infection in carcinogenesis.
Acknowledgments
The research was financed by FONACIT under the research projects G-2005000408, 2013000413 and Misión Ciencia (Science Mission) LPL-20074001088.
Conflicts of interest
The authors have no conflicts of interest to declare.
References
1. Spitz M, Fueger J, and Goepfert H (1988) Squamous cell carcinoma of the aerodigestive tract. A case comparison analysis Cancer 61 203 https://doi.org/10.1002/1097-0142(19880101)61:1<203::AID-CNCR2820610134>3.0.CO;2-6 PMID: 3334949
2. Rosenquist K, Wennerberg J, and Schildt E, et al (2005) Oral status, oral infections and some lifestyle factors as risk factors for oral and oropharyngeal squamous cell carcinoma. A population-based case-control study in southern Sweden Act Otolaryngol 125 1327–1336 https://doi.org/10.1080/00016480510012273
3. Martinez I, Wang J, and Hobson K, et al (2007) Identification of differentially expressed genes in HPV-positive and HPV-negative oropharyngeal squamous cell carcinomas Eur J Cancer 43(2) 415–432 https://doi.org/10.1016/j.ejca.2006.09.001 PMCID: 1847595
4. Badulescu F, Crisan A, and Badulescu A, et al (2010) Recent data about the role of human papillomavirus (HPV) in oncogenesis of head and neck cancer Romanian J Morphol Embryol 51(3) 437–440
5. Haigentz M (2005) Aerodigestive cancers in HIV infection Curr Opin Oncol 17 474–478 https://doi.org/10.1097/01.cco.0000174036.46785.8f PMID: 16093799
6. Gillison ML (2006) Human papillomavirus and prognosis of oropharyngeal squamous cell carcinoma: implication for clinical research in head and neck cancer J Clin Oncol 24 5623–5625 https://doi.org/10.1200/JCO.2006.07.1829 PMID: 17179099
7. Mork J, Lie A, and Glattre E, et al (2001) Human papillomavirus types in head and neck squamous-cell carcinoma of the head and neck N Engl J Med 344(15) 1125–31
8. Burns J and Maitland, N (2005) Human papillomaviruses and cancer Microbiol Tod 56 116–120
9. Koskinen W, Chen R, and Leivo I, et al (2003) Prevalence and physical status of human papillomavirus in squamous cell carcinomas of the head and neck Int J Cancer 107 401–406 https://doi.org/10.1002/ijc.11381 PMID: 14506740
10. Arias F, Villafranca E, and Dueñas MT, et al (2007) Factores de pronóstico en carcinoma de cabeza y cuello (Pamplona: Servicio de oncología, Hospital de Navarra)
11. Veitía D, Liuzzi J, and Ávila M, et al (2014) Human papilloma virus detection in head and neck squamous cell carcinoma ecancermedicalscience 8 475–482
12. Salazar C, Smith R, and Garg M, et al (2014) Human papillomavirus-associated head and neck squamous cell carcinoma survival: a comparison by tumor site and initial treatment Head Neck Pathol 8 77–87 https://doi.org/10.1007/s12105-013-0486-4 PMCID: 3950385
13. Faust H, Eldenhed E, and Roslin A, et al (2016) Prevalence of human papillomavirus types, viral load and physical status of HPV16 in head and neck squamous cell carcinoma from the South Swedish Health Care Region J Gen Virol 97 2949–2956 https://doi.org/10.1099/jgv.0.000611 PMID: 27667722
14. Harlé A, Guillet J, and Thomas J, et al (2019) HPV insertional pattern as a personalized tumor marker for the optimized tumor diagnosis and follow-up of patients with HPV-associated carcinomas: a case report BMC Cancer 19 277–284 https://doi.org/10.1186/s12885-019-5447-1 PMID: 30922253 PMCID: 6437879
15. Yang R, Klimentová J, and Göckel-Krzikalla E, et al (2019) Combined transcriptome and proteome analysis of immortalized human keratinocytes expressing human papillomavirus 16 (HPV16) oncogenes reveals novel key factors and networks in HPV-induced carcinogenesis mSphere 4 e00129–e0019. https://doi.org/10.1128/mSphere.00129-19 PMID: 30918060 PMCID: 6437273
16. Veitía D, Liuzzi J, and Ávila M, et al (2017) Prevalence of HPV and EBV infection and their relationship with the p53 and PCNA expression in oral carcinoma patients J Oral Res 6(4) 86–91 https://doi.org/10.17126/joralres.2017.026
17. Cordero A, Ruíz A, and Veitía D, et al (2017) Detección molecular de la infección por VPH en pacientes con cáncer de cabeza y cuello Rev venezolana de oncología 29(3) 194–200
18. Tao Y, Sturgis E, and Huang Z, et al (2018) TGF-β1 genetic variants predict clinical outcomes of HPV positive oropharyngeal cancer patients after definitive radiotherapy Clin Cancer Res 24(9) 2225–2233 https://doi.org/10.1158/1078-0432.CCR-17-1904 PMID: 29463556 PMCID: 5932112
19. Del Río L, Soto S, and Camargo M, et al The DNA load of six high-risk human papillomavirus types and its association with cervical lesions BMC Cancer 15 100-111
20. Betiol J, Villa L, and Sichero L (2013) Impact of HPV infection on the development of head and Neck cancer Brazilian J Med Biol Res 46 217–226 https://doi.org/10.1590/1414-431X20132703
21. Badaracco G, Venuti A, and Morello R, et al (2000) Human papillomavirus in head and neck carcinomas: prevalence, physical status and relationship with clinical/pathological parameters Anticancer Res 20 1301–1306 PMID: 10810439
22. Ringstrom E, Peters E, and Hasegawa (2002) Human papillomavirus type 16 and squamous cell carcinoma of the head and neck Clin Cancer Res 8 3187–3192 PMID: 12374687
23. Correnti M, Rivera H, and Caazza ME (2004) Detection of human papillomavirus of high oncogenic potential in oral squamous cell carcinoma in a Venezuelan population Oral Oncol Microbiol Oral Dis 10 163–166
24. Castillo A, Koriyama C, and Higashi M, et al (2011) Human papillomavirus in upper digestive tract tumors from three countries J Clin Virol 17 5295–5304
25. Kreimer A, Clifford G, and Boyle P, et al (2005) Human papillomavirus types in head and neck squamous cell carcinomas worldwide: a systematic review Cancer Epidemiol Biomarkeres Prev 14 467–475 https://doi.org/10.1158/1055-9965.EPI-04-0551
26. Speel EJ (2017) HPV integration in head and neck squamous cell carcinomas: cause and consequence Recent Results Cancer Res 206 57–72 https://doi.org/10.1007/978-3-319-43580-0_4
27. Zhang D, Zhang Q, and Zhou L, et al (2010) Comparison of prevalence, viral load, physical status and expression of human papillomavirus- 16, -18 and -58 in esophageal and cervical cancer: a case-control study BMC Cancer 10 650–658 https://doi.org/10.1186/1471-2407-10-650
28. Si H, Wah S, and Ping C, et al (2005) Physical status of HPV-16 in esophageal squamous cell carcinoma J Clin Virol 32 19–23 https://doi.org/10.1016/j.jcv.2004.04.004
29. Lillsunde G, Helenius G, and Sorbe B, et al (2014) Viral load, integration and methylation of e2bs3 and 4 in human papilloma virus (HPV) 16-positive vaginal and vulvar carcinomas PLoS One 9(11) e112839 https://doi.org/10.1371/journal.pone.0112839
30. Hatano T, Sano D, and Takahashi H, et al (2017) Identification of human papillomavirus (HPV) 16 DNA integration and the ensuing patterns of methylation in HPV-associated head and neck squamous cell carcinoma cell lines Int J Cancer 140(7) 1571–1580 https://doi.org/10.1002/ijc.30589
31. Mellin H, Dahlgren L, and Muner E, et al (2002) Human papollomavirus type 16 is an episomal and high viral load may be correlated to better prognosis in tonsillar cancer Int J Cancer 1002 152–158 https://doi.org/10.1002/ijc.10669
32. Olthof N, Straetmans J, and Snoeck R, et al (2012) Next generation treatment strategies for HPV-related head and neck squamous cell carcinomas: where do we go? Rev Med Virol 22 88–105 https://doi.org/10.1002/rmv.714
33. Olthof NC, Huebbers CU, and Kolligs J, et al (2015) Viral load, gene expression andmapping of viral integration sites in HPV16-associated HNSCC cell lines Int J Cancer 136 E207–E218 https://doi.org/10.1002/ijc.29112
34. Groves I and Coleman N (2015) Pathogenesis of human papillomavirus-associated mucosal disease J Pathol 235 527–538 https://doi.org/10.1002/path.4496 PMID: 25604863
35. Visalli G, Riso R, and Facciolà A, et al (2016) Higher levels of oxidative DNA damage in cervical cells are correlated with the grade of dysplasia and HPV infection J Med Virol 88 336–344 https://doi.org/10.1002/jmv.24327
