The Integration Pattern of Human Papillomavirus in the Host Genome as a Promising Guide for Precise Radiation Therapy in Cervical Cancer

authors:

avatar Mansoureh Gorginzadeh 1 , avatar Mahsa Motavaf 2 , *

Endometriosis Research Center, Iran University of Medical Sciences, Tehran, Iran
Functional Neurosurgery Research Center, Shohada Tajrish Neurosurgical Comprehensive Center of Excellence, Shahid Beheshti University of Medical Sciences, Tehran, Iran

how to cite: Gorginzadeh M , Motavaf M. The Integration Pattern of Human Papillomavirus in the Host Genome as a Promising Guide for Precise Radiation Therapy in Cervical Cancer. Artif Intell Precis Med. 2021;1(1):e115501. https://doi.org/10.5812/pmco.115501.

Abstract

Context:

Cervical cancer (CC) is one of the most common cancers among women worldwide. Infection with high-risk human papillomavirus (HPV) is the leading risk factor for this cancer. Currently, concurrent chemotherapy and radiotherapy are the standard of care for patients with locally advanced CC. However, the response to radiotherapy varies widely among patients, underscoring the need to identify biomarkers to predict patients’ treatment outcomes.

Evidence Acquisition:

In order to find relevant articles, we performed a systematic search using ISI Web of Science, Scopus, PubMed, and Embase, with a combination of the keywords “human papillomaviruses”, “cervical cancer”, “radiotherapy”, “biomarkers”, “viral integration”, and “episome”. All types of studies, including observational or retrospective studies, clinical trials, congress abstracts, and case reports or case series, were examined to find relevant data. No time limitation was considered for the studies.

Results:

The integration of HPV DNA into the host genome is considered to be a critical step in HPV carcinogenesis. The HPV integration pattern has been investigated as a potent prognostic factor for the tailored treatment of CC.

Conclusions:

It is indicated that the integration state of the HPV genome in host DNA is a potential prognostic or predictive biomarker for CC response to radiotherapy.

1. Context

Papillomavirus is a large family of viruses belonging to the family Papovaviridae, the genus Papillomavirus. To date, more than 200 types of human papillomavirus (HPV) have been recognized based on DNA sequence (1). About 40 types of HPV are sexually transmitted and affect fully differentiated squamous epithelium of the skin and mucosal epithelium of anogenital and non-genital sites, among which 15 types can cause cervical cancer (CC) (2). The immune system spontaneously clears most HPV infections (approximately 90%) within 1 or 2 years. ‘High-risk’ mucosal HPV types, predominantly types 16, 18, 31, 33, and 35, have been associated with most cervical, penile, vulvar, vaginal, anal, and oropharyngeal cancers and pre-cancers (3). In general, HPV is thought to be responsible for more than 90% of anal and cervical cancers, about 70% of vaginal and vulvar cancers, and more than 60% of penile cancers (4). In addition, HPV16 and HPV18 are accountable for approximately 70% of CC worldwide (5).

Structurally, HPV16 and HPV18 are small non-enveloped viruses with icosahedral capsid enclosing a circular double-stranded DNA. The total length of viral DNA is about 8000 base pairs (8 kb), encompassing eight open reading frames transcribed as polycistronic messages from a single DNA strand. Based on functional properties and location, the genome can be divided into three regions: an early region (E), a late region (L), and a long control region (LCR). The early region contains six ORFs (E1, E2, E4, E5, E6, and E7), which carry out regulatory functions. The late region encodes L1 and L2, the principal virus capsid protein (Figure 1). The L1 region, with 90% similarity, is the most conserved genome region in HPV and is the prime target for serological assay and vaccine development (6). The extended control region is a noncoding region that holds the origin of replication and several transcription factor binding sites. E1 and E2 proteins are responsible for virus amplification (replication and DNA packaging) and maintenance. Persistent HPV infections are reported to be associated with integrating the HPV genome into the host cell genome. The process of integration frequently occurs within E1 or E2 and leads to their loss. This leads to overexpression of viral oncoproteins E6 and E7 and genetic alterations known as drivers of mutation in HPV carcinogenesis.

