Widespread horse-based mobility arose around 2200 BCE in Eurasia.

Horses revolutionized human history with fast mobility1. However, the timeline between their domestication and their widespread integration as a means of transport remains contentious2-4. Here we assemble a collection of 475 ancient horse genomes to assess the period when these animals were first reshaped by human agency in Eurasia. We find that reproductive control of the modern domestic lineage emerged around 2200 BCE, through close-kin mating and shortened generation times. Reproductive control emerged following a severe domestication bottleneck starting no earlier than approximately 2700 BCE, and coincided with a sudden expansion across Eurasia that ultimately resulted in the replacement of nearly every local horse lineage. This expansion marked the rise of widespread horse-based mobility in human history, which refutes the commonly held narrative of large horse herds accompanying the massive migration of steppe peoples across Europe around 3000 BCE and earlier3,5. Finally, we detect significantly shortened generation times at Botai around 3500 BCE, a settlement from central Asia associated with corrals and a subsistence economy centred on horses6,7. This supports local horse husbandry before the rise of modern domestic bloodlines.

Animal domestication in the era of ancient genomics.

The domestication of animals led to a major shift in human subsistence patterns, from a hunter-gatherer to a sedentary agricultural lifestyle, which ultimately resulted in the development of complex societies. Over the past 15,000 years, the phenotype and genotype of multiple animal species, such as dogs, pigs, sheep, goats, cattle and horses, have been substantially altered during their adaptation to the human niche. Recent methodological innovations, such as improved ancient DNA extraction methods and next-generation sequencing, have enabled the sequencing of whole ancient genomes. These genomes have helped reconstruct the process by which animals entered into domestic relationships with humans and were subjected to novel selection pressures. Here, we discuss and update key concepts in animal domestication in light of recent contributions from ancient genomics.

A dynastic elite in monumental Neolithic society.

The nature and distribution of political power in Europe during the Neolithic era remains poorly understood1. During this period, many societies began to invest heavily in building monuments, which suggests an increase in social organization. The scale and sophistication of megalithic architecture along the Atlantic seaboard, culminating in the great passage tomb complexes, is particularly impressive2. Although co-operative ideology has often been emphasised as a driver of megalith construction1, the human expenditure required to erect the largest monuments has led some researchers to emphasize hierarchy3-of which the most extreme case is a small elite marshalling the labour of the masses. Here we present evidence that a social stratum of this type was established during the Neolithic period in Ireland. We sampled 44 whole genomes, among which we identify the adult son of a first-degree incestuous union from remains that were discovered within the most elaborate recess of the Newgrange passage tomb. Socially sanctioned matings of this nature are very rare, and are documented almost exclusively among politico-religious elites4-specifically within polygynous and patrilineal royal families that are headed by god-kings5,6. We identify relatives of this individual within two other major complexes of passage tombs 150 km to the west of Newgrange, as well as dietary differences and fine-scale haplotypic structure (which is unprecedented in resolution for a prehistoric population) between passage tomb samples and the larger dataset, which together imply hierarchy. This elite emerged against a backdrop of rapid maritime colonization that displaced a unique Mesolithic isolate population, although we also detected rare Irish hunter-gatherer introgression within the Neolithic population.

Ancient DNA at the edge of the world: Continental immigration and the persistence of Neolithic male lineages in Bronze Age Orkney.

