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.

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

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

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

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

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.

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

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

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

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

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

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

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.

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

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

Unbiased Estimation of Linkage Disequilibrium from Unphased Data.

Linkage disequilibrium (LD) is used to infer evolutionary history, to identify genomic regions under selection, and to dissect the relationship between genotype and phenotype. In each case, we require accurate estimates of LD statistics from sequencing data. Unphased data present a challenge because multilocus haplotypes cannot be inferred exactly. Widely used estimators for the common statistics r2 and D2 exhibit large and variable upward biases that complicate interpretation and comparison across cohorts. Here, we show how to find unbiased estimators for a wide range of two-locus statistics, including D2, for both single and multiple randomly mating populations. These unbiased statistics are particularly well suited to estimate effective population sizes from unlinked loci in small populations. We develop a simple inference pipeline and use it to refine estimates of recent effective population sizes of the threatened Channel Island Fox populations.

Fur colour in the Arctic fox: genetic architecture and consequences for fitness.

Genome-wide association studies provide good opportunities for studying the genetic basis of adaptive traits in wild populations. Yet, previous studies often failed to identify major effect genes. In this study, we used high-density single nucleotide polymorphism and individual fitness data from a wild non-model species. Using a whole-genome approach, we identified the MC1R gene as the sole causal gene underlying Arctic fox Vulpes lagopus fur colour. Further, we showed the adaptive importance of fur colour genotypes through measures of fitness that link ecological and evolutionary processes. We found a tendency for blue foxes that are heterozygous at the fur colour locus to have higher fitness than homozygous white foxes. The effect of genotype on fitness was independent of winter duration but varied with prey availability, with the strongest effect in years of increasing rodent populations. MC1R is located in a genomic region with high gene density, and we discuss the potential for indirect selection through linkage and pleiotropy. Our study shows that whole-genome analyses can be successfully applied to wild species and identify major effect genes underlying adaptive traits. Furthermore, we show how this approach can be used to identify knowledge gaps in our understanding of interactions between ecology and evolution.

Comparison of intra- and inter-host genetic diversity in rabies virus during experimental cross-species transmission.

The development of high-throughput genome sequencing enables accurate measurements of levels of sub-consensus intra-host virus genetic diversity and analysis of the role played by natural selection during cross-species transmission. We analysed the natural and experimental evolution of rabies virus (RABV), an important example of a virus that is able to make multiple host jumps. In particular, we (i) analyzed RABV evolution during experimental host switching with the goal of identifying possible genetic markers of host adaptation, (ii) compared the mutational changes observed during passage with those observed in natura, and (iii) determined whether the colonization of new hosts or tissues requires adaptive evolution in the virus. To address these aims, animal infection models (dog and fox) and primary cell culture models (embryo brain cells of dog and fox) were developed and viral variation was studied in detail through deep genome sequencing. Our analysis revealed a strong unidirectional host evolutionary effect, as dog-adapted rabies virus was able to replicate in fox and fox cells relatively easily, while dogs or neuronal dog cells were not easily susceptible to fox adapted-RABV. This suggests that dog RABV may be able to adapt to some hosts more easily than other host variants, or that when RABV switched from dogs to red foxes it lost its ability to adapt easily to other species. Although no difference in patterns of mutation variation between different host organs was observed, mutations were common following both in vitro and in vivo passage. However, only a small number of these mutations also appeared in natura, suggesting that adaptation during successful cross-species virus transmission is a complex, multifactorial evolutionary process.

The biological characteristics of the canine adenovirus type 1 from fox and the transcriptome analysis of the infected MDCK cell.

Canine adenovirus type 1 (CAdV-1) is the etiologic agent of fox encephalitis, and a virus strain from fox encephalitis is isolated and related research are conducted. In this experiment, the results showed that the F1301 strain was confirmed to be the CAdV-1. The whole genome of the CAdV-1 F1301 strain isolated from fox was 30,535 bp and had higher homology to the other reported CAdV-1 strains. After 0, 12, and 36 h of CAdV-1 infection, the difference gene of the 592 long noncoding RNA and 11,215 microRNA were involved in cell responses to CAdV-1 infection through the PI3K-AKT, Wnt, Herpes simplex, hepatitis C, and Epstein-Barr virus infection pathway in Madin-Darby canine kidney cell line (MDCK). The results indicate that the biological characterization of the CAdV-1 and the MDCK cell-CAdV-1 interaction are clarified.

Molecular detection of coccidian Apicomplexa Parasites isolated from wild crab-eating and pampas foxes through novel TaqMan™ probes: a contribution to their molecular epidemiology.

