Phiclust:一种用于单细胞转录组学的聚类可行性度量,揭示表型亚群

Phiclust: a clusterability measure for single-cell transcriptomics reveals phenotypic subpopulations

作者信息Maria Mircea, Mazène Hochane, Xueying Fan, Susana M Chuva de Sousa Lopes, Diego Garlaschelli, Stefan Semrau
PMID35012604
发布时间2022-01-10
DOI10.1186/s13059-021-02590-x

实验完整度

包含模拟数据、实验数据(PBMC、RNA混合、BMNC、胎儿肾脏)等多个证据层级,并进行了功能验证(如免疫染色)和机制验证(如g-phiclust基因分析)。

主要模型

PBMC数据 RNA混合数据(3种细胞系) 骨髓单核细胞(BMNC) 胎儿人类肾脏组织

重点核对

数据预处理:归一化、对数转换、基因-wise和细胞-wise标准化 混淆因子回归:总转录本数、核糖体基因表达、线粒体基因表达、细胞周期评分、应激相关基因 通过UMAP进行聚类可视化,检查聚类结构与phiclust值的对应关系 通过免疫染色验证胎儿肾脏样本中亚群的体内存在 通过差异表达分析验证基因-奇异向量识别的驱动基因

摘要

通过单细胞转录组的无监督聚类发现新的细胞表型的能力已经彻底改变了生物学。目前,还没有一种原则性的方法来判断一个细胞簇是否包含应进一步解析的有意义的亚群。在这里,我们提出了phiclust(ϕclust),一种源自随机矩阵理论的聚类可行性度量,可用于识别具有非随机子结构的细胞簇,并可测试性地导致发现以前被忽视的表型。

实验结论

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

研究问题
如何确定单细胞RNA测序数据中的一个细胞簇是否包含有意义的亚群,而不仅仅是随机噪声?
核心机制
phiclust基于随机矩阵理论和低秩扰动理论,通过衡量测量表达矩阵与无噪声信号矩阵之间奇异向量的夹角来量化聚类可行性,且不受噪声分布类型限制。
主要证据
模拟数据、PBMC数据、RNA混合数据以及BMNC和胎儿肾脏数据集上的验证显示,phiclust与理论可达的调整兰德指数(tARI)强相关,并能识别出通过后续实验验证的亚群(如胎儿肾脏中的尿上皮细胞、壁层上皮细胞和间质细胞亚型)。
研究意义
phiclust提供了一种无参数、基于第一性原理的聚类可行性度量,有助于区分有意义的亚群与随机噪声,提高单细胞RNA-seq分析的可重复性和稳健性,并可指导后续实验验证。

研究路径

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

1

推导phiclust的理论基础

从随机矩阵理论和低秩扰动理论出发,定义phiclust作为聚类可行性的度量。

将测量基因表达矩阵分解为随机噪声矩阵和信号矩阵,利用奇异值分解,通过低秩扰动理论计算信号矩阵的奇异值,进而得出phiclust。

2

在模拟数据上验证phiclust与tARI的相关性

验证phiclust是否能作为聚类可行性(以tARI衡量)的代理指标。

使用Splatter模拟双簇数据,改变差异表达基因数量或log fold change,计算phiclust和tARI,并比较两者相关性。

3

在PBMC数据上验证phiclust与tARI的相关性

在真实实验数据上验证phiclust与tARI的相关性,并控制真实噪声结构。

从PBMC数据中提取B细胞和CD8 T细胞,通过加权平均构建合成簇,改变权重α,比较phiclust和tARI。

4

在RNA混合数据上评估phiclust的性能

测试phiclust在实验数据中检测子结构的能力,并与其他方法(ROGUE,轮廓系数)比较。

使用RNA混合数据集,计算phiclust,并应用混淆因子回归去除输入RNA量和总计数的影响。

5

在BMNC数据上应用phiclust发现亚群

展示phiclust在复杂组织(骨髓单核细胞)中优先选择有意义的簇进行亚聚类。

计算BMNC数据中各个簇的phiclust,对高phiclust的簇(红细胞祖细胞、树突状细胞祖细胞)进行亚聚类,并识别差异基因。

6

在胎儿肾脏数据上应用phiclust并实验验证新亚群

验证phiclust在复杂组织(胎儿肾脏)中发现的亚群是否具有生物学意义。

分析胎儿肾脏数据,对高phiclust的簇(UBCD、SSBpr、ICa)进行亚聚类,并对新发现的亚群进行免疫荧光染色验证。

研究方法

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

产品清单

实验环节名称品牌货号
抗UPK1A抗体Atlas AntibodiesHPA049879
抗KRT7抗体Thermo Fisher ScientificMA5-11986
抗CDH1抗体Santa CruzSC-7870
抗CLDN1抗体Thermo Fisher Scientific717800
抗CAV2抗体R&D SystemsAF5788-SP
抗AKAP12抗体Santa Cruzsc-376740
抗CLDN11抗体SIGMA AldrichHPA013166
抗POSTN抗体Santa Cruzsc-398631
抗SULT1E1抗体R&D SystemsAF5545-SP
Alexa Fluor 594 驴抗小鼠二抗InvitrogenA21203
Alexa Fluor 594 驴抗兔二抗InvitrogenA21207
Alexa Fluor 647 驴抗小鼠二抗InvitrogenA31571
Alexa Fluor 647 驴抗兔二抗InvitrogenA31573
Alexa Fluor 647 驴抗山羊二抗InvitrogenA21447
尼康Ti-Eclipse落射荧光显微镜Nikon--
Andor iXON Ultra 888 EMCCD相机Andor--

关键环节

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

环节核对要点
数据预处理
归一化和log转换,以及基因-wise和细胞-wise标准化
阅读提示:Methods - Preprocessing
phiclust计算
奇异值分解,并计算phiclust值
阅读提示:Methods - Phiclust
混淆因子回归
用于回归的混淆因子列表,如总转录本数、核糖体基因表达、线粒体基因表达、细胞周期评分、应激相关基因
阅读提示:Methods - Confounder regression
免疫荧光染色
使用的抗体及其稀释度、供应商和货号
阅读提示:Methods - Staining