Ferroptosis-related gene signature correlates with the tumor immune features and predicts the prognosis of glioma patients

Background: Glioma is a malignant intracranial tumor and the most fatal cancer. The role of ferroptosis in the clinical progression of gliomas is unclear.Method: Univariate and least absolute shrinkage and selection operator (Lasso) Cox regression methods were used to develop a ferroptosis-related signature (FRSig) using a cohort of glioma patients from the Chinese Glioma Genome Atlas (CGGA), and was validated using an independent cohort of glioma patients from The Cancer Genome Atlas (TCGA). A single-sample gene set enrichment analysis (ssGSEA) was used to calculate levels of the immune infiltration. Multivariate Cox regression was used to determine the independent prognostic role of clinicopathological factors and to establish a nomogram model for clinical application.Results: We analyzed the correlations between the clinicopathological features and ferroptosis-related gene (FRG) expression and established an FRSig to calculate the risk score for individual glioma patients. Patients were stratified into two subgroups with distinct clinical outcomes. Immune cell infiltration in the glioma microenvironment and immune-related indexes were identified that significantly correlated with the FRSig, the tumor mutation burden (TMB), copy number alteration (CNA), and immune checkpoint expression was also significantly positively correlated with the FRSig score. Ultimately, an FRSig-based nomogram model was constructed using the independent prognostic factors age, World Health Organization (WHO) grade, and FRSig score.Conclusion: We established the FRSig to assess the prognosis of glioma patients. The FRSig also represented the glioma microenvironment status. Our FRSig will contribute to improve patient management and individualized therapy by offering a molecular biomarker signature for precise treatment.     

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.     

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.     

LAMB3在头颈鳞癌中的表达及预后价值分析

摘要 目的：研究层黏连蛋白亚基β3（LAMB3）在头颈鳞癌（HNSCC）中的表达、预后价值、免疫相关性及作用机制。方法：利用基于癌症基因组图谱（TCGA）的多个在线数据库,分析LAMB3在HNSCC肿瘤组织中的差异表达、预后价值以及对免疫细胞浸润的影响;通过基因共表达及蛋白-蛋白相互作用网络分析,探究LAMB3参与HNSCC的生物过程以及相关信号通路;预测调控LAMB3基因的上游靶miRNA。结果：LAMB3在HNSCC中高表达,LAMB3表达水平与年龄分布、肿瘤病理分级、淋巴结转移状态和HPV病毒感染状态相关;预后分析显示LAMB3在HNSCC中高表达预示患者不良预后,在肿瘤分期为Stage3、白色人种、肿瘤病理分期为Grade1和Grade2以及肿瘤新生抗原高表达HNSCC患者中总生存期更短;LAMB3表达水平影响HNSCC肿瘤微环境,LAMB3高表达能够降低B细胞、CD8⁺T 细胞浸润程度,升高CD4⁺T细胞浸润程度;基因共表达和功能相关分析表明,LAMB3高表达与基底膜形成、细胞连接、细胞迁移等生物过程有关;LAMB3潜在靶miRNA是hsa-miR-24-3p。结论：LAMB3在 HNSCC中高表达与免疫细胞浸润和患者不良预后相关,可作为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.     

Identification of hub genes associated with neutrophils infiltration in colorectal cancer

Colorectal cancer (CRC) is the leading cause of cancer-related mortality in the world. Accumulating evidence indicate that tumour infiltrating immune cells participated in cancer progression. Among them, tumour infiltrating neutrophils (TINs) are reported to play crucial role in various cancers. In this study, we used CIBERSORTx, a digital cytometry tool to evaluate the neutrophils infiltration in CRC based on gene expression data of CRC tissues from GSE39582 data set and The Cancer Genome Atlas data set (TCGA-COAD and TCGA-READ). Weighted gene co-expression network analysis (WGCNA) was conducted in GSE39582 data set to identify hub genes associated with neutrophil infiltration. The association of hub gene and neutrophils was then validated in TCGA cohorts and an independent RJ cohort. Functional analysis was performed to investigate the molecular mechanisms of the interested hub gene. We found that neutrophil infiltration is elevated in CRC tissues, and it is related to a poorer prognosis. A total of 18 gene modules are identified by WGCNA in GSE39582 data set, among which lightcyan module is significantly correlated with neutrophils infiltration. Furthermore, Superoxide Dismutase 2 (SOD2) in lightcyan module was proved to correlated with neutrophils infiltration in various cancer types. In addition, SOD2 expression is highly associated with several chemokines, including CXCL8, a neutrophils-related attractant, and functional analysis revealed that SOD2 is involved in neutrophils recruitment biological process. These results indicate that an ‘SOD2-CXCL8-neutrophil recruitment’ axis plays a potential role in colorectal cancer progression.     

