MBE Advance Access originally published online on March 29, 2007
Molecular Biology and Evolution 2007 24(6):1416-1422; doi:10.1093/molbev/msm062
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Articles |
Bayesian Estimation of Sequence Damage in Ancient DNA
* Department of Zoology, University of Oxford, Oxford, United Kingdom
Department of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom
E-mail: beth.shapiro{at}zoo.ox.ac.uk.
Accepted for publication March 19, 2007.
DNA extracted from archaeological and paleontological remains is usually damaged by biochemical processes postmortem. Some of these processes lead to changes in the structure of the DNA molecule, which can result in the incorporation of incorrect nucleotides during polymerase chain reaction. These base misincorporations, or miscoding lesions, can lead to the inclusion of spurious additional mutations in ancient DNA (aDNA) data sets. This has the potential to affect the outcome of phylogenetic and population genetic analyses, including estimates of mutation rates and genetic diversity. We present a novel model, termed the delta model, which estimates the amount of damage in DNA data and accounts for its effects in a Bayesian phylogenetic framework. The ability of the delta model to estimate damage is first investigated using a simulation study. The model is then applied to 13 aDNA data sets. The amount of damage in these data sets is shown to be significant but low (about 1 damaged base per 750 nt), suggesting that precautions for limiting the influence of damaged sites, such as cloning and enzymatic treatment, are worthwhile. The results also suggest that relatively high rates of mutation previously estimated from aDNA data are not entirely an artifact of sequence damage and are likely to be due to other factors such as the persistence of transient polymorphisms. The delta model appears to be particularly useful for placing upper credibility limits on the amount of sequence damage in an alignment, and this capacity might be beneficial for future aDNA studies or for the estimation of sequencing errors in modern DNA.
Key Words: miscoding lesions • delta model • sequence errors • mutation rate
Connie Mulligan, Associate Editor
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  A. Rambaut, S. Y.W. Ho, A. J. Drummond, and B. Shapiro Accommodating the Effect of Ancient DNA Damage on Inferences of Demographic Histories Mol. Biol. Evol., February 1, 2009; 26(2): 245 - 248. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
E. Axelsson, E. Willerslev, M. T. P. Gilbert, and R. Nielsen The Effect of Ancient DNA Damage on Inferences of Demographic Histories Mol. Biol. Evol., October 1, 2008; 25(10): 2181 - 2187. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
L. M. Mateiu and B. H. Rannala Bayesian Inference of Errors in Ancient DNA Caused by Postmortem Degradation Mol. Biol. Evol., July 1, 2008; 25(7): 1503 - 1511. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 M. T. P. Gilbert, D. I. Drautz, A. M. Lesk, S. Y. W. Ho, J. Qi, A. Ratan, C.-H. Hsu, A. Sher, L. Dalen, A. Gotherstrom, et al. Intraspecific phylogenetic analysis of Siberian woolly mammoths using complete mitochondrial genomes PNAS, June 17, 2008; 105(24): 8327 - 8332. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
C. P. Burridge, D. Craw, D. Fletcher, and J. M. Waters Geological Dates and Molecular Rates: Fish DNA Sheds Light on Time Dependency Mol. Biol. Evol., April 1, 2008; 25(4): 624 - 633. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 S. Y. W. Ho, B. Shapiro, M. J. Phillips, A. Cooper, and A. J. Drummond Evidence for Time Dependency of Molecular Rate Estimates Syst Biol, June 1, 2007; 56(3): 515 - 522. [Full Text] [PDF]
|
|