Identification of a prognostic gene signature based on an immunogenomic landscape analysis of bladder cancer

Cancer immune plays a critical role in cancer progression. Tumour immunology and immunotherapy are one of the exciting areas in bladder cancer research. In this study, we aimed to develop an immune‐related gene signature to improve the prognostic prediction of bladder cancer. Firstly, we identified 392 differentially expressed immune‐related genes (IRGs) based on TCGA and ImmPort databases. Functional enrichment analysis revealed that these genes were enriched in inflammatory and immune‐related pathways, including in ‘regulation of signaling receptor activity’, ‘cytokine‐cytokine receptor interaction’ and ‘GPCR ligand binding’. Then, we separated all samples in TCGA data set into the training cohort and the testing cohort in a ratio of 3:1 randomly. Data set {"type":"entrez-geo","attrs":{"text":"GSE13507","term_id":"13507"}}GSE13507 was set as the validation cohort. We constructed a prognostic six‐IRG signature with LASSO Cox regression in the training cohort, including AHNAK, OAS1, APOBEC3H, SCG2, CTSE and KIR2DS4. Six IRGs reflected the microenvironment of bladder cancer, especially immune cell infiltration. The prognostic value of six‐IRG signature was further validated in the testing cohort and the validation cohort. The results of multivariable Cox regression and subgroup analysis revealed that six‐IRG signature was a clinically independent prognostic factor for bladder cancer patients. Further, we constructed a nomogram based on six‐IRG signature and other clinicopathological risk factors, and it performed well in predict patients'' survival. Finally, we found six‐IRG signature showed significant difference in different molecular subtypes of bladder cancer. In conclusions, our research provided a novel immune‐related gene signature to estimate prognosis for patients'' survival with bladder cancer.     

Identification of a four immune-related genes signature based on an immunogenomic landscape analysis of clear cell renal cell carcinoma

Renal clear cell carcinoma (ccRCC) is the most common type of renal cell carcinoma, which has strong immunogenicity. A comprehensive study of the role of immune-related genes (IRGs) in ccRCC is of great significance in finding ccRCC treatment targets and improving patient prognosis. In this study, we comprehensively analyzed the expression of IRGs in ccRCC based on The Cancer Genome Atlas datasets. The mechanism of differentially expressed IRGs in ccRCC was analyzed by bioinformatics. In addition, Cox regression analysis was used to screen prognostic related IRGs from differentially expressed IRGs. We also identified a four IRGs signature consisting of four IRGs (CXCL2, SEMA3G, PDGFD, and UCN) through lasso regression and multivariate Cox regression analysis. Further analysis results showed that the four IRGs signature could effectively predict the prognosis of patients with ccRCC, and its predictive power is independent of other clinical factors. In addition, the correlation analysis of immune cell infiltration showed that this four IRGs signature could effectively reflect the level of immune cell infiltration of ccRCC. We also found that the expression of immune checkpoint genes CTLA-4, LAG3, and PD-1 in the high-risk group was higher than that in the low-risk group. Our research revealed the role of IRGs in ccRCC, and developed a four IRGs signature that could be used to evaluate the prognosis of patients with ccRCC, which will help to develop personalized treatment strategies for patients with ccRCC and improve their prognosis. In addition, these four IRGs may be effective therapeutic targets for ccRCC.     

Coordinated loading of IRG resistance GTPases on to the Toxoplasma gondii parasitophorous vacuole

The immunity‐related GTPases (IRGs) constitute an interferon‐induced intracellular resistance mechanism in mice against Toxoplasma gondii. IRG proteins accumulate on the parasitophorous vacuole membrane (PVM), leading to its disruption and to death of the parasite. How IRGs target the PVM is unknown. We show that accumulation of IRGs on the PVM begins minutes after parasite invasion and increases for about 1 h. Targeting occurs independently of several signalling pathways and the microtubule network, suggesting that IRG transport is diffusion‐driven. The intensity of IRG accumulation on the PVM, however, is reduced in absence of the autophagy regulator, Atg5. In wild‐type cells IRG proteins accumulate cooperatively on PVMs in a definite order reflecting a temporal hierarchy, with Irgb6 and Irgb10 apparently acting as pioneers. Loading of IRG proteins onto the vacuoles of virulent Toxoplasma strains is attenuated and the two pioneer IRGs are the most affected. The polymorphic rhoptry kinases, ROP16, ROP18 and the catalytically inactive proteins, ROP5A–D, are not individually responsible for this effect. Thus IRG proteins protect mice against avirulent strains of Toxoplasma but fail against virulent strains. The complex cooperative behaviour of IRG proteins in resisting Toxoplasma may hint at undiscovered complexity also in virulence mechanisms.     