Orkney was a major cultural center during the Neolithic, 3800 to 2500 BC. Farming flourished, permanent stone settlements and chambered tombs were constructed, and long-range contacts were sustained. From ∼3200 BC, the number, density, and extravagance of settlements increased, and new ceremonial monuments and ceramic styles, possibly originating in Orkney, spread across Britain and Ireland. By ∼2800 BC, this phenomenon was waning, although Neolithic traditions persisted to at least 2500 BC. Unlike elsewhere in Britain, there is little material evidence to suggest a Beaker presence, suggesting that Orkney may have developed along an insular trajectory during the second millennium BC. We tested this by comparing new genomic evidence from 22 Bronze Age and 3 Iron Age burials in northwest Orkney with Neolithic burials from across the archipelago. We identified signals of inward migration on a scale unsuspected from the archaeological record: As elsewhere in Bronze Age Britain, much of the population displayed significant genome-wide ancestry deriving ultimately from the Pontic-Caspian Steppe. However, uniquely in northern and central Europe, most of the male lineages were inherited from the local Neolithic. This suggests that some male descendants of Neolithic Orkney may have remained distinct well into the Bronze Age, although there are signs that this had dwindled by the Iron Age. Furthermore, although the majority of mitochondrial DNA lineages evidently arrived afresh with the Bronze Age, we also find evidence for continuity in the female line of descent from Mesolithic Britain into the Bronze Age and even to the present day.

Herded and hunted goat genomes from the dawn of domestication in the Zagros Mountains.

The Aceramic Neolithic (∼9600 to 7000 cal BC) period in the Zagros Mountains, western Iran, provides some of the earliest archaeological evidence of goat (Capra hircus) management and husbandry by circa 8200 cal BC, with detectable morphological change appearing ∼1,000 y later. To examine the genomic imprint of initial management and its implications for the goat domestication process, we analyzed 14 novel nuclear genomes (mean coverage 1.13X) and 32 mitochondrial (mtDNA) genomes (mean coverage 143X) from two such sites, Ganj Dareh and Tepe Abdul Hosein. These genomes show two distinct clusters: those with domestic affinity and a minority group with stronger wild affinity, indicating that managed goats were genetically distinct from wild goats at this early horizon. This genetic duality, the presence of long runs of homozygosity, shared ancestry with later Neolithic populations, a sex bias in archaeozoological remains, and demographic profiles from across all layers of Ganj Dareh support management of genetically domestic goat by circa 8200 cal BC, and represent the oldest to-this-date reported livestock genomes. In these sites a combination of high autosomal and mtDNA diversity, contrasting limited Y chromosomal lineage diversity, an absence of reported selection signatures for pigmentation, and the wild morphology of bone remains illustrates domestication as an extended process lacking a strong initial bottleneck, beginning with spatial control, demographic manipulation via biased male culling, captive breeding, and subsequently phenotypic and genomic selection.

Modern Siberian dog ancestry was shaped by several thousand years of Eurasian-wide trade and human dispersal.

Dogs have been essential to life in the Siberian Arctic for over 9,500 y, and this tight link between people and dogs continues in Siberian communities. Although Arctic Siberian groups such as the Nenets received limited gene flow from neighboring groups, archaeological evidence suggests that metallurgy and new subsistence strategies emerged in Northwest Siberia around 2,000 y ago. It is unclear if the Siberian Arctic dog population was as continuous as the people of the region or if instead admixture occurred, possibly in relation to the influx of material culture from other parts of Eurasia. To address this question, we sequenced and analyzed the genomes of 20 ancient and historical Siberian and Eurasian Steppe dogs. Our analyses indicate that while Siberian dogs were genetically homogenous between 9,500 to 7,000 y ago, later introduction of dogs from the Eurasian Steppe and Europe led to substantial admixture. This is clearly the case in the Iamal-Nenets region (Northwestern Siberia) where dogs from the Iron Age period (∼2,000 y ago) possess substantially less ancestry related to European and Steppe dogs than dogs from the medieval period (∼1,000 y ago). Combined with findings of nonlocal materials recovered from these archaeological sites, including glass beads and metal items, these results indicate that Northwest Siberian communities were connected to a larger trade network through which they acquired genetically distinctive dogs from other regions. These exchanges were part of a series of major societal changes, including the rise of large-scale reindeer pastoralism ∼800 y ago.

Exploring the phylogeography and population dynamics of the giant deer (Megaloceros giganteus) using Late Quaternary mitogenomes.

