US20020042051A1 - Automated analysis of real-time nucleic acid amplification - Google Patents

Automated analysis of real-time nucleic acid amplification Download PDF

Info

Publication number
US20020042051A1
US20020042051A1 US09/391,811 US39181199A US2002042051A1 US 20020042051 A1 US20020042051 A1 US 20020042051A1 US 39181199 A US39181199 A US 39181199A US 2002042051 A1 US2002042051 A1 US 2002042051A1
Authority
US
United States
Prior art keywords
fluorescence
amplification
cycle
nucleic acid
cycles
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.)
Granted
Application number
US09/391,811
Other versions
US6387621B1 (en
Inventor
Carl T. Wittwer
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Utah Research Foundation UURF
Original Assignee
University of Utah Research Foundation UURF
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 University of Utah Research Foundation UURF filed Critical University of Utah Research Foundation UURF
Priority to US09/391,811 priority Critical patent/US6387621B1/en
Priority to AT00303519T priority patent/ATE366821T1/en
Priority to EP00303519A priority patent/EP1059523B1/en
Priority to DE60035459T priority patent/DE60035459T2/en
Priority to JP2000128560A priority patent/JP4689004B2/en
Assigned to UTAH RESEARCH FOUNDATION, UNIVERSITY OF reassignment UTAH RESEARCH FOUNDATION, UNIVERSITY OF ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: WITTWER, CARL T.
Priority to US10/074,169 priority patent/US20020123062A1/en
Publication of US20020042051A1 publication Critical patent/US20020042051A1/en
Application granted granted Critical
Publication of US6387621B1 publication Critical patent/US6387621B1/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/6428Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
    • 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]

