The genomic history of the Iberian Peninsula over the past 8000 years

Genetic adaptation to pathogens and increased risk of inflammatory disorders in post-Neolithic Europe

Ancient genomics can directly detect human genetic adaptation to environmental cues. However, it remains unclear how pathogens have exerted selective pressures on human genome diversity across different epochs and affected present-day inflammatory disease risk. Here, we use an ancestry-aware approximate Bayesian computation framework to estimate the nature, strength, and time of onset of selection acting on 2,879 ancient and modern European genomes from the last 10,000 years. We found that the bulk of genetic adaptation occurred after the start of the Bronze Age, <4,500 years ago, and was enriched in genes relating to host-pathogen interactions. Furthermore, we detected directional selection acting on specific leukocytic lineages and experimentally demonstrated that the strongest negatively selected candidate variant in immunity genes, lipopolysaccharide-binding protein (LBP) D283G, is hypomorphic. Finally, our analyses suggest that the risk of inflammatory disorders has increased in post-Neolithic Europeans, possibly because of antagonistic pleiotropy following genetic adaptation to pathogens.

The power of geohistorical boundaries for modeling the genetic background of human populations: The case of the rural catalan Pyrenees

The genetic variation of the European population at a macro-geographic scale follows genetic gradients which reflect main migration events. However, less is known about factors affecting mating patterns at a micro-geographic scale. In this study we have analyzed 726,718 autosomal single nucleotide variants in 435 individuals from the catalan Pyrenees covering around 200 km of a vast and abrupt region in the north of the Iberian Peninsula, for which we have information about the geographic origin of all grand-parents and parents. At a macro-geographic scale, our analyses recapitulate the genetic gradient observed in Spain. However, we also identified the presence of micro-population substructure among the sampled individuals. Such micro-population substructure does not correlate with geographic barriers such as the expected by the orography of the considered region, but by the bishoprics present in the covered geographic area. These results support that, on top of main human migrations, long ongoing socio-cultural factors have also shaped the genetic diversity observed at rural populations.

Characterization of Danube Swabian population samples on a high-resolution genome-wide basis

BACKGROUND

German-derived ethnicities are one of the largest ethnic groups in Hungary, dating back to the formation of the Kingdom of Hungary, which took place at the beginning of the 11th century. Germans arrived in Hungary in many waves. The most significant immigration wave took place following the collapse of the Ottoman Empire in East-Central Europe which closed the 150 year long Ottoman occupation. To date, there are no comprehensive genome-wide studies investigating the genetic makeup of the Danube Swabians. Here we analyzed 47 Danube Swabian samples collected from elderly Swabian individuals living in the Dunaszekcső-Bár area, in Danube side villages of Southwest Hungary. These Swabians, according to self-declaration, did not admix with other ethnic groups for 3-6 succeeding generations. Using Illumina Infinium 720 K Beadchip genotype data, we applied allele frequency-based and haplotype-based genome-wide marker data analyses to investigate the ancestry and genetic composition of the collected Danube Swabian samples.

RESULTS

Haplotype-based analyses like identity by descent segment analysis show that the investigated Danube Swabians possess significant German and other West European ancestry, but their Hungarian ancestry is also prominent. Our results suggest that their main source of ancestry can be traced back to Western Europe, presumably to the region of Germany.

CONCLUSION

This is the first analysis of Danube Swabian population samples based on genome-wide autosomal data. Our results establish the basis for conducting further comprehensive research on Danube Swabians and on other German ethnicities of the Carpathian basin, which can help reconstruct their origin, and identify their major archaic genomic patterns.

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

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

The genetic history of Scandinavia from the Roman Iron Age to the present

We investigate a 2,000-year genetic transect through Scandinavia spanning the Iron Age to the present, based on 48 new and 249 published ancient genomes and genotypes from 16,638 modern individuals. We find regional variation in the timing and magnitude of gene flow from three sources: the eastern Baltic, the British-Irish Isles, and southern Europe. British-Irish ancestry was widespread in Scandinavia from the Viking period, whereas eastern Baltic ancestry is more localized to Gotland and central Sweden. In some regions, a drop in current levels of external ancestry suggests that ancient immigrants contributed proportionately less to the modern Scandinavian gene pool than indicated by the ancestry of genomes from the Viking and Medieval periods. Finally, we show that a north-south genetic cline that characterizes modern Scandinavians is mainly due to the differential levels of Uralic ancestry and that this cline existed in the Viking Age and possibly earlier.

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.

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.

Kinship practices in the early state El Argar society from Bronze Age Iberia

The Early Bronze Age in Europe is characterized by social and genetic transformations, starting in the early 3rd millennium BCE. New settlement and funerary structures, artifacts and techniques indicate times of change with increasing economic asymmetries and political hierarchization. Technological advances in metallurgy also played an important role, facilitating trade and exchange networks, which became tangible in higher levels of mobility and connectedness. Archeogenetic studies have revealed a substantial transformation of the genetic ancestry around this time, ultimately linked to the expansion of steppe- and forest steppe pastoralists from Eastern Europe. Evidence for emerging infectious diseases such as Yersinia pestis adds further complexity to these tumultuous and transformative times. The El Argar complex in southern Iberia marks the genetic turnover in southwestern Europe ~ 2200 BCE that accompanies profound changes in the socio-economic structure of the region. To answer the question of who was buried in the emblematic double burials of the El Argar site La Almoloya, we integrated results from biological relatedness analyses and archaeological funerary contexts and refined radiocarbon-based chronologies from 68 individuals. We find that the El Argar society was virilocally and patrilineally organized and practiced reciprocal female exogamy, supported by pedigrees that extend up to five generations along the paternal line. Synchronously dated adult males and females from double tombs were found to be unrelated mating partners, whereby the incoming females reflect socio-political alliances among El Argar groups. In three cases these unions had common offspring, while paternal half-siblings also indicate serial monogamy or polygyny.

de Pancorbo M, Calafell F. Y-chromosome target enrichment reveals rapid expansion of haplogroup R1b-DF27 in Iberia during the Bronze Age transition

