Chromosome-level assembly of the gray fox (Urocyon cinereoargenteus) confirms the basal loss of PRDM9 in Canidae.

Reference genome assemblies have been created from multiple lineages within the Canidae family; however, despite its phylogenetic relevance as a basal genus within the clade, there is currently no reference genome for the gray fox (Urocyon cinereoargenteus). Here, we present a chromosome-level assembly for the gray fox (U. cinereoargenteus), which represents the most contiguous, non-domestic canid reference genome available to date, with 90% of the genome contained in just 34 scaffolds and a contig N50 and scaffold N50 of 59.4 and 72.9 Megabases, respectively. Repeat analyses identified an increased number of simple repeats relative to other canids. Based on mitochondrial DNA, our Vermont sample clusters with other gray fox samples from the northeastern United States and contains slightly lower levels of heterozygosity than gray foxes on the west coast of California. This new assembly lays the groundwork for future studies to describe past and present population dynamics, including the delineation of evolutionarily significant units of management relevance. Importantly, the phylogenetic position of Urocyon allows us to verify the loss of PRDM9 functionality in the basal canid lineage, confirming that pseudogenization occurred at least 10 million years ago.

Molecular epidemiology of Enterocytozoon bieneusi from foxes and raccoon dogs in the Henan and Hebei provinces in China.

BACKGROUND: Enterocytozoon bieneusi is a zoonotic pathogen widely distributed in animals and humans. It can cause diarrhea and even death in immunocompromised hosts. Approximately 800 internal transcribed spacer (ITS) genotypes have been identified in E. bieneusi. Farmed foxes and raccoon dogs are closely associated to humans and might be the reservoir of E. bieneusi which is known to have zoonotic potential. However, there are only a few studies about E. bieneusi genotype identification and epidemiological survey in foxes and raccoon dogs in Henan and Hebei province. Thus, the present study investigated the infection rates and genotypes of E. bieneusi in farmed foxes and raccoon dogs in the Henan and Hebei provinces. RESULT: A total of 704 and 884 fecal specimens were collected from foxes and raccoon dogs, respectively. Nested PCR was conducted based on ITS of ribosomal RNA (rRNA), and then multilocus sequence typing (MLST) was conducted to analyze the genotypes. The result showed that infection rates of E. bieneusi in foxes and raccoon dogs were 18.32% and 5.54%, respectively. Ten E. bieneusi genotypes with zoonotic potential (NCF2, NCF3, D, EbpC, CHN-DC1, SCF2, CHN-F1, Type IV, BEB4, and BEB6) were identified in foxes and raccoon dogs. Totally 178 ITS-positive DNA specimens were identified from foxes and raccoon dogs and these specimens were then subjected to MLST analysis. In the MLST analysis, 12, 2, 7 and 8 genotypes were identified in at the mini-/ micro-satellite loci MS1, MS3, MS4 and MS7, respectively. A total of 14 multilocus genotypes were generated using ClustalX 2.1 software. Overall, the present study evaluated the infection of E. bieneusi in foxes and raccoon dogs in the Henan and Hebei province, and investigated the zoonotic potential of the E. bieneusi in foxes and raccoon dogs. CONCLUSIONS: These findings expand the geographic distribution information of E. bieneusi' host in China and was helpful in preventing against the infection of E. bieneusi with zoonotic potential in foxes and raccoon dogs.

Serological and molecular survey of canine distemper virus in red foxes (Vulpes vulpes): Exploring cut-off values and the use of protein A in ELISA tests.

The wide distribution and ecological plasticity of the red fox (Vulpes vulpes) make it a potential reservoir for many infectious diseases shared with domestic and wild carnivores. One of such diseases is canine distemper, which is caused by an RNA virus and its main domestic reservoir is the dog. However, other carnivores can also participate in its maintenance, as shown by the recent upsurge of reported cases in wildlife in many parts of the world, and by the fact that red foxes may act as true reservoirs for canine distemper virus (CDV). The lack of validated serological tests for wildlife or other non-target species may be a handicap for monitoring this virus. In this study, serological assays were compared in 147 red fox sera using a commercial ELISA validated for its use in dogs and a non-specific modified ELISA with Protein A peroxidase conjugate to detect bound antibodies. In addition, the presence of CDV RNA in brain, spleen, lung, and liver samples from 144 foxes was investigated by a RT-qPCR. Through the comparison of the results of both ELISAs and the use of a finite mixture model of the optical density values obtained by both techniques, we adjusted the cut-off point of the commercial ELISA to obtain the seroprevalence in foxes. The overall seroprevalence detected was 53.7% (79/147) and 57.1% (84/147) by the commercial and modified ELISA, respectively, with a moderate agreement according to Cohen's Kappa statistic (κ = 0.491, z = 5.97, p < 0.0001). CDV RNA was detected in 30 out of 144 foxes, which resulted in 20.8% of CDV-infected foxes. At individual level, the results obtained by relating the serological status and the presence/absence of RNA in different organs were explained in terms of the pathogenesis of the infection. Our results highlight the convenience of adjusting the cut-off point when using an ELISA assay developed in domestic dogs for its use in foxes. Moreover, Protein A is confirmed to be a good alternative to be used in red foxes, presenting a good reactivity towards its IgG.

