哺乳動(dòng)物生長發(fā)育過程一直是一個(gè)經(jīng)典的話題，隨著單細(xì)胞技術(shù)發(fā)展，單細(xì)胞轉(zhuǎn)錄組研究揭示哺乳動(dòng)物大腦皮層細(xì)胞的異質(zhì)性，但是對于表觀調(diào)控研究還未深入。為了在單細(xì)胞水平研究轉(zhuǎn)錄組和表觀調(diào)控的關(guān)系，百邁客可以提供單細(xì)胞多組學(xué)ATAC&基因表達(dá)檢測（scATAC&GEX-seq），真實(shí)反映同一個(gè)細(xì)胞內(nèi)表觀調(diào)控與基因表達(dá)內(nèi)在聯(lián)系，確定假定的調(diào)控元件（啟動(dòng)子，調(diào)節(jié)因子）和它們的靶基因之間的聯(lián)系。此外百邁客業(yè)務(wù)還包括單細(xì)胞轉(zhuǎn)錄組、空間轉(zhuǎn)錄組、、單細(xì)胞免疫組庫等，有興趣的老師可以聯(lián)系當(dāng)?shù)劁N售。接下來我們看看Nature發(fā)表單細(xì)胞水平表觀基因組學(xué)研究思路，文章詳情如下：

發(fā)表期刊：Nature

影響因子：49.962

原文鏈接：https://international.biocloud.net/zh/article/detail/34616060

發(fā)表日期：2021-10-06

 

摘要

在哺乳動(dòng)物發(fā)育過程中，染色質(zhì)狀態(tài)的差異與細(xì)胞分化相一致，反映了基因調(diào)控圖譜的變化。在大腦發(fā)育過程中，細(xì)胞命運(yùn)和圖譜特征對于定義細(xì)胞特征和賦予神經(jīng)發(fā)育障礙選擇脆弱性非常重要。為了確定發(fā)育中人腦的特定細(xì)胞類型的染色質(zhì)可及性模式，作者使用了一種單細(xì)胞轉(zhuǎn)座酶可及性測序分析方法(scATAC-seq)，對來自人前腦的原始組織樣本進(jìn)行了測序分析。應(yīng)用無偏估計(jì)來識別在神經(jīng)形成過程中經(jīng)歷了廣泛的細(xì)胞類型和大腦區(qū)域特異性可及性變化的基因組位點(diǎn)，并應(yīng)用綜合分析來預(yù)測特定細(xì)胞類型的候選調(diào)控元件。結(jié)果發(fā)現(xiàn)，腦類器官包含了大多數(shù)假定的細(xì)胞類型特異性增強(qiáng)子可及性模式，但缺乏許多體內(nèi)發(fā)現(xiàn)的細(xì)胞類型特異性開放染色質(zhì)區(qū)域。對腦區(qū)域染色質(zhì)可及性進(jìn)行系統(tǒng)比較，發(fā)現(xiàn)了大腦皮層神經(jīng)前體細(xì)胞之間出人意料的差異，表明視黃酸信號在前額葉皮質(zhì)神經(jīng)元譜系確定中起著重要作用?？傊?，結(jié)果揭示了染色質(zhì)狀態(tài)對細(xì)胞類型多樣性和細(xì)胞命運(yùn)特異性的發(fā)展模式具有重要作用，并評估了腦類器官作為皮層發(fā)育模型的√確度和穩(wěn)定性。

 

背景介紹