The Y chromosome can yield a unique perspective into the study of human demographic history. However, due to the repetitive nature of part of its sequence, only a small set of regions are suitable for variant calling and discovery from short-read sequencing data. These regions combined represent 8.9 Mbp or 0.14% of a diploid human genome. Consequently, investing in whole-genome sequencing to resolve Y-chromosome questions is poorly efficient. Here we use, as an alternative, target enrichment technology to greatly increase sequencing effectiveness, validating and applying the technique to 181 males, for 162 of whom we obtained a positive result. Additionally, 75 samples sequenced for the whole genome were also included, for a total sample size of 237. These samples were chosen for their Y chromosome haplogroup: R1b-DF27. In the context of European populations, and particularly in Iberia, this haplogroup stands out for its high frequency and its demographic history. Current evidence indicates that the diffusion of this haplogroup is related to the population movements that mark the cultural Bronze Age transition, making it remarkably interesting for population geneticists. The results of this study show the effects of the rapid radiation of the haplogroup in Spain, as even with the higher discriminating power of whole sequences, most haplotypes still fall within the R1b-DF27* paragroup rather than in the main derived branches. However, we were able to refine the ISOGG 2019-2020 phylogeny, and its two main subbranches, namely L176.2 and Z272, which present geographical differentiation between the Atlantic and Mediterranean coasts of Iberia.

Estimation of site frequency spectra from low-coverage sequencing data using stochastic EM reduces overfitting, runtime, and memory usage

The site frequency spectrum is an important summary statistic in population genetics used for inference on demographic history and selection. However, estimation of the site frequency spectrum from called genotypes introduces bias when working with low-coverage sequencing data. Methods exist for addressing this issue but sometimes suffer from 2 problems. First, they can have very high computational demands, to the point that it may not be possible to run estimation for genome-scale data. Second, existing methods are prone to overfitting, especially for multidimensional site frequency spectrum estimation. In this article, we present a stochastic expectation-maximization algorithm for inferring the site frequency spectrum from NGS data that address these challenges. We show that this algorithm greatly reduces runtime and enables estimation with constant, trivial RAM usage. Furthermore, the algorithm reduces overfitting and thereby improves downstream inference. An implementation is available at github.com/malthesr/winsfs.

Genomic ancestry, diet and microbiomes of Upper Palaeolithic hunter-gatherers from San Teodoro cave

Recent improvements in the analysis of ancient biomolecules from human remains and associated dental calculus have provided new insights into the prehistoric diet and genetic diversity of our species. Here we present a multi-omics study, integrating metagenomic and proteomic analyses of dental calculus, and human ancient DNA analysis of the petrous bones of two post-Last Glacial Maximum (LGM) individuals from San Teodoro cave (Italy), to reconstruct their lifestyle and the post-LGM resettlement of Europe. Our analyses show genetic homogeneity in Sicily during the Palaeolithic, representing a hitherto unknown Italian genetic lineage within the previously identified Villabruna cluster. We argue that this lineage took refuge in Italy during the LGM, followed by a subsequent spread to central-western Europe. Analysis of dental calculus showed a diet rich in animal proteins which is also reflected on the oral microbiome composition. Our results demonstrate the power of this approach in the study of prehistoric humans and will enable future research to reach a more holistic understanding of the population dynamics and ecology.

The arrival of millets to the Atlantic coast of northern Iberia

Despite being one of the most important crops in the recent prehistory of Eurasia, the arrival and exploitation of millets in the westernmost part of Europe are still largely underexplored. Here and for the first time, we report multipronged biomolecular evidence of millet consumption along the Atlantic façade of northern Iberia through a combination of radiocarbon dating, stable isotopes, and dental calculus analyses on the human individuals found in the burial site of El Espinoso cave (Asturias, Spain). The high-resolution chronological framework established for individuals placed the burials between 1235 and 1099 cal. BC. The discovery of high δ¹³C values on their bone collagen and the identification of polyhedral starch grains within their dental plaque underline the relevance of C₄ plants in their diet and highlights the timing of the systematic consumption of millets in the Late Bronze Age. Our data support previous regional archaeobotanical evidence and establish a more precise chronology of the dispersal of millets into northern Iberia during the Bronze Age, becoming an essential crop until the arrival of maize from America after AD 1492. This study emphasizes the importance of multidisciplinary methods to ascertain the origin and development of agricultural practices during recent prehistory.

Neolithic genomic data from southern France showcase intensified interactions with hunter-gatherer communities

Archaeological research shows that the dispersal of the Neolithic took a more complex turn when reaching western Europe, painting a contrasted picture of interactions between autochthonous hunter-gatherers (HGs) and incoming farmers. In order to clarify the mode, the intensity, and the regional variability of biological exchanges implied in these processes, we report new palaeogenomic data from Occitanie, a key region in Southern France. Genomic data from 28 individuals originating from six sites spanning from c. 5,500 to c. 2,500 BCE allow us to characterize regional patterns of ancestries throughout the Neolithic period. Results highlight major differences between the Mediterranean and Continental Neolithic expansion routes regarding both migration and interaction processes. High proportions of HG ancestry in both Early and Late Neolithic groups in Southern France support multiple pulses of inter-group gene flow throughout time and space and confirm the need for regional studies to address the complexity of the processes involved.

•

Genome-wide data from 28 individuals from Southern France (∼5,500–∼2,500 BCE)

•

Small groups associated with the Neolithic expansion along the Mediterranean

•

Early admixture between hunter-gatherers and Neolithic farmers in Southern France

•

Multiple pulses of HG legacy introgression in Western Europe throughout Neolithic

The diverse genetic origins of a Classical period Greek army

Trade and colonization caused an unprecedented increase in Mediterranean human mobility in the first millennium BCE. Often seen as a dividing force, warfare is in fact another catalyst of culture contact. We provide insight into the demographic dynamics of ancient warfare by reporting genome-wide data from fifth-century soldiers who fought for the army of the Greek Sicilian colony of Himera, along with representatives of the civilian population, nearby indigenous settlements, and 96 present-day individuals from Italy and Greece. Unlike the rest of the sample, many soldiers had ancestral origins in northern Europe, the Steppe, and the Caucasus. Integrating genetic, archaeological, isotopic, and historical data, these results illustrate the significant role mercenaries played in ancient Greek armies and highlight how participation in war contributed to continental-scale human mobility in the Classical world.

Genetic structure and differentiation from early bronze age in the mediterranean island of sicily: Insights from ancient mitochondrial genomes

