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.

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

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

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

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

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

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

DamageProfiler: Fast damage pattern calculation for ancient DNA

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

Modeling Kinetic Rate Variation in Third Generation DNA Sequencing Data to Detect Putative Modifications to DNA Bases

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

gargammel: a sequence simulator for ancient DNA

Read Length Dominates Phylogenetic Placement Accuracy of Ancient DNA Reads

Unravelling reference bias in ancient DNA datasets

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

AntCaller: an Accurate Variant Caller Incorporating Ancient DNA Damage

PCR‐Based Analysis of Mitochondrial DNA Copy Number, Mitochondrial DNA Damage, and Nuclear DNA Damage

Benchmarking software tools for trimming adapters and merging next-generation sequencing data for ancient DNA

Detecting DNA Modifications from SMRT Sequencing Data by Modeling Sequence Context Dependence of Polymerase Kinetic.

Improved detection of methylation in ancient DNA

Estimation of DNA Degradation in Archaeological Human Remains

NGSNGS: next-generation simulator for next-generation sequencing data