Anastasiya V. Vladimirova

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.

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

BackgroundTwo 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.ResultscDNA 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.ConclusionsTranscriptome 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.

Genotyping-By-Sequencing (GBS) Detects Genetic Structure and Confirms Behavioral QTL in Tame and Aggressive Foxes (Vulpes vulpes)

The silver fox (Vulpes vulpes) offers a novel model for studying the genetics of social behavior and animal domestication. Selection of foxes, separately, for tame and for aggressive behavior has yielded two strains with markedly different, genetically determined, behavioral phenotypes. Tame strain foxes are eager to establish human contact while foxes from the aggressive strain are aggressive and difficult to handle. These strains have been maintained as separate outbred lines for over 40 generations but their genetic structure has not been previously investigated. We applied a genotyping-by-sequencing (GBS) approach to provide insights into the genetic composition of these fox populations. Sequence analysis of EcoT22I genomic libraries of tame and aggressive foxes identified 48,294 high quality SNPs. Population structure analysis revealed genetic divergence between the two strains and more diversity in the aggressive strain than in the tame one. Significant differences in allele frequency between the strains were identified for 68 SNPs. Three of these SNPs were located on fox chromosome 14 within an interval of a previously identified behavioral QTL, further supporting the importance of this region for behavior. The GBS SNP data confirmed that significant genetic diversity has been preserved in both fox populations despite many years of selective breeding. Analysis of SNP allele frequencies in the two populations identified several regions of genetic divergence between the tame and aggressive foxes, some of which may represent targets of selection for behavior. The GBS protocol used in this study significantly expanded genomic resources for the fox, and can be adapted for SNP discovery and genotyping in other canid species.

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.

Platinum coat color in red fox (Vulpes vulpes) is caused by a mutation in an autosomal copy of KIT

The red fox (Vulpes vulpes) demonstrates a variety of coat colors including platinum, a common phenotype maintained in farm-bred fox populations. Foxes heterozygous for the platinum allele have a light silver coat and extensive white spotting, whereas homozygosity is embryonic lethal. Two KIT transcripts were identified in skin cDNA from platinum foxes. The long transcript was identical to the KIT transcript of silver foxes, whereas the short transcript, which lacks exon 17, was specific to platinum. The KIT gene has several copies in the fox genome: an autosomal copy on chromosome 2 and additional copies on the B chromosomes. To identify the platinum-specific KIT sequence, the genomes of one platinum and one silver fox were sequenced. A single nucleotide polymorphism (SNP) was identified at the first nucleotide of KIT intron 17 in the platinum fox. In platinum foxes, the A allele of the SNP disrupts the donor splice site and causes exon 17, which is part of a segment that encodes a conserved tyrosine kinase domain, to be skipped. Complete cosegregation of the A allele with the platinum phenotype was confirmed by linkage mapping (LOD 25.59). All genotyped farm-bred platinum foxes from Russia and the US were heterozygous for the SNP (A/G), whereas foxes with different coat colors were homozygous for the G allele. Identification of the platinum mutation suggests that other fox white-spotting phenotypes, which are allelic to platinum, would also be caused by mutations in the KIT gene.

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.Long-term selective breeding has produced strains of the red fox (Vulpes vulpes) with different behaviours. Here, the authors sequence the genomes of tame and aggressive strains to uncover the genetic regions that have responded to selection for behaviour.

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.

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

Significance The behavior of domesticated animals differs dramatically from that of wild relatives, and the Russian tame fox experiment demonstrated clearly that these changes can occur in just a few generations of selection. Analysis of gene expression in the brains of tame and aggressive foxes from this experiment allows us to ask what brain pathways have been altered by this recent, strong selection. Pathways that impact the function of both serotonergic and glutaminergic neurons were clearly modulated by selection, consistent with the roles of these neurons in learning and memory. Both allele frequency and gene expression changes also implicate genes important in neural crest cell function, supporting a possible role of neural crest cells in the domestication syndrome. 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.