Seurat Weekly NO.17 || 一個R包的CNS之路

seurat發(fā)展歷史田轧,cns之路

在單細(xì)胞數(shù)據(jù)分析中舔糖,特別是對初學(xué)者來說帮毁,Seurat是一個教科書級別的工具酌摇。至于這個用來分析單細(xì)胞多組學(xué)的工具為什么叫這個名字瀑志,我想汁果,與它問世的第一篇文獻(xiàn)不無關(guān)系:

這篇文章的第一作者Rahul Satija和作者Aviv Regev成為日后單細(xì)胞領(lǐng)域人們耳熟能詳?shù)拿峙韧鳶eurat是一位畫家的名字:

喬治·修拉(Georges Seurat驶赏,1859-1891)耕驰,1859年12月2日生于巴黎爷辱,受過完整的美術(shù)學(xué)院教育,曾師從安格爾的學(xué)生亨利·萊曼(Henri Lehmann)學(xué)習(xí)古典主義繪畫朦肘,后來又研究過 盧浮宮中的大師作品饭弓,對光學(xué)和色彩理論特別關(guān)注并為之做了大量的實驗。他的畫作風(fēng)格相當(dāng)與眾不同媒抠,Seurat的畫充滿了細(xì)膩繽紛的小點弟断,當(dāng)你靠近看,每一個點都充滿著理性的筆觸领舰,與梵高的狂野夫嗓,還有塞尚的色塊都大為不同迟螺。把文藝復(fù)興傳統(tǒng)的古典結(jié)構(gòu)和印象主義的色彩試驗結(jié)合起來。Seurat擅長畫都市中的風(fēng)景畫舍咖,也擅長將色彩理論套用到畫作當(dāng)中矩父。把最新的繪畫空間概念、傳統(tǒng)的幻象透視空間排霉、以及在色彩和光線的知覺方面的最新科學(xué)發(fā)現(xiàn)結(jié)合起來窍株。

而這也是一篇描繪空間分析的單細(xì)胞文章NBT。我們知道攻柠,空間定位是細(xì)胞命運和行為的關(guān)鍵決定因素球订,但在復(fù)雜組織中,空間解析瑰钮、轉(zhuǎn)錄組范圍內(nèi)基因表達(dá)譜的方法還缺乏冒滩。RNA染色方法只檢測少量轉(zhuǎn)錄本,而測量基因表達(dá)的scRNA-seq將細(xì)胞從其固有的空間環(huán)境中分離出來浪谴。正式在這篇文章中开睡,Seurat以一個附件的形式面世了,Seurat一種通過整合單細(xì)胞RNA-seq數(shù)據(jù)和原位RNA模式來推斷細(xì)胞定位的計算策略苟耻。


這里的幾個附件都值得讀一讀篇恒,如Supplementary Text and Figures中講述了空間分析的大部分關(guān)鍵點也給出了第一套Seurat教程:

在這個意義上,Serurat是為空間分析而生的凶杖。然后為了我們更好地理解Seurat胁艰,第一版的源代碼我們還是需要拜讀一下的:

Seurat_Source\seurat.R
Seurat_Source\seuratFxns.R
Seurat_Source\zfRenderSeurat.R

在2021年的Seurat的更新日志中,我們中不到SeuratV1的信息智蝠,這個日志是從Seurat 2.0.0開始記錄的腾么。V2是Seurat功能和形象的一大轉(zhuǎn)折和重新定型,盡管保留了空間分析的功能杈湾,但更多地是轉(zhuǎn)向單細(xì)胞多組學(xué)數(shù)據(jù)的整合分析哮翘。

隨著單細(xì)胞測序通量和緯度的提升,識別跨多個數(shù)據(jù)集(實驗條件毛秘、技術(shù)平臺饭寺、不同物種)的細(xì)胞亞群成為一個挑戰(zhàn)。Seurat V2基于公共變異源(common sources of variation)集成scRNA-seq數(shù)據(jù)集的分析策略叫挟,允許跨數(shù)據(jù)集識別共享種群并進(jìn)行下游比較分析艰匙。這即是日后人們在做數(shù)據(jù)整合(Integration)時說到的CCA的方法實現(xiàn)途徑。

  • Stuart, Butler, et al., Cell 2019 [Seurat V3]
    僅僅一年之后抹恳,在整合海量單細(xì)胞數(shù)據(jù)這條路上员凝,Seurat越走越遠(yuǎn),隨著單細(xì)胞數(shù)據(jù)解析的組織類型越來越多奋献,僅僅用Seurat的CCA還是有很多數(shù)據(jù)集無法整合到一起健霹。于是旺上,整合力度更大的算法出現(xiàn)了,策略也有所變化:直接以一個作為reference另一個作為query糖埋。這個概念像極了基因組數(shù)據(jù)分析中常用的mapping宣吱、alignment、blast瞳别。

