JP6621802B2 - 遺伝的多様体を検出する方法 - Google Patents

Description

発明の分野
本発明は、ＤＮＡ試料中の遺伝的多様体を検出する方法に関する。

発明の背景
遺伝的多様体は、所与の領域の参照ＤＮＡ配列と異なる１つまたは複数のヌクレオチドである。例えば遺伝的多様体は、１つまたは複数のヌクレオチドの欠失、置換または挿入を含み得る。

目的の領域のＤＮＡ配列を決定すること、および決定された配列を参照配列と比較すること、によって、ＤＮＡ試料を既知の遺伝的多様体について分析でき、または目的の領域にある過去に知られていない遺伝的多様体を発見できる。

ＤＮＡシーケンシングは、Ｍｅｔｚｋｅｒ， Ｍ．Ｌ．， Ｎａｔ Ｒｅｖ Ｇｅｎｅｔ ２０１０ Ｊａｎ；１１（１）： ３１−４６によって論評された、古典的な連鎖停止法、または複数のハイスループット次世代シーケンシング（ＮＧＳ）法の１つなど、様々な技術を利用して実施できる。

Ｉｌｌｕｍｉｎａシーケンシング、４５４パイロシーケンシング、Ｈｅｌｉｓｃｏｐｅ一分子シーケンシング、一分子リアルタイム（ＳＭＲＴ）シーケンシングおよびイオン半導体シーケンシングプラットフォームは、「合成によるシーケンシング」の原理に基づくＤＮＡシーケンシング法の例である。これらの方法において、ヌクレオチド塩基が新たに合成された相補鎖に組み入れられると放出されるシグナルの検出を通して、ＤＮＡのテンプレート鎖の配列が決定される。

ＤＮＡシーケンシングプラットフォームは、誤り率を有する。例えば時折、増幅反応に用いられるポリメラーゼは、合成される相補鎖に誤ったヌクレオチド塩基を組み入れて、ＤＮＡテンプレート内のその位置でのヌクレオチドの決定を間違えることがある。ＮＧＳ法の検出限界は、ライブラリーの調製（通常はＰＣＲによる増幅を含む）およびシーケンシングそのもの、の２段階での誤りによって定義される。

このことは、低頻度で、例えば用いられるシーケンシング法の誤り率に近い頻度またはそれ未満の頻度でＤＮＡ試料中に存在するに過ぎない遺伝的多様体の検出の場合には特に問題になる。そのような状況下では、同定された遺伝的多様体が事実であるか（即ち、ＤＮＡテンプレート分子中に実際に存在する）、または誤りであるかを決定することが、困難または不可能である。

Ｉｌｌｕｍｉｎａシーケンシングの場合、バックグランド誤り率が、異なる遺伝的多様体およびゲノム位置によって変動し、大きな変動性を有する。それゆえ、約１％以下の頻度でＤＮＡ試料中に存在する突然変異を検出することが、問題になる。

既存のＤＮＡシーケンシング法および遺伝的多様体同定法は、特に少量のＤＮＡを有する試料では、複数の領域で希少な新規多様体を検出することに関して限界を有する。

方法は、用いられた方法の誤り率（即ち、バックグランドノイズ）よりも低い、またはそれと類似の頻度で生じる突然変異を同定するのは、典型的には不可能である。

デジタルＰＣＲ（ｄＰＣＲ； Ｖｏｇｅｌｓｔｅｉｎ ｂ．， Ｋｉｎｚｌｅｒ Ｋ．Ｗ． Ｐｒｏｃ． Ｎａｔｌ． Ａｃａｄ． Ｓｃｉ． Ｕ．Ｓ．Ａ． １９９９ ９６（１６）：９２３６−４１； Ｓｙｋｅｓ， Ｐ．Ｊ． ｅｔ ａｌ．， ＢｉｏＴｅｃｈｎｉｑｕｅｓ １９９２ １３（３）： ４４４−９）は、特定の多様体を検出するように設計されたプライマーおよびアッセイの利用を必要とするため、新規な（即ち、過去に同定されていない）遺伝的多様体の同定には有用でない。その上、ｄＰＣＲは、特にＤＮＡ試料が限定される場合に、目的の複数の領域を並行して分析するのに限定された範囲を有する。

Ｓａｆｅ−ＳｅｑＳおよび一分子での分子反転プローブなど、ＤＮＡの単一プールからの単一ＤＮＡ分子を標識するための他の複雑な方法が、存在する（Ｋｉｎｄｅ Ｉ， ｅｔ ａｌ．， Ｐｒｏｃ Ｎａｔｌ Ａｃａｄ Ｓｃｉ Ｕ Ｓ Ａ． ２０１１ １０８（２３）： ９５３０−５； Ｈｉａｔｔ ＪＢ， ｅｔ ａｌ．， Ｇｅｎｏｍｅ Ｒｅｓ． ２０１３ ２３（５）：８４３−５４．）。これらの方法は、複数の遺伝子（即ち、該当する複数の領域）の同時分析およびＤＮＡが限られる場合には適さない。

複数の試験によって、癌ＤＮＡの非侵襲性検出が実証された（Ｄａｗｓｏｎ ＳＪ， ｅｔ ａｌ．， Ｎ Ｅｎｇｌ Ｊ Ｍｅｄ． ２０１３ ３６８（１３）：１１９９−２０９； Ｆｏｒｓｈｅｗ Ｔ， ｅｔ ａｌ．， Ｓｃｉ Ｔｒａｎｓｌ Ｍｅｄ． ２０１２ ４（１３６）：１３６ｒａ６８； Ｍｕｒｔａｚａ Ｍ， ｅｔ ａｌ． Ｎａｔｕｒｅ． ２０１３ ４９７（７４４７）：１０８−１２）。しかし、（ａ）関連の癌突然変異を検出するための十分な塩基のゲノムのスクリーニング、（ｂ）そのような突然変異に関する少量の断片化ＤＮＡのスクリーニング、および（ｃ）多くの「野生型」分子のうちの低頻度の突然変異体腫瘍ＤＮＡ分子の検出など、この分野には大きな難題が依然として存在する。

例えば、Ｆｏｒｓｈｅｗ Ｔ， ｅｔ ａｌ．， Ｓｃｉ Ｔｒａｎｓｌ Ｍｅｄ． ２０１２ ４（１３６）：１３６ｒａ６８には、血中の癌突然変異に関する大きな領域のゲノムをスクリーニングすることが記載されているが、この方法の検出限界は、約１％〜２％のアレル頻度（ＡＦ）であった。

発明の概要
本発明は、上記の問題の解決策を提供する。実際の遺伝的多様体は、方法の誤り（即ち、ＤＮＡ増幅およびシーケンシングプラットフォームの誤り率）を考慮して各位置での遺伝的多様体のバックグランド頻度を決定することにより同定される。

第一の態様において、本発明は、ＤＮＡ試料において目的の領域中の遺伝的多様体を検出するための方法であって、
（ｉ）所与のシーケンシングプラットフォーム、シーケンシング工程およびシーケンシング深度について、増幅およびシーケンシング誤差により増幅反応物のシーケンシング結果で観察されると予測される遺伝的多様体または複数の多様体を裏づけるリード数の分布（リードカウント分布）を決定すること；
（ｉｉ）ステップ（ｉ）において決定されたリードカウント分布に基づいて、それ以上で増幅反応物のシーケンシング結果において遺伝的多様体が観察されるべき閾値頻度を確定して所与の増幅反応物における遺伝的多様体の存在について陽性判定を割り付けること；
（ｉｉｉ）反復増幅反応物における目的の領域の増幅可能なテンプレート分子の平均数がステップ（ｉｉ）において決定された閾値頻度の逆数よりも少なくなるようにＤＮＡ試料を複数の反復増幅反応物に分配すること；
（ｉｖ）ステップ（ｉｉｉ）の増幅反応を実施すること、および増幅反応の産物をシーケンシングすること；
（ｖ）ステップ（ｉｉ）およびステップ（ｉｖ）の結果に基づいて、各反復増幅反応物における遺伝的多様体の有無を決定すること；ならびに
（ｖｉ）（ｖ）の結果を統合してＤＮＡ試料における目的の領域の遺伝的多様体の有無を決定すること、
を含む、方法を提供する。

本発明は、ＤＮＡ試料において目的の領域の遺伝的多様体を検出する方法であって、
（ｉ）所与のシーケンシングプラットフォームについて、増幅およびシーケンシング誤差により増幅反応物のシーケンシング結果で遺伝的多様体または複数の遺伝的多様体が観察されると予測される平均頻度および頻度の分散を決定すること；
（ｉｉ）ステップ（ｉ）で決定された平均頻度および頻度の分散に基づいて、それ以上で増幅反応物のシーケンシング結果において遺伝的多様体が観察されるべき閾値頻度を確定して所与の増幅反応物における遺伝的多様体の存在について陽性判定を割り付けること；
（ｉｉｉ）反復増幅反応物における目的の領域の増幅可能なテンプレート分子の平均数がステップ（ｉｉ）において決定された閾値頻度の逆数よりも少なくなるようにＤＮＡ試料を複数の反復増幅反応物に分配すること；
（ｉｖ）ステップ（ｉｉｉ）の増幅反応を実施すること、および増幅反応の産物をシーケンシングすること；
（ｖ）ステップ（ｉｉ）およびステップ（ｉｖ）の結果に基づいて、各反復増幅反応物における遺伝的多様体の有無を決定すること；ならびに
（ｖｉ）（ｖ）の結果を統合してＤＮＡ試料における目的の領域の遺伝的多様体の有無を決定すること、
を含む、方法を提供する。

有利には、この方法は、ＤＮＡ試料において非常に低頻度の遺伝的多様体の検出を可能にする。したがって方法は、過去に知られた突然変異検出方法と比較して、遺伝的突然変異の存在（即ち、突然変異がより低頻度で存在する場合）のより早期の同定（例えば、癌関連の突然変異など疾患病態に関連して）を可能にする。それゆえ方法は、腫瘍の再出現、腫瘍の成長、薬物耐性の発達／腫瘍の進展に関するスクリーニングおよび治療的に実行可能な突然変異の同定など、様々な適用例に用途を見出す。この方法はまた、改善された統計学的信頼性を有する低頻度の遺伝的突然変異の存在の同定を可能にする。

本発明の方法の重要な特色は、平均で、反復増幅反応物に存在する目的の領域の増幅可能なテンプレート分子の数が、陽性判定の閾値頻度の逆数よりも少なくなるように、ＤＮＡ試料を複数の反復増幅反応物に分配すること（例えば、希釈およびウェルへの分取により）である。こうして、遺伝的多様体が存在する増幅可能なテンプレート領域を有する各増幅反応物について、多様体の存在を決定するための閾値頻度よりも大きな頻度で多様体が観察される確率が高くなる。それによりこの方法は、ＤＮＡ試料中に非常に低頻度で存在する場合であっても、ＤＮＡ試料における遺伝的多様体の存在の検出を可能にする。例えばこの方法は、１％未満の頻度の検出を可能にする。

その上、本発明の方法は、遺伝的多様体の発見に有用である。即ち、目的の領域について参照ＤＮＡ配列からの３つの潜在的塩基置換のそれぞれのバックグランドレベルを即座に決定でき、バックグランドレベル（特定のＤＮＡ試料の）を超える遺伝的多様体の頻度を、この方法を用いて続いて同定できる。したがって、任意の特定の遺伝的多様体についてスクリーニングせずに、目的の領域の任意の位置にある新たな遺伝的多様体が、同定され得る。この特色は、診断、予知、突然変異発見、処置への応答のモニタリングおよび薬物耐性など、広範囲の分野における明確な利点に関連する。

