Novel prognostic model established for patients with head and neck squamous cell carcinoma based on pyroptosis-related genes

We aimed at establishing a risk – score model using pyroptosis-related genes to predict the prognosis of patients with head and neck squamous cell carcinoma (HNSCC). A total of 33 pyroptosis-related genes were selected. We then evaluated the data of 502 HNSCC patients and 44 normal patients from TCGA database. Gene expression was then profiled to detect differentially expressed genes (DEGs). Using the univariate, the least absolute shrinkage and selection operator (LASSO) Cox regression analyses, we generated a risk – score model. Tissue samples from neoplastic and normal sites of 44 HNSCC patients were collected. qRT-PCR were employed to analyze the mRNA level of the samples. Kaplan-Meier method was used to evaluate the overall survival rate (OS). Enrichment analysis was performed to elucidate the underlying mechanism of HNSCC patient's differentially survival status from the perspective of tumor immunology. 17 genes were categorized as DEGs. GSDME, IL-6, CASP8, CASP6, NLRP1 and NLRP6 were used to establish the risk – score model. Each patient's risk score in the TCGA cohort was calculated using the risk – score formula. The risk score was able to independently predict the OS of the HNSCC patients (P = 0.02). The OS analysis showed that the risk score model (P < 0.0001) was more reliable than single gene, a phenomenon verified by practical patient cohort. Additionally, enrichment analysis indicated more active immune activities in low-risk group than high-risk group. In conclusion, our risk – score model has provided novel strategy for the prediction of HNSCC patients’ prognosis.     

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.     

A robust six-miRNA prognostic signature for head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) remains a major health problem worldwide. We aimed to identify a robust microRNA (miRNA)-based signature for predicting HNSCC prognosis. The miRNA expression profiles of HNSCC were obtained from The Cancer Genome Atlas (TCGA) database. The TCGA HNSCC cohort was randomly divided into the discovery and validation cohort. A miRNA-based prognostic signature was built up based on TGCA discovery cohort, and then further validated. The downstream targets of prognostic miRNAs were subjected to functional enrichment analyses. The role of miR-1229-3p, a prognosis-related miRNA, in tumorigenesis of HNSCC was further evaluated. A total of 305 significantly differentially expressed miRNAs were found between HNSCC samples and normal tissues. A six-miRNA prognostic signature was constructed, which exhibited a strong association with overall survival (OS) in the TCGA discovery cohort. In addition, these findings were successfully confirmed in TCGA validation cohort and our own independent cohort. The miRNA-based signature was demonstrated as an independent prognostic indicator for HNSCC. A risk signature-based nomogram model was constructed and showed good performance for predicting the OS for HNSCC. The functional analyses revealed that the downstream targets of these prognostic miRNAs were closely linked to cancer progression. Mechanistically, in vitro analysis revealed that miR-1229-3p played a tumor promoting role in HNSCC. In conclusion, our study has developed a robust miRNA-based signature for predicting the prognosis of HNSCC with high accuracy, which will contribute to improve the therapeutic outcome.     

Identification of potential core genes and pathways predicting pathogenesis in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is the most common subtype of head and neck cancer; however, its pathogenesis and potential therapeutic targets remain largely unknown. In the present study, we analyzed three gene expression profiles and screened differentially expressed genes (DEGs) between HNSCC and normal tissues. The DEGs were subjected to gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG), protein–protein interaction (PPI), and survival analyses, while the connectivity map (CMap) database was used to predict candidate small molecules that may reverse the biological state of HNSCC. Finally, we measured the expression of the most relevant core gene in vitro and examined the effect of the top predicted potential drug against the proliferation of HNSCC cell lines. Among the 208 DEGs and ten hub genes identified, CDK1 and CDC45 were associated with unfavorable HNSCC prognosis, and three potential small molecule drugs for treating HNSCC were identified. Increased CDK1 expression was confirmed in HNSCC cells, and menadione, the top predicted potential drug, exerted significant inhibitory effects against HNSCC cell proliferation and markedly reversed CDK1 expression. Together, the findings of the present study suggest that the ten hub genes and pathways identified may be closely related to HNSCC pathogenesis. In particular, CDK1 and CDC45 overexpression could be reliable biomarkers for predicting unfavorable prognosis in patients with HNSCC, while the new candidate small molecules identified by CMap analysis provide new avenues for the development of potential drugs to treat HNSCC.     

