MTaxi: A comparative tool for taxon identification of ultra low coverage ancient genomes

2023-01-01
Atağ, Gözde
VURAL, KIVILCIM BAŞAK
KAPTAN, DAMLA
Özkan, Mustafa
Koptekin, Dilek
Sağlıcan, Ekin
Doğramacı, Sevcan
Köz, Mevlüt
Yılmaz, Ardan
Söylev, Arda
Togan, İnci Zehra
Somel, Mehmet
Özer, Füsun
A major challenge in zooarchaeology is to morphologically distinguish closely related species’ remains, especially using small bone fragments. Shotgun sequencing aDNA from archeological remains and comparative alignment to the candidate species’ reference genomes will only apply when reference nuclear genomes of comparable quality are available, and may still fail when coverages are low. Here, we propose an alternative method, MTaxi, that uses highly accessible mitochondrial DNA (mtDNA) to distinguish between pairs of closely related species from ancient DNA sequences. MTaxi utilises mtDNA transversion-type substitutions between pairs of candidate species, assigns reads to either species, and performs a binomial test to determine the sample taxon. We tested MTaxi on sheep/goat and horse/donkey data, between which zooarchaeological classification can be challenging in ways that epitomise our case. The method performed efficiently on simulated ancient genomes down to 0.3x mitochondrial coverage for both sheep/goat and horse/donkey, with no false positives. Trials on n=18 ancient sheep/goat samples and n=10 horse/donkey samples of known species identity also yielded 100% accuracy. Overall, MTaxi provides a straightforward approach to classify closely related species that are difficult to distinguish through zooarchaeological methods using low coverage aDNA data, especially when similar quality reference genomes are unavailable. MTaxi is freely available at https://github.com/goztag/MTaxi.
Open Research Europe
Citation Formats
G. Atağ et al., “MTaxi: A comparative tool for taxon identification of ultra low coverage ancient genomes,” Open Research Europe, vol. 2, pp. 0–0, 2023, Accessed: 00, 2023. [Online]. Available: https://www.scopus.com/inward/record.uri?partnerID=HzOxMe3b&scp=85164722199&origin=inward.