大腦皮層的細(xì)胞類型通常是根據(jù)少數(shù)形態(tài)學(xué)、解剖學(xué)和生理學(xué)特征進(jìn)行分類的(圖1)。新興的單細(xì)胞轉(zhuǎn)錄學(xué)，包括single-cell RNA sequencing (scRNA-seq)，已經(jīng)能夠大規(guī)模平行分析單個(gè)細(xì)胞的數(shù)千個(gè)分子特征，并揭示了密切相關(guān)的細(xì)胞類型之間的差異，例如位于大腦皮層不同區(qū)域的興奮性神經(jīng)元。盡管取得了一些進(jìn)展，但人大腦皮層中出現(xiàn)不同神經(jīng)元譜系的發(fā)育機(jī)制仍然基本未知。
染色質(zhì)狀態(tài)通過調(diào)節(jié)基因調(diào)控元件(如增強(qiáng)子)的可及性來確定基因組的功能結(jié)構(gòu)，這些元件是轉(zhuǎn)錄調(diào)節(jié)因子的結(jié)合位點(diǎn)。在發(fā)育過程中，轉(zhuǎn)錄因子的順序級聯(lián)逐漸重塑，并且不同細(xì)胞類型呈現(xiàn)出染色質(zhì)可及性的差異化。盡管染色質(zhì)的可及性是細(xì)胞特性的一個(gè)基本特征，但關(guān)于大腦發(fā)育過程中染色質(zhì)狀態(tài)變化的研究相對較少。目前單細(xì)胞基因組學(xué)已經(jīng)能夠使用scATAC-seq以細(xì)胞分辨率對染色質(zhì)狀態(tài)進(jìn)行可伸縮分析。在發(fā)育中的小鼠大腦中，scATAC-seq揭示了作為神經(jīng)發(fā)育過程基礎(chǔ)的染色質(zhì)可及性的高度動(dòng)態(tài)變化。有必要將這些研究擴(kuò)展到人類初級組織，以便更好地了解非編碼調(diào)控元件(包括人類特有的神經(jīng)發(fā)育增強(qiáng)子)的突變?nèi)绾胃蓴_正常的發(fā)育過程，并導(dǎo)致精神神經(jīng)發(fā)育障礙的遺傳基礎(chǔ)。

圖1：發(fā)育中皮質(zhì)的橫截面示意圖，顯示了主要細(xì)胞類型。

發(fā)育中大腦的染色質(zhì)狀態(tài)

為了在單細(xì)胞分辨率下觀察發(fā)育中人腦的染色質(zhì)狀態(tài)，作者對6例妊娠中期的人前腦進(jìn)行了scATAC-seq，包括背外側(cè)前額葉皮層(PFC)、初級視覺皮質(zhì)(V1)、初級運(yùn)動(dòng)皮質(zhì)(M1)、初級軀體感覺皮層、背外側(cè)頂葉皮質(zhì)、顳葉皮質(zhì)、島葉皮質(zhì)和內(nèi)側(cè)神經(jīng)節(jié)隆起(MGE)。共有77,354個(gè)細(xì)胞通過質(zhì)控進(jìn)行下游分析。scATAC數(shù)據(jù)庫的聚類結(jié)果與并行產(chǎn)生的混合ATAC-seq測序的數(shù)據(jù)高度相關(guān)，并且生物學(xué)重復(fù)之間的相關(guān)性系數(shù)也很高。對數(shù)據(jù)進(jìn)行了降維處理與批次效應(yīng)校正。使用Leiden群落檢測算法識別了25個(gè)不同的clusters。有力地鑒定了皮質(zhì)和皮質(zhì)下(MGE)兩種不同的細(xì)胞類型。

圖2實(shí)驗(yàn)流程與scATAC-seq細(xì)胞的細(xì)胞類型

使用基因活性分?jǐn)?shù)來區(qū)分細(xì)胞類型，鑒定了多種主要的細(xì)胞類型，包括放射狀膠質(zhì)細(xì)胞(RGS)、中間前體細(xì)胞(IPC)、深層(皮質(zhì)V-VI層)興奮性神經(jīng)元(dlENs)、上層(皮質(zhì)II-IV層)興奮性神經(jīng)元(UlENs)、MGE和CGE來源的皮質(zhì)中間神經(jīng)元(分別為IN-MGES和IN-CGES)、胰島神經(jīng)元、MGE的前體細(xì)胞、小神經(jīng)膠質(zhì)細(xì)胞、少突膠質(zhì)前體細(xì)胞（OPCs）、內(nèi)皮細(xì)胞以及壁細(xì)胞。此外，結(jié)合scRNA-seq數(shù)據(jù)，使用CellWalker將scATAC-seq細(xì)胞與細(xì)胞類型匹配。CellWalker以更精細(xì)的分辨率識別細(xì)胞類型和更廣泛的細(xì)胞亞型。此外，作者可以在兩種亞型(TRGS和ORGS)之間識別出不同的開放區(qū)，這兩種亞型分別表達(dá)CRYAB和HOPX，這表明scATAC-seq能夠在高分辨率下區(qū)分細(xì)胞亞型，具有更高的靈敏度。

