Colour moult phenology and camouflage mismatch in polymorphic populations of Arctic foxes.

Species that seasonally moult from brown to white to match snowy backgrounds become conspicuous and experience increased predation risk as snow cover duration declines. Long-term adaptation to camouflage mismatch in a changing climate might occur through phenotypic plasticity in colour moult phenology and or evolutionary shifts in moult rate or timing. Also, adaptation may include evolutionary shifts towards winter brown phenotypes that forgo the winter white moult. Most studies of these processes have occurred in winter white populations, with little attention to polymorphic populations with sympatric winter brown and winter white morphs. Here, we used remote camera traps to record moult phenology and mismatch in two polymorphic populations of Arctic foxes in Sweden over 2 years. We found that the colder, more northern population moulted earlier in the autumn and later in the spring. Next, foxes moulted earlier in the autumn and later in the spring during colder and snowier years. Finally, white foxes experienced relatively low camouflage mismatch while blue foxes were mismatched against snowy backgrounds most of the autumn through the spring. Because the brown-on-white mismatch imposes no evident costs, we predict that as snow duration decreases, increasing blue morph frequencies might help facilitate species persistence.

Resources and predation: drivers of sociality in a cyclic mesopredator.

In socially flexible species, the tendency to live in groups is expected to vary through a trade-off between costs and benefits, determined by ecological conditions. The Resource Dispersion Hypothesis predicts that group size changes in response to patterns in resource availability. An additional dimension is described in Hersteinsson's model positing that sociality is further affected by a cost-benefit trade-off related to predation pressure. In the arctic fox (Vulpes lagopus), group-living follows a regional trade-off in resources' availability and intra-guild predation pressure. However, the effect of local fluctuations is poorly known, but offers an unusual opportunity to test predictions that differ between the two hypotheses in systems where prey availability is linked to intra-guild predation. Based on 17-year monitoring of arctic fox and cyclic rodent prey populations, we addressed the Resource Dispersion Hypothesis and discuss the results in relation to the impact of predation in Hersteinsson's model. Group-living increased with prey density, from 7.7% (low density) to 28% (high density). However, it remained high (44%) despite a rodent crash and this could be explained by increased benefits from cooperative defence against prey switching by top predators. We conclude that both resource abundance and predation pressure are factors underpinning the formation of social groups in fluctuating ecosystems.

Fur colour in the Arctic fox: genetic architecture and consequences for fitness.

Genome-wide association studies provide good opportunities for studying the genetic basis of adaptive traits in wild populations. Yet, previous studies often failed to identify major effect genes. In this study, we used high-density single nucleotide polymorphism and individual fitness data from a wild non-model species. Using a whole-genome approach, we identified the MC1R gene as the sole causal gene underlying Arctic fox Vulpes lagopus fur colour. Further, we showed the adaptive importance of fur colour genotypes through measures of fitness that link ecological and evolutionary processes. We found a tendency for blue foxes that are heterozygous at the fur colour locus to have higher fitness than homozygous white foxes. The effect of genotype on fitness was independent of winter duration but varied with prey availability, with the strongest effect in years of increasing rodent populations. MC1R is located in a genomic region with high gene density, and we discuss the potential for indirect selection through linkage and pleiotropy. Our study shows that whole-genome analyses can be successfully applied to wild species and identify major effect genes underlying adaptive traits. Furthermore, we show how this approach can be used to identify knowledge gaps in our understanding of interactions between ecology and evolution.

Genomic and fitness consequences of inbreeding in an endangered carnivore.

Reduced fitness through genetic drift and inbreeding is a major threat to small and isolated populations. Although previous studies have generally used genetically verified pedigrees to document effects of inbreeding and gene flow, these often fail to capture the whole inbreeding history of the species. By assembling a draft arctic fox (Vulpes lagopus) genome and resequencing complete genomes of 23 additional foxes born before and after a well-documented immigration event in Scandinavia, we here look into the genomic consequences of inbreeding and genetic rescue. We found a difference in genome-wide diversity, with 18% higher heterozygosity and 81% lower FROH in immigrant F1 compared to native individuals. However, more distant descendants of immigrants (F2, F3) did not show the same pattern. We also found that foxes with lower inbreeding had higher probability to survive their first year of life. Our results demonstrate the important link between genetic variation and fitness as well as the transient nature of genetic rescue. Moreover, our results have implications in conservation biology as they demonstrate that inbreeding depression can effectively be detected in the wild by a genomic approach.

