單細胞RNA/蛋白組測序鑒定小鼠胰島細胞亞群對高脂飲食代謝適應性調(diào)控

圖1

2. Ace2在α細胞中的異質(zhì)表達
為了分析α細胞的異質(zhì)性，使用scRNA-seq數(shù)據(jù)將α細胞聚成8個簇（圖2a-b），并分析了8個聚類中的marker基因表達（圖2c）。本研究發(fā)現(xiàn)與嚴重急性呼吸綜合征冠狀病毒2 (SARS-CoV-2)感染相關的關鍵基因Ace2在α細胞群中顯著變化。根據(jù)Ace2的表達水平，將8個α細胞簇分為Ace2^low(簇3和5)和Ace2^high(簇0、1、2、4、6和7)亞群。ScRNA-seq和免疫熒光證實了在cd和hfd小鼠胰島中α細胞亞群的存在（圖f,g,j）。為了探索ace2 α細胞異質(zhì)性的分子特征，我們在Ace2^low和Ace2^high亞群之間鑒定出379個差異基因（圖2h），并對差異基因進行了KEGG富集分析（圖2i），表明Ace2的表達可以區(qū)分具有不同功能和生存狀態(tài)的α細胞。我們還發(fā)現(xiàn)，在Ace2^low組中，α細胞特異性轉錄因子（Mafb和Irx1）和α細胞相關基因（Ttr、Chgb、Fev和Pyy）顯著下調(diào)。相反，Ace2^low α細胞開始表達具有β細胞特征的基因（如Ucn3、Ero1lb、和Slc2a2）。其中一些結果通過qPCR（圖2m-q）和western印跡（圖2k）對分選的α細胞群進行了證實。這表明，Ace2^low亞群可能代表了未成熟的α細胞，似乎具有β細胞樣特征。
 

圖2

3. 表明α細胞亞群功能的單細胞蛋白質(zhì)組學
用胰島單細胞懸液的 FACS將α細胞（胰高血糖素[GCG]+）細胞分為ACE2^high組和ACE2^low組。值得注意的是，HFD小鼠中ACE2^high α細胞的比例（65.3%）明顯低于CD小鼠（80.0%）,這表明基于ace2的α細胞異質(zhì)性可能參與了糖尿病的病理生理過程（圖3a,b）。從FACS分選的的CD和HFD小鼠（每個樣本2000個細胞）中收集ACE2^low和ACE2^high α細胞亞群進行單細胞蛋白質(zhì)組學分析。檢測到的蛋白數(shù)量在3042~4021之間，重復性良好。PCA也能明顯區(qū)分ACE2^low細胞和ACE2^high細胞（圖3c）。在CD小鼠中，與ACE2^high α細胞相比，ACE2^low α細胞中有29個上調(diào)，80個DEPs下調(diào)(圖3d)，受影響的KEGG通路包括Ras信號通路和有絲分裂活化蛋白激酶（MAPK）信號通路（圖3 e）。在HFD小鼠中，與ACE2^high α細胞相比，ACE2^low α細胞中的DEPs明顯增加，其中17個上調(diào)，301個DEPs下調(diào)（圖3 f）。KEGG分析發(fā)現(xiàn)與營養(yǎng)和能量代謝相關的途徑在ACE2低組均顯著下調(diào)（圖3g）。這些發(fā)現(xiàn)強烈暗示了ACE2^low α細胞的代謝和其他生物活性減弱，特別是在HFD小鼠中。
 

圖3 

4. CD81在β細胞中的異質(zhì)表達
為了分析β細胞的異質(zhì)性，我們將這些細胞聚為9個簇，并確定了前10個標記基因（圖4a-c）。根據(jù)最近的一項研究，發(fā)現(xiàn)Cd81在細胞簇中有差異表達。研究將9個β細胞簇分為Cd81^low（簇0-4和7）和Cd81^high（簇5、6和8）亞群（4d-g），發(fā)現(xiàn)CD81在CD和HFD小鼠的β細胞中分布的明顯異質(zhì)性，同時也經(jīng)過了免疫組化的證實（圖4j）。接下來分析了Cd81^low相對于Cd81 ^highβ細胞的差異基因和KEGG通路分析，細胞代謝及其相關功能通路以及一系列應激相關和炎癥通路被顯著富集。而未成熟的細胞標記物（Chga、Chgb、Mafb和Rbp4）在Cd81^high亞群中高表達（圖4l）。這表明Cd81^high亞群可能代表了一組功能較弱和未成熟的細胞。其中一些結果通過qPCR（圖4m-p）和westernblot（圖4k）對分選的細胞亞群進行證實。
 

