Noninvasive Sampling Reveals Short-Term Genetic Rescue in an Insular Red Fox Population

Genetic rescue often leads to higher fitness as a result of increased heterozygosity across animal taxa

Biodiversity loss has reached critical levels partly due to anthropogenic habitat loss and degradation. These landscape changes are damaging as they can fragment species distributions into small, isolated populations, resulting in limited gene flow, population declines and reduced adaptive potential. Genetic rescue, the translocation of individuals to increase genetic diversity and ultimately fitness, has produced promising results for fragmented populations but remains underutilized due to a lack of long-term data and monitoring. To promote a better understanding of genetic rescue and its potential risks and benefits over the short-term, we reviewed and analysed published genetic rescue attempts to identify whether genetic diversity increases following translocation, and if this change is associated with increased fitness. Our review identified 19 studies that provided genetic and fitness data from before and after the translocation; the majority of these were on mammals, and included experimental, natural and conservation-motivated translocations. Using a Bayesian meta-analytical approach, we found that on average, genetic diversity and fitness increased in populations post translocations, although there were some exceptions to this trend. Overall, genetic diversity was a positive predictor of population fitness, and in some cases this relationship extended three generations post-rescue. These data suggest a single translocation can have lasting fitness benefits, and support translocation as another tool to facilitate conservation success. Given the limited number of studies with long-term data, we echo the need for genetic monitoring of populations post-translocation to understand whether genetic rescue can also limit the loss of adaptive potential in the long-term.

Whole Genomes Inform Genetic Rescue Strategy for Montane Red Foxes in North America

A few iconic examples have proven the value of facilitated gene flow for counteracting inbreeding depression and staving off extinction; yet, the practice is often not implemented for fear of causing outbreeding depression. Using genomic sequencing, climatic niche modeling, and demographic reconstruction, we sought to assess the risks and benefits of using translocations as a tool for recovery of endangered montane red fox (Vulpes vulpes) populations in the western United States. We demonstrated elevated inbreeding and homozygosity of deleterious alleles across all populations, but especially those isolated in the Cascade and Sierra Nevada ranges. Consequently, translocations would be expected to increase population growth by masking deleterious recessive alleles. Demographic reconstructions further indicated shallow divergences of less than a few thousand years among montane populations, suggesting low risk of outbreeding depression. These genomic-guided findings set the stage for future management, the documentation of which will provide a roadmap for recovery of other data-deficient taxa.

Missing history of a modern domesticate: Historical demographics and genetic diversity in farm-bred red fox populations

The first record of captive-bred red foxes (Vulpes vulpes) dates to 1896 when a breeding enterprise emerged in the provinces of Atlantic Canada. Because its domestication happened during recent history, the red fox offers a unique opportunity to examine the genetic diversity of an emerging domesticated species in the context of documented historical and economic influences. In particular, the historical record suggests that North American and Eurasian farm-bred populations likely experienced different demographic trajectories. Here, we focus on the likely impacts of founder effects and genetic drift given historical trends in fox farming on North American and Eurasian farms. A total of 15 mitochondrial haplotypes were identified in 369 foxes from 10 farm populations that we genotyped (n = 161) or that were previously published. All haplotypes are endemic to North America. Although most haplotypes were consistent with eastern Canadian ancestry, a small number of foxes carried haplotypes typically found in Alaska and other regions of western North America. The presence of these haplotypes supports historical reports of wild foxes outside of Atlantic Canada being introduced into the breeding stock. These putative Alaskan and Western haplotypes were more frequently identified in Eurasian farms compared to North American farms, consistent with historical documentation suggesting that Eurasian economic and breeding practices were likely to maintain low-frequency haplotypes more effectively than in North America. Contextualizing inter- vs. intra-farm genetic diversity alongside the historical record is critical to understanding the origins of this emerging domesticate and the relationships between wild and farm-bred fox populations.

Maintenance of a narrow hybrid zone between native and introduced red foxes (Vulpes vulpes) despite conspecificity and high dispersal capabilities

