A Digital Atlas of the Dog Brain

Abstract

There is a long history and a growing interest in the canine as a subject of study in neuroscience research and in translational neurology. In the last few years, anatomical and functional magnetic resonance imaging (MRI) studies of awake and anesthetized dogs have been reported. Such efforts can be enhanced by a population atlas of canine brain anatomy to implement group analyses. Here we present a canine brain atlas derived as the diffeomorphic average of a population of fifteen mesaticephalic dogs. The atlas includes: 1) A brain template derived from in-vivo, T1-weighted imaging at 1 mm isotropic resolution at 3 Tesla (with and without the soft tissues of the head); 2) A co-registered, high-resolution (0.33 mm isotropic) template created from imaging of ex-vivo brains at 7 Tesla; 3) A surface representation of the gray matter/white matter boundary of the high-resolution atlas (including labeling of gyral and sulcal features). The properties of the atlas are considered in relation to historical nomenclature and the evolutionary taxonomy of the Canini tribe. The atlas is available for download (https://cfn.upenn.edu/aguirre/wiki/public:data_plosone_2012_datta).

About the Authors:

Ritobrato Datta

Affiliation: Department of Neurology, School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

Jongho Lee

Affiliation: Department of Radiology, School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

Jeffrey Duda

Affiliation: Department of Radiology, School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

Brian B. Avants

Affiliation: Department of Radiology, School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

Charles H. Vite

Affiliation: Section of Neurology, Department of Clinical Studies, School of Veterinary Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

Ben Tseng

Affiliation: Department of Neurology, School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

James C. Gee

Affiliation: Department of Radiology, School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

Gustavo D. Aguirre

Affiliation: Section of Ophthalmology, Department of Clinical Studies, School of Veterinary Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

Geoffrey K. Aguirre

* E-mail: [email protected]

Affiliation: Department of Neurology, School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America

Introduction

The domestic dog has served as an experimental model in neuroscience experiments and translational neurology for several centuries. Some of the earliest evidence for specific localization of brain function derived from experiments on dogs by Gustav Fritsch and Eduard Hitzig, who electrically stimulated small regions of the exposed cortex in awake animals [1], [2]. Using similar techniques, Sir David Ferrier identified multiple cortical areas related to the precise control of movement and translated these findings to map the “eloquent” cortex of patients with tumors undergoing neurosurgical procedures [3]. One of the earliest localizations of visual cortex was in the dog, identified using focal lesions [4], and in the early 20th century, the dog was used as a model of traumatic brain injury from missile wounds [5]. Perhaps the most celebrated use of dogs in neuroscience and psychology was the work of Ivan Pavlov that characterized conditioned reflexes [6], [7].

There has been a recent revival of interest in the canine as a model of ophthalmologic and neurologic disease. The dog has become an important model system for inherited retinal disease [8], and gene therapeutic treatment of these disorders (e.g. [9]–[11]). Dogs suffer from age-related cognitive dysfunction, and the associated neuropathology resembles human Alzheimer’s Disease [12]–[14]. The dog has also become a valuable model of inherited leukodystrophies [15], [16], and potential gene therapeutic treatment of lysosomal enzyme deficiencies [17], [18].

The dog continues to be an essential model of social cognition. Recent behavioral work in canines has examined the extent (and variability) of cognitive skills in different dog breeds, such as tracking cues [19], pointing gestures [20], and even “word learning” [21]. This interest in behavior and sensory function has led to a small but growing number of studies using functional and anatomical magnetic resonance imaging to study the canine brain. A set of early studies showed that visual stimulation in the anesthetized dog could produce measurable changes in blood oxygen level dependent (BOLD) fMRI signal from the canine visual cortex [22], [23]. Subsequently, the recovery of cortical responses following treatment of retinal disease by gene therapy was studied in the dog model [24]. A recent fMRI study has examined the neural correlates of reward mechanisms in the awake dog [25].

As the number of MRI-based studies of the canine brain grows, so does the need for a standard MRI-based template of the dog brain. Such a template allows data from across animals to be registered to a common space to be combined and compared, and facilitates quantitative comparisons of anatomical features. Here we present an atlas of the canine brain that is well suited for this purpose. This atlas is a diffeomorphic [26], population-based average that is composed of a low-resolution brain volume to be used for automated registration and skull-stripping, and a co-registered, high-resolution volume for data display and referencing of effects to the cortical surface.

Materials and Methods

Animals

A total of 15 dogs with mesocephalic (mesaticephalic) conformation were studied; all were purpose-bred, mixed-breed dogs that originated from several breeds having various forms of inherited retinal degenerations [27]. Thirteen of these animals were homozygous RPE65-mutants, resulting in severe retinal photoreceptor dysfunction present at birth, later treated successfully with subretinal injections of an adeno-associated viral vector carrying wild-type RPE65 [24]; the remaining 2 animals were non-affected carriers of the mutation with normal retinal function. Each animal was studied to obtain MRI images of brain anatomy either for the low-resolution, “in-vivo” atlas, or for the high-resolution, “ex-vivo” atlas.

Ethics Statement.

This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The protocol was approved by the Institutional Animal Care and Use Committee of the University of Pennsylvania (IACUC Protocol #s 803269 and 801870). All procedures were carried out under anesthesia and all efforts were made to minimize discomfort.

In-vivo Low Resolution Template

Magnetic resonance imaging.

T1-weighted images from seven dogs were used in the creation of the in-vivo template. Two, 15-minute MPRAGE images (1 mm isotropic) were acquired for each animal on a 3 Tesla Siemens Trio (Erlangen, Germany) using a transmit–receive, quadrature volume head coil (USA instruments, Aurora, Ohio). Each dog was anesthetized, and, during image acquisition, also paralyzed and ventilated (for details of anesthesia protocol, see [24]). The two MPRAGE images were subjected to 6-parameter realignment with least-squares minimization and then averaged (Figure 1A).

[Figure omitted. See PDF.]

Figure 1. In-vivo low resolution and ex-vivo high resolution templates.

(A) In-vivo T1-weighted images from three individual canines, obtained with 1 mm isotropic resolution. (B) Individual in-vivo brains following manual brain extraction. (C) Templates in the in-vivo space. Top is the diffeomorphic average, low-resolution brain including the soft-tissues of the head. Middle is the average, low-resolution, skull-stripped brain. Bottom is the high-resolution, ex-vivo diffeomorphic average following warping to the in-vivo space. (D) The high-resolution, diffeomorphic average of the ex-vivo brains, in the ex-vivo space. (E) Examples of high-resolution, ex-vivo brains scanned at 7 Tesla with 0.33 mm isotropic resolution.

https://doi.org/10.1371/journal.pone.0052140.g001

Creation of low-resolution, in-vivo anatomical template.

The anatomical image for each animal was segmented to separate the brain from the skull (Figure 1B) using semi-automated, open-source methods available in ITK SNAP (http://www.itksnap.org/pmwiki/pmwiki.php) [28]. A cascade of transformations were then applied to generate an unbiased shape and intensity template using the Advanced Normalization Tools (ANTs) (http://www.picsl.upenn.edu/ANTS/). The initial distribution of the skull-stripped brains was estimated by directly averaging their intensities. The second step involved rigidly registering each image to the intensity distribution using ITK mutual information as a similarity metric. Then a trimmed average of the rigidly registered intensities was used to tighten the distribution [29]. Next, an elastic registration model was applied to find a sharper trimmed average intensity image and a small deformation shape average. The final step used large deformation diffeomorphic image registration and shape averaging to bring all structures of the brains into exact correspondence which yielded the final unbiased intensity average and the final optimal shape anatomy [26]. The middle panel of Figure 1C (labeled low-res) is the un-biased, diffeomorphic template average derived from the low-resolution (1 mm isotropic), in-vivo set of brains.

We then created a template in the same space as the low-res atlas, which retained the skull and head. To do so, individual skull-striped canine brains were diffeomorphically registered to the low-res atlas. The transform generated while warping each individual skull-striped brain to the low-res atlas was then used to warp the respective individual MPRAGE image containing the entire head to the low-res atlas space. A diffeomorphic image registration and shape averaging was then performed on the set of registered, whole-brain images to create the in-vivo atlas that includes both the brain and soft tissues of the head and skull (Figure 1C; labeled low-res w/ skull).

A description of the creation of the low-res atlas has been presented previously in abstract form [30].

Ex-vivo High Resolution Template

Brain collection.

Brains were collected from a separate group of eight animals following euthanasia (intravenously administered Euthasol; Virbac Animal Health, Ft. Worth, TX) after completion of gene-therapy studies. The skin and muscles overlying the skull were removed, and an oscillating saw used to cut the calvarium which was elevated bluntly with a scalpel handle and removed. The overlying meninges were removed, and then the head was rotated 180° to provide better exposure of the ventral aspect of the brain. All cranial nerves were cut with fine, blunt scissors, and the brain removed with minimal damage. To prevent compression during fixation, the brain was placed in 10% buffered formalin, and 37% formalin stock solution was added until the brain floated just below the fluid surface; paper towels soaked in formalin were placed over the surface to prevent drying from exposure. After 48–72 hours the brains were transferred to 10% formalin solution where they were stored for 2 to 3 months prior to the MRI studies. Prior to MRI scanning, the brains were transferred to phosphate buffered saline (PBS), and the fluid changed every 3–4 days for 3 changes.

Magnetic resonance imaging.

MRI images were acquired on a 7 Tesla whole body MRI system (Siemens, Erlangen, Germany) with a 32 channel phased-array head coil (Nova Medical, Wilmington, PA). The brain was stored in a cylinder filled with PBS and placed at the bottom of the coil to improve SNR. T1-weighted MPRAGE images were acquired. The brain was covered by two sequentially acquired slabs with the middle area overlapped in both slabs. The resolution was 0.33 mm isotropic, FOV  = 84×84×34.3 mm in each slab, matrix size  = 256×256×104 in each slab, TR  = 3 sec, TI  = 550 ms, TE  = 3.4 ms, flip angle  = 12°, pixel bandwidth  = 370 Hz/pixel, and total scan time  = 12:48 min for each slab. The acquisition was repeated 6 times for signal averaging. While there was no veridical movement of the studied tissues, drift of the image within the field of view can occur with warming of the gradient coils. Therefore, the six MPRAGE images were subjected to 6-parameter realignment with least-squares minimization and then averaged (Figure 1E).

Creation of high-resolution, ex-vivo anatomical template.

Creation of the high-res atlas proceeded in a manner similar to that used for the low-res atlas. Each T1-weighted structural MRI of the ex-vivo canine brain was mapped using a cross-correlation registration metric to an optimal template space, defined as the population-specific, unbiased average shape and appearance image derived from a representative population which in this study are the individual ex-vivo brains [29].

Without the confinement of the skull, the ex-vivo brains relax into wider left-right conformation. This is reflected in the diffeomorphic average (Figure 1D). The template ex-vivo brain was therefore mapped using a mutual information registration metric to the in-vivo, low-res template. This yielded the final, high-res atlas which resides within the in-vivo template space.

Canine brain inflated surface.

The high-res template was then processed using a modified version of the automatic anatomical surface reconstruction pipeline of the FreeSurfer toolkit (http://surfer.nmr.mgh.harvard.edu/) [31], [32]. After automatic tissue segmentation, the images were manually inspected to identify errors in the gray / white matter boundary definition. The gray / white matter intensity differences were very large for the majority of the areas in the canine brain and the corresponding tissue segmentation was generally accurate. However, some regions contained partial volume effects and required manual intervention to demarcate the boundary. Control points for white matter voxels were manually defined in areas of problematic tissue segmentation, and the FreeSurfer pipeline re-invoked to estimate the gray white matter boundary. This was performed iteratively until the segmentation was judged sufficiently accurate. Topological holes in the white matter segmentation produced by the ventricles and hippocampi were manually filled. The two hemispheres of the filled white matter volume were then separated. Separate three-dimensional rendering of each white matter hemispheric volume was created. The surface of the three dimensional rendering is the gray matter / white matter boundary which was then smoothed with a surface smoothing kernel to create a smoothed white matter surface. The white matter surface was then inflated using the standard tools available in FreeSurfer.

Results and Applications

The complete canine atlas set is composed of co-registered, low-res and high-res volumetric templates, including a low-res template that includes the skull and soft tissues of the head. Further, the high-res volumetric template serves as the basis of a cortical surface reconstruction of the canine brain. Below we examine the availability of sub-cortical and cortical detail in the template brain; the pattern, nomenclature, and evolutionary history of canine cortical surface topology; and suggested processing approaches for use of the canine atlas.

Detail

In the high resolution (0.33 mm isotropic), ex-vivo MRI images collected with 2½ hours of scanning at 7 Tesla, excellent contrast was available between gray and white matter. While individual differences in anatomy would be expected to induce smoothing of the high-resolution images when averaged across subjects, shape-based diffeomorphic registration matches tissue types prior to averaging. Consequently, the fine detail present in individual, high-resolution brain images is well preserved in the average atlas. Figure 2 illustrates some of the subcortical and brainstem anatomical features that are visible in the atlas. Notable is the clear appearance of the claustrum (Figure 2A, label b), which is a thin strip of gray matter located between the external and extreme capsules, and the preservation of the folia of the cerebellum. For some of these structures (e.g., the thalamus seen in axial view in Figure 2A, label d) further structure is readily apparent.

[Figure omitted. See PDF.]

Figure 2. Representative slices of the high resolution ex-vivo template demonstrating labeled cortical and subcortical structures.

(A) Dorsal plane (horizontal) slice through the basal ganglia and thalamus. The fine structure of both the lateral geniculate nucleus and head of the hippocampus can be seen in this population average image. (B) Expanded and contrast-enhanced coronal slice through the brainstem, and illustration of white matter tissue segmentation.

https://doi.org/10.1371/journal.pone.0052140.g002

In the brainstem, differentiation between gray matter and white matter structures can also be seen. We have labeled some of these anatomical features (Figure 2B, following the [33], [34]). Our goal in doing so is not to provide a comprehensive atlas of all brainstem structures (for an excellent reference for this purpose, see Palazzi [34]) but instead to illustrate that the diffeomorphic average contains sufficient anatomical detail to support such efforts.

The imaging contrast between the gray and white matter enabled the definition of a white matter tissue segmentation (Figure 2B, right). This was iteratively edited to ensure that the white matter volume in each hemisphere was a continuous volume without topological defects, and thus may be expressed as a continuous cortical surface within FreeSurfer.

Canine Cortical Surface Topology

We produced a surface rendering of the canine brain from the high-res atlas (Figure 3A). Next, following segmentation of the topologically corrected white matter, a surface based reconstruction of the canine hemisphere was performed within FreeSurfer. The resulting inflated view of the cortical surface (Figure 3B) allows the continuous cortical sheet to be seen, including cortex normally obscured within the sulcal depths. We labeled the sulci and gyri on the cortical surface, generally following the nomenclature of Miller et al [35], [36]. It should be noted that the olfactory bulb, which is a prominent feature in the canine brain, is absent in this ex-vivo atlas as this structure was transected in removal of the brain.

[Figure omitted. See PDF.]

Figure 3. Surfaces and labels.

(A) Reconstruction of the canine brain surface from the high-res atlas, viewed in three orientations. The location of the prorean gyrus, coronal sulcus (red) and ansate sulcus (blue) is indicated on the dorsal view. (B) Lateral and medial views of the inflated white matter surface with sulci and gyri labeled. The dark gray structures are the sulci and the light gray regions are the gyri. On the labeled surfaces, the sulci are colored in less saturated colors and gyri in saturated colors.

https://doi.org/10.1371/journal.pone.0052140.g003

There are disagreements of nomenclature for the canine cortical surface. A prominent variation regards the most medial, dorsal, caudal gyrus. We have adopted the label “marginal gyrus” [37]–[39] instead of the alternate “lateral gyrus” [40]–[44] for this structure, which is the location of primary visual cortex in the dog (by homology to the cat [45] and sheep [46]. Cohn and Papez [47] adopted the term “calcarine sulcus” for the structure that we have labeled the “splenial sulcus” and further described a “posterior calcarine fissue”, extending caudally from the splenial sulcus in approximately half the dogs they studied. We did not observe this structure in our atlas, although we cannot exclude that this is a breed specific difference.

The pattern of cortical folding has been used to provide taxonomic organization of the family Canidae (which includes dogs, wolves, and foxes, among many other extant and extinct species). Based upon the observed size of the prorean gyrus (also termed the proreal gyrus; Figure 3A), Huxley [48], and subsequently Radinsky [49], divided canids into dogs and foxes. Lyras [50], [51] studied endocasts of 29 different living and extinct species of the Canini subfamily (which includes wolves and dogs, but not foxes), and suggested that the overall conformation of the coronal and ansate sulci (Figure 3A) assumes one of four basic patterns which distinguishes among the genera within Canini. The appearance of these sulci in our atlas best corresponds to the “orthogonal” arrangement. This, along with the well formed and elongated prorean gyrus in the atlas, is as expected in the genus Canis, of which the domestic dog subspecies is a member.

Application Approaches

In practice, if an experimenter has collected T1-weighted anatomical images and functional data, then the workflow described below can be used to warp individual animal data to the template space for group analyses.

First, the individual subject brain with the skull is registered to the low-res w/ skull template (Figure 1C, top) using, e.g., the diffeomorphic warping tool available in ANTS. The resulting transformation matrix is then used to project a binary tissue mask (brain vs. not-brain) from low-res w/ skull template back to the individual brain space; the individual brain in template space can be extracted from its skull by applying this binary mask. Manual editing of the segmentation mask may be performed at this step if necessary.

Once a satisfactory result is achieved, the extracted brain in the original, individual animal space is registered to the low-res template. The resulting transformation matrix may be applied to raw functional data or statistical maps in the original space as well.

We expect the low-res template to serve as the best target for registration of anatomical images collected in-vivo in individual animals, as these will be the most similar in contrast and detail. Because the low-res and hi-res atlases are in register, data may be referred to the hi-res atlas following registration to the template space, and further displayed on the cortical surface reconstruction.

Discussion

The canine atlas was created using diffeomorphic registration of a population of brains, initialized with an intensity average. This approach provides the benefit of shape-based averaging, which guarantees that tissues are in correspondence prior to averaging, while avoiding the bias of using one individual from the set as a template norm [26]. Consequently, the canine digital atlas has two important properties. First, the image appearance is not driven by any specific anatomical structure (as no manual landmarking was required); and second, the image shape is independent of any individual’s anatomical coordinate system [52]. The automated image registration methods used for the generation of the templates assume that the structural correspondences are correct between the images of the different canine brains. This assumption seems well justified as all the animals studied, although mixed-breed, were of a common genetic background [27], and thus would be expected to have similar cortical structure.

The high-resolution anatomical images were obtained from ex-vivo brains. Because the brain changes shape when freed from the confines of the skull, the ex-vivo atlas is an imperfect target for registration of in-vivo images. To mitigate this limitation, we created a volumetric brain template from lower-resolution, in-vivo brain images, and then transformed the ex-vivo average to the in-vivo space. Presuming that the alteration of brain shape produced by skull removal is well captured by the plastic deformation used for co-registration of the low and high-res atlases, this distortion should be fully corrected in the our atlas.

An edge artifact is present in the high-res images acquired at 7 Tesla, consisting of a non-uniform, T1 hyperintense band at the external edges of the gray matter ribbon (best seen in the axial and coronal slices in Figure 2). A possible cause of this band is a long MPRAGE readout that induces different T1 weighting for high spatial frequencies at the gray matter edge. This theory is supported by the finding of two separate bands in the gray matter when a still longer readout was used. An alternative explanation is the effect of chemical fixative [53]. Regardless of the cause, we believe that the hyperintense rim is properly segmented as gray matter, based upon comparison to the low-res images obtained at 3 Tesla. As this artifact is restricted to the outer edge of the cortical sheet, it did not compromise the construction of the white matter label.

We anticipate several applications of the atlas to the analysis of canine neuroimaging data. The atlas may be used to register functional data from different animals to a common anatomical space, allowing group-level inferences (such as was conducted in [24]). Individual differences in brain structure as assessed by different imaging modalities, such as cortical thickness, or diffusion tensor imaging [54], may be related to normal variations of behavior or one of many disease states. Indeed, the canine is a valuable model system for many neurological diseases, including epilepsy [55]; cortical malformations such as lissencephaly [56], [57] and polymicrogyria [58], [59]; dementia [60], [61]; and focal lesions [62], [63]. Given the good registration of high-resolution anatomy with a head model, the atlas can be used to guide source localization of EEG recording in the dog [64].

An important feature of the canine atlas, and a potential limit to these applications, is that it was derived from dogs with mesaticephalic skull conformation, meaning that the skull is “medium” shaped, as opposed to elongated (dolichocephalic; e.g., Greyhounds) or shortened (brachycephalic; e.g., Boxers). Head shape may influence cortical folding in a manner more complex than a simple affine transformation of brain size. Differences in canine skull shape are associated with different sensory and behavioral profiles [65]. For example, the distribution of retinal ganglion cells differs between breeds based on muzzle length [66], presumably related to differences in the extent of binocular vision. Other studies have noted an association between head shape and biomechanical function, with brachycephalic breeds being used as guards and fighters and and dolichocephalic breeds as runners [67]. At the very least, there is evidence that human directed breeding of dogs has produced systematic differences in canine cerebral organization, for example the position of the olfactory lobe [68]. Given breed-specific differences in behavior and skull shape, it is quite possible that cortical topology also differs between breeds and in turn relates to behavioral diversity. Therefore, the canine template offered here should be used with caution in the analysis of data obtained from non-mesaticephalic animals. Conversely, our atlas can serve as a starting point to test for the existence of such differences. For example, brain anatomical images from different breeds may be registered, perhaps using explicit sulcal topology to the surface template, and systematic differences in surface deformation assessed [69].

Another potential limitation of the atlas is that it was derived from animals with a visual impairment. Thirteen of 15 animals used to construct the atlas were born with a severe, congenital form of retinal blindness. This canine analog of Leber’s Congenital Amaurosis is a rod-cone dysfunction caused by mutations in the RPE65 gene. We might consider that there exists systematically different brain structure in either RPE65 mutants specifically or animals with congenital blindness generally. We believe the effect upon our atlas is minimal, however. First, RPE65 is not expressed within the central nervous system [70]. Further, the neural retina is not altered in the disease [71], thus preserving the optic nerves and post-chiasmatic anatomy [24]. Second, while congenital blindness (where it has been studied in the human) can alter brain structure, these effects are generally subtle, such as changes in gray matter thickness [72] or the surface area [73] of the striate cortex. In both the current study, and our previous examination of controls and RPE65-mutants [24], no qualitative structural differences were observed.

While the dog brain is gyrencephalic (characterized by the development of sulci and gyri) [74]–as opposed to the lisencephalic (smooth) brains of some mammals and birds–the shape and appearance of the sulci and gyri is generally uniform across individuals [75]; although, see [47]. This implies that diffeomorphic volume based registration approaches, as used here in template construction, are sufficient to register individual brains to the template without the need for surface-based approaches. However, if the degree of variability of cortical folding pattern is subsequently found to be greater between or within breeds, then a surface based registration would be preferred over volume based approaches [76]–[77].

An additional variation of cortical surface topology that remains to be examined is hemispheric asymmetry. Cerebral asymmetry may be a fundamental feature of vertebrates [78], and there is some evidence of anatomical [79], [80] and behavioral laterality [81], [82] in the dog. For example, the right hemisphere is larger than the left hemisphere [79], although specific features, such as the ectosylvian gyrus, is larger on the left [80]. The canine brain atlas we have presented here can serve as the basis of quantitative hemispheric comparisons, using volumetric techniques following mirror reversal [83], [84], or surface-based approaches following the creation of a pseudo-hemisphere that has surface topology intermediate between the left and right [85].

Digital atlases of the healthy (http://vanat.cvm.umn.edu/mriBrainAtlas/) and diseased [39], [86] canine brain have been offered previously. A notable previous effort is the work of Tapp and colleagues [87], which constructed a dog brain template from the average of 192 animals and then used the template in a voxel-based morphometry study of the effect of aging upon the canine brain. The atlas we have created is an advance on these previous efforts in several respects. In addition to basic improvements in voxel resolution and imaging contrast, our atlas has the important property of representing the diffeomorphic central tendency of a group of animals, having both low and high-resolution versions to support an image analysis pipeline, and a surface-based implementation. The ex-vivo images at 7 Tesla provided good spatial and contrast resolution for the identification of anatomical features, which can be difficult to obtain in-vivo using clinical scanners operating at lower field strengths (although see [88], [89] for progress in obtaining in-vivo canine measures at high field). Our atlas is free to use (with appropriate attribution) for academic or commercial purposes, although the atlas may not be distributed for commercial gain. It may be downloaded from our website (https://cfn.upenn.edu/aguirre/wiki/public:data_plosone_2012_datta).

Acknowledgments

The authors are grateful to Mr. Joe Magrane for superb technical skills in isolating the brains for the high resolution, ex-vivo studies.

Author Contributions

Conceived and designed the experiments: JL GDA GKA. Performed the experiments: JL GDA GKA. Analyzed the data: RD JD BBA CHV BT JCG GKA. Contributed reagents/materials/analysis tools: BBA JCG. Wrote the paper: RD GKA.

Citation: Datta R, Lee J, Duda J, Avants BB, Vite CH, Tseng B, et al. (2012) A Digital Atlas of the Dog Brain. PLoS ONE 7(12): e52140. https://doi.org/10.1371/journal.pone.0052140

References

1. Fritsch G, Hitzig E (1870) Ueber die elektrische Erregbarkeit des Grosshirns. Archiv für Anatomie, Physiologie und wissenschaftliche Medicin: 300–332.

2. Bonin von G (1960) Some Papers on the Cerebral Cortex. Translated from the French and German by Gerhardt von Bonin.: Charles C Thomas, Publisher, Springfield, Illinois. xxiv, 396 pp.

3. Ferrier D (1886) The functions of the brain. 2nd ed. G.P. Putnam’s Sons. 498 pp.

4. Munk H (1878) Weitere Mittheilungen zur Physiologie der Grosshirnrinde. Archiv fur Anatomie und Physiologie, 2 161–178.

5. Horsley V (1915) Remarks on Gunshot Wounds of the Head: Made in Opening a Discussion at the Medical Society of London on February 8th, 1915. Br Med J 1: 321–323.

6. Razran HS, Warden CJ (1929) The sensory capacities of the dog as studied by the conditioned reflex method (Russian schools). Psychological Bulletin 26: 202.

7. Pavlov IP (1995) Lectures on Conditioned Reflexes. (Twenty-Five Years of Objective Study of the Higher Nervous Activity (Behaviour) of Animals. Limited. Classics of Medicine Library.

8. Miyadera K, Acland GM, Aguirre GD (2011) Genetic and phenotypic variations of inherited retinal diseases in dogs: the power of within- and across-breed studies. Mamm Genome: –. doi:10.1007/s00335-011-9361-3.

9. Acland GM, Aguirre GD, Ray J, Zhang Q, Aleman TS, et al. (2001) Gene therapy restores vision in a canine model of childhood blindness. Nat Genet 28: 92–95 .

10. Komáromy AM, Alexander JJ, Rowlan JS, Garcia MM, Chiodo VA, et al. (2010) Gene therapy rescues cone function in congenital achromatopsia. Hum Mol Genet 19: 2581–2593 .

11. Beltran WA, Cideciyan AV, Lewin AS, Iwabe S, Khanna H, et al. (2012) Gene therapy rescues photoreceptor blindness in dogs and paves the way for treating human X-linked retinitis pigmentosa. Proc Natl Acad Sci USA 109: 2132–2137 .

12. González-Martínez Á, Rosado B, Pesini P, Suárez M-L, Santamarina G, et al. (2011) Plasma β-amyloid peptides in canine aging and cognitive dysfunction as a model of Alzheimer’s disease. Exp Gerontol 46: 590–596 .

13. Sarasa MM, Pesini PP (2009) Natural non-trasgenic animal models for research in Alzheimer’s disease. CORD Conference Proceedings 6: 171–178.

14. Cotman CW, Head E (2008) The canine (dog) model of human aging and disease: dietary, environmental and immunotherapy approaches. J Alzheimers Dis 15: 685–707.

15. Vite CHC, McGowan JCJ (2001) Magnetization transfer imaging of the canine brain: a review. Vet Radiol Ultrasound 42: 5–8 .

16. McGowan JC, Haskins M, Wenger DA, Vite C (2000) Investigating demyelination in the brain in a canine model of globoid cell leukodystrophy (Krabbe disease) using magnetization transfer contrast: preliminary results. J Comput Assist Tomogr 24: 316–321.

17. Ellinwood NM, Vite CH, Haskins ME (2004) Gene therapy for lysosomal storage diseases: the lessons and promise of animal models. J Gene Med 6: 481–506 .

18. Haskins MM (2009) Gene therapy for lysosomal storage diseases (LSDs) in large animal models. ILAR J 50: 112–121.

19. Téglás E, Gergely A, Kupán K, Miklósi A, Topál J (2012) Dogs’ gaze following is tuned to human communicative signals. Curr Biol 22: 209–212 .

20. Miklósi A, Soproni K (2006) A comparative analysis of animals’ understanding of the human pointing gesture. Anim Cogn 9: 81–93 .

21. Kaminski J, Call J, Fischer J (2004) Word learning in a domestic dog: evidence for “fast mapping”. Science 304: 1682–1683 .

22. Willis CKC, Quinn RPR, McDonell WMW, Gati JJ, Partlow GG, et al. (2001) Functional MRI activity in the thalamus and occipital cortex of anesthetized dogs induced by monocular and binocular stimulation. Can J Vet Res 65: 188–195.

23. Willis CK, Quinn RP, McDonell WM, Gati J, Parent J, et al. (2001) Functional MRI as a tool to assess vision in dogs: the optimal anesthetic. Vet Ophthalmol 4: 243–253 .

24. Aguirre GK, Komáromy AM, Cideciyan AV, Brainard DH, Alemán TS, et al. (2007) Canine and human visual cortex intact and responsive despite early retinal blindness from RPE65 mutation. PLoS Med 4: e230 .

25. Berns GS, Brooks AM, Spivak M (2012) Functional MRI in awake unrestrained dogs. PLoS ONE 7: e38027 .

26. Avants BB, Epstein CL, Grossman M, Gee JC (2008) Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain. Med Image Anal 12: 26–41 .

27. Aguirre GD, Acland GM (2007) Models, mutants and man: searching for unique phenotypes and genes in the dog model of inherited retinal degeneration. In Ostrander EA, et al. The Dog and Its Genome. CSHL Press, 584 pp.

28. Yushkevich PA, Piven J, Hazlett HC, Smith RG, Ho S, et al. (2006) User-guided 3D active contour segmentation of anatomical structures: Significantly improved efficiency and reliability. Neuroimage 31: 1116–1128 .

29. Avants B, Gee JC (2004) Geodesic estimation for large deformation anatomical shape averaging and interpolation. Neuroimage 23 Suppl 1S139–S150 .

30. Avants B, Aguirre G, Walker J, Gee JC (2008) Unbiased Diffeomorphic Shape and Intensity Atlas Creation: Application to Canine Brain. ISMRM Abstract.

31. Dale AM, Fischl B, Sereno MI (1999) Cortical surface-based analysis. I. Segmentation and surface reconstruction. Neuroimage 9: 179–194 .

32. Fischl B, Sereno MI, Dale AM (1999) Cortical surface-based analysis. II: Inflation, flattening, and a surface-based coordinate system. Neuroimage 9: 195–207 .

33. Buxton DF, Compton RW (1986) The canine brain: Basic atlas for an auto-tutorial approach to the central nervous system. S.I.

34. Palazzi X (2011) The Beagle Brain in Stereotaxic Coordinates. Springer

35. Miller ME, Christensen G, Evans H (1965) ANATOMY OF THE DOG. Reprint. W. B. Saunders Co. 941 pp.

36. Evans HE (1993) Miller’s Anatomy of the Dog, 3rd ed. Saunders.

37. Ofri R (1993) Electrophysiological and Histological Mapping of the Cortical Area of Central Vision in the Dog. Doctoral Dissertation at University of Florida, 1993 - 218 pp.

38. Leigh EJ, Mackillop E, Robertson ID, Hudson LC (2008) Clinical Anatomy of the Canine Brain Using Magnetic Resonance Imaging. Vet Radiol Ultrasound 49: 113–121 .

39. Mogicato G, Conchou F, Raharison F, Sautet J (2011) Normal canine brain: comparison between magnetic resonance imaging and cross-sectional anatomy. Rev Med Vet-Toulouse 162: 400–405.

40. Langley JN (1883) The Structure of the Dog’s Brain. J Physiol (Lond) 4: 248–326.3.

41. Lim RKS, Chan-Nao Liu, Moffitt RL (1960) A stereotaxic atlas of the dog’s brain. Thomas, 93 pp.

42. Black AH, Young GA (1972) Electrical activity of the hippocampus and cortex in dogs operantly trained to move and to hold still. Journal of comparative and physiological psychology 79: 128.

43. Hoerlein BF (1978) Canine Neurology: Diagnosis and Treatment. 3rd ed. W.B. Saunders Company. pp.

44. Stein L, Roy K, Lei L, Kaushal S (2011) Clinical gene therapy for the treatment of RPE65-associated Leber congenital amaurosis. Expert Opin Biol Ther 11: 429–439 .

45. Sereno MI, Allman JM (1991) Cortical visual areas in mammals. The neural basis of visual function. The Neural Basis of Visual Function. Macmillan, 160–172 pp.

46. Clarke PG, Whitteridge D (1976) The cortical visual areas of the sheep. J Physiol (Lond) 256: 497–508.

47. Cohn HA, Papez JW (1933) The posterior calcarine fissure in the dog. The Journal of Comparative Neurology 58: 593–602.

48. Huxley TH (1880) 3. On the Cranial and Dental Characters of the Canidæ. Proceedings of the Zoological Society of London 48: 238–288 .

49. Radinsky L (1973) Evolution of Canid Brain. Brain Behav Evol 7: 169–202.

50. Lyras GA, Van Der Geer AAE (2003) External brain anatomy in relation to the phylogeny of Caninae (Carnivora: Canidae). Zoological Journal of the Linnean Society 138: 505–522.

51. Lyras GA (2009) The evolution of the brain in Canidae (Mammalia: Carnivora). Scripta Geologica, 139: 1–93 pp.

52. Ashburner J, Friston KJ (2009) Computing average shaped tissue probability templates. Neuroimage 45: 333–341 .

53. Cahill LS, Laliberté CL, Ellegood J, Spring S, Gleave JA, et al. (2012) Preparation of fixed mouse brains for MRI. Neuroimage 60: 933–939 .

54. Wang P, Zhu JM (2010) Quantitative Diffusion Tensor Imaging of White Matter Microstructure in Dog Brain at 7 T. The Open Medical Imaging Journal, 4, 1–5.

55. Kuwabara T, Hasegawa D, Kobayashi M, Fujita M, Orima H (2010) Clinical magnetic resonance volumetry of the hippocampus in 58 epileptic dogs. Vet Radiol Ultrasound 51: 485–490 .

56. Lee K-I, Lim C-Y, Kang B-T, Park H-M (2011) Clinical and MRI findings of lissencephaly in a mixed breed dog. J Vet Med Sci 73: 1385–1388.

57. MacKillop EE (2011) Magnetic resonance imaging of intracranial malformations in dogs and cats. Vet Radiol Ultrasound 52: S42–S51 .

58. Vanwinkle TJ, Fyfe JC, Dayrellhart B, Aguirre GD, Acland GM, et al. (1994) Blindness Due to Polymicrogyria and Asymmetrical Dilation of the Lateral Ventricles in Standard Poodles. Prog Vet Neurol 5: 66–71.

59. Jurney C, Haddad J, Crawford N, Miller AD, Van Winkle TJ, et al. (2009) Polymicrogyria in standard poodles. J Vet Intern Med 23: 871–874 .

60. Hasegawa DD, Yayoshi NN, Fujita YY, Fujita MM, Orima HH (2005) Measurement of interthalamic adhesion thickness as a criteria for brain atrophy in dogs with and without cognitive dysfunction (dementia). Vet Radiol Ultrasound 46: 452–457 .

61. Hasegawa D, Yamato O, Nakamoto Y, Ozawa T, Yabuki A, et al. (2012) Serial MRI Features of Canine GM1 Gangliosidosis: A Possible Imaging Biomarker for Diagnosis and Progression of the Disease. The Scientific World Journal 2012: 1–10 .

62. Vite CHC, Cross JRJ (2011) Correlating magnetic resonance findings with neuropathology and clinical signs in dogs and cats. Vet Radiol Ultrasound 52: S23–S31 .

63. Wolff CA, Holmes SP, Young BD, Chen AV, Kent M, et al. (2012) Magnetic resonance imaging for the differentiation of neoplastic, inflammatory, and cerebrovascular brain disease in dogs. J Vet Intern Med 26: 589–597 .

64. Pellegrino F (2004) Canine electroencephalographic recording technique: findings in normal and epileptic dogs. Clinical Neurophysiology 115: 477–487 .

65. Helton WS (2009) Cephalic index and perceived dog trainability. Behav Processes 82: 355–358 .

66. McGreevy P, Grassi TD, Harman AM (2004) A strong correlation exists between the distribution of retinal ganglion cells and nose length in the dog. Brain Behav Evol 63: 13–22 .

67. Ellis JL, Thomason J, Kebreab E, Zubair K, France J (2009) Cranial dimensions and forces of biting in the domestic dog. J Anat 214: 362–373 .

68. Roberts T, McGreevy P, Valenzuela M (2010) Human induced rotation and reorganization of the brain of domestic dogs. PLoS ONE 5: e11946 .

69. Fischl B, Rajendran N, Busa E, Augustinack J, Hinds O, et al. (2008) Cortical folding patterns and predicting cytoarchitecture. Cerebral Cortex 18: 1973–1980 .

70. Hamel CP, Tsilou E, Harris E, Pfeffer BA, Hooks JJ, et al. (1993) A developmentally regulated microsomal protein specific for the pigment epithelium of the vertebrate retina. J Neurosci Res 34: 414–425 .

71. Hernández MM, Pearce-Kelling SES, Rodriguez FDF, Aguirre GDG, Vecino EE (2010) Altered expression of retinal molecular markers in the canine RPE65 model of Leber congenital amaurosis. Invest Ophthalmol Vis Sci 51: 6793–6802 .

72. Jiang J, Zhu W, Shi F, Liu Y, Li J, et al. (2009) Thick Visual Cortex in the Early Blind. Journal of Neuroscience 29: 2205–2211 .

73. Park H-J, Lee JD, Kim EY, Park B, Oh M-K, et al. (2009) Morphological alterations in the congenital blind based on the analysis of cortical thickness and surface area. Neuroimage 47: 98–106 .

74. Budras KD (2007) Anatomy of the Dog. Manson Publishing, 224 pp.

75. Atkins DL (1978) Evolution and morphology of the coyote brain. In: Bekoff, M. (ed.), Coyotes: Behavior, Biology, and Management. Academic Press, New York. 17–35 pp.

76. Klein A, Ghosh SS, Avants B, Yeo BTT, Fischl B, et al. (2010) Evaluation of volume-based and surface-based brain image registration methods. Neuroimage 51: 214–220 .

77. Ghosh SS, Kakunoori S, Augustinack J, Nieto-Castanon A, Kovelman I, et al. (2010) Evaluating the validity of volume-based and surface-based brain image registration for developmental cognitive neuroscience studies in children 4 to 11 years of age. Neuroimage 53: 85–93 .

78. Bisazza AA, Rogers LJL, Vallortigara GG (1998) The origins of cerebral asymmetry: a review of evidence of behavioural and brain lateralization in fishes, reptiles and amphibians. Neurosci Biobehav Rev 22: 411–426.

79. Tan U, Calişkan S (1987) Asymmetries in the cerebral dimensions and fissures of the dog. Int J Neurosci 32: 943–952.

80. Natchev S, Altanjiiski I (2000) Right-Left Asymmetry of the Temporoparieta! Area (Tpt) in Dog Brain: I. Cytoarchitectonic and Quantitative Study. Comptes Rendus de l’Academie Bulgare des Sciencesi, vol. 53 1: 121.

81. Wells DL (2003) Lateralised behaviour in the domestic dog, Canis familiaris. Behav Processes 61: 27–35.

82. McGreevy PD, Brueckner A, Thomson PC, Branson NJ (2010) Motor laterality in 4 breeds of dog. Journal of Veterinary Behavior: Clinical Applications and Research 5: 318–323 .

83. Good CD, Johnsrude I, Ashburner J, Henson RN, Friston KJ, et al. (2001) Cerebral Asymmetry and the Effects of Sex and Handedness on Brain Structure: A Voxel-Based Morphometric Analysis of 465 Normal Adult Human Brains. Neuroimage 14: 685–700 .

84. Luders E, Gaser C, Jancke L, Schlaug G (2004) A voxel-based approach to gray matter asymmetries. Neuroimage 22: 656–664 .

85. Greve DN, Sabuncu MR, Shafee R, Schmansky N, Buckner RL, Fischl B (2011) Automatic Surface-based Interhemispheric Registration with FreeSurfer. OHBM Abstract.

86. Kraft SL, Gavin PR, Wendling LR, Reddy VK (2005) Canine Brain Anatomy on Magnetic Resonance Images. Veterinary Radiology 30: 147–158 .

87. Tapp PD, Head K, Head E, Milgram NW, Muggenburg BA, et al. (2006) Application of an automated voxel-based morphometry technique to assess regional gray and white matter brain atrophy in a canine model of aging. Neuroimage 29: 234–244 .

88. Kang B-TB, Ko K-JK, Jang D-PD, Han J-YJ, Lim C-YC, et al. (2009) Magnetic resonance imaging of the canine brain at 7 T. Vet Radiol Ultrasound 50: 615–621 .

89. Martín-Vaquero PP, Da Costa RCR, Echandi RLR, Tosti CLC, Knopp MVM, et al. (2011) Magnetic resonance imaging of the canine brain at 3 and 7 T. Vet Radiol Ultrasound 52: 25–32.