Phylogeography of the North American red fox: vicariance in Pleistocene forest refugia

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.

Phylogeographic History of Endangered Hokuriku Salamander, Hynobius takedai (Amphibia: Caudata)

Knowledge of the phylogeographic history of organisms is valuable for understanding their evolutionary processes. To the best of our knowledge, the phylogeographic structure of Hokuriku salamander, Hynobius takedai, an endangered species, remains unclear. This study aimed to elucidate the phylogeographic history of H. takedai, which is expected to be strongly influenced by paleogeographic events. Phylogenetic analysis based on partial sequences of the mitochondrial DNA cytochrome b gene confirmed the genetic independence of H. takedai, and the divergence time with closely related species was estimated to be from the Late Pliocene to the Early Pleistocene. In the phylogenetic tree, two clades were identified within H. takedai, and their haplotypes were found in samples collected from the west and east of the distribution range. These intraspecific divergences were strongly influenced by geohistorical subdivisions of the current major distribution areas in the Middle Pleistocene. One clade was further subdivided and its formation may have been influenced by sea level changes in the Late Pleistocene.

Detection of Asian genetic components in autochthonous human Echinococcus multilocularis infections from endemic Warmia-Masuria (north-eastern Poland)

Alveolar echinococcosis is a life-threatening zoonotic disease caused by the larval stage of the cestode Echinococcus multilocularis. People are aberrant intermediate hosts accidentally infected with the parasite eggs via faecal-oral route, usually by the consumption of unwashed fruit and vegetable or direct contact with definitive hosts. The recently reported presence of Asian admixture in E. multilocularis tapeworms from Polish red foxes prompted the question of metacestode descent in the human population. In this study, a Maximum Likelihood tree based on partial sequences of E. multilocularis mitochondrial genes cox1, cob, and nad2 coupled with a hierarchical clustering analysis of microsatellite EmsB profiles and supplemented by Sammon's nonlinear mapping with k-means clustering revealed Asian genetic components, to date associated only with the sylvatic cycle, in two autochthonous samples from alveolar echinococcosis patients living in endemic Warmia-Masuria, north-eastern Poland. The red fox is the most likely source of contamination in the environment shared by people and wildlife that led to these infections. Our results confirm that Asian genetic variants participate in the synanthropic cycle in north-eastern Poland and indicate that they may be present in the human population in other areas where Asian genetic variants were detected in red foxes.

Dental and temporomandibular joint pathology of the red fox (Vulpes vulpes)

Skulls from 305 red foxes (Vulpes vulpes) were examined using predefined criteria to assess for the presence and severity of dental and temporomandibular joint (TMJ) pathology. Of the 305 specimens examined, 131 were male (42.9%), 103 were female (33.4%), 71 were of unknown sex (23.3%), 238 were adults (78.03%) and 67 were young adults (21.97%), with juveniles excluded from the study. Of a maximum of 12,810 possible teeth, 12,355 (96.5%) were present for examination, 72 (0.6%) were absent congenitally, 97 (0.8%) were acquired tooth loss and 280 (2.2%) were absent artefactually. Eight teeth (0.06%) in five specimens (1.6%) had abnormal tooth form and four (0.03%) supernumerary teeth were present across four (1.3%) specimens. Ten persistent deciduous teeth were present in seven (2.3%) specimens and six (0.05%) unerupted teeth were found in four (1.1%) specimens. Root number variation was present in 51 (0.4%) teeth, predominantly premolar teeth. Twenty-one (6.9%) specimens had lesions consistent with enamel hypoplasia and 42 (13.8%) showed root fenestrations in the maxillary alveolar bone. Periodontitis and attrition/abrasion were present in most specimens (56.7% and 96.1%, respectively). Three-hundred and eighty (3.1%) teeth were fractured, with the most common fracture type being root fractures. Eleven periapical lesions were found in nine (0.03%) specimens. Six-hundred and eight TMJs were evaluated and many specimens (11.5%) had evidence of low-grade TMJ osteoarthritis. Overall, these red foxes share similar dental pathology to other foxes, but had a higher prevalence of congenital tooth absence, persistent deciduous teeth, TMJ pathology and enamel hypoplasia.

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.

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.

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.

Phylogenomic systematics of the spotted skunks (Carnivora, Mephitidae, Spilogale): Additional species diversity and Pleistocene climate change as a major driver of diversification

Four species of spotted skunks (Carnivora, Mephitidae, Spilogale) are currently recognized: Spilogale angustifrons, S. gracilis, S. putorius, and S. pygmaea. Understanding species boundaries within this group is critical for effective conservation given that regional populations or subspecies (e.g., S. p. interrupta) have experienced significant population declines. Further, there may be currently unrecognized diversity within this genus as some taxa (e.g., S. angustifrons) and geographic regions (e.g., Central America) never have been assessed using DNA sequence data. We analyzed species limits and diversification patterns in spotted skunks using multilocus nuclear (ultraconserved elements) and mitochondrial (whole mitogenomes and single gene analysis) data sets from broad geographic sampling representing all currently recognized species and subspecies. We found a high degree of genetic divergence among Spilogale that reflects seven distinct species and eight unique mitochondrial lineages. Initial divergence between S. pygmaea and all other Spilogale occurred in the Early Pliocene (∼ 5.0 million years ago). Subsequent diversification of the remaining Spilogale into an "eastern" and a "western" lineage occurred during the Early Pleistocene (∼1.5 million years ago). These two lineages experienced temporally coincident patterns of diversification at ∼0.66 and ∼0.35 million years ago into two and ultimately three distinct evolutionary units, respectively. Diversification was confined almost entirely within the Pleistocene during a timeframe characterized by alternating glacial-interglacial cycles, with the origin of this diversity occurring in northeastern Mexico and the southwestern United States of America. Mitochondrial-nuclear discordance was recovered across three lineages in geographic regions consistent with secondary contact, including a distinct mitochondrial lineage confined to the Sonoran Desert. Our results have direct consequences for conservation of threatened populations, or species, as well as for our understanding of the evolution of delayed implantation in this enigmatic group of small carnivores.