Human-facilitated introductions of nonnative populations can lead to secondary contact between allopatric lineages, resulting in lineage homogenisation or the formation of stable hybrid zones maintained by reproductive barriers. We investigated patterns of gene flow between the native Sacramento Valley red fox (Vulpes vulpes patwin) and introduced conspecifics of captive-bred origin in California's Central Valley. Considering their recent divergence (20-70 kya), we hypothesised that any observed barriers to gene flow were primarily driven by pre-zygotic (e.g. behavioural differences) rather than post-zygotic (e.g. reduced hybrid fitness) barriers. We also explored whether nonnative genes could confer higher fitness in the human-dominated landscape resulting in selective introgression into the native population. Genetic analysis of red foxes (n = 682) at both mitochondrial (cytochrome b + D-loop) and nuclear (19,051 SNPs) loci revealed narrower cline widths than expected under a simulated model of unrestricted gene flow, consistent with the existence of reproductive barriers. We identified several loci with reduced introgression that were previously linked to behavioural divergence in captive-bred and domestic canids, supporting pre-zygotic, yet possibly hereditary, barriers as a mechanism driving the narrowness and stability of the hybrid zone. Several loci with elevated gene flow from the nonnative into the native population were linked to genes associated with domestication and adaptation to human-dominated landscapes. This study contributes to our understanding of hybridisation dynamics in vertebrates, particularly in the context of species introductions and landscape changes, underscoring the importance of considering how multiple mechanisms may be maintaining lineages at the species and subspecies level.

Nutritional ecology of a prototypical generalist predator, the red fox (Vulpes vulpes)

Generalist species, which exploit a wide range of food resources, are expected to be able to combine available resources as to attain their specific macronutrient ratio (percentage of caloric intake of protein, lipids and carbohydrates). Among mammalian predators, the red fox Vulpes vulpes is a widespread, opportunistic forager: its diet has been largely studied, outlining wide variation according to geographic and climatic factors. We aimed to check if, throughout the species' European range, diets vary widely in macronutrient composition or foxes can combine complementary foods to gain the same nutrient intake. First, we assessed fox's intake target in the framework of nutritional geometry. Secondly, we aimed to highlight the effects of unbalanced diets on fox density, which was assumed as a proxy for Darwinian fitness, as assessed in five areas of the western Italian Alps. Unexpectedly, the target macronutrient ratio of the fox (52.4% protein-, 38.7% lipid- and 8.9% carbohydrate energy) was consistent with that of hypercarnivores, such as wolves and felids, except for carbohydrate intakes in urban and rural habitats. The inverse relation between density and the deviation of observed macronutrient ratios from the intake target suggests that fox capability of surviving in a wide range of habitats may not be exempt from fitness costs and that nutrient availability should be regarded among the biotic factors affecting animal abundance and distribution.

A hierarchical modeling approach to predict the distribution and density of Sierra Nevada Red Fox (Vulpes vulpes necator)

Carnivores play critical roles in ecosystems, yet many species are declining worldwide. The Sierra Nevada Red Fox (Vulpes vulpes necator; SNRF) is a rare and endangered subspecies of red fox limited to upper montane forests, subalpine, and alpine environments of California and Oregon, United States. Having experienced significant distribution contractions and population declines in the last century, the subspecies is listed as at-risk by relevant federal and state agencies. Updated information on its contemporary distribution and density is needed to guide and evaluate conservation and management actions. We combined 12 years (2009-2020) of detection and nondetection data collected throughout California and Oregon to model the potential distribution and density of SNRFs throughout their historical and contemporary ranges. We used an integrated species distribution and density modeling approach, which predicted SNRF density in sampled locations based on observed relationships between environmental covariates and detection frequencies, and then projected those predictions to unsampled locations based on the estimated correlations with environmental covariates. This approach provided predictions that serve as density estimates in sampled regions and projections in unsampled areas. Our model predicted a density of 1.06 (95% credible interval = 0.8-1.36) foxes per 100 km2 distributed throughout 22,926 km2 in three distinct regions of California and Oregon-Sierra Nevada, Lassen Peak, and Oregon Cascades. SNRFs were most likely to be found in areas with low minimum temperatures and high snow water equivalent. Our results provide a contemporary baseline to inform the development and evaluation of conservation and management actions, and guide future survey efforts.

