scholarly journals Glycolysis Define Two Prognostic Subgroups of Lung Adenocarcinoma With Different Mutation Characteristics and Immune Infiltration Signatures

Author(s):  
Chen Huo ◽  
Meng-Yu Zhang ◽  
Rui Li ◽  
Ting-Ting Liu ◽  
Jian-Ping Li ◽  
...  

Increasing studies have proved that malignant tumors are associated with energy metabolism. This study was aimed to explore biological variables that impact the prognosis of patients in the glycolysis-related subgroups of lung adenocarcinoma (LUAD). The mRNA expression profiling and mutation data in large LUAD samples were collected from the Cancer Genome Atlas (TCGA) database. Then, we identified the expression level and prognostic value of glycolysis-related genes, as well as the fractions of 22 immune cells in the tumor microenvironment. The differences between glycolysis activity, mutation, and immune infiltrates were discussed in these groups, respectively. Two hundred fifty-five glycolysis-related genes were identified from gene set enrichment analysis (GSEA), of which 43 genes had prognostic values (p < 0.05). Next, we constructed a glycolysis-related competing endogenous RNA (ceRNA) network which related to the survival of LUAD. Then, two subgroups of LUAD (clusters 1 and 2) were identified by applying unsupervised consensus clustering to 43 glycolysis-related genes. The survival analysis showed that the cluster 1 patients had a worse prognosis (p < 0.001), and upregulated differentially expressed genes (DEGs) are interestingly enriched in malignancy-related biological processes. The differences between the two subgroups are SPTA1, KEAP1, USH2A, and KRAS among top 10 mutated signatures, which may be the underlying mechanism of grouping. Combined high tumor mutational burden (TMB) with tumor subgroups preferably predicts the prognosis of LUAD patients. The CIBERSORT algorithm results revealed that low TMB samples were concerned with increased infiltration level of memory resting CD4+ T cell (p = 0.03), resting mast cells (p = 0.044), and neutrophils (p = 0.002) in cluster 1 and high TMB samples were concerned with increased infiltration level of memory B cells, plasma cells, CD4 memory-activated T cells, macrophages M1, and activated mast cells in cluster 2, while reduced infiltration of monocytes, resting dendritic cells, and resting mast cells was captured in cluster 2. In conclusion, significant different gene expression characteristics were pooled according to the two subgroups of LUAD. The combination of subgroups, TMB and tumor-infiltrating immune cell signature, might be a novel prognostic biomarker in LUAD.

2021 ◽  
Vol 11 ◽  
Author(s):  
Wencheng Zhang ◽  
Zhouyong Gao ◽  
Mingxiu Guan ◽  
Ning Liu ◽  
Fanjie Meng ◽  
...  

Anti-silencing function 1B histone chaperone (ASF1B) is known to be an important modulator of oncogenic processes, yet its role in lung adenocarcinoma (LUAD) remains to be defined. In this study, an integrated assessment of The Cancer Genome Atlas (TCGA) and genotype-tissue expression (GTEx) datasets revealed the overexpression of ASF1B in all analyzed cancer types other than LAML. Genetic, epigenetic, microsatellite instability (MSI), and tumor mutational burden (TMB) analysis showed that ASF1B was regulated by single or multiple factors. Kaplan-Meier survival curves suggested that elevated ASF1B expression was associated with better or worse survival in a cancer type-dependent manner. The CIBERSORT algorithm was used to evaluate immune microenvironment composition, and distinct correlations between ASF1B expression and immune cell infiltration were evident when comparing tumor and normal tissue samples. Gene set enrichment analysis (GSEA) indicated that ASF1B was associated with proliferation- and immunity-related pathways. Knocking down ASF1B impaired the proliferation, affected cell cycle distribution, and induced cell apoptosis in LUAD cell lines. In contrast, ASF1B overexpression had no impact on the malignant characteristics of LUAD cells. At the mechanistic level, ASF1B served as an indirect regulator of DNA Polymerase Epsilon 3, Accessory Subunit (POLE3), CDC28 protein kinase regulatory subunit 1(CKS1B), Dihydrofolate reductase (DHFR), as established through proteomic profiling and Immunoprecipitation-Mass Spectrometry (IP-MS) analyses. Overall, these data suggested that ASF1B serves as a tumor promoter and potential target for cancer therapy and provided us with clues to better understand the importance of ASF1B in many types of cancer.


2021 ◽  
Vol 10 ◽  
Author(s):  
Liang Zhao ◽  
Jiayue Zhang ◽  
Zhiyuan Liu ◽  
Yu Wang ◽  
Shurui Xuan ◽  
...  

