Ancient mitogenomes from Pre-Pottery Neolithic Central Anatolia and the effects of a Late Neolithic bottleneck in sheep (Ovis aries).

Occupied between ~10,300 and 9300 years ago, the Pre-Pottery Neolithic site of Asikli Höyük in Central Anatolia went through early phases of sheep domestication. Analysis of 629 mitochondrial genomes from this and numerous sites in Anatolia, southwest Asia, Europe, and Africa produced a phylogenetic tree with excessive coalescences (nodes) around the Neolithic, a potential signature of a domestication bottleneck. This is consistent with archeological evidence of sheep management at Asikli Höyük which transitioned from residential stabling to open pasturing over a millennium of site occupation. However, unexpectedly, we detected high genetic diversity throughout Asikli Höyük's occupation rather than a bottleneck. Instead, we detected a tenfold demographic bottleneck later in the Neolithic, which caused the fixation of mitochondrial haplogroup B in southwestern Anatolia. The mitochondrial genetic makeup that emerged was carried from the core region of early Neolithic sheep management into Europe and dominates the matrilineal diversity of both its ancient and the billion-strong modern sheep populations.

Biological and substitute parents in Beaker period adult-child graves.

Joint inhumations of adults and children are an intriguing aspect of the shift from collective to single burial rites in third millennium BC Western Eurasia. Here, we revisit two exceptional Beaker period adult-child graves using ancient DNA: Altwies in Luxembourg and Dunstable Downs in Britain. Ancestry modelling and patterns of shared IBD segments between the individuals examined, and contemporary genomes from Central and Northwest Europe, highlight the continental connections of British Beakers. Although simultaneous burials may involve individuals with no social or biological ties, we present evidence that close blood relations played a role in shaping third millennium BC social systems and burial practices, for example a biological mother and her son buried together at Altwies. Extended family, such as a paternal aunt at Dunstable Downs, could also act as 'substitute parents' in the grave. Hypotheses are explored to explain such simultaneous inhumations. Whilst intercommunity violence, infectious disease and epidemics may be considered as explanations, they fail to account for both the specific, codified nature of this particular form of inhumation, and its pervasiveness, as evidenced by a representative sample of 131 adult-child graves from 88 sites across Eurasia, all dating to the third and second millennia BC.

Descent, marriage, and residence practices of a 3,800-year-old pastoral community in Central Eurasia.

Our understanding of prehistoric societal organization at the family level is still limited. Here, we generated genome data from 32 individuals from an approximately 3,800-y-old burial mound attributed to the Bronze Age Srubnaya-Alakul cultural tradition at the site of Nepluyevsky, located in the Southern Ural region of Central Eurasia. We found that life expectancy was generally very low, with adult males living on average 8 y longer than females. A total of 35 first-degree, 40 second-degree, and 48 third-degree biological relationships connected 23 of the studied individuals, allowing us to propose a family tree spanning three generations with six brothers at its center. The oldest of these brothers had eight children with two women and the most children overall, whereas the other relationships were monogamous. Notably, related female children above the age of five were completely absent from the site, and adult females were more genetically diverse than males. These results suggest that biological relationships between male siblings played a structural role in society and that descent group membership was based on patrilineality. Women originated from a larger mating network and moved to join the men, with whom they were buried. Finally, the oldest brother likely held a higher social position, which was expressed in terms of fertility.

Between Raetia Secunda and the dutchy of Bavaria: Exploring patterns of human movement and diet.

During the transition from Late Antiquity to the Middle Ages, the Roman Empire dissolved in the West and medieval empires were founded. There has been much discussion about the role that migration played in this transition. This is especially true for the formation of the Baiuvariian tribe and the founding of this tribal dukedom, which took place from the 5th to the 6th century in what is now Southern Bavaria (Germany). In this study, we aimed to determine the extent of immigration during the beginning of this transformation and to shed further light on its character. To achieve this goal, we analyzed stable isotope values of strontium, carbon, and nitrogen from the teeth and bones of over 150 human remains from Southern Germany, dating from around 500 AD. This group of individuals included women with cranial modifications (ACD) which can be found sporadically in the burial grounds of this period. Our results showed an above-average migration rate for both men and women in the second half of the 5th century. They also indicate that a foreign background may also be assumed for the women with ACD. The demonstrably different origins of the immigrants from isotopically diverse regions, and the identification of local differences in detectable migration rate, as well as indication for different timing of residential changes, highlight the complexity of immigration processes and the need for more studies at the regional level.

