scholarly journals Comprehensive transcriptomic characterization reveals core genes and module associated with immunological changes via 1619 samples of brain glioma

2021 ◽  
Vol 12 (12) ◽  
Author(s):  
Ying Zhang ◽  
Wenping Ma ◽  
Wenhua Fan ◽  
Changyuan Ren ◽  
Jianbao Xu ◽  
...  

AbstractGlioma is the most common primary malignant brain tumor with limited treatment options and poor prognosis. To investigate the potential relationships between transcriptional characteristics and clinical phenotypes, we applied weighted gene co-expression network analysis (WGCNA) to construct a free-scale gene co-expression network yielding four modules in gliomas. Turquoise and yellow modules were positively correlated with the most malignant glioma subtype (IDH-wildtype glioblastomas). Of them, genes in turquoise module were mainly involved in immune-related terms and were regulated by NFKB1, RELA, SP1, STAT1 and STAT3. Meanwhile, genes in yellow module mainly participated in cell-cycle and division processes and were regulated by E2F1, TP53, E2F4, YBX1 and E2F3. Furthermore, 14 genes in turquoise module were screened as hub genes. Among them, five prognostic hub genes (TNFRSF1B, LAIR1, TYROBP, VAMP8, and FCGR2A) were selected to construct a prognostic risk score model via LASSO method. The risk score of this immune-related gene signature is associated with clinical features, malignant phenotype, and somatic alterations. Moreover, this signature showed an accurate prediction of prognosis across different clinical and pathological subgroups in three independent datasets including 1619 samples. Our results showed that the high-risk group was characterized by active immune-related activities while the low-risk group enriched in neurophysiological-related pathway. Importantly, the high-risk score of our immune signature predicts an enrichment of glioma-associated microglia/macrophages and less response to immune checkpoint blockade (ICB) therapy in gliomas. This study not only provides new insights into the molecular pathogenesis of glioma, but may also help optimize the immunotherapies for glioma patients.

2021 ◽  
Vol 12 ◽  
Author(s):  
Susu Zheng ◽  
Xiaoying Xie ◽  
Xinkun Guo ◽  
Yanfang Wu ◽  
Guobin Chen ◽  
...  

Pyroptosis is a novel kind of cellular necrosis and shown to be involved in cancer progression. However, the diverse expression, prognosis and associations with immune status of pyroptosis-related genes in Hepatocellular carcinoma (HCC) have yet to be analyzed. Herein, the expression profiles and corresponding clinical characteristics of HCC samples were collected from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Then a pyroptosis-related gene signature was built by applying the least absolute shrinkage and selection operator (LASSO) Cox regression model from the TCGA cohort, while the GEO datasets were applied for verification. Twenty-four pyroptosis-related genes were found to be differentially expressed between HCC and normal samples. A five pyroptosis-related gene signature (GSDME, CASP8, SCAF11, NOD2, CASP6) was constructed according to LASSO Cox regression model. Patients in the low-risk group had better survival rates than those in the high-risk group. The risk score was proved to be an independent prognostic factor for overall survival (OS). The risk score correlated with immune infiltrations and immunotherapy responses. GSEA indicated that endocytosis, ubiquitin mediated proteolysis and regulation of autophagy were enriched in the high-risk group, while drug metabolism cytochrome P450 and tryptophan metabolism were enriched in the low-risk group. In conclusion, our pyroptosis-related gene signature can be used for survival prediction and may also predict the response of immunotherapy.


2021 ◽  
Vol 8 ◽  
Author(s):  
Shirui Chen ◽  
Tielong Wang ◽  
Tao Luo ◽  
Shujiao He ◽  
Changjun Huang ◽  
...  