Can demographic histories explain long-term isolation and recent pulses of asymmetric gene flow between highly divergent grey fox lineages?

Secondary contact zones between deeply divergent, yet interfertile, lineages provide windows into the speciation process. North American grey foxes (Urocyon cinereoargenteus) are divided into western and eastern lineages that diverged approximately 1 million years ago. These ancient lineages currently hybridize in a relatively narrow zone of contact in the southern Great Plains, a pattern more commonly observed in smaller-bodied taxa, which suggests relatively recent contact after a long period of allopatry. Based on local ancestry inference with whole-genome sequencing (n = 43), we identified two distinct Holocene pulses of admixture. The older pulse (500-3500 YBP) reflected unidirectional gene flow from east to west, whereas the more recent pulse (70-200 YBP) of admixture was bi-directional. Augmented with genotyping-by-sequencing data from 216 additional foxes, demographic analyses indicated that the eastern lineage declined precipitously after divergence, remaining small throughout most of the late Pleistocene, and expanding only during the Holocene. Genetic diversity in the eastern lineage was highest in the southeast and lowest near the contact zone, consistent with a westward expansion. Concordantly, distribution modelling indicated that during their isolation, the most suitable habitat occurred far east of today's contact zone or west of the Great Plains. Thus, long-term isolation was likely caused by the small, distant location of the eastern refugium, with recent contact reflecting a large increase in suitable habitat and corresponding demographic expansion from the eastern refugium. Ultimately, long-term isolation in grey foxes may reflect their specialized bio-climatic niche. This system presents an opportunity for future investigation of potential pre- and post-zygotic isolating mechanisms.

High-Quality Chromosome-Level Genome Assembly of the Corsac Fox (Vulpes corsac) Reveals Adaptation to Semiarid and Harsh Environments.

The Corsac fox (Vulpes corsac) is a species of fox distributed in the arid prairie regions of Central and Northern Asia, with distinct adaptations to dry environments. Here, we applied Oxford-Nanopore sequencing and a chromosome structure capture technique to assemble the first Corsac fox genome, which was then assembled into chromosome fragments. The genome assembly has a total length of 2.2 Gb with a contig N50 of 41.62 Mb and a scaffold N50 of 132.2 Mb over 18 pseudo-chromosomal scaffolds. The genome contained approximately 32.67% of repeat sequences. A total of 20,511 protein-coding genes were predicted, of which 88.9% were functionally annotated. Phylogenetic analyses indicated a close relation to the Red fox (Vulpes vulpes) with an estimated divergence time of ~3.7 million years ago (MYA). We performed separate enrichment analyses of species-unique genes, the expanded and contracted gene families, and positively selected genes. The results suggest an enrichment of pathways related to protein synthesis and response and an evolutionary mechanism by which cells respond to protein denaturation in response to heat stress. The enrichment of pathways related to lipid and glucose metabolism, potentially preventing stress from dehydration, and positive selection of genes related to vision, as well as stress responses in harsh environments, may reveal adaptive evolutionary mechanisms in the Corsac fox under harsh drought conditions. Additional detection of positive selection for genes associated with gustatory receptors may reveal a unique desert diet strategy for the species. This high-quality genome provides a valuable resource for studying mammalian drought adaptation and evolution in the genus Vulpes.

Low levels of genetic differentiation and structure in red fox populations in Eastern Canada and implications for Arctic fox rabies propagation potential.

