一区IF9.9全套复刻模板|单细胞+空间转录+硅内扰动,肝癌肿瘤进化与靶向靶点完整闭环

细胞癌（HCC）异质性极强，门静脉癌栓、远处淋巴结转移预后极差，现有单细胞研究仅局限单一组织样本，缺少从正常肝→原发灶→癌栓→转移灶完整多维度细胞进化图谱。本研究整合多队列单细胞、10X Visium空间转录、bulk临床大队列，搭载cNMF恶性细胞转录程序拆解、多类基质亚群分群、细胞通讯、空间生态位、Geneformer硅内敲除、药敏与生存验证六层生信分析体系；系统定义4类肝癌核心恶性meta-program，阐明内皮/成纤维/巨噬构成促转移基质生态位，筛选HSP90B1为跨细胞通用耐药靶点，全程依托公共测序数据集，仅纯计算生信无细胞动物实验，消化/肝脏肿瘤一区高分标准可直接复用完整范式！

大规模单细胞分析与计算机模拟扰动揭示了肝细胞癌从起始到治疗靶向的动态演变过程

研究思路解构

多阶段单细胞整合：正常肝、原发HCC、PVTT门静脉癌栓、淋巴结转移四组样本Harmon批次校正，全局细胞分群注释；
恶性细胞cNMF分解：CopyKAT区分肿瘤/正常细胞，cNMF聚类得到4类核心恶性转录meta-program；
基质细胞精细分型：内皮、成纤维、髓系巨噬亚群细分，对比原发/转移细胞丰度差异；
细胞互作分析：CellChat全局配受体，对比原发、转移灶通讯强度变化；
空间转录验证：Cell2location细胞空间反卷积，COMMOT空间通讯定位基质生态位；
硅内扰动筛选：Geneformer虚拟敲除三类关键细胞，筛选共有核心致病基因； 7. 临床转化验证：TCGA/索拉非尼治疗队列验证靶基因表达、空间分布、预后价值。

主要结果

1.全流程单细胞研究设计+全局细胞分群图谱

左侧示意图完整展示样本分组、测序平台、全套分析算法流水线，清晰说明本研究覆盖正常、原发、癌栓、转移四类病变组织，实现肝癌进展全阶段细胞图谱绘制。右侧UMAP图完成7大类细胞无监督分群，配套标记基因气泡图区分肝细胞、内皮、成纤维、髓系、淋巴细胞谱系。下方堆叠柱状图对比四类组织细胞占比，髓系、恶性肝细胞在转移样本显著扩增，淋巴细胞在肿瘤组织大幅减少，直观体现进展期免疫抑制微环境特征。

2.cNMF解析4类肝癌恶性细胞核心转录Meta-program

左上热图对15个基础基因程序分层聚类，最终收敛为4套稳定肿瘤转录模块。四幅UMAP分色图分别展示Diff-Metabolic、Prolif-Stress、MYC-Biosynth-Immune、EMT-Inflammatory四类细胞在样本内分布。下方柱状图显示癌栓、转移灶富集增殖、EM型恶性细胞，原发灶分化型细胞占比更高；配套marker基因气泡图区分四类细胞特征分子，为后续预后分层提供分子基础。拟时序轨迹可见EMT细胞为肿瘤起始状态，可双向分化为代谢型、增殖型细胞，证明肝癌细胞具备高度可塑性。

3.四类恶性程序转录调控、预后与药物敏感性分析

上方pySCENIC转录因子热图，四类恶性细胞拥有独立主控转录因子，仅ATF3、FOSB少量调控因子共有，解释细胞状态转换的分子基础。中间四组KM曲线基于TCGA队列，Diff-Metabolic组生存期显著更长，增殖、MYC、EM亚型均提示不良预后，完美匹配四类细胞恶性表型。下方scTherapy药敏预测图展示四类细胞药物偏好，分化型对索拉非尼天然耐药，EM亚型对MEK、BCL2抑制剂敏感，为分层精准用药提供计算依据。

4.髓系巨噬细胞亚群异质性与免疫功能解析

