Scientific Reports

Sci Rep. 7: 7146

Earliest expansion of animal husbandry beyond the Mediterranean zone in the sixth millennium BC

Jonathan Ethier¹, Eszter Bánffy², Jasna Vuković³, Krassimir Leshtakov⁴, Krum Bacvarov⁵, Mélanie Roffet-Salque⁶, Richard P. Evershed⁶, Maria Ivanova

ivanova@uni-heidelberg.de

¹

1. 0000 0001 2190 4373grid.7700.0Institut für Ur- und Frühgeschichte und Vorderasiatische Archäologie, Universität Heidelberg, Marstallhof 4, 69117 Heidelberg, Germany

2. Römisch-Germanische Kommission des Deutschen Archäologischen Instituts, Palmengartenstr. 10-12, 60325 Frankfurt, Germany

3. 0000 0001 2166 9385grid.7149.bDepartment of Archaeology, Faculty of Philosophy, University of Belgrade, Čika Ljubina 18-20, 11000 Belgrade, Serbia

4. 0000 0001 2192 3275grid.11355.33Department of Archaeology, Faculty of History, Sofia University St. Kliment Ohridski, 15 Tzar Osvoboditel Boulevard, 1504 Sofia, Bulgaria

5. 0000 0001 2097 3094grid.410344.6National Institute of Archaeology and Museum, Bulgarian Academy of Sciences, 2 Saborna St., 1000 Sofia, Bulgaria

6. 0000 0004 1936 7603grid.5337.2Organic Geochemistry Unit, School of Chemistry, University of Bristol, Cantock’s Close, Bristol, BS8 1TS United Kingdom

© The Author(s) 2017

DOI: 10.1038/s41598-017-07427-x

Published online: 2 August 2017

Abstract

Since their domestication in the Mediterranean zone of Southwest Asia in the eighth millennium BC, sheep, goats, pigs and cattle have been remarkably successful in colonizing a broad variety of environments. The initial steps in this process can be traced back to the dispersal of farming groups into the interior of the Balkans in the early sixth millennium BC, who were the first to introduce Mediterranean livestock beyond its natural climatic range. Here, we combine analysis of biomolecular and isotopic compositions of lipids preserved in prehistoric pottery with faunal analyses of taxonomic composition from the earliest farming sites in southeast Europe to reconstruct this pivotal event in the early history of animal husbandry. We observe a marked divergence between the (sub)Mediterranean and temperate regions of Southeast Europe, and in particular a significant increase of dairying in the biochemical record coupled with a shift to cattle and wild fauna at most sites north of the Balkan mountain range. The findings strongly suggest that dairying was crucial for the expansion of the earliest farming system beyond its native bioclimatic zone.

Introduction

The wild progenitors of the main domestic animals in the Old World are endemic to regions with Mediterranean climate and are adapted to withstand prolonged hot summer droughts and mild but wet winters¹, ². Since domestication, farmers have brought sheep, goats, pigs and cattle to an enormous variety of environments, from semi-deserts to sub-arctic regions. Their present-day distribution, pushed out to the boundaries of the world inhabitable by humans, was mediated through human protection and breeding of animals that thrive under conditions often not tolerated by their wild ancestors. The first steps in this process can be traced back to the farming pioneers of the Balkans who penetrated beyond the borders of the sub-Mediterranean zone of Europe in the early centuries of the sixth millennium BC. How the first Balkan herders succeeded in extending the habitat of their livestock is an intriguing question with pivotal importance for the early history of human-animal relationships.

During the later seventh and early sixth millennia BC permanent farming settlements, similar to contemporary sites in the core areas of domestication in southwest Asia, became established in the (sub-)Mediterranean southern Balkans³, ⁴. The spread of farming economy into the temperate northern parts of the peninsula, however, was accompanied by pronounced changes, including higher (probably seasonal) residential mobility, smaller community sizes and a loss of sophistication in architecture and material culture⁵–⁸, a phenomenon which has been designated in the archaeological literature as “the First Temperate Neolithic”. Although the adaptation of herding economy to new bioclimatic conditions has been recognised as a major component of this phenomenon⁸–¹⁴, the human strategies which promoted it have remained uncertain.

Since early farmers kept domestic animals primarily for food and were dependent on them for survival, we can assume that their use of animal foodstuffs was finely adjusted to the animals’ productivity under different ecological conditions. Differences in the use of meat and milk at different latitudes in the Balkans are therefore likely to reflect regional adaptations in the course of the northward expansion of herding. Herein, we combine organic residue analysis of ceramic vessels with faunal data on taxonomic composition to investigate variation in the use of animal products across the different bioclimatic zones of the Balkans.

Results

Faunal remains are widely used in archaeology as a proxy for past interactions between people and animals¹⁵–¹⁷. Datasets with taxonomic frequencies are available from the majority of recently excavated early farming sites in Southeast Europe. These datasets are generally dominated by food species and provide information on the relative proportions of consumed animals, but do not directly reflect differences in the used products (milk, meat, and fat). Investigations of fat residues in archaeological pottery can directly verify the use of particular animal foodstuffs, and differentiate between dairy, ruminant carcass and non-ruminant carcass fats, as well as a range of other commodities, such as aquatic resources or beeswax¹⁸–²². Lipid residues and faunal remains are thus complementary strands of evidence for recovering past practices of animal exploitation²³.

Regional patterns in lipid residues

From the 290 pottery samples analysed in this study, 22% (n = 64) provided identifiable lipid residues with concentrations up to 5679 μg g⁻¹ (mean 802 μg g⁻¹, for the sherds containing lipids). This preservation rate is comparable with previous studies of prehistoric pottery from the Mediterranean²⁴, but lower than the proportions of over 50% of sherds with extracted lipids from Central Europe²⁵. All sites were nearly equally represented with ca. 20% to 35% of the samples containing identifiable lipids of archaeological origin, the only exception being Divostin, where only 7% of sherds contained observable lipids (Table 1). The lipid extracts were dominated by fatty acids derived from animal fat, i.e. palmitic acid (n-hexadecanoic acid; C16:0) and stearic acid (n-octadecanoic acid; C18:0; Fig. 1, Supplementary Table S1). Biomarkers characteristic of plant material were not detected.

