2,777
Views
27
CrossRef citations to date
0
Altmetric
Articles

The influence of varying proportions of terrestrial and marine dietary protein on the stable carbon-isotope compositions of pig tissues from a controlled feeding experiment

, , , , , , , & show all
Pages 28-44 | Received 03 Oct 2016, Accepted 16 Dec 2016, Published online: 16 Mar 2017

ABSTRACT

In recent years, it has become evident that limitations exist in our ability to meaningfully assess palaeodiet using stable isotope compositions. These limitations in part arise because many of the fundamental assumptions about tissue-diet relationships are poorly understood. In order to redress this deficiency, a controlled feeding experiment was undertaken to define the impact of terrestrial- vs. marine-derived dietary protein consumption on consumer tissue carbon isotopic compositions (δ13C). Two generations of pigs were raised on one of five feeds with varying proportions of terrestrial (soy) and marine (fish meal) protein. A comprehensive range of tissues and fluids from 49 pigs was submitted for δ13C analysis.

The observed tissue–whole diet and tissue–dietary protein carbon isotopic offsets were found to be highly dependent on the percentage of marine protein in diet. We suggest that the trend in δ13C offsets most likely derives from the increased routing of non-essential amino acids, especially glycine, with the increasing proportion of marine protein in the diet. These findings demonstrate that solely using bulk δ13C compositions not only masks considerable information about diet, but may also lead to erroneous representations of marine and terrestrial resource consumption in the past.

Introduction

Of central importance in many archaeological investigations are patterns of resource use through time and across geographic regions, particularly during environmentally- or culturally-mediated periods of dietary change, such as the postulated shift away from marine resource consumption during the Mesolithic-Neolithic transition across Europe (Richards et al., Citation2003; Tauber et al., Citation1981). Stable carbon-isotope (δ13C) analysis of human and faunal remains has been used for over three decades to investigate palaeodiet. These determinations make valuable contributions to our understanding of resource exploitation in the past where differential access and exploitation can have important implications for socio-political and economic behaviour. In ecologically-complex areas where numerous resources are available for exploitation, stable isotopic compositions are similarly used to evaluate relative consumption of different kinds of food, for example, Mesoamerica (e.g. Emery et al., Citation2000; Scherer et al., Citation2007; Warinner et al., Citation2013; White, Citation2005; White et al., Citation1993) or the western coast of South America (e.g. Finucane et al., Citation2006; Gil et al., Citation2011; Knudson et al., Citation2015; Tomczak, Citation2003; Webb et al., Citation2013; Yesner et al., Citation2003).

Early stable isotope research (DeNiro & Epstein, Citation1976, Citation1978) determined that there was a strong correlation between the stable isotopic compositions of various body tissues and diet. It was hypothesised that bone mineral δ13C values (bioapatite; δ13Csc) should reflect the δ13C values of dietary carbohydrates, lipids, and, to a lesser extent, protein, whereas bone protein δ13C values (collagen; δ13Ccol) should most strongly reflect dietary protein (Krueger & Sullivan, Citation1984). Numerous subsequent feeding experiments explored the influence of varying the δ13C values of different dietary macronutrients, as well as the impact of different relative proportions of dietary protein, carbohydrates and lipids on the relationships among δ13Csc, δ13Ccol, δ13Cwhole diet and δ13Cdietary protein (e.g. Ambrose & Norr, Citation1993; Howland et al., Citation2003; Jim et al., Citation2006; Krueger & Sullivan, Citation1984; Passey et al., Citation2005; Tieszen & Fagre, Citation1993). These studies further investigated the process of carbon routing from dietary macronutrients, and proposed that a reasonably straightforward linear mixing model best explained the relationship between δ13Csc and δ13Cwhole diet values. Ultimately, they also confirmed that collagen δ13C values reflect the δ13C of dietary protein when diet is protein-sufficient (Ambrose & Norr, Citation1993; Chisholm et al., 1982; Tieszan & Fagre, Citation1993).

In recent years, however, it has become increasingly evident that there are limitations associated with our ability to meaningfully assess palaeodiet using stable carbon-isotope compositions. In part this is because of the complexity of questions now being investigated, but largely because many of the fundamental assumptions about tissue-diet relationships are poorly understood (Hedges, Citation2004). Beyond constraints on interpretation imposed by the archaeological record, attempts to refine the predictive ability of isotopic compositions have included the applications of an index to estimate degree of meat consumption (Δ13Capcol), using multiple light isotope proxies, including nitrogen and sulphur (e.g. Druker & Bocherens, Citation2004; Hedges & Reynard, Citation2007; Müldner & Richards, Citation2005; Nehlich et al., Citation2010, Citation2011; Reynard & Tuross, Citation2015; Richards et al., Citation2003; Sayle et al., Citation2013). New approaches using stable strontium (δ88Sr), calcium (δ44Ca), iron (δ56Fe), copper (δ65Cu) and zinc (δ66Zn) isotopes from mineralized tissues have also been applied (e.g. Jaouen et al., Citation2013; Knudson et al., Citation2010; Reynard et al., 2007). Meta-analyses of published isotopic datasets, incorporating isotopic data from archaeological humans and fauna, as well as controlled feeding studies, have enabled reassessment of many of the earlier dietary models and the outcomes of feeding experiments. Drawing on these data, increasingly sophisticated mathematical models have been developed and applied to the problem of archaeological human palaeodietary reconstruction (e.g. Fernandes et al., Citation2012, Citation2014; Froehle et al., Citation2012; Kellner & Schoeninger, Citation2007; Newsome et al., Citation2004). Models, although useful, are unfortunately limited by a lack of knowledge about dietary routing, tissue-diet relationships, and fractional contributions of different macronutrients (e.g. dietary amino acids, fatty acids, etc.) to different tissues. Moreover, many of these models have been adapted from ecological research in which ecosystems and ecological niches are well-constrained in terms of natural isotopic variability and dietary resource exploitation. However, models appropriate for use in these contexts are unlikely to be uncritically applicable to human palaeodietary studies. Crucially, human food choice is more than a function of the relative quantities of various edible resources in the environment; indeed, edibility is understood to be culturally-defined, as are food preparation and consumption methods (e.g. cooking or assembling meals), which can impact the quality of dietary macronutrients and how they are absorbed. Similarly, dietary heterogeneity resulting from socioeconomic factors and differential food access (e.g. status, economic role, trade and exchange, age, or gender) further complicates palaeodietary reconstruction. This suite of potential environmental, social and economic influences on human dietary choice creates a complex context in which to reconstruct palaeodiet. A thorough understanding of the biological and metabolic factors influencing the relationship between consumed food and the biochemical/metabolic record archived in archaeological and faunal remains will thus add much-needed rigour to palaeodietary reconstruction, and mitigate many of the uncertainties associated with elucidating resource exploitation in the past.

With these considerations in mind, a controlled feeding study was undertaken led by the University of Bristol. Two successive generations of pigs were fed one of five feeds containing varying proportions of terrestrial (soy) and marine (fish meal) protein. Pigs were chosen since they are widely recognised as an excellent physiological analogue for humans although there are differences in gestation time, average number of offspring per pregnancy and time to reach physical and sexual maturity (inter alia Heinritz et al., Citation2013; Litten-Brown et al., Citation2010; Sullivan et al., Citation2001; Swindle et al., Citation2012). The overarching goal of this study is to improve the understanding of consumer tissue isotopic composition for palaeodietary reconstruction, particularly the exploitation of marine resources. Here, our objectives are to elucidate tissue-tissue and tissue-diet carbon isotopic discrimination under different conditions of dietary protein consumption ranging from 100% terrestrial protein to 100% marine-derived protein. A comprehensive range of tissues and fluids including bone collagen (femoral and rib), muscle (femoral and loin), liver, blood, plasma, milk, hair, and faeces from 49 pigs from the first and second generation of the study were analysed for their δ13C compositions and a preliminary set of collagen and muscle tissues have been analysed for their individual amino acid δ13C (δ13CAA) values. The δ13C compositions of pig feeds, archived throughout the study, provide a defined dietary carbon isotopic baseline.

Theoretical considerations for palaeodietary reconstruction

