CN103806111A - 高通量测序文库的构建方法及其应用 - Google Patents
高通量测序文库的构建方法及其应用 Download PDFInfo
- Publication number
- CN103806111A CN103806111A CN201210460712.2A CN201210460712A CN103806111A CN 103806111 A CN103806111 A CN 103806111A CN 201210460712 A CN201210460712 A CN 201210460712A CN 103806111 A CN103806111 A CN 103806111A
- Authority
- CN
- China
- Prior art keywords
- region
- optionally
- specific
- genome
- methylation
- 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.)
- Pending
Links
- 238000012163 sequencing technique Methods 0.000 title claims abstract description 98
- 238000010276 construction Methods 0.000 title abstract description 6
- 239000000523 sample Substances 0.000 claims abstract description 179
- 230000011987 methylation Effects 0.000 claims abstract description 121
- 238000007069 methylation reaction Methods 0.000 claims abstract description 121
- 108020004414 DNA Proteins 0.000 claims abstract description 91
- 239000012634 fragment Substances 0.000 claims abstract description 85
- 238000000034 method Methods 0.000 claims abstract description 51
- LSNNMFCWUKXFEE-UHFFFAOYSA-M Bisulfite Chemical compound OS([O-])=O LSNNMFCWUKXFEE-UHFFFAOYSA-M 0.000 claims abstract description 19
- 230000003321 amplification Effects 0.000 claims abstract description 14
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 14
- ISAKRJDGNUQOIC-UHFFFAOYSA-N Uracil Chemical compound O=C1C=CNC(=O)N1 ISAKRJDGNUQOIC-UHFFFAOYSA-N 0.000 claims abstract description 6
- 229940035893 uracil Drugs 0.000 claims abstract description 3
- 108090000623 proteins and genes Proteins 0.000 claims description 75
- 238000012165 high-throughput sequencing Methods 0.000 claims description 73
- 108700039691 Genetic Promoter Regions Proteins 0.000 claims description 29
- 108091029523 CpG island Proteins 0.000 claims description 24
- 238000002360 preparation method Methods 0.000 claims description 19
- 239000011324 bead Substances 0.000 claims description 18
- 238000000746 purification Methods 0.000 claims description 17
- 108091026890 Coding region Proteins 0.000 claims description 15
- 238000007405 data analysis Methods 0.000 claims description 15
- 230000008439 repair process Effects 0.000 claims description 12
- 230000000903 blocking effect Effects 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 10
- UCSJYZPVAKXKNQ-HZYVHMACSA-N streptomycin Chemical compound CN[C@H]1[C@H](O)[C@@H](O)[C@H](CO)O[C@H]1O[C@@H]1[C@](C=O)(O)[C@H](C)O[C@H]1O[C@@H]1[C@@H](NC(N)=N)[C@H](O)[C@@H](NC(N)=N)[C@H](O)[C@H]1O UCSJYZPVAKXKNQ-HZYVHMACSA-N 0.000 claims description 10
- 238000012408 PCR amplification Methods 0.000 claims description 8
- 238000002156 mixing Methods 0.000 claims description 8
- 238000011144 upstream manufacturing Methods 0.000 claims description 7
- 108700009124 Transcription Initiation Site Proteins 0.000 claims description 6
- 229960005322 streptomycin Drugs 0.000 claims description 5
- 241000124008 Mammalia Species 0.000 claims description 4
- 108010006785 Taq Polymerase Proteins 0.000 claims description 4
- 238000000246 agarose gel electrophoresis Methods 0.000 claims description 4
- 239000010931 gold Substances 0.000 claims description 4
- 229910052737 gold Inorganic materials 0.000 claims description 4
- 108010017826 DNA Polymerase I Proteins 0.000 claims description 3
- 102000004594 DNA Polymerase I Human genes 0.000 claims description 3
- 230000001035 methylating effect Effects 0.000 claims description 3
- 210000003819 peripheral blood mononuclear cell Anatomy 0.000 claims description 3
- 102000012410 DNA Ligases Human genes 0.000 claims description 2
- 108010061982 DNA Ligases Proteins 0.000 claims description 2
- 108010021757 Polynucleotide 5'-Hydroxyl-Kinase Proteins 0.000 claims description 2
- 102000008422 Polynucleotide 5'-hydroxyl-kinase Human genes 0.000 claims description 2
- 210000004369 blood Anatomy 0.000 claims description 2
- 239000008280 blood Substances 0.000 claims description 2
- 244000005700 microbiome Species 0.000 claims description 2
- 238000000926 separation method Methods 0.000 claims description 2
- 108010068698 spleen exonuclease Proteins 0.000 claims description 2
- 239000003054 catalyst Substances 0.000 claims 1
- 102000053602 DNA Human genes 0.000 abstract description 6
- OPTASPLRGRRNAP-UHFFFAOYSA-N cytosine Chemical class NC=1C=CNC(=O)N=1 OPTASPLRGRRNAP-UHFFFAOYSA-N 0.000 abstract description 2
- 238000003752 polymerase chain reaction Methods 0.000 abstract 2
- 239000000047 product Substances 0.000 description 41
- 238000006243 chemical reaction Methods 0.000 description 24
- 238000009396 hybridization Methods 0.000 description 21
- 238000001514 detection method Methods 0.000 description 19
- 239000011534 wash buffer Substances 0.000 description 19
- 238000009826 distribution Methods 0.000 description 16
- 108010077544 Chromatin Proteins 0.000 description 15
- 210000003483 chromatin Anatomy 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 14
- 239000007788 liquid Substances 0.000 description 14
- 230000001965 increasing effect Effects 0.000 description 12
- 239000006228 supernatant Substances 0.000 description 12
- 238000013461 design Methods 0.000 description 10
- 239000000203 mixture Substances 0.000 description 10
- 239000012149 elution buffer Substances 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 239000000243 solution Substances 0.000 description 7
- 239000003153 chemical reaction reagent Substances 0.000 description 6
- 230000000295 complement effect Effects 0.000 description 6
- 108091081062 Repeated sequence (DNA) Proteins 0.000 description 5
- 230000008901 benefit Effects 0.000 description 5
- 239000000872 buffer Substances 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 108091029865 Exogenous DNA Proteins 0.000 description 4
- 210000000349 chromosome Anatomy 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000007899 nucleic acid hybridization Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000007067 DNA methylation Effects 0.000 description 3
- 108010090804 Streptavidin Proteins 0.000 description 3
- 239000012148 binding buffer Substances 0.000 description 3
- 210000004027 cell Anatomy 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000035772 mutation Effects 0.000 description 3
- 239000002773 nucleotide Substances 0.000 description 3
- 125000003729 nucleotide group Chemical group 0.000 description 3
- 230000003252 repetitive effect Effects 0.000 description 3
- 238000013518 transcription Methods 0.000 description 3
- 230000035897 transcription Effects 0.000 description 3
- YBJHBAHKTGYVGT-ZKWXMUAHSA-N (+)-Biotin Chemical compound N1C(=O)N[C@@H]2[C@H](CCCCC(=O)O)SC[C@@H]21 YBJHBAHKTGYVGT-ZKWXMUAHSA-N 0.000 description 2
- 230000030933 DNA methylation on cytosine Effects 0.000 description 2
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 2
- 238000011278 co-treatment Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 239000013024 dilution buffer Substances 0.000 description 2
- 238000004090 dissolution Methods 0.000 description 2
- 238000001962 electrophoresis Methods 0.000 description 2
- 238000013467 fragmentation Methods 0.000 description 2
- 238000006062 fragmentation reaction Methods 0.000 description 2
- 238000003780 insertion Methods 0.000 description 2
- 230000037431 insertion Effects 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 239000007791 liquid phase Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 108020004707 nucleic acids Proteins 0.000 description 2
- 102000039446 nucleic acids Human genes 0.000 description 2
- 150000007523 nucleic acids Chemical class 0.000 description 2
- 239000012264 purified product Substances 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 239000007790 solid phase Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000005406 washing Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 229920000936 Agarose Polymers 0.000 description 1
- 208000005623 Carcinogenesis Diseases 0.000 description 1
- 108020003215 DNA Probes Proteins 0.000 description 1
- 239000003298 DNA probe Substances 0.000 description 1
- 108010033040 Histones Proteins 0.000 description 1
- 108010010677 Phosphodiesterase I Proteins 0.000 description 1
- 238000001190 Q-PCR Methods 0.000 description 1
- 238000012181 QIAquick gel extraction kit Methods 0.000 description 1
- 108091061939 Selfish DNA Proteins 0.000 description 1
- 238000012300 Sequence Analysis Methods 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 229960002685 biotin Drugs 0.000 description 1
- 235000020958 biotin Nutrition 0.000 description 1
- 239000011616 biotin Substances 0.000 description 1
- 239000007853 buffer solution Substances 0.000 description 1
- 230000036952 cancer formation Effects 0.000 description 1
- 231100000504 carcinogenesis Toxicity 0.000 description 1
- 230000003915 cell function Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000005119 centrifugation Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005202 decontamination Methods 0.000 description 1
- 230000003588 decontaminative effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 239000005547 deoxyribonucleotide Substances 0.000 description 1
- 125000002637 deoxyribonucleotide group Chemical group 0.000 description 1
- 230000019975 dosage compensation by inactivation of X chromosome Effects 0.000 description 1
- 238000010828 elution Methods 0.000 description 1
- 230000013020 embryo development Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000001973 epigenetic effect Effects 0.000 description 1
- 230000007608 epigenetic mechanism Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011534 incubation Methods 0.000 description 1
- 210000004969 inflammatory cell Anatomy 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012164 methylation sequencing Methods 0.000 description 1
- 238000010369 molecular cloning Methods 0.000 description 1
- 230000004879 molecular function Effects 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 238000011176 pooling Methods 0.000 description 1
- 238000012802 pre-warming Methods 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 239000002096 quantum dot Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000008247 solid mixture Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了提供了高通量测序文库的构建方法及其应用。其中，构建高通量测序文库的方法包括：将基因组DNA片段化；将DNA片段进行末端修复；在经过末端修复的DNA片段的3’末端添加碱基A；将具有粘性末端A的DNA片段与甲基化接头相连；利用特异性探针对连接产物进行杂交捕获，以便获得目的片段；将该目的片段进行重亚硫酸盐处理，以便将该目的片段中非甲基化的胞嘧啶转换为尿嘧啶；将经过转换的目的片段进行PCR扩增；分离纯化扩增产物，该扩增产物构成高通量测序文库。采用本发明的构建高通量测序文库的方法及其应用，能够方便有效地构建样本的基因组特定区域的高通量测序文库。
Description
技术领域
本发明涉及生物技术领域。具体地，涉及DNA甲基化检测技术，特别是涉及基因组特定区域的甲基化检测。更具体地，本发明提供了一种构建高通量测序文库的方法、一种确定样本的基因组特定区域的甲基化信息的方法、一种用于确定样本的基因组特定区域的甲基化信息的装置以及一种用于构建样本的基因组特定区域高通量测序文库的试剂盒。
背景技术
DNA甲基化是研究最为深入的表观遗传学机制，DNA甲基化在维持正常细胞功能、抑制寄生DNA成分对基因组完整性的损害、染色质结构修饰、X染色体失活、基因组印记、胚胎发育以及人类肿瘤发生中起着重要作用，是目前新的研究热点之一。
然而，目前对基因组特定区域如启动子区域、CpG岛区域、CpG岛外区域以及印记基因区域的甲基化检测的研究，仍有待改进。
发明内容
本发明旨在解决现有技术问题的至少之一。由此，为了代表性检测基因组上特定区域的甲基化信息，本发明提供了高通量测序文库的构建方法及其应用。
根据本发明的一个方面，本发明提供了一种构建高通量测序文库的方法。根据本发明的实施例，该方法包括以下步骤：将基因组DNA片段化，以便获得DNA片段；将该DNA片段进行末端修复，以便获得经过末端修复的DNA片段；在该经过末端修复的DNA片段的3’末端添加碱基A，以便获得具有粘性末端A的DNA片段；将该具有粘性末端A的DNA片段与甲基化接头相连，以便获得连接产物；利用特异性探针对该连接产物进行杂交捕获，以便获得目的片段；将该目的片段进行PCR扩增，以便获得扩增产物；以及分离纯化所述扩增产物，该扩增产物构成该高通量测序文库。根据本发明的实施例，可以采用的特异性探针是对已知甲基化位点特异性的，例如，该特异性探针是基于采用人类基因组作为参考序列，并且采用已知具有甲基化位点的基因区域作为靶序列而设计的，其中，该已知具有甲基化位点的基因区域可以为下表I中所列出基因的至少一种的编码区和启动子区。
表I-1分子功能相关基因
表I-2生物学加工相关基因
表I-3细胞组分相关基因
表I-4染色质相关基因
利用根据本发明实施例的构建高通量测序文库的方法，能够有效地构建基因组DNA样品的高通量测序文库，特别是能够有效地构建基因组DNA样品的已知甲基化位点的特定区域的高通量测序文库，从而能够有效、充分地应用于高通量测序技术，通过对文库的测序，然后基于对测序结果的数据分析，能够有效地获得基因组特定区域的甲基化位点信息，实现对基因组DNA样品的基因组特定区域的甲基化检测。
根据本发明的另一方面，本发明提供了一种确定样本的基因组特定区域的甲基化信息的方法。根据本发明的实施例，该方法包括下列步骤：根据前面所述构建高通量测序文库的方法，构建该样本的基因组特定区域的高通量测序文库；对该样本的基因组特定区域的高通量测序文库进行测序，以便得到测序结果；以及对该测序结果进行数据分析，以便确定该样本的基因组特定区域的甲基化信息。
利用根据本发明实施例的确定样本的基因组特定区域的甲基化信息的方法，能够准确地确定样本的基因组特定区域的甲基化信息，从而实现对样本的基因组特定区域的甲基化检测。
根据本发明的再一方面，本发明提供了一种用于确定样本的基因组特定区域的甲基化信息的装置。根据本发明的实施例，该装置包括：文库制备单元，该文库制备单元用于制备样本的基因组特定区域的高通量测序文库，该文库制备单元内设置有特异性探针；测序单元，该测序单元与该文库制备单元相连，并且从该文库制备单元接收该样本的基因组特定区域的高通量测序文库，以便用于对该样本的基因组特定区域的高通量测序文库进行测序，获得测序结果；以及数据分析单元，该数据分析单元与该测序单元相连，并且从该测序单元接收该测序结果，以便对该测序结果进行数据分析，确定该样本的基因组特定区域的甲基化信息。
利用根据本发明实施例的用于确定样本的基因组特定区域的甲基化信息的装置，能够方便准确地确定样本的基因组特定区域的甲基化信息，可以应用于多种针对基因组特定区域的甲基化的研究。
根据本发明的又一方面，本发明提供了一种用于构建样本的基因组特定区域的高通量测序文库的试剂盒。根据本发明的实施例，该试剂盒包括：特异性探针，该特异性探针是对已知甲基化位点特异性的。利用根据本发明实施例的用于构建样本的基因组特定区域的高通量测序文库的试剂盒，能够方便有效地构建样本的基因组特定区域的高通量测序文库。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1：显示了根据本发明一个实施例的构建高通量测序文库的方法的流程示意图；
图2：显示了根据本发明一个实施例的方法确定基因组特定区域甲基化信息时，在不同覆盖深度下（覆盖深度≥1及覆盖深度＞5），每条染色质上的捕获区域占探针靶区域的百分比图；
图3：显示了根据本发明一个实施例的方法确定基因组特定区域甲基化信息时，在不同覆盖深度下，各条染色质中检测到甲基化信息的启动子占该染色质的总启动子的百分比图；
图4：显示了根据本发明一个实施例的方法确定基因组特定区域甲基化信息时，基因组上启动子区域、CpG岛、CpG岛外（在本文中指为CGI shore）及印记基因区域的甲基化水平分布结果，其中，
(a)显示了根据本发明一个实施例的确定的样本的基因组CpG岛、CGI shore区域的甲基化水平分布图，
(b)显示了根据本发明一个实施例的确定的样本的基因组启动子区域的甲基化水平分布图，
(c)显示了样本的基因组特定区域的原始分布和根据本发明一个实施例的确定的样本的基因组特定区域的高通量测序文库的reads分布及启动子、CpG岛区域的甲基化水平分布图；
图5：显示了根据本发明一个实施例的用于确定样本的基因组特定区域的甲基化信息的装置的示意图；
图6：显示了根据本发明一个实施例的测序读段（在本文中也成为“reads”）的插入片段长度分布；
图7：显示了根据本发明一个实施例的捕获区域的各碱基的测序深度和累积测序深度统计图；
图8：显示了根据本发明一个实施例的在读段比对过程中各位置碱基错配发生频数的示意图；
图9：显示了根据本发明一个实施例的炎黄细胞系捕获的各表观相关基因在不同深度下的覆盖度和累计覆盖度的示意图；以及
图10：显示了根据本发明一个实施例的个别捕获基因的测序覆盖的直观展示结果。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
构建高通量测序文库的方法
根据本发明的一个方面，本发明提供了一种构建高通量测序文库的方法。参考图1，根据本发明的实施例，该方法包括以下步骤：
首先，将基因组DNA片段化，以便获得DNA片段。在本发明中所使用的术语“DNA”可以是任何包含脱氧核糖核苷酸的聚合物，包括但不限于经过修饰的或者未经修饰的DNA。本领域的技术人员可以理解，基因组DNA的来源不受特别限制，可以从任何可能的途径获得，可以是通过市售直接获得，也可以是从其他实验室直接获取，还可以是直接从样本中提取。根据本发明的实施例，可以从样本中提取获得基因组DNA。根据本发明的一个实施例，构建高通量测序文库的方法可以进一步包括从样本中提取基因组DNA的步骤。根据本发明的一些具体示例，样本可以来源于哺乳动物、植物、和微生物的至少一种。根据本发明的一些实施例，哺乳动物可以为人和小鼠的至少一种。根据本发明的一个实施例，基因组DNA可以为人类全血基因组DNA，优选为外周血单核细胞基因组DNA。发明人发现，当采用YH cell基因组DNA构建高通量测序文库时，从样本中提取基因组DNA的操作方便易行，且获得的DNA质量好、甲基化信息完整，由其构建的样本的基因组特定区域的高通量测序文库能够方便地应用于高通量测序技术，从而基于对测序结果的数据分析就能方便有效地获得样本的基因组特定区域的甲基化信息。根据本发明的实施例，基因组DNA的量不受特别限制，根据本发明的具体示例，优选基因组DNA的量为2μg。发明人惊奇地发现，当基因组DNA的量为2μg时，根据本发明实施例的构建高通量测序文库的方法构建的样本的基因组特定区域的高通量测序文库，能够非常方便地应用于高通量测序技术，如Solexa测序技术，且文库测序结果准确，可重复性好，包含的特定区域的甲基化信息准确、甲基化位点覆盖率高。
其次，将DNA片段进行末端修复，以便获得经过末端修复的DNA片段。根据本发明的一个实施例，在将DNA片段进行末端修复前，可以进一步包括纯化DNA片段的步骤，由此，使得后续的末端修复易于进行。根据本发明的实施例，将DNA片段进行末端修复可以利用Klenow片段、T4DNA聚合酶和T4多核苷酸激酶进行，其中，所述Klenow片段具有5’→3’聚合酶活性和3’→5’聚合酶活性，但缺少5’→3’外切酶活性。由此，能够方便准确地对DNA片段进行末端修复。根据本发明的实施例，还可以进一步包括对经过末端修复的DNA片段进行纯化的步骤，由此能够方便地进行后续处理。
接下来，在经过末端修复的DNA片段的3’末端添加碱基A，以便获得具有粘性末端A的DNA片段。根据本发明的一个实施例，可以利用Klenow (3’-5’exo-)，即具有3’→5’外切酶活性的Klenow，在经过末端修复的DNA片段的3’末端添加碱基A。由此，能够方便准确地将碱基A添加到经过末端修复的DNA片段的3’末端。根据本发明的实施例，还可以进一步包括对具有粘性末端A的DNA片段进行纯化的步骤，由此能够方便地进行后续处理。
接着，将具有粘性末端A的DNA片段与甲基化接头相连，以便获得连接产物。本发明中所使用的术语“甲基化接头”是指这样的一种接头，在其核苷酸序列中，所有C位点均被甲基化修饰。根据本发明的一个实施例，在将具有粘性末端A的DNA片段与甲基化接头相连前，可以进一步包括对常规测序所使用的接头进行甲基化的步骤。由此，能够有效避免测序接头对后续重亚硫酸盐处理等操作带来的干扰，例如重亚硫酸盐处理处理过程中接头序列可能会被改变。本领域的技术人员可以理解，对接头进行甲基化的方法不受特别限制，可以利用本领域已知的任何方法对测序接头进行甲基化。
根据本发明的一些实施例，甲基化接头中还可以进一步包含标签，由此可以方便地同时构建多种样本的基因组特定区域的高通量测序文库，并能够有效地应用于高通量测序平台，从而在对测序结果进行数据分析后，基于标签的序列信息，就能够准确地区分多种样本的基因组特定区域的高通量测序文库的序列信息以及样本的基因组特定区域的甲基化信息，由此，能够充分地利用高通量测序平台，且能够节省时间、降低成本。
根据本发明的实施例，标签的长度为6bp，在文库制备时不同样品分别连接不同的标签接头，捕获前将多个不同文库混合到一起形成一个新文库，使用该新文库进行探针捕获并测序。测序所得数据是多个样品混合在一起的，此时根据每段reads（测序仪所读出的序列结果）中不同的标签序列将不同样品区分开。此方法可以大大降低成本、时间、人力。根据本发明的实施例，
根据本发明的一个实施例，将具有粘性末端A的DNA片段与甲基化接头相连是利用T4DNA连接酶进行的，由此可以方便地获得连接产物。根据本发明的实施例，还可以进一步包括对连接产物进行纯化的步骤，由此能够方便地进行后续处理。
然后，利用特异性探针对所述连接产物进行杂交捕获，以便获得目的片段。根据本发明的实施例，这里的术语“特异性探针”是指探针是对已知甲基化位点特异性的。根据本发明的具体示例，特异性探针是基于采用人类基因组作为参考序列，并且采用基因组上已知具有甲基化位点的特定基因区域作为靶序列而设计的，具体地，已知具有甲基化位点的基因区域包括选自启动子区域、CpG岛区域、CpG岛外区域以及印记基因区域的至少一种，由此，利用根据本发明实施例的特异性探针进行杂交捕获，能够有效地捕获样本中与靶序列互补的序列，即样本中已知具有甲基化位点的基因区域（在本说明书中，有时也称为“基因组特定区域”）。
根据本发明的实施例，可以用于设计特异性探针的已知具有甲基化位点的基因区域为表I中所列出基因的至少一种的编码区和启动子区。根据本发明的实施例，所述编码区为外显子区域序列，所述启动子区为基因转录起始位点的上游2200bp到下游500bp的区域。根据本发明的实施例，所述特异性探针是采用eArray系统设计的。根据本发明的实施例，任选地，所述探针的长度为12mer。
如前面表I所示的基因是发明人利用gene ontology等数据库资源，经过大量试验筛选获得的。发明人意外发现采用这些基因制备的探针能够最有效地捕获所期望的目的片段，并且有助于后续的研究。根据本发明的实施例，可以采用基因TSS（转录起始位点）上游2200bp到下游500bp作为promoter（启动子）区域，以及exon（外显子）区域序列作为基因的编码区，并针对这些区域的序列信息进行捕获探针的设计。另外，令人吃惊地，利用这些基因座位靶筛选得到的测序结果，对于各染色体的覆盖不存在偏差。
根据核酸的互补配对原则，单链状态的捕获探针可以与单链状态的目的序列互补结合，从而成功地将目标区域捕获。根据本发明的实施例，探针设计可选择固相捕获芯片（探针固定在固体载体上）或液相捕获探针（探针游离在液体中），然而固相捕获芯片因探针长度、探针密度、价格偏高等诸多因素限制，液相捕获即作为首选。
根据本发明的实施例，采用安捷伦公司（Agilent）的探针设计系统eArray设计探针，探针长度120mer，探针可覆盖长度范围大，从小于200kb到24Mb甚至更长。eArray探针设计系统可以方便地使用生物信息学工具window masker（窗口序列屏蔽）和repeat masker（重复序列屏蔽）对目标区域分析并进行屏蔽，由此，可以避免对这些区域进行探针设计，非常有效地减少实验中的捕获干扰以及后续序列分析时发生的比对干扰；并且缩短覆盖长度可以在一定程度上减少成本。
根据本发明的实施例，核酸中高CG（CG碱基含量高于60%）序列由于碱基C、G的分子结构不同会导致捕获效率低于常规（ATCG碱基平均含量各25%）序列。对于重点研究对象CGI（CpG Island，CpG岛）区域，可以通过提高探针设计量，更好更多地得到CGI数据。
此外，根据本发明的一个实施例，在杂交捕获前，可以进一步包括利用诸如cot-1DNA和接头封闭序列的单链寡核苷酸对连接产物（尤其是连接产物的基因组序列中的重复区域）和连接产物上的甲基化接头进行杂交封闭的步骤。发明人惊奇地发现，当使用cot-1DNA和接头封闭序列分别对连接产物（尤其是连接产物的基因组序列中的重复区域）和连接产物上的甲基化接头进行杂交封闭后，能够显著地增强对连接产物的杂交捕获。根据本发明的实施例，cot-1DNA的使用量不受特别限制，根据具体的示例，优选采用过量的cot-1DNA对连接产物的基因组序列中的重复区域进行杂交封闭。其中，这里所使用的术语“过量”是指cot-1DNA的量远大于待进行杂交捕获的连接产物的量，即采用cot-1DNA的量可以是待进行杂交捕获的连接产物的量的2倍以上。根据本发明的具体示例，优选，采用cot-1DNA的量为待进行杂交捕获的连接产物的量的5倍。根据本发明的一些实施例，采用cot-1DNA的量小于待进行杂交捕获的连接产物的量的5倍，则封闭杂交不彻底，重复序列的非特异性强杂交背景信号干扰强烈，严重影响核酸杂交的效率；而采用cot-1DNA的量大于待进行杂交捕获的连接产物的量的5倍，则过多的cot-1DNA会影响探针与连接产物的结合，同样会影响核酸杂交的效率。由此，采用待进行杂交捕获的连接产物的量的5倍的cot-1DNA对连接产物的基因组区域重复序列进行杂交封闭，能够方便、有效地进行封闭，以去掉重复序列DNA，从而在后续的核酸杂交过程中，能够有效避免重复序列产生的非特异性强杂交背景信号的干扰，显著提高核酸杂交的效率，增强杂交效果。根据本发明的实施例，接头封闭序列包括选自Block1和Block2的至少一种，由此，能够有效地对连接产物上的甲基化接头进行封闭。根据本发明的实施例，可以采用1μg的连接产物进行所述杂交捕获，由此能够提高杂交捕获的效率。根据本发明的具体示例，利用特异性探针对所述连接产物进行杂交捕获，可以进一步包括利用链霉素磁珠捕获目的片段，由此，能够高效地捕获目的片段。
然后，将目的片段进行PCR扩增，以便获得扩增产物。根据本发明的实施例，可以使用热启动taq DNA聚合酶对经过转换的目的片段进行PCR扩增。根据本发明的实施例，热启动taq DNA聚合酶的种类不受特别限制，根据本发明的具体示例，热启动taq DNA聚合酶可以为r-taq聚合酶，由此PCR扩增效率高、用时少。
最后，分离纯化扩增产物，所得到的扩增产物构成全基因组甲基化高通量测序文库。根据本发明的实施例，分离纯化扩增产物的方法不受特别限制，根据本发明的具体示例，可以通过选自磁珠纯化、纯化柱纯化和2%的琼脂糖凝胶电泳的至少一种进行，优选通过2%的琼脂糖凝胶电泳进行。根据本发明的一些具体示例，高通量测序文库的文库片段长度为300-450bp，由此，高通量测序文库能够方便有效地应用于高通量测序平台如Solexa测序平台，且可重复性好，测序结果真实可靠，包含特异性探针所针对的基因组特定区域的甲基化信息较完整。
根据本发明的实施例，在得到目的片段之后，可以将目的片段进行重亚硫酸盐处理，以便将目的片段中非甲基化的胞嘧啶转换为尿嘧啶，获得经过转换的目的片段。根据本发明的实施例，在将目的片段进行重亚硫酸盐处理之前，可以进一步包括将目的片段与片段化的λ-DNA混合。发明人发现，通过添加外源DNA（λ-DNA），即将目的片段与外源DNA混合，然后进行重亚硫酸盐高效共处理，对目标DNA片段能够起到保护作用，最大限度地降低重亚硫酸盐对微量DNA的破坏，可以进一步提高检测精度，使得较少量的基因组DNA，甚至纳克级，例如50-150ng基因组的甲基化检测成为现实。根据本发明的实施例，片段化的λ-DNA的添加量不受特别限制，根据具体的示例，优选片段化的λ-DNA的量为200-400ng，更优选为200ng。本领域技术人员能够理解，可以通过本领域已知的任意方法制备这些片段化的λ-DNA，例如可以随同前面的DNA片段化处理一起进行制备。
重亚硫酸盐处理可以通过本领域已知的任何方法进行，根据本发明的具体示例，可以采用商品化的试剂盒进行，优选地采用EZ DNA Methylation-Gold KitTM（ZYMO）进行。发明人惊奇地发现，采用EZ DNA Methylation-Gold KitTM（ZYMO）对目的片段进行重亚硫酸盐处理时，方便快捷，且处理效果好，目的片段中非甲基化的胞嘧啶能够高效准确地转换为尿嘧啶，并且利于后续处理。
由此，根据本发明的实施例，可以在捕获后直接测序，进而，基于测序结果可以进行基因的单核苷酸多态性（SNPs）、核苷酸突变、插入、缺失（indel）或拷贝数变化（CNVs）分析。也可以在捕获后经重亚硫酸盐处理及测序进行DNA甲基化状态分析，由此，可以进行如甲基化密度、不同元件甲基化水平、胞嘧啶甲基化分析和差异性甲基化区域（DMRs）分析等。
利用根据本发明实施例的构建高通量测序文库的方法，能够有效地构建样本的基因组特定区域的高通量测序文库，从而能够有效、充分地应用于高通量测序技术，通过对高通量测序文库的测序，然后基于对测序结果的数据分析，就能够有效地获得样本的基因组特定区域的甲基化信息，实现对样本的基因组特定区域的甲基化检测。
确定样本的基因组特定区域的甲基化信息的方法和装置
根据本发明的另一方面，本发明提供了一种确定样本的基因组特定区域的甲基化信息的方法。根据本发明的实施例，该方法包括下列步骤：根据本发明实施例的构建高通量测序文库的方法构建样本的基因组特定区域的高通量测序文库；对该样本的基因组特定区域的高通量测序文库进行测序，以便得到测序结果；以及对测序结果进行数据分析，以便确定样本的基因组特定区域的甲基化信息。
根据本发明的一些实施例，测序是利用高通量测序技术进行的。本领域的技术人员可以理解，可以通过本领域已知的任何高通量测序技术进行测序，根据本发明的具体示例，优选地利用Hiseq2000测序仪进行测序。发明人发现，利用Hiseq2000测序仪对样本的基因组特定区域的高通量测序文库进行测序，能够有效地获得测序结果，且测序用时少、效率高、测序结果准确，可重复性好。
利用根据本发明实施例的确定样本的基因组特定区域的甲基化信息的方法，能够有效地构建样本的基因组特定区域的高通量测序文库，并且能够通过高通量测序技术如Solexa测序技术实现对文库的准确测序，基于对测序结果的数据分析，就能够准确地确定样本的基因组特定区域的甲基化信息，从而实现对样本的基因组特定区域的甲基化检测，且特定区域的甲基化位点覆盖多，获得甲基化信息完整。
根据本发明的再一方面，本发明提供了一种用于确定样本的基因组特定区域的甲基化信息的装置。参考图5，根据本发明的一个实施例，该装置1000包括：文库制备单元100、测序单元200以及数据分析单元300。
根据本发明的实施例，文库制备单元100用于制备样本的基因组特定区域的高通量测序文库，其中，文库制备单元100内设置有特异性探针。根据本发明的实施例，特异性探针是对已知甲基化位点特异性的。根据本发明的具体示例，特异性探针是基于采用人类基因组作为参考序列，并且采用基因组上已知具有甲基化位点的特定基因区域作为靶序列而设计的，具体地，已知具有甲基化位点的基因区域包括选自启动子区域、CpG岛区域、CpG岛外区域以及印记基因区域的至少一种，由此，利用根据本发明实施例的特异性探针进行杂交捕获，能够有效地捕获样本中与靶序列互补的序列，即样本中已知具有甲基化位点的基因区域。由此，文库制备单元100可以适于实施前面所述的高通量测序文库构建方法。根据本发明的实施例，可以用于设计特异性探针的已知具有甲基化位点的基因区域为表I中所列出基因的至少一种的编码区和启动子区。根据本发明的实施例，所述编码区为外显子区域序列，所述启动子区为基因转录起始位点的上游2200bp到下游500bp的区域。根据本发明的实施例，所述特异性探针是采用eArray系统设计的。根据本发明的实施例，任选地，所述探针的长度为12mer。关于探针，前面已经进行了详细描述，在此，不再赘述。
测序单元200与文库制备单元100相连，可以从文库制备单元100接收所制备的样本的基因组特定区域的高通量测序文库，并对所接收的样本的基因组特定区域的高通量测序文库进行测序，从而可以获得测序结果。
数据分析单元300与测序单元200相连，可以从测序单元200接收所获得的测序结果，并且能够进一步对测序结果进行数据分析，从而基于分析结果确定样本的基因组特定区域的甲基化信息，最终实现对样本的基因组特定区域的甲基化检测。
本领域技术人员能够理解的是，可以采用本领域中已知的任何适于进行上述操作的装置作为上述各个单元的组成部件。在本文中所使用的术语“相连”应作广义理解，可以是直接相连，也可以通过中间媒介间接相连，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。
利用根据本发明实施例的用于确定样本的基因组特定区域的甲基化信息的装置，能够方便准确地确定样本的基因组特定区域的甲基化信息，从而可以应用于多种针对基因组特定区域，如已知甲基化位点的基因组区域的甲基化的研究，例如可以用于对基因组特定区域的甲基化异常进行检测。
试剂盒
根据本发明的另一方面，本发明提供了一种用于构建样本的基因组特定区域的高通量测序文库的试剂盒。根据本发明的实施例，该试剂盒包括：特异性探针，该特异性探针是对已知甲基化位点特异性的。根据本发明的一些具体示例，特异性探针是基于采用人类基因组作为参考序列，并且采用基因组上已知具有甲基化位点的特定基因区域作为靶序列而设计的，具体地，已知具有甲基化位点的基因区域包括选自启动子区域、CpG岛区域、CpG岛外区域以及印记基因区域的至少一种，由此，利用根据本发明实施例的特异性探针进行杂交捕获，能够有效地捕获样本中与靶序列互补的序列，即样本中已知具有甲基化位点的基因区域。根据本发明的实施例，可以用于设计特异性探针的已知具有甲基化位点的基因区域为表I中所列出基因的至少一种的编码区和启动子区。根据本发明的实施例，所述编码区为外显子区域序列，所述启动子区为基因转录起始位点的上游2200bp到下游500bp的区域。根据本发明的实施例，所述特异性探针是采用eArray系统设计的。根据本发明的实施例，任选地，所述探针的长度为12mer。关于探针，前面已经进行了详细描述，在此，不再赘述。
本领域的技术人员可以理解，试剂盒中还可以进一步包括构建样本的基因组特定区域的高通量测序文库所需的任何其他组分，在此不再赘述。利用根据本发明实施例的用于构建样本的基因组特定区域的高通量测序文库的试剂盒，能够方便有效地构建样本的基因组特定区域的高通量测序文库。
需要说明的是，根据本发明实施例的构建样本的基因组特定区域的高通量测序文库的方法及其应用，是本申请的发明人经过艰苦的创造性劳动和优化工作完成的。
下面将结合实施例对本发明的实施方案进行详细描述，但是本领域技术人员将会理解，下列实施例仅用于说明本发明，而不应视为限定本发明的范围。实施例中未注明具体技术或条件的，按照本领域内的文献所描述的技术或条件（例如参考J.萨姆布鲁克等著，黄培堂等译的《分子克隆实验指南》，第三版，科学出版社）或者按照产品说明书进行。所用试剂或仪器未注明生产厂商者，均为可以通过市购获得的常规产品。
实施例1：
本实施例以2μg的人类外周血单核细胞基因组DNA为样本，按照下列步骤实施。
一、基因组DNA片段化：
利用covaris-S2打断仪，按照下表设置的参数，将样本基因组DNA进行片段化处理，以便获得DNA片段。
| 处理1 | 负载比（%） | 10 |
| 强度 | 5 | |
| 循环/脉冲 | 200 | |
| 时间（s） | 50 | |
| 处理2 | 时间（s） | 0 |
| 处理3 | 时间（s） | 0 |
| 处理4 | 时间（s） | 0 |
| 循环 | 3 |
将获得的DNA片段进行电泳检测，要求DNA片段主带集中在150-300bp之间，无蛋白、RNA污染。利用QIAquick PCR纯化试剂盒（Qiagen）或磁珠纯化，将检测合格的DNA片段纯化回溶到32μl的洗脱缓冲液中，备用。
用同样的方法制备200-400ng的片段化的λ-DNA，其中λ-DNA为外源非甲基化的。
二、末端修复：
1)将上一步获得的DNA片段按照下表在1.5mL的离心管中配制末端修复反应体系：
2)将上述反应体系置于20℃的Thermomixer(Eppendorf)上，进行反应30min。反应完后用QIAquick PCR纯化试剂盒（Qiagen）进行纯化，最后将纯化产物溶于34μl洗脱缓冲液。
三、添加碱基A：
1)将上一步得到的DNA按下表在1.5mL的离心管中配制添加碱基A的反应体系：
2)将上述反应体系置于37℃的Thermomixer(Eppendorf)上，进行反应30min。反应完后用MiniElute PCR纯化试剂盒（Qiagen）进行纯化，最后将纯化产物溶于20μL洗脱缓冲液。
四．连接甲基化接头：
1)将上一步得到的DNA按下表配制甲基化接头（有时也称为“甲基化标签接头”）的连接反应体系：
*甲基化接头序列为：
接头1：5’Phos/GATCGGAAGAGCACACGTCTGAACTCCAGTCAC
接头2：5’TACACTCTTTCCCTACACGACGCTCTTCCGATCT
或采用以下标签接头进行混合杂交（Pooling杂交）：
接头1’：5’Phos/NNNNNNAGATCGGAAGAGCACACGTCTGAACTCCAGTCAC
接头2’：5’TACACTCTTTCCCTACACGACGCTCTTCCGATCTNNNNNNT
接头1和接头2或接头1’和接头2’序列中的C均进行了甲基化修饰保护。
2)将上述反应体系置于20℃的Thermomixer(Eppendorf)上，进行反应15min，获得连接产物。反应完后用MiniElute PCR纯化试剂盒（Qiagen）进行纯化，最后将纯化的连接产物溶于22μL洗脱缓冲液。
五、杂交捕获目的片段：
1、设计获得特异性探针：通过SSAHA算法，设计确定一套仅由唯一序列组成的特异性探针。具体地，以人类基因组hg19为参考序列，选取全基因组已知甲基化位点的约10,000个启动子、28,000个CpG岛、28,000个CpG岛外和61个印记基因的区域作为靶序列设计探针，其中，将少于200bp的区域通过填充补足修改为200bp的长度，并且去除重迭区域，且要求探针的序列不存在重迭序列，所有探针要满足在允许高达5个插入、缺失或错配下仍为唯一序列的条件，并且合成的每一条DNA探针序列上分别通过偶联生物素作为后续的捕获标记。然后，通过Roche NimbleGen生产获得设计好的特异性探针，备用。
表1显示了根据本发明一个实施例的特异性探针对靶区域覆盖度的评估结果。由表1可知，该探针几乎覆盖了基因组全部的启动子区域、绝大部分的印记基因和CpG岛、CGI shore区域。发明人发现，未覆盖到的区域大多是一些存在一定重复性的短序列区域，若将其增加到探针可捕获范围内，不仅会增加许多非靶区域的数据信息，同时重复序列的存在还可能影响到对其它区域的捕获效果，而这些区域的甲基化信息含量较少，且不会显著影响甲基化的整体水平，因此不将这些区域作为探针的靶序列。
2、杂交
1)将上一步得到的连接产物，按下表配制杂交反应体系：
用Qubit fluorometer和相应的Quant-iT dsDNA HS检测试剂盒(Invitrogen)对所得到的连接产物进行定量，然后取1μg的连接产物于一个新的1.5mL的EP管中，并加入10μL 1mg/mL的cot-1DNA和各1nmol的接头封闭序列，然后置于SpeedVac中于60℃下蒸干，备用。然后，在蒸干的EP管中分别加入2×SC杂交缓冲液和SC杂交组合物A。
*接头封闭序列为：
Block1：5’GTGACTGGAGTTCAGACGTGTGCTCTTCCGATC
Block2：5’AGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGTA
或采用以下接头封闭序列进行混合杂交：
Block1’：5’GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCTNNNNNN
Block2’：5’ANNNNNNAGATCGGAAGAGCGTCGTGTAGGGAAAGAGTGTA
其中，Block1’和Block2’中的碱基NNNNNN分别和接头1’和接头2’序列中的碱基N互补配对。
将上述反应体系混匀后置于离心机上全速离心10秒，然后将其转移至heat block中于95℃下孵育10分钟，使DNA变性。
4）取出上述样品，震荡混匀后于下室温全速离心10秒，然后将其转移至一个0.2mL的PCR管中或96孔PCR板中，并加入4.5μL前面所述的探针文库，震荡混匀，然后置于离心机上全速离心10秒，再将该PCR管或96孔PCR板置于PCR仪上，于47℃下杂交64h-72h，其中PCR仪的热盖温度设置为57℃。
3、序列捕获
1）准备洗涤缓冲液
a)将以下四种洗涤缓冲液进行稀释：10×SC洗涤缓冲液I、10×SC洗涤缓冲液II、10×SC洗涤缓冲液III和2×Stringent洗涤缓冲液，稀释到1×溶液后贮存，其中贮存时间不宜超过2周。
b)于47℃下预热配制好的1mL Stringent洗涤缓冲液和1mL SC洗涤缓冲液I两种溶液
2）准备链霉素磁珠
a)从冰箱中取出Dynabeads M-280Streptavidin(Invitrogen)磁珠，充分混匀后取100μL于一个新的1.5mL的EP管中；
b)将EP管置于磁力架上至澄清，用移液器小心的去除上清，然后加入200μL的Streptavidin Dynabead结合洗涤缓冲液；
c)将上述反应体系震荡（Vortex）10s，混匀后，将EP管重新放回磁力架至液体澄清，用移液器小心的去除上清；
d)将上一步得到的反应体系进行两次洗涤；
e)用100μL的Streptavidin Dynabead结合洗涤缓冲液将上述反应体系的磁珠进行重悬，并将其转入0.2mL的小管中；
f)用磁力架结合磁珠（将小管靠到磁力架上），直到液体澄清，用移液器小心地去除上清。
3）利用链霉素磁珠捕获目的片段
将杂交混合物吸出来(记录杂交后剩余体积)加入到准备好的磁珠中，吹打混匀10次后将小管放在PCR仪上，47℃下孵育45min(PCR仪热盖温度应设为57℃，每隔15min，取出震荡3s，以防止磁珠沉淀)。
4）洗涤结合了捕获DNA的链霉素磁珠
a)孵育45min后，将混合物从0.2mL的小管转入1.5mL的EP管中，将EP管置于磁力架上至液体澄清，小心的去除上清；
b)加入100μL预热到47℃的1×洗涤缓冲液I，震荡混匀10s，将EP管置于磁力架上至液体澄清，小心的去除上清；
c)从磁力架上取下EP管，加入200μL预热到47℃的1×Stringent洗涤缓冲液，移液器吹打混匀10次（该操作应迅速以便管中的液体不低于47℃）；
d)47℃下孵育5min后，将EP管置于磁力架上至液体澄清，小心的去除上清；
e)重复步骤c)-d)，共用1×Stringent洗涤缓冲液洗两次；
f)加200μL室温放置的1×洗涤缓冲液I，震荡混匀2min，如果液体溅到管盖上，用手指轻弹EP管使其集中到管底，将EP管置于磁力架上至液体澄清，小心地去除上清；
g)加200μL室温放置的1×洗涤缓冲液II，震荡混匀1min，将EP管置于磁力架上至液体澄清，小心的去除上清；
h)加200μL室温放置的1×洗涤缓冲液III，震荡混匀30s，将EP管置于磁力架上至液体澄清，小心的去除上清。
4）洗脱结合了目的片段的链霉素磁珠
a）在以上洗涤好的磁珠中加入50μLSureSelect洗脱缓冲液，振荡5s，重悬磁珠；
b）将上述反应体系于室温下孵育10min，将EP管置于磁力架上5-10min至液体澄清；
c）用移液器将上清液转移到一个新的1.5mL离心管中（此时的上清液中含有捕获的DNA，磁珠可丢弃）；
d）向上清液中加入50μLSureSelect中和缓冲液，混匀；
e）用MiniElute PCR纯化试剂盒（Qiagen）纯化后，溶于22μL的洗脱缓冲液中。
六、重亚硫酸盐处理：
以前面制备好的200-400ng片段化的λ-DNA作为外源DNA，采用EZ DNAMethylation-Gold Kit(ZYMO)，将杂交捕获的目的片段DNA和外源DNA一起进行重亚硫酸盐共处理，具体步骤如下：
1）制备CT转换试剂（CT Conversion Reagent）溶液：从试剂盒中取出CT转换试剂（固体混合物），分别加入900μL的水、50μL的M-溶解缓冲液（M-Dissolving Buffer）和300μL的M-稀释缓冲液（M-Dilution Buffer），室温下溶解并且震荡10分钟或在摇床上摇动10分钟。
2）M-洗涤缓冲液的制备：向M-洗涤缓冲液中添加24mL 100%的乙醇，备用。
3）将待转换的目标片段DNA与λ-DNA混合物加入同一PCR管中，若不足20μL的则用水补足。
4）在PCR管中加入130μL的CT转换试剂溶液，轻弹或移液器吹悬混合样品。
5）将样品管放到PCR仪上按以下步骤操作:
98℃下持续5分钟
64℃下持续2.5小时
完成上述操作后，立刻进行下一步操作或者在4℃下存储(最多20小时)备用。
6）将Zymo-Spin ICTM Column放入收集管（Collection Tube）中，并加入600μL的M-结合缓冲液（M-Binding Buffer）。
7）将重亚硫酸盐处理的样品加入到含M-结合缓冲液的Zymo-Spin IC Column中，盖上盖子颠倒混匀。
8）全速(>10,000xg)离心30秒，弃收集管中的收集液。
9）向柱中加入100μL的M-洗涤缓冲液，全速(>10,000xg)离心30秒，弃收集管中的液体。
10）向柱中添加200μL的M-Desulphonation Buffer，室温放置15min，全速(>10,000xg)离心30s，弃收集管中的液体。
11）向柱中添加200μL的M-洗涤缓冲液，全速(>10,000xg)离心30s，弃收集管中的液体，并再重复此步骤1次。
12）将Zymo-Spin ICTM Column置于新的1.5mL EP管中，加入12μL的M-洗脱缓冲液r到柱基质中，室温放置2min，全速(>10,000xg)离心洗脱目的片段DNA。
七、PCR扩增及扩增产物分离纯化：
1）将上一步得到的目的片段DNA按以下体系配制PCR反应体系：
*其中，P1公用引物的序列为：
AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT
**标签N的序列为：
CAAGCAGAAGACGGCATACGAGATNNNNNNGTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT，其中碱基N为A、T、C、G四个碱基的任意组合，作为区别标识。
PCR反应条件：
2）将PCR扩增产物经2%琼脂糖电泳后，使用QIAquick凝胶提取试剂盒(Qiagen)回收纯化300-450bp片段的文库，备用。
八、文库检测：
使用Bioanalyzer分析系统(Agilent,Santa Clara,USA)检测文库插入片段的大小及含量；并利用Q-PCR精确定量文库的浓度。
由此，构建的样本的基因组特定区域的高通量测序文库经检测合格，备用。
实施例2：
利用Hiseq2000测序仪，按照双末端90个碱基的读长，将实施例1构建的样本的基因组特定区域的高通量测序文库进行测序，以便获得测序结果。
在上述测序后，直接获得的是原始数据，通过对原始数据进行基本分析可以获得上述测序结果，其中，该基本分析过程包括以下主要步骤：首先，通过接头或PCR引物上的序列标签区分不同样本的文库数据；然后，对测序所得的原始数据进行去污染、去接头和去低质量过滤；最后，将经过前述处理的数据进行碱基转化，具体地，将正链的C全部转化成T，互补链的G全部转化成A，由此，获得实施例1构建的样本的基因组特定区域的高通量测序文库的测序结果。
将获得的测序结果进行数据分析，以便确定样本的基因组特定区域的甲基化信息。其中，数据分析包括：使用SOAP2.01软件，将测序结果中的reads与参考基因组进行比对，其中容许错配率设置为2，以便确定唯一比对上的reads，基于对这些reads的分析，获得样本的基因组特定区域的序列信息及该区域的甲基化信息。
本实施例的一个方面，以非CpG处的单个的C作为标准，计算实施例1的重亚硫酸盐处理的转换效率；以及，基于测序结果，进行测序深度和覆盖度的分析，其中，在本实施例中，是对hg19全基因组所有启动子区域、CpG岛、CGI shore和印记基因区域的覆盖度，及不同区域的覆盖深度的分析，由此确定不同覆盖区域的甲基化水平。
此外，本实施例由测序结果确定了实施例1中的特异性探针的捕获效率。图2显示了根据本发明一个实施例的方法确定基因组特定区域甲基化信息时，在不同覆盖深度下（覆盖深度≥1及覆盖深度＞5），每条染色质上的捕获区域占探针靶区域的百分比图。图2所依据的测序数据为：测序原始数据序列为25.5M，比对率为75.27%，唯一比对测序序列约为14.9M，唯一比对率为57.78%。由图2可知，在覆盖深度≥1的条件下，99%以上的探针均可检测到其捕获区域的甲基化信息，而当覆盖深度＞5的时候，也有约90%的探针可以检测到其捕获区域的甲基化信息，（≥1是指≥1的所有，其包含了＞5）同样表明，可以通过适当的增加测序数据量来进一步提高探针捕获的实际检测范围，这表明根据本发明实施例的探针可以稳定可靠的捕获靶区域，结合重亚硫酸盐处理即可准确地进行甲基化检测。
本实施例还分析了检测到甲基化信息的不同元件的实验数据占基因组该区域的百分率，分析结果如图3和表2所示。基因组经过杂交捕获和重亚硫酸盐处理之后，依据以下测序数据进行分析及作图，从而得到图3和表2：测序原始数据序列为25.5M，比对率为75.27%，唯一比对测序序列约为14.9M，唯一比对率为57.78%。图3显示了根据本发明一个实施例的方法确定基因组特定区域甲基化信息时，在不同覆盖深度下，各条染色质中检测到甲基化信息的启动子占该染色质的总启动子的百分比图。由图3可知，当覆盖深度大于5，各条染色质中可检测到甲基化信息的启动子占该染色质上总启动子的百分比均大于70%，与理论值接近，而且覆盖深度大于10时（大于5以上的所有，包含了大于10请发明人解释一下，大于5的具体范围是多少，因为大于5也包含了大于10），各条染色质中可检测到甲基化信息的启动子仍可占该染色质上总启动子的60%以上。表2显示了根据本发明一个实施例的方法确定基因组特定区域甲基化信息时，各条染色质中可检测到的印记基因在该染色质上的分布分析结果。由表2可知，当覆盖深度≥1的时候，97.6%的印记基因的甲基化信息均可被检测到，而保持测序数据量不变的条件下，随着覆盖深度过滤的增加，检测到的基因个数明显降低，这表明，在高测序深度下分析印记基因的甲基化信息时，应该加大现有的测序数据量，提高每一个印记基因的覆盖深度。
另外，本实施例还分析了基因组上启动子区域、CpG岛、CGI shore及印记基因区域的甲基化水平分布，分析结果如图4所示。
表1 设计探针在全基因组各靶区域的覆盖信息
表2 检测到的印记基因在每一条染色质上的分布分析
其中，图4(a)显示了根据本发明一个实施例的确定的样本的基因组CpG岛、CGI shore区域的甲基化水平分布图。由图4(a)可知，高CG含量的CpG岛处于低甲基化水平，而CGIshore区域的甲基化水平，相对于CpG岛甲基化水平显著增高。图4(b)显示了根据本发明一个实施例的确定的样本的基因组启动子区域的甲基化水平分布图。由图4(b)可知，启动子区域中，其转录起始位点处的甲基化水平处于低水平，所有结果与理论相符。图4(c)显示了样本的基因组特定区域的原始分布和根据本发明一个实施例的样本的基因组特定区域的高通量测序文库的reads分布及确定的启动子、CpG岛区域的甲基化水平分布图。由图4(c)可知，根据本发明实施例的确定样本的基因组特定区域的甲基化信息的方法，能有效地捕获每一个特定区域，并能够准确检测到该区域的甲基化信息。
实施例3：
采用炎黄细胞系样品（Jun Wang et al.2008），重复实施例1，只是，用于设计特异性探针的已知具有甲基化位点的基因区域为表I中所列出基因的编码区和启动子区（合并重复基因后共867个基因），采用eArray系统设计，由安捷伦公司制备的，探针的长度为12mer。另外，对于重测序和非甲基化测序文库不需要进行重亚硫酸盐处理步骤。
采用混合标签测序，读段长度49bp，标签长度6bp，下机序列片段数量2.67Mb对，测试数据量产出约240Mb。使用bwa比对程序，将过滤了低质量和污染接头的测序片段比对到人类全基因组上。并对比对结果做了初步的分析。
检测结果：
表3给出了炎黄细胞系样品具体的下机的数据总量、通过过滤和筛选之后得到的数据量，以及通过比对之后最终能够比对到人类基因组上的序列总数，以及比对率和芯片捕获效率等。
表3 数据产出以及比对基本情况统计
表4 给出了各染色体和基因元件统计的目标区域的深度和覆盖度，从整体上看，所捕获到的数据在各个染色体的覆盖情况不存在差别。
表4 比对序列在各个染色体上的分布情况
图6示出了测序下机数据插入片段长度的分布，从图中可以看出，虽然未经过片段选择，但插入片段大小在160bp左右。图7示出了目标区域各个碱基的测序深度的分布情况，从图中可以看出，大部分碱基（约75%）的覆盖深度在20X以上，如果继续加大测序量，则能保障测序覆盖深度的一个要求。图8展示了在数据比对过程中，所有比对上的Reads各个碱基位置发生错配的频率，根据测序原理，测序质量随读长的增加而降低，测序错误率随读长的增加而增加。从图8中也可以看出，末端Reads的错配数比较多，因此在后续变异检测中，应考虑末端测序质量问题。图9展示了目标基因的测序覆盖情况。从图9中可以看出，深度10X以上，覆盖度达到60%以上的基因达到了80%，说明探针对基因的捕获无误，继续加大测序量，则可满足在一定深度上，100%覆盖所有的涉及到的基因。图10展示了探针芯片对组蛋白基因HIST2H3A基因及其启动子的捕获情况。
至此，本实施例说明了利用芯片捕获表观基因组基因启动子区和外显子区域的可行性，可以用于后续的变异检测分析和胞嘧啶甲基化检测的分析。为能提高检测的准确性，建议提高上机样品量，增加测序深度。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。
Claims (10)
1.一种构建高通量测序文库的方法，其特征在于，包括以下步骤：
将基因组DNA片段化，以便获得DNA片段；
将所述DNA片段进行末端修复，以便获得经过末端修复的DNA片段；
在所述经过末端修复的DNA片段的3’末端添加碱基A，以便获得具有粘性末端A的DNA片段；
将所述具有粘性末端A的DNA片段与甲基化接头相连，以便获得连接产物；
利用特异性探针对所述连接产物进行杂交捕获，以便获得目的片段；
将所述目的片段进行PCR扩增，以便获得扩增产物；以及
分离纯化所述扩增产物，所述扩增产物构成所述高通量测序文库，
其中，
所述特异性探针是对已知甲基化位点特异性的，所述特异性探针是基于采用人类基因组作为参考序列，并且采用已知具有甲基化位点的基因区域作为靶序列而设计的，所述已知具有甲基化位点的基因区域为表I中所列出基因的至少一种的编码区和启动子区，
任选地，在进行杂交捕获之后，在进行PCR扩增之前，将所述目的片段进行重亚硫酸盐处理，以便将所述目的片段中非甲基化的胞嘧啶转换为尿嘧啶，
任选地，进一步包括从样本中提取基因组DNA的步骤，优选所述样本来源于哺乳动物、植物、和微生物的至少一种，更优选所述哺乳动物为人和小鼠的至少一种，优选所述基因组DNA为人类全血基因组DNA，更优选所述基因组DNA为外周血单核细胞基因组DNA，
优选，所述基因组DNA的量为2μg，
任选地，利用covaris-S2打断仪将基因组DNA片段化，
任选地，所述DNA片段的长度为约150-300bp，优选200-300bp，
任选地，在将所述DNA片段进行末端修复前，进一步包括纯化DNA片段的步骤，
任选地，将所述DNA片段进行末端修复是利用Klenow片段、T4DNA聚合酶和T4多核苷酸激酶进行的，其中，所述Klenow片段具有5’→3’聚合酶活性和3’→5’聚合酶活性，但缺少5’→3’外切酶活性，
任选地，将所述经过末端修复的DNA片段的3’末端添加碱基A是利用Klenow(3’-5’exo-)进行的，
任选地，所述甲基化接头中包含标签，
任选地，将所述具有粘性末端A的DNA片段与甲基化接头相连前，进一步包括对接头进行甲基化的步骤，
任选地，将所述具有粘性末端A的DNA片段与甲基化接头相连是利用T4DNA连接酶进行的，
任选地，在获得连接产物后，进一步包括对连接产物进行纯化的步骤，
任选地，所述已知具有甲基化位点的基因区域包括选自启动子区域、CpG岛区域、CpG岛外区域以及印记基因区域的至少一种，
任选地，所述编码区为外显子区域序列，所述启动子区为基因转录起始位点的上游2200bp到下游500bp的区域，
任选地，所述特异性探针是采用eArray系统设计的，
任选地，所述探针的长度为12mer，
任选地，在所述杂交捕获前，进一步包括利用cot-1DNA和接头封闭序列分别对所述连接产物和所述连接产物上的甲基化接头进行杂交封闭的步骤，优选采用过量的cot-1DNA对所述连接产物进行杂交封闭，可选地，所述接头封闭序列包括选自Block1和Block2的至少一种，
任选地，采用1μg的连接产物进行所述杂交捕获，
任选地，利用特异性探针对所述连接产物进行杂交捕获进一步包括利用链霉素磁珠捕获所述目的片段，
任选地，在将所述目的片段进行重亚硫酸盐处理之前，进一步包括将所述目的片段与片段化的λ-DNA混合，优选所述片段化的λ-DNA的量为200-400ng，更优选200ng，
任选地，将所述目的片段进行重亚硫酸盐处理是采用EZ DNA Methylation-Gold KitTM（ZYMO）进行的，
任选地，使用热启动taq DNA聚合酶进行所述PCR扩增，
任选地，分离纯化所述扩增产物是通过选自磁珠纯化、纯化柱纯化和2%的琼脂糖凝胶电泳的至少一种进行的，优选通过2%的琼脂糖凝胶电泳进行，
任选地，所述高通量测序文库的文库片段长度为300-450bp。
2.一种确定样本的基因组特定区域的甲基化信息的方法，其特征在于，包括下列步骤：
根据权利要求1所述的方法构建所述样本的基因组特定区域的高通量测序文库；
对所述样本的基因组特定区域的高通量测序文库进行测序，以便得到测序结果；以及
对所述测序结果进行数据分析，以便确定所述样本的基因组特定区域的甲基化信息。
3.根据权利要求2所述的方法，其特征在于，所述测序是利用高通量测序技术进行的。
4.根据权利要求3所述的方法，其特征在于，所述测序是利用Hiseq2000测序仪进行的。
5.一种用于确定样本的基因组特定区域的甲基化信息的装置，其特征在于，包括：
文库制备单元，所述文库制备单元用于制备样本的基因组特定区域的高通量测序文库，所述文库制备单元内设置有特异性探针；
测序单元，所述测序单元与所述文库制备单元相连，并且从所述文库制备单元接收所述样本的基因组特定区域的高通量测序文库，以便用于对所述样本的基因组特定区域的高通量测序文库进行测序，获得测序结果；以及
数据分析单元，所述数据分析单元与所述测序单元相连，并且从所述测序单元接收所述测序结果，以便对所述测序结果进行数据分析，确定所述样本的基因组特定区域的甲基化信息。
6.根据权利要求5所述的装置，其特征在于，所述特异性探针是对已知甲基化位点特异性的。
7.根据权利要求6所述的装置，其特征在于，所述特异性探针是基于采用人类基因组作为参考序列，并且采用已知具有甲基化位点的基因区域作为靶序列而设计的。
8.根据权利要求7所述的装置，其特征在于，所述已知具有甲基化位点的基因区域包括选自启动子区域、CpG岛区域、CpG岛外区域以及印记基因区域的至少一种，
任选地，所述已知具有甲基化位点的基因区域为表I中所列出基因的至少一种的编码区和启动子区，
任选地，所述编码区为外显子区域序列，所述启动子区为基因转录起始位点的上游2200bp到下游500bp的区域，
任选地，所述特异性探针是采用eArray系统设计的，
任选地，所述探针的长度为12mer。
9.一种用于构建样本的基因组特定区域的高通量测序文库的试剂盒，其特征在于，包括：
特异性探针，所述特异性探针是对已知甲基化位点特异性的。
10.根据权利要求9所述的试剂盒，其特征在于，所述特异性探针是基于采用人类基因组作为参考序列，并且采用已知具有甲基化位点的基因区域作为靶序列而设计的，可选地，所述已知具有甲基化位点的基因区域包括选自启动子区域、CpG岛区域、CpG岛外区域以及印记基因区域的至少一种，
任选地，所述已知具有甲基化位点的基因区域为表I中所列出基因的至少一种的编码区和启动子区，
任选地，所述编码区为外显子区域序列，所述启动子区为基因转录起始位点的上游2200bp到下游500bp的区域，
任选地，所述特异性探针是采用eArray系统设计的，
任选地，所述探针的长度为12mer。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210460712.2A CN103806111A (zh) | 2012-11-15 | 2012-11-15 | 高通量测序文库的构建方法及其应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210460712.2A CN103806111A (zh) | 2012-11-15 | 2012-11-15 | 高通量测序文库的构建方法及其应用 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103806111A true CN103806111A (zh) | 2014-05-21 |
Family
ID=50703476
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210460712.2A Pending CN103806111A (zh) | 2012-11-15 | 2012-11-15 | 高通量测序文库的构建方法及其应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103806111A (zh) |
Cited By (30)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104264231A (zh) * | 2014-09-30 | 2015-01-07 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN104293940A (zh) * | 2014-09-30 | 2015-01-21 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN104294370A (zh) * | 2014-09-18 | 2015-01-21 | 北京众科生物科技服务有限公司 | 一种用单条染色体构建文库方法 |
| CN104293938A (zh) * | 2014-09-30 | 2015-01-21 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN104294371A (zh) * | 2014-09-30 | 2015-01-21 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN104357918A (zh) * | 2014-11-25 | 2015-02-18 | 北京阅微基因技术有限公司 | 血浆游离dna文库的构建方法 |
| CN105297145A (zh) * | 2015-11-06 | 2016-02-03 | 艾吉泰康生物科技(北京)有限公司 | 一种遗传代谢病的筛查方法和试剂盒 |
| CN105506109A (zh) * | 2015-12-31 | 2016-04-20 | 安诺优达基因科技(北京)有限公司 | 一种单细胞简化代表性重亚硫酸氢盐测序的方法及试剂盒 |
| CN105603052A (zh) * | 2014-11-11 | 2016-05-25 | 武汉华大医学检验所有限公司 | 探针及其用途 |
| WO2016082057A1 (zh) * | 2014-11-25 | 2016-06-02 | 深圳华大基因研究院 | 游离dna测序文库的构建方法 |
| WO2016082130A1 (zh) * | 2014-11-26 | 2016-06-02 | 深圳华大基因研究院 | 一种核酸的双接头单链环状文库的构建方法和试剂 |
| WO2016082129A1 (zh) * | 2014-11-26 | 2016-06-02 | 深圳华大基因研究院 | 一种核酸的双接头单链环状文库的构建方法和试剂 |
| CN105780129A (zh) * | 2014-12-15 | 2016-07-20 | 天津华大基因科技有限公司 | 目标区域测序文库构建方法 |
| CN105861700A (zh) * | 2016-05-17 | 2016-08-17 | 湖州昂朴医学检验有限公司 | 一种针对神经肌肉病的高通量检测方法 |
| CN106191246A (zh) * | 2016-07-12 | 2016-12-07 | 集美大学 | 一种利用液相捕获进行大黄鱼基因组基因分型的方法 |
| CN106497920A (zh) * | 2016-11-21 | 2017-03-15 | 深圳华大基因研究院 | 一种用于非小细胞肺癌基因突变检测的文库构建方法及试剂盒 |
| CN106636344A (zh) * | 2016-10-28 | 2017-05-10 | 上海阅尔基因技术有限公司 | 一种基于二代高通量测序技术的地中海贫血症的基因检测试剂盒 |
| CN107109460A (zh) * | 2014-10-22 | 2017-08-29 | 艾比斯生物科学公司 | 细菌表观基因组分析 |
| WO2017193833A1 (zh) * | 2016-05-10 | 2017-11-16 | 广州嘉检医学检测有限公司 | 一种富集4000人类致病靶基因的方法及试剂盒 |
| CN109536579A (zh) * | 2018-11-05 | 2019-03-29 | 深圳市艾斯基因科技有限公司 | 单链测序文库的构建方法及其应用 |
| CN109593839A (zh) * | 2017-09-29 | 2019-04-09 | 上海交通大学 | 一种dna突变及甲基化状况检测方法 |
| CN110970086A (zh) * | 2018-09-30 | 2020-04-07 | 深圳华大生命科学研究院 | 从古dna数据中过滤现代dna污染的方法及其应用 |
| CN111041069A (zh) * | 2019-12-26 | 2020-04-21 | 人和未来生物科技（长沙）有限公司 | 一种低起始量dna样本的高通量测序文库构建方法及其应用 |
| CN111133135A (zh) * | 2017-07-17 | 2020-05-08 | 希昆斯生物科学公司 | 用于高通量测序的快速文库构建 |
| WO2020106906A1 (en) * | 2018-11-21 | 2020-05-28 | Avida Biomed, Inc. | Methods for targeted nucleic acid library formation |
| CN111690988A (zh) * | 2020-06-22 | 2020-09-22 | 上海韦翰斯生物医药科技有限公司 | 一种捕获文库构建方法及应用 |
| WO2020224159A1 (zh) * | 2019-05-06 | 2020-11-12 | 臻和精准医学检验实验室无锡有限公司 | 基于二代测序用于脑胶质瘤的检测panel、检测试剂盒、检测方法及其应用 |
| CN112877403A (zh) * | 2019-11-29 | 2021-06-01 | 福建和瑞精创基因技术有限公司 | 构建靶序列的测序文库的方法 |
| WO2021169875A1 (zh) * | 2020-02-25 | 2021-09-02 | 博尔诚（北京）科技有限公司 | 一种癌症基因甲基化检测系统和在该系统在中执行的癌症体外检测方法 |
| CN113544282A (zh) * | 2019-06-20 | 2021-10-22 | 深圳华大智造科技股份有限公司 | 基于dna样本构建测序文库的方法及应用 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101365803A (zh) * | 2005-12-22 | 2009-02-11 | 关键基因股份有限公司 | 采用高通量测序技术的改进的转录谱描述策略 |
| CN101802223A (zh) * | 2007-08-15 | 2010-08-11 | 香港大学 | 用于高通量亚硫酸氢盐dna-测序的方法和组合物及其用途 |
| CN102409408A (zh) * | 2010-09-21 | 2012-04-11 | 深圳华大基因科技有限公司 | 一种利用微量基因组dna进行全基因组甲基化位点精确检测的方法 |
-
2012
- 2012-11-15 CN CN201210460712.2A patent/CN103806111A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101365803A (zh) * | 2005-12-22 | 2009-02-11 | 关键基因股份有限公司 | 采用高通量测序技术的改进的转录谱描述策略 |
| CN101802223A (zh) * | 2007-08-15 | 2010-08-11 | 香港大学 | 用于高通量亚硫酸氢盐dna-测序的方法和组合物及其用途 |
| CN102409408A (zh) * | 2010-09-21 | 2012-04-11 | 深圳华大基因科技有限公司 | 一种利用微量基因组dna进行全基因组甲基化位点精确检测的方法 |
Non-Patent Citations (4)
| Title |
|---|
| JAMIE K. TEER等: "Systematic comparison of three genomic enrichment methods for massively parallel DNA sequencing", 《GENOME RESEARCH》 * |
| SERGEY I. NIKOLAEV ET AL.: "Detection of Genomic Variation by Selection of a 9 Mb DNA Region and High Throughput Sequencing", 《PLOS ONE》 * |
| WANG ET AL.: "High resolution profiling of human exon methylation by liquid hybridization capture-based bisulfite sequencing", 《BMC GENOMICS》 * |
| ZACHARY D. SMITH等: "High-throughput bisulfite sequencing in mammalian genomes", 《METHODS》 * |
Cited By (49)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104294370A (zh) * | 2014-09-18 | 2015-01-21 | 北京众科生物科技服务有限公司 | 一种用单条染色体构建文库方法 |
| CN104294370B (zh) * | 2014-09-18 | 2020-06-19 | 博伟至诚科技服务(北京)有限公司 | 一种用单条染色体构建文库方法 |
| CN104264231A (zh) * | 2014-09-30 | 2015-01-07 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN104293940A (zh) * | 2014-09-30 | 2015-01-21 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN104293938A (zh) * | 2014-09-30 | 2015-01-21 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN104294371A (zh) * | 2014-09-30 | 2015-01-21 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN104293940B (zh) * | 2014-09-30 | 2017-07-28 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN104294371B (zh) * | 2014-09-30 | 2017-07-04 | 天津华大基因科技有限公司 | 构建测序文库的方法及其应用 |
| CN107109460A (zh) * | 2014-10-22 | 2017-08-29 | 艾比斯生物科学公司 | 细菌表观基因组分析 |
| CN105603052A (zh) * | 2014-11-11 | 2016-05-25 | 武汉华大医学检验所有限公司 | 探针及其用途 |
| CN105603052B (zh) * | 2014-11-11 | 2021-03-19 | 武汉华大医学检验所有限公司 | 探针及其用途 |
| CN104357918A (zh) * | 2014-11-25 | 2015-02-18 | 北京阅微基因技术有限公司 | 血浆游离dna文库的构建方法 |
| CN104357918B (zh) * | 2014-11-25 | 2016-08-17 | 北京阅微基因技术有限公司 | 血浆游离dna文库的构建方法 |
| WO2016082057A1 (zh) * | 2014-11-25 | 2016-06-02 | 深圳华大基因研究院 | 游离dna测序文库的构建方法 |
| CN107075732A (zh) * | 2014-11-25 | 2017-08-18 | 深圳华大基因研究院 | 游离dna测序文库的构建方法 |
| WO2016082129A1 (zh) * | 2014-11-26 | 2016-06-02 | 深圳华大基因研究院 | 一种核酸的双接头单链环状文库的构建方法和试剂 |
| WO2016082130A1 (zh) * | 2014-11-26 | 2016-06-02 | 深圳华大基因研究院 | 一种核酸的双接头单链环状文库的构建方法和试剂 |
| US10479991B2 (en) | 2014-11-26 | 2019-11-19 | Mgi Tech Co., Ltd | Method and reagent for constructing nucleic acid double-linker single-strand cyclical library |
| CN107002291B (zh) * | 2014-11-26 | 2019-03-26 | 深圳华大智造科技有限公司 | 一种核酸的双接头单链环状文库的构建方法和试剂 |
| CN107002292A (zh) * | 2014-11-26 | 2017-08-01 | 深圳华大基因研究院 | 一种核酸的双接头单链环状文库的构建方法和试剂 |
| CN107002291A (zh) * | 2014-11-26 | 2017-08-01 | 深圳华大基因研究院 | 一种核酸的双接头单链环状文库的构建方法和试剂 |
| CN105780129A (zh) * | 2014-12-15 | 2016-07-20 | 天津华大基因科技有限公司 | 目标区域测序文库构建方法 |
| CN105297145A (zh) * | 2015-11-06 | 2016-02-03 | 艾吉泰康生物科技(北京)有限公司 | 一种遗传代谢病的筛查方法和试剂盒 |
| CN105506109A (zh) * | 2015-12-31 | 2016-04-20 | 安诺优达基因科技(北京)有限公司 | 一种单细胞简化代表性重亚硫酸氢盐测序的方法及试剂盒 |
| WO2017193833A1 (zh) * | 2016-05-10 | 2017-11-16 | 广州嘉检医学检测有限公司 | 一种富集4000人类致病靶基因的方法及试剂盒 |
| CN105861700A (zh) * | 2016-05-17 | 2016-08-17 | 湖州昂朴医学检验有限公司 | 一种针对神经肌肉病的高通量检测方法 |
| CN105861700B (zh) * | 2016-05-17 | 2019-07-30 | 上海昂朴生物科技有限公司 | 一种针对神经肌肉病的高通量检测方法 |
| CN106191246B (zh) * | 2016-07-12 | 2019-12-31 | 集美大学 | 一种利用液相捕获进行大黄鱼基因组基因分型的方法 |
| CN106191246A (zh) * | 2016-07-12 | 2016-12-07 | 集美大学 | 一种利用液相捕获进行大黄鱼基因组基因分型的方法 |
| CN106636344A (zh) * | 2016-10-28 | 2017-05-10 | 上海阅尔基因技术有限公司 | 一种基于二代高通量测序技术的地中海贫血症的基因检测试剂盒 |
| CN106497920A (zh) * | 2016-11-21 | 2017-03-15 | 深圳华大基因研究院 | 一种用于非小细胞肺癌基因突变检测的文库构建方法及试剂盒 |
| CN111133135B (zh) * | 2017-07-17 | 2024-04-16 | 希昆斯生物科学公司 | 用于高通量测序的快速文库构建 |
| CN111133135A (zh) * | 2017-07-17 | 2020-05-08 | 希昆斯生物科学公司 | 用于高通量测序的快速文库构建 |
| CN109593839A (zh) * | 2017-09-29 | 2019-04-09 | 上海交通大学 | 一种dna突变及甲基化状况检测方法 |
| CN110970086A (zh) * | 2018-09-30 | 2020-04-07 | 深圳华大生命科学研究院 | 从古dna数据中过滤现代dna污染的方法及其应用 |
| CN110970086B (zh) * | 2018-09-30 | 2023-08-15 | 深圳华大三生园科技有限公司 | 从古dna数据中过滤现代dna污染的方法及其应用 |
| CN109536579A (zh) * | 2018-11-05 | 2019-03-29 | 深圳市艾斯基因科技有限公司 | 单链测序文库的构建方法及其应用 |
| GB2595076A (en) * | 2018-11-21 | 2021-11-17 | Avida Biomed Inc | Methods for targeted nucleic acid library formation |
| GB2595076B (en) * | 2018-11-21 | 2023-07-26 | Avida Biomed Inc | Methods for targeted nucleic acid library formation |
| WO2020106906A1 (en) * | 2018-11-21 | 2020-05-28 | Avida Biomed, Inc. | Methods for targeted nucleic acid library formation |
| WO2020224159A1 (zh) * | 2019-05-06 | 2020-11-12 | 臻和精准医学检验实验室无锡有限公司 | 基于二代测序用于脑胶质瘤的检测panel、检测试剂盒、检测方法及其应用 |
| CN113544282A (zh) * | 2019-06-20 | 2021-10-22 | 深圳华大智造科技股份有限公司 | 基于dna样本构建测序文库的方法及应用 |
| CN113544282B (zh) * | 2019-06-20 | 2024-05-14 | 深圳华大智造科技股份有限公司 | 基于dna样本构建测序文库的方法及应用 |
| CN112877403A (zh) * | 2019-11-29 | 2021-06-01 | 福建和瑞精创基因技术有限公司 | 构建靶序列的测序文库的方法 |
| CN112877403B (zh) * | 2019-11-29 | 2023-11-03 | 福建和瑞基因科技有限公司 | 构建靶序列的测序文库的方法 |
| CN111041069B (zh) * | 2019-12-26 | 2021-01-19 | 人和未来生物科技（长沙）有限公司 | 一种低起始量dna样本的高通量测序文库构建方法及其应用 |
| CN111041069A (zh) * | 2019-12-26 | 2020-04-21 | 人和未来生物科技（长沙）有限公司 | 一种低起始量dna样本的高通量测序文库构建方法及其应用 |
| WO2021169875A1 (zh) * | 2020-02-25 | 2021-09-02 | 博尔诚（北京）科技有限公司 | 一种癌症基因甲基化检测系统和在该系统在中执行的癌症体外检测方法 |
| CN111690988A (zh) * | 2020-06-22 | 2020-09-22 | 上海韦翰斯生物医药科技有限公司 | 一种捕获文库构建方法及应用 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103806111A (zh) | 高通量测序文库的构建方法及其应用 | |
| CN103103624B (zh) | 高通量测序文库的构建方法及其应用 | |
| CN109797436B (zh) | 一种测序文库构建方法 | |
| CN106591441B (zh) | 基于全基因捕获测序的α和/或β-地中海贫血突变的检测探针、方法、芯片及应用 | |
| CN110036117B (zh) | 通过多联短dna片段增加单分子测序的处理量的方法 | |
| DK2002017T3 (en) | High-capacity detection of molecular markers based on restriction fragments | |
| CN103088433B (zh) | 全基因组甲基化高通量测序文库的构建方法及其应用 | |
| US11965157B2 (en) | Compositions and methods for library construction and sequence analysis | |
| Blumenstiel et al. | Targeted exon sequencing by in‐solution hybrid selection | |
| US20100029498A1 (en) | Selection of nucleic acids by solution hybridization to oligonucleotide baits | |
| US20130244885A1 (en) | High-throughput sequencing method for methylated dna and use thereof | |
| CN109536579B (zh) | 单链测序文库的构建方法及其应用 | |
| CN102409047B (zh) | 一种构建杂交测序文库的方法 | |
| CN113308562B (zh) | 棉花全基因组40k单核苷酸位点及其在棉花基因分型中的应用 | |
| CN102796808A (zh) | 甲基化高通量检测方法 | |
| CN109576346B (zh) | 高通量测序文库的构建方法及其应用 | |
| WO2012068919A1 (zh) | DNA文库及其制备方法、以及检测SNPs的方法和装置 | |
| JP2020516281A (ja) | 試料核酸にアダプターを付着する方法 | |
| US20160215331A1 (en) | Flexible and scalable genotyping-by-sequencing methods for population studies | |
| WO2012000150A1 (zh) | Hla-a，b基因分型用pcr引物及其使用方法 | |
| CN112410331A (zh) | 带分子标签和样本标签的接头及其单链建库方法 | |
| CN114729349A (zh) | 条码化核酸用于检测和测序的方法 | |
| CN113166809B (zh) | 一种dna甲基化检测的方法、试剂盒、装置和应用 | |
| CN112259165A (zh) | 用于检测微卫星不稳定性状态的方法及系统 | |
| CN105603052B (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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140521 |