A western route of prehistoric human migration from Africa into the Iberian Peninsula.

Being at the western fringe of Europe, Iberia had a peculiar prehistory and a complex pattern of Neolithization. A few studies, all based on modern populations, reported the presence of DNA of likely African origin in this region, generally concluding it was the result of recent gene flow, probably during the Islamic period. Here, we provide evidence of much older gene flow from Africa to Iberia by sequencing whole genomes from four human remains from northern Portugal and southern Spain dated around 4000 years BP (from the Middle Neolithic to the Bronze Age). We found one of them to carry an unequivocal sub-Saharan mitogenome of most probably West or West-Central African origin, to our knowledge never reported before in prehistoric remains outside Africa. Our analyses of ancient nuclear genomes show small but significant levels of sub-Saharan African affinity in several ancient Iberian samples, which indicates that what we detected was not an occasional individual phenomenon, but an admixture event recognizable at the population level. We interpret this result as evidence of an early migration process from Africa into the Iberian Peninsula through a western route, possibly across the Strait of Gibraltar.

Demographic history and adaptation account for clock gene diversity in humans.

Circadian clocks give rise to daily oscillations in behavior and physiological functions that often anticipate upcoming environmental changes generated by the Earth rotation. In model organisms a relationship exists between several genes affecting the circadian rhythms and latitude. We investigated the allele distributions at 116 000 single-nucleotide polymorphisms (SNPs) of 25 human clock and clock-related genes from the 1000Genomes Project, and at a reference data set of putatively neutral polymorphisms. The global genetic structure at the clock genes did not differ from that observed at the reference data set. We then tested for evidence of local adaptation searching for FST outliers under both an island and a hierarchical model, and for significant association between allele frequencies and environmental variables by a Bayesian approach. A total of 230 SNPs in 23 genes, or 84 SNPs in 19 genes, depending on the significance thresholds chosen, showed signs of local adaptation, whereas a maximum of 190 SNPs in 23 genes had significant covariance with one or more environmental variables. Only two SNPs from two genes (NPAS2 and AANAT) exhibit both elevated population differentiation and covariance with at least one environmental variable. We then checked whether the SNPs emerging from these analyses fall within a set of candidate SNPs associated with different chronotypes or sleep disorders. Correlation of five such SNPs with environmental variables supports a selective role of latitude or photoperiod, but certainly not a major one.

Nine things to remember about human genome diversity.

Understanding how and why humans are biologically different is indispensable to get oriented in the ever-growing body of genomic data. Here we discuss the evidence based on which we can confidently state that humans are the least genetically variable primate, both when individuals and when populations are compared, and that each individual genome can be regarded as a mosaic of fragments of different origins. Each population is somewhat different from any other population, and there are geographical patterns in that variation. These patterns clearly indicate an African origin for our species, and keep a record of the main demographic changes accompanying the peopling of the whole planet. However, only a minimal fraction of alleles, and a small fraction of combinations of alleles along the chromosome, is restricted to a single geographical region (and even less so to a single population), and diversity between members of the same population is very large. The small genomic differences between populations and the extensive allele sharing across continents explain why historical attempts to identify, once and for good, major biological groups in humans have always failed. Nevertheless, racial categorization is all but gone, especially in clinical studies. We argue that racial labels may not only obscure important differences between patients but also that they have become positively useless now that cheap and reliable methods for genotyping are making it possible to pursue the development of truly personalized medicine.

Comparing models on the genealogical relationships among Neandertal, Cro-Magnoid and modern Europeans by serial coalescent simulations.

