Single-stranded DNA library preparation from highly degraded DNA using T4 DNA ligase

Ancient Plasmodium genomes shed light on the history of human malaria

Malaria-causing protozoa of the genus Plasmodium have exerted one of the strongest selective pressures on the human genome, and resistance alleles provide biomolecular footprints that outline the historical reach of these species1. Nevertheless, debate persists over when and how malaria parasites emerged as human pathogens and spread around the globe1,2. To address these questions, we generated high-coverage ancient mitochondrial and nuclear genome-wide data from P. falciparum, P. vivax and P. malariae from 16 countries spanning around 5,500 years of human history. We identified P. vivax and P. falciparum across geographically disparate regions of Eurasia from as early as the fourth and first millennia BCE, respectively; for P. vivax, this evidence pre-dates textual references by several millennia3. Genomic analysis supports distinct disease histories for P. falciparum and P. vivax in the Americas: similarities between now-eliminated European and peri-contact South American strains indicate that European colonizers were the source of American P. vivax, whereas the trans-Atlantic slave trade probably introduced P. falciparum into the Americas. Our data underscore the role of cross-cultural contacts in the dissemination of malaria, laying the biomolecular foundation for future palaeo-epidemiological research into the impact of Plasmodium parasites on human history. Finally, our unexpected discovery of P. falciparum in the high-altitude Himalayas provides a rare case study in which individual mobility can be inferred from infection status, adding to our knowledge of cross-cultural connectivity in the region nearly three millennia ago.

The Amphibian Genomics Consortium: advancing genomic and genetic resources for amphibian research and conservation

Amphibians represent a diverse group of tetrapods, marked by deep divergence times between their three systematic orders and families. Studying amphibian biology through the genomics lens increases our understanding of the features of this animal class and that of other terrestrial vertebrates. The need for amphibian genomics resources is more urgent than ever due to the increasing threats to this group. Amphibians are one of the most imperiled taxonomic groups, with approximately 41% of species threatened with extinction due to habitat loss, changes in land use patterns, disease, climate change, and their synergistic effects. Amphibian genomics resources have provided a better understanding of ontogenetic diversity, tissue regeneration, diverse life history and reproductive modes, antipredator strategies, and resilience and adaptive responses. They also serve as critical models for understanding widespread genomic characteristics, including evolutionary genome expansions and contractions given they have the largest range in genome sizes of any animal taxon and multiple mechanisms of genetic sex determination. Despite these features, genome sequencing of amphibians has significantly lagged behind that of other vertebrates, primarily due to the challenges of assembling their large, repeat-rich genomes and the relative lack of societal support. The advent of long-read sequencing technologies, along with computational techniques that enhance scaffolding capabilities and streamline computational workload is now enabling the ability to overcome some of these challenges. To promote and accelerate the production and use of amphibian genomics research through international coordination and collaboration, we launched the Amphibian Genomics Consortium (AGC) in early 2023. This burgeoning community already has more than 282 members from 41 countries (6 in Africa, 131 in the Americas, 27 in Asia, 29 in Australasia, and 89 in Europe). The AGC aims to leverage the diverse capabilities of its members to advance genomic resources for amphibians and bridge the implementation gap between biologists, bioinformaticians, and conservation practitioners. Here we evaluate the state of the field of amphibian genomics, highlight previous studies, present challenges to overcome, and outline how the AGC can enable amphibian genomics research to "leap" to the next level.

Coevolution with hosts underpins speciation in brood-parasitic cuckoos

Coevolution between interacting species is thought to increase biodiversity, but evidence linking microevolutionary processes to macroevolutionary patterns is scarce. We leveraged two decades of behavioral research coupled with historical DNA analysis to reveal that coevolution with hosts underpins speciation in brood-parasitic bronze-cuckoos. At a macroevolutionary scale, we show that highly virulent brood-parasitic taxa have higher speciation rates and are more likely to speciate in sympatry than less-virulent and nonparasitic relatives. We reveal the microevolutionary process underlying speciation: Hosts reject cuckoo nestlings, which selects for mimetic cuckoo nestling morphology. Where cuckoos exploit multiple hosts, selection for mimicry drives genetic and phenotypic divergence corresponding to host preference, even in sympatry. Our work elucidates perhaps the most common, but poorly characterized, evolutionary process driving biological diversification.

Network of large pedigrees reveals social practices of Avar communities

From AD 567-568, at the onset of the Avar period, populations from the Eurasian Steppe settled in the Carpathian Basin for approximately 250 years1. Extensive sampling for archaeogenomics (424 individuals) and isotopes, combined with archaeological, anthropological and historical contextualization of four Avar-period cemeteries, allowed for a detailed description of the genomic structure of these communities and their kinship and social practices. We present a set of large pedigrees, reconstructed using ancient DNA, spanning nine generations and comprising around 300 individuals. We uncover a strict patrilineal kinship system, in which patrilocality and female exogamy were the norm and multiple reproductive partnering and levirate unions were common. The absence of consanguinity indicates that this society maintained a detailed memory of ancestry over generations. These kinship practices correspond with previous evidence from historical sources and anthropological research on Eurasian Steppe societies2. Network analyses of identity-by-descent DNA connections suggest that social cohesion between communities was maintained via female exogamy. Finally, despite the absence of major ancestry shifts, the level of resolution of our analyses allowed us to detect genetic discontinuity caused by the replacement of a community at one of the sites. This was paralleled with changes in the archaeological record and was probably a result of local political realignment.

A review on the impact of single-stranded library preparation on plasma cell-free diversity for cancer detection

In clinical oncology, cell-free DNA (cfDNA) has shown immense potential in its ability to noninvasively detect cancer at various stages and monitor the progression of therapy. Despite the rapid improvements in cfDNA liquid biopsy approaches, achieving the required sensitivity to detect rare tumor-derived cfDNA still remains a challenge. For next-generation sequencing, the perceived presentation of cfDNA is strongly linked to the extraction and library preparation protocols. Conventional double-stranded DNA library preparation (dsDNA-LP) focuses on assessing ~167bp double-stranded mononucleosomal (mncfDNA) and its other oligonucleosomal cell-free DNA counterparts in plasma. However, dsDNA-LP methods fail to include short, single-stranded, or nicked DNA in the final library preparation, biasing the representation of the actual cfDNA populations in plasma. The emergence of single-stranded library preparation (ssDNA-LP) strategies over the past decade has now allowed these other populations of cfDNA to be studied from plasma. With the use of ssDNA-LP, single-stranded, nicked, and ultrashort cfDNA can be comprehensively assessed for its molecular characteristics and clinical potential. In this review, we overview the current literature on applications of ssDNA-LP on plasma cfDNA from a potential cancer liquid biopsy perspective. To this end, we discuss the molecular principles of single-stranded DNA adapter ligation, how library preparation contributes to the understanding of native cfDNA characteristics, and the potential for ssDNA-LP to improve the sensitivity of circulating tumor DNA detection. Additionally, we review the current literature on the newly reported species of plasma ultrashort single-stranded cell-free DNA plasma, which appear biologically distinct from mncfDNA. We conclude with a discussion of future perspectives of ssDNA-LP for liquid biopsy endeavors.

Targeted enrichment of whole-genome SNPs from highly burned skeletal remains

