Skip to main content

Y-chromosome phylogeographic analysis of the Greek-Cypriot population reveals elements consistent with Neolithic and Bronze Age settlements

Abstract

Background

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.

Results

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.

Conclusions

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.

Background

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 [13] 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 [4] 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 [46]. 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 [4]. 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 [8] and recoveries associated with the Assyrian, Phoenician, Hellenistic, Roman, and Ottoman periods [9].

Some extant Cypriot Y-chromosome data exists but it is either restricted to particular haplogroups [10] 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 [14]. 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?

Methods

Sample collection

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 [15] 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

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) [16]. All J1-M497-derived chromosomes were categorized on the basis of DYS388 repeat length with ‘long’ being ≥15 repeats [17]. All J2a-M530-derived chromosomes were fractionated on the basis of the DYS445 6, 10, and 11 tandem repeat alleles [18]. 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 [19].

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 4.6.1.1 (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 [20] 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 [21] between Cypriot Y-STR haplotypes were estimated with Arlequin [22] and plotted by multidimensional scaling (MDS) using R [23]. 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 [24]. 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 [25]. 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.

Last, Y-chromosome haplogroup frequencies were used to test four historical models through an analysis of molecular variance (AMOVA) [22]. In order to test if the proportion of variance between the different district populations under study was asymmetrically distributed in Cyprus, we started with an AMOVA analysis considering the six populations as one group. We then tested four groupings according to three models of settlement of Cyprus: coastal versus inland (i.e., more than 5 km from the seashore); earliest (Nicosia, Pafos, Limassol) versus Bronze Age occupation sites (Ammochostos, Larnaka, Kyreneia); and two submodels of the arrival of the Philia phase (4400–3700 y BP) with its distinctive pottery style, Red Polished (Ammochostos, Larnaka versus the rest; and Ammochostos, Larnaka versus Kyreneia and the rest) (Table 1 and Additional file 1: Table S1). In order to assess the archeological geographic subdivision between the Pottery Neolithic and Early Bronze Age settlements (densely populating Northwest/West/South Cyprus) from the Late Bronze Age settlements (East/Northeast Cyprus), we a priori divided the island into two regions: (1) Kyreneia/Nicosia/Pafos/Limassol versus (2) Ammochostos/Larnaka and performed 2 × 2 table chi-square comparisons of frequency distributions of specific Y-chromosome haplogroups.

Table 1 SAMOVA and AMOVA among the six district of Cyprus based on Y-STR haplotypes and haplogroup frequencies (Hgp freq)

Population relationships

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, 2632]. 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 [33]. 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, 3540] and performed ADMIX95 [41] 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 [42]. 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.

Date estimates

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 [43] and implemented in POWERMARKER [44]. 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 [47] and pedigree mutation rate of 0.0021 with 95 % confidence interval limit (CIL) of 0.0006–0.0049 % [48] to set lower and upper bounds. We assumed a generation time of 25 years.

Results

The complex genetic structure of Cyprus

Figure 1 presents the phylogenetic relationships and frequencies of the Y-chromosome lineages detected in the six districts of Cyprus. Like other populations in Anatolia and Lebanon, Cyprus exhibits a high level of haplogroup J2-M172 related diversity. J2a-M410 is the dominant Y-chromosome lineage, constituting 26.0 % of the overall Cypriot samples. J2b-M12/M102 splits into mainly J2b-M205 (5.9 %), frequent in Southern Levant (Additional file 5: Figure S2), and J2b-M241 (0.6 %), most frequent in Greece and the Balkans [20, 35]. Overall, the E-M35 haplogroup totals to 23.1 % and contains various E-M78 sub-haplogroups including E-V13 (7.3 %) that is common in Greece [10, 18, 35] and E-V22 (3.5 %), that is frequent in Egypt [10] and Sudan [49]. Another E-M35 related haplogroup, E-M34, previously reported in Asia Minor [31], Southern Levant [50, 51], and the Balkans [35] also was observed in Cyprus (10.3 %).

Fig. 1
figure 1

