[go: up one dir, main page]
Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Letter
  • Published:

Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis

Nature volume 514, pages 494–497 (2014)Cite this article

Subjects

This article has been updated

Abstract

Modern strains of Mycobacterium tuberculosis from the Americas are closely related to those from Europe, supporting the assumption that human tuberculosis was introduced post-contact1. This notion, however, is incompatible with archaeological evidence of pre-contact tuberculosis in the New World2. Comparative genomics of modern isolates suggests that M. tuberculosis attained its worldwide distribution following human dispersals out of Africa during the Pleistocene epoch3, although this has yet to be confirmed with ancient calibration points. Here we present three 1,000-year-old mycobacterial genomes from Peruvian human skeletons, revealing that a member of the M. tuberculosis complex caused human disease before contact. The ancient strains are distinct from known human-adapted forms and are most closely related to those adapted to seals and sea lions. Two independent dating approaches suggest a most recent common ancestor for the M. tuberculosis complex less than 6,000 years ago, which supports a Holocene dispersal of the disease. Our results implicate sea mammals as having played a role in transmitting the disease to humans across the ocean.

This is a preview of subscription content, access via your institution

Access options

Buy this article

  • Purchase on SpringerLink
  • Instant access to full article PDF
Buy now

Prices may be subject to local taxes which are calculated during checkout

Figure 1: Archaeological description of the skeletal samples.
Figure 2: Coverage plots for three ancient genomes.
Figure 3: Phylogenetic analysis.

Similar content being viewed by others

Accession codes

Primary accessions

Sequence Read Archive

Data deposits

Raw sequencing data have been deposited in the National Center for Biotechnology Information Sequence Read Archive under accession numbers SRP041177 for the ancient Peruvian samples and SRP041181 for the M. pinnipedii strains.

Change history

  • 23 October 2014

    Minor changes were made to the author list, Fig. 3 and ED Figs 1 and 8.

