Prognostic value of metabolic genes in lung adenocarcinoma via integrative analyses

BackgroundLung adenocarcinoma (LUAD) is the most common malignant lung tumor. Metabolic pathway reprogramming is an important hallmark of physiologic changes in cancers. However, the mechanisms through which these metabolic genes and pathways function in LUAD as well as their prognostic values have not been fully established.MethodsFour publicly available datasets from GEO and TCGA were used to identify differentially expressed genes (DEGs) in LUAD, which were then subjected to GO and KEGG pathway enrichment analysis. Associations between metabolic gene expressions with overall survival, tumor stage, TP53 mutation status, and infiltrated immune cells were investigated. Protein-protein interactions were evaluated using GeneMANIA and Metascape.ResultsBy integrating four public datasets, 247 DEGs were identified in LUAD. These DEGs were significantly enriched in regulation of chromosome segregation, centromeric region, and histone kinase activity GO terms, as well as in cell cycle, p53 signaling pathway, metabolic pathways, and other KEGG pathways. Elevated expressions of ten metabolic genes in LUAD were significantly associated with poor survival outcomes. These metabolic genes were highly expressed in more advanced tumor stage and TP53 mutated patients. Moreover, expression levels were significantly correlated with tumor-infiltrating immune cells. PPI interaction analysis revealed that the top 20 genes interacting with each metabolic gene were significantly enriched in DNA replication, response to radiation, and central carbon metabolism in cancer.ConclusionThis study elucidates on molecular changes in metabolic genes in LUAD, which may inform the development of genetically oriented diagnostic approaches and effective treatment options.     

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.     

Comprehensive molecular analysis of a four-pyroptosis-gene signature with prognosis and immune landscape in lung adenocarcinoma

Pyroptosis plays an important role in tumor immunity. However, the biological behavior and prognostic significance of pyroptosis remain unclear. We identified 41 pyroptosis regulators differently expressed in lung adenocarcinoma (LUAD). All cases of LUAD can be classified into two molecular subtypes using unsupervised clustering algorithm. Using multiple analyses, a four-pyroptosis-gene signature was constructed, and all LUAD patients were categorized as low-risk or high-risk with a longer overall survival (OS) time in the low-risk group(P < 0.001). This signature had power prognosis and stratification which was validated by six independent datasets and clinical subtypes. Besides, this signature showed distinct clinical outcomes, immune landscapes in different risk groups. Moreover, the low-risk group had a higher response against immunotherapy with a lower TIDE score. Importantly, this signature surpassed other biomarkers (TIDE, TMB, PD-L1) in predicting prognosis. Overall, the current study might help with precise prognostic prediction and crucial treatment strategies, eventually promoting tailored therapy for LUAD patients.     

Systematic profiling of invasion-related gene signature predicts prognostic features of lung adenocarcinoma

Due to the high heterogeneity of lung adenocarcinoma (LUAD), molecular subtype based on gene expression profiles is of great significance for diagnosis and prognosis prediction in patients with LUAD. Invasion-related genes were obtained from the CancerSEA database, and LUAD expression profiles were downloaded from The Cancer Genome Atlas. The ConsensusClusterPlus was used to obtain molecular subtypes based on invasion-related genes. The limma software package was used to identify differentially expressed genes (DEGs). A multi-gene risk model was constructed by Lasso-Cox analysis. A nomogram was also constructed based on risk scores and meaningful clinical features. 3 subtypes (C1, C2 and C3) based on the expression of 97 invasion-related genes were obtained. C3 had the worst prognosis. A total of 669 DEGs were identified among the subtypes. Pathway enrichment analysis results showed that the DEGs were mainly enriched in the cell cycle, DNA replication, the p53 signalling pathway and other tumour-related pathways. A 5-gene signature (KRT6A, MELTF, IRX5, MS4A1 and CRTAC1) was identified by using Lasso-Cox analysis. The training, validation and external independent cohorts proved that the model was robust and had better prediction ability than other lung cancer models. The gene expression results showed that the expression levels of MS4A1 and KRT6A in tumour tissues were higher than in normal tissues, while CRTAC1 expression in tumour tissues was lower than in normal tissues. The 5-gene signature prognostic stratification system based on invasion-related genes could be used to assess prognostic risk in patients with LUAD.     

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.     