See full tableTable 1. Sample information and summary of organic residue analysis of pottery

Figure 1View larger versionFigure 1. Gas chromatograms of TLEs from pottery from the three regions of study showing high concentration of palmitic acid (C16:0) and stearic acid (C18:0): (a) Ecsegfalva 23, sample 16, (b) Grivac, sample 17, and (c) Yabalkovo, sample 38. IS – internal standard: C₃₄ n-tetratriacontane.

The identification of the various commodities was based on comparison of the carbon isotope compositions of archaeological fatty acids to those of modern reference animal fats. To eliminate differences in diets between the modern reference and the prehistoric animals caused by environmental and dietary factors, the Δ¹³C values (δ¹³C18:0 − δ¹³C16:0) of the archaeological samples were calculated and compared to that of modern reference fats²⁵–²⁹. All archaeological lipids derived from ruminant or non-ruminant animals. Potential contributors of non-ruminant fats include domestic pig, wild boar, fowl and domestic dog. No biomarkers of aquatic resources, such as ω-(o-alkylphenyl)alkanoic acids, were detected.

The difference between the southern and northern part of the Balkans in respect to the distribution of commodities among the lipid residues is highly significant (p = .002, Fisher’s exact test). Over 70% (n = 17) of all animal fat residues detected in pottery from the southern part of the Balkans in the present study plot within the range of ruminant adipose fats (Fig. 2b). At Yabalkovo, all three sources of animal fat were identified, although non-ruminant and dairy residues were represented only in low proportions (20% each). The pattern is even more pronounced in the distribution of residues from Nova Nadezhda, where dairy fats were not detected and the proportion of ruminant adipose fats reached 89%. These findings agree well with previously published data from early farming sites in other areas with Mediterranean climate, e.g. the Levant, Syria, Anatolia, Greece and Italy²⁴, ²⁷, ³⁰, ³¹. We can denote this type of early livestock use, reflected in strong predominance of ruminant carcass fat and low incidence or absence of dairy fat among the residues in ceramic vessels, the “Mediterranean type”.

Figure 2View larger versionFigure 2. Regional differences in Δ¹³C values ( = δ¹³C18:0 − δ¹³C16:0) of archaeological lipid residues in pottery. (a) proportions of lipid residues according to commodity type; (b) Δ¹³C values of archaeological lipid residues from sites in (b) the southern Balkans, (c) northern Balkans (d) the Carpathian basin (including published data for Ecsegfalva³²). Each data point represents an individual vessel.

The lipid residues from sites in the northern part of the Balkans and in the Hungarian Plain, characterized by transitional and sub-continental climate, demonstrate a different pattern (Fig. 2c,d). Notable similarities exist between distant sites in present-day Serbia and Hungary (with the only exception of Divostin, which is not representative because of the very small number of sherds with residues, n = 3). Residues classified as non-ruminant fats are nearly absent, which agrees with the very low percentages of pigs in the faunal assemblages. With 60% (n = 24), dairy fats are by far the most abundant class of animal fat residues. The results obtained in the current study represent the first substantial dataset from the earliest farming sites in the temperate zone. They are comparable to the results of a previous study of pottery from Ecsegfalva 23 with n = 7 sherds with residues³². The type of animal product exploitation reflected in a residue record with a strong emphasis on dairy fats and nearly absent non-ruminant residues can be provisionally assigned as the “Temperate type”.

Faunal data on species preferences

The taxonomic compositions of the faunal assemblages from early farming sites in the interior of the southern Balkans are similar to these from the Mediterranean littoral of the peninsula, with a somewhat lower proportion of caprines and more cattle compared to the Aegean (Fig. 3). The faunal data from the northern Balkans, in contrast, suggests a clearly different strategy of animal foodstuff supply. The proportion of caprines decreased compared to the south and cattle assumed a major role among livestock in the northern Balkans, clearly predominating as a supplier of meat when their body mass is considered. Together with the increasing importance of large game, this change in taxonomic frequencies highlights a shift to ‘‘local’’ food species (animals which are native to or have close endemic relatives in temperate Europe). Surprising is the nearly complete absence of domestic pigs (represented by only a few percent of NISP), since pigs must have thrived in the temperate forests and marshlands. One possible explanation for the scarcity of pigs in the northern Balkans is seasonal mobility of herds: although pig transhumance is known from the ethnographic record³³, pigs are in general untypical for the seasonally relocated herds. It has been previously argued that the preference for mobile fauna in the northern Balkans relates to other mobility indicators, such as smaller community size, lower settlement permanence, less substantial or semi-subterranean dwellings, and possibly decreasing reliance on plant cultivation³⁴. Seasonal mobility is also tentatively suggested by the age distribution of caprines and by the absence of migratory species among birds in the faunal assemblage from Blagotin, indicating cold season (late autumn to late spring) presence at this location³⁵. However, other causes for the decline in pig husbandry also need to be considered.

Figure 3View larger versionFigure 3. General trends in the exploitation of the most abundant mammalian taxa. Faunal data from early farming sites in the Aegean, the Balkans and the Carpathian basin, symbol size corresponds to relative body size of species, each symbol represents 10% of NISP. See Supplementary Table S2 for the complete dataset.