Definitions

  • the present invention relates to a method of analyzing a sample for the presence of a nucleic acid. More particularly, the present invention is directed to an automated method for detecting and reporting the presence of a predetermined nucleic acid in a sample using polymerase chain reaction and a fluorescent detecting entity.
  • PCR polymerase chain reaction
  • PCR is an invaluable molecular biology tool
  • the practical implementation of real time PCR techniques has lagged behind the conceptual promise.
  • Currently available instrumentation does not actually analyze data during PCR; it simply acquires the data for later analysis. After PCR has been completed, multiple manual steps are necessary to analyze the acquired data, and human judgment is typically required to provide the analysis result.
  • What is needed is a system for automating data acquisition and analysis so that no user intervention is required for reporting the analytical results.
  • the system software is automatically triggered and the results, for example, the presence or absence of a given pathogen, is immediately displayed on screen. Algorithms for detection, quantification, and genotyping are needed.
  • initiation of the analysis algorithm can be implemented prior to completion of temperature cycling.
  • Data processing can occur during amplification and concomitant analysis results can be used to modify temperature cycling and to acquire additional data during the latter stages of the amplification procedure to optimize amplification protocol and data quality.
  • a major problem in automating PCR data analysis is identification of baseline fluorescence. Background fluorescence varies from reaction to reaction. Moreover, baseline drift, wherein fluorescence increases or decreases without relation to amplification of nucleic acids in the sample, is a common occurrence.
  • Prior attempts to automate amplification data analysis involved setting the baseline fluorescence as that measured at one or more predetermined early cycle numbers. This technique accounts for the variation in background fluorescence, but it does not compensate for baseline drift. Without compensation for baseline drift, automated amplification data analysis can easily provide both false negative and false positive results.
  • one aspect of the present invention is directed to a method of determining the presence of a nucleic acid in a sample by using a fluorescent entity capable of detecting the nucleic acid and amplifying the nucleic acid in the presence of the fluorescent entity.
  • a baseline fluorescence region is determined by analyzing the fluorescence measurements of a number of amplification cycles, and the fluorescence measurements during specific amplification cycles are compared to the baseline fluorescence region to determine the presence or absence of the nucleic acid.
  • the baseline fluorescence region is determined by calculating the slope of a fluorescence intensity verses amplification cycle plot at each of the amplification cycles and choosing the fluorescence measurement of the amplification cycle with the slope having an absolute value closest to zero.
  • the baseline fluorescence region is generated as the square root of the mean square error.
  • the baseline fluorescence region is determined and the fluorescence values are compared thereto after each amplification cycle.
  • the presence of the nucleic acid sequence can be determined more quickly in samples containing higher copy number.
  • the remaining cycles may be used to acquire other information concerning the nucleic acid sample, such as initial copy number and allelic data.
  • the process is automated so that a user can prepare a biological sample and simply place it in a thermal cycler having a sensor for reporting fluorescence values as a function of cycle number and a processor programmed with an algorithm capable of processing the values and reporting a positive or negative result.
  • FIG. 1 is a comparison of three fluorescence monitoring schemes, (A) dsDNA dye, (B) exonuclease probe, and (C) hybridization probe, for PCR amplification, wherein each scheme is illustrated (I) before amplification and (II) after amplification, and fluorescence values are shown (III) once during each cycle of PCR and (IV) continuously during PCR.
  • FIG. 2 is a graph illustrating logistic growth.
  • FIG. 3 shows a comparison of various cycle-verses-fluorescence curve types.
  • FIG. 4 illustrates a sliding window analysis for determining the slope of the fluorescence-verses-cycle number graph at each cycle.
  • FIG. 5 shows typical fluorescence verses amplification cycle graphs for (A) a negative sample and (B) a positive sample.
  • FIG. 6 also shows typical amplification graphs wherein (A) shows fluorescence verses amplification cycle, (B) is the first derivative of fluorescence verses amplification cycle, and (C) is the second derivative of fluorescence verses amplification cycle.
  • FIGS. 7 - 11 show the results for various samples wherein open white circles represent the fluorescence measurement at each cycle, open black circles represent the first derivatives, closed black circles represent second derivatives, large black circles connected by lines represent the points contributing to the baseline calculation, and the horizontal lines illustrate the baseline region.
  • FIGS. 7 and 8 illustrate positive results
  • FIGS. 9 - 11 illustrate negative results.
  • nucleic acid As used herein, “nucleic acid,” “DNA,” and similar terms also include nucleic acid analogs, i.e. analogs having other than a phosphodiester backbone.
  • nucleic acid analogs i.e. analogs having other than a phosphodiester backbone.
  • peptide nucleic acids which are known in the art and have peptide bonds instead of phosphodiester bonds in the backbone, are considered within the scope of the present invention.
  • fluorescence resonance energy transfer pair refers to a pair of fluorophores comprising a donor fluorophore and acceptor fluorophore, wherein the donor fluorophore is capable of transferring resonance energy to the acceptor fluorophore.
  • the emission spectrum of the donor fluorophore overlaps the absorption spectrum of the acceptor fluorophore.
  • the absorption spectrum of the donor fluorophore does not substantially overlap the absorption spectrum of the acceptor fluorophore.
  • FRET oligonucleotide pair refers to a pair of oligonucleotides, each labeled with a member of a fluorescent resonance energy transfer pair, wherein hybridization to complementary target nucleic acid sequences brings the fluorescent entities into a fluorescence resonance energy transfer relationship.
  • the present invention is directed to a method of analyzing a sample for the presence of a nucleic acid wherein the sample is amplified, preferably using PCR, in the presence of a fluorescent probe capable of detecting the presence of the nucleic acid sample.
  • a baseline region is determined by comparing the fluorescence at various amplification cycles, and the fluorescence at each of various amplification cycles is compared to the baseline region to determine whether the fluorescence measurements fall outside of that baseline region.
  • dsDNA double stranded DNA
  • SYBRTM Green I SYBRTM Green I
  • product specificity can be increased by analysis of melting curves or by acquiring fluorescence at a high temperature where nonspecific products have melted.
  • Ririe K M, Rasmussen R P and C T Wittwer Product differentiation by analysis of DNA melting curves during the polymerase chain reaction, Anal. Biochem. 245-154-160, 1997; Morrison T B, J&J Weis and C T Wittwer, Quantification of low copy transcripts by continuous SYBR Green I monitoring during amplification,BioTechniques 24:954-962, 1998.
  • Oligonucleotide probes can also be covalently labeled with fluorescent molecules.
  • Hairpin primers (SunriseTM primers), hairpin probes (Molecular BeaconsTM) and exonuclease probes (TaqManTM) are dual-labeled oligonucleotides that can be monitored during PCR. These probes depend on fluorescence quenching of a fluorophore by a quencher on the same oligonucleotide. Fluorescence increases when hybridization or exonuclease hydrolysis occurs.
  • a preferred design employs two oligonucleotides, each labeled with a fluorescent probe. Hybridization of these oligonucleotides to a target nucleic acid brings the two fluorescent probes close together to allow resonance energy transfer to occur. Wittwer C T, M G Herrmann, A A Moss and R P Rasmussen, Continuous fluorescence monitoring of rapid cycle DNA amplification,BioTechniques 22:130-138, 1997. These hybridization probes require only a single fluorescent label per probe and are easier to design and synthesize than dual labeled probes.
  • Acceptable fluorophore pairs for use as fluorescent resonance energy transfer pairs are well known to those skilled in the art and include, but are not limited to, fluorescein/rhodamine, phycoerythrin/Cy7, fluorescein/Cy5, fluorescein/Cy5.5, fluorescein/LC Red 640, and fluorescein/LC Red 705.
  • SYBRTM Green I, exonuclease probe and hybridization probe designs are shown in FIG. 1.
  • schematics both before (I) and after (II) amplification are shown, as well as cycle verses fluorescence amplification plots of positive and negative controls (III), and temperature verses fluorescence plots from continuous monitoring (IV).
  • SYBR Green I fluorescence increases as more dsDNA is made (FIG. 1A). Because the dye is not sequence specific, a negative control also increases in fluorescence during later cycles as primer dimers are formed. In FIG.
  • FIG. 1B shows dual-labeled fluorescein/rhodarnine probes are cleaved during polymerase extension by 5′-exonuclease activity, separating the fluorophores and increasing the fluorescein emission. The signal generated is cumulative and the fluorescence continues to increase even after the amount of product has reached a plateau.
  • FIG. 1C shows use of a FRET oligonucleotide pair wherein two probes hybridize next to each other, one labeled 3′ with fluorescein and the other labeled 5′ with Cy5. As product accumulates during PCR, fluorescence energy transfer to Cy5 increases. The fluorescence of hybridization probes decreases at high cycle number because of probe/product competition.
  • Standard instruments for PCR complete 30 cycles in about two to four hours.
  • a preferred system is a rapid thermal cycling device using capillary tubes and hot air temperature control. See, for example, U.S. Pat. No. 5,455,175, herein incorporated by reference. Because of the low heat capacity of air and the thin walls and high surface area of capillary tubes, small volume samples could be cycled quickly. The total amplification time for 30 cycles is reduced to 15 minutes with excellent results.
  • capillaries with forced air heating allows precise control of sample temperature at a speed not possible with other designs.
  • sample temperature verses time plots in capillaries show sharp spikes at denaturation and annealing temperatures, whereas several seconds are required for all of the sample to reach equilibrium in conical plastic tubes.
  • Rapid cycling for cycle sequencing reduces sequencing artifacts and minimizes “shadow banding” in dinucleotide repeat amplifications.
  • Swerdlow H, K Dew-Jager and R F Gesteland Rapid cycle sequencing in an air thermal cycler, BioTechniques 15:512-519, 1993; Odelberg S J and R White, A method for accurate amplification of polymorphic CA-repeat sequences, PCR Meth. Appl. 3:7-12, 1993.
  • yield is improved when the sample is exposed as little as possible to high denaturation temperatures.
  • Gustafson C E, R A Alm and T J Trust Effect of heat denaturation of target DNA on the PCR amplification.
  • the RapidCyclerTM developed by Idaho Technology, is an example of a rapid thermal cycling device.
  • the LightCyclerTM is a rapid temperature cycler with a fluorimeter, wherein light emitting diodes are used for excitation and photodiodes are used for detection.
  • the present invention is directed to methods for automating detection nucleic acids with real time PCR. While these algorithms may be applied to any amplification system, it is preferred to integrate these algorithms into the LightCyclerTM platform. These analysis routines are triggered by the completion of rapid thermal cycling for “hands off” amplification, analysis, and final results presentation in a total of less than 15 min. The analysis routines take from ⁇ 1 second for detection and quantification to ⁇ 10 seconds for genotyping. Lab View (National Instruments, Austin, Tex.), a graphical programming language, is preferred for LightCyclerTM instrument control. The LightCyclerTM is a PC-based instrument.
  • FIG. 3 displays various types of amplification curves, all of which have been observed in LightCyclerTM runs.
  • FIGS. 3A and B show curves from samples that are negative with no template present.
  • the fluorescence scales in FIGS. 3A and B are magnified (compared to FIGS. 3 C-F) to demonstrate the baseline drift and to provide algorithms capable of being independent of the fluorescence intensity.
  • This baseline drift of negative reactions must be distinguished from positive reactions of either low copy numbers (FIG. 3C) or high copy numbers (FIG. 3D) of starting template.
  • the method needs to work with various probe designs, including exonuclease (FIG. 3E) and hybridization (FIG. 3F) probes.
  • the baseline is determined as a function of measured fluorescence at a fixed range of cycles near the beginning of amplification.
  • selection of a fixed range of cycles is not adequate because both downward drift (FIG. 3A) and high copy (FIG. 3D) amplifications may be incorrectly called.
  • the background is identified by analyzing the fluorescent measurements over a wide range of amplification cycles.
  • the background is identified by selecting the sliding window (FIG. 4) with the shallowest slope. That is, calculate the slope at each cycle by linear regression of the local neighborhood (for example, a 7 point sliding window).
  • the window with the slope of lowest absolute value (least difference from zero) defines the background region.
  • the variation of these background points about their regression line (the square root of the mean square error) is multiplied by a constant to determine a confidence band. This confidence band will have a slope near zero and is extrapolated across all cycles. If the fluorescence of the last cycle is within the confidence band it is negative, if it is outside the band it is positive.
  • FIG. 5 demonstrates both cases.
  • the maximum slope (first derivative) is obtained from the sliding window analysis already performed for background identification.
  • the second derivatives are calculated by a 3-point sliding window linear regression of the first derivatives. If the curve shape is well behaved (that is, if looking at a graph of FIG. 6, and reading from lowest to highest cycle number, the features occur in the order listed above), then the background is only selected from sliding windows centered at cycle numbers less than the second derivative maximum. This solves the potential analysis problem with FIG. 3D.
  • cycle numbers less than the first derivative maximum or cycle numbers less than the second derivative minimum may be used. It will be further understood that any cycle number between the second derivative maximum and the second derivative minimum is a suitable cutoff cycle for use with this technique and is within the scope of this invention.
  • Another method is to compare the cycle with the greatest fluorescence (which is not necessarily the last cycle) to the confidence band. This is especially suited for hybridization probes that may decrease in fluorescence with extensive cycling, such as seen in FIG. 3F.
  • the cycle with the greatest fluorescence only should be used when the curve shape is well behaved, in order to prevent false positive calls with downward drifts, such as shown in FIG. 3A.
  • the variables to optimize for automatic detection are: 1) the window size for the first derivative estimate, 2) the window size for the second derivative estimate, and 3) the confidence band factor.
  • a reasonable value for the first derivative window size is 7, although 3, 5, 9, and 11 are also quite useful.
  • the preferred window size is 3, but 5, and 7 have also proven to be useful values.
  • a preferred confidence band factor is 20. As the first derivative window size increases the variance estimate is more accurate, but the edge cycles (beginning and ending) are lost.
  • A the number of fluorescence values used to determine the first derivatives. It is convenient to use odd numbers, so that the first derivatives correspond to integer cycle numbers. As discussed above, reasonable values include 3, 5, 7, 9, and 11. Preferably, 7 is used as the first derivative window size.
  • B the number of first derivative values used to determine the second derivatives. Again, it is convenient to use odd numbers, so that the second derivative values also correspond to integer cycle numbers. Reasonable values include 3, 5, and 7, with 3 being the preferred value.
  • C the confidence band factor. This factor determines the confidence band by multiplying it by a variance measure, preferably the square root of the mean square error.
  • the first step is to calculate the first and second derivatives.
  • a preferred method is to determine the first derivatives as the slope of a linear regression line through A points, and assigning the value to the central cycle number. Some cycles on either edge cannot be assigned first derivatives, but first derivatives can be provided for cycles (A+1)/2 through N-(A ⁇ 1)/2.
  • the second derivatives are calculated as the slope of the first derivative points and assigned to cycles (A+1)/2+(B ⁇ 1)/2 through [N ⁇ (A ⁇ 1)/2] ⁇ (B ⁇ 1)/2. Calculation of the first and second derivatives provide arrays Y′i and Y′′i, with some edge values missing. In FIG. 7, the first and second derivatives are displayed as open black circles and closed black circles, respectively.
  • the next step is to determine whether the fluorescence curve has a well-behaved shape.
  • the well-behaved shape occurs when the cycles with minimum fluorescence, maximum second derivative, maximum first derivative, minimum second derivative, and maximum fluorescence occur in that order, from low to high cycle number.
  • the baseline is then determined. If the fluorescence curve does not have the expected shape, the cycle whose first derivative is closest to zero is used. If the fluorescence curve has a well-behaved shape, the cycle whose first derivative is closest to zero chosen from among all cycles prior to the cycle with the maximum second derivative (again, any cycle between the maximum second derivative and the minimum second derivative may also be used as the cutoff cycle number).
  • the baseline is drawn through the fluorescence value of the chosen cycle with a slope of its first derivative. In FIG. 7, the A points contributing to the first derivative calculation for the baseline are displayed as large black dots connected by a line.
  • the next step is to determine the test point cycle, that is, the cycle used to compare against the baseline for determining a positive or negative result. It the curve is not well-behaved, the test point is the last cycle. If the fluorescence curve is well-behaved, the test point is the cycle with fluorescence farthest from the baseline. The test point fluorescence of a negative sample can be predicted as the intersection of the baseline with the test point cycle.
  • a confidence interval can be determined about the predicted negative test point. Preferably, this is done by finding the square root of the mean square error about the baseline of A points used to determine the baseline. This is multiplied by C. The product is added to the predicted negative test point to get the upper fluorescence limit of the confidence interval and is subtracted from the predicted negative test point to get the lower limit of the confidence interval. These limits are shown on FIG. 7 as two solid horizontal lines.
  • the final step is to declare the sample positive or negative. If the test point fluorescence is outside of the confidence interval, the sample is positive. If it is within the interval, the sample is negative.
  • FIGS. 7 and 8 are samples which are positive, while FIGS. 9 - 11 are negative samples.

