Search

Environmental DNA (eDNA) is the fastest growing biomonitoring tool fuelled by two key features: time efficiency and sensitivity. Technological advancements allow rapid biodiversity detection at both species and community levels with increasing accuracy. Concurrently, there has been a global demand to standardise eDNA methods, but this is only possible with an in-depth overview of the technological advancements and a discussion of the pros and cons of available methods. We therefore conducted a systematic literature review of 407 peer-reviewed papers on aquatic eDNA published between 2012 and 2021. We observed a gradual increase in the annual number of publications from four (2012) to 28 (2018), followed by a rapid growth to 124 publications in 2021. This was mirrored by a tremendous diversification of methods in all aspects of the eDNA workflow. For example, in 2012 only freezing was applied to preserve filter samples, whereas we recorded 12 different preservation methods in the 2021 literature. Despite an ongoing standardisation debate in the eDNA community, the field is seemingly moving fast in the opposite direction and we discuss the reasons and implications. Moreover, by compiling the largest PCR-primer database to date, we provide information on 522 and 141 published species-specific and metabarcoding primers targeting a wide range of aquatic organisms. This works as a user-friendly ‘distillation’ of primer information that was hitherto scattered across hundreds of papers, but the list also reflects which taxa are commonly studied with eDNA technology in aquatic environments such as fish and amphibians, and reveals that groups such as corals, plankton and algae are under-studied. Efforts to improve sampling and extraction methods, primer specificity and reference databases are crucial to capture these ecologically important taxa in future eDNA biomonitoring surveys. In a rapidly diversifying field, this review synthetises aquatic eDNA procedures and can guide eDNA users towards best practice.

Gigantopithecus blacki was a giant hominid that inhabited densely forested environments of Southeast Asia during the Pleistocene epoch1. Its evolutionary relationships to other great ape species, and the divergence of these species during the Middle and Late Miocene epoch (16-5.3 million years ago), remain unclear2,3. Hypotheses regarding the relationships between Gigantopithecus and extinct and extant hominids are wide ranging but difficult to substantiate because of its highly derived dentognathic morphology, the absence of cranial and post-cranial remains1,3-6, and the lack of independent molecular validation. We retrieved dental enamel proteome sequences from a 1.9-million-year-old G. blacki molar found in Chuifeng Cave, China7,8. The thermal age of these protein sequences is approximately five times greater than that of any previously published mammalian proteome or genome. We demonstrate that Gigantopithecus is a sister clade to orangutans (genus Pongo) with a common ancestor about 12-10 million years ago, implying that the divergence of Gigantopithecus from Pongo forms part of the Miocene radiation of great apes. In addition, we hypothesize that the expression of alpha-2-HS-glycoprotein, which has not been previously observed in enamel proteomes, had a role in the biomineralization of the thick enamel crowns that characterize the large molars in Gigantopithecus9,10. The survival of an Early Pleistocene dental enamel proteome in the subtropics further expands the scope of palaeoproteomic analysis into geographical areas and time periods previously considered incompatible with the preservation of substantial amounts of genetic information. E.C. and F.W. are supported by VILLUM FONDEN (17649) and by the European Commission through a Marie Skłodowska-Curie (MSCA) Individual Fellowship (795569). T.M.-B. is supported by BFU2017-86471-P (MINECO/FEDER, UE), NIHM grant U01 MH106874, Howard Hughes International Early Career grant, Obra Social ‘La Caixa’ and Secretaria d’Universitats i Recerca and CERCA Programme del Departament d’Economia i Coneixement de la Generalitat de Catalunya (GRC 2017 SGR 880). E.C., J.C., J.V.O., D.S. and P.G. are supported by the Marie Skłodowska-Curie European Training Network (ETN) TEMPERA, a project funded by the European Union’s EU Framework Program for Research and Innovation Horizon 2020 under grant agreement no. 722606. M.J.C. and M.M. are supported by the Danish National Research Foundation award PROTEIOS (DNRF128). Work at the Novo Nordisk Foundation Center for Protein Research is funded in part by a donation from the Novo Nordisk Foundation (NNF14CC0001). Research at Chuifeng Cave is made possible by support from the National Natural Science Foundation of China (41572023) and by a grant from the Bagui Scholar of Guangxi. M.K. was supported by a Deutsche Forschungsgemeinschaft (DFG) fellowship (KU 3467/1-1) and the Postdoctoral Junior Leader Fellowship Programme from ‘la Caixa’ Banking Foundation (LCF/BQ/PR19/11700002). M.E.A. is supported by the Independent Research Fund Denmark (7027-00147B). We thank E. Willerslev for critical reading of the manuscript, scientific support and guidance.

