Genetic footprints reveal geographic patterns of expansion in Fennoscandian red foxes

Population expansions of boreal species are among the most substantial ecological consequences of climate change, potentially transforming both structure and processes of northern ecosystems. Despite their importance, little is known about expansion dynamics of boreal species. Red foxes (Vulpes vulpes) are forecasted to become a keystone species in northern Europe, a process stemming from population expansions that began in the 19th century. To identify the relative roles of geographic and demographic factors and the sources of northern European red fox population expansion, we genotyped 21 microsatellite loci in modern and historical (1835–1941) Fennoscandian red foxes. Using Bayesian clustering and Bayesian inference of migration rates, we identified high connectivity and asymmetric migration rates across the region, consistent with source‐sink dynamics, whereby more recently colonized sampling regions received immigrants from multiple sources. There were no clear clines in allele frequency or genetic diversity as would be expected from a unidirectional range expansion from south to north. Instead, migration inferences, demographic models and comparison to historical red fox genotypes suggested that the population expansion of the red fox is a consequence of dispersal from multiple sources, as well as in situ demographic growth. Together, these findings provide a rare glimpse into the anatomy of a boreal range expansion and enable informed predictions about future changes in boreal communities.     

Ecology of wildlife rabies in Europe

• 1 The number of wildlife rabies cases has increased in Europe in recent years. We review the epizootiology of wildlife rabies in Europe, paying special attention to recent changes to the situation of two important vector species: the red fox and the raccoon dog. Red fox Vulpes vulpes has been the main vector of rabies since 1945, but the number and proportion of raccoon dog Nyctereutes procyonoides cases has rapidly increased during the past few years, particularly in north‐eastern Europe.
• 2 The transmission rate (average number of susceptible animals infected by each rabid animal) is critical for rabies spread and is partly determined by population density. Both raccoon dogs and foxes live in pairs. Foxes also live in family groups. Pairs and groups share their territories. Home range size usually correlates negatively with population density. Fox home ranges are 50–1500 ha, those of raccoon dogs 150–700 ha. The threshold value for rabies spread among foxes is estimated to be 0.63 individuals/km². Although fox density in eastern and northern Europe may be lower than this, the pooled density of foxes and raccoon dogs exceeds the threshold density.
• 3 Animal movements, especially dispersal of young, pose a risk for rabies spread. Although the likelihood of an epizootic is highest where fox and raccoon dog densities are highest, rabies may spread fastest where population densities are lower, because dispersal distances tend to correlate negatively with population density.
• 4 Oral vaccinations have been more effective in rabies control than culling foxes. Where two vector species exist, vaccination should be conducted twice a year, because most raccoon dogs disperse in autumn but some foxes do not disperse before mid‐ or late winter.
• 5 New rabies models, based on two vector species and their interaction, and which take into account the hibernation period of raccoon dogs, are needed for north‐eastern Europe.

     

Landscape genetics of the nonnative red fox of California

Invasive mammalian carnivores contribute disproportionately to declines in global biodiversity. In California, nonnative red foxes (Vulpes vulpes) have significantly impacted endangered ground‐nesting birds and native canids. These foxes derive primarily from captive‐reared animals associated with the fur‐farming industry. Over the past five decades, the cumulative area occupied by nonnative red fox increased to cover much of central and southern California. We used a landscape‐genetic approach involving mitochondrial DNA (mtDNA) sequences and 13 microsatellites of 402 nonnative red foxes removed in predator control programs to investigate source populations, contemporary connectivity, and metapopulation dynamics. Both markers indicated high population structuring consistent with origins from multiple introductions and low subsequent gene flow. Landscape‐genetic modeling indicated that population connectivity was especially low among coastal sampling sites surrounded by mountainous wildlands but somewhat higher through topographically flat, urban and agricultural landscapes. The genetic composition of populations tended to be stable for multiple generations, indicating a degree of demographic resilience to predator removal programs. However, in two sites where intensive predator control reduced fox abundance, we observed increases in immigration, suggesting potential for recolonization to counter eradication attempts. These findings, along with continued genetic monitoring, can help guide localized management of foxes by identifying points of introductions and routes of spread and evaluating the relative importance of reproduction and immigration in maintaining populations. More generally, the study illustrates the utility of a landscape‐genetic approach for understanding invasion dynamics and metapopulation structure of one of the world's most destructive invasive mammals, the red fox.     

