The Beaker phenomenon and the genomic transformation of northwest Europe

Bioarchaeology aids the cultural understanding of six characters in search of their agency (Tarquinia, ninth-seventh century BC, central Italy)

Etruria contained one of the great early urban civilisations in the Italian peninsula during the first millennium BC, much studied from a cultural, humanities-based, perspective, but relatively little with scientific data, and rarely in combination. We have addressed the unusual location of twenty inhumations found in the sacred heart of the Etruscan city of Tarquinia, focusing on six of these as illustrative, contrasting with the typical contemporary cremations found in cemeteries on the edge of the city. The cultural evidence suggests that the six skeletons were also distinctive in their ritualization and memorialisation. Focusing on the six, as a representative sample, the scientific evidence of osteoarchaeology, isotopic compositions, and ancient DNA has established that these appear to show mobility, diversity and violence through an integrated bioarchaeological approach. The combination of multiple lines of evidence makes major strides towards a deeper understanding of the role of these extraordinary individuals in the life of the early city of Etruria.

Phylogenetic insights into the genetic legacies of Hungarian-speaking communities in the Carpathian Basin

This study focuses on exploring the uniparental genetic lineages of Hungarian-speaking minorities residing in rural villages of Baranja (Croatia) and the Zobor region (Slovakia). We aimed to identify ancestral lineages by examining genetic markers distributed across the entire mitogenome and on the Y-chromosome. This allowed us to discern disparities in regional genetic structures within these communities. By integrating our newly acquired genetic data from a total of 168 participants with pre-existing Eurasian and ancient DNA datasets, our goal was to enrich the understanding of the genetic history trajectories of Carpathian Basin populations. Our findings suggest that while population-based analyses may not be sufficiently robust to detect fine-scale uniparental genetic patterns with the sample sizes at hand, phylogenetic analysis of well-characterized Y-chromosomal Short Tandem Repeat (STR) data and entire mitogenome sequences did uncover multiple lineage ties to far-flung regions and eras. While the predominant portions of both paternal and maternal DNA align with the East-Central European spectrum, rarer subhaplogroups and lineages have unveiled ancient ties to both prehistoric and historic populations spanning Europe and Eastern Eurasia. This research augments the expansive field of phylogenetics, offering critical perspectives on the genetic constitution and heritage of the communities in East-Central Europe.

A novel expectation-maximization approach to infer general diploid selection from time-series genetic data

Detecting and quantifying the strength of selection is a main objective in population genetics. Since selection acts over multiple generations, many approaches have been developed to detect and quantify selection using genetic data sampled at multiple points in time. Such time series genetic data is commonly analyzed using Hidden Markov Models, but in most cases, under the assumption of additive selection. However, many examples of genetic variation exhibiting non-additive mechanisms exist, making it critical to develop methods that can characterize selection in more general scenarios. Thus, we extend a previously introduced expectation-maximization algorithm for the inference of additive selection coefficients to the case of general diploid selection, in which heterozygote and homozygote fitnesses are parameterized independently. We furthermore introduce a framework to identify bespoke modes of diploid selection from given data, as well as a procedure for aggregating data across linked loci to increase power and robustness. Using extensive simulation studies, we find that our method accurately and efficiently estimates selection coefficients for different modes of diploid selection across a wide range of scenarios; however, power to classify the mode of selection is low unless selection is very strong. We apply our method to ancient DNA samples from Great Britain in the last 4,450 years, and detect evidence for selection in six genomic regions, including the well-characterized LCT locus. Our work is the first genome-wide scan characterizing signals of general diploid selection.

The Genetic Origin of the Indo-Europeans

The Yamnaya archaeological complex appeared around 3300BCE across the steppes north of the Black and Caspian Seas, and by 3000BCE reached its maximal extent from Hungary in the west to Kazakhstan in the east. To localize the ancestral and geographical origins of the Yamnaya among the diverse Eneolithic people that preceded them, we studied ancient DNA data from 428 individuals of which 299 are reported for the first time, demonstrating three previously unknown Eneolithic genetic clines. First, a "Caucasus-Lower Volga" (CLV) Cline suffused with Caucasus hunter-gatherer (CHG) ancestry extended between a Caucasus Neolithic southern end in Neolithic Armenia, and a steppe northern end in Berezhnovka in the Lower Volga. Bidirectional gene flow across the CLV cline created admixed intermediate populations in both the north Caucasus, such as the Maikop people, and on the steppe, such as those at the site of Remontnoye north of the Manych depression. CLV people also helped form two major riverine clines by admixing with distinct groups of European hunter-gatherers. A "Volga Cline" was formed as Lower Volga people mixed with upriver populations that had more Eastern hunter-gatherer (EHG) ancestry, creating genetically hyper-variable populations as at Khvalynsk in the Middle Volga. A "Dnipro Cline" was formed as CLV people bearing both Caucasus Neolithic and Lower Volga ancestry moved west and acquired Ukraine Neolithic hunter-gatherer (UNHG) ancestry to establish the population of the Serednii Stih culture from which the direct ancestors of the Yamnaya themselves were formed around 4000BCE. This population grew rapidly after 3750-3350BCE, precipitating the expansion of people of the Yamnaya culture who totally displaced previous groups on the Volga and further east, while admixing with more sedentary groups in the west. CLV cline people with Lower Volga ancestry contributed four fifths of the ancestry of the Yamnaya, but also, entering Anatolia from the east, contributed at least a tenth of the ancestry of Bronze Age Central Anatolians, where the Hittite language, related to the Indo-European languages spread by the Yamnaya, was spoken. We thus propose that the final unity of the speakers of the "Proto-Indo-Anatolian" ancestral language of both Anatolian and Indo-European languages can be traced to CLV cline people sometime between 4400-4000 BCE.

Accurate Bayesian inference of sex chromosome karyotypes and sex-linked scaffolds from low-depth sequencing data

The identification of sex-linked scaffolds and the genetic sex of individuals, i.e. their sex karyotype, is a fundamental step in population genomic studies. If sex-linked scaffolds are known, single individuals may be sexed based on read counts of next-generation sequencing data. If both sex-linked scaffolds as well as sex karyotypes are unknown, as is often the case for non-model organisms, they have to be jointly inferred. For both cases, current methods rely on arbitrary thresholds, which limits their power for low-depth data. In addition, most current methods are limited to euploid sex karyotypes (XX and XY). Here we develop BeXY, a fully Bayesian method to jointly infer the posterior probabilities for each scaffold to be autosomal, X- or Y-linked and for each individual to be any of the sex karyotypes XX, XY, X0, XXX, XXY, XYY and XXYY. If the sex-linked scaffolds are known, it also identifies autosomal trisomies and estimates the sex karyotype posterior probabilities for single individuals. As we show with downsampling experiments, BeXY has higher power than all existing methods. It accurately infers the sex karyotype of ancient human samples with as few as 20,000 reads and accurately infers sex-linked scaffolds from data sets of just a handful of samples or with highly imbalanced sex ratios, also in the case of low-quality reference assemblies. We illustrate the power of BeXY by applying it to both whole-genome shotgun and target enrichment sequencing data of ancient and modern humans, as well as several non-model organisms.

Predicting Functional Consequences of Recent Natural Selection in Britain