圖3細(xì)胞類型標(biāo)記基因的特異性

識別特定細(xì)胞類型的增強(qiáng)子

為了確定候選的基因調(diào)控元件，從每種細(xì)胞類別中找出總的單個(gè)細(xì)胞的開放區(qū)。將重疊的開放區(qū)合并，總計(jì)459,953個(gè)開放區(qū)。對基因組中開放區(qū)進(jìn)行注釋，其在基因內(nèi)含子和遠(yuǎn)端基因間隔區(qū)以及轉(zhuǎn)錄起始點(diǎn)附近區(qū)域富集。將開放區(qū)與Roadmap Epicenomics 25-state chromatin model相交，發(fā)現(xiàn)啟動(dòng)子和增強(qiáng)子state強(qiáng)烈富集，以及轉(zhuǎn)錄、異染色質(zhì)和靜止state的耗盡。每種細(xì)胞類型確定了特定于細(xì)胞類型的差異可及的開放區(qū)，共有265,123個(gè)開放區(qū)，大多數(shù)細(xì)胞類型具有數(shù)千個(gè)特異的開放區(qū)。此外，還確定了八個(gè)大腦區(qū)域之間可區(qū)別訪問的開放區(qū)信號值。

為了鑒定候選增強(qiáng)子，整合了ATAC-seq、CUT&Tag、Hic和基因表達(dá)數(shù)據(jù)，并使用接觸活性算法來預(yù)測所有類型的皮質(zhì)細(xì)胞的增強(qiáng)子-基因相互作用。共預(yù)測了25,659個(gè)基因連鎖增強(qiáng)子。其中在發(fā)育中的人皮層分離的細(xì)胞上進(jìn)行了染色質(zhì)免疫沉淀和測序(PLAC-seq)，發(fā)現(xiàn)了67,493個(gè)開放區(qū)以及10,050個(gè)預(yù)測的增強(qiáng)子，這些開放區(qū)和預(yù)測的增強(qiáng)子與預(yù)測的細(xì)胞類型特異性增強(qiáng)子相關(guān)的基因被富集。為了進(jìn)一步驗(yàn)證注釋的z確性，將開放區(qū)與從人類皮質(zhì)組織樣本生成的公開可用數(shù)據(jù)進(jìn)行關(guān)聯(lián)分析。發(fā)現(xiàn)scATAC-seq鑒定到了混合測序數(shù)據(jù)注釋的大多數(shù)開放區(qū)，并且還鑒定了許多在混合數(shù)據(jù)中未預(yù)測的細(xì)胞類型特異性peaks，特別是富集在諸如小神經(jīng)膠質(zhì)細(xì)胞和內(nèi)皮細(xì)胞的稀有細(xì)胞群體中的開放區(qū)。在功能驗(yàn)證的前腦增強(qiáng)子中，大多數(shù)(319個(gè)中的304個(gè))與染色質(zhì)可及性開放區(qū)重疊，但只有67個(gè)重疊的增強(qiáng)子是通過激活接觸方法進(jìn)行預(yù)測?？傊?，這些分析表明，scATAC-seq是一種從異質(zhì)組織樣本中檢測染色質(zhì)可及性模式的可靠方法。然而，預(yù)測的增強(qiáng)子與以前發(fā)表的研究的有限重疊表明，迫切需要更好地理解預(yù)測調(diào)控潛力的計(jì)算機(jī)算法的相關(guān)特征。

為了確定細(xì)胞類型的調(diào)節(jié)“語法”，計(jì)算了已知的轉(zhuǎn)錄因子結(jié)合基序在細(xì)胞類型特異的開放區(qū)中的富集。轉(zhuǎn)錄因子基序的富集與標(biāo)記基因富集的細(xì)胞類型注釋密切相關(guān)。為了在單細(xì)胞水平上檢查轉(zhuǎn)錄因子基序的富集，作者還使用了ChromVAR(方法)，發(fā)現(xiàn)每個(gè)cluster的排列靠前的基序富集基本一致。這些發(fā)現(xiàn)確定了scATAC-seq在不同細(xì)胞類型中識別出與已知轉(zhuǎn)錄因子表達(dá)模式一致的染色質(zhì)可及性模式，并為發(fā)現(xiàn)細(xì)胞圖譜和細(xì)胞命運(yùn)的轉(zhuǎn)錄因子編碼提供了路徑。