Genetic assessment of highly incinerated and/or degraded human skeletal material is a persistent challenge in forensic DNA analysis, including identifying victims of mass disasters. Few studies have investigated the impact of thermal degradation on whole-genome single-nucleotide polymorphism (SNP) quality and quantity using next-generation sequencing (NGS). We present whole-genome SNP data obtained from the bones and teeth of 27 fire victims using two DNA extraction techniques. Extracts were converted to double-stranded DNA libraries then enriched for whole-genome SNPs using unpublished biotinylated RNA baits and sequenced on an Illumina NextSeq 550 platform. Raw reads were processed using the EAGER (Efficient Ancient Genome Reconstruction) pipeline, and the SNPs filtered and called using FreeBayes and GATK (v. 3.8). Mixed-effects modeling of the data suggest that SNP variability and preservation is predominantly determined by skeletal element and burn category, and not by extraction type. Whole-genome SNP data suggest that selecting long bones, hand and foot bones, and teeth subjected to temperatures <350°C are the most likely sources for higher genomic DNA yields. Furthermore, we observed an inverse correlation between the number of captured SNPs and the extent to which samples were burned, as well as a significant decrease in the total number of SNPs measured for samples subjected to temperatures >350°C. Our data complement previous analyses of burned human remains that compare extraction methods for downstream forensic applications and support the idea of adopting a modified Dabney extraction technique when traditional forensic methods fail to produce DNA yields sufficient for genetic identification.

Ancient mitochondrial genomes unveil the origins and evolutionary history of New Zealand's enigmatic takahē and moho

Many avian species endemic to Aotearoa New Zealand were driven to extinction or reduced to relict populations following successive waves of human arrival, due to hunting, habitat destruction and the introduction of mammalian predators. Among the affected species were the large flightless South Island takahē (Porphyrio hochstetteri) and the moho (North Island takahē; P. mantelli), with the latter rendered extinct and the former reduced to a single relictual population. Little is known about the evolutionary history of these species prior to their decline and/or extinction. Here we sequenced mitochondrial genomes from takahē and moho subfossils (12 takahē and 4 moho) and retrieved comparable sequence data from takahē museum skins (n = 5) and contemporary individuals (n = 17) to examine the phylogeny and recent evolutionary history of these species. Our analyses suggest that prehistoric takahē populations lacked deep phylogeographic structure, in contrast to moho, which exhibited significant spatial genetic structure, albeit based on limited sample sizes (n = 4). Temporal genetic comparisons show that takahē have lost much of their mitochondrial genetic diversity, likely due to a sudden demographic decline soon after human arrival (~750 years ago). Time-calibrated phylogenetic analyses strongly support a sister species relationship between takahē and moho, suggesting these flightless taxa diverged around 1.5 million years ago, following a single colonisation of New Zealand by a flighted Porphyrio ancestor approximately 4 million years ago. This study highlights the utility of palaeogenetic approaches for informing the conservation and systematic understanding of endangered species whose ranges have been severely restricted by anthropogenic impacts.

Diet in phenotypically divergent sympatric species of African weakly electric fish (genus: Campylomormyrus)-A hybrid capture/HTS metabarcoding approach

Ecological speciation within the mormyrid genus Campylomormyrus resulted in sympatric species exhibiting divergence in their feeding apparatus and electric organ discharge (EOD). This study documents the overall diet of the genus Campylomormyrus and examines the hypothesis that the Campylomormyrus radiation is caused by adaptation to different food sources. We performed diet assessment of five sympatric Campylomormyrus species (C. alces, C. compressirostris, C. curvirostris, C. tshokwe, C. numenius) and their sister taxon Gnathonemus petersii with markedly different snout morphologies and EODs using hybrid capture/HTS DNA metabarcoding of their stomach contents. Our approach allowed for high taxonomic resolution of prey items, including benthic invertebrates, allochthonous invertebrates and vegetation. Comparisons of the diet compositions using quantitative measures and diet overlap indices revealed that all species are able to exploit multiple food niches in their habitats, that is fauna at the bottom, the water surface and the water column. A major part of the diet is larvae of aquatic insects, such as dipterans, coleopterans and trichopterans, known to occur in holes and interstitial spaces of the substrate. The results indicate that different snout morphologies and the associated divergence in the EOD could translate into different prey spectra. This suggests that the diversification in EOD and/or morphology of the feeding apparatus could be under functional adaptation.

Ultra-short cell-free DNA fragments enhance cancer early detection in a multi-analyte blood test combining mutation, protein and fragmentomics

OBJECTIVES

Cancer morbidity and mortality can be reduced if the cancer is detected early. Cell-free DNA (cfDNA) fragmentomics emerged as a novel epigenetic biomarker for early cancer detection, however, it is still at its infancy and requires technical improvement. We sought to apply a single-strand DNA sequencing technology, for measuring genetic and fragmentomic features of cfDNA and evaluate the performance in detecting multiple cancers.

METHODS

Blood samples of 364 patients from six cancer types (colorectal, esophageal, gastric, liver, lung, and ovarian cancers) and 675 healthy individuals were included in this study. Circulating tumor DNA mutations, cfDNA fragmentomic features and a set of protein biomarkers were assayed. Sensitivity and specificity were reported by cancer types and stages.

RESULTS

Circular Ligation Amplification and sequencing (CLAmp-seq), a single-strand DNA sequencing technology, yielded a population of ultra-short fragments (<100 bp) than double-strand DNA preparation protocols and reveals a more significant size difference between cancer and healthy cfDNA fragments (25.84 bp vs. 16.05 bp). Analysis of the subnucleosomal peaks in ultra-short cfDNA fragments indicates that these peaks are regulatory element "footprints" and correlates with gene expression and cancer stages. At 98 % specificity, a prediction model using ctDNA mutations alone showed an overall sensitivity of 46 %; sensitivity reaches 60 % when protein is added, sensitivity further increases to 66 % when fragmentomics is also integrated. More improvements observed for samples representing earlier cancer stages than later ones.

CONCLUSIONS

These results suggest synergistic properties of protein, genetic and fragmentomics features in the identification of early-stage cancers.

Wide-spread dispersal in a deep-sea brooding polychaete: the role of natural history collections in assessing the distribution in quill worms (Onuphidae, Annelida)

BACKGROUND

Modern integrative taxonomy-based annelid species descriptions are detailed combining morphological data and, since the last decades, also molecular information. Historic species descriptions are often comparatively brief lacking such detail. Adoptions of species names from western literature in the past led to the assumption of cosmopolitan ranges for many species, which, in many cases, were later found to include cryptic or pseudocryptic lineages with subtle morphological differences. Natural history collections and databases can aid in assessing the geographic ranges of species but depend on correct species identification. Obtaining DNA sequence information from wet-collection museum specimens of marine annelids is often impeded by the use of  formaldehyde and/or long-term storage in ethanol resulting in DNA degradation and cross-linking.

RESULTS

The application of ancient DNA extraction methodology in combination with single-stranded DNA library preparation and target gene capture resulted in successful sequencing of a 110-year-old collection specimen of quill worms. Furthermore, a 40-year-old specimen of quill worms was successfully sequenced using a standard extraction protocol for modern samples, PCR and Sanger sequencing. Our study presents the first molecular analysis of Hyalinoecia species including the previously known species Hyalinoecia robusta, H. tubicloa, H. artifex, and H. longibranchiata, and a potentially undescribed species from equatorial western Africa morphologically indistinguishable from H. tubicola. The study also investigates the distribution of these five Hyalinoecia species. Reassessing the distribution of H. robusta reveals a geographical range covering both the Atlantic and the Indian Oceans as indicated by molecular data obtained from recent and historical specimens.

CONCLUSION

Our results represent an example of a very wide geographical distribution of a brooding deep-sea annelid with a complex reproduction strategy and seemingly very limited dispersal capacity of its offspring, and highlights the importance of molecular information from museum specimens for integrative annelid taxonomy and biogeography.