Sicily is one of the main islands of the Mediterranean Sea, and it is characterized by a variety of archaeological records, material culture and traditions, reflecting the history of migrations and populations’ interaction since its first colonization, during the Paleolithic. These deep and complex demographic and cultural dynamics should have affected the genomic landscape of Sicily at different levels; however, the relative impact of these migrations on the genomic structure and differentiation within the island remains largely unknown. The available Sicilian modern genetic data gave a picture of the current genetic structure, but the paucity of ancient data did not allow so far to make predictions about the level of historical variation. In this work, we sequenced and analyzed the complete mitochondrial genomes of 36 individuals from five different locations in Sicily, spanning from Early Bronze Age to Iron Age, and with different cultural backgrounds. The comparison with coeval groups from the Mediterranean Basin highlighted structured genetic variation in Sicily since Early Bronze Age, thus supporting a demic impact of the cultural transitions within the Island. Explicit model testing through Approximate Bayesian Computation allowed us to make predictions about the origin of Sicanians, one of the three indigenous peoples of Sicily, whose foreign origin from Spain, historically attributed, was not confirmed by our analysis of genetic data. Sicilian modern mitochondrial data show a different, more homogeneous, genetic composition, calling for a recent genetic replacement in the Island of pre-Iron Age populations, that should be further investigated.

Performance of the Use of Genetic Information to Assess the Risk of Colorectal Cancer in the Basque Population

Although the genetic contribution to colorectal cancer (CRC) has been studied in various populations, studies on the applicability of available genetic information in the Basque population are scarce. In total, 835 CRC cases and 940 controls from the Basque population were genotyped and genome-wide association studies were carried out. Mendelian Randomization analyses were used to discover the effect of modifiable risk factors and microbiota on CRC. In total, 25 polygenic risk score models were evaluated to assess their performance in CRC risk calculation. Moreover, 492 inflammatory bowel disease cases were used to assess whether that genetic information would not confuse both conditions. Five suggestive (p < 5 × 10−6) loci were associated with CRC risk, where genes previously associated with CRC were located (e.g., ABCA12, ATIC or ERBB4). Moreover, the analyses of CRC locations detected additional genes consistent with the biology of CRC. The possible contribution of cholesterol, BMI, Firmicutes and Cyanobacteria to CRC risk was detected by Mendelian Randomization. Finally, although polygenic risk score models showed variable performance, the best model performed correctly regardless of the location and did not misclassify inflammatory bowel disease cases. Our results are consistent with CRC biology and genetic risk models and could be applied to assess CRC risk in the Basque population.

Principal Component Analyses (PCA)-based findings in population genetic studies are highly biased and must be reevaluated

Principal Component Analysis (PCA) is a multivariate analysis that reduces the complexity of datasets while preserving data covariance. The outcome can be visualized on colorful scatterplots, ideally with only a minimal loss of information. PCA applications, implemented in well-cited packages like EIGENSOFT and PLINK, are extensively used as the foremost analyses in population genetics and related fields (e.g., animal and plant or medical genetics). PCA outcomes are used to shape study design, identify, and characterize individuals and populations, and draw historical and ethnobiological conclusions on origins, evolution, dispersion, and relatedness. The replicability crisis in science has prompted us to evaluate whether PCA results are reliable, robust, and replicable. We analyzed twelve common test cases using an intuitive color-based model alongside human population data. We demonstrate that PCA results can be artifacts of the data and can be easily manipulated to generate desired outcomes. PCA adjustment also yielded unfavorable outcomes in association studies. PCA results may not be reliable, robust, or replicable as the field assumes. Our findings raise concerns about the validity of results reported in the population genetics literature and related fields that place a disproportionate reliance upon PCA outcomes and the insights derived from them. We conclude that PCA may have a biasing role in genetic investigations and that 32,000-216,000 genetic studies should be reevaluated. An alternative mixed-admixture population genetic model is discussed.

1,000 ancient genomes uncover 10,000 years of natural selection in Europe

Ancient DNA has revolutionized our understanding of human population history. However, its potential to examine how rapid cultural evolution to new lifestyles may have driven biological adaptation has not been met, largely due to limited sample sizes. We assembled genome-wide data from 1,291 individuals from Europe over 10,000 years, providing a dataset that is large enough to resolve the timing of selection into the Neolithic, Bronze Age, and Historical periods. We identified 25 genetic loci with rapid changes in frequency during these periods, a majority of which were previously undetected. Signals specific to the Neolithic transition are associated with body weight, diet, and lipid metabolism-related phenotypes. They also include immune phenotypes, most notably a locus that confers immunity to Salmonella infection at a time when ancient Salmonella genomes have been shown to adapt to human hosts, thus providing a possible example of human-pathogen co-evolution. In the Bronze Age, selection signals are enriched near genes involved in pigmentation and immune-related traits, including at a key human protein interactor of SARS-CoV-2. Only in the Historical period do the selection candidates we detect largely mirror previously-reported signals, highlighting how the statistical power of previous studies was limited to the last few millennia. The Historical period also has multiple signals associated with vitamin D binding, providing evidence that lactase persistence may have been part of an oligogenic adaptation for efficient calcium uptake and challenging the theory that its adaptive value lies only in facilitating caloric supplementation during times of scarcity. Finally, we detect selection on complex traits in all three periods, including selection favoring variants that reduce body weight in the Neolithic. In the Historical period, we detect selection favoring variants that increase risk for cardiovascular disease plausibly reflecting selection for a more active inflammatory response that would have been adaptive in the face of increased infectious disease exposure. Our results provide an evolutionary rationale for the high prevalence of these deadly diseases in modern societies today and highlight the unique power of ancient DNA in elucidating biological change that accompanied the profound cultural transformations of recent human history.

The genetic history of the Southern Arc: A bridge between West Asia and Europe

By sequencing 727 ancient individuals from the Southern Arc (Anatolia and its neighbors in Southeastern Europe and West Asia) over 10,000 years, we contextualize its Chalcolithic period and Bronze Age (about 5000 to 1000 BCE), when extensive gene flow entangled it with the Eurasian steppe. Two streams of migration transmitted Caucasus and Anatolian/Levantine ancestry northward, and the Yamnaya pastoralists, formed on the steppe, then spread southward into the Balkans and across the Caucasus into Armenia, where they left numerous patrilineal descendants. Anatolia was transformed by intra-West Asian gene flow, with negligible impact of the later Yamnaya migrations. This contrasts with all other regions where Indo-European languages were spoken, suggesting that the homeland of the Indo-Anatolian language family was in West Asia, with only secondary dispersals of non-Anatolian Indo-Europeans from the steppe.

A genetic probe into the ancient and medieval history of Southern Europe and West Asia

