轉(zhuǎn)錄控制和染色質(zhì)活性總體上涉及到調(diào)節(jié)性RNA招募特定的RNA結(jié)合蛋白（RBP）的作用。盡管其中涉及多個RBP，但仍不清楚RBP如何直接作用于染色質(zhì)。作者進行大規(guī)模的RBP ChIP-seq分析，揭示了人類基因組中活躍的染色質(zhì)區(qū)域中廣泛存在的RBP。像轉(zhuǎn)錄因子（TFs）一樣，RBP對基因組中的熱點，尤其是基因啟動子，也表現(xiàn)出強烈的偏好，因為它們的關(guān)聯(lián)經(jīng)常與轉(zhuǎn)錄輸出相關(guān)。聚類揭示了TF和RBP之間廣泛的關(guān)聯(lián)，例如YY1（已知的RNA依賴性TF）和RBM25（參與剪接調(diào)控的RBP）。最后通過HI-C和CHIA-PET分析RBM25的消耗對所有YY1依賴的活性，包括染色質(zhì)結(jié)合，DNA環(huán)化和轉(zhuǎn)錄。

圖1.核RNA結(jié)合蛋白普遍存在于基因啟動子中，許多通過與特定轉(zhuǎn)錄因子的功能相互作用直接參與轉(zhuǎn)錄。

    為了廣泛研究RBP在染色質(zhì)水平上的潛在功能，作者根據(jù)其細胞核定位、先前的研究結(jié)果、結(jié)合結(jié)構(gòu)域和功能類別等等條件，分別對HepG2和K562細胞中的58和45個RBP進行了系統(tǒng)的ChIP-seq分析，根據(jù)分析發(fā)現(xiàn)大部分的RBP表現(xiàn)出與染色質(zhì)的廣泛而特異性的相互作用（圖1A）。在大多數(shù)情況下，RBP在HepG2和K562細胞中均表現(xiàn)出強結(jié)合力，結(jié)合ENCODE注釋的染色質(zhì)狀態(tài)和DNase I超敏位點發(fā)現(xiàn)RBP偏愛開放的染色質(zhì)區(qū)域，并且通常為CTCF結(jié)合位點和CpG島，所有與染色質(zhì)相關(guān)的RBP普遍偏愛基因啟動子，而不同的RBP表現(xiàn)出對不同啟動子的偏好，一些特定的啟動子被多個RBP結(jié)合（圖1B）。同時RBP-染色質(zhì)互作傾向于與活性組蛋白修飾（例如，H3K27ac，H3K4me3和H3K36me3）呈正相關(guān)，而與抑制性H3K9me3標記負相關(guān)（圖1C）?？傮w而言，在顯示出可檢測的染色質(zhì)結(jié)合的RBP中30％–40％的具有生化活性的染色質(zhì)區(qū)域至少與一個RBP結(jié)合（圖1D）。

    為了全面表征單個RBP的結(jié)合偏好，作者將每個RBP的ChIP-seq峰分配給了7個帶有ENCODE注釋的基因組片段，并評估了這些片段中每個RBP的峰的相對分布（圖1E）。這些數(shù)據(jù)清楚地表明，活性啟動子是RBP的熱點。各個RBP對每種基因組分割類型的相對貢獻如圖1E相對氣泡大小所示，相對于其他RBP，HNRNPK是阻遏區(qū)的主要RBP。相反，AGO2在HepG2細胞的轉(zhuǎn)錄區(qū)域富集。有趣的是，三個RBP XRCC5，HNRNPL和RBM25似乎比其他RBP更普遍地與啟動子和增強子相連。

圖2.染色質(zhì)相關(guān)RBP的一般特征

（A）通過ChIP-seq在HepG2和K562細胞中調(diào)查的RBP概述。深藍色，符合ENCODE標準的高質(zhì)量數(shù)據(jù)；淺藍色，ChIP-seq數(shù)據(jù)，符合所有其他ENCODE標準；灰色，IP后無信號富集，盡管通過蛋白質(zhì)印跡檢測到有效IP；白色，未經(jīng)調(diào)查。

（B）一個典型的基因組區(qū)域，紅色和橙色突出顯示的ChromHMM片段分別對應于啟動子和增強子。

（C）Circos圖顯示了集體RBP-染色質(zhì)相互作用，DNase I超敏檢測到的開放染色質(zhì)區(qū)域和HepG2細胞中關(guān)鍵組蛋白修飾事件之間的關(guān)系。

