Sep 07, 2023

Public workspaceA method for the temperature-controlled extraction of DNA from ancient bones V.2

  • 1Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Deutscher Platz 6, D-04103 Leipzig, Germany;
  • 2EMBL-EBI, Wellcome Genome Campus, Hinxton, Cambridgeshire, CB10 1SD, UK
Icon indicating open access to content
QR code linking to this content
Protocol CitationElena Essel, Petra Korlevic, Matthias Meyer 2023. A method for the temperature-controlled extraction of DNA from ancient bones . protocols.io https://dx.doi.org/10.17504/protocols.io.rm7vz3pb2gx1/v2Version created by Elena Essel
Manuscript citation:
A method for the temperature-controlled extraction of DNA from ancient bones Elena Essel, Petra Korlević, and Matthias Meyer BioTechniques 2021 71:1, 382-386 https://doi.org/10.2144/btn-2021-0025
License: This is an open access protocol distributed under the terms of the Creative Commons Attribution License,  which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Protocol status: Working
We use this protocol and it's working
Created: September 07, 2023
Last Modified: April 10, 2024
Protocol Integer ID: 87493
Keywords: Ancient DNA, sequential DNA extraction, contamination removal, endogenous DNA, archaeological material
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.
Materials
Reagents
ReagentSodium phosphate, 0.5M buffer soln., pH 7.0 Thermo ScientificCatalog #AAJ63791AP
ReagentWater for HPLCSigma AldrichCatalog #270733
ReagentEDTA solution pH 8.0 (0.5 M) for molecular biologyAppliChemCatalog #A4892,1000
ReagentTris buffer pH 8.0 (1 M) for molecular biologyAppliChemCatalog #A4577,0500
ReagentProteinase K 100 mgSigma-aldrichCatalog #3115879001
ReagentTWEEN® 20Sigma AldrichCatalog #T2700-100ML

Consumables and equipment
ReagentDNA LoBind Tubes 2.0 mLEppendorfCatalog #0030108078
ReagentDNA LoBind Tubes 2.0 mLEppendorfCatalog #0030108078
ReagentCeramic beads 2.8 mm VWR InternationalCatalog #432-0292
Reagent50 ml CELLSTAR® Polypropylene Tube 30/115 MM Conical Bottom Blue screw cap sterile skirtgreiner bio-oneCatalog #210261
ReagentParafilm M 10 cm widneoLabCatalog #3-1012


Equipment
Thermomixer
NAME
HLC
BRAND
52 82 00133
SKU

Equipment
Incubator
NAME
Memmert
BRAND
Incubator IN55
SKU

Equipment
Tube rotator
NAME
VWR
BRAND
444-0500
SKU

Equipment
UV cross-linker
NAME
Vilber
BRAND
Bio-Link BLX 254
SKU

Equipment
Vortex mixer
NAME
Scientific Industries
BRAND
SI-0236
SKU

Equipment
Centrifuge
NAME
Bench centrifuge
TYPE
Eppendorf
BRAND
5424
SKU


Protocol materials
ReagentTWEEN® 20Merck MilliporeSigma (Sigma-Aldrich)Catalog #T2700-100ML
In Materials and 3 steps
ReagentDNA LoBind Tubes 2.0 mLEppendorfCatalog #0030108078
In Materials, Materials
ReagentProteinase K 100 mgMerck MilliporeSigma (Sigma-Aldrich)Catalog #3115879001
Materials, Step 4
Reagent50 ml CELLSTAR® Polypropylene Tube 30/115 MM Conical Bottom Blue screw cap sterile skirtgreiner bio-oneCatalog #210261
Materials
ReagentTris buffer pH 8.0 (1 M) for molecular biologyAppliChemCatalog #A4577,0500
Materials, Step 3
ReagentCeramic beads 2.8 mm VWR InternationalCatalog #432-0292
Materials, Step 6
ReagentSodium phosphate, 0.5M buffer soln., pH 7.0 Thermo ScientificCatalog #AAJ63791AP
Materials, Step 2
ReagentWater for HPLCMerck MilliporeSigma (Sigma-Aldrich)Catalog #270733
In Materials and 2 steps
ReagentParafilm M 10 cm widneoLabCatalog #3-1012
Materials
ReagentEDTA solution pH 8.0 (0.5 M) for molecular biologyAppliChemCatalog #A4892,1000
Materials, Step 4
Buffer preparation
Buffer preparation

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


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

