Introduction

The Linearbandkeramik or Linear Pottery culture (LBK) played a key role in the Neolithization of central Europe. Culturally, economically, and genetically, the LBK had its ultimate roots in western Anatolia, but it also displayed distinct features of autochthonous European Mesolithic hunter-gatherer societies. Several models for the origins of the LBK culture have been proposed over the years1.

The Indigenist model suggests the LBK was founded through the adaptation of elements of the West Asian Neolithic Package by indigenous Mesolithic populations exclusively through frontier contact and cultural diffusion. The Integrationist model views the formation of LBK as the integration of Mesolithic hunter-gatherers into an agro-pastoral lifeway through mechanisms such as leapfrog colonization, frontier mobility and contact. According to this model, small groups associated with the Starčevo-Körös-Criş (SKC) culture, the likely LBK predecessors in Europe, left their homelands in the Balkans (where most of their own ancestors had arrived earlier from Anatolia), and settled new areas to the northwest. Contacts with local Mesolithic groups and exchange of products would have resulted in the co-optation of hunter-gatherers into farming communities, where they would have adopted farming practices1. Evidence of such interactions exists at the Tiszaszőlős-Domaháza site in northeastern Hungary, containing interments of individuals of mostly hunter gatherer genetic ancestry buried in a clearly SKC context2,3.

The Migrationist model suggests that a sparsely populated territory of Mesolithic central Europe was taken over by pioneering agro-pastoral groups associated with the SKC culture, which gradually displaced indigenous hunting-gathering populations, who did not significantly influence the arriving Starčevo colonizers. According to this model, newcomers would have replicated their ancestral material culture in the newly settled territory without incorporating the material culture features of the local indigenous populations. Some variation, due to innovation and adaptation to the new environment and sources, would have involved changes in technology such as pottery and building material as well as lithic tool sources. At the same time, symbolic systems, such as decorative designs and cultural objects, would have remained unchanged. This model appeared at the end of the 1950s4 and gained wide support in the second half of 20th century5,6,7,8.

To date, ancient DNA (aDNA) studies have convincingly shown that Neolithic European farming populations were primarily genetic descendants of central and western Anatolian Neolithic farmers (ANFs)9,10,11. Their genetic signature is clearly distinct from autochthonous Mesolithic European hunter gatherers (HGs) of central Europe (WHGs) at the level of uniparental markers such as mitochondrial DNA (mtDNA) and Y chromosome as well as genome-wide. Nevertheless, the extent to which the newcomers interacted both culturally and genetically with local hunter gatherers remains unclear; that is, it remains unclear to what extent an Integrationist or Migrationist model is accurate. Genetically, Neolithic central European farmers carried a minor proportion of genetic ancestry characteristic to WHG populations, but the extent and the timing of the WHG admixture in the gene pool of the European Neolithic descendants of Anatolian farmers varies across central Europe3. While the amount of WHG ancestry in European Neolithic farmers had been observed to increase throughout the Neolithic in the present-day territories of Hungary, Germany and other regions of Europe3,9,12,13,14,15, the initial degree of exchange remains unresolved, in part due to a scarcity of human remains contemporaneous with the earliest stages of the Neolithic farming migration.

The Brunn 2 archaeological site, part of the Brunn am Gebirge, Wolfholz archaeological complex south of Vienna, Austria (Fig. 1), is the oldest Neolithic site known in Austria and one of the oldest in all of central Europe. It belongs to the earliest stage of the development of LBK, called the Formative phase. Radiocarbon dates obtained for Brunn 2 time the site to about 5670–5350 cal BCE16,17,18,19,20. The main characteristic of the settlements of the Formative phase is the absence of fine pottery and the use of coarse pottery with clear Starčevo features. The leading role of Anatolian migrants in the formation of cultural attributes of the earliest farmers of Europe is evident through the comparative typological analysis of material culture artifacts from the Brunn 2 site20. In addition to rich trove of culture artifacts, Brunn 2 yielded four human burials. The initial radiocarbon dating of the remains confirmed these to be contemporaneous with the earliest phase of the Brunn am Gebirge complex20 and, thus, to represent some of the earliest central European Neolithic farmers. We set out to perform a bioarchaeological analysis of these individuals to examine genetic ancestry as well as diet and mobility at the dawn of the European Neolithization, in an effort to refine the model of the establishment of farming in the Neolithic central Europe.

