Seurat V5 | 最火的单细胞测序分析工具重磅更新啦!~

1写在前面

现在最炙手可热的单细胞分析包,Seurat重磅跟新啦!~🤩

Seurat最初是由纽约大学的Rafael A. IrizarrySatija等人于2015年开发。😘

该工具基于R语言编写,使用了许多先进的统计学和机器学习算法,可以对scRNA-seq数据进行细胞聚类、细胞亚群鉴定、基因差异表达分析等。🥰

自推出以来,Seurat不断更新和改进,增加了许多新功能,包括对空间转录组数据的分析和可视化,以及对多组学数据的整合。

每次发文都是顶刊,也是让人着实羡慕啊!~🤩(左右滑动哦,亲)

alt
alt
alt
alt
alt

本次更新至V5也是做了很多新的升级:👇

  • 空间转录组分析和多组学分析;

  • 分析流程进一步简化和集成;

  • 处理超大数据;

  • 数据结构优化;

alt

2如何安装Seurat V5

Seurat V5暂时还不能通过CRAN安装,通过下面这段代码安装吧。🫶

remotes::install_github("satijalab/seurat", "seurat5", quiet = T)
library(Seurat)

还有一些包会在后面的教程中会用到。

remotes::install_github("satijalab/seurat-data", "seurat5", quiet = T)
remotes::install_github("satijalab/azimuth", "seurat5", quiet = T)
remotes::install_github("satijalab/seurat-wrappers", "seurat5", quiet = T)
remotes::install_github("stuart-lab/signac", "seurat5", quiet = T)
remotes::install_github("bnprks/BPCells", quiet = T)

3基于测序和成像的空间数据集分析

基于测序(即VisiumSLIDE-seq等)和基于成像(MERFISH/VizgenXeniumCosMX等)的空间数据集分析技术都具有独特的优势,需要根据实际情况选择分析方法。😜

本次更新,Seurat v5支持空间数据格式,并支持scRNA-seq整合、反卷积和生态位识别的分析技术。


alt

alt

4分析流程简化

Seurat v5提供一行代码解决整合数据问题:👇

4.1 整合数据

obj <- IntegrateLayers(
object = obj, method = CCAIntegration,
orig.reduction = "pca", new.reduction = "integrated.cca",
verbose = FALSE
)

4.2 聚类

obj <- FindNeighbors(obj, reduction = "integrated.cca", dims = 1:30)
obj <- FindClusters(obj, resolution = 2, cluster.name = "cca_clusters")

4.3 可视化

obj <- RunUMAP(obj, reduction = "integrated.cca", dims = 1:30, reduction.name = "umap.cca")
p1 <- DimPlot(
obj,
reduction = "umap.cca",
group.by = c("Method", "predicted.celltype.l2", "cca_clusters"),
combine = FALSE
)

obj <- FindNeighbors(obj, reduction = "integrated.scvi", dims = 1:10)
obj <- FindClusters(obj, resolution = 2, cluster.name = "scvi_clusters")

obj <- RunUMAP(obj, reduction = "integrated.scvi", dims = 1:10, reduction.name = "umap.scvi")
p2 <- DimPlot(
obj,
reduction = "umap.scvi",
group.by = c("Method", "predicted.celltype.l2", "scvi_clusters"),
combine = FALSE
)

wrap_plots(c(p1, p2), ncol = 2)
alt

5多组学整合

Seurat v5提供了一个新的技术,叫bridge integration,用于整合不同组学的实验结果,例如单独的scRNA-seqscATAC-seq数据集。🥰

这种方法利用了单独的多组学数据集作为分子“桥梁”,在低维空间中进行整合,并返回一个目标降维(例如integrated.rpca),旨在将不同批次中的共享细胞类型共嵌入。🧐

alt

整理一下目前提供的整合算法:👇

  • CCA:基于规范相关性分析的整合,适用于具有 相似细胞类型组成的数据集。
  • RPCA:基于稀疏主成分分析的整合,适用于具有 不同细胞类型组成的数据集。
  • MNN:基于互近邻的整合,适用于具有 相似细胞类型组成的数据集。
  • LIGER:基于因子分析的整合,适用于具有 不同细胞类型组成的数据集。
  • Bridge:基于字典学习的整合,适用于 不同组学的数据集。

