研究單位：西安交通大學(xué)基礎(chǔ)醫(yī)學(xué)院院長張保軍教授團(tuán)隊(duì)

期刊名：cellular & molecular immunology

影響因子：21.8

樣本類型：8-12周齡小鼠

文章采用技術(shù)：空間轉(zhuǎn)錄組測序、單細(xì)胞RNA測序（scRNA-seq）

引言

在對抗感染和腫瘤的免疫戰(zhàn)場上，記憶?CD8+?T?細(xì)胞是守護(hù)機(jī)體的「長效衛(wèi)士」，?它們能快速識別再次入侵的病原體或癌細(xì)胞，發(fā)動強(qiáng)力攻擊。長久以來，樹突狀細(xì)胞（DC）被認(rèn)為是主導(dǎo)記憶?T?細(xì)胞分化的核心「指揮官」。研究人員通過整合單細(xì)胞RNA測序（scRNA-seq）和空間轉(zhuǎn)錄（BMKMANU?S1000）技術(shù)，揭示了單核細(xì)胞才是調(diào)控記憶?CD8+?T?細(xì)胞分化的關(guān)鍵「操盤手」，其奧秘藏在「細(xì)胞接觸」與「分子信號」的雙重機(jī)制中。其中BMKMANU?S1000空間轉(zhuǎn)錄組學(xué)平臺大放異彩，成為揭CCR2+單核細(xì)胞促進(jìn)記憶CD8+?T細(xì)胞分化的關(guān)鍵工具。

研究背景

CD8+?T細(xì)胞在適應(yīng)性免疫中扮演著關(guān)鍵角色，能夠保護(hù)機(jī)體免受病原體感染和消除惡性細(xì)胞。當(dāng)CD8+?T細(xì)胞識別到抗原后，會迅速增殖并分化為效應(yīng)CD8+?T細(xì)胞和記憶CD8+?T細(xì)胞。效應(yīng)CD8+?T細(xì)胞負(fù)責(zé)即時清除感染，而記憶CD8+?T細(xì)胞則提供長期保護(hù)，能夠在再次遇到相同抗原時迅速啟動免疫反應(yīng)。然而，單核細(xì)胞（monocytes）作為重要的髓系免疫細(xì)胞，雖在炎癥遷移中作用明確，但其是否直接參與T細(xì)胞分化尚無定論。

研究策略

動物模型與樣本制備

實(shí)驗(yàn)動物：優(yōu)先選擇8-12周齡的CD45.1+/CD45.2+及OT-I轉(zhuǎn)基因小鼠

樣本制備：脾臟單細(xì)胞懸液通過機(jī)械分離和紅細(xì)胞裂解（ACK緩沖液）制備

技術(shù)手段

• 流式細(xì)胞術(shù)分析

細(xì)胞分選：從感染LM-WT的小鼠脾臟中分選出cDCs、moDCs和單核細(xì)胞等抗原呈遞細(xì)胞（APCs）

細(xì)胞刺激與培養(yǎng)：將分選出的APCs與OVA肽段體外共孵育后，過繼轉(zhuǎn)移到已接受初始OT-I+?CD8+?T細(xì)胞的WT受體小鼠體內(nèi)，或與體外刺激的初始CD8+?T細(xì)胞共培養(yǎng)

表型檢測：通過流式細(xì)胞術(shù)檢測CD8+?T細(xì)胞的分化表型，包括記憶相關(guān)標(biāo)記物（如cKit、Sca1、Bcl6、TCF1和Eomes）的表達(dá)

• 空間分辨轉(zhuǎn)錄組學(xué)

樣本處理：對感染第0天和第5天的小鼠脾臟組織進(jìn)行切片，并進(jìn)行H&E染色以確定組織形態(tài)學(xué)特征

空間轉(zhuǎn)錄組學(xué)分析：利用BMKMANU?S1000空間轉(zhuǎn)錄組學(xué)技術(shù)對脾臟切片進(jìn)行分析，揭示不同免疫細(xì)胞亞群在脾臟中的空間分布

數(shù)據(jù)整合：將scRNA-seq數(shù)據(jù)與空間轉(zhuǎn)錄組學(xué)數(shù)據(jù)進(jìn)行整合，通過Seurat等工具的FindTransferAnchors和TransferData函數(shù)，將scRNA-seq中識別的細(xì)胞類型映射到空間轉(zhuǎn)錄組學(xué)數(shù)據(jù)中，實(shí)現(xiàn)細(xì)胞類型的空間定位

