JP2019537440A - 配列変異体の検出 - Google Patents

Abstract

核酸配列変異体を検出するための、システム、方法、及び装置が提供される。他の実施形態は、本明細書に記載の方法と関連するシステム及びコンピュータ可読媒体を対象とする。【選択図】図２（ａ）

Description

本開示は、概して核酸増幅によって配列変異体を検出するシステム及び方法に関する。

すべての生物（例えば、動物、植物、及び微生物）の遺伝情報はデオキシリボ核酸（ＤＮＡ）にコードされている。ヒトでは、全ゲノムは２４本の染色体上に約１００，０００個の遺伝子を含有する（Ｔｈｅ Ｈｕｍａｎ Ｇｅｎｏｍｅ， Ｔ． Ｓｔｒａｃｈａｎ， ＢＩＯＳ Ｓｃｉｅｎｔｉｆｉｃ Ｐｕｂｌｉｓｈｅｒｓ， １９９２）。各遺伝子は特定のたんぱく質をコードしており、そのたんぱく質は転写及び翻訳によって発現した後、生細胞内で特定の生化学的機能を果たす。遺伝暗号の変化又は変異は、ｍＲＮＡの配列又は発現レベルを変化させ、潜在的にはｍＲＮＡにコードされたたんぱく質における配列又は発現レベルを変化させる可能性がある。多型又は突然変異として知られるこれらの変化は、疾患をもたらすｍＲＮＡ又はたんぱく質の生物学的活性に重大な有害作用を及ぼし得る。突然変異には、ヌクレオチドの欠失、挿入、置換、又は他の変化（すなわち、点突然変異）が含まれる。

遺伝子多型に起因する多くの疾患が知られており、血友病、サラセミア、デュシェンヌ筋ジストロフィー（ＤＭＤ）、ハンチントン病（ＨＤ）、アルツハイマー病、及び嚢胞性線維症（ＣＦ）が含まれる（Ｈｕｍａｎ Ｇｅｎｏｍｅ Ｍｕｔａｔｉｏｎｓ， Ｄ． Ｎ． Ｃｏｏｐｅｒ ａｎｄ Ｍ． Ｋｒａｗｃｚａｋ， ＢＩＯＳ Ｐｕｂｌｉｓｈｅｒｓ， １９９３）。これらのような遺伝性疾患は、特定の遺伝子を形成するＤＮＡ中における、単一ヌクレオチドの単一の付加、置換、又は欠失から生じ得る。特定の先天異常には、遺伝性疾患を引き起こす変異遺伝子に加え、２１トリソミー（ダウン症候群）、１３トリソミー（パトー症候群）、１８トリソミー（エドワード症候群）、Ｘモノソミー（ターナー症候群）といった染色体異常、及びクラインフェルター症候群（ＸＸＹ）といった他の性染色体異数体が原因となるものがある。さらに、特定のＤＮＡ配列が、糖尿病、動脈硬化症、肥満、種々の自己免疫疾患、及びがん（例えば、結腸直腸、乳房、卵巣、肺）などの多くの疾患のいずれかの素因となり得るという証拠が増加している。

ある種の多型は、個人によっては疾患の素因となるか、又はある種の疾患の罹病率と関連していると考えられている。アテローム性動脈硬化症、肥満、糖尿病、自己免疫障害、及びがんは、多型と相関すると考えられるこのような疾患のいくつかである。疾患との相関に加えて、多型はまた、疾患を治療するために投与される治療薬に対する患者の反応においても役割を果たすと考えられている。例えば、多型は、薬物、放射線療法、及び他の治療形態に反応する患者の能力において役割を果たすと考えられている。

多型の同定によって、特定の疾患のより良い理解、及びそのような疾患に対する潜在的により効果的な治療を導くことができる。実際、患者の同定された多型に基づいて個別化された治療レジメンにより、救命医療介入をもたらすことができる。多型が同定され単離されると、特異的な多型の産物と相互作用するような新規の薬物又は化合物が発見され得る。ウイルス、バクテリア、プリオン、及び真菌を含む感染性微生物の同定もまた、多型に基づいて達成することができ、且つ適切な治療応答を感染宿主へ投与することができる。

約１６ヌクレオチドの配列は、ヒトゲノムのサイズに対してさえも統計的根拠に基づいて特異的であるため、比較的短い核酸配列を用いて、高等生物における正常及び欠陥遺伝子を検出し、感染性微生物（例えば、バクテリア、真菌、原生生物、及び酵母）並びにウイルスを検出することができる。ＤＮＡ配列はまた、同一種内の異なる個体を検出するためのフィンガープリントとしても機能する（Ｔｈｏｍｐｓｏｎ， Ｊ． Ｓ． ａｎｄ Ｍ． Ｗ． Ｔｈｏｍｐｓｏｎ， ｅｄｓ．， Ｇｅｎｅｔｉｃｓ ｉｎ Ｍｅｄｉｃｉｎｅ， Ｗ．Ｂ． Ｓａｕｎｄｅｒｓ Ｃｏ．， Ｐｈｉｌａｄｅｌｐｈｉａ， Ｐａ． （１９９１）を参照されたい）。

核酸配列を検出及び分析するための多数の方法が開発されている。例えば、核酸配列は、増幅された核酸分子の移動度をゲル電気泳動によって既知の標準と比較するか、又は同定される配列に相補的なプローブとのハイブリダイゼーションによって、同定することができる。核酸配列を分析する当技術分野で既知の方法は、規定された核酸配列を酵素的に合成又は増幅するためのインビトロ法であるポリメラーゼ連鎖反応（ＰＣＲ）を用いる。この反応は、典型的には、反対の鎖にハイブリダイズする２つのオリゴヌクレオチドプライマーを使用し、増幅すべき鋳型又は標的ＤＮＡ配列に隣接する。プライマーの伸長は熱安定性ＤＮＡポリメラーゼによって触媒される。鋳型変性、プライマーアニーリング、及びポリメラーゼによるアニーリングされたプライマーの伸長を含む一連の反復サイクルは、特異的ＤＮＡ断片の指数関数的蓄積を提供する。蛍光プローブ又はマーカーは、典型的には、増幅プロセスの検出及び定量を容易にするためにプロセスにおいて使用される。

典型的な動態ＰＣＲ分析において、蛍光強度値は、典型的なＰＣＲプロセスについてサイクル数に対してプロットされ、ＰＣＲ産物の形成はＰＣＲプロセスの各サイクルでモニターされる。増幅は通常、増幅反応中の蛍光シグナルを測定するための成分及びデバイスを含むサーモサイクラーで測定される。そのようなサーモサイクラーの例は、Ｒｏｃｈｅ Ｄｉａｇｎｏｓｔｉｃｓ ＬｉｇｈｔＣｙｃｌｅｒ（カタログ番号第２０１１０４６８番）である。増幅産物はまた、例えば、標的核酸に結合したときのみ蛍光シグナルを発する蛍光標識ハイブリダイゼーションプローブによって、又は特定の場合には、二本鎖ＤＮＡに結合する蛍光色素によっても検出される。

しかしながら、増幅反応の結果として得られるデータに誤差が生じる可能性があり、それによって、データの分析における誤差が引き起こされ、且つ分析から得られるような物理的特性（例えば、出発物質の同定及び／又は定量）が生じる。さらに、核酸の大きな集団における比較的小さな配列変動を分析する場合、これらの誤差は複合され得る。したがって、核酸増幅方法の結果を用いて配列変異体情報を分析するための改善されたプロセスを同定することが有益であろう。

