zotero/storage/TDFCZDXH/.zotero-ft-cache

The half-life of DNA in bone: measuring decay kinetics in 158 dated fossils
Morten E. Allentoft1,2,3,*, Matthew Collins4, David Harker4,
James Haile1, Charlotte L. Oskam1, Marie L. Hale2, Paula
F. Campos3,5, Jose A. Samaniego3, M. Thomas P. Gilbert1,3,
Eske Willerslev3, Guojie Zhang6, R. Paul Scofield7,
Richard N. Holdaway2,8 and Michael Bunce1,*
1Ancient DNA Laboratory, School of Biological Sciences and Biotechnology, Murdoch University, 90 South Street, Perth, Western Australia 6150, Australia 2School of Biological Sciences, University of Canterbury, Private Bag 4800, Christchurch, New Zealand 3Centre for GeoGenetics, Natural History Museum of Denmark, University of Copenhagen, Øster Voldgade 5-7, 1350 Copenhagen K, Denmark 4Department of Archaeology, University of York, PO Box 373, York, UK 5Museu da Cieˆncia, University of Coimbra, Laboratorio Chimico, Largo Marqueˆs de Pombal, 3000-272 Coimbra, Portugal 6Beijing Genomics Institute-Shenzhen, Shenzhen 518083, People’s Republic of China 7Canterbury Museum, Rolleston Avenue, Christchurch 8050, New Zealand 8Palaecol Research Ltd, 167 Springs Road, Hornby, Christchurch 8042, New Zealand
Claims of extreme survival of DNA have emphasized the need for reliable models of DNA degradation through time. By analysing mitochondrial DNA (mtDNA) from 158 radiocarbon-dated bones of the extinct New Zealand moa, we confirm empirically a long-hypothesized exponential decay relationship. The average DNA half-life within this geographically constrained fossil assemblage was estimated to be 521 years for a 242 bp mtDNA sequence, corresponding to a per nucleotide fragmentation rate (k) of 5.50  10 – 6 per year. With an effective burial temperature of 13.18C, the rate is almost 400 times slower than predicted from published kinetic data of in vitro DNA depurination at pH 5. Although best described by an exponential model (R2 1⁄4 0.39), considerable sample-to-sample variance in DNA preservation could not be accounted for by geologic age. This variation likely derives from differences in taphonomy and bone diagenesis, which have confounded previous, less spatially constrained attempts to study DNA decay kinetics. Lastly, by calculating DNA fragmentation rates on Illumina HiSeq data, we show that nuclear DNA has degraded at least twice as fast as mtDNA. These results provide a baseline for predicting long-term DNA survival in bone.
Keywords: DNA degradation; aDNA; decay kinetics; DNA half-life
1. INTRODUCTION
Although early-1990s claims of DNA recovered from million year-old fossils [1– 4] are now widely regarded as modern contaminants [5– 13], the kinetics of longterm post-mortem DNA decay is still poorly understood. There is a lack of empirical data on which to estimate the rate of in situ DNA fragmentation. Because the field of ancient DNA (aDNA) has recently entered the era of whole-genome profiling [14 – 17], which is dependent on samples of exceptional preservation, understanding the nature and rate of DNA decay is as pertinent as ever, for both predictive and authentication purposes. DNA has limited chemical stability and decays without the enzymatic repair mechanisms of living cells [18].
Following cell death, nucleases start to cleave the DNA into fragments [19] and during decomposition, the DNA is digested by micro-organisms [18,20]. In determining long-term DNA decay, it is believed that hydrolysis of amino groups accelerates the loss of purine residues (depurination), resulting in strand cleavage [21,22]. This random DNA fragmentation generates a characteristic negative exponential correlation between DNA fragment length and number of molecules (figure 1a) [23,24,26,27]. Although other types of DNA damage occur, such as cross-links [28,29], the role of depurination in DNA fragmentation is well documented, and damage patterns from high-throughput sequencing data show that purines are overrepresented adjacent to strand breaks in ancient samples [30 – 33]. However, it is not known whether the rate of fragmentation can be regarded as constant through time, nor to what extent it varies between specimens from similar depositional environments. Indeed, a constant rate would imply that DNA decay follows first-order kinetics
* Authors for correspondence (morten.allentoft@gmail.com; m. bunce@murdoch.edu.au).
Electronic supplementary material is available at http://dx.doi.org/ 10.1098/rspb.2012.1745 or via http://rspb.royalsocietypublishing.org.
Proc. R. Soc. B (2012) 279, 4724–4733 doi:10.1098/rspb.2012.1745 Published online 10 October 2012
Received 26 July 2012 Accepted 14 September 2012 4724 This journal is q 2012 The Royal Society


