孟昶 朱磊 田雪文 王清路

摘 要：DNA甲基化調控基因的表達，在肥胖和2型糖尿病等代謝類疾病的發(fā)病機制中起重要作用。針對 DNA 甲基化檢測研究方法包括全基因組亞硫酸鹽測序（BS-seq）、甲基化 DNA 免疫共沉淀測序（Me DIP-seq）和亞硫酸氫鹽測序（RRBS）。研究表明低氧與EPAS1甲基化呈負相關，與LINE-1甲基化呈正相關，低氧運動可以使CpG和DMR甲基化顯著。綜述DNA甲基化分子機制及相關技術，探討低氧運動與DNA甲基化的分子機制及生理學意義，展望低氧運動與DNA甲基化的研究前景，為防治代謝類疾病提供更多依據(jù)。

關鍵詞：低氧運動;DNA甲基化;亞硫酸氫鹽測序

中圖分類號：G804.7? 文獻標識碼：A? 文章編號：1009-9840（2019）05-0066-04

Abstract：DNA methylation is one of the well-known epigenetic markers that regulate gene expression and play an important role in the pathogenesis of metabolic diseases such as obesity and type 2 diabetes. Methods for DNA methylation detection include genome-wide sulfite sequencing （BS-seq）， methylated DNA immunoprecipitation sequencing （Me DIP-seq）， and bisulfite sequencing （RRBS）. Studies have shown that hypoxia is negatively correlated with EPAS1 methylation and positively correlated with LINE-1 methylation. In this paper， we review the molecular mechanism and physiological significance of hypoxia and DNA methylation by reviewing the molecular mechanisms and related technologies of DNA methylation， and prospectively study the hypoxia movement and DNA methylation research， and provide more for the prevention and treatment of metabolic diseases.

Key words：hypoxic exercise;DNA methylation;RRBS

隨著社會的進步，人民生活水平的提高，生活習慣和生活環(huán)境發(fā)生了翻天覆地的變化，伴隨而來的是肥胖、三高癥以及2型糖尿病等代謝類疾病發(fā)病率不斷提高，并且患病人群具有年輕化趨勢，嚴重威脅人類的健康。表觀遺傳學被認為是基因表達的可遺傳變化，而不是由于DNA序列的任何改變[1]。表觀遺傳學的重要性在于對一些現(xiàn)象提供了部分的解釋，而這些現(xiàn)象單靠經典遺傳學是無法解釋的，因此它已引起越來越多的關注[2]。DNA甲基化是最著名的表觀遺傳學標記之一，在基因表達修飾中起到不可忽視的作用。以往研究證實了低氧條件下對EPAS1甲基化呈負相關，與LINE-1甲基化呈正相關，而現(xiàn)階段初步研究發(fā)現(xiàn)，低氧運動與常氧運動相比甲基化升高或降低大于25%的單個CpG甲基化位點共計675 個，所有 CpG 平均甲基化，這很可能是調控肥胖等代謝性疾病的一種重要干預手段[3-4]。本文通過綜述DNA甲基化分子機制及相關技術，探討低氧運動與DNA甲基化的分子機制及生理學意義，展望研究前景，為運動有益于健康提供更多重要依據(jù)。

1 DNA甲基化

DNA甲基化是一種DNA化學修飾形式，是將甲基加入胞嘧啶或腺嘌呤的過程，還可以在不改變序列的情況下改變DNA分子活性。5-甲基胞嘧啶的甲基化在原核生物與真核生物中是廣泛存在的，它是一種非常重要的表觀遺傳修飾，可以調控基因的活性，影響基因組印跡、細胞分化、轉錄調控、染色質重塑等一系列關鍵過程。通過基因組分析DNA甲基化對理解甲基化在正常生物學和代謝類疾病中的影響至關重要。有研究表明，哺乳動物基因組中有60%的CpG二核苷酸是甲基化的[5]。CpG二核苷酸通常存在于CpG中，通常是甲基化的基因啟動子[6]。而近期在哺乳動物基因組中發(fā)現(xiàn)5-羥甲基胞嘧啶、5-甲酰胞嘧啶、5-羧胞嘧啶等5-甲基胞嘧啶氧化衍生物，進一步拓展了人們對甲基化調控的認識[7]。

