Cervidae antlers exploited to manufacture prehistoric tools and hunting implements as a reliable source of ancient DNA.

Antler is one of the primary animal raw materials exploited for technical purposes by the hunter-gatherer groups of the Eurasian Upper Palaeolithic (UP) all over the ecological range of deers, and beyond. It was exhaustively employed to produce one of the most critical tools for the survival of the UP societies: hunting weapons. However, antler implements can be made from diverse deer taxa, with different ecological requirements and ethological behaviours. Identifying the antler's origin at a taxonomic level is thus essential in improving our knowledge of humans' functional, practical and symbolic choices, as well as the human-animal interface during Prehistoric times. Nevertheless, palaeogenetics analyses have focused mainly on bone and teeth, with genetic studies of antler generally focused on modern deer conservation. Here we present the results of the first whole mitochondrial genome ancient DNA (aDNA) analysis by means of in-solution hybridisation capture of antlers from pre-Holocene archaeological contexts. We analysed a set of 50 Palaeolithic and Neolithic (c. 34-8ka) antler and osseous objects from South-Western Europe, Central Europe, South-Western Asia and the Caucasus. We successfully obtained aDNA, allowing us to identify the exploited taxa and demonstrate the archaeological relevance of those finds. Moreover, as most of the antlers were sampled using a minimally-invasive method, further analyses (morphometric, technical, genetic, radiometric and more) remain possible on these objects.

Link between Monkeypox Virus Genomes from Museum Specimens and 1965 Zoo Outbreak.

We used pathogen genomics to test orangutan specimens from a museum in Bonn, Germany, to identify the origin of the animals and the circumstances of their death. We found monkeypox virus genomes in the samples and determined that they represent cases from a 1965 outbreak at Rotterdam Zoo in Rotterdam, the Netherlands.

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

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

The population genomic legacy of the second plague pandemic.

Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.1 It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).2 Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17th-19th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large FST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.

Northeastern Asian and Jomon-related genetic structure in the Three Kingdoms period of Gimhae, Korea.

The genetic history of prehistoric and protohistoric Korean populations is not well understood because only a small number of ancient genomes are available. Here, we report the first paleogenomic data from the Korean Three Kingdoms period, a crucial point in the cultural and historic formation of Korea. These data comprise eight shotgun-sequenced genomes from ancient Korea (0.7×-6.1× coverage). They were derived from two archeological sites in Gimhae: the Yuha-ri shell mound and the Daesung-dong tumuli, the latter being the most important funerary complex of the Gaya confederacy. All individuals are from between the 4th and 5th century CE and are best modeled as an admixture between a northern China Bronze Age genetic source and a source of Jomon-related ancestry that shares similarities with the present-day genomes from Japan. The observed substructure and proportion of Jomon-related ancestry suggest the presence of two genetic groups within the population and diversity among the Gaya population. We could not correlate the genomic differences between these two groups with either social status or sex. All the ancient individuals' genomic profiles, including phenotypically relevant SNPs associated with hair and eye color, facial morphology, and myopia, imply strong genetic and phenotypic continuity with modern Koreans for the last 1,700 years.

Genomes from Verteba cave suggest diversity within the Trypillians in Ukraine.

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

Human-specific changes in two functional enhancers of FOXP2.

FOXP2 is a gene involved in language development and function. Neanderthals and humans share the same coding region of the gene, although the formers are thought to have exhibited less sophisticated language abilities. In this paper, we report on several human-specific changes in two functional enhancers of FOXP2. Two of these variants are located within the binding sites for the transcription factors POLR2A and SMARCC1, respectively. Interestingly, SMARCC1 is involved in brain development and vitamin D metabolism. We hypothesize that the human specific change in this position might have resulted in a different regulation pattern of FOXP2 expression in our species compared to extinct hominins, with a potential impact on our language abilities.

TKGWV2: an ancient DNA relatedness pipeline for ultra-low coverage whole genome shotgun data.