Late Quaternary climatic fluctuations in the Northern Hemisphere had drastic effects on large mammal species, leading to the extinction of a substantial number of them. The giant deer (Megaloceros giganteus) was one of the species that became extinct in the Holocene, around 7660 calendar years before present. In the Late Pleistocene, the species ranged from western Europe to central Asia. However, during the Holocene, its range contracted to eastern Europe and western Siberia, where the last populations of the species occurred. Here, we generated 35 Late Pleistocene and Holocene giant deer mitogenomes to explore the genetics of the demise of this iconic species. Bayesian phylogenetic analyses of the mitogenomes suggested five main clades for the species: three pre-Last Glacial Maximum clades that did not appear in the post-Last Glacial Maximum genetic pool, and two clades that showed continuity into the Holocene. Our study also identified a decrease in genetic diversity starting in Marine Isotope Stage 3 and accelerating during the Last Glacial Maximum. This reduction in genetic diversity during the Last Glacial Maximum, coupled with a major contraction of fossil occurrences, suggests that climate was a major driver in the dynamics of the giant deer.

Insular Celtic population structure and genomic footprints of migration.

Previous studies of the genetic landscape of Ireland have suggested homogeneity, with population substructure undetectable using single-marker methods. Here we have harnessed the haplotype-based method fineSTRUCTURE in an Irish genome-wide SNP dataset, identifying 23 discrete genetic clusters which segregate with geographical provenance. Cluster diversity is pronounced in the west of Ireland but reduced in the east where older structure has been eroded by historical migrations. Accordingly, when populations from the neighbouring island of Britain are included, a west-east cline of Celtic-British ancestry is revealed along with a particularly striking correlation between haplotypes and geography across both islands. A strong relationship is revealed between subsets of Northern Irish and Scottish populations, where discordant genetic and geographic affinities reflect major migrations in recent centuries. Additionally, Irish genetic proximity of all Scottish samples likely reflects older strata of communication across the narrowest inter-island crossing. Using GLOBETROTTER we detected Irish admixture signals from Britain and Europe and estimated dates for events consistent with the historical migrations of the Norse-Vikings, the Anglo-Normans and the British Plantations. The influence of the former is greater than previously estimated from Y chromosome haplotypes. In all, we paint a new picture of the genetic landscape of Ireland, revealing structure which should be considered in the design of studies examining rare genetic variation and its association with traits.

Origins and genetic legacies of the Caribbean Taino.

The Caribbean was one of the last parts of the Americas to be settled by humans, but how and when the islands were first occupied remains a matter of debate. Ancient DNA can help answering these questions, but the work has been hampered by poor DNA preservation. We report the genome sequence of a 1,000-year-old Lucayan Taino individual recovered from the site of Preacher's Cave in the Bahamas. We sequenced her genome to 12.4-fold coverage and show that she is genetically most closely related to present-day Arawakan speakers from northern South America, suggesting that the ancestors of the Lucayans originated there. Further, we find no evidence for recent inbreeding or isolation in the ancient genome, suggesting that the Lucayans had a relatively large effective population size. Finally, we show that the native American components in some present-day Caribbean genomes are closely related to the ancient Taino, demonstrating an element of continuity between precontact populations and present-day Latino populations in the Caribbean.

The population genomics of archaeological transition in west Iberia: Investigation of ancient substructure using imputation and haplotype-based methods.

We analyse new genomic data (0.05-2.95x) from 14 ancient individuals from Portugal distributed from the Middle Neolithic (4200-3500 BC) to the Middle Bronze Age (1740-1430 BC) and impute genomewide diploid genotypes in these together with published ancient Eurasians. While discontinuity is evident in the transition to agriculture across the region, sensitive haplotype-based analyses suggest a significant degree of local hunter-gatherer contribution to later Iberian Neolithic populations. A more subtle genetic influx is also apparent in the Bronze Age, detectable from analyses including haplotype sharing with both ancient and modern genomes, D-statistics and Y-chromosome lineages. However, the limited nature of this introgression contrasts with the major Steppe migration turnovers within third Millennium northern Europe and echoes the survival of non-Indo-European language in Iberia. Changes in genomic estimates of individual height across Europe are also associated with these major cultural transitions, and ancestral components continue to correlate with modern differences in stature.

