Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Stable isotope analysis of human bone from Ganj Dareh, Iran, ca. 10,100 calBP

  • Deborah C. Merrett ,

    Roles Conceptualization, Data curation, Formal analysis, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    dcmerret@sfu.ca (DCM); christina.t.t.cheung@googlemail.com (CC)

    Affiliation Department of Archaeology, Simon Fraser University, Burnaby, BC, Canada

  • Christina Cheung ,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Project administration, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    dcmerret@sfu.ca (DCM); christina.t.t.cheung@googlemail.com (CC)

    Affiliations EA–Eco-anthropologie (UMR 7206), Muséum National d’Histoire Naturelle, CNRS, Université Paris Diderot, Paris, France, UMR 7269, LAMPEA, Aix-Marseille Université, CNRS, Minist Culture, Aix-en-Provence, France

  • Christopher Meiklejohn,

    Roles Validation, Writing – original draft, Writing – review & editing

    Affiliation Department of Anthropology, University of Winnipeg, Winnipeg, MB, Canada

  • Michael P. Richards

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Software, Supervision, Validation, Writing – review & editing

    Affiliation Department of Archaeology, Simon Fraser University, Burnaby, BC, Canada

Abstract

We report here on stable carbon, nitrogen, and sulfur isotope values from bone collagen of human (n = 20) and faunal (n = 11) remains from the Early Neolithic site of Ganj Dareh, Iran, dating to ca. 10,100 cal. BP. Our focus explores how isotope values of human bone vary by age and sex, and evaluates dietary practices at this site. It also provides a baseline for future studies of subsistence in the early Holocene Central Zagros Mountains, from the site with the first evidence for human ovicaprid management in the Near East. Human remains include individuals of all age groups for dietary reconstruction, as well two Ottoman intrusive burials for temporal and cultural comparison. All analyzed individuals exhibited δ13C and δ15N values consistent with a diet based heavily on C3 terrestrial sources. There is no statistically significant difference between the isotopic compositions of the two sexes, though males appear to show larger variations compared to females. Interesting patterns in the isotopic compositions of the subadults suggested weaning children may be fed with supplements with distinctive δ13C values. Significant difference in sulfur isotope values between humans and fauna could be the earliest evidence of transhumance and could identify one older adult male as a possible transhumant shepherd. Both Ottoman individuals had distinctively different δ13C, δ15N, and δ34S values compared to the Neolithic individuals. This is the first large scale analysis of human stable isotopes from the eastern end of the early Holocene Fertile Crescent. It provides a baseline for future intersite exploration of stable isotopes and insight into the lifeways, health, and processes of neolithisation associated with the origins of goat domestication at Ganj Dareh and the surrounding Central Zagros.

Introduction

Although the Levant and Anatolia have received substantial attention concerning the origins of agriculture, much less has been reported on the easternmost Fertile Crescent, the Central Zagros. The ‘accepted’ story for this region revolves around one of the key processes of neolithisation: domestication of ovicaprids, and especially goat, a transformation suggested to have been unwittingly orchestrated by altered interactions between humans and local ungulate populations [14].

While recent excavations at the slightly earlier site of Sheik-e Abad address origins of herding through identification of penning and dung [5, 6], previous research was formulated using ovicaprid age and sex profiles as a proxy for human-goat interactions [2, 4, 7]. The first evidence for human control of goats comes from Ganj Dareh, where Hesse saw a predominance of young male goats/sheep in the deposits, i.e. slaughter pattern analysis, as supporting the hypothesis that the animals were under human control, managed rather than hunted [4, 8, 9]. Goat hoof impressions in mud bricks [10] also suggested that, though morphologically wild, the Ganj Dareh goats lived in close proximity to the human settlement and, at least behaviourally, were on the way to domestication [4, 8, 9]. Recent radiocarbon dates confirm contemporaneity of the human population [11] and ovicaprids [4] at the site, emphasizing the need to further explore life in the Early Neolithic of the Central Zagros. This study provides a new avenue for elucidating human lifeways at the site, with the first stable isotope analysis of the Ganj Dareh human remains. Published data from nearby contemporaneous sites are also included to provide a more comprehensive picture of dietary practices in the region.

Introduction to the site of Ganj Dareh

Ganj Dareh Tepe, a small mound or Tepe in the High Zagros of Kermanshah Province in Western Iran, is one of several in the area (Fig 1). Lying in the Gamas-Ab Valley at an altitude of ~1400 m [12], it measures ~40 m in diameter with 7 to 8 m of cultural deposits. The initial work identified five levels, A to E, with level E at the base and the first permanent architecture in level D, well preserved by an extensive fire. New research based on re-examination of the original field notes of Philip Smith, shows the stratigraphy of the site to be more complex than previously published [13] but does not alter any of the conclusions of this paper. During excavations sponsored by the University of Toronto and the Royal Ontario Museum, Philip E.L. Smith excavated roughly 20 percent of the mound in four seasons between 1967 and 1974 [10, 1416]. Current evidence places site occupation at 10,170–9,700 cal. BP, based on goat remains [4] and ca. 10,100 cal. BP based on human remains [11], confirming the relatively short occupation period of the site (100–200 years), and the contemporaneity of the humans and ovicaprids. Analysis of the full radiocarbon record also lays to rest the earlier idea of a major hiatus between levels D and E at the base of the site [4, 11].