(i.e. DNA has a half-life; figure 1b), and hence that the age of a specimen can potentially work as a proxy for its DNA preservation and vice versa. Despite many attempts [24,31,34 –41], it has proved extremely difficult to demonstrate a general association between age and DNA preservation, probably because of variation in physical, chemical and biological factors such as taphonomy, fossil storage, oxygenation, microbial diagenesis, pH and ionic strength, and the presence of cations, humics and humates. Moreover, most studies are limited by small sample sizes and lack of individually dated specimens, leading to potential problems of mixed taphonomies. The absence of clear temporal trends in previous work suggests that the rate of DNA decay either fluctuates widely (hence not a rate per se), or that it takes a large, homogenous, and accurately dated sample, to overcome the ‘noise’ introduced by the aforementioned factors. In the light of many inconclusive results, a common refrain is that DNA decay rates cannot be predicted [42 – 44]. A few studies, however, have indicated a relationship between DNA preservation and post-mortem age. Bar et al. [45] showed that the length of amplifiable DNA fragments from human tissues (blood, muscle and spleen) decreased exponentially during the 20 days after a person’s death, and Campos et al. [46] showed a rapid decrease in DNA content in cow bone after burial. Moreover, Adler et al. [26] showed a weak exponential correlation between DNA preservation and age, on the basis of 37 ancient human samples (mainly teeth). These samples had not all been individually dated, but the result is encouraging and provides impetus for more detailed investigations of long-term DNA fragmentation. Establishing an association between age and preservation is not the only challenge. The rate of depurination is influenced by temperature, among other factors [22], which explains why the most extreme survival of DNA was documented in approximately 450–800 kyr ice cores [47]. Smith et al. [11,12] argued for the possibility of using the temperature-dependence of DNA fragmentation to normalize the samples and predict DNA survival. Such a relationship has been demonstrated for collagen, the most abundant protein in bone [48]. However, to establish
a thermal age model for DNA, the first step is to confirm that long-term in situ DNA degradation can be described by a rate kinetic. In an attempt to document a correlation between sample age and DNA preservation, we use a quantitative real-time PCR (qPCR) design to measure relative copy numbers of mitochondrial DNA (mtDNA) fragments from bones of the extinct New Zealand moa (Aves: Dinornithiformes). Aside from a large dataset (qPCR data from 158 radiocarbon-dated bones), our study differs from previous attempts by being highly localized. Using Holocene fossil tibiotarsi recovered from anoxic, limestone-buffered sediments at three adjacent (less than 5 km apart) fossil deposits (figure 2), we endeavour to provide a high level of sample homogeneity and minimize variables of climate and taphonomy. Also, by sampling only natural bone accumulations, we expect to minimize other complicating factors (e.g. butchering and cooking) that could impact on taphonomic processes [27,49]. Because the three fossil sites were excavated at different times (figure 2), the influence of storage times on DNA recovery can also be assessed. Pruvost et al. [50] argued that DNA degradation intensifies when a bone is removed from its deposition environment. The moa fossils analysed here were recovered at well-documented dates over an approximately 70 year period, allowing us to test the importance of post-excavation storage time. The three main objectives of this study were: (i) to test whether long-term DNA decay follows first-order kinetics, thereby confirming the foundation for a predictive model; (ii) to estimate the long-term decay rate in bone at a given burial temperature and compare this rate with the predicted depurination rate from DNA in solution [21,22]; and (iii) to estimate the relative importance of storage time on DNA preservation in bone. Lastly, following Deagle et al.’s [23] demonstration that frequency distributions of DNA fragment lengths from a sample of known age can be used to estimate the DNA decay rate, we argue that histograms based on the thousands-to-millions of reads generated from highthroughput sequencing platforms could represent a ‘by-product’ to investigate DNA fragmentation. We therefore apply the method to Illumina HiSeq 2000 data
10
20
30
40
50
60
70
80
90
100
time (years after cell death)
0
0 00
10
20
30
40
50
60
70
80
90
100
(a) (b)
copy number (relative)
fragment length (bp)
50 100 150 200 250 300
f(T) = f(T0)*e–k50*T
damage rate (k) = 0.02 per site per year
50 100 150 200 250 300
f(L) = f(L0)*e–l*L damage fraction (l) = 0.02 per site
50 bp copy number
(relative)
Figure 1. DNA fragmentation theory. (a) The exponential relationship caused by random fragmentation of DNA. Post-mortem, the template fragment length (L) distribution follows an exponential decline determined by the proportion of damaged sites (l). This relationship has been described from both modern and ancient samples [23– 25]. Here, a fragment size distribution representing l 1⁄4 0.02 (2% of the bonds in the DNA backbone are broken). (b) A hypothetical signal of temporal DNA decay, which has, prior to this study, been extremely difficult to demonstrate. The model assumes that the observed damage fraction (l) can be converted to a rate of decay (k) when the age (T ) of a sample is known. It implies that the number of DNA copies of a given length (L) will decline exponentially with time—hence the notion that DNA has a half-life. Here, the theoretical decay kinetics of a 50 bp DNA fragment, assuming a k of 2% per site per year. k is converted to a 50 bp decay rate (k50), according to
a Poisson distribution as: k50 1⁄4 1 – (e20.02*50).
The half-life of DNA M. E. Allentoft et al. 4725
Proc. R. Soc. B (2012)


generated from two moa specimens, in order to achieve independent DNA decay-rate estimates, comparable with the qPCR results.
2. MATERIAL AND METHODS
(a) Sampling, dating and storage
Sampling of 158 moa left tibiotarsi was conducted as described previously [51]. The material covered three adjacent sites, located within a 5 km radius in North Canterbury, South Island, New Zealand; Pyramid Valley (PV; n 1⁄4 103), Bell Hill Vineyard (BHV; n 1⁄4 47) and Rosslea (n 1⁄4 8) (figure 2). Further information on samples and sites has been given previously [52,53] and is also listed in electronic supplementary material, table S1. The PV and BHV specimens were sampled from museum collections, whereas the Rosslea material was sampled less than 14 days after excavation in 2008 [53]. Bone gelatin was extracted by Isolytix Ltd, Dunedin, New Zealand, using standard protocols. Gelatin sample splits were sent for stable isotopic analysis (Environmental Isotopes Ltd, Sydney, Australia) and accelerator mass spectrometry (AMS) radiocarbon analysis. Radiocarbon ages were measured at the Rafter Radiocarbon Laboratory, GNS Science. Lower Hutt, New Zealand. Ages are reported as conventional radiocarbon ages, and as median calibrated ages and 95% probability distributions generated by the OXCAL v. 4.1 (Oxford Radiocarbon Accelerator Unit) program. Sample quality was assessed using the C : N molar mass ratio, and stable isotopic ratios (d15N, d13C) in comparison with archive data, and was uniformly high (see the electronic supplementary material, table S1). The moa bones used in this study had been excavated at different times (figure 2). The PV elements were recovered during excavations from 1939 to 1973 [54], BHV was excavated in 2001 [52] and Rosslea in 2008 [53]. To assess the
potential DNA degradation effect of storage time, a multiple regression analysis (SPSS v. 19) was conducted between DNA preservation, calibrated 14C age and storage time. This analysis was conducted on the entire dataset (n 1⁄4 158) and exclusively on 103 elements from PV that were excavated over a 35-year time frame, with excavation dates recorded for each bone.
(b) DNA extraction and qPCR
DNA was isolated from all samples at dedicated aDNA facilities (Murdoch University, Perth, Australia), using a standardized method described previously [52] and detailed in the electronic supplementary material, text S1. All DNA extractions were undertaken from 200 mg uniformly sampled and powdered bone. The DNA isolation protocol (molecular weight cut-off (MWCO) and Silica purification) was selected owing to its ability to generate DNA largely free of inhibitors. DNA amplifications of the 262F/441R fragment (mtDNA control region, 242 bp including primers—see Bunce et al. [55]) were performed using SYBR-green detection chemistry on a BioRad My-IQ thermocycler. With a standardized level of detection, the recorded cycle threshold (CT) values represented the relative DNA preservation in each extract. To minimize effects of post-extraction storage time, qPCR was conducted immediately following DNA extraction. The quantitative response and general applicability of our qPCR set-up has been validated before based partly on the same extracts [51,56–58]. In addition, we conducted melt curve analyses for all PCR products, (see the electronic supplementary material, text S1 and figure S1), as well as ‘spike-in’ (see the electronic supplementary material, table S2), dilution(see the electronic supplementary material, figure S3 and table S3) and PCR efficiency (see the electronic supplementary material, figure S3) experiments. This was done to ensure that amplification of the correct target was being
10 km
Rosslea Bell Hill Vineyard
Pyramid Valley
Pyramid Valley (n = 103): shallow seasonal lake excavated: 1939–1973 age: 602–4136 BP (m = 1970 BP) mean CT: 34.1 (28.5–43.8)
excavated: 2001 age: 851–2876 BP (m = 1923 BP)
prehistoric stream
Bell Hill Vineyard (n = 47):
mean CT: 31.7 (24.8–40.4)
excavated: 2008
spring hole
Rosslea (n = 8):
age: 1482–7839 BP (m = 5280 BP) mean CT: 39.6 (30.5–46.9)
Figure 2. Study site. The three fossil deposits, PV (42858022.000 S, 172835049.000 E), BHV (42858019.3600 S, 172839056.1500 E) and Rosslea (42857053.8300 S, 172839022.3900 E), from which 158 radiocarbon-dated moa fossils were characterized for DNA decay kinetics. The sites in North Canterbury, South Island, New Zealand are located within a 5 km radius in the eastern rain shadows of the Southern Alps. Most of the area is more than 200 m a.s.l. and consists of flat alluvial plains and rolling downlands. Information on calibrated radiocarbon ages and DNA preservation (CT values) are shown for each site. m, mean age.
4726 M. E. Allentoft et al. The half-life of DNA
Proc. R. Soc. B (2012)