Bioinformatics analysis of the prognostic value of NEK8 and its effects on immune cell infiltration in glioma

Glioma is the most common malignancy of the nervous system with high rates of recurrence and mortality, even after surgery. The 5-year survival rate is only about 5%. NEK8 is involved in multiple biological processes in a variety of cancers; however, its role in glioma is still not clear. In the current study, we evaluated the prognostic value of NEK8, as well as its role in the pathogenesis of glioma. Using a bioinformatics approach and RNA-seq data from public databases, we found that NEK8 expression is elevated in glioma tissues; we further verified this result by RT-PCR, Western blotting and immunochemistry using clinical samples. Functional enrichment analyses of genes with correlated expression indicated that elevated NEK8 expression is associated with increased immune cell infiltration in glioma and may affect the tumour microenvironment via the regulation of DNA damage/repair. Survival analyses revealed that high levels of NEK8 are associated with a poorer prognosis; higher WHO grade, IDH status, 1p/19q codeletion, age and NEK8 were identified as an independent prognostic factor. These findings support the crucial role of NEK8 in the progression of glioma via effects on immune cell infiltration and suggest that it is a new prognostic biomarker.     

Prognostic power of a lipid metabolism gene panel for diffuse gliomas

Lipid metabolism reprogramming plays important role in cell growth, proliferation, angiogenesis and invasion in cancers. However, the diverse lipid metabolism programmes and prognostic value during glioma progression remain unclear. Here, the lipid metabolism‐related genes were profiled using RNA sequencing data from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) database. Gene ontology (GO) and gene set enrichment analysis (GSEA) found that glioblastoma (GBM) mainly exhibited enrichment of glycosphingolipid metabolic progress, whereas lower grade gliomas (LGGs) showed enrichment of phosphatidylinositol metabolic progress. According to the differential genes of lipid metabolism between LGG and GBM, we developed a nine‐gene set using Cox proportional hazards model with elastic net penalty, and the CGGA cohort was used for validation data set. Survival analysis revealed that the obtained gene set could differentiate the outcome of low‐ and high‐risk patients in both cohorts. Meanwhile, multivariate Cox regression analysis indicated that this signature was a significantly independent prognostic factor in diffuse gliomas. Gene ontology and GSEA showed that high‐risk cases were associated with phenotypes of cell division and immune response. Collectively, our findings provided a new sight on lipid metabolism in diffuse gliomas.     

Establishment and experimental validation of an immune miRNA signature for assessing prognosis and immune landscape of patients with colorectal cancer

As essential regulators of gene expression, miRNAs are engaged in the initiation and progression of colorectal cancer (CRC), including antitumour immune response. In this study, we proposed an integrated algorithm, ImmuMiRNA, for identifying miRNA modulators of immune-associated pathways. Based on these immune-associated miRNAs, we applied the LASSO algorithm to develop a reliable and individualized signature for evaluating overall survival (OS) and inflammatory landscape of CRC patients. An external public data set and qRT-PCR data from 40 samples were further utilized to validate this signature. As a result, an immune-associated miRNA prognostic signature (IAMIPS) consisting of three miRNAs (miR-194-3P, miR-216a-5p and miR-3677-3p) was established and validated. Patients in the high-risk group possessed worse OS. After stratification for clinical factors, the signature remained a powerful independent predictor for OS. IAMIPS displayed much better accuracy than the traditional clinical stage in assessing the prognosis of CRC. Further analysis revealed that patients in the high-risk group were characterized by inflammatory response, abundance immune cell infiltration, and higher immune checkpoint profiles and tumour mutation burden (TMB). In conclusion, the IAMIPS is highly predictive of OS in patients with CRC, which may serve as a powerful prognostic tool to further optimize immunotherapies for cancer.     

Construction and validation of a seven-gene signature for predicting overall survival in patients with kidney renal clear cell carcinoma via an integrated bioinformatics analysis

ABSTRACT

