Dataset Information

Prognosis and Immunological Characteristics of PGK1 in Lung Adenocarcinoma: A Systematic Analysis.

ABSTRACT: Background: Aerobic glycolysis plays a key role in tumor metabolic reprogramming to reshape the immune microenvironment. The phosphoglycerate kinase 1 (PGK1) gene codes a glycolytic enzyme that converts 1,3-diphosphoglycerate to 3-phosphoglycerate. However, in lung adenocarcinoma (LUAD), the role of PGK1 in altering the tumor microenvironment (TME) has not yet been determined. Methods: Raw data, including bulk DNA and mRNA-seq data, methylation modification data, single-cell RNA-seq data, proteomics data, clinical case characteristics survival, immunotherapy data, and so on, were obtained from multiple independent public data sets. These data were reanalyzed to uncover the prognosis and immunological characteristics of PGK1 in LUAD. Results: We found that PGK1 mRNA and protein were considerably over-expressed in LUAD compared to normal tissue and that high PGK1 expression is associated with poorer prognostic outcomes in LUAD. The enrichment analysis of PGK1 co-expressed genes in lung adenocarcinoma revealed that PGK1 may be involved in hypoxia, metabolism, DNA synthesis, cell cycle, PI3K/AKT, and various immune and inflammatory signaling pathways. Furthermore, PGK1 is also linked to the recruitment of numerous immune cells, including aDC (dendritic cells), macrophages, and neutrophils. More importantly, PGK1 was highly expressed in immunosuppressive cells, including M2 macrophages, Tregs, and exhausted T cells, among others. Finally, higher PGK1 expression indicated significant correlations to immune checkpoints, TMB (tumor mutation burden), and high response to immunotherapy. Conclusions: The presented findings imply that PGK1, as a glycolysis core gene, may be important for the modification of the immune microenvironment by interacting with the tumor metabolism. The results of this study provide clues for a potential immunometabolic combination therapy strategy in LUAD, for which more experimental and clinical translational research is needed.

SUBMITTER: Yang Y

PROVIDER: S-EPMC9653683 | biostudies-literature | 2022 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Prognosis and Immunological Characteristics of PGK1 in Lung Adenocarcinoma: A Systematic Analysis.

Yang Yuechao Y Cui Huanhuan H Li Deheng D Gao Yang Y Chen Lei L Zhou Changshuai C Feng Mingtao M Tu Wenjing W Li Sen S Chen Xin X Hao Bin B Li Liangdong L Cao Yiqun Y

Cancers 20221025 21

<b>Background</b>: Aerobic glycolysis plays a key role in tumor metabolic reprogramming to reshape the immune microenvironment. The phosphoglycerate kinase 1 (PGK1) gene codes a glycolytic enzyme that converts 1,3-diphosphoglycerate to 3-phosphoglycerate. However, in lung adenocarcinoma (LUAD), the role of PGK1 in altering the tumor microenvironment (TME) has not yet been determined. <b>Methods</b>: Raw data, including bulk DNA and mRNA-seq data, methylation modification data, single-cell RNA-se ...[more]

PMID: 36358653

Similar Datasets