Red fox takeover of arctic fox breeding den: an observation from Yamal Peninsula, Russia

Here, we report from the first direct observation of a red fox (Vulpes vulpes) intrusion on an arctic fox (Vulpes lagopus) breeding den from the southern Arctic tundra of Yamal Peninsula, Russia in 2007. At the same time, as a current range retraction of the original inhabitant of the circumpolar tundra zone the arctic fox is going on, the red fox is expanding their range from the south into arctic habitats. Thus, within large parts of the northern tundra areas the two species are sympatric which gives opportunities for direct interactions including interference competition. However, direct first-hand observations of such interactions are rare, especially in the Russian Arctic. In the present study, we observed one red fox taking over an arctic fox breeding den which resulted in den abandonment by the arctic fox. On July 19, eight arctic fox pups were observed on the den before the red fox was observed on the same den July 22. The pups were never seen at the den or elsewhere after the red fox was observed on the den for as long as we stayed in the area (until August 10). Our observation supports the view that direct interference with red fox on breeding dens may contribute to the range retraction of arctic foxes from the southern limits of the Arctic tundra in Russia.     

Natal den selection by sympatric arctic and red foxes on Herschel Island,Yukon, Canada

In the twentieth century, red fox (Vulpes vulpes) expanded into the Canadian Arctic, where it competes with arctic fox (Vulpes lagopus) for food and shelter. Red fox dominates in physical interactions with the smaller arctic fox, but little is known about competition between them on the tundra. On Hershel Island, north Yukon, where these foxes are sympatric, we focused on natal den choice, a critical aspect of habitat selection. We tested the hypothesis that red fox displaces arctic fox from dens in prey-rich habitats. We applied an approach based on model comparisons to analyse a 10-year data set and identify factors important to den selection. Red fox selected dens in habitats that were more prey-rich in spring. When red foxes reproduced, arctic fox selected dens with good springtime access, notably many burrows unblocked by ice and snow. These provided the best refuge early in the reproductive season. In the absence of red foxes, arctic foxes selected dens offering good shelter (i.e. large isolated dens). Proximity to prey-rich habitats was consistently less important than the physical aspects of dens for arctic fox. Our study shows for the first time that red foxes in the tundra select dens associated primarily with prey-rich areas, while sympatric arctic foxes do not. These results fit a model of red fox competitively interfering with arctic fox, the first detailed study of such competition in a true arctic setting.     

Oral vaccination of captive arctic foxes with lyophilized SAG2 rabies vaccine

Arctic foxes (Alopex lagopus) were immunized with lyophilized SAG2 oral rabies vaccine. The effectiveness of this vaccine was determined by serologic response and survival to challenge by rabies virus isolated from a red fox from Alaska (USA). No vaccine virus was found in saliva 1-72 hr after ingestion. At 2 wk after vaccination, all foxes had seroconverted, with rabies virus neutralizing antibody levels of 0.2-3.1 IU ml(-1). All vaccinated foxes survived to week 17 after challenge, and hippocampus, pons, and cerebellum were free of rabies virus as determined by direct immunofluorescence testing after death. One of four nonvaccinated foxes survived challenge and was free of rabies virus in neural tissue, and no rabies virus neutralizing antibody was detected in blood. Our results suggest that the lyophilized SAG2 oral rabies vaccine could be effective in arctic and subarctic regions, where freezing air and ground temperatures probably would not reduce its immunogenicity.     

Development of a genotype‐by‐sequencing immunogenetic assay as exemplified by screening for variation in red fox with and without endemic rabies exposure

Pathogens are recognized as major drivers of local adaptation in wildlife systems. By determining which gene variants are favored in local interactions among populations with and without disease, spatially explicit adaptive responses to pathogens can be elucidated. Much of our current understanding of host responses to disease comes from a small number of genes associated with an immune response. High‐throughput sequencing (HTS) technologies, such as genotype‐by‐sequencing (GBS), facilitate expanded explorations of genomic variation among populations. Hybridization‐based GBS techniques can be leveraged in systems not well characterized for specific variants associated with disease outcome to “capture” specific genes and regulatory regions known to influence expression and disease outcome. We developed a multiplexed, sequence capture assay for red foxes to simultaneously assess ~300‐kbp of genomic sequence from 116 adaptive, intrinsic, and innate immunity genes of predicted adaptive significance and their putative upstream regulatory regions along with 23 neutral microsatellite regions to control for demographic effects. The assay was applied to 45 fox DNA samples from Alaska, where three arctic rabies strains are geographically restricted and endemic to coastal tundra regions, yet absent from the boreal interior. The assay provided 61.5% on‐target enrichment with relatively even sequence coverage across all targeted loci and samples (mean = 50×), which allowed us to elucidate genetic variation across introns, exons, and potential regulatory regions (4,819 SNPs). Challenges remained in accurately describing microsatellite variation using this technique; however, longer‐read HTS technologies should overcome these issues. We used these data to conduct preliminary analyses and detected genetic structure in a subset of red fox immune‐related genes between regions with and without endemic arctic rabies. This assay provides a template to assess immunogenetic variation in wildlife disease systems.     

