Decoding per- and polyfluoroalkyl substances (PFAS) in hepatocellular carcinoma: a multi-omics and computational toxicology approach

Abstract Background Per- and polyfluoroalkyl substances (PFAS), particularly perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS), are synthetic chemicals known for their widespread use and environmental persistence. These compounds have been increasingly linked to hepatotoxicity and the development of hepatocellular carcinoma (HCC). However, the molecular mechanisms by which PFAS contribute to HCC remain underexplored. Methods This study employs a multi-omics approach that combines network toxicology, integrated machine learning, single-cell RNA sequencing, spatial transcriptomics, experimental validation, and molecular docking simulations to uncover the mechanisms through which PFAS exposure drives HCC. We analyzed publicly available transcriptomic data from several HCC cohorts and used differential gene expression analysis to identify targets associated with both PFAS exposure and HCC. We constructed a protein–protein interaction (PPI) network and a survival risk model, the PFAS-related HCC signature (PFASRHSig), based on integrated machine learning to identify prognostic biomarkers, with the goal of identifying core targets of PFAS in HCC progression and prognosis. RT-qPCR and immunohistochemical (IHC) staining were used to validate the expression levels of the targets in both tumor and normal tissues. Molecular docking simulations were conducted to assess the binding affinities between PFAS compounds and selected target proteins. Results Functional enrichment studies revealed that PFAS targets were associated with metabolic signaling pathways, which are actively involved in lipid, glucose, drug metabolism, etc. Through integrated machine learning and PPI network analysis, we identified six genes, APOA1, ESR1, IGF1, PPARGC1A, SERPINE1, and PON1, that serve as core targets of PFAS in both HCC progression and prognosis. These targets were further validated via bulk RNA-seq, single-cell RNA-seq, and spatial transcriptomics, which revealed differential expression patterns across various cell types in the HCC tumor microenvironment. The results of RT-qPCR and IHC staining were consistent with the in silico findings. Molecular docking simulations revealed strong binding affinities between PFAS compounds and these core targets, supporting their potential roles in PFAS-induced hepatocarcinogenesis. Conclusions Our study highlights key molecular targets and pathways involved in PFAS-induced liver carcinogenesis and proposes a robust survival risk model (PFASRHSig) for HCC. These findings provide new insights into PFAS toxicity mechanisms and offer potential therapeutic targets for mitigating the health risks associated with PFAS exposure. Collectively, our findings help in advancing clinical applications by providing insights into disease mechanisms and potential therapeutic interventions.