Identification of novel methylation markers in HPV-associated oropharyngeal cancer: genome-wide discovery, tissue verification and validation testing in ctDNA

The molecular characteristics of recurrent/metastatic HPV-positive head and neck squamous cell carcinoma: A systematic review of the literature

OBJECTIVES

About 17% of patients with human papillomavirus (HPV)-positive head and neck squamous cell carcinoma (HNSCC), which is mainly comprised of oropharyngeal SCC (OPSCC), will experience disease recurrence, which is often considered incurable when manifested at a metastatic and/or recurrent stage. We conducted a critical qualitative systematic review. Our objectives were to provide an overview of the molecular landscape of recurrent/metastatic HPV-positive HNSCC as well as novel molecular biomarkers.

DESIGN

A literature review was conducted to identify studies reporting on the molecular characteristics of recurrent/metastatic HPV-positive HNSCC, novel molecular biomarkers and treatment options. The reviews of abstracts, full articles, and revision of the included studies, followed by data extraction and quality assessment were performed by three independent assessors. All primary literature, such as retrospective, prospective, and clinical trials as well as basic research studies were considered, and the final search was conducted at the end of February 2023. The level of evidence was rated using the guidelines published by the Oxford Centre for Evidence-based Medicine and quality was assessed using the Newcastle-Ottawa Scale criteria.

RESULTS AND CONCLUSIONS

The literature search resulted in the identification of 1991 articles. A total of 181 full articles were screened, and 66 articles were included in this analysis. Several studies reported that recurrent/metastatic HPV-positive HNSCC had higher rates of TP53 mutation and were genomically similar to HPV-negative HNSCC. The detection of circulating tumour tissue-modified HPV DNA (ctHPVDNA) as a specific biomarker has shown promising results for monitoring treatment response and recurrence in the subset of HPV-positive HNSCC. In addition, evidence for targeted therapy in recurrent/metastatic HPV-positive HNSCC has emerged, including agents that inhibit overexpressed EGFR. Studies of combination immunotherapy are also underway. Our review outlines the latest evidence on the distinct molecular profiles of recurrent/metastatic HPV-positive HNSCC as well as the clinical potential of ctHPVDNA testing in routine practice. More controlled and longitudinal studies are needed to identify additional molecular targets and to assess the performance and benefits of novel molecular biomarkers in clinical practice.

Research advances of MAL family members in tumorigenesis and tumor progression (Review)

The myelin and lymphocyte protein (MAL) family is a novel gene family first identified and characterized in 2002. This family is comprised of seven members, including MAL, MAL2, plasmolipin, MALL, myeloid differentiation‑associated marker (MYADM), MYADML2 and CMTM8, which are located on different chromosomes. In addition to exhibiting extensive activity during transcytosis, the MAL family plays a vital role in the neurological, digestive, respiratory, genitourinary and other physiological systems. Furthermore, the intimate association between MAL and the pathogenesis, progression and metastasis of malignancies, attributable to several mechanisms such as DNA methylation has also been elucidated. In the present review, an overview of the structural and functional properties of the MAL family and the latest research findings regarding the relationship between several MAL members and various cancers is provided. Furthermore, the potential clinical and scientific significance of MAL is discussed and directions for future research are summarized.

Construction of a 12-gene prognostic model for colorectal cancer based on heat shock protein-related genes

Some heat shock proteins (HSPs) have been shown to influence tumor prognosis, but their prognostic significance in colorectal cancer (CRC) remains unclear. This study explored the prognostic significance of HSP-related genes in CRC. Transcriptional data and clinical information of CRC patients were obtained from The Cancer Genome Atlas (TCGA) database, and a literature search was conducted to identify HSP-related genes. Using Least Absolute Selection and Shrinkage Operator (LASSO) regression and univariate/multivariate Cox regression analyses, 12 HSP-related genes demonstrating significant associations with CRC survival were successfully identified and employed to formulate a predictive risk score model. The efficacy and precision of this model were validated utilizing TCGA and Gene Expression Omnibus (GEO) datasets, demonstrating its reliability in CRC prognosis prediction. gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed significant disparities between high- and low-risk groups in chromatin remodeling biological functions and neutrophil extracellular trap formation pathways. Single sample gene set enrichment analysis (ssGSEA) further revealed differences in immune cell types and immune functional status between the two risk groups. Differential analysis showed higher expression of immune checkpoints within the low-risk group, while the high-risk group exhibited notably higher Tumor Immune Dysfunction and Exclusion (TIDE) scores. Additionally, we predicted the sensitivity of different prognosis risk patients to various drugs, providing potential drug choices for tailored treatment. Combined, our study successfully crafted a novel CRC prognostic model that can effectively predict patient survival, immune landscape, and treatment response, providing important support and guidance for CRC patient prognosis.

Liquid biopsies based on cell-free DNA as a potential biomarker in head and neck cancer

In the era of 'precision medicine', liquid biopsies based on cell-free DNA (cfDNA) have emerged as a promising tool in the oncology field. cfDNA from cancer patients is a mixture of tumoral (ctDNA) and non-tumoral DNA originated from healthy, cancer and tumor microenvironmental cells. Apoptosis, necrosis, and active secretion from extracellular vesicles represent the main mechanisms of cfDNA release into the physiological body fluids. Focused on HNC, two main types of cfDNA can be identified: the circulating cfDNA (ccfDNA) and the salivary cfDNA (scfDNA). Numerous studies have reported on the potential of cfDNA analysis as potential diagnostic, prognostic, and monitoring biomarker for HNC. Thus, ctDNA has emerged as an attractive strategy to detect cancer specific genetic and epigenetic alterations including DNA somatic mutations and DNA methylation patterns. This review aims to provide an overview of the up-to-date studies evaluating the value of the analysis of total cfDNA, cfDNA fragment length, and ctDNA analysis at DNA mutation and methylation level in HNC patients.

A paradigm for post-embryonic Oct4 re-expression: E7-induced hydroxymethylation regulates Oct4 expression in cervical cancer