圖4 不同細(xì)胞類型、原始scATAC-seq細(xì)胞中每個(gè)區(qū)域的特異性peaks以及預(yù)測增強(qiáng)子-基因相互作用

調(diào)控區(qū)突變引發(fā)的疾病風(fēng)險(xiǎn)

不同發(fā)育階段和分化狀態(tài)的染色質(zhì)狀態(tài)的細(xì)胞分辨率數(shù)據(jù)可能提供突變和發(fā)育中的不同細(xì)胞類型選擇性脆弱性之間的聯(lián)系。為此，研究者將特定細(xì)胞類型的ATAC-seq peaks和預(yù)測的增強(qiáng)子與疾病相關(guān)的常見和罕見的非編碼突變體記性關(guān)聯(lián)分析。首先分析了細(xì)胞類型特異性的peaks sets、預(yù)測的增強(qiáng)子、啟動(dòng)子相互作用區(qū)域與基因組區(qū)域重疊的peaks值，這些區(qū)域在發(fā)育遲緩的個(gè)體中富集了拷貝數(shù)變異；結(jié)果發(fā)現(xiàn)dlEN、內(nèi)皮/壁層和小神經(jīng)膠質(zhì)細(xì)胞特異性peaks顯著富集，以及中間神經(jīng)元中重疊啟動(dòng)子相互作用區(qū)的peaks顯著富集。

由于這些區(qū)域不提供針對單個(gè)調(diào)控元件或基因的特異性，所以又測試了特定細(xì)胞類型的peaks、預(yù)測的增強(qiáng)子以及與ASD、神經(jīng)發(fā)育延遲(NDD)相關(guān)基因附近區(qū)域的啟動(dòng)子相互作用區(qū)域重疊peaks的富集情況，并發(fā)現(xiàn)在大多數(shù)細(xì)胞類型中這些區(qū)域顯著富集或缺失。還對ASD和NDD患者的細(xì)胞類型特異性peaks和預(yù)測的DNMs的增強(qiáng)子進(jìn)行了比較分析；但與同胞家系相比，在先證者中注釋的DNMs沒有顯著富集。此外，將預(yù)測的增強(qiáng)子與包含神經(jīng)發(fā)育疾病相關(guān)基因的拓?fù)浣Y(jié)構(gòu)域(TADS)相交，發(fā)現(xiàn)在幾種細(xì)胞類型的TADS中有顯著的共定位。

最后，研究者評估了在預(yù)測的每種細(xì)胞類型的增強(qiáng)子中與神經(jīng)精神疾病風(fēng)險(xiǎn)相關(guān)的常見變異的富集程度。對精神分裂癥、ASD、重度抑郁癥和雙相情感障礙進(jìn)行了部分遺傳連鎖不平衡(LD)回歸分析。結(jié)果發(fā)現(xiàn)，與精神分裂癥相關(guān)的常見突變的興奮性和抑制性神經(jīng)元增強(qiáng)子富集?？偠灾?，細(xì)胞類型特異性染色質(zhì)狀態(tài)數(shù)據(jù)有可能在皮質(zhì)發(fā)育過程中確定特定的調(diào)控程序，這些程序很有可能導(dǎo)致神經(jīng)發(fā)育障礙，特別是提高的疾病相關(guān)的突變。

圖5：scATAC-seq peaks中疾病相關(guān)變異的富集和耗盡

神經(jīng)發(fā)生中染色質(zhì)的動(dòng)態(tài)變化

為了更好地理解轉(zhuǎn)錄組和表觀基因組變化如何在神經(jīng)發(fā)生過程調(diào)節(jié)細(xì)胞命運(yùn)，聯(lián)合分析了來自視覺皮質(zhì)的相關(guān)細(xì)胞類型的scRNA-seq和scATAC-seq數(shù)據(jù)?；虮磉_(dá)和基基因活性在共嵌空間的預(yù)測顯示，不同細(xì)胞類型分群不受特征模式的影響。為了確定興奮性神經(jīng)元分化和成熟后的染色質(zhì)可及性軌跡，在共嵌空間對細(xì)胞進(jìn)行了擬時(shí)間排序，恢復(fù)了已知的興奮性神經(jīng)元分化的細(xì)胞類型的發(fā)育序列。發(fā)現(xiàn)了超過25000個(gè)開放區(qū)，在擬時(shí)間內(nèi)具有瞬時(shí)可及性，包括超過5000個(gè)預(yù)測的增強(qiáng)子，其中許多與細(xì)胞類型識別相關(guān)的基因相互作用。