Ancient DNA can directly reveal the contribution of natural selection to human genomic variation. However, while the analysis of ancient DNA has been successful at identifying genomic signals of selection, inferring the phenotypic consequences of that selection has been more difficult. Most trait-associated variants are noncoding, so we expect that a large proportion of the phenotypic effects of selection will also act through noncoding variation. Since we cannot measure gene expression directly in ancient individuals, we used an approach (Joint-Tissue Imputation [JTI]) developed to predict gene expression from genotype data. We tested for changes in the predicted expression of 17,384 protein coding genes over a time transect of 4,500 years using 91 present-day and 616 ancient individuals from Britain. We identified 28 genes at seven genomic loci with significant (false discovery rate [FDR] < 0.05) changes in predicted expression levels in this time period. We compared the results from our transcriptome-wide scan to a genome-wide scan based on estimating per-single nucleotide polymorphism (SNP) selection coefficients from time series data. At five previously identified loci, our approach allowed us to highlight small numbers of genes with evidence for significant shifts in expression from peaks that in some cases span tens of genes. At two novel loci (SLC44A5 and NUP85), we identify selection on gene expression not captured by scans based on genomic signatures of selection. Finally, we show how classical selection statistics (iHS and SDS) can be combined with JTI models to incorporate functional information into scans that use present-day data alone. These results demonstrate the potential of this type of information to explore both the causes and consequences of natural selection.

Reconstructing DNA methylation maps of ancient populations

Studying premortem DNA methylation from ancient DNA (aDNA) provides a proxy for ancient gene activity patterns, and hence valuable information on evolutionary changes in gene regulation. Due to statistical limitations, current methods to reconstruct aDNA methylation maps are constrained to high-coverage shotgun samples, which comprise a small minority of available ancient samples. Most samples are sequenced using in-situ hybridization capture sequencing which targets a predefined set of genomic positions. Here, we develop methods to reconstruct aDNA methylation maps of samples that were not sequenced using high-coverage shotgun sequencing, by way of pooling together individuals to obtain a DNA methylation map that is characteristic of a population. We show that the resulting DNA methylation maps capture meaningful biological information and allow for the detection of differential methylation across populations. We offer guidelines on how to carry out comparative studies involving ancient populations, and how to control the rate of falsely discovered differentially methylated regions. The ability to reconstruct DNA methylation maps of past populations allows for the development of a whole new frontier in paleoepigenetic research, tracing DNA methylation changes throughout human history, using data from thousands of ancient samples.

The contribution of gene flow, selection, and genetic drift to five thousand years of human allele frequency change

Genomic time series from experimental evolution studies and ancient DNA datasets offer us a chance to directly observe the interplay of various evolutionary forces. We show how the genome-wide variance in allele frequency change between two time points can be decomposed into the contributions of gene flow, genetic drift, and linked selection. In closed populations, the contribution of linked selection is identifiable because it creates covariances between time intervals, and genetic drift does not. However, repeated gene flow between populations can also produce directionality in allele frequency change, creating covariances. We show how to accurately separate the fraction of variance in allele frequency change due to admixture and linked selection in a population receiving gene flow. We use two human ancient DNA datasets, spanning around 5,000 y, as time transects to quantify the contributions to the genome-wide variance in allele frequency change. We find that a large fraction of genome-wide change is due to gene flow. In both cases, after correcting for known major gene flow events, we do not observe a signal of genome-wide linked selection. Thus despite the known role of selection in shaping long-term polymorphism levels, and an increasing number of examples of strong selection on single loci and polygenic scores from ancient DNA, it appears to be gene flow and drift, and not selection, that are the main determinants of recent genome-wide allele frequency change. Our approach should be applicable to the growing number of contemporary and ancient temporal population genomics datasets.

Lead isotopes of prehistoric copper tools define metallurgical phases in Late Neolithic and Eneolithic Italy

The diffuse presence of small copper ore deposits in the Alpine area, mostly exploited since Late Medieval times, led most scholars to assume that these deposits may actually be active much earlier and that many of the circulating prehistoric metal objects found in the area were produced with local copper sources. This assumption was recently validated for the Recent Bronze Age through the use of lead isotope tracers, and well supported by the archaeometallurgical evidences found in the South-Eastern Alps. However, the scarcity of available lead isotope data for pre-Bronze Age metals precluded to date the reconstruction of the metal flow through the Late Neolithic and Eneolithic (or Copper Age). Based on 49 new analyses of important archaeological objects from the Alpine region, the Po River Valley and Central Italy, mostly axes dated from the Late Neolithic to the Late Eneolithic, here we show that the diffusion of copper in Northern Italy (approximately 4500-2200 BC) includes three major periods of metal use and/or production, each related to specific ore sources. The South Alpine copper was massively used only starting from the middle of the 3rd millennium BC, in connection or slightly earlier than the Beaker event.

The Allen Ancient DNA Resource (AADR) a curated compendium of ancient human genomes

More than two hundred papers have reported genome-wide data from ancient humans. While the raw data for the vast majority are fully publicly available testifying to the commitment of the paleogenomics community to open data, formats for both raw data and meta-data differ. There is thus a need for uniform curation and a centralized, version-controlled compendium that researchers can download, analyze, and reference. Since 2019, we have been maintaining the Allen Ancient DNA Resource (AADR), which aims to provide an up-to-date, curated version of the world's published ancient human DNA data, represented at more than a million single nucleotide polymorphisms (SNPs) at which almost all ancient individuals have been assayed. The AADR has gone through six public releases at the time of writing and review of this manuscript, and crossed the threshold of >10,000 individuals with published genome-wide ancient DNA data at the end of 2022. This note is intended as a citable descriptor of the AADR.

Evolutionary Trends of Polygenic Scores in European Populations From the Paleolithic to Modern Times

This study examines the temporal and geographical evolution of polygenic scores (PGSs) across cognitive measures (Educational Attainment [EA], Intelligence Quotient [IQ]), Socioeconomic Status (SES), and psychiatric conditions (Autism Spectrum Disorder [ASD], schizophrenia [SCZ]) in various populations. Our findings indicate positive directional selection for EA, IQ, and SES traits over the past 12,000 years. Schizophrenia and autism, while similar, showed different temporal patterns, aligning with theories suggesting they are psychological opposites. We observed a decline in PGS for neuroticism and depression, likely due to their genetic correlations and pleiotropic effects on intelligence. Significant PGS shifts from the Upper Paleolithic to the Neolithic periods suggest lifestyle and cognitive demand changes, particularly during the Neolithic Revolution. The study supports a mild hypothesis of Gregory Clark's model, showing a noticeable rise in genetic propensities for intelligence, academic achievement and professional status across Europe from the Middle Ages to the present. While latitude strongly influenced height, its impact on schizophrenia and autism was smaller and varied. Contrary to the cold winters theory, the study found no significant correlation between latitude and intelligence.

Stable population structure in Europe since the Iron Age, despite high mobility

Ancient DNA research in the past decade has revealed that European population structure changed dramatically in the prehistoric period (14,000-3000 years before present, YBP), reflecting the widespread introduction of Neolithic farmer and Bronze Age Steppe ancestries. However, little is known about how population structure changed from the historical period onward (3000 YBP - present). To address this, we collected whole genomes from 204 individuals from Europe and the Mediterranean, many of which are the first historical period genomes from their region (e.g. Armenia and France). We found that most regions show remarkable inter-individual heterogeneity. At least 7% of historical individuals carry ancestry uncommon in the region where they were sampled, some indicating cross-Mediterranean contacts. Despite this high level of mobility, overall population structure across western Eurasia is relatively stable through the historical period up to the present, mirroring geography. We show that, under standard population genetics models with local panmixia, the observed level of dispersal would lead to a collapse of population structure. Persistent population structure thus suggests a lower effective migration rate than indicated by the observed dispersal. We hypothesize that this phenomenon can be explained by extensive transient dispersal arising from drastically improved transportation networks and the Roman Empire's mobilization of people for trade, labor, and military. This work highlights the utility of ancient DNA in elucidating finer scale human population dynamics in recent history.

Genetic history of Cambridgeshire before and after the Black Death

