基因甲基化检测方法

㈠ 如何检测 DNA 甲基化

需要去专门的检测机构或者医院检查。

DNA甲基化（DNA methylation）为DNA化学修饰的一种形式，能够在不改变DNA序列的前提下，改变遗传表现。所谓DNA甲基化是指在DNA甲基化转移酶的作用下，在基因组CpG二核苷酸的胞嘧啶5号碳位共价键结合一个甲基基团。大量研究表明，DNA甲基化能引起染色质结构、DNA构象、DNA稳定性及DNA与蛋白质相互作用方式的改变，从而控制基因表达。

机制

由于Dnmtl和Dnmt3基因家族没有针对CpG二核苷酸序列的特异性，人们因此提出了DNA甲基化转移酶发现靶位点的机制。

首先，甲基化转移酶并不是同等地接近所有染色体区域。具有染色体重构和DNA螺旋酶活性的蛋白质能调节哺乳动物细胞内DNA甲基化，如SNF2家族2个成员ATRX和Lsh;其次，附件因子（蛋白质、RNA等）能召集DNA甲基化转移酶到特定基因组序列或染色体结构中，如pRB蛋白等能够与Dnmtl作用，在S期晚期将它召集到高度甲基化的异染色质区。

以上内容参考：网络-DNA甲基化

㈡ 易基因｜一文读懂：十大DNA甲基化研究核心问题

DNA甲基化是最早被发现、也是研究最深入的表观遗传调控机制之一，近年来关于DNA甲基化的研究成果屡屡见刊。我翻阅各类文献，为大家总结了十大DNA甲基化研究核心问题，包括什么是DNA甲基化，DNA甲基化的主要形式、DNA甲基化与去甲基化、植物中的DNA甲基化、DNA甲基化的主要功能、DNA甲基化作为生物标志物的潜力、DNA甲基化的主要研究方向、DNA甲基化检测方法、样本不同如何选择DNA甲基化检测技术、DNA甲基化数据挖掘等，让您一文读懂DNA甲基化。

1、什么是DNA甲基化

DNA甲基化（DNA methylation）为DNA化学修饰的一种形式，是指DNA分子在DNA甲基转移酶的作用下将甲基选择性地添加到特定碱基上的过程。DNA甲基化能够在不改变DNA序列的前提下，改变遗传表现，是最重要的表观遗传调控方式之一。

2、DNA甲基化的主要形式

5-甲基胞嘧啶（5-mC）：最重要的一种DNA甲基化修饰，广泛存在于植物、动物等真核生物基因组中, 称誉为“第五碱基”。

5-羟甲基胞嘧啶世衫腊（5-hmC）：哺乳动物的“第六碱基”。

N6-甲基腺嘌呤（N6-mA）：在细菌、藻类及动植物基因组中存在。

7-甲基鸟嘌呤（7-mG）

3、DNA甲基化与去甲基化（5mC）：

DNA甲基化反应分为2种类型：一种是2条链均未甲基化的DNA被甲基化，称为从头甲基化（denovo methylation）；另一种是双链DNA的其中一条链已存在甲基化，另一条未甲基化的链被甲基化，这种类型称为保留甲基化（maintenance methylation）。

甲基化的DNA可以发生去甲基化。DNA的去甲基化由基因内部的片段及与其结合的因子所调控，包括：

主动去甲基化(Active demethyaltion)：哺乳动物TET酶主动去甲基化，5mC经过TET作用转化成5hmC。

被动物甲基化(Passive demethylation)：DNA通过不断复制丢失/稀释甲基化。

4、植物中的DNA甲基化

对于植物而言，面对生长环境塌颂的改变，表观遗传的变异会改变植物DNA的构象，从而改变染色质和蛋白质的结构，达到调节基因组的作用。研究发现，当植物面临生物胁迫和非生物胁迫时，植物基因组中DNA甲基化会发生改变，并且这些改变会遗传给后代。所以DNA甲基化的改变能够丰富植物物种的多样性，加强植物的环境适应性。

不同于哺乳动物基因组只有CG甲基化，植物基因组甲基化有CG，CHG，CHH（H代表任何非G的碱基）甲基化。而维持这三种不同的DNA甲基化的分子机制非常复杂。

5、DNA甲基化的主要功能

DNA甲基化能引起染色体结构、DNA构象、DNA稳定性及DNA与蛋白质相互作用方式改变，从而控制基因的表达。

保持基因组遗传物质的稳定性（TE的高甲基化）

调控基因的表达（顺式作用元件的动态甲基化，如Promoter/Enhancer等）

DNA甲基化参与基因转录调控、细胞分化、胚胎发育、X染色体失活、基因印记和肿瘤的发生等过程。

