JP2017506908A - 対象シークエンスを有する核酸の量の減少または増加の検出方法 - Google Patents

対象シークエンスを有する核酸の量の減少または増加の検出方法 Download PDF

Info

Publication number
JP2017506908A
JP2017506908A JP2016556745A JP2016556745A JP2017506908A JP 2017506908 A JP2017506908 A JP 2017506908A JP 2016556745 A JP2016556745 A JP 2016556745A JP 2016556745 A JP2016556745 A JP 2016556745A JP 2017506908 A JP2017506908 A JP 2017506908A
Authority
JP
Japan
Prior art keywords
nucleic acid
dna
sequence
sample
amount
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP2016556745A
Other languages
English (en)
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.)
Oxford Gene Technology (operations) Ltd
Oxford Gene Technology Operations Ltd
Original Assignee
Oxford Gene Technology (operations) Ltd
Oxford Gene Technology Operations 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 Oxford Gene Technology (operations) Ltd, Oxford Gene Technology Operations Ltd filed Critical Oxford Gene Technology (operations) Ltd
Publication of JP2017506908A publication Critical patent/JP2017506908A/ja
Pending legal-status Critical Current

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/6806Preparing nucleic acids for analysis, e.g. for polymerase chain reaction [PCR] assay
    • 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/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • 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/6844Nucleic acid amplification reactions
    • C12Q1/6851Quantitative amplification
    • 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/6844Nucleic acid amplification reactions
    • C12Q1/6853Nucleic acid amplification reactions using modified primers or templates
    • C12Q1/6855Ligating adaptors
    • 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/6844Nucleic acid amplification reactions
    • C12Q1/6858Allele-specific amplification
    • 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/6844Nucleic acid amplification reactions
    • C12Q1/686Polymerase chain reaction [PCR]
    • 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/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • C12Q1/6886Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material for cancer
    • 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
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/156Polymorphic or mutational markers

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Organic Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Genetics & Genomics (AREA)
  • Immunology (AREA)
  • Microbiology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Biotechnology (AREA)
  • General Engineering & Computer Science (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Pathology (AREA)
  • Hospice & Palliative Care (AREA)
  • Oncology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

循環セルフリーDNAのサンプル内の胎児あるいは腫瘍DNAのコピー数変化を検知する方法であって、長いDNA分子の増幅が差別れるようにDNAを増幅する方法。これは胎児あるいは腫瘍DNAの選択的な増幅をもたし、たとえばアレイ比較ゲノムハイブリダイゼーションを使用して解析される。

Description

本発明は、正常なサンプルと比較して、サンプル中の対象シークエンスを有する核酸の量の減少または増加を決定する方法に関する。特には、核酸は胎児の核酸あるいは腫瘍に由来してもよい。
核酸の分析は多くの分野の中で使用される;たとえば、伝染因子、およびトランスフォームされた癌細胞の同定、および特に発達障害の原因である人の染色体の異常の出生前診断において使用される。歴史的に、そのような分析は侵襲性の技術、たとえば疑わしい腫瘍の生検、および流産の比較的高い危険をもたらす羊水穿刺のような技術によってのみ可能だった。
血漿中での循環セルフリーDNA(circulating cell-free DNA)の発見は、前述の分析のために非侵襲性の技術の可能性を開いた。循環セルフリーDNAは、典型的に患者の血流中で少量見つかる。それは、ソースの混合物からのシークエンスを含み、しばしば2つのサイズから構成される。小さいサイズのフラクションは主に細胞のアポトーシスに由来し、より大きなサイズのフラクションは壊死細胞に由来すると考えられている。
循環セルフリーDNAは、少なくとも2つの領域において診断薬の可能性を持っている。第1は、腫瘍細胞に由来した循環セルフリーDNAは、分子的に腫瘍DNAをプロファイルするために「リキッドバイオプシー(liquid biopsy)」として使用されることができる。
第2に、循環セルフリーDNAは発達障害の原因である人の染色体の異常を発見するために使用されることができる。生まれていない胎児内の染色体異常の非侵襲性の検知についての増大する興味があり、これは胎児からのDNAが母親の血流の中にあるという観察によって可能になった。よく知られている例はダウン症候群に関連した染色体21トリソミーである。多くの方法がそのような染色体異常の同定および分析に利用可能である。出生後の欠失と増幅の検知のための最も広く用いられている方法は、マイクロアレイに基づいた比較ゲノムハイブリダイゼーション(aCGH)である。近年、いくつかの方法が、高度なシーケンシング方法を使用して、妊娠している女性からの循環セルフリーDNAの分析によって胎児の染色体異常を検知するために開発されている。
アレイCGHは、ゲノムの異なる部分をハイブリダイズするために設計された個々のDNAが多くのスポットで配置された固体表面を含むマイクロアレイを使用する。アレイCGHは2つのDNAプローブで標識を付ける。1番目は染色体異常の可能性のある対象のサンプルからのDNAである。2番目は正常なDNAであると考えられる参照DNAである。サンプルDNAは1つの蛍光染料で標識付けされる。参照DNAはスペクトルが別個の染料で標識付けされる。DNAプローブは組み合わされ、次にDNAマイクロアレイとコ−ハイブリダイズ(co-hybridise)される。未結合の物質を取り除くために洗った後に、マイクロアレイはDNAの個々のスポットの各々に存在する蛍光染料の量を測定するために走査される。サンプルの蛍光シグナルと参照の間の信号の比率の変化は、染色体の利得またはロスを示す。従来では、aCGHの中で使用される参照DNAは異なる個体からのDNAであることができた。しかしながら、これはコピー数変化を適切に考慮に入れないので、試験サンプルの測定をゆがませ、誤りの結果に結びつく場合がある。
血漿からの循環セルフリーDNAを分析する際に、問題がある。第1に、血流中の循環セルフリーDNAの量は、細胞ソース内のそれよりはるかに低い。第2に、出生前診断の場合、診断が可能になる時(つまり妊娠の初期段階で)には、母親自身のDNAの量と比較して胎児のDNAの量は少ない。同様に、正常な非腫瘍DNAと比較して、腫瘍DNAの量は少ない。第3に、胎児に由来した循環セルフリーDNAは小さいフラグメントであり、典型的にはアポトーシスの細胞からのフラグメントであり、一方、母親のDNAのほとんどは大きい。妊娠している女性では、例えば、胎児のDNAはおよそ140−180bpの寸法範囲内にあるが、フラグメントで見つかるほとんどの母親のDNAは約1kb以上である (Ellinger et al. (2008) Int. J. Cancer 122, 138-143)。循環セルフリーDNAのより小さなサイズは、より大きなサイズのDNAフラクションに比較して、いくつかの癌(例えば前立腺癌)ではより多くの腫瘍DNAを含むという示唆がある(Ellinger et al. (2008) Int. J. Cancer 122, 138-143)。
上に示されるように、現在の技術では、標準オリゴヌクレオチドアレイCGHはいくつかの理由で循環セルフリーDNAを分析するために使用することができない。アレイCGHはおよそ200ナノグラムのDNAの下限を有するが、5mlの血液サンプル中の循環セルフリーDNAの総量は、典型的には25ナノグラムである。増幅方法はアレイCGHと共に使用することができる。しかしながら、アレイCGHのための現在の増幅方法(例えばシグマによって販売されるRubicon GenomicsからのGenomePlex(登録商標)DNA増幅)は、循環セルフリーDNAの分析に使用される場合には、対象のDNA(例えば腫瘍DNAあるいは胎児のDNA)ではなく、長い背景の母親の/非腫瘍のDNAの増幅に好ましい。したがって、ターゲットとされていないフラクションの増幅が、ターゲットとされたフラクションの貢献を圧倒し、テストを敏感でなくしてしまう。
GenomePlex(登録商標)増幅は二段法である。第一段階は、PCRプライマー結合部位に融合(サンプルDNAに拘束すること)したランダムシークエンスを有するプライマーを使用したサンプルDNAのプライマーエクステンションである。プライマーエクステンションは鎖侵襲性のポリメラーゼを使用して実行される。増幅の第2段階は、ブライマー結合部位に相補的なプライマーを使用するPCR増幅から成る。
GenomePlex(登録商標)および同様の増幅方法は、アレイCGHと共に長年成功裡に使用されていた。しかしながら、これらの方法は腫瘍DNAの検出または非侵襲性の出生前診断において対象の循環セルフリーDNAを成功裡に増幅しない。したがって、最先端技術の増幅方法を使用するために、可能な限り物理的にサンプルから背景の母性DNAを取り除かなければならない。これは、たとえば、ゲル電気泳動法によって長いフラクションから短いフラクションを分離し、短いDNAを含んでいるバンドを削除して、次いでゲルからDNAを精製し、増幅することにより行うことができるかもしれない(Li et al. (2004))。EP 1524 321は、DNAの短いフラクションが、増幅に先立ってサンプル中のより長い核酸シークエンスから物理的に分けられる方法を開示する。WO 2009/097511は、母の血液サンプルから得られたDNAのDNase−処理により背景の母性DNAを取り除き、GenomePlex(登録商標)増幅を行う方法を開示する。
核酸の混合物内の特定の核酸の量を測定するための別法は、サンプルの組成物を分析するために塩基シークエンス決定を使用する。1つの方法では、サンプルは、循環セルフリーDNAのマッシブリィ パラレル ショットガン シークエンス法(massively parallel shotgun sequencing:MPSS)を使用して、増幅されシークエンスが決定される。この方法の1つの例では、循環セルフリーDNAは、高スループットショットガンシーケンシング技術を使用して、妊婦の血漿から直接シークエンスが決定された。1人の患者当たり平均で500万のシークエンス・タグが得られた。「リード」とも呼ばれる各シークエンス・タグ(短いDNAシークエンス断片)は、その元の染色体に対してマップされ、1つの染色体当たりのフラグメントが列挙される。各染色体についてマップされたシークエンス・タグの数を数えることによって、染色体の過剰または過小発現を検知することができる。ダウン症候群の場合には、各血漿サンプル中のDNAリードの総量内の染色体21リードのパーセンテージが計算される。胎児のDNA含量は通常血漿中のDNAの合計の4〜20%であるので、胎児によって寄与されたどんな増大も母体の貢献によって薄められる。この方法は、胎児のものと母体のDNAの間の分化を要求しない(Fan et al. (2008), Proc. Natl. Acad. Sci. USA 105, 16266-71 )。しかしながら、この方法は、多数のサンプルのバッチ処理を要求し、高価で、遅く、感度が低い。
循環セルフリーDNAで見い出されるように、小さな分子を含むサンプルを分析するために使用された時、フラグメントの端部へのアダプターをリゲートした後にPCRによってDNAを増幅することは標準的技法である。先行技術における、サンプルの少量のDNAを発見する次世代シーケンシング方法では、生成物の量を最大限にするために60〜90秒程度のアニーリング/エクステンション時間を使用することが標準的技法である。(Ehrich et al. (2011 ) Am. J. Obs. Gyn. 204, 205 e1 -e11 ; Bianchi et al. (2012) Obs. Gyn. 119, 890-901 ; Fan et al. (2010) Clin. Chem. 56, 1279-86)。
少量の核酸を含むサンプル内の特定のフラクションの量の測定に適用されることができるDNA分析の多くの他の方法がある。
WO 98/39474は、妊娠している女性の血流内の循環胎児DNAを利用する、非侵襲性の出生前診断の方法を開示する。血流から得られたサンプル中の母体のDNAおよび胎児のDNAを識別するために、この方法は、父親から継承される(したがって母体のDNAの中に存在しない)、胎児のDNAの内の特異的シークエンスの分析に依存する。
多くの先行技術方法がサンプル中のDNAの中の特定の前もって決定されたシークエンスを増幅することを含んでいる(たとえばWO 03/062441、米国公開特許2011/0312503、米国公開特許2013/0178371、米国公開特許2013/0143213およびWO 2007/147063を参照)。したがって、これらの方法は、増幅されるシークエンスについての知識を要求する。
他の例は、新しいシーケンシング方法による分析を含む(Chiu et al. (2011 ) Brit. Med. J., c7401 )か、あるいは量的PCRを使用してメチル化された胎児のDNAを母親のメチル化されていないDNAを識別し、母親と胎児の間の差別的にメチル化されたDNA領域をテストすること含んでいる(Papageorgiou et al. (2010) Nat. Med. 17, 510-13)。しかしながら、これらの既存の方法のすべてにおいて実用化の際には増幅することを必要とする。血流の中で循環するセルフリーの核酸の量は少なすぎるので増幅なしで分析することができないからである。WO2011/082386は、さらに胎児のDNAと比較して、母体のDNAのメチル化レベルの差に依存する方法を開示する。メチル化感受性の酵素が、提供メチル化の胎児のDNAを消化し、その後、PCRプライマー結合部位を含むリンカーがリゲートされてもよい。その後、リンカー−メディエイテッドPCR(Linker-mediated PCR)は優先的にメチル化されていないDNA(つまり消化され、リンカーを含んでいるもの)を増幅する。その後、再増幅はアレイハイブリダイゼーションに十分な生成物を供給するように要求される。
本発明は、循環セルフリー核酸のサンプル内の、対象のシークエンスを有する核酸の量が増加したかまたは減少したかを判断する改良方法を提供することを目的とする。本発明の態様によれば、循環セルフリー核酸のサンプル内の第1のソースから導かれた対象のシークエンスを有する核酸の量が増加したかまたは減少したかを判断する方法であって、該サンプルは第1のソースからの核酸と第2のソースからの核酸を含み、循環セルフリー核酸のサンプルはより短いサイズ範囲とより大きいサイズ範囲の2つのサイズ範囲内にある核酸を含み、最初のソースからの核酸はより短いサイズ範囲である方法であって、
a) 少なくともより短いサイズ範囲の核酸から増幅用テンプレートを形成すること;
b)増幅生産物を形成するためにテンプレートを増幅することにより、より短いサイズ範囲の核酸のためのサンプルを豊化すること;および
c)サンプルから導かれたコントロール増幅生成物の量と比較することにより、対象シークエンスを有する核酸の増幅生成物の相対量を決定する事を含み、ここでコントロール増幅生成物はコントロール・シークエンスを有し、該コントロール・シークエンスは対象シークエンスとは異なるものである方法が提供される。
この方法は、異なるサイズの核酸断片の異なる増幅を許可する。具体的には、増幅はサンプル中のより長いフラグメントの増幅を差別する。より短いサイズ範囲の核酸がより長いサイズ範囲にある核酸より効率的に(したがって優先的に)増幅される(したがって、その増幅は差別される)ように、ステップb)の増幅が実行される。その結果、より小さなフラグメントはより効率的に豊化され、分析の手段が何であっても、より長いフラグメントの背景に対するそれらの範囲内の特異的シークエンスを検知し測定することをより簡単にする。
対象シークエンスは、一般に分析されている(異常の部位)状態と異なると予想されるシークエンスである。コントロール・シークエンスは、一般に対象シークエンスの量が比較される試験サンプル中のシークエンスである。
理想的には、第1のソースからの核酸の本質的にすべてが増幅されるだろう。重要なことには、増幅は非特異的に行われる。したがって、増幅は、対象シークエンスを有する核酸のシークエンスに依存しない。増幅の効率は、DNA塩基シークエンスに依存せず、最初のテンプレートのサイズ、あるいは増幅の第一ラウンドの後に形成されたテンプレートのサイズに依存する。したがって、増幅ステップのためにサンプル内の特定核酸シークエンスの予備的知識を持っていることは本発明の方法では必要ではない。
第2のソースからのいくつかの背景核酸も、いずれかが小分子として存在する場合には増幅されていてもよいし、サンプル中の長い核酸断片は一般に効率的に増幅されないだろう。したがって、本発明方法の増幅ステップは、対象の小さな核酸についてサンプルを豊化するが、サンプル中に元々多量に存在する背景核酸を薄める。有利には、増幅に先立って長い核酸から物理的に短い核酸を分けるステップは必要ではない。多くの実施態様では、長い核酸も短い核酸も、ステップb)の増幅の最初のサンプルの中に存在する。いくつかの実施態様では、少なくともより長いサイズの範囲内にある核酸のうちのいくつかは、ステップbの増幅の最初のサンプルの中に存在する。本発明の方法で、より短いサイズの範囲内にある核酸だけが増幅されるように、ステップa)で短い核酸分子から形成されたテンプレートが選択的に増幅される。
増幅のためのテンプレートおよび/または増幅生産物は、たとえばおよそ750bp未満あるいはおよそ250bp未満のサイズを有する。特別の実施態様の中では、増幅のためのテンプレートおよび/または増幅生産物は、およそ140−180bpのサイズを有する。
いくつかの実施態様では、より短いサイズの範囲内にある最長の核酸は、より長いサイズの範囲内にある最短の核酸の1/2〜1/3である。好ましい実施態様では、ステップa)はサンプル中の核酸の端部へのアダプター分子をリゲートすることを含み、該アダプター分子は特異性のプライマー結合部位を含む。ステップb)はアダプター・プライマー結合部位に相補的なプライマーをアニーリングし、ポリメラーゼ連鎖反応を使用してテンプレートを増幅することを含んでもよい。
アダプターPCRの使用は、短いフラグメントの増幅に好適であり、増幅はサンプル中の核酸の特異的シークエンスに依存しない。
多くの実施態様では、増幅ステップの間にアニーリング/エクステンションに使用されたエクステンション時間は、当業者が通常使用しようとするものと比較して短くされる。組み合わされたアニーリング/エクステンション時間は30秒以下、20秒以下、10秒以下、または5秒以下である場合がある。第1のソースに由来した核酸のパーセンテージが低い場合(例えば4%以下)では、短くされたエクステンション時間は特に有用である。例えば、胎児のDNAの中の異常の検知の場合には、胎児のフラクションは変わる場合がある。これは、たとえば妊娠期間などの、母親の体質量指数に関連づけられることができる。胎児のフラクションは20から4%未満までと異なることある。エクステンション時間の低減によって、より短いDNA分子と比較して長いDNA分子の低減した増幅が見られる。結果として、胎児DNAを含まないものの検出限界は低下する。また(相対的に)少量の胎児DNAを含まないものを検知することができる。
コントロール・シークエンスは、好ましくは対象の核酸シークエンスと異なる染色体である。好ましい実施態様では、対象シークエンスは多形性(polymorphic)の疑いがなく、および/またはコントロール・シークエンスは多形性の疑いがない。
好ましい実施態様では、ステップc)は増幅生成物の量を、第2のソースからの核酸に対応する1セットの参照核酸から得られた増幅生成物と比較することを含んでおり、参照核酸は対象核酸シークエンスに対応するシークエンスを有する核酸を含む;ここで参照核酸は、対応するシークエンスを有する核酸の増加したか減少した量を含んでいない。好ましくは、コントロール・シークエンスと比較した対象のシークエンスからの増幅生成物の量の比率が得られ、これが参照核酸のセットから得られた同等の比率(equivalent ratio)と比較される。
参照核酸は、一般に同じサンプルから得られず、第1のソースからの核酸を含んでいないサンプルから得られる。「対応するシークエンス(corresponding sequence)」は、対象シークエンスと同一、または同一でないオルソロガスなシークエンス(orthologous sequence)をいう。好ましい実施態様では、参照核酸中の対応するシークエンスが対象シークエンスから成るプローブにストリンジェントな状態(stringent conditions)の下で結合することができるだろう。最も好ましい実施態様では、対象シークエンスは多形性ではない。したがって、対応するシークエンスが対象シークエンスと同一であると予想されるだろう。
実施態様では、ステップc)は以下を含んでいる:
i) 増幅生産物を、対象核酸シークエンスに本質的に相補的なテスト・プローブと少なくとも1つのコントロールプローブとハイブリダイズすること、ここでコントロールプローブはコントロールシークエンスに本質的に相補的である;
ii) 参照核酸を、テスト・プローブおよびコントロールプローブとコ−ハイブリダイズすること;
iii) サンプル/参考試験ハイブリダイゼーション比率およびサンプル/参照コントロール・ハイブリダイゼーション比率を得るために、テスト・プローブおよびコントロール・プローブとの、増幅生産物とのハイブリダイゼーションの量と参照核酸のハイブリダイゼーションの量とを比較すること;および
iv) 第1のソース内の対象の核酸シークエンスを有する核酸の量が増加したか減少したかを判断するために、テスト・ハイブリダイゼーション比率をコントロール・ハイブリダイゼーション比率と比較すること。
プローブは、たとえばチップあるいは複数のビーズのような固体担上に好ましくは提供される。固体担体が複数のビーズである場合、ビーズはそれぞれ単一のプローブ・シークエンスを含んでもよい。固体担体は多孔性表面を有することができる。
実施態様では、増幅生産物の量はマイクロアレイ・ハイブリダイゼーションを含んでいない方法を使用して比較される。例えば、本発明の方法は、ディジタルPCR方法を含む次世代シーケンシング方法あるいは定量PCR法を使用して、対象の核酸シークエンスの増幅生成物の相対量を決定することを含むことができる。
本発明の別の態様によれば、被検者の異常の可能性を予測する方法であって、上記の方法を使用して対象シークエンスを有する核酸の量が増加したか減少したかを決定することを含み、対象シークエンスを有する核酸の量の増加または減少がある場合には、被検者は異常の増加した危険を持つと予測される。
この方法のステップc)は、増幅生成物の量を、第2のソースからの核酸に対応する1セットの参照核酸から得られた増幅生成物と比較することを含んでいてもよく、そこでは参照核酸は対象の核酸シークエンスに対応するシークエンスを有する核酸を含む;ここで、参照核酸は、対応するシークエンスを有する増加したか減少した量の核酸を含んでいない。
参照核酸は被検者から得られてもよい。好ましくは、参照核酸は、第1のソースからの核酸を含むとは予想されない被検者からの細胞から得られる。例えば、参照核酸は、被検者の白血球あるいは被検者の頬側細胞から得られてもよい。実施態様では、参照核酸はヌクレオソームのDNAから抽出された。
方法の1つの実施においては、被検者は妊娠している女性である。また、異常はコピー数欠損のある胎児である。第1のソースからの核酸は胎児の核酸である。また、第2のソースからの核酸は母体の核酸である。
対象の核酸シークエンスは、例えば染色体13、染色体18、染色体21および/または、X染色体であることができる。コントロール・シークエンスは、染色体13、染色体18あるいは染色体21でない常染色体からであることができる。方法は、パトー症候群、エドワーズ症候群、ダウン症候群あるいはターナー症候群の出生前診断の非侵襲性の方法であることができる。
方法の別の実施において、被検者は腫瘍を持っている疑いがあり、第1のソースは腫瘍であり、第2のソースは非腫瘍細胞である。コントロール・シークエンスは好ましくは既知の良性のコピー数変化を含まないように選択される。実施態様では、ステップc)の中で得られたデータは、サンプル内の特異的プローブ・セットを、同じサンプルに由来したコントロールと比較するためにウィルコクソン検定を使用して分析される。
実施態様では、染色体用の単一のp値はストゥファー(Stouffer)法を使用して、p値を組み合わすことにより得られる。本発明の別の態様では、第1のソースからの核酸に対する循環セルフリー核酸サンプルの豊化方法であって、該サンプルは第1のソースからの核酸と第2のソースからの核酸を含み、セルフリー核酸のサンプルは、より短いサイズ範囲およびより長いサイズ範囲の2つのサイズ範囲の核酸を含み、第1のソースからの核酸はより短いサイズ範囲であり、
a) 少なくともより短いサイズ範囲の核酸から増幅用テンプレートを形成すること;および
b)テンプレートを増幅することにより、より短いサイズ範囲の核酸のサンプルを豊化し、対象シークエンスを含む核酸のヌクレオチドシークエンスと無関係の方法で増幅生産物を形成すること、を含む方法を提供する。
本発明の別の態様では、被検者から得られた血漿からの循環セルフリー核酸サンプル内の第1のソースからの対象のシークエンスを有する核酸の増加または減少を決定する方法であって、該サンプルが第1のソースからの核酸と第2のソースからの核酸を含み:
対象シークエンスを有する核酸の、第1のソースからのコントロール核酸との相対量を決定すること、ここでコントロール核酸はコントロール・シークエンスを有し、コントロール・シークエンスは対象シークエンスとは異なっている;および
対象のシークエンスがある核酸の相対量と、対象核酸シークエンスを有する第2のソースからの核酸に対応する参照核酸の相対量を比較すること;を含み、
参照核酸は被検者から得られる、方法を提供する。
本発明の好ましい実施態様は、核酸サンプル中の微量成分の検知および分析のためのaCGHの使用を可能にする。本発明の好ましい実施態様は、添付の図面を参照し、例示のみを目的として説明される。
図1は、本発明の方法の一般的な概観を説明するフローチャートである。 図2は、本発明の方法の実施態様のステップを説明するフローチャートである。 図3は、本発明の方法の実施態様のステップを説明するフローチャートである。 図4は、トリソミー21を検知するためにアダプター−PCRを使用する方法で得られた結果を説明する棒グラフである。 図5は、トリソミー21を検知するためにアダプター−PCRを使用する方法で得られた結果を説明する棒グラフである。 図6は、トリソミー21の検知のための異なるエクステンション時間を比較する棒グラフである。 図7は、ゲル電気泳動法によって分離され、より短いフラグメントの豊化を実証するDNA断片の写真である。 図8は、アダプターPCRの中でダイソミーDNAおよびトリソミー21 DNA生成物の産出を示すゲル・イメージである。 図9は図8の結果のグラフ表示である。 図10は、異なるエクステンション時間でのダイソミーとトリソミーのサンプル用染色体21にアラインした合計のリードのパーセンテージを示す棒グラフである。 図11は、異なるエクステンション時間でのダイソミーとトリソミーのサンプル用染色体21にアラインした合計のリードのパーセンテージを示す棒グラフである。 図12は、トリソミー/ダイソミー・サンプルの比率を異なるエクステンション時間で示す棒グラフである。 図13は、トリソミー/ダイソミー・サンプルの比率を異なるエクステンション時間で示す棒グラフである。 図14は、常染色体のすべてにわたる、トリソミー/ダイソミー・サンプルおよびダイソミー/ダイソミー・サンプルの平均比率を示す棒グラフである。
図1は、発明の好ましい実施態様に従うDNAを分析するために使用されてもよい典型的な方法の概観を説明するフローチャートである。第1のステップでは、循環セルフリーDNAが、被検者(例えば妊娠している女性あるいは腫瘍を持っている疑いのある人)の血流から得られたサンプルから抽出される。サンプルは対象のDNA(例えば胎児または腫瘍DNA)およびさらに背景DNA(母体または非腫瘍DNA)を含んでいる。DNAはサンプルから抽出される。DNAは下記に述べられた方法に従って増幅される。その後、サンプルDNAは分析される。これはいくつかの方法で行われてもよい。例えば、aCGH(背景DNAに対応する参照DNAを備えた)の中のマイクロアレイへのコ−ハイブリダイズすること、次世代シーケンシング方法を使用すること、あるいは定量PCRを使用すること。サンプル中の対象シークエンスに由来した増幅生成物の量は、コントロール・シークエンスに由来した増幅生成物と比較され、異常な量の対象DNA塩基シークエンス(期待される量より多いか少ない)があるかどうか確認する。
方法のステップは、aCGHを使用する方法に関して、および特に図2を参照して詳細に以下に説明される。
DNA抽出
方法の第1の実施態様では、DNAサンプルは妊娠している女性の血流から得られる。また、方法はトリソミー21の出生前診断の非侵襲性の方法であることができる。別の実施態様では、サンプルは、前立腺癌の疑いのある患者の血流から得ることができ、また方法は腫瘍DNAを検知する方法であることができる。方法は、トリソミー13またはトリソミー18のような、コピー数の変化によって引き起こされた他の先天異常の検知、あるいは他の腫瘍の検知に等しく適用可能である。どんな場合も、DNAは、任意の侵襲性のサンプリングの必要なしで、容易な方法で被検者の血液サンプルから簡単に得ることができる。DNAは、当該技術分野において有名な標準法(例えば、Qiagen QiAmp)を使用して、血液サンプルから抽出される。サンプルから抽出されたDNAは少量の対象のDNA(例えば胎児のDNAあるいは腫瘍DNA)を含んでいる。それは小断片の形態で、大量の背景DNA(例えば母体のDNAあるいは非腫瘍DNA)と存在し、それらの大多数はより長いフラグメントの形をしている。
参照DNAも準備される。ゲノム変化の規格化された個体群を表わすために、近年の公知の方法の中で使用される参照サンプルはしばしば多くの個体からのDNAの混合物からできている。出生前診断のための本発明の改良方法では、参照サンプルは、母体のDNA、好ましくは妊娠している女性自身のDNAへ非常にマッチするように好ましくは選ばれる。ソースの混合物を含む試験サンプル内の特定のフラクションの量を測定したい場合、問題がある。例えば、妊婦からの血液サンプルの血漿中のDNAの量は、約25ナノグラムである。また、典型的に、4−20%は胎児に由来する。したがって、妊娠している母親のそれにできるだけ、シークエンスが近似しているDNAを参照として使用することは重要である。この理由は、コピー数変化として知られている、ある染色体領域における量の変化がすべての人間において存在するからである。同様の関係は循環腫瘍DNA分析に当てはまる。参照におけるこれらの変化の存在は、もちろん、試験サンプルの測定をゆがめ、誤りの結果に結びつくだろう。比率(以下を参照)の最も正確な測定については、サンプルDNAに参照DNAのサイズおよび質を一致させることが重要である。
さらに、参照として使用されたDNAのサイズ範囲が、胎児または腫瘍に由来するものとマッチすることは重要である。胎児のフラクションは、アポトーシスによって発生するように見える;それは狭いサイズ範囲を有し、約150bpである。循環セル−フリーDNAを分析するための参照DNAのいくつかの潜在的なソースがある。
トリソミー21の出生前診断の方法用の参照DNAは、トリソミー21を含まないことが知られている母体のDNAであることができる。例えば、非トリソミー21胎児を妊娠している異なる女性からのDNAを使用することができる。胎児が1つは父親からの染色体と他方は母親からの染色体との染色体対を持っているので、すべてではないが、CNVの潜在的な貢献を、参照として母親のDNAを使用することにより幾分回避することができる。好ましくは、参照DNAは対象の妊娠している女性からのDNAを含み、参照DNAの中に胎児のDNAが存在する機会が最小限にされるように得られる。例えば、母体の白血球から参照DNAを得ることができる。妊娠している女性からの頬側スワブから参照DNAを得ることができる。このDNAは胎児のDNAを含んでいるべきでないが、母親のゲノムにおける任意のコピー数変化の存在が、データ分析ステップの間に緩和されることを可能にするだろう。したがって、これは胎児のDNAのより正確な分析に帰着する。
方法が腫瘍DNAを識別する方法である場合、参照DNAは腫瘍DNAを含んでいる疑いをかけられないDNAであることができる。再び、これは被検者からの頬側スワブから得られたDNAかもしれない。ヌクレオソームからDNAを抽出することにより、胎児のDNAのサイズ範囲とマッチさせることは可能であり、好ましくは母親の循環する白血球、あるいは上に述べられるような頬側細胞ソースから得られる。特に好ましい実施態様では、参照DNAの前述のソースのすべてについて、参照DNAはヌクレオソーマルDNA(例えば細胞からEZヌクレオソームDNA準備キット(ZymoResearch)を使用して、DNAを抽出することによる)から抽出されてもよい。それがアポトーシスのDNAに由来したDNAとサイズにおいて等しいので、これは小さいDNAフラクションとぴったりマッチするべきである。
同様に、循環腫瘍DNAのための参照として使用するのに適切なフラグメントは、患者の正常でトランスフォームされていない細胞から得ることができ、組織ソースから採取され、超音波処理あるいは他の手段によって適切なサイズ範囲のフラグメントに分解することができる。
好ましい実施態様では、参照DNAはサンプルDNAと同じ方法で抽出され、並行して処理される。
増幅
上に示されるように、DNAサンプルの中にある対象のDNAの量は少なすぎるので、従来のaCGHで直接使用することができない。
1つの実施態様では、DNAはアダプター−PCR法を使用して増幅される。ここでは特異性PCRプライマー結合部位を有するアダプターはサンプル中のDNAフラクションの端部にリゲートされる。この方法のステップは図2に述べられる。アダプターは、ヒトゲノムの任意の部分への低いシークエンス相同性のPCRプライマー結合部位を含ことを可能とする長さおよびGC含量のダブルストランドオリゴであり、A−テイルダブルストランドゲノムフラグメントへのリゲーションを助するT−オーバーハングを含む。プライマーの好ましい特徴は、それらのアニール温度がポリメラーゼの最適温度と同じであるということだろう。その結果、アニーリングと伸長は単一ステップとして実行されることができる。好ましいアダプター−PCR方法では、フラグメントは、当業者に知られている方法に従って、ギャップリペア、燐酸化、アデニル化および次にT−テイルアダプターのリゲーションにさらされる。
その後、PCR増幅が、短いエクステンション時間で好ましくは行なわれる。驚くべきことに、出願人は、開示した方法に最適な、当業者によって期待されるより短いエクステンション時間を決定した。好ましい方法では、アニーリングとエクステンションは同じ温度で単一ステップの中で実行されるだろう。出生前診断の方法では、好ましいアニーリングおよびエクステンション時間は循環胎児DNAについて150bp以上のフラグメントの増幅を達成するために必要とされるものを超過するべきでない。例えば、アニーリング/エクステンション・ステップの正確な長さは、たとえば使用されている増幅プロトコルおよびポリメラーゼの性質に依存する。当業者は、このゴールを達成するためにエクステンション条件を最適化する方法を容易に決定することができる。たとえば、典型的な条件の下で、Pfu DNAポリメラーゼのためのアニーリング/エクステンション時間は従来90秒台である。しかし本発明においては、アニーリング/エクステンション時間はおよそ30秒以下であることができる。30秒で、アニーリング/エクステンションの組み合わされた時間の正確なコーリングが達成された。出願人は時間を30秒未満にすると比率が改善され、したがって検知の感度が改善される。これは約700塩基を超えるフラグメントが増幅されないからである(下記の例3〜5を参照)。
短くされたエクステンション時間は、長いフラグメントの端部に達する前にエクステンション反応が終了することを保証する。長いフラグメントの端のプライマー結合部位はこのようにコピーされない。また、長いフラグメントのコピーは、第2の増幅サイクルのテンプレートとして利用可能ではない。したがって、この増幅方法は、サンプル中の長いフラグメントの増幅に失敗し、結果として短いフラグメントのサンプルの富化をもたらし、それは対象のDNAを含んでいる。このアプローチは、大きなDNA断片が容易に増幅され、そのため短いフラクション内に見い出される対象のDNAをさらに薄め、したがって出生前の診断または腫瘍同定に使用することができない、標準シークエンスオリゴヌクレオチドCGHと共に使用される先行技術の増幅方法の問題を解決する。次世代シーケンシング方法と共に使用された時、本発明の増幅方法は、さらに胎児のフラクション中の低量のDNAの検知を可能とする(以下を参照)。本発明の増幅方法は大きなフラグメントの増幅を差別するので、異なるサイズのフラクションが増幅に先立って物理的に差別的に分離される清浄化ステップの要求を回避する。
典型的に、aCGHおよびシークエンス分析の他の方法に関する当業者は、マイクロアレイ上で分析された時よい信号を与えるのに十分な大量の生成物を与える増幅の条件の最適化を求めている。本出願の増幅の好ましい方法は、より長いフラグメントの増幅の回避により増幅生成物の総量を低減するように設計されている。これは当業者の直観に反している。本発明の増幅方法の結果として、マイクロアレイ(以下を参照)上の信号はより低くなる。しかしながら、より小さな量の非目的信号に起因して、より低い背景ノイズレベル、および目的信号の改善された測定を提供する。
標識付け
増幅中または増幅後に、サンプルDNAは蛍光染料分子で標識付けされる、その後励起により特定波長の蛍光を放射する。参照DNAも標識付けされるが、異なる染料で標識付けされ異なる波長を放射する。染料は例えばCy3(緑)/Cy5(赤)であることができる。当業者は、他の適切な染料を使用することができることを認識するだろう。標準のアレイCGH標識付けは、標識が付けられたdNTPを組込むエキソヌクレアーゼを含まないクレノーDNAポリメラーゼ(Klenow DNA polymerase)のような、鎖侵襲性のポリメラーゼを使用する、プライマーエクステンションから成る。それはより短いフラグメントより、より長いフラグメントに好適であるが、これはうまくいく。例えばKreatechによって商業化されたULS標識付方法のような、サイズ・バイアスのより少ない化学リゲーションアプローチのような異なる標識付けアプローチを使用することができるかもしれない。
増幅中に標識を付けることができ、たとえばアダプターPCR工程中にPCRプライマーを蛍光性に標識付けすることができる。あるいは、蛍光性に標識が付けられたdNTPを加えることにより、PCRの間に蛍光性のdNTPを組込むことができる。アミノ−アリル−dNTPあるいはビオチン−dNTPを加えることのような非直接的な標識付けと、引き続く蛍光性に標識が付けられたNHSエステルまたは蛍光性に標識が付けられたストレプトアビジンによる標識付けを、それぞれ代替えとして使用することができる。
当業者に利用可能な他の多くの適切な標識付方法がある。
マイクロアレイへのコ−ハイブリダイゼーション(Co-hybridisation)
標識が付けられた増幅されたサンプルと参照DNAを、対象シークエンス(正常なソースと比較して増加したか減少した量で存在することが疑われるシークエンス)のプローブを含むマイクロアレイにコ−ハイブリダイゼーションし、プローブ(増加したか減少した量で存在するとは予想されない核酸シークエンスに相補的である)をコントロールする。
核酸アレイは複数のハイブリダイゼーションプローブであり、目標ヌクレオチドシークエンスがハイブリダイズされる固体表面上で不動化された。これは、単一反応コンパートメント内で同時に複数のプローブとサンプルが接触されることを可能にする。核酸アレイの調製および使用は当該技術において標準のものである。
アレイCGHは、ゲノムについて非常に低解像度レベルでコピー数変化を分析する高生産性方法である。コピー数における変化は位置特異的プローブ(locus-specific probes)を含むaCGHマイクロアレイで試験DNAおよび参照サンプルDNAをハイブリダイズすることによって測定される。例えば、胎児の染色体21のトリソミーを発見するのにふさわしいマイクロアレイは、染色体21上で見つかったシークエンスに特異的な少なくとも1つのプローブ(好ましくは少なくとも10、少なくとも100、または少なくとも1000プローブ)、および正常な2つのコピーで胎児の中にあると予想されるシークエンスに特異的なさらに少なくとも1つのコントロールプローブ(好ましくは少なくとも10、少なくとも100、または少なくとも1000プローブ)を含んでいるだろう。プローブは、一般に14塩基から350kb(例えば14塩基から500塩基、500塩基から150kb、あるいは150kbから350kb)までの範囲にあるだろう。
プローブ用の好ましい長さは60塩基であるかもしれない。
市販のマイクロアレイ(Oxford Gene TechnologyによるCytoSure ISCA v2(4x180k)アレイのような)は、先天異常(たとえば胎児のDNAの中のトリソミー21)の検知に、本発明の方法のために使用することができる。このマイクロアレイは、染色体21に特異的なプローブ、またさらに胎児のDNAの中の2つのコピー以外での存在を疑われないシークエンスに特異的なプローブを含む。CytoSure ISCA v2(4x180k)アレイは、単一のスライド上の4つの約154kプローブアレイを備えたスライドの形をしている。これらは、4つの別個の実験として個々にハイブリダイズすることができる。プローブの標的部位のうちのほとんどはゲノムに均一に存在し(平均のプローブの標的部位は25kb離れている)、より高い密度で細胞遺伝学者に興味のある約200のゲノム領域の小さい数の目標部位を有している。染色体21上のプローブ標的部位の数は1945である;染色体18においては3628である。また、染色体13においては4922である。残りの常染色体については、125026のプローブ標的部位がある。
腫瘍DNAの分析については、対象の核酸シークエンスへのプローブは癌に関連した染色体異常に由来するだろう。コントロール・プローブは癌に関連しない染色体領域に由来するだろう。CytoSure ISCA v2(4x180k)アレイのようなゲノムアレイの全体を使用することができるかもしれない。しかしながら、特定の癌に対するより大きな解像度で分析されることが重要なゲノムの部位がある場合、当業者はこれらの対象領域でプローブ密度を増加させたアレイを設計することができる。
マイクロアレイ上のプローブは、異なる最適化パラメーターで設計することができる。本発明の方法についてはこれは重要な点である。なぜなら、比較がプローブの緑および赤のチャンネルの間だけでなく、1つの染色体上のプローブについての蛍光比率を、異なる染色体上の同様に作動するプローブと比較する必要があるからである(以下を参照)。「同様に作動するプローブ」を構成するものがなんであるのかが議論のポイントである。性能に影響を及ぼす要因の中には、融解温度(おおむねGC含量に関連する)、染料バイアス(信号強度とともに変化してもよい)他がある。目標が理想的に等しい測定でゲノム領域を測定することであるので、設計はたとえば一塩基変異多型、コピー数変化、繰り返し領域、既知のブレークポイントおよび転座部位、高度の可変領域およびスプライス部位などの既知の領域を回避するために実行することができる。他方では、高度に保存された部位は、プローブ選択には有益である。
市販のアレイは本発明の方法の用途に適しているが、マイクロアレイ分析の騒音レベルを下げるさらなる改良をなすことができる。例えば、先天異常または癌を引き起こす形質転換をもたらす、ゲノムの領域をターゲットとするプローブの数は、それに適用可能なように増加されることができる。先天異常の場合には、chr 21(ダウン症候群)、chr 13(パトー症候群)および/またはchr 18(エドワーズ症候群)のプローブの数が多いことが好ましい。他の例では、4p(ウォルフヒルシュボン(Wolf-Hirschhorn)症候群)、22q11.2(ディ・ジョージ(diGeorge)症候群)および5p(ネコナキ(Cri du Chat)症候群)では染色体中の高められたプローブカバレッジが好ましいかもしれない。多くのそのような症候群が知られている。プローブの数を多くすることによって、分析の統計的検出力は増加される。また、従って、母体のDNAと共に血漿DNAの中にある胎児のDNAの中のこれらの症候群の検知に対する信頼性は高められる。本出願人は、過剰の他のDNAの内にあるターゲットDNAの与えられたレベルを検知するために、アレイ上に好ましく含まれるプローブの数を求めるために、プローブの数、データ中のノイズレベルおよびテストの感度の数値的関係を分析した。使用される具体的な実験・分析手法のために、信頼できる結果を得るために必要なプローブの数は調査することができる。1つの特別の例において、およそ2800のプローブはトリソミー21の信頼できるコーリングのために十分だった。一方、同じプローブ(5分の1に低減)のランダムサンプルは偽陰性の結果に帰着する不十分な精度を示した。ノンパラメトリック検定法(Wilcoxon)が使用された場合には、一般化されたモデルは生産することができない。しかし、当業者は、与えられた用途にふさわしいプローブセットを設計するために本出願の開示を利用することができる。
極端な塩基組成(例えば、高いGC含量、またはシークエンス繰り返し)、および貧弱なハイブリダイゼーション信号を備えたプローブの除去は、母体または非腫瘍DNAの背景の胎児または腫瘍DNAの検知に対する高められた信頼性に帰着する。これは、最適なプローブを選ぶことが分析(しばしばDLRSメトリックによってaCGHの中で測定された)中のノイズを減らし、比率測定の正確さおよび精度を改善するからである(以下を参照)。さらなる改良は、コントロールDNA領域のプローブのGC含量と目標DNA領域でのプローブのGC含量を合わせることにより達成され、これは比率の正確さおよび精度における増加に帰着する。
好ましい実施態様では、マイクロアレイのデザインは、既知の一塩基変異多型の領域およびコピー数変化を選択するのではなく、回避することにより、従来のaCGHデザインと比較して改善される。言いかえれば、アレイ上で好ましくは使用されるプローブは、すべての個体において同じ事が期待されるシークエンス(テスト・プローブおよびコントロール・プローブの両方)およびコピー数数(コントロール・プローブについて)むシークエンスを含んでいる。したがって、正確にマイクロアレイ上のプローブと一致するシークエンスだけが検知されるように、厳格なハイブリダイゼーション条件をこれによって使用することができる。一塩基変異多型およびコピー数変化のこの回避は、予想されるように単一性により近い信号の比率に結びつく(以下を参照)。従来のaCGHでは、対照的に、サンプル内の染色体のインバランスを発見するために、信号の比率が単一性とは異なるゲノムの領域の検知がゴールである。aCGHの中の一塩基変異多型を回避する原理は、サンプル中のゲノムの異常の検知の目的で直観に反している。先行技術方法(WO 98/39474に開示したような)は、サンプル中の母体と胎児のDNAを識別するために、胎児によって母方および父方から相続したシークエンスの違いを利用する。本発明のアプローチは、そのようなゲノム領域を回避し、胎児の中のより大規模な染色体のインバランスの検知のための新規で、改善された、非自明のプローブ選択を示す。
標識が付けられたサンプルDNAおよび参照DNAは、対象のDNAのためのプローブと複数のコントロールプローブを含むCGHマイクロアレイにコ−ハイブリダイゼーションされる。その後、マイクロアレイは洗われる。
洗浄の厳格性はテストプローブとコントロールプローブとの間で許容されるべきミスマッチのレベル、およびサンプル内で期待されるシークエンスにより決定される。上に示されるように、理想的には、プローブはどんなミスマッチもないように設計されている。それには厳格なハイブリダイゼーションおよび洗浄条件が使用されるべきである。洗浄の後、マイクロアレイは市販の装置、たとえばAgilentマイクロアレイ・スキャナで走査される。その後、個々のスポット上の蛍光の量は特別なソフトウェア(例えばAgilentの特徴抽出ソフトウェア)を使用して計算される。
アレイ上の特徴の走査および計量の後、分析されるゲノムDNA領域(例えばダウン症候群中の染色体21)についてサンプルDNAと参照DNAの間の比率が測定され、コントロール領域(例えば染色体14)で測定された比率と比較される。
比率のどんな差も、コピー数変化を含んでいるテストDNAの中の領域を示す。
従来のaCGHでは、コピー数の超過あるいは欠乏のいずれかを発見するために、各プローブ上の蛍光測定値は参照DNAのそれと比較される。
サンプルの蛍光の比率:
テスト・プローブ用の参照DNAの蛍光測定比率がこのようにして得られる。本発明の方法では、追加のステップが好ましくは行なわれ、試験サンプル:参照比率は、正常であると予想される染色体領域(つまりコントロール・プローブ)と同等の比率と比較され、テスト:コントロール比率が得られる。正常なサンプルについては、この比率の比率は単一性のものに近似しているべきである;単一性からのどんな有意の差も欠乏または超過を示す。
他の実施態様では、他の固体表面を使用することができる。固体表面は、ビーズの形であることができる。固体表面はその表面積を増加させるために多孔性であることができる。
データ分析
上に記述されるように、異なる色でテストサンプルと参照サンプルに標識を付けてこれらの信号の強度を測定することによって、参照サンプルと比較して、試験サンプルのDNAの中の倍数性レベルを定量することができる。この定量化ステップはコピー数変化が0のまわりで対称となるように、マイクロアレイの各プローブについて参照とサンプルのlog2比率を得ることを含む。理想的状態では、正常(コピー中立)のlog2比率は、log2(2/2)=0である。単一のコピー・ロスでは、log2(1/2)=−1である、また、単一のコピーゲインでは、log2(3/2)=0.58である。実際の用途では、測定誤差を計算した後でさえ、log2比率は理論数値と相当に異なることがある。
ほとんどプラットフォームに依存するaCGHデータの適切な規格化に続いて(例えばNeuvial et al. (2006) BMC Bioinformatics 7, 264 for spatial normalisationを参照)、行なわれた第一ステップは、システムの複雑な技術的なノイズを減らすためにデータを滑らかにすることである。異なる平滑法、たとえば四分位平滑法(quantile smoothing)およびウエーブレット(wavelets)を使用する異なる方法が提案されている (Eilers & de Menezes (2005) Bioinformatics 21, 1146-53; Hsu et al. (2005) Biostatistics 6, 211 -26)。分析の最も重要な部分は、異常なコピー数を示す領域(ゲノムのセグメント)を確実に検知する能力に依存する。主な目標は、ノイジーなシステムでの参照と比較して、サンプルの離散レベルによるコピー数変化を示す領域の境界を識別することである。ブレークポイント検出方法のセグメント化は、染色体のセクションのコピー数状態を検知し、異なるコピー数状態のセクション間の転移の位置を見つける。最も広く使用されたセグメント化方法は、循環的な二成分のセグメンテーションアルゴリズム(Olshen et al. (2004) Biostatistics 5, 557-72)である。一旦、セグメント化の結果が決定されたならば、各セグメントのコピー数の値または見込み、およびセグメントの重要なランが割り当てられる。
多くのアプローチがこの問題に取り組むために考案された。例えば、Hupe et al. (2004) (Bioinformatics 20, 3413-22) (2004)は、アダクティブウェイトスムージング方法(Adaptive Weights Smoothing procedure)に基づいた定数機能(piecewise constant function)をモデル化する、ガウスに基づいたアプローチ(Gaussian-based approach)を遂行する。コピー数が隠された状態にあるHidden Markovモデルを、ゲノムへの接近を考慮に入れて拡張される。Fridlyand et al. (2004) (Journal of Multivariate Analysis 90, 132-53)およびMarioni et al. (2006) (Bioinformatics 22, 1144-46)を参照。
上記のデータ解析方法は公知である。我々が本発明の中で取り組む問題は完全に異なるものである。最先端技術のデータ分析への変更は必要であり、それはここで詳述される。
ここに、染色体全体の多くの部分(恐らく全染色体でさえ)は、正常な参照サンプルに期待された値の若干上の緑−赤比率を持っていてもよい。しかし、この値は演繹的に知られておらず、本質的に過剰の対象DNAの濃度と関係のある連続的な変数である。log2比率の変化はlog2(1.05)=0.07のオーダーである場合がある。それはaCGHで典型的に考慮される値より(先の議論と比較して)相当に小さい。従って、多くのプローブからの結果を組み合わすことは、特定されていない実験に由来する信号中のノイズに結びつく実験の不確実性を克服するためにに有益である。値の分布に基づいた比率の真値の評価は、平均値、メジアン値、調整平均値およびM−エスティマー(M-estimators)のようなものと組み合わせてもよい。さらに、統計的手法は特別の個体群(たとえば母体の血液からの胎児のサンプルの染色体の異数性の検知の場合には、胎児がいるかどうか、前記異数性によって影響を受けないか)にサンプルが属するかどうかを判断するために使用されてもよい。これらの統計的手法は2つのクラスの方法を含むが、これらに制限されない。
統計分類および仮説検証
統計分類では、1つ以上のカテゴリーへのメンバーシップの問題は、多くの既知のサンプル(「トレーニング・セット」)の測定された特性を評価し、未知のサンプルのための分類関数を導き出すことにより解決される。回帰分析、クラスター分析、サポートベクターマシーンなどの当業者に公知の多くの方法があり、それら自身は指導されるかまたは指導されない学習方法を代表する。他方では、仮説検証は、仮定(帰無仮説)に対する利用可能なデータを評価し、データが、ランダム変化が観察と帰無仮説の期待値との間の差について説明することができないという仮説とデータが十分に異なるか否かの結論を求める。モデルに基づいた仮説検証方法(ステューデントt検定あるいはZ検定)および非母数試験方法(ウィルコクソン(Wilcoxon)ランク・テスト、マン・ホイットニーU検定他)がある。そのようなフリークエンティスト仮説テスト(frequentist hypothesis tests)に加えて、ベイズ推定が使用されてもよい。解決されるべき特定の問題によって、当業者は、上の1つの方法を選択してもよい。
好ましい実施では、ウィルコクソン検定を使用する非母数試験が、データ分析に使用された。ウィルコクソン検定はマイクロアレイ実験に使用された分析パッケージの一部であるが、その使用は以下の例に述べられているものとは異なることを指摘することは重要である。より具体的には、ウィルコクソン検定を含む方法を使用する統計テストは、先行技術の中では、1つのサンプルがコントロールとして働く複数のサンプル間の、特定のプローブまたはプローブの組の有意の差を決定する。一方本発明の方法では有意差試験が、同じサンプルに属するコントロールに対するサンプル内の特異的プローブ・セットのテストに適用される。
修正方法では、サンプルDNAから得られた増幅生産物はaCGHの代わりにシーケンシング方法を使用して分析される。例示の方法のアウトラインは図3のフローチャートで説明される。シーケンシング方法を使用する場合、DNAフラグメントにリゲートされたアダプター・シークエンスは、次世代シーケンシング技術によってフラグメントがシークエンスされることが可能または部分的に可能となる。優先的に短いDNA断片を増幅する短くされたエクステンション時間と共に、上に記述されるように、増幅ステップが行なわれる。この増幅ステップはさらにシーケンシングステップに必要とされてもよい追加のDNAシークエンスを加えるために使用されてもよい。増幅は、シークエンサーのフロー・セル上で直接起こることができる。あるいは、フロー・セル上でない場合には、DNAは次世代シークエンサー上へロードされてもよい。シーケンシングに続いて、対象フラグメントの数が、コントロールDNAからのDNA分子の数と比較される。
さらなる修正では、qPCRはサンプル中のDNAの量を決定するために使用される。背景DNA(母体/非腫瘍)と比較して、テストされている領域が対象DNA(胎児/腫瘍)DNAの中でよりメチル化され、その後増幅の中で短くされたエクステンション時間を使用する場合、長いフラクション(より母体の/非腫瘍DNA)と比較された、短いフラクション(より胎児の/腫瘍DNA)の増幅の差があるべきであり、そのため、より低い胎児の/腫瘍フラクションでトリソミー/腫瘍の検出限界を増加させる可能性に帰着する。
上記方法は、たとえば妊婦の循環セルフリーDNAの中の胎児のフラクションに見つかるような、ターゲットとされたシークエンスがより大きなターゲットとされていないフラグメントの背景内の小断片の分析を可能にする。本発明の方法では、DNAを切る1ステップは不必要である。本発明の方法はアダプターPCRのための短くされたアニーリング/エクステンション時間を有し、その結果、より長いフラグメントはより短いフラグメントより少ない量で増幅される。これらの増幅方法がうまくいくことは予期しえないものであった。小さな分子数のサンプルのシークエンス分析の当業者は、切断と長いPCRサイクルがサンプルから分子のより高い収率およびより一定のカバレッジを与えることを期待するだろう。
実施例
例1 − GenomePlex(登録商標)およびアダプターPCR増幅の比較
方法
増幅の差を説明するために、モデルシステムがセット・アップされた。切断されていない材料(長い)あるいは切断された材料(メーカーの指示に従ってCovarisソニケーターを使用して約180塩基対に切り、Qiagen PCR浄化カラムを使用して清浄にした)から成る、正常なゲノムDNA(Promega)が準備された。細胞ライン(Coriell)からのトリソミー21 DNA(T21)はCovarisソニケーターを使用して切断された。
実験は以下のようにセット・アップされた:
サンプルDNA:
・ 1ナノグラムの切断されたT21 + 25ナノグラムの切断されていないPromega DNA
・ 1ナノグラムの切断されたT21 + 25ナノグラムの切断されたPromega DNA
・ 25ナノグラムのPromega DNA(ネガティブコントロール)
参照DNAは切断されたPromega DNAを使用して準備された。サンプルおよび参照DNAはサプライヤーの指示に従ってGenomePlex増幅キット(シグマ)を使用して増幅された。ライブラリー調製溶液の2μlが加えられ、次いでライブラリー安定化溶液の1μlが続けて加えられた。DNAは攪拌され、2分間95℃のサーマルサイクラーに置かれた。氷の上でのインキュベーションの後、ライブラリー調製酵素の1μlが加えられた。また、DNAは次の条件の下でインキュベートされた:20分間15℃、20分間24℃、20分間37℃、5分間75℃。4℃でのインキュベーションの後、DNA調製物の15μlが次のステップに使用された。DNAの15μlに、10xの増幅のマスターミックスの7.5μl、水の47.5μlおよびWGAポリメラーゼの5μlが加えられた。DNAは3分間95℃でインキュベートされ、次に、15秒94℃および5分間65℃を14サイクルさせた。
その後、DNAはメーカーの指示を使用して、Qiagen PCR浄化カラムを使用して清掃化された。
その後、DNAの量はメーカーの指示に従って、NanoDropUV分光測光器を使用して定量された。その後、CytoSure標識付キット(Oxford Gene Technology)で1マイクログラムのサンプルDNAはCy3を使用して、1μgの参照DNAはCy5を使用して標識付けされた。DNA調製物は水を使用して18μlまで増やされ、Random Primerの10μlおよびReaction Bufferの10μlを加え、38μlの全容積とした。DNAは3分間95℃で変性された。氷の上の5分のインキュベーションに続き、以下の試薬が加えられ、2時間DNAがインキュベートされた:10μlのdCTP Labelling Mix、1μlのCy−dCTP 、および1μlのエクソフリーなクレノー・ポリメラーゼ(Klenow polymerase)。10分65℃のステップの後、サプライヤーの指示に従ってCytoSure浄化カラム(Oxford Gene Technology)を使用して、DNAは清浄化された。DNAは、4x180k CytoSure ISCAマイクロアレイ(Oxford Gene Technology)とのハイブリダイゼーションのために、スピードバック(speedvac)により乾燥され、40μlの水中で再懸濁され、コット1(Cotl)(Kreatech)の5μl、11μlのブロッキング薬(Agilent)および55μl 2xハイブリダイゼーション・バッファー(Agilent)を加えた。その後、DNAはサプライヤーの指示に従って、20rpmで回転するSurehybオーブン(Agilent)の中で24時間65℃でSurehybチャンバー(Agilent)のアレイ内でハイブリダイズされた。ハイブリダイゼーションに続き、マイクロアレイ・サンドイッチがCGHバッファー1(Agilent)下で取り外され、室温で5分間洗われた。その後、アレイはCGH2 Wash2(Agilent)で37℃で1分間洗われた。ついで、2μm解像度、16ビット、100%pmtでAgilentマイクロアレイ・スキャナで走査された。その後、TIFファイルはAgilentの特徴抽出ソフトウェアを使用して、特徴が抽出された。また、データはマイクロソフトエクセルを使用して分析された。アレイの各スポットのCy3とCy5の信号は検査され、350以下の強度の信号は削除された。Cy3/Cy5比率が計算された。比率の比率は以下のように計算した。比率の比率=(chr21の上のプローブの平均Cv3/Cv5比率)/(chr1〜20の上のプローブの平均Cy3/Cy5比率)
対照的に、実験は、GenomePlex(登録商標)ランダムプライマー−PCR方法の代替増幅方法を使用して実行された。
増幅の差を示すために、モデルシステムがセット・アップされた。切断されていない材料(長い)あるいは切断された材料(メーカーの指示に従ってCovarisソニケーターを使用して、約180 bpに切断され、Qiagen PCR浄化カラムを使用して清浄化した)から成る、正常なゲノムDNA(Promega)が準備された。細胞ライン(Coriell)からのトリソミー21 DNA(T21)はCovarisソニケーターを使用して切断された。実験は以下のようにセット・アップされた:
サンプルDNA:
・ 1ナノグラムの切断されたT21 + 25ナノグラムの切断されていないPromega DNA
・ 1ナノグラムの切断されたT21 + 25ナノグラムの切断されたPromega DNA
・ 25ナノグラムの切断されたPromega DNA(ネガティブコントロール)
参照DNAは切断されたPromega DNAを使用して準備された。サンプルおよび参照DNAはアダプター−PCR条件を使用して増幅された。DNAの量はメーカーの指示に従って、NanoDroUV分光測光器を使用して定量された。DNAサンプルおよび参照サンプルはOGT次世代シークエンシングセイブラリー調製を使用して増殖された。サンプル(混合サイズのヒトDNA、5%の切断されたT21 DNAとPromega、Coriell細胞容器(Coriell Cell Repositories))および参照DNA(切断されたヒトDNA、Promega)は、市販のDNAポリヌクレオチドキナーゼおよびDNAポリメラーゼを使用して、従来のプロトコルに従って端部変性されA−テイルされた。DNAアダプター(例えばAgilent SureSelect XT試薬キット)が、各アダプターハイブリッドおよび1200U T4のRapid DNA Ligase(Enzymatics)の1.9 pmolを直ちに加えて15分間20℃でインキュベートすることにより、端部を変性され、A−テイルされたDNAにリゲートされた。生成物はPCR Qiaquickカラムを使用して清浄化されて、30μlのQiagen溶離バッファー中に溶出した。その後、リゲートされたDNAの半分は1x Accuzyme反応緩衝液(Bioline)、100 pmolのアダプター・シークエンスをターゲットにしたそれぞれのPCRプライマー(例えばAgilent SureSelect XT試薬キット)、50nmolのdNTPsl(Bioline)および5U のAccuzyme DNAポリメラーゼ(Bioline)の中で使用された。サイクリング条件は以下の通りだった:30分の98℃、引き続く98℃30秒の変性、引き続く30秒65℃のアニーリングおよび1分72℃のエクステンションの15サイクル。PCR反応はQiaQuickカラム(Qiagen)を使用して、清掃され、30μlの中でQiagen溶離バッファーに溶出した。その後、DNAの量はメーカーの指示にしたがってNanoDrop UV分光測光器を使用して定量された。CytoSure標識付キット(Oxford Gene Technology)で1マイクログラムのサンプルDNAはCy3を使用して、1μgの参照DNAはCy5を使用して標識付けされた。DNA調製物は水を使用して18μlまで増やされ、Random Primerの10μlおよびReaction Bufferの10μlを加え、38μlの全容積とした。DNAは3分間95℃で変性された。氷の上の5分のインキュベーションに続き、以下の試薬が加えられ、2時間DNAがインキュベートされた:10μlのdCTP Labelling Mix、1μlのCy−dCTP、および1μlのエクソフリーなクレノー・ポリメラーゼ(Klenow polymerase)。10分65℃のステップの後、サプライヤーの指示に従ってCytoSure浄化カラム(Oxford Gene Technology)を使用して、DNAは清浄化された。DNAは、4x180k CytoSure ISCAマイクロアレイ(Oxford Gene Technology)とのハイブリダイゼーションのために、スピードバック(speedvac)により乾燥され、40μlの水中で再懸濁され、コット1(Cotl)(Kreatech)の5μl、11μlブロッキング薬(Agilent)および55μl 2xハイブリダイゼーション・バッファー(Agilent)を加えた。その後、DNAはサプライヤーの指示に従って、20rpmで回転するSurehybオーブン(Agilent)の中で24時間65℃でSurehybチャンバー(Agilent)のアレイ内でハイブリダイズされた。ハイブリダイゼーションに続き、マイクロアレイ・サンドイッチがCGHバッファー1(Agilent)下で取り外され、室温で5分間洗われた。その後、アレイはCGH2 バッファー1(Agilent)で37℃で1分間洗われた。ついで、2μm解像度、16ビット、100%pmtでAgilentマイクロアレイ・スキャナで走査された。その後、TIFファイルはAgilentの特徴抽出ソフトウェアを使用して、特徴が抽出された。また、データはマイクロソフトエクセルを使用して分析された。アレイの各スポットのCy3とCy5の信号は検査され、350以下の強度の信号は削除された。Cy3/Cy5比率が計算された。比率の比率は以下のように計算した。比率の比率=(chr21の上のプローブの平均Cv3/Cv5比率)/(chr 1〜20の上のプローブの平均Cy3/Cy5比率)
GenomePlex増幅およびアダプター−PCRを比較した結果は、以下の表1に示される:
結果は、GenomePlex増幅を使用する時、4%の切断されたT21と25ナノグラムへ切断されたゲノムDNAで得られた比率の比率は、1.020の予期された理論的な比率に接近している。対照的に、25ナノグラムの長いゲノムDNA内に、4%の切断されたT21がスパイクされた時に得られた比率の比率は1.000に近かった。この結果は、GenomePlex増幅方法が長いゲノムDNAの増幅に好ましいことを示す。これは、対象のDNAが循環セルフリーDNAの短いフラクション内に見つかる場合には、aCGHの使用のための既知の増幅方法が適切ではないことを実証する。
長いゲノムDNAに4%の切断されたT21DNAをスパイクしたアダプターPCR増幅方法を使用する時、理論的な比率よりも有意に大きな比率を与えた。これは、アダプターPCR増幅方法が短いDNAの増幅に好適であり、混合されたサイズのDNAを含むサンプルの短いDNAフラグメントを豊化するために使用することができることを示す。
出願人は、たとえばGenomePlexのような標準PCR方法の使用は十分にうまくいかないことを示した。これは、長いDNA(例えば母体のDNA)からの背景の信号が、短いDNA(例えば胎児のDNA)からの信号を圧倒するからである。
例2 −トリソミー21の検知用のアダプター−PCR
Wessex National Genetics Reference Laboratory (NGRL)により提供された、ダイソミーあるいはT21胎児のいずれかを有していると知られている妊婦からの冷凍のセルフリープラズマ・サンプルのメーカーの指示に従って、QIAmp Circulating Nucleic Acid (Qiagen)キットを使用して、1つのサンプル当たりの5mlのプラズマからのDNA抽出を、1つのトリソミー・サンプルおよび3つのダイソミー・サンプルに対して実行した。
DNAは10Ομlヌクレアーゼフリーの水中で溶出され、SpeedVacを使用して乾燥し、15μlヌクレアーゼフリーの水で再構成された。各サンプルの全体は従来のプロトコルに従って端部変性されA−テイルされた。DNAアダプター(例えばAgilent SureSelectXT試薬キット)は、市販のDNAリガーゼ(たとえば、 T4 Rapid DNA Ligase, Enzymatics)を使用し、メーカーの指示に従って、端部変性されA−テイルされたDNAにリゲートされた。QiaQuickカラム(Qiagen)を使用して、反応物は清浄化されて、30μlのQiagen溶離バッファー中に溶出した。その後、リゲートされたDNAの半分は、1x Accuzyme反応緩衝液(Bioline)、アダプターシークエンスをターゲットとするPCRプライマー(例えばAgilent SureSelectXT試薬キット)の各々の50−200pmol、50nmolのdNTPおよび5U Accuzyme DNAポリメラーゼ(バイオライン)がPCR反応の中で使用された。サイクリング条件は以下のとおりだった:3分間の98℃、続いて30秒間98℃の変性、30秒間65℃のアニーリング、および1分間72℃のエロンゲーションを15サイクル。PCR反応物はQiaQuickカラム(Qiagen)を使用して清浄化され、30μlのQiagen溶離バッファー中に溶出された。
その後、DNAの量はメーカーの指示に従ってNanoDrop UV分光測光器を使用して定量された。その後、サンプルDNAの1マイクログラムはCy3を使用して標識付けされた。CytoSure標識付キット(Oxford Gene Technology)で1マイクログラムのサンプルDNAはCy3を使用して、1μgの参照DNAはCy5を使用して標識付けされた。DNA調製物は水を使用して18μlまで増やされ、Random Primerの10μlおよびReaction Bufferの10μlを加えた。DNAは3分間95℃で変性された。氷の上の5分のインキュベーションに続き、以下の試薬が加えられ、2時間DNAがインキュベートされた:10μlのdCTP Labelling Mix、1μlのCy−dCTP、および1μlのエクソフリーなクレノー。10分65℃のステップの後、サプライヤーの指示に従ってCytoSure浄化カラム(Oxford Gene Technology)を使用して、DNAは清浄化された。DNAは、4x180k CytoSure ISCAマイクロアレイ(Oxford Gene Technology)とのハイブリダイゼーションのために、スピードバック(speedvac)により乾燥され、40μlの水中で再懸濁され、コット1(Cotl)(Kreatech)の5μl、11μlブロッキング薬(Agilent)および55μl 2xハイブリダイゼーション・バッファー(Agilent)を加えた。その後、DNAはサプライヤーの指示に従って、20rpmで回転するSurehybオーブン(Agilent)の中で22時間65℃でSurehybチャンバー(Agilent)のアレイ内でハイブリダイズされた。ハイブリダイゼーションに続き、マイクロアレイ・サンドイッチがCGHバッファー1(Agilent)下で取り外され、室温で5分間洗われた。その後、アレイはCGH2 バッファー2(Agilent)で37℃で1分間洗われた。ついで、2μm解像度、16ビット、100%pmtでAgilentマイクロアレイ・スキャナで走査された。その後、TIFファイルはAgilentの特徴抽出ソフトウェアを使用して、特徴が抽出された。また、データはマイクロソフトエクセルを使用して分析された。
データを分析するために、アレイの各スポットのCy3およびCy5信号が検査され、人の染色体にハイブリダイズするようには設計されなかったすべてのプローブ、および350より下の緑信号のプローブは削除された。緑/赤(Cy3/Cy5)比率はすべてのプローブについて計算された。また、平均の緑/赤比率は染色体21に対応するプローブ、およびXおよびY(参照プローブ)以外の他のすべての染色体に対応するプローブについて計算された。その後、比率(Chr21/参照)の比率は各サンプルセットについて計算された。結果は図4に示され、ダイソミー21コントロールと比較して、染色体21 DNAの増加分を検知するために本発明の方法を使用することができることを実証する。
現在までテストしたすべての臨床のサンプルについては、図5に示される結果が得られた。箱のプロットは、すべての読みの平均値を表わす中央のラインと共に、各グループの比率の4分位数の1番目と3番目の数の間の範囲を表わす。エラーバーは、各コホート(cohort)で得られる最大値と最小値を示す。これは、アダプターリゲーションに続く増幅の本発明方法が、前述のサンプルの中で母体のプラズマを使用することにより、胎児の中のトリソミーとダイソミーの21の正確な非侵襲的検知を可能にすることを実証する。
例3 − 短くされたエクステンション時間でのアダプター−PCR
サンプル(混合サイズのヒトDNA、5%の切断されたT21 DNAとPromega、コリエル セル レポジトリーズ(Coriell Cell Repositories))および参照DNA(切断されたヒトDNA、Promega)は、例2に記述されたように端部変性(end-repair)、A−テーリング(A-tailing)およびアダプターリゲーションにより処理された。清浄化されたリゲーション生成物は、異なるサンクリングパラメーターで、以前に記述されたようなPCR反応の中で使用された。サイクリング条件は以下のとおりだった:3分間の98℃、続いて30秒間98℃の変性、30秒間65℃のアニーリング、および30秒間72℃のエロンゲーションを15サイクル。1分間のエクステンション時間は、1秒エクステンション時間に対してテストされた。PCR反応は、1.8X AMPure XP磁気ビーズ(Beckman Coulter)を使用して清浄化され、30のμlヌクレアーゼフリー水中に溶出された。
DNAの量はメーカーの指示に従って、NanoDropUV分光測光器を使用して定量された。その後、CytoSure標識付キット(Oxford Gene Technology)で1マイクログラムのサンプルDNAはCy3を使用して、1μgの参照DNAはCy5を使用して標識付けされた。DNA調製物は水を使用して18μlまで増やされ、Random Primerの10μlおよびReaction Bufferの10μlを加えた。DNAは3分間95℃で変性された。氷の上の5分のインキュベーションに続き、以下の試薬が加えられ、2時間DNAがインキュベートされた:10μlのdCTP Labelling Mix、1μlのCy−dCTP、および1μlのエクソフリーなクレノー。10分65℃のステップの後、サプライヤーの指示に従ってCytoSure浄化カラム(Oxford Gene Technology)を使用して、DNAは清浄化された。DNAは、4x180k CytoSure ISCAマイクロアレイ(Oxford Gene Technology)とのハイブリダイゼーションのために、スピードバック(speedvac)により乾燥され、40μlの水中で再懸濁され、コット1(Cotl)(Kreatech)の5μl、11μlのブロッキング薬(Agilent)および55μl 2xハイブリダイゼーション・バッファー(Agilent)を加えた。
その後、DNAはサプライヤーの指示に従って、20rpmで回転するSurehybオーブン(Agilent)の中で22時間65℃でSurehybチャンバー(Agilent)のアレイ内でハイブリダイズされた。ハイブリダイゼーションに続き、マイクロアレイ・サンドイッチがCGHバッファー1(Agilent)下で取り外され、室温で5分間洗われた。その後、アレイはCGH2 バッファー2(Agilent)で37℃で1分間洗われた。ついで、2μm解像度、16ビット、100%pmtでAgilentマイクロアレイ・スキャナで走査された。その後、TIFファイルはAgilentの特徴抽出ソフトウェアを使用して、特徴が抽出された。また、データはマイクロソフトエクセルを使用して分析された。データを分析するために、アレイの各スポットのCy3およびCy5信号が検査され、人の染色体にハイブリダイズするようには設計されなかったすべてのプローブ、および350より下の緑信号のプローブは削除された。緑/赤(Cy3/Cy5)比率はすべてのプローブについて計算された。また、平均の緑/赤比率は染色体21に対応するプローブ、およびXおよびY(参照プローブ)以外の他のすべての染色体に対応するプローブについて計算された。その後、比率(Chr21/参照)の比率は各サンプル・セットについて計算され、表2に示された結果が得られた。
これらの結果も図6の棒グラフ中で説明される。結果は、60秒から1秒のエクステンション時間のトリソミー・サンプルでの比率の著しい増加、およびダイソミーサンプルでの比率の著しい減少を示す。これは短いエクステンション時間が使用される場合、トリソミーとダイソミーの比率間の違いが増加することを導く。
例4 −エクステンション時間の長短の比較 I
Hyperladder(登録商標)(Bioline)のサンプルが使用され、上記のように端部変性、Aテーリング処理され、リゲーションされた。異なるアニーリング/エクステンション時間で、上に記述されるようにPCRが実行された。また、サンプルはD1Kテープステーション(Agilentテクノロジーズ)上で実行された。
結果は図7に示される。それはTapeStation分析ソフトウェア(Agilentテクノロジーズ)を使用して得られたゲル・イメージである。増幅標準状態(90秒のアニーリング/エクステンション:30秒の65℃、60秒の72℃)の下で、サイズ範囲<350ntの生成物のフラクションはおよそ36%である。対照的に、より短いアニーリング/エクステンションレジーム(10秒のアニーリング/エクステンション:10秒の65℃)の下では、約61%はこのフラクションにある。これは、より短いアニーリング/エクステンションレジームでのより長いフラグメントの増幅に対するバイアスによる。具体的には、より長いアニーリング/エクステンションでは36%は>700ntであり、より短いアニーリング/エクステンションでは存在しない。
例5 −エクステンション時間の長短の比較 II
混合サイズ・ダイソミーヒトDNAサンプル(Promega)および、4%の切断されたT21DNA(Coriell Cell Repositories)を含む混合サイズ・ダイソミーヒトDNAサンプル(Promega)、および参照DNA(切断されたヒトDNA、Promega)が、上に記述されるように端部変性、Aテーリングおよびアダプターリゲーション処理された。清浄化されたリゲーション生成物は異なるアニーリング/エクステンション条件で、先に述べたようなPCR反応で使用された。サイクリング条件は以下のとおりだった:3分間98℃、次いで30秒98℃の変性、65℃アニーリング/エクステンションの15サイクル。テストされたアニーリング/エクステンション時間は、30秒、20秒、10秒および5秒だった。PCR反応生成物は、AMPure XP磁気ビーズ(Beckman Coulter)で清浄化され、標識が付けられ、アレイに適用され、先に述べたように信号データは分析された。結果として得られた比率は、以下の表3に示される。
この実験は、アニーリング/エクステンション時間が減少するにつれて、Chr21/Ref比率が著しく改善されることを示す。
例6 − プローブ最適化
マイクロアレイのためのプローブ選択の間の1つの目的は、比較の結果に敏感であると知られている要因が、比較プローブ選択の間に一定に保たれるということでありえる。実際上、我々は最適化プローブセットを備えたマイクロアレイにハイブリダイゼーション実験を実行した。この最適化プローブセットは最終設計のために必要であるものよりはるかに多くのプローブを含んでいた。その後、この最適化プローブセットからの選択規則は、対象の染色体(例えば1セットのコントロール染色体と染色体13、18および21)にわたる変数(例えばGC含量)に関して同様の分布を維持し、同時にいくつかの最適化実験に基づいたプローブの平均の赤/緑信号のための制約を反映するために選ばれた。
例えば、GC含量は厳重にセット・ウィンドウ(例えば20%<GC含量<60%)から選ぶことができる。あるいは、複数のそのようなウィンドウをセットすることができ、各ウィンドウ内のプローブの割合をすべての染色体組にわたり維持することができる。同様のアプローチは信号の強度にも行い、最適化実験に信号のしきい値をセットするか、またはそれぞれのプローブ(例えば0−500、500−1000、1000−1500など)の平均またはメジアン強度の1つ以上のウィンドウを設け、その間異なる染色体あるいは対象の領域からのプローブの割合が維持されることができる。当業者はこれらの例から、パフォーマンス基準(例えばGC含量、または信号レベル)に基づいてプローブセットを一グループにすることができ、プローブのバランスのとれた設計を許容し、対象のプローブセットおよび参照セットが個々のそのようなグループの内でプローブの同様の数を有し、さらなるダウンストリームにおいてバランスのとれたセットの分析を許容する。
先行技術とは対照的に、本発明の方法では、我々は、プローブごとにサンプルと参照の間の差を示す領域を捜す。
例7 − データ分析
単純なデータ分析はここで示された実験的方法で行なわれてもよい。当該技術分野において周知のように、データ分析は抽出されたマイクロアレイ・データの特徴把握からスタートする。与えられたサンプルについては、下記手続きが母体の血液サンプルからの胎児のトリソミーおよびダイソミー・サンプルの間の差を見つけるために使用された:
(a)プローブごとに、すべてのプローブについて緑と赤の比率を計算する;
(b)実験条件(例えば500以上の任意の蛍光性のユニット(a.f.u)の平均の赤および緑信号)に適切な信号のしきい値を課することによって、実験的な基準によってデータセットをフィルターする;
(c)プローブを、それらのGC含量(これがアレイに対するプローブ選択の間にまだ行われていない場合)によってフィルターする;
(d)対象の染色体(例えば染色体21)について、フィルターされたG/R比率の平均およびメジアン値を計算し、同様に得られた値により規格化し、あるいは、固定された参照染色体(たとえば染色体14)、あるいは1セットの適切に選ばれた染色体(例えば異数性(つまり、好ましくは染色体13および18を除外して)の重要な危険のないすべての常染色体)により規格化する。
これらの「比率の比率」は、調査されるトレーニング・セットについて、明瞭な分離を示す。それは特別の方法の分類されたしきい値を知るために使用することができる。例として、ダイソミー(「正常な」場合)とトリソミーの間の分類しきい値が約1.03に選ばれるかもしれない場合、図5に示されるデータはこのしきい値を説明するために役立つ。
例8 − ウィルコクソン検定を使用するデータ分析
同じデータセットのデータ分析の代替例は、ウィルコクソン検定に基づく。ここでは、下記手続きが使用された:
(a)プローブごとのベースで、すべてのプローブについて緑と赤の比率を計算する;
(b)実験条件により、例えば500以上の任意の蛍光性ユニット(a.f.u)の平均の赤および緑信号)に適切な信号のしきい値を課することによってデータセットをフィルターする;
(c)アレイに対するプローブ選択の間にまだ行われていない場合は、それらのGC含量によってフィルターする;
(d) 参考について、
(e)対象の染色体(例えば染色体21)についてのフィルターされた緑と赤の比率を規格化するためにステップ(d)で得られた値を使用する;
(f)ダイソミーを持った胎児のための期待値とデータが異ならない帰無仮説に対する(e)で得られた値の群についてウィルコクソン検定を行う;理論上、この値は1である;しかしながら、ある状況で、例えばダイソミー妊娠からのサンプルのトレーニングデータに対して実行された前の実験が示す場合に、期待値が1とは異なる状況が発生してもよい;その場合、前の経験に基づいた帰無仮説の値への必要な調節を行うことは許容可能である;
(g)タイプ1エラー(偽陽性)に対する適切なリスクを決定するためにp値のしきい値を選ぶ。
そのようなしきい値の正確な値は実験の詳細に依存し、テストの要求性能によって選ばれるだろう。当業者は、そのようなしきい値を決定するためにROC曲線を評価することができるだろう。例において、しきい値は10−20を選ぶことができる。それは典型的な有意差試験の中で使用されるp値カットオフより本質的に小さい。例は以下の表4に示され、患者のサンプル用の比率およびp値を示す。
好ましい実施態様において、ウィルコクソン検定に対してステューデントt検定が先に記述した分析法において使用される時、偽陽性の数および偽陰性の数はウィルコクソン検定方法と比較してより多くなる。
例9 − トリソミー21/ダイソミー 21比率に対するアニーリング/エクステンション時間の影響の比較
この実験で、切断されたダイソミーDNAあるいはT21 DNA(胎児のDNAのモデル)のいずれかを10%のスパイク・インとして、(2.5ナノグラム)の25ナノグラムのサイズ・ミックスDNA(母体のDNAのモデル)に加えた。サイズ・ミックスは1:1:1重量比の高分子量:約1000 bp:180 bpだった。その後、サンプルはPCR1とPCR2の30秒または5秒を用いて、全ゲノム次世代シークエンスによりNGS分析のために準備された。より短いエクステンション時間がより短いシークエンスの増幅に好ましく(あるいは実際のサンプルでは、胎児由来のフラグメントの、増加したT21/D21比率に帰着することは本発明の出願人によって提案されている。
方法
端部変性/A−テーリング反応ミックス中の量は、以下の表8に与えられる。
1つのT21および1つのダイソミー・サンプルが準備され、30分間20℃、次いで30分間72℃でインキュベートされた。サンプルは氷の上に置かれた。リゲーション混合物は、2μlのリガーゼ・バッファーを、4μlのAgilentアダプター(1:5)および4μlのT4 Rapid Ligaseに加えることによりセット・アップされた。リゲーション混合物の5マイクロリットルが各サンプルに加えられた。
生成物は30分インキュベートされ、1.8 Xのビーズで清浄化されて、次いで30μl水中に溶出された。サンプルは半分に分割され、Agilentプライマーを使用してSurecyclerで標準のランピング速度で30秒または5秒で10サイクル増幅された。(過剰に増幅された生成物が生成されないように、サイクル数が最適化された。)サイクルパラメーターは、3分の0の98℃と次いで98℃で30秒および65℃での30秒または5秒の10サイクルだった。
Agilent SureSelect PCR1試薬の5X マスターミックスは、以下により35μlにされ、15μlのライブラリーに加えられた(下記の表9を参照)。
1.8 X アンピュアを使用してサンプルは清浄化されて、Qubit broad range DS DNAキットを使用して定量された。PCR2については、10ナノグラムのPCR1生成物の二重反復試験が、上記の方法を使用して、6サイクル(8および10サイクルが過剰な増幅に帰着した)で増幅された。またlllumina MiSeqデスクトップの次世代シークエンサーの多重処理のためのインデックスを付けられたバーコードが組込まれた。MiSeqシーケンシング用のプールは、以下の表10に示される。
MiSeqプール中の目標の最終濃度は4nMだった。したがって、サンプルはそれぞれ4/8=0.5nMだった。各フラグメントの両方の鎖の100の塩基がシークエンスされた。シークエンスされたフラグメントあるいは「リード」は、ヒトゲノム中の染色体に整列された。各染色体のリードの合計が計算された。染色体21に整列されたリードの数は合計の整列されたリードの割合として計算された。
結果
表11はPCR1(QubitブロードレンジDS DNAキット)の後のサンプルの濃度を示す。
Qubit上で測定された時のDNA濃度は、30秒のエクステンションからの産出が5秒のエクステンションのほとんど3倍であることを示す。表12、図8および9は、PCR2(QubitブロードレンジDS DNAキット)の後のサンプルの結果を示す。
結果は図8および9にも示される。ここでAOはラダー(Ladder)であり、A1がT21、5秒、インデックス1,B1がT21、5秒、インデックス2,C1がダイソミー、5秒、インデックス3、D1はダイソミー、5秒、インデックス4、E1はT21、30秒、インデックス5、F1はT21、30秒、インデックス6、G1はダイソミー、30秒、インデックス7およびH1はダイソミー、30秒、インデックス8である。
すべてのサンプルの最初の量が10ナノグラムであるにもかかわらず、PCR2の30秒のエクステンションの後には、5秒のエクステンションのほぼ3倍の大きな産出を示す。5秒および30秒のエクステンションのピーク目標長さに本質的な違いはなかった。しかしながら、より短いエクステンション時間では高い分子量生成物に大幅な減少があった。
シークエンスデータは以下の表13および14にまとめられる。表13は各染色体に整列した列挙されたフラグメントを示す。表14は各染色体に整列した合計のフラグメントのパーセンテージを示す。
各サンプルの整列したリード(XとY染色体以外の)の総数は、イタリック体で示される。範囲は1つのサンプル当たり2〜400万のフラグメントあるいはリードである。
シークエンスレーン中の整列したリードの総数は21,716,666である。
染色体21は、リードの合計のちょうど1%以上(染色体のサイズから期待され、他のものによって報告されたように)を占める。T21スパイクがD21よりわずかに高い値を持っていることが見られる。これは、図10と11にグラフとして示される。これは5秒または30秒のエクステンション時間でのダイソミーとトリソミーのサンプルについて染色体21に整列した合計のリードのパーセンテージを示す。
図10の左側は、5秒のエクステンション時間のパーセンテージを示す。また、右は30秒のエクステンション時間の場合を示す。実験は2回行なわれた。図11は、個々のサンプルの平均および2回の間での違いを示す。
どちらも10%のT21スパイク−インが、D21スパイク−インより整列したリードの高い割合を有することが見られる。図11は、ダイソミーとトリソミーの間の差が30秒より5秒のエクステンション時間の方がわずかに大きいことを示す。これは、T21リード/D21リードの比率の計算により示される(下記の表15を参照、それは5秒か30秒のエクステンション・ステップでのトリソミー/ダイソミー・サンプルの比率の比率を示す)。
ダイソミーとトリソミーのサンプルがそれぞれのエクステンション時間で2回実行されるので、それぞれのエクステンション時間についてトリ/ダイ用の4つの比率と1つのダイ/ダイ比率が得られる。表15は、5秒のエクステンション時間の平均が30秒のものより高いことを示す。この結果の可能な1つの解釈は、より長いフラグメントのエクステンションに対して、より短いエクステンション時間がバイアスを有することを示す。それぞれの結果および平均は、図12および13にグラフとしてに示される。図14は、染色体21の比率を他の染色体と比較し、常染色体のすべてについてトリソミー/ダイソミー・サンプルおよびダイソミー/ダイソミー・サンプルの平均比率を示す。図14は、常染色体のすべてについて、1つの染色体当たり4つのトリソミー/ダイソミー比率のデータポイントの平均と、1つのダイソミー/ダイソミー比率のデータポイントを示す。染色体21用のデータは黒いバーで示される。データは、1と本質的に異なるただ一つの比率は、両方のエクステンション時間(点線)における染色体21用のトリソミー/ダイソミー比率であることを明白に示す。
結論
10%のトリソミー21スパイクインは、5秒および30秒のエクステンション時間の両方で先に述べたアダプターPCR増幅方法を使用して、ゲノムワイド次世代シークエンシングによって検知された。2つのダイソミーおよび2つのトリソミー・サンプルは、4つのトリソミー/ダイソミー比率の計算を許容した。オーバーラップがあるが、4つのT21/D21比率の平均は、30秒のエクステンションより5秒のエクステンションにおいて13%高い。この結果の可能な1つの解釈は、より長いフラグメントの増幅に対して、より短いエクステンション時間の方がバイアスがかかるということである。常染色体がすべて分析される場合、データが1と本質的に異なるただ一つの比率が、染色体21についてのトリソミー/ダイソミー比率であることを示す。
例10 − ストゥファー法を使用するデータ分析
107のサンプル(コントロールを含む)が集められた。そのうち24はサンプルの明白な溶血により除外された。試験サンプルは、トリソミー21胎児を持っている危険の高い妊娠している女性からの血液サンプルだった(コントロールは妊娠していない人々からの血漿だった)。その後、DNAはQIAmp Circulating Nucleic Acidキット(Qiagen)を使用してサンプルから抽出された。また、DNAは例2に述べられているようなアダプター−PCR方法を使用して増幅された。増幅に続いて、DNAは標識付けされ、標識が付けられた参照DNA(妊娠していない女性からのゲノムDNA)と一緒にマイクロアレイでハイブリダイズされた。その後、アレイは洗われ走査された。また、各スポットでの蛍光の量は例2に述べられているような特徴抽出ソフトウェアを使用して計られた。その後、Cy3/Cy5比率は以下のように計算され分析された。
各サンプルについては、ペアでウィルコクソンの順位和検定が、ボンフェローニ(Bonferroni)複合テスト補正を使用して、染色体位置(染色体Xを除き)によってグループ化されたサンプルログ−比率(sample log-ratios)間で行なわれた(Hollander & Wolfe (1999) Non Parametric Statistical Methods. New York: John Wiley & Sons)。その後、染色体21について得られたp−値はストゥファー法(Stouffer et al. (1949) The American Soldier, Volume 1: Adjustment During Army Life. Princeton University Press, Princeton)を使用して組み合わされた。これは単一のp−値に帰着する。残りの83のサンプル(11トリソミー21妊娠を含む)のうち、65は真の陰性として呼ばれた。9は真の陽性だった。また、9は高いDLR(比率における高い変化)により判断できなかった。単一のp値の典型的なカットオフ値はサンプルまたはアレイの質スコアと組み合わされて0.01だった。
上記の例は発明の好ましい実施態様を説明するために提供される。当業者は、修正が請求項の範囲から外れずになされてもよいことを認識するだろう。この特許出願の記述の長さを合理的にするために、特徴のあらゆる組み合わせの徹底的な記述は提供されていない。
記述された実施態様および従属クレームのオプションの好ましい特徴および修正はすべて、ここに教示された発明のすべての態様において使用可能である。更に、すべてのオプションの好ましい特徴および記述された実施態様の修正と同様に従属クレームの個々の特徴は、互いに結合可能で交換可能である。本用途の主題の前出生の診断方法に関して記述された特徴はすべて腫瘍の識別のために適用されてもよく、また逆の場合もある。当業者は、明示的にここに言及されない他の腫瘍および胎児の条件に当てはまるために開示をここに適応させることができるだろう。本出願は優先権を主張する英国特許出願1404063.8、および本出願に伴う要約での開示は、参照されここに組込まれる。

Claims (39)

  1. 循環セルフリー核酸のサンプル内の第1のソースに由来する対象シークエンスを有する核酸の量の増減の有無を決定する方法であって、該サンプルは第1のソースからの核酸と第2のソースからの核酸を含み、
    セルフリー核酸のサンプルはより短いサイズ範囲およびより長いサイズ範囲の2つのサイズ範囲の核酸を含み、第1のソースからの核酸はより短いサイズ範囲にある、以下を含む方法:
    a)少なくともより短いサイズ範囲にある核酸から増幅用テンプレートを形成すること;
    b)テンプレートを増幅することにより、より短いサイズ範囲の核酸についてサンプルを豊化し、対象シークエンスを有する核酸のヌクレオチドシークエンスと無関係の方法で増幅生産物を形成すること;および
    c)サンプルに由来するコントロール増幅生成物の量と比較することにより、対象シークエンスを有する核酸の増幅生成物の相対量を決定すること、ここで該コントロール増幅生成物はコントロール・シークエンスを有し、コントロール・シークエンスは対象シークエンスとは異なっている。
  2. 増幅生産物は、約750bp未満のサイズを持っている、請求項1記載の方法。
  3. 増幅生産物は、約250bp未満のサイズを持っている、請求項1記載の方法。
  4. 増幅生産物は、約140−180bpのサイズを持っている、請求項1から3のいずれか1項記載の方法。
  5. より短いサイズの範囲内にある最長の核酸が、より長いサイズ範囲内にある最短の核酸の1/2〜1/3倍である、請求項1から4のいずれか1項記載の方法。
  6. ステップa)はサンプル中の核酸の端部へ、プライマー結合部位を含むアダプター分子をリゲートすることを含む、請求項1から5のいずれか1項記載の方法。
  7. ステップb)はアダプター・プライマー結合部位に相補的なプライマーをアニーリングし、ポリメラーゼ連鎖反応を使用してテンプレートを増幅することを含む、請求項6記載の方法。
  8. アニーリング/エクステンション時間が30秒以下である、請求項7記載の方法。
  9. アニーリング/エクステンション時間が30秒未満である、請求項7記載の方法。
  10. アニーリング/エクステンション時間が約20秒または20秒未満である、請求項7記載の方法。
  11. アニーリング/エクステンション時間が約10秒または10秒未満である、請求項7記載の方法。
  12. アニーリング/エクステンション時間が約5秒または5秒未満である、請求項7記載の方法。
  13. コントロール・シークエンスは対象の核酸シークエンスとは異なる染色体からのものである、請求項1から12のいずれか1項記載の方法。
  14. 対象シークエンスは多形性の疑いがない、請求項1から13のいずれか1項記載の方法。
  15. コントロール・シークエンスは多形性の疑いがない、請求項1から14のいずれか1項記載の方法。
  16. ステップc)は前記増幅生成物の量と、第2のソースからの核酸に対応する1セットの参照核酸から得られた増幅生成物と比較することを含み、該参照核酸は対象の核酸シークエンスに対応するシークエンスを有する核酸を含み;
    該参照核酸は、対応するシークエンスを有する核酸の増加したかまたは減少した量を含んでいない、請求項1から15のいずれか1項記載の方法。
  17. コントロール・シークエンスと比較して、対象シークエンスから得られた増幅生成物の量の比率を決定することを含み、該比率を参照核酸のセットから得られた当量比と比較する、請求項16記載の方法。
  18. ステップc)は以下を含む、請求項16または17記載の方法:
    i) 増幅生産物を、対象核酸シークエンスに本質的に相補的なテスト・プローブと少なくとも1つのコントロールプローブとハイブリダイズすること、ここでコントロールプローブはコントロールシークエンスに本質的に相補的である;
    ii) 参照核酸を、テスト・プローブおよびコントロールプローブとコ−ハイブリダイズすること;
    iii) テスト・プローブおよびコントロール・プローブとの、増幅生産物とのハイブリダイゼーションの量と参照核酸のハイブリダイゼーションの量とを比較し、サンプル/参考試験ハイブリダイゼーション比率およびサンプル/参照コントロール・ハイブリダイゼーション比率を得ること;および
    iv) テスト・ハイブリダイゼーション比率をコントロール・ハイブリダイゼーション比率と比較し、第1のソース内の対象の核酸シークエンスを有する核酸の量が増加したか減少したかを判断すること。
  19. プローブは固体担体上に提供される、請求項18記載の方法。
  20. 固体担体はチップである、請求項19記載の方法。
  21. 固体担体は複数のビーズである、請求項19記載の方法。
  22. ビーズはそれぞれ単一のプローブ・シークエンスを含んでいる、請求項21記載の方法。
  23. 固体担体は多孔性表面を有する、請求項19から21のいずれか1項記載の方法。
  24. 次世代シーケンシング方法を使用して、対象の核酸シークエンスの増幅生成物の相対量を決定することを含む、請求項1から15のいずれか1項記載の方法。
  25. 定量PCR方法を使用して、対象の核酸シークエンスの増幅生成物の相対量を決定することを含んでいる、請求項1から15のいずれか1項記載の方法。
  26. 定量PCR方法はディジタルPCRを含んでいる、請求項25記載の方法。
  27. 被検者の異常の可能性を予測する方法であって、請求項1から26のいずれか1項記載の方法を使用して、対象シークエンスを有する核酸の量の増減の有無を決定することを含み、対象シークエンスを有する核酸の量の増減があめ場合には、被検者は異常の増加した危険を有すると予測される方法。
  28. ステップc)は、前記増幅生成物の量と、第2のソースからの核酸に対応する1セットの参照核酸から得られた増幅生成物と比較することを含み、該参照核酸は対象の核酸シークエンスに対応するシークエンスを有する核酸を含み、
    該参照核酸は、対応するシークエンスを有する核酸の増減した量を含んでいない、
  29. 参照核酸は被検者から得られた、請求項28記載の方法。
  30. 参照核酸は、第1のソースからの核酸を含むことが予想されない被検者からの細胞から得られた、請求項28または29記載の方法。
  31. 参照核酸は、被検者の白血球あるいは被検者の頬側細胞から得られた、請求項28から30のいずれか1項記載の方法。
  32. 参照核酸はヌクレオソーマルDNAから抽出された、請求項28から31のいずれか1項記載の方法。
  33. 被検者は妊娠している女性であり、異常は胎児のコピー数欠損であり、第1のソースからの核酸は胎児の核酸であり、第2のソースからの核酸は母体の核酸である、請求項27から32のいずれか1項記載の方法。
  34. 対象の核酸シークエンスは、染色体13、染色体18 、染色体21および/またはX染色体である、請求項1から33のいずれか1項記載の方法。
  35. コントロール・シークエンスは染色体13、染色体18あるいは染色体21でない常染色体からのものである、請求項1から34のいずれか1項記載の方法。
  36. 方法は、パトー症候群、エドワーズ症候群あるいはダウン症候群あるいはターナー症候群の非侵襲性の出生前診断方法である、請求項1から35のいずれか1項記載の方法。
  37. 被検者は腫瘍を有する疑いがあり、第1のソースは腫瘍細胞であり、第2のソースは非腫瘍細胞である、請求項1から36のいずれか1項記載の方法。
  38. コントロール・シークエンスが既知の良性のコピー数変化を含まないように選択される、請求項1から37のいずれか1項記載の方法。
  39. ステップc)で得られるデータは、サンプル内の特異的プローブ・セットと同じサンプルに由来するコントロールと比較するためにウィルコクソン検定を使用して分析される、請求項1から38のいずれか1項記載の方法。
