UCASpatial算法构建与模拟数据基准测试
评估算法在模拟空间转录组数据中反卷积细胞亚群的准确性和对低丰度、转录异质性亚群的识别能力。
从CRC scRNA-seq生成750个模拟spot，按细胞组成熵分低/中/高复杂度，比较UCASpatial与SPOTlight、RCTD、cell2location、CARD、stereoscope、Spotiphy的RMSE、PCC、F1；并分析不同细胞数和不同亚群数下的稳健性。
关于丁香通
公司信息
个人用户
企业机构
无忧采购 轻松科研
- 移动端
超精度空间转录组反卷积解码组织中的免疫异质性与命运决定程序
Ultra-precision deconvolution of spatial transcriptomics decodes immune heterogeneity and fate-defining programs in tissues.
实验完整度
高含模拟数据基准、人CRC空间转录组、TCGA/MSK队列分析及小鼠耳损伤模型，并有mIHC、流式及抗体阻断功能验证
主要模型
人结直肠癌10x Visium空间转录组样本 人结直肠癌scRNA-seq参考数据 小鼠耳穿孔伤口愈合模型（MRL与C57BL/6） TCGA-COAD与MSK-IMPACT队列
重点核对
模拟ST数据：750个spot，按细胞组成熵分低/中/高复杂度各250个，每spot含5–20个细胞，20个参考细胞亚群 人CRC发现队列：7例患者8张肿瘤切片，34个细胞亚群scRNA-seq参考 验证队列CRC8的ST与连续切片mIHC对比，12个免疫亚群比例低于10% 小鼠耳伤口模型：MRL与B6，2 mm耳穿孔，取day 0、3、7、15样本进行ST和scRNA-seq 抗体阻断：anti-IL11RA 2.5 µg/wound于伤后day 3局部注射；rosiglitazone 3 mg/kg隔日注射
摘要
阐明复杂组织内免疫细胞的空间组织和功能特化仍具挑战性。我们提出 UCASpatial，一种利用基于熵的加权来精确映射细胞亚群的超精度空间转录组反卷积算法。基准测试证实其在识别低丰度细胞亚群和区分转录异质性细胞亚群方面的优越性。将 UCASpatial 应用于人结直肠癌，我们揭示单个癌细胞克隆中的染色体 20q 获得 orchestrate 一个 T 细胞排斥的微环境，与 HERV-H 沉默和受损的 I 型干扰素反应相关。在小鼠伤口愈合模型中，我们揭示区分瘢痕性和再生性表型的时空动态。具体而言，我们在瘢痕愈合小鼠（C57BL/6）中鉴定出一个由 Igfbp5+ 软骨细胞、Cd36+Gpnmb+Il1b- 巨噬细胞和 Fmod+ 成纤维细胞组成的促纤维化群落。我们进一步证明该三联体中的 IL11-IL11RA 信号驱动促纤维化群落的形成并限制再生。总之，UCASpatial 可作为破译精细细胞景观和探索复杂动态微环境中细胞间机制的多功能工具。
实验结论
提炼研究问题、关键发现与证据，快速把握文章的核心贡献。
- 研究问题
- 如何实现复杂组织中低丰度、高异质性免疫细胞亚群的超精度空间反卷积，并据此解析肿瘤免疫微环境与伤口愈合中的细胞群落及命运决定程序。
- 核心机制
- chr20q-gain的CRC肿瘤细胞伴随HERV-H表达下调及I型干扰素反应受损，与T细胞排斥性克隆免疫微环境相关；小鼠中IL11-IL11RA轴驱动Igfbp5+软骨细胞、Cd36+Gpnmb+Il1b−巨噬细胞与Fmod+成纤维细胞促纤维化群落的形成并限制再生。
- 主要证据
- 模拟ST数据与Visium HD伪spot基准显示UCASpatial在RMSE、PCC、F1上优于SPOTlight、RCTD、cell2location、CARD、stereoscope、Spotiphy；人CRC ST结合mIHC验证低丰度T细胞亚群注释；TCGA-COAD与MSK-IMPACT支持chr20q-gain与HERV-H降低及ICB后生存相关；小鼠耳伤口模型中anti-IL11RA与rosiglitazone处理减少促纤维化群落并改善愈合。
- 研究意义
- UCASpatial为破译精细细胞景观和探索复杂动态微环境中的细胞间机制提供通用工具；chr20q-gain被提出作为T细胞排斥的基因组标志物及预测ICB反应的免疫基因组参数。
研究路径
按研究推进顺序梳理实验设计、验证步骤与关键观察。
多平台与Visium HD伪spot验证
检验UCASpatial在不同分辨率平台和更接近真实空间相关结构的合成数据上的适用性。
应用于人CRC的10x Visium HD数据（8 μm×8 μm bin）并与RCTD比较；分析小鼠海马Slide-seq V2数据；由三张Visium HD CRC切片聚合49 μm×49 μm bin生成39,813个伪spot进行基准测试。
人CRC空间免疫微环境解析与mIHC验证
验证UCASpatial在真实复杂肿瘤组织中注释低丰度免疫亚群的准确性，并刻画CD8+ T细胞亚群空间分布。
对7例患者8张CRC切片进行ST，使用含34个细胞亚群的scRNA-seq参考反卷积；在CRC8连续切片进行mIHC验证12个免疫亚群；在CRC9和CRC10中用mIHC量化CD8+ T细胞亚群与肿瘤距离及浸润比例。
克隆分辨率免疫微环境分型与多细胞群落鉴定
在克隆分辨率上关联肿瘤基因组改变与T细胞富集或排斥的免疫微环境，并识别相关多细胞群落。
构建UCASpatial-clonalTIME平台，通过空间CNV推断肿瘤克隆并估计每个克隆的免疫微环境组成；按CD4+和CD8+ T细胞亚群比例将13个克隆分为T细胞富集（3个）和T细胞排斥（10个）；计算22个免疫亚群和7个非免疫亚群邻域平均比例的PCC并进行层次聚类，鉴定3个多细胞群落。
chr20q-gain与HERV-H沉默及I型干扰素反应关联分析
探究与克隆免疫表型转变相关的肿瘤遗传改变及其潜在分子机制。
比较克隆间CNV，发现T细胞排斥克隆富集chr20q-gain、chr20p-gain和chr7p-gain；在CRC4中比较clone3与clone1/2的差异基因和转座元件表达；利用TCGA-COAD 471例样本分析CNV与免疫亚群特征、HERV-H表达及I型干扰素通路的关系；用MSK-IMPACT 71例ICB治疗患者分析chr20q状态与生存。
小鼠耳伤口愈合时空动态与巨噬细胞亚群空间隔离分析
在真实伤口愈合ST数据中验证UCASpatial对免疫细胞亚群时空动态和空间分布的解析能力。
采用MRL和B6小鼠2 mm耳穿孔模型，收集day 0、3、7、15耳组织进行10x Visium ST，使用day 7耳组织scRNA-seq（23个细胞亚群）作为参考；用mIHC验证6种主要细胞类型的分布和比例；分析巨噬细胞亚群的空间隔离模式。
促纤维化三细胞群落的鉴定与功能干预
确定与B6小鼠瘢痕愈合相关的三细胞群落，并验证其对愈合表型的因果贡献。
计算23个细胞亚群PCC并层次聚类，鉴定由Igfbp5+软骨细胞、Fmod+成纤维细胞和Cd36+Gpnmb+Il1b−巨噬细胞组成的促纤维化群落；用rosiglitazone阻断Fmod+成纤维细胞分化，或用anti-IL11RA抗体阻断IL11-IL11RA轴，检测群落形成和耳孔剩余面积。
研究方法
按研究目的归类文中使用的方法，便于定位所需技术。
产品清单
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 组织运输与保存 | RPMI 1640培养基 | Gibco | -- |
| HEPES缓冲液 | Thermo Fisher | 15630080 | |
| 青霉素-链霉素 | Thermo Fisher | 15140122 | |
| 两性霉素B | Sigma | S1636 | |
| 硫酸庆大霉素 | Gibco | 15710-064 | |
| 组织包埋 | OCT包埋剂 | Sakura Tissue-Tek | -- |
| 组织消化 | 人肿瘤解离试剂盒 | Miltenyi Biotec | 130-095-929 |
| gentleMACS组织解离仪 | Miltenyi Biotec | 130-093-235 | |
| 细胞过滤 | 70微米细胞筛网 | Biologix | 15-1070 |
| 流式分选 | CD45-FITC抗体 | BioLegend | 368508 |
| 单细胞文库构建 | Chromium Next GEM单细胞5'文库与凝胶珠试剂盒v1.1 | 10x Genomics | PN-1000165 |
| 空间转录组 | Visium空间组织优化试剂盒 | 10x Genomics | PN-1000193 |
| Visium空间基因表达试剂盒 | 10x Genomics | PN-1000187 | |
| 多重免疫组化 | Opal七色手动免疫组化试剂盒 | AKOYA BIOSCIENCES | NEL81 1001KT |
| 内源性过氧化物酶封闭缓冲液 | Beyotime | P0100A | |
| CD14抗体 | Proteintech | 60253-1-Ig | |
| C1qC抗体 | Proteintech | 16889-1-AP | |
| IGF1抗体 | Proteintech | 28530-1-AP | |
| SPP1抗体 | Proteintech | 22952-1-AP | |
| PDPN抗体 | Proteintech | 11629-1-AP | |
| ANXA1抗体 | Abcam | ab214486 | |
| CXCL13抗体 | Abcam | ab246518 | |
| PD-1抗体 | Abcam | ab137132 | |
| CD4抗体 | Abcam | ab133616 | |
| FOXP3抗体 | Abcam | ab20034 | |
| IL-17a抗体 | Proteintech | 26163-1-AP | |
| HLA-DR抗体 | Abcam | ab92511 | |
| CD20抗体 | Abcam | ab78237 | |
| CLEC9A抗体 | Proteintech | 55451-1-AP | |
| CD1c抗体 | Abcam | ab246520 | |
| CD208抗体 | Abcam | ab271053 | |
| GZMK抗体 | Abcam | ab282703 | |
| CD8抗体 | Abcam | ab17147 | |
| CD103抗体 | Abcam | ab224202 | |
| CTHRC1抗体 | Abcam | ab256458 | |
| TCF1抗体 | Cell Signaling Technology | 2203 T | |
| 聚合物HRP小鼠+兔二抗 | -- | ARH1001EA | |
| DAPI核染液 | Invitrogen | D1306 | |
| 多重免疫组化成像 | Vectra Polaris成像系统 | PerkinElmer | -- |
| 小鼠处理 | 罗格列酮 | Selleck | S2556 |
| 抗IL11RA抗体 | Invitrogen | PA5-47667 | |
| 小鼠流式分析 | Zombie Aqua可固定活力检测试剂盒 | BioLegend | 423102 |
| Fc受体封闭剂 | BioXCell | BE0307 | |
| CD45-BUV395抗体 | BD Biosciences | 564279 | |
| CD31-FITC抗体 | Invitrogen | 11-0311-81 | |
| CD29-PE抗体 | Invitrogen | 12-0291-81 | |
| Sca-1-APC/Cy7抗体 | BioLegend | 108125 | |
| CD34-Pacific Blue抗体 | BioLegend | 119321 | |
| CD55-PE/Cy7抗体 | BioLegend | 131813 | |
| CD90.2-BV711抗体 | BioLegend | 105349 | |
| 小鼠多重免疫组化 | CD11b抗体 | Abcam | ab133357 |
| IL-1β抗体 | Invitrogen | P420B | |
| CD163抗体 | Abcam | ab182422 | |
| CD68抗体 | Proteintech | 28058-1-AP | |
| CD74抗体 | Abcam | ab289885 | |
| GPNMB抗体 | Proteintech | 66926-1-Ig | |
| α-SMA抗体 | Abcam | ab124964 | |
| PDGFRα抗体 | Abcam | ab96569 | |
| CD34抗体 | Invitrogen | MA5-29674 | |
| OGN抗体 | Proteintech | 12755-1-AP | |
| SOX9抗体 | CST | 82630T | |
| II型胶原抗体 | Proteintech | 28459-1-AP | |
| IGFBP5抗体 | Proteintech | 55205-1-AP | |
| PRG4抗体 | Proteintech | 25732-1-AP | |
| K14抗体 | Santa Cruz | sc-53253 | |
| CD31抗体 | Abcam | ab281583 | |
| FcεRIα抗体 | Proteintech | 10980-1-AP | |
| Ly6G抗体 | Abclonal | A22270 | |
| IL11RA抗体 | Proteintech | 10264-1-AP | |
| 骨桥蛋白抗体 | Proteintech | 22952-1-AP | |
| 小鼠单细胞分选 | 钙黄绿素 | Invitrogen | C3099 |
| 小鼠单细胞文库构建 | Chromium Next GEM单细胞3'GEM文库与凝胶珠试剂盒V3 | 10x Genomics | PN-1000075 |
关键环节
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 模拟ST数据基准测试 | 模拟spot数量、每spot细胞数、参考亚群数、复杂度分组方式及统计检验 阅读提示：Methods中Simulation strategy与Evaluation metrics小节，Fig. 2A–G图注 |
| 人CRC空间转录组与mIHC验证 | 发现队列患者/切片数、验证队列样本、scRNA-seq参考亚群数、mIHC抗体组合及定量阈值 阅读提示：Methods中Human specimens与Multi-color immunohistochemistry小节，Fig. 3图注 |
| 克隆免疫微环境分型 | 克隆数、T细胞富集/排斥分组标准、Ro/e阈值、群落聚类方法 阅读提示：Methods中Tumor clonal inference、ClonalTIME phenotyping与Multi-cellular community analysis小节，Fig. 3H–N图注 |
| chr20q-gain与HERV-H分析 | CNV判定阈值、TCGA-COAD样本量与筛选、HERV-H位点定义、MSK-IMPACT患者数与ICB分组 阅读提示：Methods中TCGA-COAD data analysis与MSK-IMPACT data analysis小节，Fig. 4图注 |
| 小鼠耳伤口愈合ST与scRNA-seq | 小鼠品系与周龄、耳穿孔尺寸、取样时间点、scRNA-seq参考亚群数、ST渗透时间 阅读提示：Methods中Ear wound creation and treatment、Spatial gene expression assay小节，Fig. 5图注 |
| 促纤维化群落功能干预 | rosiglitazone剂量与给药频率、anti-IL11RA剂量与给药时间、对照组设置、每组样本量及终点指标 阅读提示：Methods中Ear wound creation and treatment小节，Fig. 6L–M与Supplementary Fig. S7S–W图注 |