RunPrestoAll:Seurat差异基因分析加速利器!

\,\,\,\,\,\, 单细胞转录组测序的差异基因表达分析是一种用于比较不同细胞之间基因表达水平的方法。在单细胞RNA测序数据中,每个细胞的基因表达水平都被测量,并以基因计数矩阵的形式进行存储。差异基因表达分析旨在识别不同细胞群体之间的基因表达水平差异,从而揭示细胞群体之间的功能和特征。Seurat 提供了一个 FindAllMarkers 方法用于在单细胞RNA测序数据中寻找差异表达基因。然而,对于大型数据集的DE分析,使用Seurat软件包的FindAllMarkers方法 在数据集的全部细胞上执行DE搜索将变得非常缓慢。即使开启了future并行,速度也并不会显著提升,因为多线程模式只能在某些特定的计算环节中发挥作用,而在其他计算环节中仍然需要使用单线程模式。在加速搜索差异基因, FindAllMarkers 函数本身提供了对基因或细胞进行预过滤的参数,特别是设置 max.cells.per.ident 对每个类别进行降采样,使其中用于比较的细胞数量不超过设置的阈值。虽然通常会出现一定的性能损失,但速度增加会比较显,并且差异表达最高的基因可能仍会排列在最前面。

如果计算平台允许,并不推荐使用 FindAllMarkers(max.cells.per.ident = *) 通过降采样加速来执行DEG 搜索,而是推荐基于 Rcpp 实现的 SeuratWrappers::RunPrestoAll 方法,具有与 Seurat::FindAllMarkers 一样的功能,并提供更迅速的搜索。

1、FindAllMarkers 降采样DEG计算性能损失评估

降采样差异基因搜索代码示例:

还有 79% 的精彩内容
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。
禁止转载,如需转载请通过简信或评论联系作者。
支付 ¥1.00 继续阅读

相关阅读更多精彩内容

友情链接更多精彩内容