人胰腺癌空间和单细胞转录组数据的整合分析揭示与不良预后相关的中间态癌细胞群体

Integrative analysis of spatial and single-cell transcriptome data from human pancreatic cancer reveals an intermediate cancer cell population associated with poor prognosis

作者信息Seongryong Kim, Galam Leem, Junjeong Choi, Yongjun Koh, Suho Lee, Sang-Hee Nam, Jin Su Kim, Chan Hee Park, Ho Kyoung Hwang, Kyoung Il Min, Jung Hyun Jo, Hee Seung Lee, Moon Jae Chung, Jeong Youp Park, Seung Woo Park, Si Young Song, Eui-Cheol Shin, Chang Moo Kang, Seungmin Bang, Jong-Eun Park
PMID38297291
发布时间2024-01-31
DOI10.1186/s13073-024-01287-7

实验完整度

包含单细胞RNA测序、空间转录组、RNA原位杂交、免疫组化、公共数据集生存分析及转录因子网络分析等多个独立证据层级,并包含功能验证和机制探讨。

主要模型

人胰腺癌组织(17例,含PDAC和IPMN) 非免疫细胞(CD45阴性)单细胞RNA测序 对应组织切片的空间转录组(7例)

重点核对

患者来源:Severance Hospital,2019年1月至2020年7月,共17例 单细胞测序:CD45阴性细胞,Chromium Next GEM Single Cell 5' v1.1,目标细胞数10,000-15,000 空间转录组:Visium Spatial Gene Expression for FFPE,7例患者 生存分析:TCGA和ICGC队列,基于优化截断值划分高低表达组 RNA原位杂交及免疫组化验证关键标记基因表达

摘要

背景:近期利用单细胞转录组分析的研究报道了胰腺癌肿瘤微环境中多个不同的肿瘤性上皮细胞和癌症相关成纤维细胞簇。然而,由于肿瘤内和肿瘤间的异质性,其分子特征和生物学意义尚未被明确阐明。方法:我们利用来自17例胰腺肿瘤组织(16例胰腺癌和1例高度异型增生)的富集非免疫细胞群体进行了单细胞RNA测序,并生成了来自7例患者样本的配对空间转录组数据。结果:我们鉴定出五个不同的功能性胰腺癌细胞亚簇和六个不同的癌症相关成纤维细胞亚簇。我们深入分析了它们的特征,并发现这些亚簇成功地将胰腺癌 bulk 转录组分析中提出的大多数特征去卷积。在这些亚簇中,我们鉴定了一个新的癌细胞亚簇 Ep_VGLL1,尽管其具有预后价值,但显示出介于 basal-like 和 classical 二分法两个极端之间的中间特性。Ep_VGLL1 的分子特征提示其具有 basal-like 和 classical 亚型之间的过渡特性,空间转录组数据支持了这一观点。结论:这项整合分析不仅提供了胰腺癌和成纤维细胞群体的全面图谱,还为胰腺癌细胞的动态状态提供了新的见解,并揭示了潜在的治疗靶点。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
胰腺癌肿瘤微环境中癌细胞和成纤维细胞的异质性及其分子特征和生物学意义是什么?
核心机制
鉴定出新的癌细胞亚簇 Ep_VGLL1,其具有介于经典和基底样亚型之间的过渡特性,可能通过SOX4和KLF5转录因子活性参与SMAD4缺陷的TGF-β驱动的致癌过程,并与Wnt非依赖性相关。
主要证据
单细胞RNA测序鉴定出五个癌细胞亚簇,通过转录因子活性分析、RNA原位杂交、空间转录组和公共队列生存分析,发现Ep_VGLL1与不良预后相关,且其空间位置介于经典和基底样亚型之间。
研究意义
研究为胰腺癌细胞的动态状态提供新见解,并指出靶向Ep_VGLL1以中断经典到基底样亚型转变可能成为克服化疗耐药的新治疗策略。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

单细胞转录组数据生成与预处理

获取胰腺癌组织中非免疫细胞的单细胞转录组数据,为后续亚群分析提供基础。

从17例手术切除的胰腺肿瘤组织中制备单细胞悬液,磁珠分选去除CD45阳性细胞,构建单细胞RNA文库并测序。

2

上皮细胞和成纤维细胞亚群鉴定与注释

