SearcHPV: A novel approach to identify and assemble human papillomavirus-host genomic integration events in cancer

Virusplot: a web server for viral integration analysis and visualization

The integration of viral DNA into the human genome is a critical event in the pathogenesis of various cancers. This process leads to genomic instability, disrupts cellular regulatory mechanisms, and activates oncogenes or inactivates tumor suppressor genes. Despite significant advancements in genome sequencing technologies, there remains a notable lack of computational tools, particularly web-based applications, specifically designed for viral integration analysis and visualization. To address this gap, we present virusPlot, a web server with the following functional modules: (i) automatic retrieval of virus genome sequences and their annotation; (ii) visualization of virus integration locations and read counts through a graphical representation that links viral and host genome integration sites, facilitating the interpretation of integration patterns; (iii) analysis of virus integration hotspots using Fisher's exact test; and (iv) integration of various functions into an interactive web platform via shinyapp. VirusPlot efficiently processes and visualizes integration data from viruses and host genomes, providing researchers with an intuitive and user-friendly analytical tool that simplifies the complexity of virus integration analysis.

The role of HR-HPV integration in the progression of premalignant lesions into different cancer types

High-risk human papillomavirus (HR-HPV) is associated with the development of different types of cancer, such as cervical, head and neck (including oral, laryngeal, and oropharyngeal), vulvar, vaginal, penile, and anal cancers. The progression of premalignant lesions to cancer depends on factors associated with the host cell and the different epithelia infected by HPV, such as basal cells of the flat epithelium and the cells of the squamocolumnar transformation zone (STZ) found in the uterine cervix and the anal canal, which is rich in heparan sulfate proteoglycans and integrin-like receptors. On the other hand, factors associated with the viral genotype, infection with multiple viruses, viral load, viral persistence, and type of integration determine the viral breakage pattern and the sites at which the virus integrates into the host cell genome (introns, exons, intergenic regions), inducing the loss of function of tumor suppressor genes and increasing oncogene expression. This review describes the role of viral integration and the molecular mechanisms induced by HR-HPV in different types of tissues. The purpose of this review is to identify the common factors associated with the role of integration events in the progression of premalignant lesions in different types of cancer.

Identification of transcriptionally-active human papillomavirus integrants through nanopore sequencing reveals viable targets for gene therapy against cervical cancer

Integration of the human papillomavirus (HPV) genome into the cellular genome is a key event that leads to constitutive expression of viral oncoprotein E6/E7 and drives the progression of cervical cancer. However, HPV integration patterns differ on a case-by-case basis among related malignancies. Next-generation sequencing technologies still face challenges for interrogating HPV integration sites. In this study, utilizing Nanopore long-read sequencing, we identified 452 and 108 potential integration sites from the cervical cancer cell lines (CaSki and HeLa) and five tissue samples, respectively. Based on long Nanopore chimeric reads, we were able to analyze the methylation status of the HPV long control region (LCR), which controls oncogene E6/E7 expression, and to identify transcriptionally-active integrants among the numerous integrants. As a proof of concept, we identified an active HPV integrant in between RUNX2 and CLIC5 on chromosome 6 in the CaSki cell line, which was supported by ATAC-seq, H3K27Ac ChIP-seq, and RNA-seq analysis. Knockout of the active HPV integrant, by the CRISPR/Cas9 system, dramatically crippled cell proliferation and induced cell senescence. In conclusion, identifying transcriptionally-active HPV integrants with Nanopore sequencing can provide viable targets for gene therapy against HPV-associated cancers.

Transdifferentiation of cervical squamous cell carcinoma with ERBB2 amplification to adenocarcinoma: whole genome sequence analysis and successful control by anti-HER2 therapy

Cancer cells sometimes transdifferentiate into different histological type(s) and tumors with multiple histological types can share a common ancestor cell. However, diagnosis of the origin of multiple tumor lesions with different histological features remains a clinical challenge. A 45-year-old woman with a history of cervical squamous cell carcinoma (CeSq) presented with abdominal pain and vomiting. A surgical operation revealed an ileal tumor and a peritoneal nodule with a small amount of ascites. A histological examination of the ileal tumor demonstrated squamous cell carcinoma, which was consistent with metastasis of cervical cancer, while that of the nodule and ascites showed adenocarcinoma. Whole genome sequencing (WGS) of the CeSq, ileal squamous cell carcinoma (SiSq), and peritoneal adenocarcinoma (PeAd) demonstrated that ERBB2 was commonly amplified in all lesions. Additionally, HPV-16 genome sequences were identified at identical genomic loci in these lesions. A trajectory analysis corroborated that SiSq and PeAd had a shared origin and developed simultaneously at each metastatic site. These results indicate that a subpopulation of the CeSq had transdifferentiated into adenocarcinoma in our patient. Anti-HER2 therapy showed marked effects on the recurrent disease. Our case demonstrates the plasticity of tumor cells and reinforces the potential roles of WGS in the implementation of precision oncology.