Molecular epidemiology of terrestrial rabies in the former Soviet Union

Fifty-five rabies virus isolates originating from different regions of the former Soviet Union (FSU) were compared with isolates originating from Eurasia, Africa, and North America according to complete or partial nucleoprotein (N) gene sequences. The FSU isolates formed five distinct groups. Group A represented viruses originating from the Arctic, which were similar to viruses from Alaska and Canada. Group B consisted of "Arctic-like" viruses, originating from the south of East Siberia and the Far East. Group C consisted of viruses circulating in the steppe and forest-steppe territories from the European part of Russia to Tuva and in Kazakhstan. These three phylogenetic groups were clearly different from the European cluster. Viruses of group D circulate near the western border of Russia. Their phylogenetic position is intermediate between group C and the European cluster. Group E consisted of viruses originating from the northwestern part of Russia and comprised a "northeastern Europe" group described earlier from the Baltic region. According to surveillance data, a specific host can be defined clearly only for group A (arctic fox; Alopex lagopus) and for the Far Eastern part of the group B distribution area (raccoon dog; Nyctereutes procyonoides). For other territories and rabies virus variants, the red fox (Vulpes vulpes) is the main virus reservoir. However, the steppe fox (Vulpes corsac), wolf (Canis lupus), and raccoon dog are also involved in virus circulation, depending on host population density. These molecular data, joined with surveillance information, demonstrate that the current fox rabies epizootic in the territory of the FSU developed independently of central and western Europe. No evidence of positive selection was found in the N genes of the isolates. In the glycoprotein gene, evidence of positive selection was strongly suggested in codons 156, 160, and 183. At these sites, no link between amino acid substitutions and phylogenetic placement or specific host species was detected.     

Elimination of rabies from red foxes in eastern Ontario

The province of Ontario (Canada) reported more laboratory confirmed rabid animals than any other state or province in Canada or the USA from 1958-91, with the exception of 1960-62. More than 95% of those cases occurred in the southern 10% of Ontario (approximately 100,000 km2), the region with the highest human population density and greatest agricultural activity. Rabies posed an expensive threat to human health and significant costs to the agricultural economy. The rabies variant originated in arctic foxes: the main vector in southern Ontario was the red fox (Vulpes vulpes), with lesser involvement of the striped skunk (Mephitis mephitis). The Ontario Ministry of Natural Resources began a 5 yr experiment in 1989 to eliminate terrestrial rabies from a approximately 30,000 km2 study area in the eastern end of southern Ontario. Baits containing oral rabies vaccine were dropped annually in the study area at a density of 20 baits/km2 from 1989-95. That continued 2 yr beyond the original 5 yr plan. The experiment was successful in eliminating the arctic fox variant of rabies from the whole area. In the 1980's, an average of 235 rabid foxes per year were reported in the study area. None have been reported since 1993. Cases of fox rabies in other species also disappeared. In 1995, the last bovine and companion animal cases were reported and in 1996 the last rabid skunk occurred. Only bat variants of rabies were present until 1999, when the raccoon variant entered from New York (USA). The success of this experiment led to an expansion of the program to all of southern Ontario in 1994. Persistence of terrestrial rabies, and ease of elimination, appeared to vary geographically, and probably over time. Ecological factors which enhance or reduce the long term survival of rabies in wild foxes are poorly understood.     

Arctic fox versus red fox in the warming Arctic: four decades of den surveys in north Yukon

