視網(wǎng)膜基因研究：Hi-C技術(shù)繪制高分辨率人類(lèi)視網(wǎng)膜3D基因組拓撲結構圖譜

在過(guò)去視網(wǎng)膜疾病的全基因組關(guān)聯(lián)研究（GWAS）發(fā)現常見(jiàn)的變異主要在基因組非編碼區域，但由于先導變異區域的局部連鎖不平衡（LD），阻礙了特定關(guān)聯(lián)基因和變異的識別，因此闡明與性狀相關(guān)的遺傳變異如何影響調節基因表達網(wǎng)絡(luò )，并因此影響表型，亟需了解相關(guān)組織和細胞類(lèi)型的3D基因組拓撲結構。

近日，來(lái)自美國NIH國家眼科研究所Anand Swaroop團隊通過(guò)Hi-C技術(shù)繪制了人類(lèi)視網(wǎng)膜細胞染色質(zhì)三維結構圖譜，揭示了人類(lèi)視網(wǎng)膜中超級增強子（SEs）的3D染色質(zhì)組織的特性，并將得到的視網(wǎng)膜基因組調控網(wǎng)絡(luò )與通過(guò)GWAS識別的老年性黃斑變性（AMD）和青光眼相關(guān)的數量性狀基因座（quantitativetraitloci, QTL）和遺傳變異相結合，揭示了順式調控元件(CREs)的靶標，表明了人類(lèi)視網(wǎng)膜中SEs的3D染色質(zhì)組織的特性。此次人類(lèi)視網(wǎng)膜高分辨率基因組結構的研究為組織特異性功能的遺傳調控提供了新見(jiàn)解，并能夠解剖常見(jiàn)的致盲疾病表型。

研究背景

在不同組織和細胞類(lèi)型中人類(lèi)基因組染色質(zhì)結構增強子-啟動(dòng)子相互作用建立了獨特的時(shí)空基因表達模式。非編碼遺傳變異可以通過(guò)染色質(zhì)空間構象從而調控基因表達來(lái)影響細胞表型。為了闡明人類(lèi)視網(wǎng)膜的基因組調控，利用Hi-C技術(shù)以高分辨率揭示染色質(zhì)并與超級增強子（SE）、組蛋白標記、CTCF結合和選擇轉錄因子相互作用，表明與邊緣SE的TAD相比，具有中心SE的拓撲相關(guān)結構域（TAD）表現出更強的絕緣性，并與視網(wǎng)膜基因的接觸增強。全基因組表達數量性狀位點(diǎn)（eQTL）分析與拓撲圖相結合，揭示了100個(gè)eQTL和與視網(wǎng)膜神經(jīng)變性相關(guān)的相應eGene之間的聯(lián)系，并發(fā)現了與年齡相關(guān)的黃斑變性和青光眼相關(guān)的易感變異的候選基因。對人類(lèi)視網(wǎng)膜高分辨率基因組結構的研究為組織特異性功能的遺傳控制提供了見(jiàn)解，提出了缺失遺傳性的模式，Hi-C組學(xué)技術(shù)為識別視網(wǎng)膜相關(guān)疾病表型提供有利工具。

研究設計

材料：一共取材5名患者捐贈視網(wǎng)膜組織（2名女性，3名男性，年齡65-77歲）;Hi-C互作取3名男性和1名女性視網(wǎng)膜組織以及HCT116細胞系；CUT&RUN取1名75歲的女性外周視網(wǎng)膜組織；ATAC-seq取2名女性和3名男性視網(wǎng)膜組織。

方法：Hi-C互作(5kb)、ATAC-seq、Cut&Run、ChIP-seq（數據集：GSE13731134）、RNA-seq（數據集）

文章亮點(diǎn)

1）本文整合表觀(guān)調控全面的技術(shù)，利用Hi-C互作、ATAC-seq、Cut&Run、ChIP-seq等多組學(xué)不僅重構了視網(wǎng)膜基因組拓撲結構，還深入剖析了三維結構對基因調控的網(wǎng)絡(luò )；

2）視網(wǎng)膜調控基因組拓撲結構與AMD和青光眼相關(guān)遺傳變異的首次詳細整合；

3）高分辨率的基因組結構圖像(5kb)，深入研究常染色質(zhì)A/異染色質(zhì)B區域、TAD邊界、Loop形成的變化；

