Survival of mineral-bound peptides into the Miocene
Abstract
Previously, we showed that authentic peptide sequences could be obtained from 3.8-Ma-old ostrich eggshell (OES) from the site of Laetoli, Tanzania (Demarchi et al., 2016). Here, we show that the same sequences survive in a >6.5 Ma OES recovered from a palaeosteppe setting in northwestern China. The eggshell is thicker than those observed in extant species and consistent with the Liushu Struthio sp. ootaxon. These findings push the preservation of ancient proteins back to the Miocene and highlight their potential for paleontology, paleoecology, and evolutionary biology.
Editor's evaluation
This fundamental study substantially pushes the known preservation of protein sequences bound to mineral surfaces. The successful recovery of these sequences from late Miocene fossil eggshell also has important implications for taxonomic classification. This solid work encourages future paleoproteomic research on paleontological remains from deep antiquity and across various taxa. The paper will be of great interest to a wide range of paleoscientists.
https://doi.org/10.7554/eLife.82849.sa0Introduction
The oldest authenticated peptide sequences to date were reported in 2016 from 3.8-Ma-old ostrich eggshell (OES) from the site of Laetoli, Tanzania (Demarchi et al., 2016). This finding had great scientific impact since it integrated computational chemistry (molecular dynamics simulations) as well as experimental data to propose a mechanism of preservation, concluding that mineral binding ensures the survival of protein sequences. Importantly, this study demonstrated that peptide-bound amino acids could survive into deep time even in hot environments. The effect of temperature on the kinetics of protein diagenesis has been described by several authors, both on the basis of actualistic experiments and of the quantification of the extent of degradation in ancient samples of known ages (Crisp et al., 2013; Demarchi et al., 2013; Hendy et al., 2012; Johnson et al., 1997; Kaufman, 2006; Kaufman, 2003; Schroeder and Bada, 1976; Wehmiller, 2013; Wehmiller, 1977; Wehmiller and Belknap, 1978). This discovery has fuelled the analysis of ancient proteins from other mineral matrices, namely tooth enamel (Cappellini et al., 2019; Welker et al., 2020; Welker et al., 2019) in order to reconstruct phylogenies.
The recovery of peptides has so far been limited to biomineral samples of Plio-Pleistocene age: our attempts at retrieving intact peptides from Cretaceous eggshell were unsuccessful, despite the fact that a genuine intracrystalline fraction of amino acids was preserved in the same sample (Saitta et al., 2020). Here, we show that the same peptide sequences that exhibit strong binding to the calcite surface—and that were recovered from the Laetoli OES—also persist into the late Miocene, in a >6.5-Ma-old eggshell sample from the Linxia Basin, northeastern Tibetan Plateau, China, Liushu Formation.
Results and discussion
The OES (IVPP V26107) was recovered from an area between the towns of Xinji and Songming, in Hezheng County, close to the border with Guanghe County, from mudstone facies of the Liushu Formation (Figure 1).
Both eggshell and skeletal remains are previously known from the Late Miocene Liushu Formation, Linxia Basin, of Gansu Province, where the mean annual temperature is ∼11◦C (Hou et al., 2005; Liu et al., 2016; Li et al., 2021; Wang, 2008). Age control on the highly fossiliferous units of the Lishu Formation is good (Deng et al., 2019; Deng et al., 2013; Zhang et al., 2012), with detailed correlations across China as well as with Neogene deposits globally. These units have been assessed via new magnetostratigraphic data; most of the Hipparion fauna is estimated to be earlier than 6.4–6.5 Ma, with a lower Liushu transition to a Platybelodon fauna yielding age estimates of 11.1–12.5 Ma (Zhang et al., 2012). These dates agree with the prior estimates from biostratigraphic and magnetostratigraphic data (Deng et al., 2019; Deng et al., 2013). The recovery site is closer to the southern part of the basin near the Heilinding section of Zhang et al., 2012, which had an estimated minimum age of 6.4–6.5 Ma (Chron 3An). In contrast, the candidate stratotype section, Guoniguo, that exposes lower parts of the Liushu, is north of the recovery area (Deng et al., 2019; Zhang et al., 2012).
Ostrich (Struthio) remains to have a long history of recovery from the late Miocene of northwest China (Buffetaut and Angst, 2021; Hou et al., 2005; Li et al., 2021; Lowe, 1931; Mikhailov and Zelenkov, 2020). The first reported ostrich fossils from China were from units of the ‘Hipparion Clay’ (Red Clay) dated between 6.54 and 7.18 Ma (Lowe, 1931; Zhu et al., 2008). Mikhailov and Zelenkov, 2020 concluded that the Wang, 2008 Liushu taxon is referable to Struthio and not presently supported as a distinct but related genus; it is proposed to be from a taxon larger than extant Struthio species with a thickness of ~2.4 mm consistent with the eggshell sampled here (Figure 2C).
Given the extreme antiquity of the samples, protein extraction was performed in an ultra-clean facility at the University of Copenhagen in order to minimize any chance of contamination, following Hendy et al., 2018. Furthermore, OES powders were bleached extensively to isolate the intracrystalline fraction only (as described by Demarchi et al., 2016), samples were analysed by LC-MS/MS using a new LC column, and six blanks (analytical and procedural) were included in the run, flanking the sample. Data analysis was equally cautious: raw tandem mass spectrometry data were used to reconstruct potential peptide sequences starting from raw product ion spectra; only de novo peptides with an ALC (Average Local Confidence) score ≥80% were considered. This is the most stringent threshold that can be applied, and it signifies that all potential peptide sequences with lower confidence are discarded. The de novo peptides were searched against the Uniprot/Swissprot database (containing 565,254 manually annotated and reviewed protein sequences). A database of common laboratory contaminants (common Repository of Adventitious Protein [cRAP]) was included in the search. Eleven unique peptide sequences were found to match the sequence of struthiocalcin-1 (Table 1), all containing the typical Asp-rich motif ‘DDDD’ (Figures 3 and 4), which had been shown in our previous paper to be the mineral-binding peptide belonging to the sequence of struthiocalcin-1 (SCA-1). Annotated tandem mass spectra are shown in Figure 3 (ALDDDDYPKG) and Figure 3—figure supplement 1 (ALDDDYPK), Figure 3—figure supplement 2 (SALDDDDYPKG), Figure 3—figure supplement 3 (DDDDYPKGKH), Figure 3—figure supplement 4 (LDDDDYPKGK), Figure 3—figure supplement 5 (SALDDDDYPK), Figure 3—figure supplement 6 (DDDDYPKGK), Figure 3—figure supplement 7 (LDDDDYPKG), Figure 3—figure supplement 8 (DDDYPKGK), Figure 3—figure supplement 9 (DDDDYPK), and Figure 3—figure supplement 10 (DDYPKGK).
Peptide sequences identified in sample IVPP V26107.
Scan | Precursor mass | z | m/z | Peptide sequence | Peptide mass | Error (ppm) | Peptide –10lgP |
---|---|---|---|---|---|---|---|
5401 | 1107.472 | 2 | 555.246 | ALDDDDYPKG | 1107.472 | 5.7 | 27.24 |
5268 | 1050.457 | 2 | 526.236 | ALDDDDYPK | 1050.451 | 6.1 | 24.83 |
5051 | 1194.517 | 2 | 598.266 | SALDDDDYPKG | 1194.504 | 10.9 | 21.11 |
4984 | 1137.494 | 2 | 569.754 | SALDDDDYPK | 1137.483 | 10.1 | 29.94 |
4727 | 1137.494 | 2 | 569.755 | SALDDDDYPK | 1137.483 | 9.7 | 16.76 |
4565 | 1107.483 | 2 | 554.749 | ALDDDDYPKG | 1107.472 | 9.8 | 32.57 |
4565 | 1107.484 | 2 | 554.749 | ALDDDDYPKG | 1107.472 | 10.6 | 28.7 |
4443 | 1050.459 | 2 | 526.237 | ALDDDDYPK | 1050.451 | 7.6 | 17.62 |
4313 | 1050.461 | 2 | 526.235 | ALDDDDYPK | 1050.451 | 9.9 | 15.13 |
4199 | 1050.458 | 2 | 526.236 | ALDDDDYPK | 1050.451 | 7.2 | 16.35 |
4037 | 1036.442 | 2 | 519.230 | LDDDDYPKG | 1036.435 | 7.2 | 16.57 |
2724 | 1164.539 | 2 | 583.277 | LDDDDYPKGK | 1164.53 | 8.2 | 30.91 |
2502 | 866.335 | 2 | 434.176 | DDDDYPK | 866.335 | 6.9 | 18.46 |
2176 | 866.335 | 2 | 434.175 | DDDDYPK | 866.335 | 7.1 | 21.18 |
2026 | 866.338 | 2 | 434.176 | DDDDYPK | 866.329 | 10.1 | 19.07 |
1396 | 1051.454 | 2 | 526.734 | DDDDYPKGK | 1051.446 | 7.8 | 28.62 |
1138 | 936.427 | 2 | 469.221 | DDDYPKGK | 936.419 | 8.4 | 27.8 |
1046 | 1188.513 | 3 | 397.178 | DDDDYPKGKH | 1188.505 | 7.2 | 30.65 |
1042 | 821.398 | 2 | 411.706 | DDYPKGK | 821.392 | 7.6 | 17.39 |