Mitochondrial Genomes of the United States Distribution of Gray Fox (Urocyon cinereoargenteus) Reveal a Major Phylogeographic Break at the Great Plains Suture Zone

We examined phylogeographic structure in gray fox (Urocyon cinereoargenteus) across the United States to identify the location of secondary contact zone(s) between eastern and western lineages and investigate the possibility of additional cryptic intraspecific divergences. We generated and analyzed complete mitochondrial genome sequence data from 75 samples and partial control region mitochondrial DNA sequences from 378 samples to investigate levels of genetic diversity and structure through population- and individual-based analyses including estimates of divergence (FST and SAMOVA), median joining networks, and phylogenies. We used complete mitochondrial genomes to infer phylogenetic relationships and date divergence times of major lineages of Urocyon in the United States. Despite broad-scale sampling, we did not recover additional major lineages of Urocyon within the United States, but identified a deep east-west split (∼0.8 million years) with secondary contact at the Great Plains Suture Zone and confirmed the Channel Island fox (Urocyon littoralis) is nested within U. cinereoargenteus. Genetic diversity declined at northern latitudes in the eastern United States, a pattern concordant with post-glacial recolonization and range expansion. Beyond the east-west divergence, morphologically-based subspecies did not form monophyletic groups, though unique haplotypes were often geographically limited. Gray foxes in the United States displayed a deep, cryptic divergence suggesting taxonomic revision is needed. Secondary contact at a common phylogeographic break, the Great Plains Suture Zone, where environmental variables show a sharp cline, suggests ongoing evolutionary processes may reinforce this divergence. Follow-up study with nuclear markers should investigate whether hybridization is occurring along the suture zone and characterize contemporary population structure to help identify conservation units. Comparative work on other wide-ranging carnivores in the region should test whether similar evolutionary patterns and processes are occurring.

Insight into the Genetic Population Structure of Wild Red Foxes in Poland Reveals Low Risk of Genetic Introgression from Escaped Farm Red Foxes

In this study we assessed the level of genetic introgression between red foxes bred on fur farms in Poland and the native wild population. We also evaluated the impact of a geographic barrier and isolation by distance on gene flow between two isolated subpopulations of the native red fox and their genetic differentiation. Nuclear and mitochondrial DNA was collected from a total of 308 individuals (200 farm and 108 wild red foxes) to study non-native allele flow from farm into wild red fox populations. Genetic structure analyses performed using 24 autosomal microsatellites showed two genetic clusters as being the most probable number of distinct populations. No strong admixture signals between farm and wild red foxes were detected, and significant genetic differentiation was identified between the two groups. This was also apparent from the mtDNA analysis. None of the concatenated haplotypes detected in farm foxes was found in wild animals. The consequence of this was that the haplotype network displayed two genetically distinct groups: farm foxes were completely separated from native ones. Neither the River Vistula nor isolation by distance had a significant impact on gene flow between the separated wild red fox subpopulations. The results of our research indicate a low probability of genetic introgression between farm and native red foxes, and no threat to the genetic integrity of this species.

The role of a mechanistic host in maintaining arctic rabies variant distributions: Assessment of functional genetic diversity in Alaskan red fox (Vulpes vulpes)

Populations are exposed to different types and strains of pathogens across heterogeneous landscapes, where local interactions between host and pathogen may present reciprocal selective forces leading to correlated patterns of spatial genetic structure. Understanding these coevolutionary patterns provides insight into mechanisms of disease spread and maintenance. Arctic rabies (AR) is a lethal disease with viral variants that occupy distinct geographic distributions across North America and Europe. Red fox (Vulpes vulpes) are a highly susceptible AR host, whose range overlaps both geographically distinct AR strains and regions where AR is absent. It is unclear if genetic structure exists among red fox populations relative to the presence/absence of AR or the spatial distribution of AR variants. Acquiring these data may enhance our understanding of the role of red fox in AR maintenance/spread and inform disease control strategies. Using a genotyping-by-sequencing assay targeting 116 genomic regions of immunogenetic relevance, we screened for sequence variation among red fox populations from Alaska and an outgroup from Ontario, including areas with different AR variants, and regions where the disease was absent. Presumed neutral SNP data from the assay found negligible levels of neutral genetic structure among Alaskan populations. The immunogenetically-associated data identified 30 outlier SNPs supporting weak to moderate genetic structure between regions with and without AR in Alaska. The outliers included SNPs with the potential to cause missense mutations within several toll-like receptor genes that have been associated with AR outcome. In contrast, there was a lack of genetic structure between regions with different AR variants. Combined, we interpret these data to suggest red fox populations respond differently to the presence of AR, but not AR variants. This research increases our understanding of AR dynamics in the Arctic, where host/disease patterns are undergoing flux in a rapidly changing Arctic landscape, including the continued northward expansion of red fox into regions previously predominated by the arctic fox (Vulpes lagopus).

Relationships between fox populations and rabies virus spread in northern Canada

Rabies spreads in both Arctic (Vulpes lagopus) and red foxes (Vulpes vulpes) throughout the Canadian Arctic but limited wildlife disease surveillance, due to the extensive landmass of the Canadian north and its small widely scattered human population, undermines our knowledge of disease transmission patterns. This study has explored genetic population structure in both the rabies virus and its fox hosts to better understand factors that impact rabies spread. Phylogenetic analysis of 278 samples of the Arctic lineage of rabies virus recovered over 40 years identified four sub-lineages, A1 to A4. The A1 lineage has been restricted to southern regions of the Canadian province of Ontario. The A2 lineage, which predominates in Siberia, has also spread to northern Alaska while the A4 lineage was recovered from southern Alaska only. The A3 sub-lineage, which was also found in northern Alaska, has been responsible for virtually all cases across northern Canada and Greenland, where it further differentiated into 18 groups which have systematically evolved from a common predecessor since 1975. In areas of Arctic and red fox sympatry, viral groups appear to circulate in both hosts, but both mitochondrial DNA control region sequences and 9-locus microsatellite genotypes revealed contrasting phylogeographic patterns for the two fox species. Among 157 Arctic foxes, 33 mitochondrial control region haplotypes were identified but little genetic structure differentiating localities was detected. Among 162 red foxes, 18 control region haplotypes delineated three groups which discriminated among the Churchill region of Manitoba, northern Quebec and Labrador populations, and the coastal Labrador locality of Cartwright. Microsatellite analyses demonstrated some genetic heterogeneity among sampling localities of Arctic foxes but no obvious pattern, while two or three clusters of red foxes suggested some admixture between the Churchill and Quebec-Labrador regions but uniqueness of the Cartwright group. The limited population structure of Arctic foxes is consistent with the rapid spread of rabies virus subtypes throughout the north, while red fox population substructure suggests that disease spread in this host moves most readily down certain independent corridors such as the northeastern coast of Canada and the central interior. Interestingly the evidence suggests that these red fox populations have limited capacity to maintain the virus over the long term, but they may contribute to viral persistence in areas of red and Arctic fox sympatry.

