# A Predominantly Neolithic Origin for European Paternal Lineages

Patricia Balaresque, Georgina R Bowden, Susan M Adams, Ho-Yee Leung, Turi E King, Zoë H Rosser, Jane Goodwin, Jean-Paul Moisan, Christelle Richard, Ann Millward, et al.

► To cite this version:

Patricia Balaresque, Georgina R Bowden, Susan M Adams, Ho-Yee Leung, Turi E King, et al. A Predominantly Neolithic Origin for European Paternal Lineages. PLoS Biology, 2010, 8 (1), pp.e1000285. 10.1371/journal.pbio.1000285. hal-05203138

HAL Id: hal-05203138

https://hal.science/hal-05203138v1

Submitted on 7 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# A Predominantly Neolithic Origin for European Paternal Lineages

Patricia Balaresque¹, Georgina R. Bowden¹, Susan M. Adams¹, Ho-Yee Leung¹, Turi E. King¹, Zoë H. Rosser¹, Jane Goodwin², Jean-Paul Moisan³, Christelle Richard³, Ann Millward⁴, Andrew G. Demaine⁴, Guido Barbujani⁵, Carlo Previderè⁶, Ian J. Wilson⁷, Chris Tyler-Smith⁸, Mark A. Jobling¹*

¹Department of Genetics, University of Leicester, Leicester, United Kingdom, ²Ty Celyn, Maeshafod, Blaina, Gwent, United Kingdom, ³Laboratoire d'Etude du Polymorphisme de l'ADN, Faculté de Médecine, Nantes, France, ⁴Molecular Medicine Research Group, Peninsula Medical School, Universities of Exeter and Plymouth, Plymouth, United Kingdom, ⁵Dipartimento di Biologia ed Evoluzione, Università di Ferrara, Ferrara, Italy, ⁶Dipartimento di Medicina Legale e Sanità Pubblica, Università di Pavia, Pavia, Italy, ⁷Institute of Human Genetics, Newcastle University, Newcastle upon Tyne, United Kingdom, ⁸The Wellcome Trust Sanger Institute, Wellcome Trust Genome Campus, Hinxton, United Kingdom

## Abstract

The relative contributions to modern European populations of Paleolithic hunter-gatherers and Neolithic farmers from the Near East have been intensely debated. Haplogroup R1b1b2 (R-M269) is the commonest European Y-chromosomal lineage, increasing in frequency from east to west, and carried by 110 million European men. Previous studies suggested a Paleolithic origin, but here we show that the geographical distribution of its microsatellite diversity is best explained by spread from a single source in the Near East via Anatolia during the Neolithic. Taken with evidence on the origins of other haplogroups, this indicates that most European Y chromosomes originate in the Neolithic expansion. This reinterpretation makes Europe a prime example of how technological and cultural change is linked with the expansion of a Y-chromosomal lineage, and the contrast of this pattern with that shown by maternally inherited mitochondrial DNA suggests a unique role for males in the transition.

**Citation:** Balaresque P, Bowden GR, Adams SM, Leung H-Y, King TE, et al. (2010) A Predominantly Neolithic Origin for European Paternal Lineages. PLoS Biol 8(1): e1000285. doi:10.1371/journal.pbio.1000285

**Academic Editor:** David Penny, Massey University, New Zealand

**Received** May 8, 2009; **Accepted** December 10, 2009; **Published** January 19, 2010

**Copyright:** © 2010 Balaresque et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** MAJ was supported by a Wellcome Trust Senior Fellowship in Basic Biomedical Science (grant number 057559); PB, GRB, SMA, ZHR, and CTS were supported by the Wellcome Trust (www.wellcome.ac.uk). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

**Abbreviations:** hg, haplogroup; KYA, thousand years ago; mtDNA, mitochondrial DNA; TMRCA, time to the most recent common ancestor.

* E-mail: maj4@le.ac.uk

## Introduction

Events underlying the distribution of genetic diversity among modern European populations have been the subject of intense debate since the first genetic data became available [1]. Anatomically modern humans, originating in East Africa, colonized Europe from the Near East ~40 thousand years ago (KYA), then during the last glacial maximum populations retreated into the peninsulas of Iberia, Italy, and the Balkans, followed by northward recolonization from these refugia ~14 KYA. The most important cultural transition was the adoption of agriculture originating in the Fertile Crescent in the Near East at the start of the Neolithic, ~10 KYA [2]. It spread rapidly westwards via Anatolia [3] (Figure 1A), reaching Ireland by 6 KYA, accompanied by the development of sedentary populations and demographic expansion. Debate has focused on whether this spread was due to the movement and expansion of Near-Eastern farmers (demic diffusion), or to the transmission of cultural innovation to existing populations (acculturation), who then themselves expanded.