鉴定胰腺癌细胞和成纤维细胞的亚簇,并确定其分子特征。

对上皮细胞和成纤维细胞进行聚类分析,使用标记基因命名亚簇,并通过拷贝数变异推断、细胞周期分析等区分癌细胞和非癌细胞。

3

患者分层与公共数据验证

基于细胞亚簇组成对患者进行分层,并验证亚簇的预后价值。

基于癌细胞和成纤维细胞亚簇比例进行层次聚类,得到三个患者群,并在TCGA和ICGC队列中评估亚簇标记基因的预后价值。

4

转录因子网络分析

揭示各癌细胞亚簇的转录因子调控网络及其特征。

使用SCENIC推断转录因子活性,构建相关网络,并识别各亚簇的特异性转录因子群。

5

空间转录组定位

解析癌细胞亚簇和成纤维细胞在组织中的空间分布及相互作用。

对7例患者样本进行空间转录组测序,利用单细胞数据作为参考进行去卷积,并分析细胞类型间的邻域富集关系。

6

体外功能轨迹分析

验证Ep_VGLL1在经典到基底样转变中的过渡性质。

利用已发表的FOLFIRINOX处理的肿瘤球单细胞数据,追踪亚簇标记随时间的变化。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
gentleMACS C管Miltenyi Biotec--
人肿瘤解离试剂盒Miltenyi Biotec--
gentleMACS Octo解离器Miltenyi Biotec--
死细胞去除试剂盒Miltenyi Biotec--
人CD45磁珠Miltenyi Biotec--
Chromium Next GEM单细胞5'文库v1.110x Genomics--
NovaSeq 6000测序系统Illumina--
Visium空间基因表达试剂盒(FFPE)10x Genomics--
RNAscope多重荧光试剂盒----
Pannoramic SCAN II3D Histech--
Benchmark XTVentana Medical Systems--
COL9A1抗体LSBioLS-C98645
KRT6B抗体Abcamab154313
SPRR3抗体BiossBS-11163R
TRIM54抗体OrigeneTA803871
Infinium亚洲筛查芯片-24 v1.0Illumina--
VGLL1探针--AD44673
TRIM54探针--AD555211-C3
KRT6B探针--AD805641-C4
KRT19探针--AD310221-C2
PDGFRA探针--AD604481-C3
COCH探针--AD1104401-C1
TSLP探针--AD403541-C2
PI16探针--AD569181
CD34探针--AD560821-C2

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
患者样本收集
患者纳入标准:新诊断、经病理确认的胰腺癌;样本数量:17例;组织类型:肿瘤核心病变;收集时间:2019年1月至2020年7月。
阅读提示:Methods 中的 'Patient information and sample collection method'
单细胞悬液制备
酶消化时间和温度;使用gentleMACS C-tubes和酶混合物(Enzyme H, R, A);细胞过滤:40µm细胞筛。
阅读提示:Methods 中的 'Patient information and sample collection method'
细胞分选
CD45阴性细胞富集;死细胞去除步骤;磁珠分选方法。
阅读提示:Methods 中的 'Cell sorting'
单细胞RNA文库构建
文库类型:Chromium Next GEM Single Cell 5' v1.1;目标细胞数:10,000-15,000;测序平台:NovaSeq 6000;测序深度:约50,000 reads/细胞。
阅读提示:Methods 中的 'Single-cell RNA library construction and sequencing'
单细胞数据预处理
过滤标准(UMI>2000,基因数500-7000,线粒体基因<10%,Scrublet singlet);样本去多重化(Souporcell,k=3或4);批次校正(Harmony,BBKNN)。
阅读提示:Methods 中的 'Data preprocessing procedures'
空间转录组
FFPE组织块;Visium试剂盒;6x6mm区域;测序平台。
阅读提示:Methods 中的 'Spatial transcriptomic data analysis' 和 'Spatial transcriptome library construction and sequencing'
RNA原位杂交
组织切片厚度5µm;脱蜡处理;探针序列(列表见方法);检测系统(Pannoramic SCAN II)。
阅读提示:Methods 中的 'RNA in situ hybridization'
免疫组化
抗体及稀释比(COL9A1 1:200, KRT6B 1:500, SPRR3 1:200, TRIM54 1:150);自动染色设备。
阅读提示:Methods 中的 'Immunohistochemistry'