The extent of the devastation of the Black Death pandemic (1346-1353) on European populations is known from documentary sources and its bacterial source illuminated by studies of ancient pathogen DNA. What has remained less understood is the effect of the pandemic on human mobility and genetic diversity at the local scale. Here, we report 275 ancient genomes, including 109 with coverage >0.1×, from later medieval and postmedieval Cambridgeshire of individuals buried before and after the Black Death. Consistent with the function of the institutions, we found a lack of close relatives among the friars and the inmates of the hospital in contrast to their abundance in general urban and rural parish communities. While we detect long-term shifts in local genetic ancestry in Cambridgeshire, we find no evidence of major changes in genetic ancestry nor higher differentiation of immune loci between cohorts living before and after the Black Death.

The contribution of gene flow, selection, and genetic drift to five thousand years of human allele frequency change

Genomic time series from experimental evolution studies and ancient DNA datasets offer us a chance to directly observe the interplay of various evolutionary forces. We show how the genome-wide variance in allele frequency change between two time points can be decomposed into the contributions of gene flow, genetic drift, and linked selection. In closed populations, the contribution of linked selection is identifiable because it creates covariances between time intervals, and genetic drift does not. However, repeated gene flow between populations can also produce directionality in allele frequency change, creating covariances. We show how to accurately separate the fraction of variance in allele frequency change due to admixture and linked selection in a population receiving gene flow. We use two human ancient DNA datasets, spanning around 5,000 years, as time transects to quantify the contributions to the genome-wide variance in allele frequency change. We find that a large fraction of genome-wide change is due to gene flow. In both cases, after correcting for known major gene flow events, we do not observe a signal of genome-wide linked selection. Thus despite the known role of selection in shaping long-term polymorphism levels, and an increasing number of examples of strong selection on single loci and polygenic scores from ancient DNA, it appears to be gene flow and drift, and not selection, that are the main determinants of recent genome-wide allele frequency change. Our approach should be applicable to the growing number of contemporary and ancient temporal population genomics datasets.

An individual with Sarmatian-related ancestry in Roman Britain

In the second century CE the Roman Empire had increasing contact with Sarmatians, nomadic Iranian speakers occupying an area stretching from the Pontic-Caspian steppe to the Carpathian mountains, both in the Caucasus and in the Danubian borders of the empire.1,2,3 In 175 CE, following their defeat in the Marcomannic Wars, emperor Marcus Aurelius drafted Sarmatian cavalry into Roman legions and deployed 5,500 Sarmatian soldiers to Britain, as recorded by contemporary historian Cassius Dio.4,5 Little is known about where the Sarmatian cavalry were stationed, and no individuals connected with this historically attested event have been identified to date, leaving its impact on Britain largely unknown. Here we document Caucasus- and Sarmatian-related ancestry in the whole genome of a Roman-period individual (126-228 calibrated [cal.] CE)-an outlier without traceable ancestry related to local populations in Britain-recovered from a farmstead site in present-day Cambridgeshire, UK. Stable isotopes support a life history of mobility during childhood. Although several scenarios are possible, the historical deployment of Sarmatians to Britain provides a parsimonious explanation for this individual's extraordinary life history. Regardless of the factors behind his migrations, these results highlight how long-range mobility facilitated by the Roman Empire impacted provincial locations outside of urban centers.

Accurate detection of identity-by-descent segments in human ancient DNA

Long DNA segments shared between two individuals, known as identity-by-descent (IBD), reveal recent genealogical connections. Here we introduce ancIBD, a method for identifying IBD segments in ancient human DNA (aDNA) using a hidden Markov model and imputed genotype probabilities. We demonstrate that ancIBD accurately identifies IBD segments >8 cM for aDNA data with an average depth of >0.25× for whole-genome sequencing or >1× for 1240k single nucleotide polymorphism capture data. Applying ancIBD to 4,248 ancient Eurasian individuals, we identify relatives up to the sixth degree and genealogical connections between archaeological groups. Notably, we reveal long IBD sharing between Corded Ware and Yamnaya groups, indicating that the Yamnaya herders of the Pontic-Caspian Steppe and the Steppe-related ancestry in various European Corded Ware groups share substantial co-ancestry within only a few hundred years. These results show that detecting IBD segments can generate powerful insights into the growing aDNA record, both on a small scale relevant to life stories and on a large scale relevant to major cultural-historical events.

The selection landscape and genetic legacy of ancient Eurasians

The Holocene (beginning around 12,000 years ago) encompassed some of the most significant changes in human evolution, with far-reaching consequences for the dietary, physical and mental health of present-day populations. Using a dataset of more than 1,600 imputed ancient genomes1, we modelled the selection landscape during the transition from hunting and gathering, to farming and pastoralism across West Eurasia. We identify key selection signals related to metabolism, including that selection at the FADS cluster began earlier than previously reported and that selection near the LCT locus predates the emergence of the lactase persistence allele by thousands of years. We also find strong selection in the HLA region, possibly due to increased exposure to pathogens during the Bronze Age. Using ancient individuals to infer local ancestry tracts in over 400,000 samples from the UK Biobank, we identify widespread differences in the distribution of Mesolithic, Neolithic and Bronze Age ancestries across Eurasia. By calculating ancestry-specific polygenic risk scores, we show that height differences between Northern and Southern Europe are associated with differential Steppe ancestry, rather than selection, and that risk alleles for mood-related phenotypes are enriched for Neolithic farmer ancestry, whereas risk alleles for diabetes and Alzheimer's disease are enriched for Western hunter-gatherer ancestry. Our results indicate that ancient selection and migration were large contributors to the distribution of phenotypic diversity in present-day Europeans.

Temporal Variation in Introgressed Segments' Length Statistics Computed from a Limited Number of Ancient Genomes Sheds Light on Past Admixture Pulses

Hybridization is recognized as an important evolutionary force, but identifying and timing admixture events between divergent lineages remain a major aim of evolutionary biology. While this has traditionally been done using inferential tools on contemporary genomes, the latest advances in paleogenomics have provided a growing wealth of temporally distributed genomic data. Here, we used individual-based simulations to generate chromosome-level genomic data for a 2-population system and described temporal neutral introgression patterns under a single- and 2-pulse admixture model. We computed 6 summary statistics aiming to inform the timing and number of admixture pulses between interbreeding entities: lengths of introgressed sequences and their variance within genomes, as well as genome-wide introgression proportions and related measures. The first 2 statistics could confidently be used to infer interlineage hybridization history, peaking at the beginning and shortly after an admixture pulse. Temporal variation in introgression proportions and related statistics provided more limited insights, particularly when considering their application to ancient genomes still scant in number. Lastly, we computed these statistics on Homo sapiens paleogenomes and successfully inferred the hybridization pulse from Neanderthal that occurred approximately 40 to 60 kya. The scarce number of genomes dating from this period prevented more precise inferences, but the accumulation of paleogenomic data opens promising perspectives as our approach only requires a limited number of ancient genomes.

Allelic bias when performing in-solution enrichment of ancient human DNA

In-solution hybridisation enrichment of genetic variation is a valuable methodology in human paleogenomics. It allows enrichment of endogenous DNA by targeting genetic markers that are comparable between sequencing libraries. Many studies have used the 1240k reagent-which enriches 1,237,207 genome-wide SNPs-since 2015, though access was restricted. In 2021, Twist Biosciences and Daicel Arbor Biosciences independently released commercial kits that enabled all researchers to perform enrichments for the same 1240 k SNPs. We used the Daicel Arbor Biosciences Prime Plus kit to enrich 132 ancient samples from three continents. We identified a systematic assay bias that increases genetic similarity between enriched samples and that cannot be explained by batch effects. We present the impact of the bias on population genetics inferences (e.g. Principal Components Analysis, ƒ-statistics) and genetic relatedness (READ). We compare the Prime Plus bias to that previously reported of the legacy 1240k enrichment assay. In ƒ-statistics, we find that all Prime-Plus-generated data exhibit artefactual excess shared drift, such that within-continent relationships cannot be correctly determined. The bias is more subtle in READ, though interpretation of the results can still be misleading in specific contexts. We expect the bias may affect analyses we have not yet tested. Our observations support previously reported concerns for the integration of different data types in paleogenomics. We also caution that technological solutions to generate 1240k data necessitate a thorough validation process before their adoption in the paleogenomic community.

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.