單細(xì)胞轉(zhuǎn)錄組學(xué)已經(jīng)改變了我們描述細(xì)胞狀態(tài)的能力征候,但深入的生物學(xué)理解需要的不僅僅是亞群的分類。隨著測量不同細(xì)胞模式技術(shù)的出現(xiàn)祟敛,一個關(guān)鍵的分析挑戰(zhàn)是整合這些數(shù)據(jù)集疤坝,以更好地了解細(xì)胞的身份和功能。Seurat V3將不同的數(shù)據(jù)集“錨定”在一起馆铁,使我們不僅能夠集成跨scRNA-seq技術(shù)的單細(xì)胞測量跑揉,而且能夠集成跨不同模式的單細(xì)胞測量。

  • Hao, Hao, et al., bioRxiv 2020 [Seurat V4]
    2020年10月埠巨,紐約基因組所( New York Genome Center )的 Satija Lab團(tuán)隊開發(fā)的單細(xì)胞分析工具Seurat升級到了4.0畔裕,其特性之一是開發(fā)了加權(quán)鄰近(weighted-nearest neighbor,WNN)算法來分析多模態(tài)數(shù)據(jù)乖订。所謂多模態(tài)(multimodal )是指同時在一個細(xì)胞內(nèi)測量分屬于(廣義?)中心法則不同過程的特征(如具练,RNA乍构,ATAC,膜蛋白等)扛点。
seurat引用情況哥遮,技術(shù)快風(fēng)

與2017年以來的Seurat的快速更新同時發(fā)生的是,單細(xì)胞技術(shù)的快速應(yīng)用與普及陵究。簡單的幾行命令可以看出眠饮,2021年以來,大部分時間铜邮,一日之內(nèi)就有上千人下載和和安裝Seurat:

# install.packages("cranlogs")
library(cranlogs)
tail(Seuratdl <- cran_downloads(from = "2017-08-01",  to= "2021-06-01", packages = c("Seurat")))
          date count package
1396 2021-05-27  1441  Seurat
1397 2021-05-28  1101  Seurat
1398 2021-05-29   499  Seurat
1399 2021-05-30   487  Seurat
1400 2021-05-31  1022  Seurat
1401 2021-06-01  1316  Seurat
plot(Seuratdl$date,Seuratdl$count)
seurat開發(fā)

Seurat的廣泛應(yīng)用既反映了單細(xì)胞組學(xué)的熱度仪召,也說明Seurat團(tuán)隊在開源這件是上做的很好。github上每天都有人提問交流松蒜,開發(fā)者也是其中的積極參與者扔茅。讓我們簡單回顧一下這個R包的發(fā)展史,感受一下時間的力量秸苗。

Seurat Changelog

正式因為對單細(xì)胞數(shù)據(jù)分析有著深刻的理解和與用戶的廣泛交流召娜,Seurat單細(xì)胞數(shù)據(jù)分析工具箱日漸完善。成為大部分入門單細(xì)胞數(shù)據(jù)分析的首選工具惊楼。Seurat不僅教程豐富玖瘸,還附有大量可以直接練習(xí)分析的數(shù)據(jù)集和對其他工具的友好接口秸讹,如Velocity、LIGER雅倒、Monocle3璃诀、Harmony、alevin屯断、Nebulosa等文虏。

seurat數(shù)據(jù)結(jié)構(gòu),可擴(kuò)展殖演,數(shù)據(jù)管理

既然定位于初學(xué)者氧秘,或者說是單細(xì)胞探索性數(shù)據(jù)分析工具,那我們就來看看Seurat在除了質(zhì)控趴久,降維聚類丸相,差異分析,整合分析之外一個十分重要的功能:數(shù)據(jù)管理彼棍。