• 免疫熒光染色

樣本處理：對感染第0天和第5天的小鼠脾臟組織進(jìn)行切片，并進(jìn)行免疫熒光染色

抗體選擇：使用特異性抗體標(biāo)記CD8+?T細(xì)胞、單核細(xì)胞（CCR2+）、樹突狀細(xì)胞（CD11c+）等關(guān)鍵細(xì)胞類型

結(jié)果分析：通過顯微鏡觀察并拍攝染色切片，分析不同細(xì)胞類型在脾臟中的共定位關(guān)系

• 生物信息學(xué)分析

差異基因分析：利用Wilcoxon秩和檢驗(yàn)等統(tǒng)計方法，識別不同細(xì)胞亞群間的差異表達(dá)基因

通路富集分析：通過ClusterProfiler等工具對差異表達(dá)基因進(jìn)行通路富集分析，揭示潛在的生物學(xué)過程

偽時間分析：利用Monocle3等工具構(gòu)建細(xì)胞分化軌跡的偽時間軸，分析CD8+?T細(xì)胞亞群在分化過程中的基因表達(dá)變化

研究結(jié)果

研究結(jié)果一：通過空間分辨轉(zhuǎn)錄組學(xué)繪制急性感染后脾組織的圖譜

為了研究免疫反應(yīng)過程中不同免疫細(xì)胞的空間分布和CD8+?T細(xì)胞分化的調(diào)節(jié)機(jī)制，研究人員通過整合單細(xì)胞RNA測序（scRNA-seq）和BMKMANU?S1000空間轉(zhuǎn)錄組（ST）技術(shù)，分析了急性感染模型中脾臟免疫細(xì)胞的空間分布與CD8??T細(xì)胞分化機(jī)制。實(shí)驗(yàn)采用OT-I轉(zhuǎn)基因小鼠模型，通過李斯特菌（LM-OVA）感染誘導(dǎo)免疫反應(yīng)，關(guān)鍵發(fā)現(xiàn)包括：1）脾臟形成7個功能集群（B細(xì)胞、T細(xì)胞、單核/巨噬細(xì)胞等）；2）感染后T/B細(xì)胞區(qū)擴(kuò)大，APCs（B細(xì)胞、DC、巨噬細(xì)胞）在T細(xì)胞區(qū)外圍聚集；3）空間重組與T細(xì)胞分化（MP/TE亞群）顯著相關(guān)。多組學(xué)整合揭示了免疫微環(huán)境的動態(tài)調(diào)控機(jī)制，為理解感染中細(xì)胞互作提供了多組學(xué)視角。

研究結(jié)果二：急性感染期間效應(yīng)和memoryCD8?T細(xì)胞分化軌跡不同

為了探索CD8+?T細(xì)胞分化的空間決定因素，特別是與近端細(xì)胞的相互作用，對于急性感染模型研究人員進(jìn)行了亞群分析、分化路徑分析和功能驗(yàn)證，結(jié)果表明在急性感染中，IFN反應(yīng)型CD8+?T細(xì)胞和CD8+MP?細(xì)胞之間分化軌跡的不同意味著，效應(yīng)和記憶CD8+?T細(xì)胞的命運(yùn)決定于免疫反應(yīng)的初始階段，而IFN應(yīng)答和MP?CD8+?T細(xì)胞分別是效應(yīng)和記憶CD8+?T細(xì)胞分化途徑的初始階段。

研究結(jié)果三：跨組織區(qū)域細(xì)胞亞群的鑒定和空間圖譜

為了研究脾臟免疫細(xì)胞在感染前后的動態(tài)變化，研究人員通過整合單細(xì)胞RNA測序（scRNA-seq）和BMKMANU?S1000空間轉(zhuǎn)錄組（ST）技術(shù)，進(jìn)一步鑒定了供體和受體脾臟中的免疫細(xì)胞群。結(jié)果顯示：①細(xì)胞組成：從1.6萬細(xì)胞中鑒定出19類，包括CD8??T細(xì)胞亞群（如記憶前體和IFN反應(yīng)性細(xì)胞）及髓系細(xì)胞。②動態(tài)變化：感染后第5天，T細(xì)胞區(qū)從以初始CD8+?T細(xì)胞為主轉(zhuǎn)為記憶前體和效應(yīng)細(xì)胞主導(dǎo)，且單核細(xì)胞取代DC成為主要浸潤的髓系細(xì)胞。③功能提示：不同的APC和CD8+?T細(xì)胞之間的相互作用在免疫反應(yīng)和CD8+?T細(xì)胞的分化中起著重要作用。