Similar taxonomic abundances are characteristic of the faunal distributions from the Hungarian Plain. The small ruminants seem more numerous here compared to the northern Balkans, though this picture is biased by several sites in the northern periphery of the region which focused intensively on caprines (70 to 80% of NISP; Fig. 4). Since the swampy wetland environments of the Danube and Tisza plains are adverse to sheep husbandry, these assemblages have often been described as puzzling³⁶, ³⁷ and the preference for caprines has been attributed to emotional motives of the early farmers, who adhered to the pastoral traditions of their Mediterranean ancestors against all odds³⁸. However, a large area between the Aegean and the Carpathian basin where ‘‘local’’ species (cattle and large game) prevailed in the faunal remains indicates that continuity of herding practices is not a completely satisfactory explanation for the sites with high proportions of caprines in eastern Hungary (Fig. 4). Rather, the latter can be considered more anomalous than typical for the temperate zone. A likely pragmatic reason for the prevalence of caprines at some sites in Hungary, despite the unfavourable environment, are problems with livestock overwintering due to snow cover and shorter plant growing season. “Caprine” sites are situated in grassland areas, where large game hunting was not as rewarding as on the forest margins³⁹. In view of their higher reproductive rate, opting for caprines (rather than cattle) may be an advantage when herders are forced to slaughter many of their animals before or during winter because of limited fodder, since the annual growth rate in number of heads in sheep herds and goat herds is higher than that of cattle⁴⁰. Moreover, in competition with pigs, caprines possess an advantage which becomes crucial under conditions of decreasing animal productivity – milk exploitation yields up to four times higher yield of food energy compared to hunting or growing animals for meat⁴¹.

Figure 4View larger versionFigure 4. Summary of climatic conditions in Southeast Europe. Dots indicate the location of faunal assemblages included in this study (see Supplementary Table S2); black dots – sites where caprines were represented with >60% of NISP of the large- and middle-sized domestic and wild mammalian taxa, white dots – caprines <60%. The dotted lines indicate regional boundaries between the Aegean/Adriatic littoral, southern Balkans, northern Balkans and Carpathian Basin. Numbers show location of sites where organic residue analyses were performed 1. Yabalkovo, 2. Nova Nadezhda, 3. Blagotin, 4. Divostin, 5. Grivac, 6. Alsónyék, 7. Ecsegfalva 23. (Climate data from: Hijmans, R. J., S. E. Cameron, J. L. Parra, P. G. Jones and A. Jarvis. 2005. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25(15): 1965-1978; http://worldclim.org/, Layer BIO6 = Min Temperature of Coldest Month, interpolations of observed data, representative of 1960–1990). The map was created using QGIS Version 2.8.9-Wien (http://qgis.org/de/site/).

Discussion

The interior of the Balkan Peninsula and the adjacent Hungarian Plain occupy an intermediate position between the Mediterranean and the Central European bioclimatic regions. The present-day distribution of thermophilous mixed deciduous broadleaved forests indicates fluid transitions of temperature and precipitation, from mild winter temperatures in the south and occasional winter frosts in the northern parts of the peninsula to regular periods of winter frost and snow in the Hungarian plain⁴². Reconstructions of warmest month, coldest month and mean annual surface air temperatures, based on fossil pollen data, suggest that winter temperatures in southeast Europe have generally remained at the same level from 6000 calBC through to the present (Fig. 4)⁴³.

Farmers and their domestic animals spread fast, in c. 10 human generations, across these biogeographic zones from sub-Mediterranean Macedonia to the northern limits of the temperate Carpathian Basin¹⁰, ⁴⁴–⁴⁶. The increasingly unfavourable climate, and particularly the long frosty winters, had certainly an impact on Mediterranean livestock. Climate affects animals in a direct way through thermoregulation, and indirectly, through trophic links (e.g. pasture quality and biomass)⁴⁷. Climate-driven selective pressure brings forward changes in the genome, such as the adaptation to feed scarcity through accumulation of fat as energy reserve, or through a seasonal decline in appetite and the ability to reduce metabolism in winter⁴⁸, ⁴⁹. Genetically-determined improvements of feed efficiency, pigmentation, coat and body size are further common adaptations to harsh environments.

However, since genetic changes do not occur immediately, farmers spreading into higher latitudes in the Balkans are likely to have initially experienced difficulties in maintaining the productivity of all domesticates, and in particular of sheep, which do not have native wild relatives in Europe and, being grazers, are not adapted to feed in the forested environments of the Balkan interior. The most serious hurdle to overcome was the provision of feed under conditions of a shorter period of plant growth and deep snow impeding access to grass for weeks in winter. Spring temperatures have also been recognized for their effect on animal productivity. For example, 1 °C decrease in mean spring temperature has been estimated to reduce lamb autumn body mass by nearly 0.4 kg⁵⁰. One further climate-related parameter, directly associated with animal productivity, is the seasonality of mating and birth, which is determined by day-length through a genetically controlled mechanism for sheep and through feed availability for cattle⁵¹, ⁵².

Before genetically adapted breeds became established, these obstacles must have been overcome by techniques of herd management, such as control of mating, strategic culling, adjustments in the species composition, interbreeding with wild populations, and labour-intensive feeding practices (stalling, foddering, change of pastures). Equally important must have been, however, adjustments in the use of animal foodstuffs and in the techniques of food processing. For example dairying, a practice familiar to the early farmers of southwest Asia and the Mediterranean²⁴, ²⁷, ⁵³, is highly efficient in increasing the yields of animal husbandry. Milk pastoralism can provide up to four times higher yields of food energy in comparison to growing animals for meat⁴¹. T. Ingold lucidly explains the higher productivity of milking in comparison to meat exploitation by the respective position of the herders in the food chain. Pastoralist who milk their animals intercept the chain of conversions from grass to milk, and from milk to meat, at an earlier point, thus avoiding the net loss of energy that takes place at the latter stage of this conversion⁵⁴.

Our investigation of animal lipids preserved in pottery and of taxonomic frequencies of large- and medium-sized mammalian species shows differences in the use of animal products between the Aegean and the southern Balkans, on the one hand, and the northern Balkans and the Hungarian Plain on the other hand. We observe a significant increase of dairy fats in the biochemical record coupled with a decline in caprines at most sites north of the Balkan mountain range. Our study indicates that dairying, being a highly efficient way to counterbalance reduced productivity of animal husbandry, was instrumental to the spread of the Mediterranean livestock system into the temperate areas of the Balkans, and from there across the entire European continent. The general distinction of a “Mediterranean” and a “Temperate” type of early animal husbandry outlined in the present study does not rule out the existence of small-scale regional variation. Previous studies of lipids in early ceramics in the Mediterranean, for example, suggest that a high proportion of dairy residues might be characteristic for early farmers’ habitations in caves and rock-shelters²⁴, ⁵⁵, ⁵⁶, possibly pointing at the presence of specialized seasonal herding sites.

Finally, dairying represents not only an efficient economic strategy, but also provides an important dietary component. Milk was one of the main potential sources of dietary fat for the first Balkan farmers, together with lard, oily fish and probably waterfowl. Having in mind that pigs were barely represented in the faunal assemblages from the northern Balkans and the Carpathian Basin, and that non-ruminant lipids were very rare among the residues in pottery, we can regard dairying as a key strategy of obtaining fat in this geographic region. Fresh milk in small-scale non-industrial animal husbandry is a seasonal product. Peak lactation in temperate climates generally coincides with the warm season, when pastures can support the high nutritional requirements of lactating stock⁵¹. Being a perfect medium for pathogenic microbial growth, in the hot months milk spoils to an ill-smelling liquid within hours, if not consumed or heated and processed in a controlled manner (for example through controlled souring by lactic bacteria). The low frequency or absence of the −13,910*T allele, associated with continued production of lactase into adulthood (and thus with the ability to digest fresh milk), in hunter-gatherer and early farming populations in Europe⁵⁷–⁶¹, suggests that high-lactose fresh milk could not have been consumed in large quantities by the first farmers. Hence, most of the milk must have been processed without delay both to reduce its lactose content and to convert it into a storable year-round staple food. The increase in dairying observed in this study thus represents a remarkable example of complex cultural behavior employed to circumvent the limitations of environment and of human and animal biology.

Methods

Organic residue analysis

Sampling

In this study, we analysed a total of 290 potsherds from seven pottery assemblages dating to c. 6100–5700 calBC (Table 1). The sherds derive from well-documented and published excavations in Bulgaria (Yabalkovo and Nova Nadezhda), Serbia (Blagotin, Divostin and Grivac), and Hungary (Alsónyék and Ecsegfalva 23), which have provided the earliest Neolithic assemblages in their respective geographic area³⁵, ⁶²–⁶⁷ (Supplementary materials). Sherds were selected following two main criteria – an undisturbed stratigraphic context and a representative spectrum of wares and shapes (including diverse types of plates, bowls and jars); they were mostly diagnostic rim fragments (73%).

Lipid extraction and derivatisation

Lipid extractions were performed using direct acidified methanol extraction, according to the protocol of Correa-Ascencio and Evershed⁶⁸. Using a modelling drill, a small area of each fragment was cleaned to remove modern contaminants (from soil burial, excavation and post-excavation handling) and a sample of ca. 2 g was taken. The sherd fragment was crushed into a fine powder, weighed and an internal standard of n-tetratriacontane (20 μl, 1.0 mg ml⁻¹ solution) was added. Lipids were extracted from the powder by adding 5 ml of a H₂SO₄-MeOH (4% sulphuric acid/methanol) solution (δ¹³C known), heated for 1 h at 70 °C, whirl mixed every 10 min to ensure contact between solvent and powdered sherd. In order to help observe the solvent layers, the H₂SO₄-MeOH was transferred to a clean culture tube, where 2 ml of DCM extracted-double-distilled water was added. To ensure the recovery of lipids not fully solubilized by the methanol solution, hexane (2 ml portions) washing was performed 2x for the powdered sherd and 3x for the diluted methanol extract. All hexane extracts were combined and evaporated to dryness under a gentle stream of nitrogen. Prior to instrumental analysis one fourth aliquots of total lipid extracts (TLEs) diluted in hexane were blown down to dryness with nitrogen and derivatised with 20 μl of N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA containing 1% trimethylsilyl chloride) for 1 h at 70 °C to trimethylsilylate free alcohols.