References

  1. Hershberg, R. et al. High functional diversity in Mycobacterium tuberculosis driven by genetic drift and human demography. PLoS Biol. 6, e311 (2008)

    Article  Google Scholar 

  2. Roberts, C. A. & Buikstra, J. E. The Bioarchaeology of Tuberculosis. A Global View on a Reemerging Disease 187–213 (Univ. Press of Florida, 2003)

  3. Comas, I. et al. Out-of-Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with modern humans. Nature Genet. 45, 1176–1182 (2013)

    Article  CAS  Google Scholar 

  4. Wirth, T. et al. Origin, spread, and demography of the Mycobacterium tuberculosis complex. PLoS Pathogens 4, e1000160 (2008)

    Article  Google Scholar 

  5. Cockburn, A. The Evolution and Eradication of Infectious Diseases (Johns Hopkins Press, 1963)

    Google Scholar 

  6. Brosch, R. et al. A new evolutionary scenario for the Mycobacterium tuberculosis complex. Proc. Natl Acad. Sci. USA 99, 3684–3689 (2002)

    Article  ADS  CAS  Google Scholar 

  7. Gagneux, S. & Small, P. M. Global phylogeography of Mycobacterium tuberculosis and implications for tuberculosis product development. Lancet Infect. Dis. 7, 328–337 (2007)

    Article  Google Scholar 

  8. Gagneux, S. et al. Variable host-pathogen compatibility in Mycobacterium tuberculosis. Proc. Natl Acad. Sci. USA 103, 2869–2873 (2006)

    Article  ADS  CAS  Google Scholar 

  9. Comas, I. et al. Human T-cell epitopes of Mycobacterium tuberculosis are evolutionarily hyperconserved. Nature Genet. 42, 498–503 (2010)

    Article  CAS  Google Scholar 

  10. Pepperell, C. S. et al. The role of selection in shaping diversity of natural M. tuberculosis populations. PLoS Pathogens 9, e1003543 (2013)

    Article  CAS  Google Scholar 

  11. Shapiro, B. & Gilbert, M. P. T. No proof that typhoid fever caused the Plague of Athens (a reply to Papagrigorakis et al.). Int. J. Infect. Dis. 10, 334–340 (2006)

    Article  Google Scholar 

  12. Bos, K. I. et al. A draft genome of Yersinia pestis from victims of the Black Death. Nature 478, 506–510 (2011)

    Article  ADS  CAS  Google Scholar 

  13. Bouwman, A. et al. Genotype of a historic strain of Mycobacterium tuberculosis. Proc. Natl Acad. Sci. USA 109, 18511–18516 (2012)

    Article  ADS  CAS  Google Scholar 

  14. Chan, J. Z.-M. et al. Metagenomic analysis of tuberculosis in a mummy. N. Engl. J. Med. 369, 289–290 (2013)

    Article  CAS  Google Scholar 

  15. Briggs, A. W. et al. Patterns of damage in genomic DNA sequences from a Neandertal. Proc. Natl Acad. Sci. USA 104, 14616–14621 (2007)

    Article  ADS  CAS  Google Scholar 

  16. Coscolla, M. et al. Novel Mycobacterium tuberculosis complex from a wild chimpanzee. Emerg. Infect. Dis. 19, 969–976 (2013)

    Article  CAS  Google Scholar 

  17. Bastida, R. et al. Tuberculosis in a wild subantarctic fur seal from Argentina. J. Wildl. Dis. 35, 796–798 (1999)

    Article  CAS  Google Scholar 

  18. Comas, I. et al. Whole-genome sequencing of rifampicin-resistant Mycobacterium tuberculosis strains identifies compensatory mutations in RNA polymerase genes. Nature Genet. 18, 106–110 (2012)

    Article  Google Scholar 

  19. Botella, H. et al. Metallobiology of host-pathogen interactions: an intoxicating new insight. Trends Microbiol. 20, 106–112 (2012)

    Article  CAS  Google Scholar 

  20. Bryant, J. M. et al. Inferring patient to patient transmission of Mycobacterium tuberculosis from whole genome sequencing data. BMC Infect. Dis. 13, 110 (2013)

    Article  Google Scholar 

  21. Elias, S. A., Short, S. K., Nelson, C. H. & Birks, H. H. Life and times of the Bering land bridge. Nature 382, 60–63 (1996)

    Article  ADS  CAS  Google Scholar 

  22. Patrucco, R. et al. Parasitological studies of coprolites of pre-Hispanic Peruvian populations. Curr. Anthropol. 24, 393–394 (1983)

    Article  Google Scholar 

  23. Bastida, R. et al. La tuberculosis en grupos de cazadores recolectores de Patagonia y Tierra del Fuego: nuevas alternativas de contagio a través de la fauna silvestre. Rev. Arg. Antropol. Biol. 13, 83–95 (2011)

    Google Scholar 

  24. Salo, W. et al. Identification of Mycobacterium tuberculosis DNA in a pre-Columbian Peruvian mummy. Proc. Natl Acad. Sci. USA 91, 2091–2094 (1994)

    Article  ADS  CAS  Google Scholar 

  25. Arriaza, R. T. et al. Pre-Columbian tuberculosis in northern Chile: molecular and skeletal evidence. Am. J. Phys. Anthropol. 98, 37–45 (1995)

    Article  CAS  Google Scholar 

  26. Anawalt, P. R. Cultural contacts between Ecuador, west Mexico, and the American Southwest: clothing similarities. Lat. Am. Antiq. 3, 114–129 (1992)

    Article  Google Scholar 

  27. Herring, D. A. & Sattenspiel, L. Social contexts, syndemics, and infectious disease in northern Aboriginal populations. Am. J. Hum. Biol. 19, 190–202 (2007)

    Article  Google Scholar 

  28. Jurczynski, K. et al. Pinniped tuberculosis in Malayan tapirs (Tapirus indicus) and its transmission to other terrestrial mammals. J. Zoo Wildl. Med. 42, 222–227 (2011)

    Article  Google Scholar 

  29. Berg, S. et al. The burden of mycobacterial disease in Ethiopian cattle: implications for public health. PLoS ONE 4, e5068 (2009)

    Article  ADS  Google Scholar 

  30. Cui, Y. et al. Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis. Proc. Natl Acad. Sci. USA 110, 577–582 (2012)

    Article  ADS  Google Scholar 

  31. Drummond, A. J. & Rambaut, A. (2007). BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 7, 214

    Article  Google Scholar 

Download references

Acknowledgements

