《npj Systems Biology and Applications》：3D tissue model to patient atlas: dynamic network topology reveals PDAC chemoresistance bottlenecks
编辑推荐：
胰腺导管腺癌（PDAC）通过动态细胞状态转变产生化疗耐药，但调控耐药耐受状态的调控网络架构仍表征不足，因为静态表达快照无法捕捉相互作用拓扑如何沿治疗诱导轨迹演变。研究人员在脱细胞猪小肠基质上建立3D PANC-1组织模型，该模型达到类患者增殖指数和组织微环境特
胰腺导管腺癌（PDAC）通过动态细胞状态转变产生化疗耐药，但调控耐药耐受状态的调控网络架构仍表征不足，因为静态表达快照无法捕捉相互作用拓扑如何沿治疗诱导轨迹演变。研究人员在脱细胞猪小肠基质上建立3D PANC-1组织模型，该模型达到类患者增殖指数和组织微环境特异性的药物反应，并在四种条件下进行单细胞RNA测序（scRNA-seq）：对照、吉西他滨（GEM）、转化生长因子-β1（TGF-β1）及其联合。将拟时序轨迹推断与动态蛋白-蛋白相互作用（PPI）网络分析耦合，研究人员沿细胞状态转变连续追踪介数中心性（betweenness centrality），而非比较离散状态。GEM存活细胞启用增殖程序（E2F1、CDK1、AURKA、TPX2、TOP2A、BIRC5），而TGF-β1诱导独立的上皮-间质转化（EMT）吸引子（SPOCK1、MBOAT2、COL5A1、THBS1、FN1）。二者汇聚处出现一个涌现瓶颈：CDK1介数中心性选择性激增，受CDKN1A（p21）和WEE1约束，定义了一个富含GEM存活细胞的S期持续状态（S-phase-persistence state）。在来自231例患者的726,107个细胞的PDAC图谱中交叉检验显示，三阳性细胞（triple-positive cells）的转移富集为8.7倍。该框架将受控3D组织实验、单细胞系统分析和临床规模患者数据联系起来，适用于跨癌种耐药研究。
论文解读
研究背景与问题：胰腺导管腺癌（PDAC）以极高化疗耐药和死亡率为特征，晚期诊断限制手术，GEM单药或联合氟嘧啶类方案的中位总生存约一年，辅助FOLFIRINOX疗效优于GEM单药但安全性较差。GEM应用超过25年，耐药涉及致密间质导致药物渗透差、核苷转运体介导摄取受损及细胞周期阻滞不完全等；即使撤药后检查点仍延迟，提示药物压力下动态细胞状态转变而非静态表达差异决定存活。解析这些转变需要生理相关实验平台和能连续追踪调控架构的分析框架，而传统2D培养和标准差异表达难以实现。研究人员据此建立3D PANC-1组织模型并进行四条件scRNA-seq，耦合拟时序与动态PPI网络拓扑分析，并在患者单细胞图谱中交叉验证，以识别驱动GEM耐药的调控决策点并同时验证模型与方法。
关键技术方法（不超过250字）：研究人员在脱细胞猪小肠黏膜下层（SISmuc，保留黏膜的小肠黏膜下层）上建立3D PANC-1组织模型，设置对照、GEM、TGF-β1及GEM+TGF-β1四条件，采用10x Genomics 3′ v3进行scRNA-seq；以Scanpy聚类、scFates/Palantir拟时序、OmniPath PPI网络、介数中心性、Personalized PageRank及随机敲除模拟分析网络拓扑；在Loveless et al.提供的PDAC单细胞图谱（726,107个细胞、231例患者、12项独立研究）中，采用PHATE、签名评分、Fisher精确检验和ROC分析进行交叉检验。
研究结果：3D PANC-1组织模型揭示GEM耐药的空间与转录异质性。H&E、Ki-67和MTT显示GEM在隐窝龛优先清除敏感细胞，存活增殖细胞经TGF-β1共刺激进一步富集，3D群体整体高耐药。单细胞聚类识别由细胞周期时相和TGF-β1状态定义的四类细胞状态。UMAP与Leiden聚类得到四个稳定亚群，G1期与S/G2/M期富集群分离；GEM存活细胞不形成独立群而分布于既有群，cluster 2（S期+EMT）富集最强（48%），共享复制应激耐受和动粒调控基因签名。拟时序轨迹推断解析连续TGF-β1诱导状态转变并识别GEM相关基因程序。两条平行轨迹从G1期对照亚群走向TGF-β1相关状态，TGF-β1信号和EMT持续上调、上皮基因逐步下调；三组时序基因模块分别对应TGF-β1早期/晚期诱导、G1→S细胞周期变化及二者重叠且富集于GEM存活细胞的程序，并与患者生存相关。轨迹耦合网络拓扑揭示联合S期+TGF-β1状态中涌现的CDK1–CDKN1A瓶颈。连续追踪介数中心性发现CDK1与CDKN1A在联合状态出现协同尖峰，为连接S期/有丝分裂模块与EMT/TGF-β1/SMAD模块的桥梁；复合影响评分和敲除模拟显示CDKN1A、CSNK2A1/CK2α等可卸载该瓶颈，表达逻辑为CDK1由S期推进诱导、CDKN1A由TGF-β1/SMAD3驱动。CDK1–CDKN1A–WEE1瓶颈在转移性PDAC中富集并在726,107细胞患者图谱中预测疾病进展。L2瓶颈签名定位于EMT与增殖分支交界，三阳性细胞转移OR=8.71，CDK1
+/CDKN1A
+共表达从正常上皮0.1%升至转移细胞6.3%，GEM存活者以G1阻滞和WEE1
+富集为主，FOLFIRINOX存活者保留更多cycling细胞。CDKN1A功能耦合EMT与细胞周期程序。在CDK1
+/CDKN1A
+细胞中，EMT与细胞周期评分由负相关转为正相关，CDK1
+/CDKN1A
+增殖细胞上调EMT效应基因，DUAL-program细胞最具转移倾向（OR=3.54）。瓶颈签名优于通用基因集。27基因bridge签名和19基因Leiden 2签名区分转移与初治细胞的AUC分别约0.630和0.620，优于通用增殖、通用EMT和既有耐药panel，三基因轴二元共检测OR=8.71优于连续评分。药物特异性生存策略。GEM存活者CDK1
+和DUAL-program细胞减少而WEE1
+富集，FOLFIRINOX存活者CDK1
+、DUAL-program和WEE1
+均富集，提示GEM以WEE1介导检查点阻滞为主、FOLFIRINOX维持经瓶颈的细胞周期运行；治疗臂来自单一研究且无匹配前样本，属观察性结果。
讨论部分总结：该研究整合生理相关3D组织平台、沿拟时序连续追踪PPI网络拓扑的系统框架及726,107细胞患者图谱交叉验证，识别CDK1–CDKN1A–WEE1瓶颈为人类PDAC复发性特征，并与转移倾向和药物特异性生存策略相关。2D培养增殖指数过高掩盖选择动力学；SISmuc模型约35% Ki-67阳性接近临床均值，GEM后Ki-67升高反映耐药增殖亚群选择，与Striefler观察一致。GEM存活细胞富集S/G2/M群及CDK1、BIRC5、TOP2A、TPX2、AURKA、RRM2；TGF-β1诱导FN1、THBS1、COL5A1、SPOCK1、MBOAT2并丢失KRT18/KRT19，S100A2、FN1、MBOAT2具预后意义。单独靶向TGF-β1或EMT获益有限，关键在于两程序共发生。动态网络框架揭示CDK1/CDKN1A协同尖峰仅见于联合状态，S期经E2F1/MYBL2/CCNA2驱动CDK1，TGF-β1/SMAD3以p53非依赖方式诱导p21，形成S期持续并由RAD51AP1、BRIP1、RRM2、CLSPN→CHK1维持。患者图谱中三阳性细胞定位EMT–增殖分支交界，CDKN1A共表达逆转EMT–细胞周期负相关，Werba结果亦支持化疗富集EMT程序。治疗上，GEM存活者呈WEE1主导检查点阻滞，FOLFIRINOX存活者维持cycling，提示ATR–CHK1–WEE1轴为状态依赖性脆弱点；CDK1–CDKN1A–WEE1复合签名可作为候选伴随诊断。局限包括单一PANC-1细胞系、拟时序为推断、网络依赖OmniPath curation完整性及图谱治疗比较为观察性。
研究结论部分翻译：研究人员建立了一个生理相关的3D PDAC组织平台，使药物诱导耐药轨迹的单细胞映射成为可能，并引入了一个系统生物学框架，即拟时序轨迹推断与动态PPI网络拓扑分析整合，能够解析静态方法不可见的涌现调控瓶颈。在726,107细胞患者图谱中对组织模型发现进行交叉检验，同时支持该实验平台和分析框架：在四种确定条件下受控3D系统中识别的签名，在群体规模上重现了人类PDAC最具侵袭性的生物学。CDK1–CDKN1A–WEE1轴定义了一个S期持续状态，其中增殖驱动与TGF-β1诱导的约束汇聚，产生可能 underlying 复发和最小残留病（minimal residual disease）的化疗耐受亚群。三阳性细胞显示8.7倍转移富集，CDKN1A共表达逆转EMT–细胞周期负相关，药物特异性生存模式揭示两种可采用不同治疗靶向的耐药策略。该框架可直接应用于具备单细胞轨迹数据的其他癌种耐药研究。