Landscapes

  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Immunology (AREA)
  • Organic Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Analytical Chemistry (AREA)
  • Molecular Biology (AREA)
  • General Physics & Mathematics (AREA)
  • Biotechnology (AREA)
  • Microbiology (AREA)
  • Biophysics (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • Genetics & Genomics (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
  • Automatic Analysis And Handling Materials Therefor (AREA)

Abstract

A method is described for analyzing a sample for the presence of a nucleic acid wherein the sample is amplified, illustratively using PCR, in the presence of a fluorescent probe capable of detecting the presence of the nucleic acid sample. A baseline region is determined by comparing the fluorescence at various amplification cycles, and the fluorescence at a selected amplification cycle is compared to the baseline region to determine whether the fluorescence measurement falls outside of that baseline region, indicating the presence of the nucleic acid.

Description

    CROSS-REFERENCE TO RELATED APPLICATION
  • This application claims priority under 35 U.S.C. §119 (e) to U.S. Provisional Application No. 60/131,256, filed Apr. 27, 1999, which is expressly incorporated by reference herein.[0001]
  • FIELD OF THE INVENTION
  • The present invention relates to a method of analyzing a sample for the presence of a nucleic acid. More particularly, the present invention is directed to an automated method for detecting and reporting the presence of a predetermined nucleic acid in a sample using polymerase chain reaction and a fluorescent detecting entity. [0002]
  • BACKGROUND AND SUMMARY OF THE INVENTION
  • Amplification of DNA by polymerase chain reaction (PCR) is a technique fundamental to molecular biology. Nucleic acid analysis by PCR requires sample preparation, amplification, and product analysis. Although these steps are usually performed sequentially, amplification and analysis can occur simultaneously. DNA dyes or fluorescent probes can be added to the PCR mixture before amplification and used to analyze PCR products during amplification. Sample analysis occurs concurrently with amplification in the same tube within the same instrument. This combined approach decreases sample handling, saves time, and greatly reduces the risk of product contamination for subsequent reactions, as there is no need to remove the samples from their closed containers for further analysis. The concept of combining amplification with product analysis has become known as “real time” PCR. See, for example, WO/9746707A2, WO/9746712A2, WO/9746714A1, all published Dec. 11, 1997, incorporated herein by reference. [0003]
  • Monitoring fluorescence each cycle of PCR initially involved the use of ethidium bromide. Higuchi R, G Dollinger, P S Walsh and R. Griffith, Simultaneous amplification and detection of specific DNA sequences, Bio/Technology 10:413-417, 1992; Higuchi R, C Fockler G Dollinger and R Watson, Kinetic PCR analysis: real time monitoring of DNA amplification reactions, Bio/Technology 11:1026-1030, 1993. In that system fluorescence is measured once per cycle as a relative measure of product concentration. Ethidium bromide detects double stranded DNA; if template is present fluorescence intensity increases with temperature cycling. Furthermore, the cycle number where an increase in fluorescence is first detected increases inversely proportionally to the log of the initial template concentration. Other fluorescent systems have been developed that are capable of providing additional data concerning the nucleic acid concentration and sequence. [0004]
  • While PCR is an invaluable molecular biology tool, the practical implementation of real time PCR techniques has lagged behind the conceptual promise. Currently available instrumentation does not actually analyze data during PCR; it simply acquires the data for later analysis. After PCR has been completed, multiple manual steps are necessary to analyze the acquired data, and human judgment is typically required to provide the analysis result. What is needed is a system for automating data acquisition and analysis so that no user intervention is required for reporting the analytical results. Thus, when the temperature cycling in a polymerase chain reaction amplification is complete, the system software is automatically triggered and the results, for example, the presence or absence of a given pathogen, is immediately displayed on screen. Algorithms for detection, quantification, and genotyping are needed. Moreover, initiation of the analysis algorithm can be implemented prior to completion of temperature cycling. Data processing can occur during amplification and concomitant analysis results can be used to modify temperature cycling and to acquire additional data during the latter stages of the amplification procedure to optimize amplification protocol and data quality. [0005]
  • A major problem in automating PCR data analysis is identification of baseline fluorescence. Background fluorescence varies from reaction to reaction. Moreover, baseline drift, wherein fluorescence increases or decreases without relation to amplification of nucleic acids in the sample, is a common occurrence. Prior attempts to automate amplification data analysis involved setting the baseline fluorescence as that measured at one or more predetermined early cycle numbers. This technique accounts for the variation in background fluorescence, but it does not compensate for baseline drift. Without compensation for baseline drift, automated amplification data analysis can easily provide both false negative and false positive results. [0006]
  • Thus, one aspect of the present invention is directed to a method of determining the presence of a nucleic acid in a sample by using a fluorescent entity capable of detecting the nucleic acid and amplifying the nucleic acid in the presence of the fluorescent entity. A baseline fluorescence region is determined by analyzing the fluorescence measurements of a number of amplification cycles, and the fluorescence measurements during specific amplification cycles are compared to the baseline fluorescence region to determine the presence or absence of the nucleic acid. In a preferred embodiment, the baseline fluorescence region is determined by calculating the slope of a fluorescence intensity verses amplification cycle plot at each of the amplification cycles and choosing the fluorescence measurement of the amplification cycle with the slope having an absolute value closest to zero. Preferably, the baseline fluorescence region is generated as the square root of the mean square error. [0007]
  • In another embodiment, the baseline fluorescence region is determined and the fluorescence values are compared thereto after each amplification cycle. Thus, the presence of the nucleic acid sequence can be determined more quickly in samples containing higher copy number. Furthermore, the remaining cycles may be used to acquire other information concerning the nucleic acid sample, such as initial copy number and allelic data. [0008]
  • In an additional embodiment, the process is automated so that a user can prepare a biological sample and simply place it in a thermal cycler having a sensor for reporting fluorescence values as a function of cycle number and a processor programmed with an algorithm capable of processing the values and reporting a positive or negative result.[0009]
  • Additional features of the present invention will become apparent to those skilled in the art upon consideration of the following detailed description of preferred embodiments exemplifying the best mode of carrying out the invention as presently perceived. [0010]
  • BRIEF DESCRIPTION OF THE DRAWINGS
  • FIG. 1 is a comparison of three fluorescence monitoring schemes, (A) dsDNA dye, (B) exonuclease probe, and (C) hybridization probe, for PCR amplification, wherein each scheme is illustrated (I) before amplification and (II) after amplification, and fluorescence values are shown (III) once during each cycle of PCR and (IV) continuously during PCR. [0011]
  • FIG. 2 is a graph illustrating logistic growth. [0012]
  • FIG. 3 shows a comparison of various cycle-verses-fluorescence curve types. [0013]
  • FIG. 4 illustrates a sliding window analysis for determining the slope of the fluorescence-verses-cycle number graph at each cycle. [0014]
  • FIG. 5 shows typical fluorescence verses amplification cycle graphs for (A) a negative sample and (B) a positive sample. [0015]
  • FIG. 6 also shows typical amplification graphs wherein (A) shows fluorescence verses amplification cycle, (B) is the first derivative of fluorescence verses amplification cycle, and (C) is the second derivative of fluorescence verses amplification cycle. [0016]
  • FIGS. [0017] 7-11 show the results for various samples wherein open white circles represent the fluorescence measurement at each cycle, open black circles represent the first derivatives, closed black circles represent second derivatives, large black circles connected by lines represent the points contributing to the baseline calculation, and the horizontal lines illustrate the baseline region. FIGS. 7 and 8 illustrate positive results, while FIGS. 9-11 illustrate negative results.
  • DETAILED DESCRIPTION OF THE INVENTION
  • In describing and claiming the invention, the following terminology will be used in accordance with the definitions set forth below. [0018]
  • As used herein, “nucleic acid,” “DNA,” and similar terms also include nucleic acid analogs, i.e. analogs having other than a phosphodiester backbone. For example, the so-called “peptide nucleic acids,” which are known in the art and have peptide bonds instead of phosphodiester bonds in the backbone, are considered within the scope of the present invention. [0019]
  • As used herein, “fluorescence resonance energy transfer pair” or “FRET pair” refers to a pair of fluorophores comprising a donor fluorophore and acceptor fluorophore, wherein the donor fluorophore is capable of transferring resonance energy to the acceptor fluorophore. In other words the emission spectrum of the donor fluorophore overlaps the absorption spectrum of the acceptor fluorophore. In preferred fluorescence resonance energy transfer pairs, the absorption spectrum of the donor fluorophore does not substantially overlap the absorption spectrum of the acceptor fluorophore. [0020]
  • As used herein, “FRET oligonucleotide pair” refers to a pair of oligonucleotides, each labeled with a member of a fluorescent resonance energy transfer pair, wherein hybridization to complementary target nucleic acid sequences brings the fluorescent entities into a fluorescence resonance energy transfer relationship. [0021]
  • The present invention is directed to a method of analyzing a sample for the presence of a nucleic acid wherein the sample is amplified, preferably using PCR, in the presence of a fluorescent probe capable of detecting the presence of the nucleic acid sample. A baseline region is determined by comparing the fluorescence at various amplification cycles, and the fluorescence at each of various amplification cycles is compared to the baseline region to determine whether the fluorescence measurements fall outside of that baseline region. [0022]
  • Many different probes have recently become available for monitoring PCR. Although not sequence specific, double stranded DNA (dsDNA) specific dyes can be used in any amplification without the need for probe synthesis. Such dyes include ethidium bromide and SYBR™ Green I. With dsDNA dyes, product specificity can be increased by analysis of melting curves or by acquiring fluorescence at a high temperature where nonspecific products have melted. Ririe K M, Rasmussen R P and C T Wittwer, Product differentiation by analysis of DNA melting curves during the polymerase chain reaction, Anal. Biochem. 245-154-160, 1997; Morrison T B, J&J Weis and C T Wittwer, Quantification of low copy transcripts by continuous SYBR Green I monitoring during amplification,BioTechniques 24:954-962, 1998. [0023]
  • Oligonucleotide probes can also be covalently labeled with fluorescent molecules. Hairpin primers (Sunrise™ primers), hairpin probes (Molecular Beacons™) and exonuclease probes (TaqMan™) are dual-labeled oligonucleotides that can be monitored during PCR. These probes depend on fluorescence quenching of a fluorophore by a quencher on the same oligonucleotide. Fluorescence increases when hybridization or exonuclease hydrolysis occurs. [0024]
  • A preferred design employs two oligonucleotides, each labeled with a fluorescent probe. Hybridization of these oligonucleotides to a target nucleic acid brings the two fluorescent probes close together to allow resonance energy transfer to occur. Wittwer C T, M G Herrmann, A A Moss and R P Rasmussen, Continuous fluorescence monitoring of rapid cycle DNA amplification,BioTechniques 22:130-138, 1997. These hybridization probes require only a single fluorescent label per probe and are easier to design and synthesize than dual labeled probes. Acceptable fluorophore pairs for use as fluorescent resonance energy transfer pairs are well known to those skilled in the art and include, but are not limited to, fluorescein/rhodamine, phycoerythrin/Cy7, fluorescein/Cy5, fluorescein/Cy5.5, fluorescein/LC Red 640, and fluorescein/LC Red 705. [0025]
  • SYBR™ Green I, exonuclease probe and hybridization probe designs are shown in FIG. 1. For each design, schematics both before (I) and after (II) amplification are shown, as well as cycle verses fluorescence amplification plots of positive and negative controls (III), and temperature verses fluorescence plots from continuous monitoring (IV). SYBR Green I fluorescence increases as more dsDNA is made (FIG. 1A). Because the dye is not sequence specific, a negative control also increases in fluorescence during later cycles as primer dimers are formed. In FIG. 1B, dual-labeled fluorescein/rhodarnine probes are cleaved during polymerase extension by 5′-exonuclease activity, separating the fluorophores and increasing the fluorescein emission. The signal generated is cumulative and the fluorescence continues to increase even after the amount of product has reached a plateau. FIG. 1C shows use of a FRET oligonucleotide pair wherein two probes hybridize next to each other, one labeled 3′ with fluorescein and the other labeled 5′ with Cy5. As product accumulates during PCR, fluorescence energy transfer to Cy5 increases. The fluorescence of hybridization probes decreases at high cycle number because of probe/product competition. [0026]
  • Standard instruments for PCR complete 30 cycles in about two to four hours. A preferred system is a rapid thermal cycling device using capillary tubes and hot air temperature control. See, for example, U.S. Pat. No. 5,455,175, herein incorporated by reference. Because of the low heat capacity of air and the thin walls and high surface area of capillary tubes, small volume samples could be cycled quickly. The total amplification time for 30 cycles is reduced to 15 minutes with excellent results. [0027]
  • The use of capillaries with forced air heating allows precise control of sample temperature at a speed not possible with other designs. For example, sample temperature verses time plots in capillaries show sharp spikes at denaturation and annealing temperatures, whereas several seconds are required for all of the sample to reach equilibrium in conical plastic tubes. Wittwer, C T, G B Reed and K M Ririe, Rapid cycle DNA amplification, in K Mullis, F Ferre, and R Gibbs (Eds.), The polymerase chain reaction, Springer-Verlag, Deerfield Beach, Fla. pp. 174-181, 1994; Wittwer, C T, B C Marshall, G B Reed, and J L Cherry, Rapid cycle allele-specific amplification: studies with the cystic fibrosis delta F508 locus, Clin. Chem., 39:804-809, 1993. Rapid temperature cycling with minimal annealing and denaturation times improves quantitative PCR and increases the discrimination of allele specific amplification. Weis, J H, S S Tan, B K Martin, and C T Wittwer, Detection of rare mRNA species via quantitative RT-PCR, Trends in Genetics, 8:263-4, 1992; Tan S T and J H Weis, Development of a sensitive reverse transcriptase PCR assay, RT-RPCR, utilizing rapid cycle times, PCR Meth. and Appl. 2:137-143, 1992. Rapid cycling for cycle sequencing reduces sequencing artifacts and minimizes “shadow banding” in dinucleotide repeat amplifications. Swerdlow H, K Dew-Jager and R F Gesteland, Rapid cycle sequencing in an air thermal cycler, BioTechniques 15:512-519, 1993; Odelberg S J and R White, A method for accurate amplification of polymorphic CA-repeat sequences, PCR Meth. Appl. 3:7-12, 1993. For long PCR, yield is improved when the sample is exposed as little as possible to high denaturation temperatures. Gustafson C E, R A Alm and T J Trust, Effect of heat denaturation of target DNA on the PCR amplification. Gene 23:241-244, 1993. The RapidCycler™, developed by Idaho Technology, is an example of a rapid thermal cycling device. The LightCycler™ is a rapid temperature cycler with a fluorimeter, wherein light emitting diodes are used for excitation and photodiodes are used for detection. [0028]
  • The present invention is directed to methods for automating detection nucleic acids with real time PCR. While these algorithms may be applied to any amplification system, it is preferred to integrate these algorithms into the LightCycler™ platform. These analysis routines are triggered by the completion of rapid thermal cycling for “hands off” amplification, analysis, and final results presentation in a total of less than 15 min. The analysis routines take from <1 second for detection and quantification to <10 seconds for genotyping. Lab View (National Instruments, Austin, Tex.), a graphical programming language, is preferred for LightCycler™ instrument control. The LightCycler™ is a PC-based instrument. [0029]
  • Perhaps the most basic analysis of real time PCR data is a judgement of whether a targeted nucleic acid is present. If the nucleic acid is present, further quantification and genotyping may take place. In many cases, a yes/no judgement is all that is needed. For example, one may want to determine whether [0030] E. coli 0157:H7 is in a sample of hamburger, whether anthrax is on a swab from a soldier; or whether hepatitis C is in a unit of blood. Real time PCR can improve yes/no detection over end point PCR assays because fluorescence is acquired at each cycle.
  • Inspection of cycle verses fluorescence data from positive and negative real time PCR runs (see FIG. 1BIII and [0031] 1CIII) suggests that discrimination is simple. The positive samples increase with cycle number while the negative samples remain at baseline. A trained observer expects positive samples to follow an S-shape curve, beginning with a baseline, followed by an exponential segment, and finishing with a plateau. The expected curve is similar to the logistic model for population growth, where the rate of growth is proportional to both the population size y and to the difference L-y, where L is the maximum population that can be supported. For small y, growth is exponential, but as y nears L the growth rate approaches zero. An example of logistic growth is shown in FIG. 2.
  • Although intuitively simple, accurately discriminating between positive and negative samples is not easy in practice. The simplest approach is to set a horizontal fluorescence threshold as a discriminator between positive and negative samples. This works best with a stable baseline (between and within samples) and a known fluorescence intensity that correlates with “positive.” Although this method will work on obvious samples (e.g. FIG. 1BIII and [0032] 1CIII), a more robust algorithm is desired that will work under a wider variety of conditions. For example, the baseline may drift and the fluorescence intensity may vary greatly between different samples and probe techniques. Thus, the present invention is directed to a method that will: (1) automatically identify the baseline, (2) use the baseline variance to establish a confidence region, and (3) call each sample positive or negative based on the relationship of the confidence region to the fluorescence data.
  • FIG. 3 displays various types of amplification curves, all of which have been observed in LightCycler™ runs. FIGS. 3A and B show curves from samples that are negative with no template present. The fluorescence scales in FIGS. 3A and B are magnified (compared to FIGS. [0033] 3C-F) to demonstrate the baseline drift and to provide algorithms capable of being independent of the fluorescence intensity. There is always some baseline drift during cycling. This drift usually is greatest at the beginning of cycling but later levels off, and may be either downward (FIG. 3A) or upward (FIG. 3B). This baseline drift of negative reactions must be distinguished from positive reactions of either low copy numbers (FIG. 3C) or high copy numbers (FIG. 3D) of starting template. The method needs to work with various probe designs, including exonuclease (FIG. 3E) and hybridization (FIG. 3F) probes.
  • Automatic identification of the background is surprisingly difficult. In prior art methods, the baseline is determined as a function of measured fluorescence at a fixed range of cycles near the beginning of amplification. However, selection of a fixed range of cycles is not adequate because both downward drift (FIG. 3A) and high copy (FIG. 3D) amplifications may be incorrectly called. [0034]
  • In the present invention, the background is identified by analyzing the fluorescent measurements over a wide range of amplification cycles. Preferably, the background is identified by selecting the sliding window (FIG. 4) with the shallowest slope. That is, calculate the slope at each cycle by linear regression of the local neighborhood (for example, a 7 point sliding window). The window with the slope of lowest absolute value (least difference from zero) defines the background region. Once the background region has been identified, the variation of these background points about their regression line (the square root of the mean square error) is multiplied by a constant to determine a confidence band. This confidence band will have a slope near zero and is extrapolated across all cycles. If the fluorescence of the last cycle is within the confidence band it is negative, if it is outside the band it is positive. FIG. 5 demonstrates both cases. [0035]
  • This algorithm should work well in most cases. However, with the high copy fluorescence curve type (FIG. 3D), the shallowest slope might be found at early cycles (resulting in a correct positive call) or at late cycles (resulting in an incorrect negative call). This exception may be handled by analyzing the curve shape. In a well-behaved amplification, the expected amplification curve shape is ordered by cycle number as follows: [0036]
  • 1. Minimum fluorescence [0037]
  • 2. Maximum second derivative (F″) [0038]
  • 3. Maximum first derivative (F′) [0039]
  • 4. Minimum second derivative (F′) [0040]
  • 5. Maximum fluorescence [0041]
  • This gives the characteristic S-curve shape we expect during PCR (FIG. 6A). The maximum slope (first derivative) is obtained from the sliding window analysis already performed for background identification. Preferably, the second derivatives are calculated by a 3-point sliding window linear regression of the first derivatives. If the curve shape is well behaved (that is, if looking at a graph of FIG. 6, and reading from lowest to highest cycle number, the features occur in the order listed above), then the background is only selected from sliding windows centered at cycle numbers less than the second derivative maximum. This solves the potential analysis problem with FIG. 3D. In other preferred embodiments, cycle numbers less than the first derivative maximum or cycle numbers less than the second derivative minimum may be used. It will be further understood that any cycle number between the second derivative maximum and the second derivative minimum is a suitable cutoff cycle for use with this technique and is within the scope of this invention. [0042]
  • Another method is to compare the cycle with the greatest fluorescence (which is not necessarily the last cycle) to the confidence band. This is especially suited for hybridization probes that may decrease in fluorescence with extensive cycling, such as seen in FIG. 3F. The cycle with the greatest fluorescence only should be used when the curve shape is well behaved, in order to prevent false positive calls with downward drifts, such as shown in FIG. 3A. [0043]
  • The variables to optimize for automatic detection are: 1) the window size for the first derivative estimate, 2) the window size for the second derivative estimate, and 3) the confidence band factor. A reasonable value for the first derivative window size is 7, although 3, 5, 9, and 11 are also quite useful. For the second derivative the preferred window size is 3, but 5, and 7 have also proven to be useful values. A preferred confidence band factor is 20. As the first derivative window size increases the variance estimate is more accurate, but the edge cycles (beginning and ending) are lost. [0044]
  • This algorithm is best understood by referring to the fluorescence verses cycle test result plot shown in FIGS. [0045] 7-11. The input data consist of one fluorescence value for each cycle of amplification, shown as the closed white circles. Let this equal array Yi, where i is the cycle number and N is the total number of cycles. The detection criteria are:
  • A=the number of fluorescence values used to determine the first derivatives. It is convenient to use odd numbers, so that the first derivatives correspond to integer cycle numbers. As discussed above, reasonable values include 3, 5, 7, 9, and 11. Preferably, 7 is used as the first derivative window size. [0046]
  • B=the number of first derivative values used to determine the second derivatives. Again, it is convenient to use odd numbers, so that the second derivative values also correspond to integer cycle numbers. Reasonable values include 3, 5, and 7, with 3 being the preferred value. [0047]
  • C=the confidence band factor. This factor determines the confidence band by multiplying it by a variance measure, preferably the square root of the mean square error. [0048]
  • The first step is to calculate the first and second derivatives. Although there are many ways to accomplish this, a preferred method is to determine the first derivatives as the slope of a linear regression line through A points, and assigning the value to the central cycle number. Some cycles on either edge cannot be assigned first derivatives, but first derivatives can be provided for cycles (A+1)/2 through N-(A−1)/2. Similarly, the second derivatives are calculated as the slope of the first derivative points and assigned to cycles (A+1)/2+(B−1)/2 through [N−(A−1)/2]−(B−1)/2. Calculation of the first and second derivatives provide arrays Y′i and Y″i, with some edge values missing. In FIG. 7, the first and second derivatives are displayed as open black circles and closed black circles, respectively. [0049]
  • The next step is to determine whether the fluorescence curve has a well-behaved shape. As discussed above, the well-behaved shape occurs when the cycles with minimum fluorescence, maximum second derivative, maximum first derivative, minimum second derivative, and maximum fluorescence occur in that order, from low to high cycle number. [0050]
  • The baseline is then determined. If the fluorescence curve does not have the expected shape, the cycle whose first derivative is closest to zero is used. If the fluorescence curve has a well-behaved shape, the cycle whose first derivative is closest to zero chosen from among all cycles prior to the cycle with the maximum second derivative (again, any cycle between the maximum second derivative and the minimum second derivative may also be used as the cutoff cycle number). The baseline is drawn through the fluorescence value of the chosen cycle with a slope of its first derivative. In FIG. 7, the A points contributing to the first derivative calculation for the baseline are displayed as large black dots connected by a line. [0051]
  • The next step is to determine the test point cycle, that is, the cycle used to compare against the baseline for determining a positive or negative result. It the curve is not well-behaved, the test point is the last cycle. If the fluorescence curve is well-behaved, the test point is the cycle with fluorescence farthest from the baseline. The test point fluorescence of a negative sample can be predicted as the intersection of the baseline with the test point cycle. [0052]
  • Next, a confidence interval can be determined about the predicted negative test point. Preferably, this is done by finding the square root of the mean square error about the baseline of A points used to determine the baseline. This is multiplied by C. The product is added to the predicted negative test point to get the upper fluorescence limit of the confidence interval and is subtracted from the predicted negative test point to get the lower limit of the confidence interval. These limits are shown on FIG. 7 as two solid horizontal lines. [0053]
  • The final step is to declare the sample positive or negative. If the test point fluorescence is outside of the confidence interval, the sample is positive. If it is within the interval, the sample is negative. FIGS. 7 and 8 are samples which are positive, while FIGS. [0054] 9-11 are negative samples.
  • Although the invention has been described in detail with reference to preferred embodiments, variations and modifications exist within the scope and spirit of the invention as described and defined in the following claims. [0055]