Gas chromatography (GC), GC/mass spectrometry (GC/MS) and GC-combustion-isotope ratio MS (GC-C-IRMS)

Derivatised TLEs were blown to dryness then diluted with hexane and screened by GC to make preliminary identifications of components and quantify lipid concentrations in potsherds by reference to the n-tetratriacontane internal standard. Two different GC systems were used, employing either manual or autosampling injections: (1) Agilent Technologies 7820 series 2 with a silica capillary column (50 m × 0.32 mm) coated with a dimethyl polysiloxane stationary phase (J&W Scientific; CP-Sil 5 CB, 0.12 μm film thickness) or (2) Agilent 7890a with a capillary column (15 m × 0.32 mm) coated with a dimethyl polysiloxane stationary phase (J&W Scientific; DB1-HT, 0.1 μm film thickness). On both systems, the carrier gas was helium and a flame ionization detector (FID) was used to monitor the column effluent. The temperature programs began with an initial hold at 50 °C for 2 min followed by temperature programming at 10 °C min⁻¹ and a final hold for 10 min at 350 °C (Agilent 7890a) or 300 °C (Agilent 7820). Collection of the data and quantification were performed using a Agilent Chemstation software.

Where structures of components were not identifiable based on GC alone, GC/MS analyses were performed on a ThermoFinnigan Trace MS quadrupole mass spectrometer coupled to a Trace GC. Using either manual or auto-injections samples were injected via the splitless mode onto a fused-silica capillary column (60 m × 0.32 mm) coated with a polymethylpolysiloxane stationary phase (Phenomenex; ZB1, 0.1 μm film thickness). The temperature programs for the GC began with an initial temperature held at 50 °C for 2 min and followed by a constant increase at a rate of 10 °C min⁻¹ to 300 °C then held for 10 min. The MS was operated in electron ionisation (EI) mode (70 eV) with a GC interface temperature of 300 °C and a source temperature of 200 °C. The emission current was 150 μA and the data collection was between m/z 50–650 at 1.3 scans s⁻¹. Xcalibur v3.0 software and a NIST spectra library were used to acquire, analyse and identify separated compounds. The same samples were screened for ω-(o- alkylphenyl)alkanoic acids using GC/MS selected ion monitoring (SIM; m/z 105, 262, 290, 318 and 346). The extracts were injected in the same GC as above equipped with a fused-silica capillary column (60 m × 0.32 mm) coated with a 50% cyanopropyl-methylpolysiloxane stationary phase (Varian, Factor Four; VF-23 ms, 0.15 μm film thickness). The GC oven temperature was programmed from 50 °C, following an isothermal hold for 2 min, to 100 °C at 10 °C min⁻¹, and then to 240 °C at 4 °C min⁻¹, following by an isothermal hold at 240 °C for 15 min. The extracts dominated by C16:0 and C18:0 fatty acids and free from plant lipids (e.g. long-chain n-alkanols and long-chain n-alkanoic acids), were analysed in duplicate by GC-C-IRMS to obtain their individual δ¹³C values. GC-C-IRMS analyses used an Agilent 6890 GC coupled to a Finnegan MAT Deltaplus XL mass spectrometer via a Finnigan MAT GCCIII interface. Samples were injected using a CTC A200S autosampler in splitless mode onto a Varian silica capillary column (CP-Sil 5 CB, 100% dimethylpolysiloxane with 0.12 μm film thickness, 50 m × 0.32 mm). The temperature program for the GC began with an initial temperature held at 50 °C for 2 min, followed by a constant increase in oven temperature at10 °C min⁻¹ to 250 °C, followed by a 15 min hold.

