A stochastic model of the processes in PCR based amplification of STR DNA in forensic applications.

In forensic DNA profiling use is made of the well-known technique of PCR. When the amount of DNA is high, generally unambiguous profiles can be obtained, but for low copy number DNA stochastic effects can play a major role. In order to shed light on these stochastic effects, we present a simple model for the amplification process. According to the model, three possible things can happen to an individual single DNA strand in each complete cycle: successful amplification, no amplification, or amplification with the introduction of stutter. The model is developed in mathematical terms using a recursive approach: given the numbers of chains at a given cycle, the numbers in the next can be described using a multinomial probability distribution. A full set of recursive relations is derived for the expectations and (co)variances of the number of amplicon chains with no, 1 or 2 stutters. The exact mathematical solutions of this set are given, revealing the development of the expectations and (co)variances as function of the cycle number. The equations reveal that the expected number of amplicon chains without stutter grows exponentially with the cycle number, but for the chains with stutter the relation is more complex. The relative standard deviation on the numbers of chains (coefficient of variation) is inversely proportional to the square root of the expected number of DNA strands entering the amplification. As such, for high copy number DNA the stochastic effects can be ignored, but they play an important role at low concentrations. For the allelic peak, the coefficient of variation rapidly stabilizes after a few cycles, but for the chains with stutter the decrease is more slowly. Further, the ratio of the expected intensity of the stutter peak over that of the allelic peak increases linearly with the number of cycles. Stochastic models, like the one developed in the current paper, can be important in further developing interpretation rules in a Bayesian context.

[1]  Fengzhu Sun,et al.  Taq DNA polymerase slippage mutation rates measured by PCR and quasi-likelihood analysis: (CA/GT)n and (A/T)n microsatellites. , 2003, Nucleic acids research.

[2]  D. Ubelaker.  Review of: The Scientific Study of Mummies , 2004 .

[3]  S L Lauritzen,et al.  Probabilistic expert systems for handling artifacts in complex DNA mixtures. , 2011, Forensic science international. Genetics.

[4]  Robert G Cowell,et al.  Validation of an STR peak area model. , 2009, Forensic science international. Genetics.

[5]  P. Walsh,et al.  Sequence analysis and characterization of stutter products at the tetranucleotide repeat locus vWA. , 1996, Nucleic acids research.

[6]  Dennis Y. Wang,et al.  Development and Validation of the AmpFℓSTR® MiniFilerTM PCR Amplification Kit: A MiniSTR Multiplex for the Analysis of Degraded and/or PCR Inhibited DNA * , 2008, Journal of forensic sciences.

[7]  J Buckleton,et al.  An investigation of the rigor of interpretation rules for STRs derived from less than 100 pg of DNA. , 2000, Forensic science international.

[8]  John Buckleton,et al.  Validation and development of interpretation guidelines for low copy number (LCN) DNA profiling in New Zealand using the AmpFlSTR SGM Plus multiplex. , 2010, Forensic science international. Genetics.

[9]  Andrew J Gibb,et al.  Characterisation of forward stutter in the AmpFlSTR SGM Plus PCR. , 2009, Science & justice : journal of the Forensic Science Society.

[10]  Benoît Leclair,et al.  Systematic analysis of stutter percentages and allele peak height and peak area ratios at heterozygous STR loci for forensic casework and database samples. , 2004, Journal of forensic sciences.

[11]  James Curran,et al.  A graphical simulation model of the entire DNA process associated with the analysis of short tandem repeat loci , 2005, Nucleic acids research.

[12]  Maximizing DNA profiling success from sub-optimal quantities of DNA: a staged approach. , 2009, Forensic science international. Genetics.

[13]  W R Mayr,et al.  DNA commission of the International Society of Forensic Genetics: Recommendations on the interpretation of mixtures. , 2006, Forensic science international.

[14]  P Gill,et al.  A comparison of the characteristics of profiles produced with the AMPFlSTR SGM Plus multiplex system for both standard and low copy number (LCN) STR DNA analysis. , 2001, Forensic science international.

[15]  Pieter A W M Wouters,et al.  Stochastic processes defining sensitivity and variability of internally calibrated quantitative NASBA-based viral load assays. , 2002, Nucleic acids research.

[16]  P Gill,et al.  Report of the European Network of Forensic Science Institutes (ENSFI): formulation and testing of principles to evaluate STR multiplexes. , 2000, Forensic science international.