During the last century, the red fox (Vulpes vulpes) has expanded its distribution into the Arctic, where it competes with the arctic fox (Vulpes lagopus), an ecologically similar tundra predator. The red fox expansion correlates with climate warming, and the ultimate determinant of the outcome of the competition between the two species is hypothesized to be climate. We conducted aerial and ground fox den surveys in the northern Yukon (Herschel Island and the coastal mainland) to investigate the relative abundance of red and arctic foxes over the last four decades. This region has undergone the most intense warming observed in North America, and we hypothesized that this climate change led to increasing dominance of red fox over arctic fox. Results of recent surveys fall within the range of previous ones, indicating little change in the relative abundance of the two species. North Yukon fox dens are mostly occupied by arctic fox, with active red fox dens occurring sympatrically. While vegetation changes have been reported, there is no indication that secondary productivity and food abundance for foxes have increased. Our study shows that in the western Arctic of North America, where climate warming was intense, the competitive balance between red and arctic foxes changed little in 40?years. Our results challenge the hypotheses linking climate to red fox expansion, and we discuss how climate warming’s negative effects on predators may be overriding positive effects of milder temperatures and longer growing seasons.     

Range‐wide multilocus phylogeography of the red fox reveals ancient continental divergence,minimal genomic exchange and distinct demographic histories

Widely distributed taxa provide an opportunity to compare biogeographic responses to climatic fluctuations on multiple continents and to investigate speciation. We conducted the most geographically and genomically comprehensive study to date of the red fox (Vulpes vulpes), the world's most widely distributed wild terrestrial carnivore. Analyses of 697 bp of mitochondrial sequence in ~1000 individuals suggested an ancient Middle Eastern origin for all extant red foxes and a 400 kya (SD = 139 kya) origin of the primary North American (Nearctic) clade. Demographic analyses indicated a major expansion in Eurasia during the last glaciation (~50 kya), coinciding with a previously described secondary transfer of a single matriline (Holarctic) to North America. In contrast, North American matrilines (including the transferred portion of Holarctic clade) exhibited no signatures of expansion until the end of the Pleistocene (~12 kya). Analyses of 11 autosomal loci from a subset of foxes supported the colonization time frame suggested by mtDNA (and the fossil record) but, in contrast, reflected no detectable secondary transfer, resulting in the most fundamental genomic division of red foxes at the Bering Strait. Endemic continental Y‐chromosome clades further supported this pattern. Thus, intercontinental genomic exchange was overall very limited, consistent with long‐term reproductive isolation since the initial colonization of North America. Based on continental divergence times in other carnivoran species pairs, our findings support a model of peripatric speciation and are consistent with the previous classification of the North American red fox as a distinct species, V. fulva.     

Rabies and canine distemper in an arctic fox population in Alaska

Two oil field workers were attacked by a rabid arctic fox (Alopex lagopus) in the Prudhoe Bay oil field (Alaska, USA) prompting officials to reduce the local fox population. Ninety-nine foxes were killed during winter 1994. We tested foxes for prevalence of rabies and canine distemper. Exposure to rabies was detected in five of 99 foxes. Of the five, only one fox had rabies virus in neural tissue as determined by the direct fluorescent antibody test. The other four foxes had been exposed to rabies, but had apparently produced antibodies and did not have an active infection. No evidence of canine distemper was detected as determined by the absence of distemper antibodies in serum and distemper virus in neural tissue.     

How ecological neighbourhoods influence the structure of the scavenger guild in low arctic tundra

Aim How the ecological neighbourhoods of coast and forest affect arctic tundra ecosystems is a pressing question as the circumpolar tundra belt is shrinking under global warming. Mobile facultative scavengers are likely to negatively impact tundra biodiversity as dominant competitors or predators, if they spill over into tundra. Here, we provide the first quantitative assessments of the structure of a scavenger guild in low arctic tundra with emphasis on how it changes along spatial gradients from neighbouring ecosystems (i.e. forest and coast) and with altitude (i.e. productivity gradients). We also assess the likelihood of interactions between guild members that may negatively impact vulnerable tundra species. Location North‐eastern part of Norway. Methods Extensive records of scavenger prevalence were obtained by deploying automatic digital cameras at experimental carcasses in tundra regions covering several thousand square kilometres and three winters in northern Norway. Main conclusions We found short‐range neighbourhood effects of forest and coast within the tundra scavenger guild. Species richness declined steeply with decreasing distance from the neighbouring ecosystems, in particular subarctic forest, and with increasing altitude. Bird species with strongholds in forest (golden eagle Aquila chrysaetos and hooded crow Corvus cornix) or along the coast (white‐tailed eagle Haliaeetus albicilla) were mostly responsible for short‐range neighbourhood effects on guild structure. However, the two most abundant guild members, the common raven Corvus corax and the red fox Vulpes vulpes, exhibited no spatial patterns within the range of neighbourhoods and altitudes examined. There was a clear diurnal segregation in the use of carcasses between birds and mammals reducing the likelihood of direct interactions between these two taxa. Presence of red fox appeared to exclude the arctic fox Vulpes lagopus, the only endemic tundra species within the guild, from carcasses.     