The Octamer-binding transcription factor-4 (Oct4) is upregulated in different malignancies, yet a paradigm for mechanisms of Oct4 post-embryonic re-expression is inadequately understood. In cervical cancer, Oct4 expression is higher in human papillomavirus (HPV)-related than HPV-unrelated cervical cancers and this upregulation correlates with the expression of the E7 oncogene. We have reported that E7 affects the Oct4-transcriptional output and Oct4-related phenotypes in cervical cancer, however, the underlying mechanism remains elusive. Here, we characterize the Oct4-protein interactions in cervical cancer cells via computational analyses and Mass Spectrometry and reveal that Methyl-binding proteins (MBD2 and MBD3), are determinants of Oct4-driven transcription. E7 triggers MBD2 downregulation and TET1 upregulation, thereby disrupting the methylation status of the Oct4 gene. This coincides with an increase in the total DNA hydroxymethylation leading to the re-expression of Oct4 in cervical cancer and likely affecting broader transcriptional patterns. Our findings reveal a previously unreported mechanism by which the E7 oncogene can regulate Oct4 re-expression and global transcriptional patterns by increasing DNA hydroxymethylation and lowering the barrier to cellular plasticity during carcinogenesis.

Identification of lactylation gene CALML5 and its correlated lncRNAs in cutaneous melanoma by machine learning

As a product of glycolysis, lactate contributes to cancer proliferation, immunosuppression, and metastasis via histone lactylation. However, the relationship between cutaneous melanoma (CM) and lactylation-associated genes and lncRNAs has remained unclear. In this study, 4 mechanism learning algorithms and integrated bioinformatic analyses were employed to identify the core lactylation-associated genes and lncRNAs. Subsequently, 2 risk signatures based on the hub lactylation-associated genes and lncRNAs were constructed for CM patients. As a result, CALML5 was identified as a core lactylation-associated gene in CM, and its expression was found to be associated with patients survival and immune infiltration, suggesting its relevance as a potential therapeutic target. Additionally, this study provided clarification on hub CALML5-associated lncRNAs in CM, offering insights into their roles in the disease. Meanwhile, 2 identified risk signatures were both strongly linked to the prognosis and cancer growth of CM, underscoring their potential as valuable prognostic indicators. Furthermore, mechanistic analyses suggested a significant association between the risk signature and the immune microenvironment in CM, highlighting potential immune-related implications in disease progression. In conclusion, we propose that lactylation-associated genes and lncRNAs hold promise as potential targets in CM. Moreover, our findings revealed a significant correlation between lactylation and the immune microenvironment, providing crucial insights for guiding individualized treatment strategies in CM.

A Network Landscape of HPVOPC Reveals Methylation Alterations as Significant Drivers of Gene Expression via an Immune-Mediated GPCR Signal

HPV-associated oropharynx carcinoma (HPVOPC) tumors have a relatively low mutational burden. Elucidating the relative contributions of other tumor alterations, such as DNA methylation alterations, alternative splicing events (ASE), and copy number variation (CNV), could provide a deeper understanding of carcinogenesis drivers in this disease. We applied network propagation analysis to multiple classes of tumor alterations in a discovery cohort of 46 primary HPVOPC tumors and 25 cancer-unaffected controls and validated our findings with TCGA data. We identified significant overlap between differential gene expression networks and all alteration classes, and this association was highest for methylation and lowest for CNV. Significant overlap was seen for gene clusters of G protein-coupled receptor (GPCR) pathways. HPV16-human protein interaction analysis identified an enriched cluster defined by an immune-mediated GPCR signal, including CXCR3 cytokines CXCL9, CXCL10, and CXCL11. CXCR3 was found to be expressed in primary HPVOPC, and scRNA-seq analysis demonstrated CXCR3 ligands to be highly expressed in M2 macrophages. In vivo models demonstrated decreased tumor growth with antagonism of the CXCR3 receptor in immunodeficient but not immunocompetent mice, suggesting that the CXCR3 axis can drive tumor proliferation in an autocrine fashion, but the effect is tempered by an intact immune system. In conclusion, methylation, ASE, and SNV alterations are highly associated with network gene expression changes in HPVOPC, suggesting that ASE and methylation alterations have an important role in driving the oncogenic phenotype. Network analysis identifies GPCR networks, specifically the CXCR3 chemokine axis, as modulators of tumor-immune interactions that may have proliferative effects on primary tumors as well as a role for immunosurveillance; however, CXCR3 inhibition should be used with caution, as these agents may both inhibit and stimulate tumor growth considering the competing effects of this cytokine axis. Further investigation is needed to explore opportunities for targeted therapy in this setting.

Circulating Tumor DNA in Human Papillomavirus-Mediated Oropharynx Cancer: Leveraging Early Data to Inform Future Directions

ABSTRACT

Circulating tumor DNA (ctDNA) has become an area of intense study in many solid malignancies including head and neck cancer. This is of particular interest for human papillomavirus-mediated oropharyngeal squamous cell carcinoma as this cohort of patients has excellent survival and is undergoing current clinical trials aimed at treatment de-escalation. Recent studies have demonstrated the prognostic implications of pretreatment ctDNA and the utility of monitoring ctDNA during and posttreatment; however, there is a need for a more critical understanding of ctDNA as it is beginning to be incorporated into clinical trials. This review discusses the current state of ctDNA in oropharynx cancer focusing on ctDNA kinetics and minimal residual disease detection and ends with a discussion of future applications.

Gene Expression Profiles of Methyltransferases and Demethylases Associated with Metastasis, Tumor Invasion, CpG73 Methylation, and HPV Status in Head and Neck Squamous Cell Carcinoma

Epigenetic studies on the role of DNA-modifying enzymes in HNSCC tumorigenesis have focused on a single enzyme or a group of enzymes. To acquire a more comprehensive insight into the expression profile of methyltransferases and demethylases, in the present study, we examined the mRNA expression of the DNA methyltransferases DNMT1, DNMT3A, and DNMT3B, the DNA demethylases TET1, TET2, TET3, and TDG, and the RNA methyltransferase TRDMT1 by RT-qPCR in paired tumor-normal tissue samples from HNSCC patients. We characterized their expression patterns in relation to regional lymph node metastasis, invasion, HPV16 infection, and CpG73 methylation. Here, we show that tumors with regional lymph node metastases (pN+) exhibited decreased expression of DNMT1, 3A and 3B, and TET1 and 3 compared to non-metastatic tumors (pN0), suggesting that metastasis requires a distinct expression profile of DNA methyltransferases/demethylases in solid tumors. Furthermore, we identified the effect of perivascular invasion and HPV16 on DNMT3B expression in HNSCC. Finally, the expression of TET2 and TDG was inversely correlated with the hypermethylation of CpG73, which has previously been associated with poorer survival in HNSCC. Our study further confirms the importance of DNA methyltransferases and demethylases as potential prognostic biomarkers as well as molecular therapeutic targets for HNSCC.

Identification of a novel glycolysis-related prognosis risk signature in triple-negative breast cancer

Introduction