Project description:Lung adenocarcinoma (LUAD) usually contains heterogeneous histological subtypes, among which the micropapillary (MIP) subtype was associated with poor prognosis while the lepidic (LEP) subtype possessed the most favorable outcome. However, the genomic features of the MIP subtype responsible for its malignant behaviors are substantially unknown. In this study, eight FFPE samples from LUAD patients were micro-dissected to isolate MIP and LEP components, then sequenced by whole-exome sequencing. More comprehensive analyses involving our samples and public validation cohorts on the two subtypes were performed to better decipher the key biological and evolutionary mechanisms. As expected, the LEP and MIP subtypes exhibited the largest disease-free survival (DFS) differences in our patients. EGFR was found with the highest mutation frequency. Additionally, shared mutations were observed between paired LEP and MIP components from single patients, and recurrent mutations were verified in the Lung-Broad, Lung-OncoSG, and TCGA-LUAD cohorts. Distinct biological processes or pathways were involved in the evolution of the two components. Besides, analyses of copy number variation (CNV) and intratumor heterogeneity (ITH) further discovered the possible immunosurveillance escape, the discrepancy between mutation and CNV level, ITH, and the pervasive DNA damage response and WNT pathway gene alternations in the MIP component. Phylogenetic analysis of five pairs of LEP and MIP components further confirmed the presence of ancestral EGFR mutations. Through comprehensive analyses combining our samples and public cohorts, PTP4A3, NAPRT, and RECQL4 were identified to be co-amplified. Multi-omics data also demonstrated the immunosuppression prevalence in the MIP component. Our results uncovered the evolutionary pattern of the concomitant LEP and MIP components from the same patient that they were derived from the same initiation cells and the pathway-specific mutations acquired after EGFR clonal mutation could shape the subtype-specificity. We also confirmed the immunosuppression prevalence in the MIP subtype by multi-omics data analyses, which may have resulted in its unfavorable prognosis.

Project description:Lung adenocarcinoma (LUAD) is the most common type of lung cancer, accounting for approximately 85% of pulmonary malignancies. Emerging evidence has demonstrated that ferroptosis plays a central role in both immunities as well as tumor proliferation. However, the clinical significance, immunological function, and upstream modulatory mechanism of ferroptosis-related genes in LUAD remain unclear. Here, we utilized various bioinformatics data to identify differentially expressed (DEGs) and prognosis-related ferroptosis (FRGs) genes in LUAD. Based upon identified DEGs, FRG, and ceRNA modulatory networks were constructed. Pearson's correlation analysis was used to evaluate the correlation between FRGs and the tumor mutational burden, microsatellite instability, tumor-infiltrating immunity, cellular checkpoint control, and drug sensitivity in LUAD. A loss-of-function analysis was performed to verify the function of CISD1 in LUAD progression. Our findings revealed that certain FRGs (CISD1, ATP5MC3, PGD, SLC7A11, ACSL3, and FANCD2) are significantly upregulated in LUAD and that their elevated expression is associated with both advanced tumor stage and unfavorable prognosis. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results revealed these FRGs to be primarily involved in ferroptosis and glutathione metabolism in LUAD. We constructed a prognostic FRG-based model capable of accurately predicting LUAD patient overall survival with high specificity. The upstream lncRNA GSEC/miRNA-101-3p regulatory axis involving CISD1, ATP5MC3, and PGD was identified to be relevant in tumor progression. We also found GSEC, CISD1, ATP5MC3, and PGD to be upregulated, with miRNA-101-3p downregulated, in the setting of LUAD. Immunohistochemical analysis revealed CISD1, ATP5MC3, and PGD overexpression in LUAD tissue samples; CISD1 knockdown was noted to significantly inhibit LUAD proliferation and migration. In summary, this study characterizes relevant functional roles of the lncRNA GSEC/miR-101-3p axis in the setting of LUAD and suggests diagnostic and therapeutic biomarkers potentially useful in the clinical management of this illness.

