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.

The brain of the silver fox (Vulpes vulpes): a neuroanatomical reference of cell-stained histological and MRI images.

Although the silver fox (Vulpes vulpes) has been largely overlooked by neuroscientists, it has the potential to serve as a powerful model for the investigation of brain-behavior relationships. The silver fox is a melanistic variant of the red fox. Within this species, the long-running Russian farm-fox experiment has resulted in different strains bred to show divergent behavior. Strains bred for tameness, aggression, or without selection on behavior present an excellent opportunity to investigate neuroanatomical changes underlying behavioral characteristics. Here, we present a histological and MRI neuroanatomical reference of a fox from the conventional strain, which is bred without behavioral selection. This can provide an anatomical basis for future studies of the brains of foxes from this particular experiment, as well as contribute to an understanding of fox brains in general. In addition, this can serve as a resource for comparative neuroscience and investigations into neuroanatomical variation among the family Canidae, the order Carnivora, and mammals more broadly.

Defining hippocampal area CA2 in the fox (Vulpes vulpes) brain.

Since 1959, the Russian Farm-Fox study has bred foxes to be either tame or, more recently, aggressive, and scientists have used them to gain insight into the brain structures associated with these behavioral features. In mice, hippocampal area CA2 has emerged as one of the essential regulators of social aggression, and so to eventually determine whether we could identify differences in CA2 between tame and aggressive foxes, we first sought to identify CA2 in foxes (Vulpes vulpes). As no clearly defined area of CA2 has been described in species such as cats, dogs, or pigs, it was not at all clear whether CA2 could be identified in foxes. In this study, we cut sections of temporal lobes from male and female red foxes, perpendicular to the long axis of the hippocampus, and stained them with markers of CA2 pyramidal cells commonly used in tissue from rats and mice. We observed that antibodies against Purkinje cell protein 4 best stained the pyramidal cells in the area spanning the end of the mossy fibers and the beginning of the pyramidal cells lacking mossy fibers, resembling the pattern seen in rats and mice. Our findings indicate that foxes do have a "molecularly defined" CA2, and further, they suggest that other carnivores like dogs and cats might as well. With this being the case, these foxes could be useful in future studies looking at CA2 as it relates to aggression.

Neuromorphological changes following selection for tameness and aggression in the Russian fox-farm experiment.

The Russian fox-farm experiment is an unusually long-running and well-controlled study designed to replicate wolf-to-dog domestication. As such, it offers an unprecedented window onto the neural mechanisms governing the evolution of behavior. Here we report evolved changes to gray matter morphology resulting from selection for tameness vs. aggressive responses toward humans in a sample of 30 male fox brains. Contrasting with standing ideas on the effects of domestication on brain size, tame foxes did not show reduced brain volume. Rather, gray matter volume in both the tame and aggressive strains was increased relative to conventional farm foxes bred without deliberate selection on behavior. Furthermore, tame- and aggressive-enlarged regions overlapped substantially, including portions of motor, somatosensory, and prefrontal cortex, amygdala, hippocampus, and cerebellum. We also observed differential morphological covariation across distributed gray matter networks. In one prefrontal-cerebellum network, this covariation differentiated the three populations along the tame-aggressive behavioral axis. Surprisingly, a prefrontal-hypothalamic network differentiated the tame and aggressive foxes together from the conventional strain. These findings indicate that selection for opposite behaviors can influence brain morphology in a similar way.SIGNIFICANCE STATEMENTDomestication represents one of the largest and most rapid evolutionary shifts of life on earth. However, its neural correlates are largely unknown. Here we report the neuroanatomical consequences of selective breeding for tameness or aggression in the seminal Russian fox-farm experiment. Compared to a population of conventional farm-bred control foxes, tame foxes show neuroanatomical changes in the prefrontal cortex and hypothalamus, paralleling wolf-to-dog shifts. Surprisingly, though, aggressive foxes also show similar changes. Moreover, both strains show increased gray matter volume relative to controls. These results indicate that similar brain adaptations can result from selection for opposite behavior, that existing ideas of brain changes in domestication may need revision, and that significant neuroanatomical change can evolve very quickly - within the span of less than a hundred generations.

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.

Hypothalamic transcriptome of tame and aggressive silver foxes (Vulpes vulpes) identifies gene expression differences shared across brain regions.