Triple-negative breast cancer (TNBC) is a particularly aggressive cluster of breast cancer characterized by significant molecular heterogeneity. Glycolysis is a metabolic pathway that is significantly associated with cancer progression, metastasis, recurrence and chemoresistance. However, the potential roles of glycolysis-related genes in TNBC remain unclear.

Methods

In the present study, we identified 108 glycolysis-related differentially expressed genes (DEGs) between breast cancer (BRCA) tumor tissues and normal tissues, and we divided patients into two different clusters with significantly distinct molecular characteristics, clinicopathological features, prognosis, immune cell infiltration and mutation burden. We then constructed a 10-gene signature that classified all TNBCs into low- and high-risk groups.

Results

The high-risk group had significantly lower survival than the low-risk group, which implied that the risk score was an independent prognostic indicator for TNBC patients. Consequently, we constructed and validated a prognostic nomogram, which accurately predicted individual overall survival (OS) of TNBC. Moreover, the risk score predicted the drug sensitivity of chemotherapeutic agents and immunotherapy for TNBC patients.

Discussion

The present comprehensive analysis of glycolysis-related DEGs in TNBC provides new methods for prognosis prediction and more effective treatment strategies.

Multi-omics integration analysis of GPCRs in pan-cancer to uncover inter-omics relationships and potential driver genes

G protein-coupled receptors (GPCRs) are the largest drug target family. Unfortunately, applications of GPCRs in cancer therapy are scarce due to very limited knowledge regarding their correlations with cancers. Multi-omics data enables systematic investigations of GPCRs, yet their effective integration remains a challenge due to the complexity of the data. Here, we adopt two types of integration strategies, multi-staged and meta-dimensional approaches, to fully characterize somatic mutations, somatic copy number alterations (SCNAs), DNA methylations, and mRNA expressions of GPCRs in 33 cancers. Results from the multi-staged integration reveal that GPCR mutations cannot well predict expression dysregulation. The correlations between expressions and SCNAs are primarily positive, while correlations of the methylations with expressions and SCNAs are bimodal with negative correlations predominating. Based on these correlations, 32 and 144 potential cancer-related GPCRs driven by aberrant SCNA and methylation are identified, respectively. In addition, the meta-dimensional integration analysis is carried out by using deep learning models, which predict more than one hundred GPCRs as potential oncogenes. When comparing results between the two integration strategies, 165 cancer-related GPCRs are common in both, suggesting that they should be prioritized in future studies. However, 172 GPCRs emerge in only one, indicating that the two integration strategies should be considered concurrently to complement the information missed by the other such that obtain a more comprehensive understanding. Finally, correlation analysis further reveals that GPCRs, in particular for the class A and adhesion receptors, are generally immune-related. In a whole, the work is for the first time to reveal the associations between different omics layers and highlight the necessity of combing the two strategies in identifying cancer-related GPCRs.

CALML5 is a novel diagnostic marker for differentiating thymic squamous cell carcinoma from type B3 thymoma

BACKGROUND

Thymic squamous cell carcinoma and type B3 thymoma are primary neoplasms of the anterior mediastinum that are sometimes difficult to differentiate from one another histologically. However, only a few immunohistochemical markers are available for the differential diagnosis. The purpose of this study was to discover a novel marker for differentiating between thymic squamous cell carcinoma and type B3 thymoma.

METHODS

We used histological samples of thymic carcinomas (n = 26) and type B3 thymomas (n = 38) which were resected between 1986 and 2017. To search for candidates of differential markers, gene expression levels were evaluated in samples using promoter analysis by cap analysis of gene expression (CAGE) sequencing.

RESULTS

Promoter level expression of CALML5 genes was significantly higher in thymic carcinomas than in type B3 thymomas. We further validated the results of the CAGE analysis in all 26 thymic carcinomas and 38 type B3 thymomas by immunohistochemistry (IHC). CALML5 was strongly expressed in the cytoplasm in 19 of 26 cases with thymic carcinoma, whereas positivity at the protein level was shown in two of 38 type B3 thymomas. Thus, the sensitivity (73.1%) and specificity (94.7%) of CALML5 as markers for immunohistochemical diagnosis of thymic carcinoma were extremely high.

CONCLUSION

We identified CALML5 as a potential marker for differentiating thymic squamous cell carcinoma from type B3 thymoma. It is assumed that future clinical use of CALML5 may improve the diagnostic accuracy of differentiating between these two diseases.

DNA Methylation as a Diagnostic, Prognostic, and Predictive Biomarker in Head and Neck Cancer

Head and neck squamous cell carcinoma (HNSCC) is a term collectively used to describe all cancers that develop in the oral and nasal cavities, the paranasal sinuses, the salivary glands, the pharynx, and the larynx. The majority (75%) of all newly diagnosed cases are observed in patients with locally advanced and aggressive disease, associated with significant relapse rates (30%) and poor prognostic outcomes, despite advances in multimodal treatment. Consequently, there is an unmet need for the identification and application of tools that would enable diagnosis at the earliest possible stage, accurately predict prognostic outcomes, contribute to the timely detection of relapses, and aid in the decision for therapy selection. Recent evidence suggests that DNA methylation can alter the expression of genes in a way that it favors tumorigenesis and tumor progression in HNSCC, and therefore represents a potential source for biomarker identification. This study summarizes the current knowledge on how abnormally methylated DNA profiles in HNSCC patients may contribute to the pathogenesis of HNSCC and designate the methylation patterns that have the potential to constitute clinically valuable biomarkers for achieving significant advances in the management of the disease and for improving survival outcomes in these patients.

Using Single-Cell RNA Sequencing and MicroRNA Targeting Data to Improve Colorectal Cancer Survival Prediction

Colorectal cancer has proven to be difficult to treat as it is the second leading cause of cancer death for both men and women worldwide. Recent work has shown the importance of microRNA (miRNA) in the progression and metastasis of colorectal cancer. Here, we develop a metric based on miRNA-gene target interactions, previously validated to be associated with colorectal cancer. We use this metric with a regularized Cox model to produce a small set of top-performing genes related to colon cancer. We show that using the miRNA metric and a Cox model led to a meaningful improvement in colon cancer survival prediction and correct patient risk stratification. We show that our approach outperforms existing methods and that the top genes identified by our process are implicated in NOTCH3 signaling and general metabolism pathways, which are essential to colon cancer progression.

Dalpiciclib partially abrogates ER signaling activation induced by pyrotinib in HER2+HR+ breast cancer

