Establishing a Novel Gene Signature Related to Histone Modifications for Predicting Prognosis in Lung Adenocarcinoma

Epigenetic modifications have been revealed to play an important role in tumorigenesis and tumor development. This study aims to analyze the role of histone modifications and the prognostic values of histone modifications in lung adenocarcinoma (LUAD). The promoters and enhancers of protein encoding genes (PCGs) were the regions of enriched histone modifications.

Expression profiles and clinical information of LUAD samples were downloaded from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Histone modification data of LUAD cell lines were downloaded from Encyclopedia of DNA Elements (ENCODE) database. Limma R package was used to identify differentially expressed PCGs. To identify molecular subtypes, consensus clustering was conducted based on the expression of dysregulated PCGs with abnormal histone modifications. Univariate Cox regression analysis and stepwise Akaike information criterion (stepAIC) were utilized to establish a prognostic model.

We identified a total of 699 epigenetic dysregulated genes with 122 of them significantly correlating with LUAD prognosis. We constructed three molecular subtypes (C1, C2, and C3) based on the 122 prognostic genes. C2 had the longest overall survival while C1 had the worst prognosis. In addition, three subtypes had differential immune infiltration and the response to immune checkpoint inhibitors. Moreover, we identified a risk model containing 5 epi-PCGs that had favorable performance to predict prognosis in different datasets.

This study further supported the critical histone modifications in LUAD development. Three subtypes may provide guidance for the immunotherapy of LUAD patients. Importantly, the prognostic model had great potential to predict LUAD prognosis.