The observation of southeast-northwest frequency clines for “classical” genetic markers [1,4], autosomal DNA markers [5,6], and Y-chromosomal markers [7,8] (though not for mitochondrial

DNA [mtDNA] [9]) has been used to support the demic diffusion model. No dates can be automatically attached to these clines, however, and some [1], detected by principal component analysis, may simply reflect isolation by distance [10]. The direction of movement underlying a cline can also be ambiguous: the high-frequency pole could indicate the area of preexisting substrate least affected by a migration originating far away, or the final destination of a wave of migration into thinly populated territory, where expansion and drift have had their greatest effects [11].

The origins of a frequency cline of a lineage can be illuminated by analysing the diversity within it. For Y-chromosomal lineages defined by binary markers (haplogroups), this can be done using multiple microsatellites. This approach has been applied to haplogroups E, J [12], and I [13] within Europe, but the major western European lineage has not yet been focused upon. The frequency of the major western European lineage, haplogroup (hg) R1b1b2, follows a cline from 12% in Eastern Turkey to 85% in Ireland (Figure 1B), and is currently carried by some 110 million European men. Previous studies of lineages approximately equivalent to hgR1b1b2 [7,8] suggested that it has a Paleolithic origin, based simply on its high frequency in the west. Here, in contrast, we show that the geographical distribution of diversity within the haplogroup is best explained by its spread from a single

<figure><img src="image_2.png" /></figure>
## Author Summary

Arguably the most important cultural transition in the history of modern humans was the development of farming, since it heralded the population growth that culminated in our current massive population size. The genetic diversity of modern populations retains the traces of such past events, and can therefore be studied to illuminate the demographic processes involved in past events. Much debate has focused on the origins of agriculture in Europe some 10,000 years ago, and in particular whether its westerly spread from the Near East was driven by farmers themselves migrating, or by the transmission of ideas and technologies to indigenous hunter-gatherers. This study examines the diversity of the paternally inherited Y chromosome, focusing on the commonest lineage in Europe. The distribution of this lineage, the diversity within it, and estimates of its age all suggest that it spread with farming from the Near East. Taken with evidence on the origins of other lineages, this indicates that most European Y chromosomes descend from Near Eastern farmers. In contrast, most maternal lineages descend from hunter-gatherers, suggesting a reproductive advantage for farming males over indigenous hunter-gatherer males during the cultural transition from hunting-gathering to farming.

source from the Near East via Anatolia during the Neolithic. Taken together with the evidence on the origins of many other European haplogroups, this indicates that the great majority of the Y chromosomes of Europeans have their origins in the Neolithic expansion.

## Results

To investigate the origins of hgR1b1b2, we assembled a dataset of 840 chromosomes from this haplogroup with associated nine-locus microsatellite haplotypes (Table 1; Table S1). The diversity of the lineage within each population (measured by mean microsatellite variance) should reflect its age: under a hypothesis of recolonization from southern refugia, we expect a gradient of diversity correlating with latitude, whereas Neolithic expansion from Anatolia predicts a correlation primarily with longitude. Figure 1C shows the geographical distribution of mean microsatellite variance, and Figure 2 shows that although there is no evidence for correlation with latitude (R² = 0.06; *p* = 0.268), the correlation with longitude is significant (R² = 0.358; *p* = 0.004), with greatest diversity in the east (strongly influenced by highly diverse samples within Turkey), thus providing support for the Neolithic colonization hypothesis.

The two hypotheses also make different predictions for the number of sources of diversity within hgR1b1b2: under the postglacial recolonization model, we expect multiple sources, whereas under the Neolithic expansion model, we expect only one. We can test this by examining the phylogenetic relationships among microsatellite haplotypes. A reduced median network of 859 haplotypes (Figure 3) shows a simple star-like structure indicative of expansion from one source: 74 haplotypes (8.6%) lie in its central node, and this node plus its single-step mutational neighbours together comprise 214 haplotypes (24.9%). Haplotypes belonging to populations from all three refugia are present in the core of the network. This pattern seems incompatible with recolonization from differentiated refugial populations, and in terms of the history of hgR1b1b2, the refugia possess no special status. The core of the network also contains haplotypes from

Turkey (Anatolia), which is compatible with a subpopulation from this region acting as a source for the westwards-expanding lineage.

Does the time to the most recent common ancestor (TMRCA) of the hgR1b1b2 chromosomes support a Paleolithic origin? Mean estimates for individual populations vary (Table 2), but the oldest value is in Central Turkey (7,989 y [95% confidence interval (CI): 5,661–11,014]), and the youngest in Cornwall (5,460 y [3,764–7,777]). The mean estimate for the entire dataset is 6,512 y (95% CI: 4,577–9,063 years), with a growth rate of 1.95% (1.02%–3.30%). Thus, we see clear evidence of rapid expansion, which cannot have begun before the Neolithic period.

The similarity between the isochron map of Neolithic sites (Figure 1A; [3]) and those of hgR1b1b2 frequency (Figure 1B) and diversity (Figure 1C) is striking. Further support for the association of the expansion of hgR1b1b2 with that of farming comes from a statistical comparison of the variables. The frequency of hgR1b1b2 at different points in Europe is significantly negatively correlated (R² = 0.390; *p* = 0.0005) with the dates of local Neolithic sites (Figure 4A). For the local variance of the microsatellite haplotypes within hgR1b1b2, the correlation with Neolithic dates is significantly positive (R² = 0.331; *p* = 0.0124; Figure 4B).

