# A Bayesian method for calculating real-time quantitative PCR calibration curves using absolute plasmid DNA standards

- Mano Sivaganesan
^{1}Email author, - Shawn Seifring
^{2}, - Manju Varma
^{2}, - Richard A Haugland
^{2}and - Orin C Shanks
^{1}

**9**:120

**DOI: **10.1186/1471-2105-9-120

© Sivaganesan et al; licensee BioMed Central Ltd. 2008

**Received: **15 August 2007

**Accepted: **25 February 2008

**Published: **25 February 2008

## Abstract

### Background

In real-time quantitative PCR studies using absolute plasmid DNA standards, a calibration curve is developed to estimate an unknown DNA concentration. However, potential differences in the amplification performance of plasmid DNA compared to genomic DNA standards are often ignored in calibration calculations and in some cases impossible to characterize. A flexible statistical method that can account for uncertainty between plasmid and genomic DNA targets, replicate testing, and experiment-to-experiment variability is needed to estimate calibration curve parameters such as intercept and slope. Here we report the use of a Bayesian approach to generate calibration curves for the enumeration of target DNA from genomic DNA samples using absolute plasmid DNA standards.

### Results

Instead of the two traditional methods (classical and inverse), a Monte Carlo Markov Chain (MCMC) estimation was used to generate single, master, and modified calibration curves. The mean and the percentiles of the posterior distribution were used as point and interval estimates of unknown parameters such as intercepts, slopes and DNA concentrations. The software WinBUGS was used to perform all simulations and to generate the posterior distributions of all the unknown parameters of interest.

### Conclusion

The Bayesian approach defined in this study allowed for the estimation of DNA concentrations from environmental samples using absolute standard curves generated by real-time qPCR. The approach accounted for uncertainty from multiple sources such as experiment-to-experiment variation, variability between replicate measurements, as well as uncertainty introduced when employing calibration curves generated from absolute plasmid DNA standards.

## Background

The goal for many real-time quantitative PCR (qPCR) assays with clinical, forensic, or environmental applications is to develop a standardized method that can be implemented on an inter-laboratory scale. Real-time qPCR assays are ideal for such applications due to high levels of precision, specificity, and sensitivity. Real-time PCR allows for the continuous monitoring of PCR product production as the reaction occurs. Under ideal conditions these products accumulate exponentially in the reactions, i.e. their quantities double with each thermal cycle. Thus, real-time qPCR can be applied to determine a fixed threshold where the accumulation of PCR product is first significantly detectable over a real-time measurement background signal [for review see [1]]. The fractional cycle number where PCR product accumulation passes this fixed threshold is called the threshold cycle (C_{T}) [2]. qPCR is based on the theoretical premise that there is a log-linear relationship between the starting amount of DNA target in the reaction and the C_{T} value that is obtained. The C_{T} value can then be used to estimate the initial concentration of a DNA target from an unknown sample.

### Relative and Absolute Quantification with Real-Time qPCR

Two general strategies are often used to estimate DNA concentration from C_{T} values including relative and absolute approaches [3]. A relative quantification approach measures the change in target DNA concentration relative to another reference sample. This approach is ideal in gene expression studies where the goal is to measure the regulation of a gene in response to a particular treatment. However, a relative approach can be limiting for qPCR applications designed to quantify DNA targets with no clear connections to a reference target such as assays where the DNA target is from an uncharacterized microorganism. Relative quantification based qPCR methods can also be difficult to apply on an inter-laboratory scale for the enumeration of DNA targets from highly variable, complex, and poorly described sample matrices such as gastrointestinal and environmental samples [4].

