Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Cover ArticleFeatured ArticleResearch Articles, Behavioral/Cognitive

Neuromorphological Changes following Selection for Tameness and Aggression in the Russian Farm-Fox experiment

Erin E. Hecht, Anna V. Kukekova, David A. Gutman, Gregory M. Acland, Todd M. Preuss and Lyudmila N. Trut
Journal of Neuroscience 14 July 2021, 41 (28) 6144-6156; DOI: https://doi.org/10.1523/JNEUROSCI.3114-20.2021
Erin E. Hecht
1Department of Human Evolutionary Biology, Harvard University, Cambridge, MA 02138
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Erin E. Hecht
Anna V. Kukekova
2Department of Animal Sciences, College of Agriculture, Consumer, and Environmental Sciences, University of IL Urbana-Champaign, Urbana, IL 61801
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David A. Gutman
3School of Medicine, Emory University, Atlanta, GA 30322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gregory M. Acland
4Baker Institute for Animal Health, College of Veterinary Medicine, Cornell University, Ithaca, New York, 14853
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Todd M. Preuss
5Yerkes National Primate Research Center, Emory University, Atlanta, GA 30322
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Todd M. Preuss
Lyudmila N. Trut
6Institute of Cytology and Genetics, Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia, 630090
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

The Russian farm-fox 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 versus 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 morphologic 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 STATEMENT Domestication 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 farm-fox experiment. Compared with a population of conventional farm-bred control foxes, tame foxes show neuroanatomical changes in the PFC 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 <100 generations.

  • aggression
  • canids
  • domestication
  • farm-fox experiment
  • neuroimaging
  • social behavior

Introduction

Domestication refers to the process of animal adaptation to the human niche. It represents one of the largest and most rapid evolutionary shifts in life on Earth: the biomass of domesticated animals has increased an estimated 3.5-fold in the last 100 years and now outweighs the biomass of other terrestrial mammals by a factor of ∼25 (Smil, 2011). Correspondingly, the neural changes associated with domestication constitute a major event in the history of brain evolution. Moreover, self-domestication is hypothesized to have played a role in the evolution of our own species (Brune, 2007; Hare, 2017; Wrangham, 2018). However, surprisingly little is known about the neural correlates of domestication.

Perhaps the most well-known effect of domestication on the brain is a reduction in size (Kruska, 2005). Dogs are the oldest and perhaps the archetypal domesticate, having split from wolves an estimated 10,000-30,000 years ago (Skoglund et al., 2011; MacHugh et al., 2017). Past research on the neural correlates of wolf-to-dog domestication has implicated PFC and the limbic system, particularly the hypothalamo-pituitary-adrenal axis. For example, gene expression in the hypothalamus is conserved between wolves and coyotes but diverges significantly in dogs (Saetre et al., 2004). Similarly, genes showing high differentiation between wolves and Chinese native dogs show high expression bias for the brain, particularly those expressed in PFC (Li et al., 2013). In an MRI study of 8 wild carnivore species and 13 domestic dogs, the allometric scaling of corpus callosum size to total brain size was constant across species, except in the rostral component, which interconnects prefrontal cortices (Spocter et al., 2018). Moreover, the enlargement of one component of carnivore PFC, the prorean gyrus, which has extensive connections with other limbic system components (Cavada and Reinoso-Suárez, 1985; Markow-Rajkowska and Kosmal, 1987), has been implicated in the emergence of complex social behavior in canid evolution (Radinsky, 1969).

The domestication of wolves into dogs is paralleled by the well-known and long-running Russian fox experiment (Trut et al., 2009). The experiment does not, of course, perfectly recapitulate the “naturally” occurring domestication of any species, and indeed domestication trajectories differ in important ways across species and can involve selection pressures beyond tameness (Zeder, 2012). Nonetheless, the Russian fox experiment offers a singularly well-controlled window on neurobiological shifts associated with tameness (Zeder, 2012). Since 1959, researchers at the Institute of Cytology and Genetics in Novosibirsk have been breeding conventional farm foxes on the basis of their behavioral response to human social contact. The tame strain is selected for high social approach behavior toward humans and produces dog-like behaviors, such as licking and tail wagging. In a parallel experiment, the aggressive strain is selected for the opposite behavior and reacts with defensive aggression when faced with human contact. A third conventional strain is kept on the farm but bred without deliberate selection on behavior (Statham et al., 2011). Behavioral differences between tame and aggressive foxes map to a locus on fox chromosome 12, which is homologous to the locus implicated in wolf-to-dog domestication (Kukekova et al., 2011a, 2012) and the whole-genome sequencing of tame, aggressive, and conventional foxes identifies 103 genomic regions differentiating the three populations including 46 regions, which are syntenic to canine candidate domestication regions (Kukekova et al., 2018). Differential gene expression across strains has been established in PFC, basal forebrain, hypothalamus, and anterior pituitary (Kukekova et al., 2011b; Hekman et al., 2018; Wang et al., 2018; Rosenfeld et al., 2020). Additionally, tame foxes show increased adult neurogenesis in the hippocampus in comparison to conventional foxes (Huang et al., 2015). However, the brain-wide neuroanatomical consequences of selection on behavior in the fox model are as yet unknown. In this study, we addressed this question using high-resolution, T2-weighted, whole-brain ex vivo neuroimaging in 10 tame, 10 aggressive, and 10 conventional foxes. These were the same individuals used in previous transcriptomic studies (Kukekova et al., 2011a, 2012), and their behavior was tested and analyzed as previously described (Kukekova et al., 2008). We conducted three types of analyses: a comparison of overall gray matter, white matter, and total brain volumes; a comparison of differences in regional gray matter volumes across strains; and a comparison of strain-wise differences in anatomic covariation across regions, which can identify morphologically coevolving structural networks and link them to individual variation in behavior.

Materials and Methods

Brain specimens

In the current study, we examined the brains of 10 tame, 10 conventional, and 10 aggressive foxes. All foxes were housed with littermates and mother from birth until weaning at 6 weeks, then housed with littermates for a further 4 weeks, and then singly housed in a shed for young foxes with members of the same strain until ∼7 months of age. After that, all 30 foxes used in this experiment were singly housed in the same shed for adult foxes. At all stages of rearing and maintenance, fox interaction with humans was limited to feeding, veterinary care, and testing behavior in home cages at 5.5-6 months of age. All foxes were male, sexually naive, and ∼1.5 years old at the time of brain extraction in August of 2010. Right hemispheres were preserved for gene expression studies; we report analyses in left hemispheres here. Left hemispheres were formalin-fixed and were maintained in formalin solution until scanning in 2016.

Neuroimaging data acquisition

For imaging, brains were placed in a waterproof plastic container, which was packed with polyethylene beads for stabilization. The container was then pumped full of Fluorinert FC-770 (3M). Fluorinert is a fluorocarbon; it is analogous to a hydrocarbon, but with fluorine taking the place of hydrogen. It thus produces no signal in (typical) MRI, which is tuned to the resonant frequency of hydrogen nuclei. Images were acquired on a 9.4 T/20 cm horizontal bore Bruker magnet, interfaced to an Avance console, with Paravision 5.1 software (Bruker). A 7.2-cm-diameter volume radio frequency coil was used for transmission and reception with a RARE T2 sequence (2 averages, 13 ms TE, 2500 ms TR, rare factor 8). Image resolution was 300 µm3 with a matrix size of 256 × 100 × 88.

Image analysis

Image preprocessing was accomplished using the FSL software package (Smith et al., 2004; Woolrich et al., 2009; Jenkinson et al., 2012). Images underwent bias correction and segmentation into white matter and gray matter using FAST (Zhang et al., 2001). In order to provide a common spatial framework for morphometric analysis, an unbiased nonlinear template was built from the 10 conventional foxes' T2-weighted images using the ANTS software package (Avants et al., 2009). This template represents the group average morphology across the conventional fox brain specimens. All subject's T2-weighted images were nonlinearly aligned to this template. We computed the Jacobian determinant image of each deformation field; this represents a spatial map of where and how much each individual subject's scan had to deform to come into alignment with the template. The Jacobian determinant images were then masked with the gray matter segmentation images to produce representations of each subject's gray matter deviation from the template, and were then smoothed with a 0.6 mm (2 voxel) Gaussian kernel. Individual foxes' gray-matter masked, smoothed Jacobian determinant images became the input for each of two independent, complementary statistical morphometric analyses. The average of all foxes' masked Jacobian images is shown in Figure 1B.

Voxel-based morphometry (VBM) (Ashburner and Friston, 2000) is an inherently hypothesis-based approach that performs GLM at each voxel in the image to determine whether morphology is significantly related to explanatory variables (e.g., group). VBM analysis was accomplished using FSL's randomize tool for voxel-wise Monte Carlo permutation testing of GLMs, which permutes explanatory variables across cases to build up a null distribution, and then tests whether observed associations to explanatory variables significantly differ from this random null distribution (Winkler et al., 2014). On the other hand, source-based morphometry (SBM) (Xu et al., 2009) is a data-driven, model-free approach that identifies patterns of significant morphologic correlation across subjects. In other words, it determines which regions of the brain significantly covary with each other across the entire dataset, while remaining agnostic to putative differences across subjects (e.g., group differences). Post hoc tests can then be used to determine whether these networks show significant associations with variables of interest. This was accomplished using the GIFT SBM toolbox for MATLAB (Xu et al., 2009; http://mialab.mrn.org/software/gift/index.html; for an in-depth discussion of the method and its applications, see Gupta et al., 2019). Multiple linear regression was then used to assess the relationship between factor loadings for covarying gray matter networks and factor loadings for a principal components analysis of fox behavioral traits described below. Past research has used SBM mainly in the context of identifying neuroanatomical differences between patient and control populations (e.g., schizophrenia) (Xu et al., 2009; Turner et al., 2012; Rodrigue et al., 2020), but several studies have also used SBM in animal studies focused on evolved traits: for example, chimpanzee tool use (Hopkins et al., 2019) and eye gaze (Hopkins et al., 2020); and breed-specialized skills in domestic dogs (Hecht et al., 2019). Covariance in morphology across regions is thought to reflect both shared genetic factors and similar patterns of development and functional specialization; these structurally covarying networks appear to be similar to functionally covarying resting state networks (Alexander-Bloch et al., 2013; Evans, 2013; Gupta et al., 2019). Although SBM identifies networks of structurally covarying regions, it does not assess white matter connections linking those regions.

Assignment of fox behavioral phenotypes

Fox behavior was tested at 5.5-6 months of age in the standard test described by Kukekova et al. (2008). Fox behavior during the test was scored from the video records for 98 recordable observations (Table 1). The matrix, including scores for 1003 foxes from Kukekova et al. (2011a) and scores for 30 foxes used in this study, were subjected to principal component analysis (PCA) in R using function: prcomp. The distribution of PC1-PC3 values for 30 foxes whose brains are analyzed are shown in Figure 4.

View this table:
  • View inline
  • View popup
Table 1.

Trait loadings for PCA of behavior

Results

The T2-weighted conventional fox brain template and labels for some anatomic regions are shown in Figure 1A. A map of variation in brain anatomy across all strains is shown in Figure 1B; a graph of gray matter, white matter, and total brain volumes for each strain are shown in Figure 1C. One-way ANOVAs revealed a strain-wise difference in total gray matter volume (F(2,27) = 6.855, p = 0.004). Post hoc Tukey's HSD tests indicated that both tame and aggressive foxes had significantly higher gray matter volume than conventional foxes (tame > conventional: t(18) = 3.370, p = 0.003; aggressive > conventional: t(18) = 3.556, p = 0.002), but did not differ significantly from each other (t(18) = 0.010, p = 0.992). Strain-wise differences in white matter volume and total brain volume did not reach significance (white matter: F(2,27) = 2.186, p = 0.132; total brain volume: F(2,27) = 1.504, p = 0.240).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Brain volume shifts in the Russian farm-fox experiment. A, Group-average template with anatomic labels and approximate cortical functional localizations from the dog (Evans and de Lahunta, 2013). olf, Olfactory; pfc, prefrontal; pmc, premotor; mot, motor; ss, somatosensory; aud, auditory; vis, visual. B, Variation in regional morphology across the dataset. Warmer colors represent more variation. Lateral views of 3D surface renderings and coronal cross-sections are shown. C, Strain-wise differences in gray matter, white matter, and total brain volume. Error bars indicate SEM.

Voxel-based morphometry identified a number of differences between strains (Fig. 2). Carnivore cortical mapping studies have been most extensive in cats and ferrets, but relatively less work has been done in dogs, which are more closely related to foxes. Thus, putative functions of anatomic regions are based on known dog anatomy here but should be considered tentative. Relative to aggressive foxes, tame foxes show expansion in portions of the sylvian gyrus and ectosylvian sulcus (temporal regions that include auditory cortex and other regions, potentially higher-order visual or multisensory cortex; Fig. 2A) (Kosmal, 2000). In contrast, compared with tame foxes, aggressive foxes show expansion in portions of the rostral composite gyrus and presylvian sulcus (potentially somatosensory and/or premotor-prefrontal transition cortex) (Kosmal et al., 1984) and in the ectosylvian and sylvian gyri and sylvian sulcus (auditory and association cortex in dogs; Fig. 2B) (Kosmal, 2000). Surprisingly, relative to the conventional strain, both tame and aggressive foxes show expansion in similar regions, including portions of the prorean, orbital, frontal, precruciate, and rostral composite gyri (prefrontal, premotor, and motor cortex in dogs) (Kosmal et al., 1984), amygdala, hippocampus, and cerebellum (Fig. 2C,D; overlap shown in Fig. 2E). No voxels showed reduced volume in tame or aggressive foxes relative to conventional foxes. Anatomical locations, maximum t statistics, number of voxels, and volume for each cluster are shown in Table 2.

View this table:
  • View inline
  • View popup
Table 2.

Clusters resulting from VBM analysisa

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Differences in regional gray matter volume between strains. A, Tame > aggressive. B, Aggressive > tame. C, Tame > conventional. D, Aggressive > conventional. E, Tame and aggressive both > conventional. Color code: red represents frontal cortex; green represents amygdala; magenta represents hippocampus; blue and cyan represents cerebellum.

Notably, the hypothalamus was one of the regions with the highest volumetric variation across the entire dataset (Fig. 1B), but direct group-wise comparisons did not identify significant volumetric differences between strains (Fig. 2). Past studies on the farm-fox experiment have implicated the hypothalamo-pituitary-adrenal axis generally (for review, see Trut et al., 2009) and gene expression in the hypothalamus specifically (Rosenfeld et al., 2020). Neuroanatomical consequences of rapid selection on behavior can be visible not only in changes to relative volume of brain regions, but also in the degree of morphologic correlation across brain regions, as we recently documented in domestic dog breeds (Hecht et al., 2019). To probe this possibility, we used source-based morphometry, a model-free, ICA-based approach (Xu et al., 2009), to identify structurally covarying multiregion networks across the entire dataset (but white matter connectivity is not assessed with this method). This identified four significantly covarying brain networks (Fig. 3; Table 3).

View this table:
  • View inline
  • View popup
Table 3.

Clusters resulting from SBM analysis

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Regionally covarying structural brain networks and relationship to strain membership and individual behavior scores. Red-yellow and blue-green components of each network are anticorrelated; only significant voxels are shown for each network. Bar graphs represent factor loadings for each strain; positive loadings indicate larger volumes in red-yellow regions (represented with red dot), whereas negative loadings indicate larger volumes for blue-green regions (represented with blue dot). Partial regression plots represent significant relationships between factor loadings for Networks 2-4 and behavior measurements PC1-PC3 in each fox. A, Network 1 factor loadings did not differentiate between strains and were not significantly related to behavior scores. B, Network 2 factor loadings differentiated between tame and aggressive strains and were significantly related to PC1 behavior scores. C, Network 3 factor loadings differentiated the tame and aggressive strains from the conventional strain and were significantly related to PC2 behavior scores. D, Network 4 factor loadings did not differentiate between strains and were marginally related to PC3 behavior scores.

We then investigated the extent to which these networks were related to strain membership and to individual behavior scores. Behavior scores were subjected to a PCA analysis, resulting in three components, each of which explains 29.0%, 8.0%, and 6.6% of variance, respectively. We then used multiple linear regression to probe the relationship between factor loadings for these three behavior components (Table 1; Fig. 4) and our four morphometry components (Table 3).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Population distributions for the first three PCs (PC1-PC3) of fox behavior. Aggr, “Aggressive” population; Conv, “conventional farm-bred population”; Tame, “tame” population. Black dots represent individual data points. Horizontal bars within each box represent the population median. The whiskers represent the range of data up to 1.5 times the interquartile range.

Network 1 contained clusters in the thalamus, caudate, NAc, cerebellum, and other regions (Fig. 3A). Factor loadings were in opposite directions for the tame and aggressive strains but were centered near zero with wide variance in the conventional strain; however, strain-wise differences did not reach significance (F(2,27) = 0.475, p = 0.627). Multiple linear regression indicated that factor loadings for Network 1 accounted for 11.6% of the variance in behavior PCA scores, but the overall regression model did not reach significance (F(3,26) = 1.137, p = 0.352), nor did any of the partial correlations between brain and behavior factor loadings.

Network 2 was mainly comprised of one cluster that covered most of the hypothalamus, plus several other clusters scattered throughout the cerebellum. Factor loadings strongly differentiated the tame and aggressive strains from each other and were centered near zero for the conventional strain. Strain-wise differences were significant (F(2,27) = 4.495, p = 0.021; Fig. 3B). Post hoc tests using Tukey correction for multiple comparisons confirmed that factor loadings significantly differentiated the tame and aggressive strains from each other (p = 0.016). Loading coefficients for Network 2 accounted for 31.0% of the variance in behavior scores. The overall regression model reached significance (F(3,26) = 3.908, p = 0.020). In examining partial correlations with individual behavior categories, we applied the Benjamini-Hochberg method for FDR correction. This revealed a significant partial correlation with PC1 behavior scores (β = 0.166, t = 2.964, p =0.006, which exceeded the critical value of 0.017). Generally, behavioral traits with the highest positive loadings for PC1 describe tame behavior and proximity to a human approacher, whereas traits with the most negative loadings describe aggressive behavior and greater distance from a human approacher (Table 1; Fig. 4). Examination of the partial correlation scatter plot reveals that tame foxes cluster together with higher scores for both behavior and neuroanatomy, whereas aggressive and conventional foxes cluster together on the lower end of both axes (Fig. 33).

Network 3 contained two smaller, discrete clusters in the hypothalamus, plus a large cluster covering much of prefrontal and premotor cortex, including portions of the prorean, orbital, frontal, precruciate, and rostral composite gyri. For this network, factor loadings were similar for the tame and aggressive strains, and both were strongly opposite to the conventional foxes; the strain-wise difference was significant (F(2,27) = 14.795, p < 0.001; Fig. 3C). Post hoc Tukey tests revealed that this network significantly differentiated both the tame strain from the conventional strain (p = 0.001) and the aggressive strain from the conventional strain (p = 0.000). Regression analyses revealed that behavior scores accounted for 35.3% of the variance in morphometry loading coefficients. The overall regression model reached significance (F(3,26) = 4.720, p = 0.009). Examination of partial correlations using Benjamini-Hochberg FDR correction revealed that there was a significant partial correlation with PC2 (β = 0.483, t = 3.251, p = 0.003, which exceeded the critical value of 0.017). For PC2, traits with the highest positive loadings describe active aggressive response and fox position in the front part of the cage, whereas traits with highest negative loadings describe a passive tame response: tolerance of human tactile contact and fox position in the back of the cage (Table 1; Fig. 4). Examination of the partial correlation scatter plot reveals that aggressive foxes tended to score high for both the behavior and neuroanatomical measures, whereas conventional foxes tended to score low on both, and tame foxes were tightly clustered in the intermediate range (Fig. 3C).

Network 4 contained a large cluster that spanned regions of the thalamus, most of the hypothalamus, and portions of the NAc/ventral forebrain and caudate (Fig. 3D). There were also several additional clusters located in cortex and cerebellum. Factor loadings were again opposite for the tame and aggressive strains, and centered near zero for the conventional strain; strain-wise differences did not reach significance (F(2,27) = 1.198, p = 0.312). Loading coefficients for Network 4 accounted for 14.6% of the variance in behavior scores. The overall regression model did not reach significance (F(3,26) = 1.413, p = 0.242). None of the partial correlations between morphometry and behavior scores reached significance, but there was a marginally positive relationship with PC3 (β = 0.492, t = 1.780, p = 0.087, which did not exceed the critical value of 0.017). In PC3, traits with the highest positive loadings describe neutral exploratory behavior (e.g., “ears are vertical”), whereas traits with the most negative loadings describe prosocial greeting behavior and fox position in the front part of the cage (see Table 1). Tame, aggressive, and conventional foxes showed marked overlap in both behavior and neural scores in this analysis (Fig. 3D).

Discussion

The Russian farm-fox experiment is perhaps the longest-running, best-controlled, and most well-known artificial selection study bearing on the evolution of mammalian behavior. As such, it enables a uniquely powerful window on the neural mechanisms governing behavioral adaptation. This study used high-resolution MRI to examine the brains of these foxes. We found that selection on social behavior has altered the anatomy of distributed gray matter networks, which included, among other regions, PFC, hippocampus, amygdala, caudate, NAc, cerebellum, and hypothalamus. These regions have also been implicated in past studies in these foxes (Kukekova et al., 2011b; Huang et al., 2015; Hekman et al., 2018; Wang et al., 2018; Rosenfeld et al., 2020) and in past studies on wolf-to-dog domestication (Nikulina, 1991; Saetre et al., 2004; Natt et al., 2012; Li et al., 2013; Ruan and Zhang, 2016; Spocter et al., 2018; Oshchepkov et al., 2019). Interestingly, some of these regions are also affected by domestication in rabbits (Brusini et al., 2018).

Notably, portions of fox PFC appeared in several of our results (Figs. 2E, 3A, 3C). How can the functional relevance of this region be interpreted? Although some of these prefrontal results are located on the dorsolateral surface of the brain (i.e., the prorean gyrus), this region is likely not homologous to the granular dorsolateral PFC of humans and macaque monkeys, as that region is thought to be unique to primates (Preuss, 1995; Passingham and Wise, 2012). Rather, the connectivity and cytoarchitecture of carnivore PFC (Narkiewicz and Brutkowski, 1967; Kosmal and Dabrowska, 1980; Kosmal, 1981a,b; Stepniewska and Kosmal, 1986; Markow-Rajkowska and Kosmal, 1987; Rajkowska and Kosmal, 1988) are more similar to that of the dysgranular and agranular portions of primate orbitofrontal and ventromedial PFC, regions that function to integrate external multisensory and internal visceromotor information with limbic reward and threat signals to compare the value of potential behavioral choices (Carmichael and Price, 1996; Rudebeck and Murray, 2011). Putatively, this circuitry has been altered in tame foxes to bias behavioral decisions toward the reward value of social contact with humans, and in aggressive foxes toward the opposite. This interpretation fits generally with conceptualizations of partially shared circuitry between social approach and social avoidance processing in humans and other species (Aupperle and Paulus, 2010).

However, the current study also produced some findings which were unexpected and suggest revision of existing thinking about domestication. One of these was that total gray matter volume in both the tame strain and the aggressive strain is increased compared with the conventional strain. These findings are in contrast to a number of prior studies which have reported that domestication reduces brain size in diverse species including Atlantic cod (Mayer et al., 2011), guppies (Burns et al., 2009), rainbow trout (Marchetti and Nevitt, 2003), mallard ducks (Ebinger, 1995), rats, mice, gerbils, guinea pigs, rabbits, pigs, sheep, llamas, horses, ferrets, cats, and dogs (for review, see Kruska, 2005). Brain size should be interpreted in the context of body size, given that the two covary. While body size measurements were not available for the foxes in our study, Huang et al. (2015) report a 15.9% reduction in body weight in tame compared with conventional foxes, suggesting that, if anything, differences in brain:body size ratios are more pronounced than the brain measurements alone would indicate. What might cause this increase in brain size? One possible explanation might involve neuroplastic changes resulting from differential lifetime experience. The link between naturalistic, enriched environments and larger, more complex brains has been noted for decades (van Praag et al., 2000; Lambert et al., 2019). Tame, aggressive, and conventional foxes are all housed in identical conditions and undergo identical treatment with minimally necessary human interaction. Nonetheless, tame foxes, because of their innate predisposition toward prosocial interaction with humans, may effectively experience this environment to be more naturalistic and enriched.

However, this potential interpretation does not explain why the aggressive strain also shows increased gray matter volume relative to controls. Perhaps their constant drive to avoid human contact functions as a sort of enrichment, but we instead favor an alternative explanation. Both the tame and aggressive strains have been subject to intense, sustained selection on behavior, while the conventional strain undergoes no such intentional selection. Thus, it is possible that fast evolution of behavior, at least initially, may generally proceed via increases in gray matter. Several potential mechanisms might underlie this effect. First, because of developmental linkages, selection pressure for increased size in one region may “drag along” enlargement in others to some extent (Finlay et al., 2001), although it seems apparent that the current dataset also reflects at least some degree of modular evolution (i.e., focused selection pressure on specific regions and networks) (Striedter, 2020). Second, selection pressure for expansion of a particular brain region might favor genomic variants leading to changes in expression of genes with potential pleiotropic effects (e.g., transcription factors). Consequently, these gene expression changes could produce expansion not only in that one region but also in others. Normally, trade-offs with metabolic and life history constraints (Isler and van Schaik, 2009) might favor further fine-tuning of the activity of these genes so that the effect occurs more specifically only in the “targeted” brain region, but because captivity reduces the pressure for optimization, this fine-tuning might fail to occur. Third, because of the complex interdependencies between brain systems, adaptive solutions that modify one existing system may stand a high likelihood of producing deleterious effects in others unless additional compensatory adaptations also occur (De Vries, 2004). In the context of extreme selection pressure, these compensatory changes might be more likely to occur via the addition of new neural material rather than via volumetrically net-zero alterations to existing circuits, which putatively would require more fine-tuned genetic changes. In the wild, there would be constant pressure against these “easy, wasteful” solutions, but in the context of the farm-fox experiment, these constraints may have been lessened or removed, as these foxes are not bred for increased meat yields, more rapid maturation, or reduced feeding costs. In general, it is important to note that strong selection for a single dimension of behavior in captivity, as occurred in the farm-fox experiment, is expected to differ in important ways from selection in the wild, which occurs on many behaviors at once. Because of this, experimentally applied selection pressure might be more likely to produce “nonoptimal” brain network changes, in contrast to in the wild, where such shifts might produce detrimental effects on other behaviors, and, importantly, on animal survival. There are multiple potential cellular-level causes for increased gray matter, including increased neuron count, increased neuron size, and/or dendritic changes, such as increased arborization or spine density (Zatorre et al., 2012; Keifer et al., 2015). Future research, including single-cell transcriptomics and histologic work, will be required to differentiate among these possibilities.

Notably, a recent skeletal morphology study did not find a significant difference in endocranial volume between tame and conventional foxes (Kistner et al., 2021). The foxes that became the progenitor population for Belyaev's study had existed in a farmed state in Russia for ∼50 years (Kochergin, 1936; Vahrameyev and Belyaev, 1948); this original Russian farmed population was itself drawn largely from Eastern Canadian foxes (Statham et al., 2011; Rando et al., 2017) bred at fur farms on Prince Edward Island (Forester and Forester, 1982). Kistner et al. (2021) compared endocranial volumes of modern Russian farm-fox skulls to skulls from wild foxes, which had been collected in Canada east of Quebec between 1884 and 1952, with 70% collected between 1894 and 1900. This revealed that endocranial volume in the modern, conventional Russian farm-raised foxes (i.e., those bred without selection on behavior) was significantly reduced compared with the archived Canadian wild skulls. However, no significant difference was observed between modern tame and conventional strains, but notably, like the current study, Kistner et al. (2021) were unable to incorporate body size measurements into their analyses. We also do not know whether the skull morphology of conventional foxes started to show differences at the time when selection of conventional foxes for tameness in Russia began (1959) or whether it became pronounced later when the conventional population was evolving in parallel with the tame population. Together with these findings, the current study hints that the brain volume reduction associated with domestication might occur relatively quickly, before the onset of intentional selection on behavior. Further research is required to establish whether this is true.

A second surprising result from this study was that both tame and aggressive foxes showed enlargement in substantially overlapping gray matter regions, including PFC, amygdala, hippocampus, and cerebellum (Fig. 2E). In addition to these volumetric effects, we also observed strain-wise differences in the degree of morphologic covariation across distributed, multiregion networks. These latter measurements revealed links with individual variation in behavior, including in brain regions that did not show volumetric differences, notably, the hypothalamus. Additionally, Network 3 consisted primarily of the hypothalamus and PFC, two regions strongly implicated in both fox and dog domestication (Nikulina, 1991; Saetre et al., 2004; Kukekova et al., 2011b; Natt et al., 2012; Li et al., 2013; Ruan and Zhang, 2016; Spocter et al., 2018; Wang et al., 2018; Oshchepkov et al., 2019; Rosenfeld et al., 2020). In this network, factor loadings did not differentiate the tame from aggressive strains; rather, the selectively bred strains together were differentiated from the conventional strain (Fig. 3C). Together, these results indicate that selection for opposite behavioral responses (docility vs aggression) can produce similar evolved changes in the brain. This has important implications for attempts to evaluate the hypothesis that humans are ourselves self-domesticated (Sanchez-Villagra and van Schaik, 2019), given that similar neuroanatomical patterns of change could now be interpreted to support either selection for increased or decreased aggression in our lineage. Notably, Trut et al. (2017) reported some elements of the “domestication syndrome” in aggressive foxes, including altered timing of annual reproductive activity, facial foreshortening, and facial “star” marking; they thus propose that these traits might result pleiotropically from selection on either tameness or aggression (i.e., as a result of selection on social behavior generally). This calls for additional research on the neural correlates of evolved differences in both docility and aggression at the cellular and genomic levels.

The Russian farm-fox experiment represents a uniquely well-controlled opportunity to study the effects of specific, sustained selection on behavior. Thus, apart from questions of domestication, an additional implication of these results concerns brain evolution on a more general level. We found that intense selection on behavior can produce gross changes in distributed brain morphology extremely rapidly (within the span of well under 100 generations). This suggests that the brains of other animals on this planet, including Homo sapiens, may have undergone similarly precipitous morphologic shifts any time steep selection on behavior was experienced.

aWithin each network, the negative components (shown as blue in Fig. 3) and positive components (shown as red in Fig. 3) are morphologically anticorrelated. Each component contains multiple clusters, indicated by Cluster index. Voxel coordinates are for maximum value in cluster.