We thank the following people and institutions for assistance and/or permission for sampling: Museo Contisuyo, Centro Mallqui, S. Guillen, Instituto Nacional de Cultura, Peru, G. Cock, C. Gaither, M. Murphy, M. C. Lozada, S. Burgess, D. Blom, B. Owen, A. Oquiche Hernani, P. Palacios Filinich, S. Williams, B. Vargas, D. Rice, and H. Klaus, S. Pfieffer, the University of Toronto, the Upper Mississippi Valley Archaeological Research Foundation, L. Conrad, Indiana University, G. Milner, the Pennsylvania State University, the American Museum of Natural History, A. Stodder, B. Brier, I. Tattersall, K. Mowbray, the National Museum of Natural History (Smithsonian Institution), B. Billeck, the Smithsonian Institution, N. Tuross, L.-A. Pfister, Rochester Museum & Science Center, L. P. Saunders, C. Grivas, G. Housman, and M. Nieves-Colon. We thank H. Poinar for discussions about capture regions for M. tuberculosis screening. We thank B. Coombes and B. Krause-Kyora for providing modern tuberculosis for bait manufacture. The Huron Wendat Nation is aware of the sampling of Uxbridge bone and is a recipient of information from this study. We acknowledge the following sources of funding: European Research Council starting grant APGREID (to J.K.), the National Science Foundation (to A.C.S. and J.E.B.) for NSF BCS-1063939, NSF-REU BCS-0612222, and NSF BCS-0612222, the George E. Burch Fellow in Theoretic Medicine and Affiliated Sciences at the Smithsonian Institution (2003–2007, to J.E.B.), Social Sciences and Humanities Research Council of Canada postdoctoral fellowship grant 756-2011-501 (to K.I.B.), National Science Foundation Graduate Research Fellowship DGE-1311230 and Jacob K. Javits Fellowship (to K.M.H.), Ramón y Cajal Spanish research grant RYC-2012-10627 (to I.C.), Swiss National Science Foundation PP0033_119205 (to S.G.), National Institutes of Health AI090928 (to S.G.), European Research Council 309540 (to S.G.), PICT0575 Argentina (to R.A.G.), Wadsworth Fellowship from the Wenner-Gren Foundation (to T.J.C.), Wellcome Trust 098051 (to J.P., J.M.B. and S.R.H.), and funding from the Medical Research Council (to J.M.B.).

Author information

Author notes

  1. Kirsten I. Bos, Kelly M. Harkins, Alexander Herbig and Mireia Coscolla: These authors contributed equally to this work.

Authors and Affiliations

  1. Department of Archaeological Sciences, University of Tübingen, Ruemelinstraße 23, 72070 Tübingen, Germany,

    Kirsten I. Bos, Alexander Herbig, Stephen A. Forrest, Verena J. Schuenemann, Kerttu Majander & Johannes Krause

  2. School of Human Evolution and Social Change, Arizona State University, PO Box 872402, Tempe, Arizona 85287-2402, USA,

    Kelly M. Harkins, Alicia K. Wilbur, Jane E. Buikstra & Anne C. Stone

  3. Center for Bioinformatics, University of Tübingen, Sand 14, 72076 Tübingen, Germany,

    Alexander Herbig, Nico Weber, Daniel Huson & Kay Nieselt

  4. Department of Medical Parasitology and Infection Biology, Swiss Tropical and Public Health Institute, Socinstrasse 57, 4002 Basel, Switzerland,

    Mireia Coscolla & Sebastien Gagneux

  5. University of Basel, Petersplatz 1, CH-4003 Basel, Switzerland,

    Mireia Coscolla & Sebastien Gagneux

  6. Genomics and Health Unit, FISABIO-Public Health, Avenida Cataluña 21, 46020 Valencia, Spain,

    Iñaki Comas

  7. CIBER (Centros de Investigación Biomédica en Red) in Epidemiology and Public Health, Instituto de Salud Carlos III, C/ Monforte de Lemos 3-5, Pabellón 11, Planta 0, 28029 Madrid, Spain,

    Iñaki Comas

  8. Pathogen Genomics, The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SA, UK,

    Josephine M. Bryant, Simon R. Harris & Julian Parkhill

  9. Department of Archaeology, University of Cape Town, Private Bag X1, Rondebosch, 7701, South Africa,

    Tessa J. Campbell

  10. Departamento de Biología (FCEyN, CONICET, Laboratorio de Ecología Evolutiva Humana (FACSO, UNCPBA), UNMDP), Calle 508 No. 881 (7631), Quequen, Argentina,

    Ricardo A. Guichon

  11. Department of Anthropology, University of Tennessee, 250 South Stadium Hall, Knoxville, Tennessee 37996, USA,

    Dawnie L. Wolfe Steadman

  12. Department of Anthropology, Indiana University, 701 East Kirkwood Avenue, Bloomington, Indiana 47405-7100, USA,

    Della Collins Cook

  13. Molecular Mycobacteriology, Forschungszentrum Borstel, Parkallee 1, 23845 Borstel, Germany,

    Stefan Niemann

  14. German Center for Infection Research, Forschungszentrum Borstel, Parkallee 1, 23845 Borstel, Germany,

    Stefan Niemann

  15. McGill International TB Centre, McGill University, 1650 Cedar Avenue, Montreal H3G 1A4, Canada,

    Marcel A. Behr

  16. Biotechnology Institute, CICVyA-INTA Castelar, Dr. Nicolás Repetto y De Los Reseros S/N, (B1686IGC) Hurlingham, Buenos Aires, Argentina,

    Martin Zumarraga

  17. Instituto de Investigaciones Marinas y Costeras (CONICET-UNMdP), Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Mar del Plata, San Luis 1722, Mar del Plata 7600, Argentina,

    Ricardo Bastida

  18. Department of Medicine, Imperial College, London W2 1PG, UK

    Douglas Young

  19. Division of Mycobacterial Research, MRC National Institute for Medical Research, Mill Hill, London NW7 1AA, UK,

    Douglas Young

  20. Senckenberg Centre for Human Evolution and Palaeoenvironment, University of Tübingen, Tübingen 72070, Germany,

    Johannes Krause

  21. Max Planck Institute for Science and History, Khalaische Straße 10, 07745 Jena, Germany,

    Johannes Krause

