鑒定好細(xì)胞類型后臊恋，我們還可以對(duì)其進(jìn)行一些功能性分析衣洁，如擬時(shí)序分析和SCENIC分析，對(duì)細(xì)胞的生物學(xué)意義進(jìn)行一定程度上的挖掘抖仅。本期小編主要為大家介紹擬時(shí)序分析（Pseudotime Analysis）坊夫。

1****、Q**** & ****A環(huán)節(jié)

Q****1****：我們?yōu)槭裁匆M(jìn)行擬時(shí)序分析撤卢？环凿？

l 機(jī)體為響應(yīng)各種應(yīng)激，其細(xì)胞會(huì)從一種功能“狀態(tài)”轉(zhuǎn)變?yōu)榱硪环N功能“狀態(tài)”放吩；

l 當(dāng)細(xì)胞在不同狀態(tài)之間轉(zhuǎn)變時(shí)智听，往往會(huì)經(jīng)歷轉(zhuǎn)錄重組，導(dǎo)致一些基因被沉默屎慢，一些基因被重新激活瞭稼，但純化這些瞬態(tài)細(xì)胞進(jìn)行研究是很困難或不可能的；

l ScRNA-seq擬時(shí)序分析可以讓我們?cè)诓恍枰兓那闆r下查看這些細(xì)胞狀態(tài)腻惠。

Q****2****：Sc****RNA-seq****擬時(shí)序分析是什么环肘？？

擬時(shí)序分析集灌，即根據(jù)不同細(xì)胞亞群基因表達(dá)量隨時(shí)間的變化情況悔雹，構(gòu)建細(xì)胞譜系發(fā)育，但這里的時(shí)間并不是真時(shí)間欣喧，而是一個(gè)虛擬的時(shí)間腌零，是指的細(xì)胞與細(xì)胞之間的轉(zhuǎn)化和演替的順序和軌跡。

l 即使在同一個(gè)樣本中唆阿，也會(huì)存在多種不同的細(xì)胞形態(tài)益涧。因此，不論測(cè)定了多少樣本驯鳖，我們都可以采用擬時(shí)序分析對(duì)樣本中的細(xì)胞轉(zhuǎn)化和變化進(jìn)行描述闲询。

Q****3****：擬時(shí)序分析用的工具是啥久免？它能進(jìn)行哪些分析？

l Monocle是用于ScRNA-seq擬時(shí)序分析的經(jīng)典工具（R包）扭弧，目前已更新至3版本阎姥；

l Monocle是使用算法來(lái)學(xué)習(xí)細(xì)胞狀態(tài)轉(zhuǎn)變過(guò)程中每個(gè)細(xì)胞必須經(jīng)歷的基因表達(dá)變化序列，一旦了解了基因表達(dá)變化的整體“軌跡”鸽捻，Monocle就可以將每個(gè)細(xì)胞放置在軌跡中的適當(dāng)位置呼巴。

l 我們可以使用Monocle中的相關(guān)分析工具做Beam分析（Branch Expression Analysis Modeling），揭示重要的基因和細(xì)胞御蒲。

（算法參考文獻(xiàn)：Reversed graph embedding resolves complex single-cell trajectories.）

2衣赶、擬時(shí)序分析工作流搭建 略

3、擬時(shí)序分析結(jié)果解讀

對(duì)于擬時(shí)序分析來(lái)說(shuō)删咱，研究某一特定細(xì)胞類型的轉(zhuǎn)化屑埋，如M1/M2型巨噬細(xì)胞極化、CD8⁺ T細(xì)胞激活和耗竭等等痰滋，往往具有一定的生物學(xué)意義摘能。這里小編通過(guò)對(duì)上期T細(xì)胞進(jìn)行亞群鑒定，以CD8⁺ T細(xì)胞亞型為例進(jìn)行擬時(shí)序分析（圖1）敲街。

image.png

