Bone samples
Herbivorous Edmontosaurus sp. (Hadrosauridae) sacrum bone fossils were excavated from the Upper Cretaceous zone of the Hell Creek Formation in Harding County, South Dakota, USA. A 22 kg sample from this duck-billed dinosaur fossil, together with samples of the sediment in which the fossil was found, was donated to and accessioned at the repository of the Victoria Gallery & Museum of the University of Liverpool under UOL GEO.1.
The modern bone (control) was from a common turkey (Meleagris gallopavo). Small bone segments were dried in an oven at 60°C for several hours in preparation for crushing (powderisation), and analysis following the same protocols as used for the Edmontosaurus samples.
Motion Photogrammetry
Motion photogrammetry procedures were used to capture a digital 3D model of the Edmontosaurus sp. bone fossils prior to sampling.
The model was reconstructed using 1061 images in Agisoft Metashape software version 1.7.2. Sony alpha ILCE-6400 cameras were used with Sony 50 mm Macro F2.8 FE Lens. Even diffuse lighting was used, and for the smaller bones light tents were used. Images were saved in both RAW and fine JPG format, and JPGs were used for the reconstruction of the models. Due to the level of granular detail in the objects, and the controlled conditions, high resolution (100% of image resolution) was used for the camera alignment process. The surfaces of the specimens were highly matte. This permitted the generation of high-resolution depth maps without concern for the potentially erroneous influence of specular highlights. Textures were also generated at high resolution (8k JPGs) using the software mosaic algorithm. As the specimens were not overly complex (i.e., there were no appendages and few occluding surfaces), only a single JPG texture file was generated for each skeletal fragment assigned to this specimen. The final models varied between 2–5 million faces and have been retained for future detailed measurements. The models were down-sampled to 64 thousand faces for online loading speeds while retaining adequate accuracy for general inspection (see Supplementary Table S1).
The samples were ground bone shards (cross-sections of 1–3 mm width, Supplementary Figure S1) prepared using a mortar and pestle. These were cleaned with powdered bicarbonate and hot-tap water (~ 50°C) before final rinsing with de-ionised water. A 50-micron stackable zooplankton sieve was used to filter the particles onto a freshly cut piece of aluminium foil for transfer into new vials, ready for use with FTIR or further sample preparation for LC-MS.
FTIR
FTIR was performed using an Attenuated Total Reflectance accessory (ATR) with a germanium window on a Bruker Vertex 70© equipped with Deuterated Lanthanum α Alanine doped TriGlycine Sulphate (DLaTGS) detector. Each spectrum combined an average of 32 scans, with a resolution of between 2–4 cm⁻¹ in the range 4,500 to 650 cm⁻¹. Spectra were collected and analysed with OPUS software and compared with authentic Ca₃PO₄ from the library (©Nicodom, 2014). Peak intensities correspond to the moiety abundance in the sample absorbing the energy at a certain frequency.
For this measurement samples were ground to powder with particle sizes of no more than 50 µm [44]. This provides a sufficient contact area between the sample and the crystal.
XPol
Thin sections of UOL GEO.1 were prepared and imaged using a Motic Polarizing Microscope BA310 with a Sony ILCE-7RM4 detector. Images from several focal planes were collected then stacked using Photoshop 24.5.
LC-MS bottom-up proteomics
Initial analyses were completed in the Materials Innovation Factory (MIF) laboratories, University of Liverpool, using high resolution LC-MS. The preliminary results led to further analyses at the Centre for Proteome Research, also at the University of Liverpool, which are reported herein.
Bone samples were ground by hand in a mortar & pestle. Twenty milligrams each of Edmontosaurus bone, turkey bone, and pure collagen (Bovine tendon collagen, lyophilized fibrous powder, Sigma-Aldrich product #5162) was dispensed into separate polypropylene microcentrifuge tubes. Each sample was treated with aqueous ammonium bicarbonate (AmBic, 80 µl, 25 mM) and RapiGest SF Surfactant solution (1% RapiGest solution in AmBic, 5 µl, Waters) with continuous gentle shaking (450 rpm, 80˚C, 10 min.). Cysteine reduction was then performed by the addition of dithiothreitol (DTT, 11.1 mg/ml in 25 mM AmBic, 5 µl). After mixing and incubation (60˚C, 10 min.) alkylation of free thiols was performed using iodoacetamide (46.6 mg/ml in 25 mM Ambic, 5 µl, 30 min. in the dark). Excess iodoacetamide was quenched with DTT (4.7 µl as above), and samples were acidified (neat trifluoroacetic acid, 2 µl) to a pH of 2 or less (checked with pH indicator paper). Digestion was carried out with trypsin (Promega sequencing grade, 0.2 µg/µl in 50 mM aqueous acetic acid) with incubation (37˚C, 16 hr). Following centrifugation (13000g, 15 min., 4˚C) the supernatants were transferred to clean microcentrifuge tubes and stored frozen until analysis.
Samples were analysed using nanobore reversed-phase chromatography (Ultimate™ 3000 RSLC, Thermo Scientific™, Hemel Hempstead) coupled to hybrid linear quadrupole/orbitrap mass spectrometry (Q Exactive™ HF Quadrupole-Orbitrap™, Thermo Scientific™) equipped with a nanospray ionization source. Samples (2 µl) were loaded onto the trapping column (Thermo Scientific, PepMap100, C18, 300 µm X 5 mm) equilibrated in aqueous formic acid (0.1%, v/v) using partial loop injection over seven minutes at a flow rate of 12 µl/min. After the direction of eluent flow was reversed components were transferred to and resolved on an analytical column (Easy-Spray C18 75 µm x 500 mm, 2 µm particle size) equilibrated in 96.2% eluent A (water/formic acid, 100/0.1, v/v) and 3.8% eluent B (acetonitrile/water/formic acid, 79.95/19.95/0.1, v/v/v) and eluted (0.3 µl/min.) with a linear increasing concentration of eluant B (min./% B; 0/3.8, 30/50). The mass spectrometer was operated in a data-dependant positive ion mode (FWHM 60,000 orbitrap full-scan, automatic gain control (AGC) set to 3e⁶ ions, maximum fill time (MFT) of 100 ms). The seven most abundant peaks per full scan were selected for high energy collisional dissociation (HCD, 30,000 FWHM resolution, AGC 1e⁵, MFT 300 ms) with an ion selection window of 2 m/z and a normalised collision energy of 30%. Ion selection excluded singularly charged ions and ions with ≥ + 6 charge state. A 60 s dynamic exclusion window was used to avoid repeated selection of the same ion for fragmentation.
Survey analyses of each sample were first used to determine the amount of each sample, calculated by extrapolation, needed to give a full orbitrap scan base peak intensity (BPI) of 1–2e⁹. These analyses were performed with a compacted 15 min. gradient (Supplementary Table S2). Based on these BPI results, 2 µl of neat Edmontosaurus and Hadrosaurus samples were used for the full analysis. The modern turkey sample was diluted 1:100 and Bovine collagen sample was diluted 1:1000 in water/acetonitrile/trifluoroacetic acid (97/3/0.1, v/v/v).
Particular attention was given to the order of sample injections to avoid carry-over from analysis to analysis. Four 30 min. blank analyses (injection solvent only) were performed after the turkey and bovine samples to minimize carry over, then the three fossilized samples were analysed on the 1 hr program.
Database searches
Analysis of the LC-MS/MS data was performed using MASCOT PD using the Swissprot, UniCow, UniTurkey and UniChick databases. The search parameters included cysteine carbamidomethylation, variable methionine oxidation, a precursor mass tolerance of 10 ppm, a product mass tolerance of 0.01 Da, and a maximum of one missed cleavage.
The data files were also imported onto Peaks (Peaks 11 Software) for further database searching. In this case the search parameters included cysteine carbamidomethylation, methionine oxidation, variable lysine and proline oxidation, a precursor mass tolerance of 10 ppm, a product mass tolerance of 0.01 Da, and a maximum of one missed cleavage. This software permits database searching for multiple post translation modifications (PTMs) as opposed to MASCOT PD that cannot be searched with more than two PTMs. The ‘SPIDER’ feature uses an algorithm that is specially designed to detect peptide mutations and performs cross-species homology search. The Edmontosaurus sample was searched against the Swissprot database, bovine collagen (96%) was searched against UniCow, and the modern turkey sample was searched against both UniChick and UniTurkey databases. The contaminants (cRAP) database was also included in each search.
LC-MS/MS of hydroxyproline
After being frozen with liquid nitrogen, both fossilized and turkey bone samples were manually crushed to a fine powder with a mortar and pestle. One-gram portions of the powdered bone or sediment were dispensed into polypropylene microcentrifuge tubes, suspended in water (1 ml), mixed vigorously, sonicated in a bath sonicator (30 min.), centrifuged (2000g, 15 min.), and the supernatants transferred to new tubes. The extraction procedure was repeated on the pellet by adding methanol (1 ml) and the samples were mixed, sonicated, and centrifuged as above. The supernatants were pooled and reserved for future bottom-up proteomics. The pellets were then treated with HCl (2 ml, 6 N) and incubated (2 hrs, 60°C) before the samples were dried in a vacuum centrifuge. The HCl treatment was repeated until the samples ceased effervescing after HCl addition, each time with drying in a vacuum centrifuge between acid treatments. Generally, it took two or three such treatments before effervescing ceased. The final dried samples were treated again with HCl (500 µl, 6 N) and incubated (12 hrs, 120°C) to effect protein hydrolysis. The samples were dried overnight in a vacuum centrifuge and then treated with n-butanolic HCl (300 µl, 3 N), incubated (2 hrs, 60°C) to make the butyl esters, and dried again in a vacuum centrifuge. Lastly, the samples were reconstituted in water (200 µl), mixed vigorously and centrifuged (5 min., 16,000g, room temperature). The supernatants were transferred to HPLC vials and aliquots (typically 10 µl), injected onto a reversed-phase HPLC column (Phenomenex Kinetex, 2.6 µm Polar C18, 100 Å, 100 x 2.1 mm), equilibrated in eluant A, and eluted (100 µl/min) with a step-wise linearly increasing concentration of eluant C (acetonitrile/formic acid, 100/0.1, v/v; min/%C, 0/1, 5/1, 20/25, 22/1, 60/1). The effluent from the column was passed through an electrospray ionization source connected to a hybrid linear ion trap/orbitrap mass spectrometer (Thermo Scientific Orbitrap LTQ XL) scanning either in the positive linear ion trap (for quantitative measurements) or the positive orbitrap (for accurate m/z measurements) modes.
For quantitative measurements, the mass spectrometer was set to fragment preselected parent ions under standard MS/MS fragmentation conditions for the butyl esters of hydroxyproline (Hypbe, MH⁺ at m/z 188). For the accurate m/z measurements the orbitrap was scanned (FWHM of 100,000 at m/z 400) immediately after calibration with LTQ ESI Positive Ion calibration solution mixture. Data were collected and interrogated with instrument manufacturer-supplied software (Xcalibur 2.05).
Several control samples were included with each batch of bone samples. These were negative control samples devoid of added bone extracts (in triplicate), sediment in which the fossil was sitting at the excavation site (20 mg/sample, in triplicate), pure collagen (Sigma, 20 mg/sample, in triplicate), and authentic Hyp standard in a range of amounts (typically 0, 2, 10, 20 and 50 nmoles/sample, in duplicate). These samples were prepared and processed with each batch of bone samples.
The order in which samples were analysed was carefully arranged. Injections of water (solvent blanks) were used at the start of the analysis of each batch of samples to check that there were no peaks for Hypbe resulting from carry-over from the analysis of previous sample batches. After verification that the LC/MS system was clean a typical order of sample analysis was: negative control samples 1–3; water blank #1; sediment samples 1–3, water blank #2, Edmontosaurus fossilized bone samples 1–3, water blanks #3 and #4, turkey bone samples 1–3, water blanks #5 and #6, collagen samples 1–3, water blanks #7 and #8, Hypbe standards 1–10, water blanks #9–#11.
The data from the standard Hypbe samples were processed by plotting the known amount of Hyp per sample against the measured chromatographic peak areas corresponding to the Hypbe peak. The trendline equation was then used to interpolate or extrapolate the amount of Hyp in each sample.