JP6602374B2 - 試料の分析方法 - Google Patents

試料の分析方法 Download PDF

Info

Publication number
JP6602374B2
JP6602374B2 JP2017517690A JP2017517690A JP6602374B2 JP 6602374 B2 JP6602374 B2 JP 6602374B2 JP 2017517690 A JP2017517690 A JP 2017517690A JP 2017517690 A JP2017517690 A JP 2017517690A JP 6602374 B2 JP6602374 B2 JP 6602374B2
Authority
JP
Japan
Prior art keywords
cycle
baseline
threshold
amplification
data set
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2017517690A
Other languages
English (en)
Other versions
JP2018502552A (ja
Inventor
ユーン・チュン ジョン
チョ・リ ヨン
Original Assignee
シージーン アイエヌシー
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by シージーン アイエヌシー filed Critical シージーン アイエヌシー
Priority claimed from PCT/KR2015/010308 external-priority patent/WO2016052991A1/en
Publication of JP2018502552A publication Critical patent/JP2018502552A/ja
Application granted granted Critical
Publication of JP6602374B2 publication Critical patent/JP6602374B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6844Nucleic acid amplification reactions
    • C12Q1/686Polymerase chain reaction [PCR]
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6844Nucleic acid amplification reactions
    • C12Q1/6851Quantitative amplification
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • G16B20/20Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • G16B25/20Polymerase chain reaction [PCR]; Primer or probe design; Probe optimisation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2537/00Reactions characterised by the reaction format or use of a specific feature
    • C12Q2537/10Reactions characterised by the reaction format or use of a specific feature the purpose or use of
    • C12Q2537/165Mathematical modelling, e.g. logarithm, ratio
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2561/00Nucleic acid detection characterised by assay method
    • C12Q2561/113Real time assay

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Organic Chemistry (AREA)
  • Biotechnology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Genetics & Genomics (AREA)
  • Analytical Chemistry (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Evolutionary Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Microbiology (AREA)
  • Epidemiology (AREA)
  • Evolutionary Computation (AREA)
  • Bioethics (AREA)
  • Public Health (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Apparatus Associated With Microorganisms And Enzymes (AREA)

Description

本発明は、試料の分析方法に関する。特に、本発明は、試料を分析する方法及び増幅反応の生データセットを補正する方法に関する。
試料の分析は、多様な技術分野で非常に重要な技術である。試料の分析は、特定の特性の観点で試料を詳細に説明するか、描写するか、あるいは特徴を確認するために行う。
試料の分析は、生命工学分野でさらに重要な意味を有する。特に、試料の分析は、一般に分析物質(analyte)の存在または不存在、結合親和度、酵素活性、遺伝子発現水準及びアミノ酸やヌクレオチド配列のような特定特性に対する情報を提供するために行う。代表的に、試料の分析のために免疫分析法(immunoassay)及び遺伝子分析法(genetic analysis)が広く使用されてきて、生体試料(biosample)を分析する技術に係る米国特許第6,516,276号、第6,228,593号、第7,349,809号、第7,115,229号及び第6,816,790号に公開されている。
標的核酸増幅反応は、標的核酸分子を検出する大部分の技術に広く関連している。核酸増幅は、多様な増幅技術が提案されているほど、分子生物学分野の多様な方法のうち、最も核心的な技術である。核酸増幅技術として最も広く使用される技術は、重合酵素連鎖反応(Polynucleotide Chain Reaction: PCR)であって、この技術は、二本鎖DNAの変性、DNA鋳型へのオリゴヌクレオチドプライマーのアニーリング及びDNA重合酵素によるプライマー伸長の繰り返されたサイクル過程に基づいている(Mullisら, 米国特許第4,683,195号, 第4,683,202号及び第4,800,159号; Saiki et al., (1985) Science 230, 1350−1354)。
リアルタイムPCRは、試料から標的核酸分子を実時間で検出するためのPCRベースの技術である[Logan J et al., (2009). Real−Time PCR: Current Technology and Applications. Caister Academic Press]。特定標的核酸分子を検出するために、標的核酸の量に比例して検出可能な蛍光信号を放出する信号発生手段を利用する。検出可能な蛍光信号の放出は、例えば、二本鎖DNA間に結合したら信号を放出するインターカレーター(intercalator)を使用するか、蛍光レポーター分子及びクエンチャー分子を有するオリゴヌクレオチドを使用する方法がある。標的核酸の量に比例する蛍光信号が各増幅サイクル毎に検出されて、これを増幅サイクルに対してプロットし、増幅曲線(amplification curve)または増幅プロファイル曲線(amplification profile curve)を得るようになる。
一般に、リアルタイムPCRによる増幅曲線は、ベースライン領域(baseline region)、指数領域(exponential phase)、直線領域(linear phase)及び渋滞領域(plateau phase)に分かれる。指数領域は、PCR増幅産物の増加に比例し、放出される蛍光信号が増加する領域である。直線領域では、蛍光信号増加が実質的に減少し、実質的に線形方式で行動するようになり、渋滞領域は、PCR増幅産物の増加及び蛍光信号の放出が飽和状態に至り、それ以上蛍光信号の増加が現れない領域を意味する。
ベースライン領域は、PCR反応初期サイクルの間に、蛍光信号がほとんど変化せず一定に維持される領域を意味する。この領域では、PCR産物の水準が、検出できるほど十分ではないため、この領域の蛍光信号の大部分は、反応試料自体の蛍光信号及び測定システム自体の蛍光信号を含む背景信号(background signal)に起因する。
リアルタイムPCRデータの正確且つ再現性のある分析のために、得られた増幅曲線の補正(correction)または正規化(normalization)を行わなければならない。増幅曲線の補正は、ベースライン領域を決定し、決定されたベースライン領域における背景信号(background signal)を除去することにより行うことができる。
背景信号は、PCR反応時、反応条件及び環境の変化を反映するため、反応毎に変わって、これにより、試料内の核酸分子の量とは無関係なベースライン移動(baseline drift)が頻繁に起こる。このようなベースライン移動は、反応間の増幅曲線の比較を難しくして、標的核酸検出において、偽陽性または偽陰性結果を招来しやすい。したがって、PCRデータ分析において、このようなベースライン領域の適切な設定及び設定されたベースラインに基づいたPCR実験データの修正が必要である。
増幅曲線を補正するための従来の一つの方法では、PCR初期サイクルの任意のサイクル領域(例えば、3〜15サイクル)をベースラインと決定して、他の方法では、実験により増幅曲線が決定された後、増幅信号が十分増加する前にサイクルを決定して、ベースラインを設定する。米国特許第8,219,324号は、増幅曲線の2次導関数(second derivative)を求めて、一定な特徴を有するデータポイントをベースラインのエンドポイントサイクルとし、ベースライン領域を設定する方法を開示する。
このような従来の方法は、幾つかの重要な短所を有する。
前述の従来の一つの方法において、ベースラインを一つの初期任意サイクル領域に決定する場合、反応毎に相異なる背景信号の変化を補正することはできるが、この方法は、ベースライン移動(baseline drift)を補正するものではない。また、増幅領域の始まりは、初期試料に存在する標的核酸の量によって変わるため、ベースライン領域が予め決定されている場合、これを多様なサンプルに一律的に適用することは困難である。また、ベースライン領域を実験者が任意で決定する場合、同一増幅曲線といっても、分析する実験者によってベースライン領域設定が相異なる場合があり、これにより、実際増幅産物の量が変わってしまうため、信頼できる結果を導出することが困難である。
または、米国特許第8,219,324号に開示された方法によると、ベースラインは、複雑なアルゴリズムにより決定されるが、この場合、前記アルゴリズムには、明確に定義されていない多くのパラメータが必要になり、最適化が難しくなる。
閾値を使用する多様な試料分析方法において、ノイズ信号または非典型的なパターンの信号(例えば、陰の傾きパターン)の発生により、偽陽性または偽陰性結果が引き起こされる。このような分析誤りは、従来の閾値を使用する方法によっては除去することが難しい。
したがって、当業界では、各試料(またはPCR反応)に対し、新しい方式で閾値を設定するか、正確なベースライン領域を設定することにより、試料分析方法を改善するための新しい接近法(例えば、増幅曲線の補正)の開発が切実に必要であり、これは、より正確で且つ信頼性のある分析結果に寄与する。
本明細書全体にかけて多数の引用文献及び特許文献が参照され、その引用が表示されている。引用された文献及び特許の開示内容は、その全体が本明細書に参照として取り込まれ、本発明の属する技術分野の水準及び本発明の内容がより明確に説明される。
本発明の発明者は、信号発生反応を通じて収得したデータセットの加工を通じて、さらに正確で且つ信頼性のある信号発生反応の結果を得て、これにより、さらに正確で且つ信頼性のある試料の分析結果を提供するために鋭意研究した。その結果、信号発生反応による信号の有意性を決定するに邪魔になる要素、または信号発生反応による実際信号の増加を示すものではない不正確な信号を効果的に除去するために、サイクルの少なくとも二つのサイクルの閾値が互いに異なっている可変閾値が、信号発生反応のサイクルに適用できることを見出した。また、本発明が増幅反応の生データセットの補正に適用できることを見出した。
したがって、本発明の目的は、可変閾値を使用した試料分析方法を提供することである。
本発明の他の目的は、信号発生手段を利用した増幅反応の生データセットを補正する方法を提供することである。
本発明のまた他の目的は、試料分析方法を実行させるためのプロセッサを具現する指示を含む、コンピュータ読み取り可能な記録媒体を提供することである。
本発明のまた他の目的は、信号発生手段を利用した増幅反応の生データセットを補正する方法を実行するためのプロセッサを具現する指示を含む、コンピュータ読み取り可能な記録媒体を提供することである。
本発明のまた他の目的は、試料分析のための装置を提供することである。
本発明のまた他の目的は、信号発生手段を利用した増幅反応の生データセットを補正する装置を提供することである。
本発明のまた他の目的は、試料分析方法を実行させるためのプロセッサを具現する、コンピュータ読み取り可能な記録媒体に格納されたコンピュータプログラムを提供することである。
本発明のまた他の目的は、信号発生手段を利用した増幅反応の生データセットを補正する方法を実行するためのプロセッサを具現する、コンピュータ読み取り可能な記録媒体に格納されたコンピュータプログラムを提供することである。
本発明の他の目的及び利点は、発明の詳細な説明、請求の範囲及び図面により、さらに明確にされる。
本明細書の発明の詳細な説明において、本発明と直接的な関連がなく、技術分野によく知られた一般技術及び知識に関する説明は、詳細な説明が本発明をさらによく説明するために省かれる。また、以下説明されたセクション間共通する内容は、明細書の複雑性を避けるために省かれる。
I.可変閾値を利用した試料の分析
本発明の一様態によると、本発明は、(a)サイクルにおける信号値を提供するために、前記試料を使用する信号発生反応のそれぞれのサイクルにおける信号値を収得する段階と、
(b)それぞれのサイクルに閾値を適用し、複数の閾値をサイクルに適用する段階(前記サイクルの少なくとも二つのサイクルの閾値は、互いに異なっている)と、
(c)閾値のそれぞれにより定められる閾値基準を満足する一つ以上のサイクルを確認する段階と、
(d)前記段階(c)で確認されたサイクルを使用して試料を分析する段階と、
を含む、試料の分析方法を提供する。
本発明は、信号発生反応で得たデータセット(例えば、信号値とサイクル)を使用して、これを加工して試料を分析する方法に関する。
本発明の一具現例によると、試料の分析は、試料内分析物質の存在または不存在を決定することである。用語‘試料内分析物質の存在または不存在の決定’は、試料内分析物質の存在または不存在を定量的、定性的に決定することを意味する。
本発明の発明者は、信号発生反応を通じて収得したデータセットの加工を通じて、さらに正確で且つ信頼性のある信号発生反応の結果を得て、これにより、さらに正確で且つ信頼性のある試料の分析結果を提供するために鋭意研究した。その結果、信号発生反応による信号の有意性を決定するに邪魔になる要素、または信号発生反応による実際信号の増加を示すものではない不正確な信号を効果的に除去するために、サイクルの少なくとも二つのサイクルの閾値が互いに異なっている可変閾値が、信号発生反応のサイクルに適用できることを見出した。また、本発明が増幅反応の生データセットの補正に適用できることを見出した。
サイクルのうち、少なくとも二つのサイクルの閾値が互いに異なる値を有する閾値が信号発生反応のサイクルに適用される本方法は、知られていない。
本発明の特徴と利点を要約すると以下の通りである。
(a)試料分析のための本発明は、多重反応及び工程でよく現れる、誤り信号に基づいたサイクルの決定を防止し、これにより、試料分析のための情報をより正確に得ることができる。
(b)試料分析のための本発明において、閾値は、各サイクルに適用され、複数の閾値がサイクルに独立して適用され、これにより、試料の分析において、非正常的な信号の影響を除去する。従来の方法は、試料分析において、信号そのものを分析して、非正常的な信号を除去する。したがって、本発明の方法は、従来の方法とは異なるアルゴリズムで行われるものであり、したがって、従来の方法と共に使用することができ、これにより、試料分析の正確性を大幅に向上させることができる。
(c)本発明は、各試料のより正確なベースライン領域の設定を通じて増幅曲線を補正するため、増幅反応の結果は、より正確で、信頼性高く分析される。
(d)本発明は、簡潔な過程(アルゴリズム)によって増幅曲線を補正するため、分析対象および測定装備による最適化が非常に容易である。
図1Aは、本発明の試料分析方法の一具現例に対する流れ図である。 図1Bは、本発明の増幅曲線補正方法の一具現例に対する流れ図である。 図2は、信号発生手段としてTaqManプローブを使用して行った50サイクル以上リアルタイムPCRの結果を示す補正前増幅曲線を示す。リアルタイムPCRは、信号発生手段の存在下で行われて、増幅サイクル番号と測定された信号を含む生データセットを収得する。生データセットは、プロットされる。RFUは、相対的蛍光単位を表示する。 図3は、 図2の補正前増幅曲線の傾き曲線を示す。Y軸は、生データセットを利用して計算された各増幅サイクルにおける傾きを示す。傾きは、最小二乗法により計算された。Sは、ベースライン領域のスタートポイントサイクルであり、Eは、ベースライン領域のエンドポイントサイクルである。初期サイクルの一つがスタートポイントサイクル(S)と決定されて、ベースライン閾値と傾き曲線の一番目のクロスポイントサイクルがエンドポイントサイクル(E)と決定された。 図4Aは、図2の補正前増幅曲線(生データセット)のベースライン領域と、最小二乗法により収得した最適ライン関数により書かれた最適ライン(線形回帰線)のベースライン領域の拡大図である(Sは、ベースライン領域のスタートポイントサイクルであり、Eは、ベースライン領域のエンドポイントである)。 図4Bは、生データセットの測定された信号値から最適ライン関数の値を差し引いて収得した、補正された増幅曲線のベースライン領域の拡大図である(Sは、ベースライン領域のスタートポイントサイクルであり、Eは、ベースライン領域のエンドポイントである)。 図5は、図2の補正された増幅曲線を収得する過程を示している。図2の生データセットの蛍光信号強度から最適ライン(線形回帰線)関数値が差し引かれ、補正されたデータセットが収得されて、補正されたデータセットをプロットし、補正された増幅曲線を収得した。 図6は、本発明の試料分析用プログラムを備えたリアルタイムPCRシステムの一具現例を図式的に示している。 図7は、ベースライン領域のエンドポイントを決定するために、固定ベースライン閾値(‘300’または‘30’)を全体増幅サイクルに使用する場合、ベースライン領域決定に誤りが発生することを示す。 図8Aは、高濃度試料においてベースライン領域を決定する具現例を示す。エンドポイントサイクルは、傾き曲線とベースライン閾値とのクロスポイントサイクルのうち、最小ベースラインエンドポイントサイクル(MBEC)以後サイクルに決定された。 図8Bは、低濃度試料においてベースライン領域を決定する具現例を示す。エンドポイントサイクルは、傾き曲線とベースライン閾値とのクロスポイントサイクルのうち、最小ベースラインエンドポイントサイクル(MBEC)以後サイクルに決定された。 図8Cは、エンドポイントサイクルの決定に、MBECを使用するか、使用せずに決定されたベースライン領域を使用して、増幅曲線を補正した結果を示す。 図9は、サイクルによって変われるベースライン閾値が増幅曲線(または傾き曲線)に適用された多数の具現例を図式的に示している。太線は、ベースライン閾値を示す。BTCCは、ベースライン閾値変化サイクルを示す。 図10Aは、高濃度試料の傾き曲線にVBT(可変ベースライン閾値、Variable Baseline Threshold)を適用した結果を示している。エンドポイントサイクルは、傾き曲線とBTCC(ベースライン閾値変化サイクル)によって相異して適用されるベースライン閾値と傾き曲線間のクロスポイントサイクルに決定される。 図10Bは、低濃度試料の傾き曲線にVBT(可変ベースライン閾値、Variable Baseline Threshold)を適用した結果を示している。エンドポイントサイクルは、傾き曲線とBTCC(ベースライン閾値変化サイクル)によって相異して適用されるベースライン閾値と傾き曲線間のクロスポイントサイクルに決定される。 図10Cは、エンドポイントサイクル決定にVBT(可変ベースライン閾値)を使用するか、使用せずに求められたベースライン領域を使用して増幅曲線を補正した結果を示している。 図11Aは、C値を決定するために、増幅曲線にVST(可変信号閾値)を適用した結果を示している。FST(固定信号閾値)方法は、公知技術を意味する。 図11Bは、C値を決定するために、順に希釈されたFlu Aのgenomic RNA(10−3、10−4、10−5、10−6及び10−7希釈)の増幅曲線にVST(可変信号閾値)を適用した結果を示している。 図11Cは、C値を決定するために、順に希釈されたFlu Aのgenomic RNA(10−3、10−4、10−5、10−6及び10−7希釈)の増幅曲線にVST(可変信号閾値)を適用した結果を示している。
図1Aは、本発明の試料分析方法の一具現例に対する流れ図である。本発明をより詳細に説明すると、以下のようである。
段階(a):信号値収得段階(S10)
第一、サイクルにおける信号値を提供するために、前記試料を使用する信号発生反応のそれぞれのサイクルにおける信号値を収得する。
本発明は、信号発生反応で得たデータセット(例えば、信号値とサイクル)を使用して試料を分析する方法に関する。段階(a)は、(i)試料を利用した信号発生反応のサイクル及び(ii)サイクルにおける信号発生反応の信号値を含むデータセットを収得する段階とも説明でできる。
本発明で用語‘信号発生反応’は、試料内分析物質の存在に依存的に信号を発生できる反応を意味する。
信号発生反応は、信号の変化を伴う。
本発明の一具現例によると、信号発生反応は、信号増幅反応である。
本発明で用語‘信号’は、測定可能な産出物を意味する。
信号変化は、分析物質の存在または不存在を定量的または定性的に示す指標役割をすることができる。
有用な指標の例としては、蛍光強度、発光強度、化学発光強度、生発光強度、燐光強度、電化伝達、電圧、電流、電力、エネルギー、温度、粘度、光散乱、放射線強度、反射力、透過度、吸光度がある。最も広く使われる指標としては、蛍光強度がある。
本発明の一具現例によると、信号発生反応は、増幅曲線を提供する反応である。特に、増幅曲線は、信号増幅曲線である。
このような信号発生反応は、生物学的または化学的反応を含む。生物学的反応は、PCR、リアルタイムPCR、マイクロアレイ、インベーダーアッセイ、免疫アッセイ反応及び細菌培養分析のような遺伝的分析反応を含む。特に、信号発生反応を遺伝的分析反応を含む。化学的反応は、化学物質の生産、変化または分解を含む化学的分析方法を含む。
分析物質(analyte)は、核酸(例えば、DNAまたはRNA)、タンパク質、ペプチド、炭水化物、脂質、アミノ酸、生物学的化学物、ホルモン、抗体、抗原、代謝産物、細胞などである。択一的に、分析物質は、化学物質のような非生物学的物質である。
本発明の一具現例によると、分析物質は、標的核酸分子である。用語‘標的核酸分子’は、分析または検出対象核酸分子を意味する。
本発明の一具現例によると、信号発生反応は、分析物質の増幅を伴うかあるいは伴わない反応である。
特に、信号発生反応は、分析物質、具体的に標的核酸分子の増幅を伴う反応である。さらに具体的には、信号発生反応は、標的核酸分子の増幅を伴い、標的核酸分子の増幅によって信号が増加または減少(好ましくは、増加)できる反応である。
本発明で用語‘信号発生’は、信号の発生または消滅及び信号の増加または減少を全部含む。特に、用語‘信号発生’は、信号の増加を意味する。
本発明の一具現例によると、前記信号発生反応は、信号発生手段の存在下で行われる。
本発明で用語‘信号発生手段’は、分析物質(例えば、標的核酸分子)の存在を示す信号の発生に使用される物質を意味する。例えば、オリゴヌクレオチド、標識及び酵素である。択一的に、本発明で用語‘信号発生手段’は、信号発生のための物質を使用する方法を意味するものとして使用できる。
多様な種類の信号発生手段が当業者に知られている。信号発生手段は、標識そのもの及び標識の付いたオリゴヌクレオチドを全て含む。標識は、蛍光標識、発光標識、化学発光標識、電気化学的標識及び金属標識を含む。標識は、インターカレーティング染料(intercalating dye)のように、標識そのものが信号発生手段として使用できる。択一的に、単一標識または供与分子(donor molecule)及び受容分子(acceptor molecule)を含む相互作用的二重標識が、一つ以上のオリゴヌクレオチドに結合された状態で使用できる。
信号発生手段は、核酸分解酵素(例えば、5’−ヌクレアーゼまたは/及び3’−ヌクレアーゼ)のような信号を発生させるための要素を追加的に含むことができる。
本発明の方法は、標的核酸分子の存在または不存在を決定するために使用されるため、信号発生反応は、当業界に知られた多様な方法によって行うことができる。代表的な例としては、TaqManTM probe方法(米国第5,210,015号)、Molecular Beacon方法(Tyagiら, Nature Biotechnology, 14 (3):303(1996)),スコルピオン(Scorpion)方法(Whitcombeら, Nature Biotechnology 17:804−807(1999)),サンライズ(Sunrise or Amplifluor)方法(Nazarenkoら, Nucleic Acids Research, 25(12):2516−2521(1997),及び米国特許第6,117,635号),Lux方法(米国特許第7,537,886号),CPT(Duck Pら, Biotechniques, 9:142−148(1990)),LNA方法(米国特許第6,977,295号),Plexor方法(Sherrill CBら, Journal of the American Chemical Society, 126:4550−4556(2004)),HybeaconsTM(D. J. Frenchら, Molecular and Cellular Probes(2001) 13, 363−374及び米国特許第7,348,141号),Dual−labeled, self−quenched probe(米国特許第5,876,930号),Hybridization probe(Bernard PSら, Clin Chem 2000, 46, 147−148),PTOCE(PTO cleavage and extension)方法(WO 2012/096523),PCE−SH(PTO Cleavage and Extension−Dependent Signaling Oligonucleotide Hybridization)方法(WO 2013/115442),PCE−NH(PTO Cleavage and Extension−Dependent Non−Hybridization)方法(PCT/KR2013/012312)、CER方法(WO 2011/037306)がある。
本発明で用語‘増幅’または‘増幅反応’は、信号を増加または減少させるための反応を意味する。前記信号の増加または減少は、信号発生手段により起こる。
本発明の一具現例によると、信号発生手段による信号は、分析物質(例えば、標的核酸分子)の存在に依存的に発生して、前記信号の強度は、増幅反応によって増加または減少する。
本発明の一具現例によると、増幅反応は、分析物質(例えば、標的核酸分子)の存在に依存的に信号発生手段による信号を増幅するための反応を意味する。
本発明の一具現例によると、増幅曲線は増幅反応によって得られる。
本発明において、用語‘サイクル’は、状態の変化を伴う複数の測定において状態の変化の単位を意味する。状態の変化は、例えば温度の変化、反応時間の変化、反応回数の変化、濃度の変化、pHの変化および/または測定対象(例えば、標的核酸分子)の複製回数の変化を含む。したがって、サイクルは、時間またはプロセスサイクル、単位運営サイクル及び複製サイクルでありえる。
例えば、酵素の基質分解能が基質の濃度によって分析される場合、酵素の分解能に対する複数の測定が様々な基質濃度で行われる。この場合、基質濃度の増加が状態の変化に該当し、その増加の単位がサイクルに該当する。
他の例として、等温増幅反応は、反応時間に応じて、同一温度条件で複数の測定が行われ、この場合、反応時間が状態の変化に該当し、反応時間単位がサイクルに該当する。
特に、一連の反応の繰り返しまたは一定時間間隔を有する反応の繰り返しの場合、用語‘サイクル’は、前記繰り返しの単位を意味する。
例えば、重合酵素連鎖反応(PCR)において、サイクルは、標的分子の変性、標的分子とプライマーのアニーリング(混成化)及びプライマーの伸長を含む反応単位を意味する。この場合、反応の繰り返しの増加が状態の変化に該当し、前記繰り返し単位がサイクルに該当する。
本発明の一具現例によると、標的核酸分子が試料に存在する場合、測定される信号値(例えば、信号強度)は、増幅反応のサイクルの増加によって増加または減少する。
本発明の一具現例によると、標的核酸分子の存在を示す信号の増幅のための増幅反応は、信号が標的核酸分子の増幅と同時に増幅される方法(例えば、リアルタイムPCR)で行われる。択一的に、前記増幅反応は、信号が標的核酸分子の増幅を伴わずに増幅される方法[例えば、CPT方法(Duck P、et al.、Biotechniques、9:142−148(1990))、Invaderアッセイ(米国特許第6,358,691号及び第6,194,149号)]で行うことができる。
標的核酸分子を増幅する様々な方法が知られている。これらに限定されないが、例えば、重合酵素連鎖反応(PCR、polymerase chain reaction)、リガーアゼ連鎖反応(LCR、ligase chain reaction、Wiedmann Mら、’Ligase chain reaction(LCR)−overview and applications. ‘PCR Methods and Applications、1994 Feb; 3(4):S51−64)、GLCR(gap filling LCR、WO 90/01069、EP 439182、およびWO 93/00447)、Q−beta(Q−beta replicase amplification、Cahill P 、et al.、Clin Chem。、37(9):1482−5(1991)、米国特許第5,556,751号)、鎖置換増幅(SDA、strand displacement amplification、GT Walkerら、Nucleic Acids Res.20(7) :16911696(1992)、EP 497272)、NASBA(nucleic acid sequence−based amplification、Compton、J. Nature 350(6313):912(1991))、転写媒介増幅(TMA、Transcription−Mediated Amplification、Hofmann WPなど、 J Clin Virol.32(4):289−93(2005);米国特許第5,888,779号)またはローリングサークル増幅(RCA、Rolling Circle Amplification、Hutchison CAなど、Proc。Natl Acad。Sci。USA。102:1733217336(2005))である。
本発明の一具現例によると、本発明で使用する標識は、蛍光標識であり、より具体的に、単一蛍光標識または蛍光レポーター分子とクエンチャー分子を含む相互作用的二重標識である。本発明の一具現例によると、本発明で使用する増幅反応は、標的核酸分子の増幅を伴って信号を増幅させる。本発明の一具現例によると、本発明の増幅反応は、PCR方法によって行われる。
信号発生反応は、試料分析のためのデータセット(信号値及びサイクル値)を提供する。
本発明において、用語‘信号値’は、信号発生反応のサイクルで実際に測定された信号値(例えば、増幅反応による実際蛍光強度)またはその変形を意味する。前記変形は、測定された信号値(信号強度)の数学的に加工された値でする。測定された信号値の数学的に加工された値は、例えば、測定された信号値の対数値(logarithmic values)または微分値(derivatives)である。前記測定された信号値の微分値は、多重微分を含む。
本発明において、用語‘データポイント’は、サイクルと、前記サイクルにおける信号値を含む座標値を意味する。信号発生手段を用いた増幅反応によって得られたデータポイントは、直交座標系の座標値で表示することができる。前記直交座標系において、X軸は、増幅反応のサイクルを示し、Y軸は、該サイクルで前記信号発生手段からの信号値を示す(図2参照)。
本発明において、用語‘データセット’は、データポイントの集合を意味する。前記データセットは、生データセットと変形されたデータセットとを含む。
生データセットは、本発明の分析のための最初のデータセットである。前記生データセットは、試料分析のための信号発生反応(例えば、増幅反応)を介して直接収得したデータポイントのセットである。
例えば、本発明は、増幅反応の生データセットを補正するために使用される場合、生データセットは、増幅反応から直接収得したデータポイントの集合を含む(図2参照)。
変形されたデータセットは、生データセットが数学的に加工されたデータセットを含む。前記変形されたデータセットは、補正されたデータセット及び傾きデータセットを含む。補正されたデータセットは、生データセットの補正を介して収得したデータポイントのセットである。
本明細書では、生データセットと変形されたデータセットは、相対的な意味を有する。たとえば、生データセットは、あるデータの変形になる前のデータセットを意味し、変形されたデータセットは、データの変形後に収得したデータセットである。
本発明で使用されたデータセットは、信号発生反応で得られたデータポイントの部分または全部または補正されたデータポイントの全部又は一部を含むことができる。
本発明の一具現例によると、前記信号発生反応は、標的核酸分子の増幅を伴う反応である。より具体的に、前記標的核酸分子の増幅を伴う反応は、リアルタイム重合酵素連鎖反応である(リアルタイムPCR)。
本発明の一具現例によると、前記信号値は、信号発生反応から発生した信号値であるか、または信号発生反応から発生した信号の数学的に変形された値である。
本発明の一具現例によると、前記信号発生反応は、リアルタイムPCRであり、前記信号値は、リアルタイムPCRから発生された信号の数学的に変形された値であって、前記数学的に変形された値は、信号値をサイクルに対して微分した値である(図3参照)。前記サイクルに対して微分された信号値は、生データセットの導関数を含む。
図2は、信号発生反応として、リアルタイムPCRを通じて収得したデータセットの具体的な例を示しており、リアルタイムPCRの増幅曲線に該当する。図2に示されたデータセットは、信号発生反応から直接的に収得した生データセットである。生データセットは、リアルタイムPCRの増幅サイクル及び前記増幅サイクルで測定された信号強度(例えば、RFU)を含む。
図3は、生データセットの変形の一つを示しているが、これは、増幅サイクルで計算された傾きに対するデータポイントを含む。図3の曲線は、図2の生データセットのに導関数に該当する。
増幅反応を示す増幅曲線は、信号値を増幅サイクルに対してプロットして得ることができる。本明細書の増幅曲線は、データセットをプロットして得た曲線をいう。
補正前増幅曲線は、測定された信号値またはその変形に基づいた、各増幅サイクルで測定された信号値またはその変形された値を示す補正前増幅曲線を意味する。前記補正前増幅曲線は、増幅サイクルに対して測定された信号強度をプロットして得られる。特に、補正前増幅曲線は、生データセットをプロットして得ることができる。
補正後増幅曲線は、補正前増幅曲線をもとに補正された増幅曲線を意味する。前記補正された増幅曲線は、補正されたデータセットをプロットして得ることができる。
本発明において、用語‘標的核酸’または‘標的核酸分子’は、検出または定量する関心核酸分子を意味する。前記標的核酸分子は、二本鎖だけでなく、一本鎖の配列をすべて含む。前記標的核酸分子は、反応で新たに生成された配列だけでなく、初期から核酸試料中に存在する配列を含む。
前記標的核酸分子は、あらゆるDNA(gDNA、cDNA)、RNA分子またはこれらの交雑体(chimera nucleic acid)を全て含む。前記核酸分子は、二本鎖形態または一本鎖形態である。出発物質として核酸が二本鎖である場合、二本の鎖を、一本鎖または一部一本鎖形態に変形することが好ましい。鎖を分離する既知の方法は、、これに限定されず、例えば、熱処理、アルカリ、ホルムアミド(formamide)、尿素(urea)及びグリオキサール(glyoxal)処理、酵素的処理(例えば、ヘリカーゼ活性)、および結合タンパク質がある。例えば、鎖分離は80℃乃至105℃の熱処理を通じて行うことができる。このような処理を実現する一般的な方法は、Joseph Sambrookら、Molecular Cloning、A Laboratory Manual、Cold Spring Harbor Laboratory Press、Cold Spring Harbor、N.Y.(2001)で提供される。
前記標的核酸分子は、自然に発生した原核、真核(例えば、原生動物(protozoans)と寄生動物(parasites)、カビ(fungi)、酵母(yeast)、高等植物、下等動物、哺乳類や人間を含む高等動物)、ウイルス(例えば、ヘルペスウイルス(Herpes viruses)、HIV、インフルエンザウイルス(influenza virus)、エプスタイン−バールウイルス(Epstein−Barr virus)、肝炎ウイルス(hepatitis virus)、ポリオウイルス(polio virus)など)またはウイロイド(viroid)の核酸を含む。また、核酸分子は、組換え的に生産または化学的に合成された、あるいは合成できる核酸分子も含む。したがって、前記核酸配列は、自然から発見されるか、発見されないものである。標的核酸分子は、知られているか知られていない配列をすべて含む。
本発明の用語’試料’は、細胞、組織、生物由来の流体(fluid)または本発明によって評価できるその他の中間物質を意味し、前記試料は、例えば、ウイルス、バクテリア、組織、細胞、血液、血清、血漿、リンパ、喀痰、スワップ(swap)、吸引液(aspirate)、気管支肺胞洗浄液(bronchoalveolar lavage fluid)、乳、尿、糞便、眼球液、唾液、精液、脳抽出物、脊髄液、盲腸、脾臓および扁桃腺抽出液、羊膜液、腹水液および非生物学的試料(例えば、食品または水)である。また、前記試料は、化学反応の液相、固相の物質である。さらに、前記試料は、生物学的原料から抽出した天然核酸と合成された核酸分子の両方を含む。
段階(b):各サイクルに閾値適用段階(S20)
閾値がそれぞれのサイクルに適用され、複数の閾値がサイクルに適用される。前記サイクルの少なくとも二つのサイクルの閾値は、互いに異なっている。つまり、複数の閾値は、全体がまたは一部が異なる値を有する。
本発明の最も目立つ特徴は、複数の閾値を信号発生反応のサイクルに適用し、この際、前記サイクルの少なくとも二つのサイクルの閾値は、互いに異なっていることである。
各サイクルには、一つの個別閾値が割り当てられる。例えば、信号発生反応のサイクルが30個である場合、30個の閾値が個別に割り当てられる。前記割り当てられた閾値は、同一か、相異なっている。本発明の最も顕著な特徴は、前記割り当てられた閾値の少なくとも二つは、互いに異なっていることである。
閾値の適用は、閾値によって決定される閾値基準を満足するデータポイントを選択するために行われる。伝統的に、信号発生反応における信号値を評価するために、一つの閾値が採択される。つまり、従来の技術は、信号発生反応の信号値を評価するために、すべてのサイクルに渡って同一な閾値を使用する固定閾値方法を導入した。
これとは異なって、本発明は、サイクルの少なくとも二つのサイクルの閾値は、互いに異なる可変閾値を使用して、これにより、試料を分析する。
サイクルの少なくとも二つのサイクルの閾値は、互いに異なっている。いわば、複数の閾値は、全体が異なっているか、あるいは、一部が互いに異なる値を有する。
サイクルに対する閾値をプロットして得られたグラフをTCグラフ(threshold cycle graph)という。このTCグラフは、閾値セットをプロットして得られる。前記閾値セットは、閾値ポイントの集合である。閾値ポイントとは、サイクルと前記サイクルにおける閾値を含む座標値をいう。
ベースラインを求めるために、データセットに適用される閾値は、ベースライン閾値と命名され、サイクルに対してベースライン閾値をプロットして得られたグラフをBTグラフという。
本発明の方法が増幅反応の生データセットを補正するために使用される場合(eg、baselining)、TCグラフの一種のBTグラフは、ベースライン閾値をサイクルに対してプロットして得られる(図9参照)。
本発明の一具現例によると、サイクルの少なくとも二つ以上のサイクルは、互いに異なる閾値を有し、したがって、試料分析のための情報をとびきり正確に得ることができる。本発明ではこれを‘可変閾値(VT)’方法と命名する。前記VT方法は、可変ベースライン閾値方法と可変信号閾値方法とを両方とも含む。
本発明の一具現例によると、前記閾値は、閾値−変化サイクル(threshold−changed cycle、TCC)を基準に、前記TCC以前のサイクル及びこれらに適用される閾値の集合によって形成された関数と、前記TCC以後のサイクル及びこれらに適用される閾値の集合によって形成された関数とが互いに異なっているように決定される。
本発明の一具現例によると、TCCには、TCC以前のサイクルの関数またはTCC以後のサイクルの関数のいずれも適用することができる。
本発明において、用語‘閾値−変化サイクル’(TCC)は、閾値変化のパターンがサイクルを超えながら変更される基準サイクルを意味する。特に、本発明において、用語‘閾値−変化サイクル’(TCC)は、閾値が変わる基準サイクルを意味する。前記TCCは、一つであるか、または複数である。本発明において、用語‘TCC以前のサイクル’とは、閾値−変化サイクルより前のサイクルを意味し、‘TCC以後のサイクル’とは、閾値−変化サイクルより後のサイクルを意味する。
TCC以前のサイクルとこれらのTCC以前のサイクルに適用される閾値によって形成された関数と、TCC以後のサイクルとこれらのTCC以後のサイクルに適用される閾値によって形成された関数とが互いに異なっていることは、TCC以前のサイクルとTCC以後のサイクルのTCグラフが互いに異なるパターンを示すことを意味する。この例示は、図9に示されており、セクションIIで後述する。
前記TCCは、一つの反応に一つまたはそれ以上設定することができる。
本発明の一具現例によると、前記TCCの数は1,2,3,4,5,6,7,8,9,10,20,20,40または50個である。前記TCCの数は、70,60,50,40または30以下である。具体的に、前記TCCの数は、1〜2または1〜3個である。
本発明の一具現例によると、サイクルは、少なくとも一つのTCC(閾値−変化サイクル)を基準に、少なくとも二つの異なるグループに分類される。同一なグループに分類されたサイクルは連続的であり、同一な閾値を有する。互いに隣接したグループに属するサイクルは、互いに異なる閾値を有する。したがって、互いに離れたグループに属するサイクルは、相異なる、あるいは同一な閾値を有する。TCCはTCC前のサイクルまたはTCC後のサイクルと同一な閾値を有することができる。TCCは、一つのデータセットにおいて、一つが設定可能であり、この場合、データセットのサイクルは、二つのグループに分類される。択一的に、TCCは、一つのデータセットにおいて、二つ以上設定可能であり、この場合、データセットのサイクルは、3つ以上のグループに分類される。
表現‘互いに隣接したグループに属するサイクル’の例は、以下の通りである。信号発生反応は、合計40個のサイクルを含み、グループ1は、1〜10サイクル、グループ2は、11〜20サイクル、グループ3は、21〜30サイクルであって、グループ4は、31〜40サイクルである。この場合、互いに隣接したグループは、グループ1とグループ2、グループ2とグループ3またはグループ3とグループ4である。互いに離れたグループは、グループ1とグループ4、グループ1とグループ3またはグループ2とグループ4である。
本発明の一具現例によると、前記互いに隣接したグループの関数は、互いに異なっており、前記互いに離れたグループの関数は、互いに同一であるか、異なっている。
本発明の一具現例によると、前記の段階(b)は、サイクルの少なくとも一つのサイクルに追加閾値を適用する段階をさらに含む。
本発明の一具現例によると、付加的な閾値セットがデータセットに適用される。
本発明の一具現例によると、閾値セットは、サイクルの少なくとも二つのサイクルに適用される閾値が互いに異なっている、可変閾値セットである。
本発明の一具現例によると、閾値セットは、サイクルに適用された閾値が同一な、固定閾値セットである。
本発明の一具現例によると、多重閾値セットは、可変閾値セット及び固定閾値セットで構成された群から選択される少なくとも二つの閾値セットを含む。
本発明の一具現例によると、多重閾値セットは、少なくとも一つの固定閾値セットを含む。
本発明の一具現例によると、多重閾値セットは、少なくとも一つの可変閾値セットを含む。
本発明の一具現例によると、多重閾値セットは、少なくとも一つの固定閾値セットと、少なくとも一つの可変閾値セットを含む。
本発明の一具現例によると、多重閾値セットのそれぞれは、対応する閾値基準を有する。
本発明の一具現例によると、多重閾値セットは、データセットの分析に同時に適用される。複数のサイクルが閾値基準を満たしていることが観察されると、これらのすべてが分析に使用されるか、またはその一部が分析に使用できる。
本発明の一具現例によると、多重閾値セットは、データセットの分析に順次適用される。最初に適用された閾値セットの閾値基準を満たしているサイクルがない場合、別の閾値セットが適用される。
本発明の一具現例によると、閾値セットは、データセットに適用され、追加閾値は、サイクルの少なくとも一つのサイクルに適用できる。
追加閾値の適用については、後述するセクションIIの(3)多重閾値セット方法の記載を参照して説明することができる。
段階(c):閾値基準を満足するサイクルの確認(S30)
各サイクルに閾値を適用した後、各閾値によって決定された閾値基準(criterion)を満足する1つ以上のサイクルが確認される。
本発明において、用語‘閾値基準’は、一定の特徴を有するサイクルを確認するための基準をいい、前記一定の基準は、各閾値によって決定される。
本発明の一具現例によると、閾値基準は、閾値と関心サイクルにおける信号値を含む、ある参照(reference)または基準(benchmark)である。
本発明の一具現例によると、前記閾値基準は、サイクルの信号値と該サイクルに適用された閾値間の大きさ関係である。
本発明の一具現例によると、前記閾値基準は、各サイクルの信号値と、各増幅サイクルの閾値との比較である。
特に、前記閾値基準は、信号値が閾値以上であるか、以下であるかによって定義される。
本発明の一具現例によると、前記閾値基準は、閾値と同一であるか、閾値より大きい信号値を有するものである。
たとえば、閾値が1〜10サイクルでは、5に設定されて、11〜20サイクルでは、2に設定される場合、閾値基準は、信号値が1〜10サイクルでは、5以上、11〜20サイクルでは、2以上と定義される。
段階(d):確認されたサイクルを使用して試料を分析する段階(S40)
試料は、段階(c)で確認されたサイクルを使用して分析される。
本発明の一具現例によると、前記試料の分析は、試料内の標的核酸分子の存在を決定することであり、前記閾値基準を満足する1つ以上のサイクルの確認は、増幅曲線のC値を決定することである。この場合、閾値基準は、閾値と同一な信号値を有することができる。確認されるサイクルの数は、一つでありえる。
本発明の方法が、標的核酸分子の存在を決定するために、リアルタイムPCRでC値を決定するに適用される場合、本発明の方法は、以下のような段階を含む;
(a)サイクルにおける信号値を提供するために、前記試料を使用するリアルタイムPCRのそれぞれのサイクルにおける信号値を収得する段階と、
(b)それぞれのサイクルに信号閾値を適用して、複数の信号閾値をサイクルに適用する段階(前記サイクルの少なくとも二つのサイクルの閾値は、互いに異なっている)と、
(c)信号閾値のそれぞれによって決定される閾値基準を満足する一つ以上のサイクルを確認する段階と、
(d)前記段階(c)で確認されたサイクルを使用して、リアルタイムPCRにおけるC値を決定する段階。
リアルタイムPCRにおけるC値決定に適用する例は、実施例3及び図11A乃至11Cに例示されている。前記閾値は、C値決定においては、信号閾値と命名される。C値を決定する本方法は、VST(可変信号閾値)方法と命名される。
実施例3及び図11A乃至11Cに記載されているように、可変信号閾値を使用する本方法は、初期ノイズ信号が発生される初期増幅サイクルのデータポイントが、標的核酸分子が存在すると決定される誤りを除去するすることができる。さらに、本方法は、後期増幅サイクルにおいて、信号値が増加するスタートポイントをより正確に決定することができ、これによって、C値の決定において、誤りを除去することができる。
本発明の一具現例によると、試料の分析は、試料から標的核酸分子の存在を決定することであり、閾値基準を満足する一つ以上のサイクルを確認することは、リアルタイムPCRの増幅曲線のベースライン領域のエンドポイントサイクルを決定することができる。この場合、閾値基準は、ベースライン閾値と同一な信号値を有するものである。確認されるサイクルの数は、一つでありえる。
本方法が、標的核酸分子の存在を決定するための増幅曲線のベースライン領域のエンドポイントサイクルを決定するために適用される場合、本方法は、以下の段階を含む;
(a)サイクルにおける信号値を提供するために、試料を使用するリアルタイムPCRのそれぞれのサイクルにおける信号値を収得し、(i)リアルタイムPCRの増幅サイクル、および(ii)前記増幅サイクルにおける信号値を含む生データセットを収得する段階と、
(b)前記生データセットを使用してベースライン領域のスタートポイントサイクルとエンドポイントサイクルを決定して、ベースライン領域を決定する段階(前記エンドポイントサイクルは、下記段階によって決定される)と、
(b1)ベースライン閾値を各増幅サイクルに適用し、複数のベースライン閾値を増幅サイクルに適用する段階(サイクルの少なくとも二つのサイクルのベースラインの閾値は、互いに異なっている)と、
(b2)それぞれのベースライン閾値によって決定されるベースライン閾値基準を満足する1つ以上のサイクルを確認する段階と、
(b3)前記段階(b2)で確認されたサイクを使用してベースライン領域のエンドポイントサイクルを決定する段階、
(c)前記決定されたベースライン領域内生データセットの二つ以上のデータポイントからベースライン領域の最適ライン関数を算出する段階と、
(d)前記生データセットの信号値から前記ベースライン領域の最適ライン関数によって計算された関数値を差し引き、補正されたデータセットを収得する段階(前記補正されたデータセットは、(i)前記リアルタイムPCRの増幅サイクル及び(ii)前記差し引いた結果を含む)。
さらに具体的に、前記段階(b2)における確認は、生データセットを使用し、それぞれの増幅サイクルに対して計算された傾きと、それぞれの増幅サイクルに対するベースライン閾値とを比較して行われる。
セクションIIに説明された方法がセクションIの代表的な実施例であるため、その共通する内容は、明細書の複雑性を避けるために省く。
本方法は、増加パターン(例えば、増幅反応による信号変化)及び減少パターンの信号変化を含む如何なるパターンの信号変化にも適用することができる。
II.増幅反応の生データセットの補正
本発明の他の様態によると、本発明は、以下の段階を含む信号発生手段の存在下で行われた増幅反応の生データセットを補正する方法を提供する;
(a)(i)前記増幅反応における増幅サイクル及び(ii)前記信号発生手段から得た前記増幅サイクルにおける信号値を含む生データセットを収得する段階と、
(b)前記生データセットからベースラインスタートポイントサイクルとエンドポイントサイクルを決定して、ベースライン領域を決定する段階と、
(c)前記決定されたベースライン領域内生データセットの二つ以上のデータポイントからベースライン領域の最適ライン関数を算出する段階と、
(d)前記生データセットの信号値から前記ベースライン領域の最適ライン関数により計算された関数値を差し引き、補正されたデータセットを収得する段階(前記補正されたデータセットは、(i)前記増幅反応の増幅サイクル及び(ii)前記差し引いた結果を含む)。
増幅反応の生データセットを補正する方法である本方法は、試料を分析する方法の特定具現例であるため、これらの間に共通する記載は、明細書の複雑性を避けるために省く。
本発明の発明者は、増幅反応を通じて収得した生データセットを補正して、これにより、増幅反応に対する、より正確で信頼性の高い分析結果を定性的、定量的に提供する新しい方法について鋭意研究した。その結果、増幅反応の生データセットを補正するための新しい方法を発見した。特に、本発明の新しい方法によって、ベースライン領域のエンドポイントサイクルが決定でき、本発明は、増幅反応のより正確で信頼性の高い分析結果を収得するに寄与することを発見した。
図1Bは、本発明の増幅曲線補正方法の一具現例に対する流れ図を示す。本発明をより詳細に説明すると、以下のようである。
段階(a):生データセットの収得(S110)
最初に、生データセットが収得される。前記生データセットは、(i)増幅反応の増幅サイクル及び(ii)信号発生手段から得られた増幅サイクルにおける信号値を含む。
生データセットは、信号発生手段の存在下で増幅反応を行って得られる。段階(a)は、択一的に、信号発生手段の存在下で増幅反応を行って、生データセットを収得する段階と表現することができる。
本発明の一具現例によると、信号発生手段による信号は、標的核酸分子の存在に依存的に発生し、その強度(intensity)は、増幅反応の進行によって増加または減少する。
本発明の一具現例によると、増幅反応は、標的核酸分子の存在に依存的に、信号発生手段の存在下で信号を増幅する反応を意味する。
特に、一連の反応の繰り返しまたは一定時間間隔を有する反応の繰り返しの場合、用語‘サイクル’は、前記の繰り返しの単位を意味する。
例えば、重合酵素連鎖反応(PCR)において、サイクルは、標的分子の変性、標的分子とプライマーのアニーリング(混成化)とプライマーの伸長を含む反応単位を意味する。この場合、反応の繰り返しの増加が状態変化に該当し、前記繰り返し単位がサイクルに当たる。他の例として、LAMP(Loop−mediated isothermal amplification)及びNASBA(Nucleic acid sequence−based amplification)のような等温核酸増幅反応の場合、サイクルは、時間間隔を意味する。
本発明の一具現例によると、標的核酸分子が試料内に存在する場合、測定された信号値(例えば、信号強度)は、増幅サイクル数の増加によって増加するか、減少する。
前記生データセットは、(i)増幅サイクル及び(ii)信号発生手段から得られた増幅サイクルにおける信号値を含む。
本発明の用語‘信号値’は、増幅サイクルで実際測定された信号値またはその変形を意味する。前記変形は、測定された信号値(信号強度)の数学的に加工された値を含むことができる。数学的に加工された値は、例えば、測定された信号値の対数値(logarithmic values)または微分値(derivatives)を含む。前記測定された信号値の微分値は、多重微分を含むことができる。
上述のように、生データセットを収得する段階(S110)において、増幅サイクル及び前記増幅サイクルにおける信号発生手段からの信号値を含むデータセットは、増幅反応の実行によって得られ、補正前の増幅曲線(最初の増幅曲線)を提供するためにプロットされる。
(i)増幅サイクル及び(ii)前記増幅サイクルにおける信号値を含む生データセットは、増幅反応によって得られ、図2に示したような補正前の増幅曲線を提供するためにプロットされる。図2においてRFUは、相対的蛍光単位(relative fluorescence unit)である。
上述のように、増幅曲線は、ベースライン領域(baseline region)、指数領域(exponential phase)、直線領域(linear phase)と渋滞領域(plateau phase)に区分される。
ベースライン領域では、初期増幅サイクルの間に蛍光信号の変化がほとんどない。指数領域は、増幅産物の増加に比例して、蛍光信号が増加することを示す。直線領域では、蛍光信号の増加が実質的に減少し、実質的に線形に現れる。渋滞領域では、増幅産物と蛍光信号の飽和により、蛍光信号の増加はほとんどない。
ベースライン領域の信号の大部分を占める背景信号が、試料内標的核酸分子の量と無関係なベースライン移動(drift)の原因となるため、ベースライン領域の決定及び増幅曲線の補正がなされる必要がある。
段階(b):ベースライン領域の決定(S120)
次に、生データセットを用いたベースライン領域のスタートポイントサイクル及びエンドポイントサイクルの決定によって、ベースライン領域が決定される。
ベースライン領域決定に関する‘生データセットの利用’は、前記生データセットの直接的、間接的使用をすべて含む。生データセットの間接的利用は、生データセットの変形されたデータセットの利用を含む。
本発明の一具現例によると、スタートポイントサイクル及びエンドポイントサイクルは、生データセットから直接決定されるか、あるいは生データセットが数学的に加工されたデータセットから決定される。
例えば、スタートポイントサイクルは、一定の値以上の信号値を有する最初のサイクルに決定する方法で、生データセットから直接的に決定することができる。択一的に、エンドポイントサイクルは、一定の値以上の傾き値を有する最初のサイクルに決定する方法で、数学的に加工されたデータセットから決定することができる。
本発明において、用語‘スタートポイントサイクル’は、ベースライン領域の開始に当たるサイクルを意味する。
ベースラインスタートポイントサイクル(S)は、ユーザーによって任意に決定することができる。一般には、スタートポイントサイクルは、増幅反応の初期の典型的な変動形態が現れた後のサイクルに決定する。例えば、スタートポイントサイクルは、1〜10サイクルから決定でき、例えば2〜10、2〜8、2〜6または2〜4サイクルから決定できる。
択一的に、ベースライン領域のスタートポイントサイクル(S)は、一定条件を満足するサイクルを考慮して決定することができる。
例えば、スタートポイントサイクルは、以前サイクルと異なる傾き傾向を示す一番目のサイクルに決定することができる。前記以前サイクルと異なる傾き傾向を示すサイクルは、例えば、以前サイクルより大きい傾きを有するが、初期サイクルの傾きより10%以下の傾きを有するサイクル、以前サイクルの傾きが負のとき正の傾きを示すサイクルまたは以前サイクルの傾きが正のとき、負の傾きを有するサイクルである。
本発明の一具現例によると、各サイクルにおける信号値の変化率が計算され、これをベースラインのスタートポイントサイクルまたはエンドポイントサイクルの決定に利用することができる。本明細書で特に言及しない限り、用語‘傾き’は、選択されたサイクルで信号値の変化率を意味する。本発明の一具現例によると、傾きは、サイクルに対してプロットされ、傾き曲線を提供する。
本発明において、用語‘エンドポイントサイクル’は、ベースライン領域の終了に当たるサイクルを意味する。ベースラインエンドポイントサイクル(E)がベースライン領域の終わりを決定するため、ベースラインエンドポイントサイクルは、信号増幅が発生する前のサイクルで決定される。
前記ベースライン領域のエンドポイントサイクルは、生データセットまたはその変形されたデータセットから決定することができる。
ベースライン領域のエンドポイントサイクルは、さまざまな方法で決定できる。
たとえば、エンドポイントサイクルは、2次微分値が最大となるサイクルに決定することができる。または、データセットの傾き曲線の特徴(位置及び大きさ)がエンドポイントサイクルの決定のために分析される。例えば、傾き曲線は、指数領域のピークが最も大きいはずなので、ベースライン領域のエンドポイントサイクルは、指数領域のピークが始まるサイクルに決定できる。
また、エンドポイントサイクルは、以前サイクルと比較し、変動係数(coefficient of variation)が急激に増加することを示す一番目のサイクルで決定することができる。または、エンドポイントサイクルは、予め指定した値以上の変動係数を有するサイクルで決定することができる。前記変動係数は、標準偏差を算術平均で割ったものとして定義される。前記変動係数は、変動係数が算出されるサイクルとその一定数前、後のサイクルを選択して、これらの選択されたサイクルの標準偏差と算術平均を計算して変動係数を求めることができる。前記選択されたサイクルの一定数は、1,2,3,4または5であり、好ましくは、1または2である。
本発明の一具現例によると、前記エンドポイントサイクルは、前記スタートポイントサイクルと、データセットで最大の傾きを有するデータポイントのサイクルとの間で決定される。
本発明の一具現例によると、前記エンドポイントサイクルは、超過サイクル(not−less−than cycle)の中から選択することができる。または、前記エンドポイントサイクルは、超過サイクルの中で選択されたサイクルを、エンドポイントサイクルを決定するための基準サイクルとして決定することができる。
本発明において、用語‘超過サイクル’(not−less−than cycle(s))は、ベースライン閾値と同一またはそれ以上の値を有するデータポイントのサイクルを意味する。つまり、超過サイクルは、ベースライン閾値以上の傾きを有するサイクルである。
増幅曲線が減少するパターンの場合、エンドポイントサイクルは、未満サイクル(not−more−than cycle)の中から選択されたサイクルによって決定できる。
本発明の一具現例によると、前記エンドポイントサイクルは、クロスポイントサイクルの中で選択されたサイクルによって決定することができる。または前記エンドポイントサイクルは、クロスポイントサイクルの中で選択されたサイクルをエンドポイントサイクルの決定のための基準サイクルとして決定することができる。
本発明の用語‘クロスポイントサイクル’(cross−point cycle)は、ベースライン閾値と同一な値を有するデータポイントのサイクルを意味する。
前記クロスポイントサイクルは、データポイントのサイクルの中で決定できる。
前記クロスポイントサイクルは、ベースライン閾値とデータポイントを使用して数学的に計算されたサイクルによって決定することができる。
本発明の一具現例によると、エンドポイントサイクルは、ベースライン閾値と同一な傾きを有するデータポイントのサイクル、ベースライン閾値を超過する傾きを有する一番目のデータポイントのサイクルまたはベースライン閾値未満の一次傾きを有する一番目のデータポイントのサイクルによって決定される。または、特定値の傾きを有するデータポイントのサイクルまたは特定値を超過する一番目のデータポイントのサイクルを基準サイクルとして使用し、エンドポイントサイクルを決定するための数学的方程式に代入する方法でエンドポイントサイクルを決定することができる。前記数学的方程式の例は、‘エンドポイントサイクル=基準サイクル−(1,2,3または4サイクル)’;‘エンドポイントサイクル=基準サイクル+(1,2,3または4サイクル)’;‘エンドポイントサイクル=[基準サイクル×0.9]’;及び‘エンドポイントサイクル=[基準サイクル−(ベースライン閾値×0.1)]’を含む。[X]は、X値と同一か小さい最大の整数を表示する。
前記傾きと共に記載された特定値は、ベースライン閾値(BT)である。ベースライン閾値は、一般にベースラインを決定するために使用され、本発明では、ベースライン領域のエンドポイントサイクルを決定するために使用される。前記ベースライン閾値は、測定対象および/または測定機器に応じて予め決定(または入力)でき、またはユーザーによって任意に決定することができる。
本発明の一具現例によると、エンドポイントサイクルは、ベースライン閾値と傾き曲線の一番目のクロスポイントのサイクルまたは(前記)一番目のクロスポイントのデータポイントのサイクルに最も近いサイクルによって決定できる。
クロスポイントのサイクルまたはクロスポイントのデータポイントのサイクルをクロスポイントサイクルという。前記クロスポイントサイクルは、各サイクルに対してベースライン閾値をプロットして得られたグラフ(以下、ベースライン閾値グラフまたはBTグラフという)と傾き曲線が交互するポイントのサイクルである。
クロスポイントサイクルの数は、傾き曲線および/またはBTグラフの形に応じて、一つまたはそれ以上である。エンドポイントサイクルを決定するためのクロスポイントサイクルは、例えば、最初のクロスポイントサイクルまたは最後のクロスポイントサイクルのように、予め決定された特定クロスポイントサイクルによって決定できる。またはクロスポイントサイクルの数が二つ以上である場合、最も少ないサイクル数を有するクロスポイントサイクルがエンドポイントサイクルとして決定される。
クロスポイントのサイクルの値は、整数ではないことがある。実際の実験で、サイクルは整数で表示されるため、エンドポイントサイクルは、整数であることが望ましい。したがって、前記クロスポイントサイクルを超過する一番目の整数サイクル(integral cycle)または前記一番目の整数サイクルの1,2,3または4番目前または後のサイクルがエンドポイントサイクルとして決定できる。またはエンドポイントサイクルは、クロスポイントサイクル未満の最も大きい整数サイクルまたは前記最も大きい整数サイクルの1,2,3または4番目前または後のサイクルがエンドポイントサイクルとして決定される。
本発明の一具現例によると、ベースライン領域のエンドポイントサイクルは、ベースライン閾値と傾き曲線の一番目のクロスポイントのサイクル又は前記一番目のクロスポイントのデータポイントのサイクルの1,2,3または4サイクル前または後のサイクルによって決定される。
本発明の一具現例によると、ベースライン閾値は、増幅領域が現れる前の初期サイクルの間、背景信号によって干渉されないように設定できる。
本発明の一具現例によると、ベースライン閾値は、様々な試料の分析結果によって選択された適切な値に設定することができる。
図3は、2〜4サイクルのいずれかのサイクルがスタートポイントサイクルに決定されて、ベースライン閾値と傾き曲線間の一番目のクロスポイントサイクルがエンドポイントサイクルに決定されたベースライン領域を示す。
本発明では、前記決定されたスタートポイントサイクルとエンドポイントサイクル間の相対的距離が一定値以下である場合、前記相対的距離は、適切なベースライン領域を確保するために、追加に調整することができる。
前記相対的距離は、エンドポイントサイクルからスタートポイントサイクルを差し引いて計算することができる。前記追加調整が必要な相対的距離の一定値は、0,1,2,3,4,5,6,7,8,9または10サイクルであり、好ましくは、0,1,2,3または4サイクルである。
前記追加調整の方法は、スタートポイントサイクルまたはエンドポイントサイクル、またはその両方のサイクルを、前記相対的距離が一定値以上になるように任意に設定する方法である。または、前記追加調整方法は、スタートポイントサイクルまたはエンドポイントサイクルのいずれかを、上述のスタートポイントサイクルまたはエンドポイントサイクル決定方法によって決定し、前記決定されたスタートまたはエンドサイクルから前記相対的距離が一定値以上になるサイクルの範囲で、他の一つを、上述のスタートポイントサイクルまたはエンドポイントサイクル決定方法によって決定する方法で決定することができる。
本発明の一具現例によると、データポイントの傾きは、前記データポイントのサイクルにおける信号値の変化を示す。
上述のように、本発明の増幅反応は、サイクルに応じて信号減少を示す反応を含む。
前記傾きは、微分法など、様々な方法で計算することができる。
前記傾きは、特定サイクルのデータポイント及び前記特定サイクルの前および/または後のサイクルの少なくとも一つのデータポイントを使用して、最小二乗法またはLMS(least mean square)アルゴリズムなどにより計算することができる。
以下、線形回帰分析の代表的な例として、最小二乗法(least square method)を挙げて説明するが、本発明はこれに限定されない。
最小二乗法で傾きを計算するのに使用されるデータポイントの数は、二つ以上でありえる。例えば、最小二乗法で傾きを計算するために使用されるデータの数は、2,3,4,5,6,7,8,9,10,11,12,13,14または15個以下である。特に、前記データポイントの数は、2〜3、2〜15、3〜11、3〜9、3〜7、3〜5または5〜7個である。
最小二乗法で傾きを計算するのに使用されるデータポイントは、隣接したサイクルであるか、離れているサイクルである。
たとえば、傾きは、特定サイクルのデータポイント及び特定サイクル以前および/または以後のサイクルに対応する少なくとも一つのデータポイントを使用して、最小二乗法によって計算することができる。
他の例としては、前記傾きは、前記データセットのうち、特定サイクルに対応するデータポイント、特定サイクル以前のサイクルのデータポイント及び特定サイクル以後のサイクルのデータポイントを使用して、最小二乗法によって計算することができる。
最小二乗法によって傾き計算に使用されるデータポイントの数は、サイクルに応じて変わる。たとえば、データポイントの傾きは、互いに隣接した二つまたは三つのデータポイントを使用して、最小二乗法によって計算することができる。たとえば、一番目のサイクル以前にはサイクルが存在しないため、一番目のサイクルの傾きは、一番目のサイクルとその次のサイクルのデータポイントを使用して、最小二乗法によって計算することができる。最後のサイクル以後には、サイクルがないため、最後のサイクルの傾きは、最後のサイクルと隣接した直前のサイクルを使用して、最小二乗法によって計算することができる。他のサイクルに対しては、前記傾きは、該サイクルのデータポイント、前記該サイクルの直前のサイクルのデータポイント及び前記該サイクルの直後のサイクルのデータポイントを使用して、最小二乗法によって計算することができる。
本発明の一具現例によると、前記最小二乗法は、下記数学式1で表すことができる。
(前記数学式1において、
であり、
Iは、傾きを求めようとするデータポイントのサイクルであり、mは、I番目サイクルのデータポイントの傾きであって、xは、i番目サイクルのサイクル番号であり、yは、i番目サイクルで測定された信号値であって、nは、a+b+1であり、a及びbは、独立的に0以上10以下の整数(但し、aは、I未満であり、a+b+1は、2以上、生データセットのデータポイントの数以下、I+bは、生データセットのデータポイントの数以下)である。)
‘a+b+1’は、I番目のサイクルの傾きを計算するために使用されるデータポイントの数であり、LSMR(最小二乗法範囲; Least Squares Method Range)と称する。‘a’は、I番目のサイクルの傾きを計算するために使用されるデータポイントのうち、最小サイクルを計算するための値である。‘b’は、最大サイクルを計算するための値である。データポイントの数は、全体反応で収得したデータポイントの数を意味し、増幅曲線の最大サイクル値に当たる。
‘a’及び‘b’は、独立して、0乃至10の整数であり、好ましくは、1乃至5の整数、さらに好ましくは、1乃至3の整数である。
前記‘a’及び‘b’は、同一な値であることが有利であるが、これらは、測定対象、測定環境及びサイクルによって変わる。
前記‘a’及び‘b’は、I−aが1未満であるか、I+bが全体データポイントの数を超過する場合を除いて、反応の全体データポイントに同一に適用されることが有利である。しかし、信号値の変化側面で特別な属性または区間特性を示す特定データポイントまたは特定範囲のデータポイントの傾きは、他の前記‘a’及び‘b’を適用して計算することができる。
たとえある反応のすべてのデータポイントに一定な‘a’及び‘b’が適用される場合でも、I−aが1未満であるか、I+bが全体データポイントの数を超過するデータポイントの傾きを計算するためには、‘a’と‘b’の値は、他のデータポイントに適用されたものとは異なる値を有することができる。たとえば、I−aが1未満であるデータポイントでは、‘a’は‘I−a’が1になるように調整することができる。この場合、‘b’は、一定に維持されるか、‘a’の調整によって変わる。
I+bが全体データポイントの数を超過するデータポイントの場合、‘b’は、‘I+b’が全体データポイントの数と同一になるように調整される。この場合、‘a’は、一定に維持されるか、‘b’の調整によって変わる。
LSMRの値、‘a’及び‘b’は、測定対象および/または測定装備に応じて、あらかじめ決定(入力)することができ、またはユーザーによって任意に決定することができる。
図3は、数学式1で表現される最小二乗法によって計算された傾き曲線を示す。Y軸は、最小二乗法により各増幅サイクルで計算された蛍光信号強度(または相対的蛍光単位)の傾きを示す。
ベースライン領域決定のエンドポイントサイクル
上述のように、初期サイクルにおいて、背景信号による干渉を避けながら、ベースライン領域のエンドポイントサイクルは、生データセットから各増幅サイクルで計算された傾きと各増幅サイクルに対するベースライン閾値とを比較して決定することができる。
ベースラインの閾値は、ベースライン領域のエンドポイントサイクルを決定するために設定される。前記ベースライン閾値は、測定対象および/または測定装備に応じて、あらかじめ決定(入力)することができ、またはユーザーによって任意に決定することができる。
背景信号またはノイズによる増幅分析の誤りは、エンドポイントサイクルの誤った決定によるものでありえる。核酸増幅反応のような増幅反応において、非正常的な蛍光信号は、初期サイクルで頻繁に現れて、記録される。非正常的な蛍光信号の検出は、標的核酸分子の量を反映しない蛍光信号を検出することをいう。
ベースライン閾値が過度に低く設定される場合、非正常的な蛍光信号の傾き値がエンドポイントサイクルの決定に影響を与えるようになる。もしそう決定されたエンドポイントサイクルを使用してベースライン領域が決定され、増幅曲線が補正されると、補正された増幅曲線は、偽陽性結果を示すか、増幅産物の量を反映しない。もしエンドポイントサイクルの決定に非正常的な蛍光信号の傾き値が関与することを避けるために、ベースライン閾値が過度に高く設定される場合、ベースライン閾値と傾き曲線間の交差が、初期サイクルではなく、後のサイクルで起こるか、傾き曲線のピークが低い場合は、ベースライン閾値と傾き曲線間の交差が起こらないことがある。
本発明の一具現例によると、上述の問題を解決するために、特定サイクル以前のサイクルは、エンドポイントサイクルの決定のプロセスから除く。または初期背景信号がエンドポイントサイクルの決定に関与しないよう、ベースライン閾値が調整される。
本発明の一具現例によると、ベースライン閾値は、増幅領域が観察される前の初期のサイクルにおいて、背景信号によって干渉されないように設定することができる。
(1)特定サイクル以前のサイクルを除外する方法
本発明の一具現例によると、エンドポイントサイクルを決定するために、特定サイクル以前のサイクルを除外する方法において、前記特定サイクルは、最小ベースラインエンドポイントサイクル(MBEC、Minimum Baseline End−point Cycle)である。この方法において、エンドポイントサイクルは、MBEC以上のサイクルの中から決定される。MBEC適用によって、増幅反応初期サイクルの間の背景信号またはノイズ信号によって、エンドポイントサイクルが増幅領域の始まりより過度に離れて最初に決定されることを防止することができる。このようにエンドポイントサイクルの決定にMBECを設定する方法を‘MBEC方法’という。
前記MBECは、測定機器、個々の装備の特性、分析試料、試薬の影響による背景信号またはノイズ信号のパターンに応じて多様に設定することができる。MBECは、背景信号またはノイズ信号による誤りが防止される限り、設定サイクル範囲に制限がない。例えば、MBECは、1〜50サイクルの範囲で決定することができ、好ましくは、1〜10、1〜15、1〜20、1〜25、1〜30、1〜35、1〜40、5〜10、5〜15、5〜20、5〜25、5〜30サイクル、より好ましくは、5〜15サイクルの範囲で決定することができる。
図8A及び図8Bに示すように、傾き曲線の二つのクロスポイントサイクル(E、E)のうち、MBECが適用されると、MBECを超過するEがエンドポイントサイクルとして決定される。
また、MBECを超過するクロスポイントサイクルが二つ以上である場合、より小さいクロスポイントサイクルがエンドポイントサイクルに決定される。上述のように、MBECを超過する最も小さいクロスポイントサイクルの値が整数ではない場合、前記最も小さいクロスポイントサイクルを超過する一番目の整数サイクル値または前記一番目の整数サイクル値の1,2,3または4サイクル前または後のサイクルがエンドポイントサイクルとして決定される。または、前記最も小さいクロスポイントサイクル値未満の最も大きい整数サイクルまたは前記最も大きい整数サイクルの1,2,3または4サイクル前または後のサイクルがエンドポイントサイクルとして決定される。
本発明の一具現例によると、前記の段階(b)において、ベースライン領域のエンドポイントサイクルは、最小ベースラインエンドポイントサイクル(minimum baseline end−point cycle、MBEC)以上のサイクルで決定され、前記MBECは、増幅反応の前または後に決定される。
MBECが決定された後は、ベースライン領域のエンドポイントサイクルは、さまざまな方法で容易に決定することができる。例えば、クロスポイントサイクルまたは超過サイクルが先に確認され、これらがエンドポイントサイクルに決定できるかどうかを評価するために、MBECと比較する。または前記エンドポイントサイクルを単純にMBEC以後のサイクルの中から決定することができる。
もしMBEC以上のクロスポイントサイクルが存在しない場合は、傾き曲線の最後のサイクルがエンドポイントサイクルに決定される。
本発明の一具現例によると、前記ベースライン領域のエンドポイントサイクルは、以下の段階を含む方法によって決定することができる;
(i)それぞれの増幅サイクルに対して計算された傾きを収得する段階と、
(ii)前記傾きをそれぞれの増幅サイクルに対するベースライン閾値と比較し、ベースライン領域の候補エンドポイントサイクルを収得する段階と、
(iii)前記候補エンドポイントサイクルをMBECと比較する段階(前記候補エンドポイントサイクルがMBECより大きい場合、前記候補エンドポイントサイクルをエンドポイントサイクルに決定する)。
候補エンドポイントがMBEC未満の場合、前記候補は除去され、前記段階(i)および(ii)が、新しい候補エンドポイントを検索するために繰り返される。もしMBEC以上で候補エンドポイントがない場合、傾き曲線の最後のサイクルがエンドポイントサイクルに決定される。
本発明のMBEC方法は、好ましくない増幅反応の初期サイクルでエンドポイントサイクルが決定されることを防止し、結果的に、より正確に増幅曲線を補正することができる。
(2)可変ベースライン閾値方法
本発明の一具現例によると、ベースライン閾値(BT)は、初期背景信号がエンドポイントサイクルの決定に関与しないように調整される。
上述のように、生データセットを利用して、各増幅サイクルで計算された傾きと各増幅サイクルに対するベースライン閾値とを比較して、エンドポイントサイクルを決定するために、固定ベースライン閾値(FBT)は、すべての増幅反応のサイクルに適用されることが一般的である。
本発明によると、ベースライン閾値は、各サイクルに対して同一であるか、または相異に適用(割り当て)される。またはベースライン閾値は、すべてのサイクルにおいて相異に適用される。
本発明の一具現例によると、サイクルの少なくとも二つのサイクルは、互いに異なるベースライン閾値を有する。したがって、複数のベースライン閾値は、全体が又は一部が互いに異なっている。複数のベースライン閾値全部または一部が互いに異なる値を有するように、サイクルに適用されるベースライン閾値を調整するこの方法は、‘可変ベースライン閾値’(VBT、variable baseline threshold)方法と命名する。
本発明の一具現例によると、前記ベースライン領域のエンドポイントサイクルは、以下の段階を含む方法によって決定することができる;
(b1)それぞれの増幅サイクルにベースライン閾値を適用して、複数のベースライン閾値をサイクルに適用する段階(前記サイクルの少なくとも二つのサイクルのベースライン閾値は、互いに異なっている)、
(b2)ベースライン閾値のそれぞれによって決定される閾値基準を満足する一つ以上のサイクルを確認する段階と、
(b3)前記段階(b2)で確認されたサイクルを使用してベースライン領域のエンドポイントサイクルを決定する段階。
ベースラインの閾値をサイクルに対してプロットして得られたグラフをBTグラフという。
単一ベースライン閾値がすべてのサイクルに割り当てられる場合、前記BTグラフは、x軸と平行な直線を有する。
VBT方法によると、図9に示されたように、互いに異なるベースライン閾値による多様なBTグラフを収得することができる。
本発明の一具現例によると、増幅反応のサイクルは、少なくとも二つの相異なるグループに分類されて、同一なグループに属するサイクルは、互いに同一なベースライン閾値を有し、互いに異なるグループに属するサイクルは、相異なるベースライン閾値を有する。この場合、BTグラフは、x軸と平行な直線を有する(図9の(a)と(b)を参照)。
本発明の一具現例によると、全体または一部サイクルのベースライン閾値は、サイクル数の増加によって一定比率で増加するか、減少する。この場合、BTグラフは、1次関数の形態で表れる(図9の(c)、(d)、(g)及び(h)参照)。
本発明の一具現例によると、全体または一部サイクルのベースライン閾値は、サイクル数の増加によって変化する比率で増加するか、減少する。この場合、BTグラフは、曲線関数(例えば、2次関数)形態で表れる(図9の(f)参照)。
本発明の一具現例によると、増幅サイクルは、少なくとも二つのグループに分類されて、同一グループに分類されたサイクルは、同一なベースライン閾値を有して、異なるグループに分類されたサイクルは、異なるベースライン閾値を有する。このような場合、BTグラフは、少なくとも二つの関数によって表現することができる。前記BTグラフは、連続的または不連続的にプロットできる。
本発明の一具現例によると、ベースライン閾値変化サイクル(BTCC)が設定され、前記設定されたBTCCの前、後のサイクルには、相異なるベースライン閾値がそれぞれ適用される。
本発明において、用語‘ベースライン閾値変化サイクル’(BTCC; baseline threshold−changed cycle)は、ベースライン閾値の変化パターンがサイクルを超えながら変更される基準サイクルをいう。特に用語‘ベースライン閾値変化サイクル ’(BTCC)は、ベースライン閾値が変化する基準サイクルを意味する。前記BTCCは、一つの増幅反応に対して一つまたはそれ以上設定できる。前記BTCCは、増幅反応前、増幅反応中または増幅反応後に設定できる。
本発明の一具現例によると、増幅サイクルに対するベースライン閾値は、ベースライン閾値変化サイクル(BTCC)を基準に、前記BTCC以前のサイクルとこれらに適用されるベースライン閾値の集合によって形成された第1関数と、前記BTCC以後のサイクルとこれらに適用されるベースライン閾値の集合によって形成された第2関数とが互いに異なっているように設定される。
本発明の一具現例によると、BTCCには、BTCC以前のサイクルに対する関数またはBTCC以後のサイクルに対する関数が適用される。
一つのサイクルに一つのベースライン閾値が適用されるため、サイクルに対するベースライン閾値の関数が形成できる。前記関数は、全体サイクルに対する閾値の集合から形成することができ、一部サイクルに対する閾値の集合から形成することができる。
BTCC以前のサイクルとこれらに適用されるベースライン閾値の集合によって形成された第1関数と、前記BTCC以後のサイクルとこれらに適用されるベースライン閾値の集合によって形成された第2関数とが互いに異なるように設定されるということは、BTCC以前のサイクルの第1関数のBTグラフとBTCC以後のサイクルの第2関数のBTグラフとが互いに異なるパターンを示すことを意味する。
たとえば、BTCC以前のサイクルに適用されたベースライン閾値は、1次関数で表すことができ、BTCC以後のサイクルに適用されたベースライン閾値は、定数関数で表すことができる(図9の(c)参照)。図9の(c)において、BTCC以前のサイクルに適用されたベースライン閾値は一定に減少し、BTCC以後のサイクルに適用されたベースライン閾値は、固定された値を有する。
他の例として、二つのベースライン閾値を有する場合、一番目のBTCC(BTCC1)以前のサイクルに対するベースライン閾値と二番目のBTCC(BTCC2)以後のサイクルに対するベースライン閾値は、定数関数で表し、BTCC1とBTCC1間のベースライン閾値は、前記二つの定数関数を連結する1次関数(図9の(d)参照)、または2次関数または他の関数(図9の(f)参照)で表すことができる。
BTCC前のBTグラフとBTCC後のBTグラフは、BTCCを起点に不連続的な形態でありえる(図9の(e)と(h)参照)。このような場合、BTCCを、BTCC前の関数とBTCC後の関数にそれぞれ適用すると、互いに異なる関数値を算出する。
本発明の一具現例によると、前記増幅サイクルは、ベースライン閾値変化サイクル(BTCC)を境に、2以上の異なるグループに分類される。この場合、同一なグループに属するサイクルは、連続的で、同一なグループに属するサイクルは、互いに同一なベースライン閾値を有し、互いに隣接した別のグループに属するサイクルは、相異なるベースライン閾値を有する。したがって、離隔された別のグループに分類されたサイクルは、互いに異なるか、同一な閾値を有する。BTCCのベースライン閾値は、BTCC以前のサイクルまたは以後のサイクルのベースライン閾値と同一でありえる。BTCCは、一つのデータセットに一つ設定でき、この場合、データセットのサイクルは、二つのグループに分類される。または前記BTCCは、一つのデータセットに二つ以上設定することができて、この場合、データセットのサイクルは、3つ以上のグループに分類される。
本発明の一具現例によると、増幅サイクルは、少なくとも1つ以上のBTCCを基準に、少なくとも二つの相異なるグループに分類され、同一なグループに属するサイクルは、同一なベースライン閾値を有することができる。非特異的信号またはノイズ信号が激しく表れるサイクル区間には、高いまたは低いベースライン閾値が適用されて、非特異的またはノイズ信号が正常信号として検出されないようにすることができる。また、その他のサイクル区間には、一般的なベースライン閾値が適用され、正常信号を検出および分析する。
好ましくは、増幅サイクルは、BTCCを基準に、互いに異なる二つのグループに分類され、同一なグループは、同一なベースライン閾値を有して、異なるグループに分類されたグループは、異なるベースライン閾値を有する。BTCCは、BTCCの前または後のサイクルと同一なベースライン閾値を有することができる。
本発明の一実施例では、前記VBT方法を、初期サイクルが非正常的に高い傾き値を有する増幅結果に使用した。BTCCを設定した後、増幅曲線の補正のために、高いベースライン閾値をBTCC以前のサイクルに適用して、低ベースライン閾値をBTCC以後のサイクルに適用した。VBTを適用した場合が、さらに正確に増幅曲線を補正できることを確認した。
前記BTCCは、測定機器、試料、試薬の特性に影響を受ける背景信号またはノイズ信号のパターンによって変わる。BTCCは、背景信号またはノイズ信号による誤りが防止される限り、その決定範囲に制限がない。たとえば、BTCCは、70,60,50,40,30,29,28,27,26,25,24,23,22,21,10または15サイクル以下でありえる。前記BTCCサイクルは、1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,25,30,35または40サイクル以上でありえる。BTCCは、1乃至70サイクル範囲内で決定することができ、好ましいサイクル範囲は、1〜60、1〜50、1〜40、1〜30、5〜60、5〜50、5〜40、10〜40、10〜35、15〜35、15〜30、15〜25であり、さらに好ましくは、15〜25サイクルである。
図10A及び図10Bに示されたように、BTCCを基準に、第1BTグラフ(1st BT)と第2BTグラフ(2nd BT)で表示される二つの関数が定数関数であると、傾き曲線との一番目のクロスポイントサイクルがベースライン領域のエンドポイントサイクル(E2)に決定される。
MBEC方法とVBT方法は、選択的に適用することができる。また、前記二つの方法を両方とも増幅曲線の補正に適用することができる。
実施例2では、インフルエンザAウイルス(Flu A)の核酸を収得して、増幅し、第1増幅曲線を確保した。ベースライン領域のエンドポイントサイクルをMBEC方法またはVBT方法をそれぞれ適用するか、両方とも適用せずに、それぞれ収得した。以後、得られたエンドポイントサイクルに決定されたベースライン領域を使用し、第1増幅曲線を補正した。その結果、MBEC方法またはVBT方法を適用した補正は、前記方法を適用していない補正に比べて、増幅産物の量をより正確に反映する、補正された増幅曲線を提供することができることを確認した。
(3)多重ベースライン閾値セット方法
本発明の一具現例によると、追加的なベースライン閾値が、サイクルの少なくとも一つのサイクルに適用できる。
本発明の一具現例によると、追加的なベースライン閾値セットが、データセットに適用できる。
一般的な核酸増幅では、標的核酸分子がない場合、増幅曲線は、平らな形状(flat shape)であるか、非特異的結合及び増幅によって少し増加するパターンである。可能性が低いが、非正常的な陰性増幅反応は、サイクルの増加に応じて、信号値の減少を示す(すなわち、陰性傾きパターン)。
このような場合、ベースライン領域の決定および生データセットの補正において、誤りが発生しやすいが、これは、ベースライン閾値セットと増幅曲線(または傾き曲線)との間のクロスポイントが発生せず、したがって、ベースライン領域のエンドポイントサイクルが設定されないからである。本発明の追加ベースライン閾値のセットの適用は、このような誤りのある分析を防止することができる。例えば、追加ベースライン閾値セットのベースライン閾値が負に適用され、したがって、負の傾き値を有するベースライン領域のエンドポイントサイクルを決定することができる。
本発明の一具現例によると、少なくとも二つの閾値セットがデータセットに適用される。増幅曲線の分析において誤り防止のために、少なくとも二つの閾値セットをデータセットに適用する方法は、‘多重ベースライン閾値セット’方法と命名される。ベースライン閾値セットは、ベースライン閾値ポイントの集合である。前記ベースライン閾値ポイントは、サイクルと、前記サイクルにおけるベースライン閾値を含む座標値をいう。
本発明の一具現例によると、ベースライン閾値セットは、サイクルの少なくとも二つのサイクルに適用されたベースライン閾値が互いに異なる、可変ベースライン閾値セットである。
本発明の一具現例によると、ベースライン閾値セットは、すべてのサイクルに適用されるベースライン閾値が同一な、固定ベースライン閾値セットである。
本発明の一具現例によると、多重ベースライン閾値セットは、可変ベースライン閾値セットと固定ベースライン閾値セットとを含むグループから選択された少なくとも二つのベースライン閾値セットを含む。
本発明の一具現例によると、多重ベースライン閾値セットは、少なくとも一つの固定ベースライン閾値セットを含む。
本発明の一具現例によると、多重ベースライン閾値セットは、少なくとも一つの可変ベースライン閾値セットを含む。
本発明の一具現例によると、多重ベースライン閾値セットは、少なくとも一つの固定ベースライン閾値セット及び少なくとも一つの可変ベースライン閾値セットを含む。
本発明の一具現例によると、多重ベースライン閾値セットは、少なくとも二つのベースライン閾値セットを含み、前記少なくとも二つのベースライン閾値セットは、すべて固定ベースライン閾値セットである。
本発明の一具現例によると、多重ベースライン閾値セットの各ベースライン閾値セットは、対応する閾値基準を有する。
本発明の一具現例によると、多重ベースライン閾値セットは、同時に分析対象データセットに適用することができる。閾値基準を満足する複数のサイクルが観察される場合、すべてがまたはその一部が分析に使用できる。
本発明の一具現例によると、多重ベースライン閾値セットは、順次に分析対象データセットに適用することができる。最初に適用されたベースライン閾値セットの閾値基準を満足するサイクルがない場合、また他のベースライン閾値セットが適用される。
本発明の一具現例によると、ベースライン閾値セットがデータセットに適用されて、追加的なベースライン閾値セットがサイクルの少なくとも一つのサイクルに適用できる。
本発明の一具現例によると、ベースライン領域のエンドポイントサイクルは、以下の段階を含む方法によって決定できる;
(a)各サイクルにベースライン閾値を適用する段階と、
(b)少なくとも一つのサイクルに追加ベースライン閾値を適用する段階と、
(c)各ベースライン閾値によって決定されたベースライン閾値基準を満足する1つ以上のサイクルを確認する段階と、
(d)前記段階(c)で確認されたサイクルを使用してベースライン領域のエンドポイントサイクルを決定する段階。
段階(a)のベースラインの閾値は、正常的な増幅結果のエンドポイント決定のために適用される。
段階(b)の追加的なベースライン閾値は、段階(a)のベースライン閾値によって決定された主閾値基準(main threshold criterion)を満足するサイクルがない、非正常的な増幅結果に適用される。
さらに具体的に、閾値基準を満足するサイクルは、以下の方法で確認される;n番目サイクルの差引結果の符号が(n−1)番目サイクルの差引結果の符号と相異なる場合、n番目サイクルは、閾値基準を満足するものである。前記差引結果は、信号値から閾値を差し引いた結果である。
エンドポイントサイクルを決定するために、さまざまなデータセットが使用できる。適切な閾値及び決定方法は、データセットの種類に応じて選択できる。当業者は、本発明の方法を、上述の傾きデータセットの具現例及び実施例の指示および方向に基づいて、多様なデータセットの分析に利用することができる。
段階(c):最適ライン関数の算出(S130)
ベースライン領域の決定に続いて、ベースライン領域の最適ライン関数が、ベースライン領域内生データセットの二つ以上のデータポイントから計算される。
最適ライン関数は、データにあらわれた傾向を最もよく表す関数を意味する。最適ラインは、前記最適ライン関数をプロットして得られたグラフを意味する。
最適ライン関数は、ベースライン領域内の少なくとも二つのデータポイントを使用して求めることができる。例えば、ベースライン領域の一部またはすべての地域のデータポイントを使用して求めることができる。
最適ライン関数は、様々な方法で求めることができる。例えば、ベースライン領域内のデータポイントを使用した線形回帰分析(linear regression analysis)またはLMS(least mean square)アルゴリズムによって求めることができる。
特に、線形回帰の1次方程式、‘y=mx+b’で表示される最適ライン関数は、スタートポイントサイクルからエンドポイントサイクルまでのデータポイントを使用して求めることができる。
図4Aに示されたように、最適ライン傾きの‘m’及び最適ラインのy軸切片の‘b’は、下記数学式2及び3によって計算することができる。
(前記数学式2及び3において、
であり、mは、最適ラインの傾きであって、bは、y軸切片であり、Xは、i番目サイクルのサイクル番号であって、yは、i番目サイクルで測定された信号値であり、Sは、ベースライン領域のスタートポイントサイクルであって、Eは、ベースライン領域のエンドポイントサイクルであり、nは、E−S+1である。)
段階(d):補正されたデータセットの収得(S140)
補正されたデータセットは、生データセットの信号値から最適信号関数の値を差し引いて得られる。前記得られたデータセットは、(i)増幅反応の増幅信号値及び(ii)減算された結果を含む。
本発明の一具現例によると、前記段階(a)は、第1増幅曲線を得るために、生データセットをプロットする段階をさらに含み、前記段階(d)は、補正された増幅曲線を得るために、補正されたデータセットをプロットする段階をさらに含む。
補正された増幅曲線(第2増幅曲線)は、補正前の増幅曲線(第1増幅曲線)から最適ラインを差し引いて得ることができる。図5に示したように、生データセットの信号値から、最適ライン関数から計算された値を差し引いて、補正されたデータセット及び補正された増幅曲線を収得する。
本発明の一具現例によると、生データセットの補正は、最適ライン関数で算出された値を、ベースライン領域だけでなく、ベースライン領域ではないその他の領域でも差し引くことを含む。例えば、すべてのサイクルに対して生データセットを得て、生データセットの信号値から最適ライン関数の値を差し引いて補正されたデータセットを求め、続いて前記補正されたデータセットをプロットし、補正された増幅曲線を得ることができる。または、ベースライン領域に属するサイクルの生データセットを得て、最適ライン関数及び補正されたデータセットを得て、他のサイクルに対しては、サイクルから信号発生の度に生データセットを得て、前記生データセットから最適ライン関数による該サイクルに対応する値を差し引いて、補正されたデータセットを得て、これをプロットし、補正された増幅曲線を得ることができる。
III.記録媒体、装置及びコンピュータプログラム
以下、記載された本発明の記録媒体、装置およびコンピュータプログラムは、本発明の方法をコンピュータで実行するためのものであるため、その共通する内容は、明細書の複雑性を避けるために省く。
本発明の一様態によると、本発明は、以下の段階を含む、試料を分析する方法を実行させるためのプロセッサを具現する指示を含む、コンピュータ読み取り可能な記録媒体を提供する;
(a)サイクルにおける信号値を提供するために、前記試料を使用する信号発生反応のそれぞれのサイクルにおける信号値を収得する段階と、
(b)それぞれのサイクルに閾値を適用し、複数の閾値をサイクルに適用する段階(前記サイクルの少なくとも二つのサイクルの閾値は、互いに異なっている)と、
(c)閾値のそれぞれにより定められる閾値基準を満足する一つ以上のサイクルを確認する段階と、
(d)前記段階(c)で確認されたサイクルを使用して試料を分析する段階。
本発明の一具現例によると、前記信号発生反応は、試料内の分析物質の存在に依存的に信号を発生する。この際、前記分析物質は、標的核酸分子であり、前記信号発生反応は、標的核酸分子の増幅を伴うか、あるいは伴わない反応であって、前記標的核酸分子の増幅を伴う反応は、リアルタイムPCR(real−time polymerase chain reaction)である。
本発明の一具現例によると、前記信号発生反応は、リアルタイムPCRであって、信号値は、前記リアルタイムPCRにより発生された信号値が数学的に変形された値であり、前記数学的に変形された値は、サイクルに対して信号値を微分した値である。
本発明の一具現例によると、前記閾値は、閾値変化サイクル(threshold−changed cycle、TCC)を基準に、前記TCC以前のサイクル及びこれらに適用される閾値の集合によって形成された関数と、前記TCC以後のサイクル及びこれらに適用される閾値の集合によって形成された関数とが互いに異なるように決定される。
本発明の他の一様態によると、本発明は、以下の段階を含む、信号発生手段の存在下で行われた増幅反応から得た生データセット(raw data set)を補正する方法を実行させるためのプロセッサを具現する指示を含む、コンピュータ読み取り可能な記録媒体を提供する;
(a)(i)前記増幅反応における増幅サイクル及び(ii)前記信号発生手段から得た前記増幅サイクルにおける信号値を含む生データセットを収得する段階と、
(b)前記生データセットからベースラインスタートポイントサイクルとエンドポイントサイクルを決定して、ベースライン領域を決定する段階と、
(c)前記決定されたベースライン領域内生データセットの二つ以上のデータポイントからベースライン領域の最適ライン関数を算出する段階と、
(d)前記生データセットの信号値から前記ベースライン領域の最適ライン関数により計算された関数値を差し引き、補正されたデータセットを収得する段階(前記補正されたデータセットは、(i)増幅反応の増幅サイクル及び(ii)前記差し引いた結果を含む)。
本発明の一具現例によると、前記段階(b)においてベースライン領域のエンドポイントサイクルは、以下の段階を含む方法により決定される;
(b1)それぞれの増幅サイクルにベースライン閾値を適用して、複数の閾値をサイクルに適用する段階(前記サイクルの少なくとも二つのサイクルのベースライン閾値は、互いに異なっている)、(b2)ベースライン閾値のそれぞれによって決定される閾値基準を満足する一つ以上のサイクルを確認する段階と、(b3)前記段階(b2)で確認されたサイクルを使用してベースライン領域のエンドポイントサイクルを決定する段階。
本発明の一具現例によると、前記段階(b1)の傾きは、特定サイクルにおけるデータポイント及び前記特定サイクルの以前および/または以後のサイクルの少なくとも一つのデータポイントを利用した最小二乗法によって計算することができる。
本発明の一具現例によると、本発明は、以下の段階を含む、試料を分析する方法を実行するためのプロセッサを具現する、コンピュータ読み取り可能な記録媒体に格納されたコンピュータプログラムを提供する;
(a)サイクルにおける信号値を提供するために、前記試料を使用する信号発生反応のそれぞれのサイクルにおける信号値を収得する段階と、
(b)それぞれのサイクルに閾値を適用し、複数の閾値をサイクルに適用する段階(前記サイクルの少なくとも二つのサイクルの閾値は、互いに異なっている)と、
(c)閾値のそれぞれによって定められる閾値基準を満足する一つ以上のサイクルを確認する段階と、
(d)前記段階(c)で確認されたサイクルを使用して試料を分析する段階。
本発明の一具現例によると、前記信号発生反応は、試料内の分析物質の存在に依存的に信号を発生させ、前記の分析物質は、標的核酸分子であり、前記信号発生反応は、標的核酸分子の増幅を伴うかあるいは伴わない反応であって、前記標的核酸分子の増幅を伴う反応は、リアルタイムPCR(real−time polymerase chain reaction)である。
本発明の一具現例によると、前記信号発生反応は、リアルタイムPCRであり、信号値は、前記リアルタイムPCRにより発生された信号値が数学的に変形された値であって、前記数学的に変形された値は、サイクルに対して信号値を微分した値である。
本発明の一具現例によると、前記閾値は、閾値変化サイクル(threshold−changed cycle、TCC)を基準に、前記TCC以前のサイクル及びこれらに適用される閾値の集合によって形成された関数と、前記TCC以後のサイクル及びこれらに適用される閾値の集合によって形成された関数とが互いに異なるように決定される。
本発明の他の一様態によると、本発明は、以下の段階を含む、信号発生手段の存在下で行われた増幅反応から得た生データセットを補正する方法を実行するためのプロセッサを具現する、コンピュータ読み取り可能な記録媒体に格納されたコンピュータプログラムを提供する;
(a)(i)前記増幅反応における増幅サイクル及び(ii)前記信号発生手段から得た前記増幅サイクルにおける信号値を含む生データセットを収得する段階と、
(b)前記生データセットからベースラインスタートポイントサイクルとエンドポイントサイクルを決定して、ベースライン領域を決定する段階と、
(c)前記決定されたベースライン領域内生データセットの二つ以上のデータポイントからベースライン領域の最適ライン関数を算出する段階と、
(d)前記生データセットの信号値から前記ベースライン領域の最適ライン関数により計算された関数値を差し引き、補正されたデータセットを収得する段階(前記補正されたデータセットは、(i)増幅反応の増幅サイクル及び(ii)前記差し引いた結果を含む)。
本発明の一具現例によると、前記信号発生手段の存在下で行われた増幅反応から得た生データセットを補正する方法を実行するためのプロセッサを具現する、コンピュータ読み取り可能な記録媒体に格納されたコンピュータプログラムを提供する。
本発明の一具現例によると、前記段階(b)においてベースライン領域のエンドポイントサイクルは、以下の段階を含む方法により決定される;
(b1)それぞれの増幅サイクルにベースライン閾値を適用して、複数の閾値をサイクルに適用する段階(前記サイクルの少なくとも二つのサイクルのベースライン閾値は、互いに異なっている)、(b2)ベースライン閾値のそれぞれによって決定される閾値基準を満足する一つ以上のサイクルを確認する段階と、(b3)前記段階(b2)で確認されたサイクルを使用してベースライン領域のエンドポイントサイクルを決定する段階。

本発明の一具現例によると、前記段階(b1)の傾きは、特定サイクルにおけるデータポイント及び前記特定サイクルの以前および/または以後のサイクルの少なくとも一つのデータポイントを利用した最小二乗法によって計算することができる。
前記プログラム指示は、プロセッサによって実行される場合、プロセッサが上述の本発明の方法を実行するように運営できる。試料分析方法を実行するためのプログラム指示は、サイクルにおける信号値を提供するために、前記試料を使用する信号発生反応のそれぞれのサイクルにおける信号値を収得する指示と、それぞれのサイクルに閾値を適用して、閾値のそれぞれによって決定される閾値基準を満足する一つ以上のサイクルを確認する指示と、前記確認されたサイクルを使用して試料を分析する指示とを含む。増幅反応から得られた生データセットを補正する方法を実行するためのプログラム指示は、生データセットを収得する指示と、前記生データセットから、ベースラインスタートポイントサイクルとエンドポイントサイクルを決定してベースライン領域を決定し、前記ベースライン領域の最適ライン関数を算出する指示と、前記生データセットの信号値から、前記ベースライン領域の最適ライン関数によって計算された関数値を差し引いて、補正されたデータセットを収得する指示とを含む。
本発明の方法は、プロセッサで実行され、前記プロセッサは、スタンドアロンコンピュータ、ネットワーク付きコンピュータまたはリアルタイムPCR装置のようなデータ収得装置のプロセッサである。
コンピュータ読み取り可能な記録媒体は、当業界に公知された様々な記憶媒体、例えば、CD−R、CD−ROM、DVD、フラッシュメモリ、フロッピー(登録商標)ディスク、ハードディスクドライブ、ポータブルHDD、USB、磁気テープ、MINIDISC、不揮発性のメモリカード、EEPROM、光学ディスク、光学記憶媒体、RAM、ROM、システムメモリ、およびウェブサーバーを含むが、これらに限定されるものではない。
データポイント(例えば、信号強度及び増幅サイクル)は、多様な方式でで収集することができる。たとえば、データポイントは、PCRデータ収集装置にあるプロセッサによって収集できる。データポイントは、収得され次第実時間でプロセッサに提供できて、または、メモリユニットまたはバッファに格納され、実験終了後、プロセッサに提供することができる。同様に、データセットは、前記収集装置とのネットワーク接続(例えば、LAN、VPN、インターネットおよびイントラネット)または直接接続(例えば、USBまたは他の直接有線接続または無線接続)によってデスクトップコンピュータシステムなどの別のシステムに提供でき、またはCD、DVD、フロッピー(登録商標)ディスクおよびポータブルHDDのようなポータブルメディア上に提供することができる。同様に、データセットは、ノートブックまたはデスクトップコンピュータシステムなどのクライアントにネットワーク接続(例えば、LAN、VPN、インターネット、イントラネットおよび無線通信ネットワーク)を介してサーバシステムに提供することができる。
データポイントが受信されるか得られた後、試料を分析するか、増幅反応の補正されたデータセットを得るためのデータ分析プロセスが進行される。例えば、試料分析のためのプロセッサは、各閾値によって決定された閾値基準を満足する一つまたはそれ以上のサイクルを確認するために、受信されたデータポイントを分析する。増幅反応の補正されたデータセットを得るためのプロセッサは、ベースライン領域を決定し、ベースライン領域の最適ライン関数を算出して、生データセットの信号値から、最適ライン関数によって計算された関数値を差し引き、補正されたデータセットを収得するために、受信されたデータポイントを分析する。
本発明を実行するためのプロセッサを具現する指示は、ロジックシステムに含まれる。前記指示は、たとえポータブルHDD、USB、フロッピー(登録商標)ディスク、CD及びDVDのようなソフトウェア記録媒体に提供できるが、メモリモジュール(例えば、ハードドライブ、またはローカルまたは付着RAMまたはROMのような他のメモリ)にダウンロード及び格納できる。本発明を実行するコンピュータコードは、C、C++、Java(登録商標)、Visual Basic、VBScript、JavaScript(登録商標)、Perl及びXMLのような多様なコーディング言語で実行できる。また、多様な言語及びプロトコルは、本発明によるデータと命令の外部及び内部記憶と伝達に利用できる。
本発明のまた他の様態によると、本発明は、(a)コンピュータプロセッサ、および(b)前記コンピュータプロセッサにカップリングされた前記本発明のコンピュータ読み取り可能な記録媒体を含む、試料を分析する装置を提供する。
本発明のまた他の様態によると、本発明は、(a)コンピュータプロセッサ、および(b)前記コンピュータプロセッサにカップリングされた前記本発明のコンピュータ読み取り可能な記録媒体を含む、信号発生手段の存在下で行われた増幅反応の生データセットを補正する装置を提供する。
本発明の一具現例によると、前記装置は、試料及び信号発生手段を収容できる反応容器、前記反応容器の温度を調節する温度調節手段及び/または増幅サイクルにおける信号を検出する検出器をさらに含むことができる。
本発明の一具現例によると、前記コンピュータのプロセッサは、サイクルにおける信号値を受信するだけでなく、試料を分析するか、増幅反応の補正されたデータセットを得る。
前記プロセッサは、一つのプロセッサが下記二つのパフォーマンスを実行するように構築することができる:データポイントの受信及び試料の分析や補正されたデータセットを収得。
または、プロセッサユニットは、二つのプロセッサが前記二つのパフォーマンスをそれぞれ実行するように構築することができる。
本発明の一具現例によると、プロセッサは、標的分析物質(例えば、標的核酸分子)の検出に用いられる従来の装置(例えば、リアルタイムPCR装置)にソフトウェアをインストールして具現することができる。
図6は、増幅反応の生データセットを補正するための本発明の一実施例に最適化されたリアルタイムPCRシステムを示している。前記システムは、リアルタイムPCR増幅を実行するためのリアルタイムPCR装置110、生データセットを補正し、補正された結果をディスプレイするために、ケーブル130を介して前記リアルタイムPCR装置110に連結されたロジックシステムのコンピュータシステム120を含む。前記コンピュータシステム120は、補正結果をユーザーのニーズに応じてグラフ、テーブル(表)及びワードなど、様々な形でディスプレイすることができる。前記コンピュータシステムは、増幅反応の増幅曲線を補正するための本発明の方法を実行するための、コンピュータ読み取り可能な記録媒体を含む。前記リアルタイムPCR装置110及び前記コンピュータシステム120は、システムに統合することができる。
増幅曲線に関するデータポイント(例えば、信号強度及び増幅サイクル)は、多様な方法で収集できる。たとえば、データポイントは、リアルタイムPCR装置110のデータ収集装置内のプロセッサによって受信し収集することができる。データセットは、実時間でプロセッサに提供されるか、またはメモリユニットまたはバッファに格納されて、実験終了後、プロセッサに提供される。
同様に、データセットは、ネットワーク接続(例えば、LAN、VPN、インターネット及びイントラネット)または直接接続(例えば、USBまたは他の直接有線接続または無線接続)によって、またはCD、DVD、フロッピー(登録商標)ディスクおよびポータブルHDDのようなポータブルメディアを介して、リアルタイムPCR装置110からデスクトップコンピュータシステムのようなコンピュータシステム120に提供できる。同様に、データセットは、ノートブックまたはデスクトップコンピュータシステムのようなクライアントにネットワーク接続(例えば、LAN、VPN、インターネット、イントラネット、および無線通信ネットワーク)を介してサーバシステムに提供できる。
データセットが受信されて得られた後、データ分析プロセッサは、補正された増幅曲線を反映するデータセットを準備する。
増幅曲線の補正は、コンピュータシステム120にインストールされているアプリケーション(つまり、プログラム)によって実行できる。択一的に、増幅曲線の補正は、コンピュータシステム120の運営システムにおいて、前記アプリケーションが運営できるコンピュータシステム120にアプリケーションストアサーバーまたはアプリケーションプロバイダサーバを介して直接インストールされたアプリケーションによって作られる。前記運営システムは、Window、Macintoshおよびスマートフォンや、タブレットPCのようなモバイル端末にインストールされるiOS及びアンドロイドのようなモバイル運営システムを含む。
上述のように、本発明の増幅曲線補正方法は、供給者がインストールした、またはユーザーが直接インストールしたアプリケーション(即ち、プログラム)によってコンピュータシステム120に具現でき、コンピュータ読み取り可能な記録媒体122に記録されえる。
本発明の増幅曲線を補正するための方法が具現されたコンピュータプログラム124は、前記補正のためのすべての機能を具現することができる。前記コンピュータプログラム124は、本発明の方法を実行するように具現されたコンピュータ読み取り可能な記録媒体に格納されたプログラム指示を含めたプログラムである。
コンピュータプログラム124は、C、C++、JAVA(登録商標)、Visual basic、VBScript、JavaScript(登録商標)、Perl、XML及び機械語などの適切なコンピュータ言語を使用してコードすることができる。前記プログラムコードは、前記言及された数学的関数に対する関数コード及びコンピュータシステム120のプロセッサによって順番にプロセスを移行するための制御コードを含む。
前記コードは、プロセッサによって、上述した機能の移行に必要な追加情報またはメディアが、コンピュータシステム120の内蔵または外付けメモリのどの位置(アドレス)で参照されるかに対するメモリ参照コードをさらに含む。
コンピュータシステム120がプロセッサの機能の実行のために、他のコンピュータまたはサーバーとリモートで通信が必要な場合、前記コードは、プロセッサがどのように通信モジュール(例えば、有線および/または無線通信モジュール)を使用して他のコンピュータまたはサーバーとリモートで通信するか、または通信時、どのような情報またはメディアを伝達するかに対する通信関連コードをさらに含むことができる。
本発明の具現のための機能的プログラム及びコード(コードセグメント)は、コンピュータが読み取るストレージメディア及び実行プログラムなど、システム環境を考慮して、当業界のプログラマによって容易に推論または修正できる。
記録媒体122が配布されて、コンピュータで読み取り可能なコードが記憶されて、配布形式で実行される。このような場合、複数の配布されたコンピュータのうち、少なくとも一つのコンピュータに、前記記録媒体122は、ネットワークで接続されたコンピュータシステムに分散でき、噴射方式で、コンピュータ読み取り可能なコードが保存されて、実行できる。このような場合、複数の分散されたコンピュータのうち、少なくとも一つのコンピュータは、前記機能の一部を実行して、その実行結果を少なくとも一つのコンピュータに伝達し、結果が伝達されたコンピュータも前記機能の一部を実行して、その結果を少なくとも一つの他のコンピュータに提供する。
本発明の実行のためのアプリケーション(つまり、プログラム)の記録された記録媒体122は、アプリケーションストアサーバーまたはアプリケーションプロバイダサーバーに含まれた記録媒体(例えば、ハードディスク)、またはアプリケーションプロバイダサーバそのものであるか、プログラムを記録した別のコンピュータ及びその記録媒体でありえる。
記録媒体122を読み取り可能なコンピュータシステム120は、一般的なデスクトップやノート型パソコンなどの一般PCだけではなく、スマートフォン、タブレットPC、PCA(Personal Digital Assistants)及び移動通信端末などのモバイル端末を含み、コンピューティング可能なあらゆるデバイスを含む。
以下、実施例を通じて本発明をさらに詳細に説明するが、これら実施例は、本発明をより具体的に説明するためのものであって、本発明の範囲がこれら実施例に限定されないことは、本発明の属する技術分野で通常の知識を有する者にとっては自明なことであろう。
実施例1:増幅曲線の補正
図6に示されたリアルタイムPCRシステムを使用して、増幅曲線の傾き曲線及びベースライン閾値から得たベースライン領域の最適ラインによって増幅曲線を補正することができるかを確認した。
生データセット(補正前の増幅曲線)の収得(S110)
アップストリームプライマーとダウンストリームプライマーの伸長およびTaqManプローブの切断のために、5’ヌクレアーゼ活性を有するTaq DNA重合酵素を使用した。Neisseria gonorrhoeae(NG)のゲノムDNAを標的核酸配列として使用した。
NGの検出のために、TaqManリアルタイムPCRを使用した。標的核酸が存在すると、TaqManプローブが切断されて、標識された断片(fragment)が放出される。増幅曲線は、標識された断片からの信号を測定して得られる。
NGに対するTaqManプローブは、5’末端に蛍光レポーター分子(Cal Fluor Red 610)を標識して、3’末端にクエンチャー分子(BHQ−2)を標識した(配列番号:3)。
本実施例で使用されてアップストリームプライマー、ダウンストリームプライマー、およびプローブの配列は、以下の通りである:
NG−F 5’−TACGCCTGCTACTTTCACGCTIIIIIGTAATCAGATG−3’(配列番号:1)
NG−R 5’−CAATGGATCGGTATCACTCGCIIIIICGAGCAAGAAC−3’(配列番号:2)
NG−P 5’−[Cal Fluor Red 610]TGCCCCTCATTGGCGTGTTTCG[BHQ−2]−3’(配列番号:3)
(I:Deoxyinosine、BHQ−2:Black hole quencher−2)
リアルタイムPCRは、標的核酸(10pg、1pg、100fg、10fg、または1fgのNGゲノムDNA)、NG標的増幅のための5pmoleのアップストリームプライマー(配列番号:1)と5pmoleのダウンストリームプライマー(配列番号:2)、3pmoleのTaqManプローブ(配列番号:3)、および5μlの4×マスターミックス(Master Mix)[最終、200μM dNTPs、2mM MgCl、2UのTaq DNA重合酵素]を含む最終的な容量20μlの反応混合液で行った。前記反応混合液を含むチューブをリアルタイムPCR装置(CFX96TM、Bio−Rad)で、50℃で5分間放置して、次いで95℃で15分間変性させた後、95℃で30秒間、60℃で60秒間および72℃で30秒間行われるサイクルを50回繰り返した。信号検出は、各サイクルごとに60℃で行った。
生データセットは、リアルタイムPCR増幅によって得て、補正前の増幅曲線は、生データセットを使用してプロットした(図2参照)。
補正前の増幅曲線は、以下のように補正した。
ベースライン領域の決定(S120)
増幅反応の3番目のサイクルをベースライン領域のスタートポイントサイクル(5)に設定した。
ベースライン領域のエンドポイントサイクル(E)を決定するために、数学式1で表される最小二乗法(Least Square Method)を使用した線形回帰分析(Linear Regression Analysis、LRA)を生データに対して行い、傾き曲線を得た。この際、i番目のサイクルにおける傾き値を算出するために、i−1、i、i+1の3つのサイクルに対応するデータを使用した(つまり、a=1及びb=1)。
傾き曲線でピークが生成される前の初期サイクルにおける背景信号(Background signal)と交差しないように、ベースライン閾値(Baseline threshold)を‘20’に設定した。ベースライン閾値と傾き曲線との最初のクロスポイント(cross−point、CP)のサイクルをエンドポイントサイクル(E)に決定して、最終的にベースライン領域を決定した(図3参照)。
ベースライン領域の最適ライン関数の算出(S130)
前記決定されたベースライン領域のスタートポイントサイクル(S)からエンドポイントサイクル(E)まで各サイクルのデータを利用し、最小二乗法により回帰直線の1次線形方程式形態の最適ライン関数を算出した(図4A参照)。
算出される線形回帰線1次線形方程式の一般式は、‘y=mx+b’であり、ここで傾きを示す‘m’は、上述の数学式2によって算出し、y軸切片を示す‘b’は、上述の数学式3によって算出した。その結果、算出された線形回帰線の1次線形方程式形態の最適ライン関数は、‘y=2.512x+2396.4’である。
補正されたデータセット(補正された増幅曲線)の収得(S140)
補正前の増幅曲線全体1〜50サイクルから前記算出した最適ライン関数値を差し引いて、補正されたデータセットを収得し、補正された増幅曲線を得するために、前記補正されたデータセットをプロット(plotting)した。
図4Bは、図2のベースライン領域の生データセットの各サイクルのデータ(測定された信号値)からベースライン領域の最適ライン関数によるデータ(関数値)を差し引いて得られた、補正されたデータセットを示す。
図5は、前記補正されたデータセットをプロットして得られた、補正された増幅曲線を示す。
本発明の方法は、単純なアルゴリズムを使用して増幅曲線の補正を行うことができるため、本発明は、測定された試料及び測定機器に対する特定条件をさらに簡単な方法で最適化することが可能である。
実施例:増幅曲線の補正(II)
リアルタイム重合酵素連鎖反応を行って得られた増幅曲線を補正するためのベースライン領域決定における誤りを除去することができるかどうか確認した。
アップストリームプライマーとダウンストリームプライマーの伸長およびTaqManプローブの切断のために、5’ヌクレアーゼ活性を有するTaq DNA重合酵素を使用した。インフルエンザAウイルス(Flu A)のゲノムDNAを標的核酸配列として使用した。
TaqManリアルタイムPCRを使用してFlu Aを検出した。標的核酸が存在すると、TaqManプローブが切断されて、標識された断片が放出される。増幅曲線は、標識された断片からの信号を測定して得ることができる。
Flu Aに対するTaqManプローブは、5’末端に蛍光レポーター分子(FAM)を標識して、3’末端にクエンチャー分子(BHQ−1)を標識した(配列番号:6)。
本実施例で使用されたアップストリームプライマー、ダウンストリームプライマー、およびプローブの配列は、以下の通りである:
Flu A−F5’−TGGAATGGCTAAAGACAAGACCIIIIITGTCACCTCT−3’(配列番号:4)
Flu A−R5’−CATCCTGTTGTATATGAGGCCCATIIIICTGGCAAG−3’(配列番号:5)
Flu A−P5’−[FAM]CTCACTGGGCACGGTGAGCGTGA[BHQ−1]−3’(配列番号:6)
(I:Deoxyinosine、BHQ−1:Black hole quencher−1)
リアルタイムPCRは、標的核酸(抽出したFlu AゲノムRNAの10−3、10−4、10−5、または10−6希釈)、Flu A標的増幅のための5pmoleのアップストリームプライマー(配列番号:4)及び5pmoleのダウンストリームプライマー(配列番号:5)、3pmoleのTaqManプローブ(配列番号:6)、及び5μlの5×RT−PCR緩衝液[75mM Tris−HCl(pH8.3)、50mM KCl、2.5mM MgCl、0.2mM dNTP]及び2μlの酵素混合物[最終、3.5UのTaq DNA重合酵素、25UのMMLV逆転写酵素、5UのRNase抑制剤]を含む最終容量25μlの反応混合液で行った。前記反応混合液を含むチューブをリアルタイムPCR装置(CFX96TM、Bio−Rad)で、50℃で20分間放置し、次いで95℃で15分間変性させた後、95℃で10秒間、60℃で60秒間および72℃で10秒間行われるサイクルを45回繰り返した。信号検出は、各サイクルごとに60℃で行った。
ベースライン領域の決定段階(S120)での誤りの確認
ベースライン領域のエンドポイントサイクル(E)は、増幅反応で、実際標的信号強度の増加が始まるサイクルまたはこのサイクルに近接したサイクルに決定することができる。
ベースライン領域のエンドポイントサイクルは、各サイクルで計算される傾き及び各サイクルにおける閾値を全部考慮して決定することができる。この場合、ベースライン閾値は、実施例1のように、すべての増幅サイクルに対して適用することができる。しかし、このような接近法は、ベースライン領域を決定する際、誤りを発生させるおそれがある(図7参照)。
例えば、ベースライン閾値を‘30’のように低く設定する場合、実際標的信号強度が増加し始めるポイントではなく、初期ノイズ信号が発生したポイントがベースライン領域のエンドポイントサイクルに決定されることがあって、ベースライン領域決定誤りが発生する可能性がある。一方、このような誤りを回避するために、ベースライン閾値を‘300’のように高く設定すると、低濃度の標的配列を含む試料(すなわち、傾き値が低い試料)からは、実際標的信号強度が増加し始めるポイントが検出されず、これも同様にベースライン領域決定に誤りが発生するおそれがある。
このように、ベースライン領域決定誤りにより、正確な実際増幅産物の量が反映されていない、補正された増幅曲線が得られる可能性がある。
MBEC方法によるエンドポイントのサイクルの決定
増幅反応の3番目のサイクルをベースライン領域のスタートポイントサイクル(S)に設定した。
図8A及び図8Bに示すように、傾き曲線でエンドポイントサイクルは、最小ベースラインエンドポイントサイクル(minimum baseline end−point cycle、MBEC)以後のサイクルに決定できる。実施例2において、10番目のサイクルがMBECに決定された。
図8A(高濃度試料)及び図8B(低濃度試料)で示されたように、MBECを適用しないと、エンドポイントサイクルがベースライン閾値と傾き曲線間の最初のクロスポイント(CP)のサイクルに指定されて、ベースライン領域(B)は、3〜7サイクル領域(高濃度試料)または3〜4サイクル領域(低濃度試料)に決定される。一方、MBECを適用すると、MBECを超える領域における最初のクロスポイント(CP)のサイクルがエンドポイントサイクルに指定されて、ベースライン領域(B)は、3〜29サイクル領域(高濃度試料)または3〜38サイクル領域(低濃度試料)に決定される。
図8Cに示されたように、ベースライン閾値を‘30’に決定して、MBECを適用しないと、補正された増幅曲線から増幅産物の量を不正確に反映するか、偽陰性結果と解釈される。しかし、ベースライン閾値を‘30’に決定して、MBECを適用すると、補正された増幅曲線は、増幅産物の量を正確に反映するようになる。
このように、増幅反応の初期サイクルで頻繁に発見されるノイズ信号により、このような初期ノイズ信号が発生するポイントがベースライン領域のエンドポイントサイクルに誤って決定されることがあり、ベースライン決定誤りが発生する。
前記結果によって、ベースライン領域の決定(S120)誤りは、本発明の方法により成功的に除去できることを確認した。
VBT(Variable Baseline Threshold)方法によるエンドポイントサイクルの決定
実施例2で使用したVBT(可変ベースライン閾値)方法は、ベースライン閾値−変化サイクル(baseline threshold−changed cycle、BTCC)を決定し、互いに異なるベースライン閾値を前記BTCC前と後のサイクルにそれぞれ適用する方式で行う。
BTCCを20サイクルに決定し、第1のBT(ベースライン閾値)を1〜20サイクルに適用し、第2のBTを21〜45サイクルにそれぞれ相異に適用した。第1のBTは‘300’に決定し、第2のBTは‘30’に決定した。
図10A(高濃度試料)及び図10B(低濃度試料)に示されたように、VBTを適用せずに、固定ベースライン閾値を‘30’にして、全てのサイクルに適用したとき、高濃度試料の場合、エンドポイントサイクル(E)は、7サイクルに決定されるか、または低濃度試料の場合、エンドポイントサイクル(E)は、4サイクルに決定された。
増幅反応の3番目のサイクルをベースライン領域のスタートポイントサイクル(S)に決定すると、ベースライン領域(B)は、3〜7サイクル領域(高濃度試料)または3〜4サイクル領域(低濃度試料)に決定された。一方、VBTを適用すると、エンドポイントサイクル(E)は、高濃度試料の場合、29サイクルに決定されるか、または低濃度試料の場合、38サイクルに決定された。したがって、ベースライン領域(B)は、3〜29サイクル(高濃度試料)または3〜38サイクル(低濃度試料)に決定された。
図10Cに示されたように、VBTを適用しなかった場合、補正された増幅曲線は、増幅産物の量を不正確に反映するか、偽陰性結果を示した。しかし、VBTを適用したとき、補正された増幅曲線は、増幅産物の量を正確に反映し、偽陰性結果は表れなかった。
したがって、前記結果によると、可変的なベースライン閾値を異なって設定するVBT方法を通じて、初期増幅サイクルで初期ノイズ信号を発生させるポイントをベースライン領域のエンドポイントサイクル(E)に誤って設定する誤りを除去することができることを確認した。また、VBT方法は、増幅サイクルの後半部において、信号増加が始まるポイントをより正確に決定することができ、ベースライン領域の決定(S120)誤りを除去することができた。
実施例2の結果によると、ベースライン領域は、初期サイクルにおける背景信号の干渉無しに決定できることを確認した。
上述のように、本発明は、ベースライン領域決定誤りを除去して増幅曲線を補正することにより、より信頼性の高い正確な方式で増幅曲線を分析することができる。
実施例3:正確なCt値決定による標的核酸の検出及び定量
増幅曲線からCt値を決定する段階における誤りを除去することができるかどうかを確認した。
アップストリームプライマーとダウンストリームプライマーの伸長およびTaqManプローブの切断のために、5’ヌクレアーゼ活性を有するTaq DNA重合酵素を使用した。インフルエンザAウイルス(Flu A)のゲノムDNAを標的核酸配列として使用した。
Flu A検出のために、TaqManリアルタイムPCRを使用した。標的核酸が存在すると、TaqManプローブが切断されて、標識された断片が放出される。増幅曲線は、標識された断片からの信号を測定して得ることができる。
Flu Aに対するTaqManプローブは、5’末端に蛍光レポーター分子(FAM)を標識して、3’末端にクエンチャー分子(BHQ−1)を標識した(配列番号:6)。
本実施例で使用されてアップストリームプライマー、ダウンストリームプライマー、およびプローブの配列は、以下の通りである:
Flu A−F5’−TGGAATGGCTAAAGACAAGACCIIIIITGTCACCTCT−3’(配列番号:4)
Flu A−R5’−CATCCTGTTGTATATGAGGCCCATIIIICTGGCAAG−3’(配列番号:5)
Flu A−P5’−[FAM]CTCACTGGGCACGGTGAGCGTGA[BHQ−1]−3’(配列番号:6)
(I:Deoxyinosine、BHQ−1:Black hole quencher−1)
リアルタイムPCRは、標的核酸(抽出したFlu AゲノムRNAの10−3、10−4、10−5、10−6、または10−7希釈)、Flu A標的増幅のための5pmoleのアップストリームプライマー(配列番号:4)及び5pmoleのダウンストリームプライマー(配列番号:5)、3pmoleのTaqManプローブ(配列番号:6)、及び5μlの5×RT−PCR緩衝液[75mM Tris−HCl(pH8.3)、50mM KCl、2.5mM MgCl、0.2mM dNTP]及び2μlの酵素混合物[最終、3.5UのTaq DNA重合酵素、25UのMMLV逆転写酵素、5UのRNase抑制剤]を含む最終容量25μlの反応混合液で行った。前記反応混合液を含むチューブをリアルタイムPCR装置(CFX96TM、Bio−Rad)で、50℃で20分間放置し、次いで95℃で15分間変性させた後、95℃で10秒間、60℃で60秒間および72℃で10秒間行われるサイクルを45回繰り返した。信号検出は、各サイクルごとに60℃で行った。
Ct値決定誤りの確認
増幅曲線から標的核酸の正確な量を得るための伝統的なサイクル閾値(Ct)方法は、信号閾値を使用する。Ct値は、蛍光反応が背景信号水準以上に増加して、あらかじめ決定された信号閾値と交差する、増幅曲線の指数期内のポイントに基づいて決定される。このような場合、固定信号閾値(fixed signal threshold、FST)の使用は、Ct値の決定に誤りを発生させるおそれがある。
図11Aは、抽出したFlu AゲノムRNAの希釈(10−3)から得た、補正された増幅曲線を示す。
図11Aに示されたように、FST値が‘200’RFUのように低い値に決定されると、実際標的信号が指数的に増加するポイントではなく、初期ノイズ信号を発生させるポイントがCt値に決定されることがあり、これにによって、標的核酸濃度の決定における誤りが発生するようになる。
図11Bは、抽出されたFlu AゲノムRNAの希釈(10−7〜10−3)から得た、補正された増幅曲線を示す。図11Bに示されたように、FST値が‘500’RFUのように高い値に決定されると、低い濃度の標的配列を含む試料は、実際標的核酸の量より少ない量に決定されることがある。
このように、固定信号閾値を設定する伝統的なCt方法は、標的核酸の濃度決定における誤りを除去することができないことを確認した。
VST(Variable Signal Threshold)方法によるCt値の決定
実施例3で使用したVST(可変信号閾値)方法は、信号閾値 − 変化サイクル(signal threshold−changed cycle、STCC)を決定し、互いに異なる信号ベースライン閾値を、前記STCC前及び後のサイクルにそれぞれ適用する方式で行う。
図11Aにおいて、STCCを10サイクルに決定し、第1のST(信号閾値)を1〜10サイクルに、第2のSTを11〜45サイクルにそれぞれ相異に適用した。図11B及び図11Cにおいて、STCCを38サイクルに決定し、第1のST(信号閾値)を1〜38サイクルに適用して、第2のSTを39〜45サイクルにそれぞれ相異に適用した。前記第1のSTは、‘500’RFUに決定し、第2のSTは、‘200’RFUにそれぞれ決定した。
図11Aに示されたように、FSTを‘200’RFUで適用したとき、FSTと増幅曲線間の最初のクロスポイント(CP)におけるCt値は、1.24に決定された。一方、VSTを適用したとき、VSTに対する最初のクロスポイント(CP)におけるCt値は、32.02に決定された。
図11B及び図11Cに示されたように、標的核酸の存在または不存在を区分するカット − オフ値は、Ct<40に設定した。FSTを‘500’RFUに設定して適用したとき、10−6希釈されたRNAの結果は、Ct 43.41を示したが、この値は、標的核酸の量を不正確に反映したものである。
また、標的核酸の存在または不存在を区分するカット−オフ値がCt<40に設定されたため、Ct 43.41は、標的核酸の不存在を示すが、これは、偽陰性結果に当たる。VSTを適用する場合、10−6希釈されたRNAの結果は、Ct 39.82を示すが、これは、標的核酸の量を正確に反映したものであって、陽性の結果を示す。
したがって、可変的な信号閾値を互いに異なるように設定するVST方法は、初期増幅曲線における初期ノイズ信号を発生させるポイントをCt値として誤って解釈する誤りを除去することができる。また、VST方法は、Ct値をより正確に決定することができ、標的核酸の濃度決定誤りを除去することができる。
上述のように、本発明は、適切な信号閾値の設定を通じて、正確なCt値を決定することにより、さらに信頼性の高い正確な方式で標的核酸を検出して、定量することができる。
以上、本発明の特定な部分を詳細に記述したが、当業界の通常の知識を有する者にとっては、このような具体的な記述はただ望ましい具現例に過ぎなく、これに本発明の範囲が限定されないことは明らかである。従って、本発明の実質的な範囲は、添付の請求項とその等価物により定義されると言える。

Claims (21)

  1. (a)試料を使用する信号発生反応のそれぞれのサイクルにおける信号値を収得し、サイクルにおける信号値を提供する段階と、
    (b)それぞれのサイクルに一つの閾値を適用し、複数の閾値をサイクルに適用する段階であって、前記サイクルの少なくとも二つのサイクルの閾値は、互いに異なっている段階と、
    (c)閾値のそれぞれにより定められる閾値基準を満足する一つ以上のサイクルを確認する段階と、
    (d)前記段階(c)で確認されたサイクルを使用して試料を分析する段階と、
    を含
    前記閾値は、閾値変化サイクル(threshold−changed cycle, TCC)を基準に、前記TCC以前のサイクル及びこれに対応する閾値の集合により形成された関数と、前記TCC以後のサイクル及びこれに対応する閾値の集合により形成された関数とが互いに異なっているように設定される試料の分析方法。
  2. 前記試料の分析は、試料内分析物質の存在または不存在を定量的または定性的に決定するものであることを特徴とする、請求項1に記載の方法。
  3. 前記信号発生反応は、増幅曲線を提供する反応であることを特徴とする、請求項1に記載の方法。
  4. 前記信号値は、信号発生反応から発生した信号値であるか、あるいは信号発生反応から発生した信号の数学的に変形された値であることを特徴とする、請求項1に記載の方法。
  5. 前記閾値基準は、閾値と同一であるか、閾値より大きい信号値を有することを特徴とする、請求項1に記載の方法。
  6. 前記試料の分析は、試料内標的核酸分子の存在を決定するものであり、前記閾値基準を満足する一つ以上のサイクルの確認は、増幅曲線のCt値を決定するものであることを特徴とする、請求項3に記載の方法。
  7. 前記試料の分析は、試料内標的核酸分子の存在を決定するものであり、前記閾値基準を満足する一つ以上のサイクルの確認は、増幅曲線のベースライン領域のエンドポイントサイクルを決定するものであることを特徴とする、請求項3に記載の方法。
  8. (a)(i)増幅反応における増幅サイクル及び(ii)信号発生手段から得た前記増幅サイクルにおける信号値を含む生データセットを収得する段階と、
    (b)前記生データセットからベースラインスタートポイントサイクルとエンドポイントサイクルを決定して、ベースライン領域を決定する段階と、
    (c)前記決定されたベースライン領域内生データセットの二つ以上のデータポイントからベースライン領域の最適ライン関数を算出する段階と、
    (d)前記生データセットの信号値から前記ベースライン領域の最適ライン関数により計算された関数値を差し引き、補正されたデータセットを収得する段階であって、前記補正されたデータセットは、(i)増幅反応の増幅サイクル及び(ii)前記差し引いた結果を含む段階と、
    を含む、信号発生手段の存在下で行われた増幅反応から収得した生データセットを補正する方法であって、
    前記段階(b)において、ベースライン領域のエンドポイントサイクルは、
    (b1)それぞれの増幅サイクルに一つのベースライン閾値を適用し、複数の閾値をサイクルに適用する段階であって、前記サイクルの少なくとも二つのサイクルのベースライン閾値は、互いに異なっている段階と、
    (b2)ベースライン閾値のそれぞれにより決定される閾値基準を満足する一つ以上のサイクルを確認する段階と、
    (b3)前記段階(b2)で確認されたサイクルを使用してベースライン領域のエンドポイントサイクルを決定する段階と、
    を含む方法により決定され、
    増幅サイクルに対するベースライン閾値は、ベースライン閾値変化サイクル(baseline threshold−changed cycle, BTCC)を基準に、前記BTCC以前のサイクル及びこれに対応するベースライン閾値の集合により形成された第1関数と、前記BTCC以後のサイクル及びこれに対応するベースライン閾値の集合により形成された第2関数とが互いに異なっているように設定される、方法。
  9. 前記段階(a)は、前記生データセットをプロットして増幅曲線を収得する段階をさらに含み、前記段階(d)は、前記補正されたデータセットをプロットして補正された増幅曲線を収得する段階をさらに含むことを特徴とする、請求項8に記載の方法。
  10. 前記増幅サイクルは、ベースライン閾値変化サイクルを境界に2以上のグループに分類されて、同一なグループに属するサイクルは、連続的であり、同一なグループに属するサイクルは、互いに同一なベースライン閾値を有して、互いに隣接したグループに属するサイクルは、互いに異なるベースライン閾値を有することを特徴とする、請求項8に記載の方法。
  11. 前記段階(b2)における確認は、生データセットを使用して、それぞれの増幅サイクルに対して計算された傾きと、それぞれの増幅サイクルに対するベースライン閾値を比較して行われることを特徴とする、請求項8に記載の方法。
  12. 前記傾きは、特定サイクルにおけるデータポイント及び前記特定サイクルの以前及び/または以後のサイクルの少なくとも一つのデータポイントを利用した最小二乗法により計算されることを特徴とする、請求項11に記載の方法。
  13. 前記段階(b)において、ベースラインのエンドポイントサイクルは、最小ベースラインエンドポイントサイクル(minimum baseline end−point cycle, MBEC)以上のサイクルで決定されることを特徴とする、請求項8に記載の方法。
  14. 前記ベースラインエンドポイントサイクルは、
    (i)それぞれの増幅サイクルに対して計算された傾きを収得する段階と、
    (ii)前記傾きをそれぞれの増幅サイクルに対するベースライン閾値と比較し、ベースライン領域の候補エンドポイントサイクルを収得する段階と、
    (iii)前記候補エンドポイントサイクルをMBECと比較する段階であって、前記候補エンドポイントサイクルがMBECより大きい場合、前記候補エンドポイントサイクルをエンドポイントサイクルに決定する段階と、
    を含む方法により決定することを特徴とする、請求項13に記載の方法。
  15. 前記ベースライン領域の最適ライン関数を算出する段階は、ベースライン領域内の少なくとも二つのデータポイントを使用した線形回帰分析により行うことを特徴とする、請求項8に記載の方法。
  16. 試料を分析する方法を実行するためのプロセッサを具現する指示を含む、コンピュータ読み取り可能な記録媒体であって、前記試料を分析する方法は、
    (a)前記試料を使用する信号発生反応のそれぞれのサイクルにおける信号値を収得し、サイクルにおける信号値を提供する段階と、
    (b)それぞれのサイクルに一つの閾値を適用し、複数の閾値をサイクルに適用する段階であって、前記サイクルの少なくとも二つのサイクルの閾値は、互いに異なっている段階と、
    (c)閾値のそれぞれにより定められる閾値基準を満足する一つ以上のサイクルを確認する段階と、
    (d)前記段階(c)で確認されたサイクルを使用して試料を分析する段階と、
    を含み、
    前記閾値は、閾値変化サイクル(threshold−changed cycle, TCC)を基準に、前記TCC以前のサイクル及びこれに対応する閾値の集合により形成された関数と、前記TCC以後のサイクル及びこれに対応する閾値の集合により形成された関数とが互いに異なっているように設定されることを特徴とする、記録媒体。
  17. 信号発生手段の存在下で行われた増幅反応から収得した生データセットを補正する方法を実行するためのプロセッサを具現する指示を含む、コンピュータ読み取り可能な記録媒体であって、前記生データセットを補正する方法は、
    (a)(i)前記増幅反応における増幅サイクル及び(ii)前記信号発生手段から得た前記増幅サイクルにおける信号値を含む生データセットを収得する段階と、
    (b)前記生データセットからベースラインスタートポイントサイクルとエンドポイントサイクルを決定して、ベースライン領域を決定する段階と、
    (c)前記決定されたベースライン領域内生データセットの二つ以上のデータポイントからベースライン領域の最適ライン関数を算出する段階と、
    (d)前記生データセットの信号値から前記ベースライン領域の最適ライン関数により計算された関数値を差し引き、補正されたデータセットを収得する段階であって、前記補正されたデータセットは、(i)増幅反応の増幅サイクル及び(ii)前記差し引いた結果を含む段階と、
    を含み、
    前記段階(b)において、ベースライン領域のエンドポイントサイクルは、
    (b1)それぞれの増幅サイクルに一つのベースライン閾値を適用し、複数の閾値をサイクルに適用する段階であって、前記サイクルの少なくとも二つのサイクルのベースライン閾値は、互いに異なっている段階と、
    (b2)ベースライン閾値のそれぞれにより決定される閾値基準を満足する一つ以上のサイクルを確認する段階と、
    (b3)前記段階(b2)で確認されたサイクルを使用してベースライン領域のエンドポイントサイクルを決定する段階と、
    を含む方法により決定され、
    増幅サイクルに対するベースライン閾値は、ベースライン閾値変化サイクル(baseline threshold−changed cycle, BTCC)を基準に、前記BTCC以前のサイクル及びこれに対応するベースライン閾値の集合により形成された第1関数と、前記BTCC以後のサイクル及びこれに対応するベースライン閾値の集合により形成された第2関数とが互いに異なっているように設定されることを特徴とする、記録媒体。
  18. (a)コンピュータプロセッサと、(b)前記コンピュータープロセッサと連結された、請求項16に記載のコンピュータ読み取り可能な記録媒体とを含む、試料の分析のための装置。
  19. (a)コンピュータプロセッサと、(b)前記コンピュータープロセッサと連結された、請求項17に記載のコンピュータ読み取り可能な記録媒体とを含む、信号発生手段の存在下で行われた増幅反応から得た生データセットを補正する方法のための装置。
  20. 試料を分析する方法を実行するためのプロセッサを具現する、コンピュータ読み取り可能な記録媒体に格納されたコンピュータプログラムであって、前記試料を分析する方法は、
    (a)前記試料を使用する信号発生反応のそれぞれのサイクルにおける信号値を収得し、サイクルにおける信号値を提供する段階と、
    (b)それぞれのサイクルに一つの閾値を適用し、複数の閾値をサイクルに適用する段階であって、前記サイクルの少なくとも二つのサイクルの閾値は、互いに異なっている段階と、
    (c)閾値のそれぞれにより定められる閾値基準を満足する一つ以上のサイクルを確認する段階と、
    (d)前記段階(c)で確認されたサイクルを使用して試料を分析する段階と、
    を含み、
    前記閾値は、閾値変化サイクル(threshold−changed cycle, TCC)を基準に、前記TCC以前のサイクル及びこれに対応する閾値の集合により形成された関数と、前記TCC以後のサイクル及びこれに対応する閾値の集合により形成された関数とが互いに異なっているように設定されることを特徴とする、コンピュータプログラム。
  21. 信号発生手段の存在下で行われた増幅反応から収得した生データセットを補正する方法を実行するためのプロセッサを具現する、コンピュータ読み取り可能な記録媒体に格納されたコンピュータプログラムであって、前記生データセットを補正する方法は、
    (a)(i)前記増幅反応における増幅サイクル及び(ii)前記信号発生手段から得た前記増幅サイクルにおける信号値を含む生データセットを収得する段階と、
    (b)前記生データセットからベースラインスタートポイントサイクルとエンドポイントサイクルを決定して、ベースライン領域を決定する段階と、
    (c)前記決定されたベースライン領域内生データセットの二つ以上のデータポイントからベースライン領域の最適ライン関数を算出する段階と、
    (d)前記生データセットの信号値から前記ベースライン領域の最適ライン関数により計算された関数値を差し引き、補正されたデータセットを収得する段階であって、前記補正されたデータセットは、(i)増幅反応の増幅サイクル及び(ii)前記差し引いた結果を含む段階と、
    を含み、
    前記段階(b)において、ベースライン領域のエンドポイントサイクルは、
    (b1)それぞれの増幅サイクルに一つのベースライン閾値を適用し、複数の閾値をサイクルに適用する段階であって、前記サイクルの少なくとも二つのサイクルのベースライン閾値は、互いに異なっている段階と、
    (b2)ベースライン閾値のそれぞれにより決定される閾値基準を満足する一つ以上のサイクルを確認する段階と、
    (b3)前記段階(b2)で確認されたサイクルを使用してベースライン領域のエンドポイントサイクルを決定する段階と、
    を含む方法により決定され、
    増幅サイクルに対するベースライン閾値は、ベースライン閾値変化サイクル(baseline threshold−changed cycle, BTCC)を基準に、前記BTCC以前のサイクル及びこれに対応するベースライン閾値の集合により形成された第1関数と、前記BTCC以後のサイクル及びこれに対応するベースライン閾値の集合により形成された第2関数とが互いに異なっているように設定されることを特徴とする、コンピュータプログラム。
JP2017517690A 2014-10-01 2015-09-30 試料の分析方法 Active JP6602374B2 (ja)

Applications Claiming Priority (7)

Application Number Priority Date Filing Date Title
KR10-2014-0132229 2014-10-01
KR20140132229 2014-10-01
KR20140191924 2014-12-29
KR10-2014-0191924 2014-12-29
KR10-2015-0051080 2015-04-10
KR1020150051080A KR20160039529A (ko) 2014-10-01 2015-04-10 증폭곡선을 보정하는 방법 및 시료 분석을 위한 정보를 획득하는 방법, 그 기록매체, 그 컴퓨터 프로그램, 컴퓨터 장치
PCT/KR2015/010308 WO2016052991A1 (en) 2014-10-01 2015-09-30 Methods for analyzing samples

Publications (2)

Publication Number Publication Date
JP2018502552A JP2018502552A (ja) 2018-02-01
JP6602374B2 true JP6602374B2 (ja) 2019-11-06

Family

ID=55800641

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017517690A Active JP6602374B2 (ja) 2014-10-01 2015-09-30 試料の分析方法

Country Status (4)

Country Link
US (2) US20170226563A1 (ja)
EP (1) EP3201359B1 (ja)
JP (1) JP6602374B2 (ja)
KR (2) KR20160039529A (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017209563A1 (en) 2016-06-02 2017-12-07 Seegene, Inc. . Method for detecting a target analyte in a sample using a signal change-amount data set
EP3690886A4 (en) * 2017-09-28 2021-06-23 Seegene, Inc. AMPLIFICATION DATA DISPLAY METHOD AND DEVICE
US20210317515A1 (en) * 2018-07-10 2021-10-14 Gen-Probe Incorporated Methods and systems for detecting and quantifying nucleic acids
US11150165B2 (en) * 2019-03-01 2021-10-19 Dell Products, L.P. System and method for configuration drift detection and remediation
EP4421814A1 (en) * 2021-10-21 2024-08-28 Seegene, Inc. Device and method for reading positive/negative with respect to target assay substance in sample
WO2024025390A1 (ko) * 2022-07-29 2024-02-01 주식회사 씨젠 파라미터의 임계값 설정용 인터페이스 제공 장치 및 방법
CN115985396B (zh) * 2022-12-16 2023-12-12 苏州思迈德生物科技有限公司 一种实时荧光定量pcr扩增数据的分析处理方法及装置
WO2024147568A1 (ko) * 2023-01-05 2024-07-11 주식회사 씨젠 분자진단 분석 결과의 획득 방법, 분자진단 분석 결과를 추정하는 모델의 획득 방법 및 이를 수행하는 컴퓨터 장치
CN116564411B (zh) * 2023-07-12 2023-09-29 广州普世君安生物科技有限公司 等温扩增核酸的结果判定方法、系统、设备及存储介质

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7188030B2 (en) 2001-08-21 2007-03-06 Applera Corporation Automatic threshold setting for quantitative polymerase chain reaction
US7363168B2 (en) * 2001-10-02 2008-04-22 Stratagene California Adaptive baseline algorithm for quantitative PCR
US7228237B2 (en) 2002-02-07 2007-06-05 Applera Corporation Automatic threshold setting and baseline determination for real-time PCR
US7720611B2 (en) 2005-05-13 2010-05-18 Bio-Rad Laboratories, Inc. Baselining amplification data
US7844403B2 (en) * 2005-12-20 2010-11-30 Roche Molecular Systems, Inc. Temperature step correction with double sigmoid Levenberg-Marquardt and robust linear regression
EP1804172B1 (en) * 2005-12-20 2021-08-11 Roche Diagnostics GmbH PCR elbow determination using curvature analysis of a double sigmoid
WO2011130184A2 (en) * 2010-04-11 2011-10-20 Life Technologies Corporation SYSTEMS AND METHODS FOR MODEL-BASED qPCR
KR20120103496A (ko) * 2011-03-09 2012-09-19 삼성테크윈 주식회사 실시간 중합효소 연쇄 반응에 있어서 표준 곡선 및 임계값을 획득하는 방법

Also Published As

Publication number Publication date
KR20160039529A (ko) 2016-04-11
EP3201359A1 (en) 2017-08-09
KR102336732B1 (ko) 2021-12-08
US20210238652A1 (en) 2021-08-05
EP3201359A4 (en) 2018-05-02
JP2018502552A (ja) 2018-02-01
KR20170051539A (ko) 2017-05-11
EP3201359B1 (en) 2019-12-04
US20170226563A1 (en) 2017-08-10

Similar Documents

Publication Publication Date Title
JP6602374B2 (ja) 試料の分析方法
KR102385959B1 (ko) 신호 변화량 데이터 세트를 이용한 샘플 내 타겟 분석물질 검출 방법
JP6431076B2 (ja) リアルタイムpcr信号におけるジャンプの検出及び補正
JP6685396B2 (ja) 標的分析物質に対するデータセットの補正方法
JP7379548B2 (ja) 試料内ターゲット分析物を分析する方法及び装置
US10978173B2 (en) Method for reducing noise level of data set for a target analyte
WO2016052991A1 (en) Methods for analyzing samples
US11781180B2 (en) Detection of abnormal signal using two or more datasets
KR102110985B1 (ko) 타겟 핵산 서열에 대한 시그널 추출
KR102189356B1 (ko) 데이터 세트에서 비정상적인 시그널의 검출
KR102270849B1 (ko) 타겟 핵산 서열의 존재를 결정하기 위한 분석 시그널

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170417

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170417

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180417

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20180711

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20180912

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181017

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20190207

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190319

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20190618

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190712

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20190917

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20191008

R150 Certificate of patent or registration of utility model

Ref document number: 6602374

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250