JP2005504543A

JP2005504543A - 定量的ｐｃｒのための適応ベースラインアルゴリズム

Info

Publication number: JP2005504543A
Application number: JP2003533071A
Authority: JP
Inventors: ロジャーエイチ．テイラー、; アランアール．オール、
Original assignee: Stratagene California
Current assignee: Stratagene California
Priority date: 2001-10-02
Filing date: 2002-10-02
Publication date: 2005-02-17
Also published as: WO2003029924A2; US20030148332A1; WO2003029924A3; US7363168B2; EP1438430A2; EP1438430A4; CA2462823A1; AU2002327800B2

Abstract

本発明は、リアルタイムＰＣＲ増幅中のチューブ間およびサイクル間の変動を低減するために開発されたベースライン減算アルゴリズムに関する。特に、本発明は、増幅された核酸物の最初の信頼性検出に対する閾値サイクルを決定するためのアルゴリズムに関する。本発明は、該アルゴリズムを含むコンピュータプログラム、ならびに該アルゴリズムを実行する方法および／または装置にも関する。

Description

【技術分野】
【０００１】
本発明は、ベースライン減算アルゴリズム、該アルゴリズムを用いた方法、および該アルゴリズムを実行するための装置およびコンピュータプログラムに関する。
【背景技術】
【０００２】
ポリメラーゼ連鎖反応（ＰＣＲ）は、標的核酸配列をたとえ１コピーであっても特異的に増幅および検出するための、今日広く研究室で用いられている有力な技術である。出発標的配列の量と、任意の所定のサイクルにおけるＰＣＲ産物の量の間には定量的な関係があることから、ＰＣＲは核酸配列の定量化のためにも用いられる。
【０００３】
終点ＰＣＲは、核酸鋳型を増幅する用途で広く用いられている。終点ＰＣＲにおいて、ＰＣＲ反応の開始時に鋳型を加え、反応は、通常は２０〜５０サイクルの複数のサイクルで実施する。検出および／または定量されるのは増幅反応の最終産物である。これに対し、リアルタイム定量的ＰＣＲ（ＱＰＣＲ）は、ＰＣＲ増幅の進行をその最中にモニタする。リアルタイムＱＰＣＲ技術において、（一般的に蛍光性の）シグナルは、該シグナル発生中にモニタされる。所望のレベルの蛍光を得るために必要なサイクルの数は、ＰＣＲプロセス開始時における増幅可能な標的の濃度によって直接的に変化することから、蛍光強度の測定値によって試料中の標的ＤＮＡの量の測定値を与えることができる。蛍光強度は、各ＰＣＲサイクルのアニーリング／伸長期間の間に検出され、この検出の出力はプロセッサに送られて格納およびデータ操作される。終点ＰＣＲの精度はＱＰＣＲよりも劣るが、これは、測定がＰＣＲにおいてより遅くに行われるために、より多くの変数が結果に影響を与える可能性があるためである。例えば、反応成分が消耗されるにつれ、増幅は減少する。終点におけるシグナルレベルの違いは、出発標的よりも制限試薬の僅かな差異に起因しているかも知れない。この制限はＰＣＲ反応の早い段階の方が小さい。
【０００４】
増幅中に得られるデータを、プロセッサによって正規化する。このプロセッサは、バックグラウンドシグナル（標的核酸の非存在下におけるＰＣＲチューブ内の予想シグナル）のベースラインを同定し、測定されたシグナルからバックグラウンドシグナルを除去する。ベースライン減算を行った結果、試料中の標的核酸の量をより正確に反映するシグナル強度の測定値が得られる。ベースライン減算では一般に、増幅前のサイクル中のチューブ内で観察されるバックグラウンドシグナルをベースラインとして設定する。増幅前に行うサイクルの範囲は、ユーザによって設定される終点（例えば開始および終了サイクル）によって定義される。開始および終了サイクルを選択するための典型的な指針は、最初の数サイクルにおける典型的な変動が減じた後に開始サイクルを選択し、どのチューブでも増幅が起こる前に終了サイクルを選択することを規定している。典型的には、これらのサイクルから得られるシグナルを、少なくとも平均二乗アルゴリズムを用いてラインに適合する。最良の適合線を用いて、ＱＰＣＲ実験の全てのサイクルに対するバックグラウンドを予想し、これを、（例えば、米国特許第５，９２８，９０７号に記載のようにして）評価中の各試料より生じたデータから差し引く。
【０００５】
最良適合線算出のための範囲を設定するこの装置には次のようないくつかの制限がある。１）出発量が異なる試料を用いる場合、最後の終点に対する最適な選択が、試料によって異なる。２）ユーザによって（あるいは同じユーザでも時によって）、異なる終点を選択する。３）異なる終点で解析された同一の実験が異なる結果を与えることになる。
【０００６】
ベースライン減算のばらつきは、未知の出発コピーの最終決定においてより大きな誤差を生む原因となる。
【発明の開示】
【発明が解決しようとする課題】
【０００７】
当該技術分野においては、定量的ＰＣＲ増幅中に起こるチューブ間の変動の効果を低減する、適応ベースライン減算が必要とされている。
【課題を解決するための手段】
【０００８】
１つの態様において、本発明は、ＰＣＲ反応のサイクルを循環する試料ブロック中の複数の試料ウェルに対する適応ベースラインの算出法を提供する。本方法は、各ウェルにおいて観測される実光学シグナルの強度を、該ウェルに対するサイクル数の関数としてプロットして、第１のプロットを得る工程と、該ウェルに対する開始サイクルと終了サイクルとを決定する工程と、開始サイクルと終了サイクルの間の両サイクルを含めた点を直線に適合することにより、該点と直線間の最良適合を決定する工程とを含む。しかしながら、１つの態様においては、最良適合を決定するために非線形関数を用いる。ウェルに対して得られる最良適合を用いて、該ウェルに対するベースラインを確立する。強度をプロットする工程、開始および終了サイクルを決定する工程、最良適合を決定する工程、およびベースラインを確立する工程は、複数のウェルのそれぞれに対して個別に実施される。
【０００９】
１つの態様において、開始サイクルと終了サイクルとを決定する工程は、サイクル２から始まる各サイクルに対して得られる実シグナルの強度の３点移動平均を取得することと、前記３点移動平均をサイクル２から始まるサイクル数の関数としてプロットして、第２のプロットを得ることとを含む。
【００１０】
別の態様において、サイクルＸは、第１または第２のプロットから特定される。サイクルＸは、直前のサイクルからの勾配の変化が最大となるともに、直後にそれぞれ直前のサイクルと比べて勾配が増大する４つのサイクルを有している。１つの態様において、サイクルＸは、ＰＣＲ反応の最初の８サイクル以内で特定され、終了サイクルは、サイクルＸの３サイクル、２サイクルまたは１サイクル前のものとして特定される。しかしながら、他の態様において、サイクルＸは、ＰＣＲ反応の最初の８サイクル以内では特定されない。この態様において、サイクルＸは、終了サイクルとして割り当てられる。
【００１１】
さらなる態様において、開始サイクルは、直前のサイクルとは異なる勾配傾向を示す第１のサイクルと、直前のサイクルの勾配よりも大きな正の勾配を有するサイクルと、初期勾配の１０％未満の勾配しか有しないサイクルとから成る群より選ばれる。
【００１２】
１つの態様において、最良適合線は、以下の式で表される、標準最小二乗平均平方根誤差アルゴリズムを用いて算出される。
【数１】