Literary and archaeological sources have preserved a rich history of Southern Europe and West Asia since the Bronze Age that can be complemented by genetics. Mycenaean period elites in Greece did not differ from the general population and included both people with some steppe ancestry and others, like the Griffin Warrior, without it. Similarly, people in the central area of the Urartian Kingdom around Lake Van lacked the steppe ancestry characteristic of the kingdom's northern provinces. Anatolia exhibited extraordinary continuity down to the Roman and Byzantine periods, with its people serving as the demographic core of much of the Roman Empire, including the city of Rome itself. During medieval times, migrations associated with Slavic and Turkic speakers profoundly affected the region.

Dairying, diseases and the evolution of lactase persistence in Europe

In European and many African, Middle Eastern and southern Asian populations, lactase persistence (LP) is the most strongly selected monogenic trait to have evolved over the past 10,000 years1. Although the selection of LP and the consumption of prehistoric milk must be linked, considerable uncertainty remains concerning their spatiotemporal configuration and specific interactions2,3. Here we provide detailed distributions of milk exploitation across Europe over the past 9,000 years using around 7,000 pottery fat residues from more than 550 archaeological sites. European milk use was widespread from the Neolithic period onwards but varied spatially and temporally in intensity. Notably, LP selection varying with levels of prehistoric milk exploitation is no better at explaining LP allele frequency trajectories than uniform selection since the Neolithic period. In the UK Biobank4,5 cohort of 500,000 contemporary Europeans, LP genotype was only weakly associated with milk consumption and did not show consistent associations with improved fitness or health indicators. This suggests that other reasons for the beneficial effects of LP should be considered for its rapid frequency increase. We propose that lactase non-persistent individuals consumed milk when it became available but, under conditions of famine and/or increased pathogen exposure, this was disadvantageous, driving LP selection in prehistoric Europe. Comparison of model likelihoods indicates that population fluctuations, settlement density and wild animal exploitation-proxies for these drivers-provide better explanations of LP selection than the extent of milk exploitation. These findings offer new perspectives on prehistoric milk exploitation and LP evolution.

Ancient DNA reveals five streams of migration into Micronesia and matrilocality in early Pacific seafarers

Micronesia began to be peopled earlier than other parts of Remote Oceania, but the origins of its inhabitants remain unclear. We generated genome-wide data from 164 ancient and 112 modern individuals. Analysis reveals five migratory streams into Micronesia. Three are East Asian related, one is Polynesian, and a fifth is a Papuan source related to mainland New Guineans that is different from the New Britain-related Papuan source for southwest Pacific populations but is similarly derived from male migrants ~2500 to 2000 years ago. People of the Mariana Archipelago may derive all of their precolonial ancestry from East Asian sources, making them the only Remote Oceanians without Papuan ancestry. Female-inherited mitochondrial DNA was highly differentiated across early Remote Oceanian communities but homogeneous within, implying matrilocal practices whereby women almost never raised their children in communities different from the ones in which they grew up.

Ancient Maltese genomes and the genetic geography of Neolithic Europe

Archaeological consideration of maritime connectivity has ranged from a biogeographical perspective that considers the sea as a barrier to a view of seaways as ancient highways that facilitate exchange. Our results illustrate the former. We report three Late Neolithic human genomes from the Mediterranean island of Malta that are markedly enriched for runs of homozygosity, indicating inbreeding in their ancestry and an effective population size of only hundreds, a striking illustration of maritime isolation in this agricultural society. In the Late Neolithic, communities across mainland Europe experienced a resurgence of hunter-gatherer ancestry, pointing toward the persistence of different ancestral strands that subsequently admixed. This is absent in the Maltese genomes, giving a further indication of their genomic insularity. Imputation of genome-wide genotypes in our new and 258 published ancient individuals allowed shared identity-by-descent segment analysis, giving a fine-grained genetic geography of Neolithic Europe. This highlights the differentiating effects of seafaring Mediterranean expansion and also island colonization, including that of Ireland, Britain, and Orkney. These maritime effects contrast profoundly with a lack of migratory barriers in the establishment of Central European farming populations from Anatolia and the Balkans.

hapCon: Estimating Contamination of Ancient Genomes by Copying from Reference Haplotypes

Motivation

Human ancient DNA (aDNA) studies have surged in recent years, revolutionizing the study of the human past. Typically, aDNA is preserved poorly, making such data prone to contamination from other human DNA. Therefore, it is important to rule out substantial contamination before proceeding to downstream analysis. As most aDNA samples can only be sequenced to low coverages (<1× average depth), computational methods that can robustly estimate contamination in the low coverage regime are needed. However, the ultra low-coverage regime (0.1× and below) remains a challenging task for existing approaches.

Results

We present a new method to estimate contamination in aDNA for male modern humans. It utilizes a Li&Stephens haplotype copying model for haploid X chromosomes, with mismatches modeled as errors or contamination. We assessed this new approach, hapCon, on simulated and down-sampled empirical aDNA data. Our experiments demonstrate that hapCon outperforms a commonly used tool for estimating male X contamination (ANGSD), with substantially lower variance and narrower confidence intervals, especially in the low coverage regime. We found that hapCon provides useful contamination estimates for coverages as low as 0.1× for SNP capture data (1240k) and 0.02× for whole genome sequencing data, substantially extending the coverage limit of previous male X chromosome-based contamination estimation methods. Our experiments demonstrate that hapCon has little bias for contamination up to 25–30% as long as the contaminating source is specified within continental genetic variation, and that its application range extends to human aDNA as old as ∼45 000 and various global ancestries.

Availability and implementation