The underlying neurological events accompanying dog domestication remain elusive. To reconstruct the domestication process in an experimental setting, silver foxes (Vulpes vulpes) have been deliberately bred for tame vs aggressive behaviors for more than 50 generations at the Institute for Cytology and Genetics in Novosibirsk, Russia. The hypothalamus is an essential part of the hypothalamic-pituitary-adrenal axis and regulates the fight-or-flight response, and thus, we hypothesized that selective breeding for tameness/aggressiveness has shaped the hypothalamic transcriptomic profile. RNA-seq analysis identified 70 differentially expressed genes (DEGs). Seven of these genes, DKKL1, FBLN7, NPL, PRIMPOL, PTGRN, SHCBP1L and SKIV2L, showed the same direction expression differences in the hypothalamus, basal forebrain and prefrontal cortex. The genes differentially expressed across the three tissues are involved in cell division, differentiation, adhesion and carbohydrate processing, suggesting an association of these processes with selective breeding. Additionally, 159 transcripts from the hypothalamus demonstrated differences in the abundance of alternative spliced forms between the tame and aggressive foxes. Weighted gene coexpression network analyses also suggested that gene modules in hypothalamus were significantly associated with tame vs aggressive behavior. Pathways associated with these modules include signal transduction, interleukin signaling, cytokine-cytokine receptor interaction and peptide ligand-binding receptors (eg, G-protein coupled receptor [GPCR] ligand binding). Current studies show the selection for tameness vs aggressiveness in foxes is associated with unique hypothalamic gene profiles partly shared with other brain regions and highlight DEGs involved in biological processes such as development, differentiation and immunological responses. The role of these processes in fox and dog domestication remains to be determined.

The Red Fox Y-Chromosome in Comparative Context.

While the number of mammalian genome assemblies has proliferated, Y-chromosome assemblies have lagged behind. This discrepancy is caused by biological features of the Y-chromosome, such as its high repeat content, that present challenges to assembly with short-read, next-generation sequencing technologies. Partial Y-chromosome assemblies have been developed for the cat (Feliscatus), dog (Canislupusfamiliaris), and grey wolf (Canislupuslupus), providing the opportunity to examine the red fox (Vulpesvulpes) Y-chromosome in the context of closely related species. Here we present a data-driven approach to identifying Y-chromosome sequence among the scaffolds that comprise the short-read assembled red fox genome. First, scaffolds containing genes found on the Y-chromosomes of cats, dogs, and wolves were identified. Next, analysis of the resequenced genomes of 15 male and 15 female foxes revealed scaffolds containing male-specific k-mers and patterns of inter-sex copy number variation consistent with the heterogametic chromosome. Analyzing variation across these two metrics revealed 171 scaffolds containing 3.37 Mbp of putative Y-chromosome sequence. The gene content of these scaffolds is consistent overall with that of the Y-chromosome in other carnivore species, though the red fox Y-chromosome carries more copies of BCORY2 and UBE1Y than has been reported in related species and fewer copies of SRY than in other canids. The assignment of these scaffolds to the Y-chromosome serves to further characterize the content of the red fox draft genome while providing resources for future analyses of canid Y-chromosome evolution.

Evolution of gene regulation in ruminants differs between evolutionary breakpoint regions and homologous synteny blocks.

The role of chromosome rearrangements in driving evolution has been a long-standing question of evolutionary biology. Here we focused on ruminants as a model to assess how rearrangements may have contributed to the evolution of gene regulation. Using reconstructed ancestral karyotypes of Cetartiodactyls, Ruminants, Pecorans, and Bovids, we traced patterns of gross chromosome changes. We found that the lineage leading to the ruminant ancestor after the split from other cetartiodactyls was characterized by mostly intrachromosomal changes, whereas the lineage leading to the pecoran ancestor (including all livestock ruminants) included multiple interchromosomal changes. We observed that the liver cell putative enhancers in the ruminant evolutionary breakpoint regions are highly enriched for DNA sequences under selective constraint acting on lineage-specific transposable elements (TEs) and a set of 25 specific transcription factor (TF) binding motifs associated with recently active TEs. Coupled with gene expression data, we found that genes near ruminant breakpoint regions exhibit more divergent expression profiles among species, particularly in cattle, which is consistent with the phylogenetic origin of these breakpoint regions. This divergence was significantly greater in genes with enhancers that contain at least one of the 25 specific TF binding motifs and located near bovidae-to-cattle lineage breakpoint regions. Taken together, by combining ancestral karyotype reconstructions with analysis of cis regulatory element and gene expression evolution, our work demonstrated that lineage-specific regulatory elements colocalized with gross chromosome rearrangements may have provided valuable functional modifications that helped to shape ruminant evolution.

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.