Faunal analyses

Taxonomic abundances in 39 early farming faunal assemblages from the Aegean, the Balkans and the Carpathian Basin were calculated based on NISP (number of identifiable specimens, identified to species level, with the exception of indistinguishable sheep/goat remains that represent the caprine subfamily; Supplementary Table S2). To limit bias through taphonomic history, recovery techniques and laboratory procedures, only the most ubiquitous large- and medium-sized mammalian taxa were considered. Because small assemblages of less than 500 NISP tend not to represent reliably the relative importance of animals³⁸, ³⁹, the assemblages included in this study exceed 890 NISP (with the exceptions of Nea Nikomedia n = 450 and Alsónyék n = 428).

Animal bones from excavations represent butchering and cooking remains and are thus directly related to the consumption of animal products. However, straightforward estimates of dietary contribution from taxonomic abundances based on NISP are not possible, since NISP neither corresponds to meat quantity, nor to the use of animals for meat and milk products. Other proxies, such as age-at-death (AtD) data, are available only for a handful of sites (Prodromos, Blagotin, Ecsegfalva 23, Miercurea Sibiului-Petriş, Foeni Salaș and Foeni Gaz)³⁷, ⁶⁹–⁷¹. The AtD data were considered in the interpretation of the organic residues in pottery; they are, however, insufficient for inter-regional comparisons. Although no calculation of dietary contributions was attempted in this study, we take into consideration the relative contributions of large- and medium-sized livestock and large game in the discussion of their importance in diet. We adopted a modified version of Clason’s “live-weight factor” approach to express the live weight of taxa in terms of sheep equivalents (sheep = 1, cattle = 28, pig = 1.2)⁷², ⁷³. For large game, we use a factor of 22, based on the average of the aurochs (36), equid (24) and cervid (7) factors. These factors are only rough approximations of the relative meat contribution of taxa and do not represent actual body sizes or meat yields.

Supplementary Material

Open In Web Browser

Footnotes

Electronic supplementary material

Supplementary information accompanies this paper at doi:10.1038/s41598-017-07427-x

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

The organic residue analysis of pottery was conducted as part of the collaborative project ‘‘Food cultures: Interdisciplinary studies of early farming, food technology and palaeodiet in Southeastern Europe’’, funded by an award of the German Research Council, DFG-IV101/5-1 (to M.I.). Branka Zorbic (Kragujevac), Tibor Marton and Krisztián Oross (Budapest), Nadezhda Todorova and Vanya Petrova (Sofia) are thanked for their support in selecting the pottery samples from Divostin, Grivac, Alsónyék, Ecsegfalva 23, Nova Nadezhda and Yabalkovo. We are grateful to Haskel Greenfield, Laszlo Bartosiewicz and Marie Balasse for their comments on this work. We also thank Elena Marinova-Wolff for advice on Fig. 4.

Author Contributions

J.E., M.I. and R.P.E. conceived and designed the study. M.I. and J.E. wrote the paper. J.E. collected the faunal data. J.E. and M.R.-S. performed analytical work on pottery and data analysis. M.I., E.B., J.V., K.L. and K.B. either conducted sampling of pottery for lipid residue analyses or directed excavations. All authors reviewed and approved the paper.

Notes

Competing Interests

The authors declare that they have no competing interests.

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group

PMC Copyright Notice

The articles available from the PMC site are protected by copyright, even though access is free. Copyright is held by the respective authors or publishers who provide these articles to PMC. Users of PMC are responsible for complying with the terms and conditions defined by the copyright holder.

Users should assume that standard copyright protection applies to articles in PMC, unless an article contains an explicit license statement that gives a user additional reuse or redistribution rights. PMC does not allow automated/bulk downloading of articles that have standard copyright protection.

See the copyright notice on the PMC site, https://www.ncbi.nlm.nih.gov/pmc/about/copyright/, for further details and specific exceptions.

References

Vigne JD, author. The origins of animal domestication and husbandry: A major change in the history of humanity and the biosphere. Comptes Rendus Biologies. 2011;334:171–181. [Cross Ref] [PubMed]

