Northwest African Neolithic initiated by migrants from Iberia and Levant.

In northwestern Africa, lifestyle transitioned from foraging to food production around 7,400 years ago but what sparked that change remains unclear. Archaeological data support conflicting views: (1) that migrant European Neolithic farmers brought the new way of life to North Africa1-3 or (2) that local hunter-gatherers adopted technological innovations4,5. The latter view is also supported by archaeogenetic data6. Here we fill key chronological and archaeogenetic gaps for the Maghreb, from Epipalaeolithic to Middle Neolithic, by sequencing the genomes of nine individuals (to between 45.8- and 0.2-fold genome coverage). Notably, we trace 8,000 years of population continuity and isolation from the Upper Palaeolithic, via the Epipaleolithic, to some Maghrebi Neolithic farming groups. However, remains from the earliest Neolithic contexts showed mostly European Neolithic ancestry. We suggest that farming was introduced by European migrants and was then rapidly adopted by local groups. During the Middle Neolithic a new ancestry from the Levant appears in the Maghreb, coinciding with the arrival of pastoralism in the region, and all three ancestries blend together during the Late Neolithic. Our results show ancestry shifts in the Neolithization of northwestern Africa that probably mirrored a heterogeneous economic and cultural landscape, in a more multifaceted process than observed in other regions.

Genomic ancestry and social dynamics of the last hunter-gatherers of Atlantic France.

Since the early Holocene, western and central Europe was inhabited by a genetically distinct group of Western Hunter-Gatherers (WHGs). This group was eventually replaced and assimilated by the incoming Neolithic farmers. The western Atlantic façade was home to some of the last Mesolithic sites of mainland Europe, represented by the iconic open-air sites at Hoedic and Téviec in southern Brittany, France. These sites are known for the unusually well-preserved and rich burials. Genomic studies of Mesolithic European hunter-gatherers have been limited to single or a few individuals per site and our understanding of the social dynamics of the last Mesolithic hunter-gatherers of Europe and their interactions with incoming farmers is limited. We sequenced and analyzed the complete genomes of 10 individuals from the Late Mesolithic sites of Hoedic, Téviec, and Champigny, in France, four of which sequenced to between 23- and 8-times genome coverage. The analysis of genomic, chronological and dietary data revealed that the Late Mesolithic populations in Brittany maintained distinct social units within a network of exchanging mates. This resulted in low intra-group biological relatedness that prevented consanguineous mating, despite the small population size of the Late Mesolithic groups. We found no genetic ancestry from Neolithic farmers in the analyzed hunter-gatherers, even though some of them may have coexisted with the first farming groups in neighboring regions. Hence, contrary to previous conclusions based on stable isotope data from the same sites, the Late Mesolithic forager community was limited in mate-exchange to neighboring hunter-gatherer groups, to the exclusion of Neolithic farmers.

Imputed genomes and haplotype-based analyses of the Picts of early medieval Scotland reveal fine-scale relatedness between Iron Age, early medieval and the modern people of the UK.

There are longstanding questions about the origins and ancestry of the Picts of early medieval Scotland (ca. 300-900 CE), prompted in part by exotic medieval origin myths, their enigmatic symbols and inscriptions, and the meagre textual evidence. The Picts, first mentioned in the late 3rd century CE resisted the Romans and went on to form a powerful kingdom that ruled over a large territory in northern Britain. In the 9th and 10th centuries Gaelic language, culture and identity became dominant, transforming the Pictish realm into Alba, the precursor to the medieval kingdom of Scotland. To date, no comprehensive analysis of Pictish genomes has been published, and questions about their biological relationships to other cultural groups living in Britain remain unanswered. Here we present two high-quality Pictish genomes (2.4 and 16.5X coverage) from central and northern Scotland dated from the 5th-7th century which we impute and co-analyse with >8,300 previously published ancient and modern genomes. Using allele frequency and haplotype-based approaches, we can firmly place the genomes within the Iron Age gene pool in Britain and demonstrate regional biological affinity. We also demonstrate the presence of population structure within Pictish groups, with Orcadian Picts being genetically distinct from their mainland contemporaries. When investigating Identity-By-Descent (IBD) with present-day genomes, we observe broad affinities between the mainland Pictish genomes and the present-day people living in western Scotland, Wales, Northern Ireland and Northumbria, but less with the rest of England, the Orkney islands and eastern Scotland-where the political centres of Pictland were located. The pre-Viking Age Orcadian Picts evidence a high degree of IBD sharing across modern Scotland, Wales, Northern Ireland, and the Orkney islands, demonstrating substantial genetic continuity in Orkney for the last ~2,000 years. Analysis of mitochondrial DNA diversity at the Pictish cemetery of Lundin Links (n = 7) reveals absence of direct common female ancestors, with implications for broader social organisation. Overall, our study provides novel insights into the genetic affinities and population structure of the Picts and direct relationships between ancient and present-day groups of the UK.