Absolute quantification is another widely used strategy. Absolute quantification is achieved by using a standard curve, constructed by amplifying known amounts of target DNA in a parallel set of reactions [5]. Absolute quantification requires that the exact quantity of a standard is determined by independent means using spectrophotometry or an intercalating dye such as PicoGreen^{®} [6]. For bacterial DNA targets, genomic DNA from pure cell cultures is preferred. Cultivated bacterial cells can be isolated and counted to provide a conversion factor between mass of genomic DNA, copies of target DNA, and number of cells. However, this practice imposes a substantial restriction on the development of real-time qPCR methods targeting bacterial genes because an estimated 99% of the microbial diversity on the planet has not been cultivated [7–10]. When a DNA target originates from an uncultivated microorganism, plasmid DNA standards are often used. Plasmid preparations are advantageous because these preparations generate high quality, pure, and concentrated standards that can be independently quantified and converted to number of copies of target DNA. For absolute quantification approaches, an assumption must be made that plasmid and genomic DNA amplify with the same efficiency. Factors such as DNA stability, base composition, secondary structure, and presence of complex mixtures of non-target DNA could significantly alter amplification performance. A limited number of strategies have been used in an attempt to equilibrate these two types of DNA for real-time qPCR applications such as treating genomic DNA with a cocktail of restriction enzymes and DNA ultrasonication [11]. However, many studies simply assume that there are no differences.

In addition to the uncertainty associated with amplification of plasmid versus genomic DNA targets, there are a number of other sources of variability to consider when generating a calibration curve from absolute standards. Uncertainty can arise within and between experiments from numerous sources such as inconsistencies in quality of reagents, pipet calibration, as well as dilution preparation and storage of standards. Any of these factors could significantly alter C_{T} measurements from experiment to experiment. Therefore, estimation of uncertainty becomes critical to account for sources of variability and make reasonable estimates of calibration curve parameters.

### Estimating DNA Concentrations from C_{T} Values and Propagation of Uncertainty

Simple linear regression is commonly used to estimate DNA concentration from an unknown sample where the standard calibration model is developed with a DNA concentration (ie. plasmid copy number) and associated C_{T} measurements. Typically four to five known DNA concentrations are selected and then triplicate C_{T} measurements are taken at each DNA concentration to fit a calibration curve. The fitted curve is then used to estimate the mean DNA concentrations of unknown samples.

_{T}measurement as the dependent variable. Usually each experiment is repeated three to four times, with three replicates within each experiment. Even though triplicate C

_{T}measurements are taken at each DNA concentration of each experiment, the average of the C

_{T}measurements is commonly used to fit the calibration curve [12]. The corresponding regression model is given by:

_{i}is the average of the C

_{T}measurements at the ith DNA concentration, X

_{i}is the corresponding DNA concentration,

*α*and

*β*are regression coefficients and

*σ*

^{2}is the random error variance. For an unknown mean value of log

_{10}(X), say log

_{10}(X

_{0}), a Y value, say Y

_{0}is observed. The classical method uses Y

_{0}to estimate log

_{10}(X

_{0}) by:

*α*and

*β*, respectively. Finding the standard deviation of log

_{10}(${\widehat{X}}_{0}$) is not a simple statistical problem as it is a non-linear function of the estimated intercept and slope parameters. Thus for given X, a 100(1-

*α*)% confidence interval is constructed for $\widehat{Y}(=\widehat{\alpha}+\widehat{\beta}{\mathrm{log}}_{10}(X))$ first, as it is a linear function of intercept and slope parameters. The formula for this interval is given by:

where Z_{i} = log_{10}(X_{i}). Then the corresponding fiducial interval is reported as the confidence interval for X (given Y).

_{T}measurements from one experiment to obtain the calibration curve [13]. The corresponding regression model for replicated data is then given by:

where, Y_{ij} is the jth C_{T} measurement of ith DNA concentration. Except for more data points, the above regression model is same as the model given by Equation (1). The same least squares method is used to estimate the model parameters and then equation (2) is used to estimate unknown concentrations.

_{0}is given by:

*δ*

_{0}and

*δ*

_{1}. An approximate 100(1-

*α*)% confidence interval is given by :

An alternative approach to the classical and the inverse approaches is a Bayesian method using a Monte Carlo Markov Chain (MCMC) simulation technique. A detailed description of this method to generate a master calibration curve is discussed in the results and discussion section. Bayesian approaches have been employed in many molecular applications and have been particularly useful for microarray data analyses to account for multiple sources of uncertainty arising from experimental variation, background noise, and the use of multiple hybridization probes with different lengths and base pair compositions [15, 16]. Bayesian principles have also been used to model PCR amplification curves [17] and characterize the relationship between fluorescence chemistry and determination of C_{T} values during real-time detection [18].

