Search

DNA obtained from environmental samples such as sediments, ice or water (environmental DNA, eDNA), represents an important source of information on past and present biodiversity. It has revealed an ancient forest in Greenland, extended by several thousand years the survival dates for mainland woolly mammoth in Alaska, and pushed back the dates for spruce survival in Scandinavian ice-free refugia during the last glaciation. More recently, eDNA was used to uncover the past 50 000 years of vegetation history in the Arctic, revealing massive vegetation turnover at the Pleistocene/Holocene transition, with implications for the extinction of megafauna. Furthermore, eDNA can reflect the biodiversity of extant flora and fauna, both qualitatively and quantitatively, allowing detection of rare species. As such, trace studies of plant and vertebrate DNAin the environment have revolutionized our knowledge of biogeography. However, the approach remains marred by biases related toDNAbehaviour in environmental settings, incomplete reference databases and false positive results due to contamination.We provide a review of the field. [ABSTRACT FROM AUTHOR]

Summary: Ancient DNA has emerged as a remarkable tool to infer the history of extinct species and past populations. However, many of its characteristics, such as extensive fragmentation, damage and contamination, can influence downstream analyses. To help investigators measure how these could impact their analyses in silico, we have developed gargammel, a package that simulates ancient DNA fragments given a set of known reference genomes. Our package simulates the entire molecular process from post-mortem DNA fragmentation and DNA damage to experimental sequencing errors, and reproduces most common bias observed in ancient DNA datasets. [ABSTRACT FROM AUTHOR]

The origins of the First Americans remain contentious. Although Native Americans seem to be genetically most closely related to east Asians, there is no consensus with regard to which specific Old World populations they are closest to. Here we sequence the draft genome of an approximately 24,000-year-old individual (MA-1), from Mal'ta in south-central Siberia, to an average depth of 1×. To our knowledge this is the oldest anatomically modern human genome reported to date. The MA-1 mitochondrial genome belongs to haplogroup U, which has also been found at high frequency among Upper Palaeolithic and Mesolithic European hunter-gatherers, and the Y chromosome of MA-1 is basal to modern-day western Eurasians and near the root of most Native American lineages. Similarly, we find autosomal evidence that MA-1 is basal to modern-day western Eurasians and genetically closely related to modern-day Native Americans, with no close affinity to east Asians. This suggests that populations related to contemporary western Eurasians had a more north-easterly distribution 24,000 years ago than commonly thought. Furthermore, we estimate that 14 to 38% of Native American ancestry may originate through gene flow from this ancient population. This is likely to have occurred after the divergence of Native American ancestors from east Asian ancestors, but before the diversification of Native American populations in the New World. Gene flow from the MA-1 lineage into Native American ancestors could explain why several crania from the First Americans have been reported as bearing morphological characteristics that do not resemble those of east Asians. Sequencing of another south-central Siberian, Afontova Gora-2 dating to approximately 17,000 years ago, revealed similar autosomal genetic signatures as MA-1, suggesting that the region was continuously occupied by humans throughout the Last Glacial Maximum. Our findings reveal that western Eurasian genetic signatures in modern-day Native Americans derive not only from post-Columbian admixture, as commonly thought, but also from a mixed ancestry of the First Americans. [ABSTRACT FROM AUTHOR]

Freshwater ecosystems are among the most endangered habitats on Earth, with thousands of animal species known to be threatened or already extinct. Reliable monitoring of threatened organisms is crucial for data-driven conservation actions but remains a challenge owing to nonstandardized methods that depend on practical and taxonomic expertise, which is rapidly declining. Here, we show that a diversity of rare and threatened freshwater animals-representing amphibians, fish, mammals, insects and crustaceans-can be detected and quantified based on DNA obtained directly from small water samples of lakes, ponds and streams. We successfully validate our findings in a controlled mesocosm experiment and show that DNA becomes undetectable within 2 weeks after removal of animals, indicating that DNA traces are near contemporary with presence of the species. We further demonstrate that entire faunas of amphibians and fish can be detected by high-throughput sequencing of DNA extracted from pond water. Our findings underpin the ubiquitous nature of DNA traces in the environment and establish environmental DNA as a tool for monitoring rare and threatened species across a wide range of taxonomic groups. [ABSTRACT FROM AUTHOR]

