在之前的一個視頻學(xué)習(xí)系列筆記里蓝晒，有一篇筆記提到了RNA velocity（Single cell RNA-seq data analysis with R視頻學(xué)習(xí)筆記（八）），但當(dāng)時沒有仔細的去對RNA velocity進行詳細的了解除呵。這篇文獻學(xué)習(xí)筆記就是對這個模型進行更詳細的了解。
文獻：RNA velocity of single cells
發(fā)表于2018年，Nature壁却。
如果想下載這篇文獻瘤礁，而又沒Nature官網(wǎng)權(quán)限的童鞋可以在這里下載：here
（筆記里很多補充圖沒有放上來阳懂，想看的可以從上面鏈接下載文章看）
【摘要】
RNA豐度是單個細胞狀態(tài)的一個強大的指標(biāo)。單細胞RNA測序可以精確的柜思、高通量的揭示RNA豐度岩调。然而，這種方法只能捕捉某一個時間點的“快照”（snapshot）赡盘，對于分析胚胎形成或者組織重建的這種需要觀測連續(xù)時間的“事件”就比較困難了号枕。這里我們提出了“RNA velocity”（RNA速度），可以直接測定在單細胞測序結(jié)果中的“未成熟”mRNA和“成熟”mRNA陨享。RNA velocity是一個高維度載體葱淳，可以預(yù)測單個細胞在時間軸上的“未來的”幾個小時的狀態(tài)。我們利用神經(jīng)嵴系驗證了RNA velocity的精確度抛姑，證明了它可以用在多種測序平臺上赞厕，利用它揭示了發(fā)育中的小鼠海馬體的分支譜系樹，并測定了人類腦胚胎的轉(zhuǎn)錄動力學(xué)定硝。
【正文】
在發(fā)育過程中皿桑，分化過程發(fā)生在時間軸上的幾小時到數(shù)天不等。而未成熟的mRNA（unspliced）和成熟的mRNA（spliced）的相對豐度可以評估基因剪切和降解的速率。我們推斷這種信號可以在單細胞測序數(shù)據(jù)中被檢測到诲侮，從而可以揭示在動態(tài)的過程中镀虐，整個轉(zhuǎn)錄的速率和變化。
所有的單細胞RNA測序方法都是依賴于oligo-dT引物浆西，通過富集帶有polyA的mRNA分子來得到測序材料粉私。盡管如此，現(xiàn)有的測序平臺（SMART-seq2, STRT/
C1, inDrop and 10x Genomics）的測序結(jié)果中近零，都含有15-25%的reads包含未成熟帶有內(nèi)含子的序列（圖1a）诺核。在常規(guī)bulk RNA-seq里含有14.6%的未成熟mRNA reads，單細胞測序中含有20%的未成熟mRNA reads久信。在10x Genomics文庫里窖杀，這一類mRNA的reads主要是由于PCR擴增時候產(chǎn)生的。
為了量化時間依賴的mRNA前提和成熟的mRNA裙士，我們給轉(zhuǎn)錄動力學(xué)假定了一個簡單的模型入客，成熟的mRNA豐度由未成熟mRNA產(chǎn)生的spliced mRNA和降解的速度來決定（圖1b）。在這個動態(tài)的過程中腿椎，轉(zhuǎn)錄速率a上升桌硫，導(dǎo)致了unspliced的mRNA增加，從而導(dǎo)致成熟的mRNA的增多（圖1c）啃炸，直到一個新的穩(wěn)定狀態(tài)铆隘。反之亦然。在基因表達的誘導(dǎo)過程里南用，未成熟的mRNA的存在超出了預(yù)期（圖1d）膀钠。這兩種mRNA的平衡是在未來狀態(tài)里成熟mRNA豐度的一個重要指標(biāo)，也是意味著決定了細胞未來的狀態(tài)裹虫。
為了證明這個模型可以用來推斷在“未來時間里”成熟的mRNA豐度肿嘲，我們用bulk-RNA-seq方法檢測了小鼠肝臟的日周期的一個時間過程。未成熟的mRNA水平在每一個時間點都與隨后一個時間點上成熟的mRNA水平相似（圖1e）筑公，許多與節(jié)律相關(guān)的基因也展示出了unspliced mRNA的上調(diào)的趨勢（圖1f,g）雳窟。這一模型可以讓我們在每一個周期里進行測定，準(zhǔn)確捕捉晝夜節(jié)律周期發(fā)展的預(yù)期方向(圖1h)匣屡。
接下來涩拙，為了證明在單細胞水平上，這個模型也可以預(yù)測轉(zhuǎn)錄動力學(xué)耸采，我們分析了近期發(fā)表的老鼠嗜鉻細胞的數(shù)據(jù)，是SMARTseq2平臺測序的（圖2a）工育。在發(fā)育過程中虾宇，相當(dāng)比例的嗜鉻細胞（即腎上腺髓質(zhì)的神經(jīng)內(nèi)分泌細胞）源自Schwann細胞前體，為我們的模型提供了一個簡單的測試如绸，并可以通過譜系跟蹤來驗證嘱朽。圖2b和c里綠色/灰色的圖代表著基因的表達情況旭贬，綠色代表表達水平高；而紅色/藍色圖則表達未成熟mRNA和成熟mRNA搪泳，說明Serpin2這個基因是處在一個被抑制的狀態(tài)稀轨，而Chga基因處于被激活的狀態(tài)。RNA速度對單個細胞進行評估岸军，準(zhǔn)確地再現(xiàn)了該數(shù)據(jù)的轉(zhuǎn)錄動力學(xué)奋刽，包括細胞的分化，比如細胞向嗜鉻細胞方向發(fā)展(圖2d)艰赞。RNA速度還捕獲了嗜鉻的細胞-周期動力學(xué)佣谐，包括PCA分析和細胞周期相關(guān)基因的分析。
有很多的技術(shù)可以對velocity在低維里進行可視化方妖。細胞狀態(tài)的推測可以嵌入一個低維空間里（比如PCA狭魂，圖2d）或者tSNE（圖2h）。在一個很大的數(shù)據(jù)集了党觅，RNA velocity更容易被可視化（圖Fig. 2i）雌澄。因為細胞可以同時擁有多個獨立的RNA velocity存在，比如分化杯瞻、成熟和增殖镐牺。
細胞特異性的RNA velocity評估提供了一個自然的基礎(chǔ)對細胞命運進行定量建模。有效的可推測的時間軸又兵，取決于你要研究的生物過程任柜。比如用EdU標(biāo)記的嗜鉻祖細胞，我們可以推測的“未來”時間是2.5-3.8小時（圖Fig. 2f, g）沛厨。由于這種推測是線性的宙地，所以推測時間軸也依賴于基因表達軌跡的形狀。
我們接著利用RNA velocity對小鼠的海馬體的風(fēng)雨譜系的分支進行分析逆皮。去除血管宅粥、免疫細胞、GABAergic 和卡哈爾Retzius神經(jīng)元后电谣，t-SNE圖揭示了一個復(fù)雜的分支（圖3a）秽梅。我們使用已知的標(biāo)記基因鑒定了星形細胞、少突細胞前體 (OPCs)剿牺、齒狀回顆粒神經(jīng)元和錐體神經(jīng)元五個細胞群：the subiculum, CA1, CA2, CA3企垦，hilus。單個基因的相位圖展示了特異的基因表達的誘導(dǎo)和抑制（圖3b）晒来。采用馬爾可夫速度隨機游走模型可以自動識別分支的終端和root(圖3c)钞诡，展示了RNA速度可以在沒有已知的發(fā)育過程的前提下，確定譜系的方向。
在微觀層面上荧降，“命運”選擇是一個傾向于非決定性過程接箫，這個過程涉及到基因的表達，一旦轉(zhuǎn)錄因子的反饋loop建立后朵诫，細胞的“命運”便被“鎖定”了辛友。比較一個pre-OPC群里的細胞和一個處于過度期（narrow passage）的細胞的未來狀態(tài)的可能性分布，顯然后者更可能成為一個OPC細胞（圖 3d）剪返。在圖3e里废累，單細胞分支譜系里，如果一個細胞的Prox1基因被激活随夸，那么變形成granule神經(jīng)元九默，如果Prox1基因被抑制，那么細胞會分化成neuroblasts宾毒。
為了證明RNA速度可以在人體胚胎內(nèi)被檢測到驼修，我們進行了基于液滴的單細胞RNA-seq的實驗，分析了在懷孕十周后人類前腦的發(fā)育過程诈铛，主要研究谷氨酸能神經(jīng)元譜系(圖4a)乙各。我們發(fā)現(xiàn)由增殖的前體細胞狀態(tài)(radial glia)通過一系列中間成神經(jīng)細胞的狀態(tài)向更成熟的分化的谷氨酸能神經(jīng)元（表達SLC17A7）產(chǎn)生了一個很強的velocity模式。我們利用多重原位雜交實驗幢竹，驗證了皮質(zhì)神經(jīng)元發(fā)育的已知和新的標(biāo)記物耳峦，(圖4b, c)。這些在組織中分層的基因(圖4c)與單細胞RNA-seq數(shù)據(jù)中基因的偽時間分布密切相關(guān)（圖4b）焕毫。
我們利用PCA來把細胞根據(jù)分化的“擬時間”進行排序蹲坷，確定了未成熟mRNA始終先出現(xiàn)于成熟的mRNA（圖4d）。從4d里還可以看到不同的動力學(xué)特征邑飒。比如RNASEH2B基因?qū)儆趂ast動力學(xué)循签，因為它的unspliced和spliced RNA差不多。相反的疙咸，DCX县匠，ELAVL4 和STMN2基因在剛開始有一個快速的轉(zhuǎn)錄，隨后的轉(zhuǎn)錄速度下降（藍色線下降）撒轮，成熟的mRNA有一個明顯的滯后軌跡乞旦。
由于RNA速度是基于真實的轉(zhuǎn)錄動力學(xué)，這方法有望為我們的研究帶來更堅實的定量基礎(chǔ)题山，幫助我們了解細胞在基因表達空間中的動態(tài)分化兰粉。RNA速度已經(jīng)可以用來進行詳細的對整個有機體的動態(tài)過程的研究，為譜系分析提供很大的便利顶瞳，特別是在人類胚胎中亲桦。
