==== Front Front Oncol Front Oncol Front. Oncol. Frontiers in Oncology 2234-943X Frontiers Media S.A. 10.3389/fonc.2023.1134456 Oncology Original Research LncRNA AL161431.1 predicts prognosis and drug response in head and neck squamous cell carcinoma Zhou Mingzhu 1 † Mao Mingyu 2 † Yang Fan 1 † Zhou Tao 1 * Zhou Liuqing 1 * Li Yuncheng 1 * 1 Department of Otorhinolaryngology, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China 2 Department of Neurology, Shanghai Municipal Hospital of Traditional Chinese Medicine, Shanghai University of Traditional Chinese Medicine, Shanghai, China Edited by: Xuan Zhou, Tianjin Medical University Cancer Institute and Hospital, China Reviewed by: Wei Cao, Shanghai Jiao Tong University, China; Sinan Wang, Tianjin Medical University General Hospital, China; Cheng Wang, Sun Yat-sen University, China *Correspondence: Tao Zhou, entzt2013@sina.com; Liuqing Zhou, 190042422@qq.com; Yuncheng Li, whxhent30@163.com †These authors have contributed equally to this work 16 6 2023 2023 13 113445630 12 2022 01 6 2023 Copyright © 2023 Zhou, Mao, Yang, Zhou, Zhou and Li 2023 Zhou, Mao, Yang, Zhou, Zhou and Li https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms. Background Long non-coding RNAs (lncRNAs) are increasingly recognized as essential players in various biological processes due to their interactions with DNA, RNA, and protein. Emerging studies have demonstrated lncRNAs as prognostic biomarkers in multiple cancers. However, the prognostic effect of lncRNA AL161431.1 in head and neck squamous cell carcinoma (HNSCC) patients has not been reported. Methods In the present study, we conducted a series of analyses to identify and validate the prognostic value of lncRNA AL161431.1 in HNSCC, which included differential lncRNAs screening, survival analysis, Cox regression analysis, time ROCanalysis, nomogram prediction, enrichment analysis, tumor infiltration of immune cells, drug sensitivity analysis, and quantitative real-time polymerase chain reaction (qRT-PCR). Results In this study, we performed a comprehensive survival and predictive analysis and demonstrated that AL161431.1 was an independent prognostic factor of HNSCC, for which a high AL161431.1 level indicated poor survival in HNSCC. Functional enrichment analyses found that cell growth and immune-related pathways were significantly enriched in HNSCC, suggesting that AL161431.1 may play a role in tumor development and tumor microenvironment (TME). AL161431.1-related immune cells infiltration analysis demonstrated that AL161431.1 expression is significantly positively associated with M0 macrophages in HNSCC (P<0.001). Using "OncoPredict", we recognized chemotherapy drugs sensitive to the high expression group. Quantitative real-time polymerase chain reaction (qRT-PCR) was performed to identify the expression level of AL161431.1 in HNSCC, and the results further validated our findings. Conclusions Our findings suggest that AL161431.1 is a reliable prognostic marker for HNSCC and can potentially be an effective therapeutic target. AL161431.1 long non-coding RNAs head and neck squamous cell carcinoma cell growth immune cell infiltration drug sensitivity prognosis section-in-acceptanceHead and Neck Cancer ==== Body pmc1 Introduction Head and neck squamous cell carcinoma (HNSCC), which arises from the mucosal epithelium of the oral cavity, nasopharynx, oropharynx, hypopharynx, and larynx, is the eighth most commonly occurring form of cancer in the world, responsible for approximately 745,000 new cases and 364,340 deaths in 2020 (1). The available clinical treatment modalities for HNSCC include surgery, radiotherapy, chemotherapy, and the latest emerging immunotherapy (2). However, despite a modest improvement in HNSCC survival over the past three decades, the 5-year survival rate still hovers at 60% (3). Therefore, there is an urgent need to find a new molecular biomarker capable of predicting survival, identifying new intervention targets, and predicting response to therapeutic agents. Long non-coding RNAs (lncRNAs), defined as RNAs longer than 200 nucleotides with minimal coding potential, are becoming a new biomarker for cancer diagnosis and prognosis (4–6). Its abnormal expression is closely associated with the occurrence and progression of HNSCC (7–9). MYOSLID, ZFAS1, KTN1-AS1, MIR31HG, UCA1, NEAT1, EGFR-AS1, and other lncRNAs have been identified as HNSCC-associated prognostic biomarkers (8, 10–15). At the same time, it has been reported that lncRNA AL161431.1 is involved in the occurrence and development of various tumors, including pancreatic cancer (16), endometrial cancer (17), lung squamous cell carcinoma (18), and non-small cell lung cancer (19). However, the prognostic effect of lncRNA AL161431.1 in HNSCC has never been studied. Moreover, the correlation of AL161431.1 with immune infiltration and drug sensitivity in HNSCC has not been determined. In this study, we conducted a series of analyses, which included differential lncRNAs screening, survival analysis, Cox regression analysis, time ROC analysis, nomogram prediction, enrichment analysis, tumor infiltration of immune cells, drug sensitivity analysis, and quantitative real-time polymerase chain reaction (qRT-PCR). First, we screened the prognosis-related lncRNA AL161431.1 of HNSCC. Then, we confirmed that AL161431.1 was an independent prognostic factor of HNSCC. High expression of AL161431.1 indicated poor survival in HNSCC. Our findings suggest that AL161431.1 has a high potential in the diagnosis, prognosis, and targeted therapy of HNSCC. 2 Materials and methods 2.1 Data collection and processing A total of 271 RNA-sequence (17 normal and 254 tumor tissues) data were acquired from the TCGA-HNSCC database. After removing mRNA data, lncRNAs were used for further analysis. In addition, clinical information of 313 HNSCC patients, including age, gender, grade, stage, T-N-M stage, survival time, and survival status, was extracted for subsequent analyses (https://portal. http://gdc.cancer.gov). 2.2 Analysis of differentially expressed lncRNAs The differential lncRNA expression between normal and HNSCC tissues was assessed using the “Limma” package of R software (version 4.0.5) (20). The filtering criteria for lncRNA differential expression were false discovery rate (FDR) < 0.05 and | Log fold change| >1. In addition, the R package “heatmap” was used to display the 100 differentially expressed lncRNAs. 2.3 Identification of prognosis-related lncRNAs Based on the previously obtained differentially expressed lncRNAs, the Kaplan-Meier analysis (KM) and the Cox regression analysis screened the prognosis-related lncRNAs. By merging expression and clinical data, the lncRNAs with P-values < 0.05 were regarded as independent prognostic lncRNAs and were analyzed further. Subsequently, time-dependent receiver operating characteristics (ROC) analysis was conducted to screen out lncRNAs with high accuracy in predicting overall survival (OS) using the “survival”, “survminer”, and “timeROC” R packages. 2.4 Differential expression and survival analysis LncRNA AL161431.1 was selected from the above prognostic lncRNAs for subsequent analyses. The differential analysis used the R package “limma” to investigate the differential expression of AL161431.1 between normal and tumor subgroups. The pan-cancer analysis compared expression values between cancer tissue and adjacent normal tissue by the Wilcoxon rank-sum test (21). More importantly, the expression values of AL161431.1 in HNSCC samples and normal samples were compared. The above results are presented by boxplots. According to the median level, HNSCC tissues were divided into high- and low- AL161431.1 expression groups. Finally, KM analysis compared the OS between the two sets using the “survival” and “survminer” packages in R. 2.5 Univariate, multivariate Cox regression and time ROC analysis The univariate Cox analysis was used to identify the prognostic indicators. At the same time, the multivariate Cox analysis was used to determine whether AL161431.1 is an independent risk factor for OS in HNSCC. First, the ROC curve assessed the diagnostic value of AL161431.1 in HNSCC by the “pROC” R package. Then, the “timeROC” R package was employed to score the prediction accuracy by AL161431.1 using time-dependent ROC curve analysis. 2.6 Establishment and evaluation of nomograms A nomogram encompassing the expression value of AL161431.1 attributes was created using the R packages “rms” and “regplot” to predict HNSCC patient OS (1 year, 2 years, and 3 years). The nomogram’s predictive power was evaluated using calibration curves. 2.7 Co-expression and enrichment analyses The correlation analysis was used to identify the co-expression gene of AL161431.1 in the TCGA-HNSCC cohort to better understand the important biological function that AL161431.1 plays in HNSCC. Then, all co-expression genes that positively and negatively correlated with AL161431.1 were selected for subsequent enrichment analysis. For correlation analysis, the Pearson correlation test was employed. The corresponding coefficient threshold values were set at >0.4. The adjusted P-value was <0.001. The above results were manifested as scatter plots using the “ggplot2” package and R software (version 4.0.5). Gene Ontology (GO) analyses were carried out to identify the biological effects of AL161431.1 co-expressed genes. In addition, Gene Set enrichment analyses (GSEA) were conducted to identify significant enrichment signaling pathways in the high- and low- AL161431.1 expression groups. The R package “clusterProfiler” was used to evaluate GO analyses and GSEA (22). 2.8 Immune cell infiltration The Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm (23) evaluated immune infiltration in HNSCC tissue in 22 subsets of immune cells, including M0 macrophages, activated dendritic cells, activated mast cells, eosinophils, resting NK cells, CD4 memory resting T cells, neutrophils, M2 macrophages, memory B cells, CD4 naïve T cells, plasma cells, gamma delta T cells, activated NK cells, M1 macrophages, monocytes, CD4 memory activated T cells, naïve B cells, resting dendritic cells, T follicular helper cells (Tfhs), and regulatory T cells (Tregs). Spearman correlation analysis investigated the relationship between AL161431.1 expression and the infiltrated immune cells. Lollipop plots were used to visualize the correlation coefficients of the results (P-value<0.05). 2.9 OncoPredict for drug sensitivity analysis To predict the chemotherapeutic drug sensitivity of the high- and low- AL161431.1 expression groups, we used the “OncoPredict” R package and the Genomics of Drug Sensitivity in Cancer (GDSC; https://www.cancerrxgene.org/) database to predict drug responses in cancer patients (24). First, we analyzed the difference in the activity of drugs between the two groups using the Wilcoxon test. A total of 198 drugs were calculated, and the significance level was set at P <0.001. The “ggplot2” and “ggpubr” functions of R were used to create the box plots. 2.10 Validation of the expression level of AL161431.1 in HNSCC by qRT-PCR A total of 8 HNSCC tissues and 8 adjacent normal mucosal tissues, diagnosed with primary HNSCC after surgical resection in the Wuhan Union Hospital, were collected for subsequent study. There was no history of other tumors, and no patients received radiotherapy, chemotherapy, or other treatments. All patients had signed informed consent before surgery. The collected tissue samples were frozen in liquid nitrogen and stored at -80°C for further RNA extraction. First, total RNA was extracted from samples using TRIzol reagents (Invitrogen, Waltham, MA). Then, reverse transcription and PCR reactions were performed using the PrimeScriptTM RT-PCR kit (Takara, Shiga, Japan). The primer sequence was as follows: AL161431.1: Forward 5 -TACACTGTTTCTCCAAGCCATCA -3′ Reverse 5 -GAATTGGGAGGATCTAGGACATCTA -3′ 2.11 Statistical analysis All analyses were performed using R version 4.0.5, and a P-value<0.0 was considered statistically significant. 3 Results 3.1 Differential analysis of lncRNA in HNSCC To identify the expression difference of AL161431.1 in tumors and normal tissues, we analyzed data from 271 patients (17 normal and 254 tumors) extracted from the TCGA-HNSCC database. Following differential expression analysis, 1180 differentially expressed lncRNAs were identified. The study design of this project is shown in Figure 1 . In addition, the differentially expressed lncRNAs have been presented as a heatmap, including only 100 differentially expressed lncRNAs ( Figure 2A ). Figure 1 The flow chart of our study. Figure 2 (A) Heatmaps of 100 differential lncRNA in normal and tumor samples. (B) The scatter diagram indicates the comparative expression of AL161431.1 between pan-cancers and corresponding normal tissues. (C) Boxplots present the differential expression of LncRNA AL161431.1 between normal and HNSCC subtypes. (D) Kaplan–Meier survival curves show that patients in the AL161431.1 low-level group survived dramatically longer than those in the high-level group. (E) Evaluation of the expression of AL161431.1 in HNSCC tissues and adjacent normal tissues (n = 16). The statistical significance is shown as ns, P > 0.05; *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001. 3.2 Screening and identification of prognosis-related lncRNAs in HNSCC To identify lncRNAs associated with the prognosis of HNSCC, we used the Cox regression and KM survival analysis to preliminarily screen out 137 prognostic-related lncRNAs (P <0.05). Then, 81 independent prognostic lncRNAs were identified by independent prognostic analysis (P<0.05). Finally, we used the ROC curve to select four highly accurate lncRNAs in predicting patient survival, namely SNHG26, AL161431.1, LINC00460, and AL358334.2 (AUC>0.65, Table 1 ). Table 1 Four independent prognostic lncRNAs. lncRNA SNHG26 AL161431.1 LINC00460 AL358334.2 AUC 0.674476301869521 0.672933769930221 0.67979996009358 0.678265369849573 3.3 Expression of AL161431.1 was upregulated in HNSCC and validated by qRT-PCR To explore the differential expression of AL161431.1 in the normal and tumor samples, we first analyzed its expression in pan-cancer, which included HNSCC. The results demonstrated that compared with normal tissue, the AL161431.1 expression was significantly upregulated in 11 different types of cancer, namely, urothelial bladder carcinoma (BLCA), breast invasive carcinoma (BRCA), cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), esophageal carcinoma (ESCA), HNSCC, lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), rectum adenocarcinoma(READ), stomach adenocarcinoma (STAD), and uterine corpus endometrial carcinoma (UCEC) (all P<0.01, Figure 2B ). In contrast, the AL161431.1 expression was significantly downregulated in kidney Chromophobe (KICH). Further investigation revealed that HNSCC tissues had significantly higher levels of AL161431.1 than normal tissues (P < 0.001, Figure 2C ). At the same time, KM survival analysis found that high AL161431.1 expression was closely associated with worse OS in HNSCC (P< 0.001, Figure 2D ). In addition, the expression level of AL161431.1 in HNSCC tissues was verified by qRT-PCR. qRT-PCR analysis showed that AL161431.1 was upregulated in HNSCC compared to normal tissues ( Figure 2E ). 3.4 AL161431.1 was an independent prognostic indicator for HNSCC We conducted the univariate Cox regression analysis to examine the predictive value of AL161431.1 in HNSCC patients. As illustrated in Figure 3A , the result revealed age [HR: 1.026 (1.003-1.051), P = 0.030], gender [HR: 0.585 (0.353-0.970), P = 0.038], stage [HR: 1.569 (1.101-2.235), P =0.013], and AL161431.1 expression [hazard ratio (HR): 1.021 (1.002-1.041), P = 0.030] as the significant predictors of OS in HNSCC patients. After multivariate adjustment, only AL161431.1 [HR: 1.022 (1.002-1.042), P = 0.028], gender [HR: 0.545 (0.322-0.924), P = 0.024], and stage [HR: 1.661(1.162-2.376), P = 0.005] remained significant, as displayed in Figure 3B . The results above manifested that AL161431.1 is an independent prognostic factor for HNSCC. As shown in Figure 3C , the AUC of AL161431.1 was 0.749, indicating a high diagnostic value of AL161431.1. Furthermore, time-dependent ROC curve analysis revealed that AL161431.1 had a high prognostic value, which the 1 -, 3 -, and 5-year survival rates predicted by the AL161431.1 (AUC) were 0.559, 0.674, and 0.673, respectively ( Figure 3D ). Figure 3 (A, B) Univariate and multivariate Cox analyses evaluated the independent prognostic value of AL161431.1 risk signature in HNSCC patients. (C) ROC curves show the diagnosis accuracy of AL161431.1. The area under the curve is 0.749. (D) ROC curve for AL161431.1 for predicting 1-, 3- and 5-year survival. (E) The prediction of 1-, 3-, and 5-year survival for HNSCC patients based on the prognostic nomogram derived from the expression of AL161431.1 and other clinicopathologic features. (F) Calibration curve of the nomogram. 3.5 Construction and assessment of the prognosis nomogram To precisely estimate the 1-, 3-, and 5-year survival rates, we constructed a nomogram by combining the expression value of AL161431.1 with multiple clinicopathological factors (age, grade, M stage, gender, AJCC stage, N stage, and T stage) ( Figure 3E ). The calibration curve analysis revealed that the predicted and actual 1-, 3- and 5-year survival times were comparable ( Figures 3F ). These findings suggested that the nomogram containing the expression level of AL161431.1 is accurate and reliable. 3.6 Co-expression and enrichment analyses reveal that AL161431.1 may involve in the regulation of cell growth We first identified 443 co-expressed genes significantly associated with AL161431.1 to predict the functions and pathways that AL161431.1 may affect ( Figures 4A, B ). Then, we selected these genes to perform GO enrichment analyses, and results revealed AL161431.1 to be mainly associated with cell growth in the biological process (BP) category, the chromosomal region in the cellular component (CC) category, and transcription coregulator activity in molecular function (MF) category ( Figures 5A–C ). In addition, GSEA was conducted to search pathways AL161431.1 might affect. The GSEA results showed that allograft rejection, antigen procession and presentation, autoimmune thyroid disease, the intestinal immune network for IgA production, and type I diabetes mellitus were significantly enriched in the low-AL161431.1 expression group ( Figure 5D ). Figure 4 (A, B) Scatter plots show the correlation of AL161431.1 expression with the expression of the Co-expression gene. AL161431.1 expression was positively (A) or negatively (B) correlated with the expression of the Co-expression gene. Figure 5 (A, B) The bar plot (A) and dot plot (B) of GO analysis for the top ten enrichments. (C) GO analysis was performed to detect biological processes involved in AL161431.1 regulated genes. (D) GSEA results show significant enrichment signaling pathways. 3.7 The relationship between AL161431.1 expression and the level of immune infiltration in HNSCC Immune cells that play an essential role in resisting or accelerating tumor growth are tightly linked to the initiation, development, and progression of HNSCC (25). Therefore, we evaluated immune cells significantly associated with AL161431.1 expression and their immune infiltration level in HNSCC. Correlation analysis revealed that the level of infiltration of M0 macrophages, activated dendritic cells, activated mast cells, and eosinophils were significantly and positively correlated with AL161431.1 expression. However, AL161431.1 expression was significantly negatively correlated with the infiltration level of regulatory T cells (Tregs), T follicular helper cells (Tfhs), resting dendritic cells, CD8 T cells, resting mast cells, and naive B cells (all P< 0.05, Figure 6 ). Figure 6 Lollipop charts show the correlation between AL161431.1 and immune cell infiltration in HNSCC samples. 3.8 Drug responses of high- and low- AL161431.1 expression groups in HNSCC There were statistically significant response differences between high- and low-AL161431.1 expression groups for 59 of these drugs (P < 0.001). Apart from BI−2536 and PD0325901, several chemotherapeutic drugs, including AZD1332, AZD3759, Dasatinib, Erlotinib, Gefitinib, Lapatinib, Obatoclax Mesylate, Sapitinib, SCH772984, and VSP34_8731, showed a lower score in the high-AL161431.1 expression group (all P< 0.001, Figure 7 ). This suggested that HNSCC patients in the high-AL161431.1 expression group were more likely to respond to these chemotherapeutic drugs. Figure 7 Boxplots present the differences in chemotherapeutic sensitivity between the AL161431.1 low-level and high-level groups. 4 Discussion In the past, the major difficulties in treating HNSCC were late detection, poor response to treatment, and a high recurrence rate. In addition, the molecular heterogeneity of HNSCC has hindered the identification of specific targets and the development of targeted therapy. EGFR inhibitors, as the only approved targeted drugs, have limited efficacy and face the problem of tumor drug resistance (2, 26, 27). In recent years, immunotherapy has brought new hope to patients with HNSCC (28). However, HNSCC remains difficult to treat, with high mortality and poor survival. As a result, determining how to diagnose the disease early and predict the patient’s prognosis is a main priority. With the deepening of research, more and more lncRNAs have been proven to be novel diagnostic and prognostic markers and molecularly targeted therapeutic targets (29–31). This prompted us to search for lncRNAs related to the prognosis of HNSCC to predict patients’ survival and treatment better. This study comprehensively elucidated the prognostic value of AL161431.1 in patients with HNSCC by bioinformatics methods. We confirmed that the AL161431.1 expression was upregulated in HNSCC neoplastic tissues compared to normal tissues. In addition, we elucidated that AL161431.1 is an independent prognostic indicator for deteriorative OS in HNSCC patients using KM survival analysis combined with Cox regression analysis. Besides, gender and stage were also independent prognostic factors for HNSCC. We used differentially expressed lncRNA, Kaplan-Meier analysis, and the Cox regression analysis to screen the prognosis-related lncRNAs, which differs from Cao W et al.’s methods, such as the orthogonal partial least squares discrimination analysis (OPLS-DA) and 1.5-fold expression change criterion methods (32). In our study, we first screened out lncRNAs with different expressions in normal and tumor samples. Then, the Kaplan-Meier and Cox regression analyses were performed to screen the prognosis-related lncRNA. However, Cao W et al. first screened out lncRNAs associated with the OS of the patients. Then, the patients were divided into good and poor survival groups based on survival time. Next, the OPLS-DA analysis evaluated the lncRNA expression profile differences between the groups. Finally, the lncRNAs with significant differences between groups were further screened by the 1.5x expression variation criterion. As a non-parametric estimation method, Kaplan-Meier analysis is currently the most commonly used method for survival analysis. It can analyze the impact of univariate and categorical variables on overall survival and provide a visual representation of survival function (33). Cox regression analysis can handle multiple predictors and confounding variables and provide hazard ratios, which help quantify the effect of predictors on survival (33, 34). Combining these two methods, we established a reliable method for screening the prognosis-related lncRNAs. Many studies have used these two methods to screen prognostic-related lncRNAs and verify the prognostic capability of models. For example, Zhang et al. used differentially expressed lncRNA screening, univariate survival analysis, and multivariate Cox regression analysis to identify a 4-lncRNA signature, which performed well in predicting the prognosis of laryngeal cancer (35). Using univariate and multivariate Cox regression analyses, Wu et al. screened prognosis-related lncRNAs, and an eight-immune-related lncRNA prognostic signature was acquired. Then, Kaplan-Meier analysis and ROC analysis verified the prognostic capability of models (36). Similarly, OPLS-DA and 1.5-fold expression change criterion methods are innovative and effective methods for screening prognostic lncRNAs, which can remove irrelevant variables and effectively separate samples between groups. We used different methods to screen for prognosis-related lncRNAs, but both yielded reliable results. In previous studies, numerous lncRNAs were aberrantly expressed in various tumors (37, 38). Studies have proved that AL161431.1 is implicated as an oncogene in various tumor types. For example, AL161431.1 was overexpressed in endometrial carcinoma and promoted endometrial carcinoma cell multiplication and migration by activating miR-1252-5p expression (17). Furthermore, LncRNA AL161431.1 could influence the invasion and metastasis of pancreatic cancer cells by promoting the epithelial-mesenchymal transition (EMT) process (16). Studies have also shown that AL161431.1 is an autophagy-related LncRNA, strongly correlated with the prognosis and tumor immune microenvironment of non-small cell lung cancer (19). Shao et al. found that hypoxia-related lncRNA AL161431.1 was highly expressed in LUAD. Inhibition of AL161431.1 can block the proliferation of LUAD cells. The expression level of AL161431.1 was upregulated under hypoxia (39). However, the precise role of AL161431.1 in HNSCC is unknown. In the present study, through GO enrichment analysis, we found that AL161431.1 is involved in a variety of biological functions and processes, such as cell growth, negative regulation of cell growth, and extrinsic apoptotic signaling pathway. In light of these data, AL161431.1 may affect the growth and development of HNSCC, which may be the reason for the poor prognosis of patients with high expression of AL161431.1. However, its specific role and molecular mechanism remain to be studied. HNSCC is an immunosuppressive disease with high immune infiltration and poor antigen-presenting function (28, 40). The GSEA results in this study showed that several immune-related pathways are enriched in the low-AL161431.1 expression group. Based on these results, we concluded that immunotherapy might be more effective for patients in the low-AL161431.1 expression group than in the high-AL161431.1 expression group. Immune cells are the cellular basis of immunotherapy. Infiltrating immune cells in tumors is closely related to clinical outcomes and can be used as drug targets for cancer treatment (41). Our correlation analysis of AL161431.1 expression and immune cell infiltration revealed that AL161431.1 expression correlated significantly and positively with M0 macrophages, eosinophils, activated mast cells, and activated dendritic cells. In contrast, regulatory T cells (Tregs), T follicular helper cells (Tfhs), naive B cells, resting mast cells, CD8 T cells, and resting dendritic cells were significantly negatively correlated with it. Studies have shown that lncRNA can activate the expression of immune-cell-related genes, thus leading to tumor immune cell infiltration (42). Meanwhile, immune cells, such as macrophages, dendritic cells, regulatory T cells, and mast cells, are essential to the tumor microenvironment (TME) that can positively or negatively regulate cancer development (43). Macrophages within the TME are called tumor-associated macrophages (TAMs) (44). The latest research has shown that TAMs infiltration often predicts a poor prognosis (45, 46). On the contrary, T cell infiltration, especially CD8 T cells, predicts a favorable prognosis (47). Our study results support these conclusions. Tregs are the most important type of cell in the TME. Typically, Treg cells are dedicated to regulating immune responses to prevent excessive reactivity of the immune system. However, Tregs’ function within the TME is highly complex, as they can promote cancer progression by suppressing the immune response against cancer cells (48). Considerable evidence suggests that Tregs increase during the development of HNSCC, but there is no consensus on the prognostic value of Tregs in HNSCC (49). Several studies strongly support that high levels of Tregs are negatively correlated with the prognosis of HNSCC (50, 51). However, other studies describe a completely different situation, where high Treg numbers are associated with improved overall survival (51, 52). Difficulties in correctly identifying Treg cells may cause these conflicting results. Similarly, this also explains our results that AL161431.1 expression correlated significantly and negatively with Treg cells. Based on these findings, we can better evaluate patients’ prognosis and immunotherapeutic efficacy based on the type and status of infiltrating immune cells. Hypoxia is a common cause of cancer cell proliferation, metastasis, and recurrence. In addition, hypoxia is closely associated with increased drug resistance in tumors, severely limiting the therapeutic efficacy of HNSCC (53, 54). Increasing data have shown that hypoxic-induced activation of hypoxia-inducible factors (HIFs) can regulate the expression of lncRNA and participate in tumor proliferation, metastasis, and drug resistance (55–57). AL161431.1, a hypoxia-related lncRNA, is upregulated under hypoxia conditions, associated with poor tumor patient survival (39). To improve treatment outcomes, we identified potential chemotherapeutic-sensitive drugs for HNSCC patients with high expression of AL161431.1 by analyzing the differences in sensitivity to chemotherapy drugs between the two sets. We found that patients in the high-AL161431.1 expression group may respond better to chemotherapeutic agents, such as BI−2536, PD0325901, and AZD1332, than those in the low-AL161431.1 expression group. These results may guide the chemotherapy treatment of HNSCC. However, the specific role of AL161431.1 in the difference in chemotherapeutic drug sensitivity in HNSCC warranted further analysis. Finally, we validated the expression of AL161431.1 in clinical specimens by qRT-PCR. Consistent with our previous bioinformatics analysis, the expression of AL161431.1 was upregulated in HNSCC. To our knowledge, this is the first study on AL161431.1 and HNSCC prognosis, including differential and survival analyses. Our findings suggest that AL161431.1 has a high potential in the diagnosis, prognosis, and targeted therapy of HNSCC. However, our study has several limitations. First, HNSCC has been defined as a highly heterogenous tumor (58). Therefore, the clinical impact of AL161431.1 in other specific subtypes of head and neck cancer remains unclear. Second, we only analyzed the relationship between the AL161431.1 expression and the level of immune infiltration in HNSCC. However, its relationship with the immunosuppressive microenvironment of HNSCC remains to be further studied. Finally, additional experimentation and clinical research are required to further validate the clinical value and potential mechanism of AL161431.1. 5 Conclusion In conclusion, our findings demonstrated that AL161431.1 is highly expressed and associated with a poor prognosis in HNSCC, which might make it a new therapeutic target for HNSCC patients. Data availability statement The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material. Author contributions MZ, FY, and LZ conceived and performed bioinformatics analysis. MZ and LZ collected clinical samples. MZ, MM, and FY co-wrote the manuscript. LZ, TZ, and YL undertook a manuscript review. All authors contributed to the article and approved the submitted version. Acknowledgments Thanks to all authors for their contributions to this study, especially to LZ for her guidance and assistance. Conflict of interest The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. Publisher’s note All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher. ==== Refs References 1 Hyuna S Jacques F Siegel RL Mathieu L Isabelle S Ahmedin J . Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: Cancer J Clin (2021) 71 (3 ). doi: 10.3322/caac.21660 2 Alsahafi E Begg K Amelio I Raulf N Lucarelli P Sauter T . Clinical update on head and neck cancer: molecular biology and ongoing challenges. Cell Death Dis (2019) 10 (8 ):540. doi: 10.1038/s41419-019-1769-9 31308358 3 Pulte D Brenner H . Changes in survival in head and neck cancers in the late 20th and early 21st century: a period analysis. Oncologist (2010) 15 (9 ):994–1001. doi: 10.1634/theoncologist.2009-0289 20798198 4 Statello L Guo CJ Chen LL Huarte M . Gene regulation by long non-coding RNAs and its biological functions. Nat Rev Mol Cell Biol (2021) 22 (2 ):96–118. doi: 10.1038/s41580-020-00315-9 33353982 5 Bhan A Soleimani M Mandal SS . Long non-coding RNA and cancer: a new paradigm. Cancer Res (2017) 77 (15 ):3965–81. doi: 10.1158/0008-5472.CAN-16-2634 6 Hou P Meng S Li M Lin T Chu S Li Z . LINC00460/DHX9/IGF2BP2 complex promotes colorectal cancer proliferation and metastasis by mediating HMGA1 mRNA stability depending on m6A modification. J Exp Clin Cancer Res (2021) 40 (1 ):52. doi: 10.1186/s13046-021-01857-2 33526059 7 Jiang Y Cao W Wu K Qin X Wang X Li Y . LncRNA LINC00460 promotes EMT in head and neck squamous cell carcinoma by facilitating peroxiredoxin-1 into the nucleus. J Exp Clin Cancer Res (2019) 38 (1 ):365. doi: 10.1186/s13046-019-1364-z 31429766 8 Wang R Ma Z Feng L Yang Y Tan C Shi Q . LncRNA MIR31HG targets HIF1A and P21 to facilitate head and neck cancer cell proliferation and tumorigenesis by promoting cell-cycle progression. Mol Cancer (2018) 17 (1 ):162. doi: 10.1186/s12943-018-0916-8 30458787 9 Yang Y Chen D Liu H Yang K . Increased expression of lncRNA CASC9 promotes tumor progression by suppressing autophagy-mediated cell apoptosis via the AKT/mTOR pathway in oral squamous cell carcinoma. Cell Death Dis (2019) 10 (2 ):41. doi: 10.1038/s41419-018-1280-8 30674868 10 Xiong HG Li H Xiao Y Yang QC Yang LL Chen L . Long non-coding RNA MYOSLID promotes invasion and metastasis by modulating the partial epithelial-mesenchymal transition program in head and neck squamous cell carcinoma. J Exp Clin Cancer Res (2019) 38 (1 ):278. doi: 10.1186/s13046-019-1254-4 31238980 11 Kolenda T Guglas K Kopczyńska M Teresiak A Bliźniak R Mackiewicz A . Oncogenic role of ZFAS1 lncRNA in head and neck squamous cell carcinomas. Cells (2019) 8 (4 ). doi: 10.3390/cells8040366 12 Jiang Y Wu K Cao W Xu Q Wang X Qin X . Long non-coding RNA KTN1-AS1 promotes head and neck squamous cell carcinoma cell epithelial-mesenchymal transition by targeting miR-153-3p. Epigenomics (2020) 12 (6 ):487–505. doi: 10.2217/epi-2019-0173 32267161 13 Qian Y Liu D Cao S Tao Y Wei D Li W . Upregulation of the long non-coding RNA UCA1 affects the proliferation, invasion, and survival of hypopharyngeal carcinoma. Mol Cancer (2017) 16 (1 ):68. doi: 10.1186/s12943-017-0635-6 28327194 14 Wang P Wu T Zhou H Jin Q He G Yu H . Long non-coding RNA NEAT1 promotes laryngeal squamous cell cancer through regulating miR-107/CDK6 pathway. J Exp Clin Cancer Res (2016) 35 :22. doi: 10.1186/s13046-016-0297-z 26822763 15 Tan DSW Chong FT Leong HS Toh SY Lau DP Kwang XL . Long non-coding RNA EGFR-AS1 mediates epidermal growth factor receptor addiction and modulates treatment response in squamous cell carcinoma. Nat Med (2017) 23 (10 ):1167–75. doi: 10.1038/nm.4401 16 Ma G Li G Fan W Xu Y Song S Guo K . The role of long non-coding RNA AL161431.1 in the development and progression of pancreatic cancer. Front Oncol (2021) 11 :666313. doi: 10.3389/fonc.2021.666313 34395245 17 Gu ZR Liu W . The LncRNA AL161431.1 targets miR-1252-5p and facilitates cellular proliferation and migration via MAPK signaling in endometrial carcinoma. Eur Rev Med Pharmacol Sci (2020) 24 (5 ):2294–302. doi: 10.26355/eurrev_202003_20495 18 Ju Q Zhao YJ Ma S Li XM Zhang H Zhang SQ . Genome-wide analysis of prognostic-related lncRNAs, miRNAs and mRNAs forming a competing endogenous RNA network in lung squamous cell carcinoma. J Cancer Res Clin Oncol (2020) 146 (7 ):1711–23. doi: 10.1007/s00432-020-03224-8 19 Jiang H Xu A Li M Han R Wang E Wu D . Seven autophagy-related lncRNAs are associated with the tumor immune microenvironment in predicting survival risk of non-small cell lung cancer. Brief Funct Genomics (2022) 21 (3 ):177–87. doi: 10.1093/bfgp/elab043 20 Ritchie ME Phipson B Wu D Hu Y Law CW Shi W . Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res (2015) 43 (7 ):e47. doi: 10.1093/nar/gkv007 25605792 21 Li Y Ge X Peng F Li W Li JJ . Exaggerated false positives by popular differential expression methods when analyzing human population samples. Genome Biol (2022) 23 (1 ):79. doi: 10.1186/s13059-022-02648-4 35292087 22 Yu G Wang LG Han Y He QY . clusterProfiler: an r package for comparing biological themes among gene clusters. Omics (2012) 16 (5 ):284–7. doi: 10.1089/omi.2011.0118 23 Newman AM Liu CL Green MR Gentles AJ Feng W Xu Y . Robust enumeration of cell subsets from tissue expression profiles. Nat Methods (2015) 12 (5 ):453–7. doi: 10.1038/nmeth.3337 24 Maeser D Gruener RF Huang RS . oncoPredict: an r package for predicting in vivo or cancer patient drug response and biomarkers from cell line screening data. Brief Bioinform (2021) 22 (6 ). doi: 10.1093/bib/bbab260 25 Elmusrati A Wang J Wang CY . Tumor microenvironment and immune evasion in head and neck squamous cell carcinoma. Int J Oral Sci (2021) 13 (1 ):24. doi: 10.1038/s41368-021-00131-7 34341329 26 Braakhuis BJ Brakenhoff RH Leemans CR . Treatment choice for locally advanced head and neck cancers on the basis of risk factors: biological risk factors. Ann Oncol (2012) 23 Suppl 10 :x173–7. doi: 10.1093/annonc/mds299 27 Sacco AG Worden FP . Molecularly targeted therapy for the treatment of head and neck cancer: a review of the ErbB family inhibitors. Onco Targets Ther (2016) 9 :1927–43. doi: 10.2147/ott.S93720 28 Ferris RL . Immunology and immunotherapy of head and neck cancer. J Clin Oncol (2015) 33 (29 ):3293–304. doi: 10.1200/JCO.2015.61.1509 29 Kolenda T Guglas K Ryś M Bogaczyńska M Teresiak A Bliźniak R . Biological role of long non-coding RNA in head and neck cancers. Rep Pract Oncol Radiother (2017) 22 (5 ):378–88. doi: 10.1016/j.rpor.2017.07.001 30 Chandra Gupta S Nandan Tripathi Y . Potential of long non-coding RNAs in cancer patients: from biomarkers to therapeutic targets. Int J Cancer (2017) 140 (9 ):1955–67. doi: 10.1002/ijc.30546 31 Meng X Wang ZF Lou QY Rankine AN Zheng WX Zhang ZH . Long non-coding RNAs in head and neck squamous cell carcinoma: diagnostic biomarkers, targeted therapies, and prognostic roles. Eur J Pharmacol (2021) 902 :174114. doi: 10.1016/j.ejphar.2021.174114 33901464 32 Cao W Liu JN Liu Z Wang X Han ZG Ji T . A three-lncRNA signature derived from the atlas of ncRNA in cancer (TANRIC) database predicts the survival of patients with head and neck squamous cell carcinoma. Oral Oncol (2017) 65 :94–101. doi: 10.1016/j.oraloncology.2016.12.017 28109476 33 George B Seals S Aban I . Survival analysis and regression models. J Nucl Cardiol (2014) 21 (4 ):686–94. doi: 10.1007/s12350-014-9908-2 34 Delgado J Pereira A Villamor N López-Guillermo A Rozman C . Survival analysis in hematologic malignancies: recommendations for clinicians. Haematologica (2014) 99 (9 ):1410–20. doi: 10.3324/haematol.2013.100784 35 Zhang G Fan E Zhong Q Feng G Shuai Y Wu M . Identification and potential mechanisms of a 4-lncRNA signature that predicts prognosis in patients with laryngeal cancer. Hum Genomics (2019) 13 (1 ):36. doi: 10.1186/s40246-019-0230-6 31416476 36 Wu G Wang Q Zhu T Fu L Li Z Wu Y . Identification and validation of immune-related LncRNA prognostic signature for lung adenocarcinoma. Front Genet (2021) 12 :681277. doi: 10.3389/fgene.2021.681277 34306024 37 Maruyama R Suzuki H . Long non-coding RNA involvement in cancer. BMB Rep (2012) 45 (11 ):604–11. doi: 10.5483/BMBRep.2012.45.11.227 38 Jiang Q Wang Z Qi Q Li J Xin Y Qiu J . lncRNA SNHG26 promoted the growth, metastasis, and cisplatin resistance of tongue squamous cell carcinoma through PGK1/Akt/mTOR signal pathway. Mol Ther Oncolytics (2022) 24 :355–70. doi: 10.1016/j.omto.2021.12.021 39 Shao J Zhang B Kuai L Li Q . Integrated analysis of hypoxia-associated lncRNA signature to predict prognosis and immune microenvironment of lung adenocarcinoma patients. Bioengineered (2021) 12 (1 ):6186–200. doi: 10.1080/21655979.2021.1973874 40 López-Albaitero A Nayak JV Ogino T Machandia A Gooding W DeLeo AB . Role of antigen-processing machinery in the in vitro resistance of squamous cell carcinoma of the head and neck cells to recognition by CTL. J Immunol (2006) 176 (6 ):3402–9. doi: 10.4049/jimmunol.176.6.3402 41 Zhang Y Zhang Z . The history and advances in cancer immunotherapy: understanding the characteristics of tumor-infiltrating immune cells and their therapeutic implications. Cell Mol Immunol (2020) 17 (8 ):807–21. doi: 10.1038/s41423-020-0488-6 42 Chen YG Satpathy AT Chang HY . Gene regulation in the immune system by long non-coding RNAs. Nat Immunol (2017) 18 (9 ):962–72. doi: 10.1038/ni.3771 43 Peltanova B Raudenska M Masarik M . Effect of tumor microenvironment on pathogenesis of the head and neck squamous cell carcinoma: a systematic review. Mol Cancer (2019) 18 (1 ):63. doi: 10.1186/s12943-019-0983-5 30927923 44 Xia Y Rao L Yao H Wang Z Ning P Chen X . Engineering macrophages for cancer immunotherapy and drug delivery. Adv Mater (2020) 32 (40 ):e2002054. doi: 10.1002/adma.202002054 32856350 45 Hu Y He MY Zhu LF Yang CC Zhou ML Wang Q . Tumor-associated macrophages correlate with the clinicopathological features and poor outcomes via inducing epithelial to mesenchymal transition in oral squamous cell carcinoma. J Exp Clin Cancer Res (2016) 35 :12. doi: 10.1186/s13046-015-0281-z 26769084 46 Seminerio I Kindt N Descamps G Bellier J Lechien JR Mat Q . High infiltration of CD68+ macrophages is associated with poor prognoses of head and neck squamous cell carcinoma patients and is influenced by human papillomavirus. Oncotarget (2018) 9 (13 ):11046–59. doi: 10.18632/oncotarget.24306 47 Wuerdemann N Pütz K Eckel H Jain R Wittekindt C Huebbers CU . LAG-3, TIM-3 and VISTA expression on tumor-infiltrating lymphocytes in oropharyngeal squamous cell carcinoma-potential biomarkers for targeted therapy concepts. Int J Mol Sci (2020) 22 (1 ). doi: 10.3390/ijms22010379 48 Whiteside TL Jackson EK . Adenosine and prostaglandin e2 production by human inducible regulatory T cells in health and disease. Front Immunol (2013) 4 :212. doi: 10.3389/fimmu.2013.00212 23898333 49 Maggioni D Pignataro L Garavello W . T-Helper and T-regulatory cells modulation in head and neck squamous cell carcinoma. Oncoimmunology (2017) 6 (7 ):e1325066. doi: 10.1080/2162402X.2017.1325066 28811959 50 Strauss L Bergmann C Gooding W Johnson JT Whiteside TL . The frequency and suppressor function of CD4+CD25highFoxp3+ T cells in the circulation of patients with squamous cell carcinoma of the head and neck. Clin Cancer Res (2007) 13 (21 ):6301–11. doi: 10.1158/1078-0432.CCR-07-1403 51 Liang YJ Liu HC Su YX Zhang TH Chu M Liang LZ . Foxp3 expressed by tongue squamous cell carcinoma cells correlates with clinicopathologic features and overall survival in tongue squamous cell carcinoma patients. Oral Oncol (2011) 47 (7 ):566–70. doi: 10.1016/j.oraloncology.2011.04.017 52 Bron L Jandus C Andrejevic-Blant S Speiser DE Monnier P Romero P . Prognostic value of arginase-II expression and regulatory T-cell infiltration in head and neck squamous cell carcinoma. Int J Cancer (2013) 132 (3 ):E85–93. doi: 10.1002/ijc.27728 53 Lu X Kang Y . Hypoxia and hypoxia-inducible factors: master regulators of metastasis. Clin Cancer Res (2010) 16 (24 ):5928–35. doi: 10.1158/1078-0432.CCR-10-1360 54 Wicks EE Semenza GL . Hypoxia-inducible factors: cancer progression and clinical translation. J Clin Invest (2022) 132 (11 ). doi: 10.1172/JCI159839 55 Choudhry H Harris AL . Advances in hypoxia-inducible factor biology. Cell Metab (2018) 27 (2 ):281–98. doi: 10.1016/j.cmet.2017.10.005 56 Niu Y Bao L Chen Y Wang C Luo M Zhang B . HIF2-induced long non-coding RNA RAB11B-AS1 promotes hypoxia-mediated angiogenesis and breast cancer metastasis. Cancer Res (2020) 80 (5 ):964–75. doi: 10.1158/0008-5472.CAN-19-1532 57 Zhou L Jiang J Huang Z Jin P Peng L Luo M . Hypoxia-induced lncRNA STEAP3-AS1 activates wnt/β-catenin signaling to promote colorectal cancer progression by preventing m(6)A-mediated degradation of STEAP3 mRNA. Mol Cancer (2022) 21 (1 ):168. doi: 10.1186/s12943-022-01638-1 35986274 58 Senovilla L Vacchelli E Galon J Adjemian S Eggermont A Fridman WH . Trial watch: prognostic and predictive value of the immune infiltrate in cancer. Oncoimmunology (2012) 1 (8 ):1323–43. doi: 10.4161/onci.22009