We make hapCon available as part of a python package (hapROH), which is available at the Python Package Index (https://pypi.org/project/hapROH) and can be installed via pip. The documentation provides example use cases as blueprints for custom applications (https://haproh.readthedocs.io/en/latest/hapCon.html). The program can analyze either BAM files or pileup files produced with samtools. An implementation of our software (hapCon) using Python and C is deposited at https://github.com/hyl317/hapROH.

Supplementary information

Supplementary data are available at Bioinformatics online.

At the beginnings of the funerary Megalithism in Iberia at Campo de Hockey necropolis

The excavations undertaken at the Campo de Hockey site in 2008 led to the identification of a major Neolithic necropolis in the former Island of San Fernando (Bay of Cádiz). This work presents the results of the latest studies, which indicate that the site stands as one of the oldest megalithic necropolises in the Iberian Peninsula. The main aim of this work is to present with precision the chronology of this necropolis through a Bayesian statistical model that confirms that the necropolis was in use from c. 4300 to 3800 cal BC. The presence of prestige grave goods in the earliest and most monumental graves suggest that the Megalithism phenomenon emerged in relation to maritime routes linked to the distribution of exotic products. We also aim to examine funerary practices in these early megalithic communities, and especially their way of life and the social reproduction system. As such, in addition to the chronological information and the Bayesian statistics, we provide the results of a comprehensive interdisciplinary study, including anthropological, archaeometric and genetic data.

The spatiotemporal patterns of major human admixture events during the European Holocene

Recent studies have shown that admixture has been pervasive throughout human history. While several methods exist for dating admixture in contemporary populations, they are not suitable for sparse, low coverage ancient genomic data. Thus, we developed DATES (Distribution of Ancestry Tracts of Evolutionary Signals) that leverages ancestry covariance patterns across the genome of a single individual to infer the timing of admixture. DATES provides reliable estimates under various demographic scenarios and outperforms available methods for ancient DNA applications. Using DATES on~1100 ancient genomes from sixteen regions in Europe and west Asia, we reconstruct the chronology of the formation of the ancestral populations and the fine-scale details of the spread of Neolithic farming and Steppe pastoralist-related ancestry across Europe. By studying the genetic formation of Anatolian farmers, we infer that gene flow related to Iranian Neolithic farmers occurred before 9600 BCE, predating the advent of agriculture in Anatolia. Contrary to the archaeological evidence, we estimate that early Steppe pastoralist groups (Yamnaya and Afanasievo) were genetically formed more than a millennium before the start of Steppe pastoralism. Our analyses provide new insights on the origins and spread of farming and Indo-European languages, highlighting the power of genomic dating methods to elucidate the legacy of human migrations.

Bioarchaeological and palaeogenomic portrait of two Pompeians that died during the eruption of Vesuvius in 79 AD

The archaeological site of Pompeii is one of the 54 UNESCO World Heritage sites in Italy, thanks to its uniqueness: the town was completely destroyed and buried by a Vesuvius’ eruption in 79 AD. In this work, we present a multidisciplinary approach with bioarchaeological and palaeogenomic analyses of two Pompeian human remains from the Casa del Fabbro. We have been able to characterize the genetic profile of the first Pompeian’ genome, which has strong affinities with the surrounding central Italian population from the Roman Imperial Age. Our findings suggest that, despite the extensive connection between Rome and other Mediterranean populations, a noticeable degree of genetic homogeneity exists in the Italian peninsula at that time. Moreover, palaeopathological analyses identified the presence of spinal tuberculosis and we further investigated the presence of ancient DNA from Mycobacterium tuberculosis. In conclusion, our study demonstrates the power of a combined approach to investigate ancient humans and confirms the possibility to retrieve ancient DNA from Pompeii human remains. Our initial findings provide a foundation to promote an intensive and extensive paleogenetic analysis in order to reconstruct the genetic history of population from Pompeii, a unique archaeological site.

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.

Salorno—Dos de la Forca (Adige Valley, Northern Italy): A unique cremation site of the Late Bronze Age

The archaeological site of Salorno—Dos de la Forca (Bozen, Alto Adige) provides one of the rarest and most significant documentations of cremated human remains preserved from an ancient cremation platform (ustrinum). The pyre area, located along the upper Adige valley, is dated to the Late Bronze Age (ca. 1150–950 BCE) and has yielded an unprecedented quantity of cremated human remains (about 63.5 kg), along with burnt animal bone fragments, shards of pottery, and other grave goods made in bronze and animal bone/antler. This study focuses on the bioanthropological analysis of the human remains and discusses the formation of the unusual burnt deposits at Salorno through comparisons with modern practices and protohistoric and contemporaneous archaeological deposits. The patterning of bone fragmentation and commingling was investigated using spatial data recorded during excavation which, along with the bioanthropological and archaeological data, are used to model and test two hypotheses: Salorno—Dos de la Forca would be the result of A) repeated primary cremations left in situ; or B) of residual material remaining after select elements were removed for internment in urns or burials to unknown depositional sites. By modelling bone weight and demographic data borrowed from regional affine contexts, the authors suggest that this cremation site may have been used over several generations by a small community–perhaps a local elite. With a quantity of human remains that exceeds that of any other coeval contexts interpreted as ustrina, Salorno may be the product of a complex series of rituals in which the human cremains did not receive individual burial, but were left in situ, in a collective/communal place of primary combustion, defining an area of repeated funeral ceremonies involving offerings and libations across a few generations. This would represent a new typological and functional category that adds to the variability of mortuary customs at the end of the Bronze Age in the Alpine are, at a time in which “globalising” social trends may have stimulated the definition of more private identities.

Genomes from Verteba cave suggest diversity within the Trypillians in Ukraine

The transition to agriculture occurred relatively late in Eastern Europe, leading researchers to debate whether it was a gradual, interactive process or a colonisation event. In the forest and forest-steppe regions of Ukraine, farming appeared during the fifth millennium BCE, associated with the Cucuteni-Trypillia cultural complex (CTCC, ~ 5000–3000 BCE). Across Europe, the Neolithisation process was highly variable across space and over time. Here, we investigate the population dynamics of early agriculturalists from the eastern forest-steppe region based on the analyses of 20 ancient genomes from the site of Verteba Cave (3935–825 cal BCE). Results reveal that the CTCC individuals’ ancestry is related to both western hunter-gatherers and Near Eastern farmers, has no local ancestry associated with Ukrainian Neolithic hunter-gatherers and has steppe ancestry. An Early Bronze Age individual has an ancestry profile related to the Yamnaya expansions but with 20% of ancestry related to the other Trypillian individuals, which suggests admixture between the Trypillians and the incoming populations carrying steppe-related ancestry. A Late Bronze Age individual dated to 980–825 cal BCE has a genetic profile indicating affinity to Beaker-related populations, detected close to 1000 years after the end of the Bell Beaker phenomenon during the third millennium BCE.

Genomic and dietary discontinuities during the Mesolithic and Neolithic in Sicily