Summary The Eurasian Holocene (beginning c. 12 thousand years ago) encompassed some of the most significant changes in human evolution, with far-reaching consequences for the dietary, physical and mental health of present-day populations. Using an imputed dataset of >1600 complete ancient genome sequences, and new computational methods for locating selection in time and space, we reconstructed the selection landscape of the transition from hunting and gathering, to farming and pastoralism across West Eurasia. We identify major selection signals related to metabolism, possibly associated with the dietary shift occurring in this period. We show that the selection on loci such as the FADS cluster, associated with fatty acid metabolism, and the lactase persistence locus, began earlier than previously thought. A substantial amount of selection is also found in the HLA region and other loci associated with immunity, possibly due to the increased exposure to pathogens during the Neolithic, which may explain the current high prevalence of auto-immune disease, such as psoriasis, due to genetic trade-offs. By using ancient populations to infer local ancestry tracks in hundreds of thousands of samples from the UK Biobank, we find strong genetic differentiation among ancient Europeans in loci associated with anthropometric traits and susceptibility to several diseases that contribute to present-day disease burden. These were previously thought to be caused by local selection, but in fact can be attributed to differential genetic contributions from various source populations that are ancestral to present-day Europeans. Thus, alleles associated with increased height seem to have increased in frequency following the Yamnaya migration into northwestern Europe around 5,000 years ago. Alleles associated with increased risk of some mood-related phenotypes are overrepresented in the farmer ancestry component entering Europe from Anatolia around 11,000 years ago, while western hunter-gatherers show a strikingly high contribution of alleles conferring risk of traits related to diabetes. Our results paint a picture of the combined contributions of migration and selection in shaping the phenotypic landscape of present-day Europeans that suggests a combination of ancient selection and migration, rather than recent local selection, is the primary driver of present-day phenotypic differences in Europe.

Due to postmortem DNA degradation, most ancient genomes sequenced to date have low depth of coverage, preventing the true underlying genotypes from being recovered. Genotype imputation has been put forward to improve genotyping accuracy for low-coverage genomes. However, it is unknown to what extent imputation of ancient genomes produces accurate genotypes and whether imputation introduces bias to downstream analyses. To address these questions, we downsampled 43 ancient genomes, 42 of which are high-coverage (above 10x) and three constitute a trio (mother, father and son), from different times and continents to simulate data with coverage in the range of 0.1x-2.0x and imputed these using state-of-the-art methods and reference panels. We assessed imputation accuracy across ancestries and depths of coverage. We found that ancient and modern DNA imputation accuracies were comparable. We imputed most of the 42 high-coverage genomes downsampled to 1x with low error rates (below 5%) and estimated higher error rates for African genomes, which are underrepresented in the reference panel. We used the ancient trio data to validate imputation and phasing results using an orthogonal approach based on Mendel’s rules of inheritance. This resulted in imputation and switch error rates of 1.9% and 2.0%, respectively, for 1x genomes. We further compared the results of downstream analyses between imputed and high-coverage genomes, notably principal component analysis (PCA), genetic clustering, and runs of homozygosity (ROH). For these three approaches, we observed similar results between imputed and high-coverage genomes using depths of coverage of at least 0.5x, except for African genomes, for which the decreased imputation accuracy impacted ROH estimates. Altogether, these results suggest that, for most populations and depths of coverage as low as 0.5x, imputation is a reliable method with potential to expand and improve ancient DNA studies.