研究結(jié)果四：單核細(xì)胞和CD8+?MP細(xì)胞的抗原依賴性共定位

為了探索不同細(xì)胞類型之間潛在的細(xì)胞相互作用，研究人員通過評估特異性免疫細(xì)胞標(biāo)記物的表達(dá)、免疫熒光染色及空間轉(zhuǎn)錄組分析，檢查了第0天和第5天脾臟中CD8+?T細(xì)胞和單核細(xì)胞的空間分布。結(jié)果顯示單核細(xì)胞和CD8+?MP細(xì)胞以抗原刺激依賴性的方式在空間上共定位，表明單核細(xì)胞可能對CD8+?MP細(xì)胞的分化有關(guān)鍵影響。

研究總結(jié)

CD8+T?細(xì)胞是適應(yīng)性免疫的重要執(zhí)行者；特別是記憶CD8+?T細(xì)胞對強(qiáng)效和長期保護(hù)至關(guān)重要。CD8+?T細(xì)胞分化在轉(zhuǎn)錄水平上的調(diào)控機(jī)制已被廣泛研究。然而，人們對感染過程中效應(yīng)和記憶CD8+?T細(xì)胞分化的空間要求仍然知之甚少。研究人員通過整合BMKMANU?S1000空間轉(zhuǎn)錄組（ST）技術(shù)和scRNA-seq技術(shù)，發(fā)現(xiàn)了記憶CD8+?T細(xì)胞分化的新機(jī)制，該機(jī)制涉及單核細(xì)胞和CD8+?MP細(xì)胞之間的解剖學(xué)鄰近性和TGF-β信號傳導(dǎo)。該研究結(jié)果闡述了記憶CD8+?T細(xì)胞命運(yùn)決定的新機(jī)制，并強(qiáng)調(diào)單核細(xì)胞作為關(guān)鍵的APC群體，通過細(xì)胞間接觸依賴的方式在感染過程中促進(jìn)記憶CD8+?T細(xì)胞的分化。

英文題目：Identifification of epilepsy-associated neuronal subtypes and gene expression underlying epileptogenesis

中文題目：鑒定癲癇相關(guān)的神經(jīng)元亞型和誘發(fā)癲癇的基因表達(dá)

發(fā)表時間：2020年10月7號

影響因子：12.121

研究背景

癲癇是一種神經(jīng)系統(tǒng)疾病，具有自發(fā)性和反復(fù)發(fā)作的特征，其主要發(fā)生在海馬體或大腦皮層區(qū)域。人們對癲癇的病理生理學(xué)仍然知之甚少。雖然有一些動物模型研究顯示某些神經(jīng)元亞型對癲癇發(fā)作的產(chǎn)生和傳播有一定的促進(jìn)作用，但在人類癲癇患者的相應(yīng)數(shù)據(jù)卻很少。這是由于癲癇發(fā)生過程中復(fù)雜的神經(jīng)網(wǎng)絡(luò)現(xiàn)象。

實(shí)驗(yàn)方法

10x snRNA-seq，Smart-seq2?snRNA-seq，RNA-seq，ISH等

結(jié)論

該項(xiàng)研究揭示了癲癇對神經(jīng)元轉(zhuǎn)錄的不同影響—即雖然許多亞型表現(xiàn)出輕微的基因表達(dá)變化，但其對一些特定的SUB型主要神經(jīng)元和GABA能間神經(jīng)元產(chǎn)生了很大的影響。

結(jié)果

一、癲癇皮質(zhì)中與疾病相關(guān)的神經(jīng)元亞型