Low Persistence of Genetic Rescue Across Generations in the Arctic Fox (Vulpes lagopus).

Genetic rescue can facilitate the recovery of small and isolated populations suffering from inbreeding depression. Long-term effects are however complex, and examples spanning over multiple generations under natural conditions are scarce. The aim of this study was to test for long-term effects of natural genetic rescue in a small population of Scandinavian Arctic foxes (Vulpes lagopus). By combining a genetically verified pedigree covering almost 20 years with a long-term dataset on individual fitness (n = 837 individuals), we found no evidence for elevated fitness in immigrant F2 and F3 compared to native inbred foxes. Population inbreeding levels showed a fluctuating increasing trend and emergence of inbreeding within immigrant lineages shortly after immigration. Between 0-5 and 6-9 years post immigration, the average number of breeding adults decreased by almost 22% and the average proportion of immigrant ancestry rose from 14% to 27%. Y chromosome analysis revealed that 2 out of 3 native male lineages were lost from the gene pool, but all founders represented at the time of immigration were still contributing to the population at the end of the study period through female descendants. The results highlight the complexity of genetic rescue and suggest that beneficial effects can be brief. Continuous gene flow may be needed for small and threatened populations to recover and persist in a longer time perspective.

Genetic rescue in an inbred Arctic fox (Vulpes lagopus) population.

Isolation of small populations can reduce fitness through inbreeding depression and impede population growth. Outcrossing with only a few unrelated individuals can increase demographic and genetic viability substantially, but few studies have documented such genetic rescue in natural mammal populations. We investigate the effects of immigration in a subpopulation of the endangered Scandinavian arctic fox (Vulpes lagopus), founded by six individuals and isolated for 9 years at an extremely small population size. Based on a long-term pedigree (105 litters, 543 individuals) combined with individual fitness traits, we found evidence for genetic rescue. Natural immigration and gene flow of three outbred males in 2010 resulted in a reduction in population average inbreeding coefficient (f), from 0.14 to 0.08 within 5 years. Genetic rescue was further supported by 1.9 times higher juvenile survival and 1.3 times higher breeding success in immigrant first-generation offspring compared with inbred offspring. Five years after immigration, the population had more than doubled in size and allelic richness increased by 41%. This is one of few studies that has documented genetic rescue in a natural mammal population suffering from inbreeding depression and contributes to a growing body of data demonstrating the vital connection between genetics and individual fitness.

The influence of spatiotemporal conditions and personality on survival in reintroductions-evolutionary implications.

Personality exists in non-human animals and can impact fitness. There is, however, a shortage of empirical studies in certain areas within the field, and fundamental evolutionary theory on personality remains largely untested. For example, little is known on how variation in personality is maintained over evolutionary time. Theory suggests that fluctuating selection pressures due to spatiotemporal variation in conditions, e.g. food availability, is a possible mechanism and a few studies have shown that the success of different personality types varies with spatiotemporal conditions. However, it remains unknown whether different mechanisms can maintain personality within a species. Here we use a reintroduction programme for the critically endangered European mink (Mustela lutreola) to test whether multiple personality trait domains (boldness, exploration and sociability) affected survival in two different years and islands. This was done through pre-release personality tests and post-release radio-tracking monitoring. Survival was positively correlated with boldness, whereas the relationship with exploration was either negative or positive depending on year/island. The results show a complex relationship between personality and survival and suggest that exploration can be maintained over evolutionary time via spatiotemporal variation in conditions. However, in contrast to exploration, boldness did not vary spatiotemporally and sociability had no impact on survival. This indicates that different personality trait domains might be maintained by different mechanisms. To date, personality has been studied primarily within behavioural sciences, but through empirical findings we highlight the importance of personality also in ecology and conservation biology.

Inbreeding depression in a critically endangered carnivore.

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.

On the origin of the Norwegian lemming.

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.

