ML Estimation of DNA Initial Copy Number in Polymerase Chain Reaction (PCR) Processes

Estimation of DNA copy number in a given biological sample is an extremely important problem in genomics. This problem is especially challenging when the number of the DNA strands is minuscule, which is often the case in applications such as pathogen and genetic mutation detection. A recently developed technique, real-time polymerase chain reaction (PCR), amplifies the number of initial target molecules by replicating them through a series of thermal cycles. Ideally, the number of target molecules doubles at the end of each cycle. However, in practice, due to biochemical noise the efficiency of the PCR reaction, defined as the fraction of target molecules which are successfully copied during a cycle, is always less than 1. In this paper, we formulate the problem of joint maximum-likelihood estimation of the PCR efficiency and the initial DNA copy number. As indicated by simulation studies, the performance of the proposed estimator is superior with respect to competing statistical approaches. Moreover, we compute the Cramer-Rao lower bound on the mean-square estimation error.

[1]  K. Mullis,et al.  Specific synthesis of DNA in vitro via a polymerase-catalyzed chain reaction. , 1987, Methods in enzymology.

[2]  Arjang Hassibi,et al.  A STOCHASTIC MODEL AND SIMULATION ALGORITHM FOR POLYMERASE CHAIN REACTION ( PCR ) SYSTEMS , 2004 .

[3]  G Stolovitzky,et al.  Efficiency of DNA replication in the polymerase chain reaction. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[4]  J. Peccoud,et al.  Estimation of the parameters of a branching process from migrating binomial observations , 1998, Advances in Applied Probability.

[5]  Cheng Zhao,et al.  Estimation of the Mutation Rate During Error-prone Polymerase Chain Reaction , 2000, J. Comput. Biol..

[6]  R. Abramson,et al.  Detection of specific polymerase chain reaction product by utilizing the 5'----3' exonuclease activity of Thermus aquaticus DNA polymerase. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Nadia Lalam,et al.  Modelling the PCR amplification process by a size-dependent branching process and estimation of the efficiency , 2004, Advances in Applied Probability.

[8]  Babak Hassibi,et al.  ML Estimation of DNA Initial Copy Number in Polymerase Chain Reaction (PCR) Processes , 2007, 2007 IEEE International Conference on Acoustics, Speech and Signal Processing - ICASSP '07.

[9]  Estimation of the mean and the initial probabilities of a branching process , 1974 .

[10]  Marek Kimmel,et al.  Branching processes in biology , 2002 .

[11]  J. Sninsky,et al.  PCR Applications: Protocols for Functional Genomics , 1999 .