Profiles of immune infiltration in lung adenocarcinoma and their clinical significant: A gene-expression–based retrospective study

Lung cancer is one of the fatal tumors. The tumor microenvironment plays a key role in regulating tumor progression. To figure out the role of tumor microenvironment in lung adenocarcinoma (LUAD), ESTIMATE algorithm was used to evaluate the immune scores in LUAD. Patients with low immune scores had a worse overall survival (OS) compared with high immune scores. Using RNA-Seq data of 489 patients in The Cancer Genome Atlas (TCGA), differentially expressed genes (DEGs) were identified between high- and low-immune score groups. Based on the DEGs, nine-gene signature was constructed by the least absolute shrinkage and selection operator Cox regression model in TCGA set. The signature demonstrated significant prognostic value in both TCGA and Gene Expression Omnibus database. Multivariate Cox regression analyses indicated that nine-genes signature was an independent prognostic factor. Subgroup analysis also revealed a robust prognostic ability of nine-gene signature. A nomogram with a C-index of 0.722 had a favorable power for predicting 3-, 5-, and 10-year survival for clinical use by integrating nine-gene signature and other clinical features. Co-expression and functional enrichment analysis showed that nine-gene signature was significantly associated with immune response and provided potential profound molecules for revealing the mechanism of tumor initiation and progression. In conclusion, we revealed the significance of immune infiltration and built a novel nine-gene signature as a reliable prognostic factor for patients with LUAD.     

DNA methylation biomarkers for head and neck squamous cell carcinoma

DNA methylation plays an important role in the etiology and pathogenesis of head and neck squamous cell carcinoma (HNSCC). The current study aimed to identify aberrantly methylated-differentially expressed genes (DEGs) by a comprehensive bioinformatics analysis. In addition, we screened for DEGs affected by DNA methylation modification and further investigated their prognostic values for HNSCC. We included microarray data of DNA methylation (GSE25093 and GSE33202) and gene expression (GSE23036 and GSE58911) from Gene Expression Omnibus. Aberrantly methylated-DEGs were analyzed with R software. The Cancer Genome Atlas (TCGA) RNA sequencing and DNA methylation (Illumina HumanMethylation450) databases were utilized for validation. In total, 27 aberrantly methylated genes accompanied by altered expression were identified. After confirmation by The Cancer Genome Atlas (TCGA) database, 2 hypermethylated-low-expression genes (FAM135B and ZNF610) and 2 hypomethylated-high-expression genes (HOXA9 and DCC) were identified. A receiver operating characteristic (ROC) curve confirmed the diagnostic value of these four methylated genes for HNSCC. Multivariate Cox proportional hazards analysis showed that FAM135B methylation was a favorable independent prognostic biomarker for overall survival of HNSCC patients.     

Identification of pyroptosis-related lncRNAs for constructing a prognostic model and their correlation with immune infiltration in breast cancer

The inflammasome-dependent cell death, which is denoted as pyroptosis, might be abnormally regulated during oncogenesis and tumour progression. Long non-coding RNAs (LncRNAs) are pivotal orchestrators in breast cancer (BC), which have the potential to be a biomarker for BC diagnosis and therapy. The present study aims to explore the correlation between pyroptosis-related lncRNAs and BC prognosis. In this study, a profile of 8 differentially expressed lncRNAs was screened in the TCGA database and used to construct a prognostic model. The BC patients were divided into high- and low-risk groups dependent on the median cutoff of the risk score in the model. Interestingly, the risk model significantly distinguished the clinical characteristics of BC patients between high- and low-risk groups. Then, the risk score of the model was identified to be an excellent independent prognostic factor. Notably, the GO, KEGG, GSEA and ssGSEA analyses revealed the different immune statuses between the high- and low-risk groups. Particularly, the 8 lncRNAs expressed differentially in BC tissues between two risk subgroups in vitro validation. Collectively, this constructed well-validated model is of high effectiveness to predict the prognosis of BC, which will provide novel means that is applicable for BC prognosis recognition.     