A novel noninvasive genetic survey technique for small mammals

Noninvasive genetic surveys, often conducted by collecting fecal samples, have become a popular tool for surveying wildlife, but have primarily been applied to species with large and conspicuous scat. Although many small mammals are threatened, endangered, or data deficient, noninvasive genetic surveys have rarely been applied due to the challenges of detecting their inconspicuous fecal pellets. As part of a broader study of the endangered salt marsh harvest mouse (Reithrodontomys raviventris), we developed a noninvasive genetic survey technique for the community of small mammals in their putative range. We designed bait stations to passively collect fecal samples from rodents, and developed a multiplex primer set that amplified unique fragment sizes for salt marsh harvest mice and four other sympatric species. We tested the primer set on positive controls and on fecal pellets collected from bait stations at two regularly monitored field sites known to have very different densities of salt marsh harvest mice. The multiplex amplified DNA from all five species, even when all five species were present in a single sample. A positive species identification was made for all field-collected samples, and 43% of these field-collected samples had multispecies detections. The combination of bait stations and genetic species identification proved to be an effective means of noninvasively surveying small mammals in potential salt marsh harvest mouse habitat. The sampling technique should be applicable to a wide variety of small mammals in other systems.

Estimation of breeding population size using DNA-based pedigree reconstruction in brown bears

Robust estimates of demographic parameters are critical for effective wildlife conservation and management but are difficult to obtain for elusive species. We estimated the breeding and adult population sizes, as well as the minimum population size, in a high-density brown bear population on the Shiretoko Peninsula, in Hokkaido, Japan, using DNA-based pedigree reconstruction. A total of 1288 individuals, collected in and around the Shiretoko Peninsula between 1998 and 2020, were genotyped at 21 microsatellite loci. Among them, 499 individuals were identified by intensive genetic sampling conducted in two consecutive years (2019 and 2020) mainly by noninvasive methods (e.g., hair and fecal DNA). Among them, both parents were assigned for 330 bears, and either maternity or paternity was assigned to 47 and 76 individuals, respectively. The subsequent pedigree reconstruction indicated a range of breeding and adult (≥4 years old) population sizes: 128-173 for female breeders and 66-91 male breeders, and 155-200 for female adults and 84-109 male adults. The minimum population size was estimated to be 449 (252 females and 197 males) in 2019. Long-term continuous genetic sampling prior to a short-term intensive survey would enable parentage to be identified in a population with a high probability, thus enabling reliable estimates of breeding population size for elusive species.

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.

Genotyping-in-thousands by sequencing (GT-seq) of non-invasive fecal and degraded samples: a new panel to enable ongoing monitoring of Canadian polar bear populations

Genetic monitoring using noninvasive samples provides a complement or alternative to traditional population monitoring methods. However, next‐generation sequencing approaches to monitoring typically require high quality DNA and the use of noninvasive samples (e.g., scat) is often challenged by poor DNA quality and contamination by nontarget species. One promising solution is a highly multiplexed sequencing approach called genotyping‐in‐thousands by sequencing (GT‐seq), which can enable cost‐efficient genomics‐based monitoring for populations based on noninvasively collected samples. Here, we develop and validate a GT‐seq panel of 324 single nucleotide polymorphisms (SNPs) optimized for genotyping of polar bears based on DNA from noninvasively collected faecal samples. We demonstrate (1) successful GT‐seq genotyping of DNA from a range of sample sources, including successful genotyping (>50% loci) of 62.9% of noninvasively collected faecal samples determined to contain polar bear DNA; and (2) that we can reliably differentiate individuals, ascertain sex, assess relatedness, and resolve population structure of Canadian polar bear subpopulations based on a GT‐seq panel of 324 SNPs. Our GT‐seq data reveal spatial‐genetic patterns similar to previous polar bear studies but at lesser cost per sample and through use of noninvasively collected samples, indicating the potential of this approach for population monitoring. This GT‐seq panel provides the foundation for a noninvasive toolkit for polar bear monitoring and can contribute to community‐based programmes – a framework which may serve as a model for wildlife conservation and management for species worldwide.

Vital rates of two small populations of brown bears in Canada and range-wide relationship between population size and trend

