Hit or miss - A metagenomic evaluation of intra-bone variability of host pathogen load in tuberculosis-infected human remains

Many sampling protocols have been established to successfully retrieve human DNA from archaeological remains, however the systematic detection of ancient pathogens remains challenging. Here, we present a first assessment of the intra-bone variability of metagenomic composition in human skeletal remains and its effect on the sampling success for Mycobacterium tuberculosis (MTB) and human endogenous DNA. For this purpose, four bone samples from published peer-reviewed studies with PCR-based evidence for ancient MTB DNA were selected. Two bone samples of a Neolithic individual from Halberstadt, Germany and two ribs of two 18th-century Hungarian church mummies were sampled at multiple locations for equal amounts, followed by DNA extraction and library construction. Shotgun sequencing data was generated for taxonomic profiling as well as quantitative and qualitative evaluation of MTB and human endogenous DNA. Despite low variance in microbial diversity within and across samples, intra-bone variability of up to 36.45- and 62.88-fold for authentic ancient MTB and human reads, respectively, was detected. This study demonstrates the variable sampling success for MTB and human endogenous DNA within single skeletal samples despite relatively consistent microbial composition and highlights how a multisampling approach can facilitate the detection of hotspots with highly concentrated pathogen and human endogenous DNA.