Smoking and HPV infection are known causes for the vast majority of head and neck squamous cell carcinomas (HNSCC) due to their likelihood of causing gene dysregulation and genomic alterations. Enhancer RNAs (eRNAs) are non-coding RNAs that are known to increase nearby and target gene expression, and activity that has been suggested to be affected by genetic and epigenetic alterations. Here we sought to identify the effects of smoking and HPV status on eRNA expression in HNSCC tumors. We focused on four patient cohorts including smoking/HPV+, smoking/HPV-, non-smoking/HPV+, and non-smoking/HPV- patients. We used TCGA RNA-seq data from cancer tumors and adjacent normal tissue, extracted eRNA read counts, and correlated these to survival, clinical variables, immune infiltration, cancer pathways, and genomic alterations. We found a large number of differentially expressed eRNA in each patient cohort. We also found several dysregulated eRNA correlated to patient survival, clinical variables, immune pathways, and genomic alterations. Additionally, we were able to find dysregulated eRNA nearby seven key HNSCC-related oncogenes. For example, we found eRNA chr14:103272042-103272430 (eRNA-24036), which is located close to the TRAF3 gene to be differentially expressed and correlated with the pathologic N stage and immune cell populations. Using a separate validation dataset, we performed differential expression and immune infiltration analysis to validate our results from the TCGA data. Our findings may explain the association between eRNA expression, enhancer activity, and nearby gene dysregulation.
Keywords: HNSCC; HPV; eRNA; smoking.