Estimation of genetically related individuals is playing an increasingly important role in the ancient DNA field. In recent years, the numbers of sequenced individuals from single sites have been increasing, reflecting a growing interest in understanding the familial and social organisation of ancient populations. Although a few different methods have been specifically developed for ancient DNA, namely to tackle issues such as low-coverage homozygous data, they require a 0.1-1× minimum average genomic coverage per analysed pair of individuals. Here we present an updated version of a method that enables estimates of 1st and 2nd-degrees of relatedness with as little as 0.026× average coverage, or around 18,000 SNPs from 1.3 million aligned reads per sample with average length of 62 bp-four times less data than 0.1× coverage at similar read lengths. By using simulated data to estimate false positive error rates, we further show that a threshold even as low as 0.012×, or around 4000 SNPs from 600,000 reads, will always show 1st-degree relationships as related. Lastly, by applying this method to published data, we are able to identify previously undocumented relationships using individuals that had been excluded from prior kinship analysis due to their very low coverage. This methodological improvement has the potential to enable relatedness estimation on ancient whole genome shotgun data during routine low-coverage screening, and therefore improve project management when decisions need to be made on which individuals are to be further sequenced.

The epidemiology of Plasmodium vivax among adults in the Democratic Republic of the Congo.

Reports of P. vivax infections among Duffy-negative hosts have accumulated throughout sub-Saharan Africa. Despite this growing body of evidence, no nationally representative epidemiological surveys of P. vivax in sub-Saharan Africa have been performed. To overcome this gap in knowledge, we screened over 17,000 adults in the Democratic Republic of the Congo (DRC) for P. vivax using samples from the 2013-2014 Demographic Health Survey. Overall, we found a 2.97% (95% CI: 2.28%, 3.65%) prevalence of P. vivax infections across the DRC. Infections were associated with few risk-factors and demonstrated a relatively flat distribution of prevalence across space with focal regions of relatively higher prevalence in the north and northeast. Mitochondrial genomes suggested that DRC P. vivax were distinct from circulating non-human ape strains and an ancestral European P. vivax strain, and instead may be part of a separate contemporary clade. Our findings suggest P. vivax is diffusely spread across the DRC at a low prevalence, which may be associated with long-term carriage of low parasitemia, frequent relapses, or a general pool of infections with limited forward propagation.

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.

Has COVID-19 halted winter-spring wildfires in the Mediterranean? Insights for wildfire science under a pandemic context.

Wildfires in the Mediterranean are strongly tied to human activities. Given their particular link with humans, which act as both initiators and suppressors, wildfire hazard is highly sensitive to socioeconomic changes and patterns. Many researchers have prompted the perils of sustaining the current management policy, the so-called 'total fire exclusion'. This policy, coupled to increasingly fire-prone weather conditions, may lead to more hazardous fires in the mid-long run. Under this framework, the irruption of the COVID-19 pandemic adds to the ongoing situation. Facing the lack of an effective treatment, the only alternative was the implementation of strict lockdown strategies. The virtual halt of the system undoubtedly affected economic and social behavior, triggering cascading effects such as the drop in winter-spring wildfire activity. In this work, we discuss the main impacts, challenges and consequences that wildfire science may experience due to the pandemic situation, and identify potential opportunities for wildfire management. We investigate the recent evolution of burned area (retrieved from the MCD64A1 v006 MODIS product) in the EU Mediterranean region (Portugal, Spain, France, Italy and Greece) to ascertain to what extent the 2020 winter-spring season was impacted by the public health response to COVID-19 (curfews and lockdowns). We accounted for weather conditions (characterized using the 6-month Standardized Precipitation Evapotranspiration Index; SPEI6) to disregard possible weather effects mediating fire activity. Our results suggest that, under similar drought-related circumstances (SPEI6 ≈ -0.7), the expected burned area in 2020 during the lockdown period in the EU (March-May) would lay somewhere within the range of 38,800 ha ± 18,379 ha. Instead, the affected area stands one order of magnitude below average (3325 ha). This stresses the need of considering the social dimension in the analysis of current and future wildfire impacts in the Mediterranean region.

