DNA Methylation and HPV-Associated Head and Neck Cancer

Global DNA Methylation Level in Tumour and Margin Samples in Relation to Human Papilloma Virus and Epstein-Barr Virus in Patients with Oropharyngeal and Oral Squamous Cell Carcinomas

BACKGROUND

Aberrant DNA methylation is a common epigenetic modification in cancers, including oropharyngeal squamous cell carcinoma (OPSCC) and oral squamous cell carcinoma (OSCC). Therefore, the analysis of methylation levels appears necessary to improve cancer therapy and prognosis.

METHODS

The enzyme-linked immunosorbent assay (ELISA) was used to analyse global DNA methylation levels in OPSCC and OSCC tumours and the margin samples after DNA isolation. HPV detection was conducted by hybridisation using GenoFlow HPV Array Test Kits (DiagCor Bioscience Inc., Hong Kong, China). EBV detection was performed using real-time PCR with an EBV PCR Kit (EBV/ISEX/100, GeneProof, Brno, Czech Republic).

RESULTS

OPSCC tumour samples obtained from women showed lower global DNA methylation levels than those from men (1.3% vs. 3.5%, p = 0.049). The margin samples from OPSCC patients with HPV and EBV coinfection showed global DNA methylation lower than those without coinfection (p = 0.042). G3 tumours from OSCC patients had significantly lower levels of global DNA methylation than G2 tumours (0.98% ± 0.74% vs. 3.77% ± 4.97%, p = 0.010). Additionally, tumours from HPV-positive OSCC patients had significantly lower global DNA methylation levels than those from HPV-negative patients (p = 0.013). In the margin samples, we observed a significant negative correlation between global DNA methylation and the N stage of OSCC patients (rS = -0.33, p = 0.039). HPV-positive OPSCC patients had higher global DNA methylation levels than HPV-positive OSCC patients (p = 0.015).

CONCLUSION

We confirmed that methylation could be changed in relation to viral factors, such as HPV and EBV, as well as clinical and demographical parameters.

Promoter DNA methylation patterns in oral, laryngeal and oropharyngeal anatomical regions are associated with tumor differentiation, nodal involvement and survival

Differentially methylated regions (DMRs) can be used as head and neck squamous cell carcinoma (HNSCC) diagnostic, prognostic and therapeutic targets in precision medicine workflows. DNA from 21 HNSCC and 10 healthy oral tissue samples was hybridized to a genome-wide tiling array to identify DMRs in a discovery cohort. Downstream analyses identified differences in promoter DNA methylation patterns in oral, laryngeal and oropharyngeal anatomical regions associated with tumor differentiation, nodal involvement and survival. Genome-wide DMR analysis showed 2,565 DMRs common to the three subsites. A total of 738 DMRs were unique to laryngeal cancer (n=7), 889 DMRs were unique to oral cavity cancer (n=10) and 363 DMRs were unique to pharyngeal cancer (n=6). Based on the genome-wide analysis and a Gene Ontology analysis, 10 candidate genes were selected to test for prognostic value and association with clinicopathological features. TIMP3 was associated with tumor differentiation in oral cavity cancer (P=0.039), DAPK1 was associated with nodal involvement in pharyngeal cancer (P=0.017) and PAX1 was associated with tumor differentiation in laryngeal cancer (P=0.040). A total of five candidate genes were selected, DAPK1, CDH1, PAX1, CALCA and TIMP3, for a prevalence study in a larger validation cohort: Oral cavity cancer samples (n=42), pharyngeal cancer tissues (n=25) and laryngeal cancer samples (n=52). PAX1 hypermethylation differed across HNSCC anatomic subsites (P=0.029), and was predominantly detected in laryngeal cancer. Kaplan-Meier survival analysis (P=0.043) and Cox regression analysis of overall survival (P=0.001) showed that DAPK1 methylation is associated with better prognosis in HNSCC. The findings of the present study showed that the HNSCC subsites oral cavity, pharynx and larynx display substantial differences in aberrant DNA methylation patterns, which may serve as prognostic biomarkers and therapeutic targets.

Human Papillomavirus 16 DNA Methylation Patterns and Investigation of Integration Status in Head and Neck Cancer Cases