6、DNA甲基化作为生物标志物的潜力

相比基因组，DNA甲基化能够反应环境的影响

DNA甲基化不像基因组那么一成不变，也不像转录组和蛋白组那么搜滑不稳定。

DNA甲基化处于动态变化中，能够像年轮一样记录环境因素的影响。

DNA甲基化标志物是最有应用前景的表观遗传标志物 。

7、DNA甲基化的主要研究方向

8、DNA甲基化检测方法

目前研究中常用的DNA甲基化测序方法包括全基因组（WGBS、oxWGBS等）、简化基因组（dRRBS、RRBS、XRBS等）、靶向基因组（液相捕获）、靶向基因（TBS）和850K芯片等，适用于多种不同应用场景。

WGBS和RRBS用于基因组范围内研究探索，并筛选目标基因（候选DNA甲基化标记物）

TBS用于后续目标基因的甲基化验证

9、样本不同，如何选择DNA甲基化检测技术

易基因结合十余年DNA甲基化研究经验，全面总结了不同样本类型对于不同DNA甲基化检测工具选择的不同，技术推荐标准可以分为三点：

最有力方案（金标准）：WGBS/微量WGBS/scWGBS 

最具性价比方案：RRBS/dRRBS/XRBS（动物）

大规模临床转化/目标区域甲基化验证：TBS 

10、DNA甲基化数据挖掘

DNA甲基化一般遵循三个步骤进行数据挖掘。首先，进行整体全基因组甲基化变化的分析，包括平均甲基化水平变化、甲基化水平分布变化、降维分析、聚类分析、相关性分析等。其次，进行甲基化差异水平分析，筛选具体差异基因，包括DMC/DMR/DMG鉴定、DMC/DMR在基因组元件上的分布、DMC/DMR的TF结合分析、时序甲基化数据的分析策略、DMG的功能分析等。最后，将甲基化组学&转录组学关联分析，包括Meta genes整体关联、DMG-DEG对应关联、网络关联等。

㈢ DNA甲基化测序方法介绍

内容来自于  生信人社区  欢迎关注。

DNA 甲基化是表观遗传学(Epigenetics)的重要组成部分，在维持正常细胞功能、遗传印记、胚胎发育以及人类肿瘤发生中起着重要作用，是目前新的研究热点之一。

DNA 甲基化及CpG岛

DNA 甲基化是最早发现的基因表观修饰方式之一,可能存在于所有高等生物中。DNA 甲基化能关闭某些基因的活性，去甲基化则诱导了基因的重新活化和表达。甲基化的主要形式有5-甲基胞嘧啶,N6-甲基腺嘌呤和7-甲基鸟嘌呤。原核生物中CCA/TGG和GATC常被甲基化，而真核生物中甲基化仅发生于胞嘧啶。DNA 的甲基化是在DNA 甲基化转移酶（DNMTs）的作用下使CpG二核苷酸5'端的胞嘧啶转变为5'甲基胞嘧啶。这种DNA 修饰方式并没有改变基因序列，但是它调控了基因的表达。脊椎动物基因的甲基化状态有三种：持续的低甲基化状态,如管家基因；去甲基化状态,如发育阶段中的一些基因；高度甲基化状态,如女性的一条失活的X染色体。

哺乳动物中，CpG序列在基因组中出现的频率仅有1%，远低于基因组中的其它双核苷酸序列。但在基因组的某些区域中，CpG序列密度很高，可以达均值的5倍以上，成为鸟嘌呤和胞嘧啶的富集区，形成所谓的CpG岛。通常，CpG岛大约含有500多个碱基。在哺乳动物基因组中约有4万个CpG岛，而且只有CpG岛的胞嘧啶能够被甲基化，CpG岛通常位于基因的启动子区或是第一个外显子区。健康人中，CpG岛中的CpG位点通常是处于非甲基化状态，而在CpG岛外的CpG位点则通常是甲基化的。这种甲基化的形式在细胞分裂的过程中能够稳定的保留。当肿瘤发生时，抑癌基因CpG岛以外的CpG序列非甲基化程度增加，而CpG岛中的CpG则呈高度甲基化状态，以致于染色体螺旋程度增加及抑癌基因表达的丢失。

随着高通量测序技术（NGS）技术的发展，使我们能够从全基因组水平来分析5’甲基胞嘧啶及组蛋白修饰等事件，由此能够发现很多传统的基因组学研究所不能发现的东西，这就是所谓的“DNA甲基化测序”！

DNA甲基化测序方法按原理可以分成三大类：

一、重亚硫酸盐测序；

二、基于限制性内切酶的测序；

三、靶向富集甲基化位点测序；