The presence and impact of reference bias on population genomic studies of prehistoric human populations.

Haploid high quality reference genomes are an important resource in genomic research projects. A consequence is that DNA fragments carrying the reference allele will be more likely to map successfully, or receive higher quality scores. This reference bias can have effects on downstream population genomic analysis when heterozygous sites are falsely considered homozygous for the reference allele. In palaeogenomic studies of human populations, mapping against the human reference genome is used to identify endogenous human sequences. Ancient DNA studies usually operate with low sequencing coverages and fragmentation of DNA molecules causes a large proportion of the sequenced fragments to be shorter than 50 bp-reducing the amount of accepted mismatches, and increasing the probability of multiple matching sites in the genome. These ancient DNA specific properties are potentially exacerbating the impact of reference bias on downstream analyses, especially since most studies of ancient human populations use pseudo-haploid data, i.e. they randomly sample only one sequencing read per site. We show that reference bias is pervasive in published ancient DNA sequence data of prehistoric humans with some differences between individual genomic regions. We illustrate that the strength of reference bias is negatively correlated with fragment length. Most genomic regions we investigated show little to no mapping bias but even a small proportion of sites with bias can impact analyses of those particular loci or slightly skew genome-wide estimates. Therefore, reference bias has the potential to cause minor but significant differences in the results of downstream analyses such as population allele sharing, heterozygosity estimates and estimates of archaic ancestry. These spurious results highlight how important it is to be aware of these technical artifacts and that we need strategies to mitigate the effect. Therefore, we suggest some post-mapping filtering strategies to resolve reference bias which help to reduce its impact substantially.

Khoe-San Genomes Reveal Unique Variation and Confirm the Deepest Population Divergence in Homo sapiens.

The southern African indigenous Khoe-San populations harbor the most divergent lineages of all living peoples. Exploring their genomes is key to understanding deep human history. We sequenced 25 full genomes from five Khoe-San populations, revealing many novel variants, that 25% of variants are unique to the Khoe-San, and that the Khoe-San group harbors the greatest level of diversity across the globe. In line with previous studies, we found several gene regions with extreme values in genome-wide scans for selection, potentially caused by natural selection in the lineage leading to Homo sapiens and more recent in time. These gene regions included immunity-, sperm-, brain-, diet-, and muscle-related genes. When accounting for recent admixture, all Khoe-San groups display genetic diversity approaching the levels in other African groups and a reduction in effective population size starting around 100,000 years ago. Hence, all human groups show a reduction in effective population size commencing around the time of the Out-of-Africa migrations, which coincides with changes in the paleoclimate records, changes that potentially impacted all humans at the time.

Population genomics of Mesolithic Scandinavia: Investigating early postglacial migration routes and high-latitude adaptation.

Scandinavia was one of the last geographic areas in Europe to become habitable for humans after the Last Glacial Maximum (LGM). However, the routes and genetic composition of these postglacial migrants remain unclear. We sequenced the genomes, up to 57× coverage, of seven hunter-gatherers excavated across Scandinavia and dated from 9,500-6,000 years before present (BP). Surprisingly, among the Scandinavian Mesolithic individuals, the genetic data display an east-west genetic gradient that opposes the pattern seen in other parts of Mesolithic Europe. Our results suggest two different early postglacial migrations into Scandinavia: initially from the south, and later, from the northeast. The latter followed the ice-free Norwegian north Atlantic coast, along which novel and advanced pressure-blade stone-tool techniques may have spread. These two groups met and mixed in Scandinavia, creating a genetically diverse population, which shows patterns of genetic adaptation to high latitude environments. These potential adaptations include high frequencies of low pigmentation variants and a gene region associated with physical performance, which shows strong continuity into modern-day northern Europeans.