Here we report the use of a Bayesian approach to generate calibration curves for the enumeration of target DNA from genomic DNA samples using absolute plasmid DNA standards. Calibration curves were generated from three independent real-time qPCR assays (Btheta, Entero1 and Entero2) using both genomic and plasmid DNA standards to test the assumption that both DNA types generate similar calibration curves. Finally, a calibration curve was generated for an additional real-time qPCR assay (HF183) where only a plasmid absolute standard was available. To account for potential differences in amplification performance between the plasmid standards and genomic DNA target from unknown samples, MCMC simulations were used to estimate the mean difference in slope and intercept from fitted curve equations for plasmid and genomic DNA produced from assays Btheta, Entero1, and Entero2. Using the same MCMC approach, these differences were applied to the plasmid DNA derived calibration curve for HF183. The modified calibration curve was then used to estimate DNA concentration from several unknown samples. The MCMC approach was ideal because it not only accounted for observed mean differences in plasmid and genomic DNA standards, but also propagated intra- and inter-assay variation.

## Results and discussion

### Bayesian Simulation Method

The Bayesian approach to statistical modeling is based on the premise that the uncertainty about unknown quantities, such as the parameters in a model, is described by a probability distribution; more precisely by a conditional probability distribution given all that is known, including the data as it becomes available. Initially, i.e., prior to obtaining the data, the uncertainty about the parameters are described by what is known as the prior distribution of the parameters, which probabilistically summarizes any available prior information about the parameters. Once the data is obtained and a suitable model for the observed data is chosen, the likelihood function of the parameters summarizing the information in the data can be mathematically expressed. The prior distribution is then combined with the likelihood via Bayes theorem, to obtain what is known as the posterior distribution of the parameters. The posterior distribution is a probabilistic expression of the (remaining) uncertainty about the parameters, after incorporating the available prior information and the information contained in the data. It is therefore the posterior distribution that forms the basis for Bayesian inference about the unknown parameters.

Typically, summaries of the posterior distribution such as the mean and the percentiles are used as point and interval estimates of an unknown parameter. In this paper, we use the term Bayesian credible interval (BCI) to refer to the interval with equal tail probability on either side under the posterior distribution. Closed form solutions for these quantities are usually not available, but, in most cases, MCMC methods [19–21] can be used to numerically compute the desired summaries of the posterior distribution. MCMC methods first use an iterative algorithm to generate a sequence of draws from a suitable Markov chain. Drawing a sufficiently long sequence, referred to as the burn-in phase, typically ensures convergence. Convergence is needed for the estimates of unknown model parameters. Examining the trace plots of the sample values of a model parameter provides evidence of when the simulation appears to have stabilized. Subsequent draws, after the burn-in phase, is a (Monte Carlo) sample from the posterior distribution, which can be used to calculate desired summaries of the posterior distribution.

The MCMC calculations in this study were done using the publicly available software WinBUGS [22]. Often, prior information about an unknown parameter may not be available. In such cases, standard non-informative prior distributions, i.e., probability distributions which contain little or no information about the parameters, are used, resulting in posterior distributions that are dominated by the likelihood. Some of the advantages of the Bayesian approach via MCMC are that it is capable of fitting models accounting for different sources of variability, and it allows for the appropriate processing of uncertainty when inference about complex functions of the model parameters are of interest. In such cases, the traditional methods tend to use approximations based on the basic summary values, i.e., estimates of model parameters and their standard errors, to obtain the inference, whereas the Bayesian approach via MCMC accurately evaluates the inference using the joint posterior distribution of the parameters. The Bayesian approach, however, also requires the specification of distributions of additional quantities in the models, as well as extensive simulation to fit them.

### Developing a Calibration Curve from a single qPCR experiment

A Bayesian approach was used to estimate the calibration curve parameters. To estimate X_{0}, we use all the triplicate C_{T} measurements from a single experiment to fit the calibration curve. The simple linear regression model given by the equation (3) was used here to fit the data. As no prior information is assumed for the model parameters *α*, *β* and *σ*^{2}, the following diffused prior distributions are used to estimate these model parameters:*α*, *β* ~ N (0, 10^{6})*σ*^{2} ~ Inv. Gamma(.0001,.0001).