Plankton community changes during the last 124 000 years in the subarctic Bering Sea derived from sedimentary ancient DNA

Current global warming results in rising sea-water temperatures, and the loss of sea ice in Arctic and subarctic oceans impacts the community composition of primary producers with cascading effects on the food web and potentially on carbon export rates. This study analyzes metagenomic shotgun and diatom rbcL amplicon sequencing data from sedimentary ancient DNA of the subarctic western Bering Sea that records phyto- and zooplankton community changes over the last glacial-interglacial cycles, including the last interglacial period (Eemian). Our data show that interglacial and glacial plankton communities differ, with distinct Eemian and Holocene plankton communities. The generally warm Holocene period is dominated by picosized cyanobacteria and bacteria-feeding heterotrophic protists, while the Eemian period is dominated by eukaryotic picosized chlorophytes and Triparmaceae. By contrast, the glacial period is characterized by microsized phototrophic protists, including sea ice-associated diatoms in the family Bacillariaceae and co-occurring diatom-feeding crustaceous zooplankton. Our deep-time record of plankton community changes reveals a long-term decrease in phytoplankton cell size coeval with increasing temperatures, resembling community changes in the currently warming Bering Sea. The phytoplankton community in the warmer-than-present Eemian period is distinct from modern communities and limits the use of the Eemian as an analog for future climate scenarios. However, under enhanced future warming, the expected shift toward the dominance of small-sized phytoplankton and heterotrophic protists might result in an increased productivity, whereas the community's potential of carbon export will be decreased, thereby weakening the subarctic Bering Sea's function as an effective carbon sink.

Genome sequences of 36,000- to 37,000-year-old modern humans at Buran-Kaya III in Crimea

Populations genetically related to present-day Europeans first appeared in Europe at some point after 38,000-40,000 years ago, following a cold period of severe climatic disruption. These new migrants would eventually replace the pre-existing modern human ancestries in Europe, but initial interactions between these groups are unclear due to the lack of genomic evidence from the earliest periods of the migration. Here we describe the genomes of two 36,000-37,000-year-old individuals from Buran-Kaya III in Crimea as belonging to this newer migration. Both genomes share the highest similarity to Gravettian-associated individuals found several thousand years later in southwestern Europe. These genomes also revealed that the population turnover in Europe after 40,000 years ago was accompanied by admixture with pre-existing modern human populations. European ancestry before 40,000 years ago persisted not only at Buran-Kaya III but is also found in later Gravettian-associated populations of western Europe and Mesolithic Caucasus populations.

Ancient DNA is preserved in fish fossils from tropical lake sediments

Tropical freshwater lakes are well known for their high biodiversity, and particularly the East African Great Lakes are renowned for their adaptive radiation of cichlid fishes. While comparative phylogenetic analyses of extant species flocks have revealed patterns and processes of their diversification, little is known about evolutionary trajectories within lineages, the impacts of environmental drivers, or the scope and nature of now-extinct diversity. Time-structured palaeodata from geologically young fossil records, such as fossil counts and particularly ancient DNA (aDNA) data, would help fill this large knowledge gap. High ambient temperatures can be detrimental to the preservation of DNA, but refined methodology now allows data generation even from very poorly preserved samples. Here, we show for the first time that fish fossils from tropical lake sediments yield endogenous aDNA. Despite generally low endogenous content and high sample dropout, the application of high-throughput sequencing and, in some cases, sequence capture allowed taxonomic assignment and phylogenetic placement of 17% of analysed fish fossils to family or tribe level, including remains which are up to 2700 years old or weigh less than 1 mg. The relationship between aDNA degradation and the thermal age of samples is similar to that described for terrestrial samples from cold environments when adjusted for elevated temperature. Success rates and aDNA preservation differed between the investigated lakes Chala, Kivu and Victoria, possibly caused by differences in bottom water oxygenation. Our study demonstrates that the sediment records of tropical lakes can preserve genetic information on rapidly diversifying fish taxa over time scales of millennia.

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.

Gray versus yellow ventral coloration: Identity, distribution, color polymorphism and molecular relationships of the microhylid frog Platypelis mavomavo Andreone, Fenolio & Walvoord, 2003

The Malagasy frog Platypelis mavomavo from Ambolokopatrika in the North East of Madagascar was originally diagnosed based on its bright yellow venter, but only limited information on this species has become available after its initial description in 2003. Several Platypelis specimens with yellow ventral color have been erroneously assigned to this species due to a lack of DNA sequences from the P. mavomavo type series. On the other hand, the candidate species Platypelis sp. Ca10 from Andranomapanga in the Northern Central East of Madagascar with gray ventral color has been defined based on its genetic differentiation from other nominal Platypelis species. Here we study the genetic variation of P. mavomavo and P. sp. Ca10 based on mitochondrial (16S rRNA) and nuclear-encoded (RAG-1) genes, including a newly determined sequence from the P. mavomavo holotype, which was studied using a museomics approach. We find only limited genetic variation among the samples studied, and this variation is unlinked to ventral coloration but instead reflects geographic distribution. We, therefore, conclude that P. sp. Ca10 is a gray-colored variant of P. mavomavo, and that P. mavomavo is rather widespread in the North East and Northern Central East of Madagascar, with populations in areas bordering the North West (Ambohitantely) and Sambirano (Ampotsidy) geographic regions, and the yellow-bellied morph restricted to the North East (Makira, Ambolokopatrika). Due to the range extension of P. mavomavo, the conservation status of the species requires re-assessment.

A highly efficient scheme for library preparation from single-stranded DNA

Although methods for sequencing library preparation from double-stranded DNA are well established, those from single-stranded DNA (ssDNA) have not been well studied. Further, the existing methods have limitations in efficiency and yield. Therefore, we developed a highly efficient procedure for sequencing library preparation from ssDNA. In this method, the first adaptor tagging of ssDNA is performed using terminal deoxyribonucleotidyl transferase (TdT)-assisted adenylate connector-mediated ssDNA (TACS) ligation, which we reported recently. After complementary strand synthesis using the adaptor-tagged ssDNA, second adaptor tagging via Vaccinia virus topoisomerase I (VTopoI or TOPO)-based adaptor ligation is performed. With additional steps for degradation, repression, and removal of the adaptor dimer, the proposed TACS-TOPO scheme realizes adaptor dimer-free sequencing library preparation from ssDNA samples of 24 pg. The TACS-TOPO scheme was successfully applied to cell-free DNA analysis with amplification-free library preparation from 50 µL of human serum. A modified TACS-TOPO scheme was also applied to DNA extracted from ancient human bones, bringing two to eight times more library yields than those using a conventional library preparation protocol. The procedures for preparing VTopoI and its complex with a double-stranded oligonucleotide adaptor are also described. Overall, the proposed TACS-TOPO scheme can facilitate practical and sensitive sequencing analysis of ssDNA.

Library preparation method and DNA source influence endogenous DNA recovery from 100-year-old avian museum specimens