Author Correction: The dental proteome of Homo antecessor.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Geospatial Modeling of Containment Probability for Escaped Wildfires in a Mediterranean Region.

Despite escalating expenditures in firefighting, extreme fire events continue to pose a major threat to ecosystem services and human communities in Mediterranean areas. Developing a safe and effective fire response is paramount to efficiently restrict fire spread, reduce negative effects to natural values, prevent residential housing losses, and avoid causalties. Though current fire policies in most countries demand full suppression, few studies have attempted to identify the strategic locations where firefighting efforts would likely contain catastrophic fire events. The success in containing those fires that escape initial attack is determined by diverse structural factors such as ground accessibility, airborne support, barriers to surface fire spread, and vegetation impedance. In this study, we predicted the success in fire containment across Catalonia (northeastern Spain) using a model generated with random forest from detailed geospatial data and a set of 73 fire perimeters for the period 2008-2016. The model attained a high predictive performance (AUC = 0.88), and the results were provided at fine resolution (25 m) for the entire study area (32,108 km2 ). The highest success rates were found in agricultural plains along the nonburnable barriers such as major road corridors and largest rivers. Low levels of containment likelihood were predicted for dense forest lands and steep-relief mountainous areas. The results can assist in suppression resource pre-positioning and extended attack decision making, but also in strategic fuels management oriented at creating defensive locations and fragmenting the landscape in operational firefighting areas. Our modeling workflow and methods may serve as a baseline to generate locally adapted models in fire-prone areas elsewhere.

The dental proteome of Homo antecessor.

The phylogenetic relationships between hominins of the Early Pleistocene epoch in Eurasia, such as Homo antecessor, and hominins that appear later in the fossil record during the Middle Pleistocene epoch, such as Homo sapiens, are highly debated1-5. For the oldest remains, the molecular study of these relationships is hindered by the degradation of ancient DNA. However, recent research has demonstrated that the analysis of ancient proteins can address this challenge6-8. Here we present the dental enamel proteomes of H. antecessor from Atapuerca (Spain)9,10 and Homo erectus from Dmanisi (Georgia)1, two key fossil assemblages that have a central role in models of Pleistocene hominin morphology, dispersal and divergence. We provide evidence that H. antecessor is a close sister lineage to subsequent Middle and Late Pleistocene hominins, including modern humans, Neanderthals and Denisovans. This placement implies that the modern-like face of H. antecessor-that is, similar to that of modern humans-may have a considerably deep ancestry in the genus Homo, and that the cranial morphology of Neanderthals represents a derived form. By recovering AMELY-specific peptide sequences, we also conclude that the H. antecessor molar fragment from Atapuerca that we analysed belonged to a male individual. Finally, these H. antecessor and H. erectus fossils preserve evidence of enamel proteome phosphorylation and proteolytic digestion that occurred in vivo during tooth formation. Our results provide important insights into the evolutionary relationships between H. antecessor and other hominin groups, and pave the way for future studies using enamel proteomes to investigate hominin biology across the existence of the genus Homo.

Mapping co-ancestry connections between the genome of a Medieval individual and modern Europeans.

Historical genetic links among similar populations can be difficult to establish. Identity by descent (IBD) analyses find genomic blocks that represent direct genealogical relationships among individuals. However, this method has rarely been applied to ancient genomes because IBD stretches are progressively fragmented by recombination and thus not recognizable after few tens of generations. To explore such genealogical relationships, we estimated long IBD blocks among modern Europeans, generating networks to uncover the genetic structures. We found that Basques, Sardinians, Icelanders and Orcadians form, each of them, highly intraconnected sub-clusters in a European network, indicating dense genealogical links within small, isolated populations. We also exposed individual genealogical links -such as the connection between one Basque and one Icelandic individual- that cannot be uncovered with other, widely used population genetics methods such as PCA or ADMIXTURE. Moreover, using ancient DNA technology we sequenced a Late Medieval individual (Barcelona, Spain) to high genomic coverage and identified IBD blocks shared between her and modern Europeans. The Medieval IBD blocks are statistically overrepresented only in modern Spaniards, which is the geographically closest population. This approach can be used to produce a fine-scale reflection of shared ancestry across different populations of the world, offering a direct genetic link from the past to the present.