Fox colors in relation to colors in mice and sheep

Color inheritance in foxes is explained in terms of homology between color loci in foxes, mice, and sheep. The hypothesis presented suggests that the loci A (agouti), B (black/chocolate brown pigment) and E (extension of eumelanin vs. phaeomelanin) all occur in foxes, both the red fox, Vulpes vulpes, and the arctic fox, Alopex lagopus. Two alleles are postulated at each locus in each species. At the A locus, the (top) dominant allele in the red fox, Ar, produces red color and the corresponding allele in the arctic fox, Aw, produces the winter-white color. The bottom recessive allele in both species is a, which results in the black color of the silver fox and a rare black color in the Icelandic arctic fox when homozygous. The B alleles are assumed to be similar in both species: B, dominant, producing black eumelanin, and b, recessive, producing chocolate brown eumelanin when homozygous. The recessive E allele at the E locus in homozygous form has no effect on the phenotype determined by alleles at the A locus, while Ed, the dominant allele is epistatic to the A alleles and results in Alaska black in the red fox and the dark phase in the arctic fox. Genetic formulae of various color forms of red and arctic fox and their hybrids are presented.     

A restricted hybrid zone between native and introduced red fox (Vulpes vulpes) populations suggests reproductive barriers and competitive exclusion

Introduced species can threaten native taxa in multiple ways, including competition and hybridization, which can reduce fitness, alter ecological niches or swamp native genomes. Encroachment and hybridization by introduced species also provide opportunities to study the dynamics of invasiveness and hybridization during early stages following contact. We used 33 microsatellites, 51 single nucleotide polymorphisms and a mtDNA marker to characterize the extent and spatial pattern of encroachment and hybridization between a native, endemic subspecies of red fox (Vulpes vulpes patwin) and an introduced red fox population composed of highly admixed, phylogenetically divergent stock, resulting from a century of domestication. Both nuclear and mtDNA markers indicated that hybridization was primarily restricted to a narrow zone where the two populations came into contact. Although a few introgressed genotypes were detected in the interior of the native range, we found no immigrant foxes or F(1) or F(2) hybrids there, suggesting native foxes excluded introduced individuals. We speculate that the observed interbreeding at the periphery was facilitated by low densities. In total, 98% of mtDNA haplotypes in the native range were native and 96% of the nuclear ancestry was estimated to be native. Although the introduced range had expanded fivefold over the past four decades, native and non-native haplotypes from museum samples collected in and near the native range three decades earlier showed a similar geographic distribution as today, suggesting that the native range and hybrid zone were relatively stable. We hypothesize that the monogamous mating system of red foxes and other wild canids may enhance their resistance to hybridization because of greater fitness consequences associated with mate discrimination.     

The fluctuating world of a tundra predator guild: bottom‐up constraints overrule top‐down species interactions in winter

Global warming is predicted to change ecosystem functioning and structure in Arctic ecosystems by strengthening top‐down species interactions, i.e. predation pressure on small herbivores and interference between predators. Yet, previous research is biased towards the summer season. Due to greater abiotic constraints, Arctic ecosystem characteristics might be more pronounced in winter. Here we test the hypothesis that top‐down species interactions prevail over bottom‐up effects in Scandinavian mountain tundra (Northern Sweden) where effects of climate warming have been observed and top‐down interactions are expected to strengthen. But we test this ‘a priori’ hypothesis in winter and throughout the 3–4 yr rodent cycle, which imposes additional pulsed resource constraints. We used snowtracking data recorded in 12 winters (2004–2015) to analyse the spatial patterns of a tundra predator guild (arctic fox Vulpes lagopus, red fox Vulpes vulpes, wolverine Gulo gulo) and small prey (ptarmigan, Lagopus spp). The a priori top‐down hypothesis was then tested through structural equation modelling, for each phase of the rodent cycle. There was weak support for this hypothesis, with top‐down effects only discerned on arctic fox (weakly, by wolverine) and ptarmigan (by arctic fox) at intermediate and high rodent availability respectively. Overall, bottom‐up constraints appeared more influential on the winter community structure. Cold specialist predators (arctic fox and wolverine) showed variable landscape associations, while the boreal predator (red fox) appeared strongly dependent on productive habitats and ptarmigan abundance. Thus, we suggest that the unpredictability of food resources determines the winter ecology of the cold specialist predators, while the boreal predator relies on resource‐rich habitats. The constraints imposed by winters and temporary resource lows should therefore counteract productivity‐driven ecosystem change and have a stabilising effect on community structure. Hence, the interplay between summer and winter conditions should determine the rate of Arctic ecosystem change in the context of global warming.     

