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.

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.

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.

Sequence comparison of prefrontal cortical brain transcriptome from a tame and an aggressive silver fox (Vulpes vulpes).

BACKGROUND: Two strains of the silver fox (Vulpes vulpes), with markedly different behavioral phenotypes, have been developed by long-term selection for behavior. Foxes from the tame strain exhibit friendly behavior towards humans, paralleling the sociability of canine puppies, whereas foxes from the aggressive strain are defensive and exhibit aggression to humans. To understand the genetic differences underlying these behavioral phenotypes fox-specific genomic resources are needed. RESULTS: cDNA from mRNA from pre-frontal cortex of a tame and an aggressive fox was sequenced using the Roche 454 FLX Titanium platform (> 2.5 million reads & 0.9 Gbase of tame fox sequence; >3.3 million reads & 1.2 Gbase of aggressive fox sequence). Over 80% of the fox reads were assembled into contigs. Mapping fox reads against the fox transcriptome assembly and the dog genome identified over 30,000 high confidence fox-specific SNPs. Fox transcripts for approximately 14,000 genes were identified using SwissProt and the dog RefSeq databases. An at least 2-fold expression difference between the two samples (p < 0.05) was observed for 335 genes, fewer than 3% of the total number of genes identified in the fox transcriptome. CONCLUSIONS: Transcriptome sequencing significantly expanded genomic resources available for the fox, a species without a sequenced genome. In a very cost efficient manner this yielded a large number of fox-specific SNP markers for genetic studies and provided significant insights into the gene expression profile of the fox pre-frontal cortex; expression differences between the two fox samples; and a catalogue of potentially important gene-specific sequence variants. This result demonstrates the utility of this approach for developing genomic resources in species with limited genomic information.

Mapping Loci for fox domestication: deconstruction/reconstruction of a behavioral phenotype.

During the second part of the twentieth century, Belyaev selected tame and aggressive foxes (Vulpes vulpes), in an effort known as the "farm-fox experiment", to recapitulate the process of animal domestication. Using these tame and aggressive foxes as founders of segregant backcross and intercross populations we have employed interval mapping to identify a locus for tame behavior on fox chromosome VVU12. This locus is orthologous to, and therefore validates, a genomic region recently implicated in canine domestication. The tame versus aggressive behavioral phenotype was characterized as the first principal component (PC) of a PC matrix made up of many distinct behavioral traits (e.g. wags tail; comes to the front of the cage; allows head to be touched; holds observer's hand with its mouth; etc.). Mean values of this PC for F1, backcross and intercross populations defined a linear gradient of heritable behavior ranging from tame to aggressive. The second PC did not follow such a gradient, but also mapped to VVU12, and distinguished between active and passive behaviors. These data suggest that (1) there are at least two VVU12 loci associated with behavior; (2) expression of these loci is dependent on interactions with other parts of the genome (the genome context) and therefore varies from one crossbred population to another depending on the individual parents that participated in the cross.

A meiotic linkage map of the silver fox, aligned and compared to the canine genome.

A meiotic linkage map is essential for mapping traits of interest and is often the first step toward understanding a cryptic genome. Specific strains of silver fox (a variant of the red fox, Vulpes vulpes), which segregate behavioral and morphological phenotypes, create a need for such a map. One such strain, selected for docility, exhibits friendly dog-like responses to humans, in contrast to another strain selected for aggression. Development of a fox map is facilitated by the known cytogenetic homologies between the dog and fox, and by the availability of high resolution canine genome maps and sequence data. Furthermore, the high genomic sequence identity between dog and fox allows adaptation of canine microsatellites for genotyping and meiotic mapping in foxes. Using 320 such markers, we have constructed the first meiotic linkage map of the fox genome. The resulting sex-averaged map covers 16 fox autosomes and the X chromosome with an average inter-marker distance of 7.5 cM. The total map length corresponds to 1480.2 cM. From comparison of sex-averaged meiotic linkage maps of the fox and dog genomes, suppression of recombination in pericentromeric regions of the metacentric fox chromosomes was apparent, relative to the corresponding segments of acrocentric dog chromosomes. Alignment of the fox meiotic map against the 7.6x canine genome sequence revealed high conservation of marker order between homologous regions of the two species. The fox meiotic map provides a critical tool for genetic studies in foxes and identification of genetic loci and genes implicated in fox domestication.