圖X 調(diào)控通路模型。在血清細(xì)胞中，UHRF1通過(guò)H3K9me2染色質(zhì)標(biāo)記與半甲基化DNA結(jié)合而募集Dnmt1到復(fù)制灶維持DNA甲基化。血清到2i條件轉(zhuǎn)換導(dǎo)致5mC的被動(dòng)損失蛋白質(zhì)水平的UHRF1下調(diào)和H3K9me2的喪失導(dǎo)致復(fù)制灶的募集受損。

    DNA去甲基化動(dòng)力學(xué)可歸因于三種主要途徑:( 1）維持DNA甲基化或復(fù)制依賴性被動(dòng)稀釋，（2）從頭DNA甲基化，和（3）活性DNA去甲基化，主要通過(guò)DNA羥甲基化（圖1A）。在從血清到2i ESC的轉(zhuǎn)化期間，維持甲基化組分Dnmt1和Uhrf1，十一十一易位（TET）酶（Tet1，Tet2，和低水平Tet3），Aicda和Tdg，這些都以相似的轉(zhuǎn)錄水平表達(dá)涉及在活性DNA去甲基化中。相反，從頭甲基酶Dnmt3a / b及其調(diào)節(jié)劑Dnmt3在2i狀態(tài)下被抑制（圖1B）。為了進(jìn)一步了解從血清到2i ESC的轉(zhuǎn)變動(dòng)力學(xué)，作者在幾個(gè)時(shí)間點(diǎn)確定了它們的DNA甲基化狀態(tài)。首先，通過(guò)液相色譜-質(zhì)譜（LC-MS）（圖1C）以及還原代表亞硫酸氫鹽測(cè)序（RRBS）（圖1D）等方法量化5-甲基胞嘧啶（5mC）和5-羥甲基胞嘧啶（5hmC）的總體水平，DNA的去甲基化在培養(yǎng)基替換后（32小時(shí);約兩輪復(fù)制）迅速發(fā)生，然后逐漸繼續(xù)，達(dá)到穩(wěn)態(tài)14天后水平（圖1C）。而在72小時(shí)內(nèi)觀察到5hmC水平的適度增加表明存在TET活性。

圖1 小鼠胚胎干細(xì)胞血清轉(zhuǎn)2i轉(zhuǎn)化過(guò)程中5mC和5hmC的動(dòng)態(tài)調(diào)節(jié)

（A）胞嘧啶甲基化/去甲基化循環(huán)的示意圖。指出了不同形式的修飾的CpG二元體和相應(yīng)的方法。

（B）在血清到2i轉(zhuǎn)變期間參與DNA甲基化機(jī)制的基因的表達(dá)水平。

（C）在血清至2i轉(zhuǎn)化期間通過(guò)LC MS在E14中測(cè)量的5mC（上）和5hmC（下）的百分比。

（D）在血清至2i轉(zhuǎn)化期間通過(guò)EBS中的RRBS測(cè)量的5mC甲基化的百分比。 水平條代表中值。

（E）數(shù)學(xué)建模的圖形表示。針對(duì)從頭甲基化（p1）和維持甲基化（p2）的不同值繪制最小均方誤差（MMSE）。

（F）在血清到2i轉(zhuǎn)化期間，在E14中從LC-MS獲得的實(shí)際測(cè)量（紅線）覆蓋數(shù)學(xué)模型預(yù)測(cè)（虛線）。該表總結(jié)了數(shù)學(xué)建模的結(jié)果，顯示了血清穩(wěn)態(tài)和2i中p1，p2和p3的估計(jì)值。p1，每平均細(xì)胞分裂成為半甲基化的未甲基化CpG的比例; p2，每平均細(xì)胞分裂完全甲基化的半甲基化CpG的比例;和p3,每平均細(xì)胞分裂變成羥甲基化的甲基化CpG的比例

    為了剖析這三種途徑的作用和相對(duì)貢獻(xiàn)以及所涉及的各種調(diào)節(jié)因子，作者使用數(shù)學(xué)模型來(lái)預(yù)測(cè)整個(gè)時(shí)間過(guò)程中的DNA去甲基化。計(jì)算了完全甲基化的CpG二元組（mCpG/GpCm），半甲基化的CpG二元組（mCpG/GpC）和未甲基化的CpG二元組（CpG/GpC）的百分比以及羥甲基化的CpGs的水平。來(lái)自TAB-seq數(shù)據(jù)和RRBS測(cè)序，這些輸入值以及來(lái)自LC-MS的全局5mC值用于估計(jì)以下參數(shù)，這些參數(shù)與酶豐度和/或活性成正比并反映轉(zhuǎn)化為產(chǎn)物的底物量：p1，動(dòng)態(tài)從頭甲基化的比例值; p2，維持甲基化的比例常數(shù);和p3，活性去甲基化（羥甲基化）的比例常數(shù)。通過(guò)將數(shù)學(xué)模型擬合到5mC數(shù)據(jù)的若干次迭代，能夠估計(jì)具有最低最小均方誤差（MMSE）的E14ESC中血清到2i轉(zhuǎn)變的常數(shù)值（圖1E）。p1，p2和p3反映了三種途徑對(duì)觀察到的DNA甲基化動(dòng)力學(xué)的個(gè)體活性和總體貢獻(xiàn)，并預(yù)測(cè)維持甲基化顯著受損并且觀察到DNA去甲基化的主要驅(qū)動(dòng)因素（圖1F）。

    為了驗(yàn)證并充分了解單個(gè)DNA甲基化和去甲基化酶在全基因組表觀遺傳重編程中對(duì)血清向2i ESC過(guò)渡的貢獻(xiàn)，作者檢測(cè)了對(duì)各DNA甲基化機(jī)制的成分進(jìn)行敲除后的小鼠胚胎干細(xì)胞在血清中幾個(gè)時(shí)間點(diǎn)和從血清到2i ESC過(guò)渡期間的DNA甲基化狀態(tài)（圖2A），并將其與模型的預(yù)測(cè)進(jìn)行比較。在血清→2i轉(zhuǎn)換過(guò)程特別比較了缺乏Dnmt1或Uhrf1的胚胎干細(xì)胞與對(duì)照胚胎干細(xì)胞的去甲基化率，并觀察到去甲基化率增加（圖2B），表明DNA甲基化維持的喪失導(dǎo)致去甲基化率增加。這支持了數(shù)學(xué)模型（彩色框中的虛線）的預(yù)測(cè)，并暗示了2i ESC中DNA甲基化維持的丟失，盡管不是完全丟失。接下來(lái)，比較了缺乏Dnmt3a和Dnmt3b的胚胎干細(xì)胞的去甲基化動(dòng)力學(xué)，血清培養(yǎng)的胚胎干細(xì)胞中Dnmt3a/b的缺失導(dǎo)致5mC的基因組水平僅略微下降（圖2C），并且DNA去甲基化的動(dòng)力學(xué)未改變。在血清到2i的轉(zhuǎn)化中的檢測(cè)表明從頭甲基化的缺失不是DNA甲基化全局喪失的原因（圖2D）。最后，作者評(píng)估了參與活性去甲基化途徑的酶在血清到2i轉(zhuǎn)化中的貢獻(xiàn)。如模型所預(yù)測(cè)的，缺乏Tet1 / 2/3，Tdg或Aicda的ESC與其野生型對(duì)照顯示出非常相似的去甲基化動(dòng)力學(xué)（圖2E）。由于先前已經(jīng)提出TET驅(qū)動(dòng)的氧化作為血清轉(zhuǎn)2i轉(zhuǎn)化中去甲基化的潛在驅(qū)動(dòng)因素，作者證實(shí)了Tet1 / 2/3敲除（KO）細(xì)胞中5hmC的損失（圖2F）。這表明TET酶在血清到2i轉(zhuǎn)化過(guò)程中活躍地氧化5mC，但對(duì)于全局DNA去甲基化來(lái)說(shuō)既不充分也不必要。

圖2 DNA甲基化機(jī)制突變體的全球去甲基化動(dòng)力學(xué)

（A）實(shí)驗(yàn)設(shè)計(jì)的示意圖。

（B）通過(guò)質(zhì)譜法測(cè)量的在血清至2i轉(zhuǎn)變期間Uhrf1和Dnmt1KO ESC中5mC的水平。

（C）通過(guò)質(zhì)譜法測(cè)量的血清培養(yǎng)基中缺失誘導(dǎo)后誘導(dǎo)型Dnmt3a/b KO中5mC的水平。

（D）在血清至2i轉(zhuǎn)換期間誘導(dǎo)型Dnmt3a/b KO中5mC的水平。

（E）通過(guò)質(zhì)譜法測(cè)量的血清至2i轉(zhuǎn)變期間Tet1/2/3 KO和對(duì)照中5mC的水平。

（F）通過(guò)質(zhì)譜法測(cè)量的在血清至2i轉(zhuǎn)變期間Tet1/2/3 KO中和對(duì)照的5hmC水平。

    有研究表明血清-2i轉(zhuǎn)換期間5hmC水平升高以及維生素C（vitC）處理的2i胚胎干細(xì)胞中進(jìn)一步顯示TET依賴性低甲基化，提出了TET蛋白可能有助于去甲基化的可能性。作者在血清到2i轉(zhuǎn)化期間的幾個(gè)時(shí)間點(diǎn)測(cè)定了vitC存在或不存在的情況下 5mC和5hmC的水平（圖3A），并觀察到vitC處理后5hmC顯著增加，這導(dǎo)致加速率去甲基化和進(jìn)一步的全局低甲基化。為了證實(shí)這種效應(yīng)是TET依賴性的，檢測(cè)了在vitC存在或不存在的情況下Tet1/2 KO細(xì)胞中5mC和5hmC的水平（圖3B），表明去甲基化（和羥甲基化）的增加速率取決于TET蛋白的活性。隨后分析了5mC和5hmC在功能不同的基因組區(qū)域的分布（圖3C），血清中富含5hmC，并將啟動(dòng)子分為高、中、低CpG密度啟動(dòng)子（分別為HCP，ICP和LCP），在沒(méi)有vitC的情況下，5mC/5hmC水平相對(duì)于任何三類啟動(dòng)子的變化很小。然而，在存在vitC的情況下，早在4小時(shí)就已經(jīng)明顯轉(zhuǎn)化5mC至5hm，并且5mC損失的速度與CpG密度相關(guān)并且對(duì)于每個(gè)類別是不同的。增強(qiáng)子，定義為重疊H3K4me1，H3K27ac和DNaseI超敏反應(yīng)的元素，在血清胚胎干細(xì)胞中顯示10％-40％5mC水平，與之前的研究結(jié)果一致，增強(qiáng)子的5mC豐度低，并且遵循類似于LCP的動(dòng)力學(xué)?；蚪M的最大部分表現(xiàn)與LCP相似;轉(zhuǎn)化為5hmC和隨后的擦除非常緩慢（圖3C）。最后，分析了Tet1/2/3 KO細(xì)胞中的甲基化,使用k-means聚類對(duì)基因組中平均DNA甲基化分析（圖3D），這些區(qū)域在Tet1/2/3 KO細(xì)胞中具有更高的5mC水平，但無(wú)法鑒定與他們相關(guān)的任何顯著功能的富集。

圖3 血清-2i胚胎干細(xì)胞重編程過(guò)程中TET依賴的去甲基化動(dòng)力學(xué)

（A）在不存在（藍(lán)線）或存在vitC（紅線）的血清至2i轉(zhuǎn)變期間，通過(guò)LC-MS在E14中測(cè)量的5mC（上）和5hmC（下）的百分比。

（B）在不存在（灰線）或存在vitC（綠線）的情況下，在血清到2i轉(zhuǎn)變期間通過(guò)LC-MS在Tet1/2KO中測(cè)量的5mC（上）和5hmC（下）的百分比。

（C）通過(guò)WGBS和TAB-seq在高CpG下測(cè)量的從血清（灰色）到2i（綠色）或2i+vitC（橙色）重編程的前32小時(shí)期間5mC（上）和5hmC水平（下）的高CpG啟動(dòng)子（HCP），中CpG啟動(dòng)子（ICP），低CpG啟動(dòng)子（LCP）和增強(qiáng)子。對(duì)于每類基因組元件，僅考慮在血清中具有5hmC富集的子集。

（D）Tet1/2/3 KO細(xì)胞和相應(yīng)對(duì)照細(xì)胞（Tet WT）中的平均甲基化水平。使用k均值聚類鑒定Tet依賴性區(qū)域。

    通過(guò)數(shù)學(xué)模型預(yù)測(cè)維持甲基化顯著受損并且觀察到DNA去甲基化的主要驅(qū)動(dòng)因素（圖1F）。隨后，證實(shí)了Dnmt1或Uhrf1缺失后的快速去甲基化（圖2B），并且還顯示從頭甲基化的喪失不能解釋觀察到的去甲基化動(dòng)力學(xué)（圖2C和2D）。為了理解這種損傷的機(jī)制調(diào)節(jié)，作者關(guān)注維持性甲基化機(jī)制的各個(gè)組成部分，特別是Uhrf1的作用。據(jù)報(bào)道在原始生殖細(xì)胞中，維持甲基化部分受到UHRF1核排除的影響。作者分析了UHRF1和DNMT1在血清和2i中生長(zhǎng)的ESC中的亞細(xì)胞定位（圖4A），沒(méi)有檢測(cè)到ESC中UHRF1或DNMT1的任何核排斥。然而，觀察到2i樣品中UHRF1蛋白的信號(hào)強(qiáng)度降低。在對(duì)血清到2i轉(zhuǎn)換的細(xì)胞的初始轉(zhuǎn)錄組學(xué)分析顯示Dnmt3a和Dnmt3b以及催化失活的調(diào)節(jié)同種型Dnmt3l的表達(dá)在2i中顯著降低，而Uhrf1 mRNA水平未改變（圖1B）。相反，在2i ESC中，UHRF1蛋白水平顯著降低（通過(guò)蛋白質(zhì)印跡3倍;通過(guò)定量質(zhì)譜法2倍）（圖4B-4D），同時(shí)發(fā)現(xiàn)UHRF1在血清和2i ESC中的異質(zhì)表達(dá)（圖4C），可歸因于UHRF1的細(xì)胞周期依賴性調(diào)節(jié)。為了進(jìn)一步證實(shí)UHRF1在蛋白質(zhì)水平受到調(diào)節(jié)的觀察，我們構(gòu)建了具有UHRF1-GFP融合蛋白的組成型過(guò)表達(dá)的ESC系，與對(duì)內(nèi)源性UHRF1的觀察相似，UHRF1-GFP細(xì)胞系的熒光激活細(xì)胞分選（FACS）分析顯示UHRF1-GFP在血清ESC中表達(dá)，但在血清至2i轉(zhuǎn)化后迅速喪失（圖4E）。

    2i條件下UHRF1蛋白的快速但不完全的丟失，那么其是否損害了復(fù)制灶的DNA維持甲基化機(jī)制的募集？在缺乏Uhrf1的細(xì)胞中的先前實(shí)驗(yàn)證明，UHMF1的存在對(duì)于將DNMT1募集到復(fù)制灶是絕對(duì)必需的。作者通過(guò)在UHRF1和EdU的共染色計(jì)算具有EdU陽(yáng)性病灶的細(xì)胞數(shù)量，用于UHRF1與復(fù)制灶的共定位。雖然UHRF1與血清條件下的大多數(shù)細(xì)胞中的EdU陽(yáng)性復(fù)制灶共定位，但在2i條件下該數(shù)量顯著降低（圖4F），表明血清到2i轉(zhuǎn)化中的DNA去甲基化是由于UHRF1將DNMT1募集到復(fù)制叉受損而導(dǎo)致DNA甲基化維持受損的結(jié)果。

圖4 UHRF1在2i中的蛋白質(zhì)水平受到調(diào)節(jié)

（A）UHRF1（綠色）和DNMT1（紅色）的細(xì)胞定位。DAPI（藍(lán)色）標(biāo)記核。

（B）Western印跡分析。

（C）血清和LT 2i ESC中UHRF1（綠色），DNMT1（紅色）和DAPI（藍(lán)色）的免疫熒光染色。

（D）相對(duì)于血清ESC，在血清至-2i轉(zhuǎn)化期間的不同時(shí)間點(diǎn)通過(guò)蛋白質(zhì)質(zhì)譜法（無(wú)標(biāo)記定量[LFQ]）檢測(cè)的UHRF1的蛋白質(zhì)水平。

（E）血清-2i轉(zhuǎn)化期間ESC中UHRF1-GFP融合蛋白的FACS分析。直方圖顯示轉(zhuǎn)換期間不同時(shí)間點(diǎn)的GFP信號(hào)強(qiáng)度。黑色箭頭表示用于量化GFP +細(xì)胞百分比的閾值，如折線圖所示。

（F）在血清，24小時(shí)2i和LT2i ESC中的共定位的復(fù)制灶的百分比。來(lái)自血清ESC的圖像顯示兩種定量狀態(tài)的實(shí)例;箭頭表示EdU +復(fù)制焦點(diǎn)

    UHRF1的募集依賴于通過(guò)其SRA結(jié)構(gòu)域識(shí)別半甲基化的CpG和通過(guò)其Tudor結(jié)構(gòu)域識(shí)別組蛋白3賴氨酸9的甲基化（H3K9me2）上。H3K9甲基轉(zhuǎn)移酶G9a的缺失導(dǎo)致全球DNA甲基化的喪失，H3K9甲基化的喪失特別是H3K9me2，是導(dǎo)致全基因組去甲基化的另一關(guān)鍵步驟。檢測(cè)發(fā)現(xiàn)血清到2i轉(zhuǎn)化后H3K9me2水平顯著降低（2倍減少）（圖5A，5B）。有趣的是，H3K9me2的丟失很快，并且檢測(cè)到H3K9me2的全局水平與血清ESC中的DNA甲基化之間的IF之間的弱相關(guān)性（圖5C），這與H3K9me2依賴性UHRF1募集在其中起關(guān)鍵作用的觀點(diǎn)一致。為了解H3K9me2丟失的基因座特異性影響，進(jìn)行了H3K9me2的ChIP-seq分析。為了將H3K9me2的存在與5mC的水平相關(guān)聯(lián)，進(jìn)行H3K9me2染色質(zhì)免疫沉淀，然后進(jìn)行亞硫酸氫鹽測(cè)序（ChIP-BS-seq）。與相應(yīng)的輸入樣品相比，H3K9me2結(jié)合的DNA中血清ESC和2i ESC中的5mC水平顯著增加（圖5D）。將基因組分成具有高或低H3K9me2富集的區(qū)域，發(fā)現(xiàn)H3K9me2高區(qū)域的5mC水平顯著增加（圖5E）。接下來(lái)在Uhrf1與相應(yīng)的野生型對(duì)照ESC的RRBS文庫(kù)在H3K9me2區(qū)域上測(cè)量5mC水平發(fā)現(xiàn)這些區(qū)域以UHRF1依賴性方式在2i ESC中保留高水平的DNA甲基化（圖5F）。血清到2i轉(zhuǎn)化過(guò)程中H3K9me2水平是如何調(diào)節(jié)的？作者通過(guò)質(zhì)譜法的蛋白質(zhì)水平（圖5G）量化關(guān)鍵H3K9甲基化酶和去甲基化酶的mRNA水平，發(fā)現(xiàn)幾種H3K9去甲基化酶，包括KDM3A，KDM3B和KDM4A，在2i條件下特異性地在蛋白質(zhì)水平上調(diào)。相反，許多H3K9甲基化酶，例如EHMT1（GLP）和EHMT2（G9A）被下調(diào)。

圖5.H2K9me2在2i ESC中迅速減少

（A）血清和2i ESC中H3K9me2的Western印跡分析。

（B）血清和2i ESC中H3K9me2（綠色），5mC（紅色）和DAPI（藍(lán)色）的免疫熒光染色。

（C）血清，24小時(shí)2i和LT 2i ESC的免疫熒光染色中H3K9me2和5mC的定量信號(hào)的相關(guān)性。

（D）通過(guò)H3K9me2-ChIP-BS-seq測(cè)量H3K9me2結(jié)合的DNA中的5mC水平

（E）在血清和2i ESC中具有高或低H3K9me富集的區(qū)域中的5mC水平。

（F）在E14，Uhrf1 WT和Uhrf1 KO ESC中的血清-2i轉(zhuǎn)變期間由RRBS在富含H3K9me2的區(qū)域測(cè)定的5mC的水平。

（G）在血清至2i轉(zhuǎn)化期間的不同時(shí)間點(diǎn)通過(guò)蛋白質(zhì)質(zhì)譜法檢測(cè)的已知H3K9修飾物的相對(duì)蛋白質(zhì)水平。