Search

Human intervention, pre-human climate change (or a combination of both), as well as genetic effects, contribute to species extinctions. While many species from oceanic islands have gone extinct due to direct human impacts, the effects of pre-human climate change and human settlement on the genomic diversity of insular species and the role that loss of genomic diversity played in their extinctions remains largely unexplored. To address this question, we sequenced whole genomes of two extinct New Zealand passerines, the huia (Heteralocha acutirostris) and South Island kōkako (Callaeas cinereus). Both species showed similar demographic trajectories throughout the Pleistocene. However, the South Island kōkako continued to decline after the last glaciation, while the huia experienced some recovery. Moreover, there was no indication of inbreeding resulting from recent mating among closely related individuals in either species. This latter result indicates that population fragmentation associated with forest clearing by Maōri may not have been strong enough to lead to an increase in inbreeding and exposure to genomic erosion. While genomic erosion may not have directly contributed to their extinctions, further habitat fragmentation and the introduction of mammalian predators by Europeans may have been an important driver of extinction in huia and South Island kōkako.

Arid environments provide ideal ground for investigating the mechanisms of adaptive evolution. High temperatures and low water availability are relentless stressors for many endotherms, including birds; yet birds persist in deserts. While physiological adaptation probably involves metabolic phenotypes, the underlying mechanisms (plasticity, genetics) are largely uncharacterized. To explore this, we took an intraspecific approach that focused on a species that is resident over a mesic to arid gradient, the Karoo scrub-robin (Cercotrichas coryphaeus). Specifically, we integrated environmental (climatic and primary productivity), physiological (metabolic rates: a measure of energy expenditure), genotypic (genetic variation underlying the machinery of energy production) and microbiome (involved in processing food from where energy is retrieved) data, to infer the mechanism of physiological adaptation. We that found the variation in energetic physiology phenotypes and gut microbiome composition are associated with environmental features as well as with variation in genes underlying energy metabolic pathways. Specifically, we identified a small list of candidate adaptive genes, some of them with known ties to relevant physiology phenotypes. Together our results suggest that selective pressures on energetic physiology mediated by genes related to energy homeostasis and possibly microbiota composition may facilitate adaptation to local conditions and provide an explanation to the high avian intraspecific divergence observed in harsh environments., (© 2019 John Wiley & Sons Ltd.)

The current phylogeographic pattern of European brown bears ( Ursus arctos ) has commonly been explained by postglacial recolonization out of geographically distinct refugia in southern Europe, a pattern well in accordance with the expansion/contraction model. Studies of ancient DNA from brown bear remains have questioned this pattern, but have failed to explain the glacial distribution of mitochondrial brown bear clades and their subsequent expansion across the European continent. We here present 136 new mitochondrial sequences generated from 346 remains from Europe, ranging in age between the Late Pleistocene and historical times. The genetic data show a high Late Pleistocene diversity across the continent and challenge the strict confinement of bears to traditional southern refugia during the last glacial maximum (LGM). The mitochondrial data further suggest a genetic turnover just before this time, as well as a steep demographic decline starting in the mid-Holocene. Levels of stable nitrogen isotopes from the remains confirm a previously proposed shift toward increasing herbivory around the LGM in Europe. Overall, these results suggest that in addition to climate, anthropogenic impact and inter-specific competition may have had more important effects on the brown bear's ecology, demography, and genetic structure than previously thought., Competing Interests: None declared.

Environmental gradients are caused by gradual changes in abiotic factors, which affect species abundances and distributions, and are important for the spatial distribution of biodiversity. One prominent environmental gradient is the altitude gradient. Understanding ecological processes associated with altitude gradients may help us to understand the possible effects climate change could have on species communities. We quantified vegetation cover, species richness, species evenness, beta diversity, and spatial patterns of community structure of vascular plants along altitude gradients in a subarctic mountain tundra in northern Sweden. Vascular plant cover and plant species richness showed unimodal relationships with altitude. However, species evenness did not change with altitude, suggesting that no individual species became dominant when species richness declined. Beta diversity also showed a unimodal relationship with altitude, but only for an intermediate spatial scale of 1 km. A lack of relationships with altitude for either patch or landscape scales suggests that any altitude effects on plant spatial heterogeneity occurred on scales larger than individual patches but were not effective across the whole landscape. We observed both nested and modular patterns of community structures, but only the modular patterns corresponded with altitude. Our observations point to biotic regulations of plant communities at high altitudes, but we found both scale dependencies and inconsistent magnitude of the effects of altitude on different diversity components. We urge for further studies evaluating how different factors influence plant communities in high altitude and high latitude environments, as well as studies identifying scale and context dependencies in any such influences., Competing Interests: None declared.

Many endangered species have experienced severe population declines within the last centuries [1, 2]. However, despite concerns about negative fitness effects resulting from increased genetic drift and inbreeding, there is a lack of empirical data on genomic changes in conjunction with such declines [3-7]. Here, we use whole genomes recovered from century-old historical museum specimens to quantify the genomic consequences of small population size in the critically endangered Grauer's and endangered mountain gorillas. We find a reduction of genetic diversity and increase in inbreeding and genetic load in the Grauer's gorilla, which experienced severe population declines in recent decades. In contrast, the small but relatively stable mountain gorilla population has experienced little genomic change during the last century. These results suggest that species histories as well as the rate of demographic change may influence how population declines affect genome diversity., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