Rabies is a lethal zoonosis present in most parts of the world which can be transmitted to humans through the bite from an infected mammalian reservoir host. The Arctic rabies virus variant (ARVV) persists mainly in populations of Arctic foxes (Vulpes lagopus), and to a lesser extent in red fox populations (Vulpes vulpes). Red foxes are thought to be responsible for sporadic southward movement waves of the ARVV outside the enzootic area of northern Canada. In this study, we wanted to investigate whether red foxes displayed notable levels of genetic structure across the Quebec-Labrador Peninsula, which includes portions of the provinces of Quebec and Newfoundland-Labrador in Canada, and is a region with a history of southward ARVV movement waves. We combined two datasets that were collected and genotyped using different protocols, totalling 675 red fox individuals across the whole region and genotyped across 13 microsatellite markers. We found two genetic clusters across the region, reflecting a latitudinal gradient, and characterized by low genetic differentiation. We also observed weak but significant isolation by distance, which seems to be marginally more important for females than for males. These findings suggest a general lack of resistance to movement in red fox populations across the Quebec-Labrador Peninsula, regardless of sex. Implications of these findings include additional support for the hypothesis of long-distance southward ARVV propagation through its red fox reservoir host.

North African fox genomes show signatures of repeated introgression and adaptation to life in deserts.

Elucidating the evolutionary process of animal adaptation to deserts is key to understanding adaptive responses to climate change. Here we generated 82 individual whole genomes of four fox species (genus Vulpes) inhabiting the Sahara Desert at different evolutionary times. We show that adaptation of new colonizing species to a hot arid environment has probably been facilitated by introgression and trans-species polymorphisms shared with older desert resident species, including a putatively adaptive 25 Mb genomic region. Scans for signatures of selection implicated genes affecting temperature perception, non-renal water loss and heat production in the recent adaptation of North African red foxes (Vulpes vulpes), after divergence from Eurasian populations approximately 78 thousand years ago. In the extreme desert specialists, Rueppell's fox (V. rueppellii) and fennec (V. zerda), we identified repeated signatures of selection in genes affecting renal water homeostasis supported by gene expression and physiological differences. Our study provides insights into the mechanisms and genetic underpinnings of a natural experiment of repeated adaptation to extreme conditions.

Genomic recovery lags behind demographic recovery in bottlenecked populations of the Channel Island fox, Urocyon littoralis.

With continued global change, recovery of species listed under the Endangered Species Act is increasingly challenging. One rare success was the recovery and delisting of the Channel Island fox (Urocyon littoralis) after 90%-99% population declines in the 1990s. While their demographic recovery was marked, less is known about their genetic recovery. To address genetic changes, we conducted the first multi-individual and population-level direct genetic comparison of samples collected before and after the recent bottlenecks. Using whole-exome sequencing, we found that already genetically depauperate populations were further degraded by the 1990s declines and remain low, particularly on San Miguel and Santa Rosa Islands, which underwent the most severe bottlenecks. The two other islands that experienced recent bottlenecks (Santa Cruz, and Santa Catalina islands) showed mixed results based on multiple metrics of genetic diversity. Previous island fox genomics studies showed low genetic diversity before the declines and no change after the demographic recovery, thus this is the first study to show a decrease in genetic diversity over time in U. littoralis. Additionally, we found that divergence between populations consistently increased over time, complicating prospects for using inter-island translocation as a conservation tool. The Santa Catalina subspecies is now federally listed as threatened, yet other de-listed subspecies are still recovering genetic variation which may limit their ability to adapt to changing environmental conditions. This study further demonstrates that species conservation is more complex than population size and that some island fox populations are not yet 'out of the woods'.

Synanthropic and Wild Animals as Sentinels of Zoonotic Agents: A Study of Leptospira Genotypes Circulating in Northeastern Italy.

Leptospirosis is an infectious disease widely reported in veterinary practice and a worldwide zoonosis. In Northeastern Italy, different serogroups and genotypes of Leptospira have been described in ill dogs, the most commonly detected being Icterohaemorragiae (ICT) ST 17, Australis (AUS) ST 24 and ST 198, Pomona (POM) ST 117 and ST 289, and Sejroe (SEJ) ST 155. However, there is little information available on the environmental exposure to Leptospira of wild and synanthropic animals. The aim of this study was to identify the circulating genotypes in potential reservoirs to fill this gap of knowledge. Between 2015 and 2022, 681 animal carcasses collected by the Public Veterinary Service were analyzed for Leptospira with a real-time PCR-based screening test, while positive samples were genotyped by multi-locus sequence typing analysis. To carry out our study, we tested 330 hedgehogs, 105 red foxes, 108 Norway rats, 79 mice, 22 coypus, 10 bank voles, 13 grey wolves, 5 common shrews and 9 greater mouse-eared bats. Five sequence types (STs) common in dogs were also found in wild animals: ST 24, ST 198, ST 17 and ST 155 in hedgehogs, ST 17 and ST 24 in foxes, ST 17 in rats, ST 17 and ST 155 in mice, and ST 117 in a wolf. In addition, to the best of the authors' knowledge, this is the first Italian report of SEJ ST 197 in a bank vole. Furthermore, this study described a previous survey conducted in 2009 on coypus (30 animals from the province of Trento and 41 from the province of Padua), referring to a serological positivity (L. Bratislava) without any molecular detection of Leptospira. This study on Leptospira in synanthropic and wild animals highlighted the importance of increasing our epidemiological knowledge of leptospirosis and its zoonotic risks.