Authors
  1. Kirsten I. Bos
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Kelly M. Harkins
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Alexander Herbig
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Mireia Coscolla
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Nico Weber
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Iñaki Comas
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Stephen A. Forrest
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Josephine M. Bryant
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Simon R. Harris
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Verena J. Schuenemann
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Tessa J. Campbell
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Kerttu Majander
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Alicia K. Wilbur
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Ricardo A. Guichon
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Dawnie L. Wolfe Steadman
    View author publications

    You can also search for this author in PubMed Google Scholar

  16. Della Collins Cook
    View author publications

    You can also search for this author in PubMed Google Scholar

  17. Stefan Niemann
    View author publications

    You can also search for this author in PubMed Google Scholar

  18. Marcel A. Behr
    View author publications

    You can also search for this author in PubMed Google Scholar

  19. Martin Zumarraga
    View author publications

    You can also search for this author in PubMed Google Scholar

  20. Ricardo Bastida
    View author publications

    You can also search for this author in PubMed Google Scholar

  21. Daniel Huson
    View author publications

    You can also search for this author in PubMed Google Scholar

  22. Kay Nieselt
    View author publications

    You can also search for this author in PubMed Google Scholar

  23. Douglas Young
    View author publications

    You can also search for this author in PubMed Google Scholar

  24. Julian Parkhill
    View author publications

    You can also search for this author in PubMed Google Scholar

  25. Jane E. Buikstra
    View author publications

    You can also search for this author in PubMed Google Scholar

  26. Sebastien Gagneux
    View author publications

    You can also search for this author in PubMed Google Scholar

  27. Anne C. Stone
    View author publications

    You can also search for this author in PubMed Google Scholar

  28. Johannes Krause
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

A.C.S., J.E.B., J.K., K.I.B., and K.M.H. conceived the investigation. J.K., K.I.B., A.C.S., S.A.F., N.W., and A.K.W. designed experiments. J.P., R.A.G., D.L.W.S., D.C.C., S.N., M.A.B., M.Z., and R.B. provided samples for analysis. K.I.B., K.M.H., V.J.S., T.J.C., and A.K.W. performed laboratory work. A.H., J.K., S.G., M.C., N.W., K.I.B., I.C., D.Y., J.P., J.M.B., S.R.H., D.H., K.N., A.C.S., K.M.H., J.E.B., T.J.C., D.C.C., and D.L.W.S. performed analyses. K.I.B. wrote the manuscript with contributions from all co-authors.

Corresponding authors

Correspondence to Kirsten I. Bos, Anne C. Stone or Johannes Krause.

Ethics declarations

Competing interests

The authors declare no competing financial interests.

Extended data figures and tables