【００１３】
（式中、Ｎは前記開始サイクルと終了サイクル間のサイクル数であり、Ｓは前記開始サイクルであり、Ｅは前記終了サイクルであり、γ_ｉは第ｉサイクルのシグナルであり、ｍおよびｂは前記最良適合線のパラメータである）。
【００１４】
別の態様において、試料ブロック中の試料ウェルの少なくとも１つが、鋳型核酸を含む可能性のある試料を含み、光学シグナルは、該核酸に特異的な標識（例えば蛍光標識など）に由来する。さらなる態様において、本方法は、工程（ｄ）として、少なくとも１つのウェルに対してプロットされた第１のプロットから、少なくとも１つのウェルに対して算出された適応ベースラインを差し引いて、適応シグナルプロットを作成する工程を含む。好ましくは、適応シグナルプロットは、既知量の鋳型濃度を含むＰＣＲ反応から作成された第１の標準プロットと比較される。より好ましくは、標準プロットは、ベースラインを除去するように適応されている。この比較を用いて試料中の鋳型の量を決定することができる。
【００１５】
本発明は、上述の方法を実施するためのプログラム命令を含むコンピュータプログラム製品も提供する。
【００１６】
本発明は、ＰＣＲ反応に対する適応ベースラインの算出のための装置をさらも提供する。装置は、上述のコンピュータプログラムのプログラム命令を実行するためのメモリを含む解析装置も含む。１つの態様において、装置は、試料ブロックを受容するためのサーマルサイクラー、試料ブロックの試料ウェルから光学シグナルを収集するための入力装置（例えば検出器）、および、前記解析装置によって取得または作成されたデータを表示するための出力装置（例えばコンピュータ）をさらに含んでいる。
【発明を実施するための最良の形態】
【００１７】
本発明の目的および特徴は、以下の詳細な説明および図面を参照することにより、よりよく理解されるであろう。
【００１８】
本発明は、ＰＣＲ反応の適応ベースラインの算出のためのアルゴリズム、該アルゴリズムを含むプログラム製品、該アルゴリズムを実行するための装置、およびその使用方法に関する。本アルゴリズムは、増幅反応において評価される１つ以上の標識から得られるデータに対する適応ベースラインを算出する。１つの態様において、本方法は、ＰＣＲサイクル中に得られる実シグナルを測定し、１つ以上の標識のそれぞれに対して第１のプロットを作成し、１つ以上の標識に対する開始サイクルと終了サイクルとを決定し、各試料または標識に対する適応ベースラインを算出する。適応ベースラインは、ＰＣＲ増幅サイクル中に観察される実シグナルから差し引いて、ＰＣＲ反応中の核酸鋳型の量の算出のために用いることもできる。
【００１９】
定義
次に、以下の記述において用いる具体的な用語について定義する。
【００２０】
「核酸」は共有結合されたヌクレオチド配列（すなわちＲＮＡに対してはリボヌクレオチド、ＤＮＡに対してはデオキシリボヌクレオチド）であり、ここで１つのヌクレオチドのペントースの３’位が、ホスホジエステル基によって、次のペントースの５’位に結合されている。「核酸」には、限定はされないが、一本鎖および二本鎖の核酸が含まれる。本明細書中で用いる「核酸」という用語は、化学的、酵素的、または代謝的に修飾された核酸の形態も含意する。
【００２１】
本明細書中で用いる用語「リアルタイム標的鋳型合成」または「リアルタイム合成」は、合成産物を、その生成の最中に、該産物のそれ以後の合成に影響を及ぼすことなく解析するこのできる合成プロセスのことをいう。
【００２２】
本明細書中で用いる、「ポリメラーゼ連鎖反応」または「ＰＣＲ」とは、標的鋳型配列を増幅するためのインビトロ法のことをいう。ＰＣＲ反応は、反復する一連の温度サイクルを伴い、典型的には５０〜１００μＬの体積で実施される。反応混合物は、ｄＮＴＰ（４つのデオキシヌクレオチドｄＡＴＰ，ｄＣＴＰ，ｄＧＴＰおよびｄＴＴＰ）と、プライマーと、バッファーと、ＤＮＡポリメラーゼと、ポリヌクレオチド鋳型とを含む。１つのＰＣＲ反応は、変性とポリヌクレオチド分子合成とからなる「サイクル」が５〜１００回で構成されてもよい。
【００２３】
本明細書中で用いる「サイクル」という用語は、選択された時間にわたる一連の温度工程のことをいう。１つの態様において、サイクルは、核酸を変性させる工程と、核酸にプライマーをアニーリングする工程と、アニールされたプライマーを伸長する工程とからなる３工程を含む。「工程」は、均一な温度（±１℃以内）の期間のことをいいう。しかしながら、ＰＣＲサイクルを循環させられる試料ウェルは、変性や、プライマーへのアニーリングなどを受けうる核酸を有している必要はなく、例えば、ウェルはバッファだけ、またはヌクレオチドだけ、および／またはＰＣＲ反応混合物の他の様々な成分を含んでいてもよい。試料ウェルは、変性、アニーリングおよび伸長に適した同一の時間量のあいだ同一の温度に曝されていさえすれば「循環されている」と言える。１つの態様において、ＰＣＲ反応は、複数のサイクルを含み、各サイクルが、９０〜１００℃（好ましくは９４℃）で３０秒間〜１分間（好ましくは３０秒間）の工程と、３７〜６０℃（好ましくは５５℃〜５７℃）で１〜２分間（好ましくは１分間）のアニーリング工程と、７０〜７５℃で３０秒間〜１分間（好ましくは３０秒間）の伸長工程とを含む。
【００２４】
本明細書中で用いる「終了サイクル」とは、定量的ＰＣＲ反応（「ＱＰＣＲ」）の間に標的鋳型の増幅が検出可能になるまえのサイクルのことをいう。「終了サイクル」は、試料中の標的鋳型の増幅を検出するために用いる標識に対して決定してもよい。
【００２５】
本明細書中で用いる「開始サイクル」とは、初期の増幅の変動が弱まった後であって、終了サイクルよりも前のサイクルのことをいう。１つの態様において、開始サイクルは、直前のサイクルとは異なる勾配傾向を示す第１のサイクルと、直前のサイクルの勾配よりも大きな正の勾配を有するサイクルと、最初の勾配の１０％未満の勾配しか有しないサイクルとから成る群より選ばれる。
【００２６】
本明細書中で用いる「検出可能な」は、アルゴリズムが、あるサイクルにおけるシグナルがバックグラウンドを超えていると決定したこと、すなわち、所与の試料（例えば、試料中の所与のラベル）に対して決定されたシグナルが、シグナルのバックグラウンドレベルを超えて、正の増加（例えば、少なくとも１０％、２０％、３０％、４０％、５０％またはそれ以上、あるいは２倍、３倍、４倍またはそれ以上の増加）を示していることを意味する。
【００２７】
本明細書中で用いる「バックグラウンドシグナル」とは、終了サイクルより前のシグナル、あるいは、鋳型を含まない対照試料において発生されたシグナルのことをいう。
【００２８】
本明細書中で用いる「正の増加」は、終了サイクルを、バックグラウンドシグナルが変動するサイクルから区別するために用いられる。シグナルの「正の増加」は、例えば、ランダムなバックグラウンドノイズ変動ではなく、標識から発生されるシグナルによって示されるような、試料中の増幅産物の増加を反映する。典型的には、正の増加がある場合には、蛍光シグナルは少なくとも２（または３、または４、あるいはそれ以上）の連続したサイクルに亘って継続的に増加している。
【００２９】
本明細書中で用いる用語「サイクルＸ」は、ＱＰＣＲ反応から発生されるシグナルが初めて閾値を超える、すなわち、標的核酸配列の増幅を初めて信頼性よく検出できるようになるサイクル数のことを言う。「信頼性ある」とは、シグナルが、ＱＰＣＲの最中に増幅された産物の検出可能なレベルを反映していることを意味する。サイクルＸは一般には、標的核酸の未知量の出発量に相関し、例えば標的の量が少ないほど、サイクルＸは遅くなる。
【００３０】
本明細書中で用いる用語「変動」とは、試料中の標的核酸の量には直接関係しない、増幅されたＰＣＲ産物の品質または量の変化のことをいう。変動は、試料の取扱い時のムラや使用する反応容器の違いによることもある。例えば、チューブを用いた場合、生シグナルはチューブによって変動する（例えば、「チューブ間変動」）。本発明に従う変動は、増幅された核酸物から発生するシグナルの有無に関係なく、各試料から得られる生データに終始変動を与える機器ドリフトに起因することもある（例えば、「サイクル間変動」）。いずれのタイプの変動も、サイクルＸの測定の精度を低下させる。
【００３１】
本明細書中で用いる「サイクル間変動」は、リアルタイム合成の任意の２つ以上のサイクル間での増幅または合成の変動のことをいう。
【００３２】
「成長カーブ」という用語は、ＰＣＲ反応中の各サイクルの伸長部分の終点または終点付近で存在する、増幅された核酸物（例えばｄｓＤＮＡ）の測定値の組を意味する。
【００３３】
本明細書中で用いる「実シグナル」とは、標識によって発生される検出可能なシグナルのことをいい、この場合シグナルは、シグナル検出器によって直接測定することができる。好ましくは、検出可能なシグナルは光学シグナル（例えば、蛍光標識または化学発光標識）であり、シグナル検出器は光学シグナル検出器である。
【００３４】
本明細書中で用いる「標識」とは、検出可能なシグナル、例えば光学シグナルを発生する分子のことをいう。本発明に従う有用な標識タイプとしては、限定はされないが、蛍光、化学発光、比色または酵素標識が挙げられる。
【００３５】
本明細書中で用いる「標識されたプローブ」とは、増幅産物に直接的または間接的に取り込まれるか、結合することのできる分子のことをいう。例えば、標識されたプローブは、増幅産物または標識されたプライマー中に取り込まれた標識されたヌクレオチドであってもよいし、あるいは、標識されたプローブは、それ自体の増幅産物に結合する標識された分子であってもよい（例えば、増幅産物中の内部配列に対するプローブなど）。
【００３６】
本明細書中で用いる用語「蛍光」とは、所与の波長の光を照射すると、分子が励起してそれよりも長い波長の光を発するという分子の特性をいう。
【００３７】
本明細書中で用いる用語「フルオロフォア」とは、蛍光分子のことをいう。フルオロフォアの蛍光特性を記述する多数のパラメータが存在する。これらの中には、例えば、励起および発光の最大波長、励起および発光に対するピークの幅、励起および発光の最大間の差異（ストークスシフト）、蛍光強度、量子収量、および吸光係数が含まれる。生物学的または生化学的用途に対しては、ストークスシフトは一般に短いものよりも長いが好ましい。
【００３８】
本発明に従う「標的鋳型」は、試験試料においてその正体（たとえば）または量を決定すべき鋳型のことをいう。「標的鋳型」は、複製および／または検出すべきポリヌクレオチド鋳型の一領域であってもよい。１つの態様において、「標的鋳型」は、ＰＣＲ増幅のために用いられる２つのプライマー配列間に位置する。「標的鋳型」は、ＰＣＲ反応において生成した増幅された産物であってもよい。しかしながら、試験試料は、正体または量を決定する必要のある２つ以上の標的鋳型を含んでいてもよく、そのような標的鋳型のそれぞれが、解析される際に標的鋳型であってもよい。本発明に従えば、試料中に複数の標的鋳型が存在する場合、それぞれが特異的標的鋳型に特異的に結合し、したがって特異的標的鋳型に対する「マッチした標識」と称される複数のラベルであってもよい。
【００３９】
本明細書中で用いる「標的鋳型を検出する」とは、試料中の所与の標的ポリヌクレオチド配列の有無を決定することをいう。
【００４０】
本明細書中で用いる「標的鋳型を測定する」とは、試料中の所与の標的ポリヌクレオチド配列の量を決定することをいう。測定または検出可能な標的ポリヌクレオチド配列の量は、好ましくは約１分子から１０^２０分子であり、より好ましくは約１００分子から１０^１７分子であり、最も好ましくは約１０００分子から１０^１４分子またはそれ以上である。
【００４１】
本明細書中で用いる「オリゴヌクレオチドプライマー」とは、ポリヌクレオチド鋳型にハイブリダイズすることができ、第２のポリヌクレオチド鎖の合成（例えば、酵素的合成など）を開始する一本鎖ＤＮＡまたはＲＮＡ分子のことをいう。本発明に従って有用なオリゴヌクレオチドプリマーは、約６〜１００、好ましくは１０〜１００ヌクレオチド長であり、より好ましくは約１７〜５０ヌクレオチド長であり、最も好ましくは約１７４５ヌクレオチド長である。
【００４２】
「標準プロット」とは、既知量の標的鋳型の関数として、光学シグナルの強度（例えば蛍光強度）をプロットすることによって作成されるプロットのことをいう。「標準プロット」は、試料中の未知量の標的鋳型を定量するための基準となる。
【００４３】
本明細書中で用いる「試料」または「試験試料」とは、対象とする標的鋳型を含む任意の物質（例えば、標的ポリヌクレオチドまたは標的ポリペプチド）のことをいう。「試料」という用語は、細胞、組織、体液（限定はされないが、血漿；血清；髄液；リンパ液；滑液；尿；涙；大便；皮膚、呼吸器、腸および泌尿生殖器管の外分泌物；唾液；および血液）、腫瘍、器官、生物（例えば微生物）、インビトロ細胞培養構成分の試料、環境試料（例えば、湖、貯水池、土壌試料など）あるいは産業試料（例えば、市販食品製品、産業廃棄物など）のなかに見いだされる、ポリヌクレオチド（ゲノムＤＮＡ、ｃＤＮＡ、ＲＮＡ）および／またはポリペプチドの試料を含んでいてもよい。
【００４４】
ＰＣＲ手順の概略
本発明の装置の基礎は、新奇な適応ベースライン減算アルゴリズムを使用することにあり、該アルゴリズムは、例えばリアルタイムＱＰＣＲ中などの増幅反応中に、評価中の各試料の特性に基づいて増幅に先だつ開始およびお終了サイクルを決定する。アルゴリズムは、全ての試料に対する終点に基づいてベースラインを算出する従来技術とは違い、各試料に対して個々に決定した終点からベースラインを算出する。本発明に従う適応ベースラインアルゴリズムは、増幅産物の最初の信頼性ある検出を反映する閾値サイクル（サイクルＸ）のより正確および／または精密な測定値を与えるために用いることができる。
【００４５】
アルゴリズムは、ＰＣＲに加えて様々なタイプの増幅スキームの解析において用いることができる。一般に、アルゴリズムは、増幅された核酸物（例えばｄｓＤＮＡ）の集団が増加する増幅反応中に、核酸ポリメラーゼを使用する必要のある増幅スキームを評価するために用いられる。増幅スキームの例としては、限定はされないが、ＰＣＲ；リガーゼに基づく増幅スキーム、例えばリガーゼ連鎖反応（ＬＣＲ）；Ｑ−β−レプリカーゼに基づく増幅スキーム；鎖交換増幅（ＳＤＡ）スキーム（Ｗａｌｋｅｒら、ＮｕｃｌｅｉｃＡｃｉｄｓＲｅｓｅａｒｃｈ，第２０巻、１６９１〜１６９６頁（１９９２年）に記載されているような）等が挙げられる。核酸増幅スキームの総括的な記載は、ＫｅｌｌｅｒおよびＭａｎａｋ、ＩｎＤＮＡＰｒｏｂｅｓ，第２版（ＳｏｃｋｔｏｎＰｒｅｓｓ，ニューヨーク、１９９３年）に提示されている。
【００４６】
一般に、アルゴリズムは、反応混合物中に標識を含めることによってモニタされる増幅反応に対するベースラインの算出のために用いられる。標識は、（例えば、増幅産物中に取り込まれるか、あるいは増幅産物と結合することにより）増幅産物と相互作用することができ、例えば光学シグナルなどの検出可能なシグナルを発生することができる。このシグナルの測定値は、反応中の増幅された産物の量と相関させることができる。本発明の有用な標識タイプとしては、限定はされないが、蛍光、化学発光、比色、または酵素標識が挙げられる。
【００４７】
有用な蛍光ラベルの例としては、フルオレセイン、ローダミン、フィコエリトリン、フィコシアニン、アロフィコシアニン、およびフルオレサミンが挙げられる。有用な化学発光標識タイプとしては、例えば、ルミナル、イソルミナル、芳香族アクリジニウムエステル、イミダゾール、アクリジニウム塩、シュウ酸エステル、ルシフェリン、ルシフェラーゼ、およびエクオリンが挙げられる。有用な酵素としては、例えば、マレイン酸ヒドロゲナーゼ、ブドウ球菌デヒドロゲナーゼ、デルタ−５−ステロイドイソメラーゼ、アルコールデヒドロゲナーゼ、アルファ−グリセロールリン酸デヒドロゲナーゼ、トリオースリン酸イソメラーゼ、ペルオキシダーゼ、アルカリホスファターゼ、アスパラギナーゼ、グルコースオキシダーゼ、β−ガラクトシダーゼ、リボヌクレアーゼ、ウレアーゼ、カタラーゼ、グルコース−６−リン酸デヒドロゲナーゼ、およびグルコアミラーゼ、アセチルコリンエステラーゼが挙げられる。
【００４８】
適切な標識は、当業者によって知られる標準的な技術によってポリヌクレオチド、ポリペプチド、抗体または抗体断片に結合（例えば共有結合）させるか、取り込ませることができる。例えば、Ｋｅｎｎｅｄｙら、Ｃｌｉｎ．Ｃｈｉｍ．Ａｃｔａ、第７０巻、１〜３１頁（１９７６年）、およびＳｃｈｕｒｓら、Ｃｌｉｎ．Ｃｈｉｍ．Ａｃｔａ、第８１巻、１〜４０頁（１９７７年）参照のこと。
【００４９】
本発明の１つの実施形態において、ＰＣＲ中の産物の増幅をモニタするために蛍光標識された分子を用いる。蛍光標識された分子は、増幅反応（例えばＰＣＲ反応）によって生じた増幅産物と相互作用するかこれに取り込まれた場合、（例えば、少なくともＰＣＲ反応の各サイクルの伸長部分の間に）適当な刺激が与えられると検出可能なシグナルを発生する。増幅反応を十分な時間行い（例えば、ＰＣＲ反応の間に十分なサイクル数の間反応混合物に対して熱サイクルを行う）、所望の最終濃度の増幅核酸物を確立する。
【００５０】
増幅反応の任意の時点で観察される生シグナル（例えば蛍光強度）を、本発明のアルゴリズムを用いて相殺し、シグナルが観察されている特定の試料に独特なバックグラウンドシグナルを除去する。得られたシグナルは、試料中で増幅される標的分子のより実際に近い測定値を与える。１つの態様において、バックグラウンドシグナルまたはベースラインは、標的核酸の非存在下における、試料を含むチューブまたはウェル中の予測されるシグナルとして決定される。シグナルベースラインが評価中のすべての試料に対して決定される従来の方法とは対照的に、本発明はアルゴリズムを用いて、増幅をモニタ中の個々の試料それぞれに対して特異的なベースラインを決定する。
【００５１】
１つの態様において、アルゴリズムはＱＰＣＲ反応の最中に実施される。典型的なＱＰＣＲのプロセスには、ＱＰＣＲ中に増幅された産物の量に対応するシグナルを発する少なくとも１つの標識存在下において、標的鋳型を増幅することと、発生されたシグナルを収集することと、収集されたシグナルをプロセッサに保存することが伴う。収集されたシグナルを解析し、続いて前記プロセッサと連絡するユーザ装置のスクリーンに表示する。
【００５２】
１つの実施形態において、本発明のアルゴリズムを実行するための方法は、選択された複数の開始サイクル（増幅がどのチューブ内でも起こる前）のそれぞれの少なくとも伸長部分の間に、試料中のシグナルの強度を検出および測定する工程と、該強度を増幅産物（例えばｄｓＤＮＡ）のモル濃度値に変換する工程と、各サイクルの各伸長部分に対するモル濃度値をプロセッサに保存する工程と、保存された濃度値に基づいて、増幅産物（例えばｄｓＤＮＡ）のモル濃度を表す測定曲線を、サイクル数の関数として作成する工程とを含む。
【００５３】
シグナル強度測定値からモル濃度への変換は、当該技術分野における常套的方法によって実施することができる。例えば、ＰＣＲ試料を、ゼロ（例えばベースライン）から試料中の未知量の標的核酸を検出するためのＰＣＲ反応によって産生されそうな濃度範囲の既知濃度の核酸鋳型を含むこと以外は、すべての同一の反応材料を同一の濃度で含んだＤＮＡの標準溶液（例えば光度標準）に置き換えることができる。そして、これらの既知の光度標準上で測定されたシグナルからベースラインシグナルを減算することにより、ベースラインバックグラウンドに対して補正された測定シグナルと、標準鋳型ｄｓＤＮＡの濃度との関連を示した作業曲線が作成される。ベースライン減算は、シグナル強度測定値をモル濃度に変換する前または後に実施することができる。１つの実施形態において、ベースライン減算は、シグナル強度測定値をモル濃度に変換する前に実施される。別の実施形態において、ベースライン減算は、シグナル強度ををモル濃度に変換した後に実施される。この作業曲線を準備するためにＰＣＲを実施する必要はない。既知のＤＮＡモル濃度を有する標準溶液を用いてもよい。
【００５４】
好ましい実施形態において、強度をモル濃度値に変換することなく、保存された濃度値に基づいて、増幅された産物のシグナル強度（例えばｄｓＤＮＡ）を、サイクル数の関数として表した測定曲線を作成した。
【００５５】
１つの態様において、未知量の核酸を含む試料中の標的核酸の出発モル濃度を、試験試料のアリコートを得、アリコート上の同一の増幅反応を、既知量の標的核酸を含む標準試料に対して用いられたものと同一の条件下において実施することによって決定する。増幅反応における各サイクルの伸長部分の最後またはその付近で存在する増幅された核酸物（例えば、ｄｓＤＮＡ）の測定値のセットを得ることにより、アリコートのそれぞれに対して成長曲線を記録する。アリコートの１つの成長曲線と、標準に対する作業曲線との間の最良適合は、例えば、一変数逐次近似法などによって算出される。作業曲線が、未知量の核酸を含む試料のアリコートの成長曲線と最良適合する標準の濃度を、未知試料のアリコートの濃度として取りだす。ベースライン減算は、未知量の標的核酸を含む試料に対して発生した測定シグナルに対して、以下により詳細に記載するように、シグナル強度測定値をモル濃度に変換する前または後に適用してもよい。
【００５６】
適応ベースラインアルゴリズムを用いた方法
リアルタイムＰＣＲの間に、リアルタイムＰＣＲ装置（Ｓｔｒａｔａｇｅｎｅ（商標）のＭｘ４００サーマルサイクラーなど）に挿入した複数の試料ウェルホルダの所与の試料ウェルのそれぞれに対する、実シグナルまたは生データを、ＰＣＲ増幅反応のすべてのサイクルにおいて検出する。核酸鋳型を除きＰＣＲ増幅に用いた全ての試薬を含めることにより、１つ以上のウェルを、対照ウェルとして用いてもよい。１つ以上の試料ウェル中の１つ以上の標識からの各サイクルにおける実シグナルをデータセット中に収集することができる。データセットは、１つ以上の標識のそれぞれの平均（すなわち算術平均）シグナルからなってもよいし、あるいは、各サイクルに対して標識から発生した最小または最大のシグナルからなるものであってもよい。
【００５７】
従来のＰＣＲベースライン減算アルゴリズムは、同一の開始および終了サイクルをＰＣＲ反応中のすべてのウェルに対して割り当てている。しかしながら、ＰＣＲ増幅がサイクル毎およびウェル毎で変わることは不可避である。したがっていくつかのウェルに対して、終了サイクルは大きすぎることになる。この場合、増幅プロットは、負の大きな勾配で始まり、予定より遅い増幅を示し、増幅におけるより遅いサイクルに亘って著しく垂下する（図１に示されるように）。他のいくつかのウェルに対しては、終了サイクルが小さすぎるため、増幅プロットが正の大きな勾配で始まり、予定よりも早い増幅を示す（図１参照）ということも起こる。終了サイクルが小さすぎると、初期勾配があまり影響を受けないが、測定された増幅、具体的かつ最も重要なことにはサイクルＸの精度がはるかに落ちてしまう。
【００５８】
本発明において、開始および終了サイクルは、各試料ウェルおよび／または評価中の各染料に対して個別かつ独自に決定される。各試料ウェルおよび／または染料に対して、試料ウェル中の標識によって発生される実シグナルをプロットすることにより、サイクル数の関数としての第１のプロットを作成する。好ましい実施形態において、サイクル数の関数として、実シグナルの３点移動平均（「ＭＡ」）をプロットした第２のプロットを作成し、サイクル数の関数として３サイクルごとに得られた一連の平均を得る。３点ＭＡの算出方法を表１に示す。
【表１】