Topological relationship of Y-chromosome binary markers and their observed haplogroup frequencies (absolute and relative) in the six districts of Cyprus. Nomenclature used in that recommended by [89]. Common names of markers are shown along the branches whose lengths are uninformative with respect to time. The asterisk refers to the unresolved status of the paragroup beyond the specific polymorphism. Six markers shown in italic font were not genotyped but provide context. The following 18 binary markers (with their haplogroup affiliation) were also genotyped but displayed no derived alleles: (E) V42, V6, V92, V257, M81, (G) P20, P16.1, P16.2, (I) M26, (J2a) M318, M419, M322, (O) M175, (Q) M25, M3, M378, (R2) M479, and (R1b) U106

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 [30] and attributed to reflect the early Neolithic settlement of Sardinia [52]. 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 [30] 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 [18]. 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 [30]. Lastly, we detected G2a-CTS342 (0.6 %), a lineage that has been reported in Sardinia [52] and as well in ancient DNA from Asia Minor [53].

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) [56]. 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 [57]. Lastly, traces of geographically remote B-M60, C-M130, L-M20, and Q-M346 lineages were detected, mainly in the Nicosia district.

Island substructure

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.

To assess the degree of possible local patterns of genetic diversity shaped by recent demographic forces, we conducted a SAMOVA using the distributions of 574 haplotypes within a uniform grid across the island. From K to K + 1, single areas stood out from the rest of Cyprus (Additional file 8: Table S4). Wherever we detected statistically significant zones of reduced variance (red dots in Fig. 2, Additional file 8: Table S4), the amount of within-group variance made them split at K + 1. Small geographical patterns of reduced genetic variation remained in the Southwest and in the East, attributable to recent demography processes associated with genetic drift. Before we tested for any possible signals of correspondences between settlement history and predefined groups of districts, we first confirmed that male lineage distribution among the six Cypriot districts (Table 1) showed an uneven percentage of variation (0.61 %, p value <0.001). The inland/coast partition model did not significantly pull apart the variance (0.09 %, p value >0.05). Neither the Central South West versus East partition (Pottery Neolithic to Early Bronze Age versus Late Bronze Age occupation site model) nor the East versus the rest of Cyprus (one of the Philia submodels) were found significant (0.18 and 0.19, respectively, p > 0.05). However, the best clustering with a district was found for one of the two submodels concerning the arrival of the Philia phase, namely the Pottery Neolithic to Early Bronze versus Late Bronze Age occupation site model, when Kyreneia is considered separately (0.39 %, p value <0.05).

Fig. 2
figure 2

Location of Y-STR haplotypes grouped by SAMOVA and their evolution from K groups to K + 1. Proportion of variance between groups ranged from 10.39 to 15.87 (p < 0.0001)

PCA shown in Fig. 3 depicts the main genetic relationships between Cyprus and surrounding populations based on Y-chromosome haplogroup frequencies. Axes 1 and 2 contribute 31.5 % of the total variance. The distributions of the populations trend with geography (R 2 Pearson’s coefficient of correlation with latitude and longitude above 0.092, p < 0.05). On axis 1, Bulgarian, Czech, Balkan, Hungarian, and Greek groups stand apart from Caucasus and Near Eastern populations, while axis 2 separates the Italian and Near Eastern groups. Notably, Cyprus and Crete occupy a central position. Caucasus, southern Italian, Crete, and neighbor Arab-speaking Egyptian and Iraqi populations show closer genetic relationships with our Cypriot sample set. A recent autosomal survey also revealed genetic affinity between Cypriot and Caucasus individuals, probably dating back from the early Bronze Age [13]. Noticeably, the long-term affiliated closest population, mainland Greeks, is genetically more distant than the aforementioned populations, clustering rather with the Balkan and Bulgarian groups. Such similarity between Greece and northern Balkans could trace back to the emergence of the Starcevo culture [58] in early Neolithic (8500 y BP).

Fig. 3
figure 3

Plot of the two first coordinates from principal component analysis depicting the relationship of Cypriot with 35 other regional populations based on frequencies of 25 comparable Y-chromosome haplogroups as described in the “Methods” section. The separate vector image reflects the role of various haplogroups with the population distributions

Ancestry of male Cypriot lineages