Claims (17)

1. A method for determining the presence of a nucleic acid comprising the steps of
providing a fluorescent entity capable of indicating the presence of the nucleic acid and capable of providing a signal related to the quantity of the nucleic acid,
amplifying the nucleic acid through a plurality of amplification cycles in the presence of the fluorescent entity,
measuring the fluorescence of the fluorescent entity during each of the plurality of amplification cycles,
analyzing the measured fluorescence to identify amplification cycles for use in establishing a baseline fluorescence region, and
ascertaining whether the fluorescence measurement during any of the plurality of amplification cycles is outside the baseline fluorescence region.
2. The method of claim 1 wherein the cycles used in establishing the baseline fluorescence region are identified by calculating slopes of a fluorescence verses amplification cycle plot at the plurality of amplification cycles.
3. The method of claim 2 wherein the slopes are calculated by linear regression of each local neighborhood.
4. The method of claim 2 wherein the cycles used in establishing the baseline fluorescence region comprise an interval of cycles that includes the amplification cycle with the slope having an absolute value closest to zero
5. The method of claim 4 wherein the amplification cycle with the slope having the absolute value closest to zero comprises a central cycle of the interval of cycles used to establish the baseline fluorescence region.
6. The method of claim 4 wherein the baseline fluorescence is centered on a line passing through the measured fluorescence of the amplification cycle with the slope having the absolute value closest to zero, the line further comprising a slope equal to the slope having the absolute value closest to zero, and the baseline fluorescence region including an area on at least one side of the line defined by a variance measure of the fluorescence values.
7. The method of claim 6 wherein the variance measure is calculated from the amplification cycles used in establishing the baseline fluorescence.
8. The method of claim 1 wherein the fluorescent entity comprises a FRET oligonucleotide pair.
9. The method of claim 1 further comprising the steps of
generating a well-behaved amplification curve,
determining a cutoff cycle,
calculating the slope of the fluorescence verses cycle plot for each of the of amplification cycles prior to the cutoff cycle, and
choosing the fluorescence measurement of the amplification cycle with the slope having an absolute value closest to zero.
10. The method of claim 9 wherein the cutoff cycle is selected from the group consisting of a maximum second derivative, a maximum first derivative, and a minimum second derivative of a fluorescence verses cycle curve.
11. The method of claim 1 further comprising the steps of
generating a well-behaved amplification curve,
determining the cycle having the maximum fluorescence, and
determining whether the maximum fluorescence is outside the baseline fluorescence region.
12. The method of claim 1 further comprising the steps of
generating an amplification curve which is not well behaved, and
determining if the fluorescence of a last cycle tested is outside the baseline fluorescence region.
13. A method for determining the presence of a nucleic acid comprising the steps of
providing a fluorescent entity capable of indicating the presence of the nucleic acid and capable of providing a signal related to the quantity of the nucleic acid,
initially amplifying the nucleic acid through a set number of amplification cycles in the presence of the fluorescent entity,
measuring the fluorescence of the fluorescent entity during each of the amplification cycles,
analyzing the measured fluorescence to identify amplification cycles for use in establishing a baseline fluorescence region,
ascertaining whether a positive result is indicated by determining whether the fluorescence measurement during any of the of amplification cycles is outside the baseline fluorescence region, and
continuing amplification through additional amplification cycles and repeating the measuring, determining, and ascertaining steps after each additional amplification cycle until either the positive result is obtained or a maximum cycle number is reached.
14. The method of claim 13 further comprising the steps of
performing additional amplification cycles after the positive result is obtained, and
analyzing the nucleic acid for additional information.
15. The method of claim 13 wherein the nucleic acid is further analyzed to determine the presence of a particular allele.
16. The method of claim 13 wherein the measured fluorescence is further analyzed to determine initial concentration.
17. An automated method for determining the presence of a nucleic acid comprising the steps of
placing a sample into a container containing a fluorescent entity capable of indicating the presence of the nucleic acid and capable of providing a signal related to the quantity of the nucleic acid,
placing the container into a device for amplifying the nucleic acid through a plurality of amplification cycles in the presence of the fluorescent entity, measuring the fluorescence of the fluorescent entity during each of the plurality of amplification cycles, determining a baseline fluorescence region by analyzing the fluorescence at various amplification cycles, and outputting a positive result if the fluorescence measurement during any of the plurality of amplification cycles is outside the baseline fluorescence region.
US09/391,811 1999-04-27 1999-09-08 Automated analysis of real-time nucleic acid amplification Expired - Lifetime US6387621B1 (en)