Neospora caninum, Toxoplasma gondii and Hammondia spp. are coccidian parasites similar in morphology. Molecular techniques are necessary to detect parasite DNA isolated from stool samples in wild canids because they were reported as definitive hosts of N. caninum life cycle. The objective of this study was to develop a highly sensitive and accurate molecular method for the identification of coccidian Apicomplexa parasites in crab-eating fox (Cerdocyon thous) and pampas fox (Lycalopex gymnocercus). Tissue samples from road-killed animals (pampas fox = 46, crab-eating fox = 55) and feces (pampas fox = 84, crab-eating fox = 2) were collected, and species were diagnosed through molecular assay. PCR was used for the amplification of a fragment of the coccidian Apicomplexa nss-rRNA gene. Additionally, we developed a novel real-time PCR TaqMan™ probe approach to detect T. gondii- Hammondia spp. and N. caninum. This is the first report of N. caninum DNA in pampas fox feces (n = 1), thus it was also detected from pampas fox tissues (n = 1). Meanwhile, T. gondii was found in tissues of pampas (n = 1) and crab-eating (n = 1) foxes and H. triffittae in one crab-eating fox tissue. Despite the low percentage (2.5%) of positive samples, the molecular method developed in this study proved to be highly sensitive and accurate allowing to conduct an extensive monitoring analysis for these parasites in wildlife.

Low Persistence of Genetic Rescue Across Generations in the Arctic Fox (Vulpes lagopus).

Genetic rescue can facilitate the recovery of small and isolated populations suffering from inbreeding depression. Long-term effects are however complex, and examples spanning over multiple generations under natural conditions are scarce. The aim of this study was to test for long-term effects of natural genetic rescue in a small population of Scandinavian Arctic foxes (Vulpes lagopus). By combining a genetically verified pedigree covering almost 20 years with a long-term dataset on individual fitness (n = 837 individuals), we found no evidence for elevated fitness in immigrant F2 and F3 compared to native inbred foxes. Population inbreeding levels showed a fluctuating increasing trend and emergence of inbreeding within immigrant lineages shortly after immigration. Between 0-5 and 6-9 years post immigration, the average number of breeding adults decreased by almost 22% and the average proportion of immigrant ancestry rose from 14% to 27%. Y chromosome analysis revealed that 2 out of 3 native male lineages were lost from the gene pool, but all founders represented at the time of immigration were still contributing to the population at the end of the study period through female descendants. The results highlight the complexity of genetic rescue and suggest that beneficial effects can be brief. Continuous gene flow may be needed for small and threatened populations to recover and persist in a longer time perspective.

Genomic responses to selection for tame/aggressive behaviors in the silver fox (Vulpes vulpes).

Animal domestication efforts have led to a shared spectrum of striking behavioral and morphological changes. To recapitulate this process, silver foxes have been selectively bred for tame and aggressive behaviors for more than 50 generations at the Institute for Cytology and Genetics in Novosibirsk, Russia. To understand the genetic basis and molecular mechanisms underlying the phenotypic changes, we profiled gene expression levels and coding SNP allele frequencies in two brain tissue specimens from 12 aggressive foxes and 12 tame foxes. Expression analysis revealed 146 genes in the prefrontal cortex and 33 genes in the basal forebrain that were differentially expressed, with a 5% false discovery rate (FDR). These candidates include genes in key pathways known to be critical to neurologic processing, including the serotonin and glutamate receptor pathways. In addition, 295 of the 31,000 exonic SNPs show significant allele frequency differences between the tame and aggressive populations (1% FDR), including genes with a role in neural crest cell fate determination.

Beyond the landscape: Resistance modelling infers physical and behavioural gene flow barriers to a mobile carnivore across a metropolitan area.

Urbanization affects key aspects of wildlife ecology. Dispersal in urban wildlife species may be impacted by geographical barriers but also by a species' inherent behavioural variability. There are no functional connectivity analyses using continuous individual-based sampling across an urban-rural continuum that would allow a thorough assessment of the relative importance of physical and behavioural dispersal barriers. We used 16 microsatellite loci to genotype 374 red foxes (Vulpes vulpes) from the city of Berlin and surrounding rural regions in Brandenburg in order to study genetic structure and dispersal behaviour of a mobile carnivore across the urban-rural landscape. We assessed functional connectivity by applying an individual-based landscape genetic optimization procedure. Three commonly used genetic distance measures yielded different model selection results, with only the results of an eigenvector-based multivariate analysis reasonably explaining genetic differentiation patterns. Genetic clustering methods and landscape resistance modelling supported the presence of an urban population with reduced dispersal across the city border. Artificial structures (railways, motorways) served as main dispersal corridors within the cityscape, yet urban foxes avoided densely built-up areas. We show that despite their ubiquitous presence in urban areas, their mobility and behavioural plasticity, foxes were affected in their dispersal by anthropogenic presence. Distinguishing between man-made structures and sites of human activity, rather than between natural and artificial structures, is thus essential for better understanding urban fox dispersal. This differentiation may also help to understand dispersal of other urban wildlife and to predict how behaviour can shape population genetic structure beyond physical barriers.