一態様において、試料中の少なくとも２つの核酸配列変異体を検出する方法であって、以下のステップ：（ａ）少なくとも２つの配列変異体を増幅して、第１の変異体アンプリコン及び第２の変異体アンプリコンを含む少なくとも２つの増幅産物を製造し、増幅ステップは第１及び第２のアンプリコンについてそれぞれ第１及び第２の増幅曲線を生じるステップ；（ｂ）プロセッサを用いて、第１の増幅曲線の直線性からの第１の偏差及び第２の増幅曲線の直線性からの第２の偏差を作り、且つ直線性からの第１及び第２の偏差を比較して直線性比から偏差を作ることによって、第１及び第２の増幅曲線の直線性からの相対偏差を分析するステップ；（ｃ）プロセッサを用いて、閾値マトリックスに対する直線性からの相対偏差を比較するステップ；並びに、（ｄ）プロセッサを用いて、比較ステップ（ｃ）の結果に基づく２つ以上の配列変異体を同定するステップ；を含む、方法が提供される。いくつかの実施形態において、該分析ステップ（ｂ）は、第１の増幅曲線についてデルタＢ１を算出すること、及び第２の増幅曲線についてデルタＢ２を算出することを含み、
且つ、直線性からの相対偏差＝デルタＢ１／デルタＢ２である。いくつかの実施形態において、直線性からの相対偏差が、ｌｏｇ１０（デルタＢ１／デルタＢ２）を含む。いくつかの実施形態において、
である。本明細書において、いくつかの実施形態では、ベースライン切片は、ｙ中央値、幾何平均、調和平均、又は一般化平均を含む。ある実施形態において、ベースライン切片は、ｙ中央値を含む。いくつかの実施形態において、
である。いくつかの実施形態において、直線性からの偏差はサイクルｍ〜ｐより評価され、０＜ｍ＜ｐである。ある実施形態において、ｍ＞２である。ある実施形態において、３＜ｍ＜７である、ある実施形態において、ｍ＝６である。

別の態様において、試料中の２つ以上の核酸配列変異体を検出するシステムであって、システムは、メモリ、プロセッサ、及びディスプレイと動作可能に接続された核酸増幅モジュールを含み、プロセッサは、２つ以上の配列変異体を検出するコンピュータ実施方法を実行するように構成され、方法は、以下：（ａ）核酸増幅モジュールを用いて、少なくとも２つの配列変異体を増幅して、第１の変異体アンプリコン及び第２の変異体アンプリコンを含む少なくとも２つの増幅産物を製造するステップ；（ｂ）プロセッサを用いて、第１及び第２の増幅曲線を、それぞれ第１及び第２のアンプリコンについて作成するステップ；（ｃ）プロセッサを用いて、第１の増幅曲線の直線性からの第１の偏差及び第２の増幅曲線の直線性からの第２の偏差を作り、且つ直線性からの第１及び第２の偏差を比較して直線性からの偏差比を作ることによって、第１及び第２の増幅曲線の直線性からの相対偏差を分析するステップ；（ｄ）プロセッサを用いて、閾値マトリックスに対する直線性からの相対偏差を比較するステップ；並びに、（ｅ）プロセッサを用いて、該比較ステップ（ｄ）の結果に基づく２つ以上の配列変異体を同定するステップ；を含む、システムが提供される。いくつかの実施形態において、該分析ステップ（ｃ）は、第１の増幅曲線についてデルタＢ１を算出すること、及び第２の増幅曲線についてデルタＢ２を算出することを含み、
且つ、直線性からの相対偏差＝デルタＢ１／デルタＢ２である。いくつかの実施形態において、直線性からの相対偏差は、ｌｏｇ１０（デルタＢ１／デルタＢ２）を含む。いくつかの実施形態では、
。本明細書において、いくつかの実施形態では、ベースライン切片は、ｙ中央値、幾何平均、調和平均、又は一般化平均を含む。ある実施形態において、ベースライン切片は、ｙ中央値を含む。いくつかの実施形態において、
である。いくつかの実施形態において、直線性からの偏差はサイクルｍ〜ｐより評価され、０＜ｍ＜ｐである。ある実施形態において、ｍ＞２である。ある実施形態において、３＜ｍ＜７である。ある実施形態において、ｍ＝６である。

さらに、本開示は、試料中の２つ以上の核酸配列変異体を検出するシステムを提供し、このシステムは、メモリ、プロセッサ、及びディスプレイと動作可能に接続された核酸増幅モジュールを含む。システムは、核酸増幅モジュールを用いて、少なくとも２つの配列変異体を増幅して、第１の変異体アンプリコン及び第２の変異体アンプリコンを含む少なくとも２つの増幅産物を製造し、プロセッサを用いて、第１及び第２の増幅曲線を、それぞれ第１及び第２のアンプリコンについて作成することによって、２つ以上の配列変異体を検出する方法を実行するように構成される。増幅曲線を、プロセッサを用いて分析し、第１の増幅曲線の直線性からの第１の偏差及び第２の増幅曲線の直線性からの第２の偏差を作り、且つ第１及び第２の偏差を比較して直線性からの偏差比を作ることによって、第１及び第２の増幅曲線の直線性からの相対偏差を決定する。プロセッサを用いて、本システムは、次いで閾値マトリックスに対する直線性からの相対偏差を比較し、それによって２つ以上の配列変異体を同定する。

試料中の２つ以上の核酸配列変異体を検出するためのシステムであって、メモリ、プロセッサ、及びディスプレイと動作可能に接続された核酸増幅モジュールを含むシステムもまた、提供される。核酸増幅モジュールは、少なくとも２つの配列変異体を増幅して、第１の変異体アンプリコン及び第２の変異体アンプリコンを含む少なくとも２つの増幅産物を製造し、プロセッサを用いて、第１及び第２の増幅曲線を、それぞれ第１及び第２のアンプリコンについて作成する。その後、プロセッサは、第１の核酸変異体の増幅のための増幅曲線を表す第１のデータセットを受け取る；さらなる核酸変異体の増幅のための増幅曲線を表す第２のデータセットを受け取り、ここで、該さらなる核酸変異体は、第１の核酸変異体に対する多型部位における該少なくとも１つの配列改変を含んでなる；該第１のデータセットに適合する第１の曲線及び該第２のデータセットに適合する第２の曲線を作成する；第１の曲線の直線性からの第１の偏差及び第２の曲線の直線性からの第２の偏差を生成する；直線性からの第１の偏差と直線性からの第２の偏差とを比較して、直線性からの偏差比を同定する；直線性からの偏差比を閾値マトリックスと比較する；並びに、閾値マトリックスとの比較に基づいて２つ以上の配列変異体を同定する；ように構成される。

さらに、試料中の少なくとも２つの核酸配列変異体を検出する方法であって、以下のステップ：少なくとも２つの配列変異体を増幅して、第１の変異体アンプリコン及び第２の変異体アンプリコンを含む少なくとも２つの増幅産物を製造し、増幅ステップは第１及び第２のアンプリコンについてそれぞれ第１及び第２の増幅曲線を生じるステップ；プロセッサを用いて、第１の増幅曲線の直線性からの第１の偏差及び第２の増幅曲線の直線性からの第２の偏差を作り、且つ直線性からの第１及び第２の偏差を比較して直線性比から偏差を作ることによって、第１及び第２の増幅曲線の直線性からの相対偏差を分析するステップ；プロセッサを用いて、閾値マトリックスに対する直線性からの相対偏差を比較するステップ；並びに、プロセッサを用いて、比較ステップの結果に基づく２つ以上の配列変異体を同定するステップ；を含む、方法が提供される。

また、２つ以上の配列変異体を検出するコンピュータ実施方法であって、コンピュータシステムにおいて、以下：第１の核酸変異体の増幅に関する増幅曲線を表す第１のデータセットを受け取ること；さらなる核酸変異体の増幅に関する増幅曲線を表す第２のデータセットを受け取ることであって、ここで、該さらなる核酸変異体は、第１の核酸変異体に関する多型部位における該少なくとも１つの配列改変を含んでなる；該第１のデータセットに適合する第１の曲線及び該第２のデータセットに適合する第２の曲線を作成すること；第１の曲線の直線性からの第１の偏差及び第２の曲線の直線性からの第２の偏差を作成すること；直線性からの第１の偏差を直線性からの第２の偏差と比較して、直線性からの偏差比を同定すること；直線性からの偏差比を閾値マトリックスと比較すること；並びに、該比較ステップに基づいて２つ以上の配列変異体を同定することを含む、方法も提供される。

本開示の性質及び利点のより良い理解は、以下の詳細な説明及び添付の図面を参照して得ることができる。