DNA甲基化是由DNA甲基轉移酶引發(fā)的。DNA甲基轉移酶-1（DNMT1）優(yōu)先在半甲基化DNA位點進行甲基化，并且在細胞復制過程中保持甲基化模式[8]。甲基轉移酶與DNA結合，從DNA螺旋中“翻轉”出胞嘧啶，并將來自S-腺苷甲硫氨酸的甲基連接到胞嘧啶的5'位置[9-10]。DNMT的N-末端調節(jié)結構域包含蛋白質與DNA相互作用，蛋白質與蛋白質相互作用和核定位的基序。C末端區(qū)域含有限定酶活性位點的保守基序。DNMT3L的N-末端區(qū)域與DNMT3A和DNMT3B的N-末端區(qū)域具有相似性，但在其C-末端結構域中缺乏催化活性位點[11]。DNA甲基化的特征是在胞嘧啶嘧啶環(huán)的C5位或腺嘌呤環(huán)的N6位上添加一個甲基。雖然DNA甲基化是在DNA的結構被分解之前發(fā)現(xiàn)的，但對其在生物體中的多種功能的研究仍在繼續(xù)[12]。后者可以歸因于與DNA甲基化相關的過程的復雜性，以及DNA甲基化很難分析的事實[13]。一致地，DNA甲基化可以表示為表觀遺傳修飾，因為它將不包含在DNA序列本身中的相關信息傳遞給下一代。在人類中，甲基化僅限于胞嘧啶殘基，主要存在于CpG中。這種二核苷酸在整個基因組中都沒有得到充分的表達，但是在0.33 kb的片段中，這種二核苷酸的表達率卻在增加，這就是所謂的CpG島[14]。約40%的人類基因在啟動子區(qū)中含有CpG島，其甲基化通常與沉默表達有關[15]。DNA甲基化是由DNA甲基轉移酶（DNMTs）建立和維持的，DNMT3A和DNMT3B介導新的甲基化，DNMT1維持已建立的DNA甲基化。DNA甲基化在人類中具有調控特定的基因[16]：1）特定基因的調控。2）基因組印記。3）X染色體失活。4）基因組防御。

2 DNA甲基化檢測技術

DNA甲基化檢測研究方法主要包括BS-seq、MeDIP-seq和RRBS三種技術，其中最經濟有效的方法是RRBS技術。

2.1 BS-seq

BS-Seq是在單核苷酸分辨率下分析DNA胞嘧啶甲基化基因組的首選方法，被譽為是全基因組測序的“金標準”[17]。BS-Seq技術的基礎是亞硫酸氫鹽可將未甲基化的胞嘧啶脫胺成尿嘧啶，而甲基化的胞嘧啶對治療具有耐藥性。通過對亞硫酸氫鹽轉化的基因組DNA進行深度測序，可以測量基因組中每個可映射胞嘧啶位置的甲基化水平[18]。BS-Seq是單核苷酸分辨率下甲基化狀態(tài)分析最有效的方法。由于Sanger測序和大多數(shù)第二代測序方法不能區(qū)分甲基化和非甲基化胞嘧啶，測序前需要亞硫酸氫鹽處理將非甲基化胞嘧啶轉化為尿嘧啶。然而，像PacBio這樣的第三代測序能夠根據(jù)聚合酶遇到修飾堿基時的停頓時間來檢測DNA修飾，因此最終可能使BS-Seq方法過時[19]。

2.2 MeDIP-seq

MeDIP-seq技術由Weber于2005年首次提出[20]，這是基于抗體富集原理進行的全基因組的甲基化檢測技術，該方法是將甲基化 DNA 免疫沉淀和高通量 DNA 測序技術相結合，即采用甲基化 DNA 免疫共沉淀技術，利用5′-甲基胞嘧啶抗體特異性富集基因組上已被甲基化的 DNA 片段，然后采用高通量測序，獲得基因組DNA甲基化發(fā)生區(qū)域及甲基化水平等[21]。該方法的優(yōu)點在于精確度較高，覆蓋范圍廣，不需要較多的測序通量，并且成本較低，但只能確定一段DNA 區(qū)域的甲基化程度[22]。

2.3 RRBS