Background: Early allograft dysfunction (EAD) is correlated with poor patient or graft survival in liver transplantation. However, the power of distinct definitions of EAD in prediction of graft survival is unclear.Methods: This retrospective, single-center study reviewed data of 677 recipients undergoing orthotopic liver transplant between July 2015 and June 2020. The following EAD definitions were compared: liver graft assessment following transplantation (L-GrAFT) risk score model, early allograft failure simplified estimation score (EASE), model for early allograft function (MEAF) scoring, and Olthoff criteria. Risk factors for L-GrAFT7 high risk group were evaluated with univariate and multivariable logistic regression analysis.Results: L-GrAFT7 had a satisfied C-statistic of 0.87 in predicting a 3-month graft survival which significantly outperformed MEAF (C-statistic = 0.78, P = 0.01) and EAD (C-statistic = 0.75, P < 0.001), respectively. L-GrAFT10, EASE was similar to L-GrAFT7, and they had no statistical significance in predicting survival. Laboratory model for end-stage liver disease score and cold ischemia time are risk factors of L-GrAFT7 high-risk group.Conclusion: L-GrAFT7 risk score is capable for better predicting the 3-month graft survival than the MEAF and EAD in a Chinese cohort, which might standardize assessment of early graft function and serve as a surrogate endpoint in clinical trial.


2021 ◽  
Vol 11 ◽  
Author(s):  
Xiaoxia Tong ◽  
Xiaofei Qu ◽  
Mengyun Wang

BackgroundCutaneous melanoma (CM) is one of the most aggressive cancers with highly metastatic ability. To make things worse, there are limited effective therapies to treat advanced CM. Our study aimed to investigate new biomarkers for CM prognosis and establish a novel risk score system in CM.MethodsGene expression data of CM from Gene Expression Omnibus (GEO) datasets were downloaded and analyzed to identify differentially expressed genes (DEGs). The overlapped DEGs were then verified for prognosis analysis by univariate and multivariate COX regression in The Cancer Genome Atlas (TCGA) datasets. Based on the gene signature of multiple survival associated DEGs, a risk score model was established, and its prognostic and predictive role was estimated through Kaplan-Meier (K-M) analysis and log-rank test. Furthermore, the correlations between prognosis related genes expression and immune infiltrates were analyzed via Tumor Immune Estimation Resource (TIMER) site.ResultsA total of 103 DEGs were obtained based on GEO cohorts, and four genes were verified in TCGA datasets. Subsequently, four genes (ADAMDEC1, GNLY, HSPA13, and TRIM29) model was developed by univariate and multivariate Cox regression analyses. The K-M plots showed that the high-risk group was associated with shortened survival than that in the low-risk group (P < 0.0001). Multivariate analysis suggested that the model was an independent prognostic factor (high-risk vs. low-risk, HR= 2.06, P < 0.001). Meanwhile, the high-risk group was prone to have larger breslow depth (P< 0.001) and ulceration (P< 0.001).ConclusionsThe four-gene risk score model functions well in predicting the prognosis and treatment response in CM and will be useful for guiding therapeutic strategies for CM patients. Additional clinical trials are needed to verify our findings.


2021 ◽  
Vol 8 ◽  
Author(s):  
Lingling Guo ◽  
Yu Jing

Background: Breast cancer is one of the most common malignancies in women worldwide. The purpose of this study was to identify the hub genes and construct prognostic signature that could predict the survival of patients with breast cancer (BC).Methods: We identified differentially expressed genes between the responder group and non-responder group based on the GEO cohort. Drug-resistance hub genes were identified by weighted gene co-expression network analysis, and a multigene risk model was constructed by univariate and multivariate Cox regression analysis based on the TCGA cohort. Immune cell infiltration and mutation characteristics were analyzed.Results: A 5-gene signature (GP6, MAK, DCTN2, TMEM156, and FKBP14) was constructed as a prognostic risk model. The 5-gene signature demonstrated favorable prediction performance in different cohorts, and it has been confirmed that the signature was an independent risk indicater. The nomogram comprising 5-gene signature showed better performance compared with other clinical features, Further, in the high-risk group, high M2 macrophage scores were related with bad prognosis, and the frequency of TP53 mutations was greater in the high-risk group than in the low-risk group. In the low-risk group, high CD8+ T cell scores were associated with a good prognosis, and the frequency of CDH1 mutations was greater in the low-risk group than that in the high-risk group. At the same time, patients in the low risk group have a good response to immunotherapy in terms of immunotherapy. The results of immunohistochemistry showed that MAK, GP6, and TEMEM156 were significantly highly expressed in tumor tissues, and DCTN2 was highly expressed in normal tissues.Conclusions: Our study may find potential new targets against breast cancer, and provide new insight into the underlying mechanisms.