Alternative splicing (AS) of pre-mRNA has been widely reported to be associated with the progression of malignant tumors. However, a systematic investigation into the prognostic value of AS events in glioblastoma (GBM) is urgently required. The gene expression profile and matched AS events data of GBM patients were obtained from The Cancer Genome Atlas Project (TCGA) and TCGA SpliceSeq database, respectively. 775 AS events were identified as prognostic factors using univariate Cox regression analysis. The least absolute shrinkage and selection operator (LASSO) cox model was performed to narrow down candidate AS events, and a risk score model based on several AS events were developed subsequently. The risk score-based signature was proved as an efficient predictor of overall survival and was closely related to the tumor purity and immunosuppression in GBM. Combined similarity network fusion and consensus clustering (SNF-CC) analysis revealed two distinct GBM subtypes based on the prognostic AS events, and the associations between this novel molecular classification and clinicopathological factors, immune cell infiltration, as well as immunogenic features were further explored. We also constructed a regulatory network to depict the potential mechanisms that how prognostic splicing factors (SFs) regulate splicing patterns in GBM. Finally, a nomogram incorporating AS events signature and other clinical-relevant covariates was built for clinical application. This comprehensive analysis highlights the potential implications for predicting prognosis and clinical management in GBM.


2021 ◽  
Vol 12 ◽  
Author(s):  
Yingxi Du ◽  
Yarui Ma ◽  
Qing Zhu ◽  
Tongzheng Liu ◽  
Yuchen Jiao ◽  
...  

Background: N6-methyladenosine (m6A) is related to the progression of multiple cancers. However, the underlying influences of m6A-associated genes on the tumor immune microenvironment in hepatocellular carcinoma (HCC) remain poorly understood. Therefore, we sought to construct a survival prediction model using m6A-associated genes to clarify the molecular and immune characteristics of HCC.Methods: HCC case data were downloaded from The Cancer Genome Atlas (TCGA). Then, by applying consensus clustering, we identified two distinct HCC clusters. Next, four m6A-related genes were identified to construct a prognostic model, which we validated with Gene Expression Omnibus (GEO) and International Cancer Genome Consortium (ICGC) datasets. Additionally, the molecular and immune characteristics in different subgroups were analyzed.Results: m6A RNA methylation regulators were differentially expressed between HCC and normal samples and linked with immune checkpoint expression. Using consensus clustering, we divided HCC samples into two subtypes with distinct clinical features. Cluster 2 was associated with unfavorable prognosis, higher immune checkpoint expression and immune cell infiltration levels. In addition, the immune and carcinogenic signaling pathways were enriched in cluster 2. Furthermore, we constructed a risk model using four m6A-associated genes. Patients with different risk scores had distinct survival times, expression levels of immunotherapy biomarkers, TP53 mutation rates, and sensitivities to chemotherapy and targeted therapy. Similarly, the model exhibited an identical impact on overall survival in the validation cohorts.Conclusion: The constructed m6A-based signature may be promising as a biomarker for prognostics and to distinguish immune characteristics in HCC.


2021 ◽  
Vol 12 ◽  
Author(s):  
Jixin Wang ◽  
Xiangjun Yin ◽  
Yin-Qiang Zhang ◽  
Xuming Ji

Lung adenocarcinoma (LUAD) is a major subtype of lung cancer, the prognosis of patients with which is associated with both lncRNAs and cancer immunity. In this study, we collected gene expression data of 585 LUAD patients from The Cancer Genome Atlas (TCGA) database and 605 subjects from the Gene Expression Omnibus (GEO) database. LUAD patients were divided into high and low immune-cell-infiltrated groups according to the single sample gene set enrichment analysis (ssGSEA) algorithm to identify differentially expressed genes (DEGs). Based on the 49 immune-related DE lncRNAs, a four-lncRNA prognostic signature was constructed by applying least absolute shrinkage and selection operator (LASSO) regression, univariate Cox regression, and stepwise multivariate Cox regression in sequence. Kaplan–Meier curve, ROC analysis, and the testing GEO datasets verified the effectiveness of the signature in predicting overall survival (OS). Univariate Cox regression and multivariate Cox regression suggested that the signature was an independent prognostic factor. The correlation analysis revealed that the infiltration immune cell subtypes were related to these lncRNAs.


BMC Cancer ◽  
2020 ◽  
Vol 20 (1) ◽  
Author(s):  
Shibai Yan ◽  
Juntao Fang ◽  
Yongcai Chen ◽  
Yong Xie ◽  
Siyou Zhang ◽  
...  