Zeder MA, author. Domestication and early agriculture in the Mediterranean Basin: origins, diffusion, and impact. PNAS. 2008;105:11597–11604. [Cross Ref] [PubMed]

Perles, C. The Early Neolithic in Greece. The First Farming Communities in Europe (Cambridge University Press, 2001).

Nikolov, V. In Neolithic of Southeastern Europe and its Near Eastern connections (Varia Archaeologica Hungarica 2) (eds S. Bökönyi & P. Raczky) 191–200 (Institute of Archaeology of the Hungarian Academy of Sciences, 1989).

Nandris, J. In Problems in economic and social archaeology (eds G. Sieveking, I. H. Longworth, & K. Wilson) 549–556 (Duckworth, 1976).

Whittle, A. W. R. Europe in the Neolithic: the creation of new worlds (Cambridge University Press, 1996).

Tringham, R. In Europe’s First Farmers (eds T. D. Price) 19–56 (Cambridge University Press, 2000).

Nandris, J. In A Short Walk through the Balkans: The First Farmers of the Carpathian Basin and Adjacent Regions (eds M. Spataro & P. Biagi) 11–23 (2007).

Halstead, P. In The beginnings of agriculture (eds A. Milles, D. Williams, & N. Gardner) 23–53 (British Archaeological Reports, 1989).

Krauß, R., Marinova, E., Brue, H. D. & Weninger, B. The rapid spread of early farming from the Aegean into the Balkans via the Sub-Mediterranean-Aegean Vegetation Zone. Quaternary International, doi:10.1016/j.quaint.2017.01.019 (2017).

Manning, K. et al. In The origins and spread of domestic animals in Southwest Asia and Europe (eds S. Colledge et al.) 237–252 (Left Coast, 2013).

McClure SB, author. Domesticated animals and biodiversity: early agriculture at the gates of Europe and long-term ecological consequences. Anthropocene. 2013;4:57–68. [Cross Ref]

McClure SB, author. The pastoral effect: niche construction, domesticated animals and the spread of farming in Europe. Current Anthropology. 2015;56:901–910. [Cross Ref]

Orton D, author. Herding, settlement, and chronology in the Balkan Neolithic. European Journal of Archaeology. 2012;15:5–40. [Cross Ref]

Arbuckle BS, et al., authors. Data sharing reveals complexity in the westward spread of domestic animals across Neolithic Turkey. PloS one. 2014;9:e99845 [Cross Ref] [PubMed]

Colledge, S., Conolly, J., Dobney, K., Manning, K. & Shennan, S. The origins and spread of domestic animals in Southwest Asia and Europe (Left Coast Press, 2013).

Russell, N. Social zooarchaeology. Humans and animals in prehistory (Cambridge University Press, 2012).

Evershed RP, author. Organic residue analysis in archaeology: the archaeological biomarker revolution. Archaeometry. 2008;50:895–924. [Cross Ref]

Dudd SN, Evershed RP, authors. Direct demonstration of milk as an element of archaeological economies. Science. 1998;282:1478–1481. [Cross Ref] [PubMed]

Evershed RP, et al., authors. Lipids as carriers of anthropogenic signals from prehistory. Philosophical Transactions of the Royal Society. 1999;354:19–31. [Cross Ref]

Evershed RP, et al., authors. Chemistry of archaeological animal fats. Accounts of Chemical Research. 2001;35:660–668. [Cross Ref]

Evershed, R. P. In Organic Mass Spectrometry in Art and Archaeology (eds M. Perla Columnbini & F. Modugono) 391–432 (Chichester, 2009).

Roffet-Salque, M. et al. From the inside out: Upscaling organic residue analyses of archaeological ceramics. Journal of Archaeological Science: Reports (2016).

Debono Spiteri, C. et al. Regional asynchronicity in dairy production and processing in early farming communities of the northern Mediterranean. PNAS 13594–13599 (2016).

Salque M, et al., authors. Earliest evidence for cheese making in the sixth millennium B.C. in northern Europe. Nature. 2013;493:522–525. [Cross Ref] [PubMed]

Copley MS, et al., authors. Dairying in Antiquity: III - Evidence from absorbed lipid residues dating to the British Neolithic. Journal of Archaeological Science. 2005;32:523–546. [Cross Ref]

Evershed RP, et al., authors. Earliest date for milk use in the Near East and southeastern Europe linked to cattle herding. Nature. 2008;455:528–531. [Cross Ref] [PubMed]

Mukherjee, A. J., Copley, M. S., Berstan, R., Clark, K. A. & Evershed, R. P. In The Zooarchaeology of Fats, Oils, Milk and Dairying - 9th ICAZ Conference (eds J. Mulville & A. K. Outram) 77–93 (Oxbow Books, 2005).

Dunne J, et al., authors. First dairying in green Saharan Africa in the fifth millennium BC. Nature. 2012;486:390–394. [Cross Ref] [PubMed]

Nieuwenhuyse O, Roffet-Salque M, Evershed R, Akkermans PM, Russell A, authors. Tracing pottery use and the emergence of secondary product exploitation through lipid residue analysis at Late Neolithic Tell Sabi Abyad (Syria). Journal of Archaeological Science. 2015;64:54–66. [Cross Ref]

Salque M, et al., authors. New insights into the Early Neolithic economy and management of animals in Southern and Central Europe revealed using lipid residue analyses of pottery vessels. Anthropozoologica. 2012;47:46–61. [Cross Ref]

Craig OE, et al., authors. Did the first farmers of central and eastern Europe produce dairy foods? Antiquity. 2005;79:882–894. [Cross Ref]

Albarella, U., Manconi, F. & Trentacoste, A. In EthnoZooArchaeology: The Present Past of Human-Animal Relationships (ed U. Albarella) 143–159 (Oxbow, 2011).

Greenfield, H. J. & Jongsma, T. L. In Living Well Together? Settlement and Materiality in the Neolithic of South-East and Central Europe (eds D. W. Bailey, A. Whittle, & D. Hofmann) 108–130 (Oxbow Books, 2008).