monitored and to test for PCR inhibition (see the electronic supplementary material, text S1). In most instances (greater than 70%), the actual qPCR products were directly sequenced, confirming their species origin (data not shown) and specificity of the amplicons. Extracts displaying deviating melt curve diagnostics (owing to non-specific amplifications and/or dimer accumulation) were omitted from the analysis (ca 20 samples). We acknowledge that the relatively large (242 bp) amplicons, with primers that had to work on multiple species and across a large dynamic range of DNA yields, is not ideal in a qPCR set-up. All of our tests indicated, however, that the assay is sufficiently robust and reproducable for the intended purpose. Stochastic variation between qPCR reactions is expected, in particular when amplifying low copy number DNA at CT values greater than 35. Results from 50 extracts in repeated qPCRs (same extract, but independent PCR mastermix and thermocycling run) were used to isolate this effect and estimate the variation introduced by stochastic qPCR effects. To test for consistency across DNA extracts, one moa bone was sampled and extracted three times, showing only minor deviations in CT (less than 1 cycle) between extracts at multiple levels of dilution (see the electronic supplementary material, table S4). Preliminary data suggested that DNA from one of the four North Canterbury moa species (Pachyornis elephantopus) amplified less efficiently in the qPCR assay. We suspect that primer-binding-site mutations are the cause, and to prevent a systematic bias, data from this species were excluded.
(c) Estimating the rate of decay from qPCR data
The number of DNA copies in a PCR reaction follows an exponential increase per cycle. If one DNA extract produces a CT three cycles higher than another, then the former has eight times (23) less DNA (in an ideal situation with 100% amplification efficiency). To test for an association between DNA preservation and age, we calculated the relative number of template molecules in each sample according to the ‘best’ sample (Dinornis robustus, S39955, CT 1⁄4 24.8), which was set to 100 per cent. Although some previous aDNA studies convert CT values into absolute starting copy numbers [24,59], others do not [46] as the accuracy and meaningfulness can sometimes be questioned [58]. For our purpose, the use of relative CT values was sufficient, as we were interested only in rates and not in absolute numbers. Following Adler et al. [26], we tested (in SPSS v. 19) whether the observed DNA degradation through time would be better described by a linear, inverse, S-shaped or exponential decay model. An exponential decay relation could be estimated as: Nt 1⁄4 N0  ek242t (Nt and N0 being the quantity at time t and 0, respectively, and k242 the decay constant for the entire 242 bp fragment). The average molecular half-life (t1/2 1⁄4 (ln 2)/k242) could be calculated from the entire dataset and for each fossil site, respectively. A model of random fragmentation predicts that the relationship between decay rate and fragment length is a Poisson distribution (figure 1) [23]. The average per nucleotide rate of decay is therefore calculated by isolating k in: k242 1⁄4 1 – (e2k*242).
(d) Estimating the rate of decay from Illumina HiSeq data
DNA extracts of two South Island giant moa (D. robustus) fossils were selected for Illumina HiSeq 2000 sequencing (see the electronic supplementary material, text S1 for details). Only a fraction of ‘shot-gun’ sequencing data from an ancient bone
is expected to represent DNA of the target species. Microbial DNA will constitute a proportion of the reads [60,61]. To retrieve sequences of endogenous moa DNA, only reads that mapped against an unpublished draft of the ostrich (Struthio camelus) genome (source tissue provided by San Diego Zoo and advance access to genomic data provided courtesy of Oliver Ryder, San Diego Zoo; Genome 10K; BGI, China) and the mtDNA genome of a D. robustus moa (GenBank: AY016013.1) [62] were used in the downstream analyses. The software BWA [63] was used in the mapping, using default settings. The resulting fragment length histograms of all mapped reads were used to assess DNA decay rates. According to a model of random fragmentation, the amount of amplifiable template should decline exponentially with increasing fragment size [23] (figure 1a). Log-transformed copy numbers therefore have a linear relationship with amplicon length, with the slope of the decline (l) describing the probability of a bond in the DNA backbone being cleaved [23]. Lambda can be converted to a damage rate when sample age is known (figure 1b).
3. RESULTS AND DISCUSSION
(a) Estimating the half-life of mtDNA in bone
The calibrated radiocarbon ages of the 158 fossils documented a Mid- to Late Holocene accumulation, ranging from 602 to 7839 BP, calibrated calendar years before present (‘present’ defined as CE 1950; figure 2; electronic supplementary material, table S1). When correlating relative DNA copy number with age, a linear (R2 1⁄4 0.002, p 1⁄4 0.604), inverse (R2 1⁄4 0.002, p 1⁄4 0.598) and S-shaped (R2 1⁄4 0.229, p , 0.001) model proved inferior to an exponential relationship (R2 1⁄4 0.386, p , 0.001; figure 3). Hence, the qPCR results demonstrated that the number of template mtDNA molecules in the fossils was best described with first-order decay kinetics (see the electronic supplementary material, table S4 for individual CT values). The equation for the fitted regression was:
Nt 1⁄4 3:61  e0:0013t; ð3:1Þ
with Nt being the proportion of molecules left at the time t (years BP), compared with the best sample. Solving the equation yielded a molecular half-life of 521 years for the targeted 242 bp mtDNA fragment, and assuming a Poisson distribution [23], the per site decay rate (k) was 5.50  10 – 6 per year. As noted, several studies have attempted to investigate DNA degradation kinetics using samples of different ages. Here, we show that by using a controlled experimental set-up, with samples normalized for ambient temperature, an exponential relationship between sample age and DNA preservation can be demonstrated. This supports the view that random DNA fragmentation, likely caused by depurination, applies not only to short-term degradation of DNA in a solution [21,22] but also to ancient DNA in fossil bone.
(b) Comparison with Lindahl & Nyberg’s [22] results
The rate constant (k) depends on absolute temperature (T ) according to the relationship:
ln k 1⁄4 ln(A)  Ea
R  (1=T ) ; ð3:2Þ
The half-life of DNA M. E. Allentoft et al. 4727
Proc. R. Soc. B (2012)


where A is the pre-exponential factor, Ea is the activation energy and R the gas constant. On the basis of the Arrhenius plot in Lindahl & Nyberg [22], showing the temperature-dependency of depurination at pH 5, we estimated k to 2.11  10–3 per site per year (see the electronic supplementary material, text S2), which is 384 times faster than the rate estimated from the moa data (figure 4). Assuming that long-term DNA fragmentation happens primarily as a result of depurination, this discrepancy is probably related to pH and access to free water. Even if the burial environment is slightly acidic, sacrificial dissolution of microcrystalline carbonated bioapatite will likely act as buffer within the bone. On the basis of the measurements in Berna et al. [64], we used pH 7.5 as an appropriate value for bone (see the electronic supplementary material, Text S2). From Lindahl & Nyberg [22], we find that k at a pH of 7.5 is ca 73 times slower than the same at pH 5 (see the electronic supplementary material, Text S2), but still ca five times faster than the rate observed for moa mtDNA (figure 4). Given that the same authors also reported a twofold rate decrease in the presence of apatite, as in bone (cited in [18]), the rate expected from Lindahl & Nyberg’s [22] in vitro results is perhaps only roughly twice as fast, as the rate we observe for the moa qPCR estimate. From equation (2.2), it is clear that both the preexponential factor (A) and the activation energy (Ea) affect the decay rate (k), but Lindahl & Nyberg [22] showed similar values of Ea at pH 5 and 6. We therefore assume that the estimated Ea of 127 kJ mol21 also applies to pH 7.5. With this Ea value, and an observed average k for moa mtDNA of 5.5  10 –6 per site per year, and an estimated fossil burial temperature T 1⁄4 13.18C or
286.25 K (see the electronic supplementary material, text S3), the pre-exponential factor A is determined to e41.2. The average decay rate (per site per year) for mtDNA in moa bone is therefore related to temperature as follows (figure 4):
ln k 1⁄4 41:2  15267:6  1
T : ð3:3Þ
This is arguably the best available approximation for the temperature-dependency of long-term decay of mtDNA in bone, and it should be used as a guide for future work on DNA decay kinetics. We note that Lindahl & Nyberg [22] tested only the temperature dependency of depurination between 458C and 808C (318.15 –353.15 K), and it remains to be confirmed that the relationship in their Arrhenius plot can be projected to include more typical burial temperatures, as seen in table 1. Furthermore, we note that our model of mtDNA decay in bone does not account for a potential initial post-mortem phase of faster DNA decay, facilitated by nuclease activity rather than depurination [23,31,65]. With these caveats in place, table 1 shows our predicted half-lives, and average fragment lengths of mtDNA in bone at various temperatures.
(c) Variance in DNA preservation
Only 38.6 per cent (figure 3) of the variation in DNA preservation could be explained by the age of the fossils. Despite efforts to minimize variation (see §2), part of this variance can be ascribed to qPCR stochasticity. Results from qPCR repeats of 50 DNA extractions showed an average of 0.9 cycle shift in CT value (see the electronic supplementary material, table S5). In 36 of the 50 repeats, the two qPCRs deviated with less than one cycle, confirming the general consistency in the set
(a) 100
(b)
all sites combined, n = 158
R2 = 0.386
R2 = 0.359 R2 = 0.058
PV, n = 103 BHV, n = 47 Rosslea, n = 8
R2 = 0.752
10
0 2000 4000 years BP
relative copy number
6000
0 5000
8000 10 000
0 5000 10 000 10 000 0 5000 10 000
1
0.1
0.01
0.001
0.0001
0.00001
0.00001
0.0001
0.001
0.01
0.11
10
100
Figure 3. Correlations between age and DNA preservation. Relative mtDNA copy numbers (determined by qPCR) in moa bone plotted against age for all 158 fossils (a), and for each of the three deposits respectively (b). The exponential correlations are significant (p , 0.005) except for the BHV data (p 1⁄4 0.1) when tested alone. Although a faster average decay is observed at Rosslea, the decay rates (slopes) did not differ significantly from each other.
4728 M. E. Allentoft et al. The half-life of DNA
Proc. R. Soc. B (2012)


time in years
time in years
1.0
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
2000 4000 6000 8000 10 000 12 000 14 000 16 000 18 000 20 000
mtDNA qPCR (k = 5.50 × 10–6)
predicted pH 5 (k = 2.11 × 10–3)
predicted pH 7.5 (k = 2.88 × 10–5)
mtDNA HiSeq 1 (k = 8.64 × 10–6)
mtDNA HiSeq 2 (k = 1.35 × 10–5)
nuDNA HiSeq 1 (k = 2.13 × 10–5)
nuDNA HiSeq 2 (k = 2.71 × 10–5)
(b) 1.0
(a)
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
500 1000 1500 2000 2500 3000 3500 4000
0
average mtDNA half-life measured to 521 years
242 bp survival bond survival
Figure 4. Observed and predicted rates of DNA decay. The predicted survival of DNA in bone through time, measured as intact phosphodiester bonds in the DNA backbone (a), and survival of a 242 bp fragment (b). The depicted survival rates are based on: (i) the average mtDNA decay rate measured directly from qPCR of 158 moa bones; (ii) the rate of depurination measured from DNA in solution at pH 5 in Lindahl & Nyberg [22] but adjusted to 13.18C to allow comparison with the moa data; (iii) the same rate adjusted further to pH 7.5, as expected inside a bone; (iv) mtDNA and nuDNA decay rates calculated based on Illumina HiSeq data from two moa samples (HiSeq 1 from sample S40114, HiSeq 2 from sample S39946-3). The estimated decay rate (k, per site per year) is listed for each of the seven datasets.
Table 1. Predictions of decay rates (k) of mtDNA in bone at various temperatures (based on equation (3.3)). Estimates of mtDNA half-lives for three fragment lengths are indicated as well as the expected average read length (1/l, where l is damage fraction) after 10 000 years. The decay rates do not account for the potential initial post-mortem phase of rapid DNA decay governed by nucleases. Still, the results indicate that under the right conditions of preservation, short fragments of DNA should be retrievable from very old bone (e.g. greater than 1 Myr). However, even under the best preservation conditions at 258C, our model predicts that no intact bonds (average length 1⁄4 1 bp) will remain in the DNA ‘strand’ after 6.8 Myr. This displays the extreme improbability of being able to amplify a 174 bp DNA fragment from an 80–85 Myr old Cretaceous bone [1].
temperature
k per site per year
half-life (years), 30 bp
half-life (years), 100 bp
half-life (years), 500 bp
average length at 10 kyr
time (years) until average length 1⁄4 1 bp
258C 4.5  10 –5 500 150 30 2 bp 22 000 158C 7.6  10 –6 3000 900 180 13 bp 131 000 58C 1.1  10 –6 20 000 6000 1200 88 bp 882 000 258C 1.5  10 –7 158 000 47 000 9500 683 bp 6 830 000
The half-life of DNA M. E. Allentoft et al. 4729
Proc. R. Soc. B (2012)