Sequencing Red Fox Y Chromosome Fragments to Develop Phylogenetically Informative SNP Markers and Glimpse Male-Specific Trans-Pacific Phylogeography

The red fox (Vulpes vulpes) has a wide global distribution with many ecotypes and has been bred in captivity for various traits, making it a useful evolutionary model system. The Y chromosome represents one of the most informative markers of phylogeography, yet it has not been well-studied in the red fox due to a lack of the necessary genomic resources. We used a target capture approach to sequence a portion of the red fox Y chromosome in a geographically diverse red fox sample, along with other canid species, to develop single nucleotide polymorphism (SNP) markers, 13 of which we validated for use in subsequent studies. Phylogenetic analyses of the Y chromosome sequences, including calibration to outgroups, confirmed previous estimates of the timing of two intercontinental exchanges of red foxes, the initial colonization of North America from Eurasia approximately half a million years ago and a subsequent continental exchange before the last Pleistocene glaciation (~100,000 years ago). However, in contrast to mtDNA, which showed unidirectional transfer from Eurasia to North America prior to the last glaciation, the Y chromosome appears to have been transferred from North America to Eurasia during this period. Additional sampling is needed to confirm this pattern and to further clarify red fox Y chromosome phylogeography.

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

Genetic factors in the decline of small populations are extremely difficult to study in nature. We leveraged a natural experiment to investigate evidence of inbreeding depression and genetic rescue in a remnant population of subalpine-specialized Sierra Nevada red foxes (Vulpes vulpes necator) using noninvasive genetic monitoring during 2010-2017. Only 7 individuals were detected in the first 2 years. These individuals assigned genetically to the historical population and exhibited genetic hallmarks of inbreeding and no evidence of reproduction. Two years into the study, we detected 2 first-generation immigrant males from a recently expanding population of red foxes in the Great Basin Desert. Through annual resampling of individuals (634 red fox DNA samples, 41 individuals) and molecular reconstruction of pedigrees, we documented 1-3 litters/year for 5 years, all descended directly or indirectly from matings involving immigrant foxes. The observed heterozygosity and allelic richness of the population nearly doubled in 2 years. Abundance increased, indicative of a rapidly expanding population. Throughout the study, adult survival was high. Restoration of gene flow apparently improved the demographic trajectory of this population in the short term. Whether these benefits continue in the longer term could depend on numerous factors, such as maintenance of any locally adapted alleles. This study highlights the value of noninvasive genetic monitoring to assess rapidly shifting conditions in small populations. Uncertainties about the longer-term trajectory of this population underscore the need to continue monitoring and to research potential for both negative and positive aspects of continued genetic infusion.

Phylogeography of moose in western North America

Subspecies designations within temperate species’ ranges often reflect populations that were isolated by past continental glaciation, and glacial vicariance is believed to be a primary mechanism behind the diversification of several subspecies of North American cervids. We used genetics and the fossil record to study the phylogeography of three moose subspecies (Alces alces andersoni, A. a. gigas, and A. a. shirasi) in western North America. We sequenced the complete mitochondrial genome (16,341 base pairs; n = 60 moose) and genotyped 13 nuclear microsatellites (n = 253) to evaluate genetic variation among moose samples. We also reviewed the fossil record for detections of all North American cervids to comparatively assess the evidence for the existence of a southern refugial population of moose corresponding to A. a. shirasi during the last glacial maximum of the Pleistocene. Analysis of mtDNA molecular variance did not support distinct clades of moose corresponding to currently recognized subspecies, and mitogenomic haplotype phylogenies did not consistently distinguish individuals according to subspecies groupings. Analysis of population structure using microsatellite loci showed support for two to five clusters of moose, including the consistent distinction of a southern group of moose within the range of A. a. shirasi. We hypothesize that these microsatellite results reflect recent, not deep, divergence and may be confounded by a significant effect of geographic distance on gene flow across the region. Review of the fossil record showed no evidence of moose south of the Wisconsin ice age glaciers ≥ 15,000 years ago. We encourage the integration of our results with complementary analyses of phenotype data, such as morphometrics, originally used to delineate moose subspecies, for further evaluation of subspecies designations for North American moose.

Nuclear Genetic Analysis of the Red Fox Across its Trans-Pacific Range

The red fox (Vulpes vulpes) occurs on multiple continents in diverse habitats, making it an informative system for evolutionary genomic research. However, its phylogeography remains unclear. Previously, mitochondrial DNA and small numbers of nuclear loci provided discordant views. Both markers indicated deep divergence (~ 0.5 million years [MY]) between Eurasian and southern North American populations but differed in the apparent continental affinity of Alaskan red foxes, implying some degree of gene exchange during secondary contact (~0.1 MY). We assayed >173000 nuclear genomic sites in 52 red foxes, along with 2 Rueppell's foxes (Vulpes rueppellii) and a gray wolf (Canis lupus) using the Illumina CanineHD BeadChip. We obtained 5107 single nucleotide polymorphisms (SNPs) in the foxes. Consistent with the Afro-Eurasian origins of red foxes, genetic diversity was higher in Eurasian than North American samples. Phylogenetic trees indicated that Alaskan and southern North American red foxes formed a monophyletic group nested within the Eurasian clade. However, admixture models suggested Alaskan red foxes contained up to 40% Eurasian ancestry. We hypothesize that North American red foxes either hybridized with Eurasian foxes in Beringia at the start of the last glaciation or merged with a Beringian population after the last glaciation. Future work is needed to test between these scenarios and assess speciation.

