一、写在前面文章Human pancreatic cancer single cell atlas reveals association of CXCL10 fibroblasts and basal subtype tumor cells**IF10.4发布于《Clinical Cancer Research》。**文章通过汇总大多数来自人类胰腺的scRNAseq研究(n229)克服了小样本限制。单细胞RNA测序和空间转录组学加深了对胰腺导管腺癌(PDAC)肿瘤微环境(TME)的认识利用scRNAseq图谱发现了一种新的基底肿瘤细胞与CXCL10肌成纤维细胞的浓缩和空间关联。并且证明了利用来自大量患者的基因标志在细胞系和患者衍生的有机化合物中去卷积药物治疗的价值该计算方法可以应用于未来患者在治疗前和治疗后的批量RNAseq活组织检查以揭示治疗耐药的潜在机制。了解更多单细胞和空间转录组课程如果需要单细胞数据分析教学、生信热点全文复现、自测数据个性化分析辅导、实验科研服务和常态化实验学习欢迎联系客服微信[Biomamba_zhushou]。二、主要内容1. 导管亚型揭示了9个转录上不同的群体结合了来自12个不同研究的229个胰腺scRNAseq样本包含172个原发肿瘤(PT)、25个转移灶(Met)、26个邻近正常(AdjN)和6个健康供体胰腺的更大的患者队列以揭示关于PDAC亚型和相关TME的新生物学使用Louvain算法以一种无偏的方式对细胞进行聚类。证实了推测的九个簇的CNS值比从正常邻近组织和供体组织分离的细胞更高簇2和8在供体和AdjN样品中显著富含而簇1、4和7在PT和Met中富含。虽然一些肿瘤主要由单一导管亚型组成通常是群体1或2但大多数肿瘤是异质性的包含多个导管群体这与最近的发现一致。这些发现突出了先前研究中确定的肿瘤内的异质性。2.ScRNAseq肿瘤亚型确定了两个转录上不同的基础群体使用ReactomePA我们研究了在每个簇中丰富的途径以确定潜在的功能差异。Ductal 1和Ductal 3都富含与金属转运蛋白相关的途径Ductal 3与基质金属蛋白酶活性有关Ductal 9丰富了涉及细胞代谢的途径。Ductal 6的基因特征与内质网应激反应一致Ductal 5表现出干扰素信号特征。Ductal 7丰富了与细胞周期调控有关的途径。在Ductal 4中观察到与细胞外基质组织和c-Met信号相关的途径丰富。TP63作为基底细胞状态的驱动因素在Ductal 4和7中富含。**3.两个转录上不同的scRNAseq基础基因特征与PDAC总体生存率较差相关使用未调整的Cox比例风险模型Ductal 4和7基础亚群基因特征与较差的OS一致相关Ductal 5与更差的OS相关。在Puleo数据集中Ductal 1、8和9与更好的OS一致相关而在TCGA-PAAD中Ductal 1与更好的OS相关。Ductal 2与TCGA-PAAD中较好的生存相关。总体而言这些研究表明基础浓缩样品的OS更差而经典浓缩样品的OS更好。4. 与PDAC相比在正常胰腺中补体基因丰富的成纤维细胞增多而单细胞图谱来源的肌成纤维细胞基因特征与较差的存活率相关。炎性(ICAF)、肌成纤维细胞(MyCAF)和抗原提呈(ApCAF)三种转录上不同的成纤维细胞群体具有不同的功能和空间分布。分离成纤维细胞并将它们单独聚集在一起确定了8个转录上不同的簇其中五组肌成纤维细胞和三组炎性成纤维细胞。经路径分析IFB 1和IFB 2被富含补体信号IFB-3群体显示与补体和ECM途径都相关而且唯一地表达了高水平的VEGFA。myFb1-5群体高度丰富了与ECM信号和胶原组装相关的通路MMP11和COL11A1在myFb 1中的表达最高MyFb2细胞富含与ATP合成、柠檬酸循环和低氧信号有关的通路提示这些细胞可能是先前描述的代谢群体myFb 3富含CXCL14一种先前被认为促进肿瘤细胞侵袭特性的趋化因子。除此之外确定了一组独特的表达myFbs(MyFb5)的CXCL10高度富含干扰素伽马反应基因(IFNG)。研究了成纤维细胞评分与总存活率的关系myFb集群与较差的OS有关正常/邻近的胰腺中“炎症性”成纤维细胞丰富而肌成纤维细胞在肿瘤形成过程中增多。IFb和myFb群体之间共享的免疫调节功能存在显著重叠。肌成纤维细胞基因特征的高表达可能与PDAC的整体预后较差相关。除此之外鉴定了几个主要的髓系细胞树突状细胞通常与较好的OS相关而中性粒细胞一般与较差的OS相关单核细胞和巨噬细胞评分与存活率没有明显的相关性可能是由于抗原提呈树突状细胞启动了更有效的T细胞反应而中性粒细胞在PDAC中具有高度的免疫抑制作用。5. 将单细胞图谱 scRNAseq基因特征投影到批量RNAseq数据集上显示CXCL10myFb与基底肿瘤细胞一致共浓缩基底细胞丰富的PDAC患者中观察到的更糟糕的OS对四个PDAC批量RNA-seq数据集进行评分以估计细胞组成。对所有细胞类型进行了成对的皮尔逊相关然后进行了分级聚类2号管与1号管、3号管和9号管一致相关这表明它也可能是一个经典亚型两个杂交星团之一的Ductal 5也与Ductal 4/7有关增强了该群体处于基础状态或从基础状态过渡的可能性。Ductal 4和7一致地与肌成纤维细胞聚集在一起特别是在所有四个批量数据集中的myFb5。MyFb5s之所以有别于其他Fb群体部分原因是它表达趋化因子CXCL10在scRNA-seq数据中验证了基础簇和表达肌成纤维细胞簇的CXCL10之间的正相关性。CXCL10受体CXCR3在髓系TNK细胞以及肿瘤细胞中表达还在除Ductal 8外所有肿瘤亚群中表达。此外调节CXCL10表达的IFNG在T和NK细胞中高表达而其受体IFNGR1和IFNGR2在成纤维细胞上都高表达。6. 空间转录最近邻分析显示基底细胞与肌成纤维细胞非常接近22个空间转录PDAC患者肿瘤切片进行了去卷积在所有空间PDAC组织中确定了scRNAseq中每种细胞类型的区域。强调了肿瘤内的异质性也证实了scRNAseq鉴定的导管种群在完整的组织中被发现。使用去卷积ST分析对于每个样本对每种导管细胞类型的邻居进行了计数发现Ductal 4与myFbs显著相关特别是那些高表达MMP11和CXCL10的细胞。高Ductal 4/7含量与CD8T细胞的存在呈负相关这表明在基底丰富的肿瘤中潜在的细胞毒性T细胞的渗透减少这可能是这些患者预后较差的部分原因。7. PDAC体外2D和3D模型的基因签名分析显示使用图谱获得的签名具有Ductal 4占优势的特征应用CIBERSORTx反卷积来估计Ductal簇的比例。发现簇4在所有2D细胞系中都高度表达反卷积结果表明在KRAS突变的PDAC中ERK抑制或KRAS沉默细胞系67降低了Ductal 4的含量增加了Ductal 3和Ductal 5的含量除了Panc-1和SW1990外所有的细胞系都从基底样表型转变为经典表型。进一步评估了患者来源的有机类物质批量RNAseq数据集并鉴定出一个样本带有Ductal 9一个样本带有Ductal 7这两个样本都不存在于2D细胞系数据中突出了2D和3D细胞培养之间的差异。此外在正常样本中一致地鉴定了Ductal 8反映了供体和邻近正常样本的富集性。数据与空间反卷积一致表明Ductal 4是主要的亚型。且Ductal 8可能是反应性导管细胞和早期Panin的细胞成分。8. 人PDAC中scRNAseq和空间相关性的验证显示原位有不同的过渡性基底细胞免疫荧光对组织进行Panck染色以定位导管/上皮细胞并联合染色与导管4(LAMC2)或导管7(S100A2)相关的标记。观察到Ductal 4亚型的高患病率。与肿瘤细胞相比正常导管上皮中LAMC2和S100A2的表达存在显著差异并且细胞形态向提示EMT/基础状态的更间充质纺锤形转变与在scRNAseq的UMAP可视化中类似可以看到过渡状态。对PDAC组织进行多种染色发现Ductal 1和CXCL10成纤维细胞之间缺乏亲和力观察到CXCL10附近的Ductal 4和Ductal 7区域表达SMACAFmyFb和Ductal 4含量较高的区域也显示出肿瘤细胞内的CXCL10增强了局部TME对IFNG信号反应的可能性。三、总结与展望一、研究核心总结本文通过整合多源单细胞 RNA 测序scRNAseq数据构建了包含 726,107 个高质量细胞的 “人类胰腺癌细胞单细胞图谱”并结合空间转录组与实验验证完成了三方面关键工作1. 细胞分型与亚型定位将胰腺癌细胞细分为 9 个亚群识别出高表达 KRT17/KRT5 的 “基底亚型肿瘤细胞”证实其与患者低生存率显著相关同时将基质细胞分为 5 个成纤维细胞亚群与 2 个内皮细胞亚群发现高表达 CXCL10 的 “CXCL10 成纤维细胞”。2. 空间关联验证通过单细胞共定位分析、空间转录组技术确认基因表达区域重叠及多重免疫荧光 / RNA 原位杂交实验明确 CXCL10 成纤维细胞与基底亚型肿瘤细胞在胰腺癌组织中 “近距离共存”。3. 机制初步探索预测并提出 “CXCL10-CXCR3 信号轴” 是两者相互作用的关键介导通路推测该轴可通过激活肿瘤细胞增殖通路、招募免疫抑制细胞双重促进肿瘤进展。二、研究创新与意义一核心创新点1.数据整合与图谱价值首次整合 5 个公开数据集、39 例多类型样本原发灶、转移灶、正常组织的 scRNAseq 数据构建的 “胰腺癌细胞单细胞图谱” 规模大、覆盖全为后续 PDAC 单细胞研究提供了 “基准参考地图”。2.“单细胞 空间” 的双重验证突破单一 scRNAseq “脱离空间位置” 的局限通过空间转录组与实验技术结合直接证实细胞亚群的空间关联让 “细胞互作” 从 “数据推测” 变为 “可视化证据”。3.精准定位关键细胞互作对首次发现 “CXCL10 成纤维细胞 - 基底亚型肿瘤细胞” 这一特定互作组合并关联到患者预后为 PDAC “细胞亚群层面的精准研究” 提供了新靶点。二临床与科研意义1.预后诊断价值明确 “基底亚型肿瘤细胞占比” 可作为 PDAC 患者预后判断的标志物为临床通过检测相关基因如 KRT17/KRT5评估患者生存风险提供依据。2.治疗靶点拓展提出的 “CXCL10-CXCR3 信号轴” 为 PDAC 治疗提供新方向 —— 可通过抑制该轴如阻断 CXCL10 与 CXCR3 结合同时抑制肿瘤细胞增殖与免疫抑制微环境形成为联合治疗提供思路。3.研究方法示范展示了 “scRNAseq 构建图谱→细分细胞亚群→空间转录组验证关联→信号通路推测” 的完整研究框架为其他实体瘤的 “单细胞 空间” 研究提供了可借鉴的范式。三、研究局限与未来方向一主要局限1.数据与样本的局限性依赖公开数据集整合缺乏自主采集的大样本验证尤其是不同种族、不同治疗背景的患者样本可能存在数据偏倚且未纳入胰腺癌术前治疗如化疗、靶向治疗后的样本无法分析治疗对 “CXCL10 成纤维细胞 - 基底亚型肿瘤细胞” 互作的影响。2.机制验证的初步性对 “CXCL10-CXCR3 轴” 的功能作用仅停留在 “预测与推测”未通过细胞实验如敲降 CXCL10 或 CXCR3 后观察肿瘤细胞行为变化或动物模型验证该轴对肿瘤进展的实际影响机制证据不够充分。3.空间分辨率的不足所用空间转录组技术的分辨率有限无法精确到 “单个细胞” 层面难以判断 CXCL10 成纤维细胞与基底亚型肿瘤细胞是否存在 “直接接触” 及具体的接触方式如通过细胞黏附分子还是分泌因子。二未来研究方向1.样本拓展与验证扩大样本量纳入不同临床背景如不同分期、不同治疗方案的患者样本验证 “基底亚型预后价值” 与 “CXCL10-CXCR3 轴作用” 的普适性开展多中心研究减少数据偏倚。2.机制深度验证通过体外细胞共培养实验将 CXCL10 成纤维细胞与基底亚型肿瘤细胞共培养、体内动物模型如构建过表达 CXCL10 的 PDAC 小鼠模型验证 “CXCL10-CXCR3 轴” 对肿瘤增殖、侵袭及免疫微环境的调控作用进一步探索该轴下游的具体信号分子明确调控网络。3.技术升级与多组学整合采用更高分辨率的空间转录组技术如单细胞空间转录组或结合空间蛋白组技术精准解析两种细胞的空间接触模式整合基因组、表观基因组数据分析 “CXCL10 成纤维细胞” 与 “基底亚型肿瘤细胞” 的基因变异差异探索细胞亚群形成的分子根源。4.临床转化探索基于 “CXCL10-CXCR3 轴” 开发靶向药物如 CXCR3 拮抗剂开展早期临床试验评估其对 PDAC 患者的治疗效果探索将 “基底亚型标志物” 与 “CXCL10 表达水平” 结合开发联合诊断试剂盒用于 PDAC 的预后分层与治疗方案选择。