Comparing the entire set of Y-chromosome haplogroups with those from regional populations surrounding Cyprus revealed a high Anatolian influence (mY = 66 %), followed by the Levant (mY = 24 %) then the Balkan regions (mY = 13 %, Table 2). A putative Roman contribution to Cyprus using data from Italy and Sicily [32] also showed negative values of mY (data not shown). A closer look at specific proxy lineages permitted us to dissect these results further (Table 2). Anatolia would have generated up to 83 % to the Cypriot G-P15 and up to a quarter (range, 22–25 %) of Cypriot J2a-M67, J2b-M12, and R1b-M269 related lineages. Danube Balkans would have provided most of the Cypriot J2b-M12 (67 %) and all Cypriot I2-M423 (99 %). Although, when using the entire set of Y-chromosome haplogroup frequencies, the composition of Cyprus can be explained by contributions from Anatolia, Balkans, and Levant, the actual Greek contribution stood out for the Cypriot E-V13 (87 %), J2a-M67 (74 %), R1b-M269 (48 %), and G-P15 (17 %) components. Lastly, Levant contributed up to 30 % of the Cypriot R1b-M269 and to a lesser extent regarding the Cypriot J lineages (3–8 %).

Table 2 Bootstrapped (100 iterations) mY and m estimators of admixture proportions of four parental populations to Cyprus

The pattern of structural variation in Cyprus points towards a model comprising two stages of expansion: an earlier expansion of G2a-P15, J2a-M67, and R1b-M269 (range, 11,600–13,800 y BP with a slow YSTR mutation rate ω; 3800–4500 BP with a fast ω), subsequently followed later by the expansion of E-V13, I2-M423, and J2b-M12 (slow ω, 4400–6600 y BP; fast ω 1500–4500 y BP) (Table 3). However, times of divergence of these lineages from current Anatolian, Danubian, Greek, and Levantine Y-STRs appeared more recent. Pre-historical divergence was observed for Cypriot G-P15 with Greece (3600 y BP), I2-M423 with Anatolia (4200 y BP) and Levant (9400 y BP) and J2b-M12 with Danube Balkans (3,500 y BP) and Levant (5100 y BP). Divergence of E-V13, J2a-M67, and R1b-M269 would have taken place in modern times (range, 300–2.200 y BP). Note that both M67 and M269 encapsulate high genetic variance, respectively, 0.352 and 0.320, but little genetic differentiation with Anatolia, Danube Balkans, Greece, and Levant. This suggests either an arrival of multiple diverse founders during the Bronze Age period or alternatively several subsequent flows from these regions. In addition, the use of a fast YSTRS mutation rate reduced the divergence to the historical era. While these YSTR-based estimates are consistent with two stages of settlement, the chronological framework regarding such expansions remains in doubt due to uncertainty regarding the appropriate average YSTR mutation rate for the 11 loci that compose the haplotype.

Table 3 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

Discussion

Early and subsequent haplogroup dispersals

The strong correspondence between geography and Y-chromosome binary haplogroups is well known [59]. 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., [60]) 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 [20], 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 [61], 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 [6367]. While, for reasons summarized in [63], 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., [69]). 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 [7072] including Mediterranean islands that served as both way stations and terminal settlements [3, 14, 7274]. 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 [35]. 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) [7780].

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., [80]). 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 [14]. 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 [6]. 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 [18]. Haplogroup J2b-M12 splits into J2b-M205 and J2b-M241. Since J2b-M241 frequency distribution was already well characterized [35], 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 [5]. 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 [18] to represent a Neolithic expansion from Anatolia to Greece and Italy, a pattern similar to G2a-M406 [30], 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 [86]. J2a-M67, proposed to represent both the Neolithic of Central Anatolia and the expansion of the Troia Maritime Culture in Northwestern Anatolia (13.5 %) [31], 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 [79]. 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 [20]. The remaining 31 samples overlap with the Anatolian Greek E-V13 lineages. Given that the highest frequency of I2 is in the Balkans [28], 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 [30]. 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 [56]. 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 [58]. 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 [64] 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.

Conclusions

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.

Support for the perspective that present-day male genetic picture of Cyprus is consistent with an early arrival of pre-historical lineages covered by layers of Y-chromosomes in more recent times is summarized by the following:

  1. (i)

    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.

  2. (ii)

    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 [53]. 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 [80] 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 [88] 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.

Abbreviations

ca. :

circa

SNP:

single-nucleotide polymorphism

STR:

short tandem repeat or microsatellite

y BP:

years before present

