(57)【特許請求の範囲】
【請求項１】
ポリメラーゼ連鎖反応（ＰＣＲ）プロセスについてのデータが、有意な増幅を示すかを決定する方法であって、当該方法が以下の：
ポリメラーゼ連鎖反応プロセスを示すデータセットを受容し、ここで当該データセットは複数のデータポイントを含み各データポイントが座標値の対を有し；
上記データセットを、以下の式：
【数１】
0006009724.tif 000026
で表される二次多項式にフィットさせることにより、曲線を計算し；
上記曲線について統計的に有意な値を決定し、ここで当該統計的に有意な値が、Ｒ2値であり；
上記有意な値が閾値の約０.９０を超えるかを決定し；
閾値を超えない場合、上記データセットをさらに処理し；そして
閾値を超える場合、上記データセットが有意な増殖を有しないと示すか、及び／又は当該データセットを破棄する
を含む、前記方法。
【請求項２】
前記データセットを処理することが、ＰＣＲデータセットのサイクル閾値(Ｃｔ)を決定することをさらに含む、請求項１に記載の方法。
【請求項３】
前記Ｃｔ値の決定が、以下の：
Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ(ＬＭ)回帰プロセスをダブルシグモイド関数に適用して当該関数のパラメーターを決定することにより、前記データセットにフィットする曲線の近似を計算し；
上記決定されたパラメーターを用いて曲線を正規化して、正規化された曲線をもたらし；そして
上記正規化された曲線を処理して、最大曲率の点を決定し、ここで当該最大曲率の点が、ＰＣＲ曲線のＣｔ値を表す
を含む、請求項２に記載の方法。
【請求項４】
前記データセットにフィットする曲線を計算する前に、当該データセットを正規化することをさらに含む、請求項１〜３のいずれか一項に記載の方法。
【請求項５】
ポリメラーゼ連鎖反応（ＰＣＲ）プロセスについてのデータが、有意な増幅を示すかを決定するためにプロセッサーを制御するコードを含むコンピューターが読み取り可能な媒体であって、当該コードが以下の指令：
ポリメラーゼ連鎖反応（ＰＣＲ）プロセスを示すデータセットを受容し、当該データセットは複数のデータポイントを含み、各データポイントが、座標値の対を有し；
上記データセットを以下の式：
【数２】
0006009724.tif 000027
で表される二次多項式にフィットさせることにより、曲線を計算し；
上記曲線について統計的に有意な値を決定し、ここで当該統計的に有意な値が、Ｒ2値であり；
上記有意な値が閾値の約０.９０を超えるかを決定し；そして
閾値を超えない場合、上記データセットをさらに処理し；そして
閾値を超える場合、上記データセットが有意な増殖を有さないということを示すか及び／又は当該データセットを破棄する
を含む、前記媒体。
【請求項６】
動的ポリメラーゼ連鎖反応(ＰＣＲ)システムであって、以下の：
動的ＰＣＲ増幅曲線を表すＰＣＲデータセットを作成する動的ＰＣＲ分析モジュール、ここで当該データセットは、複数のデータポイントを含み、各々は座標値の対を有し；及び
以下の：
上記ＰＣＲデータセットを以下の式：
【数３】
0006009724.tif 000028
で表される二次多項式にフィットさせることにより、曲線を計算し、
上記曲線について統計的に有意な値を決定し、ここで当該統計的に有意な値が、Ｒ2値であり；
上記値が閾値の約０.９０を超えるかを決定し；そして
閾値を超えない場合、上記ＰＣＲデータセットをさらに処理し；そして
閾値を超える場合、上記ＰＣＲデータセットが有意な増殖を有さないということを示すか、及び/又は当該ＰＣＲデータセットを破棄する
ことにより、ＰＣＲデータセットを処理して、ＰＣＲデータセットが有意な増殖を示すかを決定するように適用されたインテリジェンス・モジュール
を含む、動的ポリメラーゼ連鎖反応(ＰＣＲ)システム。
【請求項７】
前記データセットが、ＰＣＲデータセットのサイクル閾値(Ｃｔ)値を決定することをさらに含む、請求項６に記載のＰＣＲシステム。
【請求項８】
前記Ｃｔ値の決定が、以下の：
Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ(ＬＭ)回帰プロセスを適用することによりデータセットをダブルシグモイド関数にフィットさせて、当該関数のパラメーターを決定する曲線の近似を計算し；
上記決定されたパラメーターを用いて曲線を正規化して、正規化曲線をもたらし；そして
上記正規化された曲線を処理して、最大曲率の点を決定し、ここで最大曲率の点が、ＰＣＲ曲線のＣｔ値を表す
を含む、請求項７記載のＰＣＲシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一般に、シグモイド又は増殖曲線を示すデータを処理するシステム及び方法に関する。特に、本発明は、増殖曲線のデータが、有効又は有意な増殖を表すか又は示すかを決定すること、並びにそのような場合、シグモイド又はポリメラーゼ連鎖反応曲線などの増殖型の曲線におけるエルボー値などの特徴的な転移値(transition value)を決定することに関する。
【背景技術】
【０００２】
ポリメラーゼ連鎖反応(ＰＣＲ)とは、規定の核酸配列を酵素的に合成又は増幅するin vitro法である。この反応は、通常、反対の鎖とハイブリダイズし、そして増幅対象のテンプレート又は標的ＤＮＡ配列の側面に位置する２つのオリゴヌクレオチドプライマーを使用する。プライマーの伸張は、熱安定ＤＮＡポリメラーゼにより触媒される。テンプレートの変性、プライマーのアニーリング、及びポリメラーゼによるアニーリングされたプライマーの伸長を伴う反復的な一連のサイクルにより、結果的に、特定のＤＮＡ断片の指数的増殖がもたらされる。増幅プロセスの検出と定量化を促進するために、蛍光プローブ又はマーカーが一般的にこのプロセスにおいて使用されている。
【０００３】
代表的なリアルタイムＰＣＲ曲線が図１に示されており、ここで、通常のＰＣＲプロセスについて蛍光強度値がサイクル数に対してプロットされている。この場合、ＰＣＲ産物の形成が、ＰＣＲプロセスの各サイクルにおいてモニタリングされる。増幅は、通常、増幅反応の間に蛍光シグナルを計測する構成及び装置を含むサーモサイクラーで計測される。このようなサーモサイクラーの一例が、Ｒｏｃｈｅ Ｄｉａｇｎｏｓｔｉｃｓ ＬｉｇｈｔＣｙｃｌｅｒ(カタログ番号：２０１１０４６８)である。増幅産物は、例えば、標的核酸に結合した場合にのみ蛍光シグナルを発する蛍光標識されたハイブリダイゼーションプローブを用いて検出されるか、又はある場合、二本鎖ＤＮＡに結合する蛍光色素を用いて検出される。
【０００４】
一般的なＰＣＲ曲線では、エルボー値又はサイクル閾値(Ｃｔ)と一般的に呼ばれるベースライン領域の末端で転移点を同定することは、ＰＣＲ増幅プロセスの特徴を理解するためにかなり有用である。Ｃｔ値は、ＰＣＲプロセスの効率の尺度として使用することができる。例えば、一般的に、解析対象のすべての反応について、規定のシグナル閾値を決定し、この閾値に到達するのに必要とされるサイクル数(Ｃｔ)が、標的核酸、並びに標準又はハウスキーピング遺伝子などの参照核酸について決定される。標的核酸及び参照核酸について得られたＣｔ値に基づいて、標的分子の絶対数又は相対的なコピー数を決定できる(Gibsonら、Genome Research 6: 995-1001; Biecheら、Cancer Research 59:2759-2765, 1999; WO 97/46707; WO 97/46712; WO 97/46714)。図１におけるベースライン領域の末端における領域２０のエルボー値は、サイクル数が３０の領域内に位置することになろう。
【０００５】
ＰＣＲ曲線におけるエルボー値は、いくつかの既存の方法を使用することにより決定することができる。例えば、様々な最新の方法は、エルボーの実際の値を、ＡＦＬ(任意の蛍光値)と呼ばれる既定のレベルに蛍光強度が到達する値として決定する。その他の最新の方法は、サイクル数に対する蛍光の２次導関数が最大値に達するサイクル数を使用することができよう。しかしながら、これらの方法は、いずれも欠点を有している。例えば、いくつかの方法は、異常値(ノイズ)データに対してかなり影響を受け易く、そして当該ＡＦＬ値のアプローチは、ベースラインが高いデータセットについて上手く機能しない。図１に示される増殖曲線についてのベースラインストップ(又はベースラインの末端)を決定する従来の方法は、特に高タイターの状況下では上手く機能しないことがある。さらに、これらのアルゴリズムは、通常、多数の(例えば、５０以上)のパラメータを具備しており、これらのパラメータは、明確に定義されておらず、線形依存性を有しており、そして多くの場合に、最適化が不可能であるまではいかなくとも非常に困難である。
【発明の開示】
【０００６】
従って、上記問題及び他の問題を克服するシグモイド型の曲線、及びＰＣＲ曲線などの曲線におけるエルボー値を決定するシステム及び方法を提供することが望ましい。曲線が有効な増殖を示しているか、又は処理資源を消費する前に、データを破棄すべきかを最初のうちに決定することも望ましい。
【０００７】
本発明は、増殖曲線についてのデータが、有効又は有意な増殖を表すか又は示すかを決定し、そしてそのような場合、シグモイド又は増殖型の曲線におけるエルボー値などの特徴的な転移値を決定する新規の効率的なシステム及び方法を提供する。ある実施においては、本発明のシステム及び方法は、ＰＣＲ増幅曲線におけるサイクル閾値(Ｃｔ)を決定するのに特に有用である。
【０００８】
本発明の特定の態様では、シグモイド又は増殖型曲線を示すデータセットを処理して、当該データが有意又は有効な増殖を示すかを決定するために処理される。ある態様では、当該データにフィットする一次又は二次の多項式曲線を決定し、そして近似された曲線について統計的に有意な値が測定される。有意な値が、有意閾値(significance threshold)を超える場合、当該データは、有意又は有効ではない増殖を表していないとみなされる。当該データが、有意又は有効な増殖を表さない場合、当該データセットは破棄されてもよい。当該有意な値が、有意閾値を超えない場合、当該データは、有意又は有効な増殖を表すとみなされる。当該データセットが正確な増殖を表すと決定される場合、当該データはさらに処理されて、シグモイド又は増殖曲線における転移値、例えばベースライン領域の末端又はＰＣＲ増幅曲線のＣｔ値若しくはエルボー値を決定する。ある態様では、増殖プロセスを表すデータ曲線が有意閾値を超え、かつ有効な増殖を表すと判断されるように決定される場合、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ(ＬＭ)回帰プロセスにより決定されるパラメーターを有するダブルシグモイド関数を用いて、当該データセットにフィットする曲線への近似を発見した。パラメーターを決定すると、当該曲線は、決定されたパラメーターのうちの１以上を用いて正規化することができる。正規化した後に、正規化された曲線を処理して、曲線に沿った幾つか又は全ての点で曲線の曲率を決定し、例えばＰＣＲデータセットについてサイクル数に対する曲率を表すデータセット又はプロットをもたらした。最大曲率が生じるサイクル数は、ＰＣＲデータセットのＣｔ値に対応する。次に曲率及び/又はＣｔ値が返され、そして表示されるか又はそうでない場合さらなる処理に用いられることもある。
【０００９】
本発明の１の態様では、コンピューターにより実行される方法が、増殖プロセスのデータが有意な増殖を示すかの決定を提供する。当該方法は典型的に、増殖プロセスを表すデータセットを受容することを含み、ここで、当該データセットは複数のデータポイントを含み、各データポイントが座標値を有し、そして当該データセットにフィットする曲線を計算し、当該曲線が一次又は二次多項式のうちの１つを含む。当該方法は、一般的に、当該曲線についての統計的に有意な値を決定し、当該有意な値が閾値を超えるかを決定し、そして閾値を超えない場合、当該データセットをさらに処理し、そして閾値を超える場合、当該データセットが有意な増殖を有さないということを示し、及び/又は当該データセットを破棄することを含む。１の態様では、当該増殖プロセスは、ポリメラーゼ連鎖反応(ＰＣＲ)プロセスである。ここで、当該データセットを処理することは、さらに、ＰＣＲデータセットのサイクル閾値(Ｃｔ)を決定することを含んでもよい。ある態様では、当該曲線は、動的ポリメラーゼ連鎖反応(kinetic polymerase chain reaction)( ＰＣＲ)プロセスの増幅曲線であり、そしてベースライン領域の末端での点は、動的ＰＣＲ曲線のエルボー又はサイクル閾値(Ｃｔ)を表すベースライン領域の末端の位置である。１の態様では、当該曲線を処理して、曲線に沿った幾つか又は全ての点で曲率を決定するように処理され、ここで最大曲率を有する点は、Ｃｔ値を表す。ある態様では、受容されたデータセットは、１以上の異常値又はスパイクポイントを除くように処理されたデータセットを含む。ある態様では、統計的に有意な値はＲ2値であり、そしてその閾値は、約０.９０よりも高い。１の態様では、統計的に有意な値は、Ｒ2値であり、そしてその閾値は約０.９９である。ある実施態様では、当該方法はさらに、データセットにフィットする曲線を計算する前に、データセットを正規化することを含む。さらに別の実施態様では、Ｃｔ値を決定する方法は、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ(ＬＭ)回帰プロセスをダブルシグモイド関数に適用して、関数のパラメーターを測定し、規定のパラメーターを用いて当該曲線を正規化して正規化曲線をもたらし、そして当該正規化曲線を処理して、最大曲率の点を決定することによって、データセットにフィットする曲線の近似を計算することを含み、ここで当該最大曲率の点は、ＰＣＲ曲線のＣｔ値を表す。
【００１０】
本発明の別の態様によると、増殖プロセスについてのデータが有意な増殖を示すかを決定するためにプロセッサーを制御するためのコードを含むコンピューターが読み取り可能な媒体が提供される。当該コードは、一般的に、増殖プロセスを表すデータセットであって、当該データセットが複数のデータポイントを含み、各データポイントが座標値の対を有する、上記データセットを受容し、そして当該データセットにフィットする曲線であって、当該曲線が一次又は二次の多項式のうちの一つを含む、上記曲線を計算する指示を含む。当該コードは、また一般的に、当該曲線について統計的に有意な値を決定し、当該有意な値が閾値を超えるかを決定し、そして当該閾値を超えない場合、当該データセットをさらに処理し、そして当該閾値を超える場合、当該データセットが有意な増殖を有さないということを示し、及び/又は当該データセットを破棄する指示を含む。１の態様では、当該増殖プロセスは、ポリメラーゼ連鎖反応(ＰＣＲ)プロセスである。ここで、コンピューターが読み取り可能な媒体は、さらに、ＰＣＲデータセットのサイクル閾値(Ｃｔ)を決定する指示をさらに含むことがある。別の態様では、当該曲線は動的ポリメラーゼ連鎖反応(ＰＣＲ)プロセスの増幅曲線であり、そしてベースライン領域の末端でのポイントは、動的ＰＣＲ曲線のサイクル閾値(Ｃｔ)又はエルボー値を表す。１の態様では、当該曲線を処理して、当該曲線に沿った幾つか又は全てのポイントで曲率を計算し、ここで最大曲率を有する点が、Ｃｔ値を表す。ある態様では、当該統計的に有意な値は、Ｒ2値であり、そして閾値は、約０.９０よりも大きい。１の態様では、統計的に有意な値はＲ2値であり、そして当該閾値は、約０.９９である。コンピューターが読み取り可能な媒体のある実施態様では、Ｃｔ値を決定するための指示は、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ(ＬＭ)回帰プロセスをダブルシグモイド関数に適用して、関数のパラメーターを決定することにより、データセットにフィットする曲線の近似を計算し；当該決定されたパラメーターを用いて関数を正規化して正規化された曲線をもたらし；そして正規化された曲線を処理して、最大曲率の点(当該最大曲率の点が、ＰＣＲ曲線のＣｔ値を表す)を決定する指示を含む。当該コンピューターが読み取り可能な媒体の別の実施態様において、当該コードはさらに、当該データセットにフィットする曲線を計算する前にデータセットを正規化する指示をさらに含む。さらに別の実施態様では、当該コードは、増殖プロセスがＰＣＲプロセスである場合、Ｃｔ値を表すデータを出力する指示をさらに含む。
【００１１】
本発明のさらに別の態様では、動的ポリメラーゼ連鎖反応(ＰＣＲ)システムが提供される。当該システムは、一般的に、動的ＰＣＲ増幅曲線を示すＰＣＲデータセットを作成する動的ＰＣＲ分析モジュール（ここで、当該データセットが複数のデータポイントを含み、その各々が座標値の対を有し、ここで当該データセットがサイクル閾値(Ｃｔ)を含む目的の領域においてデータポイントを含む）、並びにＰＣＲデータセットが有意な増殖を示すかを決定するために適用されるインテリジェンス・モジュールを含む。インテリジェンス・モジュールは、一般的に、ＰＣＲデータセットにフィットし、そして一次又は二次多項式のうちの一つを含む曲線を計算し、そして当該曲線について統計的に有意な値を決定することにより、ＰＣＲデータセットを処理する。インテリジェンス・モジュールはまた、当該有意な値が閾値を超えるかを決定し、そして当該閾値を超えない場合、当該ＰＣＲデータセットをさらに処理し、そして当該閾値を超える場合、当該ＰＣＲデータセットが有意な増殖を有さないということを示すか及び/又は当該ＰＣＲデータセットを破棄することにより、ＰＣＲデータセットを処理する。１の態様では、当該データセットを処理することはさらに、ＰＣＲデータセットのサイクル閾値(Ｃｔ)値を決定することをさらに含む。ある態様では、当該曲線は、それにより処理されて、曲線に沿った幾つか又は全ての点で曲率を決定し、ここで最大曲率を有する点はＣｔ値を表す。特定の実施態様では、Ｃｔ値を決定するＰＣＲシステムは、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ(ＬＭ)回帰プロセスを、ダブルシグモイド関数に適用して関数のパラメーターを決定することにより、データセットにフィットする曲線の近似を計算し；決定されたパラメーターを用いて曲線を正規化して、正規化曲線をもたらし；そして正規化された曲線を処理して最大曲率の点(ここで、当該最大曲率の点は、ＰＣＲ曲線のＣｔ値を表す)を決定することを含む。ある態様では、統計的に有意な値は、Ｒ2値であり、そして閾値は約０.９０よりも大きい。１の態様では、統計的に有意な値は、Ｒ2値であり、そして閾値は約０.９９である。ある態様では、当該インテリジェンス・モジュールは、さらにデータセットにフィットする曲線を計算する前に、データセットを正規化するためにさらに適用される。ＰＣＲシステムの特定の実施態様では、動的ＰＣＲ分析モジュールは、動的サーモサイクラー装置に組み込まれており、そして分析モジュールに通信可能なように結合されたプロセッサーを含む。ＰＣＲシステムのある実施態様では、インテリジェンス・モジュールは、ネットワーク接続又は直接接続の内の一つにより、分析モジュールにつながれたコンピューターシステム内に存在するプロセッサーを含む。
【００１２】
図面及び特許請求の範囲を含めて、本明細書の残りの部分を参照すれば、本発明の他の特徴及び利点が認識できよう。本発明のさらなる特徴及び利点、並びに本発明の様々な実施態様の構造及び操作が、添付の図面に関して以下に詳細に記載されている。図面において、参照番号は、同一又は機能的に類似した要素を指す。
【発明を実施するための最良の形態】
【００１３】
本発明は、シグモイド又は増殖型曲線を示すデータが、有意な増殖を示しているかを決定するシステム及び方法を提供する。ある態様では、当該データにフィットする一次又は二次多項式曲線が決定され、そしてフィットされた曲線について統計的に有意な値が決定される。当該有意な値が、有意な閾値を超える場合、当該データは、有意又は有効な増殖を示さないとみなされる。当該データが有意又は有効な増殖を示さない場合、当該データセットは、破棄されてもよい。有意な値が有意な閾値を超えることがない場合、当該データは有意又は有効な増殖を示すとみなされる。当該データセットが有効な増殖を示すと決定されたならば、データをさらに処理して、シグモイド又は増殖曲線における転移値、例えばベースライン領域の末端又はＰＣＲ増殖曲線のエルボー値若しくはＣｔ値を決定する。ある態様では、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ(ＬＭ)回帰プロセスにより決定されるパラメーターを有するダブルシグモイド関数を用いて、当該曲線の近似を発見した。パラメーターが決定されると、決定されたパラメーターの１以上を用いて当該曲線を正規化することができる。正規化後、正規化された曲線を処理して、曲線に沿って幾つか又は全ての点で曲線の曲率を決定し、例えば曲率 対 サイクル数を表すデータセット又はプロットをもたらす。最大曲率が生じるサイクル数は、Ｃｔ値に一致する。当該Ｃｔ値は次に返され、そして表示されるか又はそうでない場合さらなる処理に用いられることがある。
【００１４】
本発明の代表的な実施態様では、当該方法は、非限定的にデータセットを入力する入力装置、例えばキーボード、マウスなど；曲線領域の目的の特異的点を表示する表示装置、例えばモニター；当該方法の各ステップを実施するために必須の処理装置、例えばＣＰＵ；モデムなどのネットワークインターフェース、データセットを格納するデータ貯蔵装置、プロセッサー上で作動するコンピューターコードなどを含む慣用されているパーソナルコンピューターシステムを用いることによって、実行され得る。さらに、本方法は、ＰＣＲ装置中で実行されることもある。
【００１５】
図２４は、ソフトウェア資源とハードウェア資源との間の関係を説明する一般的なブロック図を示す。当該システムは、サーモサイクラー装置内に位置する動的ＰＣＲ分析モジュール、並びにコンピューターシステムの一部であるインテリジェンス・モジュールを含むことがある。ここでは、当該データセット(ＰＣＲデータセット)は、ネットワーク接続又は直接接続を介して分析モジュールからインテリジェンス・モジュールへと移されるか、又はその逆である。図２又は９に表示されるように当該方法に従って、プロセッサー上で作動するコンピューターコードにより、当該データセットが処理され、そしてインテリジェンス・モジュールの貯蔵装置上に格納され、そして処理後に分析モジュールの貯蔵装置に戻され、ここで当該改変データは、表示装置上にディスプレーされてもよい。
【００１６】
有効な増殖を示すＰＣＲデータについてのＣｔ決定
ＰＣＲプロセスの点で、増殖又は増幅曲線１０の一例が図１に示される。示されるように、曲線１０は、ラグフェーズ領域１５、及び指数フェーズ領域２５を含む。ラグフェーズ領域１５は、一般的にベースライン又はベースライン領域と呼ばれる。このような曲線１０は、ラグフェーズ領域と指数フェーズ領域を繋ぐ目的の遷移領域２０を含む。領域２０は、一般的に、エルボー又はエルボー領域と呼ばれる。当該エルボー領域は、一般的に、ベースラインの末端、並びに基礎プロセス(underlying process)の増殖又は増幅率の変化を定義する。一般的なＰＣＲ曲線において、エルボー値又はサイクル閾値(Ｃｔ)値と呼ばれる遷移点を同定することは、ＰＣＲプロセスの効率特徴を理解するのに有用である。
【００１７】
同様のシグモイド又は増殖曲線を提供し得る他のプロセスは、細菌プロセス、酵素プロセス、及び結合プロセスを含む。細菌増殖曲線では、例えば、目的の転移点が、ラグフェーズ領域における時間、シータと呼ばれた。本発明に従って分析され得るデータ曲線をもたらす他の具体的プロセスは、標準置換増幅(ＳＤＡ)プロセス、核酸配列に基づく増幅(ＮＡＳＢＡ)プロセス及び遷移媒介性増幅(ＴＭＡ)プロセスを含む。ＳＤＡ及びＮＡＳＢＡプロセスの例及びデータ曲線は、Wang, Sha-Shaら、「Homogeneous Real Time Detection of Single Nucleotide Polymorphisms by Standard Displacement Amplification on the BD ProbeTec ET System」Clin Chem 2003 49 (10): 1599, 及びWeusten, Jos J. A. M.ら、「Principles of Quantitation of viral loads Using Nucleic Acid Sequence-Based Amplification in Combination With Homogeneous Detection Using Molecular Beacons」, Nucleic Acids Research, 2002 30 (6): 26に見出すことができる。こうして、本明細書の残りは、本発明の実施態様及び態様を、ＰＣＲ曲線への適用の点で論じるが、本発明が他のプロセスに関するデータ曲線に適用され得るということが認められるべきである。
【００１８】
図１に示される様に、一般的なＰＣＲ増殖曲線についてのデータは、例えばＰＣＲサイクル数をｘ軸とし、そして増幅されたポリヌクレオチドの指標をｙ軸とした二次元座標系で表すことができる。一般的に、蛍光マーカーの使用が恐らく最も広く使用される標識スキームであるので、図１に示される様に増幅の指標は蛍光強度値である。しかしながら、使用される特定の標識及び/又は検出スキームに応じて、他の指標を使用することも可能であることを理解されたい。増幅されたシグナルの他の有用な指標の例として、発光強度、化学発光強度、生物発光強度、リン光強度、荷電移動、電圧、電流、電力、エネルギー、温度、粘度、光散乱、放射活性強度、反射性、透過率、及び吸光度が上げられる。サイクルの定義には、時間、プロセスサイクル、単位動作サイクル、及び再生サイクルを含めることも可能である。
【００１９】
一般的プロセスの概説
本発明に従って、単一のシグモイド曲線における転移値、例えば動的ＰＣＲ増幅曲線のエルボー値又はＣｔ値などの決定するためのプロセス１００の１の実施態様は、図２を参照して簡潔に記載することができる。ステップ１１０では、曲線を表す実験データセットを受容し、またはそうでない場合取得する。プロットされたＰＣＲデータセットを図１に示し、ここでｙ軸及びｘ軸は、それぞれＰＣＲ曲線についての蛍光強度とサイクル数を表す。ある態様では、当該データセットは、連続的であり、そして軸に沿って均等に配置されたデータを含むべきである。
【００２０】
プロセス１００が、ＰＣＲデータ取得装置、例えばサーモサイクラー内に位置するインテリジェンス・モジュール(例えば、指示を実行するプロセッサー)で実行される場合、データセットは、データが取得されるとリアルタイムでインテリジェンス・モジュールに提供されるか、又はデータセットは記憶ユニット又はバッファーに格納され、そして実験が完了した後にインテリジェンス・モジュールに提供されることもある。同様に、当該データセットは、取得装置に対するネットワーク接続(例えば、ＬＡＮ、ＶＰＮ、イントラネット、インターネットなど)又は直接接続(例えば、ＵＳＢ又は他のダイレクトワイヤー又はワイヤレス接続)を介して、デスクトップコンピュータシステム又は他のコンピューターシステムなどの分離したシステムに提供されるか、又はＣＤ、ＤＶＤ、フロッピー（登録商標）ディスクなどのポータブルメディアで提供されることがある。
【００２１】
ある態様では、データは、座標値の対を有するデータポイント(又は２次元ベクトル)を含む。ＰＣＲデータでは、座標値の対は、一般的に、サイクル数及び蛍光強度値を一般的に表す。ステップ１１０でデータセットを受容するか又は取得した後に、データセットを分析して、ベースライン領域の末端を決定することもできる。
【００２２】
ステップ１２０では、近似曲線が計算される。このステップの間に、１の実施態様では、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ(ＬＭ)回帰プロセス又は他の回帰プロセスにより決定されたパラメーターを有するダブルシグモイド関数を使用して、当該データセットを表す近似曲線を見出した。異常値又はスパイクポイントが曲線フィットの品質に最小の効果しか与えないので、近似は「強健（robust)」であると言われる。以下に論じられる図１３は、受容されたデータセットのプロットと、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ回帰プロセスを用いて本発明に従ったダブルシグモイド関数のパラメーターを決定することにより決定されたデータセットの強健な近似のプロットの例を示す。
【００２３】
ある態様では、データセット中の異常値又はスパイクポイントが取り除かれるか、又はデータセットを処理する前に置換して、ベースライン領域の末端を決定した。スパイクの除去は、データセットをステップ１１０で取得する前又は後で行なわれる。図３は、ＰＣＲ又は他の増殖曲線を表すデータセットにおいてスパイクポイントを同定しそして置換するための処理の流れを示す。スパイクポイントを決定し、そして除去又は置換するためのプロセスのより詳細な記載は、ＵＳ２００７-０１４８６３２に見出すことができる。
【００２４】
ステップ１３０では、ステップ１２０で決定されたパラメーターは、例えばベースラインの傾きを取り除くために、曲線を正規化するために使用される。この様式における正規化は、曲線のベースライン領域の末端又はベースラインストップ位置を決定又は特定する必要なく、Ｃｔ値を決定することを可能にする。ステップ１４０では、当該正規化された曲線を処理して、以下により詳細に記載されるように、Ｃｔ値を決定した。
【００２５】
ＬＭ回帰プロセス
以下に論じられるように、図３のステップ５０２〜５２４は、データセットの曲線を近似し、そしてフィット関数のパラメーターを決定するためのプロセスフローを説明する(ステップ１２０)。これらのパラメーターは、曲線を正規化するのに、例えば本発明の１の実施態様に従って、シグモイド又はＰＣＲ曲線などの増殖型の曲線を表すデータセットのベースライン傾きを変更又は除去するのに使用することができる(ステップ１３０)。当該データセットが、スパイクポイントを除去又は取り除いて、変更済みデータセットを得るように処理された場合、変更されたスパイクを伴わないデータセットは、ステップ５０２〜５２４に従って処理されて、フィット関数のパラメーターを同定することができる。
【００２６】
示される様に１の実施態様では、データセットの強健な曲線近似を計算するために、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ(ＬＭ)法が用いられる。ＬＭ法は、非線形回帰プロセスであり；非線形関数とデータセット間の距離を最小化する反復的な技法である。当該プロセスは、最急降下法とガウス-ニュートン法の組み合わせたもののように動作する(現在の近似が良好にフィットしていない場合には、最急降下法のように動作し(低速であるが、信頼性の高い収束)、現在の近似が正確性を増すに伴って、Ｇａｕｓｓ-Ｎｅｗｔｏｎプロセスのように動作することになる(高速であるが、信頼性の低い収束))。ＬＭ回帰法は、非線形の回帰問題を解決するために広く使用されている。
【００２７】
一般に、ＬＭ回帰法は、様々な入力を必要とし、且つ、出力を提供するアルゴリズムを含んでいる。ある態様においては、入力は、処理対象のデータセット、データのフィットに使用する関数、及び関数のパラメータ又は変数の初期推定値を含んでいる。出力は、関数とデータセット間の距離を最小化する関数のパラメータの組を含んでいる。
【００２８】
ある態様によれば、フィット関数は、次のような形態のダブルシグモイドである：
【数１】
0006009724.tif 000002
【００２９】
この式は、一般的なＰＣＲ曲線及びその他の増殖曲線が取り得る様々な曲線形状にフィットするその柔軟性及び能力に基づいて、フィット関数として選択されたものである。当業者であれば、必要に応じて、その他のフィット関数を使用可能であることを理解するであろう。
【００３０】
ダブルシグモイド式(１)は、(ａ)、(ｂ)、(ｃ)、(ｄ)、(ｅ)、(ｆ)及び(ｇ)という７つのパラメータを具備している。この式は、定数、傾き、及びダブルシグモイドの集まりに分けることが可能である。ダブルシグモイド自体は、２つのシグモイドの乗算(multiplication)である。図４は、ダブルシグモイド式(１)を分解したものを示している。パラメータ(ｄ)、(ｅ)、(ｆ)、及び(ｇ)が２つのシグモイドの形状を決定している。最終的な曲線に対するこれらの影響を示すために、次の単一のシグモイドを検討する：
【数２】
0006009724.tif 000003
【００３１】
式中、パラメータ(ｄ)は、曲線の「鋭さ」を決定しており、パラメータ(ｅ)は、変曲点のｘ値を決定している。図５は、曲線に対するパラメータ(ｄ)の影響と変曲点のｘ値の位置に対するパラメータ(ｅ)の影響を示している。以下の表１は、ダブルシグモイド曲線に対するパラメータの影響を説明している。
【００３２】
【表１】
0006009724.tif 000004
【００３３】
ある態様においては、曲線が非現実的な形状になることを防止するべく、ダブルシグモイド式の「鋭さ」パラメータ(ｄ)及び(ｆ)を制約する必要がある。このため、ある態様においては、ｄ＜-１又はｄ＞１.１或いは、ｆ＜-１又はｆ＞１.１の反復を不成功であると見なしている。その他の態様においては、パラメータ(ｄ)及び(ｆ)に対して異なる制約を使用可能である。
【００３４】
Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔアルゴリズムは、反復的なアルゴリズムであるため、通常、フィットのために関数のパラメータの初期推定値が必要である。初期推定値が良好なものであるほど、近似も良好なものになり、アルゴリズムが局所的な最小値に収束する可能性が低くなる。ダブルシグモイド関数の複雑性とＰＣＲ曲線又はその他の増殖曲線の様々な形状に起因し、すべてのパラメータに１つの初期推定値を付与する方法は、アルゴリズムが頻繁に局所的な最小値に収束することを防止するのに十分ではないであろう。従って、ある態様においては、多数の(例えば、３つ以上)初期パラメータの組を入力し、最良の結果を保持している。ある態様においては、使用する多数のパラメータの組にわたって大部分のパラメータを一定に維持しており、多数のパラメータの組のそれぞれについて異なることができるのは、パラメータ(ｃ)、(ｄ)、及び(ｆ)のみである。図６は、異なるパラメータの組の３つの曲線形状の例を示している。これらの３つのパラメータの組の選択は、ＰＣＲデータを表す曲線の３つの可能な異なる形状を示している。尚、３つを上回る数のパラメータの組を処理して最良の結果を保持することも可能であることを理解されたい。
【００３５】
図３に示されているように、ステップ５１０において、ＬＭ法の初期入力パラメータが同定されている。これらのパラメータは、操作者が入力するか、或いは、算出することも可能である。態様によっては、パラメータは、後述するように、段階５０２、５０４、及び５０６に従って決定又は設定される。
【００３６】
初期パラメータ(ａ)の算出：
パラメータ(ａ)は、ベースラインの高さであり、この値は、初期バラメータのすべての組について同一である。ある態様においては、ステップ５０４において、例えば、蛍光値などの３番目に小さなｙ軸値がデータセットからパラメータａに割り当てられる。この結果、強健な計算が得られる。当然のことながら、その他の態様においては、必要に応じて、最小のｙ軸値や２番目に小さな値などの任意のその他の蛍光値をパラメータａに割り当て可能である。
【００３７】
初期パラメータ(ｂ)の算出：
パラメータ(ｂ)は、ベースラインとプラトーの傾きである。この値は、すべての初期パラメータの組について同一である。全く傾きが存在しないことが理想的であるため、ある態様においては、ステップ５０２において、０.０１の静的な値が(ｂ)に割り当てられる。その他の態様においては、パラメータ(ｂ)には、例えば、０〜約０.５の範囲の値などの異なる値を割り当て可能である。
【００３８】
初期パラメータ(ｃ)の算出：
パラメータ(ｃ)は、「プラトーの高さ」 - 「ベースラインの高さ」を表しており、これは絶対蛍光強度、つまりＡＦＩを指す。１の態様では、パラメーターの第一セットでは、ｃ＝ＡＦＩ+２であり、一方最後の２つのパラメーターについては、ｃ＝ＡＦＩである。これは、図６に示されており、ここでパラメーターの最後の２つのセットでは、ｃ＝ＡＦＩである。パラメーターの第一セットでは、ｃ＝ＡＦＩ+２である。この変化は、第一セットのパラメーターによりモデル化された曲線の形状に起因しており、この曲線形状は、プラトーを具備していない。
【００３９】
パラメータ(ｄ)及び(ｆ)の算出：
パラメータ(ｄ)及び(ｆ)は、２つのシグモイドの鋭さを定義している。曲線に基づいてこれらのパラメータの近似を付与する方法は存在しないため、ある態様においては、３つの静的な代表値をステップ５０２において使用する。その他の静的な又は非静的な値をパラメータ(ｄ)及び／又は(ｆ)に使用することも可能であることを理解されたい。これらのペアは、遭遇するＰＣＲ曲線の最も共通的な形状をモデル化している。次の表２は、図６に示されている異なるパラメータの異なるセットについて(ｄ)及び(ｆ)の値を示している。
【００４０】
【表２】
0006009724.tif 000005
【００４１】
パラメータ(ｅ)及び(ｇ)の算出：
ステップ５０６において、パラメータ(ｅ)及び(ｇ)を決定する。パラメータ(ｅ)及び(ｇ)は、２つのシグモイドの変曲点を定義している。１の態様においては、これらは、いずれも、すべての初期パラメータセットにおいて同一の値を取っている。パラメータ(ｅ)及び(ｇ)は、同一の又は異なる値を具備可能である。近似を見出すべく、ある態様においては、例えば、蛍光などの強度の平均値を上回る最初の点のｘ値を使用している(これは、スパイクではない)。この態様に従って(ｅ)及び(ｇ)の値を決定するプロセスが図７に示されており、これについては後述する。パラメータ(ｅ)及び(ｇ)並びにその他のパラメーターの値を決定するプロセスの更に詳細な説明は、ＵＳ２００７-０１４８６３２に見出すことができる。
【００４２】
図７を参照すれば、まず、曲線(例えば、蛍光強度)の平均値が決定される。次に、平均値を上回る最初のデータ点を同定する。次いで、
(ａ)その点が、曲線の開始点の近傍、例えば最初の５つのサイクル内に位置していないかどうか、
(ｂ)その点が、曲線の終了点の近傍、例えば最後の５つのサイクル内に位置していないかどうか、及び
(ｃ)その点の周辺(例えば、その周囲の２点の範囲における)導関数が符号の変化を示していないかどうか(示している場合には、その点は、スパイクである可能性が高く、従って、取り除く必要がある)、
が決定される。
【００４３】
次の表３は、１の態様に従って図６において使用されている初期パラメータの例を示す。
【００４４】
【表３】
0006009724.tif 000006
【００４５】
図３を再度参照すれば、ステップ５１０においてすべてのパラメータの設定が完了したら、入力されたデータセット、関数、及びパラメータを使用して、ＬＭプロセス５２０を実行する。伝統的に、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ法は、非線形の最小二乗問題を解決するべく使用されている。従来のＬＭ法においては、近似曲線とデータセット間の誤差の二乗の合計として定義された距離の尺度を算出している。しかしながら、二乗の合計を最小化する場合には、異常値の距離が非スパイクデータ点の距離よりも大きいため、異常値に対して大きな重みが付与され、この結果、しばしば、不適切な曲線又は望ましくない曲線が得られることになる。従って、本発明の１の態様によれば、絶対誤差の合計を極小化することにより、近似とデータセット間の距離を演算している(この場合には、異常値に対して、それほど大きな重みが付与されない)。この態様においては、近似とデータ間の距離は、次式によって与えられる：
【数３】
0006009724.tif 000007
【００４６】
前述のように、１の態様においては、多数(例えば、３つの)の初期パラメータのセットの各々を入力及び処理し、そしてステップ５２２及び５２４に示されるように最良の結果が保持される。ここで、最良の結果とは、式(３)において、最も小さい又は最小の距離を提供するパラメータである。１の態様においては、パラメータの多くは、パラメーターの複数のセットにわたって一定に保持しいる。(ｃ)、(ｄ)、及び(ｆ)のみが、パラメータの各セットについて異なってもよい。任意の数の初期パラメータセットが使用可能であることを理解されたい。
【００４７】
図８は、本発明によるパラメータのセットについてのＬＭプロセス５２０のプロセスフローを示している。前述のように、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ法は、最急降下法又はＧａｕｓｓ-Ｎｅｗｔｏｎ法のように動作可能である。この動作は、減衰係数λによって左右される。λが大きいほど、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔアルゴリズムは、最急降下プロセスのように動作することになる。一方、λが小さいほど、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔアルゴリズムは、Ｇａｕｓｓ-Ｎｅｗｔｏｎプロセスのように動作することになる。１の態様においては、０.００１でλが開始している。λが、約０.０００００１〜約１.０などのその他の値において開始することも可能であることを理解されたい。
【００４８】
前述のように、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ法は反復的な技法である。図８に示されているように、１の態様によれば、各反復において以下の内容が実行される。
１.先行する近似のヘッセ行列(Ｈ)を算出する。
２.先行する近似の転置ヤコビ行列(ＪT)を算出する。
３.先行する近似の距離ベクトル(ｄ)を算出する。
４.現在の減衰係数λによりヘッセ行列の対角を拡大させる：
【数４】
0006009724.tif 000008
５.拡大した式を解く：
【数５】
0006009724.tif 000009
６.拡大した式の解ｘを関数のパラメータに加算する。
７.新しい近似と曲線間の距離を算出する。
８.この新しいパラメータセットとの距離が、以前のパラメータセットとの距離を下回っている場合には：
・その反復を成功であると見なし、
・新しいパラメータの組を保持又は保存し、
・例えば、係数１０だけ減衰係数λを低減し、
この新しいパラメータセットとの距離が、以前のパラメータセットとの距離を上回っている場合には、
・その反復を不成功であると見なし、
・新しいパラメータの組を破棄し、
・例えば、係数１０だけ、減衰係数λを増大させる。
【００４９】
１の態様においては、図８のＬＭプロセスは、次の基準の中の１つが実現されるまで反復される：
１.指定の反復回数Ｎだけ既に反復されている。この第１の基準により、アルゴリズムが無限に反復されることを防止している。例えば、図１０に示されているある態様においては、既定の反復値Ｎは１００である。収束が可能な場合には、１００回の反復は、アルゴリズムが収束するのに十分なものであろう。一般に、Ｎは、１０未満〜１００超の範囲をとることができる。
２.２つの成功した反復の間の距離の差が、閾値(例えば、０.０００１)を下回っている。差が非常に小さくなっている場合には、既に望ましい精度が実現された状態にあり、解は、それ以上、大幅に良好なものにはならないため、反復を継続するのは無意味である。
３.減衰係数λが、指定の値を上回っている(例えば、１０２０を上回っている)。λが非常に大きくなった場合には、アルゴリズムは、現在の解よりも良好には収束せず、従って、反復を継続するのは無意味である。一般に、指定値は、１０２０よりも十分に小さい又は大きいものであってもよい。
【００５０】
正規化
パラメーターを決定した後に、１の実施態様では、曲線は、１以上の決定されたパラメーターを用いて正規化される(ステップ１３０)。例えば、１の態様では、曲線は、正規化されるか、又は曲線の線形増殖部分を減算することにより、ゼロベースラインの傾きを有するように調節される。数学的には、これは以下のように示される：
【数６】
0006009724.tif 000010
式中、ｄａｔａＮｅｗ(ＢＬＳ)は、ベースラインを引いた後に正規化されたシグナル、例えば線形増殖又はベースラインの傾きが引き算されるか又は除かれたデータセット(データ)である。パラメーターａ及びｂの値は、曲線に回帰させるためにＬＭ式を用いることにより決定された値であり、そしてｘはサイクル数である。こうして、ｘ軸にそった各データ値について、定数ａ＋傾きｂ×ｘ値が、当該データから引かれて、ゼロベースラインの傾きを具備したデータ曲線をもたらす。ある態様では、スパイク点は、ＬＭ回帰プロセスをデータセットに適用して正規化パラメーターを決定する前に、データセットから取り除かれる。
【００５１】
別の態様では、曲線は、以下の式に従ったゼロ傾きとなるように正規化又は調節されてもよい：
【数７】
0006009724.tif 000011
式中、ｄａｔａＮＥＷ(ＢＬＳＤ)は、ベースラインを引いた後の正規化されたシグナルを割ったものであり、例えば、線形増殖又はベースラインの傾きを引き算するか又は取り除き、そして結果をａで割ったデータセット(データ)である。パラメーターａ及びｂの値は、ＬＭ式を用いて曲線を回帰させることにより決定される値であり、そしてｘはサイクル数である。こうして、ｘ軸にそった全てのデータ値について、定数ａ及び傾きｂ×ｘ値が当該データから引き算され、そして結果をパラメーターａの値で割って、ゼロベースライン傾きを具備するデータ曲線をもたらした。１の態様では、式(７ａ)は、パラメーター「ａ」≧１の場合妥当であり；パラメーター「ａ」＜１の場合、以下の式を用いる：
【数８】
0006009724.tif 000012
【００５２】
ある態様では、データセットにＬＭ回帰プロセスを適用して、正規化パラメーターを決定する前に、スパイクポイントはデータから取り除かれる。
【００５３】
さらに別の態様では、曲線は、以下の式：
【数９】
0006009724.tif 000013
にしたがって正規化又は調節することもできる。式中、ｄａｔａＮＥＷ(ＢＬＤ)は、ベースラインで割り算した後の正規化シグナル、例えばパラメーターａで割られたデータセット(データ)である。パラメーターａ及びｂの値は、ＬＭ式を用いて曲線へと回帰させることにより決定され値であり、そしてｘはサイクル数である。１の態様では、式(８ａ)は、パラメーター「ａ」≧１の場合妥当であり、パラメーター「ａ」＜１の場合以下の式：
【数１０】
0006009724.tif 000014
を使用する。
【００５４】
ある態様では、ＬＭ回帰プロセスをデータセットに適用して、正規化パラメーターを決定する前に、データセットからスパイクポイントを取り除く。
【００５５】
さらにべつの態様では、曲線は、以下の式：
【数１１】
0006009724.tif 000015
に従って正規化又は調節することもできる。式中、ｄａｔａＮｅｗ(ＰＧＴ)は、ベースラインを減算して割算したのちの正規化シグナル、例えば線形増殖又はベースラインの傾きを引き算するか又は取り除き、そして結果をｃで割ったデータセット(データ)である。パラメーターａ、ｂ及びｃの値は、ＬＭ式を用いて曲線へと回帰させることにより決定した値であり、そしてｘはサイクル数である。こうして、ｘ軸に沿った各データ値について、定数ａと傾きｂ×ｘ値がデータから引き算され、そして結果がパラメータｃの値により割られて、ゼロベースラインの傾きを具備するデータ曲線を与えた。１の態様では、式(９ａ)はパラメータ「ｃ」≧１の場合妥当であり；パラメータ「ｃ」＜１及び「ｃ」≧０の場合、以下の式：
【数１２】
0006009724.tif 000016
を用いる。
【００５６】
ある態様では、ＬＭ回帰プロセスをデータセットに適用して正規化パラメータを決定する前に、データセットからスパイクポイントが取り除かれる。
【００５７】
Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ又は他の回帰プロセスにより決定されるパラメーターを用いてベースラインを正規化及び/又は変更するために、他の正規化式を使用することができるということが、当業者に認められよう。
【００５８】
曲率決定
式(６)、(７)、(８)、又は(９)のうちの１つ、或いは他の正規化式を用いて、曲線を正規化した後に、Ｃｔ値を決定することができる。１の実施態様では、曲率決定プロセス又は方法を、図９を参照して説明される正規化曲線に適用する。図９は、動的ＰＣＲ曲線におけるエルボー値又はＣｔ値を決定するためのプロセスフローを示す。ステップ９１０では、当該データセットを取得する。サーモサイクラーなどのＰＣＲデータ取得装置内に配置されたインテリジェンス・モジュール(例えば、指令を実行するプロセッサー)中で決定プロセスが実行する場合、当該データセットは、データを取得するとリアルタイムでインテリジェンス・モジュールにデータセットが提供されるか、又はデータセットはメモリーユニット又はバッファーに格納されて、そして実験が完了した後に当該モジュールへと提供されてもよい。同様に、データセットは、取得装置に対するネットワーク接続(例えば、ＬＡＮ、ＶＰＮ、イントラネット、インターネットなど)又は直接接続(例えば、ＵＳＢやその他の直接的な有線又は無線接続)を介して、デスクトップコンピュータシステムなどの別個のシステムに提供することも可能であり、或いは、ＣＤ、ＤＶＤ、フロッピー（登録商標）ディスク、又はこれらに類似したものなどの携帯型の媒体上において提供することも可能であろう。
【００５９】
データセットを受信又は取得した後に、ステップ９２０において、曲線への近似を決定する。このステップの間、１の実施態様では、Ｌｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔプロセスにより決定されるパラメーターを有するダブルシグモイド関数を用いて、データセットを表す曲線の近似を発見した。さらに、スパイクポイントは、図３を参照して記載されるように、ステップ９２０の前に、データセットから取り除くことができる。例えば、ステップ９１０において取得されるデータセットは、事前にスパイクが除去されているデータセットであってもよい。ステップ９３０において、曲線を正規化する。ある態様では、上記式(６)、(７)、(８)又は(９)のうちの一つを用いて曲線を正規化する。例えば、上記式(６)に記載されるようにベースライン傾きを引き算するためにステップ９２０で決定されるダブルシグモイド式のパラメーターを用いて、ベースラインをゼロスロープに設定することができる。ステップ９４０では、プロセスを正規化曲線へと適用して、正規化曲線に沿った点で曲率を決定する。曲率 対 サイクル数のプロットを返すことができ、及び/又は表示することができる。最大曲率の点は、エルボー又はＣｔ値に対応する。ステップ９５０では、分析を行なうシステムへと、又は当該分析を要求する別のシステムへと結果を返す。ステップ９６０では、Ｃｔ値が表示される。全体のデータセット又は曲線近似などの追加データを表示することもできる。図９の分析を行なったシステムとつながれたモニター、スクリーン又はプリンターなどの表示装置に、グラフ表示が映されるか、又はデータは、表示装置上に映すための別のシステムへと提供されることもある。
【００６０】
１の実施態様では、この曲線についてＣｔ値を得るために、最大曲率が決定される。１の態様では、曲率は正規化曲線上の幾つか又は全ての点について決定される。曲率対サイクル数のプロットが表示することができる。曲線の曲率は、以下の式：
【数１３】
0006009724.tif 000017
により与えられる。
【００６１】
円の半径ａを考慮すると、以下の式：
【数１４】
0006009724.tif 000018
により与えられる。
【００６２】
式(１１)の曲率は、ｋａｐｐａ(ｘ)＝-(１／ａ)である。こうして、曲率の半径は、曲率の逆数の負の値である。円の半径ａが定数であるので、その曲率は、-(１／ ａ)で与えられる。図１０ａのＰＣＲデータセットのフィットについての曲率のプロットである図１０ｂを検討する。Ｃｔ値は、最大曲率の位置で生じるとみなすことができ、これはサイクル数Ｃｔ＝２１.８４で生じる。好ましくはこのＣｔ値は、図１０ａに示されるＰＣＲ増殖曲線と比較される。
【００６３】
(２１.８４のＣｔ値に対応する)最大曲率における曲率の半径は：半径＝１／０.２８１８＝３.５５サイクルである。図１０ａのＰＣＲ増殖曲線に重ねあわされた半径の円が図１１に示される。図１１が、最大曲率に対応する半径の円は、ＰＣＲ曲線の増殖領域の開示時点で重ねあわすことができる一方、曲線への接線が残る。曲率の小さい(最大)半径を有する曲線は、急勾配の増殖曲線を有する一方、曲率の大きい(最大)半径を有する曲線は、緩やかな増殖曲線を有することもある。曲率の半径がかなり大きいならば、明確なシグナルがない曲線、例えば僅かな増殖又は妥当でない増殖を指し示しうる。１の実施態様では、しかしながら、以下に詳細に記載されるように、データセットが有意又は妥当な増殖を示しているかを決定するために、増殖妥当性試験が提供される。データセットが統計的に有意な増殖を有すると発見されるならば、曲率分析アルゴリズムは、Ｃｔ値を決定するために適用することができる。そうでない場合、データセットが破棄されてもよく、及び／又は妥当でない増殖の表示が返される。
【００６４】
計算に必要とされる式(１)のダブルシグモイドの一次導関数及び二次導関数は以下に示される：
一次導関数：
【数１５】
0006009724.tif 000019
式(１３)二次導関数：
【数１６】
0006009724.tif 000020
【実施例】
【００６５】
図１２ａは、増殖曲線について生データの例を示す。ダブルシグモイド／ＬＭ法を図１２ｂに示される生データプロットに適用することは、以下の表４に示される式(１)の７個のパラメーターの値を与える。
【表４】
0006009724.tif 000021
【００６６】
図１２に示されるデータにフィットされたダブルシグモイドは、図１３に示され、データポイントのかなり正確な評価を示す。次に、これらのデータを式(６)(ベースライン引き算)に従って正規化して、図１４に示されるグラフをもたらす。図１４に示される実線は、式(６)に従って正規化されたデータセットに対する式(Ｉ)のダブルシグモイド／ＬＭの適用を示す。図１５は、図１４の正規化された曲線についての、曲率 対 サイクル数のプロットを示す。曲率の最大値は、曲率０.１３７８、サイクル数３４.４２で生じる。こうして、最大曲率でのサイクル数に基づいてＣｔ＝３４.４２であり、そして曲率の半径＝１／０.１３７８＝７.２５である。当該曲率の半径を有する円と、正規化されたデータセットの重ね合わせを図１６に示す。
【００６７】
「緩やかに増加する」データセットの例を図１７に示す。このデータセットに対するダブルシグモイドのフィットは及びベースラインの減算(式６)を用いた正規化は、図１８に示されるフィットの結果を与えた。対応する曲率プロットを図１９に示す。サイクル数２５.９０、曲率＝０.００１０９２７４(曲率半径＝９１５に対応する)で最大曲率が生じる。この曲率の大きな半径は、緩やかな増殖データセットであるということを伝える。
【００６８】
別の例では、図２０に示されるＰＣＲ増殖曲線のセットが考慮される。現在の方法(閾値)を用いて得られるＣｔ値  対  ベースラインを引き算し割算した(ＢＬＳＤ-式７)後の曲率法を用いて取得されたＣｔ値の比較が、以下の表５に示される。
【００６９】
【表５】
0006009724.tif 000022
【００７０】
表５は、Ｃｔ値を計算する曲率法(この場合、ＢＬＳＤで正規化した後である)は、現剤の閾値法よりも小さいＣｖ(バリエーション係数)を与える。さらに、曲率法を用いて計算された曲率半径(ＲＯＣ)は、曲線が線形曲線であるか、又は真の増殖曲線であるかを示唆する単純な方法を提供する。
【００７１】
増殖妥当性試験
曲率が存在するためには、ＰＣＲシグナルは、高次多項式(一般的には、上記のように７次以上（a power of 7 or higher))により表すことができなければならない。そうでなければ、シグナルは、以下の式：
【数１７】
0006009724.tif 000023
の一次又は二次の多項式により表すことができ(例えば、Ｒ2≧０.９０)、次にこのようなシグナルについて曲率が、１の態様では、以下のように決定される：
(１)式(１４)からベースラインの引き算を行い、以下の式(１５)：
【数１８】
0006009724.tif 000024
をもたらし：
(２)式(１５)についての曲率を、以下の式(１６)：
【数１９】
0006009724.tif 000025
として与え；
(３)この曲率関数(式１６)は、ｘ＝０で最大値を有し、その結果二次関数への優れた曲線フィットを有するＰＣＲシグナルについての規定のエルボー値が存在しないということを示唆する。こうして、１の実施態様では、データセットが、統計的に有意なマージンで１次又は２次多項式にフィットするならば、当該データセットは有意な増殖を欠いていると決定される。
【００７２】
１の実施態様では、増殖プロセスについてのデータセットは、当該データが有意な増殖を示すかを決定するために処理される。最初は、当該データにフィットする１次又は２次多項式曲線が計算され(例えば、式(１４))、そして次に統計的に有意な値が曲線フィットについて測定される。ある態様では、統計的に有意な値はＲ2値である。統計的に有意な値が、閾値を超えない場合、当該データセットは、統計的に有意であるか又は有効な増殖を示すことが判断され、そして当該データセットはさらに処理されて、例えばＣｔ値を決定する。１の態様では、Ｒ2閾値は約０.９０であり；Ｒ2が０.９０を超えている場合、当該データセットは、有効ではない(例えば、有意な増殖ではない)と判断される。別の態様では、Ｒ2閾値は０.９９である。Ｒ2閾値が、約０.９０〜０.９９で設定されるか、又は閾値が０.９９を超えるか、又は０.９０未満であることがあるということが認められるべきである。統計的に有意な値は、閾値を超える場合、データセットは有意ではなく、又は有効ではない増殖を示すと判断される。当該データセットが有意な増殖を有さないということを示すメッセージが返されてもよいし、及び／又は当該データセットが破棄されてもよい。
【００７３】
図２１は、標的を含まないリアルタイムＰＣＲデータシグナルを示し、そしてこれは許容される範囲でベースラインインターセプト、傾き及びＡＦＩ値を具備する。式（１０)、(１２)、及び(１３)の曲率アルゴリズムは、当該Ｃｔ値が１２.９４であることを指し、そして(最大)曲率半径(ＲＯＣ)が４８１であるということを指す。増殖妥当性試験が適用された場合、当該データは、不十分な増殖又は不十分な曲率を有する決定され、当該シグナルが、統計的に意味ある値が閾値を超える、例えばＲ2＞０.９０である１次又は２次多項式にフィットすることを意味する。
【００７４】
図２２は、４８１のＲＯＣを有する別のリアルタイムＰＣＲデータシグナルを示す。Ｒ2値は、閾値、例えば０.９９よりもずっと小さく、当該プロセスはＣｔ値の計算を続ける。式(１０)、(１２)、及び(１３)の曲率アルゴリズムは、曲率の最大半径、及びＣｔ値が、サイクル３８.７で生じることを示す。図２１を図２２と比較すると、ＲＯＣ値のみについての知識では、曲線が有効な増殖を示すかを同定するには不十分であるということが明らかである。ここで、両方のシグナルは、同じ最大ＲＯＣを有しているが、一方は有効な増殖を有するが、もう一方は有効な増殖を有さない。
【００７５】
図２３は、別のリアルタイムＰＣＲシグナルを示す。ＲＯＣアルゴリズムを適用して、Ｃｔ値を決定することにより、サイクル３０.３にて(最大)ＲＯＣ７１でＣｔ値を与える。増殖妥当性試験を適用することにより、有意でなく又は有効ではない増殖であるということが示される。こうして、このかなり低い(最大)ＲＯＣでシグナルは妥当ではなく、低い(最大)ＲＯＣがそれ自身、曲線が有効でないと宣言するために不十分であるということを示す。
【００７６】
曲線フィット及び曲率決定プロセスを含む増殖妥当性試験及びＣｔ決定プロセスが、コンピューターシステムのプロセッサ上で走らされるコンピューターコードで実行されることがある。当該コードは、増殖妥当性Ｃｔ決定プロセスのさまざまな態様及びステップを実行するためにプロセッサーを制御する指令を含む。コードは典型的にハードディスク、ＲＡＭ、又はポータブルメディア、例えばＣＤ、ＤＶＤなどに一般的に格納されている。同様に、当該プロセスは、プロセッサーにつながれたメモリーユニットに格納された指令を実行するプロセッサーを含むサーモサイクラーなどのＰＣＲ装置中で実行される。このような指令を含むコードは、ネットワーク接続又はコードソースへの直接接続を介して、又は周知のようにポータブルメディアを用いてＰＣＲ装置のメモリーユニットにダウンロードされてもよい。
【００７７】
本発明のエルボー決定プロセスが、Ｃ、Ｃ＋＋、Ｃ＃、Ｆｏｒｔｒａｎ、ＶｉｓｕａｌＢａｓｉｃなどの様々なプログラミング言語、並びにデータの可視化及び分析に有用なプレパッケージングルーチン、関数、及び分析を提供する数学などの適用を用いてコードすることができる。後者の別の例は、ＭＡＴＬＡＢ（登録商標）である。
【００７８】
本発明が、例示の方法で記載され、そして具体的実施態様の点で記載されてきたが、本発明が開示された実施態様に限定されるものではないということを理解されたい。反対に、当業者に明らかであるさまざまな変更及び同等のアレンジをカバーすることが意図される。その結果、添付の特許請求の範囲は、このような変更及び同等のアレンジのすべてを包含するように、最も広い解釈に従うべきである。
【図面の簡単な説明】
【００７９】
【図１】図１は、蛍光強度 対 サイクル数としてプロットされた一般的なＰＣＲ増殖曲線の例を指す。
【図２】図２は、増殖曲線のベースライン領域の末端、又はＰＣＲ曲線のＣｔ値を決定するためのプロセスフローを示す。
【図３ａ】図３は、本発明の１の実施態様に従ったスパイク同定及び置換プロセスの詳細なプロセスフローを示す。
【図３ｂ】図３は、本発明の１の実施態様に従ったスパイク同定及び置換プロセスの詳細なプロセスフローを示す。
【図３ｃ】図３は、本発明の１の実施態様に従ったスパイク同定及び置換プロセスの詳細なプロセスフローを示す。
【図３ｄ】図３は、本発明の１の実施態様に従ったスパイク同定及び置換プロセスの詳細なプロセスフローを示す。
【図３ｅ】図３は、本発明の１の実施態様に従ったスパイク同定及び置換プロセスの詳細なプロセスフローを示す。
【図３ｆ】図３は、本発明の１の実施態様に従ったスパイク同定及び置換プロセスの詳細なプロセスフローを示す。
【図４】図４は、パラメーター(ａ)-(ｇ)を含むダブルシグモイド方程式の分離を示す。
【図５】図５は、曲線についてのパラメーター(ｄ)の影響、及び位置(ｅ)、変曲点のｘ値を示す。図５の全ての曲線は、パラメーター(ｄ)を除く同じパラメーター値を有する。
【図６】図６は、異なるパラメーターセットについての３個の曲線形状の例を示す。
【図７ａ】図７は、１の態様に従ったダブルシグモイド方程式パラメーター(ｅ)及び(ｇ)の値を決定するプロセスを示す。
【図７ｂ】図７は、１の態様に従ったダブルシグモイド方程式パラメーター(ｅ)及び(ｇ)の値を決定するプロセスを示す。
【図８ａ】図８は、パラメーターの初期セットについてのＬｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ回帰プロセスのプロセスフローを示す。
【図８ｂ】図８は、パラメーターの初期セットについてのＬｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ回帰プロセスのプロセスフローを示す。
【図８ｃ】図８は、パラメーターの初期セットについてのＬｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ回帰プロセスのプロセスフローを示す。
【図８ｄ】図８は、パラメーターの初期セットについてのＬｅｖｅｎｂｅｒｇ-Ｍａｒｑｕａｒｄｔ回帰プロセスのプロセスフローを示す。
【図９】図９は、１の実施態様に従ったＰＣＲプロセスについてのエルボー値を決定するためのより詳細なプロセスフローを示す。
【図１０】図１０ａは、ダブルシグモイドを用いて実験データにフィットさせた典型的な増殖曲線を示し、図１０ｂは、図１０ａのダブルシグモイド曲線の曲率のプロットを示す。
【図１１】図１１は、図１０ａの増殖曲線において重ねあわされた最大曲率点で接する円を示す。
【図１２ａ】図１２ａは、増殖曲線について記載されるデータセットの例を示す。
【図１２ｂ】図１２ｂは、図１２ａのデータセットのプロットを示す。
【図１３】図１３は、図１２のデータセットにフィットするダブルシグモイドを示す。
【図１４】図１４は、方程式(６)のベースライン・サブトラクション法を用いて正規化したのちの、図１２(図１３)のデータセット(及びダブルシグモイドフィット)を示す。
【図１５】図１５は、図１４の正規化データセットについての曲率 対 サイクル数のプロットを示す。
【図１６】図１６は、曲率の最大半径を有する円と、図１５の正規化されたデータセットの重ね合わせを示す。
【図１７】図１７は、「ゆっくり-増殖(slow glower)」データセットの例を示す。
【図１８】図１８は、図１７のデータセット、並びに方程式(６)のベースラインサブトラクション法を用いて正規化した後のダブルシグモイドフィットを示す。
【図１９】図１９は、図１８の正規化データセットの曲率 対 サイクル数のプロットを示す。
【図２０】図２０は、複製実行及び陰性サンプルを含む、ＰＣＲ増殖曲線のセットのプロットを示す。
【図２１】図２１は、標的を含まないリアルタイムＰＣＲデータシグナルを示し、そしてこれは許容される範囲で、ベースラインインターセプト、傾き、及びＡＦＩ値を有する。
【図２２】図２２は、図２１のシグナルと同じ(最大)半径の曲率を有するリアルタイムＰＣＲデータシグナルを示す。
【図２３】図２３は、小さい(最大)半径の曲率を有するリアルタイムＰＣＲシグナルを示す。
【図２４】図２４は、ソフトウェアとハードウェア資源との間の関係を説明する一般的なブロック図を示す。
【図１】
0006009724.tif 000029
【図２】
0006009724.tif 000030
【図３ａ】
0006009724.tif 000031
【図３ｂ】
0006009724.tif 000032
【図３ｃ】
0006009724.tif 000033
【図３ｄ】
0006009724.tif 000034
【図３ｅ】
0006009724.tif 000035
【図３ｆ】
0006009724.tif 000036
【図４】
0006009724.tif 000037
【図５】
0006009724.tif 000038
【図６】
0006009724.tif 000039
【図７ａ】
0006009724.tif 000040
【図７ｂ】
0006009724.tif 000041
【図８ａ】
0006009724.tif 000042
【図８ｂ】
0006009724.tif 000043
【図８ｃ】
0006009724.tif 000044
【図８ｄ】
0006009724.tif 000045
【図９】
0006009724.tif 000046
【図１０】
0006009724.tif 000047
【図１１】
0006009724.tif 000048
【図１２ａ】
0006009724.tif 000049
【図１２ｂ】
0006009724.tif 000050
【図１３】
0006009724.tif 000051
【図１４】
0006009724.tif 000052
【図１５】
0006009724.tif 000053
【図１６】
0006009724.tif 000054
【図１７】
0006009724.tif 000055
【図１８】
0006009724.tif 000056
【図１９】
0006009724.tif 000057
【図２０】
0006009724.tif 000058
【図２１】
0006009724.tif 000059
【図２２】
0006009724.tif 000060
【図２３】
0006009724.tif 000061
【図２４】
0006009724.tif 000062