Identification of a five-gene signature with prognostic value in colorectal cancer

Colorectal cancer (CRC) ranks as one of the most commonly diagnosed malignancies worldwide. Although mortality rates have been decreasing, the prognosis of CRC patients is still highly dependent on the individual. Therefore, identifying and understanding novel biomarkers for CRC prognosis remains crucial. The gene expression profiles of five-gene expression omnibus (GEO) data sets of CRC were first downloaded. A total of 352 consistent differentially expressed genes (DEGs) were identified for CRC and paired with normal tissues. Functional analysis including gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment revealed that these DEGs were related to metabolic pathways, tight junctions, and the cell cycle. Ten hub DEGs were identified based on the search tool for the retrieval of interacting genes database and protein–protein interaction networks. By using univariate Cox proportional hazard regression analysis, we found 11 survival-related genes among these DEGs. We finally established a five-gene signature (kinesin family member 15, N-acetyltransferase 2, glutathione peroxidase 3, secretogranin II, and chloride channel accessory 1) with prognostic value in CRC by step multivariate Cox regression analysis. Based on this risk scoring system, patients in the high-risk group had significantly poorer survival results compared with those in the low-risk group (log-rank test, p < 0.0001). Finally, we validated our gene signature scoring system in two independent GEO cohorts (GSE17536 and GSE33113). We found all five of the signature genes to be DEGs in The Cancer Genome Atlas database. In conclusion, our findings suggest that our five DEG-based signature can provide a novel biomarker with useful applications in CRC prognosis.     

An eight-long noncoding RNA expression signature for colorectal cancer patients’ prognosis

Long noncoding RNAs (lncRNAs) have recently emerged as important biomarkers of cancer progression. Here, we proposed to develop a lncRNA-based signature with a prognostic value for colorectal cancer (CRC) overall survival (OS). Through mining microarray datasets, we analyzed the lncRNA expression profiles of 122 patients with CRC from Gene Expression Omnibus. Associations between lncRNA and CRC OS were firstly evaluated through univariate Cox regression analysis. A random survival forest method was applied for further screening of the lncRNA signature, which resulted in eight lncRNAs, including PEG3-AS1, LOC100505715, MINCR, DBH-AS1, LINC00664, FAM224A, LOC642852, and LINC00662. Combination of the eight lncRNAs weighted by their multivariate Cox regression coefficients formed a prognostic signature, through which, we could divide the 122 patients with CRC into two subgroups with significantly different OS. Good robustness of the lncRNA signature's prognostic value was verified through an independent data set consisting of 55 patients with CRC. In addition, gene set enrichment analysis indicated the potential association between high prognostic value and oxygen metabolism-related processes. This result should indicate that lncRNAs could be a useful signature for CRC prognosis.     

Profiles of immune-related genes and immune cell infiltration in the tumor microenvironment of diffuse lower-grade gliomas

The tumor microenvironment is highly correlated with tumor occurrence, progress, and prognosis. We aimed to investigate the immune-related gene (IRG) expression and immune infiltration pattern in the tumor microenvironment of lower-grade glioma (LGG). We employed the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm to calculate immune and stromal scores and identify prognostic IRG based on The Cancer Genome Atlas data set. The potential molecular functions of these genes were explored with the help of functional enrichment analysis and the protein–protein interaction network. Remarkably, three cohorts that were downloaded from the Chinese Glioma Genome Atlas database were analyzed to further verify the prognostic values of these genes. Moreover, the Tumor IMmune Estimation Resource (TIMER) algorithm was used to estimate the abundance of infiltrating immune cells and explore the immune infiltration pattern in LGG. And unsupervised cluster analysis determined three clusters of the immune infiltration pattern and indicated that CD8⁺ T cells and macrophages were significantly associated with LGG outcomes. Altogether, our study identified a list of prognostic IRGs and provided a perspective to explore the immune infiltration pattern in LGG.     

A dual immune signature of CD8+ T cells and MMP9 improves the survival of patients with hepatocellular carcinoma