亞硫酸氫鹽處理DNA可使未甲基化的胞嘧啶脫氨成尿嘧啶，不受甲基化胞嘧啶的影響[23]，從而提供了一種甲基依賴的DNA突變，這種技術被稱為亞硫酸氫鹽測序（RRBS）[24]。RRBS方法是在Rudolf Jaenisch的實驗室中開發(fā)的，由于其成本效益和基因組中CpG位點的相對高覆蓋率而變得流行[25]。該方法特別適用于細胞數(shù)量有限的樣品，如生殖細胞或卵母細胞。RRBS方法基于用MspI（5個CCGG識別位點）進行DNA消化，然后進行末端修復，接頭連接，DNA亞硫酸氫鹽轉化和文庫擴增。該方法提供了位于兩個MspI位點之間的DNA片段的信息，其距離350 bp，可以代表人類基因組中大約82%的CpG島和大約70%的啟動子。近年來，以12個細胞和單個細胞為起始材料，優(yōu)化了DNA甲基化組譜的還原表征RRBS方法，以研究胚胎發(fā)育過程中的甲基化組景觀[26]。然而，為了從DC或SC樣品中制備足夠數(shù)量的下一代測序文庫，進行了過量的PCR擴增（>40個循環(huán)），這可能導致大量的數(shù)據(jù)重復。由于RRBS的測序庫是基于酶消化的，酶消化可以從同一位置染色體的多個拷貝中產生相同的序列，因此沒有任何可用的策略可以區(qū)分這些序列是來自相同片段的不同拷貝還是來自PCR的重復誘導[27-28]。

3 低氧運動對DNA甲基化的影響

人類已適應高海拔生活的三個主要地理區(qū)域包括青藏高原、埃塞俄比亞高原和安第斯高原[29]。來自這些區(qū)域中的每個區(qū)域的人群對于低海拔人群和彼此相對地表現(xiàn)出獨特的循環(huán)、呼吸和血液適應性。這些適應性包括動脈血氧飽和度、血紅蛋白濃度、靜息通氣和低氧通氣反應等差異[30]。然而，單靠遺傳學并不能完全解釋所觀察到的低氧環(huán)境適應性變異程度。在高海拔地區(qū)的生長發(fā)育也被證明有助于高海拔地區(qū)的適應性特征。

表觀遺傳調控與遺傳過程密切相關，使生物體為其生存的外部環(huán)境做好準備，并有助于提供一種潛在的機制，使其能夠適應不斷變化的環(huán)境條件[31-32]。DNA甲基化是理解最透徹、研究最廣泛的表觀遺傳修飾。這種修飾最常見的發(fā)生在胞嘧啶之后是鳥嘌呤，其排列方式稱為胞嘧啶-磷酸鹽-鳥嘌呤（GpG）位點[33]。DNA甲基化是DNMT進行的，它選擇性靶向特定的DNA序列，要么從頭開始，要么維持甲基化[34-36]。內皮PAS結構域蛋白1（EPAS1）是一種缺氧激活基因，它在低氧條件下表達較高[37]，它的甲基化水平以前已經被證明與結直腸癌有關，隨著甲基化的增加，EPAS1蛋白的表達降低?？傊?，這些發(fā)現(xiàn)使得EPAS1成為研究高海拔適應性表觀遺傳學的一個自然的先天候選。Ainash等人[17]分析了秘魯海拔4388米和海拔0米的572名克丘亞人的EPAS1啟動子區(qū)和LINE-1重復元件的DNA甲基化。在高海拔地區(qū)招募的參與者EPAS1 DNA甲基化水平較低，而LINE-1甲基化水平較高。出生的海拔與較高的1系甲基化有關，而與EPAS1甲基化無關。生活在高海拔地區(qū)的年數(shù)與EPAS1甲基化呈負相關，與LINE-1甲基化呈正相關。發(fā)現(xiàn)了4個單碳代謝SNP（MTHFD1 rs2236225， TYMS rs502396， FOLH1 rs202676， GLDC rs10975681），它們累積解釋了平均LINE-1甲基化變異的11.29%。并確定了1號線甲基化和全基因組SNP主成分之間的關聯(lián)，當前和終身暴露于高海拔缺氧環(huán)境中都會對EPAS1和LINE-1甲基化產生影響，這表明表觀遺傳修飾可能在高海拔適應性中發(fā)揮作用。而田雪文博士詮釋了低氧運動與DNA甲基化之間的聯(lián)系，通過低氧運動和常氧運動訓練肥胖大鼠，從表觀遺傳學的角度，分析比較常氧運動肥胖大鼠和低氧運動肥胖大鼠Wnt/β-catenin信號通路 DNA甲基化的差異顯著基因。初步證實了低氧運動肥胖大鼠減脂和Wnt/β-catenin 信號通路的關系。比較分析了常氧運動肥胖大鼠和低氧運動肥胖大鼠全基因組 DNA的CpG島發(fā)生甲基化的基因，具有213 個 DMR 差異區(qū)域以及 675 個發(fā)生甲基化具有顯著差異性的基因。但是，目前低氧運動與DNA甲基化的研究還處于起步階段，其分子機制還不明確，還需進一步研究。