Identifying mechanisms of population change is fundamental for conserving small and declining populations and determining effective management strategies. Few studies, however, have measured the demographic components of population change for small populations of mammals (<50 individuals). We estimated vital rates and trends in two adjacent but genetically distinct, threatened brown bear (Ursus arctos) populations in British Columbia, Canada, following the cessation of hunting. One population had approximately 45 resident bears but had some genetic and geographic connectivity to neighboring populations, while the other population had <25 individuals and was isolated. We estimated population-specific vital rates by monitoring survival and reproduction of telemetered female bears and their dependent offspring from 2005 to 2018. In the larger, connected population, independent female survival was 1.00 (95% CI: 0.96-1.00) and the survival of cubs in their first year was 0.85 (95% CI: 0.62-0.95). In the smaller, isolated population, independent female survival was 0.81 (95% CI: 0.64-0.93) and first-year cub survival was 0.33 (95% CI: 0.11-0.67). Reproductive rates did not differ between populations. The large differences in age-specific survival estimates resulted in a projected population increase in the larger population (λ = 1.09; 95% CI: 1.04-1.13) and population decrease in the smaller population (λ = 0.84; 95% CI: 0.72-0.95). Low female survival in the smaller population was the result of both continued human-caused mortality and an unusually high rate of natural mortality. Low cub survival may have been due to inbreeding and the loss of genetic diversity common in small populations, or to limited resources. In a systematic literature review, we compared our population trend estimates with those reported for other small populations (<300 individuals) of brown bears. Results suggest that once brown bear populations become small and isolated, populations rarely increase and, even with intensive management, recovery remains challenging.

Moving far, staying close: red fox dispersal patterns revealed by SNP genotyping

The genetic structure of a population can provide important insights into animal movements at varying geographical scales. Individual and social behaviors, such as philopatry and dispersal, affect patterns of relatedness, age and sex structure, shaping the local genetic structure of populations. However, these fine scale patterns may not be detected within broader population genetic structure. Using SNP genotyping for pairwise relatedness estimates, we investigated the spatial and genetic structuring of 141 red foxes within south-central Sweden at two scales. First, we looked at broad scale population structuring among red foxes at the regional level. We then estimated pairwise relatedness values to evaluate the spatial and genetic structure of male, female and mixed sex pairs for patterns of philopatry and dispersal at a more localized scale. We found limited genetic differentiation at the regional scale. However, local investigations revealed patterns of female philopatry and male biased dispersal. There were significant differences in pairwise geographic distances between highly related same sex pairs with the average distance between related males, 37.8 km, being six times farther than that of related females, averaging 6.3 km. In summary, the low levels of genetic differentiation found in this study illustrates the mobility and dispersal ability of red foxes across scales. However, relatedness plays a strong role in the spatial organization of red foxes locally, ultimately contributing to male biased dispersal patterns.

Multiscale landscape genetics of American marten at their southern range periphery

American marten (Martes americana) are a conservation priority in many forested regions of North America. Populations are fragmented at the southern edge of their distribution due to suboptimal habitat conditions. Facilitating gene flow may improve population resilience through genetic and demographic rescue. We used a multiscale approach to estimate the relationship between genetic connectivity and landscape characteristics among individuals at three scales in the northeastern United States: regional, subregional, and local. We integrated multiple modeling techniques and identified top models based on consensus. Top models were used to parameterize resistance surfaces at each scale, and circuit theory was used to identify potential movement corridors. Regional gene flow was affected by forest cover, elevation, developed land cover, and slope. At subregional and local scales, the effects were site specific and included subsets of temperature, elevation, developed land cover, and slope. Developed land cover significantly affected gene flow at each scale. At finer scales, lack of variance in forest cover may have limited the ability to detect a relationship with gene flow. The effect of slope on gene flow was positive or negative, depending on the site examined. Occupancy probability was a relatively poor predictor, and we caution its use as a proxy for landscape resistance. Our results underscore the importance of replication and multiscale approaches in landscape genetics. Climate warming and landscape conversion may reduce the genetic connectivity of marten populations in the northeastern United States, and represent the primary challenges to marten conservation at the southern periphery of their range.