Ancient genome-wide DNA from France highlights the complexity of interactions between Mesolithic hunter-gatherers and Neolithic farmers
Abstract
Starting from 12,000 years ago in the Middle East, the Neolithic lifestyle spread across Europe via separate continental and Mediterranean routes. Genomes from early European farmers have shown a clear Near Eastern/Anatolian genetic affinity with limited contribution from hunter-gatherers. However, no genomic data are available from modern-day France, where both routes converged, as evidenced by a mosaic cultural pattern. Here, we present genome-wide data from 101 individuals from 12 sites covering today’s France and Germany from the Mesolithic (N = 3) to the Neolithic (N = 98) (7000–3000 BCE). Using the genetic substructure observed in European hunter-gatherers, we characterize diverse patterns of admixture in different regions, consistent with both routes of expansion. Early western European farmers show a higher proportion of distinctly western hunter-gatherer ancestry compared to central/southeastern farmers. Our data highlight the complexity of the biological interactions during the Neolithic expansion by revealing major regional variations.
INTRODUCTION
The Neolithic transition, which broadly describes the shift from foraging to farming, is one of the most important events in human history. In western Eurasia, the Neolithic way of life has been shown to spread westward from the northern Levant and potentially bypassing inner Anatolia from the seventh millennium BCE onward (1). Farming expanded along two main, archaeologically well-defined streams: along the Danube River to central Europe and along the Mediterranean coastline to Iberia (2, 3).
Recent large-scale genomic studies have shown that the spread of farming was mediated by early farmers through demic diffusion (4–11), while studies with a regional focus hinted at more complex regional processes of admixture between incoming farmers and local indigenous hunter-gatherers (HG) (6, 7, 12). The continental route is relatively well documented from southeastern to central Europe and shows a rapid expansion with very limited initial biological interaction, especially during the Neolithic Linear Pottery culture [Linearbandkeramik (LBK)] (5–8), followed by continued coexistence and evidence for cultural exchange with local HG for more than one millennium (13, 14). The Neolithic expansion in southwestern Europe is linked to the Mediterranean route (3), where a similar pattern of increasing admixture at least one millennium after the initial settlement of first farmer groups is observed in Iberia (5, 6, 11, 12, 15). The bigger picture emerging from these studies suggests little admixture between first incoming farmers and local indigenous HG in all targeted regions, followed by an increase of HG ancestry during an advanced phase of the Neolithic (5, 6). The contributing HG component detected in later Neolithic phases has been shown to be of local origin for the Carpathian Basin and Iberia (6, 12, 15).
From an archaeological perspective, it is not trivial to define and demonstrate the various modes of interaction between late HG and early farmers, and it requires high-resolution data and precise analytical methods to generate a clear and accurate spatiotemporal framework. Along the continental route, archaeological evidence of contact between Mesolithic and Neolithic groups increases toward the West (2). Overall, the lithic assemblages of the last HG are technically and stylistically rather similar to those of the early farmers. Signals of contacts with HG are mostly microlith types with regional traditions in the late Mesolithic blade-and-trapeze horizon. Archaeological signals of contact between immigrant farmers and HG have been reported not only from the earliest LBK sites in Hesse but also from Vaihingen in southwestern Germany. Here, coexistence between farmers, HG, and groups with ceramics such as the debated La Hoguette and Limburg wares is documented within LBK settlements (2, 16). While evidence of coexistence within settlements is no longer observable, particular HG objects of adornment have been observed from burials in western central Europe, suggesting a continuous tradition of coexistence of these different populations throughout the Middle Neolithic (17).
Looking at the Western Mediterranean route, the main uncertainties revolve around the origin of the dispersal of the Impresso-Cardial complex (ICC) in Italy. The late HG settlements are concentrated in the northeast of the peninsula, while the earliest farmers appear in the south, with little geographic overlap (18). In contrast, the discontinuity between the latest HG and earliest farmers is more pronounced in southern France, due to the marked stylistic differences of the latest HG blade and trapezes industries and the earliest Neolithic toolkits. These discrepancies reasonably argue for a cultural and population shift. However, only a few stratigraphic sequences provide a Mesolithic-Neolithic succession in southern France or Iberia and, in most cases, a clear stratigraphic gap is recorded, which is hard to reconcile with putative local interactions (19–21). Successive mixed assemblages, or rare Mesolithic reoccurrences within Neolithic sequences (21), appear only in the southern Alps, at least three centuries after the pioneer colonization of the Mediterranean coasts.
When reaching the northwestern parts of Europe, the dispersal of the Neolithic way of life became more complex, as shown by archaeological research (2, 18), resulting in a highly differentiated picture of interaction and exchange from the Mesolithic to the Neolithic (including regional variations). There is an even more pronounced chronological gap between the colonization of southeastern France by the ICC pioneer groups (terminus post quem 5850 BCE) (22) and of northeastern France by the LBK (5300 BCE) (23). It was suggested that both the concurrence of the two main streams (24, 25) of the west European Neolithic expansion up to the Atlantic coasts and the degree of interactions with autochthonous Mesolithic societies (3, 18) may have created the mosaic pattern of diversity seen in the material culture across the region in the subsequent centuries (text S1). This pattern has been well described at the cultural level, but no genomic data have been available to date, in particular from today’s France. A mitochondrial DNA study conducted on farmers from the Paris basin has highlighted a higher amount of characteristic HG mitochondrial haplogroups (notably U5) than described for regions in central or southern Europe (26), suggesting different processes in action.
Our study aims to cover the key geographic region of modern-day France and neighboring regions in Germany to unravel the complexity and variability in cultural and biological interactions between human groups during the earliest stages of the Neolithic period. This study region is uniquely suited as it encompasses both (i) the convergence of early farming communities from central Europe (Danubian route) with those from southern France (Mediterranean route) and (ii) variable forms of interaction with indigenous late HG.
RESULTS
Here, we present newly typed genome-wide data from 101 individuals from 12 sites from modern-day France and Germany (3 Late Mesolithic and 98 Neolithic, 7000–3000 cal BCE; table S1). We enriched next-generation ancient DNA (aDNA) libraries for ~1.2 million single-nucleotide polymorphisms (SNPs) using targeted in-solution capture (8), as well as an independent in-house capture array for the complete mitogenome [(27) modified following (6)], and sequenced these on Illumina platforms to an average depth per site of 0.6× for 1240K and 249× for mitochondrial capture. Libraries were built with a partial uracil–DNA–glycosylase (UDG) treatment, allowing the assessment of postmortem deamination patterns (2 to 29%) expected for aDNA data. We estimated potential contamination in mitochondrial capture DNA using (28) (table S3) and the nuclear genome by testing for heterozygosity of polymorphic sites on the X chromosome in males (table S2) (29). We estimated kinship to exclude first-degree relatives from downstream analyses (text S6) (30). We coanalyzed our new data with published ancient (n = 629) and modern (n = 2583) individuals from a panel of modern-day worldwide populations genotyped on the Affymetrix Human Origins (HO) panel (31, 32) as well as with ancient individuals genotyped on the 1240K panel. We excluded two individuals with low SNP coverage on the HO panel (<10,000) from downstream genome-wide analyses (table S1). To support the interpretation of our genetic results with a robust chronological context, we report new direct radiocarbon dates for 30 individuals (text S2 and table S1). A small set of 13 SNPs associated with phenotypes of interest were investigated. Results are presented in text S7 and show expected genotypes for Mesolithic and Neolithic European individuals.
We explored our new data qualitatively using principal components analysis (PCA) by projecting the ancient samples onto the genetic variation of a HO set of west Eurasians (Fig. 1C). Newly sequenced Late Mesolithic individuals from Bad Dürrenberg (BDB001) and Bottendorf (BOT004 and BOT005; see text S6 and fig. S17) fall within the variability of western HG (WHG). Individual TGM009, attributed to the Late Neolithic site of Tangermünde from an Elb-Havel context [Trichterbecherkultur (TRB); texts 1 and 2], occupies an intermediate position between WHG and Neolithic farmers. Newly sequenced Neolithic individuals cluster with previously published data, forming two regional subgroups: the first including individuals from Neolithic central and southeastern Europe and the second individuals from Neolithic western Europe (Iberian Peninsula, France, and British Isles), which is shifted slightly toward WHG on PC1. Formal f-statistics of the form f4(Mbuti, European_HG; test, Anatolia_Neolithic), where test represents Neolithic groups, confirm these observations via different degrees of affinity with European HG (table S9).
Fig. 1 Spatial, temporal, and genetic structure of individuals in this study.
(A) Geographic map showing sample locations. (B) Chronological timeline: Maximal chronological range according to calibrated radiocarbon date ranges (2-sigma) for each site/individual. (C) Principal components analysis (zoom in). Published ancient (no outlined symbols) and newly reported (black outlined symbols) individuals projected onto 777 present-day west Eurasians (gray circles).
Early Neolithic individuals from Stuttgart-Mühlhausen (SMH), Schwetzingen (SCH), and Halberstadt (HBS) in Germany [44 published and 42 newly reported, of which 5 low-coverage individuals (<18,000 SNPs mapping on the HO panel) are shifted toward the centroid] representing the LBK horizon (17 sites from Hungary, Austria, and Germany) form a homogeneous genetic group of early farmer populations in central Europe. Neolithic groups from France cluster with western Neolithic individuals, while southern French sites Pendimoun (PEN) and Les Bréguières (LBR), attributed to the ICC, are shifted even further toward WHG. This suggests a higher HG component than in any other early Neolithic individuals, and thus a different history of admixture compared to contemporaneous western early farmer groups from Iberia. For both of the sites PEN and LBR, we split the individuals with significantly different HG ancestry contributions into subgroups A (less) and B (more HG ancestry; see table S7). Individuals from southern France do not group with ICC individuals from the Adriatic region, who fall within the subgroup of southeastern and central European farmers (Fig. 1C).
Both central and western groups are still present during the fifth and early fourth millennium BCE, suggesting a geographic border along the Rhine River for regions that are situated north of the Alps. Individuals from the Middle Neolithic French sites of Gurgy (GRG), Prissé-la-Charrière (PRI), and Fleury-sur-Orne (FLR) in the northern half of France appear homogeneous (one PRI individual outside of the cloud has low SNP coverage), while the individuals from the site Obernai (OBN) form three groups: one with other contemporary western individuals (OBN A), one with a stronger HG component (OBN B), and a third one with central European farmers (OBN C), despite sharing similar cultural and chronological backgrounds (Fig. 1C and text S2).
Quantifying universal HG ancestry in European farmer groups
We used qpAdm (ADMIXTOOLS) (5) to estimate the proportion of European_HG ancestry in all new and published Neolithic individuals with directly associated radiocarbon dates to track changes through time (MODEL A; Fig. 2, text S9, and table S10). We observe a notably similar pattern with negligible amounts of HG ancestry during the first periods of attested farming in each region of central and southwestern Europe, followed by a steady increase after centuries of established farming. The last individuals without any traceable HG ancestry disappear around 3800–3700 cal BCE. Southeastern Europe shows a specific pattern in individuals linked to the Iron Gates region, and the British Isles show a sudden and constant higher HG component at the time of the arrival of the Neolithic in the region. Our new western LBK groups (east of the Rhine) confirm the previous estimates from the Middle-Elbe-Saale region and Transdanubia, which is in contrast to southeastern European sites such as Malak Preslavets in Bulgaria, for which higher and more variable estimates were reported (6). However, the situation is different in today’s France, where we observe not only the highest HG ancestry proportion overall when compared to other regions in Europe but it is also found in the oldest individuals from the southern sites of PEN and LBR. This observation is also supported by uniparental markers. Y chromosome lineages in western early farmers in the southern region are exclusively derived from HG (I2a; table S5 and text S5). In contrast, mitochondrial DNA results show a more universal Neolithic diversity profile, as previously reported [i.e., (4, 5)], with only two haplotypes (U5 and U8) that are potentially of HG origin (table S4 and text S4).
Fig. 2 HG ancestry proportions over time.
(A) Overall timeline: Results of qpAdm (MODEL A) modeling of European_HG (represented by Loschbour, La Braña, and KO1; y axis) and Anatolia_Neolithic ancestry for each individual with a direct radiocarbon date ranging between 6200 and 2800 BCE (x axis). See table S10 for further details. (B) Regional timelines: Plots of the data shown in (A), separated by geographic location, to show regional signals.
Two individuals from the Middle Neolithic OBN site located immediately west of the Rhine (Fig. 1A) also show a high proportion of the HG ancestry component (OBN B), in contrast to LBK sites east of the Rhine. To quantify the extra-HG proportion, we modeled LBK from Germany as a mixture of Anatolia_Neolithic and European_HG, and subsequently OBN subgroups as a more proximal mixture model of LBK and European_HG. We find this model well supported (table S11), and obtained up to 31.8% of excess HG ancestry for OBN B. Given the absence of a strong HG contribution in LBK groups east of the Rhine, we assume this to be a local contribution during the centuries following the arrival of first farmers. Moreover, male individuals from OBN carry exclusively the Y chromosome haplogroups I2a1a2 and C1a2b, attributed to HG groups (text S5), providing further evidence for a greater amount of the HG contribution in this region.
Genetic structure in European HG
The shared recent ancestry of early European farmers with Anatolian Neolithic farmers (32) and the rapid expansion across mainland Europe render it challenging to distinguish between different Neolithic routes at the genomic level as proposed by archaeology. Although we observe an increase in HG ancestry from east to west in published regions, this might simply be the result of geographic distance from the farming origin in line with Lipson et al. (6). To test for signals of different routes of expansion, we therefore made use of the increasingly emerging geographic structure in European HG ancestry (6, 7, 15, 33). Here, post-glacial European HG ancestry can be described by three main clusters (Fig. 3 and fig. S1) with clines of admixture between each: (i) GoyetQ2-like ancestry, (ii) Villabruna-related WHG ancestry, and (iii) eastern HG (EHG) ancestry. In addition, two clines between the clusters can be observed: one between GoyetQ2 and WHG formed by Iberian HG, and the other between EHG and WHG formed by southeastern European, Scandinavian, and Baltic HG (7, 12, 15).
Fig. 3 Multidimensional scaling plot HG individuals.
Multidimensional scaling plot of genetic distances based on f-statistics of the form f3(Mbuti; test, test) between Eurasian HG individuals. Newly reported individuals have a black outline.
Following this observation, we also performed f4-statistics of the form f4(Mbuti, test; Villabruna, EHG) testing all HG individuals (fig. S2A). Here, negative f4-values for all European groups indicate shared ancestry with the Villabruna individual and thus the WHG cluster, whereas positive values show more attraction to EHG. Using qpAdm, we modeled HG individuals as a mixture of the distal sources Villabruna, EHG, and GoyetQ2 to establish the best-fitting subclades of European HG individuals as sources of ancestry (MODEL B; text S9, fig. S4, and table S12). We also added Anatolia_Neolithic to fit some of the HG groups who show signs of admixture with early farmers. We chose the most parsimonious model, i.e., fitting the data with the minimum number of sources, for each set of results. In addition to the GoyetQ2 component present in Iberian individuals, as shown recently (15), we confirm that individuals from the EHG-WHG cline, including southeastern European and Scandinavian HG, can be distinguished from WHG individuals, as they carry a substantial proportion of EHG ancestry (7). Figure 3 shows that the newly reported Mesolithic individuals from Germany (BDB and BOT) also form part of the WHG cluster.
In the following, we leverage the observation that, at the time of the arrival of farming, this broad ancestry cline between EHG and WHG covers key geographical regions in central/southeastern Europe (spanning today’s Germany, Hungary, Serbia, Romania, and Ukraine), which overlap with the proposed Danubian route of the Neolithic expansion. Consequently, we hypothesize that Early Neolithic farmers of the Mediterranean route associated with ICC Ware and especially those arriving in France and Spain show less ancestry from the EHG side of the EHG-WHG cline and instead carry an admixture signal dominated by Villabruna-related WHG and GoyetQ2-related ancestries.
Tracing different HG ancestries in farmer groups
Analogous to the test with HG individuals, we calculated f-statistics of the form f4(Mbuti, test; Villabruna, EHG) for all Neolithic groups (fig. S5A) to estimate the affinity of Neolithic groups (test) to the distal HG representatives of the east-west cline. Like the HG groups (fig. S2A), the f4-values from the Neolithic groups are mostly negative, indicating shared ancestry with Villabruna, but form a wide gradient reaching 0, which implies equal amounts of shared ancestry with EHG for some groups located in southeastern Europe. We used qpAdm to quantify these sources in MODEL B (text S9) for all Neolithic European groups. However, using these distal sources due to the very small amount of HG component in many early Neolithic groups, the model fails to reliably detect an EHG component (table S12).
We therefore restricted the analysis to more proximal sources. We chose the deeply covered Mesolithic individual Loschbour from Luxembourg as representative of WHG ancestry (west of the Rhine) and therefore as a proxy for the Mediterranean route (34). We selected the Hungarian individual KO1, which looks genetically like a HG, though found in a farming context, for the HG ancestry of the continental route (10). We first replicated the f-statistics of the form f4(Mbuti, test; Loschbour, KO1) (fig. S5B). Here, positive f4-values indicate shared ancestry with KO1 as a proxy for the EHG-WHG cline. In turn, negative f4-values indicate excess shared ancestry with Loschbour and thus the WHG cluster. The results show a tendency for farmer groups east of the Rhine to share more ancestry with KO1 and farmer groups west of the Rhine with Loschbour. Notable exceptions are individuals from the Middle Neolithic Blätterhöhle group (6, 13), individual N22 from Poland (35), and our individual TGM009, all of which show a strong affinity to WHG-related individuals, although they are located east of the Rhine. On the basis of our geographic rationale, we then used qpAdm to model HG components in the Neolithic groups, including both HG sources that could be encountered en route for the proposed expansions (MODEL C; text S9, fig. S6, and table S13). We chose the best-fitting model given by qpAdm in cases where both HG sources were supported. Figure 4 summarizes these results, illustrating the diverse HG ancestries in time windows across major regions in Europe, adapting the models used to the targeted groups (table S14).
Fig. 4 Maps of variable sources of HG and Anatolian Neolithic ancestries through time.
(A) Proportion of distal sources of Villabruna, EHG, Goyet_Q2, and Anatolia_Neolithic of post-LGM HG individuals (14000–4000 cal BCE) estimated according to MODEL B (table S12). The following five panels (B to F) show Neolithic farmer groups (6000–3500 cal BCE) modeled with the proximal sources Anatolia_Neolithic, KO1, and Loschbour according to MODEL C, in time increments of 500 years each. Transparent colors indicate individuals or groups not sufficiently supported by the models (P < 0.05). Note that not only N22 from the Polish BKG, TGM009, and Blätterhöhle, which carry a substantial proportion of HG ancestry, but also Anatolian_Neolithic ancestry have been modeled with MODEL B and were added to (E) and (F), respectively. See table S14 for model details.
We caution that because the amount of HG ancestry is very low in many Neolithic groups (<10%), it remains difficult to characterize the ultimate source reliably. Nevertheless, our admixture patterns from supported models show clear geographic signals. Neolithic groups associated with the LBK in central Europe (Hungary, Austria, and Germany) carry a small HG proportion, which was likely derived from admixture with HG individuals of the EHG-WHG cline and could have occurred in southeastern Europe during a preceding phase of the Neolithic expansion around 6000–5400 BCE. When using f-statistics of the form f4(Mbuti, test; BDB001, KO1) with our new Mesolithic genome (BDB001) from the Middle-Elbe-Saale region in central Germany as a geographically local HG proxy (instead of using Loschbour, which is located west of the Rhine), we do not find support for a local attraction for LBK groups, but the same pattern as for Loschbour (fig. S5C). This suggests that additional gene flow from neighboring Loschbour-like HG such as BDB001 in central Europe was negligible in the first Neolithic groups. However, the German Baalberge group (4000–3500 BCE) shows a marked increase of such HG ancestry, as well as individuals from the Blätterhöhle group, as has been suggested (5, 6), compared to a combination of both KO1-like and Loschbour-like ancestries for LBK groups (6). We can now show that this increase in WHG ancestry (up to 21.3 ± 1.5% in Baalberge; table S13) is driven by either local Loschbour-like ancestry or an expansion of farming groups from the west carrying this signal during the fifth millennium BCE, as suggested by archaeological data (36). For all studied Neolithic groups west of the Rhine, we observe a different pattern with a higher HG ancestry proportion, even for earliest groups that appears to be of a local (Loschbour-like) HG origin, consistent with archaeological data (2).
Late survival of HG ancestry in central Europe
In contrast to Middle-German Neolithic groups, the Late Neolithic individual TGM009 (~3300 BCE) found in an TRB/Elb-Havel group context shows a different pattern of HG ancestry. When modeled with distal sources, this individual carries 63.6 ± 5.2% Villabruna-related ancestry (see MODEL B; table S12 and Fig. 4). When we use proximal sources (MODEL C; table S13), we estimate the HG component to be split into 48.1 ± 6.4% KO1-related ancestry and 25.8 ± 6.1% Loschbour-like ancestry (P = 0.12). To test whether the subtle EHG-related signal observed in TGM009 (through KO1) could also come from Scandinavian HG, as suggested by the archaeological records of regional late contacts with the Scandinavian Mesolithic (37), we applied f-statistics of the form f4(Mbuti, test; Hungary_KO1, Sweden_Motala_HG). The f4 value is significantly negative (f4 = −0.0035, Z score = −6.383), indicating that TGM009 shares more ancestry with southeastern Europe HG than with Scandinavian HG (table S16). Consequently, qpAdm models with Sweden_Motala_HG in place of Hungary_KO1 received a poor fit (P = 1.32−04; table S16).
In the specific case of TGM009, we also tested whether individuals of the Pitted Ware Culture (PWC) group (~3200–2300 BCE) (38), considered to be “Neolithic HG,” would be a suitable contemporaneous proximal source with MODEL C set of outgroups (table S16). We find that the tests still support the three-way model with Anatolia_Neolithic, Luxembourg_Loschbour, and Hungary_KO1 as a third source (P = 0.12). However, when we add Sweden_Motala_HG to the outgroups, we find the best support for a four-way mixture model (P = 0.087) of 21.7 ± 2.4% Anatolia_Neolithic, 24.4 ± 6.2% Luxembourg_Loschbour, 12.6 ± 4.4% Sweden_PWC, and 41.3 ± 7.3% Hungary_KO1. The small but stable contribution of PWC groups at the northern fringes of the loess belt adds additional complexity to the modes of interaction between Late Neolithic and the last HG groups in central Europe.
Estimating admixture dates between HG and early farmers
To gain further insights into the timing of the admixture of the HG ancestry proportion, we used the software package DATES to estimate the timing of admixture events between early farmers and HG, given in generations in the past (Fig. 5, table S17, and text S10). The resulting date estimates, when ordered temporally and geographically, echo patterns of both groups that are also visible in the PCA (east versus west of the Rhine). It confirms that this pattern is dependent not only on the amount of HG ancestry but also on the qualitative signature of the HG ancestry, as supported by fig. S7. Together, the date estimates for the four southern French groups from sites PEN A and B and LBR B suggest admixture with local HG relatively soon (100 to 300 years) after the arrival of Neolithic farmers about 5850 cal BCE, with an admixture date about 5740–5450 cal BCE for all groups (table S17). These date estimates agree with local archaeological data for the establishment of early farming in this region (18), although we cannot exclude a scenario in which admixture occurred on the Italian Peninsula shortly before.
Fig. 5 Admixture dating per group.
Admixture date estimation according to DATES software, obtained in generations. Radiocarbon date intervals (given with 2-sigma) for each site are black lines; blue diamonds are the estimated admixture date (SE = 1). Admixture is estimated with two sources: European_HG and Anatolia_Neolithic. Number of years calculated on the basis of 28 years for one generation. Admixture time calculated according to the oldest date of the radiocarbon interval for each group (table S17 and text S10).
In the rest of France, post-LBK groups show an older admixture date, which, together with a strong local HG component, place the admixture event with WHG-related HG during the first phase of the local Neolithic (Fig. 4). In contrast, the site OBN reveals a heterogeneous genetic signal forming three distinguishable groups with different HG proportions and corresponding admixture dates (tables S7, S13, and S17). The differences between the date estimates for the three groups suggest that the OBN group shows a genetic substructure rather than a recent ongoing admixture event.
Tracing GoyetQ2 ancestry in France
We also used qpAdm with GoyetQ2, Villabruna, and Anatolia_Neolithic as sources to track Magdalenian-associated GoyetQ2-like ancestry in the new western Neolithic groups (MODEL D; text S9 and table S15) (15). The GoyetQ2-related component is observed in the sites of PRI (6 ± 3%) in western France and GRG (3.7 ± 1.3%) from the Paris basin. Similar to Iberian Neolithic groups, the PRI group (4300–4200 BCE) can be modeled successfully with GoyetQ2 as one of the HG components (fig. S8 and table S15), suggesting either admixture with local HG in western France who retained this post-Last Glacial Maximum component or genetic contacts with Neolithic Iberia in later phases. Of note, individuals from the site PRI on the Atlantic coast indicate a slightly more recent admixture date with European_HG (about 5200 BCE), which is consistent with a later arrival of Neolithic groups in westernmost France (table S17). We plotted these results in fig. S9 to complement the summary map with the GoyetQ2-related component in Neolithic groups, which adds another layer of complexity to the overall genetic picture of Neolithic Europe.
Connections with Britain and Ireland
To address questions at broader scale on the western fringe of Europe, we also investigated the relationship between Britain and Ireland and the European mainland during the Neolithic. Here, we performed f4-statistics of the form f4(Mbuti, Britain/Ireland groups; PRI, test), which measures whether Neolithic groups from Britain and Ireland share more genetic drift with the French Atlantic group PRI than with a test population (see fig. S10). As previously described (39, 40), we confirm that British Neolithic groups share affinities with the Mediterranean Neolithic (LBR_A, France_MN, and Iberia_MN), which is also visible in the WHG-rich HG proportion (Fig. 4). However, on the basis of the results from our French Neolithic sites, we suggest that English, Welsh, and Scottish groups are connected to the Mediterranean Neolithic sphere not only via the Atlantic coast but more plausibly also via Normandy (FLR), the Paris Basin (GRG), and southern France (LBR A, post-ICC groups).
DISCUSSION
The expansion of Anatolian Neolithic farmer–related ancestry across Europe has been described for many geographic regions (5–7, 9–11, 39, 40). A recurrent pattern of increasing HG gene flow in individuals associated with a farming lifestyle many centuries after initial contact and tentatively termed “HG resurgence” (5) has hitherto been observed in the Iberian Peninsula, northern/central Europe, the Carpathian Basin, and the Balkans (6). However, a comparison of genomic data in the light of the proposed main routes of Neolithic expansion has not yet been attempted. Our new results provide important insights for the region of modern-day France where both routes had intersected. This region has not been documented so far, but is ideally positioned to address these questions. The different Mesolithic genetic substratum in Europe (6, 7, 15, 33) enables us to track both Neolithic expansion routes in the form of the quantity and quality of the admixed HG component observed in Neolithic groups.
Neolithic groups from southern France, which are a part of the Mediterranean/ICC route of expansion, show a different genetic profile when compared to early periods of expansion in other regions, with a substantially higher HG component than groups associated with the continental route (up to 56 ± 2.9% of HG component for PEN B; table S10). Although the sites PEN and LBR date ~400 years (i.e., ~16 generations) later than the first local coastal settlements in Liguria, Provence, and Mediterranean Languedoc, they suggest a recent local admixture event (between three and six generations earlier; Fig. 5 and table S17). Archaeological research has argued for increased interaction between incoming farmers and indigenous HG in the western Mediterranean during a second stage of the Neolithization process and especially in areas with higher HG population densities, e.g., the Tosco-Emilian Apennine and Po plain (18). We are now able to confirm that these contacts left a traceable biological signal during the Neolithic expansion in southern France. From an archaeological perspective, this suggests that HG have contributed to the clear changes observed within the material culture postdating the pioneer phase. Note that ICC individuals from the eastern Adriatic coast have only a very small amount of HG ancestry with a greater affinity to central European groups (see table S8). This fits with the hypothesis of a differentiation of technical traditions within material cultures observed from both sides of Apennine Mountains in Italy: an Adriatic tradition connected to the Balkans and a Tyrrhenian one whose origin is still unknown (41). It is tempting to associate such a strong HG component on the Tyrrhenian side with the characteristic/specific pottery traditions observed in this same region and to consider these original traditions the result of a HG reinterpretation (41). However, the scarcity of genomic data available from central and southern Italy currently does not allow this hypothesis to be tested directly. Moreover, ICC individuals from the Iberian Peninsula also carry less HG ancestry. Together, this rejects the hypothesis that ICC-associated individuals represent a uniform genetic horizon per se and argues for more regionally nuanced scenarios of interaction.
Early central European farmers carried a very small percentage of HG component (about 5% on average), which is most likely nonlocal and instead derives from admixture in Transdanubia during an early stage of expansion. This not only matches general observations in the archaeological record (2) but could also explain why the earliest LBK lithic assemblages resemble those of the Late Mesolithic Blade and Trapeze Complex. It also confirms previous aDNA studies (5, 6), which argued for a fast spread of the first farmers across the German loess regions. LBK groups from southwestern and eastern Germany share more affinities with KO1 than with Loschbour (fig. S5). The admixture date estimates for southern German sites, SMH and SCH, 19.2 ± 3.8 and 12.3 ± 8.2 generations, respectively, are younger or contemporaneous to admixture date estimates from the Carpathian Basin and Austria (Fig. 5 and table S17). The temporal lag, and the subtle increase of shared KO1-like HG ancestry, allows us to trace LBK groups chronologically and matches well with the proposed model of LBK expansion from the core region in Transdanubia as suggested by archaeological research (14). However, the present picture does not explain the increasing amount of archaeological evidence of contact with HG and groups with southern influences in more western sites, particularly during the early LBK, such as the debated La Hoguette and Limburg phenomena (16).
In contrast to the situation in central Europe, regions west of the Rhine show a different profile during the following fifth millennium BCE. Here, first farmers carry a higher local Loschbour-related HG component, later increasing to up to 33.3 ± 3% in several individuals in Alsace (OBN B). Mitochondrial data support this finding with a higher average proportion of HG-affiliated haplotypes (U5 and U8) in all French groups from the fifth millennium BCE (15.5%; table S4 and text S4) compared to LBK groups east of the Rhine (1.4%). Although we do not have genomes from LBK individuals from northern France, we can infer from this HG component that admixture processes happened locally after the first farmers arrived. When modeling the European_HG component in the three groups from OBN with qpAdm, we observe an increase in the HG component of between 4.3 and 31.8%, which we attribute to local Loschbour-related sources (table S11). However, this approach cannot be directly applied to other contemporaneous French groups as potential additional movements linked to the Mediterranean route of Neolithic expansion that are supposed to have followed the first LBK farmers’ arrival in northern France could complicate the simple assumption of a two-way HG/farmer interaction (26, 42). The current genomic data do not allow us to distinguish whether the detected signal is confounded by a potential southern contribution. However, the estimated admixture dates for GRG and FLR in northern France indicate older admixture events that occurred more than 30 generations (840 to 930 years) earlier (Fig. 5 and table S17). In accordance with the established chronology of first Neolithic settlements in the French territory, the overlapping/synchronous date estimates obtained for southern ICC sites are consistent with the signal of a first HG contribution in the south of France, followed by a subsequent northward expansion of groups carrying this HG legacy (42).
The GoyetQ2-like HG component detected in individuals from the sites PRI and GRG suggests connections with the Iberian Peninsula where this post-LGM residual component is found, either by an admixture with HG (which could also be local) carrying this component or by contact with first farming communities from Iberia or because of exchanges with Neolithic groups from Iberia during the subsequent fifth millennium. To date, archaeological data are compatible with all three hypotheses (25, 43). However, while we can show affinities between PRI on the Atlantic coast and the Iberian Peninsula, the previously described affinities between the British Isles and Mediterranean Neolithic (39, 40) are currently best explained by sites mainly in Normandy and the Mediterranean area via the Paris Basin. English, Scottish, and Welsh groups show more genetic affinities to northern and southern France and Iberia than to western France (fig. S10). In contrast, Neolithic Ireland shows less affinities with the northern French coastline and the Mediterranean area than other British groups, which could be explained within an Atlantic framework. This whole pattern is in line with the hypothesis of two different phenomena and speeds of Neolithic expansion to the western and eastern British Isles, as proposed by archaeological data (44).
In summary, our study highlights a diverse pattern of cultural and biological interactions between first farmers and indigenous HG along the western Mediterranean coastline and west of the Rhine, which confirms a high variability of processes during the Neolithic expansion, distinct in the proportion of the HG component as well as processes and duration. The genetic structure among HG groups allowed us to track local admixture in early farmers, which not only is higher west of the Rhine compared to central and southeastern Europe but also is largely attributed to local and distinctly WHG-related sources. We show that nuanced sampling as well as increasing cohort numbers of both HG and early farmers can help to unravel the regional dynamics of the Neolithic transition and aid in refining our understanding of the underlying processes and developments over time, at both micro- and macroregional scales. On the basis of our observations, we find that broad-brush models are increasingly less likely to reconcile the full spectrum and details of the farmer-forager interactions and thus advocate the use of models with more regional focus in future studies.
MATERIALS AND METHODS
Study design
Archaeological samples
We processed and analyzed samples from 101 individuals, 3 Mesolithic and 98 Neolithic (7000–3300 cal BCE), from 12 sites in modern-day France and Germany. A detailed description of each site is provided in the Supplementary Materials, including contextual data and available radiocarbon dates (text S2 and table S1).
We sampled, extracted, and prepared DNA for next-generation sequencing in three different dedicated aDNA facilities. Individuals from PRI and PEN were sampled in the clean room of the Laboratoire PACEA, Bordeaux University, France. Samples coded “XN” from the site SMH were processed at the Interfaculty Centre for Archaeology at the University of Tübingen, Germany, and all other samples were processed at the Max Planck Institute for the Science of Human History (MPI-SHH) in Jena, Germany.
Sampling
We sampled mainly petrous bones and tooth and, for a few individuals, also other type of skeletal remains (femur and phalanges; table S1). Samples were first irradiated with ultraviolet light for 30 min on all sides. All teeth were cleaned with a low-concentration bleach solution (3%). Teeth prepared in Tübingen and Jena were cut along the cementum/enamel junction, and powder was collected by drilling into the pulp chamber. Teeth prepared in Bordeaux were ground to fine powder completely. Petrous bones were either cut in half and powder drilled from the denser regions around the cochlea (GRG) or drilled from the outside (FLR, LBR, OBN, and PEN). Before that, a layer of the bone surface was mechanically removed.
DNA extraction
From the extraction step onwards, all samples were processed at MPI-SHH in Jena, except the samples labeled XN, which were processed at the University of Tübingen. At MPI-SHH, DNA was extracted following the protocol described by Dabney et al. (45), with the High Pure Viral Nucleic Acid Kit (Roche). The same protocol was used in Tübingen with MinElute columns (Qiagen).
Library construction
One hundred and four double-stranded libraries were built with unique index pairs (46, 47). We applied the partial UDG (UDG half) protocol to remove most of the aDNA damage while preserving the characteristic damage pattern in the terminal nucleotides (48). A separate set of non-UDG libraries prepared in Tübingen was used for the mito-capture of SMH samples (table S3).
Shotgun screening and capture (1240K and mitochondrial)
We first screened all indexed libraries via shotgun sequencing of 5 million reads on an Illumina HiSeq 4000 sequencer using either a single [1 × 75–base pair (bp) reads] or double end (2 × 50–bp reads) kit. We used EAGER (49) to process the raw data and to select libraries with >0.1% endogenous human DNA and those showing characteristic damage aDNA patterns for downstream SNP capture. Selected libraries were hybridized in-solution to different oligonucleotide probe sets synthesized by Agilent Technologies to enrich for 1,196,358 informative nuclear SNP markers (8) and an in-house capture for the complete mitogenome following Maricic et al. (27) and modified after Haak et al. (5).
aDNA data processing
Read processing and aDNA damage
After demultiplexing, raw sequence data were processed using EAGER. This included clipping adaptors with AdapterRemoval (50), mapping with BWA (Burrows-Wheeler Aligner) v0.7.12 (51) against the Human Reference Genome hs37d5, and removing duplicate reads with the same orientation and start and end positions. After using mapDamage v.2.0.6 to observe characteristic aDNA damage patterns, we used BamUtil (https://genome.sph.umich.edu/wiki/BamUtil:_trimBam) to clip two bases at the ends of each read for each sample to remove residual deaminations and 10 bases for reads from the non-UDG libraries performed for mito-capture on SMH samples (tables S2 and S3).
Sex determination
Following Mittnik et al. (52), we determined the genetic sex by calculating the number of reads mapping to each of the sex chromosomes with respect to the autosomes. We set a threshold of Y ratio < 0.05 for a female and Y ratio > 0.4 for a male (text S3 and table S2).
Contamination estimation
We used the ANGSD (Analysis of Next Generation Sequencing Data) package to test for heterozygosity of polymorphic sites on the X chromosome in male individuals, applying a contamination threshold of 5% (table S2) (29). For mito-captured samples, we estimated contamination levels using contamMix 1.0.10 (53) by comparing the consensus mitogenome of the ancient sample to a panel of 311 worldwide mitogenomes as a potential contamination source (table S3). Fourteen individuals from FLR and SMH gave a rate of contamination over 5%, but while checking the read assembly visually, we found no consistent pattern of contaminating mitochondrial DNA lineages. On the basis of low contamination estimates for males on the X chromosome—when available—we opted to include these individuals in the downstream analyses. Libraries used for mito-capture for SMH samples are non-UDG–treated, and hence, damage is also biasing our contamination estimates (53).
Genotyping
We genotyped our bam files with pileupCaller (https://github.com/stschiff/sequenceTools/blob/master/src/SequenceTools) by randomly calling one allele per position considering the human genome as pseudo-haploid genome. We called the SNPs according to the Affymetrix HO panel (~600K SNPs) (31, 32) and the 1240K panel (8). Numbers of SNPs covered at least once are given in table S1.
Mitochondrial and Y chromosome haplogroup assignment
To process mitochondrial DNA data, we extracted reads from mito-capture data when available or from 1240K data using samtools v1.3.1 (54) and mapped these to the revised Cambridge reference sequence. We called consensus sequences using Geneious R8.1.974 (55) and used HaploGrep 2 to determine mitochondrial haplotypes (text S4 and table S4) (56).
We used pileup from the Rsamtools package (57) to call the Y chromosome SNPs of the 1240K SNP panel from all male individuals (mapping quality ≥30 and base quality ≥30). We manually assigned Y chromosome haplogroups using pileups of Y-SNPs included in the 1240K panel that overlap with SNPs included on the ISOGG SNP index v.14.07 (text S5 and table S5).
Kinship analysis
We estimated the degree of genetic relatedness between our individuals by applying Relationship Estimation from Ancient DNA (READ) (text S6) (30).
Phenotypic and functional analysis
We investigated individual genotypes for 13 SNPs associated with phenotypes of interest (8) by calculating genotype likelihoods based on the number of reads from our bam files (quality <30) for each specific position to determine the presence of the ancestral (non-effect) or derived (effect) alleles. We restricted the analysis to individuals with >200,000 SNPs mapping to the HO panel (text S7).
Population genetic analysis
Merging new dataset with published data
Two datasets were created for genome-wide analysis. We first merged our new data with published ancient data to the HO panel (~600K SNPs) (31, 32). This dataset was only used for PCA (see below). We then merged our data with published ancient data to the 1240K SNP panel (8) including 300 present-day individuals from 142 populations sequenced to high coverage (58). This second dataset, restricted to the autosomes, was used for all other population genetic analyses.
Group labels
Table S6 describes analysis groups and labels, for new, as well as published ancient samples.
Principal components analysis
We performed PCA using the “smartpca” program v10210 (EIGENSOFT) on the HO dataset (Fig. 1C) (59). We computed principal components from 777 present-day west Eurasians on which ancient individuals were then projected using the options lsqproject: YES and shrinkmode: YES. We excluded individuals with less than 10,000 covered SNPs. A zoomed-out version of the PCA is presented in text S8.
f-Statistics
Outgroup f3-statistics were calculated using qp3Pop and f4-statistics using qpDstat with the f4 mode from ADMIXTOOLS (31) on the 1240K SNP panel. To investigate HG diversity, we performed outgroup f3-statistics of the form f3(test, test; outgroup) to create a similarity matrix, which was then used to generate the heatmap using the heatmap.2 function of the R-package gplots (fig. S1) (60). By calculating the values of 1-f3, we created a dissimilarity matrix from our f3-statistics. We performed multidimensional scaling on two dimensions for the dissimilarity matrix and then plotted the principal dimensions (Fig. 3). We computed SEs using the default block jackknife approach. We report and plot three standard errors for the f4-statistics test in supplementary figures (figs. S2, S3, S5, and S9).
qpAdm modeling
We used qpAdm and the 1240K SNP panel to estimate ancestry proportions (ADMIXTOOLS) (5) with respect to a basic set of 11 outgroups: Mbuti, Papuan, Onge, Han, Karitiana, Mota, Ust_Ishim, MA1, Czech_Vestonice, Caucasus_HG, and Israel_Natufian. We then added additional outgroups according to the tests performed (see text S9 and tables S10 to S13 and S15 for further details).
Admixture date estimation
We used the method DATES v.753 to leverage patterns of ancestry covariance to estimate the date of admixture between European_HG and Anatolia_Neolithic (see text S11 and table S17 for further details).
Acknowledgments
We thank F. Aron, M. Burri, A. Mittnik, K. Nägele, R. Radzeviciute, E. Späth, R. Stahl, A. Wissgott, and G. Brandt for technical support in the DNA analyses. We thank the teams at MPI-SHH-Archaeogenetics and PACEA, University of Bordeaux for continued support and discussion, especially S. Clayton, J. Y. Fellows, T. Lamnidis, L. Papac, and V. Villalba-Mouco. We thank S. Nakagome for a script for allele frequency estimation. We thank H. Duday, D. Hoffmann, B. Gehlen, and A. Powell for critical discussions. We thank H.-C. Strien for the discussion of SMH radiocarbon dates. We thank the project Counter Culture: Investigating Neolithic social diversity for analyses on isotopic data for OBN and analyses carried out at the BioArCh facility, Department of Archaeology, University of York. Funding: M.R. was supported by the Fyssen Foundation (post-doctoral stipend, 2017–2018). C.J. was supported by the New Faculty Startup Fund from Seoul National University. The study was funded by the Max Planck Society, and the French (ANR) and German (DFG) Research Foundations, INTERACT project, ANR-17-FRAL-0010, DFG-HA-5407/4-1, 2018-2021 (M.-F.D., W.H., and M.R.), and the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program under grant agreement no. 771234 – PALEoRIDER (W.H. and A.B.R.). We thank the DFG, grant no. RO 4148/1 (E.R.), for funding SMH radiocarbon dates. Author contributions: M.R., M.-F.D., and W.H. conceived the study. K.W.A., D.B., S.F., E.G., D.G., L.L., P.L., H.M., H.R., E.R., S.R., C.S., L.S., and J.W. provided archaeological material and advised on the archaeological background and interpretation. M.R., İ.K., and M.-H.P. performed the laboratory work. M.R. performed the data analyses, with C.J., S.S., A.B.R., J.K., and W.H. providing guidance and methodologies. M.R., D.G., D.B., M.-F.D., and W.H. wrote the manuscript with input from all coauthors. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Individual bam files (aligned human reads) are available through the European Nucleotide Archive under accession number PRJEB36208.
Supplementary Material
File (aaz5344_sm.pdf)
File (aaz5344_tables_s1_to_s5.xlsx)
File (aaz5344_tables_s6_to_s17.xlsx)
REFERENCES AND NOTES
1
B. Horejs, M. Milic, F. Ostmann, U. Thanheiser, B. Weininger, A. Galik, The Aegean in the Early 7th millennium BC: Maritime networks and colonization. J. World Prehist. 28, 289–330 (2015).
2
D. Gronenborn, The persistence of hunting and gathering: Neolithic western temperate and Central Europe, in The Oxford Handbook of the Archaeology and Anthropology of Hunter-Gatherers, V. Cummings, P. Jordan, M. Zvelebil, Eds. (Oxford University Press, 2014), pp. 787–804.
3
J. Guilaine, Aspects de la néolithisation en Méditerranée et en France, in The Widening Harvest. The Neolithic Transition in Europe: Looking Back, Looking Forward, A. Ammerman, P. Biagi, Eds. (Archaeological Institute of America, 2003), pp. 189–206.
4
B. Bramanti, M. G. Thomas, W. Haak, M. Unterlaender, P. Jores, K. Tambets, I. Antanaitis-Jacobs, M. N. Haidle, R. Jankauskas, C.-J. Kind, F. Lueth, T. Terberger, J. Hiller, S. Matsumura, P. Forster, J. Burger, Genetic discontinuity between local hunter-gatherers and central Europe’s first farmers. Science 326, 137–140 (2009).
5
W. Haak, I. Lazaridis, N. Patterson, N. Rohland, S. Mallick, B. Llamas, G. Brandt, S. Nordenfelt, E. Harney, K. Stewardson, Q. Fu, A. Mittnik, E. Bánffy, C. Economou, M. Francken, S. Friederich, R. Garrido Pena, F. Hallgren, V. Khartanovich, A. Khokhlov, M. Kunst, P. Kuznetsov, H. Meller, O. Mochalov, V. Moiseyev, N. Nicklisch, S. L. Pichler, R. Risch, M. A. Rojo Guerra, C. Roth, A. Szécsényi-Nagy, J. Wahl, M. Meyer, J. Krause, D. Brown, D. Anthony, A. Cooper, K. W. Alt, D. Reich, Massive migration from the steppe was a source for Indo-European languages in Europe. Nature 552, 207–211 (2015).
6
M. Lipson, A. Szécsényi-Nagy, S. Mallick, A. Pósa, B. Stégmár, V. Keerl, N. Rohland, K. Stewardson, M. Ferry, M. Michel, J. Oppenheimer, N. Broomandkhoshbacht, E. Harney, S. Nordenfelt, B. Llamas, B. G. Mende, K. Köhler, K. Oross, M. Bondár, T. Marton, A. Osztás, J. Jakucs, T. Paluch, F. Horváth, P. Csengeri, J. Koós, K. Sebők, A. Anders, P. Raczky, J. Regenye, J. P. Barna, S. Fábián, G. Serlegi, Z. Toldi, E. G. Nagy, J. Dani, E. Molnár, G. Pálfi, L. Márk, B. Melegh, Z. Bánfai, L. Domboróczki, J. Fernández-Eraso, J. A. Mujika-Alustiza, C. A. Fernández, J. J. Echevarría, R. Bollongino, J. Orschiedt, K. Schierhold, H. Meller, A. Cooper, J. Burger, E. Bánffy, K. W. Alt, C. Lalueza-Fox, W. Haak, D. Reich, Parallel palaeogenomic transects reveal complex genetic history of early European farmers. Nature 551, 368–372 (2017).
7
I. Mathieson, S. Alpaslan-Roodenberg, C. Posth, A. Szécsényi-Nagy, N. Rohland, S. Mallick, I. Olalde, N. Broomandkhoshbacht, F. Candilio, O. Cheronet, D. Fernandes, M. Ferry, B. Gamarra, G. G. Fortes, W. Haak, E. Harney, E. Jones, D. Keating, B. Krause-Kyora, I. Kucukkalipci, M. Michel, A. Mittnik, K. Nägele, M. Novak, J. Oppenheimer, N. Patterson, S. Pfrengle, K. Sirak, K. Stewardson, S. Vai, S. Alexandrov, K. W. Alt, R. Andreescu, D. Antonović, A. Ash, N. Atanassova, K. Bacvarov, M. B. Gusztáv, H. Bocherens, M. Bolus, A. Boroneanţ, Y. Boyadzhiev, A. Budnik, J. Burmaz, S. Chohadzhiev, N. J. Conard, R. Cottiaux, M. Čuka, C. Cupillard, D. G. Drucker, N. Elenski, M. Francken, B. Galabova, G. Ganetsovski, B. Gély, T. Hajdu, V. Handzhyiska, K. Harvati, T. Higham, S. Iliev, I. Janković, I. Karavanić, D. J. Kennett, D. Komšo, A. Kozak, D. Labuda, M. Lari, C. Lazar, M. Leppek, K. Leshtakov, D. L. Vetro, D. Los, I. Lozanov, M. Malina, F. Martini, K. McSweeney, H. Meller, M. Menđušić, P. Mirea, V. Moiseyev, V. Petrova, T. D. Price, A. Simalcsik, L. Sineo, M. Šlaus, V. Slavchev, P. Stanev, A. Starović, T. Szeniczey, S. Talamo, M. Teschler-Nicola, C. Thevenet, I. Valchev, F. Valentin, S. Vasilyev, F. Veljanovska, S. Venelinova, E. Veselovskaya, B. Viola, C. Virag, J. Zaninović, S. Zäuner, P. W. Stockhammer, G. Catalano, R. Krauß, D. Caramelli, G. Zariņa, B. Gaydarska, M. Lillie, A. G. Nikitin, I. Potekhina, A. Papathanasiou, D. Borić, C. Bonsall, J. Krause, R. Pinhasi, D. Reich, The genomic history of southeastern Europe. Nature 555, 197–203 (2018).
8
I. Mathieson, I. Lazaridis, N. Rohland, S. Mallick, N. Patterson, S. A. Roodenberg, E. Harney, K. Stewardson, D. Fernandes, M. Novak, K. Sirak, C. Gamba, E. R. Jones, B. Llamas, S. Dryomov, J. Pickrell, J. L. Arsuaga, J. M. Bermúdez de Castro, E. Carbonell, F. Gerritsen, A. Khokhlov, P. Kuznetsov, M. Lozano, H. Meller, O. Mochalov, V. Moiseyev, M. A. Rojo Guerra, J. Roodenberg, J. M. Vergès, J. Krause, A. Cooper, K. W. Alt, D. Brown, D. Anthony, C. Lalueza-Fox, W. Haak, R. Pinhasi, D. Reich, Genome-wide patterns of selection in 230 ancient Eurasians. Nature 528, 499–503 (2015).
9
Z. Hofmanová, S. Kreutzer, G. Hellenthal, C. Sell, Y. Diekmann, D. Diez-Del-Molino, L. van Dorp, S. López, A. Kousathanas, V. Link, K. Kirsanow, L. M. Cassidy, R. Martiniano, M. Strobel, A. Scheu, K. Kotsakis, P. Halstead, S. Triantaphyllou, N. Kyparissi-Apostolika, D.-C. Urem-Kotsou, C. Ziota, F. Adaktylou, S. Gopalan, D. M. Bobo, L. Winkelbach, J. Blöcher, M. Unterländer, C. Leuenberger, Ç. Çilingiroğlu, B. Horejs, F. Gerritsen, S. J. Shennan, D. G. Bradley, M. Currat, K. R. Veeramah, D. Wegmann, M. G. Thomas, C. Papageorgopoulou, J. Burger, Early farmers from across Europe directly descended from Neolithic Aegeans. Proc. Natl. Acad. Sci. U.S.A. 113, 6886–6891 (2016).
10
C. Gamba, E. R. Jones, M. D. Teasdale, R. L. McLaughlin, G. Gonzalez-Fortes, V. Mattiangeli, L. Domboróczki, I. Kővári, I. Pap, A. Anders, A. Whittle, J. Dani, P. Raczky, T. F. G. Higham, M. Hofreiter, D. G. Bradley, R. Pinhasi, Genome flux and stasis in a five millennium transect of European prehistory. Nat. Commun. 5, 5257 (2014).
11
C. Valdiosera, T. Günther, J. C. Vera-Rodríguez, I. Ureña, E. Iriarte, R. Rodríguez-Varela, L. G. Simões, R. M. Martínez-Sánchez, E. M. Svensson, H. Malmström, L. Rodríguez, J.-M. Bermúdez de Castro, E. Carbonell, A. Alday, J. A. Hernández Vera, A. Götherström, J.-M. Carretero, J. L. Arsuaga, C. I. Smith, M. Jakobsson, Four millennia of Iberian biomolecular prehistory illustrate the impact of prehistoric migrations at the far end of Eurasia. Proc. Natl. Acad. Sci. U.S.A. 115, 3428–3433 (2018).
12
I. Olalde, S. Mallick, N. Patterson, N. Rohland, V. Villalba-Mouco, M. Silva, K. Dulias, C. J. Edwards, F. Gandini, M. Pala, P. Soares, M. Ferrando-Bernal, N. Adamski, N. Broomandkhoshbacht, O. Cheronet, B. J. Culleton, D. Fernandes, A. M. Lawson, M. Mah, J. Oppenheimer, K. Stewardson, Z. Zhang, J. M. Jiménez Arenas, I. J. Toro Moyano, D. C. Salazar-García, P. Castanyer, M. Santos, J. Tremoleda, M. Lozano, P. García Borja, J. Fernández-Eraso, J. A. Mujika-Alustiza, C. Barroso, F. J. Bermúdez, E. Viguera Mínguez, J. Burch, N. Coromina, D. Vivó, A. Cebrià, J. Maria Fullola, O. García-Puchol, J. I. Morales, F. X. Oms, T. Majó, J. M. Vergès, A. Díaz-Carvajal, I. Ollich-Castanyer, F. J. López-Cachero, A. M. Silva, C. Alonso-Fernández, G. Delibes de Castro, J. Jiménez Echevarría, A. Moreno-Márquez, G. P. Berlanga, P. Ramos-García, J. Ramos-Muñoz, E. Vijande Vila, G. Aguilella Arzo, Á. Esparza Arroyo, K. T. Lillios, J. Mack, J. Velasco-Vázquez, A. Waterman, L. Benitez de Lugo Enrich, M. Benito Sánchez, B. Agustí, F. Codina, G. de Prado, A. Estalrrich, Á. Fernández Flores, C. Finlayson, G. Finlayson, S. Finlayson, F. Giles-Guzmán, A. Rosas, V. Barciela González, G. García Atiénzar, M. S. Hernández Pérez, A. Llanos, Y. Carrión Marco, I. Collado Beneyto, D. López-Serrano, M. S. Tormo, A. C. Valera, C. Blasco, C. Liesau, P. Ríos, J. Daura, M. Jesús de Pedro Michó, A. A. Diez-Castillo, R. Flores Fernández, J. Francès Farré, R. Garrido-Pena, V. S. Gonçalves, E. Guerra-Doce, A. Merino Herrero-Corral, J. Juan-Cabanilles, D. López-Reyes, S. B. McClure, M. M. Pérez, A. O. Foix, M. Sanz Borràs, A. C. Sousa, J. M. Vidal Encinas, D. J. Kennett, M. B. Richards, K. W. Alt, W. Haak, R. Pinhasi, C. Lalueza-Fox, D. Reich, The genomic history of the Iberian Peninsula over the past 8000 years. Science 363, 1230–1234 (2019).
13
R. Bollongino, O. Nehlich, M. P. Richards, J. Orschiedt, M. G. Thomas, C. Sell, Z. Fajkošová, A. Powell, J. Burger, 2000 years of parallel societies in Stone Age Central Europe. Science 342, 479–481 (2013).
14
E. Bánffy, The 6th Millennium BC Boundary in Western Transdanubia and Its Role in the Central European Neolithic Transition: The Szentgyörgyvölgy-Pityerdomb Settlement (Archaeological Institute of the Hungarian Academy of Sciences, 2004).
15
V. Villalba-Mouco, M. S. van de Loosdrecht, C. Posth, R. Mora, J. Martínez-Moreno, M. Rojo-Guerra, D. C. Salazar-García, J. I. Royo-Guillén, M. Kunst, H. Rougier, I. Crevecoeur, H. Arcusa-Magallón, C. Tejedor-Rodríguez, I. G.-M. de Lagrán, R. Garrido-Pena, K. W. Alt, C. Jeong, S. Schiffels, P. Utrilla, J. Krause, W. Haak, Survival of Late Pleistocene Hunter-Gatherer Ancestry in the Iberian Peninsula. Curr. Biol. 29, 1169–1177.e7 (2019).
16
D. Hofmann, The changing role of La Hoguette pottery in an LBK context, in Something Out of the Ordinary? Interpreting Diversity in the Early Neolithic Linearbandkeramik and Beyond, L. Amkreutz, F. Haack, D. Hofmann, I. van Wijk, Eds. (CSP, 2016).
17
B. Gehlen, No future? No past? Mesolithic heritage in neolithic burials, in Mesolithic Burials—Rites, Symbols and Social Organisation of Early Postglacial Communities. Mesolithische Bestattungen—Riten, Symbole Und Soziale Organisation Früher Postglazialer Gemeinschaften, J. M. Grünberg, B. Gramsch, L. Larsson, J. Orschiedt, H. Meller, Eds. (Tagungen des Landesmuseums für Vorgeschichte Halle 13. Halle/Saale: Landesamt f. Denkmalpflege u. Archäologie Sachsen-Anhalt, 2017), pp. 827–849.
18
D. Binder, Mesolithic and Neolithic interaction in southern France and northern Italy: New data and current hypotheses, in Europe's First Farmers, T. D. Price, Ed. (Cambridge Univ. Press, 2000), pp. 117–143.
19
D. Binder, J. Battentier, C. Delhon, I. Sénépart, In pursuit of a missing transition : The Mesolithic and Neolithic radiocarbon chronology at la Font-aux-Pigeons rock shelter. Antiquity 91, 605–620 (2017).
20
O. García-Puchol, S. B. McClure, J. Juan-Cabanilles, A. A. Diez-Castillo, J. Bernabeu-Aubán, B. Martí-Oliver, S. Pardo-Gordó, J. L. Pascual-Benito, M. Pérez-Ripoll, L. Molina-Balaguer, D. J. Kennett, Cocina cave revisited: Bayesian radiocarbon chronology for the last hunter-gatherers and first farmers in Eastern Iberia. Quat. Int. 472, 259–271 (2018).
21
P.-Y. Nicod, T. Perrin, F.-X. Le Bourdonnec, S. Philibert, C. Oberlin, M. Besse, First Obsidian in the Northern French Alps during the Early Neolithic. J. Field Archaeol. 44, 180–194 (2019).
22
D. Binder, P. Lanos, L. Angeli, L. Gomart, J. Guilaine, C. Manen, R. Maggi, I. M. Muntoni, C. Panelli, G. Radi, C. Tozzi, D. Arobba, J. Battentier, M. Brandaglia, L. Bouby, F. Briois, A. Carré, C. Delhon, L. Gourichon, P. Marinval, R. Nisbet, S. Rossi, P. Rowley-Conwy, S. Thiébault, Modelling the earliest north-western dispersal of Mediterranean Impressed Wares: New dates and Bayesian chronological model. Documenta Praehistorica 44, 54–77 (2017).
23
A. Denaire, P. Lefranc, J. Wahl, C. Bronk Ramsey, E. Dunbar, T. Goslar, A. Bayliss, N. Beaven, P. Bickle, A. Whittle, The cultural project : Formal chronological modelling of the Early and Middle Neolithic sequence in Lower Alsace. J. Archaeol. Method Theory 24, 1072–1149 (2017).
24
C. Jeunesse, S. van Willigen, Westmediterranes Frühneolithikum und westliche Linearbandkeramik: Impulse, Interaktionen, Mischkulturen, in Die Neolithisierung Mitteleuropas—The Spread of the Neolithic to Central Europe, D. Gronenborn, J. Petrasch, Eds. (RGZM-Tagungen, 2010), vol. 4, pp. 569–606.
25
L. Laporte, Néolithisation de la façade atlantique du Centre-Ouest de la France, in Unité et Diversité des Processus de Néolithisation (Mémoires XXXVI de la Société Préhistorique Française, Actes du Colloque de Nantes, 2005), pp. 99–125.
26
M. Rivollat, F. Mendisco, M.-H. Pemonge, A. Safi, D. Saint-Marc, A. Brémond, C. Couture-Veschambre, S. Rottier, M.-F. Deguilloux, When the Waves of European Neolithization Met: First Paleogenetic Evidence from Early Farmers in the Southern Paris Basin. PLOS ONE 10, e0125521 (2015).
27
T. Maricic, M. Whitten, S. Pääbo, Multiplexed DNA Sequence Capture of Mitochondrial Genomes Using PCR Products. PLOS ONE 5, e14004 (2010).
28
Q. Fu, A. Mittnik, P. L. F. Johnson, K. Bos, M. Lari, R. Bollongino, C. Sun, L. Giemsch, R. Schmitz, J. Burger, A. M. Ronchitelli, F. Martini, R. G. Cremonesi, J. Svoboda, P. Bauer, D. Caramelli, S. Castellano, D. Reich, S. Pääbo, J. Krause, A Revised Timescale for Human Evolution Based on Ancient Mitochondrial Genomes. Curr. Biol. 23, 553–559 (2013).
29
T. S. Korneliussen, A. Albrechtsen, R. Nielsen, ANGSD: Analysis of Next Generation Sequencing Data. BMC Bioinform. 15, 356 (2014).
30
J. M. Monroy Kuhn, M. Jakobsson, T. Günther, Estimating genetic kin relationships in prehistoric populations. PLOS ONE 13, e0195491 (2018).
31
N. Patterson, P. Moorjani, Y. Luo, S. Mallick, N. Rohland, Y. Zhan, T. Genschoreck, T. Webster, D. Reich, Ancient admixture in human history. Genetics 192, 1065–1093 (2012).
32
I. Lazaridis, D. Nadel, G. Rollefson, D. C. Merrett, N. Rohland, S. Mallick, D. Fernandes, M. Novak, B. Gamarra, K. Sirak, S. Connell, K. Stewardson, E. Harney, Q. Fu, G. Gonzalez-Fortes, E. R. Jones, S. A. Roodenberg, G. Lengyel, F. Bocquentin, B. Gasparian, J. M. Monge, M. Gregg, V. Eshed, A.-S. Mizrahi, C. Meiklejohn, F. Gerritsen, L. Bejenaru, M. Blüher, A. Campbell, G. Cavalleri, D. Comas, P. Froguel, E. Gilbert, S. M. Kerr, P. Kovacs, J. Krause, D. McGettigan, M. Merrigan, D. A. Merriwether, S. O'Reilly, M. B. Richards, O. Semino, M. Shamoon-Pour, G. Stefanescu, M. Stumvoll, A. Tönjes, A. Torroni, J. F. Wilson, L. Yengo, N. A. Hovhannisyan, N. Patterson, R. Pinhasi, D. Reich, Genomic insights into the origin of farming in the ancient Near East. Nature 536, 419–424 (2016).
33
Q. Fu, C. Posth, M. Hajdinjak, M. Petr, S. Mallick, D. Fernandes, A. Furtwängler, W. Haak, M. Meyer, A. Mittnik, B. Nickel, A. Peltzer, N. Rohland, V. Slon, S. Talamo, I. Lazaridis, M. Lipson, I. Mathieson, S. Schiffels, P. Skoglund, A. P. Derevianko, N. Drozdov, V. Slavinsky, A. Tsybankov, R. Grifoni Cremonesi, F. Mallegni, B. Gély, E. Vacca, M. R. Gonzáles Morales, L. G. Straus, C. Neugebauer-Maresch, M. Teschler-Nicola, S. Constantin, O. T. Moldovan, S. Benazzi, M. Peresani, D. Coppola, M. Lari, S. Ricci, A. Ronchitelli, F. Valentin, C. Thevenet, K. Wehrberger, D. Grigorescu, H. Rougier, I. Crevecoeur, D. Flas, P. Semal, M. A. Mannino, C. Cupillard, H. Bocherens, N. J. Conard, K. Harvati, V. Moiseyev, D. G. Drucker, J. Svoboda, M. P. Richards, D. Caramelli, R. Pinhasi, J. Kelso, N. Patterson, J. Krause, S. Pääbo, D. Reich, The genetic history of Ice Age Europe. Nature 534, 200–205 (2016).
34
I. Lazaridis, N. Patterson, A. Mittnik, G. Renaud, S. Mallick, K. Kirsanow, P. H. Sudmant, J. G. Schraiber, S. Castellano, M. Lipson, B. Berger, C. Economou, R. Bollongino, Q. Fu, K. I. Bos, S. Nordenfelt, H. Li, C. de Filippo, K. Prüfer, S. Sawyer, C. Posth, W. Haak, F. Hallgren, E. Fornander, N. Rohland, D. Delsate, M. Francken, J.-M. Guinet, J. Wahl, G. Ayodo, H. A. Babiker, G. Bailliet, E. Balanovska, O. Balanovsky, R. Barrantes, G. Bedoya, H. Ben-Ami, J. Bene, F. Berrada, C. M. Bravi, F. Brisighelli, G. B. J. Busby, F. Cali, M. Churnosov, D. E. C. Cole, D. Corach, L. Damba, G. van Driem, S. Dryomov, J.-M. Dugoujon, S. A. Fedorova, I. Gallego Romero, M. Gubina, M. Hammer, B. M. Henn, T. Hervig, U. Hodoglugil, A. R. Jha, S. Karachanak-Yankova, R. Khusainova, E. Khusnutdinova, R. Kittles, T. Kivisild, W. Klitz, V. Kučinskas, A. Kushniarevich, L. Laredj, S. Litvinov, T. Loukidis, R. W. Mahley, B. Melegh, E. Metspalu, J. Molina, J. Mountain, K. Näkkäläjärvi, D. Nesheva, T. Nyambo, L. Osipova, J. Parik, F. Platonov, O. Posukh, V. Romano, F. Rothhammer, I. Rudan, R. Ruizbakiev, H. Sahakyan, A. Sajantila, A. Salas, E. B. Starikovskaya, A. Tarekegn, D. Toncheva, S. Turdikulova, I. Uktveryte, O. Utevska, R. Vasquez, M. Villena, M. Voevoda, C. A. Winkler, L. Yepiskoposyan, P. Zalloua, T. Zemunik, A. Cooper, C. Capelli, M. G. Thomas, A. Ruiz-Linares, S. A. Tishkoff, L. Singh, K. Thangaraj, R. Villems, D. Comas, R. Sukernik, M. Metspalu, M. Meyer, E. E. Eichler, J. Burger, M. Slatkin, S. Pääbo, J. Kelso, D. Reich, J. Krause, Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature 513, 409–413 (2014).
35
D. M. Fernandes, D. Strapagiel, P. Borówka, B. Marciniak, E. Żądzińska, K. Sirak, V. Siska, R. Grygiel, J. Carlsson, A. Manica, W. Lorkiewicz, R. Pinhasi, A genomic Neolithic time transect of hunter-farmer admixture in central Poland. Sci. Rep. 8, 14879 (2018).
36
C. Jeunesse, Pour une origine occidentale de la culture de Michelsberg? Materialhefte zur Archäologie in Baden-Württemberg 43, 29–45 (1998).
37
S. Hartz, H. Lübke, New evidence for a chronostratigraphic division of the Ertebølle culture and the earliest funnel beaker culture on the Southern Mecklenburg Bay, in After The Ice Age. Settlements, Subsistence and Social Development in the Mesolithic of Central Europe (Konrad Theiss Verlag, 2006) pp. 59–74.
38
P. Skoglund, H. Malmström, A. Omrak, M. Raghavan, C. Valdiosera, T. Günther, P. Hall, K. Tambets, J. Parik, K.-G. Sjögren, J. Apel, E. Willerslev, J. Storå, A. Götherström, M. Jakobsson, Genomic Diversity and Admixture Differs for Stone-Age Scandinavian Foragers and Farmers. Science 344, 747–750 (2014).
39
I. Olalde, S. Brace, M. E. Allentoft, I. Armit, K. Kristiansen, T. Booth, N. Rohland, S. Mallick, A. Szécsényi-Nagy, A. Mittnik, E. Altena, M. Lipson, I. Lazaridis, T. K. Harper, N. Patterson, N. Broomandkhoshbacht, Y. Diekmann, Z. Faltyskova, D. Fernandes, M. Ferry, E. Harney, P. de Knijff, M. Michel, J. Oppenheimer, K. Stewardson, A. Barclay, K. W. Alt, C. Liesau, P. Ríos, C. Blasco, J. V. Miguel, R. M. García, A. A. Fernández, E. Bánffy, M. Bernabò-Brea, D. Billoin, C. Bonsall, L. Bonsall, T. Allen, L. Büster, S. Carver, L. C. Navarro, O. E. Craig, G. T. Cook, B. Cunliffe, A. Denaire, K. E. Dinwiddy, N. Dodwell, M. Ernée, C. Evans, M. Kuchařík, J. F. Farré, C. Fowler, M. Gazenbeek, R. G. Pena, M. Haber-Uriarte, E. Haduch, G. Hey, N. Jowett, T. Knowles, K. Massy, S. Pfrengle, P. Lefranc, O. Lemercier, A. Lefebvre, C. H. Martínez, V. G. Olmo, A. B. Ramírez, J. L. Maurandi, T. Majó, J. I. McKinley, K. McSweeney, B. G. Mende, A. Modi, G. Kulcsár, V. Kiss, A. Czene, R. Patay, A. Endrődi, K. Köhler, T. Hajdu, T. Szeniczey, J. Dani, Z. Bernert, M. Hoole, O. Cheronet, D. Keating, P. Velemínský, M. Dobeš, F. Candilio, F. Brown, R. F. Fernández, A.-M. Herrero-Corral, S. Tusa, E. Carnieri, L. Lentini, A. Valenti, A. Zanini, C. Waddington, G. Delibes, E. Guerra-Doce, B. Neil, M. Brittain, M. Luke, R. Mortimer, J. Desideri, M. Besse, G. Brücken, M. Furmanek, A. Hałuszko, M. Mackiewicz, A. Rapiński, S. Leach, I. Soriano, K. T. Lillios, J. L. Cardoso, M. P. Pearson, P. Włodarczak, T. D. Price, P. Prieto, P.-J. Rey, R. Risch, M. A. Rojo Guerra, A. Schmitt, J. Serralongue, A. M. Silva, V. Smrčka, L. Vergnaud, J. Zilhão, D. Caramelli, T. Higham, M. G. Thomas, D. J. Kennett, H. Fokkens, V. Heyd, A. Sheridan, K.-G. Sjögren, P. W. Stockhammer, J. Krause, R. Pinhasi, W. Haak, I. Barnes, C. Lalueza-Fox, D. Reich, The Beaker phenomenon and the genomic transformation of northwest Europe. Nature 555, 190–196 (2018).
40
S. Brace, Y. Diekmann, T. J. Booth, L. van Dorp, Z. Faltyskova, N. Rohland, S. Mallick, I. Olalde, M. Ferry, M. Michel, J. Oppenheimer, N. Broomandkhoshbacht, K. Stewardson, R. Martiniano, S. Walsh, M. Kayser, S. Charlton, G. Hellenthal, I. Armit, R. Schulting, O. E. Craig, A. Sheridan, M. Parker Pearson, C. Stringer, D. Reich, M. G. Thomas, I. Barnes, Ancient genomes indicate population replacement in Early Neolithic Britain. Nat. Ecol. Evol. 3, 765–771 (2019).
41
L. Gomart, A. Weiner, M. Gabriele, G. Durrenmath, S. Sorin, L. Angeli, M. Colombo, C. Fabbri, R. Maggi, C. Panelli, D. F. Pisani, G. Radi, C. Tozzi, D. Binder, Spiralled patchwork in pottery manufacture and the introduction of farming to Southern Europe. Antiquity 91, 1501–1514 (2017).
42
M. Lichardus-Itten, Premières influences méditerranéennes dans le Néolithique du Bassin parisien, contribution au débat, in Le Néolithique de la France. Hommage à Gérard Bailloud, J.-P. Demoule, J. Guilaine, Eds. (Picard, 1986), pp. 147–160.
43
G. Marchand, Neolithic fragrances: Mesolithic-Neolithic interactions in western France, in Proceedings-British Academy (Oxford Univ. Press Inc., 2007), vol. 144, p. 225.
44
A. Sheridan, The neolithisation of Britain and Ireland: The big picture, in Landscapes in Transition, B. Finlayson, G. Warren, Eds. (Oxbow, 2010), pp. 89–105.
45
J. Dabney, M. Knapp, I. Glocke, M.-T. Gansauge, A. Weihmann, B. Nickel, C. Valdiosera, N. García, S. Pääbo, J.-L. Arsuaga, M. Meyer, Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. Proc. Natl. Acad. Sci. U.S.A. 110, 15758–15763 (2013).
46
M. Meyer, M. Kircher, Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb. Protoc. 6, pdb.prot5448 (2010).
47
M. Kircher, S. Sawyer, M. Meyer, Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform. Nucleic Acids Res. 40, e3 (2012).
48
N. Rohland, E. Harney, S. Mallick, S. Nordenfelt, D. Reich, Partial uracil–DNA–glycosylase treatment for screening of ancient DNA. Philos. Trans. R. Soc. B Biol. Sci. 370, 20130624 (2015).
49
A. Peltzer, G. Jäger, A. Herbig, A. Seitz, C. Kniep, J. Krause, K. Nieselt, EAGER: Efficient ancient genome reconstruction. Genome Biol. 17, 60 (2016).
50
M. Schubert, S. Lindgreen, L. Orlando, AdapterRemoval v2: Rapid adapter trimming, identification, and read merging. BMC. Res. Notes 9, 88 (2016).
51
H. Li, R. Durbin, Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics 26, 589–595 (2010).
52
A. Mittnik, C.-C. Wang, J. Svoboda, J. Krause, A molecular approach to the sexing of the triple burial at the Upper Paleolithic site of Dolní Věstonice. PLOS ONE 11, e0163019 (2016).
53
Q. Fu, H. Li, P. Moorjani, F. Jay, S. M. Slepchenko, A. A. Bondarev, P. L. F. Johnson, A. Aximu-Petri, K. Prüfer, C. de Filippo, M. Meyer, N. Zwyns, D. C. Salazar-García, Y. V. Kuzmin, S. G. Keates, P. A. Kosintsev, D. I. Razhev, M. P. Richards, N. V. Peristov, M. Lachmann, K. Douka, T. F. G. Higham, M. Slatkin, J.-J. Hublin, D. Reich, J. Kelso, T. B. Viola, S. Pääbo, Genome sequence of a 45,000-year-old modern human from western Siberia. Nature 514, 445–449 (2014).
54
H. Li, B. Handsaker, A. Wysoker, T. Fennell, J. Ruan, N. Homer, G. Marth, G. Abecasis, R. Durbin; 1000 Genome Project Data Processing Subgroup, The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
55
M. Kearse, R. Moir, A. Wilson, S. Stones-Havas, M. Cheung, S. Sturrock, S. Buxton, A. Cooper, S. Markowitz, C. Duran, T. Thierer, B. Ashton, P. Meintjes, A. Drummond, Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28, 1647–1649 (2012).
56
H. Weissensteiner, D. Pacher, A. Kloss-Brandstätter, L. Forer, G. Specht, H.-J. Bandelt, F. Kronenberg, A. Salas, S. Schönherr, HaploGrep 2: Mitochondrial haplogroup classification in the era of high-throughput sequencing. Nucleic Acids Res. 44, W58–W63 (2016).
57
M. Morgan, H. Pagès, V. Obenchain, N. Hayden, Rsamtools: Binary alignment (BAM), FASTA, variant call (BCF), and tabix file import, R package version 1.34.1 (2019).
58
S. Mallick, H. Li, M. Lipson, I. Mathieson, M. Gymrek, F. Racimo, M. Zhao, N. Chennagiri, S. Nordenfelt, A. Tandon, P. Skoglund, I. Lazaridis, S. Sankararaman, Q. Fu, N. Rohland, G. Renaud, Y. Erlich, T. Willems, C. Gallo, J. P. Spence, Y. S. Song, G. Poletti, F. Balloux, G. van Driem, P. de Knijff, I. G. Romero, A. R. Jha, D. M. Behar, C. M. Bravi, C. Capelli, T. Hervig, A. Moreno-Estrada, O. L. Posukh, E. Balanovska, O. Balanovsky, S. Karachanak-Yankova, H. Sahakyan, D. Toncheva, L. Yepiskoposyan, C. Tyler-Smith, Y. Xue, M. S. Abdullah, A. Ruiz-Linares, C. M. Beall, A. Di Rienzo, C. Jeong, E. B. Starikovskaya, E. Metspalu, J. Parik, R. Villems, B. M. Henn, U. Hodoglugil, R. Mahley, A. Sajantila, G. Stamatoyannopoulos, J. T. S. Wee, R. Khusainova, E. Khusnutdinova, S. Litvinov, G. Ayodo, D. Comas, M. F. Hammer, T. Kivisild, W. Klitz, C. A. Winkler, D. Labuda, M. Bamshad, L. B. Jorde, S. A. Tishkoff, W. S. Watkins, M. Metspalu, S. Dryomov, R. Sukernik, L. Singh, K. Thangaraj, S. Pääbo, J. Kelso, N. Patterson, D. Reich, The Simons Genome Diversity Project: 300 genomes from 142 diverse populations. Nature 538, 201–206 (2016).
59
N. Patterson, A. L. Price, D. Reich, Population Structure and Eigenanalysis. PLOS Genet. 2, e190 (2006).
60
G. R. Warnes, B. Bolker, L. Bonebakker, R. Gentleman, W. Huber, A. Liaw, T. Lumley, M. Maechler, A. Magnusson, S. Moeller, M. Schwartz, B. Venables, gplots: Various R programming tools for plotting data, R package version 2.4.1 (2019).
61
D. Gronenborn, A variation on a basic theme: The transition to farming in southern central Europe. J. World Prehist. 13, 123–210 (1999).
62
G. Marchand, Préhistoire Atlantique. Fonctionnement et Évolution des Sociétés du Paléolithique au Néolithique (Errance, 2014), pp. 520.
63
T. Perrin, D. Binder, Le Mésolithique à trapèzes et la néolithisation de l'Europe sud-occidentale, in La transition Néolithique en Méditerranée. Actes du Colloque Transitions en Méditerranée, ou Comment des Chasseurs Devinrent Agriculteurs, Muséum de Toulouse (Errance et Archives d'Ecologie Préhistorique, 2014), pp. 271–281.
64
G. Marchand, T. Perrin, Why this revolution? Explaining the major technical shift in Southwestern Europe during the 7th millennium cal. BC. Quatern. Int. 428, 73–85 (2017).
65
V. G. Childe, The Dawn of European Civilization (Kegan Paul, 1925), vol. 674, 328 pp.
66
P. Biagi, D. Kiosak, The Mesolithic of the northwestern Pontic region. New AMS dates for the origin and spread of the blade and trapeze industries in southeastern Europe. Eurasia Antiqua 16, 21–41 (2010).
67
C. Perlès, A. Quiles, H. Valladas, Early seventh-millennium AMS dates from domestic seeds in the Initial Neolithic at Franchthi Cave (Argolid, Greece). Antiquity 87, 1001–1015 (2013).
68
N. Efstratiou, A. Karetsou, M. Ntinou, The Neolithic Settlement of Knossos in Crete: New Evidence for the Early Occupation of Crete and the Aegean Islands (INSTAP Academic Press, 2013).
69
A. Reingruber, Preceramic, aceramic or early ceramic? The radiocarbon dated beginning of the Neolithic in the Aegean. Documenta Praehistorica 42, 147–158 (2015).
70
C. Perlès, The Early Neolithic In Greece: The First Farming Communities In Europe (Cambridge Univ. Press, 2001).
71
R. Duru, The Neolithic of the Lakes Region. Hacılar, Kuruçay Höyük, Höyücek, Bademağacı Höyük, in The Neolithic in Turkey. New Excavations and New Research, M. Özdoğan, N. Başgelen, P. Kuniholm, Eds. (Archaeology and Art Publications, 2012), vol. 4, pp. 1–65.
72
F. Gerritsen, R. Özbal, L. Thissen, Barcın Höyük. The beginnings of farming in the Marmara region, in The Neolithic in Turkey. New Excavations and New Research, M. Özdoğan, N. Başgelen, P. Kuniholm, Eds. (Archaeology and Art Publications, 2013), pp. 93–112.
73
F. Klopffleisch, Charakteristik und Zeitfolge der Keramik Mitteldeutschlands, in Die Grabhügel von Leubingen, Sömmerda und Nienstedt, F. Klopffleisch, Ed. (Heft II. Vorgeschichtliche Altertümer der Provinz Sachsen, Erste Abteilung, 1884).
74
C. Jeunesse, La céramique de La Hoguette: un nouvel élément non- rubané du Néolithique ancien de l’Europe du nord-ouest. Cahiers Alsaciens d’Archéologie, d’Art et d’Histoire 30, 3–33 (1987).
75
D. Binder, Mésolithique et Néolithique ancien en Italie et dans le sud-est de la France entre 7000 et 5500 BCE cal: Questions ouvertes sur les dynamiques culturelles et les procès d’interaction, in Transitions, Ruptures et Continuité Durant la Préhistoire, Actes du XXVIIe Congrès Préhistorique de France (Bordeaux - Les Eysies, 2010), T. Perrin, C. Manen, G. Marchand, P. Allard, D. Binder, M. Ilett, Eds. (Société Préhistorique Française, 2013), pp. 341–355.
76
J.-F. Berger, G. Metallinou, J. Guilaine, Vers une révision de la transition méso-néolithique sur le site de Sidari (Corfou, Grèce): Nouvelles données géoarchéologiques et radiocarbone, évaluation des processus post-dépositionnels, in La transition Néolithique en Méditerranée (Actes du Colloque Transitions en Méditerranée ou Comment des Chasseurs Devinrent Agriculteurs, Muséum de Toulouse, 14–15 Avril 2011), C. Manen, T. Perrin, J. Guilaine, Eds. (Errance & Archives d'Ecologie Préhistorique, 2014), pp. 213–232.
77
C. Manen, T. Perrin, A. Raux, D. Binder, F.-X. Le Bourdonnec, F. Briois, F. Convertini, S. Dubernet, G. Escallon, L. Gomart, J. Guilaine, C. Hamon, S. Philibert, A. Queffelec, Le sommet de l’iceberg ? Colonisation pionnière et néolithisation de la France méditerranéenne. Bull. Soc. Préhist. Fr. 116, 317–361 (2019).
78
C. Collina, Le Néolithique Ancien en Italie du Sud, Evolution des Industries Lithiques entre VIIe et VIe Millénaire (Archaeopress, 2015).
79
O. García-Puchol, D. C. Salazar-García, Times of Neolithic Transition Along the Western Mediterranean (Springer, 2017).
80
C. Manen, T. Perrin, J. Guilaine, L. Bouby, S. Bréhard, F. Briois, F. Durand, P. Marinval, J.-D. Vigne, The Neolithic transition in the Western Mediterranean: a complex and non-linear diffusion process. The radiocarbon record revisited. Radiocarbon 61, 531–571 (2018).
81
J. Roussot-Larroque, Rubané et Cardial: Le poids de l'Ouest, in Rubané et Cardial, D. Cahen, M. Otte, Eds. (E.R.A.U.L., 1990), vol. 39, pp. 315–360.
82
L. Laporte, C. Picq, Les occupations néolithiques du vallon des Ouchettes à Plassay. Gallia Préhist. (2002).
83
J. Jakucs, E. Bánffy, K. Oross, V. Voicsek, C. B. Ramsey, E. Dunbar, B. Kromer, A. Bayliss, D. Hofmann, P. Marshall, A. W. Whittle, Between the Vinča and Linearbandkeramik worlds: The diversity of practices and identities in the 54th–53rd centuries cal BC in southwest Hungary and beyond. J. World Prehist. 29, 267–336 (2016).
84
D. Gronenborn, Beyond the models:Neolithisation'in Central Europe. Proc. Br. Acad. 144, 73–98 (2007).
85
J. Preuss, Das Neolithikum in Mitteleuropa (Beier & Beran, 1992).
86
J. Müller, Megaliths and Funnel Beakers: Societies in Change, 4100-2700 BC (Stichting Nederlands Museum voor Anthropologie en Praehistorie, 2011).
87
C. Constantin, Fin du Rubané, Céramique du Limbourg et Post-Rubané: Le Néolithique le Plus Ancien en Bassin Parisien et en Hainaut (BAR, 1985), vol. 273.
88
C. Thevenet, “Des faits aux gestes... des gestes aux sens?: Pratiques funéraires et société durant le Néolithique ancien en Bassin parisien,” thesis (2010).
89
E. Ghesquière, L. Juhel, C. Marcigny, Entre Néolithique ancien et Néolithique moyen en Normandie et dans les îles anglonormandes: Parcours chronologique, in Premiers Néolithiques de l'Ouest. Culture, Réseaux, Échanges des Premières Sociétés Néolithiques à Leur Expansion, C. Billard, M. Legris, Eds. (Presses Universitaires de Rennes, 2007), pp. 117–162.
90
C. Constantin, D. Mordant, D. Simonin, La Culture de Cerny. Nouvelle Économie, Nouvelle Société au Néolithique. Actes du Colloque International de Nemours, 9–10-11 Mai 1994 (Mémoires du Musée de Préhistoire d'Île-de-France, 1997), vol. 6, pp. 740.
91
I. Sidéra, Animaux domestiques, bêtes sauvages et objets en matières animales du Rubané au Michelsberg: De l'économie aux symboles, des techniques à la culture. Gallia Préhist. 42, 107–194 (2000).
92
A. Augereau, L'industrie du Silex du Ve au IVe Millénaire Dans le Sud-Est du Bassin Parisien. Rubané, Villeneuve-Saint-Germain, Cerny et Groupe de Noyen, Documents d'Archéologie Française, Ed. (de la Maison des Sciences de l'Homme, 2004), vol. 97, p. 220.
93
G. Vernet, Le Néolithique ancien (épicardial) et le site Néolithique moyen/ancien sur le site du Brézet à Clermont-Ferrand (Massif central, France). Rencontres Méridionales de Préhistoire Récente, 85–90 (2004).
94
S. van Willigen, C. Jeunesse, Le vase à décor rubané de la Grande Grotte à Cheval-Blanc (Vaucluse): Un objet danubien dans le Néolithique ancien du Midi de la France. Bull. Soc. Préhist. Fr., 603–608 (2006).
95
S. Bonnardin, La parure funéraire au Néolithique ancien dans les bassins parisien et rhénan: Rubané, Hinkelstein et Villeneuve-Saint-Germain (Société Préhistorique Française, 2009), vol. XLIX, p. 322.
96
C. Constantin, D. Vachard, Anneaux d'origine méridionale dans le Rubané récent du Bassin parisien. Bull. Soc. Préhist. Fr. 101, 75–83 (2004).
97
P. Biagi, E. Starnini, Some aspects of the neolithization of the Adriatic region, in Atti della Societa Preistoria a Protostoria Regione Friuli-Venezia Giulia (Trieste, 1999).
98
P. Mazzieri, I siti di Via Spezia (Beneficio e Via Guidorossi) et Pontetaro a Parma nel quadro della cultura diei Vasi a Bocca Quadrata in Emila Occidentale (Università di Pisa, 2010).
99
P. Pétrequin, S. Cassen, M. Errera, L. Klassen, A. Sheridan, A.-M. Pétrequin, Jade. Grandes Haches Alpines du Néolithique Européen. Ve et IVe Millénaires av. J.-C. (Presses Universitaires de Franche-Comté & Centre de Recherche Archéologique de la Vallée de l’Ain, 2012), p. 1524.
100
R. Joussaume, Le Néolithique de l'Aunis et du Poitou Occidental Dans son Cadre Atlantique (Travaux du Laboratoire Anthropologie-Préhistoire-Protohistoire-Quaternaire Armoricains, 1981), p. 625.
101
L. Laporte, L. Jallot, M. Sohn, Mégalithismes en France. Nouveaux acquis et nouvelles perspectives de recherche. Gallia Préhist. 53, 289–334 (2011).
102
M. Porr, K. W. Alt, The burial of Bad Dürrenberg, central Germany: Osteopathology and osteoarchaeology of a Late Mesolithic shaman's grave. Int. J. Osteoarchaeol. 16, 395–406 (2006).
103
J. Grünberg, The Mesolithic burials of the Middle Elbe-Saale region, in Mesolithic Burials—Rites, Symbols and Social Organisation of Early Postglacial Communities, International Conference Halle (Saale), Germany, 18 to 21 September 2013.
104
M. Stecher, J. M. Grünberg, K. W. Alt, Bioarchaeology of the Mesolithic individuals from Bottendorf (Thuringia, Germany), in Mesolithic Burials—Rites, Symbols and Social Organisation of Early Postglacial Communities, International Conference Halle (Saale), Germany, 18 to 21 September 2013.
105
G. Le Bras-Goude, D. Binder, V. Formicola, H. Duday, C. Couture-Veschambre, J.-J. Hublin, M. Richard, Stratégies de subsistance et analyse culturelle de populations néolithiques de Ligurie : Approche par l'étude isotopique (δ13 C et δ15 N) des restes osseux. Bull. Mém. Soc. Anthropol. Paris 18, 45–55 (2006).
106
D. Binder, E. J. Brochier, H. Duday, D. Helmer, P. Marinval, S. Thiebault, J. Wattez, L'abri Pendimoun à Castellar (Alpes-Maritimes). Nouvelles données sur le complexe culturel de la céramique imprimée méditerranéenne dans son contexte stratigraphique. Gallia Préhist. 35, 177–251 (1993).
107
A. Zemour, D. Binder, A. Coppa, H. Duday, La sépulture au début du Néolithique (6e et 1ère moitié du 5e millénaire cal-BC) en France et en Italie: de l'idée d'une" simple" fosse à une typologie architecturale. Bull. Mem. Soc. Anthropol. Paris, 1–18 (2016).
108
S. Provost, D. Binder, H. Duday, G. Durrenmath, G. Goude, L. Gourichon, C. Delhon, I. Gentile, M. Vuillien, A. Zemour, Une sépulture collective à la transition des VIe et Ve millénaires BCE: Mougins–Les Bréguières (Alpes-Maritimes, France). Fouilles Maurice Sechter 1966–1967. Gallia Préhist., 289–336 (2017).
109
T. Price, J. Wahl, C. Knipper, E. Burger-Heinrich, G. Kurtz, R. Bentley, Das bandkeramische Gräberfeld vom „Viesenhäuser Hof“ bei Stuttgart-Mühlhausen: Neue Untersuchungsergebnisse zum Migrationsverhalten im frühen Neolithikum, in Fundberichte aus Baden-Württemberg, D. Funda, Ed. (Kommissionsverlag Konrad Theiss Verlag, 2003), pp. 23–58.
110
H. C. Strien, Untersuchungen zur Bandkeramik in Württemberg (Habelt, 2000).
111
A. Hujić, Das Kind in uns unter der Lupe der Isotopie, Allometrie und Pathologie. Zusammenhang zwischen δ15N und δ13C als Eiweißproxy und dem Knochenwachstum bei linienbandkeramischen Individuen aus Stuttgart-Mühlhausen‚ "Viesenhäuser Hof" und Schwetzingen unter Berücksichtigung verschiedener Indikatoren für Nährstoffversorgung (Freie Universität, 2016).
112
C. Gerling, M. Francken, Das linearbandkeramische Gräberfeld von Schwetzingen. Bull. Soc. Suisse Anthropol. 13, 43–50 (2007).
113
R. A. Bentley, T. D. Price, J. Lning, D. Gronenborn, J. Wahl, P. Fullagar, Prehistoric migration in Europe: Strontium isotope analysis of early Neolithic skeletons. Curr. Anthropol. 43, 799–804 (2002).
114
R. A. Bentley, P. Bickle, L. Fibiger, G. M. Nowell, C. W. Dale, R. E. M. Hedges, J. Hamilton, J. Wahl, M. Francken, G. Grupe, E. Lenneis, M. Teschler-Nicola, R.-M. Arbogast, D. Hofmann, A. Whittle, Community differentiation and kinship among Europe’s first farmers. Proc. Natl. Acad. Sci. U.S.A. 109, 9326–9330 (2012).
115
S. Griffiths, Radiocarbon dates from Nitra, Schwetzingen and Vedrovice, in The First Farmers of Central Europe: Diversity In Lbk Lifeways, P. Bickle, A. Whittle, Ed. (Oxbow Books, 2013).
116
B. Fritsch, E. Claßen, U. Müller, V. Dresely, Die linienbandkeramischen Gräberfelder von Derenburg “Meerenstieg II” und Halberstadt “Sonntagsfeld”. Jahressch. Mitteldt. Vorgesch. 92, 25–229 (2011).
117
C. Meyer, C. Knipper, N. Nicklisch, A. Münster, O. Kürbis, V. Dresely, H. Meller, K. W. Alt, Early Neolithic executions indicated by clustered cranial trauma in the mass grave of Halberstadt. Nat. Commun. 9, 2472 (2018).
118
S. Rottier, C. Mordant, P. Chambon, C. Thevenet, Découverte de plus d'une centaine de sépultures du Néolithique moyen à Gurgy, les Noisats (Yonne). Bull. Soc. Préhist. Fr. 102, 641–645 (2005).
119
M. Rivollat, S. Rottier, C. Couture, M.-H. Pemonge, F. Mendisco, M. G. Thomas, M.-F. Deguilloux, P. Gerbault, Investigating mitochondrial DNA relationships in Neolithic Western Europe through serial coalescent simulations. Eur. J. Hum. Genet. 25, 388–392 (2016).
120
E. Ghesquière, P. Chambon, D. Giazzon, L. Hachem, C. Thevenet, A. Thomas, Monumental cemeteries of the 5th millennium BC: The Fleury-sur-Orne contribution, in Proceedings of the International Conference “Megaliths – Societies – Landscapes. Early Monumentality and Social Differentiation in Neolithic Europe” (Kiel, Germany, 16 to 20 June 2015), in J. Müller, M. Hinz, M. Wunderlich, Eds. (Bonn, 2019), pp. 177–191.
121
P. Lefranc, F. Bachellerie, F. Chenal, A. Denaire, C. Féliu, H. Réveillas, N. Schneider, La nécropole néolithique moyen d'Obernai "Neuen Brunnen" (Bas-Rhin) : Rites funéraires de la première moitié du 5e millénaire dans le sud de la plaine du Rhin supérieur (Grossgartach, Planig-Friedberg, Roessen). Rev. Archéol. Est 67, 5–57 (2018).
122
A. Denaire, Chronologie absolue de la séquence Hinkelstein-Grossgartach-Roessen-Bischheim dans le sud de la plaine du Rhin supérieur et le nord de la Franche-Comté à la lumière des dernières données, in Nécropole et Enceintes Danubiennes du Ve Millénaire Dans le Nord-Est de la France et le Sud-Ouest de l’Allemagne, A. Denaire, C. Jeunesse, P. Lefranc, Eds. (Presses universitaires de Strasbourg, 2011), pp. 9–30.
123
C. Jeunesse, Les groupes régionaux occidentaux du Rubané (Rhin et Bassin parisien) à travers les pratiques funéraires. Gallia Préhist. 37, 115–154 (1995).
124
M. Rivollat, H. Réveillas, F. Mendisco, M.-H. Pemonge, P. Justeau, C. Couture, P. Lefranc, C. Féliu, M.-F. Deguilloux, Ancient mitochondrial DNA from the middle neolithic necropolis of Obernai extends the genetic influence of the LBK to west of the Rhine. Am. J. Phys. Anthropol. 161, 522–529 (2016).
125
L. Laporte, R. Joussaume, C. Scarre, Le tumulus C de Péré à Prissé-la-Charrière (Deux-Sèvres). Gallia Préhist. 44, 167–214 (2002).
126
C. Scarre, L. Laporte, R. Jousseaume, Long mounds and megalithic origins in western France: recent excavations at Prissé-la-Charrière. Proc. Prehistoric Soc. 69, 235–251 (2003).
127
M.-F. Deguilloux, L. Soler, M.-H. Pemonge, C. Scarre, R. Joussaume, L. Laporte, News from the west: Ancient DNA from a French megalithic burial chamber. Am. J. Phys. Anthropol. 144, 108–118 (2011).
128
J. Preuss, Das jungsteinzeitliche Koerpergraeberfeld von Tangermuende, Kr.Stendal. Wiss Z Karl Marx Univ Leipzig Ges Sprachwiss 3, 415–482 (1954).
129
F. Lanfranco, A. Kamischke, M. Zitzmann, E. Nieschlag, Klinefelter's syndrome. Lancet 364, 273–283 (2004).
130
S. S. Ebenesersdóttir, M. Sandoval-Velasco, E. D. Gunnarsdóttir, A. Jagadeesan, V. B. Guðmundsdóttir, E. L. Thordardóttir, M. S. Einarsdóttir, K. H. S. Moore, Á. Sigurðsson, D. N. Magnúsdóttir, H. Jónsson, S. Snorradóttir, E. Hovig, P. Møller, I. Kockum, T. Olsson, L. Alfredsson, T. F. Hansen, T. Werge, G. L. Cavalleri, E. Gilbert, C. Lalueza-Fox, J. W. Walser III, S. Kristjánsdóttir, S. Gopalakrishnan, L. Árnadóttir, Ó. Þ. Magnússon, M. T. P. Gilbert, K. Stefánsson, A. Helgason, Ancient genomes from Iceland reveal the making of a human population. Science 360, 1028–1032 (2018).
131
A. Mittnik, C.-C. Wang, S. Pfrengle, M. Daubaras, G. Zariņa, F. Hallgren, R. Allmäe, V. Khartanovich, V. Moiseyev, M. Tõrv, A. Furtwängler, A. Andrades Valtueña, M. Feldman, C. Economou, M. Oinonen, A. Vasks, E. Balanovska, D. Reich, R. Jankauskas, W. Haak, S. Schiffels, J. Krause, The genetic prehistory of the Baltic Sea region. Nat. Commun. 9, 442 (2018).
132
C. Posth, G. Renaud, A. Mittnik, D. G. Drucker, H. Rougier, C. Cupillard, F. Valentin, C. Thevenet, A. Furtwängler, C. Wißing, M. Francken, M. Malina, M. Bolus, M. Lari, E. Gigli, G. Capecchi, I. Crevecoeur, C. Beauval, D. Flas, M. Germonpré, J. van der Plicht, R. Cottiaux, B. Gély, A. Ronchitelli, K. Wehrberger, D. Grigorescu, J. Svoboda, P. Semal, D. Caramelli, H. Bocherens, K. Harvati, N. J. Conard, W. Haak, A. Powell, J. Krause, Pleistocene mitochondrial genomes suggest a single major dispersal of non-Africans and a late glacial population turnover in Europe. Curr. Biol. 26, 827–833 (2016).
133
G. Brandt, W. Haak, C. J. Adler, C. Roth, A. Szécsényi-Nagy, S. Karimnia, S. Möller-Rieker, H. Meller, R. Ganslmeier, S. Friederich, V. Dresely, N. Nicklisch, J. K. Pickrell, F. Sirocko, D. Reich, A. Cooper, K. W. Alt; The Genographic Consortium, Ancient DNA reveals key stages in the formation of central european mitochondrial genetic diversity. Science 342, 257–261 (2013).
134
W. Haak, P. Forster, B. Bramanti, S. Matsumura, G. Brandt, M. Tänzer, R. Villems, C. Renfrew, D. Gronenborn, K. W. Alt, J. Burger, Ancient DNA from the first european farmers in 7500-year-old neolithic sites. Science 310, 1016–1018 (2005).
135
E. R. Jones, G. Gonzalez-Fortes, S. Connell, V. Siska, A. Eriksson, R. Martiniano, R. L. McLaughlin, M. G. Llorente, L. M. Cassidy, C. Gamba, T. Meshveliani, O. Bar-Yosef, W. Müller, A. Belfer-Cohen, Z. Matskevich, N. Jakeli, T. F. G. Higham, M. Currat, D. Lordkipanidze, M. Hofreiter, A. Manica, R. Pinhasi, D. G. Bradley, Upper Palaeolithic genomes reveal deep roots of modern Eurasians. Nat. Commun. 6, 8912 (2015).
136
M. Feldman, E. Fernández-Domínguez, L. Reynolds, D. Baird, J. Pearson, I. Hershkovitz, H. May, N. Goring-Morris, M. Benz, J. Gresky, R. A. Bianco, A. Fairbairn, G. Mustafaoğlu, P. W. Stockhammer, C. Posth, W. Haak, C. Jeong, J. Krause, Late Pleistocene human genome suggests a local origin for the first farmers of central Anatolia. Nat. Commun. 10, 1218 (2019).
137
R. Martiniano, L. M. Cassidy, R. Ó'Maoldúin, R. McLaughlin, N. M. Silva, L. Manco, D. Fidalgo, T. Pereira, M. J. Coelho, M. Serra, J. Burger, R. Parreira, E. Moran, A. C. Valera, E. Porfirio, R. Boaventura, A. M. Silva, D. G. Bradley, The population genomics of archaeological transition in west Iberia: Investigation of ancient substructure using imputation and haplotype-based methods. PLOS Genet. 13, e1006852 (2017).
138
M. van de Loosdrecht, A. Bouzouggar, L. Humphrey, C. Posth, N. Barton, A. Aximu-Petri, B. Nickel, S. Nagel, E. H. Talbi, M. A. El Hajraoui, S. Amzazi, J.-J. Hublin, S. Pääbo, S. Schiffels, M. Meyer, W. Haak, C. Jeong, J. Krause, Pleistocene North African genomes link Near Eastern and sub-Saharan African human populations. Science 360, 548–552 (2018).
139
H. L. Norton, R. A. Kittles, E. Parra, P. McKeigue, X. Mao, K. Cheng, V. A. Canfield, D. G. Bradley, B. McEvoy, M. D. Shriver, Genetic Evidence for the Convergent Evolution of Light Skin in Europeans and East Asians. Mol. Biol. Evol. 24, 710–722 (2006).
140
M. E. Allentoft, M. Sikora, K.-G. Sjögren, S. Rasmussen, M. Rasmussen, J. Stenderup, P. B. Damgaard, H. Schroeder, T. Ahlström, L. Vinner, A.-S. Malaspinas, A. Margaryan, T. Higham, D. Chivall, N. Lynnerup, L. Harvig, J. Baron, P. Della Casa, P. Dąbrowski, P. R. Duffy, A. V. Ebel, A. Epimakhov, K. Frei, M. Furmanek, T. Gralak, A. Gromov, S. Gronkiewicz, G. Grupe, T. Hajdu, R. Jarysz, V. Khartanovich, A. Khokhlov, V. Kiss, J. Kolář, A. Kriiska, I. Lasak, C. Longhi, G. McGlynn, A. Merkevicius, I. Merkyte, M. Metspalu, R. Mkrtchyan, V. Moiseyev, L. Paja, G. Pálfi, D. Pokutta, Ł. Pospieszny, T. D. Price, L. Saag, M. Sablin, N. Shishlina, V. Smrčka, V. I. Soenov, V. Szeverényi, G. Tóth, S. V. Trifanova, L. Varul, M. Vicze, L. Yepiskoposyan, V. Zhitenev, L. Orlando, T. Sicheritz-Pontén, S. Brunak, R. Nielsen, K. Kristiansen, E. Willerslev, Population genomics of Bronze Age Eurasia. Nature 522, 167–172 (2015).
141
J. Burger, M. Kirchner, B. Bramanti, W. Haak, M. G. Thomas, Absence of the lactase-persistence-associated allele in early Neolithic Europeans. Proc. Natl. Acad. Sci. U.S.A. 104, 3736–3741 (2007).
142
T. M. Teslovich, K. Musunuru, A. V. Smith, A. C. Edmondson, I. M. Stylianou, M. Koseki, J. P. Pirruccello, S. Ripatti, D. I. Chasman, C. J. Willer, C. T. Johansen, S. W. Fouchier, A. Isaacs, G. M. Peloso, M. Barbalic, S. L. Ricketts, J. C. Bis, Y. S. Aulchenko, G. Thorleifsson, M. F. Feitosa, J. Chambers, M. Orho-Melander, O. Melander, T. Johnson, X. Li, X. Guo, M. Li, Y. S. Cho, M. J. Go, Y. J. Kim, J.-Y. Lee, T. Park, K. Kim, X. Sim, R. T.-H. Ong, D. C. Croteau-Chonka, L. A. Lange, J. D. Smith, K. Song, J. H. Zhao, X. Yuan, J. A. Luan, C. Lamina, A. Ziegler, W. Zhang, R. Y. L. Zee, A. F. Wright, J. C. M. Witteman, J. F. Wilson, G. Willemsen, H. E. Wichmann, J. B. Whitfield, D. M. Waterworth, N. J. Wareham, G. Waeber, P. Vollenweider, B. F. Voight, V. Vitart, A. G. Uitterlinden, M. Uda, J. Tuomilehto, J. R. Thompson, T. Tanaka, I. Surakka, H. M. Stringham, T. D. Spector, N. Soranzo, J. H. Smit, J. Sinisalo, K. Silander, E. J. G. Sijbrands, A. Scuteri, J. Scott, D. Schlessinger, S. Sanna, V. Salomaa, J. Saharinen, C. Sabatti, A. Ruokonen, I. Rudan, L. M. Rose, R. Roberts, M. Rieder, B. M. Psaty, P. P. Pramstaller, I. Pichler, M. Perola, B. W. J. H. Penninx, N. L. Pedersen, C. Pattaro, A. N. Parker, G. Pare, B. A. Oostra, C. J. O’Donnell, M. S. Nieminen, D. A. Nickerson, G. W. Montgomery, T. Meitinger, R. McPherson, M. I. McCarthy, W. McArdle, D. Masson, N. G. Martin, F. Marroni, M. Mangino, P. K. E. Magnusson, G. Lucas, R. Luben, R. J. F. Loos, M.-L. Lokki, G. Lettre, C. Langenberg, L. J. Launer, E. G. Lakatta, R. Laaksonen, K. O. Kyvik, F. Kronenberg, I. R. König, K.-T. Khaw, J. Kaprio, L. M. Kaplan, Å. Johansson, M.-R. Jarvelin, A. C. J. W. Janssens, E. Ingelsson, W. Igl, G. K. Hovingh, J.-J. Hottenga, A. Hofman, A. A. Hicks, C. Hengstenberg, I. M. Heid, C. Hayward, A. S. Havulinna, N. D. Hastie, T. B. Harris, T. Haritunians, A. S. Hall, U. Gyllensten, C. Guiducci, L. C. Groop, E. Gonzalez, C. Gieger, N. B. Freimer, L. Ferrucci, J. Erdmann, P. Elliott, K. G. Ejebe, A. Döring, A. F. Dominiczak, S. Demissie, P. Deloukas, E. J. C. de Geus, U. de Faire, G. Crawford, F. S. Collins, Y.-D. I. Chen, M. J. Caulfield, H. Campbell, N. P. Burtt, L. L. Bonnycastle, D. I. Boomsma, S. M. Boekholdt, R. N. Bergman, I. Barroso, S. Bandinelli, C. M. Ballantyne, T. L. Assimes, T. Quertermous, D. Altshuler, M. Seielstad, T. Y. Wong, E.-S. Tai, A. B. Feranil, C. W. Kuzawa, L. S. Adair, H. A. Taylor Jr., I. B. Borecki, S. B. Gabriel, J. G. Wilson, H. Holm, U. Thorsteinsdottir, V. Gudnason, R. M. Krauss, K. L. Mohlke, J. M. Ordovas, P. B. Munroe, J. S. Kooner, A. R. Tall, R. A. Hegele, J. J. P. Kastelein, E. E. Schadt, J. I. Rotter, E. Boerwinkle, D. P. Strachan, V. Mooser, K. Stefansson, M. P. Reilly, N. J. Samani, H. Schunkert, L. A. Cupples, M. S. Sandhu, P. M. Ridker, D. J. Rader, C. M. van Duijn, L. Peltonen, G. R. Abecasis, M. Boehnke, S. Kathiresan, Biological, clinical and population relevance of 95 loci for blood lipids. Nature 466, 707–713 (2010).
143
C. D. Huff, D. J. Witherspoon, Y. Zhang, C. Gatenbee, L. A. Denson, S. Kugathasan, H. Hakonarson, A. Whiting, C. T. Davis, W. Wu, J. Xing, W. S. Watkins, M. J. Bamshad, J. P. Bradfield, K. Bulayeva, T. S. Simonson, L. B. Jorde, S. L. Guthery, Crohn’s Disease and Genetic Hitchhiking at IBD5. Mol. Biol. Evol. 29, 101–111 (2011).
144
L. B. Barreiro, M. Ben-Ali, H. Quach, G. Laval, E. Patin, J. K. Pickrell, C. Bouchier, M. Tichit, O. Neyrolles, B. Gicquel, J. R. Kidd, K. K. Kidd, A. Alcaïs, J. Ragimbeau, S. Pellegrini, L. Abel, J.-L. Casanova, L. Quintana-Murci, Evolutionary dynamics of human toll-like receptors and their different contributions to host defense. PLOS Genet. 5, e1000562 (2009).
145
R. Sherva, J. P. Rice, R. J. Neuman, N. Rochberg, N. L. Saccone, L. J. Bierut, Associations and interactions between snps in the alcohol metabolizing genes and alcoholism phenotypes in european americans. Alcohol. Clin. Exp. Res. 33, 848–857 (2009).
146
V. M. Narasimhan, N. Patterson, P. Moorjani, N. Rohland, R. Bernardos, S. Mallick, I. Lazaridis, N. Nakatsuka, I. Olalde, M. Lipson, A. M. Kim, L. M. Olivieri, A. Coppa, M. Vidale, J. Mallory, V. Moiseyev, E. Kitov, J. Monge, N. Adamski, N. Alex, N. Broomandkhoshbacht, F. Candilio, K. Callan, O. Cheronet, B. J. Culleton, M. Ferry, D. Fernandes, S. Freilich, B. Gamarra, D. Gaudio, M. Hajdinjak, É. Harney, T. K. Harper, D. Keating, A. M. Lawson, M. Mah, K. Mandl, M. Michel, M. Novak, J. Oppenheimer, N. Rai, K. Sirak, V. Slon, K. Stewardson, F. Zalzala, Z. Zhang, G. Akhatov, A. N. Bagashev, A. Bagnera, B. Baitanayev, J. Bendezu-Sarmiento, A. A. Bissembaev, G. L. Bonora, T. T. Chargynov, T. Chikisheva, P. K. Dashkovskiy, A. Derevianko, M. Dobeš, K. Douka, N. Dubova, M. N. Duisengali, D. Enshin, A. Epimakhov, A. V. Fribus, D. Fuller, A. Goryachev, A. Gromov, S. P. Grushin, B. Hanks, M. Judd, E. Kazizov, A. Khokhlov, A. P. Krygin, E. Kupriyanova, P. Kuznetsov, D. Luiselli, F. Maksudov, A. M. Mamedov, T. B. Mamirov, C. Meiklejohn, D. C. Merrett, R. Micheli, O. Mochalov, S. Mustafokulov, A. Nayak, D. Pettener, R. Potts, D. Razhev, M. Rykun, S. Sarno, T. M. Savenkova, K. Sikhymbaeva, S. M. Slepchenko, O. A. Soltobaev, N. Stepanova, S. Svyatko, K. Tabaldiev, M. Teschler-Nicola, A. A. Tishkin, V. V. Tkachev, S. Vasilyev, P. Velemínský, D. Voyakin, A. Yermolayeva, M. Zahir, V. S. Zubkov, A. Zubova, V. S. Shinde, C. Lalueza-Fox, M. Meyer, D. Anthony, N. Boivin, K. Thangaraj, D. J. Kennett, M. Frachetti, R. Pinhasi, D. Reich, The formation of human populations in South and Central Asia. Science 365, eaat7487 (2019).
147
G. E. P. Box, D. R. Cox, An analysis of transformations. J. R. Stat. Soc. B. 26, 211–243 (1964).
148
B. D. Ripley, W. N. Venables, Modern Applied Statistics with S (Springer, 2002), vol. 537.
149
F. Leisch, Flexmix: A general framework for finite mixture models and latent class regression In R (2004).
(0)eLetters
eLetters is a forum for ongoing peer review. eLetters are not edited, proofread, or indexed, but they are screened. eLetters should provide substantive and scholarly commentary on the article. Neither embedded figures nor equations with special characters can be submitted, and we discourage the use of figures and equations within eLetters in general. If a figure or equation is essential, please include within the text of the eLetter a link to the figure, equation, or full text with special characters at a public repository with versioning, such as Zenodo. Please read our Terms of Service before submitting an eLetter.
Log In to Submit a ResponseNo eLetters have been published for this article yet.
Information & Authors
Information
Published In
Science Advances
Volume 6 | Issue 22
May 2020
May 2020
Copyright
Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
Submission history
Received: 17 September 2019
Accepted: 23 March 2020
Acknowledgments
We thank F. Aron, M. Burri, A. Mittnik, K. Nägele, R. Radzeviciute, E. Späth, R. Stahl, A. Wissgott, and G. Brandt for technical support in the DNA analyses. We thank the teams at MPI-SHH-Archaeogenetics and PACEA, University of Bordeaux for continued support and discussion, especially S. Clayton, J. Y. Fellows, T. Lamnidis, L. Papac, and V. Villalba-Mouco. We thank S. Nakagome for a script for allele frequency estimation. We thank H. Duday, D. Hoffmann, B. Gehlen, and A. Powell for critical discussions. We thank H.-C. Strien for the discussion of SMH radiocarbon dates. We thank the project Counter Culture: Investigating Neolithic social diversity for analyses on isotopic data for OBN and analyses carried out at the BioArCh facility, Department of Archaeology, University of York. Funding: M.R. was supported by the Fyssen Foundation (post-doctoral stipend, 2017–2018). C.J. was supported by the New Faculty Startup Fund from Seoul National University. The study was funded by the Max Planck Society, and the French (ANR) and German (DFG) Research Foundations, INTERACT project, ANR-17-FRAL-0010, DFG-HA-5407/4-1, 2018-2021 (M.-F.D., W.H., and M.R.), and the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program under grant agreement no. 771234 – PALEoRIDER (W.H. and A.B.R.). We thank the DFG, grant no. RO 4148/1 (E.R.), for funding SMH radiocarbon dates. Author contributions: M.R., M.-F.D., and W.H. conceived the study. K.W.A., D.B., S.F., E.G., D.G., L.L., P.L., H.M., H.R., E.R., S.R., C.S., L.S., and J.W. provided archaeological material and advised on the archaeological background and interpretation. M.R., İ.K., and M.-H.P. performed the laboratory work. M.R. performed the data analyses, with C.J., S.S., A.B.R., J.K., and W.H. providing guidance and methodologies. M.R., D.G., D.B., M.-F.D., and W.H. wrote the manuscript with input from all coauthors. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Individual bam files (aligned human reads) are available through the European Nucleotide Archive under accession number PRJEB36208.
Authors
Funding Information
H2020 European Research Council: 771234
Deutsche Forschungsgemeinschaft: DFG-HA-5407/4-1
Deutsche Forschungsgemeinschaft: RO 4148/1
Agence Nationale de la Recherche: ANR-17-FRAL-0010
Seoul National University: New Faculty Startup Fund
Max Planck Society
Metrics & Citations
Metrics
Article Usage
- 134 citation in Crossref
- 119 citation in Web of Science
Altmetrics
Citations
Cite as
- Maïté Rivollat et al.
Export citation
Select the format you want to export the citation of this publication.
Cited by
- Social dynamics in past societies, Encyclopedia of Measurement in Social Sciences, (600-612), (2027).https://doi.org/10.1016/B978-0-443-26629-4.00189-1
- Local ancestry inference identifies robust evidence of selection in Neolithic Europe, Molecular Biology and Evolution, 43, 9, (2026).https://doi.org/10.1093/molbev/msag236
- The Oral Microbiome of King Richard III of England , American Journal of Biological Anthropology, 191, 1, (2026).https://doi.org/10.1002/ajpa.70350
- Ancient DNA Reveals Limited Hunter‐Gatherer Ancestry in the Cerny Culture (Middle Neolithic, France) Despite Hunting‐Associated Practices , American Journal of Biological Anthropology, 191, 1, (2026).https://doi.org/10.1002/ajpa.70351
- Diachronic reuse of Neolithic burial monuments by Bronze Age newcomers in Southwestern Britain, Scientific Reports, 16, 1, (2026).https://doi.org/10.1038/s41598-026-66094-z
- The Evolution of Molecular Clocks: Concepts, Calibrations, and Challenges, Biology, 15, 15, (1329), (2026).https://doi.org/10.3390/biology15151329
- Genes, Identity, and Material Culture: In Search of a Unifying Theory, Culture Evolving, (115-143), (2026).https://doi.org/10.1007/978-3-032-16124-6_5
- Ancient DNA reveals pervasive directional selection across West Eurasia, Nature, 654, 8118, (419-428), (2026).https://doi.org/10.1038/s41586-026-10358-1
- Population discontinuity in the Paris Basin linked to evidence of the Neolithic decline, Nature Ecology & Evolution, 10, 4, (677-688), (2026).https://doi.org/10.1038/s41559-026-03027-z
- Origins and interactions of Neolithic populations from Liguria (NW Italy): a craniofacial morphometric study, Archaeological and Anthropological Sciences, 18, 4, (2026).https://doi.org/10.1007/s12520-026-02440-6
- See more
Loading...
View Options
View options
PDF format
Download this article as a PDF file
Download PDFLog in to view the full text
AAAS login provides access to Science for AAAS Members, and access to other journals in the Science family to users who have purchased individual subscriptions.
More options
Register for free to read this article
As a service to the community, this article is available for free. Login or register for free to read this article.