Figure 1
figure 1

The location of Brunn am Gebirge, Wolfholz site on the map of Europe. Map image is from https://en.wikipedia.org/wiki/File:Europe_blank_map.png (image is in the Public Domain).

Materials and Methods

Brunn 2 archaeological site

The Brunn am Gebirge - Wolfholz excavation site (48.120396, 16.291722) is located near Vienna, Austria. Six sites uncovered at the location have shown a development of the LBK from the Formative phase until the Musical Note Pottery (Notenkopfkeramik) phase. Sixteen radiocarbon dates obtained for Brunn 2 date the oldest part of this site to about 5,670–5,450 calBCE16,17,18,19,20, which places this site along with Szentgyörgyvölgy-Pityerdomb21 and Zalaegerszeg-Andráshida22 in Hungary within the Formative phase of LBK. The main characteristics of these settlements is the absence of fine pottery and the use of only coarse pottery with clear Starčevo features. In total, 3,620 pottery vessels have been recovered from Brunn am Gebirge along with 100,000 pottery sherds. The ceramic collection from Brunn 2 includes 11 narrow-neck vessels (“amphorae”) and fragments of seven clay figurines (“idols”). In addition, over 15,000 lithic artifacts were found at Brunn am Gebirge, most of them at Brunn 2, which is unusual for a central European Neolithic settlement. The presence of anthropomorphic forms, amphora, remains of musical instruments (clay flutes), as well as large quantities of lithic artifacts indicates that Brunn 2 was part of a “central settlement”, a site dedicated to the ritual activities of large LBK communities23.

Brunn 2 burials

Burial 1 was found in a clay extraction pit in the southern part of the Brunn am Gebirge site 2. Burials 2 and 4 were found in the long ditches of houses no longer in use at the time of burial. Burial 3 was not associated with an above-ground structure. All four individuals were buried in a typical LBK position on their left side and had a head orientation to the northeast facing east. The teeth from all four seem to have been heavily worn, suggestive of a plant-based diet20. One of the skeletons (Individual 2) was found with six trapezes made from radiolarite sourced some 200 km southeast at Bakony-Szentgál, near Lake Balaton in western Hungary19. A detailed description of the four individuals interred at Brunn 2 can be found in20.

Radiocarbon dating and stable isotope analysis

Stable isotope analysis of carbon and nitrogen is routinely applied in archaeological studies to investigate diet. Carbon isotope values (δ13C) measured in bone collagen and dentine allows us to distinguish between dietary proteins obtained from marine, terrestrial, and freshwater resources24,25. Nitrogen isotope values (δ15N) allow us to investigate the trophic level of an organism within a food web26,27. Dietary isotope studies of LBK groups performed to date have produced δ13C values of ca. 20‰ and δ15N values of ca. 10‰, which reflect dietary pathways that are reliant on C3 terrestrial proteins (e.g. terrestrial plants and fauna)28,29,30.

Bulk stable isotope values from bone collagen reflect a mixed signal of dietary proteins consumed in the last ca. 10 years of an individuals’ life31. In teeth, dentine is formed in layers, with each layer representing a discrete chronological period and containing a record of corresponding dietary patterns32,33. For this project, we analyze the dentine of second pre-molars, which represent a bulk measurement of diet from the age of approximately 3 to 15 years34,35 (and references therein).

The main limitation associated with this approach to sampling is the issue of nursing and weaning practices. A number of archaeological studies have utilized stable isotope analysis to examine infant diet and feeding practices36,37. Their conclusions are based on the observation that a child consuming its mother’s milk is effectively one trophic level above the mother. Previous stable isotope investigations into weaning processes at LBK sites suggests an average weaning age of approximately 3 years of age30,38. Therefore, it is a distinct possibility that the dentine samples analyzed in this study will include some of the nursing signal (e.g. 3‰ higher in δ15N than the adult values, and ca. 1‰ more enriched for δ13C).