The 5-year survival of hepatocellular carcinoma (HCC) is difficult due to the high recurrence rate and metastasis. Tumor infiltrating immune cells (TICs) and immune-related genes (IRGs) bring hope to improve survival and treatment of HCC patients. However, there are problems in predicting immune signatures and identifying novel therapeutic targets. In the study, the CIBERSORT algorithm was used to evaluate 22 immune cell infiltration patterns in gene expression omnibus (GEO) and the cancer genome atlas (TCGA) data. Eight immune cells were found to have significant infiltration differences between the tumor and normal groups. The CD8+ T cells immune signature was constructed by least absolute shrinkage and selection operator (LASSO) algorithm. The high infiltration level of CD8+ T cells could significantly improve survival of patients. The weighted gene co-expression network analysis (WGCNA) algorithm identified MMP9 was closely related to the overall survival of HCC patients. K-M survival and tROC analysis confirmed that MMP9 had an excellent prognostic prediction. Cox regression showed that a dual immune signature of CD8+ T cells and MMP9 was independent survival factor in HCC. Therefore, a dual prognostic immune signature could improve the survival of patient and may provide a new strategy for the immunotherapy of HCC.     

Development of a novel gene signature in patients without Helicobacter pylori infection gastric cancer

Gastric cancer (GC) is one of the most fatal common cancers in worldwide. Helicobacter pylori (H. pylori) infection is closely related to the development of GC, although the mechanism is still unclear. In our study, we aim to develop a robust messenger RNA (mRNA) signature associated with H. pylori (-) GC that can sensitively and efficiently predict the prognostic. The RNA-seq expression profile and corresponding clinical data of 598 gastric cancer samples and 63 normal samples obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus database. Using gene set enrichment analysis H. pylori (+) GC and H. pylori (-) GC patients and normal samples to select certain genes for further analysis. Using univariate and multivariate Cox regression model to establish a gene signature for predicting the overall survival (OS). Finally, we identified G2/M related seven-mRNA signature (TGFB1, EGF, MKI67, ILF3, INCENP, TNPO2, and CHAF1A) closely related to the prognosis of patients with H. pylori (-) GC. The seven-mRNA signature was identified to act as an independent prognostic biomarker by stratified analysis and multivariate Cox regression analysis. It was also validated on two test groups from TCGA and GSE15460 and shown that patients with high-risk scores based on the expression of the seven mRNAs had significantly shorter survival times compared to patients with low-risk scores (P < .0001). In this study, we developed a seven-mRNA signature related to G2/M checkpoint from H. pylori (-) GCs that as an independent biomarker potentially with a good performance in predicting OS and might be valuable for the clinical management for patients with GC.     

In vitro and clinical data analysis of Osteopontin as a prognostic indicator in colorectal cancer

Osteopontin (OPN) has been shown to promote colorectal cancer (CRC) progression; however, the mechanism of OPN‐induced CRC progression is largely unknown. In this study, we found that OPN overexpression led to enhanced anchorage‐independent growth, cell migration and invasion in KRAS gene mutant cells but to a lesser extent in KRAS wild‐type cells. OPN overexpression also induced PI3K signalling, expression of Snail and Matrix metallopeptidase 9 (MMP9), and suppressed the expression of E‐cadherin in KRAS mutant cells. In human CRC specimens, a high‐level expression of OPN significantly predicted poorer survival in CRC patients and OPN expression was positively correlated with MMP9 expression, and negatively correlated with E‐cadherin expression. Furthermore, we have found that 15 genes were co‐upregulated in OPN highly expression CRC and a list of candidate drugs that may have potential to reverse the secreted phosphoprotein 1 (SPP1) gene signature by connectivity mapping. In summary, OPN is a potential prognostic indicator and therapeutic target for colon cancer.     

Construction of an N6-methyladenosine lncRNA- and immune cell infiltration-related prognostic model in colorectal cancer

