RESUMO
Importance: Molecular subtypes of HPV-associated Head and Neck Squamous Cell Carcinoma (HNSCC), named IMU (immune strong) and KRT (highly keratinized), are well-recognized and have been shown to have distinct mechanisms of carcinogenesis, clinical outcomes, and potentially differing optimal treatment strategies. Currently, no standardized method exists to subtype a new HPV+ HNSCC tumor. Our paper introduces a machine learning-based classifier and webtool to reliably subtype HPV+ HNSCC tumors using the IMU/KRT paradigm and highlights the importance of subtype in HPV+ HNSCC. Objective: To develop a robust, accurate machine learning-based classification tool that standardizes the process of subtyping HPV+ HNSCC, and to investigate the clinical, demographic, and molecular features associated with subtype in a meta-analysis of four patient cohorts. Data Sources: We conducted RNA-seq on 67 HNSCC FFPE blocks from University of Michigan hospital. Combining this with three publicly available datasets, we utilized a total of 229 HPV+ HNSCC RNA-seq samples. All participants were HPV+ according to RNA expression. An ensemble machine learning approach with five algorithms and three different input training gene sets were developed, with final subtype determined by majority vote. Several additional steps were taken to ensure rigor and reproducibility throughout. Study Selection: The classifier was trained and tested using 84 subtype-labeled HPV+ RNA-seq samples from two cohorts: University of Michigan (UM; n=18) and TCGA-HNC (n=66). The classifier robustness was validated with two independent cohorts: 83 samples from the HPV Virome Consortium and 62 additional samples from UM. We revealed 24 of 39 tested clinicodemographic and molecular variables significantly associated with subtype. Results: The classifier achieved 100% accuracy in the test set. Validation on two additional cohorts demonstrated successful separation by known features of the subtypes. Investigating the relationship between subtype and 39 molecular and clinicodemographic variables revealed IMU is associated with epithelial-mesenchymal transition (p=2.25×10-4), various immune cell types, and lower radiation resistance (p=0.0050), while KRT is more highly keratinized (p=2.53×10-8), and more likely female than IMU (p=0.0082). Conclusions and Relevance: This study provides a reliable classifier for subtyping HPV+ HNSCC tumors as either IMU or KRT based on bulk RNA-seq data, and additionally, improves our understanding of the HPV+ HNSCC subtypes.
RESUMO
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.