推导phiclust的理论基础
从随机矩阵理论和低秩扰动理论出发,定义phiclust作为聚类可行性的度量。
将测量基因表达矩阵分解为随机噪声矩阵和信号矩阵,利用奇异值分解,通过低秩扰动理论计算信号矩阵的奇异值,进而得出phiclust。
Phiclust: a clusterability measure for single-cell transcriptomics reveals phenotypic subpopulations
包含模拟数据、实验数据(PBMC、RNA混合、BMNC、胎儿肾脏)等多个证据层级,并进行了功能验证(如免疫染色)和机制验证(如g-phiclust基因分析)。
PBMC数据 RNA混合数据(3种细胞系) 骨髓单核细胞(BMNC) 胎儿人类肾脏组织
数据预处理:归一化、对数转换、基因-wise和细胞-wise标准化 混淆因子回归:总转录本数、核糖体基因表达、线粒体基因表达、细胞周期评分、应激相关基因 通过UMAP进行聚类可视化,检查聚类结构与phiclust值的对应关系 通过免疫染色验证胎儿肾脏样本中亚群的体内存在 通过差异表达分析验证基因-奇异向量识别的驱动基因
提炼研究问题、关键发现与证据,快速把握文章的核心贡献。
按研究推进顺序梳理实验设计、验证步骤与关键观察。
从随机矩阵理论和低秩扰动理论出发,定义phiclust作为聚类可行性的度量。
将测量基因表达矩阵分解为随机噪声矩阵和信号矩阵,利用奇异值分解,通过低秩扰动理论计算信号矩阵的奇异值,进而得出phiclust。
验证phiclust是否能作为聚类可行性(以tARI衡量)的代理指标。
使用Splatter模拟双簇数据,改变差异表达基因数量或log fold change,计算phiclust和tARI,并比较两者相关性。
在真实实验数据上验证phiclust与tARI的相关性,并控制真实噪声结构。
从PBMC数据中提取B细胞和CD8 T细胞,通过加权平均构建合成簇,改变权重α,比较phiclust和tARI。
测试phiclust在实验数据中检测子结构的能力,并与其他方法(ROGUE,轮廓系数)比较。
使用RNA混合数据集,计算phiclust,并应用混淆因子回归去除输入RNA量和总计数的影响。
展示phiclust在复杂组织(骨髓单核细胞)中优先选择有意义的簇进行亚聚类。
计算BMNC数据中各个簇的phiclust,对高phiclust的簇(红细胞祖细胞、树突状细胞祖细胞)进行亚聚类,并识别差异基因。
验证phiclust在复杂组织(胎儿肾脏)中发现的亚群是否具有生物学意义。
分析胎儿肾脏数据,对高phiclust的簇(UBCD、SSBpr、ICa)进行亚聚类,并对新发现的亚群进行免疫荧光染色验证。
按研究目的归类文中使用的方法,便于定位所需技术。
| 实验环节 | 名称 | 品牌 | 货号 |
|---|---|---|---|
| 免疫荧光染色 | 抗UPK1A抗体 | Atlas Antibodies | HPA049879 |
| 抗KRT7抗体 | Thermo Fisher Scientific | MA5-11986 | |
| 抗CDH1抗体 | Santa Cruz | SC-7870 | |
| 抗CLDN1抗体 | Thermo Fisher Scientific | 717800 | |
| 抗CAV2抗体 | R&D Systems | AF5788-SP | |
| 抗AKAP12抗体 | Santa Cruz | sc-376740 | |
| 抗CLDN11抗体 | SIGMA Aldrich | HPA013166 | |
| 抗POSTN抗体 | Santa Cruz | sc-398631 | |
| 抗SULT1E1抗体 | R&D Systems | AF5545-SP | |
| Alexa Fluor 594 驴抗小鼠二抗 | Invitrogen | A21203 | |
| Alexa Fluor 594 驴抗兔二抗 | Invitrogen | A21207 | |
| Alexa Fluor 647 驴抗小鼠二抗 | Invitrogen | A31571 | |
| Alexa Fluor 647 驴抗兔二抗 | Invitrogen | A31573 | |
| Alexa Fluor 647 驴抗山羊二抗 | Invitrogen | A21447 | |
| 成像 | 尼康Ti-Eclipse落射荧光显微镜 | Nikon | -- |
| Andor iXON Ultra 888 EMCCD相机 | Andor | -- |
汇总复现实验时建议重点确认的条件及原文阅读提示。
| 环节 | 核对要点 |
|---|---|
| 数据预处理 | 归一化和log转换,以及基因-wise和细胞-wise标准化 阅读提示:Methods - Preprocessing |
| phiclust计算 | 奇异值分解,并计算phiclust值 阅读提示:Methods - Phiclust |
| 混淆因子回归 | 用于回归的混淆因子列表,如总转录本数、核糖体基因表达、线粒体基因表达、细胞周期评分、应激相关基因 阅读提示:Methods - Confounder regression |
| 免疫荧光染色 | 使用的抗体及其稀释度、供应商和货号 阅读提示:Methods - Staining |