（D）HepG2和K562細胞中與染色質(zhì)相關(guān)的RBP覆蓋單個組蛋白修飾事件，以及與至少一種生化活性相關(guān)的所有染色質(zhì)區(qū)域的累積覆蓋率（紅線）。

（E）RBP在HepG2細胞中特定狀態(tài)的ENCODE注釋的基因組分段中的占有率。R，抑制區(qū)域； PF，啟動子側(cè)翼區(qū)； T，轉(zhuǎn)錄區(qū)； CTCF，CTCF結(jié)合位點； WE，弱效增強劑； E，增強劑； TSS，轉(zhuǎn)錄起始位點/啟動子。

    啟動子是RBP-染色質(zhì)相互作用的主要位置，作者進一步研究單個RBP的啟動子結(jié)合情況。發(fā)現(xiàn)RBP共同優(yōu)先結(jié)合用H3K4me3和H3K27me3標記的二價啟動子以及單獨H3K4me3修飾的活性啟動子和CpG，不同的是在HepG2細胞中，二價基因比僅H3K4me3基因更為普遍，在K562細胞中則相反（圖2A），這表明染色質(zhì)相關(guān)的RBP可能積極參與細胞類型特異性基因表達程序。關(guān)于在不同類別基因的啟動子，相對于蛋白質(zhì)編碼和lncRNA基因的啟動子，四個RBP（即RBM22，PRPF4，HNRNPUL1和SNRNP70）在小RNA基因的啟動子上表現(xiàn)出額外的富集（圖2B）。此外，PRPF4在tRNA基因啟動子上高度富集，SNRNP70在小核仁RNA（snoRNA）基因啟動子上高富集（圖2C）。

   圖4D表示了RBP ChIP-seq峰相對于轉(zhuǎn)錄起始位點（TSS）的位置，盡管所有RBP在TSS的兩側(cè)均顯示結(jié)合，但出現(xiàn)三類明顯的RBP-染色質(zhì)相互作用：（1）上游TSS（即RBM25），（2）以TSS為中心（即，GTF2F1）和（3）下游TSS（即RBFOX2），第三類代表兩種細胞類型中的大多數(shù)RBP（圖2D）。 GTF2F1以TSS為中心的結(jié)合模式很可能反映出其作為核心轉(zhuǎn)錄機制的一部分的功能，暗示具有類似關(guān)聯(lián)模式的其他RBP可能具有相關(guān)功能。大多數(shù)RBP與TSS下游序列的關(guān)聯(lián)表明，許多RBP可能以新生的RNA依賴性方式與染色質(zhì)相互作用。

圖3.不同啟動子類別上的RBP-染色質(zhì)互作模式

（A）通過序列背景（帶有或不帶有CpG島）或組蛋白修飾特征（例如以H3K4me3和H3K27me3信號標記的二價啟動子）分離的啟動子亞組的總體RBP偏好性，僅包含H3K4me3或H3K27me3信號的啟動子。正確：相對豐富的關(guān)鍵。

（B）在不同類別的基因啟動子上各個RBP的相對占用頻率。

（C）相對于背景分布，六類小RNA基因啟動子之間的RBP ChIP-seq峰分布。

（D）TSS周圍的復合RBP綁定信號。 RBP根據(jù)信號最大值與TSS的相對位置進行排序。

由于RBP與基因啟動子的聯(lián)系緊密，作者接下來通過檢測每個RBP-啟動子互作與基因轉(zhuǎn)錄之間的關(guān)系，通過RNA-seq比較了每個RBP敲除之前和之后基因表達的變化。在敲除之前，大多數(shù)RBP-啟動子互作與靶基因的轉(zhuǎn)錄活性相關(guān)（圖3A，左），在敲除后，幾乎所有RBP都影響基因表達（圖3A，中間），這些RBP中至少有六個對轉(zhuǎn)錄有直接影響（圖3A，右）。為了確定誘導的基因表達是否與其啟動子結(jié)合相關(guān)，作者進一步計算了每個RBP的優(yōu)勢比，這六個RBP均顯示出明顯的優(yōu)勢比（> 1，p <0.05），表明這些RBP直接參與轉(zhuǎn)錄控制（圖3B）。另外，作者在GRO-seq和RNA-seq信號之間進行了比較，觀察到在穩(wěn)態(tài)下轉(zhuǎn)錄相對于RNA水平的總體變化的三種不同模式：對于HNRNPLL，在敲除后觀察到RNA-seq和GRO-seq之間幾乎沒有一致性（圖3C，頂部），表明該RBP可能獨立調(diào)節(jié)轉(zhuǎn)錄和轉(zhuǎn)錄后事件。對于RBM25有正相關(guān)，表明RBM25介導的轉(zhuǎn)錄在穩(wěn)態(tài)下可能有助于RNA水平（圖3C，中間）。相反，檢測到與XRCC5呈負相關(guān)（圖3C，底部）。