Persistent high-risk human papillomavirus (HPV) infection is a pivotal factor in the progression of cervical cancer. In recent years, an increasing interest has emerged in comprehending the influence of HPV on head and neck squamous cell carcinoma (HNSCC). Notably, it is well established that HPV-associated HNSCC show cases with distinct molecular and clinical attributes compared to HPV-negative cases. The present study delves into the epigenetic landscape of HPV16, specifically its L1 gene and untranslated region (UTR), through pyrosequencing, while the HPV16 DNA physical status was evaluated using E2/E6 ratio analysis in HPV16-positive HNSCC FFPE biopsies. Our findings reveal substantial methylation across six sites within the HPV16 L1 gene and seven sites in the UTR. Specifically, methylation percentages of two L1 CpG sites (7136, 7145) exhibit significant associations with tumor histological grade (p < 0.01), while proving concurrent methylation across multiple sites. The HPV16 DNA physical status was not correlated with the methylation of viral genome or tumor characteristics. This is the first study that examines epigenetic modifications and the HPV16 DNA physical status in Greek HNSCC patients. Our findings suggest an orchestrated epigenetic modulation among specific sites, impacting viral gene expression and intricate virus-host interactions.

Role of DNA Methylation Profiles as Potential Biomarkers and Novel Therapeutic Targets in Head and Neck Cancer

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common cancer worldwide and is associated with high mortality. The main reasons for treatment failure are a low rate of early diagnosis, high relapse rates, and distant metastasis with poor outcomes. These are largely due to a lack of diagnostic, prognostic, and predictive biomarkers in HNSCC. DNA methylation has been demonstrated to play an important role in the pathogenesis of HNSCC, and recent studies have also valued DNA methylation as a potential biomarker in HNSCC. This review summarizes the current knowledge on DNA methylation profiles in HPV-positive and HPV-negative HNSCC and how these may contribute to the pathogenesis of HNSCC. It also summarizes the potential value of DNA methylation as a biomarker in the diagnosis, prognosis, and prediction of the response to therapy. With the recent immunotherapy era in head and neck treatment, new strategies to improve immune responses by modulating TIMEs have been intensely investigated in early-phase trials. Therefore, this study additionally summarizes the role of DNA methylation in the regulation of TIMEs and potential predictive immunotherapy response biomarkers. Finally, this study reviews ongoing clinical trials using DNA methylation inhibitors in HNSCC.

HPV16 Impacts NHERF2 Expression in Oropharyngeal Cancers

Infection with human papillomaviruses (HPVs), in particular with HPV type 16, is now considered to be a key risk factor for the development of a subset of oropharyngeal squamous cell carcinomas (OPSCC) that show different epidemiological, clinical, and prognostic characteristics from HPV-negative (HPV-) OPSCCs. So far, extensive research efforts aiming to distinguish these two distinct entities have not identified specific biomarkers, nor led to different therapies. Previous research has shown that HPV16 E6 oncoprotein binds NHERF2, inducing its proteasomal degradation, and consequently increasing cell proliferation; we therefore aimed to investigate how this might be reflected in human histological samples. We analyzed NHERF2 expression patterns in HPV16-positive (HPV16+) and HPV- OPSCC samples, to investigate any potential differences in NHERF2 pattern. Interestingly, we observed a statistically significant decrease in NHERF2 levels in HPV16+ and poorly differentiated HPV- OPSCCs, compared with healthy tissue. Furthermore, we observed a significant reduction in the percentage of NHERF2 immunoreactive cancer cells in HPV16+ tumors, compared with well and moderately differentiated HPV- OPSCCs, suggesting the importance of 16E6's targeting of NHERF2 in HPV-driven oncogenesis in the head and neck area.

Whole-genome CpG-resolution DNA Methylation Profiling of HNSCC Reveals Distinct Mechanisms of Carcinogenesis for Fine-scale HPV+ Cancer Subtypes