Museum specimens collected prior to cryogenic tissue storage are increasingly being used as genetic resources, and though high-throughput sequencing is becoming more cost-efficient, whole genome sequencing (WGS) of historical DNA (hDNA) remains inefficient and costly due to its short fragment sizes and high loads of exogenous DNA, among other factors. It is also unclear how sequencing efficiency is influenced by DNA sources. We aimed to identify the most efficient method and DNA source for collecting WGS data from avian museum specimens. We analyzed low-coverage WGS from 60 DNA libraries prepared from four American Robin (Turdus migratorius) and four Abyssinian Thrush (Turdus abyssinicus) specimens collected in the 1920s. We compared DNA source (toepad versus incision-line skin clip) and three library preparation methods: (1) double-stranded DNA (dsDNA), single tube (KAPA); (2) single-stranded DNA (ssDNA), multi-tube (IDT); and (3) ssDNA, single tube (Claret Bioscience). We found that the ssDNA, multi-tube method resulted in significantly greater endogenous DNA content, average read length, and sequencing efficiency than the other tested methods. We also tested whether a predigestion step reduced exogenous DNA in libraries from one specimen per species and found promising results that warrant further study. The ~10% increase in average sequencing efficiency of the best-performing method over a commonly implemented dsDNA library preparation method has the potential to significantly increase WGS coverage of hDNA from bird specimens. Future work should evaluate the threshold for specimen age at which these results hold and how the combination of library preparation method and DNA source influence WGS in other taxa.

DEtail-seq is an ultra-efficient and convenient method for meiotic DNA break profiling in multiple organisms

Programmed DNA double-strand break (DSB) formation is a crucial step in meiotic recombination, yet techniques for high-efficiency and precise mapping of the 3' ends of DSBs are still in their infancy. Here, we report a novel technique, named DNA End tailing and sequencing (DEtail-seq), which can directly and ultra-efficiently characterize the 3' ends of meiotic DSBs with near single-nucleotide resolution in a variety of species, including yeast, mouse, and human. We find that the 3' ends of meiotic DSBs are stable without significant resection in budding yeast. Meiotic DSBs are strongly enriched in de novo H3K4me3 peaks in the mouse genome at leptotene stage. We also profile meiotic DSBs in human and find DSB hotspots are enriched near the common fragile sites during human meiosis, especially at CCCTC-binding factor (CTCF)-associated enhancers. Therefore, DEtail-seq provides a powerful method to detect DSB ends in various species, and our results provide new insights into the distribution and regulation of meiotic DSB hotspots.

Marine ecosystem shifts with deglacial sea-ice loss inferred from ancient DNA shotgun sequencing

Sea ice is a key factor for the functioning and services provided by polar marine ecosystems. However, ecosystem responses to sea-ice loss are largely unknown because time-series data are lacking. Here, we use shotgun metagenomics of marine sedimentary ancient DNA off Kamchatka (Western Bering Sea) covering the last ~20,000 years. We traced shifts from a sea ice-adapted late-glacial ecosystem, characterized by diatoms, copepods, and codfish to an ice-free Holocene characterized by cyanobacteria, salmon, and herring. By providing information about marine ecosystem dynamics across a broad taxonomic spectrum, our data show that ancient DNA will be an important new tool in identifying long-term ecosystem responses to climate transitions for improvements of ocean and cryosphere risk assessments. We conclude that continuing sea-ice decline on the northern Bering Sea shelf might impact on carbon export and disrupt benthic food supply and could allow for a northward expansion of salmon and Pacific herring.

Advancing Grain Legumes Domestication and Evolution Studies with Genomics

Grain legumes were domesticated in parallel with cereals in several regions of the world and formed the economic basis of early farming cultures. Since then, legumes have played a vital role in human and animal diets and in fostering agrobiodiversity. Increasing grain legume cultivation will be crucial to safeguard nutritional security and the resilience of agricultural ecosystems across the globe. A better understanding of the molecular underpinnings of domestication and crop evolution of grain legumes may be translated into practical approaches in modern breeding programs to stabilize yield, which is threatened by evolving pathogens and changing climates. During recent decades, domestication research in all crops has greatly benefited from the fast progress in genomic technologies. Yet still, many questions surrounding the domestication and diversification of legumes remain unanswered. In this review, we assess the potential of genomic approaches in grain legume research. We describe the centers of origin and the crucial domestication traits of grain legumes. In addition, we survey the effect of domestication on both above-ground and below-ground traits that have economic importance. Finally, we discuss open questions in grain legume domestication and diversification and outline how to bridge the gap between the preservation of historic crop diversity and their utilization in modern plant breeding.

Blue Turns to Gray: Paleogenomic Insights into the Evolutionary History and Extinction of the Blue Antelope (Hippotragus leucophaeus)

The blue antelope (Hippotragus leucophaeus) is the only large African mammal species to have become extinct in historical times, yet no nuclear genomic information is available for this species. A recent study showed that many alleged blue antelope museum specimens are either roan (Hippotragus equinus) or sable (Hippotragus niger) antelopes, further reducing the possibilities for obtaining genomic information for this extinct species. While the blue antelope has a rich fossil record from South Africa, climatic conditions in the region are generally unfavorable to the preservation of ancient DNA. Nevertheless, we recovered two blue antelope draft genomes, one at 3.4× mean coverage from a historical specimen (∼200 years old) and one at 2.1× mean coverage from a fossil specimen dating to 9,800-9,300 cal years BP, making it currently the oldest paleogenome from Africa. Phylogenomic analyses show that blue and sable antelope are sister species, confirming previous mitogenomic results, and demonstrate ancient gene flow from roan into blue antelope. We show that blue antelope genomic diversity was much lower than in roan and sable antelope, indicative of a low population size since at least the early Holocene. This supports observations from the fossil record documenting major decreases in the abundance of blue antelope after the Pleistocene-Holocene transition. Finally, the persistence of this species throughout the Holocene despite low population size suggests that colonial-era human impact was likely the decisive factor in the blue antelope's extinction.

The mitochondrial genome of the red icefish (Channichthys rugosus) casts doubt on its species status

Antarctic notothenioid fishes are recognised as one of the rare examples of adaptive radiation in the marine system. Withstanding the freezing temperatures of Antarctic waters, these fishes have diversified into over 100 species within no more than 10–20 million years. However, the exact species richness of the radiation remains contested. In the genus Channichthys, between one and nine species are recognised by different authors. To resolve the number of Channichthys species, genetic information would be highly valuable; however, so far, only sequences of a single species, C. rhinoceratus, are available. Here, we present the nearly complete sequence of the mitochondrial genome of C. rugosus, obtained from a formalin-fixed museum specimen sampled in 1974. This sequence differs from the mitochondrial genome of C. rhinoceratus in no more than 27 positions, suggesting that the two species may be synonymous.

Extended longevity of DNA preservation in Levantine Paleolithic sediments, Sefunim Cave, Israel

Paleogenomic research can elucidate the evolutionary history of human and faunal populations. Although the Levant is a key land-bridge between Africa and Eurasia, thus far, relatively little ancient DNA data has been generated from this region, since DNA degrades faster in warm climates. As sediments can be a source of ancient DNA, we analyzed 33 sediment samples from different sedimentological contexts in the Paleolithic layers of Sefunim Cave (Israel). Four contained traces of ancient Cervidae and Hyaenidae mitochondrial DNA. Dating by optical luminescence and radiocarbon indicates that the DNA comes from layers between 30,000 and 70,000 years old, surpassing theoretical expectations regarding the longevity of DNA deposited in such a warm environment. Both identified taxa are present in the zooarchaeological record of the site but have since gone extinct from the region, and a geoarchaeological study suggests little movement of the sediments after their deposition, lending further support to our findings. We provide details on the local conditions in the cave, which we hypothesize were particularly conducive to the long-term preservation of DNA-information that will be pertinent for future endeavors aimed at recovering ancient DNA from the Levant and other similarly challenging contexts.

A comparative analysis of extraction protocol performance on degraded mammalian museum specimens