為了識別哪些神經(jīng)元亞型受癲癇影響或?qū)е掳d癇，文章中比較了來自癲癇性顳皮質(zhì)和非癲癇性顳皮質(zhì)的snRNASeq數(shù)據(jù)，研究結(jié)果表明，在一些位置，例如，皮層上層L2_Cux2_Lamp5和L2-3_Cux2_Frem3，在癲癇和非癲癇的皮層產(chǎn)生的神經(jīng)元之間存在明顯的轉(zhuǎn)錄組轉(zhuǎn)移(圖1a)。此外，盡管癲癇病和非癲癇病樣本中每個亞型的神經(jīng)元數(shù)量大致相似，但觀察到幾種亞型中已被鑒定的細(xì)胞核數(shù)量顯著減少(圖1b)。在癲癇患者中，L2/3亞型的數(shù)量減少了，當(dāng)對每種情況下的序列神經(jīng)元總數(shù)進(jìn)行歸一化時，這種減少更加明顯(圖1c)。對于中間神經(jīng)元，Pvalb_Sulf1亞型的神經(jīng)元數(shù)量減少較多(圖1c)。

通過基因表達(dá)相關(guān)性分析，發(fā)現(xiàn)癲癇和非癲癇神經(jīng)亞型之間的巨大差異(圖1d)，這可能表明癲癇對其轉(zhuǎn)錄組有特定的影響。對于主要神經(jīng)元，L5-6_Fezf2_Tle4_Abo、L5-6_Themis_Ntng2、L2_
Cux2_Lamp5和L2-3_Cux2_Frem3亞型在癲癇和非癲癇皮質(zhì)中轉(zhuǎn)錄組差異大(圖1d)。對于gaba能間神經(jīng)元，在癲癇發(fā)作的顳葉皮層中，Vip_Cbln1、Sst_Tac1、Pvalb_Sulf1、Pvalb_Nos1和Id2_Lamp5_Nos1的轉(zhuǎn)錄組改變大。同時也證實(shí)了癲癇相關(guān)DE基因在具有較大轉(zhuǎn)錄組差異的神經(jīng)元亞型中普遍存在(圖1e、f)。

圖1癲癇和非癲癇數(shù)據(jù)集的整合和疾病相關(guān)神經(jīng)元亞型的鑒定

二、共有和亞型特異性癲癇相關(guān)途徑

作者計算了每種已鑒定亞型GO中DE基因的富集程度，一些神經(jīng)元亞型在癲癇中表現(xiàn)出了大的轉(zhuǎn)錄組改變(>100豐富的GO)，但多數(shù)神經(jīng)元亞型中只有少數(shù)GO，特別是對于生物過程(BP)來說，這應(yīng)該與通路的生物學(xué)功能相關(guān)(圖2a)。這些受影響較小的亞型中，大多數(shù)在癲癇細(xì)胞和非癲癇細(xì)胞之間也具有高的基因表達(dá)相關(guān)性(圖1d)，從而證實(shí)癲癇和非癲癇之間轉(zhuǎn)錄組相關(guān)性高的亞型中，癲癇改變的信號通路也較少。

根據(jù)GO分析，發(fā)現(xiàn)一些Sst亞型、Vip_Cbln1和Id2_Lamp5亞型顯示DO高度富集，如局灶性癲癇(DOID:2234)、癲癇綜合征(DOID:1826)和顳葉癲癇(DOID:3328)。此外，還發(fā)在L3_Cux2_Prss12、L5-6_Fezf2_Lrrk1_Sema3e和L5-6_Fezf2_Tle4_Abo主要神經(jīng)元亞型中富集的與癲癇相關(guān)的DO相同。這表明，與癲癇相關(guān)的DO在不同的神經(jīng)元亞型中富集差異很大，與谷氨酰胺能細(xì)胞相比，在GABAergic和Sst亞型Vip_Cbln1和Id2_Lamp5中，癲癇相關(guān)的DO在GABAergic和Sst亞型Vip_Cbln1和Id2_Lamp5中富集更廣泛，顯示出對癲癇的特殊易感性。

圖2神經(jīng)元亞型中癲癇相關(guān)通路和轉(zhuǎn)錄轉(zhuǎn)移的鑒定

三、癲癇活動的信號通路

谷氨酸受體，尤其是AMPA受體被證明是癲癇發(fā)作的主要驅(qū)動力之一，在篩選出缺乏表達(dá)的基因后，尋找富含與奮性相關(guān)轉(zhuǎn)錄組學(xué)變化的亞型，并從這些GO項(xiàng)中繪制DE基因（圖3a–d）。CKAMP44是貫穿整個顳葉皮質(zhì)的主要AMPA受體輔助亞基，這一事實(shí)也強(qiáng)調(diào)了CKAMP44上調(diào)的重要性(圖3c, d)。癲癇對AMPA受體輔助亞基編碼基因的上調(diào)有相當(dāng)普遍的影響。

