Abstract:1. Motivation Under favourable conditions DNA molecules can persist for hundreds of thousands of years. Such genetic remains make up invaluable resources to study past assemblages, populations, and even the evolution of species. However, DNA is subject to degradation, and hence over time decrease to ultra low concentrations which makes it highly prone to contamination by modern sources. Strict precautions are therefore necessary to ensure that DNA from modern sources does not appear in the final data is authenticated as ancient. The most generally accepted and widely applied authenticity for ancient DNA studies is to test for elevated deaminated cytosine residues towards the termini of the molecules: DNA damage. To date, this has primarily been used for single organisms and recently for read assemblies, however, these methods are not applicable for estimating DNA damage for ancient metagenomes with tens and even hundreds of thousands of species. 2. Methods We present metaDMG, a novel framework and toolkit that allows for the estimation, quantification and visualization of postmortem damage for single reads, single genomes and even metagenomic environmental DNA by utilizing the taxonomic branching structure. It bypasses any need for initial classification, splitting reads by individual organisms, and realignment. We have implemented a Bayesian approach that combines a modified geometric damage profile with a beta-binomial model to fit the entire model to the individual misincorporations at all taxonomic levels. 3. Results We evaluated the performance using both simulated and published environmental DNA datasets and compared to existing methods when relevant. We find metaDMG to be an order of magnitude faster than previous methods and more accurate – even for complex metagenomes. Our simulations show that metaDMG can estimate DNA damage at taxonomic levels down to 100 reads, that the estimated uncertainties decrease with increased number of reads and that the estimates are more significant with increased number of C to T misincorporations. 4. Conclusion metaDMG is a state-of-the-art program for aDNA damage estimation and allows for the computation of nucleotide misincorporation, GC-content, and DNA fragmentation for both simple and complex ancient genomic datasets, making it a complete package for ancient DNA damage authentication. ### Competing Interest Statement The authors have declared no competing interest.

AuthentiCT: a model of ancient DNA damage to estimate the proportion of present-day DNA contamination

PyDamage: automated ancient damage identification and estimation for contigs in ancient DNA de novo assembly

Assessing the impact of post-mortem damage and contamination on imputation performance in ancient DNA

DamageProfiler: Fast damage pattern calculation for ancient DNA

Estimation of DNA Degradation in Archaeological Human Remains

mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters

A likelihood method for estimating present-day human contamination in ancient DNA samples using low-depth haploid chromosome data

A likelihood method for estimating present-day human contamination in ancient male samples using low-depth X-chromosome data

Schmutzi: estimation of contamination and endogenous mitochondrial consensus calling for ancient DNA

metaDMG – A Fast and Accurate Ancient DNA Damage Toolkit for Metagenomic Data

Revisiting the Briggs Ancient DNA Damage Model: A Fast Maximum Likelihood Method to Estimate Post-Mortem Damage.

New Insights Into Mitochondrial DNA Reconstruction and Variant Detection in Ancient Samples

Revisiting the Briggs ancient DNA damage model: a fast regression method to estimate postmortem damage

Reducing microbial and human contamination in DNA extractions from ancient bones and teeth.

Separating endogenous ancient DNA from modern day contamination in a Siberian Neandertal

Competitive mapping allows for the identification and exclusion of human DNA contamination in ancient faunal genomic datasets

Read Length Dominates Phylogenetic Placement Accuracy of Ancient DNA Reads

An open-sourced bioinformatic pipeline for the processing of Next-Generation Sequencing derived nucleotide reads: Identification and authentication of ancient metagenomic DNA

Mitochondrial DNA Consensus Calling and Quality Filtering for Constructing Ancient Human Mitogenomes: Comparison of Two Widely Applied Methods

AntCaller: an Accurate Variant Caller Incorporating Ancient DNA Damage

Improved detection of methylation in ancient DNA