2021 ◽  
Vol 11 ◽  
Author(s):  
Shenbin Xu ◽  
Zefeng Wang ◽  
Juan Ye ◽  
Shuhao Mei ◽  
Jianmin Zhang

Lower-grade glioma (LGG) is characterized by genetic and transcriptional heterogeneity, and a dismal prognosis. Iron metabolism is considered central for glioma tumorigenesis, tumor progression and tumor microenvironment, although key iron metabolism-related genes are unclear. Here we developed and validated an iron metabolism-related gene signature LGG prognosis. RNA-sequence and clinicopathological data from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA) were downloaded. Prognostic iron metabolism-related genes were screened and used to construct a risk-score model via differential gene expression analysis, univariate Cox analysis, and the Least Absolute Shrinkage and Selection Operator (LASSO)-regression algorithm. All LGG patients were stratified into high- and low-risk groups, based on the risk score. The prognostic significance of the risk-score model in the TCGA and CGGA cohorts was evaluated with Kaplan-Meier (KM) survival and receiver operating characteristic (ROC) curve analysis. Risk- score distributions in subgroups were stratified by age, gender, the World Health Organization (WHO) grade, isocitrate dehydrogenase 1 (IDH1) mutation status, the O6‐methylguanine‐DNA methyl‐transferase (MGMT) promoter-methylation status, and the 1p/19q co-deletion status. Furthermore, a nomogram model with a risk score was developed, and its predictive performance was validated with the TCGA and CGGA cohorts. Additionally, the gene set enrichment analysis (GSEA) identified signaling pathways and pathological processes enriched in the high-risk group. Finally, immune infiltration and immune checkpoint analysis were utilized to investigate the tumor microenvironment characteristics related to the risk score. We identified a prognostic 15-gene iron metabolism-related signature and constructed a risk-score model. High risk scores were associated with an age of > 40, wild-type IDH1, a WHO grade of III, an unmethylated MGMT promoter, and 1p/19q non-codeletion. ROC analysis indicated that the risk-score model accurately predicted 1-, 3-, and 5-year overall survival rates of LGG patients in the both TCGA and CGGA cohorts. KM analysis showed that the high-risk group had a much lower overall survival than the low-risk group (P < 0.0001). The nomogram model showed a strong ability to predict the overall survival of LGG patients in the TCGA and CGGA cohorts. GSEA analysis indicated that inflammatory responses, tumor-associated pathways, and pathological processes were enriched in high-risk group. Moreover, a high risk score correlated with the infiltration immune cells (dendritic cells, macrophages, CD4+ T cells, and B cells) and expression of immune checkpoint (PD1, PDL1, TIM3, and CD48). Our prognostic model was based on iron metabolism-related genes in LGG, can potentially aid in LGG prognosis, and provides potential targets against gliomas.


2020 ◽  
Author(s):  
Zhenning Zou ◽  
Shuguang Liu ◽  
Yanping Ha ◽  
Bowan Huang

Abstract Increasing evidence indicated that aberrant expression of long noncoding RNAs (lncRNAs) are involved in tumorigenesis of nasopharyngeal carcinoma (NPC). The purpose of this study was to construct a lncRNA-mediated ceRNA network based on weighted correlation network analysis (WGCNA). First, modules with highly correlated genes were identified from GSE102349 via WGCNA, and the preservation of the modules was evaluated by GSE68799. Then, the differentially expressed lncRNAs and mRNAs identified from GSE12452 which belonged to the same WGCNA modules and the differentially expressed miRNAs identified from GSE32960 were used to construct a ceRNA network. The prognostic value of the network was evaluated by survival analysis. Furthermore, a risk score model for predicting progression-free survival (PFS) of NPC patients was established via lasso penalized Cox regression, and the differences in the expression of the lncRNAs between high- and low-risk groups were investigated. Finally, 14 stable modules were identified, and a ceRNA network composed of 11 lncRNAs, 15 miRNAs and 40 mRNAs was established. The lncRNAs and mRNAs in the network belonged to the turquoise and salmon modules. Survival analysis indicated that ZNF667-AS1, LDHA, LMNB2, TPI1, UNG and hsa-miR-142-3p were significantly correlated with the prognosis of NPC. Gene set enrichment analysis indicated that the up-regulation of ZNF667-AS1 was associated with some immune-related pathways. Besides, a risk score model consisting of 12 genes was constructed and showed a good performance in predicting PFS for NPC patients. Among the 11 lncRNAs in the ceRNA network, SNHG16, SNHG17, and THAP9-AS1 were up-regulated in the high-risk group of NPC, while ZNF667-AS1 was down-regulated in the high-risk group of NPC. These results will promote our understanding of the crosstalk among lncRNAs, miRNAs, and mRNAs in the tumorigenesis and progression of NPC.


PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e12547
Author(s):  
Zihao Yan ◽  
Siwen Chu ◽  
Chen Zhu ◽  
Yunhe Han ◽  
Qingyu Liang ◽  
...  

Background Despite the rise in the use of immune checkpoint blockade drugs (ICBs) in recent years, there are no ICB drugs that are currently approved or under large-scale clinical trials for glioblastoma (GBM). T-cells, which mainly mediate adaptive immunity, are an important part of the tumor immune microenvironment. The activation of T-cells in tumors plays a key role in evaluating the sensitivity of patients to immunotherapy. Therefore, we applied bioinformatics approaches to construct a T-cell activation related risk score to study the effect of the activation of T-cells on the prognosis and ICB response of patients with GBM. Materials and Methods This study collected TCGA, CGGA, and GSE16011 glioma cohorts, as well as the IMvigor210 immunotherapy dataset, with complete mRNA expression profiles and clinical information. GraphPad Prism 8 and R 3.6.3 were used for bioinformatics analysis and plotting. Results The activation of T-cells in patients with GBM is characterized by obvious heterogeneity. We established a T-cell activation-related risk score based on five univariate Cox regression prognostic genes (CD276, IL15, SLC11A1, TNFSF4, and TREML2) in GBM. The risk score was an independent risk factor for poor prognosis. The overall survival time of patients in the high-risk group was significantly lower than in the low-risk group. Moreover, the high-risk score was accompanied by a stronger immune response and a more complex tumor immune microenvironment. “Hot tumors” were mainly enriched in the high-risk group, and high-risk group patients highly expressed inhibitory immune checkpoints (PD1, PD-L1, TIM3 etc.). By combining the risk and priming scores we obtained the immunotherapy score, which was shown to be a good evaluation index for sensitivity to GBM immunotherapy. Conclusions As an independent risk factor for poor prognosis, the T-cell activation-related risk score, combined with other clinical characteristics, could efficiently evaluate the survival of patients with GBM. The immunotherapy score obtained by combining the risk and priming scores could evaluate the ICB response of patients with GBM, providing treatment opportunities.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Yang Peng ◽  
Haochen Yu ◽  
Yingzi Zhang ◽  
Fanli Qu ◽  
Zhenrong Tang ◽  
...  

AbstractFerroptosis is a new form of regulated cell death (RCD), and its emergence has provided a new approach to the progression and drug resistance of breast cancer (BRCA). However, there is still a great gap in the study of ferroptosis-related genes in BRCA, especially luminal-type BRCA patients. We downloaded the mRNA expression profiles and corresponding clinical data of BRCA patients from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA) databases. Then, we built a prognostic multigene signature with ferroptosis-related differentially expressed genes (DEGs) in the METABRIC cohort and validated it in the TCGA cohort. The predictive value of this signature was investigated in terms of the immune microenvironment and the probability of a response to immunotherapy and chemotherapy. The patients were divided into a high-risk group and a low-risk group according to the ferroptosis-associated gene signature, and the high-risk group had a worse overall survival (OS). The risk score based on the 10 ferroptosis-related gene-based signature was determined to be an independent prognostic predictor in both the METABRIC and TCGA cohorts (HR, 1.41, 95% CI, 1.14–1.76, P = 0.002; HR, 2.19, 95% CI, 1.13–4.26, P = 0.02). Gene set enrichment analysis indicated that the term “cytokine-cytokine receptor interaction” was enriched in the high-risk score subgroup. Moreover, the immune infiltration scores of most immune cells were significantly different between the two groups, the low-risk group was much more sensitive to immunotherapy, and six drugs might have potential therapeutic implications in the high-risk group. Finally, a nomogram incorporating a classifier based on the 10 ferroptosis-related genes, tumor stage, age and histologic grade was established. This nomogram showed favorable discriminative ability and could help guide clinical decision-making for luminal-type breast carcinoma.


