PTM-Related Signatures Predict Lymph Node Metastasis and Shape Microenvironment in Hepatocellular Carcinoma: A
Zhu Ai1, Yuying Liang1, Yuzhi Cai2
1Department of Radiology, The Affiliated Panyu Central Hospital, Guangzhou Medical University, Guangzhou, 511400, People's Republic of China.
Background:
Hepatocellular carcinoma (HCC) has a poor prognosis with lymph node metastasis (LNM), severely impacting survival. Protein post-translational modifications (PTMs) regulate tumor microenvironment (TME) dynamics, but their role in HCC LNM is unknown.
Methods:
This study systematically analyzed single-cell RNA sequencing (scRNA-seq) and bulk transcriptomic data from TCGA-LIHC and validation cohorts. Key cellular subpopulations and differentially expressed genes (DEGs) were identified. PTM-related candidate genes were identified through intersection analysis. Prognostic genes were determined via Cox regression and integrated into a random survival forest (RSF) model. Patients were stratified to assess associations with TME remodeling, genomic alterations, immune checkpoints, and drug sensitivity. Pseudotime trajectory and transcription factor activity analyses elucidated molecular dynamics.
Results:
Hepatocyte populations were significantly expanded in both LNM and non-LNM HCC groups. The intersection of scRNA-seq DEGs, bulk DEGs, and PTM-related genes identified five candidate genes, among which PJA1 and RFPL4B identified as potential prognostic risk factors and were incorporated into a robust RSF model. High-risk patients exhibited poor survival outcomes, increased tumor mutational burden, and pronounced tumor microenvironment alterations. PJA1 expression was negatively correlated with hepatocyte infiltration. Pseudotime analysis revealed that LNM-associated hepatocytes were arrested in early differentiation states, coinciding with dysregulated transcription factor activity.
Conclusion:
This study computationally identified PJA1 and RFPL4B as post-translational modification-related candidate prognostic genes potentially associated with HCC LNM. The proposed risk model stratified patients into subgroups characterized by distinct tumor microenvironment features and genomic instability. Computational analysis suggested that arrest of hepatocyte differentiation and dysregulation of transcription factor networks may be implicated in metastatic progression, thereby offering novel therapeutic insights.


