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:

Early Pleistocene enamel proteome from Dmanisi resolves Stephanorhinus phylogeny

Abstract

The sequencing of ancient DNA has enabled the reconstruction of speciation, migration and admixture events for extinct taxa1. However, the irreversible post-mortem degradation2 of ancient DNA has so far limited its recovery—outside permafrost areas—to specimens that are not older than approximately 0.5 million years (Myr)3. By contrast, tandem mass spectrometry has enabled the sequencing of approximately 1.5-Myr-old collagen type I4, and suggested the presence of protein residues in fossils of the Cretaceous period5—although with limited phylogenetic use6. In the absence of molecular evidence, the speciation of several extinct species of the Early and Middle Pleistocene epoch remains contentious. Here we address the phylogenetic relationships of the Eurasian Rhinocerotidae of the Pleistocene epoch7,8,9, using the proteome of dental enamel from a Stephanorhinus tooth that is approximately 1.77-Myr old, recovered from the archaeological site of Dmanisi (South Caucasus, Georgia)10. Molecular phylogenetic analyses place this Stephanorhinus as a sister group to the clade formed by the woolly rhinoceros (Coelodonta antiquitatis) and Merck’s rhinoceros (Stephanorhinus kirchbergensis). We show that Coelodonta evolved from an early Stephanorhinus lineage, and that this latter genus includes at least two distinct evolutionary lines. The genus Stephanorhinus is therefore currently paraphyletic, and its systematic revision is needed. We demonstrate that sequencing the proteome of Early Pleistocene dental enamel overcomes the limitations of phylogenetic inference based on ancient collagen or DNA. Our approach also provides additional information about the sex and taxonomic assignment of other specimens from Dmanisi. Our findings reveal that proteomic investigation of ancient dental enamel—which is the hardest tissue in vertebrates11, and is highly abundant in the fossil record—can push the reconstruction of molecular evolution further back into the Early Pleistocene epoch, beyond the currently known limits of ancient DNA preservation.

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

Access options

Buy this article

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

Fig. 1: Location of Dmanisi, stratigraphy, and specimen Dm.5/157–16635.
Fig. 2: Degradation of the enamel proteome.
Fig. 3: Sequence motif analysis of phosphorylation sites in the proteome of ancient enamel.
Fig. 4: Phylogenetic relationships between the comparative dataset of enamel proteomes and specimen Dm.5/157–16635.

Similar content being viewed by others

Data availability

