scRNA-seq and bulk RNA-seq reveal the characteristics of macrophage copper metabolism and establish a risk signature in hepatocellular carcinoma.

BACKGROUND

Hepatocellular carcinoma (HCC) is a prevalent malignancy with an urgent need for improved prognostic stratification and treatment-response prediction. This study aimed to explore a macrophage copper metabolism-associated prognostic model and to investigate the relationship between this risk model and the tumor immune microenvironment.

METHODS

The FindClusters function was used to analyze cell clusters, and CellChat and CellPhoneDB/LIANA were employed for cell-cell communication analysis. Copper metabolism-related genes were sourced from the MSigDB database. A prognostic risk model was established using least absolute shrinkage and selection operator (LASSO) analysis and multivariate Cox regression analysis, and a nomogram was constructed by integrating the prognostic model with clinicopathological factors. Additional analyses were performed to map the seven model genes in single-cell data, assess model uncertainty and robustness, evaluate macrophage/copper/cuproptosis-related transcriptional programs, and examine the correlations between risk score, immune infiltration and predicted drug sensitivity.

RESULTS

Using single-cell RNA sequencing (scRNA-seq) data, we identified four macrophage subpopulations. Macrophages with high SPP1 expression showed close interaction with T cell populations and were associated with copper ion metabolism. By incorporating 141 copper metabolism-related genes and using The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) cohort, we constructed a seven-gene risk prediction model. Additional single-cell mapping showed that the model genes were detectable in the HCC single-cell dataset and showed a macrophage-associated expression pattern. The model showed moderate prognostic discrimination in TCGA-LIHC, whereas its external performance was heterogeneous and remained evaluable across external cohorts, with performance varying among datasets. Immune and mechanism-related analyses suggested that the risk signature was associated with macrophage-related infiltration, copper metabolism and cuproptosis-related transcriptional programs. Drug sensitivity analysis nominated Daporinad as a computationally predicted candidate compound, supporting Daporinad as a pharmacogenomic candidate for follow-up investigation.

CONCLUSIONS

By integrating scRNA-seq and bulk RNA sequencing (RNA-seq) data, we constructed a macrophage copper metabolism-associated prognostic signature for HCC. The risk score was associated with survival, immune microenvironment features and predicted drug response, providing a transcriptomic framework for risk stratification and therapeutic hypothesis generation.