圖6：人腦皮層神經(jīng)發(fā)生過程中染色質(zhì)可及性的動(dòng)態(tài)變化

區(qū)域特異性染色質(zhì)狀態(tài)

在神經(jīng)發(fā)生的早期，出現(xiàn)了特定區(qū)域類型的皮質(zhì)興奮性神經(jīng)元，但在不同區(qū)域的前體細(xì)胞之間只發(fā)現(xiàn)了有限的轉(zhuǎn)錄差異。調(diào)控元件的可及性變化通常先于基因表達(dá)。研究者比較了頭-尾軸兩端、PFC和V1的興奮性譜系細(xì)胞的scRNA-seq和scATAC-seq圖譜。在每種模式中，對細(xì)胞進(jìn)行擬時(shí)間分析，確定分化軌跡，并沿著該軌跡確定PFC和V1譜系之間轉(zhuǎn)錄或染色質(zhì)狀態(tài)差異變得明顯的“分支點(diǎn)”。與轉(zhuǎn)錄組數(shù)據(jù)揭示了興奮性神經(jīng)元的區(qū)域特異性分群不同的是，染色質(zhì)狀態(tài)特征表明PFC和V1中間前體細(xì)胞群體之間存在顯著差異。在轉(zhuǎn)錄水平上，PFC和V1IPCs只有少數(shù)幾個(gè)基因有差異表達(dá)，而染色質(zhì)可及性分析發(fā)現(xiàn)這些細(xì)胞類型之間有1800多個(gè)差異可及區(qū)。

為了確定可能導(dǎo)致PFC和V1譜系分化的調(diào)控程序，又對PFC和V1細(xì)胞之間可區(qū)分的開放區(qū)進(jìn)行了轉(zhuǎn)錄因子結(jié)合位點(diǎn)富集分析。該分析鑒定了幾個(gè)轉(zhuǎn)錄組研究預(yù)測的轉(zhuǎn)錄因子。該分析還在PFC細(xì)胞中發(fā)現(xiàn)了視黃酸(RA)信號通路的組成部分——在妊娠中期的PFC中RA活性增加。

圖7：祖細(xì)胞染色質(zhì)狀態(tài)的區(qū)域差異預(yù)示著區(qū)域特異性興奮神經(jīng)元的出現(xiàn)

視黃酸在大腦皮層區(qū)域化中的作用

在哺乳動(dòng)物大腦發(fā)育的過程中，視黃酸信號在神經(jīng)組織中發(fā)揮著重要的作用。為了測試RA是否促進(jìn)人類PFC的分化，在含有或不含有維生素A(RA合成的前體)的情況下培養(yǎng)皮質(zhì)類器官。同時(shí)利用DEAB(一種有效的RA合成抑制劑)處理與維生素A一起培養(yǎng)的腦類器官。在分化的第10周，此時(shí)深層神經(jīng)正在發(fā)生，利用scRNA-seq分析了腦類器官。發(fā)現(xiàn)在維生素A存在下培養(yǎng)的興奮性前腦神經(jīng)元與不含維生素A或在DEAB存在下培養(yǎng)的腦類器官的神經(jīng)元是分開聚類的。還發(fā)現(xiàn)了區(qū)分PFC和V1皮層神經(jīng)元的信號。利用目前的分類方法注釋類器官神經(jīng)元中的PFC和V1神經(jīng)元的特征；與未使用維生素A培養(yǎng)或使用DEAB處理的神經(jīng)元相比，維生素A培養(yǎng)的類器官中PFC類似神經(jīng)元的比例始終較高。免疫染色證實(shí)了基因產(chǎn)物的差異表達(dá)，并發(fā)現(xiàn)與含有維生素A培養(yǎng)的腦類器官的PFC特性一致。這些發(fā)現(xiàn)表明，RA信號通路有助于人類皮質(zhì)發(fā)育過程中PFC神經(jīng)元譜系的發(fā)育。

