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 IUCN. 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 to 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 (Ovis gmelini). Studying diversity indices, we detect a considerable load of short runs of homozygosity (ROH) blocks (<2 Mb) in both Anatolian and Cyprian mouflons, reflecting small effective population size (Ne). Meanwhile, Ne as well as mutation load estimates are lower in Cyprian compared to 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 IUCN conservation status. Factors such as recent inbreeding, introgression, and unique population dynamics may have contributed to the observed disparities.

Benchmarking kinship estimation tools for ancient genomes using pedigree simulations.

There is growing interest in uncovering genetic kinship patterns in past societies using low-coverage palaeogenomes. Here, we benchmark four tools for kinship estimation with such data: lcMLkin, NgsRelate, KIN, and READ, which differ in their input, IBD estimation methods, and statistical approaches. We used pedigree and ancient genome sequence simulations to evaluate these tools when only a limited number (1 to 50 K, with minor allele frequency ≥0.01) of shared SNPs are available. The performance of all four tools was comparable using ≥20 K SNPs. We found that first-degree related pairs can be accurately classified even with 1 K SNPs, with 85% F1 scores using READ and 96% using NgsRelate or lcMLkin. Distinguishing third-degree relatives from unrelated pairs or second-degree relatives was also possible with high accuracy (F1 > 90%) with 5 K SNPs using NgsRelate and lcMLkin, while READ and KIN showed lower success (69 and 79% respectively). Meanwhile, noise in population allele frequencies and inbreeding (first-cousin mating) led to deviations in kinship coefficients, with different sensitivities across tools. We conclude that using multiple tools in parallel might be an effective approach to achieve robust estimates on ultra-low-coverage genomes.

aMeta: an accurate and memory-efficient ancient metagenomic profiling workflow.

Analysis of microbial data from archaeological samples is a growing field with great potential for understanding ancient environments, lifestyles, and diseases. However, high error rates have been a challenge in ancient metagenomics, and the availability of computational frameworks that meet the demands of the field is limited. Here, we propose aMeta, an accurate metagenomic profiling workflow for ancient DNA designed to minimize the amount of false discoveries and computer memory requirements. Using simulated data, we benchmark aMeta against a current state-of-the-art workflow and demonstrate its superiority in microbial detection and authentication, as well as substantially lower usage of computer memory.

Spatial and temporal heterogeneity in human mobility patterns in Holocene Southwest Asia and the East Mediterranean.

We present a spatiotemporal picture of human genetic diversity in Anatolia, Iran, Levant, South Caucasus, and the Aegean, a broad region that experienced the earliest Neolithic transition and the emergence of complex hierarchical societies. Combining 35 new ancient shotgun genomes with 382 ancient and 23 present-day published genomes, we found that genetic diversity within each region steadily increased through the Holocene. We further observed that the inferred sources of gene flow shifted in time. In the first half of the Holocene, Southwest Asian and the East Mediterranean populations homogenized among themselves. Starting with the Bronze Age, however, regional populations diverged from each other, most likely driven by gene flow from external sources, which we term "the expanding mobility model." Interestingly, this increase in inter-regional divergence can be captured by outgroup-f3-based genetic distances, but not by the commonly used FST statistic, due to the sensitivity of FST, but not outgroup-f3, to within-population diversity. Finally, we report a temporal trend of increasing male bias in admixture events through the Holocene.

A genomic snapshot of demographic and cultural dynamism in Upper Mesopotamia during the Neolithic Transition.

Upper Mesopotamia played a key role in the Neolithic Transition in Southwest Asia through marked innovations in symbolism, technology, and diet. We present 13 ancient genomes (c. 8500 to 7500 cal BCE) from Pre-Pottery Neolithic Çayönü in the Tigris basin together with bioarchaeological and material culture data. Our findings reveal that Çayönü was a genetically diverse population, carrying mixed ancestry from western and eastern Fertile Crescent, and that the community received immigrants. Our results further suggest that the community was organized along biological family lines. We document bodily interventions such as head shaping and cauterization among the individuals examined, reflecting Çayönü's cultural ingenuity. Last, we identify Upper Mesopotamia as the likely source of eastern gene flow into Neolithic Anatolia, in line with material culture evidence. We hypothesize that Upper Mesopotamia's cultural dynamism during the Neolithic Transition was the product not only of its fertile lands but also of its interregional demographic connections.

