Next Article in Journal
Role of Gut Microbiota through Gut–Brain Axis in Epileptogenesis: A Systematic Review of Human and Veterinary Medicine
Next Article in Special Issue
Crystallographic Texture of the Mineral Matter in the Bivalve Shells of Gryphaea dilatata Sowerby, 1816
Previous Article in Journal
Comparison of Different Dietary Fatty Acids Supplement on the Immune Response of Hybrid Grouper (Epinephelus fuscoguttatus × Epinephelus lanceolatus) Challenged with Vibrio vulnificus
Previous Article in Special Issue
Evolution of the Family Equidae, Subfamily Equinae, in North, Central and South America, Eurasia and Africa during the Plio-Pleistocene
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fossil Biomarkers and Biosignatures Preserved in Coprolites Reveal Carnivorous Diets in the Carboniferous Mazon Creek Ecosystem

1
Western Australian Organic and Isotope Geochemistry Centre, The Institute for Geoscience Research, School of Earth and Planetary Sciences, Curtin University, Kent Street, Bentley, WA 6102, Australia
2
Department of Earth & Planetary Sciences, Yale University, 210 Whitney Avenue, New Haven, CT 06511, USA
3
Division of Geological and Planetary Sciences, California Institute of Technology, 1200 E. California Blvd., Pasadena, CA 91125, USA
4
Research School of Earth Sciences, The Australian National University, Canberra, ACT 2601, Australia
5
The Field Museum, 1400 S Lake Shore Dr., Chicago, IL 60605, USA
6
Organic Geochemistry Unit, Institute of Geoscience, Christian-Albrechts-University, 24118 Kiel, Germany
*
Authors to whom correspondence should be addressed.
Biology 2022, 11(9), 1289; https://doi.org/10.3390/biology11091289
Submission received: 2 July 2022 / Revised: 25 August 2022 / Accepted: 26 August 2022 / Published: 30 August 2022
(This article belongs to the Special Issue Paleontology in the 21st Century)

Abstract

:

Simple Summary

Coprolites (fossilised faeces) can preserve important dietary information through geological time, offering insights into extinct animal diets. When digestion of dietary items leaves no unambiguous morphology to reconstruct the food spectrum of a coprolite producer, preserved biomolecular information can offer unique perspectives into the individual dietary composition and trophic relationships in ancient ecosystems. In this study we combine a uniquely diverse array of chemical techniques to demonstrate that biomarkers and macromolecular biosignatures from Carboniferous coprolites can reveal the dietary spectrum and trophic position of their extinct producers: an overwhelming abundance of cholesteroids, biomarkers of animal cholesterol, and an animal-affinity of the preserved macromolecular phase revealed by the statistical analysis of in situ Raman spectra, indicate a likely carnivorous diet for the coprolite producer. The presence of intact primary metabolites, such as sterols, and informative fossilization products of biopolymers, demonstrates the significance of siderite (iron carbonate) concretions in the exceptional preservation of biomolecular information in deep time, facilitated by the rapid encapsulation and remineralisation of organic matter within days to months.

Abstract

The reconstruction of ancient trophic networks is pivotal to our understanding of ecosystem function and change through time. However, inferring dietary relationships in enigmatic ecosystems dominated by organisms without modern analogues, such as the Carboniferous Mazon Creek fauna, has previously been considered challenging: preserved coprolites often do not retain sufficient morphology to identify the dietary composition. Here, we analysed n = 3 Mazon Creek coprolites in concretions for dietary signals in preserved biomarkers, stable carbon isotope data, and macromolecular composition. Cholesteroids, metazoan markers of cholesterol, show an increased abundance in the sampled coprolites (86 to 99% of the total steranes) compared to the surrounding sediment, indicating an endogenous nature of preserved organics. Presence of unaltered 5α-cholestan-3β-ol and coprostanol underline the exceptional molecular preservation of the coprolites, and reveal a carnivorous diet for the coprolite producer. Statistical analyses of in situ Raman spectra targeting coprolite carbonaceous remains support a metazoan affinity of the digested fossil remains, and suggest a high trophic level for the coprolite producer. These currently oldest, intact dietary stanols, combined with exquisitely preserved macromolecular biosignatures in Carboniferous fossils offer a novel source of trophic information. Molecular and biosignature preservation is facilitated by rapid sedimentary encapsulation of the coprolites within days to months after egestion.

1. Introduction

The reconstruction of ancient trophic networks is pivotal to our understanding of ecosystem function and change through time. However, inferring dietary relationships in enigmatic ecosystems dominated by organisms without modern analogues, such as the Carboniferous Mazon Creek fauna, has previously been considered challenging. Coprolites, fossil faecal materials, can offer unique insights into the diets and trophic relationships of extinct life forms in deep time. However, fossil faecal matter can be difficult to interpret, due to the digestion and thus substantial degradation of organismal morphologies. Coprolites can be linked to a producer based on their shape, mineralogy, and geological context of the specimen as well as the presence of any identifiable remains, e.g., [1,2,3,4,5,6,7]. Recently, methods of identification have been expanded to include 13C/12C, 15N/14N and DNA analysis [8,9]. However, many coprolite specimens still remain of ambiguous origin and composition and may not contain DNA remains, especially in samples from deep time.
The preservation of biomolecules in deep time is primarily dictated by the early diagenetic chemo-environment, e.g., [10,11,12]. Delicate details in soft tissues are only known from Konservat Lagerstätten [13,14,15,16]; Lagerstätten that preserve soft tissues offer generally more representative insights into extinct biodiversity than fossil sites biased towards the preservation of only hard tissues. The conditions, which result in soft-tissue preservation, are often conducive to the preservation of detailed molecular information, which is otherwise lost during heterotrophic reworking prior to preservation. Carbonate concretions are known to frequently contain carbonaceous fossil soft tissues, e.g., [17], and are thus inferred to form through rapid authigenic mineralisation, which halts significant destruction of the specimen [14,18,19]. During fossilization, lipids are transformed into stable derivatives preserving their original hydrocarbon skeleton. The resulting biomarkers are commonly used to identify original sources of organic matter input into sediments [20].
Biomolecules such as sterols are present in most eukaryotes and contain structural features specific to their function within a group of organisms [21]. Lipids such as sterols tend to be better preserved in sediments compared to other classes of biomolecules (e.g., carbohydrates, proteins) which are often completely degraded under anoxic conditions [22]. Through early diagenesis, biological steroids are defunctionalised as a result of the activity of microbes and clay catalysts, while typically retaining most of their isomeric characteristics [23]. These compounds can include stanols, sterenes, diasterenes and A/B/C-ring monoaromatic steroids. Further diagenetic and catagenetic reactions will result in isomerisation and aromatisation to form more thermodynamically stable configurations [23,24,25], primarily steranes and triaromatic steroids, which can reside in the sediments for up to hundreds of millions of years. A diagenetic continuum of such steroid hydrocarbons was identified in a study of a Devonian calcite concretion from the Gogo Formation, an exceptional Konservat Lagerstätte in the Canning Basin of Western Australia, by Melendez et al. [26]. Here, microbially mediated eogenetic processes were determined to have resulted in the parallel preservation of stenols and their diagenetic products encompassing steranes and triaromatic steroids. Such instances of preservation offer the opportunity to study not only the characteristic biomolecules from the fossilised specimen, but also the post-depositional processes, which transform these biomolecules and the broader taphonomic history of coprolites.
In contrast, structural biomolecules shared among all organisms, such as proteins, tend to crosslink with oxidation products of lipids and sugars via advanced glycoxidation and lipoxidation reaction schemes to form insoluble complex organic matter composed of heteroatom-rich polymers. The resulting insoluble organic matter retains chemically altered, but not unrecognisable evidence of original biosignatures.
This compositional complexity and potential for chemical transformation through digestive and diagenetic processes requires the combination of complementary chemical analyses. Coprolites have been known to preserve lipid biomarkers, e.g., [27,28,29,30,31], which can be used to reconstruct various molecular inputs including direct dietary information and the processes, which alter these molecules through the digestive tract of the producer [28,32]. In studies of faecal samples, 5β-stanols are common reduced products of dietary sterols including cholesterol, campesterol, sitosterol and stigmasterol [27]. There is a predominance of 5β-cholestan-3β-ol (coprostanol) in human faeces [33], while herbivorous mammal faeces contain mostly 5β-campestanol and 5β-stigmastanol [27,33].
In addition, insoluble organic matter resulting from the diagenetic crosslinking of structural biomolecules has been shown to preserve biologically informative heterogeneities for fossils from the Mazon Creek locality [34]. Extracting dietary information from molecular biomarkers and biosignatures preserved in coprolites opens up new opportunities for tracing trophic networks through time. For example, the presence of high amounts of cholestane in coprolite samples has been attributed to a predominantly carnivorous or omnivorous diet, e.g., [28,29], due to ubiquity of cholesterol in animal tissue. In contrast, the presence of an array of phytosterols and long-chain n-alkanes derived from leaf waxes can indicate a herbivorous diet, e.g., [29,33,35]. A fraction of the organic material is chemically transformed into insoluble complex organic matter following ingestion, digestion and egestion, potentially preserving dietary tissue information in form of informative heterogeneities. On the other hand, primary dietary information can also be obtained from indigestible components, which pass through the gastrointestinal system without chemical modification.
The accuracy and detail of the molecular dietary reconstruction is reliant on the degree of biomolecule and biomarker preservation, and diagenetic transformation. The processes, which are responsible for preservation of soft tissue (e.g., rapid burial, rapid mineral growth [14]) are also those which have been demonstrated to preserve biomarkers of low maturity and intact biomolecules, e.g., [11,12,26]. Therefore, exceptionally preserved coprolites that were sealed from environmental influence during earliest diagenesis, such as those from the Mazon Creek Lagerstätte, are likely to preserve biomolecular information alongside soft tissues and thus are ideal candidates for ancient dietary reconstructions.

Palaeoenvironmental Setting

The Mazon Creek Konservat Lagerstätte is one of the most productive fossil Lagerstätten worldwide, with over 350 species of plants and 465 of animals identified [36,37]. The site is renowned for its preservation of delicate soft tissue fossils within iron carbonate (siderite) concretions, of which there are many hundreds of thousands collected. Palynological and palaeobotanical studies of the Lagerstätte, e.g., [36,38,39] have determined its age to be Middle Pennsylvanian (Westphalian D) (306–311 Ma). The siderite concretions are found in the lower 3–8 m of the Francis Creek Shale Member of the Lower Carbondale Formation, overlying the Colchester Coal (No. 2) Member [40,41,42,43,44]. The Mazon Creek site was interpreted as representing a river delta system, e.g., [41,42,43,44,45], wherein some large-scale events such as storms or flash flooding caused the rapid burial of massive amounts of organisms, e.g., [42]. More recently the palaeoenvironment has been re-evaluated to represent a low energy, brackish marine environment where under anaerobic conditions input from peaty forests provided a means of rapid burial of organisms [43]. The abundance of siderite concretions and its occasional co-occurrence with pyrite indicates minimal marine input and low sulfate concentrations, where iron initially reacted with H2S (where available) to form pyrite. Under anoxic conditions the presence and activity of methanogenic archaea and methanotrophic bacteria initiated the precipitation of iron carbonate, e.g., [42,46,47,48]. This has also been supported in geochemical studies of 34S/32S, 13C/12C and 18O/16O isotopic compositions of the Mazon Creek concretions [43]. Preservation occurs during early microbial decay of deposited organisms, wherein a ‘proto-concretion’ is formed by the breakdown of degrading organic matter into fatty acids, resulting in mineral precipitation of primarily siderite in an iron-rich environment with limited sulfate [43]. Growth and decay experiments, e.g., [49,50,51,52] have demonstrated that this process is initiated within the weeks after deposition. Studies [53,54] have identified that these processes occurred rapidly (within weeks to months) in carbonate concretions formed around decaying soft tissues of tusk-shells.
Here, we analyze n = 3 Mazon Creek coprolites preserved in carbonate concretions for dietary signals in preserved biomarkers, stable carbon isotope data, insoluble fossil organic matter composition, and preserved morphology, to provide comprehensive and complementary insights into the trophic structure of the enigmatic Mazon Creek ecosystem and improve our understanding of the role of concretion formation in the preservation of biomolecules in deep time.

