HPV status is associated with altered PIWI-interacting RNA expression pattern in head and neck cancer

Abstract

As HPV-induced cases of oral malignancy increase, it is important to understand the molecular differences between HPV positive and negative head and neck squamous cell carcinoma (HNSCC). PIWI-interacting RNAs (piRNAs) are a class of small non-coding RNAs aberrantly expressed in cancer. We analyzed global piRNA expression patterns to define the HNSCC piRNA transcriptome and assess whether HPV infection status associates with changes in piRNA levels. A total of 498 HNSCC small RNA sequencing libraries were acquired from the Cancer Genomics Hub (cgHUB) Data Repository and a custom sequence analysis pipeline was developed to deduce piRNA expression from raw sequencing data. Expression matrices were aligned to clinicopathological features in order to analyze piRNA expression patterns across different HNSCC groups. The association of a piRNA signature with HPV-positive patient survival was evaluated using a Cox proportional hazard model. Analysis of piRNA levels between HNSCC and non-malignant tissues revealed distinct expression patterns, with 87 piRNAs exclusively expressed in tumor samples. HPV infection status affected the expression of 41 of these piRNAs. Eleven (26.8%) piRNAs were significantly downregulated in HPV16/18 tumors compared to other HPV types. Remarkably, expression of a combination of five-piRNAs in HPV-positive HNSCC tumors was associated with worse overall survival. The expression of specific piRNAs is deregulated in HNSCC, and changes with both HPV status and type. Importantly, a five-piRNA signature is able to delineate a subset of HPV-positive HNSCC patients with poor outcome, highlighting the potential utility of piRNAs in patient management.