Intratumor microbiome derived glycolysis-lactate signatures depicts immune heterogeneity in lung adenocarcinoma by integration of microbiomic, transcriptomic, proteomic and single-cell data.

Abstract

Microbiome plays roles in lung adenocarcinoma (LUAD) development and anti-tumor treatment efficacy. Aberrant glycolysis in tumor might promote lactate production that alter tumor microenvironment, affecting microbiome, cancer cells and immune cells. We aimed to construct intratumor microbiome score to predict prognosis of LUAD patients and thoroughly investigate glycolysis and lactate signature's association with LUAD immune cell infiltration. The Cancer Genome Atlas-LUAD (TCGA-LUAD) microbiome data was downloaded from cBioPortal and analyzed to examine its association with overall survival to create a prognostic scoring model. Gene Set Enrichment Analysis (GSEA) was used to find each group's major mechanisms involved. Our study then investigated the glycolysis and lactate pattern in LUAD patients based on 19 genes, which were correlated with the tumor microenvironment (TME) phenotypes and immunotherapy outcomes. We developed a glycolysis-lactate risk score and signature to accurately predict TME phenotypes, prognosis, and response to immunotherapy. Using the univariate Cox regression analysis, the abundance of 38 genera were identified with prognostic values and a lung-resident microbial score (LMS) was then developed from the TCGA-LUAD-microbiome dataset. Glycolysis hallmark pathway was significantly enriched in high-LMS group and three distinct glycolysis-lactate patterns were generated. Patients in Cluster1 exhibited unfavorable outcomes and might be insensitive to immunotherapy. Glycolysis-lactate score was constructed for predicting prognosis with high accuracy and validated in external cohorts. Gene signature was developed and this signature was elevated in epithelial cells especially in tumor mass on single-cell level. Finally, we found that the glycolysis-lactate signature levels were consistent with the malignancy of histological subtypes. Our study demonstrated that an 18-microbe prognostic score and a 19-gene glycolysis-lactate signature for predicting prognosis of LUAD patients. Our LMS, glycolysis-lactate score and glycolysis-lactate signature have potential roles in precision therapy of LUAD patients.