基于以上原理又有数种不同的测序方法，下面，就介绍10种DNA甲基化测序的常用方法及参考文献：

1) 重亚硫酸盐测序

该方法可以从单个碱基水平分析基因组中甲基化的胞嘧啶。首先，利用重烟硫酸盐对基因组DNA进行处理，将未发生甲基化的胞嘧啶脱氨基变成尿嘧啶。而发生了甲基化的胞嘧啶未发生脱氨基，因而，可以基于此将经重亚硫酸盐处理的和未处理的测序样本进行比较来发现甲基化的位点。

【相关文献】

Shotgun bisulphite sequencing of theArabidopsis genome reveals DNA methylation patterning

Highly integrated single-base resolutionmaps of the epigenome in Arabidopsis

2）重亚硫酸盐处理后接头标记技术(PBAT)

为了避免重亚硫酸盐处理时模板的丢失，通常会在重亚硫酸盐处理后进行接头连接和随机引物的扩增。

【相关文献】

Amplification-free whole-genome bisulfitesequencing by post-bisulfite adaptor tagging

3）限制性内切酶-重亚硫酸盐靶向测序(RRBS)

该技术是指对基因组上CpG岛或CpG甲基化较密集的区域进行靶向测序。样本首先经几种限制酶进行消化处理，然后经重亚硫酸盐处理，最后再测序。这种方法可以发现单个核苷酸水平的甲基化。

【相关文献】

Reced representation bisulfite sequencingfor comparative high-resolution DNA methylation analysis

4）氧化-重亚硫酸盐测序(oxBS-Seq)

5’羟甲基胞嘧啶(5’hmC)是5’甲基胞嘧啶脱甲基成胞嘧啶过程的中间产物，重亚硫酸盐测序无法对二者进行区分。通过氧化-重亚硫酸盐测序，5’甲基胞嘧啶保留，而5’羟甲基胞嘧啶(5’hmC)被氧化，进而脱氨基成尿嘧啶。通过将经过氧化处理和未处理的样本进行测序比较，即可从单个碱基水平分辨5’羟甲基胞嘧啶(5’hmC)和5’甲基胞嘧啶。

【相关文献】

Quantitative sequencing of 5-formylcytosinein DNA at single-base resolution.

5）TET辅助的重亚硫酸盐测序(TAB-seq)

TAB-seq采用葡萄糖亚胺与5’羟甲基胞嘧啶(5’hmC)作用来保护免受TET蛋白的氧化。5’甲基胞嘧啶和未甲基化的胞嘧啶被脱氨基成尿嘧啶，进而可以从单个碱基水平鉴定5’羟甲基胞嘧啶(5’hmC)。

【相关文献】

Base-resolution analysis of5-hydroxymethylcytosine in the Mammalian genome

6）甲基化敏感性的限制酶测序(MRE-Seq)

MRE-Seq将甲基化作用的敏感性和限制酶的特异性结合起来进而鉴定CpG岛的甲基化状态。

【相关文献】

Genome-scale DNA methylation analysis

7）HELP-Seq

HELP-Seq采用HpaII及其甲基化不敏感的限制性内切酶MSPI处理，来对基因组内及基因组间的甲基化位点进行比较，进而实现甲基化测序。

【相关文献】

Comparative isoschizomer profiling ofcytosine methylation: the HELP assay

8）甲基化DNA免疫共沉淀测序(MeDIP)

MeDIP是一种采用抗体或甲基化DNA结合蛋白来捕获富集甲基化DNA的技术，这种技术可以发现基因组中高度甲基化的区域，如CpG岛，但不能进行单个碱基水平的分析。

【相关文献】

Chromosome-wide and promoter-specificanalyses identify sites of differential DNA methylation in normal andtransformed human cells

9)甲基化结合域捕获技术(MBD-CAP)

MBD-CAP技术利用甲基化DNA能够结合蛋白MeCP2,MBD1-2 和 MBD3LI来对甲基化的DNA进行免疫沉淀。与MeDIP技术相似，该技术也是可以发现基因组中高度甲基化的区域，不能从单个碱基水平分析甲基化。

【相关文献】

High-resolution mapping of DNAhypermethylation and hypomethylation in lung cancer

10）基于探针的靶向富集技术

甲基化测序靶向富集技术采用合成寡核苷酸探针来捕获CpG岛、基因启动子区域以及其他一些显着性甲基化的区域。目前，Agilent 和 Roche Nimblegen公司已有这种商品化的试剂盒。

最后，Pacific Biosciences（Pacbio）公司的这项SMRT DNA测序技术采用动力学原理来直接检测甲基化的胞嘧啶。

其中三代测序技术针对真核中的5mc检测需要的深度较深，大概要250x,但是利用特殊的试剂会降低其深度。