【００５９】
算出された３点移動平均は、適応ベースラインのための終点を決定するために用いてもよく（例えば、実際の増幅が起こる前の開始および終了サイクルが、評価中の個別のウェル中で起こる）、各データセット、すなわち各試料ウェルに対して算出される。
【００６０】
好ましい実施形態において、開始サイクル（増幅反応の初期変動が消失した後、検出可能な増幅が始まる前のサイクル）は、以下の１）〜３）のいずれかに１つとして割り当てられる。１）直前のサイクルとは異なる勾配傾向を示す第１のサイクル。例えば、直前の勾配が負であったのに正となる、あるいは直前の勾配が正であったのに負となるようなサイクル。２）直前の勾配よりも大きな正の勾配を有するサイクル。３）第１または第２のプロットの勾配である初期勾配の１０％未満の勾配しか有しないサイクル。どのサイクルも上記１）〜３）の基準を満たさない場合は、開始サイクルはサイクル０として割り当てられる。
【００６１】
他の好ましい実施形態において、開始サイクルは、最初のサイクルと８番目のサイクルとの間のサイクルに限定される。この限定により、データの広がりを少なくし、および／または人工的疑似結果の発生を低減することができるかもしれない。
【００６２】
別の実施形態において、アルゴリズムは、一次微分および二次微分によってよりむしろデータ内のノイズレベルに基づいている。この場合、アルゴリズムは、提案されたベースライン付近の点の標準偏差を算出する。そして、サイクル６までの点をこの値と比較する。サイクル６までの全てのポイントがこの値の外にあれば、サイクル７が開始サイクルである。もしそうでなければ、アルゴリズムはサイクル５までの点を参照する。アルゴリズムは条件が合うまでこれを続けるが、サイクル１はアルゴリズムによっては解析されない。
【００６３】
評価中の個々の試料に対する終了サイクルを決定するために、サイクルＸを決定することが重要である。サイクルＸは、増幅が信頼性よく検出され始めるサイクルのことをいう。サイクルＸは、直前のサイクルからの勾配の変化が最大となるサイクルのことをいい、サイクルＸの直後には、それぞれに直前のサイクルと比較して勾配が増大する４サイクルが続く。サイクルＸを決定するための工程は、サイクル回数に対する蛍光強度のプロット上の点を分類することと、プロット上の異なる点における一次微分および二次微分を確認することを含んでいてもよい。一次微分は、あるサイクルから次のサイクルまでに観察される蛍光強度の変化に相当する勾配を表す。二次微分は、増幅反応の間の変化を表す。
【００６４】
サイクルを表す任意の点において、一次微分が正であり、次の４つの点が増加勾配を有する場合、このサイクルはサイクルＸとして認められる。もしこの基準が満たされなければ、そのサイクルは排斥され、サイクルに対応する次の点を評価して、そのサイクルがサイクルＸに対する基準に合うかどうかを判定する。
【００６５】
１つの実施形態において、サイクルＸは、ＰＣＲ反応の最初の８サイクルを表す点を評価することによって決定される。サイクルＸが最初の８サイクル以内で合格すれば、終了サイクル（標的鋳型の増幅前におこるサイクル）は、サイクルＸから一定の任意のサイクル数を差し引くことによって割り当てられる。この差し引かれる一定の数は、実験毎に変わるものであってもよいし、後述するように、評価中の試料に対して適応ベースラインがどの程度よく適合しているかに従って決定することもできる。好ましい実施形態において、終了サイクルはサイクルＸの３、または２、または１サイクル前である。
【００６６】
別の実施形態において、サイクルＸは、ＰＣＲ反応の最初の８サイクル以内で決定される。この実施形態において、サイクルＸ自身が終了サイクルとして割り当てられる。
【００６７】
別の実施形態において、サイクルＸは、増幅曲線の最大二次微分に基づいて決定される。サイクルＸは、最大二次微分におけるシグナルの固定分数として設定される。さらに別の実施形態において、サイクルＸは、例えば、最大二次微分と最大勾配、あるいは最大二次微分と最小二次微分などの２つのシグナルレベルに基づいて決定される。さらに別の実施形態において、サイクルＸは、適応ベースラインによって決定される終了サイクルに基づいて決定される。
【００６８】
１つの態様において、サイクルＸは、複製試料（同じように調製され、複製であると機器によって確認されている試料）の群に対するサイクルＸの範囲に基づいて決定される。これらの試料が同一であるので、これらは同一の閾値サイクルを与えるに違いない。この場合、サイクルＸは、シグナルレベルが、全ての複製に対する閾値サイクルの範囲が最小になるサイクルとして決定される。
【００６９】
終了サイクルが正確に選択されれば、データが指示範囲から出る範囲もはるかに小さくなるであろう。終了サイクルを決定する適応ベースラインアルゴリズムは、増幅曲線の勾配（二次微分）の最大変化に基づいている。１つの実施形態において、アルゴリズムは３つのサイクルを差し引くことにより、終了サイクルを作成する。他の実施形態において、アルゴリズムは、３サイクルではなく４サイクルを差し引くことにより、終了サイクルを算出する。疑似終点を除去するために、様々な検査を適用することができる。例えば、シグナルではなくノイズが最大二次微分を与える場合などである。１つの実施形態において、最大二次微分を有するサイクル上の終了サイクルは、一次微分の兆候と、そのサイクル付近の一次微分と二次微分の傾向をチェックすることによって実施した。
【００７０】
個々のウェルに対して決定される開始および終了サイクルは、ベースラインが、機器ドリフトおよび／またはサンプリングのばらつきによるＰＣＲ反応の変動に、より正確に適合するように、さらに調整してもよい。１つの実施形態において、算出された成長曲線は、鋳型に対して作成される。この成長曲線は、既知量の鋳型を用いた各サイクルにおいて発生されるシグナルを反映している。開始および終了サイクルは、既知または未知量の鋳型を使用して、上述のようにＰＣＲ増幅反応に対して決定してもよい。増幅の実シグナルが測定される。測定された成長曲線は、上述のようにして決定された開始および終了サイクルに基づいて誘導された適応ベースラインを用いて実シグナルから作成してもよい。そして、測定成長曲線を、理論成長曲線と比較して、測定成長曲線と理論成長曲線の間の適合の質に従って、試料に対する開始および終了サイクルを調整できるようにする。
【００７１】
測定成長曲線と理論成長曲線の間の適合の質は、差異の二乗の平均をどれほど小さくできるかによって測ることができる。有用な手段としては、適合の二乗平均平方根誤差、すなわち「ｒｍｓｅ」である。ｒｍｓｅは、各サイクルにおける測定シグナル強度と予測（例えば理論）値の差異を二乗することによって算出される。適合の範囲（例えば、適合に含まれるサイクル数）に亘ってこれらの二乗の平均を算出し、平方根を求めると、適合に対するｒｍｓｅが得られる。
【００７２】
１つの態様において、ｒｍｓｅは、下記の式に従って決定される。
【数２】