Differential DNA methylation of vocal and facial anatomy genes in modern humans.

Changes in potential regulatory elements are thought to be key drivers of phenotypic divergence. However, identifying changes to regulatory elements that underlie human-specific traits has proven very challenging. Here, we use 63 reconstructed and experimentally measured DNA methylation maps of ancient and present-day humans, as well as of six chimpanzees, to detect differentially methylated regions that likely emerged in modern humans after the split from Neanderthals and Denisovans. We show that genes associated with face and vocal tract anatomy went through particularly extensive methylation changes. Specifically, we identify widespread hypermethylation in a network of face- and voice-associated genes (SOX9, ACAN, COL2A1, NFIX and XYLT1). We propose that these repression patterns appeared after the split from Neanderthals and Denisovans, and that they might have played a key role in shaping the modern human face and vocal tract.

Evolutionary History, Genomic Adaptation to Toxic Diet, and Extinction of the Carolina Parakeet.

As the only endemic neotropical parrot to have recently lived in the northern hemisphere, the Carolina parakeet (Conuropsis carolinensis) was an iconic North American bird. The last surviving specimen died in the Cincinnati Zoo in 1918 [1]. The cause of its extinction remains contentious: besides excessive mortality associated to habitat destruction and active hunting, their survival could have been negatively affected by its range having become increasingly patchy [2] or by the exposure to poultry pathogens [3, 4]. In addition, the Carolina parakeet showed a predilection for cockleburs, an herbaceous plant that contains a powerful toxin, carboxyatractyloside, or CAT [5], which did not seem to affect them but made the birds notoriously toxic to most predators [3]. To explore the demographic history of this bird, we generated the complete genomic sequence of a preserved specimen held in a private collection in Espinelves (Girona, Spain), as well as of a close extant relative, Aratinga solstitialis. We identified two non-synonymous genetic changes in two highly conserved proteins known to interact with CAT that could underlie a specific dietary adaptation to this toxin. Our genomic analyses did not reveal evidence of a dramatic past demographic decline in the Carolina parakeet; also, its genome did not exhibit the long runs of homozygosity that are signals of recent inbreeding and are typically found in endangered species. As such, our results suggest its extinction was an abrupt process and thus likely solely attributable to human causes.

Plasmodium vivax Malaria Viewed through the Lens of an Eradicated European Strain.

The protozoan Plasmodium vivax is responsible for 42% of all cases of malaria outside Africa. The parasite is currently largely restricted to tropical and subtropical latitudes in Asia, Oceania, and the Americas. Though, it was historically present in most of Europe before being finally eradicated during the second half of the 20th century. The lack of genomic information on the extinct European lineage has prevented a clear understanding of historical population structuring and past migrations of P. vivax. We used medical microscope slides prepared in 1944 from malaria-affected patients from the Ebro Delta in Spain, one of the last footholds of malaria in Europe, to generate a genome of a European P. vivax strain. Population genetics and phylogenetic analyses placed this strain basal to a cluster including samples from the Americas. This genome allowed us to calibrate a genomic mutation rate for P. vivax, and to estimate the mean age of the last common ancestor between European and American strains to the 15th century. This date points to an introduction of the parasite during the European colonization of the Americas. In addition, we found that some known variants for resistance to antimalarial drugs, including Chloroquine and Sulfadoxine, were already present in this European strain, predating their use. Our results shed light on the evolution of an important human pathogen and illustrate the value of antique medical collections as a resource for retrieving genomic information on pathogens from the past.