A palaeogenomic investigation of overharvest implications in an endemic wild reindeer subspecies

Overharvest can severely reduce the abundance and distribution of a species and thereby impact its genetic diversity and threaten its future viability. Overharvest remains an ongoing issue for Arctic mammals, which due to climate change now also confront one of the fastest changing environments on Earth. The high-arctic Svalbard reindeer (Rangifer tarandus platyrhynchus), endemic to Svalbard, experienced a harvest-induced demographic bottleneck that occurred during the 17-20th centuries. Here, we investigate changes in genetic diversity, population structure, and gene-specific differentiation during and after this overharvesting event. Using whole-genome shotgun sequencing, we generated the first ancient and historical nuclear (n = 11) and mitochondrial (n = 18) genomes from Svalbard reindeer (up to 4000 BP) and integrated these data with a large collection of modern genome sequences (n = 90) to infer temporal changes. We show that hunting resulted in major genetic changes and restructuring in reindeer populations. Near-extirpation followed by pronounced genetic drift has altered the allele frequencies of important genes contributing to diverse biological functions. Median heterozygosity was reduced by 26%, while the mitochondrial genetic diversity was reduced only to a limited extent, likely due to already low pre-harvest diversity and a complex post-harvest recolonization process. Such genomic erosion and genetic isolation of populations due to past anthropogenic disturbance will likely play a major role in metapopulation dynamics (i.e., extirpation, recolonization) under further climate change. Our results from a high-arctic case study therefore emphasize the need to understand the long-term interplay of past, current, and future stressors in wildlife conservation.

Contrasting genomic consequences of anthropogenic reintroduction and natural recolonization in high-arctic wild reindeer

Anthropogenic reintroduction can supplement natural recolonization in reestablishing a species' distribution and abundance. However, both reintroductions and recolonizations can give rise to founder effects that reduce genetic diversity and increase inbreeding, potentially causing the accumulation of genetic load and reduced fitness. Most current populations of the endemic high-arctic Svalbard reindeer (Rangifer tarandus platyrhynchus) originate from recent reintroductions or recolonizations following regional extirpations due to past overharvesting. We investigated and compared the genomic consequences of these two paths to reestablishment using whole-genome shotgun sequencing of 100 Svalbard reindeer across their range. We found little admixture between reintroduced and natural populations. Two reintroduced populations, each founded by 12 individuals around four decades (i.e. 8 reindeer generations) ago, formed two distinct genetic clusters. Compared to the source population, these populations showed only small decreases in genome-wide heterozygosity and increases in inbreeding and lengths of runs of homozygosity. In contrast, the two naturally recolonized populations without admixture possessed much lower heterozygosity, higher inbreeding and longer runs of homozygosity, possibly caused by serial population founder effects and/or fewer or more genetically related founders than in the reintroduction events. Naturally recolonized populations can thus be more vulnerable to the accumulation of genetic load than reintroduced populations. This suggests that in some organisms even small-scale reintroduction programs based on genetically diverse source populations can be more effective than natural recolonization in establishing genetically diverse populations. These findings warrant particular attention in the conservation and management of populations and species threatened by habitat fragmentation and loss.

Circumpolar status of Arctic ptarmigan: Population dynamics and trends

Rock ptarmigan (Lagopus muta) and willow ptarmigan (L. lagopus) are Arctic birds with a circumpolar distribution but there is limited knowledge about their status and trends across their circumpolar distribution. Here, we compiled information from 90 ptarmigan study sites from 7 Arctic countries, where almost half of the sites are still monitored. Rock ptarmigan showed an overall negative trend on Iceland and Greenland, while Svalbard and Newfoundland had positive trends, and no significant trends in Alaska. For willow ptarmigan, there was a negative trend in mid-Sweden and eastern Russia, while northern Fennoscandia, North America and Newfoundland had no significant trends. Both species displayed some periods with population cycles (short 3-6 years and long 9-12 years), but cyclicity changed through time for both species. We propose that simple, cost-efficient systematic surveys that capture the main feature of ptarmigan population dynamics can form the basis for citizen science efforts in order to fill knowledge gaps for the many regions that lack systematic ptarmigan monitoring programs.

Spatial heterogeneity in climate change effects decouples the long-term dynamics of wild reindeer populations in the high Arctic

The 'Moran effect' predicts that dynamics of populations of a species are synchronized over similar distances as their environmental drivers. Strong population synchrony reduces species viability, but spatial heterogeneity in density dependence, the environment, or its ecological responses may decouple dynamics in space, preventing extinctions. How such heterogeneity buffers impacts of global change on large-scale population dynamics is not well studied. Here, we show that spatially autocorrelated fluctuations in annual winter weather synchronize wild reindeer dynamics across high-Arctic Svalbard, while, paradoxically, spatial variation in winter climate trends contribute to diverging local population trajectories. Warmer summers have improved the carrying capacity and apparently led to increased total reindeer abundance. However, fluctuations in population size seem mainly driven by negative effects of stochastic winter rain-on-snow (ROS) events causing icing, with strongest effects at high densities. Count data for 10 reindeer populations 8-324 km apart suggested that density-dependent ROS effects contributed to synchrony in population dynamics, mainly through spatially autocorrelated mortality. By comparing one coastal and one 'continental' reindeer population over four decades, we show that locally contrasting abundance trends can arise from spatial differences in climate change and responses to weather. The coastal population experienced a larger increase in ROS, and a stronger density-dependent ROS effect on population growth rates, than the continental population. In contrast, the latter experienced stronger summer warming and showed the strongest positive response to summer temperatures. Accordingly, contrasting net effects of a recent climate regime shift-with increased ROS and harsher winters, yet higher summer temperatures and improved carrying capacity-led to negative and positive abundance trends in the coastal and continental population respectively. Thus, synchronized population fluctuations by climatic drivers can be buffered by spatial heterogeneity in the same drivers, as well as in the ecological responses, averaging out climate change effects at larger spatial scales.