The present paper aims to shed light on the influence of N6-methyladenosine (m6A) long non-coding RNAs (lncRNAs) and immune cell infiltration on colorectal cancer (CRC). We downloaded workflow-type data and xml-format clinical data on CRC from The Cancer Genome Atlas project. The relationship between lncRNA and m6A was identified by using Perl and R software. Kyoto Encyclopedia of Genes and Genomes enrichment analysis was performed. Lasso regression was utilized to construct a prognostic model. Survival analysis was used to explore the relationship between clusters of m6A lncRNAs and clinical survival data. Differential analysis of the tumor microenvironment and an immune correlation analysis were used to determine immune cell infiltration levels in different clusters and their correlation with clinical prognosis. The expression of lncRNA was tightly associated with m6A. The univariate Cox regression analysis showed that lncRNA was a risk factor for the prognosis. Differential expression analysis demonstrated that m6A lncRNAs were partially highly expressed in tumor tissue. m6A lncRNA-related prognostic model could predict the prognosis of CRC independently. “ECM_RECEPTOR_INTERACTION” was the most significantly enriched gene set. PARP8 was overexpressed in tumor tissue and high-risk cluster. CD4 memory T cells, activated resting NK cells, and memory B cells were highly clustered in the high-risk cluster. All of the scores were higher in the low-risk group. m6A lncRNA is closely related to the occurrence and progression of CRC. The corresponding prognostic model can be utilized to evaluate the prognosis of CRC. m6A lncRNA and related immune cell infiltration in the tumor microenvironment can provide novel therapeutic targets for further research.

     

Identification of a novel cell cycle-related gene signature predicting survival in patients with gastric cancer

Gastric cancer (GC) is one of the most fatal cancers in the world. Thousands of biomarkers have been explored that might be related to survival and prognosis via database mining. However, the prediction effect of single gene biomarkers is not specific enough. Increasing evidence suggests that gene signatures are emerging as a possible better alternative. We aimed to develop a novel gene signature to improve the prognosis prediction of GC. Using the messenger RNA (mRNA)-mining approach, we performed mRNA expression profiling in a large GC cohort (n = 375) from The Cancer Genome Atlas (TCGA) database. Gene Set Enrichment Analysis (GSEA) was performed, and we recovered genes related to the G2/M checkpoint, which we identified with a Cox proportional regression model. We identified a set of five genes (MARCKS, CCNF, MAPK14, INCENP, and CHAF1A), which were significantly associated with overall survival (OS) in the test series. Based on this five-gene signature, the test series patients could be classified into high-risk or low-risk subgroups. Multivariate Cox regression analysis indicated that the prognostic power of this five-gene signature was independent of clinical features. In conclusion, we developed a five-gene signature related to the cell cycle that can predict survival for GC. Our findings provide novel insight that is useful for understanding cell cycle mechanisms and for identifying patients with GC with poor prognoses.     

Identification of m6A methyltransferase-related lncRNA signature for predicting immunotherapy and prognosis in patients with hepatocellular carcinoma

N6-methyladenosine (m6A) methyltransferase has been shown to be an oncogene in a variety of cancers. Nevertheless, the relationship between the long non-coding RNAs (lncRNAs) and hepatocellular carcinoma (HCC) remains elusive. We integrated the gene expression data of 371 HCC and 50 normal tissues from The Cancer Genome Atlas (TCGA) database. Differentially expressed protein-coding genes (DE-PCGs)/lncRNAs (DE-lncRs) analysis and univariate regression and Kaplan–Meier (K–M) analysis were performed to identify m6A methyltransferase-related lncRNAs. Three prognostic lncRNAs were selected by univariate and LASSO Cox regression analyses to construct the m6A methyltransferase-related lncRNA signature. Multivariate Cox regression analyses illustrated that this signature was an independent prognostic factor for overall survival (OS) prediction. The Gene Set Enrichment Analysis (GSEA) suggested that the m6A methyltransferase-related lncRNAs were involved in the immune-related biological processes (BPs) and pathways. Besides, we discovered that the lncRNAs signature was correlated with the tumor microenvironment (TME) and the expression of critical immune checkpoints. Tumor Immune Dysfunction and Exclusion (TIDE) analysis revealed that the lncRNAs could predict the clinical response to immunotherapy. Our study had originated a prognostic signature for HCC based on the potential prognostic m6A methyltransferase-related lncRNAs. The present study had deepened the understanding of the TME status of HCC patients and laid a theoretical foundation for the choice of immunotherapy.     

Five genes as a novel signature for predicting the prognosis of patients with laryngeal cancer

