单细胞分析:聚类流程(六)
1. 学习目标
描述质量控制步骤后scRNA-seq
分析的工作流程(聚类)。
现在有了高质量的细胞,可以继续工作流程。最终,希望对细胞进行聚类并识别不同的潜在细胞类型,但是在那之前需要完成几个步骤。下面的工作流程示意图中的绿色框对应于QC
后采取的步骤,共同构成了聚类工作流程。
2. 聚类流程
对于具有信息性的事物,它需要表现出变化,但并非所有变化都是信息性的。聚类分析的目的是保留数据集中应该定义的细胞类型的主要变异来源,同时消减由于无意义的变异来源(测序深度、细胞周期差异、线粒体表达、批次效应等)引起的变异。然后,为了确定存在的细胞类型,将使用特异的基因进行聚类分析,以定义数据集的主要来源。
要识别clusters
,将执行以下步骤:
2.1. 变异溯源
探索 unwanted variation
的来源
工作流程的第一步是查看数据是否包含任何不需要的变异。在scRNA-seq
数据中评估的最常见的生物学效应是细胞周期对转录组的影响。另一种已知的生物学效应是线粒体基因表达,它被解释为细胞应激的一种迹象。工作流程的这一步涉及探索数据以确定想要回归的协变量。
2.2. 回归
归一化和回归 unwanted variation
来源
Seurat
最近推出了一种名为 sctransform
的新方法,该方法对 scRNA-seq
数据执行多个处理步骤。需要标准化以缩放原始计数数据以获得细胞之间正确的相对基因表达丰度。sctransform
函数实现了数据的高级归一化和方差稳定。 sctransform
函数还回归了数据中不需要的变化的来源。在上一步中,已经确定了这些可变性来源,在这里指定了这些协变量是什么。
2.3. 整合
整合(可选)
通常使用 scRNA-seq
,会处理对应于不同样本组、多个实验或不同模式的多个样本。如果想最终比较组之间的细胞类型表达,建议整合数据。整合是一种强大的方法,它使用这些最大变化的共享源来识别跨条件或数据集的共享子群。在Seurat
中执行整合涉及几个步骤。完成后,使用可视化方法来确保良好的整合,然后再对单元进行聚类。
2.4. 细胞聚类
细胞聚类
通过基于其基因表达谱的相似性对细胞进行分组来获得细胞簇。表达轮廓相似性是通过距离度量来确定的,它通常以降维表示作为输入。Seurat
根据从整合的最多变基因的表达中获得的 PCA
分数将细胞分配到簇中。
2.5. 质量评估
聚类质量评估
数据中确定的clusters
代表可能属于相似细胞类型的细胞组。在确认一组成员单元格的单元格类型之前,需要执行以下步骤:
检查 clusters
是否不受无意义变异来源的影响。
检查主要的主成分是否在驱动不同的 clusters
。
通过查看 clusters
中已知标记的表达式来探索细胞类型身份。
本文由 mdnice 多平台发布
单细胞分析:聚类流程(六)相关推荐
- 生信入门(六)——单细胞分析(Seurat)
生信入门(六)--单细胞分析(Seurat) 文章目录 生信入门(六)--单细胞分析(Seurat) 一.数据导入 1.数据来源 2.数据导入 二.标准预处理 1.QC和选择细胞进行进一步分析 2.规 ...
- 中科院单细胞分析算法开发博士带你做单细胞转录组分析
" 福利公告:为了响应学员的学习需求,经过易生信培训团队的讨论筹备,现决定安排扩增子16S分析.宏基因组和Python课程的线上直播课.报名参加线上直播课的老师可在1年内选择参加同课程的一次 ...
- 如何使用Bioconductor进行单细胞分析?
最近的技术进步使得能够在单个细胞中分析全基因组特征.但是,单细胞数据为分析提出了独特的挑战,需要开发专用的方法和数据架构才能成功解析数据背后的生物问题.Bioconductor项目托管了社区开发的开源 ...
- 复现Cell附图 |类器官的单细胞分析
类器官的单细胞分析 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单细胞测序分析 (重 ...
- 单细胞分析Seurat使用相关的10个问题答疑精选!
NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单细胞测序分析 (重磅综述:三万字长文读 ...
- 收藏|北大生信平台单细胞分析、染色质分析视频和PPT分享
生物信息学习的正确姿势 NGS系列文章包括NGS基础.在线绘图.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单细胞 ...
- 单细胞cellranger处理流程
单细胞cellranger处理流程 单细胞cellranger处理代码 aggr.csv的制作 单细胞cellranger处理代码 一套流程直接跑通! 首先你要cd切换到工作目录,需要更改你的refe ...
- 单细胞分析可视化工具盘点
论文阅读 单细胞分析可视化工具盘点 首先:大概知道单细胞分析的十几种格式和可视化分析工具. 第二:输入格式中(下图表绿色标记):csv/txt格式是最常被接受的格式,有8个工具支持.更专业的格式,如h ...
- 熟练掌握R语言的Meta分析全流程和不确定性分析,并结合机器学习等方法讲解Meta分析在文献大数据的延伸应用
Meta分析是针对某一科研问题,根据明确的搜索策略.选择筛选文献标准.采用严格的评价方法,对来源不同的研究成果进行收集.合并及定量统计分析的方法,最早出现于"循证医学",现已广泛应 ...
最新文章
- redis主从复制、高可用和集群
- 互联网+连接教育 全息时代即将来临
- python输出去空格_Python3基础 print(,end=) 输出内容的末尾加入空格
- Netflix推荐系统(Part two)-系统架构
- Android----获取activity上所有的控件
- 项目实践中Linux集群的总结和思考
- 网站渗透测试原理及详细过程
- linux共享文件系统sy,Linux使用Samba实现文件共享
- iOS开发--XMPPFramework--环境的配置(一)
- 深度强化学习算法的未来——样本效率研究
- videoview实现视频引导页及从assets文件中读取文件到本地SD卡
- Redis 使用 scan 命令代替 keys
- Steam版XCOM: Enemy Within(内部敌人)不能启动的问题
- Sencha学习笔记4: Creating your First App - 官方创建您的第一个Sencha Touch应用指导
- 苹果手表Apple Watch 无法检查更新 iOS14.5 14.6 14.7
- 牛牛的猜球游戏(前缀和+逆交换)
- PyQt(Python+Qt)学习随笔:Qt Designer中部件的调色板palette属性和字体font属性设置
- supervisor 进程托管程序的使用示例
- 利用腾讯公有云建设灾备中心的解决方案
- RationalDMIS 7.1 程序示例