thumbnail
Fig 1. Location of Ganj Dareh (1) and other sites of Central Zagros mentioned in the text: Abdul Hosein (2), Sheikh-e-Abad (3), Jani (4), and Wezmeh Cave (5).

Inset: Location of Central Zagros (CZ) and Çatalhöyük (ÇH). Scale = 30 km; inset scale = 150 km. Modified from https://www.visibleearth.nasa.gov/images/106461/the-middle-east.

https://doi.org/10.1371/journal.pone.0247569.g001

Carbonized macrobotanical remains were isolated during excavation [17], with two-row hulled barley (Hordeum spontaneum), a C3 crop, the only cereal isolated. However, general poor preservation made it clear that presence of low levels of non-shattering rachis morphology alone was not sufficient to identify domestication [18]; at Early Neolithic sites of the Euphrates Valley with excellent plant preservation, non-shattering internode appearance was identical to that of modern wild populations. Thus, while there is insufficient evidence for domestication of barley, presence of barley and grinding stones suggests, at a minimum, the use of wild plants [19].

As with the botanical remains [17], faunal proportions did not vary through the occupation sequence, suggesting that subsistence did not change during use of the site [2, 4, 11]. Though a broad faunal spectrum was in use, with 49 identifiable taxa; goats, sheep, fox and partridge predominated and, of identifiable mammalian bone, 89.9% were ovicaprids, with goat predominating [7]. All levels displayed slaughter patterns consistent with goat management rather than hunting [2].

Principals of stable isotope analysis

Stable carbon and nitrogen isotope analysis of bone collagen, expressed as δ13C and δ15N, is in routine use for dietary and mobility reconstruction [2023]. More recently, analysis of the stable sulfur isotope (δ34S) has become a significant tool in reconstructing ancient lifeways [24]. The stable carbon isotope ratio (13C/12C) distinguishes organisms exploiting marine, terrestrial, and plant based ecosystems, with their different photosynthetic pathways (e.g. C3 vs. C4) [2527]. The stable nitrogen isotope ratio (15N/14N) identifies trophic level effects in the food chain, with a 3–5‰ enrichment from one trophic level to another, making δ15N values particularly useful for estimating animal protein intake [28]. However, they are also impacted by physiological and environmental phenomena, such as breastfeeding [29], long-term starvation [30], manuring [31, 32], and aridity [33, 34]. The stable sulfur isotope ratio (34S/32S) is primarily influenced by geochemical conditions [35] and, while highly variable and locally specific, can be used to identify migrants from regions with distinctive δ34S baselines. The primary sources of sulfur in bone are dietary methionine, and its internal recycling. Methionine is an essential amino acid. Thus, δ34S values reflect methionine from dietary protein [24], and may also contribute to health assessment by providing insight into nitrogen metabolism and stress exposure.

Materials