4）利用GWAS、eQTL等疾病性狀關(guān)聯(lián)分析建立與視網(wǎng)膜神經(jīng)變性相關(guān)的相應eGene之間的聯(lián)系。

主要研究?jì)热?/h3>

1、Hi-C測序以5kb高分辨率識別人類(lèi)視網(wǎng)膜中的染色質(zhì)空間構象

利用Hi-C技術(shù)檢測了4名眼科相關(guān)疾病患者的視網(wǎng)膜組織，以超高的分辨率（5kb）的重新繪制了人類(lèi)視網(wǎng)膜細胞染色質(zhì)的組織（每個(gè)視網(wǎng)膜樣本的測序深度約為3億對reads，HCT116的深度約為5000萬(wàn)對read）。分析表明Hi-C數據在樣本間表現出高度的相似性，所有常染色體的分層相關(guān)系數(SCC)>0.97。整合所有樣本的數據一共獲得了11.48億對reads,其中有效互作reads的比例高達95%，不同染色體之間的互作trans占23.6%。

Hi-C測序數據統計（補充表1）

此次研究總共獲得了7.04億個(gè)有效染色質(zhì)互作，這些染色質(zhì)互作平均分布在染色體之間和染色體之內。

Hi-C檢測每個(gè)常染色體上染色體內和染色體間的互作

本文Hi-C以3 kb的高分辨率識別出67,841個(gè)顯著(zhù)的染色質(zhì)互作和2948個(gè)拓撲關(guān)聯(lián)域（TADs）、以及60000多個(gè)loops的形成。如預測一樣，視網(wǎng)膜中活躍表達的基因在A(yíng) Component染色質(zhì)區域，而大多數低活性的基因出現在B Component中，檢測全基因組loop環(huán)結構特征發(fā)現，在A(yíng)-A互作、A-B互作、B-B互作形成loops三種模式中，A-A互作形成loops距離小于B-B互作和A-B互作形成的loops，并且A-A形成的loops都是活躍的染色質(zhì)區域，跨度平均小于1Mb，大多數都是發(fā)生在TAD內。

高分辨率Hi-C（5kb）識別人類(lèi)視網(wǎng)膜染色質(zhì)結構

2、人視網(wǎng)膜染色質(zhì)拓撲結構具有組織特異性和物種保守性

為了評估已識別的染色質(zhì)特征在人視網(wǎng)膜組織特異性和物種保守性，這里比較了人類(lèi)神經(jīng)元（ACC）、GM12878淋巴細胞（LCL）、結腸癌細胞（HCT116，clone cancer）Hi-C測序數據，發(fā)現視網(wǎng)膜特異性基因，如OTX2、EYS、PDC，只存在于視網(wǎng)膜的A component；由于視網(wǎng)膜和ACC神經(jīng)元之間的高度相關(guān)性，比較發(fā)現PAX6基因座在神經(jīng)元和視網(wǎng)膜中的相互作用與其在兩個(gè)組織中的表達相似，而編碼OTX2和CRX兩個(gè)關(guān)鍵的視網(wǎng)膜TF，在視網(wǎng)膜中表現出大量的局部相互作用，但在神經(jīng)元中沒(méi)有，表明人視網(wǎng)膜組織拓撲關(guān)聯(lián)域具有組織特異性，并在基因調節中發(fā)揮著(zhù)重要作用。

小鼠和人類(lèi)染色質(zhì)組織之間的相似性表明了跨物種的守恒，強調了染色質(zhì)組織模式與視網(wǎng)膜基因調控的相關(guān)性。發(fā)現小鼠TAD中21.8%的小鼠基因對也位于相應的人類(lèi)TAD中，超過(guò)三分之一（35.7%）的基因對在小鼠體內通過(guò)染色質(zhì)環(huán)相互作用，在人類(lèi)視網(wǎng)膜中也是如此。

在這兩個(gè)物種中，EGF和ENPEP之間都觀(guān)察到loops的形成，這些發(fā)現表明，視網(wǎng)膜表達的基因表現出保守的染色質(zhì)拓撲結構，突顯了其對調節視網(wǎng)膜基因的重要性，并可能涉及組織特異性順式調控元件CREs。

人視網(wǎng)膜染色質(zhì)相互作用具有組織特異性和物種保守性

3、視網(wǎng)膜中SE與高度表達的組織特異性基因重疊，并富集視網(wǎng)膜TF結合motfis

