Archaeogenomic analysis of the first steps of Neolithization in Anatolia and the Aegean
Abstract
The Neolithic transition in west Eurasia occurred in two main steps: the gradual development of sedentism and plant cultivation in the Near East and the subsequent spread of Neolithic cultures into the Aegean and across Europe after 7000 cal BCE. Here, we use published ancient genomes to investigate gene flow events in west Eurasia during the Neolithic transition. We confirm that the Early Neolithic central Anatolians in the ninth millennium BCE were probably descendants of local hunter–gatherers, rather than immigrants from the Levant or Iran. We further study the emergence of post-7000 cal BCE north Aegean Neolithic communities. Although Aegean farmers have frequently been assumed to be colonists originating from either central Anatolia or from the Levant, our findings raise alternative possibilities: north Aegean Neolithic populations may have been the product of multiple westward migrations, including south Anatolian emigrants, or they may have been descendants of local Aegean Mesolithic groups who adopted farming. These scenarios are consistent with the diversity of material cultures among Aegean Neolithic communities and the inheritance of local forager know-how. The demographic and cultural dynamics behind the earliest spread of Neolithic culture in the Aegean could therefore be distinct from the subsequent Neolithization of mainland Europe.
1. Introduction
The primary zone of Neolithization in western Eurasia encompassed the Levant, Taurus-Zagros ranges of Mesopotamia, central Anatolia and Cyprus [1–4]. The earliest evidence for sedentary life and food storage in this region goes back to the Natufians (ca 12 500–10 800 cal BCE) [5,6]. Sedentary communities were established across this zone during the first phase of the Pre-Pottery Neolithic (PPN, or Aceramic Neolithic, ca 10 000–8500 cal BCE), and the first indications of plant cultivation appeared [7–9]. Between ca 8500 and 7000 cal BCE, community sizes increased, architectural elaboration intensified and a subsistence economy based on agriculture gradually became the norm [10–14]. Meanwhile, portable artefacts such as figurines and stamps evolved into staples of sedentary life, and pottery production became widespread ca 7000 cal BCE [10,11]. The elements of the subsequent Pottery Neolithic culture (PN, ca 7000–5500 cal BCE), including integrated cultivation practices of domestic plants and animals, the architectural practices of sedentary life, together with portable artefacts, have been collectively described as the Near Eastern ‘Neolithic Package’ [15–18].
During the same period, there were no signs of a Neolithization in west Anatolia and the Aegean. Only after ca 7000 cal BCE did elements of the ‘Neolithic Package’ appear in these regions, eventually spreading towards Europe [19–21]. Some archaeologists suggest that the emergence of the Neolithic elements in the Aegean and in Europe without a preceding PPN development period indicates the role of demic processes (i.e. migrations from the Neolithic primary zones through land and sea routes), frequently described as a leap-frog model where migrants form enclaves in new territory [15,16,22–27]. Others, in contrast, favour a role for interaction between local foragers and primary zone Neolithic populations, including the adoption of Neolithic elements by locals and acculturation [16,18,28,29].
Recent archaeogenomic data have shown that the Neolithization of central, western and northern Europe involved migration from a single eastern source, frequently termed ‘Anatolian farmer’ [30–33], while in other regions, such as in the Baltic [33,34] and in southern Greece [33], acculturation may have played role. In most of Europe, there is limited genetic evidence for early admixture between farmers and local European Mesolithic (WHG) communities in the sixth millennium BCE, such that early European farmers studied to date (with few exceptions [35]) carry ancestry similar to farmers from northwest Anatolian Barcın [31–33] (electronic supplementary material, figure S1). In subsequent millenia, however, WHG-like ancestry appears in Middle and Late Neolithic European populations [36–38]. These observations support a leap-frog model of Neolithic spread in Europe [28]: farmers only occupied enclaves in the new territories while Mesolithic groups persisted in the same regions [39–42].
The processes behind the earliest steps of Neolithization and the Neolithic spread in the Aegean are less understood. For instance, whether Aegean Neolithic populations were recent colonists originating from areas of the primary zone of Neolithization (e.g. [27]), descendants of indigenous foragers (e.g. [20]) or admixed groups (e.g. [15,16]) is still contentious. Additionally, whether Aegeans' demographic or cultural relationships were stronger with central Anatolians [16] or with the Levantine seafaring populations [27] remains unclear. We re-analyse published ancient human genomes to answer these questions and to dissect the demographic dynamics behind the Neolithic transition in Anatolia and the Aegean.
2. Material and methods
(a) Compiling and mapping genomic data
We obtained DNA sequencing data of 99 published ancient individuals (electronic supplementary material, table S1), generated using whole-genome shotgun sequencing and/or sequencing of libraries enriched by hybridization capture [30–32,35,37,38,42–46]. We mapped sequencing reads to the human reference genome (hs37d5) using the Burrows–Wheeler Aligner (BWA, v. 0.7.12) [47], with the parameters ‘-l 16500, -n 0.01, -o 2’. We filtered PCR duplicates using FilterUniqSAMCons.py [48]. We filtered reads shorter than 35 base pairs, with greater than 10% mismatches to the reference, and less than 30 mapping quality per read.
(b) Preparation of population genetics analysis datasets
We restricted our analysis to known present-day DNA variants to minimize false positives. We used two different modern reference panels, calling genotypes of ancient individuals for SNPs overlapping with (i) the Human Origins genotype dataset [42,49] and (ii) the 1000 Genomes whole-genome sequence data [50] using SAMtools mpileup (v. 1.3) [51]. For (i) we obtained a curated version of the Human Origins panel of 594 924 autosomal SNP genotype calls for 2730 present-day individuals from Lazaridis et al. [42]. We determined the SNPs of the ancient samples overlapping with this dataset. We encoded transitions as missing to avoid confounding with cytosine deamination in ancient DNA. To prepare (ii), we obtained the BAM and VCF files for the African Yoruba individuals from 1000 Genomes Project phase 3 from McVean et al. [50]. Using vcftools [52], we extracted a total of 1 938 919 transversion SNPs with minor allele frequencies of greater than or equal to 10% in the Yoruba population to avoid false positive calls [36,40]. We determined the positions in the ancient samples overlapping with this dataset. We merged ancient genotypes with these two datasets using PLINK [53] requiring base quality greater than or equal to 30 per overlapping position. We haploidized each full dataset by randomly selecting one allele per position. The Human Origins-merged dataset, which has higher number of present-day populations, was used for principal component analysis (PCA) and for calculating f3-statistics. The 1000 Genomes-merged dataset, with a higher number of SNPs, was used for D-statistics, where we require high statistical power.
(c) Principal component analysis
We performed PCA by calculating principal components using west Eurasian populations from the Human Origins dataset using the smartpca program of EIGENSOFT [54], with the ‘numoutlieriter:0’ parameter. We projected ancient genomes onto the reference space using the ‘lsqproject:YES’ option and plotted the results using R (v. 3.3.0).
(d) D- and f3-statistics
We computed D-statistics using the qpDstat program of the ADMIXTOOLS package [49]. We assessed statistical significance by calculating standard errors using a block jacknife of 0.5 Mbp. We used the Yoruba population as outgroup for the D-statistics [32]. We computed f3-statistics (i.e. genetic affinity between pairs of populations based on an estimate of shared drift between them as their divergence from an outgroup population) using the qp3Pop program of the ADMIXTOOLS package [49]. The Human Origins dataset's African Mbuti population was used as outgroup for calculating f3-statistics [42]. We performed multiple testing correction using Benjamini–Yekutieli method for all 207 D-statistics results and reported adjusted p-values together with Z scores per each test [55]. For the pairwise f3-statistics, as genetic distance measure between a pair of populations, X and Y, we used: 1-f3(Mbuti;X,Y) [30]. These pairwise distances were summarized with the multidimensional scaling (MDS) method using the cmdscale function of R. We evaluated the goodness of fit for MDS using ‘GOF’ component obtained from cmdscale function.
(e) Heterozygosity estimates
We calculated heterozygosity as a measure of genetic diversity in a population, using genome sequence data of (i) Bon002 (from Boncuklu, central Anatolia, pre-7000 cal BCE) [30], Tep003 (Tepecik-Çiftlik, central Anatolia, post-7000 cal BCE) [30], Bar8 (Barcın, north Aegean, post-7000 cal BCE) [31] and Rev5 (Revenia, north Aegean, post-7000 cal BCE) [31]. We calculated genome coverage per sample using GenomeCoverageBed [56]. We downsampled the genome sequences of Bon002 and Bar8 to similar levels as the other two samples using SAMtools (v. 1.3) [51]. We calculated heterozygosity per sample using ANGSD [57] as ‘angsd -GL 1 -doGlf 2 -doMajorMinor 1 -sites ReferenceSNP.pos -bam bamlist -doSaf 1 -anc referencegenome.fasta’. To minimize false positives, we only considered transversions overlapping with of Yoruba individuals from 1000 Genomes Project phase 3 from McVean et al. [50].
(f) Modelling of admixture
We used the qpWave/qpAdm framework [38,58] in the ADMIXTOOLS package [49] to model populations as mixtures of two or more sources. The following worldwide set of ancient and present-day outgroups, which most probably did not experience any post-split gene flow from Anatolian/Aegean populations, was used: Mbuti, Yoruba, Ust Ishim, El Miron, Goyet Q116, Villabruna, Kostenki14, Vestonice16, Papuan, Onge, Karitiana, Mixe, Chipewyan, Oroqen, Koryak, Dai, Japanese. Adding East European hunter–gatherers (EHG) as a close outgroup to increase the resolution did not change the results.
(g) Serial coalescent simulations
We performed serial coalescent simulations using fastsimcoal [59] under four various demographic models involving Neolithic central Anatolians, Aegeans, Iranians and WHG (not including the Levantine populations, for whom we lack whole-genome data). The simulations were designed to mimick the data with respect to tree topology, divergence times and sample sizes. We then performed D-statistics on the simulated DNA and compared these with the observed data to gain understanding into the plausibility of different models. Specifically, we generated data to represent Iranian Neolithics (10 000 BP), WHGs (Loschbour: 7200 BP), central Anatolian Neolithics (Tepecik-Çiftlik: 8500 BP; Boncuklu: 10 000 BP), the Aegean Neolithics (Revenia: 8300 BP) and present-day sub-Saharan Africans (Yoruba-YRI). We launched 100 runs for each model defined in the parameter file (input.par) for testing different population histories. For all models, we sampled 30 Mb DNA sequences for: five present-day Yoruba, two Iranian Neolithics, two WHGs, four central Anatolian Neolithics (two Tepecik-Çiftlik, two Boncuklu) and two Aegean Neolithics (Revenia). We assumed a mutation rate of 1.00 × 10−9 bp yr−1, and a recombination rate of 1.00 × 10−8 bp yr−1, and assumed 25 years per generation, again following [45]. We set the effective population size (Ne) of these populations and times of divergence between Anatolian Neolithic, WHGs and Iranian Neolithic populations based on [45]. We converted all outputs (arp file) to plink format and computed D-statistics with topology of D(YRI, Test, central Anatolian N, Aegean N) to test the relationships among populations via ADMIXTOOLS [49]. Note that the tree topology involving the Anatolian/Aegean populations, Iran, WHG and the Africans, were based on the phylogenetic analysis from Broushaki et al. [45]. The Anatolian/Aegean populations were assumed to diverge simultaneously from the same source (star shaped).
3. Results
(a) Early Holocene gene pools of west Eurasia and the Anatolian/Aegean gene pool
We compiled published genome sequence data of 99 ancient individuals (sample ages: ca 11 840–4360 cal BCE) (electronic supplementary material, figure S2, table S1). Both a PCA using present-day and ancient populations (electronic supplementary material, figure S1) and an MDS analysis using only ancient genomes (figures 1 and 2; electronic supplementary material, table S2) revealed the presence of four distinct gene pools in Early Holocene west Eurasia: (a) a ‘Caucasia/Iran gene pool’, (b) a ‘Levant gene pool’, (c) a ‘European pre-Neolithic gene pool’ and (d) an ‘Anatolian/Aegean gene pool’. To objectively measure clustering in gene pools (a)–(c), we used D-statistics of the form D(Yoruba, p1; p2, p3) where ‘p’ refers to the Caucasia/Iran, the Levant or European pre-Neolithic gene pools, correcting for multiple testing. In 80% comparisons (p < 0.05; Z ≥ 3), populations belonging to the same gene pool shared more alleles with each other compared to external populations (figure 1a; electronic supplementary material, figure S3a, table S3). The only exceptions were comparisons involving a single pre-Neolithic individual from Iran for which we had relatively few SNPs and low statistical power.
We then investigated the relationships among ancient Anatolians and other west Eurasian gene pools, using the oldest Anatolian population yet sequenced: Boncuklu from central Anatolia (sample ages: ca 8300–7952 cal BCE), an Aceramic Neolithic population previously predicted to be the descendants of local Epipalaeolithic groups [30,60]. We computed D-statistics of the form D(Yoruba, p1; p2, Boncuklu), where ‘p1’ and ‘p2’ refer to populations belonging to different gene pools: Caucasia/Iran, the Levant or the European pre-Neolithic. In 56% of the comparisons (p < 0.05; Z ≥ 2.8), all three regional gene pools showed higher affinity to Boncuklu than to each other (figure 1a,b; electronic supplementary material, table S4). Using the qpWave/qpAdm algorithm [38,58], we further modelled the Boncuklu population as a mixture of CHG (59.1%), the Levant (31.4%) and WHG (9.5%) (electronic supplementary material, table S5).
We next included three post-7000 cal BCE Neolithic populations from Anatolia and Aegean in the analyses: Tepecik-Çiftlik in central Anatolia [30,61], Barcın in northwest Anatolia [31,37,62] and ‘Revenia’ in Pieria of northeast Greece [31]. We computed D-statistics of the form D(Yoruba, p1; p3, p2) where ‘p1’ and ‘p2’ are Anatolian/Aegean populations, and ‘p3’ is an external population (Caucasia/Iran, the Levant or European pre-Neolithic). In 94% of the comparisons (p < 0.05; Z ≥ 2.8) all Anatolian/Aegean populations were genetically closer to each other than to any other gene pool (electronic supplementary material, figure S3b, table S5).
Given archaeological indication that Aegean Neolithic was influenced by east Mediterranean sources [27], we further studied the genetic affinities of Aegean Neolithic people to central Anatolian Neolithics and to the Levantines. Calculating D-statistics of the form D(Yoruba, northAegean; Levant, centralAnatolia) revealed that the post-7000 cal BCE Neolithic north Aegean individuals (Barcın and Revenia) consistently share more alleles with central Anatolians compared to south Levantines, where 50% of the comparisons were significant (p < 0.05; Z ≥ 2.8) (electronic supplementary material, figure S3b, table S5).
(b) Notable genetic diversity in the Aegean
To assess demographic events in the Near East during the Neolithic transition, we studied signatures of regional admixture using diachronic populations from the same region (figure 1c). In 83% of the comparisons, pre-7000 cal BCE Neolithic populations of the Levant and of Iran were genetically closer to all post-7000 cal BCE Anatolian/Aegean populations (Tepecik, Barcın, Revenia) compared with the pre-7000 cal BCE Anatolian Boncuklu (p < 0.05; Z ≥ 3) (electronic supplementary material, figure S3c, table S6). Considering the radiocarbon dates of the investigated individuals, this is consistent with gene flow from both the Levant and from Iran into Anatolia, within a period ranging from the PPN to the PN (figure 1c, arrows ‘e’ and ‘f’). These results are also compatible with a regional increase in the levels of admixture during the Neolithic [30,32], although alternative explanations to gene flow remain plausible, such as population structure confounding the analysis results [63].
Next, to gain understanding into Aegean Neolithization, we studied the population genetic characteristics of the PN Aegean groups relative to central Anatolian groups. We first compared heterozygosity estimates among these populations. If the Aegeans were recent colonists from a single origin, due to a founder effect, one might expect lower heterozygosity in the Aegean than in central Anatolia. By contrast, Barcın and Revenia individuals had higher heterozygosity levels (mean 0.25 and 0.26, respectively) than those of Boncuklu and Tepecik (0.22 and 0.19, respectively) (electronic supplementary material, table S7).
Second, we calculated D-statistics focused on the Aegeans, which suggested higher admixture in this region than in central Anatolia:
(i) D(Yoruba, Natufian; centralAnatolia, northAegean) revealed that pre-Neolithic population of the Levant had stronger genetic affinity to the two north Aegean Neolithic populations (Barcın and Revenia, post-7000 cal BCE) than to the two central Anatolian Neolithic groups (Boncuklu and Tepecik, pre- and post-7000 cal BCE) (p < 0.05; Z ≥ 3) (figure 2a, arrow ‘a’; electronic supplementary material, figure S4a, table S8). Given the above-proposed gene flow event from the Levant into Anatolia during the Neolithic, this result might imply additional genetic interactions between Natufian-related populations and the ancestors of north Aegean populations that bypassed central Anatolia.
(ii) D(Yoruba, Caucasia/Iran; centralAnatolia, northAegean) revealed that in 50% of the comparisons CHGs and Neolithic Iran individuals shared more alleles with the two north Aegean PN populations than with the two central Anatolians (p < 0.05; Z ≥ 2.8) (figure 2a, arrows ‘b’ and ‘c’; electronic supplementary material, figure S4a, table S9).
(iii) Likewise, WHG individuals showed higher affinity to the two north Aegean PN populations than PN central Anatolian group groups (p < 0.05; Z ≥ 3) (figure 2a, arrow ‘d’; electronic supplementary material, figure S4a, table S10).
(iv) Natufians, WHGs and Iranian PPN individuals were consistently more similar to the Revenia individual than to those in Barcın (electronic supplementary material, figure S4b, table S11).
(v) Both the Boncuklu (PPN) and the Tepecik (PN) groups of central Anatolia had stronger affinity to the north Aegean PN populations, Barcın and Revenia, than to each other (figure 2b; electronic supplementary material, table S12). Likewise, all Anatolian groups (Boncuklu, Tepecik-Çiftlik and Barcın) were genetically closer to Revenia than they were to each other (figure 2b; electronic supplementary material, table S13).
(vi) All European early farmer populations examined were genetically closer to Revenia than to each other (electronic supplementary material, figure S4c, table S14)
Observations (ii), (iii) and in particular (v) are intriguing. We asked whether these could be consistent with a number of demographic scenarios, assuming a phylogenetic topology that included Iran, WHG and Aegean/Anatolian populations, estimated by [45]. We considered the following scenarios: (a) separate extreme bottlenecks in the ancestors of the two central Anatolian populations (possibly causing differentiation between the central Anatolian populations from each other, and from all other groups, (b) independent gene flow events from external sources (WHG and Iran) into the two central Anatolian groups (possibly causing differentiation between the two), (c) independent gene flow from WHG and Iran into the Aegean, and (d) independent gene flow from WHG, Iran and the two central Anatolian lineages into the Aegean. We performed serial coalescent simulations using realistic settings and compared the results with the observed D-statistics. We could only replicate the observed results under scenario (d) that describes rampant admixture in the Aegean (electronic supplementary material, figure S5).
4. Discussion
The analyses presented here highlight two points regarding the process of Neolithization. First, the observation that the two central Anatolian populations cluster together to the exclusion of Neolithic populations of south Levant or of Iran restates the conclusion that farming in central Anatolia in the PPN was established by local groups instead of immigrants, which is consistent with the described cultural continuity between central Anatolian Epipalaeolithic and Aceramic communities [9,64]. This reiterates the earlier conclusion [32] that the early Neolithization in the primary zone was largely a process of cultural interaction instead of gene flow.
The second point relates to whether Aegean Neolithization (post-7000 cal BCE) involved similar acculturation processes, or was driven by migration similar to Neolithization in mainland Europe—a long-standing debate in archaeology [16,20,22,27,28]. Here, we discuss the two scenarios based on the genetic analysis.
(a) Model 1: migration from Anatolia to the Aegean
A recent study reported that by the seventh millennium BCE the eastward border of the WHG gene pool extended to the Iron Gates (on the border between Romania and Serbia) [33]. Plausibly, during the Early Holocene, the WHG population could also have been present along the Aegean coastline, such that the border between central Anatolian and WHG gene pools ran along west Anatolia. If so, the Aegean Neolithization must have involved replacement of a local, WHG-related Mesolithic population by incoming easterners.
If migration occurred, where did it originate? Because Revenia and Barcın cluster with PPN and PN central Anatolian Neolithic groups to the exclusion of the south Levant (figure 1c; electronic supplementary material, figure S3c), the latter is unlikely to be the source, leaving central Anatolia or south Anatolia (north Levant) as potential origins.
Notably, the north Aegeans (Revenia and Barcın) show higher diversity than the central Anatolians. We had earlier shown that the highest-quality Barcın genome carries a smaller proportion of short runs of homozygosity than the highest quality Boncuklu genome [30], which also supports the notion that the ancestral effective population size of the Aegeans was larger than those of central Anatolians. Moreover, we find that the north Aegeans share more alleles with eastern, western and southern gene pools, as estimated using the D-statistic (figure 2). Although the D-statistic can be sensitive to technical biases, our result is unlikely to be a technical artefact because (a) the north Aegean data were derived from two independent studies [31,37], (b) the Barcın data were produced using two different techniques, whole-genome shotgun sequencing and SNP capture, and (c) both Barcın and Revenia display the same population genetic patterns, suggesting that the admixture signals in the Aegean individuals are reproducible. In addition, although unknown population structure can complicate interpretation of the D-statistic [63], we note that the admixture estimates are consistent with the estimated higher genetic diversity in the Aegean.
If the Revenia and Barcın individuals studied here were descendants of Anatolian Neolithic immigrants, they must have been recent settlers, as all samples analysed here date to early stages of the Aegean Neolithic (Revenia: 6438–6264 and Barcın: 6500–6200 cal BCE). Furthermore, if the migration was directly of central Anatolian origin (represented by Boncuklu and Tepecik-Çiftlik), the putative migrants must have admixed with populations carrying alleles of distinct gene pools (the Levant, Caucasus/Iran and WHG) within a few centuries, in order to explain our observations above (figure 2a).
Alternatively, the migration event could have originated from the Anatolian south coast or north Levant [27] (currently no genome data are available from these groups). This region could have hosted a hypothetical central Anatolian-related population exposed to admixture from CHG-, Iran- and the Levant-related gene pools in earlier millenia. A south Anatolian population could have been in contact with different central Anatolian populations from the Konya Plain (Boncuklu) and Cappadocia (Tepecik-Çiftlik), explaining the affinity of both Boncuklu and Tepecik-Çiftlik to Barcın. A seafaring population could also be in genetic contact with putative WHG-related populations of the Aegean. This hypothetical population could have initiated the Cyprus Neolithic in the eleventh millennium BCE and later Aegean Neolithic communities in the seventh millennium BCE [27].
One surprising observation here is the apparent absence of WHG-like ancestry in Late Neolithic/Chalcolithic Aegean genomes: admixture analysis results from two individuals from northwest Anatolia (Kumtepe, approx. 5000 BCE) [65] and four individuals from south Greece (Franchti Cave and Diros, approxi. 4000 BCE) [33] all lack noticeable WHG-like ancestry components [30,31,33]. This contrasts with WHG admixture emerging in European farmer populations in the Middle and Late Neolithic [36,38], and perhaps earlier in the Balkans [33], indicating the persistence of Mesolithic populations in Europe after Neolithic migrations. Therefore, if the Mesolithic populations of the Aegean coast had indeed been WHG-related, they must have been fully replaced by the eastern migrant farmers.
(b) Model 2: adoption of Neolithic elements by local foragers
Alternatively, the Aegean coast Mesolithic populations may have been part of the Anatolian-related gene pool that occupied the Aegean seaboard during the Early Holocene. Under this scenario, the north Aegean PN populations would be at least partial descendants of local hunter–gatherers who adopted Neolithic lifestyle post-7000 cal BCE, triggered by contacts with central Anatolian and the Levantine populations. The following events would be conceivable. (a) during the Last Glacial Maximum (LGM), the Aegean evolved into a refuge hosting a significant human population, which is in line with climatic modelling [66–69]; estimates of human population density during the Marine Isotope Stage 2 in west (but not central) Anatolia reach one of their highest levels in Europe [70,71]. The existence of an Aegean human population going back to the LGM is also consistent with mitochondrial haplogroup-based analyses [72], and that Anatolian-like mitochondrial haplogroups are found also in Mesolithic Balkan and Aegean populations [31,33]. (b) Following the LGM, Aegean emigrants dispersed into central Anatolia and established populations that eventually gave rise to the local Epipalaeolithic and later Neolithic communities, in line with the earliest direct evidence for human presence in central Anatolia ca 14 000 cal BCE [60]. This hypothetical out-of-the-Aegean event coincides with the post-LGM Near East-related migration signatures in European Mesolithic genomes [73]. (c) Between the LGM and post-7000 cal BCE Neolithization, WHG, Natufian and Caucasus/Iran-related groups admixed with north Aegeans, differentiating the latter from their central Anatolian relatives and leading to our observations in figure 2a. (d) Post-7000 cal BCE, there occurred additional, albeit limited central Anatolian gene flow back into the Aegean, giving rise to our observation in figure 2b.
(c) The archaeological evidence
Both the migration and acculturation models for Aegean Neolithization enjoy support from material culture investigations, but the overall evidence points to a complex process where Aegean societies were culturally influenced by diverse sources, including the central Anatolian Neolithic, the Levant Neolithic and possibly local Mesolithic traditions. In contrast to the relative homogeneity of European Neolithic cultures, such as the LBK and Cardial, the Aegean Neolithic is noted for its diversity [64]. Variation in Neolithic Package elements and primary zone traditions is notable across Aegean sites, among regions (e.g. east and west of Marmara), even between closely neighbouring villages [16,17,20,64,74–79]. This diversity includes, for example, obsidian, with Greek Aegean (Melos) [80,81] or mainland Anatolian (Cappadocian) [82] sources being preferred in some settlements, and yet other settlements showing no evidence of obsidian use [64]. Cultural trait diversity involves architecture, tool types, ceramics and symbolic elements (such as figurines and intramural burial), which may show partial similarities to either central Anatolia or to the Levant, or may be unique [16,64,75]. For instance, intramural burial, a common feature among primary zone sites, is also widespread in east Marmara early Neolithic villages (including Barcın), but totally absent in settlements only 200 km west [16]. Mesolithic-like lithic industries and the prominence of seafood in some settlements further imply the continuing presence of Aegean Mesolithic traditions into the Neolithic [16,20,27,83,84]. Indeed, lively seafaring activity was prevalent in the Mediterranean and the Aegean already by the eleventh millennium BCE [85,86], as evidence from Cyprus, Crete, Franchti, Cyclops Cave, Ouirakos and other Aegean island and coastal mainland Mesolithic sites demonstrates [20,27,83,86–94].
Instead of a single-sourced colonization process, the Aegean Neolithization may thus have flourished upon already existing coastal and interior interaction networks connecting Aegean foragers with the Levantine and central Anatolian PPN populations, and involved multiple cultural interaction events from its early steps onward [16,20,64,74]. This wide diversity of cultural sources and the potential role of local populations in Neolithic development may set apart Aegean Neolithization from that in mainland Europe. While Mesolithic Aegean genetic data are awaited to fully resolve this issue, researchers should be aware of the possibility that the initial emergence of the Neolithic elements in the Aegean, at least in the north Aegean, involved cultural and demographic dynamics different than those in European Neolithization.
Data accessibility
This article has no additional data.
Authors' contributions
M.S., I.T., G.M.K., A.G. and Ç.A. designed and supervised the study; G.M.K., D.K. and A.P.S. analysed data with contributions from H.M.D., R.Y., E.A., P.F.; Ç.A., D.B., A.M.B. and C.C.B. provided archaeological and ecological interpretations; M.S., I.T., G.M.K. and Ç.A. wrote the manuscript with input from all authors.
Competing interests
We declare we have no competing interests.
Funding
This work was supported by TÜBİTAK (grant no. 114Z927 to M.S.), TÜBA GEBİP (to M.S.), METU BAP (to M.S.), Knut and Alice Wallenberg Foundation: ‘1000 ancient genomes' (to A.G.). E.A. was supported by the Institutional Development Program of the University of Ostrava and by EU structural funding Operational Programme Research and Development for Innovation, project No. CZ.1.05/2.1.00/19.0388.
Acknowledgements
We thank T. Günther, A. Munters, S. C. Acan, A. Birand, Y. S. Erdal, A. Fairbairn, A. Omrak, M. Jakobsson, M. Krzewinska, F. Özer, R. Özbal, L. Excoffier and J. Stora for technical support in analysis, suggestions and/or comments on the manuscript, and two anonymous reviewers for helpful suggestions.
References
- 1
Kelly RL . 1995 The foraging spectrum: diversity in hunter-gatherer lifeways. Washington, DC: Smithsonian Institution Press. Google Scholar - 2
Zeder MA . 2009 Evolutionary biology and the emergence of agriculture: the value of co-opted models of evolution in the study of culture change. In Macroevolution in human prehistory (ed.Prentiss AM ), pp. 157–210. New York, NY: Springer New York. Crossref, Google Scholar - 3
Aiello LC . 2011 The origins of agriculture: new data, new ideas. Curr. Anthropol. 52, S161–S162. (doi:10.1086/660154) Crossref, Web of Science, Google Scholar - 4
Özdoğan M . 1997 The beginning of neolithic economies in southeastern Europe: an anatolian perspective. J. Eur. Archaeol. 5, 1–33. (doi:10.1179/096576697800660267) Crossref, Google Scholar - 5
Bar-Yosef O . 2001 From sedentary foragers to village hierarchies: the emergence of social institutions. Proc. Br. Acad. 101, 1–38. Google Scholar - 6
Kuijt I, Finlayson B . 2009 Evidence for food storage and predomestication granaries 11 000 years ago in the Jordan Valley. Proc. Natl Acad. Sci. USA 106, 10 966–10 970. (doi:10.1073/pnas.0812764106) Crossref, Web of Science, Google Scholar - 7
Zeder MA . 2012 The broad spectrum revolution at 40: resource diversity, intensification, and an alternative to optimal foraging explanations. J. Anthropol. Archaeol. 31, 241–264. (doi:10.1016/j.jaa.2012.03.003) Crossref, Web of Science, Google Scholar - 8
Weiss E, Zohary D . 2011 The Neolithic Southwest Asian founder crops. Curr. Anthropol. 52, S237–S254. (doi:10.1086/658367) Crossref, Web of Science, Google Scholar - 9
Baird D. 2012 The late epipaleolithic, Neolithic, and chalcolithic of the anatolian plateau, 13 000–4000 BC. In A companion to the archaeology of the ancient near east (ed.Potts DT ), pp. 431–465. Oxford, UK: Wiley-Blackwell. Crossref, Google Scholar - 10
Kuijt I, Goring-Morris N . 2002 Foraging, farming, and social complexity in the pre-Pottery neolithic of the southern Levant: a review and synthesis. J. World Prehistory 16, 361–440. (doi:10.1023/A:1022973114090) Crossref, Web of Science, Google Scholar - 11
Akkermans PMMG, Schwartz GM . 2003 The archaeology of Syria: from complex hunter-gatherers to early urban societies (c. 16 000–300 BC). Cambridge, UK: Cambridge University Press. Google Scholar - 12
Flannery KV . 2002 The origins of the village revisited: from nuclear to extended households. Am. Antiq. 67, 417–433. (doi:10.2307/1593820) Crossref, Web of Science, Google Scholar - 13
Asouti E, Fuller DQ . 2013 A contextual approach to the emergence of agriculture in southwest Asia. Curr. Anthropol. 54, 299–345. (doi:10.1086/670679) Crossref, Web of Science, Google Scholar - 14
Atakuman Ç . 2014 Architectural discourse and social transformation during the early Neolithic of southeast Anatolia. J. World Prehistory 27, 1–42. (doi:10.1007/s10963-014-9070-4) Crossref, Web of Science, Google Scholar - 15
Özdoğan M . 2011 Archaeological evidence on the westward expansion of farming communities from eastern Anatolia to the Aegean and the Balkans. Curr. Anthropol. 52, S415–S430. (doi:10.1086/658895) Crossref, Web of Science, Google Scholar - 16
Özdoğan M . 2014 A new look at the introduction of the Neolithic way of life in Southeastern Europe. Changing paradigms of the expansion of the Neolithic way of life. Doc. Praehist. 41, 33–49. (doi:10.4312/dp.41.2) Crossref, Google Scholar - 17
Perles C . 2005 From the near east to greece: let's reverse the focus—cultural elements that did not transfer. In How did farming reach Europe: anatolian-european relations from the second half of the 7th through the first half of the 6th millennium cal BC (ed.Lichter C ), pp. 275–290. İstanbul, Turkey: Byzas 2. Google Scholar - 18
Zvelebil M . 1989 On the transition to farming in Europe, or what was spreading with the Neolithic: a reply to Ammerman (1989). Antiquity 63, 379–383. (doi:10.1017/S0003598X00076110) Crossref, Web of Science, Google Scholar - 19
Zvelebil M . 2009 Mesolithic societies and the transition to farming: problems of time, scale and organisation. In Hunters in transition: mesolithic societies of temperate eurasia and their transition to farming (edsZvelebil M, Dennell R, Domanka L ), pp. 167–188. Cambridge, UK: Cambridge University Press. Google Scholar - 20
Reingruber A . 2011 Early Neolithic settlement patterns and exchange networks in the Aegean. Doc. Praehist. 38, 291. (doi:10.4312/dp.38.23) Crossref, Google Scholar - 21
Brami M, Zanotti A . 2015 Modelling the initial expansion of the Neolithic out of Anatolia. Doc. Praehist. 42, 103. (doi:10.4312/dp.42.6) Google Scholar - 22
Çilingiroğlu Ç, Çakırlar C . 2013 Towards configuring the neolithisation of Aegean Turkey. Doc. Praehist. 40, 21–29. (doi:10.4312/dp.40.3) Crossref, Google Scholar - 23
Ammerman AJ, Cavalli-Sforza LL, Luigi L . 1984 The neolithic transition and the genetics of populations in Europe. Princeton, NJ: Princeton University Press. Crossref, Google Scholar - 24
Renfrew C . 1990 Archaeology and language : the puzzle of indo-european origins. Cambridge, UK: Cambridge University Press. Google Scholar - 25
Özdoğan M . 2007 Von Zentralanatolien nach Europa. In Vor 12 000 jahren in anatolien: Die ältesten monumente der menschheit (ed.Lichter C ), p. 160. Stuttgart, Germany: Badisches Landesmuseum Karlsruhe. Google Scholar - 26
Özdoğan M . 2008 An alternative approach in tracing changes in demographic composition. In Neolithic demographic transition and its consequences (edsBocquet-Appel J, Bar-Yosef O ), pp. 139–178. Heidelberg, Germany: Springer. Google Scholar - 27
Horejs B, Milić B, Ostmann F, Thanheiser U, Weninger B, Galik A . 2015 The Aegean in the early 7th millennium BC: maritime networks and colonization. J. World Prehistory 28, 289–330. (doi:10.1007/s10963-015-9090-8) Crossref, PubMed, Web of Science, Google Scholar - 28
Zvelebil M . 2001 The agricultural transition and the origins of Neolithic society in Europe. Doc. Praehist. 28, 1–26. (doi:10.4312/dp.28.1) Crossref, Google Scholar - 29
Whittle AWR . 1996 Europe in the neolithic: the creation of new worlds. Cambridge, UK: Cambridge University Press. Google Scholar - 30
Kılınç GM 2016 The demographic development of the first farmers in Anatolia. Curr. Biol. 26, 2659–2666. (doi:10.1016/j.cub.2016.07.057) Crossref, PubMed, Web of Science, Google Scholar - 31
Hofmanová Z 2016 Early farmers from across Europe directly descended from Neolithic Aegeans. Proc. Natl Acad. Sci. USA 113, 6886–6891. (doi:10.1073/pnas.1523951113) Crossref, PubMed, Web of Science, Google Scholar - 32
Lazaridis I 2016 Genomic insights into the origin of farming in the ancient Near East. Nature 536, 419–424. (doi:10.1038/nature19310) Crossref, PubMed, Web of Science, Google Scholar - 33
Mathieson I 2017 The Genomic History of Southeastern Europe. BioRxiv (doi:10.1101/135616) Google Scholar - 34
Jones ER, Zarina G, Moiseyev V, Lightfoot E, Nigst PR, Manica A, Pinhasi R, Bradley DG . 2017 The Neolithic transition in the Baltic was not driven by admixture with early European Farmers. Curr. Biol. 27, 576–582. (doi:10.1016/J.CUB.2016.12.060) Crossref, PubMed, Web of Science, Google Scholar - 35
Gamba C 2014 Genome flux and stasis in a five millennium transect of European prehistory. Nat. Commun. 5, 5257. (doi:10.1038/ncomms6257) Crossref, PubMed, Web of Science, Google Scholar - 36
Günther T 2015 Ancient genomes link early farmers from Atapuerca in Spain to modern-day Basques. Proc. Natl Acad. Sci. USA 112, 11 917–11 922. (doi:10.1073/pnas.1509851112) Crossref, Web of Science, Google Scholar - 37
Mathieson I 2015 Genome-wide patterns of selection in 230 ancient Eurasians. Nature 528, 499–503. (doi:10.1038/nature16152) Crossref, PubMed, Web of Science, Google Scholar - 38
Haak W 2015 Massive migration from the steppe was a source for Indo-European languages in Europe. Nature 522, 207–211. (doi:10.1038/nature14317) Crossref, PubMed, Web of Science, Google Scholar - 39
Skoglund P, Malmström H, Raghavan M, Storå J, Hall P, Willerslev E, Gilbert MTP, Götherström A, Jakobsson M. 2012 Origins and genetic legacy of Neolithic farmers and hunter-gatherers in Europe. Science 336, 466–469. (doi:10.1126/science.1216304) Crossref, PubMed, Web of Science, Google Scholar - 40
Skoglund P 2014 Genomic diversity and admixture differs for Stone-Age Scandinavian foragers and farmers. Science 344, 747–750. (doi:10.1126/science.1253448) Crossref, PubMed, Web of Science, Google Scholar - 41
Bollongino R, Nehlich O, Richards MP, Orschiedt J, Thomas MG, Sell C, Fajkosova Z, Powell A, Burger J . 2013 2000 years of parallel societies in Stone Age central Europe. Science 342, 479–481. (doi:10.1126/science.1245049) Crossref, PubMed, Web of Science, Google Scholar - 42
Lazaridis I 2014 Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature 513, 409–413. (doi:10.1038/nature13673) Crossref, PubMed, Web of Science, Google Scholar - 43
Olalde I 2014 Derived immune and ancestral pigmentation alleles in a 7000-year-old Mesolithic European. Nature 507, 225–228. (doi:10.1038/nature12960) Crossref, PubMed, Web of Science, Google Scholar - 44
Gallego-Llorente M 2016 The genetics of an early Neolithic pastoralist from the Zagros, Iran. Sci. Rep. 6, 31326. (doi:10.1038/srep31326) Crossref, PubMed, Web of Science, Google Scholar - 45
Broushaki F 2016 Early Neolithic genomes from the eastern Fertile Crescent. Science 353, 499–503. (doi:10.1126/science.aaf7943) Crossref, PubMed, Web of Science, Google Scholar - 46
Jones ER 2015 Upper Palaeolithic genomes reveal deep roots of modern Eurasians. Nat. Commun. 6, 8912. (doi:10.1038/ncomms9912) Crossref, PubMed, Web of Science, Google Scholar - 47
Li H, Durbin R . 2009 Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760. (doi:10.1093/bioinformatics/btp324) Crossref, PubMed, Web of Science, Google Scholar - 48
Kircher M . 2012 Analysis of high-throughput ancient DNA sequencing data. In Methods in molecular biology (Clifton, N.J.) (edsShapiro B, Hofreiter M ), pp. 197–228. New York, NY: Humana Press. Google Scholar - 49
Patterson N, Moorjani P, Luo Y, Mallick S, Rohland N, Zhan Y, Genschoreck T, Webster T, Reich D . 2012 Ancient admixture in human history. Genetics 192, 1065–1093. (doi:10.1534/genetics.112.145037) Crossref, PubMed, Web of Science, Google Scholar - 50
McVean GA 2012 An integrated map of genetic variation from 1,092 human genomes. Nature 491, 56–65. (doi:10.1038/nature11632) Crossref, PubMed, Web of Science, Google Scholar - 51
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R . 2009 The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079. (doi:10.1093/bioinformatics/btp352) Crossref, PubMed, Web of Science, Google Scholar - 52
Danecek P 2011 The variant call format and VCFtools. Bioinformatics 27, 2156–2158. (doi:10.1093/bioinformatics/btr330) Crossref, PubMed, Web of Science, Google Scholar - 53
Purcell S 2007 PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. (doi:10.1086/519795) Crossref, PubMed, Web of Science, Google Scholar - 54
Patterson N, Price AL, Reich D . 2006 Population structure and eigenanalysis. PLoS Genet. 2, e190. (doi:10.1371/journal.pgen.0020190) Crossref, PubMed, Web of Science, Google Scholar - 55
Benjamini Y, Yekutieli D . 2001 The control of the false discovery rate in multiple testing under dependency. Ann. Stat. 29, 1165–1188. (doi:10.1214/aos/1013699998) Crossref, Web of Science, Google Scholar - 56
Quinlan AR, Hall IM . 2010 BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842. (doi:10.1093/bioinformatics/btq033) Crossref, PubMed, Web of Science, Google Scholar - 57
Korneliussen TS, Albrechtsen A, Nielsen R. 2014 ANGSD: analysis of next generation sequencing data. BMC Bioinformatics 15, 356. (doi:10.1186/s12859-014-0356-4) Crossref, PubMed, Web of Science, Google Scholar - 58
Reich D 2012 Reconstructing native American population history. Nature 488, 370–374. (doi:10.1038/nature11258) Crossref, PubMed, Web of Science, Google Scholar - 59
Excoffier L, Foll M . 2011 fastsimcoal: a continuous-time coalescent simulator of genomic diversity under arbitrarily complex evolutionary scenarios. Bioinformatics 27, 1332–1334. (doi:10.1093/bioinformatics/btr124) Crossref, PubMed, Web of Science, Google Scholar - 60
Baird D 2013 Juniper smoke, skulls and wolves' tails: the Epipalaeolithic of the Anatolian plateau in its South-west Asian context; insights from Pınarbaşı. Levant 45, 175–209. (doi:10.1179/0075891413Z.00000000024) Crossref, Web of Science, Google Scholar - 61
Bıçakçı E, Godon M, Çakan YG . 2012 Tepecik-Çiftlik. In Neolithic In Turkey: New excavations and new research 3–central Turkey (edsÖzdoğan M, Başgelen N, Kuniholm P ), pp. 89–134. İstanbul, Turkey: Archaeology and Art Publications. Google Scholar - 62
Roodenberg SA, Gerritsen FA, Özbal R . 2013 Neolithic burials from Barcın Höyük: the 2007–2012 excavation seasons. Anatolica 39, 93–111. Google Scholar - 63
Eriksson A, Manica A . 2012 Effect of ancient population structure on the degree of polymorphism shared between modern human populations and ancient hominins. Proc. Natl Acad. Sci. USA 109, 13 956–13 960. (doi:10.1073/pnas.1200567109) Crossref, Web of Science, Google Scholar - 64
Düring B . 2011 The prehistory of Asia minor, from complex hunter-gatherers to early urban societies. Cambridge, UK: Cambridge University Press. Google Scholar - 65
Omrak A 2016 Genomic evidence establishes Anatolia as the source of the European Neolithic gene pool. Curr. Biol. 26, 270–275. (doi:10.1016/j.cub.2015.12.019) Crossref, PubMed, Web of Science, Google Scholar - 66
Şenkul Ç, Doğan U . 2013 Vegetation and climate of Anatolia and adjacent regions during the Last Glacial period. Quat. Int. 302, 110–122. (doi:10.1016/j.quaint.2012.04.006) Crossref, Web of Science, Google Scholar - 67
Sarıkaya MA, Çiner A . 2015 Late Pleistocene glaciations and paleoclimate of Turkey. Bull. Miner. Res. Explor. 151, 107–127. (doi:10.19111/bmre.35245) Google Scholar - 68
Shumilovskikh LS, Tarasov P, Arz HW, Fleitmann D, Marret F, Nowaczyk N, Plessen B, Schlütz F, Behling H . 2012 Vegetation and environmental dynamics in the southern Black Sea region since 18kyr BP derived from the marine core 22-GC3. Palaeogeogr. Palaeoclimatol. Palaeoecol. 337, 177–193. (doi:10.1016/j.palaeo.2012.04.015) Crossref, Web of Science, Google Scholar - 69
Valsecchi V, Sanchez Goñi MF, Londeix L . 2012 Vegetation dynamics in the Northeastern Mediterranean region during the past 23 000 yr: insights from a new pollen record from the Sea of Marmara. Clim. Past 8, 1941–1956. (doi:10.5194/cp-8-1941-2012) Crossref, Web of Science, Google Scholar - 70
Timmermann A, Friedrich T . 2016 Late Pleistocene climate drivers of early human migration. Nature 538, 92–95. (doi:10.1038/nature19365) Crossref, PubMed, Web of Science, Google Scholar - 71
Giampoudakis K, Marske KA, Borregaard MK, Ugan A, Singarayer JS, Valdes PJ, Rahbek C, Nogués-Bravo D . 2017 Niche dynamics of Palaeolithic modern humans during the settlement of the Palaearctic. Glob. Ecol. Biogeogr. 26, 359–370. (doi:10.1111/geb.12543) Crossref, Google Scholar - 72
Pereira JB 2017 Reconciling evidence from ancient and contemporary genomes: a major source for the European Neolithic within Mediterranean Europe. Proc. R. Soc. B 284, 20161976. (doi:10.1098/rspb.2016.1976) Link, Web of Science, Google Scholar - 73
Fu Q 2016 The genetic history of Ice Age Europe. Nature 534, 200–205. (doi:10.1038/nature17993) Crossref, PubMed, Web of Science, Google Scholar - 74
Yakar J . 2016 The nature and extent of Neolithic Anatolia's contribution to the emergence of farming communities in the Balkans: an overview. In South east Europe and anatolia in prehistory band 263: essays in honor of vassil nikolov on His 65th anniversary (edsBacvarov K, Gleser R ), pp. 25–68. Bonn, Germany: Verlag Dr. Rudolf Habelt GmbH. Google Scholar - 75
Düring BS . 2013 Breaking the bond: investigating the neolithic expansion in Asia Minor in the seventh millennium BC. J. World Prehistory 26, 75–100. (doi:10.1007/s10963-013-9065-6) Crossref, Web of Science, Google Scholar - 76
Thissen L . 2011 The Neolithic–Chalcolithic sequence in the SW anatolian lakes region. Doc. Praehist. 37, 269. (doi:10.4312/dp.37.23) Crossref, Google Scholar - 77
Borić D . 2007 Mesolithic-Neolithic interactions in the Danube Gorges. In Mesolithic/neolithic interactions in the balkans and in the middle Danube basin (edsKozłowski JK, Nowak M ), pp. 31–45. Oxford, UK: Archaeopress. Google Scholar - 78
Fuller DQ, Willcox G, Allaby RG . 2011 Cultivation and domestication had multiple origins: arguments against the core area hypothesis for the origins of agriculture in the Near East. World Archaeol. 43, 628–652. (doi:10.1080/00438243.2011.624747) Crossref, Web of Science, Google Scholar - 79
Bailey DW . 2000 Balkan prehistory: exclusion, incorporation and identity. London, UK: Routledge. Crossref, Google Scholar - 80
Perlès C, Takaoğlu T, Gratuze B . 2011 Melian obsidian in NW Turkey: evidence for early Neolithic trade. J. Field Archaeol. 36, 42–49. (doi:10.1179/009346910X12707321242313) Crossref, Web of Science, Google Scholar - 81
Milić M . 2014 PXRF characterisation of obsidian from central Anatolia, the Aegean and central Europe. J. Archaeol. Sci. 41, 285–296. (doi:10.1016/j.jas.2013.08.002) Crossref, Web of Science, Google Scholar - 82
Çilingiroğlu A, Çilingiroğlu Ç . 2007 Ulucak. In Anadolu'da uygarlığın doğuşu avrupaya yayılımı türkiye'de neolitik dönem, yeni kazılar, yeni bulgular (edsÖzdoğan M, Başgelen N ), pp. 361–372. İstanbul, Turkey: Archaeology and Art Publications. Google Scholar - 83
Efstratiou N, Biagi P, Karkanas E, Starnini E . 2013 A Late Palaeolithic site at Ouriakos (Limnos, Greece) in the north-eastern Aegean. See http://antiquity.ac.uk/projgall/efstratiou335/. Google Scholar - 84
Yalçınkaya I 1995 Les occupations tardiglaciaires du site d’Öküzini (Sudouest de la Turquie) Résultats préliminaires. L'Anthropologie 99, 562–583. Google Scholar - 85
Broodbank C . 2006 The origins and early development of Mediterranean maritime activity. J. Mediterr. Archaeol. 19, 199–230. (doi:10.1558//jmea.2006.v19i2.199) Crossref, Google Scholar - 86
Vigne J-D, Zazzo A, Saliège J-F, Poplin F, Guilaine J, Simmons A . 2009 Pre-Neolithic wild boar management and introduction to Cyprus more than 11 400 years ago. Proc. Natl Acad. Sci. USA 106, 16 135–16 138. (doi:10.1073/pnas.0905015106) Crossref, Web of Science, Google Scholar - 87
Peltenburg EJ, Wasse A . 2004 Neolithic revolution: new perspectives on southwest Asia in light of recent discoveries on Cyprus. Oxford, UK: Oxbow Books. Google Scholar - 88
Vigne J-D 2012 First wave of cultivators spread to Cyprus at least 10 600 y ago. Proc. Natl Acad. Sci. USA 109, 8445–8449. (doi:10.1073/pnas.1201693109) Crossref, PubMed, Web of Science, Google Scholar - 89
Bar-Yosef O . 2015 Southeast Turkish Neolithic: A view from the southern Levant. In The neolithic in Turkey 6 (edsÖzdoğan M, Başgelen N, Kuniholm P ), pp. 293–320. İstanbul, Turkey: Arkeoloji ve Sanat Yayıniarı. Google Scholar - 90
Takaoğlu T, Korkut T, Erdoğu B, Işin G . 2014 Archaeological evidence for 9th and 8th millennia BC at girmeler cave near tlos in SW turkey. Doc. Praehist. 41, 111–118. (doi:10.4312/dp.41.6) Crossref, Google Scholar - 91
Simmons AH . 2007 The neolithic revolution in the near east: transforming the human landscape. Tucson, AZ: University of Arizona Press. Google Scholar - 92
Sampson A . 2005 New evidence from the early pottery production stages in the Aegean Basin from the 9th to the 7th millennium cal BCE. In How did farming reach Europe? Anatolian-European relations from the second half of the 7th through the first half of the 6th millennium cal BC (ed.Lichter C ), pp. 131–141. İstanbul: Byzas 2. Google Scholar - 93
Sampson A. 2008 The cave of the cyclops: Mesolithic and Neolithic networks in the northern Aegean, Greece. Volume I: Intra-Site analysis, local industries, and regional site distribution. Philadelphia, PA: Institute for Aegean Prehistory Academic Press. Crossref, Google Scholar - 94
Sampson A (ed.). 2011 The cave of the cyclops: mesolithic and neolithic networks in the northern aegean, Greece, volume II: bone tool industries, dietary resources and the paleoenvironment, and archeometrical studies. Philadelphia, PA: INSTAP Academic Press. Google Scholar