Molecular Tumor Subtypes of HPV-Positive Head and Neck Cancers: Biological Characteristics and Implications for Clinical Outcomes

Identification of TTC21A as a Potential Prognostic Marker in Head and Neck Squamous Cell Carcinoma: In Silico Analysis

BACKGROUND/AIM

Tetratricopeptide repeat domain 21A (TTC21A) plays a crucial role in ciliary function and has been associated with various pathogenic processes, including carcinogenesis. However, its role in head and neck squamous cell carcinoma (HNSCC) has not been elucidated.

MATERIALS AND METHODS

Based on the sequencing and microarray data of HNSCC from publicly available databases, the expression of TTC21A was compared between different subgroups based on clinical and molecular parameters. The survival analysis and regression analysis were conducted using the Kaplan-Meier method and the Cox method, respectively. Functional analysis was performed by the Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and gene set enrichment analysis (GSEA) tools. Immune infiltration analysis was performed based on the expression of TTC21A.

RESULTS

TTC21A decreased in tumor tissues and was associated with N stage, histologic grade, HPV infection, and TP53 mutation in HNSCC. TTC21A was an independent indicator of overall survival for patients with HNSCC. A high level of TTC21A expression indicated a favorable prognosis. The TTC21A expression level was involved with immune-related signaling regulation, immune-related gene expression, and immune cell infiltration. TTC21A expression was potent in predicting immunotherapeutic benefits.

CONCLUSION

TTC21A, as a potential predictor of favorable outcomes and immunotherapy response for HNSCC, is related to immune-related signaling regulation, immune-related gene expression, and immune cell infiltration.

Bioinformatics analysis of immune characteristics in tumors with alternative carcinogenesis pathways induced by human papillomaviruses

BACKGROUND

Human papillomaviruses (HPVs) induce a subset of head and neck squamous cell carcinomas (HNSCC) and anogenital cancers, particularly cervical cancer (CC). The major viral proteins that contribute to tumorigenesis are the E6 and E7 oncoproteins, whose expression is usually enhanced after the integration of viral DNA into the host genome. Recently, an alternative tumorigenesis pathway has been suggested in approximately half of HNSCC and CC cases associated with HPV infection. This pathway is characterized by extrachromosomal HPV persistence and increased expression of the viral E2, E4, and E5 genes. The E6, E7, E5, and E2 proteins have been shown to modify the expression of numerous cellular immune-related genes. The antitumor immune response is a critical factor in the prognosis of HPV-driven cancers, and its characterization may contribute to the prediction and personalization of the increasingly used cancer immunotherapy.

METHODS

We analyzed the immune characteristics of HPV-dependent tumors and their association with carcinogenesis types. Transcriptomic HNSCC and CC datasets from The Cancer Genome Atlas were used for this analysis.

RESULTS

Clustering with immune-related genes resulted in two clusters of HPV16-positive squamous cell carcinomas in both tumor types: cluster 1 had higher activation of immune responses, including stimulation of the antigen processing and presentation pathway, which was associated with higher immune cell infiltration and better overall survival, and cluster 2 was characterized by keratinization. In CC, the distribution of tumor samples into clusters 1 and 2 did not depend on the level of E2/E5 expression, but in HNSCC, most E2/E5-high tumors were localized in cluster 1 and E2/E5-low tumors in cluster 2. Further analysis did not reveal any association between the E2/E5 levels and the expression of immune-related genes.

CONCLUSIONS

Our results suggest that while the detection of immune responses associated with preserved expression of genes encoding components of antigen processing and presentation machinery in HPV-driven tumors may be markers of better prognosis and an important factor in therapy selection, the type of carcinogenesis does not seem to play a decisive role in the induction of antitumor immunity.

Clinical, morphologic and molecular heterogeneity of HPV-associated oropharyngeal cancer