染色質(zhì)互作形成的loops促進(jìn)CRE與其目標基因之間的物理相互作用，并且在控制基因表達方面發(fā)揮著(zhù)至關(guān)重要的作用。為了識別視網(wǎng)膜CRE，利用染色質(zhì)可及性檢測（ATAC-seq）確定了染色質(zhì)開(kāi)放狀態(tài)和CUT&Run檢測活性（H3K4me3，H3K27Ac，H3K4me234）和抑制性（H3K9me3）組蛋白標記互作的特征序列。通過(guò)全基因組染色質(zhì)開(kāi)放程度以及開(kāi)放區域（啟動(dòng)子、增強子、異染色質(zhì)、低豐度信號）構建了10種染色質(zhì)狀態(tài)模型，啟動(dòng)子和增強子區域富集了用于結合CTCF和關(guān)鍵視網(wǎng)膜TF，如CRX、NRL、OTX2、MEF2D、CREB和RORB，另外增強子主要在TSS上游區域富集，而啟動(dòng)子僅在TSS富集。

人視網(wǎng)膜全基因組染色質(zhì)可及性特征

利用H3K27Ac蛋白標記互作在視網(wǎng)膜順式調控元件啟動(dòng)子和活性或穩定增強子區域鑒定了1325個(gè)超級增強子（SEs），發(fā)現視網(wǎng)膜表達基因的數量與每個(gè)染色體的SE數量之間存在高度相關(guān)性。SE與轉錄組測序數據集結合分析表明，與CRE相比，SE重疊的基因的平均表達量高2.17倍，并且比100組隨機SE大小的位點(diǎn)高2.89倍，視網(wǎng)膜中關(guān)鍵的光感受器特異性基因NRL主要在活性組蛋白標記和染色質(zhì)開(kāi)放區域的大型SE中富集，與CRE或隨機區域相比，TF結合位點(diǎn)主要存在于SE中，其中編碼許多具有豐富結合基序TF的基因也與SE重疊，這表明它們在維持視網(wǎng)膜穩態(tài)方面的潛在作用。

人視網(wǎng)膜超級增強子SE的識別和特征鑒定

4、SE在TAD內部loops富集并參與基因表達調控

以往研究表明SE可能會(huì )與其他調節區域互作，形成能夠調節多個(gè)基因的大型轉錄單元。視網(wǎng)膜SE在距離較小的loops中富含，并廣泛存在與loops內的互作，幾乎所有SE參與內部互作的loops都小于1Mb，并且大多數SE主要在TAD中存在交互。

與隨機區域相比較，SE在調節相互作用中大量富集，與3,059個(gè)獨特的CRE、1個(gè)獨特的TSS和217個(gè)獨特的SE互作。在視網(wǎng)膜和/或大腦中特異性表達的基因大多與SE重疊，而與SE相互作用但不重疊的基因則是一些常規基因，也有一些SE與許多CRE相互作用，因此可能在協(xié)調基因表達模式方面發(fā)揮重要作用，這些結果均表明人視網(wǎng)膜中SE的組織特異性和調控獨特性。

人視網(wǎng)膜SE參與染色質(zhì)loops形成和互作獨特特征

5、SE的TAD內定位與靶基因的邊界絕緣和生物學(xué)功能有關(guān)

進(jìn)一步研究發(fā)現，大多數SE與桿狀細胞基因有關(guān)，反映了人類(lèi)視網(wǎng)膜樣本中桿狀細胞的優(yōu)勢。研究還證明了SE靶基因的生物學(xué)功能可能與TAD中的SE定位以及TAD邊界絕緣有關(guān)。在TADs的邊緣，SEs的存在是與較弱的絕緣有關(guān)，與其他類(lèi)型TAD相比，我們發(fā)現含SE的TAD中有更多的染色質(zhì)相互作用。值得注意的是，SE本身大多與位點(diǎn)的相互作用非常接近，表明在局部、TAD邊緣有很強的絕緣性，可能是由于CTCF結合的富集和避免附近基因的隨機激活，這與之前的研究一致。

SE靶向基因功能研究，即TSS與SE重疊或相互作用的基因，發(fā)現SE處在邊緣的TAD中都是富含應激反應基因，而視網(wǎng)膜基因主要在SE處于中心的TAD中富集，其中應激反應FOS重疊的SE位于TAD的邊緣，該TAD可能會(huì )受到超出自身TAD邊界的區域變化的影響。

