JP5926795B2 - 胎芽又は腫瘍染色体のコピー数を検出するシステム及び装置 - Google Patents

胎芽又は腫瘍染色体のコピー数を検出するシステム及び装置 Download PDF

Info

Publication number
JP5926795B2
JP5926795B2 JP2014513025A JP2014513025A JP5926795B2 JP 5926795 B2 JP5926795 B2 JP 5926795B2 JP 2014513025 A JP2014513025 A JP 2014513025A JP 2014513025 A JP2014513025 A JP 2014513025A JP 5926795 B2 JP5926795 B2 JP 5926795B2
Authority
JP
Japan
Prior art keywords
chromosome
dna
detection target
sample
sequencing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2014513025A
Other languages
English (en)
Other versions
JP2014516541A (ja
Inventor
ガオ,ヤン
チャオ,ダイシン
チャン,ジャングアン
ティエン,ファン
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Berry Genomics Co Ltd
Original Assignee
Berry Genomics Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Berry Genomics Co Ltd filed Critical Berry Genomics Co Ltd
Publication of JP2014516541A publication Critical patent/JP2014516541A/ja
Application granted granted Critical
Publication of JP5926795B2 publication Critical patent/JP5926795B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing
    • C12Q1/6874Methods for sequencing involving nucleic acid arrays, e.g. sequencing by hybridisation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • G16B25/10Gene or protein expression profiling; Expression-ratio estimation or normalisation
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6809Methods for determination or identification of nucleic acids involving differential detection
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00ICT 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)

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(10)を含み、一回で数億ないし数十億個の配列を検出できる。これらの技術で母体血漿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断片を占める割合のみに基づいて判定する方法では検出できなかった複数の仮陰性結果を検出できる。具体的な実施形態において具体的な実験で証明している。
一方、例えば、文献に開示されたように(15)、母体血漿中の胎芽DNAの大部分は100bp〜250bpの断片であって、且つ150bp〜170bpのDNAが特に多い。母体のDNAの極少ない部分が当該断片領域に分布されているが、断片が250bpを超えるDNAは基本的に母体のDNAに所属される。本発明者は以下のことを発現した:理由不明であるが、各染色体からのDNA断片が総DNA断片を占める割合は、100bp〜250bp間のいずれかの一点又はいずれかの一つの区間のDNAにおいて均一に分布され、即ち、各染色体が100bp〜250bp間のいずれかの一点、例えば110bp又は167bp(当該点のDNA量が最も多い)において、総DNAを占める割合が他の点の割合を示し、従って、各染色体が100bp〜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はヌクレオソームの形式で存在するので、大部分が100bp〜250bpの断片であって、また、各染色体からのDNA断片が総DNA断片を占める割合が100bp〜250bp間のいずれかの一点又はいずれかの一つの区間のDNAは均一に分布され、即ち、各染色体が100bp〜250bp間のいずれか一点において、総DNAを占める割合は他の点の割合を示し、従って、各染色体が100bp〜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断片を占める割合と一定の線形関係、即ち,即ち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の配列長さが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であることがさらに好ましい。
母体血漿のDNAを両端シークエンシング用のライブラリに作成した後、1%のアガロース・ゲルで電気泳動した後の画像である。左側の通路は、DNA100bp標準画像で、右側の通路は両端シークエンシング用のライブラリDNAの画像で、最も明らかなストリップは280bp前後に位置し、その中、120bpの接続プライマーを含む。 サンプルG356のDNA断片サイズの分布図である。 比較した後にサンプルG356で確定したX染色体からのDNA断片を断片サイズに応じて分布した図である。 計算して得たG356と実際に検出して得たG356のX染色体のDNA断片を断片サイズに応じて分布した図で、図において、星付きの線は、G356サンプルの総DNA配列数にX染色体の百分比を掛けて得た図形を示し、円は実際に検出したX染色体の配列数を示す。図に示すように、計算して得たG356X染色体のDNA断片の分布と実際に検出して得た分布は大体同じで、2つの線は大体合致する。 サンプル番号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含有量を示す。5A〜Fはそれぞれ、三つの染色体の標準曲線を示し、順に、第1〜18番染色体の標準曲線を示す図で、図5Gは第19〜22番染色体の標準曲線を示す図で、図5HはX染色体の標準曲線を示す図である。各染色体の標準曲線の関数は、表1と表2に示すとおりである。 サンプル番号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含有量を示す。5A〜Fはそれぞれ、三つの染色体の標準曲線を示し、順に、第1〜18番染色体の標準曲線を示す図で、図5Gは第19〜22番染色体の標準曲線を示す図で、図5HはX染色体の標準曲線を示す図である。各染色体の標準曲線の関数は、表1と表2に示すとおりである。 サンプル番号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含有量を示す。5A〜Fはそれぞれ、三つの染色体の標準曲線を示し、順に、第1〜18番染色体の標準曲線を示す図で、図5Gは第19〜22番染色体の標準曲線を示す図で、図5HはX染色体の標準曲線を示す図である。各染色体の標準曲線の関数は、表1と表2に示すとおりである。 サンプル番号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含有量を示す。5A〜Fはそれぞれ、三つの染色体の標準曲線を示し、順に、第1〜18番染色体の標準曲線を示す図で、図5Gは第19〜22番染色体の標準曲線を示す図で、図5HはX染色体の標準曲線を示す図である。各染色体の標準曲線の関数は、表1と表2に示すとおりである。 サンプル番号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含有量を示す。5A〜Fはそれぞれ、三つの染色体の標準曲線を示し、順に、第1〜18番染色体の標準曲線を示す図で、図5Gは第19〜22番染色体の標準曲線を示す図で、図5HはX染色体の標準曲線を示す図である。各染色体の標準曲線の関数は、表1と表2に示すとおりである。 サンプル番号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含有量を示す。5A〜Fはそれぞれ、三つの染色体の標準曲線を示し、順に、第1〜18番染色体の標準曲線を示す図で、図5Gは第19〜22番染色体の標準曲線を示す図で、図5HはX染色体の標準曲線を示す図である。各染色体の標準曲線の関数は、表1と表2に示すとおりである。 サンプル番号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含有量を示す。5A〜Fはそれぞれ、三つの染色体の標準曲線を示し、順に、第1〜18番染色体の標準曲線を示す図で、図5Gは第19〜22番染色体の標準曲線を示す図で、図5HはX染色体の標準曲線を示す図である。各染色体の標準曲線の関数は、表1と表2に示すとおりである。 サンプル番号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含有量を示す。5A〜Fはそれぞれ、三つの染色体の標準曲線を示し、順に、第1〜18番染色体の標準曲線を示す図で、図5Gは第19〜22番染色体の標準曲線を示す図で、図5HはX染色体の標準曲線を示す図である。各染色体の標準曲線の関数は、表1と表2に示すとおりである。 本発明による方法で検出した第13番染色体3倍体のサンプルを示す図である。ここで、X軸は13番染色体からのDNA断片の数量が総DNA断片を占める割合を示し、Y軸は13番染色体からのDNA断片のGC含有量を示す。X軸において、菱形で各サンプル中の第13番染色体からのDNA断片の数量が総DNA断片を占める割合を示し、矢印で示す矩形は、第13番染色体からのDNA断片のGC含有量を補正した後に検出した第13番染色体3倍体のサンプルを示す。本発明による検出対象染色体のDNA断片のGC含有量の補正を行わず、第13番染色体からのDNA断片が総DNA断片を占める割合から見ると、矢印で示す三つサンプルの異常を発現することができなく、そうなると、仮陰性の結果が現れる。 本発明による方法で検出した第18番染色体3倍体のサンプルを示す図である。なかでは、X軸は18番染色体からのDNA断片の数量が総DNA断片を占める割合を示し、Y軸は18番染色体からのDNA断片のGC含有量を示す。X軸で菱形で各サンプル中の第18番染色体からのDNA断片の数量が総DNA断片を占める割合を示し、矢印で示す矩形は第18番染色体からのDNA断片のGC含有量を補正した後に検出した第18番染色体3倍体のサンプルを示す。本発明に係わる検出対象染色体のDNA断片のGC含有量の補正を行わず、第18番染色体からのDNA断片が総DNA断片を占める割合のみから見ると、第18番染色体からのDNA断片が総DNA断片を占める割合が0.031以上である五つのサンプルが第18番染色体3倍体であることのみを発現でき、割合が0.03前後である2つのサンプルの異常は発現できなく、当該2つのサンプルの場合、仮陰性の結果が現れる。 本発明による方法で検出したX染色体単体又は3倍体サンプルを示す図である。そのなか、X軸はX染色体からのDNA断片の数量が総DNA断片を占める割合を示し、Y軸はX染色体からのDNA断片のGC含有量を示す。X軸で菱形で各サンプル中のX染色体からのDNA断片の数量が総DNA断片を占める割合を示し、矢印で示す標準曲線左側の矩形はX染色体からのDNA断片のGC含有量を補正した後に検出したX染色体単体のサンプルを示し、矢印で示す標準曲線右側の矩形はX染色体からのDNA断片のGC含有量を補正した後に検出したX染色体3倍体のサンプルを示す。本発明に係わる検出対象染色体のDNA断片のGC含有量の補正を行わず、X染色体からのDNA断片が総DNA断片を占める割合のみから見ると、X染色体3倍体の一つのサンプルのみ発現でき、X染色体単体のサンプルを発現しにくく、そうなると、仮陰性の結果が現れる。 本発明による方法で検出した第21番染色体3倍体のサンプルを示す図である。その中、X軸は21番染色体からのDNA断片の数量が総DNA断片を占める割合を示し、Y軸は21番染色体からのDNA断片のGC含有量を示す。X軸で菱形で各サンプル中の第21番染色体からのDNA断片の数量が総DNA断片を占める割合を示し、矢印で示す矩形で第21番染色体からのDNA断片のGC含有量を補正した後に検出した第21番染色体3倍体のサンプルを示す。本発明による検出対象染色体のDNA断片のGC含有量の補正を行わず、第21番染色体からのDNA断片が総DNA断片を占める割合のみから見ると、図9中の第21番染色体からのDNA断片が総DNA断片を占める割合が最大の四つのサンプルが第21番染色体3倍体であることのみを発現でき、割合が0.0138であるサンプルの異常は発現できなく、当該サンプルに仮陰性の結果が現れる。 本発明による方法で検出した第21番染色体3倍体のサンプルを示す図である。その中、X軸は21番染色体からのDNA断片の数量が総DNA断片を占める割合を示し、Y軸は21番染色体からのDNA断片のGC含有量を示す。X軸で菱形で各サンプル中の第21番染色体からのDNA断片の数量が総DNA断片を占める割合を示し、矢印で示す矩形で第21番染色体からのDNA断片のGC含有量を補正した後に検出した第21番染色体3倍体のサンプルを示す。本発明による検出対象染色体のDNA断片のGC含有量の補正を行わず、第21番染色体からのDNA断片が総DNA断片を占める割合のみから見ると、図9中の第21番染色体からのDNA断片が総DNA断片を占める割合が最大の四つのサンプルが第21番染色体3倍体であることのみを発現でき、割合が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鋳型DNA、PCRプライマー、PCRポリメラーゼ、自由塩基を含むPCR(Polymerase Chain Reaction)の反応物を一つのオイル・ドロップ内でPCR反応を行うことを指す。通常、一つのオイル・ドロップ内で、PCRの乳状化を行う鋳型は一つのDNA分子だけである。
GC含有量は、核酸又はデオキシリボ核酸において、グアニンとシトシンの数量が塩基総量で占める比例を指す。
実施例1:胎芽染色体コピー数の検出方法
工程1:母体血液を採取して血漿を調製する。
本実施例において、合計で14個の母体血液サンプルを取得し、サンプルの番号は、G356、G397、G426、G735、G756、G760、G763、G770、G778、G779、G780、G781、G824、G825であって、いずれも中南大学湘雅医学院の▲ウ▼玲仟教授が羊水穿刺によって染色体核型鑑定を経て正常の女性の胎芽サンプル[46,XX]であると鑑定されたものである。上記サンプルを検出したデータは標準曲線の作成に用いられ、また、上記サンプルは標準サンプルになる。血液サンプルを高速遠心処理して血細胞が除去された血漿サンプルを取得し、各サンプルの血漿量は約1mlである。
工程2:血漿DNAを抽出する。
Qiagen会社が生産したDNA抽出試薬カートリッジを用いて血漿中のDNA(製品番号57704)を抽出することができる。
工程3:血漿DNAをシークエンシング用のライブラリに作成する。
血漿DNAを双端短配列シークエンシング、単端長配列シークエンシング又は単端短配列シークエンシング用のライブラリに作成することができ、以下は双端短配列シークエンシング用のライブラリに作成する工程である。
抽出したDNA末端を充満して5’端のリン酸化を行って、DNA30μl、純水45μl、10mM ATPを有するT4 DNAリガーゼ緩衝液10μl、10mM dNTP Mix 4μl、T4 DNA ポリメラーゼ5μl、Klenow酵素1μl、T4 PNK5μl混合した後、20℃で、30分の温浴を行う(試薬は、Illuminaサンプル準備試薬カートリッジPE〜102〜1001が提供)。温浴後にQIAGEN QIAquick PCR浄化試薬カートリッジ(part # 28104)を用いてDNAを浄化した。
末端懸A:前工程で得られた製品を32μlの緩衝液に溶解し、Klenow緩衝液5μl、1mM dATP 10μl、Klenow Exo〜3μlを添加し、在37℃で30分保持し(試薬は、Illuminaサンプル準備試薬カートリッジPE〜102〜1001が提供)、製品はQIAGEN MinElute PCR浄化試薬カートリッジ(part # 28004)
接続:DNAを10μlの緩衝液に溶解し、DNA リガーゼ緩衝液2X25μl、PE Adapter Olige Mix 10μl,DNAリガーゼ5μlを添加し、20℃で15分保持する(試薬は、Illuminaサンプル準備試薬カートリッジPE〜102〜1001が提供)。温浴後にQIAGEN QIAquick PCR浄化試薬カートリッジ(part # 28104)を用いてDNAを浄化する。
図1は、サンプルのDNAが両端シークエンシング用のライブラリに作成されたもので、当該ライブラリ中のDNAは1%のアガロース・ゲルで電気泳動を行い、図1はアガロース・ゲル電気泳動の画像で、最も明確なストリップは280bp前後にある。120bpの接続プライマーを含むので、母体血漿DNAの主な断片は160bp前後に集中されている。
両端シークエンシング用のライブラリに対してPCR増幅を行うことが好ましい。つまり、1μlのDNA、22μlの純水、1μlのPE PCR プライマーPE2.0、1μlのPE PCRプライマーPE1.0、2x Phusion DNAポリメラーゼ(Finnzymes Oy)(試薬はIlluminaサンプル準備試薬カートリッジPE〜102〜1001が提供)。PCR機器でDNA増幅を行って、その工程は、98℃で30秒、98℃で40秒、65℃で30秒、72℃で30秒で、合計12個の循環を行って、72℃で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/assembly/grc/human/data/?build=37)(当該データベースはHG19とも呼ばれる)と比較し、当該両端の配列各々の染色体上の位置を確定した。当該両端配列間の距離が当該DNA断片の長さであって、同時に、当該両端配列が位置する染色体の位置によって当該DNA断片が何番の染色体からのものであるかを確定できる。
図2は、上記の方法で確定して得たサンプルG356のDNA断片サイズの分布図で、図2に示すように、母体血漿の短DNAは主に100bp〜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値に基づいて検出対象染色体のコピー数を確定する。
図5A〜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は常数で、且つaは負数である。ここで、異なる参照サンプルの場合、式中の具体的なパラメーターが若干の変化があるが、総傾向は変化しない。表1に第2、3、4、5、6、7、8、12、13、14、18、X番染色体の関数を示した。
第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は正数である。
そして、検出対象サンプルを検出して得た染色体の割合を、当該染色体特有の関数に従って、ある固定のGC値に補正し、通常、GC値の算術平均を用いる。
例えば、各染色体からのDNA断片のGC含有量と各染色体からのDNA断片が総DNA断片を占める割合を示すy=ax+bの関数から分かるように、各サンプル中の各染色体からのDNA断片が総DNA断片を占める割合xを、当該関数によってGC値で補正することができる。補正式は、
式1
である。ここで、
は検出したサンプルが当該所定の染色体断片が総DNA断片を占める割合を示し、yは検出対象サンプルの所定の染色体断片上のGC含有量を示し、
は検出した標準サンプル(既に正常のサンプルであることを把握、例えば、実施例1において検出した既に正常であることを把握した14個の標準サンプル)の当該断片の染色体における算術平均GC値を示し、aは当該曲線の傾斜度を示す。各染色体の関数は表1と表2に示す通りで、異なるサンプルのy、yの値は今回試験の検出値であって、当該2つのパラメーターは検出対象サンプルが異なると相違する。
このような補正済みの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番の染色体の検出及び補正結果である。
ここで、上記検出と計算の方法は染色体全体のコピー数の異常の検出に適すると共に、一部の染色体のコピー数の異常の検出にも適する。
以下は、検出対象サンプル中の胎芽13番の染色体のコピー数を検出する実施例である。
上述のように、合計で15個の母体血液サンプルを採取し、血液サンプルは高速遠心処理を経て血細胞が除去された血漿サンプルであって、各サンプルの血漿量は約1mlで、サンプルの番号はG352、G362、G372、G383、G397(重複して8回検出した)、G402、G409、G415、G424、G445、G440、G503、G588、G735、G783である。上記サンプルはいずれも中南大学湘雅医学院の▲ウ▼玲仟教授によって収集したものである。
存在する可能性のある13番染色体3倍体を検出する。上述のように得られた13番の染色体の標準曲線を用いて、上述の検出の結果を検証する。図6に示すように、GC補正を行っていない場合(図においてX軸上の菱形で標記)、13番染色体3倍体のサンプルG445、G352、G402(X軸上三つの円形で標記したサンプル)と正常のサンプル(X軸上で標記していないサンプル)とを区別することができない。しかし、GC補正を行った場合(図において矩形で示す)、13番染色体3倍体のサンプルG445、G352、G402(図において三つの矢印で標記したサンプル)と正常のサンプル(図において他の矩形で示すサンプル)とを明確に区別することができ、また、GC補正を行った結果は中南大学湘雅医学院の▲ウ▼玲仟教授が羊水穿刺して染色体核型鑑定した結果と一致する。従って、GC含有量の補正は、13番の染色体3倍体の検出に用いられ、仮陰性結果の出現を低減できる。
以下の表6に一部の検出対象サンプルの検出と計算結果のみを示す、他のサンプルの結果は示していない。
表に示すように、G445、G352、G402サンプルのZ値が3を超えると、13番の染色体3倍体と確定することができる。他の各サンプルのZ値はいずれも−3〜+3の間である。
以下は、検出対象サンプル中の胎芽18番染色体のコピー数を検出する実施例である。
上述のように、合計16個の母体血液サンプルを採取し、血液サンプルは高速遠心処理を経て血細胞が除去された血漿サンプルで、各サンプルの血漿量は約1mlで、サンプルの番号はG362、G372、G383、G397(重複して8回検出した)、G407、G409、G415、G424、G442、G432、G445、G440、G595、G588、G735、G783である。上記サンプルはいずれも中南大学湘雅医学院の▲ウ▼玲仟教授が収集したものである。
存在する可能性のある18番の染色体3倍体を検出する。図7に示すように、上述のように得られた18番の染色体の標準曲線を用いて、前記検出の結果を検証する。GC補正を行っていない場合(図においてX軸上の菱形で示す)、一部の18番染色体3倍体のサンプル(X軸上の黒の円形で標記したサンプル)と正常のサンプル(X軸上で菱形で示すが標記していないサンプル)とを区別することができないので、当該サンプルの仮陰性をもたらした。他の18番染色体3倍体のサンプル(図において下向き矢印で標記したX軸上のサンプル)はGC含有量の補正を行わずに検出することができる。しかし、GC補正を行った場合(図において矩形で標記)、全ての18番染色体3倍体のサンプル(図において横向き矢印で標記したサンプル)と正常のサンプル(図において他の矩形で示すサンプル)とを明確に区別することができ、また、GC補正を行った結果は中南大学湘雅医学院の▲ウ▼玲仟教授が羊水穿刺して染色体核型鑑定した結果と一致する。従って、GC含有量の補正は、18番染色体3倍体の検出に用いられ、仮陰性結果の出現を低減できる。
以下の表7に一部の検出対象サンプルの検出と計算結果のみを示し、他のサンプルの結果は示していない。
表に示すように、G424、G442、G432、G415、G595、G407サンプルのZ値が3を超えると、18番染色体3倍体であると確定することができる。他の各サンプルのZ値はいずれも−3〜+3間である。
以下は、検出対象サンプル中の胎芽21番染色体のコピー数を検出する実施例である。
上述のように、合計で14個の母体血液サンプルを採取し、血液サンプルは高速遠心処理を経て血細胞が除去された血漿サンプルであって、各サンプルの血漿量は約1mlで、サンプルの番号はG267、G387、G393、G376、G397(重複して8回検出した)、G405、G408、G409、G440、G491、G588、G641、G735、G783である。上記サンプルはいずれも中南大学湘雅医学院の▲ウ▼玲仟教授が収集したものである。
21番染色体3倍体が存在する可能性のあるダウン症候群のサンプルを検出する。図9Aに示すように、上述のように得られた21番染色体の標準曲線を用いて、前記検出の結果を検証する。GC補正を行っていない場合(図においてX軸上の菱形で示す)、21番染色体3倍体のサンプルを、総染色体を占める百分比のみに基づいて検出することができる(X軸上の黒の円形と三つの縦方向矢印で標記するサンプル)。その中の一つのサンプル(X軸上の黒の円形で標記したサンプル)と正常のサンプル(他の菱形で標記したサンプル)との区別が明確でなく、当該サンプルの仮陰性が検出される恐れがある。しかし、GC補正を行った場合(図において矩形で標記)、全ての21番染色体3倍体のサンプル(図において四つの横方向の矢印で標記したサンプル)と正常のサンプル(図において他の矩形で標記するサンプル)とを明確に区別でき、また、GC補正を行った結果は中南大学湘雅医学院の▲ウ▼玲仟教授が羊水穿刺して染色体核型鑑定した結果に一致する。GC含有量の補正の、21番染色体3倍体の検出精度の向上に対する影響を、21番染色体3倍体のサンプルと正常のサンプルとの最小距離に基づいて判断することができる。図9Bに示すように、補正後の最小距離d1は補正前の最小距離d2より大きい。従って、GC含有量の補正によると、21番染色体3倍体をさらに正確に検出し、仮陰性結果の出現を低減することができる。
以下の表8には一部の検出対象サンプルの検出と計算結果のみを示し、他のサンプルの結果は示していない。
表に示すように、G405、G387、G376、G393サンプルのZ値が3を超えると、21番染色体3倍体であると確定できる。他の各サンプルのZ値はいずれも:−3〜+3の間である。
実施例2:胎芽又は腫瘍染色体のコピー数を検出する試薬カートリッジ
実施例1の検出方法に対応し、本願の発明者は胎芽又は腫瘍染色体のコピー数を検出できる試薬カートリッジを開発した。当該試薬カートリッジは、血液の採取に用いることのできる全ての採血針、注射器等であることができる、母体又は腫瘍患者から採血する器具と、遠心分離機に適して血液を収容するための微小管又は他の分離に適する全ての容器又は器具であることができる、前記血液中の血細胞と血漿との分離に適する器具と、プロテアーゼ、飽和フェノール、クロロホルム:イソペンチルアルコール(24:1)、酢酸ナトリウム、無水アルコール、70%アルコール、TE溶液等を含むことができる前記血漿中のDNAを抽出する試薬及び器具と、ここで、Qiagen会社が生産したDNA抽出試薬カートリッジを用いて血漿中のDNA(製品番号57704)を抽出することができ、また他のDNAの抽出に適する試薬又は容器を選択することができ、前記DNAを、双端短配列シークエンシングのライブラリ、単端長配列シークエンシングのライブラリ又は単端短配列シークエンシングのライブラリであることのできるシークエンシング用のライブラリに作成する試薬及び器具と、ここで、前記DNAを双端短配列シークエンシング用のライブラリに作成する試薬及び器具は、10mM ATPを有するT4 DNAリガーゼ緩衝液と、10mM dNTP Mixと、T4 DNAポリメラーゼと、Klenow酵素と、T4PNK(以上の試薬はIlluminaサンプル準備試薬カートリッジPE―102―1001が提供)と、特定の環境でDNAと親和性を有するイオン交換樹脂と、を含み、DNAを分離させる。また、QIAGENQIAquick PCR製品である分離試薬カートリッジ(製品番号#28104)又はQIAGEN MinElute PCR製品である分離試薬カートリッジ(製品番号#28004)を選択することもでき、前記DNAに対して、双端短配列シークエンシング、単端長配列シークエンシング又は単端短配列シークエンシングであることのできるシークエンシングを行う試薬及び器具と、を含み、ここで、双端短配列シークエンシングを行う試薬及び器具は、PEPCRプライマーPE2.0と、PE PCRプライマーPE1.0と、Phusion DNAポリメラーゼ(FinnzymesOy)と、を含む(試薬はIlluminaサンプル準備試薬カートリッジPE―102―1001が提供)。
実施例3:他の胎芽又は腫瘍染色体のコピー数を検出する試薬カートリッジ
本願の発明者は、他の胎芽又は腫瘍染色体のコピー数を検出できる試薬カートリッジを開発し、それは、血液の採取に用いることのできる全ての採血針、注射器等であることができる、母体又は腫瘍患者から採血する器具と、遠心分離機に適して血液を収容するための微小管又は他の分離に適する全ての容器又は器具であることができる、前記血液中の血細胞と血漿との分離に適する器具と、プロテアーゼ、飽和フェノール、クロロホルム:イソペンチルアルコール(24:1)、酢酸ナトリウム、無水アルコール、70%アルコール、TE溶液等を含むことができる前記血漿中のDNAを抽出する試薬及び器具と、ここで、Qiagen会社が生産したDNA抽出試薬カートリッジを用いて血漿中のDNA(製品番号57704)を抽出することができ、また他のDNAの抽出に適する試薬又は容器を選択することができ、アガロース粉(Biowest 11860)、markermarker(Takara 100bp DNA marker、製品番号D505A)等を含むことができる、物理方法で前記DNAの断片をそのサイズに応じて分離する試薬及び器具と、例えば、一定の区間のアガロース・ゲルを取るカッターを含むことができる、100bp〜250bp間のいずれかの一点又はいずれかの一つの区間の全てのDNAを選択してシークエンシングを行う試薬及び器具と、を含む。
切り取ったアガロース・ゲルからDNAを回収して増幅を行ってシークエンシング用のライブラリに作成する試薬及び器具を含むことが好ましい。
実施例4:胎芽又は腫瘍染色体のコピー数を検出する装置
胎芽又は腫瘍染色体のコピー数を検出する装置は、母体の血漿サンプル又は腫瘍患者の血漿サンプル中のDNAに対してシークエンシングを行う検出モジュールと、ここで、前記シークエンシングは母体の血漿サンプル又は腫瘍患者の血漿サンプル中のDNAをシークエンシング用のライブラリに作成する工程であり;母体血漿のサンプル中のDNAに対するシークエンシングに、IlluminaのcBot機器とIlluminaのGenome Analyzer又はHiSeq2000シークエンシング器又はABI会社のSOLiDシリーズのシークエンシング器を用いることができ、前記DNAのシークエンシング結果と染色体組配列図譜とを比較して前記DNA中の各DNA断片が何番の染色体からのものであるか、及び前記各DNA断片の配列長さを確定する比較モジュールと、ここで、人類ゲノム標準配列データベースhg19を利用することができ、同一のサンプル中の検出対象染色体からの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 fortrisomies 21 and 18. N Engl J Med, 349: 1405-1413.
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 formolecular detection of fetal chromosomal aneuploidy. Proc Natl Acad Sci USA, 104: 13116-13121.
8. Bentley DR, et al. (2008) Accurate whole human genome sequencing using reversible terminator chemistry. Nature, 456: 53-59.
9. McKernan KJ, et al. (2009) Sequence and structure variation in a human genome uncovered by short-read, massively paralle lligation 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.
15. 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 (24)

  1. 胎芽又は腫瘍染色体のコピー数を検出する試薬カートリッジと胎芽又は腫瘍染色体のコピー数を検出する装置とからなる胎芽又は腫瘍染色体のコピー数を検出するシステムであって、
    前記胎芽又は腫瘍染色体のコピー数を検出する試薬カートリッジは、
    母体又は腫瘍患者体内から血液を採取する器具と、
    前記血液中の血細胞と血漿を分離させるに適する器具と、
    前記血漿中のDNAを抽出する試薬及び器具と、
    物理方法で、前記DNAの断片のサイズに応じてそれを分離する試薬及び器具と、
    100bp〜250bp間のいずれかの一点又はいずれかの一つの区間の全てのDNAを選択してシークエンシングを行う試薬及び器具と、を含み、
    前記腫瘍染色体のコピー数を検出する装置は、
    前記DNAのシークエンシング結果と染色体組配列図譜とを比較して、前記DNA中の各DNA断片が何番の染色体からのものであるか、及び前記各DNA断片の配列長さを確定する比較モジュールと、
    同一のサンプル中の検出対象染色体からのDNA断片の数量がDNA断片の総数を占める割合を計算し、前記検出対象染色体からのDNA断片のGC含有量に基づいて前記割合を補正し、また、下記の式に従って、検出対象サンプルにおける前記検出対象染色体からのDNA断片の補正後の割合の変異を計算し、前記変異レベルに基づいて検出対象染色体のコピー数を確定する計算モジュールと、
    前記検出対象染色体のコピー数を出力する出力モジュールと、を含み、
    補正後の割合の変異=検出対象サンプルの補正済みの割合−補正済みの各標準サンプル割合の平均値)/標準誤差
    検出対象染色体からのDNA断片のGC含有量がそれぞれ、各染色体からのDNA断片が総DNA断片を占める割合と一定の線形関係、即ちy=ax+bの線形関係を有し、前記計算モジュールが、当該関数に従って前記割合を補正し、但し、yは検出対象染色体からのDNA断片のGC含有量を示し、xは検出対象染色体からのDNA断片の数量が総DNAを占める割合を示し、aとbは常数である胎芽又は腫瘍染色体のコピー数を検出するシステム。
  2. 前記胎芽又は腫瘍染色体が全本の染色体又は一部の染色体である請求項1に記載のシステム
  3. 前記DNAをシークエンシング用のライブラリに作成する試薬及び器具をさらに含む請求項1に記載のシステム
  4. 前記DNAをシークエンシング用のライブラリに作成する試薬及び器具をさらに含む請求項2に記載のシステム
  5. 前記シークエンシング用のライブラリに対してPCR増幅を行う試薬及び器具をさらに含む請求項4に記載のシステム
  6. 前記100〜250bp間のいずれかの一点又はいずれかの一つの区間のDNAが150bp〜170bpのDNAである請求項1又は5に記載のシステム
  7. 前記100〜250bp間のいずれかの一点又はいずれかの一つの区間のDNAが167bpのDNAである請求項6に記載のシステム
  8. 胎芽又は腫瘍染色体のコピー数を検出する試薬カートリッジと胎芽又は腫瘍染色体のコピー数を検出する装置とからなる胎芽又は腫瘍染色体のコピー数を検出するシステムであって、
    前記胎芽又は腫瘍染色体のコピー数を検出する試薬カートリッジは、
    母体又は腫瘍患者体内から血液を採取する器具と、
    前記血液中の血細胞と血漿を分離させるに適する器具と、
    前記血漿中のDNAを抽出する試薬及び器具と、
    前記DNAをシークエンシング用のライブラリに作成するための試薬及び器具と、
    前記DNAに対してシークエンシングを行う試薬及び器具とを、含み
    前記腫瘍染色体のコピー数を検出する装置は、
    前記DNAのシークエンシング結果と染色体組配列図譜とを比較して、前記DNA中の各DNA断片が何番の染色体からのものであるか、及び前記各DNA断片の配列長さを確定する比較モジュールと、
    同一のサンプル中の検出対象染色体からのDNA断片の数量がDNA断片の総数を占める割合を計算し、前記検出対象染色体からのDNA断片のGC含有量に基づいて前記割合を補正し、また、下記の式に従って、検出対象サンプルにおける前記検出対象染色体からのDNA断片の補正後の割合の変異を計算し、前記変異レベルに基づいて検出対象染色体のコピー数を確定する計算モジュールと、
    前記検出対象染色体のコピー数を出力する出力モジュールと、を含み、
    補正後の割合の変異=検出対象サンプルの補正済みの割合−補正済みの各標準サンプル割合の平均値)/標準誤差
    検出対象染色体からのDNA断片のGC含有量がそれぞれ、各染色体からのDNA断片が総DNA断片を占める割合と一定の線形関係、即ちy=ax+bの線形関係を有し、前記計算モジュールが、当該関数に従って前記割合を補正し、但し、yは検出対象染色体からのDNA断片のGC含有量を示し、xは検出対象染色体からのDNA断片の数量が総DNAを占める割合を示し、aとbは常数である胎芽又は腫瘍染色体のコピー数を検出するシステム。
  9. 前記胎芽又は腫瘍染色体が全本の染色体又は一部の染色体である請求項8に記載のシステム
  10. 血漿から抽出した前記DNAに対してPCR増幅を行う試薬及び器具をさらに含む請求項8に記載のシステム
  11. 母体の血漿サンプル又は腫瘍患者の血漿サンプル中のDNAに対してシークエンシングを行う検出モジュールと、
    前記DNAのシークエンシング結果と染色体組配列図譜とを比較して、前記DNA中の各DNA断片が何番の染色体からのものであるか、及び前記各DNA断片の配列長さを確定する比較モジュールと、
    同一のサンプル中の検出対象染色体からのDNA断片の数量がDNA断片の総数を占める割合を計算し、前記検出対象染色体からのDNA断片のGC含有量に基づいて前記割合を補正し、そして、下記の式に従って検出対象サンプルにおける前記検出対象染色体からのDNA断片の補正後の割合の変異を計算し、前記変異レベルに基づいて検出対象染色体のコピー数を確定する計算モジュールと、
    前記検出対象染色体のコピー数を出力する出力モジュールと、を含み、
    補正後の割合の変異=検出対象サンプルの補正済みの割合−補正済みの各標準サンプル割合の平均値)/標準誤差
    検出対象染色体からのDNA断片のGC含有量がそれぞれ、各染色体からのDNA断片が総DNA断片を占める割合と一定の線形関係、即ちy=ax+bの線形関係を有し、前記計算モジュールが、当該関数に従って前記割合を補正し、但し、yは検出対象染色体からのDNA断片のGC含有量を示し、xは検出対象染色体からのDNA断片の数量が総DNAを占める割合を示し、aとbは常数である胎芽又は腫瘍染色体のコピー数を検出する装置。
  12. 前記計算モジュールが、第2、3、4、5、6、7、8、12、13、18、X番の染色体の割合を補正し、且つaは負数である請求項11に記載の装置。
  13. 前記第2、3、4、5、6、7、8、12、13、18、X番の染色体からのDNA断片のGC含有量と各染色体からのDNA断片が総DNA断片を占める割合との間の線形関係が以下のとおりである請求項12に記載の装置。
  14. 前記計算モジュールが第1、9、10、11、15、16、17、19、20、21、22番の染色体の割合を補正し、且つaは正数である請求項11に記載の装置。
  15. 前記第1、9、10、11、15、16、17、19、20、21、22番の染色体からのDNA断片のGC含有量と各染色体からのDNA断片が総DNA断片を占める割合との間の線形関係が以下の通りである請求項14に記載の装置。
  16. 前記検出モジュールが、母体の血漿サンプル又は腫瘍患者の血漿サンプル中の全てのDNAをシークエンシング用のライブラリに作成する請求項11乃至15の中のいずれかに記載の装置。
  17. 前記母体の血漿サンプル又は腫瘍患者の血漿サンプル中のDNAに対するシークエンシングが、双端短配列シークエンシング、単端長配列シークエンシング又は単端短配列シークエンシングである請求項11に記載の装置。
  18. 前記胎芽又は腫瘍染色体が全本の染色体又は一部の染色体である請求項11に記載の装置。
  19. 前記母体の血漿サンプル又は腫瘍患者の血漿サンプル中のDNAに対して、シークエンシングを行う前にPCR増幅を行う請求項17に記載の装置。
  20. 計算モジュールが、同一のサンプル中の100bp〜250bp間のいずれかの一点又はいずれかの一つの区間の全てのDNAの中の、検出対象染色体からのDNA断片の数量が100bp〜250bp間のいずれかの一点又はいずれかの一つの区間の全てのDNA断片の総数を占める割合を計算し、100bp〜250bp間のいずれかの一点又はいずれかの一つの区間の全てのDNA中の検出対象染色体からのDNA断片のGC含有量に基づいて前記割合を補正し、そして、検出対象サンプルにおける前記検出対象染色体からのDNA断片の補正後の割合の変異を計算し、前記変異レベルに基づいて検出対象染色体のコピー数を確定する請求項19に記載の装置。
  21. 前記母体の血漿サンプル又は腫瘍患者の血漿サンプル中の100bp〜250bp間のいずれかの一点又はいずれかの一つの区間の全てのDNAをシークエンシング用のライブラリに作成する請求項20に記載の装置。
  22. 前記母体の血漿サンプル又は腫瘍患者の血漿サンプル中のDNAを、シークエンシング用のライブラリに作成される前又はその後にPCR増幅を行う請求項21に記載の装置。
  23. 前記100〜250bp間のいずれかの一点又はいずれかの一つの区間のDNAが150bp〜170bpのDNAである請求項20に記載の装置。
  24. 前記100〜250bp間のいずれかの一点又はいずれかの一つの区間のDNAが167bpのDNAである請求項23に記載の装置。
