Guido Barbujani

Y-Chromosomal Diversity in Europe Is Clinal and Influenced Primarily by Geography, Rather than by Language

Clinal patterns of autosomal genetic diversity within Europe have been interpreted in previous studies in terms of a Neolithic demic diffusion model for the spread of agriculture; in contrast, studies using mtDNA have traced many founding lineages to the Paleolithic and have not shown strongly clinal variation. We have used 11 human Y-chromosomal biallelic polymorphisms, defining 10 haplogroups, to analyze a sample of 3,616 Y chromosomes belonging to 47 European and circum-European populations. Patterns of geographic differentiation are highly nonrandom, and, when they are assessed using spatial autocorrelation analysis, they show significant clines for five of six haplogroups analyzed. Clines for two haplogroups, representing 45% of the chromosomes, are continentwide and consistent with the demic diffusion hypothesis. Clines for three other haplogroups each have different foci and are more regionally restricted and are likely to reflect distinct population movements, including one from north of the Black Sea. Principal-components analysis suggests that populations are related primarily on the basis of geography, rather than on the basis of linguistic affinity. This is confirmed in Mantel tests, which show a strong and highly significant partial correlation between genetics and geography but a low, nonsignificant partial correlation between genetics and language. Genetic-barrier analysis also indicates the primacy of geography in the shaping of patterns of variation. These patterns retain a strong signal of expansion from the Near East but also suggest that the demographic history of Europe has been complex and influenced by other major population movements, as well as by linguistic and geographic heterogeneities and the effects of drift.

Y genetic data support the Neolithic demic diffusion model

There still is no general agreement on the origins of the European gene pool, even though Europe has been more thoroughly investigated than any other continent. In particular, there is continuing controversy about the relative contributions of European Palaeolithic hunter-gatherers and of migrant Near Eastern Neolithic farmers, who brought agriculture to Europe. Here, we apply a statistical framework that we have developed to obtain direct estimates of the contribution of these two groups at the time they met. We analyze a large dataset of 22 binary markers from the non-recombining region of the Y chromosome (NRY), by using a genealogical likelihood-based approach. The results reveal a significantly larger genetic contribution from Neolithic farmers than did previous indirect approaches based on the distribution of haplotypes selected by using post hoc criteria. We detect a significant decrease in admixture across the entire range between the Near East and Western Europe. We also argue that local hunter-gatherers contributed less than 30% in the original settlements. This finding leads us to reject a predominantly cultural transmission of agriculture. Instead, we argue that the demic diffusion model introduced by Ammerman and Cavalli-Sforza [Ammerman, A. J. & Cavalli-Sforza, L. L. (1984) The Neolithic Transition and the Genetics of Populations in Europe (Princeton Univ. Press, Princeton)] captures the major features of this dramatic episode in European prehistory.

Evidence for a genetic discontinuity between Neandertals and 24,000-year-old anatomically modern Europeans.

During the late Pleistocene, early anatomically modern humans coexisted in Europe with the anatomically archaic Neandertals for some thousand years. Under the recent variants of the multiregional model of human evolution, modern and archaic forms were different but related populations within a single evolving species, and both have contributed to the gene pool of current humans. Conversely, the Out-of-Africa model considers the transition between Neandertals and anatomically modern humans as the result of a demographic replacement, and hence it predicts a genetic discontinuity between them. Following the most stringent current standards for validation of ancient DNA sequences, we typed the mtDNA hypervariable region I of two anatomically modern Homo sapiens sapiens individuals of the Cro-Magnon type dated at about 23 and 25 thousand years ago. Here we show that the mtDNAs of these individuals fall well within the range of variation of todays humans, but differ sharply from the available sequences of the chronologically closer Neandertals. This discontinuity is difficult to reconcile with the hypothesis that both Neandertals and early anatomically modern humans contributed to the current European gene pool.

Geographic patterns of mtDNA diversity in Europe.

Genetic diversity in Europe has been interpreted as a reflection of phenomena occurring during the Paleolithic ( approximately 45,000 years before the present [BP]), Mesolithic ( approximately 18,000 years BP), and Neolithic ( approximately 10,000 years BP) periods. A crucial role of the Neolithic demographic transition is supported by the analysis of most nuclear loci, but the interpretation of mtDNA evidence is controversial. More than 2,600 sequences of the first hypervariable mitochondrial control region were analyzed for geographic patterns in samples from Europe, the Near East, and the Caucasus. Two autocorrelation statistics were used, one based on allele-frequency differences between samples and the other based on both sequence and frequency differences between alleles. In the global analysis, limited geographic patterning was observed, which could largely be attributed to a marked difference between the Saami and all other populations. The distribution of the zones of highest mitochondrial variation (genetic boundaries) confirmed that the Saami are sharply differentiated from an otherwise rather homogeneous set of European samples. However, an area of significant clinal variation was identified around the Mediterranean Sea (and not in the north), even though the differences between northern and southern populations were insignificant. Both a Paleolithic expansion and the Neolithic demic diffusion of farmers could have determined a longitudinal cline of mtDNA diversity. However, additional phenomena must be considered in both models, to account both for the north-south differences and for the greater geographic scope of clinical patterns at nuclear loci. Conversely, two predicted consequences of models of Mesolithic reexpansion from glacial refugia were not observed in the present study.