Sicily is a key region for understanding the agricultural transition in the Mediterranean because of its central position. Here, we present genomic and stable isotopic data for 19 prehistoric Sicilians covering the Mesolithic to Bronze Age periods (10,700–4,100 yBP). We find that Early Mesolithic hunter-gatherers (HGs) from Sicily are a highly drifted lineage of the Early Holocene western European HGs, whereas Late Mesolithic HGs carry ∼20% ancestry related to northern and (south) eastern European HGs, indicating substantial gene flow. Early Neolithic farmers are genetically most similar to farmers from the Balkans and Greece, with only ∼7% of ancestry from local Mesolithic HGs. The genetic discontinuities during the Mesolithic and Early Neolithic match the changes in material culture and diet. Three outlying individuals dated to ∼8,000 yBP; however, suggest that hunter-gatherers interacted with incoming farmers at Grotta dell’Uzzo, resulting in a mixed economy and diet for a brief interlude at the Mesolithic-Neolithic transition.

•

Genetic transition between Early Mesolithic and Late Mesolithic hunter-gatherers

•

A near-complete genetic turnover during the Mesolithic-Neolithic transition

•

Exchange of subsistence practices between hunter-gatherers and early farmers

Incidence Trends of Inflammatory Bowel Disease in a Southern European Country: A Mirror of the Western World?

INTRODUCTION

Inflammatory bowel disease (IBD) affects people from all age categories worldwide. Although the incidence of the disease is stabilizing or decreasing in most Western world countries, its prevalence is still increasing because of the rise in life expectancy and better disease management. This work intends to identify the trends related to IBD incidence nationwide, analyzing regional, sex, and age distributions.

METHODS

Data were provided by the Portuguese Shared Services of the Ministry of Health. This study consisted of a retrospective analysis of all first consultations coded for "Chronic enteritis/ulcerative colitis" (D94) in a primary healthcare setting, between 2017 and 2020, in Portugal. The primary outcome measure was the IBD incidence rate per 100,000 inhabitants. We also calculated the incidence rate per person-year and forecasted incidence until 2024.

RESULTS

Between 2017 and 2019, the incidence rate of IBD in Portugal decreased from 54.9 to 48.6 per 100,000 inhabitants. The average incidence was 20 new cases of IBD per 1,000 person-year. It was predicted that, in December 2023, IBD incidence would reach 305.4 new cases (95% Prediction Interval 156.6-454.3), a similar result to the values forecasted for December 2021 (305.4, 95% Prediction Interval 197.3-413.6).

DISCUSSION

The incidence of IBD slightly declined from 2017 to 2019, and it is posed to stabilize in the future. The presented data are of the utmost importance for the characterization of IBD in Southern European countries and the establishment of future health policies in the setting of compounding prevalence in the Western world.

Stone Age Yersinia pestis genomes shed light on the early evolution, diversity, and ecology of plague

The bacterium Yersinia pestis has caused numerous historically documented outbreaks of plague and research using ancient DNA could demonstrate that it already affected human populations during the Neolithic. However, the pathogen’s genetic diversity, geographic spread, and transmission dynamics during this early period of Y. pestis evolution are largely unexplored. Here, we describe a set of ancient plague genomes up to 5,000 y old from across Eurasia. Our data demonstrate that two genetically distinct forms of Y. pestis evolved in parallel and were both distributed across vast geographic distances, potentially occupying different ecological niches. Interpreted within the archeological context, our results suggest that the spread of plague during this period was linked to increased human mobility and intensification of animal husbandry.

The bacterial pathogen Yersinia pestis gave rise to devastating outbreaks throughout human history, and ancient DNA evidence has shown it afflicted human populations as far back as the Neolithic. Y. pestis genomes recovered from the Eurasian Late Neolithic/Early Bronze Age (LNBA) period have uncovered key evolutionary steps that led to its emergence from a Yersinia pseudotuberculosis-like progenitor; however, the number of reconstructed LNBA genomes are too few to explore its diversity during this critical period of development. Here, we present 17 Y. pestis genomes dating to 5,000 to 2,500 y BP from a wide geographic expanse across Eurasia. This increased dataset enabled us to explore correlations between temporal, geographical, and genetic distance. Our results suggest a nonflea-adapted and potentially extinct single lineage that persisted over millennia without significant parallel diversification, accompanied by rapid dispersal across continents throughout this period, a trend not observed in other pathogens for which ancient genomes are available. A stepwise pattern of gene loss provides further clues on its early evolution and potential adaptation. We also discover the presence of the flea-adapted form of Y. pestis in Bronze Age Iberia, previously only identified in in the Caucasus and the Volga regions, suggesting a much wider geographic spread of this form of Y. pestis. Together, these data reveal the dynamic nature of plague’s formative years in terms of its early evolution and ecology.

Origin and mobility of Iron Age Gaulish groups in present-day France revealed through archaeogenomics

The Iron Age period occupies an important place in French history because the Gauls are regularly presented as the direct ancestors of the extant French population. We documented here the genomic diversity of Iron Age communities originating from six French regions. The 49 acquired genomes permitted us to highlight an absence of discontinuity between Bronze Age and Iron Age groups in France, lending support to a cultural transition linked to progressive local economic changes rather than to a massive influx of allochthone groups. Genomic analyses revealed strong genetic homogeneity among the regional groups associated with distinct archaeological cultures. This genomic homogenization appears to be linked to individuals' mobility between regions and gene flow with neighbouring groups from England and Spain. Thus, the results globally support a common genomic legacy for the Iron Age population of modern-day France that could be linked to recurrent gene flow between culturally differentiated communities.

From Human Remains to Powerful Objects: Ancestor Research from a Deep-Time Perspective

Family history research has seen a surge in popularity in recent years; however, is this preoccupation with who we are and where we come from new? Archaeological evidence suggests that ancestors played crucial and ubiquitous roles in the identities and cosmologies of past societies. This paper will explore how, in the absence of genealogical websites and DNA testing, kinship structures and understandings of personhood beyond genealogy may have influenced concepts of ancestry. Case studies from later prehistoric Britain will demonstrate the ways in which monuments, objects and human remains themselves created bonds between the living and the dead, prompting us to reflect on genealogy as just one aspect of our identity in the present.

Local genetic variation of inflammatory bowel disease in Basque population and its effect in risk prediction

Inflammatory bowel disease (IBD) is characterised by chronic inflammation of the gastrointestinal tract. Although its aetiology remains unknown, environmental and genetic factors are involved in its development. Regarding genetics, more than 200 loci have been associated with IBD but the transferability of those signals to the Basque population living in Northern Spain, a population with distinctive genetic background, remains unknown. We have analysed 5,411,568 SNPs in 498 IBD cases and 935 controls from the Basque population. We found 33 suggestive loci (p < 5 × 10⁻⁶) in IBD and its subtypes, namely Crohn’s Disease (CD) and Ulcerative Colitis (UC), detecting a genome-wide significant locus located in HLA region in patients with UC. Those loci contain previously associated genes with IBD (IL23R, JAK2 or HLA genes) and new genes that could be involved in its development (AGT, BZW2 or FSTL1). The overall genetic correlation between European populations and Basque population was high in IBD and CD, while in UC was lower. Finally, the use of genetic risk scores based on previous GWAS findings reached area under the curves > 0.68. In conclusion, we report on the genetic architecture of IBD in the Basque population, and explore the performance of European-descent genetic risk scores in this population.