図１は、本明細書に記載するシステムのブロック図である。 図２（ａ）は、システムによって使用される方法の概略図である。 図２（ｂ）は、本方法のデータ処理部分の概略図である。 図３は、配列特異的閾値マトリックスの使用の具体例の説明図である。 図４は、従来の配列変異体理論の具体例を示す図である。 図５は、核酸増幅曲線に関するデルタＢのグラフ表示である。 図６は、改善された配列改変方法を示す。 図７は、新規の方法を用いて得られた結果を示す。 図８は、従来の配列変体異理論を用いた偽ＧＧ対立遺伝子の同定を示す。 図９は、従来の配列変異体理論を使用すると、特定の試料が、カットオフを下回ったＣＹ５．５又はＦＡＭ ＲＦＩ（ＡＡのみ）のために無効化されたことを示す。 図１０は、従来の配列変異体理論を使用すると、特定のＧＡ検体が、デルタＣＴ値が範囲外であるために無効化されたことを示す。 図１１は、ＨＥＸ ＲＦＩ値がカットオフ値を超えていたため、ＡＡ検体がどのように無効化されたかを示している。全てのチャネル（ＦＡＭ、ＨＥＸ、及びＣＹ５．５）ＣＴが存在したが、ΔＣＴカットオフに失敗したので、試料をＧＡ試料として処理した。

実施形態は、核酸増幅反応に基づいて配列変異体を検出することができる。従来の方法は配列変異体を不正確に同定するが、本明細書に記載する改善された方法は、従来の方法を用いて正しく検出されなかった変異体を一貫して正確に同定する機構を提供する。さらに、改善された方法は、挑戦的な試料が、例えば潜在的な阻害成分の存在下において評価される場合に、よりロバストであることが見出されている。

本明細書で特に定義しない限り、本明細書で使用される科学的用語及び技術的用語は、当業者によって一般に理解される意味を有する。さらに、文脈によって特に要求されない限り、単数形の用語は複数形を含み、且つ複数形の用語は単数形を含む。冠詞「ひとつの（ａ）」及び「ひとつの（ａｎ）」が本明細書にて使用される場合、冠詞の文法的目的の１つ、又は１つ以上（すなわち、少なくとも１つ）を指す。例として、「要素（ａｎ ｅｌｅｍｅｎｔ）」とは、１つの要素又は１つ以上の要素を意味する。

「検出する（ｄｅｔｅｃｔ）」、「検出すること（ｄｅｔｅｃｔｉｎｇ）」、「検出（ｄｅｔｅｃｔｉｏｎ）」、及び同様の用語は、本出願において、プロセス、又はある物の存在若しくは非存在、及び程度、量、又はレベル、又は発生確率を発見若しくは決定することを広範に指すために使用される。例えば、用語「検出すること（ｄｅｔｅｃｔｉｎｇ）」とは、標的核酸配列に関して使用される場合、その配列の存在、存在、レベル、若しくは量、及び確率若しくは尤度の発見又は決定を示すことがある。「存在又は非存在を検出すること（ｄｅｔｅｃｔｉｎｇ ｐｒｅｓｅｎｃｅ ｏｒ ａｂｓｅｎｃｅ）」、「存在又は非存在の検出（ｄｅｔｅｃｔｉｏｎ ｏｆ ｐｒｅｓｅｎｃｅ ｏｒ ａｂｓｅｎｃｅ）」、及び関連する表現は、定性的且つ定量的検出を含むことを理解されたい。

用語「核酸」、「ポリヌクレオチド」、及び「オリゴヌクレオチド」は、ヌクレオチド（例えば、リボヌクレオチド又はデオキシリボヌクレオチド）のポリマーを指し、天然に存在する（アデノシン、グアニジン、シトシン、ウラシル、及びチミジン）、天然に存在しない、並びに修飾された核酸を含む。この用語は、ポリマーの長さ（例えば、モノマーの数）によって制限されない。核酸は一本鎖であっても又は二本鎖であってもよく、概して５’−３’ホスホジエステル結合を含有し得るが、いくつかの場合において、ヌクレオチド類似体は他の結合を有してもよい。モノマーは、典型的にはヌクレオチドと呼ばれる。用語「非天然ヌクレオチド」又は「修飾ヌクレオチド」とは、修飾された窒素塩基、糖、若しくはリン酸基を含有するヌクレオチド、又はその構造中に非天然部分を組み込むヌクレオチドを指す。非天然ヌクレオチドの例には、ジデオキシヌクレオチド、ビオチン化ヌクレオチド、アミノ化ヌクレオチド、脱アミノ化ヌクレオチド、アルキル化ヌクレオチド、ベンジル化ヌクレオチド、及び蛍光体標識ヌクレオチドが含まれる。

本明細書に記載の標準仕様を含み得る好適な核酸増幅方法としては、これらに限定されないが、ｔｈｅ Ｌｉｇａｓｅ Ｃｈａｉｎ Ｒｅａｃｔｉｏｎ （ＬＣＲ； Ｗｕ Ｄ． Ｙ． ａｎｄ Ｗａｌｌａｃｅ Ｒ． Ｂ．， Ｇｅｎｏｍｉｃｓ ４ （１９８９） ５６０−６９； ａｎｄ Ｂａｒａｎｙ Ｆ．， Ｐｒｏｃ． Ｎａｔｌ． Ａｃａｄ． Ｓｃｉ． ＵＳＡ ８８ （１９９１）１８９−１９３）； Ｐｏｌｙｍｅｒａｓｅ Ｌｉｇａｓｅ Ｃｈａｉｎ Ｒｅａｃｔｉｏｎ （Ｂａｒａｎｙ Ｆ．， ＰＣＲ Ｍｅｔｈｏｄｓ ａｎｄ Ａｐｐｌｉｃ． １ （１９９１） ５−１６）； Ｇａｐ−ＬＣＲ （ＷＯ ９０／０１０６９）；Ｒｅｐａｉｒ Ｃｈａｉｎ Ｒｅａｃｔｉｏｎ （ＥＰ ０４３９１８２ Ａ２）， ３ＳＲ （Ｋｗｏｈ Ｄ．Ｙ． ｅｔ ａｌ．， Ｐｒｏｃ． Ｎａｔｌ． Ａｃａｄ． Ｓｃｉ． ＵＳＡ ８６ （１９８９） １１７３−１１７７；Ｇｕａｔｅｌｌｉ Ｊ．Ｃ．， ｅｔ ａｌ．， Ｐｒｏｃ． Ｎａｔｌ． Ａｃａｄ． Ｓｃｉ． ＵＳＡ ８７ （１９９０） １８７４−１８７８； ＷＯ ９２／０８８０８）、及びＮＡＳＢＡ （米国特許第５，１３０，２３８号明細書）が挙げられる。さらに、以下の方法のうち１つ以上はまた、鎖置換増幅（ＳＤＡ、ｓｔｒａｎｄ ｄｉｓｐｌａｃｅｍｅｎｔ ａｍｐｌｉｆｉｃａｔｉｏｎ）、転写媒介増幅（ＴＭＡ、ｔｒａｎｓｃｒｉｐｔｉｏｎ ｍｅｄｉａｔｅｄ ａｍｐｌｉｆｉｃａｔｉｏｎ）、及びＱｂ増幅（例えば、Ｗｈｅｌｅｎ Ａ． Ｃ． ａｎｄ Ｐｅｒｓｉｎｇ Ｄ． Ｈ．， Ａｎｎｕ． Ｒｅｖ． Ｍｉｃｒｏｂｉｏｌ． ５０（１９９６） ３４９−３７３； Ａｂｒａｍｓｏｎ Ｒ． Ｄ． ａｎｄ Ｍｙｅｒｓ Ｔ． Ｗ．， Ｃｕｒｒ Ｏｐｉｎ Ｂｉｏｔｅｃｈｎｏｌ ４ （１９９３） ４１−４７を概説のために参照されたい）を使用することもできる。