The reconstruction of palaeodiet using stable isotope analysis is based on the well-tested assumption that tissue isotopic composition reflects the isotopic composition of consumed food (Ambrose, Citation1993). There are systematic differences in isotopic composition between tissue and food, and among different tissues within the body (Ambrose, Citation1993; DeNiro & Epstein, Citation1978). Here we define differences between tissues and diet and different tissues using a capital delta notations where Δ13CX−Y = δ13Cx − δ13Cy. These tissue-diet and tissue-tissue offsets are a result of isotopic fractionation, the differential partitioning of isotopes between phases in a reaction (e.g. ingested food → consumer tissue) caused by the mass differences among isotopes of the same element (Urey, Citation1947). Isotopic data are used to assess relative contributions to diet of isotopically distinct foods by comparing tissue isotopic compositions to a food web model that describes the variability in isotopic compositions of local food resources (i.e. terrestrial and marine fauna and plant material; Ambrose, Citation1993; Kellner & Schoeninger, Citation2007).

The δ13C composition of any proteinaceous tissue is a weighted average of the δ13C compositions of both its essential and non-essential amino acids (δ13CAA). Essential amino acids (e.g. threonine, valine, methionine, isoleucine, leucine, histidine, lysine, and phenylalanine) cannot be generated by the body and therefore must be ingested in sufficient quantities. Non-essential amino acids (e.g. asparagine/aspartic acid, hydroxyproline, glutamic acid/glutamate, serine, glycine, alanine, proline and arginine) can be assimilated with minimal modification from a dietary source, or may be synthesised de novo using components drawn from the body’s biochemical pools (Ambrose, Citation1993; Schwarcz, Citation2000). As a result, tissue essential amino acid δ13C values are expected to closely approximate dietary essential amino acid δ13C values due to direct routing, that is, δ13Ctissue AAδ13Cdiet AA ≈ 0‰. Non-essential amino acid δ13C values, however, may show evidence of both direct routing and biosynthesis.

The controls of the relative proportion of direct routing vs. de novo synthesis for non-essential amino acids is dependent on both the total protein in the diet and the amount of a given amino acid in the diet. In a high protein diet direct routing is more likely to occur, since it is more energetically efficient to assimilate non-essential amino acids than it is to synthesise them and one would expect δ13Ctissue AA to approximate δ13Cdiet AA (Corr et al., Citation2005; Jim et al., Citation2006; Jones, Citation2002; Schwarcz, Citation2000; Umbarger, Citation1978). Conversely when the diet is protein poor de novo synthesis dominates and a stronger relationship to whole diet δ13C values rather than to the corresponding dietary amino acid δ13C value is expected. A final caveat is that whilst it is convenient to discuss total dietary protein the metabolic effects are driven at the amino acid level. Thus for a diet that may be poor in overall protein but rich in a particular individual non-essential amino acid, say glycine, de novo synthesis may dominate for most non-essential amino acids but glycine will be directly routed.

Although this variability in essential and non-essential amino acid metabolism and its influence on amino acid isotopic compositions is subsumed within the sampling resolution for this study, these processes underlie “bulk”Footnote1 tissue δ13C compositions, and their impact may be detectable at the bulk protein sampling resolution under controlled feeding conditions.

Based on numerous controlled feeding and well-constrained wild faunal studies, it is generally assumed that bone collagen δ13C values are elevated by ∼+5‰ relative to whole diet, and by 0 to +1‰ relative to dietary protein (Ambrose & Norr, Citation1993). The δ13C values of other proteinaceous tissues, including muscle, blood and plasma are ∼1‰ higher than whole diet, whereas liver and faeces δ13C values are generally moderately reduced relative to diet (Sponheimer et al., Citation2003, Citation2006). Many experimental studies show considerable variability in the relationships between collagen and whole diet and collagen and dietary protein. Larger tissue–diet offsets are typically attributed to experimentally-controlled isotopic differences between the “energy” (non-protein) and protein components of diet, and to variable contributions to bone collagen from dietary energy vs. dietary protein (explored in Froehle et al., Citation2010, Citation2012; Kellner & Schoeninger, Citation2007). Early work determined that, by varying the energy and protein component isotopic compositions, it was possible to induce tissue-diet offsets ranging from ∼−2 to +10‰ (e.g. Ambrose & Norr, Citation1993; Young, Citation2002). Moreover, although the δ13Ccol values do predominately reflect dietary protein δ13C values, direct routing of protein from diet to tissue can occur, indicating that non-essential amino acids can be assimilated directly rather than biosynthesised from non-protein carbon skeletons. At higher levels of protein consumption, dietary amino acids are routed more or less directly from diet to tissue with minimal modification. When protein consumption is low, however, carbon drawn from dietary carbohydrates and lipids can constitute ∼49–58% of the carbon in collagen. If that small amount of consumed protein is from a marine resource, with the remainder of the diet derived from plant foods, tissue δ13C values may be quite low, obscuring the marine carbon contribution. In contrast, protein-rich diets that incorporate both marine and terrestrial protein resources may over-estimate marine protein consumption due to preferential routing of dietary marine protein to collagen. Preliminary models suggest that a contribution of up to 20% marine protein to diet may cause tissue isotopic changes of only +0.3 to +1.8‰, depending on the overall protein content of diet and thus on the importance of routing versus biosynthesis of non-essential amino acids (Hedges, Citation2004).

Methodology

Controlled feeding study

The pigs were raised at Harper Adams University (Shropshire, UK), and a full suite of tissues, fluids and excreta was archived (). Each pig was fed one of five diets of known dietary protein source composition: (i) 100% terrestrial-derived (soy), (ii) 87.5% terrestrial/12.5% marine, (iii) 75% terrestrial/25% marine, (iv) 50% terrestrial/50% marine, and (v) 100% marine-derived (fish meal). All five diets were nutritionally equivalent, that is, there was no difference in the amount or quality of dietary protein, carbohydrates or lipids among feeds. The experiment was run over two successive generations, in the first generation, five groups of gilts were fed one of the above diets from weaning until sacrifice. All gilts were artificially inseminated and the second generation pigs were fed exclusively on one of the five diets from weaning until sacrifice in adolescence. Five tonnes of each feed was produced from the same batches of ingredients by Parnutt Foods, Ltd. (Sleaford, UK) and held in cold storage until required. Diet formulations are presented in . Four subsamples of each feed were taken as batches were released from storage to confirm homogeneity.

Figure 1. Flowchart illustrating the structure of the controlled feeding study.

Figure 1. Flowchart illustrating the structure of the controlled feeding study.

Table 1. Diet formulations and carbon-isotope compositions.

This study was explicitly designed to address limitations recognized in earlier feeding studies, particularly difficulties associated with sample size, differential tissue turnover, and nutritional stress. As such, each diet group is represented by several pigs, all of which have only consumed a single experimental diet (including sow milk from the same dietary group), and only one dietary variable – the ratio of terrestrial to marine protein – was changed. All five diets have a constant 20% protein contribution to whole diet, which eliminates variability in isotopic data and tissue-diet offsets associated with low (≤5%) or high (∼70%) protein consumption. In total, 10 sows (first generation), 19 piglets (second generation, aged four weeks) and 39 pigs (second generation, aged >160 days) were reared and sacrificed over the course of the study.

Laboratory methods

Ten first generation sows, 19 piglets, and 20 second generation pigs were selected for isotopic analysis. Bone collagen (femur and rib), muscle (femoral and loin) and liver samples were analyzed for their δ13C compositions. Milk (n = 10; sows only) and hair (n = 10; piglets only) samples were also analyzed. Plasma (sampled at 70 and 140 days of age), blood, faeces and urine samples were analyzed as available. Bone collagen was extracted using a modified Longin (Citation1971) method. For each femur or rib sample, a section of bone was taken using a hacksaw or rotary tool (Dremel tool, 3000JB with diamond cutting wheel SC545). The bone was mechanically defleshed using a scalpel, freeze-dried, and a rotary tool with a silicon carbide burr was used to remove all trabecular bone and ∼0.5 mm of surface bone. Each sample was then crushed and sieved, retaining approximately 450 mg of ≥212 µm and ≤1 mm fragments for collagen extraction. Lipids were extracted using 2:1 v/v chloroform: methanol solution (3 × 8 mL solvent solution, 3 × 20 min sonication). Collagen was extracted by soaking in 0.5 M hydrochloric acid (HCl) until bone chips were entirely soft. Extracted collagen was solubilised in 10−3 M HCl at 75°C for 48 h, filtered (E-Zee filters, 60–90 µm), and freeze-dried for ≥24 h.