Kidney renal clear cell carcinoma (KIRC) remains a significant challenge worldwide because of its poor prognosis and high mortality rate, and accurate prognostic gene signatures are urgently required for individual therapy. This study aimed to construct and validate a seven-gene signature for predicting overall survival (OS) in patients with KIRC. The mRNA expression profile and clinical data of patients with KIRC were obtained from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC). Prognosis-associated genes were identified, and a prognostic gene signature was constructed. Then, the prognostic efficiency of the gene signature was assessed. The results obtained using data from the TCGA were validated using those from the ICGC and other online databases. Gene set enrichment analyses (GSEA) were performed to explore potential molecular mechanisms. A seven-gene signature (PODXL, SLC16A12, ZIC2, ATP2B3, KRT75, C20orf141, and CHGA) was constructed, and it was found to be effective in classifying KIRC patients into high- and low-risk groups, with significantly different survival based on the TCGA and ICGC validation data set. Cox regression analysis revealed that the seven-gene signature had an independent prognostic value. Then, we established a nomogram, including the seven-gene signature, which had a significant clinical net benefit. Interestingly, the seven-gene signature had a good performance in distinguishing KIRC from normal tissues. GSEA revealed that several oncological signatures and GO terms were enriched. This study developed a novel seven-gene signature and nomogram for predicting the OS of patients with KIRC, which may be helpful for clinicians in establishing individualized treatments.     

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.     

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.     

The landscape and prognostic value of immune characteristics in uterine corpus endometrial cancer

In the present study, we explored the clinical and immunological characteristics of 575 uterine corpus endometrial carcinoma (UCEC) samples obtained from The Cancer Genome Atlas (TCGA) using the ESTIMATE and CIBERSORT algorithms. First, Kaplan–Meier and univariate Cox regression analyses indicated that the immune cell score was a prognostic factor for overall survival (OS) and recurrence-free survival (RFS). Multivariate Cox regression analysis further revealed that the immune cell score was an independent prognostic factor for UCEC patients. Second, we investigated the correlation between the infiltration levels of 22 types of immune cells and the immune score. Survival analysis based on the 22 immune cell types showed that higher levels of regulatory T cell, activated NK cell, and follicular helper T-cell infiltration were associated with longer OS, while higher levels of CD8+ T cell and naive B-cell infiltration were associated with longer RFS. Next, we performed differential expression and prognosis analyses on 1534 immune-related genes and selected five from 14 candidate genes to construct a prognostic prediction model. The area under the receiver-operating characteristic (ROC) curve (AUC) for 3- and 5-year survival were 0.711 and 0.728, respectively. Further validation using a stage I–II subgroup showed similar results, presenting AUC values for 3- and five-year survival of 0.677 and 0.692, respectively. Taken together, the present study provides not only a deeper understanding of the relationship between UCEC and the immune landscape but also guidance for the future development of UCEC immunotherapy.     

MicroRNA signature predicts survival in papillary thyroid carcinoma

Papillary thyroid cancer (PTC) accounts for the majority of malignant thyroid tumors. Recently, several microRNA (miRNA) expression profiling studies have used bioinformatics to suggest miRNA signatures as potential prognostic biomarkers in various malignancies. However, a prognostic miRNA biomarker has not yet been established for PTC. The aim of the present study was to identify miRNAs with prognostic value for the overall survival (OS) of patients with PTC by analyzing high-throughput miRNA data and their associated clinical characteristics downloaded from The Cancer Genome Atlas database. From our dataset, 150 differentially expressed miRNAs were identified between tumor and nontumor samples; of these miRNAs, 118 were upregulated and 32 were downregulated. Among the 150 differentially expressed miRNAs, a four miRNA signature was identified that reliably predicts OS in patients with PTC. This miRNA signature was able to classify patients into a high-risk group and a low-risk group with a significant difference in OS (P < .01). The prognostic value of the signature was validated in a testing set ( P < .01). The four miRNA signature was an independent prognostic predictor according to the multivariate analysis and demonstrated good performance in predicting 5-year disease survival with an area under the receiver operating characteristic curve area under the curve (AUC) score of 0.886. Thus, this signature may serve as a novel biomarker for predicting the survival of patients with PTC.     

Comprehensive pan-cancer analysis identifies centromere-associated protein E as a novel prognostic and immunological biomarker in human tumors