㈣ DNA甲基化检测有哪些方法

DNA甲基化检测: Southern Blot杂交、 Bisulfite Sequencing and McrBC-PCR

㈤ DNA甲基化检测方法

为什么要检测 DNA 甲基化？ 首先，我们为什么要检测 DNA 甲基化呢？检测 DNA 甲基化能回答什么问题？我们都知道，DNA 甲基化可以调控基因表达，因此检测 DNA 甲基化至少可以为解释基因表达的调控提供一些线索。其次，DNA 甲基化参与一系列重要生物学过程，包括早期胚胎发育、基因组印记、X 染色体失活、重复序列的沉默以及癌症的发生发展转移，DNA 甲基化也有助于阐明这些重要的生命过程背后隐藏的奥秘。此外，DNA 甲基化一个非常重要的应用，是可以作为肿瘤的生物标志物。DNA 甲基化如此重要，我们可能需要时刻想着它。

图 2. 人类基因组中 DNA 甲基化概况

位点特异性的 DNA 甲基化 刚才我们介绍了总体 DNA 甲基化水平的检测，但是这些方法都没有提供序列信息。如果我们需要知道位点特异性的甲基化信息，这个时候需要问第二个问题，我是只检测感兴趣的几个基因就行了，还是需要知道全基因组每一个基因的甲基化状态。如果我们只需要检测特定基因的 DNA 甲基化状态，紧接着，我们需要问第三个问题，我们只是定性地看看甲基化状态，还是定量地看甲基化的水平。

图 13. 四个问题 2.0 版本之第二、三个问题 Methylation-specific PCR，MSP 如果只是定性地看甲基化的状态，最简单快捷的是甲基化特异性的 PCR。MSP 可以快速检测 CpG 岛中任意一组 CpG 位点的 DNA 甲基化状态。

图 16. 重亚硫酸盐转化&克隆测序[7]

基于高通量测序的全基因组 DNA 甲基化检测方法 接下来，我们看一下第四个问题的另一种答案，我们用高通量测序来检测全基因组的 DNA 甲基化。之所以最后讲这个，是因为这个比较烧钱，而且后续的生物信息学分析比较复杂。 WGBS & RRBS 对于二代测序的方法，WGBS（Whole-Genome Bisulfite Sequencing）是目前的“金标准”。WGBS 想必大家都已经很熟悉了，但是考虑到 WGBS 成本较高，因此还有一种跟 WGBS 十分相似，但是只测基因组中的一部分的方法，即 RRBS（Reced representation bisulfite sequencing）。

图 21. WGBS & RRBS 的比较[10]两者的差别在于，WGBS 通常使用超声打断整个基因组，并对整个基因组进行重亚硫酸盐处理和测序；但是 RRBS 使用 MspI 限制性内切酶。RRBS 可以检测到人类基因组中 85% 的 CpG 岛。RRBS 大约需要 100ng-1ug 的 DNA，而对于 WGBS 而言，需要 50-100ng DNA。尽管 WGBS 被誉为检测 DNA 甲基化的“金标准”，但是依然存在很多的不足。首先，重亚硫酸盐必须转化未甲基化的 DNA，否则会造成假阳性，不过目前我们可以通过加入 spike-in 作为对照来克服这个问题；还有，目前测到的 DNA 甲基化是多个细胞的平均甲基化状态，单细胞的 DNA 甲基化检测可能克服这个问题；其次，重亚硫酸盐转化，会降低基因组的复杂度，造成后续的 mapping 率不高；重亚硫酸盐转化还可能造成 DNA 断裂，这使得扩增长片段比较困难。另外，WGBS 的成本还是非常高的。

参考资料： https://mp.weixin.qq.com/s?__biz=MzIxMjA0NzU0OA==&mid=2650983484&idx=1&sn=&scene=21#wechat_redirect

㈥ 甲基化的甲基化检测方法