由于啟動子結(jié)合可能會指導下游RNA加工事件，例如RNA穩(wěn)定性，輸出或翻譯。為此，作者進行HepG2細胞中的細胞分級測序（CeFra-seq）和RBP敲低誘導的剪接變化（通過“剪接百分比”或“ PSI），確定與RBP染色質(zhì)結(jié)合活性相關(guān)的功能性。觀察到mRNA的核保留與在基因啟動子處檢測到的RBP數(shù)量成反比（圖3D），編碼和非編碼基因均是如此。另外還發(fā)現(xiàn)，不同的RBP對這種效果的貢獻不同（圖3E）。特別是，RBM25的啟動子締合活性最能預測核保留（圖3F）。

圖4. RBP-啟動子相互作用與基因表達之間的相關(guān)性

（A）RBP與啟動子結(jié)合的概率與GRO-seq鑒定的靶基因轉(zhuǎn)錄活性之間的相關(guān)性（左）； RNA-seq（中），GRO-seq（右）對單個RBP敲低的反應。

（B）由GRO-seq確定的RBP占用的啟動子與非占用啟動子相比的轉(zhuǎn)錄反應的幾率。

（C）敲低三個代表性RBP時，RNA-seq和GRO-seq分析的基因表達變化之間的比較。

（D）其啟動子被不同的RBP占據(jù)的每組基因的RNA核保留指數(shù)（核/（核+胞質(zhì)））的分布。

（E）由隨機森林算法確定的變量重要性，以評估每個變量的預測能力，顯示了前十個RBP。

（F）具有或不具有RBM25結(jié)合在其啟動子上的證據(jù)的基因的RNA核保留分布

    為了了解RPB可能如何影響轉(zhuǎn)錄，RBP-染色質(zhì)互作如何與特定的TF相互協(xié)調(diào)。作者將RBP ChIP-seq數(shù)據(jù)與TF的ChIP-seq數(shù)據(jù)進行了整合，并通過使用新開發(fā)的非負矩陣分解（NMF）方法分析了它們的共共結(jié)合。通過分析找到了17個最佳組號（圖4A）。根據(jù)ENCODE注釋，這些因子組中的每一個都顯示出對不同基因組區(qū)域的特定偏好（圖4B）。已知CTCF，RAD21和SMC3這3種蛋白在介導長距離基因組相互作用中相互協(xié)作，分析發(fā)現(xiàn)這些蛋白質(zhì)都主要位于HepG2細胞中富含CTCF的基因組片段中（圖4B），這與已知的一致。當前分析的RBP可以其中一組（第10組）僅由RBP組成（圖4C），每組中的多個因素都在蛋白質(zhì)水平上發(fā)生物理相互作用（圖4C中的藍色邊緣）。通過根據(jù)人類基因組HOT區(qū)域中峰的比例對RBP進行排名并將其分為四個升四分位數(shù)，發(fā)現(xiàn)第4組中超過一半的RBP與各種TF一起排在了前四分位數(shù)中，包括PRPF4， SRSF4等（圖4D）。因此，與TF一樣，RBP也顯示出對人類基因組中HOT區(qū)的高度偏好。

   在互作網(wǎng)絡分析中，TF YY1分為兩個獨立的聯(lián)合組：第9組，包含多個TF和RBP NONO，是參與RNA加工，DNA修復和重組的多功能核蛋白；第13組，由YY1，XRCC5和RBM25（圖4C）組成，已知前兩者彼此互作。YY1基因組結(jié)合位點在HepG2細胞中也表現(xiàn)出RBM25結(jié)合（圖4E，上圖），鑒于最近發(fā)現(xiàn)YY1似乎以依賴RNA的方式與DNA相互作用，RBM25可能以細胞類型特異性方式在調(diào)節(jié)YY1的轉(zhuǎn)錄功能中發(fā)揮作用。為了探索這種可能性作者分別分析了RBM25和YY1結(jié)合圖譜與TSS周圍的YY1結(jié)合基序分布的關(guān)系（圖4E，底部），觀察到核心YY1結(jié)合基序在YY1和RBM25共結(jié)合的啟動子區(qū)域富集或單獨的YY1結(jié)合，但不是單獨的RBM25結(jié)合，表明核心YY1結(jié)合基序介導YY1結(jié)合。

