Abstract
Background: Although the role of tumor microenvironment in lung adenocarcinoma (LUAD) has been explored in a number of studies, the value of TME-related signatures in immunotherapy has not been comprehensively characterized. Materials and Methods: Consensus clustering was conducted to characterize TME-based molecular subtypes using transcription data of LUAD samples. The biological pathways and immune microenvironment were assessed by CIBERSORT, ESTIMATE, and gene set enrichment analysis. A TME-related risk model was established through the algorithms of least absolute shrinkage and selection operator (Lasso) and stepwise Akaike information criterion (stepAIC). Results: Four TME-based molecular subtypes including C1, C2, C3, and C4 were identified, and they showed distinct overall survival, genomic characteristics, DNA methylation pattern, immune microenvironment, and biological pathways. C1 had the worst prognosis and high tumor proliferation rate. C3 and C4 had higher enrichment of anti-tumor signatures compared to C1 and C2. C4 had evidently low enrichment of epithelial-mesenchymal transition (EMT) signature and tumor proliferation rate. C3 was predicted to be more sensitive to immunotherapy compared with other subtypes. C1 is more sensitive to chemotherapy drugs, including Docetaxel, Vinorelbine and Cisplatin, while C3 is more sensitive to Paclitaxel. A five-gene risk model was constructed, which showed a favorable performance in three independent datasets. Low-risk group showed a longer overall survival, more infiltrated immune cells, and higher response to immunotherapy than high-risk group. Conclusion: This study comprehensively characterized the molecular features of LUAD patients based on TME-related signatures, demonstrating the potential of TME-based signatures in exploring the mechanisms of LUAD development. The TME-related risk model was of clinical value to predict LUAD prognosis and guide immunotherapy.