Abstract
Lung adenocarcinoma (LUAD) is a major subtype of non-small-cell lung cancer and accompanies high mortality rates. While the role of bilirubin metabolism in cancer is recognized, its specific impact on LUAD and patient response to immunotherapy needs to be elucidated. This study aimed to develop a prognostic signature of bilirubin metabolism-associated genes (BMAGs) to predict outcomes and efficacy of immunotherapy in LUAD. We analysed gene expression data from The Cancer Genome Atlas (TCGA) to identify survival-related BMAGs and construct a prognostic model in LUAD. The prognostic efficacy of our model was corroborated by employing TCGA-LUAD and five Gene Expression Omnibus datasets, effectively stratifying patients into risk-defined cohorts with marked disparities in survival. The BMAG signature was indeed an independent prognostic determinant, outperforming established clinical parameters. The low-risk group exhibited a more favourable response to immunotherapy, highlighted by increased immune checkpoint expression and immune cell infiltration. Further, somatic mutation profiling differentiated the molecular landscapes of the risk categories. Our screening further identified potential drug candidates preferentially targeting the high-risk group. Our analysis of critical BMAGs showed the tumour-suppressive role of FBP1, highlighting its suppression in LUAD and its inhibitory effects on tumour proliferation, migration and invasion, in addition to its involvement in cell cycle and apoptosis regulation. These findings introduce a potent BMAG-based prognostic indicator and offer valuable insights for prognostication and tailored immunotherapy in LUAD.