Priority Applications (6)

Application Number Priority Date Filing Date Title
US09/391,811 US6387621B1 (en) 1999-04-27 1999-09-08 Automated analysis of real-time nucleic acid amplification
EP00303519A EP1059523B1 (en) 1999-04-27 2000-04-26 Automated real-time analysis of nucleic acid amplification
DE60035459T DE60035459T2 (en) 1999-04-27 2000-04-26 Automated real-time analysis of nucleic acid amplification
AT00303519T ATE366821T1 (en) 1999-04-27 2000-04-26 AUTOMATED REAL-TIME ANALYSIS OF NUCLEIC ACID AMPLIFICATION
JP2000128560A JP4689004B2 (en) 1999-04-27 2000-04-27 Automatic analysis of real-time nucleic acid amplification
US10/074,169 US20020123062A1 (en) 1999-04-27 2002-02-12 Automated analysis of real-time nucleic acid amplification

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US13125699P 1999-04-27 1999-04-27
US09/391,811 US6387621B1 (en) 1999-04-27 1999-09-08 Automated analysis of real-time nucleic acid amplification

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US10/074,169 Continuation-In-Part US20020123062A1 (en) 1999-04-27 2002-02-12 Automated analysis of real-time nucleic acid amplification

Publications (2)

Publication Number Publication Date
US20020042051A1 true US20020042051A1 (en) 2002-04-11
US6387621B1 US6387621B1 (en) 2002-05-14