DNA methylation is a vital early step in carcinogenesis. Most findings of aberrant DNA methylation in head and neck squamous cell carcinomas (HNSCC) are array based with limited coverage and resolution, and mainly explored by human papillomavirus (HPV) status, ignoring the high heterogeneity of this disease. In this study, we performed whole-genome bisulfite sequencing on a well-studied HNSCC cohort (n = 36) and investigated the methylation changes between fine-scaled HNSCC subtypes in relation to genomic instability, repetitive elements, gene expression, and key carcinogenic pathways. The previously observed hypermethylation phenotype in HPV-positive (HPV+) tumors compared with HPV-negative tumors was robustly present in the immune-strong (IMU) HPV+ subtype but absent in the highly keratinized (KRT) HPV+ subtype. Methylation levels of IMU tumors were significantly higher in repetitive elements, and methylation showed a significant correlation with genomic stability, consistent with the IMU subtype having more genomic stability and better prognosis. Expression quantitative trait methylation (cis-eQTM) analysis revealed extensive functionally-relevant differences, and differential methylation pathway analysis recapitulated gene expression pathway differences between subtypes. Consistent with their characteristics, KRT and HPV-negative tumors had high regulatory potential for multiple regulators of keratinocyte differentiation, which positively correlated with an expression-based keratinization score. Together, our findings revealed distinct mechanisms of carcinogenesis between subtypes in HPV+ HNSCC and uncovered previously ignored epigenomic differences and clinical implications, illustrating the importance of fine-scale subtype analysis in cancer.

Significance

This study revealed that the previously observed hypermethylation of HPV(+) HNSCC is due solely to the IMU subtype, illustrating the importance of fine-scale subtype analysis in such a heterogeneous disease. Particularly, IMU has significantly higher methylation of transposable elements, which can be tested as a prognosis biomarker in future translational studies.

erarchical tumor artificial intelligence classifier traces tissue of origin and tumor type in primary and metastasized tumors using DNA methylation

Human cancers are heterogenous by their cell composition and origination site. Cancer metastasis generates the conundrum of the unknown origin of migrated tumor cells. Tracing tissue of origin and tumor type in primary and metastasized cancer is vital for clinical significance. DNA methylation alterations play a crucial role in carcinogenesis and mark cell fate differentiation, thus can be used to trace tumor tissue of origin. In this study, we employed a novel tumor-type-specific hierarchical model using genome-scale DNA methylation data to develop a multilayer perceptron model, HiTAIC, to trace tissue of origin and tumor type in 27 cancers from 23 tissue sites in data from 7735 tumors with high resolution, accuracy, and specificity. In tracing primary cancer origin, HiTAIC accuracy was 99% in the test set and 93% in the external validation data set. Metastatic cancers were identified with a 96% accuracy in the external data set. HiTAIC is a user-friendly web-based application through https://sites.dartmouth.edu/salaslabhitaic/. In conclusion, we developed HiTAIC, a DNA methylation-based algorithm, to trace tumor tissue of origin in primary and metastasized cancers. The high accuracy and resolution of tumor tracing using HiTAIC holds promise for clinical assistance in identifying cancer of unknown origin.

Hypermethylation of nc886 in HPV-positive oropharyngeal cancer and its clinical implications: An epigenome-wide association study

The incidence of oropharyngeal squamous cell carcinoma (OPSCC) has increased rapidly in the United States, driven by rising human papillomavirus (HPV) infections in the U.S. population. HPV-positive OPSCC patients have a better prognosis than HPV-negative patients. To gain insights into the unique biology of HPV(+) OPSCC that may contribute to its clinical behaviors, we performed a multi-stage epigenome-wide methylation profiling of leukocyte and tumor DNA in OPSCC patients and compared the methylation levels of CpG sites between HPV(+) and HPV(-) OPSCC patients. We identified and validated a significantly differentially methylated region (DMR) of 1,355 bp encompassing non-coding RNA 886 (nc886) gene and its promoter region. Nc886 is hypermethylated in both leukocytes and tumor DNA of HPV(+) OPSCC patients. Homozygous knockout of nc886 by CRISPR-Cas9 in head and neck cell lines was lethal, but nc886 could be knocked out on the background of protein kinase R (PKR) knockout. Our data suggest that HPV induces nc886 hypermethylation, and nc886 acts as both a viral sensor and a tumor sensor in OPSCC patients and contribute to the better prognosis of HPV(+) OPSCC patients. Nc886 may become a therapeutic target in OPSCC.