Nunataks are isolated bedrocks protruding through ice sheets. They vary in age, but represent island environments in 'oceans' of ice through which organism dispersals and replacements can be studied over time. The J.A.D. Jensen's Nunataks at the southern Greenland ice sheet are the most isolated nunataks on the northern hemisphere - some 30 km from the nearest biological source. They constitute around 2 km2 of ice-free land that was established in the early Holocene. We have investigated the changes in plant composition at these nunataks using both the results of surveys of the flora over the last 130 years and through reconstruction of the vegetation from the end of the Holocene Thermal Maximum (5528 ± 75 cal year BP) using meta-barcoding of plant DNA recovered from the nunatak sediments ( sedaDNA). Our results show that several of the plant species detected with sedaDNA are described from earlier vegetation surveys on the nunataks (in 1878, 1967 and 2009). In 1967, a much higher biodiversity was detected than from any other of the studied periods. While this may be related to differences in sampling efforts for the oldest period, it is not the case when comparing the 1967 and 2009 levels where the botanical survey was exhaustive. As no animals and humans are found on the nunataks, this change in diversity over a period of just 42 years must relate to environmental changes probably being climate-driven. This suggests that even the flora of fairly small and isolated ice-free areas reacts quickly to a changing climate. [ABSTRACT FROM AUTHOR]

DNA molecules originating from animals and plants can be retrieved directly from sediments and have been used for reconstructing both contemporary and past ecosystems. However, the extent to which such 'dirt' DNA reflects taxonomic richness and structural diversity remains contentious. Here, we couple second generation high-throughput sequencing with 16S mitochondrial DNA (mtDNA) meta-barcoding, to explore the accuracy and sensitivity of 'dirt' DNA as an indicator of vertebrate diversity, from soil sampled at safari parks, zoological gardens and farms with known species compositions. PCR amplification was successful in the full pH range of the investigated soils (6.2 ± 0.2 to 8.3 ± 0.2), but inhibition was detected in extracts from soil of high organic content. DNA movement (leaching) through strata was evident in some sporadic cases and is influenced by soil texture and structure. We find that DNA from the soil surface reflects overall taxonomic richness and relative biomass of individual species. However, one species that was recently introduced was not detected. Furthermore, animal behaviour was shown to influence DNA deposition rates. The approach potentially provides a quick methodological alternative to classical ecological surveys of biodiversity, and most reliable results are obtained with spatial sample replicates, while relative amounts of soil processed per site is of less importance. [ABSTRACT FROM AUTHOR]

Although ancient DNA from sediments ( sedaDNA) has been used to investigate past ecosystems, the approach has never been directly compared with the traditional methods of pollen and macrofossil analysis. We conducted a comparative survey of 18 ancient permafrost samples spanning the Late Pleistocene (46-12.5 thousand years ago), from the Taymyr Peninsula in northern Siberia. The results show that pollen, macrofossils and sedaDNA are complementary rather than overlapping and, in combination, reveal more detailed information on plant palaeocommunities than can be achieved by each individual approach. SedaDNA and macrofossils share greater overlap in plant identifications than with pollen, suggesting that sedaDNA is local in origin. These two proxies also permit identification to lower taxonomic levels than pollen, enabling investigation into temporal changes in species composition and the determination of indicator species to describe environmental changes. Combining data from all three proxies reveals an area continually dominated by a mosaic vegetation of tundra-steppe, pioneer and wet-indicator plants. Such vegetational stability is unexpected, given the severe climate changes taking place in the Northern Hemisphere during this time, with changes in average annual temperatures of >22 °C. This may explain the abundance of ice-age mammals such as horse and bison in Taymyr Peninsula during the Pleistocene and why it acted as a refugium for the last mainland woolly mammoth. Our finding reveals the benefits of combining sedaDNA, pollen and macrofossil for palaeovegetational reconstruction and adds to the increasing evidence suggesting large areas of the Northern Hemisphere remained ecologically stable during the Late Pleistocene. [ABSTRACT FROM AUTHOR]

Motivation: Large-scale genotype datasets can help track the dispersal patterns of epidemiological outbreaks and predict the geographic origins of individuals. Such genetically-based geographic assignments also show a range of possible applications in forensics for profiling both victims and criminals, and in wildlife management, where poaching hotspot areas can be located. They, however, require fast and accurate statistical methods to handle the growing amount of genetic information made available from genotype arrays and next-generation sequencing technologies. Results: We introduce a novel statistical method for geopositioning individuals of unknown origin from genotypes. Our method is based on a geostatistical model trained with a dataset of georeferenced genotypes. Statistical inference under this model can be implemented within the theoretical framework of Integrated Nested Laplace Approximation, which represents one of the major recent breakthroughs in statistics, as it does not require Monte Carlo simulations. We compare the performance of our method and an alternative method for geospatial inference, SPA in a simulation framework. We highlight the accuracy and limits of continuous spatial assignment methods at various scales by analyzing genotype datasets from a diversity of species, including Florida Scrubjay birds Aphelocoma coerulescens, Arabidopsis thaliana and humans, representing 41-197,146 SNPs. Our method appears to be best suited for the analysis of medium-sized datasets (a few tens of thousands of loci), such as reduced-representation sequencing data that become increasingly available in ecology. Availability and implementation: http://www2.imm.dtu.dk∼gigu/Spasiba/. [ABSTRACT FROM AUTHOR]