Populations of anatomically archaic (Neandertal) and early modern (Cro-Magnoid) humans are jointly documented in the European fossil record, in the period between 40 000 and 25 000 years BP, but the large differences between their cultures, morphologies and DNAs suggest that the two groups were not close relatives. However, it is still unclear whether any genealogical continuity between them can be ruled out. Here, we simulated a broad range of demographic scenarios by means of a serial coalescence algorithm in which Neandertals, Cro-Magnoids and modern Europeans were either part of the same mitochondrial genealogy or of two separate genealogies. Mutation rates, population sizes, population structure and demographic growth rates varied across simulations. All models in which anatomically modern (that is, Cro-Magnoid and current) Europeans belong to a distinct genealogy performed better than any model in which the three groups were assigned to the same mitochondrial genealogy. The maximum admissible level of gene flow between Neandertals and the ancestors of current Europeans is 0.001% per generation, one order of magnitude lower than estimated in previous studies not considering genetic data on Cro-Magnoid people.

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.

Genetic characterization of the body attributed to the evangelist Luke.

Historical sources indicate that the evangelist Luke was born in Syria, died in Greece, and then his body was transferred to Constantinople, and from there to Padua, Italy. To understand whether there is any biological evidence supporting a Syrian origin of the Padua body traditionally attributed to Luke, or a replacement in Greece or Turkey, the mtDNA was extracted from two teeth and its control region was cloned and typed. The sequence determined in multiple clones is an uncommon variant of a set of alleles that are common in the Mediterranean region. We also collected and typed modern samples from Syria and Greece. By comparison with these population samples, and with samples from Anatolia that were already available in the literature, we could reject the hypothesis that the body belonged to a Greek, rather than a Syrian, individual. However, the probability of an origin in the area of modern Turkey was only insignificantly lower than the probability of a Syrian origin. The genetic evidence is therefore compatible with the possibility that the body comes from Syria, but also with its replacement in Constantinople.

A multistep process for the dispersal of a Y chromosomal lineage in the Mediterranean area.

In this work we focus on a microsatellite-defined Y-chromosomal lineage (network 1.2) identified by us and reported in previous studies, whose geographic distribution and antiquity appear to be compatible with the Neolithic spread of farmers. Here, we set network 1.2 in the Y-chromosomal phylogenetic tree, date it with respect to other lineages associated with the same movements by other authors, examine its diversity by means of tri- and tetranucleotide loci and discuss the implications in reconstructing the spread of this group of chromosomes in the Mediterranean area. Our results define a tripartite phylogeny within HG 9 (Rosser et al. 2000), with the deepest branching defined by alleles T (Haplogroup Eu10) or G (Haplogroup Eu9) at M172 (Semino et al. 2000), and a subsequent branching within Eu9 defined by network 1.2. Population distributions of HG 9 and network 1.2 show that their occurrence in the surveyed area is not due to the spread of people from a single parental population but, rather, to a process punctuated by at least two phases. Our data identify the wide area of the Balkans, Aegean and Anatolia as the possible homeland harbouring the largest variation within network 1.2. The use of recently proposed tests based on the stepwise mutation model suggests that its spread was associated to a population expansion, with a high rate of male gene flow in the Turkish-Greek area.

DNA diversity and population admixture in Anatolia.

The Turkic language was introduced in Anatolia at the start of this millennium, by nomadic Turkmen groups from Central Asia. Whether that cultural transition also had significant population-genetics consequences is not fully understood. Three nuclear microsatellite loci, the hypervariable region I of the mitochondrial genome, six microsatellite loci of the Y chromosome, and one Alu insertion (YAP) were amplified and typed in 118 individuals from four populations of Anatolia. For each locus, the number of chromosomes considered varied between 51-200. Genetic variation was large within samples, and much less so between them. The contribution of Central Asian genes to the current Anatolian gene pool was quantified using three different methods, considering for comparison populations of Mediterranean Europe, and Turkic-speaking populations of Central Asia. The most reliable estimates suggest roughly 30% Central Asian admixture for both mitochondrial and Y-chromosome loci. That (admittedly approximate) figure is compatible both with a substantial immigration accompanying the arrival of the Turkmen armies (which is not historically documented), and with continuous gene flow from Asia into Anatolia, at a rate of 1% for 40 generations. Because a military invasion is expected to more deeply affect the male gene pool, similar estimates of admixture for female- and male-transmitted traits are easier to reconcile with continuous migratory contacts between Anatolia and its Asian neighbors, perhaps facilitated by the disappearance of a linguistic barrier between them.