肺腺癌诊断标志物筛选及免疫细胞浸润分析

用生物信息学方法筛选肺腺癌(Lung adenocarcinoma,LUAD)的诊断生物标志物,并分析肺腺癌中免疫细胞浸润情况。从GEO和TCGA数据库下载肺腺癌的表达数据集,利用R软件筛选肺腺癌与正常肺组织间的差异表达基因(DEGs),使用DAVID网站对DEGs进行GO及KEGG富集分析,使用STRING及Cytoscape等工具对DEGs构建蛋白相互作用网络并筛选hub基因;利用Kaplan-Meier法对DEGs进行生存分析,并对hub基因进行ROC分析筛选诊断生物标志物,利用GSEA预测有预后价值的基因参与的信号通路;并用Cibersort软件反卷积算法分析肺腺癌中免疫细胞浸润情况。共得到肺腺癌的234个DEGs,这些基因主要参与信号转导、物质代谢、免疫反应等相关信号通路;构建PPI网络筛选出的20个hub基因中8个存在预后价值(CCNA2、DLGAP5、HMMR、MMP1、MMP9、MMP13、SPP1、TOP2A),ROC分析中DLGAP5、SPP1值分别是0.703、0.706;DLGAP5、SPP1基因表达水平与肺腺癌组织浆细胞、未活化的CD4⁺记忆细胞、调节T细胞、巨噬细胞M0、M1、M2及中性粒细胞浸润密切相关(P＜0.05)。肺腺癌中DLGAP5、SPP1具有较高诊断价值且参与肺腺癌组织免疫细胞浸润;DLGAP5、SPP1基因可作为肺腺癌诊断的生物标志物,可为肺腺癌的靶向治疗研究提供新思路。     

Construction and validation of a novel apoptosis-associated prognostic signature related to osteosarcoma metastasis and immune infiltration

BackgroundApoptosis played vital roles in the formation and progression of osteosarcoma. However, no studies elucidated the prognostic relationships between apoptosis-associated genes (AAGs) and osteosarcoma.MethodsThe differentially expressed genes associated with osteosarcoma metastasis and apoptosis were identified from GEO and MSigDB databases. The apoptosis-associated prognostic signature was established through univariate and multivariate cox regression analyses. The Kaplan–Meier (KM) survival curve, ROC curve and nomogram were constructed to investigate the predictive value of this signature. CIBERSORT algorithm and ssGSEA were used to explore the relationships between immune infiltration and AAG signature. The above results were validated in another GEO dataset and the expression of AAGs was also validated in osteosarcoma patient samples by immunohistochemistry.ResultsHSPB1 and IER3 were involved in AAG signature. In training and validation datasets, apoptosis-associated risk scores were negatively related to patient survival rates and the AAG signature was regarded as the independent prognostic factor. ROC and calibration curves demonstrated the signature and nomogram were reliable. GSEA revealed the signature related to immune-associated pathways. ssGSEA indicated that one immune cell and three immune functions were significantly dysregulated. The immunohistochemistry analyses of patients’ samples revealed that AAGs were significantly differently expressed between metastasis and non-metastasis osteosarcomas.ConclusionsThe present study identified and validated a novel apoptosis-associated prognostic signature related to osteosarcoma metastasis. It could serve as the potential biomarker and therapeutic targets for osteosarcoma in the future.     

Discovery of core gene families associated with liver metastasis in colorectal cancer and regulatory roles in tumor cell immune infiltration