For soft tissue samples, a cross section of the tissue was removed using a new scalpel blade. Soft tissues, faeces and milk were then freeze dried for more than 48 hours and then lipid extracted using a 2:1 v/v chloroform: methanol solution (3 × 8 mL solvent solution, 3 × 20 min sonication). Fluid samples were freeze dried and hair samples were cleaned of contaminant lipids and other organics by soaking in 2:1 v/v chloroform: methanol solution (2 × 8 mL for >24 h). For all samples, homogenised aliquots were then weighed into tin capsules (∼0.70 ± 0.1 mg) for isotopic analysis. Twenty homogenised feed samples were weighed into tin capsules (∼1.4 ± 0.1 mg) for isotopic analysis. Prior to use, all glassware was washed with Decon 90 and solvent-rinsed before furnacing at 450°C for four hours. Aluminium foil and disposable gloves were used to handle samples.

All bulk isotopic analyses were performed using a Flash HT elemental analyser interfaced with a Thermo Electron DeltaPlus XP mass spectrometer at the Natural Environment Research Council Life Science Mass Spectrometry Facility in East Kilbride, Scotland. For collagen, methodological reproducibility was determined through duplicate collagen preparation and analysis for 10% of samples and was ± 0.1‰. Analytical reproducibility for all samples was assessed by repeated analyses of 10% of samples, and was also ±0.1‰. Analytical precision was monitored throughout using USGS-40, and the average δ13C value obtained over all analytical sessions was −26.3 ± 0.1‰, which compares well with the accepted value of −26.4‰ (Brand et al., Citation2014).

δ13C determinations on individual amino acids from the 5 diet feeds, 10 sow collagen and 10 sow muscle samples were made by LC-IRMS using a Thermo-Finnigan Surveyor HPLC coupled to a Delta V IRMS via an LC-Isolink (Krummen et al., Citation2004). Amino acids are separated using a PrimeStep A column (250 mm × 3.2 mm, 5 μm particle size, 100 Å pore size; SIELC Technologies Ltd., Prospect Heights, IL, USA) following the method of Dunn et al. (2011). Repeated analysis of an amino acid mixture gives reproducibility within ±0.5‰ of the EA-IRMS value.

Results

Growth performance

All results are presented as average ± 1σ [range] unless otherwise noted and all growth performance data are summarised in . There were no significant health issues throughout the production stages, and there was no evidence of nutritional stress or delayed growth among diet groups, which is reflected in the average weight gain per day and total weight gain.

Table 2. Pig growth performance summary.

Sows were sacrificed on 11th February 2011. The average sow weight during pregnancy (adjusted downward to compensate for total litter weight) was 195 ± 15 kg, and was 203 ± 15 kg when the piglets were weaned. The average weight at slaughter was 205 ± 13 kg. All piglets were born between 29th August and 5th September 2010. A subset of piglets (three to five per diet group) were sacrificed at approximately four weeks of age on 30th September 2010. The average birth weight for all piglets was 1.4 ± 0.2 kg, and did not differ significantly among sow diet groups (Kruskal-Wallis, p = 0.371). Similarly, there was no difference in weaning weight (6.4 ± 1.5 kg), total gain (4.9 ± 1.4 kg) or average gain per day (183 ± 55 g day−1) among piglets from different diet groups (Kruskal-Wallis, p > 0.36 for all tests).

The 39 second generation pigs were sacrificed on either 8th February or 9th March 2011 at ∼160 or ∼190 days of age, respectively. The average birth weight for all pigs was 1.5 ± 0.3 kg [0.7 to 2.2 kg], and differed significantly across the five diets (Kruskal-Wallis, p = 0.017), but the difference in weight among diet groups was no longer apparent at weaning. The average weaning weight was 14.6 ± 2.9 kg [7.3 to 20.0 kg], and did not vary significantly with maternal diet (Kruskal-Wallis, p = 0.261). Similarly, the total weight gain from weaning until slaughter (107.8 ± 13.9 kg [59.0 to 138.2 kg]) did not differ among the five diet groups (ANOVA, p = 0.466). The average gain per day was 624 ± 90 g day−1 [316 to 752 g day−1], and also did not differ among the five diet groups (ANOVA, p = 0.274). Finally, the average slaughter weight was 109.3 ± 13.8 kg [60.6 to 139.6 kg], and did not differ among pigs consuming different diets (ANOVA, p = 0.480).

Pig feed isotopic compositions

Lots 1 through 4 of each pig feed are isotopically equivalent for all five diets (Kruskal-Wallis, p = 0.288; Table S1). The differences in δ13C compositions of diet 1 through 5 are likewise not statistically significant (Kruskal-Wallis, p = 0.275). The average δ13C values over the five diets are: diet 1: −25.0 ± 0.2‰ [0.3‰], diet 2: −24.7 ± 0.2‰ [0.4‰], diet 3: −24.8 ± 0.2‰ [0.5‰], diet 4: −24.8 ± 0.2‰ [0.4‰], and diet 5: −24.5 ± 0.3‰ [0.5‰]. Thus, although the dietary protein source was qualitatively different for each diet, all pigs were consuming feed with the same overall δ13C composition (−24.7 ± 0.2‰). The δ13C value of the fish meal and soymeal dietary protein sources were determined to be −20.5‰ and −25.6‰, respectively.

Tissue isotopic compositions

Collagen

Collagen δ13C values were determined for both femoral and rib bone collagen for 10 sows, 19 pigletsFootnote2 and 20 pigs (Tables S2-S4). The average rib collagen δ13C values are −20.8 ± 1.5‰ [4.5‰] for sows, −20.0 ± 1.4‰ [5.0‰] for piglets and −20.2 ± 1.5‰ [4.7‰] for pigs. The average femoral bone collagen δ13C values are −20.2 ± 1.4‰ [4.1‰] for sows, −19.6 ± 1.2‰ [4.2‰] for piglets and −22.6 ± 1.4‰ [4.1‰] for pigs. Within each age category, rib and femoral bone collagen are not statistically significantly different (Mann-Whitney, sows p = 0.247; piglets p = 0.327; pigs p = 0.512). Similarly, across age categories, neither rib nor femoral bone collagen δ13C values are significantly different (Kruskal-Wallis, p = 0.227 for femoral collagen and p = 0.274 for rib collagen). As such, an average collagen δ13C value was determined for each pig, and all age categories are pooled for the remainder of this paper (). We caution however, that although there are no interpretively or statistically important differences among sows, piglets and pigs at the bulk carbon-isotope level for collagen, this determination may not hold true at the compound-specific amino acid level of analysis.

Figure 2. The δ13C values of all tissues for all pigs discussed herein. For feed samples, the error bars represent one standard deviation about the average value. Note that collagen and muscle carbon-isotope compositions are average values for rib and femoral collagen and loin and femoral muscle, respectively.

Figure 2. The δ13C values of all tissues for all pigs discussed herein. For feed samples, the error bars represent one standard deviation about the average value. Note that collagen and muscle carbon-isotope compositions are average values for rib and femoral collagen and loin and femoral muscle, respectively.

Muscle and liver

Muscle δ13C values were determined for both femoral and loin muscle samples for 10 sows, 19 pigletsFootnote3 and 20 pigs. The average femoral muscle δ13C values are −22.3 ± 1.4‰ [3.9‰] for sows, −22.3 ± 1.2‰ [3.6‰] for piglets and −22.6 ± 1.4‰ [4.1‰] for pigs. The average loin muscle δ13C values are −22.3 ± 1.4‰ [3.9‰] for sows, −22.1 ± 1.0‰ [3.3‰] for piglets and −22.5 ± 1.4‰ [4.2‰] for pigs. Liver δ13C values were determined for 10 sows, 19 piglets and 20 pigs, and are −23.2 ± 1.1‰ [3.4‰], −22.9 ± 1.0‰ [3.3‰] and −23.7 ± 1.2‰ [3.9‰], respectively. Within each age category, femoral and loin muscle are not statistically significantly different (Mann-Whitney, sows p = 0.912; piglets p = 0.367; pigs p = 0.620). Across age categories, femoral and loin muscle δ13C values are similarly statistically indistinguishable (Kruskal-Wallis, p = 0.305 for femoral muscle and p = 0.336 for loin muscle). Thus, for all pigs with both femoral and loin muscle isotopic data, an average δ13C value is used for the remainder of this paper. Liver δ13C values are statistically different between piglets and pigs (Mann-Whitney, p = 0.010), but the absolute difference in δ13C is only 0.8‰ and this is not considered interpretively significant (). Other inter-age category comparisons do not reveal statistically significant differences.