The HPV Induced Cancer Resource (THInCR): a Suite of Tools for Investigating HPV-Dependent Human Carcinogenesis

Human papillomaviruses (HPVs) are highly infectious and cause the most common sexually transmitted viral infections. They induce hyperproliferation of squamous epithelial tissue, often forming warts. Virally encoded proteins reprogram gene expression and cell growth to create an optimal environment for viral replication. In addition to their normal roles in infection, functional alterations induced by viral proteins establish conditions that frequently contribute to human carcinogenesis. In fact, ~5% of human cancers are caused by HPVs, with virtually all cervical squamous cell carcinomas (CESC) and an increasing number of head and neck squamous cell carcinomas (HNSC) attributed to HPV infection. The Cancer Genome Atlas (TCGA) molecularly characterized thousands of primary human cancer samples in many cancer types, including CESC and HNSC, and created a comprehensive atlas of genomic, epigenomic, and transcriptomic data. This publicly available genome-wide information provides an unprecedented opportunity to expand the knowledge of the role that HPV plays in human carcinogenesis. While many tools exist to mine these data, few, if any, focus on the comparison of HPV-positive cancers with their HPV-negative counterparts or adjacent normal control tissue. We have constructed a suite of web-based tools, The HPV Induced Cancer Resource (THInCR), to utilize TCGA data for research related to HPV-induced CESC and HNSC. These tools allow investigators to gain greater biological and medical insights by exploring the impacts of HPV on cellular gene expression (mRNA and microRNA), altered gene methylation, and associations with patient survival and immune landscape features. These tools are accessible at https://thincr.ca/. IMPORTANCE The suite of analytical tools of THInCR provides the opportunity to investigate the roles that candidate target genes identified in cell lines or other model systems contribute to in actual HPV-dependent human cancers and is based on large-scale TCGA data sets. Expression of target genes, including both mRNA and microRNA, can be correlated with HPV gene expression, epigenetic changes in DNA methylation, patient survival, and numerous immune features, like leukocyte infiltration, interferon gamma response, T cell response, etc. Data from these analyses may immediately provide evidence to validate in vitro observations, reveal insights into mechanisms of virus-mediated alterations in cell growth, behavior, gene expression, and innate and adaptive immunity and may help hypothesis generation for further investigations.

A Pan-Cancer Analysis Reveals the Prognostic and Immunotherapeutic Value of Stanniocalcin-2 (STC2)

Background: Stanniocalcin-2 (STC2) is a secreted glycoprotein which plays an important role in regulating the homeostasis of calcium, glucose homeostasis, and phosphorus metastasis. Accumulating evidence suggests that STC2 is implicated in cancer mechanisms. However, the effects of STC2 on cancer development and progression across pan-cancer are not yet completely known.

Methods: Data were downloaded from The Cancer Genome Atlas database to obtain differentially expressed genes significantly associated with prognosis (key genes). A gene was selected for subsequent correlation studies by integrating the significance of prognosis and the time-dependent ROC curve. Gene expression of different tumor types was analyzed based on the UCSC XENA website. Furthermore, our study investigated the correlation of STC2 expression between prognosis, immune cell infiltration, immune checkpoint genes (ICGs), mismatch repair genes (MMRs), tumor mutation burden (TMB), microsatellite instability (MSI), and drug sensitivity in various malignant tumors. Gene set enrichment analysis (GSEA) was conducted for correlated genes of STC2 to explore potential mechanisms.

Results: A total of 3,429 differentially expressed genes and 397 prognosis-related genes were identified from the TCGA database. Twenty-six key genes were found by crossing the former and the latter, and the highest risk gene, STC2, was selected for subsequent correlation studies. STC2 had good diagnostic performance for HNSCC, and was closely related to the survival status and clinicopathological stage of HNSCC patients. In pan-cancer analysis, STC2 was upregulated in 20 cancers and downregulated in seven cancers. STC2 overexpression was overall negatively correlated with overall survival, disease-free survival, disease-specific survival, and progress-free survival. STC2 was profoundly correlated with the tumor immune microenvironment, including immune cell infiltration, ICGs, MMRs, TMB, and MSI. Moreover, STC2 was significantly negatively correlated with the sensitivity or resistance of multiple drugs.