up. The 50 qPCR repeats spanned a total of 14.8 cycles in CT, implying that roughly 6 per cent (0.9/14.8) of the total variation in CT could be ascribed to qPCR stochasticity. This value is the procedurally introduced background noise in our estimated decay rate. Lastly, with regard to accuracy, it is possible that low-template aDNA extracts (typically with CT . 35) could include an amplification signal of template fragments forming as a result of incomplete extension (and where fragments eventually forming a full length product in latter PCR cycles). Such fragmentary chimaeric templates forming as a qPCR artefact might subtly impact our estimates of decay kinetics on the most poorly preserved bones. The mean CT value for the Rosslea samples (39.6) was significantly higher than those from BHV (31.7; t-test, unequal variances, td.f. 1⁄4 8 1⁄4 4, p 1⁄4 0.004) and PV (34.1; t-test, unequal variances, td.f. 1⁄4 7 1⁄4 2.9, p 1⁄4 0.022). When converted to relative DNA copy numbers, PV fossils contained on average 19.6 per cent of the template molecules present in BHV fossils, and Rosslea contained just 0.3 per cent, corresponding to approximately 330 times fewer DNA templates. A ‘spike-in’ experiment showed that the poor DNA preservation at Rosslea did not result from PCR inhibition (see the electronic supplementary material, table S2). Whereas poorer DNA preservation at Rosslea may be explained by its greater age (figure 2), the difference between the largely contemporaneous PV and BHV sites is more intriguing. This difference was significant, not only when including data from the entire combined PV-BHV time frame (t-test, equal variances, td.f. 1⁄4 148 1⁄4 4.5, p , 0.001), but also within the 1882 years of temporal overlap (915–2797 BP, n 1⁄4 72). Eliminating age as a contributing factor, PV fossils still displayed significantly higher CT (i.e. poorer DNA preservation) than those from BHV (t-test, unequal variances, td.f. 1⁄4 72 1⁄4 3.8, p , 0.001). We tested whether these site-specific differences in DNA preservation could be explained by differing decay rates. After log-transformation, PV (R2 1⁄4 0.359, ANOVAd.f. 1⁄4 102: F 1⁄4 56.53, p , 0.001) and Rosslea (R2 1⁄4 0.749, ANOVAd.f. 1⁄4 7: F 1⁄4 17.91, p 1⁄4 0.005) both displayed highly significant regressions, whereas the same decay signal was not supported for BHV (R2 1⁄4 0.058, ANOVAd.f. 1⁄4 46, F 1⁄4 2.79, p 1⁄4 0.10; figure 3). The average molecular half-lives for DNA in bones from PV and Rosslea were 506 years and 389 years, respectively. However, an analysis of covariance (based on standard error of difference between slopes) showed that the decay rate was not significantly higher at Rosslea than in PV (td.f. 1⁄4 107 1⁄4 0.92, p 1⁄4 0.36) or the whole dataset (td.f. 1⁄4 162 1⁄4 0.70, p 1⁄4 0.48). It was therefore justifiable to pool the data from the three sites to achieve an average DNA decay rate for DNA in a North Canterbury bone (figure 3). We could not detect an exponential decay process in the BHV material; so the decay rate was not estimated for this site alone. Notably, whereas 74.9 per cent of the variance at Rosslea was explicable by geological age, only 5.8 per cent was explicable by age at BHV (figure 3). The discrepancy probably reflects differences in the periods of deposition at the sites (shorter at BHV), allowing geological age to dominate the relationship at Rosslea but not at BHV. We then examined the effect of storage time. The difference in average CT between PV and BHV could be an effect of the seasonal water saturation at the PV site, but it
could also be a storage effect, given that the material was collected from BHV in 2001, whereas most of the PV material had been at Canterbury Museum for more than 50 years (figure 2). However, we note that the DNA in bones from Rosslea was more degraded, despite the bones having been sampled shortly after excavation. Moreover, the range in CT (figure 2) was highest at BHV although the 47 bones from this site had identical storage times (7 years between excavation and sampling). A previous investigation showed that storage conditions can be a major contributor to variation in DNA preservation [50]. Our results nuance this observation by showing that storage time in itself does not seem to have a large effect. A multiple regression analysis showed that both variables (age and storage) had a significant effect (p , 0.001 for both), but whereas in combination, they explained 45.2 per cent of the variation in DNA preservation, only 7.7 per cent of this could be ascribed solely to storage time. Therefore, in our experimental setup, the geological age of a sample was a far better predictor of its DNA preservation than storage time.
(d) Decay rates from Illumina data: mtDNA versus nuDNA
Illumina sequencing of two DNA extracts yielded 10.1 and 24.8 million reads, respectively. When filtered and mapped against a draft of the ostrich genome (see §2), the return was 3128 and 34 495 unique, non-duplicated reads (see the electronic supplementary material, table S6). These should represent endogenous fragmented moa DNA and hence suitable for calculations of DNA decay rates. It was first confirmed that the number of fragments declined exponentially with fragment lengths (see the electronic supplementary material, table S6 and figure S3), thereby conforming to the random fragmentation model [23]. The DNA decay rate (k) was calculated as 2.13  10–5 and 2.71  10–5 per site per year, for each sample, respectively (figure 4). The majority of DNA within most vertebrate cells is derived from the nucleus, thus it was a priori hypothesized that the majority of the mapped ‘shotgun’ reads would represent nuclear DNA fragments. If there is a difference between post-mortem DNA decay rates in the nuclear and mtDNA genomes, it could explain why the rate is four to five times faster than estimated from our qPCR data (5.5  10–6 per site per year), which represent fragmentation of mtDNA. To investigate this further, the shotgun sequences were then mapped against the moa (D. robustus) mitochondrial genome, returning 878 and 4032 reads, respectively. When the decay rate was calculated based only on these mtDNA sequences, it was found to be 2 and 2.5 times slower, respectively, than that estimated from the entire dataset (figure 4, electronic supplementary material, table S6). These numbers suggest that mtDNA degrades at a slower rate than nuDNA, consistent with observations reported in the study of Schwarz et al. [24]. This could perhaps be explained by the circular structure of mtDNA, making it less accessible to exonuclease activity. The two decay rates of mtDNA measured from Illumina data appear slightly faster than the rate based on qPCR (figure 4). This could be sampling error because these are just two individual data points compared with the qPCR average of 158 samples. However, it could also be explained by the fact that the qPCR decay rate
4730 M. E. Allentoft et al. The half-life of DNA
Proc. R. Soc. B (2012)