Neolithic and Bronze Age migration to Ireland and establishment of the insular Atlantic genome.

The Neolithic and Bronze Age transitions were profound cultural shifts catalyzed in parts of Europe by migrations, first of early farmers from the Near East and then Bronze Age herders from the Pontic Steppe. However, a decades-long, unresolved controversy is whether population change or cultural adoption occurred at the Atlantic edge, within the British Isles. We address this issue by using the first whole genome data from prehistoric Irish individuals. A Neolithic woman (3343-3020 cal BC) from a megalithic burial (10.3× coverage) possessed a genome of predominantly Near Eastern origin. She had some hunter-gatherer ancestry but belonged to a population of large effective size, suggesting a substantial influx of early farmers to the island. Three Bronze Age individuals from Rathlin Island (2026-1534 cal BC), including one high coverage (10.5×) genome, showed substantial Steppe genetic heritage indicating that the European population upheavals of the third millennium manifested all of the way from southern Siberia to the western ocean. This turnover invites the possibility of accompanying introduction of Indo-European, perhaps early Celtic, language. Irish Bronze Age haplotypic similarity is strongest within modern Irish, Scottish, and Welsh populations, and several important genetic variants that today show maximal or very high frequencies in Ireland appear at this horizon. These include those coding for lactase persistence, blue eye color, Y chromosome R1b haplotypes, and the hemochromatosis C282Y allele; to our knowledge, the first detection of a known Mendelian disease variant in prehistory. These findings together suggest the establishment of central attributes of the Irish genome 4,000 y ago.

Early farmers from across Europe directly descended from Neolithic Aegeans.

Farming and sedentism first appeared in southwestern Asia during the early Holocene and later spread to neighboring regions, including Europe, along multiple dispersal routes. Conspicuous uncertainties remain about the relative roles of migration, cultural diffusion, and admixture with local foragers in the early Neolithization of Europe. Here we present paleogenomic data for five Neolithic individuals from northern Greece and northwestern Turkey spanning the time and region of the earliest spread of farming into Europe. We use a novel approach to recalibrate raw reads and call genotypes from ancient DNA and observe striking genetic similarity both among Aegean early farmers and with those from across Europe. Our study demonstrates a direct genetic link between Mediterranean and Central European early farmers and those of Greece and Anatolia, extending the European Neolithic migratory chain all the way back to southwestern Asia.

Animal origin of 13th-century uterine vellum revealed using noninvasive peptide fingerprinting.

Tissue-thin parchment made it possible to produce the first pocket Bibles: Thousands were made in the 13th century. The source of this parchment, often called "uterine vellum," has been a long-standing controversy in codicology. Use of the Latin term abortivum in many sources has led some scholars to suggest that the skin of fetal calves or sheep was used. Others have argued that it would not be possible to sustain herds if so many pocket Bibles were produced from fetal skins, arguing instead for unexpected alternatives, such as rabbit. Here, we report a simple and objective technique using standard conservation treatments to identify the animal origin of parchment. The noninvasive method is a variant on zooarchaeology by mass spectrometry (ZooMS) peptide mass fingerprinting but extracts protein from the parchment surface by using an electrostatic charge generated by gentle rubbing of a PVC eraser on the membrane surface. Using this method, we analyzed 72 pocket Bibles originating in France, England, and Italy and 293 additional parchment samples that bracket this period. We found no evidence for the use of unexpected animals; however, we did identify the use of more than one mammal species in a single manuscript, consistent with the local availability of hides. These results suggest that ultrafine vellum does not necessarily derive from the use of abortive or newborn animals with ultrathin hides, but could equally well reflect a production process that allowed the skins of maturing animals of several species to be rendered into vellum of equal quality and fineness.

