Phylogeography and post-glacial recolonization in wolverines (Gulo gulo) from across their circumpolar distribution

In Search of the Elusive North: Evolutionary History of the Arctic Fox (Vulpes lagopus) in the Palearctic from the Late Pleistocene to the Recent Inferred from Mitogenomic Data

Despite the high level of interest, the population history of arctic foxes during the Late Pleistocene and Holocene remains poorly understood. Here we aimed to fill gaps in the demographic and colonization history of the arctic fox by analyzing new ancient DNA data from fossil specimens aged from 50 to 1 thousand years from the Northern and Polar Urals, historic DNA from museum specimens from the Novaya Zemlya Archipelago and the Taymyr Peninsula and supplementing these data by previously published sequences of recent and extinct arctic foxes from other regions. This dataset was used for reconstruction of a time-calibrated phylogeny and a temporal haplotype network covering four time intervals: Late Pleistocene (ranging from 30 to 13 thousand years bp), Holocene (ranging from 4 to 1 thousand years bp), historical (approximately 150 years), and modern. Our results revealed that Late Pleistocene specimens showed no genetic similarity to either modern or historical specimens, thus supporting the earlier hypothesis on local extinction rather than habitat tracking.

Wolverines (Gulo gulo) in the Arctic: Revisiting distribution and identifying research and conservation priorities amid rapid environmental change

Wolverines (Gulo gulo) occupy most of the globe’s Arctic tundra. Given the rapidly warming climate and expanding human activity in this biome, understanding wolverine ecology, and therefore the species’ vulnerability to such changes, is increasingly important for developing research priorities and effective management strategies. Here, we review and synthesize knowledge of wolverines in the Arctic using both Western science sources and available Indigenous Knowledge (IK) to improve our understanding of wolverine ecology in the Arctic and better predict the species’ susceptibility to change. To accomplish this, we update the pan-Arctic distribution map of wolverines to account for recent observations and then discuss resulting inference and uncertainties. We use these patterns to contextualize and discuss potential underlying drivers of distribution and population dynamics, drawing upon knowledge of food habits, habitat associations, and harvest, as well as studies of wolverine ecology elsewhere. We then identify four broad areas to prioritize conservation and research efforts: (1) Monitoring trends in population abundance, demographics, and distribution and the drivers thereof, (2) Evaluating and predicting wolverines’ responses to ongoing climate change, particularly the consequences of reduced snow and sea ice, and shifts in prey availability, (3) Understanding wolverines’ response to human development, including the possible impact of wintertime over-snow travel and seismic testing to reproductive denning, as well as vulnerability to hunting and trapping associated with increased human access, and (4) Ensuring that current and future harvest are sustainable.

Dynamic landscapes in northwestern North America structured populations of wolverines (Gulo gulo)

Cyclic climatic and glacial fluctuations of the Late Quaternary produced a dynamic biogeographic history for high latitudes. To refine our understanding of this history in northwestern North America, we explored geographic structure in a wide-ranging carnivore, the wolverine (Gulo gulo). We examined genetic variation in populations across mainland Alaska, coastal Southeast Alaska, and mainland western Canada using nuclear microsatellite genotypes and sequence data from the mitochondrial DNA (mtDNA) control region and Cytochrome b (Cytb) gene. Data from maternally inherited mtDNA reflect stable populations in Northwest Alaska, suggesting the region harbored wolverine populations since at least the Last Glacial Maximum (LGM; 21 Kya), consistent with their persistence in the fossil record of Beringia. Populations in Southeast Alaska are characterized by minimal divergence, with no genetic signature of long-term refugial persistence (consistent with the lack of pre-Holocene fossil records there). The Kenai Peninsula population exhibits mixed signatures depending on marker type: mtDNA data indicate stability (i.e., historical persistence) and include a private haplotype, whereas biparentally inherited microsatellites exhibit relatively low variation and a lack of private alleles consistent with a more recent Holocene colonization of the peninsula. Our genetic work is largely consistent with the early 20th century taxonomic hypothesis that wolverines on the Kenai Peninsula belong to a distinct subspecies. Our finding of significant genetic differentiation of wolverines inhabiting the Kenai Peninsula, coupled with the peninsula’s burgeoning human population and the wolverine’s known sensitivity to anthropogenic impacts, provides valuable foundational data that can be used to inform conservation and management prescriptions for wolverines inhabiting these landscapes.

Characterizing the elusive Vancouver Island wolverine, Gulo gulo vancouverensis, using historical DNA

