dna甲基化测序文库及其构建方法和检测方法
技术领域
1.本技术涉及基因检测技术领域，具体涉及一种dna甲基化测序文库及其构建方法和检测方法。


背景技术：

2.生物体内的dna包括四种碱基，分别为腺嘌呤(a)、鸟嘌呤(g)、胸腺嘧啶(t)和胞嘧啶(c)。其中，部分胞嘧啶(c)会发生甲基化而带上甲基基团。dna甲基化是dna甲基化转移酶将甲基基团选择性地添加到胞嘧啶上以形成甲基化胞嘧啶的过程。dna甲基化的水平通常与诱导或抑制基因的表达有关。如果dna的甲基化发生异常，则可能与肿瘤的发生、发展有关，因此，通过特定基因进行测序，以获知其甲基化水平，这对于肿瘤早筛有重要的价值。
3.目前通过高通量测序平台来对特定基因进行测序，这首先需要根据特定基因来建立dna甲基化测序文库。现有的dna甲基化测序文库的建库方法包括：基于亚硫酸盐转化法的单链建库法和基于酶学转化法的双链建库法。
4.如图1所示，基于亚硫酸盐转化法的单链建库法包括如下步骤：将目标双链dna变性为单链dna；通过化学反应改变单链dna中胞嘧啶上的甲基，以产生受保护基团保护的胞嘧啶；利用转化试剂亚硫酸盐将未甲基化的胞嘧啶转化为尿嘧啶，然后建立用于供测序平台测序的dna甲基化文库。该方法是目前广泛采用的文库构建方法，性能比较稳定，但是亚硫酸盐处理对dna的损伤较大，特别是针对末端缺失的双链dna，会造成dna断链、碱基丢失(例如形成无碱基位点)、甲基化信息失真等现象。
5.如图2所示，基于酶学转化法的双链建库法包括如下步骤：将目标双链dna进行3’末端修复，并直接加接头序列，接头序列里的胞嘧啶(c)为全部甲基化(如图2中m所示)的胞嘧啶，在双链状态下通过化学反应改变单链dna中胞嘧啶上的甲基，以产生受保护基团保护的胞嘧啶；然后变性为单链dna，采用转化酶将未甲基化的胞嘧啶转化为尿嘧啶(u)，然后通过pcr扩增来建立用于供测序平台测序的dna甲基化文库。因为接头序列在最开始已经加上了，因此，转化后直接通过pcr扩增来构建文库。该方法没有采用易对dna产生损伤的亚硫酸盐，因此，该方法对dna结构的损伤小，不会造成dna断链、碱基丢失等现象，但是该方法采用未甲基化的dntp对双链dna的缺失末端进行修复。针对双链不平整，带有粘性单链末端的dna而言，修复的碱基也会被测序仪读取，因此该方法在对缺失末端进行修复时引入新的、不同的甲基化水平，这会导致甲基化信息失真。另外，由于难以了解原始的双链dna到底有多少个碱基缺失，这些失真的甲基化信息无法在后续生物信息学分析中修正。再者，测序出来的dna序列中的碱基数目多于原始的dna序列的碱基数目，这会造成数据污染。


技术实现要素：

6.本技术提供一种dna甲基化测序文库及其构建方法和检测方法，其采用的试剂性质温和，对dna结构的损伤小，不会造成dna断链、碱基丢失等现象，也不会产生甲基化信息失真和数据污染等现象。
7.为了克服上述缺点，本技术采用了以下技术方案：
8.[dna甲基化测序文库的构建方法]
[0009]
本技术提供一种dna甲基化测序文库的构建方法，其包括如下步骤：
[0010]
(1)、获得目标双链dna；
[0011]
(2)、采用甲基化保护试剂将目标双链dna的甲基化胞嘧啶氧化为受保护胞嘧啶，以得到受保护双链dna；
[0012]
(3)、将受保护双链dna变性为单链dna；
[0013]
(4)、采用转化试剂将单链dna中未甲基化的胞嘧啶转化为尿嘧啶，得到转化后单链dna；
[0014]
(5)、利用转化后单链dna构建dna甲基化测序文库。
[0015]
其中，在步骤(1)中，目标双链dna可以为血液中游离的cfdna。
[0016]
在一些实施例中，cfdna在血液中可以主要以缠绕核小体的形式存在，因为在该种情况下没有缠绕在核小体上的cfdna会很快被降解。一个核小体上缠绕的cfdna大约170bp，因此，cfdna的主要片段在170bp存在一个主峰(相当于一个核小体)，然后在340bp有一个小峰(相当于两个核小体)，以此类推。
[0017]
在另一些实施例中，cfdna的长度可以位于170bp至210bp之间(包括本数)。然而，在其它一些实施例中，cfdna的长度也可以为180bp、185bp、190bp、195bp、200bp等。
[0018]
在其它的一些实施例中，cfdna的长度可以为170bp至210bp之间的任意一个数值的整数倍(根据所缠绕的核小体的个数确定)。例如为340bp、420bp等。
[0019]
在步骤(1)中，目标双链dna也可以为ctdna。ctdna是cfdna中来自于肿瘤细胞的部分，其来自于组织细胞的非正常凋亡过程，因此，所产生的片段大小不同，其片段的均一性小于组织细胞的正常凋亡过程形成的cfdna。
[0020]
在一些实施例中，ctdna的长度位于150bp至500bp之间(包括本数)。例如，ctdna的长度也可以为210bp、220bp、230bp、250bp、270bp、280bp、300bp、320bp、350bp、380bp、400bp、500bp等。
[0021]
在步骤(1)中，目标双链dna也可以为基因组dna。在构建文库时，该基因组dna需要使用打断仪打断为100bp至400bp之间(包括本数)。在一些实施例中，其端点值可以参照上述的cfdna或ctdna取值。
[0022]
在步骤(1)中，目标双链dna的每条链的3’端可以存在1至80个碱基的缺失，也可以为5至70个碱基的缺失，也可以为10至50个碱基的缺失，可以进一步为20至30个碱基的缺失。在其它一些实施例中，目标双链dna的每条链的3’端的碱基缺失数可以为2、3、5、6、7、8、9、10、11、12、13、14、15、16、17、18、19、25、35、45、55、65、75等。以上情况除了适用于3’末端缺失的目标双链dna之外，也适用于5’端缺失的目标双链dna。
[0023]
在步骤(2)中，本技术的目标双链dna的起始量可以为3ng至200ng，也可以为5ng至100ng，还可以为20ng至50ng。起始量为建立测序文库时最初使用的目标双链dna的用量。
[0024]
在步骤(2)中，甲基化保护试剂包括tet2蛋白和氧化增强剂。tet2蛋白用于将甲基化胞嘧啶中的甲基和羟甲基经过氧化反应转换为保护基团(例如，甲酰基或羧基)。受到保护基团保护的胞嘧啶能够对抗胞嘧啶脱氨酶的脱氨作用，从而不被转换为尿嘧啶。未甲基化的胞嘧啶不存在甲基，不会连接上保护基团，从而不能够对抗胞嘧啶脱氨酶的脱氨作用，
后继会被转化为尿嘧啶。
[0025]
在步骤(2)中，甲基化胞嘧啶为5-甲基胞嘧啶或5-羟甲基胞嘧啶；
[0026]
在步骤(2)中，受保护胞嘧啶为5-甲酰基胞嘧啶或5-羧基胞嘧啶。
[0027]
在步骤(3)中，变性的温度为85℃，添加的变性剂为甲酰胺。
[0028]
在步骤(4)中，转化试剂包括apobec蛋白。apobec蛋白具有胞嘧啶脱氨酶的活性，能够使得未甲基化的胞嘧啶脱氨，从而转化为尿嘧啶。
[0029]
上述的步骤(5)具体包括如下步骤：
[0030]
(5-1)、采用衔接酶在转化后单链dna的3’端衔接上部分双链的链截短接头序列，得到单链接头产物；
[0031]
(5-2)、采用dna聚合酶使单链接头产物中转化后单链dna利用底物dntp分子以该截短接头序列为聚合起点沿5’端至3’端的方向聚合成双链dna，得到双链延伸产物；
[0032]
(5-3)、在双链延伸产物的未添加截短接头序列的另一端连接上双链接头序列，得到双链接头产物；
[0033]
(5-4)、采用dna连接酶在双链接头产物中双链dna的5’端分别添加具有双链结构的测序接头序列，进行标记pcr反应(indexing pcr)，得到dna甲基化测序文库。
[0034]
在步骤(5-4)中，标记pcr反应所使用的标记引物包括：编号为cat.no.36024的index x、编号为cat.no.38096的index d50x、编号为cat.no.38096的index d7xx、以及编号为cat.no.39096和390384的index u001-u096。
[0035]
[dna甲基化测序文库]
[0036]
本技术提供了一种dna甲基化测序文库，该dna甲基化测序文库由上述的构建方法构建而成。该dna甲基化测序文库的每个dna分子的两个平头端带有测序接头序列，能够被illumina测序仪读取，从而能够测得该dna的甲基化水平。
[0037]
[dna甲基化水平的检测方法]
[0038]
本技术提供了一种dna甲基化水平的检测方法，其包括如下步骤：
[0039]
(1)、对目标双链dna按照上述的构建方法构建dna甲基化测序文库；
[0040]
(2)、采用测序仪对dna甲基化测序文库进行测序，得到测序结果；
[0041]
(3)、对测序结果进行分析，得到dna甲基化水平。
[0042]
由于采用了上述技术方案，本技术取得了如下的技术效果：
[0043]
首先，本技术先将末端缺失或未缺失的双链dna变性为单链dna，然后利用该单链dna建立dna甲基化测序文库，因此，本技术并未对该双链dna进行末端修复(如5’末端修复和/或3’末端修复)，能够避免由于末端修复而带来的修复末端的碱基被测序仪读取的现象，从而避免由此造成的测序数据冗余现象。
[0044]
另外，本技术并未采用未甲基化的dntps进行末端修复，不会在末端引入新的甲基化水平，不会引起甲基化水平失真现象。
[0045]
再者，本技术在将未甲基化的胞嘧啶转化为尿嘧啶时采用反应条件温和的酶学转化法，而并未添加亚硫酸盐，因此，不会造成dna断链、碱基丢失等现象。
附图说明
[0046]
为了更清楚地说明本技术实施例中的技术方案，下面将对实施例描述中所需要使
用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0047]
图1为基于亚硫酸盐转化的单链建库法的步骤示意图。
[0048]
图2为基于酶学转化的双链建库法的步骤示意图。
[0049]
图3为本技术的基于酶学转化的单链建库法的步骤示意图。
[0050]
图4为本技术与现有酶学转化甲基化建库的测序结果对比结果图。
[0051]
上述各图中的m表示甲基化胞嘧啶上的甲基。p
*
表示受保护胞嘧啶上的保护基团。上述各图只显示了双链dna中其中一条链的处理方法，另一条链作同样处理，故上述各图仅仅显示其中一条链。在测序文库构建的过程中，尿嘧啶u会转换为胸腺嘧啶t，故最终的测序文库不再显示尿嘧啶。上述各图中的黑色区域表示dna的链，空白区域和/或花纹区域表示接头序列。不同的花纹区域表示不同的接头序列。
具体实施方式
[0052]
下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
[0053]
以下分别对本技术进行详细说明。需说明的是，以下实施例的描述顺序不作为对实施例优选顺序的限定。
[0054]
除非另有说明，文中涉及的试剂与材料均可商购获得，或本领域技术人员可依据公知常识自行制备。
[0055]
如图3所示，本技术提供了一种dna甲基化测序文库的构建方法，其特征在于，包括如下步骤：
[0056]
(1)、获得目标双链dna；
[0057]
(2)、采用甲基化保护试剂(又称氧化试剂)将目标双链dna的甲基化胞嘧啶氧化为受保护胞嘧啶，即相当于在胞嘧啶上连接了保护基团，并且未甲基化的胞嘧啶没有连接上保护基团。由此，得到受保护双链dna；
[0058]
(3)、将受保护双链dna变性为单链dna；
[0059]
(4)、采用转化试剂(又称脱氨试剂)将单链dna中未甲基化的胞嘧啶转化为尿嘧啶，已连接上保护基团的胞嘧啶由于受到保护基团的保护而不会转化为尿嘧啶，从而得到转化后单链dna；
[0060]
(5)、利用转化后单链dna构建dna甲基化测序文库。
[0061]
其中，在步骤(1)中，目标双链dna的获取方法可以采用常用的dna提取方法。图3仅显示了3’末端缺失的双链dna，实际上，本技术的方法除了适用于3’末端缺失的双链dna之外，也适用于5’末端缺失的双链dna，另外也适用于具有平头末端的双链dna(即其不存在末端缺失)。
[0062]
在步骤(1)中，在一些实施例中，所获得的目标双链dna可以为末端缺失的双链dna。生物体的血液中往往会出现处于游离状态的末端缺失的双链dna，例如，由组织细胞经
由正常的凋亡过程而产生的循环游离dna(circulating free dna，cfdna)，或由肿瘤细胞经由非正常的凋亡过程而产生的循环肿瘤dna(circulating tumor dna，ctdna)等。
[0063]
正常的cfdna主要是通过细胞凋亡过程中产生的，长度呈均匀分布并且位于170至210bp之间，或者为170至210bp之间的任意一个自然数值的倍数。浓度一般小于或等于100ng/ml，平均浓度约为30ng/ml。
[0064]
ctdna则是在非正常凋亡过程中产生，其特异性地来自于肿瘤细胞，能够用作对肿瘤细胞的超早期筛查。该ctdna的长度呈大小不同分布并且位于150bp至500bp之间(包括本数)，浓度按肿瘤进展期不同一般为0.01至10ng/ml，平均浓度约为0.1ng/ml。这些双链dna在进入血液中后，由于受到一些酶的作用，其3’端往往有1-80个碱基的缺失，从而出现不平整的末端。如果在建立测序文库时首先采用填补末端的方法，由于末端修复一般会采用未甲基化的dntps，因此，末端修复会导致原始双链dna的末端碱基缺失数目无法精确测得，这会造成测序所得到的碱基数目失真，并且也会在末端引入新的甲基化水平而造成甲基化信息失真。
[0065]
本技术中，如果待测样本来自于基因组dna(genomic dna，gdna)，则需要将该待测样本中的dna打断，以得到片段化的dna，然后再构建dna甲基化测序文库。如果待测样本来自于cfdna，由于cfdna本身的片段较短，故无需打断过程。组织细胞提取的dna也会存在不平整的粘性末端，但是由于组织dna一般的片段的长度较长，实验中需要通过超声波机械打断成短片段，而超声打断的“切口”一般较为整齐，所以相对来说不平整末端的影响并不那么大。
[0066]
在步骤(2)中，甲基化保护试剂包括tet2蛋白和氧化增强剂。
[0067]
在步骤(2)中，甲基化胞嘧啶为5-甲基胞嘧啶或5-羟甲基胞嘧啶，受保护胞嘧啶为5-甲酰胞嘧啶或5-羧基胞嘧啶。
[0068]
在步骤(3)中，变性的温度为85℃，添加的变性剂为甲酰胺。
[0069]
在步骤(4)中，转化试剂包括apobec蛋白。
[0070]
步骤(5)具体包括如下步骤：
[0071]
(5-1)、采用衔接酶在转化后单链dna的3’端衔接上部分双链的截短接头序列，得到单链接头产物；其中，“单链”指的是转化后单链dna，“接头”指的是具有部分双链的截短接头序列。具有部分双链的截短接头序列为swift试剂盒的truncated adaptor，试剂盒名称为swift accel-ngs methyl-seq dna library kit。该截短接头序列的一部分为单链，一部分为双链，为了简便描述起见，在图3中以单链表示，但该图并不意味着该截短接头序列整体上为单链结构或双链结构。
[0072]
(5-2)、采用dna聚合酶使单链接头产物中转化后单链dna利用底物dntp分子以截短接头序列为聚合起点沿5’端至3’端的方向聚合成双链dna，得到双链延伸产物；
[0073]
(5-3)、通过连接反应在双链延伸产物的未添加截短接头序列的另一端连接上双链接头序列，得到双链接头产物；该双链接头序列具有双链结构，能同时连接到双链延伸产物的未添加截短接头序列的另一端的5’末端和3’末端；本步骤是单链建库流程中的一个步骤，以便后继在原始链5’端方向上连接index接头序列，如果没有这一步，则在原始链5’端方向无法进行后续的接头连接；
[0074]
(5-4)、采用dna连接酶在双链接头产物的双链dna的两个5’端分别添加具有双链
结构的测序接头序列，并进行标记pcr反应(indexing pcr)，得到dna甲基化测序文库。测序接头序列也具有双链结构，因此，能同时连接到双链接头产物的两个平头末端。该测序接头序列含有index序列(图3中未显示该index序列)。
[0075]
由于测序仪进行一次测序需要处理很多的样本，为了避免不同来源的样本之间的混淆，同一来源的每个样本都有一段特定的测序接头序列，每个测序接头序列里面有一段特定的序列，该段特定的序列称为index或者barcode。本技术的测序采用illumina公司的测序仪及其测序原理。因此，index的具体序列来自于swift bioscience公司的商用index试剂盒。
[0076]
在步骤(5-4)中，标记pcr反应所使用的标记引物包括：编号为cat.no.36024的index x、编号为cat.no.38096的index d50x、编号为cat.no.38096的index d7xx、以及编号为cat.no.39096和390384的index u001-u096。
[0077]
本技术还提供了一种dna甲基化测序文库，该dna甲基化测序文库由上述任意一种构建方法或其组合构建而成。
[0078]
本技术还提供了一种dna甲基化水平的检测方法，其包括如下步骤：
[0079]
(1)、对目标双链dna按照上述的任意一种构建方法或其组合构建dna甲基化测序文库；
[0080]
(2)、采用测序仪对dna甲基化测序文库进行测序，得到测序结果；
[0081]
(3)、对测序结果进行分析，得到dna甲基化水平。
[0082]
本技术的方法没有末端修复的过程，故既适用于末端缺失(5’末端缺失或者3’末端缺失)的双链dna，又适用于具有平头端的双链dna。
[0083]
在本技术中，目标双链dna先变性为单链dna，然后在单链dna的3’端衔接上截短接头序列，以截短接头序列为起点通过5’端至3’端的延伸反应补平单链dna，从而形成平末端dna，因此，片段化的目标双链dna的5’磷酸基团并不参与延伸反应，无需去除该5’磷酸基团，故本技术无需加入能够去除dna的5
’‑
磷酸基团并且3’端磷酸酶活性缺失的t4多聚核苷酸激酶。
[0084]
在本技术中，dna样品没有经过重亚硫酸盐处理，而是经过较温和的酶处理，dna的结构不会受损伤，使得测得的甲基化水平更加接近真实值。
[0085]
本技术的dna甲基化测序文库能够适用于illumina高通量测序平台，所得到的dna甲基化水平非常精确，能够绘制单碱基分辨率的细胞游离dna的甲基化图谱。
[0086]
本技术采用酶学转化法和单链建库法相结合的dna测序文库构建方法，其先将末端缺失或末端未缺失的双链dna变性为单链dna，然后采用酶学转化法将单链dna中未甲基化的胞嘧啶转化为尿嘧啶，之后在每段单链dna的3’末端加上截短接头序列，采用pcr延伸反应以截短接头序列为起点补全互补链，得到双链延伸产物。接着在该双链延伸产物的未添加截短接头序列的另一端连接上双链接头序列，在已补全的平头双链dna的每条链的5’端加上具有双链结构的测序接头序列(含有index序列)，通过indexingpcr进行扩增，以得到dna甲基化测序文库。
[0087]
由上述可知，本技术并未对末端缺失的双链dna进行末端修复，能够避免由于末端修复所带来的测序数据冗余现象。原因如下：如果先进行末端修复，那么在末端修复时无法确定缺失的3’末端或5’末端有几个碱基被修复，并且在测序时末端修复的碱基信息也被测
序仪计算在内，导致测序仪读取到不属于原始样本中的碱基数据，从而造成数据冗余现象。
[0088]
再者，因为测序时会将修复的碱基也计算在内，导致测得的平均甲基化水平低于真实的平均甲基化水平，这也会导致甲基化水平失真。本技术并非采用末端修复和双链建库法，而是采用单链建库法，不会产生上述的因为末端修复而造成的数据冗余现象和甲基化水平失真现象，所测得的碱基数据与真实的碱基数据一致。
[0089]
以下结合实施例对本技术作进一步的说明。
[0090]
实施例一
[0091]
本实施例提供了一种目标双链dna的获得方法，其包括如下步骤：
[0092]
1、提取目标样品中的dna，检测dna的浓度，并且采用电泳检测提取的dna条带的分子量是否属于目标dna的分子量。如果属于，则进行下一步。如果不属于，则重新提取dna。
[0093]
2、超声打断提取的dna，其具体包括：根据dna的浓度和电泳条带确定dna的上样量，一般打断起始量为150ng，总体积为100μl，置于打断仪中进行打断，将打断后的dna进行2％琼脂糖凝胶电泳检测，dna片段在200bp至300bp之间时，即得含有目标双链dna的dna测序样品。打断仪的程序为打断30s，然后暂停30s，共持续9个循环。
[0094]
当从组织中提取基因组dna，由于其片段较长，故需要采用上述的打断方法。而从血液中提取到的cfdna的片段一般位于120bp至200bp之间，例如为170bp，故无需上述的打断过程。
[0095]
本技术除了适用于具有平头端的双链dna之外，还能够适用于血液中游离并且具有不平整末端的cfdna和ctdna。上述两种dna在血液环境中容易被酶消化，从而3’末端的不平整现象要比来自于细胞核的组织dna更常见。
[0096]
实施例二
[0097]
本实施例提供了一种目标双链dna的甲基化保护方法，其包括如下步骤：
[0098]
1、加入内参dna：
[0099]
取26μl实施例1所得的打断后的dna测序样品，加入1μl打断后的λdna和1μl打断后puc19 dna(含有0.01-0.02ng)，混合后得到总体积为28μl的dna混合样本。
[0100]
本实施例中，加入的puc19 dna是完全甲基化，即其上所有的胞嘧啶(c)是完全甲基化的，puc19的甲基化率用于指示实验中对甲基化胞嘧啶的保护是否有效。如果最终分析结果中，puc19甲基化率大于98％时可以视为甲基化保护成功。因此，加入的puc19 dna可以作为内参。
[0101]
本实施例中，加入的λdna是完全未甲基化的，即其上所有的胞嘧啶(c)是完全不甲基化的。λdna的甲基化率可以指示实验中对未甲基化c的转化是否有效。如果最终分析λdna甲基化率低于1.5％时可以视为试验成功。因此，加入的λdna也可以作为内参。
[0102]
2、配制甲基化保护液：
[0103]
将10μl的tet2反应缓冲液(tet2 reaction buffer)、1μl的氧化补剂(oxidation supplement)、1μl的二硫苏糖醇(dithiothreitol，dtt)、1μl的氧化增强剂(oxidation enhancer)和4μl的tet2试剂在pcr薄壁管中混匀，得到17μl的甲基化保护液。
[0104]
其中，使用的试剂盒为neb公司货号为e7125的酶法甲基化转化模块试剂盒，即enzymatic methyl-seq conversion module。
[0105]
tet2酶和氧化增强剂作为甲基化保护试剂用于保护甲基化胞嘧啶上的甲基，其能
够通过化学反应改变甲基基团，使其对脱氨试剂惰性，而未甲基化的胞嘧啶上不存在甲基，故不会被氧化而仍然以胞嘧啶形式存在于dna中。dna上的甲基化胞嘧啶包括5-甲基胞嘧啶(5-methylcytosine，5mc)或5-羟甲基胞嘧啶(5-hydroxymethylcytosine，5hmc)。5-羟甲基胞嘧啶是5-甲基胞嘧啶的羟基化形式。tet2酶(ten eleven translocation)属于双加氧酶家族，其结构包括一个保守的c端的催化区和一个n端的调节区。c端的催化区为fe
2+
和/或α-酮戊二酸(α-kg)依赖型，是氧化5mc或5hmc的催化结构域。tet2酶和氧化增强剂的共同作用能够将5-甲基胞嘧啶(5mc)氧化为5-羟甲基胞嘧啶(5hmc)，并将5-羟甲基胞嘧啶(5hmc)氧化为5-甲酰基胞嘧啶(5-formylcytosine，5fc)或5-羧基胞嘧啶(5-carboxylcytosine，5cac)，这相当于在胞嘧啶上加上了甲酰基或羧基等保护基团，从而对脱氨试剂表现出惰性，在后继的反应过程中使得胞嘧啶不会由于发生脱氨而生成尿嘧啶。与5mc和5hmc相比，5-甲酰基胞嘧啶和5-羧基胞嘧啶在哺乳动物基因组中含量极低。氧化增强剂用于增加tet2酶的活性。
[0106]
3、进行氧化反应，具体包括：
[0107]
(3-1)、将总体积为28μl的dna混合样本与17μl的预混液相混合，加入5μl二价铁稀释液，立即震荡混匀和短暂离心，置于pcr仪中，在37℃反应1h。
[0108]
(3-2)、反应结束后，将反应液置于冰盒上，加入1μl的终止试剂(stop reagent)，震荡混匀后短时离心。
[0109]
(3-3)、将反应液置于pcr仪中，在37℃反应0.5h，此时，氧化反应已经进行完毕，得到了甲基受保护的双链dna。
[0110]
其中，上述的二价铁稀释液的配制方法包括：取1μl 500mm二价铁溶液(fe ii solution)，加入到1249μl nf水中。因为二价铁溶液容易被氧化，故需要在使用前配制，现配现用。二价铁离子是tet2酶的催化区所依赖的金属离子，为了最大化tet2酶的催化活性，故需要在氧化反应体系中加入二价铁离子。
[0111]
4、在氧化反应后进行纯化，得到纯化后的dna。
[0112]
实施例三
[0113]
实施例提供了一种dna变性的方法，其包括如下步骤：
[0114]
取16μl实施例二所得的纯化后的dna置于pcr反应管中，加入4μl甲酰胺，放入预热后的pcr仪中，盖上热盖，85℃孵育10min，孵育结束后，立即取下反应管，置于冰盒上保存。
[0115]
其中，甲酰胺作为变性剂用于使得纯化后的双链dna变性为单链dna。本技术先氧化反应保护甲基化的胞嘧啶，再变性为单链转化未甲基化的胞嘧啶，这是neb公司酶学转化方法设计决定的。
[0116]
实施例四
[0117]
本实施例提供了一种单链dna的转化方法，其包括如下步骤：
[0118]
1、配制脱氨反应液
[0119]
取1支pcr薄壁管，加入10μlapobec反应缓冲液(apobec reaction buffer)、1μl牛血清白蛋白(bovine serum albumin，bsa)、1μlapobec，然后加入68μl无核酸酶水(nuclease-free water，nf water)，充分震荡混匀后，配制成80μl脱氨反应液。
[0120]
2、进行脱氨反应
[0121]
用移液器吸取80μl脱氨反应液置于实施例三所得的已变性完成的单链dna中，充
分震荡混匀，并置于预热好的pcr仪上，于37℃孵育3h。
[0122]
在该步骤中，apobec蛋白的c端具有胞嘧啶脱氨酶区，能够使胞嘧啶脱去氨基变成尿嘧啶，因此，apobec作为胞嘧啶脱氨酶将实施例三所得的单链dna中未甲基化的胞嘧啶(c)转化为尿嘧啶(u)，而受保护胞嘧啶不能被apobec蛋白脱氨，故不会转化为尿嘧啶。脱氨反应后，得到了转化后单链dna。
[0123]
3、在脱氨后进行纯化，得到能够用于建库的转化后的单链dna。
[0124]
实施例五
[0125]
本实施例提供了一种利用转化后单链dna构建dna甲基化测序文库的方法，其包括如下步骤：
[0126]
1、dna重变性：
[0127]
取6.63μl实施例四所得的转化后的单链dna至新的0.2ml pcr管中，然后置于pcr仪中，盖上热盖，95℃孵育2min。孵育完成后，立即将pcr管放置在冰上孵育2min，并立即进行下一步的衔接反应。
[0128]
dna重变性的目的是使得下一步的衔接反应建立在单链dna(single stranded dna，ssdna)模板基础上。因为实施例四虽然得到了单链dna，但是在存放的过程中单链dna会发生复性，重新变成双链dna，故需要在衔接反应之间确保待衔接的dna处于完全单链dna的状态。
[0129]
2、配制衔接酶(adaptase)反应液
[0130]
取一个pcr薄壁管，加入1μl g1缓冲液(buffer g1)、1μl g2试剂(reagent g2)、0.63μl g3试剂(reagent g3)、0.25μl g4酶(enzyme g4)、0.25μl g5酶(enzyme g5)、0.25μl g6酶(enzyme g6)，配制成3.38μl的衔接酶反应液。
[0131]
本实施例使用的建库试剂盒为swift公司的货号30024的甲基化建库试剂盒swift accel-ngs methyl-seq dna library kit。建库流程完全按照产品说明书进行。
[0132]
本步骤的衔接酶反应液仅适用于单个变性后的dna样本。如果存在多个变性后的dna样本，可以同比例放大，混匀后分别加入单个dna样本中。
[0133]
3、进行衔接反应：
[0134]
将步骤2所得的3.38μl的衔接酶反应液加入步骤1所得的含有6.63μl变性完毕的单链dna的pcr管中，混合均匀。然后置于预热好的pcr仪中，盖上热盖，热盖的温度设定为105℃。在37℃孵育15min，然后在95℃孵育2min。孵育完毕后，即在转化后单链dna的3’端衔接上截短接头序列，从而得到单链接头产物。
[0135]
截短接头序列为swift公司的甲基化建库试剂盒swift accel-ngs methyl-seq dna library kit中的成分，其结构为部分双链，一端与样本单链dna片段连接，另一端带有pcr引物。本步骤中，5’端不加截短接头序列。在单链dna延伸后变成双链dna后可以在5’端直接加双链的测序接头序列。
[0136]
衔接反应进行完毕后，将dna样本置于冰上保存。
[0137]
4、配制延伸反应液：
[0138]
取一个pcr薄壁管，加入0.5μl的试剂y1(reagent y1)和10.5μl的酶y2(enzyme y2)，震荡混匀后，配成总体积为11μl的延伸反应液。
[0139]
本步骤的延伸反应液仅适用于单个变性后的dna样本。如果存在多个变性后的dna
样本，可以同比例放大，混匀后分别加入单个dna样本中。
[0140]
5、进行延伸反应：
[0141]
将步骤4所得的11μl的延伸反应液加入到步骤3所得的含有10.01μl衔接反应完毕后的dna样本的pcr管中，吹打混匀。将该pcr管置于提前预热的pcr仪中，盖上热盖，设置热盖的温度为105℃。在98℃下孵育1min，62℃下孵育1min，65℃下孵育5min。
[0142]
在延伸反应过程中，酶y2为dna聚合酶，能够利用底物dntp分子以截短接头序列为起点从5’端向3’端进行聚合反应，并最终聚合形成双链dna(即为双链延伸产物)。
[0143]
6、进行纯化步骤：
[0144]
延伸反应完成后，使用1.2倍磁珠(25.2μl)进行纯化，使用80％乙醇进行清洗，最后溶解在4.25μl low te中，备用，
[0145]
7、配制连接反应液：
[0146]
本步骤仍然使用swift公司试剂盒。
[0147]
取1个pcr管，按照顺序加入0.75μl的b1缓冲液(buffer b1)、2.5μl的b2试剂(reagent b2)、0.5μl的b3酶(enzyme b3)，吹打混匀后，配成总体积为3.75μl的连接反应液。
[0148]
其中，b3酶需要在使用之前加入，以防止酶的失活。
[0149]
本步骤的连接反应液仅适用于单个变性后的dna样本。如果存在多个变性后的dna样本，可以同比例放大，混匀后分别加入单个dna样本中。
[0150]
8、进行连接反应：
[0151]
从步骤6所得的纯化后的双链延伸产物中吸取3.75μl置于pcr管中，将3.75μl的连接反应液加入该pcr管中，吹打混匀。将该pcr管置于pcr仪中，打开热盖，在无热盖加热的条件下于25℃孵育15min。连接反应完成后，使用磁珠进行纯化，使用80％乙醇进行清洗，最后溶解在6μl low te中，备用。
[0152]
在本步骤中，磁珠的体积需要根据dna的样本类型确定。对于cfdna，样本体积为7.5μl时(3.75μl的延伸反应产物和3.75μl的连接反应液)，用于纯化的磁珠体积为cfdna体积的1.2倍，即为9μl。对于gdna，样本体积为7.5μl时，用于纯化的磁珠体积为cfdna体积的1.0倍，即为7.5μl。
[0153]
在本步骤中，连接反应用于将双链接头序列添加到双链dna的末端，得到双链接头产物。
[0154]
在本步骤中，b3酶为t4 dna连接酶(t4dna ligase)，其能在双链dna的5’末端添加双链的双链接头序列，而并非在3’末端添加双链接头序列。
[0155]
双链接头序列为上述的swift公司的甲基化建库试剂盒，属于一种双链接头，但并非颈环结构
[0156]
9、配制标记pcr(indexing pcr)反应液
[0157]
取pcr管，加入2.5μl的low edta te，2.5μl的r1缓冲液(buffer r1)、1μl的r2试剂(reagent r2)和0.25μl的r3酶(enzyme r3)，吹打均匀后配成6.25μl的标记pcr反应液。
[0158]
其中，r3酶需要在使用之前加入，以防止酶失活。
[0159]
10、进行标记pcr反应(indexing pcr)
[0160]
(10-1)、取步骤8所得的5μl纯化后的连接产物(含有双链接头产物)，加入6.25μl
的标记pcr反应液，吹打混合均匀。
[0161]
(10-2)、针对每个纯化后的连接产物，加入以下的标记引物(index primer)：1.25μl编号为cat.no.36024的index x、0.625μl编号为cat.no.38096的index d50x、0.625μl编号为cat.no.38096的index d7xx、以及1.25μl编号为cat.no.39096和390384的index u001-u096。
[0162]
上述的引物的编号为试剂盒中厂商命名的编号，针对同一个样品而言，上述引入需要同时加入。
[0163]
index u001-096为96个不同的index(index通常为8个碱基)。每个dna单链(下一步需要高温变性)的两端分别连接一对不同的index，用于唯一标识这个样本，因为测序时会有很多样本混在一起，需要靠每个样本上不同的index区分。
[0164]
(10-3)、之后置于pcr仪中，进行pcr反应，反应条件如下：98℃持续2min，进行1个循环；98℃变性过程持续10s、60℃退火过程持续20s、68℃延伸过程持续60s，依次进行13个循环；72℃持续5min，共进行1个循环，然后4℃保存。
[0165]
(10-4)、上述的pcr反应结束后，使用磁珠进行纯化，使用80％乙醇进行清洗，磁珠体积和洗脱体积如下所示：若dna种类为cfdna，以样本体积为15μl计算，磁珠体积为样本体积的0.8倍(即10.0μl)，洗脱体积为15μl，吸取体积为14μl。若dna种类为基因组dna(genomic dna，gdna)，以样本体积为12.5μl计算，磁珠体积为样本体积的0.85倍(即10.62μl)，洗脱体积为15μl，吸取体积为14μl。经过纯化后，得到dna甲基化测序文库。将该dna甲基化测序文库置于测序仪中，可以测得目标dna的甲基化水平。
[0166]
在经过上述的indexing pcr扩增之后，u碱基的配对碱基是a，a碱基配对碱基是t。因为pcr过程中使用的寡核苷酸是a、t、c、g，不含u，在经过n个循环(n≥13)过后，含u的原始链变成了仅有1/(2
13
)，可以忽略不计。由此，双链dna中的尿嘧啶(u)被转化为胸腺嘧啶(t)。
[0167]
在标记pcr反应完毕后，双链接头产物的5’端添加了具有双链结构的测序接头序列，该测序接头序列含有用于标记样本来源的index序列，从而得到能够被illumina测序仪读取的dna甲基化测序文库。每个dna甲基化测序文库的dna分子都包括截短接头序列、双链接头序列和测序接头序列(该测序接头序列含有index序列)。
[0168]
上述步骤的序号并不视为对上述步骤的顺序的唯一限定。
[0169]
实施例六
[0170]
本技术提供了一种dna甲基化水平的检测方法，其包括如下步骤：
[0171]
(1)、对目标双链dna按照上述的构建方法构建dna甲基化测序文库；
[0172]
(2)、采用illumina测序仪对所述dna甲基化测序文库进行测序，得到测序结果；
[0173]
(3)、对测序结果进行分析，得到dna甲基化水平。
[0174]
实验例
[0175]
本实验例对本技术的测序文库与现有的商用试剂盒采用的酶学法甲基化实验构建的测序文库的测序结果进行对比，其结果如图4所示。图4是本技术的方案相对现有技术的结果对比图。现有技术包括如下步骤：先在双链状态下使用未甲基化的dntp进行末端修复，再加a，加接头，进行酶学转化，并完成pcr扩增，从而得到甲基化测序文库。
[0176]
在图4中，read2通过illumina公司的novaseq600型号测序仪读取平均甲基化率变化。纵坐标表示read2在各读取位置上的平均甲基化率，横坐标从左至右的方向表示读长2
(read 2)从p5端至p7端的方向，横坐标的值表示测序仪读取的碱基位置。图4中的上方曲线表示本技术建立的文库测得的甲基化率，下方曲线表示现有酶学法建立的文库测得的甲基化率。本实验例的样本与现有酶学法所使用的样本相同。
[0177]
进入测序仪进行测序的每个dna小片段的长度约在200至300个碱基之间，甚至可能更短。测序仪在读每个dna小片段的时候，并不能从头读到尾，只能从头读150个碱基的长度，并从尾读150个碱基的长度。即测序仪分别从dna片段的两头读取，各读150个碱基长度，故把这些测序仪读取产生的短读段称为reads，单个短读段称为read。每条文库dna单链在二代测序仪上读取的时候，是先从p5端到p7端读150碱基，为read1；再从p7端到p5端读150个碱基，为read2。p5和p7是illumina测序仪的芯片上的接头名称，p5对应5'端，p7对应3'端。测序仪的芯片上的接头与测序文库的dna分子上的接头序列有互补性。
[0178]
从图4可知，针对相同样本，现有酶学法测得的平均甲基化率在1至33个碱基范围内与本技术测得的平均甲基化率基本相同，在34个碱基到133个碱基范围内存在显著不同。理论上，因为所有待测序分子在基因组上的位置分布是随机的，因此其平均甲基化率不应该随读取位置而变化。但是从图4可以看到，现有酶学法在越接近p7端(即进行末端修复的3’端方向)时，平均甲基化率越低，说明末端修复引入的未甲基化dntp对样本原始甲基化信号造成了污染，且在越靠近原始单链的3’端时，这种污染越明显。而本技术的方法则避免了这个问题，在各个读取位置的平均甲基化率上保持了良好的一致性。这能够说明：与现有的基于酶学转化法的双链建库法相比，本技术不会在末端引入新的甲基化水平，不会产生甲基化水平失真或降低现象。
[0179]
以上对本技术进行了详细介绍，本文中应用了具体个例对本技术的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本技术的方法及其核心思想；同时，对于本领域的技术人员，依据本技术的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本技术的限制。