Samples were collected from 32 human and 25 faunal bones (Table 1). 30 human samples are Early Neolithic (ca. 10,100 cal. BP). An additional two (GD# 1150 and GD# 1151–1), surface burials of Ottoman age (450–320 years cal. BP) [11], are included for the record. The human remains were from burials beneath house floors, in the rubble that accumulated in collapsed houses and, in one case, from a test pit outside the mound.

thumbnail
Table 1. Human and faunal samples analysed showing level in the site and burial context relative to architectural features.

https://doi.org/10.1371/journal.pone.0247569.t001

In general, the human bone collection is highly fragmentary, leading to differing MNI estimates over the course of the reconstruction and analysis of the collection. Factors included identification of multiple burials over several rounds of analysis and recovery of 17 isolated single human bones/bone fragments/single teeth from the faunal sample. The current MNI is 116, with 56 catalogued as skeletons, defined as those with >4 skeletal elements [36]. Of these, 52 were identified during the excavation process [40]. For biological and metabolic consistency of bone used and minimal bone destruction, ribs were only sampled from individuals classified as skeletons (N = 32 with ribs present). No permits were required for this study, which complied with all relevant regulations. Samples are curated by one of us (DCM) at Department of Archaeology, Simon Fraser University. The faunal remains used in this study were included in the work of Hesse (7–9), then associated in the collection with the human assemblage as fauna accompanying the human burials. Although probably not intentionally associated with the human remains, they none-the-less were within the same stratigraphic level and in close enough proximity to be identified as with the human remains.

The internal site structure gave no evidence for courtyards or roads; remains were recovered from within architectural features, primarily below house floors or in bricked features [10]. A large expanse of level D was exposed to fire, preserving mud brick walls up to a metre in height, with some below-house-floor burials sufficiently exposed to heat to cause bone calcination [10]. This reduced the number of successful collagen extractions in the present study.

Methods

Osteological analysis

Ages-at-death estimations for the human skeletal remains are by Merrett [36, 37]. Non-adult ages were based on multiple methods: tooth development and eruption [41, 42], enamel cross-striation age [43], bone developmental patterns [44], bone length [4449], and patterns of epiphyseal fusion [44, 50, 51]. Adult age estimation methods included tooth wear [52, 53], pubic symphysis [54], auricular surface morphology [55], and cranial and palatal suture closure [56, 57]. Adult sex estimates were based on sexual dimorphism of cranial and pelvic bones [50, 5860]. Sex of six individuals was determined through aDNA analysis [38, 39]. Bone samples for stable isotope analysis were from rib fragments. In addition, although most human remains were from below-house floor contexts, two were identified as from midden-like sediments: one within roof rubble, the other outside the settlement (Table 1).

Faunal remains used here were excavated as ‘associated’ with the human remains, and part of the sample analysed and placed in the unidentified category by Hesse [7]. The original excavation notes give no evidence for ritual interment of fauna. Although the fauna were too fragmented for identification, all used were from the ’medium mammal long bone’ category, and were most likely ovicaprid. As a result, the faunal isotopic results can only serve as a regional baseline, and cannot be used to evaluate the dietary compositions of the Ganj Dareh residents.

Stable isotope analysis

Bone samples were prepared at the Isotope Chemistry Laboratory, Simon Fraser University; extracted collagen was sent to Isoanalytical Limited (UK) for analysis. Collagen was extracted following the modified Longin method [61] and additional ultrafiltration. Sample analysis was with a Europa ScientificTM elemental analyzer, coupled to a mass spectrometer. All carbon and nitrogen isotope values are averaged, based on duplicate analysis, and reported in per mil (‰). Results are calibrated to VPDB and AIR, respectively, using international standards IAEA-N-1, IAEA-C7, and IAEA-CH-6. Accuracy of measurements are monitored using in-house check standards IA-R068 (soy protein, δ13C = –25.22‰, δ15N = +0.99‰), IA-R038 (L-alanine, δ13C = –24.99‰, δ15N = –0.65‰), IA-R069 (tuna protein, δ13C = –18.88‰, δ15N = +11.60‰), and a mixture of IA-R046 and IAEA-C7 (ammonium sulfate and oxalic acid, δ13C = –14.48‰, δ15N = +22.04‰). Averaged measured δ13C values for IA-R068 (n = 12), IA-R038 (n = 5), IA-R069 (n = 5), and IA-R046/IAEA-C7 (n = 5) are –25.23‰± 0.04‰, –25.04‰± 0.03‰, –18.90‰± 0.02‰, and –14.53‰± 0.06‰, respectively. Averaged measured δ15N values for IA-R068 (n = 12), IA-R038 (n = 5), IA-R069 (n = 5), and IA-R046/IAEA-C7 (n = 5) are +1.00‰± 0.04‰, –0.59‰± 0.04‰, +11.71‰± 0.04‰, and +21.89‰± 0.09‰, respectively.

For sulfur isotope analysis, all measurements are reported in per mil (‰). Results are calibrated to VCDT using international standards NBS-127, IAEA-S-1, and IAEA-SO-5. Accuracy of measurements are monitored using in-house check standards IA-R061 (barium sulfate, δ34S = +20.33‰), IA-R069 (tuna protein, δ34S = +18.91‰), and NBS-1577B (bovine liver, δ34S = +7.50‰). Averaged measured δ34S values for IA-R061 (n = 6), IA-R069 (n = 3), and NBS-1577B are +20.31‰±0.08‰, +18.69‰±0.31‰, and +7.61‰±0.20‰, respectively.

Collagen quality was assessed using conventional criteria: %collagen between 0.5% and 22% by weight, %C between 15.3% and 47%, %N between 5.5% and 17.3%, %S between 13% and 35%, atomic C/N ratio between 2.9 and 3.6, atomic C/S ratio between 300 and 900, and atomic N/S ratio between 100 and 300 [6167]. Only samples with elemental compositions within these ranges were accepted for analysis.

Results

With many samples heavily calcined, only 20 of 32 human and 11 of 25 faunal samples yielded sufficient preserved collagen for C and N measurements; 10 of 20 human and 10 of 11 faunal samples for S measurements. Results of accepted measurements, elemental compositions, and sample information are listed in Table 2.

thumbnail
Table 2. Summary of sample information, elemental compositions, and bone collagen carbon, nitrogen, and sulfur isotopic compositions of all samples analysed in this study.

https://doi.org/10.1371/journal.pone.0247569.t002

Human and faunal C, N, and S measurements are plotted in Fig 2, together with published human C and N data (n = 3) from the nearby and contemporary sites Abdul Hosein (AH) and Wezmeh Cave (WC) (2 and 5 in Fig 1). All humans have δ13C values between –20.4 and –18.2‰, δ15N between +8.7 and +14.4‰, and δ34S between +9.0 and +14.2‰. All fauna have δ13C values between –20.8 and –17.9‰, δ15N between +6.9 and +10.3‰, and δ34S between +11.5 and +14.2‰. δ13C and δ15N values of the three individuals from AH and WC fall well within the Ganj Dareh range. However, the two Ottoman intrusive burials (Fig 2) have distinctive stable isotopic compositions (C, N, and S) compared to the Ganj Dareh Neolithic group. Unfortunately, the Ottoman sample size is too small for meaningful discussion. S measurements are only available for Ganj Dareh (Fig 2B). The S compositions of the humans are significantly depleted in 34S compared to those of the fauna (Wilcoxon test: diff: 2.0‰, p = 0.01).

thumbnail
Fig 2. Human and faunal stable carbon, nitrogen, and sulfur isotope data from western Iran.

Closed circles (●) correspond to Neolithic samples, open circles (○) correspond to human samples dated to the Ottoman period: a) carbon and nitrogen isotope data of all humans and fauna from western Iran; data from Adbul Hosein and Wezmeh Cave are from [68]; b) sulfur isotope data of all humans and fauna from Ganj Dareh.