The fight for evolution.

As the Turkish government intensified its attacks on the theory of evolution, the academic community rallied to push back. A researcher recounts how she decided to join them.

Indian genetic heritage in Southeast Asian populations.

The great ethnolinguistic diversity found today in mainland Southeast Asia (MSEA) reflects multiple migration waves of people in the past. Maritime trading between MSEA and India was established at the latest 300 BCE, and the formation of early states in Southeast Asia during the first millennium CE was strongly influenced by Indian culture, a cultural influence that is still prominent today. Several ancient Indian-influenced states were located in present-day Thailand, and various populations in the country are likely to be descendants of people from those states. To systematically explore Indian genetic heritage in MSEA populations, we generated genome-wide SNP data (using the Affymetrix Human Origins array) for 119 present-day individuals belonging to 10 ethnic groups from Thailand and co-analyzed them with published data using PCA, ADMIXTURE, and methods relying on f-statistics and on autosomal haplotypes. We found low levels of South Asian admixture in various MSEA populations for whom there is evidence of historical connections with the ancient Indian-influenced states but failed to find this genetic component in present-day hunter-gatherer groups and relatively isolated groups from the highlands of Northern Thailand. The results suggest that migration of Indian populations to MSEA may have been responsible for the spread of Indian culture in the region. Our results also support close genetic affinity between Kra-Dai-speaking (also known as Tai-Kadai) and Austronesian-speaking populations, which fits a linguistic hypothesis suggesting cladality of the two language families.

Archaeogenetic analysis of Neolithic sheep from Anatolia suggests a complex demographic history since domestication.

Sheep were among the first domesticated animals, but their demographic history is little understood. Here we analyzed nuclear polymorphism and mitochondrial data (mtDNA) from ancient central and west Anatolian sheep dating from Epipaleolithic to late Neolithic, comparatively with modern-day breeds and central Asian Neolithic/Bronze Age sheep (OBI). Analyzing ancient nuclear data, we found that Anatolian Neolithic sheep (ANS) are genetically closest to present-day European breeds relative to Asian breeds, a conclusion supported by mtDNA haplogroup frequencies. In contrast, OBI showed higher genetic affinity to present-day Asian breeds. These results suggest that the east-west genetic structure observed in present-day breeds had already emerged by 6000 BCE, hinting at multiple sheep domestication episodes or early wild introgression in southwest Asia. Furthermore, we found that ANS are genetically distinct from all modern breeds. Our results suggest that European and Anatolian domestic sheep gene pools have been strongly remolded since the Neolithic.

Human inbreeding has decreased in time through the Holocene.

The history of human inbreeding is controversial.1 In particular, how the development of sedentary and/or agricultural societies may have influenced overall inbreeding levels, relative to those of hunter-gatherer communities, is unclear.2-5 Here, we present an approach for reliable estimation of runs of homozygosity (ROHs) in genomes with ≥3× mean sequence coverage across >1 million SNPs and apply this to 411 ancient Eurasian genomes from the last 15,000 years.5-34 We show that the frequency of inbreeding, as measured by ROHs, has decreased over time. The strongest effect is associated with the Neolithic transition, but the trend has since continued, indicating a population size effect on inbreeding prevalence. We further show that most inbreeding in our historical sample can be attributed to small population size instead of consanguinity. Cases of high consanguinity were rare and only observed among members of farming societies in our sample. Despite the lack of evidence for common consanguinity in our ancient sample, consanguineous traditions are today prevalent in various modern-day Eurasian societies,1,35-37 suggesting that such practices may have become widespread within the last few millennia.