All of the mass spectrometry proteomics data have been deposited in the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository with the dataset identifier PXD011008. Genomic BAM files used for Rhinocerotidae protein sequence translation and protein sequence alignments used for phylogenetic reconstruction are available on Figshare (https://doi.org/10.6084/m9.figshare.7212746).

Code availability

The in-house R script used to align the peptide sequences confidently identified by the PEAKS searches is available to everyone upon request to the corresponding authors.

References

  1. Cappellini, E. et al. Ancient biomolecules and evolutionary inference. Annu. Rev. Biochem. 87, 1029–1060 (2018).

    Article  CAS  PubMed  Google Scholar 

  2. Dabney, J., Meyer, M. & Pääbo, S. Ancient DNA damage. Cold Spring Harb. Perspect. Biol. 5, a012567 (2013).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Meyer, M. et al. Nuclear DNA sequences from the Middle Pleistocene Sima de los Huesos hominins. Nature 531, 504–507 (2016).

    Article  CAS  ADS  PubMed  Google Scholar 

  4. Wadsworth, C. & Buckley, M. Proteome degradation in fossils: investigating the longevity of protein survival in ancient bone. Rapid Commun. Mass Spectrom. 28, 605–615 (2014).

    Article  CAS  ADS  PubMed  PubMed Central  Google Scholar 

  5. Schweitzer, M. H. et al. Analyses of soft tissue from Tyrannosaurus rex suggest the presence of protein. Science 316, 277–280 (2007).

    Article  CAS  ADS  PubMed  Google Scholar 

  6. Schroeter, E. R. et al. Expansion for the Brachylophosaurus canadensis collagen I sequence and additional evidence of the preservation of Cretaceous protein. J. Proteome Res. 16, 920–932 (2017).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Willerslev, E. et al. Analysis of complete mitochondrial genomes from extinct and extant rhinoceroses reveals lack of phylogenetic resolution. BMC Evol. Biol. 9, 95 (2009).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Welker, F. et al. Middle Pleistocene protein sequences from the rhinoceros genus Stephanorhinus and the phylogeny of extant and extinct Middle/Late Pleistocene Rhinocerotidae. PeerJ 5, e3033 (2017).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Kirillova, I. et al. Discovery of the skull of Stephanorhinus kirchbergensis (Jäger, 1839) above the Arctic Circle. Quat. Res. 88, 537–550 (2017).

    Article  CAS  Google Scholar 

  10. Lordkipanidze, D. et al. A complete skull from Dmanisi, Georgia, and the evolutionary biology of early Homo. Science 342, 326–331 (2013).

    Article  CAS  ADS  PubMed  Google Scholar 

  11. Eastoe, J. E. Organic matrix of tooth enamel. Nature 187, 411–412 (1960).

    Article  CAS  ADS  PubMed  Google Scholar 

  12. Orlando, L. et al. Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse. Nature 499, 74–78 (2013).

    Article  CAS  ADS  PubMed  Google Scholar 

  13. Demarchi, B. et al. Protein sequences bound to mineral surfaces persist into deep time. eLife 5, e17092 (2016).

    Article  PubMed  PubMed Central  Google Scholar 

  14. Welker, F. et al. Ancient proteins resolve the evolutionary history of Darwin’s South American ungulates. Nature 522, 81–84 (2015).

    Article  CAS  ADS  PubMed  Google Scholar 

  15. Chen, F. et al. A late Middle Pleistocene Denisovan mandible from the Tibetan Plateau. Nature 569, 409–412 (2019).

    Article  CAS  ADS  PubMed  Google Scholar 

  16. Nei, M. Molecular Evolutionary Genetics Vol. 75, 39–63 (Columbia Univ. Press, 1987).

  17. Buckley, M., Warwood, S., van Dongen, B., Kitchener, A. C. & Manning, P. L. A fossil protein chimera; difficulties in discriminating dinosaur peptide sequences from modern cross-contamination. Proc. R. Soc. Lond. B 284, 20170544 (2017).

    Article  CAS  Google Scholar 

  18. Gabunia, L. et al. Earliest Pleistocene hominid cranial remains from Dmanisi, Republic of Georgia: taxonomy, geological setting, and age. Science 288, 1019–1025 (2000).

    Google Scholar 

  19. Ferring, R. et al. Earliest human occupations at Dmanisi (Georgian Caucasus) dated to 1.85–1.78 Ma. Proc. Natl Acad. Sci. USA 108, 10432–10436 (2011).

    Article  CAS  ADS  PubMed  PubMed Central  Google Scholar 

  20. Castiblanco, G. A. et al. Identification of proteins from human permanent erupted enamel. Eur. J. Oral Sci. 123, 390–395 (2015).

    Article  CAS  PubMed  Google Scholar 

  21. Stewart, N. A. et al. The identification of peptides by nanoLC-MS/MS from human surface tooth enamel following a simple acid etch extraction. RSC Advances 6, 61673–61679 (2016).

    Article  CAS  Google Scholar 

  22. van Doorn, N. L., Wilson, J., Hollund, H., Soressi, M. & Collins, M. J. Site-specific deamidation of glutamine: a new marker of bone collagen deterioration. Rapid Commun. Mass Spectrom. 26, 2319–2327 (2012).

    Article  ADS  CAS  PubMed  Google Scholar 

  23. Catak, S., Monard, G., Aviyente, V. & Ruiz-López, M. F. Computational study on nonenzymatic peptide bond cleavage at asparagine and aspartic acid. J. Phys. Chem. A 112, 8752–8761 (2008).

    Article  CAS  PubMed  Google Scholar 

  24. Hunter, T. Why nature chose phosphate to modify proteins. Phil. Trans. R. Soc. Lond. B 367, 2513–2516 (2012).

    Article  CAS  Google Scholar 

  25. Hu, J. C. C., Yamakoshi, Y., Yamakoshi, F., Krebsbach, P. H. & Simmer, J. P. Proteomics and genetics of dental enamel. Cells Tissues Organs 181, 219–231 (2005).

    Article  CAS  PubMed  Google Scholar 

  26. Tagliabracci, V. S. et al. Secreted kinase phosphorylates extracellular proteins that regulate biomineralization. Science 336, 1150–1153 (2012).

    Article  CAS  ADS  PubMed  PubMed Central  Google Scholar 

  27. Cleland, T. P. Solid digestion of demineralized bone as a method to access potentially insoluble proteins and post-translational modifications. J. Proteome Res. 17, 536–542 (2018).

    Article  CAS  PubMed  Google Scholar 

  28. Antoine, P.-O. et al. A revision of Aceratherium blanfordi Lydekker, 1884 (Mammalia: Rhinocerotidae) from the Early Miocene of Pakistan: postcranials as a key. Zool. J. Linn. Soc. 160, 139–194 (2010).

    Article  Google Scholar 

  29. Steiner, C. C. & Ryder, O. A. Molecular phylogeny and evolution of the Perissodactyla. Zool. J. Linn. Soc. 163, 1289–1303 (2011).

    Article  Google Scholar 

  30. Hobolth, A., Dutheil, J. Y., Hawks, J., Schierup, M. H. & Mailund, T. Incomplete lineage sorting patterns among human, chimpanzee, and orangutan suggest recent orangutan speciation and widespread selection. Genome Res. 21, 349–356 (2011).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Rieseberg, L. H. Evolution: replacing genes and traits through hybridization. Curr. Biol. 19, R119–R122 (2009).

    Article  CAS  PubMed  Google Scholar 

  32. Guérin, C. Les Rhinocéros (Mammalia, Perissodactyla) du Miocène Terminal au Pleistocène Supérieur en Europe occidentale, Comparaison avec les Espèces Actuelles (Documents du Laboratoire de Geologie de la Faculte des Sciences de Lyon, volume 79) (Univ. Claude-Bernard, 1980).

  33. Deng, T. et al. Out of Tibet: Pliocene woolly rhino suggests high-plateau origin of Ice Age megaherbivores. Science 333, 1285–1288 (2011).

    Article  CAS  ADS  PubMed  Google Scholar 

  34. Orlando, L. et al. Ancient DNA analysis reveals woolly rhino evolutionary relationships. Mol. Phylogenet. Evol. 28, 485–499 (2003).

    Article  CAS  PubMed  Google Scholar 

  35. Yuan, J. et al. Ancient DNA sequences from Coelodonta antiquitatis in China reveal its divergence and phylogeny. Sci. China Earth Sci. 57, 388–396 (2014).

    Article  CAS  Google Scholar 

  36. Penkman, K. E. H., Kaufman, D. S., Maddy, D. & Collins, M. J. Closed-system behaviour of the intra-crystalline fraction of amino acids in mollusc shells. Quat. Geochronol. 3, 2–25 (2008).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Hendy, J. et al. A guide to ancient protein studies. Nat. Ecol. Evol. 2, 791–799 (2018).

    Article  PubMed  Google Scholar 

  38. Wiśniewski, J. R., Zougman, A., Nagaraj, N. & Mann, M. Universal sample preparation method for proteome analysis. Nat. Methods 6, 359–362 (2009).

    Article  CAS  PubMed  Google Scholar 

  39. Cappellini, E. et al. Resolution of the type material of the Asian elephant, Elephas maximus Linnaeus, 1758 (Proboscidea, Elephantidae). Zool. J. Linn. Soc. 170, 222–232 (2014).

    Google Scholar 

  40. Kulak, N. A., Pichler, G., Paron, I., Nagaraj, N. & Mann, M. Minimal, encapsulated proteomic-sample processing applied to copy-number estimation in eukaryotic cells. Nat. Methods 11, 319–324 (2014).

    Article  CAS  PubMed  Google Scholar 

  41. Mackie, M. et al. Palaeoproteomic profiling of conservation layers on a 14th century Italian wall painting. Angew. Chem. Int. Edn 57, 7369–7374 (2018).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Cappellini, E. et al. Proteomic analysis of a Pleistocene mammoth femur reveals more than one hundred ancient bone proteins. J. Proteome Res. 11, 917–926 (2012).

    Article  CAS  PubMed  Google Scholar 

  43. Cox, J. & Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat. Biotechnol. 26, 1367–1372 (2008).

    Article  CAS  PubMed  Google Scholar 

  44. Zhang, J. et al. PEAKS DB: de novo sequencing assisted database search for sensitive and accurate peptide identification. Mol. Cell. Proteomics 11, M111.010587 (2012).

    Article  CAS  PubMed  Google Scholar 

  45. The UniProt Consortium. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 45, D158–D169 (2017).

    Article  CAS  Google Scholar 

  46. O’Leary, N. A. et al. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 44, D733–D745 (2016).

    Article  CAS  PubMed  Google Scholar 

  47. Welker, F. et al. Palaeoproteomic evidence identifies archaic hominins associated with the Châtelperronian at the Grotte du Renne. Proc. Natl Acad. Sci. USA 113, 11162–11167 (2016).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Kearse, M. et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649 (2012).

    Article  PubMed  PubMed Central  Google Scholar 

  49. Gabriels, R., Martens, L. & Degroeve, S. Updated MS2PIP web server delivers fast and accurate MS2 peak intensity prediction for multiple fragmentation methods, instruments and labeling techniques. Nucleic Acids Res. 47, W295–W299 (2019).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Tyanova, S., Temu, T. & Cox, J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. Nat. Protocols 11, 2301–2319 (2016).

    Article  CAS  PubMed  Google Scholar 

  51. Colaert, N., Helsens, K., Martens, L., Vandekerckhove, J. & Gevaert, K. Improved visualization of protein consensus sequences by iceLogo. Nat. Methods 6, 786–787 (2009).

    Article  CAS  PubMed  Google Scholar 

  52. Korneliussen, T. S., Albrechtsen, A. & Nielsen, R. ANGSD: analysis of next generation sequencing data. BMC Bioinformatics 15, 356 (2014).

    Article  PubMed  PubMed Central  Google Scholar 

  53. Briggs, A. W. et al. Removal of deaminated cytosines and detection of in vivo methylation in ancient DNA. Nucleic Acids Res. 38, e87 (2010).

    Article  CAS  PubMed  Google Scholar 

  54. Altschul, S. F. et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402 (1997).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Sea Urchin Genome Sequencing Consortium. The genome of the sea urchin Strongylocentrotus purpuratus. Science 314, 941–952 (2006).

    Article  ADS  PubMed Central  Google Scholar 

  56. Katoh, K. & Frith, M. C. Adding unaligned sequences into an existing alignment using MAFFT and LAST. Bioinformatics 28, 3144–3146 (2012).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Schliep, K. P. phangorn: phylogenetic analysis in R. Bioinformatics 27, 592–593 (2011).

    Article  CAS  PubMed  Google Scholar 

  58. Guindon, S. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321 (2010).

    Article  CAS  PubMed  Google Scholar 

  59. Ronquist, F. et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539–542 (2012).

    Article  PubMed  PubMed Central  Google Scholar 

  60. Rohland, N. & Hofreiter, M. Comparison and optimization of ancient DNA extraction. Biotechniques 42, 343–352 (2007).

    Article  CAS  PubMed  Google Scholar 

  61. Meyer, M. & Kircher, M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb. Protoc. 2010, pdb.prot5448 (2010).

    Article  PubMed  Google Scholar 

  62. Schubert, M. et al. Characterization of ancient and modern genomes by SNP detection and phylogenomic and metagenomic analysis using PALEOMIX. Nat. Protocols 9, 1056–1082 (2014).

    Article  CAS  PubMed  Google Scholar 

  63. Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760 (2009).

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Dickinson, M. R., Lister, A. M. & Penkman, K. E. H. A new method for enamel amino acid racemization dating: a closed system approach. Quat. Geochronol. 50, 29–46 (2019).

    Article  Google Scholar 

Download references

Acknowledgements

E.C. and F.W. are supported by the VILLUM FONDEN (grant number 17649) and by the European Commission through a Marie Skłodowska Curie (MSC) Individual Fellowship (grant number 795569). E.W. is supported by the Lundbeck Foundation, the Danish National Research Foundation, the Novo Nordisk Foundation, the Carlsberg Foundation, KU2016 and the Wellcome Trust. E.C., C.K., J.V.O., P.R. and D.S. are supported by the European Commission through the MSC European Training Network ‘TEMPERA’ (grant number 722606). M.M. and R.R.J.-C. are supported by the University of Copenhagen KU2016 (UCPH Excellence Programme) grant. M.M. is also supported by the Danish National Research Foundation award PROTEIOS (DNRF128). Work at the Novo Nordisk Foundation Center for Protein Research is funded in part by a donation from the Novo Nordisk Foundation (grant number NNF14CC0001). M.R.D. is supported by a PhD DTA studentship from NERC and the Natural History Museum (NE/K500987/1 & NE/L501761/1). K.P. is supported by the Leverhulme Trust (PLP -2012-116). L.R. and L.P. are supported by the Italian Ministry for Foreign Affairs (MAECI, DGSP-VI). L.P. was also supported by the EU-SYNTHESYS project (AT-TAF-2550, DE-TAF-3049, GB-TAF-2825, HU-TAF-3593 and ES-TAF-2997) funded by the European Commission. L.D. is supported by the Swedish Research Council (grant number 2017-04647) and FORMAS (grant number 2015-676). M.T.P.G. is supported by ERC Consolidator Grant ‘Extinction genomics’ (grant number 681396). L.O. is supported by the ERC Consolidator Grant ‘PEGASUS’ (grant agreement number 681605). B.S., J.K. and P.D.H. are supported by the Gordon and Betty Moore foundation. B.M.-N. is supported by the Spanish Ministry of Sciences (grant number CGL2016-80975-P) and the Generalitat de Catalunya, Spain (grant number 2017SGR 859). J.A. is supported by the Spanish Ministry of Sciences (grant number CGL2016-80000-P). R.F. is supported by National Science Foundation (grant number 1025245). The ancient DNA analysis was carried out using the facilities of the University of Luxembourg, the Swedish Museum of Natural History and UC Santa Cruz. We acknowledge support from the Science for Life Laboratory, the National Genomics Infrastructure (Sweden) and UPPMAX for providing assistance with massive parallel sequencing and computational infrastructure. Research at Dmanisi is supported by the John Templeton Foundation (grant number 52935), and the Shota Rustaveli Science Foundation (grant number 18-27262). We thank B. Triozzi and K. Murphy Gregersen for technical support.

Author information

Authors and Affiliations

Authors

Contributions

E.C., D. Lordkipanidze and E.W. designed the study. A.K.F., M.M., R.R.J.-C., M.E.A., M.R.D., K.P. and E.C. performed laboratory experiments. M.B., M.T., R.F., E.P., T.W.S. Jr, Y.L.C., A. Götherström, S.K.S.S.N., P.D.H., J.D.K., I.K., Y.M., J.A., R.-D.K., G.K., B.M.-N., M.-H.S.S., S.L., M.S.V., B.S., L.D., M.T.P.G. and D. Lordkipanidze provided ancient samples or modern reference material. E.C., F.W., L.P., J.R.-M., D. Lyon, J.V.M.-M., D.S., C.D.K., A. Ginolhac, L.O., L.R., J.V.O., P.L.R., M.R.D. and K.P. performed analyses and data interpretation. E.C., F.W., J.R.-M., L.P. and E.W. wrote the manuscript with contributions from all authors.

Corresponding authors

Correspondence to Enrico Cappellini, Jesper V. Olsen or Eske Willerslev.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Peer review information Nature thanks Benedikt Kessler, Tina Warinner and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Extended data figures and tables

Extended Data Fig. 1 Generalized stratigraphic profiles for Dmanisi, indicating origins of the specimens.

a, Type section of the Dmanisi M5 excavation block. b, Stratigraphic profile of excavation area M6. M6 preserves a larger gully associated with the pipe-gully phase of stratigraphic–geomorphic development in stratum B1. The thickness of the stratum B1 gully fill extends to the basalt surface but includes ‘rip-ups’ of strata A1 and A2, showing that the deposits in stratum B1 post-date those of stratum A. c, Stratigraphic section of excavation area M17. Here, Stratum B1 was deposited after the erosion of stratum A deposits. The stratigraphic position of specimen Dm.5/157–16635 is highlighted with a red diamond. The Masavara basalt is about 50 cm below the base of the profile shown. d, Northern section of block 2. Following the collapse of a pipe and erosion to the basalt, the deeper part of this area was filled with local gully fill of strata B1x, B1y and B1z. Note the uniform burial of all stratum B1 deposits by strata B2, B3 and B4. The sampled specimens are indicated by the five-digit CGG numbers. Extended Data Table 1 provides both the CGG and GNM specimen numbers.

Extended Data Fig. 2 Proteome-sequence coverage for specimen Dm.5/157-16635.

a, c, e, g, i, j, Peptide–spectrum match (PSM) sequence coverage of the proteins AMBN (a), ENAM (c), AMELX (e), AMTN (g), MMP20 (i) and ALB (j). Annotations include ‘amino acid position, amino acid called in that position (number of PSMs and peptides covering that position)’ for the phylogenetically informative single-amino-acid polymorphisms within Rhinocerotidae. b, d, f, h, Frequency (per cent) of phosphorylated (green) and unphosphorylated (red) PSMs per amino acid position for AMBN (b), ENAM (d), AMELX (f) and AMTN (h). Numbers within the bars provide the PSM counts. k, Violin plot of distribution of PSM coverage for all covered sites (n = 693), and for sites of phylogenetic relevance (single-amino-acid polymorphisms, n = 30). The box plots define the range of the data, with whiskers extending to 1.5× interquartile range, boxes denoting the 25th and 75th percentiles and dots indicating the median. All panels are based only on MaxQuant search results. The Supplementary Data contains examples of MS/MS spectra, and fragment-ion series alignments for each of the marked single-amino-acid polymorphisms.

Source data

Extended Data Fig. 3 Peptide and fragment-ion coverage of AMELX isoform 1 and isoform 2 from specimen Dm.M6/7.II.296–16856.

Peptides specific to AMELX isoform 1 and isoform 2 appear in the top and bottom parts of the figure, respectively. No AMELX isoform 2 is currently reported in public databases for the Cervidae group. Accordingly, the AMELX-isoform-2-specific peptides were identified by MaxQuant spectral matching against bovine (Bos taurus) AMELX isoform 2 (UniProt accession number P02817-2). AMELX isoform 2 (also known as leucine-rich amelogenin peptide (LRAP)) is a naturally occurring isoform of AMELX from the translation product of an alternatively spliced transcript.

Extended Data Fig. 4 Amino acid racemization.

Extent of intra-crystalline racemization in enamel for the free amino acid (FAA, x axis) fraction and the total hydrolysable amino acids (THAA, y axis) fraction for four amino acids (Asp plus Asn (here denoted Asx), Glu plus Gln (here denoted Glx), Ala and Phe). Note the differences in axis scale. Intra-crystalline data from Proboscidea enamel from a range of sites in the UK64 have been shown for comparison (grey crosses). Taxa from both Dmanisi and the UK exhibit a similar relationship between FAA and THAA racemization, and R2 values have been calculated on the basis of a polynomial relationship (order = 2, all > 0.93).

Source data

Extended Data Fig. 5 Phosphorylation in the proteome of ancient enamel.

Annotated spectra including phosphorylated (here denoted ph) serine (S). a, Phosphorylation in the S-X-E motif of AMELX. b, Phosphorylation in the S-X-phosphorylated S motif of AMBN. Phosphorylation was independently observed in all three separate analyses of Dm.5/157–16635, including multiple spectra and peptides (Extended Data Fig. 2).

Extended Data Fig. 6 Phylogenetic relationships between the comparative reference dataset and specimen Dm.bXI–16857.

Consensus tree from Bayesian inference. The posterior probability of each bipartition is shown as a percentage to the left of each node.

Extended Data Fig. 7 AMELY-specific matches.

a, Specimen Dm.6/151.4.A4.12–16630. b, Specimen Dm.69/64.3.B1.53–16631. c, Specimen Dm.8/154.4.A4.22–16639. d, Specimen Dm.M6/7.II.296–16856. Note the presence of deamidated glutamine (deQ) and asparagine (deN), oxidated methionine (oxM) and phosphorylated serine (phS).

Extended Data Fig. 8 Effect of the missingness in the tree topology.

a, Maximum-likelihood phylogeny obtained using PhyML and the protein alignment that excludes Dm.5/157–16635. b, Topologies obtained from 100 random replicates of the woolly rhinoceros (Cantiquitatis). In each replicate, the number of missing sites was similar to that observed for the Dm.5/157–16635 specimen (72.4% missingness). The percentage shown for each topology indicates the number of replicates in which that particular topology was recovered. c, As in b, but for the Javan rhinoceros (R. sondaicus). d, As in b, but for the black rhinoceros (D. bicornis).

Extended Data Table 1 Genome and proteome survival in 23 specimens of fossil fauna from Dmanisi
Extended Data Table 2 Proteome composition and coverage

Supplementary information

Supplementary Information

Supplementary Materials, Methods and Results, with Figures and Tables. Detailed description, enriched with figures and tables, of: (i) the studied specimens, (ii) the experimental procedures used to generate the data, and (iii) the results, both positive and negative, supporting the conclusions reported in the main text.

Reporting Summary

Supplementary Data

Selection of automatically and manually annotated tandem MS/MS spectra, retrieved from Dmanisi ~1.77 Myr old specimens as well as synthetic peptides, supporting the identification of: (i) phylogenetically informative amino acid positions observed in specimen Dm.5/157-16635, and (ii) phosphorylated sites.

Source data

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cappellini, E., Welker, F., Pandolfi, L. et al. Early Pleistocene enamel proteome from Dmanisi resolves Stephanorhinus phylogeny. Nature 574, 103–107 (2019). https://doi.org/10.1038/s41586-019-1555-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41586-019-1555-y

This article is cited by

Search

Quick links

Nature Briefing: Translational Research

Sign up for the Nature Briefing: Translational Research newsletter — top stories in biotechnology, drug discovery and pharma.

Get what matters in translational research, free to your inbox weekly. Sign up for Nature Briefing: Translational Research