Project description:BackgroundsColorectal adenocarcinoma (COAD), accounting for the most common subtype of colorectal cancer (CRC), is a kind of malignant digestive tumor. Some cell cycle checkpoints (CCCs) have been found to contribute to CRC progression, whereas the functional roles of a lot of CCCs, especially the integrated role of checkpoint mechanism in the cell cycle, remain unclear.Materials and methodsThe Genomic Data Commons (GDC) The Cancer Genome Atlas (TCGA) COAD cohort was retrieved as the training dataset, and GSE24551 and GSE29623 were downloaded from Gene Expression Omnibus (GEO) as the validation datasets. A total of 209 CCC-related genes were derived from the Gene Ontology Consortium and were subsequently enrolled in the univariate, multivariate, and least absolute shrinkage and selection operator (LASSO) Cox regression analyses, finally defining a CCC signature. Cell proliferation and Transwell assay analyses were utilized to evaluate the functional roles of signature-related CCCs. The underlying CCC signature, molecular characteristics, immune-related features, and therapeutic response were finally estimated. The Genomics of Drug Sensitivity in Cancer (GDSC) database was employed for the evaluation of chemotherapeutic responses.ResultsThe aberrant gene expression of CCCs greatly contributed to COAD development and progression. Univariate Cox regression analysis identified 27 CCC-related genes significantly affecting the overall survival (OS) of COAD patients; subsequently, LASSO analysis determined a novel CCC signature. Noticeably, CDK5RAP2, MAD1L1, NBN, RGCC, and ZNF207 were first identified to be correlated with the prognosis of COAD, and it was proven that all of them were significantly correlated with the proliferation and invasion of HCT116 and SW480 cells. In TCGA COAD cohort, CCC signature robustly stratified COAD patients into high and low CCC score groups (median OS: 57.24 months vs. unreached, p< 0.0001), simultaneously, with the good AUC values for OS prediction at 1, 2, and 3 years were 0.74, 0.78, and 0.77. Furthermore, the prognostic capacity of the CCC signature was verified in the GSE24551 and GSE29623 datasets, and the CCC signature was independent of clinical features. Moreover, a higher CCC score always indicated worse OS, regardless of clinical features, histological subtypes, or molecular subgroups. Intriguingly, functional enrichment analysis confirmed the CCC score was markedly associated with extracellular, matrix and immune (chemokine)-related signaling, cell cycle-related signaling, and metabolisms. Impressively, a higher CCC score was positively correlated with a majority of chemokines, receptors, immunostimulators, and anticancer immunity, indicating a relatively immune-promoting microenvironment. In addition, GSE173839, GSE25066, GSE41998, and GSE194040 dataset analyses of the underlying CCC signature suggested that durvalumab with olaparib and paclitaxel, taxane-anthracycline chemotherapy, neoadjuvant cyclophosphamide/doxorubicin with ixabepilone or paclitaxel, and immunotherapeutic strategies might be suitable for COAD patients with higher CCC score. Eventually, the GDSC database analysis showed that lower CCC scores were likely to be more sensitive to 5-fluorouracil, bosutinib, gemcitabine, gefitinib, methotrexate, mitomycin C, and temozolomide, while patients with higher CCC score seemed to have a higher level of sensitivity to bortezomib and elesclomol.ConclusionThe novel CCC signature exhibited a good ability for prognosis prediction for COAD patients, and the CCC score was found to be highly correlated with molecular features, immune-related characteristics, and therapeutic responses, which would greatly promote clinical management and precision medicine for COAD.

Project description:CCL17 is an important chemokine that plays a vital immunomodulatory role in the tumor microenvironment (TME). Analysis of lung adenocarcinoma (LUAD) data in Kaplan-Meier plotter databases found that the overall survival of patients in the CCL17 high-expression group was higher than that of the low-expression group, especially for patients with early (stages I and II) LUAD, which has a more positive prognostic value. Expression of CCL17 in LUAD was positively correlated with the proportion of tumor-infiltrating lymphocytes, immunostimulators, and major histocompatibility complexes using the TISIDB databases. Based on the RNA-seq and clinical data of 491 LUAD patients obtained from the TCGA database, 1,455 differential genes were found between the CCL17 high- and low-expression groups. Using WGCNA analysis confirmed that the expression of differential genes in the blue module is negatively correlated with poor survival and clinical stages of LUAD patients, and CCL17 and CCR4 genes belong to the hub genes in the blue module. Further analysis by the ESTIMATE and CIBERSORT algorithm found that the naive B cells and CD8+ T cells in the CCL17 high-expression group have a higher distribution ratio in the early LUAD patients, and the high immune score has a positive relationship with the overall survival rate. Using somatic mutation data of TCGA-LUAD, we found that 1) the tumor mutation burden values of the CCL17 high-expression group were significantly lower than those of the CCL17 low-expression group and 2) the expression levels of CCL17 and the tumor mutation burden values were negatively correlated. Transwell chemotaxis and cytotoxicity assays confirmed that CCL17 contributes to the migration of CCR4-positive lymphocytes into the H1993 LUAD TME and enhances the specific lysis of LUAD cells. In summary, high expression of CCL17 in the LUAD TME promotes local immune cell infiltration and antitumor immune response, which may contribute to the better survival and prognosis of patients with early LUAD.