Amount49.5 mL ReagentSodium phosphate, 0.5M buffer soln., pH 7.0 Thermo ScientificCatalog #AAJ63791AP
Amount50 µL ReagentTWEEN® 20Sigma AldrichCatalog #T2700-100ML
Tris-Tween wash buffer (10 mM Tris-HCl, pH 8.0, 0.1% Tween-20) is prepared by combining the following reagents:
Amount49.5 mL ReagentWater for HPLCSigma AldrichCatalog #270733
Amount0.5 mL ReagentTris buffer pH 8.0 (1 M) for molecular biologySigma AldrichCatalog #A4577,0500
Amount50 µL ReagentTWEEN® 20Sigma AldrichCatalog #T2700-100ML
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:
Amount3.725 mL ReagentWater for HPLCSigma AldrichCatalog #270733
Amount45 mL ReagentEDTA solution pH 8.0 (0.5 M) for molecular biologySigma AldrichCatalog #A4892,1000
Amount25 µL ReagentTWEEN® 20Sigma AldrichCatalog #T2700-100ML
Amount1.25 mL 10 mg/ml proteinase K solution in water (prepared from ReagentProteinase K 100 mgSigma AldrichCatalog #3115879001 )

Note
Proteinase K is added after UV irradiation



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

To facilitate resuspension of the bone powder during the subsequent incubation and wash steps, add 3-4ReagentCeramic beads 2.8 mm Sigma AldrichCatalog #432-0292 to the sample material.

Temperature-controlled phosphate treatment
Temperature-controlled phosphate treatment
Add Amount0.5 mL 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 Shaker900 rpm, 00:15:00.

Note
Temperature-controlled phosphate treatment steps
Temperature37 °C 2 times
Temperature60 °C 2 times
Temperature90 °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).

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

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.


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.

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

Transfer tubes to a tabletop centrifuge and spin for 2 min at maximum speed (e.g., 16,400g/13,200 rpm)
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
Final digestion of sample material

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


Note
Wrap the tube with parafilm to prevent leaking.



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

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
DNA purification of phosphate fractions and final lysate

Thaw the sodium phosphate fractions (and lysates if necessary) at Temperature37 °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.

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.

Library preparation, sequencing, and data processing
Library preparation, sequencing, and data processing
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.


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.


Sequence analysis
Sequence analysis

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.

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).

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.).

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

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).

Citations
Step 18
Glocke I, Meyer M. Extending the spectrum of DNA sequences retrieved from ancient bones and teeth.
https://doi.org/10.1101/gr.219675.116
Step 19
Gansauge MT, Aximu-Petri A, Nagel S, Meyer M. 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.
https://doi.org/10.1038/s41596-020-0338-0
Step 20
Kircher M, Sawyer S, Meyer M. Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform.
https://doi.org/10.1093/nar/gkr771
Step 21
Renaud G, Stenzel U, Kelso J. leeHom: adaptor trimming and merging for Illumina sequencing reads.
https://doi.org/10.1093/nar/gku699
Step 22
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. A high-coverage genome sequence from an archaic Denisovan individual.
https://doi.org/10.1126/science.1224344
Step 22
Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform.
https://doi.org/10.1093/bioinformatics/btp698
Step 25
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genome Project Data Processing Subgroup.. The Sequence Alignment/Map format and SAMtools.
https://doi.org/10.1093/bioinformatics/btp352