The extraction of nucleic acids is one of the most routine procedures used in molecular biology laboratories, yet kit performance may influence the downstream processing of samples, particularly for samples which are degraded, and in low concentrations. Here we tested several commercial kits for specific use on commonly sampled mammalian museum specimens to evaluate the yield, size distribution, and endogenous content. Samples were weighed and had approximately equal input material for each extraction. These sample types are typical of natural history repositories ranged from 53 to 130 years old. The tested protocols spanned spin-column based extractions, magnetic bead purification, phenol/chloroform isolation, and specific modifications for ancient DNA. Diverse types of mammalian specimens were tested including adherent osteological material, bone and teeth, skin, and baleen. The concentration of DNA was quantified via fluorometry, and the size distributions of extracts visualized on an Agilent TapeStation. Overall, when DNA isolation was successful, all methods had quantifiable concentrations, albeit with variation across extracts. The length distributions varied based on the extraction protocol used. Shotgun sequencing was performed to evaluate if the extraction methods influenced the amount of endogenous versus exogenous content. The DNA content was similar across extraction methods indicating no obvious biases for DNA derived from different sources. Qiagen kits and phenol/chloroform isolation outperformed the Zymo magnetic bead isolations in these types of samples. Statistical analyses revealed that extraction method only explained 5% of the observed variation, and that specimen age explained variation (29%) more effectively.

Application of palaeogenetic techniques to historic mollusc shells reveals phylogeographic structure in a New Zealand abalone

Natural history collections worldwide contain a plethora of mollusc shells. Recent studies have detailed the sequencing of DNA extracted from shells up to thousands of years old and from various taphonomic and preservational contexts. However, previous approaches have largely addressed methodological rather than evolutionary research questions. Here we report the generation of DNA sequence data from mollusc shells using such techniques, applied to Haliotis virginea Gmelin, 1791, a New Zealand abalone, in which morphological variation has led to the recognition of several forms and subspecies. We successfully recovered near-complete mitogenomes from 22 specimens including 12 dry-preserved shells up to 60 years old. We used a combination of palaeogenetic techniques that have not previously been applied to shell, including DNA extraction optimized for ultra-short fragments and hybridization-capture of single-stranded DNA libraries. Phylogenetic analyses revealed three major, well-supported clades comprising samples from: 1) the Three Kings Islands; 2) the Auckland, Chatham and Antipodes Islands; and 3) mainland New Zealand and Campbell Island. This phylogeographic structure does not correspond to the currently recognized forms. Critically, our non-reliance on freshly collected or ethanol-preserved samples enabled inclusion of topotypes of all recognized subspecies as well as additional difficult-to-sample populations. Broader application of these comparatively cost-effective and reliable methods to modern, historical, archaeological and palaeontological shell samples has the potential to revolutionize invertebrate genetic research.

Genomic analysis of a novel Neanderthal from Mezmaiskaya Cave provides insights into the genetic relationships of Middle Palaeolithic populations

The Mezmaiskaya cave is located on the North Caucasus near the border that divides Europe and Asia. Previously, fossil remains for two Neanderthals were reported from Mezmaiskaya Cave. A tooth from the third archaic hominin specimen (Mezmaiskaya 3) was retrieved from layer 3 in Mezmaiskaya Cave. We performed genome sequencing of Mezmaiskaya 3. Analysis of partial nuclear genome sequence revealed that it belongs to a Homo sapiens neanderthalensis female. Based on a high-coverage mitochondrial genome sequence, we demonstrated that the relationships of Mezmaiskaya 3 to Mezmaiskaya 1 and Stajnia S5000 individuals were closer than those to other Neanderthals. Our data demonstrate the close genetic connections between the early Middle Palaeolithic Neanderthals that were replaced by genetically distant later group in the same geographic areas. Based on mitochondrial DNA (mtDNA) data, we suggest that Mezmaiskaya 3 was the latest Neanderthal individual from the early Neanderthal’s branches. We proposed a hierarchical nomenclature for the mtDNA haplogroups of Neanderthals. In addition, we retrieved ancestral mtDNA mutations in presumably functional sites fixed in the Neanderthal clades, and also provided the first data showing mtDNA heteroplasmy in Neanderthal specimen.

Taxonomic Identification of Two Poorly Known Lantern Shark Species Based on Mitochondrial DNA From Wet-Collection Paratypes

Etmopteridae (lantern sharks) is the most species-rich family of sharks, comprising more than 50 species. Many species are described from few individuals, and re-collection of specimens is often hindered by the remoteness of their sampling sites. For taxonomic studies, comparative morphological analysis of type specimens housed in natural history collections has been the main source of evidence. In contrast, DNA sequence information has rarely been used. Most lantern shark collection specimens, including the types, were formalin fixed before long-term storage in ethanol solutions. The DNA damage caused by both fixation and preservation of specimens has excluded these specimens from DNA sequence-based phylogenetic analyses so far. However, recent advances in the field of ancient DNA have allowed recovery of wet-collection specimen DNA sequence data. Here we analyse archival mitochondrial DNA sequences, obtained using ancient DNA approaches, of two wet-collection lantern shark paratype specimens, namely Etmopterus litvinovi and E. pycnolepis, for which the type series represent the only known individuals. Target capture of mitochondrial markers from single-stranded DNA libraries allows for phylogenetic placement of both species. Our results suggest synonymy of E. benchleyi with E. litvinovi but support the species status of E. pycnolepis. This revised taxonomy is helpful for future conservation and management efforts, as our results indicate a larger distribution range of E. litvinovi. This study further demonstrates the importance of wet-collection type specimens as genetic resource for taxonomic research.

Simultaneous Barcode Sequencing of Diverse Museum Collection Specimens Using a Mixed RNA Bait Set

A growing number of publications presenting results from sequencing natural history collection specimens reflect the importance of DNA sequence information from such samples. Ancient DNA extraction and library preparation methods in combination with target gene capture are a way of unlocking archival DNA, including from formalin-fixed wet-collection material. Here we report on an experiment, in which we used an RNA bait set containing baits from a wide taxonomic range of species for DNA hybridisation capture of nuclear and mitochondrial targets for analysing natural history collection specimens. The bait set used consists of 2,492 mitochondrial and 530 nuclear RNA baits and comprises specific barcode loci of diverse animal groups including both invertebrates and vertebrates. The baits allowed to capture DNA sequence information of target barcode loci from 84% of the 37 samples tested, with nuclear markers being captured more frequently and consensus sequences of these being more complete compared to mitochondrial markers. Samples from dry material had a higher rate of success than wet-collection specimens, although target sequence information could be captured from 50% of formalin-fixed samples. Our study illustrates how efforts to obtain barcode sequence information from natural history collection specimens may be combined and are a way of implementing barcoding inventories of scientific collection material.

Larix species range dynamics in Siberia since the Last Glacial captured from sedimentary ancient DNA

Climate change is expected to cause major shifts in boreal forests which are in vast areas of Siberia dominated by two species of the deciduous needle tree larch (Larix). The species differ markedly in their ecosystem functions, thus shifts in their respective ranges are of global relevance. However, drivers of species distribution are not well understood, in part because paleoecological data at species level are lacking. This study tracks Larix species distribution in time and space using target enrichment on sedimentary ancient DNA extracts from eight lakes across Siberia. We discovered that Larix sibirica, presently dominating in western Siberia, likely migrated to its northern distribution area only in the Holocene at around 10,000 years before present (ka BP), and had a much wider eastern distribution around 33 ka BP. Samples dated to the Last Glacial Maximum (around 21 ka BP), consistently show genotypes of L. gmelinii. Our results suggest climate as a strong determinant of species distribution in Larix and provide temporal and spatial data for species projection in a changing climate.

