Lounès Chikhi Evolution et Diversité Biologique CNRS Université Paul Sabatier, Toulouse

1. Lounès Chikhi Evolution et Diversité Biologique CNRS Université Paul Sabatier, Toulouse The Neolithic transition in Europe: different views from population genetic (a tentative discussion around some methodological questions)

2. Inference in population genetics • Data collection • Genetic typing • Description of patterns of genetic variability • Analysis and interpretation • Test (simulations)

3. Inference in population genetics • Sampling of “populations” • “Choice” of the markers (genome sampling) • mitochondrial DNA : female demography • Y Chromosome : male demography • nuclear genes (markers: allozymes, microsatellites, RFLP, AFLP, SNPs, etc.) • Description of the patterns : • Diversity within samples • Diversity between samples • Are there spatial patterns ?

4. A similarpattern with Y chromosome data Semino et al. (2000) Science What to do of the patterns ? How to interpret them ?

5. Inference in population genetics • Are the patterns, if any, compatible with hypotheses or demographic scenarios from other areas (archaeology, linguistics, etc.) ?

6. A possible scheme of population movements since Paleolithic 18,000 BP 45,000 BP 10,000 BP

7. Inference in population genetics • Is there a link between these images (archeo-genetico-linguistic) ? • Can we estimate demographic parameters ? • Population : stable ? growing ? bottleneck ? • Admixture between populations ? • Can we date these events ? • Can we detect selection ?

8. nA = 7 r = 8 nA/r = 0.88 He = 0.74 gap Allele frequency nA = 4 r = 7 nA/r = 0.57 He = 0.71 gaps Allele frequency Allele sizes (nb of repetitions) Effect of population size changes on some measures of genetic diversity • nA drops quicker than He because rare alleles are eliminated and do not contribute to • He = 1-Σpi2. • gappy allelic size distributions • range varies little (r=range) Bottleneck

9. Inference in population genetics • Thus there is some information in genetic data about ancient demographic events. • However, this information, may be qualitative rather than quantitative and does not allow us to determine whether other scenarios could have played a role (or selection).

10. Recent data from the Y chromosome have been interpreted as indicating a Neolithic contribution of 22% (Semino et al., 2000). This figure (22%) is the sum of the frequencies of 4 haplotypes called Eu4, 9, 10 and 11 Question : why should the proportion of haplotypes exhibiting a clinal distribution today represent the so-called “Neolithic” contribution?

11. There are two problems with this “estimation”: 1. Clines are only expected for alleles that were present in different frequencies in the populations when they mixed (dilution problem). Moreover, drift in the last 4000-8000 years may have blurred clines that were visible at the time. Many haplotypes are observed only 1, 2 or 3 times in each sample (i.e. no cline is going to be as visible by eye as those observed for the 4 selected haplotypes) 2. Even if it were estimated properly it would be meaningless for understanding the processes of European colonization. A single number cannot summarise a cline.

12. PN = proportion of farmers in any admixed population n= number of admixture events Geometric decrease of Neolithic contrib. from PN to PNn Ex: PN=0.9 and n=25, then PNn=0.07 Average = 100(PN +PN2+…+ PNn)/n. Same value of PN = 0.9 (90% farmers + 10% hunter-gatherers). Horizontal lines are averages: n=10: average = 62% n=25: average = 36% n=50: average = 21% Thus, a lack of pattern or a low average can correspond to a high PN value. +: n=10 Δ: n=25 O: n=50

13. Demic diffusion Cultural diffusion Average Large Genetic contribution Small Genetic contribution Two major models have been proposed (or at least structure current debate) The demic diffusion model: significant correlations between archaeological and genetic maps are explained by a movement of people entering Europe from the Levant and Anatolia during the Neolithic.We would expect a significant genetic contribution. The cultural diffusion model: the spread of agriculture in Europe involved the movement of ideas, not of people.The genetic contribution of Near East farmers to the European gene pool should be limited.

14. Inférence en génétique des populations What kind of inference ? • Qualitative versus quantitative ? • Detection versusestimation. • Models and underlying assumptions.

15. T/N1 T/Nh T/N2 Admixturemodel p1 1 – p1 Parent 1 (hunter-gatherers) Parent 2 (farmers.) Hybrid Past T Parent 1 (Basques) Hybrid (Europe) Parent 2 (Near East) Present Separates the effects of drift and admixture

16. The effect of drift is that the « hybrid » population may not even be intermediate after a limited number of generations. • In other words: • the information on admixture decreases with time. • (ii) It is risky to analyse single locus data when demographic events are ancient.

17. p1 = 0.3 Little drift 1) We simulate data according to the model (figure above) varying some parameters (here drift) 2) The outputs are given to the program implementing the method 3) One distribution is obtained for each simulation More drift

18. 1) We note that for the VERY SAMEscenario inference can be extremely different ! 2) This inference varies from one locus to the other. 3) When two loci produce different estimates, we cannot conclude that they had a different demographic history. 4) Worse : we are in an optimalsituation : we « know » the real p and the data were simulated according to the model. This NEVER happens in real life. p1 = 0.3

19. 1) One solution : multi-locus data. 2) Increasing sample sizes is NOT very useful. 3) Better to have multilocus data much later than one locus juste after: Ex: 5 loci after 100 generations versus 1 locus after 1 generation (for N=1000) 4) Don’t throw your allozyme data away.

20. What if we re-analyse Semino et al.’s data ?

21. Y chromosome data (Semino et al., 2000). p1 represents the hunter-gatherer contribution (descendants = Basques) Each curve corresponds to the analysis of a European population. Significant cline observed for (1-p1) values (i.e. Near Eastern contribution) against geog. distance calculated from the Near East.

22. After Semino et al., 2000 Science

23. As a test we can analyse the same data considering Sardinians as descendants of the hunter-gatherers. We find an extremely similar result. The « Neolithic » contribution is even slightly superior: on the order of 65% instead of 50%.

24. Model-based results (i) (mostly on Y chromosome data) • 1) There are significant clines for the parameter representing the Neolithic contribution Néolithique across Europe. • 2) This “trend” is signifcantly different from that “obtained” by Semino et al. (2000). • 3) The Neolithic contribution appears to be around 50% rather than 22%. • Re-analysis of all European populations using the Sardinian population as P1 shows very similar results with higher Neolithic contribution (average of 65%). • Conclusion: • The cultural diffusion model is unlikely to explain the patterns observed using the Y chrom. data.

25. Tests performed are partial and the model is simplistic but it is a first step towards quantification of demographic parameters clearly identified. • Qualitative approach : • Easy and useful BUT little or misleadingly precise • Quantitative approach: • Assumptions are explicit • Results can be precise (or not) BUT often complicated to interpret and (maybe) model-dependent.

26. Inference in population genetics • Data collection • Genetic typing • Description of patterns of genetic variability • Analysis and interpretation • Test (simulations)

27. Inference in population genetics In case I was not specific enough : Beware the use of any method whose assumptions you do not understand or which have not been extensively tested on simulations : • Nested Clade Analysis • Median-network Thank you AND MANY THANKS TO Mark Beaumont, Mike Bruford, Guido Barbujani, Richard Nichols, etc.