*α*,

*β*and

*σ*

^{2}, given the data

*y*

_{ 1 }, ...,

*y*

_{ n }, is proportional to the likelihood, and the probability density of the prior distribution of

*α*,

*β*and

*σ*

^{2}. The MCMC method is employed using the WinBUGS software to obtain the required summaries of the posterior distributions of

*α*,

*β*and

*σ*

^{2}. For given Y

_{0}, the posterior distribution of

can be easily used to obtain summary statistics, such as mean, median and 95% BCI, for the unknown DNA concentration log_{10} (X_{0}).

### Developing a Master Calibration Curve from Multiple qPCR Experiments

_{ijk}is the kth Ct measurement of jth copy number and ith run, X

_{ij}is the jth copy number for ith run,

*α*

_{i}and

*β*

_{i}are regression coefficients for ith run, ${\sigma}_{\text{i}}^{2}$ is the random error variance of the ith calibration curve, $\overline{\alpha}$ and $\overline{\beta}$ are the overall regression coefficients, combining information from all runs. The following diffused prior distributions are used to estimate the model parameters:

*σ*

_{a}and

*σ*

_{b}, which is based on the harmonic mean of the estimated variances of the intercepts and slopes of individual calibration curves [23]. DuMouchel priors for

*σ*

_{a}and

*σ*

_{b}are given by:

where, U stands for the standard Uniform distribution U(0,1) and var(${\widehat{\alpha}}_{i}$) and var(${\widehat{\beta}}_{i}$) are respectively the estimated variances of the least squares estimates of *α*_{i} and *β*_{i}. The results obtained using the DuMouchel and Gamma priors for *σ*_{a} and *σ*_{b} are very similar. A MCMC simulation method was used to estimate the model parameters via WinBUGS software. Convergence diagnostics of Markov Chain draws from the posterior distributions of the parameters were checked using trace plots, auto-correlation plots, and Gelman and Rubin diagnostics [24, 25], and found to be satisfactory (data not shown).

_{0}, by requesting the posterior distribution of

_{10}(X

_{0}). Replacing $\overline{\alpha}$ by

*α*

_{i}and $\overline{\beta}$ by

*β*

_{i}in equation (10), we get the posterior distribution for the ith run (see Additional file 1). The estimated mean copy number corresponding to different C

_{T}measurements are plotted in Figure 1 for Entero2 genomic type (seven independent runs). Notice that the 95% upper and lower credible bounds and the fitted curve are for the copy number (in log base 10) in Figure 1. For comparison purposes, the averaged concentration data is used to obtain a fitted master curve and 95% BCI, for mean DNA concentration, and these are given in Figure 2 along with the corresponding 95% BCI using the raw data. It is better to use the raw data (than the averaged data) as it allows accounting for the within and between run variations in constructing credible interval for DNA concentration. Allowing for these (additional) variations would lead to more realistic and wider confidence intervals. Consequently, the 95% BCI is wider for the raw data than for the averaged data.

### Fitting a Genomic DNA Calibration Curve using Three Independent qPCR Assays

In real time qPCR studies using absolute standards, usually a calibration curve is developed to estimate an unknown DNA concentration. Typically, either plasmid or genome type calibration curves can be developed for a given assay. But, there are instances where PCR assays designed to target genomic DNA sequences must rely on plasmid derived absolute DNA standards to generate calibration curves such as PCR assays targeting genes from uncultivated microorganisms. qPCR assays that rely on plasmid absolute DNA standards to estimate genomic DNA concentrations from unknown samples must either assume that there is no difference in the amplification efficiencies between these two DNA types or estimate differences and account for this uncertainty in respective calibration curve statistics. A simulation method to estimate the genomic DNA type calibration curve for the assay HF183 using both plasmid and genomic DNA type curves of Btheta, Entero1 and Entero2 assays is discussed in this section.

where the DuMouchel priors for *σ*_{a1} and *σ*_{b1} are based on the least square estimates of *α*_{il} and *β*_{il}, respectively (see equation (10)).

### Estimating DNA Concentration from a Modified Master Calibration Curve

