JP6431076B2 - リアルタイムpcr信号におけるジャンプの検出及び補正 - Google Patents
リアルタイムpcr信号におけるジャンプの検出及び補正 Download PDFInfo
- Publication number
- JP6431076B2 JP6431076B2 JP2016543190A JP2016543190A JP6431076B2 JP 6431076 B2 JP6431076 B2 JP 6431076B2 JP 2016543190 A JP2016543190 A JP 2016543190A JP 2016543190 A JP2016543190 A JP 2016543190A JP 6431076 B2 JP6431076 B2 JP 6431076B2
- Authority
- JP
- Japan
- Prior art keywords
- value
- data set
- jump
- amplification
- derivative
- 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
Links
- 238000012937 correction Methods 0.000 title claims description 37
- 238000003753 real-time PCR Methods 0.000 title description 14
- 238000001514 detection method Methods 0.000 title description 10
- 230000003321 amplification Effects 0.000 claims description 171
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 171
- 238000000034 method Methods 0.000 claims description 151
- 230000008569 process Effects 0.000 claims description 72
- 238000003752 polymerase chain reaction Methods 0.000 claims description 60
- 238000012408 PCR amplification Methods 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 5
- 230000006870 function Effects 0.000 description 33
- 108020004414 DNA Proteins 0.000 description 14
- 150000007523 nucleic acids Chemical group 0.000 description 13
- 238000004458 analytical method Methods 0.000 description 12
- 108020004707 nucleic acids Proteins 0.000 description 12
- 102000039446 nucleic acids Human genes 0.000 description 12
- 239000000523 sample Substances 0.000 description 12
- 239000000047 product Substances 0.000 description 11
- 238000006243 chemical reaction Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 6
- 239000007850 fluorescent dye Substances 0.000 description 5
- 239000013615 primer Substances 0.000 description 5
- 238000011002 quantification Methods 0.000 description 5
- 230000007704 transition Effects 0.000 description 5
- 102000053602 DNA Human genes 0.000 description 4
- 230000001580 bacterial effect Effects 0.000 description 4
- 150000001875 compounds Chemical class 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 108091033319 polynucleotide Proteins 0.000 description 4
- 102000040430 polynucleotide Human genes 0.000 description 4
- 239000002157 polynucleotide Substances 0.000 description 4
- 108091028043 Nucleic acid sequence Proteins 0.000 description 3
- 108020005187 Oligonucleotide Probes Proteins 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 3
- 230000027455 binding Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000005415 bioluminescence Methods 0.000 description 3
- 230000029918 bioluminescence Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 239000000975 dye Substances 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000004020 luminiscence type Methods 0.000 description 3
- 239000002751 oligonucleotide probe Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 238000010222 PCR analysis Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 239000012472 biological sample Substances 0.000 description 2
- 238000011088 calibration curve Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 239000013039 cover film Substances 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000002255 enzymatic effect Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000005499 meniscus Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 239000007858 starting material Substances 0.000 description 2
- 230000004544 DNA amplification Effects 0.000 description 1
- 239000003155 DNA primer Substances 0.000 description 1
- 230000004568 DNA-binding Effects 0.000 description 1
- 102000016928 DNA-directed DNA polymerase Human genes 0.000 description 1
- 108010014303 DNA-directed DNA polymerase Proteins 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108091034117 Oligonucleotide Proteins 0.000 description 1
- 241000239226 Scorpiones Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 238000000149 argon plasma sintering Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000006227 byproduct Substances 0.000 description 1
- 230000003197 catalytic effect Effects 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000004925 denaturation Methods 0.000 description 1
- 230000036425 denaturation Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 239000000539 dimer Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000002795 fluorescence method Methods 0.000 description 1
- 238000002866 fluorescence resonance energy transfer Methods 0.000 description 1
- 238000001506 fluorescence spectroscopy Methods 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 238000000338 in vitro Methods 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000013642 negative control Substances 0.000 description 1
- 230000009871 nonspecific binding Effects 0.000 description 1
- 239000002773 nucleotide Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 102000054765 polymorphisms of proteins Human genes 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000013207 serial dilution Methods 0.000 description 1
- 238000013518 transcription Methods 0.000 description 1
- 230000035897 transcription Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000003612 virological effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B50/00—ICT programming tools or database systems specially adapted for bioinformatics
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6844—Nucleic acid amplification reactions
- C12Q1/6851—Quantitative amplification
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6844—Nucleic acid amplification reactions
- C12Q1/686—Polymerase chain reaction [PCR]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
- G16B40/10—Signal processing, e.g. from mass spectrometry [MS] or from PCR
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Organic Chemistry (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioethics (AREA)
- Molecular Biology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Data Mining & Analysis (AREA)
- Epidemiology (AREA)
- Software Systems (AREA)
- Public Health (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Biochemistry (AREA)
- Genetics & Genomics (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Microbiology (AREA)
- Signal Processing (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Apparatus Associated With Microorganisms And Enzymes (AREA)
Description
図1Aには、動的PCRプロセスの環境における増幅曲線10の一例が示されている。図示のように、曲線10は、ラグフェーズ領域15と、指数フェーズ領域25と、を含む。ラグフェーズ領域15は、ベースライン又はベースライン領域と一般に呼称される。このような曲線は、ラグフェーズと指数フェーズをリンクする遷移領域20を含む。領域20は、エルボー又はエルボー領域と一般に呼称される。エルボー領域は、通常、ベースラインの末尾と、基礎をなすプロセスの成長又は増幅レートにおける遷移と、を定義している。領域20内において特定の遷移ポイントを識別することは、基礎をなすプロセスの振る舞いを分析するために有用であろう。
ジャンプエラーは、信号強度の非常に大きな増大として表れる。実施形態を使用することにより、ジャンプエラーを検出するプロセスを自動化することができる。従って、実施形態は、コンピュータシステムがジャンプエラーを自動的に検出できるようにするための技法を表すことができる。
実施形態は、増幅曲線の二次導関数を使用してジャンプエラーを識別することができる。具体的には、実施形態は、二次導関数値がそれぞれのサイクル数において取得されるように、二次導関数の離散計算を実行することができる。二次導関数は、データセットとの間の二次導関数カーネルの畳込みを介して算出することができる。二次導関数値を分析することにより、ジャンプエラーを表す振る舞いを識別することができる。
図3は、本発明の実施形態による増幅プロセスの増幅曲線においてジャンプエラーを検出する方法300のフローチャートである。方法300は、コンピュータシステムによって部分的に又は全体的に実行することができる。方法300の特定のブロックは、省略することができる。例示用の増幅プロセスは、リアルタイムPCR増幅である。その他の増幅プロセスは、バクテリアプロセス、酵素プロセス、及び結合プロセスを含む。増幅プロセスは、一連のポイントによって計測されることが可能であり、それぞれのデータポイントは、サイクル数における信号強度を提供する。
上述のように、二次導関数は、カーネルとデータセットの畳込みとして算出することができる。畳込みは、隣接するサイクルにおける信号値の間の有限な差として実装することができる。当業者には明らかなように、且つ、後述するように、様々なカーネルを利用することができる。
図4Aは、その他の点においてはフラットな信号においてジャンプを有する実際のリアルタイムPCRデータを示している。図4Bは、図4Aからのデータの二次導関数を示している。図4Bにおけるデータは、中心有限差分法における畳込みを使用して判定されている。図4Aにおいては、データは、ジャンプの後においてわずかに減少しているが、全体的には、相対的にフラットな状態に留まっている。従って、このデータセットは、増幅を示してはおらず、その理由は、信号における唯一の大きな増大がジャンプに起因するものであるからである。ジャンプが識別されなければ、データセットが、増幅を示しているとして、誤って識別されうるという点において、偽陽性が発生しうるであろう。このような不正確な分析は、生物学的サンプルにおける基礎をなす増幅プロセスの不正確な計測を結果的にもたらしうるであろう。図4Bにおいては、観察することができるように、正の二次導関数及び負の二次導関数を有する連続的なサイクルのペア(9及び10のX値として示されており、これらは、サイクル10及び11に対応している)が存在している。又、これらの値は、その他のサイクルにおける値と比較した場合に、最大の大きさを有している。
データセットは、複数のジャンプエラーを有することができよう。このような状況においては、連続的なサイクルの1つのペアが、その他のジャンプエラーのペアよりも、大きな大きさを有することになろう。すべてのジャンプエラーを検出すると共に、恐らくは、それらについて補正するべく、データを通じて複数のパスを実施することができよう。スプリアスノイズが、非常に小さな正の二次導関数値が次のサイクルにおける非常に小さな負の二次導関数値によって後続されている状態をジャンプエラーとして識別させないことを保証するべく、大きさの最小閾値を使用することができよう。
図示の例は、正のジャンプについてのものであるが、負のジャンプも出現しうるであろう。負のジャンプの場合に、ペアの第1サイクルが、負の二次導関数値を有することになり、且つ、ペアの第2サイクルが、正の二次導関数値を有することになろう。類似の方式により、分析を依然として実装することができる。
ジャンプが検出されたら、ジャンプ及びデータセットを分析することにより、データセットを無効化するのか又はデータセットを補正するのか、或いは、恐らくは、これらのいずれをも実施しないのか、を判定することができる。様々な基準を使用することにより、これらの決定を下すことができる。例えば、1つ又は複数の無効化規準の組を使用することにより、データセットを無効化するかどうかを判定することができる。1つ又は複数の補正基準の組を使用することにより、データセットを補正するかどうかを判定することができる。これらの基準は、データを分析するコンピュータシステムに入力することができる。無効化するのか又は補正するのかの判定は、データの収集時点において、或いは、後の事後処理ステップにおいて、実行することができる。
上述のように、様々な基準を使用することにより、ジャンプを示す増幅データを無効化するのか又は補正するのかを判定することができる。増幅曲線を無効化又は補正するための決定プロセスにおいては、例えば、ジャンプの高さの大きさ又は正味の増幅との関係における大きさ(例えば、合計増幅−ジャンプの高さ)、ベースライン切片との関係における大きさ、及びジャンプのサイクル数を使用してもよいであろう。
ジャンプエラーは、様々な方法で補正することができる。ジャンプについてデータを補正する1つの方法は、ジャンプの後のすべてのデータポイントからジャンプの高さを減算するというものである。別の方法は、ジャンプの前のすべてのデータポイントにジャンプの高さを加算するというものである。ジャンプ補正の後の選択肢は、補正済みのジャンプの後のすべての蛍光値を1つのサイクル数だけ左にシフトさせるというものである。これにより、なんらかの不連続性が補正済みのジャンプ曲線内に出現することを防止することができる。このようなシフトは、ジャンプが指数フェーズ(増幅領域)において発生した際に有益であろう。
図9〜図12は、本発明の実施形態による、ジャンプを有する増幅曲線と、補正済みの曲線と、を示している。これらの増幅曲線は、図4A、図5A、図6A、及び図7Aに対応している。その他のプロットと同様に、垂直軸は、強度であり、且つ、水平軸は、サイクル数である。
本明細書において言及されているコンピュータシステムの任意のものは、任意の適切な数のサブシステムを利用しうる。図14には、このようなサブシステムの例が、コンピュータ装置1400内において示されている。いくつかの実施形態においては、コンピュータシステムは、単一のコンピュータ装置を含み、この場合には、サブシステムは、コンピュータ装置のコンポーネントであってもよい。その他の実施形態においては、コンピュータシステムは、内部コンポーネントを有する複数のコンピュータ装置を含むことが可能であり、この場合には、コンピュータ装置のそれぞれがサブシステムである。
又、特定の態様においては、本発明は、PCRシステムをも提供している。図15及び図16には、例示用のPCRシステムが表示されている。図15は、本発明の方法及びシステムを実装するべく使用されうるソフトウェア及びハードウェアリソースの間における関係を説明する概略ブロックダイアグラムを示している。図16において示されているシステムは、サーモサイクラ装置内において配置されうる動的PCR分析モジュールと、コンピュータシステムの一部分であるインテリジェントモジュールと、を有する。データセット(PCRデータセット)が、ネットワーク接続又は直接接続を介して、分析モジュールからインテリジェントモジュールに、又はこの逆に、転送される。データセットは、例えば、図3に示されているフローチャートに従って処理されてもよい。このフローチャートは、例えば、図15に示されているフローチャートに従って、コンピュータシステムのハードウェア上に保存されたソフトウェアにより、便利な方式で実装されてもよい。図15を参照すれば、コンピュータシステム(2000)は、例えば、PCR反応の際に取得された蛍光データを受け取る受取り手段(2100)と、データセットとの間のカーネルの離散畳込みを算出してサイクル数におけるデータセットの二次導関数を取得することにより、前記データを処理する計算手段(2200)と、2つの連続的なサイクル数が、異なる符号を有する二次導関数の値を有するかどうかを識別する第1識別手段(2300)と、第1サイクル数及び第2サイクル数が、それぞれ、異なる符号を有する二次導関数の第1値及び第2値を有している際に、第1値及び第2値のうちの少なくとも1つの値の大きさが閾値を上回っているかどうかを判定する判定手段(2400)と、第1値及び第2値のうちの少なくとも1つの値の大きさが閾値を上回っているという判定に基づいて、増幅曲線を表すデータセットにおいてジャンプエラーが発生したと識別する第2識別手段(2500)と、を有してもよい。又、特定の実施形態においては、システムは、コンピュータ画面上において結果を表示するための表示手段を有してもよい。図16は、サーモサイクラ装置とコンピュータシステムの間におけるやり取りを示している。システムは、サーモサイクラ内において配置されうる動的PCR分析モジュールと、コンピュータシステムの一部分であるインテリジェントモジュールと、を有する。データセット(PCRデータセット)が、ネットワーク接続又は直接的接続を介して、分析モジュールからインテリジェントモジュールに、又はこの逆に、転送される。データセットは、プロセッサ上において稼働すると共にインテリジェントモジュールのストレージ装置上において保存されているコンピュータコードにより、図15に従って処理されてもよく、且つ、処理の後に、分析装置のストレージ装置に転送して戻されてもよく、そこで、変更済みのデータは、表示装置上において表示されてもよい。又、いくつかの実施形態においては、インテリジェントモジュールは、サーモサイクラ上において実装されてもよい。特定の態様においては、本発明は、ポリメラーゼ連鎖反応(PCR)システムを提供し、システムは、PCR増幅曲線を表すPCRデータセットを生成するPCRデータ取得装置であって、前記データセットは、複数のデータポイントを含み、データポイントのそれぞれは、サイクル数とそのサイクル数における増幅プロセスの信号強度の座標値のペアを有する、装置と、データセットとの間のカーネルの離散畳込みを算出してサイクル数におけるデータセットの二次導関数を取得し、2つの連続的なサイクル数が、異なる符号を有する二次導関数の値を有するかどうかを識別し、第1サイクル数及び第2サイクル数が、それぞれ、異なる符号を有する二次導関数の第1値及び第2値を有している際に、第1値及び第2値のうちの少なくとも1つの値の大きさが閾値を上回っているかどうかを判定し、第1値及び第2値のうちの少なくとも1つの値の大きさが閾値を上回っているという判定に基づいて、増幅曲線を表すデータセットにおいてジャンプエラーが発生したと識別することにより、PCRデータセットを処理してジャンプエラーを検出するように適合されたインテリジェントモジュールと、を有する。
Claims (20)
- 増幅プロセスの増幅曲線においてジャンプエラーを検出する方法であって、
−コンピュータシステムにおいて、増幅曲線を表すデータセットを受け取るステップであって、前記データセットは、複数のデータポイントを含み、それぞれのデータポイントは、サイクル数とそのサイクル数における前記増幅プロセスの信号強度の座標値のペアを有する、ステップと、
−前記コンピュータシステムにより、前記データセットとの間のカーネルの離散畳込みを算出して前記サイクル数における前記データセットの二次導関数を取得するステップと、
−2つの連続的なサイクル数が、異なる符号を有する前記二次導関数の値を有するかどうかを識別するステップと、
−第1サイクル数及び第2サイクル数が、それぞれ、異なる符号を有する前記二次導関数の第1の値及び第2の値を有している際に、前記第1値及び前記第2値のうちの少なくとも1つの値の大きさが閾値を上回っているかどうかを判定するステップと、
−前記第1値及び前記第2値のうちの少なくとも1つの値の大きさが前記閾値を上回っているという判定に基づいて、前記増幅曲線を表す前記データセットにおいてジャンプエラーが発生したと識別するステップと、
を有する方法。 - −前記第1値及び前記第2値のうちの両方の値の大きさが1つ又は複数の閾値を上回っているかどうかを判定するステップを更に有し、
ジャンプエラーが発生したという前記識別は、前記第1値及び前記第2値のうちの両方の値の大きさが前記1つ又は複数の閾値を上回っているという判定に基づいている請求項1に記載の方法。 - 前記閾値は、前記第1値及び前記第2値が、前記二次導関数のその他の値との関係において最大の大きさを有するかどうかを試験する請求項2に記載の方法。
- 前記閾値は、既定の値である請求項1乃至3のいずれか一項に記載の方法。
- 前記第1値は、正であり、且つ、前記第2値は、負である請求項1乃至4のいずれか一項に記載の方法。
- −前記ジャンプエラーの1つ又は複数の特性を判定するステップと、
−前記ジャンプエラーの前記1つ又は複数の特性が1つ又は複数の無効化規準を充足している際に、前記データセットを無効化するステップと、
を更に有する請求項1乃至5のいずれか一項に記載の方法。 - 前記1つ又は複数の無効化規準は、
−前記第1値及び/又は前記第2値の大きさと、
−前記ジャンプエラーの場所と、
のうちの少なくとも1つを含む請求項6に記載の方法。 - 前記第1値及び/又は前記第2値の前記大きさは、前記データセット内のベースライン又は正味の増幅との関係における大きさである請求項7に記載の方法。
- 前記ジャンプエラーの前記場所は、前記増幅曲線の領域の場所との関係における場所である請求項7に記載の方法。
- −前記ジャンプエラーについて前記データセットを補正するステップを更に有する請求項1乃至9のいずれか一項に記載の方法。
- 前記ジャンプエラーの1つ又は複数の特性を判定するステップであって、前記データセットは、前記ジャンプエラーの前記1つ又は複数の特性が1つ又は複数の補正基準を充足している際に、補正される、ステップを更に有する請求項10に記載の方法。
- 前記1つ又は複数の補正基準は、
前記第1値及び/又は前記第2値の大きさと、
−前記ジャンプエラーのサイクル数と、
のうちの少なくとも1つを含む請求項11に記載の方法。 - −前記補正済みのデータセットに関数をフィッティングするステップと、
−前記フィッティングされた関数を分析して前記増幅プロセスの特性を判定するステップと、
を更に有する請求項10に記載の方法。 - 前記データセットを補正するステップは、
−前記ジャンプエラーの場所Xを識別するステップと、
−前記ジャンプエラーの高さを判定するステップと、
−前記大きさを使用して前記ジャンプエラーの前のサイクル数の前記信号強度を補正するか又は前記ジャンプエラーの後のサイクル数について前記データセットを補正するステップと、
を含む請求項10に記載の方法。 - 前記データセットを補正するステップは、
−1つだけ少ないサイクル数に対応するように、X+2を上回るサイクル数の前記信号強度の値をシフトさせるステップを更に含む請求項14に記載の方法。 - −ジャンプエラーを識別した後に、2つの連続的なサイクル数が、異なる符号を有する前記二次導関数の値を有するその他の場所について前記データセットをサーチするステップを更に有する請求項1乃至15のいずれか一項に記載の方法。
- 実行された際に、増幅プロセスの増幅曲線においてジャンプエラーを検出するようにコンピュータシステムを制御する複数の命令を保存した一時的ではないコンピュータ可読媒体を有するコンピュータプロダクトであって、
前記命令は、
−増幅曲線を表すデータセットを受け取るステップであって、前記データセットは、複数のデータポイントを含み、それぞれのデータポイントは、サイクル数とそのサイクル数における前記増幅プロセスの信号強度の座標値のペアを有する、ステップと、
−前記データセットとの間のカーネルの離散畳込みを算出して前記サイクル数における前記データセットの二次導関数を取得するステップと、
−2つの連続的なサイクル数が、異なる符号を有する前記二次導関数の値を有するかどうかを識別するステップと、
−第1サイクル数及び第2サイクル数が、それぞれ、異なる符号を有する前記二次導関数の第1値及び第2値を有している際に、前記第1値及び前記第2値のうちの少なくとも1つの値の大きさが閾値を上回っているかどうかを判定するステップと、
−前記第1値及び前記第2値のうちの少なくとも1つの値の大きさが前記閾値を上回っているという判定に基づいて、前記増幅曲線を表す前記データセットにおいてジャンプエラーが発生したと識別するステップと、
を有する、コンピュータプロダクト。 - 前記命令は、
−前記ジャンプエラーの1つ又は複数の特性を判定するステップと、
−前記ジャンプエラーの前記1つ又は複数の特性が1つ又は複数の無効化規準を充足している際に、前記データセットを無効化するステップと、
前記ジャンプエラーの前記1つ又は複数の特性が1つ又は複数の補正基準を充足している際に、前記データセットを補正するステップと、
を更に有する請求項17に記載のコンピュータプロダクト。 - ポリメラーゼ連鎖反応(PCR)システムであって、
−PCR増幅曲線を表すPCRデータセットを生成するPCRデータ取得装置であって、前記データセットは、複数のデータポイントを含み、前記データポイントのそれぞれは、サイクル数とそのサイクル数における前記増幅プロセスの信号強度の座標値のペアを有する、装置と、
−前記データセットとの間のカーネルの離散畳込みを算出して前記サイクル数における前記データセットの二次導関数を取得し、
−2つの連続的なサイクル数が、異なる符号を有する前記二次導関数の値を有するかどうかを識別し、
−第1サイクル数及び第2サイクル数が、それぞれ、異なる符号を有する前記二次導関数の第1値及び第2値を有している際に、前記第1値及び前記第2値のうちの少なくとも1つの値の大きさが閾値を上回っているかどうかを判定し、
−前記第1値及び前記第2値のうちの少なくとも1つの値の大きさが前記閾値を上回っているという判定に基づいて、前記増幅曲線を表す前記データセットにおいてジャンプエラーが発生したと識別することにより、
−前記PCRデータセットを処理してジャンプエラーを検出するように適合されたインテリジェントモジュールと、
を有するPCRシステム。 - 前記インテリジェントモジュールは、
−前記ジャンプエラーの1つ又は複数の特性を判定し、且つ、
−前記ジャンプエラーの前記1つ又は複数の特性が1つ又は複数の無効化規準を充足している際に、前記データセットを無効化し、且つ、
−前記ジャンプエラーの前記1つ又は複数の特性が1つ又は複数の補正基準を充足している際に、前記データセットを補正する、
ように、更に適合されている請求項19に記載のPCRシステム。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/144,382 US9607128B2 (en) | 2013-12-30 | 2013-12-30 | Detection and correction of jumps in real-time PCR signals |
US14/144,382 | 2013-12-30 | ||
PCT/EP2014/079356 WO2015101591A1 (en) | 2013-12-30 | 2014-12-29 | Detection and correction of jumps in real-time pcr signals |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2017509040A JP2017509040A (ja) | 2017-03-30 |
JP6431076B2 true JP6431076B2 (ja) | 2018-11-28 |
Family
ID=52232210
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016543190A Active JP6431076B2 (ja) | 2013-12-30 | 2014-12-29 | リアルタイムpcr信号におけるジャンプの検出及び補正 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9607128B2 (ja) |
EP (1) | EP3090372B1 (ja) |
JP (1) | JP6431076B2 (ja) |
CN (1) | CN105849736B (ja) |
WO (1) | WO2015101591A1 (ja) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR102165931B1 (ko) * | 2015-09-24 | 2020-10-14 | 주식회사 씨젠 | 타겟 분석물질의 존재 또는 부존재를 결정하기 위한 다중 데이터 세트 분석법 |
WO2017171469A1 (en) * | 2016-04-01 | 2017-10-05 | Seegene, Inc. . | Detection of abnormal signal using two or more datasets |
WO2017196112A1 (en) * | 2016-05-11 | 2017-11-16 | Seegene, Inc. | Detection of abnormal signal in dataset |
CN105755159A (zh) * | 2016-05-13 | 2016-07-13 | 苏州国科闻普生物科技有限公司 | 自动快速判别荧光定量pcr结果的方法 |
US20190156920A1 (en) * | 2016-06-02 | 2019-05-23 | Seegene, Inc. | Method for detecting a target analyte in a sample using a signal change-amount data set |
JP2019004814A (ja) * | 2017-06-27 | 2019-01-17 | 東洋紡株式会社 | 解析用プログラムおよび解析装置 |
US11169899B2 (en) | 2019-04-15 | 2021-11-09 | Toyota Motor Engineering & Manufacturing North America, Inc. | Mitigating data offsets for machine learning |
CN113851189B (zh) * | 2021-09-24 | 2022-09-23 | 中元汇吉生物技术股份有限公司 | 一种自适应检测pcr荧光基线的方法、装置及其设备 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US2997704A (en) | 1958-02-24 | 1961-08-22 | Epsco Inc | Signal conversion apparatus |
US7680868B2 (en) * | 2005-12-20 | 2010-03-16 | Roche Molecular Systems, Inc. | PCR elbow determination by use of a double sigmoid function curve fit with the Levenburg-Marquardt algorithm and normalization |
US7668663B2 (en) | 2005-12-20 | 2010-02-23 | Roche Molecular Systems, Inc. | Levenberg-Marquardt outlier spike removal method |
US7991562B2 (en) | 2005-12-20 | 2011-08-02 | Roche Molecular Systems, Inc. | PCR elbow determination using curvature analysis of a double sigmoid |
US7587283B2 (en) * | 2006-07-01 | 2009-09-08 | Celera Corporation | Growth and baseline identification determination for amplification data |
WO2008070328A2 (en) * | 2006-10-24 | 2008-06-12 | Celera Corporation | Systems and methods for baselining and real-time pcr data analysis |
US20090119020A1 (en) | 2007-09-25 | 2009-05-07 | Roche Molecular Systems, Inc. | Pcr elbow determination using quadratic test for curvature analysis of a double sigmoid |
US8374795B2 (en) | 2008-05-13 | 2013-02-12 | Roche Molecular Systems, Inc. | Systems and methods for step discontinuity removal in real-time PCR fluorescence data |
CA2736243A1 (en) * | 2008-09-09 | 2010-03-18 | Bio-Rad Laboratories, Inc. | Multi-stage, regression-based pcr analysis system |
US8265879B2 (en) * | 2009-04-17 | 2012-09-11 | Roche Molecular Systems, Inc. | Determination of single peak melting temperature by PCR analogy and double sigmoid equation |
US8219366B2 (en) * | 2009-08-26 | 2012-07-10 | Roche Molecular Sytems, Inc. | Determination of elbow values for PCR for parabolic shaped curves |
-
2013
- 2013-12-30 US US14/144,382 patent/US9607128B2/en active Active
-
2014
- 2014-12-29 WO PCT/EP2014/079356 patent/WO2015101591A1/en active Application Filing
- 2014-12-29 CN CN201480071584.XA patent/CN105849736B/zh active Active
- 2014-12-29 JP JP2016543190A patent/JP6431076B2/ja active Active
- 2014-12-29 EP EP14820896.0A patent/EP3090372B1/en active Active
Also Published As
Publication number | Publication date |
---|---|
EP3090372B1 (en) | 2023-06-28 |
WO2015101591A1 (en) | 2015-07-09 |
CN105849736A (zh) | 2016-08-10 |
US20150186598A1 (en) | 2015-07-02 |
CN105849736B (zh) | 2019-08-13 |
JP2017509040A (ja) | 2017-03-30 |
EP3090372A1 (en) | 2016-11-09 |
US9607128B2 (en) | 2017-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6431076B2 (ja) | リアルタイムpcr信号におけるジャンプの検出及び補正 | |
US11615863B2 (en) | Universal method to determine real-time PCR cycle threshold values | |
US7680868B2 (en) | PCR elbow determination by use of a double sigmoid function curve fit with the Levenburg-Marquardt algorithm and normalization | |
JP5599986B2 (ja) | リアルタイムpcr蛍光データ内のステップ型不連続部除去用のシステム及び方法 | |
JP6602374B2 (ja) | 試料の分析方法 | |
JP2010081937A (ja) | 方程式のないアルゴリズムによるリアルタイムpcrエルボーの呼び出し | |
KR102385959B1 (ko) | 신호 변화량 데이터 세트를 이용한 샘플 내 타겟 분석물질 검출 방법 | |
KR102165933B1 (ko) | 둘 이상의 데이터 세트를 이용한 비정상적인 시그널의 검출 | |
EP2471007B1 (en) | Determination of elbow values for pcr for parabolic shaped curves | |
KR101771402B1 (ko) | 핵산 정량 방법 | |
KR102189356B1 (ko) | 데이터 세트에서 비정상적인 시그널의 검출 | |
US20180314789A1 (en) | Melt curve analysis | |
JP2017533719A (ja) | 出発試料の特性を決定する方法 | |
KR20120097215A (ko) | 실시간 중합효소 연쇄 반응의 증폭 프로파일 곡선에 있어서 베이스라인을 결정하는 방법 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20171124 |
|
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: 20181002 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20181101 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6431076 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 |