Family

ID=26829286

Family Applications (2)

Application Number Title Priority Date Filing Date
US09/391,811 Expired - Lifetime US6387621B1 (en) 1999-04-27 1999-09-08 Automated analysis of real-time nucleic acid amplification
US10/074,169 Abandoned US20020123062A1 (en) 1999-04-27 2002-02-12 Automated analysis of real-time nucleic acid amplification

Family Applications After (1)

Application Number Title Priority Date Filing Date
US10/074,169 Abandoned US20020123062A1 (en) 1999-04-27 2002-02-12 Automated analysis of real-time nucleic acid amplification

Country Status (5)

Country Link
US (2) US6387621B1 (en)
EP (1) EP1059523B1 (en)
JP (1) JP4689004B2 (en)
AT (1) ATE366821T1 (en)
DE (1) DE60035459T2 (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080176221A1 (en) * 2002-11-15 2008-07-24 Gen-Probe Incorporated Assay and compositions for detection of bacillus anthracis nucleic acid
US20090035779A1 (en) * 2007-06-29 2009-02-05 Roche Molecular Systems, Inc. Systems and methods for determining cross-talk coefficients in pcr and other data sets
US8306754B2 (en) 2007-12-26 2012-11-06 Arkray, Inc. Nucleic acid amplification determining method and nucleic acid amplification determining device
US20120282611A1 (en) * 2009-10-21 2012-11-08 Brandeis University Methods, kits and reaction mixtures for analyzing single-stranded nucleic acid sequences
US20130273547A1 (en) * 2012-04-16 2013-10-17 Samsung Techwin Co., Ltd. Method to determine and correct baseline and to characterize pcr amplification kinetics
CN111426661A (en) * 2020-04-09 2020-07-17 基蛋生物科技股份有限公司 Fluorescence data acquisition and discrimination method in nucleic acid amplification stage
CN112669910A (en) * 2020-12-30 2021-04-16 杭州博日科技股份有限公司 Amplification curve baseline determination method and device and electronic equipment

Families Citing this family (74)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6387621B1 (en) * 1999-04-27 2002-05-14 University Of Utah Research Foundation Automated analysis of real-time nucleic acid amplification
EP1290219B1 (en) * 2000-03-27 2005-06-01 Applera Corporation Improved invasion assay
WO2002089844A1 (en) * 2001-05-07 2002-11-14 Mayo Foundation For Medical Education And Research Legionella detection with fluorescent resonance energy transfer
US6979541B1 (en) 2001-07-26 2005-12-27 University Of Utah Research Foundation Methods for identifying chromosomal aneuploidy
US7188030B2 (en) * 2001-08-21 2007-03-06 Applera Corporation Automatic threshold setting for quantitative polymerase chain reaction
DE60213730T2 (en) * 2001-08-31 2007-08-16 The University Of Utah Research Foundation, Salt Lake City Real-time quantification with internal standards
US7373253B2 (en) * 2002-02-12 2008-05-13 Idaho Technology Multi-test analysis of real-time nucleic acid amplification
EP1335028B1 (en) * 2002-02-12 2008-01-16 University of Utah Research Foundation Multi-test analysis of real-time nucleic acid amplifications
EP2428588B1 (en) 2002-02-21 2016-08-10 Alere San Diego, Inc. Recombinase Polymerase Amplification
US8030000B2 (en) * 2002-02-21 2011-10-04 Alere San Diego, Inc. Recombinase polymerase amplification
US7399590B2 (en) 2002-02-21 2008-07-15 Asm Scientific, Inc. Recombinase polymerase amplification
CA2477761A1 (en) 2002-03-01 2003-09-12 Ravgen, Inc. Methods for detection of genetic disorders
US6977162B2 (en) * 2002-03-01 2005-12-20 Ravgen, Inc. Rapid analysis of variations in a genome
US7442506B2 (en) * 2002-05-08 2008-10-28 Ravgen, Inc. Methods for detection of genetic disorders
US20070178478A1 (en) * 2002-05-08 2007-08-02 Dhallan Ravinder S Methods for detection of genetic disorders
US7727720B2 (en) * 2002-05-08 2010-06-01 Ravgen, Inc. Methods for detection of genetic disorders
AUPS267802A0 (en) 2002-05-30 2002-06-20 Bio-Molecular Holdings Pty Limited Improved dna amplification apparatus and method
US7083974B2 (en) * 2002-07-12 2006-08-01 Applera Corporation Rotatable sample disk and method of loading a sample disk
EP1552001A4 (en) * 2002-09-19 2006-01-04 Applera Corp Methods and compositions for detecting targets
US7074598B2 (en) * 2002-09-25 2006-07-11 Mayo Foundation For Medical Education And Research Detection of vancomycin-resistant enterococcus spp.
EP1601785A4 (en) * 2003-02-28 2007-02-07 Ravgen Inc Methods for detection of genetic disorders
US20040185446A1 (en) * 2003-03-18 2004-09-23 Jones Alison M. Cpn60 targets for quantification of microbial species
EP2390352A1 (en) 2003-03-18 2011-11-30 Quantum Genetics Ireland Limited Systems and methods for improving protein and milk production of dairy herds
US20050053950A1 (en) * 2003-09-08 2005-03-10 Enrique Zudaire Ubani Protocol and software for multiplex real-time PCR quantification based on the different melting temperatures of amplicons
WO2005112544A2 (en) 2004-02-19 2005-12-01 The Governors Of The University Of Alberta Leptin promoter polymorphisms and uses thereof
KR100738073B1 (en) 2004-09-01 2007-07-12 삼성전자주식회사 Method for quantification of initial nucleic acid concentration from real-time nucleic acid amplification data
US9464310B2 (en) * 2005-02-10 2016-10-11 Medical Diagnostic Laboratories, Llc Integrated method for collection and maintenance of detectability of a plurality of microbiological agents in a single clinical sample and for handling a plurality of samples for reporting a sum of diagnostic results for each sample
US20060246423A1 (en) * 2005-02-10 2006-11-02 Adelson Martin E Method and kit for the collection and maintenance of the detectability of a plurality of microbiological species in a single gynecological sample
US20060178838A1 (en) * 2005-02-10 2006-08-10 Adelson Martin E Method of receiving and handling a plurality of clinical samples for reporting a sum of diagnostic results for each sample
US7794659B2 (en) 2005-03-10 2010-09-14 Gen-Probe Incorporated Signal measuring system having a movable signal measuring device
EP1897016B1 (en) * 2005-06-22 2012-08-01 Gen-Probe Incorporated Method and algorithm for quantifying polynucleotides
EP1929049B1 (en) 2005-07-25 2013-04-10 Alere San Diego, Inc. Methods for multiplexing recombinase polymerase amplification
EP1798542A1 (en) * 2005-12-19 2007-06-20 Roche Diagnostics GmbH Analytical method and instrument
US20070238093A1 (en) * 2006-04-11 2007-10-11 Espy Mark J Detection of influenza A virus
US20070238095A1 (en) * 2006-04-11 2007-10-11 Mayo Foundation For Medical Education And Research , A Minnesota Corporation Detection of Influenza A Virus
EP3088533B1 (en) 2006-05-04 2018-01-17 Alere San Diego, Inc. Recombinase polymerase amplification
JP4993260B2 (en) * 2006-05-15 2012-08-08 石川県 Gene quantification apparatus and quantification method
US7802215B2 (en) * 2006-06-06 2010-09-21 Fujitsu Limited System and method for providing an improved sliding window scheme for clock mesh analysis
US7629124B2 (en) 2006-06-30 2009-12-08 Canon U.S. Life Sciences, Inc. Real-time PCR in micro-channels
US20080178653A1 (en) * 2007-01-29 2008-07-31 Applera Corporation System and Method for Interpolative Calibration
US20090055243A1 (en) 2007-08-21 2009-02-26 Jayson Lee Lusk Systems and methods for predicting a livestock marketing method
US8386184B2 (en) * 2007-08-28 2013-02-26 Becton, Dickinson And Company Systems and methods for determining an amount of starting reagent using the polymerase chain reaction
US8380457B2 (en) 2007-08-29 2013-02-19 Canon U.S. Life Sciences, Inc. Microfluidic devices with integrated resistive heater electrodes including systems and methods for controlling and measuring the temperatures of such heater electrodes
WO2009064717A2 (en) * 2007-11-14 2009-05-22 Talecris Biotherapeutics, Inc. Method for adjusting results of a polymerase chain reaction (pcr) instrument
KR101413659B1 (en) 2007-12-06 2014-07-01 삼성전자주식회사 Method of determining initial concentration of nucleic acids in a sample using a real time amplification data
US20100055733A1 (en) * 2008-09-04 2010-03-04 Lutolf Matthias P Manufacture and uses of reactive microcontact printing of biomolecules on soft hydrogels
EP2166470A1 (en) * 2008-09-19 2010-03-24 Corbett Research Pty Ltd Analysis of melt curves, particularly dsDNA and protein melt curves
US20100119454A1 (en) * 2008-11-03 2010-05-13 Ping Shen Use of the conserved Drosophila NPFR1 system for uncovering interacting genes and pathways important in nociception and stress response
GB0905325D0 (en) 2009-03-30 2009-05-13 Selex Sensors & Airborne Sys Detection system
WO2010135310A1 (en) 2009-05-20 2010-11-25 Biosite Incorporated Dna glycosylase/lyase and ap endonuclease substrates
US9057097B2 (en) 2009-06-05 2015-06-16 Alere San Diego Inc. Recombinase polymerase amplification reagents and kits
WO2011017567A1 (en) * 2009-08-05 2011-02-10 Life Technologies Corporation Methods for the analysis of proximity binding assay data
US8023794B2 (en) * 2009-08-10 2011-09-20 The Boeing Company Apparatus and method for establishing an optical path spanning a discontinuity in an optical channel
US9046507B2 (en) 2010-07-29 2015-06-02 Gen-Probe Incorporated Method, system and apparatus for incorporating capacitive proximity sensing in an automated fluid transfer procedure
CN103403533B (en) 2011-02-24 2017-02-15 简.探针公司 Systems and methods for distinguishing optical signals of different modulation frequencies in an optical signal detector
AU2012239989B2 (en) 2011-04-07 2016-03-10 Abbott Diagnostics Scarborough, Inc. Monitoring recombinase polymerase amplification mixtures
EP2756101B1 (en) 2011-09-15 2018-05-23 David A. Shafer Probe: antiprobe compositions for high specificity dna or rna detection
EP2809767B1 (en) 2012-02-03 2017-11-08 Axxin Pty Ltd Nucleic acid amplification and detection apparatus
AU2015275310B2 (en) * 2012-02-03 2018-10-18 Axxin Pty Ltd Nucleic acid amplification and detection apparatus and method
WO2014000037A1 (en) 2012-06-26 2014-01-03 Axxin Pty Ltd Nucleic acid amplification and detection kit
EP2877952B1 (en) 2012-07-27 2023-04-12 Gen-Probe Incorporated Dual reference calibration method and system for quantifying polynucleotides
CA2880583C (en) 2012-09-10 2021-02-16 Randy P. Rasmussen Multiple amplification cycle detection
US9714447B2 (en) 2013-08-19 2017-07-25 General Electric Company Detection of nucleic acid amplification in a porous substrate
EP3218698B1 (en) 2014-11-14 2020-05-13 Axxin Pty Ltd Biological sample collection and storage assembly
AU2016295422B2 (en) 2015-07-17 2022-01-06 Axxin Pty Ltd Diagnostic test assembly, apparatus, method
US10328147B2 (en) 2017-03-24 2019-06-25 Board Of Supervisors Of Louisiana State University And Agricultural And Mechanical College Herpes simplex virus type-1(HSV-1) vaccine strain VC2 generating an anti-EHV-1 immune response
EP3655549B1 (en) 2017-07-21 2023-01-04 BioFire Diagnostics, LLC Detection using concurrent melting curves
JP7276716B2 (en) 2017-09-27 2023-05-18 アックスイン ピーティーワイ リミテッド Diagnostic test system and method
GB201812192D0 (en) 2018-07-26 2018-09-12 Ttp Plc Variable temperature reactor, heater and control circuit for the same
EP3895170A4 (en) * 2018-12-14 2022-10-12 Seegene, Inc. Method for detecting a target analyte in a sample using an s-shaped function for a slope data set
KR102407585B1 (en) * 2020-08-21 2022-06-10 한림대학교 산학협력단 Test method for negative determination of polymerase chain reaction using artificial intelligence
KR102407584B1 (en) * 2020-08-21 2022-06-10 한림대학교 산학협력단 Test method for positive determination of polymerase chain reaction using artificial intelligence
CN112816446B (en) * 2020-12-24 2022-02-01 四川长虹电器股份有限公司 Method for detecting powder decay of fluorescent wheel based on fluorescence spectrum
WO2024110392A1 (en) 2022-11-22 2024-05-30 bioMérieux New nucleic acid binding compounds and uses

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4592365A (en) * 1981-08-10 1986-06-03 Ivac Corporation Electronic sphygmomanometer
US5455175A (en) 1990-06-04 1995-10-03 University Of Utah Research Foundation Rapid thermal cycling device
ATE309388T1 (en) 1996-06-04 2005-11-15 Univ Utah Res Found DEVICE AND METHOD FOR FLUORESCENCE MEASUREMENT OF DNA AMPLIFICATION
DK1179600T3 (en) 1996-06-04 2005-09-05 Univ Utah Res Found Hybridization monitoring during PCR
US5736333A (en) * 1996-06-04 1998-04-07 The Perkin-Elmer Corporation Passive internal references for the detection of nucleic acid amplification products
US6303305B1 (en) * 1999-03-30 2001-10-16 Roche Diagnostics, Gmbh Method for quantification of an analyte
US6387621B1 (en) * 1999-04-27 2002-05-14 University Of Utah Research Foundation Automated analysis of real-time nucleic acid amplification

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080176221A1 (en) * 2002-11-15 2008-07-24 Gen-Probe Incorporated Assay and compositions for detection of bacillus anthracis nucleic acid
US7494774B2 (en) 2002-11-15 2009-02-24 Gen-Probe Incorporated Assay and compositions for detection of Bacillus anthracis nucleic acid
US20090035779A1 (en) * 2007-06-29 2009-02-05 Roche Molecular Systems, Inc. Systems and methods for determining cross-talk coefficients in pcr and other data sets
US8095322B2 (en) * 2007-06-29 2012-01-10 Roche Molecular Systems, Inc. Systems and methods for determining cross-talk coefficients in PCR and other data sets
US8306754B2 (en) 2007-12-26 2012-11-06 Arkray, Inc. Nucleic acid amplification determining method and nucleic acid amplification determining device
US20120282611A1 (en) * 2009-10-21 2012-11-08 Brandeis University Methods, kits and reaction mixtures for analyzing single-stranded nucleic acid sequences
US9353407B2 (en) * 2009-10-21 2016-05-31 Brandeis University Methods, kits and reaction mixtures for analyzing single-stranded nucleic acid sequences
US10570462B2 (en) 2009-10-21 2020-02-25 Brandeis University Kits and reaction mixtures for analyzing single-stranded nucleic acid sequences
US20130273547A1 (en) * 2012-04-16 2013-10-17 Samsung Techwin Co., Ltd. Method to determine and correct baseline and to characterize pcr amplification kinetics
CN111426661A (en) * 2020-04-09 2020-07-17 基蛋生物科技股份有限公司 Fluorescence data acquisition and discrimination method in nucleic acid amplification stage
CN112669910A (en) * 2020-12-30 2021-04-16 杭州博日科技股份有限公司 Amplification curve baseline determination method and device and electronic equipment

Also Published As

Publication number Publication date
US20020123062A1 (en) 2002-09-05
ATE366821T1 (en) 2007-08-15
JP4689004B2 (en) 2011-05-25
EP1059523A3 (en) 2003-05-21
DE60035459T2 (en) 2008-03-20
DE60035459D1 (en) 2007-08-23
US6387621B1 (en) 2002-05-14
EP1059523B1 (en) 2007-07-11
EP1059523A2 (en) 2000-12-13
JP2000333700A (en) 2000-12-05

Similar Documents

Publication Publication Date Title
US6387621B1 (en) Automated analysis of real-time nucleic acid amplification
US7373253B2 (en) Multi-test analysis of real-time nucleic acid amplification
EP0686699B1 (en) Apparatus and method for determining the concentration of the target nucleic acid in PCR
US7363168B2 (en) Adaptive baseline algorithm for quantitative PCR
US7630837B2 (en) Real-time gene quantification with internal standards
US8182995B2 (en) Real time gene expression profiling
US6232079B1 (en) PCR method for nucleic acid quantification utilizing second or third order rate constants
EP0640828B1 (en) Monitoring multiple reactions simultaneously and analyzing same
US20030219754A1 (en) Fluorescence polarization detection of nucleic acids
EP1335028B1 (en) Multi-test analysis of real-time nucleic acid amplifications
RU2691763C2 (en) Device for simultaneous monitoring in real time scale of nucleic acid amplifications

Legal Events

Date Code Title Description
AS Assignment

Owner name: UTAH RESEARCH FOUNDATION, UNIVERSITY OF, UTAH

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:WITTWER, CARL T.;REEL/FRAME:010788/0762

Effective date: 20000419

STCF Information on status: patent grant

Free format text: PATENTED CASE

FPAY Fee payment

Year of fee payment: 4

FEPP Fee payment procedure

Free format text: PAYOR NUMBER ASSIGNED (ORIGINAL EVENT CODE: ASPN); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY

FPAY Fee payment

Year of fee payment: 8

FPAY Fee payment

Year of fee payment: 12