Using ancient sedimentary DNA from up to 50 kya, dramatic distributional shifts are documented in two dominant boreal larch species, likely guided by environmental changes suggesting climate as a strong determinant of species distribution.

Ancient genomes from the Himalayas illuminate the genetic history of Tibetans and their Tibeto-Burman speaking neighbors

Present-day Tibetans have adapted both genetically and culturally to the high altitude environment of the Tibetan Plateau, but fundamental questions about their origins remain unanswered. Recent archaeological and genetic research suggests the presence of an early population on the Plateau within the past 40 thousand years, followed by the arrival of subsequent groups within the past 10 thousand years. Here, we obtain new genome-wide data for 33 ancient individuals from high elevation sites on the southern fringe of the Tibetan Plateau in Nepal, who we show are most closely related to present-day Tibetans. They derive most of their ancestry from groups related to Late Neolithic populations at the northeastern edge of the Tibetan Plateau but also harbor a minor genetic component from a distinct and deep Paleolithic Eurasian ancestry. In contrast to their Tibetan neighbors, present-day non-Tibetan Tibeto-Burman speakers living at mid-elevations along the southern and eastern margins of the Plateau form a genetic cline that reflects a distinct genetic history. Finally, a comparison between ancient and present-day highlanders confirms ongoing positive selection of high altitude adaptive alleles.

Characteristics, origin, and potential for cancer diagnostics of ultrashort plasma cell-free DNA

Current evidence suggests that plasma cell-free DNA (cfDNA) is fragmented around a mode of 166 bp. Data supporting this view has been mainly acquired through the analysis of double-stranded cfDNA. The characteristics and diagnostic potential of single-stranded and damaged double-stranded cfDNA in healthy individuals and cancer patients remain unclear. Here, through a combination of high-affinity magnetic bead-based DNA extraction and single-stranded DNA sequencing library preparation (MB-ssDNA), we report the discovery of a large proportion of cfDNA fragments centered at ∼50 bp. We show that these "ultrashort" cfDNA fragments have a greater relative abundance in plasma of healthy individuals (median = 19.1% of all sequenced cfDNA fragments, n = 28) than in plasma of patients with cancer (median = 14.2%, n = 21, P < 0.0001). The ultrashort cfDNA fragments map to accessible chromatin regions of blood cells, particularly in promoter regions with the potential to adopt G-quadruplex (G4) DNA secondary structures. G4-positive promoter chromatin accessibility is significantly enriched in ultrashort plasma cfDNA fragments from healthy individuals relative to patients with cancers (P < 0.0001), in whom G4-cfDNA enrichment is inversely associated with copy number aberration-inferred tumor fractions. Our findings redraw the landscape of cfDNA fragmentation by identifying and characterizing a novel population of ultrashort plasma cfDNA fragments. Sequencing of MB-ssDNA libraries could facilitate the characterization of gene regulatory regions and DNA secondary structures via liquid biopsy. Our data underline the diagnostic potential of ultrashort cfDNA through classification for cancer patients.

Integrated multidisciplinary ecological analysis from the Uluzzian settlement at the Uluzzo C Rock Shelter, south-eastern Italy

The Middle to Upper Palaeolithic transition, between 50 000 and 40 000 years ago, is a period of important ecological and cultural changes. In this framework, the Rock Shelter of Uluzzo C (Apulia, southern Italy) represents an important site due to Late Mousterian and Uluzzian evidence preserved in its stratigraphic sequence. Here, we present the results of a multidisciplinary analysis performed on the materials collected between 2016 and 2018 from the Uluzzian stratigraphic units (SUs) 3, 15 and 17. The analysis involved lithic technology, use-wear, zooarchaeology, ancient DNA of sediments and palaeoproteomics, completed by quartz single-grain optically stimulated luminescence dating of the cave sediments. The lithic assemblage is characterized by a volumetric production and a debitage with no or little management of the convexities (by using the bipolar technique), with the objective to produce bladelets and flakelets. The zooarchaeological study found evidence of butchery activity and of the possible exploitation of marine resources, while drawing a picture of a patchy landscape, composed of open forests and dry open environments surrounding the shelter. Ancient mitochondrial DNA from two mammalian taxa were recovered from the sediments. Preliminary zooarchaeology by mass spectrometry results are consistent with ancient DNA and zooarchaeological taxonomic information, while further palaeoproteomics investigations are ongoing. Our new data from the re-discovery of the Uluzzo C Rock Shelter represent an important contribution to better understand the meaning of the Uluzzian in the context of the Middle/Upper Palaeolithic transition in south-eastern Italy.

Predicting Physical Appearance from DNA Data-Towards Genomic Solutions

The idea of forensic DNA intelligence is to extract from genomic data any information that can help guide the investigation. The clues to the externally visible phenotype are of particular practical importance. The high heritability of the physical phenotype suggests that genetic data can be easily predicted, but this has only become possible with less polygenic traits. The forensic community has developed DNA-based predictive tools by employing a limited number of the most important markers analysed with targeted massive parallel sequencing. The complexity of the genetics of many other appearance phenotypes requires big data coupled with sophisticated machine learning methods to develop accurate genomic predictors. A significant challenge in developing universal genomic predictive methods will be the collection of sufficiently large data sets. These should be created using whole-genome sequencing technology to enable the identification of rare DNA variants implicated in phenotype determination. It is worth noting that the correctness of the forensic sketch generated from the DNA data depends on the inclusion of an age factor. This, however, can be predicted by analysing epigenetic data. An important limitation preventing whole-genome approaches from being commonly used in forensics is the slow progress in the development and implementation of high-throughput, low DNA input sequencing technologies. The example of palaeoanthropology suggests that such methods may possibly be developed in forensics.

To bee, or not to bee? One leg is the question

Temporal genomic studies that utilise museum insects are invaluable for understanding changes in ecological processes in which insects are essential, such as wild and agricultural pollination, seed dispersal, nutrient cycling, and food web architecture, to name a few. However, given such analyses come at the cost of physical damage to museum specimens required for such work, there is a natural interest in the development and/or application of methods to minimise the damage incurred. We explored the efficacy of a recently published single stranded library construction protocol, on DNA extracted from single legs taken from eight dry-preserved historic bee specimens collected 150 years ago. Specifically, the DNA was extracted using a 'minimally destructive' method that leaves the samples' exterior intact. Our sequencing data revealed not only that the endogenous DNA recovered from some of the samples was at a relatively high level (up to 58%), but that the complexity of the libraries was sufficient in the best samples to theoretically allow deeper sequencing to a predicted level of 69x genome coverage. As such, these combined protocols offer the possibility to generate sequencing data at levels that are suitable for many common evolutionary genomic analyses, while simultaneously minimising the damage conferred to the valuable dried museum bee samples. Furthermore, we anticipate that these methods may have much wider application on many other invertebrate taxa stored in a similar way. We hope that the results from this research may be able to contribute to the increased willingness of museums to loan much needed dry-preserved insects for future genomic studies.

Social stratification without genetic differentiation at the site of Kulubnarti in Christian Period Nubia