The wolverine (Gulo gulo) is a Holarctic species found in North America primarily across the boreal forest, the subarctic, and along the Pacific coast, including Vancouver Island (VI), British Columbia. While wolverines on VI are rare and possibly extirpated, they have been previously described as a unique subspecies, G. g. vancouverensis, distinct from G. g. luscus from the mainland of North America. However, the validity of the VI subspecies is contentious, with conflicting results from studies of skull morphology. Here, we used molecular analyses to characterize the genetic diversity of the VI population and resolve this taxonomic debate to assist with conservation priorities. Historical DNA of VI wolverines was obtained from museum specimens, amplified at 16 nuclear microsatellite loci, and sequenced at the mitochondrial D-loop control region to compare with wolverines from mainland British Columbia. The VI population had lower allelic richness and was fixed for a single common mtDNA haplotype. Bayesian and non-Bayesian assignments using microsatellites generally revealed admixture across populations, implying allele frequencies between the VI and mainland populations were not significantly different. Hence, both types of genetic markers showed little evolutionary divergence between VI and the mainland population. Combined, these results do not provide evidence of significant genetic distinction for VI wolverines, nor support the subspecific classification. Immediate conservation efforts should focus on estimating population size, while future conservation planning can assume VI wolverines likely are not a unique genetic population and there remains the potential for natural recolonization of wolverines to VI.

Population genetic structures at multiple spatial scales: importance of social groups in European badgers

Population viability and metapopulation dynamics are strongly affected by gene flow. Identifying ecological correlates of genetic structure and gene flow in wild populations is therefore a major issue both in evolutionary ecology and species management. Studying the genetic structure of populations also enables identification of the spatial scale at which most gene flow occurs, hence the scale of the functional connectivity, which is of paramount importance for species ecology. In this study, we examined the genetic structure of a social, continuously distributed mammal, the European badger (Meles meles), both at large spatial scales (among populations) and fine (within populations) spatial scales. The study was carried out in 11 sites across France utilizing a noninvasive hair trapping protocol at 206 monitored setts. We identified 264 badgers genotyped at 24 microsatellite DNA loci. At the large scale, we observed high and significant genetic differentiation among populations (global Fst = 0.139; range of pairwise Fst [0.046–0.231]) that was not related to the geographic distance among sites, suggesting few large-scale dispersal events. Within populations, we detected a threshold value below which badgers were genetically close (< 400 m), highlighting that sociality is the major structuring process within badger populations at the fine scale.

Population genetics of the wolverine in Finland: the road to recovery?

After decades, even centuries of persecution, large carnivore populations are widely recovering in Europe. Considering the recent recovery of the wolverine (Gulo gulo) in Finland, our aim was to evaluate genetic variation using 14 microsatellites and mtDNA control region (579 bp) in order (1) to determine whether the species is represented by a single genetic population within Finland, (2) to quantify the genetic diversity, and (3) to estimate the effective population size. We found two major genetic clusters divided between eastern and northern Finland based on microsatellites (FST = 0.100) but also a significant pattern of isolation by distance. Wolverines in western Finland had a genetic signature similar to the northern cluster, which can be explained by former translocations of wolverines from northern to western Finland. For both main clusters, most estimates of the effective population size Ne were below 50. Nevertheless, the genetic diversity was higher in the eastern cluster (HE = 0.57, AR = 4.0, AP = 0.3) than in the northern cluster (HE = 0.49, AR = 3.7, AP = 0.1). Migration between the clusters was low. Two mtDNA haplotypes were found: one common and identical to Scandinavian wolverines; the other rare and not previously detected. The rare haplotype was more prominent in the eastern genetic cluster. Combining all available data, we infer that the genetic population structure within Finland is shaped by a recent bottleneck, isolation by distance, human-aided translocations and postglacial recolonization routes.

Origins of Chinese reindeer (Rangifer tarandus) based on mitochondrial DNA analyses

The most southern population of reindeer (Rangifer tarandus) inhabits northeastern China, but the migration route and origin of this population have not been confirmed. The sequences of mitochondrial DNA control regions from domestic and wild herds from Eurasia and China were analysed. The results showed that the Chinese reindeer population originated independently from north-central Russian domestic herds, belonging to a large reindeer population that was present across Beringia during the last glacial period. Some studies have reported that the Chinese reindeer population is closely related to wild forest reindeer herds in Russia. Our results, however, indicate that wild forest reindeer herds of southeastern Russia contributed little or nothing to the Chinese reindeer herd gene pool. Chinese reindeer herds have a much greater genetic similarity to more northerly distributed tundra-type herds that inhabit open areas. The present findings will be essential for future conservation planning for Chinese reindeer.