Are vipers prototypic fear-evoking snakes? A cross-cultural comparison of Somalis and Czechs

Snakes are known as highly fear-evoking animals, eliciting preferential attention and fast detection in humans. We examined the human fear response to snakes in the context of both current and evolutionary experiences, conducting our research in the cradle of humankind, the Horn of Africa. This region is characterized by the frequent occurrence of various snake species, including deadly venomous viperids (adders) and elapids (cobras and mambas). We conducted experiments in Somaliland and compared the results with data from Czech respondents to address the still unresolved questions: To which extent is human fear of snakes affected by evolutionary or current experience and local culture? Can people of both nationalities recognize venomous snakes as a category, or are they only afraid of certain species that are most dangerous in a given area? Are respondents of both nationalities equally afraid of deadly snakes from both families (Viperidae, Elapidae)? We employed a well-established picture-sorting approach, consisting of 48 snake species belonging to four distinct groups. Our results revealed significant agreement among Somali as well as Czech respondents. We found a highly significant effect of the stimulus on perceived fear in both populations. Vipers appeared to be the most salient stimuli in both populations, as they occupied the highest positions according to the reported level of subjectively perceived fear. The position of vipers strongly contrasts with the fear ranking of deadly venomous elapids, which were in lower positions. Fear scores of vipers were significantly higher in both populations, and their best predictor was the body width of the snake. The evolutionary, cultural, and cognitive aspects of this phenomenon are discussed.

Predicting functional consequences of recent natural selection in Britain

Ancient DNA can directly reveal the contribution of natural selection to human genomic variation. However, while the analysis of ancient DNA has been successful at identifying genomic signals of selection, inferring the phenotypic consequences of that selection has been more difficult. Most trait-associated variants are non-coding, so we expect that a large proportion of the phenotypic effects of selection will also act through non-coding variation. Since we cannot measure gene expression directly in ancient individuals, we used an approach (Joint-Tissue Imputation; JTI) developed to predict gene expression from genotype data. We tested for changes in the predicted expression of 17,384 protein coding genes over a time transect of 4500 years using 91 present-day and 616 ancient individuals from Britain. We identified 28 genes at seven genomic loci with significant (FDR < 0.05) changes in predicted expression levels in this time period. We compared the results from our transcriptome-wide scan to a genome-wide scan based on estimating per-SNP selection coefficients from time series data. At five previously identified loci, our approach allowed us to highlight small numbers of genes with evidence for significant shifts in expression from peaks that in some cases span tens of genes. At two novel loci (SLC44A5 and NUP85), we identify selection on gene expression not captured by scans based on genomic signatures of selection. Finally we show how classical selection statistics (iHS and SDS) can be combined with JTI models to incorporate functional information into scans that use present-day data alone. These results demonstrate the potential of this type of information to explore both the causes and consequences of natural selection.

Multivariate assessments of activity-related skeletal changes: Interpreting Bell Beaker specialized male archery and social organization in Central Europe

OBJECTIVES

The Bell Beaker period witnessed the rise of individual inhumations with "wealthy" burial contexts containing archery-related grave goods, leading archaeologists to label the individuals in these tombs as "archers." This study looks to (1) compare the skeletons from male "archer" burials with those from male "non-archer" burials-those not having archery-related grave goods-in order to assess a possible link between burial context and physical activity, and (2) apply a biomechanics profile to evaluate whether the individuals associated with these "archer" burials practiced specialized archer activity.

MATERIALS AND METHODS

The corpus (males only) included 46 "archers" and 40 "non-archers" from Bell Beaker individual inhumations. Osteological data included measurements, scores of entheseal changes, and a diagnosis of certain pathologies. Data analyses involved visual observations, hypothesis tests, dimension reduction, and MANOVA, with approaches aimed at exploring the treatment of data missingness.

RESULTS

Measurement data revealed no differences between the two groups. Evaluations of entheseal changes found that "non-archers" had consistently more instances of bone surface modifications than "archers." Individual assessments of specialized archer occupation identified 11 possible specialized archers.

DISCUSSION

These findings indicate a possible labor differentiation represented through the presence of a probably prestigious "archer" burial context. This suggests a link between grave good presence and labor, but not between a Bell Beaker archery occupation and an "archer" burial context. Data analyses support the application of biomechanics to osteological analyses in order to assess specialized activity on the skeleton.

High-coverage genome of the Tyrolean Iceman reveals unusually high Anatolian farmer ancestry

The Tyrolean Iceman is known as one of the oldest human glacier mummies, directly dated to 3350-3120 calibrated BCE. A previously published low-coverage genome provided novel insights into European prehistory, despite high present-day DNA contamination. Here, we generate a high-coverage genome with low contamination (15.3×) to gain further insights into the genetic history and phenotype of this individual. Contrary to previous studies, we found no detectable Steppe-related ancestry in the Iceman. Instead, he retained the highest Anatolian-farmer-related ancestry among contemporaneous European populations, indicating a rather isolated Alpine population with limited gene flow from hunter-gatherer-ancestry-related populations. Phenotypic analysis revealed that the Iceman likely had darker skin than present-day Europeans and carried risk alleles associated with male-pattern baldness, type 2 diabetes, and obesity-related metabolic syndrome. These results corroborate phenotypic observations of the preserved mummified body, such as high pigmentation of his skin and the absence of hair on his head.

Interdisciplinary Analyses of Bronze Age Communities from Western Hungary Reveal Complex Population Histories

In this study, we report 21 ancient shotgun genomes from present-day Western Hungary, from previously understudied Late Copper Age Baden, and Bronze Age Somogyvár-Vinkovci, Kisapostag, and Encrusted Pottery archeological cultures (3,530-1,620 cal Bce). Our results indicate the presence of high steppe ancestry in the Somogyvár-Vinkovci culture. They were then replaced by the Kisapostag group, who exhibit an outstandingly high (up to ∼47%) Mesolithic hunter-gatherer ancestry, despite this component being thought to be highly diluted by the time of the Early Bronze Age. The Kisapostag population contributed the genetic basis for the succeeding community of the Encrusted Pottery culture. We also found an elevated hunter-gatherer component in a local Baden culture-associated individual, but no connections were proven to the Bronze Age individuals. The hunter-gatherer ancestry in Kisapostag is likely derived from two main sources, one from a Funnelbeaker or Globular Amphora culture-related population and one from a previously unrecognized source in Eastern Europe. We show that this ancestry not only appeared in various groups in Bronze Age Central Europe but also made contributions to Baltic populations. The social structure of Kisapostag and Encrusted Pottery cultures is patrilocal, similarly to most contemporaneous groups. Furthermore, we developed new methods and method standards for computational analyses of ancient DNA, implemented to our newly developed and freely available bioinformatic package. By analyzing clinical traits, we found carriers of aneuploidy and inheritable genetic diseases. Finally, based on genetic and anthropological data, we present here the first female facial reconstruction from the Bronze Age Carpathian Basin.

The genetic legacy of African Americans from Catoctin Furnace