DNA甲基化是最早发现的基因表观修饰方式之一，真核生物中的甲基化仅发生于胞嘧啶，即在DNA甲基化转移酶（DNMTs）的作用下使CpG二核苷酸5’-端的胞嘧啶转变为5’-甲基胞嘧啶。DNA甲基化通常抑制基因表达，去甲基化则诱导了基因的重新活化和表达。这种DNA修饰方式在不改变基因序列前提下实现对基因表达的调控。脊椎动物DNA的甲基化状态与生长发育调控密切相关，比如在肿瘤发生时，抑癌基因CpG岛以外的CpG序列非甲基化程度增加，CpG岛中的CpG则呈高度甲基化状态，导致抑癌基因表达的下降。
1、甲基化特异性的PCR（Methylation-specific PCR，MSP）
用亚硫酸氢盐处理基因组DNA，所有未发生甲基化的胞嘧啶被转化为尿嘧啶，而甲基化的胞嘧啶不变；随后设计针对甲基化和非甲基化序列的引物进行PCR。通过电泳检测MSP扩增产物，如果用针对处理后甲基化DNA链的引物能得到扩增片段，则说明该位点存在甲基化；反之，说明被检测的位点不存在甲基化。
2、亚硫酸氢盐测序法（Bisulfite sequencing PCR，BSP）
用亚硫酸氢盐处理基因组DNA，则未发生甲基化的胞嘧啶被转化为尿嘧啶，而甲基化的胞嘧啶不变。随后设计BSP引物进行PCR，在扩增过程中尿嘧啶全部转化为胸腺嘧啶，最后对PCR产物进行测序就可以判断CpG位点是否发生甲基化称为BSP-直接测序方法。将PCR产物克隆至载体后进行测序，可以提高测序成功率，这种方法称为BSP-克隆测序法。
3、高分辨率熔解曲线法（High Resolution Melting，HRM）
在非CpG岛位置设计一对针对亚硫酸氢盐修饰后的DNA双链的引物，这对引物中间的片段包含感兴趣的CpG岛。若这些CpG岛发生了甲基化，用亚硫酸氢盐处理后，未甲基化的胞嘧啶经PCR扩增后转变成胸腺嘧啶，而甲基化的胞嘧啶不变，样品中的GC含量发生改变，从而导致熔解温度的变化（图1）。
其中，样品要求：细胞（≥106 个）、组织（≥300mg）、血液（≥1ml）、血清（≥1.5ml）等样品材料，基因组DNA（体积≥20μl，浓度≥50 ng/μl）。

㈦ 易基因｜全基因组DNA甲基化测序分析全流程

全基因组DNA甲基化实验怎么做？从技术原理、建库测序流程、信息分析流程和研究套路等四方面详细介绍。

表观修饰不需要改变 DNA 序列便能实现对性状的改变，表观修饰的改变与基因功能乃至细胞状态段烂陆、发育、衰老、疾病等存在重要的关联。在众多的表观遗传修饰中，最为重要且研究最为广泛的修饰之一是 DNA 甲基化，而全基因组甲基化测序（WGBS-seq）无疑是最有效的研究手段。

全基因组甲基化测序利用重亚硫酸盐能够将未甲基化的胞嘧啶（C）转化为胸腺嘧啶 （T）的特性，将基因组用重亚硫酸盐处理后测序，即可根据单个 C 位点上未转化为 C 未转化为 T 的 reads 数目与所有覆盖的 reads 数目的比例，计算得到甲基化率。该技术对于全面研究胚胎发育、衰老机制、疾病发生发展的表观遗传机制，以及筛选疾病相关的表观遗传学标记位点具有重要的应用价值。

全基因组甲基化测序原理示意图入下：

样品检测——样品打断 ——文库构建——BS处理——文库质检

（一）样品检测

对DNA样品的检测主要包括2种方法：

（1）琼脂糖凝胶电泳分析DNA降解程度以及是否有污染，检测具有明显的主带，且条带清晰；

Qubit 2.0对DNA浓度进行精确定量，DNA检测总量不低于1ug。

（二）文库构建

样本检测合格后，使用Bioruptor系统将1µg样品基因组DNA与未甲基化的lambda DNA混合，然后将其片段化，平均大小约为250bp。片段化后，纯化的随机片段化DNA随后用T4 DNA聚合酶，Klenow片段和T4多核苷酸激酶的混合物进行修复，钝化和磷酸化末端。随后使用Klenow片段（3'-5'exo-）对钝的DNA片段进行3'腺苷酸化，然后与连接5'-甲基胞嘧啶而不是使用T4 DNA连接酶的胞嘧啶连接的衔接子进行连接。完成每个步骤后，使用磁珠纯化DNA。之后，根据说明使用ZYMO EZ DNA甲基化金试剂盒将未甲基化的胞嘧啶转化为尿嘧啶。最后，用JumpStart Taq DNA聚合酶进行PCR扩增，再使用磁珠对PCR产物进行纯化获得最终文库。

（三）文库质检

文库构建完成后，先使用Qubit2.0进行初步定量，稀释文库至1ng/ul，随后使用Agilent 2100对文库的insert size进行检测，insert size符合预期后，使用qPCR方握顷法对文库的有效浓度进行准确定量（文库有效浓度> 2nM），以保证文库质量。

（四）上机测序

文库检测合格后，把不同文库按照有效浓度及目标下机数据量的需求pooling后在illumina Nova平台测序，测序策略为PE150。

（一）原始下机数据质控

