《Journal of Cellular and Molecular Medicine》：Single-Cell Characterization of Terminal States and State-Specific Transcriptional Regulatory Networks in Hepatocellular Carcinoma
编辑推荐：
**第三部分：论文摘要翻译**
肝细胞癌（HCC）在肿瘤进展过程中表现出显著的细胞异质性和动态状态转变，然而肿瘤细胞状态沿命运轨迹的涌现及其转录调控机制尚不清楚。在此，研究人员利用公开可用的单细胞RNA测序（scRNA-seq）数据，建立了一个整合细胞命运动态
**第三部分：论文摘要翻译**
肝细胞癌（HCC）在肿瘤进展过程中表现出显著的细胞异质性和动态状态转变，然而肿瘤细胞状态沿命运轨迹的涌现及其转录调控机制尚不清楚。在此，研究人员利用公开可用的单细胞RNA测序（scRNA-seq）数据，建立了一个整合细胞命运动态与转录调控网络的分析框架。首先通过推断的拷贝数变异（CNV）识别恶性细胞，随后应用CellRank进行命运推断，并使用SCENIC（单细胞调控网络推断与聚类）重建调控网络。研究人员鉴定了三种HCC终末细胞状态——免疫激活型、代谢型和增殖型肝细胞——它们形成了从干细胞样到更成熟表型的分化连续体。每种状态由具有特定核心调控因子（包括IRF3、PPARA和JUN）的独特调控网络所控制。结合癌症基因组图谱-肝细胞癌（TCGA-LIHC）bulk转录组和临床数据分析发现，增殖状态来源的转录特征与更差的总体生存相关。综上，本研究提供了一个连接HCC肿瘤细胞命运轨迹、调控异质性与临床结局的机制框架。
**第四部分：论文主体内容解读**
肝细胞癌是全球范围内导致癌症相关死亡的主要原因之一，在中国其疾病负担尤为沉重。现有研究虽然利用单细胞RNA测序（scRNA-seq）技术揭示了HCC的细胞组成，但大多局限于对细胞类型或亚群的静态描述，缺乏对肿瘤细胞在连续转录空间内状态转变的系统分析，特别是不同细胞状态潜在的转录调控网络架构及其关键调控因子尚未被全面阐明。因此，深入探究HCC细胞状态异质性及其背后的调控机制，对于理解疾病发生发展至关重要。为此，研究人员整合单细胞转录组数据与细胞命运推断、转录调控网络重建技术，从动态和机制层面系统刻画了HCC的细胞状态异质性。研究发表于《Journal of Cellular and Molecular Medicine》。
本研究主要采用了以下关键技术方法：数据来源于公共数据库（GSE149614）的scRNA-seq数据，包含来自10名HCC患者的原发肿瘤、门静脉癌栓、转移淋巴结和癌旁组织的样本。分析流程包括：使用inferCNV基于推断的拷贝数变异区分恶性与非恶性细胞；应用Harmony算法对多样本数据进行批次校正；利用CellRank模型推断细胞命运轨迹并识别终末状态；采用SCENIC框架重建转录因子-基因调控网络；使用AUCell算法评估单细胞水平的分化程序活性；进行基因本体（GO）和京都基因与基因组百科全书（KEGG）通路富集分析；并利用癌症基因组图谱-肝细胞癌（TCGA-LIHC）队列数据进行临床关联性分析。
研究结果揭示了HCC细胞在单细胞分辨率下的三个具有不同分化轨迹的终末细胞状态。CellRank分析识别出三个稳定的终末状态（Fate 0, 1, 2），UMAP可视化及批次校正效果评估证实了数据整合的有效性。命运轨迹分析显示HCC细胞通过多分支路径向不同终点演化。这些终末状态并非随机分布，而是对应功能分化的亚群，其恶性程度呈现从Fate 1（免疫激活型）到Fate 0（代谢型）再到Fate 2（增殖型）的递增梯度。
终末细胞状态沿肝细胞分化和干性程序呈现连续梯度。通过评估五个肝细胞分化相关基因程序（HB1, HB2, FH1, FH2, AH）的活性发现，Fate 2细胞富集于早期肝母细胞样特征（HB1/HB2），Fate 1细胞富集于胎儿肝细胞特征（FH1/FH2），而Fate 0细胞则表现出更高的成年肝细胞（AH）成熟程序活性。这表明三个终末状态代表了肝细胞分化谱系上的不同位置，从高干性状态过渡到功能成熟的肝细胞样表型。
HCC细胞的转录调控网络组织成功能上不同的调控模块。利用SCENIC框架重建调控网络并计算连接特异性指数（CSI）后，发现调控子（regulon）呈现出显著的模块化组织结构，模块内部调控子活性高度协调，而模块间则存在功能分离。这表明HCC细胞的转录调控具有高度组织性和模块性，细胞状态异质性源于多个功能不同调控模块的协调活动。
调控模块在终末细胞状态中呈现状态特异性活性模式。分析显示，同一调控模块内的大多数调控子在特定终末命运中选择性地高表达，表明转录因子驱动的调控程序与特定细胞命运紧密关联。部分调控子的活性在命运间呈现渐进变化，反映了命运演化的连续性特征。
不同的终末细胞状态由状态特异性的转录调控模块控制。通过整合SCENIC重建的网络与CellRank推断的命运驱动基因，研究人员构建了精细化的命运特异性调控网络。Fate 0的网络主要与模块3相关，核心调控因子包括IRF3、IRF6；Fate 1的网络与模块1和3富集，核心调控因子包括PPARA、NR1I3、CEBPD；Fate 2的网络则主要集中在模块2、4和6，核心调控因子包括EGR2、ZNF467、HOXA5。这表明每个终末命运并非由单一主调控因子控制，而是由协调行动的、命运特异的转录调控模块所驱动。
终末细胞状态的调控模块表现出不同的功能特化。GO和KEGG功能富集分析表明，Fate 0细胞主要富集于免疫和炎症反应相关通路；Fate 1细胞富集于核心的肝细胞代谢功能，如氨基酸、脂质和碳水化合物代谢；Fate 2细胞则显著富集于细胞周期、DNA复制、核糖体生物发生等与增殖和转录活性相关的通路。这清晰地体现了三个终末状态各自的功能特性。
增殖型（Fate 2）调控特征在HCC中具有临床相关性。研究人员构建了基于Fate 2调控网络中26个转录因子的特征签名，并在TCGA-LIHC队列中评估其临床意义。生存分析显示，Fate 2特征评分较高的患者总体生存较差。在单变量Cox回归中，Fate 2特征评分每增加一个标准差，死亡风险增加约88%（HR=1.88）。虽然评分在不同AJCC分期中无显著差异，但该结果提示状态特异性的转录程序能提供超越传统分期的预后信息。
总结讨论部分指出，本研究整合单细胞命运推断与转录调控网络分析，建立了统一框架解析肝细胞癌的终末细胞命运。研究识别出三种沿连续分化谱系分布的独特终末命运，它们各自由特异的调控模块和核心转录因子所定义。例如，Fate 0的核心调控因子IRF3是cGAS-STING信号通路的核心组件，与免疫激活相关；Fate 1的核心调控因子如PPARA、NR1I3等支持其代谢稳定和成熟肝细胞表型；Fate 2的网络则由HOX家族等转录因子主导，强化了其干性维持和恶性增殖特征。研究结论部分强调，本研究从动态和调控角度为肝细胞癌的细胞异质性提供了机制性见解，并凸显了命运特异性调控程序作为潜在精准治疗靶点的价值。