_{10}(X

_{0}), where

was used to estimate the mean, standard deviation and 95% credible intervals for unknown DNA concentration. Estimates for four unknown samples are given in the output section of Appendix B (see Additional file 2). Even though log_{10}(X_{0}) is a non-linear function of the parameters ${\overline{\alpha}}_{1},\mathrm{...}{\overline{\alpha}}_{7};{\overline{\beta}}_{1},\mathrm{...}{\overline{\beta}}_{7}$, the Bayesian MCMC simulation method can be easily applied to estimate X_{0}. To evaluate the impact of prior distributions, Uniform prior was used for each of *σ*_{a1}and *σ*_{b1} (l = 1...7). No apparent difference was seen in the resulting mean, median or 95% BCI of the two posterior distributions of any of the model parameters (data not shown).

## Conclusion

We employed a Bayesian approach for the estimation of DNA concentrations from environmental samples using absolute standard curves generated by real-time qPCR. Our approach allowed us to account for uncertainty from multiple sources such as experiment-to-experiment variation, variability between replicate measurements, as well as uncertainty introduced when employing calibration curves generated from absolute plasmid DNA standards. The Bayesian approach also allowed for the estimation of model parameters from multiple models simultaneously unlike stepwise progression of estimates typically used in real-time PCR calibration calculations. The flexible modeling capability of the Bayesian approach was ideal for real-time qPCR assays that rely on absolute plasmid DNA standards for quantification and this method should be applicable over a wide range of study designs.

## Methods

### Sample collection and DNA extraction

Select individual fecal and recreational water samples were collected as previously described [26]. All DNA extractions were performed with the FastDNA Kit for Soils (Q-Biogene; Carlsbad, CA) [26].

### Genomic DNA standard preparations from pure bacterial cultures

