Seurat Weekly NO.07 || V4 新特性

最近(202101)安装Seurat(3.2.3)的话可能会收到一条提示信息:

Seurat v4 will be going to CRAN in the near future;
 for more details, please visit https://satijalab.org/seurat/v4_changes

我们介绍过的Seurat V4版就要通过CRAN的规范,可以在CRAN上直接安装了,这意味着V4走出了内侧和\beta版,以更成熟的方式提供给广大单细胞研究的科学家。我们在Seurat4.0|| 单细胞多模态数据分析启示录系列笔记中学习了V4整合多模态数据的过程,我们应该记得V4在预印的文章:Integrated analysis of multimodal single-cell data。

Seurat 4.0 ||单细胞数据分析工具箱有更新
Seurat 4.0 ||单细胞多模态数据整合算法WNN
Seurat 4.0 || 分析scRNA和表面抗体数据
Seurat 4.0 || WNN整合scRNA和scATAC数据
Seurat 4.0 || 单细胞PBMC多模态参考数据集
Seurat 4.0 || 单细胞BMNC多模态参考数据集

除了以上新功能的集成,Seurat还做了已有函数的小改动,大部分只是性能的提升,尽可能地保证了与Seurat V3 的兼容性。这些修改包括对默认参数设置的微小更改以及为相同任务使用性能更好的包,例如标识k最近的邻居和基于图的聚类。这些变化不会对下游结果产生负面影响。

默认参数的改变

这当然是值得注意的,在不少时候会困惑:为什么同样的代码,同样的数据,在小明那里是这样的,而我是另一个样!默认参数,随机种子都是可能的原因,或者没有仔细检查数据和代码(也就是你说的一样,其实并不一样,坏了啦)。

  • FindNeighbors

用于识别k近邻的默认方法被设置为annoy。这是一种近似最近邻方法,广泛应用于高维分析,包括单细胞分析。广泛的社区基准测试表明,annoy极大地提高了邻居发现的速度并减少了内存需求,对下游结果的影响可以忽略不计(当然会有些许不同),这与我们的分析和测试是一致的。用户可以使用nn.method="rann"切换回原来的默认设置。

?FindNeighbors 
## Default S3 method:
FindNeighbors(
  object,
  distance.matrix = FALSE,
  k.param = 20,
  compute.SNN = TRUE,
  prune.SNN = 1/15,
  nn.method = "rann",  # V3 的默认
  annoy.metric = "euclidean",
  nn.eps = 0,
  verbose = TRUE,
  force.recalc = FALSE,
  ...
)
  • FindMarkers

我们重新构造了FindMarkers()函数的代码,使其更易于理解、解释和调试。差异表达的结果保持不变。然而,在默认情况下,我们现在使用的FC(fold change )则基于log2,这在其他差异表达包中是常见的,而不是V3之前的自然对数。如果设置了默认选项,FindMarkers()的输出将包括列avg_log2FC,而不是avg_logFC。用户可以通过指定base = exp(1)来恢复以前的行为(以自然对数味为底的FC)。

  • IntegrateData/TransferData

我们对锚点加权矩阵(anchor weight matrix)的精确计算做了微小的改变。这些更改反映了改进的工作流,但不会对下游分析产生有意义的差异(例如,您可以使用Seurat v3和Seurat v4来比较)

  • SCTransform

在SCTransform()中,我们稍微修改了默认参数,以提高大型数据集参数估计的可伸缩性(scalability )。例如,在估计mu和theta之间的正则关系时,我们通过将ncells参数设置为5000来计算数据子集。sctransform v0.3中的vst()函数(可在CRAN上获得)还对正则化过程进行了微小的更改。我们对这些更改进行了广泛的测试,发现在速度和内存方面有了实质性的改进,特别是对于大型数据集,而且没有对结果产生不利影响。用户可以比较使用Seurat v3和Seurat v4计算的SCTransform vignette的结果,或者在更大的数据集中设置ncells=NULL来比较结果。

移除的函数

The following functions have been removed in Seurat v4:

  • CreateGeneActivityMatrix replaced by GeneActivity in Signac
  • RunLSI replaced by RunTFIDF and RunSVD in Signac
  • ReadAlevin and ReadAlevinCsv moved to SeuratWrappers, see details here
  • ExportToCellbrowser and StopCellbrowser moved to SeuratWrappers, see details here

英国哲学家弗朗西斯·培根在1620年出版了《新工具》一书,原名《新工具或解释自然的一些指导》,并作为《伟大的复兴》的第二部分。我们相信任何一个学科门类,工具既是技术发展的缩影也是科学发展的推手,随着Seurat的迭代,我们也感受到了单细胞基因组学的发展轨迹。



©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 213,186评论 6 492
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 90,858评论 3 387
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 158,620评论 0 348
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 56,888评论 1 285
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,009评论 6 385
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,149评论 1 291
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,204评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 37,956评论 0 268
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,385评论 1 303
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,698评论 2 327
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,863评论 1 341
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,544评论 4 335
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,185评论 3 317
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,899评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,141评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,684评论 2 362
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,750评论 2 351

推荐阅读更多精彩内容