《HUMAN MUTATION》：Single-Cell RNA Sequencing–Integrated Fibroblast Markers and Candidate Drug–Target Associations in Pancreatic Ductal Adenocarcinoma
编辑推荐：
背景：胰腺导管腺癌（PDAC）预后极差，其发生发展受富含成纤维细胞的肿瘤微环境（TME）塑造。目前对癌症相关成纤维细胞（CAFs）及其治疗相关性的细胞类型分辨、模块水平认识仍不充分。方法：研究人员整合了GSE212966中6例未经治疗的PDAC切除标本的单细胞
背景：胰腺导管腺癌（PDAC）预后极差，其发生发展受富含成纤维细胞的肿瘤微环境（TME）塑造。目前对癌症相关成纤维细胞（CAFs）及其治疗相关性的细胞类型分辨、模块水平认识仍不充分。方法：研究人员整合了GSE212966中6例未经治疗的PDAC切除标本的单细胞RNA测序（scRNA-seq）数据；公共元数据中未提供个体肿瘤分期信息。研究进行了标准化预处理、聚类和主要细胞类型注释，分析了向成纤维细胞的接触介导的细胞间信号传导。研究人员对成纤维细胞应用高维加权基因共表达网络分析（hdWGCNA）以识别共表达模块和枢纽基因。通过单因素Cox模型和Kaplan–Meier分析评估与总生存期（OS）的关联；采用CIBERSORT和ESTIMATE评估免疫浸润；通过基因集富集分析（GSEA）检查通路关联。候选化合物通过Enrichr/DSigDB筛选，并使用PyMOL–Chem3D–AutoDockTools工作流程与靶蛋白进行分子对接。进一步在PDAC细胞系中通过siRNA介导的SEMA3C敲低，随后进行qRT–PCR、CCK-8、伤口愈合和Transwell侵袭实验对SEMA3C进行检查。结果：研究人员分辨出10个细胞亚群，并观察到向成纤维细胞的广泛接触信号，包括CD99–CD99轴。hdWGCNA delineated成纤维细胞模块，并从中优先筛选出四个关键模块。五个基因——SEMA3C、PMP22、CPE、SNX9和UACA——与总生存期显著相关。这些生物标志物与特定免疫细胞组分呈正相关和负相关，但与StromalScore、ImmuneScore和ESTIMATEScore呈正相关。高表达组富集于癌症相关程序，如上皮–间充质转化（EMT）、炎症反应和TNFα–NF-κB信号传导，而若干增殖和代谢通路在CPE低表达组中富集。体外实验中，siRNA介导的SEMA3C敲低在PANC-1和Capan-2细胞中降低了细胞活力、伤口闭合和侵袭能力，支持SEMA3C在PDAC细胞中具有细胞内在的促肿瘤作用。Coumestrol成为计算优先筛选的化合物，分子对接预测其与SNX9具有潜在有利相互作用（预测结合能?7.5 kcal/mol）；该关联需要实验验证。结论：一个整合的单细胞和网络框架识别了PDAC中成纤维细胞相关共表达模块和五个生存相关基因。这些基因的预后相关性以及预测的coumestrol–SNX9关联需要进一步的独立验证和实验验证。
胰腺导管腺癌（PDAC）是胰腺癌中最常见的组织学类型，其5年生存率长期处于个位数水平，多数患者确诊时已属不可切除或转移性疾病。尽管围手术期管理、以吉西他滨或FOLFIRINOX为基础的联合化疗方案以及有限的靶向和免疫治疗策略取得了一定进展，但人群水平的生存改善仍然有限，耐药和复发普遍存在。除恶性上皮细胞固有的基因组不稳定性和多克隆演化外，以深度免疫抑制、致密间质屏障和异常血管为特征的肿瘤微环境（TME）已成为制约疗效的主要因素。传统以肿瘤细胞为中心的研究范式难以解释患者间治疗反应的显著差异和获得性耐药的快速出现。因此，从肿瘤–微环境共演化的视角重新审视PDAC，系统刻画TME的细胞组成、相互作用和分子程序，并将其与临床结局和治疗反应相关联，具有重要的科学和转化价值。在此背景下，癌症相关成纤维细胞（CAFs）在PDAC的TME中丰度高且功能关键，通过分泌和重塑细胞外基质（ECM）驱动促结缔组织增生和组织硬化，改变间质压力和液体流动形成药物递送屏障，并通过细胞因子、趋化因子和代谢物重编程血管生成和免疫浸润，促进免疫抑制和免疫逃逸。CAFs还是重要的"信号放大器"，与恶性细胞、内皮细胞和免疫细胞发生广泛的配体–受体相互作用，协调肿瘤生长、侵袭和转移。值得注意的是，CAFs并非均一群体，而是表现出显著的分子和功能异质性，单细胞和空间组学研究已 delineated 肌成纤维细胞样（myCAF）、炎症性（iCAF）和抗原呈递性（apCAF）等表型，这些亚型在空间定位、分泌程序、免疫调节和代谢串扰方面存在差异。多项研究将CAF异质性与患者生存、化疗敏感性和免疫治疗反应性相关联，但不同研究之间在亚群定义、标志物集和功能解释方面仍存在差异。然而，既往研究多聚焦于预定义的CAF亚型或单一标志物，而在PDAC中，将成纤维细胞程序与预后和免疫背景相关联的细胞类型分辨的共表达模块仍表征不足。在PDAC之外，CD99已在CAFs中被检测到，同嗜性CD99结合已被证明可调节β1-整合素/FAK依赖性细胞黏附。然而，CD99–CD99介导的接触是否影响PDAC中CAF的活化仍不清楚。为解决这些空白，研究人员整合了一个公开的PDAC单细胞数据集，并应用高维加权基因共表达网络分析（hdWGCNA）在细胞分辨率下识别成纤维细胞相关共表达模块，随后对模块相关基因进行优先级排序，并在TCGA-PDAC队列中评估其与生存、免疫浸润和功能通路的关联。同时，研究人员结合药物扰动资源和分子对接生成假设水平的候选化合物–靶点关联。研究人员假设成纤维细胞特异性共表达分析能够识别临床相关的间质程序并生成可实验检验的假设。研究的贡献包括一个以CAF为中心的模块水平分析框架、五个成纤维细胞相关生存相关基因以及需要进一步验证的初步计算药物–靶点预测。
在技术方法方面，研究使用Seurat进行单细胞RNA测序数据的质量控制、归一化和聚类分析（质控后保留49,191个细胞），采用SCTransform进行方差稳定化，Harmony校正样本间批次效应，UMAP进行二维可视化。差异表达分析采用FindAllMarkers（Wilcoxon秩和检验）。细胞间通讯分析使用CellChat及其配体–受体知识库。成纤维细胞亚群的高维加权基因共表达网络分析（hdWGCNA）用于构建共表达网络并识别模块和枢纽基因。生存分析采用单因素Cox比例风险模型和Kaplan–Meier方法（TCGA-PDAC队列，129例肿瘤样本）。免疫浸润评估使用CIBERSORT（LM22签名矩阵）和ESTIMATE算法。基因集富集分析（GSEA）使用HALLMARK基因集。体外实验使用PDAC细胞系PANC-1和Capan-2以及非恶性胰腺导管上皮细胞系HPDE6-C7，通过siRNA介导的SEMA3C敲低结合qRT–PCR、CCK-8、伤口愈合和Transwell侵袭实验进行功能验证。候选药物预测通过Enrichr/DSigDB筛选，分子对接采用PyMOL–Chem3D–AutoDockTools工作流程。
研究结果方面，在"PDAC单细胞图谱"部分，研究人员通过质控、归一化、批次校正和降维聚类，在UMAP嵌入中分辨出10个主要细胞亚群，包括T细胞、成纤维细胞、巨噬细胞、上皮细胞、B细胞、平滑肌细胞、内皮细胞、肥大细胞、雪旺细胞和MKI67+祖细胞样细胞，并通过CellMarker 2.0和文献策展的经典标志物完成注释；成纤维细胞特异性上调基因的功能注释显示ECM和基质重塑通路显著富集。在"细胞–细胞通讯分析"部分，通讯分析揭示了PDAC微环境中广泛的跨簇信号网络，在限制为"细胞–细胞接触"类别后，预测到来自T细胞、巨噬细胞和上皮细胞向成纤维细胞的CD99–CD99相互作用。在"hdWGCNA识别成纤维细胞相关模块"部分，研究人员在成纤维细胞亚群中构建加权共表达网络（软阈值功率β = 8），通过动态树切割识别出9个稳健模块（M1–M9），其中M1、M3、M4和M8在成纤维细胞中表现出较高的表达比例和水平，被指定为关键模块，并提取了这些模块的枢纽基因交互网络。在"生存关联分析"部分，通过单因素Cox模型在TCGA-PDAC队列中评估成纤维细胞相关基因的生存关联，发现SEMA3C、PMP22、CPE、SNX9和UACA五个基因与总生存期显著相关，Kaplan–Meier分析显示所有五个基因均具有显著的生存分层效应（log-rank p < 0.05）。在"生物标志物与免疫浸润的相关性"部分，基于CIBERSORT推断的22种免疫细胞类型比例，计算了五个生物标志物与免疫浸润水平的Spearman相关性，多个关联具有统计学显著性，表明与固有免疫和适应性免疫区室均有关联；ESTIMATE分析一致显示所有五个生物标志物与StromalScore、ImmuneScore和ESTIMATEScore呈显著正相关。在"GSEA"部分，以各生物标志物的中位表达量将样本二分为高表达组和低表达组进行GSEA分析，高表达组富集于癌症相关HALLMARK通路，包括上皮–间充质转化、炎症反应和TNFα–NF-κB信号传导等；相反，CPE低表达组富集于E2F靶点、MYC靶点V1/V2和氧化磷酸化通路。在"SEMA3C在PDAC细胞系中的实验验证"部分，qRT–PCR显示SEMA3C、PMP22、CPE、SNX9和UACA在PDAC细胞系PANC-1和Capan-2中的表达水平高于非恶性胰腺导管上皮细胞系HPDE6-C7，其中SEMA3C的升高幅度最大；两种独立的SEMA3C靶向siRNA均显著降低了SEMA3C mRNA表达，选择敲低效率更高的si-SEMA3C#1进行后续功能实验；SEMA3C敲低显著降低了48和72小时的细胞活力、48小时的伤口闭合以及Matrigel包被Transwell实验中的侵袭细胞数，支持SEMA3C在促进PDAC细胞增殖、迁移和侵袭中的细胞内在作用，但未确立SEMA3C的CAF特异性功能。在"药物优先级排序与分子对接"部分，经DSigDB筛选后，由于SNX9具有实验解析的结构（PDB: 3DYT）且存在兼容的候选配体，SNX9–coumestrol对被选为代表性结构评估对象，coumestrol–SNX9复合物的预测结合能为?7.5 kcal/mol，满足预设的计算阈值（≤ ?5.0 kcal/mol），预测的配体–受体接触在3.5 ?以内。
讨论部分对研究结果进行了深入阐释。CD99作为一种同嗜性细胞表面黏附分子，在PDAC之外的研究中已在CAFs中被检测到，其结合与β1-整合素/FAK依赖性黏附的调节相关。本分析预测了T细胞、巨噬细胞、上皮细胞与成纤维细胞之间的CD99–CD99相互作用，这种同嗜性接触可能改变免疫/上皮–成纤维细胞接触的持续时间或强度，进而影响与CAF状态相关的黏附、细胞骨架组织和基质重塑程序。然而，CellChat仅从转录组数据推断配体–受体兼容性，不能确定信号方向、CAF活化状态或因果关系，因此所提出的CD99依赖性机制应被视为假设生成性的。关于SEMA3C，其在PDAC中上调，促进肿瘤细胞存活并重塑免疫环境（包括调节PD-L1和与M2样巨噬细胞标志物的关联），其高表达与较差的生存相关；研究人员在PANC-1和Capan-2细胞中的siRNA介导的敲低实验表明SEMA3C缺失降低了细胞活力、迁移和侵袭，为其在PDAC中的促肿瘤作用提供了直接功能支持。关于PMP22、CPE和UACA，讨论指出这些基因不应被解释为已确立的PDAC CAF驱动因子，其纳入主要基于它们属于成纤维细胞来源的共表达模块以及与生存和间质评分的关联，因此研究结果将这些基因定位为模块相关候选标志物而非成纤维细胞特异性因果功能。在免疫浸润方面，研究人员观察到五个基因与特定免疫细胞组分之间的混合正负相关，同时StromalScore、ImmuneScore和ESTIMATEScore升高，提示免疫富集但可能功能失调的微环境而非均匀增强的抗肿瘤免疫。CIBERSORT仅估计细胞比例而不提供空间信息，因此当前的相关性并非CAF介导T细胞排斥的直接证据。关于coumestrol，近期体外和体内研究表明其在多种恶性肿瘤中表现出抗肿瘤活性，作为天然haspin抑制剂，富含芽的食物饮食在小鼠模型中减少了肠道息肉并改善了恶病质相关全身表型；PLCO的流行病学分析表明较高的总异黄酮和选定成分（包括coumestrol）摄入量与胰腺癌风险降低相关，主要在既往吸烟者中。然而，coumestrol在肿瘤学中仍处于临床前/早期证据阶段，流行病学关联并不能证明生理可达到的血浆或肿瘤组织浓度足以作用于SNX9或发挥抗肿瘤效应，因此对接结果应仅被视为计算假设而非治疗疗效的证据。研究局限性包括：单细胞分析仅纳入6例未经治疗的PDAC切除标本且个体肿瘤分期信息不可用，样本量有限可能无法捕捉CAF状态的完整患者间、分期相关和治疗相关异质性；生存分析使用基因特异性Kaplan–Meier最佳截断点和单变量Cox模型筛选个体生存相关候选基因，未构建多基因风险评分或临床预测模型，未进行多变量调整和独立外部队列验证，因此所识别的基因不应被解释为独立预后因素；免疫浸润景观通过转录组去卷积推断，缺乏空间、影像或蛋白质水平确认；五个基因与CAFs的关联主要基于转录组和网络关联，谱系特异性和因果功能需要在定义的myCAF、iCAF和apCAF群体中验证；SEMA3C敲低实验仅展示了癌细胞内在表型，缺乏CAF–肿瘤细胞共培养实验阻碍了关于旁分泌或CAF特异性SEMA3C功能的结论；分子对接仅提供预测性证据，直接结合实验、分子动力学模拟、细胞靶点结合实验、药代动力学分析以及体内疗效和安全性研究仍有必要。结论部分翻译如下：这项以单细胞为基础的研究表征了PDAC微环境，识别了9个成纤维细胞共表达模块和5个成纤维细胞相关生存相关基因——SEMA3C、PMP22、CPE、SNX9和UACA。细胞间通讯分析识别了一个推定的CD99–CD99接触轴，而免疫和通路分析将候选基因与间质富集和肿瘤相关生物学程序相关联。SEMA3C敲低在体外降低了PDAC细胞的活力、迁移和侵袭。分子对接进一步生成了初步的coumestrol–SNX9相互作用假设。总体而言，这些发现提供了候选生物标志物和可实验检验的机制，但其预后独立性、CAF特异性和药理学相关性需要进一步验证。