American Type Culture Collection (ATCC) bacterial strains were used to prepare genomic DNA calibration standards. *E. faecalis* (ATCC #29212) was cultured as previously described [27]. *B. thetaiotaomicron* (ATCC # 29741) cells were grown in chopped meat carbohydrate broth (Remel, Lenexa, KS) according to manufacturer's instructions. Both cultures were harvested by centrifugation at 8,000 × *g* for 5 min, washed twice using sterile phosphate buffered saline (Sigma, St. Louis, MO) and stored in aliquots at -40°C. Cell concentrations of each organism in the final washed suspensions were determined by bright field microscopy at 40× magnification in disposable hemocytometer chambers (Nexcelom Bioscience, #CP2-002, Lawrence, MA). DNA was isolated from the cell suspensions using a bead beating extraction approach [27] and incubated for one hour at 37°C with 0.017 *μ* g/*μ* l RNase A (Gentra Systems, USA). DNA purification was performed using a silica column adsorption kit (DNA-EZ, GeneRite, Kendall Park, NJ.). DNA concentrations of cell extracts were determined by spectrophotometric absorbance readings at 260 nm (A_{260}) and purity of the DNA preparations was determined by A_{260}/A_{280} ratios.

### Plasmid DNA standard preparation

^{®}TET labeled probe sequence flanked by PCR primer binding sites for all four qPCR assays was developed using overlap extension PCR [Figure 5, [28]]. To build the plasmid construct, long oligonucleotides (> 100 bp, Table 1) containing multiple primer sequences [29] were designed such that their 3' ends overlapped. Overlapping fragments were then combined into a single DNA molecule using a two step overlap extension PCR, i.e. the partially overlapping products of two initial overlap extension PCR experiments were combined by a second overlap extension PCR. The plasmid construct was then inserted into a pCR4

^{®}TOPO plasmid vector (Invitrogen) and the resulting recombinant plasmid was purified from transformed

*E. coli*cell cultures using a Qiagen Plasmid Purification Kit (Qiagen, Valencia, CA). Plasmid DNA was linearized by a Not1 restriction digestion (New England BioLabs, Beverly, MA), quantified with a NanoDrop ND-1000 UV spectrophotometer (NanoDrop Technologies), and diluted in 10 mM Tris, 0.1 mM EDTA, pH 8.0 to generate samples ranging from approximately 100 to 4 × 10

^{4}molecules. Dilutions were stored in TE buffer (10 mM Tris, 0.1 mM EDTA, pH 8.0) in single use aliquots.

Oligonucleotides and probe used in study.

Sequence 5' → 3' | Reference | |
---|---|---|

Btheta |
| [30] |

| ||

Entero1 |
| [31] |

| ||

Entero2 |
| [32] |

| ||

HF183 |
| |

| ||

UC1F1 | CCGTCATCCTTCACGCTACTGAGGACCGAACCCACGTACCCTTC | This Study |

AGTGCCGCAGTCGTTCCATTAGGCAGTTGGTGAGAAA | ||

UC1R1 | CCTGCCGTCTCGTGCTCCTCAAACGCTTCTTAGTCAGGTACCGT | |

CAAGTTCGTTTGGAATTTCTCACCAACTGCCTAATG | ||

UC1F2 | TGAGGAGCACGAGACGGCAGGAACCTTCCTCTCAGAACCCCAATG | |

ATGGAGGTAGAGCACTGACACGGTCCAAACTCCTA | ||

UC1R2 | GATCATGAGTTCACATGTCCGCGTCGCAGGATGTCAAGACAGTA | |

AATCCGGATAACGCTCGTAGGAGTTTGGACCGTGTCA | ||

UC1 | [TET] CCTGCCGTCTCGTGCTCCTCA [TAMRA]* | [35] |

### Quantitative real-time PCR

Four qPCR assays were used in this study including HF183, Btheta, Entero1, and Entero2 (Table 1) [30–33]. Amplification was performed in a 7900 HT Fast Real Time Sequence Detector (Applied Biosystems) with default thermal cycle conditions. Reaction mixtures (25 *μ* l) contained 1X TaqMan^{®} Universal PCR Master Mix with AmpErase^{®} uracil-N-glycosylase (UNG, Applied Biosystems), 0.2 mg/ml bovine serum albumin (Sigma), 1 *μ* M of each primer, 80 nM FAM™ or TET™ labeled TaqMan^{®} probe (Applied Biosystems), and either 2 ng genomic DNA (unknown samples) or 100 to 4 × 10^{4} target gene copies (plasmid or purified genomic DNA). All reactions were performed in triplicate. Data was initially analyzed with Sequence Detector Software (Version 2.2.2) at a threshold determination of 0.03. Threshold cycle (C_{T}) values were exported to Microsoft Excel for further statistical analysis.

### Data analysis

An analysis of Covariance (ANCOVA) model was used to compare the overall mean intercept and slope of genome standard curves with the corresponding overall mean intercept and slope of the corresponding plasmid standard curves. ANCOVA tests were performed using SAS programs (Cary, North Carolina) with the following procedure "PROC MIXED" [34]. Markov Chain Monte Carlo (MCMC) simulation method was used to obtain single, master, and modified calibration curves. Summaries of the posterior distribution such as the mean and the percentiles were used as point and interval estimates of unknown parameters of interest. The software WinBUGS versions 1.4.1 [22] was used to perform all simulations.

## Declarations

### Acknowledgements

Any opinions expressed in this paper are those of the author(s) and do not, necessarily, reflect the official positions and policies of the U.S. EPA and any mention of products or trade names does not constitute recommendation for use.

## Authors’ Affiliations

## References

- Sambrook J, Russell DW:
*Molecular Cloning: A Laboratory Manual.**Volume 1–3*. 4th edition. New York: Cold Spring Harbor Laboratory Press; 2001.Google Scholar - Higuchi R, Fockler C, Dollinger G, Watson R: Kinetic PCR analysis: real-time monitoring of DNA amplification reactions.
*Biotechnology*1993, 11: 1026–1030. 10.1038/nbt0993-1026View ArticlePubMedGoogle Scholar - ABI: Essentials of Real Time PCR.
*Applied Biosystems*2006.Google Scholar - Fierer N, Jackson JA, Vilgalys R, Jackson RB: Assessment of Soil Microbial Community Structure by Use of Taxon-Specific Quantitative PCR Assays.
*Applied and Environmental Microbiology*2005, 71(7):4117–4120. 10.1128/AEM.71.7.4117-4120.2005PubMed CentralView ArticlePubMedGoogle Scholar - ABI: Absolute Quantitation Using Standard Curve Getting Started Guide.
*Applied Biosystems*2006.Google Scholar - Singer VL, Jones LJ, Yue ST, Haugland RP: Characterization of PicoGreen Reagent and Development of a Fluorescence-Based Solution Assay for Double-Stranded DNA Quantitation.
*Analytical Biochemistry*1997, 249(2):228–238. 10.1006/abio.1997.2177View ArticlePubMedGoogle Scholar - Handlesman J: Metagenomics: Application of genomics to uncultured microorganisms.
*Microbiology and Molecular Biology Reviews*2004, 68(4):669–685. 10.1128/MMBR.68.4.669-685.2004View ArticleGoogle Scholar - Grimes DJ, Atwell RW, Brayton PR, Palmer LM, Rollins DM, Roszak DB, Singleton FL, Tamplin ML, Colwell RR: The fate of enteric pathogenic bacteria in estuarine and marine environments.
*Microbiology Science*1986, 3: 324–329.Google Scholar - Torsvik V, Goksoyer J, Daae FL: High diversity in DNA of soil bacteria.
*Applied and Environmental Microbiology*1990, 56: 782–787.PubMed CentralPubMedGoogle Scholar - Staley JT, Konopka A: Measurment of insitue activities of nonphotosynthetic microorganisms in aquatic and terrestrial habitats.
*Annual Reviews in Microbiology*1985, 39: 321–346. 10.1146/annurev.mi.39.100185.001541View ArticleGoogle Scholar - Toyota A, Akiyama H, Sugimura M, Watanabe T, Sakata K, Shiramasa Y, Kitta K, Hino A, Esaka M, Maitani T: Rapid Quantification Methods for Genetically Modified Maize Contents Using Genomic DNAs Pretreated by Sonication and Restriction Endonuclease Digestion for a Capillary-Type Real-Time PCR System with a Plasmid Reference Standard.
*Bioscience, Biotechnology, and Biochemistry*2006, 70(12):2965–2973. 10.1271/bbb.60366View ArticlePubMedGoogle Scholar - Martin B, Jofre A, Garriga M, Pla M, Aymerich T: Rapid quantitative detection of
*Lactobacillus sakei*in meat and fermented sausages by real-time PCR.*Applied and Environmental Microbiology*2006, 72(9):6040–6048. 10.1128/AEM.02852-05PubMed CentralView ArticlePubMedGoogle Scholar - Ibekwe AW, Watt PM, Grieve CM, Sharma VK, Lyons SR: Multiplex fluoregenic real-time PCR for detection and quantification of
*Escherichia coli*O157:H7 in dairy wastewater wetlands.*Applied and Environmental Microbiology*2002, 68(10):4853–4862. 10.1128/AEM.68.10.4853-4862.2002PubMed CentralView ArticlePubMedGoogle Scholar - Rutledge RG, Cote C: Mathematics of quantitative kinetic PCR and the application of standard curves.
*Nucleic Acids Research*2003, 31(16):e93. 10.1093/nar/gng093PubMed CentralView ArticlePubMedGoogle Scholar - Frigessi A, van de Wiel MA, Holden M, Svendsrud DH, Glad IK, Lyng H: Genome-wide estimation of transcript concentrations from spotted cDNA microarry data.
*Nucleic Acids Research*2005, 33(17):e143. 10.1093/nar/gni141PubMed CentralView ArticlePubMedGoogle Scholar - Conlon EM, Song JJ, Liu A: Bayesian meta-analysis models for microarray data: a comparative study.
*BMC Bioinformatics*2007, 8(80):1–21.Google Scholar - Lalam N:
*Statistical Inference for Quantitative Polymerase Chain Reaction Using a Hidden Markov Model: A Bayesian Approach.**Volume 6*. The Berkeley Electronic Press; 2007:1–33.Google Scholar - Lalam N, Jacob C: Bayesian Estimation for Quantification by Real-time Polymerase Chain Reaction Under a Branching Process Model of the DNA Molecules Amplification Process.
*Mathematical Population Studies*2007, 14(2):111–129. 10.1080/08898480701298418View ArticleGoogle Scholar - Brooks SP: Markov chain Monte Carlo method and its application.
*The Statistician*1998, 47: 69–100.Google Scholar - Gelfand AE, Smith AFM: Sampling based approaches to calculating marginal densities.
*Journal of the American Statistical Association*1990, 85: 398–409. 10.2307/2289776View ArticleGoogle Scholar - Gelman A, Carlin JC, Stern H, Rubin DB:
*Bayesian Data Analysis.*New York: Chapman & Hall; 1995.Google Scholar - The BUGS Project-Bayesian inference Using Gibbs Sampling[http://www.mrc-bsu.cam.ac.uk/bugs]
- Spiegelhalter DJ, Abrams KR, Myles JP:
*Bayesian Approaches to Clinical Trials and Health-Care Evaluation.*USA: John Wiley & Sons Inc; 2003.View ArticleGoogle Scholar - Gelman A, Rubin DB: Inference from iterative simulation using multiple sequences.
*Statistical Science*1992, 7: 457–511. 10.1214/ss/1177011136View ArticleGoogle Scholar - Cowles MK, Carlin BP: Markov Chain Monte Carlo Convergence Diagnostics: A Comparative Review.
*Journal of the American Statistical Association*1996, 91: 883–904. 10.2307/2291683View ArticleGoogle Scholar - Shanks OC, Santo Domingo JW, Lamendella R, Kelty CA, Graham JE: Competitive metagenomic DNA hybridization identifies host-specific microbial genetic markers in cow fecal samples.
*Applied and Environmental Microbiology*2006, 72(6):4054–4060. 10.1128/AEM.00023-06PubMed CentralView ArticlePubMedGoogle Scholar - Haugland RA, Siefring SC, Wymer LJ, Brenner KP, Dufour AP: Comparison of
*Enterococcus*measurements in freshwater at two recreational beaches by quantitative polymerase chain reaction and membrane filter culture analysis.*Water Research*2005, 39: 559–568. 10.1016/j.watres.2004.11.011View ArticlePubMedGoogle Scholar - Higuchi R, Krummel B, Saiki RK: A general method of in vitro preparation and specific mutagenesis of DNA fragments: study of protein and DNA interactions.
*Nucleic Acids Research*1988, 16: 7351–7367. 10.1093/nar/16.15.7351PubMed CentralView ArticlePubMedGoogle Scholar - Stocher M, Leb V, Berg J: A convenient approach to the generation of multiple internal control DNA for a panel of real-time PCR assays.
*Journal of Virology Methods*2003, 108: 1–8. 10.1016/S0166-0934(02)00266-5View ArticleGoogle Scholar - Blackwood AD, Noble RT: Development of a rapid quantitative PCR method for quantification of
*Bacteroides thetaiotaomicron*as an alternative indicator of fecal contamination.*In Preparation*2007.Google Scholar - Ludwig W, Schleifer KH: How quantitative is quantitative PCR with respect to cell counts?
*Systematic and Applied Microbiology*2000, 23(4):556–562.View ArticlePubMedGoogle Scholar - Seifring SC, Varma M, Atikovic E, Wymer LJ, Haugland RA: Improved real-time PCR assays for the detection of fecal indicator bacteria in surface waters with different instrument and reagent systems.
*Journal of Water and Health*2007, in press.Google Scholar - Bernhard AE, Field KG: A PCR Assay to Discriminate Human and Ruminant Feces on the Basis of Host Differences in
*Bacteroides-Prevotella*Genes Encoding for 16S rRNA.*Applied and Environmental Microbiology*2000, 66(10):4571–4574. 10.1128/AEM.66.10.4571-4574.2000PubMed CentralView ArticlePubMedGoogle Scholar - SAS:
*SAS/STAT User's Guide Version 6.*4th edition. Cary, USA: SAS Institute Inc; 1990.Google Scholar - Zhang Y, Zhang D, Wenquan L, Chen J, Peng Y, Cao W: A novel real-time quantitative PCR method using attached universal template probe.
*Nucleic Acids Research*2003, 31: e123. 10.1093/nar/gng123PubMed CentralView ArticlePubMedGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.