【００７３】
（式中、Ｎは、開始サイクルと終了サイクルの間のサイクル数であり、Ｓは開始サイクルであり、Ｅは終了サイクルであり、γ_ｉはｉ番目のサイクルであり、ｍおよびｂは最良適合線のパラメータである）。
【００７４】
ｒｍｓｅはシグナルがバックグラウンドを超えて増幅し始めると迅速に増加する。１つの態様において、サイクルに対する適合のｒｍｓｅが適合するように、終了サイクルを変えて、サイクルに対する適合のｒｍｓｅを適合するようにプロットする。ｒｍｓｅが所定の値を超える終了サイクルが選ばれる。等価な結果を与える任意の方法が有効であることは当業者にとって自明である。例えば、観察された差異の重大性を判定するために当該技術分野において常套的に用いられている統計的方法を用いることができる。
【００７５】
各増幅反応に対する終了サイクルを決定する別の方法は、ベースライン関数の最小平均二乗適合のｒｓｑを、適合サイクル関数としてのデータにプロットすることである。終了サイクルは、ｒｓｑ値が所定の値を超えるサイクルとして定義される。
【００７６】
測定成長曲線および理論成長曲線は、それらの多項または対数形あるいは他の適当な形態で比較してもよい。
【００７７】
試料に対する開始および終了サイクルの別の調整方法としては、終了サイクルを適宜割り当てることである。１つの実施形態において、開始サイクルと終了サイクルの間の相対距離が０である場合、例えば、開始サイクルと終了サイクルが同じサイクルである場合、終了サイクルは、適宜、増幅反応の最終サイクルに割り当てられる。開始サイクルは適宜、サイクル０に割り当てられる。
【００７８】
別の実施形態において、開始サイクルと終了サイクルの間の相対距離が０より大きく（すなわち、開始サイクルと終了サイクルが同一でない場合）、所定のサイクル数（例えば、６〜１２サイクル）よりは小さい場合、開始サイクルと終了サイクルは、開始サイクルと終了サイクルの間の距離が所定の数（例えば、６または７または８または９または１０または１１または１２）となるように適宜割り当てられる。
【００７９】
好ましい実施形態において、一定の数は適宜８に決定される。たとえば、相対距離が８サイクル未満であれば、開始および終了サイクルは８サイクル離れて設定される。このことは、終了サイクルの数に対しての開始サイクルの数を調整することによって達成される。したがって、もし終了サイクルがサイクル１０であり、当初の開始サイクルがサイクル５である場合には、開始サイクルは、終了サイクルからの距離８を維持するためにサイクル２に調整される。終了サイクルが８以下のサイクルである場合には、開始サイクルは、サイクル０に設定され、終了サイクルは、開始サイクルと終了サイクルの間の距離８を維持するために、サイクル８に設定される。
【００８０】
別の実施形態において、それぞれの相対距離に基づいて開始サイクルと終了サイクルを調整する前に、開始サイクルと終了サイクルは、測定成長曲線上の開始サイクルと終了サイクルの間の所定のサイクルにおける蛍光強度を表す点を、理論成長曲線上の点と比較することにより、標準偏差に基づいて調整される。開始サイクルと終了サイクルは、すべての点が互いに何らかの標準偏差内に収まるまで、内方（互いに向けて）調整することができる。
【００８１】
本方法の実施にあたって実行すべきサイクルの数は、少なくとも試料中の増幅産物（例えばｄｓＤＮＡ）の蛍光シグナルの強度を、目的とする最小開始濃度を有した試料に対する各サイクルの終了時において、蛍光強度測定に用いる手段（例えば、サーマルサイクラーの検出器）による検出が可能な限界を超える値まで高めるのに十分でなければならない。サイクル数が大きければ大きいほど良好な結果が得られるが、一旦、増幅産物のモル濃度の各サイクルにおける有意な増加が止まってしまうと、サイクルを追加してもその利益は減る。
【００８２】
本発明を定量的リアルタイムＰＣＲに適用するために必要な装置
本発明のアルゴリズムおよび方法を適用するために有用な装置は、サーマルサイクラー、シグナルを収集するための入力装置（例えば、検出装置）、シグナルを解析するための解析装置（例えばプロセッサ）、および解析装置によって取得あるいは生成されたデータを表示するための出力装置（例えば、解析装置と連絡するコンピュータなどのユーザ装置）とを含んでいてもよい。解析装置は、当該技術分野において周知のように入力装置を介してサーマルサイクラーに接続してもよいし、および／またはデータ表示装置の技術分野において周知の出力装置に接続されているか、該装置を含んでいてもよい。１つの実施形態において、解析装置は、例えば図１２Ａおよび１２Ｂに示すようなコンピュータである。
【００８３】
本方法を実施するために用いられるサーマルサイクラーは、標準的なＰＣＲチューブ内、またはＰＣＲ試料プレート内でサーマルサイクリングブロック中に９６個までの反応試料を保持できる（本明細書中では総括的に「試料ウェルを含む試料ブロック」と呼ぶ）従来の設計のものであってもよい。チューブを用いる場合には、透明のキャップ（例えば、ポリプロピレンまたは他の透明の熱耐性プラスチックを含む）を有することが好ましい。好適には、チューブ／ウェルの上にくるサーマルサイクラーの蓋は、検出装置の入力装置を受容および支持するように特別に適応されている。
【００８４】
１つの態様において、検出装置は、１つ以上の蛍光標識から光学シグナルを検出するのに適している。好ましくは、前記装置は、アセンブリ内の試料チューブとの光の送受を行うための１つ以上の走査用光学ファイバを含んでいる。好ましくは、複数の光学ファイバ（たとえば、検出中のシグナルまたは標識のタイプ毎に１つ）を用いて、チューブの表面を走査する（例えば行毎に）。１走査毎に複数の光学シグナル（好ましくは、少なくとも９光学シグナル）を収集し、これを平均して１走査あたりの試料チューブに対する平均シグナルをもたらす。好ましくは、検出器は、光信号を、アセンブリに連絡したプロセッサへのデジタル入力に変換するための、適当な信号増幅変換回路も含んでいる。検出器は、Ｈｉｇｕｃｈｉら（１９９３年）Ｂｉｏｔｅｃｈｎｏｌｏｇｙ、１１（９）、１０２６〜１０３０頁に記載されているようなデジタルカメラ装置を含んでいてもよい。
【００８５】
検出装置の出力（例えば、増幅反応中に生じるものに対応するシグナル）は、プロセッサに送られて、データの保存と操作が行われる。１つの実施形態において、前記装置は、複数の異なるタイプの光学信号、例えば、複数の異なるタイプの蛍光標識を検出するとともに、ミクロプレート蛍光読取器の能力も有している（例えば、個々のウェルから得られたシグナルの強度を単離および解析することができる）。
【００８６】
検出装置は、好ましくは、可視光レーザーまたは紫外線ランプまたはハロゲンランプのいずれであってもよい励起光源と、反応チューブからの蛍光を受容するためのファイバ光学系とコネクタを介して、個々の反応チューブに励起光を分配するための多重化装置と、励起光から波長によって蛍光を分離するための濾過手段と、蛍光強度を測定するための検出手段とを含む多重蛍光光度計である。
【００８７】
好ましくは、サーマルサイクラーの検出装置は、３５０ｎｍ〜８３０ｎｍの検出範囲を提供し、フルオロフォア選択に大きな自由度を与え、高い感受性と優れたシグナル−ノイズ比を与える。装置の光源は好ましくは、３５０〜７５０ｎｍの広い範囲の励起光を発生する。これにより、ユーザは、スペクトルの重複が殆どもしくは全くなく、優れた多重化のための明瞭な線状シグナルを与えることができるフルオロフォアを選ぶことができる、強度評価の対象となる各フルオロフォアの励起および発光波長を正確に一致させて、スペクトル上隣接するフルオロフォアからの不要なクロストークを阻止するために、最適化インターフェースフィルタを設けることもできる（例えば、図６に示されるような）。例えば、ＦＡＭ、ＴＥＴ，ＨＥＸ／ＪＯＥ／ＶＩＣ、ＴＡＭＵＲＡ、テキサスレッド／ＲＯＸ、Ｃｙ５およびＣｙ３フィルタセットが市販されており、他のフルオロフォアに対してはカスタムフィルタセットを作成することもできる（Ｓｔｒａｔａｇｅｎｅ、カリフォルニア州）。
【００８８】
好ましくは、リアルタイム増幅プロットは増幅の進行としてみなされる。これにより、アセンブリのユーザは、実行が終了するまで待たなくても、サーマルサイクリング中の任意の時点において実験がどのように運んでいるかを一目で判定することができる。ユーザは反応中に問題が起こった場合に実行を中断するか、実験を終了し、所望の情報が生成され次第データを保存するかのいずれかを選択することができる。
【００８９】
検出装置によって受信された光学シグナル（例えば、所与のチューブにおける所与の時点における蛍光シグナル強度に対応）は、一般にプロセッサによって演算可能な信号に変換されて、該プロセッサと連絡したユーザー装置（例えばコンピュータ）の表示装置上でユーザが見ることのできるデータを与える。好ましくは、解析装置のプログラムによって、コンピュータへの電力損失時や通信エラー時であってもデータ収集が継続するというような、サーマルサイクラー装置とユーザ装置の間の強力な通信が可能になる。通信が再開されると、実行からのデータがサーマルサイクラーの内蔵ソフトウェアからコンピュータ上のソフトウェアに送られて、確実に実験が完結しデータはしっかり保存される。
【００９０】
ユーザ装置は、ユーザインターフェースを含んでいてもよく、あるいは、キーボードとビデオモニタを備えた従来の市販のパーソナルコンピュータ（ＰＣ）装置であってもよい。ユーザ装置によって表示しうるデータとしては、例えば、増幅プロット（図８を参照）、散布図（図１０を参照）、アセンブリ中のすべてのチューブおよびすべての使用標識に対する試料値スクリーン、光学シグナル強度スクリーン（例えば、蛍光シグナル強度スクリーン）、最終案内結果、融解曲線、アニーリング範囲、レポート文書などが挙げられる。好ましい態様において、ユーザ装置は作業曲線、成長曲線、および第１および第２のプロット、ならびにこれらの曲線／プロット上での様々な演算の結果を表示することができる。例えば、１つの態様において、ユーザ装置は、解析装置／プロセッサが第１のプロットと第２のプロットの間の適合を判定した後に得られるｒｍｓｅを表示する。
【００９１】
ユーザ装置は、ユーザが解析装置またはプロセッサに、例えば、装置のサイクリングパラメータを変更する指示またはベースライン減算を実行する指示などの指示を与えることができるようなユーザインターフェースを表示することもできる。１つの態様において、ユーザ装置からの指示に応答して、解析装置またはプロセッサは、本発明の適応アルゴリズムを含むプログラムを実行する。アルゴリズムは、例えば、コンピュータによって使用されるソフトウェアの一部または、コンピュータのメモリ内またはコンピュータがアクセス可能なサーバーのメモリ内に記憶されているアプリケーションの一部としての、コンピュータによって使用可能なプログラム製品の一部であってもよい。
【００９２】
１つの態様において、解析装置は、上述した本発明の方法に従って、検出装置によって検出されたサーマルサイクラーの試料ブロック内の試料チューブまたはウェルに対する開始サイクルと終了サイクルを決定する。１つの実施形態において、解析装置は、例えば、測定成長曲線を理論成長曲線と比較することによって、開始サイクルおよび終了サイクルも調整する。開始サイクルと終了サイクルを適応ベースラインを決定するために用いるので、それぞれのチューブまたはウェルに対して独自の開始サイクルおよび終了サイクルが割り当てられ、機器ドリフトおよびサンプリングのばらつきに起因する変動を回避することができる。それぞれのチューブ／ウェルに対する測定成長曲線と理論成長曲線の比較は、各チューブおよびウェルに対するベースラインを決定する手段として開始サイクルおよび終了サイクルを調整するために、（例えば、ｒｍｓｅを決定するためにプログラムを実行することにより）コンピュータの使用者によって行われる。それぞれのチューブ／ウェルに対する独自のベースラインを、増幅反応中に得られた生シグナルから差し引いて、チューブ／ウェル内の標的試料の量をより正確に反映する測定シグナルを得る。
【実施例】
【００９３】
実施例１
ＰＣＲは、ＰＣＲ試料プレートのウェル中のサーマルサイクリングブロック内に最大で９６個までの反応試料を保持するサーマルサイクラーを用いて実施した。９６ウェルすべてに同じ試料を入れた。適応ベースラインアルゴリズムを適用してＰＣＲ増幅を算出した。開始サイクルは、１〜８サイクルの間のサイクルに限定した。これらの試験において、すべてのウェルが同一のサイクルにおいてアンプリコンを検出し、増幅プロットはすべて一致していなければならない。図１１Ａおよび１１Ｂは、典型的な均一実行からの得られた結果である。図１１Ａは、開始サイクルを１〜８サイクルの間のサイクルに限定していない場合のデータ解析結果を示し、図１１Ｂは開始サイクルの決定に１〜８サイクルの限定を適用した場合の同じデータから得られた結果を示している。１〜８サイクルの限定は、データの広がりを小さくするだけでなく、図１１Ａの２試料の除去によって見られるように、人工的な疑似結果の発生を低減する。
【００９４】
本明細書中に記載および請求する本発明の精神と範囲から逸脱しない限りにおいて、当業者によって、ここに記載されているものの変形、変更および他の実施がなされるであろう。
【００９５】
特許および特許出願を含む上記で確認した全ての参考文献は、該参考文献が、本発明の１つ以上の実施形態の実施に重要であるかもしれない構造および／または装置の基礎を記述、列挙、提供するかこれらを可能にする範囲において、参照により本願に明確に組み込む。
【図面の簡単な説明】
【００９６】
【図１】本発明の１つの実施形態に従って、一定の開始サイクルおよび終了サイクルを用いたベースラインによって算出した、４試料に対するＰＣＲ増幅反応を示すグラフ。
【図２Ａ】本発明の１つの実施形態に従うサイクル数の関数としての実増幅シグナルのプロットを示すグラフ。
【図２Ｂ】本発明の１つの実施形態に従って、図２Ａのデータを用いてＰＣＲ反応の各サイクルにおけるシグナルの測定値と予測値の間の差異の二乗平均平方根誤差（「ｒｍｓｅ」）算出の結果を示したグラフ。
【図２Ｃ】本発明の１つの実施形態に従って、一定の開始サイクルおよび終了サイクルを用いたベースラインによって算出した、４試料に対するＰＣＲ増幅反応を示すグラフ。多くの試料のデータは、シグナルが検出可能な強度に上昇するまでは、０未満である。
【図２Ｄ】本発明の１つの実施形態に従って、図２Ｃで用いたものと同じ４試料に対する適応ベースラインを用いて算出したＰＣＲ増幅を示すグラフ。終点はそれぞれの試料に対して独立して選択されており、増幅が検出できるようになるまでのゼロ未満のデータ点はより少ない。
【図３】本発明の１つの態様に従う３５０〜８３０ｎｍの範囲の放射波長での蛍光シグナルの検出。
【図４】本発明の１つの態様に従う１つの試料チューブ中での４つまでの蛍光シグナルの検出。
【図５】本発明の１つの態様に従う最適化されたフィルタセットによって測定された様々な蛍光染料の、シグナル−ノイズ比。
【図６】本発明の１つの態様に従って用いられる典型的な染料に対する蛍光染料発光曲線。
【図７】本発明の１つの態様に従う、サイクル回数の関数としての蛍光シグナルとして示されるＰＣＲの進行に伴うリアルタイム増幅プロット。
【図８】本発明の１つの態様に従う、１つの増幅プロット中に見られる４つまでの蛍光シグナル。
【図９】本発明の１つの態様に従うベースライン減算を用いて作成した増幅プロット。
【図１０】ＱＰＣＲのために用いられ、本発明の１つの態様に従うベースライン減算を実施するために適応されるサーマルサイクラーと連絡するユーザー装置のインターフェース上に設けられるデータ解析オプション。
【図１１Ａ】本発明の１つの実施形態に従って、開始サイクルが１〜８の間のサイクルに限定されていない、ＰＣＲ増幅解析の結果を示すグラフ。
【図１１Ｂ】本発明の１つの実施形態に従って、開始サイクルが１〜８の間のサイクルに限定されている、ＰＣＲ増幅解析の結果を示すグラフ。
【図１２Ａ】本発明の１つの実施形態に従う解析装置の一例。
【図１２Ｂ】本発明の１つの実施形態に従う解析装置の別の例。