Properties of vertebrate predator-prey networks in the high Arctic.

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.

Life history traits in a cyclic ecosystem: a field experiment on the arctic fox.

The reproduction of many species depends strongly on variation in food availability. The main prey of the arctic fox in Fennoscandia are cyclic small rodents, and its number of litters and litter size vary depending on the phase of the rodent cycle. In this experiment, we studied if the arctic fox adjusts its reproduction as a direct response to food abundance, in accordance with the food limitation hypothesis, or if there are additional phase-dependent trade-offs that influence its reproduction. We analysed the weaning success, i.e. proportion of arctic fox pairs established during mating that wean a litter in summer, of 422 pairs of which 361 were supplementary winter fed, as well as the weaned litter size of 203 litters of which 115 were supplementary winter fed. Females without supplementary winter food over-produced cubs in relation to food abundance in the small rodent increase phase, i.e. the litter size was equal to that in the peak phase when food was more abundant. The litter size for unfed females was 6.38 in the increase phase, 7.11 in the peak phase and 3.84 in the decrease phase. The litter size for supplementary winter-fed litters was 7.95 in the increase phase, 10.61 in the peak phase and 7.86 in the decrease phase. Thus, feeding had a positive effect on litter size, but it did not diminish the strong impact of the small rodent phase, supporting phase-dependent trade-offs in addition to food determining arctic fox reproduction.

Microsatellite genotyping reveals end-Pleistocene decline in mammoth autosomal genetic variation.

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.

Genomic Consequences of Fragmentation in the Endangered Fennoscandian Arctic Fox (Vulpes lagopus).

Accelerating climate change is causing severe habitat fragmentation in the Arctic, threatening the persistence of many cold-adapted species. The Scandinavian arctic fox (Vulpes lagopus) is highly fragmented, with a once continuous, circumpolar distribution, it struggled to recover from a demographic bottleneck in the late 19th century. The future persistence of the entire Scandinavian population is highly dependent on the northernmost Fennoscandian subpopulations (Scandinavia and the Kola Peninsula), to provide a link to the viable Siberian population. By analyzing 43 arctic fox genomes, we quantified genomic variation and inbreeding in these populations. Signatures of genome erosion increased from Siberia to northern Sweden indicating a stepping-stone model of connectivity. In northern Fennoscandia, runs of homozygosity (ROH) were on average ~1.47-fold longer than ROH found in Siberia, stretching almost entire scaffolds. Moreover, consistent with recent inbreeding, northern Fennoscandia harbored more homozygous deleterious mutations, whereas Siberia had more in heterozygous state. This study underlines the value of documenting genome erosion following population fragmentation to identify areas requiring conservation priority. With the increasing fragmentation and isolation of Arctic habitats due to global warming, understanding the genomic and demographic consequences is vital for maintaining evolutionary potential and preventing local extinctions.

Kin encounter rate and inbreeding avoidance in canids.

Mating with close kin can lead to inbreeding depression through the expression of recessive deleterious alleles and loss of heterozygosity. Mate selection may be affected by kin encounter rate, and inbreeding avoidance may not be uniform but associated with age and social system. Specifically, selection for kin recognition and inbreeding avoidance may be more developed in species that live in family groups or breed cooperatively. To test this hypothesis, we compared kin encounter rate and the proportion of related breeding pairs in noninbred and highly inbred canid populations. The chance of randomly encountering a full sib ranged between 1-8% and 20-22% in noninbred and inbred canid populations, respectively. We show that regardless of encounter rate, outside natal groups mates were selected independent of relatedness. Within natal groups, there was a significant avoidance of mating with a relative. Lack of discrimination against mating with close relatives outside packs suggests that the rate of inbreeding in canids is related to the proximity of close relatives, which could explain the high degree of inbreeding depression observed in some populations. The idea that kin encounter rate and social organization can explain the lack of inbreeding avoidance in some species is intriguing and may have implications for the management of populations at risk.

Pulses of movement across the sea ice: population connectivity and temporal genetic structure in the arctic fox.

