A method for the temperature-controlled extraction of DNA from ancient bones

Elena Essel, Matthias Meyer, Petra Korlevic

Published: 2021-09-10 DOI: 10.17504/protocols.io.bx5epq3e

Disclaimer

DISCLAIMER – FOR INFORMATIONAL PURPOSES ONLY; USE AT YOUR OWN RISK

The protocol content here is for informational purposes only and does not constitute legal, medical, clinical, or safety advice, or otherwise; content added to protocols.io is not peer reviewed and may not have undergone a formal approval of any kind. Information presented in this protocol should not substitute for independent professional judgment, advice, diagnosis, or treatment. Any action you take or refrain from taking using or relying upon the information presented here is strictly at your own risk. You agree that neither the Company nor any of the authors, contributors, administrators, or anyone else associated with protocols.io, can be held responsible for your use of the information contained in or linked to this protocol or any of our Sites/Apps and Services.

Abstract

We here provide a protocol for the decontamination of ancient bones and teeth that is based on a temperature-controlled, sequential release of DNA. DNA can be extracted from all fractions generated with this method and the fraction with the highest proportion of endogenous DNA identified for further analysis. The protocol proceeds through repeated incubation of the sample powder in phosphate buffer at 37, 60 and 90 °C, followed by the complete lysis of the residual sample powder. As DNA is denatured at high temperature, subsequent DNA extraction and library preparation has to be performed using methods optimized for single-stranded DNA.

Steps

Buffer preparation

1.

Note
All buffers are irradiated with UV-C light at a dose of 7 kJ/cm2 using a cross-linker.

2.

Sodium-phosphate buffer (0.5 M sodium phosphate, pH 7.0, 0.1 % Tween 20) is prepared by combining the following reagents:

49.5mL

50µL

3.

Tris-Tween wash buffer (10 mM Tris-HCl, pH 8.0, 0.1% Tween-20) is prepared by combining the following reagents:

49.5mL

0.5mL

50µL

4.

Lysis buffer (0.45 M EDTA, pH 8.0, 0.05% Tween-20 and 0.25 mg/ml proteinase K) is prepared by combining the following reagents:

3.725mL

45mL

25µL

1.25mL 10 mg/ml proteinase K solution in water (prepared from )

Note
Proteinase K is added after UV irradiation

Sample preparation

5.

In an ancient DNA cleanroom, remove approximately 50mg of sample powder from each specimen using a sterile dentist drill and transfer the powder to a 2.0 ml DNA LoBind tube.

6.

To facilitate resuspension of the bone powder during the subsequent incubation and wash steps, add 3-4 to the sample material.

Temperature-controlled phosphate treatment

7.

Add 0.5mL sodium phosphate buffer to the sample powder, completely resuspend the powder by thorough vortexing, and incubate the tube in a thermo block adjusted to the desired temperature 900rpm

Note
Temperature-controlled phosphate treatment steps37°C 2 times60°C 2 times90°C 2 times

Note
At least one negative control (tube without sample material) should be included in each experiment and carried through all subsequent steps).

8.

Transfer tubes to a tabletop centrifuge and spin for 2 min at maximum speed (e.g., 16,400g/13,200 rpm).

9.

Transfer supernatant to a 1.5 mL LoBind tube and store at -20 °C until the day of DNA extraction.

Note
Beads facilitate the resuspension of the sample powder after centrifugation steps, but make it harder to remove supernatant.Pipette slowly and carefully.

10.

Repeat steps 7-9 once at each temperature (for a total of 2 wash steps).

Note
For the 90 °C incubation, make sure the liquid in the tube reaches 90 °C by the end of the 15 min incubation time. If necessary, set the thermo block to a higher temperature.

11.

The temperature-controlled phosphate treatment is followed by a room-temperature wash step with 1mL Tris-Tween buffer at the end of the last temperature cycle. Completely resuspend the powder by thorough vortexing.

12.

Transfer tubes to a tabletop centrifuge and spin for 2 min at maximum speed (e.g., 16,400g/13,200 rpm)

13.

Transfer supernatant to a 1.5 mL LoBind tube and store at -20 °C until the day of DNA extraction.

Final digestion of sample material

14.

Add 1mL of lysis buffer to the sample powder, completely resuspended the powder by vortexing, and incubate overnight (8 – 16 h) with rotation at 37°C

