泛素化相关的导管-成纤维细胞串扰塑造胰腺导管腺癌的肿瘤进展与预后

《HUMAN MUTATION》：Ubiquitination-Associated Ductal–Fibroblast Crosstalk Shapes Tumor Progression and Prognosis in Pancreatic Ductal Adenocarcinoma

　　

背景：胰腺导管腺癌（Pancreatic Ductal Adenocarcinoma, PDAC）是一种以复杂肿瘤微环境（Tumor Microenvironment, TME）为特征的高度恶性肿瘤。泛素化调节关键的致癌过程和微环境重塑，然而其在PDAC微环境中的细胞类型特异性活性和空间组织仍知之甚少。材料与方法：整合了PDAC的单细胞RNA测序（single-cell RNA sequencing, scRNA-seq）、空间转录组学（spatial transcriptomics）和批量RNA测序（bulk RNA-seq）数据集进行多组学分析。利用基因集评分算法量化泛素化活性，随后进行细胞通讯和空间相互作用分析。利用批量队列构建预后模型，并通过基因沉默功能实验验证关键发现。结果：与邻近正常组织相比，PDAC组织中的泛素化活性显著升高，其中导管上皮细胞和成纤维细胞的升高最为显著。高泛素化状态与增强的导管-成纤维细胞相互作用以及与侵袭相关肿瘤区域相关的独特空间模式有关。整合泛素化相关基因特征确定了一个稳健的预后模型，强调了一种在肿瘤中持续过表达并与不良生存相关的细胞外基质相关因子。功能实验表明，抑制该因子可抑制胰腺癌细胞的增殖、迁移、侵袭和存活。结论：泛素化组织了PDAC微环境内的导管-成纤维细胞串扰，并将空间肿瘤生态学与疾病侵袭性和患者预后联系起来。靶向泛素化相关的微环境程序可能为预后分层和治疗干预提供新策略。

该论文发表于《HUMAN MUTATION》，题为“泛素化相关的导管-成纤维细胞串扰塑造胰腺导管腺癌的肿瘤进展与预后”。研究针对胰腺导管腺癌（PDAC）缺乏早期临床症状、诊断困难且复发率极高，以及泛素化在肿瘤微环境（TME）中作用机制尚不明确的问题，旨在通过多组学整合揭示泛素化如何调控细胞间通讯及影响患者预后。

研究人员采用了多项关键技术方法：首先，整合了来自GEO数据库的单细胞RNA测序数据（GSE212966, GSE197177, GSE235315）和TCGA数据库的批量RNA测序数据；其次，利用AUCell、UCell、Singscore、ssGSEA和AddModuleScore五种算法计算泛素化评分；第三，应用CellChat和mistyR进行细胞通讯与空间相互作用分析；第四，使用RCTD算法对空间转录组数据（stRNA-seq）进行细胞类型反卷积；最后，结合LASSO回归与支持向量机（survivalSVM）构建机器学习预后模型，并通过体外细胞实验验证靶点功能。

研究结果如下：

1. 细胞类型鉴定与泛素化评分

通过对89,208个细胞的分析，识别出16种细胞类型。五种算法评分显示，导管上皮细胞和成纤维细胞在肿瘤组织中的泛素化水平显著高于邻近正常组织。

2. 不同泛素化水平下的细胞通讯与拷贝数变异

基于泛素化中位值将细胞分为高、低两组。CellChat分析表明，高泛素化组的导管细胞与成纤维细胞间通讯强度和频率增加。差异基因富集分析显示，泛素介导的蛋白水解通路在两组细胞中均上调，关键基因包括导管细胞中的ANAPC3、ANAPC7等，以及成纤维细胞中的CUL4A、MDM2等。

3. 空间聚类识别PDAC微环境中的代谢活跃生态位

对两个PDAC组织切片的空间转录组数据进行聚类，发现特定的代谢活跃集群（如切片1的Cluster 2），这些集群在空间上局限于特定组织区域，表明局部代谢重编程。

4. 空间反卷积揭示跨尺度的泛素化相关细胞串扰

空间反卷积分析显示，高泛素化的导管细胞和成纤维细胞具有独特的空间分布模式。mistyR分析进一步揭示，细胞内部（intra view）及15μm半径内的邻近细胞（para_15 view）相互作用对预测结果贡献最大，高泛素化细胞群表现出更强的空间互作网络。

5. ANAPC13⁺导管细胞与HECW2⁺成纤维细胞的空间图谱与互作网络

基于ANAPC13和HECW2的表达将细胞分组。空间映射显示，HECW2⁺成纤维细胞在空间上与具有侵袭相关转录特征的肿瘤细胞区域共定位。COL1A2–ITGB1配体-受体对在肿瘤区域富集，提示其参与细胞外基质介导的信号传导。

6. 预后与免疫治疗相关性

生存分析显示，低浸润ANAPC13⁺导管细胞标记物或HECW2⁺成纤维细胞标记物的患者预后显著更好，且对免疫治疗的反应更佳。

7. 共同差异基因筛选与LAMA3的潜在作用

通过交集分析获得16个共同差异基因，经机器学习筛选出关键基因LAMA3。LAMA3在肿瘤组织中高表达，且与高级别肿瘤和更短的生存期相关，并在多个GEO数据集中得到验证。

8. LAMA3的基因组改变景观

LAMA3突变频率较低，但在LAMA3高表达与低表达组间，常见驱动基因（如TP53、KRAS）的突变谱及拷贝数变异存在差异。免疫浸润分析显示LAMA3改变状态与特定免疫细胞群的分布相关。

9. LAMA3敲低抑制癌细胞恶性表型

体外实验证实，敲低LAMA3可显著抑制SW1990和PANC-1细胞的增殖（CCK-8、克隆形成）、迁移（伤口愈合）、侵袭（Transwell）能力，并促进细胞凋亡（流式细胞术）。

讨论部分指出，该研究综合scRNA-seq、stRNA-seq和常规RNA-seq数据，揭示了泛素化在PDAC中的生物学意义。研究发现导管细胞和成纤维细胞是泛素化介导的蛋白水解通路的关键执行者，并通过空间分析阐明了ANAPC13与HECW2相关细胞亚群的空间动态及互作机制。LAMA3被确定为与不良预后相关的关键因子，其高表达与患者生存率低密切相关。

结论部分总结：本研究揭示了泛素化相关调控程序在PDAC中关键地塑造了导管-成纤维细胞的串扰，从而影响肿瘤进展、侵袭行为和临床结局。导管细胞和成纤维细胞内异常的泛素化活性有助于微环境重塑，并可能构成对治疗反应的异质性基础。值得注意的是，LAMA3相关程序反映了侵袭性肿瘤状态和基质相互作用，凸显了其用于风险分层和治疗靶向的潜在价值。这些发现提供了一个机制框架，将泛素化驱动的细胞-细胞相互作用与疾病进展联系起来，并支持其作为胰腺癌转化生物标志物开发和环境导向治疗策略候选者的相关性。