多種谷氨酸受體亞基和神經(jīng)元活性相關(guān)基因表達(dá)復(fù)雜失調(diào)，其中大部分在癲癇中未見報道。因此,雖然大多數(shù)的基因編碼谷氨酸受體亞基表達(dá)上調(diào)(GRIA1、GRIA3 GRIA4, GRIK3, GRIK4, GRIK5, GRIN2B, GRIN3A, GRM1, GRM7,和GRM8),但一些是表達(dá)下調(diào)(GRIA2、GRIA3 GRIN2A, GRM5, GRIK1,和GRIK2)(圖3?b, c)。

圖3鑒別可能導(dǎo)致癲癇發(fā)作的皮層神經(jīng)元亞型的信號通路和基因

癲癇患者皮層谷氨酸介導(dǎo)興奮相關(guān)基因顯著分層上調(diào)或下調(diào)。作者用單分子熒光原位雜交(smFISH)方法標(biāo)記了幾個高度調(diào)控基因的mRNA。如上所述，在主要神經(jīng)元的亞型中，編碼AMPA輔助亞基的許多基因分層上調(diào)。因此，在一組癲癇和非癲癇樣本的皮質(zhì)切片中標(biāo)記了CKAMP44 mRNA，通過與Rorb和DAPI共標(biāo)記，并確定了L2/3、L4和L5/6的位置(圖4b)。重要的是，證實(shí)了各層中CKAMP44的表達(dá)顯著上調(diào)(圖4b)。

圖4癲癇患者皮層谷氨酸介導(dǎo)興奮相關(guān)基因的復(fù)雜分層失調(diào)

四、受癲癇影響的基因網(wǎng)絡(luò)

在癲癇樣本中，六個與癲癇密切相關(guān)的模塊上調(diào)，六個下調(diào)。另外，在細(xì)胞水平分析中突出的幾個涉及過度興奮信號和癲癇發(fā)作的基因也是一個或多個癲癇相關(guān)基因模塊的成員。編碼谷氨酸受體亞基的基因也廣泛存在于癲癇相關(guān)模塊中。這些基因共表達(dá)分析捕獲了主要神經(jīng)元L5-6_Fezf2_Tle4_Abo和L2-3_Cux2亞型中活躍的核心轉(zhuǎn)錄網(wǎng)絡(luò)。與癲癇狀態(tài)相關(guān)的基因模塊指向與突觸和離子通道相關(guān)的轉(zhuǎn)錄變化，包括AMPA受體輔助亞基、谷氨酸受體亞基和電壓門控鈉通道，這些可能與癲癇回路異常高的興奮有關(guān)。同時觀察到主要神經(jīng)元和GABAergic間神經(jīng)元的受影響的亞型根據(jù)豐富的GO項(xiàng)的相似性聚集在一起(圖5b)。因此，L2-3_Cux2和L5-6_Fezf2亞型與Sst_Tac1和Vip_Cbln1亞型共聚，而L3_Cux2_Prss12亞型與Pvalb_Sulf1亞型共聚，這可能是癲癇轉(zhuǎn)錄組中受影響嚴(yán)重的局部神經(jīng)網(wǎng)絡(luò)的基礎(chǔ)。

圖5綜合分析癲癇易受影響的神經(jīng)元亞型

結(jié)論

總之，該研究通過單核轉(zhuǎn)錄組測序（snRNA-seq）發(fā)現(xiàn)癲癇患者的神經(jīng)轉(zhuǎn)錄組發(fā)生了大規(guī)模而復(fù)雜的變化，其中一些亞型表現(xiàn)出了明顯的癲癇驅(qū)動的基因表達(dá)異常，而其他亞型則大部分未受影響。癲癇相關(guān)的轉(zhuǎn)錄組改變可以聚集成包含多個神經(jīng)元亞型的模塊，這些神經(jīng)元亞型可能是受癲癇影響的不同神經(jīng)元聚集的基礎(chǔ)。未來需要在小鼠模型和人類組織中進(jìn)行轉(zhuǎn)化研究，以解決哪些已確定的通路導(dǎo)致癲癇的產(chǎn)生和傳播，以及哪些更能代表神經(jīng)網(wǎng)絡(luò)的穩(wěn)態(tài)可塑性。