Few African Americans have been able to trace family lineages back to ancestors who died before the 1870 United States Census, the first in which all Black people were listed by name. We analyzed 27 individuals from Maryland's Catoctin Furnace African American Cemetery (1774-1850), identifying 41,799 genetic relatives among consenting research participants in 23andMe, Inc.'s genetic database. One of the highest concentrations of close relatives is in Maryland, suggesting that descendants of the Catoctin individuals remain in the area. We find that many of the Catoctin individuals derived African ancestry from the Wolof or Kongo groups and European ancestry from Great Britain and Ireland. This study demonstrates the power of joint analysis of historical DNA and large datasets generated through direct-to-consumer ancestry testing.

Early contact between late farming and pastoralist societies in southeastern Europe

Archaeogenetic studies have described two main genetic turnover events in prehistoric western Eurasia: one associated with the spread of farming and a sedentary lifestyle starting around 7000-6000 BC (refs. 1-3) and a second with the expansion of pastoralist groups from the Eurasian steppes starting around 3300 BC (refs. 4,5). The period between these events saw new economies emerging on the basis of key innovations, including metallurgy, wheel and wagon and horse domestication6-9. However, what happened between the demise of the Copper Age settlements around 4250 BC and the expansion of pastoralists remains poorly understood. To address this question, we analysed genome-wide data from 135 ancient individuals from the contact zone between southeastern Europe and the northwestern Black Sea region spanning this critical time period. While we observe genetic continuity between Neolithic and Copper Age groups from major sites in the same region, from around 4500 BC on, groups from the northwestern Black Sea region carried varying amounts of mixed ancestries derived from Copper Age groups and those from the forest/steppe zones, indicating genetic and cultural contact over a period of around 1,000 years earlier than anticipated. We propose that the transfer of critical innovations between farmers and transitional foragers/herders from different ecogeographic zones during this early contact was integral to the formation, rise and expansion of pastoralist groups around 3300 BC.

Patrilocality and hunter-gatherer-related ancestry of populations in East-Central Europe during the Middle Bronze Age

The demographic history of East-Central Europe after the Neolithic period remains poorly explored, despite this region being on the confluence of various ecological zones and cultural entities. Here, the descendants of societies associated with steppe pastoralists form Early Bronze Age were followed by Middle Bronze Age populations displaying unique characteristics. Particularly, the predominance of collective burials, the scale of which, was previously seen only in the Neolithic. The extent to which this re-emergence of older traditions is a result of genetic shift or social changes in the MBA is a subject of debate. Here by analysing 91 newly generated genomes from Bronze Age individuals from present Poland and Ukraine, we discovered that Middle Bronze Age populations were formed by an additional admixture event involving a population with relatively high proportions of genetic component associated with European hunter-gatherers and that their social structure was based on, primarily patrilocal, multigenerational kin-groups.

Language trees with sampled ancestors support a hybrid model for the origin of Indo-European languages

The origins of the Indo-European language family are hotly disputed. Bayesian phylogenetic analyses of core vocabulary have produced conflicting results, with some supporting a farming expansion out of Anatolia ~9000 years before present (yr B.P.), while others support a spread with horse-based pastoralism out of the Pontic-Caspian Steppe ~6000 yr B.P. Here we present an extensive database of Indo-European core vocabulary that eliminates past inconsistencies in cognate coding. Ancestry-enabled phylogenetic analysis of this dataset indicates that few ancient languages are direct ancestors of modern clades and produces a root age of ~8120 yr B.P. for the family. Although this date is not consistent with the Steppe hypothesis, it does not rule out an initial homeland south of the Caucasus, with a subsequent branch northward onto the steppe and then across Europe. We reconcile this hybrid hypothesis with recently published ancient DNA evidence from the steppe and the northern Fertile Crescent.

Explaining population booms and busts in Mid-Holocene Europe

Archaeological evidence suggests that the population dynamics of Mid-Holocene (Late Mesolithic to Initial Bronze Age, ca. 7000-3000 BCE) Europe are characterized by recurrent booms and busts of regional settlement and occupation density. These boom-bust patterns are documented in the temporal distribution of 14C dates and in archaeological settlement data from regional studies. We test two competing hypotheses attempting to explain these dynamics: climate forcing and social dynamics leading to inter-group conflict. Using the framework of spatially-explicit agent-based models, we translated these hypotheses into a suite of explicit computational models, derived quantitative predictions for population fluctuations, and compared these predictions to data. We demonstrate that climate variation during the European Mid-Holocene is unable to explain the quantitative features (average periodicities and amplitudes) of observed boom-bust dynamics. In contrast, scenarios with social dynamics encompassing density-dependent conflict produce population patterns with time scales and amplitudes similar to those observed in the data. These results suggest that social processes, including violent conflict, played a crucial role in the shaping of population dynamics of European Mid-Holocene societies.

Pilot study of correlation of selected genetic factors with cribra orbitalia in individuals from a medieval population from Slovakia

OBJECTIVE

The aim of this study is to investigate the potential genetic etiology of cribra orbitalia noted on human skeletal remains.

MATERIALS

We obtained and analyzed ancient DNA of 43 individuals with cribra orbitalia. The analyzed set represented medieval individuals from two cemeteries in western Slovakia, Castle Devín (11th-12th century AD) and Cífer-Pác (8th-9th century AD).

METHODS

We performed a sequence analysis of 5 variants in 3 genes associated with anemia (HBB, G6PD, PKLR), which are the most common pathogenic variants in present day of European populations, and one variant MCM6:c.1917 + 326 C>T (rs4988235) associated with lactose intolerance.

RESULTS

DNA variants associated with anemia were not found in the samples. The allele frequency of MCM6:c.1917 + 326 C was 0.875. This frequency is higher but not statistically significant in individuals displaying cribra orbitalia compared to individuals without the lesion.

SIGNIFICANCE

This study seeks to expand our knowledge of the etiology of cribra orbitalia by exploring the potential association between the lesion and the presence of alleles linked to hereditary anemias and lactose intolerance.

LIMITATIONS

A relatively small set of individuals were analyzed, so an unequivocal conclusion cannot be drawn. Hence, although it is unlikely, a genetic form of anemia caused by rare variants cannot be ruled out.

SUGGESTIONS FOR FURTHER RESEARCH

Genetic research based on larger sample sizes and in more diverse geographical regions.

Yersinia pestis genomes reveal plague in Britain 4000 years ago

Extinct lineages of Yersinia pestis, the causative agent of the plague, have been identified in several individuals from Eurasia between 5000 and 2500 years before present (BP). One of these, termed the 'LNBA lineage' (Late Neolithic and Bronze Age), has been suggested to have spread into Europe with human groups expanding from the Eurasian steppe. Here, we show that the LNBA plague was spread to Europe's northwestern periphery by sequencing three Yersinia pestis genomes from Britain, all dating to ~4000 cal BP. Two individuals were from an unusual mass burial context in Charterhouse Warren, Somerset, and one individual was from a single burial under a ring cairn monument in Levens, Cumbria. To our knowledge, this represents the earliest evidence of LNBA plague in Britain documented to date. All three British Yersinia pestis genomes belong to a sublineage previously observed in Bronze Age individuals from Central Europe that had lost the putative virulence factor yapC. This sublineage is later found in Eastern Asia ~3200 cal BP. While the severity of the disease is currently unclear, the wide geographic distribution within a few centuries suggests substantial transmissibility.

Genetic predictors of cultural values variation between societies

Associations between the STin2 and 5-HTTLPR polymorphisms within the serotonin transporter gene, SLC6A4, and culture across societies were examined. Based on an analysis of 75 primary studies (28,726 individuals), STin2 allelic frequencies were found to vary widely across countries, ranging from 26% in Germany to 85% in Singapore. Across 53 countries, and after controlling for all major environmental influences of culture, STin2 and 5-HTTLPR were found to explain 23.6% unique variance in monumentalism but none in individualism. Our findings evidence a significant role of genetics in predicting cross-societal cultural values variation, and potentially speak to the need for and importance of incorporating both nature and nurture in theories of cultural values variation across societies.

