scINSIGHT:用于解释生物学异质性数据中单细胞基因表达的方法

scINSIGHT for interpreting single-cell gene expression from biologically heterogeneous data

作者信息Kun Qian, Shiwei Fu, Hongwei Li, Wei Vivian Li
PMID35313930
发布时间2022-03-21
DOI10.1186/s13059-022-02649-3

实验完整度

高

包含模拟数据基准测试、三个真实数据应用(黑色素瘤、COVID-19、皮肤创面)以及计算性能评估,并进行了与多种方法的比较验证。

主要模型

模拟单细胞数据(6个样本,3个时间点) 黑色素瘤患者的CD8+ T细胞(48个肿瘤活检样本) COVID-19患者外周血B细胞(14个血液样本) 小鼠皮肤创面真皮细胞(对照与Hh激活条件)

重点核对

模拟数据中,三个细胞类型(C1, C2, C3)在所有样本中存在,而C4, C5, C6仅存在于特定条件;每个样本每种细胞类型100个细胞 黑色素瘤数据中,48个肿瘤样本来自32例患者,分为31例无应答(NR)和17例应答(R) COVID-19数据中,14个血液样本分为5例健康、4例复杂、5例恢复/出院前 皮肤创面数据中,两个样本(对照和Hh激活) scINSIGHT参数:K(公共模块数)在模拟中为13,黑色素瘤为9,COVID-19为13,创面为11;λ1=λ2在模拟中为0.01,黑色素瘤和COVID-19为1,创面为0.1

摘要

单细胞RNA测序(scRNA-seq)数据的不断增加强调了对整合分析的需求,以解释单细胞样本之间的相似性和差异性。尽管已经开发了不同的批次效应去除方法,但没有一种方法适用于来自多个生物学条件的异质性单细胞样本。我们提出了一种名为scINSIGHT的方法,用于学习在不同生物学条件下共同或特异的协调基因表达模式,并识别跨单细胞样本的细胞身份和过程。我们将scINSIGHT与最新方法在模拟和真实数据上进行了比较,结果证明了其性能的改进。我们的结果表明scINSIGHT在多种生物医学和临床问题中的适用性。

实验结论

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

研究问题
如何从多个生物学异质性来源的单细胞样本中整合分析,以识别共同的细胞身份和条件特异性基因表达模式?
核心机制
scINSIGHT通过非负矩阵分解,将每个样本的基因表达分解为条件特异性基因模块和公共基因模块,从而捕捉共同和条件特异的生物学过程。
主要证据
模拟数据中,scINSIGHT实现了最高的ARI(0.99)和可靠的整合评分;在黑色素瘤数据中,scINSIGHT识别的6个聚类中有4个与免疫治疗应答显著相关(P<0.01),并通过公共模块的表达和通路富集分析揭示了T细胞状态;在COVID-19数据中,scINSIGHT识别的B细胞聚类与参考注释的ARI最高;在创面数据中,scINSIGHT实现了最高的Silhouette和整合评分,并揭示了Hh激活成纤维细胞的异质性。
研究意义
scINSIGHT能够整合来自不同生物学条件的单细胞数据,同时保留条件特异性信息,适用于广泛的生物医学和临床问题。

研究路径

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

1

模拟数据生成与基准测试

使用已知的细胞类型组成和条件特异性效应,评估scINSIGHT的聚类性能和整合能力。

使用scDesign模拟了六个单细胞样本,属于三个时间点(T1, T2, T3),每个时间点两个样本。包含六种细胞类型,其中C1-C3在所有样本中,C4-C6仅在特定条件中出现。将scINSIGHT应用于模拟数据,并与六种先进方法(Seurat、LIGER、Harmony、scMerge、scGen、Scanorama)进行比较。

2

集成CD8+ T细胞数据以识别黑色素瘤免疫治疗应答相关状态

验证scINSIGHT在真实数据中整合来自不同患者和条件的样本,并识别与免疫治疗应答相关的T细胞状态。

应用scINSIGHT到48个黑色素瘤肿瘤活检样本的CD8+ T细胞单细胞数据(6350个细胞),将NR/R作为条件因子,识别细胞聚类和公共/条件特异性基因模块。与Seurat、LIGER、Harmony、scGen、Scanorama比较。

3

整合COVID-19 B细胞数据以分析疾病阶段相关细胞类型

评估scINSIGHT在复杂临床数据中识别与疾病阶段相关的B细胞亚型和基因模块的能力。

应用scINSIGHT到14个外周血样本的9741个B细胞数据,将疾病阶段(健康、复杂、恢复/出院前)作为条件因子,识别公共和阶段特异性基因模块,并使用SingleR注释细胞类型。与六种方法比较。

4

整合小鼠皮肤创面真皮细胞数据以分析Hh激活效应

测试scINSIGHT在仅有两个样本(每个条件一个)时,识别细胞类型和公共模块的能力。

应用scINSIGHT到对照和Hh激活条件下的鼠真皮细胞单细胞数据,将条件因子设为对照/治疗,识别公共和条件特异性模块,并注释细胞类型。与六种方法比较。

5

评估运行时间和内存使用

评估scINSIGHT的计算性能,并与现有方法比较。

记录了scINSIGHT和其他六种方法在多个数据集上的运行时间和内存使用,包括模拟数据、黑色素瘤、COVID-19、创面和小鼠视网膜数据集。

研究方法

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

产品清单

实验环节名称品牌货号
scDesign----
Seurat--v3.2.3
LIGER--v0.5.0
Harmony--v0.1.0
scMerge--v1.6.0
scGen--v2.0.0
Scanorama--v1.7.1
SingleR----
cluster----
mclust----

关键环节

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

环节核对要点
模拟数据生成
模拟样本数量、细胞类型组成、每个细胞类型的细胞数、条件效应权重(0.9和0.1)
阅读提示:Methods: Simulation and analysis of synthetic data
黑色素瘤数据整合
基因特征选择(CD8+基因签名或2000个高度可变基因)、条件因子(NR/R)、参数K和λ1/λ2值(K=9,λ1=λ2=1)
阅读提示:Methods: Analysis of real data(黑色素瘤部分)
COVID-19数据整合
基因特征选择(2000个高度可变基因)、条件因子(疾病阶段)、参数K和λ1/λ2(K=13,λ1=λ2=1)
阅读提示:Methods: Analysis of real data(COVID-19部分)
皮肤创面数据整合
条件因子(对照/治疗)、参数K和λ1/λ2(K=11,λ1=λ2=0.1)
阅读提示:Methods: Analysis of real data(创面部分)
计算性能评估
运行时间、内存使用、数据集大小
阅读提示:Results: Computational time and memory usage