Conclusion: STC2 was a potential prognostic biomarker for pan-cancer and a new immunotherapy target.

Novel Systemic Treatment Modalities Including Immunotherapy and Molecular Targeted Therapy for Recurrent and Metastatic Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinomas (HNSCCs) are the sixth most common cancers worldwide. More than half of patients with HNSCC eventually experience disease recurrence and/or metastasis, which can threaten their long-term survival. HNSCCs located in the oral cavity and larynx are usually associated with tobacco and/or alcohol use, whereas human papillomavirus (HPV) infection, particularly HPV16 infection, is increasingly recognized as a cause of oropharyngeal HNSCC. Despite clinical, histologic, and molecular differences between HPV-positive and HPV-negative HNSCCs, current treatment approaches are the same. For recurrent disease, these strategies include chemotherapy, immunotherapy with PD-1-inhibitors, or a monoclonal antibody, cetuximab, that targets epidermal growth factor; these therapies can be administered either as single agents or in combination. However, these treatment strategies carry a high risk of toxic side effects; therefore, more effective and less toxic treatments are needed. The landscape of HNSCC therapy is changing significantly; numerous clinical trials are underway to test novel therapeutic options like adaptive cellular therapy, antibody-drug conjugates, new targeted therapy agents, novel immunotherapy combinations, and therapeutic vaccines. This review helps in understanding the various developments in HNSCC therapy and sheds light on the path ahead in terms of further research in this field.

A novel necroptosis-related LncRNA signature for prediction of prognosis and therapeutic responses of head and neck squamous cell carcinoma

Background: Long non-coding RNAs (lncRNAs) play an essential role in the occurrence and prognosis of tumors, and it has great potential as biomarkers of tumors. However, the roles of Necroptosis-related lncRNA (NRLs) in Head and neck squamous cell carcinoma (HNSCC) remain elusive. Methods: We comprehensively analyzed the gene expression and clinical information of 964 HNSCC in four cohorts. LASSO regression was utilized to construct a necroptosis-related lncRNA prognosis signature (NLPS). We used univariate and multivariate regression to assess the independent prognostic value of NLPS. Based on the optimal cut-off, patients were divided into high- and low-risk groups. In addition, the immune profile, multi-omics alteration, and pharmacological landscape of NLPS were further revealed. Results: A total of 21 NRLs associated with survival were identified by univariate regression in four cohorts. We constructed and validated a best prognostic model (NLPS). Compared to the low-risk group, patients in the high group demonstrated a more dismal prognosis. After adjusting for clinical features by multivariate analysis, NLPS still displayed independent prognostic value. Additionally, further analysis found that patients in the low-risk group showed more abundant immune cell infiltration and immunotherapy response. In contrast, patients in the high-risk group were more sensitive to multiple chemotherapeutic agents. Conclusion: As a promising tool, the establishment of NLPS provides guidance and assistance in the clinical management and personalized treatment of HNSCC.

Genetic and methylation status of CDKN2A (p14ARF/p16INK4A) and TP53 genes in recurrent respiratory papillomatosis

Recurrent respiratory papillomatosis (RRP) is a rare and chronic disease affecting the upper airway with papillomatous lesions caused by the human papillomavirus (HPV) infection, especially HPV-6 and/or HPV-11 types. Little is known about the genetic and epigenetic drivers in RRP pathophysiology. For this purpose, we analyzed 27 papillomatous lesions from patients with RRP to evaluate somatic mutations and methylation status in CDKN2A (p14^ARF/p16^INK4A) and TP53, which are key tumor suppressor genes for the cell cycle control. Sanger sequencing analysis revealed one somatic mutation in TP53 (c.733_734insA) and four mutations in CDKN2A (c.-30G > T, c.29_30insA, c.69delT, and c.300C > A). These mutations were observed in 10 patients, 6 of which carried double mutation. Furthermore, 50% (5/10) of these patients carrying somatic mutations had RRP severity, representing 62.5% (5/8) of the severity cases in this study, albeit no significant association was found between somatic mutations and disease severity. Methylation-specific polymerase chain reaction assays revealed p14^ARF promoter hypermethylation in 100% of cases, followed by TP53 (96.3%) and p16^INK4A (55.6%), suggesting the influence of HPV in the DNA methylation machinery. In conclusion, somatic mutations were not common events identified in patients with RRP. However, epigenetic modulation by high methylation rates, particularly for the p14^ARF/TP53 pathway, seems to be in the course of RRP development.