百邁客引進(jìn)10xGenomics單細(xì)胞測序平臺，使用Chromium系統(tǒng)采用微流控、油滴包裹和barcode標(biāo)記等技術(shù)實(shí)現(xiàn)一次性分離、高效標(biāo)記捕獲；同時具有10x 單細(xì)胞轉(zhuǎn)錄組、單細(xì)胞核轉(zhuǎn)錄組、單細(xì)胞ATAC-seq、單細(xì)胞免疫組庫、全長轉(zhuǎn)錄組測序、空間轉(zhuǎn)錄組，實(shí)現(xiàn)10x平臺全面優(yōu)質(zhì)服務(wù)；已經(jīng)具有大量單細(xì)胞分離捕獲，極低量RNA反轉(zhuǎn)錄擴(kuò)增建庫成功經(jīng)驗(yàn)；提供單細(xì)胞分離捕獲、反轉(zhuǎn)錄建庫、測序、標(biāo)準(zhǔn)分析和高級分析全套單細(xì)胞測序服務(wù)；強(qiáng)大的生信團(tuán)隊(duì)不僅提供基本分析，還提供細(xì)胞分化軌跡分析等多種高級分析；資深單細(xì)胞技術(shù)人員為您提供專業(yè)的課題方案設(shè)計，為您量身訂造專屬個性化分析。

點(diǎn)擊下方按鈕聯(lián)系我們，免費(fèi)獲取文章思路設(shè)計方案。

本次我們給大家?guī)韮善脝渭?xì)胞轉(zhuǎn)錄組技術(shù)助力動物細(xì)胞圖譜繪制的高分文章解讀，希望能給老師后續(xù)的研究提供思路。

動物細(xì)胞圖譜分析繪制案例：黑腹果蠅

中文標(biāo)題：果蠅細(xì)胞圖譜：成年果蠅的單核轉(zhuǎn)錄組圖譜

英文標(biāo)題：Fly Cell Atlas: A single-nucleus transcriptomic atlas of the adult fruit fly

期刊：Science[IF: 63.714] 2022.3.4

DOI：10.1126/science.abk2432

實(shí)驗(yàn)設(shè)計

實(shí)驗(yàn)材料：黑腹果蠅

實(shí)驗(yàn)方法：解剖了來自雌性和雄性黑腹果蠅的12個組織（觸角、體壁、脂肪體、平衡棒、心臟、腸道、腿、馬氏管、絳色細(xì)胞、喙和下顎須、氣管、翅膀）以及3個性別特異性組織（雄性生殖腺、睪丸、卵巢）。對于遍布全身的組織使用特異性GAL4驅(qū)動核-GFP蛋白，再用流式細(xì)胞熒光分選技術(shù)標(biāo)記和收集細(xì)胞核，進(jìn)行單細(xì)胞核snRNA-seq

測序策略：單細(xì)胞核轉(zhuǎn)錄組10x Genomics、Smart-seq2

研究內(nèi)容

黑腹果蠅在生物學(xué)研究方面有著豐富的歷史，以往研究探索了其不同組織中的表達(dá)模式，但缺乏細(xì)胞類型分辨率級別的數(shù)據(jù)庫。近年來單細(xì)胞技術(shù)的進(jìn)步使得同時對數(shù)千個細(xì)胞進(jìn)行轉(zhuǎn)錄組分析成為可能，促進(jìn)了全組織細(xì)胞圖譜的創(chuàng)建。然而，現(xiàn)有scRNA-seq數(shù)據(jù)集來自不同實(shí)驗(yàn)室，是基于不同遺傳背景，不同分離方法以及不同測序平臺產(chǎn)生的，阻礙了跨細(xì)胞和組織的基因表達(dá)的系統(tǒng)性分析。

