A Bayesian method for calculating real-time quantitative PCR calibration curves using absolute plasmid DNA standards
© Sivaganesan et al; licensee BioMed Central Ltd. 2008
Received: 15 August 2007
Accepted: 25 February 2008
Published: 25 February 2008
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.
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.
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.
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 ]. The fractional cycle number where PCR product accumulation passes this fixed threshold is called the threshold cycle (CT) . 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 CT value that is obtained. The CT 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 CT values including relative and absolute approaches . 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 .
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 . Absolute quantification requires that the exact quantity of a standard is determined by independent means using spectrophotometry or an intercalating dye such as PicoGreen® . 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 . 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 CT 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 CT 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 CT measurements. Typically four to five known DNA concentrations are selected and then triplicate CT 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.
where Zi = log10(Xi). Then the corresponding fiducial interval is reported as the confidence interval for X (given Y).
where, Yij is the jth CT 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.
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  and characterize the relationship between fluorescence chemistry and determination of CT values during real-time detection .
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 . 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 X0, we use all the triplicate CT 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, 106)σ2 ~ Inv. Gamma(.0001,.0001).
can be easily used to obtain summary statistics, such as mean, median and 95% BCI, for the unknown DNA concentration log10 (X0).
Developing a Master Calibration Curve from Multiple qPCR Experiments
where, U stands for the standard Uniform distribution U(0,1) and var() and var() 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).
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
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 log10(X0) is a non-linear function of the parameters , the Bayesian MCMC simulation method can be easily applied to estimate X0. To evaluate the impact of prior distributions, Uniform prior was used for each of σa1and σ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).
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.
Sample collection and DNA extraction
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 . 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  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 (A260) and purity of the DNA preparations was determined by A260/A280 ratios.
Plasmid DNA standard preparation
Oligonucleotides and probe used in study.
Sequence 5' → 3'
[TET] CCTGCCGTCTCGTGCTCCTCA [TAMRA]*
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 × 104 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 (CT) values were exported to Microsoft Excel for further statistical 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" . 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  was used to perform all simulations.
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.
- 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
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.