## Discussion

Previous observations of the east-west clinal distribution of the common Western European hgR1b1b2 (or its equivalent) [7,8] considered it to be part of a Paleolithic substrate into which farmers from the Near East had diffused. Later analyses have also considered variance, and have conformed to the Paleolithic explanation [14,15]. Here, we concur that the cline results from demic diffusion, but our evidence supports a different interpretation: that R1b1b2 was carried as a rapidly expanding lineage from the Near East via Anatolia to the western fringe of Europe during the Neolithic. Such mutations arising at the front of a wave of expansion have a high probability of surviving and being propagated, and can reach high frequencies far from their source [11]. Successive founder effects at the edge of the expansion wave can lead to a reduction in microsatellite diversity, even as the lineage increases in frequency.

The innovations in the Near East also spread along the southern shore of the Mediterranean, reflected in the expansion of hgE1b1b1b (E-M81) [16], which increases in frequency and reduces in diversity from east to west. In sub-Saharan Africa, hgE1b1a (E-M2) underwent a massive expansion associated with the Bantu expansion [17,18]. In India, the spread of agriculture has been associated with the introduction of several Y lineages [19], and in Japan, lineages within hgO spread with the Yayoi migration [20], which brought wet rice agriculture to the archipelago. On a more recent timescale, the expansion of the Han culture in China has been linked to demic diffusion [21]. In this context, the apparently low contribution of incoming Y chromosomes to the European Neolithic, despite its antiquity and impact, has appeared anomalous. Our interpretation of the history of hgR1b1b2 now makes Europe a prime example of how expansion of a Y-chromosomal lineage tends to accompany technological and cultural change.

Other lineages also show evidence of European Neolithic expansion, hgE1b1b1 (E-M35) and hgJ, in particular [12]. Indeed, hgI is the only major lineage for which a Paleolithic origin is generally accepted, but it comprises only 18% of European Y chromosomes [13]. The Basques contain only 8%–20% of this lineage, but 75%–87% hgR1b1b2 (Table S1); our findings therefore challenge their traditional “Mesolithic relict” status, and in particular, their use as a proxy for a Paleolithic parental population in admixture modelling of European Y-chromosomal prehistory [22].

<figure><img src="image_3.png" /></figure>
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

**Figure 1.** Maps showing dates of the spread of early farming in Europe, and the frequency and microsatellite variance of **haplogroup R1b1b2**. (A) Isochron map representing dates of early Neolithic sites in Europe, based on data of Pinhasi et al. (2005) [3]. KYBP, thousand years before present. (B) Geographical distribution of haplogroup frequency of hgR1b1b2, shown as an interpolated spatial frequency surface. Filled circles indicate populations for which microsatellite data and TMRCA estimates are available. Unfilled circles indicate populations included to illustrate R1b1b2 frequency only. Population codes are defined in Table 1. (C) Geographical distribution of mean microsatellite variance within hgR1b1b2, shown as an interpolated spatial frequency surface. Samples shown are those used for the calculation of variance only.
doi:10.1371/journal.pbio.1000285.g001

<figure><img src="image_7.png" /></figure>
<table><caption>**Table 1.** Frequency of haplogroup R1b1b2 in European populations, with geographical coordinates for sampled populations.</caption><thead><tr><th>Country</th><th>Area of Sampling (Population)</th><th>Abbreviation</th><th>Included in All Analyses?</th><th>Longitude West</th><th>Latitude North</th><th>N<sup>a</sup></th><th>% R1b1b2<sup>b</sup></th><th>Source</th></tr></thead><tbody><tr><td>Bosnia-Herzegovina</td><td>National</td><td>BO</td><td></td><td>17.650</td><td>43.850</td><td>256</td><td>3.9</td><td>[40]</td></tr><tr><td>Denmark</td><td>National</td><td>DK</td><td>Y</td><td>9.654</td><td>54.513</td><td><strong>56</strong></td><td>42.9</td><td>Present study</td></tr><tr><td>England</td><td>Cornwall</td><td>EN1</td><td>Y</td><td>-4.955</td><td>50.442</td><td><strong>64</strong></td><td>78.1</td><td>Present study</td></tr><tr><td>England</td><td>Leicestershire</td><td>EN2</td><td>Y</td><td>-1.130</td><td>52.637</td><td><strong>43</strong></td><td>62.0</td><td>Present study</td></tr><tr><td>France</td><td>Basques</td><td>FR1</td><td>Y</td><td>-1.305</td><td>43.384</td><td><strong>61</strong></td><td>75.4</td><td>Present study</td></tr><tr><td>France</td><td>Baie de Somme</td><td>FR2</td><td>Y</td><td>1.603</td><td>50.237</td><td><strong>43</strong></td><td>62.8</td><td>Present study</td></tr><tr><td>France</td><td>Finistère</td><td>FR3</td><td>Y</td><td>-4.264</td><td>48.233</td><td><strong>75</strong></td><td>76.0</td><td>Present study</td></tr><tr><td>France</td><td>Haute-Garonne</td><td>FR4</td><td>Y</td><td>1.443</td><td>43.604</td><td><strong>57</strong></td><td>78.9</td><td>Present study</td></tr><tr><td>France</td><td>Ile et Vilaine</td><td>FR5</td><td></td><td>-1.605</td><td>48.170</td><td><strong>82</strong></td><td>80.5</td><td>Present study</td></tr><tr><td>France</td><td>Loire-Atlantique</td><td>FR6</td><td></td><td>-1.741</td><td>47.348</td><td><strong>48</strong></td><td>77.1</td><td>Present study</td></tr><tr><td>France</td><td>Vendée</td><td>FR7</td><td>Y</td><td>-1.469</td><td>46.676</td><td><strong>50</strong></td><td>68.0</td><td>Present study</td></tr><tr><td>Germany</td><td>Bavaria</td><td>GE1</td><td>Y</td><td>11.319</td><td>48.985</td><td><strong>80</strong></td><td>32.3</td><td>Present study</td></tr><tr><td>Germany</td><td>National</td><td>GE</td><td>Y</td><td>10.451</td><td>51.165</td><td>1215</td><td>38.9</td><td>[41]<sup>c</sup></td></tr><tr><td>Greece</td><td>National</td><td>GR</td><td></td><td>21.824</td><td>39.074</td><td>171</td><td>13.5</td><td>[42]</td></tr><tr><td>Italy</td><td>North-East (Ladin)</td><td>IT1</td><td>Y</td><td>11.552</td><td>46.528</td><td><strong>79</strong></td><td>60.8</td><td>Present study</td></tr><tr><td>Italy</td><td>North-West</td><td>IT2</td><td>Y</td><td>7.912</td><td>44.875</td><td>99</td><td>45.0</td><td>Present study</td></tr><tr><td>Ireland</td><td>National</td><td>IR</td><td>Y</td><td>-8.244</td><td>53.413</td><td>796</td><td>85.4</td><td>[43]</td></tr><tr><td>Italy</td><td>Sardinia</td><td>IT3</td><td></td><td>8.948</td><td>39.991</td><td>930</td><td>17.0</td><td>[44]</td></tr><tr><td>Netherlands</td><td>National</td><td>NL</td><td>Y</td><td>5.417</td><td>52.246</td><td><strong>84</strong></td><td>42.0</td><td>Present study</td></tr><tr><td>Poland</td><td>National</td><td>PL</td><td></td><td>19.145</td><td>51.919</td><td>913</td><td>11.6</td><td>[41]<sup>c</sup></td></tr><tr><td>Portugal</td><td>South</td><td>PO</td><td></td><td>-8.176</td><td>37.750</td><td><strong>78</strong></td><td>46.2</td><td>Present study</td></tr><tr><td>Russia</td><td>Belgorod</td><td>RU1</td><td></td><td>36.480</td><td>50.780</td><td>143</td><td>2.8</td><td>[45]</td></tr><tr><td>Russia</td><td>Ostrov</td><td>RU2</td><td></td><td>28.320</td><td>57.350</td><td>75</td><td>2.7</td><td>[45]</td></tr><tr><td>Russia</td><td>Pristen</td><td>RU3</td><td></td><td>36.710</td><td>51.230</td><td>45</td><td>2.2</td><td>[45]</td></tr><tr><td>Russia</td><td>Repievka</td><td>RU4</td><td></td><td>38.650</td><td>51.080</td><td>96</td><td>5.2</td><td>[45]</td></tr><tr><td>Russia</td><td>Roslavl</td><td>RU5</td><td></td><td>32.870</td><td>53.950</td><td>107</td><td>11.2</td><td>[45]</td></tr><tr><td>Spain</td><td>Andalucia East</td><td>SP1</td><td></td><td>-3.209</td><td>37.513</td><td><strong>95</strong></td><td>72.0</td><td>Present study</td></tr><tr><td>Spain</td><td>Andalucia West</td><td>SP2</td><td>Y</td><td>-5.17</td><td>36.34</td><td><strong>72</strong></td><td>55.0</td><td>Present study</td></tr><tr><td>Spain</td><td>Basques</td><td>SP3</td><td></td><td>-2.430</td><td>42.580</td><td><strong>116</strong></td><td>87.1</td><td>Present study</td></tr><tr><td>Spain</td><td>Catalonia</td><td>SP4</td><td>Y</td><td>2.460</td><td>41.560</td><td><strong>80</strong></td><td>81.3</td><td>Present study</td></tr><tr><td>Spain</td><td>Castilla La Mancha</td><td>SP5</td><td>Y</td><td>-3.15</td><td>39.41</td><td><strong>63</strong></td><td>72.0</td><td>Present study</td></tr><tr><td>Serbia</td><td>National</td><td>SB</td><td></td><td>20.759</td><td>44.178</td><td><strong>100</strong></td><td>10.0</td><td>Present study</td></tr><tr><td>Spain</td><td>Galicia</td><td>SP6</td><td>Y</td><td>-8.150</td><td>42.510</td><td><strong>88</strong></td><td>58.0</td><td>Present study</td></tr><tr><td>Slovenia</td><td>National</td><td>SL</td><td></td><td>15.366</td><td>45.609</td><td><strong>70</strong></td><td>20.6</td><td>Present study</td></tr><tr><td>Turkey</td><td>Central</td><td>TK1</td><td>Y</td><td>34.036</td><td>38.942</td><td>152</td><td>19.1</td><td>[14]</td></tr><tr><td>Turkey</td><td>East</td><td>TK2</td><td>Y</td><td>40.110</td><td>38.921</td><td>208</td><td>12.0</td><td>[14]</td></tr><tr><td>Turkey</td><td>West</td><td>TK3</td><td>Y</td><td>28.570</td><td>39.243</td><td>163</td><td>13.5</td><td>[14]</td></tr><tr><td>Wales</td><td>National</td><td>WA</td><td></td><td>-3.793</td><td>52.170</td><td><strong>65</strong></td><td>92.3</td><td>Present study</td></tr></tbody></table>

Is the predominance of Neolithic-expansion lineages among Y chromosomes reflected in other parts of the genome? Mitochondrial DNA diversity certainly presents a different picture: no east-west cline is discernible, most lineages have a Paleolithic TMRCA [23], and hgH [24] and hgV [25] show signatures of postglacial expansion from the Iberian peninsula. Demic diffusion involves

both females and males, but the disparity between mtDNA and Y-chromosomal patterns could arise from an increased and transmitted reproductive success for male farmers compared to indigenous hunter-gatherers, without a corresponding difference between females from the two groups. This would lead to the expansion of incoming Y lineages—as suggested by the high

<figure><img src="image_8.png" /></figure>

<footnote><sup>a</sup>Figures in bold indicate samples typed for M269 in this study.</footnote>

<footnote><sup>b</sup>The number of men currently carrying hgR1b1b2 chromosomes (see Introduction) was approximated from these proportions and population census sizes given at http://www populationdata net/europe.php.</footnote>

<footnote><sup>c</sup>Chromosomes considered to belong to hgR1*(xR1a1).</footnote>

<footnote>Y, yes.</footnote>

<footnote>doi:10.1371/journal.pbio.1000285.t001</footnote>
<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 2.** Relationship of diversity among 840 R1b1b2 chromosomes with (A) longitude and (B) latitude. Population codes are defined in Table 1.
doi:10.1371/journal.pbio.1000285.g002</figcaption></figure>

<figure><img src="image_11.png" /></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 3.** Reduced median network of microsatellite haplotypes within haplogroup R1b1b2. Molecular relationships between the nine-locus microsatellite haplotypes of 849 hgR1b1b2 chromosomes, including seven Serbian and two Greek haplotypes not included in the other analyses because population sample sizes were too small. Circles represent haplotypes, with area proportional to frequency and coloured according to population. Lines between circles represent microsatellite mutational steps.
doi:10.1371/journal.pbio.1000285.g003</figcaption></figure>

growth rate observed for hgR1b1b2. Similar conclusions have been reached for the Bantu expansion (in which the current Bantu-speaking populations carry many mtDNA lineages originating from hunter-gatherers [26]), the introduction of agriculture to India [19] and the Han expansion [21].

Some studies have found evidence of east-west clines for autosomal loci [6,27]. By contrast, recent genome-wide SNP typing surveys [28–30] find a basic south-north division or gradient, including greater diversity in the south, but they provide no indication of the time-depth of the underlying events, which could in principle involve contributions from the original colonization, postglacial Paleolithic recolonization, Neolithic expansion, and later contact between Africa and southern Europe [31].

The distinction between the geographical patterns of variation of the Y chromosome and those of mtDNA suggest sex-specific factors in patterning European diversity, but the rest of the genome has yet to reveal definitive information. Detailed studies of X-chromosomal and autosomal haplotypes promise to further illuminate the roles of males and females in prehistory.

# Materials and Methods

## Ethics Statement

Males were recruited with informed consent, following ethical approval by the Leicestershire Research Ethics Committee and the ethics committees of the Universities of Ferrara, Pavia, and Exeter and Plymouth.

## DNA Samples and Haplotyping

A total of 2,574 DNA samples from European males, assigned to populations based on two generations of residence, were typed for the SNP M269 [17], defining hgR1b1b2. Following PCR amplification using the primers 5'-CTAAAGATCAGAGTATCTCCCTTTG-3' and 5'-ATTTCTAGGAGTTCACTGTATTAC-3', the T to C transition was analysed by digestion with BstNI, which cleaves M269-C-allele chromosomes only. Samples from the Iberian peninsula were typed using the SNaPshot (ABI) procedure [31]. Haplotype data were obtained for up to 20 Y-specific microsatellites [32,33]. Data from the Ysearch database (http://www.ysearch.org) for Germany (GE) and Ireland (IR) were added, together with published data for Turkey,

<figure><img src="image_13.png" /></figure>
<table><caption>**Table 2.** Estimates of TMRCA for individual populations, arranged from west to east.</caption><thead><tr><th>Population</th><th>TMRCA/y (mean [95% CI])</th></tr></thead><tbody><tr><td>Ireland</td><td>5,533 (4,094–7,391)</td></tr><tr><td>Spain–Galicia</td><td>6,584 (4,923–8,684)</td></tr><tr><td>Spain–Andalucia West</td><td>6,208 (4,476–8,463)</td></tr><tr><td>England–Cornwall</td><td>5,460 (3,764–7,777)</td></tr><tr><td>France–Finistère</td><td>6,432 (4,786–8,571)</td></tr><tr><td>Spain–Castilla La Mancha</td><td>6,706 (4,772–9,261)</td></tr><tr><td>France–Vendée</td><td>6,787 (4,575–9,853)</td></tr><tr><td>France–Basques</td><td>5,797 (4,133–8,065)</td></tr><tr><td>England–Leicestershire</td><td>5,981 (4,051–8,439)</td></tr><tr><td>France–Haute–Garonne</td><td>5,925 (4,296–8,114)</td></tr><tr><td>France–Baie de Somme</td><td>7,384 (5,259–10,131)</td></tr><tr><td>Spain–Catalonia</td><td>5,800 (4,410–7,544)</td></tr><tr><td>Netherlands</td><td>6,952 (5,051–9,410)</td></tr><tr><td>Italy–North-West</td><td>5,944 (3,718–8,842)</td></tr><tr><td>Denmark</td><td>6,555 (4,391–9,386)</td></tr><tr><td>Germany–National</td><td>6,138 (4,627–7,997)</td></tr><tr><td>Germany–Bavaria</td><td>7,282 (5,059–10,139)</td></tr><tr><td>Italy–North-East (Ladin)</td><td>6,995 (4,635–10,396)</td></tr><tr><td>Turkey–West</td><td>7,304 (5,022–10,359)</td></tr><tr><td>Turkey–Central</td><td>7,989 (5,661–11,014)</td></tr><tr><td>Turkey–East</td><td>7,000 (4,423–10,490)</td></tr></tbody></table>

doi:10.1371/journal.pbio.1000285.t002

subdivided into East, West, and Central subpopulations based on published sampling information [14]. To avoid a bias from very large samples of hgR1b1b2 (GE and IR), these were randomly subsampled to give sample sizes of 75. This allowed a comparison of nine-locus haplotypes (DYS19, DYS388, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, and DYS439) for 849 hgR1b1b2 chromosomes, subdivided into 23 populations. Greek and Serbian samples were too small for population-based analyses, but were included in Network analysis.

## Analysis

Neolithic dates, frequencies of hgR1b1b2, and local microsatellite variances were displayed using Surfer 8.02 (Golden Software) by the gridding method. Latitudes and longitudes were based on sampling centres.

Intrahaplogroup diversity was assessed for populations with hgR1b1b2 sample size ≥15 as the mean of the individual microsatellite variances [34], as has been done elsewhere (e.g., [35]); this measure is highly correlated ($R^2$ = 0.871; $p$ = 6.72×10⁻¹⁰) with a more conventional measure, average squared distance (ASD) [36]. Regression analyses were carried out in the R statistical package [37] to compare these two measures, and also to compare mean of variance with latitude and longitude.

A reduced median network [38] of microsatellite haplotypes was constructed using Network 4.5 and Network Publisher, using weighting based on the inverse of the microsatellite variances.

TMRCA and population growth rates were estimated using BATWING [39], under a model of exponential population growth and splitting. Whereas standard use of BATWING assumes a random sample from a population, we validated its use to analyse single haplogroups. Justification of this, together with other details, is given in Text S1.

To assess the correlation between the dates of Neolithic sites and the local hgR1b1b2 frequency and variance, we considered 765 sites and their associated calibrated radiocarbon dates [3]. We identified sites lying within a buffer-zone of 150-km radius around each location for which we had frequency or variance data (Figure 1B and 1C). When more than one site was identified in a given buffer-zone, we considered the mean of the dates. Regression analyses were carried out as described above.

## Supporting Information

**Table S1 Haplotype data.** Population abbreviations are as in Table 1; for each microsatellite (DYS19–DYS439), repeat unit numbers are given.

Found at: doi:10.1371/journal.pbio.1000285.s001 (0.14 MB PDF)

**Text S1 Details of application of BATWING.**

Found at: doi:10.1371/journal.pbio.1000285.s002 (0.14 MB DOC)

## Acknowledgments

We thank all DNA donors, Agnar Helgason for discussions, Nicolas Poulet for help with spatial analysis, and Jean-Michel Dugoujon and Evelyne Guitard for assistance with sample collection.

## Author Contributions

The author(s) have made the following declarations about their contributions: Conceived and designed the experiments: CTS MAJ. Performed the experiments: PB GB SA HYL ZHR. Analyzed the data: PB IJW MAJ. Contributed reagents/materials/analysis tools: TEK ZHR JG JPM AM AGD GB CP. Wrote the paper: PB CTS MAJ.

<figure><img src="image_14.png" /></figure>
<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 4.** Correlation of dates of Neolithic sites with *hgR1b1b2* (A) frequency and (B) variance. Population codes are defined in Table 1.
YBP, years before present.
doi:10.1371/journal.pbio.1000285.g004</figcaption></figure>

<figure><img src="image_17.png" /></figure>
References

1. Menozzi P, Piazza A, Cavalli-Sforza LL (1978) Synthetic maps of human gene frequencies in Europeans. Science 201: 786-792.

2. Jobling MA, Hurles ME, Tyler-Smith C (2004) Human evolutionary genetics: origins, peoples and disease. New York (New York): Garland Science. 523 p.

3. Pinhasi R, Fort J, Ammerman AJ (2005) Tracing the origin and spread of agriculture in Europe. PLoS Biol 3: e410. doi:10.1371/journal.pbio.00030410.

4. Cavalli-Sforza LL, Menozzi P, Piazza A (1994) The history and geography of human genes. Princeton (New Jersey): Princeton University Press. 518 p.

5. Chikhi L, Destro-Bisol G, Bertorelle G, Pascali V, Barbujani G (1998) Clines of nuclear DNA markers suggest a largely neolithic ancestry of the European gene pool. Proc Natl Acad Sci USA 95: 9053-9058.

6. Belle EM, Landry PA, Barbujani G (2006) Origins and evolution of the Europeans' genome: evidence from multiple microsatellite loci. Proc Biol Sci 273: 1595-1602.

7. Rosser ZH, Zerjal T, Hurles ME, Adojaan M, Alavantic D, et al. (2000) Y-chromosomal diversity within Europe is clinal and influenced primarily by geography, rather than by language. Am J Hum Genet 67: 1526-1543.

8. Semino O, Passarino G, Oefner PJ, Lin AA, Arbuzova S, et al. (2000) The genetic legacy of Paleolithic Homo sapiens sapiens in extant Europeans: a Y chromosome perspective. Science 290: 1155-1159.

9. Richards M, Corte-Real H, Forster P, Macaulay V, Wilkinson-Herbots H, et al. (1996) Paleolithic and neolithic lineages in the European mitochondrial gene pool. Am J Hum Genet 59: 185-203.

10. Novembre J, Stephens M (2008) Interpreting principal component analyses of spatial population genetic variation. Nat Genet 40: 646-649.

11. Edmonds CA, Lillie AS, Cavalli-Sforza LL (2004) Mutations arising in the wave front of an expanding population. Proc Natl Acad Sci USA 101: 975-979.

12. Semino O, Magri C, Benuzzi G, Lin AA, Al-Zahery N, et al. (2004) 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 74: 1023-1034.

13. Rootsi S, Magri C, Kivisild T, Benuzzi G, Help H, et al. (2004) Phylogeography of Y-chromosome haplogroup I reveals distinct domains of prehistoric gene flow in Europe. Am J Hum Genet 75: 128-137.

14. Cinnioglu C, King R, Kivisild T, Kalfoglu E, Atasoy S, et al. (2004) Excavating Y-chromosome haplotype strata in Anatolia. Hum Genet 114: 127-148.

15. Pericic M, Lauc LB, Klaric IM, Rootsi S, Janicijevic B, et al. (2005) High-resolution phylogenetic analysis of southeastern Europe traces major episodes of paternal gene flow among Slavic populations. Mol Biol Evol 22: 1964-1975.

16. Arredi B, Poloni ES, Paracchini S, Zerjal T, Fathallah DM, et al. (2004) A predominantly neolithic origin for Y-chromosomal DNA variation in North Africa. Am J Hum Genet 75: 338-345.

17. Cruciani F, Santolamazza P, Shen PD, Macaulay V, Moral P, et al. (2002) A back migration from Asia to sub-Saharan Africa is supported by high-resolution analysis of human Y-chromosome haplotypes. Am J Hum Genet 70: 1197-1214.

18. Beleza S, Gusmao L, Amorim A, Carracedo A, Salas A (2005) The genetic legacy of western Bantu migrations. Hum Genet 117: 366-375.

19. Cordaux R, Deepa E, Vishwanathan H, Stoneking M (2004) Genetic evidence for the demic diffusion of agriculture to India. Science 304: 1125.

20. Hammer MF, Karafet TM, Park H, Omoto K, Harihara S, et al. (2006) Dual origins of the Japanese: common ground for hunter-gatherer and farmer Y chromosomes. J Hum Genet 51: 47-58.

21. Wen B, Li H, Lu D, Song X, Zhang F, et al. (2004) Genetic evidence supports demic diffusion of Han culture. Nature 431: 302-305.

22. Chikhi L, Nichols RA, Barbujani G, Beaumont MA (2002) Y genetic data support the Neolithic demic diffusion model. Proc Natl Acad Sci USA 99: 11008-11013.

23. Richards M, Macaulay V, Hickey E, Vega E, Sykes B, et al. (2000) Tracing European founder lineages in the near eastern mtDNA pool. Am J Hum Genet 67: 1251-1276.

24. Achilli A, Rengo C, Magri C, Battaglia V, Olivieri A, et al. (2004) The molecular dissection of mtDNA haplogroup H confirms that the Franco-Cantabrian glacial

refuge was a major source for the European gene pool. Am J Hum Genet 75: 910-918.

25. Torroni A, Bandelt HJ, Macaulay V, Richards M, Cruciani F, et al. (2001) A signal, from human mtDNA, of postglacial recolonization in Europe. Am J Hum Genet 69: 844-852.

26. Wood ET, Stover DA, Ehret C, Destro-Bisol G, Spedini G, et al. (2005) Contrasting patterns of Y chromosome and mtDNA variation in Africa: evidence for sex-biased demographic processes. Eur J Hum Genet 13: 867-876.

27. Sokal RR, Oden NL, Wilson C (1991) Genetic evidence for the spread of agriculture in Europe by demic diffusion. Nature 351: 143-145.

28. Lao O, Lu TT, Nothnagel M, Junge O, Freitag-Wolf S, et al. (2008) Correlation between genetic and geographic structure in Europe. Curr Biol 18: 1241-1248.

29. Novembre J, Johnson T, Bryce K, Kutalik Z, Boyko AR, et al. (2008) Genes mirror geography within Europe. Nature 456: 98-101.

30. Auton A, Bryce K, Boyko A, Lohmueller K, Novembre J, et al. (2009) Global distribution of genomic diversity underscores rich complex history of continental human populations. Genome Res 19: 795-803.

31. Adams SM, Bosch E, Balaresque PL, Lee AC, Arroyo E, et al. (2008) The genetic legacy of religious diversity and intolerance: paternal lineages of Christians, Jews and Muslims in the Iberian Peninsula. Am J Hum Genet 83: 725-736.

32. Bosch E, Lee AC, Calafell F, Arroyo E, Henneman P, et al. (2002) High resolution Y chromosome typing: 19 STRs amplified in three multiplex reactions. Forensic Sci Int 125: 42-51.

33. Parkin EJ, Kraayenbrink T, van Driem GL, Tshering K, de Knijff P, et al. (2006) 26-locus Y-STR typing in a Bhutanese population sample. Forensic Sci Int 161: 1-7.

34. Kayser M, Krawczak M, Excoffier L, Dieltjes P, Corach D, et al. (2001) An extensive analysis of Y-chromosomal microsatellite haplotypes in globally dispersed human populations. Am J Hum Genet 68: 990-1018.

35. Sengupta S, Zhivotovsky LA, King R, Mehdi SQ, Edmonds CA, et al. (2006) 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 Hum Genet 78: 202-221.

36. Goldstein DB, Linares AR, Cavalli-Sforza LL, Feldman MW (1995) An evaluation of genetic distances for use with microsatellite loci. Genetics 139: 463-471.

37. R Development Core Team (2007) R: a language and environment for statistical computing. Vienna (Austria): R Foundation for Statistical Computing.

38. Bandelt H-J, Forster P, Röhl A (1999) Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16: 37-48.

39. Wilson IJ, Weale ME, Balding DJ (2003) Inferences from DNA data: population histories, evolutionary processes and forensic match probabilities. J Roy Statist Soc A 166: 1-33.

40. Marjanovic D, Fornarino S, Montagna S, Primorac D, Hadziselimovic R, et al. (2005) The peopling of modern Bosnia-Herzegovina: Y-chromosome haplogroups in the three main ethnic groups. Ann Hum Genet 69: 757-763.

41. Kayser M, Lao O, Anslinger K, Augustin C, Bargel G, et al. (2005) Significant genetic differentiation between Poland and Germany follows present-day political borders, as revealed by Y-chromosome analysis. Hum Genet 117: 428-443.

42. King RJ, Ozcan SS, Carter T, Kalfoglu E, Atasoy S, et al. (2008) Differential Y-chromosome Anatolian influences on the Greek and Cretan Neolithic. Ann Hum Genet 72: 205-214.

43. Moore LT, McEvoy B, Cape E, Simms K, Bradley DG (2006) A Y-chromosome signature of hegemony in Gaelic Ireland. Am J Hum Genet 78: 334-338.

44. Contu D, Morelli L, Santoni F, Foster JW, Francalacci P, et al. (2008) Y-chromosome based evidence for pre-neolithic origin of the genetically homogeneous but diverse Sardinian population: inference for association scans. PLoS ONE 3: e1430. doi:10.1371/journal.pone.0001430.

45. Balanovsky O, Rootsi S, Pshenichnov A, Kivisild T, Churnosov M, et al. (2008) Two sources of the Russian patrilineal heritage in their Eurasian context. Am J Hum Genet 82: 236-250.

<figure><img src="image_18.png" /></figure>