Blood, plasma, milk, hair and faeces

For the second generation pigs only, blood and plasma samples were collected and analysed for their δ13C composition. The average blood δ13C value was −23.7 ± 1.4‰ [4.0‰]. The average plasma δ13C value at 70 days of age was −23.0 ± 1.2‰ [3.2‰], and −23.3 ± 1.3‰ [3.8‰] at 140 days of age; these paired samples do not have statistically different δ13C values (Mann-Whitney, p = 0.567). The average sow milk δ13C value was −22.5 ± 1.0‰ [2.9‰]. Piglet hair had an average δ13C value of −20.4 ± 1.2‰ [3.2‰]. Faeces (n = 15) samples were analysed as available from second generation pigs, and are generally lower in δ13C relative to tissues. The average δ13C value was −25.1 ± 0.7‰ [2.3‰] and is very similar to that of the feed ().

Differences across diet groups for all tissues and samples

Statistical comparison of all eight tissues and samples (collagen, muscle, liver, milk, hair, faeces, blood and plasma) has determined that there are differences in δ13C values among tissues (Kruskal-Wallis, p = 0.000; ). Collagen and hair δ13C values do not differ significantly, but both tissues are significantly different from the other six tissues and samples. Faeces are similarly distinct from all other tissues and samples. Liver, muscle, milk, blood and plasma are all statistically similar.

The Δ13Ctissue − tissue, Δ13Ctissue − whole diet, and Δ13Ctissue −dietary protein isotopic offsets

Tissue–tissue, tissue–whole diet and tissue–dietary protein isotopic offsets for all pigs and tissues are presented in , and , respectively. With the exception of Δ13Ctissue − faeces, there are no consistent, significant linear trends in tissue–tissue isotopic offsets from diets 1 to 5. There are some weak linear trends for the tissue-tissue isotopic offsets, for example, in the average Δ13Cmuscle − milk isotopic offsets which range from −0.3‰ (Diet 2) to +1.0‰ (Diet 5). However, these trends only narrowly exceed the analytical reproducibility associated with the isotopic measurements and the magnitude of the fully propagated error associated with replicated pigs consuming the same diet (± 0.5‰ across all diet groups). As faeces δ13C values closely reflect whole diet δ13C values, there is a strong linear relationship between pig tissues and faeces, wherein the Δ13Ctissue − faeces isotopic offset increases with greater contributions of marine protein to total dietary protein.

Figure 3. The δ13Ctissueδ13Cwhole diet isotopic offsets for all tissues for all pigs discussed herein. The error bars represent one standard deviation about the average value.

Figure 3. The δ13Ctissue – δ13Cwhole diet isotopic offsets for all tissues for all pigs discussed herein. The error bars represent one standard deviation about the average value.

Figure 4. The δ13Ctissueδ13Cdietary protein isotopic offsets for all tissues for all pigs discussed herein. The error bars represent one standard deviation about the average value.

Figure 4. The δ13Ctissue – δ13Cdietary protein isotopic offsets for all tissues for all pigs discussed herein. The error bars represent one standard deviation about the average value.

Table 3. Tissue–tissue isotopic offsets.

As established, muscle, plasma, blood, liver and milk are not statistically different, nor are collagen and hair. In contrast, the tissue–whole diet isotopic offsets vary significantly for all tissues across diet groups for all three age categories, ranging, for example, from +3.6 ± 0.1‰ to +7.0 ± 0.4‰ for pig collagen from diet 1 to 5 (). These offsets are also significantly linearly correlated with the percentage of marine protein in diet. Spearman’s correlation tests determined that p = 0.000 for all comparisons with ρ values ranging from 0.832 to 0.991 (fully reported in Table S5). The magnitude of change in Δ13Cfaeces – whole diet offset is smaller than that observed for other tissues and samples, but there is nonetheless a linear trend apparent in the data.

The δ13C values for dietary protein for diets 1 through 5 were estimated using balance equation. The δ13Cdietary protein values are thereby estimated to be −25.6, −25.0, −24.3, −23.1 and −20.5‰ for diets 1 through 5. The linear relationships between % marine protein and Δ13Ctissue – dietary protein isotopic offsets for collagen, muscle, liver, milk, hair and faeces are statistically significant, however the δ13Ctissueδ13Cdietary protein values do vary among the tissues ().

Preliminary δ13AA results

Preliminary δ13C­ results for individual amino acids from the 5 diet feeds, 10 sow collagen and 10 sow muscle samples are presented in and supplementary table 6. δ13CAA values for individuals on replicated diets are typically within 0.75 ‰ of each other and so are combined as an average δ13CAA value for each feeding group. For essential amino acids, Δ13CAA tissue – AA whole diet ≈ 0‰ as expected, since these amino acids must be directly routed from dietary protein. Non-essential amino acids in low marine protein feed groups show an increased Δ13CAA tissue – AA whole diet reflecting their mixed direct assimilation vs. de novo synthesis origins. As the marine protein content of the diet increases Δ13CAA tissue – AA whole diet for non-essential amino acids trends towards 0‰. Across all feeds and tissues there is a trend towards increasing δ13CAA with increasing marine protein content which is indicative of the marine origin of these proteins.

Figure 5. Δ13Ctissue AA – whole diet AA for essential and non-essential amino acids from sow bone collagen and muscle samples. For essential amino acids phenylalanine (◊), lysine (□), threonine (▵), valine (○) and leucine/isoleucine (▹) are shown for non-essential amino acids aspartic acid (◊), glutamic acid (□), glycine (▵), alanine (○) and proline (▹) are shown.

Figure 5. Δ13Ctissue AA – whole diet AA for essential and non-essential amino acids from sow bone collagen and muscle samples. For essential amino acids phenylalanine (◊), lysine (□), threonine (▵), valine (○) and leucine/isoleucine (▹) are shown for non-essential amino acids aspartic acid (◊), glutamic acid (□), glycine (▵), alanine (○) and proline (▹) are shown.

Discussion

Differences among individuals and tissues

The similarity in δ13C compositions across age categories for commonly-sampled tissues (i.e. collagen, muscle and liver) is not unexpected. All pigs consumed one of the controlled feeding study diets during periods of significant growth or for a prolonged period of time (i.e. sows and pigs). Piglets, although slaughtered at a young age and likely still consuming sow milk, would also have been consuming pelleted feed. Sow milk production typically peaks 21 days post-partum and the volume of milk produced after that time is insufficient to meet the nutritional requirements of growing piglets (Patience et al., Citation1995). Thus, it is expected that piglets were consuming pelleted feed ad libitum for at least one week before slaughter, which, coupled with a rapid growth rate and high tissue turnover, would obscure any potential milk-related 13C-depletion in piglet tissues relative to sow tissues.

