一种新型多样本多片段DNA甲基化检测方法与流程

一种新型多样本多片段dna甲基化检测方法
技术领域
1.本发明属于基因检测技术领域，具体涉及一种新型多样本多片段dna甲基化检测方法。


背景技术：

2.dna甲基化可以调控基因表达，因此检测dna甲基化至少可以为解释基因表达的调控提供一些线索。其次，dna甲基化已证实参与一系列重要生物学过程，包括早期胚胎发育、基因组印记、x染色体失活、重复序列的沉默以及癌症的发生发展转移。此外，研究表明，dna甲基化一个非常重要的应用，是可以作为肿瘤的早期筛查和预后的生物标志物。
3.人类大约有30亿个碱基对，dna甲基化主要发生在cpg二核苷酸上。基因组中大约有2800万个cpg位点，这些cpg位点中，大约60～80％被甲基化，而在一些特定的区域，比如启动子，存在cpg位点富集的序列(cpg岛)，cpg岛通常未被甲基化。但是在肿瘤细胞中，整体甲基化水平降低至20～50％，与正常细胞相比，即可能有10～60％的cpg位点的甲基化状态发生了改变，特别是抑癌基因上，也就是大约280～1680万个cpg位点。
4.现有的dna甲基化检测方法主要包括全基因组的dna甲基化检测和特异性位点甲基化检测。全基因组的dna甲基化检测主要用于研究层面、寻找差异化甲基化位点等，又可分为基于高通量测序和基于芯片的全基因组dna甲基化检测方法。特异性位点甲基化检测主要用于转化医学层面，包括甲基化特异性pcr(msp)、亚硫酸氢盐测序pcr(bsp)、焦磷酸测序和质谱检测。
5.在实际临床应用和操作过程中，检测基质中的dna含量一般较少，例如在癌症早筛检测中，血浆中的ctdna(circulating tumor dna)含量，在游离dna(cfdna)的占比较低(可低至0.01％)。cfdna的浓度在健康个体血浆中，每毫升约含有0至100ng，平均浓度大约为30ng/ml。而且ctdna片段长度，一般约为166bp(死亡细胞裂解)或者150-250bp(癌细胞外囊泡分泌)。在此种条件下，传统的通过杂交的方式捕获dna后再进行高通量测序的方法灵敏性和特异性差。且基于全基因组的甲基化测序方式，流程复杂、后续分析复杂、捕获成本较高。其他通过pcr的方法，比如bsp的方法，因为引物不包含cpg位点，虽然可以测到甲基化，但有两点问题，一是由于不能包含cpg位点，一般的c经过bs转化变成t，所以在实际操作中，特别是cpg岛里面设计引物较困难，存在gc含量低，pcr实际效率差，非特异性产物过多的情况；二是针对于痕量的甲基化片段，因为pcr没有选择性，非甲基化的信号往往会将这部分信号湮没，比如一段序列可能正常细胞分泌了99％但都没有甲基化，而癌细胞就分泌了1％发生甲基化了，所以如果用bsp可能99％因为无法区分噪音还是真实信号，1％这部分信息就会丢失掉。使用msp，或芯片，仅能覆盖有限的cpg位点(一般目标只有1个位点)，灵活度低。一般cpg岛区域有较多位点会被甲基化，传统方法不能够知道一段序列内的多个位点的cpg情况。


技术实现要素：