https://doi.org/10.1371/journal.pone.0247569.g002

Sex comparison

For this comparison, all “possible male” and “possible female” individuals are treated as “male” and “female”, respectively. Only Neolithic individuals above the age of 4 are considered, so that δ15N values will not be complicated by breastfeeding effect. δ13C, δ15N, and δ34S values do not significantly differ between males and females (Wilcoxon tests: δ13C: W = 10, p = 1.0; δ15N: W = 13, p = 0.5167; and δ34S: W = 5, p = 1.0) (Fig 3).

thumbnail
Fig 3.

Stable carbon, nitrogen, and sulfur isotope data of all sexed individuals from Ganj Dareh: a) carbon and nitrogen isotope data; b) nitrogen and sulfur isotope data; c) carbon and sulfur isotope data; and d) the sulfur isotope data.

https://doi.org/10.1371/journal.pone.0247569.g003

Non-adult comparison

Fig 4 shows the age profiles of stable carbon and nitrogen isotope values in non-adults (n = 9). Because we have no genetic evidence for mother/child relationships in the sample, the means of all adult females from the site (δ13C: –19.7‰; δ15N: +10.4‰, n = 3) are used as the baseline to evaluate the effect of breastfeeding in non-adults. Both δ13C and δ15N values decrease with increasing age.

thumbnail
Fig 4.

Carbon and nitrogen isotope data for Ganj Dareh non-adults (age < 15 years): a) carbon isotope values of non-adults relative to mean Neolithic adult female value (–19.7‰); b) nitrogen isotope values of non-adults relative to mean Neolithic adult female value (+10.4‰). Dotted lines correspond to mean female adult values.

https://doi.org/10.1371/journal.pone.0247569.g004

Discussion

The overall δ13C and δ15N human and fauna values from Ganj Dareh provide an initial baseline for the Early Neolithic of the Central Zagros and are consistent with the limited data available from nearby Adbul Hosein and Wezmeh Cave (Fig 2A). However, further research is needed before evaluation of the meanings of intersite variation can proceed. We need not expect the trajectory away from hunting and gathering to be linear, occur at similar rates, nor vary in similar ways between sites or even within one region [69]. These data are however invaluable in the exploration of the pathways to goat domestication in the Central Zagros, especially given that recent genetic findings [70] indicate the relative isolation of the region from processes of neolithisation elsewhere in the Fertile Crescent to the west.

Sex comparison

While no statistically significant difference exists between the means of the two sexes, the males have much larger ranges for all three isotopes than do the females (Fig 3 and Table 3). Given the small sample size, any interpretation of higher male variability is speculative without further evidence. Some possible explanations include: i) males had access to a larger variety of foods, ii) males had more access to ‘exotic’ foods for ceremonial purposes, and iii) males had higher mobility than females. More than one of these processes may be occurring.

thumbnail
Table 3. Comparing variability in Ganj Dareh isotopic compositions between males and females (age 4+) using total range and standard deviation.

https://doi.org/10.1371/journal.pone.0247569.t003

Stable carbon and nitrogen isotope analysis

The Ganj Dareh fauna (only one specimen identified to species, dog) have δ13C (–20.8‰ to –18.8‰) and δ15N values (+6.9‰ to +10.3‰) consistent with those from contemporary Wezmeh Cave (5 in Fig 1) [66], Jani (4 in Fig 1) and Sheikh-e-Abad (3 in Fig 1) [71]. This contrasts with contemporary Çatalhöyük, at a much lower elevation on the Turkish Konya Plain and with a less continental climate. In this case, maximum faunal δ13C values are substantially higher than the Zagros values, at –13.6‰ for sheep and –12.5‰ for goats. Pearson and colleagues [72] suggest presence of C4 plants such as wet/salt tolerant grasses in the Çatalhöyük diet. The Ganj Dareh data suggest that either C4 grasses were not present or were not utilized in the early Holocene. Following the δ13C trend, minimum δ15N values from Çatalhöyük fauna are similar to those from Zagros sites while the maximum value for sheep is substantially higher (+14.4‰) [72], again illustrating different resource availability or utilization in the two regions.

Turning to the human material, adult δ13C and δ15N values from Neolithic Ganj Dareh are also in keeping with the one individual from Wezmeh Cave and the two from Abdul Hosein (Fig 2A), suggesting similar resource access and utilization in the Central Zagros. All individuals have a δ13C value range between –20.4‰ and –19.2‰ (n = 17), values consistent with a terrestrial diet based primarily on C3 resources, as noted above.

In archaeology, weaning practices in past populations can be detected using stable carbon and nitrogen isotope analysis. Previous studies have assumed that the enrichment in nitrogen isotope ratios from mother to infant is between 2–5‰ [7377]. A smaller enrichment between 0.5 and 1.4‰ is also expected in the δ13C values of breastfeeding individuals [29]. In this study, similar pattern in both δ13C and δ15N values of the subadults can be seen (Fig 4), where the highest δ13C and δ15N values are reported amongst children below the age of 4. Unfortunately, our sample size is too small to allow for a more specific estimation of weaning age at the site.