Project description:IntroductionAs the molecular features of lung adenocarcinoma (LUAD) have been evaluated as a cross-sectional study, the course of tumor characteristics has not been modeled. The temporal evolution of the tumor immune microenvironment (TIME), as well as the clinico-molecular features of LUAD, could provide a precise strategy for immunotherapy and surrogate biomarkers for the course of LUAD.MethodsA pseudotime trajectory was constructed in patients with LUAD from the Cancer Genome Atlas and non-small cell lung cancer radiogenomics datasets. Correlation analyses were performed between clinical features and pseudotime. Genes associated with pseudotime were selected, and gene ontology analysis was performed. F-18 fluorodeoxyglucose positron emission tomography images of subjects were collected, and imaging parameters, including standardized uptake value (SUV), were obtained. Correlation analyses were performed between imaging parameters and pseudotime. Correlation analyses were performed between the enrichment scores of various immune cell types and pseudotime. In addition, correlation analyses were performed between the expression of PD-L1, tumor mutation burden, and pseudotime.ResultsPseudotime trajectories of LUAD corresponded to clinical stages. Molecular profiles related to cell division and natural killer cell activity were changed along the pseudotime. The maximal SUV of LUAD tumors showed a positive correlation with pseudotime. Type 1 helper T (Th1) cells showed a positive correlation, whereas M2 macrophages showed a negative correlation with pseudotime. PD-L1 expression showed a negative correlation, whereas tumor mutation burden showed a positive correlation with pseudotime.ConclusionThe estimated pseudotime associated with the stage suggested that it could reflect the clinico-molecular evolution of LUAD. Specific immune cell types in the TIME as well as cell division and glucose metabolism were dynamically changed according to the progression of the pseudotime. As a molecular progression of LUAD, different cellular targets should be considered for immunotherapy.

Project description:BackgroundMicropapillary adenocarcinoma is one of the most aggressive histologic subtypes of lung adenocarcinoma (LADC), and even a minor proportion of micropapillary component (MPC) within the LADC could contribute to poor prognosis. Comprehensive analysis of genetic and immunological features of LADC with different percentages of MPC would help better understand cancer biology of this LADC subtype and direct future treatments.MethodsWe performed next-generation sequencing (NGS) for a discovery cohort of 43 LADC patients whose tumors were micro-dissected to separate MPC and non-MPC lesions and a reference cohort of 113 LADC patients. MPC-enriched genetic alterations that were detected in the discovery cohort were then confirmed using a validation cohort of 183 LADC patients. Immunological staining was also conducted on the MPC-containing samples in the discovery cohort.ResultsTumors with a higher percentage of MPC tended to harbor more tumor mutation burdens (TMBs) and chromosome instability (CIN). Some rare genetic events may serve as the genetic landscape to drive micropapillary tumor progression. Specifically, alterations in transcription termination factor 1 (TTF1), brain-specific angiogenesis inhibitor 3 (BAI3), mammalian target of rapamycin (MTOR), and cyclin-dependent kinase inhibitor 2A (CDKN2A) were cross-validated to be enriched in MPC-contained LADC. Additionally, tumors with a higher percentage of MPC were associated with a higher percentage of CD4+, CD8+, and PD-L1+ staining, and some genetic changes that were enriched in MPC, including MET amplification and MTOR mutation, were correlated with increased PD-L1 expression.ConclusionWe identified multiple novel MPC-enriched genetic changes that could help us understand the nature of this aggressive cancer subtype. High MPC tumors tended to have elevated levels of TMBs, T cell infiltration, and immunosuppression than low MPC tumors, implying the potential link between MPC content and sensitivity to immunotherapy.

Dataset Information

Prognosis and Immunological Characteristics of PGK1 in Lung Adenocarcinoma: A Systematic Analysis.

Publications

Prognosis and Immunological Characteristics of PGK1 in Lung Adenocarcinoma: A Systematic Analysis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets