单细胞分辨率下细胞间通讯的比较分析

Comparative analysis of cell-cell communication at single-cell resolution.

作者信息Aaron J Wilk, Alex K Shalek, Susan Holmes, Catherine A Blish
PMID37169965
发布时间2024-03
DOI10.1038/s41587-023-01782-z

实验完整度

高

含公共图谱数据集、CRISPRa扰动筛选、NK-B细胞转染共培养scRNA-seq及空间转录组多重证据与功能验证

主要模型

已发表scRNA-seq图谱数据集(SCC、麻风肉芽肿、胎肠、SARS-CoV-2感染HBEC) CRISPRa Perturb-seq激活人T细胞 转染CD40L/CD40的NK-B细胞共培养 空间转录组数据集

重点核对

NK细胞转染CD40LG mRNA、B细胞转染CD40 mRNA,设GFP-GFP、GFP-CD40、CD40L-GFP、CD40L-CD40四组,样本量4934/5665/4908/6031 NK与B细胞按1:1共培养12 h后进行流式与scRNA-seq CRISPRa Perturb-seq使用靶向15种细胞表面配体/受体的gRNA,AUC约0.93 空间转录组11个数据集比较top 1%高互作细胞对与随机置换距离 麻风肉芽肿数据集含5例LL与4例RR患者

摘要

从单细胞RNA测序数据推断细胞间通讯是揭示细胞间通讯通路的有力技术,然而现有方法在细胞类型或聚类水平进行此类分析,丢弃了单细胞水平的信息。在此,我们提出Scriabin,一个灵活且可扩展的框架,用于在单细胞分辨率下进行细胞间通讯的比较分析,无需细胞聚集或降采样。我们使用多个已发表的大规模图谱数据集、遗传扰动筛选和直接实验验证,证明Scriabin能够准确恢复预期的细胞间通讯边,并识别可能被凝聚方法掩盖的通讯网络。此外,我们使用空间转录组数据证明,Scriabin仅从解离数据即可揭示相互作用的空间特征。最后,我们展示了对纵向数据集的应用,以追踪时间点之间运行的通讯通路。我们的方法代表了一种广泛适用的策略,可揭示健康和疾病中生态位-表型关系的完整结构。

实验结论

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

研究问题
能否在单细胞分辨率下进行细胞间通讯的比较分析,避免细胞聚集或降采样丢失信息,并准确识别具有生物学意义的通讯边
核心机制
Scriabin结合配体-受体数据库、下游信号模型、基于锚定的数据整合(binning)与基因网络分析,对每对细胞的配体-受体表达取几何均值构建CCIM,并用NicheNet配体活性加权
主要证据
SCC数据中识别耗竭T细胞与CD1C+ DC通讯的CTLA4/TIGIT及CCL4/CCL5丢失;NK-B转染实验仅CD40LG-CD40边显著改变;空间转录组高互作细胞对距离显著更近;麻风肉芽肿中IFNG为最重要配体
研究意义
提供广泛适用的策略以揭示健康和疾病中生态位-表型关系的完整结构,发现被凝聚方法掩盖的通讯通路并支持纵向通讯电路分析

研究路径

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

1

Scriabin框架与CCIM构建

在单细胞分辨率推断细胞间通讯边

对每对细胞的每个配体-受体对计算表达几何均值构建CCIM,并以NicheNet配体活性加权

2

技术噪声与鲁棒性评估

评估稀疏噪声对单细胞通讯分析的影响

将高深度数据集随机降采样至inDrop覆盖度,比较原始、ALRA去噪、邻域聚集与聚类伪bulk的CCIM对GT的再现

3

空间转录组验证

验证Scriabin预测的高互作细胞对是否空间邻近

在11个空间转录组数据集上运行Scriabin,限制为膜结合配体-受体,比较top 1%互作细胞对距离与随机置换

4

CRISPRa扰动筛选验证

验证能否从通讯谱预测gRNA

分析Schmidt等激活人T细胞CRISPRa Perturb-seq数据,基于CCIM预测gRNA并计算ROC

5

NK-B转染共培养直接验证

提供Scriabin检测通讯变化的直接实验证据

NK细胞转染CD40LG mRNA、B细胞转染CD40 mRNA,四组共培养后scRNA-seq检测CCIM差异

6

麻风肉芽肿比较分析

检验Scriabin能否避免凝聚方法掩盖的通讯边

对9例肉芽肿数据binning后比较LL与RR差异互作,识别CRTAM+ T细胞IFNG与髓系细胞应答

7

胎肠相互作用程序发现

展示可扩展的共表达配体-受体模块发现

对76,592细胞胎肠图谱迭代近似TOM发现75个相互作用程序并评分

8

纵向通讯电路组装

追踪跨时间点通讯通路

对SARS-CoV-2感染HBEC每日采样3天数据binning并组装纵向通讯电路

研究方法

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

产品清单

实验环节名称品牌货号
Ficoll-Paque密度梯度分离液GE Healthcare--
苯甲酸酶EMD Millipore--
RPMI 1640培养基----
NK细胞磁珠分选试剂盒Miltenyi Biotec130-092-657
B细胞磁珠分选试剂盒Miltenyi Biotec130-101-638
电荷改变可释放转运体BDK-O7:N7:A13----
编码eGFP的mRNATriLink BioTechnologies--
编码CD40的mRNATriLink BioTechnologies--
编码CD40L的mRNATriLink BioTechnologies--
可固定活性染料eFluor 780Thermo Fisher Scientific--
Aurora流式细胞仪Cytek Biosciences--
mRNA捕获磁珠ChemGenes--
Maxima H逆转录酶Thermo Fisher ScientificEPO0753
核酸外切酶INew England Biolabs (NEB)M0293L
Klenow片段NEBM0212L
KAPA HiFi PCR预混液Kapa BiosystemsKK2602
Agencourt AMPure XP磁珠Beckman CoulterA63881
Nextera XT DNA文库制备试剂盒IlluminaFC-131-1096
NovaSeq 6000测序仪Illumina--
聚碳酸酯膜----

关键环节

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

环节核对要点
NK-B转染共培养
四组转染条件及各自样本量:GFP-GFP 4934、GFP-CD40 5665、CD40L-GFP 4908、CD40L-CD40 6031
阅读提示:Methods中Transfection and co-culture of primary NK and B cells;Fig. 3d
转染条件
CART/mRNA电荷比10:1,血清-free孵育6 h,共培养12 h
阅读提示:Methods中Transfection and co-culture of primary NK and B cells
scRNA-seq质量控制
过滤标准:UMI<750或>15000、线粒体/rRNA>20%、每100 UMI表达>75基因
阅读提示:Methods中Alignment and quality control of scRNA-seq data
麻风肉芽肿比较
5例LL与4例RR分组,bin比较采用t统计量
阅读提示:Fig. 4 legend及Results相关段落
空间转录组验证
11个独立数据集,限制膜结合配体-受体,top 1%互作距离与随机置换比较
阅读提示:Methods中Analysis of spatial transcriptomic datasets with Scriabin;Fig. 3a,b
CRISPRa扰动验证
靶向15种细胞表面配体/受体的gRNA,使用作者gRNA注释
阅读提示:Methods中Analysis of CRISPRa Perturb-seq data;Fig. 3c