Understanding natural selection in Holocene Europe using multi-locus genotype identity scans

Ancient DNA (aDNA) has been a revolutionary technology in understanding human history but has not been used extensively to study natural selection as large sample sizes to study allele frequency changes over time have thus far not been available. Here, we examined a time transect of 708 published samples over the past 7,000 years of European history using multi-locus genotype-based selection scans. As aDNA data is affected by high missingness, ascertainment bias, DNA damage, random allele calling, and is unphased, we first validated our selection scan, G 12 a n c i e n t , on simulated data resembling aDNA under a demographic model that captures broad features of the allele frequency spectrum of European genomes as well as positive controls that have been previously identified and functionally validated in modern European datasets on data from ancient individuals from time periods very close to the present time. We then applied our statistic to the aDNA time transect to detect and resolve the timing of natural selection occurring genome wide and found several candidates of selection across the different time periods that had not been picked up by selection scans using single SNP allele frequency approaches. In addition, enrichment analysis discovered multiple categories of complex traits that might be under adaptation across these periods. Our results demonstrate the utility of applying different types of selection scans to aDNA to uncover putative selection signals at loci in the ancient past that might have been masked in modern samples.

Clinical, genetic, epidemiologic, evolutionary, and functional delineation of TSPEAR-related autosomal recessive ectodermal dysplasia 14

TSPEAR variants cause autosomal recessive ectodermal dysplasia (ARED) 14. The function of TSPEAR is unknown. The clinical features, the mutation spectrum, and the underlying mechanisms of ARED14 are poorly understood. Combining data from new and previously published individuals established that ARED14 is primarily characterized by dental anomalies such as conical tooth cusps and hypodontia, like those seen in individuals with WNT10A-related odontoonychodermal dysplasia. AlphaFold-predicted structure-based analysis showed that most of the pathogenic TSPEAR missense variants likely destabilize the β-propeller of the protein. Analysis of 100000 Genomes Project (100KGP) data revealed multiple founder TSPEAR variants across different populations. Mutational and recombination clock analyses demonstrated that non-Finnish European founder variants likely originated around the end of the last ice age, a period of major climatic transition. Analysis of gnomAD data showed that the non-Finnish European population TSPEAR gene-carrier rate is ∼1/140, making it one of the commonest AREDs. Phylogenetic and AlphaFold structural analyses showed that TSPEAR is an ortholog of drosophila Closca, an extracellular matrix-dependent signaling regulator. We, therefore, hypothesized that TSPEAR could have a role in enamel knot, a structure that coordinates patterning of developing tooth cusps. Analysis of mouse single-cell RNA sequencing (scRNA-seq) data revealed highly restricted expression of Tspear in clusters representing enamel knots. A tspeara -/-;tspearb -/- double-knockout zebrafish model recapitulated the clinical features of ARED14 and fin regeneration abnormalities of wnt10a knockout fish, thus suggesting interaction between tspear and wnt10a. In summary, we provide insights into the role of TSPEAR in ectodermal development and the evolutionary history, epidemiology, mechanisms, and consequences of its loss of function variants.

A 23,000-year-old southern Iberian individual links human groups that lived in Western Europe before and after the Last Glacial Maximum

Human populations underwent range contractions during the Last Glacial Maximum (LGM) which had lasting and dramatic effects on their genetic variation. The genetic ancestry of individuals associated with the post-LGM Magdalenian technocomplex has been interpreted as being derived from groups associated with the pre-LGM Aurignacian. However, both these ancestries differ from that of central European individuals associated with the chronologically intermediate Gravettian. Thus, the genomic transition from pre- to post-LGM remains unclear also in western Europe, where we lack genomic data associated with the intermediate Solutrean, which spans the height of the LGM. Here we present genome-wide data from sites in Andalusia in southern Spain, including from a Solutrean-associated individual from Cueva del Malalmuerzo, directly dated to ~23,000 cal yr BP. The Malalmuerzo individual carried genetic ancestry that directly connects earlier Aurignacian-associated individuals with post-LGM Magdalenian-associated ancestry in western Europe. This scenario differs from Italy, where individuals associated with the transition from pre- and post-LGM carry different genetic ancestries. This suggests different dynamics in the proposed southern refugia of Ice Age Europe and posits Iberia as a potential refugium for western European pre-LGM ancestry. More, individuals from Cueva Ardales, which were thought to be of Palaeolithic origin, date younger than expected and, together with individuals from the Andalusian sites Caserones and Aguilillas, fall within the genetic variation of the Neolithic, Chalcolithic and Bronze Age individuals from southern Iberia.

Inferring biological kinship in ancient datasets: comparing the response of ancient DNA-specific software packages to low coverage data

BACKGROUND

The inference of biological relations between individuals is fundamental to understanding past human societies. Caregiving, resource sharing and sexual behaviours are often mediated by biological kinship and yet the identification and interpretation of kin relationships in prehistoric human groups is difficult. In recent years, the advent of archaeogenetic techniques have offered a fresh approach, and when combined with more traditional osteological and interpretive archaeological methods, allows for improved interpretation of the burial practices, cultural behaviours, and societal stratification in ancient societies. Although archaeogenetic techniques are developing at pace, questions remain as to their accuracy, particularly when applied to the low coverage datasets that results from the sequencing of DNA derived from highly degraded ancient material.

RESULTS

The performance of six of the most commonly used kinship identifcation software methods was explored at a range of low and ultra low genome coverages. An asymmetrical response was observed across packages, with decreased genome coverage resulting in differences in both direction and degree of change of calculated kinship scores and thus pairwise relatedness estimates are dependant on both package used and genome coverage. Methods reliant upon genotype likelihoods methods (lcMLkin, NGSrelate and NGSremix) show a decreased level of prediction at coverage below 1x, although were consistent in the particular relationships identified at these coverages when compared to the pseudohaploid reliant methods tested (READ, the Kennett 2017 method and TKGWV2.0). The three pseudohaploid methods show predictive potential at coverages as low as 0.05x, although the accuracy of the relationships identified is questionable given the increase in the number of relationships identifIed at the low coverage (type I errors).

CONCLUSION

Two pseudohaploid methods (READ and Kennett 2017) show relatively consistent inference of kin relationships at low coverage (0.5x), with READ only showing a significant performance drop off at ultralow coverages (< 0.2x). More generally, our results reveal asymmetrical kinship classifications in some software packages even at high coverages, highlighting the importance of applying multiple methods to authenticate kin relationships in ancient material, along with the continuing need to develop laboratory methods that maximise data output for downstream analyses.

The impacts of bronze age in the gene pool of Chinese: Insights from phylogeographics of Y-chromosomal haplogroup N1a2a-F1101

Objectives: Previous studies of archaeology and history suggested that the rise and prosperity of Bronze Age culture in East Asia had made essential contribution to the formation of early state and civilization in this region. However, the impacts in perspective of genetics remain ambiguous. Previous genetic researches indicated the Y-chromosome Q1a1a-M120 and N1a2a-F1101 may be the two most important paternal lineages among the Bronze Age people in ancient northwest China. Here, we investigated the 9,000-years history of haplogroup N1a2a-F1101 with revised phylogenetic tree and spatial autocorrelation analysis.

Materials and Methods: In this study, 229 sequences of N1a2a-F1101 were analyzed. We developed a highly-revised phylogenetic tree with age estimates for N1a2a-F1101. In addition, we also explored the geographical distribution of sub-lineages of N1a2a-F1101, and spatial autocorrelation analysis was conducted for each sub-branch.