A unified genealogy of modern and ancient genomes

The sequencing of modern and ancient genomes from around the world has revolutionized our understanding of human history and evolution. However, the problem of how best to characterize ancestral relationships from the totality of human genomic variation remains unsolved. Here, we address this challenge with nonparametric methods that enable us to infer a unified genealogy of modern and ancient humans. This compact representation of multiple datasets explores the challenges of missing and erroneous data and uses ancient samples to constrain and date relationships. We demonstrate the power of the method to recover relationships between individuals and populations as well as to identify descendants of ancient samples. Finally, we introduce a simple nonparametric estimator of the geographical location of ancestors that recapitulates key events in human history.

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

The Orcadian Neolithic has been intensively studied and celebrated as a major center of cultural innovation, whereas the Bronze Age is less well known and often regarded as a time of stagnation and insularity. Here, we analyze ancient genomes from the Orcadian Bronze Age in the context of the variation in Neolithic Orkney and Bronze Age Europe. We find clear evidence for Early Bronze Age immigration into Orkney, but with an extraordinary pattern: continuity from the Neolithic on the male line of descent but immigration from continental Europe on the female side, echoed in the genome-wide picture. This suggests that despite substantial immigration, indigenous male lineages persisted for at least a thousand years after the end of the Neolithic.

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

The genetic origin of Daunians and the Pan-Mediterranean southern Italian Iron Age context

The geographical location and shape of Apulia, a narrow land stretching out in the sea at the South of Italy, made this region a Mediterranean crossroads connecting Western Europe and the Balkans. Such movements culminated at the beginning of the Iron Age with the Iapygian civilization which consisted of three cultures: Peucetians, Messapians, and Daunians. Among them, the Daunians left a peculiar cultural heritage, with one-of-a-kind stelae and pottery, but, despite the extensive archaeological literature, their origin has been lost to time. In order to shed light on this and to provide a genetic picture of Iron Age Southern Italy, we collected and sequenced human remains from three archaeological sites geographically located in Northern Apulia (the area historically inhabited by Daunians) and radiocarbon dated between 1157 and 275 calBCE. We find that Iron Age Apulian samples are still distant from the genetic variability of modern-day Apulians, they show a degree of genetic heterogeneity comparable with the cosmopolitan Republican and Imperial Roman civilization, even though a few kilometers and centuries separate them, and they are well inserted into the Iron Age Pan-Mediterranean genetic landscape. Our study provides for the first time a window on the genetic make-up of pre-Roman Apulia, whose increasing connectivity within the Mediterranean landscape, would have contributed to laying the foundation for modern genetic variability. In this light, the genetic profile of Daunians may be compatible with an at least partial autochthonous origin, with plausible contributions from the Balkan peninsula.

Genetic continuity of Indo-Iranian speakers since the Iron Age in southern Central Asia

Since prehistoric times, southern Central Asia has been at the crossroads of the movement of people, culture, and goods. Today, the Central Asian populations are divided into two cultural and linguistic groups: the Indo-Iranian and the Turko-Mongolian groups. Previous genetic studies unveiled that migrations from East Asia contributed to the spread of Turko-Mongolian populations in Central Asia and the partial replacement of the Indo-Iranian populations. However, little is known about the origin of the latters. To shed light on this, we compare the genetic data on two current-day Indo-Iranian populations — Yaghnobis and Tajiks — with genome-wide data from published ancient individuals. The present Indo-Iranian populations from Central Asia display a strong genetic continuity with Iron Age samples from Turkmenistan and Tajikistan. We model Yaghnobis as a mixture of 93% Iron Age individual from Turkmenistan and 7% from Baikal. For the Tajiks, we observe a higher Baikal ancestry and an additional admixture event with a South Asian population. Our results, therefore, suggest that in addition to a complex history, Central Asia shows a remarkable genetic continuity since the Iron Age, with only limited gene flow.

Three assays for in-solution enrichment of ancient human DNA at more than a million SNPs

The strategy of in-solution enrichment for hundreds of thousands of single-nucleotide polymorphisms (SNPs) has been used to analyze >70% of individuals with genome-scale ancient DNA published to date. This approach makes it economical to study ancient samples with low proportions of human DNA and increases the rate of conversion of sampled remains into interpretable data. So far, nearly all such data have been generated using a set of bait sequences targeting about 1.24 million SNPs (the "1240k reagent"), but synthesis of the reagent has been cost-effective for only a few laboratories. In 2021, two companies, Daicel Arbor Biosciences and Twist Bioscience, made available assays that target the same core set of SNPs along with supplementary content. We test all three assays on a common set of 27 ancient DNA libraries and show that all three are effective at enriching many hundreds of thousands of SNPs. For all assays, one round of enrichment produces data that are as useful as two. In our testing, the "Twist Ancient DNA" assay produces the highest coverages, greatest uniformity on targeted positions, and almost no bias toward enriching one allele more than another relative to shotgun sequencing. We also identify hundreds of thousands of targeted SNPs for which there is minimal allelic bias when comparing 1240k data to either shotgun or Twist data. This facilitates coanalysis of the large data sets that have been generated using 1240k and Twist capture, as well as shotgun sequencing approaches.

Direct detection of natural selection in Bronze Age Britain

We developed a novel method for efficiently estimating time-varying selection coefficients from genome-wide ancient DNA data. In simulations, our method accurately recovers selective trajectories and is robust to misspecification of population size. We applied it to a large data set of ancient and present-day human genomes from Britain and identified seven loci with genome-wide significant evidence of selection in the past 4500 yr. Almost all of them can be related to increased vitamin D or calcium levels, suggesting strong selective pressure on these or related phenotypes. However, the strength of selection on individual loci varied substantially over time, suggesting that cultural or environmental factors moderated the genetic response. Of 28 complex anthropometric and metabolic traits, skin pigmentation was the only one with significant evidence of polygenic selection, further underscoring the importance of phenotypes related to vitamin D. Our approach illustrates the power of ancient DNA to characterize selection in human populations and illuminates the recent evolutionary history of Britain.