SE的TAD內部定位與TAD的生物特征有關(guān)

6、染色質(zhì)loops、SEs和CREs與視網(wǎng)膜 eQTLs的聯(lián)合分析

eQTL將特定的遺傳變異與目標基因表達的變化聯(lián)系起來(lái)（下文稱(chēng)為eGene），為了識別可能與視網(wǎng)膜基因調節相關(guān)的eQTL，結合14,859個(gè)視網(wǎng)膜eQTLs，發(fā)現大多數視網(wǎng)膜eQTL都存在于活躍的染色質(zhì)中，并駐留在同一個(gè)TAD中，這為解釋基因變體對eGene的影響提供了直接機制。此外，評估了變異與其eGene的啟動(dòng)子區域（距離TSS±2.5kb）之間的物理距離，能夠區分啟動(dòng)子eQTL、遠端eQTL和與啟動(dòng)子相互作用的遠端eQTL。分別確定了2,410個(gè)CRE（CRE-eQTL）和880個(gè)SE（SE-eQTL）特殊的eQTL，其中58.5%的CRE-eQTL和69.3%的SE-eQTL也與相關(guān)的eGene相交?？傮w而言，大多數eQTL在物理上與eGene或其啟動(dòng)子重疊或者通過(guò)染色質(zhì)loops互作接近各自的eGene。

視網(wǎng)膜eQTL和與眼病相關(guān)的eQTL的表觀(guān)遺傳變異

7、基因組拓撲結構將目標基因與AMD和青光眼相關(guān)的風(fēng)險變異進(jìn)行關(guān)聯(lián)

關(guān)于青光眼和AMD的GWAS研究已經(jīng)確定了大量的非編碼變異，包括eQTLs分析在內的其他研究提供了進(jìn)一步見(jiàn)解，但許多相關(guān)位點(diǎn)的目標基因和變異仍然難以捉摸。研究人員將染色質(zhì)互作與從GWAS研究中識別的遺傳變異數據進(jìn)行了整合，并對帶有GWAS先導和LD變體的成人視網(wǎng)膜基因組拓撲結構進(jìn)行了綜合分析，發(fā)現了幾個(gè)先前未知的可能與青光眼和AMD有關(guān)的基因。例如，研究展示了一個(gè)顯著(zhù)的遠程相互作用導致目標基因變化的青光眼相關(guān)基因TFAP2B/PKHD1。此外，小鼠中TFAP2B表達的中斷導致與青光眼一致的強烈病理表型，這些發(fā)現指向了與這些疾病相關(guān)的特定候選致病基因。研究團隊發(fā)現了AMD的候選致病基因，在CTRB2/CTRB1位點(diǎn)建立了CFDP1作為AMD相關(guān)的候選基因。

AMD和青光眼變異基因通過(guò)視網(wǎng)膜染色質(zhì)loops與靶基因作用

總結

研究團隊通過(guò)將高分辨率Hi-C數據與表觀(guān)遺傳圖譜和順式調控元件整合，為人類(lèi)視網(wǎng)膜研究提供了一個(gè)重要的數據資源，促進(jìn)了基因組調控的研究。研究團隊識別了視網(wǎng)膜病變中缺失的遺傳性，以及包括AMD和青光眼在內的常見(jiàn)致盲疾病的候選因果基因和變異，該研究為將調控變異與視網(wǎng)膜疾病表型聯(lián)系起來(lái)提供了一個(gè)框架。整合的基因組調控圖譜還將有助于評估與其他常見(jiàn)視網(wǎng)膜相關(guān)疾病(如糖尿病視網(wǎng)膜病變)相關(guān)的基因，確定缺失的遺傳性，并了解遺傳性視網(wǎng)膜和黃斑疾病的基因型-表型相關(guān)性。

這項研究的首席研究員Anand Swaroop教授說(shuō):“這是視網(wǎng)膜調控基因組拓撲與老年性黃斑變性(AMD)和青光眼相關(guān)的遺傳變異的首次詳細整合，這兩種基因變異是視力喪失和失明的兩個(gè)主要原因?！盨waroop說(shuō):“擁有如此高分辨率的基因組結構圖像，將繼續為組織特異性功能的基因控制提供見(jiàn)解?！?/p>

如果您對文中測序技術(shù)感興趣，歡迎點(diǎn)擊下方按鈕聯(lián)系我們，我們將免費為您設計文章思路方案