North America is currently home to a number of grey wolf (Canis lupus) and wolf-like canid populations, including the coyote (Canis latrans) and the taxonomically controversial red, Eastern timber and Great Lakes wolves. We explored their population structure and regional gene flow using a dataset of 40 full genome sequences that represent the extant diversity of North American wolves and wolf-like canid populations. This included 15 new genomes (13 North American grey wolves, 1 red wolf and 1 Eastern timber/Great Lakes wolf), ranging from 0.4 to 15x coverage. In addition to providing full genome support for the previously proposed coyote-wolf admixture origin for the taxonomically controversial red, Eastern timber and Great Lakes wolves, the discriminatory power offered by our dataset suggests all North American grey wolves, including the Mexican form, are monophyletic, and thus share a common ancestor to the exclusion of all other wolves. Furthermore, we identify three distinct populations in the high arctic, one being a previously unidentified "Polar wolf" population endemic to Ellesmere Island and Greenland. Genetic diversity analyses reveal particularly high inbreeding and low heterozygosity in these Polar wolves, consistent with long-term isolation from the other North American wolves., Competing Interests: The authors have declared that no competing interests exist.

For millennia, the Pontic-Caspian steppe was a connector between the Eurasian steppe and Europe. In this scene, multidirectional and sequential movements of different populations may have occurred, including those of the Eurasian steppe nomads. We sequenced 35 genomes (low to medium coverage) of Bronze Age individuals (Srubnaya-Alakulskaya) and Iron Age nomads (Cimmerians, Scythians, and Sarmatians) that represent four distinct cultural entities corresponding to the chronological sequence of cultural complexes in the region. Our results suggest that, despite genetic links among these peoples, no group can be considered a direct ancestor of the subsequent group. The nomadic populations were heterogeneous and carried genetic affinities with populations from several other regions including the Far East and the southern Urals. We found evidence of a stable shared genetic signature, making the eastern Pontic-Caspian steppe a likely source of western nomadic groups.

Although many large mammal species went extinct at the end of the Pleistocene epoch, their DNA may persist due to past episodes of interspecies admixture. However, direct empirical evidence of the persistence of ancient alleles remains scarce. Here, we present multifold coverage genomic data from four Late Pleistocene cave bears (Ursus spelaeus complex) and show that cave bears hybridized with brown bears (Ursus arctos) during the Pleistocene. We develop an approach to assess both the directionality and relative timing of gene flow. We find that segments of cave bear DNA still persist in the genomes of living brown bears, with cave bears contributing 0.9 to 2.4% of the genomes of all brown bears investigated. Our results show that even though extinction is typically considered as absolute, following admixture, fragments of the gene pool of extinct species can survive for tens of thousands of years in the genomes of extant recipient species.

Hyenas (family Hyaenidae), as the sister group to cats (family Felidae), represent a deeply diverging branch within the cat-like carnivores (Feliformia). With an estimated population size of <10,000 individuals worldwide, the brown hyena (Parahyaena brunnea) represents the rarest of the four extant hyena species and has been listed as Near Threatened by the IUCN. Here, we report a high-coverage genome from a captive bred brown hyena and both mitochondrial and low-coverage nuclear genomes of 14 wild-caught brown hyena individuals from across southern Africa. We find that brown hyena harbor extremely low genetic diversity on both the mitochondrial and nuclear level, most likely resulting from a continuous and ongoing decline in effective population size that started ∼1 Ma and dramatically accelerated towards the end of the Pleistocene. Despite the strikingly low genetic diversity, we find no evidence of inbreeding within the captive bred individual and reveal phylogeographic structure, suggesting the existence of several potential subpopulations within the species.

Species and populations are disappearing at an alarming rate as a direct result of human activities. Loss of genetic diversity associated with population decline directly impacts species' long-term survival. Therefore, preserving genetic diversity is of considerable conservation importance. However, to assist in conservation efforts, it is important to understand how genetic diversity is spatially distributed and how it changes due to anthropogenic pressures. In this study, we use historical museum and modern faecal samples of two critically endangered eastern gorilla taxa, Grauer's (Gorilla beringei graueri) and mountain gorillas (Gorilla beringei beringei), to directly infer temporal changes in genetic diversity within the last century. Using over 100 complete mitochondrial genomes, we observe a significant decline in haplotype and nucleotide diversity in Grauer's gorillas. By including historical samples from now extinct populations we show that this decline can be attributed to the loss of peripheral populations rather than a decrease in genetic diversity within the core range of the species. By directly quantifying genetic changes in the recent past, our study shows that human activities have severely impacted eastern gorilla genetic diversity within only four to five generations. This rapid loss calls for dedicated conservation actions, which should include preservation of the remaining peripheral populations.

Understanding how species respond to population declines is a central question in conservation and evolutionary biology. Population declines are often associated with loss of genetic diversity, inbreeding and accumulation of deleterious mutations, which can lead to a reduction in fitness and subsequently contribute to extinction. Using temporal approaches can help us understand the effects of population declines on genetic diversity in real time. Sequencing pre-decline as well as post-decline mitogenomes representing all the remaining mitochondrial diversity, we estimated the loss of genetic diversity in the critically endangered kākāpō ( Strigops habroptilus ). We detected a signal of population expansion coinciding with the end of the Pleistocene last glacial maximum (LGM). Also, we found some evidence for northern and southern lineages, supporting the hypothesis that the species may have been restricted to isolated northern and southern refugia during the LGM. We observed an important loss of neutral genetic diversity associated with European settlement in New Zealand but we could not exclude a population decline associated with Polynesian settlement in New Zealand. However, we did not find evidence for fixation of deleterious mutations. We argue that despite high pre-decline genetic diversity, a rapid and range-wide decline combined with the lek mating system, and life-history traits of kākāpō contributed to a rapid loss of genetic diversity following severe population declines.

