Human papillomavirus (HPV) causes an increasing number of head and neck squamous cell carcinomas (HNSCCs). Altered metabolism contributes to patient prognosis, but the impact of HPV status on HNSCC metabolism remains relatively uncharacterized. We hypothesize that metabolism-related gene expression differences unique to HPV-positive HNSCC influences patient survival. The Cancer Genome Atlas RNA-seq data from primary HNSCC patient samples were categorized as 73 HPV-positive, 442 HPV-negative, and 43 normal-adjacent control tissues. We analyzed 229 metabolic genes and identified numerous differentially expressed genes between HPV-positive and negative HNSCC patients. HPV-positive carcinomas exhibited lower expression levels of genes involved in glycolysis and higher levels of genes involved in the tricarboxylic acid cycle, oxidative phosphorylation, and β-oxidation than the HPV-negative carcinomas. Importantly, reduced expression of the metabolism-related genes SDHC, COX7A1, COX16, COX17, ELOVL6, GOT2, and SLC16A2 were correlated with improved patient survival only in the HPV-positive group. This work suggests that specific transcriptional alterations in metabolic genes may serve as predictive biomarkers of patient outcome and identifies potential targets for novel therapeutic intervention in HPV-positive head and neck cancers.

Oncoproteins from high-risk human papillomaviruses (HPV) downregulate the transcription of the class I major histocompatibility complex (MHC-I) antigen presentation apparatus in tissue culture model systems. This could allow infected or transformed cells to evade the adaptive immune response. Using data from over 800 human cervical and head & neck tumors from The Cancer Genome Atlas (TCGA), we determined the impact of HPV status on the mRNA expression of all six MHC-I heavy chain genes, and the β2 microglobulin light chain. Unexpectedly, these genes were all expressed at high levels in HPV positive (HPV+) cancers compared with normal control tissues. Indeed, many of these genes were expressed at significantly enhanced levels in HPV+ tumors. Similarly, the transcript levels of several other components of the MHC-I peptide-loading complex were also high in HPV+ cancers. The coordinated expression of high mRNA levels of the MHC-I antigen presentation apparatus could be a consequence of the higher intratumoral levels of interferon γ in HPV+ carcinomas, which correlate with signatures of increased infiltration by T- and NK-cells. These data, which were obtained from both cervical and oral tumors in large human cohorts, indicates that HPV oncoproteins do not efficiently suppress the transcription of the antigen presentation apparatus in human tumors.

Although infection with human papillomavirus (HPV) is associated with nearly all cervical cancers (CC), a small proportion are HPV-negative. Recently, it has become clear that HPV-negative CC represent a distinct disease phenotype compared to HPV-positive disease and exhibit increased mortality. In addition, variations between different HPV types associated with CC have been linked to altered molecular pathology and prognosis. We compared the immune microenvironments of CC caused by HPV α9 species (HPV16-like), HPV α7 species (HPV18-like) and HPV-negative disease. HPV-negative CC appeared distinct from other subtypes, with greatly reduced levels of lymphocyte infiltration compared to either HPV α9 or α7 CC. Besides reduced levels of markers indicative of B, T, and NK lymphocytes, the expression of T-cell effector molecules, activation/exhaustion markers, and T-cell receptor diversity were also significantly lower in HPV-negative CC. Interestingly, HPV-negative CC expressed much higher levels of potential neoantigens than HPV-positive CC. These results identify profound differences between the immune landscape of HPV-positive and HPV-negative CC as well as modest differences between HPV α9 and α7 CC. These differences may contribute to altered patient outcomes between HPV-negative and HPV-positive CC and potentially between CC associated with different HPV types.

High-risk human papillomaviruses (HPVs) are responsible for a subset of head and neck squamous cell carcinomas (HNSCC). Expression of class II major histocompatibility complex (MHC-II) is associated with antigen presenting cells (APCs). During inflammation, epithelial cells can be induced to express MHC-II and function as accessory APCs. Utilizing RNA-seq data from over 500 HNSCC patients from The Cancer Genome Atlas, we determined the impact of HPV-status on the expression of MHC-II genes and related genes involved in their regulation, antigen presentation, and T-cell co-stimulation. Expression of virtually all MHC-II genes was significantly upregulated in HPV+ carcinomas compared to HPV− or normal control tissue. Similarly, genes that encode products involved in antigen presentation were also significantly upregulated in the HPV+ cohort. In addition, the expression of CIITA and RFX5—regulators of MHC-II—were significantly upregulated in HPV+ tumors. This coordinated upregulation of MHC-II genes was correlated with higher intratumoral levels of interferon-gamma in HPV+ carcinomas. Furthermore, genes that encode various co-stimulatory molecules involved in T-cell activation and survival were also significantly upregulated in HPV+ tumors. Collectively, these results suggest a previously unappreciated role for epithelial cells in antigen presentation that functionally contributes to the highly immunogenic tumor microenvironment observed in HPV+ HNSCC.

Epstein–Barr virus (EBV) is responsible for approximately 9% of stomach adenocarcinomas. EBV-encoded microRNAs have been reported as reducing the function of the class I major histocompatibility complex (MHC-I) antigen presentation apparatus, which could allow infected cells to evade adaptive immune responses. Using data from nearly 400 human gastric carcinomas (GCs), we assessed the impact of EBV on MHC-I heavy and light chain mRNA levels, as well as multiple other components essential for antigen processing and presentation. Unexpectedly, mRNA levels of these genes were as high, or higher, in EBV-associated gastric carcinomas (EBVaGCs) compared to normal control tissues or other GC subtypes. This coordinated upregulation could have been a consequence of the higher intratumoral levels of interferon γ in EBVaGCs, which correlated with signatures of increased infiltration by T and natural killer (NK) cells. These results indicate that EBV-encoded products do not effectively reduce mRNA levels of the MHC-I antigen presentation apparatus in human GCs.