Capturing goats: documenting two hundred years of mitochondrial DNA diversity among goat populations from Britain and Ireland.

The domestic goat (Capra hircus) plays a key role in global agriculture, being especially prized in regions of marginal pasture. However, the advent of industrialized breeding has seen a dramatic reduction in genetic diversity within commercial populations, while high extinction rates among feral herds have further depleted the reservoir of genetic variation available. Here, we present the first survey of whole mitochondrial genomic variation among the modern and historical goat populations of Britain and Ireland using a combination of mtDNA enrichment and high throughput sequencing. Fifteen historical taxidermy samples, representing the indigenous 'Old Goat' populations of the islands, were sequenced alongside five modern Irish dairy goats and four feral samples from endangered populations in western Ireland. Phylogenetic and network analyses of European mitochondrial variation revealed distinct groupings dominated by historical British and Irish samples, which demonstrate a degree of maternal genetic structure between the goats of insular and continental Europe. Several Irish modern feral samples also fall within these clusters, suggesting continuity between these dwindling populations and the ancestral 'Old Goats' of Ireland and Britain.

Storytelling and story testing in domestication.

The domestication of plants and animals marks one of the most significant transitions in human, and indeed global, history. Traditionally, study of the domestication process was the exclusive domain of archaeologists and agricultural scientists; today it is an increasingly multidisciplinary enterprise that has come to involve the skills of evolutionary biologists and geneticists. Although the application of new information sources and methodologies has dramatically transformed our ability to study and understand domestication, it has also generated increasingly large and complex datasets, the interpretation of which is not straightforward. In particular, challenges of equifinality, evolutionary variance, and emergence of unexpected or counter-intuitive patterns all face researchers attempting to infer past processes directly from patterns in data. We argue that explicit modeling approaches, drawing upon emerging methodologies in statistics and population genetics, provide a powerful means of addressing these limitations. Modeling also offers an approach to analyzing datasets that avoids conclusions steered by implicit biases, and makes possible the formal integration of different data types. Here we outline some of the modeling approaches most relevant to current problems in domestication research, and demonstrate the ways in which simulation modeling is beginning to reshape our understanding of the domestication process.

Definition of the cattle killer cell Ig-like receptor gene family: comparison with aurochs and human counterparts.

Under selection pressure from pathogens, variable NK cell receptors that recognize polymorphic MHC class I evolved convergently in different species of placental mammal. Unexpectedly, diversified killer cell Ig-like receptors (KIRs) are shared by simian primates, including humans, and cattle, but not by other species. Whereas much is known of human KIR genetics and genomics, knowledge of cattle KIR is limited to nine cDNA sequences. To facilitate comparison of the cattle and human KIR gene families, we determined the genomic location, structure, and sequence of two cattle KIR haplotypes and defined KIR sequences of aurochs, the extinct wild ancestor of domestic cattle. Larger than its human counterpart, the cattle KIR locus evolved through successive duplications of a block containing ancestral KIR3DL and KIR3DX genes that existed before placental mammals. Comparison of two cattle KIR haplotypes and aurochs KIR show the KIR are polymorphic and the gene organization and content appear conserved. Of 18 genes, 8 are functional and 10 were inactivated by point mutation. Selective inactivation of KIR3DL and activating receptor genes leaves a functional cohort of one inhibitory KIR3DL, one activating KIR3DX, and six inhibitory KIR3DX. Functional KIR diversity evolved from KIR3DX in cattle and from KIR3DL in simian primates. Although independently evolved, cattle and human KIR gene families share important function-related properties, indicating that cattle KIR are NK cell receptors for cattle MHC class I. Combinations of KIR and MHC class I are the major genetic factors associated with human disease and merit investigation in cattle.

