一文讀懂全基因組甲基化測序技術(shù)原理及應用

 

技術(shù)背景

DNA甲基化是一種在原核和真核生物基因組中常見(jiàn)的復制后表觀(guān)遺傳修飾，可以在不改變 DNA分子一級結構的情況下調節基因組的功能，其在基因表達、轉座子沉默、染色質(zhì)相互作用、細胞分化以及生長(cháng)發(fā)育過(guò)程中起著(zhù)重要作用。

其形成過(guò)程主要是在DNA甲基轉移酶的作用下，將S-腺苷甲硫氨酸分子上的甲基(CH3-)轉移到DNA的CpG兩個(gè)核苷酸的胞嘧啶上進(jìn)行DNA 甲基化修飾，形成 5-甲基胞嘧啶 (5-mC)、N6-甲基腺嘌呤 (N6-mA) 以及7-甲基鳥(niǎo)嘌呤(7-mG)。

DNA 甲基化造成基因沉默的原因大多都是由于轉錄受到抑制從而影響了目的基因的表達：（1）啟動(dòng)子區域發(fā)生甲基化 ；（2）基因序列中某位點(diǎn)發(fā)生甲基化抑制蛋白的表達 ；（3）染色體結構由于受到 DNA 序列的變化從而變得緊密。

因此，闡明DNA甲基化的動(dòng)力學(xué)和遺傳學(xué)及其分子調控機制，可在基因表達層面上探究動(dòng)植物的生長(cháng)發(fā)育及抗逆性機理，為優(yōu)良新品種的培育和遺傳改良研究提供更多理論依據。

目前應用廣泛的全基因組甲基化檢測方式主要是BS-seq：用亞硫酸氫鹽處理DNA將胞嘧啶殘基（C）轉化為尿嘧啶（U），而5-甲基胞嘧啶殘基（5mC）不會(huì )發(fā)生轉變。因此可通過(guò)測序的手段識別甲基化位點(diǎn)。但是亞硫酸氫鹽處理后DNA會(huì )發(fā)生降解，測序長(cháng)度受限，使得較大的區域難易被解析。

三代測序平臺則可以在不依賴(lài)DNA實(shí)驗處理的情況下，根據信號差來(lái)識別甲基化位點(diǎn)信息。其中，PacBio 已經(jīng)實(shí)現通過(guò)HiFi測序直接檢測 DNA 樣本中 CpG 位點(diǎn)的 5mC甲基化?；贖iFi 的 5mC 檢測能夠在提供既長(cháng)又準的測序結果的同時(shí)，給出準確的單倍型甲基化信息。同時(shí)，具有甲基化信息的HiFi數據，所占存儲僅增~5%，極大的壓縮了可分析的空間，為表觀(guān)遺傳挖掘保駕護航！

技術(shù)原理

PacBio平臺主要通過(guò)DNA聚合酶與模板鏈結合，在鏈合成的過(guò)程中，不同熒光標記的堿基會(huì )釋放不同的熒光信號，根據光的波長(cháng)與峰值從而判斷堿基類(lèi)型與順序。聚合酶合成中每個(gè)堿基間都存在時(shí)間間隔，當模板堿基有修飾時(shí)，兩個(gè)相鄰的脈沖峰之間的距離就會(huì )延長(cháng)。當IPD比例（兩個(gè)相鄰的脈沖峰之間的距離和參考序列距離的比值）明顯大于1 時(shí)，可推斷這個(gè)位置存在修飾。

技術(shù)優(yōu)勢

 

分析方案

測序深度；≥10x。根據數據檢測圖可以看到，當深度達到10X，即可達到很好的檢出效果。

應用場(chǎng)景

（1）全基因組甲基化

對單個(gè)樣本進(jìn)行HiFi測序可檢測整個(gè)基因組的甲基化模式，如轉錄起始位點(diǎn)的低甲基化。不同組織部位甲基化差異，構建全基因組的甲基化圖譜。

（2）泛基因組甲基化