JP2014513025A 2011-05-31 2011-05-31 胎芽又は腫瘍染色体のコピー数を検出するシステム及び装置 Active JP5926795B2 (ja)

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
JP2014516541A JP2014516541A (ja) 2014-07-17
JP5926795B2 true JP5926795B2 (ja) 2016-05-25

Family

ID=47258277

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014513025A Active JP5926795B2 (ja) 2011-05-31 2011-05-31 胎芽又は腫瘍染色体のコピー数を検出するシステム及び装置

Country Status (11)

Country Link
US (1) US9885080B2 (ja)
EP (1) EP2716766B1 (ja)
JP (1) JP5926795B2 (ja)
CN (1) CN103080336B (ja)
DK (1) DK2716766T3 (ja)
ES (1) ES2605372T3 (ja)
HU (1) HUE031239T2 (ja)
PL (1) PL2716766T3 (ja)
PT (1) PT2716766T (ja)
RS (1) RS55518B1 (ja)
WO (1) WO2012162884A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102532991B1 (ko) * 2019-12-23 2023-05-18 주식회사 랩 지노믹스 태아의 염색체 이수성 검출방법

Families Citing this family (41)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2425240A4 (en) 2009-04-30 2012-12-12 Good Start Genetics Inc METHOD AND COMPOSITION FOR EVALUATING GENETIC MARKERS
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
US20140235474A1 (en) 2011-06-24 2014-08-21 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
US10424394B2 (en) 2011-10-06 2019-09-24 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
US10196681B2 (en) 2011-10-06 2019-02-05 Sequenom, Inc. Methods and processes for non-invasive assessment of genetic variations
WO2013058907A1 (en) 2011-10-17 2013-04-25 Good Start Genetics, Inc. Analysis methods
ES2929923T3 (es) 2012-01-20 2022-12-02 Sequenom Inc Procesos de diagnóstico que condicionan las condiciones experimentales
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
US10504613B2 (en) 2012-12-20 2019-12-10 Sequenom, Inc. Methods and processes for non-invasive assessment of genetic variations
US9920361B2 (en) 2012-05-21 2018-03-20 Sequenom, Inc. Methods and compositions for analyzing nucleic acid
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
US8778609B1 (en) 2013-03-14 2014-07-15 Good Start Genetics, Inc. Methods for analyzing nucleic acids
HUE061261T2 (hu) 2013-04-03 2023-05-28 Sequenom Inc Eljárások és folyamatok genetikai variánsok nem invazív értékelésére
WO2014190286A2 (en) 2013-05-24 2014-11-27 Sequenom, Inc. Methods and processes for non-invasive assessment of genetic variations
LT3011051T (lt) 2013-06-21 2019-05-10 Sequenom, Inc. Genetinių variacijų neinvazinis vertinimo būdas
KR102514024B1 (ko) 2013-10-04 2023-03-23 시쿼넘, 인코포레이티드 유전적 변이의 비침습 평가를 위한 방법 및 프로세스
WO2015054080A1 (en) 2013-10-07 2015-04-16 Sequenom, Inc. Methods and processes for non-invasive assessment of chromosome alterations
US10851414B2 (en) 2013-10-18 2020-12-01 Good Start Genetics, Inc. Methods for determining carrier status
WO2015061359A1 (en) 2013-10-21 2015-04-30 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
AU2015266665C1 (en) 2014-05-30 2021-12-23 Verinata Health, Inc. Detecting fetal sub-chromosomal aneuploidies and copy number variations
EP3149640B1 (en) * 2014-05-30 2019-09-04 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
US10429399B2 (en) 2014-09-24 2019-10-01 Good Start Genetics, Inc. Process control for increased robustness of genetic assays
AU2015360298B2 (en) 2014-12-12 2018-06-07 Verinata Health, Inc. Using cell-free DNA fragment size to determine copy number variations
WO2016112073A1 (en) 2015-01-06 2016-07-14 Good Start Genetics, 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
CA3207879A1 (en) 2017-01-24 2018-08-02 Sequenom, Inc. Methods and processes for assessment of genetic variations
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拷贝数的方法
CN113096728B (zh) * 2021-06-10 2021-08-20 臻和(北京)生物科技有限公司 一种微小残余病灶的检测方法、装置、存储介质及设备

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ES2228520T3 (es) * 1999-05-04 2005-04-16 Ortho-Clinical Diagnostics, Inc. Capatura rapida y eficaz de adn de una muestra sin usar reactivo de lisis celular.
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
PT2562268T (pt) * 2008-09-20 2017-03-29 Univ Leland Stanford Junior Diagnóstico não invasivo de aneuploidia fetal por sequenciação
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 赖新莲 一种治疗颈椎病的药枕
CN102108406B (zh) * 2010-12-20 2014-05-21 北京贝瑞和康生物技术有限公司 检测胚胎染色体拷贝数的试剂盒、装置和方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102532991B1 (ko) * 2019-12-23 2023-05-18 주식회사 랩 지노믹스 태아의 염색체 이수성 검출방법

