Abstract
Mitochondrial DNA (mtDNA) analysis is crucial for understanding human population structure and genetic diversity. However, post-mortem DNA damage poses challenges, that make analysis difficult. DNA preservation is affected by environmental conditions which, among other factors, complicates the differentiation of endogenous variants from artefacts in ancient mtDNA mix profiles. This study aims to develop a molecular damage profile for ancient mtDNA that can become a useful tool in analysing mtDNA from ancient remains. A dataset of 427 whole genomes or capture of mtDNA sequences from individuals representing different historical periods and climatic regions was compiled from the ENA database. Present-day and UDG-treated ancient samples were also included and used to establish levels of damaged reads. Results indicated that samples from cold regions exhibited the lowest percentage of damaged reads, followed by arid, cold, tropical and temperate regions, with significant differences observed between cold and temperate regions. A global damage profile was generated, identifying 2933 positions (25% of the positions considered) with damage in more than 23.8% of the samples analysed, deemed as damage hotspots. Notably, 2856 of these hotspots had never been reported as damage or mutational hotspots, or heteroplasmic positions. Damage hotspot frequency by position was slightly higher in the non-coding region compared with the coding region. In conclusion, this study provides a molecular damage profile for ancient mtDNA analysis that is expected to be a valuable tool in the interpretation of mtDNA variation in ancient samples.