Admixture analyses and phylogeographic relationships reveal complete genetic distinctiveness of Polish farm and wild red foxes (Vulpes vulpes) and the North American origin of farm-bred individuals

A number of studies showed that many mtDNA haplotypes were shared among contemporary farm red foxes bred on different continents and the historical wild red foxes of North American origin. Therefore, in this study, the population genetic structure and phylogeographic relationships of Polish red foxes kept on fur farms and their wild conspecifics were investigated to assess the ancestry of the farm red foxes in Poland. A total of 330 tissue samples (200 from farm foxes and 130 from wild foxes) were used for the genetic analyses. Thirty microsatellite loci and two regions of mtDNA were used to assess the level of admixture between farm- and wild red foxes, to construct haplotype networks and create a phylogenetic tree. The genetic structure analysis clearly indicated two genetic clusters as being the most probable number of genetically distinct populations. The fixation index revealed a significant genetic distance between the farm- and wild red fox populations (F_ST  = 0.27, p < 0.05). Haplotype networks based on frequencies showing relationships between concatenated haplotypes of Polish farm- and wild red foxes and the constructed phylogenetic tree clearly indicated two genetically distinct groups. The results of this study provide strong evidence confirming the North American origin of red foxes bred on Polish farms and the genetic distinctiveness of both studied populations.

Assessment of the global pattern of genetic diversity in Echinococcus multilocularis inferred by mitochondrial DNA sequences

The aim of this review was to assess our current knowledge on phylogeography and global genetic structure of Echinococcus multilocularis populations originating from rodents, wild canid hosts, and human. Six bibliographic databases were searched from 1990 to 2017, identifying a total of 110 publications. The cytochrome c oxidase subunit 1 (cox1) and cytochrome b (cytb) sequences of E. multilocularis from Asia, Europe, and North Americas were analyzed to estimate the diversity and neutrality indices, and genetic differentiation. A total of 69 (cox1, 36.7%) and 16 haplotypes (cytb, 19.2%) were grouped into various geographical clades. A parsimonious haplotype network demonstrated a star-like feature with haplo-groups Em2 (Asia: 36%), Em105 (Eastern Tibetan plateau: 4.8%), Em46 (Europe: 9.1%), Em73, (Europe: 2.7%) and Em92 (North Americas: 4.3%) as the most common haplotypes. A relatively high level of genetic diversity was detected in rodent-derived E. multilocularis isolates (Haplotype diversity: 0.944), wild canids (Hd: 0.912), and human origin (Hd: 0.704). The highest number of haplotypes (n = 59) and the highest haplotype diversity (0.969) were identified in the Asian and European populations, respectively. Cladistic phylogenetic tree indicated the European clade has a sister relationship with the Asian clade. However, some North American haplotypes were assigned to the European clade together with haplotypes from Poland. The statistically significant Fst values indicated that E. multilocularis populations of Asian-European, Asian-North American, and European-North American origins were genetically differentiated (Fst: 0.22624 to 0.43059). An occurrence of distinct parasite populations suggests that E. multilocularis derived from glacial refugia have been plausibly sustained by indigenous hosts during the Pleistocene Epoch.

Red Fox Ancestry and Connectivity Assessments Reveal Minimal Fur Farm Introgression in Greater Yellowstone Ecosystem

Rocky Mountain red foxes Vulpes vulpes macroura potentially encounter other red fox Vulpes vulpes lineages at lower elevations, which may include nonindigenous red foxes derived from fur farms. Introgression from nonindigenous red foxes could have negative evolutionary consequences for the rare Rocky Mountain red fox subspecies. Red foxes at high elevations in the Greater Yellowstone Ecosystem exhibit lighter coat colors than those at lower elevations, potentially indicating that they represent the indigenous subspecies and that gene flow across the elevational gradient is restricted. We collected tissue samples across a 1,750-m elevation range and examined mitochondrial DNA sequences and nuclear DNA microsatellite genotypes to assess the ancestry and genetic population structure of red foxes in the northern Greater Yellowstone Ecosystem. We also used reference samples from fur farm red foxes and indigenous red foxes of the western United States to assess the extent of nonindigenous introgression across the ecosystem. We found little overlap in the elevational distribution of maternally inherited mitochondrial DNA haplotypes: above 1,600 m, we only found indigenous Rocky Mountain haplotypes (n = 4), whereas below 1,600 m, we found haplotypes not indigenous to the Rocky Mountains (n = 5) that were associated with fur farms or indigenous to the Great Plains. In contrast, biparentally inherited microsatellite variation showed little population structure across the elevational gradient. Despite this evidence of nuclear gene flow across the elevational gradient, we found little fur farm introgression in the microsatellite genotypes. It is possible that long-standing nuclear (but apparently not mitochondrial) gene flow between Rocky Mountain red foxes and indigenous red foxes on the Great Plains explained the low nuclear differentiation of these populations. Importantly, our results suggested that high elevations of the northern Greater Yellowstone Ecosystem remained free of significant fur farm introgression. Mitonuclear discordance could reflect sex-biased dispersal, which we hypothesize could be the effect of elevational differences in reproductive phenology.

Predicting the Potential Distribution of the Sierra Nevada Red Fox in the Oregon Cascades

The Sierra Nevada red fox Vulpes vulpes necator is a native subspecies associated with subalpine regions in the Sierra Nevada and Cascade mountain ranges of California and Oregon. In the past century, the Sierra Nevada red fox experienced a major range contraction and decline in California. However, the number, size, and connectivity of populations extant in Oregon remain unclear. This knowledge gap impedes efficient monitoring and hinders development of a cohesive conservation strategy for the subspecies. The historical range is large and includes rugged terrain with low accessibility; therefore, a predictive model is needed to facilitate more comprehensive and systematic surveys in the future. We initiated a multiagency collaborative effort to survey portions of the range in the Oregon Cascades during 2011–2016 (verified genetic and photographic detections) and to assemble existing sighting reports dating back to 1985 (unverified), which we used to create Maxent models to predict the potential distribution of Sierra Nevada red fox within Oregon. To identify optimal levels of model complexity, we compared cross-validation accuracy of models that varied in levels of protection against overfitting (regularization). The highest-performing models utilized intermediate regularization, and included minimum January temperature and land-cover type. Regardless of regularization or data set (verified detections, all putative detections), all models agreed in predictions of a high-probability region covering approximately 3,470 km2 or 6% of the Cascade region, corresponding to the high-elevation portion of the crest. With the exception of a gap between Mount Hood and Mt. Jefferson, this core area of predicted presence was continuous along the north–south extent of the crest, suggesting a capacity for high connectivity among observed clusters of occurrence. Use of modeled potential distributions in future survey design will improve efficiency of field data collection, facilitating more precise evaluations of the distribution, abundance, and genetic integrity and connectivity of Sierra Nevada red fox in Oregon.

Y-Chromosome Markers for the Red Fox

The de novo assembly of the red fox (Vulpes vulpes) genome has facilitated the development of genomic tools for the species. Efforts to identify the population history of red foxes in North America have previously been limited by a lack of information about the red fox Y-chromosome sequence. However, a megabase of red fox Y-chromosome sequence was recently identified over 2 scaffolds in the reference genome. Here, these scaffolds were scanned for repeated motifs, revealing 194 likely microsatellites. Twenty-three of these loci were selected for primer development and, after testing, produced a panel of 11 novel markers that were analyzed alongside 2 markers previously developed for the red fox from dog Y-chromosome sequence. The markers were genotyped in 76 male red foxes from 4 populations: 7 foxes from Newfoundland (eastern Canada), 12 from Maryland (eastern United States), and 9 from the island of Great Britain, as well as 48 foxes of known North American origin maintained on an experimental farm in Novosibirsk, Russia. The full marker panel revealed 22 haplotypes among these red foxes, whereas the 2 previously known markers alone would have identified only 10 haplotypes. The haplotypes from the 4 populations clustered primarily by continent, but unidirectional gene flow from Great Britain and farm populations may influence haplotype diversity in the Maryland population. The development of new markers has increased the resolution at which red fox Y-chromosome diversity can be analyzed and provides insight into the contribution of males to red fox population diversity and patterns of phylogeography.

St. Matthew Island colonized through multiple long-distance red fox (Vulpes vulpes) dispersal events

Expansion of red fox (Vulpes vulpes (L., 1758)) into new arctic habitat and the potential for competition with arctic fox (Vulpes lagopus (L., 1758)) are of considerable conservation concern. Previous work has focused on red fox expanding into contiguous areas with few barriers to dispersal. Here, we examine mitochondrial DNA in red fox on recently colonized St. Matthew Island in the Bering Sea to determine their ultimate origin. Though limited in sample size (n = 7), we found that St. Matthew Island was colonized by North American lineages; surprisingly, despite the >400 km distance to the mainland, we found the island was colonized by at least three mitochondrial matrilines. These results suggest that even extremely isolated places may be colonized by red fox, and that the over-ice or over-ocean dispersal ability of red fox may have been previously underappreciated.

Phylogeography of a widespread small carnivore, the western spotted skunk (Spilogale gracilis) reveals temporally variable signatures of isolation across western North America

We analyzed phylogeographic patterns in the western spotted skunk, Spilogale gracilis Merriam, 1890 (Carnivora: Mephitidae) in relation to historical events associated with Pre-Pleistocene Divergence (PPD) and Quaternary climate change (QCC) using mitochondrial DNA from 97 individuals distributed across Western North America. Divergence times were generated using BEAST to estimate when isolation in putative refugia occurred. Patterns and timing of demographic expansion was performed using Bayesian skyline plot. Putative climatic refugia resulting from Quaternary climate change were identified using paleoecological niche modeling and divergence dates compared to major vicariant events associated with Pre-Pleistocene conditions. We recovered three major mitochondrial clades corresponding to western North America (California, Baja, and across the Great Basin), east-central North America (Texas, central Mexico, New Mexico), and southwestern Arizona/northwestern Mexico. Time to most recent common ancestor for S. gracilis occurred ~1.36 Ma. Divergence times for each major clade occurred between 0.25 and 0.12 Ma, with signature of population expansion occurring 0.15 and 0.10 Ma. Ecological niche models identified three potential climatic refugia during the Last Interglacial, (1) west coast of California and Oregon, (2) northwestern Mexico, and (3) southern Texas/northeastern Mexico as well as two refugia during the Last Glacial Maximum, (1) western USA and (2) southern Texas/northeastern Mexico. This study supports PPD in shaping species-level diversity compared to QCC-driven changes at the intraspecific level for Spilogale, similar to the patterns reported for other small mammals (e.g., rodents and bats). Phylogeographic patterns also appear to have been shaped by both habitat and river vicariance, especially across the desert southwest. Further, continuing climate change during the Holocene coupled with anthropogenic modifications during the Anthropocene appears to be removing both of these barriers to current dispersal of western spotted skunks.

Variation in fur farm and wild populations of the red fox, Vulpes vulpes (Carnivora: Canidae) — Part I: Morphometry

This paper demonstrates the influence of artificial selection on morphometric traits in the red fox [Vulpes vulpes (Linnaeus, 1758)]. Measurements and two proportion coefficients were analysed in 132 wild and 199 farm red foxes. The two groups differed significantly (P ≤ 0.05) on all but one of the measurements. Eight out of 11 measurements were significantly greater in the farm fox population, while only tail length, ear height, and length of the right hind limb were greater in the population of wild foxes. The opposite trend was observed when analysing variation in the measurements — the farm foxes were characterized by a greater variability only in the case of body weight, body length, and breadth of chest. When analysing the sexual dimorphism index in different sex and population groups, in almost all analysed traits, the greatest differences occurred between farm males and wild females. All of the traits examined in this study are important for survival of wild foxes. However, because importance of some traits was reduced during domestication and selective breeding (farm foxes do not have to fight for survival), the genetic relationship between them may have weakened. Other possible causes of morphological differences between the studied groups of red foxes are discussed as well.

High gene flow in the American badger overrides habitat preferences and limits broadscale genetic structure

