CN103080336A - 检测胚胎或肿瘤染色体拷贝数的试剂盒、装置和方法 - Google Patents
检测胚胎或肿瘤染色体拷贝数的试剂盒、装置和方法 Download PDFInfo
- Publication number
- CN103080336A CN103080336A CN2011800355107A CN201180035510A CN103080336A CN 103080336 A CN103080336 A CN 103080336A CN 2011800355107 A CN2011800355107 A CN 2011800355107A CN 201180035510 A CN201180035510 A CN 201180035510A CN 103080336 A CN103080336 A CN 103080336A
- Authority
- CN
- China
- Prior art keywords
- dna
- chromosome
- ratio
- measured
- dna fragmentation
- 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.)
- Granted
Links
- 210000000349 chromosome Anatomy 0.000 title claims abstract description 259
- 206010028980 Neoplasm Diseases 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 48
- 210000001161 mammalian embryo Anatomy 0.000 title claims abstract description 45
- 238000012163 sequencing technique Methods 0.000 claims abstract description 54
- 210000004369 blood Anatomy 0.000 claims abstract description 38
- 239000008280 blood Substances 0.000 claims abstract description 38
- 108020004414 DNA Proteins 0.000 claims description 337
- 238000006062 fragmentation reaction Methods 0.000 claims description 144
- 238000013467 fragmentation Methods 0.000 claims description 140
- 210000002381 plasma Anatomy 0.000 claims description 86
- 239000012634 fragment Substances 0.000 claims description 69
- 108090000623 proteins and genes Proteins 0.000 claims description 52
- 238000001514 detection method Methods 0.000 claims description 45
- 230000008774 maternal effect Effects 0.000 claims description 40
- 238000012937 correction Methods 0.000 claims description 39
- 239000003153 chemical reaction reagent Substances 0.000 claims description 37
- 210000001766 X chromosome Anatomy 0.000 claims description 31
- 238000010186 staining Methods 0.000 claims description 24
- 238000012408 PCR amplification Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 108700026220 vif Genes Proteins 0.000 claims description 8
- 108091028043 Nucleic acid sequence Proteins 0.000 claims description 7
- 238000000926 separation method Methods 0.000 claims description 6
- 238000001712 DNA sequencing Methods 0.000 claims description 5
- 238000000053 physical method Methods 0.000 claims description 3
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims 1
- 210000000601 blood cell Anatomy 0.000 abstract 1
- 208000037280 Trisomy Diseases 0.000 description 28
- 238000003752 polymerase chain reaction Methods 0.000 description 21
- 230000002759 chromosomal effect Effects 0.000 description 13
- 238000012360 testing method Methods 0.000 description 11
- 230000008859 change Effects 0.000 description 9
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 8
- 238000009826 distribution Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 8
- 238000004519 manufacturing process Methods 0.000 description 8
- 206010016825 Flushing Diseases 0.000 description 6
- 238000002669 amniocentesis Methods 0.000 description 6
- 238000011010 flushing procedure Methods 0.000 description 6
- 238000007400 DNA extraction Methods 0.000 description 5
- 230000001605 fetal effect Effects 0.000 description 5
- 230000000869 mutational effect Effects 0.000 description 5
- 239000007787 solid Substances 0.000 description 5
- HEDRZPFGACZZDS-UHFFFAOYSA-N Chloroform Chemical compound ClC(Cl)Cl HEDRZPFGACZZDS-UHFFFAOYSA-N 0.000 description 4
- 102000012410 DNA Ligases Human genes 0.000 description 4
- 108010061982 DNA Ligases Proteins 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 4
- 208000036878 aneuploidy Diseases 0.000 description 4
- 231100001075 aneuploidy Toxicity 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000007847 digital PCR Methods 0.000 description 4
- 201000010099 disease Diseases 0.000 description 4
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 4
- 210000003754 fetus Anatomy 0.000 description 4
- 238000000703 high-speed centrifugation Methods 0.000 description 4
- PHTQWCKDNZKARW-UHFFFAOYSA-N isoamylol Chemical compound CC(C)CCO PHTQWCKDNZKARW-UHFFFAOYSA-N 0.000 description 4
- 210000004881 tumor cell Anatomy 0.000 description 4
- 102000016928 DNA-directed DNA polymerase Human genes 0.000 description 3
- 108010014303 DNA-directed DNA polymerase Proteins 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 3
- 239000000499 gel Substances 0.000 description 3
- 239000000178 monomer Substances 0.000 description 3
- 239000011049 pearl Substances 0.000 description 3
- 238000002360 preparation method Methods 0.000 description 3
- 238000000746 purification Methods 0.000 description 3
- 102000053602 DNA Human genes 0.000 description 2
- 102000004190 Enzymes Human genes 0.000 description 2
- 108090000790 Enzymes Proteins 0.000 description 2
- 108091005804 Peptidases Proteins 0.000 description 2
- 239000004365 Protease Substances 0.000 description 2
- 102100037486 Reverse transcriptase/ribonuclease H Human genes 0.000 description 2
- VMHLLURERBWHNL-UHFFFAOYSA-M Sodium acetate Chemical compound [Na+].CC([O-])=O VMHLLURERBWHNL-UHFFFAOYSA-M 0.000 description 2
- 238000000246 agarose gel electrophoresis Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 239000002585 base Substances 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000013020 embryo development Effects 0.000 description 2
- 238000004945 emulsification Methods 0.000 description 2
- 235000019441 ethanol Nutrition 0.000 description 2
- UYTPUPDQBNUYGX-UHFFFAOYSA-N guanine Chemical compound O=C1NC(N)=NC2=C1N=CN2 UYTPUPDQBNUYGX-UHFFFAOYSA-N 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 150000002989 phenols Chemical class 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 235000017281 sodium acetate Nutrition 0.000 description 2
- 239000001632 sodium acetate Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 108091032973 (ribonucleotides)n+m Proteins 0.000 description 1
- 208000031404 Chromosome Aberrations Diseases 0.000 description 1
- 206010067477 Cytogenetic abnormality Diseases 0.000 description 1
- 230000004544 DNA amplification Effects 0.000 description 1
- 201000010374 Down Syndrome Diseases 0.000 description 1
- 201000006360 Edwards syndrome Diseases 0.000 description 1
- 108010033040 Histones Proteins 0.000 description 1
- 208000026350 Inborn Genetic disease Diseases 0.000 description 1
- 108010047956 Nucleosomes Proteins 0.000 description 1
- 241001494479 Pecora Species 0.000 description 1
- 208000020584 Polyploidy Diseases 0.000 description 1
- 208000007159 Trisomy 18 Syndrome Diseases 0.000 description 1
- 230000000890 antigenic effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 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 description 1
- 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 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 238000011304 droplet digital PCR Methods 0.000 description 1
- 238000004043 dyeing Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001962 electrophoresis Methods 0.000 description 1
- 230000001804 emulsifying effect Effects 0.000 description 1
- 230000001973 epigenetic effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000012458 free base Substances 0.000 description 1
- 238000001502 gel electrophoresis Methods 0.000 description 1
- 208000016361 genetic disease Diseases 0.000 description 1
- 238000012165 high-throughput sequencing Methods 0.000 description 1
- 238000005342 ion exchange Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000010977 jade Substances 0.000 description 1
- 108020004999 messenger RNA Proteins 0.000 description 1
- 239000011859 microparticle Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 108020004707 nucleic acids Proteins 0.000 description 1
- 102000039446 nucleic acids Human genes 0.000 description 1
- 150000007523 nucleic acids Chemical class 0.000 description 1
- 210000001623 nucleosome Anatomy 0.000 description 1
- 230000001717 pathogenic effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000026731 phosphorylation Effects 0.000 description 1
- 238000006366 phosphorylation reaction Methods 0.000 description 1
- 210000002826 placenta Anatomy 0.000 description 1
- 230000003169 placental effect Effects 0.000 description 1
- 238000003793 prenatal diagnosis Methods 0.000 description 1
- 239000000376 reactant Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 210000002966 serum Anatomy 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- -1 thus standard curve Substances 0.000 description 1
- 206010053884 trisomy 18 Diseases 0.000 description 1
- 231100000588 tumorigenic Toxicity 0.000 description 1
- 230000000381 tumorigenic effect Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
- C12Q1/6874—Methods for sequencing involving nucleic acid arrays, e.g. sequencing by hybridisation
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
- G16B25/10—Gene or protein expression profiling; Expression-ratio estimation or normalisation
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6809—Methods for determination or identification of nucleic acids involving differential detection
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Organic Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Engineering & Computer Science (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Genetics & Genomics (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biotechnology (AREA)
- General Engineering & Computer Science (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Microbiology (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Medical Informatics (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Apparatus Associated With Microorganisms And Enzymes (AREA)
Abstract
本发明提供了检测胚胎或肿瘤细胞染色体拷贝数的试剂盒及其装置和方法。该方法包括取母体或肿瘤患者血液，分离血细胞和血浆，从血浆中提取DNA，构建测序文库和测序。
Description
检测胚胎或肿瘤染色体拷贝数的试剂盒、 装置和方法 技术领域 本发明涉及检测胚胎或肿瘤染色体拷贝数的试剂盒、 装置以及方法。 背景技术 染色体拷贝数异常和人类的疾病有着密切的关系。 在携带遗传性疾病的 胚胎细包中以及肿瘤细包中均有染色体异常的存在。 平均每 1000 个新生儿 中， 有 9个会携带因染色体拷贝数异常而导致的疾病（1)。 因此在胎儿尚未出 生之前能够检测出染色体拷贝数是很重要的。 然而， 目前所用的诊断方法， 包括羊膜穿刺和羊绒毛膜取样， 都属于有创伤的方法， 会给孕妇和胎儿带来 一定的风险。 用血清蛋白标记物和超声波来检测胎儿是否有染色体拷贝数异 常的疾病虽无创伤， 但不是直接检测致病因素， 因而准确性和灵敏度都不太 好 (2), 也存在不能尽早地发现染色体拷贝数异常疾病的问题。 这样的现状促 使研究人员去开发一种精确而且灵敏度高的无创诊断和检测方法。 自从母体血液中的胚胎 DNA被发现之后（3), 无创伤地直接诊断和检测 胚胎染色体异常性成了一个重大的研究课题。 2007年， 卢煜明教授和他的同 事们证明可以用母体血浆 mRNA中胎盘特异基因 4的突变位点比例来断定胚 胎是否有 21号染色体 3倍体 (4)。 突变位点比例同时也被用来判断 18号染色 体是否为 3倍体 (5)。 它的局限性在于突变位点在人群中并不普遍， 所以这些 方法只适用于一部分人群。 在同一时期， 数码 PCR ( digital PCR, 或 dPCR ) 被用来检测胚胎的染色体 3倍体 (6),(7)。 数码 PCR的优势是不依赖任何突变 位点， 但它的精确度不够高， 而且需要大量血样， 加大了釆样的困难。 近年来迅猛发展的高通量测序技术解决了以上的问题。 这些技术包括 Illumina公司的 Genome Analyzer(8), Life Technologies公司的 SOLiD(9), 以 及 Helicos公司的 Heliscope(lO), 能一次性地检测出几亿乃至几十亿个序列。 用这些技术来检测母体血浆 DNA时，血浆中微量的胚胎 DNA的染色体数目 的变化能被检测出来（11),(12),(13)。但因为测序的成本太高， 目前这些技术尚 未被普遍使用。 同时， 从母体血液中检测胚胎染色体的局部拷贝数的变化还 属于未解决的难题。 用高通量测序来检测母体血浆中胚胎染色体的拷贝数变 化有一些优势， 但目前价钱昂贵， 不能普及。 而且测序的偏差（CV )比较高， 检测的精确度和稳定性也都有待改进。 测序的偏差也决定了这种方法只适合 少数几个染色体， 如 21号染色体， 18号染色体， 而且目前还不适合于染色
体局部拷贝数变化的检测。 用高通量测序来检测染色体数目的变化的高成本和难度主要因为母体血 浆中胚胎 DNA的含量比较低， 低时只占 5%, 尤其在胚胎发育的早期。 母体 血浆中的大部分 DNA还是母体 DNA。 胚胎染色体数目或局部拷贝数的变化 艮容易被母体 DNA的背景所覆盖。因此分离母体和胚胎 DNA的方法也就成 为多年来研究的课题， 但成效不大。 比较成功的方法应属 Baylor医学院发明 的针对组蛋白的分离方法（14)。 不过分离出来的 DNA量太少， 只适合于突变 位点的检测， 不宜用来检测染色体拷贝数的变化。 发明内容 针对上面检测胚胎染色体拷贝数的方法中的种种问题， 发明人设计出能 有效地氐成本检测出胚胎染色体全部或局部拷贝数的试剂盒、 装置及方法。 本发明是基于以下的事实： 发明人发现， 通过本发明的方法所测得的来 自各个染色体的 DNA片段的 GC含量分别与来自各个染色体的 DNA片段占 总 DNA片段的比值具有一定的线性关系， 上述现象可能与检测的方法相关， 该线性关系可用 y = ax + b表示，其中 y代表来自待测染色体的 DNA片段的 GC含量， X代表来自待测染色体的 DNA片段数量占总 DNA的比值， a和 b 是常数， 对于不同的染色体 a和 b可以是不同的值， 可根据所述来自待测染 色体的 DNA片段中的 GC含量对所述比值进行校正， 并计算待测样本中所 述来自待测染色体的 DNA 片段校正后的比值的变异， 居所述变异的程度 确定待测染色体的拷贝数。 通过 GC含量的校正， 能够有效地检测出原来仅 仅靠各个染色体的 DNA片段占总 DNA片段比值的判断方法检测不出来的多 个假阴性结果。 在具体实施方式部分有具体的实验为证。 另夕卜， 如在文献中 4艮道的（15), 母体血浆中的胚胎 DNA大部分为 lOObp 到 250bp的片段， 且 150bp到 170bp的 DNA尤其占多数， 虽然母体的 DNA 也有很小部分分布在这个片段区域， 但片段大于 250bp的 DNA基本上属于 母体的 DNA。 本发明人发现， 虽然理由未知， 来自各个染色体的 DNA片段 占总 DNA片段的比值在 lOObp到 250bp之间的任意一点或任意一个区间处 的 DNA是均匀分布的， 即各个染色体在 lOObp到 250bp之间的任何一点， 比如 HObp或 167bp (此位点的 DNA量最多；)， 与总 DNA的比值代表了其 他点的比值， 因此， 也就代表了各个染色体占 lOObp 到 250bp 之间的所有 DNA的比值。 如果通过对所述 DNA中 100bp-250bp之间的任意一点或任意 一个区间的 DNA进行测序， ^！夺所述 DNA中 100bp-250bp之间的任意一点或 任意一个区间的所有 DNA 的测序结果与染色体组序列图谱进行比对， 以确
定所述 DNA中 100bp-250bp之间的任意一点或任意一个区间的所有 DNA中 的每段 DNA来自几号染色体和所述每段 DNA的序列长度；并计算在同一样 本内所述 DNA中 100bp-250bp之间的任意一点或任意一个区间的所有 DNA 中来自待测染色体的 DNA片段占 100bp-250bp之间的任意一点或任意一个 区间的所有 DNA片段的比值，也就得到了各个胚胎染色体与总 DNA的比值。 这样大大降低了检测的结果。 结合上面的根据 GC校正的方法， 根据所述来 自待测染色体的 DNA片段中的 GC含量对所述比值进行校正， 并计算待测 样本中所述来自待测染色体的 DNA 片段校正后的比值的变异， 根据所述变 异的程度确定待测染色体的拷贝数。 同时， 发明人发现在肿瘤发生的过程中， 同样地， 患者的血液中出现了 与胚胎发育过程中母体血液中类似的现象， 即在肿瘤患者的血液中， 能够检 测到游离的肿瘤细胞的 DNA。通过本发明的方法所测得的来自各个染色体的 DNA片段的 GC含量分别与来自各个染色体的 DNA片段占总 DNA片段的 比值具有一定的线性关系应同样适用于肿瘤细胞的 DNA 的染色体多倍体的 检测。 而且血浆中游离肿瘤细胞的 DNA 以核小体的形式存在， 所以它们也 大部分为 lOObp到 250bp的片段，而且来自各个染色体的 DNA片段占总 DNA 片段的比值在 lOObp到 250bp之间的任意一点或任意一个区间处的 DNA是 均匀分布的， 即各个染色体在 lOObp到 250bp之间的任何一点， 与总 DNA 的比值代表了其他点的比值，因此，也就代表了各个染色体占 lOObp到 250bp 之间的所有 DNA的比值。 因此， 本发明的试剂盒、 装置和方法也同样适用 于检测肿瘤细胞染色体或局部染色体的拷贝数。 基于以上的发现， 本发明人发明了相对无创且经济方便且更准确地检测 胚胎或肿瘤染色体或局部染色体拷贝数的试剂盒、 装置和方法。 本发明提供的一种检测胚胎或肿瘤染色体拷贝数的试剂盒， 包括： 从母 体或肿瘤患者体内取血液的器械； 适合将所述血液中的血细胞和血浆分离的 器械； 抽提所述血浆中的 DNA的试剂和器械； 通过物理方法根据所述 DNA 的片段大小对其进行分离的试剂和器械； 和取 100bp-250bp之间的任意一点 或任意一个区间的所有 DNA进行测序的试剂和器械。。 优选地， 所述胚胎或肿瘤染色体是整条染色体或局部染色体。 优选地， 本发明的试剂盒还包括： 将所述 DNA制成供测序的文库的试 剂和器械。 优选地， 本发明的试剂盒还包括： 对所述从血浆中抽提的 DNA或所述 供测序的文库进行 PCR扩增的试剂和器械。
优选地， 其中所述 100-250bp 之间的任意一点或任意一个区间的 DNA 是 150bp-170bp的 DNA, 更优选地， 是 167bp的 DNA。 本发明提供的另一种检测胚胎或肿瘤染色体拷贝数的试剂盒， 包括： 从 母体或肿瘤患者体内取血液的器械； 适合将所述血液中的血细胞和血浆分离 的器械； 抽提所述血浆中的 DNA的试剂和器械； 将所述 DNA制成可供测 序的文库的试剂和器械； 和对所述 DNA进行测序的试剂和器械。 其中所述 胚胎或肿瘤染色体是整条染色体或局部染色体。 优选地，本发明的试剂盒还包括：对从血浆中抽提的所述 DNA进行 PCR 扩增的试剂和器械。 本发明提供的一种检测胚胎或肿瘤染色体拷贝数的装置， 包括： 检测模 块， 用于对母体血浆或肿瘤患者血浆样本中的 DNA进行测序， 所述测序包 括了将母体血浆或肿瘤患者血浆样本的所有 DNA制成可供测序的文库的过 程； 比对模块， 用于将所述 DNA的测序结果与染色体组序列图谱进行比对， 以确定所述 DNA中的每段 DNA来自几号染色体以及所述每段 DNA的序列 长度； 计算模块， 用于计算在同一样本内来自待测染色体的 DNA 片段的数 量占 DNA片段总数的比值， 才艮据所述来自待测染色体的 DNA片段中的 GC 含量对所述比值进行校正； 并计算待测样本中所述来自待测染色体的 DNA 片段校正后的比值的变异， 根据所述变异的程度确定待测染色体的拷贝数； 以及输出模块， 用于输出所述待测染色体的拷贝数。 优选地， 本发明的装置中的计算模块根据以下的函数对第 2、 3、 4、 5、
6、 7、 8、 12、 13、 18和 X号染色体的比值进行校正： 来自第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18和 X号染色体的 DNA片段的 GC含量分别与来自 各个染色体的 DNA片段占总 DNA片段的比值具有一定的线性关系，所述线 性关系可用 y = ax + b表示， 其中 y代表来自待测染色体的 DNA片段的 GC 含量， X代表来自待测染色体的 DNA片段数量占总 DNA的比值， a和 b是 常数， 且 a是负数。 优选地， 所述来自第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18和 X号染 色体的 DNA片段的 GC含量与来自各个染色体的 DNA片段占总 DNA片段 的比值之间的线性关系如在表 1中所示。 优选地， 本发明的装置中的计算模块根据以下的函数对第 1、 9、 10、 11、
15、 16、 17、 19、 20、 21和 22号染色体的比值进行校正： 来白第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21和 22号染色体的 DNA片段的 GC含量分别与 来自各个染色体的 DNA片段占总 DNA片段的比值具有一定的线性关系，所
述线性关系可用 y = ax + b表示，其中 y代表来自待测染色体的 DNA片段的 GC含量， X代表来自待测染色体的 DNA片段数量占总 DNA的比值， a和 b 是常数， 且 a是正数。 优选地， 所述来白 白第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21和 22 号染色体的 DNA片段的 GC含量与来自各个染色体的 DNA片段占总 DNA 片段的比值之间的线性关系如在表 2中所示。 优选地，所述测序包括了将母体血浆或肿瘤患者血浆样本中的所有 DNA 制成可供测序的文库的过程。 优选地， 所述对母体血浆或肿瘤患者血浆样本中的 DNA进行测序是双 端短序列测序、 单端长序列测序或单端短序列测序。 优选地， 其中所述胚胎或肿瘤染色体是整条染色体或局部染色体。 优选地， 其中对所述母体血浆或肿瘤患者血浆样本中的 DNA在测序前 进行了 PCR扩增。 优选地，本发明的装置中，计算模块用于计算在同一样本内 100bp-250bp 之间的任意一点或任意一个区间的所有 DNA中来自待测染色体的 DNA片段 数量占 100bp-250bp之间的任意一点或任意一个区间的所有 DNA片段总数 的比值， 才艮据 100bp-250bp之间的任意一点或任意一个区间的所有 DNA 中 来自待测染色体的 DNA片段的 GC含量对所述比值进行校正； 并计算待测 样本中所述来自待测染色体的 DNA 片段校正后的比值的变异， 根据所述变 异的程度确定待测染色体的拷贝数。。 优选地， 其中将所述母体血浆或肿瘤患者血浆样本中 100bp-250bp之间 的任意一点或任意一个区间的所有 DNA制成可供测序的文库。 优选地， 其中所述母体血浆或肿瘤患者血浆样本中的 DNA在被制作成 可供测序的文库之前或之后进行了 PCR扩增。 优选地， 其中所述 100-250bp 之间的任意一点或任意一个区间的 DNA 是 150bp-170bp的 DNA, 更优选地， 是 167bp的 DNA。 本发明提供的一种检测胚胎或肿瘤染色体拷贝数的方法，包括以下步骤： 取母体血浆或肿瘤患者血浆； 将所述血液中的血浆和血细胞分离； 将所述血 浆中的 DNA制作成可供测序的文库； 对所述 DNA测序文库进行测序； 将所 述测序的结果与染色体组序列图谱进行比对， 确定每段所述 DNA序列来自
几号染色体以及每段所述 DNA序列的长度；以及用所述 DNA的测序和比对 结果，计算在同一样本内来自待测染色体的 DNA片段的数量占 DNA片段总 数的比值， 才艮据所述来自待测染色体的 DNA片段中的 GC含量对所述比值 进行校正， 并计算待测样本中所述来自待测染色体的 DNA 片段校正后的比 值的变异， 根据所述变异的程度确定待测染色体的拷贝数。 优选地， 本发明的方法才艮据以下的函数对第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18和 X号染色体的比值进行校正： 来自第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18和 X号染色体的 DNA片段的 GC含量分别与来自各个染色 体的 DNA片段占总 DNA片段的比值具有一定的线性关系，所述线性关系可 用 = & + 1？表示， 其中 y代表来自待测染色体的 DNA片段的 GC含量， x 代表来自待测染色体的 DNA片段数量占总 DNA的比值， a和 b是常数， 且 a是负数。 优选地， 其中所述来自第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18和 X 号染色体的 DNA片段的 GC含量与来自各个染色体的 DNA片段占总 DNA 片段的比值之间的线性关系如在表 1中所示。 优选地， 本发明的方法才艮据以下的函数对第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21和 22号染色体的比值进行校正： 来白第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21和 22号染色体的 DNA片段的 GC含量分别与来自各个染色 体的 DNA片段占总 DNA片段的比值具有一定的线性关系，所述线性关系可 用 = & + 1？表示， 其中 y代表来自待测染色体的 DNA片段的 GC含量， x 代表来自待测染色体的 DNA片段数量占总 DNA的比值， a和 b是常数， 且 a是正数。 优选地， 其中所述来自第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21 和 22号染色体的 DNA片段的 GC含量与来自各个染色体的 DNA片段占总 DNA 片段的比值之间的线性关系如在表 2中所示。 优选地， 仅耳又 DNA的序列长度在 100bp-250bp之间的任意一点或任意 一个区间的所有 DNA 的测序和比对结果， 计算在同一样本内 100bp-250bp 之间的任意一点或任意一个区间的所有 DNA中来自待测染色体的 DNA片段 数量占 100bp-250bp之间的任意一点或任意一个区间的所有 DNA片段总数 的比值， 才艮据 100bp-250bp之间的任意一点或任意一个区间的所有 DNA 中 所述来自待测染色体的 DNA片段中的 GC含量对所述比值进行校正， 并计 算待测样本中所述来自待测染色体的 DNA 片段校正后的比值的变异， 才艮据 所述变异的程度确定待测染色体的拷贝数。
优选地， 对所述 DNA测序文库进行测序是双端短序列测序、 单端长序 列测序或单端短序列测序。 优选地， 所述 100-250bp 之间的任意一点或任意一个区间的 DNA 是 150bp-170bp的 DNA, 更优选地， 是 167bp的 DNA。 附图说明 图 1是将母体血浆的 DNA制成可供双端测序的文库后，用 1%的琼脂糖 凝胶电泳后的图像。 左边的通道是 DNA lOObp标准样的图像， 右边的通道是 可供双端测序的文库 DNA的图像， 最明显的条带位于 280bp左右， 其中包 含了 120bp的连接引物。 图 2是样品 G356的 DNA片段大小的分布图。 图 3是比对后样品 G356中确定来自 X染色体的 DNA片段根据片段大 小所故的分布图。 图 4是计算得出的和实际测得的 G356的 X染色体的 DNA片段根据片 段大小所做的分布图， 图中带星号的线表示 G356样品的总 DNA序列数乘以 X染色体的百分比后得到的图形，圓圏代表实际检测到的 X染色体的序列数。 如图中所示， 计算得出的 G356X染色体的 DNA片段的分布与实际测量得到 的分布基本相同， 两条线基本吻合。 图 5A-H是根据样品代号为： G356、 G397(重复测了 8次）、 G426、 G735、 G756、 G760、 G763、 G770、 G778、 G779、 G780、 G781、 G824 和 G825 中来自待测染色体的 DNA片段的 GC含量和来自待测染色体的 DNA片段数 量占总 DNA片段的比值，所故的标准曲线，样品 G356、 G397、 G426、 G735、 G756、 G760、 G763、 G770、 G778、 G779、 G780、 G781、 G824 和 G825 均由中南大学湘雅医学院通过羊膜腔穿刺经染色体核型鉴定为正常的女性胚 胎样本 [46 , XX]。其中， X轴表示来自待测染色体的 DNA片段数量占总 DNA 片段的比值， Y轴表示来自待测染色体的 DNA片段的 GC含量。 5 A-F每幅 图示出了针对三个染色体的标准曲线，依次示出了第 1-18号染色体的标准曲 线图， 图 5G示出了第 19-22号染色体的标准曲线图， 图 5H示出了 X染色 体的标准曲线图。针对每个染色体的标准曲线的函数，如表 1和表 2中所示。 图 6是用本发明的方法检测出来的第 13号染色体三体的样品的示意图。 其中， X轴表示来自 13号染色体的 DNA片段数量占总 DNA片段的比值， Y轴表示来自 13号染色体的 DNA片段的 GC含量。在 X轴上用菱形示出了
各个样品中来自第 13号染色体的 DNA片段数量占总 DNA片段的比值， 而 用箭头指示出来的矩形表示经来自第 13号染色体的 DNA片段的 GC含量校 正后， 检测出的第 13 号染色体三体的样品。 如果不釆用本发明的待测染色 体的 DNA片段的 GC含量的校正， 而仅从来自第 13号染色体的 DNA片段 占总 DNA 片段的比值来看， 根本发现不了箭头所指示的三个样品的异常， 则会出现假阴性的结果。 图 7是用本发明的方法检测出来的第 18号染色体三体的样品的示意图。 其中， X轴表示来自 18号染色体的 DNA片段数量占总 DNA片段的比值， Y轴表示来自 18号染色体的 DNA片段的 GC含量。在 X轴上用菱形示出了 各个样品中来自第 18号染色体的 DNA片段数量占总 DNA片段的比值， 而 用箭头指示出来的矩形表示经来自第 18号染色体的 DNA片段的 GC含量校 正后， 检测出的第 18 号染色体三体的样品。 如果不釆用本发明的待测染色 体的 DNA片段的 GC含量的校正， 而仅从来自第 18号染色体的 DNA片段 占总 DNA片段的比值来看， 仅仅能发现来自第 18号染色体的 DNA片段占 总 DNA片段的比值大于等于 0.031的 5个样品是第 18号染色体三体， 不能 发现比值在 0.03左右的两个样品的异常， 对于这两个样品， 就会出现假阴性 的结果。 图 8是用本发明的方法检测出来的 X染色体单体或三体样品的示意图。 其中， X轴表示来自 X染色体的 DNA片段数量占总 DNA片段的比值， Y 轴表示来自 X染色体的 DNA片段的 GC含量。 在 X轴上用菱形示出了各个 样品中来自 X染色体的 DNA片段数量占总 DNA片段的比值， 而用箭头指 示出来的在标准曲线左侧的矩形表示经来自 X染色体的 DNA片段的 GC含 量校正后， 检测出的 X染色体单体的样品， 用箭头指示出来的在标准曲线右 侧的矩形表示经来自 X染色体的 DNA片段的 GC含量校正后， 检测出的 X 染色体三体的样品。 如果不釆用本发明的待测染色体的 DNA片段的 GC含 量的校正， 而仅从来自 X染色体的 DNA片段占总 DNA片段的比值来看， 仅仅能发现 X染色体三体的一个样品，很难发现 X染色体单体的样品， 就会 出现假阴性的结果。 图 9A-B是用本发明的方法检测出来的第 21号染色体三体的样品的示意 图。 其中， X轴表示来自 21号染色体的 DNA片段数量占总 DNA片段的比 值， Y轴表示来自 21号染色体的 DNA片段的 GC含量。 在 X轴上用菱形示 出了各个样品中来自第 21号染色体的 DNA片段数量占总 DNA片段的比值， 而用箭头指示出来的矩形表示经来自第 21号染色体的 DNA片段的 GC含量 校正后， 检测出的第 21 号染色体三体的样品。 如果不釆用本发明的待测染 色体的 DNA片段的 GC含量的校正， 而仅从来自第 21号染色体的 DNA片
段占总 DNA 片段的比值来看， 仅仅能发现在图 9 中来自第 21 号染色体的 DNA片段占总 DNA片段的比值最大的 4个样品是第 21号染色体三体， 不 能发现比值在 0.0138的样品的异常， 对于该样品， 就会出现假阴性的结果。 具体实施方式 需要说明的是， 在不矛盾的情况下， 本申请中的实施例及实施例中的特 征可以相互组合。 下面将参考附图并结合实施例来详细说明本发明。 本发明 的附图及其实施例仅仅用于解释本发明， 并不构成对本发明的任何限制。 名词解释 双端短序列是指紧接着 5，端链接引物的小于 50bp的序列和紧接着 3，端 链接引物的小于 50bp的序列。 优选地， 双端短序列是指紧接着 5，端链接引 物的不大于 36bp的序列和紧接着 3'端链接引物的不大于 36bp的序列。 单端短序列是指紧接着 5，端链接引物的小于 50bp的序列或紧接着 3，端 链接引物的小于 50bp的序列。 优选地， 单端短序列是指紧接着 5，端链接引 物的不大于 36bp的序列或紧接着 3'端链接引物的不大于 36bp的序列。 单端长序列是指紧接着 5，端链接引物的大于 99bp的序列或紧接着 3，端 链接引物的大于 99bp的序列。 双端测序是指分别测试位于序列两端的序列。 单端测序是指测试位于序列一端的序列。
DNA簇指由单一 DNA分子扩增而成的位于一定固体面积的多个 DNA 分子。 在该申请的实施例中， DNA簇指由单一 DNA分子扩增而成的位于 1 平方 ^啟米内的大约 1000个 DNA分子。 乳化 PCR指的是将 PCR ( Polymerase Chain Reaction ) 的反应物， 包括 PCR模板 DNA、PCR引物、 PCR聚合酶和自由碱基置于一个油珠内进行 PCR 反应。 通常情况下， 在一个油珠内， 乳化 PCR的模板只有一个 DNA分子。 GC 含量： 是指在核酸或脱氧核糖核酸中， 鸟嘌呤和胞嘧啶的数量占所 有碱基总数量的比例。 实施例 1： 一种检测胚胎染色体拷贝数的方法 步骤 1 : 取母体血液， 制备血浆。
在本实施例中， 总共抽取了 14个母体血液样品， 样品的代号为： G356、 G397、 G426、 G735、 G756、 G760、 G763、 G770、 G778、 G779、 G780、 G781、 G824和 G825， 均由中南大学湘雅医学院的邬玲仟教授通过羊膜腔穿 刺经染色体核型鉴定为正常的女性胚胎样本 [46 , XX] , 对上述样品检测的数 据用于故标准曲线， 也将上述样品成为标准样品。 血液样品经高速离心后得 到除去了血细包的血浆样品， 每个样品的血浆量大约为 lml。 步骤 2: 提取血浆 DNA。 可以使用 Qiagen公司生产的 DNA抽提试剂盒来提取血浆中的 DNA(产 品号为 57704 )。 步 4聚 3 : 将血浆 DNA制成可供测序的文库。 可以 夺血浆 DNA制成可供双端短序列测序、 单端长序列测序或单端短 序列测序的文库， 以下是制成双端短序列测序文库的过程。 将提取的 DNA末端补平并进行 5，端磷酸化： 将 DNA30 μ1、 纯水 45 μ1、 具有 10mM ATP的 T4 DNA连接酶緩冲液 10 μ1、 10mM dNTP Mix 4 μ1、 T4 DNA 聚合酶 5 μ1、 Klenow酶 1 μ1、 Τ4 ΡΝΚ 5 μΐ混合后， 在 20°C温浴 30分 钟（试剂为 Illumina样本准备试剂盒 PE- 102- 1001提供 )„温浴后釆用 QIAGEN QIAquick PCR 纯化试剂盒 ( part # 28104 ) 纯化 DNA。 末端悬 A: 将上步的产物溶解在 32 μΐ緩冲液中， 加入 Klenow 緩冲液 5 μ1、 ImM dATP 10 μ1、 Klenow Εχο- 3 μΐ, 在 37°C保持 30分钟（试剂为 Illumina 样本准备试剂盒 PE- 102- 1001提供 ), 产物由 QIAGEN MinElute PCR纯化试 剂盒 （ part # 28004 ) 连接： DNA溶解在 ΙΟμΙ緩冲液中， 加入 DNA连接酶緩冲液 2χ 25μ1、 PE Adapter Oligo Mix ΙΟμΙ, DNA 连接酶 5μ1, 在 20°C保持 15分钟 (试剂为 Illumina样本准备试剂盒 PE- 102- 1001提供 )„温浴后釆用 QIAGEN QIAquick PCR纯化试剂盒 ( part # 28104 ) 纯化 DNA。 图 1是样品的 DNA被制作成可供双端测序的文库， 将该文库中的 DNA 用 1%的琼脂糖凝胶进行电泳，图 1是凝胶电泳的图像，最明显的条带在 280bp 左右， 由于包含了 120bp的连接引物， 因此母体血浆 DNA主要片段集中在 160bp左右。 优选地， 还可以对可供双端测序的文库进行 PCR扩增： 1 μ1 ϋΝΑ、 22 μΐ 纯水、 Ιμΐ PE PCR 引物 PE 2.0、 1 μΐ PE PCR引物 PE 1.0、 2x Phusion DNA
聚合酶 ( Finnzymes Oy ) (试剂为 Illumina样本准备试剂盒 PE- 102- 1001提 供）。通过 PCR仪器扩增 DNA,程序为 98°C 30秒； 98°C 40秒， 65 °C 30秒， 72 °C 30秒,共进行 12个循环； 72 °C 5分钟。 步骤 4： 对所述 DNA测序文库进行测序。 才艮据在步骤 3中所制备的文库的不同， 可以分别进行双端短序列测序、 单端长序列测序或单端短序列测序。 以下是进行双端短序列测序的过程。 用 Illumina的 cBot仪器将 DNA双端测序文库中单一的 DNA分子制成 DNA簇， 这一步也可由乳化 PCR将单一 DNA分子变成珠子上的多分子。 ^！夺上面生成的 DNA袭或乳化 PCR得到的 DNA珠在 Illumina 的 Genome Analyzer或 HiSeq2000测序仪进行双端测序。此过程均由仪器本身自动完成。 可替换地， 也可将上面生成的 DNA簇或乳化 PCR得到的 DNA珠在 Illumina的 Genome Analyzer、 HiSeq2000或 Life Technologies公司的 SOLiD 测序仪进行单端长序列的测序。 在 Illumina的 Genome Analyzer, HiSeq2000 和 Life Technologies公司的 SOLiD测序仪上的测序步骤和反应条件与以上描 述的一样。 步骤 5： 确定血浆中的 DNA片段来自几号染色体， 并确定待测染色体的 拷贝数是否正常。 在进行了 DNA文库的双端短序列测序后 （可替换地， 也可以是测定单 端长序列或单端短序列测序）， 当已知了一个 DNA片段两端各 36bp的碱基 序 列 后 ， 可 以 将这两 端 的序 列 与 人类基因 组标准序 列 37. 1 ( http： //www. ncbi . nlm.nih.gov/projects/genome/as sembly/gr c/human/data/?build =37 ), 该数据库也称 hg l9 比对， 确定了这两端的序列分别在染色体上的位 置， 该两端序列之间的距离就是该 DNA 片段的长度， 同时该两端序列所在 的染色***置即确定了该 DNA片段来自几号染色体。 图 2是根据上述的方法确定得到的样品 G356的 DNA片段大小的分布 图， 从图中可以看出， 母体血浆的短 DNA主要集中在 lOObp到 220bp之间。 图 3是将比对后样品 G356中确定来自 X染色体的 DNA片段根据片段 大小所 ^¾丈的分布图。 所得到的图形与图 2几乎一致。 步骤 6: 计算在同一样本内所述 DNA中来自待测染色体的 DNA片段数 量占所有 DNA片段的比值， 根据所述来自待测染色体的 DNA片段中的 GC 含量对所述比值进行校正， 并计算待测样本中所述来自待测染色体的 DNA
片段校正后的比值的变异， 根据所述变异的程度确定待测染色体的拷贝数。 优选地， 取 DNA的序列长度在 100bp-250bp之间的任意一点或任意一 个区间的所有 DNA的测序和比对结果， 计算在同一样本内 100bp-250bp之 间的任意一点或任意一个区间的所有 DNA中来自待测染色体的 DNA片段数 量占 100bp-250bp之间的任意一点或任意一个区间的所有 DNA片段的比值， 才艮据所述来自待测染色体的 DNA片段中的 GC含量对所述比值进行校正， 并计算待测样本中所述来自待测染色体的 DNA 片段校正后的比值的变异， 根据所述变异的程度确定待测染色体的拷贝数。 实验中发现， 取 DNA的序列长度在 100bp-250bp之间的任意一点或任 意一个区间的所有 DNA的测序和比对结果， 能够提高实 -险结果的准确度。 具体的算法如下： 使用标准样品， 计算在同一样本内来自待测染色体的 DNA片段占所有 DNA片段的比值， 根据测序的结果可以知道来自待测染色 体的 DNA片段中的 GC含量。并根据上述比值和 GC含量制定出 GC含量（ Y 轴） 和各个染色体或局部染色体占全部染色体百分比 （X轴） 的标准曲线。 将待测样品所测得的染色体的比值根据该染色体自身特有的函数将其校正到 某一固定的 GC值， 一般常用 GC值的算术平均数。 并计算待测样本中所述 来自待测染色体的 DNA片段校正后的比值的变异 Z值， 根据 Z值确定待测 染色体的拷贝数。 图 5 A-H是根据所测的样品中来自待测染色体的 DNA片段的 GC含量 和来自待测染色体的 DNA片段占总 DNA片段的比值， 所做的标准曲线。 从 图中可以看出， 来自第 2、 3、 4、 5、 6、 7、 8、 12、 13、 14、 18和 X号 染色体的 DNA片段的 GC含量分别与来自各个染色体的 DNA片段占总 DNA 片段的比值成反比， 所述线性关系可用 y=ax+b表示， 其中 y代表来自待测 染色体的 DNA片段的 GC含量， X代表来自待测染色体的 DNA片段数量占 总 DNA的比值， a和 b是常数， JL a是负数。 需要说明的是， 针对不同的参 考样本， 公式的具体参数可能会有微小的变化， 但总体趋势是不变的。 在表 1中示出了 4十对第 2、 3、 4、 5、 6、 7、 8、 12、 13、 14、 18和 X号染色 体的函数： 表 1
染色体编号 函数
2 Y = - ■1274.5X + 154.48
3 Υ = - ■682.68Χ + 88.11
4 Υ = - ■391.42X + 62.075
5 Υ = - ■772.25Χ + 88.517
6 Υ = - ■729.61X + 84.354
7 -2874.7Χ + 197.28
8 -1599.4X + 122.84
12 -1936.7X + 129.32
13 -827.29Χ + 67.049
14 -933.9Χ + 74.22
18 -1946.4X + 97.323
X -749.77Χ + 71.33 来自第 1、 9、 γ γ γ γ γ γ γ 10、 11、 15、 16、 17、 19、 20、 21和 22号染色体的 DNA 片段的 GC含量分别与来自各个染色体的 DNA片段占总 DNA片段的比值成 正比，所述线性关系可用 y = ax +b表示，其中 y代表来自待测染色体的 DNA 片段的 GC含量， X代表来自待测染色体的 DNA片段数量占总 DNA的比值， a和 b是常数， JL a是正数。
表 2
染色体编号 函数
1 Y = 909.17X - 29.372
9 Υ = 2211.3X - 45.632
10 Υ = 1886.8X - 51.026
11 Υ = 1255X - 15.714
15 Υ = 1775.8Χ - 10.124
16 Υ = 797.17X + 23.883
17 Υ = 560.83Χ + 31.227
19 Υ = 513.79X + 43.088
20 Υ = 1006.7X + 22.818
21 Υ = 7298Χ - 51.083
22 Υ = 596.24Χ + 43.346 然后， 将待测样品所测得的染色体的比值根据该染色体自身特有的函数 将其校正到某一固定的 GC值， 一般常用 GC值的算术平均数。 举例来说， 由来自各个染色体的 DNA片段的 GC含量与来自各个染色 体的 DNA片段占总 DNA片段的比值表示的 y=ax+b的函数可知， 每个样本 中来自各个染色体的 DNA片段占总 DNA片段的比值 X可以通过此函数根据
GC值校正。 校正公式可以是： a 。 其中， X为校正后来自各个 染色体的 DNA片段占总 DNA片段的比值； έ为检测到的样本在该指定染色 体片段占总 DNA片段的比值； y为待测样本在指定的染色体片段上的 GC含 量； γ为检测到的标准样本 （已知正常的样本， 如实施例 1 中检测的已知正 常的 14个标准样品）在这一段染色体的算术平均 GC值； a为这条曲线的斜 率。 针对每条染色体的函数见表 1和 2; 针对不同的样本， έ、 y这两个值为 该次试-险的测定值， 这两个参数随着待测样本不同而不同。
通过计算这些校正过的 X 值的平均值和标准误差， 可以得到数值 Z, Z= (待测样品校正过的比值 X-校正过的各标准样品比值的平均值 )/标准误差。 计算各待测样本所述校正过的比值 X的变异， -据所述变异的数值 Z确定待 测染色体或待测局部染色体的拷贝数。 一般情况下， 将 Z的绝对值小于 3视 为正常检测误差， Z的绝对值大于 3视为异常。 以下表 3、 4和 5是通过实施例 1的方法检测样品 G356、 G397(重复测 了 8次）、 G426、 G735、 G756、 G760、 G763、 G770、 G778、 G779、 G780、 G781、 G824和 G825计算得到的针对第 13、 18和 21号染色体的检测以及校 正结果。 表 3 : 第 13号染色体校正示意表
表 4: 第 18号染色体校正示意表
样品 X Y 校正后的 X
G356 0.028355 42.32462 0.028009403
G397L1 0.028042 42.66689 0.027871871
G397L2 0.028211 42.5936 0.028003058
G397L3 0.027866 42.92468 0.027828101
G397L4 0.028163 42.8853 0.028105316
s/uu O z-iz-oosldsoiAV
这里需要特别说明的是上述检测和计算的方法不仅适用于检测整个染色 体拷贝数的异常， 同样适用于检测局部染色体的拷贝数异常。 以下是检测待测样品中胚胎 13号染色体拷贝数的实例 如上面所述的操作， 总共抽取了 15 个母体血液样品， 血液样品经高速 离心后得到除去了血细胞的血浆样品， 每个样品的血浆量大约为 1ml, 样品 的代号为： G352、 G362, G372, G383 , G397(重复测了 8 次）、 G402、 G409, G415 , G424, G445 , G440, G503 , G588, G735和 G783。 上述样 品均由中南大学湘雅医学院的邬玲仟教 ·ί受釆集得到。 对可能存在的 13号染色体三体进行检测。 利用上面得到的针对 13号染 色体的标准曲线， 对所述检测的结果进行-险证。 如图 6所示， 在没有 GC校 正的情况下 （图中以 X轴上的菱形来标记）， 13 号染色体三体样品 G445、 G352和 G402 ( x轴上 3个用圓圏标记的样品）与正常样品（ x轴上不加标记 的样品） 无法区分开来。 但在 GC 校正的情况下 （图中以矩形来表示）， 13 号染色体三体样品 G445、 G352和 G402 (图中 3个用箭头标 ΐ己的样品）与正 常样品 （图中其他以矩形表示的样品） 就可以清楚地区分开来， 而且经 GC 校正的结果与中南大学湘雅医学院的邬玲仟教 ·ί受通过羊膜腔穿刺经染色体核 型鉴定的结果一致。 所以 GC含量的校正可以用来检测 13号染色体三体， 减 少假阴性结果的出现。 下面的表 6中仅示出了部分待测样品的检测和计算结果， 其他样品的结 果未示出。 第 13号染色体校正后 Ζ值计算示意表
从上表可以见到 G445、 G352、 G402样本 Z值超过 3 , 可以判断其为 13
号染色体三体。 其余各样本 Z值均在 -3到 +3之间。 以下是检测待测样品中胚胎 18号染色体拷贝数的实例 如上面所述的操作， 总共抽取了 16 个母体血液样品， 血液样品经高速 离心后得到除去了血细胞的血浆样品， 每个样品的血浆量大约为 1ml, 样品 的代号为： G362, G372, G383 , G397(重复测了 8 次）、 G407, G409, G415 , G424, G442, G432, G445 , G440, G595 , G588, G735和 G783。 上述样品均由中南大学湘雅医学院的邬玲仟教授釆集得到。 对可能存在的 18号染色体三体进行检测。 如图 7所示， 利用上面得到 的针对 18号染色体的标准曲线， 对所述检测的结果进行验证。 在没有 GC校 正的情况下 （图中以 X轴上的菱形表示；)， 有些 18号染色体三体样品 （X轴 上 1个用黑色圓圏标记的样品） 与正常样品 （X轴上以菱形表示但不加标记 的样品） 无法区分开来， 所以导致该样品的^ _阴性。 其他 18 号染色体三体 的样品（图中用向下箭头标记的在 X轴上的样品） 可以在没有 GC含量校正 的情况下检测出来。 但在 GC校正的情况下 （图中以矩形来标记）， 所有 18 号染色体三体样品 （图中用横向箭头标记的样品） 与正常样品 （图中其他以 矩形表示的样品） 就可以清楚地区分开来， 而且经 GC校正的结果与中南大 学湘雅医学院的邬玲仟教授通过羊膜腔穿刺经染色体核型鉴定的结果一致。 所以 GC含量的校正可以用来检测 18号染色体三体，减少假阴性结果的出现。 下面的表 7中仅示出了部分待测样品的检测和计算结果， 其他样品的结 果未示出。 第 18号染色体校正后 Z值计算示意表
从上表可以见到 G424、 G442、 G432、 G415、 G595、 G407样本 Z值超 过 3 , 可以判断其为 18号染色体三体。 其余各样本 Z值均在 -3到 +3之间。 以下是检测待测样品中胚胎 21号染色体拷贝数的实例 如上面所述的操作， 总共抽取了 14 个母体血液样品， 血液样品经高速 离心后得到除去了血细胞的血浆样品， 每个样品的血浆量大约为 1ml, 样品 的代号为： G267, G387, G393 , G376, G397(重复测了 8次）、 G405 , G408, G409, G440, G491 , G588, G641 , G735和 G783。 上述样品均由中南大 学湘雅医学院的邬玲仟教授釆集得到。 对可能存在的 21号染色体三体即唐氏综合征样品进行检测。如图 9A所 示， 利用上面得到的针对 21 号染色体的标准曲线， 对所述检测的结果进行 验证。 在没有 GC校正的情况下 （图中以 X轴上的菱形表示；)， 21号染色体 三体的样品可以单凭它占全部染色体的百分比检测出来 （X轴上 1个用黑色 圓圏和 3个用竖向箭头标记的样品）。 其中一个样品（X轴上 1个用黑色圓圏 标记的样品） 与正常样品 （其他菱形标记的样品） 的区别不是很明显， 可能 导致该样品的假阴性检测。 但在 GC校正的情况下（图中以矩形来标记；)， 所 有 21号染色体三体样品（图中 4个用横向箭头标记的样品）与正常样品（图 中其他以矩形标记的样品） 就可以清楚地区分开来， 而且经 GC校正的结果 与中南大学湘雅医学院的邬玲仟教授通过羊膜腔穿刺经染色体核型鉴定的结 果一致。 GC含量校正对 21号染色体三体检测精确度的提高可以用 21号染 色体三体样品距正常样品之间的最小距离来衡量。 如图 9B 所示， 校正后的 最小距离 d l比校正前的最小距离 d 2要大。所以 GC含量的校正可以用来更 精确地检测 21号染色体三体， 减少假阴性结果的出现。 下面的表 8中仅示出了部分待测样品的检测和计算结果， 其他样品的结 果未示出。 表 8: 第 21号染色体校正后 Z值计算示意表
第 18号染色体 校正后 X
平均值 0.013225973
标准误 0.0000762249 从上表可以见到 G405、 G387、 G376、 G393样本 Z值超过 3 , 可以判断 其为 21号染色体三体。 其余各样本 Z值均在 -3到 +3之间。 实施例 2: 检测胚胎或肿瘤染色体拷贝数的试剂盒 对应于实施例 1的检测方法， 本申请的发明人开发了一种能够用于检测 胚胎或肿瘤染色体拷贝数的试剂盒， 其包括： 从母体或肿瘤患者取血液的器械， 可以是任何能用于取血的釆血针、 注 射器等等； 适合将所述血液中的血细胞和血浆分离的器械， 可以是适合于在离心机 上用于盛装血液的微管或其他任何适合分离的容器或器械； 抽提所述血浆中的 DNA的试剂和器械， 可以包括： 蛋白酶、 饱和酚、 氯仿：异戊醇（24 ： 1 )、 醋酸钠、 无水乙醇、 70%乙醇、 TE溶液等。 也可以 选取 Qiagen公司生产的 DNA抽提试剂盒来提取血浆中的 DNA (产品号为 57704 ) 以及其它任何可以用于进行 DNA提取的试剂或容器； 将所述 DNA制成可供测序的文库的试剂和器械， 供测序的文库可以是 供双端短序列测序的文库、单端长序列测序的文库或单端短序列测序的文库。 将所述 DNA制成可供双端短序列测序的文库的试剂和器械包括：具有 10mM ATP的 T4 DNA连接酶緩冲液、 lOmM dNTP Mix, T4 DNA聚合酶、 Klenow 酶、 T4 PNK (以上这些试剂为 Illumina样本准备试剂盒 PE-102-1001提供 ) 以及在特定环境下对 DNA具有亲和性的离子交换脂以实现对 DNA的分离。 也可以选取 QIAGEN QIAquick PCR 产物分离试剂盒 （产品号 # 28104 )或 QIAGEN MinElute PCR产物分离试剂盒 （产品号 # 28004 )„ 和对所述 DNA进行测序的试剂和器械，对所述 DNA进行测序可以是双 端短序列测序、 单端长序列测序或单端短序列测序。 进行双端短序列测序的 试剂和器才成可以包括： PE PCR引物 PE 2.0、 PE PCR 引物 PE 1.0、 Phusion DNA聚合酶 ( Finnzymes Oy ) (试剂为 Illumina样本准备试剂盒 PE-102-1001 提供）。 实施例 3： 另一种检测胚胎或肿瘤染色体拷贝数的试剂盒 本申请的发明人开发了另一种能够用于检测胚胎或肿瘤染色体拷贝数的 试剂盒， 其包括：
从母体或肿瘤患者取血液的器械， 可以是任何能用于取血的釆血针、 注 射器等等； 适合将所述血液中的血细胞和血浆分离的器械， 可以是适合于在离心机 上用于盛装血液的微管或其他任何适合分离的容器或器械； 抽提所述血浆中的 DNA的试剂和器械， 可以包括： 蛋白酶、 饱和酚、 氯仿：异戊醇（24 ： 1 )、 醋酸钠、 无水乙醇、 70%乙醇、 TE溶液等。 也可以 选取 Qiagen公司生产的 DNA抽提试剂盒来提取血浆中的 DNA (产品号为 57704 ) 以及其它任何可以用于进行 DNA提取的试剂或容器； 通过物理方法根据所述 DNA 的片段大小对其进行分离的试剂和器械， 可以包括： 琼月旨糖粉 （Biowest 11860 ), marker ( Takara lOObp DNA marker, 产品号 D505A ) 等； 和取 100bp-250bp之间的任意一点或任意一个区间的所有 DNA进行测 序的试剂和器械， 可以包括： 如切取一定区间的琼脂糖凝胶的切刀。 优选地， 可以包括将从切取的琼脂糖凝胶回收 DNA进行扩增并将其制 成可供测序的文库的试剂和器才成。 实施例 4: 一种检测胚胎或肿瘤染色体拷贝数的装置
一种检测胚胎或肿瘤染色体拷贝数的装置， 包括： 检测模块， 用于对母体血浆或肿瘤患者血浆样本中 DNA进行测序， 所 述测序包括了将母体血浆或肿瘤患者血浆样本中的 DNA制成可供测序的文 库的过程， 对母体血浆样本中的 DNA进行测序， 可以包括 Illumina的 cBot 仪器和 Illumina的 Genome Analyzer或 HiSeq2000测序仪或 ABI公司的 SOLiD 系列的测序仪； 比对模块， 用于将所述 DNA的测序结果与染色体组序列图谱进行比对， 以确定所述 DNA中的每段 DNA来自几号染色体以及所述每段 DNA的序列 长度， 可以使用人类基因组标准序列数据库 hgl9; 计算模块， 用于计算在同一样本内来自待测染色体的 DNA 片段数量占 所有 DNA片段的比值，根据所述来自待测染色体的 DNA片段中的 GC含量 对所述比值进行校正； 并计算待测样本中所述来自待测染色体的 DNA 片段 校正后的比值的变异， 根据所述变异的程度确定待测染色体的拷贝数； 以及 输出模块， 用于输出所述待测染色体的拷贝数。
可选地， 检测模块可以检测样品中所有的 DNA 片段， 也可以仅仅检测 100bp-250bp之间的任意一点或任意一个区间如 150bp到 175bp的所有 DNA, 这可以通过在检测装置的上游包括一个根据所述 DNA 的片段大小， 对母体 血浆中的 DNA进行分离的试剂和器械， 可以是进行琼脂糖凝胶电泳的模块 或装置。 显然， 本领域的技术人员应该明白， 上述的本发明的一些模块或一些步 骤可以用通用的计算装置来实现， 它们可以集中在单个的计算装置上， 或者 分布在多个计算装置所组成的网络上， 可选地， 它们可以用计算装置可执行 的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行， 或者将它们分别制作成各个集成电路模块， 或者将它们中的多个模块或步骤 制作成单个集成电路模块来实现。 这样， 本发明不限制于任何特定的硬件和 软件结合。 以上仅为本发明的优选实施例而已， 并不用于限制本发明， 对于本领域 的技术人员来说， 本发明可以有各种更改和变化。 凡在本发明的 ^"神和原则 之内， 所作的任何修改、 等同替换、 改进等， 均应包含在本发明的保护范围 之内。
参考文献：
1. Cunningham F, et al. ( 2002 ) In Williams Obstretrics ( McGraw-Hill Professional, New York ) , p942.
2. Wapner R, et al. ( 2003 ) First- trimester screening for trisomies 21 and 18.
N Engl J Med , 3 9 \ 05- \ \3.
3. Lo YM, et al. ( 1997 ) Presence of fetal DAN in maternal plasma and serum. Lancet , 350: 485-487.
4. Lo YM, et al. ( 2007 ) Plasma placental RNA allelic ratio permits
noninvasive prenatal chromosomal aneuploidy detection. Nat Med, 13: 218-223.
5. Tong YK, et al. ( 2006 ) Noninvasive prenatal detection of fetal trisomy 18 by epigenetic allelic ratio analysis in maternal plasma: Theoretical and empirical considerations. Clin Chem, 52: 2194-2202.
6. Fan HC, Quake SR. ( 2007 ) Detection of aneuploidy with digital
polymerase chain reaction. Anal Chem, 79: 7576-7579.
7. Lo YM, et al. ( 2007 ) Digital PCR for molecular detection of fetal
chromosomal aneuploidy. Proc Natl Acad Sci USA, 104: 131 16-13121.
8. Bentley DR, et al. ( 2008 ) Accurate whole human genome sequencing using reversible terminator chemistry. Nature, 456: 53-59.
9. McKeman KJ, et al. ( 2009 ) Sequence and structure variation in a human genome uncovered by short-read, massively parallel ligation sequencing using two-base encoding. Genome Research, 119: 1527-1541.
10. Harris TD, et al. ( 2008 ) Single-molecule DNA sequencing of a viral genome. Science, 320: 106- 109.
11. Fan HC, et al. ( 2008 ) Noninvasive diagnosis of fetal aneuploidy by sequencing DNA from maternal blood. Proc Natl Acad Sci USA, 105:
16266-16271.
12. Chiu RWK, et al. ( 2008 ) Noninvasive prenatal diagnosis of fetal
chromosomal aneuploidy by massively parallel genomic sequencing of DNA in maternal plasma. Proc Natl Acad Sci USA, 105: 20458-20463.
13. Chiu RWK, et al. ( 2010 ) Maternal plasma DNA analysis with massively parallel sequencing by ligation for noninvasive prenatal diagnosis of trisomy 21. Chin Chem, 56:459-463.
14. Lewis DE, et al. ( 2010 ) Antigenic approach to the detection and
isolation of microparticles associated fetal DNA. PCT, US2010, #025209.
Fan HC, et al. ( 2010 ) Analysis of the size distributions of fetal and maternal cell-free DNA by paired-end sequencing. Clin Chem, 56: 1279-1286.
Claims (32)
- 权 利 要 求 书1. 一种检测胚胎或肿瘤染色体拷贝数的试剂盒， 包括：从母体或肿瘤患者体内取血液的器械；适合将所述血液中的血细胞和血浆分离的器械；抽提所述血浆中的 DNA的试剂和器械；通过物理方法根据所述 DNA 的片段大小对其进行分离的试剂和 器械； 和取 100bp-250bp之间的任意一点或任意一个区间的所有 DNA进行 测序的试剂和器械。
- 2. 根据权利要求 1所述的试剂盒， 其中所述胚胎或肿瘤染色体是整条染色 体或局部染色体。
- 3. 根据权利要求 1所述的试剂盒， 还包括： 将所述 DNA制成供测序的文 库的试剂和器械。
- 4. 根据权利要求 2所述的试剂盒， 还包括： 将所述 DNA制成供测序的文 库的试剂和器械。
- 5. 根据权利要求 4所述的试剂盒， 还包括： 对所述供测序的文库进行 PCR 扩增的试剂和器械。
- 6. 根据权利要求 1或 5所述的试剂盒， 其中所述 100-250bp之间的任意一 点或任意一个区间的 DNA是 150bp-170bp的 DNA。
- 7. 根据权利要求 6所述的试剂盒， 其中所述 100-250bp之间的任意一点或 任意一个区间的 DNA是 167bp的 DNA。
- 8. —种检测胚胎或肿瘤染色体拷贝数的试剂盒， 包括：从母体或肿瘤患者体内取血液的器械；适合将所述血液中的血细胞和血浆分离的器械；抽提所述血浆中的 DNA的试剂和器械；将所述 DNA制成可供测序的文库的试剂和器戈； 和对所述 DNA进行测序的试剂和器械。. 根据权利要求 8所述的试剂盒， 其中所述胚胎或肿瘤染色体是整条染色 体或局部染色体。
- 10. 根据权利要求 8 所述的试剂盒， 还包括： 对从血浆中抽提的所述 DNA 进行 PCR扩增的试剂和器械。
- 11. 一种检测胚胎或肿瘤染色体拷贝数的装置， 包括：检测模块，用于对母体血浆或肿瘤患者血浆样本中的 DNA进行测 序；比对模块，用于将所述 DNA的测序结果与染色体组序列图谱进行 比对， 以确定所述 DNA中的每段 DNA来自几号染色体以及所述每段 DNA的序列长度；计算模块，用于计算在同一样本内来自待测染色体的 DNA片段的 数量占 DNA片段总数的比值， 根据所述来自待测染色体的 DNA片段 中的 GC含量对所述比值进行校正； 并计算待测样本中所述来自待测 染色体的 DNA片段校正后的比值的变异，根据所述变异的程度确定待 测染色体的拷贝数； 以及输出模块， 用于输出所述待测染色体的拷贝数。
- 12. 根据权利要求 11所述的装置，其中所述计算模块根据以下的函数对第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18 和 X号染色体的比值进行校正： 来自第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18和 X号染色体的 DNA 片段的 GC含量分别与来自各个染色体的 DNA片段占总 DNA片段的比 值具有一定的线性关系， 所述线性关系可用 y=ax+b表示， 其中 y代表 来自待测染色体的 DNA片段的 GC含量， X代表来自待测染色体的 DNA 片段数量占总 DNA的比值， a和 b是常数， 且 a是负数。
- 13. 根据权利要求 12所述的装置， 其中所述来自第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18和 X号染色体的 DNA片段的 GC含量与来自各个染色 体的 DNA片段占总 DNA片段的比值之间的线性关系如以下所示染色体编号 函数2 Y = -1274.5X + 154.483 Υ = -682.68Χ + 88.114 Y = -391.42X + 62.0755 Υ = -772.25Χ + 88.5176 Υ = -729.61X + 84.3547 Υ = -2874.7Χ + 197.288 Υ = -1599.4X + 122.8412 Υ = -1936.7X + 129.32 Y = -827.29X + 67.049Y = - 1946.4X + 97.323Υ = -749.77Χ + 71.33
- 14. 根据权利要求 11所述的装置，其中所述计算模块根据以下的函数对第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21和 22号染色体的比值进行校正： 来自第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21和 22号染色体的 DNA 片段的 GC含量分别与来自各个染色体的 DNA片段占总 DNA片段的比 值具有一定的线性关系， 所述线性关系可用 y=ax+b表示， 其中 y代表 来自待测染色体的 DNA片段的 GC含量， X代表来自待测染色体的 DNA 片段数量占总 DNA的比值， a和 b是常数， 且 a是正数。
- 15. 才艮据权利要求 14所述的装置， 其中所述来自第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21和 22号染色体的 DNA片段的 GC含量与来自各个染色 体的 DNA片段占总 DNA片段的比值之间的线性关系如在以下所示： 染色体编号 函数1 Y = 909. 17X - 29.3729 Υ = 2211.3X - 45.63210 Υ = 1886.8X - 5 1.02611 Υ = 1255X - 15.71415 Υ = 1775.8Χ - 10. 12416 Υ = 797. 17X + 23.88317 Υ = 560.83Χ + 31.22719 Υ = 513.79X + 43.08820 Υ = 1006.7X + 22.81821 Υ = 7298Χ - 5 1.08322 Υ = 596.24Χ + 43.346
- 16. 根据权利要求 11- 15所述的装置， 其中所述检测模块将母体血浆或肿瘤 患者血浆样本中的所有 DNA制成可供测序的文库的过程。
- 17. 居权利要求 11所述的装置，其中所述对母体血浆或肿瘤患者血浆样本 中的 DNA进行测序是双端短序列测序、 单端长序列测序或单端短序列 测序。
- 18. 根据权利要求 11所述的装置，其中所述胚胎或肿瘤染色体是整条染色体 或局部染色体。
- 19. 根据权利要求 17所述的装置， 其中对所述母体血浆或肿瘤患者血浆样本 中的 DNA在测序前进行了 PCR扩增。
- 20. 根据权利要求 19 所述的装置， 其中计算模块用于计算在同一样本内 100bp-250bp之间的任意一点或任意一个区间的所有 DNA中来自待测染 色体的 DNA片段数量占 100bp-250bp之间的任意一点或任意一个区间 的所有 DNA片段总数的比值， 根据 100bp-250bp之间的任意一点或任 意一个区间的所有 DNA中来自待测染色体的 DNA片段的 GC含量对所 述比值进行校正； 并计算待测样本中所述来自待测染色体的 DNA 片段 校正后的比值的变异， 根据所述变异的程度确定待测染色体的拷贝数。
- 21. 居权利要求 20 所述的装置， 其中将所述母体血浆或肿瘤患者血浆样 本中 100bp-250bp之间的任意一点或任意一个区间的所有 DNA制成可 供测序的文库。
- 22. 居权利要求 21 所述的装置， 其中所述母体血浆或肿瘤患者血浆样本 中的 DNA在被制作成可供测序的文库之前或之后进行了 PCR扩增。
- 23. 根据权利要求 20所述的装置，其中所述 100-250bp之间的任意一点或任 意一个区间的 DNA是 150bp-170bp的 DNA。
- 24. 根据权利要求 23所述的装置，其中所述 100-250bp之间的任意一点或任 意一个区间是 DNA是 167bp的 DNA。
- 25. 一种检测胚胎或肿瘤染色体拷贝数的方法， 包括以下步骤：取母体血浆或肿瘤患者血浆；将所述血液中的血浆和血细月包分离；将所述血浆中的 DNA制作成可供测序的文库；对所述 DNA测序文库进行测序；将所述测序的结果与染色体组序列图谱进行比对， 确定每段所述 DNA序列来自几号染色体以及每段所述 DNA序列的长度； 以及用所述 DNA的测序和比对结果，计算在同一样本内来自待测染色 体的 DNA片段的数量占 DNA片段总数的比值， 根据所述来自待测染 色体的 DNA片段中的 GC含量对所述比值进行校正，并计算待测样本 中所述来自待测染色体的 DNA片段校正后的比值的变异，根据所述变 异的程度确定待测染色体的拷贝数。
- 26. 根据权利要求 25中任一项所述的方法， 根据以下的函数对第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18和 X号染色体的比值进行校正： 来自第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18和 X号染色体的 DNA片段的 GC 含量分别与来自各个染色体的 DNA片段占总 DNA片段的比值具有一定 的线性关系， 所述线性关系可用 y = ax + b表示， 其中 y代表来自待测 染色体的 DNA片段的 GC含量， X代表来自待测染色体的 DNA片段数 量占总 DNA的比值， a和 b是常数， 且 a是负数。
- 27. 根据权利要求 26所述的方法， 其中所述来自第 2、 3、 4、 5、 6、 7、 8、 12、 13、 18和 X号染色体的 DNA片段的 GC含量与来自各个染色 体的 DNA片段占总 DNA片段的比值之间的线性关系如在权利要求 13 中所示。
- 28. 根据权利要求 25中任一项所述的方法， 根据以下的函数对第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21和 22号染色体的比值进行校正： 来自第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21和 22号染色体的 DNA片段的 GC含量分别与来自各个染色体的 DNA片段占总 DNA片段的比值具有 一定的线性关系， 所述线性关系可用 y=ax+b表示， 其中 y代表来自待 测染色体的 DNA片段的 GC含量， X代表来自待测染色体的 DNA片段 数量占总 DNA的比值， a和 b是常数， JL a是正数。
- 29. 根据权利要求 28所述的方法， 其中所述来自第 1、 9、 10、 11、 15、 16、 17、 19、 20、 21和 22号染色体的 DNA片段的 GC含量与来自各个染色 体的 DNA片段占总 DNA片段的比值之间的线性关系如在权利要求 15 中所示。
- 30. 根据权利要求 25-29 所述的方法， 其中仅取 DNA 的序列长度在 100bp-250bp之间的任意一点或任意一个区间的所有 DNA的测序和比对 结果， 计算在同一样本内 100bp-250bp之间的任意一点或任意一个区间 的所有 DNA中来自待测染色体的 DNA片段数量占 100bp-250bp之间的 任意一点或任意一个区间的所有 DNA 片段总数的比值， 根据 100bp-250bp之间的任意一点或任意一个区间的所有 DNA中所述来自待 测染色体的 DNA片段中的 GC含量对所述比值进行校正， 并计算待测 样本中所述来自待测染色体的 DNA 片段校正后的比值的变异， 根据所 述变异的程度确定待测染色体的拷贝数。
- 31. 根据权利要求 25-29所述的方法， 其中对所述 DNA测序文库进行测序 是双端短序列测序、 单端长序列测序或单端短序列测序。
- 32. 根据权利要求 30所述的方法，其中所述 100-250bp之间的任意一点或任 意一个区间的 DNA是 150bp-170bp的 DNA。
- 33. 根据权利要求 32所述的方法，其中所述 100-250bp之间的任意一点或任 意一个区间的 DNA是 167bp的 DNA。
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/CN2011/075037 WO2012162884A1 (zh) | 2011-05-31 | 2011-05-31 | 检测胚胎或肿瘤染色体拷贝数的试剂盒、装置和方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103080336A true CN103080336A (zh) | 2013-05-01 |
| CN103080336B CN103080336B (zh) | 2014-06-04 |
Family
ID=47258277
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201180035510.7A Active CN103080336B (zh) | 2011-05-31 | 2011-05-31 | 检测胚胎或肿瘤染色体拷贝数的试剂盒、装置和方法 |
Country Status (11)
| Country | Link |
|---|---|
| US (1) | US9885080B2 (zh) |
| EP (1) | EP2716766B1 (zh) |
| JP (1) | JP5926795B2 (zh) |
| CN (1) | CN103080336B (zh) |
| DK (1) | DK2716766T3 (zh) |
| ES (1) | ES2605372T3 (zh) |
| HU (1) | HUE031239T2 (zh) |
| PL (1) | PL2716766T3 (zh) |
| PT (1) | PT2716766T (zh) |
| RS (1) | RS55518B1 (zh) |
| WO (1) | WO2012162884A1 (zh) |
Families Citing this family (42)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012525147A (ja) | 2009-04-30 | 2012-10-22 | グッド スタート ジェネティクス， インコーポレイテッド | 遺伝マーカーを評価するための方法および組成物 |
| US9163281B2 (en) | 2010-12-23 | 2015-10-20 | Good Start Genetics, Inc. | Methods for maintaining the integrity and identification of a nucleic acid template in a multiplex sequencing reaction |
| WO2014014498A1 (en) * | 2012-07-20 | 2014-01-23 | Verinata Health, Inc. | Detecting and classifying copy number variation in a fetal genome |
| WO2012177792A2 (en) | 2011-06-24 | 2012-12-27 | Sequenom, Inc. | Methods and processes for non-invasive assessment of a genetic variation |
| WO2013052907A2 (en) | 2011-10-06 | 2013-04-11 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
| US10196681B2 (en) | 2011-10-06 | 2019-02-05 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
| US9367663B2 (en) | 2011-10-06 | 2016-06-14 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
| US9984198B2 (en) | 2011-10-06 | 2018-05-29 | Sequenom, Inc. | Reducing sequence read count error in assessment of complex genetic variations |
| US10424394B2 (en) | 2011-10-06 | 2019-09-24 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
| EP2768983A4 (en) | 2011-10-17 | 2015-06-03 | Good Start Genetics Inc | METHODS OF IDENTIFYING MUTATIONS ASSOCIATED WITH DISEASES |
| EP4148739A1 (en) | 2012-01-20 | 2023-03-15 | Sequenom, Inc. | Diagnostic processes that factor experimental conditions |
| US8209130B1 (en) | 2012-04-04 | 2012-06-26 | Good Start Genetics, Inc. | Sequence assembly |
| US10227635B2 (en) | 2012-04-16 | 2019-03-12 | Molecular Loop Biosolutions, Llc | Capture reactions |
| US9920361B2 (en) | 2012-05-21 | 2018-03-20 | Sequenom, Inc. | Methods and compositions for analyzing nucleic acid |
| US10504613B2 (en) | 2012-12-20 | 2019-12-10 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
| US10497461B2 (en) | 2012-06-22 | 2019-12-03 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
| US10482994B2 (en) | 2012-10-04 | 2019-11-19 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
| US20130309666A1 (en) | 2013-01-25 | 2013-11-21 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
| WO2014152421A1 (en) | 2013-03-14 | 2014-09-25 | Good Start Genetics, Inc. | Methods for analyzing nucleic acids |
| FI2981921T3 (fi) | 2013-04-03 | 2023-03-09 | Sequenom Inc | Menetelmiä ja prosesseja geneettisten variaatioiden ei-invasiiviseen arviointiin |
| KR102665592B1 (ko) | 2013-05-24 | 2024-05-21 | 시쿼넘, 인코포레이티드 | 유전적 변이의 비침습 평가를 위한 방법 및 프로세스 |
| AU2014284180B2 (en) | 2013-06-21 | 2020-03-19 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
| MY181069A (en) | 2013-10-04 | 2020-12-17 | Sequenom Inc | Methods and processes for non-invasive assessment of genetic variations |
| US10438691B2 (en) | 2013-10-07 | 2019-10-08 | Sequenom, Inc. | Non-invasive assessment of chromosome alterations using change in subsequence mappability |
| US10851414B2 (en) | 2013-10-18 | 2020-12-01 | Good Start Genetics, Inc. | Methods for determining carrier status |
| AU2014340239B2 (en) | 2013-10-21 | 2019-11-28 | Verinata Health, Inc. | Method for improving the sensitivity of detection in determining copy number variations |
| WO2015175530A1 (en) * | 2014-05-12 | 2015-11-19 | Gore Athurva | Methods for detecting aneuploidy |
| KR102566176B1 (ko) | 2014-05-30 | 2023-08-10 | 베리나타 헬스, 인코포레이티드 | 태아 아-염색체 홀배수체 및 복제수 변이 탐지 |
| US20150347676A1 (en) * | 2014-05-30 | 2015-12-03 | Sequenom, Inc. | Chromosome representation determinations |
| WO2016019042A1 (en) | 2014-07-30 | 2016-02-04 | Sequenom, Inc. | Methods and processes for non-invasive assessment of genetic variations |
| WO2016040446A1 (en) | 2014-09-10 | 2016-03-17 | Good Start Genetics, Inc. | Methods for selectively suppressing non-target sequences |
| WO2016048829A1 (en) | 2014-09-24 | 2016-03-31 | Good Start Genetics, Inc. | Process control for increased robustness of genetic assays |
| US11072814B2 (en) | 2014-12-12 | 2021-07-27 | Verinata Health, Inc. | Using cell-free DNA fragment size to determine copy number variations |
| EP3271480B8 (en) | 2015-01-06 | 2022-09-28 | Molecular Loop Biosciences, Inc. | Screening for structural variants |
| US10095831B2 (en) | 2016-02-03 | 2018-10-09 | Verinata Health, Inc. | Using cell-free DNA fragment size to determine copy number variations |
| US11200963B2 (en) | 2016-07-27 | 2021-12-14 | Sequenom, Inc. | Genetic copy number alteration classifications |
| JP7237003B2 (ja) | 2017-01-24 | 2023-03-10 | セクエノム， インコーポレイテッド | 遺伝子片の評価のための方法およびプロセス |
| CN107217308A (zh) * | 2017-06-21 | 2017-09-29 | 北京贝瑞和康生物技术股份有限公司 | 一种用于检测染色体拷贝数变异的测序文库构建方法和试剂盒 |
| CN109112191A (zh) * | 2017-06-23 | 2019-01-01 | 沈阳精准医疗技术有限公司 | 一种肿瘤ctDNA的信息统计方法 |
| CN109136363A (zh) * | 2018-08-31 | 2019-01-04 | 领航基因科技（杭州）有限公司 | 一种检测胎儿染色体非整数倍的试剂盒和检测胎儿甲基化dna拷贝数的方法 |
| KR102532991B1 (ko) * | 2019-12-23 | 2023-05-18 | 주식회사 랩 지노믹스 | 태아의 염색체 이수성 검출방법 |
| CN113096728B (zh) * | 2021-06-10 | 2021-08-20 | 臻和(北京)生物科技有限公司 | 一种微小残余病灶的检测方法、装置、存储介质及设备 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102108406A (zh) * | 2010-12-20 | 2011-06-29 | 北京贝瑞和康生物技术有限公司 | 检测胚胎染色体拷贝数的试剂盒、装置和方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE60014399T2 (de) * | 1999-05-04 | 2005-12-08 | Ortho-Clinical Diagnostics, Inc. | Schnelle und wirksame dns-immobilisierung aus einer probe ohne die verwendung von zelllysereagenzien |
| US20100216153A1 (en) | 2004-02-27 | 2010-08-26 | Helicos Biosciences Corporation | Methods for detecting fetal nucleic acids and diagnosing fetal abnormalities |
| US20100112590A1 (en) | 2007-07-23 | 2010-05-06 | The Chinese University Of Hong Kong | Diagnosing Fetal Chromosomal Aneuploidy Using Genomic Sequencing With Enrichment |
| CA3069082C (en) * | 2008-09-20 | 2022-03-22 | The Board Of Trustees Of The Leland Stanford Junior University | Noninvasive diagnosis of fetal aneuploidy by sequencing |
| EP2494065B1 (en) | 2009-10-26 | 2015-12-23 | Lifecodexx AG | Means and methods for non-invasive diagnosis of chromosomal aneuploidy |
| CN102018406A (zh) | 2010-11-30 | 2011-04-20 | 赖新莲 | 一种治疗颈椎病的药枕 |
-
2011
- 2011-05-31 PL PL11866914T patent/PL2716766T3/pl unknown
- 2011-05-31 PT PT118669142T patent/PT2716766T/pt unknown
- 2011-05-31 RS RS20161036A patent/RS55518B1/sr unknown
- 2011-05-31 CN CN201180035510.7A patent/CN103080336B/zh active Active
- 2011-05-31 HU HUE11866914A patent/HUE031239T2/en unknown
- 2011-05-31 EP EP11866914.2A patent/EP2716766B1/en active Active
- 2011-05-31 WO PCT/CN2011/075037 patent/WO2012162884A1/zh active Application Filing
- 2011-05-31 DK DK11866914.2T patent/DK2716766T3/da active
- 2011-05-31 JP JP2014513025A patent/JP5926795B2/ja active Active
- 2011-05-31 ES ES11866914.2T patent/ES2605372T3/es active Active
- 2011-05-31 US US13/811,634 patent/US9885080B2/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102108406A (zh) * | 2010-12-20 | 2011-06-29 | 北京贝瑞和康生物技术有限公司 | 检测胚胎染色体拷贝数的试剂盒、装置和方法 |
Non-Patent Citations (2)
| Title |
|---|
| H. CHRISTINA FAN: "《Analysis of the Size Distributions of Fetal and Maternal Cell-Free DNA by Paired-End Sequencing》", 《CLINICAL CHEMISTRY》 * |
| ROSSA W.K.: "《Maternal Plasma DNA Analysis with Massively Parallel Sequencing by Ligation for Noninvasive Prenatal Diagnosis of Trisomy 21》", 《CLINICAL CHEMISTRY》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| EP2716766A1 (en) | 2014-04-09 |
| PT2716766T (pt) | 2016-11-21 |
| WO2012162884A1 (zh) | 2012-12-06 |
| HUE031239T2 (en) | 2017-07-28 |
| JP5926795B2 (ja) | 2016-05-25 |
| EP2716766A4 (en) | 2014-11-26 |
| JP2014516541A (ja) | 2014-07-17 |
| ES2605372T3 (es) | 2017-03-14 |
| US20130130921A1 (en) | 2013-05-23 |
| PL2716766T3 (pl) | 2017-09-29 |
| RS55518B1 (sr) | 2017-05-31 |
| EP2716766B1 (en) | 2016-09-28 |
| US9885080B2 (en) | 2018-02-06 |
| CN103080336B (zh) | 2014-06-04 |
| DK2716766T3 (da) | 2017-01-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103080336A (zh) | 检测胚胎或肿瘤染色体拷贝数的试剂盒、装置和方法 | |
| AU2020202153B2 (en) | Single-molecule sequencing of plasma DNA | |
| AU2015289414B2 (en) | Non-invasive prenatal diagnosis of fetal genetic condition using cellular DNA and cell free DNA | |
| CN107077537B (zh) | 用短读测序数据检测重复扩增 | |
| JP6534191B2 (ja) | コピー数変動を決定することにおける検出の感度を向上させるための方法 | |
| CN103298955A (zh) | 用于构建血浆dna测序文库的方法和试剂盒 | |
| CN107475403A (zh) | 从外周血游离dna中检测循环肿瘤dna的方法、试剂盒及其测序结果的分析方法 | |
| CN102108406A (zh) | 检测胚胎染色体拷贝数的试剂盒、装置和方法 | |
| TW201639967A (zh) | 檢測胎兒基因資訊的方法、試劑盒、裝置和系統 | |
| CN104694654A (zh) | 一种检测胎儿染色体数目变异的方法 | |
| WO2015196752A1 (en) | A method and a kit for quickly constructing a plasma dna sequencing library | |
| CN108315320A (zh) | 一种rna高通量测序文库的制备方法 | |
| CN107583052A (zh) | miR‑6734‑5p在制备Luminal型乳腺癌诊断工具中的应用 | |
| CN113755592A (zh) | 一种染色质外环状dna的检测方法 | |
| CN111690737A (zh) | 一种a-地中海贫血基因检测试剂盒及使用方法 | |
| WO2019009431A1 (ja) | 腫瘍細胞で生じた突然変異を高精度に識別する方法 | |
| CN108315824A (zh) | 一种rna高通量测序文库的构建方法 | |
| Dennis et al. | Discovery of Cell-Free Fetal DNA in Maternal Blood and Development of Noninvasive Prenatal Testing 2022 Lasker-DeBakey Clinical Medical Research Award |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20130501 Assignee: Auspicious and the health gene diagnosis Technology Co., Ltd. of Hangzhou shellfish Assignor: Beijing Berry Genomics Co., Ltd. Contract record no.: 2015990000958 Denomination of invention: Kits, devices and methods for detecting chromosome copy number of embryo or tumor Granted publication date: 20140604 License type: Common License Record date: 20151116 |
|
| LICC | Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model | ||
| C56 | Change in the name or address of the patentee | ||
| CP03 | Change of name, title or address |
Address after: 102206, No. 5, building 4, 8 life Garden Road, Changping District science and Technology Park, Beijing, 801 Patentee after: Beijing Berry Genomics Biotechnology Co., Ltd Address before: Room 2, floor 207, building C, Shuguang building, No. 5, Shun Tong Road, Chaoyang District, Beijing, China Patentee before: Beijing Berry Genomics Co., Ltd. |
|
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20130501 Assignee: Chengdu Berry and Kang medical laboratory Co., Ltd. Assignor: Beijing Berry Genomics Biotechnology Co., Ltd Contract record no.: 2016990000383 Denomination of invention: Kits, devices and methods for detecting chromosome copy number of embryo or tumor Granted publication date: 20140604 License type: Common License Record date: 20160907 |
|
| LICC | Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 102200 Beijing Changping District Science and Technology Park Life Park Road, No. 4 Courtyard, No. 5 Building, 8th Floor, 801 Patentee after: Beijing Berry Genomics Co., Ltd. Address before: 102206 8th Floor 801, No. 5 Building, No. 4 Life Garden Road, Changping District Science Park, Beijing Patentee before: Beijing Berry Genomics Biotechnology Co., Ltd |