上記のように、特定の実施形態では、とりわけ米国特許第４，６８３，２０２号明細書、同第４，６８３，１９５号明細書、同第４，８００，１５９号明細書、及び同第４，９６５，１８８号明細書に開示されているように、核酸配列はＰＣＲによって分析される。ＰＣＲは、典型的には、選択された核酸鋳型（例えば、ＤＮＡ又はＲＮＡ）に結合する二つ以上のオリゴヌクレオチドプライマーを使用する。核酸分析に有用なプライマーには、標的核酸の核酸配列内で核酸合成の開始点として作用することができるオリゴヌクレオチドが含まれる。プライマーは、従来の方法によって制限消化物から精製することができるか、又は合成的に生産することができる。プライマーは、増幅の最大効率のために好ましくは一本鎖であるが、プライマーは二本鎖であってもよい。二本鎖プライマーは、最初に変性、すなわち鎖を分離するために処理される。二本鎖核酸を変性する１つの方法は、加熱によるものである。「耐熱性ポリメラーゼ」は、熱安定性であるポリメラーゼ酵素、すなわち、鋳型に相補的なプライマー伸長産物の形成を触媒する酵素であり、二本鎖鋳型核酸の変性を行うのに必要な時間、高温に曝されたときに不可逆的に変性しない。一般に、合成は各プライマーの３’末端で開始され、鋳型鎖に沿って５’から３’方向に進行する。熱安定性ポリメラーゼは、例えば、Ｔｈｅｒｍｕｓ ｆｌａｖｕｓ、Ｔ．ｒｕｂｅｒ、Ｔ．ｔｈｅｒｍｏｐｈｉｌｕｓ、Ｔ．ａｑｕａｔｉｃｕｓ、Ｔ．ｌａｃｔｅｕｓ、Ｔ．ｒｕｂｅｎｓ、Ｂａｃｉｌｌｕｓ ｓｔｅａｒｏｔｈｅｒｍｏｐｈｉｌｕｓ、及びＭｅｔｈａｎｏｔｈｅｒｍｕｓ ｆｅｒｖｉｄｕｓから単離されている。それにもかかわらず酵素が補充されるならば、熱安定性でないポリメラーゼをＰＣＲアッセイに用いることもできる。

鋳型核酸が二本鎖である場合、ＰＣＲにおいて鋳型として使用し得る前に、二本鎖を分離する必要がある。鎖の分離は、物理的、化学的、又は酵素的手段を含むあらゆる適切な変性方法によって達成することができる。核酸鎖を分離する１つの方法は、主に変性するまで核酸を加熱することを含む。鋳型核酸を変性させるために必要な加熱条件は、例えば、緩衝塩濃度、変性される核酸の長さ及びヌクレオチド組成に依存するが、典型的には、温度及び核酸長などの反応の特徴に応じて、一度に約９０℃〜約１０５℃の範囲である。変性は典型的には約５秒〜９分間行われる。

二本鎖鋳型核酸が熱によって変性される場合、反応混合物を、各プライマーの標的核酸上のその標的配列へのアニーリングを促進する温度まで冷却させる。アニーリング温度は、好ましくは約３５℃〜約７０℃、さらに好ましくは約４５℃〜約６５℃である。さらに好ましくは約５０℃〜約６０℃、さらに好ましくは約５５〜約５８℃である。アニーリング時間は約１０秒〜約１分（例えば、約２０〜約５０秒；約３０〜約４０秒）であり得る。次いで、反応混合物を、ポリメラーゼの活性が促進又は最適化される温度、すなわち、分析すべき核酸に相補的な産物を生成するためにアニーリングされたプライマーから伸長が起こるのに充分な温度へ調節する。温度は、核酸鋳型にアニーリングされる各プライマーから伸長生成物を合成するのに充分であるべきであるが、その相補的鋳型から伸長生成物を変性させるほど高くすべきではない（例えば、伸長のための温度は、一般に約４０〜８０℃（例えば、約５０℃〜約７０℃；約６０℃）の範囲である）。伸長時間は、約１０秒〜約５分、好ましくは約１５秒〜２分、さらに好ましくは約２０秒〜約１分、さらに好ましくは約２５〜約３５秒であり得る。

新たに合成された鎖は二本鎖分子を形成し、反応の次のステップで利用できる。ストランド分離、アニーリング、及び伸長のステップは、標的核酸に対応する所望の量の増幅産物を生成するのに必要な回数だけ繰り返すことができる。反応の制限因子は、反応に存在するプライマー、熱安定性酵素、及びヌクレオシド三リン酸の量である。サイクリングステップ（すなわち、変性、アニーリング及び伸長）は、好ましくは少なくとも一回繰り返される。検出における使用のために、サイクリングステップの数は、例えば、試料の性質に依存する。試料が核酸の複雑な混合物である場合、検出に充分な標的配列を増幅するためにより多くのサイクルステップが必要とされる。一般的に、サイクリングステップは約２０回繰り返されるが、２０回未満、又は４０回、６０回、若しくは１００回でも繰り返されてもよい。さらに、アニーリング及び伸長のステップを同じステップ（１ステップＰＣＲ）で行うか、又は別々のステップ（２ステップＰＣＲ）で行うＰＣＲアッセイを行うことができる。

用語「試料」とは、個体由来の核酸を含有するか、又は含有すると推定されるあらゆる組成物を指す。この用語は、細胞、組織、又は血液及び血液成分の精製又は分離された成分、例えば、ＤＮＡ、ＲＮＡ、たんぱく質、無細胞部分、細胞溶解物、血漿、血小板、血清、軟膜、及び乾燥血液スポットを含む。試料はまた、細胞株を含む個体から得られた細胞のインビトロ培養物の成分及び成分を含み得る。試料の特定のさらなる例は、糞便、粘膜スワブ、組織吸引物、組織ホモジネート、細胞培養物及び細胞培養上清（真核細胞及び原核細胞の培養を含む）、尿、唾液、痰、及び脳脊髄試料を含む。

本願明細書に記載される方法は、核酸配列変異体を同定するために使用することができる。対立遺伝子、多型、ハプロタイプ、及び／又は遺伝子型を含むがこれらに限定されない、種々の型の配列変異体を同定することができる。用語「対立遺伝子」は、遺伝子の配列変異体を指す。１つ以上の遺伝的差異が対立遺伝子を構成することができる。用語「多型」は、特定のヌクレオチド配列の複数の変異体が集団中に見出され得る状態を指す。多型位置は、変異体を区別する多型ヌクレオチドが存在する核酸配列中の部位を指す。「一塩基多型」又はＳＮＰは、単一ヌクレオチドからなる多型部位を指す。用語「ハプロタイプ」は、個体中の同一染色体上の異なる場所（遺伝子座）にある対立遺伝子の組合せを指し、用語「遺伝子型」は、個体又は試料に含まれる遺伝子の対立遺伝子の合計を指す。

本明細書に記載される方法において、配列変異体は、１つ以上の多型部位に存在するヌクレオチドを検出することによって同定することができる。対象核酸を含有するあらゆるタイプの組織を使用することができる。配列の異常を識別するための多くの方法が当該技術分野で知られており、使用される特定の方法は本開示の重要な側面ではない。例示的な方法は、増幅された核酸の対立遺伝子特異的増幅又はプローブベースの検出を含む。

１つの実施形態において、対立遺伝子は、対立遺伝子特異的増幅又はプライマー伸長方法を用いて同定することができる。対立遺伝子は、プライマーを伸長するＤＮＡポリメラーゼの能力に対する末端プライマーミスマッチの阻害効果に基づく。対立遺伝子特異的増幅又は伸長ベースの方法を用いて対立遺伝子配列を検出するために、対立遺伝子遺伝子のそれぞれに相補的なプライマーは、３’末端ヌクレオチドが多型位置でハイブリダイズするように選択される。標的対立遺伝子の存在下では、プライマーは３’末端の標的配列と一致し、伸長される。別の対立遺伝子の存在下では、プライマーは標的配列に対して３’ミスマッチを有し、プライマー伸長は除去されるか、又は有意に減少する。対立遺伝子特異的増幅ベース方法又は伸長ベース方法は、例えば、米国特許第５，１３７，８０６号明細書、同第５，５９５，８９０号明細書、同第５，６３９，６１１号明細書、及び同第４，８５１，３３１号明細書に記載されている。対立遺伝子特異的増幅ベースの遺伝子型決定方法は、ハプロタイプの同定を容易にすることができる。本質的に、対立遺伝子特異的増幅は、ヘテロ接合試料中の２つの対立遺伝子のうちの１つから複数の多型部位を含む領域を増幅するために使用される。次いで、増幅された配列内に存在するＳＮＰ変異体を、プローブハイブリダイゼーション又は配列決定などによって同定する。