In this study, we aimed to uncover genes that drive the pathogenesis of liver metastasis in colorectal cancer (CRC), and identify effective genes that could serve as potential therapeutic targets for treating with colorectal liver metastasis patients based on two GEO datasets. Several bioinformatics approaches were implemented. First, differential expression analysis screened out key differentially expressed genes (DEGs) across the two GEO datasets. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we identified the enrichment functions and pathways of the DEGs that were associated with liver metastasis in CRC. Second, immune infiltration analysis identified key immune signature gene sets associated with CRC liver metastasis, among which two key immune gene families (CD and CCL) identified as key DEGs were filtered by protein-protein interaction (PPI) network. Some of the members in these gene families were associated with disease free survival (DFS) or overall survival (OS) in two subtypes of CRC, namely COAD and READ. Finally, functional enrichment analysis of the two gene families and their neighboring genes revealed that they were closely associated with cytokine, leukocyte proliferation and chemotaxis. These results are valuable in comprehending the pathogenesis of liver metastasis in CRC, and are of seminal importance in understanding the role of immune tumor infiltration in CRC. Our study also identified potentially effective therapeutic targets for liver metastasis in CRC including CCL20, CCL24 and CD70.     

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.     

Typical tumor immune microenvironment status determine prognosis in lung adenocarcinoma

BackgroundImmune cells, vital components of tumor microenvironment, regulate tumor survival and progression. Lung adenocarcinoma (LUAD), the tumor with the highest mortality rate worldwide, reconstitutes tumor immune microenvironment (TIME) to avoid immune destruction. Data have shown that TIME influences LUAD prognosis and predicts immunotherapeutic efficacy. The related information about the role of TIME's characteristics in LUAD is limited.MethodsWe performed unsupervised consensus clustering via machine-learning techniques to identify TIME clusters among 1906 patients and gathered survival data. The characteristics of TIME clusters of LUAD were visualized by multi-omics analysis, pseudo-time dynamic analysis, and enrichment analysis. TIME score model was constructed by principal component analysis. Comprehensive analysis and validation were conducted to test the prognostic efficacy and immunotherapeutic response of TIME score.ResultsTIME clusters (A, B and C) were constructed and exhibited different immune infiltration states. Multi-omics analyses included significant mutated genes (SMG), copy number variation (CNV) and cancer stemness that were significantly different among the three clusters. TIME cluster A had a lower SMG, lower CNV, and lower stemness but a higher immune infiltration level compared to TIME clusters B and C. TIME score showed that patients in low TIME score group had higher overall survival rates, higher immune infiltration level and high expression of immune checkpoints. In validation cohorts, low TIME score subgroup had better drug sensitivity and favorable immunotherapeutic response.ConclusionWe constructed a stable model of LUAD immune microenvironment characteristics that may improve the prognostic accuracy of patients, provide improved explanations of LUAD responses to immunotherapy, and provide new strategies for LUAD treatment.     

Comprehensive analysis of tumor mutation burden and immune microenvironment in gastric cancer

Tumor mutation burden (TMB) was a promising marker for immunotherapy. We aimed to investigate the prognostic role of TMB and its relationship with immune cells infiltration in gastric cancer (GC). We analyzed the mutation landscape of all GC cases and TMB of each GC patient was calculated and patients were divided into TMB-high and TMB-low group. Differentially expressed genes (DEGs) between the two groups were identified and pathway analysis was performed. The immune cells infiltration in each GC patient was evaluated and Kaplan–Meier analysis was performed to investigate the prognostic role of immune cells infiltration. At last, hub immune genes were identified and a TMB prognostic risk score (TMBPRS) was constructed to predict the survival outcome of GC patients. The relationships between mutants of hub immune genes and immune infiltration level in GC was investigated. We found higher TMB was correlated with better survival outcome and female patients, patients with T1-2 and N0 had higher TMB score. Altogether 816 DEGs were harvested and pathway analysis demonstrated that patients in TMB-high group were associated with neuroactive ligand–receptor interaction, cAMP signaling pathway, calcium signaling pathway. The infiltration of activated CD4⁺ memory T cells, follicular helper T cells, resting NK cells, M0 and M1 macrophages and neutrophils in TMB-high group were higher compared than that in TMB-low group and high macrophage infiltration was correlated with inferior survival outcome of GC patients. Lastly, the TMBPRS was constructed and GC patients with high TMBPRS had poor prognosis.     