幾つかの実施形態において、ステップ（ｖｉ）においてＤＮＡ試料の目的の領域に遺伝的多様体が存在すると決定されるためには、ステップ（ｖ）における遺伝的多様体の存在に関する陽性判定が１つよりも多くの復増幅反応物でなされなければならない。幾つかの実施形態において、遺伝的多様体の存在に関する陽性判定は、少なくとも３つの反復増幅反応物でなされなければならない。

幾つかの実施形態において、遺伝的多様体は、一ヌクレオチド多様体、即ち、同じ位置における１つのヌクレオチドから異なるヌクレオチドへの置換であり得る。幾つかの実施形態において、遺伝的多様体は、ヌクレオチドを付加または除去する挿入または欠失であり得る。幾つかの実施形態において、遺伝的多様体は、一ヌクレオチド多様体と挿入および／または欠失とを含む複数のイベントの組み合わせであり得る。幾つかの実施形態において、遺伝的多様体は、目的の異なる領域に存在する複数の遺伝的多様体で構成され得る。

複数の反復増幅反応物における遺伝的多様体の陽性判定を必要とすることで、ＤＮＡ試料に存在する遺伝的多様体の偽陽性判定の確率が低下する。反復増幅反応物において複数の陽性判定を必要とするこの方法はそれゆえ、遺伝的多様体の検出に関してより高い特異度を有する。

ＤＮＡ試料をシーケンシングするために用いられた方法における誤差の結果として所与の多様体が観察される（即ち、シーケンシングの結果において）平均頻度および変動係数（ＣＶ）を用いて、遺伝的多様体のバックグランドレベル（即ち、ノイズ）を決定および／またはモデル化できる。これらの値は、例えば累積分布関数（ＣＤＦ）値を決定するため、そして／またはｚスコアを計算するために、用いられ得る。

次に、遺伝的多様体のバックグランドノイズの測定および／またはモデルをその後用いて閾値頻度を確定し、それを超えて遺伝的多様体が観察されて所与の増幅反応物中に存在する（陽性判定）と決定されなければならない。陽性判定では、多様体の頻度はバックグランドレベルでの平均頻度よりも高くなければならない。

本発明の方法の幾つかの実施形態において、ステップ（ｉｉ）の閾値頻度は、二項分布、過分散二項分布、ベータ分布、正規分布、指数分布またはガンマ確率分布モデルを用いて決定される。幾つかの実施形態において、それ以上で所与の遺伝的多様体が観察されて反復増幅反応物中に存在すると決定される閾値頻度は、その遺伝的多様体の累積分布関数（ＣＤＦ）値が０．９９、０．９９５、０．９９９、０．９９９９、０．９９９９９またはそれを超える所定の閾値（ＣＤＦ＿ｔｈｒｅｓｈ）に達する頻度である。

本発明の方法の幾つかの実施形態において、ステップ（ｉｉ）の閾値頻度は、ｚスコアカットオフを利用して決定される。ここで、ステップ（ｉ）で決定された遺伝的多様体のバックグランド平均頻度および頻度の分散は、正規分布でモデル化され、突然変異をコール（ｃａｌｌ）するための閾値頻度は、バックグランド平均頻度を超える多数の標準偏差であるｚスコアでの頻度である。幾つかの実施形態において、閾値頻度は、ｚスコアが２０での頻度である。幾つかの実施形態において、閾値頻度は、ｚスコアが３０での頻度である。

幾つかの実施形態において、それ以上で増幅反応物のシーケンシング結果において遺伝的多様体が観察されるべき閾値頻度を確定して所与の増幅反応物中の遺伝的多様体の存在について陽性判定を割り付けることは、
（ａ）ステップ（ｉ）において複数の遺伝的多様体について決定されたリードカウント分布（場合により、複数の遺伝的多様体について決定された平均頻度および頻度の分散によって定義された正規分布）に基づいて、それ以上で増幅反応物のシーケンシング結果において遺伝的多様体が観察されるべき複数の閾値頻度を確定して所与の増幅反応物中の遺伝的多様体の存在について陽性判定を割り付けること、および
（ｂ）ステップ（ａ）に基づいて、それ以上で所与の増幅反応物のシーケンシング結果において遺伝的多様体が観察されるべき全体的閾値頻度を確定してその増幅反応物中の遺伝的多様体の存在について陽性判定を割り付けることであり、全体的閾値頻度が、ステップ（ａ）で決定された閾値頻度の９０％、９５％、９７．５％、９９％またはそれを超える閾値頻度がこの値未満である閾値頻度であること、
を含む。こうして、閾値頻度が目的の領域の各位置の各可能な塩基について決定される必要はなく、複数の遺伝的多様体に基づく全体的閾値が、本発明の方法で用いられ得る。

幾つかの実施形態において、遺伝的多様体が反復増幅反応物の１つの増幅可能なテンプレート分子中に存在する場合に、陽性判定がその反復測定でなされる確率が０．９以上であるように、所与の反復増幅反応物中に存在する目的の領域の増幅可能なテンプレート分子の平均数が決定される。

これによって、多様体が実際に存在する場合に、所与の反復増幅反応物が遺伝的多様体について陰性であると誤って決定される確率が最小限になる。

幾つかの実施形態において、複数の反応物は全て、試料からの同量の開始テンプレート材料を有する。幾つかの実施形態において、異なる反応物は、試料からの異なる量のテンプレート材料を有する。これらの量は、試料中の遺伝的変異体のアレル頻度を推定する場合に考慮され得る。

幾つかの実施形態において、ステップ（ｉ）は、ＤＮＡ試料を複数回シーケンシングして、遺伝的多様体または複数の遺伝的多様体についてリードカウント分布を決定することを含み、リードカウント分布は場合により、遺伝的多様体または複数の遺伝的多様体の平均頻度および頻度の分散によって定義された正規分布である。

有利には、バックグランド誤り率が、この方法で目的の領域中の遺伝的多様体について経験的に決定される場合には、バックグランド誤り率の推定がより正確であり、それによりＤＮＡ試料中の遺伝的多様体の存在を検出する感度がより大きくなり、信頼度がより大きくなる。

幾つかの実施形態において、増幅およびシーケンシング誤差により増幅反応物のシーケンシング結果において遺伝的多様体または複数の遺伝的多様体が観察されると予測されるリードカウント分布（場合により、増幅およびシーケンシング誤差により増幅反応物のシーケンシング結果において遺伝的多様体または複数の遺伝的多様体が観察されると予測される平均頻度および頻度の分散によって定義された正規分布）は、シーケンサーおよび／またはポリメラーゼの誤り率に基づいてステップ（ｉ）において決定される。幾つかの実施形態において、これは、配列の状況を考慮して決定される。

幾つかの実施形態において、ステップ（ｉ）は、データベース、チャート、表、リスト、カタログ、インデックス、ディレクトリ、またはレジスタにおいて参照値または複数の参照値を「探索すること」を含む。幾つかの実施形態において、参照値は、参照ＤＮＡ試料を複数回シーケンシングして、遺伝的多様体または複数の遺伝的多様体についてのリードカウント分布を決定することにより決定され、リードカウント分布は、場合により遺伝的多様体または複数の遺伝的多様体についての平均頻度および頻度の分散によって定義された正規分布である。幾つかの実施形態において、参照ＤＮＡ試料は、「一致した正常な」試料である。

有利には、本発明の方法は、各回でバックグランド誤り率を経験的に決定する必要なしに実施され得る。

幾つかの実施形態において、ステップ（ｉｉｉ）でのＤＮＡ試料の分配に続いて、各反復増幅反応物は、反復増幅反応物あたりに、目的の領域について１より多い増幅可能なテンプレート分子平均を有する。幾つかの実施形態において、反復増幅反応物あたりの増幅可能なテンプレート分子の平均数は、１より多く１０００未満、２より多く１０００未満、または５より多く１０００未満であろう。

有利には、ＤＮＡ試料は、テンプレート分子の分配に関連する方法の効率を実現するのに必要な量よりも多く分配する必要はない。これは、少なくともランニングコストおよび材料費を低減することに関する利点を有する。

幾つかの実施形態において、この方法は、ＤＮＡ試料の増幅可能なテンプレート分子の集団中に、２％未満、１％未満、０．５％未満、０．１％未満または０．０５％未満の頻度で存在する遺伝的多様体を検出できる。

幾つかの実施形態において、ＤＮＡ試料を複数の反復増幅反応物に分配することは、ＤＮＡ試料を希釈すること、および反復増幅反応物に分取すること、を含む。例えば増幅反応は、別々のウェル中で実施され得る。あるいは反復増幅反応物は、当業者に知られる他の手段により分配され得る。

本発明の方法の幾つかの実施形態において、反復増幅反応は、並行して実施され、増幅反応の産物のシーケンシングは、並行して実施される。

本発明の幾つかの実施形態において、この方法は、
（ｖｉｉ）ＤＮＡ試料における遺伝的多様体の頻度を決定すること、
をさらに含む。

有利には、これは、遺伝的多様体の頻度における経時的な変化（例えば、疾患経過および／または処置過程の間）および／または試料間の遺伝的多様体の頻度における差の決定を可能にする。

増幅反応は、一段階ＰＣＲまたは二段階ＰＣＲによって実施され得る。

幾つかの実施形態において、増幅反応は、目的の領域に隣接する１つまたは複数のプライマー対を使用して実施されて、試料および／または増幅反応反復測定物の特異的識別子配列を増幅産物中に組み込む。目的の結合した標的領域と共に読むと、識別子を利用して、どの試料および／または増幅反応複製物から標的配列が由来したかを同定できるように、識別子配列は、別の連続したＤＮＡ塩基と十分に異なる任意の連続したＤＮＡ塩基として定義され得る。用語「識別子配列」および「バーコード」は、本明細書では互換的に用いられる。

幾つかの実施形態において、プライマーは、配列アダプターを増幅産物に組み込む。幾つかの実施形態において、目的の領域に隣接するプライマーは、試料および／または増幅反応複製物特異的識別子配列、ならびに配列アダプターを含み、これらの識別子およびアダプター配列を一段階ＰＣＲの際に付加させる。幾つかの実施形態において、汎用の「標識」配列が、プライマー対の中に含まれる。標識配列は、ＰＣＲの次の回での標的として用いられ得る任意の連続したＤＮＡ塩基として定義され得る。複数のプライマーへの共通の標識配列の結合により、増幅後に、結合したこれらの共通の標識配列を有する複数の産物が得られる。最適な標識配列は、ゲノム配列の非特異的増幅を予防するために、目的のゲノムと十分異なる。幾つかの実施形態において、標識配列は、シーケンシングプライマーの結合部位などのさらなる特色を含み得る。幾つかの実施形態において、二回目のＰＣＲが、シーケンシングアダプターおよび場合によりさらなるバーコードを元々のＰＣＲ産物に結合させるために、標識配列、配列アダプターおよび場合によりさらなるバーコードを含むプライマーを使用して実施される。幾つかの実施形態において、ライゲーションを利用して、配列アダプターおよび場合によりさらなるバーコードを元々のＰＣＲ産物に結合させる。幾つかの実施形態において、複数の目的の領域が、並行して分析される。

それゆえ幾つかの例においてこの方法は、ハイスループット法として実行され得、遺伝的多様体に関する目的の複数の領域のスクリーニングを同時に可能にする。これは、速度、効率、ならびにランニングコストおよび材料費の低下に関して利点を有する。

本発明の方法の幾つかの実施形態において、リードカウント分布は、増幅およびシーケンシング誤差により増幅反応物のシーケンシング結果で遺伝的多様体または複数の遺伝的多様体が観察される、または観察されると予測される平均頻度および頻度の分散であるパラメータを特徴とする正規分布として定義され得る。

リードカウント分布は、任意のカウントで非参照アレルを観察する一般的確率セットとして定義され得る。

リードカウント分布は、増幅反応物のシーケンシング誤差により増幅反応物のシーケンシング結果で観察されると予測される遺伝的多様体または複数の遺伝的多様体を裏づけるリード数の分布であり得る。遺伝的多様体を裏づけるリードは、遺伝的多様体に関して陽性リードである。

例えば、リードカウント分布は、増幅およびシーケンシング誤差により、所与のシーケンシングプラットフォーム、シーケンシング工程およびシーケンシング深度について、遺伝的多様体または複数の遺伝的多様体を観察する一般的確率セットであり得る。

シーケンシング深度とも称されるリード深度は、特異的ゲノム位置（例えば、所与のヌクレオチド）がシーケンシング工程で読まれる回数として定義され得る。

本発明の方法は、広範囲の適用例に用途を見出す。実際、この方法は、任意の目的のために、任意のＤＮＡ試料中の目的の任意の領域にある遺伝的多様体の検出に有用である。

別の態様において、本発明は、試料中の腫瘍ＤＮＡを検出および／または定量するための先に記載された方法を提供する。幾つかの実施形態において、この方法は、試料中の循環する腫瘍ＤＮＡを検出および／または定量するためのものである。

本発明の任意の態様による幾つかの例において、試料は、対象から得られた生物学的試料である。幾つかの実施形態において、試料は、組織試料、例えば外科的試料である。幾つかの実施形態において、試料は、液体生検、例えば血液、血漿、尿、精液、糞便、痰、胸水、腹水、滑液、脳脊髄液、リンパ、乳頭液（ｎｉｐｐｌｅ ｆｌｕｉｄ）、または気管支洗浄液である。幾つかの実施形態において、試料は、細胞学的試料あるいは細胞材料を含有するスメアまたは流体、例えば子宮頚管スメア、鼻腔擦過物、またはスポンジ（細胞スポンジ）による食道採取物、内視鏡／胃内視鏡／結腸内視鏡生検または擦過物、子宮頚管粘膜または擦過物である。

上記試料の多くは非侵襲的に得ることができ、それゆえ対象への大きなリスクまたは不快感を伴わずに定期的に採取できる。

したがって、一態様において本発明は、診断もしくは予知または疾患のモニタリングのインビトロ方法における、先に記載された方法を提供する。

方法を利用して、ＤＮＡ試料を、易罹患性、所与の治療への耐性または応答に関連する、またはそれを予測する遺伝的多様性に関して分析できる。

したがって、別の態様において本発明は、疾患を発症する確率の高い（即ち、疾患への易罹患性が高い）対象を同定する方法であって、対象からＤＮＡ含有試料を得ること、および対象から得られたＤＮＡ含有試料で本発明の第一の態様による方法を実施すること、を含む、方法を提供する。幾つかの例において、本発明のこの態様によれば、ＤＮＡ含有試料は疾患関連の遺伝的多様体を含むと決定され、それにより疾患を発症する確率の高い対象（ＤＮＡ含有試料中で検出される疾患関連の遺伝的多様体を有さない対象と比較して）を同定する。

別の態様において、本発明は、対象からＤＮＡ含有試料を得ること、および対象から得られたＤＮＡ含有試料で本発明の第一の態様による方法を実施すること、を含む、診断の方法を提供する。幾つかの例において、本発明のこの態様によれば、ＤＮＡ含有試料は、疾患関連の遺伝的多様体を含むと決定される。

特定の易罹患性、所与の治療への耐性または応答を有し得る対象または複数の対象が、そのような易罹患性、耐性または応答を予測する遺伝的多様性について分析され得る。したがって、一態様において本発明は、易罹患性、治療への耐性または応答を予測する遺伝的多様性を同定する方法における先に記載された方法を提供する。

したがって本発明は、治療のための患者を選択する方法であって、患者からＤＮＡ含有試料を得ること、および患者から得られたＤＮＡ含有試料で本発明の第一の態様による方法を実施すること、を含む、方法を提供する。幾つかの例において、本発明のこの態様によれば、ＤＮＡ含有試料は、疾患関連の遺伝的多様体（例えば、リスクアレル）および／または治療への耐性もしくは応答を予測する遺伝的多様体を含むと決定される。この方法は、患者からのＤＮＡ含有試料が遺伝的多様体を含むという決定に基づいて、治療のための患者を選択することをさらに含み得る。幾つかの例において、この方法は、患者からのＤＮＡ含有試料が遺伝的多様体を含むという決定に基づいて、治療薬を投与するステップおよび／または治療薬の投与を推奨するステップをさらに含み得る。

本発明の方法はまた、疾患への易罹患性もしくは疾患の予後の、または治療への耐性もしくは応答の新規な（即ち、過去に同定されていない）決定子および／または予測子を同定するために用いることができる。

したがって、別の態様において本発明は、疾患を有する患者または複数の患者からＤＮＡ含有試料を得ること、および患者または複数の患者から得られたＤＮＡ含有試料で本発明の第一の態様による方法を実施すること、を含む、疾患関連の遺伝的多様体を同定する方法を提供する。幾つかの例において、本発明のこの態様によれば、ＤＮＡ含有試料または複数のＤＮＡ含有試料は、遺伝的多様体を含むと決定され、それにより遺伝的多様体が、疾患関連の遺伝的多様体であると同定される。この態様によれば、幾つかの実施形態において患者または複数の患者は、特定の臨床表現型および／または疾患予後を有し得る。したがって、幾つかの実施形態においてこの方法は、特定の臨床表現型および／または疾患予後の決定子または予測子を同定する。

別の態様において、本発明は、疾患の治療への耐性または応答を有する患者または複数の患者からＤＮＡ含有試料を得ること、および患者または複数の患者から得られたＤＮＡ含有試料で本発明の第一の態様による方法を実施すること、を含む、治療への耐性または応答に関連する遺伝的多様体を同定する方法を提供する。幾つかの例において、本発明のこの態様によれば、ＤＮＡ含有試料または複数のＤＮＡ含有試料は、遺伝的多様体を含むと決定され、それにより遺伝的多様体が、治療への耐性または応答に関連する遺伝的多様体であると同定される。

別の態様において、本発明は、治療への応答をモニタリングまたは評価する方法であって、本発明の第一の態様による方法が、治療の異なる段階で（例えば、介入前および介入時／介入後）患者から得られたＤＮＡ試料で実施される、方法を提供する。異なる時間に得られた試料中の遺伝的多様性の比較は、治療への応答の変化を明らかにし得る。例えば、治療時または治療後に得られた試料に存在し、治療前に得られた試料に存在しない遺伝的多様体、または遺伝的多様体の相対頻度は、例えば腫瘍の進展および／もしくは腫瘍量、または治療への他の生物学的応答を反映し得る。

別の態様において、遺伝的多様体の頻度、または異なる多様体の頻度の比率を利用して、治療への応答を予測、モニタリング、または評価できる。別の態様において、遺伝的多様体の存在もしくは頻度、または異なる多様体の頻度の比率を利用して、患者が処置されない場合、または患者が一連の治療の１つを与えられた場合の、患者のリスクレベルまたは予後を予測できる。

アンプリコンシーケンシング工程のバックグランドノイズ値の分布を示すグラフ。アンプリコンシーケンシングパネルで観察された異なる一ヌクレオチド多様体に関する（Ａ）平均、および（Ｂ）変動係数（ＣＶ、平均で割られた標準偏差として定義）。 複数の塩基置換についての累積分布を示すグラフであり、その累積分布関数は、データの経験的分布に対して比較した、それらの塩基置換について計算された平均およびＣＶの値を利用して適合させた、正規分布およびベータ分布に関する。（Ａ）ＥＧＦＲ＿Ｄ００１６＿００６＿Ｒ ５５２４８９５９ Ｃ＞Ｇ。（ｂ）ＥＧＦＲ＿Ｄ００１６＿００９＿Ｆ ５５２４９１５９ Ｇ＞Ｔ。（Ｃ）ＥＧＦＲ＿Ｄ００１６＿００３＿Ｒ ５５２４１７５７ Ｃ＞Ａ。（Ｄ）ＥＧＦＲ＿Ｄ００１６＿０１２＿Ｒ ５５２５９５７３ Ｇ＞Ｃ。 ＣＤＦ＿ｔｈｒｅｓｈの異なる値に関する最小検出可能頻度（ＭＤＦ）の累積値を示すグラフ。（Ｂ）は、（Ａ）の拡大図。垂直線は、９７．５パーセンタイルのポイントを示し、ここでの可能な塩基置換の９７．５％は、より低いＭＤＦ値を有する。ＣＤＦ＿ｔｈｒｅｓｈ＝０．９９９９の場合、９７．５パーセンタイルは、水平線によって示されたアレル頻度（ＡＦ）＝０．０３８２にある。 表２、表３に、細胞株の希釈系列におけるこの方法の主たる適用例を立証した結果を示す。 複数の反応での複合突然変異検出の確率を示すグラフであり、予測された分子数の関数として（材料の濃度および出発量に基づく）試料中で少なくともＮ個の分子を見出すポアソン確率を示しており、異なる線は複合コールがＮ以上の陽性反応でなされるように、異なるＮを示す。 例示的なライブラリー調製方策の概要。ＴＡｍ−Ｓｅｑプライマーは、複数の異なるバーコードを含んで作製される。各バーコードの組み合わせは、別のＰＣＲウェルに分配される。Ａ）最初のＰＣＲで目的の領域を増幅して、その反応物を、同じターゲット特異性プライマーであるが異なるバーコードを有する他の反復増幅反応物から同定する、分子識別子標識（バーコード）を付加する。Ｂ）１つの試料からのＰＣＲ産物の全てをプールして、二回目のＰＣＲでシーケンサー特異性アダプターおよび試料特異性バーコードを結合させる。 ターゲット特異性プライマーバーコードの組み合わせの代表例。各ターゲット特異性プライマーは、７つの異なるバーコードを含んで７回合成された。その後、７つのフォワードプライマーのそれぞれを７つのリバースプライマーのそれぞれと共に混合することにより、フォワードおよびリバースプライマーを組み合わせて、４９の異なるウェルのバーコード組を作製した。 コールの確率と反応物あたりの平均分子数との関連を示すグラフ。異なる値のＭＤＦ（図の凡例に示される）について、ＡＦ＞ＭＤＦの確率が、テンプレート希釈に基づく反応物あたりの平均分子数の関数としてプロットされている。水平線は、確率＝０．９０を示す。ＭＤＦ＝０．０３８２の場合、ＡＦ＞ＭＤＦの確率０．９０は、外挿されたＮ＝２０．５９（垂直線によって示される）で得られる。 ３種の細胞株からのＤＮＡの希釈系列でこの方法を実施することによる、同定された突然変異（真陽性（ＴＰ））、偽陽性コール（ＦＰ）、および見落とされた突然変異（偽陰性（ＦＮ））を示すプロット。（Ａ）Ｎが２以上のベータ分布でＣＤＦカットオフ＝０．９９９９。（Ｂ）Ｎが３以上のベータ分布でＣＤＦカットオフ＝０．９９９９。 ＡＦ予測値とＡＦ測定値の一致を示すプロット。偽陰性は、垂直軸の「検出されず（ＮＤ）」として示されている。予測値は、バツで示され、見やすいように点線で結ばれている。

発明の詳細な説明
本発明の方法は、ＤＮＡ試料を反復測定物に分配し、その中で所与の多様体が存在するならば全ＤＮＡ試料中よりも高い頻度であり、複数のウェル／反応物中の突然変異の検出と組み合わされ、それにより配列決定に用いられる方法、即ちシーケンシングプラットフォームに本質的に備わる誤り率よりも多く多様体を「コール」させることを利用する。

シーケンサーおよびポリメラーゼの誤り率は、当業者に周知の様々な方法で決定され得る（Ｔｉｎｄａｌｌ， Ｋ．Ｒ．， Ｋｕｎｋｅｌ， Ｔ．Ａ． Ｂｉｏｃｈｅｍｉｓｔｒｙ １９８８ ２７（１６）： ６００８−１３； Ｆｏｒｓｈｅｗ， Ｔ． ｅｔ ａｌ． Ｓｃｉ． Ｔｒａｎｓｌ． Ｍｅｄ． ２０１２ ４（１３６）：１３６ｒａ６８）。これらの多くは発表されており、そして／または製造業者から入手可能である。その後、誤り率を利用して、遺伝的多様体がバックグランドノイズにより観察される予測頻度を決定できる。

方法の誤差は、ＮＧＳ法の場合には、例えばライブラリーの調製の際にも導入され得る。

誤差は多くの場合、シーケンシングされる分子中の塩基の位置に応じて異なる（例えば、Ｌｏｍａｎ ＮＪ ｅｔ ａｌ．， Ｎａｔｕｒｅ Ｂｉｏｔｅｃｈｎｏｌｏｇｙ ３０： ４３４−４３９， ２０１２； Ｆｏｒｓｈｅｗ， Ｔ． ｅｔ ａｌ． Ｓｃｉ． Ｔｒａｎｓｌ． Ｍｅｄ． ２０１２ ４（１３６）：１３６ｒａ６８参照）。これに取り組むための一方法は、配列状況などのパラメータを利用してノイズをモデル化することによるが（Ｒｏｓｓ， ＭＧ ｅｔ ａｌ．， Ｇｅｎｏｍｅ Ｂｉｏｌｏｇｙ ２０１３ １４：Ｒ５１）、これは、各遺伝子座での各可能な変化についての正確な値を与えることができず、シーケンシングシステムまたは増幅工程の特性の経時的変化を説明できない。

特定の塩基での遺伝的多様体のバックグランド頻度は、参照ＤＮＡ試料の目的の領域を複数回シーケンシングすることによって、経験的に決定することもできる。例えば、参照ＤＮＡ試料は、健常対象の全血もしくは血漿から、または細胞株から得ることができる。参照試料の目的の領域を複数回シーケンシングした結果を利用して、その後、目的の領域のＤＮＡ配列の各位置での４つの可能な塩基（Ａ、Ｇ、ＴおよびＣ）のそれぞれについてバックグランド頻度を確定できる。これにより、方法の誤差（即ち、用いられるＤＮＡシーケンシングプラットフォームの誤り率）によって遺伝的多様体が同定される頻度を、目的の領域のＤＮＡ配列の各位置について決定できる。幾つかの実施形態において、参照ＤＮＡ試料は、「マッチした正常な」試料であり得る。例えば、参照ＤＮＡ試料は、健常対象からの同じ組織および／または試料のタイプから得ることができる。

本発明の方法の幾つかの実施形態において、遺伝的多様体の頻度の平均および分散は、目的の領域を含む特異的アンプリコンパネルについて決定される。

本明細書で用いられる「目的の領域」は、検討されるゲノムの一部分または複数の部分を意味する。それは、ＤＮＡの１つの配列またはＤＮＡの複数の配列であり得る。目的の領域が複数の領域である場合、これらはゲノム全体に広がり得る。即ち、「目的の領域」は、検討されるＤＮＡの複数の配列に、検討されないゲノムの部分が散在したものを包含する。幾つかの例において、目的の領域は、例えば検討される任意の臨床的問題に依存する。

例えば、幾つかの実施形態において、目的の領域は、疾患関連の遺伝的多様体または治療への応答に影響を及ぼす遺伝的多様体を含むことが知られている、または含む候補である、複数の領域であり得る。例えば、幾つかの実施形態において、目的の領域は、一連の癌関連遺伝子であり得る。

所与の位置にある所与のヌクレオチドの頻度（アレル頻度、ＡＦ）は、その位置のそのヌクレオチドの総決定数からその位置のそのヌクレオチドを同定する、その位置のリードの割合として決定される。

参照試料の各シーケンシング反応から決定された目的の領域の各位置にある各ヌクレオチドのＡＦをその後利用して、平均ＡＦ、およびＡＦの分散係数（ＣＶ）を計算し得る。ＡＦのＣＶは、平均ＡＦで割られたＡＦの標準偏差として計算される。

平均ＡＦおよびＡＦのＣＶをその後用いて、所与の遺伝的多様体についてバックグランドノイズ（即ち、シーケンシング誤差）をモデル化できる。例えば、正規分布、ベータ分布、指数分布もしくはガンマ分布、または他の関数を用い、どのモデルがデータの経験的分布と最良適合するかに応じて、遺伝的多様体のバックグランドノイズをモデル化できる。同様に、突然変異体のリード数およびシーケンシング深度を、別個の関数でモデル化できる。

確率分布およびそれらをデータにモデル化する方法は、当業者に周知である。好ましくは複数の分布をＡＦデータにモデル化して、最良適合したモデルを、本発明の方法の次の段階でのモデルとして選択する。例えば、異なるモデルを、経験的データへの適合度のコルモゴロフ・スミルノフ値について解析できる。

以下の実施例では、ベータおよび正規確率分布が用いられている。

遺伝的多様体の存在を決定するための閾値
経験的に決定されたバックグランド誤差ＡＦおよびＣＶ値に適合させた確率分布、またはシーケンサーおよびポリメラーゼの誤り率、目的の領域ならびに遺伝的多様体に基づく予測されたバックグランド誤差に基づく確率分布をその後用いて、それ以上で反復増幅反応物のシーケンシング結果において遺伝的多様体が観察されてＤＮＡ試料に存在すると決定されるべき、閾値ＡＦ（本明細書では最小検出可能頻度（ＭＤＦ）と称される）を確定し得る。

例えば、以下の実施例７において、２０および３０のｚスコアでのＡＦが選択される。即ち、所与の反復増幅反応物中のＤＮＡ試料に存在すると決定される遺伝的多様体について、そのアレルの頻度は、その特定の多様体に関するバックグランド平均ＡＦを超える２０以上、または３０以上の標準偏差でなければならなかった。

問題となる特定の遺伝的多様体、目的の領域、用いられるシーケンシングプラットフォームなどに応じて、様々なｚスコアカットオフが、本発明の方法で有用である。例えばｚスコア閾値は、１、２、３、４、５、６、７、８、９、１０、１１、１２、１３、１４、１５、１６、１７、１９、２０、２２、２４、２６、２８、３０、３５、４０、４５、５０、６０、７０、８０、９０および１００のうちの１つから選択され得る。

確率分布を利用して、多様体が実際に反復増幅反応物中に存在する場合に、特定の遺伝的多様体を所与のＡＦで観察する累積確率を決定できる。これを次に利用して、反復増幅反応物の産物においてそれ以上で遺伝的多様体が観察されて、ＤＮＡ試料に存在すると決定されるべき閾値頻度を確定する。

シーケンシングデータにおいてＭＤＦ以上の頻度で観察された推定突然変異がシーケンシング誤差の結果よりもむしろ「現実的」である確率は、非常に高い。そのような多様体はそれゆえ、高い信頼度で所与の増幅反応物中に存在すると決定され得る。

本発明の方法の幾つかの実施形態において、ステップ（ｉｉ）の閾値頻度は、遺伝的多様体のバックグランド頻度のベータ確率分布モデルを利用して決定される。この場合、目的の領域の各位置にある各可能な遺伝的多様体について、ステップ（ｉ）で決定された遺伝的多様体のバックグランド平均頻度および頻度の分散を用いて、ベータ分布を定義する。所与の遺伝的多様体が観察されて反復増幅反応物中に存在すると決定されるべき閾値頻度は、その遺伝的多様体の累積分布関数（ＣＤＦ）値が所定の閾値（ＣＤＦ＿ｔｈｒｅｓｈ）に達する頻度である。幾つかの実施形態において、ＣＤＦ＿ｔｈｒｅｓｈは、０．９９、０．９９５、０．９９９、０．９９９９、０．９９９９９またはそれを超える。

例えば、以下の実施例６において、特定の遺伝的多様体を所与のＡＦで観察する累積確率（ベータ分布に基づく）０．９９９９を、陽性判定についての閾値ＡＦ（即ち、ＭＤＦ）として用いた。即ち、遺伝的多様体が存在しない場合に、反復増幅反応物でのシーケンシング結果において所与の遺伝的多様体を観察する確率は、０．０１％以下である。幾つかの実施形態において、閾値頻度は、遺伝的多様体をその頻度（ＣＤＦ＿ｔｈｒｅｓｈ）で観察する確率に関する累積分布関数が０．８、０．８５、０．９、０．９５、０．９６、０．９７、０．９８、０．９９、０．９９５、０．９９９、０．９９９５、０．９９９９、０．９９９９５、０．９９９９９、０．９９９９９５、０．９９９９９９またはそれを超える頻度であり、それは遺伝的多様体が存在しない増幅反応物のシーケンシング結果において遺伝的多様体が観察される確率２０％、１５％、１０％、５％、４％、３％、２％、１％、０．５％、０．１％、０．０５％、０．０１％、０．００５％、０．００１％、０．０００５％、０．０００１％またはそれ未満に対応する。

幾つかの実施形態において、ＭＤＦ値が複数の可能な遺伝的多様体で決定され、続いてそれを用いて目的の所与の領域および／または目的の一連の遺伝的多様体について全体的ＭＤＦを確定する。したがって幾つかの実施形態において、それ以上で増幅反応物のシーケンシング結果において遺伝的多様体が観察されるべき閾値頻度を確定して所与の増幅反応での遺伝的多様体の存在について陽性判定を割り付けることは、
（ａ）ステップ（ｉ）で複数の遺伝的多様体について決定された平均頻度および頻度の分散に基づいて、それ以上で増幅反応物のシーケンシング結果において遺伝的多様体が観察されるべき複数の閾値頻度を確定して所与の増幅反応における遺伝的多様体の存在について陽性判定を割り付けること（即ち、ＭＤＦを複数の遺伝的多様体について確定すること）、および
（ｂ）ステップ（ａ）に基づいて、それ以上で増幅反応物のシーケンシング結果において遺伝的多様体が観察されるべき全体的閾値頻度を確定してその増幅反応における遺伝的多様体の存在について陽性判定を割り付けること（即ち、全体的ＭＤＦを確定すること）、
を含み、それは、ステップ（ａ）で決定された複数の閾値頻度の少なくとも９０％、９５％、９７．５％、または９９％がこの値未満である閾値頻度である。