別の実施形態では、プローブと複数の標的対立遺伝子を含むその対応する標的配列との間に形成されるハイブリダイゼーション二本鎖の安定性の差に依存するプローブベースの方法が使用される。充分に厳密なハイブリダイゼーション条件下では、安定な二本鎖はプローブとその標的対立遺伝子配列の間でのみ形成され、他の対立遺伝子配列とは形成されない。安定なハイブリダイゼーション二本鎖の存在は、多くの周知の方法のいずれかによって検出することができる。特定の実施形態において、多型部位を含む標的対立遺伝子由来の複数の核酸配列は、充分に厳密なハイブリダイゼーション条件下で増幅され、プローブのセットにハイブリダイズされる。存在する対立遺伝子は、増幅された標的配列へのプローブの結合パターンから推測される。この態様において、増幅は、プローブハイブリダイゼーションによる分析のために充分な核酸を提供するために行われる。したがって、プライマーは、多型部位を含む種々の対立遺伝子の領域が、試料中に存在する対立遺伝子に関係なく増幅されるように設計される。対立遺伝子非依存性増幅は、対立遺伝子の保存領域にハイブリダイズするプライマーを用いて達成される。当業者は、典型的には、増幅システムの実験的最適化が有用であることを認識するであろう。

プローブベースの遺伝子型同定は、「ＴａｑＭａｎ」又は「５’ヌクレアーゼアッセイ」を用いて、米国特許第５，２１０，０１５号明細書、同第５，４８７，９７２号明細書、及び同第５，８０４，３７５号明細書、並びにＨｏｌｌａｎｄ ｅｔ ａｌ．， １９８８， Ｐｒｏｃ． Ｎａｔｌ． Ａｃａｄ． Ｓｃｉ． ＵＳＡ， ８８：７２７６ ７２８０に記載される通り、実行することができる。ＴａｑＭａｎアッセイにおいて、増幅された領域内でハイブリダイズする標識検出プローブは、増幅反応混合物中に添加される。プローブは、ＤＮＡ合成のためのプライマーとして作用しないように修飾される。増幅は、５’から３’エキソヌクレアーゼ活性を有するＤＮＡポリメラーゼ、例えばＴｔｈ ＤＮＡポリメラーゼを用いて行われる。増幅の各合成段階の間に、伸長されるプライマーから下流の標的核酸にハイブリダイズするあらゆるプローブは、ＤＮＡポリメラーゼの５’から３’へのエキソヌクレアーゼ活性によって分解される。このように、新しい標的鎖の合成もプローブの分解をもたらし、分解産物の蓄積が標的配列の合成の尺度となる。

特定の態様において、検出プローブは２つの蛍光色素で標識され、そのうちの１つは、他の色素の蛍光を消光することができる。色素は、好ましくは５’末端に付着された１つのプローブ及び内部部位に付着され、その結果、プローブがハイブリダイズしていない状態にあるときに消光が起こり、ＤＮＡポリメラーゼの５’から３’へのエキソヌクレアーゼ活性によるプローブの切断が２つの色素の間で起こる。増幅の結果、色素間でプローブが切断され、同時に消光が除去され、最初に消光された色素から観察できる蛍光が増加する。分解生成物の蓄積は、反応蛍光の増加を測定することによってモニターされる。米国特許第５，４９１，０６３号明細書及び同第５，５７１，６７３号明細書は、増幅と同時に起こるプローブの分解を検出するための代替方法を記載している。

ＴａｑＭａｎアッセイは、プローブが増幅産物の存在を検出するためにのみ使用されるように、対立遺伝子特異的増幅プライマーと共に使用することができる。このようなアッセイは、上述した動態ＰＣＲベースの方法について記載したように実施する。或いは、ＴａｑＭａｎアッセイを標的特異的プローブと共に用いることができる。

基礎となるアッセイ法にかかわらず、本明細書中に記載されるシステム及び方法は、分析された各配列についてのＰＣＲ増幅曲線の直線性からの相対偏差を比較し、その値を、試料中に存在し得る配列変異体について規定された閾値マトリックスと比較することによって、試料中に存在する特異的な配列変異体を同定する。

図１は、核酸増幅モジュール１００、メモリ１０１、プロセッサ１０２、及びディスプレイ１０３を含む標的核酸配列中の少なくとも１つの配列変異体を検出するためのシステム１０４のブロック図を示す。核酸増幅モジュール１０１は、１つ以上の試料処理モジュール１０５を含む。各試料処理モジュール１０５は、分析のための１つ以上の核酸配列を含む試料を処理するために必要な種々のステップを実施するための１つ以上のユニット又はステーションを含む。試料処理モジュール１０５は、例えば温度サイクラーといった、反応チャンバ１０６及び熱電冷却装置１０７、並びに随意には、試料分注ステーション、分離ステーション、及び１つ以上の消耗品及び／又は試薬貯蔵ステーションのうちの１つ以上（示さず）を含む。反応チャンバは、１つ以上の核酸増幅反応ステップ中に試料を収容するように構成される。さらに、核酸増幅モジュール１０１は、１つ以上の試料処理モジュールに電気的に接続された少なくとも１つの制御ユニット１０８も含む。制御ユニット１０８はまた、検出可能なシグナルを得るために核酸を分析するように構成された分析モジュール１０９を含む。

メモリ１０２は、あらゆるタイプの揮発性又は不揮発性メモリ、例えばランダムアクセスメモリ（ＲＡＭ）、リードオンリーメモリ、例えば電気消去可能プログラマブル読み出し専用メモリ（ＥＥＰＲＯＭ）、フラッシュメモリ、ハードドライブ、ソリッドステートドライブ、光学ディスク、及び同種のものといったあらゆる組み合わせを含むことができる。簡潔にするために、メモリ１０２は、図１において単一のデバイスとして示されているが、メモリ１０２は、複数のデバイスに分散されてもよいことが理解される。プロセッサ１０３は、中央処理装置（ＣＰＵ）、グラフィックス処理装置（ＧＰＵ）、特殊目的信号又は画像プロセッサ、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）、テンソル処理装置（ＴＰＵ）、及び同種のものなど、あらゆるタイプの１つ以上のプロセッサを含むことができる。図１では、プロセッサ１０３は単一のデバイスとして示されているが、プロセッサ１０３はまた、あらゆる数のデバイスに分散されてもよい。ディスプレイ１０４は、ＬＣＤ、ＬＥＤ、ＯＬＥＤ、ＴＦＴ、プラズマなどのあらゆる好適な技術を使用して実装され得る。いくつかの実装において、ディスプレイ１０４は、タッチセンシティブディスプレイ（タッチスクリーン）であり得る。

また、システム１００は、デスクトップコンピュータ、ラップトップコンピュータ、タブレット、スマートフォン、サーバ、特定用途向けコンピューティングデバイス、又は本明細書に記載する技術及び動作を実行することができるあらゆる他のタイプの電子デバイスなどの、１つ以上のコンピューティングデバイス（示さず）に動作可能に接続することもできる。いくつかの実施形態では、システムの要素及び各要素のサブコンポーネントは、本明細書で説明する様々な機能を実現する、単一のデバイス又は２つ以上のデバイスの組み合わせとして提供することができる。例えば、核酸増幅モジュール１０１は、１つ以上のローカルエリアネットワーク及び／又はワイドエリアネットワークを介して互いに通信可能に結合された、１つ以上のサーバコンピュータ及び１つ以上のクライアントコンピュータを含むことができる。そして、システム１００はまた、１つ以上の周辺装置１１０（例えば、図１に示すようなプリンタ１１１及びキーボード１１２）を含むことができ、コンピュータサブシステムはシステムバスを介して相互接続することができる。Ｉ／Ｏコントローラに結合する周辺装置及び入出力（Ｉ／Ｏ）装置は、シリアルポートなどの当技術分野で公知の任意の手段によってシステム１００に接続することができる。例えば、シリアルポート又は外部インターフェース（例えば、Ｅｔｈｅｒｎｅｔ、Ｗｉ−Ｆｉなど）を使用して、システム１００をインターネット、マウス入力装置、又はスキャナなどの広域ネットワークに接続することができる。システムバスを介した相互接続により、中央プロセッサ１０３は、各サブシステムと通信し、システムメモリ１０２又は記憶装置（示さず）からの命令の実行、並びにサブシステム間の情報交換を制御することができる。システムメモリ及び／又は記憶装置は、コンピュータ可読媒体を具現することができる。本明細書に記載のあらゆる値は、１つのコンポーネントから別のコンポーネントに出力することができ、ユーザに出力することができる。