圖4

5. β細胞亞群功能的單細胞蛋白質(zhì)組學鑒定

為了進一步探索HFD對β細胞亞群的影響，F(xiàn)ACS顯示與CD小鼠相比，HFD小鼠的CD81^high細胞比例降低。然后，使用單細胞蛋白質(zhì)組學（每個樣本10000個細胞）分析FACS分選的CD81^low和CD81^highβ細胞（圖5a,b）。PCA顯示CD81^low和CD81^highβ細胞被聚集成兩個獨立的組（圖5c）。比較CD81^low和CD81^highβ細胞，在CD小鼠中檢測到869個上調(diào)和373個下調(diào)的DEPs，在HFD小鼠中檢測到707個上調(diào)和312個下調(diào)的DEPs（圖5d,e）。在CD和HFD小鼠中，與能量代謝相關的途徑（包括氧化磷酸化和阿爾茨海默�。┰贑D81低亞群中被下調(diào)（圖5f,g），并通過免疫組化和western blot進行了驗證。
 

圖5

6. HFD增強了成熟的Slc2a2 ^low δ細胞的功能

本研究將δ細胞分為 Slc2a2^low and Slc2a2^high 亞群，通過Slc2a2^low 相比于 Slc2a2^high差異基因和KEGG功能富集分析顯示，上調(diào)的通路與鈣離子和激素分泌有關，提示Slc2a2^low亞組可能功能更成熟（S6a,b）。然而，8周的HFD并沒有改變兩個delta細胞亞群的比例（圖6a,c）。免疫熒光證實了在小鼠和人類胰島中均存在δ細胞亞群（圖6b）。接下來采用FASC對CD和HFD小鼠的Slc2a2^low and Slc2a2^high  δ細胞亞群（每個樣本3000個細胞）進行單細胞蛋白質(zhì)組學分析。PCA數(shù)據(jù)顯示，Slc2a2^low 和Slc2a2^high  δ細胞聚集成兩個不同的群體（圖6d）。Slc2a2^low 和 Slc2a2^high    δ細胞亞群之間有548個上調(diào)和171個下調(diào)的DEPs；同時，HFD小鼠中DEPs上調(diào)較多（832），下調(diào)較少（102）（圖6e,f)。KEGG分析顯示，與β細胞相似，CD小鼠的Slc2a2^low （相對于Slc2a2^high  δ細胞）中的氧化磷酸化和阿爾茨海默病通路顯著下調(diào)(圖6g)。在CD和HFD小鼠的Slc2a2^low細胞中，多種營養(yǎng)代謝途徑也顯著上調(diào)。綜上所述，Slc2a2^low δ細胞可能比Slc2a2^high δ更具有功能和代謝活性，而HFD進一步刺激Slc2a2^low δ細胞的功能和代謝。
 

圖6

總之，我們發(fā)現(xiàn)未成熟的胰島細胞部分失去了它們的身份維持基因，并開始表現(xiàn)出其他內(nèi)分泌細胞的特征。在HFD誘導的葡萄糖耐受不良中，低容量和未成熟的ACE2^low細胞的比例隨著功能活性和成熟的CD81^low細胞比例的增加而增加。成熟的Slc2a2^low δ細胞功能增強，但比例不變。這些動態(tài)變化揭示了胰島對HFD應激反應的代謝可塑性。這種可塑性適應可以通過抑制高血糖來抵消代謝應激，恢復葡萄糖穩(wěn)態(tài)，增強細胞的低血糖作用，并增加δ對α和β細胞的調(diào)節(jié)。
 

圖7
 

參考文獻：
Fu Q, Jiang H, Qian Y, et al. Single-cell RNA sequencing combined with single-cell proteomics identifies the metabolic adaptation of islet cell subpopulations to high-fat diet in mice[J]. Diabetologia, 2022: 1-17.