Although Brazil was inhabited by more than 3,000 Indigenous populations prior to European colonization, today’s Indigenous peoples represent less than 1% of Brazil’s census population. Some of the decimated communities belonged to the so-called “Botocudos” from central-eastern Brazil. These peoples are thought to represent a case of long-standing genetic continuity bearing a strong craniometric resemblance to that of the oldest Indigenous Americans (“Paleoamericans”). Yet, little is known about their origins and genetic relationship to other Native Americans, as only two “Botocudo” genomes have been sequenced so far and those were surprisingly of Polynesian ancestry. To deepen our knowledge on the genomic history of pre-contact Indigenous Americans and the pathogens they were exposed to, we carbon-dated and sequenced 24 ancient Brazilians (including 22 “Botocudos”) whose remains were hosted at the National Museum of Rio de Janeiro and recovered prior to the tragic 2018 fire. The resulting genomes’ depth of coverage ranged from 0.001× to 24×. Their genetic ancestry was found to be Indigenous American without gene flow from external populations such as Europeans, Africans or Polynesians. Unlike Mesoamericans, the “Botocudos” and Amazonians do not seem to have experienced a population expansion once in the Americas. Moreover, remarkably, their genomes exhibit amongst the lowest levels of heterozygosity worldwide and long runs of homozygosity, which could be explained by unique social practices or a very small effective size. Finally, whole genomes of likely ancient pathogens were recovered, including lineages of Human parvovirus B19 that were possibly introduced after the European contact.Significance statementTo better understand the genetic relationship among Indigenous populations in Brazil, we sequenced the genomes of 24 ancient individuals (22 of which labelled as “Botocudos”, a term used to describe hunter-gatherer tribes) whose remains were hosted at the Museu Nacional of Rio de Janeiro prior to the tragic fire that consumed it in 2018. Unlike two previously published “Botocudo” genomes, the 22 “Botocudos” from this study have Indigenous American-related ancestry without any Polynesian-related ancestry, and they are similarly related to several Native Brazilian populations. Finally, unlike Eurasian hunter-gatherers, the “Botocudos” exhibit among the lowest heterozygosity and longest runs of homozygosity worldwide – compatible with a very small effective size and suggesting a unique social structure among hunter-gatherers in the Americas.

Monitoring of biota is pivotal for the assessment and conservation of ecosystems. Environments worldwide are being continuously and increasingly exposed to multiple adverse impacts, and the accuracy and reliability of the biomonitoring tools that can be employed shape not only the present, but more importantly, the future of entire habitats. The analysis of environmental DNA (eDNA) metabarcoding data provides a quick, affordable, and reliable molecular approach for biodiversity assessments. However, while extensively employed in aquatic and terrestrial surface environments, eDNA-based studies targeting subterranean ecosystems are still uncommon due to the lack of accessibility and the cryptic nature of these environments and their species. Recent advances in genetic and genomic analyses have established a promising framework for shedding new light on subterranean biodiversity and ecology. To address current knowledge and the future use of eDNA methods in groundwaters and caves, this review explores conceptual and technical aspects of the application and its potential in subterranean systems. We briefly introduce subterranean biota and describe the most used traditional sampling techniques. Next, eDNA characteristics, application, and limitations in the subsurface environment are outlined. Last, we provide suggestions on how to overcome caveats and delineate some of the research avenues that will likely shape this field in the near future. We advocate that eDNA analyses, when carefully conducted and ideally combined with conventional sampling techniques, will substantially increase understanding and enable crucial expansion of subterranean community characterisation. Given the importance of groundwater and cave ecosystems for nature and humans, eDNA can bring to the surface essential insights, such as study of ecosystem assemblages and rare species detection, which are critical for the preservation of life below, as well as above, the ground.

