Head and neck squamous cell carcinoma (HNSCC) is one of the ten most common cancers. Most cancer cases originate from alcohol and tobacco consumption. However, studies have demonstrated that human papillomavirus (HPV) infection, particularly HPV-16, may also significantly influence disease progression. The KRAB-ZNF family of genes is involved in epigenetic suppression, and its involvement in carcinogenesis is the subject of extensive studies. The available literature data demonstrate that they may play different roles, both as tumor suppressors and oncogenes. In this study, six ZNF genes, ZFP28, ZNF132, ZNF418, ZNF426, ZNF540, and ZNF880, were tested using several in silico approaches based on the TCGA and GEO datasets. Our analyses indicate that the expression of the analyzed ZNFs was significantly downregulated in tumor tissues and depended on tumor localization. The expression levels of ZNFs differed between HPV-positive vs. HPV-negative patients depending on the clinical-pathological parameters. More specifically, the patients with higher levels of ZNF418 and ZNF540 showed better survival rates than those with a lower expression. In addition, the level of ZNF540 expression in HPV-positive (HPV(+)) patients was higher than in HPV-negative (HPV(−)) patients (p < 0.0001) and was associated with better overall survival (OS). In conclusion, we demonstrate that ZNF540 expression highly correlates with HPV infection, which renders ZNF540 a potential biomarker for HNSCC prognosis and treatment.

Long non-coding RNAs have proven to be important molecules in carcinogenesis. Due to little knowledge about them, the molecular mechanisms of tumorigenesis are still being explored. The aim of this work was to study the effect of ionizing radiation on the expression of lncRNAs in head and neck squamous cell carcinoma (HNSCC) in patients responding and non-responding to radiotherapy. The experimental model was created using a group of patients with response (RG, n = 75) and no response (NRG, n = 75) to radiotherapy based on the cancer genome atlas (TCGA) data. Using the in silico model, statistically significant lncRNAs were defined and further validated on six HNSCC cell lines irradiated at three different doses. Based on the TCGA model, C10orf55, C3orf35, C5orf38, CASC2, MEG3, MYCNOS, SFTA1P, SNHG3, and TMEM105, with the altered expression between the RG and NRG were observed. Analysis of pathways and immune profile indicated that these lncRNAs were associated with changes in processes, such as epithelial-to-mesenchymal transition, regulation of spindle division, and the p53 pathway, and differences in immune cells score and lymphocyte infiltration signature score. However, only C10orf55, CASC2, and SFTA1P presented statistically altered expression after irradiation in the in vitro model. In conclusion, the expression of lncRNAs is affected by ionization radiation in HNSCC, and these lncRNAs are associated with pathways, which are important for radiation response and immune response. Potentially presented lncRNAs could be used as biomarkers for personalized radiotherapy in the future. However, these results need to be verified based on an in vitro experimental model to show a direct net of interactions.

Long non-coding RNAs (lncRNAs) consist of at least 200 nucleotides. Although these molecules do not code proteins, they carry many regulatory functions in normal cells, as well as in cancer cells. For instance, many of these molecules have been previously correlated with tumorigenesis of different cancers and their reaction to various stress factors, such as radiotherapy, chemotherapy, or reactive oxygen species (ROS). The lncRNAs are associated not only with dysregulation in cancers after applied treatment but also with beneficial effects that may be achieved by modulating their expression, often significantly enhancing the patients’ outcomes. A multitude of these molecules was previously considered as potential biomarkers of tumor development, progression, or cells’ response to radio- or chemotherapy. Irradiation, which is often used in treating numerous cancer types, is not always sufficient due to cells gaining resistance in multiple ways. In this review, studies considering lncRNAs and their reaction to radiotherapy were examined. These molecules were divided regarding their role in specific processes strictly related to irradiation, and their influence on this type of treatment was explained, showing how vast an impact they have on IR-supported combat with the disease. This review aims to shed some light on potential future lncRNA-based biomarkers and therapeutic targets.

HPV infection is one of the most important risk factors for head and neck squamous cell carcinoma among younger patients. YRNAs are short non-coding RNAs involved in DNA replication. YRNAs have been found to be dysregulated in many cancers, including head and neck squamous cell carcinoma (HNSCC). In this study, we investigated the role of YRNAs in HPV-positive HNSCC using publicly available gene expression datasets from HNSCC tissue, where expression patterns of YRNAs in HPV(+) and HPV(−) HNSCC samples significantly differed. Additionally, HNSCC cell lines were treated with YRNA1-overexpressing plasmid and RNA derived from these cell lines was used to perform a NGS analysis. Additionally, a deconvolution analysis was performed to determine YRNA1’s impact on immune cells. YRNA expression levels varied according to cancer pathological and clinical stages, and correlated with more aggressive subtypes. YRNAs were mostly associated with more advanced cancer stages in the HPV(+) group, and YRNA3 and YRNA1 expression levels were found to be correlated with more advanced clinical stages despite HPV infection status, showing that they may function as potential biomarkers of more advanced stages of the disease. YRNA5 was associated with less-advanced cancer stages in the HPV(−) group. Overall survival and progression-free survival analyses showed opposite results between the HPV groups. The expression of YRNAs, especially YRNA1, correlated with a vast number of proteins and cellular processes associated with viral infections and immunologic responses to viruses. HNSCC-derived cell lines overexpressing YRNA1 were then used to determine the correlation of YRNA1 and the expression of genes associated with HPV infections. Taken together, our results highlight the potential of YRNAs as possible HNSCC biomarkers and new molecular targets.

MicroRNAs and their role in cancer have been extensively studied for the past decade. Here, we analyzed the biological role and diagnostic potential of miR-154-5p and miR-154-3p in head and neck squamous cell carcinoma (HNSCC). miRNA expression analyses were performed using The Cancer Genome Atlas (TCGA) data accessed from cBioPortal, UALCAN, Santa Cruz University, and Gene Expression Omnibus (GEO). The expression data were correlated with clinicopathological parameters. The functional enrichment was assessed with Gene Set Enrichment Analysis (GSEA). The immunological profiles were assessed using the ESTIMATE tool and RNAseq data from TCGA. All statistical analyses were performed with GraphPad Prism and Statistica. The study showed that both miR-154-5p and miR-154-3p were downregulated in the HNSCC samples and their expression levels correlated with tumor localization, overall survival, cancer stage, tumor grade, and HPV p16 status. GSEA indicated that individuals with the increased levels of miR-154 had upregulated AKT-MTOR, CYCLIN D1, KRAS, EIF4E, RB, ATM, and EMT gene sets. Finally, the elevated miR-154 expression correlated with better immune response. This study showed that miR-154 is highly involved in HNSCC pathogenesis, invasion, and immune response. The implementation of miR-154 as a biomarker may improve the effectiveness of HNSCC treatment.