Four millennia of Iberian biomolecular prehistory illustrate the impact of prehistoric migrations at the far end of Eurasia.

Population genomic studies of ancient human remains have shown how modern-day European population structure has been shaped by a number of prehistoric migrations. The Neolithization of Europe has been associated with large-scale migrations from Anatolia, which was followed by migrations of herders from the Pontic steppe at the onset of the Bronze Age. Southwestern Europe was one of the last parts of the continent reached by these migrations, and modern-day populations from this region show intriguing similarities to the initial Neolithic migrants. Partly due to climatic conditions that are unfavorable for DNA preservation, regional studies on the Mediterranean remain challenging. Here, we present genome-wide sequence data from 13 individuals combined with stable isotope analysis from the north and south of Iberia covering a four-millennial temporal transect (7,500-3,500 BP). Early Iberian farmers and Early Central European farmers exhibit significant genetic differences, suggesting two independent fronts of the Neolithic expansion. The first Neolithic migrants that arrived in Iberia had low levels of genetic diversity, potentially reflecting a small number of individuals; this diversity gradually increased over time from mixing with local hunter-gatherers and potential population expansion. The impact of post-Neolithic migrations on Iberia was much smaller than for the rest of the continent, showing little external influence from the Neolithic to the Bronze Age. Paleodietary reconstruction shows that these populations have a remarkable degree of dietary homogeneity across space and time, suggesting a strong reliance on terrestrial food resources despite changing culture and genetic make-up.

The Neolithic Pitted Ware culture foragers were culturally but not genetically influenced by the Battle Axe culture herders.

OBJECTIVES: In order to understand contacts between cultural spheres in the third millennium BC, we investigated the impact of a new herder culture, the Battle Axe culture, arriving to Scandinavia on the people of the sub-Neolithic hunter-gatherer Pitted Ware culture. By investigating the genetic make-up of Pitted Ware culture people from two types of burials (typical Pitted Ware culture burials and Battle Axe culture-influenced burials), we could determine the impact of migration and the impact of cultural influences. METHODS: We sequenced and analyzed the genomes of 25 individuals from typical Pitted Ware culture burials and from Pitted Ware culture burials with Battle Axe culture influences in order to determine if the different burial types were associated with different gene-pools. RESULTS: The genomic data show that all individuals belonged to one genetic population-a population associated with the Pitted Ware culture-irrespective of the burial style. CONCLUSION: We conclude that the Pitted Ware culture communities were not impacted by gene-flow, that is, via migration or exchange of mates. These different cultural expressions in the Pitted Ware culture burials are instead a consequence of cultural exchange.

Ancient X chromosomes reveal contrasting sex bias in Neolithic and Bronze Age Eurasian migrations.

Dramatic events in human prehistory, such as the spread of agriculture to Europe from Anatolia and the late Neolithic/Bronze Age migration from the Pontic-Caspian Steppe, can be investigated using patterns of genetic variation among the people who lived in those times. In particular, studies of differing female and male demographic histories on the basis of ancient genomes can provide information about complexities of social structures and cultural interactions in prehistoric populations. We use a mechanistic admixture model to compare the sex-specifically-inherited X chromosome with the autosomes in 20 early Neolithic and 16 late Neolithic/Bronze Age human remains. Contrary to previous hypotheses suggested by the patrilocality of many agricultural populations, we find no evidence of sex-biased admixture during the migration that spread farming across Europe during the early Neolithic. For later migrations from the Pontic Steppe during the late Neolithic/Bronze Age, however, we estimate a dramatic male bias, with approximately five to 14 migrating males for every migrating female. We find evidence of ongoing, primarily male, migration from the steppe to central Europe over a period of multiple generations, with a level of sex bias that excludes a pulse migration during a single generation. The contrasting patterns of sex-specific migration during these two migrations suggest a view of differing cultural histories in which the Neolithic transition was driven by mass migration of both males and females in roughly equal numbers, perhaps whole families, whereas the later Bronze Age migration and cultural shift were instead driven by male migration, potentially connected to new technology and conquest.

Northeast African genomic variation shaped by the continuity of indigenous groups and Eurasian migrations.