Elephantids are the world's most iconic megafaunal family, yet there is no comprehensive genomic assessment of their relationships. We report a total of 14 genomes, including 2 from the American mastodon, which is an extinct elephantid relative, and 12 spanning all three extant and three extinct elephantid species including an ∼120,000-y-old straight-tusked elephant, a Columbian mammoth, and woolly mammoths. Earlier genetic studies modeled elephantid evolution via simple bifurcating trees, but here we show that interspecies hybridization has been a recurrent feature of elephantid evolution. We found that the genetic makeup of the straight-tusked elephant, previously placed as a sister group to African forest elephants based on lower coverage data, in fact comprises three major components. Most of the straight-tusked elephant's ancestry derives from a lineage related to the ancestor of African elephants while its remaining ancestry consists of a large contribution from a lineage related to forest elephants and another related to mammoths. Columbian and woolly mammoths also showed evidence of interbreeding, likely following a latitudinal cline across North America. While hybridization events have shaped elephantid history in profound ways, isolation also appears to have played an important role. Our data reveal nearly complete isolation between the ancestors of the African forest and savanna elephants for ∼500,000 y, providing compelling justification for the conservation of forest and savanna elephants as separate species., Competing Interests: The authors declare no conflict of interest.

Many species have undergone dramatic population size declines over the past centuries. Although stochastic genetic processes during and after such declines are thought to elevate the risk of extinction, comparative analyses of genomic data from several endangered species suggest little concordance between genome-wide diversity and current population sizes. This is likely because species-specific life-history traits and ancient bottlenecks overshadow the genetic effect of recent demographic declines. Therefore, we advocate that temporal sampling of genomic data provides a more accurate approach to quantify genetic threats in endangered species. Specifically, genomic data from predecline museum specimens will provide valuable baseline data that enable accurate estimation of recent decreases in genome-wide diversity, increases in inbreeding levels, and accumulation of deleterious genetic variation., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

The onset of the Holocene was associated with a global temperature increase, which led to a rise in sea levels and isolation of the last surviving population of woolly mammoths on Wrangel Island. Understanding what happened with the population's genetic diversity at the time of the isolation and during the ensuing 6000 years can help clarify the effects of bottlenecks and subsequent limited population sizes in species approaching extinction. Previous genetic studies have highlighted questions about how the Holocene Wrangel population was established and how the isolation event affected genetic diversity. Here, we generated high-quality mitogenomes from 21 radiocarbon-dated woolly mammoths to compare the ancestral large and genetically diverse Late Pleistocene Siberian population and the small Holocene Wrangel population. Our results indicate that mitogenome diversity was reduced to one single haplotype at the time of the isolation, and thus that the Holocene Wrangel Island population was established by a single maternal lineage. Moreover, we show that the ensuing small effective population size coincided with fixation of a nonsynonymous mutation, and a comparative analysis of mutation rates suggests that the evolutionary rate was accelerated in the Holocene population. These results suggest that isolation on Wrangel Island led to an increase in the frequency of deleterious genetic variation, and thus are consistent with the hypothesis that strong genetic drift in small populations leads to purifying selection being less effective in removing deleterious mutations.

While present-day taxa are valuable proxies for understanding the biology of extinct species, it is also crucial to examine physical remains in order to obtain a more comprehensive view of their behavior, social structure, and life histories [1, 2]. For example, information on demographic parameters such as age distribution and sex ratios in fossil assemblages can be used to accurately infer socioecological patterns (e.g., [3]). Here we use genomic data to determine the sex of 98 woolly mammoth (Mammuthus primigenius) specimens in order to infer social and behavioral patterns in the last 60,000 years of the species' existence. We report a significant excess of males among the identified samples (69% versus 31%; p < 0.0002). We argue that this male bias among mammoth remains is best explained by males more often being caught in natural traps that favor preservation. We hypothesize that this is a consequence of social structure in proboscideans, which is characterized by matriarchal hierarchy and sex segregation. Without the experience associated with living in a matriarchal family group, or a bachelor group with an experienced bull, young or solitary males may have been more prone to die in natural traps where good preservation is more likely., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Population-scale molecular studies of endangered and cryptic species are often limited by access to high-quality samples. The use of noninvasively collected samples or museum-preserved specimens reduces the pressure on modern populations by removing the need to capture and handle live animals. However, endogenous DNA content in such samples is low, making shotgun sequencing a financially prohibitive approach. Here, we apply a target enrichment method to retrieve mitochondrial genomes from 65 museum specimens and 56 noninvasively collected faecal samples of two endangered great ape species, Grauer's gorilla and the eastern chimpanzee. We show that the applied method is suitable for a wide range of sample types that differ in endogenous DNA content, increasing the proportion of target reads to over 300-fold. By systematically evaluating biases introduced during target enrichment of pooled museum samples, we show that capture is less efficient for fragments shorter or longer than the baits, that the proportion of human contaminating reads increases postcapture although capture efficiency is lower for human compared to gorilla fragments with a gorilla-generated bait, and that the rate of jumping PCR is considerable, but can be controlled for with a double-barcoding approach. We succeed in capturing complete mitochondrial genomes from faecal samples, but observe reduced capture efficiency as sequence divergence increases between the bait and target species. As previously shown for museum specimens, we demonstrate here that mitochondrial genome capture from field-collected faecal samples is a robust and reliable approach for population-wide studies of nonmodel organisms., (© 2017 John Wiley & Sons Ltd.)

Background: An increasing number of studies are addressing the evolutionary genomics of dog domestication, principally through resequencing dog, wolf and related canid genomes. There is, however, only one de novo assembled canid genome currently available against which to map such data - that of a boxer dog (Canis lupus familiaris). We generated the first de novo wolf genome (Canis lupus lupus) as an additional choice of reference, and explored what implications may arise when previously published dog and wolf resequencing data are remapped to this reference., Results: Reassuringly, we find that regardless of the reference genome choice, most evolutionary genomic analyses yield qualitatively similar results, including those exploring the structure between the wolves and dogs using admixture and principal component analysis. However, we do observe differences in the genomic coverage of re-mapped samples, the number of variants discovered, and heterozygosity estimates of the samples., Conclusion: In conclusion, the choice of reference is dictated by the aims of the study being undertaken; if the study focuses on the differences between the different dog breeds or the fine structure among dogs, then using the boxer reference genome is appropriate, but if the aim of the study is to look at the variation within wolves and their relationships to dogs, then there are clear benefits to using the de novo assembled wolf reference genome.

