H3K4me3-related lncRNAs signature and comprehensive analysis of H3K4me3 regulating tumor immunity in lung adenocarcinoma

Backgroud The role of epigenetic modifications in tumorigenesis has been widely reported. However, the role and mechanism of H3K4me3 modification in lung adenocarcinoma (LUAD) are rarely reported systematically. We, therefore, sought to analyze the characteristics of LUAD associated with H3K4me3 modification, build an H3K4me3-lncRNAs score model to predict the prognosis of patients with LUAD and clarify the potential value of H3K4me3 in immunotherapy of LUAD. Methods We evaluated H3K4me3-lncRNA patterns and H3K4me3-lncRNA scores of 477 LUAD samples based on 53 lncRNAs closely correlated to H3K4me3 regulators and comprehensive analyzed the role of these patterns in tumorigenesis and tumor immunity. Using Gene set variation analysis (GSVA), we systematically evaluated the H3K4me3 level of every sample and deeply analyzed the effect of H3K4me3 on the prognosis of LUAD. In addition, we included two independent immunotherapy cohorts to study the impact of high H3K4me3 score on the prognosis of patients. We also used an independent cohort with 52 matched paraffin specimens of LUAD to verify the impact of high H3K3me3 expression on the prognosis of patients. Results We identified three H3K4me3-lncRNA patterns that exhibited specific immune characteristics. Characterized by immunosuppressive and increased TGFβ-mediated epithelial-mesenchymal transition (EMT), patients with high H3K4me3-lncRNA score had a poor overall survival and decreased H3K4me3 score. H3K4me3 score was significantly positively correlated with CD4 + T-cell and CD8 + T-cell activation, programmed cell death and immune checkpoints (ICs) expression, and was negatively correlated with MYC pathway, TP53 pathway, and cell proliferation. Patients with high H3K4me3 score showed elevated expression of ICs, potentiated CD4 T-cell and CD8 T-cell activation, increased programmed cell death, and suppressed cell proliferation and TGFβ-mediated EMT. Patients with high H3K4me3 score and high expression of CTLA4, ICOS, TIGIT, PDCD1LG2, IDO1, CD274, PDCD1, LAG3, or HAVCR2 had the best survival advantage. Two independent immunotherapy cohorts verified that patients with high H3K4me3 score showed an increased inflamed tumor microenvironment (TME) phenotype and enhanced anti-PD-1/L1 immunotherapy response. Immunohistochemistry (IHC) data from 52 matched paraffin specimens of LUAD confirmed that the protein level of H3K4me3 in tumor was significantly lower than that of paracancerous tissues and H3K4me3 brought significant survival benefits to patients with LUAD. Conclusions We build an H3K4me3-lncRNAs score model to predict the prognosis of patients with LUAD. More importantly, this study revealed characteristics of H3K4me3 modification in LUAD and clarified the important potential role of H3K4me3 on tumor immunotherapy and patients’ survival.