Epigenetic Dysregulations in Merkel Cell Polyomavirus-Driven Merkel Cell Carcinoma

Merkel cell polyomavirus (MCPyV) is a small DNA virus with oncogenic potential. MCPyV is the causative agent of Merkel Cell Carcinoma (MCC), a rare but aggressive tumor of the skin. The role of epigenetic mechanisms, such as histone posttranslational modifications (HPTMs), DNA methylation, and microRNA (miRNA) regulation on MCPyV-driven MCC has recently been highlighted. In this review, we aim to describe and discuss the latest insights into HPTMs, DNA methylation, and miRNA regulation, as well as their regulative factors in the context of MCPyV-driven MCC, to provide an overview of current findings on how MCPyV is involved in the dysregulation of these epigenetic processes. The current state of the art is also described as far as potentially using epigenetic dysregulations and related factors as diagnostic and prognostic tools is concerned, in addition to targets for MCPyV-driven MCC therapy. Growing evidence suggests that the dysregulation of HPTMs, DNA methylation, and miRNA pathways plays a role in MCPyV-driven MCC etiopathogenesis, which, therefore, may potentially be clinically significant for this deadly tumor. A deeper understanding of these mechanisms and related factors may improve diagnosis, prognosis, and therapy for MCPyV-driven MCC.

Role of Virus-Induced Host Cell Epigenetic Changes in Cancer

The tumor viruses human T-lymphotropic virus 1 (HTLV-1), hepatitis C virus (HCV), Merkel cell polyomavirus (MCPyV), high-risk human papillomaviruses (HR-HPVs), Epstein-Barr virus (EBV), Kaposi’s sarcoma-associated herpes virus (KSHV) and hepatitis B virus (HBV) account for approximately 15% of all human cancers. Although the oncoproteins of these tumor viruses display no sequence similarity to one another, they use the same mechanisms to convey cancer hallmarks on the infected cell. Perturbed gene expression is one of the underlying mechanisms to induce cancer hallmarks. Epigenetic processes, including DNA methylation, histone modification and chromatin remodeling, microRNA, long noncoding RNA, and circular RNA affect gene expression without introducing changes in the DNA sequence. Increasing evidence demonstrates that oncoviruses cause epigenetic modifications, which play a pivotal role in carcinogenesis. In this review, recent advances in the role of host cell epigenetic changes in virus-induced cancers are summarized.

HPV Infection Leaves a DNA Methylation Signature in Oropharyngeal Cancer Affecting Both Coding Genes and Transposable Elements

HPV oncoproteins can modulate DNMT1 expression and activity, and previous studies have reported both gene-specific and global DNA methylation alterations according to HPV status in head and neck cancer. However, validation of these findings and a more detailed analysis of the transposable elements (TEs) are still missing. Here we performed pyrosequencing to evaluate a 5-CpG methylation signature and Line1 methylation in an oropharyngeal squamous cell carcinoma (OPSCC) cohort. We further evaluated the methylation levels of the TEs, their correlation with gene expression and their impact on overall survival (OS) using the TCGA cohort. In our dataset, the 5-CpG signature distinguished HPV-positive and HPV-negative OPSCC with 66.67% sensitivity and 84.33% specificity. Line1 methylation levels were higher in HPV-positive cases. In the TCGA cohort, Line1, Alu and long terminal repeats (LTRs) showed hypermethylation in a frequency of 60.5%, 58.9% and 92.3%, respectively. ZNF541 and CCNL1 higher expression was observed in HPV-positive OPSCC, correlated with lower methylation levels of promoter-associated Alu and LTR, respectively, and independently associated with better OS. Based on our findings, we may conclude that a 5-CpG methylation signature can discriminate OPSCC according to HPV status with high accuracy and TEs are differentially methylated and may regulate gene expression in HPV-positive OPSCC.