6更高的计算性能和更低的硬件要求

scRNA-seq数据集的大小和规模正在迅速增加,甚至超过了摩尔定律。😷

Seurat v5更新后,作者引入了新方法来分析、解释和探索跨越数百万个细胞的数据集,即使它们无法完全加载到内存中。😅

具体而言,可以理解为一种抽样的分析方法,将大型数据集的代表性子样本存储在内存中以实现快速和迭代分析,而完整数据集仍可通过磁盘存储访问。😂

此外,基于Ben ParksGreenleaf实验室开发的BPCells包实现了更高性能的分析。😜

在原文中,作者对高达130 万个脑细胞数据集和来自多项研究的150 万个细胞数据集都进行了测试。🙂

看来以后我们就可以用笔记本电脑分析数据了。😊

alt

alt

7兼容性问题

Seurat v5更新比较大,有人可能会担心兼容该问题,怕与低版本的Seurat不兼容,作者也在原文中解释了这个问题,不用担心,完美兼容。🤩

alt

alt
最后祝大家早日不卷!~

点个在看吧各位~ ✐.ɴɪᴄᴇ ᴅᴀʏ 〰

📍 往期精彩

📍 🤩 WGCNA | 值得你深入学习的生信分析方法!~
📍 🤩 ComplexHeatmap | 颜狗写的高颜值热图代码!
📍 🤥 ComplexHeatmap | 你的热图注释还挤在一起看不清吗!?
📍 🤨 Google | 谷歌翻译崩了我们怎么办!?(附完美解决方案)
📍 🤩 scRNA-seq | 吐血整理的单细胞入门教程
📍 🤣 NetworkD3 | 让我们一起画个动态的桑基图吧~
📍 🤩 RColorBrewer | 再多的配色也能轻松搞定!~
📍 🧐 rms | 批量完成你的线性回归
📍 🤩 CMplot | 完美复刻Nature上的曼哈顿图
📍 🤠 Network | 高颜值动态网络可视化工具
📍 🤗 boxjitter | 完美复刻Nature上的高颜值统计图
📍 🤫 linkET | 完美解决ggcor安装失败方案(附教程)
📍 ......

本文由 mdnice 多平台发布

  • 3
    点赞
  • 15
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
单细胞测序(single-cell sequencing)是一种高通量测序技术,可以对单个细胞的基因组或转录组进行全面的分析。而R语言是一种广泛用于统计分析和数据可视化的编程语言。 在单细胞测序实验中,通过测序技术可以获取到大量的细胞的基因表达数据,包括每个细胞中数以千计的基因的表达水平。而这些数据的处理和分析就需要使用到R语言以及相关的数据分析包和函数。 首先,我们可以使用R语言中的数据读取函数将单细胞测序的原始数据导入到R环境中,并进行数据清洗和预处理。例如,可以通过R的数据处理包如‘Seurat’对单细胞数据进行降噪、标准化和归一化等处理,以确保数据的准确性和可靠性。 接下来,我们需要使用R语言的统计分析技术对这些单细胞数据进行分析。例如,可以通过差异表达分析(DEG)来寻找在不同细胞类型或条件下差异表达的基因。也可以使用聚类算法将细胞进行分组,寻找不同细胞群体之间的差异和相似性。 同时,R语言还提供了多种数据可视化的方法,我们可以使用R语言中的绘图包如‘ggplot2’和‘pheatmap’等对单细胞测序数据进行可视化。可视化可以帮助我们更直观地展示细胞群体的分布情况、基因表达的模式等,从而更好地理解和解释实验结果。 总而言之,单细胞测序数据的R语言分析可以帮助我们深入理解细胞的表达特征和功能,发现新的细胞类型和亚群体,并为研究细胞发育、疾病机制等提供重要的生物学信息。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值