Greenfield HJ, Greenfield TLJ, Jezik S, authors. Subsistence and settlement in the Early Neolithic of temperate SE Europe: a view from Blagotin, Serbia. Archaeologia Bulgarica. 2014;18:1–33

Bartosiewicz, L. In Current Research in Animal Palaeopathology: Proceedings of the Second Animal Palaeopathology Working Group Conference. BAR International Series S1844 (eds Z. Miklíková & R. Thomas) 3–13 (Archaeopress, 2008).

Bartosiewicz, L. In The Early Neolithic on the Great Hungarian Plain: investigations of the Körös culture site of Ecsegfalva 23, Co. Bekes. (=Varia Archaeologica Hungarica 21) (ed A. Whittle) 287–325 (2007).

Bartosiewicz, L. In The First Neolithic Sites in Central/South-East European Transect. Volume III. The Körös Culture in Eastern Hungary. British Archaeological Reports, IBAR International Series 2334 (eds A. Anders & Z. Siklósi) 195-204 (Archaeopress, 2012).

Bartosiewicz, L. In (Un)settling the Neolithic (eds D. Bailey & A. Whittle) 51–63 (Oxbow Books, 2005).

Dahl, G. & Hjort, A. Having Herds: Pastoral Herd Growth and Household Economy. Stockholm Studies in Social Anthropology No. 2 (University of Stockholm, 1976).

Legge, A. J. In The Zooarchaeology of Fats, Oils, Milk and Dairying (Proceedings of the 9th ICAZ conference) (eds J. Mulville & A. Outram) 8–13 (Oxbow, 2005).

Bohn, U. & Neuhäusl, R. Karte der natürlichen Vegetation Europas/Map of the Natural Vegetation of Europe. Maßstab/Scale 1: 2 500 000. (Landwirtschaftsverlag, 2014).

Davis, B. A. S., Brewer, S., Stevenson, A. C. & Guiot, J. The temperature of Europe during the Holocene reconstructed from pollen data. Quaternary Science Reviews 1701–1716 (2003).

Domboróczki, L. In Die Neolithisierung Mitteleuropas (The spread of the Neolithic to Central Europe). Internationale Tagung Mainz, 24–26. Juni 2005 (eds D. Gronenborn & J. Petrasch) 175–187 (2010).

Krauß R, et al., authors. Beginnings of the Neolithic in Southeast Europe. The Early Neolithic sequence and absolute dates from Džuljunica-Smărdeš (Bulgaria). Documenta Praehistorica. 2014;41:51–77. [Cross Ref]

Whittle A, Bartosiewicz L, Borić D, Pettit PB, Richards MP, authors. In the beginning: new radiocarbon dates for the Early Neolithic in northern Serbia and south-east Hungary. Antaeus. 2002;25:63–117

Nielsen A, Steinheim G, Mysterud A, authors. Do different sheep breeds show equal responses to climate fluctuations? Basic and Applied Ecology. 2013;14:137–145. [Cross Ref]

Atti N, Bocquier F, Khaldi G, authors. Performance of the fat-tailed barbarine sheep in its environment: adaptive capacity to alternation of underfeeding and re-feeding periods. Anim. Res. 2004;52:165–176. [Cross Ref]

Dwyer CM, author. Welfare of sheep: providing for welfare in an extensive environment. Small Rum. Res. 2009;86:14–21. [Cross Ref]

Nielsen A, et al., authors. Are responses of herbivores to environmental variability spatially consistent in alpine ecosystems? Global Change Biology. 2012;18:3050–3062. [Cross Ref] [PubMed]

Balasse M, Tresset A, authors. Environmental constraints on reproductive activity of domestic sheep and cattle: what latitude for the herder? Anthropozoologica. 2007;42:71–88

Balasse, M. et al. Sheep birth distribution in past herds: a review for prehistoric Europe (6th – 3rd millennia BC). Animal: An International Journal of Animal Bioscience (2017).

Özbal H, et al., authors. Neolitik Batı Anadolu ve Marmara Yerleşimleri Çanak Çömleklerinde Organik Kalıntı Analizleri. Arkeometri Sonucları Toplantısı. 2013;28:105–114

Ingold, T. Hunters, pastoralists and ranchers: reindeer economies and their transformations (Cambridge Univ. Press, 1980).

Šoberl L, Gašparič AZ, Budja M, Evershed RP, authors. Early herding practices revealed through organic residue analysis of pottery from the early Neolithic rock shelter of Mala Triglavca, Slovenia. Documenta Praehistorica. 2008;35:253–260. [Cross Ref]

Šoberl L, et al., authors. Neolithic and Eneolithic activities inferred from organic residue analysis of pottery from Mala Triglavca, Moverna vas and Ajdovska jama, Slovenia. Documenta Praehistorica. 2014;XLI:149–179

Burger J, Kirchner M, Bramanti B, Haak W, Thomas MG, authors. Absence of the lactase–persistence–associated allele in Early Neolithic Europeans. PNAS. 2007;104:3736–3741. [Cross Ref] [PubMed]

Malmström H, et al., authors. High frequency of lactose intolerance in a prehistoric hunter-gatherer population in northern Europe. BMC Evol Biol. 2010;10:89 [Cross Ref] [PubMed]

Lacan M, et al., authors. Ancient DNA reveals male diffusion through the Neolithic Mediterranean route. PNAS. 2011;108:9788–9791. [Cross Ref] [PubMed]

Lacan M, et al., authors. Ancient DNA suggests the leading role played by men in the Neolithic dissemination. PNAS. 2011;108:18255–18259. [Cross Ref] [PubMed]

Sverrisdóttir OO, et al., authors. Direct estimates of natural selection in Iberia indicate calcium absorption was not the only driver of lactase persistence in Europe. Mol. Biol. Evol. 2014;31:975–983. [Cross Ref] [PubMed]