Among pigs from all age categories consuming the same diet, tissue δ13C values were very similar, with standard deviations of ±0.1 to ±0.7‰ (). These low inter-individual differences suggest that, among healthy organisms consuming an identical, nutritionally-adequate diet, individual metabolic differences do not introduce significant variability into bone collagen, muscle, liver, blood, faeces, hair, milk or plasma δ13C compositions. This outcome further suggests that small differences in δ13C values among individuals may be interpretively significant. Although variable bone turnover rates may induce differences in tissue isotopic compositions intra-skeletally if diet has changed over several years (Cox & Sealy, Citation1997; Hedges et al., Citation2007; Hill & Orth, Citation1998; Taylor et al., Citation2013), our comparison of the femoral and rib collagen data support the contention that there is no inherent δ13C difference of the protein component of different bones given a constant dietary input. Similarly, co-forming femoral and loin muscle do not differ isotopically, further supporting the hypothesis that intra-tissue differences across an organism are minimal. As a result of the comparatively fast turnover of blood, plasma and liver (i.e. several days), these samples will reflect short-term diet. Diet did not change throughout the experiment and, as expected, there is good agreement among the δ13C values of these tissues, and between replicate plasma samples taken mid-study and shortly before slaughter. Similarly, faeces δ13C also reflect short-term diet (i.e. a few days). Faeces δ13C values were remarkably consistent among pigs consuming the same diet, ranging from ±0.3 to ±0.6 ‰ across the five diet groups, and there is also good agreement between faeces and whole diet δ13C values (+0.2 ± 0.3‰). This trend is readily observable in our data because diet did not change during the pig’s lifetime, but a similar relationship between whole diet and faeces has been suggested for other species in wild contexts (Codron et al., Citation2007; Sponheimer et al., Citation2003).

As expected, different tissue types have different δ13C compositions. These differences most likely result from the impact of tissue-specific rates tissue growth and of the rate of reaction on carbon-isotope fractionation, wherein the δ13C values of fast-growing tissues would thus be expected to be more similar to those of whole diet. Further, bulk tissue isotope compositions are essentially weighted averages of the δ13C values of the individual amino acids that constitute the protein(s) making up the tissue. Individual amino acids can vary significantly in their δ13C values; for example, isotopic ranges of greater than 20‰ have been determined under controlled feeding conditions for pig collagen amino acids (Hare et al., Citation1991; Howland, Citation2003; Jones, Citation2002). Thus, the δ13C composition of the various tissues and fluids sampled here are also expected to be moderately different because the dietary δ13CAA values vary among the five different feeds due to the changing dietary protein source. Tissues are relatively 13C-enriched as follows: δ13Ccollagen ≈ δ13Chair > δ13Cmuscle ≈ δ13Cmilk ≈ δ13Cblood ≈ δ13Cplasma > δ13Cliver > δ13Cfaeces ≈ δ13Cdiet.

Relationships between tissues, whole diet, and dietary protein carbon-isotope compositions

The δ13C offsets between each tissue and the whole diet δ13C composition of the corresponding pig feed were examined. As expected, the δ13C values of all tissues and faeces increased as the proportion of marine protein in diet increased. The overall change is small from 0 to 100% marine protein contribution (e.g. ∼+3.5‰ for collagen), but, again, this is not unexpected given the relatively small difference between the terrestrial (soy meal, −25.6‰) and marine (fish meal, −20.5‰) end members for this study. Contrary to the widely-held assumption that tissue – whole diet fractionation would be similar for all diets, we determined that the Δ13Ctissue – whole diet values vary widely relative to the isotopic compositions of tissues and whole diet. For example, the overall change in δ13Ccollagen values for pigs was +4.7‰ from diet 1 to diet 5 (−21.8 to −17.1‰, Supplementary Table 4), and the tissue – whole diet offset ranged from +2.8 to +7.2‰ (Supplementary Table 5), however, whole diet δ13C remains essentially unchanged across all diets. The Δ13Ctissue – whole diet values thus appear to be highly sensitive to marine protein consumption. In the context of archaeological and ecological research, this degree of variability is large enough to make dietary reconstruction considerably less informative. Without making a priori assumptions about the likely proportions of marine vs. terrestrial protein intake, relating the δ13Ctissue values to the isotopic composition of various dietary resources in the local food web becomes problematic.

It is possible that the variability in tissue–whole diet offset and its apparent relationship to the consumption of marine protein observed here is, in fact, spurious. Re-assessment of controlled feeding study δ13C data has clearly demonstrated that when diet is not monoisotopic, that is to say, the energy and protein components have significantly different δ13C compositions, the Δ13Ccollagen – whole diet offset will depend on the specific isotopic composition of the different dietary constituents (Ambrose & Norr, Citation1993; Froehle et al., Citation2012). The offset between dietary protein and whole diet (Δ13Cdietary protein – whole diet) appears to directly influence the Δ13Ccollagen – whole diet relationship, regardless of the kind of protein consumed. Thus, when the δ13Cdietary protein value is high relative to the δ13Cwhole diet value, the Δ13Ccollagen – whole diet offset is also high and the Δ13Ccollagen – dietary protein offset is low. As a result, the isotopic data will seemingly indicate that C4 plant/marine protein in diet contributes more carbon to collagen than does C3 plant foods (Froehle et al., Citation2010). Although this assertion cannot be challenged at the bulk δ13C level, it nonetheless does not take into consideration differences at the amino acid level. By pooling C4 plant and marine protein consumers together, this explanation implicitly treats “protein” as a homogenous macronutrient. As previously discussed, different proteins have specific amino acid abundances, and derive their constituent amino acids from different resources in the natural environment. We contend that assuming that all nutritionally-adequate proteins contain the same proportion of bioavailable amino acids with the same range of isotopic variability among individual δ13CAA values underestimates the actual complexity of different dietary protein sources and their metabolism.