幾つかの実施形態において、全体的ＭＤＦは、そこで複数の遺伝的多様体について決定されたＭＤＦ値の８０％、８５％、９０％、９１％、９２％、９３％、９４％、９５％、９６％、９７％、９７．５％、９８％、９８．５％、９９％もしくは９９．５％またはそれを超える値が全体的ＭＤＦよりも小さい（即ち、全体的ＭＤＦ未満）頻度である。幾つかの実施形態において、全体的ＭＤＦは、目的の領域全体に分布される複数の遺伝的多様体について確定されたＭＤＦ値に基づいて決定される。幾つかの実施形態において、複数の遺伝的多様体は、５、１０、２０、２５、３０またはそれを超える多様体であろう。

遺伝的多様体の存在を決定するための閾値は、（所与の反復増幅反応物における目的の領域のＤＮＡテンプレートのコピー数と共に）、存在するならば所与の反復増幅反応物における遺伝的多様体の頻度がその遺伝的多様性のバックグランドレベルよりも有意に高くなるように選択される。こうして、所与の増幅反応物における遺伝的多様体の有無は、偽陽性判定の確率を最小限にして、信頼性を持って決定できる。

同様に、遺伝的多様体の存在を決定するための閾値は、偽陰性の判定数を最小限にするように選択される。

反応物あたりの多数の増幅可能なテンプレート分子
突然変異が同定され得る下限は、遺伝的多様体のバックグランド頻度（即ち、方法誤差により観察される頻度）により決定される。即ち、例えば誤った決定が５％の頻度でなされる方法では、それらの頻度が５％を実質的に超える場合に、実際の多様体が、誤りと区別され得るに過ぎない。

本発明の方法は、所与の反復増幅反応物における目的の領域の増幅可能なテンプレート分子の平均数に上限を設定することによって、この問題を克服している。

本発明の方法において、所与の反復増幅反応物における目的の領域の増幅可能なテンプレート分子の平均数（Ｍｏｌ＿ｍｅａｎ）は、以下の通り、遺伝的多様体の存在を決定するための閾値（ＭＤＦ）を用いて決定される：
Ｍｏｌ＿ｍｅａｎ ＜ １／ＭＤＦ。

したがって遺伝的多様体は、所与の反復増幅反応物中に存在する場合、所与の増幅反応において遺伝的多様体の存在を決定するための閾値頻度よりも大きな頻度で、その反復増幅反応物内で目的の領域の増幅可能なテンプレートコピーの平均総数以内で表される（分率として）と予測されよう。

こうして、多様体を有する単一分子は、ステップ（ｉｉ）において決定されるＭＤＦ以上の頻度で検出される確率が高い。

この場合、ＭＤＦは、所与の遺伝的多様体について計算されたＭＤＦ、または複数のＭＤＦを代表する包括的値であり本明細書の先に記載されている全体的ＭＤＦであり得る。

例えば、所与の遺伝的多様体が、例えば５％以上の頻度で観察された場合に（ＭＤＦによって定義された通り）ＤＮＡ試料中に存在していると信頼性をもって決定され得るのにすぎないならば、その方法は、平均で２０未満の分子が所与の反復増幅反応物中に存在することを必要とする。こうして、遺伝的多様体が存在する（そして最大１９の非多様体分子が存在する）任意の所与の反復測定物では、多様体は、５％以上の頻度で観察されるであろう。

以下に提供された実施例において、反復増幅反応物への分子の分布が、ポアソン分布でモデル化されている。所与の増幅反応において遺伝的多様体を有する増幅可能な分子についての陽性判定の確率は、ポアソン分布に基づいて決定され得る。あるいは反復増幅反応物への分子の分布は、負の二項分布でモデル化され得る。

したがって、幾つかの実施形態において、Ｍｏｌ＿ｍｅａｎは、ポアソン分布を用いて決定され、所与の反復増幅反応物中に存在する増幅可能な分子の数および増幅可能な分子の予測される平均数をモデル化する。Ｍｏｌ＿ｍｅａｎの異なる可能な値について、ポアソン分布の累積分布関数が計算される。ポアソン分布での累積分布関数値が０．８、０．８５、０．９、０．９５または０．９９などの閾値よりも大きくなるように、Ｍｏｌ＿ｍｅａｎがその後選択される。

本発明の方法の幾つかの実施形態において、反応物あたりの増幅可能なテンプレート分子の平均数が反応物における分子の予測されるポアソン分布に基づくように、ＤＮＡ試料が分配され、遺伝的多様体がもし存在するならば、９０％を超える確率で検出される。この閾値は、より高くてもより低くてもよく、反応物あたりの平均分子数が１に減少すると、定常に達する（例えば図８参照）。幾つかの実施形態において、遺伝的多様体が存在するならば、８０％を超える、８１％を超える、８２％を超える、８３％を超える、８４％を超える、８５％を超える、８６％を超える、８７％を超える、８８％を超える、８９％を超える、９０％を超える、９１％を超える、９２％を超える、９３％を超える、９４％を超える、９５％を超える、９６％を超える、９７％を超える、９８％を超えるおよび９９％を超える確率のうちの１つで検出されるように、反復増幅反応物あたりの増幅可能なテンプレート分子の平均数が選択される。

時間および資源の効率および費用効果的使用という関心事については、ＤＮＡ試料は、テンプレート分子の分配に関連する方法の効率を実現するのに必要な量よりも多く分配されないことが好ましい。このことは、少なくともランニングコストおよび材料費を低減することに関して利点を有する。もちろん増幅反応物あたりの最小分子数は、決定される遺伝的多様体、目的の領域、方法の誤差（即ち、ＤＮＡシーケンシングプラットフォームの誤り率）および計算された閾値などに依存するであろう。本発明の方法の特定の実施形態において、目的の領域において１つよりも多くの増幅可能なテンプレート分子が反復増幅反応物あたりに提供されることが、好ましい。

例えば、幾つかの実施形態において、反復増幅反応物あたりの増幅可能なテンプレート分子の平均数は、１よりも多くＭｏｌ＿ｍｅａｎ未満であろう。

幾つかの実施形態において、反復増幅反応物あたりの増幅可能なテンプレート分子の平均数は、１よりも多く１０００未満、１よりも多く７５０未満、１よりも多く５００未満、１よりも多く２５０未満、１よりも多く２００未満、１よりも多く１５０未満、１よりも多く１００未満、１よりも多く８０未満、１よりも多く６０未満、１よりも多く５０未満、１よりも多く４０未満、１よりも多く３５未満、１よりも多く３０未満、１よりも多く２９未満、１よりも多く２８未満、１よりも多く２７未満、１よりも多く２６未満、１よりも多く２５未満、１よりも多く２４未満、１よりも多く２３未満、１よりも多く２２未満、１よりも多く２１未満、１よりも多く２０未満、１よりも多く１９未満、１よりも多く１８未満、１よりも多く１７未満、１よりも多く１６未満、１よりも多く１５未満、１よりも多く１４未満、１よりも多く１３未満、１よりも多く１２未満、１よりも多く１１未満、１よりも多く１０未満、１よりも多く９未満、１よりも多く８未満、１よりも多く７未満、１よりも多く６未満、１よりも多く５未満、１よりも多く４未満、１よりも多く３未満、または１よりも多く２未満である。

幾つかの実施形態において、反復増幅反応物あたりの増幅可能なテンプレート分子の平均数は、２よりも多く１０００未満、２よりも多く１００未満、２よりも多く７５未満、２よりも多く５０未満、２よりも多く４０未満、２よりも多く３０未満、２よりも多く２９未満、２よりも多く２８未満、２よりも多く２７未満、２よりも多く２６未満、２よりも多く２５未満、２よりも多く２４未満、２よりも多く２３未満、２よりも多く２２未満、２よりも多く２１未満、２よりも多く２０未満、２よりも多く１９未満、２よりも多く１８未満、２よりも多く１７未満、２よりも多く１６未満、２よりも多く１５未満、２よりも多く１４未満、２よりも多く１３未満、２よりも多く１２未満、２よりも多く１１未満、２よりも多く１０未満、２よりも多く９未満、２よりも多く８未満、２よりも多く７未満、２よりも多く６未満、２よりも多く５未満、２よりも多く４未満、または２よりも多く３未満である。

幾つかの実施形態において、反復増幅反応物あたりの増幅可能なテンプレート分子の平均数は、５よりも多く１０００未満、５よりも多く１００未満、５よりも多く７５未満、５よりも多く５０未満、５よりも多く４０未満、５よりも多く３０未満、５よりも多く２９未満、５よりも多く２８未満、５よりも多く２７未満、５よりも多く２６未満、５よりも多く２５未満、５よりも多く２４未満、５よりも多く２３未満、５よりも多く２２未満、５よりも多く２１未満、５よりも多く２０未満、５よりも多く１９未満、５よりも多く１８未満、５よりも多く１７未満、５よりも多く１６未満、５よりも多く１５未満、５よりも多く１４未満、５よりも多く１３未満、５よりも多く１２未満、５よりも多く１１未満、５よりも多く１０未満、５よりも多く９未満、５よりも多く８未満、５よりも多く７未満、または５よりも多く６未満である。

幾つかの実施形態において、反復増幅反応物あたりの増幅可能なテンプレート分子の平均数は、２〜４０、３〜３０、４〜３０、５〜３０、５〜２５、１０〜２５、１５〜２５、または１８〜２２の範囲である。

「分配」は、反復増幅反応物を分離するのに適した任意の手段により実現され得る。例えば、分配は、（希釈された）ＤＮＡ試料を別々のウェル中に分取することにより実現され得る。別々の（即ち、分離した、個々のまたは独立した）反復増幅反応物に区分するための様々な他の方法が、当業者に公知であろう。

陽性判定のために複数の陽性反復測定を必要とすること
ＤＮＡ試料中の遺伝的多様体の存在の偽陽性判定の確率を最小限にするために、本発明の方法は、遺伝的多様体が１つよりも多くの反復測定物中に存在すると決定されることを必要とする。

反応物あたりの平均分子数、および偽陽性のコールを最小限にするのに必要とされる反復測定コール数が、必要とされる感度を得るために実行される反応物数を決定する。

全ての反応物数ＴのうちＮの反応物中で起こる偽陽性の理論的確率（Ｐ＿ｆｐＮＴ）は、１つの反応物中の偽陽性の確率（Ｐ＿ｆｐ１Ｔ、１−ＣＤＦ＿ｔｈｒｅｓｈに等しい）および反応物の数に依存し、
Ｐ＿ｆｐＮＴ＝（（Ｐ＿ｆｐ１Ｔ）＾Ｎ）^＊Ｃ＿ＮＴ＝（（１−ＣＤＦ＿ｔｈｒｅｓｈ）＾Ｎ）^＊Ｃ＿ＮＴ
に等しく、ここで、Ｃ＿ＮＴは、以下の通り、反応物の総数Ｔに依存する組み合わせの係数である（「！」は、階乗関数を示す）
Ｃ＿ＮＴ＝Ｔ！／（（Ｔ−Ｎ）！^＊Ｎ！）

この工程から予測される偽陽性の総数は、この工程で検査される多様体数を掛けた、偽陽性の確率（Ｐ＿ｆｐＮＴ）に等しい。例えば、一ヌクレオチド多様体が検査される場合、これは、３をかけた位置の数になろう（３つの可能な非参照変化の場合）。