4 小結

DNA甲基化是著名的表觀遺傳學標記之一，調控基因的表達，在肥胖和2型糖尿病等代謝類疾病的發(fā)病機制中發(fā)揮重要作用。DNA 甲基化檢測研究方法主要包括BS-seq、MeDIP-seq和RRBS技術，其中最經濟有效的方法是RRBS技術。目前初步證實低氧運動可以使CpG和DMR甲基化顯著，但是低氧運動與DNA甲基化的研究還處于起步階段，還需要進一步研究探討其分子機制。深入研究低氧運動與 DNA 甲基化的關系，為防治代謝類疾病提供更多依據(jù)。

參考文獻：

[1]Kaplow I. M.， Macisaac J. L.， Mah S. M.，et al.A pooling-based approach to mapping genetic variants associated with DNA methylation [J]. Genome Research， 2015， 25（6）： 907.

[2]Isabel Lokody.Epigenetics：Histone methyltransferase mutations promote leukaemia [J].Nature Reviews Cancer，2014，14（4）：214-5.

[3]Dirk Schübeler. Function and information content of DNA methylation [J]. Nature， 2015， 517（7534）： 321-326.

[4]Fabio Coppedè. Epigenetic biomarkers of colorectal cancer： Focus on DNA methylation [J]. Cancer Letters， 2014， 342（2）： 238-247.

[5]Bestor T H， Hellewell S B， Ingram V M. Differentiation of two mouse cell lines is associated with hypomethylation of their genomes [J]. Molecular & Cellular Biology， 1984， 4（9）： 1800-1806.

[6]Okano M.，Bell D W， Haber D A，et al.DNA methyltransferases dnmt3a and dnmt3b are essential for de novo methylation and mammalian development [J]. Cell， 1999， 99（3）： 247-257.

[7]Von Knel T， Huber A. R. DNA methylation analysis [J]. Swiss Medical Weekly， 2013， 143（5）： w13799.

[8]Sturgeon Susan R.， Pilsner J. Richard， Arcaro Kathleen F.，et al.White blood cell DNA methylation and risk of breast cancer in the prostate， lung， colorectal， and ovarian cancer screening trial （plco） [J]. Breast Cancer Research，2017，19（1）：94.

[9]Keil K. P.， Vezina C. M. DNA methylation as a dynamic regulator of development and disease processes： Spotlight on the prostate [J]. Epigenomics， 2015， 7（3）： 413-425.

[10]Kenichiro Hata， Masaki Okano， Hong Lei，et al.Dnmt3l cooperates with the dnmt3 family of de novo DNA methyltransferases to establish maternal imprints in mice [J]. Development， 2002， 129（8）： 1983-1993.

[11]Zhi-Min Zhang， Shuo Liu， Krystal Lin，et al.Crystal structure of human DNA methyltransferase 1 [J]. Journal of Molecular Biology， 2015， 427（15）： 2520-2531.

[12]Wyatt G. R. The purine and pyrimidine composition of deoxypentose nucleic acids [J]. Biochemical Journal， 1951， 48（5）： 584-590.

[13]Jeltsch Albert. Beyond watson and crick： DNA methylation and molecular enzymology of DNA methyltransferases [J]. Chembiochem， 2010， 33（22）.

[14]Josep Casadesús， David Low. Epigenetic gene regulation in the bacterial world [J]. Microbiology & Molecular Biology Reviews Mmbr， 2006， 70（3）： 830.

[15]Mehrnaz Fatemi， Pao Martha M， Shinwu Jeong，et al.Footprinting of mammalian promoters： Use of a cpg DNA methyltransferase revealing nucleosome positions at a single molecule level [J]. Nucleic Acids Research， 2005， 33（20）： e176-e176.

[16]Smith Zachary D， Alexander Meissner. DNA methylation： Roles in mammalian development [J]. Nature Reviews Genetics， 2013， 14（3）： 204-220.

[17]Barrès R， Yan J.， Egan B，et al.Acute exercise remodels promoter methylation in human skeletal muscle [J]. Cell Metabolism， 2012， 15（3）： 405-411.