Author Correction: Red fox genome assembly identifies genomic regions associated with tame and aggressive behaviours.

In the version of this Article originally published, there were some errors in the affiliations: Stephen J. O'Brien's affiliations were incorrectly listed as 8,9; they should have been 7,9. Affiliation 3 was incorrectly named the Institute of Cytology and Genetics of the Russian Academy of Sciences; it should have read Institute of Cytology and Genetics of the Siberian Branch of the Russian Academy of Sciences. Affiliation 4 was incorrectly named the Institute of Molecular and Cell Biology of the Russian Academy of Sciences; it should have read Institute of Molecular and Cellular Biology of the Siberian Branch of the Russian Academy of Sciences. These have now been corrected.

Sequencing of Supernumerary Chromosomes of Red Fox and Raccoon Dog Confirms a Non-Random Gene Acquisition by B Chromosomes.

B chromosomes (Bs) represent a variable addition to the main karyotype in some lineages of animals and plants. Bs accumulate through non-Mendelian inheritance and become widespread in populations. Despite the presence of multiple genes, most Bs lack specific phenotypic effects, although their influence on host genome epigenetic status and gene expression are recorded. Previously, using sequencing of isolated Bs of ruminants and rodents, we demonstrated that Bs originate as segmental duplications of specific genomic regions, and subsequently experience pseudogenization and repeat accumulation. Here, we used a similar approach to characterize Bs of the red fox (Vulpes vulpes L.) and the Chinese raccoon dog (Nyctereutes procyonoides procyonoides Gray). We confirm the previous findings of the KIT gene on Bs of both species, but demostrate an independent origin of Bs in these species, with two reused regions. Comparison of gene ensembles in Bs of canids, ruminants, and rodents once again indicates enrichment with cell-cycle genes, development-related genes, and genes functioning in the neuron synapse. The presence of B-chromosomal copies of genes involved in cell-cycle regulation and tissue differentiation may indicate importance of these genes for B chromosome establishment.

Red fox genome assembly identifies genomic regions associated with tame and aggressive behaviours.

Strains of red fox (Vulpes vulpes) with markedly different behavioural phenotypes have been developed in the famous long-term selective breeding programme known as the Russian farm-fox experiment. Here we sequenced and assembled the red fox genome and re-sequenced a subset of foxes from the tame, aggressive and conventional farm-bred populations to identify genomic regions associated with the response to selection for behaviour. Analysis of the re-sequenced genomes identified 103 regions with either significantly decreased heterozygosity in one of the three populations or increased divergence between the populations. A strong positional candidate gene for tame behaviour was highlighted: SorCS1, which encodes the main trafficking protein for AMPA glutamate receptors and neurexins and suggests a role for synaptic plasticity in fox domestication. Other regions identified as likely to have been under selection in foxes include genes implicated in human neurological disorders, mouse behaviour and dog domestication. The fox represents a powerful model for the genetic analysis of affiliative and aggressive behaviours that can benefit genetic studies of behaviour in dogs and other mammals, including humans.

Construction of Red Fox Chromosomal Fragments from the Short-Read Genome Assembly.

The genome of a red fox (Vulpes vulpes) was recently sequenced and assembled using next-generation sequencing (NGS). The assembly is of high quality, with 94X coverage and a scaffold N50 of 11.8 Mbp, but is split into 676,878 scaffolds, some of which are likely to contain assembly errors. Fragmentation and misassembly hinder accurate gene prediction and downstream analysis such as the identification of loci under selection. Therefore, assembly of the genome into chromosome-scale fragments was an important step towards developing this genomic model. Scaffolds from the assembly were aligned to the dog reference genome and compared to the alignment of an outgroup genome (cat) against the dog to identify syntenic sequences among species. The program Reference-Assisted Chromosome Assembly (RACA) then integrated the comparative alignment with the mapping of the raw sequencing reads generated during assembly against the fox scaffolds. The 128 sequence fragments RACA assembled were compared to the fox meiotic linkage map to guide the construction of 40 chromosomal fragments. This computational approach to assembly was facilitated by prior research in comparative mammalian genomics, and the continued improvement of the red fox genome can in turn offer insight into canid and carnivore chromosome evolution. This assembly is also necessary for advancing genetic research in foxes and other canids.