HaploBlocks: Efficient Detection of Positive Selection in Large Population Genomic Datasets.

Genomic regions under positive selection harbor variation linked for example to adaptation. Most tools for detecting positively selected variants have computational resource requirements rendering them impractical on population genomic datasets with hundreds of thousands of individuals or more. We have developed and implemented an efficient haplotype-based approach able to scan large datasets and accurately detect positive selection. We achieve this by combining a pattern matching approach based on the positional Burrows-Wheeler transform with model-based inference which only requires the evaluation of closed-form expressions. We evaluate our approach with simulations, and find it to be both sensitive and specific. The computational resource requirements quantified using UK Biobank data indicate that our implementation is scalable to population genomic datasets with millions of individuals. Our approach may serve as an algorithmic blueprint for the era of "big data" genomics: a combinatorial core coupled with statistical inference in closed form.

Thresholds for the presence of glacial megafauna in central Europe during the last 60,000 years.

Lake sediment records from Holzmaar and the infilled maar of Auel (Eifel, Germany) are used to reconstruct landscape changes and megafauna abundances. Our data document a forested landscape from 60,000 to 48,000 yr b2k and a stepwise vegetation change towards a glacial desert after 26,000 yr b2k. The Eifel landscape was continuously inhabited from 48,000 to 9000 yr b2k by large mammals, documented by the presence of spores of coprophilous fungi from Sordaria and Sporormiella fungi that grow on fecal remains of the megafauna. Megafauna reached higher numbers during cold stadial climates but was present also during the warmer interstadials. Highest abundance was at 56,500/48,500/38,500/33,000/27,000/21,000/16,200/14,000 yr b2k, i.e. under different climate regimes. Some of these dates were associated with clear human presence, which indicates that megafauna was not overkilled by humans. In contrast, human presence could quite likely have been stimulated by the abundant food supply. Megafauna presence decreased significantly when tree abundance increased during interstadials. The Megafauna disappeared finally at 11,400 yr b2k with the development of the early Holocene forest cover, which appears to be the most important threshold for megafauna presence.

Ancient mitochondrial diversity reveals population homogeneity in Neolithic Greece and identifies population dynamics along the Danubian expansion axis.

The aim of the study is to investigate mitochondrial diversity in Neolithic Greece and its relation to hunter-gatherers and farmers who populated the Danubian Neolithic expansion axis. We sequenced 42 mitochondrial palaeogenomes from Greece and analysed them together with European set of 328 mtDNA sequences dating from the Early to the Final Neolithic and 319 modern sequences. To test for population continuity through time in Greece, we use an original structured population continuity test that simulates DNA from different periods by explicitly considering the spatial and temporal dynamics of populations. We explore specific scenarios of the mode and tempo of the European Neolithic expansion along the Danubian axis applying spatially explicit simulations coupled with Approximate Bayesian Computation. We observe a striking genetic homogeneity for the maternal line throughout the Neolithic in Greece whereas population continuity is rejected between the Neolithic and present-day Greeks. Along the Danubian expansion axis, our best-fitting scenario supports a substantial decrease in mobility and an increasing local hunter-gatherer contribution to the gene-pool of farmers following the initial rapid Neolithic expansion. Οur original simulation approach models key demographic parameters rather than inferring them from fragmentary data leading to a better understanding of this important process in European prehistory.

The genomic origins of the world's first farmers.

The precise genetic origins of the first Neolithic farming populations in Europe and Southwest Asia, as well as the processes and the timing of their differentiation, remain largely unknown. Demogenomic modeling of high-quality ancient genomes reveals that the early farmers of Anatolia and Europe emerged from a multiphase mixing of a Southwest Asian population with a strongly bottlenecked western hunter-gatherer population after the last glacial maximum. Moreover, the ancestors of the first farmers of Europe and Anatolia went through a period of extreme genetic drift during their westward range expansion, contributing highly to their genetic distinctiveness. This modeling elucidates the demographic processes at the root of the Neolithic transition and leads to a spatial interpretation of the population history of Southwest Asia and Europe during the late Pleistocene and early Holocene.

