Non-genetic stratification reveals epigenetic heterogeneity and identifies vulnerabilities of glycolysis addiction in lung adenocarcinoma subtype.

Lung adenocarcinoma (LUAD) exhibits high heterogeneity and is well known for its high genetic variation. Recently, the understanding of non-genetic variation provides a new perspective to study the heterogeneity of LUAD. Little is known about whether super-enhancers (SEs) may be primarily responsible for the inter-tumor heterogeneity of LUAD. We used super-enhancer RNA (seRNA) levels of a large-scale clinical well-annotated LUAD cohort to stratify patients into three clusters with different prognosis and other malignant characteristics. Mechanistically, estrogen-related receptor alpha (ERRα) in cluster 3-like cell lines acts as a cofactor of BRD4 to assist SE-promoter loops to activate glycolysis-related target gene expression, thereby promoting glycolysis and malignant progression, which confers a therapeutic vulnerability to glycolytic inhibitors. Our study identified three groups of patients according to seRNA levels, among which patients in cluster 3 have the worst prognosis and vulnerability of glycolysis dependency. We also proposed a 3-TF index model to stratify patients with glycolysis-addicted tumors according to tumor SE stratification.