Survival and Cause-Specific Mortality of Corsac and Red Foxes in Mongolia

ABSTRACT The corsac fox (Vulpes corsac) and red fox (Vulpes vulpes) range widely across northern and central Asia and may be declining in many regions due to overhunting and other causes. However, details of the fundamental causes of survival and mortality of both species remain largely unquantified, but may be crucial for understanding interspecific relationships and developing effective conservation actions. We studied a radiomarked population of sympatric corsac and red foxes in central Mongolia to quantify survival and cause-specific mortality rates from April 2005 to April 2007. Survival probability was 0.34 for corsacs (n = 18) and 0.46 for red foxes (n = 17) and did not vary by year within or between each species. Among both foxes, mortality occurred mainly from hunting by humans, but also from predation by larger canids and unknown causes. Our results suggest that illegal human hunting represents the principal source of mortality for both species and that a recently initiated ranger patrol program in the study area did not affect fox survival. As such, more stringent protective measures will likely be necessary to halt declines of both foxes. Our results also suggest that interference competition occurs between species as red foxes killed but did not consume corsacs. Our results will be useful for developing science-based management strategies to protect foxes in Mongolia, and in understanding the competitive relationships between them.     

Prevalence of Toxoplasma gondii infection diagnosed by PCR in farmed red foxes, arctic foxes and raccoon dogs

The aim of this study was to compare Toxoplasma gondii infection in three canid species: red fox Vulpes vulpes, arctic fox Vulpes lagopus and raccoon dog Nyctereutesprocyonoides kept at the same farm. Anal swabs were taken from 24 adult and 10 juvenile red foxes, 12 adult arctic foxes, three adult and seven juvenile raccoon dogs. Additionally, muscle samples were taken from 10 juvenile red foxes. PCR was used to detect T. gondii DNA. T. gondii infection was not detected in any of the arctic foxes; 60% ofraccoon dogs were infected; the prevalence of the parasite in material from red fox swabs was intermediate between the prevalence observed in arctic foxes and raccoon dogs. It is possible that susceptibility and immune response to the parasite differ between the three investigated canid species. T. gondii DNA was detected in muscle tissue of five young foxes. The results of this study suggest that T. gondii infection is not rare in farmed canids.     

Responses of sympatric canids to human development revealed through citizen science

Measuring wildlife responses to anthropogenic activities often requires long‐term, large‐scale datasets that are difficult to collect. This is particularly true for rare or cryptic species, which includes many mammalian carnivores. Citizen science, in which members of the public participate in scientific work, can facilitate collection of large datasets while increasing public awareness of wildlife research and conservation. Hunters provide unique benefits for citizen science given their knowledge and interest in outdoor activities. We examined how anthropogenic changes to land cover impacted relative abundance of two sympatric canids, coyote (Canis latrans), and red fox (Vulpes vulpes) at a large spatial scale. In order to assess how land cover affected canids at this scale, we used citizen science data from bow hunter sighting logs collected throughout New York State, USA, during 2004–2017. We found that the two species had contrasting responses to development, with red foxes positively correlated and coyotes negatively correlated with the percentage of low‐density development. Red foxes also responded positively to agriculture, but less so when agricultural habitat was fragmented. Agriculture provides food and denning resources for red foxes, whereas coyotes may select forested areas for denning. Though coyotes and red foxes compete in areas of sympatry, we did not find a relationship between species abundance, likely a consequence of the coarse spatial resolution used. Red foxes may be able to coexist with coyotes by altering their diets and habitat use, or by maintaining territories in small areas between coyote territories. Our study shows the value of citizen science, and particularly hunters, in collection of long‐term data across large areas (i.e., the entire state of New York) that otherwise would unlikely be obtained.