Search

Predation is an important ecological process that can significantly impact the maintenance of ecosystem services. In arctic environments, the relative ecological importance of predation is thought to be increasing due to climate change, partly because of increased productivity with rising temperatures. Therefore, understanding predator–prey interactions in arctic ecosystems is vital for the sustainable management of these northern regions. Network theory provides a framework for quantifying the structures of ecological interactions. In this study, we use dietary observations on mammalian and avian predators in a high arctic region, including isolated peninsulas on Ellesmere Island and north Greenland, to construct bipartite trophic networks. We quantify the complexity, specialization, and nested as well as modular structures of these networks and also determine if these properties varied among the peninsulas. Mammal prey remains were the dominant diet item for all predators, but there was spatial variation in diet composition among peninsulas. The predator–prey networks were less complex, had more specialized interactions, and were more nested and more modular than random expectations. However, the networks displayed only moderate levels of modularity. Predator species had less specialized interactions with prey than prey had with predators. All network properties differed among the peninsulas, which highlights that ecosystems often show complex responses to environmental characteristics. We suggest that gaining knowledge about spatial variation in the characteristics of predator–prey interactions can enhance our ability to manage ecosystems exposed to environmental perturbations, particularly in high arctic environments subject to rapid environmental change.

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

The iconic Cape buffalo has experienced several documented population declines in recent history. These declines have been largely attributed to the late 19th century rinderpest pandemic. However, the effect of the rinderpest pandemic on their genetic diversity remains contentious, and other factors that have potentially affected this diversity include environmental changes during the Pleistocene, range expansions and recent human activity. Motivated by this, we present analyses of whole genome sequencing data from 59 individuals from across the Cape buffalo range to assess present-day levels of genome-wide genetic diversity and what factors have influenced these levels. We found that the Cape buffalo has high average heterozygosity overall (0.40%), with the two southernmost populations having significantly lower heterozygosity levels (0.33% and 0.29%) on par with that of the domesticated water buffalo (0.29%). Interestingly, we found that these lower levels are probably due to recent inbreeding (average fraction of runs of homozygosity 23.7% and 19.9%) rather than factors further back in time during the Pleistocene. Moreover, detailed investigations of recent demographic history show that events across the past three centuries were the main drivers of the exceptional loss of genetic diversity in the southernmost populations, coincident with the onset of colonialism in the southern extreme of the Cape buffalo range. Hence, our results add to the growing body of studies suggesting that multiple recent human-mediated impacts during the colonial period caused massive losses of large mammal abundance in southern Africa.

Non-invasive biological samples benefit studies that investigate rare, elusive, endangered, or dangerous species. Integrating genomic techniques that use non-invasive biological sampling with advances in computational approaches can benefit and inform wildlife conservation and management. Here, we used non-invasive fecal DNA samples to generate low- to medium-coverage genomes (e.g., >90% of the complete nuclear genome at six X-fold coverage) and metagenomic sequences, combining widely available and accessible DNA collection cards with commonly used DNA extraction and library building approaches. DNA preservation cards are easy to transport and can be stored non-refrigerated, avoiding cumbersome or costly sample methods. The genomic library construction and shotgun sequencing approach did not require enrichment or targeted DNA amplification. The utility and potential of the data generated was demonstrated through genome scale and metagenomic analyses of zoo and free-ranging African savanna elephants (Loxodonta africana). Fecal samples collected from free-ranging individuals contained an average of 12.41% (5.54–21.65%) endogenous elephant DNA. Clustering of these elephants with others from the same geographic region was demonstrated by a principal component analysis of genetic variation using nuclear genome-wide SNPs. Metagenomic analyses identified taxa that included Loxodonta, green plants, fungi, arthropods, bacteria, viruses and archaea, showcasing the utility of this approach for addressing complementary questions based on host-associated DNA, e.g., pathogen and parasite identification. The molecular and bioinformatic analyses presented here contributes towards the expansion and application of genomic techniques to conservation science and practice., Non-invasive biological samples benefit studies that investigate rare, elusive, endangered, or dangerous species. Integrating genomic techniques that use non-invasive biological sampling with advances in computational approaches can benefit and inform wildlife conservation and management. Here, we used non-invasive fecal DNA samples to generate low- to medium-coverage genomes (e.g., >90% of the complete nuclear genome at six X-fold coverage) and metagenomic sequences, combining widely available and accessible DNA collection cards with commonly used DNA extraction and library building approaches. DNA preservation cards are easy to transport and can be stored non-refrigerated, avoiding cumbersome or costly sample methods. The genomic library construction and shotgun sequencing approach did not require enrichment or targeted DNA amplification. The utility and potential of the data generated was demonstrated through genome scale and metagenomic analyses of zoo and free-ranging African savanna elephants (Loxodonta africana). Fecal samples collected from free-ranging individuals contained an average of 12.41% (5.54–21.65%) endogenous elephant DNA. Clustering of these elephants with others from the same geographic region was demonstrated by a principal component analysis of genetic variation using nuclear genome-wide SNPs. Metagenomic analyses identified taxa that included Loxodonta, green plants, fungi, arthropods, bacteria, viruses and archaea, showcasing the utility of this approach for addressing complementary questions based on host-associated DNA, e.g., pathogen and parasite identification. The molecular and bioinformatic analyses presented here contributes towards the expansion and application of genomic techniques to conservation science and practice.