Centromere-associated protein E (CENP-E), a core component of the kinetochore, mediates chromosome congression and spindle microtubule capture during mitosis. Partial experimental evidence has illustrated the carcinogenic effects of CENPE in tumors, but the corresponding pan-cancer analysis of CENPE still lacking. Based on public databases, including the Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), and Human Protein Atlas (HPA), we take an array of bioinformatics methods to investigate the potential oncogenic roles of CENPE. Then, we validated CENPE, cell cycle-related proteins, and immune checkpoint molecule findings expression in clinical colon cancer samples by western blot. Our results showed that CENPE was up-regulated in almost all tumors, and the expression level of CENPE was associated with worse overall survival (OS) and disease-specific survival (DSS) in patients. The strong relationship between CENPE with gene mutation and MMR has also been validated. Moreover, CENPE gene expression was positively correlated with immune checkpoint molecular, and reversely correlated with infiltration levels of most immune cells. In the human colon cancer tissues, the expression of CENPE, cell cycle-related proteins, and immune checkpoint molecules were significantly higher than in the adjacent normal tissues. Our results indicated that CENPE can function as an oncogene in various cancers, and may be regarded as a promising prognostic and diagnostic biomarker in cancer treatment.     

A six-gene-based signature for breast cancer radiotherapy sensitivity estimation

Breast cancer (BRCA) represents the most common malignancy among women worldwide with high mortality. Radiotherapy is a prevalent therapeutic for BRCA that with heterogeneous effectiveness among patients. Here, we proposed to develop a gene expression-based signature for BRCA radiotherapy sensitivity estimation. Gene expression profiles of BRCA samples from the Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) were obtained and used as training and independent testing dataset, respectively. Differential expression genes (DEGs) in BRCA samples compared with their paracancerous samples in the training set were identified by using the edgeR Bioconductor package. Univariate Cox regression analysis and LASSO Cox regression method were applied to screen optimal genes for constructing a radiotherapy sensitivity estimation signature. Nomogram combining independent prognostic factors was used to predict 1-, 3-, and 5-year OS of radiation-treated BRCA patients. Relative proportions of tumor infiltrating immune cells (TIICs) calculated by CIBERSORT and mRNA levels of key immune checkpoint receptors was adopted to explore the relation between the signature and tumor immune response. As a result, 603 DEGs were obtained in BRCA tumor samples, six of which were retained and used to construct the radiotherapy sensitivity prediction model. The signature was proved to be robust in both training and testing sets. In addition, the signature was closely related to the immune microenvironment of BRCA in the context of TIICs and immune checkpoint receptors’ mRNA levels. In conclusion, the present study obtained a radiotherapy sensitivity estimation signature for BRCA, which should shed new light in clinical and experimental research.     

A robust immune-related lncRNA signature for the prognosis of human colorectal cancer

Background: Colorectal cancer (CRC) is one of the most prevalent malignant cancers worldwide. Immune-related long non-coding RNAs (IRlncRNAs) are proved to be essential in the development and progression of carcinoma. The purpose of the present study was to develop and validate a prognostic IRlncRNA signature for CRC patients.Methods: Gene expression profiles of CRC samples were downloaded from The Cancer Genome Atlas (TCGA) database. Immune-related genes were obtained from the ImmPort database and were used to identify IRlncRNA by correlation analysis. Through LASSO Cox regression analyses, a prognostic signature was constructed. Functional enrichment analysis was performed by gene set enrichment analysis (GSEA). TIMER2.0 web server and tumor immune dysfunction and exclusion (TIDE) algorithm were employed to analyze the association between our model and tumor-infiltrating immune cells and immunotherapy response. The expression levels of IRlncRNAs in cell lines were detected by quantitative real-time PCR (qPCR).Results: A 9-IRlncRNA signature was developed by a LASSO Cox proportional regression model. Based on the signature, CRC patients were divided into high- and low-risk groups with different prognoses. GSEA results indicated that patients in high-risk group were associated with cancer-related pathways. In addition, patients in low-risk group were found to have more infiltration of anti-tumor immune cells and might show a favorable response to immunotherapy. Finally, the result of qPCR revealed that most IRlncRNAs were differently expressed between normal and tumor cell lines.Conclusion: The constructed 9-IRlncRNA signature has potential to predict the prognosis of CRC patients and may be helpful to guide personalized immunotherapy.