Genomic analyses of gray fox lineages suggest ancient divergence and secondary contact in the southern Great Plains.

The gray fox (Urocyon cinereoargenteus) lineage diverged from all other extant canids at their most basal node and is restricted to the Americas. Previous mitochondrial analysis from coastal populations identified deeply divergent (up to 1 Mya) eastern and western lineages that predate most intraspecific splits in carnivores. We conducted genotyping by sequencing and mitochondrial analysis on gray foxes sampled across North America to determine geographic concordance between nuclear and mitochondrial contact zones and divergence times. We also estimated the admixture within the contact zone between eastern and western gray foxes based on nuclear DNA. Both datasets confirmed that eastern and western lineages met in the southern Great Plains (i.e. Texas and Oklahoma), where they maintained high differentiation. Admixture was generally low, with the majority of admixed individuals carrying <10% ancestry from the other lineage. Divergence times confirmed a mid-Pleistocene split, similar to the mitochondrial estimates. Taken together, findings suggest gray fox lineages represent an ancient divergence event, far older than most intraspecific divergences in North American carnivores. Low admixture may reflect a relatively recent time since secondary contact (e.g. post-Pleistocene) or, alternatively, ecological or reproductive barriers between lineages. Though further research is needed to disentangle these factors, our genomic investigation suggests species-level divergence exists between eastern and western gray fox lineages.

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.

Population Genetics of California Gray Foxes Clarify Origins of the Island Fox.

We used mitochondrial sequences and nuclear microsatellites to investigate population structure of gray foxes (Urocyon cinereoargenteus) and the evolutionary origins of the endemic island fox (Urocyon littoralis), which first appeared in the northern Channel Islands <13,000 years ago and in the southern Channel Islands <6000 years ago. It is unclear whether island foxes evolved directly from mainland gray foxes transported to the islands one or more times or from a now-extinct mainland population, already diverged from the gray fox. Our 345 mitochondrial sequences, combined with previous data, confirmed island foxes to be monophyletic, tracing to a most recent common ancestor approximately 85,000 years ago. Our rooted nuclear DNA tree additionally indicated genome-wide monophyly of island foxes relative to western gray foxes, although we detected admixture in northern island foxes from adjacent mainland gray foxes, consistent with some historical gene flow. Southern California gray foxes also bore a genetic signature of admixture and connectivity to a desert population, consistent with partial replacement by a late-Holocene range expansion. Using our outgroup analysis to root previous nuclear sequence-based trees indicated reciprocal monophyly of northern versus southern island foxes. Results were most consistent with island fox origins through multiple introductions from a now-extirpated mainland population.

On taming the effect of transcript level intra-condition count variation during differential expression analysis: A story of dogs, foxes and wolves.

The evolution of RNA-seq technologies has yielded datasets of scientific value that are often generated as condition associated biological replicates within expression studies. With expanding data archives opportunity arises to augment replicate numbers when conditions of interest overlap. Despite correction procedures for estimating transcript abundance, a source of ambiguity is transcript level intra-condition count variation; as indicated by disjointed results between analysis tools. We present TVscript, a tool that removes reference-based transcripts associated with intra-condition count variation above specified thresholds and we explore the effects of such variation on differential expression analysis. Initially iterative differential expression analysis involving simulated counts, where levels of intra-condition variation and sets of over represented transcripts are explicitly specified, was performed. Then counts derived from inter- and intra-study data representing brain samples of dogs, wolves and foxes (wolves vs. dogs and aggressive vs. tame foxes) were used. For simulations, the sensitivity in detecting differentially expressed transcripts increased after removing hyper-variable transcripts, although at levels of intra-condition variation above 5% detection became unreliable. For real data, prior to applying TVscript, ≈20% of the transcripts identified as being differentially expressed were associated with high levels of intra-condition variation, an over representation relative to the reference set. As transcripts harbouring such variation were removed pre-analysis, a discordance from 26 to 40% in the lists of differentially expressed transcripts is observed when compared to those obtained using the non-filtered reference. The removal of transcripts possessing intra-condition variation values within (and above) the 97th and 95th percentiles, for wolves vs. dogs and aggressive vs. tame foxes, maximized the sensitivity in detecting differentially expressed transcripts as a result of alterations within gene-wise dispersion estimates. Through analysis of our real data the support for seven genes with potential for being involved with selection for tameness is provided. TVscript is available at: https://sourceforge.net/projects/tvscript/.