Stable isotope measurements were obtained at BETA Analytic, Miami, FL (BETA) using a modified version of the Longin (1971) collagen extraction method39. δ13C measurements are reported on the VPDB scale and δ15N is reported with reference to AIR. A number of quality control parameters are in place, including a collagen yield of 1% at 5 mg or above, and a C:N ratio of between 2.9 and 3.640,41.

Radiocarbon dating was accomplished at BETA Analytic Accelerated Mass Spectrometry (AMS) Laboratory, Miami, FL (BETA) and the Penn State Radiocarbon 14C Laboratory, University Park, PA (PSUAMS).

Strontium isotopes from enamel analysis

It is possible to obtain specific clues about the movement of people in the past from the chemistry of prehistoric human teeth. The basic principle for the isotopic proveniencing of human remains essentially involves the comparison of isotope ratios in human tooth enamel with local, or baseline, levels from the place of burial. Tooth enamel is a remarkable repository of childhood environment. Tooth enamel forms in the first years of life and remains unchanged through life and often for a very long period after death. A variety of studies have demonstrated that enamel is highly resistant to post-mortem contamination42,43,44. Enamel is largely a mineral, hydroxyapatite, Ca10(PO4)6(OH)2, composed primarily of calcium and oxygen. A few other elements also can be deposited in the apatite. Strontium and lead, for example, substitute for calcium during mineral formation.

We lightly abrade the surface of the enamel to be sampled using a dental drill to remove surficial dirt and calculus and the outermost enamel due to the possibility of contamination by diffusion. After abrading the surface, we remove one or more small chips from the side of the molar or drill 5 to 10 milligrams of powder from the enamel. Any remaining dentine is removed. Samples of enamel weighing 3–5 mg are dissolved in 5-molar nitric acid for strontium analysis. The strontium fraction is purified using EiChrom Sr-Spec resin and eluted with nitric acid followed by water. Isotopic compositions are obtained on the strontium fraction using a VG (Micromass) Sector 54 thermal ionization mass spectrometer (TIMS). Strontium is placed on single Re filaments and analyzed using a quintuple-collector dynamic mode of data collection. Internal precision for 87Sr/86Sr analyses is typically 0.0006 to 0.0009 percent standard error, based on 100 dynamic cycles of data collection i.e., ±0.000006. The analysis was performed at the Geo Chemistry Labs, University of North Carolina at Chapel Hill, Chapel Hill, NC.

Ancient DNA procedures

Powder was obtained from teeth of all four individuals interred at Brunn 2 in a dedicated clean room. DNA was extracted45,46 using between 69 and 82 mg of powder, and a portion of the extract was converted into DNA sequencing libraries. For each of the 4 samples an individually barcoded and UDG treated library47 was built (L1). For one sample (Individual 1, S6912) additional libraries were prepared from the same DNA extract using 4 different protocols (S6912.E1.L2 non-UDG treated double-stranded, S6912.E1.L3 non-UDG treated single-stranded48, S6912.E1.L5 UDG treated double-stranded, S6912.E1.L6 UDG treated single-stranded). Libraries were then enriched for both the mitochondrial genome49 and about 1.2 million single nucleotide polymorphisms (SNPs)50, and sequenced on an Illumina NextSeq. 500 instrument. We followed a previously described bioinformatics procedure9, merging sequences overlapping by at least 15 base pairs, mapping to the mitochondrial genome reference sequence rsrs and to the human genome reference sequence hg19 using bwa (v.0.6.1)51, and removing duplicated sequences that mapped to the same start and stop locations and had the same molecular barcodes. For mitochondrial genome analysis, we built a consensus sequence52, and for nuclear genome analysis, we represented each targeted SNP by one randomly chosen sequence passing previously reported minimum mapping and base qualities9. We evaluated ancient DNA authenticity by tabulating characteristic C-to-T damage rates at the terminal nucleotides of sequencing reads and by measuring apparent heterozygosity rates in haploid genome regions (mitochondrial genome [all >95% matching consensus] and male X chromosomes [one sample with sufficient coverage, >98% consistency])53,54. Full technical information on the data we produced for each sample is given in Supplementary Table S1.