コンピュータシステムは、例えば、外部インターフェース又は内部インターフェースによって互いに接続された、複数の同じ構成要素又はサブシステムを含むことができる。いくつかの実施形態では、コンピュータシステム、サブシステム、又は装置は、ネットワークを介して通信することができる。このような場合、１台のコンピュータをクライアント、もう１台のコンピュータをサーバと見なすことができ、各コンピュータを同じコンピュータシステムの一部とすることができる。クライアントとサーバはそれぞれ、複数のシステム、サブシステム、又はコンポーネントを含むことができる。

当然のことながら、本開示の実施形態のいずれも、ハードウェア（例えば、特定用途向け集積回路又はフィールドプログラマブルゲートアレイ）を使用する制御ロジックの形態、及び／又はモジュラ方式若しくは統合方式で一般にプログラム可能なプロセッサを備えたコンピュータソフトウェアを使用する制御ロジックの形態で実施することができる。本明細書中のユーザとして、プロセッサは、同一の集積チップ上のマルチコアプロセッサ、又は単一の回路基板若しくはネットワーク上の複数の処理ユニットを含む。本明細書で提供される開示及び教示に基づいて、当業者は、ハードウェア及びハードウェアとソフトウェアの組み合わせを使用して、本開示の実施形態を実施するための他の方法及び／又は方法を知っており、かつ理解するであろう。

本出願に記載されたソフトウェア構成要素又は機能のいずれも、例えば従来の又はオブジェクト指向技術を用いて、例えばＪａｖａ、Ｃ＋＋、又はＰｅｒｌのようなあらゆる好適なコンピュータ言語を用いてプロセッサによって実行されるソフトウェアコードとして実装することができる。ソフトウェアコードは、記憶及び／又は伝送のため、コンピュータ可読媒体上に一連の命令又はコマンドとして記憶することができ、適切な媒体には、ランダムアクセスメモリ（ＲＡＭ）、リードオンリーメモリ（ＲＯＭ）、ハードドライブ、若しくはフロッピーディスクなどの磁気媒体、又はコンパクトディスク（ＣＤ）若しくはＤＶＤ（デジタル汎用ディスク）などの光学媒体、フラッシュメモリ及び同種のものがある。コンピュータ可読媒体は、このような記憶装置又は送信装置の任意の組み合わせであってよい。

このようなプログラムはまた、インターネットを含む様々なプロトコルに準拠する有線、光、及び／又は無線ネットワークを介して送信するように構成されたキャリア信号を用いて符号化及び送信することもできる。したがって、本開示の一実施形態によるコンピュータ可読媒体は、そのようなプログラムで符号化されたデータ信号を使用して生成されてもよい。プログラムコードでコード化されたコンピュータ可読媒体は、互換性のある装置と共にパッケージ化されてもよいか、又は他の装置（例えば、インターネットからのダウンロード）とは別に提供されてもよい。このようなコンピュータ可読媒体は、単一のコンピュータプログラム製品（例えば、ハードドライブ、ＣＤ、又はコンピュータシステム全体など）上又は内部に存在してもよく、システム又はネットワーク内の異なるコンピュータプログラム製品上又は内部に存在してもよい。コンピュータシステムは、本明細書に記載するあらゆる結果をユーザに提供するためのモニタ、プリンタ、又は他の適切なディスプレイを含むことができる。

本明細書に記載の方法のいずれも、ステップを実行するように構成することができる１つ以上のプロセッサを含むコンピュータシステムを用いて、全体的に又は部分的に実行することができる。したがって、実施形態は、本明細書に記載される方法のいずれかのステップを実行するように構成されたコンピュータシステムに向けられ、場合によっては、異なる構成要素がそれぞれのステップ又はそれぞれのステップのグループを実行する。番号を付したステップとして示されているが、本明細書における方法のステップは、同時に、又は異なる順序で実施することができる。さらに、これらのステップの一部は、他の方法からの他のステップの一部と共に使用されてもよい。また、ステップの全部又は一部は随意であってもよい。さらに、あらゆる方法のあらゆるステップを、これらのステップを実行するためのモジュール、回路、又は他の手段を用いて実行することができる。

図２（ａ）は、標的核酸配列の配列変異体を検出するために図１のシステムによって使用される方法を示す。試料が受け取られ、適切に処理されると、例えば、試料中に存在する種々の核酸を、核酸増幅反応（２０１）を行うために必要な１つ以上のプライマー及びプローブ並びに付加的反応成分と組み合わせると、試料中の核酸は増幅反応に供され、各核酸は、検出可能な部分、例えば、異なる蛍光色素（２０２）を用いて同定される。特定の実施形態では、少なくとも２つの核酸、例えば、第１の変異体及び第２の変異体が試料中に存在し、それぞれは、異なる検出可能な部分を用いて増幅反応において一意的に同定可能である。随意に、この方法は、増幅反応において差次的に検出することができるようにユニークに標識された対照物質、例えばハウスキーピング遺伝子も含む。システムは、各増幅反応について増幅曲線（ＧＣ；２０３）を生成し、次いで、プロセッサは、各増幅曲線の直線性（ＤＬ；２０４）からの偏差を決定する。直線性からの偏差を互いに比較して直線性からの偏差比（ＤＬＲ；２０５）、すなわち各データセットの直線性からの相対的偏差の比較を同定し、次いで直線性からの偏差比を標的核酸配列に特異的な予め規定された閾値マトリックス（２０６）と比較して、試料中に存在する特異的配列変異体を同定する（２０７）。

ステップ２０３〜２０７は、本方法のデータ解析構成要素（２０８）を含み、これらのステップは図２（ｂ）により詳細に示される。プロセッサは、第１の変異体（Ｖ１）の第１の核酸増幅反応に関する第１のデータセット（２０９）、及び第２の変異体（２１０）の第２の核酸増幅反応に関する第２のデータセット（Ｖ２）を受け取る。各データセットは、各増幅反応（２１１及び２１２）についての増幅曲線を生成するためにプロセッサによって使用され、その後、各増幅曲線についての直線性からの偏差がプロセッサによって生成される（２１３及び２１４）。増幅曲線の直線性からの偏差を比較することによって、プロセッサは、直線性からの偏差比（２１５）を生成し、次いで、これを（２１６）配列特異的閾値マトリックス（ＴＭ）（２１７）と比較して（２１８）、試料中に存在する配列変異体を同定する。

配列特異的閾値マトリックス（２１７）は、複数の既知の特異的核酸配列の個々の増幅曲線を分析することによって決定される。例えば、２つの変異体を評価する場合、各変異体の試料のプールを分析して、これらの配列について期待される増幅曲線を同定する。次いで、他のものに対するそれぞれの可能な対象の変異体についての直線性からの偏差比が決定され、このデータが配列特異的閾値マトリックスを生成するために使用される。さらに、閾値マトリックスは、限定されるものではないが、デルタＣＴＲ又はデルタＣＴなどの、増幅曲線の他の特徴の評価も含むことができる。サイクル閾値（ＣＴ）は、検出可能な信号が閾値を超え、バックグラウンドレベルを超え、増幅が開始したことを示すのに必要なサイクル数である。ＣＴは反応に存在する核酸の量に特異的である。デルタＣＴは、評価される核酸のＣＴと参照配列のＣＴとの間の差に対応する。増幅曲線の形状は、対象とする特定の配列、プライマー、プローブ、及び熱サイクル条件に特異的であり得る。閾値マトリックスは、評価中の各変異体に対する直線性からの偏差比の許容可能な範囲のセットを含む。