The typical genomic structure of high-risk HPV types
The typical genomic structure of high-risk HPV types

The development of CC often takes many years or even decades after persistent infection with HPV. This long time window provides a golden opportunity to elucidate HPV‐host interaction, malignant transformation, and carcinogenesis. HPV genotype, viral load, physical state, integration site, and genomic change at the insertion site are suggested as key features with predictive and prognostic values for clinical management of HPV-related cancers (7).

Nearly all locally advanced HPV-related cancers, regardless of their clinical stage, tumor histology, or tumor size, are treated with a similar paradigm of surgery, external or internal radiation therapy, and systemic or regional chemotherapy (8). However, the standard treatment regimen does not reflect the biologic heterogeneity of HPV-related cancers. Therefore, to improve radiotherapy outcomes, it is necessary to optimize the intensity of radiotherapy based on the individualpatient’s responses. Given the advent of next-generation sequencing (NGS) technology and the implementation of omics approaches (genomics, transcriptomics, and proteomics) in clinical oncology, cancer treatments shift from one-size-fits-all toward more comprehensive, multidisciplinary, precision approaches. However, clinically useful biomarkers to guide radiotherapy optimization are not yet available.

In the present brief review, we aimed to summarize recent research on the prognostic value of the physical status of HPV DNA for the survival rate of patients with HPV-related cancers undergoing radiotherapy.

2. Evidence Acquisition

We performed a systematic search using ISI Web of Science, Scopus, PubMed, and Embase, with a combination of the keywords “human papillomaviruses”, “cervical cancer, “radiotherapy”, “biomarkers”, “viral integration”, and “episome”. All types of studies, including observational or retrospective studies, clinical trials, congress abstracts, and case reports or case series, were investigated to find relevant data with no time limitation.

3. Results

3.1. HPV Physical State

Following infection, the HPV virus can remain in its episomal form or become integrated into the host genome. Besides, both episomal and integrated patterns can be present simultaneously (9). Two main mechanisms of HPV integration into the genomic DNA include looping (which may result in viral-host DNA concatemers) integration and direct insertion. In addition, HPV can be integrated into the host genome as the single-copy or multi-copy of head-to-tail tandem repeats. Single copy HPV integration chiefly occurs somewhere in the E1 and E2 gene region, which normally regulates the expression of the viral oncogenes E6 and E7. This integration results in disruption of the E1 or E2 open reading frame [ORS], leading to deregulation of negative feedback control of the oncogenes E6 and E7. The E6 oncoprotein forms a physical complex with the p53 cellular tumor suppressor, resulting in P53-mediated G1/S checkpoint dysfunction and cell cycle dysregulation (10).

On the other hand, the E7 oncoprotein binds to the retinoblastoma tumor suppressor protein, pRB, releasing the transcription factor E2F from pRb (11). This results in transcriptional activation of E2F-responsive genes participating in cell-cycle control. These events contribute to HPV oncogenesis and progression of low-grade lesions to neoplastic conditions and, ultimately, carcinoma.

In contrast to one copy integration of the HPV genome, when numerous copies are being integrated, only the E1 or E2 of copies surrounding the cellular DNA will be disrupted, while the remaining copies, E1 and E2 ORFs, will stay intact; thus, they can maintain transcriptional control on E6/E7 (12).

Even though HPV oncogenesis is usually regarded as dependent on viral genomic integration, the rate of HPV DNA integration is quite variable according to HPV genotypes. For instance, the integration rate is reported to be almost 100% for HPV-18 positive CC cases but about 80% for HPV-16 cancer cases (13). The importance of integration in the carcinogenesis of CC is further supported by the fact that HPV DNA has been integrated into merely 3% of premalignant lesions (14-16).

Nevertheless, HPV integration is not the only requisite factor in the carcinogenesis of CC as purely episomal forms of the viral genome have also been detected in HPV-related cancers (15, 17-19). According to preliminary studies, the pathogenesis of HPV episome-related cancer, which is associated with numerous copies of episomal HPV, mainly results from overexpression of the viral oncogenes E6 and E7 (19, 20). However, subsequent studies have suggested some roles for genetic or epigenetic changes resulting in dysregulation of E6/E7 gene expression. For instance, the E2-mediated repression of viral oncogenesis decreased through methylation of E2 binding sites in the viral upstream regulatory region (21, 22). This refers to the impact of HPV episome on HPV-related carcinogenesis apart from E2 disruption secondary to integrating the viral genome into the host genome.