The Eastern Arc Mountains in Tanzania represent a hotspot for biological diversity of global importance. The level of endemism is high, and Eastern Arc biodiversity has been studied extensively in vertebrates and invertebrates, including millipedes. However, millipede evolution is vastly understudied at the molecular level. Therefore, we used next-generation ‘shotgun’ sequencing to obtain mitochondrial genome sequences of 26 museum specimens, representing six genera and 12 millipede species found across the Eastern Arc Mountains. Bayesian and maximum likelihood methods yielded consistent topologies with high node support, confirming a high level of congruence between molecular and morphological analyses. The only exception was a Tropostreptus sigmatospinus individual from Zanzibar, which was placed outside an otherwise monophyletic group consisting of mainland individuals of the same assumed species. For two species with a distribution across several mountain blocks (Tropostreptus sigmatospinus and Tropostreptus hamatus), each mountain population represents a distinct monophyletic lineage. In contrast, we also observe that distinct species exist sympatrically in the same montane forests, indicative of older speciation events that are not defined by current forest distribution. Our results are important for understanding speciation mechanisms in montane rain forests and highlight that ethanol-preserved invertebrates exhibit a tremendous potential for genomic analyses.

Global food production, food supply chains and food security are increasingly stressed by human population growth and loss of arable land, becoming more vulnerable to anthropogenic and environmental perturbations. Numerous mutualistic and antagonistic species are interconnected with the cultivation of crops and livestock and these can be challenging to identify on the large scales of food production systems. Accurate identifications to capture this diversity and rapid scalable monitoring are necessary to identify emerging threats (i.e. pests and pathogens), inform on ecosystem health (i.e. soil and pollinator diversity), and provide evidence for new management practices (i.e. fertiliser and pesticide applications). Increasingly, environmental DNA (eDNA) is providing rapid and accurate classifications for specific organisms and entire species assemblages in substrates ranging from soil to air. Here, we aim to discuss how eDNA is being used for monitoring of agricultural ecosystems, what current limitations exist, and how these could be managed to expand applications into the future. In a systematic review we identify that eDNA-based monitoring in food production systems accounts for only 4 % of all eDNA studies. We found that the majority of these eDNA studies target soil and plant substrates (60 %), predominantly to identify microbes and insects (60 %) and are biased towards Europe (42 %). While eDNA-based monitoring studies are uncommon in many of the world's food production systems, the trend is most pronounced in emerging economies often where food security is most at risk. We suggest that the biggest limitations to eDNA for agriculture are false negatives resulting from DNA degradation and assay biases, as well as incomplete databases and the interpretation of abundance data. These require in silico, in vitro, and in vivo approaches to carefully design, test and apply eDNA monitoring for reliable and accurate taxonomic identifications. We explore future opportunities for eDNA research which could further develop this useful tool for food production system monitoring in both emerging and developed economies, hopefully improving monitoring, and ultimately food security.

Snake venoms represent a danger to human health, but also a gold mine of bioactive proteins that can be harnessed for drug discovery purposes. The evolution of snakes and their venom has been studied for decades, particularly via traditional morphological and basic genetic methods alongside venom proteomics. However, while the field of genomics has matured rapidly over the past 2 decades, owing to the development of next-generation sequencing technologies, snake genomics remains in its infancy. Here, we provide an overview of the state of the art in snake genomics and discuss its potential implications for studying venom evolution and toxinology. On the basis of current knowledge, gene duplication and positive selection are key mechanisms in the neofunctionalization of snake venom proteins. This makes snake venoms important evolutionary drivers that explain the remarkable venom diversification and adaptive variation observed in these reptiles. Gene duplication and neofunctionalization have also generated a large number of repeat sequences in snake genomes that pose a significant challenge to DNA sequencing, resulting in the need for substantial computational resources and longer sequencing read length for high-quality genome assembly. Fortunately, owing to constantly improving sequencing technologies and computational tools, we are now able to explore the molecular mechanisms of snake venom evolution in unprecedented detail. Such novel insights have the potential to affect the design and development of antivenoms and possibly other drugs, as well as provide new fundamental knowledge on snake biology and evolution.