Y-chromosome mismatch distributions in Europe.

Ancient demographic events can be inferred from the distribution of pairwise sequence differences (or mismatches) among individuals. We analyzed a database of 3,677 Y chromosomes typed for 11 biallelic markers in 48 human populations from Europe and the Mediterranean area. Contrary to what is observed in the analysis of mitochondrial polymorphisms, Tajima's test was insignificant for most Y-chromosome samples, and in 47 populations the mismatch distributions had multiple peaks. Taken at face value, these results would suggest either (1) that the size of the male population stayed essentially constant over time, while the female population size increased, or (2) that different selective regimes have shaped mitochondrial and Y-chromosome diversity, leading to an excess of rare alleles only in the mitochondrial genome. An alternative explanation would be that the 11 variable sites of the Y chromosome do not provide sufficient statistical power, so a comparison with mitochondrial data (where more than 200 variable sites are studied in Europe) is impossible at present. To discriminate between these possibilities, we repeatedly analyzed a European mitochondrial database, each time considering only 11 variable sites, and we estimated mismatch distributions in stable and growing populations, generated by simulating coalescent processes. Along with theoretical considerations, these tests suggest that the difference between the mismatch distributions inferred from mitochondrial and Y-chromosome data are not a statistical artifact. Therefore, the observed mismatch distributions appear to reflect different underlying demographic histories and/or selective pressures for maternally and paternally transmitted loci.

Genetics and the population history of Europe.

Analysis of genetic variation among modern individuals is providing insight into prehistoric events. Comparisons of levels and patterns of genetic diversity with the predictions of models based on archeological evidence suggest that the spread of early farmers from the Levant was probably the main episode in the European population history, but that both older and more recent processes have left recognizable traces in the current gene pool.

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.

Mitochondrial DNA sequences in prehistoric human remains from the Alps.

The spread of agriculture that started in the Near East about 10 000 years ago caused a dramatic change in the European archaeological record. It is still unclear if that change was caused mostly by movement of people or by cultural transformations. In particular, there is disagreement on what proportion of the current European gene pool is derived either from the pre-agricultural, paleolithic and mesolithic people, or from neolithic farmers immigrating from the south-east. To begin to characterise the mtDNA gene pool of prehistoric Europe we examined five human remains from the Eastern Italian Alps, dated between 14 000 and 3000 years ago. Three of them yielded sufficient amount of mtDNA for analysis. DNA extracts were prepared in two independent laboratories, and PCR products from the first hypervariable segment of the mtDNA control region were cloned and sequenced. Together with the 5200 year old 'ice man', these DNA sequences show that European mtDNA diversity was already high at the beginning of the neolithic period. All the neolithic sequences have been observed in contemporary Europeans, suggesting genealogical continuity between the neolithic and present-day European mtDNA gene pool. The mtDNA sequence from a 14 000 year-old specimen was not observed in any contemporary Europeans, raising the possibility of a lack of continuity between the mesolithic and present-day European gene pools.

High carrier frequency of the 35delG deafness mutation in European populations. Genetic Analysis Consortium of GJB2 35delG.