The incidence of human papillomavirus-positive (HPV+) oropharyngeal squamous cell carcinoma (OPSCC) is rising rapidly and has exceeded cervical cancer to become the most common HPV-induced cancer in developed countries. Since patients with HPV + OPSCC respond very favorably to standard aggressive treatment, the emphasis has changed to reducing treatment intensity. However, recent multi-center clinical trials failed to show non-inferiority of de-escalation strategies on a population basis, highlighting the need to select low-risk patients likely to respond to de-intensified treatments. In contrast, there is a substantial proportion of patients who develop recurrent disease despite aggressive therapy. This supports that HPV + OPSCC is not a homogeneous disease, but comprises distinct subtypes with clinical and biological variations. The overall goal for this review is to identify biomarkers for HPV + OPSCC that may be relevant for patient stratification for personalized treatment. We discuss HPV + OPSCC as a heterogeneous disease from multifaceted perspectives including clinical behavior, tumor morphology, and molecular phenotype. Molecular profiling from bulk tumors as well as single-cell sequencing data are discussed as potential driving factors of heterogeneity between tumor subgroups. Finally, we evaluate key challenges that may impede in-depth investigations of HPV + OPSCC heterogeneity and outline potential future directions, including a section on racial and ethnic differences.

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.

Tumorigenic activation around HPV integrated sites in head and neck squamous cell carcinoma

Human papillomavirus (HPV) is causally involved in the development of head and neck squamous cell carcinoma (HNSCC). The integration of HPV drives tumorigenesis through expression of oncogenic viral genes as well as genomic alterations in surrounding regions. To elucidate involvement of epigenetic dysregulation in tumorigenesis, we here performed integrated analyses of the epigenome, transcriptome and interactome using ChIP-seq, RNA-seq and Hi-C and 4C-seq for HPV(+) HNSCCs. We analyzed clinical HNSCC using The Cancer Genome Atlas data and found that genes neighboring HPV integration sites were significantly upregulated and were correlated with oncogenic phenotypes in HPV(+) HNSCCs. While we found four HPV integration sites in HPV(+) HNSCC cell line UPCI-SCC-090 through target enrichment sequencing, 4C-seq revealed 0.5 to 40 Mb of HPV-interacting regions (HPVIRs) where host genomic regions interacted with integrated HPV genomes. While 9% of the HPVIRs were amplified and activated epigenetically forming super-enhancers, the remaining non-amplified regions were found to show a significant increase in H3K27ac levels and an upregulation of genes associated with GO terms, for example, Signaling by WNT and Cell Cycle. Among those genes, ITPR3 was significantly upregulated, involving UPCI-SCC-090-specific super-enhancer formation around the ITPR3 promoter and in the 80-kb-downstream region. The knockdown of ITPR3 by siRNA or CRISPR deletions of the distant enhancer region led to a significant suppression of cell proliferation. The epigenetic activation of HPVIRs was also confirmed in other cell lines, UM-SCC-47 and UM-SCC-104. These data indicate that epigenetic activation in HPVIRs contributes, at least partially, to genesis of HPV(+) HNSCC.

Distribution and Clinical Significance of HPV16 Variants in Head and Neck Squamous Cell Carcinomas: Data from a Portuguese Cohort and Systematic Review

INTRODUCTION

Genomic variants of the human papillomavirus type 16 (HPV16) are thought to play differential roles in the susceptibility to head and neck squamous cell carcinomas (HNSCC) and its biological behaviour. This study aimed to establish the prevalence of HPV16 variants in an HNSCC cohort and associate them with clinical pathological characteristics and patient survival.

METHODS

We retrieved samples and clinical data from 68 HNSCC patients. DNA samples were available from tumour biopsy at the time of the primary diagnosis. Targeted next-generation sequencing was used to obtain whole-genome sequences, and variants were established based on phylogenetic classification.

RESULTS

74% of samples clustered in lineage A, 5.7% in lineage B, 2.9% in lineage C, and 17.1% in lineage D. Comparative genome analysis revealed 243 single nucleotide variations. Of these, one hundred were previously reported, according to our systematic review. No significant associations with clinical pathological variables or patient survival were observed. The E6 amino acid variations E31G, L83V, and D25E and E7 N29S, associated with cervical cancer, were not observed, except for N29S in a single patient.

CONCLUSION

These results provide a comprehensive genomic map of HPV16 in HSNCC, highlighting tissue-specific characteristics which will help design tailored therapies for cancer patients.