Recent evidences from clinical trials (NCT04486911) revealed that the combination of pyrotinib, letrozole, and dalpiciclib exerted optimistic therapeutic effect in treating HER2⁺HR⁺ breast cancer; however, the underlying molecular mechanism remained elusive. Through the drug sensitivity test, the drug combination efficacy of pyrotinib, tamoxifen, and dalpiciclib to BT474 cells was tested. The underlying molecular mechanisms were investigated using immunofluorescence, Western blot analysis, immunohistochemical staining, and cell cycle analysis. Potential risk factor that may indicate the responsiveness to drug treatment in HER2⁺/HR⁺ breast cancer was identified using RNA-sequence and evaluated using immunohistochemical staining and in vivo drug susceptibility test. We found that pyrotinib combined with dalpiciclib exerted better cytotoxic efficacy than pyrotinib combined with tamoxifen in BT474 cells. Degradation of HER2 could enhance ER nuclear transportation, activating ER signaling pathway in BT474 cells, whereas dalpiciclib could partially abrogate this process. This may be the underlying mechanism by which combination of pyrotinib, tamoxifen, and dalpiciclib exerted best cytotoxic effect. Furthermore, CALML5 was revealed to be a risk factor in the treatment of HER2⁺/HR⁺ breast cancer and the usage of dalpiciclib might overcome the drug resistance to pyrotinib + tamoxifen due to CALML5 expression. Our study provided evidence that the usage of dalpiciclib in the treatment of HER2⁺/HR⁺ breast cancer could partially abrogate the estrogen signaling pathway activation caused by anti-HER2 therapy and revealed that CALML5 could serve as a risk factor in the treatment of HER2⁺/HR⁺ breast cancer.

Fluid Biomarkers in HPV and Non-HPV Related Oropharyngeal Carcinomas: From Diagnosis and Monitoring to Prognostication-A Systematic Review

Biomarkers are crucial in oncology, from detection and monitoring to guiding management and predicting treatment outcomes. Histological assessment of tissue biopsies is currently the gold standard for oropharyngeal cancers, but is technically demanding, invasive, and expensive. This systematic review aims to review current markers that are detectable in biofluids, which offer promising non-invasive alternatives in oropharyngeal carcinomas (OPCs). A total of 174 clinical trials from the PubMed search engine in the last 5 years were identified and screened by 4 independent reviewers. From these, 38 eligible clinical trials were found and subsequently reviewed. The biomarkers involved, categorized by human papillomavirus (HPV)-status, were further divided according to molecular and cellular levels. Recent trials investigating biomarkers for both HPV-positive and HPV-negative OPCs have approaches from various levels and different biofluids including plasma, oropharyngeal swabs, and oral rinse. Promising candidates have been found to aid in detection, staging, and predicting prognosis, in addition to well-established factors including HPV-status, drinking and smoking status. These studies also emphasize the possibility of enhancing prediction results and increasing statistical significance by multivariate analyses. Liquid biopsies offer promising assistance in enhancing personalized medicine for cancer treatment, from lowering barriers towards early screening, to facilitating de-escalation of treatment. However, further research is needed, and the combination of liquid biopsies with pre-existing methods, including in vivo imaging and invasive techniques such as neck dissections, could also be explored in future trials.

Clinical applications of liquid biopsy in HPV‐negative and HPV‐positive head and neck squamous cell carcinoma: advances and challenges

Head and neck squamous cell carcinomas (HNSCCs) represent the most common epithelial tumors that arise from mucosa of the oral cavity, pharynx, and larynx. The development of HNSCCs is usually associated with tobacco use, alcohol consumption, and human papillomavirus (HPV) infection. Most HNSCCs are diagnosed in advanced states, leading to a worse clinical outcome. Screening tests based on potential biomarkers associated with HNSCCs could improve this scenario. Liquid biopsy has emerged as a promising area of cancer investigation, offering a minimally invasive approach to track circulating biomarkers in body fluids that could potentially contribute to the diagnosis, predict prognosis, and monitor response to treatment. This review will discuss translational studies describing the clinical applications of liquid biopsy in HPV-negative and HPV-positive HNSCCs focused on circulating nucleic acids [cell-free DNA (cfDNA) and cell-free RNA (cfRNA)], circulating tumor cells (CTCs), and extracellular vesicles (EVs), which can be found in plasma, serum, and saliva.

Epigenetic and Transcriptomic Regulation Landscape in HPV+ Cancers: Biological and Clinical Implications

Human Papillomavirus (HPV) is an oncogenic virus that causes the highest number of viral-associated cancer cases and deaths worldwide, with more than 690,000 new cases per year and 342,000 deaths only for cervical cancer (CC). Although the incidence and mortality rates for CC are declining in countries where screening and vaccination programs have been implemented, other types of cancer in which HPV is involved, such as oropharyngeal cancer, are increasing, particularly in men. Mutational and transcriptional profiles of various HPV-associated neoplasms have been described, and accumulated evidence has shown the oncogenic capacity of E6, E7, and E5 genes of high-risk HPV. Interestingly, transcriptomic analysis has revealed that although a vast majority of the human genome is transcribed into RNAs, only 2% of transcripts are translated into proteins. The remaining transcripts lacking protein-coding potential are called non-coding RNAs. In addition to the transfer and ribosomal RNAs, there are regulatory non-coding RNAs classified according to size and structure in long non-coding RNAs (lncRNAs), circular RNAs (circRNAs), and small RNAs; such as microRNAs (miRNAs), piwi-associated RNAs (piRNAs), small nucleolar RNAs (snoRNAs) and endogenous short-interfering RNAs. Recent evidence has shown that lncRNAs, miRNAs, and circRNAs are aberrantly expressed under pathological conditions such as cancer. In addition, those transcripts are dysregulated in HPV-related neoplasms, and their expression correlates with tumor progression, metastasis, poor prognosis, and recurrence. Nuclear lncRNAs are epigenetic regulators involved in controlling gene expression at the transcriptional level through chromatin modification and remodeling. Moreover, disruption of the expression profiles of those lncRNAs affects multiple biological processes such as cell proliferation, apoptosis, and migration. This review highlights the epigenetic alterations induced by HPV, from infection to neoplastic transformation. We condense the epigenetic role of non-coding RNA alterations and their potential as biomarkers in transformation’s early stages and clinical applications. We also summarize the molecular mechanisms of action of nuclear lncRNAs to understand better their role in the epigenetic control of gene expression and how they can drive the malignant phenotype of HPV-related neoplasia. Finally, we review several chemical and epigenetic therapy options to prevent and treat HPV-associated neoplasms.

Epigenetic editing and tumor‑dependent immunosuppressive signaling in head and neck malignancies (Review)

