CN104894233B - 一种多样本多片段dna甲基化高通量测序方法 - Google Patents
一种多样本多片段dna甲基化高通量测序方法 Download PDFInfo
- Publication number
- CN104894233B CN104894233B CN201510196548.2A CN201510196548A CN104894233B CN 104894233 B CN104894233 B CN 104894233B CN 201510196548 A CN201510196548 A CN 201510196548A CN 104894233 B CN104894233 B CN 104894233B
- Authority
- CN
- China
- Prior art keywords
- product
- pcr
- sequence
- dna
- multisample
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000007067 DNA methylation Effects 0.000 title claims abstract description 23
- 238000012163 sequencing technique Methods 0.000 claims abstract description 42
- 239000012634 fragment Substances 0.000 claims abstract description 38
- 230000011987 methylation Effects 0.000 claims abstract description 13
- 238000007069 methylation reaction Methods 0.000 claims abstract description 13
- LSNNMFCWUKXFEE-UHFFFAOYSA-M Bisulfite Chemical compound OS([O-])=O LSNNMFCWUKXFEE-UHFFFAOYSA-M 0.000 claims description 20
- 108090000623 proteins and genes Proteins 0.000 claims description 19
- 238000012408 PCR amplification Methods 0.000 claims description 16
- 230000008439 repair process Effects 0.000 claims description 12
- 238000011144 upstream manufacturing Methods 0.000 claims description 12
- 102000004190 Enzymes Human genes 0.000 claims description 10
- 108090000790 Enzymes Proteins 0.000 claims description 10
- 238000002156 mixing Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 claims description 8
- 102100024154 Cadherin-13 Human genes 0.000 claims description 7
- 101000762243 Homo sapiens Cadherin-13 Proteins 0.000 claims description 7
- 101000632056 Homo sapiens Septin-9 Proteins 0.000 claims description 7
- 102100028024 Septin-9 Human genes 0.000 claims description 7
- 238000002360 preparation method Methods 0.000 claims description 7
- 238000000246 agarose gel electrophoresis Methods 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 6
- 230000002068 genetic effect Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 102000012410 DNA Ligases Human genes 0.000 claims description 4
- 108010061982 DNA Ligases Proteins 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 4
- SUYVUBYJARFZHO-RRKCRQDMSA-N dATP Chemical compound C1=NC=2C(N)=NC=NC=2N1[C@H]1C[C@H](O)[C@@H](COP(O)(=O)OP(O)(=O)OP(O)(O)=O)O1 SUYVUBYJARFZHO-RRKCRQDMSA-N 0.000 claims description 4
- SUYVUBYJARFZHO-UHFFFAOYSA-N dATP Natural products C1=NC=2C(N)=NC=NC=2N1C1CC(O)C(COP(O)(=O)OP(O)(=O)OP(O)(O)=O)O1 SUYVUBYJARFZHO-UHFFFAOYSA-N 0.000 claims description 4
- 238000004064 recycling Methods 0.000 claims description 4
- 241000238631 Hexapoda Species 0.000 claims description 3
- 108010021757 Polynucleotide 5'-Hydroxyl-Kinase Proteins 0.000 claims description 3
- 102000008422 Polynucleotide 5'-hydroxyl-kinase Human genes 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 241001465754 Metazoa Species 0.000 claims description 2
- 241000196324 Embryophyta Species 0.000 claims 1
- 239000000047 product Substances 0.000 description 44
- 108020004414 DNA Proteins 0.000 description 43
- 239000011324 bead Substances 0.000 description 13
- 239000000872 buffer Substances 0.000 description 12
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 10
- 238000006243 chemical reaction Methods 0.000 description 8
- 238000005119 centrifugation Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 7
- 239000011534 wash buffer Substances 0.000 description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 7
- ISAKRJDGNUQOIC-UHFFFAOYSA-N Uracil Chemical compound O=C1C=CNC(=O)N1 ISAKRJDGNUQOIC-UHFFFAOYSA-N 0.000 description 6
- 230000003321 amplification Effects 0.000 description 6
- 238000003199 nucleic acid amplification method Methods 0.000 description 6
- 108091093088 Amplicon Proteins 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 235000019441 ethanol Nutrition 0.000 description 4
- 239000000706 filtrate Substances 0.000 description 4
- 238000000746 purification Methods 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 239000008223 sterile water Substances 0.000 description 4
- 102000008579 Transposases Human genes 0.000 description 3
- 108010020764 Transposases Proteins 0.000 description 3
- 239000003153 chemical reaction reagent Substances 0.000 description 3
- 238000010367 cloning Methods 0.000 description 3
- 229940035893 uracil Drugs 0.000 description 3
- 108010067770 Endopeptidase K Proteins 0.000 description 2
- LSNNMFCWUKXFEE-UHFFFAOYSA-N Sulfurous acid Chemical compound OS(O)=O LSNNMFCWUKXFEE-UHFFFAOYSA-N 0.000 description 2
- IQFYYKKMVGJFEH-XLPZGREQSA-N Thymidine Chemical compound O=C1NC(=O)C(C)=CN1[C@@H]1O[C@H](CO)[C@@H](O)C1 IQFYYKKMVGJFEH-XLPZGREQSA-N 0.000 description 2
- 238000003556 assay Methods 0.000 description 2
- 239000012148 binding buffer Substances 0.000 description 2
- 238000001369 bisulfite sequencing Methods 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 238000000749 co-immunoprecipitation Methods 0.000 description 2
- OPTASPLRGRRNAP-UHFFFAOYSA-N cytosine Chemical class NC=1C=CNC(=O)N=1 OPTASPLRGRRNAP-UHFFFAOYSA-N 0.000 description 2
- 239000012149 elution buffer Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002255 enzymatic effect Effects 0.000 description 2
- 239000000499 gel Substances 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 238000001114 immunoprecipitation Methods 0.000 description 2
- 230000002779 inactivation Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 239000006228 supernatant Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 108700028369 Alleles Proteins 0.000 description 1
- 241000894006 Bacteria Species 0.000 description 1
- DWRXFEITVBNRMK-UHFFFAOYSA-N Beta-D-1-Arabinofuranosylthymine Natural products O=C1NC(=O)C(C)=CN1C1C(O)C(O)C(CO)O1 DWRXFEITVBNRMK-UHFFFAOYSA-N 0.000 description 1
- 108091029430 CpG site Proteins 0.000 description 1
- 230000006429 DNA hypomethylation Effects 0.000 description 1
- KFZMGEQAYNKOFK-UHFFFAOYSA-N Isopropanol Chemical class CC(C)O KFZMGEQAYNKOFK-UHFFFAOYSA-N 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 108091028043 Nucleic acid sequence Proteins 0.000 description 1
- 108700020796 Oncogene Proteins 0.000 description 1
- 102000001253 Protein Kinase Human genes 0.000 description 1
- 108700025716 Tumor Suppressor Genes Proteins 0.000 description 1
- 102000044209 Tumor Suppressor Genes Human genes 0.000 description 1
- 239000002253 acid Substances 0.000 description 1
- 230000000692 anti-sense effect Effects 0.000 description 1
- 230000001580 bacterial effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- IQFYYKKMVGJFEH-UHFFFAOYSA-N beta-L-thymidine Natural products O=C1NC(=O)C(C)=CN1C1OC(CO)C(O)C1 IQFYYKKMVGJFEH-UHFFFAOYSA-N 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 210000000349 chromosome Anatomy 0.000 description 1
- 238000005352 clarification Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 239000013024 dilution buffer Substances 0.000 description 1
- XPPKVPWEQAFLFU-UHFFFAOYSA-N diphosphoric acid Chemical compound OP(O)(=O)OP(O)(O)=O XPPKVPWEQAFLFU-UHFFFAOYSA-N 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000001035 drying Methods 0.000 description 1
- 238000010828 elution Methods 0.000 description 1
- 230000006718 epigenetic regulation Effects 0.000 description 1
- 125000005909 ethyl alcohol group Chemical group 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012268 genome sequencing Methods 0.000 description 1
- 239000003292 glue Substances 0.000 description 1
- 229930182470 glycoside Natural products 0.000 description 1
- 150000002338 glycosides Chemical class 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000006607 hypermethylation Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 238000011534 incubation Methods 0.000 description 1
- 239000003112 inhibitor Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 210000004962 mammalian cell Anatomy 0.000 description 1
- 238000002493 microarray Methods 0.000 description 1
- 238000010208 microarray analysis Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000002018 overexpression Effects 0.000 description 1
- XEBWQGVWTUSTLN-UHFFFAOYSA-M phenylmercury acetate Chemical compound CC(=O)O[Hg]C1=CC=CC=C1 XEBWQGVWTUSTLN-UHFFFAOYSA-M 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 239000002244 precipitate Substances 0.000 description 1
- 108060006633 protein kinase Proteins 0.000 description 1
- 229940005657 pyrophosphoric acid Drugs 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 229940104230 thymidine Drugs 0.000 description 1
- 239000003643 water by type Substances 0.000 description 1
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明的一种多样本多片段DNA甲基化高通量测序方法，能够在一个文库中同时分析多个样本、多个目标片段的DNA甲基化状态，片段长度可以达到500bp左右，并且PCR重组发生率低，该方法还具有起始样本量少、针对性强、成本低、充分利用测序通量等优点。
Description
技术领域
本发明涉及高通量基因组甲基化DNA检测领域，尤其涉及一种多样本多片段DNA甲基化高通量测序方法。
背景技术
DNA甲基化是一种重要的表观遗传调控因子。在哺乳动物细胞中，DNA甲基化通常发生在CpG二核苷酸的胞嘧啶上，形成5甲基胞嘧啶。DNA甲基化的改变，尤其是抑癌基因启动子区的高甲基化在多种肿瘤的形成和发展过程当中起着重要的作用，同时DNA低甲基化也与癌基因的过表达和染色体的不稳定相关。
当前，研究DNA甲基化的方法有基因组整体甲基化检测的甲基化DNA免疫共沉淀芯片(microarray analysis following methylated DNA immunoprecipitation,MeMeDIP-chip)、甲基化DNA免疫共沉淀测序(sequencing analysis following methylated DNAimmunoprecipitation,MeDIP-seq)、简化的表观亚硫酸氢盐测序技术(reducedrepresentation bisulfite sequencing,RRBS)、重亚硫酸盐全基因组测序(bisulfitegenomic sequencing)和特异基因位点甲基化检测的传统的重亚硫酸盐克隆测序、焦磷酸测序。其中，单个基因位点DNA甲基化检测的金标准是重亚硫酸盐克隆测序，其利用重亚硫酸盐将基因组DNA中的未甲基化的胞嘧啶(C)修饰为尿嘧啶(U)，甲基化的胞嘧啶(C)保持不变，在PCR扩增中，尿嘧啶(U)变成胸腺嘧啶(T)，然后将PCR产物连接到特定载体上，转化到感受态细菌中，挑取10个以上单克隆菌株进行一代测序，根据参考序列中胞嘧啶位置C/T的频率来判断该位点的甲基化程度。
最近有文献报道了一种新的定量检测特异片段DNA甲基化的方法，这种方法结合了重亚硫酸盐转化、特异性片段扩增、依赖转座酶的建库方式和第二代测序，称为重亚硫酸盐扩增子测序(BiSulfite Amplicon Sequencing，BSAS)。这种方法可以对特异基因组区域内的甲基化胞嘧啶进行精准的绝对定量。但是该文章中使用的是依赖转座酶的建库方式，使用Illumina公司的Nextera XT建库试剂盒，使用双向标签序列Index，每个测序芯片上可以最多同时检测96个样本。
重亚硫酸盐克隆测序方法是一种可靠性和精确度很高的方法，能明确目的片段中每一个CpG位点的甲基化状态，但是不足之处是至少要测序10个以上克隆才能获得可靠数据，过程较为繁琐，耗时耗资过多。
BSAS的缺点是对于PCR产物使用转座酶酶切以后，就不能检测PCR产物两端被酶切掉序列的甲基化状态；每个样本需要构建一个文库。而且，还有文献报道扩增子测序中会出现等位基因间PCR重组的现象。
发明内容
针对上述重亚硫酸盐扩增子测序中存在的问题，本发明提供了一种针对多样本多片段的DNA甲基化高通量检测方法，该方法能够在一个文库中同时分析多个样本、多个目标片段的DNA甲基化状态，片段长度可以达到500bp左右，并且PCR重组发生率低，该方法还具有起始样本量少、针对性强、成本低、充分利用测序通量等优点。
本发明解决上述技术问题所采用的技术方案为：
一种多样本多片段DNA甲基化高通量测序方法，其特征在于包括以下步骤：步骤一：提取基因组DNA，对基因组DNA进行重亚硫酸盐处理，然后对处理好的基因组DNA进行PCR扩增，其中针对不同基因片段使用不同的PCR引物，相同的基因片段使用标签序列区分不同的样本；
步骤二：将步骤一中不同的PCR产物进行混合，使用混合好的PCR产物进行测序文库制备；
步骤三：对测序文库进行QPCR质检定量；
步骤四：上机测序，对不同样本特异性片段甲基化状况进行高精度分析。
为优化上述技术方案，具体措施还包括：
上述样品基因组DNA可以来源于任意物种，包括但不限于动物、植物、昆虫的基因组DNA。
上述步骤1中PCR引物针对不同基因、不同的扩增区域设计，片段长度可以达到500bp左右(最长片段长度大于500bp)；同时在上下游PCR引物的5’端分别加上3个碱基以上的标签序列，含不同标签序列的上下游引物分别组合用以区分不同的样本。
上述步骤一中进行PCR扩增所使用的聚合酶包括针对重亚硫酸盐转换后的DNA的热启动taq酶、常规的r-taq酶或其它聚合酶扩增，优选的是针对重亚硫酸盐转换后的DNA的热启动taq酶；还包括对扩增产物进行1.5％的琼脂糖胶电泳纯化的步骤。
还包括将纯化回收的产物定量，然后将不同样本的产物按照相同的摩尔数混合，得到混合的PCR产物的步骤。
上述对步骤二的文库制备利用如下步骤进行：
步骤1，对混合后的PCR产物进行定量，优选荧光定量的方法，取50ng进行文库制备；
步骤2，将混合的PCR产物进行3’末端修复和3’末端添加碱基A；优选通过T4聚核苷酸激酶、Klenow片段聚合酶、Klenow片段(3’-5’exo-)聚合酶、dNTP、dATP的作用同时进行3’末端修复和3’末端添加碱基A；
步骤3，,将3’末端添加碱基A的产物连接测序接头，获得测序接头连接产物；3’末端添加上碱基A的产物，通过优选的T4DNA连接酶的作用，与测序通用的接头进行连接，得到接头连接效率高的连接产物；
步骤4，将测序接头连接产物进行PCR扩增，其中PCR循环数为6轮，获得测序文库；
还包括对文库制备中各步骤的产物进行纯化回收的步骤，纯化方法包括但不限于磁珠纯化、纯化柱纯化、2％琼脂糖胶电泳纯化，优选磁珠纯化。
本发明还提供了根据上述的方法建立的测序文库。
本发明还提供了根据上述的方法检测多个样本、多个目标片段的DNA甲基化状态。
针对重亚硫酸盐扩增子测序中存在的问题，本发明提供了一种针对多样本多片段的DNA甲基化高通量检测方法，该方法能够在一个文库中同时分析多个样本、多个目标片段的DNA甲基化状态，片段长度可以达到500bp左右，并且PCR重组发生率低，该方法还具有起始样本量少、针对性强、成本低、充分利用测序通量等优点。
附图说明
图1显示了实施例中所有样本在CDH13片段内每个胞嘧啶位点甲基化的水平。
图2显示了实施例中所有样本在SEPT9片段内每个胞嘧啶位点甲基化的水平。
具体实施方式
本发明的主要技术流程：提取基因组DNA，然后对基因组DNA进行重亚硫酸盐处理，然后对处理好的基因组DNA进行PCR扩增，其中针对不同基因片段使用不同的PCR引物，相同的基因片段使用标签(barcode)序列区分不同的样本，然后将不同的PCR产物进行混合，混合好的PCR产物进行测序文库的制备，然后对文库进行QPCR质检定量，然后上机测序，最后进行不同样本特异性片段甲基化状况的高精度分析。
根据本发明的实施例，其中样品基因组DNA可以来源于任意物种，包括但不限于动物、植物、昆虫的基因组DNA。
根据本发明的实施例，将抽提的基因组DNA进行重亚硫酸盐处理，如利用EZ DNAMethylation-Gold KitTM(ZYMO)使非甲基化胞嘧啶转化为尿嘧啶。
根据本发明的实施例，其中PCR引物首先针对不同基因、不同的扩增区域设计，片段长度最长可达500bp左右；其次在上下游PCR引物的5’端分别加上3个碱基以上的标签(barcode)序列，含不同barcode的上下游引物分别组合用来区分不同的样本。
根据本发明的实施例，将区分样本的引物组合设计好以后，以重亚硫酸酸盐转化后的DNA为模板进行PCR扩增，其中所使用的聚合酶包括针对重亚硫酸盐转换后的DNA的热启动taq酶、常规的r-taq或其它聚合酶扩增，优选的是针对重亚硫酸盐转换后的DNA的热启动taq酶；对扩增产物进行1.5％的琼脂糖胶电泳纯化。
根据本发明的实施例，将纯化回收的产物使用NanoDrop2000(Thermo Fisher)定量，然后将不同样本的产物按照相同的摩尔数混合，得到混合的PCR产物。
根据本发明的实施例，其中文库制备的步骤如下：
步骤一对混合的PCR产物进行定量，并取50ng进行如下步骤建库；
对混合后的PCR产物进行定量，优选荧光定量的方法，如2.0荧光定量仪(Invitrogen)，取50ng进行如下所述步骤建库；
步骤二将混合的PCR产物进行3’末端修复和3’末端添加碱基A；
混合后的PCR产物，通过优选的T4聚核苷酸激酶、Klenow片段聚合酶、Klenow片段(3’-5’exo-)聚合酶、dNTP、dATP的作用同时进行3’末端修复和3’末端添加碱基A；
步骤三将3’末端添加碱基A的产物连接测序接头，以便获得测序接头连接产物；
3’末端添加上碱基A的产物，通过优选的T4DNA连接酶的作用，与测序通用的接头进行连接，以便得到接头连接效率高的连接产物；
步骤四将测序接头连接产物进行PCR扩增，其中PCR循环数为6轮，获得测序文库；
以接头连接后的产物为模板，使用Illumina公司通用的建库PCR引物，包括一条通用的上游引物，和一条带有标签(Index)序列的下游引物，使用高效的PCR扩增酶进行PCR，其中PCR循环数选择6轮。使用带有标签(Index)序列的引物进行PCR以后，可以将不同来源的文库进行混合，然后上机测序。
根据本发明的实施例，为了有效提高建库各步骤中的产物纯度、减少杂质干扰、有利于后续步骤的进行，可以对文库制备中各步骤的产物进行纯化回收，纯化方法包括但不限于磁珠纯化、纯化柱纯化、2％琼脂糖胶电泳纯化，优选磁珠纯化。
根据本发明的实施例，使用QPCR的方法对测序文库进行质检定量，其中使用Illumina phix control kit v3作为标准品。
根据本发明的实施例，通过二代测序平台进行测序，优选Illumina Miseq平台，Miseq的测序读长可以达到2*300，因此我们的目标片段长度可达500bp左右。
为了进一步说明本发明的内容，下面结合本发明的一个具体实施例作进一步详细描述：
1、提取基因组DNA
样本为临床全血样本，DNA抽提采用High Pure PCR Template Preparation Kit(Roche)，具体步骤如下：
a)第一次使用新的试剂盒时：
i.用4.5ml的无菌水完全溶解蛋白酶K粉末，使用完毕后放于-20℃储存；
ii.向装有Wash-buffer的瓶子中加入80ul的无水乙醇，并混匀；
b)将Elution Buffer于70度预热；
c)取1.5ml离心管，加入200ul血样，200ul Binding Buffer，40ul Proteinase K，立即70度温浴十分钟；
d)加入100ul异丙醇并混合完全；
e)将一个High Filter Tube插入一个收集管中，用移液管移取样品到过滤管的上层缓冲液池中，8000×g.离心1min；
f)从收集管中取下过滤管，弃滤液和收集管，将过滤管和一个新的收集管组装起来，加入500ul Inhibitor Removal Buffer到过滤管的上层缓冲液池中，8000×g离心1min；
g)从收集管中取下过滤管，弃滤液和收集管，将过滤管和一个新的收集管组装起来，加入500ul Wash Buffer到过滤管的上层缓冲液池中，8000×g.离心1min；
h)从收集管中取下过滤管，弃滤液和收集管，将过滤管和一个新的收集管组装起来加入500ul Wash Buffer到过滤管的上层缓冲液池中，8000×g.离心1min；
i)弃滤液，并额外全速离心整个High Filter Tube10s，弃收集管；
j)将过滤管插入一个新的无菌的1.5ml微型离心管，加入70μl预热的ElutionBuffer至过滤管的上层缓冲液池中，8000×g.离心1min，即得到基因组DNA。
2、重亚硫酸盐处理基因组DNA
采用EZDNA Methylation-Gold KitTM(ZYMO)试剂盒，具体步骤如下：
a)CT Conversion Reagent的制备：添加900ul水，50ul的M-Dissolving Buffer，和300ul的M-Dilution Buffer到一管CT Conversion Reagent中，在室温下溶解并且在摇床上摇动10min；
b)M-WASH BUFFER的制备：加入24ml无水乙醇至M-WASH BUFFER瓶中，并在瓶盖上做好标记；
c)在PCR管中添加：
①DNA量500ng，根据浓度算好体积，用水来补至20μl
②130ul的CT Conversion Reagent
轻弹试管或移液器操作来混合样品；
d)将样品管放到循环变温器并按以下步骤操作：
①98℃放置10min
②64℃放置2.5h
e)将柱子放入Collection Tube中,并加入600ul的M-Binding Buffer，并将步骤2的样品加入到柱子中，盖上盖子将柱颠倒数次来混合样品；
f)全速(>10,000x g)离心30sec，弃废液；
g)添加200ul的M-Wash Buffer到柱中，全速离心30sec；
h)添加200ul的M-Desphonation Buffer到柱中并且在室温(20℃-30℃)下放置15-20min；
i)全速离心30sec，弃废液；
j)加200ul的M-Wash Buffer到柱中，全速离心30sec；
k)再添加200ul的M-Wash Buffer并且离心30sec；
l)直接添加10ul的M-Elution Buffer到柱基质中，将柱放置在1.5ml的管中，全速离心来洗脱DNA。
3、针对不同基因片段和样本使用不同的PCR引物扩增
本实施例中，采用两个基因片段，基因名称分别为CDH13和SEPT9，均为人源的基因，引物序列如下表：
| 基因名称 | 上游引物序列 | 下游引物序列 |
| CDH13 | TTTGGGAAGTTGGTTGGTTG | ACAACCCCTCTTCCCTACCT |
| 基因名称 | 上游引物序列 | 下游引物序列 |
| SEPT9 | GTTATGTGGATTTTTTAAGTTAA | AAAAACAAAAAACTCCATCCTCC |
本实例中，针对108个样品检测上述两个基因片段的DNA甲基化，采用3个碱基的标签序列区分样本，因此上下游分别有11个标签序列，每个样本有唯一的上下游标签序列组合。所使用的标签序列如下表：
| 上游标签序列 | 下游标签序列 |
| ACA | AGT |
| ACG | AGA |
| TAC | TGC |
| TAG | TGA |
| CAT | CTG |
| ATC | ACT |
| ATG | AGC |
| TCA | TGT |
| TCG | CGT |
| CAG | CGC |
| CTA | CAC |
| ACA | AGT |
| ACG | AGA |
PCR扩增体系为：
PCR反应条件：
PCR扩增产物采用1.5％的琼脂糖胶电泳胶回收，胶回收试剂盒采用GeneJET GelExtraction Kit(Thermo)试剂盒，最后20ul洗脱产物。
4、PCR产物混合
取1ul凝胶回收的产物用NanoDrop2000(Thermo Fisher)定量，然后将不同样本的产物按照相同的摩尔数混合，得到混合的PCR产物。
5、文库制备
步骤一对混合后的PCR产物使用2.0(Invitrogen)进行定量，取50ng建库；然后对50ng待建库的DNA使用Agencourt AMPure XP beads纯化，具体步骤如下：
a)将AMPure XP beads置于室温下至少30min；
b)向1.5mL离心管中加入50ng待建库的DNA，再加入1.5倍体积混匀的AMPure XPbeads，在漩涡混匀器上混匀并室温放置5min；
c)将离心管放在磁力架上直至溶液澄清(大约3到5min)；
d)保持离心管在磁力架上，小心弃掉离心管中的上清，注意不要碰到beads；
e)继续保持离心管在磁力架上，加入500μL的70％乙醇；
f)将管子放置1min，使所有beads沉淀，弃掉乙醇；
g)重复步骤e)和f)一次；
h)在37℃加热模块中干燥样品5min或直到残余的乙醇完全消失；
i)加入16μL无菌水，在漩涡混匀器上混匀，室温温育2min；
j)保持离心管在磁力架上，静置2-3min，直到溶液澄清；
k)将16μL左右的上清液转移至一个新的1.5mL离心管中，弃掉beads。
步骤二将混合的PCR产物进行3’末端修复和3’末端添加碱基A；
将上一步得到的DNA按下表在1.5ml离心管中配制末端修复和加A反应体系：
| 样品DNA | 16ul |
| 10X End Repair Buffer | 2.5ul |
| 10mM ATP | 2.5ul |
| 10mM dNTP | 1ul |
| End Repair E1 | 1ul |
| End Repair E2 | 1ul |
| End Repair E3 | 1ul |
| 总体系 | 25ul |
然后将离心管放在25℃水浴锅中水浴20minutes，72℃失活20minutes。反应完后使用Agencourt AMPure XP beads纯化，最后用10.5ul无菌水洗脱DNA。
步骤三将3’末端添加碱基A的产物连接测序接头，以便获得测序接头连接产物；
将上一步得到的DNA按下表在1.5ml离心管中配制连接反应体系：
| 样品DNA | 10.5ul |
| Ligation Adaptor Mix | 1ul |
| 2X Quick Ligation Buffer | 12.5ul |
| Quick Ligation Enzyme Mix | 1ul |
| 总体系 | 25ul |
然后将离心管放在25℃水浴锅中水浴20minutes，72℃失活20minutes。反应完后使用Agencourt AMPure XP beads纯化，最后用20ul无菌水洗脱DNA。
备注：接头由一下序列退火形成：
Adaptor 1:5’-ACACTCTTTCCCTACACGACGCTCTTCCGATCT-3’
Adaptor2:5’-GATCGGAAGAGCACACGTCTGAACTCCAGTCAC-3’
步骤四将测序接头连接产物进行PCR扩增，其中PCR循环数为6轮，获得测序文库；
将上一步得到的DNA按下表在0.2ml离心管中配制PCR反应体系：
| 成分 | 体积(μL) |
| 样品DNA | 20 |
| TrueSeq Universal Primer | 1.25 |
| TrueSeq Primer-Index4 | 1.25 |
| 2X HiFi PCR Master Mix | 25 |
| Water | 2.5 |
| 总体系 | 50 |
PCR反应条件：
备注：PCR引物序列如下：
TrueSeq Universal Primer:
5’-AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT-3’
TrueSeq Primer-Index4:
5’-CAAGCAGAAGACGGCATACGAGATTGGTCAGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT-3’
反应完后使用Agencourt AMPure XP beads纯化，最后用10ul无菌水洗脱DNA。
6、文库质检定量
上步得到的文库用2.0(Invitrogen)进行定量，QPCR进行质检。
7、上机测序及数据分析
将样品使用Illumina MiseqPE-300程序进行双末端测序；Miseq产出的测序结果是fastq形式的DNA序列，通过测序文库标签(Index)、样本标签序列(barcode)将测序序列对应到每个样本，然后计算每个样本所测片段每个CG位点的甲基化状态。附图1显示了所有样本在CDH13片段内每个胞嘧啶位点甲基化的水平，附图2为所有样本在SEPT9片段内每个胞嘧啶位点甲基化的水平。
由上述的实施例可以看出，本发明的多样本多片段DNA甲基化高通量测序方法，能够在一个文库中同时分析多个样本、多个目标片段的DNA甲基化状态，片段长度可以达到500bp左右，并且PCR重组发生率低，该方法还具有起始样本量少、针对性强、成本低、充分利用测序通量等优点。同时需要指出的是，本实施例以人源基因为例进行本发明方法的阐述，但本发明的方法在人源基因上的应用仅仅局限于科学研究或者获取中间数据，并不用于人类疾病的诊断和治疗。
以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。
Claims (8)
1.一种多样本多片段DNA甲基化高通量测序方法，其特征在于由以下步骤组成：
步骤一：提取基因组DNA，对基因组DNA进行重亚硫酸盐处理，然后对处理好的基因组DNA进行PCR扩增，其中针对不同基因片段使用不同的PCR引物，相同的基因片段使用标签序列区分不同的样本，其中，不同基因是CDH13和SEPT9，CDH13的上游引物序列是TTTGGGAAGTTGGTTGGTTG，CDH13的下游引物序列是ACAACCCCTCTTCCCTACCT，SEPT9的上游引物序列是GTTATGTGGATTTTTTAAGTTAA，SEPT9的下游引物序列是AAAAACAAAAAACTCCATCCTCC，同时在上下游引物的5’端分别加上如下表所示的标签序列，含不同标签序列的上下游引物分别组合用以区分不同的样本
步骤二：将步骤一中不同的PCR产物进行混合，使用混合好的PCR产物进行测序文库制备，其中文库制备利用如下步骤进行：
步骤1，对混合后的PCR产物进行定量，进行文库制备；
步骤2，将混合的PCR产物进行3’末端修复和3’末端添加碱基A；
步骤3，将3’末端添加碱基A的产物连接测序接头，获得测序接头连接产物；
步骤4，将测序接头连接产物进行PCR扩增，其中PCR循环数为6轮，获得测序文库；
步骤三：对测序文库进行QPCR质检定量；
步骤四：上机测序，对不同样本特异性片段甲基化状况进行高精度分析。
2.根据权利要求1所述的一种多样本多片段DNA甲基化高通量测序方法，其中样品基因组DNA来源于动物、植物或昆虫的基因组DNA。
3.根据权利要求1所述的一种多样本多片段DNA甲基化高通量测序方法，其中步骤一中进行PCR扩增所使用的聚合酶包括针对重亚硫酸盐转换后的DNA的热启动taq酶、常规的r-taq酶或其它聚合酶。
4.根据权利要求3所述的一种多样本多片段DNA甲基化高通量测序方法，其中步骤一中进行PCR扩增所使用的聚合酶是针对重亚硫酸盐转换后的DNA的热启动taq酶。
5.根据权利要求1所述的一种多样本多片段DNA甲基化高通量测序方法，其中步骤一还包括对扩增产物进行1.5％的琼脂糖胶电泳纯化的步骤。
6.根据权利要求5所述的一种多样本多片段DNA甲基化高通量测序方法，其中还包括将纯化回收的产物定量，然后将不同样本的产物按照相同的摩尔数混合，得到混合的PCR产物的步骤。
7.根据权利要求1所述的一种多样本多片段DNA甲基化高通量测序方法，其中对步骤二的文库制备利用如下步骤进行：
步骤1，对混合后的PCR产物进行荧光定量，进行文库制备；
步骤2，将混合的PCR产物进行3’末端修复和3’末端添加碱基A是通过T4聚核苷酸激酶、Klenow片段聚合酶、dNTP、dATP的作用同时进行3’末端修复和3’末端添加碱基A；
步骤3，3’末端添加上碱基A的产物通过T4DNA连接酶的作用，与测序通用的接头进行连接，得到接头连接效率高的连接产物；
步骤4，将测序接头连接产物进行PCR扩增，其中PCR循环数为6轮，获得测序文库。
8.根据权利要求1所述的一种多样本多片段DNA甲基化高通量测序方法，其中对步骤二的文库制备利用如下步骤进行：
步骤1，对混合后的PCR产物进行荧光定量，进行文库制备；
步骤2，将混合的PCR产物进行3’末端修复和3’末端添加碱基A是通过T4聚核苷酸激酶、Klenow片段(3’-5’exo-)聚合酶、dNTP、dATP的作用同时进行3’末端修复和3’末端添加碱基A；
步骤3，3’末端添加上碱基A的产物通过T4DNA连接酶的作用，与测序通用的接头进行连接，得到接头连接效率高的连接产物；
步骤4，将测序接头连接产物进行PCR扩增，其中PCR循环数为6轮，获得测序文库。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510196548.2A CN104894233B (zh) | 2015-04-22 | 2015-04-22 | 一种多样本多片段dna甲基化高通量测序方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510196548.2A CN104894233B (zh) | 2015-04-22 | 2015-04-22 | 一种多样本多片段dna甲基化高通量测序方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104894233A CN104894233A (zh) | 2015-09-09 |
| CN104894233B true CN104894233B (zh) | 2018-04-27 |
Family
ID=54027185
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510196548.2A Active CN104894233B (zh) | 2015-04-22 | 2015-04-22 | 一种多样本多片段dna甲基化高通量测序方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104894233B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105274629B (zh) * | 2015-11-17 | 2017-11-17 | 安诺优达基因科技(北京)有限公司 | 简化的表观重亚硫酸氢盐测序用文库的构建方法及试剂盒 |
| CN108588176A (zh) * | 2018-05-06 | 2018-09-28 | 湖南大地同年生物科技有限公司 | 一种基于单链DNA构建R-loop高通量测序文库的方法 |
| WO2020135347A1 (zh) * | 2018-12-29 | 2020-07-02 | 深圳华大生命科学研究院 | 一种dna甲基化检测的方法、试剂盒、装置和应用 |
| CN111041069B (zh) * | 2019-12-26 | 2021-01-19 | 人和未来生物科技（长沙）有限公司 | 一种低起始量dna样本的高通量测序文库构建方法及其应用 |
| CN113088562A (zh) * | 2020-01-08 | 2021-07-09 | 深圳华大基因股份有限公司 | 一种新型的低起始量dna甲基化建库方法 |
| CN113930487B (zh) * | 2020-06-29 | 2023-03-17 | 广州市基准医疗有限责任公司 | 一种新型多样本多片段dna甲基化检测方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101824478A (zh) * | 2010-04-27 | 2010-09-08 | 武汉大学 | ChREBP基因甲基化状态的检测 |
| CN102796808A (zh) * | 2011-05-23 | 2012-11-28 | 深圳华大基因科技有限公司 | 甲基化高通量检测方法 |
| CN104450872A (zh) * | 2013-09-25 | 2015-03-25 | 上海市肿瘤研究所 | 一种高通量多样本多靶点单碱基分辨率的甲基化水平检测方法 |
-
2015
- 2015-04-22 CN CN201510196548.2A patent/CN104894233B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101824478A (zh) * | 2010-04-27 | 2010-09-08 | 武汉大学 | ChREBP基因甲基化状态的检测 |
| CN102796808A (zh) * | 2011-05-23 | 2012-11-28 | 深圳华大基因科技有限公司 | 甲基化高通量检测方法 |
| CN104450872A (zh) * | 2013-09-25 | 2015-03-25 | 上海市肿瘤研究所 | 一种高通量多样本多靶点单碱基分辨率的甲基化水平检测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104894233A (zh) | 2015-09-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104894233B (zh) | 一种多样本多片段dna甲基化高通量测序方法 | |
| CN103103624B (zh) | 高通量测序文库的构建方法及其应用 | |
| US20200232016A1 (en) | Methods and Compositions for the Analysis of Biological Molecules | |
| JP2022025140A (ja) | 標的化ゲノム解析のための方法 | |
| EP3172321B1 (en) | Polynucleotide enrichment using crispr-cas systems | |
| CN108300716B (zh) | 接头元件、其应用和基于不对称多重pcr进行靶向测序文库构建的方法 | |
| CN106715713B (zh) | 试剂盒及其在核酸测序中的用途 | |
| US20180195118A1 (en) | Systems and methods for detection of genomic copy number changes | |
| EP3612641A1 (en) | Compositions and methods for library construction and sequence analysis | |
| CN101586150A (zh) | 检测探针、通用寡核苷酸芯片及核酸检测方法及其用途 | |
| EP2691544B1 (en) | Method for verifying bioassay samples | |
| CN101921840A (zh) | 一种基于dna分子标签技术和dna不完全打断策略的pcr测序方法 | |
| CN108611398A (zh) | 通过新一代测序进行基因分型 | |
| CN103806111A (zh) | 高通量测序文库的构建方法及其应用 | |
| CN105695577B (zh) | 微量DNA中甲基化CpG岛高通量测序方法 | |
| CN109797436A (zh) | 一种测序文库构建试剂盒及其使用方法和应用 | |
| CN108138365A (zh) | 一种高通量的单细胞转录组建库方法 | |
| CN107760772A (zh) | 用于核酸配对末端测序的方法、组合物、系统、仪器和试剂盒 | |
| CN113166809B (zh) | 一种dna甲基化检测的方法、试剂盒、装置和应用 | |
| CN103555846A (zh) | 一种基于液相分子杂交原理的高通量低成本snp分型方法 | |
| CN115715323A (zh) | 一种高兼容性的PCR-free建库和测序方法 | |
| CN109112217A (zh) | 一种与猪体长和乳头数显著关联的遗传标记及应用 | |
| CN104450872A (zh) | 一种高通量多样本多靶点单碱基分辨率的甲基化水平检测方法 | |
| CN111074353A (zh) | 全基因组甲基化文库单链建库方法和得到的全基因组甲基化文库 | |
| CN102559856B (zh) | 去除测序文库中的载体片段的方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Chen Jing Inventor after: Wang Xilu Inventor before: Chen Jing Inventor before: Wu Qihan Inventor before: Dou Tonghai Inventor before: Wang Xilu |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20181101 Address after: 201100 D405 fourth, level 245, new Chun Ring Road, Minhang District, Shanghai. Co-patentee after: Jiangsu ang Pu Medical Technology Co., Ltd. Patentee after: Ang Piao bio tech ltd, Shanghai Address before: 201100 D405 4, level 245, new Chun Ring Road, Minhang District, Shanghai. Patentee before: Ang Piao bio tech ltd, Shanghai |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200226 Address after: 201100 D405 room, 4 floor, No. 245, Chun Chun Road, Shanghai, Minhang District Co-patentee after: Jiangsu ang Pu Medical Technology Co., Ltd. Patentee after: Ang Piao bio tech ltd, Shanghai Co-patentee after: Shanghai Petten Biotechnology Co., Ltd Address before: 201100 D405 room, fourth floor, No. 245, Chun Chun Road, Shanghai, Minhang District Co-patentee before: Jiangsu ang Pu Medical Technology Co., Ltd. Patentee before: Ang Piao bio tech ltd, Shanghai |