There has been growing evidence that other than diet, high levels of stress exposure and malnutrition may also cause variability in the isotopic compositions of humans [29, 7880], mainly in δ15N values [79]. However, at Ganj Dareh, the corresponding patterns in both δ13C and δ15N values of the subadults suggested that diet was a stronger factor than metabolic disruption for the trend in these individuals’ isotopic compositions.

Stable sulfur isotope analysis

The difference in δ34S values between the Ganj Dareh human and faunal samples is quite unusual (Fig 2B). As mentioned earlier, δ34S values vary in different geological environments; δ34S values should reflect where food was grown and/or produced. These data suggest that, at Ganj Dareh, humans and fauna were exploiting different geographical regions. This could involve a specific animal husbandry regime, with animals raised in a different area, with a distinctive S isotope baseline compared to where plant foods were grown and/or collected. Alternative explanations include fauna being brought in from elsewhere for herd replacement following local faunal extinctions [36] and/or for funerary feasting events. A final possibility is that none of the faunal remains tested are from individuals or species that were part of the human diet with only one sample identified to species—dog, though the size profile of the faunal remains used and the frequency of various species demonstrated by Hesse [79] appears to make that unlikely.

Overall although sample size is small, similar to the patterns in δ13C and δ15N values, the male variation in δ34S values is much greater than that in females, suggesting possible sex differences. These could have origins in different food sources for men and women, both in regular diet and in ritual practices, matrilocal postmarital residence practices, or more variation in male adult travel. Consistent with the latter scenario, one adult male GD#34 has a very similar δ34S value to those of the fauna; he may have been closely associated with the animals. In the High Zagros, people practicing transhumance with goats could potentially spend much of the year in geologically different locations. Were this the case, then δ34S of GD#34 may be evidence during the early Neolithic for an early form of what later would become what is referred to as transhumance.

Conclusions

We present the first large scale stable isotope analysis on human bone collagen from Ganj Dareh, Central Zagros, Iran, ca. 10,100 cal. BP. Though limited results have been previously published for other sites in the region, sample sizes were much smaller. As a result, our data establish a much more solid data set for future inter-site comparisons of lifeways for the Early Neolithic period of this region, and to test the hypothesis of multi-linear trajectories in the transition to agricultural subsistence enunciated by Smith’s scenario of ’low-level’ food production economies [69]. The carbon isotope measurements indicate a terrestrial C3-based diet, with no evidence of C4 plant usage, in strong support for the conclusion of van Zeist et al. [17].

Although not statistically different for all isotopes analysed, increased variation of males over females points to possible difference in dietary breadth between the sexes. Two children, aged 7 months and 2.5–3 years, have elevated δ15N values. This could indicate that they had not been weaned, but other explanations are possible, including exposure to stress episodes, or use of "special" weaning supplements.

Finally, sulfur isotope analysis identified an individual with an outlier δ34S value, an older adult male. We see the adult as either spending a considerable portion of his later life in a geologically non-local region, or at least ingesting non-local foods, consistent with possible interpretation of this individual as a transhumant shepherd. If so, the origins of pastoral transhumance may be earlier than previously thought.

Within the regional context of early neolithisation and incipient goat domestication, only three individuals from two other sites were available for comparison. Carbon and nitrogen measurements from these are consistent with those from Ganj Dareh. Given the considerably larger sample size for the Ganj Dareh analysis, we see the combined results as suggesting similar subsistence and lifeways across the Early Neolithic Central Zagros. Further investigations are indeed needed. It remains possible that life in the Central Zagros was one of ’low-level’ food production, wherein the trajectory from hunting/gathering to herding and farming was neither uniform within the region nor unidirectional [69]. What is also clearly needed is a sense of how our understanding of the archaeological models for the Early Neolithic, within a broader perspective, are affected by the recent aDNA evidence suggesting that the Ganj Dareh population had little if any biological relationship with contemporary and later groups in the Tigris/Euphrates lowlands or further west in the Levant, but were tied to groups to the north and northeast such as in the Caucasus (see especially [70]).

