The genetic history of Europeans

Whole-genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.

The discovery of an in situ Neanderthal remain in the Bawa Yawan Rockshelter, West-Central Zagros Mountains, Kermanshah

Neanderthal extinction has been a matter of debate for many years. New discoveries, better chronologies and genomic evidence have done much to clarify some of the issues. This evidence suggests that Neanderthals became extinct around 40,000–37,000 years before present (BP), after a period of coexistence with Homo sapiens of several millennia, involving biological and cultural interactions between the two groups. However, the bulk of this evidence relates to Western Eurasia, and recent work in Central Asia and Siberia has shown that there is considerable local variation. Southwestern Asia, despite having a number of significant Neanderthal remains, has not played a major part in the debate over extinction. Here we report a Neanderthal deciduous canine from the site of Bawa Yawan in the West-Central Zagros Mountains of Iran. The tooth is associated with Zagros Mousterian lithics, and its context is preliminary dated to between ~43,600 and ~41,500 years ago.

Investigating Mitonuclear Genetic Interactions Through Machine Learning: A Case Study on Cold Adaptation Genes in Human Populations From Different European Climate Regions

Cold climates represent one of the major environmental challenges that anatomically modern humans faced during their dispersal out of Africa. The related adaptive traits have been achieved by modulation of thermogenesis and thermoregulation processes where nuclear (nuc) and mitochondrial (mt) genes play a major role. In human populations, mitonuclear genetic interactions are the result of both the peculiar genetic history of each human group and the different environments they have long occupied. This study aims to investigate mitonuclear genetic interactions by considering all the mitochondrial genes and 28 nuclear genes involved in brown adipose tissue metabolism, which have been previously hypothesized to be crucial for cold adaptation. For this purpose, we focused on three human populations (i.e., Finnish, British, and Central Italian people) of European ancestry from different biogeographical and climatic areas, and we used a machine learning approach to identify relevant nucDNA–mtDNA interactions that characterized each population. The obtained results are twofold: (i) at the methodological level, we demonstrated that a machine learning approach is able to detect patterns of genetic structure among human groups from different latitudes both at single genes and by considering combinations of mtDNA and nucDNA loci; (ii) at the biological level, the analysis identified population-specific nuclear genes and variants that likely play a relevant biological role in association with a mitochondrial gene (such as the “obesity gene” FTO in Finnish people). Further studies are needed to fully elucidate the evolutionary dynamics (e.g., migration, admixture, and/or local adaptation) that shaped these nucDNA–mtDNA interactions and their functional role.

The Dutch Y-chromosomal landscape

Previous studies indicated existing, albeit limited, genetic-geographic population substructure in the Dutch population based on genome-wide data and a lack of this for mitochondrial SNP based data. Despite the aforementioned studies, Y-chromosomal SNP data from the Netherlands remain scarce and do not cover the territory of the Netherlands well enough to allow a reliable investigation of genetic-geographic population substructure. Here we provide the first substantial dataset of detailed spatial Y-chromosomal haplogroup information in 2085 males collected across the Netherlands and supplemented with previously published data from northern Belgium. We found Y-chromosomal evidence for genetic–geographic population substructure, and several Y-haplogroups demonstrating significant clinal frequency distributions in different directions. By means of prediction surface maps we could visualize (complex) distribution patterns of individual Y-haplogroups in detail. These results highlight the value of a micro-geographic approach and are of great use for forensic and epidemiological investigations and our understanding of the Dutch population history. Moreover, the previously noted absence of genetic-geographic population substructure in the Netherlands based on mitochondrial DNA in contrast to our Y-chromosome results, hints at different population histories for women and men in the Netherlands.

Complex spatio-temporal distribution and genomic ancestry of mitochondrial DNA haplogroups in 24,216 Danes

Mitochondrial DNA (mtDNA) haplogroups (hgs) are evolutionarily conserved sets of mtDNA SNP-haplotypes with characteristic geographical distribution. Associations of hgs with disease and physiological characteristics have been reported, but have frequently not been reproducible. Using 418 mtDNA SNPs on the PsychChip (Illumina), we assessed the spatio-temporal distribution of mtDNA hgs in Denmark from DNA isolated from 24,642 geographically un-biased dried blood spots (DBS), collected from 1981 to 2005 through the Danish National Neonatal Screening program. ADMIXTURE was used to establish the genomic ancestry of all samples using a reference of 100K+ autosomal SNPs in 2,248 individuals from nine populations. Median-joining analysis determined that the hgs were highly variable, despite being typically Northern European in origin, suggesting multiple founder events. Furthermore, considerable heterogeneity and variation in nuclear genomic ancestry was observed. Thus, individuals with hg H exhibited 95%, and U hgs 38.2% - 92.5%, Danish ancestry. Significant clines between geographical regions and rural and metropolitan populations were found. Over 25 years, macro-hg L increased from 0.2% to 1.2% (p = 1.1*E-10), and M from 1% to 2.4% (p = 3.7*E-8). Hg U increased among the R macro-hg from 14.1% to 16.5% (p = 1.9*E-3). Genomic ancestry, geographical skewedness, and sub-hg distribution suggested that the L, M and U increases are due to immigration. The complex spatio-temporal dynamics and genomic ancestry of mtDNA in the Danish population reflect repeated migratory events and, in later years, net immigration. Such complexity may explain the often contradictory and population-specific reports of mito-genomic association with disease.

Estimating the relative importance of demic and cultural diffusion in the spread of the Neolithic in Scandinavia

Using a database of early farming sites in Scandinavia, we estimate that the spread rate of the Neolithic was in the range 0.44-0.66 km yr-1 This is substantially slower (by about 50%) than the rate in continental Europe. We interpret this result in the framework of a new mathematical model that includes horizontal cultural transmission (acculturation), vertical cultural transmission (interbreeding) and demic diffusion (reproduction and dispersal of farmers). To parametrize the model, we estimate reproduction rates of early farmers using archaeological data (sum-calibrated probabilities for the dates of early Neolithic Scandinavian sites) and use them in a wave-of-advance model for the first time. Comparing the model with the archaeological data, we find that the percentage of the spread rate due to cultural diffusion is below 50% (except for very extreme parameter values, and even for them it is below 54%). This strongly suggests that the spread of the Neolithic in Scandinavia was driven mainly by demic diffusion. This conclusion, obtained from archaeological data, agrees qualitatively with the implications of ancient genetic data, but the latter are yet too few in Scandinavia to produce any quantitative percentage for the spread rate due to cultural diffusion. We also find that, on average, fewer than eight hunter-gatherers were incorporated in the Neolithic communities by each group of 10 pioneering farmers, via horizontal and/or vertical cultural transmission.

The importance of fine-scale studies for integrating paleogenomics and archaeology