Abstract Background Ovarian cancer (OV) is one of the most common malignant tumors of gynecology oncology. The lack of effective early diagnosis methods and treatment strategies result in a low five-year survival rate. Also, immunotherapy plays an important auxiliary role in the treatment of advanced OV patient, so it is of great significance to find out effective immune-related tumor markers for the diagnosis and treatment of OV. Methods Based on the consensus clustering analysis of single-sample gene set enrichment analysis (ssGSEA) score transformed via The Cancer Genome Atlas (TCGA) mRNA profile, we obtained two groups with high and low levels of immune infiltration. Multiple machine learning methods were conducted to explore prognostic genes associated with immune infiltration. Simultaneously, the correlation between the expression of mark genes and immune cells components was explored. Results A prognostic classifier including 5 genes (CXCL11, S1PR4, TNFRSF17, FPR1 and DHRS95) was established and its robust efficacy for predicting overall survival was validated via 1129 OV samples. Some significant variations of copy number on gene loci were found between two risk groups and it showed that patients with fine chemosensitivity has lower risk score than patient with poor chemosensitivity (P = 0.013). The high and low-risk groups showed significantly different distribution (P < 0.001) of five immune cells (Monocytes, Macrophages M1, Macrophages M2, T cells CD4 menory and T cells CD8). Conclusion The present study identified five prognostic genes associated with immune infiltration of OV, which may provide some potential clinical implications for OV treatment.


2021 ◽  
Vol 8 ◽  
Author(s):  
Enchong Zhang ◽  
Fujisawa Shiori ◽  
Mo Zhang ◽  
Peng Wang ◽  
Jieqian He ◽  
...  

Prostate cancer (PCa) is the most common malignancy among men worldwide. However, its complex heterogeneity makes treatment challenging. In this study, we aimed to identify PCa subtypes and a gene signature associated with PCa prognosis. In particular, nine PCa-related pathways were evaluated in patients with PCa by a single-sample gene set enrichment analysis (ssGSEA) and an unsupervised clustering analysis (i.e., consensus clustering). We identified three subtypes with differences in prognosis (Risk_H, Risk_M, and Risk_L). Differences in the proliferation status, frequencies of known subtypes, tumor purity, immune cell composition, and genomic and transcriptomic profiles among the three subtypes were explored based on The Cancer Genome Atlas database. Our results clearly revealed that the Risk_H subtype was associated with the worst prognosis. By a weighted correlation network analysis of genes related to the Risk_H subtype and least absolute shrinkage and selection operator, we developed a 12-gene risk-predicting model. We further validated its accuracy using three public datasets. Effective drugs for high-risk patients identified using the model were predicted. The novel PCa subtypes and prognostic model developed in this study may improve clinical decision-making.


Author(s):  
Xinyu Gu ◽  
Haibo Zhou ◽  
Qingfei Chu ◽  
Qiuxian Zheng ◽  
Jing Wang ◽  
...  

Background: 5-Methylcytosine (m5C) plays essential roles in hepatocellular carcinoma (HCC), but the association between m5C regulation and immune cell infiltration in HCC has not yet been clarified.Methods: In this study, we analysed 371 patients with HCC from The Cancer Genome Atlas (TCGA) database, and the expression of 13 m5C regulators was investigated. Additionally, gene set variation analysis (GSVA), unsupervised clustering analysis, single-sample gene set enrichment analysis (ssGSEA), correlation analysis, and immunohistochemical (IHC) staining were performed.Results: Among the 371 patients, 41 had mutations in m5C regulators, the frequency of which was 11.26%. Compared with normal hepatic tissues, the expression of m5C regulators with copy number variations (CNVs) expansion was significantly higher than that in HCC tissues. Then, we identified three m5C modification patterns that had obvious tumour microenvironment (TME) cell infiltration characteristics. The prognostic analysis of the three major m5C modification subtypes showed that Cluster-2 had a clear survival advantage over the others. In addition, we found that DNMT1 was highly expressed in tumour tissues compared with normal tissues in a tissue microarray (TMA) and that it was positively correlated with many TME-infiltrating immune cells. High expression of the m5C regulator DNMT1 was related to a poor prognosis in patients with HCC. Furthermore, we developed three distinct Immu-clusters. Importantly, mRNAs related to the transcription of growth factor β (TGF-β)/EMT pathway were significantly up-regulated in Immu-cluster 2, indicating that this cluster is considered to be the immune rejection phenotype. Immu-cluster 3 showed elevated expression of mRNAs related to immune checkpoint genes.Conclusion: Our work revealed the association between m5C modification and immune regulators in the TME. These findings also suggest that DNMT1 has great potential as a prognostic biomarker and therapeutic target for HCC.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Yangming Hou ◽  
Yingjuan Xu ◽  
Dequan Wu