2021 ◽  
Vol 2021 ◽  
pp. 1-13
Author(s):  
Zemin Zhu ◽  
Caixi Tang ◽  
Tao Xu ◽  
Zhijian Zhao

Background. Tumor necrosis factor (TNF) family members play a vital role in anticancer therapy. This study aimed to screen the critical markers for the prognostic analysis of pancreatic adenocarcinoma (PAAD) by analyzing the clustering patterns of TNF family members in PAAD. Methods. In this study, the NMF clustering method was adopted to cluster samples from The Cancer Genome Atlas (TCGA) to acquire the clustering pattern of the TNF family in PAAD. Differential gene analysis was performed according to TNF family gene clusters. The support vector machine (SVM) method was conducted for further gene screening, and the risk score model of the screened genes was constructed by Lasso. The single sample gene set enrichment analysis (ssGSEA) method was adopted for immunoenrichment analysis and tumor immune cycle analysis. Genes associated with risk scores were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Results. We clustered PAAD into two groups based on TNF family genes. Nineteen TNF family genes were significantly associated with the clinical characteristics of PAAD patients. The risk score formula was composed of RHOD, UBE2C, KLHDC7b, MSLN, ADAM8, NME3, GNG2, and MCOLN3. GSE57495 and GSE62452 datasets verified that patients in the high-risk group had a worse prognosis than those in the low-risk group. The risk score-related genes analyzed by GO and KEGG were mainly involved in the modulation of chemical synaptic transmission and synaptic vesicle cycle pathway. There were significant differences in the expression of 15 immune cells between the high-risk group and the low-risk group. The risk score was positively correlated with HCK, interferon, MHC-I, and STAT1. The expression of genes relevant to chemokine, immunostimulator, MHC, and receptor was strongly associated with the risk score. Conclusion. The risk score model based on the TNF family can predict the prognosis and immune status of PAAD patients. Further research is needed to verify the clinical prognostic value of risk scores.


2022 ◽  
Vol 2022 ◽  
pp. 1-17
Author(s):  
Lei Zhang ◽  
Dahai Hu ◽  
Shuchen Huangfu ◽  
Jiaxin Zhou ◽  
Wei Wang ◽  
...  

The genomic variant features (mutations, deletions, structural variants, etc.) within gastric cancer impact its evolution and immunogenicity. The tumor has developed several coping strategies to respond to these changes by DNA repair and replication (DRR). However, the intrinsic relationship between the associated DRR-related genes and gastric cancer progression remained unknown. This study selected DRR-related genes with tumor mutation burden based on the TCGA (The Cancer Genome Atlas) database of gastric cancer transcriptome and mutation data. The prognosis model of seven genes (LAMA2, CREB3L3, SELP, ABCC9, CYP1B1, CDH2, and GAMT) was constructed by a univariate and LASSO regression analysis and divided into high-risk and low-risk groups with the median risk score. Survival analysis showed that overall survival (OS) was lower in the high-risk group than that in the low-risk group. Moreover, patients with gastric cancer in the high-risk group have worse survival in different subgroups, including age, gender, histological grade, and TNM stage. The nomogram that included risk scores for DRR-related genes could accurately foresee OS of patients with gastric cancer. Interestingly, the tumor mutation burden score was higher in the low-risk group than that in the high-risk group, and the risk score for DRR-related genes was negatively correlated with tumor mutation burden in gastric cancer. Next, we further combined the risk score and tumor mutation burden to evaluate the prognosis of gastric cancer patients. The low-risk cohort had a better prognosis than the high-risk cohort in the high tumor mutation burden subgroup. The number of mutation types in the high-risk group was lower than that in the low-risk group. In the immune microenvironment of gastric cancer, more naïve B cells, memory resting CD4+ T cells, Treg cells, monocytes cells, and resting mast cells were infiltrated in the high-risk group. At last, PD-L1 and IAP expressions were negatively correlated with the risk scores; patients with gastric cancer in the low-risk group showed better immunotherapy outcomes than those in the high-risk group. Overall, the DRR-related gene signature based on tumor mutation burden is a novel biomarker for prognostic and immunotherapy response in patients with gastric cancer.


Sign in / Sign up

Export Citation Format

Share Document