References

  1. Perlés C, editor. The Early Neolithic in Greece. Cambridge: The first farming communities in Europe; 2001.

    Google Scholar 

  2. 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.

    Article  CAS  PubMed  Google Scholar 

  3. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Knapp BA. The archaeology of Cyprus. In: From Earliest Prehistory through the Bronze Age. Cambridge: Cambridge University Press; 2013.

    Google Scholar 

  5. Colledge S, Conolly J, Shennan S. Archaeobotanical evidence for the spread of farming in the eastern Mediterranean. Curr Anthropol. 2004;45.

  6. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Webb JM, Frankel D. Characterizing the Philia facies: material culture, chronology and the Bronze Age in Cyprus. Am J Archaeol. 1999;103:3–43.

    Article  Google Scholar 

  8. Steel L. Cyprus before history. London: Bristol Classical Press; 2011.

    Google Scholar 

  9. 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 

  10. 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.

    Article  CAS  PubMed  Google Scholar 

  11. 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.

    Article  CAS  PubMed  Google Scholar 

  12. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. 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

  14. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16:1215.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Consortium TYC. A nomenclature system for the tree of human Y-chromosomal binary haplogroups. Genome Res. 2002;12:339–48.

    Article  PubMed Central  Google Scholar 

  17. 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.

    Article  Google Scholar 

  18. 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.

    Article  CAS  PubMed  Google Scholar 

  19. 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.

    Article  CAS  PubMed  Google Scholar 

  20. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  21. 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 

  22. 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.

    CAS  Google Scholar 

  23. Team RDC. R: a language and environment for statistical computing. 2150th ed. Vienna: Computing RFfS; 2011.

    Google Scholar 

  24. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Dupanloup I, Schneider S, Excoffier L. A simulated annealing approach to define the genetic structure of populations. Mol Ecol. 2002;11:2571–81.

    Article  CAS  PubMed  Google Scholar 

  26. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  27. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. 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.

    Article  PubMed  Google Scholar 

  32. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Dupanloup I, Bertorelle G. Inferring admixture proportions from molecular data: extension to any number of parental populations. Mol Biol Evol. 2001;18:672–5.

    Article  CAS  PubMed  Google Scholar 

  34. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  38. 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.

    Article  CAS  PubMed  Google Scholar 

  39. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  41. 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 

  42. 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 

  43. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Liu K, Muse SV. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics. 2005;21:2128–9.

    Article  CAS  PubMed  Google Scholar 

  45. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. 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.

    Article  CAS  PubMed  Google Scholar 

  47. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. 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.

    Article  CAS  PubMed  Google Scholar 

  49. 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.

    Article  PubMed  Google Scholar 

  50. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. 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.

    Article  PubMed  Google Scholar 

  52. 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.

    Article  CAS  PubMed  Google Scholar 

  53. 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.

  54. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  55. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  57. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  58. 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.

    Article  CAS  PubMed  Google Scholar 

  59. Jobling MA, Tyler-Smith C. The human Y chromosome: an evolutionary marker comes of age. Nat Rev Genet. 2003;4:598–612.

    Article  CAS  PubMed  Google Scholar 

  60. Keinan A, Clark AG. Recent explosive human population growth has resulted in an excess of rare genetic variants. Science. 2012;336:740–3.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  61. Schraiber JG, Akey JM. Methods and models for unravelling human evolutionary history. Nature Review Genetics. 2015;16:727–40.

    Article  CAS  Google Scholar 

  62. 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.

    Article  CAS  PubMed  Google Scholar 

  63. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  64. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  65. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  67. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  68. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  69. 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.

    Article  CAS  PubMed  Google Scholar 

  70. Bellwood P. First farmers: the origins of agricultural societies. Oxford: Blackwell Publishing; 2005.

    Google Scholar 

  71. Pinhasi R, Thomas MG, Hofreiter M, Currat M, Burger J. The genetic history of Europeans. Trends Genet. 2012;28:496–505.

    Article  CAS  PubMed  Google Scholar 

  72. Pinhasi R, Fort J, Ammerman AJ. Tracing the origin and spread of agriculture in Europe. PLoS Biol. 2005;3, e410.

    Article  PubMed Central  PubMed  Google Scholar 

  73. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  74. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  75. 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.

    Article  PubMed Central  PubMed  Google Scholar 

  76. 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.

    Article  PubMed  Google Scholar 

  77. 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.

  78. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  79. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  80. 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.

    Article  PubMed  Google Scholar 

  81. Simmons AH. Humans, island colonization and Pleistocene extinctions in the Mediterranean: the view from Akrotiri Aetokremnos. Cyprus Antiquity. 1991;65:857–69.

    Article  Google Scholar 

  82. 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 

  83. Mandel RD, Simmons AH. Geoarchaeology of the Akrotiri Aetokremnos Rockshelter, Southern Cyprus. Geoarchaeology: An International Journal. 1997;12:567–605.

    Article  Google Scholar 

  84. 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 

  85. 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.

    Article  CAS  PubMed  Google Scholar 

  86. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  87. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  88. 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.

    Article  Google Scholar 

  89. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

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 [20]. We thank Professor Ornella Semino for the M205, L293, Z2105, and M589 results for Bulgaria [28]. S.M. thanks Julien Corny for the assistance in the use of R.