Analysis of Human Papilloma Virus Content and Integration in Mucoepidermoid Carcinoma

Mucoepidermoid Carcinomas (MEC) represent the most common malignancies of salivary glands. Approximately 50% of all MEC cases are known to harbor CRTC1/3-MAML2 gene fusions, but the additional molecular drivers remain largely uncharacterized. Here, we sought to resolve controversy around the role of human papillomavirus (HPV) as a potential driver of mucoepidermoid carcinoma. Bioinformatics analysis was performed on 48 MEC transcriptomes. Subsequent targeted capture DNA sequencing was used to annotate HPV content and integration status in the host genome. HPV of any type was only identified in 1/48 (2%) of the MEC transcriptomes analyzed. Importantly, the one HPV16+ tumor expressed high levels of p16, had high expression of HPV16 oncogenes E6 and E7, and displayed a complex integration pattern that included breakpoints into 13 host genes including PIK3AP1, HIPI, OLFM4,SIRT1, ARAP2, TMEM161B-AS1, and EPS15L1 as well as 9 non-genic regions. In this cohort, HPV is a rare driver of MEC but may have a substantial etiologic role in cases that harbor the virus. Genetic mechanisms of host genome integration are similar to those observed in other head and neck cancers.

HPV Integration Site Mapping: A Rapid Method of Viral Integration Site (VIS) Analysis and Visualization Using Automated Workflows in CLC Microbial Genomics

Human papillomavirus (HPV) integration within the host genome may contribute to carcinogenesis through various disruptive mechanisms. With next-generation sequencing (NGS), identification of viral and host genomic breakpoints and chimeric sequences are now possible. However, a simple, streamlined bioinformatics workflow has been non-existent until recently. Here, we tested two new, automated workflows in CLC Microbial Genomics, i.e., Viral Hybrid Capture (VHC) Data Analysis and Viral Integration Site (VIS) Identification for software performance and efficiency. The workflows embedded with HPV and human reference genomes were used to analyze a publicly available NGS dataset derived from pre- and cancerous HPV+ cervical cytology of 21 Gabonese women. The VHC and VIS workflow median runtimes were 19 and 7 min per sample, respectively. The VIS dynamic graphical outputs included read mappings, virus-host genomic breakpoints, and virus-host integration circular plots. Key findings, including disrupted and nearby genes, were summarized in an auto-generated report. Overall, the VHC and VIS workflows proved to be a rapid and accurate means of localizing viral-host integration site(s) and identifying disrupted and neighboring human genes. Applying HPV VIS-mapping to pre- or invasive tumors will advance our understanding of viral oncogenesis and facilitate the discovery of prognostic biomarkers and therapeutic targets.

Upregulated miRNAs on the TP53 and RB1 Binding Seedless Regions in High-Risk HPV-Associated Penile Cancer

Cancer development by the human papillomavirus (HPV) infection can occur through the canonical HPV/p53/RB1 pathway mediated by the E2/E6/E7 viral oncoproteins. During the transformation process, HPV inserts its genetic material into host Integration Sites (IS), affecting coding genes and miRNAs. In penile cancer (PeCa) there is limited data on the miRNAs that regulate mRNA targets associated with HPV, such as the TP53 and RB1 genes. Considering the high frequency of HPV infection in PeCa patients in Northeast Brazil, global miRNA expression profiling was performed in high-risk HPV-associated PeCa that presented with TP53 and RB1 mRNA downregulated expression. The miRNA expression profile of 22 PeCa tissue samples and five non-tumor penile tissues showed 507 differentially expressed miRNAs: 494 downregulated and 13 upregulated (let-7a-5p, miR-130a-3p, miR-142-3p, miR-15b-5p miR-16-5p, miR-200c-3p, miR-205-5p, miR-21-5p, miR-223-3p, miR-22-3p, miR-25-3p, miR-31-5p and miR-93-5p), of which 11 were identified to be in HPV16-IS and targeting TP53 and RB1 genes. One hundred and thirty-one and 490 miRNA binding sites were observed for TP53 and RB1, respectively, most of which were in seedless regions. These findings suggest that up-regulation of miRNA expression can directly repress TP53 and RB1 expression by their binding sites in the non-canonical seedless regions.