Northeast Africa has a long history of human habitation, with fossil-finds from the earliest anatomically modern humans, and housing ancient civilizations. The region is also the gate-way out of Africa, as well as a portal for migration into Africa from Eurasia via the Middle East and the Arabian Peninsula. We investigate the population history of northeast Africa by genotyping ~3.9 million SNPs in 221 individuals from 18 populations sampled in Sudan and South Sudan and combine this data with published genome-wide data from surrounding areas. We find a strong genetic divide between the populations from the northeastern parts of the region (Nubians, central Arab populations, and the Beja) and populations towards the west and south (Nilotes, Darfur and Kordofan populations). This differentiation is mainly caused by a large Eurasian ancestry component of the northeast populations likely driven by migration of Middle Eastern groups followed by admixture that affected the local populations in a north-to-south succession of events. Genetic evidence points to an early admixture event in the Nubians, concurrent with historical contact between North Sudanese and Arab groups. We estimate the admixture in current-day Sudanese Arab populations to about 700 years ago, coinciding with the fall of Dongola in 1315/1316 AD, a wave of admixture that reached the Darfurian/Kordofanian populations some 400-200 years ago. In contrast to the northeastern populations, the current-day Nilotic populations from the south of the region display little or no admixture from Eurasian groups indicating long-term isolation and population continuity in these areas of northeast Africa.

Copy number determination of the gene for the human pancreatic polypeptide receptor NPY4R using read depth analysis and droplet digital PCR.

BACKGROUND: Copy number variation (CNV) plays an important role in human genetic diversity and has been associated with multiple complex disorders. Here we investigate a CNV on chromosome 10q11.22 that spans NPY4R, the gene for the appetite-regulating pancreatic polypeptide receptor Y4. This genomic region has been challenging to map due to multiple repeated elements and its precise organization has not yet been resolved. Previous studies using microarrays were interpreted to show that the most common copy number was 2 per genome. RESULTS: We have investigated 18 individuals from the 1000 Genomes project using the well-established method of read depth analysis and the new droplet digital PCR (ddPCR) method. We find that the most common copy number for NPY4R is 4. The estimated number of copies ranged from three to seven based on read depth analyses with Control-FREEC and CNVnator, and from four to seven based on ddPCR. We suggest that the difference between our results and those published previously can be explained by methodological differences such as reference gene choice, data normalization and method reliability. Three high-quality archaic human genomes (two Neanderthal and one Denisova) display four copies of the NPY4R gene indicating that a duplication occurred prior to the human-Neanderthal/Denisova split. CONCLUSIONS: We conclude that ddPCR is a sensitive and reliable method for CNV determination, that it can be used for read depth calibration in CNV studies based on already available whole-genome sequencing data, and that further investigation of NPY4R copy number variation and its consequences are necessary due to the role of Y4 receptor in food intake regulation.

The genomic ancestry of the Scandinavian Battle Axe Culture people and their relation to the broader Corded Ware horizon.

The Neolithic period is characterized by major cultural transformations and human migrations, with lasting effects across Europe. To understand the population dynamics in Neolithic Scandinavia and the Baltic Sea area, we investigate the genomes of individuals associated with the Battle Axe Culture (BAC), a Middle Neolithic complex in Scandinavia resembling the continental Corded Ware Culture (CWC). We sequenced 11 individuals (dated to 3330-1665 calibrated before common era (cal BCE)) from modern-day Sweden, Estonia, and Poland to 0.26-3.24× coverage. Three of the individuals were from CWC contexts and two from the central-Swedish BAC burial 'Bergsgraven'. By analysing these genomes together with the previously published data, we show that the BAC represents a group different from other Neolithic populations in Scandinavia, revealing stratification among cultural groups. Similar to continental CWC, the BAC-associated individuals display ancestry from the Pontic-Caspian steppe herders, as well as smaller components originating from hunter-gatherers and Early Neolithic farmers. Thus, the steppe ancestry seen in these Scandinavian BAC individuals can be explained only by migration into Scandinavia. Furthermore, we highlight the reuse of megalithic tombs of the earlier Funnel Beaker Culture (FBC) by people related to BAC. The BAC groups likely mixed with resident middle Neolithic farmers (e.g. FBC) without substantial contributions from Neolithic foragers.

Ancient genomes link early farmers from Atapuerca in Spain to modern-day Basques.