Genomic Assessment of Cancer Susceptibility in the Threatened Catalina Island Fox (Urocyon littoralis catalinae).

Small effective population sizes raise the probability of extinction by increasing the frequency of potentially deleterious alleles and reducing fitness. However, the extent to which cancers play a role in the fitness reduction of genetically depauperate wildlife populations is unknown. Santa Catalina island foxes (Urocyon littoralis catalinae) sampled in 2007-2008 have a high prevalence of ceruminous gland tumors, which was not detected in the population prior to a recent bottleneck caused by a canine distemper epidemic. The disease appears to be associated with inflammation from chronic ear mite (Otodectes) infections and secondary elevated levels of Staphyloccus pseudointermedius bacterial infections. However, no other environmental factors to date have been found to be associated with elevated cancer risk in this population. Here, we used whole genome sequencing of the case and control individuals from two islands to identify candidate loci associated with cancer based on genetic divergence, nucleotide diversity, allele frequency spectrum, and runs of homozygosity. We identified several candidate loci based on genomic signatures and putative gene functions, suggesting that cancer susceptibility in this population may be polygenic. Due to the efforts of a recovery program and weak fitness effects of late-onset disease, the population size has increased, which may allow selection to be more effective in removing these presumably slightly deleterious alleles. Long-term monitoring of the disease alleles, as well as overall genetic diversity, will provide crucial information for the long-term persistence of this threatened population.

Melanocortin-1 receptor mutations and pigmentation: Insights from large animals.

The melanocortin-1 receptor (MC1R) is a G protein-coupled receptor expressed in cutaneous and hair follicle melanocytes, and plays a central role in coat color determination in vertebrates. Numerous MC1R variants have been identified in diverse species. Some of these variants have been associated with specific hair and skin color phenotypes in humans as well as coat color in animals. Gain-of-function mutations of the MC1R gene cause dominant or partially dominant black/dark coat color, and loss-of-function mutations of the MC1R gene cause recessive or partially recessive red/yellow/pale coat color phenotypes. These have been well documented in a large number of mammals, including human, dog, cattle, horse, sheep, pig, and fox. Higher similarities between large mammals and humans makes them better models to understand pathogenesis of human diseases caused by MC1R mutations. High identities in MC1Rs and similar variants identified in both humans and large mammals also provide an opportunity for receptor structure and function study. In this review, we aim to summarize the naturally occurring mutations of MC1R in humans and large animals.

Contrasting genetic trajectories of endangered and expanding red fox populations in the western U.S.

As anthropogenic disturbances continue to drive habitat loss and range contractions, the maintenance of evolutionary processes will increasingly require targeting measures to the population level, even for common and widespread species. Doing so requires detailed knowledge of population genetic structure, both to identify populations of conservation need and value, as well as to evaluate suitability of potential donor populations. We conducted a range-wide analysis of the genetic structure of red foxes in the contiguous western U.S., including a federally endangered distinct population segment of the Sierra Nevada subspecies, with the objectives of contextualizing field observations of relative scarcity in the Pacific mountains and increasing abundance in the cold desert basins of the Intermountain West. Using 31 autosomal microsatellites, along with mitochondrial and Y-chromosome markers, we found that populations of the Pacific mountains were isolated from one another and genetically depauperate (e.g., estimated Ne range = 3-9). In contrast, red foxes in the Intermountain regions showed relatively high connectivity and genetic diversity. Although most Intermountain red foxes carried indigenous western matrilines (78%) and patrilines (85%), the presence of nonindigenous haplotypes at lower elevations indicated admixture with fur-farm foxes and possibly expanding midcontinent populations as well. Our findings suggest that some Pacific mountain populations could likely benefit from increased connectivity (i.e., genetic rescue) but that nonnative admixture makes expanding populations in the Intermountain basins a non-ideal source. However, our results also suggest contact between Pacific mountain and Intermountain basin populations is likely to increase regardless, warranting consideration of risks and benefits of proactive measures to mitigate against unwanted effects of Intermountain gene flow.