6.基于此，本发明提供一种新型多样本多片段dna甲基化检测方法，所述检测方法可以同时对多样本多基因的dna甲基化进行检测，整个检测流程和后续的数据分析步骤简单科学，对起始dna量要求低，非常适用于低拷贝数的甲基化片段检测，且具有灵敏度高、特异性好、检测限低以及检测成本低等优点。
7.具体技术方案为：
8.一种多样本多片段dna甲基化检测方法，包括以下步骤：
9.s1、提取样本基因组dna或游离dna(cell-free dna)，对提取的dna使用重亚硫酸盐进行处理，然后使用特异性引物对处理好的dna进行依赖甲基化多重pcr(methylation dependent multiplex pcr)；其中针对不同基因或不同甲基化位点使用不同的特异性pcr引物；
10.s2、多重pcr产物处理和纯化：对步骤s1获得的多重pcr产物用核酸内切酶或者dna片段分选磁珠进行处理，然后进行纯化；
11.s3、对步骤s2获得的纯化产物的甲基化水平进行检测。
12.本发明所述依赖甲基化多重pcr扩增(methylation dependent multiplex pcr)是指在同一个pcr扩增体系中同时针对多个感兴趣甲基化位点及其附近50-300bp区域进行的pcr扩增，其中感兴趣甲基化位点可位于引物上或扩增子内部。
13.本发明所述特异性pcr引物是指针对感兴趣的甲基化位点，或其旁边的cpg位点进行设计的引物；所述特异性pcr引物包含亚硫酸氢钠处理后成功转化的甲基化的cpg位点，以保证对甲基化dna进行特异性扩增。
14.本发明检测方法可以同时对多样本多基因的dna甲基化水平进行检测，针对不同基因或不同甲基化位点使用不同的特异性pcr引物，针对不同的样本在检测时使用特异的标签序列来进行区别。
15.在其中一些实施例中，步骤s1所述依赖甲基化多重pcr的反应条件为：98℃ 30s-5min；10-25个循环(98℃ 15s，60
±
10℃ 15s-10min，68-72℃ 15s-5min)；68-72℃ 0-15min。
16.在其中一些实施例中，步骤s1所述依赖甲基化多重pcr的反应优化条件为：98℃ 30s-5min；10-25个循环(98℃ 15s，60
±
5℃ 15s-10min，72℃ 15s-5min)；72℃ 0-15min。
17.在其中一些实施例中，步骤s1所述依赖甲基化多重pcr的反应优化条件为：98℃ 30s-5min；10-25个循环(98℃ 15s，60
±
5℃ 15s-5min，72℃ 15s-5min)；72℃ 0-15min。
18.在其中一些实施例中，步骤s1所述依赖甲基化多重pcr的扩增方法为降落pcr法。发明人经过研究法发现，在本发明检测体系中使用降落pcr法进行依赖甲基化多重pcr，能有效提高依赖甲基化多重pcr反应的特异性，降低非特异性扩增产物。
19.在其中一些实施例中，步骤s1所述依赖甲基化多重pcr的反应条件为：98℃ 30s；5-10个循环(98℃ 15s，65
±
3℃(每个循环下降0.2-0.8℃)15s，72℃ 15s)；10-25个循环(98℃ 15s，60
±
10℃ 15s-10min，68-72℃ 15s-5min)；72℃ 0-15min。
20.在其中一些实施例中，步骤s1所述依赖甲基化多重pcr的反应条件为：98℃ 30s；5-10个循环(98℃ 15s，65
±
3℃(每个循环后下降0.2-0.8℃)15s，72℃ 15s)；10-25个循环(98℃ 15s，60
±
5℃ 15s-10min，68-72℃ 15s-5min)；72℃ 0-15min。
21.在其中一些实施例中，步骤s1所述依赖甲基化多重pcr的反应条件为：98℃ 30s；5-10个循环(98℃ 15s，65
±
3℃(每个循环后下降0.2-0.8℃)15s，72℃ 15s)；10-25个循环(98℃ 15s，60
±
5℃ 15s-5min，68-72℃ 15s-5min)；72℃ 0-15min。
22.在其中一些实施例中，优选步骤s1所述依赖甲基化多重pcr的反应条件为：98℃ 30s；5-10个循环(98℃ 15s，65
±
3℃(每个循环后下降0.2-0.8℃)15s，72℃ 15s)；15-20个循环(98℃ 15s，60
±
3℃ 15s，72℃ 15s)；72℃ 15min。
23.在其中一些实施例中，更优选地，步骤s1所述依赖甲基化多重pcr的反应条件为：98℃30s；10个循环(98℃ 15s，65
±
3℃(每个循环下降0.5℃)15s，72℃ 15s)；15个循环(98℃ 15s，60
±
3℃ 15s，72℃ 15s)；72℃ 15min。发明人经过研究发现，此多重pcr扩增反应条件下能成功检测到目标基因的数目最多。
24.在其中一些实施例中，步骤s1所述依赖甲基化多重pcr的反应酶为phusion dna聚合酶、q5酶、hieff酶、kapa dna聚合酶、pfu酶、superfi酶、相应可兼容尿嘧啶的聚合酶、优选为phusionu酶、kapau酶、q5u酶、hieffu酶。
25.在其中一些实施例中，优选步骤s1所述依赖甲基化多重pcr的反应酶为phusion dna聚合酶。
26.在其中一些实施例中，更优选地，所述依赖甲基化多重pcr的反应酶为phusion u hot start dna聚合酶。使用phusion u hot start dna聚合酶能更好地提高本发明检测方法依赖甲基化多重pcr反应的特异性，降低非特异性扩增产物。
27.在其中一些实施例中，所述降落pcr法为巢式pcr法。
28.在其中一些实施例中，所述步骤s1中同时对不同基因和参考序列进行依赖甲基化多重pcr。
29.在其中一些实施例中，所述参考序列选自参考序列1～4中的至少一种：
30.参考序列1：以seq id no.25和seq id no.26为特异性pcr引物进行pcr扩增的片段在epha3基因中所对应的序列；
31.参考序列2：以seq id no.27和seq id no.28为特异性pcr引物进行pcr扩增的片段在kbtbd4基因中所对应的序列；
32.参考序列3：以seq id no.29和seq id no.30为特异性pcr引物进行pcr扩增的片段在plekhf1基因中所对应的序列；
33.参考序列4：以seq id no.31和seq id no.32为特异性pcr引物进行pcr扩增的片段在syt10基因中所对应的序列。
34.在其中一些实施例中，所述参考序列选自参考序列1～4中的至少两种。使用至少两种所述参考序列可以使检测结果更稳定，使不同样本之间的检测结果比较更准确。
35.在其中一些实施例中，所述不同基因的序列选自seq id no.1～seq id no.8所示序列中的至少一种。
36.在其中一些实施例中，针对seq id no.1～seq id no.8所示序列中的至少一种的依赖甲基化多重pcr扩增的反应条件为：98℃ 30s；10个循环(98℃ 15s，65
±
3℃(每个循环下降0.2～0.8℃)15s，72℃ 15s)；15～20个循环(98℃ 15s，60
±
3℃ 15s，72℃ 15s)；72℃ 15min。
37.在其中一些实施例中，优选针对seq id no.1～seq id no.8所示序列中的至少一
种的依赖甲基化多重pcr扩增的反应条件为：98℃ 30s；10个循环(98℃ 15s，65
±
3℃(每个循环下降0.5℃)15s，72℃ 15s)；15个循环(98℃ 15s，60
±
3℃ 15s，72℃ 15s)；72℃ 15min。
38.在其中一些实施例中，步骤s2所述核酸内切酶处理条件为：37
±
1℃处理10-15min，处理体系中核酸内切酶的终浓度为1-10u/ul。
39.在其中一些实施例中，步骤s2所述核酸内切酶处理条件为：37
±
1℃处理10-15min，处理体系中核酸内切酶的终浓度为3.5-4.5u/ul。
40.在其中一些实施例中，优选步骤s2所述核酸内切酶处理条件为：37
±
1℃处理10min，处理体系中核酸内切酶的终浓度为4u/ul。
41.在其中一些实施例中，所述核酸内切酶为t4核酸内切酶。
42.在其中一些实施例中，步骤s2所述dna片段分选磁珠为xp磁珠。
43.在其中一些实施例中，步骤s2所述纯化为使用磁珠法进行纯化。
44.在其中一些实施例中，所述磁珠法进行纯化使用的磁珠为xp磁珠。发明人经过研究发现，在本发明检测方法中，使用xp磁珠可获得更好的纯化效果。
45.在其中一些实施例中，步骤s3利用测序法对步骤s2获得的纯化产物的甲基化水平进行检测。
46.在其中一些实施例中，步骤s3所述测序法包括以下步骤：(1)将步骤s2获得的纯化产物同时进行3'末端修复和3'末端添加碱基a；(2)将步骤(1)获得的产物与测序接头连接；(3)使用接头引物对步骤(2)获得的连接产物进行pcr扩增，得到测序文库；(4)将不同样本制备所得文库按照相同的摩尔数混合，测序。
47.在其中一些实施例中，所述步骤(3)中pcr扩增的循环数为3-7。优选地，所述步骤(3)中pcr扩增的循环数为7。
48.在其中一些实施例中，所述步骤(3)中还包括对测序文库进行纯化和质检。
49.在其中一些实施例中，使用xp磁珠对测序文库进行纯化。
50.在其中一些实施例中，所述测序法包括二代测序技术与三代测序技术。二代测序技术原理为大规模平行测序(massive parallel sequencing,mps)，三代测序技术为单分子测序技术。
51.在其中一些实施例中，所述二代测序技术可以是基于dna聚合酶合成测序技术(sequencing by synthesis technology，sbs)、基于dna连接酶连接测序技术(sequencing by ligation technology,sbl)。
52.在其中一些实施例中，所述三代测序技术可以是单分子实时荧光测序技术、纳米孔测序技术、纳米门测序技术、基于dna水解测序技术。
53.本发明将上述步骤s3通过测序法对步骤s2获得的纯化产物的甲基化水平进行检测的多样本多片段dna甲基化检测方法命名为依赖甲基化的扩增和测序法(methylation-dependent amplification and sequencing,medas)。
54.在其中一些实施例中，步骤s3利用通用阵列法对步骤s2获得的纯化产物的甲基化水平进行检测。所述通用阵列(universal array)是指将已知的目标dna片段(可以是甲基化片段)和未知的核酸序列之间的一方以有序的阵列固定到特定载体上，通过序列互补杂交原理，并将杂交的结果以荧光技术和模式识别分析来检测的技术。检测时不同样本使用
特异的标签序列进行区别。
55.在其中一些实施例中，所述通用阵列可以是微阵列技术(microarray)、珠阵列技术(bead array)、xmap(multi-analyte profiling)技术、ncounter技术。
56.在其中一些实施例中，所述样本为生物体液、细胞或组织。
57.在其中一些实施例中，所述生物体液为人体在正常或病理状态下各个器官和组织分泌的液体。
58.在其中一些实施例中，所述生物体液为血液、尿液、唾液、汗液、脑脊液、胸腔积水、腹积水。
59.在其中一些实施例中，所述生物体液优选为血清、血浆、玻璃体、痰、尿、眼泪、汗液、唾液等。
60.本发明还提供了一种dna甲基化检测试剂盒。
61.具体技术方案为：
62.一种dna甲基化检测试剂盒，所述试剂盒包括针对seq id no.1～seq id no.8所示序列中的至少一种序列的特异性pcr引物，所述引物为：
63.针对seq id no.1的特异性pcr引物为seq id no.9所示的上游引物和seq id no.10所示的下游引物；
64.和/或，针对seq id no.2的特异性pcr引物为seq id no.11所示的上游引物和seq id no.12所示的下游引物；
65.和/或，针对seq id no.3的特异性pcr引物为seq id no.13所示的上游引物和seq id no.14所示的下游引物；
66.和/或，针对seq id no.4的特异性pcr引物为seq id no.15所示的上游引物和seq id no.16所示的下游引物；
67.和/或，针对seq id no.5的特异性pcr引物为seq id no.17所示的上游引物和seq id no.18所示的下游引物；
68.和/或，针对seq id no.6的特异性pcr引物为seq id no.19所示的上游引物和seq id no.20所示的下游引物；
69.和/或，针对seq id no.7的特异性pcr引物为seq id no.21所示的上游引物和seq id no.22所示的下游引物；
70.和/或，针对seq id no.8的特异性pcr引物为seq id no.23所示的上游引物和seq id no.24所示的下游引物。
71.在其中一些实施例中，所述试剂盒还包括以下特异性pcr引物中的至少一种：
72.由seq id no.25所示的上游引物和seq id no.26所示的下游引物组成的特异性pcr引物；
73.和/或，由seq id no.27所示的上游引物和seq id no.28所示的下游引物组成的特异性pcr引物；
74.和/或，由seq id no.29所示的上游引物和seq id no.30所示的下游引物组成的特异性pcr引物；
75.和/或，由seq id no.31所示的上游引物和seq id no.32所示的下游引物组成的特异性pcr引物。
76.在其中一些实施例中，所述试剂盒包括上述特异性pcr引物中的至少两种。
77.在其中一些实施例中，所述试剂盒包括上述四种特异性pcr引物。
78.在其中一些实施例中，所述试剂盒还包括phusion dna聚合酶、t4核酸内切酶、dna纯化磁珠。
79.在其中一些实施例中，优选所述phusion dna聚合酶为phusion u hot start dna聚合酶。
80.在其中一些实施例中，优选所述dna纯化磁珠为xp磁珠。
81.在其中一些实施例中，所述试剂盒还包括测序文库构建试剂。
82.与现有技术相比，本发明具有以下有益效果：
83.本发明所述检测方法引入依赖甲基化多重pcr技术，能快速对感兴趣的甲基化区域进行特异性扩增富集，并且还对依赖甲基化多重pcr技术进行了优化：先使用高保真聚合酶进行依赖甲基化多重pcr扩增，再进一步利用核酸内切酶或dna片段分选磁珠处理多重pcr产物，然后进行纯化。发明人发现，通过上述优化能有效去除依赖甲基化多重pcr扩增产物中的非特异性产物，避免非特异性产物对后续检测的影响，例如避免非特异性产物在测序过程中产生杂信号和占据测序深度，有效提高后续检测结果的准确性，降低检测限。尤其是当使用phusion u dna聚合酶进行依赖甲基化多重pcr时，可以更好地提高本发明依赖甲基化多重pcr反应的特异性。此外，本发明检测方法还使用了合适的参考序列，通过参考序列来对不同甲基化比例的样本的检测结果进行标准化，能平行比较出甲基化的差别，进一步提高检测结果的准确性和可靠性。
84.本发明通过引入依赖甲基化多重pcr技术并对依赖甲基化多重pcr扩增方法和产物处理方法进行优化，克服了现有甲基化检测方法对起始dna含量要求高，以及现有多重pcr扩增产物不适用于直接进行后续高通量检测的缺陷，具有灵敏度高、特异性好、准确性高、检测限低和检测成本低的优点，可以检测出0.05％甲基化水平，且检测流程和数据分析步骤更简单更科学，尤其适用于对低拷贝数的甲基化片段进行检测。
附图说明
85.图1本发明检测方法的流程示意图。
86.图2为实施例3中seq id no.5所示序列的两个cpg site的信息。
87.图3为实施例3中本发明检测方法对seq id no.5所示序列的两个cpg site的loq检测结果图。
88.图4为实施例3中亚硫酸氢盐修饰后测序法对seq id no.5所示序列的两个cpg site的loq检测结果图。
89.图5为实施例4中本发明检测方法和亚硫酸氢盐修饰后测序法对样本cpg位点1ch5_40681550的检测结果图。
90.图6为实施例4中本发明检测方法和亚硫酸氢盐修饰后测序法对样本cpg位点2ch5_40681569的检测结果图。
91.图7为使用q5u酶和phusion u酶进行依赖甲基化多重pcr扩增获得的产物中二聚体比例检测结果图。
92.图8为使用传统pcr法和touchdown pcr法进行依赖甲基化多重pcr扩增获得的产
物中二聚体比例和目标测序量比例检测结果图。
93.图9为使用xp beads和smart beads纯化依赖甲基化多重pcr扩增产物后纯化产物中二聚体比例检测结果图。
94.图10为使用xp beads和column纯化依赖甲基化多重pcr扩增产物后纯化产物中二聚体比例和目标测序量比例检测结果图。
95.图11为使用不同多重pcr扩增循环数进行依赖甲基化多重pcr扩增的检测效果对比图。
96.图12为实施例5中17个靶点在内的不同甲基化浓度的标品上的信号进行很好的线性拟合结果示意图。
97.图13为实施例6中的roc curve结果示意图。
98.图14为实施例7中的roc curve结果示意图。
具体实施方式
99.本发明下列实施例中未注明具体条件的实验方法，通常按照常规条件，例如sambrook等人，分子克隆：实验室手册(new york:cold spring harbor laboratory press,1989)中所述的条件，或按照制造厂商所建议的条件。实施例中所用到的各种常用化学试剂，均为市售产品。
100.除非另有定义，本发明所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不用于限制本发明。
101.本发明的术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤的过程、方法、装置、产品或设备没有限定于已列出的步骤或模块，而是可选地还包括没有列出的步骤，或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤。
102.在本发明中提及的“多个”是指两个或两个以上。“和/或”，描述关联对象的关联关系，表示可以存在三种关系，例如，a和/或b，可以表示：单独存在a，同时存在a和b，单独存在b这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。
103.本发明所述测序，包括二代测序技术与三代测序技术。二代测序技术原理为大规模平行测序(massive parallel sequencing,mps)，在其中一些实施例中，二代测序技术可以是：1.基于dna聚合酶合成测序技术(sequencing by synthesis technology，sbs)：代表公司为illumina(可逆终止测序，reversible terminator sequencing)、thermo fisher/life technologies(ion torrent)，genapsys，罗氏诊断(454焦磷酸测序)；2.基于dna连接酶连接测序技术(sequencing by ligation technology,sbl)，代表公司为华大基因/complete genomics(复合探针-锚定分子连接，cpal)、thermo fisher/applied biosystems(sequencing by oligonucleotide ligation and detection,solid)。三代测序技术为单分子测序技术，在其中一些实施例中，三代测序技术可以是：单分子实时荧光测序技术(smrt,pacific biosciences)、纳米孔测序技术[oxford nanopore technologies(ont),genia technologies and stratos genomics(罗氏诊断)]、纳米门测序技术(nanogate,quantum biosystems)、基于dna水解测序技术(sequencing by de-synthesis,
pyrophosphorolysis,base4)。
[0104]
本发明所述通用阵列(universal array)是指将已知的目标dna片段(可以是甲基化片段)和未知的核酸序列之间的一方以有序的阵列固定到特定载体上，通过序列互补杂交原理，并将杂交的结果以荧光技术和模式识别分析来检测的技术。在其中一些实施例中，通用阵列可以是：1.微阵列技术(microarray)，代表公司为affymetrix(thermo fisher)/agilent等；2.珠阵列技术(bead array)，代表公司为illumina等；3.xmap(multi-analyte profiling)技术,代表公司为luminex等；4.ncounter技术,代表公司为nanostring等。
[0105]
本发明提供的检测方法的检测流程图如图1所示。
[0106]
实施例1
[0107]
本实施例一种多样本多片段dna甲基化检测方法，包括以下步骤：
[0108]
s1、提取样本基因组dna或游离dna，对提取的dna使用采用ezdna methylation-gold(zymo)试剂盒进行处理，然后使用特异性引物对处理好的dna进行依赖甲基化多重pcr(methylation dependent multiplex pcr)；pcr反应使用的酶为高保真酶，可选用phusion dna聚合酶(例如phusion hot start ii dna polymerase，phusion u hot start dna polymerase)、q5酶(例如hot start high-fidelity dna polymerase,hot start high-fidelity dna polymerase)、hieff酶(例如hg热启动多重pcr酶)、kapa dna聚合酶(例如kapa hifi uracil kit，kapa2g快速热启动dna聚合酶)、pfu酶或superfi酶，本实施例优选phusion u hot start dna polymerase；pcr反应扩增方法为降落pcr法；其中针对不同基因或不同甲基化位点使用不同的特异性pcr引物；
[0109]
s2、多重pcr产物处理和纯化：可选方法1：对步骤s1获得的多重pcr产物用t4核酸内切酶进行处理，处理条件为：37
±
1℃处理10-15min，处理体系中酶的终浓度为1～10u/ul(优选处理条件为：37
±
1℃处理10min，处理体系中酶的终浓度为4u/ul)；将酶处理后的pcr产物进行纯化，所述纯化方法为xp磁珠纯化法；可选方法2：直接使用dna片段分选磁珠进行处理和纯化，其中磁珠为xp磁珠。
[0110]
s3、利用测序法对步骤s2获得的纯化产物的甲基化水平进行检测：(1)将步骤s2获得的纯化产物同时进行3'末端修复和3'末端添加碱基a；(2)将步骤(1)获得的产物与测序接头连接；(3)使用接头引物对步骤(2)获得的连接产物进行pcr扩增，pcr循环数为3-7轮(优选为7轮)，对扩增产物进行纯化，获得测序文库，对测序文库进行质检和定量；(4)将不同样本制备所得文库按照相同的摩尔数混合，利用illuminamiseq/miseqdx/nextseq/nextseqdx平台对文库进行测序。
[0111]
本实施例提供的检测方法命名为依赖甲基化的扩增和测序法(methylation-dependent amplification and sequencing,medas)，能同时对多样本多基因的dna甲基化进行检测，整个检测流程和后续的数据分析步骤简单科学，对起始dna量要求低，非常适用于低拷贝数的甲基化片段检测，且具有灵敏度高、特异性好、检测限低以及检测成本低等优点。表1为本发明medas法与现有的甲基化检测方法的特征比较：
[0112]
表1本发明medas法与现有的甲基化检测方法的比较
[0113][0114][0115]
表中，medas代表methylation-dependent amplification and sequencing；msp代表methylation specific pcr；bsp代表bisulfite sequencing pcr。
[0116]
从表1可知，与现有技术相比，本发明提供的medas技术具有通量高、灵敏度高、特异性好、检测限低以及检测成本低等优点。
[0117]
实施例2
[0118]
本实施例一种dna甲基化检测试剂盒，包括针对seq id no.1～seq id no.8所示序列的甲基化水平检测试剂，所述seq id no.1～seq id no.8所示序列如表2所示：
[0119]
表2 seq id no.1～seq id no.8所示序列
[0120][0121][0122]
上述试剂盒包含以下组分：
[0123]
(1)针对seq id no.1～seq id no.8所示序列的特异性pcr引物，所述特异性pcr引物具体如表3所示：
[0124]
表3针对seq id no.1～seq id no.8所示序列的特异性pcr引物
[0125][0126][0127]
(2)针对参考序列1-4(ref1-4)的特异性pcr引物：
[0128]
参考序列1：以seq id no.25和seq id no.26为特异性pcr引物进行pcr扩增的片段在epha3基因中所对应的序列；
[0129]
参考序列2：以seq id no.27和seq id no.28为特异性pcr引物进行pcr扩增的片段在kbtbd4基因中所对应的序列；
[0130]
参考序列3：以seq id no.29和seq id no.30为特异性pcr引物进行pcr扩增的片段在plekhf1基因中所对应的序列；
[0131]
参考序列4：以seq id no.31和seq id no.32为特异性pcr引物进行pcr扩增的片段在syt10基因中所对应的序列。
[0132]
具体如表4所示：
[0133]
表4参考序列1-4具体的特异性pcr引物信息
[0134][0135]
(3)phusion u hot start dna polymerase、t4核酸内切酶、agencourt ampure xp beads。
[0136]
上述试剂盒提供的针对seq id no.1～seq id no.8所示序列的特异性pcr引物是发明人经过大量研究和分析后设计得到的，所述检测引物既能保证依赖甲基化多重pcr扩增的灵敏度和特异性，又不会影响后续测序的深度，在2个拷贝～0.05％甲基化检测限上检测到甲基化片段需要的测序深度约小于现有方法(亚硫酸氢盐修饰后测序法)2-3个数量级。
[0137]
上述试剂盒提供的参考序列1-4(ref1-4)是发明人经过优化后得出的，在不同样本中表达稳定且不受甲基化水平的影响，非常适合用于对测序数据进行标准化处理，使不同样本间dna的甲基化水平具有可比性，提高比较结果的准确性和可靠性。
[0138]
实施例3
[0139]
本实施例对本发明dna甲基化检测方法(medas)和亚硫酸氢盐修饰后测序法(bsp)的检测效果进行比较。分别利用本发明dna甲基化检测方法和亚硫酸氢盐修饰后测序法(bsp)对实施例2所述seq id no.1～seq id no.8所示序列的甲基化情况进行检测。本发明dna甲基化检测方法使用实施例2所述试剂盒，具体步骤如下：
[0140]
s1、获得dna标准品，对dna标准品使用重亚硫酸盐进行处理，然后对处理好的dna标准品进行依赖甲基化多重pcr扩增
[0141]
(1)制备cfdna mock标准品
[0142]
1)genomic dna定量
‑‑
qubit hs
[0143]
a)qubit溶液配制：
[0144]
表5 qubit溶液配方
[0145]
配方体积μlqubit dsdna hs buffer(qb)199qubit dsdna hs reagent(qa)1总体积200
[0146]
b)genomic dna样品的稀释：取1μl dna样品 9μl eb，将dna样品稀释10倍，再取1μl稀释后的dna进行定量。
[0147]
2)genomic dna打断
[0148]
将genomic dna用eb稀释为约10ng/μl后，取300μl，分别向3支配套打断管(规格为130μl)中各加入50μl(重复两次)，用以下程序进行打断：
[0149]
表6 cfdna打断程序
[0150][0151][0152]
取3μl打断后dna放4℃保存，待进行浓度测定和2100bioanalyzer质检。
[0153]
3)打断后dna纯化与保存
[0154]
a)80％乙醇制备
[0155]
b)纯化步骤：原倍ampure磁珠室温平衡30min
→
取240μl与打断后的所有genomic dna于1.5ml的lobind管充分涡旋混匀，简单离心
→
室温孵育5min
→
上磁力架3-4min至溶液澄清
→
将上清转移至新的1.5ml lobind管
→
重新加入240μl原倍ampure磁珠
→
充分涡旋混匀，简单离心
→
室温孵育5min
→
上磁力架4min至溶液澄清，吸弃上清
→
加入1.5ml 80％乙醇洗涤30s，弃上清
→
加入1.5ml 80％乙醇洗涤30s，弃上清
→
离心1min，上架至液体澄清，吸余液
→
开盖晾干至磁珠表面无反光
→
加100μl eb洗脱，混匀，简单离心，室温孵育2min
→
上磁力架吸附2min至液体澄清，收集上清至新的1.5ml lobind管；重复一次洗脱步骤，dna样品终体积为200μl。
[0156]
c)保存：-20℃保存。
[0157]
4)纯化后dna定量和质检
[0158]
质检：取1μl纯化后的cfdna mock样品进行2100bioanalyzer质检；定量：取1μl进行qubit hs定量。
[0159]
(2)混合制备梯度浓度的dna标准品底物
[0160]
1)重亚硫酸盐处理基因组dna制备100％methylation control
[0161]
采用ezdna methylation-gold(zymo)试剂盒处理hg dna以及cfdna mock，具体步骤如下：
[0162]
a)相关buffer制备：
[0163]
a.ct conversion reagent的制备：添加900ul水，50ul的μ-dissolving buffer，和300ul的m-dilution buffer到一管ct conversion reagent中，在室温下溶解并且在摇床上摇动i0min；
[0164]
b.m-wash buffer的制备：加入24ml无水乙醇至m-wash buffer瓶中，并在瓶盖上做好标记；
[0165]
b)在pcr管中添加：
[0166]
a.dna量500ng，根据浓度算好体积，用水来补至20μl；
[0167]
b.130ul的ct conversion reagent：轻弹试管或移液器操作来混合样品；
[0168]
c)将样品管放到循环变温器并按以下步骤操作：98℃放置10min，64℃放置14h；
[0169]
d)将柱子放入collection tube中，并加入600ul的m-binding buffer，并将步骤2的样品加入到柱子中，盖上盖子将柱颠倒数次来混合样品；
[0170]
e)全速(》10,000x g)离心30s，弃废液；
[0171]
f)添加200ul的m-wash buffer到柱中，全速离心30s；
[0172]
g)添加200ul的m-desphonation buffer到柱中并且在室温(20-30℃)下放置15-20min；
[0173]
h)全速离心30s，弃废液；
[0174]
i)加200ul的m-wash buffer到柱中，全速离心30s；
[0175]
j)再添加200ul的m-wash buffer，并且离心30s；
[0176]
h)直接添加25ul的m-elution buffer到柱基质中，将柱放置在1.5ml的管中，全速离心来洗脱dna，得到最终转化产物(-20℃暂时储存，可保存1周)；
[0177]
i)通过测序分析，通过7-10个cpg、chg、chh位点判断标准品中甲基化转化程度。
[0178]
2)通过0％methylation control和100％methylation control混合制备梯度标准品：0％、0.05％、0.1％、0.25％、0.5％、1％、100％。
[0179]
(3)针对实施例2所述seq id no.1～seq id no.8所示序列和参考序列1-4分别进行本发明依赖甲基化多重pcr和multiplex bisulfide-specific pcr(多重bsp)。其中，本发明针对seq id no.1～seq id no.8所示序列的依赖甲基化多重pcr引物如实施例2所示，针对seq id no.1～seq id no.8所示序列的多重bsp引物如表7所示，针对参考序列1-4的引物如实施例2所示。
[0180]
表7多重bsp引物
[0181]
[0182][0183]
本发明依赖甲基化多重pcr扩增体系如表8所示，扩增程序如表9所示：
[0184]
表8本发明依赖甲基化多重pcr扩增体系
[0185][0186]
表9本发明依赖甲基化多重pcr扩增程序(touchdown法)
[0187]
[0188][0189]
bsp检测体系如表10所示，检测程序如表11所示：
[0190]
表10多重bsp检测体系
[0191][0192]
表11多重bsp检测程序
[0193][0194]
s2、pcr产物处理和纯化
[0195]
(1)对上述的依赖甲基化多重pcr产物用t4核酸内切酶进行处理，处理条件为：37
±
1℃处理10min，处理体系终t4核酸内切酶的终浓度为4u/ul；
[0196]
(2)将t4核酸内切酶处理后的不同样本的pcr产物分别使用agencourt ampure xp beads纯化，具体步骤如下：
[0197]
a)将ampure xp beads置于室温下平衡至少30min；
[0198]
b)向1.5ml离心管中加入50ng待建库的dna，再加入90ul(1.8x)ampure xp beads，在vortex上混勾并室温放置5min；
[0199]
c)将离心管放在磁力架上直至溶液澄清(大约3-5min)；
[0200]
d)保持离心管在磁力架上，小心弃掉离心管中的上清；
[0201]
e)继续保持离心管在磁力架上，加入200μl的80％乙醇；
[0202]
f)将管子放置1min，使所有beads沉淀，弃掉乙醇，重复e) f)一次；
[0203]
g)简单离心，上架至液体澄清，吸余液，在37℃加热模块中干燥样品5min或直到残余的乙醇完全消失；
[0204]
h)加55ul eb洗脱，混匀，室温孵育5min；
[0205]
j)简单离心，上架吸附3-5min至液体澄清，收集上清(共50ul，挑选个别样本2ul测qubit后再补齐50ul)。
[0206]
s3、利用测序法对步骤s2获得的纯化产物的甲基化水平进行检测
[0207]
使用neb文库制备试剂盒z1901s/l进行文库构建，不同样本使用特异的标签序列进行区别，具体步骤如下：
[0208]
(1)将步骤s2获得的纯化产物同时进行3'末端修复和3'末端添加碱基a。
[0209]
将上一步得到的dna纯化产物按下表在1.5ml离心管中配制末端修复和加a反应体系：
[0210]
表12末端修复和加a反应体系
[0211][0212]
反应程序如下：20℃,30min
→
65℃,30min
→
4℃,hold，热盖：85℃，体积60μl。
[0213]
(2)将3'末端添加碱基a的产物连接测序接头，以便获得测序接头连接产物。
[0214]
将上一步得到的产物直接按下表在1.5ml离心管中配制连接反应体系：
[0215]
表13测序接头连接反应体系
[0216][0217]
反应程序如下：20℃,15min
→
4℃,hold，关闭热盖，体积100μl。
[0218]
反应完后使用agencourt ampure xp beads纯化，最后用20ul eb洗脱dna并收集上清液(15ul)。
[0219]
(3)将测序接头连接产物进行pcr扩增，其中pcr循环数为7轮，获得测序文库，反应体系如下：
[0220]
表14 indexing pcr反应体系
[0221][0222]
反应程序如下：98℃,30s
→
7个循环【98℃,10s
→
65℃,75s】
→
65℃,5min
→
4℃,hold，热盖105℃，体积50μl。
[0223]
反应完后使用agencourt ampure xp beads纯化，最后用35ul eb洗脱dna并收集上清液(30ul)。
[0224]
定量：制备得到的文库用2.0(invitrogen)进行定量，-20℃保存。
[0225]
2100质检：用于2100检测的文库浓度范围：100pg/ul-10ng/ul，按说明书《5.2100_highsensitivitydna_qsg》进行操作。
[0226]
(4)二代高通量测序
[0227]
将不同样本制备所得文库按照相同的摩尔数混合，然后进行二代高通量测序，包括以下步骤：
[0228]
1)样品pooling
[0229]
用qubit测量文库的质量浓度(ng/ul)，通过2100确定文库的平均长度(bp)，根据以下公式进行质量浓度和摩尔浓度的换算。混合文库上机，现将文库按照上机数据量比例混合，即pooling，再计算混合文库的理论摩尔浓度。
[0230][0231]
注：yg001-047文库平均长度为350bp。
[0232]
2)将样品使用illumina miseq pe-300程序进行单末端测序；miseq产出的测序结果是fastq形式的dna序列，通过测序文库标签(index)、样本标签序列(barcode)将测序序列对应到每个样本，然后计算每个样本所测片段每个cg位点的甲基化状态。
[0233]
检测结果如下：
[0234]
(1)本发明检测方法比亚硫酸氢盐修饰后测序法具有更低的loq-hypothesis：本发明检测方法可以检测到0.05％甲基化(平均测序深度仅需100-200x，随浓度上升检测信号呈线性趋势)，而亚硫酸氢盐修饰后测序法检测限在0.5％甲基化之上(测序深度虽然上万x，但仍难以有线性关系出现)。以seq id no.5所示序列的两个cpg site的检测结果为例，示例位点信息如图2所示，本发明方法检测两个cpg site的loq情况如图3所示，测序深度平均每个梯度在100-200x；亚硫酸氢盐修饰后测序法检测相同两个cpg site的loq情况如图4所示，测序深度平均每个梯度大于9500x。
[0235]
(2)本发明检测方法与亚硫酸氢盐修饰后测序法相比具有更低的lod-hypothesis：以seq id no.7所示序列的ch1_63795447(-)和chr1_6379497(-)两个cpg位点为例，在2copies这个lod梯度，本发明检测方法在上千x测序深度可以稳定检出两个位点的甲基化情况，但是亚硫酸氢盐修饰后测序法即使是在～40000x的深度，仍然无法检出(检测出的c在测序噪音的区间)，具体如表15所示：
[0236]
表15本发明检测方法与bsp的lod-hypothesis
[0237][0238]
实施例4
[0239]
本实施例分别利用本发明所述检测方法和亚硫酸氢盐修饰后测序法对3例临床肺癌组织样本及3例对照组织样本的seq id no.1～seq id no.8所示8个序列的甲基化水平进行检测，具体检测步骤如下：
[0240]
s1、提取临床肺癌组织样本及对照组织样本基因组dna，对提取的dna使用重亚硫酸盐进行处理，然后对处理好的dna进行依赖甲基化多重pcr扩增：
[0241]
利用qiagen-qiaamp-dna-ffpe组织试剂盒(qiagen，cat#56404)提取样本基因组
dna，提取步骤根据试剂盒说明书进行。提取得到的dna按实施例3所述方法进行重亚硫酸盐处理以及依赖甲基化多重pcr扩增。
[0242]
步骤s2和s3同实施例3所述。
[0243]
检测结果表明，本发明检测方法可以在上千x的测序深度检测到肺癌与正常标本中的甲基化显著差异，而亚硫酸氢盐修饰后测序法则不可以。以seq id no.5所示序列的两个cpg site(cpg位点1ch5_40681550、cpg位点2ch5_40681569)的检测结果为例(所述两个cpg已知在肺癌组织中甲基化上调)，cpg位点1ch5_40681550的检测结果如图5所示，cpg位点2ch5_40681569的检测结果如图6所示。由图5和图6可知，使用本发明方法，在2000x的测序深度，即可发现检测位点在肺癌组织中标本中有明显的差别，而使用bsp的方法，在10倍的测序深度，即20000x深度，也无法有效检测到正常标本与肺癌标本的差别。
[0244]
对比例1
[0245]
本对比例一种dna甲基化检测方法，除了步骤s1中使用q5酶进行依赖甲基化多重pcr扩增外，其他与实施例3中本发明dna甲基化检测方法相同，其中，q5酶的用量与phusion dna聚合酶的用量相同。本发明检测方法在依赖甲基化多重pcr扩增时使用高保真酶phusion dna聚合酶(例如phusion hot start ii dna polymerase，phusion u hot start dna polymerase)、q5酶(例如hot start high-fidelity dna polymerase,hot start high-fidelity dna polymerase)、hieff酶(例如hieffhg热启动多重pcr酶)、kapa dna聚合酶(例如kapa hifi uracil kit，kapa2g快速热启动dna聚合酶)、pfu酶或superfi酶，均能取得较好的多重扩增效果；其中，phusion dna聚合酶的多重扩增效果最好，本对比例以phusion dna聚合酶(优选phusion u hot start dna polymerase)和q5酶(优选hot start high-fidelity dna polymerase)为例进行多重扩增效果的比较。
[0246]
分别利用本对比例dna甲基化检测方法和实施例3本发明dna甲基化检测方法中的依赖甲基化多重pcr法对10％甲基化cfdna mock标准品中50个不同的基因进行依赖甲基化多重pcr扩增，pcr测试统一条件为：98℃ 30s；18个循环(98℃ 15s，60℃ 2min，72℃ 1min)；72℃ 15min。
[0247]
通过2100分析，根据片段的长度范围来计算二聚体的比例，并比较两种扩增方法获得的扩增产物中二聚体的比例。
[0248]
结果如图7所示，结果表明，与使用q5酶进行依赖甲基化多重pcr扩增相比，使用phusion dna聚合酶(实施例3所述本发明检测方法中的依赖甲基化多重pcr法)进行依赖甲基化多重pcr扩增能有效降低扩增产物中二聚体的比例，使多重pcr扩增产物更利于后续的检测。
[0249]
对比例2
[0250]
本对比例一种dna甲基化检测方法，除了步骤s1中使用传统pcr法进行多重pcr扩增外，其他与实施例3中本发明dna甲基化检测方法相同，所述传统pcr法的反应程序如表16所示：
[0251]
表16传统pcr法反应程序
solution加入pcr反应产物并颠倒均匀
→
转移到dna recovery column，室温放置2min，8,000rpm离心1min
→
倒掉收集管中废液
→
重复一次
→
加入wash solution洗脱一次，12,000rpm离心1min
→
加入50ul的elution buffer，放置2min后12,000rpm离心1min。
[0261]
分别利用本对比例dna甲基化检测方法和实施例3本发明dna甲基化检测方法中的纯化方法对10％甲基化cfdna mock标准品中50个不同的基因的依赖甲基化多重pcr扩增产物进行纯化，比较两种纯化方法获得的扩增产物中二聚体的比例以及目标测序量的比例。
[0262]
结果如图10(图中column代表纯化柱纯化组)所示，结果表明，与使用纯化柱进行依赖甲基化多重pcr扩增产物的纯化相比，使用xp beads(本发明检测方法中的纯化方法)进行依赖甲基化多重pcr扩增产物的纯化能有效降低扩增产物中二聚体的比例，使纯化后的多重pcr扩增产物更利于后续的测序检测，提高目标测序量的比例。
[0263]
对比例5
[0264]
本对比例一种dna甲基化检测方法，除了步骤s1中依赖甲基化多重pcr扩增的循环数不同外，其他与实施例3中本发明dna甲基化检测方法相同，本对比例进行了2种不同条件的多重pcr扩增，分别为：(1)98℃ 30s；10个循环(98℃ 15s，65℃(每个循环下降0.5℃)15s，72℃ 15s)；17个循环(98℃ 15s，60℃ 15s，72℃ 15s)；72℃ 15min；(2)98℃ 30s；10个循环(98℃ 15s，65℃(每个循环下降0.5℃)15s，72℃ 15s)；20个循环(98℃ 15s，60℃ 15s，72℃ 15s)；72℃ 15min。本发明检测方法在依赖甲基化多重pcr扩增时扩增循环数为15～20时均能实现检测，尤其是当扩增循环数为15时能成功检测到目标基因数目是最多。
[0265]
分别利用本对比例2种dna甲基化检测方法和实施例3中本发明dna甲基化检测方法中的依赖甲基化多重pcr法对不同甲基化cfdna mock标准品(0.5％、1％、5％和10％)中的153个不同基因进行依赖甲基化多重pcr及后续测序分析比较，结果如图11所示，15个循环数在0.5％、1％、5％和10％的甲基化浓度条件下，能成功检测到目标基因数目明显多余17和20个循环数，其没检测到的目标基因数目是最少的。
[0266]
实施例5
[0267]
本实施例对本发明所述dna甲基化检测方法(medas)在不同甲基化程度的人类细胞系基因组dna(genomic dna)的检测效果进行评估，具体实施步骤如下：
[0268]
(1)制备不同甲基化梯度的基因组dna标品：0％和100％的人类基因组dna标品从zymo公司购买(cat#d5014),其中0％的标品源为hct116[dnmt1(-/-)dnmt3b(-/-)]细胞系；100％的标品源为0％的标品进行相关甲基化酶处理获得且经过相关测序验证。25％和50％的标品由0％与100％的标品按比例混合而成；
[0269]
(2)本实施例针对基因组上32个随机挑选cpg位点进行依赖甲基化多重pcr引物设计，详见下表：
[0270]
[0271][0272]
(3)之后，对0％25％50％和100％的标品按实施例3所述方法进行重亚硫酸盐处理以及依赖甲基化多重pcr扩增，步骤s2和s3同实施例3所述。
[0273]
检测结果表明，本发明检测方法对包含cg16673106，cg25381667，cg24016939，cg22101924，cg21715963，cg16712637，cg15811719，cg14603466，cg14589148，cg13119884，cg12622139，cg12180984，cg07696033，cg06080005，cg04234680，cg03556653，cg02596331 17个靶点在内的不同甲基化浓度的标品上的信号进行很好的线性拟合(r^2》0.8)，如图12所示。该结果显示，本发明所述检测方法可以在细胞系基因组dna上进行相关靶点甲基化程度的高通量、有效检测。
[0274]
实施例6
[0275]
本实施例对本发明所述dna甲基化检测方法(medas)在胃癌患者的血液标本中标志物进行检测评估，用于胃癌的早筛早诊检测，具体实施步骤如下：
[0276]
1.全血处理：本实施例共对152例胃镜无胃癌的正常人以及109例胃癌病人的血浆样本进行检测
[0277]
1.1利用edtak2抗凝真空采血管(bd，cat#367525)采集10ml全血，充分混匀，避免出现溶血，并在4-6小时内对全血进行血浆分离处理。
[0278]
1.2全血于低速离心机进行4℃，1600g，15min离心处理，小心吸取上层血浆，避免吸取中间的白膜层，所得血浆再次于高速离心机进行4℃，16000g，10min离心处理，获得所
需样品血浆。
[0279]
2.对血浆进行cfdna的提取
[0280]
具体方法：血浆dna提取具体操作步骤按照thermo fisher公司的magmax
tm cell-free dna isolation kit说明书进行。
[0281]
3.对提取的cfdna进行亚硫酸盐转化
[0282]
将提取的dna进行亚硫酸氢盐转化，使dna中未发生甲基化的胞嘧啶脱氨基转变成尿嘧啶，而甲基化的胞嘧啶保持不变，得到亚硫酸氢盐转化后的dna，转化具体操作按照zymo dna methylation-direct magprep的protocol进行，其中，投入的cfdna范围在5-20ng，本实施例优选在10ng。亚硫酸氢盐转化的产物全部用于进行多重甲基化扩增。
[0283]
4.对转化后的cfdna进行特定多个标志物的多重甲基化扩增
[0284]
转化后的产物全部进行多重甲基化扩增，其中的反应组分为通过组织筛选到的胃癌差异性103种标志物的引物组合，其中浓度在50-200nm，镁离子的浓度在2-5mm，本实施例优选在3mm，dntp mix浓度在100-600um，本实施例优选在200um，采用的酶为kapa2g fast multiplex pcr kit(roche,cat#kk5802)。具体的反应条件为：预变性，95℃ 5min,15-30个循环(变性，95℃ 15s,退火，58-66℃，本实施例优选63℃，4min。)，优选20个循环。
[0285]
多重反应体系配制如下：
[0286][0287]
5.对多重甲基化扩增产物进行文库构建
[0288]
利用nebultra
tm
ii dna library prep kit for与与multiplex oligos for(dual index primers set 1)对多重产物文库构建。相关操作见试剂盒说明书。
[0289]
6.按照等量的原则，对样本进行pooling处理，并上测序仪测序、进行相关统计分析。
[0290]
实验结果：选取全部的103个标志物利用随机森林模型进行建模分析，按照7：3的切分，进行100次重复，在99％的特异性下，测试集的1期的检测灵敏度为87.5％，2期的检测灵敏度为92.9％，3期的检测灵敏度为77.8％，4期的检测灵敏度为86.7％，对于所有的结直肠癌样本，检测的整体灵敏度为85.3％。整体的auc为0.974。表明了这些标志物在该方法下可用于胃癌的早期筛查。具体的roc curve请见图13。
[0291]
实施例7
[0292]
本实施例对本发明所述dna甲基化检测方法(medas)在乳腺癌患者的血液标本中标志物进行检测评估，用于乳腺癌的早筛早诊检测，具体实施步骤参照补充本发明实施例2，对30例正常人以及30例乳腺癌病人的血浆样本进行检测。通过组织中筛选到全部的54个
甲基化标志物利用随机森林模型进行建模分析，按照7：3的切分，进行100次重复，在99％的特异性下，测试集的11个1期样本的检测灵敏度为54.5％，15个2期样本的检测灵敏度为26.7％，4个3期样本的检测灵敏度为25％。对于所有的乳腺癌样本，检测的整体灵敏度为36.7％。整体的auc为0.948。表明了这些标志物在该方法下可用于乳腺癌的早期筛查。具体的roc curve请见图14。
[0293]
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对以上实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0294]
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。