Relatively little is known about Nubia's genetic landscape prior to the influence of the Islamic migrations that began in the late 1st millennium CE. Here, we increase the number of ancient individuals with genome-level data from the Nile Valley from three to 69, reporting data for 66 individuals from two cemeteries at the Christian Period (~650-1000 CE) site of Kulubnarti, where multiple lines of evidence suggest social stratification. The Kulubnarti Nubians had ~43% Nilotic-related ancestry (individual variation between ~36-54%) with the remaining ancestry consistent with being  introduced through Egypt and ultimately deriving from an ancestry pool like that found in the Bronze and Iron Age Levant. The Kulubnarti gene pool - shaped over a millennium - harbors disproportionately female-associated West Eurasian-related ancestry. Genetic similarity among individuals from the two cemeteries supports a hypothesis of social division without genetic distinction. Seven pairs of inter-cemetery relatives suggest fluidity between cemetery groups. Present-day Nubians are not directly descended from the Kulubnarti Nubians, attesting to additional genetic input since the Christian Period.

Genetic Evidence of Authenticity of a Hair Shaft Relic from the Portrait of Tsesarevich Alexei, Son of the Last Russian Emperor

To determine the value of an art piece, authenticity of the artwork must be verified. We demonstrate here a genetic approach to determine origin of a historic relic in the museum piece. We tested two hair shafts of unknown origin framed into a watercolor portrait of Tsesarevich Alexei Romanov, son of the last Russian Tzar Nicholas II, which is a unique item kept in the State Historical Museum. Genetic identification of the hair shafts was performed by analysis of mitochondrial DNA (mtDNA) markers using both massive parallel genomic sequencing and multiplex targeted PCR, followed by Sanger sequencing. In previous works, we reconstructed the complete mtDNA sequence inherited to Alexei Romanov through the Queen Victoria lineage [Rogaev et al. (2009) Proc. Natl. Acad. Sci. USA, 106, 5258-5263]. DNA extracts were obtained from the two thin hair shafts and used for comparative genetic analysis. Despite the very low quantity and quality of the DNA templates retrieved from the historical single hair shaft specimen, informative mtDNA sequences were determined. The mtDNA haplotype in the hair shafts corresponds to the mtDNA haplotype of Tsarevich Alexei, his sisters, and his mother, Empress Alexandra Feodorovna. This haplotype remains unique in the currently available mtDNA databases. Our results reveal that the hair relic from the portrait is associated with the family of the last Russian Emperor Nicholas II. The study is an example of first application of the genetic methodology for verification of the value of museum artwork items.

Short single-stranded DNAs with putative non-canonical structures comprise a new class of plasma cell-free DNA

Background

Cell-free DNA (cfDNA), which is extracellular DNA present in the circulating plasma and other body fluids, is currently investigated as a minimally invasive, highly informative biomarker. While nucleosome-sized cfDNA fragments have been investigated intensively, shorter DNA fragments in the plasma have not been studied due to several technical limitations.

Results

We aimed to investigate the existence of shorter cfDNA fragments in the blood. Using an improved cfDNA purification protocol and a 3′-end-labeling method, we found DNA fragments of approximately 50 nucleotides in length in the human plasma, present at a molar concentration comparable to that of nucleosome-sized fragments. Unfortunately, these short fragments cannot be recovered by widely used cfDNA isolation methods. In addition, they are composed of single-stranded DNA (ssDNA), thus escaping detection in previous studies. Therefore, we established a library-preparation protocol based on our unique ssDNA ligation technique and applied it to the isolated cfDNA. Deep sequencing of these libraries revealed that the short fragments are derived from hundreds of thousands of genomic sites in open chromatin regions and enriched with transcription factor-binding sites. Remarkably, antisense strands of putative G-quadruplex motifs occupy as much as one-third of the peaks by these short fragments.

Conclusions

We propose a new class of plasma cfDNA composed of short single-stranded fragments that potentially form non-canonical DNA structures.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01160-8.

Unlocking inaccessible historical genomes preserved in formalin

Museum specimens represent an unparalleled record of historical genomic data. However, the widespread practice of formalin preservation has thus far impeded genomic analysis of a large proportion of specimens. Limited DNA sequencing from formalin-preserved specimens has yielded low genomic coverage with unpredictable success. We set out to refine sample processing methods and to identify specimen characteristics predictive of sequencing success. With a set of taxonomically diverse specimens collected between 1962 and 2006 and ranging in preservation quality, we compared the efficacy of several end-to-end whole genome sequencing workflows alongside a k-mer-based trimming-free read alignment approach to maximize mapping of endogenous sequence. We recovered complete mitochondrial genomes and up to 3× nuclear genome coverage from formalin-preserved tissues. Hot alkaline lysis coupled with phenol-chloroform extraction out-performed proteinase K digestion in recovering DNA, while library preparation method had little impact on sequencing success. The strongest predictor of DNA yield was overall specimen condition, which additively interacts with preservation conditions to accelerate DNA degradation. Here, we demonstrate a significant advance in capability beyond limited recovery of a small number of loci via PCR or target-capture sequencing. To facilitate strategic selection of suitable specimens for genomic sequencing, we present a decision-making framework that utilizes independent and nondestructive assessment criteria. Sequencing of formalin-preserved specimens will contribute to a greater understanding of temporal trends in genetic adaptation, including those associated with a changing climate. Our work enhances the value of museum collections worldwide by unlocking genomes of specimens that have been disregarded as a valid molecular resource.

Mitochondrial DNA sequencing of a wet-collection syntype demonstrates the importance of type material as genetic resource for lantern shark taxonomy (Chondrichthyes: Etmopteridae)

After initial detection of target archival DNA of a 116-year-old syntype specimen of the smooth lantern shark, Etmopterus pusillus, in a single-stranded DNA library, we shotgun-sequenced additional 9 million reads from this same DNA library. Sequencing reads were used for extracting mitochondrial sequence information for analyses of mitochondrial DNA characteristics and reconstruction of the mitochondrial genome. The archival DNA is highly fragmented. A total of 4599 mitochondrial reads were available for the genome reconstruction using an iterative mapping approach. The resulting genome sequence has 12 times coverage and a length of 16 741 bp. All 37 vertebrate mitochondrial loci plus the control region were identified and annotated. The mitochondrial NADH2 gene was subsequently used to place the syntype haplotype in a network comprising multiple E. pusillus samples from various distant localities as well as sequences from a morphological similar species, the shortfin smooth lantern shark Etmopterus joungi. Results confirm the almost global distribution of E. pusillus and suggest E. joungi to be a junior synonym of E. pusillus. As mitochondrial DNA often represents the only available reference information in non-model organisms, this study illustrates the importance of mitochondrial DNA from an aged, wet collection type specimen for taxonomy.

Uncovering the genomic and metagenomic research potential in old ethanol-preserved snakes

Natural history museum collections worldwide represent a tremendous resource of information on past and present biodiversity. Fish, reptiles, amphibians and many invertebrate collections have often been preserved in ethanol for decades or centuries and our knowledge on the genomic and metagenomic research potential of such material is limited. Here, we use ancient DNA protocols, combined with shotgun sequencing to test the molecular preservation in liver, skin and bone tissue from five old (1842 to 1964) museum specimens of the common garter snake (Thamnophis sirtalis). When mapping reads to a T. sirtalis reference genome, we find that the DNA molecules are highly damaged with short average sequence lengths (38-64 bp) and high C-T deamination, ranging from 9% to 21% at the first position. Despite this, the samples displayed relatively high endogenous DNA content, ranging from 26% to 56%, revealing that genome-scale analyses are indeed possible from all specimens and tissues included here. Of the three tested types of tissue, bone shows marginally but significantly higher DNA quality in these metrics. Though at least one of the snakes had been exposed to formalin, neither the concentration nor the quality of the obtained DNA was affected. Lastly, we demonstrate that these specimens display a diverse and tissue-specific microbial genetic profile, thus offering authentic metagenomic data despite being submerged in ethanol for many years. Our results emphasize that historical museum collections continue to offer an invaluable source of information in the era of genomics.