Bird remains that are difficult to identify taxonomically using morphological methods, are common in the palaeontological record. Other types of challenging avian material include artefacts and food items from endangered taxa, as well as remains from aircraft strikes. We here present a DNA-based method that enables taxonomic identification of bird remains, even from material where the DNA is heavily degraded. The method is based on the amplification and sequencing of two short variable parts of the 16S region in the mitochondrial genome. To demonstrate the applicability of this approach, we evaluated the method on a set of Holocene and Late Pleistocene postcranial bird bones from several palaeontological and archaeological sites in Europe with good success., Competing Interests: The authors declare no conflict of interest. The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

The beluga whale is one of three endemic Arctic whales. The species is philopatric, and its migration patterns are passed from mother to calf. Management of the species is informed by the levels of genetic structuring among summer aggregation sites based on mitochondrial D-Loop data. To assess the levels of differentiation across the entire mitochondrial genome within belugas, we present a comparison between the first two complete mitochondrial genomes from opposite sides of their distribution range: Baffin Bay and the Russian Far East. Our analyses reveal that additional phylogenetic insights can be gained from expanding the genetic region analyzed. Further, we estimate the divergence time between the two mitochondrial genomes to be 0.469 MYA., Competing Interests: The authors report no conflicts of interest and are alone responsible for the content and writing of the paper., (© 2017 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.)

Global warming is predicted to cause substantial habitat rearrangements, with the most severe effects expected to occur in high-latitude biomes. However, one major uncertainty is whether species will be able to shift their ranges to keep pace with climate-driven environmental changes. Many recent studies on mammals have shown that past range contractions have been associated with local extinctions rather than survival by habitat tracking. Here, we have used an interdisciplinary approach that combines ancient DNA techniques, coalescent simulations and species distribution modelling, to investigate how two common cold-adapted bird species, willow and rock ptarmigan (Lagopus lagopus and Lagopus muta), respond to long-term climate warming. Contrary to previous findings in mammals, we demonstrate a genetic continuity in Europe over the last 20 millennia. Results from back-casted species distribution models suggest that this continuity may have been facilitated by uninterrupted habitat availability and potentially also the greater dispersal ability of birds. However, our predictions show that in the near future, some isolated regions will have little suitable habitat left, implying a future decrease in local populations at a scale unprecedented since the last glacial maximum., (© 2016 John Wiley & Sons Ltd.)

Near the end of the Pleistocene epoch, populations of the woolly mammoth (Mammuthus primigenius) were distributed across parts of three continents, from western Europe and northern Asia through Beringia to the Atlantic seaboard of North America. Nonetheless, questions about the connectivity and temporal continuity of mammoth populations and species remain unanswered. We use a combination of targeted enrichment and high-throughput sequencing to assemble and interpret a data set of 143 mammoth mitochondrial genomes, sampled from fossils recovered from across their Holarctic range. Our dataset includes 54 previously unpublished mitochondrial genomes and significantly increases the coverage of the Eurasian range of the species. The resulting global phylogeny confirms that the Late Pleistocene mammoth population comprised three distinct mitochondrial lineages that began to diverge ~1.0-2.0 million years ago (Ma). We also find that mammoth mitochondrial lineages were strongly geographically partitioned throughout the Pleistocene. In combination, our genetic results and the pattern of morphological variation in time and space suggest that male-mediated gene flow, rather than large-scale dispersals, was important in the Pleistocene evolutionary history of mammoths.

The Oriental Hornet ( Vespa orientalis ) is a social insect belonging to the Vespiade family (Wasps, Hornets, Yellowjackets), genus Vespa (true Hornets). The oriental hornet is a scavenger and an agricultural pest, especially to bee farmers, but is also recently described as a harvester of solar energy. Here, we report the mitochondrial genome sequence of the Oriental Hornet, Vespa orientalis F ., which may play a vital role in understanding this wasp biology, light trapping and generation of electricity. The mitochondrial genome of this hornet is 16,099 bp in length, containing 13 protein-coding genes, 21 transfer RNA genes, and 2 ribosomal RNA genes. The overall base composition of the heavy-strand is 40.3% A, 5.9% C, 13.2% G, and 40.6% T, the percentages of A and T being higher than that of G and C. The mitochondrial genome of the Oriental Hornet, Vespa orientalis F . represents the first mitogenome of a solar energy harvesting insect., Competing Interests: The author reports no conflicts of interest. The authors alone are responsible for the content and writing of the paper., (© 2017 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.)

The mammalian evolutionary tree has lost several major clades through recent human-caused extinctions. This process of historical biodiversity loss has particularly affected tropical island regions such as the Caribbean, an area of great evolutionary diversification but poor molecular preservation. The most enigmatic of the recently extinct endemic Caribbean mammals are the Nesophontidae, a family of morphologically plesiomorphic lipotyphlan insectivores with no consensus on their evolutionary affinities, and which constitute the only major recent mammal clade to lack any molecular information on their phylogenetic placement. Here, we use a palaeogenomic approach to place Nesophontidae within the phylogeny of recent Lipotyphla. We recovered the near-complete mitochondrial genome and sequences for 17 nuclear genes from a ∼750-year-old Hispaniolan Nesophontes specimen, and identify a divergence from their closest living relatives, the Solenodontidae, more than 40 million years ago. Nesophontidae is thus an older distinct lineage than many extant mammalian orders, highlighting not only the role of island systems as "museums" of diversity that preserve ancient lineages, but also the major human-caused loss of evolutionary history., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Inflow of matter and organisms may strongly affect the local density and diversity of organisms. This effect is particularly evident on shores where organisms with aquatic larval stages enter the terrestrial food web. The identities of such trophic links are not easily estimated as spiders, a dominant group of shoreline predator, have external digestion. We compared trophic links and the prey diversity of spiders on different shore types along the Baltic Sea: on open shores and on shores with a reed belt bordering the water. A priori, we hypothesized that the physical structure of the shoreline reduces the flow between ecosystem and the subsidies across the sea-land interface. To circumvent the lack of morphologically detectable remains of spider prey, we used a combination of stable isotope and molecular gut content analyses. The two tools used for diet analysis revealed complementary information on spider diets. The stable isotope analysis indicated that spiders on open shores had a marine signal of carbon isotopes, while spiders on reedy shores had a terrestrial signal. The molecular analysis revealed a diverse array of dipteran and lepidopteran prey, where spiders on open and reedy shores shared a similar diet with a comparable proportion of chironomids, the larvae of which live in the marine system. Comparing the methods suggests that differences in isotope composition of the two spider groups occurred because of differences in the chironomid diets: as larvae, chironomids of reedy shores likely fed on terrestrial detritus and acquired a terrestrial isotope signature, while chironomids of open shores utilized an algal diet and acquired a marine isotope signature. Our results illustrate how different methods of diet reconstruction may shed light on complementary aspects of nutrient transfer. Overall, they reveal that reed belts can reduce connectivity between habitats, but also function as a source of food for predators.