Motivation The alignment of sequencing reads is a critical step in the characterization of ancient genomes. However, reference bias and spurious mappings pose a significant challenge, particularly as cutting-edge wet lab methods generate datasets that push the boundaries of alignment tools. Reference bias occurs when reference alleles are favoured over alternative alleles during mapping, whereas spurious mappings stem from either contamination or when endogenous reads fail to align to their correct position. Previous work has shown that these phenomena are correlated with read length but a more thorough investigation of reference bias and spurious mappings for ancient DNA has been lacking. Here, we use a range of empirical and simulated palaeogenomic datasets to investigate the impacts of mapping tools, quality thresholds, and reference genome on mismatch rates across read lengths. Results For these analyses, we introduce AMBER, a new bioinformatics tool for assessing the quality of ancient DNA mapping directly from BAM-files and informing on reference bias, read length cut-offs and reference selection. AMBER rapidly and simultaneously computes the sequence read mapping bias in the form of the mismatch rates per read length, cytosine deamination profiles at both CpG and non-CpG sites, fragment length distributions, and genomic breadth and depth of coverage. Using AMBER, we find that mapping algorithms and quality threshold choices dictate reference bias and rates of spurious alignment at different read lengths in a predictable manner, suggesting that optimized mapping parameters for each read length will be a key step in alleviating reference bias and spurious mappings. Availability and implementation AMBER is available for noncommercial use on GitHub (https://github.com/tvandervalk/AMBER.git). Scripts used to generate and analyse simulated datasets are available on Github (https://github.com/sdolenz/refbias_scripts). [ABSTRACT FROM AUTHOR]

Background: The mammalian Major Histocompatibility Complex (MHC) is a genetic region containing highly polymorphic genes with immunological functions. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. The MHC class II sub-region contains genes expressed in antigen presenting cells. The antigen binding site is encoded by the second exon of genes encoding antigen presenting molecules. The exon 2 sequences of these MHC genes have evolved under the selective pressure of pathogens. Interspecific differences can be observed in the class II sub-region. The family Equidae includes a variety of domesticated, and free-ranging species inhabiting a range of habitats exposed to different pathogens and represents a model for studying this important part of the immunogenome. While equine MHC class II DRA and DQA loci have received attention, the genetic diversity and effects of selection on DRB and DQB loci have been largely overlooked. This study aimed to provide the first in-depth analysis of the MHC class II DRB and DQB loci in the Equidae family. Results: Three DRB and two DQB genes were identified in the genomes of all equids. The genes DRB2, DRB3 and DQB3 showed high sequence conservation, while polymorphisms were more frequent at DRB1 and DQB1 across all species analyzed. DQB2 was not found in the genome of the Asiatic asses Equus hemionus kulan and E. h. onager. The bioinformatic analysis of non-zero-coverage-bases of DRB and DQB genes in 14 equine individual genomes revealed differences among individual genes. Evidence for recombination was found for DRB1, DRB2, DQB1 and DQB2 genes. Trans-species allele sharing was identified in all genes except DRB1. Site-specific selection analysis predicted genes evolving under positive selection both at DRB and DQB loci. No selected amino acid sites were identified in DQB3. Conclusions: The organization of the MHC class II sub-region of equids is similar across all species of the family. Genomic sequences, along with phylogenetic trees suggesting effects of selection as well as trans-species polymorphism support the contention that pathogen-driven positive selection has shaped the MHC class II DRB/DQB sub-regions in the Equidae. [ABSTRACT FROM AUTHOR]

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

The horse was essential to past human societies but became a recreational animal during the twentieth century as the world became increasingly mechanized. As the author reviews here, recent studies of ancient genomes have revisited the understanding of horse domestication, from the very early stages to the most modern developments. They have uncovered several extinct lineages roaming the far ends of Eurasia some 4000 years ago. They have shown that the domestic horse has been significantly reshaped during the last millennium and experienced a sharp decline in genetic diversity within the last two centuries. At a time when no truly wild horses exist any longer, this calls for enhanced conservation in all endangered populations. These include the Przewalski's horse native to Mongolia, and the many local breeds side‐lined by the modern agenda, but yet representing the living heritage of over five millennia of horse breeding. [ABSTRACT FROM AUTHOR]

The control of transcription is crucial for homeostasis in mammals. A previous selective sweep analysis of horse racing performance revealed a 19.6 kb candidate regulatory region 50 kb downstream of the Endothelin3 (EDN3) gene. Here, the region was narrowed to a 5.5 kb span of 14 SNVs, with elite and sub-elite haplotypes analyzed for association to racing performance, blood pressure and plasma levels of EDN3 in Coldblooded trotters and Standardbreds. Comparative analysis of human HiCap data identified the span as an enhancer cluster active in endothelial cells, interacting with genes relevant to blood pressure regulation. Coldblooded trotters with the sub-elite haplotype had significantly higher blood pressure compared to horses with the elite performing haplotype during exercise. Alleles within the elite haplotype were part of the standing variation in pre-domestication horses, and have risen in frequency during the era of breed development and selection. These results advance our understanding of the molecular genetics of athletic performance and vascular traits in both horses and humans. Author summary: A previous study discovered that a genomic region close to the Endothelin3 gene was associated with harness racing performance. Here, careful phenotypic documentation of athletic performance and blood pressure measurements in horses, followed by state-of-the-art genomics, allowed us to identify a 5.5 kb regulatory region located approximately 50 kb 3' of the EDN3 gene. A comparative analysis of the region using human HiCap data supported a regulatory role as, in endothelial cells, interaction was observed between the region and multiple genes relevant to blood pressure regulation and athletic performance. Long range cis-regulatory modules are critical for cooperatively controlling multiple genes located within transcriptionally active domains. We measured blood pressure in Coldblooded trotters during exercise and demonstrated that horses with two copies of the elite-performing haplotype had lower blood pressure during exercise and better racing performance results, compared to horses with two copies of the sub-elite performing haplotype. In addition, horses with the elite-performing haplotype also had higher levels of Endothelin3 in plasma. The results reported here are important for understanding the biological mechanisms behind blood pressure regulation in relation to racing performance in both horses and humans. [ABSTRACT FROM AUTHOR]

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

Recent advances in genomic sequencing technology and computational assembly methods have allowed scientists to improve reference genome assemblies in terms of contiguity and composition. EquCab2, a reference genome for the domestic horse, was released in 2007. Although of equal or better quality compared to other first-generation Sanger assemblies, it had many of the shortcomings common to them. In 2014, the equine genomics research community began a project to improve the reference sequence for the horse, building upon the solid foundation of EquCab2 and incorporating new short-read data, long-read data, and proximity ligation data. Here, we present EquCab3. The count of non-N bases in the incorporated chromosomes is improved from 2.33 Gb in EquCab2 to 2.41 Gb in EquCab3. Contiguity has also been improved nearly 40-fold with a contig N50 of 4.5 Mb and scaffold contiguity enhanced to where all but one of the 32 chromosomes is comprised of a single scaffold. Theodore Kalbfleisch et al. present an improved genome assembly for the domestic horse by combining short- and long-read data, as well as proximity ligation data. They improve contiguity of the assembly by 40-fold, with a 10-fold reduction in gaps. [ABSTRACT FROM AUTHOR]

This paper reports a high-resolution isotopic study of medieval horse mobility, revealing their origins and in-life mobility both regionally and internationally. The animals were found in an unusual horse cemetery site found within the City of Westminster, London, England. Enamel strontium, oxygen, and carbon isotope analysis of 15 individuals provides information about likely place of birth, diet, and mobility during the first approximately 5 years of life. Results show that at least seven horses originated outside of Britain in relatively cold climates, potentially in Scandinavia or the Western Alps. Ancient DNA sexing data indicate no consistent sex-specific mobility patterning, although three of the five females came from exceptionally highly radiogenic regions. Another female with low mobility is suggested to be a sedentary broodmare. Our results provide direct and unprecedented evidence for a variety of horse movement and trading practices in the Middle Ages and highlight the importance of international trade in securing high-quality horses for medieval London elites. [ABSTRACT FROM AUTHOR]

Whole genome characterizations of crop plants based on ancient DNA have provided unique keys for a better understanding of the evolutionary origins of modern cultivars, the pace and mode of selection underlying their adaptation to new environments and the production of phenotypes of interest. Although forests are among the most biologically rich ecosystems on earth and represent a fundamental resource for human societies, no ancient genome sequences have been generated for trees. This contrasts with the generation of multiple ancient reference genomes for important crops. Here, we sequenced the first ancient tree genomes using two white oak wood remains from Germany dating to the Last Little Ice Age (15th century CE, 7.3× and 4.0×) and one from France dating to the Bronze Age (1700 BCE, 3.4×). We assessed the underlying species and identified one medieval remains as a hybrid between two common oak species (Quercus robur and Q. petraea) and the other two remains as Q. robur. We found that diversity at the global genome level had not changed over time. However, exploratory analyses suggested that a reduction of diversity took place at different time periods. Finally, we determined the timing of leaf unfolding for ancient trees for the first time. The study extends the application of ancient wood beyond the classical proxies of dendroclimatology, dendrochronology, dendroarchaeology and dendroecology, thereby enhancing resolution of inferences on the responses of forest ecosystems to past environmental changes, epidemics and silvicultural practices. see also the Perspective by Andrew L. Hipp and Desanka Lazic [ABSTRACT FROM AUTHOR]

Dynamic evolutionary processes and complex structure make the Y chromosome among the most diverse and least understood regions in mammalian genomes. Here, we present an annotated assembly of the male specific region of the horse Y chromosome (eMSY), representing the first comprehensive Y assembly in odd-toed ungulates. The eMSY comprises single-copy, equine specific multi-copy, PAR transposed, and novel ampliconic sequence classes. The eMSY gene density approaches that of autosomes with the highest number of retained X-Y gametologs recorded in eutherians, in addition to novel Y-born and transposed genes. Horse, donkey and mule testis RNAseq reveals several candidate genes for stallion fertility. A novel testis-expressed XY ampliconic sequence class, ETSTY7, is shared with the parasite Parascaris genome, providing evidence for eukaryotic horizontal transfer and inter-chromosomal mobility. Our study highlights the dynamic nature of the Y and provides a reference sequence for improved understanding of equine male development and fertility. The rapidly evolving Y chromosome accumulates male-benefit genes but is often poorly characterized in many mammals. Here, the authors assemble the male specific region of the horse Y chromosome and investigate its evolution and function. [ABSTRACT FROM AUTHOR]

Horses have been valued for their diversity of coat colour since prehistoric times; this is especially the case since their domestication in the Caspian steppe in ~3,500 BC. Although we can assume that human preferences were not constant, we have only anecdotal information about how domestic horses were influenced by humans. Our results from genotype analyses show a significant increase in spotted coats in early domestic horses (Copper Age to Iron Age). In contrast, medieval horses carried significantly fewer alleles for these phenotypes, whereas solid phenotypes (i.e., chestnut) became dominant. This shift may have been supported because of (i) pleiotropic disadvantages, (ii) a reduced need to separate domestic horses from their wild counterparts, (iii) a lower religious prestige, or (iv) novel developments in weaponry. These scenarios may have acted alone or in combination. However, the dominance of chestnut is a remarkable feature of the medieval horse population. [ABSTRACT FROM AUTHOR]