Identification of a novel glycolysis-related gene signature that can predict the survival of patients with lung adenocarcinoma

Lung cancer is one of the most malignant cancers worldwide, and lung adenocarcinoma (LUAD) is the most common histologic subtype. Thousands of biomarkers related to the survival and prognosis of patients with this cancer type have been investigated through database mining; however, the prediction effect of a single gene biomarker is not satisfactorily specific or sensitive. Thus, the present study aimed to develop a novel gene signature of prognostic values for patients with LUAD. Using a data-mining method, we performed expression profiling of 1145 mRNAs in large cohorts with LUAD (n = 511) from The Cancer Genome Atlas database. Using the Gene Set Enrichment Analysis, we selected 198 genes related to GLYCOLYSIS, which is the most important enrichment gene set. Moreover, these genes were identified using Cox proportional regression modeling. We established a risk score staging system to predict the outcome of patients with LUAD and subsequently identified four genes (AGRN, AKR1A1, DDIT4, and HMMR) that were closely related to the prognosis of patients with LUAD. The identified genes allowed us to classify patients into the high-risk group (with poor outcome) and low-risk group (with better outcome). Compared with other clinical factors, the risk score has a better performance in predicting the outcome of patients with LUAD, particularly in the early stage of LUAD. In conclusion, we developed a four-gene signature related to glycolysis by utilizing the Cox regression model and a risk staging model for LUAD, which might prove valuable for the clinical management of patients with LUAD.     

Identification of KIF4A and its effect on the progression of lung adenocarcinoma based on the bioinformatics analysis

Background: Lung adenocarcinoma (LUAD) is the most frequent histological type of lung cancer, and its incidence has displayed an upward trend in recent years. Nevertheless, little is known regarding effective biomarkers for LUAD.Methods: The robust rank aggregation method was used to mine differentially expressed genes (DEGs) from the gene expression omnibus (GEO) datasets. The Search Tool for the Retrieval of Interacting Genes (STRING) database was used to extract hub genes from the protein–protein interaction (PPI) network. The expression of the hub genes was validated using expression profiles from TCGA and Oncomine databases and was verified by real-time quantitative PCR (qRT-PCR). The module and survival analyses of the hub genes were determined using Cytoscape and Kaplan–Meier curves. The function of KIF4A as a hub gene was investigated in LUAD cell lines.Results: The PPI analysis identified seven DEGs including BIRC5, DLGAP5, CENPF, KIF4A, TOP2A, AURKA, and CCNA2, which were significantly upregulated in Oncomine and TCGA LUAD datasets, and were verified by qRT-PCR in our clinical samples. We determined the overall and disease-free survival analysis of the seven hub genes using GEPIA. We further found that CENPF, DLGAP5, and KIF4A expressions were positively correlated with clinical stage. In LUAD cell lines, proliferation and migration were inhibited and apoptosis was promoted by knocking down KIF4A expression.Conclusion: We have identified new DEGs and functional pathways involved in LUAD. KIF4A, as a hub gene, promoted the progression of LUAD and might represent a potential therapeutic target for molecular cancer therapy.     

NCF2, MYO1F,S1PR4, and FCN1 as potential noninvasive diagnostic biomarkers in patients with obstructive coronary artery: A weighted gene co-expression network analysis