JP2016556745A 2014-03-07 2015-03-09 対象シークエンスを有する核酸の量の減少または増加の検出方法 Pending JP2017506908A (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
GB1404063.8A GB2524948A (en) 2014-03-07 2014-03-07 Detecting Increase or Decrease in the Amount of a Nucleic Acid having a Sequence of Interest
GB1404063.8 2014-03-07
PCT/GB2015/050680 WO2015132614A1 (en) 2014-03-07 2015-03-09 Detecting increase or decrease in the amount of a nucleic acid having a sequence of interest

Publications (1)

Publication Number Publication Date
JP2017506908A true JP2017506908A (ja) 2017-03-16

Family

ID=50554701

Family Applications (2)

Application Number Title Priority Date Filing Date
JP2016556716A Pending JP2017506907A (ja) 2014-03-07 2015-03-09 循環セルフリー核酸のサンプルの富化
JP2016556745A Pending JP2017506908A (ja) 2014-03-07 2015-03-09 対象シークエンスを有する核酸の量の減少または増加の検出方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
JP2016556716A Pending JP2017506907A (ja) 2014-03-07 2015-03-09 循環セルフリー核酸のサンプルの富化

Country Status (5)

Country Link
US (2) US20170016054A1 (ja)
EP (2) EP3114236A1 (ja)
JP (2) JP2017506907A (ja)
GB (1) GB2524948A (ja)
WO (2) WO2015132614A1 (ja)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9892230B2 (en) 2012-03-08 2018-02-13 The Chinese University Of Hong Kong Size-based analysis of fetal or tumor DNA fraction in plasma
TWI813141B (zh) 2014-07-18 2023-08-21 香港中文大學 Dna混合物中之組織甲基化模式分析
JP6432261B2 (ja) 2014-09-30 2018-12-05 ブラザー工業株式会社 液体消費装置
US10364467B2 (en) 2015-01-13 2019-07-30 The Chinese University Of Hong Kong Using size and number aberrations in plasma DNA for detecting cancer
US11514289B1 (en) * 2016-03-09 2022-11-29 Freenome Holdings, Inc. Generating machine learning models using genetic data
US20190203200A1 (en) * 2016-09-14 2019-07-04 Toray Industries, Inc. Method of collecting cell-free dna
JP6844833B2 (ja) * 2016-11-02 2021-03-17 学校法人日本医科大学 絞扼性腸閉塞の術前診断補助方法
WO2018099418A1 (en) 2016-11-30 2018-06-07 The Chinese University Of Hong Kong Analysis of cell-free dna in urine and other samples
CN107988344A (zh) * 2017-11-23 2018-05-04 厦门思诺恩生物工程有限公司 一种基于多孔滤膜的核酸测序方法和装置
KR20200103666A (ko) 2017-12-27 2020-09-02 도레이 카부시키가이샤 핵산의 회수 방법
US11211147B2 (en) 2020-02-18 2021-12-28 Tempus Labs, Inc. Estimation of circulating tumor fraction using off-target reads of targeted-panel sequencing
US11211144B2 (en) 2020-02-18 2021-12-28 Tempus Labs, Inc. Methods and systems for refining copy number variation in a liquid biopsy assay
US11475981B2 (en) 2020-02-18 2022-10-18 Tempus Labs, Inc. Methods and systems for dynamic variant thresholding in a liquid biopsy assay

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005160470A (ja) * 2003-10-16 2005-06-23 Sinuhe Hahn 胎児遺伝形質の非侵襲的検出
JP2008522590A (ja) * 2004-12-02 2008-07-03 エピゲノミックス アクチェンゲゼルシャフト 前立腺細胞増殖障害の予後に関連する遺伝子発現分析のための方法および核酸
JP2009504168A (ja) * 2005-08-17 2009-02-05 メディクシス エス.エー. Ck19発現を確認する組成物及び方法
JP2010162033A (ja) * 2000-12-01 2010-07-29 Response Genetics Inc 上皮増殖因子受容体及びHER2−neu遺伝子発現並びにそれらのレベルと生存との相関関係を決定する方法
US20100285537A1 (en) * 2009-04-02 2010-11-11 Fluidigm Corporation Selective tagging of short nucleic acid fragments and selective protection of target sequences from degradation
JP2011521631A (ja) * 2008-05-30 2011-07-28 キアゲン ゲゼルシャフト ミット ベシュレンクテル ハフツング 核酸を単離する方法
US20110294699A1 (en) * 2007-08-29 2011-12-01 Sequenom, Inc. Methods and compositions for universal size-specific pcr
GB2488358A (en) * 2011-02-25 2012-08-29 Univ Plymouth Enrichment of foetal DNA in maternal plasma
US20130123120A1 (en) * 2010-05-18 2013-05-16 Natera, Inc. Highly Multiplex PCR Methods and Compositions
US20130237431A1 (en) * 2012-03-08 2013-09-12 The Chinese University Of Hong Kong Size-based analysis of fetal dna fraction in maternal plasma
WO2013150503A1 (en) * 2012-04-06 2013-10-10 The Chinese University Of Hong Kong Noninvasive prenatal diagnosis of fetal trisomy by allelic ratio analysis using targeted massively parallel sequencing
JP2014502845A (ja) * 2010-12-22 2014-02-06 ナテラ, インコーポレイテッド 非侵襲性出生前親子鑑定法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005514956A (ja) * 2002-01-18 2005-05-26 ジェンザイム・コーポレーション 胎児dnaの検出および対立遺伝子の定量化のための方法
AU2007260750A1 (en) * 2006-06-16 2007-12-21 Sequenom, Inc. Methods and compositions for the amplification, detection and quantification of nucleic acid from a sample
EP2245195B1 (en) * 2008-01-30 2012-06-20 Novartis AG Two stage enrichment of cell-free fetal dna in maternal plasma
WO2011082386A1 (en) * 2009-12-31 2011-07-07 The Trustees Of Columbia University In The City Of New York Specific amplification of fetal dna sequences from a mixed, fetal-maternal source
US20110312503A1 (en) * 2010-01-23 2011-12-22 Artemis Health, Inc. Methods of fetal abnormality detection
US10131947B2 (en) * 2011-01-25 2018-11-20 Ariosa Diagnostics, Inc. Noninvasive detection of fetal aneuploidy in egg donor pregnancies
US20120190020A1 (en) * 2011-01-25 2012-07-26 Aria Diagnostics, Inc. Detection of genetic abnormalities

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010162033A (ja) * 2000-12-01 2010-07-29 Response Genetics Inc 上皮増殖因子受容体及びHER2−neu遺伝子発現並びにそれらのレベルと生存との相関関係を決定する方法
JP2005160470A (ja) * 2003-10-16 2005-06-23 Sinuhe Hahn 胎児遺伝形質の非侵襲的検出
JP2008522590A (ja) * 2004-12-02 2008-07-03 エピゲノミックス アクチェンゲゼルシャフト 前立腺細胞増殖障害の予後に関連する遺伝子発現分析のための方法および核酸
JP2009504168A (ja) * 2005-08-17 2009-02-05 メディクシス エス.エー. Ck19発現を確認する組成物及び方法
US20110294699A1 (en) * 2007-08-29 2011-12-01 Sequenom, Inc. Methods and compositions for universal size-specific pcr
JP2011521631A (ja) * 2008-05-30 2011-07-28 キアゲン ゲゼルシャフト ミット ベシュレンクテル ハフツング 核酸を単離する方法
US20100285537A1 (en) * 2009-04-02 2010-11-11 Fluidigm Corporation Selective tagging of short nucleic acid fragments and selective protection of target sequences from degradation
US20130123120A1 (en) * 2010-05-18 2013-05-16 Natera, Inc. Highly Multiplex PCR Methods and Compositions
JP2014502845A (ja) * 2010-12-22 2014-02-06 ナテラ, インコーポレイテッド 非侵襲性出生前親子鑑定法
GB2488358A (en) * 2011-02-25 2012-08-29 Univ Plymouth Enrichment of foetal DNA in maternal plasma
US20130237431A1 (en) * 2012-03-08 2013-09-12 The Chinese University Of Hong Kong Size-based analysis of fetal dna fraction in maternal plasma
WO2013150503A1 (en) * 2012-04-06 2013-10-10 The Chinese University Of Hong Kong Noninvasive prenatal diagnosis of fetal trisomy by allelic ratio analysis using targeted massively parallel sequencing

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
FETAL DIAGNOSIS AND THERAPY, vol. 25, no. 5, JPN6019000216, 2008, pages 314 - 319, ISSN: 0004087679 *

Also Published As

Publication number Publication date
EP3114236A1 (en) 2017-01-11
GB2524948A (en) 2015-10-14
EP3114237A1 (en) 2017-01-11
WO2015132614A1 (en) 2015-09-11
GB201404063D0 (en) 2014-04-23
JP2017506907A (ja) 2017-03-16
WO2015132615A1 (en) 2015-09-11
US20170016047A1 (en) 2017-01-19
US20170016054A1 (en) 2017-01-19

Similar Documents

Publication Publication Date Title
JP2017506908A (ja) 対象シークエンスを有する核酸の量の減少または増加の検出方法
JP6585117B2 (ja) 胎児の染色体異数性の診断
TWI817187B (zh) 偵測突變以用於癌症篩選分析
EP3117011B1 (en) Methods and processes for non-invasive assessment of genetic variations
JP6830094B2 (ja) 染色体異常を検出するための核酸及び方法
CN112037860B (zh) 用于非入侵性性染色体非整倍性确定的统计分析
JP6328934B2 (ja) 非侵襲性出生前親子鑑定法
US20150275290A1 (en) Non-invasive method for detecting a fetal chromosomal aneuploidy
JP6793112B2 (ja) 胎児コピー数変異の統計的尤度を提供するアッセイ方法および胎児染色体異数性の尤度決定のためのアッセイ方法
US20140186827A1 (en) Assays for the detection of genotype, mutations, and/or aneuploidy
EP3018213A1 (en) Method for determining the presence of a biological condition by determining total and relative amounts of two different nucleic acids
JP2021520809A (ja) がんまたは新生物評価のための組成物および方法
TW202144583A (zh) 血漿粒線體dna分析之應用
WO2020020174A1 (en) Cell-free dna damage analysis and its clinical applications
CA2986200A1 (en) Multiplexed parallel analysis of targeted genomic regions for non-invasive prenatal testing
JP2023537215A (ja) 無細胞核酸についてのヌクレアーゼに関連する末端シグネチャ分析
CN112888783A (zh) 改善游离dna质量
JP2009518051A (ja) 異なる個体の核酸を含む生体試料から遺伝子型を決定する方法
EP3596229A1 (en) Method and system for nucleic acid sequencing
Class et al. Patent application title: ASSAYS FOR THE DETECTION OF GENOTYPE, MUTATIONS, AND/OR ANEUPLOIDY

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180227

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20181219

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190115

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20190806