Head and neck cancer (HNC) comprises a heterogeneous variety of malignant tumors, characterized by a relatively high tumor mutation burden. Previous data have revealed that immune system dysfunction appears to serve a key role in the development and progression of HNC and established immunosuppression is vital for evading the host immune response. Despite progress in chemotherapy and radiotherapy, the survival rate of patients with HNC is still low. Therefore, the present review discusses the development of novel immunotherapy approaches based on the various immune cell signaling routes that trigger drug resistance and immunosuppression. Additionally, the present review discusses the epigenetic alterations, including DNA methylation, histone modifications, chromatin remodeling and non-coding RNAs that drive and support HNC progression. Furthermore, the role of cancer-associated fibroblasts, tumor macrophages and myeloid cells in tumor-related immunosuppression are considered. Specifically, the molecular immune-related mechanisms in the tumor microenvironment, which lead to decreased drug sensitivity and tumor relapse, and strategies for reversing drug resistance and targeting immunosuppressive tumor networks are discussed. Deciphering these molecular mechanisms is essential for preclinical and clinical investigations in order to enhance therapeutic efficacy. Furthermore, an improved understanding of these immune cell signaling pathways that drive immune surveillance, immune-driven inflammation and tumor-related immunosuppression is necessary for future personalized HNC-based therapeutic approaches.

Characterization of Epigenomic Alterations in HPV16+ Head and Neck Squamous Cell Carcinomas

BACKGROUND

Epigenetic changes associated with human papillomavirus (HPV)-driven tumors have been described; however, HPV type-specific alterations are less well understood. We sought to compare HPV16-specific methylation changes with those in virus-unassociated head and neck squamous cell carcinomas (HNSCC).

METHODS

Within The Cancer Genome Atlas, 59 HPV16+ HNSCC, 238 nonviral HNSCC (no detectable HPV or other viruses), and 50 normal head and neck tissues were evaluated. Significant differentially methylated regions (DMR) were selected, and key associated genes were identified. Partial least squares models were generated to predict HPV16 status in additional independent samples.

RESULTS

HPV infection in HNSCC is associated with type-specific methylomic profiles. Multiple significant DMRs were identified between HPV16+, nonviral, and normal samples. The most significant differentially methylated genes, SYCP2, MSX2, HLTF, PITX2, and GRAMD4, demonstrated HPV16-associated methylation patterns with corresponding alterations in gene expression. Phylogenetically related HPV types (alpha-9 species; HPV31, HPV33, and HPV35) demonstrated a similar methylation profile to that of HPV16 but differed from those seen in other types, such as HPV18 and 45 (alpha-7).

CONCLUSIONS

HNSCC linked to HPV16 and types from the same alpha species are associated with a distinct methylation profile. This HPV16-associated methylation pattern is also detected in cervical cancer and testicular germ cell tumors. We present insights into both shared and unique methylation alterations associated with HPV16+ tumors and may have implications for understanding the clinical behavior of HPV-associated HNSCC.

IMPACT

HPV type-specific methylomic changes may contribute to understanding biologic mechanisms underlying differences in clinical behavior among different HPV+ and HPV- HNSCC.

Combined multimodal ctDNA analysis and radiological imaging for tumor surveillance in Non-small cell lung cancer

BACKGROUND

Radiology is the current standard for monitoring treatment responses in lung cancer. Limited sensitivity, exposure to ionizing radiations and related sequelae constitute some of its major limitation. Non-invasive and highly sensitive methods for early detection of treatment failures and resistance-associated disease progression would have additional clinical utility.

METHODS

We analyzed serially collected plasma and paired tumor samples from lung cancer patients (61 with stage IV, 48 with stages I-III disease) and 61 healthy samples by means of next-generation sequencing, radiological imaging and droplet digital polymerase chain reaction (ddPCR) mutation and methylation assays.

RESULTS

A 62% variant concordance between tumor-reported and circulating-free DNA (cfDNA) sequencing was observed between baseline liquid and tissue biopsies in stage IV patients. Interestingly, ctDNA sequencing allowed for the identification of resistance-mediating p.T790M mutations in baseline plasma samples for which no such mutation was observed in the corresponding tissue. Serial circulating tumor DNA (ctDNA) mutation analysis by means of ddPCR revealed a general decrease in ctDNA loads between baseline and first reassessment. Additionally, serial ctDNA analyses only recapitulated computed tomography (CT) -monitored tumor dynamics of some, but not all lesions within the same patient. To complement ctDNA variant analysis we devised a ctDNA methylation assay (^methcfDNA) based on methylation-sensitive restriction enzymes. cfDNA methylation showed and area under the curve (AUC) of > 0.90 in early and late stage cases. A decrease in ^methcfDNA between baseline and first reassessment was reflected by a decrease in CT-derive tumor surface area, irrespective of tumor mutational status.

CONCLUSION

Taken together, our data support the use of cfDNA sequencing for unbiased characterization of the molecular tumor architecture, highlights the impact of tumor architectural heterogeneity on ctDNA-based tumor surveillance and the added value of complementary approaches such as cfDNA methylation for early detection and monitoring.

Transcriptomic-Based Identification of the Immuno-Oncogenic Signature of Cholangiocarcinoma for HLC-018 Multi-Target Therapy Exploration

Cholangiocarcinomas (CHOLs), hepatobiliary malignancies, are characterized by high genetic heterogeneity, a rich tumor microenvironment, therapeutic resistance, difficulty diagnosing, and poor prognoses. Current knowledge of genetic alterations and known molecular markers for CHOL is insufficient, necessitating the need for further evaluation of the genome and RNA expression data in order to identify potential therapeutic targets, clarify the roles of these targets in the tumor microenvironment, and explore novel therapeutic drugs against the identified targets. Consequently, in our attempt to explore novel genetic markers associated with the carcinogenesis of CHOL, five genes (SNX15, ATP2A1, PDCD10, BET1, and HMGA2), collectively termed CHOL-hub genes, were identified via integration of differentially expressed genes (DEGs) from relatively large numbers of samples from CHOL GEO datasets. We further explored the biological functions of the CHOL-hub genes and found significant enrichment in several biological process and pathways associated with stem cell angiogenesis, cell proliferation, and cancer development, while the interaction network revealed high genetic interactions with a number of onco-functional genes. In addition, we established associations between the CHOL-hub genes and tumor progression, metastasis, tumor immune and immunosuppressive cell infiltration, dysfunctional T-cell phenotypes, poor prognoses, and therapeutic resistance in CHOL. Thus, we proposed that targeting CHOL-hub genes could be an ideal therapeutic approach for treating CHOLs, and we explored the potential of HLC-018, a novel benzamide-linked small molecule, using molecular docking of ligand-receptor interactions. To our delight, HLC-018 was well accommodated with high binding affinities to binding pockets of CHOL-hub genes; more importantly, we found specific interactions of HLC-018 with the conserved sequence of the AT-hook DNA-binding motif of HMGA2. Altogether, our study provides insights into the immune-oncogenic phenotypes of CHOL and provides valuable information for our ongoing experimental validation.