複数の反復増幅反応物中の遺伝的多様体で陽性判定を必要とすることで、ＤＮＡ試料中に存在する遺伝的多様体の偽陽性判定の確率が低下する。これは、以下の実験的実施例から明確である。

反復増幅反応物中で複数の陽性判定を必要とする方法はしたがって、遺伝的多様体の検出のより高い特異度を有する。

しかし、所与の遺伝的多様体について陽性である複数の反復測定物を必要とすることはまた、偽陰性判定がなされる確率を増加させる。図５は、少なくともＮ個のＤＮＡテンプレート分子が予測された分子数の関数として（材料の濃度および出発量に基づき）ＤＮＡ試料中に存在する、ポアソン確率を示す。

突然変異体分子の無作為な分布を前提とすれば、少なくとも必要とされる数の反復測定物が１つの突然変異を含む確率が十分に高くなるように、反応物の数、およびこれらの反応物中の多様体分子の数は、十分に高くなければならない。

それゆえ、ステップ（ｖｉ）において遺伝的多様体がＤＮＡ試料中に存在すると決定されるのに必要とされる陽性反復測定物の数は、偽陽性および偽陰性を最小限にするように選択される。即ち、必要とされる陽性反復測定物の数は、この方法の所望の感度および特異度を実現するように決定されるであろう。幾つかの実施形態において、陽性反復測定物の数は、この方法の感度および特異度を最大にする数であろう。

その数は常に、１を超えるであろう。幾つかの実施形態において、その数は、２、３、４、５、６、７、８、９、または１０の反復測定物であろう。特定の実施形態において、その数は、２または３であろう。

その数は、分析される領域の数およびサイズ、利用可能なＤＮＡの量、反復増幅反応物の総数などに応じて変動し得る。

その数はまた、参照ヌクレオチドからの多様性の特定のタイプに基づいて、例えば多様体が参照アレルからの転位であるか、または転換であるかに基づいて変動し得る

必要とされる陽性反復測定物の数は、ＤＮＡ試料についての反復増幅反応物の総数によって一部は決定され得る。例えば幾つかの実施形態において、ステップ（ｖｉ）において遺伝的多様体がＤＮＡ試料に存在すると決定される場合、多様体の存在についての陽性判定が、反復測定物の１％、２％、３％、４％、５％、６％、７％、８％、９％または１０％よりも多くなされなければならない。

所与の多様体および／または目的の領域について必要とされる陽性反復測定物の最適数は（偽陽性および偽陰性を最小限にするため）、例えば参照ＤＮＡ試料で本発明の方法を実施することによって、決定され得る。

本発明の方法は、ＤＮＡ試料中の目的の領域における遺伝的多様体の有無の最終決定のために、反復増幅反応物についての有無決定の結果の統合を必要とする。

本明細書で用いられる「統合する」は、反復増幅反応物についての有無の結果をまとめること、または合併整理することを意味する。

多様体コーリングの閾値と、増幅可能なテンプレートＤＮＡ分子の数と 反復増幅反応物の数との関連
先に、そして以下の実験的実施例に記載された通り、反復増幅反応物中の遺伝的多様体の存在を決定するための閾値ＡＦと、所与の反復増幅反応物中に存在する増幅可能なテンプレートＤＮＡ分子の数と、遺伝的多様体がＤＮＡ試料中に存在すると決定されるために「陽性」であることが必要とされる反復測定物の数と、が相互に関連することは、当業者に直ちに明白であろう。

本質的に、主要な原理は、所与の反復増幅反応物中に遺伝的多様体が存在する場合、その頻度がその遺伝的多様体のバックグランドレベルよりも有意に高くなるように、目的の領域の遺伝的多様体の存在およびＤＮＡテンプレートのコピー数を決定するための閾値が選択されなければならないこと、ならびに陽性である反復測定物の数がこの方法の感度および特異度を最適にするように選択されなければならないこと、である。

こうして、所与の増幅反応物中の遺伝的多様体の有無は、偽陽性判定の確率を最小限に抑えて、信頼性をもって決定できる。

これらのパラメータを確定する上での感度および特異度の相対的重要性は、個々の研究にも依存する。例えば、本発明の方法が新規な突然変異を同定するのに用いられている場合、主要な判断事項は、感度であり得る。逆に、ＤＮＡ試料が治療決定の情報を与えるために多様体の存在について分析されている場合、主要な判断事項は、特異度であり得る。

標識されたアンプリコンのディープシーケンシング
本発明の方法は、ハイスループットＤＮＡシーケンシング法での使用に適する。有利には、目的の複数の領域が、突然変異について並行してスクリーニングされ得る。

以下の実験的実施例において、標識されたアンプリコンのディープシーケンシング（ＴＡｍ−Ｓｅｑ）法が用いられる。この方法は、Ｆｏｒｓｈｅｗ ｅｔ ａｌ． ２０１２ Ｓｃｉ Ｔｒａｎｓｌ Ｍｅｄ ４（１３６） １３６ｒａ６８に詳細に記載される。ＴＡｍ−Ｓｅｑは、ＤＮＡ断片の個々のコピーという小さなものから数千の塩基にわたるゲノム配列の増幅およびディープシーケンシングを可能にする。

手短に述べると、プライマーは、テンプレートＤＮＡの性質（例えば、平均断片長）に基づいて選択されたサイズ範囲の断片で目的の領域を並べて表示するアンプリコンを作製するよう設計され、汎用シーケンシングアダプター配列を組み込み各反復測定物を識別子配列または「バーコード」で標識する（図６および７参照）。その後、産物をシーケンシングして、識別子配列を用いてリードをデマルチプレックスして、それらをゲノムに割り付ける。

本発明の方法の使用
この方法は、広範囲の適用例に有用であり、その適用例は、熟練者に直ちに明白となろう。本質的にこの方法は、目的の任意の試料における任意の遺伝的多様体の検出に有用である。

例えば本発明の方法は、診断および／もしくは予知方法、または微生物もしくはウイルス分析において有用である。

特に、この方法は、希少な遺伝的多様体および／または希少な突然変異の検出に有用である。この方法は、ＤＮＡ試料において１％未満の頻度で存在する多様体の検出を可能にする。さらに方法は、そのような低頻度での新規多様体を検出すること、および試料を既知の多様体について分析すること、に適している。

したがって、この方法は、癌であるまたは妊娠時などの体内に出現し得る遺伝的多様体の存在のより早期の同定を可能にする（即ち、遺伝的多様体が低頻度で存在する場合）。

本発明の任意の態様の方法による幾つかの例において、この方法は、対象から過去に得られた試料をアッセイすることを含む。試料は一般に、ＤＮＡが単離され得る任意の適切な生物学的試料であり得る。幾つかの例において、試料は、尿、唾液、血液、血清、糞便、他の生物学的流体、毛髪、細胞および組織からなる群から選択される。

本発明の任意の態様の方法による幾つかの例において、ＤＮＡは、この方法を実施する前にビスルファイトで処置され得る。ＤＮＡのビスルファイト処置は、非メチル化シトシン残基をウラシルに変換するが、メチル化シトシン残基は影響を受けないままである。こうしてこの方法は、ＤＮＡメチル化において多様体を検出するのに有用である。

低頻度の多様体の存在を検出する方法の能力は、非侵襲的な、そして／または最小限に侵襲的な手段によって得られる試料の分析を可能にする。例えば血液、血清、尿、精液、糞便、痰、胸水、腹水、滑液、脳脊髄液または気管支洗浄液の試料などの液体生検。これは、対象への大きなリスクまたは不快感を伴わずに試料を定期的に採取し得るという利点を有する。試料は、新鮮であっても、過去に貯蔵されていても（例えば、凍結）、そして／または過去に加工されていてもよい。

その上、この方法は、定量的であり、例えば臨床的相関のための循環する腫瘍ＤＮＡ（ｃｔＤＮＡ）の突然変異レベルの正確な測定、またはウイルスもしくは微生物ＤＮＡなどの他のＤＮＡの正確な測定を可能にする。

定量は、例えば多様体を有すると決定された分子の数をカウントすることによって、実施され得る。これは、多様体を有する分子の数が少ない場合に適する。定量はまた、カウントのポアソン補正によって、または観察されたアレル頻度をモデル化することによって実施され得る。定量はまた、複数の反応物それぞれの中の遺伝的多様体のアレル頻度、および各反応物中の相対的テンプレート量を考慮することによって、例えば平均もしくは材料の量を加重された平均を利用することによって、または異なる量の出発原料に関して多様体が検出される反応物の比率を考慮することによって、実施され得る。

例えばこの方法は、腫瘍量モニタリングに有用であり、その場合、この方法は、ｃｔＤＮＡ突然変異のレベルをモニタリングするため、または突然変異の再出現をモニタリングするため、腫瘍の再成長の早期指示物質として情報を与えるために用いられる。この方法は、明らかに診断的および予知的用途を有する。

この方法はまた、治療への応答において、薬物耐性および／または腫瘍進展をモニタリングするのに有用である。この方法は、多数の目的の領域を並行して評価するのに特によく適しており、薬物耐性の早期指示物質であり得る可能な耐性経路での突然変異を検出するのに用いることができる。そのような用途は、過去の方法の能力を上回る可能性がある。

さらにこの方法は、腫瘍の一次診断において有用である。この方法を利用して、早期診断のために健常患者でｃｔＤＮＡ突然変異を探索できる。例えば、重要な発癌遺伝子および腫瘍抑制遺伝子を、モニタリングできる。

この方法はまた、遺伝的多様体の日常的スクリーニングまたはテスト（即ち、モニタリング）のために、健常対象から得られた試料でも有用である。幾つかの実施形態において、健常対象は疾患から回復していてもよく、例えば対象は癌が鎮静中であってもよい。あるいは健常対象は疾患のリスクが増大していてもよく、例えば健常対象は疾患の家族歴を有していてもよい。本発明の方法を利用して、遺伝的多様体について、場合により所与の疾患に関連することが知られる遺伝的多様体について、日常的にスクリーニングし得る。こうして、健常対象から得られた試料を用いて、臨床症状発現の前に疾患関連の遺伝的多様体を検出し、それにより早期治療介入を容易にすることができる。

本発明の方法は、治療決定の情報を与えるのに有用である。即ち、易罹患性、治療への耐性または応答を予測する遺伝的多様体の同定を利用して、対象に関する適切な処置過程を選択できる。

例として、ＥＧＦＲにおけるＴ７９０Ｍ突然変異は、ＥＦＧＲ阻害剤ゲフィチニブおよびエルロチニブを用いる癌処置への耐性に関連することが公知である。それゆえ、この突然変異が本発明の方法を用いて同定される対象は、これらの阻害剤を用いる処置の不適当な候補であると同定され得る。

同様に、Ｂ−ＲａｆにおけるＶ６００Ｅ突然変異は、Ｂ−Ｒａｆ阻害剤への感受性増大に関連する。それゆえ、この突然変異が本発明の方法を用いて同定される対象は、これらの阻害剤を用いる処置の良好な候補であると同定され得る。

この方法は、異なる時点に、例えば異なる疾患段階または処置過程に、対象から得られたＤＮＡ試料で実施できる。

この方法を利用して、易罹患性、進行（即ち、予後）、所与の治療への耐性または応答に関連する、またはそれを予測する遺伝的多様体について、ＤＮＡ試料を分析できる。例えば、この方法を用いて、疾患の易罹患性もしくは予後、または治療への耐性もしくは応答の新規な（即ち、過去に同定されていない）決定子および／または予測子を同定できる。