In this study, we purpose to investigate a novel five-gene signature for predicting the prognosis of patients with laryngeal cancer. The laryngeal cancer datasets were obtained from The Cancer Genome Atlas (TCGA). Both univariate and multivariate Cox regression analysis was applied to screening for prognostic differential expressed genes (DEGs), and a novel gene signature was obtained. The performance of this Cox regression model was tested by receiver operating characteristic (ROC) curves and area under the curve (AUC). Further survival analysis for each of the five genes was carried out through the Kaplan-Meier curve and Log-rank test. Totally, 622 DEGs were screened from the TCGA datasets in this study. We construct a five-gene signature through Cox survival analysis. Patients were divided into low- and high-risk groups depending on the median risk score, and a significant difference of the 5-year overall survival was found between these two groups (P < .05). ROC curves verified that this five-gene signature had good performance to predict the prognosis of laryngeal cancer (AUC = 0.862, P < .05). In conclusion, the five-gene signature consist of EMP1, HOXB9, DPY19L2P1, MMP1, and KLHDC7B might be applied as an independent prognosis predictor of laryngeal cancer.     

Expression signature of six‐snoRNA serves as novel non‐invasive biomarker for diagnosis and prognosis prediction of renal clear cell carcinoma

Increasing evidence has verified that small nucleolar RNAs (snoRNAs) play significant roles in tumorigenesis and exhibit prognostic value in clinical practice. In the study, we analysed the expression profile and clinical relevance of snoRNAs from TCGA database including 530 ccRCC (clear cell renal cell carcinoma) and 72 control cases. By using univariate and multivariate Cox analysis, we established a six‐snoRNA signature and divided patients into high‐risk or low‐risk groups. We found patients in high‐risk group had significantly shorter overall survival and recurrence‐free survival than those in low‐risk group in test series, validation series and entire series by Kaplan‐Meier analysis. We also confirmed this signature had a great accuracy and specificity in 64 clinical tissue cases and 50 serum samples. Then, depending on receiver operating characteristic curve analysis we found the six‐snoRNA signature was an superior indicator better than conventional clinical factors (AUC = 0.732). Furthermore, combining the signature with TNM stage or Fuhrman grade were the optimal indicators (AUC = 0.792; AUC = 0.800) and processed the clinical applied value for ccRCC. Finally, we found the SNORA70B and its hose gene USP34 might directly regulate Wnt signalling pathway to promote tumorigenesis in ccRCC. In general, our study established a six‐snoRNA signature as an independent and superior diagnosis and prognosis indicator for ccRCC.     

Neurexophilin and PC-esterase domain family member 4 (NXPE4) and prostate androgen-regulated mucin-like protein 1 (PARM1) as prognostic biomarkers for colorectal cancer

Owing to the high morbidity and mortality, novel biomarkers in the occurrence and development of colorectal cancer (CRC) are needed nowadays. In this study, the CRC-related datasets were downloaded from the Gene Expression Omnibus (GEO) database and The Cancer Genome Atlas (TCGA) database. After screening the differentially expressed genes (DEGs) in R software, a total of 238 upregulated and 199 downregulated DEGs were revealed simultaneously. Then the Kaplan-Meier survival analysis and Cox regression analysis were used to reveal the prognostic function of these DEGs. Neurexophilin and PC-esterase domain family member 4 (NXPE4) and prostate androgen-regulated mucin-like protein 1 (PARM1) were two outstanding independent overall survival (OS) and relapse-free survival (RFS) prognostic genes of CRC in TCGA database. We next verified the expression of NXPE4 and PARM1 messenger RNA (mRNA) levels were significantly lower in CRC tumor tissue than in the adjacent noncancerous tissue in our clinical samples, and NXPE4 mRNA expression level was related to the tumor location and tumor size, while PARM1 was related to tumor location, lymph nodes metastasis, and tumor size. This study demonstrated that NXPE4 and PARM1 might be two potential novel prognostic biomarkers for CRC.     

A novel LncRNA-based prognostic score reveals TP53-dependent subtype of lung adenocarcinoma with poor survival