UMAP图将髓系细胞细分为巨噬、单核、树突、肥大细胞，其中巨噬进一步拆分为10个功能亚群。分组柱状图显示Macro-SPP1、Macro-TREM2抑制型巨噬在癌栓、转移灶大量富集，正常肝以MARCO库普弗细胞为主。多通路富集热图对比巨噬功能，SPP1/TREM2亚群高表达糖酵解、EMT、血管生成通路，CXCL9型巨噬富集IFN抗肿瘤信号。GSEA对比两类抑制巨噬，SPP1亚群免疫抑制、促转移通路激活程度更强，是肝癌核心免疫抑制细胞。

5.内皮、成纤维细胞亚型、预后与通路富集

A/B幅UMAP分别完成13种内皮、7种成纤维亚群分群，标记基因气泡图清晰区分各亚型特征分子。森林图筛选Endo-ESM1、Fib-POSTN为不良预后基质细胞，高表达患者生存显著缩短。下方通路堆叠图显示ESM1内皮、POSTN成纤维高富集TGFβ、血管生成、ECM重塑通路；伪时序轨迹证实正常成纤维可逐步分化为促转移基质亚群，肿瘤微环境信号驱动基质重编程。

6.原发vs转移灶细胞通讯网络差异对比

左右全局通讯图分别展示原发、门静脉癌栓/淋巴结转移的配受体互作强度，转移灶内皮、成纤维双向通讯信号大幅增强。分群差异热图量化各细胞互作增减幅度，基质细胞与恶性EM细胞、抑制巨噬的交流显著上调，与T细胞通讯全面减弱。下方配受体气泡图锁定核心轴：Endo-ESM1与Fib-POSTN通过FN1、ANGPT家族形成基质交互环路，同步激活肿瘤侵袭信号，构成转移生态位核心通讯。

7.10X Visium空间转录原位验证基质生态位

多例患者组织空间反卷积图，Prolif-Stress恶性细胞占据肿瘤核心，Endo-ESM1、Fib-CD36、POSTN成纤维集中分布于肿瘤间质交界区域。COMMOT空间信号矢量图显示ANGPTL信号集中在基质富集区，MISTy建模证实间质区域TGFβ、缺氧通路显著上调。局部ROI分层展示空间排布：中心基质核心包裹外周增殖/EM肿瘤细胞，Macro-SPP1巨噬环绕外围，形成分层促转移微环境空间架构。

8.Geneformer硅内扰动锁定核心靶点HSP90B1

左上Geneformer三套细胞状态虚拟敲除，交集得到11个共调控关键基因。DepMap数据库CRISPR依赖评分证明HSP90B1是肝癌细胞生存必需基因。下方多组箱图显示HSP90B1在肿瘤、索拉非尼免疫耐药样本显著高表达；空间分布图可见治疗无应答者基因弥漫全瘤，应答者仅局灶表达。末尾生存曲线证实HSP90B1高表达提示短生存期，从计算预测、空间分布、临床预后多层验证其作为耐药、预后双靶点价值。

易选团队总结

顶刊完整多层证据链：多阶段单细胞→恶性转录程序→基质分型→细胞通讯→空间原位→AI虚拟扰动→临床药敏预后，完全匹配npj系列一区收稿标准，无分析短板；
双重核心创新：首创四分类肝癌恶性meta-program进化模型，依托AI硅内扰动不依赖湿实验筛选可用药靶点，区别普通单细胞浅度分群；
复刻成本极低：全部单细胞、空间、TCGA公开数据集，无需临床测序、细胞动物造模，纯生信即可完成全部主体内容；
通用性极强：肝胆、胰腺、结直肠等实体瘤，仅替换单细胞样本即可完整复用“单细胞+空间+虚拟扰动”全套分析管线；
临床转化价值突出：HSP90B1可作为肝癌预后标志物与联合治疗靶点，四类恶性程序可用于术前分层指导靶向/免疫方案选择。

Seurat/cNMF/CellChat/Geneformer均为开源成熟工具，整套流程无冷门算法。不会文中生信分析的科研人员可一站式承接，全套代码、标准化绘图快速落地；胰腺癌、结直肠癌、胆管癌等实体瘤可直接替换单细胞数据集复用整套逻辑，审稿证据链完整返修极少。