基于R語言的S4結(jié)構(gòu)灭忠,Seurat構(gòu)建了一個以表達(dá)矩陣為核心的單細(xì)胞多組學(xué)數(shù)據(jù)結(jié)構(gòu)。這樣座硕,一個Seurat其實就是一個study的實驗設(shè)計弛作。

這個數(shù)據(jù)結(jié)構(gòu)也是大部分剛接觸R語言的同學(xué)最容易困惑的:

  • 為什么函數(shù)輸入的是Seurat對象,輸出的也還是Seurat呢华匾?
  • 降維的數(shù)據(jù)它是如何調(diào)用的呢映琳?數(shù)據(jù)又藏在哪里呢?
  • 更多找不到
開啟你的seurat之旅蜘拉,最短入門路徑

最短入門路徑是:

  • 安裝R語言
  • 安裝Seurat
  • 安裝SeuratData
  • 打開Seurat網(wǎng)頁
  • 練習(xí)單細(xì)胞數(shù)據(jù)分析
  • 分析自己課題組的數(shù)據(jù)
  • 發(fā)表文獻(xiàn)時引用Seurat
  • 寫自己的SeuratWrappers包

先練習(xí)再分析萨西,先用SeurtData的數(shù)據(jù)集練習(xí)教程理解單細(xì)胞分析中的概念和Seurat的數(shù)據(jù)分析結(jié)構(gòu)。


https://satijalab.org/seurat/
https://satijalab.org/seurat/news/index.html

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市持寄,隨后出現(xiàn)的幾起案子源梭,更是在濱河造成了極大的恐慌,老刑警劉巖稍味,帶你破解...
    沈念sama閱讀 206,013評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件咸产,死亡現(xiàn)場離奇詭異,居然都是意外死亡仲闽,警方通過查閱死者的電腦和手機脑溢,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,205評論 2 382
  • 文/潘曉璐 我一進(jìn)店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人屑彻,你說我怎么就攤上這事验庙。” “怎么了社牲?”我有些...
    開封第一講書人閱讀 152,370評論 0 342
  • 文/不壞的土叔 我叫張陵粪薛,是天一觀的道長。 經(jīng)常有香客問我搏恤,道長违寿,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 55,168評論 1 278
  • 正文 為了忘掉前任熟空,我火速辦了婚禮藤巢,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘息罗。我一直安慰自己掂咒,他們只是感情好,可當(dāng)我...
    茶點故事閱讀 64,153評論 5 371
  • 文/花漫 我一把揭開白布迈喉。 她就那樣靜靜地躺著绍刮,像睡著了一般。 火紅的嫁衣襯著肌膚如雪挨摸。 梳的紋絲不亂的頭發(fā)上孩革,一...
    開封第一講書人閱讀 48,954評論 1 283
  • 那天,我揣著相機與錄音得运,去河邊找鬼膝蜈。 笑死,一個胖子當(dāng)著我的面吹牛澈圈,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播帆啃,決...
    沈念sama閱讀 38,271評論 3 399
  • 文/蒼蘭香墨 我猛地睜開眼瞬女,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了努潘?” 一聲冷哼從身側(cè)響起诽偷,我...
    開封第一講書人閱讀 36,916評論 0 259
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎疯坤,沒想到半個月后报慕,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,382評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡压怠,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 35,877評論 2 323
  • 正文 我和宋清朗相戀三年眠冈,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 37,989評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡蜗顽,死狀恐怖布卡,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情雇盖,我是刑警寧澤忿等,帶...
    沈念sama閱讀 33,624評論 4 322
  • 正文 年R本政府宣布,位于F島的核電站崔挖,受9級特大地震影響贸街,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜狸相,卻給世界環(huán)境...
    茶點故事閱讀 39,209評論 3 307
  • 文/蒙蒙 一薛匪、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧卷哩,春花似錦蛋辈、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,199評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至尊浓,卻和暖如春逞频,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背栋齿。 一陣腳步聲響...
    開封第一講書人閱讀 31,418評論 1 260
  • 我被黑心中介騙來泰國打工苗胀, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人瓦堵。 一個月前我還...
    沈念sama閱讀 45,401評論 2 352
  • 正文 我出身青樓基协,卻偏偏與公主長得像,于是被迫代替她去往敵國和親菇用。 傳聞我的和親對象是個殘疾皇子澜驮,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 42,700評論 2 345

推薦閱讀更多精彩內(nèi)容