Habitat associations are a function of habitat preferences and dispersal capabilities, both of which can influence how species responded to Quaternary climatic changes and contemporary habitat heterogeneity. Predicting resultant genetic structure is not always straightforward, especially in species where high dispersal potential and habitat preferences yield opposing predictions. The American badger has high dispersal capabilities that predict widespread panmixia, but avoids closed-canopy forests and clay soils, which could restrict gene flow and create ecologically based population genetic structure. We used mitochondrial sequence and microsatellite data sets to characterize how these opposing forces contribute to genetic structure in badgers at a continent-wide scale. Our data revealed an overall lack of ecologically based population genetic structure, suggesting that high dispersal capabilities were sufficiently realized to overcome most habitat-based genetic structure. At a broadscale, badger gene flow is limited only by geographic distance (isolation by distance) and large water barriers (Lake Michigan and the Mississippi River). The absence of genetic structure in a species with strong avoidance of unsuitable habitats advances our understanding of when and how genetic structure emerges in widespread, highly mobile species.

Inconsistent phylogeographic pattern between a sperm dependent fish and its host: in situ hybridization vs dispersal

Background

Co-dispersal of sperm-dependent hybrids and their sexual relatives is expected to result in consistent spatial patterns between assemblages of hybrids and genetic structure of parental species. However, local hybridization events may blur this signal as assemblages could be organized under different connectivity constraints. This study aims at testing the hypothesis of local hybridization events by comparing the assemblage of hybrid fish Chrosomus eos-neogaeus to the genetic diversity of one of its parental species, Chrosomus eos.

Results

An extensive survey performed on a total of 132 sites located in two regions of Southern Quebec (West-Qc and East-Qc) revealed a distinct organization of hybrid lineages. One of the six hybrid lineages detected in West-Qc is widespread throughout this region resulting in a low α-diversity (1.38) and β-diversity (4.35). On the other hand, 36 hybrid lineages were detected in East-Qc and displayed narrow geographic distributions leading to a high α-diversity (2.30) and β-diversity (15.68). In addition, the C. eos multilocus haplotype of several of these hybrids is assigned to their respective sympatric C. eos population. Finally, contrasting with hybrids, the paternal species C. eos displayed a higher ρ_ST in West-Qc (0.2300) than in East-Qc (0.0734).

Conclusion

The unusually high diversity of hybrid lineages in East-Qc as well as the spatial organization and the close genetic relationship with C. eos sympatric populations support the hypothesis that multiple hybridization events occurred in situ. These findings coupled to the near absence of the maternal species Chrosomus neogeaus suggest that the decline of this species could be the trigger event at the origin of the high rates of spontaneous hybridization in this region.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0754-5) contains supplementary material, which is available to authorized users.

Spatio-temporal Analysis of the Genetic Diversity of Arctic Rabies Viruses and Their Reservoir Hosts in Greenland

There has been limited knowledge on spatio-temporal epidemiology of zoonotic arctic fox rabies among countries bordering the Arctic, in particular Greenland. Previous molecular epidemiological studies have suggested the occurrence of one particular arctic rabies virus (RABV) lineage (arctic-3), but have been limited by a low number of available samples preventing in-depth high resolution phylogenetic analysis of RABVs at that time. However, an improved knowledge of the evolution, at a molecular level, of the circulating RABVs and a better understanding of the historical perspective of the disease in Greenland is necessary for better direct control measures on the island. These issues have been addressed by investigating the spatio-temporal genetic diversity of arctic RABVs and their reservoir host, the arctic fox, in Greenland using both full and partial genome sequences. Using a unique set of 79 arctic RABV full genome sequences from Greenland, Canada, USA (Alaska) and Russia obtained between 1977 and 2014, a description of the historic context in relation to the genetic diversity of currently circulating RABV in Greenland and neighboring Canadian Northern territories has been provided. The phylogenetic analysis confirmed delineation into four major arctic RABV lineages (arctic 1–4) with viruses from Greenland exclusively grouping into the circumpolar arctic-3 lineage. High resolution analysis enabled distinction of seven geographically distinct subclades (3.I – 3.VII) with two subclades containing viruses from both Greenland and Canada. By combining analysis of full length RABV genome sequences and host derived sequences encoding mitochondrial proteins obtained simultaneously from brain tissues of 49 arctic foxes, the interaction of viruses and their hosts was explored in detail. Such an approach can serve as a blueprint for analysis of infectious disease dynamics and virus-host interdependencies. The results showed a fine-scale spatial population structure in Greenland arctic foxes based on mitochondrial sequences, but provided no evidence for independent isolated evolutionary development of RABV in different arctic fox lineages. These data are invaluable to support future initiatives for arctic fox rabies control and elimination in Greenland.

Next to dog-mediated rabies, wildlife rabies continues to pose a public health problem, particularly in the northern hemisphere. Control of this zoonosis at the animal source has been proven the most efficient route to reduction of human rabies burden. Successful elimination of red fox-mediated rabies in Western Europe and parts of North America has demonstrated the viability of wildlife rabies control strategies. In some regions, the epidemiology of wildlife rabies is well understood; this is not the case for arctic rabies, particularly in Greenland. Previous molecular epidemiological studies demonstrated the occurrence of one particular arctic rabies virus (RABV) lineage (arctic-3) but were limited by low sample numbers and limited sequence length so as to preclude generation of high resolution phylogenetic analysis. Here, a unique set comprised of 79 complete genome sequences of RABVs from Greenland, Canada, USA (Alaska) and Russia collected over the past four decades was analysed. The use of next generation sequencing (NGS) allowed simultaneous determination of host derived sequences encoding mitochondrial proteins from the same brain tissue of 49 arctic foxes. These sequence data combined with geographical and temporal information permit the study of the genetic diversity and evolution of circulating RABVs in Greenland against the background of reservoir host genetics. The results reveal the existence of a single arctic RABV lineage (arctic-3) in Greenland, which has evolved into multiple distinct variants. These analyses provide an improved knowledge of the evolution of the circulating viruses at the molecular level and a better understanding of the historical perspective of the disease in Greenland compared to other parts of the Arctic. This knowledge will support policy on rabies control in mammalian wildlife reservoirs.

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.