AbstractThe infiltration degree of immune and stromal cells has been shown clinically significant in tumor microenvironment (TME). However, the utility of stromal and immune components in Gastric cancer (GC) has not been investigated in detail. In the present study, ESTIMATE and CIBERSORT algorithms were applied to calculate the immune/stromal scores and the proportion of tumor-infiltrating immune cell (TIC) in GC cohort, including 415 cases from The Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) were screened by Cox proportional hazard regression analysis and protein–protein interaction (PPI) network construction. Then ADAMTS12 was regarded as one of the most predictive factors. Further analysis showed that ADAMTS12 expression was significantly higher in tumor samples and correlated with poor prognosis. Gene Set Enrichment Analysis (GSEA) indicated that in high ADAMTS12 expression group gene sets were mainly enriched in cancer and immune-related activities. In the low ADAMTS12 expression group, the genes were enriched in the oxidative phosphorylation pathway. CIBERSORT analysis for the proportion of TICs revealed that ADAMTS12 expression was positively correlated with Macrophages M0/M1/M2 and negatively correlated with T cells follicular helper. Therefore, ADAMTS12 might be a tumor promoter and responsible for TME status and tumor energy metabolic conversion.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Wenjie Chen ◽  
Wen Li ◽  
Zhenkun Liu ◽  
Guangzhi Ma ◽  
Yunfu Deng ◽  
...  

AbstractTo identify the prognostic biomarker of the competitive endogenous RNA (ceRNA) and explore the tumor infiltrating immune cells (TIICs) which might be the potential prognostic factors in lung adenocarcinoma. In addition, we also try to explain the crosstalk between the ceRNA and TIICs to explore the molecular mechanisms involved in lung adenocarcinoma. The transcriptome data of lung adenocarcinoma were obtained from The Cancer Genome Atlas (TCGA) database, and the hypergeometric correlation of the differently expressed miRNA-lncRNA and miRNA-mRNA were analyzed based on the starBase. In addition, the Kaplan–Meier survival and Cox regression model analysis were used to identify the prognostic ceRNA network and TIICs. Correlation analysis was performed to analysis the correlation between the ceRNA network and TIICs. In the differently expressed RNAs between tumor and normal tissue, a total of 190 miRNAs, 224 lncRNAs and 3024 mRNAs were detected, and the constructed ceRNA network contained 5 lncRNAs, 92 mRNAs and 10 miRNAs. Then, six prognostic RNAs (FKBP3, GPI, LOXL2, IL22RA1, GPR37, and has-miR-148a-3p) were viewed as the key members for constructing the prognostic prediction model in the ceRNA network, and three kinds of TIICs (Monocytes, Macrophages M1, activated mast cells) were identified to be significantly related with the prognosis in lung adenocarcinoma. Correlation analysis suggested that the FKBP3 was associated with Monocytes and Macrophages M1, and the GPI was obviously related with Monocytes and Macrophages M1. Besides, the LOXL2 was associated with Monocytes and Activated mast cells, and the IL22RA1 was significantly associated with Monocytes and Macrophages M1, while the GPR37 and Macrophages M1 was closely related. The constructed ceRNA network and identified Monocytes, Macrophages M1 and activated Mast cells are all prognostic factors for lung adenocarcinoma. Moreover, the crosstalk between the ceRNA network and TIICs might be a potential molecular mechanism involved.


2021 ◽  
Vol 20 ◽  
pp. 153303382199208
Author(s):  
Shufang Wang ◽  
Xinlong Huo

Background: Estrogen-related receptor alpha (ESRRA) was reported to play an important role in multiple biological processes of neoplastic diseases. The roles of ESRRA in endometrial cancer have not been fully investigated yet. Methods: Expression data and clinicopathological data of patients with uteri corpus endometrial carcinoma (UCEC) were obtained from The Cancer Genome Atlas (TCGA). Comprehensive bioinformatics analysis was performed, including receiver operating characteristics (ROC) curve analysis, Kaplan-Meier survival analysis, gene ontology (GO) enrichment analysis, and Gene Set Enrichment Analysis (GSEA). Immunohistochemistry was used to detect the protein expression level of ESRRA and CCK-8 assay was performed to evaluate the effect of ESRRA on the proliferation ability. Results: A total of 552 UCEC tissues and 35 normal tissues were obtained from the TCGA database. The mRNA and protein expression level of ESRRA was highly elevated in UCEC compared with normal tissues, and was closely associated with poor prognosis. ROC analysis indicated a very high diagnostic value of ESRRA for patients with UCEC. GO and GSEA functional analysis showed that ESRRA might be mainly involved in cellular metabolism processes, in turn, tumorigenesis and progression of UCEC. Knockdown of ESRRA inhibited the proliferation of UCEC cells in vitro. Further immune cell infiltration demonstrated that ESRRA enhanced the infiltration level of neutrophil cell and reduced that of T cell (CD4+ naïve), NK cell, and cancer associated fibroblast (CAF). The alteration of immune microenvironment will greatly help in developing immune checkpoint therapy for UCEC. Conclusions: Our study comprehensively analyzed the expression level, clinical value, and possible mechanisms of action of ESRRA in UCEC. These findings showed that ESRRA might be a potential diagnostic and therapeutic target.


Sign in / Sign up

Export Citation Format

Share Document