Funding

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.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Jacques Chiaroni or Constantinos Deltas.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

KV collected the samples. KV and DH analyzed the data and drafted the manuscript. SM performed most of the statistical data analyses and drafted the manuscript. DH, with the assistance and contribution of AI and CS, isolated and genotyped all DNA samples, organized all data, and helped to draft the manuscript. JDC supervised the genotyping of STRs and helped to draft the manuscript. RJK participated in the study design, integrate archeological framework, and drafted and revised the manuscript. PAU contributed to the phylogenetic analysis, the interpretation of the results and edited earlier drafts, and revised the manuscript. JC created spatial frequency maps, did network analysis, provided institutional support, helped coordinate the Marseille team, and helped to revise the manuscript. CD and KV conceived, supervised the study, coordinated the Cypriot team, and helped to revise the manuscript. All authors read and approved the final version of the manuscript.

Additional files

Additional file 1: Table S1.

Chronological list of major archaeological strata (adapted from [4]) (DOCX 12 kb)

Additional file 2: Figure S1.

Geographic sampling locations within the six districts of Cyprus. The Troodos Mountains mainly overlap Pafos, Nicosia, and Limassol districts. (PDF 2414 kb)

Additional file 3: Table S2.

SNP and STR marker specifications. First sheet presents the genotyped SNPs and indels, method of detection, assigned haplogroup, alternate common names, rs number, chromosomal position (Human Reference GRCh37/hg19 assembly), and nucleotide substitution. Second sheet presents the STR name, repeat structure, and allele range in repeat numbers. (XLSX 19 kb)

Additional file 4: Table S3.

Y-chromosome haplogroup frequencies and populations included in the PCA of Fig. 3. (XLSX 18 kb)

Additional file 5: Figure S2.

Spatial frequency distributions of haplogroups J2b-M205 and G2-L293 were generated using Surfer 10 (Golden Software). (PDF 109 kb)

Additional file 6: Figure S3.

Spatial autocorrelation analysis of Cypriot Y-STR haplotypes (PDF 413 kb)

Additional file 7: Figure S4.

MDS of the Nei genetic distances between Cypriot Y-STR haplotypes and R script employed to plot the 574 × 574 matrix of genetic distances (PDF 265 kb)

Additional file 8: Table S4.

SAMOVA results amongst 38 groups of Y-STR haplotypes in Cyprus. FCT: proportion of variance between groups. *p < 0.05; ***p < 0.0001 (DOCX 14 kb)

Additional file 9: Table S5.

Y-chromosome haplogroup specific STR haplotypes for 574 Cypriot samples and haplotype diversity by district (XLSX 101 kb)

Additional file 10: Figure S5.

E-V13 haplotype network based on 15 Y-STR loci in 46 Greek-Cypriot and 9 Anatolian Greek (Phokaia and Smyrna) samples. Circles represent microsatellite haplotypes, the areas of the circles and sectors are proportional to haplotype frequency (smallest circle corresponds to one individual). Red oval delineates the cluster of haplotypes with DYS437 = 15 repeat allele (PDF 237 kb)

Rights and permissions

Open Access This 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Voskarides, K., Mazières, S., Hadjipanagi, D. et al. Y-chromosome phylogeographic analysis of the Greek-Cypriot population reveals elements consistent with Neolithic and Bronze Age settlements. Investig Genet 7, 1 (2016). https://doi.org/10.1186/s13323-016-0032-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13323-016-0032-8

Keywords