該研究通過對果蠅多組織進(jìn)行單細(xì)胞核轉(zhuǎn)錄組測序后，得到了果蠅細(xì)胞圖譜。之后對15個組織進(jìn)行詳細(xì)的細(xì)胞注釋，在果蠅頭部，注釋到了81種主要的神經(jīng)細(xì)胞類型，在果蠅的軀體中注釋了最豐富的33種細(xì)胞類型。該圖譜還可用于進(jìn)行相同細(xì)胞的跨組織分析，作者通過分析血細(xì)胞在整個組織中的表達(dá)情況，發(fā)現(xiàn)了血細(xì)胞中最常見的細(xì)胞類型——漿細(xì)胞，并未在成蟲血細(xì)胞中觀察到葉狀血細(xì)胞。通過對肌肉細(xì)胞在不同組織類型進(jìn)行比較，發(fā)現(xiàn)它們主要在果蠅軀體、體壁和腿部有特定的富集，并發(fā)現(xiàn)內(nèi)臟肌、骨骼肌、間接飛行肌的分離。

通過轉(zhuǎn)錄因子分析，找到了500個轉(zhuǎn)錄因子和細(xì)胞類型具有高特異性，繪制了每種細(xì)胞類型中特異轉(zhuǎn)錄因子的表達(dá)熱圖。此外，作者還使用SCENIC預(yù)測了基于共表達(dá)和motif富集的基因調(diào)控網(wǎng)絡(luò)。通過對廣泛細(xì)胞類型或組織的基因進(jìn)行比較分析，發(fā)現(xiàn)了在各種細(xì)胞類型中的常見表達(dá)基因和特定表達(dá)基因。

通過對性別依賴的基因表達(dá)和性別特異性組織分析，發(fā)現(xiàn)體細(xì)胞中的主要性別決定基因doublesex（dsx）的表達(dá)在很大程度上不具有性別特異性，而許多其他基因的表達(dá)具有性別依賴。作者進(jìn)一步對精母細(xì)胞和精細(xì)胞進(jìn)行軌跡推斷，發(fā)現(xiàn)精母細(xì)胞期被轉(zhuǎn)錄的基因數(shù)量持續(xù)增加，許多強(qiáng)烈上調(diào)的基因在任何其他細(xì)胞類型中都沒有基本表達(dá)。然而，晚期精母細(xì)胞顯示了來自許多其他細(xì)胞類型的標(biāo)記基因的表達(dá)。

果蠅細(xì)胞圖譜：成年果蠅的單核轉(zhuǎn)錄組圖譜

動物細(xì)胞圖譜分析繪制案例：食蟹猴

中文標(biāo)題：食蟹猴單細(xì)胞轉(zhuǎn)錄組和調(diào)控組參考圖譜

英文標(biāo)題：A reference single-cell regulomic and transcriptomic map of cynomolgus monkeys

期刊：Nature comuiations?[IF: 17.694]?2022.07.13

DOI：10.1038/s41467-022-31770-x

實(shí)驗(yàn)設(shè)計

實(shí)驗(yàn)材料：成年食蟹猴

實(shí)驗(yàn)方法：選取16個器官（心、肝、脾、肺、腎、胃、結(jié)腸、肌肉、氣管、主動脈、脂肪、膀胱、舌頭、乳腺、子宮和睪丸）

測序策略：10x Genomics單細(xì)胞轉(zhuǎn)錄組測序（scRNA-seq）和單細(xì)胞染色質(zhì)開放性測序（scATAC-seq）

研究內(nèi)容

非人類靈長類動物（Non-human primates）在系統(tǒng)發(fā)育上與人類非常接近，在遺傳、器官發(fā)育、生理功能、病理反應(yīng)和生化代謝等諸多方面表現(xiàn)出與人類相似的特征，是生物醫(yī)學(xué)研究和藥物開發(fā)的理想實(shí)驗(yàn)動物模型。其中，最為典型的代表當(dāng)屬食蟹猴。系統(tǒng)評估食蟹猴等非人靈長類動物模型與人類的細(xì)胞組成差異、器官異質(zhì)性和基因表達(dá)時空特異性等在基礎(chǔ)研究中具有十分重要的價值。

該研究對成年食蟹猴的16個代表性器官進(jìn)行單細(xì)胞轉(zhuǎn)錄組測序和單細(xì)胞染色質(zhì)開放性測序，構(gòu)建了食蟹猴多器官的單細(xì)胞多組學(xué)參考圖譜，將25萬個細(xì)胞進(jìn)行 t-SNE降維聚類得到40余種不同的細(xì)胞亞群。