The consequences of the Neolithic transition in Europe--one of the most important cultural changes in human prehistory--is a subject of great interest. However, its effect on prehistoric and modern-day people in Iberia, the westernmost frontier of the European continent, remains unresolved. We present, to our knowledge, the first genome-wide sequence data from eight human remains, dated to between 5,500 and 3,500 years before present, excavated in the El Portalón cave at Sierra de Atapuerca, Spain. We show that these individuals emerged from the same ancestral gene pool as early farmers in other parts of Europe, suggesting that migration was the dominant mode of transferring farming practices throughout western Eurasia. In contrast to central and northern early European farmers, the Chalcolithic El Portalón individuals additionally mixed with local southwestern hunter-gatherers. The proportion of hunter-gatherer-related admixture into early farmers also increased over the course of two millennia. The Chalcolithic El Portalón individuals showed greatest genetic affinity to modern-day Basques, who have long been considered linguistic and genetic isolates linked to the Mesolithic whereas all other European early farmers show greater genetic similarity to modern-day Sardinians. These genetic links suggest that Basques and their language may be linked with the spread of agriculture during the Neolithic. Furthermore, all modern-day Iberian groups except the Basques display distinct admixture with Caucasus/Central Asian and North African groups, possibly related to historical migration events. The El Portalón genomes uncover important pieces of the demographic history of Iberia and Europe and reveal how prehistoric groups relate to modern-day people.

A female Viking warrior confirmed by genomics.

OBJECTIVES: The objective of this study has been to confirm the sex and the affinity of an individual buried in a well-furnished warrior grave (Bj 581) in the Viking Age town of Birka, Sweden. Previously, based on the material and historical records, the male sex has been associated with the gender of the warrior and such was the case with Bj 581. An earlier osteological classification of the individual as female was considered controversial in a historical and archaeological context. A genomic confirmation of the biological sex of the individual was considered necessary to solve the issue. MATERIALS AND METHODS: Genome-wide sequence data was generated in order to confirm the biological sex, to support skeletal integrity, and to investigate the genetic relationship of the individual to ancient individuals as well as modern-day groups. Additionally, a strontium isotope analysis was conducted to highlight the mobility of the individual. RESULTS: The genomic results revealed the lack of a Y-chromosome and thus a female biological sex, and the mtDNA analyses support a single-individual origin of sampled elements. The genetic affinity is close to present-day North Europeans, and within Sweden to the southern and south-central region. Nevertheless, the Sr values are not conclusive as to whether she was of local or nonlocal origin. DISCUSSION: The identification of a female Viking warrior provides a unique insight into the Viking society, social constructions, and exceptions to the norm in the Viking time-period. The results call for caution against generalizations regarding social orders in past societies.

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps.

Altitudinal gradients in mountain regions are short-range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19-29 individuals from each population. The sample includes two pairs of low- and high-altitude populations from two different valleys. High-altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low-altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low- and high-altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low- and high-altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low- and high-altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.

Mutational bias and gene conversion affect the intraspecific nitrogen stoichiometry of the Arabidopsis thaliana transcriptome.

The transcriptome and proteome of Arabidopsis thaliana are reduced in nitrogen content when compared with other taxa, which may result from ecological nitrogen limitation. We hypothesized that if the A. thaliana transcriptome is selected for a low nitrogen content, nitrogen-reducing derived alleles of single nucleotide polymorphisms (SNPs) should segregate at higher frequencies than nitrogen-increasing alleles. This pattern should be stronger in populations with a larger effective population size (N(e)) if natural selection is more efficient in large than in small populations. We analyzed variation in the nitrogen content in the transcriptome of 80 natural accessions of A. thaliana. In contrast to our expectations, derived alleles increase the nitrogen content in all accessions, and there is a positive correlation between nitrogen difference and derived allele frequency, which is strongest with nonsynonymous SNPs (nsSNPs). Also, there is a positive correlation between nitrogen difference and N(e) that was mainly caused by nsSNPs. These observations led us to reject the hypothesis that the transcriptome of A. thaliana is currently under selection to reduce nitrogen content. Instead, we show that a change in nitrogen content is a side effect of interacting evolutionary factors that influence base composition and include mutational bias, purifying selection of functionally deleterious alleles, and GC-biased gene conversion. We provide strong evidence that GC-biased gene conversion may play an important role for base composition in the highly selfing plant A. thaliana.