The Intratumor Bacterial and Fungal Microbiome Is Characterized by HPV, Smoking, and Alcohol Consumption in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) tumor phenotypes and clinical outcomes are significantly influenced by etiological agents, such as HPV infection, smoking, and alcohol consumption. Accordingly, the intratumor microbiome has been increasingly implicated in cancer progression and metastasis. However, few studies characterize the intratumor microbial landscape of HNSCC with respect to these etiological agents. In this study, we aimed to investigate the bacterial and fungal landscape of HNSCC in association with HPV infection, smoking, and alcohol consumption. RNA-sequencing data were extracted from The Cancer Genome Atlas (TCGA) regarding 449 tissue samples and 44 normal samples. Pathoscope 2.0 was used to extract the microbial reads. Microbe abundance was compared to clinical variables, oncogenic signatures, and immune-associated pathways. Our results demonstrated that a similar number of dysregulated microbes was overabundant in smokers and nonsmokers, while heavy drinkers were characterized by an underabundance of dysregulated microbes. Conversely, the majority of dysregulated microbes were overabundant in HPV+ tumor samples when compared to HPV- tumor samples. Moreover, we observed that many dysregulated microbes were associated with oncogenic and metastatic pathways, suggesting their roles in influencing carcinogenesis. These microbes provide insights regarding potential mechanisms for tumor pathogenesis and progression with respect to the three etiological agents.

p63 Directs Subtype-Specific Gene Expression in HPV+ Head and Neck Squamous Cell Carcinoma

The complex heterogeneity of head and neck squamous cell carcinoma (HNSCC) reflects a diverse underlying etiology. This heterogeneity is also apparent within Human Papillomavirus-positive (HPV+) HNSCC subtypes, which have distinct gene expression profiles and patient outcomes. One aggressive HPV+ HNSCC subtype is characterized by elevated expression of genes involved in keratinization, a process regulated by the oncogenic transcription factor ΔNp63. Furthermore, the human TP63 gene locus is a frequent HPV integration site and HPV oncoproteins drive ΔNp63 expression, suggesting an unexplored functional link between ΔNp63 and HPV+ HNSCC. Here we show that HPV+ HNSCCs can be molecularly stratified according to ΔNp63 expression levels and derive a ΔNp63-associated gene signature profile for such tumors. We leveraged RNA-seq data from p63 knockdown cells and ChIP-seq data for p63 and histone marks from two ΔNp63^high HPV+ HNSCC cell lines to identify an epigenetically refined ΔNp63 cistrome. Our integrated analyses reveal crucial ΔNp63-bound super-enhancers likely to mediate HPV+ HNSCC subtype-specific gene expression that is anchored, in part, by the PI3K-mTOR pathway. These findings implicate ΔNp63 as a key regulator of essential oncogenic pathways in a subtype of HPV+ HNSCC that can be exploited as a biomarker for patient stratification and treatment choices.

Clinical Validity of a Prognostic Gene Expression Cluster-Based Model in Human Papillomavirus-Positive Oropharyngeal Carcinoma

Under common therapeutic regimens, the prognosis of human papillomavirus (HPV)–positive squamous oropharyngeal carcinomas (OPCs) is more favorable than HPV-negative OPCs. However, the prognosis of some tumors is dismal, and validated prognostic factors are missing in clinical practice. The present work aimed to validate the prognostic significance of our published three-cluster model and to compare its prognostic value with those of the 8^th edition of the tumor-node-metastasis staging system (TNM8) and published signatures and clustering models.

Therapeutic strategies for systemic therapies of human papillomavirus-related oropharyngeal cancer

Our understanding of human papillomavirus (HPV)-related oropharyngeal cancer (OPC) and its molecular basis continues to evolve and produce important insights into customized therapeutic strategies. Novel therapeutics exploiting HPV-related targets are being evaluated in the incurable setting, while the favorable prognosis of locoregionally advanced disease has stimulated investigation into de-escalation strategies. There is much opportunity for better personalization of standard therapy according to HPV status. This review discusses both current and investigational therapeutic strategies for HPV-related OPC.