On 2018-03-21 19:13:17, user James Fellows Yates wrote:
We find this novel application of ancient dental calculus metagenomes an appealing example of how historical samples can be used to demonstrate host-associated microbial evolution. However we would like to make a few suggestions regarding the use of ancient microbial data in this study.
The following are three major recommendations that we feel will strengthen the results of this aspect of the study:
1) One major challenge in analysing ancient DNA is false positive identification of taxa due to the presence of modern contaminating environmental organisms (Warinner et al. 2017, Ann. Rev. Hu. Gen). We suggest to run your TM7-identified reads through mapDamage2 (https://ginolhac.github.io/..., Jónsson et al. 2013, Bioinformatics) to help authenticate that the TM7 DNA that you have detected in the ancient samples is indeed ‘ancient’. The program will generate plots that should show elevated cytosine to thymine deamination patterns at the termini of the fragmented DNA if the reads are truly ancient (Sawyer et al. 2012, PLoS One). This is particularly important for the data from Weyrich et al., who reported substantial soil contamination in their Neanderthal calculus samples (2017, Nature, Supplemental tables S5, S7). This would then be clearer evidence of age-related DNA damage than Supplementary Table 4, as referred to in Figure 2.
2) Related to this, it would be useful to the reader to justify why a 1% or 15% genome coverage is enough for proof of the presence of this species (assuming that ‘1% mapped’ in the section ‘Reduced genomes from Environment…’ refers to 1% genome coverage, rather than 1% of reads in the library). The high risk of false positives resulting from mismapped and/or modern contamination may play a role here, and TM7-identified reads may come from environmental DNA not originating from the individual’s oral cavity. Additionally, while there is a paucity of data for the TM7 strain from the Neanderthal sample (‘1%’), further evidence to show this strain is indeed from the oral cavity (such the presence of a marker gene) would be useful confirmation here that this strain is authentic and not derived from a soil TM7 ; such as with the tree for the medieval sample. Furthermore, is the low percentage of TM7-identified reads also the reason why genome assembly was performed only on the Warinner et al. 2014 (Nat. Genet.) data, despite the El Sidron1 sample being described as ‘well preserved’? Clarification on how the preservation of B61 and El Sidron1 was assessed, as well as justification for using one or the other, rather than both, throughout the manuscript would be welcome.
3) Finally, the suggestion that human acquisition of TM7 during animal domestication is highly speculative. This is demonstrated by an ‘increase’ of mapped reads between a single Neanderthal and a single Medieval sample; however, factors such as sequencing effort, relative abundance, and individual microbiome variation are not taken into account. Far more samples from relevant time periods would be necessary to substantiate such a claim. Thus, we feel the paper would benefit by having this statement removed entirely.
We would finally like to make a few minor comments regarding the structure of the paper:
* Make an individual section for the renaming of the TM7 groups, as this does not stand out despite being a major aim of the paper, according to the last paragraph of the introduction. Make Supplemental Table S2, which lays out the new naming convention, a main text table for this new section.<br />
* Clarify in both the results and the methods that the Warinner et al. 2014 data was used in both assembly AND mapping. In the methods section this data is only described in the mapping section, suggesting a genome bin was generated from mapping, and this may confuse the reader.<br />
* Provide more detailed information in the methods section regarding the parameters used for all software. Currently it is not possible for the reader to reproduce or assess the reliability of the analysis performed.<br />
* Roary cannot be reliably used to cluster genomes across species (Page, et al. 2015, Bioinf.), only within a single species, yet the authors have used it to cluster genomes across a phylum. The title states that the members of the Saccharibacteria phylum are highly diverse, which is countered by the fact that they clustered well by Roary. Can the authors comment on this discrepancy? These results imply that the individual organisms are actually all the same species, which means the naming convention suggested in Supplemental Table S2 will need to be adjusted.<br />
* Several of the Supplemental tables have very small text and are difficult to read and parse. Consider putting them in a standalone spreadsheet file and making the spreadsheet file available as a supplemental file. Additionally, the tree currently displayed to the side of Supplementary Table 1 is offset from the labels in the table itself, and is difficult to visualise.<br />
* In the text, it is stated that B61 has mild to severe periodontal disease; however, the original publication (Warinner et al. 2014) provides detailed oral pathology records and characterizes the individual as having moderate to severe periodontal disease with specific clinical features in the Supplementary Information.<br />
* Please correct the idiosyncratic use of capitalization and italicization when referring to species throughout the text. For example, “Nanosynbacter Lyticus”, “ecoli-like”, “streptococcus thermophilus”, “in vibrio genomes,” etc.<br />
* Citation 12 has been badly formatted
James A. Fellows Yates (fellows[at]shh.mpg.de)<br />
Christina Warinner (warinner[at]shh.mpg.de)<br />
Department of Archaeogenetics, Max Planck Institute for the Science of Human History, Germany
Irina Velsko (ivelsko[at]clemson.edu)<br />
Department of Biological Sciences, Clemson University, USA
References:<br />
Jónsson, H., Ginolhac, A., Schubert, M., Johnson, P. L. F., & Orlando, L. (2013). mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters. Bioinformatics , 29(13), 1682–1684. https://doi.org/10.1093/bio...
Page A.J., Cummins C.A., Hunt M., Wong V.K., Reuter S., Holden M.T.G., Fookes M., Falush D., Keane J.A., Parkhill, J. (2015) Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics,31(22):3691-3693. https://doi.org/10.1093/bio...
Sawyer, S., Krause, J., Guschanski, K., Savolainen, V., & Pääbo, S. (2012). Temporal patterns of nucleotide misincorporations and DNA fragmentation in ancient DNA. PloS One, 7(3), e34131. https://doi.org/10.1371/jou...
Warinner, C., Rodrigues, J. F. M., Vyas, R., Trachsel, C., Shved, N., Grossmann, J., … Cappellini, E. (2014). Pathogens and host immunity in the ancient human oral cavity. Nature Genetics, 46(4), 336–344. https://doi.org/10.1038/ng....
Warinner, C., Herbig, A., Mann, A., Fellows Yates, J. A., Weiß, C. L., Burbano, H. A., … Krause, J. (2017). A Robust Framework for Microbial Archaeology. Annual Review of Genomics and Human Genetics, 18, 321–356. https://doi.org/10.1146/ann...
Weyrich, L. S., Duchene, S., Soubrier, J., Arriola, L., Llamas, B., Breen, J., … Cooper, A. (2017). Neanderthal behaviour, diet, and disease inferred from ancient DNA in dental calculus. Nature, 544(7650), 357–361. https://doi.org/10.1038/nat...