Also Published As

Publication number Publication date
EP2716766A4 (en) 2014-11-26
JP2014516541A (ja) 2014-07-17
RS55518B1 (sr) 2017-05-31
PL2716766T3 (pl) 2017-09-29
PT2716766T (pt) 2016-11-21
DK2716766T3 (da) 2017-01-02
WO2012162884A1 (zh) 2012-12-06
CN103080336B (zh) 2014-06-04
CN103080336A (zh) 2013-05-01
ES2605372T3 (es) 2017-03-14
EP2716766B1 (en) 2016-09-28
EP2716766A1 (en) 2014-04-09
US20130130921A1 (en) 2013-05-23
US9885080B2 (en) 2018-02-06
HUE031239T2 (en) 2017-07-28

Similar Documents

Publication Publication Date Title
JP5926795B2 (ja) 胎芽又は腫瘍染色体のコピー数を検出するシステム及び装置
ES2945311T3 (es) Detección rápida de aneuploidía
CN103003447B (zh) 用于确定样品中存在或不存在不同非整倍性的方法
KR102605749B1 (ko) 자간전증에 특이적인 순환 rna 시그니처
Ray et al. Essential interpretations of bioinformatics in COVID-19 pandemic
CN102108406A (zh) 检测胚胎染色体拷贝数的试剂盒、装置和方法
CN108728522A (zh) 药物基因组检测方法
Welkers et al. Improved detection of artifactual viral minority variants in high-throughput sequencing data
CN110628898B (zh) Baz1b易感snp位点检测试剂及其制备的试剂盒
KR101907650B1 (ko) 비침습적 태아 염색체 이수성 판별 방법
Link et al. HIV-Quasipore: a suite of HIV-1-specific nanopore basecallers designed to enhance viral quasispecies detection
Hung et al. Genetic diversity and structural complexity of the killer-cell immunoglobulin-like receptor gene complex: A comprehensive analysis using human pangenome assemblies
WO2019009431A1 (ja) 腫瘍細胞で生じた突然変異を高精度に識別する方法
KR20170036649A (ko) 비침습적 태아 염색체 이수성 판별 방법
Morin et al. Genetic and epigenetic links to asthma
CN115948533B (zh) 一种用于检测mybpc3突变基因的试剂及其应用
CN109280697A (zh) 利用孕妇血浆游离dna进行胎儿基因型鉴定的方法
CN110628897B (zh) 一种kfs致病基因新突变及其应用
CN107312864B (zh) 一种用于ppd辅助诊断的wisp3基因突变及其应用
WO2020119626A1 (zh) 一种无创产前检测胎儿是否患有遗传疾病的方法
CN114134228A (zh) 评估PI3K/Art/mTOR通路相关基因突变的试剂盒、系统、储存介质及其应用
Nazario-Toole et al. Genomic Surveillance of SARS-CoV-2
CN115786492A (zh) 用于检测与肺动脉高压有关的新突变基因acvrl1的试剂及其应用
CN116042801A (zh) 用于检测肥厚型心肌病突变基因tnni3、actc1和mypn的试剂及其应用
CN116004806A (zh) 一种非综合征型骨性ⅲ类错颌畸形相关的低频/罕见突变及其检测方法

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150331

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150630

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20151104

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160303

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20160316

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160405

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160422

R150 Certificate of patent or registration of utility model

Ref document number: 5926795

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250