Automatic threshold setting and baseline determination for real-time PCR
View Patent ↗The invention discloses a system and methods for quantitating the presence of nucleic acid sequences by evaluation of amplification data generated using real-time PCR. In one aspect, the methods may be adapted to identify a threshold and threshold cycle for one or more reactions based upon evaluation of exponential and baseline regions for each amplification reaction. The methodology used in the analysis may be readily automated such that subjective user interpretation of the data is substantially reduced or eliminated.
1. A method for target sequence quantitation by a biological analysis system comprising:
acquiring signal information of an amplification of a plurality of target sequences amplified in the presence of a detectable reporter construct from a detector configured to detect the associated detectable reporter construct, the detectable reporter construct comprising of a fluorescent label molecule and a quencher molecule, and the detector comprising of a charge coupled device configured to detect an emission spectra associated with the fluorescent label molecule;
identifying, by the biological analysis system, an exponential region comprising increases in the signal information for each target sequence amplification, each exponential region comprising an upper bound and a lower bound;
identifying, by the biological analysis system, a baseline component for the signal information located prior to the lower bound of the exponential region associated with each target sequence amplification to normalize the signal information for each target sequence to compensate for non-specific fluorescence in the signal information;
determining, by the biological analysis system, an exponential region threshold for the target sequence amplifications using a comparison of the upper bounds and lower bounds of the target sequence amplifications;
determining, by the biological analysis system, a threshold cycle corresponding to an initial quantity of each target sequence using the exponential region threshold; and
determining, by the biological analysis system, the initial quantity of at least one target sequence using the determined threshold cycle,
wherein the biological analysis system comprises of a reaction module, the reaction module configured to release the quencher molecule to allow the emission spectra associated with the fluorescent label molecule to be detected by the detector.
2. The method of claim 1 , wherein determining a exponential region threshold comprises:
identifying a minimal upper bound from the upper bounds of each of the target sequence amplifications;
identifying a maximal lower bound from the lower bounds of each of the target sequence amplifications; and
identifying the intersection between the minimal upper bound and maximal lower bound as the exponential region threshold.
3. The method of claim 2 , further comprising:
identifying an exponential region for each of the target sequence amplifications based on the exponential region threshold determined from the plurality of target sequence amplifications.
4. The method of claim 3 , further comprising performing a polynomial fitting operation to fit a polynomial to at least an interval of the exponential region including the exponential region threshold for each of the target sequence amplifications.
5. The method of claim 4 , further comprising identifying a threshold cycle for each of the target sequence amplifications by calculating a root of the associated polynomial of each of the target sequence amplifications.
6. The method of claim 1 , wherein the identification of the exponential region lower bound comprises detecting incremental differences in the signal information starting at the exponential region upper bound of each of the target sequence amplifications.
7. The method of claim 6 , wherein the incremental differences comprise at least one of differences in a ratio of detected intensity of successive values of the signal information and differences in detected intensity of successive values of the signal information.
8. The method of claim 1 , wherein the normalizing comprises:
identifying a baseline region based on the exponential region lower bound for each of the target sequence amplifications;
performing a regression analysis of the baseline region to generate a characteristic equation for each of the target sequence amplifications; and
differencing the baseline component using the characteristic equation to normalize the signal information for each of the target sequence amplifications.
9. A biological analysis system for target sequence quantitation, comprising:
a detector configured to detect reporter label intensities generated by amplifying a plurality of target sequences in the presence of a detectable reporter construct, the detectable reporter construct comprising of a fluorescent label molecule and a quencher molecule, and the detector comprising of a charge coupled device configured to detect an emission spectra associated with the fluorescent label molecule;
a reaction module configured to release the quencher molecule to allow the emission spectra associated with the fluorescent label molecule to be detected by the detector;
a non-transitory computer-readable storage medium having stored thereon a set of instructions;
a data processor coupled to the detector and the storage medium, the processor configured to execute the instructions stored on the medium and, upon execution of the instructions, cause the system to perform the steps of:
identifying an exponential region associated with substantial increases in the signal information for each target sequence amplification, each exponential region comprising an upper bound and a lower bound;
identifying a baseline component for the signal information located prior to the lower bound of the exponential region associated with each target sequence amplification;
identifying the signal information for each target sequence amplification using the baseline component to compensate for non-specific fluorescence in the signal information;
determining an exponential region threshold for the target sequence amplifications using a comparison of the upper bounds and lower bounds of the target sequence amplifications;
determining a threshold cycle related to an initial quantity of each target sequence from the exponential region threshold; and
determining the initial quantity of at least one target sequence using the determined threshold cycle.
10. The system of claim 9 , wherein determining the exponential region threshold comprises:
identifying a minimal upper bound from the upper bounds of each of the target sequence amplifications;
identifying a maximal lower bound from the lower bounds of each of the target sequence amplifications; and
identifying the intersection between the minimal upper bound and maximal lower bound as the combined exponential region threshold.
11. The system of claim 9 , the data processor being further configured to:
identify an exponential region for each of the target sequence amplifications based on the exponential region threshold determined from the plurality of target sequence amplifications.
12. The system of claim 9 , the data processor being further configured to perform a polynomial fitting operation to fit a polynomial to at least an interval of the exponential region including the exponential region threshold for each of the target sequence amplifications.
13. The system of claim 12 , wherein the data processor is further configured to identify a threshold cycle for each of the target sequence amplifications by calculating a root of the associated polynomial of each of the target sequence amplifications.