Y-chromosome phylogeographic analysis of the Greek-Cypriot population reveals elements consistent with Neolithic and Bronze Age settlements
- Konstantinos Voskarides†1,
- Stéphane Mazières†2,
- Despina Hadjipanagi1,
- Julie Di Cristofaro2,
- Anastasia Ignatiou1,
- Charalambos Stefanou1,
- Roy J. King3,
- Peter A. Underhill4,
- Jacques Chiaroni2Email author and
- Constantinos Deltas1Email author
© Voskarides et al. 2016
Received: 26 November 2015
Accepted: 2 February 2016
Published: 11 February 2016
The archeological record indicates that the permanent settlement of Cyprus began with pioneering agriculturalists circa 11,000 years before present, (ca. 11,000 y BP). Subsequent colonization events followed, some recognized regionally. Here, we assess the Y-chromosome structure of Cyprus in context to regional populations and correlate it to phases of prehistoric colonization.
Analysis of haplotypes from 574 samples showed that island-wide substructure was barely significant in a spatial analysis of molecular variance (SAMOVA). However, analyses of molecular variance (AMOVA) of haplogroups using 92 binary markers genotyped in 629 Cypriots revealed that the proportion of variance among the districts was irregularly distributed. Principal component analysis (PCA) revealed potential genetic associations of Greek-Cypriots with neighbor populations. Contrasting haplogroups in the PCA were used as surrogates of parental populations. Admixture analyses suggested that the majority of G2a-P15 and R1b-M269 components were contributed by Anatolia and Levant sources, respectively, while Greece Balkans supplied the majority of E-V13 and J2a-M67. Haplotype-based expansion times were at historical levels suggestive of recent demography.
Analyses of Cypriot haplogroup data are consistent with two stages of prehistoric settlement. E-V13 and E-M34 are widespread, and PCA suggests sourcing them to the Balkans and Levant/Anatolia, respectively. The persistent pre-Greek component is represented by elements of G2-U5(xL30) haplogroups: U5*, PF3147, and L293. J2b-M205 may contribute also to the pre-Greek strata. The majority of R1b-Z2105 lineages occur in both the westernmost and easternmost districts. Distinctively, sub-haplogroup R1b- M589 occurs only in the east. The absence of R1b- M589 lineages in Crete and the Balkans and the presence in Asia Minor are compatible with Late Bronze Age influences from Anatolia rather than from Mycenaean Greeks.
KeywordsCyprus Y-chromosome Neolithic Bronze Age Migration
The island of Cyprus is located ca. 100 km from the northern Levant and Anatolia. Evidence from both site excavations and genetics support it being a threshold from which maritime colonists commenced entry to the Mediterranean basin and southeast Europe [1–3] as well as a recipient of different cultural traditions reflecting subsequent human migratory events.
Major phases of prehistoric settlement in Cyprus based on material culture are summarized in Additional file 1: Table S1. In brief, the presence of human activity commenced 13,000 years ago when Mesolithic hunter-gatherers appear in the archeological record  followed by colonists associated with Pre-Pottery Neolithic A (PPNA, 11,000–10,400 years before present (y BP)) and B (PPNB, 10,500–8800 y BP) traditions, respectively [4–6]. Subsequent influential episodes of settlement and commerce occurred during the Pottery Neolithic (7200–6000 y BP) and the Early Bronze Age or Philia Horizon (4400–3700 y BP) [4, 7, 8]. The Pottery Neolithic and Early Bronze Age settlements were found in the Northwest/West/South portions of the island . The Late Bronze Age settlements reflected maritime commerce in the Eastern Mediterranean and are concentrated in the East/Northeast regions of Cyprus. This was followed by widespread societal collapse throughout the eastern Mediterranean 3200 years ago  and recoveries associated with the Assyrian, Phoenician, Hellenistic, Roman, and Ottoman periods .
Some extant Cypriot Y-chromosome data exists but it is either restricted to particular haplogroups  or of limited phylogenetic resolution [11, 12]. Considering their language and customs, the Greek element is expected in Greek-Cypriots but little is known about the genetic constitution provided from earlier occupation periods. However, genome-wide studies indicate that the genetic affinity of Cyprus is nearest to current populations of the Levant [2, 13], and an analysis of ancient mitochondrial DNA (mtDNA) from PPNB era associated specimens from Syria identified both U* and K lineages (also present in modern Cypriots) as part of the pre-Bronze context . Here, we report a high-resolution analysis of over 600 Y-chromosomes from contemporary Greek-Cypriots throughout Cyprus, whereby we explore the hypothesis that the present-day male genetic diversity of Cyprus also retains some elements distributed prior to the Hellenic period, with the following objectives in mind: (i) How does the Cypriot population compare genetically with surrounding populations? (ii) Which Y-chromosomes may reflect the Greek versus the pre-Greek settlements of Cyprus?
All aspects of sample selection and anonymity protections were managed by the Greek-Cypriot coauthors at the Molecular Medicine Research Center (MMRC) of the University of Cyprus. The Cyprus National Bioethics Committee approved the research program and the informed consent process. Most samples were collected through organized blood donation activities island-wide. Supplement DNA samples came from the DNA Biobank at the MMRC of the University of Cyprus. Specifically, 105 DNA samples came from the Biobank resource and the remainder from volunteers recruited during blood donations. All donors gave their signed consent for analyzing their DNA anonymously, according to the procedures approved by the Cyprus National Bioethics Committee. Genomic DNA was isolated from the whole peripheral blood  from 629 healthy unrelated Greek-Cypriot adult males. We attempted to achieve adequate geographical coverage by sampling from 300 different locations of all six official districts of Cyprus with sample number proportional to population size of each district (Additional file 2: Figure S1). The six districts are Kyreneia, Nicosia, Pafos, Limassol, Larnaka, and Ammochostos (alternate name: Famagusta). Criteria for district membership were continuity of familial origin over at least two generations. Although north Cyprus has been occupied by Turkey since 1974, the northern samples shown in Additional file 2: Figure S1 reflect displaced Greek-Cypriot donors currently living in southern Cyprus but whose familial paternal ancestral heritage traces to the north.
Genotyping was conducted in both Nicosia and Marseille. Samples were screened hierarchically, using 92 biallelic markers whose specifications are given in Additional file 3: Table S2. The nomenclature used reflects the 2015 version reported in the International Society of Genetic Genealogy (ISOGG) resource, an open science marker aggregator that adheres to the mutation-based rules proposed by the Y Chromosome Consortium (YCC) . All J1-M497-derived chromosomes were categorized on the basis of DYS388 repeat length with ‘long’ being ≥15 repeats . All J2a-M530-derived chromosomes were fractionated on the basis of the DYS445 6, 10, and 11 tandem repeat alleles . The particular nucleotide sequence repeat element counted for the Y-STR loci is given in Additional file 3: Table S2, and the nomenclature used is that of .
We genotyped Y-STRs in 574 of the 629 samples, 572 of which pertain to haplogroups (E-M215, G-Page94, J-M304, R-M207) using 17 loci: DYS19, DYS385 a/b, DYS389 I, DYS389 II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635 (Y-GATA-C4), and Y-GATA-H4 in the AmpFLSTR® Yfiler® PCR Amplification Kit (Applied Biosystem ® Yfiler, Foster City, USA) following the manufacturer’s instructions. We also genotyped three additional Y-STRs: DYS388, DYS445, and DYS461 in the same 574 samples. A Y-STR based phylogenetic network of E-V13 haplotypes defined by 15 loci was constructed using the program Network 184.108.40.206 (Fluxus-Engineering) with the median joining algorithm. The 15 loci used were DYS19, DYS389 I, DYS389 II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635 (Y-GATA-C4), and Y-GATA-H4. For the multi-copy short tandem repeat or microsatellite (STR) DYS389I,II, the DYS389b value was DYS389I subtracted from DYS389II. We also included nine E-V13 haplotypes from Anatolian Greeks—Phokaia and Smyrna  for a Greek comparison population.
Population genetic structure
We examined the population structure within Cyprus by using four approaches: genetic distance, spatial autocorrelation analysis, spatial analysis of molecular variance (SAMOVA), and analysis of molecular variance (AMOVA). First, Nei and Takezaki’s genetic distance  between Cypriot Y-STR haplotypes were estimated with Arlequin  and plotted by multidimensional scaling (MDS) using R . Second, Cypriot Y-STR haplotypes and their geographical coordinates were used to detect possible patterns of isolation by distance within Cyprus through an autocorrelation spatial analysis with even distance classes of 15, 25, and 50 km with the GENALEX software . Third, Y-STR haplotypes and Y haplogroup frequencies were subject to a SAMOVA in order to examine the genetic variance among the six Cypriot districts . SAMOVA implements an objective analysis of the genetic variance to search for geographically homogeneous groups and those differentiated from each other by the highest proportion of variance. We then attempted to identify geographic locations of the sharpest gradients of genetic variation in a manner independent of the six administrative districts by means of a SAMOVA analysis of the STR haplotype data. For this, we partitioned Cyprus into 38 areas of equal sample size of 15 haplotypes. Numbers of areas per district are: Paphos: 7, Ammochostos: 7, Kyreneia: 3, LarnaKa: 4, Limassol: 11, Nicosia: 6. Mean distance between centers of two adjacent areas is 10.09 ±3.93 km. We then ran SAMOVA 2.0 from K = 2 to K = 10 groups, K = 15, and K = 20 groups.
SAMOVA and AMOVA among the six district of Cyprus based on Y-STR haplotypes and haplogroup frequencies (Hgp freq)
No. of groups
Percent of variation
Among populations within groups
Departing districts from the rest of Cyprus
Kyreneia. Ammochostos. Limassol
Kyreneia. Ammochostos. Pafos
Cyprus. 6 districts as a whole
(Nicosia. Pafos. Limassol). (Kyreneia. Ammochostos. Larnaka)
(Kyreneia). (Nicosia. Pafos. Limassol). (Ammochostos. Larnaka)
(Larnaka. Ammochostos). (Kyreneia. Limassol. Nicosia. Pafos)
To assess Cyprus in a broader regional context, we gathered Y-chromosome haplogroup frequencies from 36 populations totaling 4666 males, from relevant surrounding regions, Levant-Middle East (Lebanon, Iraq/Baghdad, Egypt, Asia Minor), Caucasus (Armenia), continental Balkans (Bulgaria, Albania, Croatia), Mediterranean northern coast (Italy, Greece), Central Europe (Czech Republic, Hungary), and Mediterranean Islands (Crete, Sicily, Sardinia, present study) [12, 18, 26–32]. Since the various studies differed in regards to haplogroup resolution, we standardized the data to 25 haplogroups as well as appropriately collapsing the Cypriot data (Additional file 4: Table S3).
We then investigated the genetic affinities among the populations by principal component analyses (PCA) using XLSTAT 7.5.2.
Admixture contribution to Cyprus
Based on the pre- and historical events in Cyprus (Additional file 1: Table S1), we attempted to measure the genetic contribution of Anatolian, Balkans, Greek, and Levantine parental populations to Cyprus. Towards this aim, we first used Y-chromosome haplogroup frequencies from Additional file 4: Table S3 to measure the mY estimators with 100 replications as implemented in ADMIX 2.0 . We then focused on E-V13, G-P15, I2-M423, J2a-M67, J2b-M12, and R1b-M269 haplogroups, that we view as the most appropriate Anatolian, Balkan (Danube basin, Greece), and Levant proxies reflecting Neolithic and Bronze Age human dispersals (e.g., [18, 30, 31, 34, 35]). We compiled lineage (i.e., coupled single-nucleotide polymorphism (SNP) and STR) data for these surrogates from [10, 12, 18, 20, 28, 31, 35–40] and performed ADMIX95  to measure the m estimators of Anatolia, Danube Balkans, Greece, and Levant to the six fore mentioned Cypriot haplogroups. The lineages were composed of 11 loci (DYS19, DYS388, DYS389I, DYS389B, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, and DYS439). All input files for both ADMIX packages were built using the AdFiT v1.7 tool . A high R 2 means that the allelic frequencies in the hybrid population can be explained by the allelic frequencies in each of the parental population.
We estimated the time of expansion (TIMEX) of E-V13, G2a-P15, I2-M423, J2a-M67, J2b-M12, and R1b-M269 lineages in Cyprus and their time of divergence (TD) from Anatolian, Balkans, Greek, and Levantine sources. TIMEX and TD date lineages assuming that the observed variance has arisen indigenously within a given population as a result of a unique founder-based migration from an outside source. In practice, this means measuring the variation between the source and sink population assuming that it has accumulated since the time of settlement.
For TIMEX, variation was measured from the mean variance of the abovementioned 11 microsatellites. TD was based on the square difference between the means of allele size incorporated in the genetic distance denoted (δμ)2  and implemented in POWERMARKER . Under the assumption of a single population splitting into two fully isolated groups, this genetic distance is supposed to increase linearly with time since divergence as (δμ)2 = 2ωτ, where ω is the mutation rate and τ is the number of generations since isolation. Note that these approaches are sensitive to both multiple founders during a particular migration as well as subsequent population gene flows, both of which inflate STR variance and the age estimate of the event. Since the choice of mutation is debatable (e.g., [45, 46]), we used two mutation rates for STR: 0.00069 per STR per generation  and pedigree mutation rate of 0.0021 with 95 % confidence interval limit (CIL) of 0.0006–0.0049 %  to set lower and upper bounds. We assumed a generation time of 25 years.
The complex genetic structure of Cyprus
In Cyprus, haplogroup G2a-U5 (12.9 %) is widely distributed. While only 0.3 % distribute to the U5* paragroup, the PF3147 component (2.5 %) includes lineages like L91 (1.3 %), also seen in Asia Minor and Crete  and attributed to reflect the early Neolithic settlement of Sardinia . The G2a-L293 lineage emerges directly from the U5 node and occurs at 2.4 %, a level similar to that observed in Anatolia and the Levant (Additional file 5: Figure S2). Within the L30 defined clade (7.5 %), all lineages with the exception of L497 also occur in Anatolia  including the L30* paragroup (0.2 %). Haplogroup G2a-M406 (2.4 %) occurs at similar levels in Asia Minor and also on Crete (1.9 %), one of the earliest (ca. 9000 y BP) known sites outside the Levant colonized by Neolithic peoples . Within the G2a-P303 portion (4.9 %), the as yet unresolved P303* paragroup is seen (3.5 %) with both the Central European and Northern Italian G2a-L497 and the G2a-M527 lineages occurring at 0.3 %. The overall frequency, the haplogroup G2a-P303, especially the U1 branch is highest (20–39 %) in populations of the southern and northwest Caucasus . Lastly, we detected G2a-CTS342 (0.6 %), a lineage that has been reported in Sardinia  and as well in ancient DNA from Asia Minor .
Unlike samples from the present day interior Levant, such as Palestinians, Bedouins, and Jordanians [17, 54], J1-M267 is less common in Cyprus at 6.5 %. Haplogroup I2a lineages, thought to have arisen from a post-Last Glacial Maximum refuge and now present in Balkans, Sardinia, and Northwestern continental Europe [35, 55], were observed at ca. 3.5 % in our sample. Overall, haplogroup R1 presence was 15.1 %. The total frequency of associated R1a-M449 and R1b-M415 sub-haplogroups was 4.5 % and 10.7 %, respectively. The paragroup R1b-M269*(xL23) lineage is present (2.5 %). Furthermore sub-branches reflecting distinctive European versus Asian divergences similarly occur in both R1a and R1b. Within R1b, the central/west Europe M412 constituent (2.2 %) is offset by the western/central/south Asia Z2105 fraction (5.4 %) that was previously reported as paragroup L23*(xM412) . Similarly in R1a, both the European Z282 component (3.0 %) and the counterpart Asian Z93 clade (1.1 %) occur. Notable is finding that none of the R1a-Z93 Cypriots carried the diagnostic Ashkenazi Levite DYS456 14 repeat YSTR allele . Lastly, traces of geographically remote B-M60, C-M130, L-M20, and Q-M346 lineages were detected, mainly in the Nicosia district.
Spatial analysis of Y-STR variance among the six Cypriot districts barely showed significant geographical structure at K = 4 with northern, eastern, and southern regions separating (4.68 %, p < 0.05, Table 1). The nearly featureless geographic structure of the haplotype data is reiterated by the non-significant spatial autocorrelation (Additional file 6: Figure S3) as well as a lack of genetic affinity with district affiliation displayed in the MDS plot (Additional file 7: Figure S4). On the contrary, dispersal of Y-chromosome haplogroups (Fig. 3) reached significance at K = 3 to 5 (range of percentage of variation, 0.39–0.66 %, p values <0.01, Table 1).
We then searched for shared STR haplotypes .While the majority of our 574 STR haplotypes are unique we observed 197 perfect matches, most of which are shared within the same district and to a lesser extent between adjacent districts, indicative of recent demographic growth. To assess the degree of possible local patterns of genetic diversity we conducted a SAMOVA using the distributions of 574 hts within a uniform grid of 38 areas across the island. Mainly single areas separated first (Additional file 8: Table S4). However at the K= 3 level (Fig. 3) we detected two clusters that separate from the remaining samples, one of which is composed of 3 coastal grid areas that is characterized by 17.8 % of J2a-M67 derived chromosomes. The other cluster, comprised of two grid areas in the center of the island, entirely lacks any J2a-M67. At higher levels, K=6 and above, significant zones of reduced variance are restricted to single grid areas indicative of recent growth.
Ancestry of male Cypriot lineages
Bootstrapped (100 iterations) mY and m estimators of admixture proportions of four parental populations to Cyprus
Variance, time of expansion (TIMEX in y BP), and time of divergence (TD in y BP) of six major male lineages in Cyprus as compared with Anatolian, Balkan, Greek, and Levantine genetic contributors
TD Danube Balkans
0.0021 [95 % CIL]
0.0021 [95 % CIL]
0.0021 [95 % CIL]
0.0021 [95 % CIL]
0.0021 [95 % CIL]
Early and subsequent haplogroup dispersals
The strong correspondence between geography and Y-chromosome binary haplogroups is well known . This feature is consistent with a link between the distribution of haplogroups and past human movements. However, the task of deconvoluting prehistoric gene flows from subsequent transformations owing to ensuing migrations, local differentiations, and recent demographic growth (e.g., ) overlaid upon previous ones is complicated. For example, more recent migrations may also contain older haplogroups. Despite rigorous geographically targeted sampling and time parallels with cultural traits , current attempts to link of modern Y-chromosome patterns to prehistoric events are preliminary and best viewed with prudence. Such interpretations will be reappraised in the future using a combination of approaches, including simulation modeling , ancient DNA (e.g., [53, 58, 62]), and assessment of haplogroups that coalesce near the time frames of interest. This latter strategy is plausible due to the development of elaborately branched SNP dense phylogenies with branches proportional to time [63–67]. While, for reasons summarized in , ambiguity currently exists regarding an established mutation rate to use for calibration, this uncertainty will narrow as additional pedigree and clan based “whole” single-copy X-degenerate sequences are analyzed [46, 68] and cross-checked by ancient DNA-based rate estimates (e.g., ). With these caveats in mind, we proceed to define the putative prehistoric roots of Cypriot male genetic diversity by: (a) identifying lineages representative of non-Greek genetic influences, (b) reporting statistical support for correlations between settlement zones and haplogroup frequencies, (c) taking guidance from preliminary temporal estimates reported in vanguard studies of Y-chromosome phylogenies with meaningful branch lengths (e.g., [63, 64, 66]), and (d) noting coherence of ancient DNA data with pertinent archeological context.
Genetic legacy/substratum of the aceramic Neolithic
The Neolithic transition has diffused a wide array of culture, economic strategies, and social changes spanning the Levant, the Caucasus, and Europe [70–72] including Mediterranean islands that served as both way stations and terminal settlements [3, 14, 72–74]. Previous Y-chromosome studies [30, 31, 34, 37, 75, 76] hypothesized that lineages from haplogroups G, J, E, and R1b-M269 would have accompanied this cultural expansion although high levels of I2 and E-V13 versus low levels of extant G and J in the Balkans raise the possibility that only vestiges of pioneering agriculturalists to southeast Europe remain . In Europe, certain sub-haplogroups of G and specifically E-V13 were detected in ancient DNA, including Linear Band Keramik (LBK) remains from Central Europe (ca. 8000 y BP), Epicardial skeletons from Iberia (7000 y BP), South of France Late Neolithic (5000 y BP), and a Tyrol specimen (5300 y BP) [77–80].
In Cyprus, the G2a-U5 assemblage has an overall frequency of 12.9 %. Notably, PF3147 and related lineages occur mostly in the insular Mediterranean and display frequencies consistent with a relic distribution (e.g., ). Haplogroup G-M406 (2.4 %) is widely distributed across the island, but highest near Khirokitia, an aceramic Neolithic site (Table 1) located on the southern coast of Larnaka district. Interestingly, the more deeply rooted sub-haplogroup G2a-L293 also occurs in Anatolia and northern Levant (Additional file 5: Figure S2), consistent with the PPNB crescent including Syrian areas whose maternal genetic legacy is coherent with the maritime movements of early farmers . In addition, recent archeological studies have demonstrated that Epipaleolithic hunter-fisher-foragers colonized the island first (11,000–13,000 y BP) [81, 82], and that these forays are located on the southern coast as well (Aetokremnos, Amathus, Klimonas, and Asprokremnos), on the opposite slope of the Troodos Mountains [6, 83, 84]. Considered together the lithic industry, chronology and locations of these first human settlements match the PPNA tradition from Levant . Cyprus would thus represent one of the first stops of this diffusion, bringing probably some G male lineages to continental and insular Europe [6, 81]. As far as J, E, and R are concerned, refined investigation of their derived lineages could match with further events, characterized by the intensification of commercial exchanges throughout the Levantine Sea and other regions discussed below.
Centripetal gene flows during pottery Neolithic
It has been hypothesized that J2b-M12 may have been associated with the Neolithic immigration of farmers to Greece . Haplogroup J2b-M12 splits into J2b-M205 and J2b-M241. Since J2b-M241 frequency distribution was already well characterized , we mapped the frequency distribution of J2b-M205 (Additional file 5: Figure S2).
Both genetic data from the literature and Syria (Chiaroni J, unpublished results) show a frequency peak of J2b-M205 in the Southern Levant in which the frequency decreases northwards with latitude (Pearson’s R 2 = 0.282, p value = 0.011). The J2b-M205 distribution coincides with the Neolithic crop package dissemination found in common between the Levant and Cyprus . Also J2b-M12 Td estimates (Table 3) coupled with the J2b-M205 distribution overlap significantly with the Pottery Neolithic to Early Bronze Age pattern of settlements in Nicosia, Pafos, Limassol, and Kyreneia (chi-square = 11.29, p < .00084). This suggests the possibility that Cyprus experienced a later (Pottery Neolithic) immigration from the Southern Levant.
Regarding J2a-M410, the most common M530 sublineages are J2a-Z387 (4.9 %) linked with the distinctive six repeat allele at the DYS445 short tandem repeat locus proposed  to represent a Neolithic expansion from Anatolia to Greece and Italy, a pattern similar to G2a-M406 , at 4.9 % and J2a-Page55*(xM67,M319,Z387) at 5.2 %. The J2a-M319 lineage, previously observed in Crete and the Levant [18, 85] is also present in Cyprus at 1.1 %. However, its Y-STR haplotype diversity (Additional file 9: Table S5) is considerably higher from that in Crete (variance, 0.279 versus 0.121).
Early Bronze Age: in search of metals
Recent insights from ancient DNA studies suggest the spread of genes during the Bronze Age involved J2a-M67 individuals who appear in the Central European plains during the Late Bronze Age . J2a-M67, proposed to represent both the Neolithic of Central Anatolia and the expansion of the Troia Maritime Culture in Northwestern Anatolia (13.5 %) , is also quite common in Cyprus (10.1 %). Testing the hypothesis that the origin of the Early Bronze Philia culture in Cyprus derives from Western/Northwest Anatolia, the distribution of J2a- M67 on Cyprus fits well (chi-square = 3.42, p < .032, one-tailed). J2a-Z489, present in Pafos and Northwest/Central Anatolia, may reflect Bronze Age immigration from Western Anatolia, the Philia phase, or mirror the Jewish population on the island from the Hellenistic/Roman Eras.
Late Bronze Age Cyprus and maritime trade
E-V13 is common in the Balkans and may mark some of the Greek demographic input to Cyprus from the Late Bronze Age through the Iron Age . Network analysis of 46 E-V13 haplotypes (Additional file 10: Figure S5) shows a discrete clustering of 15 samples suggestive of a sub-haplogroup (encircled with an oval). This cluster is characterized by DYS437 = 15 repeats not seen in the Anatolian Greek population, or in the Provence samples . The remaining 31 samples overlap with the Anatolian Greek E-V13 lineages. Given that the highest frequency of I2 is in the Balkans , we also propose that I2-M423 (1.9 %) and I2-M436 (1.3 %) lineages reflect Greek influence. Additionally, the presence of G2a-M527 and G2a-U1 is consistent with remnants of Greek heritage . E-V22 and E-M34 are common in the Southern Levant, Sicily, Algeria, and in Egypt and rare in Europe [27, 36]. These lineages, like J2b-M205, could mirror a Pottery Neolithic movement to Cyprus from the Southern Levant (Pearson R 2 coefficient of correlation of E- M34 to longitude: 0.164, p = 0.003).
Haplogroup R1b-L23 membership on Cyprus is predominantly R1b-Z2105 and parallels the M412 lineage that is prolific in central and west Europe . Whole Y-chromosome phylogeny-based estimates of the coalescent times for the M412 and Z2105 companion lineages (and/or their phylogenetic equivalents or nearest neighbors) are reported in the range of 5000–6000 years [63, 64, 87]. While these dates should be viewed as preliminary as discussed previously, all estimates postdate the earlier dates of the aceramic and early ceramic phases of Neolithic settlements. Further underpinning the inference that R1b-Z2105 lineages are plausibly associated with Bronze Age settlement is that the oldest R1b-Z2105 designated lineages detected in ancient DNA specimens occur in the steppe belt regions of Russia and radiometrically date to ca. 5000 years ago . In Cyprus, Z2105 has the opposite distribution from J2a-M67 and J2b-M205, concentrating in the East/Northeast regions of Cyprus (chi-square = 5.01, p < .0256). The Late Bronze Age cities of Enkomi, Kition, and Hala Sultan Tekke, found in the Larnaka and Ammochostos districts, may have received immigrants from Hittite/Luwian Anatolia involved in the trade of the Late Bronze Age Eastern Mediterranean.
Updating previous R1b-L23* data pertinent to Cyprus [18, 28, 56] reveals that while the presence of R1b-Z2105 in central and west Europe is minimal, conversely, it is informative in present-day Anatolia (10.2 %), Greece (7.0 %), Bulgaria (5.7 %), and Crete (3.1 %). While the emerging R1b-Z2105 substructure  has yet to be evaluated at population levels, one relevant sub-haplogroup defined by M589 is illuminating. The majority of Cypriot R1b-Z2105 lineages occur in both the westernmost and easternmost districts, Pafos and Ammochostos, respectively. Notably, of the two districts, the M589 sub-haplogroup occurs only in the east. The absence of M589 lineages in Crete and the Balkans and the presence in Asia Minor are compatible with a record of Late Bronze Age influences from Anatolia rather than from Mycenaean Greeks.
We report a comprehensive granular Y-chromosome portrait of modern-day Cyprus. Some structural elements are consistent with Neolithic settlements in the central and eastern Mediterranean. Ancient DNA surveys from continental Europe have revealed a subsequent major population replacement within the last five millennia, masking a main part of the pre-Bronze Age genetic substratum [62, 80, 87]. We also detect lineages compatible with Bronze Age communities and subsequent events on Cyprus. The unstructured character of Y-STR lineages within Cyprus, low genetic diversity of Cypriot E-V13, and little genetic differentiation with surrounding populations would support this view.
Regionally speaking, Cyprus occupies an intermediate position, both geographically and in terms of its Y-chromosome patterning, between Levant, Crete, Italian, and Anatolian/Caucasus populations. Notably, Greek populations show genetic similarities with groups from Balkans. The Greek influence while culturally and linguistically profound only represents a small number of Y-chromosomes common in the Balkans and Carpathian areas.
The pre-Greek influence is most plausibly encapsulated by the following G2a haplogroups: U5*, PF3147*, L91, L293, P303*, and CTS342. Notably, most of these lineages occur in Anatolian ancient DNA specimens over 8200 years old . In addition, some J2a lineages such as M67, M319, and J2b-M205 may also have contributed to the pre-Greek strata.
Taking a broad perspective, our results are consistent with the model that Cyprus was an early recipient of Levantine-based origin. Interestingly, it apparently remained relatively isolated until experiencing immigration during the Early Bronze Age from Anatolia by early copper metallurgists. This period of isolation and genetic stasis is reminiscent of islands like Sardinia  that seem more detached from subsequent vigorous demographic episodes experienced on the mainland.
Although our results report just contemporary patterns and more recent migrations can potentially overwhelm signals of earlier genetic strata, our data approximates perspectives from the archeological record and provides important contextual reference for future ancient DNA studies. The ongoing revolution in ancient genomics  heralds the credible opportunity to more comprehensively intersect the cultural and genetic histories of the area and underscores the need for well-attested ancient specimens at population scale sample sizes from Cyprus, Anatolia, the Levant, and southeast Europe.
Ethics approval and consent to participate
The Cyprus National Bioethics Committee approved the research program and the informed consent process. All donors gave their signed consent for analyzing their DNA anonymously, according to the procedures approved by the Cyprus National Bioethics Committee.
- ca. :
short tandem repeat or microsatellite
- y BP:
years before present
We express our gratitude to the male donors of Cyprus who participated in this study. We also thank the following individuals for their insights either during the conceptualization phase or assistance during the progression of the study:
Tristan Carter, Rolandos Katsiaounis+, Nicos Koureas, Vassos Karageorghis, Demetrius Somas, Kiki Papastavrou, Olympia Pavlidou, Petro’s Lazio, Eftychia Kyriacou, Panagiota Demosthenous, Argyro Argyrou, Ted Kandell, Ray Banks, Blood Collection Units (Nicosia, Limassol, Larnaka). We thank Alice A. Lin for genotyping L293 and M589 in relevant samples from [18, 31] and M205 from . We thank Professor Ornella Semino for the M205, L293, Z2105, and M589 results for Bulgaria . S.M. thanks Julien Corny for the assistance in the use of R.
This work was supported by internal University of Cyprus funds to CD, through chapters 3/311 and 3/346. P.A.U. was supported by SAP SP0#115016 grant to Prof. Carlos D. Bustamante.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Perlés C, editor. The Early Neolithic in Greece. Cambridge: The first farming communities in Europe; 2001.Google Scholar
- Behar DM, Yunusbayev B, Metspalu M, Metspalu E, Rosset S, Parik J, et al. The genome-wide structure of the Jewish people. Nature. 2010;466:238–42.View ArticlePubMedGoogle Scholar
- Paschou P, Drineas P, Yannaki E, Razou A, Kanaki K, Tsetsos F, et al. Maritime route of colonization of Europe. Proc Natl Acad Sci U S A. 2014;111:9211–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Knapp BA. The archaeology of Cyprus. In: From Earliest Prehistory through the Bronze Age. Cambridge: Cambridge University Press; 2013.Google Scholar
- Colledge S, Conolly J, Shennan S. Archaeobotanical evidence for the spread of farming in the eastern Mediterranean. Curr Anthropol. 2004;45.
- Vigne JD, Briois F, Zazzo A, Willcox G, Cucchi T, Thiebault S, et al. First wave of cultivators spread to Cyprus at least 10,600 y ago. Proc Natl Acad Sci U S A. 2012;109:8445–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Webb JM, Frankel D. Characterizing the Philia facies: material culture, chronology and the Bronze Age in Cyprus. Am J Archaeol. 1999;103:3–43.View ArticleGoogle Scholar
- Steel L. Cyprus before history. London: Bristol Classical Press; 2011.Google Scholar
- Siddall R. "Not a day without a line drawn”: pigments and painting techniques of Roman Artists. In: Magazine F, editor. Proceedings of the Royal Microscopical Society. 2006. p. 18–23.Google Scholar
- Cruciani F, La Fratta R, Trombetta B, Santolamazza P, Sellitto D, Colomb EB, et al. Tracing past human male movements in northern/eastern Africa and western Eurasia: new clues from Y-chromosomal haplogroups E-M78 and J-M12. Mol Biol Evol. 2007;24:1300–11.View ArticlePubMedGoogle Scholar
- Capelli C, Redhead N, Romano V, Cali F, Lefranc G, Delague V, et al. Population structure in the Mediterranean basin: a Y chromosome perspective. Ann Hum Genet. 2006;70:207–25.View ArticlePubMedGoogle Scholar
- Zalloua PA, Xue Y, Khalife J, Makhoul N, Debiane L, Platt DE, et al. Y-chromosomal diversity in Lebanon is structured by recent historical events. Am J Hum Genet. 2008;82:873–82.PubMed CentralView ArticlePubMedGoogle Scholar
- Haber M, Mezzavilla M, Xue Y, Comas D, Gasparini P, Zalloua P et al. Genetic evidence for an origin of the Armenians from Bronze Age mixing of multiple populations. Eur J Hum Genet. 2015:1–6. doi: 10.1038/ejhg.2015.206
- Fernandez E, Perez-Perez A, Gamba C, Prats E, Cuesta P, Anfruns J, et al. Ancient DNA analysis of 8000 B.C. near eastern farmers supports an early neolithic pioneer maritime colonization of Mainland Europe through Cyprus and the Aegean Islands. PLoS Genet. 2014;10:e1004401.PubMed CentralView ArticlePubMedGoogle Scholar
- Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16:1215.PubMed CentralView ArticlePubMedGoogle Scholar
- Consortium TYC. A nomenclature system for the tree of human Y-chromosomal binary haplogroups. Genome Res. 2002;12:339–48.PubMed CentralView ArticleGoogle Scholar
- Chiaroni J, King RJ, Underhill PA. Correlation of annual precipitation with human Y-chromosome diversity and the emergence of Neolithic agricultural and pastoral economies in the Fertile Crescent. Antiquity. 2008;82:281–9.View ArticleGoogle Scholar
- King RJ, Ozcan SS, Carter T, Kalfoglu E, Atasoy S, Triantaphyllidis C, et al. Differential Y-chromosome Anatolian influences on the Greek and Cretan Neolithic. Ann Hum Genet. 2008;72:205–14.View ArticlePubMedGoogle Scholar
- Gusmao L, Butler JM, Carracedo A, Gill P, Kayser M, Mayr WR, et al. DNA Commission of the International Society of Forensic Genetics (ISFG): an update of the recommendations on the use of Y-STRs in forensic analysis. Forensic Sci Int. 2006;157:187–97.View ArticlePubMedGoogle Scholar
- King RJ, Di Cristofaro J, Kouvatsi A, Triantaphyllidis C, Scheidel W, Myres NM, et al. The coming of the Greeks to Provence and Corsica: Y-chromosome models of archaic Greek colonization of the western Mediterranean. BMC Evol Biol. 2011;11:69.PubMed CentralView ArticlePubMedGoogle Scholar
- Nei M, Takezaki N. Estimation of genetic distances and phylogenetic trees from DNA analysis. Proc 5th World Cong Genet Appl Livstock Prod. 1983;21:405–12.Google Scholar
- Excoffier L, Laval G, Schneider S. Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005;1:47–50.Google Scholar
- Team RDC. R: a language and environment for statistical computing. 2150th ed. Vienna: Computing RFfS; 2011.Google Scholar
- Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012;28:2537–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Dupanloup I, Schneider S, Excoffier L. A simulated annealing approach to define the genetic structure of populations. Mol Ecol. 2002;11:2571–81.View ArticlePubMedGoogle Scholar
- Herrera KJ, Lowery RK, Hadden L, Calderon S, Chiou C, Yepiskoposyan L, et al. Neolithic patrilineal signals indicate that the Armenian plateau was repopulated by agriculturalists. Eur J Hum Genet. 2012;20:313–20.PubMed CentralView ArticlePubMedGoogle Scholar
- Luis JR, Rowold DJ, Regueiro M, Caeiro B, Cinnioglu C, Roseman C, et al. The Levant versus the Horn of Africa: evidence for bidirectional corridors of human migrations. Am J Hum Genet. 2004;74:532–44.PubMed CentralView ArticlePubMedGoogle Scholar
- Karachanak S, Grugni V, Fornarino S, Nesheva D, Al-Zahery N, Battaglia V, et al. Y-chromosome diversity in modern Bulgarians: new clues about their ancestry. PLoS ONE. 2013;8:e56779.PubMed CentralView ArticlePubMedGoogle Scholar
- Grugni V, Battaglia V, Hooshiar Kashani B, Parolo S, Al-Zahery N, Achilli A, et al. Ancient migratory events in the Middle East: new clues from the Y-chromosome variation of modern Iranians. PLoS ONE. 2012;7, e41252.PubMed CentralView ArticlePubMedGoogle Scholar
- Rootsi S, Myres NM, Lin AA, Jarve M, King RJ, Kutuev I, et al. Distinguishing the co-ancestries of haplogroup G Y-chromosomes in the populations of Europe and the Caucasus. Eur J Hum Genet. 2012;20:1275–82.PubMed CentralView ArticlePubMedGoogle Scholar
- Cinnioglu C, King R, Kivisild T, Kalfoglu E, Atasoy S, Cavalleri GL, et al. Excavating Y-chromosome haplotype strata in Anatolia. Hum Genet. 2004;114:127–48.View ArticlePubMedGoogle Scholar
- Boattini A, Martinez-Cruz B, Sarno S, Harmant C, Useli A, Sanz P, et al. Uniparental markers in Italy reveal a sex-biased genetic structure and different historical strata. PLoS ONE. 2013;8, e65441.PubMed CentralView ArticlePubMedGoogle Scholar
- Dupanloup I, Bertorelle G. Inferring admixture proportions from molecular data: extension to any number of parental populations. Mol Biol Evol. 2001;18:672–5.View ArticlePubMedGoogle Scholar
- Semino O, Magri C, Benuzzi G, Lin AA, Al-Zahery N, Battaglia V, et al. Origin, diffusion, and differentiation of Y-chromosome haplogroups E and J: inferences on the neolithization of Europe and later migratory events in the Mediterranean area. Am J Hum Genet. 2004;74:1023–34.PubMed CentralView ArticlePubMedGoogle Scholar
- Battaglia V, Fornarino S, Al-Zahery N, Olivieri A, Pala M, Myres NM, et al. Y-chromosomal evidence of the cultural diffusion of agriculture in Southeast Europe. Eur J Hum Genet. 2009;17:820–30.PubMed CentralView ArticlePubMedGoogle Scholar
- Arredi B, Poloni ES, Paracchini S, Zerjal T, Fathallah DM, Makrelouf M, et al. A predominantly neolithic origin for Y-chromosomal DNA variation in North Africa. Am J Hum Genet. 2004;75:338–45.PubMed CentralView ArticlePubMedGoogle Scholar
- Balaresque P, Bowden GR, Adams SM, Leung HY, King TE, Rosser ZH, et al. A predominantly neolithic origin for European paternal lineages. PLoS Biol. 2010;8, e1000285.PubMed CentralView ArticlePubMedGoogle Scholar
- Regueiro M, Rivera L, Damnjanovic T, Lukovic L, Milasin J, Herrera RJ. High levels of Paleolithic Y-chromosome lineages characterize Serbia. Gene. 2012;498:59–67.View ArticlePubMedGoogle Scholar
- El-Sibai M, Platt DE, Haber M, Xue Y, Youhanna SC, Wells RS, et al. Geographical structure of the Y-chromosomal genetic landscape of the Levant: a coastal-inland contrast. Ann Hum Genet. 2009;73:568–81.PubMed CentralView ArticlePubMedGoogle Scholar
- Chiaroni J, King RJ, Myres NM, Henn BM, Ducourneau A, Mitchell MJ, et al. The emergence of Y-chromosome haplogroup J1e among Arabic-speaking populations. Eur J Hum Genet. 2009;18:348–53.PubMed CentralView ArticlePubMedGoogle Scholar
- Chakraborty R. Gene identity in racial hybrids and estimation of admixture rates. In: Neel JV, Ahuja Y, editors. Genetic Microdifferentiation in Man and Others Animals. Delhi: Indian Anthropological Association, Delhi University, Anthropological Department; 1985. p. 171–80.Google Scholar
- Gourjon G, Degioanni A. AdFiT v1.7 (Admixture File Tool): input files creating tool for genetic admixture estimation software. Bulletins et Mémoires de la Société d'Anthropologie de Paris. 2009;21:223–9.Google Scholar
- Goldstein DB, Ruiz Linares A, Cavalli-Sforza LL, Feldman MW. Genetic absolute dating based on microsatellites and the origin of modern humans. Proc Natl Acad Sci U S A. 1995;92:6723–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu K, Muse SV. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005;21:2128–9.View ArticlePubMedGoogle Scholar
- Xue Y, Wang Q, Long Q, Ng BL, Swerdlow H, Burton J, et al. Human Y chromosome base-substitution mutation rate measured by direct sequencing in a deep-rooting pedigree. Curr Biol. 2009;19:1453–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Helgason A, Einarsson AW, Guethmundsdottir VB, Sigurethsson A, Gunnarsdottir ED, Jagadeesan A, et al. The Y-chromosome point mutation rate in humans. Nat Genet. 2015;47:453–7.View ArticlePubMedGoogle Scholar
- Zhivotovsky LA, Underhill PA, Cinnioglu C, Kayser M, Morar B, Kivisild T, et al. The effective mutation rate at Y chromosome short tandem repeats, with application to human population-divergence time. Am J Hum Genet. 2004;74:50–61.PubMed CentralView ArticlePubMedGoogle Scholar
- Heyer E, Puymirat J, Dieltjes P, Bakker E, de Knijff P. Estimating Y chromosome specific microsatellite mutation frequencies using deep rooting pedigrees. Hum Mol Genet. 1997;6:799–803.View ArticlePubMedGoogle Scholar
- Hassan HY, Underhill PA, Cavalli-Sforza LL, Ibrahim ME. Y-chromosome variation among Sudanese: restricted gene flow, concordance with language, geography, and history. Am J Phys Anthropol. 2008;137:316–23.View ArticlePubMedGoogle Scholar
- Hammer MF, Behar DM, Karafet TM, Mendez FL, Hallmark B, Erez T, et al. Extended Y chromosome haplotypes resolve multiple and unique lineages of the Jewish priesthood. Hum Genet. 2009;126:707–17.PubMed CentralView ArticlePubMedGoogle Scholar
- Flores C, Maca-Meyer N, Larruga JM, Cabrera VM, Karadsheh N, Gonzalez AM. Isolates in a corridor of migrations: a high-resolution analysis of Y-chromosome variation in Jordan. J Hum Genet. 2005;50:435–41.View ArticlePubMedGoogle Scholar
- Francalacci P, Morelli L, Angius A, Berutti R, Reinier F, Atzeni R, et al. Low-pass DNA sequencing of 1200 Sardinians reconstructs European Y-chromosome phylogeny. Science. 2013;341:565–9.View ArticlePubMedGoogle Scholar
- Mathieson I, Lazaridis I, Rohland N, Mallick S, Patterson N, Roodenberg SA, et al. Genome-wide patterns of selection in 230 ancient Eurasians. Nature. 2015.
- Tofanelli S, Ferri G, Bulayeva K, Caciagli L, Onofri V, Taglioli L, et al. J1-M267 Y lineage marks climate-driven pre-historical human displacements. Eur J Hum Genet. 2009;17:1520–4.PubMed CentralView ArticlePubMedGoogle Scholar
- Rootsi S, Magri C, Kivisild T, Benuzzi G, Help H, Bermisheva M, et al. Phylogeography of Y-chromosome haplogroup I reveals distinct domains of prehistoric gene flow in Europe. Am J Hum Genet. 2004;75:128–37.PubMed CentralView ArticlePubMedGoogle Scholar
- Myres NM, Rootsi S, Lin AA, Jarve M, King RJ, Kutuev I, et al. A major Y-chromosome haplogroup R1b Holocene era founder effect in Central and Western Europe. Eur J Hum Genet. 2011;19:95–101.PubMed CentralView ArticlePubMedGoogle Scholar
- Rootsi S, Behar DM, Jarve M, Lin AA, Myres NM, Passarelli B et al. Phylogenetic applications of whole Y-chromosome sequences and the Near Eastern origin of Ashkenazi Levites. Nat Commun. 2013;4:2928.PubMed CentralView ArticlePubMedGoogle Scholar
- Haak W, Lazaridis I, Patterson N, Rohland N, Mallick S, Llamas B, et al. Massive migration from the steppe was a source for Indo-European languages in Europe. Nature. 2015;522:207–11.View ArticlePubMedGoogle Scholar
- Jobling MA, Tyler-Smith C. The human Y chromosome: an evolutionary marker comes of age. Nat Rev Genet. 2003;4:598–612.View ArticlePubMedGoogle Scholar
- Keinan A, Clark AG. Recent explosive human population growth has resulted in an excess of rare genetic variants. Science. 2012;336:740–3.PubMed CentralView ArticlePubMedGoogle Scholar
- Schraiber JG, Akey JM. Methods and models for unravelling human evolutionary history. Nature Review Genetics. 2015;16:727–40.View ArticleGoogle Scholar
- Allentoft ME, Sikora M, Sjogren KG, Rasmussen S, Rasmussen M, Stenderup J, et al. Population genomics of Bronze Age Eurasia. Nature. 2015;522:167–72.View ArticlePubMedGoogle Scholar
- Hallast P, Batini C, Zadik D, Maisano Delser P, Wetton JH, Arroyo-Pardo E, et al. The Y-chromosome tree bursts into leaf: 13,000 high-confidence SNPs covering the majority of known clades. Mol Biol Evol. 2015;32:661–73.PubMed CentralView ArticlePubMedGoogle Scholar
- Karmin M, Saag L, Vicente M, Wilson Sayres MA, Jarve M, Talas UG, et al. A recent bottleneck of Y chromosome diversity coincides with a global change in culture. Genome Res. 2015;25:459–66.PubMed CentralView ArticlePubMedGoogle Scholar
- Poznik GD, Henn BM, Yee MC, Sliwerska E, Euskirchen GM, Lin AA, et al. Sequencing Y chromosomes resolves discrepancy in time to common ancestor of males versus females. Science. 2013;341:562–5.PubMed CentralView ArticlePubMedGoogle Scholar
- Scozzari R, Massaia A, Trombetta B, Bellusci G, Myres NM, Novelletto A, et al. An unbiased resource of novel SNP markers provides a new chronology for the human Y chromosome and reveals a deep phylogenetic structure in Africa. Genome Res. 2014;24:535–44.PubMed CentralView ArticlePubMedGoogle Scholar
- Wei W, Ayub Q, Chen Y, McCarthy S, Hou Y, Carbone I, et al. A calibrated human Y-chromosomal phylogeny based on resequencing. Genome Res. 2013;23:388–95.PubMed CentralView ArticlePubMedGoogle Scholar
- Balanovsky O, Zhabagin M, Agdzhoyan A, Chukhryaeva M, Zaporozhchenko V, Utevska O, et al. Deep phylogenetic analysis of haplogroup G1 provides estimates of SNP and STR mutation rates on the human Y-chromosome and reveals migrations of Iranic speakers. PLoS ONE. 2015;10, e0122968.PubMed CentralView ArticlePubMedGoogle Scholar
- Fu Q, Li H, Moorjani P, Jay F, Slepchenko SM, Bondarev AA, et al. Genome sequence of a 45,000-year-old modern human from western Siberia. Nature. 2014;514:445–9.View ArticlePubMedGoogle Scholar
- Bellwood P. First farmers: the origins of agricultural societies. Oxford: Blackwell Publishing; 2005.Google Scholar
- Pinhasi R, Thomas MG, Hofreiter M, Currat M, Burger J. The genetic history of Europeans. Trends Genet. 2012;28:496–505.View ArticlePubMedGoogle Scholar
- Pinhasi R, Fort J, Ammerman AJ. Tracing the origin and spread of agriculture in Europe. PLoS Biol. 2005;3, e410.PubMed CentralView ArticlePubMedGoogle Scholar
- Fernandez H, Hughes S, Vigne JD, Helmer D, Hodgins G, Miquel C, et al. Divergent mtDNA lineages of goats in an Early Neolithic site, far from the initial domestication areas. Proc Natl Acad Sci U S A. 2006;103:15375–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Zeder MA. Domestication and early agriculture in the Mediterranean Basin: origins, diffusion, and impact. Proc Natl Acad Sci U S A. 2008;105:11597–604.PubMed CentralView ArticlePubMedGoogle Scholar
- Hovhannisyan A, Khachatryan Z, Haber M, Hrechdakian P, Karafet T, Zalloua P, et al. Different waves and directions of Neolithic migrations in the Armenian Highland. Investig Genet. 2014;5:15.PubMed CentralView ArticlePubMedGoogle Scholar
- Brandt G, Szecsenyi-Nagy A, Roth C, Alt KW, Haak W. Human paleogenetics of Europe—the known knowns and the known unknowns. J Hum Evol. 2014;79:73–92.View ArticlePubMedGoogle Scholar
- Szécsényi-Nagy A, Brandt G, Keerl V, Jakucs Jn, Haak W, Möller-Rieker S, et al. Tracing the genetic origin of Europeans first farmers reveals insights into their social organization. Proc Biol Sci. 2015;22(282):1805.
- Lacan M, Keyser C, Ricaut FX, Brucato N, Duranthon F, Guilaine J, et al. Ancient DNA reveals male diffusion through the Neolithic Mediterranean route. Proc Natl Acad Sci U S A. 2011;108:9788–91.PubMed CentralView ArticlePubMedGoogle Scholar
- Lacan M, Keyser C, Ricaut FX, Brucato N, Tarrus J, Bosch A, et al. Ancient DNA suggests the leading role played by men in the Neolithic dissemination. Proc Natl Acad Sci U S A. 2011;108:18255–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Keller A, Graefen A, Ball M, Matzas M, Boisguerin V, Maixner F, et al. New insights into the Tyrolean Iceman's origin and phenotype as inferred by whole-genome sequencing. Nat Commun. 2012;3:698.View ArticlePubMedGoogle Scholar
- Simmons AH. Humans, island colonization and Pleistocene extinctions in the Mediterranean: the view from Akrotiri Aetokremnos. Cyprus Antiquity. 1991;65:857–69.View ArticleGoogle Scholar
- Simmons AH, editor. Akrotiri-Aetokremnos (Cyprus) 20 years later: an assessment of its significance. Island archaeology and the origins of seafaring in the eastern Mediterranean. Reggio Calabria: Proceedings of the Wenner Gren Workshop; 2012.Google Scholar
- Mandel RD, Simmons AH. Geoarchaeology of the Akrotiri Aetokremnos Rockshelter, Southern Cyprus. Geoarchaeology: An International Journal. 1997;12:567–605.View ArticleGoogle Scholar
- McCartney C, Croft P, Manning SW, Rosendhal S. Preliminary report on the 2008 EENC excavations at Agia Varvara-Asprokremnos and regional field survey. 2009.Google Scholar
- Shen P, Lavi T, Kivisild T, Chou V, Sengun D, Gefel D, et al. Reconstruction of patrilineages and matrilineages of Samaritans and other Israeli populations from Y-chromosome and mitochondrial DNA sequence variation. Hum Mutat. 2004;24:248–60.View ArticlePubMedGoogle Scholar
- Gamba C, Jones ER, Teasdale MD, McLaughlin RL, Gonzalez-Fortes G, Mattiangeli V, et al. Genome flux and stasis in a five millennium transect of European prehistory. Nat Commun. 2014;5:5257.PubMed CentralView ArticlePubMedGoogle Scholar
- Batini C, Hallast P, Zadik D, Delser PM, Benazzo A, Ghirotto S, et al. Large-scale recent expansion of European patrilineages shown by population resequencing. Nat Commun. 2015;6:7152.PubMed CentralView ArticlePubMedGoogle Scholar
- Der Sarkissian C, Allentoft ME, Avila-Arcos MC, Barnett R, Campos PF, Cappellini E, et al. Ancient genomics. Philos Trans R Soc Lond B Biol Sci. 2015;370:20130387.View ArticleGoogle Scholar
- Karafet TM, Mendez FL, Meilerman MB, Underhill PA, Zegura SL, Hammer MF. New binary polymorphisms reshape and increase resolution of the human Y chromosomal haplogroup tree. Genome Res. 2008;18:830–8.PubMed CentralView ArticlePubMedGoogle Scholar