Trio exome sequencing is highly relevant in prenatal diagnostics.

OBJECTIVE: About 3% of newborns show malformations, with about 20% of the affected having genetic causes. Clarification of genetic diseases in postnatal diagnostics was significantly improved with high-throughput sequencing, in particular through whole exome sequencing covering all protein-coding regions. Here, we aim to extend the use of this technology to prenatal diagnostics. METHOD: Between 07/2018 and 10/2020, 500 pregnancies with fetal ultrasound abnormalities were analyzed after genetic counseling as part of prenatal diagnostics using WES of the fetus and parents. RESULTS: Molecular genetic findings could explain ultrasound abnormalities in 38% of affected fetuses. In 47% of these, disease-causing de novo variants were found. Pathogenic variants in genes with autosomal recessive or X-linked inheritance were detected in more than one-third (70/189 = 37%). The latter are associated with increased probability of recurrence, making their detection important for further pregnancies. Average time from sample receipt to report was 12 days in the recent cases. CONCLUSION: Trio exome sequencing is a useful addition to prenatal diagnostics due to its high diagnostic yield and short processing time (comparable to chromosome analysis). It covers a wide spectrum of genetic changes. Comprehensive interdisciplinary counseling before and after diagnostics is indispensable.

The genomic history of the Aegean palatial civilizations.

The Cycladic, the Minoan, and the Helladic (Mycenaean) cultures define the Bronze Age (BA) of Greece. Urbanism, complex social structures, craft and agricultural specialization, and the earliest forms of writing characterize this iconic period. We sequenced six Early to Middle BA whole genomes, along with 11 mitochondrial genomes, sampled from the three BA cultures of the Aegean Sea. The Early BA (EBA) genomes are homogeneous and derive most of their ancestry from Neolithic Aegeans, contrary to earlier hypotheses that the Neolithic-EBA cultural transition was due to massive population turnover. EBA Aegeans were shaped by relatively small-scale migration from East of the Aegean, as evidenced by the Caucasus-related ancestry also detected in Anatolians. In contrast, Middle BA (MBA) individuals of northern Greece differ from EBA populations in showing ∼50% Pontic-Caspian Steppe-related ancestry, dated at ca. 2,600-2,000 BCE. Such gene flow events during the MBA contributed toward shaping present-day Greek genomes.

Reliability of genomic variants across different next-generation sequencing platforms and bioinformatic processing pipelines.

BACKGROUND: Next Generation Sequencing (NGS) is the fundament of various studies, providing insights into questions from biology and medicine. Nevertheless, integrating data from different experimental backgrounds can introduce strong biases. In order to methodically investigate the magnitude of systematic errors in single nucleotide variant calls, we performed a cross-sectional observational study on a genomic cohort of 99 subjects each sequenced via (i) Illumina HiSeq X, (ii) Illumina HiSeq, and (iii) Complete Genomics and processed with the respective bioinformatic pipeline. We also repeated variant calling for the Illumina cohorts with GATK, which allowed us to investigate the effect of the bioinformatics analysis strategy separately from the sequencing platform's impact. RESULTS: The number of detected variants/variant classes per individual was highly dependent on the experimental setup. We observed a statistically significant overrepresentation of variants uniquely called by a single setup, indicating potential systematic biases. Insertion/deletion polymorphisms (indels) were associated with decreased concordance compared to single nucleotide polymorphisms (SNPs). The discrepancies in indel absolute numbers were particularly prominent in introns, Alu elements, simple repeats, and regions with medium GC content. Notably, reprocessing sequencing data following the best practice recommendations of GATK considerably improved concordance between the respective setups. CONCLUSION: We provide empirical evidence of systematic heterogeneity in variant calls between alternative experimental and data analysis setups. Furthermore, our results demonstrate the benefit of reprocessing genomic data with harmonized pipelines when integrating data from different studies.

Low Prevalence of Lactase Persistence in Bronze Age Europe Indicates Ongoing Strong Selection over the Last 3,000 Years.