Claims

ＰＣＲ反応のサイクルを循環する試料ブロック内の複数の試料ウェルに対する適応ベースラインの算出法であって、前記方法は以下を含む。
（ａ）ウェル内で観察される実光学光学シグナルの強度を、該ウェルに対するサイクル数の関数としてプロットして、第１のプロットを得る工程と、
（ｂ）前記ウェルに対する開始サイクルと終了サイクルとを決定し、該開始サイクルと終了サイクルならびに両サイクルの間の点を直線に適合して、前記点と前記直線の間の最良適合を決定する工程と、
（ｃ）前記ウェルに対して得られた最良適合線を、該ウェルに対するベースラインとして用いる工程と、
（ｄ）前記複数のウェルのそれぞれに対して工程（ａ）〜（ｃ）を実施する工程。
前記開始サイクルと終了サイクルとを決定する工程（ｂ）は、サイクル２から始まる各サイクルに対して得られる実シグナルの強度の３点移動平均を得ることと、前記３点移動平均を、サイクル２から始まるサイクル数の関数としてプロットして第２のプロットを得ることとを含むことを特徴とする請求項１に記載の方法。
前記第１または第２のプロットからサイクルＸを同定する工程をさらに含み、前記サイクルＸは直前のサイクルからの勾配変化が最大となり、その直後に、それぞれ直前のサイクルと比較して増加勾配を有する少なくとも４つのサイクルを有することを特徴とする請求項２に記載の方法。
サイクルＸは、前記ＰＣＲ反応の最初の８サイクル以内で特定され、終了サイクルは前記サイクルＸの３サイクルまたは２サイクルまたは１サイクル前で特定されることを特徴とする請求項３に記載の方法。
サイクルＸは、ＰＣＲ反応の最初の８サイクル以内で特定され、サイクルＸは終了サイクルとして割り当てられることを特徴とする請求項３に記載の方法。
開始サイクルは、直前のサイクルとは異なる勾配傾向を示す第１のサイクルと、直前のサイクルの勾配よりも大きな正の勾配を有するサイクルと、初期勾配の１０％未満の勾配しか有しないサイクルとから成る群より選ばれることを特徴とする請求項３に記載の方法。
開始サイクルはサイクル１〜８の間のサイクルに限定されていることを特徴とする請求項３に記載の方法。
前記最良適合線は、式：