本発明の方法は、定量的であるため、疾患の異なる段階で処置への応答において、そして／または経時的に健常対象で、遺伝的多様体の相対頻度を検査するのに有用である。

この情報は今度は、治療決定についての情報を与えるために用いることができる。例えば、易罹患性または特定の治療への耐性に関連することが既知である、そして／または予測される遺伝的多様体の頻度における上昇または低下は、どの治療が所与の対象の処置に最も適するかについての決定を誘導するであろう。

実施例
実施例１ − バックグランドノイズの決定
ＴＰ５３ならびにＥＧＦＲ、ＫＲＡＳ、ＢＲＡＦおよびＰＩＫ３ＣＡのホットスポット領域の全てに及ぶように４１のアンプリコンを設計して、前方および後方に読んで、合計８２のリードファミリーを与えた（表１）。アンプリコンは、前方および後方のリードの重複を含む（プライマー配列を除く）５，０３８の塩基に及ぶ。

各遺伝子座は、３つの可能な非参照アレルのうちの１つに変化し得る。既知の多形、用いられた細胞株中で同定された多型、およびモデル化され得ない遺伝子座を除く、合計１３，２７８の可能な一塩基置換が考慮された。このアンプリコンのパネルについてのシーケンシングデータを、ＬＮＣａＰ細胞株（ヒト前立腺腺癌、アンドロゲン感受性細胞株；ＡＴＣＣ ＣＲＬ−１７４０）ＤＮＡ試料の３３６の反復増幅物から回収した。

先に除外された遺伝子座以外の各リードファミリーにおける各塩基置換（アンプリコンおよびリードの方向）について、平均アレル頻度（ＡＦ）およびＡＦの変動係数（ＣＶ）を計算した。所与のアレル（即ち、遺伝的多様体）について、ＡＦを、そのリードファミリーの全リードのうちのこのアレルを含むリードの割合として計算した。

図１Ａは、増加するＡＦ値によってビン化された平均ＡＦの異なる値の出現数を示す。垂直線は、平均＝０．０００２２を示し、それは分布の９５％パーセンタイルである。

図１Ｂは、増加するＣＶ値によってビン化されたＣＶの異なる値の出現数を示す。垂直線は、ＣＶ＝１８．１６を示し、それは分布の９５％パーセンタイルである。

実施例２ − 所与の反応における塩基置換について陽性コールに必要なＡＦの決定
観察されたＡＦがそのアレルの最小検出可能頻度（ＭＤＦ）よりも大きいと見出された場合、反応を所与の配列変化について陽性と定義した。アレルについてのＭＤＦを、その特定の塩基置換の累積分布関数（ＣＤＦ）が所定の閾値（ＣＤＦ＿ｔｈｒｅｓｈ）を交差するＡＦとして計算した。

図２は、実施例１のパネルの複数の塩基置換に関する累積分布を示す。経験的データ分布と比較した、それらの塩基置換について計算された平均ＡＦおよびＣＶの値を利用して適合させた、正規分布およびベータ分布の累積分布関数が示される。経験的データ分布に対する正規分布（ＫＳ Ｎｏｒｍａｌ）およびベータ分布（ＫＳ Ｂｅｔａ）についての適合度のコルモゴロフ・スミルノフ値は、以下の通りである。

記載された工程を利用して、複数の試料中の１３，２７８の可能な一塩基置換をスキャンした。偽陽性を低い割合で維持するために、ＡＦ＞ＭＤＦとなるＡＦおよびＣＤＦ＿ｔｈｒｅｓｈ＝０．９９９９で突然変異が存在する場合には、突然変異をコールした。

バックグランドノイズをベータ分布としてモデル化し、１３，２７８の可能な塩基置換それぞれで、対照試料における塩基置換で測定された平均ＡＦおよびＣＶを条件に（図１）、
Ｆ（ＡＦ｜Ｍｅａｎ，ＣＶ）＝ＣＤＦ＿ｔｈｒｅｓｈ
（式中、「Ｆ」は、ベータ分布のＣＤＦである）
でのＡＦである、ＭＤＦ（ＣＤＦ＿ｔｈｒｅｓｈ）を計算した。図３は、異なる値のＣＤＦ＿ｔｈｒｅｓｈ（凡例に示された）についてのＭＤＦ（ＣＤＦ＿ｔｈｒｅｓｈ）を示す。各例のデータは、得られたＭＤＦの値によって区分されている。

図３Ｂは、拡大図を示す。垂直線は、可能な塩基置換の９７．５％がより低いＡＦ値を有する９７．５パーセンタイルポイントを示す。ＣＤＦ＿ｔｈｒｅｓｈ＝０．９９９９の場合、９７．５パーセンタイルは、水平線で示されたＡＦ＝０．０３８２の位置である。

したがって、１３，２７８の可能な一塩基置換の場合、ＣＤＦ＿ｔｈｒｅｓｈ＝０．９９９９でのＭＤＦは、可能な一塩基置換の９７．５％について０．０３８２以下となろう（図３参照）。

各反応における所与の塩基置換についての偽陽性判定の確率は、このモデルによれば、１−ＣＤＦ＿ｔｈｒｅｓｈ（ＡＦ＝ＭＤＦの時のＣＤＦの値）である：
Ｐ＿ｆｐ１Ｔ＝１−Ｆ（ＡＦ＝ＭＤＦ｜Ｍｅａｎ，ＣＶ）＝１−ＣＤＦ＿ｔｈｒｅｓｈ

実施例３ − 試料中の塩基置換の存在を決定するための所与の塩基置換について陽性である反復増幅反応物の数の決定
試料中の塩基置換の存在についての偽陽性判定を最小限にするために、その試料の反復増幅反応物の総数（Ｔ）のうち、複数の陽性反復増幅反応物（Ｎ）が必要とされた。

Ｔの反応物のうちＮに出現する偽陽性の理論的確率（Ｐ＿ｆｐＮＴ）は、
Ｐ＿ｆｐＮＴ＝（Ｐ＿ｆｐ１Ｔ）＾Ｎ）^＊Ｃ＿ＮＴ＝（（１−ＣＤＦ＿ｔｈｒｅｓｈ）＾Ｎ）^＊Ｃ＿ＮＴ
であり、ここで「Ｃ＿ＮＴ」は、順序には関係なくＴのうちＮのウェルを選択する可能性の数であり、
Ｃ＿ＮＴ＝Ｔ！／（（Ｔ−Ｎ）！^＊Ｎ！）
である（ここで、「！」は、階乗関数を示す）。

例えば、４８のうちの３の反応物を選択する場合、１７，２９６の可能性が存在する。４８のうち２の反応を選択する場合、１，１２８の可能性が存在する。予測された偽陽性率を各試料について計算し、表２および３に示している（図４）。

複数の陽性反応を観察する必要性は、偽陰性の確率も上昇させ、真陽性が複数の反応物中で観察される確率を低下させる。

図５は、予測される分子数（材料の濃度および出発量に基づく）の関数としての試料中の少なくともＮの分子数を見出すポアソン確率を示す。複合コールがＮ以上の陽性反応物でなされるように、異なる線が異なるＮを示す。

（ｉ）各試料の反復増幅反応物の総数（Ｔ）が、Ｎよりもかなり多くなるように（この場合、Ｔ≧４８）、そして
（ｉｉ）突然変異体分子が存在した場合に、任意の反応物においてコールの確率が１に近似するように、
システムを設計した。

実施例４ − 各増幅反応物に存在する増幅可能な分子の数の決定
予定される希釈率および反応物あたりの対応する平均分子数では、実際の反応物は、ポアソン分布でモデル化され得る無作為な分子数を有すると予測される。

それぞれのそのようなプール中の単一突然変異分子によって、分子の数の逆数と等しい観察ＡＦが得られるであろう（リード分率が、突然変異体アレル頻度の代表であると仮定、Ｆｏｒｓｈｅｗ ｅｔ ａｌ． ２０１２ Ｓｃｉ Ｔｒａｎｓｌ Ｍｅｄ ４（１３６） １３６ｒａ６８， Ｆｉｇｕｒｅ ３ｂ参照）。

異なる値のＭＤＦについて（図の凡例に示される）、ＡＦ＞ＭＤＦの確率を、テンプレート希釈に基づいて、反応物あたりの平均分子数の関数としてプロットした（図８）。

水平線は、確率＝０．９０を示す。ＭＤＦ＝０．０３８２（塩基置換の９７．５％がＣＤＦ＿ｔｈｒｅｓｈ＝０．９９９９を通過する、図３参照）では、ＡＦ＞ＭＤＦである確率０．９０が、外挿されたＮ＝２０．５９（垂直線で示される）で得られる。

各増幅反応物が予測された２０の増幅可能なテンプレート分子を有し、それにより増幅可能な突然変異体分子の存在を有する各反応物では、突然変異体ＡＦ＞０．０３８２を有する確率が０．９を超えるように、ＤＮＡ試料の希釈率が選択された。

試料中でＮ以上の分子を観察するポアソン確率は、Ｎ以上の反応物中でバックグランド率を超える突然変異体配列の同定に基づく複合コールの確率とほぼ等しかった。偽陰性の近似的確率を各試料で計算し、表２および３に示している（図４参照）。

実施例５ − 基本実験を証明するためのＰＣＲプライマーおよび設計
４１のプライマー対を、４１のアンプリコンのパネルについて設計した。各プライマーを、５’末端での標識配列と、中央での７つの異なるバーコード配列の１つと、を有して作製した（図６）。標的特異性プライマーを、２つの最適化されたマルチプレックスＰＣＲプールに分割し、そのフォワードプライマーおよびリバースプライマーを、４９の異なるバーコード組が作成されるようにプールした（図７）。Ｆｌｕｉｄｉｇｍプラットフォームを用いて４８のプライマープールを異なるＰＣＲウェルに分配して、これをハイスループットで実施した。

ＤＮＡをデジタル方式で定量して、平均２０個の分子が各ＰＣＲに投入されるように、ＤＮＡを添加した。Ｆｌｕｉｄｉｇｍ Ａｃｃｅｓｓ Ａｒｒａｙは、６８．３２％のデッドボリュームを有し、それはＤＮＡを添加する場合に考慮された（デッドボリューム＝３３ｎｌ^＊４８／５μｌ）。

これにより、試料あたり約２０個の分子の４８の別個のプールの増幅が可能になった。低頻度アレイが予測された試料では、４８ウェルの複数組で実施した。

この最初のＰＣＲの回収に続き、２回目のＰＣＲをその後実施して、シーケンサーアダプターおよび試料バーコードを結合させた。その後、ライブラリーを洗浄して、Ｉｌｌｕｍｉｎａ ＭｉＳｅｑおよびＨｉＳｅｑ２０００の両シーケンサーでシーケンシングした。

実施例６ − 細胞株の希釈系列における原理証明
ヘテロ接合体突然変異が３〜０．０４％のＡＦで存在すると予測された、細胞株ＤＮＡの希釈系列を作製した。ＮＣＩ−Ｈ１９７５およびＶＣＡＰ細胞株は、実施例５に記載されたプライマー／アンプリコンパネルを用いて検出可能であり、ＬＮＣａＰ細胞株には存在しない５つの既知の突然変異／ＳＮＰを有する（表４）。ＶＣＡＰ、ＮＣＩ−Ｈ１９７５およびＬＮＣａＰ ＤＮＡをデジタル方式で定量して標準化し、その後、ＬＮＣａＰのＤＮＡ中に系列希釈した。

