The effect of ancient DNA damage on inferences of demographic histories.

The field of ancient DNA (aDNA) is casting new light on many evolutionary questions. However, problems associated with the postmortem instability of DNA may complicate the interpretation of aDNA data. For example, in population genetic studies, the inclusion of damaged DNA may inflate estimates of diversity. In this paper, we examine the effect of DNA damage on population genetic estimates of ancestral population size. We simulate data using standard coalescent simulations that include postmortem damage and show that estimates of effective population sizes are inflated around, or right after, the sampling time of the ancestral DNA sequences. This bias leads to estimates of increasing, and then decreasing, population sizes, as observed in several recently published studies. We reanalyze a recently published data set of DNA sequences from the Bison (Bison bison/Bison priscus) and show that the signal for a change in effective population size in this data set vanishes once the effects of putative damage are removed. Our results suggest that population genetic analyses of aDNA sequences, which do not accurately account for damage, should be interpreted with great caution.

[1]  Alexander F. Auch,et al.  Metagenomics to Paleogenomics: Large-Scale Sequencing of Mammoth DNA , 2006, Science.

[2]  C. J-F,et al.  THE COALESCENT , 1980 .

[3]  E. G. Shpaer,et al.  Coalescent estimates of HIV-1 generation time in vivo. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[4]  C. Wiuf,et al.  Assessing the Fidelity of Ancient DNA Sequences Amplified From Nuclear Genes , 2006, Genetics.

[5]  James Haile,et al.  Ancient Biomolecules from Deep Ice Cores Reveal a Forested Southern Greenland , 2007, Science.

[6]  Beth Shapiro,et al.  Rise and Fall of the Beringian Steppe Bison , 2004, Science.

[7]  Edward M. Rubin,et al.  Genomic Sequencing of Pleistocene Cave Bears , 2005, Science.

[8]  A. Rambaut,et al.  BEAST: Bayesian evolutionary analysis by sampling trees , 2007, BMC Evolutionary Biology.

[9]  S. Pääbo.  Ancient DNA: extraction, characterization, molecular cloning, and enzymatic amplification. , 1989, Proceedings of the National Academy of Sciences of the United States of America.

[10]  H. Munro,et al.  Mammalian protein metabolism , 1964 .

[11]  E. Willerslev,et al.  Review Paper. Ancient DNA , 2005, Proceedings of the Royal Society B: Biological Sciences.

[12]  S. Ho,et al.  Evidence for time dependency of molecular rate estimates. , 2007, Systematic biology.

[13]  M. Beaumont,et al.  Novel high-resolution characterization of ancient DNA reveals C > U-type base modification events as the sole cause of post mortem miscoding lesions , 2007, Nucleic acids research.

[14]  J. Rothberg,et al.  Patterns of nucleotide misincorporations during enzymatic amplification and direct large-scale sequencing of ancient DNA , 2006, Proceedings of the National Academy of Sciences.

[15]  J. Kingman.  On the genealogy of large populations , 1982, Journal of Applied Probability.

[16]  A. von Haeseler,et al.  DNA sequences from multiple amplifications reveal artifacts induced by cytosine deamination in ancient DNA. , 2001, Nucleic acids research.

[17]  Feng Chen,et al.  Sequencing and Analysis of Neanderthal Genomic DNA , 2006, Science.

[18]  Adrian W. Briggs,et al.  Analysis of one million base pairs of Neanderthal DNA , 2006, Nature.

[19]  R. Nielsen,et al.  Ancient bacteria show evidence of DNA repair , 2007, Proceedings of the National Academy of Sciences.

[20]  Philip L. F. Johnson,et al.  Patterns of damage in genomic DNA sequences from a Neandertal , 2007, Proceedings of the National Academy of Sciences.

[21]  Stephan C. Schuster,et al.  Whole-Genome Shotgun Sequencing of Mitochondria from Ancient Hair Shafts , 2007, Science.

[22]  Alexei J Drummond,et al.  Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. , 2005, Molecular biology and evolution.

[23]  C. Wiuf,et al.  Statistical evidence for miscoding lesions in ancient DNA templates. , 2001, Molecular biology and evolution.

[24]  A. Wilson,et al.  DNA sequences from the quagga, an extinct member of the horse family , 1984, Nature.

[25]  S. Ho,et al.  Bayesian estimation of sequence damage in ancient DNA. , 2007, Molecular biology and evolution.

[26]  W. Miller,et al.  Recharacterization of ancient DNA miscoding lesions: insights in the era of sequencing-by-synthesis , 2006, Nucleic acids research.

[27]  E. Willerslev,et al.  Characterization of genetic miscoding lesions caused by postmortem damage. , 2003, American journal of human genetics.

[28]  S. Jeffery.  Evolution of Protein Molecules , 1979 .

[29]  J. W. Pendleton,et al.  Surveys of Gene Families Using Polymerase Chain Reaction: PCR Selection and PCR Drift , 1994 .

[30]  S. Pääbo,et al.  Lack of phylogeography in European mammals before the last glaciation. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[31]  O. Pybus,et al.  Bayesian coalescent inference of past population dynamics from molecular sequences. , 2005, Molecular biology and evolution.

[32]  T. Jukes.  CHAPTER 24 – Evolution of Protein Molecules , 1969 .

[33]  R. Hudson.  Properties of a neutral allele model with intragenic recombination. , 1983, Theoretical population biology.