Ancient Sheep Genomes Reveal Four Millennia of North European Short-Tailed Sheep in the Baltic Sea Region.

Sheep are among the earliest domesticated livestock species, with a wide variety of breeds present today. However, it remains unclear how far back this diversity goes, with formal documentation only dating back a few centuries. North European short-tailed (NEST) breeds are often assumed to be among the oldest domestic sheep populations, even thought to represent relicts of the earliest sheep expansions during the Neolithic period reaching Scandinavia <6,000 years ago. This study sequenced the genomes (up to 11.6X) of five sheep remains from the Baltic islands of Gotland and Åland, dating from the Late Neolithic (∼4,100 cal BP) to historical times (∼1,600 CE). Our findings indicate that these ancient sheep largely possessed the genetic characteristics of modern NEST breeds, suggesting a substantial degree of long-term continuity of this sheep type in the Baltic Sea region. Despite the wide temporal spread, population genetic analyses show high levels of affinity between the ancient genomes and they also exhibit relatively high genetic diversity when compared to modern NEST breeds, implying a loss of diversity in most breeds during the last centuries associated with breed formation and recent bottlenecks. Our results shed light on the development of breeds in Northern Europe specifically as well as the development of genetic diversity in sheep breeds, and their expansion from the domestication center in general.

Population Genomic History of the Endangered Anatolian and Cyprian Mouflons in Relation to Worldwide Wild, Feral, and Domestic Sheep Lineages.

Once widespread in their homelands, the Anatolian mouflon (Ovis gmelini anatolica) and the Cyprian mouflon (Ovis gmelini ophion) were driven to near extinction during the 20th century and are currently listed as endangered populations by the International Union for Conservation of Nature. While the exact origins of these lineages remain unclear, they have been suggested to be close relatives of domestic sheep or remnants of proto-domestic sheep. Here, we study whole genome sequences of n = 5 Anatolian mouflons and n = 10 Cyprian mouflons in terms of population history and diversity, comparing them with eight other extant sheep lineages. We find reciprocal genetic affinity between Anatolian and Cyprian mouflons and domestic sheep, higher than all other studied wild sheep genomes, including the Iranian mouflon (O. gmelini). Studying diversity indices, we detect a considerable load of short runs of homozygosity blocks (<2 Mb) in both Anatolian and Cyprian mouflons, reflecting small effective population size (Ne). Meanwhile, Ne and mutation load estimates are lower in Cyprian compared with Anatolian mouflons, suggesting the purging of recessive deleterious variants in Cyprian sheep under a small long-term Ne, possibly attributable to founder effects, island isolation, introgression from domestic lineages, or differences in their bottleneck dynamics. Expanding our analyses to worldwide wild and feral Ovis genomes, we observe varying viability metrics among different lineages and a limited consistency between viability metrics and International Union for Conservation of Nature conservation status. Factors such as recent inbreeding, introgression, and unique population dynamics may have contributed to the observed disparities.

Demographic reconstruction of the Western sheep expansion from whole-genome sequences.

As one of the earliest livestock, sheep (Ovis aries) were domesticated in the Fertile Crescent about 12,000-10,000 years ago and have a nearly worldwide distribution today. Most of our knowledge about the timing of their expansions stems from archaeological data but it is unclear how the genetic diversity of modern sheep fits with these dates. We used whole-genome sequencing data of 63 domestic breeds and their wild relatives, the Asiatic mouflon (O. gmelini, previously known as O. orientalis), to explore the demographic history of sheep. On the global scale, our analysis revealed geographic structuring among breeds with unidirectional recent gene flow from domestics into Asiatic mouflons. We then selected 4 representative breeds from Spain, Morocco, the United Kingdom, and Iran to build a comprehensive demographic model of the Western sheep expansion. We inferred a single domestication event around 11,000 years ago. The subsequent westward expansion is dated to approximately 7,000 years ago, later than the original Neolithic expansion of sheep and slightly predating the Secondary Product Revolution associated with wooly sheep. We see some signals of recent gene flow from an ancestral population into Southern European breeds which could reflect admixture with feral European mouflon. Furthermore, our results indicate that many breeds experienced a reduction of their effective population size during the last centuries, probably associated with modern breed development. Our study provides insights into the complex demographic history of Western Eurasian sheep, highlighting interactions between breeds and their wild counterparts.