通過對該細(xì)胞圖譜進(jìn)行分析，鑒定到了新的細(xì)胞類型。為了解析上皮細(xì)胞異質(zhì)性，作者提取上皮細(xì)胞并進(jìn)行亞群聚類分析。根據(jù)marker基因的獨(dú)特表達(dá)模式，分析鑒定出14個上皮細(xì)胞簇，包括基底細(xì)胞、分泌細(xì)胞、纖毛細(xì)胞和非纖毛細(xì)胞。纖毛上皮細(xì)胞在各種組織中都占很大比例，為了探索纖毛上皮細(xì)胞亞型的發(fā)育和功能動態(tài)，作者選用纖毛上皮細(xì)胞相對較多的器官，使用Monocle和RNA速度分析對纖毛上皮細(xì)胞進(jìn)行軌跡分析，觀察到纖毛細(xì)胞從祖細(xì)胞狀態(tài)向成熟狀態(tài)分化的過程，沿著擬時間高表達(dá)的基因在與代謝過程、細(xì)胞對刺激的反應(yīng)和防御反應(yīng)相關(guān)的基因本體(GO)術(shù)語中順序富集。

使用CellPhoneDB細(xì)胞通訊分析，發(fā)現(xiàn)基質(zhì)細(xì)胞、上皮細(xì)胞和髓系細(xì)胞之間存在強(qiáng)烈的細(xì)胞間相互作用。通常，細(xì)胞間相互作用的強(qiáng)度和模式是器官特異性的。為了重新繪制調(diào)節(jié)細(xì)胞-細(xì)胞相互作用的分子相互作用，作者在不同器官的特定細(xì)胞亞群中繪制了配體-受體對，揭示了猴子各種器官細(xì)胞間通信的潛在分子機(jī)制。

結(jié)合表達(dá)譜和染色質(zhì)開放性數(shù)據(jù)，作者分析預(yù)測了控制不同細(xì)胞類型基因表達(dá)模式的關(guān)鍵調(diào)控子。發(fā)現(xiàn)SPIB、POU2F2、SPI1、CEBPD和IRF4是髓系細(xì)胞中的關(guān)鍵調(diào)控子；而FEV具有調(diào)節(jié)造血干細(xì)胞的潛能，該轉(zhuǎn)錄因子廣泛參與對免疫細(xì)胞和上皮細(xì)胞的基因表達(dá)調(diào)控。大多數(shù)順式調(diào)控元件(CREs或ATAC峰)來自啟動子、內(nèi)含子或遠(yuǎn)端基因間調(diào)控區(qū)。在scATAC-seq數(shù)據(jù)中預(yù)測了9種細(xì)胞類型，在RNA簇中發(fā)現(xiàn)了兩種罕見的細(xì)胞類型，單核細(xì)胞和循環(huán)B細(xì)胞，但在ATAC簇中沒有發(fā)現(xiàn)。總之，本文中的scATAC-seq數(shù)據(jù)為無偏見地發(fā)現(xiàn)食蟹猴的細(xì)胞類型和調(diào)控DNA元件提供了豐富的資源。

最后，作者進(jìn)行了跨物種比較分析，通過整合分析器官匹配的單細(xì)胞轉(zhuǎn)錄組數(shù)據(jù)，探討了人、食蟹猴和小鼠三個物種之間的細(xì)胞組成和基因表達(dá)特異性。猴子和人類在同源基因表達(dá)方面表現(xiàn)出明顯高于其他比較的細(xì)胞類型相似性。其中，免疫細(xì)胞相較于非免疫細(xì)胞，人類和猴子之間的基因表達(dá)具有更高的相似性。間質(zhì)細(xì)胞在人鼠和猴鼠比較中表現(xiàn)出最高的相似性。這些發(fā)現(xiàn)表明，猴子在免疫系統(tǒng)中與人類具有高度相似的轉(zhuǎn)錄程序，因此可能為研究對癌癥或新冠病毒COVID-19等疾病的免疫反應(yīng)提供理想的模型。

食蟹猴多器官單細(xì)胞多組學(xué)參考圖譜構(gòu)建與分析

如果您對單細(xì)胞轉(zhuǎn)錄組測序技術(shù)感興趣，歡迎點(diǎn)擊下方按鈕聯(lián)系我們，我們將免費(fèi)為您設(shè)計文章思路方案。