This study aims to explore the predictive noninvasive biomarker for obstructive coronary artery disease (CAD). By using the data set GSE90074, weighted gene co-expression network analysis (WGCNA), and protein–protein interactive network, construction of differentially expressed genes in peripheral blood mononuclear cells was conducted to identify the most significant gene clusters associated with obstructive CAD. Univariate and multivariate stepwise logistic regression analyses and receiver operating characteristic analysis were used to predicate the diagnostic accuracy of biomarker candidates in the detection of obstructive CAD. Furthermore, functional prediction of candidate gene biomarkers was further confirmed in ST-segment elevation myocardial infarction (STEMI) patients or stable CAD patients by using the datasets of GSE62646 and GSE59867. We found that the blue module discriminated by WGCNA contained 13 hub-genes that could be independent risk factors for obstructive CAD (P < .05). Among these 13 hub-genes, a four-gene signature including neutrophil cytosol factor 2 (NCF2, P = .025), myosin-If (MYO1F, P = .001), sphingosine-1-phosphate receptor 4 (S1PR4, P = .015), and ficolin-1 (FCN1, P = .012) alone or combined with two risk factors (male sex and hyperlipidemia) may represent potential diagnostic biomarkers in obstructive CAD. Furthermore, the messenger RNA levels of NCF2, MYO1F, S1PR4, and FCN1 were higher in STEMI patients than that in stable CAD patients, although S1PR4 showed no statistical difference (P > .05). This four-gene signature could also act as a prognostic biomarker to discriminate STEMI patients from stable CAD patients. These findings suggest a four-gene signature (NCF2, MYO1F, S1PR4, and FCN1) alone or combined with two risk factors (male sex and hyperlipidemia) as a promising prognostic biomarker in the diagnosis of STEMI. Well-designed cohort studies should be implemented to warrant the diagnostic value of these genes in clinical purpose.     

Identification of prognostic biomarkers for breast cancer brain metastases based on the bioinformatics analysis

PurposeThe prognosis of breast cancer (BC) patients who develop into brain metastases (BMs) is very poor. Thus, it is of great significance to explore the etiology of BMs in BC and identify the key genes involved in this process to improve the survival of BC patients with BMs.Patients and methodsThe gene expression data and the clinical information of BC patients were downloaded from TCGA and GEO database. Differentially expressed genes (DEGs) in TCGA-BRCA and GSE12276 were overlapped to find differentially expressed metastatic genes (DEMGs). The protein-protein interaction (PPI) network of DEMGs was constructed via STRING database. ClusterProfiler R package was applied to perform the gene ontology (GO) enrichment analysis of DEMGs. The univariate Cox regression analysis and the Kaplan-Meier (K-M) curves were plotted to screen DEMGs associated with the overall survival and the metastatic recurrence survival, which were identified as the key genes associated with the BMs in BC. The immune infiltration and the expressions of immune checkpoints for BC patients with brain relapses and BC patients with other relapses were analyzed respectively. The correlations among the expressions of key genes and the differently infiltrated immune cells or the differentially expressed immune checkpoints were calculated. The gene set enrichment analysis (GSEA) of each key gene was conducted to investigate the potential mechanisms of key genes involved in BC patients with BMs. Moreover, CTD database was used to predict the drug-gene interaction network of key genes.ResultsA total of 154 DEGs were identified in BC patients at M0 and M1 in TCGA database. A total of 667 DEGs were identified in BC patients with brain relapses and with other relapses. By overlapping these DEGs, 17 DEMGs were identified, which were enriched in the cell proliferation related biological processes and the immune related molecular functions. The univariate Cox regression analysis and the Kaplan-Meier curves revealed that CXCL9 and GPR171 were closely associated with the overall survival and the metastatic recurrence survival and were identified as key genes associated with BMs in BC. The analyses of immune infiltration and immune checkpoint expressions showed that there was a significant difference of the immune microenvironment between brain relapses and other relapses in BC. GSEA indicated that CXCL9 and GPR171 may regulate BMs in BC via the immune-related pathways.ConclusionOur study identified the key genes associated with BMs in BC patients and explore the underlying mechanisms involved in the etiology of BMs in BC. These findings may provide a promising approach for the treatments of BC patients with BMs.     

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.