is based on differences in DNA preservation between bones that have all been in the ground for more than 600 years. This rate does not therefore incorporate the presumed faster DNA degradation occurring in the initial post-mortem phase. It has also been demonstrated that the library preparation for high-throughout sequencing may skew the fragment length distribution towards shorter average reads [66], and such bias could potentially explain the observed rate difference between our qPCR and HiSeq data. However, it is not clear whether such bias would affect only the average read length, or whether it could also alter the slope of the frequency distribution, which we use to estimate k. Further data are required to clarify these potential biases.
4. CONCLUDING REMARKS
We have demonstrated that in situ DNA decay is described by first-order kinetics, confirming that long-term postmortem DNA fragmentation can be treated as a rate process. This closes the gap to theoretical in vitro observations made four decades earlier. We argue that equation (3.3) represents the best available approximation of the rate of mtDNA decay in fossil bone. Empirical data from a range of different depositional environments, also extending back into the Pleistocene, are now needed to further test and refine this model. Importantly, we show how high-throughput sequencing data can be used to estimate DNA decay rates, and thus the vast amount of such data currently being generated worldwide may provide an excellent resource by which to study DNA decay. Factors such as bone thickness, burial depth, surrounding pH and water saturation can be tested and modelled in future research, whereas information on season of death, speed of cadaver incorporation into sediment and climatic fluctuations will probably be less accessible. It is tempting to suggest that we can now predict the temporal limits of DNA survival, and finally refute the claims of authentic DNA from Cretaceous and Miocene specimens. This is, however, not straightforward. One needs information on the number of template molecules in living tissues, and estimates of post-mortem DNA decay rates for each tissue type. However, the half-life predictions (table 1) display the extreme improbability that an authentic 174 bp long mtDNA fragment of an 80 – 85 Myr old bone could have been amplified [1]. Our results indicate that short fragments of DNA could be present for a very long time; at –58C, the model predicts a half-life of 158 000 years for a 30 bp mtDNA fragment in bone (table 1). Even rough estimates such as this imply that sequenceable bone DNA fragments may still be present more than 1 Myr after deposition in deep frozen environments. It therefore seems reasonable to suggest that future research may identify authentic DNA that is significantly older than the current record of approximately 450–800 kyr from Greenlandic ice cores [47].
We are grateful to Malene Møhl for help with the sampling and to Emma McLay, Helen Hunt and Jayne Houston for technical assistance, and to Søren Overballe-Petersen for useful discussions. We thank the Hodgen (PV), Giesen (BHV) and Earl (Rosslea) families for their support of our research. Thanks to Oliver Ryder, San Diego Zoo for providing ostrich tissue. Financial support was provided by
the Marsden Fund of the Royal Society of New Zealand (contract 06-PAL-001-EEB to Palaecol Research Ltd). M.C. and D.H. received support from the SYNTHESYS project (http://www.synthesys.info/), financed by European Community Research Infrastructure Action under the FP7 ‘Capacities’ programme. M.T.P.G. and J.A.S. were supported by the Danish Council for Independent Research (Grant 10-081390). M.B. is funded by the Australian Research Council (Future Fellowship FT0991741).
REFERENCES
1 Woodward, S. R., Weyand, N. J. & Bunnell, M. 1994 DNA-sequence from Cretaceous period bone fragments. Science 266, 1229 –1232. (doi:10.1126/science.7973705) 2 Golenberg, E. M., Giannasi, D. E., Clegg, M. T., Smiley, C. J., Durbin, M., Henderson, D. & Zurawski, G. 1990 Chloroplast DNA-sequence from a Miocene magnolia species. Nature 344, 656–658. (doi:10.1038/344656a0) 3 Desalle, R., Gatesy, J., Wheeler, W. & Grimaldi, D. 1992 DNA-sequences from a fossil termite in Oligomiocene amber and their phylogenetic implications. Science 257, 1933–1936. (doi:10.1126/science.1411508) 4 Cano, R. J., Poinar, H. N., Pieniazek, N. J., Acra, A. & Poinar, G. O. 1993 Amplification and sequencing of DNA from a 120– 135-million-year-old weevil. Nature 363, 536 –538. (doi:10.1038/363536a0) 5 Sutlovic, D., Gamulin, S., Definis-Gojanovic, M., Gugic, D. & Andjelinovic, S. 2008 Interaction of humic acids with human DNA: proposed mechanisms and kinetics. Electrophoresis 29, 1467–1472. (doi:10.1002/Elps. 200700699) 6 Marota, I., Basile, C., Ubaldi, M. & Rollo, F. 2002 DNA decay rate in papyri and human remains from Egyptian archaeological sites. Am. J. Phys. Anthropol. 117, 310 –318. (doi:10.1002/ajpa.10045) 7 Binladen, J., Gilbert, M. T. P. & Willerslev, E. 2007 800 000 year old mammoth DNA, modern elephant DNA or PCR artefact? Biol. Lett. 3, 55–56. (doi:10. 1098/rsbl.2006.0537) 8 Pa ̈a ̈bo, S. & Wilson, A. C. 1991 Miocene DNA sequences: a dream come true? Curr. Biol. 1, 45–46. (doi:10.1016/0960-9822(91)90125-G) 9 Willerslev, E., Hansen, A. J. & Poinar, H. N. 2004 Isolation of nucleic acids and cultures from fossil ice and permafrost. Trends. Ecol. Evol. 19, 141 –147. (doi:10. 1016/J.Tree.2003.11.010) 10 Willerslev, E., Hansen, A. J., Ronn, R., Brand, T. B., Barnes, I., Wiuf, C., Gilichinsky, D., Mitchell, D. & Cooper, A. 2004 Long-term persistence of bacterial DNA. Curr. Biol. 14, R9 –R10. (doi:10.1016/j.cub. 2003.12.012) 11 Smith, C. I., Chamberlain, A. T., Riley, M. S., Stringer, C. & Collins, M. J. 2003 The thermal history of human fossils and the likelihood of successful DNA amplification. J. Hum. Evol. 45, 203–217. (doi:10.1016/s00472484(03)00106-4) 12 Smith, C. I., Chamberlain, A. T., Riley, M. S., Cooper, A., Stringer, C. B. & Collins, M. J. 2001 Neanderthal DNA: not just old but old and cold? Nature 410, 771 772. (doi:10.1038/35071177) 13 Gilbert, M. T. P., Barnes, I., Collins, M. J., Smith, C., Eklund, J., Goudsmit, J., Poinar, H. & Cooper, A. 2005 Long-term survival of ancient DNA in Egypt: response to Zink and Nerlich (2003). Am. J. Phys. Anthropol. 128, 110 –114. (doi:10.1002/Ajpa.20045) 14 Rasmussen, M. et al. 2011 An aboriginal Australian genome reveals separate human dispersals into Asia. Science 334, 94–98. (doi:10.1126/science.1211177)
The half-life of DNA M. E. Allentoft et al. 4731
Proc. R. Soc. B (2012)


15 Rasmussen, M. et al. 2010 Ancient human genome sequence of an extinct Palaeo-Eskimo. Nature 463, 757 –762. (doi:10.1038/nature08835) 16 Green, R. E. et al. 2010 A draft sequence of the Neandertal genome. Science 328, 710 –722. (doi:10.1126/science. 1188021) 17 Keller, A. et al. 2012 New insights into the Tyrolean Iceman’s origin and phenotype as inferred by whole-genome sequencing. Nat. Commun. 3, 698. (doi:10.1038/ ncomms1701) 18 Lindahl, T. 1993 Instability and decay of the primary structure of DNA. Nature 362, 709 –715. (doi:10.1038/ 362709a0) 19 Darzynkiewicz, Z., Juan, G., Li, X., Gorczyca, W., Murakami, T. & Traganos, F. 1997 Cytometry in cell necrobiology: analysis of apoptosis and accidental cell death (necrosis). Cytometry 27, 1–20. (doi:10.1002/ (SICI)1097-0320(19970101)27:1,1::AID-CYTO2.3. 0.CO;2-L) 20 Eglinton, G., Logan, G. A., Ambler, R. P., Boon, J. J. & Perizonius, W. R. K. 1991 Molecular preservation. Phil. Trans. R. Soc. Lond. B 333, 315 –328. (doi:10.1098/ rstb.1991.0081) 21 Lindahl, T. & Andersson, A. 1972 Rate of chain breakage at apurinic sites in double-stranded deoxyribonucleic acid. Biochemistry 11, 3618 –3623. (doi:10.1021/ bi00769a019) 22 Lindahl, T. & Nyberg, B. 1972 Rate of depurination of native deoxyribonucleic acid. Biochemistry 11, 3610–3618. (doi:10.1021/bi00769a018) 23 Deagle, B. E., Eveson, J. P. & Jarman, S. N. 2006 Quantification of damage in DNA recovered from highly degraded samples—a case study on DNA in faeces. Front. Zool. 3, 11. (doi:10.1186/1742-9994-3-11) 24 Schwarz, C., Debruyne, R., Kuch, M., McNally, E., Schwarcz, H., Aubrey, A. D., Bada, J. & Poinar, H. 2009 New insights from old bones: DNA preservation and degradation in permafrost preserved mammoth remains. Nucleic Acids Res. 37, 3215–3229. (doi:10. 1093/nar/gkp159) 25 Brotherton, P., Endicott, P., Sanchez, J. J., Beaumont, M., Barnett, R., Austin, J. & Cooper, A. 2007 Novel high-resolution characterization of ancient DNA reveals C . U-type base modification events as the sole cause of post mortem miscoding lesions. Nucleic Acids Res. 35, 5717–5728. (doi:10.1093/nar/gkm588) 26 Adler, C. J., Haak, W., Donlon, D., Cooper, A. & The Genographic Consortium. 2011 Survival and recovery of DNA from ancient teeth and bones. J. Arch. Sci. 38, 956 –964. (doi:10.1016/J.Jas.2010.11.010) 27 Ottoni, C., Koon, H., Collins, M. J., Penkman, K. E. H., Rickards, O. & Craig, O. E. 2009 Preservation of ancient DNA in thermally damaged archaeological bone. Naturwissenschaften 96, 267–278. (doi:10.1007/s00114-008-0478-5) 28 Hansen, A. J., Mitchell, D. L., Wiuf, C., Panikert, L., Brand, T. B., Binladen, J., Gilichinsky, D. A., Ronn, R. & Willerslev, E. 2006 Crosslinks rather than strand breaks determine access to ancient DNA sequences from frozen sediments. Genetics 173, 1175–1179. (doi:10.1534/ genetics.106.057349) 29 Sczepanski, J. T., Wong, R. S., McKnight, J. N., Bowman, G. D. & Greenberg, M. M. 2010 Rapid DNA-protein cross-linking and strand scission by an abasic site in a nucleosome core particle. Proc. Natl Acad. Sci. USA 107, 22 475 –22 480. (doi:10.1073/ pnas.1012860108) 30 Briggs, A. W. et al. 2007 Patterns of damage in genomic DNA sequences from a Neandertal. Proc. Natl Acad. Sci. USA 104, 14 616 –14 621. (doi:10.1073/pnas.070 4665104)
31 Sawyer, S., Krause, J., Guschanski, K., Savolainen, V. & P ̈a ̈abo, S. 2012 Temporal patterns of nucleotide misincorporations and DNA fragmentation in ancient DNA. PLoS ONE 7, e34131. (doi:10.1371/journal.pone. 0034131) 32 Orlando, L. et al. 2011 True single-molecule DNA sequencing of a Pleistocene horse bone. Genome Res. 21, 1705 –1719. (doi:10.1101/gr.122747.111) 33 Overballe-Petersen, S., Orlando, L. & Willerslev, E. 2012 Next-generation sequencing offers new insights into DNA degradation. Trends Biotech. 30, 364 –368. (doi:10.1016/j.tibtech.2012.03.007) 34 Dobberstein, R. C., Huppertz, J., von Wurmb-Schwark, N. & Ritz-Timme, S. 2008 Degradation of biomolecules in artificially and naturally aged teeth: implications for age estimation based on aspartic acid racemization and DNA analysis. Forensic. Sci. Int. 179, 181 –191. (doi:10. 1016/j.forsciint.2008.05.017) 35 Gilbert, M. T. P., Hansen, A. J., Willerslev, E., Rudbeck, L., Barnes, I., Lynnerup, N. & Cooper, A. 2003 Characterization of genetic miscoding lesions caused by postmortem damage. Am. J. Hum. Genet. 72, 48–61. (doi:10.1086/345379) 36 Hagelberg, E., Bell, L. S., Allen, T., Boyde, A., Jones, S. J. & Clegg, J. B. 1991 Analysis of ancient bone DNA: techniques and applications. Phil. Trans. R. Soc. Lond. B 333, 399 –407. (doi:10.1098/rstb.1991.0090) 37 Hoss, M., Jaruga, P., Zastawny, T. H., Dizdaroglu, M. & P ̈a ̈abo, S. 1996 DNA damage and DNA sequence retrieval from ancient tissues. Nucleic Acids Res. 24, 1304–1307. (doi:10.1093/nar/24.7.1304) 38 P ̈a ̈abo, S., Higuchi, R. G. & Wilson, A. C. 1989 Ancient DNA and the polymerase chain-reaction: the emerging field of molecular archaeology. J. Biol. Chem. 264, 9709 – 9712. 39 Poinar, H. N., Hoss, M., Bada, J. L. & Pa ̈a ̈bo, S. 1996 Amino acid racemization and the preservation of ancient DNA. Science 272, 864 –866. (doi:10.1126/science.272. 5263.864) 40 Haynes, S., Searle, J. B., Bertman, A. & Dobney, K. M. 2002 Bone preservation and ancient DNA: the application of screening methods for predicting DNA survival. J. Arch. Sci. 29, 585 –592. (doi:10.1006/jasc. 2001.0731) 41 Colson, I. B., Bailey, J. F., Vercauteren, M. & Sykes, B. C. 1997 The preservation of ancient DNA and bone diagenesis. Ancient Biomol. 1, 109–117. 42 Poulakakis, N., Parmakelis, A., Lymberakis, P., Mylonas, M., Zouros, E., Reese, D. S., Glaberman, S. & Caccone, A. 2007 It remains a mammoth DNA fragment: a reply to Binladen et al. (2006) and Orlando et al. (2006). Biol. Lett. 3, 60–63. (doi:10.1098/rsbl.2006.0555) 43 Zink, A. R. & Nerlich, A. G. 2005 Long-term survival of ancient DNA in Egypt: reply to Gilbert et al. Am. J. Phys. Anthropol. 128, 115– 118. (doi:10.1002/ajpa.20046) 44 Hummel, S. 2001 Encyclopedia of life sciences. Chichester, UK: John Wiley & Sons, Ltd. 45 Bar, W., Kratzer, A., Machler, M. & Schmid, W. 1988 Postmortem stability of DNA. Forensic Sci. Int. 39, 59–70. (doi:10.1016/0379-0738(88)90118-1) 46 Campos, P. F., Craig, O. E., Turner-Walker, G., Peacock, E., Willerslev, E. & Gilbert, M. T. P. 2012 DNA in ancient bone: where is it located and how should we extract it? Ann. Anat. Anatomischer Anzeiger 194, 7– 16. (doi:10.1016/j.aanat.2011.07.003) 47 Willerslev, E. et al. 2007 Ancient biomolecules from deep ice cores reveal a forested Southern Greenland. Science 317, 111 –114. (doi:10.1126/Science.1141758) 48 Buckley, M et al. 2008 Comment on ‘Protein sequences from mastodon and Tyrannosaurus rex revealed by mass
4732 M. E. Allentoft et al. The half-life of DNA
Proc. R. Soc. B (2012)


spectrometry’. Science 319, 33; author reply. (doi:10. 1126/science.1147046) 49 Jans, M., Nielsen-Marsh, C., Smith, C., Collins, M. J. & Kars, H. 2004 Characterisation of microbial attack on archaeological bone. J. Arch. Sci. 31, 87– 95. (doi:10. 1016/j.jas.2003.07.007) 50 Pruvost, M., Schwarz, R., Correia, V. B., Champlot, S., Braguier, S., Morel, N., Fernandez-Jalvo, Y., Grange, T. & Geigl, E.-M. 2007 Freshly excavated fossil bones are best for amplification of ancient DNA. Proc. Natl Acad. Sci. USA 104, 739–744. (doi:10.1073/pnas. 0610257104) 51 Allentoft, M. E. et al. 2009 Identification of microsatellites from an extinct moa species using high-throughput (454) sequence data. Biotechniques 46, 195–200. (doi:10.2144/000113086) 52 Allentoft, M. E., Bunce, M., Scofield, R. P., Hale, M. L. & Holdaway, R. N. 2010 Highly skewed sex ratios and biased fossil deposition of moa: ancient DNA provides new insight on New Zealand’s extinct megafauna. Q. Sci. Rev. 29, 753–762. (doi:10.1016/j.quascirev.2009.11.022) 53 Allentoft, M. E., Scofield, R. P., Oskam, C. L., Hale, M. L., Holdaway, R. N. & Bunce, M. 2012 A molecular characterisation of a newly discovered megafaunal fossil site in North Canterbury, South Island, New Zealand. J. R. Soc. N Z. (doi:10.1080/03036758.2011.574821) 54 Holdaway, R. N. & Worthy, T. H. 1997 A reappraisal of the late quaternary fossil vertebrates of Pyramid Valley Swamp, North Canterbury, New Zealand. N Z J. Zool. 24, 69–121. (doi:10.1080/03014223.1997.9518107) 55 Bunce, M., Worthy, T. H., Ford, T., Hoppitt, W., Willerslev, E., Drummond, A. & Cooper, A. 2003 Extreme reversed sexual size dimorphism in the extinct New Zealand moa Dinornis. Nature 425, 172–175. (doi:10.1038/nature01871) 56 Allentoft, M. E. et al. 2011 Profiling the dead: generating microsatellite data from fossil bones of extinct megafauna—protocols, problems, and prospects. PLoS ONE 6, e16670. (doi:10.1371/journal.pone.0016670)
57 Oskam, C. L. et al. 2010 Fossil avian eggshell preserves ancient DNA. Proc. R. Soc. B 277, 1991 –2000. (doi:10.1098/rspb.2009.2019) 58 Bunce, M., Oskam, C. L. & Allentoft, M. E. 2012 The use of quantitative real-time PCR in ancient DNA research. In Methods in molecular biology—ancient DNA (eds B. Shapiro & M. Hofreiter), pp. 121– 132. New York, NY: Humana Press. (doi:10.1007/978-1-61779516-9_16) 59 Rohland, N. & Hofreiter, M. 2007 Comparison and optimization of ancient DNA extraction. Biotechniques 42, 343–352. (doi:10.2144/000112383) 60 Noonan, J. P. et al. 2005 Genomic sequencing of Pleistocene cave bears. Science 309, 597 –600. (doi:10.1126/ science.1113485) 61 Poinar, H. N. et al. 2006 Metagenomics to paleogenomics: large-scale sequencing of mammoth DNA. Science 311, 392– 394. (doi:10.1126/science.1123360) 62 Cooper, A., Lalueza-Fox, C., Anderson, S., Rambaut, A., Austin, J. & Ward, R. 2001 Complete mitochondrial genome sequences of two extinct moas clarify ratite evolution. Nature 409, 704 –707. (doi:10.1038/35055536) 63 Li, H. & Durbin, R. 2009 Fast and accurate short read alignment with Burrows– Wheeler transform. Bioinformatics 25, 1754–1760. (doi:10.1093/bioinformatics/ btp324) 64 Berna, F., Matthews, A. & Weiner, S. 2004 Solubilities of bone mineral from archaeological sites: the recrystallization window. J. Archaeol. Sci. 31, 867 –882. (doi:10. 1016/j.jas.2003.12.003) 65 Molak, M. & Ho, S. 2011 Evaluating the impact of postmortem damage in ancient DNA: a theoretical approach. J. Mol. Evol. 73, 244–255. (doi:10.1007/s00239-0119474-z) 66 Dabney, J. & Meyer, M. 2012 Length and GC-biases during sequencing library amplification: a comparison of various polymerase-buffer systems with ancient and modern DNA sequencing libraries. Biotechniques 52, 87– 94. (doi:10.2144/000113809)
The half-life of DNA M. E. Allentoft et al. 4733
Proc. R. Soc. B (2012)