Tracking the transition to agriculture in Southern Europe through ancient DNA analysis of dental calculus

Archaeological dental calculus, or mineralized plaque, is a key tool to track the evolution of oral microbiota across time in response to processes that impacted our culture and biology, such as the rise of farming during the Neolithic. However, the extent to which the human oral flora changed from prehistory until present has remained elusive due to the scarcity of data on the microbiomes of prehistoric humans. Here, we present our reconstruction of oral microbiomes via shotgun metagenomics of dental calculus in 44 ancient foragers and farmers from two regions playing a pivotal role in the spread of farming across Europe-the Balkans and the Italian Peninsula. We show that the introduction of farming in Southern Europe did not alter significantly the oral microbiomes of local forager groups, and it was in particular associated with a higher abundance of the species Olsenella sp. oral taxon 807. The human oral environment in prehistory was dominated by a microbial species, Anaerolineaceae bacterium oral taxon 439, that diversified geographically. A Near Eastern lineage of this bacterial commensal dispersed with Neolithic farmers and replaced the variant present in the local foragers. Our findings also illustrate that major taxonomic shifts in human oral microbiome composition occurred after the Neolithic and that the functional profile of modern humans evolved in recent times to develop peculiar mechanisms of antibiotic resistance that were previously absent.

Extracellular Vesicle-Derived DNA vs. CfDNA as a Biomarker for the Detection of Colon Cancer

Liquid biopsy has emerged as a promising non-invasive way to diagnose tumor and monitor its progression. Different types of liquid biopsies have different advantages and limitations. In the present research, we compared the use of two types of liquid biopsy, extracellular vesicle-derived DNA (EV-DNA) and cell-free DNA (cfDNA) for identifying tumor mutations in patients with colon carcinoma. Method: DNA was extracted from the tumor tissue of 33 patients diagnosed with colon carcinoma. Targeted NGS panel, based on the hotspots panel, was used to identify tumor mutations. Pre-surgery serum and plasma were taken from the patients in which mutation was found in the tumor tissue. Extracellular vesicles were isolated from the serum followed by the extraction of EV-DNA. CfDNA was extracted from the plasma. The mutations found in the tumor were used to detect the circulating tumor DNA using ultra-deep sequencing. We compared the sensitivity of mutation detection and allele frequency obtained in EV-DNA and cfDNA. Results: The sensitivity of mutation detection in EV-DNA and cfDNA was 61.90% and 66.67%, respectively. We obtained almost identical sensitivity of mutation detection in EV-DNA and cfDNA in each of the four stages of colon carcinoma. The total DNA concentration and number mutant copies were higher in cfDNA vs. EV-DNA (p value = 0.002 and 0.003, respectively). Conclusion: Both cfDNA and EV-DNA can serve as tumor biomarkers. The use of EV-DNA did not lead to improved sensitivity or better detection of tumor DNA in the circulation.

A Fast and Efficient Single-stranded Genomic Library Preparation Method Optimized for Ancient DNA

We present a protocol to prepare extracted DNA for sequencing on the Illumina sequencing platform that has been optimized for ancient and degraded DNA. Our approach, the Santa Cruz Reaction or SCR, uses directional splinted ligation of Illumina’s P5 and P7 adapters to convert natively single-stranded DNA and heat denatured double-stranded DNA into sequencing libraries in a single enzymatic reaction. To demonstrate its efficacy in converting degraded DNA molecules, we prepare 5 ancient DNA extracts into sequencing libraries using the SCR and 2 of the most commonly used approaches for preparing degraded DNA for sequencing: BEST, which targets and converts double-stranded DNA, and ssDNA2.0, which targets and converts single-stranded DNA. We then compare the efficiency with which each approach recovers unique molecules, or library complexity, given a standard amount of DNA input. We find that the SCR consistently outperforms the BEST protocol in recovering unique molecules and, despite its relative simplicity to perform and low cost per library, has similar performance to ssDNA2.0 across a wide range of DNA inputs. The SCR is a cost- and time-efficient approach that minimizes the loss of unique molecules and makes accessible a taxonomically, geographically, and a temporally broader sample of preserved remains for genomic analysis.

Identifying the true number of specimens of the extinct blue antelope (Hippotragus leucophaeus)

Native to southern Africa, the blue antelope (Hippotragus leucophaeus) is the only large African mammal species known to have become extinct in historical times. However, it was poorly documented prior to its extinction ~ 1800 AD, and many of the small number of museum specimens attributed to it are taxonomically contentious. This places limitations on our understanding of its morphology, ecology, and the mechanisms responsible for its demise. We retrieved genetic information from ten of the sixteen putative blue antelope museum specimens using both shotgun sequencing and mitochondrial genome target capture in an attempt to resolve the uncertainty surrounding the identification of these specimens. We found that only four of the ten investigated specimens, and not a single skull, represent the blue antelope. This indicates that the true number of historical museum specimens of the blue antelope is even smaller than previously thought, and therefore hardly any reference material is available for morphometric, comparative and genetic studies. Our study highlights how genetics can be used to identify rare species in natural history collections where other methods may fail or when records are scarce. Additionally, we present an improved mitochondrial reference genome for the blue antelope as well as one complete and two partial mitochondrial genomes. A first analysis of these mitochondrial genomes indicates low levels of maternal genetic diversity in the ‘museum population’, possibly confirming previous results that blue antelope population size was already low at the time of the European colonization of South Africa.

The evolutionary history of Neanderthal and Denisovan Y chromosomes

Ancient DNA has provided new insights into many aspects of human history. However, we lack comprehensive studies of the Y chromosomes of Denisovans and Neanderthals because the majority of specimens that have been sequenced to sufficient coverage are female. Sequencing Y chromosomes from two Denisovans and three Neanderthals shows that the Y chromosomes of Denisovans split around 700 thousand years ago from a lineage shared by Neanderthals and modern human Y chromosomes, which diverged from each other around 370 thousand years ago. The phylogenetic relationships of archaic and modern human Y chromosomes differ from the population relationships inferred from the autosomal genomes and mirror mitochondrial DNA phylogenies, indicating replacement of both the mitochondrial and Y chromosomal gene pools in late Neanderthals. This replacement is plausible if the low effective population size of Neanderthals resulted in an increased genetic load in Neanderthals relative to modern humans.

Moose genomes reveal past glacial demography and the origin of modern lineages

Background

Numerous megafauna species from northern latitudes went extinct during the Pleistocene/Holocene transition as a result of climate-induced habitat changes. However, several ungulate species managed to successfully track their habitats during this period to eventually flourish and recolonise the holarctic regions. So far, the genomic impacts of these climate fluctuations on ungulates from high latitudes have been little explored. Here, we assemble a de-novo genome for the European moose (Alces alces) and analyse it together with re-sequenced nuclear genomes and ancient and modern mitogenomes from across the moose range in Eurasia and North America.

Results

We found that moose demographic history was greatly influenced by glacial cycles, with demographic responses to the Pleistocene/Holocene transition similar to other temperate ungulates. Our results further support that modern moose lineages trace their origin back to populations that inhabited distinct glacial refugia during the Last Glacial Maximum (LGM). Finally, we found that present day moose in Europe and North America show low to moderate inbreeding levels resulting from post-glacial bottlenecks and founder effects, but no evidence for recent inbreeding resulting from human-induced population declines.

Conclusions

Taken together, our results highlight the dynamic recent evolutionary history of the moose and provide an important resource for further genomic studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07208-3.