Lactase persistence (LP), the continued expression of lactase into adulthood, is the most strongly selected single gene trait over the last 10,000 years in multiple human populations. It has been posited that the primary allele causing LP among Eurasians, rs4988235-A [1], only rose to appreciable frequencies during the Bronze and Iron Ages [2, 3], long after humans started consuming milk from domesticated animals. This rapid rise has been attributed to an influx of people from the Pontic-Caspian steppe that began around 5,000 years ago [4, 5]. We investigate the spatiotemporal spread of LP through an analysis of 14 warriors from the Tollense Bronze Age battlefield in northern Germany (∼3,200 before present, BP), the oldest large-scale conflict site north of the Alps. Genetic data indicate that these individuals represent a single unstructured Central/Northern European population. We complemented these data with genotypes of 18 individuals from the Bronze Age site Mokrin in Serbia (∼4,100 to ∼3,700 BP) and 37 individuals from Eastern Europe and the Pontic-Caspian Steppe region, predating both Bronze Age sites (∼5,980 to ∼3,980 BP). We infer low LP in all three regions, i.e., in northern Germany and South-eastern and Eastern Europe, suggesting that the surge of rs4988235 in Central and Northern Europe was unlikely caused by Steppe expansions. We estimate a selection coefficient of 0.06 and conclude that the selection was ongoing in various parts of Europe over the last 3,000 years.

Screening archaeological bone for palaeogenetic and palaeoproteomic studies.

The recovery and analysis of ancient DNA and protein from archaeological bone is time-consuming and expensive to carry out, while it involves the partial or complete destruction of valuable or rare specimens. The fields of palaeogenetic and palaeoproteomic research would benefit greatly from techniques that can assess the molecular quality prior to sampling. To be relevant, such screening methods should be effective, minimally-destructive, and rapid. This study reports results based on spectroscopic (Fourier-transform infrared spectroscopy in attenuated total reflectance [FTIR-ATR]; n = 266), palaeoproteomic (collagen content; n = 226), and palaeogenetic (endogenous DNA content; n = 88) techniques. We establish thresholds for three different FTIR indices, a) the infrared splitting factor [IRSF] that assesses relative changes in bioapatite crystals' size and homogeneity; b) the carbonate-to-phosphate [C/P] ratio as a relative measure of carbonate content in bioapatite crystals; and c) the amide-to-phosphate ratio [Am/P] for assessing the relative organic content preserved in bone. These thresholds are both extremely reliable and easy to apply for the successful and rapid distinction between well- and poorly-preserved specimens. This is a milestone for choosing appropriate samples prior to genomic and collagen analyses, with important implications for biomolecular archaeology and palaeontology.

Ancient cattle genomics, origins, and rapid turnover in the Fertile Crescent.

Genome-wide analysis of 67 ancient Near Eastern cattle, Bos taurus, remains reveals regional variation that has since been obscured by admixture in modern populations. Comparisons of genomes of early domestic cattle to their aurochs progenitors identify diverse origins with separate introgressions of wild stock. A later region-wide Bronze Age shift indicates rapid and widespread introgression of zebu, Bos indicus, from the Indus Valley. This process was likely stimulated at the onset of the current geological age, ~4.2 thousand years ago, by a widespread multicentury drought. In contrast to genome-wide admixture, mitochondrial DNA stasis supports that this introgression was male-driven, suggesting that selection of arid-adapted zebu bulls enhanced herd survival. This human-mediated migration of zebu-derived genetics has continued through millennia, altering tropical herding on each continent.

Tracking Five Millennia of Horse Management with Extensive Ancient Genome Time Series.

Horse domestication revolutionized warfare and accelerated travel, trade, and the geographic expansion of languages. Here, we present the largest DNA time series for a non-human organism to date, including genome-scale data from 149 ancient animals and 129 ancient genomes (≥1-fold coverage), 87 of which are new. This extensive dataset allows us to assess the modern legacy of past equestrian civilizations. We find that two extinct horse lineages existed during early domestication, one at the far western (Iberia) and the other at the far eastern range (Siberia) of Eurasia. None of these contributed significantly to modern diversity. We show that the influence of Persian-related horse lineages increased following the Islamic conquests in Europe and Asia. Multiple alleles associated with elite-racing, including at the MSTN "speed gene," only rose in popularity within the last millennium. Finally, the development of modern breeding impacted genetic diversity more dramatically than the previous millennia of human management.