Harmful effects arising from matings between relatives (inbreeding) is a long-standing observation that is well founded in theory. Empirical evidence for inbreeding depression in natural populations is however rare because of the challenges of assembling pedigrees supplemented with fitness traits. We examined the occurrence of inbreeding and subsequent inbreeding depression using a unique data set containing a genetically verified pedigree with individual fitness traits for a critically endangered arctic fox (Vulpes lagopus) population. The study covered nine years and was comprised of 33 litters with a total of 205 individuals. We recorded that the present population was founded by only five individuals. Over the study period, the population exhibited a tenfold increase in average inbreeding coefficient with a final level corresponding to half-sib matings. Inbreeding mainly occurred between cousins, but we also observed two cases of full-sib matings. The pedigree data demonstrated clear evidence of inbreeding depression on traditional fitness traits where inbred individuals displayed reduced survival and reproduction. Fitness traits were however differently affected by the fluctuating resource abundande. Inbred individuals born at low-quality years displayed reduced first-year survival, while inbred individuals born at high-quality years were less likely to reproduce. The documentation of inbreeding depression in fundamental fitness traits suggests that inbreeding depression can limit population recovery. Introducing new genetic material to promote a genetic rescue effect may thus be necessary for population long-term persistence., (© 2016 John Wiley & Sons Ltd.)

According to the nearly-neutral theory of evolution, the relative strengths of selection and drift shift in favour of drift at small population sizes. Numerous studies have analysed the effect of bottlenecks and small population sizes on genetic diversity in the MHC, which plays a central role in pathogen recognition and immune defense and is thus considered a model example for the study of adaptive evolution. However, to understand changes in genetic diversity at loci under selection, it is necessary to compare the genetic diversity of a population before and after the bottleneck. In this study, we analyse three fragments of the MHC DQA gene in woolly mammoth samples radiocarbon dated to before and after a well-documented bottleneck that took place about ten thousand years ago. Our results indicate a decrease in observed heterozygosity and number of alleles, suggesting that genetic drift had an impact on the variation on MHC. Based on coalescent simulations, we found no evidence of balancing selection maintaining MHC diversity during the Holocene. However, strong trans-species polymorphism among mammoths and elephants points to historical effects of balancing selection on the woolly mammoth lineage.

Recent palaeogenetic studies indicate a highly dynamic history in collared lemmings (Dicrostonyx spp.), with several demographical changes linked to climatic fluctuations that took place during the last glaciation. At the western range margin of D. torquatus, these changes were characterized by a series of local extinctions and recolonizations. However, it is unclear whether this pattern represents a local phenomenon, possibly driven by ecological edge effects, or a global phenomenon that took place across large geographical scales. To address this, we explored the palaeogenetic history of the collared lemming using a next-generation sequencing approach for pooled mitochondrial DNA amplicons. Sequences were obtained from over 300 fossil remains sampled across Eurasia and two sites in North America. We identified five mitochondrial lineages of D. torquatus that succeeded each other through time across Europe and western Russia, indicating a history of repeated population extinctions and recolonizations, most likely from eastern Russia, during the last 50 000 years. The observation of repeated extinctions across such a vast geographical range indicates large-scale changes in the steppe-tundra environment in western Eurasia during the last glaciation. All Holocene samples, from across the species' entire range, belonged to only one of the five mitochondrial lineages. Thus, extant D. torquatus populations only harbour a small fraction of the total genetic diversity that existed across different stages of the Late Pleistocene. In North American samples, haplotypes belonging to both D. groenlandicus and D. richardsoni were recovered from a Late Pleistocene site in south-western Canada. This suggests that D. groenlandicus had a more southern and D. richardsoni a more northern glacial distribution than previously thought. This study provides significant insights into the population dynamics of a small mammal at a large geographical scale and reveals a rather complex demographical history, which could have had bottom-up effects in the Late Pleistocene steppe-tundra ecosystem., (© 2016 John Wiley & Sons Ltd.)

The origin of domestic dogs is poorly understood [1-15], with suggested evidence of dog-like features in fossils that predate the Last Glacial Maximum [6, 9, 10, 14, 16] conflicting with genetic estimates of a more recent divergence between dogs and worldwide wolf populations [13, 15, 17-19]. Here, we present a draft genome sequence from a 35,000-year-old wolf from the Taimyr Peninsula in northern Siberia. We find that this individual belonged to a population that diverged from the common ancestor of present-day wolves and dogs very close in time to the appearance of the domestic dog lineage. We use the directly dated ancient wolf genome to recalibrate the molecular timescale of wolves and dogs and find that the mutation rate is substantially slower than assumed by most previous studies, suggesting that the ancestors of dogs were separated from present-day wolves before the Last Glacial Maximum. We also find evidence of introgression from the archaic Taimyr wolf lineage into present-day dog breeds from northeast Siberia and Greenland, contributing between 1.4% and 27.3% of their ancestry. This demonstrates that the ancestry of present-day dogs is derived from multiple regional wolf populations., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