Anterior Pituitary Transcriptome Suggests Differences in ACTH Release in Tame and Aggressive Foxes.

Domesticated species exhibit a suite of behavioral, endocrinological, and morphological changes referred to as "domestication syndrome." These changes may include a reduction in reactivity of the hypothalamic-pituitary-adrenal (HPA) axis and specifically reduced adrenocorticotropic hormone release from the anterior pituitary. To investigate the biological mechanisms targeted during domestication, we investigated gene expression in the pituitaries of experimentally domesticated foxes (Vulpes vulpes). RNA was sequenced from the anterior pituitary of six foxes selectively bred for tameness ("tame foxes") and six foxes selectively bred for aggression ("aggressive foxes"). Expression, splicing, and network differences identified between the two lines indicated the importance of genes related to regulation of exocytosis, specifically mediated by cAMP, organization of pseudopodia, and cell motility. These findings provide new insights into biological mechanisms that may have been targeted when these lines of foxes were selected for behavior and suggest new directions for research into HPA axis regulation and the biological underpinnings of domestication.

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

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

X Chromosome Evolution in Cetartiodactyla.

The phenomenon of a remarkable conservation of the X chromosome in eutherian mammals has been first described by Susumu Ohno in 1964. A notable exception is the cetartiodactyl X chromosome, which varies widely in morphology and G-banding pattern between species. It is hypothesized that this sex chromosome has undergone multiple rearrangements that changed the centromere position and the order of syntenic segments over the last 80 million years of Cetartiodactyla speciation. To investigate its evolution we have selected 26 evolutionarily conserved bacterial artificial chromosome (BAC) clones from the cattle CHORI-240 library evenly distributed along the cattle X chromosome. High-resolution BAC maps of the X chromosome on a representative range of cetartiodactyl species from different branches: pig (Suidae), alpaca (Camelidae), gray whale (Cetacea), hippopotamus (Hippopotamidae), Java mouse-deer (Tragulidae), pronghorn (Antilocapridae), Siberian musk deer (Moschidae), and giraffe (Giraffidae) were obtained by fluorescent in situ hybridization. To trace the X chromosome evolution during fast radiation in specious families, we performed mapping in several cervids (moose, Siberian roe deer, fallow deer, and Pere David's deer) and bovid (muskox, goat, sheep, sable antelope, and cattle) species. We have identified three major conserved synteny blocks and rearrangements in different cetartiodactyl lineages and found that the recently described phenomenon of the evolutionary new centromere emergence has taken place in the X chromosome evolution of Cetartiodactyla at least five times. We propose the structure of the putative ancestral cetartiodactyl X chromosome by reconstructing the order of syntenic segments and centromere position for key groups.

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.

Genetics of Interactive Behavior in Silver Foxes (Vulpes vulpes).

Individuals involved in a social interaction exhibit different behavioral traits that, in combination, form the individual's behavioral responses. Selectively bred strains of silver foxes (Vulpes vulpes) demonstrate markedly different behaviors in their response to humans. To identify the genetic basis of these behavioral differences we constructed a large F2 population including 537 individuals by cross-breeding tame and aggressive fox strains. 98 fox behavioral traits were recorded during social interaction with a human experimenter in a standard four-step test. Patterns of fox behaviors during the test were evaluated using principal component (PC) analysis. Genetic mapping identified eight unique significant and suggestive QTL. Mapping results for the PC phenotypes from different test steps showed little overlap suggesting that different QTL are involved in regulation of behaviors exhibited in different behavioral contexts. Many individual behavioral traits mapped to the same genomic regions as PC phenotypes. This provides additional information about specific behaviors regulated by these loci. Further, three pairs of epistatic loci were also identified for PC phenotypes suggesting more complex genetic architecture of the behavioral differences between the two strains than what has previously been observed.