圖5. HepG2細胞中與染色質(zhì)相關(guān)的RBP和TF的綜合分析

（A）通過NMF推斷的系數(shù)矩陣將染色質(zhì)相關(guān)的RBP和TF分為17組。

（B）由單個NMF組（左）覆蓋和注釋總的順式調(diào)控元件（CRE，在圖1E中定義）。未注釋的區(qū)域。右側(cè)顯示了每個組所占的不同CRE的分數(shù)。

（C）代表NMF隔離的團體。藍線，由GeneMANIA注釋的每個組中成員之間的已知物理相作用。

（D）RBP與HOT區(qū)域的優(yōu)先關(guān)聯(lián)?；赗BP在HOT啟動子區(qū)域中的相對結(jié)合，它們被分為四個四分位數(shù)（灰色系）。 y軸：落入HOT區(qū)域的總峰的百分比。

（E）頂部：Rep25，XRCC5和YY1在HepG2細胞中的共定位。韋恩圖的每個部分被進一步量化為每個RBP的峰的百分比，基于此可計算出單獨的成對共定位。下圖：相對于YY1，RBM25和YY1結(jié)合峰的核心YY1結(jié)合基序的分布。

    為了追蹤RBM25在介導YY1依賴性轉(zhuǎn)錄中的潛在功能，作者首先驗證了RNA依賴性YY1在染色質(zhì)上的結(jié)合。我們從ENCODE數(shù)據(jù)中選擇了多個在HepG2細胞中顯示YY1結(jié)合峰的基因啟動子，并用5,6dichloro-1-β-D-核呋喃糖基-苯并咪唑阻止轉(zhuǎn)錄然后對YY1和RBM25進行了ChIP-qPCR，在DRB處理后，所有靶基因座上的YY1和RBM25結(jié)合都減弱了，在DRB洗脫后恢復了（圖5A，5B）。接下來通過coIP 實驗驗證YY1和RBM25形成復合物（圖5C）。siRNA介導的YY1或RBM25耗竭對每種蛋白質(zhì)的特異性敲除作用，但對非靶向蛋白質(zhì)卻沒有作用（圖5D）。細胞的GRO-seq分析結(jié)果顯示出高度可重復的模式：在每種情況下均鑒定出大量上調(diào)或下調(diào)的基因，如圖5E所示。GRO-seq信號的倍數(shù)變化（FC）之間的比較顯示YY1和RBM25敲低誘導的基因表達之間存在一致性。由于每種因子的敲低都會導致大量基因的下調(diào)和上調(diào)，因此在通常受影響的基因中，絕大多數(shù)上調(diào)和下調(diào)的基因在相同方向上受到影響（圖5G）。此外，RBM25或YY1敲低的基因上調(diào)和下調(diào)都與它們的啟動子關(guān)聯(lián)，如在所有成對比較中均具有很高的比值比（圖5H）。

圖6（A+B） YY1和RBM25在HepG2細胞（A和B）中對基因表達的共調(diào)控ChIP-qPCR分析了DRB處理后和DRB洗脫后結(jié)合在代表性靶基因啟動子上的YY1（A）和RBM25（B）。

（C）RBM25和YY1的相互coIP。 

（D）RBM25或YY1的有效敲低而不影響其他蛋白質(zhì)。

（E）在兩個代表性基因位點上的RBM25和YY1 ChIP-seq譜圖，一個通過GRO-seq確定，一個用于上調(diào)（左）基因，一個用于下調(diào)（右）基因。

（F）對RBM25和YY1敲低的轉(zhuǎn)錄反應的整體比較。標明了倍數(shù)變化（FC）的斯皮爾曼相關(guān)系數(shù)（SCC）。

（G）在耗盡了RBM25或YY1的HepG2細胞中差異表達的基因的重疊。彩色框和線型分別表示由RBM25（實線）或YY1（虛線）調(diào)節(jié)的向上（紅色框）或向下（藍色框）基因表達事件。