Results: The initial differentiation location of N1a2a-F1101 and its most closely related branch, N1a2b-P43, a major lineage of Uralic-speaking populations in northern Eurasia, is likely the west part of northeast China. After ~4 thousand years of bottleneck effect period, haplgroup N1a2a-F1101 experienced continuous expansion during the Chalcolithic age (~ 4.5 kya to 4 kya) and Bronze age (~ 4 kya to 2.5 kya) in northern China. Ancient DNA evidence supported that this haplogroup is the lineage of ruling family of Zhou Dynasty (~ 3 kya-2.2 kya) of ancient China.

Discussion: In general, we proposed that the Bronze Age people in the border area between the eastern Eurasian steppe and northern China not only played a key role in promoting the early state and civilization of China, but also left significant traces in the gene pool of Chinese people.

ancIBD - Screening for identity by descent segments in human ancient DNA

Long DNA sequences shared between two individuals, known as Identical by descent (IBD) segments, are a powerful signal for identifying close and distant biological relatives because they only arise when the pair shares a recent common ancestor. Existing methods to call IBD segments between present-day genomes cannot be straightforwardly applied to ancient DNA data (aDNA) due to typically low coverage and high genotyping error rates. We present ancIBD, a method to identify IBD segments for human aDNA data implemented as a Python package. Our approach is based on a Hidden Markov Model, using as input genotype probabilities imputed based on a modern reference panel of genomic variation. Through simulation and downsampling experiments, we demonstrate that ancIBD robustly identifies IBD segments longer than 8 centimorgan for aDNA data with at least either 0.25x average whole-genome sequencing (WGS) coverage depth or at least 1x average depth for in-solution enrichment experiments targeting a widely used aDNA SNP set ('1240k'). This application range allows us to screen a substantial fraction of the aDNA record for IBD segments and we showcase two downstream applications. First, leveraging the fact that biological relatives up to the sixth degree are expected to share multiple long IBD segments, we identify relatives between 10,156 ancient Eurasian individuals and document evidence of long-distance migration, for example by identifying a pair of two approximately fifth-degree relatives who were buried 1410km apart in Central Asia 5000 years ago. Second, by applying ancIBD, we reveal new details regarding the spread of ancestry related to Steppe pastoralists into Europe starting 5000 years ago. We find that the first individuals in Central and Northern Europe carrying high amounts of Steppe-ancestry, associated with the Corded Ware culture, share high rates of long IBD (12-25 cM) with Yamnaya herders of the Pontic-Caspian steppe, signaling a strong bottleneck and a recent biological connection on the order of only few hundred years, providing evidence that the Yamnaya themselves are a main source of Steppe ancestry in Corded Ware people. We also detect elevated sharing of long IBD segments between Corded Ware individuals and people associated with the Globular Amphora culture (GAC) from Poland and Ukraine, who were Copper Age farmers not yet carrying Steppe-like ancestry. These IBD links appear for all Corded Ware groups in our analysis, indicating that individuals related to GAC contexts must have had a major demographic impact early on in the genetic admixtures giving rise to various Corded Ware groups across Europe. These results show that detecting IBD segments in aDNA can generate new insights both on a small scale, relevant to understanding the life stories of people, and on the macroscale, relevant to large-scale cultural-historical events.

Hunter-gatherer admixture facilitated natural selection in Neolithic European farmers

Ancient DNA has revealed multiple episodes of admixture in human prehistory during geographic expansions associated with cultural innovations. One important example is the expansion of Neolithic agricultural groups out of the Near East into Europe and their consequent admixture with Mesolithic hunter-gatherers.^1,2,3,4 Ancient genomes from this period provide an opportunity to study the role of admixture in providing new genetic variation for selection to act upon, and also to identify genomic regions that resisted hunter-gatherer introgression and may thus have contributed to agricultural adaptations. We used genome-wide DNA from 677 individuals spanning Mesolithic and Neolithic Europe to infer ancestry deviations in the genomes of admixed individuals and to test for natural selection after admixture by testing for deviations from a genome-wide null distribution. We find that the region around the pigmentation-associated gene SLC24A5 shows the greatest overrepresentation of Neolithic local ancestry in the genome (|Z| = 3.46). In contrast, we find the greatest overrepresentation of Mesolithic ancestry across the major histocompatibility complex (MHC; |Z| = 4.21), a major immunity locus, which also shows allele frequency deviations indicative of selection following admixture (p = 1 × 10^-56). This could reflect negative frequency-dependent selection on MHC alleles common in Neolithic populations or that Mesolithic alleles were positively selected for and facilitated adaptation in Neolithic populations to pathogens or other environmental factors. Our study extends previous results that highlight immune function and pigmentation as targets of adaptation in more recent populations to selection processes in the Stone Age.

Estimating human mobility in Holocene Western Eurasia with large-scale ancient genomic data

The recent increase in openly available ancient human DNA samples allows for large-scale meta-analysis applications. Trans-generational past human mobility is one of the key aspects that ancient genomics can contribute to since changes in genetic ancestry-unlike cultural changes seen in the archaeological record-necessarily reflect movements of people. Here, we present an algorithm for spatiotemporal mapping of genetic profiles, which allow for direct estimates of past human mobility from large ancient genomic datasets. The key idea of the method is to derive a spatial probability surface of genetic similarity for each individual in its respective past. This is achieved by first creating an interpolated ancestry field through space and time based on multivariate statistics and Gaussian process regression and then using this field to map the ancient individuals into space according to their genetic profile. We apply this algorithm to a dataset of 3138 aDNA samples with genome-wide data from Western Eurasia in the last 10,000 y. Finally, we condense this sample-wise record with a simple summary statistic into a diachronic measure of mobility for subregions in Western, Central, and Southern Europe. For regions and periods with sufficient data coverage, our similarity surfaces and mobility estimates show general concordance with previous results and provide a meta-perspective of genetic changes and human mobility.

Nostalgia in the prehistoric archaeological record

Evidence from the prehistoric archaeological record clearly shows that ancient societies had a sense of and engaged with their own histories, be it by reusing, re-appropriating or recreating past material culture. The affective qualities of materials, places and even human remains would have enabled people to remember and connect with aspects of their immediate and more distant pasts. In some cases, this may have elicited specific emotive responses, similar to how nostalgic triggers operate today. Nostalgia is not a word commonly used by archaeologists, but through exploring the materiality and sensory affect of objects and spaces in the past, we can consider that what we deal with archaeologically may have held nostalgic qualities.

Health inequity in genomic personalized medicine in underrepresented populations: a look at the current evidence

Improvements in sequencing technology coupled with dramatic declines in the cost of genome sequencing have led to a proportional growth in the size and number of genetic datasets since the release of the human genetic sequence by The Human Genome Project (HGP) international consortium. The HGP was undeniably a significant scientific success, a turning point in human genetics and the beginning of human genomics. This burst of genetic information has led to a greater understanding of disease pathology and the potential of employing this data to deliver more precise patient care. Hence, the recognition of high-penetrance disease-causing mutations which encode drivers of disease has made the management of most diseases more specific. Nonetheless, while genetic scores are becoming more extensively used, their application in the real world is expected to be limited due to the lack of diversity in the data used to construct them. Underrepresented populations, such as racial and ethnic minorities, low-income individuals, and those living in rural areas, often experience greater health disparities and worse health outcomes compared to the general population. These disparities are often the result of systemic barriers, such as poverty, discrimination, and limited access to healthcare. Addressing health inequity in underrepresented populations requires addressing the underlying social determinants of health and implementing policies and programs which promoted health equity and reduce disparities. This can include expanding access to affordable healthcare, addressing poverty and unemployment, and promoting policies that combat discrimination and racism.