Variable kinship patterns in Neolithic Anatolia revealed by ancient genomes.

The social organization of the first fully sedentary societies that emerged during the Neolithic period in Southwest Asia remains enigmatic,1 mainly because material culture studies provide limited insight into this issue. However, because Neolithic Anatolian communities often buried their dead beneath domestic buildings,2 household composition and social structure can be studied through these human remains. Here, we describe genetic relatedness among co-burials associated with domestic buildings in Neolithic Anatolia using 59 ancient genomes, including 22 new genomes from Asikli Höyük and Çatalhöyük. We infer pedigree relationships by simultaneously analyzing multiple types of information, including autosomal and X chromosome kinship coefficients, maternal markers, and radiocarbon dating. In two early Neolithic villages dating to the 9th and 8th millennia BCE, Asikli Höyük and Boncuklu, we discover that siblings and parent-offspring pairings were frequent within domestic structures, which provides the first direct indication of close genetic relationships among co-burials. In contrast, in the 7th millennium BCE sites of Çatalhöyük and Barcin, where we study subadults interred within and around houses, we find close genetic relatives to be rare. Hence, genetic relatedness may not have played a major role in the choice of burial location at these latter two sites, at least for subadults. This supports the hypothesis that in Çatalhöyük,3-5 and possibly in some other Neolithic communities, domestic structures may have served as burial location for social units incorporating biologically unrelated individuals. Our results underscore the diversity of kin structures in Neolithic communities during this important phase of sociocultural development.

Palaeo-Eskimo genetic ancestry and the peopling of Chukotka and North America.

Much of the American Arctic was first settled 5,000 years ago, by groups of people known as Palaeo-Eskimos. They were subsequently joined and largely displaced around 1,000 years ago by ancestors of the present-day Inuit and Yup'ik1-3. The genetic relationship between Palaeo-Eskimos and Native American, Inuit, Yup'ik and Aleut populations remains uncertain4-6. Here we present genomic data for 48 ancient individuals from Chukotka, East Siberia, the Aleutian Islands, Alaska, and the Canadian Arctic. We co-analyse these data with data from present-day Alaskan Iñupiat and West Siberian populations and published genomes. Using methods based on rare-allele and haplotype sharing, as well as established techniques4,7-9, we show that Palaeo-Eskimo-related ancestry is ubiquitous among people who speak Na-Dene and Eskimo-Aleut languages. We develop a comprehensive model for the Holocene peopling events of Chukotka and North America, and show that Na-Dene-speaking peoples, people of the Aleutian Islands, and Yup'ik and Inuit across the Arctic region all share ancestry from a single Palaeo-Eskimo-related Siberian source.

Genomic study of the Ket: a Paleo-Eskimo-related ethnic group with significant ancient North Eurasian ancestry.

The Kets, an ethnic group in the Yenisei River basin, Russia, are considered the last nomadic hunter-gatherers of Siberia, and Ket language has no transparent affiliation with any language family. We investigated connections between the Kets and Siberian and North American populations, with emphasis on the Mal'ta and Paleo-Eskimo ancient genomes, using original data from 46 unrelated samples of Kets and 42 samples of their neighboring ethnic groups (Uralic-speaking Nganasans, Enets, and Selkups). We genotyped over 130,000 autosomal SNPs, identified mitochondrial and Y-chromosomal haplogroups, and performed high-coverage genome sequencing of two Ket individuals. We established that Nganasans, Kets, Selkups, and Yukaghirs form a cluster of populations most closely related to Paleo-Eskimos in Siberia (not considering indigenous populations of Chukotka and Kamchatka). Kets are closely related to modern Selkups and to some Bronze and Iron Age populations of the Altai region, with all these groups sharing a high degree of Mal'ta ancestry. Implications of these findings for the linguistic hypothesis uniting Ket and Na-Dene languages into a language macrofamily are discussed.