Extended Data Figure 1 Coverage and damage plots for the M. tuberculosis capture regions for samples 54, 58, and 64.

Extended Data Figure 2 Histograms of SNP allele frequency distributions for the ancient samples and the Hungarian mummy sample using standard mapping parameters.

The x axis denotes the frequency of reads covering a SNP position in which the SNP was detected. The y axis denotes the number of observed SNP calls with the respective frequency. All variants with a SNP allele frequency below 90% are shown.

Extended Data Figure 3 Histograms of SNP allele frequency distributions for the ancient samples, the Hungarian mummy sample, and two modern isolates using stricter mapping and filtering parameters.

The x axis denotes the frequency of reads covering a SNP position in which the SNP was detected. The y axis denotes the number of observed SNP calls with the respective frequency. All variants with a SNP allele frequency below 90% are shown.

Extended Data Figure 4 Maximum parsimony analysis.

a, Maximum parsimony tree of all 262 samples of the complete data set. Positions with missing data were excluded. b, Subtree of the full maximum parsimony tree showing the lineage 6 and animal strains. Positions with missing data were excluded. Branches are labelled with the absolute number of substitutions. Internal nodes are labelled with bootstrap statistics obtained from 1,000 replicates.

Extended Data Figure 5 Maximum likelihood analysis.

a, Maximum likelihood tree of all 262 samples of the complete data set. Positions with missing data were excluded. b, Maximum likelihood subtree showing the lineage 6 and animal strains. Positions with missing data were excluded. Internal nodes are labelled with bootstrap statistics obtained from 200 replicates.

Extended Data Figure 6 Neighbour joining analysis.

a, Neighbour joining tree of all 262 samples of the complete data set. Positions with missing data were excluded. b, Neighbour joining subtree showing the lineage 6 and animal strains. Positions with missing data were excluded. Internal nodes are labelled with bootstrap statistics obtained from 1,000 replicates.

Extended Data Figure 7 Maximum clade credibility tree of M. tuberculosis.

The tree was estimated using the uncorrelated log-normal relaxed clock model in BEAST 1.7.5 (ref. 31). The radiocarbon dates of the ancient Peruvian strains were used as temporal estimates to date the tree. Branch lengths are scaled to years. Branch colours indicate the estimated branch substitution rate on the logarithmic scale shown in the legend at the left.

Extended Data Figure 8

a, Posterior distributions of times to most recent common ancestor (TMRCA) for different MTBC branches, and exponential growth and constant size models. b, Bayesian skyline plot showing estimated effective population sizes for the human lineages. c, Bayesian skyline plot showing estimated effective population sizes for the animal lineages.

Extended Data Figure 9 Maximum likelihood phylogeny of L4 lineage including modern and ancient strains.

The mixed samples are separated out into Hungarian 1 and 2. SNPs were mapped back onto the phylogeny, and branches marked in red are those defined by variants found to be mixed in the Hungarian sample. This allowed us to determine the ancestral nodes and branches for each of the two strains on the tree. The dotted lines represent the unknown length of the terminal branches, with the stars representing the theoretical penultimate node for which age priors were determined.

Extended Data Figure 10 Maximum clade credibility tree produced using BEAST31.

Produced using TreeAnnotator from 9,000 trees. Branch lengths are scaled by age. The mean age (yr bp) of the MRCA plus 95% HPD, and the position of the separated Hungarian ancient strains, are marked on the phylogeny.

Related audio

Supplementary information

Supplementary Information

This file contains Supplementary Methods and archaeological descriptions, Supplementary Tables 2-4, 7, 10-11 and Supplementary References. (PDF 1072 kb)

Supplementary Table 1

Summary table of all samples subject to screening. (XLSX 25 kb)

Supplementary Table 5

A list of all M. tuberculosis strains used for phylogeny and dating. (XLS 53 kb)

Supplementary Table 6

This table contains mapping statistics for the three Peruvian tuberculosis strains. (XLS 31 kb)

Supplementary Table 8

This table contains coverage statistics for all genomes used in analyses. (XLS 25 kb)

Supplementary Table 9

This table contains SNPs identified in the animal cluster. (XLS 95 kb)

PowerPoint slides

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Bos, K., Harkins, K., Herbig, A. et al. Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis. Nature 514, 494–497 (2014). https://doi.org/10.1038/nature13591

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nature13591

This article is cited by

Search

Advanced search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing