《Molecular Oncology》：Network divergence analysis identifies adaptive gene modules and two orthogonal vulnerability axes in pancreatic cancer
编辑推荐：
本研究运用创新的基因-基因共表达发散框架，首次在胰腺导管腺癌（PDAC）单细胞转录组中系统性地揭示了驱动肿瘤细胞可塑性的四种动态调控模块，并进一步整合功能基因组学数据将其归纳为两个正交的适应性脆弱性轴，为理解PDAC异质性、生存机制及精准治疗提供了全新的网络生物学视角。
在复杂的肿瘤微环境中，癌细胞如何通过重塑其内部的基因调控网络，动态适应外界压力并产生治疗抵抗，是癌症生物学的一大核心问题。尤其是在凶险的胰腺导管腺癌（PDAC）中，肿瘤细胞表现出广泛的异质性和显著的转录可塑性，这不仅使其极具侵袭性，也导致了极差的治疗效果。传统研究方法，如差异表达分析和稳态共表达网络，通常只能捕捉稳定的转录状态或平均的连接模式，难以揭示那些在不同细胞状态间动态重构、驱动适应性行为的基因协调程序。为了突破这一局限，一项开创性研究采用了新颖的“网络发散分析”框架，旨在解码PDAC中隐藏的适应性逻辑。
这项研究分析了来自42个PDAC肿瘤的77,155个恶性细胞的单细胞RNA测序（scRNA-seq）数据。通过整合与批次校正后，研究者识别出10个主要的恶性细胞亚群（Leiden聚类）。研究的关键创新在于其“发散边”算法：它不是简单地比较基因的表达高低，而是计算每个基因对（gene pair）在这些不同细胞亚群中的皮尔逊相关性，并量化这些相关性在亚群间变化的“发散”程度。那些相关性在不同背景下差异最大的基因对，即“发散边”，被认为代表了动态重连的、背景依赖性的基因协调关系，是细胞适应性状态转换的潜在驱动力。
通过构建和可视化这些高发散边构成的网络，并应用社区检测，研究惊人地、可重复地识别出了四个清晰且稳定的适应性转录模块。这四个模块各自代表了一套紧密协调的基因程序，它们突破了经典PDAC状态（如经典型/基底样）的界限，揭示了更深层次的调控结构：
- 1.
集成生长-能量模块：这个模块代表了一个向增强的生物合成和代谢活性转变的转录程序。其核心基因富含核糖体蛋白和线粒体氧化磷酸化组件，功能富集于翻译、细胞呼吸、内质网蛋白加工和mRNA剪接等过程。该模块的启动与MYC和E2F转录因子网络密切相关，其标志性靶点如CDK4和LDHA也显示出与核心基因的强共表达。有趣的是，研究者还发现染色质重塑因子CHD9和生长因子NRG4作为负相关的外周调节因子，提示该生长程序受到内在反馈回路的约束。
- 2.
应激适应性转录模块：该模块定义了一个在经受代谢、蛋白质毒性或环境压力的PDAC细胞亚群中激活的协调基因表达程序。其核心特征是与整合应激反应通路紧密相关，包含ATF4、DDIT3等关键效应因子。该模块对氨基酸饥饿、自噬、缺氧等特定应激类型敏感，但与非葡萄糖饥饿通路关联较弱。模块中包含大量长链非编码RNA和一组与DNA损伤反应而非DNA复制更相关的组蛋白基因，提示了染色质水平的应激调控机制。Notch信号组件PSENEN作为负调节因子出现，暗示了Notch信号可能抑制过度的应激反应激活。
- 3.
IL-2关联的免疫逃逸模块：此模块代表了一种独特的、与细胞因子白介素-2相关的程序，在恶性PDAC细胞内部协调表达多种免疫调节基因。深入分析显示，其富集于IL-2信号通路，并伴随IL2受体亚基及下游JAK3/STAT5A/B的共表达。更重要的是，该模块核心基因的表达与多种免疫检查点分子的表达呈强正相关，包括CTLA4、LAG3、TIGIT等，同时还与免疫抑制细胞因子TGFB1相关。这表明该模块可能通过协调表达多个免疫抑制分子，帮助肿瘤细胞逃避免疫系统的攻击。
- 4.
多通路集体侵袭模块：这是四个模块中最大的一个，其程序支持一种保留细胞连接的协同侵袭模式。功能富集分析指向肌动蛋白丝组织、膜运输、细胞迁移、黏着斑形成等过程。该模块整合了多种信号通路，包括Rho GTP酶、受体酪氨酸激酶、PI3K、TGF-β和Notch，其激活剂和抑制剂呈现协调的正相关，表明其维持的是一种精细的多通路信号平衡，而非单一通路的过度激活。与单细胞/间质侵袭标志物呈负相关，而与集体侵袭/上皮标志物呈正相关，进一步证实了其功能特性。自分泌的EGFR-EREG/TGFA配体-受体对的发现，为自我强化的EGFR信号环路维持定向集体侵袭提供了机制支持。
仅仅识别这些模块并未完整揭示其生物学意义。研究的另一大突破在于，将这些从单细胞数据中衍生出的模块活动评分，与功能基因组学数据库（DepMap）中的CRISPR-Cas9基因依赖性数据和药物敏感性数据进行整合分析。这一整合从根本上重塑了对这些模块的理解，揭示出一个更高层次的组织架构：这四个适应性模块并非四个独立的平行状态，而是坍缩为两个正交的脆弱性轴。
- •
生物合成-代谢轴：主要由IGE模块构成。高IGE活性的细胞系强烈依赖于线粒体、核糖体和RNA加工相关基因，并显示出对PARP抑制剂（如他拉唑帕尼）的选择性敏感，这与MYC/E2F驱动的、依赖持续翻译和氧化磷酸化的生物合成状态所伴随的复制压力和DNA修复需求相吻合。
- •
应激-免疫-侵袭轴：由SAT、IL2和MPC三个模块汇聚而成。尽管它们在单细胞层面代表不同的生物学过程，但在功能依赖性上却共享特征。SAT模块依赖内质网蛋白折叠和蛋白质稳态调节因子；IL2模块关联自噬和细胞因子信号相关因子；MPC模块则依赖细胞极性、细胞骨架重塑和NF-κB信号相关因子。在药物反应上，这三个模块均表现出对加剧蛋白质毒性负荷或干扰细胞骨架/信号传导的化合物的相对敏感性，而对多种DNA损伤剂显示出相对抵抗。
这种“两轴”架构仅在将发散边模块映射到功能基因组学空间时才显现出来，凸显了整合多维度数据在理解肿瘤生物学中的强大力量。
这些适应性程序同样具有临床相关性。在TCGA-PAAD队列的生存分析中，研究者发现SAT模块的高活性与较差的总体生存期相关，MPC模块也显示出相似的不良趋势。当将SAT和MPC模块的活动合并考虑时，定义了一个明确的“应激-侵袭”不良预后轴。相比之下，IGE模块活动与总体风险无显著关联，但一个由最优切点定义的IGE高活性亚组却显示出生存期适度改善的趋势。IL2模块在整体队列中未显示出独立的预后分层能力，这可能与其仅局限于少数恶性细胞亚群有关。在调整了病理分期等临床病理学协变量的多变量模型中，模块评分未能保持独立的预后意义，表明这些适应性程序可能是影响疾病进展阶段的上游调控状态，而非独立的生物标志物。
这项研究通过聚焦于基因-基因协调关系的动态重连，而非静态的表达水平，成功地绘制了PDAC细胞适应性的核心调控图谱。所识别的四个适应性模块及其衍生的两个脆弱性轴，不仅加深了我们对PDAC异质性和可塑性的理解，更重要的是，它们将转录程序与实验可验证的基因依赖性和药物敏感性直接联系起来，为基于肿瘤特定适应性状态的精准治疗策略提供了全新的理论框架和潜在的靶点线索。未来的研究需要进一步在空间维度、实验扰动和纵向治疗背景下验证这些模块的因果角色和动态演变，以最终实现将这些网络生物学见解转化为改善患者预后的有效疗法。
