CURBIO13061S0960-9822(16)30850-810.1016/j.cub.2016.07.057The Author(s)Figure 1Geographical Location and Timeline of Ancient Individuals Included in This StudyMap showing the geographical distribution and timeline showing the approximate log-scaled time period (BC) of the ancient individuals used in this study. The colors and symbols for each individual are same with the principal component analysis (PCA). The regions where the Neolithic first emerged and was established are shaded. See Figure S1 for deamination patterns, sequencing efficiency using different methods for the individuals sequenced in this study, and an mtDNA haplogroup network. See also Data S1.Figure 2Genetic Structure and Diversity of Central Anatolian Neolithic Populations(A) PCA on contemporary west Eurasian populations onto which a total of 85 ancient individuals are projected from this study and previous studies. See Table S1 for number of SNPs per individual. Neighboring modern populations and ancient Anatolian populations are shown encircled. Modern population names are in italics.(B) Conditional nucleotide diversity calculated as the average pairwise mismatches between individuals. Diversities for each group were calculated using the SNPs identified in Yoruba individuals. We used two individuals per group, which yields the highest number of SNPs. Western European, eastern European, Swedish, and Caucasus hunter-gatherers are represented as WHG, EHG, SHG, and CHG, respectively. The European early Neolithic population is denoted with EN. Note that the diversities calculated for CHG and WHG are overestimates, as the individuals representing CHG are separated by three millennia and those representing WHG are separated by >1,000 km (Supplemental Experimental Procedures; Table S1). The error bars represent ±2 SEMs.(C) Distribution of runs of homozygosity (ROH) for Loschbour (European Mesolithic), Bon002 (Anatolian Aceramic), Bar8 (Anatolian Pottery Neolithic), and Stuttgart (early European Neolithic).(D) Multidimensional scaling analysis based on the Weir and Cockerham’s Fst calculated between populations using transversions overlapping with African Yoruba individuals. See Data S3 for f3 statistics, D statistics, pairwise mismatch estimates, and Fst estimates; Figure S2 for outgroup f3 statistics with present-day populations; and Figure S3 for D statistics, mean pairwise f3 statistics, and MDS analysis based on pairwise f3 statistics.Figure 3Admixture Analysis and Genetic Affinities among Neolithic/Chalcolithic Populations(A) ADMIXTURE ancestry components (K = 10) for present-day world populations and for ancient individuals. Admixture fractions are shown on map for modern-day individuals and as bar charts for ancient individuals. See Figure S4 for K = 2 to K = 20 plots with all individuals. Western European, eastern European, Swedish, and Caucasus hunter-gatherers are represented as WHG, EHG, SHG, and CHG, respectively. European early, middle, and late Neolithic populations are denoted with EN, MN, and LN, respectively.(B–D) Distributions of D statistics calculated as (B) D(Denisova,Iceman;X,Kumtepe) and D(Denisova,Remedello;X,Kumtepe), (C) D(Denisova,CHG;X,Kumtepe) and D(Denisova,CHG;X,Iceman/Remedello), and (D) D(Denisova,Tepecik;X,Kumtepe) and D(Denisova,Tepecik;X,Iceman/Remedello), where X stands for an ancient Anatolian or European early Neolithic (EN) or middle Neolithic (MN) individual, indicated on the left-hand y axis. (See Figure S3 for a plot of D statistics of comparisons of CO1, EHG, and WHG.) In brief, D < 0 indicates higher genetic affinity between the test population (name indicated on the top) and X, and D > 0 indicates higher genetic affinity between the test population and the second population (name indicated on the right-hand y axis). In each comparison, lighter-color boxplots show all D statistics calculated using all available individuals in the populations compared, and darker-color boxplots show only nominally significant D statistics with |Z| ≥ 2. The numbers in the middle indicate the percentage of comparisons where the test population resembles the population indicated on the right-hand y axis (i.e., D > 0). See Data S3 for D statistics.Table 1Summary Statistics of the Sequencing Data for Nine Ancient IndividualsSampleGenome CoveragemtDNA CoverageRead Length (Mean)mtDNA HaplogroupGenetic SexBon0010.166654.60463.208U3XYBon0026.6882,379.09069.841K1aXXBon0040.243351.23470.703N1a1a1XYBon0050.03968.61571.021N1a1a1XXTep0010.02366.81280.863K1aXYTep0020.721730.83360.814K1a12aXXTep0030.694281.96360.849N1b1aXYTep0040.473391.60861.473N1a1a1XXTep0060.267259.87983.585N1a1a1XYSee Data S2 for summary statistics for each library and for SNPs used for haplogroup classification.ReportThe Demographic Development of the First Farmers in AnatoliaGülşah MerveKılınç1210AyçaOmrak310FüsunÖzer110TorstenGünther2Ali MetinBüyükkarakaya4ErhanBıçakçı5DouglasBaird6Handan MelikeDönertaş1AyshinGhalichi1ReyhanYaka1DilekKoptekin7Sinan CanAçan1PooryaParvizi1MajaKrzewińska2Evangelia A.Daskalaki23ErenYüncü1Nihan DilşadDağtaş1AndrewFairbairn8JessicaPearson6GökhanMustafaoğlu9Yılmaz SelimErdal4Yasin GökhanÇakan5İnciTogan1∗togan@metu.edu.trMehmetSomel1∗∗msomel@metu.edu.trJanStorå2∗∗∗jan.stora@ofl.su.seMattiasJakobsson3∗∗∗∗mattias.jakobsson@ebc.uu.seAndersGötherström211∗∗∗∗∗anders.gotherstrom@arklab.su.se1Department of Biological Sciences, Middle East Technical University, 06800 Ankara, Turkey2Department of Organismal Biology, Uppsala University, Norbyvägen 18C, 75236, Uppsala, Sweden3Department of Archaeology and Classical Studies, Stockholm University, Lilla Frescativaegen 7, 114 18 Stockholm, Sweden4Department of Anthropology, Hacettepe University, Beytepe, 06800 Ankara, Turkey5Department of Prehistory, Istanbul University, Laleli, 34134 Istanbul, Turkey6Department of Archaeology, Classics, and Egyptology, University of Liverpool, Liverpool L69 7WZ, UK7Department of Health Informatics, Middle East Technical University, 06800 Ankara, Turkey8School of Social Science, University of Queensland, St Lucia, Brisbane, QLD 4072, Australia9Department of Archaeology, Bülent Ecevit University, 67100 İncivez, Zonguldak, Turkey∗Corresponding author∗∗Corresponding author∗∗∗Corresponding author∗∗∗∗Corresponding author∗∗∗∗∗Corresponding author10Co-first author11Lead ContactPublished: August 4, 2016SummaryThe archaeological documentation of the development of sedentary farming societies in Anatolia is not yet mirrored by a genetic understanding of the human populations involved, in contrast to the spread of farming in Europe [1–3]. Sedentary farming communities emerged in parts of the Fertile Crescent during the tenth millennium and early ninth millennium calibrated (cal) BC and had appeared in central Anatolia by 8300 cal BC [4]. Farming spread into west Anatolia by the early seventh millennium cal BC and quasi-synchronously into Europe, although the timing and process of this movement remain unclear. Using genome sequence data that we generated from nine central Anatolian Neolithic individuals, we studied the transition period from early Aceramic (Pre-Pottery) to the later Pottery Neolithic, when farming expanded west of the Fertile Crescent. We find that genetic diversity in the earliest farmers was conspicuously low, on a par with European foraging groups. With the advent of the Pottery Neolithic, genetic variation within societies reached levels later found in early European farmers. Our results confirm that the earliest Neolithic central Anatolians belonged to the same gene pool as the first Neolithic migrants spreading into Europe. Further, genetic affinities between later Anatolian farmers and fourth to third millennium BC Chalcolithic south Europeans suggest an additional wave of Anatolian migrants, after the initial Neolithic spread but before the Yamnaya-related migrations. We propose that the earliest farming societies demographically resembled foragers and that only after regional gene flow and rising heterogeneity did the farming population expansions into Europe occur.Highlights•Pre-pottery farmers had low genetic diversity, akin to Mesolithic hunter-gatherers•Genetic diversity levels are higher in the subsequent Pottery Neolithic•Central Anatolian farmers belonged to the same gene pool as early European farmers•Copper Age genetic affinities suggest a second wave of Anatolian gene flowKılınç et al. study ancient genomes from the earliest farmers of central Anatolia, one of the first areas where farming appears outside the Fertile Crescent. Genetic diversity increases as the Neolithic develops, indicating rising mobility. Similarities between Anatolian and European farmers suggest two gene flow events from Anatolia into Europe.Keywordsancient DNAarchaeogenomicsNeolithicAnatoliapopulation geneticsgenetic diversityResults and DiscussionThe causes, effects, and mechanisms of the transition from foraging to farming in western Eurasia are key issues in understanding the development of our species, especially in understanding the development of larger, more dense, and more socially complex populations. Over the past decade, archaeogenetic studies have largely focused on processes that drove the spread of farming practices, particularly the introduction of farming and sedentism into Europe [2, 3, 5–9]. However, the demographic aspects of the transformation of forager communities in Southwest Asia into communities practicing substantial-scale mixed farming and the full extent of the role of Anatolian populations in the spread of farming into Europe have remained unclear. Here, we investigate human remains excavated from two different Neolithic settlements in central Anatolia, Boncuklu and Tepecik-Çiftlik, between circa (ca.) 8300 and 5800 calibrated (cal) BC to explore the demographic processes during the earliest (Aceramic) phase of the Neolithic transition, as well as the later Pottery Neolithic period in Anatolia.Archaeological records show that the Neolithic era in Anatolia spanned more than 3,000 years—from around 9500 cal BC to around 6000 cal BC [4]. Farming practices were first established in the Fertile Crescent in the tenth and early ninth millennium cal BC [10] and in central Anatolia by 8300 cal BC [11, 12], or possibly earlier [12]. Between ca. 8000 cal BC and 6600 cal BC, farming spread west of central Anatolia, reaching the Aegean coast before 6600 cal BC and northwest Anatolia by 6600 at the latest [13, 14]. Debate exists as to whether this may have been a slow, steady process over those 1,400 years or relatively rapid between ca. 7000 and 6600 cal BC. Boncuklu, the earliest Anatolian site in our sample, and with evidence of very early crop cultivation in central Anatolia, is a small settlement mound dating between ca. 8300 and 7500 cal BC in the Aceramic Neolithic [11]. The excavators suggest that the Boncuklu community consisted of indigenous foragers who adopted small-scale cultivation and possibly experimented with animal herding alongside substantial traditional foraging practices [4, 11]. Tepecik-Çiftlik is a village with mixed and complex plant and animal exploitation practices, including notable elements of farming, located in the volcanic Cappadocian region of central Anatolia, dating between ca. 7500 and 5800 cal BC, from the latter Pre-Pottery Neolithic into the Pottery Neolithic [15, 16]. The evidence from Tepecik-Çiftlik indicates more substantial scale mixed farming relative to Boncuklu, although both hunting and gathering played a part in plant and animal exploitation. Both Boncuklu and Tepecik-Çiftlik show evidence of significant scale regional and inter-regional interactions, in the Tepecik-Çiftlik case especially with communities in the Fertile Crescent possibly related to the widespread distribution of obsidian [11, 15, 16]. The differences in subsistence patterns between these two settlements reflect a larger regional pattern seen in several other Aceramic and Pottery Neolithic sites in Anatolia [4, 13].We investigated a total of nine ancient individuals excavated from Boncuklu (n = 4) and Tepecik-Çiftlik (n = 5) (Data S1). We generated genome sequence data from these individuals with a mean coverage between 0.03-fold and 6-fold per individual, using a combination of whole-genome capture and direct shotgun sequencing strategies (Supplemental Experimental Procedures; Table 1; Data S2; Figures S1A and S1B). We authenticated the sequence data using multiple well-established approaches (Supplemental Experimental Procedures; Data S1; Figure S1C). Mitochondrial genome coverages were between 66- and 2,379-fold (Table 1), and all five Tepecik-Çiftlik and three Boncuklu individuals carried the haplogroups previously found in Neolithic farmers in Europe (haplogroups K and N) (Table 1; Data S2; Figure S1D) [17]. One of the Boncuklu individuals carried the haplogroup U3, which has also been observed in a later northwest Anatolian (Pottery) Neolithic site, Barcın (Figure 1), and in early Neolithic European farmers [8, 17, 18], but not among Eurasian hunter-gatherers [19]. We identified four individuals as females and the other five as males (Table 1; Data S1).We analyzed the new sequence data in the context of published ancient genetic variation (Figure 1). To discover the genetic affinities among ancient and modern-day individuals, we carried out principal component analysis (PCA). We calculated the principal components from 55 modern-day west Eurasian populations and projected the Boncuklu and Tepecik-Çiftlik individuals, as well as 85 published ancient individuals (Supplemental Experimental Procedures; Table S1), onto the first two principal components (Figure 2A). All individuals from the central Anatolian Neolithic, both the Aceramic Boncuklu group and the Pottery Neolithic Tepecik-Çiftlik group, were positioned within the genetic variation of present day southern European populations, consistent with outgroup f3 statistics (Figure S2; Data S3). Our central Anatolian Neolithic individuals (Boncuklu and Tepecik-Çiftlik), together with later (Pottery) Neolithic and Chalcolithic (Copper Age) individuals from northwest Anatolia (Barcın, Menteşe, and Kumtepe) and with early and middle Neolithic individuals from Europe, formed a distinct cluster to the exclusion of hunter-gatherers from western and eastern Europe (WHG and EHG, respectively), Sweden (SHG), and the Caucasus (CHG) (Figure 2A). Consistent with the PCA, D tests confirmed a clustering of Neolithic and Chalcolithic Anatolians to the exclusion of hunter-gatherers from Europe and the Caucasus. Hunter-gatherers from Europe and the Caucasus also share more alleles with their own groups than with Neolithic Anatolians (Figure S3A; Data S3). Interestingly, although geographically close, the Anatolian Neolithic populations from different time phases each formed discrete but proximate clusters in the PCA. Boncuklu individuals, representing the earliest phase of the Neolithic transition on the central Anatolian plateau, clustered tightly together, implying low genetic diversity within the population. In contrast, Tepecik-Çiftlik individuals, representing the later phase of the Neolithic in central Anatolia, were positioned at a peripheral position within the whole cluster and displayed high within-group diversity (Figure 2A). Pairwise f3 statistics between populations also showed significant differentiation between Boncuklu and Tepecik-Çiftlik populations (permutation test p < 0.05) (Data S3).To directly gauge levels of genetic diversity in Anatolian Neolithic populations, we calculated conditional nucleotide diversity in Boncuklu, Tepecik-Çiftlik, and Barcın, as well as in European Neolithic and hunter-gatherer populations (Data S3). Herein, we restricted the analysis to transversions identified in Yoruba as in [5] to avoid ascertainment bias, sequencing errors, and post-mortem degradation effects (Supplemental Experimental Procedures; Table S1). The Boncuklu population had remarkably low diversity relative to later ancient Anatolian populations, Tepecik-Çiftlik and Barcın, and European early Neolithic individuals from Hungary (Figure 2B). Comparison of the mean pairwise f3 statistics within populations also supported this result, with conspicuously higher genetic similarity within the Boncuklu group compared to Barcın and Tepecik-Çiftlik (Figure S3B; Data S3; 100% jackknife support). We further investigated short and intermediate runs of homozygosity (0.5–1.6 Mb); this is an indicator of historical effective population size and is expected to be influenced by geographic isolation and bottlenecks, but not recent inbreeding [20]. Our highest quality genome, Bon002 of Boncuklu, had 30% fewer such runs than the central European forager Loschbour, but 25%–40% more such runs relative to high-quality genomes from the Pottery Neolithic, Bar8 of Barcın and Stuttgart of Germany (Supplemental Experimental Procedures; Figure 2C). This supports the notion of a small ancestral population size in the Boncuklu population.We further evaluated genetic differentiation among Boncuklu, Tepecik-Çiftlik, Barcın, European Mesolithic, and Neolithic populations by calculating Fst (Supplemental Experimental Procedures; Data S3). The results were consistent with the pattern of differentiation in the PCA; particularly, Boncuklu appeared to be distinct from both Tepecik-Çiftlik and Barcın (Fst = 0.020 and 0.030, respectively; Z > 4). A multidimensional scaling (MDS) plot summarizing pairwise Fst values revealed clustering of Tepecik-Çiftlik and Barcın with European Neolithic populations, whereas Boncuklu attained a peripheral location (Figure 2C). This peripheral location is most likely due to high genetic homogeneity and drift in Boncuklu, as such a pattern was not observed in an MDS analysis of mean f3 statistics (Figure S3C).We next conducted ADMIXTURE analysis [21], inferring ancestral clusters from modern-day worldwide populations and estimating the ancestry proportions of each ancient individual based on the inferred ancestral cluster allele frequencies (Figures 3A and S4). With ten clusters (K = 10), ancestry proportions of all Anatolian (Boncuklu, Tepecik-Çiftlik, Barcın, Menteşe, and Kumtepe) and European Neolithic individuals consisted of two components, a “northern component” associated with European hunter-gatherers (WHG, SHG, and EHG) and found in modern-day northern Europe at highest frequency (orange), and a “southern component” found in the modern-day Middle East and North Africa (gray). Notably, Boncuklu displayed lower amounts of this “southern component” compared to individuals from Tepecik-Çiftlik and Barcın (Mann-Whitney U test, p < 0.001; Data S3), implying an influx of “southern component” alleles into late Aceramic and/or Pottery Neolithic settlements in Anatolia. This finding was also in line with higher genetic diversity in the later Neolithic Anatolian populations compared to Boncuklu (Figures 2B and 2C). D statistics results revealed genetic affinity between Caucasus hunter-gatherers (CHGs) and one of the individuals from Tepecik-Çiftlik, Tep003, which was greater than the rest of the individuals from Tepecik-Çiftlik and other Neolithic individuals from central Anatolia, northwest Anatolia, and Europe (Data S3). An admixture graph fitted by modeling gene flow from CHG to Tep003 using TreeMix [22] further confirmed the genetic relationship between Tep003 and CHG individuals (admixture proportion = 0.012, p = 0.002) (Figure S3D). These results show the buildup of genetic diversity during the development of the Neolithic in Anatolia.We next used our data to investigate a more recent case of possible regional migration. Previous work [6] had noted genetic affinity between Kumtepe from northwest Anatolia and the Tyrolean Iceman [23] from northern Italy. We found that the three Remedello individuals from Chalcolithic northern Italy [24], largely contemporary and possibly genetically and culturally affiliated with the Iceman, also had high affinity to Kumtepe in D statistics (Figure 3B; Data S3). A similar tendency for Kumtepe allele sharing was seen for a Chalcolithic individual from Hungary, CO1 [7], but was non-significant (Figure S3E; Data S3). Intriguingly, the Iceman/Remedello group was more similar to Kumtepe than to Boncuklu, Barcın, Tepecik-Çiftlik, or European Neolithic individuals. We further found that both Kumtepe and the Iceman/Remedello group carried more CHG alleles than other Neolithic populations (Figure 3C). This pattern of additional CHG allele sharing simultaneously observed in Iceman/Remedello and in Kumtepe is not mirrored in convergent allele sharing with other European hunter-gatherers (Figures S3F and S3G). We also found that Tepecik-Çiftlik individuals were consistently closer to Iceman/Remedello and to Kumtepe than to any other Anatolian or European early Neolithic population, including their contemporary Barcın and the neighboring Boncuklu (Figure 3D). These results point to gene flow from an eastern source into Chalcolithic Kumtepe and later into Europe, which could have crossed central Anatolia already before the Chalcolithic.Archaeogenetic studies have shown the existence of two distinct Mesolithic hunter-gatherer gene pools in west Eurasia: hunter-gatherers from Europe, ranging from Iberia to Scandinavia and to the Urals, and hunter-gatherers from the Caucasus [3, 5, 25]. The whereabouts of the so-called “early/first European farmer” gene pool [3], however, had remained unclear. Here we show that the genomes of Aceramic and Pottery Neolithic populations in central Anatolia belonged to the same group as northwestern Neolithic Anatolians and the first European farmers but were distinct from European and Caucasus foragers. The adoption of farming in central Anatolia by indigenous foragers, as suggested for Boncuklu [4, 11], would safely link the “early/first European farmer” gene pool to Anatolian foragers. However, the full geographic range of this forager population still remains to be described.The low genetic diversity of the Boncuklu population, resembling the low diversity in European hunter-gatherers [5, 25] is interesting (Figures 2B and 2C). It suggests that the population sizes at the very early stages of the Neolithic were not different from those of hunter-gatherers. This accords well with the view of indigenous forager adoption of cultivation and possible local initiation of herding in central Anatolia [4, 11]. Nearly 1,500 years later, Tepecik-Çiftlik and Barcın, fully established Neolithic populations practicing mixed farming (and within 200 km east and 400 km northwest of Boncuklu, respectively), were significantly more diverse (Figure 2B). Part of this increased genetic diversity could be linked to (1) putative southern gene flow (Figure 3A) that could be related to the Aceramic Neolithic to Pottery Neolithic transition in the Neolithic Levant or could be related to widespread interactions in the late Aceramic Neolithic between central Anatolia and the Fertile Crescent in the late Pre-Pottery Neolithic B [26]; (2) migration from the east related to similar factors of inter-regional exchanges (Figure S3D); and (3) admixture among local populations. Southern and eastern gene flow into Tepecik-Çiftlik is consistent with the site’s presumed role as an obsidian hub and its cultural links with the Levant and might have started already before the Pottery Neolithic [15, 16]. For Barcın, these results are also in line with archaeological evidence indicating cultural influx from central Anatolia [27]. This diverse Neolithic population most likely served as one of the sources for the well-documented wave of Neolithic migration to Europe [8, 9].Post-Neolithic contacts between parts of Anatolia and central Europe are a matter of discussion. Genetic affinity between a Chalcolithic group in northwest Anatolia represented by Kum6 of Kumtepe and by a group represented by the Tyrolean Iceman was earlier explained by gene flow post-dating the earlier stages of the Neolithic in Europe [6]. But it has alternatively been interpreted as the Iceman representing a relic of the first migratory event from Anatolia [9]. As we have shown in this paper, individuals of the Chalcolithic Remedello group [24] from northern Italy also share strong affinity with Kumtepe. This pattern may be explained with one out of four scenarios: (1) Iceman/Remedello representing a relict population stemming from an early farmer migratory event, (2) late-Neolithic/Chalcolithic back-migration from central Europe into Anatolia, (3) a third source-population admixing with both the population represented by Iceman/Remedello and the population represented by Kumtepe, and (4) secondary late-Neolithic/Chalcolithic migration from Anatolia. Because the Iceman/Remedello group is genetically closer to Chalcolithic Kumtepe than to earlier Anatolian Neolithic populations, including Boncuklu and Barcın, the first scenario seems unlikely. The fact that both Iceman/Remedello and Kumtepe display shared drift with Caucasus hunter-gatherers, independent of the Bronze Age Yamnaya expansions [24, 28], also argues against Iceman/Remedello being a relict population. Second, as Kumtepe predates the Iceman/Remedello group by some 1,300 years, back migration is an unlikely explanation. Finally, the Tepecik-Çiftlik population shows significant affinity to the Iceman/Remedello group and Kumtepe relative to other Anatolian and European Neolithic populations (Figure 3D); but Tepecik-Çiftlik also predates Iceman/Remedello by approximately 3,000 years. This implies gene flow events from Tepecik-Çiftlik-related populations into the Kumtepe-related west Anatolian populations, as predicted by archaeological evidence [29], and further gene flow that reached northern Italy by the fourth millennium BC. We propose an additional, yet undescribed, gene flow process from Anatolia into Europe as a better explanation than a contribution from a hypothetical third source into Neolithic central Anatolia, Chalcolithic northwest Anatolia, and Chalcolithic central Europe. Thus, Neolithic population dynamics that initiated in the Anatolian region resulted in multiple waves of expansion and admixture in west Eurasia.Experimental ProceduresDNA was isolated from petrous bone and teeth samples of nine ancient individuals. Double-stranded libraries were prepared and sequenced on Illumina HiSeq2500 and X platforms. Paired-end reads were merged, and adapters were removed. Reads were mapped to the human reference genome version hg18 and hs37d5 using BWA 0.7.12 [30]. Published ancient genomes were also mapped with the same parameters. Data was authenticated using four different methods [31–34]. Mitochondrial haplogroups were discovered using PhyloTree and Haplofind [35, 36]. Biological sex was determined using the Ry method [2, 37]. Principal component analysis was conducted using Eigensoft [38], and model-based clustering was performed using ADMIXTURE [21]. For ADMIXTURE analysis, ancestral components were determined using modern populations, and cluster memberships of each ancient individual were then inferred on the basis of these ancestral allele frequencies as in [39]. Outgroup f3 statistics were computed using popstats.py (https://github.com/pontussk/popstats). D statistics were calculated using qpDstat program of ADMIXTOOLS [40]. For computation of conditional nucleotide diversity, two approximately contemporaneous individuals with the highest quality genomes were selected to represent each group, and the average number of mismatches per each site overlapping with African Yoruba population between two individuals was calculated as in [41]. Weir and Cockerham’s Fst was calculated using popstats.py (https://github.com/pontussk/popstats). Runs of homozygosity for four high-quality genomes were calculated using PLINK [42]. See the Supplemental Experimental Procedures for details.Author ContributionsA. Götherström, M.J., J.S., İ.T., and M.S. designed and supervised the study; A.O., F. Ö., R.Y., M.K., E.D., E.Y., and N.D.D. performed the experiments; G.M.K. analyzed population-genetic data with contributions from M.S., A.O., T.G., H.M.D., A. Ghalichi, D.K., S.C.A., P.P., R.Y., E.Y., and N.D.D.; D.B., A.M.B., A.F., J.P., G.M., Y.S.E., Y.G.Ç., and E.B. excavated the samples, performed osteological assessments, and provided archaeological interpretations; and G.M.K., A.O., F.Ö., T.G., D.B., A.M.B., R.Y., A. Ghalichi, E.B., A. Götherström, M.J., J.S., İ.T., and M.S. wrote the manuscript with input from all authors.AcknowledgmentsWe thank C. Bilgin, A. Birand, R. Özbal, C. Knüsel, M. Özdoğan, F.S. Quinto, and two anonymous referees for helpful suggestions; İ.D. Toker and the Konya and Niğde Museums for support; A. Munters for assistance; and EMBO (Short-Term Fellowship to G.M.K.), Hacettepe University Scientific Research Projects Coordination Unit (project no. 13G602003 to A.M.B), İstanbul University Scientific Research Projects Unit (project no. 52349 to E.B), Australian Research Council (grant no. DP120100969 to A.F.), British Academy Research Development Award (grant no. BR100077 to D.B.), British Institute at Ankara grants 2012–2015 (to D.B.), National Geographic GEFNE 1–11 (to D.B.), Wainwright Fund University of Oxford (to D.B.), TÜBİTAK (grant no. 114Z927 to M.S.), TÜBA (GEBİP Prize to M.S), Sci. Acad. Turkey (BAGEP Prize to M.S.), METU (internal grants to İ.T. and M.S.), and ERC (starting grant no. 311413 to M.J.). Sequencing was conducted at the Uppsala University SNP&SEQ Technology Platform. Computational analyses were carried out at the Swedish National Infrastructure for Computing (SNIC-UPPMAX, projects b2013236, b2013240, b2015307, and b2015364).Accession NumbersThe accession number for the genome data produced in this study is European Nucleotide Archive: PRJEB14675.Supplemental InformationSupplemental Information includes Supplemental Experimental Procedures, four figures, one table, and three data sets and can be found with this article online at http://dx.doi.org/10.1016/j.cub.2016.07.057.Supplemental InformationDocument S1. Supplemental Experimental Procedures, Figures S1–S4, and Table S1Data S1. Overview of the Data, Related to Figure 1Data S2. Sequencing Statistics for Each Library and Haplogroup Classification, Related to Table 1Data S3. All f3 Statistics, D Statistics, Conditional Nucleotide Diversity and Fst Estimates, Mann-Whitney U Test Results, and Permutation Test Results, Related to Figures 2 and 3Document S2. Article plus Supplemental InformationReferences1B.BramantiM.G.ThomasW.HaakM.UnterlaenderP.JoresK.TambetsI.Antanaitis-JacobsM.N.HaidleR.JankauskasC.-J.KindGenetic discontinuity between local hunter-gatherers and central Europe’s first farmersScience32620091371402P.SkoglundH.MalmströmM.RaghavanJ.StoråP.HallE.WillerslevM.T.P.GilbertA.GötherströmM.JakobssonOrigins and genetic legacy of Neolithic farmers and hunter-gatherers in EuropeScience33620124664693I.LazaridisN.PattersonA.MittnikG.RenaudS.MallickK.KirsanowP.H.SudmantJ.G.SchraiberS.CastellanoM.LipsonAncient human genomes suggest three ancestral populations for present-day EuropeansNature51320144094134D.BairdThe Late Epipaleolithic, Neolithic, and Chalcolithic of the Anatolian Plateau, 13,000–4000 BCD.PottsA Companion to the Archaeology of the Ancient Near East2012Wiley-Blackwell4314665P.SkoglundH.MalmströmA.OmrakM.RaghavanC.ValdioseraT.GüntherP.HallK.TambetsJ.ParikK.-G.SjögrenGenomic diversity and admixture differs for Stone-Age Scandinavian foragers and farmersScience34420147477506A.OmrakT.GüntherC.ValdioseraE.M.SvenssonH.MalmströmH.KiesewetterW.AylwardJ.StoråM.JakobssonA.GötherströmGenomic evidence establishes Anatolia as the source of the European neolithic gene poolCurr. Biol.2620162702757C.GambaE.R.JonesM.D.TeasdaleR.L.McLaughlinG.Gonzalez-FortesV.MattiangeliL.DomboróczkiI.KőváriI.PapA.AndersGenome flux and stasis in a five millennium transect of European prehistoryNat. Commun.5201452578I.MathiesonI.LazaridisN.RohlandS.MallickN.PattersonS.A.RoodenbergE.HarneyK.StewardsonD.FernandesM.NovakGenome-wide patterns of selection in 230 ancient EurasiansNature52820154995039Z.HofmanováS.KreutzerG.HellenthalC.SellY.DiekmannD.Díez-Del-MolinoL.van DorpS.LópezA.KousathanasV.LinkEarly farmers from across Europe directly descended from Neolithic AegeansProc. Natl. Acad. Sci. USA11320166886689110E.AsoutiD.Q.FullerA contextual approach to the emergence of agriculture in Southwest Asia: reconstructing early neolithic plant-food productionCurr. Anthropol.54201329934511D.BairdA.FairbairnL.MartinC.MiddletonThe Boncuklu project: the origins of sedentism, cultivation and herding in central AnatoliaM.ÖzdoğanN.BaşgelenP.KuniholmThe Neolithic in Turkey: New Excavations and New Research, Volume 3, Central Turkey2012Archaeology and Art Publications21924412M.C.StinerH.BuitenhuisG.DuruS.L.KuhnS.M.MentzerN.D.MunroN.PöllathJ.QuadeG.TsartsidouM.ÖzbaşaranA forager-herder trade-off, from broad-spectrum hunting to sheep management at Aşıklı Höyük, TurkeyProc. Natl. Acad. Sci. USA11120148404840913M.ÖzdoğanAnatolia: from the pre-pottery neolithic to the end of the early bronze age (10,500–2000 BCE)C.RenfrewP.BahnThe Cambridge World Prehistory2014Cambridge University Press1508154414B.WeningerL.ClareF.GerritsenB.HorejsR.KraußJ.LinstädterR.ÖzbalE.J.RohlingNeolithisation of the Aegean and Southeast Europe during the 6600–6000 calBC period of rapid climate changeDoc. Praehist.XLI201413215E.BıçakcıTepecik-Çiftlik kazısı 2010 yılı çalışmalarıH.DonmezO.Otgun33. Kazı Sonucları Toplantısı 1.Cilt2011Kültür Varlıkları ve Müzeler Genel Müdürlüğü698916E.BıçakçıM.GodonY.G.ÇakanTepecik-ÇiftlikM.OzdoğanN.BaşgelenP.KuniholmNeolithic in Turkey: New Excavations and New Research, Volume 3, Central Turkey2012Archaeology and Art Publications8913417G.BrandtW.HaakC.J.AdlerC.RothA.Szécsényi-NagyS.KarimniaS.Möller-RiekerH.MellerR.GanslmeierS.FriederichGenographic ConsortiumAncient DNA reveals key stages in the formation of central European mitochondrial genetic diversityScience342201325726118A.Szécsényi-NagyG.BrandtW.HaakV.KeerlJ.JakucsS.Möller-RiekerK.KöhlerB.G.MendeK.OrossT.MartonTracing the genetic origin of Europe’s first farmers reveals insights into their social organizationProc. Biol. Sci.28220152015033919Q.FuC.PosthM.HajdinjakM.PetrS.MallickD.FernandesA.FurtwänglerW.HaakM.MeyerA.MittnikThe genetic history of Ice Age EuropeNature534201620020520T.J.PembertonD.AbsherM.W.FeldmanR.M.MyersN.A.RosenbergJ.Z.LiGenomic patterns of homozygosity in worldwide human populationsAm. J. Hum. Genet.91201227529221D.H.AlexanderJ.NovembreK.LangeFast model-based estimation of ancestry in unrelated individualsGenome Res.1920091655166422J.K.PickrellJ.K.PritchardInference of population splits and mixtures from genome-wide allele frequency dataPLoS Genet.82012e100296723A.KellerA.GraefenM.BallM.MatzasV.BoisguerinF.MaixnerP.LeidingerC.BackesR.KhairatM.ForsterNew insights into the Tyrolean Iceman’s origin and phenotype as inferred by whole-genome sequencingNat. Commun.3201269824M.E.AllentoftM.SikoraK.-G.SjögrenS.RasmussenM.RasmussenJ.StenderupP.B.DamgaardH.SchroederT.AhlströmL.VinnerPopulation genomics of Bronze Age EurasiaNature522201516717225E.R.JonesG.Gonzalez-FortesS.ConnellV.SiskaA.ErikssonR.MartinianoR.L.McLaughlinM.Gallego LlorenteL.M.CassidyC.GambaUpper Palaeolithic genomes reveal deep roots of modern EurasiansNat. Commun.62015891226M.ÖzdoğanAn alternative approach in tracing changes in demographic compositionO.Bar-YosefJ.P.Bocquet-AppelThe Neolithic Demographic Transition and Its Consequences2008Springer13917827L.ThissenH.ÖzbalA.TürkekulF.GerritsenR.ÖzbalThe land of milk? Approaching dietary preferences of Late Neolithic communities in NW AnatoliaLeiden J. Pottery Stud.26201015717228W.HaakI.LazaridisN.PattersonN.RohlandS.MallickB.LlamasG.BrandtS.NordenfeltE.HarneyK.StewardsonMassive migration from the steppe was a source for Indo-European languages in EuropeNature522201520721129M.ÖzdoğanArchaeological evidence on the westward expansion of farming communities from Eastern Anatolia to the Aegean and the BalkansCurr. Anthropol.522011S415S43030H.LiR.DurbinFast and accurate short read alignment with Burrows-Wheeler transformBioinformatics2520091754176031P.SkoglundB.H.NorthoffM.V.ShunkovA.P.DereviankoS.PääboJ.KrauseM.JakobssonSeparating endogenous ancient DNA from modern day contamination in a Siberian NeandertalProc. Natl. Acad. Sci. USA11120142229223432Q.FuM.MeyerX.GaoU.StenzelH.A.BurbanoJ.KelsoS.PääboDNA analysis of an early modern human from Tianyuan Cave, ChinaProc. Natl. Acad. Sci. USA11020132223222733R.E.GreenA.S.MalaspinasJ.KrauseA.W.BriggsP.L.F.JohnsonC.UhlerM.MeyerJ.M.GoodT.MaricicU.StenzelA complete Neandertal mitochondrial genome sequence determined by high-throughput sequencingCell134200841642634M.RasmussenX.GuoY.WangK.E.LohmuellerS.RasmussenA.AlbrechtsenL.SkotteS.LindgreenM.MetspaluT.JombartAn Aboriginal Australian genome reveals separate human dispersals into AsiaScience3342011949835D.VianelloF.SeviniG.CastellaniL.LomartireM.CapriC.FranceschiHAPLOFIND: a new method for high-throughput mtDNA haplogroup assignmentHum. Mutat.3420131189119436M.van OvenM.KayserUpdated comprehensive phylogenetic tree of global human mitochondrial DNA variationHum. Mutat.302009E386E39437P.SkoglundJ.StoråA.GötherströmM.JakobssonAccurate sex identification of ancient human remains using DNA shotgun sequencingJ. Archaeol. Sci.4020134477448238N.PattersonA.L.PriceD.ReichPopulation structure and eigenanalysisPLoS Genet.22006e19039M.SikoraM.L.CarpenterA.Moreno-EstradaB.M.HennP.A.UnderhillF.Sánchez-QuintoI.ZaraM.PitzalisC.SidoreF.BusoneroPopulation genomic analysis of ancient and modern genomes yields new insights into the genetic ancestry of the Tyrolean Iceman and the genetic structure of EuropePLoS Genet.102014e100435340N.PattersonP.MoorjaniY.LuoS.MallickN.RohlandY.ZhanT.GenschoreckT.WebsterD.ReichAncient admixture in human historyGenetics19220121065109341T.GüntherC.ValdioseraH.MalmströmI.UreñaR.Rodriguez-VarelaÓ.O.SverrisdóttirE.A.DaskalakiP.SkoglundT.NaidooE.M.SvenssonAncient genomes link early farmers from Atapuerca in Spain to modern-day BasquesProc. Natl. Acad. Sci. USA1122015119171192242S.PurcellB.NealeK.Todd-BrownL.ThomasM.A.R.FerreiraD.BenderJ.MallerP.SklarP.I.W.de BakkerM.J.DalyP.C.ShamPLINK: a tool set for whole-genome association and population-based linkage analysesAm. J. Hum. Genet.812007559575