After over a decade of developments in field collection, laboratory methods and advances in high-throughput sequencing, contamination remains a key issue in ancient DNA research. Currently, human and microbial contaminant DNA still impose challenges on cost-effective sequencing and accurate interpretation of ancient DNA data.

Rapid and cost-effective retrieval of endogenous DNA from ancient specimens remains a limiting factor in palaeogenomic research. Many methods have been developed to increase ancient DNA yield, but modifications to existing protocols are often based on personal experience rather than systematic testing. Here, we present a new silica column-based extraction protocol, where optimizations were tested in controlled experiments. Using relatively well-preserved permafrost samples, we tested the efficiency of pretreatment of bone and tooth powder with a bleach wash and a predigestion step. We also tested the recovery efficiency of MinElute and QIAquick columns, as well as Vivaspin columns with two molecular weight cut-off values. Finally, we tested the effect of uracil-treatment with two different USER enzyme concentrations. We find that neither bleach wash combined with a predigestion step, nor predigestion by itself, significantly increased sequencing efficiency. Initial results, however, suggest that MinElute columns are more efficient for ancient DNA extractions than QIAquick columns, whereas different molecular weight cut-off values in centrifugal concentrator columns did not have an effect. Uracil treatments are effective at removing DNA damage even at concentrations of 0.15 U/µL (as compared to 0.3 U/µL) of ancient DNA extracts.

To understand the causes and implications of an extinction event, detailed information is necessary. However, this can be challenging when working with poorly resolved paleontological data sets. One approach to increase the data resolution is by combining different methods. In this study, we used both radiocarbon and genetic data to reconstruct the population history and extinction dynamics of the woolly mammoth in northern Siberia. We generated 88 new radiocarbon dates and combined these with previously published dates from 626 specimens to construct Bayesian age models. These models show that mammoths disappeared on the eastern Siberian mainland before the onset of the Younger Dryas (12.9–11.7 ky cal BP). Mammoths did however persist in the northernmost parts of central and western Siberia until the early Holocene. Further genetic results of 131 high quality mitogenomes, including 22 new mitogenomes generated in this study, support the hypothesis that mammoths from, or closely related to, a central and/or west- Siberian population recolonized Wrangel Island over the now submerged northern Siberian plains. As mammoths became trapped on the island due to rising sea levels, they lived another ca. 6000 years on Wrangel Island before eventually going extinct ca. 4000 years ago.