Genome sequencing and conservation genomics in the Scandinavian wolverine population

Genetic approaches have proved valuable to the study and conservation of endangered populations, especially for monitoring programs, and there is potential for further developments in this direction by extending analyses to the genomic level. We assembled the genome of the wolverine (Gulo gulo), a mustelid that in Scandinavia has recently recovered from a significant population decline, and obtained a 2.42 Gb draft sequence representing >85% of the genome and including >21,000 protein-coding genes. We then performed whole-genome resequencing of 10 Scandinavian wolverines for population genomic and demographic analyses. Genetic diversity was among the lowest detected in a red-listed population (mean genome-wide nucleotide diversity of 0.05%). Results of the demographic analyses indicated a long-term decline of the effective population size (N_e ) from 10,000 well before the last glaciation to <500 after this period. Current N_e appeared even lower. The genome-wide F_IS level was 0.089 (possibly signaling inbreeding), but this effect was not observed when analyzing a set of highly variable SNP markers, illustrating that such markers can give a biased picture of the overall character of genetic diversity. We found significant population structure, which has implications for population connectivity and conservation. We used an integrated microfluidic circuit chip technology to develop an SNP-array consisting of 96 highly informative markers that, together with a multiplex pre-amplification step, was successfully applied to low-quality DNA from scat samples. Our findings will inform management, conservation, and genetic monitoring of wolverines and serve as a genomic roadmap that can be applied to other endangered species. The approach used here can be generally utilized in other systems, but we acknowledge the trade-off between investing in genomic resources and direct conservation actions.

A new species of Gulo from the Early Pliocene Gray Fossil Site (Eastern United States); rethinking the evolution of wolverines

The wolverine (Gulo gulo) is the largest living terrestrial member of the Mustelidae; a versatile predator formerly distributed throughout boreal regions of North America and Eurasia. Though commonly recovered from Pleistocene sites across their range, pre-Pleistocene records of the genus are exceedingly rare. Here, we describe a new species of Gulo from the Gray Fossil Site in Tennessee. Based on biostratigraphy, a revised estimate of the age of the Gray Fossil Site is Early Pliocene, near the Hemphillian—Blancan transition, between 4.9 and 4.5 Ma. This represents the earliest known occurrence of a wolverine, more than one million years earlier than any other record. The new species of wolverine described here shares similarities with previously described species of Gulo, and with early fishers (Pekania). As the earliest records of both Gulo and Pekania are known from North America, this suggests the genus may have evolved in North America and dispersed to Eurasia later in the Pliocene. Both fauna and flora at the Gray Fossil Site are characteristic of warm/humid climates, which suggests wolverines may have become ‘cold-adapted’ relatively recently. Finally, detailed comparison indicates Plesiogulo, which has often been suggested to be ancestral to Gulo, is not likely closely related to gulonines, and instead may represent convergence on a similar niche.

Colonizing the High Arctic: Mitochondrial DNA Reveals Common Origin of Eurasian Archipelagic Reindeer (Rangifer tarandus)

In light of current debates on global climate change it has become important to know more on how large, roaming species have responded to environmental change in the past. Using the highly variable mitochondrial control region, we revisit theories of Rangifer colonization and propose that the High Arctic archipelagos of Svalbard, Franz Josef Land, and Novaia Zemlia were colonized by reindeer from the Eurasian mainland after the last glacial maximum. Comparing mtDNA control region sequences from the three Arctic archipelagos showed a strong genetic connection between the populations, supporting a common origin in the past. A genetic connection between the three archipelagos and two Russian mainland populations was also found, suggesting colonization of the Eurasian high Arctic archipelagos from the Eurasian mainland. The age of the Franz Josef Land material (>2000 years before present) implies that Arctic indigenous reindeer colonized the Eurasian Arctic archipelagos through natural dispersal, before humans approached this region.

Wolf population genetics in Europe: a systematic review, meta-analysis and suggestions for conservation and management