Construction of a prognostic signature in Ewing's sarcoma: Based on metabolism-related genes

ObjectiveBy combining the expression profiles of metabolism-related genes (MRGS) with clinical information, the expression quantities of MRGS and the influence on development and prognosis were systematically analyzed, so as to provide a theoretical basis for the clinical study on the prognosis of Ewing's sarcoma.MethodsMRGs expression profiles of 64 patients with Ewing's sarcoma were obtained from GEO dataset. Univariate Cox regression analysis was used to identify metabolization-related differentially expressed genes (DEGs) related with prognosis in Ewing's sarcoma patients. Then, multivariate Cox analysis was used to calculate novel prognostic markers based on metabolism-related DEGs. Besides, We validate the model using ICGC datasets. Finally, the new prognostic index was verified on the basis of the prognostic models.ResultsMultivariate Cox regression analysis identified 74 metabolization-related DEGs, 25 of which were associated with Ewing's sarcoma patients' overall survival. Subsequently, we used 25 DEGs to construct metabolism-related prognostic signature for patients with Ewing's sarcoma. Based on the 18 DEGs regression coefficient, we propose the formula of each patient's risk score, and then divided the patients into high-risk group and low-risk group. The results indicated that the survival rate and survival time were higher in the low-risk group and lower in the high-risk group. Multivariate Cox analysis showed that risk score index was an independent prognostic factor for Ewing's sarcoma.ConclusionThe experimental results suggest that the 18 metabolism-related DEGs marker may be effective in predicting the prognosis of Ewing's sarcoma to some extent, helping to individualize treatment of patients at different risks.     

Construction of a prognostic signature in Ewing's sarcoma: Based on metabolism-related genes

ObjectiveBy combining the expression profiles of metabolism-related genes (MRGS) with clinical information, the expression quantities of MRGS and the influence on development and prognosis were systematically analyzed, so as to provide a theoretical basis for the clinical study on the prognosis of Ewing's sarcoma.MethodsMRGs expression profiles of 64 patients with Ewing's sarcoma were obtained from GEO dataset. Univariate Cox regression analysis was used to identify metabolization-related differentially expressed genes (DEGs) related with prognosis in Ewing's sarcoma patients. Then, multivariate Cox analysis was used to calculate novel prognostic markers based on metabolism-related DEGs. Besides, We validate the model using ICGC datasets. Finally, the new prognostic index was verified on the basis of the prognostic models.ResultsMultivariate Cox regression analysis identified 74 metabolization-related DEGs, 25 of which were associated with Ewing's sarcoma patients' overall survival. Subsequently, we used 25 DEGs to construct metabolism-related prognostic signature for patients with Ewing's sarcoma. Based on the 18 DEGs regression coefficient, we propose the formula of each patient's risk score, and then divided the patients into high-risk group and low-risk group. The results indicated that the survival rate and survival time were higher in the low-risk group and lower in the high-risk group. Multivariate Cox analysis showed that risk score index was an independent prognostic factor for Ewing's sarcoma.ConclusionThe experimental results suggest that the 18 metabolism-related DEGs marker may be effective in predicting the prognosis of Ewing's sarcoma to some extent, helping to individualize treatment of patients at different risks.     

急性髓系白血病铁死亡相关基因筛选及预后模型的建立

