The coming of the Greeks to Provence and Corsica: Y-chromosome models of archaic Greek colonization of the western Mediterranean
© King et al; licensee BioMed Central Ltd. 2011
Received: 28 July 2010
Accepted: 14 March 2011
Published: 14 March 2011
The process of Greek colonization of the central and western Mediterranean during the Archaic and Classical Eras has been understudied from the perspective of population genetics. To investigate the Y chromosomal demography of Greek colonization in the western Mediterranean, Y-chromosome data consisting of 29 YSNPs and 37 YSTRs were compared from 51 subjects from Provence, 58 subjects from Smyrna and 31 subjects whose paternal ancestry derives from Asia Minor Phokaia, the ancestral embarkation port to the 6th century BCE Greek colonies of Massalia (Marseilles) and Alalie (Aleria, Corsica).
19% of the Phokaian and 12% of the Smyrnian representatives were derived for haplogroup E-V13, characteristic of the Greek and Balkan mainland, while 4% of the Provencal, 4.6% of East Corsican and 1.6% of West Corsican samples were derived for E-V13. An admixture analysis estimated that 17% of the Y-chromosomes of Provence may be attributed to Greek colonization. Using the following putative Neolithic Anatolian lineages: J2a-DYS445 = 6, G2a-M406 and J2a1b1-M92, the data predict a 0% Neolithic contribution to Provence from Anatolia. Estimates of colonial Greek vs. indigenous Celto-Ligurian demography predict a maximum of a 10% Greek contribution, suggesting a Greek male elite-dominant input into the Iron Age Provence population.
Given the origin of viniculture in Provence is ascribed to Massalia, these results suggest that E-V13 may trace the demographic and socio-cultural impact of Greek colonization in Mediterranean Europe, a contribution that appears to be considerably larger than that of a Neolithic pioneer colonization.
Few genetic studies have explored the Greek contribution to the modern populations of Italy and France. A recent study of Y-chromosome haplogroups in a Sicilian population  showed a major impact of presumptive Greek immigration to the island estimated by admixture analysis to be about 37% using a localized Balkan/Greek marker E-V13. This level of admixture was higher than that predicted by classical demographic studies. Previous Y-chromosome genetic studies of Phoenician colonization have demonstrated that haplogroup J2 frequency was amplified in regions containing the Phoenician colonies of Iberia and North Africa in comparison to areas not containing Phoenician colonies . However, these studies did not address the role of either Greek colonization or early Neolithic colonization of Western Europe.
Y-chromosome studies have investigated the contribution of various Y haplogroups to the spread of farming from the Near East to Europe [8–10]. Haplogroup J2 frequency has been correlated with aspects of the symbolic material culture of the Neolithic in Europe and the Near East (painted pottery and ceramic figurines)  and sub-Haplogroups of J2 have also been associated with the Neolithic colonization of mainland Greece, Crete and southern Italy . On the other hand, E-V13 appears to have originated in Greece or the southern Balkans [13, 14] and then spread to Sicily at high frequencies with the Greek colonization of the island. E-V13 is also found at low frequencies on the Anatolian mainland  and thus may be useful in teasing apart the relative contributions of Greek colonization (E-V13) from Early Neolithic colonization (J2) to Western Europe. In this report, a sampling of individuals whose ancestry traces to the Ionian Greek city of Phokaia will be compared through Y-chromosome genotyping to samples from the Aeolian/Ionian city of Smyrna and a set of samples from Provence. These data will reveal genetic patterning characteristic of 1) the Ionian foundation of Phokaia versus the Aeolian/Ionian foundation of Smyrna. 2) the relative Y chromosome contributions of Phokaian Greeks and local Anatolian/Neolithic and/or central Anatolian populations in these two Asia Minor Greek city-states and 3) the contribution of Greek and/or Neolithic Y-chromosomes to the demographic pattern of Provence.
Populations used for AMOVA, MDS and Admixture Analyses
West Asia Minor
Anatolian Greek source
West Asia Minor
Anatolian Greek source
Anatolian Neolithic source
Anatolian Neolithic Source
Provence (surname based)
Fst analysis of haplogroup frequency
AMOVA results according to language and geography
% Var. among
% Var. among
% Var. within
The dominant haplogroup of Provence is R1b-M269 at 58.8% (Figure 2). Also found in Provence is haplogroup E-V13 (3.9%) and J2a-DYS445 = 6 (7.8%). All the V13 derived samples are from western Provence along the Rhone, while all the J2a-DYS445 = 6 are from Var in eastern Provence.
Admixture proportions (mY).
Neolithic (Central/Mediterranean Anatolian)
This study presents the first genetic data on those Greeks whose ancestry traces to western Anatolia before the 1923 exchange with Turkey. The two sites: Phokaia and Smyrna have a long established historical record and represent somewhat different Archaic Greek dialects and regions. Archaic Smyrna, a small polis of approximately 6 hectares in size, perhaps containing 700 individuals, was initially Aeolic with a subsequent immigration of Ionic Greeks from nearby Kolophon [16, 17]. Phokaia was a larger Ionic city-state (50 hectares), containing an estimated 6000 individuals including its surrounding chora, its agricultural territory [3, 16]. Smyrna, on the other hand, being a smaller polis, may show evidence of indigenous Anatolian admixture likely from neighbouring Lydia  with higher frequencies of J2a-Page55 derived chromosomes.
The frequency of J2a-DYS445 = 6 in Phokaia (6.5%) is comparable to that of central Anatolia (5.5%). Interestingly, the Anatolian Greek samples derived for J2a with DYS445 = 6 have DYS391 = 9 repeats, while samples from central Anatolia and Antalya in Mediterranean Anatolia and Crete, either are equally mixed with DYS391 = 9/10 or dominated by J2a-DYS445 = 6 with DYS391 = 10 or more repeats. The similar frequencies of J2a-DYS445 = 6 in the Greek city-state and Anatolia make the marker less useful for detecting a pure Neolithic component in other regions; however, the separation by DYS391 offers some utility in teasing apart the relevant components.
In France, Massalia was the unique initial Greek colony founded by the Phokaians circa 600 BCE . The initial colony was small, likely 12 hectares in area, but rapidly expanded during the following century to 40 hectares . Thus, its initial population may have numbered 1000 to 1500 rapidly growing to 5000 people including its small hinterland chora, later cultivated in large part with vineyards. In contrast, the departments of Var, Vaucluse and Bouche-du-Rhone contain an area of 14,000 square kilometers. During the Roman period, according to Beloch  who estimates a density of 10 individuals per square km in northern Italy, the population might have numbered 140,000. Even with an earlier 600 BCE reduction in population density around Massalia, it is probable that the indigenous Ligurians may have numbered at least 50,000. This would have yielded a maximum of 10% Greek input to Provence, much lower than the estimated 20% Y-chromosome input. However, this increase in Y-chromosome admixture from Greece is in accord with the recent results from Sicily, which estimated a 37% Greece input, in accordance with the demographic estimate of [18, 19]. We acknowledge that population history of Provence has been influenced by additional demographic events besides the Neolithic and Greek colonization events. One potential confound is the impact of the Roman Empire. However in other regions well known to have been settled by Romans, e.g. England, southern Spain, Morocco and Sardinia, the frequency of E-V13 ranges from zero to 1% . The impact of Phoenicians is minimal since the frequency of E-V13 in Lebanon is zero out of 42 samples (unpublished results, OS). Thus the presence of E-V13 in the western Mediterranean is most likely driven by Greek colonists. Interestingly the female input, estimated using mtDNA data may be minimal in Provence. One mtDNA study of Var, showed a negligible Neolithic (Near Eastern and hence Greek) component to the mtDNA distribution of Var . Results from a single locus like the Y chromosome phylogeny must be interpreted cautiously since haplogroup designation and population are not absolutely equivalent. In addition founder effects, sex-biased reproduction, sexual selection can skew the interpretation of a population's history.
The Greeks of Massalia, between 500 BCE and 300 BCE, conquered a vast nearby area and set up satellite trading posts, settlements and forts. These sites included Monoikos (Monaco), Nikaia (Nice), Antipolis (Antibes), Olbia and Tauroeis [4, 5]. The Greeks from Massalia also engaged in a major trading network along the Mediterranean coast and up the Rhone evidenced by Massaliote wine amphora and other ceramics . Our data are consistent with a male-mediated asymmetric gene flow into the indigenous Celto-Ligurian populations of Iron Age Provence due to possibly differential mating practices, elite dominance or enslavement.
The island of Corsica contains E-V13 Y-chromosomes, particularly in the eastern portion of the island at a frequency of 4.6%. Eastern Corsica was the site of a major Phokaian colony, Alalie, and the E-V13 network pattern suggests overlap among the regions studied. On the other hand, using J2a-DYS445 = 6, G-M406 and J2a-M92, we detected a Neolithic (Anatolian), impact on the demography of east Provence. This may be a slight overestimate, since no J2a-M92 or G-M406 derived chromosomes were found in the Provence samples. That said, the predominant region in which J2a-DYS445 = 6 lineages are present is Var, situated near initial Neolithic impressed ware sites . West of Var, J2a-DYS445 = 6 frequency drops off precipitously suggesting the demographic impact of Neolithic colonists from Anatolia does exceed beyond this region. The western districts of Vaucluse and Bouches du Rhônes contains Mesolithic sites and later cardial Neolithic package . The high level of indigenous Basque admixture in Provence is consistent of a model of the cultural diffusion of agriculture. The lack of Y-chromosome Neolithic markers in west Provence suggests that the subsequent cardial Neolithic may reflect a cultural adoption of farming in this area.
The Greeks from both mainland Greece and Anatolia made a major contribution to the development of western European culture through their Mediterranean colonies (Italy, France, and Spain) during the Iron Age. Haplogroup E-V13 may trace the movement of the Ionian Greeks to key areas of France and Corsica that introduced viniculture to Western Europe . Further studies will help elucidate the relative contribution of the Greek and Neolithic migrations in other areas of the western Mediterranean.
Our population samples included a total of 89 male subjects, currently living in Greece, who trace their grand-paternal ancestry to either the area near Phokaia (n = 31) or Smyrna (n = 58) prior to the 1923 Exchange of Lausanne. In addition 323 males living throughout Corsica who trace their paternal ancestry to the island, and 51 subjects from villages near Neolithic sites in Provence who trace their grand-paternal ancestry to Provence and the Principality of Monaco were also studied. A total of 23 of the subjects from Provence villages were from the western departments of Provence: Vaucluse and Bouche-du- Rhone, while 28 apportioned to the eastern departments: Var and Alps-de-Haute-Provence or to Monaco. Regarding the new samples introduced in this study, the Anatolian Greek component was approved by the IRB of Aristotle University, Thessaloniki, Greece. The French samples were approved by the French Committee for the Protection of Persons in Biomedical Research (CCPPRB) and the entire French collection were also declared to and approved by the French Ministry of Higher Education and Research. All subjects gave their informed consent to participate in the study. The location of the Anatolian Greek, mainland Greek, Turkish and Basques samples are shown in Figure 1. In addition the locations of Massalia and its trading posts and the Greek city of Alalie in Corsica are indicated. Additionally, a description of populations analyzed in this study is summarized in Table 1.
All 89 samples from Anatolian Greeks were genotyped using 29 Y-chromosome binary polymorphisms in a sequential manner using Y tree branching patterns to infer upstream haplogroup status. The following binary markers were genotyped: YAP, M35, V13, M78, M123, M34, M102, M9, M70, M74, M198, M269, M304, M497, M12, M241, M205, Page55, M67, M92, M530, M258, M253, M436, M223, P37.2, M423, M406, and M530. M497, Page55, Page94, and M530 are newly listed SNPs whose specifications are listed in Additional file 3: Supplemental Table S3. Binary marker genotyping was done by RFLP assay, DHPLC or direct sequencing. Each of the Phokaia and Smyrna samples were typed at 37 YSTRs listed in Additional file 4: Supplemental Table S4.
The 51 samples from areas near Neolithic sites in Provence had derived alleles for the following markers: V13, M34, Page94, M253, M438, M497, M530, M67, M198 and M269. In order to compare the selected 51 Neolithic samples to a larger Provence set, 368 subjects from the departments of Provence: Var (n = 68), Bouche-du-Rhone (n = 209), Vaucluse (n = 60) and Alps-de-Haute-Provence (n = 31) whose surname was blindly determined to be of French origin, were genotyped for E-V13, M406, Page94, M423, M269 and all the following J-lineages: M304, Page55, M267, M12, M410, M67, M530 and M92 as well as ten YSTRs for the J and E-V13 derived samples. DYS445 was typed in M530 derived samples (Schrack B.E., Athey T.W., Wilson J.F., 2006, The American Society of Human Genetics. Abstract). The 323 samples from Corsica were only typed for E-V13.
An AMOVA  was performed using Arlequin 2000  to test the population affinities of the two Anatolian Greek samples to three mainland Greek samples (Nea Nikomedeia, Sesklo/Dimini, Lerna/Franchthi Cave), and four regions of Turkey (western Aegean, Marmara, central Anatolia and Mediterranean Turkey) . Furthermore a Multidimensional Scaling analysis (MDS) (SPSS 18.0) was performed using the Fst measure as a distance metric across the 9 populations. An AMOVA comparing the effects of geography (Asia Minor vs. Mainland Greece) and religion/language (Christian/Greek vs. Muslim/Turkish) was also calculated using these 9 populations.
To analyze the impact of the attested Greek colonization of Provence, an admixture analysis  was conducted using a Basque population (n = 116)  as an indigenous (non-Neolithic pre-Greek) source population and the Phokaia/Smyrna data as the Greek colonizing source represented by E-V13 frequency. As a signal of putative Neolithic immigration to Provence, central Anatolian and Mediterranean Turkey data  were used. Specifically the following markers M92, M406 and J2a-(DYS445 = 6) were chosen as indicative of Neolithic ancestry. The frequencies of M92 and J2a-(DYS445 = 6) in the Basque population were estimated from their YSTR pattern . In order to assess the degree of E-V13 affinity, a 8 loci YSTR network using Phokaia, Smyrna, Provence and Corsica samples was constructed . Networks were constructed by the median joining method using Network 22.214.171.124, where ε = 0 and microsatellite loci were weighted proportionally to the inverse of the repeat variance observed in each haplogroup. Coalescent times for E-V13 based on the following 8 loci DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393 and DYS439 were computed using the methodology of Zhivotovsky et al. as modified according to Sengupta et al.. A microsatellite evolutionary effective mutation rate of 6.9 × 10-4 per 25 years was used .
We thank all the men who donated DNA samples used in this study. We wish to thank Al Aburto and Bonnie Schrack for haplogroup J2a phylogenetic insights. Professor Joachim Hallmayer provided laboratory facilities for PAU and AAL. We thank Scott R. Woodward and the Sorenson Molecular Genealogy Foundation for providing support for AAL and PAU. This research was supported by a grant to RJK from the Stanford-France Interdisciplinary Program and Grant PRIN 2007 to AN.
- Murray O: Early Greece. 1993, Cambridge: Harvard University Press, 2Google Scholar
- Brun JP: Archaeologie du Vin et de l'Huile. 2004, Edited by Editions Errance, ParisGoogle Scholar
- Hansen MH, Nielson TH: An Inventory of Classical Poleis. 2004, Oxford University PressGoogle Scholar
- Dominguez AJ: Spain and France (including Corsica). An Inventory of Classical Poleis. Edited by: Hansen MH, Nielson TH. 2004, Oxford University PressGoogle Scholar
- Dietler M: The Iron Age in Mediterranean France: Colonial Encounters, Entanglements and Transformations. Journal of World Prehistory. 1997, 11: 269-358. 10.1007/BF02221134.View ArticleGoogle Scholar
- Di Gaetano C, Cerutti N, Crobu F, Robino C, Inturri S, Gino S, Guarrera S, Underhill PA, King RJ, Romano V, et al: Differential Greek and northern African migrations to Sicily are supported by genetic evidence from the Y chromosome. Eur J Hum Genet. 2009, 17 (1): 91-9. 10.1038/ejhg.2008.120.View ArticlePubMedPubMed CentralGoogle Scholar
- Zalloua PA, Platt DE, El Sibai M, Khalife J, Makhoul N, Haber M, Xue Y, Izaabel H, Bosch E, Adams SM, Arroyo E, et al: Genographic Consortium. Identifying genetic traces of historical expansions: Phoenician footprints in the Mediterranean. Am J Hum Genet. 2008, 83 (5): 633-42. 10.1016/j.ajhg.2008.10.012.View ArticlePubMedPubMed CentralGoogle Scholar
- Di Giacomo F, Luca F, Anagnou N, Ciavarella G, Corbo RM, Cresta M, Cucci F, Di Stasi L, Agostiano V, Giparaki M, et al: Clinal patterns of human Y chromosomal diversity in continental Italy and Greece are dominated by drift and founder effects. Mol Phylogenet Evol. 2003, 28: 387-395. 10.1016/S1055-7903(03)00016-2.View ArticlePubMedGoogle Scholar
- Semino O, Passarino G, Oefner P, Lin AA, Arbuzova S, Beckman LE, De Benedictis G, Francalacci P, Kouvatsi A, Limborska S, et al: The genetic legacy of Palaeolithic Homo sapiens sapiens in extant Europeans: a Y-chromosome perspective. Science. 2000, 290: 1155-1159. 10.1126/science.290.5494.1155.View ArticlePubMedGoogle Scholar
- Semino O, Magri C, Benuzzi G, Lin AA, Al-Zahery N, Battaglia V, Maccioni L, Triantaphyllidis C, Shen P, Oefner PJ, 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 Gen. 2004, 74: 1023-1034. 10.1086/386295.View ArticleGoogle Scholar
- King R, Underhill PA: Congruent distributions of Neolithic painted pottery and ceramic figurines with Y-chromosome lineages. Antiquity. 2002, 76: 707-14.View ArticleGoogle Scholar
- King R, Ozcan S, Carter T, Kalfo011Flu E, Atasoy S, Triantiphyllidis C, Kouvatsi A, Lin A, Chow CE, Zhivotovsky L, et al: Differential Y chromosome influences on the Greek and Cretan Neolithic. Ann Hum Genet. 2008, 72: 205-214. 10.1111/j.1469-1809.2007.00414.x.View ArticlePubMedGoogle Scholar
- Cruciani F, La Fratta R, Trombetta B, Santolamazza P, Sellitto D, Colomb EB, Dugoujon JM, Crivellaro F, Benincasa T, Pascone R, 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 (6): 1300-11. 10.1093/molbev/msm049.View ArticlePubMedGoogle Scholar
- Battaglia V, Fornarino S, Al-Zahery N, Olivieri A, Pala M, Myres NM, King RJ, Rootsi S, Marjanovic D, Primorac D, et al: Y-chromosomal evidence of the cultural diffusion of agriculture in Southeast Europe. Eur J Hum Genet. 2008, 17 (6): 820-30. 10.1038/ejhg.2008.249.View ArticlePubMedPubMed CentralGoogle Scholar
- Myres NM, Rootsi S, Lin AA, Järve M, King RJ, Kutue v I, Cabrera VM, Khusnutdinova EK, Pshenichnov A, Yunusbayev B, Balanovsky O, Balanovska E, Rudan P, Marian Baldovic M, Herrera RJ, Chiaroni J, Di Cristofaro J, Villems R, Kivisild T, Underhill PA: A major Y-chromosome haplogroup R1b Holocene era founder effect in central and western Europe. Eur J Hum Genet. 2011, 19: 95-101. 10.1038/ejhg.2010.146.View ArticlePubMedPubMed CentralGoogle Scholar
- Hansen MH: The Shotgun Method: The Demography of Ancient Greek City- State Culture. 2006, Columbia: University of Missouri PressGoogle Scholar
- Herodotus: The History. Grene, D trans. 1987, Chicago: Chicago University PressGoogle Scholar
- Beloch J: Die Bevolkerung der Griechisch-Romischen Welt. 1886, Leipzig, Verlag von Duncker & HumbotGoogle Scholar
- Scheidel W: The Greek Demographic Expansion: Models and Comparisons. Journal Hellenic Studies. 2003, 123: 120-140. 10.2307/3246263.View ArticleGoogle Scholar
- Dubut V, Chollet L, Murail P, Cartaul F, Beraud-Colomb E, Serre M, Mogentale-Profizi N: mt-DNA Polymorphisms in five French groups: importance of regional sampling. Eur J Hum Genet. 2004, 12: 293-300. 10.1038/sj.ejhg.5201145.View ArticlePubMedGoogle Scholar
- Binder D: Mesolithic and Neolithic Interactions in southern France and northern Italy: New data and current hypotheses. Europe's First Farmers. Edited by: Douglas Price T. 2000, Cambridge University Press, 117-143.View ArticleGoogle Scholar
- Myles L, Boyko AR, Owens CL, Brown PJ, Grassi F, Aradhya MK, Prins B, Reynolds A, Chia JM, Ware D, Bustamante CD, and Buckler ES: Genetic structure and domestication history of the grape. PNAS. 2011Google Scholar
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: applications to human mitochondrial DNA restriction data. Genetics. 1992, 131 (2): 479-91.PubMedPubMed CentralGoogle Scholar
- Schneider S, Roessli D, Excoffier L: Arlequin ver. 2.0: A software for Population Genetics Data Analysis. 2000, Geneva, Switzerland: University of GenevaGoogle Scholar
- Cinnioğlu C, King R, Kivisild T, Kalfoğlu E, Atasoy S, Cavalleri GL, Lillie AS, Roseman CC, Lin AA, Prince K, et al: Excavating Y-chromosome haplotype strata in Anatolia. Human Genetics. 2004, 114: 127-148. 10.1007/s00439-003-1031-4.View 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-675. 10.1093/oxfordjournals.molbev.a003847.View ArticlePubMedGoogle Scholar
- Adams SM, Bosch E, Balaresque PL, Ballereau SJ, Lee AC, Arroyo E, López-Parra AM, Aler M, Grifo MS, Brion M, et al: The genetic legacy of religious diversity and intolerance: paternal lineages of Christians, Jews, and Muslims in the Iberian Peninsula. Am J Hum Genet. 2008, 83 (6): 725-36. 10.1016/j.ajhg.2008.11.007.View ArticlePubMedPubMed CentralGoogle Scholar
- Bandelt HJ, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16 (1): 37-48.View ArticlePubMedGoogle Scholar
- Qamar R, Ayub Q, Mohyuddin A, et al: Y-chromosomal DNA variation in Pakistan. Am J Hum Genet. 2002, 70: 1107-1124. 10.1086/339929.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhivotovsky LA, Underhill PA, Cinnioglu C, et al: On the effective mutation rate at Ychromosome STRs with application to human population divergence time. Am J Hum Genet. 2004, 74: 54-61. 10.1086/380911.View ArticleGoogle Scholar
- Sengupta S, Zhivotosky LA, King R, et al: Polarity and temporality of high resolution Y chromosome distributions in India identify both indigenous and exogenous expansions and reveal minor genetic influence of Central Asian pastoralists. Am J of Hum Genet. 2006, 78: 202-221. 10.1086/499411.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.