A Transcriptomic Analysis of Head and Neck Squamous Cell Carcinomas for Prognostic Indications

Survival analysis of the Cancer Genome Atlas (TCGA) dataset is a well-known method for discovering gene expression-based prognostic biomarkers of head and neck squamous cell carcinoma (HNSCC). A cutoff point is usually used in survival analysis for patient dichotomization when using continuous gene expression values. There is some optimization software for cutoff determination. However, the software's predetermined cutoffs are usually set at the medians or quantiles of gene expression values. There are also few clinicopathological features available in pre-processed datasets. We applied an in-house workflow, including data retrieving and pre-processing, feature selection, sliding-window cutoff selection, Kaplan-Meier survival analysis, and Cox proportional hazard modeling for biomarker discovery. In our approach for the TCGA HNSCC cohort, we scanned human protein-coding genes to find optimal cutoff values. After adjustments with confounders, clinical tumor stage and surgical margin involvement were found to be independent risk factors for prognosis. According to the results tables that show hazard ratios with Bonferroni-adjusted p values under the optimal cutoff, three biomarker candidates, CAMK2N1, CALML5, and FCGBP, are significantly associated with overall survival. We validated this discovery by using the another independent HNSCC dataset (GSE65858). Thus, we suggest that transcriptomic analysis could help with biomarker discovery. Moreover, the robustness of the biomarkers we identified should be ensured through several additional tests with independent datasets.

Application of Data Science in Circulating Tumor DNA Detection: A Promising Avenue Towards Liquid Biopsy

The circulating tumor DNA (ctDNA), as a promising biomarker of liquid biopsy, has potential clinical relevance on the molecular diagnosis and monitoring of cancer. However, the trace concentration level of ctDNA in the peripheral blood restricts its extensive clinical application. Recently, high-throughput-based methodologies have been leveraged to improve the sensitivity and specificity of ctDNA detection, showing a promising avenue towards liquid biopsy. This review briefly summarizes the high-throughput data features concerned by current ctDNA detection strategies and the technical obstacles, potential solutions, and clinical relevance of current ctDNA profiling technologies. We also highlight future directions improving the limit of detection of ctDNA for better clinical application. This review may serve as a reference for the crosslinks between data science and ctDNA-based liquid biopsy, benefiting clinical translation in advanced cancer diagnosis.

Salivary DNA Methylation as an Epigenetic Biomarker for Head and Neck Cancer. Part I: A Diagnostic Accuracy Meta-Analysis

DNA hypermethylation is an important epigenetic mechanism for gene expression inactivation in head and neck cancer (HNC). Saliva has emerged as a novel liquid biopsy representing a potential source of biomarkers. We performed a comprehensive meta-analysis to evaluate the overall diagnostic accuracy of salivary DNA methylation for detecting HNC. PubMed EMBASE, Web of Science, LILACS, and the Cochrane Library were searched. Study quality was assessed by the Quality Assessment for Studies of Diagnostic Accuracy-2, and sensitivity, specificity, positive likelihood ratio (PLR), negative likelihood ratio (NLR), diagnostic odds ratio (dOR), and their corresponding 95% confidence intervals (CIs) were calculated using a bivariate random-effect meta-analysis model. Meta-regression and subgroup analyses were performed to assess heterogeneity. Eighty-four study units from 18 articles with 8368 subjects were included. The pooled sensitivity and specificity of salivary DNA methylation were 0.39 and 0.87, respectively, while PLR and NLR were 3.68 and 0.63, respectively. The overall area under the curve (AUC) was 0.81 and the dOR was 8.34. The combination of methylated genes showed higher diagnostic accuracy (AUC, 0.92 and dOR, 36.97) than individual gene analysis (AUC, 0.77 and dOR, 6.02). These findings provide evidence regarding the potential clinical application of salivary DNA methylation for HNC diagnosis.

Identification of calmodulin-like protein 5 as tumor-suppressor gene silenced during early stage of carcinogenesis in squamous cell carcinoma of uterine cervix

In the course of identifying the molecular mechanism that is related to strong cell-cell adhesion in stratified structures of the squamous epithelium, calmodulin-like protein 5 (CALML5) was identified as a spinous structure-associated protein by producing monoclonal antibodies with the use of the crude intercellular portion of squamous tissue as an immunogen and by subsequent morphologic screening. By electrophoretic mobility shift assay (EMSA) and a series of mutagenesis studies, two transcription factors, ZNF750 and KLF4, by binding in line to the CALML5 gene promoter, were found to play a central role in CALML5 transcription. Knockdown of CALML5 by siRNA in the A431 cell line that expresses high levels of CALML5 resulted in the acceleration of wound confluence in a scratch assay, indicating that CALML5 functions as a tumor-suppressor in uterine cervical cancer. Immunohistochemical evaluation of squamous intraepithelial lesions, carcinoma in situ (CIS) and invasive uterine cancer, revealed a reduction in CALML5 expression during the stages of CIS through various molecular pathways including the blockage of the nuclear translocation of KLF4. Conversely, restoration of the nuclear translocation of KLF4 by inhibiting ERK-signaling reactivated CALML5 expression in ME180 cells expressing low levels of CALML5. Thus, alteration of the p63-ZNF750-KLF4 axis may result in critical functional loss of CALM-related genes during cancer progression. Although the morphological association of CALML5 with the spiny-structure in relation to cell motility is not clear, evaluation of CALML5 expression provides a useful diagnostic indicator of differentiating dysplasia, preinvasive and invasive cervical cancers.

Validation of Selected Head and Neck Cancer Prognostic Markers from the Pathology Atlas in an Oral Tongue Cancer Cohort