Is History the Same as Evolution? No. Is it Independent of Evolution? Certainly Not

History is full of violence and oppression within and between groups, and although group conflicts enhance within-group cooperation (mediated by oxytocin, which promotes parochial altruism) the hierarchy within groups ensures that spoils accrue very unevenly. Darwin suggested, and we now know, that sexual selection is as powerful as selection by mortality, and the main purpose of survival is reproduction. Male reproductive skew is greater than that among females in all societies, but the difference became much greater after the hunting-gathering era, and the rise of so-called "civilization" was everywhere a process of predatory expansion, producing kingdoms and empires where top males achieved astounding heights of reproductive success. This was shown by historical and ethnographic data now strongly confirmed by genomic science. Psychological research confirms that group identity, out-group stigmatization, leadership characterized by charisma, the will to power, narcissism, sociopathy, and cruelty, and followership characterized by hypnotic obedience, loss of individuality, and cruelty are integral parts of human nature. We can thank at least ten or twelve millennia of microevolutionary processes such as those described above, all more prominent in males than females. Followers in wars have faced a difficult risk-benefit analysis, but if they survived and won they too could increase their reproductive success through the rape and other sexual exploitation that have accompanied almost all wars. For modern leaders, social monogamy and contraception have separated autocracy from reproductive success, but only partly, and current worldwide autocratic trends still depend on the evolved will to power, obedience, and cruelty.

The spread of agriculture in Iberia through Approximate Bayesian Computation and Neolithic projectile tools

In the present article we use geometric microliths (a specific type of arrowhead) and Approximate Bayesian Computation (ABC) in order to evaluate possible origin points and expansion routes for the Neolithic in the Iberian Peninsula. In order to do so, we divide the Iberian Peninsula in four areas (Ebro river, Catalan shores, Xúquer river and Guadalquivir river) and we sample the geometric microliths existing in the sites with the oldest radiocarbon dates for each zone. On this data, we perform a partial Mantel test with three matrices: geographic distance matrix, cultural distance matrix and chronological distance matrix. After this is done, we simulate a series of partial Mantel tests where we alter the chronological matrix by using an expansion model with randomised origin points, and using the distribution of the observed partial Mantel test’s results as a summary statistic within an Approximate Bayesian Computation-Sequential Monte-Carlo (ABC-SMC) algorithm framework. Our results point clearly to a Neolithic expansion route following the Northern Mediterranean, whilst the Southern Mediterranean route could also find support and should be further discussed. The most probable origin points focus on the Xúquer river area.

Large-scale migration into Britain during the Middle to Late Bronze Age

Present-day people from England and Wales harbour more ancestry derived from Early European Farmers (EEF) than people of the Early Bronze Age¹. To understand this, we generated genome-wide data from 793 individuals, increasing data from the Middle to Late Bronze and Iron Age in Britain by 12-fold, and Western and Central Europe by 3.5-fold. Between 1000 and 875 BC, EEF ancestry increased in southern Britain (England and Wales) but not northern Britain (Scotland) due to incorporation of migrants who arrived at this time and over previous centuries, and who were genetically most similar to ancient individuals from France. These migrants contributed about half the ancestry of Iron Age people of England and Wales, thereby creating a plausible vector for the spread of early Celtic languages into Britain. These patterns are part of a broader trend of EEF ancestry becoming more similar across central and western Europe in the Middle to Late Bronze Age, coincident with archaeological evidence of intensified cultural exchange^2-6. There was comparatively less gene flow from continental Europe during the Iron Age, and Britain's independent genetic trajectory is also reflected in the rise of the allele conferring lactase persistence to ~50% by this time compared to ~7% in central Europe where it rose rapidly in frequency only a millennium later. This suggests that dairy products were used in qualitatively different ways in Britain and in central Europe over this period.

A high-resolution picture of kinship practices in an Early Neolithic tomb

To explore kinship practices at chambered tombs in Early Neolithic Britain, here we combined archaeological and genetic analyses of 35 individuals who lived about 5,700 years ago and were entombed at Hazleton North long cairn¹. Twenty-seven individuals are part of the first extended pedigree reconstructed from ancient DNA, a five-generation family whose many interrelationships provide statistical power to document kinship practices that were invisible without direct genetic data. Patrilineal descent was key in determining who was buried in the tomb, as all 15 intergenerational transmissions were through men. The presence of women who had reproduced with lineage men and the absence of adult lineage daughters suggest virilocal burial and female exogamy. We demonstrate that one male progenitor reproduced with four women: the descendants of two of those women were buried in the same half of the tomb over all generations. This suggests that maternal sub-lineages were grouped into branches whose distinctiveness was recognized during the construction of the tomb. Four men descended from non-lineage fathers and mothers who also reproduced with lineage male individuals, suggesting that some men adopted the children of their reproductive partners by other men into their patriline. Eight individuals were not close biological relatives of the main lineage, raising the possibility that kinship also encompassed social bonds independent of biological relatedness.

The genetic scenario of Mercheros: an under-represented group within the Iberian Peninsula

BACKGROUND

The general picture of human genetic variation has been vastly depicted in the last years, yet many populations remain broadly understudied. In this work, we analyze for the first time the Merchero population, a Spanish minority ethnic group that has been scarcely studied and historically persecuted. Mercheros have been roughly characterised by an itinerant history, common traditional occupations, and the usage of their own language.

RESULTS

Here, we examine the demographic history and genetic scenario of Mercheros, by using genome-wide array data, whole mitochondrial sequences, and Y chromosome STR markers from 25 individuals. These samples have been complemented with a wide-range of present-day populations from Western Eurasia and North Africa. Our results show that the genetic diversity of Mercheros is explained within the context of the Iberian Peninsula, evidencing a modest signal of Roma admixture. In addition, Mercheros present low genetic isolation and intrapopulation heterogeneity.

CONCLUSIONS

This study represents the first genetic characterisation of the Merchero population, depicting their fine-scale ancestry components and genetic scenario within the Iberian Peninsula. Since ethnicity is not only influenced by genetic ancestry but also cultural factors, other studies from multiple disciplines are needed to further explore the Merchero population. As with Mercheros, there is a considerable gap of underrepresented populations and ethnic groups in publicly available genetic data. Thus, we encourage the consideration of more ethnically diverse population panels in human genetic studies, as an attempt to improve the representation of human populations and better reconstruct their fine-scale history.