The processes leading up to species extinctions are typically characterized by prolonged declines in population size and geographic distribution, followed by a phase in which populations are very small and may be subject to intrinsic threats, including loss of genetic diversity and inbreeding. However, whether such genetic factors have had an impact on species prior to their extinction is unclear; examining this would require a detailed reconstruction of a species' demographic history as well as changes in genome-wide diversity leading up to its extinction. Here, we present high-quality complete genome sequences from two woolly mammoths (Mammuthus primigenius). The first mammoth was sequenced at 17.1-fold coverage and dates to ∼4,300 years before present, representing one of the last surviving individuals on Wrangel Island. The second mammoth, sequenced at 11.2-fold coverage, was obtained from an ∼44,800-year-old specimen from the Late Pleistocene population in northeastern Siberia. The demographic trajectories inferred from the two genomes are qualitatively similar and reveal a population bottleneck during the Middle or Early Pleistocene, and a more recent severe decline in the ancestors of the Wrangel mammoth at the end of the last glaciation. A comparison of the two genomes shows that the Wrangel mammoth has a 20% reduction in heterozygosity as well as a 28-fold increase in the fraction of the genome that comprises runs of homozygosity. We conclude that the population on Wrangel Island, which was the last surviving woolly mammoth population, was subject to reduced genetic diversity shortly before it became extinct., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

The global loss of biodiversity continues at an alarming rate. Genomic approaches have been suggested as a promising tool for conservation practice as scaling up to genome-wide data can improve traditional conservation genetic inferences and provide qualitatively novel insights. However, the generation of genomic data and subsequent analyses and interpretations remain challenging and largely confined to academic research in ecology and evolution. This generates a gap between basic research and applicable solutions for conservation managers faced with multifaceted problems. Before the real-world conservation potential of genomic research can be realized, we suggest that current infrastructures need to be modified, methods must mature, analytical pipelines need to be developed, and successful case studies must be disseminated to practitioners., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

Inference of population demographic history has vastly improved in recent years due to a number of technological and theoretical advances including the use of ancient DNA. Approximate Bayesian computation (ABC) stands among the most promising methods due to its simple theoretical fundament and exceptional flexibility. However, limited availability of user-friendly programs that perform ABC analysis renders it difficult to implement, and hence programming skills are frequently required. In addition, there is limited availability of programs able to deal with heterochronous data. Here we present the software BaySICS: Bayesian Statistical Inference of Coalescent Simulations. BaySICS provides an integrated and user-friendly platform that performs ABC analyses by means of coalescent simulations from DNA sequence data. It estimates historical demographic population parameters and performs hypothesis testing by means of Bayes factors obtained from model comparisons. Although providing specific features that improve inference from datasets with heterochronous data, BaySICS also has several capabilities making it a suitable tool for analysing contemporary genetic datasets. Those capabilities include joint analysis of independent tables, a graphical interface and the implementation of Markov-chain Monte Carlo without likelihoods.

Polar bears are uniquely adapted to life in the High Arctic and have undergone drastic physiological changes in response to Arctic climates and a hyper-lipid diet of primarily marine mammal prey. We analyzed 89 complete genomes of polar bear and brown bear using population genomic modeling and show that the species diverged only 479-343 thousand years BP. We find that genes on the polar bear lineage have been under stronger positive selection than in brown bears; nine of the top 16 genes under strong positive selection are associated with cardiomyopathy and vascular disease, implying important reorganization of the cardiovascular system. One of the genes showing the strongest evidence of selection, APOB, encodes the primary lipoprotein component of low-density lipoprotein (LDL); functional mutations in APOB may explain how polar bears are able to cope with life-long elevated LDL levels that are associated with high risk of heart disease in humans., (Copyright © 2014 Elsevier Inc. All rights reserved.)

The Pleistocene glacial cycles resulted in significant changes in species distributions, and it has been discussed whether this caused increased rates of population divergence and speciation. One species that is likely to have evolved during the Pleistocene is the Norwegian lemming (Lemmus lemmus). However, the origin of this species, both in terms of when and from what ancestral taxon it evolved, has been difficult to ascertain. Here, we use ancient DNA recovered from lemming remains from a series of Late Pleistocene and Holocene sites to explore the species' evolutionary history. The results revealed considerable genetic differentiation between glacial and contemporary samples. Moreover, the analyses provided strong support for a divergence time prior to the Last Glacial Maximum (LGM), therefore likely ruling out a postglacial colonization of Scandinavia. Consequently, it appears that the Norwegian lemming evolved from a small population that survived the LGM in an ice-free Scandinavian refugium., (© 2014 John Wiley & Sons Ltd.)

Ancient DNA analyses have provided enhanced resolution of population histories in many Pleistocene taxa. However, most studies are spatially restricted, making inference of species-level biogeographic histories difficult. Here, we analyse mitochondrial DNA (mtDNA) variation in the woolly mammoth from across its Holarctic range to reconstruct its history over the last 200 thousand years (kyr). We identify a previously undocumented major mtDNA lineage in Europe, which was replaced by another major mtDNA lineage 32-34 kyr before present (BP). Coalescent simulations provide support for demographic expansions at approximately 121 kyr BP, suggesting that the previous interglacial was an important driver for demography and intraspecific genetic divergence. Furthermore, our results suggest an expansion into Eurasia from America around 66 kyr BP, coinciding with the first exposure of the Bering Land Bridge during the Late Pleistocene. Bayesian inference indicates Late Pleistocene demographic stability until 20-15 kyr BP, when a severe population size decline occurred.