Congenital deafness accounts for about 1 in 1000 infants and approximately 80% of cases are inherited as an autosomal recessive trait. Recently, it has been demonstrated that connexin 26 (GJB2) gene is a major gene for congenital sensorineural deafness. A single mutation (named 35delG) was found in most recessive families and sporadic cases of congenital deafness, among Caucasoids, with relative frequencies ranging from 28% to 63%. We present here the analysis of the 35delG mutation in 3270 random controls from 17 European countries. We have detected a carrier frequency for 35delG of 1 in 35 in southern Europe and 1 in 79 in central and northern Europe. In addition, 35delG was detected in five out of 376 Jewish subjects of different origin, but was absent in other non-European populations. The study suggests either a single origin for 35delG somewhere in Europe or in the Middle East, and the possible presence of a carrier advantage together with a founder effect. The 35delG carrier frequency of 1 in 51 in the overall European population clearly indicates that this genetic alteration is a major mutation for autosomal recessive deafness in Caucasoids. This finding should facilitate diagnosis of congenital deafness and allow early treatment of the affected subjects.

Geographic patterns: how to identify them and why.

Geographic patterns of genetic diversity allow us to make inferences about population histories and the evolution of inherited disease. The statistical methods describing genetic variation in space, such as estimation of genetic variances, mapping of allele frequencies, and principal components analysis, have opened up the possibility to reconstruct demographic processes whose effects have been tested by a variety of approaches, including spatial autocorrelation, cladistic analyses, and simulations. These studies have significantly contributed to our understanding of human genetic variation; however, the molecular data that have accumulated since the mid-1980s have also created new complications. Reasons include the generally limited sample sizes, but, more generally, it is the nature of molecular variation itself that makes it necessary to develop and apply specific models and methods for the treatment of DNA data. The foreseeable diffusion of laboratory techniques for the rapid typing of many DNA markers will force us to change our approach to the study of human variation anyway, moving from the gene level toward the genome level. Because extensive variation among loci is the rule rather than the exception, an important practical tip is to be skeptical of inferences based on single-locus diversity.

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.

Patterns of gene flow inferred from genetic distances in the Mediterranean region.

The analysis of population structure may lead to inferences about demographic phenomena. In particular, regions of sharp genetic differentiation suggest the existence of factors that impaired gene flow and increased the evolutionary role of genetic drift. Here, we present an analysis of a data set of 10 allele frequencies in 39 populations of the Mediterranean region. As a preliminary step, we describe spatial patterns of allele frequencies using spatial autocorrelation analysis. We then construct a network connecting localities and estimate genetic distances along the edges of the network. By applying specific algorithms, we locate on the map the areas of sharpest genetic differentiation, or genetic boundaries. The main boundaries separate the northern and the southern coasts, especially in their western portions; in addition, several localities appear genetically isolated. The comparatively high genetic differentiation across the western Mediterranean, where the sea distances between localities are shorter, strongly suggests that the sea distance by itself can hardly be regarded as a major isolating factor among these populations. On the contrary, the decrease in genetic resemblance between populations of the 2 coasts as one proceeds westward may reflect an increased genetic exchange in the eastern Mediterranean basin or independent human dispersal along the 2 coasts or both.

Clinal variation in the nuclear DNA of Europeans.

Allele frequencies are clinally distributed for many protein polymorphisms in Europe, suggesting that the current populations are derived from an ancestral group that expanded from the Near East. It is not yet fully established whether that expansion took place during the Neolithic or earlier or whether the detectable protein variation faithfully reflects the underlying molecular variation. In this study we address the latter question by describing geographic patterns of genetic diversity at seven highly polymorphic DNA markers. Two of these markers are minisatellites, four are microsatellites, and the seventh is a locus of the HLA system. By analyzing a database of 304 samples, with more than 130,000 chromosomes, we found evidence for a major clinal component of genetic variation. At most loci spatially close populations resemble each other genetically, and the degree of genetic similarity, as measured by spatial autocorrelation statistics, decreases at increasing distances. The observed patterns of molecular variation do not seem to differ qualitatively from those identified for protein polymorphisms. This suggest that low levels of population structuring, described in some mitochondrial DNA studies, may reflect different evolutionary histories for nuclear and maternally inherited markers or, alternatively, that spatial patterns of mitochondrial DNA variation may need more sensitive statistical methods to be recognized.