In order to better understand the integration mechanism in the carcinogenic process, the physical state (integrated or episomal) of the viral genome should be determined. Of many different techniques being applied, including Southern blot, fluorescent in situ hybridization, and PCR-based approaches like restriction site PCR and the detection of integrated papillomavirus sequences PCR, none was able to differentiate between the two forms of the HPV genome (23, 24). Discrimination between HPV DNA patterns was also not possible when real-time PCR that can measure both E2 and E6 genes was used (25).

In a study by Joo et al. (26), both in situ hybridization [ISH] and PCR were adopted to identify the HPV integration pattern. Locally advanced cervical tumors were categorized into three groups: the episomal group, the single-copy integrated or multi-copy tandem repetition-integrated pattern group, and the undetectable HPV group. According to their results, what happens first is the disruption of the E2 region without being amplified by PCR after the full integration of the virus. Then, following HPV insertion in the form of multi-copy tandem repetition, given that unbroken E2 genes have been repeatedly put into host cells, both E2 and E6 genes can be traced. Due to the multiplicity of E2 and E6 copies in both episomal and mixed patterns, it was concluded that the sole use of PCR was not sufficient to differentiate between different HPV DNA patterns. They evaluated 204 cervical tumors using both PCR and ISH, although there was significant discordance between these methods. Therefore, since ISH can better differentiate between episomal and multi-copy tandem integrated patterns, it can be combined with PCR to improve discrimination (26).

Lately, integration events, like specifying viral integration sites into the human genome and identifying disrupted genes (esp. oncogenes, tumor suppressor genes, and DNA repair genes) can be successfully detected using the NGS technique (27).

Both targeted (the whole exome sequencing of DNA or polyA-specific RNA sequencing) and untargeted, unbiased (the whole-genome sequencing of DNA or total RNA sequencing) methods have been used in available studies. Some other novel targeted approaches, including capture sequencing technologies, have also been introduced, which are based on DNA enrichment before sequencing and allow for deeper analysis of specific regions, including genome-wide promoter assessment (28).

So far, only a few recurrent sites for HPV integration known as genomic “hotspots”, which are strongly associated with regions in host DNA active in transcription and susceptible to damage, have been identified by whole-genome association studies (23). When DNA replication 2stress ensues from HPV integration, structures called stalled replication forks are developed at regions of the host genome, which sometimes show microscopic homology of 1 to 1 base pairs to the viral genome and have high AT bases. HPV integration can also induce regulatory or mutagenic effects on adjacent genes.

Considering the importance of HPV integration into the human genome as an early key stage in malignancy, recognizing the physical state of the HPV genome in low-grade or atypical lesions can have important clinical implications in both identifying prognostic factors for disease progression and selecting more precise therapeutic approaches (29, 30).

3.2. HPV Physical State as a Prognostic Factor for Radiation Response in CC

Cervical malignancy is caused by high-risk HPV types in nearly all cases. When diagnosed at an early stage, CC can be cured using a surgical approach. As an alternative to surgery or complementary treatment, radiotherapy can also be used for some women with early-stage disease. For more advanced cases, however, a combination of both radiotherapy and chemotherapy is often needed. Given that HPV-related tumors are generally well-responsive to radiotherapy, the implementation of a precise radiotherapeutic approach improves both quality and outcome of treatment. However, even in individuals with similar clinical and morphological parameters of their malignant lesion, including anatomic region, histological type, clinical stage, and degree of cancerous cell differentiation, radiosensitivity of tumor cells is quite variable. Therefore, by identifying some unresponsive to conventional radiotherapy, the radiation therapy schedule can be optimized by elucidating the underlying molecular mechanisms. It is believed that radiosensitivity or radioresistance of tumor cells is related to the biological behavior of their genes. Hence, decoding genetic changes within these cells can aid in improving clinical outcomes (31). As the integration of HPV DNA in host somatic cells is a critical step in carcinogenesis, determining the physical status of the HPV genome and the site of integrating or targeting host genes is presumed to modify the disease-free survival rate (32, 33).