Endoparasitic fauna of red foxes (Vulpes vulpes) and golden jackals (Canis aureus) in Serbia

Wild canides have a high epizootiological - epidemiological significance, considering that they are hosts for some parasites which spread vector born diseases. Increased frequency of certain interactions between domestic and wild canides increases the risk of occurrence, spreading and maintaining the infection of parasitic etiology in domestic canides. The research was conducted in 232 wild canides (172 red foxes and 60 golden jackals). The examined material was sampled from foxes and jackals, which were hunted down between 2010 and 2014, from 8 epizootiological areas of Serbia (North-Bačka, West-Bačka, Southern-Banat, Moravički, Zlatiborski, Raški, Rasinski and Zaječarski district). On completing the parasitological dissection and the coprological diagnostics, in wild canides protozoa from the genus Isospora were identified, 3 species of trematoda (Alaria alata, Pseudamphistomum truncatum and Metagonimus yokogawai), cestods from the genus Taenia and 5 species of nematodes (Toxocara canis, Ancylostomatidae, Trichuris vulpis and Capillaria aerophila). The finding of M. yokogawai in golden jackals were, to the best of our knowledge, one of the first diagnosed cases of metagonimosis in golden jackals in Serbia. The continued monitoring of the parasitic fauna of wild canides is needed to establish the widespread of the zoonoses in different regions of Serbia, because they present the reservoirs and/or sources of these infections.

Population structure of two rabies hosts relative to the known distribution of rabies virus variants in Alaska

For pathogens that infect multiple species, the distinction between reservoir hosts and spillover hosts is often difficult. In Alaska, three variants of the arctic rabies virus exist with distinct spatial distributions. We tested the hypothesis that rabies virus variant distribution corresponds to the population structure of the primary rabies hosts in Alaska, arctic foxes (Vulpes lagopus) and red foxes (Vulpes vulpes) to possibly distinguish reservoir and spillover hosts. We used mitochondrial DNA (mtDNA) sequence and nine microsatellites to assess population structure in those two species. mtDNA structure did not correspond to rabies virus variant structure in either species. Microsatellite analyses gave varying results. Bayesian clustering found two groups of arctic foxes in the coastal tundra region, but for red foxes it identified tundra and boreal types. Spatial Bayesian clustering and spatial principal components analysis identified 3 and 4 groups of arctic foxes, respectively, closely matching the distribution of rabies virus variants in the state. Red foxes, conversely, showed eight clusters comprising two regions (boreal and tundra) with much admixture. These results run contrary to previous beliefs that arctic fox show no fine-scale spatial population structure. While we cannot rule out that the red fox is part of the maintenance host community for rabies in Alaska, the distribution of virus variants appears to be driven primarily by the arctic fox. Therefore, we show that host population genetics can be utilized to distinguish between maintenance and spillover hosts when used in conjunction with other approaches.

Fate of the other redcoat: remnants of colonial British foxes in the eastern United States

Red foxes were absent or rare in the southeastern United States until the late 1800s. Their origins potentially include natural population increase/expansion, translocations from Europe, and, eventually, 20th century fur farming. Previous studies have found no European haplotypes in North America, but few samples were sourced from the Atlantic coastal plain, closer to the source of putative introductions. Through analysis of mitochondrial DNA in 584 red foxes from this region, we identified indigenous haplotypes in ≥ 35% of foxes, 1 of 2 European haplotypes in 17% of foxes and fur farm haplotypes in ≥ 13% of foxes; another 35% of foxes had haplotypes potentially indigenous or native. In contrast, only 3 of 135 (2%) male foxes carried a single European Y chromosome haplotype. Most European and fur farm haplotypes were found near the densely human-populated coastal plain and Hudson River lowlands; most red foxes of the Appalachians and Piedmont had native eastern haplotypes. Our findings suggest that the more remote, upland populations primarily reflect indigenous red fox matrilines, whereas urban-associated populations in and around the mid-Atlantic coastal plain and Hudson lowlands reflect an admixture of native and nonnative maternal sources. Autosomal markers are needed to further elucidate the extent of European and fur farm introgression in the Appalachians and further west.

Phylogeography of the Golden Jackal (Canis aureus) in India

The golden jackal (Canis aureus) is one of the most common and widely distributed carnivores in India but phylogeographic studies on the species have been limited across its range. Recent studies have observed absence of mitochondrial (mt) DNA diversity in European populations while some North African populations of golden jackal were found to carry gray wolf (Canis lupus lupaster) mtDNA lineages. In the present study, we sequenced 440 basepairs (bp) of control region (CR) and 412 bp of cytochrome b (cyt b) gene of mtDNA from 62 golden jackals sampled from India (n = 55), Israel (n = 2) and Bulgaria (n = 5), to obtain a total of eighteen haplotypes, comprising sixteen from India and one each from Israel and Bulgaria. Except for three previously described haplotypes represented by one cyt b and one CR haplotype both from India, and one CR haplotype from Bulgaria, all haplotypes identified in this study are new. Genetic diversity was high in golden jackals compared to that reported for other canids in India. Unlike the paraphyletic status of African conspecifics with the gray wolf, the Indian (and other Eurasian) golden jackal clustered in a distinct but shallow monophyletic clade, displaying no evidence of admixture with sympatric and related gray wolf and domestic dog clades in the region. Phylogeographic analyses indicated no clear pattern of genetic structuring of the golden jackal haplotypes and the median joining network revealed a star-shaped polytomy indicative of recent expansion of the species from India. Indian haplotypes were observed to be interior and thus ancestral compared to haplotypes from Europe and Israel, which were peripheral and hence more derived. Molecular tests for demographic expansion confirmed a recent event of expansion of golden jackals in the Indian subcontinent, which can be traced back ~ 37,000 years ago during the late Pleistocene. Our results suggest that golden jackals have had a potentially longer evolutionary history in India than in other parts of the world, although further sampling from Africa, the Middle East and south-east Asia is needed to test this hypothesis.

Mitochondrial Analysis of the Most Basal Canid Reveals Deep Divergence between Eastern and Western North American Gray Foxes (Urocyon spp.) and Ancient Roots in Pleistocene California