Species worldwide are subject to contractions in both abundance and geographical range, and their persistence in a changing environment may thus depend on the ability to survive in small and fragmented populations. Despite the urgent need to understand how extinction works, our knowledge of pre-extinction genetic processes is limited because techniques allowing population and conservation genomics to be studied in wild threatened populations have become available only recently. In this thesis, I used the last surviving population of the woolly mammoth (Mammuthus primigenius) as a model for studying pre-extinction population dynamics. I used ancient DNA as a tool to study microevolutionary processes in real time, analysing genetic changes in response to environmental shifts at the end of the last Ice Age and exploring impacts of genetic drift and inbreeding as woolly mammoths became isolated on Wrangel Island and survived for 6000 years at small population size. Using mitochondrial genomes, I found evidence of a founder effect that decreased the maternal diversity to a single lineage at the time when mammoths became trapped on Wrangel Island (~10,500 years ago). Moreover, a two- to three-fold higher mitochondrial mutation rate in Holocene and a fixed, potentially detrimental mutation in the ATP6 gene encoding for one of the key enzymes of the oxidative phosphorylation pathway, is consistent with the hypothesis that selection is less effective in removing deleterious mutations in small populations. A loss of diversity was also observed in an immunity gene that belongs to the major histocompatibility complex (MHC), even though the MHC is considered to be under balancing selection. Low-coverage genomic data was analysed in order to estimate endogenous DNA content and molecular sex of the mammoth samples. The observation of a male bias (69%) in the sex ratio led to the conclusion that male mammoths were more likely to die in a way that ensured good preservation. Another, At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 4: Manuscript. Paper 5: Manuscript.

Arctic environments have relatively simple ecosystems. Yet, we still lack knowledge of the spatio-temporal dynamics of many Arctic organisms and how they are affected by local and regional processes. The Arctic hare (Lepus arcticus) is a large lagomorph endemic to high Arctic environments in Canada and Greenland. Current knowledge about this herbivore is scarce and the temporal and spatial dynamics of their populations are poorly understood. Here, we present observations on Arctic hares in two sites on north Greenland (Hall and Washington lands) and one adjacent site on Ellesmere Island (Judge Daly Promontory). We recorded a large range of group sizes from 1 to 135 individuals, as well as a substantial variation in hare densities among the three sites (Hall land: 0 animals/100 km(2), Washington land 14.5-186.7 animals/100 km(2), Judge Daly Promontory 0.18-2.95 animals/100 km(2)). However, pellet counts suggested that both Hall land and Judge Daly Promontory hosted larger populations at other times. We suggest that our results could have been caused by three spatially differentiated populations with asynchronous population fluctuations. With food limitation being a likely driver behind the observed variation, we argue that food limitation likely interacts with predation and competition in shaping the spatial dynamics of Arctic hares in this region.

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.

A large number of species have recently recovered from near-extinction events. Understanding the genetic consequences of severe population declines followed by demographic recoveries is key to predict the long-term viability of species in order to mitigate future extinction risks. Although these species have avoided the immediate extinction threat, their long-term viability remains questionable due to the genetic consequences of population declines, which are not understood on a time scale beyond a few generations. The woolly mammoth (Mammuthus primigenius) population on Wrangel Island is an excellent model system to investigate long-term genetic consequences of a population bottleneck. Mammoths became isolated on the island in the early Holocene due to rising sea levels, and persisted for over 200 generations (~6,000 years) before becoming extinct ~4,000 years ago. To study the evolutionary processes leading up to the extinction of the woolly mammoth on the island, we analysed 21 Siberian woolly mammoth genomes, including that of one of the last known mammoths. Our results show that the Wrangel Island mammoths recovered quickly from an initially severe bottleneck, and subsequently remained demographically stable during the ensuing 6 millennia. Further, we find that highly deleterious mutations were gradually purged from the population, whereas there was an accumulation of mildly deleterious mutations. The gradual purging of highly deleterious mutations suggests an ongoing inbreeding depression that lasted for hundreds of generations. This time-lag between demographic and genetic recovery has wide-ranging implications for conservation management of recently bottlenecked present-day populations.