Note
Wrap the tube with parafilm to prevent leaking.

15.

Transfer tubes to a tabletop centrifuge and spin for 2 min at maximum speed (commonly at 16,400 g/13,200 rpm).

16.

Transfer supernatant to a 1.5 mL LoBind tube and proceed to DNA extraction or store the tube at -20 °C until the day of DNA extraction.

DNA purification of phosphate fractions and final lysate

17.

Thaw the sodium phosphate fractions (and lysates if necessary) at 37°C in a thermo block with gentle shaking.

Note
Make sure the liquid is fully thawed and any crystals have completely dissolved.

Note
If desired, DNA extraction can also be performed on the Tris-Tween buffer, but DNA yields are expected to be extremely low.

18.

For the sodium phosphate fractions, purify 100 µl of the supernatant, and for the final lysate, purify 500 µl using binding buffer ‘G’ of the DNA extraction method described in Glocke and Meyer (2017). Final volume of all DNA extracts is 50 µl.

Citation
Glocke I, Meyer M 2017 Extending the spectrum of DNA sequences retrieved from ancient bones and teeth. Genome research https://doi.org/10.1101/gr.219675.116

Library preparation, sequencing, and data processing

19.

Prepare DNA libraries using 20% of the DNA extract as input, following the protocol for library preparation, quantification and indexing by Gansauge et al. (2020).

Citation
Gansauge MT, Aximu-Petri A, Nagel S, Meyer M 2020 Manual and automated preparation of single-stranded DNA libraries for the sequencing of DNA from ancient biological remains and other sources of highly degraded DNA. Nature protocols https://doi.org/10.1038/s41596-020-0338-0

20.

Perform shallow shotgun sequencing on Illumina’s MiSeq or HiSeq2500 platforms (or other Illumina platforms) using a paired-end double-index configuration (2x 76 + 2x 7 cycles).

Citation
Kircher M, Sawyer S, Meyer M 2012 Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform. Nucleic acids research https://doi.org/10.1093/nar/gkr771

Sequence analysis

21.

Trim adapters and merge overlapping paired-end reads into single-molecule sequences using leeHom.

Citation
Renaud G, Stenzel U, Kelso J 2014 leeHom: adaptor trimming and merging for Illumina sequencing reads. Nucleic acids research https://doi.org/10.1093/nar/gku699

22.

Use the Burrows-Wheeler Aligner (BWA, https://github.com/mpieva/network-aware-bwa) to align merged sequences to a suitable reference genome (e.g. turTru1.75, bosTauUMD3.1, loxAfr4) using ancient parameters (“-n 0.01 –o 2 –l 16500”) allowing more mismatches and indels.

Citation
Li H, Durbin R 2010 Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics (Oxford, England) https://doi.org/10.1093/bioinformatics/btp698

Citation
Meyer M, Kircher M, Gansauge MT, Li H, Racimo F, Mallick S, Schraiber JG, Jay F, Prüfer K, de Filippo C, Sudmant PH, Alkan C, Fu Q, Do R, Rohland N, Tandon A, Siebauer M, Green RE, Bryc K, Briggs AW, Stenzel U, Dabney J, Shendure J, Kitzman J, Hammer MF, Shunkov MV, Derevianko AP, Patterson N, Andrés AM, Eichler EE, Slatkin M, Reich D, Kelso J, Pääbo S 2012 A high-coverage genome sequence from an archaic Denisovan individual. Science (New York, N.Y.) https://doi.org/10.1126/science.1224344

23.

Restrict further analyses to sequences of length 35 bp and above to avoid spurious alignments of short sequences with random similarity to the reference genome.

24.

Merge sequences with the same start- and end-coordinate into one consensus sequence using bam-rmdup (https://github.com/mpieva/biohazard-tools).

25.

Generate summary statistics using samtools and choose the library with the highest proportion of endogenous DNA for further sequencing. Prepare additional libraries from remaining DNA extract if necessary.

Citation
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genome Project Data Processing Subgroup. 2009 The Sequence Alignment/Map format and SAMtools. Bioinformatics (Oxford, England) https://doi.org/10.1093/bioinformatics/btp352

推荐阅读

Nature Protocols
Protocols IO
Current Protocols
扫码咨询