例えば、以下の実施例においてより詳細に記載されるように、閾値マトリックスは、図３に示されるように、配列、例えば、ＧＧ、ＧＡ及びＡＡにおける異なる対立遺伝子の間を区別するために使用され得る。核酸増幅は、核酸増幅モジュールにおける３つのチャネルにおいて行われ、各チャネルは、異なる検出可能な標識を使用して異なる対立遺伝子を検出するように構成される。例えば、第１のチャネルは、例としてＣｙ５．５を用いて識別される制御を含むことができ、第２のチャネルは、Ｇの存在を評価するために使用され、ＨＥＸで検出可能に標識され、第３のチャネルは、Ａの存在を評価するために使用され、ＦＡＭで検出可能に標識される。（特異的対立遺伝子及び検出可能な標識の選択は、限定的なものではなく、例示目的のためにのみ本明細書に提供される。）閾値マトリックスは、各対立遺伝子の直線性からの偏差比に対する許容範囲のセットを含むので、この例において評価される３つの対立遺伝子について、マトリックスは、ＧＧ範囲、ＧＡ範囲、及びＡＡ範囲を含む。２つの範囲の間の接合部において、直線性からの偏差比がその範囲内にある場合、対立遺伝子の同一性をさらに確認するために１つ以上のさらなる計算を使用することができる、マトリックスの領域である緩衝帯と呼ばれる範囲があり得る。

各チャネル（３０１）における増幅の後に経験的に有効な結果が得られる場合、例えば、制御チャネルに対する増幅曲線ＣＴが特定の範囲外にないか又は存在しない場合、ＦＡＭ及びＨＥＸに対するＣＴは両方とも存在しない場合、及びＦＡＭ及びＨＥＸに対する個々のＣＴが特定の範囲外にない場合、直線性比は以下のように評価される：
●直線性比が１．０以上（３０２）であれば、対立遺伝子はＧＧである；
●直線性比が−０．２〜０．６の間で、デルタＣＴが存在する場合（３０３）、対立遺伝子はＧＡである；
●ＧＡ緩衝地帯における直線性比が０．６〜１．０の間であり、デルタＣＴが−４．０〜１．０の間である場合（３０４）、対立遺伝子はＧＡである；
●ＡＡ緩衝地帯における直線性比が−０．６〜−０．２の間にあり、デルタＣＴが３．０以上であるか、又は欠落している（３０５）場合、対立遺伝子はＡＡである；又は、
●直線性比が−０．６（３０６）以下の場合、対立遺伝子はＡＡである。

特定の実施形態では、方法は、一定範囲のサイクル数にわたる核酸増幅データからの直線
の直線性からの相対最大偏差として定義されるデルタＢの計算を含む。したがって、デルタＢを使用して、例えば以下の通りに、核酸ｎについての増幅曲線の直線性からの偏差を計算することができる；
式中、ｎは分析される核酸に対応し、ベースラインの一般化された平均は、ｙ中央値、幾何平均、調和平均、又は一般化された平均を含み、特定の実施形態では、ベースラインの一般化された平均はｙ中央値である。したがって、デルタＢは次のように表すことができる；

直線性からの偏差は、好ましくは、サイクル数の範囲、例えばサイクルｍ〜ｐまでで評価され、０＜ｍ＜ｐである。一実施形態では、ｍ≧２、特に３≦ｍ≦７、好ましくはｍ＝６である。

したがって、この実施形態では、第１及び第２の核酸の増幅を表す２つの増幅曲線の直線性からの相対的近似偏差を反映する直線性からの偏差比は、以下のように計算することができる；
直線性からの偏差比＝デルタＢ₁／デルタＢ₂。

特定の実施形態において、直線性からの偏差比は、次を含む：
ｌｏｇ１０（デルタＢ１／デルタＢ２）。

したがって、直線性からの偏差比の一例は、ＬｏｇデルタＢ率（ＬＤＢＲ、Ｌｏｇ ＤｅｌｔａＢ Ｒａｔｉｏ）である。驚くべきことに、ＬＤＢＲのような直線性からの偏差比の使用は、ＲＦＩ（相対蛍光強度）を従来のアルゴリズムと比較して低下させることを可能にし、Ｓａｎｇｅｒ配列決定との比較に基づいて従来の方法よりも高い感度で配列変異体を同定する。ＬＤＢＲは、従来のロジックを使用して遺伝子型決定エラーを導入することが分かっている、１つ又は複数の材料若しくは条件で試料を調べる場合に、無効又は不正確な結果ではなく正しい遺伝子型を同定する上でより頑健であることが分かった。

したがって、本願明細書に記載されるシステム及び方法において、方法は、以下のステップ：（ａ）第１の核酸変異体の増幅のための増幅曲線を表す第１のデータセットを受信するステップ；（ｂ）さらなる核酸変異体の増幅のための増幅曲線を表す第２のデータセットを受け取るステップあって、該さらなる核酸変異体が、第１の核酸変異体に対する多型部位における該少なくとも１つの配列改変を備えるステップ；（ｃ）該第１のデータセットに適合する第１の曲線及び該第２のデータセットに適合する第２の曲線を生成するステップ；（ｄ）第１の曲線のデルタＢ及び第２の曲線のデルタＢを生成すること；（ｅ）ＬＤＢＲを導出するために、第１及び第２の曲線についてのデルタＢの比のｌｏｇ₁₀を取ることによって、第１及び第２の曲線についてのデルタＢを比較するステップ；（ｆ）ＬＤＢＲを閾値マトリックスと比較するステップ；並びに、（ｇ）該比較ステップ（ｆ）に基づいて、２つ以上の配列変異体を同定するステップ；を含むプロセッサによって実行される。

特定の実施形態の特定の詳細は、あらゆる好適な方法で組み合わせることができる。しかしながら、本開示の他の実施形態は、各個々の態様に関する特定の実施形態、又はこれらの個々の態様の特定の組み合わせを対象とすることができる。

本開示の好ましい実施形態の上記の説明は、例示及び説明の目的で提示されている。これは、網羅的であること、又は開示を記載された正確な形態に限定することを意図しておらず、多くの修正及び変形が、上記の教示に照らして可能である。実施形態は、本開示の原理及びその実用的な用途を最もよく説明するために選択及び説明され、それによって、当業者が、考えられる特定の用途に適しているような様々な実施形態及び様々な修正で本開示を最もよく利用することを可能にした。

遺伝子型決定アッセイのためのアルゴリズム選択及びカットオフ決定
従来の遺伝子型論理はストレス条件下で不正確な遺伝子型を同定することが分かった（例えば、ＰＣＲ阻害剤を試料に添加した場合）。これらの不正確な結果は、ｃｏｂａｓ（商標）４８００システムにおけるＦＡＭ（トレオニンリンカーを有する６−カルボキシフルオレセイン）又はＨＥＸ（トレオニンリンカーを有するヘキサクロロ−フルオレセイン）チャネルのいずれかにおけるＲＦＩ（相対蛍光強度）値の変動に起因し、ＲＦＩカットオフを越える結果となった。従来の方法を使用して、ＲＦＩ値がカットオフを下回る場合、ＣＴ値はこのチャネルについて除去された。続いて、遺伝子型をチャンネルのサブセット（例えば、図４参照）のみに存在するＣＴ値で決定した。例えば、ＦＡＭ ＲＦＩ値がヘテロ接合標本ＧＡのカットオフ値を下回った場合、ＨＥＸ及びＣＹ５．５ＣＴのみが存在することになり、他のすべての値が仕様の範囲内であれば、ＧＧ遺伝子型を誤って同定することになる。

よりロバストなアルゴリズムを開発するために、他のアルゴリズムを評価した。図５に示すように、実時間ＰＣＲデータからの直線
の直線性からの相対最大偏差をサイクル数の範囲にわたって定義したデルタＢに基づいて、新しいアルゴリズムを開発した。改良されたアルゴリズムは、ＬｏｇデルタＢ率（ＬＤＢＲ）を使用し、図６にアルゴリズムを概略的に示す。改良されたアルゴリズムを使用して得られた結果を図７に示し、新しいアルゴリズムは、表１に示すように、従来のアルゴリズムに対するＲＦＩの低下を可能にする。

ＣＴは、ＲＦＩ値がカットオフ値未満の試料では存在しないように設定される。

試料調製のためにｃｏｂａｓ（商標）ｘ４８０を用いた試料スクリーニング中に、古いアルゴリズム及び新しいアルゴリズムの両方は、＞９９％の検体についてＳａｎｇｅｒ配列決定と比較して一致した結果をもたらした（表２）。

