Genetic history from the Middle Neolithic to present on the Mediterranean island of Sardinia

Urbanization and genetic homogenization in the medieval Low Countries revealed through a ten-century paleogenomic study of the city of Sint-Truiden

BACKGROUND

Processes shaping the formation of the present-day population structure in highly urbanized Northern Europe are still poorly understood. Gaps remain in our understanding of when and how currently observable regional differences emerged and what impact city growth, migration, and disease pandemics during and after the Middle Ages had on these processes.

RESULTS

We perform low-coverage sequencing of the genomes of 338 individuals spanning the eighth to the eighteenth centuries in the city of Sint-Truiden in Flanders, in the northern part of Belgium. The early/high medieval Sint-Truiden population was more heterogeneous, having received migrants from Scotland or Ireland, and displayed less genetic relatedness than observed today between individuals in present-day Flanders. We find differences in gene variants associated with high vitamin D blood levels between individuals with Gaulish or Germanic ancestry. Although we find evidence of a Yersinia pestis infection in 5 of the 58 late medieval burials, we were unable to detect a major population-scale impact of the second plague pandemic on genetic diversity or on the elevated differentiation of immunity genes.

CONCLUSIONS

This study reveals that the genetic homogenization process in a medieval city population in the Low Countries was protracted for centuries. Over time, the Sint-Truiden population became more similar to the current population of the surrounding Limburg province, likely as a result of reduced long-distance migration after the high medieval period, and the continuous process of local admixture of Germanic and Gaulish ancestries which formed the genetic cline observable today in the Low Countries.

Performance of qpAdm-based screens for genetic admixture on graph-shaped histories and stepping stone landscapes

qpAdm is a statistical tool that is often used for testing large sets of alternative admixture models for a target population. Despite its popularity, qpAdm remains untested on 2D stepping stone landscapes and in situations with low prestudy odds (low ratio of true to false models). We tested high-throughput qpAdm protocols with typical properties such as number of source combinations per target, model complexity, model feasibility criteria, etc. Those protocols were applied to admixture graph-shaped and stepping stone simulated histories sampled randomly or systematically. We demonstrate that false discovery rates of high-throughput qpAdm protocols exceed 50% for many parameter combinations since: (1) prestudy odds are low and fall rapidly with increasing model complexity; (2) complex migration networks violate the assumptions of the method; hence, there is poor correlation between qpAdm P-values and model optimality, contributing to low but nonzero false-positive rate and low power; and (3) although admixture fraction estimates between 0 and 1 are largely restricted to symmetric configurations of sources around a target, a small fraction of asymmetric highly nonoptimal models have estimates in the same interval, contributing to the false-positive rate. We also reinterpret large sets of qpAdm models from 2 studies in terms of source-target distance and symmetry and suggest improvements to qpAdm protocols: (1) temporal stratification of targets and proxy sources in the case of admixture graph-shaped histories, (2) focused exploration of few models for increasing prestudy odds; and (3) dense landscape sampling for increasing power and stringent conditions on estimated admixture fractions for decreasing the false-positive rate.

High continuity of forager ancestry in the Neolithic period of the eastern Maghreb

Ancient DNA from the Mediterranean region has revealed long-range connections and population transformations associated with the spread of food-producing economies1-6. However, in contrast to Europe, genetic data from this key transition in northern Africa are limited, and have only been available from the far western Maghreb (Morocco)1-3. Here we present genome-wide data for nine individuals from the Later Stone Age through the Neolithic period from Algeria and Tunisia. The earliest individuals cluster with pre-Neolithic people of the western Maghreb (around 15,000-7,600 years before present (BP)), showing that this 'Maghrebi' ancestry profile had a substantial geographic and temporal extent. At least one individual from Djebba (Tunisia), dating to around 8,000 years BP, harboured ancestry from European hunter-gatherers, probably reflecting movement in the Early Holocene across the Strait of Sicily. Later Neolithic people from the eastern Maghreb retained largely local forager ancestry, together with smaller contributions from European farmers (by around 7,000 years BP) and Levantine groups (by around 6,800 years BP), and were thus far less impacted by external gene flow than were populations in other parts of the Neolithic Mediterranean.

Punic people were genetically diverse with almost no Levantine ancestors

The maritime Phoenician civilization from the Levant transformed the entire Mediterranean during the first millennium BCE1-3. However, the extent of human movement between the Levantine Phoenician homeland and Phoenician-Punic settlements in the central and western Mediterranean has been unclear in the absence of comprehensive ancient DNA studies. Here, we generated genome-wide data for 210 individuals, including 196 from 14 sites traditionally identified as Phoenician and Punic in the Levant, North Africa, Iberia, Sicily, Sardinia and Ibiza, and an early Iron Age individual from Algeria. Levantine Phoenicians made little genetic contribution to Punic settlements in the central and western Mediterranean between the sixth and second centuries BCE, despite abundant archaeological evidence of cultural, historical, linguistic and religious links4. Instead, these inheritors of Levantine Phoenician culture derived most of their ancestry from a genetic profile similar to that of Sicily and the Aegean. Much of the remaining ancestry originated from North Africa, reflecting the growing influence of Carthage5. However, this was a minority contributor of ancestry in all of the sampled sites, including in Carthage itself. Different Punic sites across the central and western Mediterranean show similar patterns of high genetic diversity. We also detect genetic relationships across the Mediterranean, reflecting shared demographic processes that shaped the Punic world.

Genetic structure and demographic history of house mice in western Europe inferred using whole-genome sequences

The western house mouse, Mus musculus domesticus, is a human commensal and an outstanding model organism for studying a wide variety of traits and diseases. However, we have few genomic resources for wild mice and only a rudimentary understanding of the demographic history of house mice in Europe. Here, we sequenced 59 whole genomes of mice collected from England, Scotland, Wales, Guernsey, northern France, Italy, Portugal and Spain. We combined this dataset with 24 previously published sequences from southern France, Germany and Iran and compared patterns of population structure and inferred demographic parameters for house mice in western Europe to patterns seen in humans. Principal component and phylogenetic analyses identified three genetic clusters in western European mice. Admixture and f-branch statistics identified historical gene flow between these genetic clusters. Demographic analyses suggest a shared history of population bottlenecks prior to 20 000 years ago. Estimated divergence times between populations of house mice from western Europe ranged from 1500 to 5500 years ago, in general agreement with the zooarchaeological record. These results correspond well with key aspects of contemporary human population structure and the history of migration in western Europe, highlighting the commensal relationship of this important genetic model.

Performance of qpAdm-based screens for genetic admixture on admixture-graph-shaped histories and stepping-stone landscapes

qpAdm is a statistical tool that is often used for testing large sets of alternative admixture models for a target population. Despite its popularity, qpAdm remains untested on two-dimensional stepping-stone landscapes and in situations with low pre-study odds (low ratio of true to false models). We tested high-throughput qpAdm protocols with typical properties such as number of source combinations per target, model complexity, model feasibility criteria, etc. Those protocols were applied to admixture-graph-shaped and stepping-stone simulated histories sampled randomly or systematically. We demonstrate that false discovery rates of high-throughput qpAdm protocols exceed 50% for many parameter combinations since: 1) pre-study odds are low and fall rapidly with increasing model complexity; 2) complex migration networks violate the assumptions of the method, hence there is poor correlation between qpAdm p-values and model optimality, contributing to low but non-zero false positive rate and low power; 3) although admixture fraction estimates between 0 and 1 are largely restricted to symmetric configurations of sources around a target, a small fraction of asymmetric highly non-optimal models have estimates in the same interval, contributing to the false positive rate. We also re-interpret large sets of qpAdm models from two studies in terms of source-target distance and symmetry and suggest improvements to qpAdm protocols: 1) temporal stratification of targets and proxy sources in the case of admixture-graph-shaped histories; 2) focused exploration of few models for increasing pre-study odds; 3) dense landscape sampling for increasing power and stringent conditions on estimated admixture fractions for decreasing the false positive rate.

The extent of the hip bone sexual dimorphism in two Italian coeval modern skeletal samples

The rate of sexual dimorphism in the human hip bone is primarily due to the structural demands of childbirth. Genetic, environmental, and socio-cultural factors can also influence pelvic shape variations across populations. This study examines intra-population sex variation within the Italian population based on regional differences of 280 coxal bones from two documented human osteological collections (Bologna and Sassari) coming from different geographical areas, the northern continental and island regions. Nineteen metric variables were used to evaluate sexual dimorphism and population differences. Most of the variables showed sexual dimorphism, particularly the hip height and ischio-pubic measures within both populations, and accurately predicted sex for each population (Bologna: 100%; Sassari: 91.2%). Results show the Bologna sample have larger dimensions for most of the measurements than the Sassari one, especially when considering the longitudinal ones. Some female traits of the Bologna sample are larger than the correspondent ones in the Sassari males. The rate of sexual dimorphism between the populations shows significant differences, with better sex distinction for Bologna than Sassari. This study aims at interpreting these intra-population differences, considering the effect of environmental (physical and social milieu) and genetic factors, underscoring the importance of this local differences for accurate diagnostic criteria.

Patterns of population structure and genetic variation within the Saudi Arabian population

The Arabian Peninsula is considered the initial site of historic human migration out of Africa. The modern-day indigenous Arabians are believed to be the descendants who remained from the ancient split of the migrants into Eurasia. Here, we investigated how the population history and cultural practices such as endogamy have shaped the genetic variation of the Saudi Arabians. We genotyped 3,352 individuals and identified twelve genetic sub-clusters that corresponded to the geographical distribution of different tribal regions, differentiated by distinct components of ancestry based on comparisons to modern and ancient DNA references. These sub-clusters also showed variation across ranges of the genome covered in runs of homozygosity, as well as differences in population size changes over time. Using 25,488,981 variants found in whole genome sequencing data (WGS) from 302 individuals, we found that the Saudi tend to show proportionally more deleterious alleles than neutral alleles when compared to Africans/African Americans from gnomAD (e.g. a 13% increase of deleterious alleles annotated by AlphaMissense between 0.5 - 5% frequency in Saudi, compared to 7% decrease of the benign alleles; P < 0.001). Saudi sub-clusters with greater inbreeding and lower effective population sizes showed greater enrichment of deleterious alleles as well. Additionally, we found that approximately 10% of the variants discovered in our WGS data are not observed in gnomAD; these variants are also enriched with deleterious annotations. To accelerate studying the population-enriched deleterious alleles and their health consequences in this population, we made available the allele frequency estimates of 25,488,981 variants discovered in our samples. Taken together, our results suggest that Saudi's population history impacts its pattern of genetic variation with potential consequences to the population health. It further highlights the need to sequence diverse and unique populations so to provide a foundation on which to interpret medical- and pharmaco- genomic findings from these populations.

The cranium from the Octagon in Ephesos

During excavations in 1929, a well-preserved skeleton was discovered in a sarcophagus in the Octagon at Ephesos (Turkey). For the following century, archaeologists have speculated about the identity of this obviously notable person. Repeated claim is that the remains could represent Arsinoë IV, daughter of Ptolemy XII, and younger (half-)sister of Cleopatra VII. To address these questions we undertook state-of-the-art morphological, genetic and dating analyses of the cranium and further analyses of bone samples from a femur and a rib of the skeleton found in the same tomb. We confirm based on genetic analyses from the cranium and the femur that they derive from the same person. 14C-dating of the cranium provides a most likely time range between 205-36 BC. The connection with Arsinoë IV can be excluded because we confirmed that the individual is a male. The cranium represents an 11-14-year-old boy who suffered from significant developmental disturbances. Genetics suggest an ancestry from the Italian peninsula or Sardinia. The fate of the body of Arsinoë IV, who reportedly was killed in 41 BC in Ephesos, remains open. In contrast, investigations regarding the fate and social background of the boy from the Octagon can now proceed free of speculation.

The genomic portrait of the Picene culture provides new insights into the Italic Iron Age and the legacy of the Roman Empire in Central Italy

BACKGROUND

The Italic Iron Age is characterized by the presence of various ethnic groups partially examined from a genomic perspective. To explore the evolution of Iron Age Italic populations and the genetic impact of Romanization, we focus on the Picenes, one of the most fascinating pre-Roman civilizations, who flourished on the Middle Adriatic side of Central Italy between the 9th and the 3rd century BCE, until the Roman colonization.

RESULTS

More than 50 samples are reported, spanning more than 1000 years of history from the Iron Age to Late Antiquity. Despite cultural diversity, our analysis reveals no major differences between the Picenes and other coeval populations, suggesting a shared genetic history of the Central Italian Iron Age ethnic groups. Nevertheless, a slight genetic differentiation between populations along the Adriatic and Tyrrhenian coasts can be observed, possibly due to different population dynamics in the two sides of Italy and/or genetic contacts across the Adriatic Sea. Additionally, we identify several individuals with ancestries deviating from their general population. Lastly, in our Late Antiquity site, we observe a drastic change in the genetic landscape of the Middle Adriatic region, indicating a relevant influx from the Near East, possibly as a consequence of Romanization.

CONCLUSIONS

Our findings, consistently with archeological hypotheses, suggest genetic interactions across the Adriatic Sea during the Bronze/Iron Age and a high level of individual mobility typical of cosmopolitan societies. Finally, we highlight the role of the Roman Empire in shaping genetic and phenotypic changes that greatly impact the Italian peninsula.

Ancient DNA challenges prevailing interpretations of the Pompeii plaster casts

The eruption of Somma-Vesuvius in 79 CE buried several nearby Roman towns, killing the inhabitants and burying under pumice lapilli and ash deposits a unique set of civil and private buildings, monuments, sculptures, paintings, and mosaics that provide a rich picture of life in the empire. The eruption also preserved the forms of many of the dying as the ash compacted around their bodies. Although the soft tissue decayed, the outlines of the bodies remained and were recovered by excavators centuries later by filling the cavities with plaster. From skeletal material embedded in the casts, we generated genome-wide ancient DNA and strontium isotopic data to characterize the genetic relationships, sex, ancestry, and mobility of five individuals. We show that the individuals' sexes and family relationships do not match traditional interpretations, exemplifying how modern assumptions about gendered behaviors may not be reliable lenses through which to view data from the past. For example, an adult wearing a golden bracelet with a child on their lap-often interpreted as mother and child-is genetically an adult male biologically unrelated to the child. Similarly, a pair of individuals who were thought to have died in an embrace-often interpreted as sisters-included at least one genetic male. All Pompeiians with genome-wide data consistently derive their ancestry largely from recent immigrants from the eastern Mediterranean, as has also been seen in contemporaneous ancient genomes from the city of Rome, underscoring the cosmopolitanism of the Roman Empire in this period.

Genetic ancestry in Puerto Rican afro-descendants illustrates diverse histories of African diasporic populations

OBJECTIVES

Genetic studies of contemporary Puerto Ricans reflect a demographic history characterized by admixture between Indigenous American, African, and European peoples. While previous studies provide genetic perspectives on the general Puerto Rican population, less is known about the island's sub-populations, specifically Afro-Puerto Ricans.

MATERIALS AND METHODS

In this study, the genetic ancestry of Afro-Puerto Ricans is characterized and compared to other Caribbean populations. Thirty DNA samples collected among self-identified Puerto Ricans of African descent in Loíza (n = 2), Piñones (n = 13), San Juan (n = 2), Mayagüez (n = 9), and Ponce (n = 4), were genotyped at 750,000 loci on the National Geographic Genochip. We then applied unsupervised clustering and dimensionality-reduction methods to detect continental and subcontinental African and European genetic ancestry patterns.

RESULTS

Admixture analyses reveal that on average, the largest genetic ancestry component for Afro-Puerto Ricans is African in origin, followed by European and Indigenous American genetic ancestry components. African biogeographic origins of Afro-Puerto Ricans align most closely with contemporary peoples of Lower Guinea and the Bight of Biafra, while the European genetic ancestry component is most similar to contemporary Iberian, Italian, and Basque populations. These findings contrast with the biogeographic origins of comparative Barbadian and Puerto Rican populations.

DISCUSSION

Our results suggest that while there are similarities with regard to general patterns of genetic ancestry among African descendants in the Caribbean, there is previously unrecognized regional heterogeneity, including among Puerto Rican sub-populations. These results are also consistent with available historical sources, while providing depth absent from the documentary record, particularly with regard to African ancestry.

Analysis of Gyimes Csango population samples on a high-resolution genome-wide basis

BACKGROUND

The Csangos are an East-Central European ethnographic group living mainly in east of Transylvania in Romania. Traditionally, ethnography distinguishes three Csango subpopulations, the Moldavian, Gyimes and Burzenland Csangos. In our previous study we found that the Moldavian Csangos have East Asian/Siberian Turkic ancestry components that might be unique in the East-Central European region and might help to better understand the history of Hungarian speaking ethnic groups of the area. Since then, we obtained further Csango samples from Moldavia and from a distinct region of Gyimes, which two Csango subgroups are traditionally different since they live in a degree of isolation not only from other people but also from each other. Here we present the first genomic analysis of Gyimes Csangos, which intended to compare the genetic makeup of these two Csango subgroups using both allele-frequency and haplotype-based methods. The main goal of the study was to investigate the genetic isolation of the Csangos on a genome-wide SNP basis and to assess the isolation of Gyimes Csangos, which in contrast to the Moldavians was not yet studied.

RESULTS

Our results show that these two Csango groups show slight differences from each other. We confirmed the genetic isolation of Moldavian Csangos and revealed that Gyimes Csangos have a similar, but detectably weaker isolation. In the case of Gyimes Csangos we detected also a stronger East European or presumably Asian derived ancestry.

CONCLUSION

The Gyimes Csangos show a degree of genetic isolation comparable to that of the Moldavians. The Asian ancestry that differentiates the Moldavian Csango people from the other East-Central European populations may be present in the Gyimes Csangos in an even higher degree, since Gyimes Csango individuals show a more significant share from that ancestry component.

Pervasive findings of directional selection realize the promise of ancient DNA to elucidate human adaptation

We present a method for detecting evidence of natural selection in ancient DNA time-series data that leverages an opportunity not utilized in previous scans: testing for a consistent trend in allele frequency change over time. By applying this to 8433 West Eurasians who lived over the past 14000 years and 6510 contemporary people, we find an order of magnitude more genome-wide significant signals than previous studies: 347 independent loci with >99% probability of selection. Previous work showed that classic hard sweeps driving advantageous mutations to fixation have been rare over the broad span of human evolution, but in the last ten millennia, many hundreds of alleles have been affected by strong directional selection. Discoveries include an increase from ~0% to ~20% in 4000 years for the major risk factor for celiac disease at HLA-DQB1; a rise from ~0% to ~8% in 6000 years of blood type B; and fluctuating selection at the TYK2 tuberculosis risk allele rising from ~2% to ~9% from ~5500 to ~3000 years ago before dropping to ~3%. We identify instances of coordinated selection on alleles affecting the same trait, with the polygenic score today predictive of body fat percentage decreasing by around a standard deviation over ten millennia, consistent with the "Thrifty Gene" hypothesis that a genetic predisposition to store energy during food scarcity became disadvantageous after farming. We also identify selection for combinations of alleles that are today associated with lighter skin color, lower risk for schizophrenia and bipolar disease, slower health decline, and increased measures related to cognitive performance (scores on intelligence tests, household income, and years of schooling). These traits are measured in modern industrialized societies, so what phenotypes were adaptive in the past is unclear. We estimate selection coefficients at 9.9 million variants, enabling study of how Darwinian forces couple to allelic effects and shape the genetic architecture of complex traits.

Testing times: disentangling admixture histories in recent and complex demographies using ancient DNA

Our knowledge of human evolutionary history has been greatly advanced by paleogenomics. Since the 2020s, the study of ancient DNA has increasingly focused on reconstructing the recent past. However, the accuracy of paleogenomic methods in resolving questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation remains an open question. We evaluated the performance and behavior of two commonly used methods, qpAdm and the f3-statistic, on admixture inference under a diversity of demographic models and data conditions. We performed two complementary simulation approaches-firstly exploring a wide demographic parameter space under four simple demographic models of varying complexities and configurations using branch-length data from two chromosomes-and secondly, we analyzed a model of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudohaploidization. We observe that population differentiation is the primary factor driving qpAdm performance. Notably, while complex gene flow histories influence which models are classified as plausible, they do not reduce overall performance. Under conditions reflective of the historical period, qpAdm most frequently identifies the true model as plausible among a small candidate set of closely related populations. To increase the utility for resolving fine-scaled hypotheses, we provide a heuristic for further distinguishing between candidate models that incorporates qpAdm model P-values and f3-statistics. Finally, we demonstrate a significant performance increase for qpAdm using whole-genome branch-length f2-statistics, highlighting the potential for improved demographic inference that could be achieved with future advancements in f-statistic estimations.

Low Genetic Impact of the Roman Occupation of Britain in Rural Communities

The Roman period saw the empire expand across Europe and the Mediterranean, including much of what is today Great Britain. While there is written evidence of high mobility into and out of Britain for administrators, traders, and the military, the impact of imperialism on local, rural population structure, kinship, and mobility is invisible in the textual record. The extent of genetic change that occurred in Britain during the Roman military occupation remains underexplored. Here, using genome-wide data from 52 ancient individuals from eight sites in Cambridgeshire covering the period of Roman occupation, we show low levels of genetic ancestry differentiation between Romano-British sites and indications of larger populations than in the Bronze Age and Neolithic. We find no evidence of long-distance migration from elsewhere in the Empire, though we do find one case of possible temporary mobility within a family unit during the Late Romano-British period. We also show that the present-day patterns of genetic ancestry composition in Britain emerged after the Roman period.

READv2: advanced and user-friendly detection of biological relatedness in archaeogenomics

The advent of genome-wide ancient DNA analysis has revolutionized our understanding of prehistoric societies. However, studying biological relatedness in these groups requires tailored approaches due to the challenges of analyzing ancient DNA. READv2, an optimized Python3 implementation of the most widely used tool for this purpose, addresses these challenges while surpassing its predecessor in speed and accuracy. For sufficient amounts of data, it can classify up to third-degree relatedness and differentiate between the two types of first-degree relatedness, full siblings and parent-offspring. READv2 enables user-friendly, efficient, and nuanced analysis of biological relatedness, facilitating a deeper understanding of past social structures.

Epidemiology of seropositive myasthenia gravis in Sardinia: A population-based study in the district of Sassari

INTRODUCTION/AIMS

The global incidence and prevalence of myasthenia gravis (MG) range between 6-31/million and 10-37/100,000, respectively. Sardinia is a high-risk region for different immune-mediated disorders, but the epidemiology of MG remains unclear. We determined the epidemiology of MG with acetylcholine receptor (AChR)-immunoglobulin G (IgG) and muscle-specific tyrosine kinase (MuSK)-IgG in the district of Sassari (North-Western Sardinia; population, 325,288).

METHODS

From the laboratory of the University Hospital of Sassari (reference for AChR/MuSK-IgG testing in Sardinia since 1998) and the main neurology units in Sardinia, we retrospectively identified MG patients with (1) AChR-IgG and/or MuSK-IgG positivity by radioimmunoprecipitation assay; and (2) residency in the district of Sassari. Incidence (January 2010-December 2019) and prevalence (December 31, 2019) were calculated.

RESULTS

A total of 202 patients were included (incident, 107; prevalent, 180). Antibody specificities were AChR (n = 187 [93%]) and MuSK (n = 15 [7%]). The crude MG incidence (95% confidence interval) was 32.6 (26.8-39.2)/million, while prevalence was 55.3 (47.7-63.9)/100,000. After age-standardization to the world population, incidence decreased to 18.4 (14.3-22.5)/million, while prevalence decreased to 31.6 (26.1-37.0)/100,000. Among incident cases, age strata (years) at MG onset were: <18 (2%), 18-49 (14%), 50-64 (21%), and ≥65 (63%).

DISCUSSION

Sardinia is a high-risk region for MG, with a prevalence that exceeds the European threshold for rare disease. Identification of the environmental and genetic determinants of this risk may improve our understanding of disease pathophysiology.

Echocardiographic heart ageing patterns predict cardiovascular and non-cardiovascular events and reflect biological age: the SardiNIA study

AIMS

Age is a crucial risk factor for cardiovascular (CV) and non-CV diseases. As people age at different rates, the concept of biological age has been introduced as a personalized measure of functional deterioration. Associations of age with echocardiographic quantitative traits were analysed to assess different heart ageing rates and their ability to predict outcomes and reflect biological age.

METHODS AND RESULTS

Associations of age with left ventricular mass, geometry, diastolic function, left atrial volume, and aortic root size were measured in 2614 healthy subjects. Based on the 95% two-sided tolerance intervals of each correlation, three discrete ageing trajectories were identified and categorized as 'slow', 'normal', and 'accelerated' heart ageing patterns. The primary endpoint included fatal and non-fatal CV events, and the secondary endpoint was a composite of CV and non-CV events and all-cause death. The phenotypic age of the heart (HeartPhAge) was estimated as a proxy of biological age. The slow ageing pattern was found in 8.7% of healthy participants, the normal pattern in 76.9%, and the accelerated pattern in 14.4%. Kaplan-Meier curves of the heart ageing patterns diverged significantly (P = 0.0001) for both primary and secondary endpoints, with the event rate being lowest in the slow, intermediate in the normal, and highest in the accelerated pattern. In the Cox proportional hazards model, heart ageing patterns predicted both primary (P = 0.01) and secondary (P = 0.03 to <0.0001) endpoints, independent of chronological age and risk factors. Compared with chronological age, HeartPhAge was 9 years younger in slow, 4 years older in accelerated (both P < 0.0001), and overlapping in normal ageing patterns.

CONCLUSION

Standard Doppler echocardiography detects slow, normal, and accelerated heart ageing patterns. They predict CV and non-CV events, reflect biological age, and provide a new tool to calibrate prevention timing and intensity.

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.

Socio-cultural practices may have affected sex differences in stature in Early Neolithic Europe

The rules and structure of human culture impact health as much as genetics or environment. To study these relationships, we combine ancient DNA (n = 230), skeletal metrics (n = 391), palaeopathology (n = 606) and dietary stable isotopes (n = 873) to analyse stature variation in Early Neolithic Europeans from North Central, South Central, Balkan and Mediterranean regions. In North Central Europe, stable isotopes and linear enamel hypoplasias indicate high environmental stress across sexes, but female stature is low, despite polygenic scores identical to males, and suggests that cultural factors preferentially supported male recovery from stress. In Mediterranean populations, sexual dimorphism is reduced, indicating male vulnerability to stress and no strong cultural preference for males. Our analysis indicates that biological effects of sex-specific inequities can be linked to cultural influences at least as early as 7,000 yr ago, and culture, more than environment or genetics, drove height disparities in Early Neolithic Europe.

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.

Prevalence of epilepsy in childhood: An epidemiological study in Sardinia

BACKGROUND

The aim of this study was to investigate the frequency and characteristics of pediatric epilepsy in the geographic isolate of Sardinia island and to calculate the prevalence of active epilepsy.

METHODS

The study was retrospective, observational and involved a systematic review of medical records and computerized archives containing all clinical and EEG recordings of patients with epilepsy referred to the regional structures that could have followed patients with epilepsy in South Sardinia, during the period 2003-2021.

RESULTS

The study population included 112,912 children and adolescents (age ≤ 18 years). 618 children and adolescents (women 42.4 %) were identified. Family history of epilepsy was reported in 153 (26.1 %). Etiology was genetic in 64.5 % and structural in 26.7 % subjects. Focal seizures were reported in 51.6 % of subjects, followed by 34.7 % with generalized seizures and 10.6 % of patients experienced both type of seizures. A total of 301 subjects with active epilepsy in 2019 were identified resulting in a prevalence of 2.67 per 1000 (95 % CI 2.37-2.97). Prevalence in the age class 5-14 years was 4.21 per 1000 (95 % CI 3.72-4.76).

CONCLUSION

Compared to the previous studies in distinct geographic isolates, the present study showed a significantly low prevalence rate of active epilepsy; a high percentage of focal seizures and genetic etiology.

Testing Times: Challenges in Disentangling Admixture Histories in Recent and Complex Demographies

Paleogenomics has expanded our knowledge of human evolutionary history. Since the 2020s, the study of ancient DNA has increased its focus on reconstructing the recent past. However, the accuracy of paleogenomic methods in answering questions of historical and archaeological importance amidst the increased demographic complexity and decreased genetic differentiation within the historical period remains an open question. We used two simulation approaches to evaluate the limitations and behavior of commonly used methods, qpAdm and the f3-statistic, on admixture inference. The first is based on branch-length data simulated from four simple demographic models of varying complexities and configurations. The second, an analysis of Eurasian history composed of 59 populations using whole-genome data modified with ancient DNA conditions such as SNP ascertainment, data missingness, and pseudo-haploidization. We show that under conditions resembling historical populations, qpAdm can identify a small candidate set of true sources and populations closely related to them. However, in typical ancient DNA conditions, qpAdm is unable to further distinguish between them, limiting its utility for resolving fine-scaled hypotheses. Notably, we find that complex gene-flow histories generally lead to improvements in the performance of qpAdm and observe no bias in the estimation of admixture weights. We offer a heuristic for admixture inference that incorporates admixture weight estimate and P-values of qpAdm models, and f3-statistics to enhance the power to distinguish between multiple plausible candidates. Finally, we highlight the future potential of qpAdm through whole-genome branch-length f2-statistics, demonstrating the improved demographic inference that could be achieved with advancements in f-statistic estimations.

Population genomics identifies genetic signatures of carrot domestication and improvement and uncovers the origin of high-carotenoid orange carrots

Here an improved carrot reference genome and resequencing of 630 carrot accessions were used to investigate carrot domestication and improvement. The study demonstrated that carrot was domesticated during the Early Middle Ages in the region spanning western Asia to central Asia, and orange carrot was selected during the Renaissance period, probably in western Europe. A progressive reduction of genetic diversity accompanied this process. Genes controlling circadian clock/flowering and carotenoid accumulation were under selection during domestication and improvement. Three recessive genes, at the REC, Or and Y2 quantitative trait loci, were essential to select for the high α- and β-carotene orange phenotype. All three genes control high α- and β-carotene accumulation through molecular mechanisms that regulate the interactions between the carotenoid biosynthetic pathway, the photosynthetic system and chloroplast biogenesis. Overall, this study elucidated carrot domestication and breeding history and carotenoid genetics at a molecular level.

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.

A genetic history of continuity and mobility in the Iron Age central Mediterranean

The Iron Age was a dynamic period in central Mediterranean history, with the expansion of Greek and Phoenician colonies and the growth of Carthage into the dominant maritime power of the Mediterranean. These events were facilitated by the ease of long-distance travel following major advances in seafaring. We know from the archaeological record that trade goods and materials were moving across great distances in unprecedented quantities, but it is unclear how these patterns correlate with human mobility. Here, to investigate population mobility and interactions directly, we sequenced the genomes of 30 ancient individuals from coastal cities around the central Mediterranean, in Tunisia, Sardinia and central Italy. We observe a meaningful contribution of autochthonous populations, as well as highly heterogeneous ancestry including many individuals with non-local ancestries from other parts of the Mediterranean region. These results highlight both the role of local populations and the extreme interconnectedness of populations in the Iron Age Mediterranean. By studying these trans-Mediterranean neighbours together, we explore the complex interplay between local continuity and mobility that shaped the Iron Age societies of the central Mediterranean.

The genomic history of the indigenous people of the Canary Islands

The indigenous population of the Canary Islands, which colonized the archipelago around the 3rd century CE, provides both a window into the past of North Africa and a unique model to explore the effects of insularity. We generate genome-wide data from 40 individuals from the seven islands, dated between the 3rd-16rd centuries CE. Along with components already present in Moroccan Neolithic populations, the Canarian natives show signatures related to Bronze Age expansions in Eurasia and trans-Saharan migrations. The lack of gene flow between islands and constant or decreasing effective population sizes suggest that populations were isolated. While some island populations maintained relatively high genetic diversity, with the only detected bottleneck coinciding with the colonization time, other islands with fewer natural resources show the effects of insularity and isolation. Finally, consistent genetic differentiation between eastern and western islands points to a more complex colonization process than previously thought.

Climbing the Giara: A quantitative reassessment of movement and visibility in the Nuragic landscape of the Gesturi plateau (South-Central Sardinia, Italy)

The built landscape of Nuragic Sardinia is an exceptional case for geostatistical analysis, allowing for a discussion of long-held assumptions and ideas. The function of nuraghi (ca. 1700-1100 BCE), the most prominent settled monument of the Sardinian Bronze Age, has been addressed via a multiplicity of landscape approaches, mainly relying on intuitive assessments of their spatial properties: nuraghi were assumed as means of territorial control. The series of nuraghi crowning the mesa plateau named Giara of Gesturi (South-Central Sardinia) provides a significant case for the study of their relations with visibility and movement. Context-oriented GIS models based on viewsheds and least-cost paths have been devised as targeted tools. The results show a certain correlation between nuraghi and potential movement on the slopes, thanks to the selection of plateau morphologies such as outward crests. Anyway, nuraghi do not stand exactly at the most accessible points of the plateau. Nuraghi offered ample visual control, especially at large distances, but not specifically over the closest accessible ways. This suggests that the function of nuraghi is somehow connected to defense and visibility, but it is not explained directly by local territorial control: a role as landmarks and multifaceted monuments has likely to be envisioned.

Bioarchaeological and paleogenomic profiling of the unusual Neolithic burial from Grotta di Pietra Sant'Angelo (Calabria, Italy)

The Neolithic burial of Grotta di Pietra Sant'Angelo (CS) represents a unique archaeological finding for the prehistory of Southern Italy. The unusual placement of the inhumation at a rather high altitude and far from inhabited areas, the lack of funerary equipment and the prone deposition of the body find limited similarities in coeval Italian sites. These elements have prompted wider questions on mortuary customs during the prehistory of Southern Italy. This atypical case requires an interdisciplinary approach aimed to build an integrated bioarchaeological profile of the individual. The paleopathological investigation of the skeletal remains revealed the presence of numerous markers that could be associated with craft activities, suggesting possible interpretations of the individual's lifestyle. CT analyses, carried out on the maxillary bones, showed the presence of a peculiar type of dental wear, but also a good density of the bone matrix. Biomolecular and micromorphological analyses of dental calculus highlight the presence of a rich Neolithic-like oral microbiome, the composition of which is consistent with the presence pathologies. Finally, paleogenomic data obtained from the individual were compared with ancient and modern Mediterranean populations, including unpublished high-resolution genome-wide data for 20 modern inhabitants of the nearby village of San Lorenzo Bellizzi, which provided interesting insights into the biodemographic landscape of the Neolithic in Southern Italy.

Machine learning-based Sr isoscape of southern Sardinia: A tool for bio-geographic studies at the Phoenician-Punic site of Nora

Since prehistoric times, the island of Sardinia-in the western Mediterranean-has played a leading role in the dynamics of human population and mobility, in the circulation of raw materials and artefacts, idioms and customs, of technologies and ideas that have enriched the biological, linguistic and cultural heritage of local groups. For the Phoenician and Punic periods (from the 9th to the 3rd centuries BCE), the ancient site of Nora-in southern Sardinia-represents an emblematic case in the study of migratory phenomena that occurred on the Island from the Iron Age until the Roman conquest. Despite the importance of exploring (and characterising) such movements from a wider bio-cultural perspective, the application of bio-geochemical tools for geographical provenance to the ancient skeletal populations of Sardinia is yet scarce. The present work is the first step towards filling this gap with the development of the first isoscape of southern Sardinia using new bioavailable Sr isotope data and a machine-learning approach. From a geolithological point of view, Sardinia is rather heterogeneous and requires detailed studies to correctly assess the distribution of the isotopic signature of bioavailable Sr. The random forest model employed here to construct the Sr isoscape uses several external environmental and geological variables. The most important predictors are related to age and bedrock type, with additional input from local soil properties. A 10-fold cross-validation gives a mean square error of 0.0008 and an R-squared of 0.81, so the model correctly predicts the 87Sr/86Sr ratio of unknown areas. By using a Bayesian provenance assignment workflow, we tested the isoscape here produced to determine the geographic origin and the mobility of archaeological and modern fauna collected from the Phoenician-Punic site of Nora and the surrounding Pula Plain. Our results indicate that archaeological sheep and goats (87Sr/86Sr < 0.7090) are compatible with areas close to Nora and Pula Plain, in agreement with archaeological evidence of pastoralism in those areas. Modern wild and domesticated fauna (87Sr/86Sr > 0.7090) show compatibility with several natural and anthropogenic locations in southern Sardinia, as expected based on modern species distribution data. Finally, we discuss the large Sr isotopic variability of the Nora baseline, where human mobility studies of human cremated and inhumed individuals are currently underway.

Northwest African Neolithic initiated by migrants from Iberia and Levant

In northwestern Africa, lifestyle transitioned from foraging to food production around 7,400 years ago but what sparked that change remains unclear. Archaeological data support conflicting views: (1) that migrant European Neolithic farmers brought the new way of life to North Africa1-3 or (2) that local hunter-gatherers adopted technological innovations4,5. The latter view is also supported by archaeogenetic data6. Here we fill key chronological and archaeogenetic gaps for the Maghreb, from Epipalaeolithic to Middle Neolithic, by sequencing the genomes of nine individuals (to between 45.8- and 0.2-fold genome coverage). Notably, we trace 8,000 years of population continuity and isolation from the Upper Palaeolithic, via the Epipaleolithic, to some Maghrebi Neolithic farming groups. However, remains from the earliest Neolithic contexts showed mostly European Neolithic ancestry. We suggest that farming was introduced by European migrants and was then rapidly adopted by local groups. During the Middle Neolithic a new ancestry from the Levant appears in the Maghreb, coinciding with the arrival of pastoralism in the region, and all three ancestries blend together during the Late Neolithic. Our results show ancestry shifts in the Neolithization of northwestern Africa that probably mirrored a heterogeneous economic and cultural landscape, in a more multifaceted process than observed in other regions.

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.

Isotopic evidence for human mobility in late antique Bulla Regia (Tunisia)

This paper represents the first isotopic study on late antique human mobility in North Africa, using the urban site of Bulla Regia in Tunisia as a case study. We also present the first values for bioavailable ⁸⁷Sr/⁸⁶Sr in northern Tunisia, analysing 63 plant and snail samples, as well as a simple method for the pre-processing of plants in the field to facilitate their export. Bulla Regia was a prominent Roman and late antique town situated on an important axis of transport and communication in North Africa and is therefore an ideal site to explore mobility in the region during this time period. Strontium (⁸⁷Sr/⁸⁶Sr) and oxygen (δ¹⁸O_Carb) isotopic analysis of 22 late antique individuals from a Christian church and cemetery identified at least seven or eight non-locals, while comparative analysis of five Roman individuals from a funerary enclosure on the same site classified all but one of them as potential locals. Most non-local individuals exhibit ⁸⁷Sr/⁸⁶Sr values that match various areas of northern Tunisia, which supports regional mobility rather than long-distance migration, although when combined with the oxygen results, inter-regional mobility from an area with a warmer climate may be hypothesised for some individuals. Examination of the spatial distribution of non-local individuals in their cemetery setting reveals that they were privileged individuals, thus they may reflect the mobility of wealthier town-dwellers in late antiquity, particularly perhaps along the Carthage-Hippo route.

Phylogeographic and population genetic structure of hound-like native dogs of the Mediterranean Basin

The dog was probably the first domesticated animal. Despite extensive archaeological and genetic investigations, the origin and the evolution of the extant dogs are still being debated. Dog breeds that have over time been selected for hunting share common ancestral traits. This study represents the first comprehensive attempt to survey at the genomic and mitochondrial level eight hound-like dogs breeds indigenous to the Mediterranean Basin to determine if they share common ancient origins. Results from the microsatellite analysis indicate that all the dog populations have a low inbreeding value.The Kelb tal-Fenek has a high divergence from the current Egyptian street population, however there is not enough evidence from this study to exclude completely the potential of an ancient common relationship. Overall, the mitochondrial results indicate high frequencies of haplogroups A and B and a low representation of haplogroup C, while only one Egyptian dog could be assigned to haplogroup D. Results reveal identities and shared clades, suggesting the conservation of ancient European mitotypes in the Mediterranean hound-like breeds, especially in the Egyptian population. Although none of the dog populations/breeds participating in this study indicate to be direct descendants of the Egyptian dogs, they still have a very close morphologically resemblance to those iconic Egyptian dogs often depicted in ancient art forms and share some genetic links with the current Egyptian population. Further research is required with other markers such us complete mitogenomes and SNP panels to confirm the complex history of the Mediterranean dogs involved in this study.

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.

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

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.

Genetic Variation among Pharmacogenes in the Sardinian Population

Pharmacogenetics (PGx) aims to identify the genetic factors that determine inter-individual differences in response to drug treatment maximizing efficacy while decreasing the risk of adverse events. Estimating the prevalence of PGx variants involved in drug response, is a critical preparatory step for large-scale implementation of a personalized medicine program in a target population. Here, we profiled pharmacogenetic variation in fourteen clinically relevant genes in a representative sample set of 1577 unrelated sequenced Sardinians, an ancient island population that accounts for genetic variation in Europe as a whole, and, at the same time is enriched in genetic variants that are very rare elsewhere. To this end, we used PGxPOP, a PGx allele caller based on the guidelines created by the Clinical Pharmacogenetics Implementation Consortium (CPIC), to identify the main phenotypes associated with the PGx alleles most represented in Sardinians. We estimated that 99.43% of Sardinian individuals might potentially respond atypically to at least one drug, that on average each individual is expected to have an abnormal response to about 17 drugs, and that for 27 drugs the fraction of the population at risk of atypical responses to therapy is more than 40%. Finally, we identified 174 pharmacogenetic variants for which the minor allele frequency was at least 10% higher among Sardinians as compared to other European populations, a fact that may contribute to substantial interpopulation variability in drug response phenotypes. This study provides baseline information for further large-scale pharmacogenomic investigations in the Sardinian population and underlines the importance of PGx characterization of diverse European populations, such as Sardinians.

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.

Temporal population structure, a genetic dating method for ancient Eurasian genomes from the past 10,000 years

Radiocarbon dating is the gold standard in archeology to estimate the age of skeletons, a key to studying their origins. Many published ancient genomes lack reliable and direct dates, which results in obscure and contradictory reports. We developed the temporal population structure (TPS), a DNA-based dating method for genomes ranging from the Late Mesolithic to today, and applied it to 3,591 ancient and 1,307 modern Eurasians. TPS predictions aligned with the known dates and correctly accounted for kin relationships. TPS dating of poorly dated Eurasian samples resolved conflicting reports in the literature, as illustrated by one test case. We also demonstrated how TPS improved the ability to study phenotypic traits over time. TPS can be used when radiocarbon dating is unfeasible or uncertain or to develop alternative hypotheses for samples younger than 10,000 years ago, a limitation that may be resolved over time as ancient data accumulate.

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

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

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 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.

Assessing temporal and geographic contacts across the Adriatic Sea through the analysis of genome-wide data from Southern Italy

Southern Italy was characterised by a complex prehistory that started with different Palaeolithic cultures, later followed by the Neolithization and the demic dispersal from the Pontic-Caspian Steppe during the Bronze Age. Archaeological and historical evidences point to a link between Southern Italians and the Balkans still present in modern times. To shed light on these dynamics, we analysed around 700 South Mediterranean genomes combined with informative ancient DNAs. Our findings revealed high affinities of South-Eastern Italians with modern Eastern Peloponnesians, and a closer affinity of ancient Greek genomes with those from specific regions of South Italy than modern Greek genomes. The higher similarity could be associated with a Bronze Age component ultimately originating from the Caucasus with high Iranian and Anatolian Neolithic ancestries. Furthermore, extremely differentiated allele frequencies among Northern and Southern Italy revealed putatively adapted SNPs in genes involved in alcohol metabolism, nevi features and immunological traits.

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.

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.