Mutations in the telomerase reverse transcriptase promoter and PIK3CA gene are common events in penile squamous cell carcinoma of Italian and Ugandan patients

Penile carcinoma develops either through human papillomavirus (HPV) related or unrelated carcinogenic pathways. Genetic alterations and nucleotide changes in coding regions (ie, TP53, CDKN2A, PIK3CA and NOTCH1) are main cancer driver events either in HPV positive or in HPV negative tumours. We investigated the presence of hotspot nucleotide mutations in TERT promoter (TERTp) and PIK3CA exon 9 and their relationship with HPV status in 69 penile cancer cases from Italian and Ugandan patients. Genetic variations and viral sequences have been characterised by end‐point polymerase chain reaction (PCR) and Sanger sequencing. The mutant allele frequencies (MAFs) of TERTp −124A/−146A and PIK3CA E545K have been determined by droplet digital PCR (ddPCR) assays. The results showed that TERTp mutations are highly prevalent in penile carcinoma (53.6%) and significantly more frequent in HPV negative (67.6%) than HPV positive (32.4%) cases (P = .0482). PIK3CA mutations were similarly distributed in virus‐related and unrelated cases (25.9% and 26.7%, respectively) and coexisted with TERTp changes in 15.8% of penile carcinoma samples. Notably, MAFs of co‐occurring mutations were frequently discordant indicating that PIK3CA E545K nucleotide changes are subsequent genetic events occurring in subclones of TERTp mutated cells. The frequencies of TERTp and PIK3CA mutations were higher among Italian compared to Ugandan cases and inversely correlated with the HPV status. In conclusion, TERTp mutations are very common in penile carcinoma and their coexistence with PIK3CA in a substantial number of cases may represent a novel oncogenic synergy relevant for patient stratification and use of therapeutic strategies against new actionable targets.

Early HPV ctDNA Kinetics and Imaging Biomarkers Predict Therapeutic Response in p16+ Oropharyngeal Squamous Cell Carcinoma

PURPOSE

In locally advanced p16+ oropharyngeal squamous cell carcinoma (OPSCC), (i) to investigate kinetics of human papillomavirus (HPV) circulating tumor DNA (ctDNA) and association with tumor progression after chemoradiation, and (ii) to compare the predictive value of ctDNA to imaging biomarkers of MRI and FDG-PET.

EXPERIMENTAL DESIGN

Serial blood samples were collected from patients with AJCC8 stage III OPSCC (n = 34) enrolled on a randomized trial: pretreatment; during chemoradiation at weeks 2, 4, and 7; and posttreatment. All patients also had dynamic-contrast-enhanced and diffusion-weighted MRI, as well as FDG-PET scans pre-chemoradiation and week 2 during chemoradiation. ctDNA values were analyzed for prediction of freedom from progression (FFP), and correlations with aggressive tumor subvolumes with low blood volume (TVLBV) and low apparent diffusion coefficient (TVLADC), and metabolic tumor volume (MTV) using Cox proportional hazards model and Spearman rank correlation.

RESULTS

Low pretreatment ctDNA and an early increase in ctDNA at week 2 compared with baseline were significantly associated with superior FFP (P < 0.02 and P < 0.05, respectively). At week 4 or 7, neither ctDNA counts nor clearance were significantly predictive of progression (P = 0.8). Pretreatment ctDNA values were significantly correlated with nodal TVLBV, TVLADC, and MTV pre-chemoradiation (P < 0.03), while the ctDNA values at week 2 were correlated with these imaging metrics in primary tumor. Multivariate analysis showed that ctDNA and the imaging metrics performed comparably to predict FFP.

CONCLUSIONS

Early ctDNA kinetics during definitive chemoradiation may predict therapy response in stage III OPSCC.

Causes and Consequences of HPV Integration in Head and Neck Squamous Cell Carcinomas: State of the Art

A constantly increasing incidence in high-risk Human Papillomaviruses (HPV)s driven head and neck squamous cell carcinomas (HNSCC)s, especially of oropharyngeal origin, is being observed. During persistent infections, viral DNA integration into the host genome may occur. Studies are examining if the physical status of the virus (episomal vs. integration) affects carcinogenesis and eventually has further-reaching consequences on disease progression and outcome. Here, we review the literature of the most recent five years focusing on the impact of HPV integration in HNSCCs, covering aspects of detection techniques used (from PCR up to NGS approaches), integration loci identified, and associations with genomic and clinical data. The consequences of HPV integration in the human genome, including the methylation status and deregulation of genes involved in cell signaling pathways, immune evasion, and response to therapy, are also summarized.