Whole-genome statistical analysis

Principal component analysis (PCA) was carried out using the smartpca software55. We computed axes using 1035 present-day individuals with West Eurasian ancestry genotyped at 593,124 SNPs on the Affymetrix Human Origins array56, and projected the newly reported and previously published ancient individuals2,3,9,10,11,13,57,58,59,60,61 using the least-squares option (‘lsqproject: YES’). We inferred patterns of shared ancestry with WHGs using f-statistics as previously described3. To test for symmetry of Individual #3 to different Neolithic populations, we used the statistic f4(#3, Mbuti; Neolithic_A, Neolithic_B), with Central African hunter-gatherers as an outgroup62.

Results

Genetic analyses

We obtained genetic data passing quality control for three out of the four individuals interred at Brunn 2. No usable genetic data was obtained from Individual 4. Individuals 1–3 were males by genetic typing. The mitochondrial lineages of Individuals 1–3 were J1, U5a1, and K1b1a (Table 1), while their Y chromosomal lineages were BT, CT, and G2a2a1a, respectively. For Individual 1, we note that we ostensibly observed derived alleles at the diagnostic haplogroup P sites CTS3446 and F212, the R1 site CTS997, and the R1b1a1a2 sites PF6444 and L749 (nomenclature from the International Society of Genetic Genealogy, http://www.isogg.org), but these were mostly carried on long sequencing reads (41, 96, 74, 131, and 96 bases, respectively), none of which had evidence of ancient DNA damage, so we believe some or all of them to be due to low levels of contamination. We also observed an ancestral allele at the haplogroup R site L1225 (read length 45, likewise not damaged).

Table 1 Genetic data for Brunn 2 individuals.

SNP data obtained from whole genome sequencing were used as the basis for a PCA plot comparing to Neolithic Anatolian and early Neolithic European individuals (Fig. 2). Individual 2 (I6913) fell closest to WHGs but shifted toward EEFs/ANFs, while Individuals 1 (I6912) and 3 (I6914) grouped with Anatolian Neolithic farmers and closely related central European farming groups. In the zoomed-in plot (Supplementary Fig. S1), I6914 appears to be borderline between ANFs and ENFs, while I6912 is on the high end of WHG relatedness for early European farmers. We note though that I6912 and especially I6913 have relatively low sequencing coverage, so their exact positions in PCA should be interpreted with caution.

Figure 2
figure 2

Individuals from Brunn am Gebirge site 2 plotted onto the space of principal components defined by 1035 present-day individuals with West Eurasian ancestry and including selected published ancient individuals of European and Anatolian hunter gatherer and Anatolian farming ancestries. Sources of the published data are given in the Materials and Methods. I6912, Individual 1; I6913, Individual 2; I6914, Individual 3. The PCA computation was done with smartpca, version 16690 and the visualization was made in R, version 3.5.1.

To formalize these observations, we used f-statistics to measure asymmetries in allele sharing (see Materials and Methods). For Individuals 1–3, we estimated genome-wide WHG-related ancestry proportions of 12 ± 3%, 57 ± 8%, and <1%, respectively, consistent with the PCA results (Table 1). For Individual 1, we determined his WHG ancestry to be more closely related to western and central European hunter-gatherers than to southeastern European hunter-gatherers (f4(#1, Anatolia N; W/C WHG, SE WHG) = −2.0*10−3, |Z| = 2.4; cf. ref. 3); for Individual 2, this distinction could not be made with precision (|Z| < 1). Finally, we tested for unequal relatedness of Individual 3 (who yielded the highest sequencing coverage) to various published Neolithic farmers and found that while he is approximately symmetrically related to Neolithic Anatolians and Starčevo-associated individuals from southeastern Europe, he does share excess alleles with other LBK groups from central Europe (|Z| > 3 for Germany, |Z| > 5 for Austria). We note that this signal is not caused by WHG-related ancestry in the LBK groups, as the statistic value for WHG has the opposite sign (f4(#3, Mbuti; WHG, Anatolia N) ≪ 0, Z < −17).

Stable isotope (δ13C and δ15N) data

Two of the bone collagen samples (Individuals 1 and 2) failed the quality control parameters, with C:N ratios outside the accepted range of 2.9–3.6 (see40), likely indicating collagen depletion (Table 2). Therefore, the data from those two bone samples should be considered unreliable. The δ13C and δ15N values for the bone collagen sample which passed quality control (Individual 3), are consistent with the isotope values observed in the dentine samples (see below). Ideally, a number of wild and domestic animal bone samples from the site would have been analyzed, to form an isotopic baseline to compare and interpret the human samples against27,31. In their absence, carbon and nitrogen isotope analyses of domestic fauna from the approximately contemporary, and environmentally similar LBK site of Vedrovice, in the Czech Republic63, will be utilized. The domestic cattle samples at Vedrovice produced δ13C values of −20.2‰ ± 0.3 and δ15N values of 6.2‰ ± 1, the sheep/goat samples produced values of −19.8‰ ± 0.2 and 5.9‰ ± 0.4 for δ13C and δ15N respectively, and finally the domestic pig samples −20.4‰ ± 0.4 (δ13C) and 8.2‰ ± 0.9 (δ15N). The faunal δ13C and δ15N values from Vedrovice reflect a C3 terrestrial environment, and if we use these values as an imperfect baseline to examine human diet at Brunn, we could tentatively suggest that cattle and sheep/goat (or similar herbivores) formed the mainstay of diet. If we used the δ15N values of pig samples from Vedrovice, we would suggest the values are too enriched to contribute significantly to human diet, as we would anticipate human δ15N values of 11‰ or above. However, as the fauna do not originate from the same site, it is difficult to produce a definite conclusion.

Table 2 Radiocarbon and stable isotope data of the individuals from Brunn am Gebirge site 2, as well as mean stable isotope values (±SD) for European inland Mesolithic (EIM) and European inland Neolithic (EIN)64, central European LBK30, and Anatolian Neolithic (AN)9. Radiocarbon dates for Individuals 2 and 4 (ETH-14827 and ETH-11150) are from20.

Despite the caveats associated with using bulk dentine values, the isotope values for the dentine samples of Individuals 1, 2 and 3 (Table 2) indicate dietary pathways reliant on C3 terrestrial proteins (e.g. C3 plant species and/or the herbivores which feed on them). It is not possible to elucidate whether a nursing signal is present in the isotopic values of the dentine.

When the δ13C and δ15N values from the Brunn 2 individuals are plotted against the average δ13C and δ15N values from European LBK sites30, European Inland Neolithic farming communities (EIN) and European Inland Mesolithic HGs (EIM)64, as well as Anatolian Neolithic farmers (AN)9, all three Brunn 2 individuals place within the range of Anatolian and European Neolithic farmers (Fig. 3). The values for Individual 2 from dentin and bone collagen overlap with the lower margin for δ15N values of EIM, while the δ15N values for Individuals 1 and 2 fall outside of the EIM margin and within the EIN margin (Fig. 3).

Figure 3
figure 3

Distribution of carbon and nitrogen isotope ratios of Individuals 1, 2 and 3 from Brunn am Gebirge site 2 in relation to the mean isotope ratios (±SD) from Anatolian Neolithic (AN), European Inland Mesolithic (EIM), European Inland Neolithic (EIN), and Linearbandkeramik (LBK) populations. Sources of the published data are given in Table 2. The figure and underlying statistical analyses were generated using Microsoft Excel 2019 for Mac, version 16.25.

Strontium isotope analysis

A small-scale strontium isotope analysis of two of the burials (Individuals 1 and 2) from Brunn 2 was undertaken, limited by the amount of enamel that had survived. Strontium isotope ratios (86Sr/87Sr) in human tooth enamel reflect the nutrients consumed during tissue formation in early childhood, assumed to represent the place of birth. Local or baseline ratios at the place of burial should reflect the place of death. Difference between enamel and baseline should identify individuals of non-local birth.

Baseline values were not directly available from the site of Brunn 2 but had been measured at a nearby Longobard cemetery less than 100 m from Brunn 2 (Peter Stadler, personal communication, 2019). Measurement of modern vegetation suggests a local range of 86Sr/87Sr values between 0.7085 and 0.7098.

The results of 86Sr/87Sr measurements are presented in Table 2. If we use the Longobard cemetery values in consideration of the two human burials from Brunn 2 then Individual 1 is local and Individual 2 is non-local, falling outside the local range based on the value from the nearby cemetery.

Discussion

While the remains of early Neolithic farmers in central Europe are relatively abundant, very few remains of contemporaneous hunter gatherers are known, making the understanding of the HGs’ lifeways and their integrations with incoming Anatolian-related farming migrants difficult, particularly during the earliest steps of the Neolithization of Europe. While recent genetic studies have pointed to a limited genetic exchange between immigrant farmers and local European HGs, the flow of material goods from farmers to HGs has been documented (see references in65). The reciprocal material goods flow into the farming communities has been more difficult to identify65.

A bioarchaeological analysis of the remains of the interred at Brunn 2 presented in this study allows insight into the life history of early European Neolithic farmers who lived near the beginning of the establishment of farming economies in central Europe, revealing evidence of biological interaction between incoming Anatolians and local HGs during the earliest stages of the arrival of farming in the Neolithic Central Europe.

The mtDNA lineages of Individuals 1 and 3 belong to two of the most common mtDNA subclades found in Neolithic individuals from the Near East as well as their Neolithic European descendants3,66,67. At the same time, an individual belonging to K1b (K1b2) has been identified in a Mesolithic forager from the Baltic68. On the other hand, divisions of haplogroup U5 such as the U5a1 lineage identified in Individual 2 are generally considered to be characteristic of European hunter-gatherers69,70. At the same time, a U5-carrying individual (U5b2) has recently been identified in the Çatalhöyük population of central Anatolia71.

Individual 3 carries Y chromosome haplogroup G2a2a1a, from the larger set of G2a Y chromosomal lineages, which are characteristic of ANF and ENF populations9,15,66. Earlier studies of the LBK remains from early Neolithic sites in the Carpathian Basin and southeastern Europe that utilized uniparental genetic markers (mtDNA and Y chromosome) identified the G2a Y-chromosomal lineage to be the prevailing Y haplogroup in the early European Neolithic farmers67. The contrasting high diversity of contemporaneous mtDNA lineages suggested a reduced male-specific lineage diversity in early LBK communities67. Subsequent studies such as61 using genome-wide methods found no evidence of sex bias in the LBK. Individuals 1 and 2 have Y-chromosomes from the macro-lineages BT and CT, respectively, but due to their low-coverage data, we are not able to assign them with greater precision.

The genetic signature of Individual 3 is that of Neolithic Anatolian-related ancestry, consistent with that of most of the representatives of European Neolithic farming cultures, including LBK and Starčevo. Our analyses indicate that this individual had very little ancestry derived from European hunter-gatherers (likely zero, and no more than 1%). At the same time, Individuals 1 and 2 had WHG ancestry that was acquired after their ancestors had left Anatolia. We were unable to determine dates for this admixture, so it is possible either that it occurred locally or that the Brunn 2 migrants encountered WHGs along their journey and integrated WHG ancestry into the their predominantly ANF-derived genetic pool prior to their arrival in central Europe. It is also possible that the Brunn 2 migrants interacted with, or descended from, the Anatolia-derived farming communities (monochrome white painted pottery groups) that settled in the Balkans ca. 600 years earlier and would have also had opportunities to incorporate WHG ancestry in their gene pool since leaving Anatolia. However, the very high WHG-related proportion in Individual 2, combined with the western European affinity of the HG-related ancestry in Individual 1, points toward recent post-arrival admixture in central Europe as the most likely scenario.

Lithic artifacts from Brunn 2 indicate active interaction between early ENF and local HG population groups in the Early Neolithic. A likely explanation of the presence of 15,000 lithic artifacts at the Brunn am Gebirge - Wolfholz site is that the Neolithic farmers produced hunting implements for trade with local HGs. The absence of the evidence of violence at early LBK settlements suggests low hostility between the local HGs and LBK farmers65. The material for lithic implements from the grave of Individual 2 was sourced from around Lake Balaton20 where other settlements of the Formative LBK phase have been found. Individual 2’s strontium isotope ratio measurements indicate he was not born at the Brunn 2 settlement site and could conceivably have come from the area where the lithic material had been procured, namely Bakony-Szentgál in Hungary. The sourcing of the lithic material from distal areas is not uncommon for LBK settlements in Austria72, but the number of lithic artifacts found at Brunn 2 may have significance. A lithic production center at Brunn 2 would have needed the knowledge of dedicated craftsmen, some of them could conceivably have come from the local HG communities. This might explain the presence of individuals such as Individual 2, with high proportions of HG-related ancestry, on a permanent or a semi-permanent basis within LBK settlements, and their subsequent integration into the LBK communities.

The diet isotope δ15N ratios from dentin ranged from 9.21 for Individual 2 to 9.78 for individual 1 to 10.35 for Individual 3. Individuals 1 and 2 are roughly contemporaneous, and their δ15N variation likely reflects individual dietary specifics. However, the δ15N value for Individual 3 is somewhat elevated compared to that of Individuals 1 and 2. The δ15N value for Individual 3 is in the upper range for δ15N variation for EIN and LBK, within the lower range of δ15N variation for EIM, and within the average for δ15N of ANF (Table 2, Fig. 3). Individual 3 is also chronologically the youngest of the three. As with Individuals 1 and 2, the δ15N measurements in the Individual 3 could reflect the specifics of individual dietary patterns. At the same time, the varying nitrogen isotope ratios could be a result of oscillating environmental conditions during the Formative phase of LBK leading to the failure of initially maladapted domesticated plants from semi-arid Anatolia to thrive in the continental climate of central Europe73 and causing early European farmers to periodically rely more on animal protein rather than agricultural crops. Another explanation for the δ15N variation could be the progressive expansion of domestic animal herds in early LBK leading to an increased availability of animal protein and/or increased crop manuring, also leading to elevated δ15N values73.

The symmetric genetic relationship of Individual 3 to Starčevo and ANF individuals studied to date and his greater genetic affinity to other LBK individuals implies that the individual was from a population that had experienced a small amount of genetic drift not shared with Anatolian and southeastern European farmers studied to date. In theory, the excess relatedness of Individual 3 to other LBK-associated individuals could be due to shared WHG ancestry, but (a) we find approximately zero such ancestry in Individual 3, and (b) direct allele-sharing tests show that such a signal would in fact be in the opposite direction. In any case, our results show that the lineage that gave rise to the primary ancestry of central European LBK-associated populations was represented at Brunn 2 together with other sampled early Neolithic sites.

It is clear that the process of formation of the Starčevo and Linear Pottery cultures was more complicated than a mere immigration into a new area and the subsequent cultural deterioration during the movement. It had to also include the influence of local populations (the Early Neolithic in Bulgaria, Serbia and Croatia, and the Mesolithic in Hungary and Austria), and the adaptation to new ecological conditions, as well as new sources of stone, clay etc. We can thus conclude that the migration model of the European Neolithization involved the movement of the carriers of the agrarian economy from Anatolia, who were variably influenced by either the Mesolithic or Neolithic populations from earlier migration events already living in the Balkans, which then established the LBK culture once they arrived at Brunn 2 and other sites of the formative LBK phase. The finding of remains of a possible first generation ANF/WHG admixed individual interred at Brunn 2 points to the economic, cultural and biological integration of HGs into the early LBK farming community. The full extent of contribution of European HGs to incoming Anatolian farmers remains an important subject for future work.