References

  1. 1. Larson G, Fuller DQ. The evolution of animal domestication. Annu Rev Ecol Evol Syst. 2014;45: 115–136.
  2. 2. Zeder MA. A metrical analysis of a collection of modern goats (Capra hircus aegargus and Capra hircus hircus) from Iran and Iraq: implications for the study of caprine domestication. J Archaeol Sci. 2001;28: 61–79.
  3. 3. Zeder MA. The origins of agriculture in the Near East. Cur Anthropol. 2011;52(S4): S221–S235.
  4. 4. Zeder MA, Hesse B. The initial domestication of goats (Capra hircus) in the Zagros Mountains 10,000 years ago. Science. 2000;287: 2254–2257. pmid:10731145
  5. 5. Fazeli Nashli H, Matthews R. The neolitisation of Iran: Patterns of change and continuity. In: Matthews R, Faseli Nashli H, editors. The neolithisation of Iran: The formation of new societies. Oxford: Oxbow Books; 2013. pp. 1–13.
  6. 6. Matthews W, Shillito L-M, Elliott S, Bull ID, Williams J. Neolithic lifeways: Microstratigraphic traces within houses, animal pens, and settlements. Proc Br Acad. 2014;198: 251–279.
  7. 7. Hesse BC.1978. Evidence for husbandry from the Early Neolithic site of Ganj Dareh in western Iran. PhD dissertation, Columbia University. 1978.
  8. 8. Hesse BC.1982. Slaughter patterns and domestication: The beginnings of pastoralism in western Iran. Man (N.S.). 1982;17: 403–417.
  9. 9. Hesse BC. These are our goats: The Origins of herding in west central Iran. In: Clutton-Brock J, Grigson C, editors, Animals and archaeology: 3. Early herders and their flocks. Oxford: British Archaeological Reports International Series 202;1984. pp. 243–264.
  10. 10. Smith PEL. Architectural innovation and experimentation at Ganj Dareh, Iran. World Archaeol. 1990;21(3): 323–335.
  11. 11. Meiklejohn C, Merrett DC, Reich D, Pinhasi R. Direct dating of human skeletal material from Ganj Dareh, Early Neolithic of the Iranian Zagros. J Archaeol Sci Rep. 2017;12: 165–172.
  12. 12. Young TC Jr, Smith PEL. Research in the prehistory of central western Iran. Science. 1966;153: 386–391. pmid:17839705
  13. 13. Riel-Salvatore J, Lythe A, Uribe Albornoz A. New insights into the PPN spatial organization, stratigraphy, and human occupations of Ganj Dareh, Iran. (submitted).
  14. 14. Smith PEL. Survey of excavations in Iran during 1979–71. Iran. 1972;X: 165–169.
  15. 15. Smith PEL. Ganj Dareh Tepe. Paléorient. 1974;2: 207–209.
  16. 16. Smith PEL. Survey of excavation in Iran during 1973–74. Iran. 1978;13: 178–180.
  17. 17. van Zeist W, Smith PEL, Palfenier-Vegter RM, Suwijn M, Casparie WA. An archaeobotanical study of Ganj Dareh Tepe, Iran. Palaeohistoria. 1984;XXVI: 201–224.
  18. 18. Kislev E. Pre-domesticated cereals in the Pre-Pottery Neolithc A period. In: Hershkovitz I, editor. People and culture in change: Proceedings of the second symposium on Upper Palaeolithic, Mesolithic and Neolithic populations of Europe and the Mediterranean basin. Oxford: BAR International Series 508;1989. pp. 147–151. https://doi.org/10.1016/0022-2836(89)90238-6 pmid:2547971
  19. 19. Willcox G. Measuring grain size and identifying Near Eastern cereal domestication: evidence from the Euphrates valley. J Archaeol Sci. 2004;31: 145–150.
  20. 20. Cox G, and Sealy J. Investigating identity and life histories: Isotopic analysis and historical documentation of slave skeletons found on the Cape Town foreshore, South Africa. Int J Hist Archaeol. 1997;1(3): 207–224.
  21. 21. Katzenberg MA. Stable isotope analysis: A tool of studying past diet, demography, and life history. In: Katzenberg M, Saunders S, editors. Biological anthropology of the human skeleton. Hoboken, NJ: John Wiley & Sons, Inc.; 2008. pp. 413–441.
  22. 22. Müldner G, Richards MP. Stable isotope evidence for 1500 years of human diet at the city of York, UK. Am J Phys Anthropol. 2007;133(1): 682–697. pmid:17295296
  23. 23. Richards MP, Hedges R, Molleson T, Vogel J. Stable isotope analysis reveals variations in human diet at the Poundbury Camp Cemetery site. J Archaeol Sci. 1998;25(12): 1247–1252.
  24. 24. Nehlich O. The application of sulphur isotope analyses in archaeological research: A review. Earth Sci Rev. 2015;142: 1–17.
  25. 25. Boutton TW. Stable carbon isotope ratios of natural materials: II. Atmospheric, terrestrial, marine, and freshwater environments. In: Coleman D., editor. Carbon isotope techniques. San Diego: Academic Press: 1991. pp. 173–186.
  26. 26. Chisholm BS, DE Nelson, Schwarcz HP. Stable-carbon isotope ratios as a measure of marine versus terrestrial protein in ancient diets. Science. 1982;216(4550): 1131–1132. pmid:17808502
  27. 27. van der Merwe NJ. Carbon isotopes, photosynthesis, and archaeology: Different pathways of photosynthesis cause characteristic changes in carbon isotope ratios that make possible the study of prehistoric human diets. Am Sci. 1982;70(6): 596–606.
  28. 28. Hedges REM, Reynard LM. Nitrogen isotopes and the trophic level of humans in archaeology. J Archaeol Sci. 2007;34(8): 1240–1251.
  29. 29. Fuller BT, Fuller JL, Harris DA, Hedges REM. Detection of breastfeeding and weaning in modern human infants with carbon and nitrogen stable isotope ratios. Am J Phys Anthropol. 2006;129: 279–293. pmid:16261548
  30. 30. Mekota A, Grupe G, Ufer S, Cuntz U. Serial analysis of stable nitrogen and carbon isotopes in hair: Monitoring starvation and recovery phases of patients suffering from anorexia nervosa. Rapid Commun Mass Spectrom. 2006;20: 1604–1610. pmid:16628564
  31. 31. Fraser RA, Bogaard A, Heaton T, Charles M, Jones G, Christensen Bent T, et al. Manuring and stable nitrogen isotope ratios in cereals and pulses: Towards a new archaeobotanical approach to the inference of land use and dietary practices. J Archaeol Sci. 2011;38(10): 2790–2804.
  32. 32. Szpak P, Millaire J-F, White C, Longstaffe F. Influence of seabird guano and camelid dung fertilization on the nitrogen isotopic composition of field-grown maize (Zea Mays). J Archaeol Sci. 2012;39(12): 3721–3740.
  33. 33. Sealy J, van der Merwe N, Thorp J, Lanham J. Nitrogen isotopic ecology in southern Africa: Implications for environmental and dietary tracing. Geochim Cosmochim Acta. 1987;51(10): 2707–2717.
  34. 34. Thompson AH, Richards M, Shortland A, Zakrzewski S. Isotopic palaeodiet studies of ancient Egyptian fauna and humans. J Archaeol Sci. 2005;32(3): 451–463.
  35. 35. Zazzo A, Monahan FJ, Moloney AP, Green S, Schmidt O. Sulphur isotopes in animal hair track distance to sea. Rapid Commun Mass Spectrom. 2011;25(17): 2371–2378. pmid:21818798
  36. 36. Merrett DC. Bioarchaeology in Early Neolithic Iran: Assessment of health status and subsistence strategy. PhD dissertation, University of Manitoba. 2004.
  37. 37. Merrett DC. Infant health: Assessing adaptation and stability of transitional Early Neolithic subsistence economies in the Zagros Mountains. 75th Annual Meeting of American Association of Physical Anthropologists, Anchorage, Alaska. Am J Phys Anthropol. 2006;S41: 130.
  38. 38. Lazaridis I, Nadel D, Rollefson G, Merrett DC, Rohland N, Mallick S, et al. Genomic insights into the origin of farming in the ancient Near East. Nature. 2016;536: 419–424. pmid:27459054
  39. 39. Narasimhan VM, Patterson N, Moorjani P, Rohland N, Bernardos R, Mallick S, et al. The formation of human populations in South and Central Asia. Science. 2019;365(6457): eaat7487, 15pp. pmid:31488661
  40. 40. Meiklejohn C, Agelarakis A, Akkermans PA, Smith PEL, Solecki R. Artificial cranial deformation in the Proto-Neolithic and Neolithic Near East and its possible origin: Evidence from four sites. Paléorient. 1992;18(2): 83–97.
  41. 41. Moorees CFA, Fanning EA, Hunt EE. Formation and resorption of three deciduous teeth in children. Am J Phys Anthropol. 1963a;21: 205–213. pmid:14110696
  42. 42. Moorees CFA, Fanning EA, Hunt EE. Age formation by stages for ten permanent teeth. J Dent Res. 1963b;42: 1490–1502.
  43. 43. Smith P, Avishai G. The use of dental criteria for estimating postnatal survival in skeletal remains of infants. J Archaeol Sci. 2005;32: 83–89.
  44. 44. Scheur L, Black S. Developmental juvenile osteology. London: Academic Press; 2000.
  45. 45. Fazakas IG, Kósa F. Forensic fetal osteology. Budapest: Akadémiae Kiadó; 1978.
  46. 46. Gindhart PS. Growth standards for the tibia and radius in children aged one month through eighteen years. Am J Phys Anthropol. 1973;39: 41–48. pmid:4351576
  47. 47. Maresh MM. Measurements from roentgenograms. In: McCammon RW, editor. Human growth and development. Springfield, IL: CC Thomas; 1970. pp.157–200. pmid:5508382
  48. 48. Merchant VL, Ubelaker DH. Skeletal growth of the protohistoric Arikara. Am J Phys Anthropol. 1977;46: 61–72. pmid:319684
  49. 49. Sundick RI. Human skeletal growth and age determination. Homo. 1978;29: 228–249.
  50. 50. Buikstra JE, Ubelaker DH. Standards for data collection from human skeletal remains. Fayetteville, AK: Arkansas Archeological Survey Research Series No. 44; 1994.
  51. 51. McKern T, Stewart TD. Skeletal changes in young American males: Analyzed from the standpoint of identification. Technical Report EP-45. Natick (MA): Headquarters, Quartermaster Research and Development Command; 1957. https://doi.org/10.1002/ajpa.1330150318 pmid:13508848
  52. 52. Scott EC. Dental wear scoring technique. Am J Phys Anthropol. 1979;51: 213–218.
  53. 53. Smith BH. Patterns of molar wear in hunter-gatherers and agriculturalists. Am J Phys Anthropol. 1984;63: 39–56. pmid:6422767
  54. 54. Brooks ST, Suchey JM. Skeletal age determination based on the os pubis: A comparison of the Acsádi-Nemeskéri and Suchey-Brooks methods. Hum Evol. 1990;5: 227–238.
  55. 55. Lovejoy CO, Meindl RS, Pryzbeck TR, Mensforth RP. Chronological metamorphosis of the auricular surface of the ilium: A new method for the determination of age at death. Am J Phys Anthropol. 1985;68: 47–56. pmid:4061601
  56. 56. Mann RW, Symes SA, Bass WM. 1987. Maxillary suture obliteration: Aging the human skeleton based on intact or fragmentary maxilla. J Forensic Sci. 1987;32: 148–157. pmid:3819673
  57. 57. Meindl RS, Lovejoy CO. Ectocranial suture closure: A revised method for the determination of skeletal age at death based on the lateral-anterior sutures. Am J Phys Anthropol. 1985;68: 57–66. pmid:4061602
  58. 58. Acsádi GY, Nemeskéri J. History of human life span and mortality. Budapest: Akadémai Kiadó; 1970.
  59. 59. Phenice TW. Newly developed visual method of sexing the os pubis. Am J Phys Anthropol. 1969;30: 297–302. pmid:5772048
  60. 60. Walker PL. Greater sciatic notch morphology: Sex, age, and population differences. Am J Phys Anthropol. 2005;127: 385–391. pmid:15693026
  61. 61. Longin R. New method of collagen extraction for radiocarbon dating. Nature. 1971;230: 241–242. pmid:4926713
  62. 62. Ambrose SH. Preparation and characterization of bone and tooth collagen for isotopic analysis. J Archaeol Sci. 1990;17(4): 431–451.
  63. 63. Bocherens H, Drucker D, Taubald H. Preservation of bone collagen sulphur isotopic compositions in an early Holocene river-bank. Archaeological Site. Palaeogeogr Palaeoclimatol Palaeoecol. 2011;310: 32–38.
  64. 64. DeNiro MJ. Postmortem preservation and alteration of in vivo bone collagen isotope ratios in relation to palaeodietary reconstruction. Nature. 1985;317: 806–809.
  65. 65. Harbeck M, Grupe G. Experimental chemical degradation compared to natural diagenetic alteration of collagen: Implications for collagen quality indicators for stable isotope analysis. Archaeol Anthropol Sci. 2009;1(1): 43–57.
  66. 66. Nehlich O, Richards MP. Establishing collagen quality criteria for sulphur Isotope analysis of archaeological bone collagen. Archaeol Anthropol Sci. 2009;1(1): 59–75.
  67. 67. van Klinken GJ. Bone collagen quality indicators for palaeodietary and radiocarbon measurements. J Archaeol Sci. 1999;26(6): 687–695.
  68. 68. Broushaki F, Thomas MG, Link V, López S, van Dorp L, Kirsanow K, et al. 2016. Early Neolithic genomes from the eastern Fertile Crescent. Science Jul 29 2016;353(6298): 499–503. pmid:27417496
  69. 69. Smith BD. Low-level food production. J Archaeol Res. 2001;9(1): 1–43.
  70. 70. Gallego-Llorente M, Connell S, Jones ER, Merrett DC, Jeon Y, Eriksson A, et al. The genetics of an early Neolithic pastoralist from the Zagros, Iran. Nature Scientific Reports 2016;6: 31326. pmid:27502179
  71. 71. Müldner G. Isotope analysis of animal bone. In: Matthews R, Matthews W, Mohammadifar Y, editors. The earliest Neolithic of Iran, 2008 excavations at Sheikh-E Abad and Jani. Oxford: Oxbow Books; 2013. pp. 321–325. https://doi.org/10.1002/ajpa.22420 pmid:24226894
  72. 72. Pearson JA, Bogaard A, Charles M, Hillson SW, Larsen CS, Russell N, et al. Stable carbon and nitrogen isotope analysis at Neolithic Çatalhöyük: Evidence for human and animal diet and their relationship to households. J Archaeol Sci. 2015;57: 69–79.
  73. 73. Dupras TL, Tocheri MW. Reconstructing infant weaning histories at Roman period Kellis, Egypt using stable isotope analysis of dentition. Am J Phys Anthropol. 2007;115: 204–212. pmid:17568441
  74. 74. Fogel ML, Tuross N, Johnson BJ, Miller GH. Biogeochemical record of ancient humans. Org Geochem. 1997;27(5–6): 275–287.
  75. 75. Herrscher E, Goude G, Letz L. Longitudinal study of stable isotope compositions of maternal milk and implications for the palaeodiet of infants. Bull Mem Soc Anthropol Paris. 2017;29(3–4): 31–139.
  76. 76. Perkins MJ, McDonald RA, van Veen FJF, Kelly SD, Rees G, Bearhop S. Application of nitrogen and carbon stable isotopes (δ15N and δ13C) to quantify food chain length and trophic structure. PLoS One. 2014;9(3) e93281. pmid:24676331
  77. 77. Stantis C, Schutkowski H, Sołtysiak A. Reconstructing breastfeeding and weaning practices in the Bronze Age Near East using stable nitrogen isotopes. Am J Phys Anthropol. 2019;172(1): 58–69. pmid:31797366
  78. 78. Burt NM, Garvie-Lok S. A new method of dentine microsampling of deciduous teeth for stable isotope analysis. J Archaeol Sci. 2013;40: 3854–3864.
  79. 79. D’Ortenzio L, Brickley M, Schwarcz H, Prowse T. You are not what you eat during physiological stress: Isotopic evaluation of human hair. Am J Phys Anthropol. 2015;157: 374–388. pmid:25711625
  80. 80. Schurr MR. Exploring ideas about isotopic variation in breastfeeding and weaning within and between populations: Case studies from the American midcontinent. Int J Osteoarchaeol. 2018;2018: 1–13.