A genome-wide association study for genetic susceptibility to Mycobacterium bovis infection in dairy cattle identifies a susceptibility QTL on chromosome 23.

BACKGROUND: Bovine tuberculosis (bTB) infection in cattle is a significant economic concern in many countries, with annual costs to the UK and Irish governments of approximately €190 million and €63 million, respectively, for bTB control. The existence of host additive and non-additive genetic components to bTB susceptibility has been established. METHODS: Two approaches i.e. single-SNP (single nucleotide polymorphism) regression and a Bayesian method were applied to genome-wide association studies (GWAS) using high-density SNP genotypes (n = 597,144 SNPs) from 841 dairy artificial insemination (AI) sires. Deregressed estimated breeding values for bTB susceptibility were used as the quantitative dependent variable. Network analysis was performed using the quantitative trait loci (QTL) that were identified as significant in the single-SNP regression and Bayesian analyses separately. In addition, an identity-by-descent analysis was performed on a subset of the most prolific sires in the dataset that showed contrasting prevalences of bTB infection in daughters. RESULTS: A significant QTL region was identified on BTA23 (P value >1 × 10(-5), Bayes factor >10) across all analyses. Sires with the minor allele (minor allele frequency = 0.136) for this QTL on BTA23 had estimated breeding values that conferred a greater susceptibility to bTB infection than those that were homozygous for the major allele. Imputation of the regions that flank this QTL on BTA23 to full sequence indicated that the most significant associations were located within introns of the FKBP5 gene. CONCLUSIONS: A genomic region on BTA23 that is strongly associated with host susceptibility to bTB infection was identified. This region contained FKBP5, a gene involved in the TNFα/NFκ-B signalling pathway, which is a major biological pathway associated with immune response. Although there is no study that validates this region in the literature, our approach represents one of the most powerful studies for the analysis of bTB susceptibility to date.

Homozygosity mapping in an Irish ALS case-control cohort describes local demographic phenomena and points towards potential recessive risk loci.

Runs of homozygosity are common in European populations and are indicative of consanguinity, restricted population size and recessively inherited traits. Here, we map runs of homozygosity (ROHs) in an Irish case-control cohort for amyotrophic lateral sclerosis (ALS), a devastating neurological condition with high heritability yet only partially established genetic cause. We compare the extent of homozygosity in the Irish cohort with a large British cohort and observe that ROHs are longer and more frequent in the Irish population than in the British, and that extent of ROHs is correlated with demographic factors within the island of Ireland. ROHs are also longer and more frequent in ALS cases compared to population-matched controls, supporting the hypothesis that recessively inherited loci play a pathogenic role in ALS. Comparing homozygous haplotypes between cases and controls reveals several potential recessive risk loci for ALS, including a genomic interval spanning ARHGEF1, a compelling ALS candidate gene.

Comprehensive analysis of the gene encoding filaggrin uncovers prevalent and rare mutations in ichthyosis vulgaris and atopic eczema.

We recently reported two common filaggrin (FLG) null mutations that cause ichthyosis vulgaris and predispose to eczema and secondary allergic diseases. We show here that these common European mutations are ancestral variants carried on conserved haplotypes. To facilitate comprehensive analysis of other populations, we report a strategy for full sequencing of this large, highly repetitive gene, and we describe 15 variants, including seven that are prevalent. All the variants are either nonsense or frameshift mutations that, in representative cases, resulted in loss of filaggrin production in the epidermis. In an Irish case-control study, the five most common European mutations showed a strong association with moderate-to-severe childhood eczema (chi2 test: P = 2.12 x 10(-51); Fisher's exact test: heterozygote odds ratio (OR) = 7.44 (95% confidence interval (c.i.) = 4.9-11.3), and homozygote OR = 151 (95% c.i. = 20-1,136)). We found three additional rare null mutations in this case series, suggesting that the genetic architecture of filaggrin-related atopic dermatitis consists of both prevalent and rare risk alleles.