The Late Pleistocene global extinction of many terrestrial mammal species has been a subject of intensive scientific study for over a century, yet the relative contributions of environmental changes and the global expansion of humans remain unresolved. A defining component of these extinctions is a bias toward large species, with the majority of small-mammal taxa apparently surviving into the present. Here, we investigate the population-level history of a key tundra-specialist small mammal, the collared lemming (Dicrostonyx torquatus), to explore whether events during the Late Pleistocene had a discernible effect beyond the large mammal fauna. Using ancient DNA techniques to sample across three sites in North-West Europe, we observe a dramatic reduction in genetic diversity in this species over the last 50,000 y. We further identify a series of extinction-recolonization events, indicating a previously unrecognized instability in Late Pleistocene small-mammal populations, which we link with climatic fluctuations. Our results reveal climate-associated, repeated regional extinctions in a keystone prey species across the Late Pleistocene, a pattern likely to have had an impact on the wider steppe-tundra community, and one that is concordant with environmental change as a major force in structuring Late Pleistocene biodiversity.

Previous studies have suggested that the presence of sea ice is an important factor in facilitating migration and determining the degree of genetic isolation among contemporary arctic fox populations. Because the extent of sea ice is dependent upon global temperatures, periods of significant cooling would have had a major impact on fox population connectivity and genetic variation. We tested this hypothesis by extracting and sequencing mitochondrial control region sequences from 17 arctic foxes excavated from two late-ninth-century to twelfth-century AD archaeological sites in northeast Iceland, both of which predate the Little Ice Age (approx. sixteenth to nineteenth century). Despite the fact that five haplotypes have been observed in modern Icelandic foxes, a single haplotype was shared among all of the ancient individuals. Results from simulations within an approximate Bayesian computation framework suggest that the rapid increase in Icelandic arctic fox haplotype diversity can only be explained by sea-ice-mediated fox immigration facilitated by the Little Ice Age.

The last glaciation was a dynamic period with strong impact on the demography of many species and populations. In recent years, mitochondrial DNA sequences retrieved from radiocarbon-dated remains have provided novel insights into the history of Late Pleistocene populations. However, genotyping of loci from the nuclear genome may provide enhanced resolution of population-level changes. Here, we use four autosomal microsatellite DNA markers to investigate the demographic history of woolly mammoths (Mammuthus primigenius) in north-eastern Siberia from before 60 000 years ago up until the species' final disappearance c.4000 years ago. We identified two genetic groups, implying a marked temporal genetic differentiation between samples with radiocarbon ages older than 12 thousand radiocarbon years before present (ka) and those younger than 9ka. Simulation-based analysis indicates that this dramatic change in genetic composition, which included a decrease in individual heterozygosity of approximately 30%, was due to a multifold reduction in effective population size. A corresponding reduction in genetic variation was also detected in the mitochondrial DNA, where about 65% of the diversity was lost. We observed no further loss in genetic variation during the Holocene, which suggests a rapid final extinction event., (© 2012 Blackwell Publishing Ltd.)

Low genetic diversity in the endangered Iberian lynx, including lack of mitochondrial control region variation, is thought to result from historical or Pleistocene/Holocene population bottlenecks, and to indicate poor long-term viability. We find no variability in control region sequences from 19 Iberian lynx remains from across the Iberian Peninsula and spanning the last 50,000 years. This is best explained by continuously small female effective population size through time. We conclude that low genetic variability in the Iberian lynx is not in itself a threat to long-term viability, and so should not preclude conservation efforts., (© 2011 Blackwell Publishing Ltd.)

During the Late Pleistocene, the woolly mammoth (Mammuthus primigenius) experienced a series of local extinctions generally attributed to human predation or environmental change. Some small and isolated populations did however survive far into the Holocene. Here, we investigated the genetic consequences of the isolation of the last remaining mammoth population on Wrangel Island. We analysed 741 bp of the mitochondrial DNA and found a loss of genetic variation in relation to the isolation event, probably caused by a demographic bottleneck or a founder event. However, in spite of ca 5000 years of isolation, we did not detect any further loss of genetic variation. Together with the relatively high number of mitochondrial haplotypes on Wrangel Island near the final disappearance, this suggests a sudden extinction of a rather stable population.

Climate change in the past has led to significant changes in species' distributions. However, how individual species respond to climate change depends largely on their adaptations and environmental tolerances. In the Quaternary, temperate-adapted taxa are in general confined to refugia during glacials while cold-adapted taxa are in refugia during interglacials. In the Northern Hemisphere, evidence appears to be mounting that in addition to traditional southern refugia for temperate species, cryptic refugia existed in the North during glacials. Equivalent cryptic southern refugia, to the south of the more conventional high-latitude polar refugia, exist in montane areas during periods of warm climate, such as the current interglacial. There is also a continental/oceanic longitudinal gradient, which should be included in a more complete consideration of the interaction between species ranges and climates. Overall, it seems clear that there is large variation in both the size of refugia and the duration during which species are confined to them. This has implications for the role of refugia in the evolution of species and their genetic diversity.

Background: The Neotropical ovenbird-woodcreeper family (Furnariidae) is an avian group characterized by exceptionally diverse ecomorphological adaptations. For instance, members of the family are known to construct nests of a remarkable variety. This offers a unique opportunity to examine whether changes in nest design, accompanied by expansions into new habitats, facilitates diversification. We present a multi-gene phylogeny and age estimates for the ovenbird-woodcreeper family and use these results to estimate the degree of convergent evolution in both phenotype and habitat utilisation. Furthermore, we discuss whether variation in species richness among ovenbird clades could be explained by differences in clade-specific diversification rates, and whether these rates differ among lineages with different nesting habits. In addition, the systematic positions of some enigmatic ovenbird taxa and the postulated monophyly of some species-rich genera are evaluated., Results: The phylogenetic results reveal new examples of convergent evolution and show that ovenbirds have independently colonized open habitats at least six times. The calculated age estimates suggest that the ovenbird-woodcreeper family started to diverge at ca 33 Mya, and that the timing of habitat shifts into open environments may be correlated with the aridification of South America during the last 15 My. The results also show that observed large differences in species richness among clades can be explained by a substantial variation in net diversification rates. The synallaxines, which generally are adapted to dry habitats and build exposed vegetative nests, had the highest diversification rate of all major furnariid clades., Conclusion: Several key features may have played an important role for the radiation and evolution of convergent phenotypes in the ovenbird-woodcreeper family. Our results suggest that changes in nest building strategy and adaptation to novel habitats may have played an important role in a diversification that included multiple radiations into more open and bushy environments. The synallaxines were found to have had a particularly high diversification rate, which may be explained by their ability to build exposed vegetative nests and thus to expand into a variety of novel habitats that emerged during a period of cooling and aridification in South America.