Alternatively, we propose that the relationship between the amount of marine protein in the pig feed and the Δ13Ctissue – whole diet and Δ13Ctissue – dietary protein values results from preferential routing of non-essential amino acids, particularly glycine, from diet to tissue with increasing fish meal consumption (sensu Corr et al., Citation2005). The biosynthesis of non-essential amino acids occurs via enzyme-catalysed reactions and, in some cases, as many as three enzymes are required for synthesis. Biosynthesis can, however, be inhibited by an over-abundance of the end-product, largely because it is more energetically-efficient to make use of available amino acids than it is to make new ones. The enzyme utilised in the first irreversible reaction of glycine synthesis will not be produced in quantity if there is a high concentration of glycine readily available for tissue growth and repair (Schwarcz, Citation2000; Umbarger, Citation1978). The bioavailability in the dietary protein source of a particular non-essential amino acid could thus influence the dominance of biosynthesis vs. routing. Glycine constitutes approximately 17% of collagen-carbon and ∼one-third of all amino acid residues in collagen (Herring, Citation1972). There is a higher concentration of glycine in marine protein than in most kinds of terrestrial protein, including soy,Footnote4 and the δ13C values of marine glycine are typically higher than those of C3 or C4 plants and their consumers (Corr et al., Citation2005; Fantle et al., Citation1999; Hare et al., Citation1991; Howland, Citation2003; Keil & Fogel, Citation2001; Young, Citation2002). Preferential routing of non-essential amino acids from diet is thought to occur simply under conditions of higher protein consumption (e.g. as recently discussed in Fernandes et al. (Citation2012), regardless of the type of protein.

Stable isotopic analysis of amino acids from archaeological human and faunal and controlled feeding study samples provides some support for the preferential routing explanation (e.g. Corr et al., Citation2005; Howland, Citation2003; Howland et al., Citation2003; Jim et al., Citation2006; Jones, Citation2002; Young, Citation2002). For rats fed on diets containing different energy and protein resources (e.g. C3 carbohydrates and lipids with C4 protein), Jim et al. (Citation2006) determined that, when dietary protein is present in excess of basic requirements, a non-trivial proportion of non-essential amino acids are routed with minimal fractionation from the protein component of diet (e.g. glycine ∼43%, aspartate ∼28%). It would therefore be expected that tissue non-essential amino acid carbon isotopic compositions would approach the corresponding feed isotopic compositions with higher protein content diets, and this did indeed occur for the controlled feeding study with rats. Crucially, however, the kind of protein (milk casein) was the same for all diets (i.e. only the carbon-isotope compositions were different), so the results of this controlled feeding study (Jim et al., Citation2006) are not directly applicable here. A subset of carbon isotopic data from the controlled feeding study conducted by Harvard University and the United States Department of Agriculture (Young, Citation2002; see also Howland, Citation2003; Howland et al., Citation2003) is more relevant to elucidating the mechanism that underlies our findings. A comparison of the δ13CAA values from two of the diets in that study, specifically, the maize-soybean and maize-fish meal diets, reveals distinct differences in the Δ13Ctissue – whole diet isotopic offsets for glycine. For the exclusively terrestrial protein diet, the Δ13CGly tissue – Gly whole diet offset was +5.2 ± 2.5‰, and was +0.8 ± 2.6‰ for the marine protein diet. Preliminary compound-specific carbon-isotope compositions for muscle and collagen from the 10 sows in this controlled feeding study further demonstrate that the relative predominance of glycine biosynthesis vs. direct routing gradually and linearly increases as the proportion of marine dietary protein increases (). We contend that this change in carbon isotopic discrimination indicates a shift from glycine biosynthesis to direct routing caused by the change in the nature of dietary protein, rather than protein quantity or quality, and that it becomes significant after approximately 50% of the dietary protein is marine-derived. The Δ13Ctissue AA – AA whole diet relationships of other non-essential amino acids (e.g. proline and glutamic acid) also appear to be linearly correlated to the proportion of marine protein in diet. Glycine is likely to have a greater impact on bulk tissue δ13C values, however, because both the Δ13CGly tissue – Gly whole diet offsets and the absolute δ13CGly values change significantly with increasing marine protein consumption. Future compound-specific isotopic analysis of the University of Bristol controlled feeding study pigs will provide further insight into the nature of this proposed change in glycine metabolism through the analysis of multiple tissues and larger numbers of pigs per diet group.

Conclusion

Stable carbon-isotope analysis of archaeological human and animal remains is routinely used to investigate palaeodiet, and is a particularly valuable tool in ecologically-complex regions where many different classes of dietary resources may have been consumed. In recent years, however, it has become increasingly evident that many of the fundamental assumptions about the biochemical relationships between consumer tissues and dietary intake are poorly understood. In this study, we investigated the impact of marine protein consumption on consumer tissue isotopic compositions through a controlled pig feeding study. This study has shown that, when a mixed diet incorporating marine resources is consumed, bulk tissue δ13C values may be misleading. Further, without considerable a priori knowledge about the relative proportions of protein sources consumed, it is difficult to accurately relate tissue carbon-isotope compositions to food resource isotopic data from the natural environment.

We have quantitatively demonstrated that bulk carbon-isotope values must be interpreted with caution when reconstructing palaeodiet and that any attempt at quantifying the proportion of marine/terrestrial resource consumption beyond the broadest level must be very carefully considered indeed. We have shown that even a small change in the amount of marine protein consumed can induce large changes in the relationships between tissue and whole diet and tissue and dietary protein, with comparatively small changes in the δ13Ctissue value. Tissue–whole diet carbon-isotope discrimination changed significantly and was correlated with dietary protein source, specifically with increasing marine protein consumption. This outcome implies that our current understanding of the underlying amino acid isotopic composition and how it relates to the more commonly-used bulk isotopic compositions is poor, and that solely using bulk isotopic compositions not only masks considerable information about diet but may, in fact, present an erroneous picture of resource consumption. Using preliminary amino acid carbon-isotope compositions, we argue that increased routing of non-essential amino acids, especially isotopically-heavy amino acids like glycine, is the primary cause of this dynamic tissue–whole diet offset when marine protein is consumed. Ongoing stable isotope analysis of individual amino acids will determine if this is indeed the mechanism driving the increasing Δ13Ctissue – whole diet isotopic offset.

Supplemental material

Supplementary_Tables.docx

Download MS Word (34.6 KB)

Acknowledgements

This work was conducted as part of a responsive mode grant from the UK Natural Environment Research Council (NERC: NE/D004535/1). The authors would also like to thank Miss Helen Fewlass (Department of Archaeology and Anthropology, University of Bristol) and Dr. Jason Newton (NERC Life Sciences and Mass Spectrometry Facility, East Kilbride node).

Notes

1. As opposed to compound-specific isotopic compositions, e.g. of individual amino acids.

2. Excluding piglet W15 for which only femoral bone collagen could be extracted.

3. Loin muscle samples were only available for eight piglets.

4. Fish meal contains approximately 2.3 times as much glycine as soymeal, i.e. 10.2g/16g N vs. 4.5g/16g N, where 16g N ≈ 100g of protein (Jorgensen et al., Citation1984).

References

  • Ambrose, S.K. 1993. Isotopic Analysis of Paleodiets: Methodological and Interpretive Considerations. In: M.K. Sandford, ed. Investigations of Ancient Human Tissue: Chemical Analyses in Anthropology. Langhorne, PA: Gordon and Breach Science Publishers, pp. 59–130.
  • Ambrose, S.K. & Norr, L 1993. Experimental Evidence for the Relationship of the Carbon Isotope Ratios of Whole Diet and Dietary Protein to Those of Bone Collagen and Carbonate. In: J.B. Lambert and G. Grupe, eds. Prehistory Human Bone: Archaeology at the Molecular Level. London: Springer-Verlag, pp. 1–38.
  • Ambrose, S.H., Butler, B.M., Hanson, D.B., Hunter-Anderson, R.L. & Krueger H.W. 1997. Stable Isotopic Analysis of Human Diet in the Marianas Archipelago, Western Pacific. American Journal of Physical Anthropology, 104: 343–61. doi: 10.1002/(SICI)1096-8644(199711)104:3<343::AID-AJPA5>3.0.CO;2-W
  • Brand, W.A., Coplen, T.B., Vogl, J., Rosner, M. & Prohaska, T. 2014. Assessment of International Reference Materials for Isotope-Ratio Analysis (IUPAC Technical Report). Pure and Applied Chemistry, 86: 425–67. doi: 10.1515/pac-2013-1023
  • Codron, D., Codron, J., Lee-Thorp, J.A., Sponheimer, M., de Ruiter, D., Sealy, J., Grant, R. & Fourie, N. 2007. Diets of Savannah Ungulates From Stable Carbon Isotope Composition of Faeces. Journal of Zoology, 273: 21–29. doi: 10.1111/j.1469-7998.2007.00292.x
  • Corr, L.T, Sealy, J.C, Horton, M.C. & Evershed, R.P. 2005. A Novel Marine Indicator Utilizing Compound-Specific Bone Collagen Amino Acid δ13C Values of Ancient Humans. Journal of Archaeological Science, 32: 321–30. doi: 10.1016/j.jas.2004.10.002
  • Cox, G. & Sealy, J. 1997. Investigating Identity and Life History: Isotopic Analysis and Historical Documentation of Slave Skeletons Found in the Cape Town Foreshore, South Africa. International Journal of Historical Archaeology, 1: 207–24. doi: 10.1023/A:1027349115474
  • DeNiro, M.J. & Epstein, S. 1976. You are What you eat (Plus A few ‰): the Carbon Isotope Cycle in Food Chains. Geological Society of America Abstracts with Programs, 8: 834–5.
  • DeNiro, M.J. & Epstein, S. 1978. Influence of Diet on the Distribution of Carbon Isotopes in Animals. Geochimica et Cosmochimica Acta, 42: 495–506. doi: 10.1016/0016-7037(78)90199-0
  • Druker, D. & Bocherens, H. 2004. Carbon and Nitrogen Stable Isotopes as Tracers of Change in Diet Breadth During Middle and Upper Palaeolithic in Europe. International Journal of Osteoarchaeology, 14: 162–77. doi: 10.1002/oa.753
  • Dufour, E., Bocherens, H. & Mariotti, A. 1999. Palaeodietary Implications of Isotopic Variability in Eurasian Lacustrine Fish. Journal of Archaeological Science, 26: 617–27. doi: 10.1006/jasc.1998.0379
  • Emery, K.F., Wright, L.E. & Schwarcz, H. 2000. Isotopic Analysis of Ancient Deer Bone: Biotic Stability in Collapse Period Maya Land-use. Journal of Archaeological Science, 27: 537–50. doi: 10.1006/jasc.1999.0491
  • Fantle, M.S., Dittel, A.I., Schwalm, S.M., Epifanio, C. E. & Fogel, M. L. 1999. A food web analysis of the juvenile blue crab, Callinectes sapidus , using stable isotopes in whole animals and individual amino acids. Oecologia, 120: 416–26. doi: 10.1007/s004420050874
  • Fernandes, R., Nadeau, M.J. & Grootes, P.M. 2012. Macronutrient-based Model for Dietary Carbon Routing in Bone Collagen and Bioapatite. Archaeological and Anthropological Sciences, 4: 291–301. doi: 10.1007/s12520-012-0102-7
  • Fernandes, R., Millard, A.R., Brabec, M., Nadeau, M.J. & Grootes, P. 2014. Food Reconstruction Using Isotopic Transferred Signals (FRUITS): A Bayesian Model for Diet Reconstruction. PLoS ONE, 9: e87436. doi: 10.1371/journal.pone.0087436
  • Finucane, B., Agurto, P.M. & Isbell, W.H. 2006. Human and Animal Diet at Conchopata, Peru: Stable Isotope Evidence for Maize Agriculture and Animal Management Practices During the Middle Horizon. Journal of Archaeological Science, 33: 1766–76. doi: 10.1016/j.jas.2006.03.012
  • Froehle, A.W., Kellner, C.M. & Schoeninger, M.J. 2010. FOCUS: Effect of Diet and Protein Source on Carbon Stable Isotope Ratios in Collagen: Follow up to Warinner and Tuross (2009). Journal of Archaeological Science, 37: 2662–70. doi: 10.1016/j.jas.2010.06.003
  • Froehle, A.W., Kellner, C.M. & Schoeninger, M.J. 2012. Multivariate Carbon and Nitrogen Stable Isotope Model for the Reconstruction of Prehistoric Human Diet. American Journal of Physical Anthropology, 147: 352–69. doi: 10.1002/ajpa.21651
  • Gil, A.F., Neme, G.A. & Tykot, R.H. 2011. Stable Isotopes and Human Diet in Central Western Argentina. Journal of Archaeological Science, 38: 1395–404. doi: 10.1016/j.jas.2011.01.010
  • Hare, P., Fogel, M., Stafford, T., Mitchell, A. & Hoering, T. 1991. The Isotopic Composition of Carbon and Nitrogen in Individual Amino Acids Isolated From Modern and Fossil Proteins. Journal of Archaeological Science, 18: 211–92.
  • Hedges, R.E.M. 2004. Isotopes and red Herrings: Comments on Milner et al. and Lidén et al. Antiquity, 78: 34–37. doi: 10.1017/S0003598X00092905
  • Hedges, R.E.M. & Reynard, L.M. 2007. Nitrogen Isotopes and the Trophic Level of Humans in Archaeology. Journal of Archaeological Science, 34: 1240–51. doi: 10.1016/j.jas.2006.10.015
  • Hedges, R.E., Clement, J.G., Thomas, C.D.L. & O'Connell, T.C. 2007. Collagen Turnover in the Adult Femoral mid-Shaft: Modelled From Anthropogenic Radiocarbon Tracer Measurements. American Journal of Physical Anthropology, 133: 808–16. doi: 10.1002/ajpa.20598
  • Heinritz, S.N., Mosenthin, R. & Weiss, E. 2013. Use of Pigs as a Potential Model for Research Into Dietary Modulation of the Human gut Microbiota. Nutrition Research Reviews, 26: 191–209. doi: 10.1017/S0954422413000152
  • Herring, G.M. 1972. The Organic Matrix of Bone. In: G.H. Bourne, ed. The Biochemistry and Physiology of Bone. Academic Press, pp. 127–89.
  • Hill, P.A. & Orth, M. 1998. Bone Remodelling. British Journal of Orthodontics, 25: 101–7. doi: 10.1093/ortho/25.2.101
  • Howland, M.R. 2003. Compound-specific Stable Isotope Investigations of the Influence of Diet on the Stable Isotope Composition of Body Tissues. PhD dissertation, University of Bristol.
  • Howland, M.R., Corr, L.T., Young, S.M.M., Jones, V., Jim, S., van der Merwe, N.J., Mitchell, A.D. & Evershed, R.P. 2003. Expression of the Dietary Isotope Signal in the Compound-Specific δ13C Values of pig Bone Lipids and Amino Acids. International Journal of Osteoarchaeology, 13: 54–65. doi: 10.1002/oa.658
  • Jauoen, K., Pons, M.L. & Balter, V. 2013. Iron, Copper and Zinc Isotopic Fractionation up Mammal Trophic Chains. Earth and Planetary Science Letters 374: 164–72. doi: 10.1016/j.epsl.2013.05.037
  • Jim, S., Jones, V., Ambrose, S.H., Evershed, R.P. 2006. Quantifying Dietary Macronutrient Sources of Carbon for Bone Collagen Biosynthesis Using Natural Abundance Stable Isotope Analysis. British Journal of Nutrition, 95: 1055–62. doi: 10.1079/BJN20051685
  • Jones, V. 2002. Investigating the Routing and Synthesis of Amino Acids Between Diet and Bone Collagen via Feeding Experiments and Applications to Palaeodietary Reconstruction. PhD dissertation, School of Chemistry, University of Bristol.
  • Jorgensen, H., Sauer, W.C. & Thacker, P.A. 1984. Amino Acid Availabilities in Soybean Meal, Sunflower Meal, Fish Meal and Meat and Bone Meal fed to Growing Pigs. Journal of Animal Science, 58: 926–34. doi: 10.2527/jas1984.584926x
  • Keil, R.G. & Fogel, M.L. 2001. Reworking of Amino Acid in Marine Sediments: Stable Carbon Isotopic Composition of Amino Acids in Sediments Along the Washington Coast. Limnology and Oceanography, 46: 14–26. doi: 10.4319/lo.2001.46.1.0014
  • Kellner, C.M. & Schoeninger, M.J. 2007. A Simple Carbon Isotope Model for Reconstructing Prehistoric Human Diet. American Journal of Physical Anthropology, 133: 1112–27. doi: 10.1002/ajpa.20618
  • Knudson, K.J., Peters, A.H. & Cagigao, E.T. 2015. Paleodiet in the Paracas Necropolis of Wari Kayan: Carbon and Nitrogen Isotope Analysis of Keratin Samples From the South Coast of Peru. Journal of Archaeological Science, 55: 231–43. doi: 10.1016/j.jas.2015.01.011
  • Knudson, K.J., Williams, H.M., Buikstra, J.E., Tomczak, P.D., Gordon, G.W. & Anbar, A.D. 2010. Introducing δ88/86Sr Analysis in Archaeology: A Demonstration of the Utility of Strontium Isotope Fractionation in Palaeodietary Studies. Journal of Archaeological Science, 37: 2352–64. doi: 10.1016/j.jas.2010.04.009
  • Krueger, H.W. & Sullivan, C.H. 1984. Models for Carbon Isotope Fractionation Between Diet and Bone. In: J. Turnlund and P.E. Johnson, eds. Stable Isotopes in Nutrition, ACS Symposium Series. Washington, DC: American Chemical Society, pp. 205–22.
  • Krummen, M., Hilkert, A.W., Juchelka D., Duhr, A., Schlüter, H.-J. & Pesch., R. 2004. A new Concept for Isotope Ratio-Monitoring Liquid Chromatography/Mass Spectrometry. Rapid Communications in Mass Spectrometry, 18: 2260. doi: 10.1002/rcm.1620
  • Litten-Brown, J.C., Corson, A.M. & Clarke, L. 2010. Porcine Models for the Metabolic Syndrome, Digestive and Bone Disorders: A General Overview. Animal, 4: 899–920. doi: 10.1017/S1751731110000200
  • Longin, R. 1971. New Method of Collagen Extraction for Radiocarbon Dating. Nature, 230: 241–6. doi: 10.1038/230241a0
  • Marino, B.D. & McElroy, M.B. 1991. Isotopic Composition of Atmospheric CO2 Inferred From Carbon in C4 Plant Cellulose. Nature, 349: 127–31. doi: 10.1038/349127a0
  • Müldner, G. & Richards, M.P. 2005. Fast or Feast: Reconstructing Diet in Later Medieval England by Stable Isotope Analysis. Journal of Archaeological Science, 32: 39–48. doi: 10.1016/j.jas.2004.05.007
  • Nehlich, O., Borić, D., Stefanović, S. & Richards, M.P. 2010. Sulphur Isotope Evidence for Freshwater Fish Consumption: A Case Study From the Danube Gorges, SE Europe. Journal of Archaeological Science, 37: 1131–9. doi: 10.1016/j.jas.2009.12.013
  • Nehlich, O., Fuller, B.T., Jay, M., Mora, A., Nicholson, R.A., Smith, C.I. & Richards, M.P. 2011. Application of Sulphur Isotope Ratios to Examine Weaning Patterns and Freshwater Fish Consumption in Roman Oxfordshire, UK. Geochimica et Cosmochimica Acta, 75: 4963–77. doi: 10.1016/j.gca.2011.06.009
  • Newsome, S.D., Phillips, D.L., Culleton, B.J., Guilderson, T.P. & Koch, P.L. 2004. Dietary Reconstruction of an Early to Middle Holocene Human Population From the Central California Coast: Insights From Advanced Stable Isotope Mixing Models. Journal of Archaeological Science, 31: 1101–15. doi: 10.1016/j.jas.2004.02.001
  • Passey, B.H., Robinson, T.F., Aylifee, L.K., Cerling, T.E., Sponheimer, M., Dearing, M.D., Roeder, B.L. & Ehleringer, J.R. 2005. Carbon Isotope Fractionation Between Diet, Breath CO2, and Bioapatite in Different Mammals. Journal of Archaeological Science, 32: 1459–70. doi: 10.1016/j.jas.2005.03.015
  • Patience, J.F., Thacker, P.A. & de Lange, C.F.M. 1995. Swine Nutrition Guide, 2nd ed. Saskatoon: Prairie Swine Centre, Inc.
  • Reynard, L.M., Henderson, G.M. & Hedges, R.E.M. 2010. Calcium Isotope Ratios in Animal and Human Bone. Geochimica et Cosmochimica Acta, 74: 3735–50. doi: 10.1016/j.gca.2010.04.002
  • Reynard, L.M. & Tuross, N. 2015. The Known, the Unknown and the Unknowable: Weaning Times From Archaeological Bones Using Nitrogen Isotope Ratios. Journal of Archaeological Science, 53: 618–25. doi: 10.1016/j.jas.2014.11.018
  • Richards, M.P., Schulting, R.J. & Hedges, R.E.M. 2003. Sharp Shift in Diet at Onset of Neolithic. Nature 425: 366. doi: 10.1038/425366a
  • Sayle, K.L., Cook, G.T., Ascough, P.L., Hastie, H.R., Einarsson, A., McGovern, T.H., Hicks, M.T., Edwald, A. & Friðriksson, A. 2013. Application of 34S Analysis for Elucidating Terrestrial, Marine and Freshwater Ecosystems: Evidence of Animal Movement/Husbandry Practices in an Early Viking Community Around Lake Mývatn, Iceland. Geochimica et Cosmochimica Acta, 120: 531–44. doi: 10.1016/j.gca.2013.07.008
  • Scherer, A.K., Wright, L.E. & Yoder, C.J. 2007. Bioarchaeological Evidence for Social and Temporal Differences in Diet at Piedras Negras, Guatemala. Latin American Antiquity, 18: 85–104. doi: 10.2307/25063087
  • Schwarcz, H.P. 2000. Some Biochemical Aspects of Carbon Isotopic Paleodiet Studies. In: S.H. Ambrose and M.A. Katzenberg, eds. Biogeochemical Approaches to Palaeodietary Analysis. London: Kluwer Academic Publishers, pp. 189–209.
  • Sponheimer, M., Robinson, T., Ayliffe, L., Passey, B., Roeder, B., Shipley, L., Lopez, E., Cerling, T., Dearing, D. & Ehleringer, J. 2003. An Experimental Study of Carbon-Isotope Fractionation Between Diet, Hair, and Faeces of Mammalian Herbivores. Canadian Journal of Zoology, 81: 871–6. doi: 10.1139/z03-066
  • Sponheimer, M., Robinson, T.F., Cerling, T.E., Tegland, L., Roeder, B.L., Ayliffe, L., Dearing, M.D. & Ehleringer, J.R. 2006. Turnover of Stable Carbon Isotopes in the Muscle, Liver, and Breath CO2 of Alpacas (Lama Pacos). Rapid Communications in Mass Spectrometry, 20: 1395–9. doi: 10.1002/rcm.2454
  • Sullivan, T.P., Eaglstein, W.H., Davis, S.C. & Mertz, P. 2001. The pig as A Model for Human Wound Healing. Wound Repair and Regeneration, 9: 66–76. doi: 10.1046/j.1524-475x.2001.00066.x
  • Swindle, M.M., Makin, A., Herron, A.J., Clubb Jr, F.J. & Frazier, K.S. 2012. Swine as Models in Biomedical Research and Toxicology Testing. Veterinary Pathology, 49: 344–56. doi: 10.1177/0300985811402846
  • Taylor, G., Tucker,K., Butler, R., Pike, A.W.G., Lewis, J., Roffey, S., Marter, P., Lee, O., Wu, H.H.T., Minnikin, D.E., Besra, G.S., Sing, P., Cole, S.T. & Stewart, G.R., 2013. Detection and Strain-Typing of Ancient Mycobacterium Leprae From A Medieval Leprosy Hospital. PLoS One, 8 (4): e62406. doi: 10.1371/journal.pone.0062406
  • Tauber, H. 1981. 13C Evidence for Dietary Habits of Prehistoric man in Denmark. Nature, 292: 332–3. doi: 10.1038/292332a0
  • Tieszan, L.L. & Fagre, T. 1993. Effect of Diet Quality and Composition on the Isotopic Composition of Respiratory CO2, Bone Collagen, Bioapatite, and Soft Tissues. In: J. Lambert and G. Grupe, eds. Molecular Archaeology of Prehistoric Human Bone. Berlin: Springer-Verlag, pp. 123–35.
  • Tomczak, P.D. 2003. Prehistoric Diet and Socioeconomic Relationships Within the Osmore Valley of Southern Peru. Journal of Anthropological Archaeology, 22: 262–78. doi: 10.1016/S0278-4165(03)00039-4
  • Turner, B.L., Kamenov, G.D., Kingston, J.D. & Armelagos, G.J. 2009. Insights Into Immigration and Social Class at Machu Picchu, Peru Based on Oxygen, Strontium, and Lead Isotopic Analysis. Journal of Archaeological Science, 36: 317–32. doi: 10.1016/j.jas.2008.09.018
  • Urey, H., 1947. The Thermodynamic Properties of Isotopic Substances. Journal of the Chemical Society: 562–81.
  • Umbarger, H.E. 1978. Amino Acid Biosynthesis and its Regulation. Annual Review of Biochemistry, 47: 533–606. doi: 10.1146/annurev.bi.47.070178.002533
  • Van der Merwe, N.J. & Vogel, J.C. 1978. 13C Content of Human Collagen as A Measure of Prehistoric Diet in Woodland North America. Nature, 276: 815–6. doi: 10.1038/276815a0
  • Warinner, C., Garcia, N.R. & Tuross, N. 2013. Maize, Beans and the Floral Isotopic Diversity of Highland Oaxaca, Mexico. Journal of Archaeological Science, 40: 868–73. doi: 10.1016/j.jas.2012.07.003
  • Webb, E., White, C. & Longstaffe, F. 2013. Dietary Shifting in the Nasca Region as Inferred From the Carbon- and Nitrogen-Isotope Compositions of Archaeological Hair and Bone. Journal of Archaeological Science, 40: 129–39. doi: 10.1016/j.jas.2012.08.020
  • White, C.D. 2005. Gendered Food Behaviour among the Maya: Time, Place, Status and Ritual. Journal of Social Archaeology, 5: 356–82. doi: 10.1177/1469605305057572
  • White, C.D., Healy, P.F. & Schwarcz, H.P. 1993. Intensive Agriculture, Social Status and Maya Diet at Pacbitun, Belize. Journal of Anthropological Research, 49: 347–75. doi: 10.1086/jar.49.4.3630154
  • Yesner, D.R., Torres, M.J.F., Guichon, R.A. & Borrero. L.A. 2003. Stable Isotope Analysis of Human Bone and Ethnohistoric Subsistence Patterns in Tierra del Fuego. Journal of Anthropological Archaeology, 22: 279–91. doi: 10.1016/S0278-4165(03)00040-0
  • Young, S.M.M. 2002. Metabolic mechanisms and the isotopic reconstruction of ancient diets with an application on remains from Cuello, Belize. PhD dissertation, Harvard University.