An efficient method for simultaneous species, individual, and sex identification via in-solution single nucleotide polymorphism capture from low-quality scat samples.

Understanding predator population dynamics is important for conservation management because of the critical roles predators play within ecosystems. Noninvasive genetic sampling methods are useful for the study of predators like canids that can be difficult to capture or directly observe. Here, we introduce the FAECES* method (Fast and Accurate Enrichment of Canid Excrement for Species* and other analyses) which expands the toolbox for canid researchers and conservationists by using in-solution hybridization sequence capture to produce single nucleotide polymorphism (SNP) genotypes for multiple canid species from scat-derived DNA using a single enrichment. We designed a set of hybridization probes to genotype both coyotes (Canis latrans) and kit foxes (Vulpes macrotis) at hundreds of polymorphic SNP loci and we tested the probes on both tissues and field-collected scat samples. We enriched and genotyped by sequencing 52 coyote and 70 kit fox scats collected in and around a conservation easement in the Nevada Mojave Desert. We demonstrate that the FAECES* method produces genotypes capable of differentiating coyotes and kit foxes, identifying individuals and their sex, and estimating genetic diversity and effective population sizes, even using highly degraded, low-quantity DNA extracted from scat. We found that the study area harbours a large and diverse population of kit foxes and a relatively smaller population of coyotes. By replicating our methods in the future, conservationists can assess the impacts of management decisions on canid populations. The method can also be adapted and applied more broadly to enrich and sequence multiple loci from any species of interest using scat or other noninvasive genetic samples.

Estimating red fox density using non-invasive genetic sampling and spatial capture-recapture modelling.

Spatial capture-recapture modelling (SCR) is a powerful tool for estimating density, population size, and space use of elusive animals. Here, we applied SCR modelling to non-invasive genetic sampling (NGS) data to estimate red fox (Vulpes vulpes) densities in two areas of boreal forest in central (2016-2018) and southern Norway (2017-2018). Estimated densities were overall lower in the central study area (mean = 0.04 foxes per km2 in 2016, 0.10 in 2017, and 0.06 in 2018) compared to the southern study area (0.16 in 2017 and 0.09 in 2018). We found a positive effect of forest cover on density in the central, but not the southern study area. The absence of an effect in the southern area may reflect a paucity of evidence caused by low variation in forest cover. Estimated mean home-range size in the central study area was 45 km2 [95%CI 34-60] for females and 88 km2 [69-113] for males. Mean home-range sizes were smaller in the southern study area (26 km2 [16-42] for females and 56 km2 [35-91] for males). In both study areas, detection probability was session-dependent and affected by sampling effort. This study highlights how SCR modelling in combination with NGS can be used to efficiently monitor red fox populations, and simultaneously incorporate ecological factors and estimate their effects on population density and space use.

Next-generation phylogeography resolves post-glacial colonization patterns in a widespread carnivore, the red fox (Vulpes vulpes), in Europe.

Carnivores tend to exhibit a lack of (or less pronounced) genetic structure at continental scales in both a geographic and temporal sense and this can confound the identification of post-glacial colonization patterns in this group. In this study we used genome-wide data (using genotyping by sequencing [GBS]) to reconstruct the phylogeographic history of a widespread carnivore, the red fox (Vulpes vulpes), by investigating broad-scale patterns of genomic variation, differentiation and admixture amongst contemporary populations in Europe. Using 15,003 single nucleotide polymorphisms (SNPs) from 524 individuals allowed us to identify the importance of refugial regions for the red fox in terms of endemism (e.g., Iberia). In addition, we tested multiple post-glacial recolonization scenarios of previously glaciated regions during the Last Glacial Maximum using an Approximate Bayesian Computation (ABC) approach that were unresolved from previous studies. This allowed us to identify the role of admixture from multiple source population post-Younger Dryas in the case of Scandinavia and ancient land-bridges in the colonization of the British Isles. A natural colonization of Ireland was deemed more likely than an ancient human-mediated introduction as has previously been proposed and potentially points to a larger mammalian community on the island in the early post-glacial period. Using genome-wide data has allowed us to tease apart broad-scale patterns of structure and diversity in a widespread carnivore in Europe that was not evident from using more limited marker sets and provides a foundation for next-generation phylogeographic studies in other non-model species.