The prognostic signatures play an essential role in the era of personalised therapy for cancer patients including lung adenocarcinoma (LUAD). Long noncoding RNA (LncRNA), a relatively novel class of RNA, has shown to play a crucial role in all the areas of cancer biology. Here, we developed and validated a robust LncRNA-based prognostic signature for LUAD patients using three different cohorts. In the discovery cohort, four LncRNAs were identified with 10% false discovery rate and a hazard ratio of >10 using univariate Cox regression analysis. A risk score, generated from the four LncRNAs’ expression, was found to be a significant predictor of survival in the discovery and validation cohort (p = 9.97 × 10 ⁻⁸ and 1.41 × 10 ⁻³, respectively). Further optimisation of four LncRNAs signature in the validation cohort, generated a three LncRNAs prognostic score (LPS), which was found to be an independent predictor of survival in both the cohorts ( p = 1.00 × 10 ⁻⁶ and 7.27 × 10 ⁻⁴, respectively). The LPS also significantly divided survival in clinically important subsets, including Stage I ( p = 9.00 × 10 ⁻⁴ and 4.40 × 10 ⁻², respectively), KRAS wild-type (WT), KRAS mutant ( p = 4.00 × 10 ⁻³ and 4.30 × 10 ⁻², respectively) and EGFR WT ( p = 2.00 × 10 ⁻⁴). In multivariate analysis LPS outperformed, eight previous prognosticators. Further, individual members of LPS showed a significant correlation with survival in microarray data sets. Mutation analysis showed that high-LPS patients have a higher mutation rate and inactivation of the TP53 pathway. In summary, we identified and validated a novel LncRNA signature LPS for LUAD.     

A novel prognostic immunoscore based on The Cancer Genome Atlas to predict overall survival in colorectal cancer patients

Colorectal cancer (CRC) is highly prevalent worldwide. The relationship between the infiltration of immunocytes in CRC and clinical outcome has been investigated in recent years. The present study aims to construct a new prognostic signature using an immunocyte panel. Our novel prognostic immunoscore included 13 types of immunocytes, which were identified by least absolute shrinkage and selection operator (LASSO)-Cox regression. The time-dependent receiver operating characteristic (ROC) curve and Kaplan–Meier survival estimates were applied to evaluate the prognostic ability. Compared with the signature based on a single immune marker (i.e., CD8 mRNA expression and CD8⁺ expressing T cells), the novel prognostic immunoscore possessed better specificity and sensitivity of prognosis (area under the curves (AUCs) are 0.852, 0.856, and 0.774 for 1-, 2-, and 3-year survival times, respectively). Significant differences were identified between the high and low immunoscore groups in overall survival and disease-free survival in training and validation cohorts. Combining the immunoscore with clinical information may provide a more accurate prognosis for CRC. The immunoscore can identify patients with poor outcomes in the high Tumor Mutational Burden (TMB) group, who may benefit the most from immunotherapy. The immunoscore was also closely related to two immune checkpoints (i.e., PD-L1 and PD-1, r = 0.3087 and r = 0.3341, respectively). Collectively, our study demonstrates that the novel prognostic immunoscore reported here may be useful in distinguishing different prognoses and may improve the clinical management of patients with CRC.     

Identification and validation of autophagy-related prognostic signature for head and neck squamous cell carcinoma

BackgroundMany studies have demonstrated that autophagy plays a significant role in regulating tumor growth and progression. However, the effect of autophagy-related genes (ARGs) on the prognosis have rarely been analyzed in head and neck squamous cell carcinoma (HNSCC).MethodsWe obtained differentially expressed ARGs from HNSCC mRNA data in The Cancer Genome Atlas (TCGA) database. And then we performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses to explore the autophagy-related biological functions. The overall survival (OS)-related and disease specific survival (DSS)-related ARGs were identified by univariate Cox regression analyses. With these genes, we established OS-related and DSS-related risk signature by LASSO regression method, respectively. We validated the reliability of the risk signature with receiver operating characteristic (ROC) analysis, Kaplan-Meier survival curves, clinical correlation analysis, and nomogram. Then we analyzed relationships between risk signature and immune cell infiltration.ResultsWe established the prognostic signatures based on 14 ARGs for OS and 12 ARGs for DSS. The ROC curves, survival analysis, and nomogram validated the predictive accuracy of the models. Clinic correlation analysis showed that the risk group was closely related to Stage, pathological T stage, pathological N stage and human papilloma virus (HPV) subtype. Cox regression demonstrated that the risk score was an independent predictor for the prognosis of HNSCC patients. Furthermore, patients in low-risk score group exhibited higher immunescore and distinct immune cell infiltration than high-risk score group. And we further analysis revealed that the copy number alterations (CNAs) of ARGs-based signature affected the abundance of tumor-infiltrating immune cells.ConclusionIn this study, we identified novel autophagy-related signature for the prediction of OS and DSS in patients with HNSCC. Meanwhile, our study provides a novel sight to understand the role of autophagy and elucidate the important role of autophagy in tumor immune microenvironment (TIME) of HNSCC.