We report the first two complete mitochondrial genome sequences of the thylacine (Thylacinus cynocephalus), or so-called Tasmanian tiger, extinct since 1936. The thylacine's phylogenetic position within australidelphian marsupials has long been debated, and here we provide strong support for the thylacine's basal position in Dasyuromorphia, aided by mitochondrial genome sequence that we generated from the extant numbat (Myrmecobius fasciatus). Surprisingly, both of our thylacine sequences differ by 11%-15% from putative thylacine mitochondrial genes in GenBank, with one of our samples originating from a direct offspring of the previously sequenced individual. Our data sample each mitochondrial nucleotide an average of 50 times, thereby providing the first high-fidelity reference sequence for thylacine population genetics. Our two sequences differ in only five nucleotides out of 15,452, hinting at a very low genetic diversity shortly before extinction. Despite the samples' heavy contamination with bacterial and human DNA and their temperate storage history, we estimate that as much as one-third of the total DNA in each sample is from the thylacine. The microbial content of the two thylacine samples was subjected to metagenomic analysis, and showed striking differences between a wild-captured individual and a born-in-captivity one. This study therefore adds to the growing evidence that extensive sequencing of museum collections is both feasible and desirable, and can yield complete genomes.

We report five new complete mitochondrial DNA (mtDNA) genomes of Siberian woolly mammoth (Mammuthus primigenius), sequenced with up to 73-fold coverage from DNA extracted from hair shaft material. Three of the sequences present the first complete mtDNA genomes of mammoth clade II. Analysis of these and 13 recently published mtDNA genomes demonstrates the existence of two apparently sympatric mtDNA clades that exhibit high interclade divergence. The analytical power afforded by the analysis of the complete mtDNA genomes reveals a surprisingly ancient coalescence age of the two clades, approximately 1-2 million years, depending on the calibration technique. Furthermore, statistical analysis of the temporal distribution of the (14)C ages of these and previously identified members of the two mammoth clades suggests that clade II went extinct before clade I. Modeling of protein structures failed to indicate any important functional difference between genomes belonging to the two clades, suggesting that the loss of clade II more likely is due to genetic drift than a selective sweep.

Models for the development of species distribution in Europe typically invoke restriction in three temperate Mediterranean refugia during glaciations, from where recolonization of central and northern Europe occurred. The brown bear, Ursus arctos, is one of the taxa from which this model is derived. Sequence data generated from brown bear fossils show a complex phylogeographical history for western European populations. Long-term isolation in separate refugia is not required to explain our data when considering the palaeontological distribution of brown bears. We propose continuous gene flow across southern Europe, from which brown bear populations expanded after the last glaciation.

Unlike Oceanic islands, the islands of the Arctic Sea are not completely isolated from migration by terrestrial vertebrates. The pack ice connects many Arctic Sea islands to the mainland during winter months. The Arctic fox (Alopex lagopus), which has a circumpolar distribution, populates numerous islands in the Arctic Sea. In this study, we used genetic data from 20 different populations, spanning the entire distribution of the Arctic fox, to identify barriers to dispersal. Specifically, we considered geographical distance, occurrence of sea ice, winter temperature, ecotype, and the presence of red fox and polar bear as nonexclusive factors that influence the dispersal behaviour of individuals. Using distance-based redundancy analysis and the BIOENV procedure, we showed that occurrence of sea ice is the key predictor and explained 40-60% of the genetic distance among populations. In addition, our analysis identified the Commander and Pribilof Islands Arctic populations as genetically unique suggesting they deserve special attention from a conservation perspective.

Although the application of sequencing-by-synthesis techniques to DNA extracted from bones has revolutionized the study of ancient DNA, it has been plagued by large fractions of contaminating environmental DNA. The genetic analyses of hair shafts could be a solution: We present 10 previously unexamined Siberian mammoth (Mammuthus primigenius) mitochondrial genomes, sequenced with up to 48-fold coverage. The observed levels of damage-derived sequencing errors were lower than those observed in previously published frozen bone samples, even though one of the specimens was >50,000 14C years old and another had been stored for 200 years at room temperature. The method therefore sets the stage for molecular-genetic analysis of museum collections.

How species respond to an increased availability of habitat, for example at the end of the last glaciation, has been well established. In contrast, little is known about the opposite process, when the amount of habitat decreases. The hypothesis of habitat tracking predicts that species should be able to track both increases and decreases in habitat availability. The alternative hypothesis is that populations outside refugia become extinct during periods of unsuitable climate. To test these hypotheses, we used ancient DNA techniques to examine genetic variation in the arctic fox (Alopex lagopus) through an expansion/contraction cycle. The results show that the arctic fox in midlatitude Europe became extinct at the end of the Pleistocene and did not track the habitat when it shifted to the north. Instead, a high genetic similarity between the extant populations in Scandinavia and Siberia suggests an eastern origin for the Scandinavian population at the end of the last glaciation. These results provide new insights into how species respond to climate change, since they suggest that populations are unable to track decreases in habitat avaliability. This implies that arctic species may be particularly vulnerable to increases in global temperatures.