Footnotes

  • E.E.H., T.M.P., D.A.G., and neuroimaging scan costs were supported by National Science Foundation IOS #1457291. A.V.K., sample collection, and behavioral data analysis were supported by National Institutes of Health Grant GM120782. We thank Jaekeun Park and Orion Keifer for assistance with scan acquisition; Olivia Zarella and Jeromy Dooyema for assistance with sample preparation; Yury E. Herbeck and Anastasiya V. Kharlamova for insightful discussions and help in preparation of the experiments at the experimental farm of the Institute of Cytology and Genetics (Novosibirsk, Russia); and Irina V. Pivovarova, Anastasiya V. Vladimirova, Tatyana I. Semenova, Eugene A. Martinov, and all the animal keepers at the Institute of Cytology and Genetics experimental farm for research assistance. L.N.T., the fox colony, sample collection, and behavioral data analysis were supported by the Institute of Cytology and Genetics of the Siberian Branch of the Russian Academy of Sciences grant 0259-2021-0016

  • The authors declare no competing financial interests.

  • Correspondence should be addressed to Erin E. Hecht at erin_hecht{at}fas.harvard.edu

SfN exclusive license.

References

  1. ↵
    1. Alexander-Bloch A,
    2. Giedd JN,
    3. Bullmore E
    (2013) Imaging structural co-variance between human brain regions. Nat Rev Neurosci 14:322–336. doi:10.1038/nrn3465 pmid:23531697
    OpenUrlCrossRefPubMed
  2. ↵
    1. Ashburner J,
    2. Friston KJ
    (2000) Voxel-based morphometry: the methods. Neuroimage 11:805–821. doi:10.1006/nimg.2000.0582 pmid:10860804
    OpenUrlCrossRefPubMed
  3. ↵
    1. Aupperle RL,
    2. Paulus MP
    (2010) Neural systems underlying approach and avoidance in anxiety disorders. Dialogues Clin Neurosci 12:517–531. pmid:21319496
    OpenUrlCrossRefPubMed
  4. ↵
    1. Avants B,
    2. Tustison N,
    3. Song G,
    4. Gee J
    (2009) ANTS: Advanced Open-Source Normalization Tools for Neuroanatomy. Philadelphia: Penn Image Computing and Science Laboratory.
  5. ↵
    1. Brune M
    (2007) On human self-domestication, psychiatry, and eugenics. Philos Ethics Humanit Med 2:21.
    OpenUrlCrossRefPubMed
  6. ↵
    1. Brusini I,
    2. Carneiro M,
    3. Wang C,
    4. Rubin CJ,
    5. Ring H,
    6. Afonso S,
    7. Blanco-Aguiar JA,
    8. Ferrand N,
    9. Rafati N,
    10. Villafuerte R,
    11. Smedby O,
    12. Damberg P,
    13. Hallbook F,
    14. Fredrikson M,
    15. Andersson L
    (2018) Changes in brain architecture are consistent with altered fear processing in domestic rabbits. Proc Natl Acad Sci USA 115:7380–7385. doi:10.1073/pnas.1801024115 pmid:29941556
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Burns J,
    2. Saravanan A,
    3. Rodd FH
    (2009) Rearing environment affects the brain size of guppies: lab-reared guppies have smaller brains than wild-caught guppies. Ethology 115:122–133. doi:10.1111/j.1439-0310.2008.01585.x
    OpenUrlCrossRef
  8. ↵
    1. Carmichael ST,
    2. Price JL
    (1996) Connectional networks within the orbital and medial prefrontal cortex of macaque monkeys. J Comp Neurol 371:179–207. doi:10.1002/(SICI)1096-9861(19960722)371:2<179::AID-CNE1>3.0.CO;2-#
    OpenUrlCrossRefPubMed
  9. ↵
    1. Cavada C,
    2. Reinoso-Suárez F
    (1985) Topographical organization of the cortical afferent connections of the prefrontal cortex in the cat. J Comp Neurol 242:293–324. doi:10.1002/cne.902420302 pmid:2418073
    OpenUrlCrossRefPubMed
  10. ↵
    1. De Vries GJ
    (2004) Minireview: sex differences in adult and developing brains: compensation, compensation, compensation. Endocrinology 145:1063–1068. doi:10.1210/en.2003-1504 pmid:14670982
    OpenUrlCrossRefPubMed
  11. ↵
    1. Ebinger P
    (1995) Domestication and plasticity of brain organization in mallards (Anas platyrhynchos). Brain Behav Evol 45:286–300. doi:10.1159/000113557 pmid:7620877
    OpenUrlCrossRefPubMed
  12. ↵
    1. Evans AC
    (2013) Networks of anatomical covariance. Neuroimage 80:489–504. doi:10.1016/j.neuroimage.2013.05.054 pmid:23711536
    OpenUrlCrossRefPubMed
  13. ↵
    1. Evans H,
    2. de Lahunta A
    (2013) Miller's anatomy of the dog. New York: Elsevier Saunders.
  14. ↵
    1. Finlay BL,
    2. Darlington RB,
    3. Nicastro N
    (2001) Developmental structure in brain evolution. Behav Brain Sci 24:263–278; discussion 278-308. doi:10.1017/S0140525X01003958 pmid:11530543
    OpenUrlCrossRefPubMed
  15. ↵
    1. Forester JE,
    2. Forester AD
    (1982) Silver fox odyssey: history of the Canadian silver fox industry. Canada: Canadian Silver Fox Breeders Association.
  16. ↵
    1. Gupta CN,
    2. Turner JA,
    3. Calhoun VD
    (2019) Source-based morphometry: a decade of covarying structural brain patterns. Brain Struct Funct 224:3031–3044. doi:10.1007/s00429-019-01969-8 pmid:31701266
    OpenUrlCrossRefPubMed
  17. ↵
    1. Hare B
    (2017) Survival of the friendliest: homo sapiens evolved via selection for prosociality. Annu Rev Psychol 68:155–186. doi:10.1146/annurev-psych-010416-044201 pmid:27732802
    OpenUrlCrossRefPubMed
  18. ↵
    1. Hecht EE,
    2. Smaers JB,
    3. Dunn WD,
    4. Kent M,
    5. Preuss TM,
    6. Gutman DA
    (2019) Significant neuroanatomical variation among domestic dog breeds. J Neurosci 39:7748–7758. doi:10.1523/JNEUROSCI.0303-19.2019 pmid:31477568
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. Hekman JP,
    2. Johnson JL,
    3. Edwards W,
    4. Vladimirova AV,
    5. Gulevich RG,
    6. Ford AL,
    7. Kharlamova AV,
    8. Herbeck Y,
    9. Acland GM,
    10. Raetzman LT,
    11. Trut LN,
    12. Kukekova AV
    (2018) Anterior pituitary transcriptome suggests differences in ACTH release in tame and aggressive foxes. G3 (Bethesda) 8:859–873. doi:10.1534/g3.117.300508 pmid:29378821
    OpenUrlAbstract/FREE Full Text
  20. ↵
    1. Hopkins WD,
    2. Latzman RD,
    3. Mareno MC,
    4. Schapiro SJ,
    5. Gomez-Robles A,
    6. Sherwood CC
    (2019) Heritability of gray matter structural covariation and tool use skills in chimpanzees (Pan troglodytes): a source-based morphometry and quantitative genetic analysis. Cereb Cortex 29:3702–3711. doi:10.1093/cercor/bhy250 pmid:30307488
    OpenUrlCrossRefPubMed
  21. ↵
    1. Hopkins WD,
    2. Mulholland MM,
    3. Reamer LA,
    4. Mareno MC,
    5. Schapiro SJ
    (2020) The role of early social rearing, neurological, and genetic factors on individual differences in mutual eye gaze among captive chimpanzees. Sci Rep 10:7412. doi:10.1038/s41598-020-64051-y pmid:32366881
    OpenUrlCrossRefPubMed
  22. ↵
    1. Huang S,
    2. Slomianka L,
    3. Farmer AJ,
    4. Kharlamova AV,
    5. Gulevich RG,
    6. Herbeck YE,
    7. Trut LN,
    8. Wolfer DP,
    9. Amrein I
    (2015) Selection for tameness, a key behavioral trait of domestication, increases adult hippocampal neurogenesis in foxes. Hippocampus 25:963–975. doi:10.1002/hipo.22420 pmid:25616112
    OpenUrlCrossRefPubMed
  23. ↵
    1. Isler K,
    2. van Schaik CP
    (2009) The Expensive Brain: a framework for explaining evolutionary changes in brain size. J Hum Evol 57:392–400. doi:10.1016/j.jhevol.2009.04.009 pmid:19732937
    OpenUrlCrossRefPubMed
  24. ↵
    1. Jenkinson M,
    2. Beckmann CF,
    3. Behrens TE,
    4. Woolrich MW,
    5. Smith SM
    (2012) Fsl. Neuroimage 62:782–790. doi:10.1016/j.neuroimage.2011.09.015 pmid:21979382
    OpenUrlCrossRefPubMed
  25. ↵
    1. Keifer OP Jr.,
    2. Hurt RC,
    3. Gutman DA,
    4. Keilholz SD,
    5. Gourley SL,
    6. Ressler KJ
    (2015) Voxel-based morphometry predicts shifts in dendritic spine density and morphology with auditory fear conditioning. Nat Commun 6:7582. doi:10.1038/ncomms8582 pmid:26151911
    OpenUrlCrossRefPubMed
  26. ↵
    1. Kistner T,
    2. Zink KD,
    3. Worthington S,
    4. Lieberman DE
    (2021) Geometric morphometric investigation of craniofacial morphological change in domesticated silver foxes. Sci Rep 11:2582.
    OpenUrl
  27. ↵
    1. Kochergin AF
    (1936) Breeding and maintenance of silver-black foxes and racoon dogs. Novosibirsk: Western-Siberian region publisher.
  28. ↵
    1. Kosmal A
    (1981a) Subcortical connections of the prefrontal cortex in dogs: afferents to the proreal gyrus. Acta Neurobiol Exp (Wars) 41:69–85. pmid:6267904
    OpenUrlPubMed
  29. ↵
    1. Kosmal A
    (1981b) Subcortical connections of the prefrontal cortex in dogs: afferents to the medial cortex. Acta Neurobiol Exp (Wars) 41:339–356. pmid:7325014
    OpenUrlPubMed
  30. ↵
    1. Kosmal A
    (2000) Organization of connections underlying the processing of auditory information in the dog. Prog Neuropsychopharmacol Biol Psychiatry 24:825–854. doi:10.1016/S0278-5846(00)00109-3 pmid:11191716
    OpenUrlCrossRefPubMed
  31. ↵
    1. Kosmal A,
    2. Dabrowska J
    (1980) Subcortical connections of the prefrontal cortex in dogs: afferents to the orbital gyrus. Acta Neurobiol Exp (Wars) 40:593–608. pmid:6254345
    OpenUrlPubMed
  32. ↵
    1. Kosmal A,
    2. Markow G,
    3. Stepniewska I
    (1984) The presylvian cortex as a transitional prefronto-motor zone in dog. Acta Neurobiol Exp (Wars) 44:273–287. pmid:6536200
    OpenUrlPubMed
  33. ↵
    1. Kruska DC
    (2005) On the evolutionary significance of encephalization in some eutherian mammals: effects of adaptive radiation, domestication, and feralization. Brain Behav Evol 65:73–108. doi:10.1159/000082979 pmid:15627722
    OpenUrlCrossRefPubMed
  34. ↵
    1. Kukekova AV,
    2. Trut LN,
    3. Chase K,
    4. Shepeleva DV,
    5. Vladimirova AV,
    6. Kharlamova AV,
    7. Oskina IN,
    8. Stepika A,
    9. Klebanov S,
    10. Erb HN,
    11. Acland GM
    (2008) Measurement of segregating behaviors in experimental silver fox pedigrees. Behav Genet 38:185–194. doi:10.1007/s10519-007-9180-1 pmid:18030612
    OpenUrlCrossRefPubMed
  35. ↵
    1. Kukekova AV,
    2. Trut LN,
    3. Chase K,
    4. Kharlamova AV,
    5. Johnson JL,
    6. Temnykh SV,
    7. Oskina IN,
    8. Gulevich RG,
    9. Vladimirova AV,
    10. Klebanov S,
    11. Shepeleva DV,
    12. Shikhevich SG,
    13. Acland GM,
    14. Lark KG
    (2011a) Mapping loci for fox domestication: deconstruction/reconstruction of a behavioral phenotype. Behav Genet 41:593–606. doi:10.1007/s10519-010-9418-1 pmid:21153916
    OpenUrlCrossRefPubMed
  36. ↵
    1. Kukekova AV,
    2. Johnson JL,
    3. Teiling C,
    4. Li L,
    5. Oskina IN,
    6. Kharlamova AV,
    7. Gulevich RG,
    8. Padte R,
    9. Dubreuil MM,
    10. Vladimirova AV,
    11. Shepeleva DV,
    12. Shikhevich SG,
    13. Sun Q,
    14. Ponnala L,
    15. Temnykh SV,
    16. Trut LN,
    17. Acland GM
    (2011b) Sequence comparison of prefrontal cortical brain transcriptome from a tame and an aggressive silver fox (Vulpes vulpes). BMC Genomics 12:482. doi:10.1186/1471-2164-12-482
    OpenUrlCrossRefPubMed
  37. ↵
    1. Kukekova AV,
    2. Temnykh SV,
    3. Johnson JL,
    4. Trut LN,
    5. Acland GM
    (2012) Genetics of behavior in the silver fox. Mamm Genome 23:164–177. doi:10.1007/s00335-011-9373-z pmid:22108806
    OpenUrlCrossRefPubMed
  38. ↵
    1. Kukekova AV,
    2. Johnson JL,
    3. Xiang X,
    4. Feng S,
    5. Liu S,
    6. Rando HM,
    7. Kharlamova AV,
    8. Herbeck Y,
    9. Serdyukova NA,
    10. Xiong Z,
    11. Beklemischeva V,
    12. Koepfli KP,
    13. Gulevich RG,
    14. Vladimirova AV,
    15. Hekman JP,
    16. Perelman PL,
    17. Graphodatsky AS,
    18. O'Brien SJ,
    19. Wang X,
    20. Clark AG, et al
    . (2018) Red fox genome assembly identifies genomic regions associated with tame and aggressive behaviours. Nat Ecol Evol 2:1479–1491. doi:10.1038/s41559-018-0611-6 pmid:30082739
    OpenUrlCrossRefPubMed
  39. ↵
    1. Lambert K,
    2. Eisch AJ,
    3. Galea LA,
    4. Kempermann G,
    5. Merzenich M
    (2019) Optimizing brain performance: identifying mechanisms of adaptive neurobiological plasticity. Neurosci Biobehav Rev 105:60–71. doi:10.1016/j.neubiorev.2019.06.033 pmid:31356835
    OpenUrlCrossRefPubMed
  40. ↵
    1. Li Y,
    2. Vonholdt BM,
    3. Reynolds A,
    4. Boyko AR,
    5. Wayne RK,
    6. Wu DD,
    7. Zhang YP
    (2013) Artificial selection on brain-expressed genes during the domestication of dog. Mol Biol Evol 30:1867–1876. doi:10.1093/molbev/mst088 pmid:23660689
    OpenUrlCrossRefPubMed
  41. ↵
    1. MacHugh DE,
    2. Larson G,
    3. Orlando L
    (2017) Taming the past: ancient DNA and the study of animal domestication. Annu Rev Anim Biosci 5:329–351. doi:10.1146/annurev-animal-022516-022747 pmid:27813680
    OpenUrlCrossRefPubMed
  42. ↵
    1. Marchetti MP,
    2. Nevitt GA
    (2003) Effects of hatchery rearing on brain structures of rainbow trout, Oncorhynchus mykiss. Environ Biol Fish 66:9–14. doi:10.1023/A:1023269221678
    OpenUrlCrossRef
  43. ↵
    1. Markow-Rajkowska G,
    2. Kosmal A
    (1987) Organization of cortical afferents to the frontal association cortex in dogs. Acta Neurobiol Exp (Wars) 47:137–161. pmid:3442268
    OpenUrlPubMed
  44. ↵
    1. Mayer I,
    2. Meager J,
    3. Skjæraasen JE,
    4. Rodewald P,
    5. Sverdrup G,
    6. Fernö A
    (2011) Domestication causes rapid changes in heart and brain morphology in Atlantic cod (Gadus morhua). Environ Biol Fish 92:181–186. doi:10.1007/s10641-011-9831-1
    OpenUrlCrossRef
  45. ↵
    1. Narkiewicz O,
    2. Brutkowski S
    (1967) The organization of projections from the thalamic mediodorsal nucleus to the prefrontal cortex of the dog. J Comp Neurol 129:361–374. doi:10.1002/cne.901290406 pmid:6033987
    OpenUrlCrossRefPubMed
  46. ↵
    1. Natt D,
    2. Rubin CJ,
    3. Wright D,
    4. Johnsson M,
    5. Belteky J,
    6. Andersson L,
    7. Jensen P
    (2012) Heritable genome-wide variation of gene expression and promoter methylation between wild and domesticated chickens. BMC Genomics 13:59. doi:10.1186/1471-2164-13-59 pmid:22305654
    OpenUrlCrossRefPubMed
  47. ↵
    1. Nikulina EM
    (1991) Neural control of predatory aggression in wild and domesticated animals. Neurosci Biobehav Rev 15:545–547. doi:10.1016/s0149-7634(05)80146-0 pmid:1792017
    OpenUrlCrossRefPubMed
  48. ↵
    1. Oshchepkov D,
    2. Ponomarenko M,
    3. Klimova N,
    4. Chadaeva I,
    5. Bragin A,
    6. Sharypova E,
    7. Shikhevich S,
    8. Kozhemyakina R
    (2019) A rat model of human behavior provides evidence of natural selection against underexpression of aggressiveness-related genes in humans. Front Genet 10:1267. doi:10.3389/fgene.2019.01267 pmid:31921305
    OpenUrlCrossRefPubMed
  49. ↵
    1. Passingham RE,
    2. Wise SP
    (2012) The neurobiology of the prefrontal cortex: anatomy, evolution, and the origin of insight. Oxford: Oxford UP.
  50. ↵
    1. Preuss TM
    (1995) Do rats have prefrontal cortex? The Rose-Woolsey-Akert program reconsidered. J Cogn Neurosci 7:1–24. doi:10.1162/jocn.1995.7.1.1 pmid:23961750
    OpenUrlCrossRefPubMed
  51. ↵
    1. Radinsky L
    (1969) Outlines of canid and felid brain evolution. Ann NY Acad Sci 167:277–288. doi:10.1111/j.1749-6632.1969.tb20450.x
    OpenUrlCrossRef
  52. ↵
    1. Rajkowska G,
    2. Kosmal A
    (1988) Intrinsic connections and cytoarchitectonic data of the frontal association cortex in the dog. Acta Neurobiol Exp (Wars) 48:169–192. pmid:3188998
    OpenUrlPubMed
  53. ↵
    1. Rando HM,
    2. Stutchman JT,
    3. Bastounes ER,
    4. Johnson JL,
    5. Driscoll CA,
    6. Barr CS,
    7. Trut LN,
    8. Sacks BN,
    9. Kukekova AV
    (2017) Y-chromosome markers for the red fox. J Hered 108:678–685. doi:10.1093/jhered/esx066 pmid:28821189
    OpenUrlCrossRefPubMed
  54. ↵
    1. Rodrigue AL,
    2. Alexander-Bloch AF,
    3. Knowles EE,
    4. Mathias SR,
    5. Mollon J,
    6. Koenis MM,
    7. Perrone-Bizzozero NI,
    8. Almasy L,
    9. Turner JA,
    10. Calhoun VD,
    11. Glahn DC
    (2020) Genetic contributions to multivariate data-driven brain networks constructed via source-based morphometry. Cereb Cortex 30:4899–4913. doi:10.1093/cercor/bhaa082 pmid:32318716
    OpenUrlCrossRefPubMed
  55. ↵
    1. Rosenfeld CS,
    2. Hekman JP,
    3. Johnson JL,
    4. Lyu Z,
    5. Ortega MT,
    6. Joshi T,
    7. Mao J,
    8. Vladimirova AV,
    9. Gulevich RG,
    10. Kharlamova AV,
    11. Acland GM,
    12. Hecht EE,
    13. Wang X,
    14. Clark AG,
    15. Trut LN,
    16. Behura SK,
    17. Kukekova AV
    (2020) Hypothalamic transcriptome of tame and aggressive silver foxes (Vulpes vulpes) identifies gene expression differences shared across brain regions. Genes Brain Behav 19:e12614.
    OpenUrlCrossRef
  56. ↵
    1. Ruan C,
    2. Zhang Z
    (2016) Laboratory domestication changed the expression patterns of oxytocin and vasopressin in brains of rats and mice. Anat Sci Int 91:358–370. doi:10.1007/s12565-015-0311-0 pmid:26553093
    OpenUrlCrossRefPubMed
  57. ↵
    1. Rudebeck PH,
    2. Murray EA
    (2011) Balkanizing the primate orbitofrontal cortex: distinct subregions for comparing and contrasting values. Ann NY Acad Sci 1239:1–13. doi:10.1111/j.1749-6632.2011.06267.x pmid:22145870
    OpenUrlCrossRefPubMed
  58. ↵
    1. Saetre P,
    2. Lindberg J,
    3. Leonard JA,
    4. Olsson K,
    5. Pettersson U,
    6. Ellegren H,
    7. Bergstrom TF,
    8. Vila C,
    9. Jazin E
    (2004) From wild wolf to domestic dog: gene expression changes in the brain. Brain Res Mol Brain Res 126:198–206. doi:10.1016/j.molbrainres.2004.05.003 pmid:15249144
    OpenUrlCrossRefPubMed
  59. ↵
    1. Sanchez-Villagra MR,
    2. van Schaik CP
    (2019) Evaluating the self-domestication hypothesis of human evolution. Evol Anthropol 28:133–143. doi:10.1002/evan.21777 pmid:30938920
    OpenUrlCrossRefPubMed
  60. ↵
    1. Skoglund P,
    2. Gotherstrom A,
    3. Jakobsson M
    (2011) Estimation of population divergence times from non-overlapping genomic sequences: examples from dogs and wolves. Mol Biol Evol 28:1505–1517. doi:10.1093/molbev/msq342 pmid:21177316
    OpenUrlCrossRefPubMed
  61. ↵
    1. Smil V
    (2011) Harvesting the biosphere: the human impact. Popul Dev Rev 37:613–636. doi:10.1111/j.1728-4457.2011.00450.x pmid:22319767
    OpenUrlCrossRefPubMed
  62. ↵
    1. Smith SM,
    2. Jenkinson M,
    3. Woolrich MW,
    4. Beckmann CF,
    5. Behrens TE,
    6. Johansen-Berg H,
    7. Bannister PR,
    8. De Luca M,
    9. Drobnjak I,
    10. Flitney DE,
    11. Niazy RK,
    12. Saunders J,
    13. Vickers J,
    14. Zhang Y,
    15. De Stefano N,
    16. Brady JM,
    17. Matthews PM
    (2004) Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23 Suppl 1:S208–S219. doi:10.1016/j.neuroimage.2004.07.051 pmid:15501092
    OpenUrlCrossRefPubMed
  63. ↵
    1. Spocter MA,
    2. Uddin A,
    3. Ng JC,
    4. Wong E,
    5. Wang VX,
    6. Tang C,
    7. Wicinski B,
    8. Haas J,
    9. Bitterman K,
    10. Raghanti MA,
    11. Dunn R,
    12. Hof PR,
    13. Sherwood CC,
    14. Jovanovik J,
    15. Rusbridge C,
    16. Manger PR
    (2018) Scaling of the corpus callosum in wild and domestic canids: insights into the domesticated brain. J Comp Neurol 526:2341–2359. doi:10.1002/cne.24486 pmid:30025159
    OpenUrlCrossRefPubMed
  64. ↵
    1. Statham MJ,
    2. Trut LN,
    3. Sacks BN,
    4. Kharlamova AV,
    5. Oskina IN,
    6. Gulevich RG,
    7. Johnson JL,
    8. Temnykh SV,
    9. Acland GM,
    10. Kukekova AV
    (2011) On the origin of a domesticated species: identifying the parent population of Russian silver foxes (Vulpes vulpes). Biol J Linn Soc Lond 103:168–175. doi:10.1111/j.1095-8312.2011.01629.x pmid:21625363
    OpenUrlCrossRefPubMed
  65. ↵
    1. Stepniewska I,
    2. Kosmal A
    (1986) Distribution of mediodorsal thalamic nucleus afferents originating in the prefrontal association cortex of the dog. Acta Neurobiol Exp (Wars) 46:311–322. pmid:3565103
    OpenUrlPubMed
  66. ↵
    1. Striedter GF
    (2020) A history of ideas in evolutionary neuroscience. In: Evolutionary neuroscience (Kaas J, ed), pp 3–16. San Diego: Academic.
  67. ↵
    1. Trut LN,
    2. Kharlamova AV,
    3. Vladimirova AV,
    4. Herbeck YE
    (2017) On selection of foxes for enhanced aggressiveness and its correlated implications. Vestn VOGiS 21:392–401. doi:10.18699/VJ17.257
    OpenUrlCrossRef
  68. ↵
    1. Trut L,
    2. Oskina I,
    3. Kharlamova A
    (2009) Animal evolution during domestication: the domesticated fox as a model. Bioessays 31:349–360. doi:10.1002/bies.200800070 pmid:19260016
    OpenUrlCrossRefPubMed
  69. ↵
    1. Turner JA,
    2. Calhoun VD,
    3. Michael A,
    4. van Erp TG,
    5. Ehrlich S,
    6. Segall JM,
    7. Gollub RL,
    8. Csernansky J,
    9. Potkin SG,
    10. Ho BC,
    11. Bustillo J,
    12. Schulz SC
    (2012) Heritability of multivariate gray matter measures in schizophrenia. Twin Res Hum Genet 15:324–335. doi:10.1017/thg.2012.1
    OpenUrlCrossRefPubMed
  70. ↵
    1. Vahrameyev KA,
    2. Belyaev DK
    (1948) Guide for fox breeding (In Russian). Moscow: International Book.
  71. ↵
    1. van Praag H,
    2. Kempermann G,
    3. Gage FH
    (2000) Neural consequences of environmental enrichment. Nat Rev Neurosci 1:191–198. doi:10.1038/35044558 pmid:11257907
    OpenUrlCrossRefPubMed
  72. ↵
    1. Wang X,
    2. Pipes L,
    3. Trut LN,
    4. Herbeck Y,
    5. Vladimirova AV,
    6. Gulevich RG,
    7. Kharlamova AV,
    8. Johnson JL,
    9. Acland GM,
    10. Kukekova AV,
    11. Clark AG
    (2018) Genomic responses to selection for tame/aggressive behaviors in the silver fox (Vulpes vulpes). Proc Natl Acad Sci USA 115:10398–10403. doi:10.1073/pnas.1800889115 pmid:30228118
    OpenUrlAbstract/FREE Full Text
  73. ↵
    1. Winkler AM,
    2. Ridgway GR,
    3. Webster MA,
    4. Smith SM,
    5. Nichols TE
    (2014) Permutation inference for the general linear model. Neuroimage 92:381–397. doi:10.1016/j.neuroimage.2014.01.060 pmid:24530839
    OpenUrlCrossRefPubMed
  74. ↵
    1. Woolrich MW,
    2. Jbabdi S,
    3. Patenaude B,
    4. Chappell M,
    5. Makni S,
    6. Behrens T,
    7. Beckmann C,
    8. Jenkinson M,
    9. Smith SM
    (2009) Bayesian analysis of neuroimaging data in FSL. Neuroimage 45:S173–S186. doi:10.1016/j.neuroimage.2008.10.055 pmid:19059349
    OpenUrlCrossRefPubMed
  75. ↵
    1. Wrangham RW
    (2018) Two types of aggression in human evolution. Proc Natl Acad Sci USA 115:245–253. doi:10.1073/pnas.1713611115 pmid:29279379
    OpenUrlAbstract/FREE Full Text
  76. ↵
    1. Xu L,
    2. Groth KM,
    3. Pearlson G,
    4. Schretlen DJ,
    5. Calhoun VD
    (2009) Source-based morphometry: the use of independent component analysis to identify gray matter differences with application to schizophrenia. Hum Brain Mapp 30:711–724. doi:10.1002/hbm.20540 pmid:18266214
    OpenUrlCrossRefPubMed
  77. ↵
    1. Zatorre RJ,
    2. Fields RD,
    3. Johansen-Berg H
    (2012) Plasticity in gray and white: neuroimaging changes in brain structure during learning. Nat Neurosci 15:528–536. doi:10.1038/nn.3045 pmid:22426254
    OpenUrlCrossRefPubMed
  78. ↵
    1. Zeder MA
    (2012) Pathways to animal domestication. In: Biodiversity in agriculture: domestication, evolution, and sustainability, pp 227–259. Cambridge: Cambridge University Press.
  79. ↵
    1. Zhang Y,
    2. Brady M,
    3. Smith S
    (2001) Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans Med Imaging 20:45–57. doi:10.1109/42.906424
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 41 (28)
Journal of Neuroscience
Vol. 41, Issue 28
14 Jul 2021
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Neuromorphological Changes following Selection for Tameness and Aggression in the Russian Farm-Fox experiment
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Neuromorphological Changes following Selection for Tameness and Aggression in the Russian Farm-Fox experiment
Erin E. Hecht, Anna V. Kukekova, David A. Gutman, Gregory M. Acland, Todd M. Preuss, Lyudmila N. Trut
Journal of Neuroscience 14 July 2021, 41 (28) 6144-6156; DOI: 10.1523/JNEUROSCI.3114-20.2021

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Neuromorphological Changes following Selection for Tameness and Aggression in the Russian Farm-Fox experiment
Erin E. Hecht, Anna V. Kukekova, David A. Gutman, Gregory M. Acland, Todd M. Preuss, Lyudmila N. Trut
Journal of Neuroscience 14 July 2021, 41 (28) 6144-6156; DOI: 10.1523/JNEUROSCI.3114-20.2021
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • aggression
  • canids
  • domestication
  • farm-fox experiment
  • neuroimaging
  • social behavior

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Lumateperone Normalizes Pathological Levels of Acute Inflammation through Important Pathways Known to Be Involved in Mood Regulation
  • Estrous Cycle Mediates Midbrain Neuron Excitability Altering Social Behavior upon Stress
  • Chronic Ethanol Exposure Modulates Periaqueductal Gray to Extended Amygdala Dopamine Circuit
Show more Research Articles

Behavioral/Cognitive

  • Pulvinar Response Profiles and Connectivity Patterns to Object Domains
  • Subtype-Specific Roles of Ellipsoid Body Ring Neurons in Sleep Regulation in Drosophila
  • Proprioceptive and Visual Feedback Responses in Macaques Exploit Goal Redundancy
Show more Behavioral/Cognitive
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.