Disruption of KCNQ1 prevents methylation of the ICR2 and supports the hypothesis that its transcription is necessary for imprint establishment.

Beckwith-Wiedemann syndrome (BWS; OMIM #130650) is an imprinting disorder caused by genetic or epigenetic alterations of one or both imprinting control regions on chromosome 11p15.5. Hypomethylation of the centromeric imprinting control region (KCNQ1OT1:TSS-DMR, ICR2) is the most common molecular cause of BWS and is present in about half of the cases. Based on a BWS family with a maternal deletion of the 5' part of KCNQ1 we have recently hypothesised that transcription of KCNQ1 is a prerequisite for the establishment of methylation at the KCNQ1OT1:TSS-DMR in the oocyte. Further evidence for this hypothesis came from a mouse model where methylation failed to be established when a poly(A) truncation cassette was inserted into this locus to prevent transcription through the DMR. Here we report on a family where a balanced translocation disrupts the KCNQ1 gene in intron 9. Maternal inheritance of this translocation is associated with hypomethylation of the KCNQ1OT1:TSS-DMR and BWS. This finding strongly supports our previous hypothesis that transcription of KCNQ1 is required for establishing the maternal methylation imprint at the KCNQ1OT1:TSS-DMR.

Ancient goat genomes reveal mosaic domestication in the Fertile Crescent.

Current genetic data are equivocal as to whether goat domestication occurred multiple times or was a singular process. We generated genomic data from 83 ancient goats (51 with genome-wide coverage) from Paleolithic to Medieval contexts throughout the Near East. Our findings demonstrate that multiple divergent ancient wild goat sources were domesticated in a dispersed process that resulted in genetically and geographically distinct Neolithic goat populations, echoing contemporaneous human divergence across the region. These early goat populations contributed differently to modern goats in Asia, Africa, and Europe. We also detect early selection for pigmentation, stature, reproduction, milking, and response to dietary change, providing 8000-year-old evidence for human agency in molding genome variation within a partner species.

Population genomic analysis of elongated skulls reveals extensive female-biased immigration in Early Medieval Bavaria.

Modern European genetic structure demonstrates strong correlations with geography, while genetic analysis of prehistoric humans has indicated at least two major waves of immigration from outside the continent during periods of cultural change. However, population-level genome data that could shed light on the demographic processes occurring during the intervening periods have been absent. Therefore, we generated genomic data from 41 individuals dating mostly to the late 5th/early 6th century AD from present-day Bavaria in southern Germany, including 11 whole genomes (mean depth 5.56×). In addition we developed a capture array to sequence neutral regions spanning a total of 5 Mb and 486 functional polymorphic sites to high depth (mean 72×) in all individuals. Our data indicate that while men generally had ancestry that closely resembles modern northern and central Europeans, women exhibit a very high genetic heterogeneity; this includes signals of genetic ancestry ranging from western Europe to East Asia. Particularly striking are women with artificial skull deformations; the analysis of their collective genetic ancestry suggests an origin in southeastern Europe. In addition, functional variants indicate that they also differed in visible characteristics. This example of female-biased migration indicates that complex demographic processes during the Early Medieval period may have contributed in an unexpected way to shape the modern European genetic landscape. Examination of the panel of functional loci also revealed that many alleles associated with recent positive selection were already at modern-like frequencies in European populations ∼1,500 years ago.

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.

Ancient European dog genomes reveal continuity since the Early Neolithic.

Europe has played a major role in dog evolution, harbouring the oldest uncontested Palaeolithic remains and having been the centre of modern dog breed creation. Here we sequence the genomes of an Early and End Neolithic dog from Germany, including a sample associated with an early European farming community. Both dogs demonstrate continuity with each other and predominantly share ancestry with modern European dogs, contradicting a previously suggested Late Neolithic population replacement. We find no genetic evidence to support the recent hypothesis proposing dual origins of dog domestication. By calibrating the mutation rate using our oldest dog, we narrow the timing of dog domestication to 20,000-40,000 years ago. Interestingly, we do not observe the extreme copy number expansion of the AMY2B gene characteristic of modern dogs that has previously been proposed as an adaptation to a starch-rich diet driven by the widespread adoption of agriculture in the Neolithic.