The Pathology Atlas is an open-access database that reports the prognostic value of protein-coding transcripts in 17 cancers, including head and neck cancer. However, cancers of the various head and neck anatomical sites are specific biological entities. Thus, the aim of the present study was to validate promising prognostic markers for head and neck cancer reported in the Pathology Atlas in oral tongue squamous cell carcinoma (OTSCC). We selected three promising markers from the Pathology Atlas (CALML5, CD59, LIMA1), and analyzed their prognostic value in a Norwegian OTSCC cohort comprising 121 patients. We correlated target protein and mRNA expression in formalin-fixed, paraffin-embedded cancer tissue to five-year disease-specific survival (DSS) in univariate and multivariate analyses. Protein expression of CALML5 and LIMA1 were significantly associated with five-year DSS in the OTSCC cohort in univariate analyses (p = 0.016 and p = 0.043, respectively). In multivariate analyses, lymph node metastases, tumor differentiation, and CALML5 were independent prognosticators. The prognostic role of the other selected markers for head and neck cancer patients identified through unbiased approaches could not be validated in our OTSCC cohort. This underlines the need for subsite-specific analyses for head and neck cancer.

The MAL Protein, an Integral Component of Specialized Membranes, in Normal Cells and Cancer

The MAL gene encodes a 17-kDa protein containing four putative transmembrane segments whose expression is restricted to human T cells, polarized epithelial cells and myelin-forming cells. The MAL protein has two unusual biochemical features. First, it has lipid-like properties that qualify it as a member of the group of proteolipid proteins. Second, it partitions selectively into detergent-insoluble membranes, which are known to be enriched in condensed cell membranes, consistent with MAL being distributed in highly ordered membranes in the cell. Since its original description more than thirty years ago, a large body of evidence has accumulated supporting a role of MAL in specialized membranes in all the cell types in which it is expressed. Here, we review the structure, expression and biochemical characteristics of MAL, and discuss the association of MAL with raft membranes and the function of MAL in polarized epithelial cells, T lymphocytes, and myelin-forming cells. The evidence that MAL is a putative receptor of the epsilon toxin of Clostridium perfringens, the expression of MAL in lymphomas, the hypermethylation of the MAL gene and subsequent loss of MAL expression in carcinomas are also presented. We propose a model of MAL as the organizer of specialized condensed membranes to make them functional, discuss the role of MAL as a tumor suppressor in carcinomas, consider its potential use as a cancer biomarker, and summarize the directions for future research.

DNA Methylation and HPV-Associated Head and Neck Cancer

Head and neck squamous cell carcinoma (HNSCC), especially oropharyngeal squamous cell carcinoma (OPSCC), has recently been found to be significantly associated with human papillomavirus (HPV) infection. The incidence of OPSCC has been increasing and surpassed the number of cervical cancer cases in the United States. Although HPV-associated OPSCC has a relatively better prognosis than HPV-negative cancer, approximately 20% of HPV-associated HNSCC patients show a poor prognosis or therapeutic response, and the molecular mechanism behind this outcome in the intermediate-risk group is yet to be elucidated. These biological differences between HPV-associated HNSCC and HPV-negative HNSCC are partly explained by the differences in mutation patterns. However, recent reports have revealed that epigenetic dysregulation, such as dysregulated DNA methylation, is a strikingly common pathological feature of human malignancy. Notably, viral infections can induce aberrant DNA methylation, leading to carcinogenesis, and HPV-associated HNSCC cases tend to harbor a higher amount of aberrantly methylated DNA than HPV-negative HNSCC cases. Furthermore, recent comprehensive genome-wide DNA-methylation analyses with large cohorts have revealed that a sub-group of HPV-associated HNSCC correlates with increased DNA methylation. Accordingly, in this review, we provide an overview of the relationship between DNA methylation and HPV-associated HNSCC.

Screening of gene markers related to the prognosis of metastatic skin cutaneous melanoma based on Logit regression and survival analysis

BACKGROUND

Bioinformatics was used to analyze the skin cutaneous melanoma (SKCM) gene expression profile to provide a theoretical basis for further studying the mechanism underlying metastatic SKCM and the clinical prognosis.

METHODS

We downloaded the gene expression profiles of 358 metastatic and 102 primary (nonmetastatic) CM samples from The Cancer Genome Atlas (TCGA) database as a training dataset and the GSE65904 dataset from the National Center for Biotechnology Information database as a validation dataset. Differentially expressed genes (DEGs) were screened using the limma package of R3.4.1, and prognosis-related feature DEGs were screened using Logit regression (LR) and survival analyses. We also used the STRING online database, Cytoscape software, and Database for Annotation, Visualization and Integrated Discovery software for protein-protein interaction network, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses based on the screened DEGs.

RESULTS

Of the 876 DEGs selected, 11 (ZNF750, NLRP6, TGM3, KRTDAP, CAMSAP3, KRT6C, CALML5, SPRR2E, CD3G, RTP5, and FAM83C) were screened using LR analysis. The survival prognosis of nonmetastatic group was better compared to the metastatic group between the TCGA training and validation datasets. The 11 DEGs were involved in 9 KEGG signaling pathways, and of these 11 DEGs, CALML5 was a feature DEG involved in the melanogenesis pathway, 12 targets of which were collected.

CONCLUSION

The feature DEGs screened, such as CALML5, are related to the prognosis of metastatic CM according to LR. Our results provide new ideas for exploring the molecular mechanism underlying CM metastasis and finding new diagnostic prognostic markers.

Telomere shortening in head and neck cancer: association between DNA demethylation and survival

A growing body of evidence indicates that telomere dysfunction is a biological marker of progression in several types of cancer. However, the association between head and neck squamous cell carcinoma (HNSCC) and telomere length (TL) remains unknown. We measured the absolute TL levels in a well-characterised dataset of 211 tumoral vs normal tissues obtained from the same patients by quantitative polymerase chain reaction assay. Normalised TL levels were significantly lower in tumour samples than in normal tissue (P < 0.001) and there was a positive correlation between tumour tissue and normal mucosal tissue (R² = 0.176, P < 0.001). We were able to distinguish two classes, one with a tumour/normal TL ratio ≤ 0.3 (38.4%), which showed clear telomere erosion, and the other with a tumour/normal TL ratio > 0.3 (61.6%), in which the TL was slightly shorter or longer than that in normal tissue. Notably, the tumour/normal TL ratio was correlated with the likelihood of disease recurrence (P = 0.002), the 5-hydroxymethylcytosine level (P = 0.043), and expression of the ten-eleven translocation (TET) gene (P = 0.043). Our findings show that TL shortening and subsequent low levels of 5-hydroxymethylcytosine and TET expression may contribute to development of HNSCC.

Surveillance and Monitoring Techniques for HPV-Related Head and Neck Squamous Cell Carcinoma: Circulating Tumor DNA

OPINION STATEMENT