（H）靶基因啟動子處的RBM25和YY1結(jié)合與通過比值比確定的誘導基因表達正相關(guān)。

    為了了解YY1和RBM25協(xié)調(diào)結(jié)合和調(diào)控轉(zhuǎn)錄的機制，首先在以YY1和RBM25為靶標的啟動子的同一面板上進行ChIP-qPCR，發(fā)現(xiàn)RBM25敲低顯著降低了所有六個靶基因上的YY1結(jié)合（圖6A），但反之則沒有（圖6B）。這些結(jié)果表明RBM25首先作用于那些靶啟動子，這是隨后的YY1結(jié)合所必需的。為了獲得對此的整體證據(jù)，作者分別在耗盡RBM25之前和之后對YY1進行了ChIP-seq，在耗盡YY1之前和之后對RBM25進行了ChIP-seq，大多數(shù)（即使不是全部）YY1結(jié)合事件在RBM25耗盡后減少，但是RBM25 ChIP-seq信號在YY1耗盡后基本上保持不變（圖6C）。對啟動子，增強子和CTCF結(jié)合位點上的RBM25 ChIP-seq信號的元基因分析表明，RBM25 ChIP-seq信號對YY1耗盡基本不敏感（圖6D）。YY1 ChIP-seq信號在RBM25相關(guān)的基因組位點更為普遍，而RBM25的耗盡顯著減弱了所有三類DNA元件上的YY1結(jié)合（圖6E）。統(tǒng)計分析表明，與沒有RBM25關(guān)聯(lián)的證據(jù)相比，YY1結(jié)合中的FC對RBM25占據(jù)的基因組位點的影響更大（圖6F）。這些數(shù)據(jù)表明，RBM25是在HepG2細胞中有效靶向全基因組的YY1所必需的。

    那么RBM25如何調(diào)節(jié)此類YY1介導的基因組互作。作者在敲除RBM25之前和之后在HepG2細胞上進行了BL-Hi-C測序和雙末端標簽測序（ChIA- PET）。發(fā)現(xiàn)，RBM25敲減誘導的啟動子錨定PET的整體變化略有減少（圖6G），這可能反映了RBM25耗竭對轉(zhuǎn)錄的直接和間接影響。使用類似于然后在處理來自YY1缺失的mESC的HiChIP數(shù)據(jù)時的嚴格性，并著眼于顯著差異的啟動子-增強子相互作用的變化。發(fā)現(xiàn)，與啟動子/增強子（PE）連接的PET在人類基因組中所有帶注釋的啟動子之間在兩個方向上都顯示出相似的變化（圖6H）。與共結(jié)合的啟動子和增強子相比，單獨與YY1結(jié)合的啟動子和增強子顯示出更弱的變化（圖6H）?？傮w而言，這些數(shù)據(jù)通過調(diào)節(jié)YY1募集和YY1介導的基因組相互作用證明了RBM25在YY1依賴性轉(zhuǎn)錄中的作用。

圖7.人類基因組中的RBM25依賴性YY1結(jié)合

（A和B）敲除其他因子后，YY1（A）和RBM25（B）與代表性靶基因啟動子結(jié)合的ChIP-qPCR分析。 

（C）具有代表性的基因組位點，顯示下調(diào)的YY1結(jié)合（ChIP-seq）和YY1介導的染色質(zhì)環(huán)化（BL-Hi-C），以及下（上）-/上（下）調(diào)控的基因表達數(shù)字：敲除RBM25時通過GRO-seq在個別條件下測得的表達水平。陰影區(qū)域突出顯示了靶基因啟動子（綠色）與其增強子（藍色）之間相互作用的降低。 

（D和E）RBM25（D）或YY1（E）啟動子（左），增強子（中間）和CTCF結(jié)合位點（右）上ChIP-seq信號的有基因分析，有或沒有YY1（D）或RBM25（E）綁定以及YY1（D）或RBM25（E）敲除之前（藍色）或之后（紅色）。

（F）RBM25耗竭對所有基因組位點上有或沒有RBM25占用證據(jù)的YY1結(jié)合的影響的統(tǒng)計分析。

 （G）使用所有帶注釋的啟動子作為錨點，將RBM25敲低后歸一化PETs頻率的log2倍變化相對于對照樣品中歸一化PETs頻率作圖。突出顯示了與YY1結(jié)合的啟動子，其相互作用顯著增加（紅色）或減少（藍色）。

（H）基因啟動子分為四個不同的類別，并分別分析了在RBM25敲除之前和之后，它們與活性增強子的相互作用。歸一化PET頻率的log2倍變化的分布顯示為箱線圖（請參見“方法詳細信息”）。 P，啟動子； E，增強劑。核內(nèi)RBP控制染色體活性和基因的轉(zhuǎn)錄