（式中、Ｎは前記開始サイクルと終了サイクルの間のサイクル数であり、Ｓは前記開始サイクルであり、Ｅは前記終了サイクルであり、γ_ｉは第ｉサイクルのシグナルであり、ｍおよびｂは前記最良適合線のパラメータである。）
で表される標準最小二乗平均平方根誤差アルゴリズムを用いて算出することを特徴とする請求項１に記載の方法。
前記ウェルの少なくとも１つは、鋳型核酸を含む可能性のある試料を含み、前記光学シグナルは、前記核酸に特異的な標識されたプローブに由来することを特徴とする請求項１に記載の方法。
前記ウェルの前記少なくとも１つに対して算出された前記適応ベースラインを、前記少なくとも１つのウェルに対してプロットされた第１のプロットから差し引いて、適応シグナルプロットを作成する工程（ｄ）をさらに含む請求項９に記載の方法。
前記適応シグナルプロットを、前記標準プロットがベースラインを除去するように調整されている、既知量の前記鋳型濃度を含むＰＣＲ反応から作成した第１の標準プロットと比較する工程をさらに含む請求項１０に記載の方法。
前記比較は、前記鋳型の量を決定するために用いられることを特徴とする請求項１１に記載の方法。
前記標識は蛍光標識であることを特徴とする請求項９に記載の方法。
前記試料は前記ウェル内のチューブ内にあることを特徴とする請求項９に記載の方法。
請求項１に記載の方法を実施するためのプログラム命令を含むコンピュータプログラム製品。
ＰＣＲ反応の適応ベースラインの算出のための装置であって、前記装置は、請求項１５に記載のコンピュータプログラム製品のプログラム命令を実行するためのメモリを含む解析装置を含んでいることを特徴とする装置。
前記試料ブロックを受容するためのサーマルサイクラーと、前記試料ウェルからの光学シグナルを収集するための入力装置と、解析装置によって取得または発生されたデータを表示するための出力装置とをさらに含むことを特徴とする請求項１６に記載の装置。
前記出力装置はコンピュータであることを特徴とする請求項１７に記載の装置。