The grey wolf (Canis lupus) is an iconic large carnivore that has increasingly been recognized as an apex predator with intrinsic value and a keystone species. However, wolves have also long represented a primary source of human-carnivore conflict, which has led to long-term persecution of wolves, resulting in a significant decrease in their numbers, genetic diversity and gene flow between populations. For more effective protection and management of wolf populations in Europe, robust scientific evidence is crucial. This review serves as an analytical summary of the main findings from wolf population genetic studies in Europe, covering major studies from the 'pre-genomic era' and the first insights of the 'genomics era'. We analyse, summarize and discuss findings derived from analyses of three compartments of the mammalian genome with different inheritance modes: maternal (mitochondrial DNA), paternal (Y chromosome) and biparental [autosomal microsatellites and single nucleotide polymorphisms (SNPs)]. To describe large-scale trends and patterns of genetic variation in European wolf populations, we conducted a meta-analysis based on the results of previous microsatellite studies and also included new data, covering all 19 European countries for which wolf genetic information is available: Norway, Sweden, Finland, Estonia, Latvia, Lithuania, Poland, Czech Republic, Slovakia, Germany, Belarus, Russia, Italy, Croatia, Bulgaria, Bosnia and Herzegovina, Greece, Spain and Portugal. We compared different indices of genetic diversity in wolf populations and found a significant spatial trend in heterozygosity across Europe from south-west (lowest genetic diversity) to north-east (highest). The range of spatial autocorrelation calculated on the basis of three characteristics of genetic diversity was 650-850 km, suggesting that the genetic diversity of a given wolf population can be influenced by populations up to 850 km away. As an important outcome of this synthesis, we discuss the most pressing issues threatening wolf populations in Europe, highlight important gaps in current knowledge, suggest solutions to overcome these limitations, and provide recommendations for science-based wolf conservation and management at regional and Europe-wide scales.

A test of somatic mosaicism in the androgen receptor gene of Canada lynx (Lynx canadensis)

Background

The androgen receptor, an X-linked gene, has been widely studied in human populations because it contains highly polymorphic trinucleotide repeat motifs that have been associated with a number of adverse human health and behavioral effects. A previous study on the androgen receptor gene in carnivores reported somatic mosaicism in the tissues of a number of species including Eurasian lynx (Lynx lynx). We investigated this claim in a closely related species, Canada lynx (Lynx canadensis). The presence of somatic mosaicism in lynx tissues could have implications for the future study of exonic trinucleotide repeats in landscape genomic studies, in which the accurate reporting of genotypes would be highly problematic.

Methods

To determine whether mosaicism occurs in Canada lynx, two lynx individuals were sampled for a variety of tissue types (lynx 1) and tissue locations (lynx 1 and 2), and 1,672 individuals of known sex were genotyped to further rule out mosaicism.

Results

We found no evidence of mosaicism in tissues from the two necropsied individuals, or any of our genotyped samples.

Conclusions

Our results indicate that mosaicism does not manifest in Canada lynx. Therefore, the use of hide samples for further work involving trinucleotide repeat polymorphisms in Canada lynx is warranted.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0284-y) contains supplementary material, which is available to authorized users.

Lack of Spatial Immunogenetic Structure among Wolverine (Gulo gulo) Populations Suggestive of Broad Scale Balancing Selection

Elucidating the adaptive genetic potential of wildlife populations to environmental selective pressures is fundamental for species conservation. Genes of the major histocompatibility complex (MHC) are highly polymorphic, and play a key role in the adaptive immune response against pathogens. MHC polymorphism has been linked to balancing selection or heterogeneous selection promoting local adaptation. However, spatial patterns of MHC polymorphism are also influenced by gene flow and drift. Wolverines are highly vagile, inhabiting varied ecoregions that include boreal forest, taiga, tundra, and high alpine ecosystems. Here, we investigated the immunogenetic variation of wolverines in Canada as a surrogate for identifying local adaptation by contrasting the genetic structure at MHC relative to the structure at 11 neutral microsatellites to account for gene flow and drift. Evidence of historical positive selection was detected at MHC using maximum likelihood codon-based methods. Bayesian and multivariate cluster analyses revealed weaker population genetic differentiation at MHC relative to the increasing microsatellite genetic structure towards the eastern wolverine distribution. Mantel correlations of MHC against geographical distances showed no pattern of isolation by distance (IBD: r = -0.03, p = 0.9), whereas for microsatellites we found a relatively strong and significant IBD (r = 0.54, p = 0.01). Moreover, we found a significant correlation between microsatellite allelic richness and the mean number of MHC alleles, but we did not observe low MHC diversity in small populations. Overall these results suggest that MHC polymorphism has been influenced primarily by balancing selection and to a lesser extent by neutral processes such as genetic drift, with no clear evidence for local adaptation. This study contributes to our understanding of how vulnerable populations of wolverines may respond to selective pressures across their range.