In an early study by Unger et al. (34), the ISH approach was applied to determine HPV-18 DNA (episomal, integrated, or mixed) in primary CC tumors and evaluate the correlation of overall disease-free survival rates with the HPV DNA integration rate. According to their results, the integration of HPV DNA in the primary tumor was strongly associated with treatment failure (34). In another complementary study, the same authors evaluated the association between the integration of disrupted E2 and the overall disease-free survival in patients with HPV 16-positive primary CC tumors (32). They thoroughly examined the whole E1/E2 region by amplifying eight separate fragments spanning the region using ISH. Based on their results, the amount of disrupted E2 in sample cells was significantly correlated with reduced disease-free survival compared to other patients with intact E2 (32). In another study, Kalantari et al. (35) applied the PCR method to map E1 and E2 reading frames in HPV+ CC samples. They demonstrated that the 5’ end disruption was the most common pattern among different disruption/deletion patterns in the E1-E2 area of the HPV 16 genome in cervical carcinoma patients. The survival rate was poorer in patients diagnosed with carcinomas showing disruptions in E1/E2 than those without such changes. Concerning prognosis, E1 disruptions were considered critical (35). Contrary to the aforementioned studies in which combinations of therapeutic approaches were adopted, Lindel et al.’s study (33) was the first to include patients with HPV 16-positive uterine cervix cancer who were supposed to only receive radiotherapy as a curative modality. The results of their study are consistent with other studies demonstrating that disease-free survival is better in patients with intact E2 than in patients with disrupted E2.

In addition to the E2 gene status, another significant prognostic factor regarding the post-radiotherapy survival rate is the proportion of episomal to integrated copies of HPV DNA (29). According to a study by Shin et al. (29), cervical tumors with more episomal components are associated with better post-radiotherapy survival rates compared to purely integrated HPV DNA cells. In this study, ISH and the PCR technique were applied to differentiate between HPV DNA integration patterns. A multivariate model was used to categorize patients according to their different risk profiles. The best prognosis was found in tumors mainly composed of episomal HPV DNA, whereas the worst prognosis was detected in tumors with undetectable episomal DNA.

In a recent study by Zamulaeva et al. (36), radioresistance of CC cells was found to be dependent on both integration of HPV DNA into the host cell genome and amount of this integration. This study demonstrated that the percentage of cancer stem cells (CSCs) increased after irradiation in tumors with fully integrated HPV DNA whereas the proportion of CSCs decreased in tumors with episomal or partly integrated forms. This increase in the CSC pool after irradiation results from the subsequent proliferation of radio-resistant cancer cells and the differentiation of non-stem cells to stem cells as it occurs in radiation-induced epithelial-mesenchymal transition. Viral oncoproteins whose expression increases following HPV DNA integration into the host cell genome are known to enhance this transition (37).

What is fascinating is that the HPV DNA integration site, apart from the integration rate, appears to have a prognostic significance. Although HPV DNA integration is thought to be a random event, several studies have pointed out some hotspots that may affect cancer evolution (38-40).

In a study by Yang et al. (41), single-cell sequencing was employed to analyze mutations before and after radiotherapy in cervical tumors. They concluded that genetic characteristics were different in tumor cells pre- and post-radiotherapy. They concluded that genetic characteristics were different in tumor cells pre- and post-radiotherapy. Based on their results, tumor cells with HPV integration in POU5F1B, a well-known cancer susceptibility locus (42, 43), and a hotspot for HPV integration (44) are more likely to survive the following radiotherapy. This is assumed to be one of the possible underlying genetic mechanisms related to radiotherapy resistance.

4. Conclusions

Currently, radiotherapy alone or in combination with chemotherapy is the primary treatment for CC. However, the clinical outcome of treated patients is not the same as an endpoint survival rate after radiotherapy. Thus, it is necessary to identify potential biomarkers useful in personalizing therapy and optimizing radiotherapy outcomes. Several available reports have suggested that HPV integration patterns are an important and potentially useful clinical biomarker to guide high precision radiation therapy.