原始下机数据为FASTQ格式，是高通量测序的标准格式。FASTQ文件每四行为一个单位，包含一条测序序列（read）的信息。该单位第一行为read的历逗ID，一般以@符号开头；第二行为测序的序列，也就是reads的序列；第三行一般是一个+号，或者与第一行的信息相同；第四行是碱基质量值，是对第二行序列的碱基的准确性的描述，一个碱基会对应一个碱基质量值，所以这一行和第二行的长度相同。以下为一条read信息的示例：

原始下机数据包含建库时引进的接头序列以及质量过低的碱基，这些因素会导致后续比对到基因组的reads较少，从而导致得到的信息较少，因此需要进行过滤。利用trim_galore软件对原始数据进行去除接头序列及低质量碱基等质控步骤。

（二）序列比对

经过质控的reads需要根据与参考基因组的序列相似度比对到参考基因组上。相比于常规基因组及转录组测序，WGBS测序方法产生的数据的特点决定其在比对时存在三大困难：

（1）DNA片段正链和负链经过重亚硫酸盐转化后将不再反向互补，再经过PCR，便会产生四条不同的序列，这将大大增加比对时的计算量。

（2）经过重亚硫酸盐转化后，DNA序列大部分C碱基被转化成T碱基，因此序列含大量T而缺乏C；经过PCR后，产生的互补链则含有大量A而缺乏G。这样便导致序列的复杂度降低（即序列的组成特征更单一），从而增加比对的难度。

（3）C和T的比对是不对称的。经过重亚硫酸盐转化后，序列中非甲基化的C碱基（占大部分）被转化为T，这将导致测序序列与参考基因组不匹配，T既可能应该比对到T上，有可能应该比对到C上；而C则只能比对到C上。这也增加了比对的难度。

利用BSMAP软件进行比对。BSMAP进行比对时，先以参考基因组上C碱基的位置作为指导，将reads中对应参考基因组C碱基位置的T标记为C，其他T保持不变，从而使reads可以直接比对到参考基因组。

（三）甲基化水平计算

甲基化水平可根据未转化为 T 的 C 与转化为 T 的 C 的 reads 的比例计算得到，即：

Beta-value = C-reads / (C-reads + T-reads) * 100%

其中，Beta-value 即为该胞嘧啶的甲基化水平，C-reads 为覆盖该位点的支持甲基化的reads 数目（测得该位点为 C 的 reads），T-reads 为覆盖该位点的不支持甲基化的 reads 数目（测得该位点为 T 的 reads）。 计算原理示意图如下：

利用BSMAP统计甲基化水平。

（四）差异甲基化区域（DMR）鉴定及统计

DMR检测使用权威期刊发表的metilene软件。该软件先将基因组进行预分段，以排除较长序列中不包含CG位点的片段。随后，利用二元分隔算法，递归缩小检测范围，以搜索得到组间累积平均甲基化差异最大的区域，作为可能的DMR；最后，结合双重统计学检验（MWU-test和2D KS-test），得到准确的DMR。检测原理如下图所示：

本分析检测DMR的标准如下：

（1）区域平均甲基化差异不小于0.1；

（2）CpG位点数不少于5个；

（3）区域长度不小于50 bp；

（4）甲基化水平差异统计检验的校正P值小于0.05；

（5）2D KS-test检验P值小于0.05。

（五）信息分析流程示意图

DNA甲基化组学研究的核心内容在于对DNA甲基化数据的挖掘。DNA甲基化一般遵循三个步骤进行数据挖掘。

首先，进行整体全基因组甲基化变化的分析，包括平均甲基化水平变化、甲基化水平分布变化、降维分析、聚类分析、相关性分析等。

其次，进行甲基化差异水平分析，筛选具体差异基因，包括DMC/DMR/DMG鉴定、DMC/DMR在基因组元件上的分布、DMC/DMR的TF结合分析、时序甲基化数据的分析策略、DMG的功能分析等。

最后，将甲基化组学&转录组学关联分析，包括Meta genes整体关联、DMG-DEG对应关联、网络关联等。

Whole-Genome Bisulfite Sequencing of Two Distinct Interconvertible DNA Methylomes of Mouse Embryonic Stem Cells. 两种状态的小鼠胚胎干细胞的甲基化组学研究

1、背景

小鼠胚胎干细胞一般生长在含有血清的基质中，被称作血清干细胞(serum ESCs)；加两种激酶抑制因子使胚胎干细胞在无血清的情况下更能保持多能性的基态，这种干细胞称为2i干细胞(2i ESCs)；这两种状态的胚胎干细胞可以互相转化。以前这方面的甲基化研究大多基于质谱，覆盖度和研究结果有限，尚缺乏2i胚胎干细胞的甲基化组学研究。