通過(guò)對多個(gè)樣本進(jìn)行測序來(lái)確定物種品系間的差異甲基化，從而揭示品系間功能差異的調控機制。

本研究構建了青鳉（Oryzias latipes）近親繁殖系中12尾材料的基因組，并通過(guò)比對構建了青鳉的圖型泛基因組。結果顯示泛基因組總長(cháng)度為1.3 Gb，大約是其參考基因組的1.8倍。通過(guò)泛基因組展示了更大、更復雜的結構變異-“暗變異”，為進(jìn)一步研究結構變異（SV）如何影響感興趣的表型提供了重要的研究基礎。作者通過(guò)構建基因組時(shí)的三代測序數據，進(jìn)一步對12尾近親繁殖系的大腦樣本中的CpG甲基化進(jìn)行了差異研究。結果顯示，甲基化模式在多代中表現一致，說(shuō)明甲基化是青鳉魚(yú)的一種遺傳性狀，能夠持續多代，并可能影響表型性狀。最后結合RNA等數據，表明galm基因甲基化可能性與肝臟表達之間存在顯著(zhù)相關(guān)性。

（3）單體型甲基化

HiFi測序能夠同時(shí)分階段讀取母系和父系單倍型，并檢測甲基化。這直接揭示了印記的模式，其中染色體的甲基化狀態(tài)取決于它是從母親還是父親遺傳的。

日本田中魚(yú)（Oryzias latipes）作為脊椎動(dòng)物模式生物有著(zhù)悠久的歷史。我們利用其異常高的近親繁殖耐受性建立了Medaka近親繁殖Kiyosu Karlsruhe（MIKK）小組：來(lái)自單個(gè)野生種群的近等基因脊椎動(dòng)物系的集合[3]。

在與這篇文章一起發(fā)表的配套文章中，我們根據Illumina短讀序列與最近、完全組裝的參考基因組（日本南部medaka自交系HdrR）的比對，提供了80個(gè)個(gè)體MIKK面板系[4]的詳細遺傳特征。

盡管這使我們能夠發(fā)現與HdrR相關(guān)的MIKK面板中的許多遺傳變異，但該方法不可避免地隱藏了某些變異，包括更大、更復雜的結構變異-“暗變異”-這可能會(huì )對每個(gè)品系產(chǎn)生功能性后果。在這里，我們描述了我們如何使用牛津納米孔技術(shù)（ONT）長(cháng)讀測序來(lái)揭示12個(gè)MIKK面板系中的一些暗變異，從而使我們對其基因組變異進(jìn)行更完整的評估，并為進(jìn)一步研究結構變異（SV）如何影響感興趣的表型鋪平了道路。

 

DNA甲基化分析：

• MIKK面板樣本中所有顯著(zhù)CpG島差異甲基化的熱圖。
• CpG島按基因組位置和X軸排序。樣本通過(guò)Y軸上的分層聚類(lèi)排序，并進(jìn)行顏色編碼，以便用相同的顏色表示兄弟線(xiàn)。
• 顏色方案根據來(lái)自的中值對數似然比的值?3（藍色）到3（紅色）之間的值不明確?1和1為白色。
• B紅色顯著(zhù)差異甲基化CpG（n=4249）和藍色所有非顯著(zhù)島（n=8727）的最接近基因TSS的距離分布。
• C 24條HdrR染色體基因組窗口中顯著(zhù)CpG島的數量，從0（白色）到6（深紅色）。
• D示例100kb區域包含紅色的顯著(zhù)CpG島（chr15:1565040-1565987），而非白色的非顯著(zhù)CpG島。
• E面板D中突出顯示的CpG島的CpG水平對數似然比核密度圖。通過(guò)降低中值llr，在Y軸上對樣本進(jìn)行排序。單個(gè)CpG值由點(diǎn)表示。
• 圖D中突出顯示的CpG島的對數似然比熱圖，按樣本分層聚類(lèi)。X軸上是按基因組位置排序的單個(gè)CpG

 

如果您對我們提供的服務(wù)感興趣，歡迎聯(lián)系我們。

聯(lián)系我們