Abstract
Humans have bred different lineages of domestic dogs for different tasks such as hunting, herding, guarding, or companionship. These behavioral differences must be the result of underlying neural differences, but surprisingly, this topic has gone largely unexplored. The current study examined whether and how selective breeding by humans has altered the gross organization of the brain in dogs. We assessed regional volumetric variation in MRI studies of 62 male and female dogs of 33 breeds. Neuroanatomical variation is plainly visible across breeds. This variation is distributed nonrandomly across the brain. A whole-brain, data-driven independent components analysis established that specific regional subnetworks covary significantly with each other. Variation in these networks is not simply the result of variation in total brain size, total body size, or skull shape. Furthermore, the anatomy of these networks correlates significantly with different behavioral specialization(s) such as sight hunting, scent hunting, guarding, and companionship. Importantly, a phylogenetic analysis revealed that most change has occurred in the terminal branches of the dog phylogenetic tree, indicating strong, recent selection in individual breeds. Together, these results establish that brain anatomy varies significantly in dogs, likely due to human-applied selection for behavior.
SIGNIFICANCE STATEMENT Dog breeds are known to vary in cognition, temperament, and behavior, but the neural origins of this variation are unknown. In an MRI-based analysis, we found that brain anatomy covaries significantly with behavioral specializations such as sight hunting, scent hunting, guarding, and companionship. Neuroanatomical variation is not simply driven by brain size, body size, or skull shape, and is focused in specific networks of regions. Nearly all of the identified variation occurs in the terminal branches of the dog phylogenetic tree, indicating strong, recent selection in individual breeds. These results indicate that through selective breeding, humans have significantly altered the brains of different lineages of domestic dogs in different ways.
Introduction
A major goal of modern neuroscience is to understand how variation in behavior, cognition, and emotion relates to underlying neural mechanisms. A massive “natural experiment” in this arena has been right under our noses: domestic dogs. Humans have selectively bred dogs for different, specialized abilities—herding or protecting livestock, hunting by sight or smell, guarding property, or providing companionship. Significant breed differences in temperament, trainability, and social behavior are readily appreciable by the casual observer, and have also been documented quantitatively (Serpell and Hsu, 2005; Tonoike et al., 2015). Furthermore, recent genetic research indicates that this behavioral variation is highly heritable (MacLean et al., 2019).
This panoply of behavioral specializations must rely on underlying neural specializations. A small number of studies have investigated neural variation in dogs, including, for example, the effects of skull shape on brain morphology (Carreira and Ferreira, 2015; Pilegaard et al., 2017) and anatomical correlates of aggression (Jacobs et al., 2007; Våge et al., 2010). However, the neural underpinnings of behavioral differences between breeds remain largely unknown.
Most modern dog breeds were developed in an intentional, goal-driven manner relatively recently in evolutionary time; estimates for the origins of the various modern breeds vary between the past few thousand to the past few hundred years (Larson et al., 2012). This strong selection pressure suggests that brain differences between breeds may be closely tied to behavior. However, selection also occurred for outward physical appearance, including craniofacial morphology. This may have placed constraints on the internal dimensions of the skull, which in turn may have had secondary effects on brain morphology. There is substantial diversification of skull shape across dog breeds, and this has been linked to behavioral differences (Drake and Klingenberg, 2010; McGreevy et al., 2013). Alternatively, neuroanatomical variation may be explained primarily by body size rather than breed membership, with different breeds' brains representing minor, random, scaled-up or scaled-down variants of a basic species-wide pattern.
Any attempt to determine whether breeding for behavior has altered dog brains would have to be able to differentiate between these competing (and potentially interacting) hypotheses. A simple comparison of regional volumes would be insufficient for several reasons. First, a significant difference in the volume of, for example, the amygdala in pit bulls versus golden retrievers might seem intuitively meaningful, but to ascertain whether such a difference was truly the result of selection pressure on behavior, the phylogenetic structure of the dog family tree needs to be taken in to account to partition variance attributable to inheritance, and equal statistical priority needs to be given to the alternative hypotheses that observed variation in morphology. Second, and perhaps most importantly, a priori comparisons of regional gray matter volumes presuppose that experimenters can identify meaningful borders between regions. For highly conserved structures with clear anatomical boundaries, like the amygdala, this task is surmountable, but very little is known about the organization of higher-order cortical regions in dogs, and some complex behaviors that are the focus of selective breeding, like herding or interspecies communication, almost certainly rely on some of these areas. For this reason, even creating the regional outlines for a simple ROI analysis would be problematic.
Therefore, the current study took a totally data-driven, whole-brain, agnostic approach to assessing morphological variation across dog brains. Our goal was to determine whether significant nonrandom variation in brain anatomy exists across dogs and, if so, to differentiate between the competing and possibly interacting explanations for this variation.
Materials and Methods
Subjects
The dataset included T2-weighted MRI scans from 62 purebred dogs of 33 different breeds. These were grouped into 10 different breed groups as defined by American Kennel Club (AKC), which ostensibly represent groupings that were developed for similar behavioral specializations, such as herding or hunting. Table 1 lists the breed, breed group, and other data for all dogs included in the study.
Image acquisition and preprocessing
T2-weighted MRI images were acquired on a 3.0 T GE Healthcare HDx MRI unit with a GE Healthcare 5147137–2 3.0T HD T/R Quad Extremity Coil. Images were opportunistically collected at the Veterinary Teaching Hospital at the University of Georgia at Athens from dogs that were referred for neurological examination but were not found to have any neuroanatomical abnormalities. All scans were re-reviewed by a board-certified veterinary neurologist before inclusion.
The preprocessing pipeline was implemented using the NiPype workflow engine (Gorgolewski et al., 2011). Both transverse-acquired and sagittally acquired images were available for each dog. Transverse-acquired images ranged from 0.234 mm2 in-plane resolution and 2.699 mm slice distance to 0.352 mm2 in-plane resolution and 3.499 mm slice distance. Sagittally acquired images ranged from 0.273 mm2 in-plane resolution and 3.200 mm slice distance to 0.430 mm2 in-plane resolution and 3.200 mm slice distance. To maximize the use of all available anatomical information, the transverse and sagittal images were combined as follows. First, we manually performed skull-stripping on the transverse image. Next, we determined the smallest ROI that completely covered the brain from the brain mask image. The transverse image and transverse brain mask were then cropped using the computed ROI coordinates. Then, the transverse images were resampled to produce isotropic voxels in all three dimensions, the sagittal image was resliced so that it was in the same orientation as the transverse images, and a rigid registration was computed from the sagittally acquired image to the original transverse image. The region containing the brain was then cropped in the sagittal image, and we then registered the smaller cropped sagittal image to the isotropically resampled transverse brain image using a rigid registration. Finally, the cropped transverse and sagittal images were then rescaled so that the robust mean intensity of both images was 100, the images were averaged together, and then the brain mask applied to this combined image. A general diagram illustrating the overall processing pipeline is included in Fig. 1-1 and a detailed NiPype registration workflow is included in Fig. 1-2. The accompanying registration code is available at https://gist.github.com/dgutman/a0e05028fab9c6509a997f703a1c7413.
Template creation
We produced a study-specific template representing the average brain morphology across the entire group, equally unbiased toward any particular image. This was accomplished using the buildtemplateparallel.sh script in the ANTS software package (Avants et al., 2009), which nonlinearly registers each image into a common spatial framework.
Experimental design and statistical analyses
Morphological analyses.
During nonlinear registration, a warp-field is produced that represents the mapping from the original image to the target image. The Jacobian of the warp-field represents the degree of warping that had to occur in each original image to bring it into alignment with the target image. To localize significant variation in gray matter morphology, we applied a one-sample t test on the demeaned log Jacobian determinant images. This was accomplished using FSL's randomize, a tool for Monte Carlo permutation testing on general linear models (Winkler et al., 2014). This analysis permutes the sign of the log Jacobian and tests the null hypothesis that variation from the mean is random and therefore symmetrically distributed and centered around zero. The resultant t-statistic image was thresholded at p < 0.05, after multiple-comparisons correction was performed using threshold-free cluster enhancement (Smith and Nichols, 2009).
To calculate neurocephalic index, we identified maximally distant points on the left–right, rostral–caudal, and dorsal–ventral axes; neurocephalic index was computed the ratio of brain width to brain length × 100.
Cephalic index is defined as the ratio of skull width to skull length × 100. For many scans in our database, the exterior of the skull was not visible, but a large database of skull measurements is publicly available (Stone et al., 2016). We computed male and female average cephalic indices separately for each breed and used these sex-specific, breed-average measures in our analyses.
To identify regional covariation in gray matter morphology, we used GIFT, a software package for MATLAB (Calhoun et al., 2001). GIFT's toolbox for source-based morphometry (SBM) (Xu et al., 2009) is a multivariate alternative to voxel-based morphometry (VBM). It uses independent components analysis to identify spatially distinct, distributed networks of regions that covary across individuals, and computes their statistical relationship to other categorical or continuous variables. T2-weighted images underwent bias field correction using ANTS's Atropos N4 tool (Avants et al., 2011) and segmentation into gray matter, white matter, and CSF using FSL's FAST tool (Zhang et al., 2001). Gray matter segmentations were warped to the study-specific template and modulated by their log Jacobian determinants to produce per-subject maps of the degree of morphological divergence from the study-specific group-average template. In other words, the input to SBM consisted of gray matter maps for each subject, where intensity at each voxel corresponded to the degree of deformation required to come into alignment with the template (i.e., the demeaned log Jacobians). The number of sources was estimated using Akaike's information criterion (AIC) (Akaike, 1974); the application of AIC in SBM is described in Xu et al. (2009). This procedure identified six components, each of which were thresholded at Z scores >1.96 or below −1.96. Multiple regression and ANOVA analyses were then used to compute the relationship of each component to AKC-defined breed groups, with the statistical threshold set at p < 0.05 after multiple comparisons correction.
Phylogenetic statistics.
Because comparative data may be nonindependent due to shared phylogenetic history, the assumptions of standard statistical methods may be violated (Harvey and Pagel, 1991). We therefore used phylogenetic comparative methods that account for phylogenetic nonindependence by including expected phylogenetic variance-covariance among species into the error term of generalized least-squares [phylogenetic generalized least squares (pGLS)] linear models (Rohlf, 2001). When quantifying linear models we additionally included a lambda parameter to account for phylogenetic signal (Pagel, 1997). To test for differences in statistical fit among linear models that include different parameters (e.g., the inclusion of grouping variables to test for differences among breed groups), we used least-squares phylogenetic analysis of covariance (pANCOVA) (Smaers and Rohlf, 2016; Smaers and Mongle, 2018). It should be noted that “phylogenetic” approaches such as pGLS and pANCOVA are interpreted in the same way as standard least-squares approaches. The only difference between standard and phylogenetic least-squares approaches is that the phylogenetic approaches weight data points according to phylogenetic relatedness (Rohlf, 2001).
We further investigated the relationship between morphological components and the phylogenetic tree by estimating the amount of change that occurs on each lineage using a multiple variance Brownian motion approach (Smaers et al., 2016; Smaers and Mongle, 2018). This approach estimates phenotypic change along individual lineages of a tree and has been shown to provide more accurate estimates than traditional ancestral estimation methods (Smaers and Mongle, 2017).
Last, we use multiregime Ornstein-Uhlenbeck (OU) approaches to estimate phylogenetic shifts in mean value directly from the data. This approach has become a standard approach in comparative biology to model trait change across a phylogeny. Specifically, this approach quantifies the evolution of a continuous trait X as dX(t) = α[θ − X(t)]dt + σdB(t) where σ captures the stochastic evolution of Brownian motion, α determines the rate of adaptive evolution toward an optimum trait value θ (90). This standard OU model has been modified into multiple-regime OU models allowing optima to vary across the phylogeny (Butler and King, 2004). Such multiregime OU models allow modeling trait evolution toward different “regimes” that each display a different mean trait value. In other words, these approaches allow estimating directly from the data where in a phylogeny a shift in mean value of a trait has occurred. To overcome inherent difficulties with optimizing OU parameters (Ho and Ané, 2014), several algorithmic improvements have been proposed. Here, we use the approach proposed by Khabbazian et al. (2016).
Results
Neuromorphological variation is plainly visible across breeds. Midline sagittal images from the raw, native-space scans of selected dogs are shown in Figure 1A. To provide a common spatial reference for measuring this variation, we created an unbiased, diffeomorphic template using the ANTS software package (Avants et al., 2009). This template represents the average brain for the entire dataset and is shown in Figure 1B.
Figure 1-1
Figure 1-2
To visualize morphological variation in a more standardized manner, we nonlinearly warped the template to each dog's native-space image. This allowed us to examine breed variation in brain morphology and size with invariant contrast and resolution. We also additionally rescaled these images to have constant rostral-caudal lengths. This allowed us to more clearly visualize variation in morphology independent from variation in size. Both sets of scaled template images are shown in Figure 1A.
To carry out quantitative assessments of regional variation in gray matter morphology, we used the Jacobian determinants of the native-space-to-template spatial deformation fields to produce a variation intensity map. These fields represent a map of where and how much each dog's scan had to adjust to become aligned to the group-average template. The SD of these maps thus indexes the extent to which brain anatomy varies across individuals and is shown in Figure 1C.
To determine whether this variation was randomly distributed across the brain or focused in specific areas, we applied Monte Carlo permutation testing on the demeaned Jacobian determinant images. Importantly, this revealed that a large proportion of the brain shows significant gray matter morphological variation across subjects, as illustrated in Figure 1D.
Given these results, we next sought to determine what accounts for this variation by probing the extent to which it is related to body size, head shape, and/or breed group membership.
Figure 2A shows the relationship between brain volume and body mass. The scaling coefficient of this relationship [pGLS; b = 0.231, 95% confidence interval (CI) = 0.26–0.36] is significantly lower than that observed across most mammals (∼0.67), indicating the occurrence of more variation in body size relative to variation in brain size than would be expected. Importantly, using the tree structure from a recent large-scale genomic analysis (Parker et al., 2017), we were able to determine that the phylogenetic signal of the brain-body allometry is negative; that is, that variation present at the tree's terminal branches is not predicted by the deeper structure of the tree. If grade shifts in the brain–body allometry exist, then these would putatively show differences among different breeds. We tested this hypothesis by estimating putative grade shifts in the brain to body allometry directly from the data using an OU modeling approach (Khabbazian et al., 2016). This analysis revealed no grade shifts, thereby indicating that a one-grade allometry is the best explanation of the bivariate brain-to-body relationship.
In mammals, head shape is commonly measured using cephalic index (also known as skull index), calculated as maximum head width divided by maximum head length. We were interested in the possibility that human-driven selection on external craniofacial morphology may have had on the internal dimensions of the skull. To assess this, we computed an analogous neurocephalic index for each dog (maximum internal cranial cavity length divided by maximum internal cranial cavity width). Figure 2B shows the relationship between neurocephalic and cephalic index. Cephalic index is a significant predictor of neurocephalic index (pGLS: b = 0.37, t = 3.70, p < 0.01). Also, here we questioned whether grade shifts in this allometry exist, putatively showing differences among breeds. This analysis revealed that the neurocephalic–cephalic allometry was thus best explained by a two-grade model (F = 31.19, p < 0.001). The breeds on the higher grade, with a greater neurocephalic index for a given cephalic index, were as follows: Basset hound, beagle, German short-haired pointer, dachshund, cavalier King Charles spaniel, springer spaniel, west highland white terrier, silky terrier, bichon frise, and maltese. Importantly, this grade difference in the neurocephalic to cephalic index aligns with a significant difference in body size (pANOVA: F = 9.73, p < 0.01; average body size 11 kg vs 23 kg in other breeds). Smaller-bodied dogs hereby have a higher neurocephalic index (more spherical brains) for a given cephalic index (external head shape).
If variation in dog brain anatomy is unrelated to behavior, then variation should be randomly distributed across regions. Alternatively, if this variation represents heritable adaptations for behavior, then significant covariance should exist in separable, independent subnetworks of regions. To assess this, we performed source-based morphometry, a multivariate alternative to voxel-based morphometry which makes use of independent components analysis. This was accomplished using the GIFT software package (Xu et al., 2009). Results revealed six networks where regional volume covaried significantly across individuals. Figure 3 shows these networks, along with factor loadings for each breed group. Major anatomical constituents of each network are labeled. Additional research is needed to definitively link the function of each network to its adaptive role in response to behavior selection. However, we note putative roles that may serve as initial hypotheses for future research.
Network 1 includes the nucleus accumbens, dorsal and ventral caudate, cingulate gyrus, olfactory peduncle, and gyrus rectus (medial prefrontal cortex). These regions are part of or connected to the mesolimbic reward system, a network implicated in reward signaling related to reinforcement learning, incentive salience, and motivation broadly across species (Alcaro et al., 2007; O'Connell and Hofmann, 2011); in dogs, the caudate nucleus activates for both food reward and human social reward (Cook et al., 2016). Tentatively, this network might be relevant for social bonding to humans, training, and skill learning.
Network 2 involves brain regions involved in olfaction and gustation, including the piriform lobe, which contains olfactory cortex, and the insula and pseudosylvian sylcus, where the cortical representation of taste is located (Evans and de Lahunta, 2013). This component also involves regions of medial frontal cortex, which is involved in downstream or higher-order processing of chemosensation and shows activation in response to olfactory stimulation in awake but not sedated dogs (Jia et al., 2014). We propose that this network might support volitional (as opposed to instinctive) responses to olfactory and gustatory stimuli.
Network 3 includes a distributed network of subcortical regions that are involved movement, eye movement, vision, and spatial navigation, including the lateral geniculate nucleus, pulvinar, hippocampus, cerebellum, oculomotor nucleus, interpeduncular nucleus, ventral tegmental area, and substantia nigra. It also involves cortical regions, including the medial part of the frontal gyrus (supplementary motor area) and the lateral gyrus (visual cortex). Tentatively, this network may reflect a circuit involved in moving through the physical environment.
Network 4 involves higher-order cortical regions that may be involved in social action and interaction. The precruciate and prorean gyri house premotor and prefrontal cortex, respectively, while the gyrus rectus is part of medial prefrontal cortex. The expansion of frontal cortex has been linked to increased sociality in extant hyena species (Holekamp et al., 2007) and, notably, the prorean gyrus has been linked to the emergence of pack structure in canid evolution (Radinsky, 1969). The sylvian, ectosylvian, and suprasylvian gyri represent regions of lateral sensory cortex situated between gustatory, auditory, and somatosensory cortex (Evans and de Lahunta, 2013) and likely contain higher-order association areas related to sensation and perception. In domestic dog fMRI studies, multisensory activation in these regions has been observed during the presentation of dog and human faces and vocalizations (Cuaya et al., 2016; Andics et al., 2017; Thompkins et al., 2018).
Network 5 includes limbic regions that have a well established role in fear, stress, and anxiety, including the hypothalamus, amygdala, and hippocampus and adjacent dentate gyrus (for review, see Tovote et al., 2015). These regions are involved in the HPA axis, which regulates behavioral and endocrine responses to environmental stressors and threats. Some of these regions are also involved in other affective and instinctual processes, including mating, memory, and aggression (O'Connell and Hofmann, 2011).
Network 6 includes early sensory processing regions for olfaction and vision, including the olfactory peduncle and part of the lateral gyrus, which is the location of primary visual cortex (Evans and de Lahunta, 2013).
Next, we investigated the relationship between these components, total brain size, and skull morphology. A significant relationship with total brain volume was present for all but component 6, where it was marginal but did not meet significance (component 1: t = 3.663, p = 0.001; component 2: t = −2.608, p = 0.014; component 3: t = 6.219, p < 0.001; component 4: t = −6.325, p < 0.001; component 5: t = 3.938, p < 0.001; component 6: t = 1.845, p = 0.076). Components 3, 4, and 6 showed significant relationships with cephalic index, whereas component 1 was marginal (component 1: t = −1.945, p = 0.064; component 3: t = −2.165, p = 0.041; component 4: t = 2.411, p = 0.024; component 6: t = −2.171, p = 0.041; pGLS). Components 1, 3, 4, and 6 showed significant relationships with neurocephalic index (component 1: t = −2.258, p = 0.032; component 3: t = −3.823, p = 0.001; component 4: t = 7.066; p < 0.001; component 6: t = −2.890, p = 0.007, pGLS).
We also investigated the relationship between these covarying morphological components and the phylogenetic tree. If variation in brain organization mainly reflects the deep ancestry of the tree, with little relationship to recent behavioral specializations, then brain morphometry should be highly statistically dependent on phylogenetic structure (i.e., high phylogenetic signal). Conversely, if brain organization is strongly tied to selective breeding for behavioral traits, then morphological traits should be divorced from the structure of the tree (i.e., low phylogenetic signal). We observed the latter (Fig. 4). The majority of changes that occur in these components take place on the terminal branches of the phylogenetic tree.
Finally, we investigated whether these regionally covarying morphological networks were related to behavior. The AKC groups individual breeds into breed groups, but these breed groups change periodically and some groups contain breeds with disparate behavioral functions: for example, the nonsporting group includes both poodles and Shar-Peis. Therefore, rather using AKC breed groups, we identified each individual breed's ostensible behavioral specialization(s) as noted on the AKC website (www.akc.org). These were entered into in a multiple regression analysis using the GIFT Source Based Morphometry toolbox. Each of the six components showed significant correlation with at least one behavioral specialization (Fig. 5). The behavioral specialization associated with the most components (four of six) was explicit companionship, and the component associated with the most behavioral specializations (six of 10) was component 4, which involves regions involved in social action and interaction. Specific associations between associated brain networks and behavioral specializations are also apparent. For example, component 3, which involves regions involved in movement, eye movement, and spatial navigation, showed a significant correlation with sight hunting, whereas Network 2, which involves regions involved in olfaction and gustation, showed a significant correlation with scent hunting.
Discussion
The current study took a comprehensive, data-driven, agnostic approach to investigating neuroanatomical variation in domestic dogs. We first questioned whether significant variation in dog brain morphology even exists. The answer is a clear “yes”: differences in gross brain anatomy are readily appreciable (Fig. 1A). This observation was further confirmed by a whole-brain, multiple-comparison-corrected, voxelwise statistical analysis (Fig. 1C,D). Having established this basic finding, we then went on to probe the relationship between multiple, potentially interacting factors that might be linked to this variation: the total size of the body or brain, the external and internal morphology of the skull, the structure of the dog phylogenetic tree, and the organization of internal brain networks.
Dogs show intraspecific variation in morphology to a degree rarely seen in nature. There is a 100-fold difference between the body mass of a Chihuahua (∼1 kg) and the body mass of a Great Dane (∼100 kg) (Sutter et al., 2008). However, we found that dog brain sizes do not scale commensurately to dog body sizes, as indicated by a relatively low scaling coefficient for the relationship between brain size and body mass. To appreciate this effect, consider the adjacent dachshund and golden retriever images in Figure 1A: the dachshund's brain takes up most of the available endocranial space, whereas the golden retriever shows noticeably larger sinuses. A phylogenetic analysis revealed that changes in relative brain size are not predicated by relatedness and are more likely the result of selection on specific terminal branches of the phylogenetic tree (i.e., individual breeds).
In comparative animal cognition research, total brain size is often used as a gross index of cognitive capacity. Several previous studies have investigated the relationship between dog body size and cognition or behavior, with apparently contradictory results (Helton and Helton, 2010; Stone et al., 2016; cf. Broadway et al., 2017). Additionally, a study that used a single scaling metric across breeds found that larger-brained (i.e., larger-bodied) dogs performed better on tests of executive function (Horschler et al., 2019). We found that larger dogs do tend to have larger brains, but that the brain to body allometry across breeds is low, indicating high variability in brain to body ratio across breeds (Fig. 2A). Furthermore, we found that a substantial amount of variation in internal dog brain morphology is related to total brain size, suggesting that evolutionary increases or decreases in relative brain volume may be driven by changes in specific groups of regions. Moreover, we found that these networks differed across breed groups. Therefore, shifts in relative brain size may be related to expansion or contraction of specific networks, potentially leading to the presence or absence of correlations between body size and behavior depending on the specific breeds or behaviors being studied.
We also found that selection for smaller body size has significantly influenced the internal morphology of the cranial cavity. For a given cephalic index, or exterior skull shape, smaller-bodied dogs have more spherical brains (Fig. 2B). This is consistent with a previous analysis linking foreshortening of the skull to ventral pitching of the brain and olfactory bulb, resulting in a more spherical brain (Roberts et al., 2010). We assessed the extent to which internal and exterior skull morphology were related to the covarying morphometric networks we identified. More networks showed a significant relationship with neurocephalic index than with cephalic index, suggesting that variation in brain morphology appears to be more tied to the internal morphology of the cranial cavity than to external craniofacial morphology, which is perhaps not surprising. Our results indicate that skull morphology is linked to the underlying anatomy of specific, different networks of brain regions; it is possible that this could underlie the reported associations between behavior and head shape (Gácsi et al., 2009; Helton, 2009; McGreevy et al., 2013). Not all networks showed a significant relationship with either cephalic index or neurocephalic index, indicating that variation in dog brain morphology is partially but not totally dependent on variation in skull morphology. Importantly, we cannot say from the current analyses whether variation in skull morphology drives variation in brain morphology, the reverse, or both.
In addition to these analyses of the gross external shape and size of the brain and skull, we also investigated internal brain organization. This was accomplished using source-based morphometry to identify maximally independent networks that explain the variation present in the dataset. We identified six such networks (Fig. 3). In the case of circuitry that is highly conserved across species, such as circuitry for reward and motivation or fear and anxiety, it is a safe bet that research on other species is a good indicator of the functional role of these systems in dogs. This cannot be assumed to be the case for circuits that involve higher-order cortical association areas. Particularly in the case of our network 4, it may be tempting to jump to conclusions about parallels with human cortical regions that are located in approximately the same location and are involved in similar tasks; for example, the fusiform face area, Wernicke's area, or the mirror system. However, it is important to remember that primates and carnivores diverged further back in time than primates and rodents: humans are more closely related to mice than to dogs. Our last common ancestor with dogs likely had a fairly smooth, simple brain (Kaas, 2011), and higher-order cortical association areas, along with whatever complex perceptual and cognitive abilities they support, have evolved independently in dogs and humans. Therefore, we stress that the functional roles of these networks, and their relationship to selection on behavior in specific breeds, should at this point still be considered an open question.
Having identified these six networks, we then investigated their relationship to the dog phylogenetic tree. We found that the majority of changes that occur in these components take place in the tree's terminal branches (i.e., individual breeds). This suggests that brain evolution in domestic dog breeds follows a “late burst model,” with directional changes in brain organization being primarily lineage specific. We also assessed whether these networks were related to selective breeding, as evidenced by the ostensible behavioral specialization(s) of each breed as noted by the AKC. In all six of the regionally covarying networks that we found, significant correlations were found with at least one behavioral specialization. Associations between brain networks and related behavioral specializations are apparent. For example, network 2, which involves regions that support higher-order olfactory processing, shows a significant correlation with scent hunting, whereas network 3, which involves regions that support movement, eye movement, and spatial navigation, shows a significant correlation with sight hunting. These findings strongly suggest that humans have altered the brains of different breeds of dogs in different ways through selective breeding.
It is important to note that the current study was performed on opportunistically acquired data. The dataset included different numbers of dogs from different breeds, and some breeds are not represented at all. We used permutation testing for statistical hypothesis testing, which is a nonparametric approach appropriate for differing group sizes, but it is still possible that different patterns of variation may have been obtained with a different sample makeup. Nonetheless, we expect the basic finding that this variation exists would remain.
Additionally, it should be noted that as dogs are increasingly bred to be house pets rather than working animals, selection on behavior is relaxing; significant behavioral differences have been found between working, show, and pet animals within a breed (Lofgren et al., 2014). To our knowledge, the dogs in the current study were all house pets. Therefore, the findings reported here should be taken as representative of the innate breed-typical adaptations to brain organization that emerge without the input of specific experience and may actually reflect relaxed or reduced versions of these adaptations. This might be akin to studying language circuitry in a lineage of language-deprived humans: humans almost certainly have some specialized “hard-wired” adaptations to this circuitry, but experience is required for the anatomical phenotype to fully emerge, and indeed it is difficult to consider language-related neural adaptations divorced from the context of language exposure and learning. Thus, future studies on purpose-bred dogs that are actively performing the tasks for which they are presumably adapted might expect to find additional or more pronounced neuroanatomical effects than we observed here.
These findings have relevance to both basic and applied science. First and foremost, our findings introduce neural variation in domestic dog breeds as a novel opportunity for studying the evolution of brain–behavior relationships. Dogs represent a “natural experiment” in behavioral selection that has been ongoing for thousands of years; it seems remarkable that attempts to observe the neurological results of this experiment have so far been fairly minimal. Our findings also have implications for the current proliferation of fMRI studies in pet dogs, which nearly always group together dogs of varying breeds. The current study suggests that this approach might not be ideal because there may be evolved breed differences in, for example, functional responses to stimuli or anatomical distribution of receptors. Consistent with this possibility, one study has already found that border collies and Siberian huskies respond significantly differently to intranasal oxytocin (Kovács et al., 2016). Additionally, on a practical level, our findings open the door to brain-based assessment of the utility of different dogs for different tasks. It might be possible, for example, to identify neural features that are linked to different breeds' specializations for specific behaviors, and to selectively breed or train dogs for enhanced expression of those neural features. Finally, on a philosophical level, these results tell us something fundamental about our own place in the larger animal kingdom: we have been systematically shaping the brains of another species.
Footnotes
This work was supported by the National Science Foundation–Division of Integrative Organismal Systems (Grant NSF-IOS 1457291). Additional support was provided by NIH OD P51OD11132 to the Yerkes National Primate Research Center. We appreciate the contributions of the veterinary and imaging staff at the UGA Veterinary Teaching Hospital.
The authors declare no competing financial interests.
- Correspondence should be addressed to Erin E. Hecht at erin_hecht{at}fas.harvard.edu