2、方法

利用全基因组重亚硫酸盐甲基化测序（WGBS），对这两种可互相转换的小鼠胚胎干细胞进行甲基化组学研究

3、结论

全面准确的检测了两种小鼠胚胎干细胞的DNA甲基化修饰并进行了系统的比较；同serum ESCs相比，雄性2iESCs全局低甲基化；在血清中，雌性ESCs跟雄性2i ESCs类似呈现全局低甲基化，而在2i ESCs状态下，甲基化水平会进一步降低。

以上就是关于全基因组甲基化测序实验流程和分析思路的介绍。

参考文献：

[1] Ashburner, M. and C. A. Ball, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet, 2000, 25 (1): 25-9.

[2] Dirk Schübeler. Function and information content of DNA methylation. Nature, 2015, 517: 321–326.

[3] Frank Jühling et al. metilene: Fast and sensitive calling of differentially methylated regions from bisulfite sequencing data. Genome Research, 2016, 26: 256-262.

[4] Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic acids research, 2000,28(1): 27-30.

[5] Tadafumi Kato Kazuya Iwamoto. Comprehensive DNA methylation and hydroxymethylation analysis in the human brain and its implication in mental disorders. Neuropharmacology, 2014, 80: 133-139.

[6] Xiaojing Yang et al. Gene Body Methylation Can Alter Gene Expression and Is a Therapeutic Target in Cancer. Cancer Cell 26, 577–590.

[7] Yuanxin Xi et al. BSMAP: whole genome bisulfite sequence MAPping program. BMC Bioinformatics, 2009, 10:232.

[8] Gao F, et al. De novo DNA methylation ring monkey pre-implantation embryogenesis. Cell Res. 2017 Apr;27(4):526-539. pii: cr201725.

㈧ 如何检测甲基化的程度

先用中亚硫酸盐处理,使DNA中未发生甲基化得胞嘧啶脱氨基转化裤卖段为尿配正嘧啶而甲基化的胞嘧啶保持不变,PCR扩增得到所需片段,片段中尿嘧啶全部转化成T,然胡誉后对PCR产物测序,并与未处理的序列比较,判断CpG位点是否发生甲基化.我们研究院对甲基化检测很专业,感兴趣的话,可以看看：中关村华康基因研究院

㈨ 谁知道测DNA甲基化的方法

http://www.ebiotrade.com/emgzf/genenews200404/gene5.htm您自己看吧！好多图片粘贴不过来，还不如你自己看。

一种全新的DNA甲基化研究方法——Pyrosequencing技术

DNA甲基化是一种表观遗传修饰，它是由DNA甲基转移酶（DNA methyl-transferase, Dnmt）催化S-腺苷甲硫氨酸作为甲基供体，将胞嘧啶转变为5-甲基胞嘧啶（mC）的一种反应，在真核生物DNA中，5-甲基胞嘧啶是唯一存在的化学性修饰碱基。CG二核苷酸是最主要的甲基化位点，它在基因组中呈不均匀分布，存在高甲基化、低甲基化和非甲基化的区域，在哺乳动物中mC约占C总量的2－7％。一般说来，DNA的甲基化会抑制基因的表达。DNA的甲基化对维持染色体的结构、X染色体的失活、基因印记和肿瘤的发生发展都起重要的作用。

CpG双核苷酸在人类基因组中的分布很不均一，而在基因组的某些区段，CpG保持或高于正常概率，这些区段被称作CpG岛。CpG岛主要位于基因的启动子和第一外显子区域，约有60％以上基因的启动子含有CpG岛。 CpG甲基化的研究在肿瘤的研究中有着非常主要的地位。通过基因启动子区及附近区域CpG岛胞嘧啶的甲基化可以在转录水平调节基因的表达，从而引起相应基因沉默，去甲基化又可恢复其表达。DNA甲基化在生理情况下就参与了控制基因的时空表达，在肿瘤发生时，肿瘤细胞全基因组低甲基化是一个重要特征。肿瘤细胞基因组甲基化的程度与正常细胞相比仅为20-60% , 同时伴有局部区域基因的高甲基化,包括肿瘤抑制基因、抑制肿瘤转移和浸润的基因、细胞周期调节基因、DNA修复基因、血管形成抑制基因等。但是目前研究手段的局限，限制了DNA甲基化的广泛研究。

近年来，研究者不断探索定性及定量检测单个或多个甲基化位点的方法，但由于甲基化多态性区域存在的密度很高，所以对于延伸反应其引物的位置很难设计。Pyrosequencing技术作为一种新的序列分析技术，能够快速地检测甲基化的频率，对样品中的甲基化位点进行定性及定量检测，为甲基化研究提供了新的途径。