異なる試料について異なる数の反応物を使用し、合計５種の異なる試料で、２種の細胞株の系列希釈物で、この方法を実施した（表２および３（図４）参照）。

２種の細胞株は、合計２５の（フォワードおよびリバースを独立してコールする場合には５０の）陽性について、含まれるゲノムの領域中の５つの一塩基多型によって、ＬＮＣａＰと異なっていた。

バックグランド率をモデル化するために、目的の領域中の各位置での平均ＡＦおよびＣＶを用いて、パラメータを推定し、各遺伝子座および各可能な塩基置換に特異的なベータ分布を適合させた。

先のモデルならびに表２および３（図４）の計算に基づき、ＣＤＦ＿ｔｈｒｅｓｈ＝０．９９９９および複数（Ｎ）の陽性反応物を必要とする複合コールを用いて、以下の結果を得た。

Ｎ≧２
全体的感度０．９６、３６の偽陽性複合コール。表２（図４Ａ）は、異なる希釈率での感度を示す。

Ｎ≧３
全体的感度０．９０、１の偽陽性複合コール。表３（図４Ｂ）は、異なる希釈率での感度を示す。

２種の細胞株からのＤＮＡの希釈系列でこの方法を実施することによる、同定された突然変異（真陽性（ＴＰ））、偽陽性コール（ＦＰ）、および見落とされた突然変異（偽陰性（ＦＮ））を、表４および図９に示す。図９のデータポイントは、測定されたＡＦの増加の順に並べられ、予測ＡＦを示すＦＮ突然変異を除き、各突然変異に関する測定ＡＦを示している。

図１０は、ＣＤＦ＿ｔｈｒｅｓｈ＝０．９９９９およびＮ≧３を用いる複合突然変異コールに基づく予測ＡＦと測定ＡＦとの一致を示す。偽陰性は、垂直軸上の「検出されず」（ＮＤ）として示されている。予測値は、バツで示され、見やすいように点線で結ばれている。

実施例７ − 突然変異コーリングの別の方法による原理証明
この方法を、実施例６に記載された希釈系列で実施した。

ＶＣＡＰおよびＮＣＩ−Ｈ１９７５のＤＮＡを、表６に示される通り、ＬＮＣａＰのＤＮＡ中に系列希釈した。

最初の実験では、７（ｘ４８）ＬＮＣａＰ、１（ｘ４８）ＶＣＡＰ、１（ｘ４８）ＮＣＩ−Ｈ１９７５、そしてその後１つまたは複数の希釈物をシーケンシングした（表６、「繰り返し数」）。

突然変異コーリング
３３６のＬＮＣａＰ ＰＣＲ反応を用いて、各塩基でのバックグランドＡＦを決定した。正規分布を利用して、目的の領域の全ての位置での非参照塩基から可能な全ての塩基についてバックグランドをモデル化した。

各反応物を、特定のｚスコアおよび深度によってバックグランドと異なる各塩基での置換についてその後スクリーニングした。このスコアを超える特定回数が同定された塩基置換が、突然変異を有すると決定された。

結果
２０のｚスコアカットオフおよび３つの陽性ウェルを利用して、突然変異コーリングを最初に実施した。表７に、突然変異検出結果を示している。

突然変異／ＳＮＰの全てが、０．３３％に至るまで検出された。０．３３％希釈でスクリーニングされた推定１，９２０個の分子のうち（２試料×４８ウェル×２０個ＤＮＡ分子）、４と１２の間の陽性ウェル／分子が検出され、それはほぼ厳密に０．３３％に相当した（表７、Ｄｉｌ０．３３）。

０．１１％希釈では、５つの置換のうち４つが検出された。２つのウェルのみが陽性であったため、１つは見落とされた（表７、Ｄｉｌ０．１１）。最後に、０．０３７％希釈では、５つの突然変異のうち１つが完全に見落とされ、２つの重複アンプリコンうちの１つで、もう１つ（ｃｈｒ１７：７５７７６４４ Ｃ＞Ｇ）が見落とされた。これらの結果は、これらの突然変異体分子の無作為分布と適合する（表７、Ｄｉｌ０．０３７）。

２つの偽陽性反応は、早期のライブラリー増幅回中のポリメラーゼ誤差による可能性が最も高かった。そのような誤差は、典型的には通常の変化よりも低い頻度であるはずである。

ｚスコアカットオフを３０に増加させることによって、実際の変化の全てが保持され、２つの偽陽性が排除された（表７）。

Claims (16)

1. ＤＮＡ試料中の目的の領域中の遺伝的多様体を検出するための方法であって、
   （ｉ）所与のシーケンシングプラットフォーム、シーケンシング工程およびシーケンシング深度について、それ以上で増幅反応物のシーケンシング結果において前記遺伝的多様体が観察されるべき閾値頻度を確定して所与の増幅反応物中の前記遺伝的多様体の存在について陽性判定を割り付けることであって、増幅およびシーケンシング誤差により増幅反応物のシーケンシング結果で観察されると予測される遺伝的多様体を裏づけるリード数の分布（リードカウント分布）に前記閾値頻度が基づくこと；
   （ｉｉ）前記遺伝的多様体について、反復増幅反応物中の目的の領域の増幅可能なテンプレート分子の平均数がステップ（ｉ）において決定された閾値頻度の逆数よりも少なくなるように、前記ＤＮＡ試料を複数の反復増幅反応物に分配すること；
   （ｉｉｉ）ステップ（ｉｉ）の増幅反応を実施することおよび増幅反応の産物をシーケンシングすること；
   （ｉｖ）ステップ（ｉ）およびステップ（ｉｉｉ）の結果に基づいて、各反復増幅反応物中の前記遺伝的多様体の有無を決定すること；ならびに
   （ｖ）（ｉｖ）の結果を統合して前記ＤＮＡ試料中の目的の領域中の前記遺伝的多様体の有無を決定すること、
   を含む、方法。
2. ステップ（ｖ）において前記遺伝的多様体が前記ＤＮＡ試料中の目的の領域中に存在すると決定されるために、ステップ（ｉｖ）において前記遺伝的多様体の存在についての陽性判定が１つよりも多くの反復増幅反応物においてなされなければならない、請求項１に記載の方法であって、必要に応じて、ステップ（ｉｖ）において前記遺伝的多様体の存在について陽性判定が少なくとも３つの反復増幅反応物においてなされなければならない、方法。
3. ステップ（ｉ）の閾値頻度が、その頻度で前記遺伝的多様体を観察する確率についての累積分布関数値が０．９９以上、０．９９５以上、０．９９９以上、０．９９９９以上または０．９９９９９以上である頻度である、請求項１または２に記載の方法。
4. ステップ（ｉ）の閾値頻度が、二項分布、過分散二項分布、ベータ分布、正規分布、指数分布またはガンマ確率分布のモデルを用いて決定される、請求項１〜３のいずれか１項に記載の方法。
5. ステップ（ｉ）が、
   （ａ）複数の遺伝的多様体についてのリードカウント分布に基づいて、それ以上で増幅反応物のシーケンシング結果において前記遺伝的多様体が観察されるべき複数の閾値頻度を確定して所与の増幅反応物中の前記遺伝的多様体の存在について陽性判定を割り付けること、および
   （ｂ）ステップ（ａ）に基づいて、それ以上で所与の増幅反応物のシーケンシング結果において遺伝的多様体が観察されるべき全体的閾値頻度を確定して前記増幅反応物中の前記遺伝的多様体の存在について陽性判定を割り付けることであって、前記全体的閾値頻度が、ステップ（ａ）において決定された複数の閾値頻度の少なくとも９０％、９５％、９７．５％、９９％またはそれを超える閾値頻度がこの値未満である閾値頻度であること、
   を含む、請求項１〜４のいずれか１項に記載の方法。
6. 反復増幅反応物中に存在する増幅可能なテンプレート分子の平均数が、反復増幅反応物の１つの増幅可能なテンプレート分子中に前記遺伝的多様体が存在する場合に、その反復測定物について陽性判定がなされる確率が０．９以上である数である、請求項１〜５のいずれか１項に記載の方法。
7. 所与のシーケンシングプラットフォーム、シーケンシング工程およびシーケンシング深度について、増幅およびシーケンシング誤差により増幅反応物のシーケンシング結果で観察されると予測される遺伝的多様体または複数の遺伝的多様体を裏づけるリード数の分布（リードカウント分布）を決定するステップを含み、必要に応じて、リードカウント分布を決定する方法が、ＤＮＡ試料を複数回シーケンシングして遺伝的多様体または複数の遺伝的多様体について前記リードカウント分布を決定することをさらに含む、請求項１〜６のいずれか１項に記載の方法。
8. 前記リードカウント分布が、場合により配列状況を考慮して、シーケンサーおよび／またはポリメラーゼの誤り率に基づいている、請求項１〜６のいずれか１項に記載の方法。
9. データベース、チャート、表、リスト、カタログ、インデックス、ディレクトリ、またはレジスタにおいて遺伝的多様体または複数の遺伝的多様体について前記リードカウント分布のための参照値または複数の参照値を「探索すること」を含み、必要に応じて、前記参照値または複数の参照値が、遺伝的多様体または複数の遺伝的多様体について前記リードカウント分布を決定するために、ＤＮＡ試料を複数回シーケンシングすることにより決定されたものである、請求項１〜５のいずれか１項に記載の方法。
10. ステップ（ｉｉ）における前記ＤＮＡ試料の分配に続いて、各反復増幅反応物が、反復増幅反応物あたり目的の領域について１より多い増幅可能なテンプレート分子を有し、必要に応じて、反復増幅反応物あたりの増幅可能なテンプレート分子の平均数が、１より多く１０００未満、２より多く１０００未満、または５より多く１０００未満である、請求項１〜９のいずれか１項に記載の方法。
11. ２％未満、１％未満、０．５％未満、０．２％未満、０．１％未満または０．０５％未満の頻度でＤＮＡ試料の増幅可能なテンプレート分子の集団内に存在する遺伝的多様体を検出することが可能である、請求項１〜１０のいずれか１項に記載の方法。
12. 前記ＤＮＡ試料を複数の反復増幅反応物に分配することが、前記ＤＮＡ試料を希釈することおよび反復増幅反応物に分取することを含む、請求項１〜１１のいずれか１項に記載の方法。
13. 反復増幅反応が、並行して実施され、かつ増幅反応の産物のシーケンシングが、並行して実施され、
    必要に応じて、増幅反応が、目的の領域に隣接する１つまたは複数のプライマー対を使用して実施されて試料および／または増幅反応反復測定物の特異的識別子配列を増幅の産物中に組み込み、
    さらに必要に応じて、前記プライマーが、配列アダプターを増幅の産物中に組み込む、
    請求項１〜１２のいずれか１項に記載の方法。
14. 複数の目的の領域が、並行して分析される、請求項１〜１３のいずれか１項に記載の方法。
15. （ｖｉ）前記ＤＮＡ試料中の遺伝的多様体の頻度を決定すること、
    をさらに含む、請求項１〜１４のいずれか１項に記載の方法。
16. 対象から得られた試料中の循環する腫瘍ＤＮＡを検出および／または定量するための、
    疾患のモニタリングのインビトロ方法における、
    易罹患性、治療への耐性または応答を予測する遺伝的多様性を同定する方法における、あるいは
    治療への応答をモニタリングするための、
    請求項１〜１５のいずれか１項に記載の方法の利用。
    

Images