A Genome-Wide Perspective on the Persistence of Red Wolf Ancestry in Southeastern Canids.

The red wolf (Canis rufus), a legally recognized and critically endangered wolf, is known to interbreed with coyotes (Canis latrans). Declared extirpated in the wild in 1980, red wolves were reintroduced to northeastern North Carolina nearly a decade later. Interbreeding with coyotes was thought to be restricted to a narrow geographic region adjacent to the reintroduced population and largely believed to threaten red wolf recovery. However, red wolf ancestry was recently discovered in canids along the American Gulf Coast, igniting a broader survey of ancestry in southeastern canid populations. Here, we examine geographic and temporal patterns of genome-wide red wolf ancestry in 260 canids across the southeastern United States at over 164 000 SNP loci. We found that red wolf ancestry was most prevalent in canids sampled from Texas in the mid-1970s, although non-trivial amounts of red wolf ancestry persist in this region today. Further, red wolf ancestry was also observed in a subset of coyotes inhabiting North Carolina, despite management efforts to limit the occurrence of hybridization events. Lastly, we found no evidence of substantial red wolf ancestry in southeastern canids outside of these 2 admixture zones. Overall, this study provides a genome-wide survey of red wolf ancestry in canids across the southeastern United States, which may ultimately inform future red wolf restoration efforts.

Comparative Study of the Bush Dog (Speothos venaticus) Karyotype and Analysis of Satellite DNA Sequences and Their Chromosome Distribution in Six Species of Canidae.

The bush dog (Speothos venaticus, 2n = 74) is a near threatened species taxonomically classified among South American canids. We revised the bush dog karyotype and performed a comparative sequence analysis of satellite and satellite-like DNAs in 6 canids: the bush dog, domestic dog (Canis familiaris, 2n = 78), grey wolf (C. lupus, 2n = 78), Chinese raccoon dog (Nyctereutes procyonoides procyonoides, 2n = 54+B), red fox (Vulpes vulpes, 2n = 34+B), and arctic fox (V. lagopus, 2n = 48-50) to specify the species position among Canidae. Using FISH with painting and BAC probes, we found that the distribution of canid evolutionarily conserved chromosome segments in the bush dog karyotype is similar to that of the domestic dog and grey wolf. The bush dog karyotype differs by 2 acrocentric chromosome pairs formed by tandem fusions of the canine (29;34) and (26;35) orthologues. An interstitial signal of the telomeric probe was observed in the (26;35) fusion site in the bush dog indicating a recent evolutionary origin of this rearrangement. Sequences and hybridisation patterns of satellite DNAs were compared, and a phylogenetic tree of the 6 canid species was constructed which confirmed the bush dog position close to the wolf-like canids, and apart from the raccoon dog and foxes.

Using Ancestry-Informative SNPs to Quantify Introgression of European Alleles into North American Red Foxes.

A recent study demonstrated that British red foxes introduced to the mid-Atlantic coastal plain (ACP) of the eastern United States during the late 18th century successfully interbred with indigenous American red foxes despite half a million year's divergence. However, a large disparity in frequency of European mitochondria (27%) versus Y chromosomes (1%) left unclear the magnitude of genetic exchange. We sought to quantify genomic introgression using 35 autosomal and 5 X-chromosome ancestry-informative markers (AIMs) in conjunction with diagnostic Y chromosome single nucleotide polymorphism (Y-SNP) markers to characterize the modern state of red foxes in the eastern United States and to gain insight into the potential role of reproductive barriers. European admixture was highest in the ACP and apparently restricted to the central eastern United States. We estimated only slightly (and nonsignificantly) European ancestry in autosomal than X-chromosome markers. European ancestry from autosomal and X-chromosome markers (36.4%) was higher than the corresponding mitochondrial (mt) DNA estimate (26.4%) in the ACP. Only 1 of 124 males (<1%) in the ACP had European Y chromosomes, which was similar to the neighboring regions, in which 2 of 99 (2%) males carried a European Y chromosome (the same haplotype). Although we could not rule out drift as the cause of low European Y-chromosome frequency, results were also consistent with F1 male infertility. In the future, more extensive genomic sequencing will enable a more thorough investigation of possible barrier genes on the X chromosome as well as throughout the genome.