从原理上来看，Pyrosequencing是一种通过合成方法进行序列分析的方法，它通过核苷酸和模板结合后释放的焦磷酸引发酶级联反应，促使荧光素发光并进行检测。这项技术曾经被用作单核苷酸多态性（SNP）的基因型和单倍型的检测，以及细菌和病毒的鉴定和分型研究。这项技术的一个主要特点是在Pyrogarm™软件上显示的峰值高度来自于序列分析的原始数据，通过峰值的高度可以精确的检测混合DNA模板中等位基因的频率。

目前甲基化研究方面，很多甲基化定量分析的报道采用亚硫酸氢盐处理甲基化样本，并用混合的PCR产物
作为校正。其主要原理是：亚硫酸氢盐可以将没有甲基化的胞嘧啶转化为尿嘧啶，而在适当的实验条件下甲基化的胞嘧啶保持不变。因而，用它处理样本后，再进行PCR扩增，甲基化的位点可以被当作一个C/T的SNP来处理，它的基因频率为0－100％。在此，我们给大家介绍一个研究人员使用Pyrosequencing技术分析并精确定量DNA甲基化水平的例子。

研究者在一个Pyrosequencing反应中同时检测了6个甲基化位点。这种方法同样可以用于石蜡包埋的组织，并且具有较高的重复性和精确性。实验选择谷胱甘肽-S-转移酶π（GSTP1）转录启动位点的CpG岛进行检测。这些位点在正常前列腺组织中是非甲基化的，而在肿瘤样本中高甲基化。通过PCR扩增一个包含17个甲基化多态位点140bp的片段，并用4个测序引物研究其中15个位点（Table 1）。使用在线的SNP测序引物设计软件（Pyrosequencing AB）设计测序引物，其中一些甲基化多态性位点用最可能的碱基所代替，以减少计算的数量。再通过人工检测测序引物可能存在的错配。此外，同时在PSQ 96MA DNA分析仪上运行空白对照，扣除由测序引物、生物素标记的引物或是模板引起的背景。

PCR引物设计完全与模板相匹配，不覆盖任何甲基化多态性区域。使用10ng亚硫酸氢盐转化的DNA样本或是10 fmol纯化的PCR产物，10 pmol 正向（5’-GTGATTTAGTATTGG-3’）和反向（5’-biotin-AACTCTAAACCCCATC-3’）引物扩增GSTP1转录启动位点的基因片段，扩增片段长度为144bp。反应体系为60 mM Tris-SO4, pH 8.9, 18 mM (NH4)2SO4, 1 mM MgSO4, 200 μM dNTPs，以及3 U Platinum Taq DNA高保真聚合酶，终体积为50μL。PCR循环设置：首先在95℃下变性4分钟，然后在95℃ 30S，50℃ 45S以及72℃ 20S条件下重复50个循环，最后一步延伸步骤在72℃下4分钟，中止反应。PCR反应在Eppendorf的Mastercycler 96

哺乳动物基因组中，DNA甲基化是指CpG二核苷酸中的胞嘧啶第5位碳原子被甲基化. DNA甲基化是一种基因外修饰，不改变DNA的一级结构; 他在细胞正常发育、基因表达模式以及基因组稳定性中起着至关重要的作用. 全基因组低甲基化，维持甲基化模式酶的调节失控和正常非甲基化CpG岛的高甲基化是人类肿瘤中普遍存在的现象. DNA高甲基化是导致抑癌基因失活的又一个机制.
http://www.wjgnet.com/1009-3079/abstract_cn.asp?f=1420&v=11

㈩ 甲基化测序

         WGBS全称Whole Genome Bisulfite Seuqneicng： 即全基因组重亚硫酸盐测序。该方法通过Bisulfite处理，将原基因组中未发生甲基化的C碱基转换成U的同时，保留所有甲基化C的碱基不发生转变，从而帮助科研人员识别发生甲基化的CpG位点。该种测序技术适用于绘制单碱基分辨率的全基因组DNA甲基化图谱。

        誉凯仿 RRBS全称Reced Representation Bisulfite Sequencing： 即简化代表性重亚硫酸盐测序孙雀。该方法在Bisulfite处理前，使用MspI（该酶的酶切位点为CCGG）酶切对样本进行处理， 去除低CG含量DNA片段 ，从而使用较小的数据量富集到尽可能多的包含CpG位点的DNA片段。

        相比于WGBS技术，RRBS是一种准确、高效且经济的DNA甲基化研究方法，庆纤通过酶切，并进行Bisulfite测序，该方法在保证DNA甲基化状态检测的高分辨率的同时提升测序数据的高利用率。

甲基化数据处理流程：