Genomic perspectives on human dispersals during the Holocene

Nearly 20 y ago, Jared Diamond and Peter Bellwood reviewed the evidence for the associated spread of farming and large language families by the demographic expansions of farmers. Since then, advances in obtaining and analyzing genomic data from modern and ancient populations have transformed our knowledge of human dispersals during the Holocene. Here, we provide an overview of Holocene dispersals in the light of genomic evidence and conclude that they have a complex history. Even when there is a demonstrated connection between a demographic expansion of people, the spread of agriculture, and the spread of a particular language family, the outcome in the results of contact between expanding and resident groups is highly variable. Further research is needed to identify the factors and social circumstances that have influenced this variation and complex history.

Genetic admixture and language shift in the medieval Volga-Oka interfluve

The Volga-Oka interfluve in northwestern Russia has an intriguing history of population influx and language shift during the Common Era. Today, most inhabitants of the region speak Russian, but until medieval times, northwestern Russia was inhabited by Uralic-speaking peoples.^1,2,3 A gradual shift to Slavic languages started in the second half of the first millennium with the expansion of Slavic tribes, which led to the foundation of the Kievan Rus' state in the late 9^th century CE. The medieval Rus' was multicultural and multilingual-historical records suggest that its northern regions comprised Slavic and Uralic peoples ruled by Scandinavian settlers.^4,5,6 In the 10^th-11^th centuries, the introduction of Christianity and Cyrillic literature raised the prestige status of Slavic, driving a language shift from Uralic to Slavic.³ This eventually led to the disappearance of the Uralic languages from northwestern Russia. Here, we study a 1,500-year time transect of 30 ancient genomes and stable isotope values from the Suzdal region in the Volga-Oka interfluve. We describe a previously unsampled local Iron Age population and a gradual genetic turnover in the following centuries. Our time transect captures the population shift associated with the spread of Slavic languages and illustrates the ethnically mixed state of medieval Suzdal principality, eventually leading to the formation of the admixed but fully Slavic-speaking population that inhabits the area today. We also observe genetic outliers that highlight the importance of the Suzdal region in medieval times as a hub of long-reaching contacts via trade and warfare.

Neolithic expansion and the 17q21.31 inversion in Iberia: an evolutionary approach to H2 haplotype distribution in the Near East and Europe

The chromosomal region 17q21.31 harbors a 900 kb inversion polymorphism named after the microtubule-associated protein tau (MAPT) gene. Since no recombination occurs, two haplotypes are recognized: a directly oriented variant (H1) and an inverted variant (H2). The H2 haplotype features a distribution pattern with high frequencies in the Near East and Europe, medium levels in South Asia and North Africa, and low levels elsewhere. Studies of this genomic region are relevant owing to its likely association with numerous neurodegenerative diseases. However, the causes underlying the geographic distribution of the haplotype frequencies remain a bone of contention among researchers. With this work, we have intended to outline a plausible hypothesis on the origin of the high European H2 frequencies. To that end, we have analyzed an extensive population database (including three new Iberian populations) to explore potential clinal variations of H2 frequencies. We found a sigmoidal frequency cline with an upward trend from South Asia to Europe. The maximum value was detected in the Basques from Gipuzkoa province (0.494) with the curve's inflection point in the Near East. From our results, we suggest that the most likely scenario for high H2 frequencies in Europe would be a founding event in the Near East during the late Paleolithic or early Neolithic. Subsequently, such H2 overrepresentation would have reached Europe with the arrival of the first Neolithic farmers. The current frequencies and geographic distribution of the 17q21.31 inversion suggest that the founding events mainly affected the H2D subhaplotype.

Sex bias in Neolithic megalithic burials

Objectives

A statistical study comparing osteological and ancient DNA determinations of sex was conducted in order to investigate whether there are sex biases in United Kingdom and Irish Neolithic megalithic burials.

Materials and Methods

Genetic and osteological information from human individuals from 32 megalithic sites in the UK and Ireland dating from 4000 to 2500 cal. BCE was collected and statistically analyzed to test whether there is a true over‐representation of males at these sites. The published dataset from the study by Sánchez‐Quinto et al. in 2019 was initially analyzed before being refined and included in a larger dataset. Osteological analysis of sex bias was limited to adults with available sex estimations, and genetic analysis limited to published data

Results

Two sites consistently returned significant p ‐values suggesting a potential over‐representation in osteological males at one site (Knowe of Midhowe, Orkney) and genetic males in the other (Primrose Grange, Ireland). Cumulative statistical analyses point towards a male bias in the representation of sexes in Neolithic megalithic burials, but these results do not reflect the site‐by‐site and regional variation found in this study.

Discussion

The interpretation of sex bias, that is, the over‐representation of one sex over another ‐ depends on other socio‐cultural variables (e.g., kinship) and the emphasis placed on statistical significance. The trend towards males being over‐represented in Neolithic megalithic burials is not as clear as previously thought, and requires further testing and data collection to uncover.

Ancient DNA sequencing: telling the tale of human history and evolution

Standfirst: In the last decade, ancient DNA research has provided invaluable insights into the lives of ancient populations, adding detail and enriching the story of human evolution and development. [Formula: see text].

Ancient DNA reveals admixture history and endogamy in the prehistoric Aegean

The Neolithic and Bronze Ages were highly transformative periods for the genetic history of Europe but for the Aegean-a region fundamental to Europe's prehistory-the biological dimensions of cultural transitions have been elucidated only to a limited extent so far. We have analysed newly generated genome-wide data from 102 ancient individuals from Crete, the Greek mainland and the Aegean Islands, spanning from the Neolithic to the Iron Age. We found that the early farmers from Crete shared the same ancestry as other contemporaneous Neolithic Aegeans. In contrast, the end of the Neolithic period and the following Early Bronze Age were marked by 'eastern' gene flow, which was predominantly of Anatolian origin in Crete. Confirming previous findings for additional Central/Eastern European ancestry in the Greek mainland by the Middle Bronze Age, we additionally show that such genetic signatures appeared in Crete gradually from the seventeenth to twelfth centuries BC, a period when the influence of the mainland over the island intensified. Biological and cultural connectedness within the Aegean is also supported by the finding of consanguineous endogamy practiced at high frequencies, unprecedented in the global ancient DNA record. Our results highlight the potential of archaeogenomic approaches in the Aegean for unravelling the interplay of genetic admixture, marital and other cultural practices.

Genomic analyses reveal range-wide devastation of sea otter populations

The genetic consequences of species-wide declines are rarely quantified because the timing and extent of the decline varies across the species' range. The sea otter (Enhydra lutris) is a unique model in this regard. Their dramatic decline from thousands to fewer than 100 individuals per population occurred range-wide and nearly simultaneously due to the 18th-19th century fur trade. Consequently, each sea otter population represents an independent natural experiment of recovery after extreme population decline. We designed sequence capture probes for 50 Mb of sea otter exonic and neutral genomic regions. We sequenced 107 sea otters from five populations that span the species range to high coverage (18-76×) and three historical Californian samples from ~1500 and ~200 years ago to low coverage (1.5-3.5×). We observe distinct population structure and find that sea otters in California are the last survivors of a divergent lineage isolated for thousands of years and therefore warrant special conservation concern. We detect signals of extreme population decline in every surviving sea otter population and use this demographic history to design forward-in-time simulations of coding sequence. Our simulations indicate that this decline could lower the fitness of recovering populations for generations. However, the simulations also demonstrate how historically low effective population sizes prior to the fur trade may have mitigated the effects of population decline on genetic health. Our comprehensive approach shows how demographic inference from genomic data, coupled with simulations, allows assessment of extinction risk and different models of recovery.