圖8：視黃酸在發(fā)育中的皮質(zhì)PFC-V1分裂的作用

基準(zhǔn)測試腦類腦類器官

由于人體原代組織比較稀缺，對人類神經(jīng)發(fā)育的研究需要合適的體外模型，如腦類器官。之前單細(xì)胞轉(zhuǎn)錄組學(xué)和整體表觀基因組學(xué)強(qiáng)調(diào)了大腦類器官細(xì)胞與體內(nèi)同類細(xì)胞之間的相似性。利用三個(gè)遺傳正常的樣本在三個(gè)分化時(shí)間點(diǎn)的23,555個(gè)皮質(zhì)類器官細(xì)胞進(jìn)行了scATAC-seq。為了驗(yàn)證所產(chǎn)生的細(xì)胞系，從來自同一細(xì)胞系的類器官中生成了scRNA-seq數(shù)據(jù)，并進(jìn)行了平行培養(yǎng)，結(jié)果表明所有細(xì)胞系都表達(dá)FOXG1和主要細(xì)胞類型的標(biāo)記。scATAC-seq數(shù)據(jù)中確定了主要的細(xì)胞類型，但是單個(gè)分群包含的細(xì)胞類型特異性開放區(qū)比來自原代細(xì)胞少。隨后，通過從原代細(xì)胞的開放區(qū)量化了類器官細(xì)胞之間的染色質(zhì)可及性。在原代細(xì)胞中共確定了459,953個(gè)開放區(qū)，然而在類器官只確定了239,661個(gè)開放區(qū)。隨后又發(fā)現(xiàn)在類器官中未檢測到包括小膠質(zhì)細(xì)胞、內(nèi)皮細(xì)胞、星形膠質(zhì)細(xì)胞和OPCs等細(xì)胞類型的特定開放區(qū)。當(dāng)從分析中去除這些開放區(qū)的peaks信號后，在類器官開放區(qū)中尚未發(fā)現(xiàn)V1興奮神經(jīng)元富集。然而，除了小膠質(zhì)細(xì)胞增強(qiáng)子外，大多數(shù)(超過80%)的預(yù)測增強(qiáng)子都存在于類器官中。還鑒定了在原始細(xì)胞中沒有發(fā)現(xiàn)的109,960個(gè)類器官開放區(qū)。

圖9：腦類器官和原代組織的peaks比較揭示染色質(zhì)圖譜的顯著差異

討論

在這項(xiàng)研究中，分析了發(fā)育中的人腦中單個(gè)細(xì)胞的染色質(zhì)狀態(tài)，并發(fā)現(xiàn)了數(shù)千個(gè)追蹤神經(jīng)元分化的瞬時(shí)可及位點(diǎn)。這些狀態(tài)可能揭示了在神經(jīng)發(fā)生過程中控制細(xì)胞命運(yùn)建立的機(jī)制，將它們與來自成人大腦的數(shù)據(jù)關(guān)聯(lián)分析，可能使表觀基因組神經(jīng)發(fā)育軌跡的完全重建成為可能。與以前的研究一致，與疾病變異相關(guān)的染色質(zhì)與神經(jīng)精神障礙發(fā)病中的有絲分裂期后發(fā)育的皮質(zhì)興奮性神經(jīng)元相關(guān)。未來需要探索這些調(diào)控區(qū)域中的疾病相關(guān)變異如何改變發(fā)育中的皮層中的細(xì)胞命運(yùn)。通過比較發(fā)育過程中不同皮層區(qū)域(V1和PFC)的調(diào)控情況，發(fā)現(xiàn)了不同的轉(zhuǎn)錄因子結(jié)合位點(diǎn)，它們在這兩個(gè)譜系之間存在差異富集。該結(jié)果擴(kuò)展了RA信號在前腦發(fā)育中的作用，表明RA信號有助于人類PFC興奮性神經(jīng)元的發(fā)育。RA信號調(diào)節(jié)失調(diào)與一系列神經(jīng)發(fā)育和精神疾病有關(guān)，因此該研究結(jié)果可能促進(jìn)對這些疾病的研究。

如果對單細(xì)胞空間轉(zhuǎn)錄組系列感興趣的老師，歡迎來電咨詢，百邁客為您提供更詳盡的技術(shù)服咨詢和測序服務(wù)。