[18]Chen Yun Ru， Yu Sheng， Zhong Silin. Profiling DNA methylation using bisulfite sequencing （bs-seq） [J]. Methods in Molecular Biology，2018（1675）：31.

[19]Frommer M.，Mcdonald L E， Millar D S，et al.A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands [J]. Proceedings of the National Academy of Sciences of the United States of America， 1992.

[20]Oluwatosin Taiwo， Wilson Gareth A， Tiffany Morris， et al.Methylome analysis using medip-seq with low DNA concentrations [J]. Nature Protocols， 2012， 7（4）： 617-36.

[21]馬浪浪， 梁振娟， 先新，等.基于高通量測序技術檢測全基因組dna甲基化水平的方法 [J]. 生物技術通報， 2015， 31（7）： 45-50.

[22]Jeoung N. H.， Harris R. A. Pyruvate dehydrogenase kinase-4 deficiency lowers blood glucose and improves glucose tolerance in diet-induced obese mice [J]. Am J Physiol Endocrinol Metab， 2008， 295（1）： E46.

[23]David Sims， Ilott Nicholas E， Sansom Stephen N，et al.Cgat： Computational genomics analysis toolkit [J]. Bioinformatics， 2014， 30（9）： 1290-1291.

[24]Yin Danqing， Ritchie Matthew E.， Jabbari Jafar S.， et al.High concordance between illumina hiseq2500 and nextseq500 for reduced representation bisulfite sequencing （rrbs） [J]. Genomics Data， 2016， 10（C）： 97-100.

[25]Alexander Meissner， Andreas Gnirke， Bell George W，et al.Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis [J].Nucleic Acids Research， 2005， 33（18）： 5868-5877.

[26]Hongshan Guo， Ping Zhu， Xinglong Wu，et al.Single-cell methylome landscapes of mouse embryonic stem cells and early embryos analyzed using reduced representation bisulfite sequencing [J]. Genome Research， 2013， 23（12）： 2126-2135.

[27]Seung-Gi Jin， Swati Kadam， Pfeifer Gerd P. Examination of the specificity of DNA methylation profiling techniques towards 5-methylcytosine and 5-hydroxymethylcytosine [J]. Nucleic Acids Research，38，11（2010-4-5）， 2010， 38（11）： e125.

[28]Ivanov M， Ingelman-Sundberg M， Kals M，et al.Targeted profiling of 5-（hydroxy）methylcytosine in genomic DNA from human livers： Next-generation sequencing of target enriched DNA reveals unexpectedly high interindividual variability of cytosine methylation and hydroxymethylation [C].IEEE International Conference on Bioinformatics & Biomedicine，2014.

[29]Moore L. G.Human genetic adaptation to high altitude [J]. High Altitude Medicine & Biology， 2001， 2（2）： 257.

[30]Bigham Abigail W， Julian Colleen G， Wilson Megan J，et al.Maternal prkaa1 and ednra genotypes are associated with birth weight， and prkaa1 with uterine artery diameter and metabolic homeostasis at high altitude [J]. Physiological Genomics， 2014， 46（18）： 687-697.

[31]Jingwen， Yongjuan， Wenhao，et al.The epigenetic switches for neural development and psychiatric disorders [J]. Journal of Genetics & Genomics， 2013， 40（7）： 339-346.

[32]Mohn F， Schübeler D. Genetics and epigenetics： Stability and plasticity during cellular differentiation [J]. Trends in Genetics， 2009， 25（3）： 129-136.

[33]Tollefsbol Trygve O. Advances in epigenetic technology [M]. 2011.

[34]Fraga Mario F， Manel Esteller. DNA methylation： A profile of methods and applications [J]. Biotechniques， 2002， 33（3）： 636-649.

[35]Beall Cynthia M，Cavalleri Gianpiero L， Libin Deng， et al.Natural selection on epas1 （hif2alpha） associated with low hemoglobin concentration in tibetan highlanders [J]. Proc Natl Acad Sci U S A， 2010， 107（25）： 11459-11464.

[36]Matthieu Foll， Gaggiotti Oscar E， Daub Josephine T，et al.Widespread signals of convergent adaptation to high altitude in asia and america [J]. American Journal of Human Genetics， 2014， 95（4）： 394-407.

[37]Nanduri J， Semenza G. L.， Prabhakar N. R. Epigenetic changes by DNA methylation in chronic and intermittent hypoxia [J]. Am J Physiol Lung Cell Mol Physiol， 2017， 313（6）： L1096.