１つの検体は、Ｓａｎｇｅｒ配列決定を有するＧＧと比較して、遺伝子型不一致ＧＡ／無効（新旧アルゴリズム）であり、１つの検体は、Ｓａｎｇｅｒ配列決定を有するＧＡと比較して、遺伝子型無効／ＧＡ（新旧アルゴリズム）であった。これらの検体の両方は、他の対立遺伝子のいくつかの背景のみを示し、典型的にはホモ接合性又はヘテロ接合性ではなかったので、Ｓａｎｇｅｒ遺伝子型トレースについては決定的ではなかった。さらなる調査は、これら２つの検体が検体収集中に交差汚染された可能性を示した。従来の方法を使用すると、そのような決定的でない結果は、使用者に実験を繰り返すことを強いることになり、検体収集中に検体が交差汚染された場合、使用者は、実験を繰り返すために新鮮な試料を収集しなければならない可能性がある。

異なる試料調製をＳａｎｇｅｒ配列決定と比較すると、古いアルゴリズムと新しいアルゴリズムは、Ｓａｎｇｅｒ配列決定（表３）と１００％の一致した結果を示した。しかし、従来のアルゴリズムでは、ＦＡＭ ＣＴＲ値がＤＮＡ収率の高い試料調製法のカットオフ範囲外であるため、３つの検体を無効としていた。

内部統制（集積回路）は、ＤＮＡの質と量の統制として機能した。貧弱な試料調製をシミュレートするために、ＰＣＲ阻害剤のスパイクイン、試料希釈、異なる溶出緩衝液のシミュレーションなどを含む一連の挑戦的条件を実行し、これらのストレス条件下で、新しいアルゴリズムはＳａｎｇｅｒ配列決定と比較して１００％の一致した結果を示したか、又は貧弱なＤＮＡ量又は質（表４）のために試料を無効化した。しかし、古いアルゴリズムはサンガー配列決定法と比較して９７％の一致した結果しか示さず、新しいアルゴリズムと比較して４６以上の試料（又は１７％）が無効と呼ばれた（表４）。

ＰＣＲ阻害剤（血液、エタノール）のスパイクイン時に、１９のＧＡ試料は古いアルゴリズム（表４）で誤ってＧＧと呼ばれた。これは、ＦＡＭ ＲＦＩ値がＲＦＩ_最小カットオフを下回り、ＦＡＭチャネルからＣＴ値が削除されたためである。このように、試料はＨＥＸ及びＣＹ５．５ＣＴ値のみを示し、ＧＧ（図８）として遺伝子型決定した。新しいアルゴリズムを用いて、これらの試料の遺伝子型を正確に決定した。

さらに多くの資料が、以下の種々の理由により旧遺伝子型アルゴリズムでは無効となった：
●ＲＦＩ_最小を高くすると、ＣＹ５．５又はＦＡＭ（ＡＡのみ）がカットオフ値を下回ったために一部の試料が無効になった（図９）；
●ΔＣＴ値が範囲外であったためＧＡ検体を無効とした（図１０）；
●ＨＥＸ ＲＦＩ値がカットオフ値以上であったため、ＡＡ検体は無効となった。すべてのチャネル（ＦＡＭ、ＨＥＸ、ＣＹ５．５）ＣＴが存在したが、ΔＣＴカットオフに失敗したので、試料をＧＡ試料として処理した（図１１）。

したがって、ＬＤＢＲを用いた新しい遺伝子型アルゴリズムは、ＰＣＲ阻害剤、交差汚染、高及び低ＤＮＡ濃度、及びある種のマスターミックス成分変化（ここに示されていないデータ）で試料を攻撃するとき、無効又は不正確な結果の代わりに正しい遺伝子型を同定する上でよりロバストであることが分かった。古いアルゴリズムを使用すると無効な結果となり、ユーザは実験を繰り返す必要があるが、これはしばしば不可能である。しかしながら、新しいアルゴリズムは、従来の方法を用いて無効であることが見出されたであろう試料を正しく遺伝子型決定した。

Claims (16)

1. 試料中の少なくとも２つの核酸配列変異体を検出する方法であって、以下のステップ：
   （ａ）少なくとも２つの配列変異体を増幅して、第１の変異体アンプリコン及び第２の変異体アンプリコンを含む少なくとも２つの増幅産物を製造し、増幅ステップは第１及び第２のアンプリコンについてそれぞれ第１及び第２の増幅曲線を生じるステップ；
   （ｂ）プロセッサを用いて、第１の増幅曲線の直線性からの第１の偏差及び第２の増幅曲線の直線性からの第２の偏差を作り、且つ直線性からの第１及び第２の偏差を比較して直線性からの偏差の比を作ることによって、第１及び第２の増幅曲線の直線性からの相対偏差を分析するステップ；
   （ｃ）プロセッサを用いて、閾値マトリックスに対する直線性からの相対偏差を比較するステップ；並びに、
   （ｄ）プロセッサを用いて、前記比較ステップ（ｃ）の結果に基づく２つ以上の配列変異体を同定するステップ；
   を含む、方法。
2. 前記分析ステップ（ｂ）が、前記第１の増幅曲線についてデルタＢ₁を算出すること、及び前記第２の増幅曲線についてデルタＢ₂を算出することを含み、
   且つ、前記直線性からの相対偏差＝デルタＢ₁／デルタＢ₂
   である、請求項１に記載の方法。
3. 前記直線性からの相対偏差が、ｌｏｇ₁₀（デルタＢ₁／デルタＢ₂）を含む、請求項１〜２のいずれか一項に記載の方法。
4. である、請求項２〜３のいずれか一項に記載の方法。
5. 前記ベースライン切片が、ｙ中央値、幾何平均、調和平均、又は一般化平均を含む、請求項４に記載の方法。
6. 前記ベースライン切片が、ｙ中央値を含む、請求項４〜５のいずれか一項に記載の方法。
7. である、請求項２〜３のいずれか一項に記載の方法。
8. 直線性からの偏差がサイクルｍ〜ｐより評価され、０＜ｍ＜ｐである、請求項１〜７のいずれか一項に記載の方法。
9. 試料中の２つ以上の核酸配列変異体を検出するシステムであって、システムは、メモリ、プロセッサ、及びディスプレイと動作可能に接続された核酸増幅モジュールを含み、プロセッサは、２つ以上の配列変異体を検出するコンピュータ実施方法を実行するように構成され、方法は、以下：
   （ａ）核酸増幅モジュールを用いて、少なくとも２つの配列変異体を増幅して、第１の変異体アンプリコン及び第２の変異体アンプリコンを含む少なくとも２つの増幅産物を製造するステップ；
   （ｂ）プロセッサを用いて、第１及び第２の増幅曲線を、それぞれ第１及び第２のアンプリコンについて作成するステップ；
   （ｃ）プロセッサを用いて、第１の増幅曲線の直線性からの第１の偏差及び第２の増幅曲線の直線性からの第２の偏差を作り、且つ直線性からの第１及び第２の偏差を比較して直線性からの偏差比を作ることによって、第１及び第２の増幅曲線の直線性からの相対偏差を分析するステップ；
   （ｄ）プロセッサを用いて、閾値マトリックスに対する直線性からの相対偏差を比較するステップ；並びに、
   （ｅ）プロセッサを用いて、前記比較ステップ（ｄ）の結果に基づく２つ以上の配列変異体を同定するステップ；
   を含む、システム。
10. 前記分析ステップ（ｃ）が、第１の曲線についてデルタＢ１を算出すること、及び第２の曲線についてデルタＢ２を算出することを含み、
    且つ、直線性からの相対偏差＝デルタＢ１／デルタＢ２
    である、請求項９に記載のシステム。
11. 前記直線性からの相対偏差が、ｌｏｇ１０（デルタＢ１／デルタＢ２）を含む、請求項９〜１０のいずれか一項に記載のシステム。
12. である、請求項１０〜１１のいずれか一項に記載のシステム。
13. 前記ベースライン切片が、ｙ中央値、幾何平均、調和平均、又は一般化平均を含む、請求項１２に記載のシステム。
14. 前記ベースライン切片が、ｙ中央値を含む、請求項１２〜１３のいずれか一項に記載のシステム。
15. である、請求項１０〜１１のいずれか一項に記載のシステム。
16. 直線性からの偏差がサイクルｍ〜ｐより評価され、０＜ｍ＜ｐである、請求項９〜１５のいずれか一項に記載のシステム。