References

  • 1.

    Winer RL, Gheit T, Feng Q, Stern JE, Lin J, Cherne S, et al. Prevalence and Correlates of beta- and gamma-Human Papillomavirus Detection in Oral Samples From Mid-Adult Women. J Infect Dis. 2019;219(7):1067-75. [PubMed ID: 30395247]. [PubMed Central ID: PMC6420166]. https://doi.org/10.1093/infdis/jiy632.

  • 2.

    Adebamowo SN, Adeyemo AA, Rotimi CN, Olaniyan O, Offiong R, Adebamowo CA, et al. Genome-wide association study of prevalent and persistent cervical high-risk human papillomavirus (HPV) infection. BMC Med Genet. 2020;21(1):231. [PubMed ID: 33225922]. [PubMed Central ID: PMC7682060]. https://doi.org/10.1186/s12881-020-01156-1.

  • 3.

    Smith JS, Lindsay L, Hoots B, Keys J, Franceschi S, Winer R, et al. Human papillomavirus type distribution in invasive cervical cancer and high-grade cervical lesions: a meta-analysis update. Int J Cancer. 2007;121(3):621-32. [PubMed ID: 17405118]. https://doi.org/10.1002/ijc.22527.

  • 4.

    Chaturvedi AK. Beyond cervical cancer: burden of other HPV-related cancers among men and women. J Adolesc Health. 2010;46(4 Suppl):S20-6. [PubMed ID: 20307840]. https://doi.org/10.1016/j.jadohealth.2010.01.016.

  • 5.

    Forman D, de Martel C, Lacey CJ, Soerjomataram I, Lortet-Tieulent J, Bruni L, et al. Global burden of human papillomavirus and related diseases. Vaccine. 2012;30 Suppl 5:F12-23. [PubMed ID: 23199955]. https://doi.org/10.1016/j.vaccine.2012.07.055.

  • 6.

    Woodman CB, Collins SI, Young LS. The natural history of cervical HPV infection: unresolved issues. Nat Rev Cancer. 2007;7(1):11-22. [PubMed ID: 17186016]. https://doi.org/10.1038/nrc2050.

  • 7.

    Harle A, Guillet J, Thomas J, Demange J, Dolivet G, Peiffert D, et al. 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. 2019;19(1):277. [PubMed ID: 30922253]. [PubMed Central ID: PMC6437879]. https://doi.org/10.1186/s12885-019-5447-1.

  • 8.

    Miles B, Safran HP, Monk BJ. Therapeutic options for treatment of human papillomavirus-associated cancers - novel immunologic vaccines: ADXS11-001. Gynecol Oncol Res Pract. 2017;4:10. [PubMed ID: 28725449]. [PubMed Central ID: PMC5512733]. https://doi.org/10.1186/s40661-017-0047-8.

  • 9.

    Oyervides-Munoz MA, Perez-Maya AA, Rodriguez-Gutierrez HF, Gomez-Macias GS, Fajardo-Ramirez OR, Trevino V, et al. Understanding the HPV integration and its progression to cervical cancer. Infect Genet Evol. 2018;61:134-44. [PubMed ID: 29518579]. https://doi.org/10.1016/j.meegid.2018.03.003.

  • 10.

    Thomas M, Pim D, Banks L. The role of the E6-p53 interaction in the molecular pathogenesis of HPV. Oncogene. 1999;18(53):7690-700. [PubMed ID: 10618709]. https://doi.org/10.1038/sj.onc.1202953.

  • 11.

    Berezutskaya E, Yu B, Morozov A, Raychaudhuri P, Bagchi S. Differential regulation of the pocket domains of the retinoblastoma family proteins by the HPV16 E7 oncoprotein. Cell Growth Differ. 1997;8(12):1277-86. [PubMed ID: 9419416].

  • 12.

    Wagatsuma M, Hashimoto K, Matsukura T. Analysis of integrated human papillomavirus type 16 DNA in cervical cancers: amplification of viral sequences together with cellular flanking sequences. J Virol. 1990;64(2):813-21. [PubMed ID: 2153245]. [PubMed Central ID: PMC249176]. https://doi.org/10.1128/JVI.64.2.813-821.1990.

  • 13.

    Cullen AP, Reid R, Campion M, Lorincz AT. Analysis of the physical state of different human papillomavirus DNAs in intraepithelial and invasive cervical neoplasm. J Virol. 1991;65(2):606-12. [PubMed ID: 1846186]. [PubMed Central ID: PMC239798]. https://doi.org/10.1128/JVI.65.2.606-612.1991.

  • 14.

    Li H, Yang Y, Zhang R, Cai Y, Yang X, Wang Z, et al. Preferential sites for the integration and disruption of human papillomavirus 16 in cervical lesions. J Clin Virol. 2013;56(4):342-7. [PubMed ID: 23290390]. https://doi.org/10.1016/j.jcv.2012.12.014.

  • 15.

    Dutta S, Chakraborty C, Dutta AK, Mandal RK, Roychoudhury S, Basu P, et al. Physical and methylation status of human papillomavirus 16 in asymptomatic cervical infections changes with malignant transformation. J Clin Pathol. 2015;68(3):206-11. [PubMed ID: 25563334]. https://doi.org/10.1136/jclinpath-2014-202611.

  • 16.

    Mazumder Indra D, Singh RK, Mitra S, Dutta S, Chakraborty C, Basu PS, et al. Genetic and epigenetic changes of HPV16 in cervical cancer differentially regulate E6/E7 expression and associate with disease progression. Gynecol Oncol. 2011;123(3):597-604. [PubMed ID: 21911249]. https://doi.org/10.1016/j.ygyno.2011.08.004.

  • 17.

    Cheung JL, Cheung TH, Yu MY, Chan PK. Virological characteristics of cervical cancers carrying pure episomal form of HPV16 genome. Gynecol Oncol. 2013;131(2):374-9. [PubMed ID: 24012799]. https://doi.org/10.1016/j.ygyno.2013.08.026.

  • 18.

    Li W, Wang W, Si M, Han L, Gao Q, Luo A, et al. The physical state of HPV16 infection and its clinical significance in cancer precursor lesion and cervical carcinoma. J Cancer Res Clin Oncol. 2008;134(12):1355-61. [PubMed ID: 18478264]. https://doi.org/10.1007/s00432-008-0413-3.

  • 19.

    Hong D, Liu J, Hu Y, Lu X, Li B, Li Y, et al. Viral E6 is overexpressed via high viral load in invasive cervical cancer with episomal HPV16. BMC Cancer. 2017;17(1):136. [PubMed ID: 28202002]. [PubMed Central ID: PMC5311840]. https://doi.org/10.1186/s12885-017-3124-9.

  • 20.

    Karbalaie Niya MH, Keyvani H, Safarnezhad Tameshkel F, Salehi-Vaziri M, Teaghinezhad SS, Bokharaei Salim F, et al. Human Papillomavirus Type 16 Integration Analysis by Real-time PCR Assay in Associated Cancers. Transl Oncol. 2018;11(3):593-8. [PubMed ID: 29547758]. [PubMed Central ID: PMC5854915]. https://doi.org/10.1016/j.tranon.2018.02.017.

  • 21.

    Reuschenbach M, Huebbers CU, Prigge ES, Bermejo JL, Kalteis MS, Preuss SF, et al. Methylation status of HPV16 E2-binding sites classifies subtypes of HPV-associated oropharyngeal cancers. Cancer. 2015;121(12):1966-76. [PubMed ID: 25731880]. https://doi.org/10.1002/cncr.29315.

  • 22.

    Bhattacharjee B, Sengupta S. CpG methylation of HPV 16 LCR at E2 binding site proximal to P97 is associated with cervical cancer in presence of intact E2. Virology. 2006;354(2):280-5. [PubMed ID: 16905170]. https://doi.org/10.1016/j.virol.2006.06.018.

  • 23.

    Thorland EC, Myers SL, Persing DH, Sarkar G, McGovern RM, Gostout BS, et al. Human papillomavirus type 16 integrations in cervical tumors frequently occur in common fragile sites. Cancer Res. 2000;60(21):5916-21. [PubMed ID: 11085503].

  • 24.

    Luft F, Klaes R, Nees M, Durst M, Heilmann V, Melsheimer P, et al. Detection of integrated papillomavirus sequences by ligation-mediated PCR (DIPS-PCR) and molecular characterization in cervical cancer cells. Int J Cancer. 2001;92(1):9-17. [PubMed ID: 11279600].

  • 25.

    Arias-Pulido H, Peyton CL, Joste NE, Vargas H, Wheeler CM. Human papillomavirus type 16 integration in cervical carcinoma in situ and in invasive cervical cancer. J Clin Microbiol. 2006;44(5):1755-62. [PubMed ID: 16672403]. [PubMed Central ID: PMC1479176]. https://doi.org/10.1128/JCM.44.5.1755-1762.2006.

  • 26.

    Joo J, Shin HJ, Park B, Park SY, Yoo CW, Yoon KA, et al. Integration Pattern of Human Papillomavirus Is a Strong Prognostic Factor for Disease-Free Survival After Radiation Therapy in Cervical Cancer Patients. Int J Radiat Oncol Biol Phys. 2017;98(3):654-61. [PubMed ID: 28581408]. https://doi.org/10.1016/j.ijrobp.2017.02.226.

  • 27.

    Tuna M, Amos CI. Next generation sequencing and its applications in HPV-associated cancers. Oncotarget. 2017;8(5):8877-89. [PubMed ID: 27784002]. [PubMed Central ID: PMC5352450]. https://doi.org/10.18632/oncotarget.12830.

  • 28.

    Bewicke-Copley F, Arjun Kumar E, Palladino G, Korfi K, Wang J. Applications and analysis of targeted genomic sequencing in cancer studies. Comput Struct Biotechnol J. 2019;17:1348-59. [PubMed ID: 31762958]. [PubMed Central ID: PMC6861594]. https://doi.org/10.1016/j.csbj.2019.10.004.

  • 29.

    Shin HJ, Joo J, Yoon JH, Yoo CW, Kim JY. Physical status of human papillomavirus integration in cervical cancer is associated with treatment outcome of the patients treated with radiotherapy. PLoS One. 2014;9(1). e78995. [PubMed ID: 24427262]. [PubMed Central ID: PMC3888442]. https://doi.org/10.1371/journal.pone.0078995.

  • 30.

    Vozenin MC, Lord HK, Hartl D, Deutsch E. Unravelling the biology of human papillomavirus (HPV) related tumours to enhance their radiosensitivity. Cancer Treat Rev. 2010;36(8):629-36. [PubMed ID: 20417037]. https://doi.org/10.1016/j.ctrv.2010.03.010.

  • 31.

    Wong YF, Selvanayagam ZE, Wei N, Porter J, Vittal R, Hu R, et al. Expression genomics of cervical cancer: molecular classification and prediction of radiotherapy response by DNA microarray. Clin Cancer Res. 2003;9(15):5486-92. [PubMed ID: 14654527].

  • 32.

    Vernon SD, Unger ER, Miller DL, Lee DR, Reeves WC. Association of human papillomavirus type 16 integration in the E2 gene with poor disease-free survival from cervical cancer. Int J Cancer. 1997;74(1):50-6. [PubMed ID: 9036869]. https://doi.org/10.1002/(sici)1097-0215(19970220)74:1<50::aid-ijc9>3.0.co;2-#.

  • 33.

    Lindel K, de Villiers EM, Burri P, Studer U, Altermatt HJ, Greiner RH, et al. Impact of viral E2-gene status on outcome after radiotherapy for patients with human papillomavirus 16-positive cancer of the uterine cervix. Int J Radiat Oncol Biol Phys. 2006;65(3):760-5. [PubMed ID: 16682151]. https://doi.org/10.1016/j.ijrobp.2006.01.023.

  • 34.

    Unger ER, Vernon SD, Thoms WW, Nisenbaum R, Spann CO, Horowitz IR, et al. Human papillomavirus and disease-free survival in FIGO stage Ib cervical cancer. J Infect Dis. 1995;172(5):1184-90. [PubMed ID: 7594652]. https://doi.org/10.1093/infdis/172.5.1184.

  • 35.

    Kalantari M, Karlsen F, Kristensen G, Holm R, Hagmar B, Johansson B. Disruption of the E1 and E2 reading frames of HPV 16 in cervical carcinoma is associated with poor prognosis. Int J Gynecol Pathol. 1998;17(2):146-53. [PubMed ID: 9553812]. https://doi.org/10.1097/00004347-199804000-00009.

  • 36.

    Zamulaeva I, Selivanova E, Matchuk O, Kiseleva V, Mkrtchyan L, Krikunova L. Radiation Response of Cervical Cancer Stem Cells Is Associated with Pretreatment Proportion of These Cells and Physical Status of HPV DNA. Int J Mol Sci. 2021;22(3). [PubMed ID: 33535561]. [PubMed Central ID: PMC7867083]. https://doi.org/10.3390/ijms22031445.

  • 37.

    Zamulaeva IA, Selivanova EI, Kiseleva VI, Matchuk ON, Krikunova LI, Mkrtchyan LS, et al. Correlation of Radiation Response of Cervical Cancer Stem Cells with Their Initial Number before Treatment and Molecular Genetic Features of Papillomavirus Infection. Bull Exp Biol Med. 2020;170(2):241-5. [PubMed ID: 33263849]. https://doi.org/10.1007/s10517-020-05043-w.

  • 38.

    Zhang R, Shen C, Zhao L, Wang J, McCrae M, Chen X, et al. Dysregulation of host cellular genes targeted by human papillomavirus (HPV) integration contributes to HPV-related cervical carcinogenesis. Int J Cancer. 2016;138(5):1163-74. [PubMed ID: 26417997]. [PubMed Central ID: PMC5057319]. https://doi.org/10.1002/ijc.29872.

  • 39.

    Zhang Y, Koneva LA, Virani S, Arthur AE, Virani A, Hall PB, et al. Subtypes of HPV-Positive Head and Neck Cancers Are Associated with HPV Characteristics, Copy Number Alterations, PIK3CA Mutation, and Pathway Signatures. Clin Cancer Res. 2016;22(18):4735-45. [PubMed ID: 27091409]. [PubMed Central ID: PMC5026546]. https://doi.org/10.1158/1078-0432.CCR-16-0323.

  • 40.

    Koneva LA, Zhang Y, Virani S, Hall PB, McHugh JB, Chepeha DB, et al. HPV Integration in HNSCC Correlates with Survival Outcomes, Immune Response Signatures, and Candidate Drivers. Mol Cancer Res. 2018;16(1):90-102. [PubMed ID: 28928286]. [PubMed Central ID: PMC5752568]. https://doi.org/10.1158/1541-7786.MCR-17-0153.

  • 41.

    Yang D, Zhang W, Liu Y, Liang J, Zhang T, Bai Y, et al. Single-cell whole-genome sequencing identifies human papillomavirus integration in cervical tumour cells prior to and following radiotherapy. Oncol Lett. 2018;15(6):9633-40. [PubMed ID: 29928338]. [PubMed Central ID: PMC6004640]. https://doi.org/10.3892/ol.2018.8567.

  • 42.

    Pomerantz MM, Ahmadiyeh N, Jia L, Herman P, Verzi MP, Doddapaneni H, et al. The 8q24 cancer risk variant rs6983267 shows long-range interaction with MYC in colorectal cancer. Nat Genet. 2009;41(8):882-4. [PubMed ID: 19561607]. [PubMed Central ID: PMC2763485]. https://doi.org/10.1038/ng.403.

  • 43.

    Yeager M, Orr N, Hayes RB, Jacobs KB, Kraft P, Wacholder S, et al. Genome-wide association study of prostate cancer identifies a second risk locus at 8q24. Nat Genet. 2007;39(5):645-9. [PubMed ID: 17401363]. https://doi.org/10.1038/ng2022.

  • 44.

    Hu Z, Zhu D, Wang W, Li W, Jia W, Zeng X, et al. Genome-wide profiling of HPV integration in cervical cancer identifies clustered genomic hot spots and a potential microhomology-mediated integration mechanism. Nat Genet. 2015;47(2):158-63. [PubMed ID: 25581428]. https://doi.org/10.1038/ng.3178.