Pleistocene aridification in central North America caused many temperate forest-associated vertebrates to split into eastern and western lineages. Such divisions can be cryptic when Holocene expansions have closed the gaps between once-disjunct ranges or when local morphological variation obscures deeper regional divergences. We investigated such cryptic divergence in the gray fox (Urocyon cinereoargenteus), the most basal extant canid in the world. We also investigated the phylogeography of this species and its diminutive relative, the island fox (U. littoralis), in California. The California Floristic Province was a significant source of Pleistocene diversification for a wide range of taxa and, we hypothesized, for the gray fox as well. Alternatively, gray foxes in California potentially reflected a recent Holocene expansion from further south. We sequenced mitochondrial DNA from 169 gray foxes from the southeastern and southwestern United States and 11 island foxes from three of the Channel Islands. We estimated a 1.3% sequence divergence in the cytochrome b gene between eastern and western foxes and used coalescent simulations to date the divergence to approximately 500,000 years before present (YBP), which is comparable to that between recognized sister species within the Canidae. Gray fox samples collected from throughout California exhibited high haplotype diversity, phylogeographic structure, and genetic signatures of a late-Holocene population decline. Bayesian skyline analysis also indicated an earlier population increase dating to the early Wisconsin glaciation (~70,000 YBP) and a root height extending back to the previous interglacial (~100,000 YBP). Together these findings support California's role as a long-term Pleistocene refugium for western Urocyon. Lastly, based both on our results and re-interpretation of those of another study, we conclude that island foxes of the Channel Islands trace their origins to at least 3 distinct female founders from the mainland rather than to a single matriline, as previously suggested.

Beringian sub-refugia revealed in blackfish (Dallia): implications for understanding the effects of Pleistocene glaciations on Beringian taxa and other Arctic aquatic fauna

Background

Pleistocene climatic instability had profound and diverse effects on the distribution and abundance of Arctic organisms revealed by variation in phylogeographic patterns documented in extant Arctic populations. To better understand the effects of geography and paleoclimate on Beringian freshwater fishes, we examined genetic variability in the genus Dallia (blackfish: Esociformes: Esocidae). The genus Dallia groups between one and three nominal species of small, cold- and hypoxia-tolerant freshwater fishes restricted entirely in distribution to Beringia from the Yukon River basin near Fairbanks, Alaska westward including the Kuskokwim River basin and low-lying areas of Western Alaska to the Amguema River on the north side of the Chukotka Peninsula and Mechigmen Bay on the south side of the Chukotka Peninsula. The genus has a non-continuous distribution divided by the Bering Strait and the Brooks Range. We examined the distribution of genetic variation across this range to determine the number and location of potential sub-refugia within the greater Beringian refugium as well as the roles of the Bering land bridge, Brooks Range, and large rivers within Beringia in shaping the current distribution of populations of Dallia. Our analyses were based on DNA sequence data from two nuclear gene introns (S7 and RAG1) and two mitochondrial genome fragments from nineteen sampling locations. These data were examined under genetic clustering and coalescent frameworks to identify sub-refugia within the greater Beringia refugium and to infer the demographic history of different populations of Dallia.

Results

We identified up to five distinct genetic clusters of Dallia. Four of these genetic clusters are present in Alaska: (1) Arctic Coastal Plain genetic cluster found north of the Brooks Range, (2) interior Alaska genetic cluster placed in upstream locations in the Kuskokwim and Yukon river basins, (3) a genetic cluster found on the Seward Peninsula, and (4) a coastal Alaska genetic cluster encompassing downstream Kuskokwim River and Yukon River basin sample locations and samples from Southwest Alaska not in either of these drainages. The Chukotka samples are assigned to their own genetic cluster (5) similar to the coastal Alaska genetic cluster. The clustering and ordination analyses implemented in Discriminant Analysis of Principal Components (DAPC) and STRUCTURE showed mostly concordant groupings and a high degree of differentiation among groups. The groups of sampling locations identified as genetic clusters correspond to geographic areas divided by likely biogeographic barriers including the Brooks Range and the Bering Strait. Estimates of sequence diversity (θ) are highest in the Yukon River and Kuskokwim River drainages near the Bering Sea. We also infer asymmetric migration rates between genetic clusters. The isolation of Dallia on the Arctic Coastal Plain of Alaska is associated with very low estimated migration rates between the coastal Alaska genetic cluster and the Arctic Coastal Plain genetic cluster.

Conclusions

Our results support a scenario with multiple aquatic sub-refugia in Beringia during the Pleistocene and the preservation of that structure in extant populations of Dallia. An inferred historical presence of Dallia across the Bering land bridge explains the similarities in the genetic composition of Dallia in West Beringia and western coastal Alaska. In contrast, historic and contemporary isolation across the Brooks Range shaped the distinctiveness of present day Arctic Coastal Plain Dallia. Overall this study uncovered a high degree of genetic structuring among populations of Dallia supporting the idea of multiple Beringian sub-refugia during the Pleistocene and which appears to be maintained to the present due to the strictly freshwater nature and low dispersal ability of this genus.

Electronic supplementary material

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

Hair and bone as predictors of tissular mercury concentration in the western Alaska red fox, Vulpes vulpes

We evaluated if total mercury (THg) concentrations of keratin-based and bone-based tissues can predict THg concentrations in skeletal muscle, renal medulla, renal cortex, and liver. The THg concentration in matched tissues of 65 red foxes, Vulpes vulpes, from western Alaska was determined. Hair THg concentration had a significant positive correlation with liver, renal medulla, renal cortex, and muscle. The THg concentration for males and females is moderately predictive of THg concentration in the renal cortex and liver for these foxes based on R(2) values (R(2)=0.61 and 0.63, respectively). Bone is weakly predictive of THg concentration in muscle (R(2)=0.40), but not a reliable tissue to predict THg concentration in liver (R(2)=0.24), renal cortex (R(2)=0.35), or renal medulla (R(2)=0.25). These results confirm the potential use of trapped animals, specifically foxes, as useful Arctic sentinel species to inform researchers about patterns in THg levels over time as industrialization of the Arctic continues.