A Novel Immune-Related lncRNA Signature predict the Immune Landscape and Prognosis for Head and Neck Squamous Cell Carcinoma

A Novel Immune-Related lncRNA Signature predict the Immune Landscape and Prognosis for Head and Neck Squamous Cell Carcinoma | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Primary research A Novel Immune-Related lncRNA Signature predict the Immune Landscape and Prognosis for Head and Neck Squamous Cell Carcinoma Zhimou Cai, Lin Chen, Yu Lin, Wenbin Lei This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-214493/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background: The tumor immune microenvironment is known to play an important role in head and neck squamous cell carcinomas (HNSCC). Reliable prognostic signatures that could accurately predict immune landscape and survival in HNSCC patients are vitally needed to promote a better individualized and effective treatment. Methods: HNSCC transcriptome data and clinical data in The Cancer Genome Atlas (TCGA) were embedded in our study. The differentially expressed irlncRNAs were identified by differential co-expression analysis, and recognized differently expressed irlncRNA (DEirlncRNA) pairs using univariate analysis. Cox and lasso regression analysis was used to identify DEirlncRNA pairs related to overall survival (OS) and build the prediction model. Then, we compared the areas under curve (AUC) with other published lncRNA signatures, counted the akaike information criterion (AIC) values of 3-year receiver operating characteristic curve, and identified the cut-off point to set up an optimal model for distinguishing the high- or low- risk groups among patients with HNSCC. Receiver operating characteristic (ROC)curves and Kaplan–Meier plot curves were used to validate the prediction model. Besides, We then reevaluated them from the viewpoints of clinical factor, tumor-infiltrating immune cells, chemotherapeutics efficacy, and immunosuppressed biomarkers. Results: We built a risk score model based on 18 DEirlncRNA pairs. The risk model is closely related to the OS of HNSCC patients. The hazard ratio (HR) is 1.376 [95% CI (confidence interval) 1.302-1.453] and log-rank P -value < 0.0001. Compared with two recently published lncRNA signatures, DEirLncRNA pairs signature has higher AUC score which showed the better prognostic performance. Additionally, the signature score showed a positive correlation with aggressive outcomes of HNSCC, such as low immunity score, significantly reduced CD8+ T cell infiltration and lowly expressed immunosuppressed biomarkers. However, high-risk groups of patients may have high chemotherapeutics sensitivity. Conclusions: The signature established by paring irlncRNA regardless of expression levels showed a promising clinical prediction value and revealed tumor immune microenvironment in HNSCC patients which might help in distinguishing those who could benefit from anti-tumor immunotherapy. Cancer Biology Head and neck squamous cell carcinoma Long non-coding RNAs Tumor-infiltrating immune cells Prognosis TCGA Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Figure 9 Background Head and neck tumors are the sixth most common malignancy in the world, head and neck squamous cell carcinoma (HNSCC) is the most common histological type [ 1 ]. HNSCC encompasses a heterogeneous group of malignancies that arise in the oral cavity, larynx and pharynx and are mainly associated with tobacco, alcohol consumption and human papillomavirus (HPV) infection [ 2 – 4 ]. However, due to the early symptoms of HNSCC are not obvious, the 5-year survival rate is still below 50%. Furthermore, HNSCC is characterized by high lymph node metastasis rate, local recurrence and metastasis can reduce the survival rate to 35% [ 5 ]. Currently, HNSCC is estimated to have an annual incidence of 900,000 and causes 450,000 death every year [ 1 ]. HNSCC is a disease characterized by profound immunosuppression [ 6 ]. immune checkpoint inhibitors (ICIs) therapy is becoming a new hope for HNSCC treatment. In recent years, considerable progress has been made in ICIs therapy for HNSCC. However, the response rate of recurrent or metastatic HNSCC for PD-1/PD-L1 inhibitors was only found to be 13.3–22% in previous clinical trials [ 6 ]. Therefore, biomarkers which could effectively predict the efficacy of ICIs are crucial for patient selection. Generally, PD-L1 expression level is used as the representative marker for predicting the efficacy of ICIs. However, for most tumors, only the detection of PD-L1 alone cannot screen the dominant population [ 7 – 9 ]. Long non-coding RNAs (lncRNAs) refer to functional RNA molecules with a length more than 200 nucleotides and have minimal or no function to encode proteins [ 10 ]. It has become increasingly recognized as lncRNAs can interact with DNA, RNA and proteins which regulate gene expression at the transcription and post-transcription levels and is associated with various types of cancer [ 10 , 11 ]. Recently, several lncRNA signatures have been established for prognosis prediction in some cancers, including HNSCC. Mao et al. [ 12 ] identified an eight lncRNA signature is an independent prognostic factor of HNSCC patients. Bao et al. [ 13 ] identify genome instability-associated lncRNAs and constructed a signature to evaluate the prognosis of patients with breast cancer. It is know that the diagnosis, evaluation, and treatment of cancer is closely related to the tumor immune microenvironment, especially to infiltrating immune cells [ 14 , 15 , 16 ]. Now research confirms that, lncRNAs are critical regulators of gene expression in the immune system. Research shows that lncRNAs direct the expression of genes related to the activation of immune cells that directing the development of diverse immune cells so resulting in tumor immune cell infiltration and altering the immune microenvironment [ 17 , 18 ]. Therefore, in this study, we tried to identify immune-related lncRNAs (irlncRNAs) to investigate the relationship between irlncRNAs and immune infiltration, construct a powerful prognostic model for risk assessment of HNSCC and predict the response of patients to immunotherapy. Different from simple genes in the past, we utilized a novel modeling algorithm, paring, and iteration to construct an irlncRNA signature, which adopted two-biomarker combinations will help improve the accuracy of the cancer prognostic model. Methods Data collectin HNSCC RNAseq expression profile data were downloaded from The Cancer Genome Atlas (TCGA) database ( https://portal.gdc.cancer.gov/ ). The database contained 501 cases tumor tissue sequencing information with HNSCC patients and 44 cases of normal tissue. Clinical and follow-up data of patients with HNSCC were retrieved from the HNSC project of TCGA, valid data were extracted after eliminating the data with follow-up time of 30 days and repeated data. “Ensembl_Stable_id” is converted to "Gene Symbol" in the RNAseq expression profile according to the GTF file which downloaded from Ensembl ( http://asia.ensembl.org ), and annotation to distinguish the mRNAs and lncRNAs for subsequent analysis. The list of immune-related genes (ir-genes) was downloaded from the ImmPort database ( http://www.immport.org ), and correlation analysis was conducted between ir-genes and all lncRNAs to screen for irlncRNAs. Immune gene correlation coefficients greater than 0.4 or less than − 0.4 and p value less than 0.001 were considered as irlncRNAs. The differentially expressed irlncRNAs (DEirlncRNAs) were identified through the “limma” R package, FDR 1.5 were defined as DEirlncRNAs. Pairing DEirlncRNAs The DEirlncRNAs were cyclically singly paired, and construct a 0-or-1 matrix. If the expression level of the former lncRNA was higher than that of the latter, it was defined as 1; otherwise, it was defined as 0. Then the constructed 0-or-1 matrix is further screened. When the amounts of DEirlncRNA pairs of which expression quantity was 0 or 1 accounted for more than 20% and less than 80% of total pairs, it was considered a valid match. Building the prognostic pairing DEirlncRNAs-based signature Univariate Cox analysis was performed by “survival” R package to screen DEirlncRNA pairs with prognostic values and P < 0.01 DEirlncRNA pairs is considered to be significantly related to the prognosis of HNSCC patients. Next, Lasso regression analysis performed by “glmnet” and “survival” R package to select the minimum error value to eliminate over-fitting data, the optimal penalty parameter "lambda" value was calculated via a 1,000 cross-validations. Then multivariate analyses were performed using the Cox regression model to identify the prognostic-associated DEirlncRNA pairs factors in HNSCC. Based on the coefficients from the multivariate regression analysis and the status of DEirlncRNA pairs, we constructed a DEirlncRNA pairs signature(DEirLnc-PSig) for outcome prediction as follows: where β i indicates the coefficient for each DEirlncRNA pairs and Pi indicates the DEirlncRNA pairs status of 0-or-1 matrix. Validation of the constructed risk model A time-dependent receiver operating characteristic (ROC) curve and area under curve(AUC) was estimated using the “ survivalROC” R package to predict the accuracy of DEirLnc-PSig for HNSCC patients. The 1-, 2-, and 3-year ROC curves and AUC of the model were plotted. The AIC values of every points of the 3-year ROC curve were evaluated to identify the maximum inflection point, which was considered as the cut-off point to classify patients into the high-risk group with high DEirLnc-PSig or low-risk group with low DEirLnc-PSig. The Kaplan–Meier method was used to construct survival curves, and the log-rank test was used to assess the difference in survival between the high- and low-risk groups with a significant level of 0.05. The R packages utilized in these steps included “survival” and “survminer”. To verify the clinical application value of the DEirLnc-PSig, multivariate Cox regression and stratified analysis were used to assess the independence of DEirLnc-PSig from other key clinical factors. A forest map was used to demonstrate the results. The R packages utilized in these operations were “survival” and “survminer”. In addition, we compared our defined DEirLnc-PSig with two previously published lncRNA signatures to further compare the prognostic model performance. Investigation of tumor-infiltrating immune cells To explore the differences in immune cell infiltration between different risk groups, we used the TIMER, CIBERSORT, XCELL, QUANTISEQ, MCPcounter, EPIC, and CIBERSORT methods to assess the samples from the TCGA project of the HNSCC dataset. Differences in the number of different types of immune infiltrating cell between high- and low-risk groups of the constructed model were analyzed using the Wilcoxon's signed-rank test and the results are shown in a box chart. The correlation between the risk scores and immune cell infiltration were evaluated using Spearman's correlation analysis. The correlation coefficients of the results were shown in a lollipop diagram. The significance threshold was set as p < 0.05. The R packages utilized in these operations were “limma”, “scales”, “ggplot2”, “ggtext” and “ggpubr”. Exploration of the significance of the model in the clinical treatment We calculated the half inhibitory centration (IC 50 ) of common administrating chemotherapeutic drugs in the TCGA project of the HNSCC dataset by " pRRophetic" R package to evaluate the model in the clinic for head and neck carcinoma treatment. Antitumor drugs such as cisplatin, docetaxeland, methotrexate and paclitaxel are recommended for HNSCC treatment by NCCN guidelines. The difference in the IC 50 between the high- and low-risk groups was compared by Wilcoxon’s signed-rank test and the results are shown as box drawings. Analyses of the immunosuppressive molecules expressing related to ICIs The relationship between the model and the expression level of genes related to ICIs were performed by "limma" and “ggpubr” R package and violin plot visualization. Results Identification of irlncRNAs and differentially expressed irlncRNAs A study process flow chart is given in Fig. 1 . A total of 805 irlncRNAs were identified by co-expression analysis of ir-genes and lncRNAs. The differential expression analysis was subsequentially carried out using “limma” R package, 132 DEirLncRNAs (75 upregulated and 16 downregulated) were identified between the tumor and normal tissues (Fig. 2 ). Establishment of DEirlncRNA pairs and a DEirlncRNA pairs signature Using an iteration loop and a 0-or-1 matrix screening among 132 DEirlncRNAs and 6,686 valid DEirlncRNA pairs were identified. By conducting the univariate Cox regression analysis, 500 survival-related DEirlncRNA pairs were selected (p < 0.01). To further avoid result overfitting, the Lasso regression analysis was employed to screen out 33 DEirlncRNA pairs showing the highest correlation with survival (Fig. 3 a, b). Those 33 DEirlncRNA pairs were further subjected to multivariable Cox regression analysis, finalizing the most significant 18 DEirlncRNA pairs was developed as a DEirlncRNA pairs signature for HNSCC patients by the stepwise elimination method and obtaining corresponding coefficients (Fig. 3 c). We then applied the 1-, 2-, and 3-year ROC curves to verify the predictive capacity of the DEirLnc-Psig; the AUCs of the ROC were almost all ~ 0.75, suggesting a satisfactory predictive performance (Fig. 4 a). In addition, we compared the 3-year ROC curve with other clinical features and the superiority of DEirLnc-PSig have been further demonstrated (Fig. 4 b). Clinical evaluation by risk assessment model and independence of the DEirLnc-PSig from other clinical factors We recognized the maximum inflection point as the cut-off point on the 3-year ROC curve using the Akaike information criterion (AIC) values and according to the cut-off point (Fig. 4 c, d), HNSCC patients from the TCGA cohort was divided into a low-risk group(n = 255) and a high-risk group(n = 235). We then compared the survival of low-risk group with the high-risk group and found that the DEirLnc-PSig possessed powerful capacity to predict the prognosis of HNSCC patients between the low- and high-risk group, in which the low-risk group are significantly better than the high-risk group (p < 0.001, log-rank test) (Fig. 5 a, b, c). To evaluate whether the prognostic value of the DEirLnc-PSig was independent of common clinical variables, univariate and multivariate Cox regression analyses were performed on age, gender, pathologic grade, clinical stage and our DEirLnc-PSig-based prognostic risk score model (Fig. 6 a, b). The results demonstrated that risk score showed statistical differences by univariate and multivariate Cox regression analysis which presented as an independent prognostic predictor. In addition to the DEirLnc-PSig, there were other two clinical factors, age and clinical stage, that were observed to be significant in the multivariate analysis. Therefore, stratification analysis was performed to determine whether the DEirLnc-PSig possessed a prognostic value that was independent of the age and clinical stage. Patients in the TCGA set were stratified into a young-patient group ( 65 years, n = 253). Using the DEirLnc-PSig, patients in each age group could be further divided into high-risk or low-risk group. There was a significant difference in OS between the high- and low-risk groups in the young-patient group (p < 0.001, log-rank test) as was in the old-patient group (p < 0.001, log-rank test) (Fig. 6 c, d). Next, all HNSCC patients were then stratified into the early-stage group (stage I and II, n = 70) and the late-stage group (stage III and IV, n = 307). Using the DEirLnc-PSig, patients in each stage group could be further divided into high-risk or low-risk group. There was a significant difference in OS between the high- and low-risk groups in the early-stage group (p = 0.013, log-rank test) as was in the late-stage group (p < 0.001, log-rank test) (Fig. 6 e, f). These results indicated that the DEirLnc-PSig was an independent prognostic factor associated with OS in HNSCC patients. Performance comparison of the DEirlncRNA pairs signature with existing lncRNA-related signatures in survival prediction We further compared the prediction performance of the DEirLnc-PSig with two recently published lncRNA signatures: 3-lncRNA signature derived from Jiang’s study (hereinafter referred to as JiangLncSig) [ 19 ] and 5-lncRNA signature derived from Liu’s study (hereinafter referred to as LiuLncSig) [ 20 ] using the same TCGA patient cohort. As shown in Fig. 7 , the AUC at 3 years of OS for the DEirLnc-PSig is 0.809, which is significantly higher than that of JiangLncSig (AUC = 0.632) and LiuLncSig (AUC = 0.5). These results demonstrated the better prognostic performance of the DEirLnc-PSig in predicting survival than two recently published lncRNA signatures. Estimation of tumor-infiltrating immune cells and immunosuppressive molecules with risk assessment model Since the model was established by irlncRNAs, we consequently investigated whether the model was related to the tumor immune microenvironment. Correlation analysis of immunocyte infiltration revealed that the high-risk group was more positively associated with CD4 + T cells, cancer associated fibroblast, monocyte and macrophages, whereas were negatively associated with CD8 + T cells and T cell regulatory (Tregs), as revealed by the Wilcoxon signed-rank test (see supplementary file, Table S1 and Figure S1). In addition, the immune scores in the low-risk group were significantly higher than those in the high-risk group (p = 0.0038) in the analysis with xCell, which suggesting that the tumor cells in the low-risk group had more immune cell infiltration than those in the high-risk group. A detailed Spearman correlation analysis was conducted, and the resulting diagram exhibited a lollipop shape, as shown in Fig. 8 a. Immunotherapy with ICIs has changed HNSCC treatment strategies dramatically, improve survival of HNSCC patients. We investigated whether the risk model was related to ICI-related biomarkers and discovered that high risk scores were positively correlated with low expression of PDCD1 (p < 0.001, Fig. 8 b), LAG3 (p < 0.01, Fig. 8 c) and CTLA4 (p < 0.001, Fig. 8 d). Analysis of the correlation between the risk model and chemotherapeutics As common therapy for HNSCC, we attempted to identify associations between risk and the efficacy of chemotherapeutics in treating head and neck squamous cell carcinoma in the TCGA project of the HNSCC dataset. We showed that a high risk score was associated with a lower IC 50 of chemotherapeutics such as docetaxel (p < 0.001), paclitaxel (p = 0.02), whereas it was associated with a higher IC 50 for methotrexate (p < 0.001), which indicated that the model acted as a potential predictor for chemosensitivity (Fig. 9 ). Discussion Despite the effectiveness of targeted therapies and immunotherapy agents in HNSCC, not all patients respond to treatment and development of resistance is almost universal after a time. In clinical practice, ICIs therapy does not appear to improve survival in HNSCC patients. Precision medicine approaches that seek to individualise therapy through the use of predictive biomarkers and stratification strategies, have been the key to improve therapeutic success in HNSCC [ 21 , 22 ]. In recent years, many studies have shown that lncRNA is an important and significant factor in the current and future prognosis of patients with HNSCC [ 22 ]. For this reason, many prognostic lncRNA expression signatures have been applied to diagnose and predict prognosis of HNSCC. However, it should be pointed out that HNSCC are a type of immunosuppressive disease [ 23 ]. In the case of HNSCC, a deeper understanding of the heterogeneity and complex of the tumor immune microenvironment is essential to develop sensitive predictive biomarkers for providing information of tumor microenvironment and help to improve the efficacy of tumor immunotherapy. An ideal signature would not only to predict the prognosis of cancer patients, but also reflect the characteristics of the tumour and immune microenvironment that have an impact on disease progression and treatment response. Therefore, we established a prediction model of immune-related lncRNAs and evaluated the relationships between the risk score draw from DEirLnc-PSig, immune checkpoints, tumor-infiltrating immune cell abundances and the response to chemotherapy. Unlike most signatures, which are based on quantified transcript expression levels, we attempted to construct a reasonable model with two-lncRNA combinations and did not adopt their exact expression levels in the signature. Thus, pairs with higher or lower expression only had to be detected instead of examining specific expression values of every lncRNA. The results of this study showed that this novel model had a good predictive performance and achieved a better prediction performance with an AUC value of 0.809, which is higher than other two expression based lncRNA signature, risk score derived from the DEirLnc-PSig is reliability and independent for predicting HNSCC prognosis. In this study, some of the irlncRNAs in the process of modeling that have been already identified play an important role in malignant phenotypes of various cancer types, such as LINC01063 [ 24 , 25 ], HOXC-AS1 [ 26 , 27 ], LINC02454 [ 28 , 29 ] and SCAT1 [ 30 ]. Furthermore, previous reports have shown that LINC00460 [ 31 ], C5orf66-AS1 [ 32 ] and TMPO-AS1 [ 33 ] are associated with prognosis and treatment outcomes in HNSCC. Jiang et al. [ 31 ] revealed that LINC00460 promotes EMT in HNSCC by facilitating peroxiredoxin-1 into the nucleus. Lu et al. [ 32 ] first report that C5orf66-AS1 was able to prevent OSCC progression by inhibiting OSCC cell growth and metastasis via regulating CYC1 expression. Xing et al. [ 33 ] disclosed that a novel TMPO-AS1/miR-320a/SOX4 pathway associated with nasopharyngeal carcinoma(NPC) progression, suggesting that TMPO-AS1 may be a potential therapeutic target for NPC. However, there are also some lncRNAs that have been disclosed for the first time in this model which suggested that the proposed model can identify novel biomarkers for further research. Tumor-infiltrating immune cells which are integral to the tumor immune microenvironment (TIME) are recognized as highly important for prognosis prediction, while the responsiveness to immune checkpoint blockade may be subject to the features of TIME. To explore the relationship between risk scores and tumor-infiltrating immune cells, we used seven common acceptable methods to estimate the immune infiltrating cell, including xCell [ 34 , 35 ], TIMER [ 36 , 37 ], QUANTISEQ [ 38 , 39 ], MCPcounter [ 40 ], EPIC [ 41 ], CIBERSORT-abs [ 42 ], and CIBERSORT [ 43 ]. Multiple studies have indicated that greater infiltration of CD8 T cells and high immune scores correlated with better response to standard chemotherapy and to immune checkpoint blockade therapy in tumor patients [ 44 , 45 ]. By integrating analyzing, our results revealed that the high-risk group identified by DEirLnc-PSig was negatively associated with CD8 + T cell and immune scores. Correspondingly, the high-risk group was also negatively correlated with checkpoint-related biomarkers such as PDCD1, LAG3 and CTLA4. Our results demonstrated that patients in high-risk groups identified by our model may not benefit from immunotherapy. But these conclusions need to be verified through further studies especially supported by clinical trials. On the other hand, we checked several more drugs like cisplatin, paclitaxel, doxorubicin and methotrexate, we observed from the estimated IC 50 of these chemotherapeutic drugs, that the high-risk group could be more sensitive to chemotherapies than those in low-risk group. It would be more helpful to identify patients at high clinical risk and give sensitive chemotherapy rather than immunotherapy. However, more benchwork and clinical studies are needed to further validation. Conclusions In summary, our study revealed a comprehensive landscape of the tumor immune microenvironment in HNSCC, and provides a promising prognostic signature for HNSCC. Patients with the low risk scores could benefit more from ICIs therapy while patients with the high risk scores may benefit more from chemotherapy. However, further investigations are needed to verify the accuracy in estimating prognoses and to test its clinical utility in patient management. Abbreviations HNSCC: Head and neck squamous cell carcinomas; TCGA: The Cancer Genome Atlas; DEirlncRNA: Differently expressed irlncRNA; OS: Overall survival; AUC: Areas under curve; AIC: Akaike information criterion; ROC: Receiver operating characteristic; HR: hazard ratio; HPV: human papillomavirus; ICIs: immune checkpoint inhibitors; lncRNAs: Long non-coding RNAs; irlncRNAs: Immune-related lncRNAs; ir-genes: immune-related genes; DEirlncRNAs: Differentially expressed irlncRNAs; IC 50: Inhibitory centration; Tregs: T cell regulatory; TIME: Tumor immune microenvironment; Declarations Ethics approval and consent to participate Not applicable. Consent for publication Not applicable. Availability of data and materials The datasets analyzed during the current study are available in the TCGA repository, https://portal.gdc.cancer.gov/. Competing interests The authors report no competing interest in this study. Funding This work was supported by the National Key R&D Program of China (No 2020YFC1316903), the 5010 Clinical Research Program of Sun Yat-sen University (No. 2017004) and the National Natural Science Foundation of China (No. 81972528). Authors’ contributions ZMC designed the research, collected the data, analyzed the date and wrote the manuscript. LC took part in designing the research. YL collected the data, analyzed the date. WBL took part in designing the research and solved the disagreements between ZMC and LC. All authors read and approved the final manuscript. Acknowledgements None. Author information ZMC, LC and YL contributed equally to this work. References Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018 Nov;68(6):394–424. Katiyar SK. Emerging Phytochemicals for the Prevention and Treatment of Head and Neck Cancer. Molecules. 2016 Nov 24;21(12):1610. Fleming JC, Woo J, Moutasim K, Mellone M, Frampton SJ, Mead A, et al. HPV, tumour metabolism and novel target identification in head and neck squamous cell carcinoma. Br J Cancer. 2019 Feb;120(3):356–67. Nelson HH, Pawlita M, Michaud DS, McClean M, Langevin SM, Eliot MN, et al. Immune response to HPV16 E6 and E7 proteins and patient outcomes in head and neck cancer. JAMA Oncol. 2016;3:178–85. Posner M, Vermorken JB. Induction therapy in the modern era of combined-modality therapy for locally advanced head and neck cancer. Semin Oncol. 2008 Jun;35(3):221–8. Romano E, Romero P. The therapeutic promise of disrupting the PD-1/PD-L1 immune checkpoint in cancer: Unleashing the CD8 T cell mediated anti-tumor activity results in significant, unprecedented clinical efficacy in various solid tumors. J Immunother Cancer. 2015;3:15. Wang SM, Jiang B, Deng Y, Huang SL, Fang MZ, Wang Y. Clinical significance of MLH1/MSH2 for stage II/III sporadic colorectal cancer. World J Gastrointest Oncol. 2019 Nov;15(11):1065–80. 11(. Schepisi G, Brighi N, Cursano MC, Gurioli G, Ravaglia G, Altavilla A, et al. Inflammatory Biomarkers as Predictors of Response to Immunotherapy in Urological Tumors. J Oncol. 2019 Sep;19:2019:7317964. Liu Y, Zhou S, Du Y, Sun L, Jiang H, Zhang B, et al. Efficacy and safety of programmed death 1 inhibitors in patients with advanced non-small cell lung cancer: a meta-analysis. Cancer Manag Res. 2019 May;21:11:4619–30. Kopp F, Mendell JT. Functional Classification and Experimental Dissection of Long Noncoding RNAs. Cell. 2018 Jan 25;172(3):393–407. Schmitt AM, Chang HY. Long Noncoding RNAs: At the Intersection of Cancer and Chromatin Biology. Cold Spring Harb Perspect Med. 2017 Jul;5(7):a026492. 7(. Mao R, Chen Y, Xiong L, Liu Y, Zhang T. Identification of a nomogram based on an 8-lncRNA signature as a novel diagnostic biomarker for head and neck squamous cell carcinoma. Aging (Albany NY). 2020 Oct 22;12(20):20778–20800. Bao S, Zhao H, Yuan J, Fan D, Zhang Z, Su J, Zhou M. Computational identification of mutator-derived lncRNA signatures of genome instability for improving the clinical outcome of cancers: a case study in breast cancer. Brief Bioinform. 2020 Sep 25;21(5):1742–1755. Martinez LM, Robila V, Clark NM, Du W, Idowu MO, Rutkowski MR, et al. Regulatory T cells control the switch from in situ to invasive breast cancer. Front Immunol. 2019;10:1942. Peng D, Wang L, Li H, Cai C, Tan Y, Xu B, et al. An immune infiltration signature to predict the overall survival of patients with colon cancer. IUBMB Life. 2019 Nov;71(11):1760–70. Li AL, Zhu YM, Gao LQ, Wei SY, Wang MT, Ma Q, Zheng YY, et al. Exploration of the Immune-Related Signatures and Immune Infiltration Analysis in Melanoma. Anal Cell Pathol (Amst). 2021 Jan;2021:16. 4743971. Atianand MK, Caffrey DR, Fitzgerald KA. Immunobiology of Long Noncoding RNAs. Annu Rev Immunol. 2017 Apr 26;35:177–198. Chen YG, Satpathy AT, Chang HY. Gene regulation in the immune system by long noncoding RNAs. Nat Immunol. 2017 Aug;22(9):962–72. 18(. Jiang H, Ma B, Xu W, Luo Y, Wang X, Wen S, et al. A Novel Three-lncRNA Signature Predicts the Overall Survival of HNSCC Patients. Ann Surg Oncol. 2020 Oct 23. Liu G, Zheng J, Zhuang L, Lv Y, Zhu G, Pi L, et al. A Prognostic 5-lncRNA Expression Signature for Head and Neck Squamous Cell Carcinoma. Sci Rep. 2018 Oct;15(1):15250. 8(. Lee TW, Lai A, Harms JK, Singleton DC, Dickson BD, Macann AMJ, et al. Patient-Derived Xenograft and Organoid Models for Precision Medicine Targeting of the Tumour Microenvironment in Head and Neck Cancer. Cancers (Basel). 2020 Dec 12;12(12):3743. Wang Y, Wang S, Ren Y, Zhou X. The Role of lncRNA Crosstalk in Leading Cancer Metastasis of Head and Neck Squamous Cell Carcinoma. Front Oncol. 2020 Oct 2;10:561833. Lee TW, Lai A, Harms JK, Singleton DC, Dickson BD, Macann AMJ, et al. Patient-Derived Xenograft and Organoid Models for Precision Medicine Targeting of the Tumour Microenvironment in Head and Neck Cancer. Cancers (Basel). 2020 Dec 12;12(12):3743. Deng X, Bi Q, Chen S, Chen X, Li S, Zhong Z, et al. Identification of a Five-Autophagy-Related-lncRNA Signature as a Novel Prognostic Biomarker for Hepatocellular Carcinoma. Front Mol Biosci. 2021 Jan;11:7:611626. Cheng L, Han T, Zhang Z, Yi P, Zhang C, Zhang S, et al. Identification and Validation of Six Autophagy-related Long Non-coding RNAs as Prognostic Signature in Colorectal Cancer. Int J Med Sci. 2021 Jan 1;18(1):88–98. Zhou C, An N, Cao C, Wang G. lncRNA. HOXC-AS1 promotes gastric cancer via binding eIF4AIII by activating Wnt/β-catenin signaling. J Gene Med. 2020 Sep;22(9):e3202. Dong Y, Li X, Lin Z, Zou W, Liu Y, Qian H, et al. HOXC-AS1-MYC regulatory loop contributes to the growth and metastasis in gastric cancer. J Exp Clin Cancer Res. 2019 Dec;23(1):502. 38(. Cai WY, Chen X, Chen LP, Li Q, Du XJ, Zhou YY. Role of differentially expressed genes and long non-coding RNAs in papillary thyroid carcinoma diagnosis, progression, and prognosis. J Cell Biochem. 2018 Nov;119(10):8249–59. Tan J, Liu L, Zuo Z, Song B, Cai T, Ding D, et al. Overexpression of novel long intergenic noncoding RNA LINC02454 is associated with a poor prognosis in papillary thyroid cancer. Oncol Rep. 2020 Oct;44(4):1489–501. Zhang C, Zhang Z, Zhang G, Xue L, Yang H, Luo Y, et al. A three-lncRNA signature of pretreatment biopsies predicts pathological response and outcome in esophageal squamous cell carcinoma with neoadjuvant chemoradiotherapy. Clin Transl Med. 2020 Aug;10(4):e156. Jiang Y, Cao W, Wu K, Qin X, Wang X, Li Y, et al. LncRNA LINC00460 promotes EMT in head and neck squamous cell carcinoma by facilitating peroxiredoxin-1 into the nucleus. J Exp Clin Cancer Res. 2019 Aug 20;38(1):365. Lu T, Liu H, You G. Long non-coding RNA C5orf66-AS1 prevents oral squamous cell carcinoma through inhibiting cell growth and metastasis. Int J Mol Med. 2018 Dec;42(6):3291–9. Xing B, Qiao XF, Qiu YH, Li X. TMPO-AS1 Regulates the Aggressiveness-Associated Traits of Nasopharyngeal Carcinoma Cells Through Sponging miR-320a. Cancer Manag Res. 2021 Jan 15;13:415–425. Aran D. Cell-Type Enrichment Analysis of Bulk Transcriptomes Using xCell. Methods Mol Biol. 2020;2120:263–76. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017 Nov 15;18(1):220. Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 2017 Nov 1;77(21):e108-e110. Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020 Jul;48(W1)(2):W509–14. Finotello F, Mayer C, Plattner C, Laschober G, Rieder D, Hackl H, et al. Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data. Genome Med. 2019 May 24;11(1):34. doi: 10.1186/s13073-019-0638-6 . Erratum in: Genome Med. 2019 Jul 29;11(1):50. PMID: 31126321; PMCID: PMC6534875. Plattner C, Finotello F, Rieder D. Deconvoluting tumor-infiltrating immune cells from RNA-seq data using quanTIsEq. Methods Enzymol. 2020;636:261–85. Becht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016 Oct 20;17(1):218. doi: 10.1186/s13059-016-1070-5 . Erratum in: Genome Biol. 2016 Dec 1;17 (1):249. Racle J, Gfeller D. EPIC: A Tool to Estimate the Proportions of Different Cell Types from Bulk Gene Expression Data. Methods Mol Biol. 2020;2120:233–48. Wu YY, Wang SH, Wu CH, Yen LC, Lai HF, Ho C, et al. In silico immune infiltration profiling combined with functional enrichment analysis reveals a potential role for naïve B cells as a trigger for severe immune responses in the lungs of COVID-19 patients. PLoS One. 2020 Dec 2;15(12):e0242900. Chen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling Tumor Infiltrating Immune Cells with CIBERSORT. Methods Mol Biol. 2018;1711:243–59. Chen Y, Chen H, Mao B, Zhou Y, Shi X, Tang L, et al. Transcriptional Characterization Of The Tumor Immune Microenvironment And Its Prognostic Value For Locally Advanced Lung Adenocarcinoma In A Chinese Population. Cancer Manag Res. 2019 Oct 31;11:9165–9173. Ge PL, Li SF, Wang WW, Li CB, Fu YB, Feng ZK, et al. Prognostic values of immune scores and immune microenvironment-related genes for hepatocellular carcinoma. Aging. 2020 Mar;25(6):5479–99. 12(. Supplementary Files SupplementalInformation.docx Additional file 1: Table S1. The detail comparison results of correlation ship between tumour infiltrating immune cells and risk sore. Figure S1. The representative results of the evaluation of tumor infiltrating immune cells with risk assessment model. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-214493","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Primary research","associatedPublications":[],"authors":[{"id":11848305,"identity":"798eecaa-773a-480c-9048-2650accdf651","order_by":0,"name":"Zhimou Cai","email":"","orcid":"","institution":"Sun Yat-sen University First Affiliated Hospital","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Zhimou","middleName":"","lastName":"Cai","suffix":""},{"id":11848306,"identity":"1dd086c3-1843-421d-aa80-d477e0cb6570","order_by":1,"name":"Lin Chen","email":"","orcid":"","institution":"Sun Yat-sen University First Affiliated Hospital","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Lin","middleName":"","lastName":"Chen","suffix":""},{"id":11848307,"identity":"de81c4b6-f458-4cf9-896b-dc1fc68857db","order_by":2,"name":"Yu Lin","email":"","orcid":"","institution":"Sun Yat-sen University First Affiliated Hospital","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Yu","middleName":"","lastName":"Lin","suffix":""},{"id":11848308,"identity":"dac178dd-1db3-4b4f-a9be-1b1c0336134f","order_by":3,"name":"Wenbin Lei","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA8UlEQVRIiWNgGAWjYBACPgbGBmaGCgYGNna4WAJ+LWxgLWeADGbitTAwMDO2gUiitbA3t0kXztsmz8fMwPyZ589hBn72HAOGnzvwaOE52CY9c9ttwzZmBjZp3rbDDJI9bwwYe8/g0SKR2Habd9ttRpAWZt6GwwwGN3IMIE7FpUX+IVDLnNv2bTCH2RPUIsEI1NJwOxGohUGahw1oiwQhLTyJ7b9nHLud3AZUJjm3LZ1H4syzgoO9eLTwsx9/bFxQc9t2fnvz4Q9v/ljL8bcnb3zwE48WJMDYwMTDwMADYh4gSgNY0w+ilY6CUTAKRsFIAgDG90gdiE1SQgAAAABJRU5ErkJggg==","orcid":"https://orcid.org/0000-0002-9720-1997","institution":"Sun Yat-sen University First Affiliated Hospital","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Wenbin","middleName":"","lastName":"Lei","suffix":""}],"badges":[],"createdAt":"2021-02-06 07:03:54","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-214493/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-214493/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":6224456,"identity":"7fa89ee8-d5e1-4116-a069-01c370333d33","added_by":"auto","created_at":"2021-02-22 21:06:52","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":31512,"visible":true,"origin":"","legend":"Flow chart of this study","description":"","filename":"OnlineFig.1.png","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/dbfce809381544f0e40c4485.png"},{"id":6224457,"identity":"887f6f5f-8247-4e4d-a96a-4b12ef4fcb3c","added_by":"auto","created_at":"2021-02-22 21:06:52","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":104199,"visible":true,"origin":"","legend":"The heatmap and volcano plot of differently expressed immune-related IncRNAs. a The heatmap. b The volcano plot","description":"","filename":"OnlineFig.2.png","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/99f2af624c5f249c65628202.png"},{"id":6224248,"identity":"b96f5df3-3435-4077-b9c0-5713daacafb0","added_by":"auto","created_at":"2021-02-22 21:03:52","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":40073,"visible":true,"origin":"","legend":"Diagrams of lasso regression analysis and multivariable Cox regression analysis. a lasso regression analysis. b Penalty coefficient was utilized to minimize the mean square error of the models. c A forest graph of the 18 DEirlncRNA pairs identified by multivariable Cox regression analysis","description":"","filename":"OnlineFig.3.png","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/96d177fa156ec6641ab31a19.png"},{"id":6224461,"identity":"afd4a08c-9b20-4a9f-89c0-df992eef9009","added_by":"auto","created_at":"2021-02-22 21:06:53","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":32656,"visible":true,"origin":"","legend":"The results of ROC curves. a The 1-, 2-, and 3-year ROC of the DEirLnc-Psig. b A comparison of 3-year ROC curves with other common clinical characteristics showed the superiority of the risk score. c-d The maximum inflection point obtained by AIC","description":"","filename":"OnlineFig.4.png","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/5c93f0fb92ac649af1c21775.png"},{"id":6224247,"identity":"50512545-b6ed-4459-bd9d-c4c5b19b4b6c","added_by":"auto","created_at":"2021-02-22 21:03:52","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":100340,"visible":true,"origin":"","legend":"Survival prediction of HNSCC patients with the risk score model. a-b Patients were scored and grouped into the low-risk group (green) and high-risk group (red). Then scatter diagrams of patient’s survival were plotted based on the risk scores from low to high, with green referring to survival patients and red referring to deaths. c Kaplan–Meier survival curves","description":"","filename":"Fig.5.png","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/ccdac1d2032479dfe7dc7e8c.png"},{"id":6224249,"identity":"bc236c90-52a2-4acd-b7d2-16e677ce03df","added_by":"auto","created_at":"2021-02-22 21:03:52","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":52958,"visible":true,"origin":"","legend":"Univariate and multivariate Cox analysis to identify the risk factors and stratification analyses by age and stage. a Univariate. b Mutivariate. c-d KM curve analysis of OS in high- and low-risk groups for young patients and old patients. e-f KM curve analysis of OS in high- and low-risk groups for early-stage patients and late-stage patients. Statistical analysis was performed using the log-rank test.","description":"","filename":"OnlineFig.6.png","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/ff3acdc29e9c89179e6cb1c3.png"},{"id":6224496,"identity":"8ec1f52b-a6eb-4445-b642-236d18e60756","added_by":"auto","created_at":"2021-02-22 21:09:53","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":18227,"visible":true,"origin":"","legend":"The ROC analysis at 3-year of OS for DEirLnc-Psig, JiangLncSig, and LiuLncSig\n\n","description":"","filename":"OnlineFig.7.png","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/be7ca3080d0639dedba68579.png"},{"id":6224254,"identity":"32727b3e-23b7-4b99-9642-8ecefe526de7","added_by":"auto","created_at":"2021-02-22 21:03:53","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":55434,"visible":true,"origin":"","legend":"Estimation of tumor-infiltrating cells and immunosuppressed molecules by the risk assessment model. a Lollipop diagram of the correlation between the risk scores and immune cell infiltration. b-d Violin plot showed the expression of genes related to ICIs between low- and high-risk group","description":"","filename":"OnlineFig.8.png","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/8946be4fd07083957af37ce0.png"},{"id":6224545,"identity":"7dda65ca-8f7c-409e-978c-f64947c96bb8","added_by":"auto","created_at":"2021-02-22 21:12:53","extension":"png","order_by":9,"title":"Figure 9","display":"","copyAsset":false,"role":"figure","size":16288,"visible":true,"origin":"","legend":"Boxplots for the difference in the IC50 or chemotherapeutic drugs such as cisplatin, docetaxeland, methotrexate and paclitaxel between the high- and low-risk groups","description":"","filename":"OnlineFig.9.png","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/6fc12d7a277df91daa5c5604.png"},{"id":13669218,"identity":"6b90eb78-bc98-432b-b75c-b2402d8462d4","added_by":"auto","created_at":"2021-09-17 10:59:42","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1044577,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/88bd8502-369f-4dac-ba6a-6a539ac06b48.pdf"},{"id":6224458,"identity":"e69b480d-0db4-4424-9320-6edff9a6db77","added_by":"auto","created_at":"2021-02-22 21:06:53","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":416836,"visible":true,"origin":"","legend":"Additional file 1: Table S1. The detail comparison results of correlation ship between tumour infiltrating immune cells and risk sore. Figure S1. The representative results of the evaluation of tumor infiltrating immune cells with risk assessment model.","description":"","filename":"SupplementalInformation.docx","url":"https://assets-eu.researchsquare.com/files/rs-214493/v1/acec4aedc8f54e9a0ea29ce8.docx"}],"financialInterests":"","formattedTitle":"\u003cp\u003eA Novel Immune-Related lncRNA Signature predict the Immune Landscape and Prognosis for Head and Neck Squamous Cell Carcinoma\u003c/p\u003e","fulltext":[{"header":"Background","content":" \u003cp\u003eHead and neck tumors are the sixth most common malignancy in the world, head and neck squamous cell carcinoma (HNSCC) is the most common histological type [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]. HNSCC encompasses a heterogeneous group of malignancies that arise in the oral cavity, larynx and pharynx and are mainly associated with tobacco, alcohol consumption and human papillomavirus (HPV) infection [\u003cspan additionalcitationids=\"CR3\" citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e]. However, due to the early symptoms of HNSCC are not obvious, the 5-year survival rate is still below 50%. Furthermore, HNSCC is characterized by high lymph node metastasis rate, local recurrence and metastasis can reduce the survival rate to 35% [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. Currently, HNSCC is estimated to have an annual incidence of 900,000 and causes 450,000 death every year [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eHNSCC is a disease characterized by profound immunosuppression [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. immune checkpoint inhibitors (ICIs) therapy is becoming a new hope for HNSCC treatment. In recent years, considerable progress has been made in ICIs therapy for HNSCC. However, the response rate of recurrent or metastatic HNSCC for PD-1/PD-L1 inhibitors was only found to be 13.3\u0026ndash;22% in previous clinical trials [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. Therefore, biomarkers which could effectively predict the efficacy of ICIs are crucial for patient selection. Generally, PD-L1 expression level is used as the representative marker for predicting the efficacy of ICIs. However, for most tumors, only the detection of PD-L1 alone cannot screen the dominant population [\u003cspan additionalcitationids=\"CR8\" citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eLong non-coding RNAs (lncRNAs) refer to functional RNA molecules with a length more than 200 nucleotides and have minimal or no function to encode proteins [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]. It has become increasingly recognized as lncRNAs can interact with DNA, RNA and proteins which regulate gene expression at the transcription and post-transcription levels and is associated with various types of cancer [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e, \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. Recently, several lncRNA signatures have been established for prognosis prediction in some cancers, including HNSCC. Mao et al. [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e] identified an eight lncRNA signature is an independent prognostic factor of HNSCC patients. Bao et al. [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e] identify genome instability-associated lncRNAs and constructed a signature to evaluate the prognosis of patients with breast cancer.\u003c/p\u003e \u003cp\u003eIt is know that the diagnosis, evaluation, and treatment of cancer is closely related to the tumor immune microenvironment, especially to infiltrating immune cells [\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e, \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e, \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. Now research confirms that, lncRNAs are critical regulators of gene expression in the immune system. Research shows that lncRNAs direct the expression of genes related to the activation of immune cells that directing the development of diverse immune cells so resulting in tumor immune cell infiltration and altering the immune microenvironment [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e, \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. Therefore, in this study, we tried to identify immune-related lncRNAs (irlncRNAs) to investigate the relationship between irlncRNAs and immune infiltration, construct a powerful prognostic model for risk assessment of HNSCC and predict the response of patients to immunotherapy. Different from simple genes in the past, we utilized a novel modeling algorithm, paring, and iteration to construct an irlncRNA signature, which adopted two-biomarker combinations will help improve the accuracy of the cancer prognostic model.\u003c/p\u003e "},{"header":"Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e\n\u003ch2\u003eData collectin\u003c/h2\u003e\n\u003cp\u003eHNSCC RNAseq expression profile data were downloaded from The Cancer Genome Atlas (TCGA) database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://portal.gdc.cancer.gov/\u003c/span\u003e\u003c/span\u003e). The database contained 501 cases tumor tissue sequencing information with HNSCC patients and 44 cases of normal tissue. Clinical and follow-up data of patients with HNSCC were retrieved from the HNSC project of TCGA, valid data were extracted after eliminating the data with follow-up time of 30 days and repeated data. \u0026ldquo;Ensembl_Stable_id\u0026rdquo; is converted to \"Gene Symbol\" in the RNAseq expression profile according to the GTF file which downloaded from Ensembl (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://asia.ensembl.org\u003c/span\u003e\u003c/span\u003e), and annotation to distinguish the mRNAs and lncRNAs for subsequent analysis. The list of immune-related genes (ir-genes) was downloaded from the ImmPort database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.immport.org\u003c/span\u003e\u003c/span\u003e), and correlation analysis was conducted between ir-genes and all lncRNAs to screen for irlncRNAs. Immune gene correlation coefficients greater than 0.4 or less than \u0026minus;\u0026thinsp;0.4 and p value less than 0.001 were considered as irlncRNAs. The differentially expressed irlncRNAs (DEirlncRNAs) were identified through the \u0026ldquo;limma\u0026rdquo; R package, FDR\u0026thinsp;\u0026lt;\u0026thinsp;0.05 and |log\u003csub\u003e2\u003c/sub\u003eFC|\u0026gt;1.5 were defined as DEirlncRNAs.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec4\" class=\"Section2\"\u003e\n\u003ch2\u003ePairing DEirlncRNAs\u003c/h2\u003e\n\u003cp\u003eThe DEirlncRNAs were cyclically singly paired, and construct a 0-or-1 matrix. If the expression level of the former lncRNA was higher than that of the latter, it was defined as 1; otherwise, it was defined as 0. Then the constructed 0-or-1 matrix is further screened. When the amounts of DEirlncRNA pairs of which expression quantity was 0 or 1 accounted for more than 20% and less than 80% of total pairs, it was considered a valid match.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec5\" class=\"Section2\"\u003e\n\u003ch2\u003eBuilding the prognostic pairing DEirlncRNAs-based signature\u003c/h2\u003e\n\u003cp\u003eUnivariate Cox analysis was performed by \u0026ldquo;survival\u0026rdquo; R package to screen DEirlncRNA pairs with prognostic values and P\u0026thinsp;\u0026lt;\u0026thinsp;0.01 DEirlncRNA pairs is considered to be significantly related to the prognosis of HNSCC patients. Next, Lasso regression analysis performed by \u0026ldquo;glmnet\u0026rdquo; and \u0026ldquo;survival\u0026rdquo; R package to select the minimum error value to eliminate over-fitting data, the optimal penalty parameter \"lambda\" value was calculated via a 1,000 cross-validations. Then multivariate analyses were performed using the Cox regression model to identify the prognostic-associated DEirlncRNA pairs factors in HNSCC. Based on the coefficients from the multivariate regression analysis and the status of DEirlncRNA pairs, we constructed a DEirlncRNA pairs signature(DEirLnc-PSig) for outcome prediction as follows:\u003c/p\u003e\n\u003cp\u003e\u003cimg src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAUgAAAA5CAYAAABUHbLjAAAKWUlEQVR4Ae2bjZHUOhCELwViIAVyIARiIAUyIAMyIAIiIAESIIPLYanvqvrRzJO8ttf/26ryyZZG89MjtWWv7+WWEgSCQBAIAk0EXpqtaQwCQSAIBIFbCDKTIAgEgSDQQSAE2QEmzUEgCASBEGTmQBAIAkGgg0AIsgNMmoNAEAgCIcjMgSAQBIJAB4EQZAeYNAeBIBAEQpCZA0EgCASBDgIhyA4waQ4CSyLw+/fvGwfl169f/50vaSO6lkcgBLk8ptEYBP5B4Pv377f379/fPn78ePv8+fPtw4cPt3fv3v0jk4tjIhCCPGZe4tVOCHz79u328vKyyPH169e3KNg5/vz5840U2T2+vr6+keROIcbsBARCkBPAiuj1ERB5QZLs8rieUny8CJLxnH/58uVNFTtKnU/RHdntEQhBbo95LB4cAXZ52kXyWDy1sGNkvBMken78+PGmisdsJ8yp+iO/HQIhyO2wjqUTIeCP2pxPLbxzrAQpHewgP336NHl3qvGpt0MgBLkd1rF0MgQgMe0k2VVOKewYnSCnjI3scRAIQR4nF/HkYAjwPpGdICRJPfV95MHCiTszEAhBzgAtQ54HAX8fybvDlOdCIAT5XPlOtDMQ4FFZj9q8PzxakW+qj+bfmf0JQZ45e/F9MwR4pwgB8emP/iNmM+MDhvCpllZblcn1OAT+j+64cZE6MALscuZ8nnLEkPg0hh9L+NB6z8L7R8gR8uE/YY5SWmTYajuKv2fzYzRB6mU14OtgEdZHDq7VX2tftJz7dQ843gFpYuLD0Usrfvzn/ZW/5Gfhs9CEESTgv5TSPudXUOxs9a5sbA4fzRk7NrCqc+1RvVPHQ9LK15zcTLU3Rh5/vNRr7/Nz5iKy+jbT++acCxevW58y8YH8kW4w92L9F90BaU0O3cmZtHo3Uxdka+EwuZ0QWzI987LT6z9ae41NpMmEobSwhER90THR/HpMjMg7xmPGjJHRh89VtsZZ++dcM5dacWsH5zeROfofHcMCFwloLTyq85Hx8kX1WF2ag0vhCdHig14/oF8bA/eJOaN14O1HPZ9NkAqIyQwwfndfeuHIhmwevW7FTxs4UYin3kX5GLlFDGNjXZNANPnH+vKIHLj0cICcwHHvgo/kEgLwp4I9/NKckm2ua5v65tTkfsy/RZKz+oTXWgdzfNhzzMMEqa267yIrMNwx/G6rHZVPdj3CI6e7tPoBfyjp2KYfOQ4mLoeTNiBDQrKD/JjEz0lOjR8dtGkCKf7qn2x5PGrjTo8O/PZDGKFL+jWGWouZfuEKNvW/Q1gIknXsGMe12wRj2mmTfWwxF/Bd8v6Ihf9q91gYL5KRfbflsTCOPu1SvE/nmiuuo3Uu+Tk19hXL3rshYmuVXjuymkdD85/5wcF6BFPmh6/xahOdro+cgpHa0KM8KN9VxxGv2+g2PFWA1LUQuC8UJUCAUNdxyPgY6adNyVC/Jn2169fYgCBkh4lLglTQybX60Sn9klmqRq/rFplgU0VkhRy+1UI8Lk9smqAsUK41+RhLn9t0feiiT7Ez8WmDcChco0/EU7FTblwn5+h0m/gAybEAOOiTz8grj/hNvwgP+yqM8bjVrhq/W3ipf6saH/CFY8/Ss99rl6/038OReau5QB57xEa762MeaZzmFHbJK/PjTGV0drVItMg8SMDxhcK5XwNWHVdlpL+VNC0st1nP8cEXliaw5EiME4raWzW+oW/o8Piqjjoe2+6b5IkZbLBD7aXGw7VjiD73gfOWDXRWXbRxs+jJ084YFeVG16qxKR/qIkFGu0/JV720V2yG4kC+FYv0b1mz8Icw3MoXz5Pb7LUjg+/09whPesgX+VFOevJaa+jUwXzWDVj6INmxa1Bj9q7/roI7nmiR+CJliBaGLzZfOD21Vaann/GthVX1khj3QfokV/vVvkZdY7tng4nEYvOdVPWXHR6TDrw5fEeJfk3ilq2qqyWPbfRzIM+hUrFUu8cpGY31WvKtPFa/67XGqkav51ntqmXD7bfOJT+3FnHMHb/UOMXm+mgbKty48H+oMB84yCuYUvvTgI+F9JiP9woyrQ3QvXF79g8jaZ5pAVB7AUQS4ncLXziS5a5Fu0qV6elHXpNeY6khCU8yPvjCkT6NITner/Y16hpbtUF/nSi0uX+teCBR2jnqIw/je3fnqgt/nJDRBTGqVLwrlpLzOLUrqfNDstRVL2017nrt4zlvxVJl1r4mDvAj5r0LeFCoddzzacpOjnnam1eyU2/Wavda84N1e6YymyAJWBO+AugLR2CQFNpVqowWYWuByY7GUmPT9dWFo22/EoIOEqlJjZ3qt+t/5By/nLyrLvW7Lyw47uwqNR70tbCRfMVX7dTowqZj4QucPhEkPnHNGBXPDTdC4YYchwrn+Km46uJq5RF56UOP6yQm6ZIN/Ko3F/VtUQsLz9UWdns2PE89mdquzQIxPIql8CC3QwVb+MocJK9nKX9XwR2PAZUA/WAyV4AFhMvpXIvJZdiBUqSfheu7UT1+SofX0gfgakfex6CXosRIjjEijDeBhf54bL2JAGb+KKsJKxdqPLQTq3xX7WSEXcUqPaqRRxZsde55Y6z6ICv6kBNpun38ADePUzkEd4/LMfacSK/HiT4KtvGFQ22KQ4uxkqb6166JG4zl/9r2xugnT1OLcnSP1Mboxb4O8tMr9JFT5iFz4SxlOrpniexCfjKhmFhO6JzTpknONROwNUmZwJI7MyyQt26Ke8QBsUCQnoc9/HCbcwjSx+d8GIEQ5DA+h+iF3FiYfueFCGnzXRZyLQK5AkGya+QG4BhsmRztmFs3oC39iK1tEQhBbov3LGvsWHgchSD0OMPuUY+2rhQ5DhU9TjHOyVT9Z6i1W97Lf5HzFXbhZ8j3kXwMQR4pGwv5ApG0dpILqd9UDTs2SH7PnRs3I46U50MgBPl8OU/EExBg18jOfa8fhia4GtEVEAhBrgBqVF4DAXatvJpovcq4RoSJ4h4CIch7CKX/KRHgvecjn/RAqvU1B++Gp/7IxOuSvd69PmXiS9AhyAJILoMACPDek0fruZ/0ML4SpL9H1SdLyPghAqWGUP1TrmRmewRCkNtjHosHR4AdG4/WTmhTXNajuX9NAAn6fwxBgMjVoxIy70Dz6/kU9JeVDUEui2e0nRwBfdKjz6keqZ3YIEu/HgtTCHIsUuvIhSDXwTVaT4qAfzf6CDky1gkRgvQd6b1HbMEXghQS+9QhyH1wj9UnQ4AffCj6XGjsIzZE6kT7ZLDtHm4IcvcUxIFnQIAdJT+46EeYezEjx3tLxvFjEecp2yMQgtwe81h8QgQgvPoDzBAMyNYfcIbk07cOAiHIdXCN1iAQBC6AQAjyAklMCEEgCKyDQAhyHVyjNQgEgQsgEIK8QBITQhAIAusgEIJcB9doDQJB4AIIhCAvkMSEEASCwDoIhCDXwTVag0AQuAACIcgLJDEhBIEgsA4CIch1cI3WIBAELoBACPICSUwIQSAIrINACHIdXKM1CASBCyAQgrxAEhNCEAgC6yDwBxeAHT78N3yoAAAAAElFTkSuQmCC\" alt=\"\" /\u003e\u003c/p\u003e\n\u003cp\u003ewhere \u0026beta;\u003cem\u003ei\u003c/em\u003e indicates the coefficient for each DEirlncRNA pairs and \u003cem\u003ePi\u003c/em\u003e indicates the DEirlncRNA pairs status of 0-or-1 matrix.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec6\" class=\"Section2\"\u003e\n\u003ch2\u003eValidation of the constructed risk model\u003c/h2\u003e\n\u003cp\u003eA time-dependent receiver operating characteristic (ROC) curve and area under curve(AUC) was estimated using the \u0026ldquo; survivalROC\u0026rdquo; R package to predict the accuracy of DEirLnc-PSig for HNSCC patients. The 1-, 2-, and 3-year ROC curves and AUC of the model were plotted. The AIC values of every points of the 3-year ROC curve were evaluated to identify the maximum inflection point, which was considered as the cut-off point to classify patients into the high-risk group with high DEirLnc-PSig or low-risk group with low DEirLnc-PSig.\u003c/p\u003e\n\u003cp\u003eThe Kaplan\u0026ndash;Meier method was used to construct survival curves, and the log-rank test was used to assess the difference in survival between the high- and low-risk groups with a significant level of 0.05. The R packages utilized in these steps included \u0026ldquo;survival\u0026rdquo; and \u0026ldquo;survminer\u0026rdquo;. To verify the clinical application value of the DEirLnc-PSig, multivariate Cox regression and stratified analysis were used to assess the independence of DEirLnc-PSig from other key clinical factors. A forest map was used to demonstrate the results. The R packages utilized in these operations were \u0026ldquo;survival\u0026rdquo; and \u0026ldquo;survminer\u0026rdquo;. In addition, we compared our defined DEirLnc-PSig with two previously published lncRNA signatures to further compare the prognostic model performance.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec7\" class=\"Section2\"\u003e\n\u003ch2\u003eInvestigation of tumor-infiltrating immune cells\u003c/h2\u003e\n\u003cp\u003eTo explore the differences in immune cell infiltration between different risk groups, we used the TIMER, CIBERSORT, XCELL, QUANTISEQ, MCPcounter, EPIC, and CIBERSORT methods to assess the samples from the TCGA project of the HNSCC dataset. Differences in the number of different types of immune infiltrating cell between high- and low-risk groups of the constructed model were analyzed using the Wilcoxon's signed-rank test and the results are shown in a box chart. The correlation between the risk scores and immune cell infiltration were evaluated using Spearman's correlation analysis. The correlation coefficients of the results were shown in a lollipop diagram. The significance threshold was set as p\u0026thinsp;\u0026lt;\u0026thinsp;0.05. The R packages utilized in these operations were \u0026ldquo;limma\u0026rdquo;, \u0026ldquo;scales\u0026rdquo;, \u0026ldquo;ggplot2\u0026rdquo;, \u0026ldquo;ggtext\u0026rdquo; and \u0026ldquo;ggpubr\u0026rdquo;.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec8\" class=\"Section2\"\u003e\n\u003ch2\u003eExploration of the significance of the model in the clinical treatment\u003c/h2\u003e\n\u003cp\u003eWe calculated the half inhibitory centration (IC\u003csub\u003e50\u003c/sub\u003e) of common administrating chemotherapeutic drugs in the TCGA project of the HNSCC dataset by \" pRRophetic\" R package to evaluate the model in the clinic for head and neck carcinoma treatment. Antitumor drugs such as cisplatin, docetaxeland, methotrexate and paclitaxel are recommended for HNSCC treatment by NCCN guidelines. The difference in the IC\u003csub\u003e50\u003c/sub\u003e between the high- and low-risk groups was compared by Wilcoxon\u0026rsquo;s signed-rank test and the results are shown as box drawings.\u003c/p\u003e\n\u003c/div\u003e\n\u003cdiv id=\"Sec9\" class=\"Section2\"\u003e\n\u003ch2\u003eAnalyses of the immunosuppressive molecules expressing related to ICIs\u003c/h2\u003e\n\u003cp\u003eThe relationship between the model and the expression level of genes related to ICIs were performed by \"limma\" and \u0026ldquo;ggpubr\u0026rdquo; R package and violin plot visualization.\u003c/p\u003e\n\u003c/div\u003e"},{"header":"Results","content":" \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eIdentification of irlncRNAs and differentially expressed irlncRNAs\u003c/h2\u003e \u003cp\u003eA study process flow chart is given in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e. A total of 805 irlncRNAs were identified by co-expression analysis of ir-genes and lncRNAs. The differential expression analysis was subsequentially carried out using \u0026ldquo;limma\u0026rdquo; R package, 132 DEirLncRNAs (75 upregulated and 16 downregulated) were identified between the tumor and normal tissues (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003eEstablishment of DEirlncRNA pairs and a DEirlncRNA pairs signature\u003c/h2\u003e \u003cp\u003eUsing an iteration loop and a 0-or-1 matrix screening among 132 DEirlncRNAs and 6,686 valid DEirlncRNA pairs were identified. By conducting the univariate Cox regression analysis, 500 survival-related DEirlncRNA pairs were selected (p\u0026thinsp;\u0026lt;\u0026thinsp;0.01). To further avoid result overfitting, the Lasso regression analysis was employed to screen out 33 DEirlncRNA pairs showing the highest correlation with survival (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ea, b). Those 33 DEirlncRNA pairs were further subjected to multivariable Cox regression analysis, finalizing the most significant 18 DEirlncRNA pairs was developed as a DEirlncRNA pairs signature for HNSCC patients by the stepwise elimination method and obtaining corresponding coefficients (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ec). We then applied the 1-, 2-, and 3-year ROC curves to verify the predictive capacity of the DEirLnc-Psig; the AUCs of the ROC were almost all ~\u0026thinsp;0.75, suggesting a satisfactory predictive performance (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea). In addition, we compared the 3-year ROC curve with other clinical features and the superiority of DEirLnc-PSig have been further demonstrated (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eb).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003eClinical evaluation by risk assessment model and independence of the DEirLnc-PSig from other clinical factors\u003c/h2\u003e \u003cp\u003eWe recognized the maximum inflection point as the cut-off point on the 3-year ROC curve using the Akaike information criterion (AIC) values and according to the cut-off point (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ec, d), HNSCC patients from the TCGA cohort was divided into a low-risk group(n\u0026thinsp;=\u0026thinsp;255) and a high-risk group(n\u0026thinsp;=\u0026thinsp;235). We then compared the survival of low-risk group with the high-risk group and found that the DEirLnc-PSig possessed powerful capacity to predict the prognosis of HNSCC patients between the low- and high-risk group, in which the low-risk group are significantly better than the high-risk group (p\u0026thinsp;\u0026lt;\u0026thinsp;0.001, log-rank test) (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ea, b, c).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eTo evaluate whether the prognostic value of the DEirLnc-PSig was independent of common clinical variables, univariate and multivariate Cox regression analyses were performed on age, gender, pathologic grade, clinical stage and our DEirLnc-PSig-based prognostic risk score model (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ea, b). The results demonstrated that risk score showed statistical differences by univariate and multivariate Cox regression analysis which presented as an independent prognostic predictor. In addition to the DEirLnc-PSig, there were other two clinical factors, age and clinical stage, that were observed to be significant in the multivariate analysis. Therefore, stratification analysis was performed to determine whether the DEirLnc-PSig possessed a prognostic value that was independent of the age and clinical stage. Patients in the TCGA set were stratified into a young-patient group (\u0026thinsp;\u0026lt;\u0026thinsp;=\u0026thinsp;65years, n\u0026thinsp;=\u0026thinsp;124) and an old-patient group (\u0026gt;\u0026thinsp;65 years, n\u0026thinsp;=\u0026thinsp;253). Using the DEirLnc-PSig, patients in each age group could be further divided into high-risk or low-risk group. There was a significant difference in OS between the high- and low-risk groups in the young-patient group (p\u0026thinsp;\u0026lt;\u0026thinsp;0.001, log-rank test) as was in the old-patient group (p\u0026thinsp;\u0026lt;\u0026thinsp;0.001, log-rank test) (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ec, d). Next, all HNSCC patients were then stratified into the early-stage group (stage I and II, n\u0026thinsp;=\u0026thinsp;70) and the late-stage group (stage III and IV, n\u0026thinsp;=\u0026thinsp;307). Using the DEirLnc-PSig, patients in each stage group could be further divided into high-risk or low-risk group. There was a significant difference in OS between the high- and low-risk groups in the early-stage group (p\u0026thinsp;=\u0026thinsp;0.013, log-rank test) as was in the late-stage group (p\u0026thinsp;\u0026lt;\u0026thinsp;0.001, log-rank test) (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ee, f). These results indicated that the DEirLnc-PSig was an independent prognostic factor associated with OS in HNSCC patients.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003ePerformance comparison of the DEirlncRNA pairs signature with existing lncRNA-related signatures in survival prediction\u003c/h2\u003e \u003cp\u003eWe further compared the prediction performance of the DEirLnc-PSig with two recently published lncRNA signatures: 3-lncRNA signature derived from Jiang\u0026rsquo;s study (hereinafter referred to as JiangLncSig) [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e] and 5-lncRNA signature derived from Liu\u0026rsquo;s study (hereinafter referred to as LiuLncSig) [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e] using the same TCGA patient cohort. As shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e, the AUC at 3 years of OS for the DEirLnc-PSig is 0.809, which is significantly higher than that of JiangLncSig (AUC\u0026thinsp;=\u0026thinsp;0.632) and LiuLncSig (AUC\u0026thinsp;=\u0026thinsp;0.5). These results demonstrated the better prognostic performance of the DEirLnc-PSig in predicting survival than two recently published lncRNA signatures.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003eEstimation of tumor-infiltrating immune cells and immunosuppressive molecules with risk assessment model\u003c/h2\u003e \u003cp\u003eSince the model was established by irlncRNAs, we consequently investigated whether the model was related to the tumor immune microenvironment. Correlation analysis of immunocyte infiltration revealed that the high-risk group was more positively associated with CD4\u0026thinsp;+\u0026thinsp;T cells, cancer associated fibroblast, monocyte and macrophages, whereas were negatively associated with CD8\u0026thinsp;+\u0026thinsp;T cells and T cell regulatory (Tregs), as revealed by the Wilcoxon signed-rank test (see supplementary file, Table S1 and Figure S1). In addition, the immune scores in the low-risk group were significantly higher than those in the high-risk group (p\u0026thinsp;=\u0026thinsp;0.0038) in the analysis with xCell, which suggesting that the tumor cells in the low-risk group had more immune cell infiltration than those in the high-risk group. A detailed Spearman correlation analysis was conducted, and the resulting diagram exhibited a lollipop shape, as shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003ea. Immunotherapy with ICIs has changed HNSCC treatment strategies dramatically, improve survival of HNSCC patients. We investigated whether the risk model was related to ICI-related biomarkers and discovered that high risk scores were positively correlated with low expression of PDCD1 (p\u0026thinsp;\u0026lt;\u0026thinsp;0.001, Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003eb), LAG3 (p\u0026thinsp;\u0026lt;\u0026thinsp;0.01, Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003ec) and CTLA4 (p\u0026thinsp;\u0026lt;\u0026thinsp;0.001, Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003ed).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003eAnalysis of the correlation between the risk model and chemotherapeutics\u003c/h2\u003e \u003cp\u003eAs common therapy for HNSCC, we attempted to identify associations between risk and the efficacy of chemotherapeutics in treating head and neck squamous cell carcinoma in the TCGA project of the HNSCC dataset. We showed that a high risk score was associated with a lower IC\u003csub\u003e50\u003c/sub\u003e of chemotherapeutics such as docetaxel (p\u0026thinsp;\u0026lt;\u0026thinsp;0.001), paclitaxel (p\u0026thinsp;=\u0026thinsp;0.02), whereas it was associated with a higher IC\u003csub\u003e50\u003c/sub\u003e for methotrexate (p\u0026thinsp;\u0026lt;\u0026thinsp;0.001), which indicated that the model acted as a potential predictor for chemosensitivity (Fig.\u0026nbsp;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e9\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e "},{"header":"Discussion","content":" \u003cp\u003eDespite the effectiveness of targeted therapies and immunotherapy agents in HNSCC, not all patients respond to treatment and development of resistance is almost universal after a time. In clinical practice, ICIs therapy does not appear to improve survival in HNSCC patients. Precision medicine approaches that seek to individualise therapy through the use of predictive biomarkers and stratification strategies, have been the key to improve therapeutic success in HNSCC [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e, \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eIn recent years, many studies have shown that lncRNA is an important and significant factor in the current and future prognosis of patients with HNSCC [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]. For this reason, many prognostic lncRNA expression signatures have been applied to diagnose and predict prognosis of HNSCC. However, it should be pointed out that HNSCC are a type of immunosuppressive disease [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e]. In the case of HNSCC, a deeper understanding of the heterogeneity and complex of the tumor immune microenvironment is essential to develop sensitive predictive biomarkers for providing information of tumor microenvironment and help to improve the efficacy of tumor immunotherapy. An ideal signature would not only to predict the prognosis of cancer patients, but also reflect the characteristics of the tumour and immune microenvironment that have an impact on disease progression and treatment response. Therefore, we established a prediction model of immune-related lncRNAs and evaluated the relationships between the risk score draw from DEirLnc-PSig, immune checkpoints, tumor-infiltrating immune cell abundances and the response to chemotherapy.\u003c/p\u003e \u003cp\u003eUnlike most signatures, which are based on quantified transcript expression levels, we attempted to construct a reasonable model with two-lncRNA combinations and did not adopt their exact expression levels in the signature. Thus, pairs with higher or lower expression only had to be detected instead of examining specific expression values of every lncRNA. The results of this study showed that this novel model had a good predictive performance and achieved a better prediction performance with an AUC value of 0.809, which is higher than other two expression based lncRNA signature, risk score derived from the DEirLnc-PSig is reliability and independent for predicting HNSCC prognosis.\u003c/p\u003e \u003cp\u003eIn this study, some of the irlncRNAs in the process of modeling that have been already identified play an important role in malignant phenotypes of various cancer types, such as LINC01063 [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e, \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e], HOXC-AS1 [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e, \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e], LINC02454 [\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e, \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e] and SCAT1 [\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e]. Furthermore, previous reports have shown that LINC00460 [\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e], C5orf66-AS1 [\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e] and TMPO-AS1 [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e] are associated with prognosis and treatment outcomes in HNSCC. Jiang et al. [\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e] revealed that LINC00460 promotes EMT in HNSCC by facilitating peroxiredoxin-1 into the nucleus. Lu et al. [\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e] first report that C5orf66-AS1 was able to prevent OSCC progression by inhibiting OSCC cell growth and metastasis via regulating CYC1 expression. Xing et al. [\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e] disclosed that a novel TMPO-AS1/miR-320a/SOX4 pathway associated with nasopharyngeal carcinoma(NPC) progression, suggesting that TMPO-AS1 may be a potential therapeutic target for NPC. However, there are also some lncRNAs that have been disclosed for the first time in this model which suggested that the proposed model can identify novel biomarkers for further research.\u003c/p\u003e \u003cp\u003eTumor-infiltrating immune cells which are integral to the tumor immune microenvironment (TIME) are recognized as highly important for prognosis prediction, while the responsiveness to immune checkpoint blockade may be subject to the features of TIME. To explore the relationship between risk scores and tumor-infiltrating immune cells, we used seven common acceptable methods to estimate the immune infiltrating cell, including xCell [\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e, \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e], TIMER [\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e, \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e], QUANTISEQ [\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e, \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e], MCPcounter [\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e], EPIC [\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e], CIBERSORT-abs [\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e], and CIBERSORT [\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e]. Multiple studies have indicated that greater infiltration of CD8 T cells and high immune scores correlated with better response to standard chemotherapy and to immune checkpoint blockade therapy in tumor patients [\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e, \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e]. By integrating analyzing, our results revealed that the high-risk group identified by DEirLnc-PSig was negatively associated with CD8\u0026thinsp;+\u0026thinsp;T cell and immune scores. Correspondingly, the high-risk group was also negatively correlated with checkpoint-related biomarkers such as PDCD1, LAG3 and CTLA4. Our results demonstrated that patients in high-risk groups identified by our model may not benefit from immunotherapy. But these conclusions need to be verified through further studies especially supported by clinical trials. On the other hand, we checked several more drugs like cisplatin, paclitaxel, doxorubicin and methotrexate, we observed from the estimated IC\u003csub\u003e50\u003c/sub\u003e of these chemotherapeutic drugs, that the high-risk group could be more sensitive to chemotherapies than those in low-risk group. It would be more helpful to identify patients at high clinical risk and give sensitive chemotherapy rather than immunotherapy. However, more benchwork and clinical studies are needed to further validation.\u003c/p\u003e "},{"header":"Conclusions","content":" \u003cp\u003eIn summary, our study revealed a comprehensive landscape of the tumor immune microenvironment in HNSCC, and provides a promising prognostic signature for HNSCC. Patients with the low risk scores could benefit more from ICIs therapy while patients with the high risk scores may benefit more from chemotherapy. However, further investigations are needed to verify the accuracy in estimating prognoses and to test its clinical utility in patient management.\u003c/p\u003e "},{"header":"Abbreviations","content":"\u003cp\u003eHNSCC: Head and neck squamous cell carcinomas; TCGA: The Cancer Genome Atlas; DEirlncRNA: Differently expressed irlncRNA; OS: Overall survival; AUC: Areas under curve; AIC: Akaike information criterion; ROC: Receiver operating characteristic; HR: hazard ratio; HPV: human papillomavirus; ICIs: immune checkpoint inhibitors; lncRNAs: Long non-coding RNAs; irlncRNAs: Immune-related lncRNAs; ir-genes: immune-related genes; DEirlncRNAs: Differentially expressed irlncRNAs; IC\u003csub\u003e50: \u003c/sub\u003eInhibitory centration; Tregs: T cell regulatory; TIME: Tumor immune microenvironment;\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe datasets analyzed during the current study are available in the TCGA repository, https://portal.gdc.cancer.gov/.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors report no competing interest in this study.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis work was supported by the National Key R\u0026amp;D Program of China (No 2020YFC1316903), the 5010 Clinical Research Program of Sun Yat-sen University (No. 2017004) and the National Natural Science Foundation of China (No. 81972528).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026rsquo; contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eZMC designed the research, collected the data, analyzed the date and wrote the manuscript. LC took part in designing the research. YL collected the data, analyzed the date. WBL took part in designing the research and solved the disagreements between ZMC and LC. All authors read and approved the final manuscript.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNone.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor information\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eZMC, LC and YL contributed equally to this work.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eBray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018 Nov;68(6):394\u0026ndash;424.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKatiyar SK. Emerging Phytochemicals for the Prevention and Treatment of Head and Neck Cancer. Molecules. 2016 Nov 24;21(12):1610.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFleming JC, Woo J, Moutasim K, Mellone M, Frampton SJ, Mead A, et al. HPV, tumour metabolism and novel target identification in head and neck squamous cell carcinoma. Br J Cancer. 2019 Feb;120(3):356\u0026ndash;67.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eNelson HH, Pawlita M, Michaud DS, McClean M, Langevin SM, Eliot MN, et al. Immune response to HPV16 E6 and E7 proteins and patient outcomes in head and neck cancer. JAMA Oncol. 2016;3:178\u0026ndash;85.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePosner M, Vermorken JB. Induction therapy in the modern era of combined-modality therapy for locally advanced head and neck cancer. Semin Oncol. 2008 Jun;35(3):221\u0026ndash;8.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRomano E, Romero P. The therapeutic promise of disrupting the PD-1/PD-L1 immune checkpoint in cancer: Unleashing the CD8 T cell mediated anti-tumor activity results in significant, unprecedented clinical efficacy in various solid tumors. J Immunother Cancer. 2015;3:15.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang SM, Jiang B, Deng Y, Huang SL, Fang MZ, Wang Y. Clinical significance of MLH1/MSH2 for stage II/III sporadic colorectal cancer. World J Gastrointest Oncol. 2019 Nov;15(11):1065\u0026ndash;80. 11(.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchepisi G, Brighi N, Cursano MC, Gurioli G, Ravaglia G, Altavilla A, et al. Inflammatory Biomarkers as Predictors of Response to Immunotherapy in Urological Tumors. J Oncol. 2019 Sep;19:2019:7317964.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu Y, Zhou S, Du Y, Sun L, Jiang H, Zhang B, et al. Efficacy and safety of programmed death 1 inhibitors in patients with advanced non-small cell lung cancer: a meta-analysis. Cancer Manag Res. 2019 May;21:11:4619\u0026ndash;30.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKopp F, Mendell JT. Functional Classification and Experimental Dissection of Long Noncoding RNAs. Cell. 2018 Jan 25;172(3):393\u0026ndash;407.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSchmitt AM, Chang HY. Long Noncoding RNAs: At the Intersection of Cancer and Chromatin Biology. Cold Spring Harb Perspect Med. 2017 Jul;5(7):a026492. 7(.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMao R, Chen Y, Xiong L, Liu Y, Zhang T. Identification of a nomogram based on an 8-lncRNA signature as a novel diagnostic biomarker for head and neck squamous cell carcinoma. Aging (Albany NY). 2020 Oct 22;12(20):20778\u0026ndash;20800.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBao S, Zhao H, Yuan J, Fan D, Zhang Z, Su J, Zhou M. Computational identification of mutator-derived lncRNA signatures of genome instability for improving the clinical outcome of cancers: a case study in breast cancer. Brief Bioinform. 2020 Sep 25;21(5):1742\u0026ndash;1755.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMartinez LM, Robila V, Clark NM, Du W, Idowu MO, Rutkowski MR, et al. Regulatory T cells control the switch from \u003cem\u003ein situ\u003c/em\u003e to invasive breast cancer. Front Immunol. 2019;10:1942.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePeng D, Wang L, Li H, Cai C, Tan Y, Xu B, et al. An immune infiltration signature to predict the overall survival of patients with colon cancer. IUBMB Life. 2019 Nov;71(11):1760\u0026ndash;70.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi AL, Zhu YM, Gao LQ, Wei SY, Wang MT, Ma Q, Zheng YY, et al. Exploration of the Immune-Related Signatures and Immune Infiltration Analysis in Melanoma. Anal Cell Pathol (Amst). 2021 Jan;2021:16. 4743971.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAtianand MK, Caffrey DR, Fitzgerald KA. Immunobiology of Long Noncoding RNAs. Annu Rev Immunol. 2017 Apr 26;35:177\u0026ndash;198.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen YG, Satpathy AT, Chang HY. Gene regulation in the immune system by long noncoding RNAs. Nat Immunol. 2017 Aug;22(9):962\u0026ndash;72. 18(.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJiang H, Ma B, Xu W, Luo Y, Wang X, Wen S, et al. A Novel Three-lncRNA Signature Predicts the Overall Survival of HNSCC Patients. Ann Surg Oncol. 2020 Oct 23.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLiu G, Zheng J, Zhuang L, Lv Y, Zhu G, Pi L, et al. A Prognostic 5-lncRNA Expression Signature for Head and Neck Squamous Cell Carcinoma. Sci Rep. 2018 Oct;15(1):15250. 8(.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLee TW, Lai A, Harms JK, Singleton DC, Dickson BD, Macann AMJ, et al. Patient-Derived Xenograft and Organoid Models for Precision Medicine Targeting of the Tumour Microenvironment in Head and Neck Cancer. Cancers (Basel). 2020 Dec 12;12(12):3743.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang Y, Wang S, Ren Y, Zhou X. The Role of lncRNA Crosstalk in Leading Cancer Metastasis of Head and Neck Squamous Cell Carcinoma. Front Oncol. 2020 Oct 2;10:561833.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLee TW, Lai A, Harms JK, Singleton DC, Dickson BD, Macann AMJ, et al. Patient-Derived Xenograft and Organoid Models for Precision Medicine Targeting of the Tumour Microenvironment in Head and Neck Cancer. Cancers (Basel). 2020 Dec 12;12(12):3743.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDeng X, Bi Q, Chen S, Chen X, Li S, Zhong Z, et al. Identification of a Five-Autophagy-Related-lncRNA Signature as a Novel Prognostic Biomarker for Hepatocellular Carcinoma. Front Mol Biosci. 2021 Jan;11:7:611626.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCheng L, Han T, Zhang Z, Yi P, Zhang C, Zhang S, et al. Identification and Validation of Six Autophagy-related Long Non-coding RNAs as Prognostic Signature in Colorectal Cancer. Int J Med Sci. 2021 Jan 1;18(1):88\u0026ndash;98.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhou C, An N, Cao C, Wang G. lncRNA. HOXC-AS1 promotes gastric cancer via binding eIF4AIII by activating Wnt/β-catenin signaling. J Gene Med. 2020 Sep;22(9):e3202.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDong Y, Li X, Lin Z, Zou W, Liu Y, Qian H, et al. HOXC-AS1-MYC regulatory loop contributes to the growth and metastasis in gastric cancer. J Exp Clin Cancer Res. 2019 Dec;23(1):502. 38(.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCai WY, Chen X, Chen LP, Li Q, Du XJ, Zhou YY. Role of differentially expressed genes and long non-coding RNAs in papillary thyroid carcinoma diagnosis, progression, and prognosis. J Cell Biochem. 2018 Nov;119(10):8249\u0026ndash;59.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTan J, Liu L, Zuo Z, Song B, Cai T, Ding D, et al. Overexpression of novel long intergenic noncoding RNA LINC02454 is associated with a poor prognosis in papillary thyroid cancer. Oncol Rep. 2020 Oct;44(4):1489\u0026ndash;501.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang C, Zhang Z, Zhang G, Xue L, Yang H, Luo Y, et al. A three-lncRNA signature of pretreatment biopsies predicts pathological response and outcome in esophageal squamous cell carcinoma with neoadjuvant chemoradiotherapy. Clin Transl Med. 2020 Aug;10(4):e156.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJiang Y, Cao W, Wu K, Qin X, Wang X, Li Y, et al. LncRNA LINC00460 promotes EMT in head and neck squamous cell carcinoma by facilitating peroxiredoxin-1 into the nucleus. J Exp Clin Cancer Res. 2019 Aug 20;38(1):365.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLu T, Liu H, You G. Long non-coding RNA C5orf66-AS1 prevents oral squamous cell carcinoma through inhibiting cell growth and metastasis. Int J Mol Med. 2018 Dec;42(6):3291\u0026ndash;9.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXing B, Qiao XF, Qiu YH, Li X. TMPO-AS1 Regulates the Aggressiveness-Associated Traits of Nasopharyngeal Carcinoma Cells Through Sponging miR-320a. Cancer Manag Res. 2021 Jan 15;13:415\u0026ndash;425.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAran D. Cell-Type Enrichment Analysis of Bulk Transcriptomes Using xCell. Methods Mol Biol. 2020;2120:263\u0026ndash;76.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017 Nov 15;18(1):220.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 2017 Nov 1;77(21):e108-e110.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi T, Fu J, Zeng Z, Cohen D, Li J, Chen Q, et al. TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020 Jul;48(W1)(2):W509\u0026ndash;14.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFinotello F, Mayer C, Plattner C, Laschober G, Rieder D, Hackl H, et al. Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq data. Genome Med. 2019 May 24;11(1):34. doi: \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1186/s13073-019-0638-6\u003c/span\u003e\u003c/span\u003e. Erratum in: Genome Med. 2019 Jul 29;11(1):50. PMID: 31126321; PMCID: PMC6534875.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePlattner C, Finotello F, Rieder D. Deconvoluting tumor-infiltrating immune cells from RNA-seq data using quanTIsEq. Methods Enzymol. 2020;636:261\u0026ndash;85.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBecht E, Giraldo NA, Lacroix L, Buttard B, Elarouci N, Petitprez F, et al. Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016 Oct 20;17(1):218. doi: \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1186/s13059-016-1070-5\u003c/span\u003e\u003c/span\u003e. Erratum in: Genome Biol. 2016 Dec 1;17 (1):249.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRacle J, Gfeller D. EPIC: A Tool to Estimate the Proportions of Different Cell Types from Bulk Gene Expression Data. Methods Mol Biol. 2020;2120:233\u0026ndash;48.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWu YY, Wang SH, Wu CH, Yen LC, Lai HF, Ho C, et al. In silico immune infiltration profiling combined with functional enrichment analysis reveals a potential role for na\u0026iuml;ve B cells as a trigger for severe immune responses in the lungs of COVID-19 patients. PLoS One. 2020 Dec 2;15(12):e0242900.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen B, Khodadoust MS, Liu CL, Newman AM, Alizadeh AA. Profiling Tumor Infiltrating Immune Cells with CIBERSORT. Methods Mol Biol. 2018;1711:243\u0026ndash;59.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen Y, Chen H, Mao B, Zhou Y, Shi X, Tang L, et al. Transcriptional Characterization Of The Tumor Immune Microenvironment And Its Prognostic Value For Locally Advanced Lung Adenocarcinoma In A Chinese Population. Cancer Manag Res. 2019 Oct 31;11:9165\u0026ndash;9173.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGe PL, Li SF, Wang WW, Li CB, Fu YB, Feng ZK, et al. Prognostic values of immune scores and immune microenvironment-related genes for hepatocellular carcinoma. Aging. 2020 Mar;25(6):5479\u0026ndash;99. 12(.\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Head and neck squamous cell carcinoma, Long non-coding RNAs, Tumor-infiltrating immune cells, Prognosis, TCGA","lastPublishedDoi":"10.21203/rs.3.rs-214493/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-214493/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground:\u003c/strong\u003e The tumor immune microenvironment is known to play an important role in head and neck squamous cell carcinomas (HNSCC). Reliable prognostic signatures that could accurately predict immune landscape and survival in HNSCC patients are vitally needed to promote a better individualized and effective treatment.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eMethods:\u003c/strong\u003e HNSCC transcriptome data and clinical data in The Cancer Genome Atlas (TCGA) were embedded in our study.\u0026nbsp;The differentially expressed irlncRNAs were identified\u0026nbsp;by differential co-expression analysis, and recognized differently expressed irlncRNA (DEirlncRNA) pairs using univariate analysis. Cox and lasso regression analysis was used to identify DEirlncRNA pairs related to overall survival\u0026nbsp;(OS) and build the prediction model.\u0026nbsp;Then, we compared the areas under curve (AUC) with other published lncRNA signatures, counted the akaike information criterion (AIC) values of 3-year receiver operating characteristic curve, and identified the cut-off point to set up an optimal model for distinguishing the high- or low- risk groups among patients with HNSCC.\u0026nbsp;Receiver operating characteristic (ROC)curves and Kaplan–Meier plot curves were used to validate the prediction model. Besides, We then reevaluated them from the viewpoints of clinical factor, tumor-infiltrating immune cells, chemotherapeutics efficacy, and immunosuppressed biomarkers.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eResults:\u003c/strong\u003e We built a risk score model based on 18 DEirlncRNA pairs. The risk model is closely related to the OS of HNSCC patients. The hazard ratio (HR) is 1.376 [95% CI (confidence interval) 1.302-1.453] and log-rank \u003cem\u003eP\u003c/em\u003e-value \u0026lt; 0.0001. Compared with two recently published lncRNA signatures,\u0026nbsp;DEirLncRNA pairs signature has higher AUC score which showed\u0026nbsp;the better prognostic performance. Additionally, the signature score showed a positive correlation with aggressive outcomes of HNSCC, such as low immunity score, significantly reduced CD8+ T cell infiltration and lowly expressed immunosuppressed biomarkers. However, high-risk groups of patients may have high chemotherapeutics sensitivity.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eConclusions:\u003c/strong\u003e The signature established by paring irlncRNA regardless of expression levels showed a promising clinical prediction value and revealed tumor immune microenvironment in HNSCC patients which might help in distinguishing those who could benefit from anti-tumor immunotherapy.\u003c/p\u003e","manuscriptTitle":"A Novel Immune-Related lncRNA Signature predict the Immune Landscape and Prognosis for Head and Neck Squamous Cell Carcinoma","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2021-02-22 21:03:50","doi":"10.21203/rs.3.rs-214493/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"c87f4a62-3540-4bad-b3b8-fffc24a60d3b","owner":[],"postedDate":"February 22nd, 2021","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":2538353,"name":"Cancer Biology"}],"tags":[],"updatedAt":"2021-02-22T21:03:52+00:00","versionOfRecord":[],"versionCreatedAt":"2021-02-22 21:03:50","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-214493","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-214493","identity":"rs-214493","version":["v1"]},"buildId":"WrCJVZZCHTDjtuVLN7oU0","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}