2. Materials and Methods

2.1. Sample Preparation and Extraction

Three coprolite fossils in siderite concretions were prepared for analysis (Figure 1). The samples used in this study were obtained from the Chicago Field Museum. Samples had previously been collected from the Pit 11 coal strip mine in Kankakee County, Illinois, from the Francis Creek Shale Member of the Carbondale Formation.
One half of each concretion was cut using a handheld Dremel rotary tool with a diamond blade (previously cleaned by sonicating in a mixture of dichloromethane (DCM) and methanol (MeOH) (9:1 v/v) for 15 min intervals, separating the central coprolite fossil (referred to herein as the ‘fossil’) from the surrounding rock (‘matrix’) for each sample. The different sections of each concretion were washed by repeated sonication in 15 min intervals to trace removal of external contamination in a mixture of DCM:MeOH (9:1 v/v), before being ground using pre-annealed (450 °C for 3 h) ceramic mortars and pestles. In between sample treatments, pre-annealed sand was ground to clean mortars.
The ground sample material was Soxhlet extracted in individual pre-extracted cellulose thimbles (Soxhlet extracted three times for 24 h using a mixture of DCM:MeOH (9:1 v/v)) for 72 h. A procedural blank of a pre-extracted thimble was run alongside each extraction. The samples were filtered through activated copper powder to remove elemental sulfur. Small scale column chromatography (5 cm × 0.5cm i.d.) using silica gel activated at 160 °C for 24 h was used to separate the total lipid extracts into aliphatic (4 mL n-hexane), aromatic (4 m n-hexane:DCM (9:1 v/v)), porphyrin (4 mL n-hexane:DCM (7:3 v/v)) and polar (4 mL DCM:MeOH (1:1)) fractions for analysis.

2.2. Gas Chromatography-Mass Spectrometry (GC-MS)

Full scan gas chromatography-mass spectrometry analysis (GC-MS) was performed on the aliphatic fractions using an Agilent 7890B GC with a DB-1MS UI capillary column (J and W Scientific, 60 m, 0.25 mm i.d., 0.25 μm film thickness) coupled to an Agilent 5977B MSD. Aromatic fractions were analysed on an Agilent 6890N GC with a DB-5MS UI capillary column (J and W Scientific, 60 m, 0.25 mm i.d., 0.25 μm film thickness) coupled to an Agilent 5975B MSD. The GC oven was ramped from 40 °C to 325 °C at a rate of 3 °C/min with initial and final hold times of 1 min and 30 min, respectively.
Saturated and aromatic steroids and hopanoids were quantified using GC-MS analyses on an Agilent 6890 GC coupled to a Micromass Autospec Premier double sector MS (Waters Corporation, Milford, MA, USA). The GC was equipped with a 60 m DB-5 capillary column (0.25 mm i.d., 0.25 μm film thickness; Agilent J&W Scientific, Agilent Technologies, Santa Clara, CA, USA), and helium was used as the carrier gas at a constant flow of 1 mL/min. Samples were injected in splitless mode into a Gerstel PTV injector at 60 °C (held for 0.1 min) and heated at 260 °C min−1 to 300 °C. The MS source was operated at 260 °C in EI mode at 70 eV ionization energy and 8000 V acceleration voltage. All samples were injected in n-hexane to avoid deterioration of chromatographic signals by FeCl2 build-up in the MS ion source through use of halogenated solvents [55]. The GC oven was programmed from 60 °C (held for 4 min) to 315 °C at 4 °C min−1, with a total run time of 100 min. Saturated steranes and hopanes were quantified using metastable reaction monitoring (MRM) in M+ → 217 and M+ → 191 transitions, respectively. Mono- and triaromatic steroids were detected using selected ion recording (SIR) under magnet control of base ions m/z 253 and 231, respectively. All ratios and abundance proportions are reported uncorrected for differences in MS-response. Saturated and aromatic steroid hydrocarbons as measured using MRM are shown in Figure 2 and Figure 3.

2.3. Gas Chromatography-Isotope Ratio-Mass Spectrometry (GC-irMS)

Stable isotope ratio mass spectrometric analyses of individual compounds were performed on the aliphatic fractions of each sample to determine δ13C, using a Thermo Scientific Trace GC Ultra coupled to a Thermo Scientific Delta V Advantage mass spectrometer via a GC isolink and a Conflo IV. The reactors consisted of a combustion interface containing a ceramic tube lined with NiO and filled with NiO and CuO, held at 1000 °C. The programs used for the GC column, carrier gas, injector conditions and oven temperatures were identical to those used for GC-MS analysis as described above. The gas chromatography-isotope ratio-mass spectrometer (GC-irMS) measures the δ13C values by monitoring the CO2 produced by the sample and measuring the response of ions of m/z 44, m/z 45 and m/z 46, relative to the reference gas of known δ13C content.

2.4. Bulk Stable Carbon Isotopes

Ground residue from lipid biomarker extractions were treated with hydrochloric acid (4M) to remove carbonate mineral via repeated addition of fresh acid solution, stirring and heating at 50 °C until gas production ceased. Samples were subsequently washed with Milli-Q water and freeze-dried to remove excess water. δ13C analyses were performed using a Thermo Flash 2000 HT elemental analyser (EA) connected to a Delta V Advantage isotope ratio monitoring mass spectrometer (irMS) via a Conflo IV. Samples were weighed (approximately 6 mg) in triplicate into tin cups (SerCon) and combusted to CO2 in the nitrogen-carbon reactor (1020 °C). CO2 passed through the Conflo IV interface into the irMS, which measured m/z 44, 45 and 46. δ13C values were calculated by Thermo Isodat software and normalised to the international VPDB scale by multi-point normalisation using the standard reference materials NBS 19 (+1.95‰) and L-SVEC (−46.60‰) [56]. The standard reference material IAEA-600 was measured to evaluate the accuracy of the normalization. The normalized δ13C values of IAEA-600 from these measurements were within ±0.1‰ of the reported value of −27.77‰ [56].

2.5. Polar Compound Analysis

Polar fractions were analysed at Leeder Analytical (Victoria, Australia). Fractions were dried and internal standard (13C-cholesterol) added, and were then combined with N,O-bis(trimethylsilyl)fluoroacetamide and trimethylchlorosilane (99:1) and heated (60 °C) for 20 min. Samples were dissolved in toluene (500 μL) before analysis. Gas chromatography-tandem mass spectrometry (GC-MS/MS) was performed on an Agilent 7890B Gas Chromatograph with a DB-5MS UI column (30 m × 0.25 mm 0.25 um film) coupled to an Agilent 7000D Triple Quadruple Mass Spectrometer. Results were quantified against sterol trimethylsilyl-derivative standards.

2.6. In Situ Raman Microspectroscopy and ChemoSpace Analysis of Spectral Data

The set of analysed coprolites was microscopically screened for evidence of preserved carbonaceous matter characterized by a dark brown-to-black colouration (Figure 1), and two coprolites (FMNH PE 52316, FMNH PE 52336) with suitable preservation were identified. A total of n = 35 carbonaceous vertebrate, annelid, non-annelid invertebrate, and plant fossils from the Mazon Creek locality (Supplementary Table S1: specimen details) and the two carbonaceous coprolites (FMNH PE 52316, FMNH PE 52336) were surface-cleaned with 70% EtOH, and subjected to high-resolution in situ Raman microspectroscopy in the Department of Earth and Planetary Sciences at Yale University. Raman microspectroscopy was performed using a Horiba LabRam HR800 with 532 nm excitation (holographic notch filter; 20 mW at the sample surface). The spectra were obtained in LabSpec 5 software, and the instant processing included only a standard spike removal. Raman scattering was detected by an electron multiplying charge-coupled device (EM-CCD) following dispersal with an 1800 grooves/mm grating and passing through a 200 μm slit (hole size 300 μm). The spectrometer was calibrated using the first order Si band at 520.7 cm−1. Ten spectra were accumulated in the 500–1800 cm−1 region, also known as ‘organic fingerprint region’, for 10 s exposure time each, at 32× magnification. All spectra were analyzed in an identical fashion in SpectraGryph 1.24 spectroscopic software: A conservative adaptive baseline (30%) was fitted, no baseline offset was imposed, and all spectra were normalized to the common highest peak. Relative intensities (n = 53; arbitrary units) at pre-selected informative band positions ([34,57]; listed in the Supplementary Information) were exported using the ‘Multicursor’ function in SpectraGryph 1.2. The resulting variance-covariance matrix was exported into PAST 3 (file available as Source Data), and two separate sample identification strategies were applied: One data matrix contains the extracted spectroscopic signals and binary characters identifying carbonaceous plant, vertebrate, annelid, non-annelid invertebrate remains and coprolites as separate tissue types. A second matrix uses a more agnostic approach and contains, in addition to the spectroscopic signals, only binary characters identifying samples as plants, annelids, non-annelid invertebrates, and vertebrates; in this data matrix coprolite tissue affinity was coded as ‘unknown’. Both data matrices are available as Source Data. Endogeneity of organic matter in the carbonaceous films associated with fossil morphology has previously been demonstrated [57], and is here separately assessed using lipid biomarkers. A Canonical Correspondence Analysis of the first data matrix revealed the diagnostic molecular features distinguishing the fossil coprolites from other fossil soft tissues from the Mazon Creek locality (Figure 4A). A second Canonical Correspondence Analysis allowed the coprolite samples to locate in the compositional space (ChemoSpace) based on the affinity of the digested, fossil tissues (Figure 4B). Canonical Correspondence Analysis is a discriminant multivariate analysis that distinguishes previously identified groups of samples in a ChemoSpace, and, if samples of unknown affinity are included, reveals their affinity (which is here translated into the dietary spectrum of a coprolite producer). Due to the discriminant, comparative nature of the analysis and spectroscopic data, axis loadings do not have a dedicated unit. The impact of all n = 53 extracted relative intensities is represented by colour-coded (see caption) ChemoSpace vector arrows in Figure 5 (Figure 5B corresponding to the ChemoSpace shown in Figure 4A, and Figure 5C corresponding to the ChemoSpace shown in Figure 4B). Functional groups were identified using Lambert et al. [58], and have been previously published [34,57]; select functional groups that are characteristically enriched in coprolites and reveal the affinity of their digested tissues are labelled in Figure 5.

2.7. X-ray Diffraction

Powdered samples were analysed using a Bruker-AXS (Karlsruhe, Germany) D8 Advance Powder Diffractometer with a CuKα radiation source (40 kV, 40 mA) and a LynxEye detector. The scan ranged from 5° to 90° 2θ with a step size of 0.015° and a collection time of 0.7 s per step. Crystalline phases were identified by using the Search/Match algorithm, DIFFRAC.EVA 5.2 (Bruker-AXS, Karlsruhe, Germany) to search the International Center for Diffraction Data (ICDD) Powder Diffraction File (PDF4+ 2021 edition).

2.8. Elemental Analysis

Elemental analyses were performed at Source Certain International (Western Australia). Samples (0.25 g) were accurately weighed and digested in nitric acid (16 mL, 65 wt%), perchloric acid (4 mL, 70 wt%) and hydrofluoric acid (10 mL, 50 wt%) at approximately 180 °C for a minimum of 16 h under reflux. The acids were removed by evaporation at approximately 220 °C. Once the residue reached incipient dryness it was dissolved in hydrochloric acid (0.75 mL, 32 wt%), nitric acid (0.25 mL, 65 wt%) and DI-water (20 mL, >16.4 MΩ cm). The solution was suitably diluted for the instrument. Samples were analysed using an Agilent 5110 ICP-AES and an Agilent 7700 ICP-MS.

2.9. Total Organic Carbon

The rock samples were ground to a fine powder and digested in acid (HCl) to remove the carbonate minerals. The remaining residues were analysed using a LECO Carbon-Sulfur Analyser (CS-230). The CO2 produced was measured with an infra-red detector, and values calculated according to standard calibration.

3. Results and Discussion

Results are discussed in general terms of all three samples, unless a particular sample is specified. Samples subject to the same analysis showed generally consistent results as reflected in biomarker ratios presented in Table 1, Table 2, Table 3 and Table 4. Figures present the sample which best demonstrates compositional features.

3.1. Inorganic Composition

The coprolite components of samples PE 52316 and PE 52336 are approximately 2 cm and 3.5 cm in length, respectively and appear to be similar in composition (Figure 1). Three primary mineral regions were identified by X-ray diffraction (XRD). The coprolites are preserved in three dimensions as calcium phosphate with cracks filled with sphalerite. The concretionary material is siderite with minor amounts of quartz. Total organic carbon ranged from 0.31 to 0.78 wt% of rock in the fossils and 0.32 to 0.51 wt% in the matrices.
Phosphatic preservation is a characteristic component of carnivore coprolites, e.g., [1,4,59]. It has been suggested that rapid precipitation (within weeks—e.g., [14,19,60]) of dietary calcium phosphate can result in the preservation of fine morphological information [4,61]. Much like the formation of carbonate concretions, the remineralisation of calcium phosphate from faecal material has been demonstrated as autochthonous [7], occurring rapidly after deposition and prior to diagenetic permineralisation.
Elemental analyses revealed enrichment of rare earth elements (REE) in the coprolite fossils compared to the siderite concretion. Phosphatic minerals such as apatite are able to incorporate REE during early fossilisation [62] via substitution for calcium, e.g., [63]. REE are commonly used in palaeontological, palaeoenvironmental and palaeoredox studies, particularly those focused on vertebrate bones, e.g., [64,65,66,67,68] as well as on coprolites fossilised as apatite, e.g., [69,70].

3.2. Lipid Biomarkers of Coprolites

The coprolite fossil regions are characterised by a predominance of cholesterol-derived steroidal hydrocarbons (e.g., Figure 2 and Figure 3, Table 1). C27 cholestanes make up 86 to 99% of the total steranes composition of the aliphatic fraction (Table 1). Similarly, aromatic cholesteroids make up the majority of the monoaromatic and triaromatic steroid distributions of each of the fossils (Figure 3, Table 1). Cholesterol and its derivatives have been considered in past studies of faecal material as indicators of an animal diet, e.g., [28,29]. Cholesterol is generally known to be synthesised by animals, while ergosteroids can be found in fungi and some groups of algae [71,72] and sitosterol and stigmasterol are commonly made by higher plants [71]. As such, the relative proportions of related sterane biomarkers can help distinguish the contributions of these inputs in sedimentary organic matter input [73]. However, cholesteroids can also be produced by herbivores via modification of other sterol analogues, particularly phytosterols [74] or synthesised de novo, e.g., [74,75,76] and are also present in almost all eukaryotic cell walls [77]. While the exceptionally high abundance of cholesteroid biomarkers suggests a primarily animal diet for the coprolite producer, minor contributions from plant material as a dietary component or due to occasional grazing cannot be completely ruled out.
The diagenetic fate of biomarkers is typically controlled by factors such as burial depth, redox conditions, mineral and porewater chemistry and geological time [20]. The low diagenetic transformation of the cholesteroids in the coprolite specimens is likely related to favourable preservation conditions during aromatisation. Similarly immature biomolecular signals have been found in previous studies of carbonate concretions where rapid encapsulation, burial and mineralisation supported the preservation of soft tissues and intact biomolecules, e.g., [10,12,26].
Diagenetic processes promote rearrangement of ααα 20R stereoisomers to the more thermodynamically stable ααα 20S and αββ 20R + S compounds [20,25]. Within the cholestanes distribution the biologically derived cholestane 20R isomer is most dominant (e.g., Figure 2A). Additionally, present is C27 ααα 20S cholestane coeluting with coprostane (C27 βαα 20R). Here, the predominance of the C27 ααα 20R over the 20S isomer, (i.e., C27 ααα 20S/(20S + 20R) < 0.2), support a low diagenetic conversion of the steroids inside the fossil.
In contrast, the side chain isomerisation of C27 to C29 steranes as well as the transformation of ααα to αββ-steranes progressed much further in the concretion matrix (Figure 2). This indicates that diagenetic conversion of labile steroids was hampered by early cementation and the lack of clay mineral catalysis, e.g., [11,78,79] in the coprolite, allowing the labile steroids to persevere.
Due to coelution, the proportion of βαα 20R compared to ααα 20S was not precisely determined; however, it is presumed that a large area of this peak comprises C27 βαα 20R based on the low diagenetic conversion of steroids within the fossil and generally low maturity. This was confirmed by the mass spectrum which showed the presence of two coeluting peaks, one with m/z 149 fragment (C27 αββ 20S) and the other with m/z 151 fragment (C27 βαα 20R). Steroids with βαα stereochemistry can be produced by reduction of sterols in the intestine of mammals depending on the primary steroid components of their diet [27] or may also form in sediment through microbial reduction of ∆5 sterols [80].
2-sterenes and ∆3,5-steradienes are products formed in early diagenesis via dehydration of 5α(H)-stanols and ∆5-sterols, respectively [23,81]. ∆2-sterenes can also undergo isomerisation and rearrangement to ∆4- and ∆5-sterenes and subsequently to diasterenes [23,24,82,83]. A cholestadiene compound has been tentatively identified by comparison of its mass spectrum with that of a 3,5-cholestadiene [81], based on predominant fragments at m/z 213 and 368, in sample PE 52336. This compound is present in a higher relative abundance in the fossil than in the matrix, suggesting that it is likely derived from the original cholesterol content in the coprolite sample. Diasterenes in low abundance were also identified, although in the fossil of PE 52336 only, and are shown in the Supplementary Materials (Figure S1).
A- and B-ring monoaromatic steroids are derived directly from sterol compounds during early diagenesis, e.g., [84,85,86] while C-ring monoaromatic steroids have previously been attributed to later-stage diagenesis, e.g., [87], which can then be aromatised to form triaromatic hydrocarbons via loss of a methyl group, with increasing thermal maturity [88]. However, it has been suggested, e.g., [26,89] that aromatisation of sterene hydrocarbons to C-ring monoaromatic and triaromatic steroids can occur during early diagenesis via microbially mediated processes. Both the monoaromatic and triaromatic steroids show abundant cholesteroids compared to ergosteroids and stigmasteroids, with cholesteroids comprising between 55 to 82% of the total monoaromatic steroids distribution and 43 to 94% of total triaromatic steranes (Table 1). A predominance of the cholesteroid analogues in all steroid classes demonstrates that diagenetic transformation, while altering the cholesterol predominance somewhat, does not eliminate the primary dietary composition of steroids. Based on conventional models of triaromatic steroid formation, presence of triaromatic steroids in the fossil suggests a catagenetic history, which is inconsistent with the immaturity of the sample. It is therefore plausible that triaromatic steroids were also formed by microbially mediated reaction mechanisms, e.g., [26,89].
5α-cholestan-3β-ol was identified in the fossil of PE 52336 and PE 52316 (Table 2). This stanol may be either from a direct biological input or from early diagenesis by reduction of cholesterol [23,90]. The presence of 5α-cholestan-3β-ol in several of the fossils and its absence in any of the matrices supports that this is derived from the fossil coprolite itself. Coprostan-3-ol (5β-cholestan-3β-ol) was identified, albeit in low abundance in the PE 52315 fossil only (Table 2). This is a possible precursor sterol of coprostane that can be formed via reduction of cholesterol in the gut of many higher mammals, e.g., [27] as well as through microbial reduction in sediment [80]. Cholesterol was also identified in the fossil and matrix of all samples; however, this must be considered with caution, as it is also present in comparable concentrations in the procedural blank (Table 2).

3.3. Early Diagenetic Transformation of Dietary Sterols

A series of diagenetic products derived from cholesterol, such as cholestane and triaromatic steroids are preserved in all fossil coprolites. These components support a primarily carnivorous dietary source as they would support a high amount of cholesterol in the original faecal material, which has been rearranged and partially preserved. The cholesteroid compounds identified are summarised in Figure 6. Cholesterol can undergo reduction to 5α- and 5β-stanols [23], such as 5α-cholestan-3β-ol as present in the fossil coprolites. Further reduction would yield cholestane; specifically, the 5α-cholestane, which is the abundant sterane identified in each of the fossils. Cholesterol can also undergo reduction directly to Δ2-sterenes and Δ3,5-steradienes [23]. While sterenes were not observed, a single cholestadiene was identified in the PE 52336 fossil, which is formed from cholesterol and may be further converted into monoaromatic (Figure 3A) and triaromatic (Figure 3B) steroids.

3.4. Raman ChemoSpace

The Raman ChemoSpace is here used to complement compositional insights from the analysis of biomarkers with additional information on the molecular makeup and preserved biosignatures in the macromolecular fraction of Mazon Creek coprolites (FMNH PE 52316, FMNH PE 52336). Focus of the analysis are in situ Raman spectra collected for carbonaceous fossil remains. The discriminant analysis (CCA) of spectral data collected for all n = 37 samples, shown in Figure 4A, identifies key compositional characteristics of fossil vertebrate, non-annelid invertebrate, annelid, plant, and coprolite samples from the Mazon Creek locality. Among these samples, only minimal overlap between the clusters of the diverse array of fossil metazoan tissues is observed; the cluster of carbonaceous plant tissues is separated from the metazoan tissue clusters (Figure 4A). Using the spread of samples in the compositional space (=ChemoSpace), the circular fractions including all metazoan (orange) and plant (green) tissues were plotted in the CCA (Figure 4A,B). It is the nature of the CCA that separates different tissue categories (black vector arrows) radiating from the origin into the ChemoSpace, therefore allowing to constrain the distribution of samples through circular fractions. The trajectories of all tissue vectors and their associated, characteristic molecular features are plotted in Figure 5. Figure 4A locates the two coprolite samples within the circular fraction occupied by metazoan tissues, and reveals as characteristic molecular features a relatively increased abundance of alcohol (C-OH) and carbonyl groups (C=O), aromatics, inorganic phosphate (P-O, PO43-), and organo-sulfur moieties (C-S; Figure 5B). The majority of these characteristic molecular features that discriminate fossil coprolites from undigested fossil tissues coincide with the key chemical modifications of macromolecules experienced prior to fossilisation, during in vivo enzymatic digestion (Figure 5A): enzyme- and acid-catalysed (HCl(aq), stomach acid) hydrolysis cleaves primarily peptide bonds (amides), glucosidic bonds, esters (including phosphoesters), and breaks down lipids into their aliphatic and aromatic building blocks (please see steroid and hopanoid analyses in this study). Hydrolytic cleavage introduces alcohol groups, which can be converted into carbonyls during subsequent alteration. The molecular fingerprint of hydrolytic digestion preserves in the sampled Mazon Creek coprolites even after oxidative alteration (crosslinking reactions) acting during fossilisation [34,57]. Digestion moves data points in the CCA ChemoSpace along the vectors associated with functional groups introduced during hydrolytic cleavage (colored vectors in Figure 5B,C). The relative increase in organo-sulfur moieties is not related to digestive hydrolysis and thus relates either to the gastro-intestinal concentration processes or post-egestion alteration.
Figure 4B shows a complementary CCA ChemoSpace which, contrary to the analysis in Figure 4A, did not treat coprolites as a separate tissue category, but instead allowed them to group with the other n = 4 tissue categories based on the affinity of the contained, digested remains. Both coprolite samples plot within the circular fraction occupied by metazoan samples, and one tangentially overlaps with the vertebrate tissue cluster. Considering that digestion shifts samples in the ChemoSpace CCE (shown in Figure 4A and Figure 5B,C) towards the right, we thus infer a metazoan affinity of the digested faecal matter and identify the coprolite producer as a carnivore rather than omnivore or herbivore (no primary contribution of plants to the diet).
The carnivorous diet, potentially based on vertebrates (as suggested by the proximity to the vertebrate cluster in Figure 4A,B), suggests that the producer of the coprolites represented a high trophic level in the ancient Mazon Creek ecosystem and was potentially an apex predator. Future experimental work on the ChemoSpace effects of different metazoan digestive processes has the potential to allow for the detailed quantification of the contributions of vertebrate, plant, non-annelid invertebrate, and annelid food items to the diet of coprolite producers.
The independent data collected for the soluble organic phase (biomarker analysis, stable isotope fractionation) and the insoluble organic phase (ChemoSpace of carbonaceous fossils) converge in their result: a carnivorous diet is inferred for the producers of the analysed coprolites.

3.5. Lipid Biomarkers of Matrix (Palaeoenvironmental Signal)

Terpane and hopane compounds were more dominant relative to steranes in the concretion matrix compared to the coprolite, as reflected in the regular sterane/hopane ratios (1.9–80 in the fossils compared to 0.54–1.1 in matrices) (Table 3). The Ts/(Ts + Tm) ratio, traditionally considered to be indicative of thermal maturity and clay catalysis, indicates the end of the diagenetic or the onset of the catagenetic stage of thermal transformation for both the coprolite and concretion matrix, which stands in contrast to the presence of a thermally labile steradiene and stanol (Table 3). C29 20S/(20S + 20R) steranes ratios (Table 3) are also low overall but consistently higher in the fossil than in the matrix. The Ts/(Ts + Tm) ratios are consistent between the fossil and matrix while the diasterane/regular steranes ratios vary, supporting that different mechanisms are responsible for formation of trisnorneohopane than for the isomerization of steranes. Dia/regular sterane ratios for all steroid homologs show a higher proportions of rearranged steroids in the matrix than in the fossil. In general, even in the matrix the C27 dia/regular sterane ratio indicates lower maturity; however, the C27 sterane ratios are likely strongly influenced by source input.

3.6. Stable Carbon Isotopes

Carbon isotopes of the bulk organic matter were measured on the residue of extracted sediment from samples PE 52316 and PE 52336 after treatment with hydrochloric acid to remove carbonate. Both fossil and matrix of each sample showed δ13C values around −23.8‰ (Table 4).
Compound specific carbon isotopes were measured for the n-alkanes of these two samples, as well as for the C27 ααα 20S sterane in the fossil portions (Table 3). δ13C values of cholestane (−32.8‰) were depleted by approximately 2.0‰ compared to the average isotopic values of n-C20–25 (−30.8‰) and n-C18 (−31.1‰), and by approximately 9.0‰ compared to the bulk organic matter δ13C values.
Sterane hydrocarbons are typically enriched compared to linear hydrocarbons from the same source by up to 8‰ [93]. Cholestane is here interpreted as derived from dietary cholesterol from the coprolite, while slightly more 13C enriched n-alkanes represent input from external sources, most likely palaeoenvironmental signals. The depletion of δ13C value of cholestane in the two coprolites by approximately 2‰ compared to the average of the abundant mid-chain n-alkanes is therefore consistent with input from two different sources. The isotopic consistency (approximately −30.8‰) of the n-alkanes in the fossil and the matrix is indicative of being derived from a common source, while the depletion of cholestane suggests it has a source different to the n-alkanes. The n-alkanes in the samples are therefore likely to originate from freshwater producers (e.g., algae, aquatic macrophytes, mosses), while cholesteroids intrinsic to the coprolite derive from animal sources.
Phytol is also typically 13C enriched with respect to fatty acids synthesised in phytoplanktonic cells by 2–5‰ [93]. A similar 2–5‰ distinction between phytane and straight-chain hydrocarbons (and pristane) would therefore be expected of phytane derived from phytol, synthesised by chlorophyll a [93]. In these samples phytane was depleted in 13C by approximately 2.5–3.5‰ compared to n-C18 and 2.5–4‰ compared to mid-chain n-alkanes (Table 3) which is more typical of an origin within a methane cycle, wherein phytane is derived from ether lipids of methanotrophs and not phytol, e.g., [94,95]. Methanogenic archaea and methanotrophs are considered important components of the microbial growth mechanisms of siderite concretions, e.g., [95].

4. Conclusions

This study demonstrates that:
(1)
Cholesteroids including intact 5α-cholestan-3β-ol and coprostanol have been preserved in siderite concretions hosting 306 million-year-old coprolites.
(2)
The molecular data obtained by GC-MS, GC-MRM, GC-irMS and Raman microspectroscopy supports a primarily carnivorous diet and suggest an elevated trophic position for the coprolite producer.
(3)
The preservation of intact dietary sterols and macromolecular biosignatures is attributed to rapid encapsulation of the coprolites within days to months after egestion.
(4)
Siderite (FeCO3) concretions (of Carboniferous age) seem to preserve intact and modified, but not unrecognisable, biomolecules much like calcium carbonate concretions of Jurassic and Devonian ages.
The results demonstrate that molecular information preserved within fossils can provide important ancient dietary insights either alongside or independent of traditional mineralogical or morphological studies of coprolites. Intact dietary sterols present in fossils support the significance of rapid encapsulation and organo-templated growth of carbonate concretions in the preservation of biomolecules in geological time. Carbonate concretions, which are host to soft-tissue fossil preservation are evidently important samples for molecular studies and represent each a unique opportunity to study extinct species and past environments.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology11091289/s1, Figure S1: GC-MS m/z 257 chromatogram of aliphatic fraction showing presence of diasterenes. A: 10α, Δ13(17) diacholestene 20S; B: 10α, Δ13(17) diacholestene 20R; C: 10α, Δ13(17) 24-methyldiacholestene 20S (*Tentatively identified based on location of R-isomer and mass spectrum); D:10α, Δ13(17) 24-methyldiacholestene 20R; E: 10α, Δ13(17) 24-ethyldiacholestene 20S; F: 10α, Δ13(17) 24-ethyldiacholestene 20R. Compounds were tentatively identified based on comparison of retention time and peak patterns with sample where compounds had been previously identified (e.g., [26]), confirmed using mass spectrum of peaks.; Table S1: A list of all specimens analyzed for the Figure 4 and Figure 5, including catalogue numbers, institutional and museum identifiers, and sampled tissue types. In addition, raw ChemoSpace Canonical Correspondence Analyses (CCA) are plotted (in PAST 3) for the Figure 4 and Figure 5 (separate plots: (1) specimen data points and tissue trajectories for Figure 4, and (2) corresponding discriminant functional group vectors in the CCA space for Figure 5).

Author Contributions

M.T. and K.G. designed the experiments and overall project concept. M.T. performed all laboratory on the fossil concretions except Raman microspectroscopy and GC-MRM. J.W. performed Raman microspectroscopy and J.B. performed GC-MRM. M.T. wrote the manuscript with contributions from all co-authors. K.G. supervised PhD scholar M.T. and P.M. provided Mazon Creek samples and hosted K.G. and M.T. at The Field Museum for this study. L.S. contributed with interpretive work and written contributions. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Australian Research Council (ARC) for an ARC-Laureate Fellowship grant (FL210100103) and ARC infrastructure grants (LE110100119; LE100100041; LE0882836; LE0668345; LE0775551).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Source data for ChemoSpace analyses (Figure 4 and Figure 5) are available in the Supplementary Materials. All other data can be made available upon request from corresponding authors.

Acknowledgments

We thank Peter Hopper, Alex Holman and Janet Hope for their technical support with GC-MS and stable isotope analyses. Tripp thanks Curtin University for a Research Training Postgraduate award, The Australian Institute of Nuclear Science and Engineering for an AINSE Postgraduate Research Award and support from The Institute for Geoscience Research for Star-fish Soxhlet extraction apparatus. Part of this research was undertaken using the XRD instrumentation at the John de Laeter Centre, Curtin University. Tripp thanks Veronica Avery and Matthew Rowles for technical assistance with XRD analysis. We thank Scott Lidgard (Field Museum, Chicago) for providing samples. We thank the three anonymous reviewers for their constructive comments that helped improve this manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Hollocher, K.T.; Alcober, O.A.; Colombi, C.E.; Hollocher, T.C. Carnivore Coprolites from the Upper Triassic Ischigualasto Formation, Argentina: Chemistry, Mineralogy, and Evidence for Rapid Initial Mineralization. Palaios 2005, 20, 51–63. [Google Scholar] [CrossRef]
  2. Chin, K. The Paleobiological Implications of Herbivorous Dinosaur Coprolites from the Upper Cretaceous Two Medicine Formation of Montana: Why Eat Wood? Palaios 2007, 22, 554–566. [Google Scholar] [CrossRef]
  3. Gill, F.L.; Crump, M.P.; Schouten, R.; Bull, I.D. Lipid analysis of a ground sloth coprolite. Quat. Res. 2009, 72, 284–288. [Google Scholar] [CrossRef]
  4. Hollocher, K.T.; Hollocher, T.C.; Rigby, J.K. A Phosphatic Coprolite Lacking Diagenetic Permineralization from the Upper Cretaceous Hell Creek Formation, Northeastern Montana: Importance of Dietary Calcium Phosphate in Preservation. Palaios 2010, 25, 132–140. [Google Scholar] [CrossRef]
  5. Khosla, A.; Chin, K.; Alimohammadin, H.; Dutta, D. Ostracods, plant tissues, and other inclusions in coprolites from the Late Cretaceous Lameta Formation at Pisdura, India: Taphonomical and palaeoecological implications. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2015, 418, 90–100. [Google Scholar] [CrossRef]
  6. Northwood, C. Early Triassic coprolites from Australia and their palaeobiological significance. Palaeontology 2005, 48, 49–68. [Google Scholar] [CrossRef]
  7. Bajdek, P.; Owocki, K.; Niedźwiedzki, G. Putative dicynodont coprolites from the Upper Triassic of Poland. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2014, 411, 1–17. [Google Scholar] [CrossRef]
  8. Boast, A.P.; Weyrich, L.S.; Wood, J.R.; Metcalf, J.L.; Knight, R.; Cooper, A. Coprolites reveal ecological interactions lost with the extinction of New Zealand birds. Proc. Natl. Acad. Sci. USA 2018, 115, 1546–1551. [Google Scholar] [CrossRef]
  9. Witt, K.E.; Yarlagadda, K.; Allen, J.M.; Bader, A.C.; Simon, M.L.; Kuehn, S.R.; Swanson, K.S.; Cross, T.-W.L.; Hedman, K.M.; Ambrose, S.H.; et al. Integrative analysis of DNA, macroscopic remains and stable isotopes of dog coprolites to reconstruct community diet. Sci. Rep. 2021, 11, 3113. [Google Scholar] [CrossRef]
  10. Lengger, S.K.; Melendez, I.M.; Summons, R.E.; Grice, K. Mudstones and embedded concretions show differences in lithology-related, but not source-related biomarker distributions. Org. Geochem. 2017, 113, 67–74. [Google Scholar] [CrossRef] [Green Version]
  11. Melendez, I.; Grice, K.; Trinajstic, K.; Ladjavardi, M.; Greenwood, P.; Thompson, K. Biomarkers reveal the role of photic zone euxinia in exceptional fossil preservation: An organic geochemical perspective. Geology 2013, 41, 123–126. [Google Scholar] [CrossRef]
  12. Plet, C.; Grice, K.; Pagès, A.; Verrall, M.; Coolen, M.J.L.; Ruebsam, W.; Rickard, W.D.A.; Schwark, L. Palaeobiology of red and white blood cell-like structures, collagen and cholesterol in an ichthyosaur bone. Sci. Rep. 2017, 7, 13776. [Google Scholar] [CrossRef] [PubMed]
  13. Allison, P.A. Konservat-Lagerstatten: Cause and Classification. Paleobiology 1988, 14, 331–344. [Google Scholar] [CrossRef]
  14. Briggs, D.E.G. The role of decay and mineralization in the preservation of soft-bodied fossils. Annu. Rev. Earth Planet. Sci. 2003, 31, 275–301. [Google Scholar] [CrossRef]
  15. Parry, L.A.; Smithwick, F.; Nordén, K.K.; Saitta, E.T.; Lozano-Fernandez, J.; Tanner, A.R.; Caron, J.B.; Edgecombe, G.D.; Briggs, D.E.; Vinther, J. Soft-Bodied Fossils Are Not Simply Rotten Carcasses—Toward a Holistic Understanding of Exceptional Fossil Preservation. Bioessays 2018, 40, 1700167. [Google Scholar] [CrossRef]
  16. Seilacher, A. Begriff und bedeutung der Fossil-Lagerstätten: Neues Jarhbuch für Geologie und Paläontologie. Neues Jahrb. Geol. Paläontologie Mon. 1970, 34–39. [Google Scholar]
  17. Grice, K.; Holman, A.I.; Plet, C.; Tripp, M. Fossilised Biomolecules and Biomarkers in Carbonate Concretions from Konservat-Lagerstätten. Minerals 2019, 9, 158. [Google Scholar] [CrossRef]
  18. Mccoy, V.E. Concretions as agents of soft-tissue preservation: A review. Paleontol. Soc. Pap. 2014, 20, 147–162. [Google Scholar] [CrossRef]
  19. Sagemann, J.; Bale, S.J.; Briggs, D.E.; Parkes, R.J. Controls on the formation of authigenic minerals in association with decaying organic matter: An experimental approach. Geochim. Cosmochim. Acta 1999, 63, 1083–1095. [Google Scholar]
  20. Peters, K.E.; Walters, C.C.; Moldowan, J.M. Volume I: Biomarkers and Isotopes in the Environment and Human History. In The Biomarker Guide; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  21. Volkman, J.K. A review of sterol markers for marine and terrigenous organic matter. Org. Geochem. 1986, 9, 83–99. [Google Scholar] [CrossRef]
  22. Huang, W.Y.; Meinschein, W.G. Sterols as source indicators of organic materials in sediments. Geochim. Cosmochim. Acta 1976, 40, 323–330. [Google Scholar]
  23. Mackenzie, A.S.; Brassel, S.C.; Eglinton, G.; Maxwell, J.R. Chemical fossils: The geological fate of steroids. Science 1982, 217, 491–504. [Google Scholar] [CrossRef] [PubMed]
  24. Brassell, S.C.; McEvoy, J.; Hoffmann, C.F.; Lamb, N.A.; Peakman, T.M.; Maxwell, J.R. Isomerisation, rearrangement and aromatisation of steroids in distinguishing early stages of diagenesis. Org. Geochem. 1984, 6, 11–23. [Google Scholar] [CrossRef]
  25. Mackenzie, A.S.; Lamb, N.A.; Maxwell, J.R. Steroid hydrocarbons and the thermal history of sediments. Nature 1982, 295, 223–226. [Google Scholar] [CrossRef]
  26. Melendez, I.; Grice, K.; Schwark, L. Exceptional preservation of Palaeozoic steroids in a diagenetic continuum. Sci. Rep. 2013, 3, 2768. [Google Scholar] [CrossRef]
  27. Bull, I.D.; Lockheart, M.J.; Elhmmali, M.M.; Roberts, D.J.; Evershed, R.P. The origin of faeces by means of biomarker detection. Environ. Int. 2002, 27, 647–654. [Google Scholar] [CrossRef]
  28. Gill, F.L.; Bull, I.D. Lipid analysis of vertebrate coprolites. N. M. Mus. Nat. Hist. Sci. Bull. 2012, 57, 93–98. [Google Scholar]
  29. Umamaheswaran, R.; Prasad, G.V.R.; Rudra, A.; Dutta, S. Biomarker Signatures in Triassic Coprolites. Palaios 2019, 34, 458–467. [Google Scholar] [CrossRef]
  30. Weber, D.J.; Lawler, G.C. Lipid Components of the Coprolites. Geology, Paleontology and Paleoecology of a Late Triassic Lake, Western New Mexico. Brigh. Young Univ. Geol. Stud. 1978, 25, 75–87. [Google Scholar]
  31. Zatoń, M.; Niedźwiedzki, G.; Marynowski, L.; Benzerara, K.; Pott, C.; Cosmidis, J.; Krzykawski, T.; Filipiak, P. Coprolites of Late Triassic carnivorous vertebrates from Poland: An integrative approach. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2015, 430, 21–46. [Google Scholar] [CrossRef]
  32. Harrault, L.; Milek, K.; Jardé, E.; Jeanneau, L.; Derrien, M.; Anderson, D.G. Faecal biomarkers can distinguish specific mammalian species in modern and past environments. PLoS ONE 2019, 14, e0211119. [Google Scholar] [CrossRef] [PubMed]
  33. Leeming, R.; Ball, A.; Ashbolt, N.; Nichols, P. Using faecal sterols from humans and animals to distinguish faecal pollution in receiving waters. Water Res. 1996, 30, 2893–2900. [Google Scholar] [CrossRef]
  34. McCoy, V.E.; Wiemann, J.; Lamsdell, J.C.; Whalen, C.D.; Lidgard, S.; Mayer, P.; Petermann, H.; Briggs, D.E. Chemical signatures of soft tissues distinguish between vertebrates and invertebrates from the Carboniferous Mazon Creek Lagerstatte of Illinois. Geobiology 2020, 18, 560–565. [Google Scholar] [CrossRef]
  35. Sistiaga, A.; Berna, F.; Laursen, R.; Goldberg, P. Steroidal biomarker analysis of a 14,000 years old putative human coprolite from Paisley Cave, Oregon. J. Archaeol. Sci. 2014, 41, 813–817. [Google Scholar] [CrossRef]
  36. Clements, T.; Purnell, M.; Gabbott, S. The Mazon Creek Lagerstätte: A diverse late Paleozoic ecosystem entombed within siderite concretions. J. Geol. Soc. 2018, 176, 1–11. [Google Scholar] [CrossRef] [Green Version]
  37. Wittry, J. The Mazon Creek Fossil Fauna; ESCONI: Downers Grove, IL, USA, 2012. [Google Scholar]
  38. Peppers, R.A. Palynological Correlation of Major Pennsylvanian (Middle and Upper Carboniferous) Chronostratigraphic Boundaries in the Illinois and Other Coal Basins; The Geological Society of America, Inc.: Boulder, CO, USA, 1996. [Google Scholar]
  39. Pfefferkorn, H.W. High diversity and stratigraphic age of the Mazon Creek flora. In Mazon Creek Fossils; Nitecki, M.H., Ed.; Academic Press, Inc.: New York, NY, USA, 1979; pp. 129–142. [Google Scholar]
  40. Baird, G.C. Lithology and fossil distribution, Francis Creek Shale in northeastern Illinois. In Mazon Creek Fossils; Nitecki, M.H., Ed.; Academic Press, Inc.: New York, NY, USA, 1979; pp. 41–67. [Google Scholar]
  41. Baird, G.C.; Shabica, C.W.; Anderson, J.L.; Richardson, E.S. Biota of a Pennsylvanian muddy coast: Habitats within the Mazonian delta complex, northeast Illinois. J. Paleontol. 1985, 59, 253–281. [Google Scholar]
  42. Baird, G.C.; Sroka, S.D.; Shabica, C.W.; Keucher, G.J. Taphonomy of Middle Pennsylvanian Mazon Creek area fossil localities, Northeast Illinois; significance of exceptional fossil preservation in syngenetic concretions. Palaios 1986, 1, 271–285. [Google Scholar] [CrossRef]
  43. Cotroneo, S.; Schiffbauer, J.D.; McCoy, V.E.; Wortmann, U.G.; Darroch, S.A.F.; Peng, Y.; Laflamme, M. A new model of the formation of Pennsylvanian iron carbonate concretions hosting exceptional soft-bodied fossils in Mazon Creek, Illinois. Geobiology 2016, 14, 543–555. [Google Scholar] [CrossRef]
  44. Shabica, C.W. Pennsylvanian sedimentation in Northern Illinois: Examination of delta models. In Mazon Creek Fossils; Nitecki, M.H., Ed.; Academic Press, Inc.: New York, NY, USA, 1979. [Google Scholar]
  45. Baird, G.C.; Sroka, S.D.; Shabica, C.W.; Beard, T.L. Mazon Creek-type fossil assemblages in the U.S. midcontinent Pennsylvanian: Their recurrent character and palaeoenvironmental significance. Philos. Trans. R. Soc. Lond. B Biol. Sci. 1985, 311, 87–99. [Google Scholar]
  46. Berner, R.A. A new geochemical classification of sedimentary environments. J. Sediment. Petrol. 1981, 51, 359–365. [Google Scholar]
  47. Berner, R.A. Sulphate reduction, organic matter decomposition and pyrite formation. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci. 1985, 315, 25–38. [Google Scholar]
  48. Claypool, G.E.; Kaplan, I.R. The Origin and Distribution of Methane in Marine Sediments. In Natural Gases in Marine Sediments; Maine Science; Kaplan, I.R., Ed.; Springer: Boston, MA, USA, 1974; Volume 3, pp. 99–139. [Google Scholar]
  49. Allison, P.A. The role of anoxia in the decay and mineralization of proteinaceous macro-fossils. Paleobiology 1988, 14, 139–154. [Google Scholar] [CrossRef]
  50. Berner, R.A. Rate of concretion growth. Geochim. Cosmochim. Acta 1968, 32, 477–483. [Google Scholar] [CrossRef]
  51. Briggs, D.E.; Kear, A.J. Decay and preservation of polychaetes: Taphonomic thresholds in soft-bodied organisms. Paleobiology 1993, 19, 107–135. [Google Scholar] [CrossRef]
  52. McCoy, V.E.; Young, R.T.; Briggs, D.E.G. Sediment Permeability and the Preservation of Soft-Tissues in Concretions: An Experimental Study. Palaios 2015, 30, 608–612. [Google Scholar] [CrossRef]
  53. Yoshida, H.; Ujihara, A.; Minami, M.; Asahara, Y.; Katsuta, N.; Yamamoto, K.; Sirono, S.; Maruyama, I.; Nishimoto, S.; Metcalfe, R. Early post-mortem formation of carbonate concretions around tusk-shells over week-month timescales. Sci. Rep. 2015, 5, 14123. [Google Scholar] [CrossRef]
  54. Yoshida, H.; Yamamoto, K.; Minami, M.; Katsuta, N.; Sin-ichi, S.; Metcalfe, R. Generalized conditions of spherical carbonate formation around decaying organic matter in early diagenesis. Sci. Rep. 2018, 8, 6308. [Google Scholar] [CrossRef]
  55. Brocks, J.J.; Hope, J.M. Tailing of chromatographic peaks in GC-MS caused by interaction of halogenated solvents with the ion source. J. Chromatogr. Sci. 2014, 52, 471–475. [Google Scholar] [CrossRef]
  56. Coplen, T.B.; Brand, W.A.; Gehre, M.; Gröning, M.; Meijer, H.A.J.; Toman, B.; Verkouteren, R.M. New Guidelines for δ13C Measurements. Anal. Chem. 2006, 78, 2439–2441. [Google Scholar] [CrossRef] [Green Version]
  57. Wiemann, J.; Crawford, J.M.; Briggs, D.E.G. Phylogenetic and physiological signals in metazoan fossil biomolecules. Sci. Adv. 2020, 6, eaba6883. [Google Scholar] [CrossRef] [PubMed]
  58. Lambert, J.B.; Shurvell, H.F.; Lightner, D.A.; Cooks, R.G.; Stout, G.H. Introduction to Organic Spectroscopy; Macmillan: New York, NY, USA, 1987. [Google Scholar]
  59. Chin, K.; Tokaryk, T.T.; Erickson, G.M.; Calk, L.C. A king-sized theropod coprolite. Nature 1998, 393, 680–682. [Google Scholar] [CrossRef]
  60. Briggs, D.E.; Wilby, P.R. The role of the calcium carbonate-calcium phosphate switch in the mineralization of soft-bodied fossils. J. Geol. Soc. 1996, 153, 665–668. [Google Scholar] [CrossRef]
  61. Briggs, D.E.; Kear, A.J.; Martill, D.M.; Wilby, P.R. Phosphatization of soft-tissue in experiments and fossils. J. Geol. Soc. 1993, 150, 1035–1038. [Google Scholar] [CrossRef]
  62. Trueman, C.N. Rare Earth Element Geochemistry and Taphonomy of Terrestrial Vertebrate Assemblages. Palaios 1999, 14, 555–568. [Google Scholar] [CrossRef]
  63. Fleet, M.E.; Pan, Y. Site preference of rare earth elements in fluorapatite. Am. Mineral. 1995, 80, 329–335. [Google Scholar] [CrossRef]
  64. Suarez, C.A.; Macpherson, G.L.; González, L.A.; Grandstaff, D.E. Heterogeneous rare earth element (REE) patterns and concentrations in a fossil bone: Implications for the use of REE in vertebrate taphonomy and fossilization history. Geochim. Cosmochim. Acta 2010, 74, 2970–2988. [Google Scholar] [CrossRef]
  65. Trueman, C.N.; Benton, M.J.; Palmer, M.R. Geochemical taphonomy of shallow marine vertebrate assemblages. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2003, 197, 151–169. [Google Scholar] [CrossRef]
  66. Trueman, C.N.; Behrensmeyer, A.K.; Potts, R.; Tuross, N. High-resolution records of location and stratigraphic provenance from the rare earth element composition of fossil bones. Geochim. Cosmochim. Acta 2006, 70, 4343–4355. [Google Scholar] [CrossRef]
  67. Trueman, C.N.; Palmer, M.R.; Field, J.; Privat, K.; Ludgate, N.; Chavagnac, V.; Eberth, D.A.; Cifelli, R.; Rogers, R.R. Comparing rates of recrystallisation and the potential for preservation of biomolecules from the distribution of trace elements in fossil bones. Comptes Rendus Palevol 2008, 7, 145–158. [Google Scholar] [CrossRef]
  68. Trueman, C.N.; Tuross, N. Trace Elements in Recent and Fossil Bone Apatite. Rev. Mineral. Geochem. 2002, 48, 489–521. [Google Scholar] [CrossRef]
  69. Niedźwiedzki, G.; Bajdek, P.; Owocki, K.; Kear, B.P. An Early Triassic polar predator ecosystem revealed by vertebrate coprolites from the Bulgo Sandstone (Sydney Basin) of southeastern Australia. Palaeogeogr. Palaeoclimatol. Palaeoecol. 2016, 464, 5–15. [Google Scholar] [CrossRef]
  70. Owocki, K.; Niedźwiedzki, G.; Sennikov, A.G.; Golubev, V.K.; Janiszewska, K.; Sulej, T. Upper Permian vertebrate coprolites from Vyazniki and Gorokhovets, Vyatkian Regional Stage, Russian Platform. Palaios 2012, 27, 867–877. [Google Scholar] [CrossRef]
  71. Brocks, J.J.; Grice, K. Biomarkers (Molecular Fossils). In Encyclopedia of Geobiology; Reitner, J., Thiel, V., Eds.; Springer: Dordrecht, The Netherlands, 2011; pp. 147–167. [Google Scholar]
  72. Volkman, J.K. Sterols in microorganisms. Appl. Microbiol. Biotechnol. 2003, 60, 495–506. [Google Scholar] [CrossRef] [PubMed]
  73. Peters, K.E.; Walters, C.C.; Moldowan, J.M. Volume 2: Biomarkers and Isotopes in Petroleum Exploration and Earth History. In The Biomarker Guide; Cambridge University Press: Cambridge, UK, 2005. [Google Scholar]
  74. Grice, K.; Klein Breteler, W.C.M.; Schouten, S.; Grossi, V.; de Leeuw, J.W.; Damsté, J.S.S. Effects of zooplankton herbivory on biomarker proxy records. Paleoceanography 1998, 13, 686–693. [Google Scholar] [CrossRef]
  75. Goad, L.J. Sterol biosynthesis and metabolism in marine invertebrates. Pure Appl. Chem. 1981, 53, 837–852. [Google Scholar] [CrossRef]
  76. Kanazawa, A. Sterols in marine invertabrates. Fish. Sci. 2001, 67, 997–1007. [Google Scholar] [CrossRef]
  77. Prost, K.; Birk, J.J.; Lehndorff, E.; Gerlach, R.; Amelung, W. Steroid Biomarkers Revisited—Improved Source Identification of Faecal Remains in Archaeological Soil Material. PLoS ONE 2017, 12, e0164882. [Google Scholar] [CrossRef]
  78. Nabbefeld, B.; Grice, K.; Schimmelmann, A.; Summons, R.E.; Troitzsch, U.; Twitchett, R.J. A comparison of thermal maturity parameters between freely extracted hydrocarbons (Bitumen I) and a second extract (Bitumen II) from within the kerogen matrix of Permian and Triassic sedimentary rocks. Org. Geochem. 2010, 41, 78–87. [Google Scholar] [CrossRef]
  79. Van Kaam-Peters, H.M.E.; Köster, J.; Van Der Gaast, S.J.; Dekker, M.; De Leeuw, J.W.; Sinninge Damsté, J.S. The effect of clay minerals on diasterane/sterane ratios. Geochim. Cosmochim. Acta 1998, 62, 2923–2929. [Google Scholar] [CrossRef]
  80. Peakman, T.M.; De Leeuw, J.W.; Rijpstra, C. Identification and origin of Δ8(14)5α- and Δ145α-sterenes and related hydrocarbons in an immature bitumen from the Monterey Formation, California. Geochim. Cosmochim. Acta 1992, 56, 1223–1230. [Google Scholar] [CrossRef]
  81. Gagosian, R.B.; Farrington, J.W. Sterenes in surface sediments from the southwest African shelf and slope. Geochim. Cosmochim. Acta 1978, 42, 1091–1101. [Google Scholar] [CrossRef]
  82. Brassell, S.C.; Murchison, D.G.; Mason, R.; Durand, B.; Eglinton, G.; Comet, P.A.; Curtis, C.D.; Bada, J.; de Leeuw, J.W. Molecular Changes in Sediment Lipids as Indicators of Systematic Early Diagenesis [and Discussion]. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Sci. 1985, 315, 57–75. [Google Scholar]
  83. Rubinstein, I.; Sieskind, O.; Albrecht, P. Rearranged sterenes in a shale: Occurrence and simulated formation. J. Chem. Soc. Perkin Trans. 1 1975, 19, 1833–1836. [Google Scholar] [CrossRef]
  84. Hussler, G.; Albrecht, P. C27–C29 Monoaromatic anthrasteroid hydrocarbons in Cretaceous black shales. Nature 1983, 304, 262–263. [Google Scholar] [CrossRef]
  85. Hussler, G.; Chappe, B.; Wehrung, P.; Albrecht, P. C27–C29 ring A monoaromatic steroids in Cretaceous black shales. Nature 1981, 294, 556–558. [Google Scholar] [CrossRef]
  86. Schüpfer, P.; Finck, Y.; Houot, F.; Gülaçar, F.O. Acid catalysed backbone rearrangement of cholesta-2,4,6-triene: On the origin of ring A and ring B aromatic steroids in recent sediments. Org. Geochem. 2007, 38, 671–681. [Google Scholar] [CrossRef]
  87. Riolo, J.; Hussler, G.; Albrecht, P.; Connan, J. Distribution of aromatic steroids in geological samples: Their evaluation as geochemical parameters. Org. Geochem. 1986, 10, 981–990. [Google Scholar] [CrossRef]
  88. Mackenzie, A.S. Organic Reactions as Indicators of the Burial and Temperature Histories of Sedimentary Sequences. Clay Miner. 1984, 19, 271–286. [Google Scholar] [CrossRef]
  89. Ling, Y.C.; Moreau, J.; Berwick, L.; Tulipani, S.; Grice, K.; Bush, R. Distribution of iron- and sulfate-reducing bacteria across a coastal acid sulfate soil (CASS) environment: Implications for passive bioremediation by tidal inundation. Front. Microbiol. 2015, 6, 624. [Google Scholar] [CrossRef]
  90. Gaskell, S.J.; Eglinton, G. Sterols of a contemporary lacustrine sediment. Geochim. Cosmochim. Acta 1976, 40, 1221–1228. [Google Scholar] [CrossRef]
  91. Cranwell, P.A. Lipid geochemistry of sediments from Upton Broad, a small productive lake. Org. Geochem. 1984, 7, 25–37. [Google Scholar] [CrossRef]
  92. Ficken, K.J.; Li, B.; Swain, D.L.; Eglinton, G. An n-alkane proxy for the sedimentary input of submerged/floating freshwater aquatic macrophytes. Org. Geochem. 2000, 31, 745–749. [Google Scholar] [CrossRef]
  93. Schouten, S.; Klein Breteler, W.C.M.; Blokker, P.; Schogt, N.; Rijpstra, W.I.C.; Grice, K.; Baas, M.; Damsté, J.S.S. Biosynthetic effects on the stable carbon isotopic compositions of algal lipids: Implications for deciphering the carbon isotopic biomarker record. Geochim. Cosmochim. Acta 1998, 62, 1397–1406. [Google Scholar] [CrossRef]
  94. Summons, R.E.; Jahnke, L.L.; Roksandic, Z. Carbon isotopic fractionation in lipids from methanotrophic bacteria: Relevance for interpretation of the geochemical record of biomarkers. Geochim. Cosmochim. Acta 1994, 58, 2853–2863. [Google Scholar] [CrossRef]
  95. Bojanowski, M.J.; Clarkson, E.N.K. Origin of Siderite Concretions in Microenvironments of Methanogenesis Developed in a Sulfate Reduction Zone: An Exception or a Rule? J. Sediment. Res. 2012, 82, 585–598. [Google Scholar] [CrossRef]
Figure 1. Images of the three samples used in this study. Regions are defined as the ‘Fossil’ and the ‘Matrix’ corresponding to the coprolite fossil region and the concretionary region, respectively, as demonstrated by white dashed lines. Subsamples of each region were used for geochemical analysis. ‘Matrix’ is used throughout to refer to the concretionary region, which is assumed to consist of primarily concretionary material plus potential external organic matter inputs.
Figure 1. Images of the three samples used in this study. Regions are defined as the ‘Fossil’ and the ‘Matrix’ corresponding to the coprolite fossil region and the concretionary region, respectively, as demonstrated by white dashed lines. Subsamples of each region were used for geochemical analysis. ‘Matrix’ is used throughout to refer to the concretionary region, which is assumed to consist of primarily concretionary material plus potential external organic matter inputs.
Biology 11 01289 g001
Figure 2. Metastable Reaction Monitoring (MRM) chromatograms of M+ → 217 precursor to product transitions of C27–29 steranes of PE 52316 Fossil (A) and Matrix (B). Peaks are coloured according to the number of carbon atoms in the sterane. α and β nomenclature refers to the stereochemistry of hydrogen at C-5, C-14 and C-17 for regular steranes and C-13 and C-17 for diasteranes, while S/R refers to stereochemistry at the C-20 position. Percentages represent the relative abundance of the most abundant peak in each transition.
Figure 2. Metastable Reaction Monitoring (MRM) chromatograms of M+ → 217 precursor to product transitions of C27–29 steranes of PE 52316 Fossil (A) and Matrix (B). Peaks are coloured according to the number of carbon atoms in the sterane. α and β nomenclature refers to the stereochemistry of hydrogen at C-5, C-14 and C-17 for regular steranes and C-13 and C-17 for diasteranes, while S/R refers to stereochemistry at the C-20 position. Percentages represent the relative abundance of the most abundant peak in each transition.
Biology 11 01289 g002
Figure 3. (A) Selected Ion Recording (SIR) m/z = 253 chromatogram of C27–29 monoaromatic steroids of PE 52336 Fossil (i) and Matrix (ii). I = 5β(H),10β(CH3); II = 5α(H),10β(CH3); V = 5β(CH3),10β(H); VII = 5α(CH3),10α(H). Percentages represent the relative abundance of the most abundant peak in each transition. (B) SIM m/z = 231 chromatogram of C26–28 triaromatic steroids of PE 52336 Fossil (i) and Matrix (ii). S/R refers to the stereochemistry at C-20.
Figure 3. (A) Selected Ion Recording (SIR) m/z = 253 chromatogram of C27–29 monoaromatic steroids of PE 52336 Fossil (i) and Matrix (ii). I = 5β(H),10β(CH3); II = 5α(H),10β(CH3); V = 5β(CH3),10β(H); VII = 5α(CH3),10α(H). Percentages represent the relative abundance of the most abundant peak in each transition. (B) SIM m/z = 231 chromatogram of C26–28 triaromatic steroids of PE 52336 Fossil (i) and Matrix (ii). S/R refers to the stereochemistry at C-20.
Biology 11 01289 g003
Figure 4. Dietary ChemoSpace analysis of tissue type signals preserved in carbonaceous coprolites from the Mazon Creek locality. A total of n = 37 Mazon Creek carbonaceous fossils and coprolites were spectroscopically fingerprinted (n = 1 biological replicate and n = 10 mean-averaged technical replicates per data point), identified, and analysed by means of a Canonical Correspondence Analysis (CCA; discriminant analysis). (A) ChemoSpace resulting from a CCA that treated coprolites as a separate tissue category (black arrows); corresponding trajectories of functional groups in the ChemoSpace are plotted in Figure 5B. (B) ChemoSpace resulting from a CCA that treated the tissue affinity of coprolites as unknown; corresponding trajectories of functional groups in the ChemoSpace are plotted in Figure 5C. ChemoSpaces in both, (A) (clustered in n = 5 categories, each contained in a dotted, convex hull; except from annelids and coprolites due to small sample size) and (B) (clustered in n = 4 categories, each contained in a dotted, convex hull, except from annelids), reveal that fossil coprolites are compositionally distinct from undigested fossil soft tissues due to digestive hydrolysis of macromolecules prior to fossilisation (detailed differences are plotted in Figure 5B), and contain predominantly digested tissues of vertebrate prey items. Thus, carnivory (orange circle fraction, contrary to omnivory/herbivory = green circle fraction) can be inferred for both coprolite producers. Source Data are available for the CCAs in (A,B).
Figure 4. Dietary ChemoSpace analysis of tissue type signals preserved in carbonaceous coprolites from the Mazon Creek locality. A total of n = 37 Mazon Creek carbonaceous fossils and coprolites were spectroscopically fingerprinted (n = 1 biological replicate and n = 10 mean-averaged technical replicates per data point), identified, and analysed by means of a Canonical Correspondence Analysis (CCA; discriminant analysis). (A) ChemoSpace resulting from a CCA that treated coprolites as a separate tissue category (black arrows); corresponding trajectories of functional groups in the ChemoSpace are plotted in Figure 5B. (B) ChemoSpace resulting from a CCA that treated the tissue affinity of coprolites as unknown; corresponding trajectories of functional groups in the ChemoSpace are plotted in Figure 5C. ChemoSpaces in both, (A) (clustered in n = 5 categories, each contained in a dotted, convex hull; except from annelids and coprolites due to small sample size) and (B) (clustered in n = 4 categories, each contained in a dotted, convex hull, except from annelids), reveal that fossil coprolites are compositionally distinct from undigested fossil soft tissues due to digestive hydrolysis of macromolecules prior to fossilisation (detailed differences are plotted in Figure 5B), and contain predominantly digested tissues of vertebrate prey items. Thus, carnivory (orange circle fraction, contrary to omnivory/herbivory = green circle fraction) can be inferred for both coprolite producers. Source Data are available for the CCAs in (A,B).
Biology 11 01289 g004
Figure 5. In vivo digestive hydrolysis reactions of key macromolecules and functional group trajectories in the compositional spaces plotted in Figure 4. (A) Functional group chemistry during digestive hydrolysis reactions of proteins, polysaccharides, adenosine triphosphate (ATP), and lipids: stomach acid and digestive enzymes catalyse the hydrolytic cleavage of amide, glycosidic, and ester bonds and yield a relative increase in the relative abundance of alcohols and carbonyls (red), phosphate (teal), and aromatic compounds (dark blue); Depending on the dietary source, faecal matter can be enriched in organo-sulfur moieties (orange; S-bearing amino acids: cystein, methionine). Compositional differences associated with digestive hydrolysis survive (as shown in Figure 4) biomolecule fossilization through oxidative crosslinking. ADP = adenosine diphosphate; ATP = adenosine triphosphate. (B) Discriminant vector arrows for tissue categories (black, n = 5) and functional groups (n = 53 extracted relative intensities) in the Canonical Correspondence Analysis (CCA) shown in Figure 4A. Functional groups are colour-coded (corresponding to functional group labels in (A)) for vectors that explain the ChemoSpace placement of the two analyzed coprolites. This CCA treated coprolites (n = 2) as a separate tissue category to reveal the unique compositional features distinguishing them from other types of carbonaceous soft tissues from the Mazon Creek locality (n = 35). (C) Discriminant vector arrows for tissue categories (black, n = 4) and functional groups (n = 53 extracted relative intensities) in the CCA shown in Figure 4B. Functional groups are colour-coded (corresponding to functional group labels in (A)) for vectors that explain the ChemoSpace placement of the two analysed coprolites (see Figure 4B). This CCA treated coprolites (n = 2) as samples with unknown tissue affinity to identify the primary source of fossil faecal matter (n = 35).
Figure 5. In vivo digestive hydrolysis reactions of key macromolecules and functional group trajectories in the compositional spaces plotted in Figure 4. (A) Functional group chemistry during digestive hydrolysis reactions of proteins, polysaccharides, adenosine triphosphate (ATP), and lipids: stomach acid and digestive enzymes catalyse the hydrolytic cleavage of amide, glycosidic, and ester bonds and yield a relative increase in the relative abundance of alcohols and carbonyls (red), phosphate (teal), and aromatic compounds (dark blue); Depending on the dietary source, faecal matter can be enriched in organo-sulfur moieties (orange; S-bearing amino acids: cystein, methionine). Compositional differences associated with digestive hydrolysis survive (as shown in Figure 4) biomolecule fossilization through oxidative crosslinking. ADP = adenosine diphosphate; ATP = adenosine triphosphate. (B) Discriminant vector arrows for tissue categories (black, n = 5) and functional groups (n = 53 extracted relative intensities) in the Canonical Correspondence Analysis (CCA) shown in Figure 4A. Functional groups are colour-coded (corresponding to functional group labels in (A)) for vectors that explain the ChemoSpace placement of the two analyzed coprolites. This CCA treated coprolites (n = 2) as a separate tissue category to reveal the unique compositional features distinguishing them from other types of carbonaceous soft tissues from the Mazon Creek locality (n = 35). (C) Discriminant vector arrows for tissue categories (black, n = 4) and functional groups (n = 53 extracted relative intensities) in the CCA shown in Figure 4B. Functional groups are colour-coded (corresponding to functional group labels in (A)) for vectors that explain the ChemoSpace placement of the two analysed coprolites (see Figure 4B). This CCA treated coprolites (n = 2) as samples with unknown tissue affinity to identify the primary source of fossil faecal matter (n = 35).
Biology 11 01289 g005
Figure 6. General schematic demonstrating diagenetic rearrangements of cholesteroid compounds, as identified in fossil coprolites in the current study. Dashed lines represent multiple possible intermediates or pathways.In both fossil and matrix of all samples, n-alkanes are abundant and consistently range from n-C15 to n-C26 with a maximum at n-C23 but with no odd or even carbon number preference (Figure 7). The predominance of mid-chain (C20-C25) n-alkanes and lack of high-molecular-weight n-alkanes supports input from aquatic organic material in a limnic or deltaic environment receiving a minimal input of land plant material, e.g., [91,92]. The freshwater-influenced n-alkane distribution is thus in agreement with a geologically inferred large delta setting, e.g., [41,42,43,44,45]. The regular isoprenoids pristane (Pr) and phytane (Ph) are also present with Pr/Ph ratios ranging from 0.85–1.65 (Table 2), ratios supporting fluctuating redox conditions within the environment [73].
Figure 6. General schematic demonstrating diagenetic rearrangements of cholesteroid compounds, as identified in fossil coprolites in the current study. Dashed lines represent multiple possible intermediates or pathways.In both fossil and matrix of all samples, n-alkanes are abundant and consistently range from n-C15 to n-C26 with a maximum at n-C23 but with no odd or even carbon number preference (Figure 7). The predominance of mid-chain (C20-C25) n-alkanes and lack of high-molecular-weight n-alkanes supports input from aquatic organic material in a limnic or deltaic environment receiving a minimal input of land plant material, e.g., [91,92]. The freshwater-influenced n-alkane distribution is thus in agreement with a geologically inferred large delta setting, e.g., [41,42,43,44,45]. The regular isoprenoids pristane (Pr) and phytane (Ph) are also present with Pr/Ph ratios ranging from 0.85–1.65 (Table 2), ratios supporting fluctuating redox conditions within the environment [73].
Biology 11 01289 g006
Figure 7. GC-MS total ion chromatogram of PE 52336 showing differences in n-alkanes distribution and cholestanes abundances in the Fossil (A) versus the Matrix (B). Numbers represent the carbon number of the hydrocarbon chain. Pr = pristane, Ph = phytane.
Figure 7. GC-MS total ion chromatogram of PE 52336 showing differences in n-alkanes distribution and cholestanes abundances in the Fossil (A) versus the Matrix (B). Numbers represent the carbon number of the hydrocarbon chain. Pr = pristane, Ph = phytane.
Biology 11 01289 g007
Table 1. Steroid distributions in extracted organic matter.
Table 1. Steroid distributions in extracted organic matter.
PE52315PE52316PE52336
FossilMatrixFossilMatrixFossilMatrix
1 Steranes (%)C2786.269.192.850.199.078.0
C287.316.33.525.00.69.8
C296.514.63.724.90.412.2
2 Reg steranes/hopanes1.90.94.80.580.11.1
3 Monoaromatic steroids (%)C2754.841.868.338.982.256.1
C2831.727.725.125.915.222.7
C2913.530.56.635.12.521.2
4 Triaromatic steroids (%)C2743.626.784.641.494.868.4
C2819.427.26.119.92.711.5
C2937.046.19.338.72.520.1
1 Distributions of regular steranes were computed using the sum of regular steranes and diasteranes in MRM M+ → 217 transitions. Steranes: C27–29 ααα- and αββ-20(S + R)-steranes and C27 βαα 20R; ααα = 5α(H),14α(H),17α(H); αββ = 5α(H),14β(H),17β(H); βαα = 5β(H),14 α(H),17 α(H); diasteranes: βα-22(S + R)-diasteranes; 13β(H),17α(H). 2 Reg steranes/hopanes = [∑(C27–29 steranes)]/[∑(C27–35 hopanes)], steranes: as above; hopanes: C27 Ts, Tm 17α, Tm 17β; C29 αβ, Ts, βα; C30 αβ, βα; C31–35 αβ-22(S + R); αβ = 17α(H),21β(H); βα = 17β(H),21α(H); C27 Ts = 18α-22,29,30-trisnorneohopane; C27 Tm = 17α-22,29,30-trisnorhopane; C29 Ts = 18α-30-Norneohopane. 3 Distributions of monoaromatic steroids (MAS) were computed using the sum of S and R I + V isomers of MAS C27–29 homologs in m/z 253 mass chromatograms. 4 Triaromatic steroid distributions (TAS) were computed using the sum of C26–28 homologs in the m/z 231 mass chromatograms.
Table 2. Values of sterol compounds reported as total amount in micrograms (μg) total in each fraction, quantified against external sterol standards in 500 mL, as an average of two measurements.
Table 2. Values of sterol compounds reported as total amount in micrograms (μg) total in each fraction, quantified against external sterol standards in 500 mL, as an average of two measurements.
PE 52315PE 52316PE 52336Procedural Blank
FossilMatrixFossilMatrixFossilMatrix
5α-Cholestan-3β-ol--0.43-0.16--
Cholesterol0.070.060.050.170.070.140.07
Coprostan-3-ol0.03------
Table 3. Biomarker parameters in extracts of coprolite samples.
Table 3. Biomarker parameters in extracts of coprolite samples.
PE 52315PE 52316PE 52336
FossilMatrixFossilMatrixFossilMatrix
1 Dia/reg C27 steranes0.641.030.301.500.320.76
2 Dia/reg C28 steranes1.612.350.741.870.442.09
3 Dia/reg C29 steranes1.251.520.891.47*1.46
C29 20S/(20S + 20R) steranes0.280.300.180.270.210.25
4 Ts/(Ts + Tm)0.520.550.470.520.500.44
βα/(βα + αβ) C30 hopane0.110.100.080.090.090.11
5 Pr/Ph1.000.990.851.031.651.42
1 Dia/Reg C27 steranes = [∑(βα-22(S + R)-diacholestane)]/[∑(ααα- and αββ-20(S + R)-cholestane)]; 2 Dia/Reg C28 steranes = [∑(βα-22(S + R)-diaergostane)]/[∑(ααα- and αββ-20(S + R)-ergostane)]; 3 Dia/Reg C29 steranes = [∑(βα-22(S + R)-diastigmastane)]/[∑(ααα- and αββ-20(S + R)-stigmastane)]; 4 Ts = 18α-22,29,30-trisnorneohopane; Tm = 17α-22,29,30-trisnorhopane. Steranes and hopanes identified using the MRM M+ → 217 and M+ → 191 transitions, respectively. 5 Pristane and phytane integrated using total ion chromatogram. * Value omitted—coelution with abundant C27 steroid hydrocarbons was such that value was not able to be determined.
Table 4. Average δ13C values (bulk residue and compound specific) for PE 52316 and PE 52336, given in per mil (‰).
Table 4. Average δ13C values (bulk residue and compound specific) for PE 52316 and PE 52336, given in per mil (‰).
PE 52316PE 52336
FossilMatrixFossilMatrix
δ13Corg−23.6 (0.07)3−23.8 (0.03)2−23.9 (0.12)3−23.7 (0.03)3
δ13C17−28.4 (0.65 *)3−30.0 (0.24)3−29.4 (0.54 *)3−29.0 (0.62 *)3
δ13CPristane−28.8 (0.82 *)2−29.2 (0.30)3−28.4 (0.53 *)3−29.3 (0.51 *)3
δ13C18−30.4 (0.41 *)3−31.7 (0.31)3−31.5 (0.13)3−30.7 (0.40)3
δ13CPhytane−33.0 (0.11)2−29.5 (0.23)3−35.1 (0.39)3−33.1 (0.31)3
1δ13C20–25−30.5 −31.0 −30.8 −30.7
δ13Ccholestane−32.9 (0.39)3-−32.6 (0.29)3-
* Number in brackets indicates standard deviation; superscript refers to number of analyses used in average. Values with standard deviations greater than 0.4‰ are marked with an asterisk (*). 1δ13C values for C20–25 n-alkanes represent an average of the δ13C values for each n-alkane from C20 to C25, which were each determined from the average of three isotopic measurements. Standard deviations ranged from 0.02–0.60.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tripp, M.; Wiemann, J.; Brocks, J.; Mayer, P.; Schwark, L.; Grice, K. Fossil Biomarkers and Biosignatures Preserved in Coprolites Reveal Carnivorous Diets in the Carboniferous Mazon Creek Ecosystem. Biology 2022, 11, 1289. https://doi.org/10.3390/biology11091289

AMA Style

Tripp M, Wiemann J, Brocks J, Mayer P, Schwark L, Grice K. Fossil Biomarkers and Biosignatures Preserved in Coprolites Reveal Carnivorous Diets in the Carboniferous Mazon Creek Ecosystem. Biology. 2022; 11(9):1289. https://doi.org/10.3390/biology11091289

Chicago/Turabian Style

Tripp, Madison, Jasmina Wiemann, Jochen Brocks, Paul Mayer, Lorenz Schwark, and Kliti Grice. 2022. "Fossil Biomarkers and Biosignatures Preserved in Coprolites Reveal Carnivorous Diets in the Carboniferous Mazon Creek Ecosystem" Biology 11, no. 9: 1289. https://doi.org/10.3390/biology11091289

APA Style

Tripp, M., Wiemann, J., Brocks, J., Mayer, P., Schwark, L., & Grice, K. (2022). Fossil Biomarkers and Biosignatures Preserved in Coprolites Reveal Carnivorous Diets in the Carboniferous Mazon Creek Ecosystem. Biology, 11(9), 1289. https://doi.org/10.3390/biology11091289

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop