Pleistocene Mitochondrial Genomes Suggest a Single Major Dispersal of Non-Africans and a Late Glacial Population Turnover in Europe

Parallels and discrepancies between non-native species introductions and human migration

Biological invasions and human migrations have increased globally due to socio-economic drivers and environmental factors that have enhanced cultural, economic, and geographic connectivity. Both processes involve the movement, establishment, and spread of species, yet unfold within fundamentally different philosophical, social and biological contexts. Hence, studying biological invasions (invasion science) and human migration (migration studies) presents complex parallels that are potentially fruitful to explore. Here, we examined nuanced parallels and differences between these two phenomena, integrating historical, socio-political, and ethical perspectives. Our review underscores the need for context-specific approaches in policymaking and governance to address effectively the challenges and opportunities of human migration and harm from biological invasions. We suggest that approaches to studying the drivers of biological invasions and human migration provide an excellent opportunity for transdisciplinary research; one that acknowledges the complexities and potential insights from both fields of study. Ultimately, integrating natural and social sciences offers a promising avenue for enriching the understanding of invasion biology and migration dynamics while pursuing just, equitable, and sustainable solutions. However, while human migration is a clear driver of biological invasions, drawing on principles from biological invasions to understand past and current human migration risks oversimplification and the potential for harmful generalisations that disregard the intrinsic rights and cultural dynamics of human migrations. By doing so, we provide insights and frameworks to support the development of context-specific policies that respect human dignity, foster cultural diversity, and address migration challenges in ways that promote global cooperation and justice. This interdisciplinary approach highlights the potential for transdisciplinary research that acknowledges complexities in both fields, ultimately enriching our understanding of invasion biology and migration dynamics while pursuing equitable and sustainable solutions.

Why humans evolved blue eyes

A surprising number of humans are equipped with a subpar eye model-featuring pale, colorful irides that are nowhere as good as the original dark ones at guarding the retina from sunlight and do, in fact, raise one's risk of eye disease. Here I apply evolutionary theory to understand why. I propose that the allele for human blue eyes, which arose just once, managed to spread from one individual to millions at an astonishing speed because it is a greenbeard. "Greenbeards"-imaginary genes, or groups of genes, that produce both a green beard and a behavior that favors other bearers of a green beard-have been deemed exceedingly unlikely to show up in the real world. And yet, as individuals who prefer blue eyes are more inclined to mate with blue-eyed partners and invest in blue-eyed offspring, any blue-eye preference (whether random or arising from the bias for colorful stimuli shared by all recognition systems) becomes rapidly linked to the blue-eye trait. Thus, blue eyes gain an edge by working like a peacock's colorful tail and a nestling's colorful mouth: twice self-reinforcing, "double runaway" evolution via sexual and parental selection. The blue-eye ornament gene, by binding to a behavior that favors other bearers of the blue-eye ornament gene, is ultimately recognizing and helping copies of itself in both kin and strangers-and greatly prospering, just like theory predicts.

Potential Adaptive Introgression From Dogs in Iberian Grey Wolves (Canis lupus)

Invading species along with increased anthropogenization may lead to hybridization events between wild species and closely related domesticates. As a consequence, wild species may carry introgressed alleles from domestic species, which is generally assumed to yield adverse effects in wild populations. The opposite evolutionary consequence, adaptive introgression, where introgressed genes are positively selected in the wild species, is possible but has rarely been documented. Grey wolves (Canis lupus) are widely distributed across the Holarctic and frequently coexist with their close relative, the domestic dog (C. familiaris). Despite ample opportunity, hybridization rarely occurs in most populations. Here we studied the geographically isolated grey wolves of the Iberian Peninsula, who have coexisted with a large population of loosely controlled dogs for thousands of years in a human-modified landscape. We assessed the extent and impact of dog introgression on the current Iberian grey wolf population by analysing 150 whole genomes of Iberian and other Eurasian grey wolves as well as dogs originating from across Europe and western Siberia. We identified almost no recent introgression and a small (< 5%) overall ancient dog ancestry. Using a combination of single scan statistics and ancestry enrichment estimates, we identified positive selection on six genes (DAPP1, NSMCE4A, MPPED2, PCDH9, MBTPS1, and CDH13) for which wild Iberian wolves carry alleles introgressed from dogs. The genes with introgressed and positively selected alleles include functions in immune response and brain functions, which may explain some of the unique behavioural phenotypes in Iberian wolves such as their reduced dispersal compared to other wolf populations.

Persistence of Ancestral KhoeSan Mitochondrial Patterns in Contemporary South African Populations

INTRODUCTION

Southern Africa has been inhabited by hunter-gatherers for at least 20,000 years and has received diverse immigration flows in the last 2000 years. The original inhabitants have interacted with the pastoralist migrants from Eastern Africa (∼2000 ybp), followed by the southern Bantu migration arriving some 1000 ybp, and more recently with the European and Asian settlers after the 17th century. Many of the original Khoekhoe and San inhabitants have either become extinct or have disappeared through admixture in South Africa (SA), in a sex-biased manner involving KhoeSan women.

METHODS

In this study, we generated mitochondrial DNA (mtDNA) control region (CR) sequences for 247 South African individuals. The sampling effort was concentrated in regions and populations with historical links to the KhoeSan population groups: admixed (Coloured, Griqua), Nama (Khoekhoe) and Bantu in three provinces. Here we evaluate the composition and extent of connectivity between population groups and regions, and to assess the distribution of haplotypes for the practical application of mtDNA CR data in forensic identifications.

RESULTS

The analysis of the newly generated sequences revealed 142 distinct haplotypes, of which 122 were unique. Haplogroup L0 was predominant (overall 71.7%). A high-frequency L0d2a haplotype dominated the pool of the admixed groups with 10%-12.5% incidence overall or per region. Comparative analysis with 545 extant mtDNA CR sequences from South African KhoeSan and admixed descendants revealed extensive population structure and high within-group haplotype sharing.

CONCLUSION

The observed population and regional variations, combined with the prevalence of high-frequency haplotypes, align with patterns of matrilocality. These findings highlight the limitations of using mtDNA control region analysis for forensic applications in South Africa.

A sedimentary ancient DNA perspective on human and carnivore persistence through the Late Pleistocene in El Mirón Cave, Spain

Caves are primary sites for studying human and animal subsistence patterns and genetic ancestry throughout the Palaeolithic. Iberia served as a critical human and animal refugium in Europe during the Last Glacial Maximum (LGM), 26.5 to 19 thousand years before the present (cal kya). Therefore, it is a key location for understanding human and animal population dynamics during this event. We recover and analyse sedimentary ancient DNA (sedaDNA) data from the lower archaeological stratigraphic sequence of El Mirón Cave (Cantabria, Spain), encompassing the (1) Late Mousterian period, associated with Neanderthals, and (2) the Gravettian (c. 31.5 cal kya), Solutrean (c. 24.5-22 cal kya), and Initial Magdalenian (d. 21-20.5 cal kya) periods, associated with anatomically modern humans. We identify 28 animal taxa including humans. Fifteen of these taxa had not been identified from the archaeozoological (i.e., faunal) record, including the presence of hyenas in the Magdalenian. Additionally, we provide phylogenetic analyses on 70 sedaDNA mtDNA genomes of fauna including the densest Iberian Pleistocene sampling of C. lupus. Finally, we recover three human mtDNA sequences from the Solutrean levels. These sequences, along with published data, suggest mtDNA haplogroup continuity in Iberia throughout the Solutrean/Last Glacial Maximum period.

Biomolecular analysis of the Epigravettian human remains from Riparo Tagliente in northern Italy

The Epigravettian human remains from Riparo Tagliente in northern Italy represent some of the earliest evidence of human occupation in the southern Alpine slopes after the Last Glacial Maximum. Genomic analyses of the 17,000-year-old Tagliente 2 mandible revealed the oldest presence of a genetic profile with affinities to the Near East in the Italian peninsula, which later became the most widespread hunter-gatherer ancestry across Europe. However, a comparable biomolecular characterization of the Tagliente 1 burial remains unavailable, preventing us from defining its biological relationships with Tagliente 2. Here, we apply paleogenomic, isotopic, and radiocarbon dating analyses on a femur fragment of Tagliente 1 and compare the reconstructed data with previously reported results from Tagliente 2. Despite their different isotopic signatures and non-overlapping radiocarbon dates, we reveal that the two human remains belong to the same male individual. We determine that the distinct isotopic values can be explained by different dietary practices during lifetime, whereas the non-overlapping radiocarbon dates can be caused by minimal radiocarbon contamination, possibly deriving from chemical treatments for conservation purposes. These findings highlight the importance of interdisciplinary biomolecular studies in offering new perspectives on the Palaeolithic fossil record and addressing long-standing bioarchaeological questions.

Life history and ancestry of the late Upper Palaeolithic infant from Grotta delle Mura, Italy

The biological aspects of infancy within late Upper Palaeolithic populations and the role of southern refugia at the end of the Last Glacial Maximum are not yet fully understood. This study presents a multidisciplinary, high temporal resolution investigation of an Upper Palaeolithic infant from Grotta delle Mura (Apulia, southern Italy) combining palaeogenomics, dental palaeohistology, spatially-resolved geochemical analyses, direct radiocarbon dating, and traditional anthropological studies. The skeletal remains of the infant - Le Mura 1 - were directly dated to 17,320-16,910 cal BP. The results portray a biological history of the infant's development, early life, health and death (estimated at ~72 weeks). They identify, several phenotypic traits and a potential congenital disease in the infant, the mother's low mobility during gestation, and a high level of endogamy. Furthermore, the genomic data indicates an early spread of the Villabruna-like components along the Italian peninsula, confirming a population turnover around the time of the Last Glacial Maximum, and highlighting a general reduction in genetic variability from northern to southern Italy. Overall, Le Mura 1 contributes to our better understanding of the early stages of life and the genetic puzzle in the Italian peninsula at the end of the Last Glacial Maximum.

Long genetic and social isolation in Neanderthals before their extinction

Neanderthal genomes have been recovered from sites across Eurasia, painting an increasingly complex picture of their populations' structure that mostly indicates that late European Neanderthals belonged to a single metapopulation with no significant evidence of population structure. Here, we report the discovery of a late Neanderthal individual, nicknamed "Thorin," from Grotte Mandrin in Mediterranean France, and his genome. These dentognathic fossils, including a rare example of distomolars, are associated with a rich archeological record of Neanderthal final technological traditions in this region ∼50-42 thousand years ago. Thorin's genome reveals a relatively early divergence of ∼105 ka with other late Neanderthals. Thorin belonged to a population with a small group size that showed no genetic introgression with other known late European Neanderthals, revealing some 50 ka of genetic isolation of his lineage despite them living in neighboring regions. These results have important implications for resolving competing hypotheses about causes of the disappearance of the Neanderthals.

Human population dynamics in Upper Paleolithic Europe inferred from fossil dental phenotypes

Despite extensive archaeological research, our knowledge of the human population history of Upper Paleolithic Europe remains limited, primarily due to the scarce availability and poor molecular preservation of fossil remains. As teeth dominate the fossil record and preserve genetic signatures in their morphology, we compiled a large dataset of 450 dentitions dating between ~47 and 7 thousand years ago (ka), outnumbering existing skeletal and paleogenetic datasets. We tested a range of competing demographic scenarios using a coalescent-based machine learning Approximate Bayesian Computation (ABC) framework that we modified for use with phenotypic data. Mostly in agreement with but also challenging some of the hitherto available evidence, we identified a population turnover in western Europe at ~28 ka, isolates in western and eastern refugia between ~28 and 14.7 ka, and bottlenecks during the Last Glacial Maximum. Methodologically, this study marks the pioneering application of ABC to skeletal phenotypes, paving the way for exciting future research avenues.

Ancient DNA sheds light on the funerary practices of late Neolithic collective burial in southern France

The Aven de la Boucle (Corconne, Gard, southern France) is a karst shaft used as a collective burial between 3600 and 2800 cal BCE. The site encompasses the skeletal remains of approximately 75 individuals comprising a large majority of adult individuals, represented by scattered and commingled remains. To date, few studies have explored the potential of ancient DNA to tackle the documentation of Neolithic collective burials, and the funerary selection rules within such structures remain largely debated. In this study, we combine genomic analysis of 37 individuals with archaeo-anthropological data and Bayesian modelling of radiocarbon dates. Through this multidisciplinary approach, we aim to characterize the identity of the deceased and their relationships, as well as untangle the genetic diversity and funerary dynamics of this community. Genomic results identify 76% of male Neolithic individuals, suggesting a marked sex-biased selection. Available data emphasize the importance of biological relatedness and a male-mediated transmission of social status, as the affiliation to a specific male-lineage appears as a preponderant selection factor. The genomic results argue in favour of 'continuous' deposits between 3600 and 2800 BCE, carried out by the same community, despite cultural changes reflected by the ceramic material.

Polygenic Indices (a.k.a. Polygenic Scores) in Social Science: A Guide for Interpretation and Evaluation

Polygenic indices (PGI)-the new recommended label for polygenic scores (PGS) in social science-are genetic summary scales often used to represent an individual's liability for a disease, trait, or behavior based on the additive effects of measured genetic variants. Enthusiasm for linking genetic data with social outcomes and the inclusion of premade PGIs in social science datasets have facilitated increased uptake of PGIs in social science research-a trend that will likely continue. Yet, most social scientists lack the expertise to interpret and evaluate PGIs in social science research. Here, we provide a primer on PGIs for social scientists focusing on key concepts, unique statistical genetic considerations, and best practices in calculation, estimation, reporting, and interpretation. We summarize our recommended best practices as a checklist to aid social scientists in evaluating and interpreting studies with PGIs. We conclude by discussing the similarities between PGIs and standard social science scales and unique interpretative considerations.

Convergent Mutations and Single Nucleotide Variants in Mitochondrial Genomes of Modern Humans and Neanderthals

The genetic contributions of Neanderthals to the modern human genome have been evidenced by the comparison of present-day human genomes with paleogenomes. Neanderthal signatures in extant human genomes are attributed to intercrosses between Neanderthals and archaic anatomically modern humans (AMHs). Although Neanderthal signatures are well documented in the nuclear genome, it has been proposed that there is no contribution of Neanderthal mitochondrial DNA to contemporary human genomes. Here we show that modern human mitochondrial genomes contain 66 potential Neanderthal signatures, or Neanderthal single nucleotide variants (N-SNVs), of which 36 lie in coding regions and 7 result in nonsynonymous changes. Seven N-SNVs are associated with traits such as cycling vomiting syndrome, Alzheimer's disease and Parkinson's disease, and two N-SNVs are associated with intelligence quotient. Based on recombination tests, principal component analysis (PCA) and the complete absence of these N-SNVs in 41 archaic AMH mitogenomes, we conclude that convergent evolution, and not recombination, explains the presence of N-SNVs in present-day human mitogenomes.

The Persian plateau served as hub for Homo sapiens after the main out of Africa dispersal

A combination of evidence, based on genetic, fossil and archaeological findings, indicates that Homo sapiens spread out of Africa between ~70-60 thousand years ago (kya). However, it appears that once outside of Africa, human populations did not expand across all of Eurasia until ~45 kya. The geographic whereabouts of these early settlers in the timeframe between ~70-60 to 45 kya has been difficult to reconcile. Here we combine genetic evidence and palaeoecological models to infer the geographic location that acted as the Hub for our species during the early phases of colonisation of Eurasia. Leveraging on available genomic evidence we show that populations from the Persian Plateau carry an ancestry component that closely matches the population that settled the Hub outside Africa. With the paleoclimatic data available to date, we built ecological models showing that the Persian Plateau was suitable for human occupation and that it could sustain a larger population compared to other West Asian regions, strengthening this claim.

Population genomics of post-glacial western Eurasia

Western Eurasia witnessed several large-scale human migrations during the Holocene1-5. Here, to investigate the cross-continental effects of these migrations, we shotgun-sequenced 317 genomes-mainly from the Mesolithic and Neolithic periods-from across northern and western Eurasia. These were imputed alongside published data to obtain diploid genotypes from more than 1,600 ancient humans. Our analyses revealed a 'great divide' genomic boundary extending from the Black Sea to the Baltic. Mesolithic hunter-gatherers were highly genetically differentiated east and west of this zone, and the effect of the neolithization was equally disparate. Large-scale ancestry shifts occurred in the west as farming was introduced, including near-total replacement of hunter-gatherers in many areas, whereas no substantial ancestry shifts happened east of the zone during the same period. Similarly, relatedness decreased in the west from the Neolithic transition onwards, whereas, east of the Urals, relatedness remained high until around 4,000 BP, consistent with the persistence of localized groups of hunter-gatherers. The boundary dissolved when Yamnaya-related ancestry spread across western Eurasia around 5,000 BP, resulting in a second major turnover that reached most parts of Europe within a 1,000-year span. The genetic origin and fate of the Yamnaya have remained elusive, but we show that hunter-gatherers from the Middle Don region contributed ancestry to them. Yamnaya groups later admixed with individuals associated with the Globular Amphora culture before expanding into Europe. Similar turnovers occurred in western Siberia, where we report new genomic data from a 'Neolithic steppe' cline spanning the Siberian forest steppe to Lake Baikal. These prehistoric migrations had profound and lasting effects on the genetic diversity of Eurasian populations.

Upper Paleolithic tuberculosis: A probable case illustrated by paleoimaging methods (Azilian site of Les Iboussières, France)

The aim of this paper is to present the results of μCT-scan and 3D imaging analyses of two skeletal lesions observed on human remains of one of the last European hunter-gatherers from the late Paleolithic (Azilian period): a sacroiliac osteoarthritis and a femoral lesion suggestive of a soft tissue abscess imprint. These two skeletal elements (fused left sacrum and coxal bone, and right femur) displayed osteometric criteria indicating that they belonged to the same individual. These two associated lesions are consistent with a low-grade osteoarticular infection, and suggest a diagnosis of pelvic tuberculosis with a cold abscess of the thigh. If molecular confirmation is obtained, this case would be the first evidence of tuberculosis among Upper Palaeolithic populations.

The paleopathology and paleoepidemiology of Upper paleolithic tuberculosis: Review of evidence and hypotheses

Molecular phylogeny work has shown that tuberculosis is ancient human-adapted infection predating the Neolithic period. They also show that the Upper Paleolithic is a key period of emergence of the MTB complex strains, contemporary with the exit of modern man from Africa. Despite the richness of Upper Paleolithic sites in Eurasia and the relative abundance of human remains, the only proven case of Paleolithic tuberculosis has been described so far date from the Azilian, a culture of the European Final Paleolithic, which is more recent than the ancient Neolithic sites of the Near East, area that currently hold the record for the oldest paleopathological evidence of tuberculosis. The purpose of this review is to present evidence for the existence of tuberculosis in the Paleolithic and to list hypotheses explaining the weak demonstrative contribution of paleopathology for pre-Neolithic periods.

Unraveling the mitochondrial phylogenetic landscape of Thailand reveals complex admixture and demographic dynamics

The evolutionary dynamics of mitochondrial DNA within the Thai population were comprehensively explored with a specific focus on the influence of South Asian admixture. A total of 166 samples were collected through randomized sampling, ensuring a diverse representation. Our findings unveil substantial genetic and haplogroup diversity within the Thai population. We have identified 164 haplotypes categorized into 97 haplogroups, with a notable inclusion of 20 novel haplogroups. The distribution of haplogroups exhibited variations across different populations and countries. The central Thai population displayed a high diversity of haplogroups from both the M and N clades. Maternal lineage affinities were discerned between several Mainland Southeast Asia (MSEA) and South Asian populations, implying ancestral genetic connections and a substantial influence of South Asian women in establishing these relationships. f4-statistics indicates the presence of a Tibeto-Burman genetic component within the Mon population from Thailand. New findings demonstrate two phases of population expansion occurring 22,000-26,000 and 2500-3800 years ago, coinciding with the Last Glacial Maximum, and Neolithic demographic transition, respectively. This research significantly enhances our understanding of the maternal genetic history of Thailand and MSEA, emphasizing the influence of South Asian admixture. Moreover, it underscores the critical role of prior information, such as mutation rates, within the Bayesian framework for accurate estimation of coalescence times and inferring demographic history.

Deep-time paleogenomics and the limits of DNA survival

Although most ancient DNA studies have focused on the last 50,000 years, paleogenomic approaches can now reach into the early Pleistocene, an epoch of repeated environmental changes that shaped present-day biodiversity. Emerging deep-time genomic transects, including from DNA preserved in sediments, will enable inference of adaptive evolution, discovery of unrecognized species, and exploration of how glaciations, volcanism, and paleomagnetic reversals shaped demography and community composition. In this Review, we explore the state-of-the-art in paleogenomics and discuss key challenges, including technical limitations, evolutionary divergence and associated biases, and the need for more precise dating of remains and sediments. We conclude that with improvements in laboratory and computational methods, the emerging field of deep-time paleogenomics will expand the range of questions addressable using ancient DNA.

A pan-European dataset revealing variability in lithic technology, toolkits, and artefact shapes ~15-11 kya

Comparative macro-archaeological investigations of the human deep past rely on the availability of unified, quality-checked datasets integrating different layers of observation. Information on the durable and ubiquitous record of Paleolithic stone artefacts and technological choices are especially pertinent to this endeavour. We here present a large expert-sourced collaborative dataset for the study of stone tool technology and artefact shape evolution across Europe between ~15.000 and 11.000 years before present. The dataset contains a compendium of key sites from the study period, and data on lithic technology and toolkit composition at the level of the cultural taxa represented by those sites. The dataset further encompasses 2D shapes of selected lithic artefact groups (armatures, endscrapers, and borers/perforators) shared between cultural taxa. These data offer novel possibilities to explore between-regional patterns of material culture change to reveal scale-dependent processes of long-term technological evolution in mobile hunter-gatherer societies at the end of the Pleistocene. Our dataset facilitates state-of-the-art quantitative analyses and showcases the benefits of collaborative data collation and synthesis.

Early presence of Homo sapiens in Southeast Asia by 86-68 kyr at Tam Pà Ling, Northern Laos

The timing of the first arrival of Homo sapiens in East Asia from Africa and the degree to which they interbred with or replaced local archaic populations is controversial. Previous discoveries from Tam Pà Ling cave (Laos) identified H. sapiens in Southeast Asia by at least 46 kyr. We report on a recently discovered frontal bone (TPL 6) and tibial fragment (TPL 7) found in the deepest layers of TPL. Bayesian modeling of luminescence dating of sediments and U-series and combined U-series-ESR dating of mammalian teeth reveals a depositional sequence spanning ~86 kyr. TPL 6 confirms the presence of H. sapiens by 70 ± 3 kyr, and TPL 7 extends this range to 77 ± 9 kyr, supporting an early dispersal of H. sapiens into Southeast Asia. Geometric morphometric analyses of TPL 6 suggest descent from a gracile immigrant population rather than evolution from or admixture with local archaic populations.

A previously uncharacterized Factor Associated with Metabolism and Energy (FAME/C14orf105/CCDC198/1700011H14Rik) is related to evolutionary adaptation, energy balance, and kidney physiology

In this study we use comparative genomics to uncover a gene with uncharacterized function (1700011H14Rik/C14orf105/CCDC198), which we hereby name FAME (Factor Associated with Metabolism and Energy). We observe that FAME shows an unusually high evolutionary divergence in birds and mammals. Through the comparison of single nucleotide polymorphisms, we identify gene flow of FAME from Neandertals into modern humans. We conduct knockout experiments on animals and observe altered body weight and decreased energy expenditure in Fame knockout animals, corresponding to genome-wide association studies linking FAME with higher body mass index in humans. Gene expression and subcellular localization analyses reveal that FAME is a membrane-bound protein enriched in the kidneys. Although the gene knockout results in structurally normal kidneys, we detect higher albumin in urine and lowered ferritin in the blood. Through experimental validation, we confirm interactions between FAME and ferritin and show co-localization in vesicular and plasma membranes.

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

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

Palaeogenomics of Upper Palaeolithic to Neolithic European hunter-gatherers

Modern humans have populated Europe for more than 45,000 years1,2. Our knowledge of the genetic relatedness and structure of ancient hunter-gatherers is however limited, owing to the scarceness and poor molecular preservation of human remains from that period3. Here we analyse 356 ancient hunter-gatherer genomes, including new genomic data for 116 individuals from 14 countries in western and central Eurasia, spanning between 35,000 and 5,000 years ago. We identify a genetic ancestry profile in individuals associated with Upper Palaeolithic Gravettian assemblages from western Europe that is distinct from contemporaneous groups related to this archaeological culture in central and southern Europe4, but resembles that of preceding individuals associated with the Aurignacian culture. This ancestry profile survived during the Last Glacial Maximum (25,000 to 19,000 years ago) in human populations from southwestern Europe associated with the Solutrean culture, and with the following Magdalenian culture that re-expanded northeastward after the Last Glacial Maximum. Conversely, we reveal a genetic turnover in southern Europe suggesting a local replacement of human groups around the time of the Last Glacial Maximum, accompanied by a north-to-south dispersal of populations associated with the Epigravettian culture. From at least 14,000 years ago, an ancestry related to this culture spread from the south across the rest of Europe, largely replacing the Magdalenian-associated gene pool. After a period of limited admixture that spanned the beginning of the Mesolithic, we find genetic interactions between western and eastern European hunter-gatherers, who were also characterized by marked differences in phenotypically relevant variants.

Western visitors at the Blätterhöhle (city of Hagen, southern Westphalia) during the Younger Dryas? A new final palaeolithic assemblage type in western Germany

Until now, it was considered certain that the last reindeer hunters of the Ahrensburgian (tanged point groups) existed exclusively in northwestern Central Europe during the Younger Dryas Cold Period (~ Greenland Stadial 1). The excavations carried out since 2006 on the forecourt (Vorplatz) of the small Blätterhöhle in Hagen on the northern edge of the Sauerland uplands of southern Westphalia (North Rhine-Westphalia, western Germany) have now changed this view. Beneath a surprisingly extensive sequence of Mesolithic find horizons, Pleistocene sediments could be reached whose excavations yielded a Final Palaeolithic lithic ensemble of the Younger Dryas, unusual for the region and beyond. It is characterised by numerous backed lithic projectile points of high variability. Comparisons suggest a typological-technological connection with the Western European Laborian / Late Laborian. Neither in the nearer nor in the wider surroundings has a comparable lithic find ensemble been found so far. In addition, there is a lack of clear evidence for the reindeer in the fauna. Surprisingly, the vast majority of radiocarbon dates of bones and charcoals from the investigated archaeological horizon of the Final Pleistocene proved to be significantly older than expected from their stratigraphic position. This phenomenon has not yet been clarified.

Post-LGM environments and foragers on the move: New data from the lower Altmühl Valley (Franconian Jura, SE Germany)

After the Last Glacial Maximum, the Swabian and Franconian Jura (in SW and SE Germany, respectively) were repopulated by Magdalenian hunter-gatherers within the same communication network. However, while the Magdalenian settlement of the Swabian Jura dates to 17-14 ka cal BP, permanent Magdalenian occupations in the Franconian Jura date to 15-14 ka cal BP. In comparison with its western counterpart, the Franconian Jura was mostly excavated in the early days of archaeological research. Does this different chronology reflect the different history of research? Why did Magdalenian foragers establish permanent occupation in the Franconian Jura nearly 2 millennia after settling in Swabia, despite the fact these regions are only 150 km apart? To address these questions, we reinvestigated two sites in the Altmühl Valley with micromorphology and luminescence dating, namely Felsenhäusl-Kellerhöhle and Klausennische. Our data show that both sites have intact Pleistocene deposits. Among these, we identified sediments dating between 17 and 15 ka that show only rare lithic artifacts and microfeatures indicative of cold and arid conditions. Our work and published data suggest that the steady settlement of Magdalenian foragers in the Altmühl Valley starting 15 ka cal BP coincides with the end of this harsh period and the onset of cool and wetter environments. Data from the Swabian Jura demonstrated that in the Lone Valley, similar environments and Magdalenian occupations commenced earlier, starting 17 ka cal BP. Therefore, we propose that regional environments acted as a barrier against the dispersal of foragers in the Franconian Jura and determined its later Magdalenian occupation. Our research highlighted that different environments, taphonomic processes, and site uses probably coexisted across the German Jura. Therefore, it remains fundamental to expand the multisite data set proposed in this article to further test hypotheses about human/environment interaction in this region.

Ancient oral microbiomes support gradual Neolithic dietary shifts towards agriculture

The human microbiome has recently become a valuable source of information about host life and health. To date little is known about how it may have evolved during key phases along our history, such as the Neolithic transition towards agriculture. Here, we shed light on the evolution experienced by the oral microbiome during this transition, comparing Palaeolithic hunter-gatherers with Neolithic and Copper Age farmers that populated a same restricted area in Italy. We integrate the analysis of 76 dental calculus oral microbiomes with the dietary information derived from the identification of embedded plant remains. We detect a stronger deviation from the hunter-gatherer microbiome composition in the last part of the Neolithic, while to a lesser extent in the early phases of the transition. Our findings demonstrate that the introduction of agriculture affected host microbiome, supporting the hypothesis of a gradual transition within the investigated populations.

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.

Following the niche: the differential impact of the last glacial maximum on four European ungulates

Predicting the effects of future global changes on species requires a better understanding of the ecological niche dynamics in response to climate; the large climatic fluctuations of the last 50,000 years can be used as a natural experiment to that aim. Here we test whether the realized niche of horse, aurochs, red deer, and wild boar changed between 47,000 and 7500 years ago using paleoecological modelling over an extensive archaeological database. We show that they all changed their niche, with species-specific responses to climate fluctuations. We also suggest that they survived the climatic turnovers thanks to their flexibility and by expanding their niche in response to the extinction of competitors and predators. Irrespective of the mechanism behind such processes, the fact that species with long generation times can change their niche over thousands of years cautions against assuming it to stay constant both when reconstructing the past and predicting the future.

European megafaunal ungulates living in open habitats over the last 50,000 years showed evidence for niche change, possibly driven by climatic change and extinction of competitors and predators

Phylogeographic analysis reveals an ancient East African origin of human herpes simplex virus 2 dispersal out-of-Africa

Human herpes simplex virus 2 (HSV-2) is a ubiquitous, slowly evolving DNA virus. HSV-2 has two primary lineages, one found in West and Central Africa and the other found worldwide. Competing hypotheses have been proposed to explain how HSV-2 migrated out-of-Africa (i)HSV-2 followed human migration out-of-Africa 50-100 thousand years ago, or (ii)HSV-2 migrated via the trans-Atlantic slave trade 150-500 years ago. Limited geographic sampling and lack of molecular clock signal has precluded robust comparison. Here, we analyze newly sequenced HSV-2 genomes from Africa to resolve geography and timing of divergence events within HSV-2. Phylogeographic analysis consistently places the ancestor of worldwide dispersal in East Africa, though molecular clock is too slow to be detected using available data. Rates 4.2 × 10-8-5.6 × 10-8 substitutions/site/year, consistent with previous age estimates, suggest a worldwide dispersal 22-29 thousand years ago. Thus, HSV-2 likely migrated with humans from East Africa and dispersed after the Last Glacial Maximum.

Aurignacian dynamics in Southeastern Europe based on spatial analysis, sediment geochemistry, raw materials, lithic analysis, and use-wear from Românești-Dumbrăvița

The Aurignacian is one of the first cultural-technological traditions commonly associated with the expansion of Homo sapiens in Europe. Early Homo sapiens demographics across the continent are therefore typically inferred using the distribution of Aurignacian assemblages. Western Romania has been used as a tie-point to connect the well-researched lithic assemblages from the eastern Mediterranean and Western Europe through its early Homo sapiens fossils. However, Romania’s archeological record remains underexplored thereby hindering our ability to directly connect better understood regions through time and space. Here we report on excavations from the open-air Middle/Upper Paleolithic site of Românești-Dumbrăvița I in southwestern Romania. Three stratified Paleolithic assemblages were extensively excavated within a 1-m-thick eolian-deposited sequence. Spatial, geochemical, raw material, techno-typological, and use-wear analysis of the site reveal patterns of artifact configuration, resource exploitation, fire history, knapping objectives, and functionality. Taken together, Românești-Dumbrăvița I is the first well-contextualized archeological site in close spatiotemporal proximity to many early, well-preserved human fossils and in East-Central Europe.

Sulfonylureas or biguanides is associated with a lower risk of rheumatoid arthritis in patients with diabetes: A nationwide cohort study

Objective

Diabetes mellitus (DM) is associated with immune dysregulation, while sulfonylureas or biguanides have been linked to anti-inflammatory mechanisms. In this study, we aimed to examine the occurrence rate of rheumatoid arthritis (RA) among DM patients and its incidence rate between different treatments.

Methods

This cohort study used the Taiwan National Health Insurance Research Database between 1997 and 2013 to evaluate the primary outcomes of the preventive role of sulfonylureas or biguanides in the development of RA. We used the Chi-square test for categorical variables and Cox proportional hazard regression and log-rank test to explore the time for development of RA in DM patients. Logistic regression was adopted to estimate the odds ratio of RA in different dosages of medication exposure.

Results

Our cohort study included 94,141 DM cases. The risk of RA development of non-sulfonylureas/biguanides users among the DM group in each analysis was set as the reference, and the adjusted hazard ratio of RA in DM patients who were using sulfonylureas or biguanides was 0.73 (95% confidence interval 0.60–0.90). Within 1 year before the index date, compared with no-biguanides users, patients with more than 180 days of prescription of biguanides had a significantly lower RA risk. Similarly, the significantly lower risk of RA was still observed in DM patients who had more than 365 days of prescription of sulfonylurea within 2 or 3 years before the index date of first RA visit (all p &lt; 0.05).

Conclusion

Our data suggest that sulfonylureas or biguanides are associated with a lower rate of RA development in patients with DM; the effect of biguanides appeared more rapid than that of sulfonylureas, but the sulfonylureas might have a longer effect on lowering RA development incidence.

Variation in the substitution rates among the human mitochondrial haplogroup U sublineages

Resolving the absolute timescale of phylogenetic trees stipulates reliable estimates for the rate of DNA sequence evolution. For this end, various calibration methods have been developed and studied intensively. Intraspecific rate variation among distinct genetic lineages, however, has gained less attention. Here, we have assessed lineage-specific molecular rates of human mitochondrial DNA (mtDNA) by performing tip-calibrated Bayesian phylogenetic analyses. Tip-calibration, as opposed to traditional nodal time stamps from dated fossil evidence or geological events, is based on sample ages and becoming ever more feasible as ancient DNA data from radiocarbon-dated samples accumulate. We focus on subhaplogroups U2, U4, U5a, and U5b, the data including ancient mtDNA genomes from ¹⁴C-dated samples (n = 234), contemporary genomes (n = 301), and two outgroup sequences from haplogroup R. The obtained molecular rates depended on the data sets (with or without contemporary sequences), suggesting time-dependency. More notable was the rate variation between haplogroups: U4 and U5a stand out having a substantially higher rate than U5b. This is also reflected in the divergence times obtained (U5a: 17,700 years and U5b: 29,700 years), a disparity not reported previously. After ruling out various alternative causes (e.g., selection, sampling, and sequence quality), we propose that the substitution rates have been influenced by demographic histories, widely different among populations where U4/U5a or U5b are frequent. As with the Y-chromosomal subhaplogroup R1b, the mitochondrial U4 and U5a have been associated with remarkable range extensions of the Yamnaya culture in the Bronze Age.

Tracing the mobility of a Late Epigravettian (~ 13 ka) male infant from Grotte di Pradis (Northeastern Italian Prealps) at high-temporal resolution

We present the results of a multi-disciplinary investigation on a deciduous human tooth (Pradis 1), recently recovered from the Epigravettian layers of the Grotte di Pradis archaeological site (Northeastern Italian Prealps). Pradis 1 is an exfoliated deciduous molar (Rdm₂), lost during life by an 11–12-year-old child. A direct radiocarbon date provided an age of 13,088–12,897 cal BP (95% probability, IntCal20). Amelogenin peptides extracted from tooth enamel and analysed through LC–MS/MS indicate that Pradis 1 likely belonged to a male. Time-resolved ⁸⁷Sr/⁸⁶Sr analyses by laser ablation mass spectrometry (LA-MC-ICPMS), combined with dental histology, were able to resolve his movements during the first year of life (i.e. the enamel mineralization interval). Specifically, the Sr isotope ratio of the tooth enamel differs from the local baseline value, suggesting that the child likely spent his first year of life far from Grotte di Pradis. Sr isotopes are also suggestive of a cyclical/seasonal mobility pattern exploited by the Epigravettian human group. The exploitation of Grotte di Pradis on a seasonal, i.e. summer, basis is also indicated by the faunal spectra. Indeed, the nearly 100% occurrence of marmot remains in the entire archaeozoological collection indicates the use of Pradis as a specialized marmot hunting or butchering site. This work represents the first direct assessment of sub-annual movements observed in an Epigravettian hunter-gatherer group from Northern Italy.

Evolution and dispersal of mitochondrial DNA haplogroup U5 in Northern Europe: insights from an unsupervised learning approach to phylogeography

Background

We combined an unsupervised learning methodology for analyzing mitogenome sequences with maximum likelihood (ML) phylogenetics to make detailed inferences about the evolution and diversification of mitochondrial DNA (mtDNA) haplogroup U5, which appears at high frequencies in northern Europe.

Methods

Haplogroup U5 mitogenome sequences were gathered from GenBank. The hierarchal Bayesian Analysis of Population Structure (hierBAPS) method was used to generate groups of sequences that were then projected onto a rooted maximum likelihood (ML) phylogenetic tree to visualize the pattern of clustering. The haplogroup statuses of the individual sequences were assessed using Haplogrep2.

Results

A total of 23 hierBAPS groups were identified, all of which corresponded to subclades defined in Phylotree, v.17. The hierBAPS groups projected onto the ML phylogeny accurately clustered all haplotypes belonging to a specific haplogroup in accordance with Haplogrep2. By incorporating the geographic source of each sequence and subclade age estimates into this framework, inferences about the diversification of U5 mtDNAs were made. Haplogroup U5 has been present in northern Europe since the Mesolithic, and spread in both eastern and western directions, undergoing significant diversification within Scandinavia. A review of historical and archeological evidence attests to some of the population interactions contributing to this pattern.

Conclusions

The hierBAPS algorithm accurately grouped mitogenome sequences into subclades in a phylogenetically robust manner. This analysis provided new insights into the phylogeographic structure of haplogroup U5 diversity in northern Europe, revealing a detailed perspective on the diversity of subclades in this region and their distribution in Scandinavian populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08572-y.

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

Mitochondrial DNA Consensus Calling and Quality Filtering for Constructing Ancient Human Mitogenomes: Comparison of Two Widely Applied Methods

Retrieving high-quality endogenous ancient DNA (aDNA) poses several challenges, including low molecular copy number, high rates of fragmentation, damage at read termini, and potential presence of exogenous contaminant DNA. All these factors complicate a reliable reconstruction of consensus aDNA sequences in reads from high-throughput sequencing platforms. Here, we report findings from a thorough evaluation of two alternative tools (ANGSD and schmutzi) aimed at overcoming these issues and constructing high-quality ancient mitogenomes. Raw genomic data (BAM/FASTQ) from a total of 17 previously published whole ancient human genomes ranging from the 14th to the 7th millennium BCE were retrieved and mitochondrial consensus sequences were reconstructed using different quality filters, with their accuracy measured and compared. Moreover, the influence of different sequence parameters (number of reads, sequenced bases, mean coverage, and rate of deamination and contamination) as predictors of derived sequence quality was evaluated. Complete mitogenomes were successfully reconstructed for all ancient samples, and for the majority of them, filtering substantially improved mtDNA consensus calling and haplogroup prediction. Overall, the schmutzi pipeline, which estimates and takes into consideration exogenous contamination, appeared to have the edge over the much faster and user-friendly alternative method (ANGSD) in moderate to high coverage samples (>1,000,000 reads). ANGSD, however, through its read termini trimming filter, showed better capabilities in calling the consensus sequence from low-quality samples. Among all the predictors of overall sample quality examined, the strongest correlation was found for the available number of sequence reads and bases. In the process, we report a previously unassigned haplogroup (U3b) for an Early Chalcolithic individual from Southern Anatolia/Northern Levant.

Mitochondrial DNA Footprints from Western Eurasia in Modern Mongolia

Mongolia is located in a strategic position at the eastern edge of the Eurasian Steppe. Nomadic populations moved across this wide area for millennia before developing more sedentary communities, extended empires, and complex trading networks, which connected western Eurasia and eastern Asia until the late Medieval period. We provided a fine-grained portrait of the mitochondrial DNA (mtDNA) variation observed in present-day Mongolians and capable of revealing gene flows and other demographic processes that took place in Inner Asia, as well as in western Eurasia. The analyses of a novel dataset (N = 2,420) of mtDNAs highlighted a clear matrilineal differentiation within the country due to a mixture of haplotypes with eastern Asian (EAs) and western Eurasian (WEu) origins, which were differentially lost and preserved. In a wider genetic context, the prevalent EAs contribution, larger in eastern and central Mongolian regions, revealed continuous connections with neighboring Asian populations until recent times, as attested by the geographically restricted haplotype-sharing likely facilitated by the Genghis Khan’s so-called Pax Mongolica. The genetic history beyond the WEu haplogroups, notably detectable on both sides of Mongolia, was more difficult to explain. For this reason, we moved to the analysis of entire mitogenomes (N = 147). Although it was not completely possible to identify specific lineages that evolved in situ, two major changes in the effective (female) population size were reconstructed. The more recent one, which began during the late Pleistocene glacial period and became steeper in the early Holocene, was probably the outcome of demographic events connected to western Eurasia. The Neolithic growth could be easily explained by the diffusion of dairy pastoralism, as already proposed, while the late glacial increase indicates, for the first time, a genetic connection with western Eurasian refuges, as supported by the unusual high frequency and internal sub-structure in Mongolia of haplogroup H1, a well-known post-glacial marker in Europe. Bronze Age events, without a significant demographic impact, might explain the age of some mtDNA haplogroups. Finally, a diachronic comparison with available ancient mtDNAs made it possible to link six mitochondrial lineages of present-day Mongolians to the timeframe and geographic path of the Silk Route.

Modern human teeth unearthed from below the ∼128,000-year-old level at Punung, Java: A case highlighting the problem of recent intrusion in cave sediments

The emergence of modern humans in the eastern edge of the Eurasian Continent is debated between two major models: early (∼130-70 ka) and late (∼50 ka) dispersal models. The former view is grounded mainly on the claims that several cave sites in Southeast Asia and southern China yielded modern human fossils of those early ages, but such reports have been disputed for the lack of direct dating of the human remains and insufficient documentation of stratigraphy and taphonomy. By tracing possible burial process and conducting direct dating for an early Late Pleistocene paleontological site of Punung III, East Java, we here report a case that demonstrates how unexpected intrusion of recent human remains into older stratigraphic levels could occur in cave sediments. This further highlights the need of direct dating and taphonomic assessment before accepting either model. We also emphasize that the state of fossilization of bones and teeth is a useful guide for initial screening of recent intrusion and should be reported particularly when direct dating is unavailable. Additionally, we provide a revised stratigraphy and faunal list of Punung III, a key site that defines the tropical rainforest Punung Fauna during the early Late Pleistocene of the region.

An infant burial from Arma Veirana in northwestern Italy provides insights into funerary practices and female personhood in early Mesolithic Europe

The evolution and development of human mortuary behaviors is of enormous cultural significance. Here we report a richly-decorated young infant burial (AVH-1) from Arma Veirana (Liguria, northwestern Italy) that is directly dated to 10,211-9910 cal BP (95.4% probability), placing it within the early Holocene and therefore attributable to the early Mesolithic, a cultural period from which well-documented burials are exceedingly rare. Virtual dental histology, proteomics, and aDNA indicate that the infant was a 40-50 days old female. Associated artifacts indicate significant material and emotional investment in the child's interment. The detailed biological profile of AVH-1 establishes the child as the earliest European near-neonate documented to be female. The Arma Veirana burial thus provides insight into sex/gender-based social status, funerary treatment, and the attribution of personhood to the youngest individuals among prehistoric hunter-gatherer groups and adds substantially to the scant data on mortuary practices from an important period in prehistory shortly following the end of the last Ice Age.

A multi-proxy approach to exploring Homo sapiens' arrival, environments and adaptations in Southeast Asia

The capability of Pleistocene hominins to successfully adapt to different types of tropical forested environments has long been debated. In order to investigate environmental changes in Southeast Asia during a critical period for the turnover of hominin species, we analysed palaeoenvironmental proxies from five late Middle to Late Pleistocene faunas. Human teeth discoveries have been reported at Duoi U'Oi, Vietnam (70-60 ka) and Nam Lot, Laos (86-72 ka). However, the use of palaeoproteomics allowed us to discard the latter, and, to date, no human remains older than ~ 70 ka are documented in the area. Our findings indicate that tropical rainforests were highly sensitive to climatic changes over that period, with significant fluctuations of the canopy forests. Locally, large-bodied faunas were resilient to these fluctuations until the cooling period of the Marine Isotope Stage 4 (MIS 4; 74-59 ka) that transformed the overall biotope. Then, under strong selective pressures, populations with new phenotypic characteristics emerged while some other species disappeared. We argue that this climate-driven shift offered new foraging opportunities for hominins in a novel rainforest environment and was most likely a key factor in the settlement and dispersal of our species during MIS 4 in SE Asia.

A benchmarking of human mitochondrial DNA haplogroup classifiers from whole-genome and whole-exome sequence data

The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. We also assessed the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short-read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

Insights into human history from the first decade of ancient human genomics

[Figure: see text].

The genetic and cultural impact of the Steppe migration into Europe

BACKGROUND

During the early 3rd millennium BCE migration from Pontic Steppe, mainly related to Yamnaya culture, has affected European populations both culturally and genetically, however, it has long been debated to what extent this migration was male-driven, and how this replacement process took place which eliminated partially/largely Neolithic male lines over time.

AIM

This paper aims to evaluate the influence of the Steppe migration on European Bronze Age populations by calculating both male and female genetic contributions of the Steppe-related ancestry to the European Bronze Age populations. With this approach, we will be able to clarify the hypotheses on whether it was male-biased migration or not.

SUBJECTS AND METHODS

To evaluate the genetic impact and the proportion of the Steppe-related ancestry to the European Bronze Age populations, we performed PCA and qpAdm analyses by using published genome-wide data. In addition, we quantified male and female genetic contribution into Europe by using the analysis of uniparental markers and the X-chromosome.

RESULTS

The Steppe migration had a considerable impact on the genetic makeup of the Bronze Age European populations. The data suggest that the Steppe-related ancestry arriving into Central Europe was male-driven, dominantly in the Corded Ware culture populations and lesser in the Bell Beaker populations. In fact, there is no evidence that this migration had a significant input on the mitochondrial genetic pool of all European Bronze Age populations.

CONCLUSIONS

Our analyses suggest that the Steppe-related ancestry had genetic impact on mainly Central-Eastern Europe. Moreover, this migration was male-driven for most of the Central European populations belonging to the Corded Ware groups, and to a lesser extent for the Bell Beaker groups.

More data on ancient human mitogenome variability in Italy: new mitochondrial genome sequences from three Upper Palaeolithic burials

BACKGROUND

Recently, the study of mitochondrial variability in ancient humans has allowed the definition of population dynamics that characterised Europe in the Late Pleistocene and Early Holocene. Despite the abundance of sites and skeletal remains few data are available for Italy.

AIM

We reconstructed the mitochondrial genomes of three Upper Palaeolithic individuals for some of the most important Italian archaeological contexts: Paglicci (South-Eastern Italy), San Teodoro (South-Western Italy) and Arene Candide (North-Western Italy) caves.

SUBJECTS AND METHODS

We explored the phylogenetic relationships of the three mitogenomes in the context of Western Eurasian ancient and modern variability.

RESULTS

Paglicci 12 belongs to sub-haplogroup U8c, described in only two other Gravettian individuals; San Teodoro 2 harbours a U2'3'4'7'8'9 sequence, the only lineage found in Sicily during the Late Pleistocene and Early Holocene; Arene Candide 16 displays an ancestral U5b1 haplotype already detected in other Late Pleistocene hunter-gatherers from Central Europe.

CONCLUSION

Regional genetic continuity is highlighted in the Gravettian groups that succeeded in Paglicci. Data from one of the oldest human remains from Sicily reinforce the hypothesis that Epigravettian groups carrying U2'3'4'7'8'9 could be the first inhabitants of the island. The first pre-Neolithic mitogenome from North-Western Italy, sequenced here, shows more affinity with continental Europe than with the Italian peninsula.

Through 40,000 years of human presence in Southern Europe: the Italian case study

The Italian Peninsula, a natural pier across the Mediterranean Sea, witnessed intricate population events since the very beginning of the human occupation in Europe. In the last few years, an increasing number of modern and ancient genomes from the area have been published by the international research community. This genomic perspective started unveiling the relevance of Italy to understand the post-Last Glacial Maximum (LGM) re-peopling of Europe, the earlier phase of the Neolithic westward migrations, and its linking role between Eastern and Western Mediterranean areas after the Iron Age. However, many open questions are still waiting for more data to be addressed in full. With this review, we summarize the current knowledge emerging from the available ancient Italian individuals and, by re-analysing them all at once, we try to shed light on the avenues future research in the area should cover. In particular, open questions concern (1) the fate of pre-Villabruna Europeans and to what extent their genomic components were absorbed by the post-LGM hunter-gatherers; (2) the role of Sicily and Sardinia before LGM; (3) to what degree the documented genetic structure within the Early Neolithic settlers can be described as two separate migrations; (4) what are the population events behind the marked presence of an Iranian Neolithic-like component in Bronze Age and Iron Age Italian and Southern European samples.

Using Y-chromosome capture enrichment to resolve haplogroup H2 shows new evidence for a two-path Neolithic expansion to Western Europe

Uniparentally-inherited markers on mitochondrial DNA (mtDNA) and the non-recombining regions of the Y chromosome (NRY), have been used for the past 30 years to investigate the history of humans from a maternal and paternal perspective. Researchers have preferred mtDNA due to its abundance in the cells, and comparatively high substitution rate. Conversely, the NRY is less susceptible to back mutations and saturation, and is potentially more informative than mtDNA owing to its longer sequence length. However, due to comparatively poor NRY coverage via shotgun sequencing, and the relatively low and biased representation of Y-chromosome variants on capture assays such as the 1240 k, ancient DNA studies often fail to utilize the unique perspective that the NRY can yield. Here we introduce a new DNA enrichment assay, coined YMCA (Y-mappable capture assay), that targets the "mappable" regions of the NRY. We show that compared to low-coverage shotgun sequencing and 1240 k capture, YMCA significantly improves the mean coverage and number of sites covered on the NRY, increasing the number of Y-haplogroup informative SNPs, and allowing for the identification of previously undiscovered variants. To illustrate the power of YMCA, we show that the analysis of ancient Y-chromosome lineages can help to resolve Y-chromosomal haplogroups. As a case study, we focus on H2, a haplogroup associated with a critical event in European human history: the Neolithic transition. By disentangling the evolutionary history of this haplogroup, we further elucidate the two separate paths by which early farmers expanded from Anatolia and the Near East to western Europe.

Genome-scale sequencing and analysis of human, wolf, and bison DNA from 25,000-year-old sediment

Cave sediments have been shown to preserve ancient DNA but so far have not yielded the genome-scale information of skeletal remains. We retrieved and analyzed human and mammalian nuclear and mitochondrial environmental "shotgun" genomes from a single 25,000-year-old Upper Paleolithic sediment sample from Satsurblia cave, western Georgia:first, a human environmental genome with substantial basal Eurasian ancestry, which was an ancestral component of the majority of post-Ice Age people in the Near East, North Africa, and parts of Europe; second, a wolf environmental genome that is basal to extant Eurasian wolves and dogs and represents a previously unknown, likely extinct, Caucasian lineage; and third, a European bison environmental genome that is basal to present-day populations, suggesting that population structure has been substantially reshaped since the Last Glacial Maximum. Our results provide new insights into the Late Pleistocene genetic histories of these three species and demonstrate that direct shotgun sequencing of sediment DNA, without target enrichment methods, can yield genome-wide data informative of ancestry and phylogenetic relationships.

The evolution of human skin pigmentation involved the interactions of genetic, environmental, and cultural variables

The primary biological role of human skin pigmentation is as a mediator of penetration of ultraviolet radiation (UVR) into the deep layers of skin and the cutaneous circulation. Since the origin of Homo sapiens, dark, protective constitutive pigmentation and strong tanning abilities have been favored under conditions of high UVR and represent the baseline condition for modern humans. The evolution of partly depigmented skin and variable tanning abilities has occurred multiple times in prehistory, as populations have dispersed into environments with lower and more seasonal UVR regimes, with unique complements of genes and cultural practices. The evolution of extremes of dark pigmentation and depigmentation has been rare and occurred only under conditions of extremely high or low environmental UVR, promoted by positive selection on variant pigmentation genes followed by limited gene flow. Over time, the evolution of human skin pigmentation has been influenced by the nature and course of human dispersals and modifications of cultural practices, which have modified the nature and actions of skin pigmentation genes. Throughout most of prehistory and history, the evolution of human skin pigmentation has been a contingent and non-deterministic process.