Mercer Family Plot at Bonaventure Cemetery. Photo © Ron Cogswell, Creative Commons
The analysis of the human remains from the megalithic tomb at Alto de Reinoso represents the widest integrative study of a Neolithic collective burial in Spain. Combining archaeology, osteology, molecular genetics and stable isotope analysis (87Sr/86Sr, δ15N, δ13C) it provides a wealth of information on the minimum number of individuals, age, sex, body height, pathologies, mitochondrial DNA profiles, kinship relations, mobility, and diet. The grave was in use for approximately one hundred years around 3700 cal BC, thus dating from the Late Neolithic of the Iberian chronology. At the bottom of the collective tomb, six complete and six partial skeletons lay in anatomically correct positions. Above them, further bodies represented a subsequent and different use of the tomb, with almost all of the skeletons exhibiting signs of manipulation such as missing skeletal parts, especially skulls. The megalithic monument comprised at least 47 individuals, including males, females, and subadults, although children aged 0–6 years were underrepresented. The skeletal remains exhibited a moderate number of pathologies, such as degenerative joint diseases, healed fractures, cranial trauma, and a low intensity of caries. The mitochondrial DNA profiles revealed a pattern pointing to a closely related local community with matrilineal kinship patterns. In some cases adjacent individuals in the bottom layer showed familial relationships. According to their strontium isotope ratios, only a few individuals were likely to have spent their early childhood in a different geological environment, whilst the majority of individuals grew up locally. Carbon and nitrogen isotope analysis, which was undertaken to reconstruct the dietary habits, indicated that this was a homogeneous group with egalitarian access to food. Cereals and small ruminants were the principal sources of nutrition. These data fit in well with a lifestyle typical of sedentary farming populations in the Spanish Meseta during this period of the Neolithic.
Once the Neolithic had been established in the Middle East the new way of life spread from the 7th millennium BC onwards via different routes into Europe. The earliest evidence of farming and animal husbandry on the Iberian Peninsula dates from as early as the 6th millennium BC and is found in the eastern Mediterranean coastal regions, southern Spain and to a lesser extent in the interior [1–5]. Towards the end of the 6th millennium farming had already spread throughout almost the entire peninsula . The Neolithic on the Iberian Peninsula can be divided into an early (5700−4500 cal BC), a middle (4500−4000 cal BC) and a late (4000−3000 cal BC) phase . The earliest evidence of pulses and cereals found so far was unearthed in the interior of Iberia at La Lámpara, La Revilla del Campo, and La Vaquera [2,4]. This early phase of agriculture already comprised a wide range of farming produce: wheat, barley, pulses (peas, lentils, beans) as well as vetch, flax, and poppy. The considerable variety of agricultural crops that have different requirements with regard to their cultivation, processing and use shows that these early farmers already had extensive knowledge of farming .
The Neolithic also introduced new funerary rites into Western Europe. The first so-called megalithic tombs were built in Spain, Portugal and France using large blocks and slabs of stone. Soon thereafter similar megalithic tombs appeared along the European Atlantic coast, in Britain and Ireland, Scandinavia and the European interior . In the interior of the Iberian Peninsula these tombs began to appear at the end of the 5th millennium BC  and spread throughout the entire region over the course of the 4th millennium BC . The funerary monuments were extremely varied in their construction . A characteristic feature of megalithic tombs is that these were frequented and used permanently and over an extended period of time as collective burial places for the members of a community whilst also serving as venues for ritual acts. Complex patterns of “treating” and “reburying” the skeletal remains have been identified in some megalithic tombs .
Skeletal remains are valuable biogenic archives for the study of the human past [10–11]. In this investigation we chose an integrative approach that combined archaeological and anthropological methods to explore the lifestyles and living conditions of the burial community. The osteological findings provided information about the minimum number of individuals (MNI), age, sex, body height, stress markers, diseases, traces of violence, and activity patterns [12–14]. Molecular genetic (mtDNA) data made biological relationships visible and revealed the genetic profile and relations of the individuals buried at Alto de Reinoso with other Iberian and European Neolithic populations [15–16]. Strontium isotope ratios of enamel and bone are geological tracers and were analysed as an indicator of mobility and change of residence . The 87Sr/86Sr ratios of the trace element strontium in bedrock vary according to the kinds and ages of geological units. The biologically available strontium is transferred through the food chain, stored in enamel during tooth development, and remains afterwards unchanged. Under the precondition of consumption of locally produced food, strontium isotope ratios in enamel that differ from the local baselines can identify non-local individuals in burial communities. The light stable isotope ratios of carbon (δ13C) and nitrogen (δ15N) reflect dietary habits . Isotope fractionation allows for estimations of the proportions of C3 and C4 plants in the human diet as well as the contribution of meat, dairy products, or fish. Other studies with a similarly comprehensive approach to reconstructing the past have been published in recent years for different regions and epochs [19–22]. However, similarly extensive presentations of the research results in collective graves are still generally quite rare and completely lacking for Spain [24–25].
The analysis of the skeletal remains from the megalithic tomb at Alto de Reinoso is therefore the most wide-ranging study of the Neolithic in Spain. First of all, it aims to characterise the population group who buried their dead in the megalithic tomb. Did it comprise several generations of a village-like community with a restricted territory and indications of kinship relations, or do the genetic, isotopic, and osteological characteristics point to the presence of people originating from various communities throughout the wider area and/or their mobile lifestyles? What were the specific living conditions of the population, including dietary habits as well as stress due to disadvantageous health conditions or interpersonal violence? Moreover, do the mtDNA data reflect population-dynamic processes over the course of the establishment and consolidation of the Neolithic on the Iberian Peninsula as well as long-distance contact and migration?
The Late Neolithic Site at Alto de Reinoso, Spain
The megalithic tomb at Alto de Reinoso is situated in sprawling hills between the villages of Monasterio de Rodilla to the north and Fresno de Rodilla to the south (Province Burgos, northern Spain). Located on a narrow platform in the Paramo highlands, the burial site dominates the entire surrounding area and is placed in a striking topographical location (Fig 1). The tomb is surrounded by the Northern Meseta landscape, a vast plateau in the interior of the Iberian Peninsula. This plateau has an average elevation of 600 to 700 metres and is flanked by mountain ranges that separate it from the coastal regions: in the north by the Cantabrian Mountains (Cordillera Cantabrica), in the north-east by the Iberian Mountains (Sistema Iberico) and in the south by the Central System. Highly diverse landscapes with strong contrasts in altitude within relatively short distances are characteristic of the Iberian Peninsula .
Fig 1. Map showing the geographical location of the Alto de Reinoso megalithic tomb in the Burgos province (Castilla y León, Spain) (graphic: Héctor Arcusa Magallón; Copyright by the authors). (Click Image to Enlarge)
The collective tomb had a diameter of approximately three metres and was excavated between 2006 and 2007 by archaeologists from the University of Valladolid, Spain. The primary aim of the rescue excavation was to prevent further destruction of the tomb. Only a few of the stones used in the construction of the mound had been preserved in situ. All the findings made during the excavation suggested that the tomb had originally been some kind of hut constructed in vegetal materials (house of the dead) that was eventually dismantled, closed off and finally turned into a monumental structure by erecting a stone mound above the ossuary (Fig 2).
Fig 2. Reconstruction of the original appearance of the megalithic mound of Alto de Reinoso (graphic: Héctor Arcusa Magallón).
The excavators identified three distinct layers in the area of the tomb. The uppermost Bronze Age layer, which had been disturbed by the plough during long-term agricultural activity around the grave in more recent times, contained numerous widely scattered stone remains as well as severely disturbed human bones and grave goods. The skeletal remains of two individuals dating from around 1700–1500 cal BC were found in this layer. These two individuals were not included in the study presented here.
The actual Neolithic ossuary containing a surprisingly high density of human bones in two burial layers (top and bottom) of very different states of preservation was discovered beneath this layer. The skeletal remains of the top layer were highly disturbed and included the upper body parts and skull of a male adult (Rein 30) and numerous concentrations of disarticulated human bones: groups of skulls (e.g. the isolated skulls Rein 17 and Rein 23 arranged next to the skull of individual Rein 8; S1 Fig), assemblages of long bones or even skulls intentionally framed by long bones. Most assemblages of skulls and postcranial bones appeared to have been arranged (e.g. Rein 29; see S2 Fig). Furthermore, a huge amount of disarticulated postcranial bones and a number of isolated skeletal parts were found in this layer. Most likely, these were the result of human interference in the context of rituals or were perhaps the result of bone dislocation whilst new bodies were placed in the megalithic tomb.
The bottom layer, on the other hand, revealed complete primary burials in situ, which showed hardly any signs of disturbance. It lay directly on the natural bedrock and contained six relatively complete (Rein 1, 2, 3, 4, 5, and 10) and six half skeletons (Rein 6, 7, 8, 9, 11 and 12) in anatomical order as primary burials (S1 Table and S3 Fig). The individuals had been interred in a crouched position, as it was typical in many Neolithic archaeological contexts : six individuals lay on their left side (Rein 1, 2, 4, 5, 8 and 11) whilst two were found lying on their right side (Rein 3 and 10). However, it was not possible to reconstruct the exact burial sequence of the bodies in the tomb. The skeletal remains had lain relatively undisturbed until being covered by other burials, which provided them with additional protection against relocation. Due to the collective utilisation of the tomb characterised by successive burials whereby older corpses were displaced in favour of more recent ones, individualisation of the human remains was only possible for some of the skeletons. Hence, the reconstruction of the demographic profile of the burial community proved challenging.
The archaeological material recovered from inside the ossuary was not overly abundant and included both personal adornments (stone necklace beads, wild boar tusk pendants) and grave goods (polished stone axes, flint blades and microliths, bone spatulae) (S4 Fig). Some had been intentionally deposited beneath or near the skulls (flint blades and stone polished axes). Two bone spatulae, one of which was decorated, and a stone bead had been placed beneath a skull (Rein 29), and it had also been surrounded by a triangular frame consisting of three long bones (see S2 Fig). Small artefacts such as microliths or beads were easily removed and repositioned along with the bones; one exception to this rule was a lignite bead found on the neck of individual Rein 1. Remarkably, there was a distinct lack of ceramic vessels in megalithic tombs, at least during the earlier period of their use. Researchers assume that this was due to ritual reasons . The archaeological material suggested a date for the megalithic tomb at Alto de Reinoso in the Late Neolithic. Overall the finds exhibited strong chronological coherence and were all artefacts typically found in megalithic tombs throughout the Meseta during the 4th millennium BC.
Radiocarbon dates were obtained from three human bone samples recovered from different layers in the ossuary, giving a maximum chronological range of between 3770 and 3539 (2 sigma cal BC) (S2 Table). The OxCal v4.2.4 software  allowed us to narrow down the actual time span represented by the radiocarbon dates to a period of 60 to 80 years, or roughly three generations (S2 Table and S5 Fig). Both the archaeological finds and the radiocarbon dates are characteristic of collective burials from the Late Neolithic .
In principle, we may state that collective burials were abundant at the beginning of the 4th millennium BC and exhibited diverse architecture ranging from simple non megalithic barrows, that is to say lacking colossal stone structures, to large passage graves as simple dolmens and allées couvertes . Alto de Reinoso is a good example of a non-megalithic barrow, as the result of the monumentalization of a more basic structure, probably a burial chamber made of wood, mud and other organic materials, the actual shape of which is difficult to ascertain. This structure was finally dismantled as part of a closure ritual, and then marked and monumentalised by erecting a small stone barrow covering the remains of the ossuary . Where a tomb had once stood, a monument to the memory of the ancestors was now erected as a permanent reference to their presence in the symbolic and cultural landscape. Similar examples of this type of monument have been discovered at La Tarayuela (Ambrona Valley, Soria) and at El Castillejo .
Materials and Methods
This study deals with archaeological skeletal material and adheres to the Helsinki Declaration on Ethical Principles for Medical Research Involving Human Subjects. The excavation licenses were issued by the General Director of Cultural Heritage of the Junta de Castilla y Léon (Spain) and are stored in its archives. Following excavation campaigns 2006/2007, the skeletons were transferred to the former Bioarchaeology group at the University of Mainz, Germany for investigation. At present all archaeological remains, including the human bones, recovered at the site of Alto de Reinoso, are stored in the Prehistory and Archaeology lab of the University of Valladolid. The inventory list of the bones in their current storage place is identical to the list supplied in S1 Table of the supplemental information.
A basic concern of anthropological examination of human remains from megalithic monuments is the dense occupation of the tombs both by completely preserved individuals and partial burials as well as isolated bones. One of the preconditions of investigating the demographic structure of a burial community is the analysis of the minimum number of individuals (MNI) which serves as a starting point for all further studies including osteological determination of age and sex, reconstruction of body height, investigation of diseases, trauma and stress markers as well as molecular genetics and isotope analysis. The techniques used in this interdisciplinary integrative study were based on established methods and these are described in the supporting material (S1 Text). A varying number of individuals were available for each section of the examination, depending on the part of the skeleton required to carry out the analysis (Table 1 and S1 Table).
Table 1. Results of osteological, mtDNA and isotopic analysis for the human individuals from Alto de Reinoso. (Click Image to Enlarge)
In order to reproduce the aDNA results, two or three samples were taken from each of the available individuals (n = 27) for molecular genetics analyses. We generated mitochondrial profiles by analysing the hypervariable regions I and II (HVR I, HVR II) and coding region single nucleotide polymorphisms (SNPs) of the mitochondrial genome (S7, S8 and S9 Tables). Protocols for sample preparations, extraction, amplification setups, and contamination control were applied as described previously [15, 16, 22, 32]. During genetic sample preparation, the tooth enamel—which is not required for molecular analyses—was separated from the dentine and passed on to undergo strontium isotope analysis. This investigation was based on 36 enamel samples from the permanent dentition of 19 individuals. In most cases both an early and a late mineralising tooth were sampled from each individual in order to gain information on potential changes of residence during childhood. Individuals whose bones were too poorly preserved to determine their sex and age at death were excluded from the sampling process. The analytical procedure included the mechanical abrasion of all surfaces and remaining dentine, chemical cleaning with buffered acetic acid, ashing, strontium separation with Eichrom Sr-Spec resin and determination of Sr concentrations with a Quadrupole-Inductively Coupled Plasma-Mass Spectrometer (Q-ICP-MS) and the 87Sr/86Sr ratios with a Multi Collector-ICP-MS (VG Axiom) at the Curt Engelhorn Centre Archaeometry gGmbH in Mannheim, Germany following the procedures described in Knipper et al. . Details on sample preparation and analyses are given in S1 Text. In addition to the tooth samples, three human bones from the grave and two samples from domestic cattle (Bos taurus) also derived from the grave provided comparative data for the interpretation of the 87Sr/86Sr ratios. During the osteological examination bone samples were also taken for collagen extraction with the aim of analysing the carbon and nitrogen isotope ratios (δ13C and δ15N) in both humans (n = 29) and animals, including cattle (Bos taurus) and a rabbit (Oryctolagus cuniculus) (n = 4). The taphonomic alterations of the rabbit bone suggest that it belonged to the same burial context and time period as the human remains. The analytical procedure included mechanical cleaning, demineralization, removal of humic acids, gelatinization, filtration, concentration of the long-chained collagen by ultrafiltration and lyophilisation . The carbon and nitrogen contents were determined by an elemental analyzer (vario EL III, Elementar Analytical Systems) and the isotope ratios by an IsoPrime High Performance isotope ratio mass spectrometer (IRMS; VG Instruments) at the Department of Organic Chemistry at the University of Mainz. More information on sample preparation and analysis is given in S1 Text. The sampled skeletal remains were given laboratory assignations ranging from Rein 1 to Rein 33 (Table 1 and S1 Table). In cases where sampling was not possible the bones were only listed with an inventory number (Inv.). Isolated postcranial bones are not listed, but recorded under its inventory number.
Results and Discussion
The skeletal remains near the surface in particular bore numerous traces of taphonomic processes including high fragmentation rates, evidence of severe weathering (bleaching, surface changes, micro-cracks, and calcification) as well as root etching. Apart from one exception the absence of gnawing or bite marks on the skeletal remains proves that the dead were quickly placed in the megalithic tomb and that it was inaccessible to carnivores.
Minimum number of individuals (MNI).
The minimum number of individuals (MNI) was determined by counting the most abundant parts of the skeletons recovered from the site: crania, pelvis, and femoral diaphysis . Based on right femora and elements of the crania an MNI of 47 can be estimated (S3 Table). However, only 38 near-complete crania or mandibles with teeth were found (21 adults and 17 subadults) which corresponded to a deficit of at least nine adult crania in the burial context compared to the MNI of 30 adults. One possible explanation is that the missing crania were removed during the phase of active utilisation of the burial mound, perhaps in the context of ancestor worship.
Sex and age estimation.
Sex estimation of human remains is commonly based on morphological traits of the pelvis and cranium or on a metric evaluation following  and . At Alto de Reinoso these skeletal elements were only present in individuals from the largely undisturbed bottom layer, so that numerous individuals from the top layer could not be sexed due to poor preservation, and/or due to significant displacement of the bones. The sex estimation of 27 (26 adults, one juvenile) pelvic bones attested to a balanced distribution of both sexes (11 males [40.7%], 12 females [44.4%] and four undetermined adults [14.8%]) (S4 Table). In contrast, apart from a few instances in the bottom layer, the sex estimation of most adult crania was not possible because of poor preservation and high fragmentation. There were more males (n = 9 [42.9%]) than females (n = 5 [24%]) in the sample, although the adults of undetermined sex (n = 7 [33.3%]) could of course include more women than men (S5 Table). Long bones which could be sexed due to a clear affiliation to a sex estimated individual exhibited remarkable sexual dimorphism. It is worth noting that sex determination based solely on crania, without taking into account the pelvis results in an imbalance in the sexes.
In complete skeletons biological age can be expressed either as the dental age or the skeletal age [35–42]. In the case of Alto de Reinoso, at least 30 adults (>20 y) and 17 subadults (≤20 y) had originally been buried in the collective grave [S3 Table]. Crania and mandibles with teeth formed the basis for more precise individual age examinations [S6 Table].
The classification of individuals in different age groups followed Martins’ system . Based on 21 adult and 17 subadult individuals for whom ectocranial suture closure  and tooth wear [37;40] or tooth mineralisation  could be assessed, twelve individuals [31.5%] were estimated to have been between 21 and 40 years old (age group adult) and four individuals [10.5%] were classified as “mature” (41–60 y). Five further non-infant individuals [13.2%] could not be firmly assigned to either of the age groups. Subadults were represented by four individuals of the age group infans I (0–6 y, [10.5%]), five infans II (7–12 y, [13.2%]), and eight juveniles (13–20 y, [21.1%]). It must be noted in this context that the deposition below ground has a much more significant impact on the preservation of the skeletal remains of infants than on the more robust bones of adults , which results in a situation where the bones of some of the children are either not preserved or can at least no longer be individualised. Therefore, the burial community was certainly larger than the MNI would suggest, especially since the individuals of the infans I age group usually comprise the largest proportion of deaths in a prehistoric sample. The deficiency in infants may also have resulted from special burial rites which excluded them from collective burials while older children were interred with the adults .
The variation in the adult skeletal morphology and metrics is determined by several factors including health and nutrition during a person’s growth and development, as well as genetics, sex, geography, climate and social conditions. In principle, numerous bones are suitable for the determination of body height [14,46–47]. The long bones of the legs and arms offer the best correlation with a person’s real stature, especially if more than one bone is available for evaluation. Due to the differences in body size between the sexes and between different ethnic groups and populations very different formulae are used in calculating body height [48–52]. Unfortunately, as seen in the other examinations of the skeletons, the state of preservation of the long bones was not very suitable for body height estimation. Taking into account only long bones of clearly sexed individuals, the body height of the males (n = 6) was 159.4 ± 2.0 cm on average (ranging between 156 and 163 cm; ). The females (n = 3) exhibited an average body height of 150.0 ± 2.4 cm (ranging between 147 and 152 cm), i.e. nine centimetres shorter than the males. When using the formulae proposed by Breitinger  for males and by Bach  for females the mean values for both sexes were significantly higher: 162.7 ± 4.8 cm (ranging between 160 and 165 cm) for males and 157.9 ± 4.1 cm (ranging between 154 and 160 cm) for females. Such differences between the sexes are quite common in homogenous populations—seen as larger groups of individuals without subpopulations—where most men are taller than the tallest woman. The low variation within both sexes also suggests that Alto de Reinoso was a closely related community . However, this result must be regarded as an approximation because the total sample size is very small; in seven individuals only one length measurement, mainly of the femur, could be taken into consideration.
Palaeopathology of bones and teeth.
Substantial surface alterations in the highly fragmented skeletal material were challenging for the diagnosis of diseases, traumata and stress markers. The following description of the pathological findings is therefore not complete with regard to the frequency or the spectrum of pathologies. Indications of pathologies were recorded regardless of whether or not they could be individualised; in cases where it was possible to individualise them, they were listed by inventory number. Overall, the palaeopathological findings were typical of Neolithic farming populations [53–55]. The joints of the spine showed an entire range of degenerative alterations: osteoarthritis (Inv. 3809; Inv. 3589), Schmorl’s nodes (Rein 4), spondylarthrosis and spondylosis (Rein 6, 8, 10, 12, 30) as well as spondylitis (Rein 10). With two fractured ribs (Rein 5, 22), a healed femur fracture (Inv. 1246–47) and a tibia fracture with pseudo-arthrosis (Inv. 1526–27) as well as a hint of blunt force trauma to the cranium (Rein 30), the number of commonly found traumata and violence-related cranial traumata was within the expected low range for the populations of that period . Exceptional enthesiopathies were only diagnosed on the leg bones of Rein 11 and were interpreted as an indicator for more physical activity . Unspecific inflammations were also found in some cases (Inv. 2136–37; Inv. 3383–94; Rein 19; Rein 29; Rein 33). This observation coincides with a general trend showing that the frequency of such lesions increases at the beginning of the Neolithic period .
Due to the limitations with regard to individualisation and the large number of isolated teeth it was not possible to determine the frequency of caries (number of affected individuals compared to the total number of individuals). Caries intensity of 14.5% (number of decayed teeth compared to the number of preserved teeth including antemortem tooth loss) among children and adult individuals (loose teeth of unknown individualization are included) was quite low compared to Middle Eastern and European Neolithic populations [55,59–60], and especially compared to other Late Neolithic collective burials [25,54,61].
If compared to the contemporary sites at San Juan ante Portam Latinam  and Cova del Pantá de Foix , the caries intensity of only the permanent dentition at Alto de Reinoso, which was accounted for 14.3% (including antemortem tooth loss), fits in well with 13.2% at San Juan Portam Latinam and 15.3% at Cova del Pantá de Foix, although it is unclear whether antemortem tooth loss was taken into account for Cova del Pantá de Foix or not.
However, the unusually high caries intensity in deciduous teeth (17.8%) suggests that nearly all subadults suffered from tooth decay (exclusively molars) in their primary dentition. With 13.4% affected deciduous teeth similar values exist at San Juan ante Portam Latinam  while the collective grave of Benzingerode, Germany, did not contain a single carious deciduous tooth . Several datasets on caries intensity in deciduous teeth in earlier Linearbandkeramik communities in Germany exist: 5.4% at Aiterhofen  and 2.7% at Wandersleben . The reasons for the rather high percentages of caries in deciduous teeth seem to be associated less with the socioeconomic status of the population but rather with specific dietary habits, especially the consumption of carbohydrate-rich sticky food such as porridge but also due to prolonged breastfeeding [66,67].
With regard to stress markers, enamel hypoplasia [66,68] was observed in 11.8% of the permanent teeth (canines in particular) and in 10.2% of the deciduous teeth (molars in particular), although the poor state of preservation did not allow for an individualised analysis. When comparing the incidence of enamel hypoplasia in the permanent dentition at Alto de Reinoso with 36.5% at Cova del Pantá de Foix, significantly fewer teeth showed enamel hypoplasia at Reinoso. A similar comparison between the deciduous teeth was not possible. Studies from the Portuguese collective tombs at Feteira II, Cova da Moura, Paimogo I and Cabeço da Arruda I , which dated from the Neolithic to the Chalcolithic, revealed less enamel hypoplasia than was observed at Alto de Reinoso but were based on small sample sizes. Cribra orbitalia could not be analysed due to considerably altered cranial bone surfaces and severe crushing of the crania.
The HVR I and coding region single nucleotide polymorphisms (SNPs) of the mitochondrial genome (GenoCoRe22; ) were successfully reproduced for 26 out of the 27 individuals analysed (www.ncbi.nlm.nih.gov/genbank; accession numbers KT868907-KT868932), resulting in an extremely high amplification success for ancient DNA analyses especially in relation to studies from south-western Europe (S1 and S7 Tables). Furthermore, H –plex analyses  defined the specific subhaplogroup of an individual belonging to haplogroup H (S8 Table). The individuals buried in the megalithic tomb at Alto de Reinoso belonged to mtDNA (sub-) haplogroups H3 (n = 1; 3.7%), J (n = 1; 3.7%), K (n = 8; 29.6%), T2 (n = 4; 14.8%), U4 (n = 1; 3.7%), U5b (n = 7; 25.9%), V (n = 2; 7.4%), and X (n = 3; 11.1%). The dominance of haplogroups K, T2 and U5b in particular can be ascribed to the large number of shared lineages among the individuals buried in the tomb. More than half of all individuals analysed have characteristics in their mitochondrial profile that match one or two of the other individuals, pointing to maternal relations within the burial community. Shared lineages for all markers analysed could be found in individuals Rein 17, 20 and probably 24 [K], in Rein 5, 18, and 19 [T2b], between Rein 1, 15, and potentially 2 [U5b3], and also between Rein 12 and 32 [K1a], Rein 22 and 31 [U5b], and Rein 6 and 9 [V] (S1 Table). With reference to the layers in the tomb, shared lineages can be observed between individuals from the same layer (e.g., Rein 6 and 9; Rein 2 and 1 [bottom layer]; Rein 17 and 24; Rein 22 and 31; Rein 18 and 19 [top layer]), but also in individuals found in different layers (e.g., Rein 1 and 2 with Rein 15, Rein 5 with Rein 18 and 19 and Rein 17 and 24 with Rein 20; Fig 3). All in all 17 mitochondrial links were identified among the 26 individuals.
Fig 3. Superposition of different layers of the Neolithic ossuary indicating the individuals with the same genetic profile (graphic: Héctor Arcusa Magallón). (Click Image to Enlarge)
In order to investigate genetic connections between the Alto de Reinoso site and other 18 prehistoric European populations, comparisons were made with published data from elsewhere in West Europe on one hand [71–80], and from nine prehistoric Central European groups on the other [16,32,81–89] by means of a principle component analysis (PCA) and Ward type hierarchical clustering (Fig 4a and 4b). The first and second principal components (PCs) account for approximately 38% of the overall variability of the dataset on the PCA. The closest groups to Alto de Reinoso are the Early/Middle Neolithic Bernburg (BEC) populations in Central Europe, the Late Neolithic south French Treilles group (TRE) and the Chalcolithic El Mirador community (Fig 4b). The north Spanish and the central Portugal Neolithic groups together with the south-western European hunter-gatherers (HGSW) are separated from the main cluster along the N, U, and H haplogroups. The Central European hunter-gatherer metapopulation (HGC) is clearly separated from the Iberian groups and Central European Neolithic too, indicating population discontinuity during Neolithisation of Central Europe [16,86,89].
Fig 4. Principal component analysis (A) and Ward type hierarchical clustering (B) with 19 Western and Central European populations. / On both image parts south-western European hunter-gatherers [HGSW] and Central European hunter-gatherers [HGC] are symbolised by grey hexagons, Early Neolithic Iberian groups (Central Portugal [CP_EN], north Spain [NS_EN], and northeast Spain [NE_EN]) by red circles and Middle/Late Neolithic groups by pink star and circles (Alto de Reinoso [Reinoso], north Spain [NS_MLN], northeast Spain [NES_MLN], south French Treilles group [TRE]). Magenta-coloured circles indicate the chalcolithic El Mirador cave [El_Mirador] and Late Neolithic-Chalcolithic Central Portugal [CP_LN_CH]. Early Neolithic groups in Central Europe and Carpathian Basin are represented by brown circles (Starcevo culture [STA], Linear pottery culture in Transdanubia [LBKT], Linear Pottery culture in Central Europe [LBK], Rössen culture [RSC], Schöningen group [SCG]), Middle Neolithic groups in Central Europe (Salzmünde Culture [SMC], Baalberge Culture [BAC], Bernburg Culture [BEC]) are in orange. (B) Cluster significance is given as AU (Approximately Unbiased) p-value in %, which is computed by multi-scale bootstrap resampling with 10,000 replicates (graphic: A. Szécsényi-Nagy). / Click Image to Enlarge
Alto de Reinoso and the Early Neolithic population in north-eastern Spain appear to consist of a mixture of haplogroups typical of other Neolithic populations in Spain and Portugal that are mainly characterised by haplogroups H, N and U5b, and haplogroups V, J, K, N1a, T2, and X that are typical for Early and Middle Neolithic groups in Central Europe. However, the PCA reveals a separation between the Early/Middle Central European cluster (including Alto de Reinoso) and the cluster composed of the Neolithic-Chalcolithic Portugal, Neolithic north Spain, and the north-eastern Spanish Middle Neolithic population, where the latter cluster shows marked mixture with populations of hunter-gatherer legacy. Furthermore the PCA indicates common traits between the Central European Early/Middle Neolithic genetic substrate, the Early Neolithic population of northeastern Spain, the Late Neolithic Reinoso and Treilles community and the population of the Chalcolithic El Mirador site, probably originated from a common Near Eastern ancestry.
The 87Sr/86Sr ratios of a total of 36 human enamel samples (s) from 19 individuals (n) varied between 0.70893 and 0.70948 (mean: 0.70920 ± 0.00027; 2 σ) (Fig 5 and S1 Table). The ranges of the adults of both sexes and the subadults largely overlapped. Most likely due to the small sample size the data of the adult females (0.70899–0.70933, n = 3, s = 6) appeared less variable than those of the subadults (0.70891–0.70951, n = 8, s = 14) and the adult males (0.70897–0.70957, n = 7, s = 14). Indications pointing to local baseline values at the site itself came from three human bones (0.70924–0.70952) and two samples of domestic cattle enamel (0.70918–0.70927) (S10 Table). Their average, at two standard deviations, exhibited a range of between 0.70912 and 0.70952. The teeth of a number of individuals had 87Sr/86Sr ratios just below the lower limit (Rein 5, 10, 14, 15, 18, 23, 29, 30) or above the upper limit of this range (Rein 3). However, because the number of comparative samples was very small and may lead us to underestimate the variation of the local baseline values, these numbers provide no clear indication with regard to non-local individuals. Moreover, the overall data distribution of the human enamel was continuous and lacked any results that clearly diverged from the main data cluster.
Fig 5. 87Sr/86Sr ratios of human and faunal enamel and bone samples from Alto de Reinoso and from Algría-Dulantzi in comparison. / The data from Reinoso are grouped by age and sex and sorted by the 87Sr/86Sr ratios of deciduous teeth and first molars within each group. The dashed lines illustrate the isotope ratios of the local biologically available strontium in human bone and cattle enamel as calculated from the average plus and minus two standard deviations. (Data for Algría-Dulantzi ) (graphic: C. Knipper). / Click Image to Enlarge
Nevertheless it is worth mentioning that in addition to the 87Sr/86Sr ratios of their first molars beyond the local baseline range, the differences between the early and the later forming teeth of Rein 23 (Δ87Sr/86SrM1-M3 = 0.00042), Rein 3 (Δ87Sr/86SrM1-M3 = 0.00031), Rein 30 (Δ87Sr/86SrM1-M2 = 0.00026), and Rein 5 (Δ87Sr/86SrM1-M2 = 0.00019) were larger than those of all individuals with both measurements within the local range. As these differences are larger than the analytical error they may point to a shift in diet catchment areas during childhood, even though it would have involved localities with very similar geological conditions to those at Alto de Reinoso. The site is surrounded by calcarenites and conglomerates of the Lower Neogene , which were also used in the construction of the tomb itself and are widely distributed around the site (S6 Fig).
Only a few studies on bioavailable 87Sr/86Sr isotope ratios from archaeological sites on the Iberian Peninsula have been published so far [91–93]. A rich comparative dataset from the area, which includes human enamel, bone, dentine, and faunal samples was gathered at the early medieval cemetery of San Martín de Dulantzi (Algría-Dulantzi, Álava) some 90 km north-east of Reinoso, where Cretaceous limestone and marls prevail . This study illustrates the existence of localities with less radiogenic baseline values than those observed at Reinoso and the possibility of identifying distinctly non-local individuals in northern Spanish burial contexts.
The results of the stable carbon and nitrogen isotope analyses of the human and animal bones are presented in Table 1, S1 and S10 Tables. Most of the human bones and four animal samples (three cattle, one rabbit) analysed met the quality indicators for collagen. Four human samples yielded less than 1% of collagen. Because the other quality criteria (C %, N %, and C:N ratio) of these samples were in the accepted ranges [94–96] the samples were included in our considerations. Their lower collagen yields are likely to have resulted from a substantial loss of higher degraded collagen during the ultra-filtration step . The human stable isotope ratios ranged from -20.2 to -19.2 ‰ (average -19.6 ± 0.3 ‰) for δ13C and from 9.1 to 11 ‰ (average 9.9 ± 0.4 ‰) for δ15N (Figs 6 and 7). The average values of all adults were -19.7 ± 0.3 ‰ in δ13C and 9.8 ± 0.4 ‰ in δ15N (n = 19). The δ13C mean ratios for males were -19.7 ± 0.3 ‰ (n = 9) and those for females -19.5 ± 0.2 ‰ (n = 3). The δ15N mean values of the males (10.1 ± 0.5 ‰) were slightly higher than those of the females (9.7 ± 0.1 ‰). The stable isotope ratios from the faunal samples varied from -21.5 to -20.7 ‰ for δ13C (average -21.2 ± 0.4 ‰) and from 4.5 to 6.1 ‰ for δ15N (average 5.1 ± 0.7 ‰) (Figs 6 and 7 and S10 Table). The data for the different human age groups are presented in Fig 7. All the groups tested were normally distributed (Kolmogorov-Smirnov-test for normal distribution) and an analysis using the ANOVA and Post-Hoc Tuckey HSD-tests showed no significant differences between the sexes or between the ages.
Fig 6. Carbon and nitrogen isotope ratios of the human and faunal samples from Alto de Reinoso in comparison to the isotope ratios of potential major food sources. / Comparative data are taken from [107–108, 135, 144–146] (graphic: C. Knipper). / Click Image to Enlarge
Fig 7. Scatter plot of δ13C and δ15N values and faunal collagen data from the Alto de Reinoso collective tomb and data from comparable Mesolithic and Neolithic sites. / 1: Valencia: Santa Maria, Penya del Camtador, Clingel des Mans Nou (Spain ), 2: El Collado (Spain [137), 3: Costamar (Spain ), 4: Cova de la Pastore (Spain ), 5: Mean values from Casa da Moura, Gruta do Caldeirao, Gruta da Feteira, Gruta da Fontainhas, Roche Forte II, Cerca do Zambujal, Melides (Portugal ); 6: Fateira II (Portugal ), 7: Paimogo I (Portugal ), 8: Cabeco da Arruda (Portugal ), 9: Languedoc area, Les Crés, Montou, (France ), 10: Garonne area, Cugnaux, Villeneuve-Tolosane, Narbons (France ) (graphic: C. Knipper). / Click Image to Enlarge
The stable isotope data of the humans from Alto de Reinoso showed values that are typical of a mixed diet in a food web based on C3 plants. There is no indication for substantial millet or marine diet intake, which would both have resulted in higher δ13C values. The latter is not surprising considering the location of the site in the Meseta .
To reveal better information on the protein sources we calculated the difference between the faunal and human values. The average offsets between cattle and humans were +1.7 ‰ in δ13C and +5.2 ‰ in δ15N. Both the differences in the carbon and in the nitrogen values were higher than the enrichment typical of one trophic level, which is 1–2 ‰ for carbon and 3–5 ‰ for nitrogen [99–103]. This does not exclude the consumption of beef, but suggests that foodstuffs with higher δ15N values have contributed substantially to the human diet. The differences between the isotope ratios of the rabbit collagen and the average of the humans were +1.1 ‰ in δ13C and +3.7 ‰ in δ15N. These numbers are closer to the normal trophic level values even though the overall context makes it unlikely that rabbits contributed largely to the animal protein supply.
Because cattle are unlikely a primary source of animal protein, other animals with higher δ13C and δ15N values such as pig, goat and sheep and their dairy products as well as possibly also freshwater fish should be considered. However, samples of these species were not available in the megalithic tomb at Alto de Reinoso. Therefore, we drew on analyses from other sites on the Iberian Peninsula, even though their explanatory power for our data was perhaps limited due to spatial and chronological differences. Collagen of wild boars from the Mesolithic site at Clingel del Mas Nou showed δ13C values of -19.9 ± 1.1 ‰ and δ15N values of 6.6 ± 1.0 ‰ (n = 6) . Pigs from Late Neolithic sites in Portugal with -20.5 ± 0.6 ‰ δ13C and 7.1 ± 2.3 ‰ δ15N (n = 8) revealed even higher nitrogen isotope values. The ovicaprine collagen from Late Neolithic Portuguese sites showed mean values of -20.7 ± 0.3 ‰ δ13C and 5.6 ± 0.4 ‰ δ15N . The data ranges of these samples fit in better with the isotope ratios of animals that could have contributed protein to the human diet consumed at Reinoso. Moreover, archaeological findings suggest that goat and sheep were the most important domestic animals in Spain during the Early and Middle Neolithic [105–106]. Goats and sheep could have been among the most essential basic foods consumed at Reinoso. Moreover, consumption of freshwater fish is possible, but due to their often very variable stable isotope data [107–108] hard to estimate without contemporary comparative samples from the Meseta.
Durum wheat, hulled and naked barley and faba beans were probably the main sources of plant protein [109–110]. Aguilera and colleagues  analysed the nitrogen isotope compositions of naked wheat and naked barley recovered from the Los Castillejos site which dated from between 4000 and 2500 BC. The nitrogen isotope values of barley of between 4 and 6 ‰ remained comparatively stable over time, while the wheat varied more, from 2 to 7 ‰. The nitrogen isotope values of plants can be strongly influenced by fertilising and by arid climate conditions . The data published by Aguilera and colleagues  clearly indicate that it is essential to consider plant protein data when examining the dietary habits of the Alto de Reinoso population. Crop plants provided a mean nitrogen isotope value of 4.8 ‰ which is highly similar to the mean nitrogen values of the cattle. It is therefore hardly possible to distinguish between plant and animal derived protein in the diet.
Integrative Analysis and Interpretation
Because of its remote elevated location (Meseta) the economy at Alto de Reinoso would have been characterised by small-scale farming and herding. The economic subsistence and the remote geographical location would obviously be reflected in the results. Moreover, the short period of use of the tomb and the preservation of the skeletal remains, which was of a sufficient quality for most of the analyses carried out, provided an opportunity to generate scientific data that would allow us to ascertain whether the subsistence conditions had had an impact on the social structure within the community, what it could have looked like and what dynamics may have developed. Besides the standard osteological methods to reconstruct the fundamental “biographical” data, this approach requires the use of modern methods of analysis, several of which were in fact employed in this study.
How prehistoric communities dealt with death is reflected in predetermined communal rituals. These also incorporated the place where the dead were laid to rest. The Neolithic megalithic tombs in Europe represent a particular type of collective burial. Compared to individual graves in a cemetery, the communal use of a “house of the dead” more clearly illustrates the bond within a community where the individual is not at the forefront, in life or in death . It was not until the Chalcolithic period and the Bronze Age that social differentiation increased and individuals and groups came to the fore more and more clearly .
According to the archaeological finds and features as well as the radiocarbon dates the collective tomb was used by a small farming community for a period of three generations at most (generation = 25 years). This rather short period of use appears typical of the megalithic ossuaries in the Iberian interior, as seen at La Tarayuela, La Velilla, La Sima (especially phase II), La Vega and Los Zumacales ; all these tombs were used by only a few generations (no more than four). Given the estimated period of use, the MNI calculated at 47 for Alto de Reinoso suggests that an average of one or two people died every year; this number would include infants, which were largely absent. The MNI is within the range typically observed at Neolithic megalithic tombs situated relatively close to Alto de Reinoso such as La Peña de la Abuela , La Sima  and La Tarayuela , and also at other Spanish sites further afield such as San Juan ante Portam Latinam  in the Basque country, Cerro Virtud in Almeria  and Cueva de Malalmuerzo in Andalucia  and at Late Neolithic sites in Portugal . All the data and findings generated at Alto de Reinoso in relation to the tomb, its size, the number of burials and its period of use attest to relatively uniform funerary rites at that time in the Meseta and (in some cases far) beyond.
Is the demographic composition of the community representative?
From the point of view of palaeodemographics a skeletal sample is considered representative if the adults exhibit a natural stratification regarding age and sex distribution and if the proportion of infants is not below 20% due to the high mortality rates of children under 1 year [120–121]. Because the data gathered at Alto de Reinoso did not fulfil these criteria, the sample would at first glance not be deemed representative. The percentages of adults and subadults at Alto de Reinoso were offset in favour of the adults because the infants and small children in particular were clearly underrepresented (S3 Table). The lack of bones of the age group infans I can generally be attributed to taphonomic processes, particularly in collective burials, where a lot of disturbance occurs over their period of use. At Alto de Reinoso the group was represented by individual bones but could not be accounted for as a whole. In a comparative study on the demographic structures in megalithic graves in northern Spain Fernández-Crespo and de-la-Rúa  presented results similar to those found at Alto de Reinoso with a shortage of children less than five years, although they also pointed to examples of a lack of mature and senile adults. However, compared to Cerro Virtud , La Peña de la Abuela , La Tarayuela  and La Sima  the tomb at Alto de Reinoso contained more children and adolescents, whilst yielding fewer compared to San Juan ante Portam Latinam  and Cueva de Malalmuerzo .
The age determinations yielded mainly adolescents (n = 8) and young adults (n = 12), whilst older adults (> 41 y) were underrepresented. However, taking into account the adults whose ages could not be ascertained (n = 7) could potentially increase the number of mature individuals. The ratio of sexes which favoured the males (determination based on crania) was put into perspective by the examinations of the pelvic bones which indicated a rather more equal distribution (S4 Table). Therefore, the predominance of males in megalithic graves as has been suggested due to comparable data gathered at Cerro Virtud , La Peña de la Abuela , La Tarayuela  and San Juan ante Portam Latinam  was not observed at Alto de Reinoso. The idea of the predominance of males was already discredited elsewhere . Based on the Las Arnillas dolmen where more males than females were identified Delibes de Castro argued for a possible social relegation of women . He believed that, besides problems linked to pregnancy and childbirth, women may have had a lower life expectancy due to social discrimination, which could perhaps even have impacted on their diet.
The striking lack of skulls can again be explained by taphonomic processes. Some of them were found as part of groupings with isolated long bones elsewhere in the tomb which, however, does not explain the shortfall. We would suspect that skulls were removed from the tomb, perhaps as part of ritual acts within the context of ancestor worship. The fact that children and young adults would not be considered ancestors could explain the relative underrepresentation of older individuals in the tomb. Numerous examples of ancestor worship using skulls in pars pro toto ceremonies are known from the Neolithic period [124–128].
The structure and dynamics of the community at Alto de Reinoso
The best way to answer questions regarding the social structure and potential dynamics within the Alto de Reinoso community, beyond the demographic data, is to study the molecular profiles. Taking into account the MNI, the number of burials in the two layers and its period of use, the bottom layer probably contained one generation, whilst the top layer is more likely to have included members of two generations. This means that all individuals from the same layer were separated by one generation at most. Reproducible genetic profiles were generated for 9 individuals from the bottom layer and 17 individuals from the top layer. The fact that, according to the molecular genetic analyses, 15 of the individuals buried at Alto de Reinoso shared their mtDNA haplotypes with at least one other person is not surprising given the limited size of the community, nor should it come as a surprise that pairs of individuals with identical haplotypes were found buried directly beside each other in the less disturbed part of the tomb, which suggests very close genetic relations (Fig 3). The fact that shared maternal lineages generally occurred between individuals within both layers and also between individuals from different layers is population-biological evidence for a certain degree of continuity within the group. On the other hand, some of the genetic profiles only occurred in one of the layers.
One could argue that the lineages that were especially frequent at Alto de Reinoso—mainly of haplogroups K and T2—were also identified at other Neolithic sites in Spain (Cova d´Avellaner, Can Sadurní, San Pau de Camp, and Chaves; [74,77]) and therefore perhaps simply reflect common Neolithic lineages instead of mitochondrial kinship. However, these lineages were not identified as frequently by any of the prehistoric aDNA studies carried out elsewhere and the rCRS lineage of haplogroup V (Rein 6 and 9) was exclusively found at Alto de Reinoso. The fact that a total of six different mitochondrial lineages occur in two or three individuals argues against random haplotype matches and for the existence of maternal kinship among the individuals buried at Alto de Reinoso. Given the period of use of the tomb one would not expect more than two or three members of any particular nuclear family to have been buried in the tomb, and these may represent successive generations. In summary, the overall picture is that of a closely related community that did not experience any noteworthy changes over the period of time the collective tomb was in use.
The absence of dynamic changes within the group is further confirmed by the 87Sr/86Sr isotope ratios. There was no clear evidence pointing to non-local individuals from distinctly different geological areas. In examining the dataset as a whole, one of the striking features is that three of four individuals with a slightly divergent Sr isotope pattern compared to the remainder of the population (Rein 3, 23 and 30), had genetic profiles that did not have any parallels in the collective burial (Fig 5), despite the fact that numerous haplotypes were represented at the site. The lack of evidence regarding to biological relations between these three individuals and the other members of the community supports the assumption that they had come from elsewhere. The isotope ratios in the late mineralising tooth crowns of Rein 3 and Rein 30 were equivalent to the local baseline, which means that if there had been a change of location it would have occurred during their childhood.
The individuals with identical haplotypes, on the other hand, largely yielded Sr isotope ratios typical of the locality and very slight differences between the 87Sr/86Sr ratios in the earlier and later mineralised teeth, which both argues against residential changes. This further confirms the notion of a closely related community, several generations of which sourced their food from geologically uniform economic areas. One possible exception was the 16–20 year-old male individual Rein 5 (hg T2b) with 87Sr/86Sr ratios typical of the locality in the first permanent molar, whilst the second permanent molar yielded the lowest Sr isotope ratio of the entire dataset. Whilst it is quite possible that the individual was born locally, he may have lived with a different community for some of his childhood. However, due to the only slight difference in the 87Sr/86Sr ratios of both teeth and the small number of comparable samples it must be emphasised that this is only one possible interpretation. The Sr isotope ratios overall paint a picture of the funerary site at Alto de Reinoso as a largely closed community within which (long-distance) mobility probably played a rather insignificant role. Because the sample only included a small number of females, no definitive conclusions could be drawn with regard to the rules of residence (matri- or patrilocality).
Alto de Reinoso in view of the Western and Central European Neolithic population and settlement history
It is widely accepted that the Neolithic spread to the Iberian Peninsula from the Near East via the Mediterranean route . Dependent on the Mesolithic background, there is, however, considerable variation with regard to the acceptance of pioneer farmers and the assimilation of autochthonous hunter-gatherers in different areas of Spain and Portugal. Because Western and Central Europe differ both in the beginning of the Neolithic and in the admixture of the first farmers, it is essential to analyze sites from various regions . General statements concerning the settlement and population history of Europe can only be reached in comparison.
When compared to other Neolithic groups on the Iberian Peninsula, the haplogroup composition of the Alto de Reinoso site showed close similarity to the Early and Middle Neolithic of Central Europe. (S11 and S12 Tables). Most interestingly, Middle to Late Neolithic groups in north Spain (NS_MLN) located in the Upper Ebro and Ambrona Valley [76,86] and Alto de Reinoso—though situated in close geographical proximity—showed essential differences at the mitochondrial level (see PCA, Fig 4). The haplogroup composition in north Spain (NS_EN and NS_MLN) appears to have remained relatively unchanged over the course of the Neolithic and a great many haplogroups typical of south-western hunter-gatherer populations (H, U, U5b) still existed in the area. On the other hand, a much greater number of “Neolithic haplogroups” (e.g. J, K, T2, and X) were present in the Early Neolithic north-east Spain, at Alto de Reinoso, and at the nearby Chalcolithic El Mirador site  arguing for an influx of Near Eastern Neolithic lineages from north-eastern Spain into the northern Meseta during the Early and Middle Neolithic. The results point to different population-dynamic processes in these regions in the latter stages of the Neolithic period. This stands in contrast to genetic data obtained in Central Europe, where populations seem to have remained rather uniform in their mitochondrial composition during the Early and Middle Neolithic [16,21,86].
The subsistence and dietary data confirm the impression of a closely related community
Overall, the stable carbon and nitrogen values show that all the members of the community consumed a very similar diet, again confirming the impression of a largely homogenous community. The evidence supporting this assumption is low standard deviations of 0.3 ‰ in δ13C and 0.4 ‰ in δ15N . The dietary evidence points exclusively to C3 plants and a slightly higher consumption of meat and dairy products or increased fertilisation of the farmland compared to contemporaneous cultural groups in Central Europe [131–133]. Increased fertilisation is less likely to occur in the keeping of smaller ruminants due to their usually higher mobility than in cattle and pig farming. Whilst at that time cattle and pigs were predominant in Central Europe, sheep and goats were the preferred sources of food in south-western Europe, where pigs and cattle were rarely part of the diet up to the Late Neolithic [105–106].
Furthermore there were no relevant differences in the consumption of food between the sexes or the different age groups within the community. Due to the lack of grave goods, which also underlines the impression of an egalitarian group, any social differentiation between the individuals could be excluded. Though not statistically verified, the females at Alto de Reinoso tended to yield lower nitrogen isotope values than the males (with values more than 0.2 ‰ lower being outside of the range which could be explained by analytical error; Fig 6). Moreover, the data of seven individuals from the top layer were somewhat set apart from the main data cluster due to more negative carbon isotope values (Rein 28, 28/2, 29–33). The lower carbon values could have been the result of different carbon sources in the diet, for example due to the consumption of animals and plants from wooded areas with a distinct canopy [134–135].
It is worth noting that all these individuals were from the final phase of use of the megalithic tomb. Besides climate change, which is difficult to identify in short-term variation, other possible reasons for the lower carbon values could have been changes in the subsistence economy such as an increase in legumes whose cultivation requires rather more shaded locations. The data of two individuals (Rein 29 and 32) differ from those of the others by lower carbon and higher nitrogen isotope values. The latter could have been as a result of an increased consumption of animal protein including freshwater fish .
Compared to most other sites in Spain, Portugal, and the South of France for which nutritional data are available the nitrogen isotope values at Alto de Reinoso were higher (Fig 6). The only site showing higher δ15N value was the Mesolithic site at El Collado near Valencia. However, this site was located on the Spanish Mediterranean coast and the high nitrogen isotope values were probably a result of the consumption of marine fish . The values of Alto de Reinoso are better comparable to late Neolithic data from the Garonne area in the South of France . In this case researchers interpreted the relatively high nitrogen isotope values as the result of a high consumption of animal protein or as a consequence of fertilising practices, and this may also have applied to Alto de Reinoso. Another reason for increased nitrogen values at Alto de Reinoso could have been the aridity in the Meseta compared to other regions on the Iberian Peninsula and in the South of France . The carbon isotope values measured at Alto de Reinoso were similar to those from the Late Neolithic site of Cova de la Pastora  in Spain, from Neolithic sites in Portugal (Casa da Moura, Gruta do Caldeirao, Gruta da Feteira, Gruta da Fontainhas, Roche Forte II, Cerca do Zambujal, Melides; ) and Late Neolithic sites in the Languedoc area of France .
Methodologically comprehensive and impressively supported by a multitude of results, the study is the first of its kind to provide a more in-depth understanding of the living environment of a burial community such as this. The bioarchaeological data generated have thus allowed us to go beyond the most important individual data such as age and sex and to trace biographical details including familial kinship relations, to reconstruct dietary habits and to highlight patterns of mobility and migration within the region and beyond. Some 50 deceased found their final resting place in the tomb. Based on radiocarbon dates we may state that they covered two or three generations and some of them exhibited close genetic similarities. As shown by the mtDNA profiles the community was composed of a series of independent familial groups whose members generally grew up in the vicinity of Alto de Reinoso. The demographic and genetic composition of the burial community and the reconstruction of its lifestyles and dietary habits, supported by the funerary rites, paint a picture of a close-knit local community whose subsistence was based on farming. With regard to their origins and ancestry the farmers buried at Alto de Reinoso exhibited mainly early farming genetic component, introduced during the Neolithisation of the Iberian Peninsula by immigrants from the Near East, as it is typically seen on the north-eastern part of Spain and in Central Europe. However, a smaller proportion of hunter-gatherer genetic legacy is also reflected in the maternal genetic composition of the population buried at Alto de Reinoso.
The number of bioarchaeological analyses has increased considerably in recent years [141–142], although it is still quite rare to find a study that combines a multitude of available scientific methods [20,22–23]. The often inadequate social-historical interpretation of such data must, however, be optimised at a meta-level, which will require an even more intensive discourse with archaeological research. For instance the fact that, as today, the societal conditions and social relationships were not necessarily always stable in the past, is often not taken into account. Therefore, we should generally assume that biographies included times of crisis, upheaval and social dislocation and both individuals and communities as a whole were continuously forced to adapt to changing social circumstances. From the point of view of the data generated, however, no particularly unusual individual situations or substantial group-dynamic breaks could be identified that would have impacted on the population of Alto de Reinoso.
We thank Willi Dindorf and Sabine Moeller-Rieker (University of Mainz), Sigrid Klaus, Michael Brauns and Bernd Höppner (Curt-Engelhorn-Centre Archaeometry gGmbH Mannheim) as well as Bernd Kromer, Susanne Lindauer, Robin van Gyseghem and Ute Blach (Klaus-Tschira Archaeometry Centre at the Curt-Engelhorn-Centre Archaeometry Mannheim) for technical assistance for the preparation of samples and for conducting measurements, Sandy Haemmerle (prehistrans) for language editing, and Sandra Pichler (Integrative Prehistory and Archaeological Science, Basel, Switzerland) and the anonymous reviewers for comments that greatly improved the manuscript.
Conceived and designed the experiments: KWA MARG RGP. Performed the experiments: KWA CK PH CR DN. Analyzed the data: KWA ASN CR PH CK SZ MARG RGP HAM. Contributed reagents/materials/analysis tools: KWA MARG RGP CTR IGML HAM. Wrote the paper: KWA CK MARG RGP PH ASN CR SZ. Participated in the excavation: MARG CTR IGML. Provided ancient samples, contextual information, and access to unpublished archaeological data: MARG RGP CTR IGML HAM. Provided input on the manuscript: CTR IGML HAM.
- Gibaja-Bao JF. Neolithic societies of Northeast Iberia (Catalonia): burials, grave goods and lithic tools. Curr Anthropol. 2004;45: 679–684. doi: 10.1086/425634
- Zapata L, Peña-Chocarro L, Pérez-Jordá G, Stika H-P. Early Neolithic agriculture in the Iberian Peninsula. Journal of World Prehistory. 2004;18: 283–325. doi: 10.1007/s10963-004-5621-4
- Rojo-Guerra MA, Kunst M, Garrido-Pena R, García-Martínez de Lagrán I, Morán Dauchez G. Paisajes de la Memoria. Asentamientos del Neolítico antiguo en el Valle de Ambrona (Soria, España). Valladolid: Servicio de Publicaciones de la Universidad de Valladolid; 2008.
- Pérez-Obiol R, Jalut G, Julià R, Pèlachs A, Iriarte MJ, Otto T, et al. Mid-Holocene vegetation and climatic history of the Iberian Peninsula. Holocene. 2011;21: 75–93. doi: 10.1177/0959683610384161
- Rojo-Guerra MA, Garrido-Pena R. From pits to megaliths: Neolithic burials in the interior of Iberia. In: Gibaja JF, Carvalho AF, Chambon P, editors. Funerary practices in the Iberian Peninsula from the Mesolithic to the Chalcolithic. BAR International Series 2417. Oxford: Archaeopress; 2012. pp. 21–28.
- Rojo-Guerra MA, Kunst M, Garrido-Pena R, García-Martínez de Lagrán I. El Neolítico en la Península Ibérica y su contexto europeo. Madrid: Editorial Cátedra; 2012.
- Darvill T, Malone C. Megaliths from antiquity. Cambridge: Antiquity publications; 2003.
- Bueno P, Barroso R, Balbín R, Campo M, Etxebarria F, Gonzalez A, et al. Áreas habitacionales y funerarias en el Neolítico de la Cuenca Interior del Tajo: la provincia de Toledo. Trabajos de Prehistoria. 2002;59: 65–79. doi: 10.3989/tp.2002.v59.i2.198
- Scarre C, Arias P, Burenhult G, Fano MA, Oosterbeek L, Schulting R, et al. Megalithic chronologies. In: Burenhult G, Westergaard S, editors. Stones and bones: formal disposal of dead in Atlantic Europe during the Mesolithic-Neolithic interface 6000–3000 BC. BAR International Series 1201. Oxford: Archaeopress; 2003. pp. 65–111.
- Steckel RH, Rose JC. The backbone of history: health and nutrition in the Western Hemisphere. Cambridge: Cambridge University Press; 2002.
- Buikstra JE, Beck LA. Bioarchaeology: the contextual analysis of human remains. Amsterdam: Academic Press; 2006.
- Ortner DJ. Identification of pathological conditions in human skeletal remains. 2nd ed. San Diego: Academic Press; 2003.
- Roberts C, Manchester K. The archaeology of disease. 3. ed. Ithaca: Cornell University Press; 2005.
- White TD, Folkens PA. Human osteology. 2nd ed. San Diego: Academic Press; 2000.
- Haak W, Brandt G, de Jong HN, Meyer C, Ganslmeier R, Heyd V, et al. Ancient DNA, strontium isotopes, and osteological analyses shed light on social and kinship organization of the Later Stone Age. Proc Natl Acad Sci U S A. 2008;105: 18226–18231. doi: 10.1073/pnas.0807592105. pmid:19015520
- Brandt G, Haak W, Adler C J, Roth C, Szécsényi-Nagy A, Karimnia S, et al. Ancient DNA reveals key stages in the formation of Central European mitochondrial genetic diversity. Science. 2013;342: 257–261. doi: 10.1126/science.1241844. pmid:24115443
- Bentley RA. Strontium isotopes from the earth to the archaeological skeleton: a review. Journal of Archaeological Method and Theory. 2006;13: 135–187. doi: 10.1007/s10816-006-9009-x
- Reitsema LJ. Beyond diet reconstruction: stable isotope applications to human physiology, health, and nutrition. Am J Hum Biol. 2013;25: 445–456. doi: 10.1002/ajhb.22398. pmid:23784719
- Alt KW, Knipper C, Held P, Nicklisch N, Fecher M, Roth C, et al. Königin Editha—ein Indizienbeweis zur Identifikation einer historischen Persönlichkeit aus dem Magdeburger Dom. In: Meller H, Schenkluhn W, Schmuhl BEH, editors. Königin Editha und ihre Grablegen in Magdeburg. Archäologie in Sachsen-Anhalt, Sonderband 18. Halle: Landesamt für Denkmalpflege und Archäologie Sachsen-Anhalt. Landesmuseum für Vorgeschichte; 2012. pp. 105–156.
- Alt KW, Knipper C, Peters D, Müller W, Maurer AF, Kollig I, et al. Lombards on the move—an integrative study of the Migration Period cemetery at Szólád, Hungary. PLoS One. 2014;9: doi: 10.1371/journal.pone.0110793.
- Brandt G, Knipper C, Nicklisch N, Ganslmeier R, Klamm M, Alt KW. Settlement burials at the Karsdorf LBK site, Saxony-Anhalt, Germany—biological ties and residential mobility. In: Whittle A, Bickle P, editors. Early farmers. The view from archaeology und Science. Proc Br Acad. 198; 2014. pp. 95–114. doi: 10.5871/bacad/9780197265758.003.0006
- Knipper C, Meyer C, Jacobi F, Roth C, Fecher M, Stephan E, et al. Social differentiation and land use at an Early Iron Age “princely seat”: bioarchaeological investigations at the Glauberg (Germany). J Archaeol Sci. 2014;41: 818–835. doi: 10.1016/j.jas.2013.09.019
- Doppler T, editor. Spreitenbach-Moosweg (Aargau, Schweiz): ein Kollektivgrab um 2500 v. Chr. Antiqua 51. Basel: Archäologie Schweiz; 2012.
- Lee EJ, Renneberg R, Harder M, Krause-Kyora B, Rinne C, Müller J, et al. Collective burials among agro-pastoral societies in later Neolithic Germany: perspectives from ancient DNA. J Archaeol Sci. 2014;51: 174–180. doi: 10.1016/j.jas.2012.08.037
- Meyer C, Kranzbühler J, Drings S, Bramanti B, Nehlich O, Richards MP, et al. Die menschlichen Skelettfunde aus der neolithischen Totenhütte von Benzingerode. Anthropologische Untersuchungen an den Bestattungen eines Kollektivgrabs der Bernburger Kultur. In: Meller H, editor. Die Totenhütte von Benzingerode. Archäologie und Anthropologie. Archäologie in Sachsen-Anhalt, Sonderband 7. Halle: Landesamt für Denkmalpflege und Archäologie Sachsen-Anhalt. Landesmuseum für Vorgeschichte; 2008. pp. 107–160.
- Arias P. Neighbours but diverse: social change in north-west Iberia during the transition from Mesolithic to the Neolithic (5500–4000 ca. BC). In: Whittle A, Cummings V, editors. Going over: the Mesolithic-Neolithic transition in north-west Europe. Proc Br Acad. 144;2007: 53–71. doi: 10.5871/bacad/9780197264140.003.0004
- Oross K, Marton T. Neolithic burials of the Linearbandkeramik settlement at Balatonszárszó and their European context. Acta Archaeologica Academiae Scientiarum Hungariae. 2012;63: 257–299. doi: 10.1556/aarch.63.2012.2.1
- Delibes G, Rojo Guerra M. Reflexiones sobre el trasfondo cultural del polimorfismo megalitico en la Lora burgalesa. Archivo Español de Arqueología. 75;2002: 1–15. doi: 10.3989/aespa.2002.v75.126
- Rojo Guerra M, Kunst M, Garrido Pena R, García Martínez-de-Lagrán I, Morán Dauchez G. Un desafio a la eternidad. Tumbas monumentales del Valle de Ambrona. Junta de Castilla y León. Memorias, Arqueología en Castilla y León, 14. Conserjería de Cultura y Turismo; 2005.
- Bueno P, de Balbín R, Barroso R, Rojas JM, Villa R, Félix R, et al. Neolítico y Calcolítico en Huecas (Toledo): el túmulo de Castillejo. Campaña de 1998. Trabajos de Prehistoria. 1999;56: 141–160. doi: 10.3989/tp.1999.v56.i2.280
- Haak W, Balanovsky O, Sanchez JJ, Koshel S, Zaporozhchenko V, Adler CJ, et al. Ancient DNA from European early Neolithic farmers reveals their Near Eastern affinities. PLoS Biol. 2010;8: doi: 10.1371/journal.pbio.1000536.
- Knipper C, Maurer A-F, Peters D, Meyer C, Brauns M, Galer SJG, et al. Mobility in Thuringia or mobile Thuringians: a strontium isotope study from Early Medieval Central Germany. In: Kaiser E, Burger J, Schier W, editors. Population dynamics in prehistory and early history. New approaches using stable isotopes and genetic. Berlin, Boston: de Gruyter; 2012. pp. 287–310.
- Marshall F, Pilgram T. NISP vs. MNI in quantification of body-part representation. Am Antiq. 1993;58: 261–269. doi: 10.2307/281968
- Ferembach D, Schwidetzky I, Stloukal M. Recommendations for age and sex diagnoses of skeletons. J Hum Evol. 1980;9: 517–549.
- Buikstra JE, Ubelaker DH. Standards for data collection from human skeletal remains. Fayetteville: Arkansas Archaeological Survey Series 44, 1994.
- Miles AEW. The dentition in the assessment of individual age in skeletal material. In: Brothwell DR, editor. Dental Anthropology. Symposia of the Society for the Study of Human Biology 5. Oxford: Symposium Publications Division, Pergamon Press; 1963. pp. 191–209.
- Stloukal M, Hanáková H. Die Länge der Längsknochen altslawischer Bevölkerungen—Unter besonderer Berücksichtigung von Wachstumsfragen. Homo. 1978;29: 53–69.
- 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 doi: 10.1002/ajpa.1330680106
- Lovejoy CO. Dental wear in the Libben population: its functional pattern and role in the determination of adult skeletal age at death. Am J Phys Anthropol. 1985;68: 47–56. pmid:4061601 doi: 10.1002/ajpa.1330680105
- Ubelaker DH. Human skeletal remains: excavation, analysis, interpretation. 2nd ed. Washington, DC: Taraxacum; 1989.
- Albert AM, Maples WR. Stages of epiphyseal union for thoracal and lumbar vertebral centra as a method of age determination for teenage and young adult skeletons. J Forensic Sci. 1995;40: 623–33. pmid:7595300
- Martin R. Lehrbuch der Anthropologie in systematischer Darstellung mit besonderer Berücksichtigung der anthropologischen Methoden. 2nd ed. Jena: Fischer; 1928.
- Lewis ME. The bioarchaeology of children. Perspectives from biological and forensic anthropology. Cambridge: Cambridge University Press; 2007.
- Waterman A, Thomas JT. When the bough breaks: childhood mortality and burial practice in Late Neolithic Atlantic Europe. Oxford Journal of Archaeology. 2011; 30 (2): doi: 10.1111/j.1468-0092.2011.00363.x.
- Rösing FW. Körperhöhenrekonstruktion aus Skelettmaßen. In: Knußmann R, editor. Anthropologie. Handbuch der vergleichenden Biologie des Menschen Vol I, 1. Stuttgart: Fischer; 1988. pp. 586–599.
- Formicola V, Franceschi M. Regression equations for estimating stature from long bones of early Holocene European samples. Am J Phys Anthropol. 1996;100: 83–88. pmid:8859956 doi: 10.1002/(sici)1096-8644(199605)100:1<83::aid-ajpa8>3.3.co;2-0
- Pearson K. Mathematical contributions to the theory of evolution. V. On the reconstruction of the stature of prehistoric races. Philos Trans R Soc Lond A. 1899;192: 169–245. doi: 10.1098/rsta.1899.0004
- Trotter M. Estimation of stature from intact long limb bones. In: Stewart TD, editor. Personal identification in mass disasters. Washington, DC: National Museum of Natural History, Smithsonian Institution; 1970, pp. 71–83.
- Sjøvold T. Estimation of stature from long bones utilizing the line of organic correlation. Hum Evol. 1990;5: 431–447. doi: 10.1007/bf02435593
- Breitinger E. Zur Berechnung der Körperhöhe aus den langen Gliedmaßenknochen. Anthropol Anz. 1938;14: 249–274.
- Bach H. Zur Berechnung der Körperhöhe aus den langen Gliedmaßenknochen weiblicher Skelette. Anthropol Anz. 1965;29: 12–21.
- Eshed V, Gopher A, Hershkovitz I. Tooth wear and dental pathology at the advent of agriculture: new evidence from the Levant. Am J Phys Anthropol. 2006;130: 145–159. pmid:16353225 doi: 10.1002/ajpa.20362
- Nicklisch N. Paläodemografische und epidemiologische Untersuchungen an neolithischen und frühbronzezeitlichen Bestattungen aus dem Mittelelbe-Saale-Gebiet im Kontext populationsdynamischer Prozesse. Ph.D. Thesis, The University of Mainz. 2014.
- Papathanasiou A. Health status of the Neolithic population of Alepotrypa Cave, Greece. Am J Phys Anthropol. 2005;126: 377–390. pmid:15386282 doi: 10.1002/ajpa.20140
- Schulting R, Fibiger L. Sticks, stones, and broken bones: Neolithic violence in a European perspective. Oxford: Oxford University Press; 2012.
- Molnar E. Tracing prehistoric activities: Musculoskeletal stress marker analysis of a stone-age population on the Island of Gotland in the Baltic sea. Am J Phys Anthropol. 2006;129: 12–23. pmid:16161142 doi: 10.1002/ajpa.20234
- Eshed V, Gopher A, Pinhasi R, Hershkovitz I. Paleopathology and the origin of agriculture in the Levant. Am J Phys Anthropol. 2010;143: 121–133. doi: 10.1002/ajpa.21301. pmid:20564538
- Caselitz P. Caries—ancient plague of humankind. In: Alt KW, Rösing FW, Teschler-Nicola M, editors. Dental anthropology. Fundamentals, limits, and prospects. Wien, New York: Springer; 1998. pp. 203–226.
- Meiklejohn C, Schentag C, Venema A, Key P. Socioeconomic change and patterns of pathology and variation in the Mesolithic and Neolithic of western Europe: some suggestions. In: Cohen MN, Armelagos GJ, editors. Paleopathology at the origins of agriculture. New York: Academic Press; 1984. pp. 75–100.
- Meyer C, Alt KW. Anthropologische Untersuchung der menschlichen Skelettfunde aus dem endneolithischen Kollektivgrab von Spreitenbach: osteologischer Individualbefund und Populationscharakteristik. In: Doppler T, editor. Spreitenbach-Moosweg (Aargau, Schweiz): ein Kollektivgrab um 2500 v. Chr. Antiqua 51. Basel: Archäologie Schweiz; 2012. pp. 104–157.
- Fernández Crespo T. Final Neolithic multiple burials in the Upper Ebro Valley: the case of San Juan Ante Portam Latinam (Basque Country, Spain). EAA Summer School eBook 1. 2007; pp.55–63. Available: http://eaa.elte.hu/FERNANDEZ.pdf.
- Subirà ME, López-Onaindia D, Yll R. Cultural changes in funeral rites during the Neolithic in the northeast in the Iberian Peninsula? The Cave of Pantà de Foix (Barcelona). Int J Osteoarchaeol. 2014; doi: 10.1002/oa.2400.
- Carlie-Thiele P. Spuren von Mangelerkrankungen an steinzeitlichen Kinderskeletten. Göttingen: Goltze; 1996.
- Bach A. Neolithische Populationen im Mittelelbe-Saale-Gebiet: Zur Anthropologie des Neolithikums unter besonderer Berücksichtigung der Bandkeramiker. Weimarer Monographien zur Ur- und Frühgeschichte 1. Weimar: Museum für Ur- und Frühgeschichte Thüringens; 1978.
- Larsen CS. Bioarchaeology. Interpreting behavior from the human skeleton. Cambridge: Cambridge University Press; 1997.
- Valaitis R, Hesch R, Passarelli C, Sheehan D, Sinton J. A systematic review of the relationship between breastfeeding and early childhood caries. Can J Public Health. 2000;91: 411–417. pmid:11200729
- Hillson S. Dental pathology. In: Katzenberg MA, Saunders SR, editors. Biological anthropology of the human skeleton. New Jersey: Wiley & Sons; 2008. pp. 301–340.
- Waterman AJ. Marked in life and death: identifying biological markers of social differentiation in Late Prehistoric Portugal. Ph. D. Thesis, The University of Iowa. 2012.
- Martínez-Cruz B, Harmant C, Platt DE, Haak W, Manry J, Ramos-Luis E, et al. Evidence of pre-Roman tribal genetic structure in Basques from uniparentally inherited markers. Mol Biol Evol. 2012;29: 2211–2222. pmid:22411853 doi: 10.1093/molbev/mss091
- Chandler HC. Using ancient DNA to link culture and biology in human populations. Ph.D. Thesis. The University of Oxford. 2003.
- Chandler H, Sykes B, Zilhão J. Using ancient DNA to examine genetic continuity at the Mesolithic-Neolithic transition in Portugal. In: Arias P, Ontañón R, García-Moncó C, editors. Actas del III Congreso del Neolítico en la Península Ibérica, Santander, Monografías del Instituto internacional de Investigaciones Prehistóricas de Cantabria 1. Santander: Servicio de Publicaciones, Universidad de Cantabria; 2005. pp. 781–786.
- de-la-Rua C, Izagirre N, Alonso S, Hervella M. Ancient DNA in the Cantabrian fringe populations: a mtDNA study from Prehistory to Late Antiquity. Quat Int. 2015;364: 306–311. doi: 10.1016/j.quaint.2015.01.035
- Gamba C, Fernández E, Tirado M, Deguilloux MF, Pemonge MH, Utrilla P, et al. Ancient DNA from an Early Neolithic Iberian population supports a pioneer colonization by first farmers. Mol Ecol. 2011;21: 45–56. doi: 10.1111/j.1365-294X.2011.05361.x. pmid:22117930
- Gómez-Sánchez D, Olalde I, Pierini F, Matas-Lalueza L, Gigli E, Lari M, et al. Mitochondrial DNA from El Mirador cave (Atapuerca, Spain) reveals the heterogeneity of Chalcolithic populations. PLoS One. 2014;9: doi: 10.1371/journal.pone.0105105.
- Hervella M, Izagirre N, Alonso S, Fregel R, Alonso A, Cabrera VM, et al. Ancient DNA from hunter-gatherer and farmer groups from Northern Spain supports a random dispersion model for the Neolithic expansion into Europe. PloS One. 2012;7: doi: 10.1371/journal.pone.0034417.
- Lacan M, Keyser C, Ricaut F-X, Brucato N, Duranthon F, Guilaine J, et al. Ancient DNA reveals male diffusion through the Neolithic Mediterranean route. Proc Natl Acad Sci U S A. 2011a;108: 9788–9791. doi: 10.1073/pnas.1100723108
- Lacan M, Keyser C, Ricaut F-X, Brucato N, Tarrús J, Bosch A, et al. Ancient DNA suggests the leading role played by men in the Neolithic dissemination. Proc Natl Acad Sci U S A., 2011b;108: 18255–18259. doi: 10.1073/pnas.1113061108
- Sampietro ML, Lao O, Caramelli D, Lari M, Pou R, Martí M, et al. Palaeogenetic evidence supports a dual model of Neolithic spreading into Europe. Proc R Soc B. 2007;274: 2161–2167. pmid:17609193 doi: 10.1098/rspb.2007.0465
- Sánchez-Quinto F, Schroeder H, Ramirez O, Avila-Arcos MC, Pybus M, Olalde I, et al. Genomic affinities of two 7,000-year-old Iberian hunter-gatherers. Current Biology. 2012;22: 1494–99. doi: 10.1016/j.cub.2012.06.005. pmid:22748318
- Bollongino R, Nehlich O, Richards MP, Orschiedt J, Thomas MG, Sell C, et al. 2000 years of parallel societies in Stone Age Central Europe. Science. 2013;342: 479–481. doi: 10.1126/science.1245049. pmid:24114781
- Bramanti B, Thomas MG, Haak W, Unterländer M, Jores P, Tambets K,et al. Genetic discontinuity between local hunter-gatherers and central Europe’s first farmers. Science. 2009;326: 137–140. doi: 10.1126/science.1176869. pmid:19729620
- Fu Q, Mittnik A, Johnson PLF, Bos K, Lari M, Bollongino R, et al. A revised timescale for human evolution based on ancient mitochondrial genomes. Current Biology. 2013;23: 553–559. doi: 10.1016/j.cub.2013.02.044. pmid:23523248
- Gamba C, Jones ER, Teasdale MD, McLaughlin RL, Gonzalez-Fortes G, Mattiangeli V, et al. Genome flux and stasis in a five millennium transect of European prehistory. Nat Commun. 2014;5: doi: 10.1038/ncomms6257.
- Haak W, Forster P, Bramanti B, Matsumura S, Brandt G, Tänzer M, et al. Ancient DNA from the first European farmers in 7500-year-old Neolithic sites. Science. 2005;310: 1016–1018. pmid:16284177
- Haak W, Lazaridis I, Patterson N, Rohland N, Mallick S, Llamas B, et al. Massive migration from the steppe was a source for Indo-European languages in Europe. Nature. 2015; 522: doi: 10.1038/nature14317.
- Lazaridis I, Patterson N, Mittnik A, Renaud G, Mallick S, Kirsanow K, et al. Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature. 2014;513: 409–413. doi: 10.1038/nature13673. pmid:25230663
- Lee EJ, Krause-Kyora B, Rinne C, Schütt R, Harder M, Müller J, et al. Ancient DNA insights from the Middle Neolithic in Germany. Archaeol Anthropol Sci. 2013;6: 199–204. doi: 10.1007/s12520-013-0173-0
- Szécsényi-Nagy A, Brandt G, Haak W, Keerl V, Jakucs J, Möller-Rieker S, et al. Tracing the genetic origin of Europe’s first farmers reveals insights into their social organization. Proc R Soc B. 2015;282: doi: 10.1098/rspb.2015.0339.
- Schönenberger R, Neugebauer J. Einführung in die Geologie Europas. Freiburg im Breisgau: Rombach; 1997.
- Prevedorou EA, Díaz-Zorita Bonilla M, Romero A, de Miguel Ibañez MP, Buikstra J E, Knudson K J. Residential mobility and dental decoration in Early Medieval Spain: results from eighth centuary Site of Plaza del Castillo, Pamplona. Dental Anthropology. 2010;23: 42–52.
- Ortega LA, Guede I, Cruz Zuluaga M, Alonso-Olazabal A, Murelaga X, Niso J, et al. Strontium isotopes of humain remains from the San Martín de Dulantzi graveyard (Alegría-Dulantzi, Álava) and population mobility in the Early Middle Ages. Quat Int. 2013;303: 54–63. doi: 10.1016/j.quaint.2013.02.008
- Waterman AJ, Peate DW, Silva AM, Thomas JT. In search of homelands: using strontium isotopes to identify biological markers of mobility in late prehistoric Portugal. J Archaeol Sci. 2014;42: 119–127. doi: 10.1016/j.jas.2013.11.004
- DeNiro MJ, Postmortem preservation and alteration of in vivo bone collagen isotope ratios in relation to palaeodietary reconstruction. Nature. 1985;317: 806–809. doi: 10.1038/317806a0
- Van Klinken GJ, Bone collagen quality indicators for palaeodietary and radiocarbon measurements. J Archaeol Sci. 1999;26: 687–695. doi: 10.1006/jasc.1998.0385
- Ambrose SH. Preparation and characterization of bone and tooth collagen for isotopic analysis. J Archaeol Sci. 1990;17: 431–451. doi: 10.1016/0305-4403(90)90007-r
- Jørkov MLS. Evaluating bone collagen extraction methods for stable isotope analysis in dietary studies. J Archaeol Sci. 2007;34: 1824–1829. doi: 10.1016/j.jas.2006.12.020
- Katzenberg MA. Stable isotope analysis: a tool for studying past diet, demography, and life history. In: Katzenberg AM, Saunders SR, editors. Biological anthropology of the human skeleton. New York, Chichester: Wiley-Liss; 2000. pp. 305–327.
- Schoeninger MJ, DeNiro MJ. Nitrogen and carbon isotopic composition of bone collagen from marine and terrestrial animals. Geochim Cosmochim Acta. 1984;48: 625–639. doi: 10.1016/0016-7037(84)90091-7
- DeNiro MJ, Epstein S. Influence of diet on the distribution of nitrogen isotopes in animals. Geochim Cosmochim Acta. 1981;45: 341–351. doi: 10.1016/0016-7037(81)90244-1
- Lee-Thorp JA. On isotopes and old bones. Archaeometry. 2008;50: 925–950. doi: 10.1111/j.1475-4754.2008.00441.x
- Minigawa M, Wada E. Stepwise enrichment of 15N along the food chains: further evidence and the relation between δ15N and animal age. Geochim Cosmochim Acta. 1984;48: 1135–1140. doi: 10.1016/0016-7037(84)90204-7
- van der Merwe NJ. Carbon isotope ecology of herbivores and carnivores. Palaeoecology of Africa and the surrounding islands. 1986;17: 123–131.
- Salazar-García DC, Aura JE, Olària CR, Talamo S, Morales JV, Richards MP. Isotopic evidence for the use of marine resources in the Eastern Iberian Mesolithic. J Archaeol Sci. 2014;42: 231–240. doi: 10.1016/j.jas.2013.11.006
- Castaños Ugarte P. Estudio arqueozoológico de los mesomamíferos del Neolítico de la cueva de Chaves (Huesca). Salduie. 2004;4: 125–171.
- Pérez Ripoll M. La explotación ganadera durente el III milénio a.C. en la Península Ibérica. In: Bernabeu Aubán J, Orozco Köhler T, editors. Actes del II Congrèso del Neolítico a la Península Ibèrica, València. PUV-Service de Publicacions de la Universitat de València, Saguntum-Extra 2; 1999. pp. 95–103.
- Häberle S, Fuller BT, Nehlich O, van Neer W, Schibler J, Hüster Plogmann H. Inter- and intraspecies variability in stable isotope ratio values of archaeological freshwater fish remains from Switzerland (11th–19th centuries AD). Environmental Archaeology; 2014; In press; doi: 10.1179/1749631414Y.0000000042.
- Robson HK, Andersen SH, Clarke L, Craig OE, Gron KJ, Jones AKG et al. Carbon and nitrogen stable isotope values in freshwater, brackish and marine fish bone collagen from Mesolithic and Neolithic sites in central and northern Europe. Environmental Archaeology; 2015; In press; doi: 10.1179/1749631415Y.0000000014.
- Aguilera M, Araus JL, Voltas J, Rodríguez-Ariza MO, Molina F, Rovira N, et al. Stable carbon and nitrogen isotopes and quality traits of fossil grains provide clues on substainability at the beginnings of Mediterranean agriculture. Rapid Commun Mass Spectrom. 2008;22: 1653–1663. doi: 10.1002/rcm.3501. pmid:18438779
- Araus JL, Febrero A, Buxó R, Rodríguez-Ariza MO, Molina F, Camalich MD, et al. Identification of ancient irrigation practices based on the carbon isotope discrimination of plant seeds: a case study from the South-East Iberian Peninsula. J Archaeol Sci. 1997;24: 729–740. doi: 10.1006/jasc.1997.0154
- Van Klinken GJ, Richards MP, Hedges REM. An overview of causes for stable isotopic variations in past European human populations: environmental, ecophysiological, and cultural effects. In: Ambrose SH, Katzenberg MA, editors. Biogeochemical approaches to paleodietary analysis. Advances in Archaeological and Museum Science 5. New York: Kluwer Academic publishers; 2000. pp. 39–63.
- Furholt M, Hinz M, Mischka D. As time goes by—meanings, memories and monuments. In: Furholt M, Hinz M, Mischka D, editors. As time goes by? Monumentality, landscapes and the temporal perspective. Proceedings of the international workshop “Socio-environmental dynamics over the last 12,000 years: the creation of landscapes II”, Kiel 2. Universitätsforschungen zur Prähistorischen Archäologie 206. Bonn: Verlag Dr. Rudolf Habelt; 2012. pp. 13–20.
- Hansen S, Müller J, editors. Sozialarchäologische Perspektiven: Gesellschaftlicher Wandel 5000–1500 v. Chr. Zwischen Atlantik und Kaukasus, Archäologie in Eurasien 24. Darmstadt: Verlag Philipp von Zabern; 2011.
- Lohrke B, Wiedmann B, Alt KW. Die anthropologische Bestimmung der menschlichen Skelettreste aus La Peña de la Abuela, Ambrona (Prov. Soria). In: Rojo-Guerra MA, Kunst M, editors. Sobre el significado del fuego en los rituales funerarios del Neolitico. Studia Archaeologica 91. Valladolid: Universidad de Valladolid; 2002: pp. 83–98.
- Nicklisch N, Uerpmann A, Alt KW. Estudio de los restos humanos neolíticos de la tumba de La Sima, Miño de Medinaceli. In: Rojo-Guerra MA, Kunst M, Garrido-Pena R, García-Martínez de Lagrán I, Dauchez GM, editors. Un desafio a la eternidad: tumbas monumentales del Valle de Ambrona. Arqueologia en Castilla y León 14. Valladolid: Junta de Castilla y León, Consejería de Educación y Cultura; 2005. pp. 323–339.
- Velasco Vázquez X. Estudio bioantropológio de los restos óseos humanos del túmulo de La Tarayuela. In: Rojo-Guerra MA, Kunst M, Garrido-Pena R, García-Martínez de Lagrán I, Dauchez GM, editors. Un desafio a la eternidad: tumbas monumentales del Valle de Ambrona. Arqueologia en Castilla y León 14. Valladolid: Junta de Castilla y León, Consejería de Educación y Cultura; 2005. pp. 340–361.
- Montero I, Rihuete C, Ruíz A. Precisiones sobre el enterramiento colectivo neolítico de Cerro Virtud (Cuevas de Almanzora, Almería). Trabajos de Prehistorica. 1999;56: 119–130. doi: 10.3989/tp.1999.v56.i1.294
- Solari Giachino A. Identificación de huellas de manipulación intencional en restos óseos humanos de origin arqueológico. Ph.D. Thesis, The University of Granada. 2010.
- Boaventura R, Ferreira MT, Neves MJ, Silva AM. Funerary practices and anthropology during Middle-Late Neolithic (4th and 3rd Millennia BCE) in Portugal: old bones, new insight. Anthropologie. 2014;12: 183–205.
- Donat P, Ullrich H. Einwohnerzahlen und Siedlungsgröße der Merowingerzeit. Zeitschrift für Archäologie. 1971;5: 234–265.
- Bocquet-Appel J-P. Recent advances in palaeodemography. Data, techniques, patterns. Dordrecht: Springer; 2008.
- Fernández-Crespo T, de-la-Rúa C. Demografic evidence of selective burial in megalithic graves of northern Spain. J Archaeol Sci. 2015;53: 604–617. doi: 10.1016/j.jas.2014.11.015
- Delibes de Castro G. Ritos funerarios, demografía y estructura social entre las comunidades neolíticas de la Submeseta Norte. In: Fábregas Valcarce R, Pérez Losada F, Fernández Ibáñez C, editors. Arqueoloxia da Morte na Peninsula Ibérica desde os Orixes ata o Medievo, Xinzo de Limia: Excmo. Concello de Xinzo de Limia. Serie cursos e congresos 3; 1995. pp. 61–94.
- Bienert HD. Skull cult in the Prehistoric Near East. Journal of Prehistoric Religion. 1991;5: 9–23.
- Bonogofsky M. Cranial modelling and Neolithic bone modification at `Ain Ghazal: new interpretations. Paléorient. 2001;27: 141–146. doi: 10.3406/paleo.2001.4734
- Meller H, Schunke T. Die Ahnen schützen den heiligen Ort—Belege für Kopf- und Schädelkult in Salzmünde. In: Meller H, editor. 3300 BC—Mysteriöse Steinzeittote und ihre Welt. Halle (Saale): Landesamt für Denkmalpflege und Archäologie Sachsen-Anhalt, Landesmuseum für Vorgeschichte; 2013. pp. 349–361.
- Veit U. Studien zum Problem der Siedlungsbestattungen im europäischen Neolithikum. Tübinger Studien zur Ur- und Frühgeschichtlichen Archäologie 1. Münster, New York: Waxmann; 1996.
- Zeeb-Lanz A. Besondere Schädel und mehr—der rätselhafte bandkeramische Ritualplatz von Herxheim (Pfalz). In: Wieczorek A, Rosendahl W, editors. Schädelkult. Kopf und Schädel in der Kulturgeschichte des Menschen. Regensburg: Curt-Engelhorn-Stiftung für die Reiss-Engelhorn-Museen, Schnell und Steiner; 2011. pp. 47–51.
- Zilhao J. From the Mesolithic to the Neolithic in the Iberian Peninsula. In: Price TD, editor. Europe`s first farmers. Cambridge University Press; 2000. pp. 144–182.
- Lovell NC, Nelson DE, Schwarcz HP. Carbon isotope ratios in palaeodiet: lack of age or sex effect. Archaeometry. 1986;28: 51–55. doi: 10.1111/j.1475-4754.1986.tb00373.x
- Knipper C, Fragata M, Brauns M, Alt KW. Isotopenanalysen an den menschlichen Skeletten aus dem endneolithischen Kollektivgrab von Spreitenbach: Studien zur Ernährung und Mobilität. In: Doppler T, editor. Spreitenbach-Moosweg (Aargau, Schweiz): ein Kollektivgrab um 2500 v. Chr. Antiqua 51. Basel: Archäologie Schweiz; 2012. pp. 188–219.
- Meyer C, Ramsthaler F, Brandt G, Knipper C, Schlenker B, Meller H, et al. Eine Neunfachbestattung gibt Rätsel auf. Archäologie in Deutschland. 2013;5: 30–33.
- Oelze VM, Siebert A, Nicklisch N, Meller H, Dresely V, Alt KW. Early Neolithic diet and animal husbandry: stable isotope evidence from three Linearbandkeramik (LBK) sites in Central Germany. J Archaeol Sci. 2011;38: 270–279. doi: 10.1016/j.jas.2010.08.027
- Drucker DG, Bridault A, Hobsen KA, Szuma E, Bocherens H. Can carbon-13 in large herbivores reflect the canopy effect in temperate and boreal ecosystems? Evidence from modern and ancient ungulates. Palaeogeogr Palaeoclimatol Palaeoecol. 2008;266: 69–82. doi: 10.1016/j.palaeo.2008.03.020
- Bogaard A, Fraser R, Heaton THE, Wallace M, Vaiglova P, Charles M, et al. Crop manuring and intensive land management by Europe’s first farmers. Proc Natl Acad Sci U S A 2013;110: 12589–12594. doi: 10.1073/pnas.1305918110. pmid:23858458
- Hedges REM, Reynard LM. Nitrogen isotopes and the trophic level of humans in archaeology. J Archaeol Sci. 2007;34: 1240–1251. doi: 10.1016/j.jas.2006.10.015
- Garcia-Guixé E, Richards MP, Subirà ME. Palaeodiets of humans and fauna at the Spanish Mesolithic site of El Collado. Curr Anthropol. 2006;47: 549–557. doi: 10.1086/504170
- Herrscher E, Le Bras-Goude G. Southern French Neolithic populations: isotopic evidence for regional specificities in environment and diet. Am J Phys Anthropol. 2010;141: 259–272. doi: 10.1002/ajpa.21141. pmid:19672846
- McClure SB, García O, de Togores CR, Culleton BJ, Kennett DJ. Osteological and paleodietary investigation of burials from Cova de la Pastora, Alicante, Spain. J Archaeol Sci. 2011;38: 420–428. doi: 10.1016/j.jas.2010.09.023
- Lubell D, Jackes M, Schwarcz H, Knyf M, Meiklejohn C. The Mesolitic-Neolithic transition in Portugal: isotopic den dental evidence of diet. J Archaeol Sci. 1994;21: 201–216. doi: 10.1006/jasc.1994.1022
- Alt KW. Bioarchäologie—die Methoden. In: Meller H, editor. 3300 BC—Mysteriöse Steinzeittote und ihre Welt. Halle (Saale): Landesamt für Denkmalpflege und Archäologie Sachsen-Anhalt, Landesmuseum für Vorgeschichte; 2013. pp. 38–41.
- Agarwal SC, Glencross BA. Social Bioarchaeology. West Sussex: Wiley & Sons; 2011.
- Salazar-García DC. Estudio de la dieta en la población neolítica de Costamar. Resultados preliminares de análisis de isótopos estables de C y N, Torre la Sal (Ribera de Cabanes, Castellón). Evolución del paisaje antrópico desde la prehistoria hasta el medioevo (E. Flors, coord.), Monografies de Prehistòria i Arqueologia Castellonenques 8; 2009. pp. 411–420.
- Doppler T, Gerling C, Heyd V, Knipper C, Kuhn T, Lehmann MF et al. Landscape opening and herding strategies: Carbon isotope analyses of herbivore bone collagen from the Neolithic and Bronze Age lakeshore site of Zürich-Mozartstrasse, Switzerland. Quaternary International; 2015; In press; doi: 10.1016/j.quaint.2015.09.007
- Fischer A, Olsen J, Richards M, Heinemeier J, Sveinbjornsdottir AE, Bennike P. Coast-inland mobility and diet in the Danish Mesolithic and Neolithic: evidence from stable isotope values of humans and dogs. J Archaeol Sci. 2007;34: 2125–2150. doi: 10.1016/j.jas.2007.02.028
- An C-B, Dong W, Li H, Zhang P, Zhao Y, Zhao X et al. Variability of the stable carbon isotope ratio in modern and archaeological millets: evidence from northern China. J Archaeol Sci. 2015;53: 316–322. doi: 10.1016/j.jas.2014.11.001
- Bronk Ramsey C, Lee S. Recent and planned developments of the program OxCal. Radiocarbon. 2013;55: 720–730. doi: 10.2458/azu_js_rc.55.16215
- Reimer PJ, Bard E, Bayliss A, Beck JW, Blackwell PG, Bronk Ramsey C et al. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. Radiocarbon. 2013;55: 1869–1887. doi: 10.2458/azu_js_rc.55.16947
- Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet 1999;23: 147. pmid:10508508 doi: 10.1038/13779