Oross K, et al., authors. The early days of Neolithic Alsónyék: the Starčevo occupation. Berichte der Römisch-Germanischen Kommission. 2016;94:93–121

Roodenberg, J., Leshtakov, K. & Petrova, V. Yabalkovo, Volume 1 (Sofia University, 2014).

Whittle, A. The Early Neolithic on the Great Hungarian Plain. Investigations of the Körös Culture site of Ecsegfalva 23, County Békés. Varia Archaeologica Hungarica. (2007).

Vuković, J. In Beginnings - new research in the appearance of the Neolithic between Northwest Anatolia and the Carpathian Basin (eds R. Krauß) 205–212 (Verlag Marie Leidorf, 2011).

McPherron, A. & Srejović, D. Divostin and the Neolithic of Central Serbia (1988).

Bogdanović, M. Grivac. Naselja Protostarčevačke i Vinčanske kulture (2004).

Correa-Ascencio M, Evershed RP, authors. High throughput screening of organic residues in archaeological potsherds using direct methanolic acid extraction. Analytical Methods. 2014;6:1330–1340. [Cross Ref]

Greenfield H, Arnold E, authors. ‘Go(a)t milk?’ New perspectives on the zooarchaeological evidence for the earliest intensification of dairying in south eastern Europe. World Archaeology. 2015;47:792–818. [Cross Ref]

Halstead P, Jones G, authors. Early Neolithic economy in Thessaly: some evidence from excavations at Prodromos. Anthropologika (Athen). 1980;1:93–117

El Susi, G. In The First Neolithic Sites in Central/South- East European Transect Volume II: Early Neolithic (Starčevo-Criş) sites on the territory of Romania. British Archaeological Reports International Series 218 (eds S.A. Luca & C. Suciu) 47–56 (Archaeopress, 2011).

Clason, A. T. In Domestikationsforschung und Geschichte der Haustiere: Internationales Symposion in Budapest, 1971 (ed J. Matolcsi) 205–212 (Akadémiai Kiadó, 1973).

Russell, N. & Martin, L. In Inhabiting Çatalhöyük: Reports from the 1995-1999 Seasons McDonald Institute Monographs (eds I. Hodder) 33–98 (McDonald Institute for Archaeological Research, 2005).

[Back]

Table 1.

Sample information and summary of organic residue analysis of pottery

Sites | Dating | Potsherds analysed | Residues Detected | Non-ruminant adipose fats | Ruminant adipose fats | Ruminant dairy fats
n | % | n | % | n | % | n | %
C. Hungarian plain
Alsonyek | 5800-5730 BC | 43 | 14 | 32 | — | — | 5 | 36 | 9 | 64
Ecsegfalva 23 | 5800-5650 BC | 41 | 8 | 20 | 1 | 12 | 2 | 25 | 5 | 63
B. Northern Balkans
Blagotin | 6200-6000 BC | 43 | 9 | 21 | — | — | 2 | 22 | 7 | 78
Grivac | 6200-6000 BC | 33 | 6 | 18 | — | — | 3 | 50 | 3 | 50
Divostin | 6000-5700 BC | 43 | 3 | 7 | 2 | 67 | 0 | 0 | 1 | 33
A. Southern Balkans
Yabalkovo | 6000-5700 BC | 42 | 15 | 35 | 3 | 20 | 9 | 60 | 4 | 20
Nova Nadezhda | 6000-5600 BC | 45 | 9 | 20 | 1 | 11 | 8 | 89 | — | —
Total |  | 290 | 64 | 22 | 7 | — | 30 | — | 28 | —

(descriptions of sites, radiocarbon dates, stratigraphic contexts, and sample descriptions are given in the Supplementary Information).

[Back]

[Back]

Figure 1.

Gas chromatograms of TLEs from pottery from the three regions of study showing high concentration of palmitic acid (C16:0) and stearic acid (C18:0): (a) Ecsegfalva 23, sample 16, (b) Grivac, sample 17, and (c) Yabalkovo, sample 38. IS – internal standard: C₃₄ n-tetratriacontane.

[Back]

[Back]

Figure 2.

Regional differences in Δ¹³C values ( = δ¹³C18:0 − δ¹³C16:0) of archaeological lipid residues in pottery. (a) proportions of lipid residues according to commodity type; (b) Δ¹³C values of archaeological lipid residues from sites in (b) the southern Balkans, (c) northern Balkans (d) the Carpathian basin (including published data for Ecsegfalva³²). Each data point represents an individual vessel.

[Back]

[Back]

Figure 3.

General trends in the exploitation of the most abundant mammalian taxa. Faunal data from early farming sites in the Aegean, the Balkans and the Carpathian basin, symbol size corresponds to relative body size of species, each symbol represents 10% of NISP. See Supplementary Table S2 for the complete dataset.

[Back]

[Back]

Figure 4.

Summary of climatic conditions in Southeast Europe. Dots indicate the location of faunal assemblages included in this study (see Supplementary Table S2); black dots – sites where caprines were represented with >60% of NISP of the large- and middle-sized domestic and wild mammalian taxa, white dots – caprines <60%. The dotted lines indicate regional boundaries between the Aegean/Adriatic littoral, southern Balkans, northern Balkans and Carpathian Basin. Numbers show location of sites where organic residue analyses were performed 1. Yabalkovo, 2. Nova Nadezhda, 3. Blagotin, 4. Divostin, 5. Grivac, 6. Alsónyék, 7. Ecsegfalva 23. (Climate data from: Hijmans, R. J., S. E. Cameron, J. L. Parra, P. G. Jones and A. Jarvis. 2005. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25(15): 1965-1978; http://worldclim.org/, Layer BIO6 = Min Temperature of Coldest Month, interpolations of observed data, representative of 1960–1990). The map was created using QGIS Version 2.8.9-Wien (http://qgis.org/de/site/).

[Back]