Annotated product ion spectrum of peptide ALDDDDYPKG, m/z 554.749, −10lgP=32.57.
Figure created using http://www.interactivepeptidespectralannotator.com/PeptideAnnotator.html (Brademan et al., 2019).

Coverage of struthiocalcin-1 (SCA-1) in Miocene OES specimen IVPP V26107.
Numbers indicate spectrum matches for each position in the sequence.
Conclusions
We present the first evidence for peptide survival into the Miocene, confirming our previous (Pliocene) data and providing further support to the mechanism of preservation based on the binding of Asx-rich peptides to calcite surfaces (Demarchi et al., 2016). The sequence recovered from the Chinese specimen is identical to the peptides found in the 3.8 Ma Laetoli OES, thus also supporting the attribution of the Liushu ootaxon to genus Struthio. While the four Asp residues are conserved across several avian taxa, in some species, including other ratites, Asp can be substituted by Glu and some of the flanking residues are also variable (Demarchi et al., 2022, fig. S1); therefore, variability within this sequence could be informative for evolutionary relationships between extinct and extant taxa. The absolute age of the Chinese sample is greater than that of the African material (Demarchi et al., 2016), but, given the latitudinal differences between the two sites, it is likely to have experienced lower temperatures throughout its burial history. This suggests that sequences of even greater antiquity may be recovered from biominerals harboring a closed system of proteins, particularly from sites in cold environments (e.g., high altitude and/or latitude).
Materials and methods
Sample preparation and analysis
Request a detailed protocolA subsample of OES specimen IVPP V26107 was prepared in the ultra-clean facility at the University of Copenhagen, following the protocol of Demarchi et al., 2016; Demarchi et al., 2022 and omitting the digestion step. In brief, the fragment was powdered, bleached for 72 hr (NaOCl, 15% w/v) and demineralized in cold weak hydrochloric acid (0.6 M HCl). The acid was added in 200 µl increments until it stopped reacting, for a total volume of 1400 µl 0.6 M HCl. The extracts were exchanged in ammonium bicarbonate buffer (pH=7.8) using 3 kDa MWCO ultrafilters and the peptides were purified and concentrated using C18 Stage Tips (Cappellini et al., 2019; Rappsilber et al., 2007).
Prepared StageTips of the procedural blank and sample were eluted with 30 µl of 40% acetonitrile (ACN) 0.1% formic acid into a 96-well plate prior to LC-MS/MS analysis. To remove the ACN, the plate was vacuum centrifuged until approximately 5 µl remained. Samples were then resuspended with 6 µl of 0.1% trifluoroacetic acid (TFA) 5% ACN. Based on protein concentration results at 205 nm (NanoDrop, Thermo Fisher Scientific), 5 µl of each sample and procedural blank was then separated over a 77 min gradient by an EASY-nLC 1200 (Proxeon, Odense, Denmark) attached to a Q-Exactive HF-X mass spectrometer (Thermo Fisher Scientific, Germany) using a 15 cm column. The column (75 μm inner diameter) was made in-house, laser pulled, and packed with 1.9 μm C18 beads (Dr. Maisch, Germany). Parameters were the same as those already published for historical samples (Mackie et al., 2018). In short, MS1: 120k resolution, maximum injection time (IT) 25 ms, scan target 3E6. MS2: 60k resolution, top 10 mode, maximum IT 118 ms, minimum scan target 3E3, normalized collision energy of 28, dynamic exclusion 20 s, and isolation window of 1.2 m/z. Wash-blanks consisting of 0.1% TFA 5% ACN were also run in order to hinder cross-contamination. The LC-MS/MS run included, in this order: two wash blanks, one procedural blank, one wash blank, the OES sample, and two wash blanks. The data sets have been deposited to the ProteomeXchange Consortium via the Proteomics Identifications Database (PRIDE) partner repository with the identifier PXD035872.
Data analysis
Request a detailed protocolBioinformatic analysis was carried out using PEAKS Studio 8.5 (Bioinformatics Solutions Inc; Zhang et al., 2012). The Uniprot_swissprot database (downloaded 11/08/2021) was used for carrying out the searches and common contaminants were included (cRAP: http://www.thegpm.org/crap/). No enzyme was specified for the digestion and the tolerance was set to 10 ppm on the precursor and 0.05 Da on the fragments. The thresholds for peptide and protein identification were set as follows: peptide score −10lgP≥15, protein score −10lgP≥20, and de novo sequences scores (ALC%)≥80.
Data availability
Tandem mass spectra supporting peptide sequence identification are reported in Figure 3 and Figure 3 - Supplement 1 to 10. Raw mass spectrometry data and results of bioinformatics analysis are available via ProteomeXchange with identifier PXD035872.
-
PRIDEID PXD035872. Ostrich eggshell peptides survival into the Miocene (6-9 Ma).
References
-
Interactive peptide spectral annotator: a versatile web-based tool for proteomic applicationsMolecular & Cellular Proteomics 18:S193–S201.https://doi.org/10.1074/mcp.TIR118.001209
-
BookLate Cenozoic Biostratigraphy of the Linxia Basin, Northwestern China. Fossil Mammals of AsiaColumbia University Press.
-
Neogene integrative stratigraphy and timescale of ChinaScience China Earth Sciences 62:310–323.https://doi.org/10.1007/s11430-017-9155-4
-
Assessing amino acid racemization variability in coral intra-crystalline protein for geochronological applicationsGeochimica Cosmochimica Acta 86:338–353.https://doi.org/10.1016/j.gca.2012.02.020
-
A guide to ancient protein studiesNature Ecology & Evolution 2:791–799.https://doi.org/10.1038/s41559-018-0510-x
-
A Miocene ostrich fossil from Gansu province, northwest ChinaChinese Science Bulletin 50:1808.https://doi.org/10.1360/982005-575
-
The determination of late quaternary paleoenvironments at Equus cave, South Africa, using stable isotopes and amino acid racemization in ostrich eggshellPalaeogeography, Palaeoclimatology, Palaeoecology 136:121–137.https://doi.org/10.1016/S0031-0182(97)00043-6
-
Amino acid paleothermometry of Quaternary ostracodes from the Bonneville basin, UtahQuaternary Science Reviews 22:899–914.https://doi.org/10.1016/S0277-3791(03)00006-4
-
Exceptional preservation of an extinct ostrich from the Late Miocene Linxia basin of ChinaVertebrata PalAsiatica 59:229–244.
-
BookStruthious remains from northern China and Mongolia: with descriptions of Struthio wimani, Struthio anderssoni and Struthio mongolicus spp. novGeological Survey of China.
-
A review of the geochemical applications of the amino acid racemization reactionEarth-Science Reviews 12:347–391.https://doi.org/10.1016/0012-8252(76)90011-8
-
Reexamination of taxonomic assignment of Struthio linxiaensis Hou et al., 2005Acta Palaeontol Sin 47:362–368.
-
Amino acid studies of the Del Mar, California, midden site: apparent rate constants, ground temperature models, and chronological implicationsEarth and Planetary Science Letters 37:184–196.https://doi.org/10.1016/0012-821X(77)90163-7
-
PEAKS DB: de novo sequencing assisted database search for sensitive and accurate peptide identificationMolecular & Cellular Proteomics 11:M111.https://doi.org/10.1074/mcp.M111.010587
-
A new magnetostratigraphic framework for Late Neogene Hipparion Red Clay in the eastern loess plateau of ChinaPalaeogeography, Palaeoclimatology, Palaeoecology 268:47–57.https://doi.org/10.1016/j.palaeo.2008.08.001
Decision letter
-
Yonatan SahleReviewing Editor; University of Cape Town, South Africa
-
George H PerrySenior Editor; Pennsylvania State University, United States
Our editorial process produces two outputs: (i) public reviews designed to be posted alongside the preprint for the benefit of readers; (ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.
Decision letter after peer review:
Thank you for submitting your article "Survival of mineral-bound peptides into the Miocene" for consideration by eLife. Your article has been reviewed by 2 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and George Perry as the Senior Editor. The reviewers have opted to remain anonymous.
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
1) Both reviewers found merit in your work and appreciated the significance of it for the field of paleoproteomics. While your successful identification and recovery of authentic peptide sequences from a late Miocene fossil eggshell pushes the known time range of ancient protein preservation, both reviewers questioned the relevance of the reported thermal age owing to its inherent limitations. Both reviewers highlighted that thermal age provides unreliable estimates, erring on the generous side. Because of this limitation, and the fact that there are more reliable chronometric dates for the analyzed eggshell sample, the relevance, and validity of the reported thermal age is unwarranted. As such, both reviewers recommended that the reported thermal age be omitted from the paper altogether, owing to the inherent limitations of it and the availability of more accurate geological dates for your sample.
2) Reviewer #2 strongly recommended that your MS data be reported, which I agree with.
Both reviewers have provided additional constructive comments, which I encourage you to address in your revision.
Reviewer #1 (Recommendations for the authors):
"The oldest authenticated peptide sequences to date were reported in 2016 from 3.8 Ma old ostrich eggshell (OES) from the site of Laetoli, Tanzania (Demarchi et al., 2016), with a thermal age of the peptides estimated to be equivalent to ~16 Ma at a constant 10°C."
Please remove all reference to thermal age from the entire manuscript. Thermal age isn't real and its usage obscures patterns of preservation/protein breakdown. Also, it artificially inflates the age of samples that are already old. The previously discussed 3.8 Ma ostrich egg shells are actually 3.8 Ma and not 16 Ma. The arbitrariness of thermal age is especially evident in the last paragraph before the acknowledgements considering local temperature across the entire 6.4-6.5 Ma cannot be known and the thermal age calculation is described as 'crude'. Additionally, this sample is actually 2.6-2.7 Ma older than the previously characterized eggshells regardless of whatever the thermal age model outputs.
"25 sequences were found with an ALC (Average Local Confidence) score {greater than or equal to} 80%. This is the most stringent threshold that can be applied" If this was only from the de novo search, please include a table listing the 25 de novo peptides. It may be helpful to indicate if they match the database search as well.
"individual annotated tandem mass spectra are shown in Figure 3 and Figure 3—figure supplement 1, Figure 3—figure supplement 2, Figure 3—figure supplement 3, Figure 3—figure supplement 4, Figure 3—figure supplement 5, Figure 3—figure supplement 6, Figure 3—figure supplement 7, Figure 3—figure supplement 8, Figure 3—figure supplement 9, Figure 3—figure supplement 10, Figure 3—figure supplement 11" Should it be Figure 4? Also, please list each peptide sequence here as well for ease of connecting to the supplemental figures.
Figure 4: Please include a table in the main manuscript with the 12 peptides. I appreciate that they are annotated in the supplement, but there is no reason to not include them. Alternatively, put the aligned peptides in figure 4 in place of the gray circles. I don't know if it is possible with the software used, but please try and make this figure colorblind accessible. Red and green together are not great. Please also include details on how this figure was produced.
"demineralised in just enough cold weak hydrochloric acid" Please quantify how much acid this is.
Reviewer #2 (Recommendations for the authors):
The weakest part of this manuscript is the inclusion of thermal age data, which is a metric I do not find reliable or informative in its own right, but moreover, its inclusion does not seem to add anything of value to the MS data that is presented. The thermal age estimates at the end of this paper constitute using multiple assumptions about ancient temperature and humidity values and their variability (e.g., from a "compilation of proxies") to return what the authors themselves describe as a "very crude" estimate that is nearly identical to the specimen's actual age. For what purpose? Whether "thermal age" as a metric is meant to be a proxy for determining the likelihood of the persistence of peptides (i.e., to assess that a fossil is good for analyses because it is likely less degraded than its age suggests) or to use for authentication of ancient peptides (i.e., to show that very old fossils are in fact 'thermally younger' than their geologic age because they've undergone less heating and therefore the preservation of proteins in them beyond the hypothetical limit of preservation is not problematic) both this paper and the earlier paper giving a thermal age of 16 Ma to the 3.8 Ma Laetoli eggshell show that neither usage is valid. In fact, if thermal age were to be used in this paper as it has been used in the past, it would suggest that the peptides the authors report here represent contamination. Of course, I think that would be absurd, and that the annotated spectra the authors produced in a sterile environment, with no trace of contamination in the blanks, confirm the authenticity of their peptides. However, both things can't be true. If the peptides are real--and I think they are--then thermal age is not a useful metric.
I strongly recommend the authors simply report their MS data, which I think are extremely impactful and worthy of publication on their own.
Technical Comments:
1. It is unclear to me how the authors came to the specific age range of 6.5 Ma to 9 Ma. The text is clear that 6.5 Ma is the minimum age estimate for their recovery site, and they also state that previously reported ostrich fossils were dated between 6.54-7.18 Ma, and that Platybelodon had been dated from lower in the Liushu at 11.1 – 12.5 Ma. So where does the upper bound estimate come from? Normally I wouldn't nitpick this point, but as the authors claim to have the oldest recovered peptides, the rationale for that upper bound is important and should be justified in the text. Otherwise, the authors should drop the range and stick with the minimum estimate of 6.5 Ma.
2. Figure 4 is visually difficult to interpret. I strongly suggest that you change this figure to depict the letters signifying the residues as equally sized (to make them readable) and place proportionately sized bars above them (like a column chart) to show their relative coverage.
3. Materials and methods: What is "just enough" HCl? Please give a volume.
4. It's unclear to me why the authors chose to use a fragment error tolerance of 0.5 Da, given that 0.05 Da or less is standard for ancient protein LC-MS/MS data, and that data from the Q-Exactive is certainly high enough resolution to set a higher tolerance. Additionally, the annotated peptides all appear to be within the 0.05 Da range. Is this a typo? If so, please correct it. If not, I'm curious about how a low fragment error tolerance coupled with a large search space (since it's a no-digest search) impacts the FDR rates. Please report the %FDR represented by a peptide score of −10lgP {greater than or equal to} 30 for the ostrich data.
[Editors' note: further revisions were suggested prior to acceptance, as described below.]
Thank you for resubmitting your work entitled "Survival of mineral-bound peptides into the Miocene" for further consideration by eLife. Your revised article has been evaluated by George Perry (Senior Editor) and a Reviewing Editor.
The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:
I [George Perry] am writing this on behalf of Yonathan Sahle, who has departed for remote, international research, but we corresponded before he left so that I could incorporate his feedback into this decision letter here and get this back to you more efficiently than otherwise.
Specifically, while you have addressed most of the reviewers' concerns, the maintained position regarding the relevance of the thermal age clashes with the strong recommendations by both reviewers (a position the reviewing editor agrees with).
The geological age of the specimens should always be used as the true age, given that the thermal age is sensitive to many factors and has been repeatedly shown to artificially inflate the true age of paleontological remains. In lines 36-43 of the revised manuscript, the thermal age estimate for the Laetoli (Tanzania) specimen is invoked, despite the very large discrepancy of >10 Ma between the true age (3.6 Ma) and the thermal age (~16 Ma) of the Tanzanian sample.
True temperature affects protein preservation and there are several works that have tried to quantify the extent of temperature-induced protein degradation. However, (a) there are several unknowns (estimates) that affect such assessments (e.g., past temperatures) , (b) more samples and contexts are required for a better understanding of these temperature-protein diagenesis relationships, and (c) regardless of the level of confidence of the estimated/modelled thermal age, its relevance for the current paper on ancient proteins from geologically dated OES from Linxia Basin in China is unclear (see also the disconnect in Lines 171-173).
If the thermal age data are omitted altogether, or limited to a very brief and cautious mention in results/discussion, then we would likely be ready to move forward favorably with the paper. But this remaining point is considered essential.
https://doi.org/10.7554/eLife.82849.sa1Author response
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
1) Both reviewers found merit in your work and appreciated the significance of it for the field of paleoproteomics. While your successful identification and recovery of authentic peptide sequences from a late Miocene fossil eggshell pushes the known time range of ancient protein preservation, both reviewers questioned the relevance of the reported thermal age owing to its inherent limitations. Both reviewers highlighted that thermal age provides unreliable estimates, erring on the generous side. Because of this limitation, and the fact that there are more reliable chronometric dates for the analyzed eggshell sample, the relevance, and validity of the reported thermal age is unwarranted. As such, both reviewers recommended that the reported thermal age be omitted from the paper altogether, owing to the inherent limitations of it and the availability of more accurate geological dates for your sample.
We thank the Reviewers and the Editors for their positive assessment of our manuscript. We agree that, in the present case, thermal ages are affected by significant errors (because of the lack of paleotemperature models for this area, and for ages that span all the way back to the Miocene), and therefore only "crude estimates" can be provided. As such, we removed the calculated thermal ages as requested. However, we wish to clarify that thermal ages were not used to provide a chronometric assessment of the samples (neither here nor in our previous paper), but rather as a “diagenesis” index, which allows us to compare the extent of degradation in samples from very different environments. Indeed, protein diagenesis is strongly affected by temperature (Crisp et al., 2013; Demarchi et al., 2013; Hendy et al., 2012; Johnson et al., 1997; Kaufman, 2006, 2003; Schroeder and Bada, 1976; Wehmiller, 2013, 1977; Wehmiller and Belknap, 1978). This effect cannot be ignored in palaeoproteomics – proteins extracted from identical samples from sites in the Tropics vs the Arctic will not show the same extent of preservation even if the samples were of exactly the same age and composition. In order to draw the reader's attention to this important concept we added the following sentence to the Introduction:
“The effect of temperature on the kinetics of protein diagenesis has been described by several authors, both on the basis of actualistic experiments and by quantifying the extent of degradation in ancient samples of known ages (Crisp et al., 2013; Demarchi et al., 2013; Hendy et al., 2012; Johnson et al., 1997; Kaufman, 2006, 2003; Schroeder and Bada, 1976; Wehmiller, 2013, 1977; Wehmiller and Belknap, 1978). By modelling the combined effect of temperature and time on protein diagenesis in a closed system, the peptides recovered from the Laetoli OES were estimated to be equivalent to ~16 Ma at a constant 10°C.”
2) Reviewer #2 strongly recommended that your MS data be reported, which I agree with.
In addition to the tandem mass spectra, we now also present a table with all the identified peptides.
Both reviewers have provided additional constructive comments, which I encourage you to address in your revision.
Reviewer #1 (Recommendations for the authors):
"The oldest authenticated peptide sequences to date were reported in 2016 from 3.8 Ma old ostrich eggshell (OES) from the site of Laetoli, Tanzania (Demarchi et al., 2016), with a thermal age of the peptides estimated to be equivalent to ~16 Ma at a constant 10°C."
Please remove all reference to thermal age from the entire manuscript. Thermal age isn't real and its usage obscures patterns of preservation/protein breakdown. Also, it artificially inflates the age of samples that are already old. The previously discussed 3.8 Ma ostrich egg shells are actually 3.8 Ma and not 16 Ma. The arbitrariness of thermal age is especially evident in the last paragraph before the acknowledgements considering local temperature across the entire 6.4-6.5 Ma cannot be known and the thermal age calculation is described as 'crude'. Additionally, this sample is actually 2.6-2.7 Ma older than the previously characterized eggshells regardless of whatever the thermal age model outputs.
As requested, we have removed our crude estimate of thermal age from the text, as well as the whole final paragraph. However, we respectfully point out that theoretical kinetic models of racemization have been used for half a century to correlate sites of equivalent geological age but with different thermal histories and thus different degrees of protein degradation (e.g. Figure 6 of Wehmiller, J. F. 2013. United States Quaternary coastal sequences and molluscan racemization geochronology – What have they meant for each other over the past 45 years? Quaternary Geochronology, 16, 3–20).
Further, it is demonstrable that when eggshell experiences differences in temperature at a single site (e.g. over glacial/interglacial cycles), this results in corresponding fluctuations in rates of protein degradation (e.g. Figure 9 of Miller, G. H., Fogel, M. L., Magee, J. W., and Gagan, M. K. (2016). Disentangling the impacts of climate and human colonization on the flora and fauna of the Australian arid zone over the past 100 ka using stable isotopes in avian eggshell. Quaternary Science Reviews, 151, 27–57).
Protein degradation in closed systems is influenced by temperature (which thermal age attempts to normalize). Despite the greater antiquity of the samples reported here, our (admittedly crude) calculations suggest that their persistence is in fact not as remarkable as the sequences recovered from Africa. The take home message from the reviewer that “this sample is actually 2.6-2.7 Ma older than the previously characterized eggshells regardless…” therefore minimizes the significance of the African findings, whose thermal age predicted the persistence of sequences of greater antiquity in higher latitudes and consequently even greater antiquity at still higher altitudes and latitudes.
We have added the following sentence to the conclusions: The absolute age of the Chinese sample is greater than that of the African material (Demarchi et al., 2016), but it is estimated to have experienced lower temperatures throughout its burial history. This suggests that sequences of even greater antiquity will be recovered from closed systems buried in even colder sites (e.g. high altitude and/or latitude).
"25 sequences were found with an ALC (Average Local Confidence) score {greater than or equal to} 80%. This is the most stringent threshold that can be applied" If this was only from the de novo search, please include a table listing the 25 de novo peptides. It may be helpful to indicate if they match the database search as well.
We now include a detailed table with all peptide sequences matched to SCA-1, including the score for each peptide, therefore we replaced the statement as follows: “Eleven unique peptide sequences were found to match the sequence of struthiocalcin-1 (Table 1), all containing the typical Asp-rich motif “DDDD” (Figures 3, 4), which had been shown in our previous paper to be the mineral-binding peptide belonging to the sequence of struthiocalcin-1 (SCA-1).”
"individual annotated tandem mass spectra are shown in Figure 3 and Figure 3—figure supplement 1, Figure 3—figure supplement 2, Figure 3—figure supplement 3, Figure 3—figure supplement 4, Figure 3—figure supplement 5, Figure 3—figure supplement 6, Figure 3—figure supplement 7, Figure 3—figure supplement 8, Figure 3—figure supplement 9, Figure 3—figure supplement 10, Figure 3—figure supplement 11" Should it be Figure 4? Also, please list each peptide sequence here as well for ease of connecting to the supplemental figures.
The numbering has been kept as in the original, as figure 3 shows one of the tandem mass spectra. We have now added peptide sequences in parenthesis within the text.
Figure 4: Please include a table in the main manuscript with the 12 peptides. I appreciate that they are annotated in the supplement, but there is no reason to not include them. Alternatively, put the aligned peptides in figure 4 in place of the gray circles.
We have now added a table with the peptide sequences and their scores.
I don't know if it is possible with the software used, but please try and make this figure colorblind accessible. Red and green together are not great. Please also include details on how this figure was produced.
We have redrawn the figure, which now is a simple histogram. We hope this has improved its readability.
"demineralised in just enough cold weak hydrochloric acid" Please quantify how much acid this is.
This information was added: “The acid was added in 200 µl increments until it stopped reacting, for a total volume of 1400 µl 0.6 M HCl.”
Reviewer #2 (Recommendations for the authors):
The weakest part of this manuscript is the inclusion of thermal age data, which is a metric I do not find reliable or informative in its own right, but moreover, its inclusion does not seem to add anything of value to the MS data that is presented. The thermal age estimates at the end of this paper constitute using multiple assumptions about ancient temperature and humidity values and their variability (e.g., from a "compilation of proxies") to return what the authors themselves describe as a "very crude" estimate that is nearly identical to the specimen's actual age. For what purpose? Whether "thermal age" as a metric is meant to be a proxy for determining the likelihood of the persistence of peptides (i.e., to assess that a fossil is good for analyses because it is likely less degraded than its age suggests) or to use for authentication of ancient peptides (i.e., to show that very old fossils are in fact 'thermally younger' than their geologic age because they've undergone less heating and therefore the preservation of proteins in them beyond the hypothetical limit of preservation is not problematic) both this paper and the earlier paper giving a thermal age of 16 Ma to the 3.8 Ma Laetoli eggshell show that neither usage is valid.
The thermal age was included to provide a comparative index; we have not used it to provide numerical estimates of age or as a way to make our data look “better” – in this instance the thermal age actually made our data less impressive because, coming from a colder environment, it is thermally younger than our earlier published findings. As explained above, we do however maintain that the combined effect of temperature and time cannot be ignored, indeed they must be considered. In the absence of better models for temperature variations over time an attempt of normalization, such as that provided by thermal age, provides the only sensible way for comparing samples from widely different geographic regions – as many other authors studying amino acid racemization and other forms of protein diagenesis have convincingly shown in the past.
In fact, if thermal age were to be used in this paper as it has been used in the past, it would suggest that the peptides the authors report here represent contamination. Of course, I think that would be absurd, and that the annotated spectra the authors produced in a sterile environment, with no trace of contamination in the blanks, confirm the authenticity of their peptides. However, both things can't be true. If the peptides are real--and I think they are--then thermal age is not a useful metric.
We thank the reviewer for the positive comments on our data. We however respectfully point out that the reality of the peptide and their “younger” thermal age compared to the African OES are two sides of the same coin and do not contradict each other.
The first side of the coin concerns the authenticity of the peptides recovered: we obtained these from a closed system environment – the intracrystalline fraction within ostrich eggshell; we checked for the absence of contamination. The mechanism of preservation was put forward in 2016 and is supported experimentally and independently here, by the fact that both the Laeotoli and Chines samples yielded the same peptides.
The second side of the coin concerns the way proteins "age": certainly the absolute age of the specimens is crucial, but the local environment is equally important, because diagenesis is a complex network of chemical reactions, which are by definition influenced by temperature. The Chinese samples were preserved in a colder environment than the Laetoli OES, hence the rates of degradation (particularly, peptide bond hydrolysis) must have been slower than the rates of degradation at Laetoli. This is perfectly in line with the fact that peptides were recovered from this older specimen.
I strongly recommend the authors simply report their MS data, which I think are extremely impactful and worthy of publication on their own.
Technical Comments:
1. It is unclear to me how the authors came to the specific age range of 6.5 Ma to 9 Ma. The text is clear that 6.5 Ma is the minimum age estimate for their recovery site, and they also state that previously reported ostrich fossils were dated between 6.54-7.18 Ma, and that Platybelodon had been dated from lower in the Liushu at 11.1 – 12.5 Ma. So where does the upper bound estimate come from? Normally I wouldn't nitpick this point, but as the authors claim to have the oldest recovered peptides, the rationale for that upper bound is important and should be justified in the text. Otherwise, the authors should drop the range and stick with the minimum estimate of 6.5 Ma.
We have changed the range to the minimum age of 6.5 Ma.
2. Figure 4 is visually difficult to interpret. I strongly suggest that you change this figure to depict the letters signifying the residues as equally sized (to make them readable) and place proportionately sized bars above them (like a column chart) to show their relative coverage.
We have modified the figure following the Reviewer’s suggestions
3. Materials and methods: What is "just enough" HCl? Please give a volume.
This has been amended following both reviewers’ comment
4. It's unclear to me why the authors chose to use a fragment error tolerance of 0.5 Da, given that 0.05 Da or less is standard for ancient protein LC-MS/MS data, and that data from the Q-Exactive is certainly high enough resolution to set a higher tolerance. Additionally, the annotated peptides all appear to be within the 0.05 Da range. Is this a typo? If so, please correct it. If not, I'm curious about how a low fragment error tolerance coupled with a large search space (since it's a no-digest search) impacts the FDR rates. Please report the %FDR represented by a peptide score of −10lgP {greater than or equal to} 30 for the ostrich data.
We thank the reviewer for spotting this – it was a typo, which we have corrected. Following the Reviewers’ comments, we now report peptide matches to SCA-1 sequence rather than de novo sequences; the thresholds for peptide and protein identification were set as follows: peptide score −10lgP ≥ 15, protein score −10lgP ≥ 20, de novo sequences scores (ALC%) ≥ 80. While the FDR values are high – due to the low number of peptides, a peptide score of −10lgP ≥ 15 is standard and appears adequate here (see Author response image 1). Most peptides shown in table 1 yielded scores > 20.
References
Crisp M, Demarchi B, Collins M, Morgan-Williams M, Pilgrim E, Penkman K. 2013. Isolation of the intra-crystalline proteins and kinetic studies in Struthio camelus (ostrich) eggshell for amino acid geochronology. Quat Geochronol 16:110–128.
Demarchi B, Collins MJ, Tomiak PJ, Davies BJ, Penkman KEH. 2013. Intra-crystalline protein diagenesis (IcPD) in Patella vulgata. Part II: Breakdown and temperature sensitivity. Quat Geochronol 16:158–172.
Hendy EJ, Tomiak PJ, Collins MJ, Hellstrom J, Tudhope AW, Lough JM, Penkman KEH. 2012.
Assessing amino acid racemization variability in coral intra-crystalline protein for geochronological applications. Geochim Cosmochim Acta 86:338–353.
Johnson BJ, Miller GH, Fogel ML, Beaumont PB. 1997. The determination of late Quaternary paleoenvironments at Equus Cave, South Africa, using stable isotopes and amino acid racemization in ostrich eggshell. Palaeogeogr Palaeoclimatol Palaeoecol 136:121–137.
Kaufman DS. 2006. Temperature sensitivity of aspartic and glutamic acid racemization in the foraminifera Pulleniatina. Quat Geochronol 1:188–207.
Kaufman DS. 2003. Amino acid paleothermometry of Quaternary ostracodes from the Bonneville Basin, Utah. Quat Sci Rev 22:899–914.
Schroeder RA, Bada JL. 1976. A review of the geochemical applications of the amino acid racemization reaction. Earth-Sci Rev 12:347–391.
Wehmiller JF. 2013. United States Quaternary coastal sequences and molluscan racemization geochronology – What have they meant for each other over the past 45 years? Quat Geochronol 16:3–20.
Wehmiller JF. 1977. Amino acid studies of the Del Mar, California, midden site: Apparent rate constants, ground temperature models, and chronological implications. Earth Planet Sci Lett 37:184–196.
Wehmiller JF, Belknap DF. 1978. Alternative Kinetic Models for the Interpretation of Amino Acid Enantiomeric ratios in Pleistocene Mollusks: Examples from California, Washington, and Florida. Quat Res 9:330–348.
[Editors' note: further revisions were suggested prior to acceptance, as described below.]
Specifically, while you have addressed most of the reviewers' concerns, the maintained position regarding the relevance of the thermal age clashes with the strong recommendations by both reviewers (a position the reviewing editor agrees with).
The geological age of the specimens should always be used as the true age, given that the thermal age is sensitive to many factors and has been repeatedly shown to artificially inflate the true age of paleontological remains. In lines 36-43 of the revised manuscript, the thermal age estimate for the Laetoli (Tanzania) specimen is invoked, despite the very large discrepancy of >10 Ma between the true age (3.6 Ma) and the thermal age (~16 Ma) of the Tanzanian sample.
True temperature affects protein preservation and there are several works that have tried to quantify the extent of temperature-induced protein degradation. However, (a) there are several unknowns (estimates) that affect such assessments (e.g., past temperatures) , (b) more samples and contexts are required for a better understanding of these temperature-protein diagenesis relationships, and (c) regardless of the level of confidence of the estimated/modelled thermal age, its relevance for the current paper on ancient proteins from geologically dated OES from Linxia Basin in China is unclear (see also the disconnect in Lines 171-173).
If the thermal age data are omitted altogether, or limited to a very brief and cautious mention in results/discussion, then we would likely be ready to move forward favorably with the paper. But this remaining point is considered essential.
Following the request to remove all mentions of thermal age and/or thermal age estimates, we amended the text as follows:
1) Lines 41-43. This sentence was deleted entirely: “By modelling the combined effect of temperature and time on protein diagenesis in a closed system, the peptides recovered from the Laetoli OES were estimated to be equivalent to ~16 Ma at a constant 10°C”
2) Lines 170-173. The sentence: “The absolute age of the Chinese sample is greater than that of the African material (Demarchi et al., 2016), but it is estimated to have experienced lower temperatures throughout its burial history. This suggests that sequences of even greater antiquity will be recovered from closed systems buried in even colder sites (e.g. high altitude and/or latitude)“ was reworded more cautiously to “The absolute age of the Chinese sample is greater than that of the African material (Demarchi et al., 2016), but, given the latitudinal differences between the two sites, it is likely to have experienced lower temperatures throughout its burial history. This suggests that sequences of even greater antiquity may be recovered from biominerals harbouring a closed system of proteins, particularly from sites in cold environments (e.g. high altitude and/or latitude).”
https://doi.org/10.7554/eLife.82849.sa2Article and author information
Author details
Funding
Ministry of University and Research - Italy (Young Researchers - Rita Levi Montalcini)
- Beatrice Demarchi
Danish National Research Foundation (PROTEIOS (DNRF128))
- Meaghan Mackie
- Matthew J Collins
Alexander von Humboldt Foundation
- Julia Clarke
University of Texas at Austin
- Julia Clarke
National Natural Science Foundation of China
- Zhiheng Li
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
BD’s work was supported by a “Young Researchers - Rita Levi Montalcini” grant from the Italian Ministry of University and Research. MM and MJC were supported by the Danish National Research Foundation Award PROTEIOS (DNRF128). JC acknowledges the Alexander von Humboldt Foundation and the Jackson School of Geosciences and ZL and DT the Chinese National Science Foundation. The authors thank Prof. Jesper Velgaard Olsen at the Novo Nordisk Center for Protein Research for providing access and resources, which were also funded in part by a donation from the Novo Nordisk Foundation (Grant no. NNF14CC0001). For the purpose of open access, MJC is part of the Rights Retention Pilot at the University of Cambridge and has applied a Creative Commons Attribution (CC BY) licence to any Author Accepted Manuscript version arising from this submission. The authors wish to thank the reviewers and editors for their feedback on this manuscript.
Senior Editor
- George H Perry, Pennsylvania State University, United States
Reviewing Editor
- Yonatan Sahle, University of Cape Town, South Africa
Publication history
- Preprint posted: August 19, 2022 (view preprint)
- Received: August 24, 2022
- Accepted: December 16, 2022
- Accepted Manuscript published: December 19, 2022 (version 1)
- Version of Record published: December 30, 2022 (version 2)
Copyright
© 2022, Demarchi et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 545
- Page views
-
- 81
- Downloads
-
- 1
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Allosteric modulation of G protein-coupled receptors (GPCRs) is a major paradigm in drug discovery. Despite decades of research, a molecular-level understanding of the general principles that govern the myriad pharmacological effects exerted by GPCR allosteric modulators remains limited. The M4 muscarinic acetylcholine receptor (M4 mAChR) is a validated and clinically relevant allosteric drug target for several major psychiatric and cognitive disorders. In this study, we rigorously quantified the affinity, efficacy, and magnitude of modulation of two different positive allosteric modulators, LY2033298 (LY298) and VU0467154 (VU154), combined with the endogenous agonist acetylcholine (ACh) or the high-affinity agonist iperoxo (Ipx), at the human M4 mAChR. By determining the cryo-electron microscopy structures of the M4 mAChR, bound to a cognate Gi1 protein and in complex with ACh, Ipx, LY298-Ipx, and VU154-Ipx, and applying molecular dynamics simulations, we determine key molecular mechanisms underlying allosteric pharmacology. In addition to delineating the contribution of spatially distinct binding sites on observed pharmacology, our findings also revealed a vital role for orthosteric and allosteric ligand–receptor–transducer complex stability, mediated by conformational dynamics between these sites, in the ultimate determination of affinity, efficacy, cooperativity, probe dependence, and species variability. There results provide a holistic framework for further GPCR mechanistic studies and can aid in the discovery and design of future allosteric drugs.
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Aurora B, together with IN-box, the C-terminal part of INCENP, forms an enzymatic complex that ensures faithful cell division. The [Aurora B/IN-box] complex is activated by autophosphorylation in the Aurora B activation loop and in IN-box, but it is not clear how these phosphorylations activate the enzyme. We used a combination of experimental and computational studies to investigate the effects of phosphorylation on the molecular dynamics and structure of [Aurora B/IN-box]. In addition, we generated partially phosphorylated intermediates to analyze the contribution of each phosphorylation independently. We found that the dynamics of Aurora and IN-box are interconnected, and IN-box plays both positive and negative regulatory roles depending on the phosphorylation status of the enzyme complex. Phosphorylation in the activation loop of Aurora B occurs intramolecularly and prepares the enzyme complex for activation, but two phosphorylated sites are synergistically responsible for full enzyme activity.