A Predominantly Neolithic Origin for European Paternal Lineages

Most present-day European men inherited their Y chromosomes from the farmers who spread from the Near East 10,000 years ago, rather than from the hunter-gatherers of the Paleolithic.

A latitudinal cline in a Drosophila clock gene

The clock gene period determines biological rhythmicity in Drosophila melanogaster and encodes a protein characterized by an alternating series of threonine-glycine pairs. The minisatellite region encoding the threonine-glycine repeat is polymorphic in length in natural Drosophila melanogaster populations. In this paper we report the geographical analysis of this polymorphism within Europe and North Africa. A robust clinal pattern is observed along a north-south axis. We suggest the possibility that the length polymorphism could be maintained by thermal selection because the threonine-glycine region has been shown to provide thermostability to the circadian phenotype.

DETECTING REGIONS OF ABRUPT CHANGE IN MAPS OF BIOLOGICAL VARIABLES

The areas where the rates of change of biological variables across space are partic- ularly high may correspond to either steep ecological gradients or regions of limited admixture among demes. A method for detecting such biological boundaries was proposed by Womble (1951), who suggested averaging the absolute values of the derivatives of the functions describing biological variation in space at various locations. We present here algorithms that quantify both the mean magnitude and the mean direction of change in surfaces representing distributions of biological measures (such as gene frequencies, measures of quantitative traits, etc.). Infer- ences on the microevolutionary processes affecting the populations can be made by comparing the boundaries detected with the distribution of environmental characteristics, or with the location of factors that may have prevented population admixture. Examples of the application of this method to both simulated data and gene frequencies of two natural populations are given. (Boundary detection; gene frequencies; Wombling)

Y chromosomal haplogroup J as a signature of the post-neolithic colonization of Europe

In order to attain a finer reconstruction of the peopling of southern and central-eastern Europe from the Levant, we determined the frequencies of eight lineages internal to the Y chromosomal haplogroup J, defined by biallelic markers, in 22 population samples obtained with a fine-grained sampling scheme. Our results partially resolve a major multifurcation of lineages within the haplogroup. Analyses of molecular variance show that the area covered by haplogroup J dispersal is characterized by a significant degree of molecular radiation for unique event polymorphisms within the haplogroup, with a higher incidence of the most derived sub-haplogroups on the northern Mediterranean coast, from Turkey westward; here, J diversity is not simply a subset of that present in the area in which this haplogroup first originated. Dating estimates, based on simple tandem repeat loci (STR) diversity within each lineage, confirmed the presence of a major population structuring at the time of spread of haplogroup J in Europe and a punctuation in the peopling of this continent in the post-Neolithic, compatible with the expansion of the Greek world. We also present here, for the first time, a novel method for comparative dating of lineages, free of assumptions of STR mutation rates.

Evidence for Paleolithic and Neolithic gene flow in Europe.

We thank Hans Bandelt, Peter Forster, Luca Cavalli-Sforza, and Eric Minch for giving us access to their unpublished letters and for discussing them with us. We also thank Italo Barrai for fruitful discussion. Of course, these individuals do not share all the views expressed here.

A recent shift from polygyny to monogamy in humans is suggested by the analysis of worldwide Y-chromosome diversity.

Molecular genetic data contain information on the history of populations. Evidence of prehistoric demographic expansions has been detected in the mitochondrial diversity of most human populations and in a Y-chromosome STR analysis, but not in a previous study of 11 Y-chromosome SNPs in Europeans. In this paper, we show that mismatch distributions and tests of mutation/drift equilibrium based on up to 166 Y-chromosome SNPs, in 46 samples from all continents, also fail to support an increase of the male effective population size. Computer simulations show that the low nuclear versus mitochondrial mutation rates cannot explain these results. However, ascertainment bias, i.e., when only highly variable SNP sites are typed, may be concealing any Y SNPs evidence for a recent, but not an ancient, increase in male effective population sizes. The results of our SNP analyses can be reconciled with the expansion of male effective population sizes inferred from STR loci, and with mitochondrial evidence, by admitting that humans were essentially polygynous during much of their history. As a consequence, until recently only a few men may have contributed a large fraction of the Y-chromosome pool at every generation. The number of breeding males may have increased, and the variance of their reproductive success may have decreased, through a recent shift from polygyny to monogamy, which is supported by ethnological data and possibly accompanied the shift from mobile to sedentary communities.