There has been an undercurrent of intellectual tension between geneticists studying human population history and archaeologists for almost 40 years. The rapid development of paleogenomics, with geneticists working on the very material discovered by archaeologists, appears to have recently heightened this tension. The relationship between these two fields thus far has largely been of a multidisciplinary nature, with archaeologists providing the raw materials for sequencing, as well as a scaffold of hypotheses based on interpretation of archaeological cultures from which the geneticists can ground their inferences from the genomic data. Much of this work has taken place in the context of western Eurasia, which is acting as testing ground for the interaction between the disciplines. Perhaps the major finding has not been any particular historical episode, but rather the apparent pervasiveness of migration events, some apparently of substantial scale, over the past ∼5000 years, challenging the prevailing view of archaeology that largely dismissed migration as a driving force of cultural change in the 1960s. However, while the genetic evidence for `migration' is generally statistically sound, the description of these events as structured behaviours is lacking, which, coupled with often over simplistic archaeological definitions, prevents the use of this information by archaeologists for studying the social processes they are interested in. In order to integrate paleogenomics and archaeology in a truly interdisciplinary manner, it will be necessary to focus less on grand narratives over space and time, and instead integrate genomic data with other form of archaeological information at the level of individual communities to understand the internal social dynamics, which can then be connected amongst communities to model migration at a regional level. A smattering of recent studies have begun to follow this approach, resulting in inferences that are not only helping ask questions that are currently relevant to archaeologists, but also potentially opening up new avenues of research.

A finely resolved phylogeny of Y chromosome Hg J illuminates the processes of Phoenician and Greek colonizations in the Mediterranean

In order to improve the phylogeography of the male-specific genetic traces of Greek and Phoenician colonizations on the Northern coasts of the Mediterranean, we performed a geographically structured sampling of seven subclades of haplogroup J in Turkey, Greece and Italy. We resequenced 4.4 Mb of Y-chromosome in 58 subjects, obtaining 1079 high quality variants. We did not find a preferential coalescence of Turkish samples to ancestral nodes, contradicting the simplistic idea of a dispersal and radiation of Hg J as a whole from the Middle East. Upon calibration with an ancient Hg J chromosome, we confirmed that signs of Holocenic Hg J radiations are subtle and date mainly to the Bronze Age. We pinpointed seven variants which could potentially unveil star clusters of sequences, indicative of local expansions. By directly genotyping these variants in Hg J carriers and complementing with published resequenced chromosomes (893 subjects), we provide strong temporal and distributional evidence for markers of the Greek settlement of Magna Graecia (J2a-L397) and Phoenician migrations (rs760148062). Our work generated a minimal but robust list of evolutionarily stable markers to elucidate the demographic dynamics and spatial domains of male-mediated movements across and around the Mediterranean, in the last 6,000 years.

A mosaic genetic structure of the human population living in the South Baltic region during the Iron Age

Despite the increase in our knowledge about the factors that shaped the genetic structure of the human population in Europe, the demographic processes that occurred during and after the Early Bronze Age (EBA) in Central-East Europe remain unclear. To fill the gap, we isolated and sequenced DNAs of 60 individuals from Kowalewko, a bi-ritual cemetery of the Iron Age (IA) Wielbark culture, located between the Oder and Vistula rivers (Kow-OVIA population). The collected data revealed high genetic diversity of Kow-OVIA, suggesting that it was not a small isolated population. Analyses of mtDNA haplogroup frequencies and genetic distances performed for Kow-OVIA and other ancient European populations showed that Kow-OVIA was most closely linked to the Jutland Iron Age (JIA) population. However, the relationship of both populations to the preceding Late Neolithic (LN) and EBA populations were different. We found that this phenomenon is most likely the consequence of the distinct genetic history observed for Kow-OVIA women and men. Females were related to the Early-Middle Neolithic farmers, whereas males were related to JIA and LN Bell Beakers. In general, our findings disclose the mechanisms that could underlie the formation of the local genetic substructures in the South Baltic region during the IA.

Reconciling evidence from ancient and contemporary genomes: a major source for the European Neolithic within Mediterranean Europe

Important gaps remain in our understanding of the spread of farming into Europe, due partly to apparent contradictions between studies of contemporary genetic variation and ancient DNA. It seems clear that farming was introduced into central, northern, and eastern Europe from the south by pioneer colonization. It is often argued that these dispersals originated in the Near East, where the potential source genetic pool resembles that of the early European farmers, but clear ancient DNA evidence from Mediterranean Europe is lacking, and there are suggestions that Mediterranean Europe may have resembled the Near East more than the rest of Europe in the Mesolithic. Here, we test this proposal by dating mitogenome founder lineages from the Near East in different regions of Europe. We find that whereas the lineages date mainly to the Neolithic in central Europe and Iberia, they largely date to the Late Glacial period in central/eastern Mediterranean Europe. This supports a scenario in which the genetic pool of Mediterranean Europe was partly a result of Late Glacial expansions from a Near Eastern refuge, and that this formed an important source pool for subsequent Neolithic expansions into the rest of Europe.

Estimating mobility using sparse data: Application to human genetic variation

Mobility is one of the most important processes shaping spatiotemporal patterns of variation in genetic, morphological, and cultural traits. However, current approaches for inferring past migration episodes in the fields of archaeology and population genetics lack either temporal resolution or formal quantification of the underlying mobility, are poorly suited to spatially and temporally sparsely sampled data, and permit only limited systematic comparison between different time periods or geographic regions. Here we present an estimator of past mobility that addresses these issues by explicitly linking trait differentiation in space and time. We demonstrate the efficacy of this estimator using spatiotemporally explicit simulations and apply it to a large set of ancient genomic data from Western Eurasia. We identify a sequence of changes in human mobility from the Late Pleistocene to the Iron Age. We find that mobility among European Holocene farmers was significantly higher than among European hunter-gatherers both pre- and postdating the Last Glacial Maximum. We also infer that this Holocene rise in mobility occurred in at least three distinct stages: the first centering on the well-known population expansion at the beginning of the Neolithic, and the second and third centering on the beginning of the Bronze Age and the late Iron Age, respectively. These findings suggest a strong link between technological change and human mobility in Holocene Western Eurasia and demonstrate the utility of this framework for exploring changes in mobility through space and time.

Mitogenome Diversity in Sardinians: A Genetic Window onto an Island's Past

Sardinians are "outliers" in the European genetic landscape and, according to paleogenomic nuclear data, the closest to early European Neolithic farmers. To learn more about their genetic ancestry, we analyzed 3,491 modern and 21 ancient mitogenomes from Sardinia. We observed that 78.4% of modern mitogenomes cluster into 89 haplogroups that most likely arose in situ. For each Sardinian-specific haplogroup (SSH), we also identified the upstream node in the phylogeny, from which non-Sardinian mitogenomes radiate. This provided minimum and maximum time estimates for the presence of each SSH on the island. In agreement with demographic evidence, almost all SSHs coalesce in the post-Nuragic, Nuragic and Neolithic-Copper Age periods. For some rare SSHs, however, we could not dismiss the possibility that they might have been on the island prior to the Neolithic, a scenario that would be in agreement with archeological evidence of a Mesolithic occupation of Sardinia.

Population resequencing of European mitochondrial genomes highlights sex-bias in Bronze Age demographic expansions

Interpretations of genetic data concerning the prehistory of Europe have long been a subject of great debate, but increasing amounts of ancient and modern DNA data are now providing new and more informative evidence. Y-chromosome resequencing studies in Europe have highlighted the prevalence of recent expansions of male lineages, and focused interest on the Bronze Age as a period of cultural and demographic change. These findings contrast with phylogeographic studies based on mitochondrial DNA (mtDNA), which have been interpreted as supporting expansions from glacial refugia. Here we have undertaken a population-based resequencing of complete mitochondrial genomes in Europe and the Middle East, in 340 samples from 17 populations for which Y-chromosome sequence data are also available. Demographic reconstructions show no signal of Bronze Age expansion, but evidence of Paleolithic expansions in all populations except the Saami, and with an absence of detectable geographical pattern. In agreement with previous inference from modern and ancient DNA data, the unbiased comparison between the mtDNA and Y-chromosome population datasets emphasizes the sex-biased nature of recent demographic transitions in Europe.

Identification and analysis of mtDNA genomes attributed to Finns reveal long-stagnant demographic trends obscured in the total diversity

In Europe, modern mitochondrial diversity is relatively homogeneous and suggests an ubiquitous rapid population growth since the Neolithic revolution. Similar patterns also have been observed in mitochondrial control region data in Finland, which contrasts with the distinctive autosomal and Y-chromosomal diversity among Finns. A different picture emerges from the 843 whole mitochondrial genomes from modern Finns analyzed here. Up to one third of the subhaplogroups can be considered as Finn-characteristic, i.e. rather common in Finland but virtually absent or rare elsewhere in Europe. Bayesian phylogenetic analyses suggest that most of these attributed Finnish lineages date back to around 3,000–5,000 years, coinciding with the arrival of Corded Ware culture and agriculture into Finland. Bayesian estimation of past effective population sizes reveals two differing demographic histories: 1) the ‘local’ Finnish mtDNA haplotypes yielding small and dwindling size estimates for most of the past; and 2) the ‘immigrant’ haplotypes showing growth typical of most European populations. The results based on the local diversity are more in line with that known about Finns from other studies, e.g., Y-chromosome analyses and archaeology findings. The mitochondrial gene pool thus may contain signals of local population history that cannot be readily deduced from the total diversity.

Detecting past male-mediated expansions using the Y chromosome

Males and females display biological differences that lead to a higher variance of offspring number in males, and this is frequently exacerbated in human societies by mating practices, and possibly by past socio-cultural circumstances. This implies that the genetic record might contain the imprint of past male-mediated expansions, which can be investigated by analysing the male-specific region of the Y chromosome (MSY). Here, we review studies that have used MSY data to infer such expansions. Sets of short-tandem repeats define haplotypes of very low average frequencies, but in a few cases, high-frequency haplotypes are observed, forming the core of descent clusters. Estimates of the ages of such clusters, together with geographical information, have been used to propose powerful historical founders, including Genghis Khan, although without direct supporting evidence. Resequencing of multi-megabase segments of MSY has allowed the construction of detailed phylogenies in which branch lengths are proportional to time, leading to the identification of lineage expansions in the last few millennia as well as the more distant past. Comparisons with maternally-inherited mitochondrial DNA sequence data allow the male specificity of some of these expansions to be demonstrated. These include expansions in Europe in the last ~5000 years that may be associated with a cultural shift during the Bronze Age, as well as expansions elsewhere in the world for which explanations from archaeological evidence are not yet clear.

Complete mitochondrial sequences from Mesolithic Sardinia

Little is known about the genetic prehistory of Sardinia because of the scarcity of pre-Neolithic human remains. From a genetic perspective, modern Sardinians are known as genetic outliers in Europe, showing unusually high levels of internal diversity and a close relationship to early European Neolithic farmers. However, how far this peculiar genetic structure extends and how it originated was to date impossible to test. Here we present the first and oldest complete mitochondrial sequences from Sardinia, dated back to 10,000 yBP. These two individuals, while confirming a Mesolithic occupation of the island, belong to rare mtDNA lineages, which have never been found before in Mesolithic samples and that are currently present at low frequencies not only in Sardinia, but in the whole Europe. Preliminary Approximate Bayesian Computations, restricted by biased reference samples for Mesolithic Sardinia (the two typed samples) and Neolithic Europe (limited to central and north European sequences), suggest that the first inhabitants of the island have had a small or negligible contribution to the present-day Sardinian population, which mainly derives its genetic diversity from continental migration into the island by Neolithic times.

[aDNA Research From a Historical Perspective]

aDNA studies are a cooperative field of research with a broad range of applications including evolutionary biology, genetics, anthropology and archaeology. Scientists are using ancient molecules as source material for historical questions. Colleagues from the humanities are observing this with both interest and concern because aDNA research is affecting academic identities and both concepts of history and historiography. aDNA research developed in a way that can be described as a Hype Cycle (Chackie Fenn). Technological triggers such as Sanger Sequencing and the Polymerase Chain Reaction kicked off a multitude of experiments with ancient DNA during the 1980s and 1990s. Geneticists, microbiologists, anthropologists and many more euphorically joined a "molecule hunt". aDNA was promoted as a time machine. Media attention was enormous. As experiments and implementations began to fail and contamination was discovered to be a tremendous problem, media interest waned and many labs lost their interest. Some turned their disillusionment into systematic research into methodology and painstakingly established lab routines. The authenticity problem was first addressed by control oriented measures but later approached from a more cognitive theoretical perspective as the pitfalls and limits of aDNA became clearer. By the end of the 2000s the field reached its current plateau of productivity. Cross-disciplinary debates, conflicts and collaborations are increasing critical reflection among all participants. Historians should consider joining the field in a kind of critical friendship to both make the most of its possibilities and give an input from a constructivist perspective.

Spatially Explicit Models to Investigate Geographic Patterns in the Distribution of Forensic STRs: Application to the North-Eastern Mediterranean

Human forensic STRs used for individual identification have been reported to have little power for inter-population analyses. Several methods have been developed which incorporate information on the spatial distribution of individuals to arrive at a description of the arrangement of diversity. We genotyped at 16 forensic STRs a large population sample obtained from many locations in Italy, Greece and Turkey, i.e. three countries crucial to the understanding of discontinuities at the European/Asian junction and the genetic legacy of ancient migrations, but seldom represented together in previous studies. Using spatial PCA on the full dataset, we detected patterns of population affinities in the area. Additionally, we devised objective criteria to reduce the overall complexity into reduced datasets. Independent spatially explicit methods applied to these latter datasets converged in showing that the extraction of information on long- to medium-range geographical trends and structuring from the overall diversity is possible. All analyses returned the picture of a background clinal variation, with regional discontinuities captured by each of the reduced datasets. Several aspects of our results are confirmed on external STR datasets and replicate those of genome-wide SNP typings. High levels of gene flow were inferred within the main continental areas by coalescent simulations. These results are promising from a microevolutionary perspective, in view of the fast pace at which forensic data are being accumulated for many locales. It is foreseeable that this will allow the exploitation of an invaluable genotypic resource, assembled for other (forensic) purposes, to clarify important aspects in the formation of local gene pools.

Demography of the Early Neolithic Population in Central Balkans: Population Dynamics Reconstruction Using Summed Radiocarbon Probability Distributions

The Central Balkans region is of great importance for understanding the spread of the Neolithic in Europe but the Early Neolithic population dynamics of the region is unknown. In this study we apply the method of summed calibrated probability distributions to a set of published radiocarbon dates from the Republic of Serbia in order to reconstruct population dynamics in the Early Neolithic in this part of the Central Balkans. The results indicate that there was a significant population growth after ~6200 calBC, when the Neolithic was introduced into the region, followed by a bust at the end of the Early Neolithic phase (~5400 calBC). These results are broadly consistent with the predictions of the Neolithic Demographic Transition theory and the patterns of population booms and busts detected in other regions of Europe. These results suggest that the cultural process that underlies the patterns observed in Central and Western Europe was also in operation in the Central Balkan Neolithic and that the population increase component of this process can be considered as an important factor for the spread of the Neolithic as envisioned in the demic diffusion hypothesis.

Population structure from NOS genes correlates with geographical differences in coronary incidence across Europe

OBJECTIVES

The population analysis of cardiovascular risk and non-risk genetic variation can help to identify adaptive or random demographic processes that shaped coronary incidence variation across geography.

MATERIAL AND METHODS

In this study, 114 single nucleotide polymorphisms and 17 tandem repeat polymorphisms from Nitric Oxide Synthases (NOS) regions were analyzed in 1686 individuals from 35 populations from Europe, North Africa, and the Middle East. NOS genes encode for key enzymes on nitric oxide availability, which is involved in several cardiovascular processes. These genetic variations were used to test for selection and to infer the population structure of NOS regions. Moreover, we tested whether the variation in the incidence of coronary events and in the levels of classical risk factors in 11 of these European populations could be explained by the population structure estimates.

RESULTS

Our results supported, first, the absence of clear signs of selection for NOS genetic variants associated with cardiovascular diseases, and second, the presence of a continuous genetic pattern of variation across European and North African populations without a Mediterranean barrier for gene flow. Finally, population structure estimates from NOS regions are closely correlated with coronary event rates and classical risk parameters (explaining 39-98%) among European populations.

CONCLUSION

Our results reinforce the hypothesis that genetic bases of cardiovascular diseases and associated complex phenotypes could be geographically shaped by random demographic processes.

Maternal Genetic Composition of a Medieval Population from a Hungarian-Slavic Contact Zone in Central Europe

The genetic composition of the medieval populations of Central Europe has been poorly investigated to date. In particular, the region of modern-day Slovakia is a blank spot in archaeogenetic research. This paper reports the study of mitochondrial DNA (mtDNA) in ancient samples from the 9^th–12^th centuries originating from the cemeteries discovered in Nitra-Šindolka and Čakajovce, located in western Slovakia (Central Europe). This geographical region is interesting to study because its medieval multi-ethnic population lived in the so-called contact zone of the territory of the Great Moravian and later Hungarian state formations. We described 16 different mtDNA haplotypes in 19 individuals, which belong to the most widespread European mtDNA haplogroups: H, J, T, U and R0. Using comparative statistical and population genetic analyses, we showed the differentiation of the European gene pool in the medieval period. We also demonstrated the heterogeneous genetic characteristics of the investigated population and its affinity to the populations of modern Europe.

Y-chromosome phylogeographic analysis of the Greek-Cypriot population reveals elements consistent with Neolithic and Bronze Age settlements

Background

The archeological record indicates that the permanent settlement of Cyprus began with pioneering agriculturalists circa 11,000 years before present, (ca. 11,000 y BP). Subsequent colonization events followed, some recognized regionally. Here, we assess the Y-chromosome structure of Cyprus in context to regional populations and correlate it to phases of prehistoric colonization.

Results

Analysis of haplotypes from 574 samples showed that island-wide substructure was barely significant in a spatial analysis of molecular variance (SAMOVA). However, analyses of molecular variance (AMOVA) of haplogroups using 92 binary markers genotyped in 629 Cypriots revealed that the proportion of variance among the districts was irregularly distributed. Principal component analysis (PCA) revealed potential genetic associations of Greek-Cypriots with neighbor populations. Contrasting haplogroups in the PCA were used as surrogates of parental populations. Admixture analyses suggested that the majority of G2a-P15 and R1b-M269 components were contributed by Anatolia and Levant sources, respectively, while Greece Balkans supplied the majority of E-V13 and J2a-M67. Haplotype-based expansion times were at historical levels suggestive of recent demography.

Conclusions

Analyses of Cypriot haplogroup data are consistent with two stages of prehistoric settlement. E-V13 and E-M34 are widespread, and PCA suggests sourcing them to the Balkans and Levant/Anatolia, respectively. The persistent pre-Greek component is represented by elements of G2-U5(xL30) haplogroups: U5*, PF3147, and L293. J2b-M205 may contribute also to the pre-Greek strata. The majority of R1b-Z2105 lineages occur in both the westernmost and easternmost districts. Distinctively, sub-haplogroup R1b- M589 occurs only in the east. The absence of R1b- M589 lineages in Crete and the Balkans and the presence in Asia Minor are compatible with Late Bronze Age influences from Anatolia rather than from Mycenaean Greeks.

Electronic supplementary material

The online version of this article (doi:10.1186/s13323-016-0032-8) contains supplementary material, which is available to authorized users.

Estimating time to the most recent common ancestor (TMRCA): comparison and application of eight methods

Investigating how an ancestral population diverges to give rise to distinct subpopulations remains a fundamental pursuit in population genetics. There is broad consensus for the 'Out-of-Africa' hypothesis that states that modern humans arose ∼200 000 years ago in Africa and spread throughout the continent ∼100 000 years ago. This was followed by several waves of major population dispersals across the globe, although the exact nature of the population divergence remains debatable. Existing methods to estimate population divergence time differ in their methodological frameworks and demographic assumptions, and require different types of genetic data as input. These fundamental differences often result in the methods producing inconsistent estimates of the population divergence time, further confounding attempts to robustly uncover the history of human migration, especially when most population genetic studies do not employ multiple methods to estimate the time to the most recent common ancestor (TMRCA). Here, we chose eight popular methods for estimating TMRCA and evaluated their robustness and accuracy in correctly identifying the true TMRCA through a series of simulations that mimicked different evolutionary scenarios. We subsequently applied all eight methods to estimate the population divergence time between Southeast Asian Malays and South Asian Indians using deep whole-genome sequencing data.

Fine Dissection of Human Mitochondrial DNA Haplogroup HV Lineages Reveals Paleolithic Signatures from European Glacial Refugia

Genetic signatures from the Paleolithic inhabitants of Eurasia can be traced from the early divergent mitochondrial DNA lineages still present in contemporary human populations. Previous studies already suggested a pre-Neolithic diffusion of mitochondrial haplogroup HV*(xH,V) lineages, a relatively rare class of mtDNA types that includes parallel branches mainly distributed across Europe and West Asia with a certain degree of structure. Up till now, variation within haplogroup HV was addressed mainly by analyzing sequence data from the mtDNA control region, except for specific sub-branches, such as HV4 or the widely distributed haplogroups H and V. In this study, we present a revised HV topology based on full mtDNA genome data, and we include a comprehensive dataset consisting of 316 complete mtDNA sequences including 60 new samples from the Italian peninsula, a previously underrepresented geographic area. We highlight points of instability in the particular topology of this haplogroup, reconstructed with BEAST-generated trees and networks. We also confirm a major lineage expansion that probably followed the Late Glacial Maximum and preceded Neolithic population movements. We finally observe that Italy harbors a reservoir of mtDNA diversity, with deep-rooting HV lineages often related to sequences present in the Caucasus and the Middle East. The resulting hypothesis of a glacial refugium in Southern Italy has implications for the understanding of late Paleolithic population movements and is discussed within the archaeological cultural shifts occurred over the entire continent.

Long-Distance Dispersal Shaped Patterns of Human Genetic Diversity in Eurasia

Most previous attempts at reconstructing the past history of human populations did not explicitly take geography into account or considered very simple scenarios of migration and ignored environmental information. However, it is likely that the last glacial maximum (LGM) affected the demography and the range of many species, including our own. Moreover, long-distance dispersal (LDD) may have been an important component of human migrations, allowing fast colonization of new territories and preserving high levels of genetic diversity. Here, we use a high-quality microsatellite data set genotyped in 22 populations to estimate the posterior probabilities of several scenarios for the settlement of the Old World by modern humans. We considered models ranging from a simple spatial expansion to others including LDD and a LGM-induced range contraction, as well as Neolithic demographic expansions. We find that scenarios with LDD are much better supported by data than models without LDD. Nevertheless, we show evidence that LDD events to empty habitats were strongly prevented during the settlement of Eurasia. This unexpected absence of LDD ahead of the colonization wave front could have been caused by an Allee effect, either due to intrinsic causes such as an inbreeding depression built during the expansion or due to extrinsic causes such as direct competition with archaic humans. Overall, our results suggest only a relatively limited effect of the LGM contraction on current patterns of human diversity. This is in clear contrast with the major role of LDD migrations, which have potentially contributed to the intermingled genetic structure of Eurasian populations.

Upper Palaeolithic genomes reveal deep roots of modern Eurasians

We extend the scope of European palaeogenomics by sequencing the genomes of Late Upper Palaeolithic (13,300 years old, 1.4-fold coverage) and Mesolithic (9,700 years old, 15.4-fold) males from western Georgia in the Caucasus and a Late Upper Palaeolithic (13,700 years old, 9.5-fold) male from Switzerland. While we detect Late Palaeolithic-Mesolithic genomic continuity in both regions, we find that Caucasus hunter-gatherers (CHG) belong to a distinct ancient clade that split from western hunter-gatherers ∼45 kya, shortly after the expansion of anatomically modern humans into Europe and from the ancestors of Neolithic farmers ∼25 kya, around the Last Glacial Maximum. CHG genomes significantly contributed to the Yamnaya steppe herders who migrated into Europe ∼3,000 BC, supporting a formative Caucasus influence on this important Early Bronze age culture. CHG left their imprint on modern populations from the Caucasus and also central and south Asia possibly marking the arrival of Indo-Aryan languages.

Admixture facilitates genetic adaptations to high altitude in Tibet

Admixture is recognized as a widespread feature of human populations, renewing interest in the possibility that genetic exchange can facilitate adaptations to new environments. Studies of Tibetans revealed candidates for high-altitude adaptations in the EGLN1 and EPAS1 genes, associated with lower haemoglobin concentration. However, the history of these variants or that of Tibetans remains poorly understood. Here we analyse genotype data for the Nepalese Sherpa, and find that Tibetans are a mixture of ancestral populations related to the Sherpa and Han Chinese. EGLN1 and EPAS1 genes show a striking enrichment of high-altitude ancestry in the Tibetan genome, indicating that migrants from low altitude acquired adaptive alleles from the highlanders. Accordingly, the Sherpa and Tibetans share adaptive haemoglobin traits. This admixture-mediated adaptation shares important features with adaptive introgression. Therefore, we identify a novel mechanism, beyond selection on new mutations or on standing variation, through which populations can adapt to local environments.

Ancient mitochondrial DNA from the northern fringe of the Neolithic farming expansion in Europe sheds light on the dispersion process

The European Neolithization process started around 12 000 years ago in the Near East. The introduction of agriculture spread north and west throughout Europe and a key question has been if this was brought about by migrating individuals, by an exchange of ideas or a by a mixture of these. The earliest farming evidence in Scandinavia is found within the Funnel Beaker Culture complex (Trichterbecherkultur, TRB) which represents the northernmost extension of Neolithic farmers in Europe. The TRB coexisted for almost a millennium with hunter-gatherers of the Pitted Ware Cultural complex (PWC). If migration was a substantial part of the Neolithization, even the northerly TRB community would display a closer genetic affinity to other farmer populations than to hunter-gatherer populations. We deep-sequenced the mitochondrial hypervariable region 1 from seven farmers (six TRB and one Battle Axe complex, BAC) and 13 hunter-gatherers (PWC) and authenticated the sequences using postmortem DNA damage patterns. A comparison with 124 previously published sequences from prehistoric Europe shows that the TRB individuals share a close affinity to Central European farmer populations, and that they are distinct from hunter-gatherer groups, including the geographically close and partially contemporary PWC that show a close affinity to the European Mesolithic hunter-gatherers.

Vestiges of an Ancient Border in the Contemporary Genetic Diversity of North-Eastern Europe

It has previously been demonstrated that the advance of the Neolithic Revolution from the Near East through Europe was decelerated in the northernmost confines of the continent, possibly as a result of space and resource competition with lingering Mesolithic populations. Finland was among the last domains to adopt a farming lifestyle, and is characterized by substructuring in the form of a distinct genetic border dividing the northeastern and southwestern regions of the country. To explore the origins of this divergence, the geographical patterns of mitochondrial and Y-chromosomal haplogroups of Neolithic and Mesolithic ancestry were assessed in Finnish populations. The distribution of these uniparental markers revealed a northeastern bias for hunter-gatherer haplogroups, while haplogroups associated with the farming lifestyle clustered in the southwest. In addition, a correlation could be observed between more ancient mitochondrial haplogroup age and eastern concentration. These results coupled with prior archeological evidence suggest the genetic northeast/southwest division observed in contemporary Finland represents an ancient vestigial border between Mesolithic and Neolithic populations undetectable in most other regions of Europe.

New clues to the evolutionary history of the main European paternal lineage M269: dissection of the Y-SNP S116 in Atlantic Europe and Iberia

The dissection of S116 in more than 1500 individuals from Atlantic Europe and the Iberian Peninsula has provided important clues about the controversial evolutionary history of M269. First, the results do not point to an origin of M269 in the Franco-Cantabrian refuge, owing to the lack of sublineage diversity within M269, which supports the new theories proposing its origin in Eastern Europe. Second, S116 shows frequency peaks and spatial distribution that differ from those previously proposed, indicating an origin farther west, and it also shows a high frequency in the Atlantic coastline. Third, an outstanding frequency of the DF27 sublineage has been found in Iberia, with a restricted distribution pattern inside this peninsula and a frequency maximum in the area of the Franco-Cantabrian refuge. This entire panorama indicates an old arrival of M269 into Western Europe, because it has generated at least two episodes of expansion in the Franco-Cantabrian area. This study demonstrates the importance of continuing the dissection of the M269 lineage in different European populations because the discovery and study of new sublineages can adjust or even completely revise the theories about European peopling, as has been the case for the place of origin of M269.

The maternal perspective for five Slovenian regions: The importance of regional sampling

BACKGROUND

The Slovenian territory is geographically positioned between the Alps, Adriatic Sea, Pannonian basin and the Dinaric Mountains and, as such, has served as a passageway for various populations in different periods of time. Turbulent historic events and diverse geography of the region have produced a diverse contemporary population whose genetic analysis could provide insight into past demographic events.

AIM

The aims of this study were to characterize the Slovenian mitochondrial gene pool at the micro-geographic level and to compare it with surrounding populations.

SUBJECTS AND METHODS

A total of 402 individuals from five Slovenian regions were analysed in this study by typing HVR I, HVR II and coding region polymorphisms of mtDNA.

RESULTS

Analysis revealed 47 haplogroups and sub-haplogroups, the most common of which were H*, H1, J1c, T2 and U5a. Intra-population comparisons revealed a sharp gradient of the J1c haplogroup between Slovenian regions, with a peak frequency of 24.5% being observed in the population of the Littoral Region.

CONCLUSION

The sharp gradient of the J1c haplogroup between Slovenian regions is in line with the archaeological horizon known as Impressed Ware culture and could, therefore, represent a genetic trace of the early Neolithic expansion route along the East Adriatic coastal region.

Ancient DNA from South-East Europe Reveals Different Events during Early and Middle Neolithic Influencing the European Genetic Heritage

The importance of the process of Neolithization for the genetic make-up of European populations has been hotly debated, with shifting hypotheses from a demic diffusion (DD) to a cultural diffusion (CD) model. In this regard, ancient DNA data from the Balkan Peninsula, which is an important source of information to assess the process of Neolithization in Europe, is however missing. In the present study we show genetic information on ancient populations of the South-East of Europe. We assessed mtDNA from ten sites from the current territory of Romania, spanning a time-period from the Early Neolithic to the Late Bronze Age. mtDNA data from Early Neolithic farmers of the Starčevo Criş culture in Romania (Cârcea, Gura Baciului and Negrileşti sites), confirm their genetic relationship with those of the LBK culture (Linienbandkeramik Kultur) in Central Europe, and they show little genetic continuity with modern European populations. On the other hand, populations of the Middle-Late Neolithic (Boian, Zau and Gumelniţa cultures), supposedly a second wave of Neolithic migration from Anatolia, had a much stronger effect on the genetic heritage of the European populations. In contrast, we find a smaller contribution of Late Bronze Age migrations to the genetic composition of Europeans. Based on these findings, we propose that permeation of mtDNA lineages from a second wave of Middle-Late Neolithic migration from North-West Anatolia into the Balkan Peninsula and Central Europe represent an important contribution to the genetic shift between Early and Late Neolithic populations in Europe, and consequently to the genetic make-up of modern European populations.

The genetic prehistory of domesticated cattle from their origin to the spread across Europe

Background

Cattle domestication started in the 9^th millennium BC in Southwest Asia. Domesticated cattle were then introduced into Europe during the Neolithic transition. However, the scarcity of palaeogenetic data from the first European domesticated cattle still inhibits the accurate reconstruction of their early demography. In this study, mitochondrial DNA from 193 ancient and 597 modern domesticated cattle (Bos taurus) from sites across Europe, Western Anatolia and Iran were analysed to provide insight into the Neolithic dispersal process and the role of the local European aurochs population during cattle domestication.

Results

Using descriptive summary statistics and serial coalescent simulations paired with approximate Bayesian computation we find: (i) decreasing genetic diversity in a southeast to northwest direction, (ii) strong correlation of genetic and geographical distances, iii) an estimated effective size of the Near Eastern female founder population of 81, iv) that the expansion of cattle from the Near East and Anatolia into Europe does not appear to constitute a significant bottleneck, and that v) there is evidence for gene-flow between the Near Eastern/Anatolian and European cattle populations in the early phases of the European Neolithic, but that it is restricted after 5,000 BCE.

Conclusions

The most plausible scenario to explain these results is a single and regionally restricted domestication process of cattle in the Near East with subsequent migration into Europe during the Neolithic transition without significant maternal interbreeding with the endogenous wild stock. Evidence for gene-flow between cattle populations from Southwestern Asia and Europe during the earlier phases of the European Neolithic points towards intercontinental trade connections between Neolithic farmers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0203-2) contains supplementary material, which is available to authorized users.

Mitochondrial DNA perspective of Serbian genetic diversity

Although south-Slavic populations have been studied to date from various aspects, the population of Serbia, occupying the central part of the Balkan Peninsula, is still genetically understudied at least at the level of mitochondrial DNA (mtDNA) variation. We analyzed polymorphisms of the first and the second mtDNA hypervariable segments (HVS-I and HVS-II) and informative coding-region markers in 139 Serbians to shed more light on their mtDNA variability, and used available data on other Slavic and neighboring non-Slavic populations to assess their interrelations in a broader European context. The contemporary Serbian mtDNA profile is consistent with the general European maternal landscape having a substantial proportion of shared haplotypes with eastern, central, and southern European populations. Serbian population was characterized as an important link between easternmost and westernmost south-Slavic populations due to the observed lack of genetic differentiation with all other south-Slavic populations and its geographical positioning within the Balkan Peninsula. An increased heterogeneity of south Slavs, most likely mirroring turbulent demographic events within the Balkan Peninsula over time (i.e., frequent admixture and differential introgression of various gene pools), and a marked geographical stratification of Slavs to south-, east-, and west-Slavic groups, were also found. A phylogeographic analyses of 20 completely sequenced Serbian mitochondrial genomes revealed not only the presence of mtDNA lineages predominantly found within the Slavic gene pool (U4a2a*, U4a2a1, U4a2c, U4a2g, HV10), supporting a common Slavic origin, but also lineages that may have originated within the southern Europe (H5*, H5e1, H5a1v) and the Balkan Peninsula in particular (H6a2b and L2a1k).

Human paleogenetics of Europe--the known knowns and the known unknowns

The number of ancient human DNA studies has drastically increased in recent years. This results in a substantial record of mitochondrial sequences available from many prehistoric sites across Western Eurasia, but also growing Y-chromosome and autosomal sequence data. We review the current state of research with specific emphasis on the Holocene population events that likely have shaped the present-day genetic variation in Europe. We reconcile observations from the genetic data with hypotheses about the peopling and settlement history from anthropology and archaeology for various key regions, and also discuss the data in light of evidence from related disciplines, such as modern human genetics, climatology and linguistics.

Toward a new history and geography of human genes informed by ancient DNA

Genetic information contains a record of the history of our species, and technological advances have transformed our ability to access this record. Many studies have used genome-wide data from populations today to learn about the peopling of the globe and subsequent adaptation to local conditions. Implicit in this research is the assumption that the geographic locations of people today are informative about the geographic locations of their ancestors in the distant past. However, it is now clear that long-range migration, admixture, and population replacement subsequent to the initial out-of-Africa expansion have altered the genetic structure of most of the world's human populations. In light of this we argue that it is time to critically reevaluate current models of the peopling of the globe, as well as the importance of natural selection in determining the geographic distribution of phenotypes. We specifically highlight the transformative potential of ancient DNA. By accessing the genetic make-up of populations living at archaeologically known times and places, ancient DNA makes it possible to directly track migrations and responses to natural selection.

Hypothesis: Possible respiratory advantages for heterozygote carriers of cystic fibrosis linked mutations during dusty climate of last glaciation

This paper puts forward a new hypothesis to interpret the high carrier frequency of CFTR mutations in individuals of European descent. The proposed heterozygote advantage factor is related to the specific climate conditions in Europe during the last 50 ky that might have heavily compromised the respiratory function of our ancestors in Eurasia. A large part of the last 50 ky was cold, and the coldest period was the Last Glacial Maximum (LGM) (26.5 to 19 kya). The global climate was dry with a dust-laden atmosphere (20 to 25 times more dust than the present level). High levels of atmospheric dust started more than 40 kya and ended less than 10 kya. Secretion of airway fluid is usually related to the submucosal tissue hydration, while salt reabsorption relies on activation of CFTRs that allow ENaCs to absorb salt and water. The water loss by evaporation depends on the air humidity and flow rate. Salt accumulation in the mucus is normally prevented by reabsorption of Na(+) and Cl(-) by epithelial cells if the presence of functional CFTRs is normal. If one gene for CFTR is mutated, the number of functional CFTRs is reduced and this limits the capacity of salt reabsorption by epithelial cells. This means that evaporation makes the airway fluid more hypertonic, and osmotic forces bring more water from the interstitial space, thus leading to a new balance in mucosal fluid traffic. Increased osmolarity and volume of airway fluid can be more moveable in cases when evaporation and dust exposure is increased. If both CFTR genes are mutated, low number of functional CFTRs diminishes salt resorption of epithelial cells. Salt accumulated in the mucous fluid within respiratory ducts, as previously described. The hypertonic ductal content forces more water and some electrolytes to enter the airway fluid from the interstitial fluid, and evaporation leads to further concentration of thick immobile mucus. The proposed interpretation is that CFTR mutations have spread among our ancestors that roamed the central Eurasia after the LGM. The heterozygote individuals might have benefitted from the limited water resorption in their respiratory mucosa that allowed improved airway cleansing.

Mitochondrial DNA from El Mirador cave (Atapuerca, Spain) reveals the heterogeneity of Chalcolithic populations

Previous mitochondrial DNA analyses on ancient European remains have suggested that the current distribution of haplogroup H was modeled by the expansion of the Bell Beaker culture (ca 4,500–4,050 years BP) out of Iberia during the Chalcolithic period. However, little is known on the genetic composition of contemporaneous Iberian populations that do not carry the archaeological tool kit defining this culture. Here we have retrieved mitochondrial DNA (mtDNA) sequences from 19 individuals from a Chalcolithic sample from El Mirador cave in Spain, dated to 4,760–4,200 years BP and we have analyzed the haplogroup composition in the context of modern and ancient populations. Regarding extant African, Asian and European populations, El Mirador shows affinities with Near Eastern groups. In different analyses with other ancient samples, El Mirador clusters with Middle and Late Neolithic populations from Germany, belonging to the Rössen, the Salzmünde and the Baalberge archaeological cultures but not with contemporaneous Bell Beakers. Our analyses support the existence of a common genetic signal between Western and Central Europe during the Middle and Late Neolithic and points to a heterogeneous genetic landscape among Chalcolithic groups.

Craniometric analysis of European Upper Palaeolithic and Mesolithic samples supports discontinuity at the Last Glacial Maximum

The Last Glacial Maximum (LGM) represents the most significant climatic event since the emergence of anatomically modern humans (AMH). In Europe, the LGM may have played a role in changing morphological features as a result of adaptive and stochastic processes. We use craniometric data to examine morphological diversity in pre- and post-LGM specimens. Craniometric variation is assessed across four periods--pre-LGM, late glacial, Early Holocene and Middle Holocene--using a large, well-dated, data set. Our results show significant differences across the four periods, using a MANOVA on size-adjusted cranial measurements. A discriminant function analysis shows separation between pre-LGM and later groups. Analyses repeated on a subsample, controlled for time and location, yield similar results. The results are largely influenced by facial measurements and are most consistent with neutral demographic processes. These findings suggest that the LGM had a major impact on AMH populations in Europe prior to the Neolithic.

Improved phylogenetic resolution and rapid diversification of Y-chromosome haplogroup K-M526 in Southeast Asia

The highly structured distribution of Y-chromosome haplogroups suggests that current patterns of variation may be informative of past population processes. However, limited phylogenetic resolution, particularly of subclades within haplogroup K, has obscured the relationships of lineages that are common across Eurasia. Here we genotype 13 new highly informative single-nucleotide polymorphisms in a worldwide sample of 4413 males that carry the derived allele at M526, and reconstruct an NRY haplogroup tree with significantly higher resolution for the major clade within haplogroup K, K-M526. Although K-M526 was previously characterized by a single polytomy of eight major branches, the phylogenetic structure of haplogroup K-M526 is now resolved into four major subclades (K2a-d). The largest of these subclades, K2b, is divided into two clusters: K2b1 and K2b2. K2b1 combines the previously known haplogroups M, S, K-P60 and K-P79, whereas K2b2 comprises haplogroups P and its subhaplogroups Q and R. Interestingly, the monophyletic group formed by haplogroups R and Q, which make up the majority of paternal lineages in Europe, Central Asia and the Americas, represents the only subclade with K2b that is not geographically restricted to Southeast Asia and Oceania. Estimates of the interval times for the branching events between M9 and P295 point to an initial rapid diversification process of K-M526 that likely occurred in Southeast Asia, with subsequent westward expansions of the ancestors of haplogroups R and Q.

Population genomic analysis of ancient and modern genomes yields new insights into the genetic ancestry of the Tyrolean Iceman and the genetic structure of Europe

Genome sequencing of the 5,300-year-old mummy of the Tyrolean Iceman, found in 1991 on a glacier near the border of Italy and Austria, has yielded new insights into his origin and relationship to modern European populations. A key finding of that study was an apparent recent common ancestry with individuals from Sardinia, based largely on the Y chromosome haplogroup and common autosomal SNP variation. Here, we compiled and analyzed genomic datasets from both modern and ancient Europeans, including genome sequence data from over 400 Sardinians and two ancient Thracians from Bulgaria, to investigate this result in greater detail and determine its implications for the genetic structure of Neolithic Europe. Using whole-genome sequencing data, we confirm that the Iceman is, indeed, most closely related to Sardinians. Furthermore, we show that this relationship extends to other individuals from cultural contexts associated with the spread of agriculture during the Neolithic transition, in contrast to individuals from a hunter-gatherer context. We hypothesize that this genetic affinity of ancient samples from different parts of Europe with Sardinians represents a common genetic component that was geographically widespread across Europe during the Neolithic, likely related to migrations and population expansions associated with the spread of agriculture.

The analysis of the genome of the Tyrolean Iceman, a 5,300 year old mummy from Central Europe, revealed a surprising recent common ancestry with modern Sardinians for this ancient genome. However, this study was limited both by the availability of data from Sardinians and by a lack of genomic data from other ancient European samples. Here, we use genomic data from modern Sardinians and from ancient European individuals from different geographic regions and cultural contexts, to demonstrate that this ancestry component is shared among individuals associated with the onset of agriculture in Europe. Our results thus suggest that the Iceman's Sardinian ancestry actually reflects a more widespread genetic component related to the migration of people during the Neolithic transition in Central Europe.

Regional population collapse followed initial agriculture booms in mid-Holocene Europe

Following its initial arrival in SE Europe 8,500 years ago agriculture spread throughout the continent, changing food production and consumption patterns and increasing population densities. Here we show that, in contrast to the steady population growth usually assumed, the introduction of agriculture into Europe was followed by a boom-and-bust pattern in the density of regional populations. We demonstrate that summed calibrated radiocarbon date distributions and simulation can be used to test the significance of these demographic booms and busts in the context of uncertainty in the radiocarbon date calibration curve and archaeological sampling. We report these results for Central and Northwest Europe between 8,000 and 4,000 cal. BP and investigate the relationship between these patterns and climate. However, we find no evidence to support a relationship. Our results thus suggest that the demographic patterns may have arisen from endogenous causes, although this remains speculative.

The phylogenetic and geographic structure of Y-chromosome haplogroup R1a

R1a-M420 is one of the most widely spread Y-chromosome haplogroups; however, its substructure within Europe and Asia has remained poorly characterized. Using a panel of 16 244 male subjects from 126 populations sampled across Eurasia, we identified 2923 R1a-M420 Y-chromosomes and analyzed them to a highly granular phylogeographic resolution. Whole Y-chromosome sequence analysis of eight R1a and five R1b individuals suggests a divergence time of ∼25,000 (95% CI: 21,300-29,000) years ago and a coalescence time within R1a-M417 of ∼5800 (95% CI: 4800-6800) years. The spatial frequency distributions of R1a sub-haplogroups conclusively indicate two major groups, one found primarily in Europe and the other confined to Central and South Asia. Beyond the major European versus Asian dichotomy, we describe several younger sub-haplogroups. Based on spatial distributions and diversity patterns within the R1a-M420 clade, particularly rare basal branches detected primarily within Iran and eastern Turkey, we conclude that the initial episodes of haplogroup R1a diversification likely occurred in the vicinity of present-day Iran.

The impact of whole-genome sequencing on the reconstruction of human population history

Examining patterns of molecular genetic variation in both modern-day and ancient humans has proved to be a powerful approach to learn about our origins. Rapid advances in DNA sequencing technology have allowed us to characterize increasing amounts of genomic information. Although this clearly provides unprecedented power for inference, it also introduces more complexity into the way we use and interpret such data. Here, we review ongoing debates that have been influenced by improvements in our ability to sequence DNA and discuss some of the analytical challenges that need to be overcome in order to fully exploit the rich historical information that is contained in the entirety of the human genome.

Direct estimates of natural selection in Iberia indicate calcium absorption was not the only driver of lactase persistence in Europe

Lactase persistence (LP) is a genetically determined trait whereby the enzyme lactase is expressed throughout adult life. Lactase is necessary for the digestion of lactose--the main carbohydrate in milk--and its production is downregulated after the weaning period in most humans and all other mammals studied. Several sources of evidence indicate that LP has evolved independently, in different parts of the world over the last 10,000 years, and has been subject to strong natural selection in dairying populations. In Europeans, LP is strongly associated with, and probably caused by, a single C to T mutation 13,910 bp upstream of the lactase (LCT) gene (-13,910*T). Despite a considerable body of research, the reasons why LP should provide such a strong selective advantage remain poorly understood. In this study, we examine one of the most widely cited hypotheses for selection on LP--that fresh milk consumption supplemented the poor vitamin D and calcium status of northern Europe's early farmers (the calcium assimilation hypothesis). We do this by testing for natural selection on -13,910*T using ancient DNA data from the skeletal remains of eight late Neolithic Iberian individuals, whom we would not expect to have poor vitamin D and calcium status because of relatively high incident UVB light levels. None of the eight samples successfully typed in the study had the derived T-allele. In addition, we reanalyze published data from French Neolithic remains to both test for population continuity and further examine the evolution of LP in the region. Using simulations that accommodate genetic drift, natural selection, uncertainty in calibrated radiocarbon dates, and sampling error, we find that natural selection is still required to explain the observed increase in allele frequency. We conclude that the calcium assimilation hypothesis is insufficient to explain the spread of LP in Europe.

Neolithic mitochondrial haplogroup H genomes and the genetic origins of Europeans

Haplogroup H dominates present-day Western European mitochondrial DNA variability (>40%), yet was less common (~19%) among Early Neolithic farmers (~5450 BC) and virtually absent in Mesolithic hunter-gatherers. Here we investigate this major component of the maternal population history of modern Europeans and sequence 39 complete haplogroup H mitochondrial genomes from ancient human remains. We then compare this 'real-time' genetic data with cultural changes taking place between the Early Neolithic (~5450 BC) and Bronze Age (~2200 BC) in Central Europe. Our results reveal that the current diversity and distribution of haplogroup H were largely established by the Mid Neolithic (~4000 BC), but with substantial genetic contributions from subsequent pan-European cultures such as the Bell Beakers expanding out of Iberia in the Late Neolithic (~2800 BC). Dated haplogroup H genomes allow us to reconstruct the recent evolutionary history of haplogroup H and reveal a mutation rate 45% higher than current estimates for human mitochondria.

Ancient DNA reveals key stages in the formation of central European mitochondrial genetic diversity

The processes that shaped modern European mitochondrial DNA (mtDNA) variation remain unclear. The initial peopling by Palaeolithic hunter-gatherers ~42,000 years ago and the immigration of Neolithic farmers into Europe ~8000 years ago appear to have played important roles but do not explain present-day mtDNA diversity. We generated mtDNA profiles of 364 individuals from prehistoric cultures in Central Europe to perform a chronological study, spanning the Early Neolithic to the Early Bronze Age (5500 to 1550 calibrated years before the common era). We used this transect through time to identify four marked shifts in genetic composition during the Neolithic period, revealing a key role for Late Neolithic cultures in shaping modern Central European genetic diversity.