圖1 CD8⁺ T細(xì)胞亞群鑒定流程

（1）CD****8****⁺****T細(xì)胞樹形結(jié)構(gòu)軌跡

image.png

圖2 CD8⁺ T細(xì)胞樹形結(jié)構(gòu)軌跡

圖中每個(gè)點(diǎn)代表一個(gè)細(xì)胞团搞，具有相似狀態(tài)的細(xì)胞被聚到一起，每個(gè)分支點(diǎn)代表一個(gè)可能的細(xì)胞生物學(xué)過(guò)程決策點(diǎn)（該例圖中只有1個(gè)分支點(diǎn)）多艇。圖2展示了CD8⁺ T細(xì)胞每個(gè)cluster在偽時(shí)間軸上的分布逻恐，不同顏色代表不同cluster。我們這邊通過(guò)對(duì)每個(gè)cluster進(jìn)行細(xì)胞類型鑒定峻黍，初步將cluster 2/3/6鑒定為記憶T細(xì)胞复隆，cluster 7鑒定為Na?ve T細(xì)胞。

（2）CD****8****⁺****T細(xì)胞所處分化時(shí)間軌跡

Monocle基于離默認(rèn)偽時(shí)間軸起點(diǎn)的遠(yuǎn)近姆涩，將細(xì)胞排布如圖3所示挽拂。在該圖中顏色越深代表默認(rèn)的起點(diǎn)，顏色越淺表示離偽時(shí)間軸起點(diǎn)越遠(yuǎn)骨饿。注意：這邊的起點(diǎn)是monocle算法計(jì)算出來(lái)的亏栈，而非實(shí)際真正的起點(diǎn)。

image.png

圖3 CD8⁺ T細(xì)胞所處分化時(shí)間的軌跡圖

以該例圖為例宏赘，我們結(jié)合每個(gè)分支（共3個(gè)分支）中的cluster類群是否高表達(dá)CD8⁺ Na?ve T細(xì)胞相關(guān)Marker（如LEF1绒北、SELL和CCR7等）, 耗竭性T細(xì)胞相關(guān)Marker（如LAG3、HAVCR2察署、CD27闷游、TIGIT等），以及細(xì)胞增殖相關(guān)Marker（如MKI67、TOP2A脐往、CCNB1等）俱济，將起始狀態(tài)鑒定為Na?ve T細(xì)胞，終末狀態(tài)為耗竭性（Cell fate 1）和增殖型（Cell fate 2）T細(xì)胞钙勃。

（3）Beam分析

確定了分化起點(diǎn)后，Monocle可以模擬出每個(gè)細(xì)胞所處的分化時(shí)間聂喇，并尋找隨著分化時(shí)間逐漸升高或降低的基因表達(dá)熱圖和分布圖辖源，即Beam分析。圖4將Top200差異基因聚成4類希太，展示了pre-branch向cell fate1和cell fate2狀態(tài)分化相關(guān)的命運(yùn)決定基因克饶。圖5展示了與細(xì)胞增殖相關(guān)的基因，如MKI67和TOP2A在cluster 5中高表達(dá)誊辉，提示cluster 5作為CD8⁺ Na?ve T細(xì)胞終末分化狀態(tài)是一種增殖型細(xì)胞（Proliferating）矾湃。除此之外，如果有特別關(guān)注的基因也可以基于關(guān)注的基因進(jìn)行繪制堕澄。

image.png

圖4 擬時(shí)序Beam分析基因表達(dá)熱圖

image.png

圖5 擬時(shí)序Beam分析基因表達(dá)分布圖

總的來(lái)說(shuō)邀跃，通過(guò)對(duì)CD8+T細(xì)胞亞型進(jìn)行擬時(shí)序分析得到了CD8+ T細(xì)胞的分化軌跡，進(jìn)行Beam分析可以幫助得到與狀態(tài)分化相關(guān)的命運(yùn)決定基因