Protein extractionAnimal fossils<\/p>\n
Powder was drilled from the dentin (Supplementary Data\u00a01<\/a>), then mixed with 1.5\u2009ml of 0.6\u2009M HCl for decalcification. The precipitate was washed with ultrapure water at least 3 times until the pH reached approximately 7. Then, 200\u2009\u00b5l of 50\u2009mM NH4HCO3 was added, and the mixture was incubated at 65\u2009\u00b0C for 3\u2009h to extract soluble proteins. The supernatant containing the proteins was then transferred into a new centrifuge tube, and 1\u2009\u00b5g of trypsin (Promega) was added. The mixture was incubated at 37\u2009\u00b0C for 18\u2009h. A 2.5% trifluoroacetic acid (TFA) solution (final concentration 0.1%) was added to stop the reaction. Subsequently, the trypsinized peptides were desalted and purified using C18 ZipTips57<\/a>. For enamel samples, the powder was mixed with 1\u2009ml of 5% HCl37<\/a> for decalcification, and the acid was replaced daily until the reaction ceased. The acid solution containing dissolved enamel peptides was concentrated using a vacuum concentrator, and the peptides were then desalted and purified using C18 ZipTips58<\/a>.<\/p>\n Peptides were eluted into a solution of 80% acetonitrile with 0.1% TFA for MALDI-TOF mass spectrometry analysis. All sample preparations were performed in the dedicated clean room at the Molecular Paleontology Laboratory, IVPP of the Chinese Academy of Sciences in Beijing.<\/p>\n Hominin fossils<\/p>\n An acid etching method was used to extract protein from tooth enamel, modified from the process described in ref. 37<\/a>. Disposable toothbrushes were used to remove surface contaminants from a small area of the enamel for etching. At the same time, the remaining teeth were wrapped with parafilm to prevent contact with any liquids. Before etching, the small enamel area was initially washed with 3% H2O2 for 30\u2009s, followed by a rinse with ultrapure water. Approximately 100\u2009\u00b5l of 5% (v\/v) HCl was placed in the cap of a 1.5-ml microcentrifuge tube. A 2-min etch was performed by immersing the etching region in the HCl solution, and the initial etch solution was discarded. A second etch, lasting 15\u2009min, was carried out in the cap of another separate microcentrifuge tube, and the etch solution was retained. This second etch was repeated, and the etch solutions were combined. After etching, the etched area was treated with 100\u2009\u00b5l of 50\u2009mM ammonium bicarbonate solution for 1\u2009min to neutralize the acid. It was then rinsed with ultrapure water for 30\u2009s and dried. The combined etch solution was then desalted using C18 ZipTips (Thermo Fisher Scientific) and eluted into a solution of 0.1% TFA and 80% acetonitrile (ACN). The peptide mixture (50\u2009\u00b5l) was further divided into three aliquots; one aliquot was composed of 16\u2009\u00b5l, among which 3\u2009\u03bcl was used for the MALDI-TOF mass spectrometry test and 13\u2009\u00b5l was retained as backup in our laboratory; two aliquots (each composed of 17\u2009\u00b5l) were dried for LC\u2013MS\/MS analysis in two independent laboratories. All sample preparation for the experiment was conducted in the dedicated clean room at the Molecular Paleontology Laboratory, IVPP of the Chinese Academy of Sciences in Beijing.<\/p>\n MALDI-TOF mass spectrometry analysis<\/p>\n The peptide mixture was analysed on a Bruker autoflex maX MALDI-TOF mass spectrometer. In detail, 1\u2009\u00b5l of peptide mixture was spotted onto a MTP384 Bruker ground-steel MALDI target plate, and 1\u2009\u00b5l of \u03b1-cyano-4-hydroxycinnamic acid matrix solution (1% in 50% ACN\/0.1% TFA (v\/v\/v)) was added on top. They were mixed, dried, and analysed on the mass spectrometer with a m\/z range of 700\u20133,500. Each sample was analysed in triplicate. The raw data files were processed by mMass (v5.5.0)59<\/a>.<\/p>\n LC\u2013MS\/MS analysis<\/p>\n For each hominin fossil, the eluted peptides from the\u00a0enamel extraction were analysed under DDA mode in triplicate, including two runs on an Orbitrap Fusion Lumos mass spectrometer (Thermo Fisher Scientific) at Capital Medical University, Beijing, and one run on an Orbitrap Exploris 480 mass spectrometer (Thermo Fisher Scientific) at Fudan University. Both devices were coupled to an Easy nLC 1200 HPLC system (Thermo Fisher Scientific).<\/p>\n For the Orbitrap Fusion Lumos at Capital Medical University, the peptides were initially loaded onto a 100\u2009\u03bcm internal diameter\u2009\u00d7\u20092\u2009cm trap column and then separated on a 150\u2009\u03bcm internal diameter\u2009\u00d7\u200915\u2009cm analytical column. Both columns were packed in-house using 3\u2009\u03bcm reversed-phase silica (Reprosil-Pur C18 AQ, Dr. Maisch). The peptides were eluted using a 120\u2009min linear gradient programme (0\u20138\u2009min, 7\u201311% B; 8\u201396\u2009min, 11\u201328% B; 96\u2013108\u2009min, 28\u201340% B; 108\u2013113\u2009min, 40\u201390% B; 113\u2013120\u2009min, 90% B) at a flow rate of 500\u2009nl\u2009min\u22121. Buffer A was 0.1% formic acid in water, and buffer B was 80% acetonitrile and 0.1% formic acid. The MS1 data were acquired across 375\u20131,400\u2009m\/z, with a resolution of 120k at m\/z 200, a 250% AGC target, and a maximum injection time of 50\u2009ms. The MS2 scans were performed with a resolution of 15k at m\/z 200, an AGC target of 100%, a 35% normalized collision energy, and a maximum injection time of 22\u2009ms.<\/p>\n For the Orbitrap Exploris 480 at Fudan University, Shanghai, the peptides were separated on a 75\u2009\u03bcm internal diameter\u2009\u00d7\u200925\u2009cm analytical column, which was packed in-house using reversed-phase silica of 1.9 \u03bcm (Reprosil-Pur C18 AQ, Dr. Maisch). Buffer A was 0.1% formic acid in water, and buffer B was 80% acetonitrile and 0.1% formic acid. An 80\u2009min gradient was used with the following profile: 5\u20138% B, 2\u2009min, at a flow rate of 200\u2009nl\u2009min\u22121; 8\u201344% B, 38\u2009min, 200\u2009nl\u2009min\u22121; 44\u201370% B, 8\u2009min, 200\u2009nl\u2009min\u22121; 70\u2013100% B, 2\u2009min, 200\u2009nl\u2009min\u22121; 100% B, 10\u2009min, 200\u2009nl\u2009min\u22121; 100\u20135% B, 2\u2009min, 200\u2009nl\u2009min\u22121; 5% B, 2\u2009min, 300\u2009nl\u2009min\u22121; 5\u2013100% B, 6\u2009min, 300\u2009nl\u2009min\u22121; 100% B, 10\u2009min, 300\u2009nl\u2009min\u22121. Full mass spectrometry scans were acquired for the first 65\u2009min, after which the column was washed and re-equilibrated for 15\u2009min without data acquisition. The full mass spectrometry data acquisition was conducted across the range of m\/z 350\u20131,600, with a resolution of 60k at m\/z 200. The AGC target was set to \u2018standard\u2019, and the maximum injection time mode was set to \u2018auto\u2019. The MS\/MS spectra were acquired with a resolution of 15k at m\/z 200, a maximum injection time of 30\u2009ms and a normalized collision energy of 30%. The AGC target was also set to standard.<\/p>\n For each animal fossil, the eluted peptides from enamel extractions were analysed for one run under DDA mode on the Orbitrap Fusion Lumos mass spectrometer (Thermo Fisher Scientific), either at Capital Medical University or Fudan University. For the Orbitrap Fusion Lumos at Capital Medical University, the liquid chromatography gradient and mass spectrometry parameters were the same as those of the hominin fossils. The Orbitrap Fusion Lumos at Fudan University was also interfaced with an Easy nLC 1200 HPLC system (Thermo Scientific). The peptides were separated on a 75\u2009\u03bcm internal diameter\u2009\u00d7\u200920\u2009cm analytical column packed with 1.9\u2009\u03bcm reversed-phase silica. Mobile phase A consisted of 0.1% formic acid, and mobile phase B consisted of 80% acetonitrile and 0.1% formic acid. An 80\u2009min gradient was used with the following profile: 2\u20135% B, 3\u2009min, at a flow rate of 200\u2009nl\u2009min\u22121; 5\u201335% B, 40\u2009min, 200\u2009nl\u2009min\u22121; 35\u201344% B, 5\u2009min, 200\u2009nl\u2009min\u22121; 44\u2013100% B, 2\u2009min, 200\u2009nl\u2009min\u22121; 100% B, 10\u2009min, 200\u2009nl\u2009min\u22121; 100\u20135% B, 2\u2009min, 200\u2009nl\u2009min\u22121; 5% B, 2\u2009min, 300\u2009nl\u2009min\u22121; 5\u2013100% B, 6\u2009min, 300\u2009nl\u2009min\u22121; 100% B, 10\u2009min, 300\u2009nl\u2009min\u22121. Full mass spectrometry scans were acquired for the first 65\u2009min, after which the column was washed and re-equilibrated for 15\u2009min without data acquisition. The full mass spectrometry data acquisition was conducted across the m\/z range 350\u20131,600, with a resolution of 60k at m\/z 200, a 100% AGC target, and a maximum injection time of 50\u2009ms. The MS\/MS spectra were acquired with a resolution of 15k at m\/z 200, a 100% AGC target, a 30% normalized collision energy, and a maximum injection time of 30\u2009ms. Blank extractions were processed concurrently to monitor the exogenous contaminants during the procedure.<\/p>\n Data search strategy<\/p>\n MaxQuant (v2.6.0.0)32<\/a>, PEAKS Online (v12)33<\/a> and pFind (v3.2.1)34<\/a> were used to search the raw data. The H. erectus raw files were searched with the corresponding laboratory blanks and modern H. sapiens raw files43<\/a> against the \u2018Hominidae enamel database\u2019, supplemented with the contaminant database. Unspecific digestion was selected in each software. The animal raw files were searched against the \u2018mammal enamel database\u2019.<\/p>\n Database composition<\/p>\n The \u2018Hominidae enamel database\u2019 comprised 13 selected enamel proteins from Hominidae. Besides the commonly used 12 proteins (AHSG, ALB, AMBN, AMELX, AMELY, AMTN, COL17A1, ENAM, KLK4, MMP20, ODAM and TUFT1), SERPINC1 was also added because we identified this protein in Harbin with abundant peptides and elevated deamidation rates (nearly 100%), and this protein was also reported in modern enamel60<\/a>,61<\/a>. The sequences were retrieved from UniProt, downloaded from the \u2018Hominid Palaeoproteomic Reference Dataset\u2019 (https:\/\/zenodo.org\/records\/7728060<\/a>), translated from the genomes of public projects30<\/a>, and specific sequences from published palaeoproteomes11<\/a>. The \u2018mammal enamel database\u2019 was composed of mammal sequences of the same enamel proteins above, retrieved from UniProt using the gene names and \u2018Mammalia (mammals) [40674]\u2019. The contaminant database was composed of a previously published contaminant database62<\/a>, the cRAP database (https:\/\/www.thegpm.org\/crap\/<\/a>), and the contaminant database from MaxQuant (v2.6.0.0). The \u2018Hominidae enamel database\u2019 was accessible through the ProteomeXchange Consortium (Data availability).<\/p>\n PEAKS search<\/p>\n The precursor ion (MS1) mass tolerance was set to 10 ppm, with a fragment ion (MS2) mass tolerance of 0.02\u2009Da for all PEAKS searches, with unspecific digestion. Variable modifications included deamidation (NQ), oxidation (M), hydroxylation (P), phosphorylation (STY), N-terminal pyro-Glu from E, and N-terminal pyro-Glu from Q, with no fixed modifications and up to three modifications allowed per peptide. The peptide length was set to 6\u201345. PSMs were filtered using a false discovery rate (FDR) of 1%, and proteins were filtered with criteria of \u221210logP\u2009\u2265\u200920 and average local confidence (ALC) \u2009\u2265\u200950% (de novo only).<\/p>\n After our initial search with PEAKS, additional variable modifications were included in the second-round search for Hexian H. erectus samples: chlorination and dichlorination of tyrosine residues, dehydration, dioxidation (W), carbonyl E, dioxidation (M), oxidation (HW), ornithine derived from arginine, tryptophan oxidation to kynurenine, tryptophan oxidation to oxolactone, and proline oxidation to pyroglutamic acid. Up to five modifications per peptide were permitted. The peptide length range was set to 6\u201330.<\/p>\n After removing peptides from the contaminant database and those detected in the extraction blank samples or matching multiple genes, the deamidation rates of glutamine (Q) and asparagine (N) were calculated for each sample based on PSM counts. We used auxiliary tools for PSM prediction with PEAKS.<\/p>\n pFind search<\/p>\n We included Deamidation [N], Deamidation [Q], Oxidation [M], Oxidation [P], Oxidation [W], Gln->pyro-Glu[AnyN-termQ], Glu->pyro-Glu[AnyN-termE], Pro->pyro-Glu[P], Phospho[S], Phospho[T], Phospho[Y], Dehydrated[S], Dehydrated[T], Dehydrated[Y], Arg->Orn[R], Dioxidation[M], Dioxidation[W], Thiazolidine[W], Trp->Kynurenin[W], Trp->Oxolactone[W], Ammonia-loss[N], His->Asp[H], Pro->HAVA[P], and Amidated[AnyC-term] as variable modifications, with no fixed modifications included. Spectra FDR was set at 1%, and protein FDR was 10%. The mass range of each peptide was set from 350 to 4,000\u2009Da. Open search was enabled. Within pFind\u2019s modification configuration, the default mass \u2018X\u2019 is preset to that of isoleucine\/leucine. To avoid incorrect identification of the X residue, we set its mass to 6,228.71\u2009Da (Sm41), which substantially exceeds the mass of any natural amino acid. Therefore, any in silico peptide with an X produces abnormal theoretical precursor and fragment ion masses, effectively preventing its matching to experimental spectra during database searches and thus reducing false positives from these ambiguous sequence regions. The other parameters were the same as the settings in PEAKS.<\/p>\n After our first-round search with pFind, some additional variable modifications were selected for inclusion in the second-round search for Hexian H. erectus samples: Chlorination[Y], dichlorination[Y], Carbonyl[E], and Dioxidation[P](Pro->Glu[P]).<\/p>\n MaxQuant search<\/p>\n No fixed modifications were specified. Variable modifications included Deamidation (NQ), Phosphorylation (STY), Gln to N-terminal pyro-Glu, Glu to N-terminal pyro-Glu, Dioxidation (MW), Oxidation (M), Oxidation (P), and Oxidation (W). PSM FDR was set at 1% for all 7 samples. The search also enabled the identification of dependent and secondary peptides. The remaining parameters were set to their default values. After the search, the deamidation rates of N and Q were calculated63<\/a> for each sample.<\/p>\n In addition, we determined the extra variable modification for Hexian H. erectus samples with PEAKS and pFind results. Cl(Y) and diCl (Y) were added for HX-S1 (both post-translational modifications were self-made in Configuration-Modifications), and Cl (Y) and Dehydrated (STY) were included for HX-S2.<\/p>\n Considering the high resolution of the instrumentation used, we also reduced the precursor ion (MS1) mass tolerance to 5 ppm and repeated all analyses in PEAKS and pFind. The results were similar, and the two main SAPs (AMBN 253 and AMBN 273) within the Homo genus were consistently identified (Supplementary Table 8<\/a>). As the final results do not change significantly, we focus mainly on the results from the commonly used 10 ppm search to make our results more comparable to previous studies.<\/p>\n Construction of the consensus protein sequences and phylogenetic analysis<\/p>\n Consensus sequences of endogenous proteins were reconstructed for phylogenetic analysis (Supplementary Data\u00a01<\/a>). Peptides shorter than eight amino acids or with abnormal or artificial post-translational modifications were excluded for the consensus sequence reconstruction of each endogenous protein. Only alleles with a PSM count \u22652, a PSM ratio \u226510%, and an intensity ratio \u226510% were considered reliable and retained. For heterozygous sites, at least two peptides were required per allele. The heterozygous variant site, observed only in the Harbin specimen, was identified with 79 peptides for the V allele and 38 for the M allele at position 273 in AMBN (Table 1<\/a>). We include both alleles of the Harbin consensus in Supplementary Data\u00a02<\/a>. We performed additional AMELY sequence correction as described in Supplementary Note\u00a04<\/a>.<\/p>\n The protein data in the phylogenetic tree include Denisova 3 (ref. 52<\/a>), two Neanderthals (Altai, Vindija33.19)45<\/a>,53<\/a>, two modern humans from the HGDP project (San, HGDP0987; Bougainville, HGDP01027)30<\/a>, a chimpanzee (Pan)29<\/a>, and a Gorilla29<\/a>. We used PartitionFinder2 (v2.1.1)64<\/a> to identify the best partitioning schemes and amino acid substitution models for our dataset. Then, we built a consensus Bayesian phylogenetic tree using the software MrBayes (v3.2.6)65<\/a>, running 8 Markov Chain Monte Carlo (MCMC) chains for 1 million iterations in 2 independent runs. Sampling was done every 500 generations, and the first 200,000 iterations were discarded as burn-in. The tree was plotted with FigTree v1.4.4 (http:\/\/tree.bio.ed.ac.uk\/software\/figtree<\/a>).<\/p>\n DNA analysisSliding window analysis<\/p>\n Variant sites were identified where either Denisova 3 or two Neanderthals (Vindija33.19 and Altai) differed from two modern African diploid genomes (S_Khomani_San-1.DG, S_Mandenka-2.DG) from the phased Simons Genome Diversity Project panel (SGDP) (https:\/\/sharehost.hms.harvard.edu\/genetics\/reich_lab\/sgdp\/phased_data2021<\/a>)66<\/a>,67<\/a>. For each of 129 overlapping 20\u2009kb windows (with 2\u2009kb steps), pairwise matching rates were calculated for variants between each African haploid and each archaic artificial haploid of either Denisova 3 or the Altai Neanderthal (where each variant is randomly assigned to a haploid to preserve all variants), covering in total 276\u2009kb of DNA sequence surrounding the rs564905233 SNP. For each window, statistical significance of the differences between the African\u2013Denisovan and African\u2013Neanderthal matching rates were assessed using a Wilcoxon rank-sum test, and a permutation test (n\u2009=\u20091,000 permutations). The Wilcoxon rank-sum test was used to compare distributions of pairwise matching rates, with W representing the sum of ranks in the African\u2013Denisovan group. The permutation test was used to evaluate differences in mean matching rates. No adjustments were made for multiple comparisons, since each window was analysed independently to identify localized regions of divergence. All statistical tests were performed as two-sided tests. Full results, including W, exact P values, and the number of comparisons per window, are provided in Supplementary Data\u00a03<\/a>.<\/p>\n Ethics statement<\/p>\n Permission to test for ancient proteins in the human specimens from this study was granted by the collection room of the IVPP, the Hexian Culture, Tourism, and Sports Bureau, and the Luanchuan County Culture, Radio, Television, and Tourism Bureau. The work was conducted in collaboration with local researchers, who are co-authors because of their contributions to assembling archaeological materials and\/or discussions that informed the study.<\/p>\n Reporting summary<\/p>\n