基于急性髓系白血病(Acute Myeloid Leukemia,AML)临床大数据及多组学数据库探讨铁死亡相关基因在AML中的作用,并建立铁死亡基因表达相关预后模型。整合TCGA数据库中151例AML患者和GTEx数据库中337例正常人外周血的临床和转录组数据。将Wilcoxon检验和单因素Cox分析结果取交集,筛选出预后相关差异表达基因(Differential Expression Genes, DEGs),使用Lasso回归建立基因标志物预后模型,利用受试者工作特征曲线(Receiver Operating Characteristic Curve,ROC曲线)评价预测价值,Kaplan-Meier法进行生存分析,对AML患者临床数据进行单因素和多因素Cox回归分析,使用差异基因表达分析等方法比较高、低风险患者间的组学差异,最后,利用BeatAML数据库对基因标志物进行验证。将差异基因表达分析和单因素分析结果取交集,得到13个预后相关DEGs。构建了8个基因标志物的预后评分模型,并将患者分为高、低风险两组;ROC曲线分析证实了模型良好的预测性能;生存分析提示高、低风险组患者的生存率具有显著差异;单因素分析显示年龄和风险评分与患者整体生存显著相关,多因素分析显示,年龄和风险评分是独立预后指标。在2个风险组之间筛选出384个DEGs,GO富集分析结果显示,富集的基因大多与中性粒细胞和白细胞的趋化与迁移等免疫相关分子和通路显著相关,KEGG富集通路主要与TNF信号通路、细胞因子与细胞因子受体相互作用相关。BeatAML数据库验证结果显示,5个基因与预后显著相关。铁死亡相关基因在AML中显著表达,且高风险患者预后较差,该研究对AML铁死亡相关潜在生物标志物的发现和应用奠定了一定的基础。     

Correlation of zinc finger protein 2, a prognostic biomarker,with immune infiltrates in liver cancer

Purpose: The expression and clinical value of zinc finger protein 2 gene (ZIC2) in hepatocellular carcinoma (HCC) were analyzed by mining gene information from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases.Methods: Gene chip data sets were retrieved from GEO and TCGA and screened for differentially expressed genes in HCC. Gene expression profile interaction analysis (GEPIA) and Kaplan–Meier curves were used to analyze the relationship between differentially expressed genes (DEGs) and survival and prognosis in patients with HCC. Moreover, the Genecards database was used to extract ZIC2-related proteins and to analyze the physiological process of protein enrichment. Furthermore, the relationships between ZIC2 gene and tumor cell immune invasion and that between immune cell infiltration and the 5-year survival rate were studied using the tumor immune evaluation resource (TIMER) database.Results: Datasets from GEO and TCGA revealed that ZIC2 was differentially expressed in HCC tissues and normal tissues (P<0.05). High ZIC2 expression was associated with overall survival (OS) and progress-free survival in HCC patients. Overall, 25 ZIC2 related proteins, including Gli3, PRKDC, and rnf180 were identified and protein enrichment analysis indicated these were associated with four types of cell components, six types of cell functions, and eight types of biological processes. ZIC2 was positively correlated with immune infiltration cells in patients with HCC, and higher expression of ZIC2 mRNA CD4+T cells is associated with a better 5-year survival.Conclusion: ZIC2 gene may be used as an immune response marker in liver cancer to predict the prognosis of HCC.     

Identification and Validation of a Prognostic Model Based on Three MVI-Related Genes in Hepatocellular Carcinoma

MVI has significant clinical value for treatment selection and prognosis evaluation in hepatocellular carcinoma (HCC). We aimed to construct a model based on MVI-Related Genes (MVIRGs) for risk assessment and prognosis prediction in patients with HCC. This study utilized various statistical analysis methods for prognostic model construction and validation in the Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) cohorts, respectively. In addition, immunohistochemistry and qRT-PCR were used to analyze and identify the value of the model in our cohort. After the analyses, 153 differentially expressed MVIRGs were identified, and three key genes were selected to construct a prognostic model. The high-risk group showed significantly lower overall survival (OS), and this trend was observed in all subgroups: different age groups, genders, stages, and grades. Risk score was a risk factor independent of age, gender, stage, and grade. Moreover, the ICGC cohort validated the prognostic value of the model corresponding to the TCGA. In our cohort, qRT-PCR and immunohistochemistry showed that all three genes had higher expression levels in HCC samples than in normal controls. High expression levels of genes and high-risk scores showed significantly lower recurrence-free survival (RFS) and OS, especially in MVI-positive HCC samples. Therefore, the prognostic model constructed by three MVIRGs can reliably predict the RFS and OS of patients with HCC and is valuable for guiding clinical treatment selection and prognostic assessment of HCC.     

Identification of a six-gene prognostic signature for oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is one of the most common types of malignancies worldwide, and its morbidity and mortality have increased in the near term. Consequently, the purpose of the present study was to identify the notable differentially expressed genes (DEGs) involved in their pathogenesis to obtain new biomarkers or potential therapeutic targets for OSCC. The gene expression profiles of the microarray datasets GSE85195, GSE23558, and GSE10121 were obtained from the Gene Expression Omnibus (GEO) database. After screening the DEGs in each GEO dataset, 249 DEGs in OSCC tissues were obtained. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathway enrichment analysis was employed to explore the biological functions and pathways of the above DEGs. A protein–protein interaction network was constructed to obtain a central gene. The corresponding total survival information was analyzed in patients with oral cancer from The Cancer Genome Atlas (TCGA). A total of six candidate genes (CXCL10, OAS2, IFIT1, CCL5, LRRK2, and PLAUR) closely related to the survival rate of patients with oral cancer were identified, and expression verification and overall survival analysis of six genes were performed based on TCGA database. Time-dependent receiver operating characteristic curve analysis yields predictive accuracy of the patient's overall survival. At the same time, the six genes were further verified by quantitative real-time polymerase chain reaction using samples obtained from the patients recruited to the present study. In conclusion, the present study identified the prognostic signature of six genes in OSCC for the first time via comprehensive bioinformatics analysis, which could become potential prognostic markers for OCSS and may provide potential therapeutic targets for tumors.     

Integrated analysis of gene expression and DNA methylation datasets identified key genes and a 6-gene prognostic signature for primary lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the main subtype of non-small cell lung cancer with a poor survival prognosis. In our study, gene expression, DNA methylation, and clinicopathological data of primary LUAD were utilized to identify potential prognostic markers for LUAD, which were recruited from The Cancer Genome Atlas (TCGA) database. Univariate regression analysis showed that there were 21 methylation-associated DEGs related to overall survival (OS), including 9 down- and 12 up-regulated genes. The 12 up-regulated genes with hypomethylation may be risky genes, whereas the other 9 down-regulated genes with hypermethylation might be protective genes. By using the Step-wise multivariate Cox analysis, a methylation-associated 6-gene (consisting of CCL20, F2, GNPNAT1, NT5E, B3GALT2, and VSIG2) prognostic signature was constructed and the risk score based on this gene signature classified patients into high- or low-risk groups. Patients of the high-risk group had shorter OS than those of the low-risk group in both the training and validation cohort. Multivariate Cox analysis and the stratified analysis revealed that the risk score was an independent prognostic factor for LUAD patients. The methylation-associated gene signature may serve as a prognostic factor for LUAD patients and the represent hypermethylated or hypomethylated genes might be potential targets for LUAD therapy.     

Identification of immune-related hub genes and construction of an immune-related gene prognostic index for low-grade glioma

Low-grade glioma (LGG) poses significant management challenges and has a dismal prognosis. While immunotherapy has shown significant promise in cancer treatment, its progress in glioma has confronted with challenges. In our study, we aimed to develop an immune-related gene prognostic index (IRGPI) which could be used to evaluate the response and efficacy of LGG patients with immunotherapy. We included a total of 529 LGG samples from TCGA database and 1152 normal brain tissue samples from the GTEx database. Immune-related differentially expressed genes (DEGs) were screened. Then, we used weighted gene co-expression network analysis (WGCNA) to identify immune-related hub genes in LGG patients and performed Cox regression analysis to construct an IRGPI. The median IRGPI was used as the cut-off value to categorize LGG patients into IRGPI-high and low subgroups, and the molecular and immune mechanism in IRGPI-defined subgroups were analysed. Finally, we explored the relationship between IRGPI-defined subgroups and immunotherapy related indicators in patients after immunotherapy. Three genes (RHOA, NFKBIA and CCL3) were selected to construct the IRGPI. In a survival analysis using TCGA cohort as a training set, patients in the IRGPI-low subgroup had a better OS than those in IRGPI-high subgroup, consistent with the results in CGGA cohort. The comprehensive results showed that IRGPI-low subgroup had a more abundant activated immune cell population and lower TIDE score, higher MSI, higher TMB score, lower T cell dysfunction score, more likely benefit from ICIs therapy. IRGPI is a promising biomarker in the field of LGG ICIs therapy to distinguish the prognosis, the molecular and immunological characteristics of patients.     

Survival prediction in patients with colon adenocarcinoma via multiomics data integration using a deep learning algorithm

The present study proposed a deep learning (DL) algorithm to predict survival in patients with colon adenocarcinoma (COAD) based on multiomics integration. The survival-sensitive model was constructed using an autoencoder for DL implementation based on The Cancer Genome Atlas (TCGA) data of patients with COAD. The autoencoder framework was compared with PCA, NMF, t-SNE, and univariable Cox-PH model for identifying survival-related features. The prognostic robustness of the inferred survival risk groups was validated using three independent confirmation cohorts. Differential expression analysis, Pearson’s correlation analysis, construction of miRNA–target gene network, and function enrichment analysis were performed. Two risk groups with significant survival differences were identified in TCGA set using the autoencoder-based model (log-rank P-value = 5.51e⁻⁰⁷). The autoencoder framework showed superior performance compared with PCA, NMF, t-SNE, and the univariable Cox-PH model based on the C-index, log-rank P-value, and Brier score. The robustness of the classification model was successfully verified in three independent validation sets. There were 1271 differentially expressed genes, 10 differentially expressed miRNAs, and 12 hypermethylated genes between the survival risk groups. Among these, miR-133b and its target genes (GNB4, PTPRZ1, RUNX1T1, EPHA7, GPM6A, BICC1, and ADAMTS5) were used to construct a network. These genes were significantly enriched in ECM–receptor interaction, focal adhesion, PI3K–Akt signaling pathway, and glucose metabolism-related pathways. The risk subgroups obtained through a multiomics data integration pipeline using the DL algorithm had good robustness. miR-133b and its target genes could be potential diagnostic markers. The results would assist in elucidating the possible pathogenesis of COAD.     

基于肿瘤免疫微环境鉴定IDH突变型弥漫性胶质瘤的预后标志物

异柠檬酸脱氢酶(IDH)突变存在于大多数低级别胶质细胞瘤中,免疫逃逸是肿瘤标志性特征之一,免疫治疗在胶质瘤的治疗中的作用越来越重要。利用生信手段分析TCGA、CGGA、GEO数据集中IDH突变胶质瘤的770个免疫相关基因及其临床相关数据,从而获得每个患者的免疫风险评分(IMRS);结合IMRS和临床信息,筛选出6个差异表达基因(TRAF3、ATG10、BID、TAB1、MAP3K1、RPS6)组成IMRS模型并生成诺莫图对患者预后进行评估,发现低风险组患者的总生存期(OS)较高风险组均明显延长。此外,签名相关免疫细胞浸润分析发现肿瘤相关巨噬细胞浸润评分(TAM)与肿瘤相关T细胞浸润评分(TIS)呈明显的负相关,表明高IMRS富集了促肿瘤免疫浸润,而低IMRS则富集了相对较多的抗肿瘤免疫浸润。     

Outlier Analysis Defines Zinc Finger Gene Family DNA Methylation in Tumors and Saliva of Head and Neck Cancer Patients

Head and Neck Squamous Cell Carcinoma (HNSCC) is the fifth most common cancer, annually affecting over half a million people worldwide. Presently, there are no accepted biomarkers for clinical detection and surveillance of HNSCC. In this work, a comprehensive genome-wide analysis of epigenetic alterations in primary HNSCC tumors was employed in conjunction with cancer-specific outlier statistics to define novel biomarker genes which are differentially methylated in HNSCC. The 37 identified biomarker candidates were top-scoring outlier genes with prominent differential methylation in tumors, but with no signal in normal tissues. These putative candidates were validated in independent HNSCC cohorts from our institution and TCGA (The Cancer Genome Atlas). Using the top candidates, ZNF14, ZNF160, and ZNF420, an assay was developed for detection of HNSCC cancer in primary tissue and saliva samples with 100% specificity when compared to normal control samples. Given the high detection specificity, the analysis of ZNF DNA methylation in combination with other DNA methylation biomarkers may be useful in the clinical setting for HNSCC detection and surveillance, particularly in high-risk patients. Several additional candidates identified through this work can be further investigated toward future development of a multi-gene panel of biomarkers for the surveillance and detection of HNSCC.