Lemmings are involved in several important functions in the Arctic ecosystem. The Arctic fox (Vulpes lagopus) can be divided into two discrete ecotypes: "lemming foxes" and "coastal foxes". Crashes in lemming abundance can result in pulses of "lemming fox" movement across the Arctic sea ice and immigration into coastal habitats in search for food. These pulses can influence the genetic structure of the receiving population. We have tested the impact of immigration on the genetic structure of the "coastal fox" population in Svalbard by recording microsatellite variation in seven loci for 162 Arctic foxes sampled during the summer and winter over a 5-year period. Genetic heterogeneity and temporal genetic shifts, as inferred by STRUCTURE simulations and deviations from Hardy-Weinberg proportions, respectively, were recorded. Maximum likelihood estimates of movement as well as STRUCTURE simulations suggested that both immigration and genetic mixture are higher in Svalbard than in the neighbouring "lemming fox" populations. The STRUCTURE simulations and AMOVA revealed there are differences in genetic composition of the population between summer and winter seasons, indicating that immigrants are not present in the reproductive portion of the Svalbard population. Based on these results, we conclude that Arctic fox population structure varies with time and is influenced by immigration from neighbouring populations. The lemming cycle is likely an important factor shaping Arctic fox movement across sea ice and the subsequent population genetic structure, but is also likely to influence local adaptation to the coastal habitat and the prevalence of diseases.

Behavioural responses of breeding arctic sandpipers to ground-surface temperature and primary productivity.

Most birds incubate their eggs, which requires time and energy at the expense of other activities. Birds generally have two incubation strategies: biparental where both mates cooperate in incubating eggs, and uniparental where a single parent incubates. In harsh and unpredictable environments, incubation is challenging due to high energetic demands and variable resource availability. We studied the relationships between the incubation behaviour of sandpipers (genus Calidris) and two environmental variables: temperature and a proxy of primary productivity (i.e. NDVI). We investigated how these relationships vary between incubation strategies and across species among strategies. We also studied how the relationship between current temperature and incubation behaviour varies with previous day's temperature. We monitored the incubation behaviour of nine sandpiper species using thermologgers at 15 arctic sites between 2016 and 2019. We also used thermologgers to record the ground surface temperature at conspecific nest sites and extracted NDVI values from a remote sensing product. We found no relationship between either environmental variables and biparental incubation behaviour. Conversely, as ground-surface temperature increased, uniparental species decreased total duration of recesses (TDR) and mean duration of recesses (MDR), but increased number of recesses (NR). Moreover, small species showed stronger relationships with ground-surface temperature than large species. When all uniparental species were combined, an increase in NDVI was correlated with higher mean duration, total duration and number of recesses, but relationships varied widely across species. Finally, some uniparental species showed a lag effect with a higher nest attentiveness after a warm day while more recesses occurred after a cold day than was predicted based on current temperatures. We demonstrate the complex interplay between shorebird incubation strategies, incubation behaviour, and environmental conditions. Understanding how species respond to changes in their environment during incubation helps predict their future reproductive success.

Temporal genetic change in the last remaining population of woolly mammoth.

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.

Documenting lemming population change in the Arctic: Can we detect trends?

Lemmings are a key component of tundra food webs and changes in their dynamics can affect the whole ecosystem. We present a comprehensive overview of lemming monitoring and research activities, and assess recent trends in lemming abundance across the circumpolar Arctic. Since 2000, lemmings have been monitored at 49 sites of which 38 are still active. The sites were not evenly distributed with notably Russia and high Arctic Canada underrepresented. Abundance was monitored at all sites, but methods and levels of precision varied greatly. Other important attributes such as health, genetic diversity and potential drivers of population change, were often not monitored. There was no evidence that lemming populations were decreasing in general, although a negative trend was detected for low arctic populations sympatric with voles. To keep the pace of arctic change, we recommend maintaining long-term programmes while harmonizing methods, improving spatial coverage and integrating an ecosystem perspective.

Correction to: Documenting lemming population change in the Arctic: Can we detect trends?

In the original published article, some of the symbols in figure 1A were modified incorrectly during the typesetting and publication process. The correct version of the figure is provided in this correction.

Altitude effects on spatial components of vascular plant diversity in a subarctic mountain tundra.

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.