Human papilloma virus (HPV) related head and neck cancer is rising in prevalence, preferentially affecting young patients and imparting long term toxicities. Despite this, there are no screening tests or clinical biomarkers for treatment monitoring. HPV circulating tumor DNA (HPV ctDNA) represents a novel circulating biomarker which may provide real-time assessment of tumor response to therapy and recurrence. Early work suggests the promise of this assay as a predictive biomarker in numerous clinical settings, namely risk of recurrence after chemoradiation in locally advanced disease. Advancement of these findings to the clinic will require a collaborative effort in the field, including technical harmonization of assay testing characteristics, understanding of the normal kinetics in patients being treated with standard of care therapies, and appropriately designed phase III trials prior to implementation in the clinic. If successful, HPV ctDNA has the potential to revolutionize clinical trial treatment paradigms and transform patient care.

Head and Neck Squamous Cell Carcinoma: Epigenetic Landscape

Head and neck squamous carcinoma (HNSCC) constitutes the sixth most prevalent cancer worldwide. The molecular pathogenesis of HNSCC includes disorders in cell cycle, intercellular signaling, proliferation, squamous cell differentiation and apoptosis. In addition to the genetic mutations, changes in HNSCC are also characterized by the accumulation of epigenetic alterations such as DNA methylation, histone modifications, non-coding RNA activity and RNA methylation. In fact, some of them may promote cancer formation and progression by controlling the gene expression machinery, hence, they could be used as biomarkers in the clinical surveillance of HNSCC or as targets for therapeutic strategies. In this review, we focus on the current knowledge regarding epigenetic modifications observed in HNSCC and its predictive value for cancer development.

Long interspersed nuclear element 1 hypomethylation has novel prognostic value and potential utility in liquid biopsy for oral cavity cancer

Background

New biomarkers are urgently needed to improve personalized treatment approaches for head and neck squamous cell carcinoma (HNSCC). Global DNA hypomethylation has wide-ranging functions in multistep carcinogenesis, and the hypomethylation of long interspersed nucleotide element-1 (LINE-1) is related to increased retrotransposon activity and induced genome instability. However, little information is available regarding LINE-1 hypomethylation and its prognostic implications in HNSCC.

Methods

In this study, we analyzed LINE-1 hypomethylation levels in a well-characterized dataset of 317 primary HNSCC tissues and 225 matched pairs of normal mucosa tissues, along with five oral cavity cancer (OCC) circulating tumor DNA (ctDNA) samples using quantitative real-time methylation and unmethylation PCR. The analysis was performed according to various clinical characteristics and prognostic implications.

Results

The results demonstrated that LINE-1 hypomethylation levels were significantly higher in the HNSCC tissues than in corresponding normal tissues from the same individuals (P < 0.001). Univariate analysis revealed that high levels of LINE-1 hypomethylation were correlated with poor disease-free survival (DFS; log-rank test, P = 0.038), whereas multivariate analysis demonstrated that they were significant independent prognostic factor for DFS (hazard ratio: 2.10, 95% confidence interval: 1.02–4.36; P = 0.045). Moreover, samples with high LINE-1 hypomethylation levels exhibited the greatest decrease in 5-hydroxymethylcytosine (5-hmC) levels and increase in tumor-suppressor gene methylation index (P = 0.006 and P < 0.001, respectively). Further, ctDNA studies also showed that LINE-1 hypomethylation had high predictive ability in OCC.

Conclusions

LINE-1 hypomethylation is associated with a higher risk of early OCC relapse, and is hence, a potential predictive biomarker for OCC. Furthermore, 5-hmC levels also exhibited predictive potential in OCC, based on their inverse correlation with LINE-1 hypomethylation levels. LINE-1 hypomethylation analysis, therefore, has applications in determining patient prognosis and real-time surveillance of disease recurrence, and could serve as an alternative method for OCC screening.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s40364-020-00235-y.

Exploration of Feasible Immune Biomarkers for Immune Checkpoint Inhibitors in Head and Neck Squamous Cell Carcinoma Treatment in Real World Clinical Practice

Recurrent locally advanced or metastatic head and neck squamous cell carcinoma (HNSCC) is associated with dismal prognosis because of its highly invasive behavior and resistance to conventional intensive chemotherapy. The combination of targeted therapy and conventional chemotherapy has significantly improved clinical outcomes. In recent years, the development of immunotherapies, such as immune checkpoint inhibitors (ICIs), has further increased treatment responses and prolonged survival. However, the limited response rate, risk of immunotherapy-related adverse effects and high cost of immunotherapy make the identification of predictive markers to optimize treatment efficacy a critical issue. Biomarkers are biological molecules that have been widely utilized to predict treatment response to certain treatments and clinical outcomes or to detect disease. An ideal biomarker should exhibit good predictive ability, which can guide healthcare professionals to achieve optimal treatment goals and bring clinical benefit to patients. In this review, we summarized the results of recent and important studies focused on HNSCC ICI immunotherapy and discussed potential biomarkers including their strengths and limitations, aiming to gain more insight into HNSCC immunotherapy in real world clinical practice.

G Protein-Coupled Receptor Genes, PTGDR1, PTGDR2, and PTGIR, Are Candidate Epigenetic Biomarkers and Predictors for Treated Patients with HPV-Associated Oropharyngeal Cancer

Differences in the biology of human papillomavirus (HPV)-associated oropharyngeal cancers (OPCs) and HPV-negative OPCs may have implications in patient management. Early detection is imperative to reduce HPV-associated OPC mortality. Circulating tumor DNA (ctDNA) can potentially serve as a biomarker for monitoring clinically relevant cancer-related genetic and epigenetic modifications. We analyzed the methylation status of 24 G protein-coupled receptor (GPCR) genes in verification (85 OPC primary samples) and validation (8 OPC ctDNA samples) studies using quantitative methylation-specific polymerase chain reaction (Q-MSP). The Q-MSP-based verification study with 85 OPC primary samples revealed the GPCR genes that were significantly associated with recurrence in high methylation groups (≥14 methylated genes) with OPC and HPV-associated OPC (p < 0.001). In the Kaplan–Meier estimate and multivariate Cox proportional hazard analyses, 13 GPCR genes were significantly related to increased recurrence in the methylation group. Furthermore, the validation study on ctDNA showed that three of these genes (Prostaglandin D2 receptor 1: PTGDR1, Prostaglandin D2 receptor 2: PTGDR2, and Prostaglandin I2 Receptor: PTGIR) had a prediction performance as emerging biomarkers. We characterized the relationship between the methylation status of GPCR genes and outcomes in HPV-associated OPC. Our results highlight the potential utility of ctDNA methylation-based detection for the clinical management of HPV-associated OPC.