Creation of the WHS rat brain atlas
The Waxholm Space (WHS) rat brain atlas consist of four parts: a reference data set, a coordinate system, an annotation set, and a terminology (Fig. 5a). The reference data set is volumetric and has a coordinate system based on internal landmarks, 222 three-dimensional annotations representing brain regions manually drawn. The terminology is largely consistent with existing reference atlases but also includes additional expert neuroanatomist opinions. In the following, we present the methods for the generation of the different parts of the atlas.
Reference data and coordinate system
The reference data used as a basis for the WHS rat brain atlas comprise high-resolution, contrast-enhanced T2*-weighted gradient recalled echo anatomic images with an isotropic spatial resolution of 39 μm (sMRI), and colour-coded principal diffusion direction maps (DTI) with a resolution of 78 μm (Fig. 5a1). These data were acquired ex vivo from the head of an adult male Sprague Dawley rat (age 80 days, weight 397.6 g, Charles River, Wilmington, MS, USA), perfusion-fixed using a mixture of formalin and a gadolinium-based MRI contrast agent (Johnson et al., 2002). Imaging was performed using a 7T small animal MRI system (Magnex Scientific, Yarnton, Oxford, UK) at the Duke Centre for In Vivo Microscopy (Durham, NC, USA). Technical details are provided in Papp et al. (2014, 2015).
The WHS is a continuous, three-dimensional Cartesian coordinate system with the origin defined at the decussation of the anterior commissure (Fig. 5a2). The origin of the WHS, as first described by Papp and colleagues (page 376)16: “at the intersection of: a) the mid-sagittal plane, b) a coronal plane passing midway (rostro-caudal) through the decussation of the anterior and posterior part of the anterior commissure, and c) a horizontal plane passing midway through the most dorsal and ventral aspect of the decussation of the anterior commissure”. The stereotaxic orientation of the WHS is defined in the flat skull position with the lambda and bregma landmarks at the same height. The WHS reference data deviates 4 degrees from the flat skull position.
Terminology and structural hierarchy
The terminology used in the first version of the WHS rat brain atlas16 is based on the standard rat brain atlases from Paxinos and Watson (version 6) and Swanson (version 3)29,61. Subsequent revisions of the atlas incorporated domain specific nomenclatures for the annotations of the hippocampal formation and parahippocampal region62, and sought consensus among several sources for structures of the auditory system29,63,64. The terminology used for new structures in cortical, thalamic, striatopallidal, and midbrain regions (Supplementary File 1) is largely compatible with the terminology of Paxinos and Watson29,30, but have at several points been adapted to be more compatible with the terminology commonly used by domain experts. This particularly concerns the orbitofrontal, posterior parietal and insular cortex, which reflect the terminology used in seminal papers on the prefrontal cortex in several species65–67. For striatopallidal and midbrain dopaminergic regions, the terminology used by Paxinos and Watson has changed considerably across atlas versions, and our terminology is mostly consistent with their most recent MRI-based atlas37. For the subdivisions of the thalamus most of our terms are compatible with Paxinos and Watson29,30, with addition of some terms used by domain experts for structures in the posterior thalamus68, where the Paxinos and Watson atlas is less detailed. The hierarchical organization was based on the five secondary brain vesicles, with cortical areas, brain regions, and subregions for each vesicle sorted into a hierarchy of structures that largely followed the terminologies used in other atlases7,30,31, in particular the the Allen Mouse brain Common Coordinate Framework7. The hierarchy is provided as a structured .ilf file.
Annotation set and delineation procedure
All anatomical delineation (Fig. 5a3) are manually drawn based on a combination of features in the reference data, histological material, and information available in the literature. Below, we describe the general procedure of the delineation process, before giving a brief summary of the annotations from previously published versions of the atlas. Lastly, we summarize the annotations created for version 4 of the atlas.
General delineation procedure. Annotations in the WHS rat brain atlas were made using the ITK-SNAP software (version 3.6.0; http://itksnap.org39). Delineations were made in the principal orthogonalplanes (coronal, sagittal, and horizontal) using three approaches (Fig. 5b): 1) interpretation of sMRI/DTI contrast in the reference data; 2) inspection of histological reference material showing cyto-, chemo- or myeloarchitecture; and 3) consultation of literature, expert knowledge, and other brain atlases. These three approaches have been used for all previous versions of the atlas, as well as for the new version 4, and are briefly described below.
When interpreting sMRI/DTI contrast in the reference data, distinct contrast between coherent axonal fibre bundles and surrounding cell-rich areas forms a starting point for the delineation process by distinguishing between white and grey matter. The T2*-weighted sMRI greyscale maps show high grey to white matter contrast and highlights several cytoarchitectonic features. In the DTI maps, the prevailing orientation and magnitude of water diffusion in each DTI voxel is represented by the red, green, and blue colours, each signifying diffusion in one of the three principal directions (red = mediolateral, green = rostrocaudal, blue = dorsoventral). Intermediate diffusion orientations are indicated by intermediate colours according to the RGB model. The brightness of the colours represents the magnitude of oriented diffusion in each voxel, as determined by fractional anisotropy (FA) values69. Fluid-containing spaces appear white in T2*-weighted sMRI and dark in DTI, while air-filled regions are black in both modalities. Combined interpretation of sMRI and DTI maps viewed in coronal, sagittal, and horizontal planes allows identification of a large number of regional and subregional borders, aided by interactive inspection of surface rendering of structures to ensure smooth and 3D coherent surfaces.
At locations where the sMRI/DTI images is difficult to interpret, images of histological images showing myeloarchitecture and cytoarchitecture40,41,70 aid more detailed interpretations of the sMRI/DTI images, and refinement of border delineation. These public images have been spatially registered to the WHS rat brain atlas, and are available for inspection in an online image viewer, together with reference atlas overlay images40–42. Thus, observed features in the reference data can be interpreted in comparison with cyto- and myelin-stained histological sections, as previously described18.
Lastly, the delineation process is aided by consulting the literature, standard reference atlases29,61, and neuroanatomy experts. Key literature sources include studies of Sprague Dawley rats describing cyto-, chemo, or myelo-architectural characteristics of brain regions, with specific delineation criteria and clear descriptions of borders between regions71–73. For additional validation of the delineations, we used datasets containing 2D atlas plates from three reference atlases29,31,37 spatially registered to our volumetric reference data set74–76. This allowed us to directly compare spatial correspondences of anatomical landmarks and structural delineation across atlases in the standard planes provided by the 2D atlas plates.
Annotations from previously published versions. The annotation set in version 1 of the atlas included major anatomical structures that were identified from readily distinguishable differences in sMRI and DTI signals16. For version 2 of the atlas, detailed subdivisions in the hippocampus and parahippocampal region were identified based on observed sMRI/DTI contrast, corresponding to histologically defined borders17. For version 3 of the atlas, the complete ascending auditory system was delineated based on the interpretation of image features in the WHS rat brain sMRI/DTI data, and validated in relation to spatially corresponding images cell- and myelin-stained histological sections18.
New annotations in version 4. To create the new annotations, we used the following files from the atlas home page on NITRC (https://www.nitrc.org/projects/whs-sd-atlas/):
- WHS_SD_rat_T2star_v1.01.nii.gz; Image template showing anatomical MRI, T2*- weighted gradient echo image at 39 μm original resolution (1024x512x512 voxels).
- WHS_SD_rat_DWI_v1.01.nii.gz; Image template showing diffusion weighted image map (DWI), resampled to 39 µm resolution (1024x512x512 voxels).
- WHS_SD_rat_FA_color_v1.01.nii.gz; Image template of diffusion tensor (DTI) showing colour Fractional Anisotropy (FA) map resampled to 39 μm resolution (1024x512x512 voxels).
- WHS_SD_rat_atlas_v3.nii.gz; Volumetric atlas file containing 118 anatomical structures.
- WHS_SD_rat_atlas_v3.label; Label file specifying the ID, colour code, and name of each anatomical structure.
We delineated new regions across the brain, focusing on the cerebral cortex, striatopallidal region, thalamus, and midbrain dopaminergic regions, using the criteria outlined below.
Cerebral cortex.We added 36 new structures in the cerebral cortex. The cerebral cortex appears highly homogeneous in the sMRI/DTI data (Fig. 5), which has insufficient spatial resolution to resolve finer details of cortical layers and subtle cytoarchitectonic borders. We therefore used the cortical delineation of the 6th edition of the atlas by Paxinos and Watson29, the 4th edition of Swanson’s atlas31, and the MRI based atlas by Paxinos and colleagues37 as starting points. Diagrams from these atlases have been spatially registered to the WHS reference data (v1.01)74–76, allowing co-visualization and comparison of the different cortical delineation from these atlases to anatomical features visible in the sMRI/DTI data. This allowed us to identify the approximate stereotaxic location of cortical areas in coronal, sagittal, and horizontal levels in the sMRI/DTI data, and create the initial delineation. For the somatomotor and somatosensory areas, the locations and shapes of borders were adjusted to stereotaxic positions transferred from published maps defined by electrophysiological measurements33–36,77,78. For areas of the orbitofrontal, insular and posterior parietal cortex, the delineation were adapted in three planes on the basis of descriptions from previous anatomical studies79–81. Visual areas were adjusted in accordance with additional literature73,82, abutting the boundaries of the parahippocampal region17 (v2) and the auditory cortex18 (v3). In this way, the new delineations of cortical areas represent a composite of several sources, where most of the delineation are largely compatible with other reference atlases and published cortical maps. Adjacent to the cerebral cortex we defined three new structures, the claustrum, endopiriform nucleus, and a collective area called the amygdaloid area, unspecified, based on observed grey-white matter contrast.
Striatopallidal regions.We added 7 new and revised 3 existing structures in the striatopallidal regions. These regions were only coarsely delineated in previous versions. The annotation of the striatum from previous versions16 was subdivided into the caudate putamen as well as the core and shell of the nucleus accumbens. The dorsal boundary of the caudate-putamen complex is clearly demarcated by white matter, and core / shell region of the nucleus accumbens were defined by sMRI signal intensity contrast and their location relative to the anterior commissure. Furthermore, we delineated an area located directly caudal to the nucleus accumbens and ventral to the caudate putamen. This ventral striatal region can be identified by its sMRI signal intensity, and corresponds partly to regions referred to as the fundus of striatum and interstitial limb of the posterior part of the anterior commissure (IPAC) in other atlases29,31. However, we did not find adequate information in the template or in our reference data to subdivide the fundus of striatum or IPAC separately, and thus termed it as an ‘unspecified’ ventral striatal region. This area was included in the striatum annotation in previous versions of the WHS atlas. The ventral boundary of the ventral striatal region and shell of the nucleus accumbens towards the basal forebrain region was defined by using principal diffusion orientation differences visible in the DTI data. In the pallidum we subdivided the delineation of the external globus pallidus (referred to as globus pallidus in previous versions) into a medial and a lateral part. Furthermore, we delineated the ventral pallidum, which extends ventrally and rostrally from the external globus pallidus. The ventral pallidum is highly clustered and was identified by highly oriented (anisotropic) diffusion signal related to the anterioposteriorly oriented fibres from the olfactory tract and neighbouring basal forebrain areas. Lastly, in the pallidum, we also revised the delineation of the entopeduncular nucleus.
Thalamus.We added 55 new subregions in the thalamus, which was delineated as a single structure in previously published versions. Within the thalamus, white matter bands such as the medial lemniscus, the external and internal medullary lamina, and the superior cerebellar peduncle, are distinctly visible in the sMRI/DTI data and provide a framework for delineating subregions of the pre-, epi- and dorsal thalamus. In the prethalamus, we delineated the reticular (pre)thalamic nucleus (only the auditory segment was delineated as part of the previous version18, then referred to as the auditory segment of the reticular thalamic nucleus). We furthermore delineated rostral, dorsal, ventral and caudal parts of the zona incerta, the A13 and A11 dopamine cell groups, and the Fields of Forel.In the epithalamus, we delineated the habenula and subdivided it into lateral and medial parts. In the dorsal thalamus, we delineated 43 areas belonging to nine regional groups: the anterior nuclei; the dorsal-caudal midline group; the ventral midline group; the mediodorsal nucleus; the ventral nuclei; the intralaminar nuclei; the posterior complex; the lateral posterior (pulvinar) complex; the laterodorsal thalamic nuclei; and the medial geniculate complex. The medial geniculate complex was subdivided in the previously published version 3 of the atlas18, but all subregions except the marginal zone were revised and have therefore been assigned new region IDs in version 4. In this version of the atlas, the medial geniculate complex thus includes the dorsal, ventral, and medial divisions, as well as the marginal zone and the suprageniculate nucleus. In general, the delineation of thalamic subregions largely follows the annotations of Paxinos and Watson29,30. However, in the posterior thalamus, the atlases by Paxinos and Watson29,30 are incomplete, with large areas remaining unspecified. The posterior complex of the thalamus in the WHS rat brain atlas includes the posterior thalamic nucleus and the posterior thalamic nuclear group, triangular part.
Midbrain dopaminergic regions. We added 4 new structures to the midbrain dopaminergic regions. We delineated the subregions of the substantia nigra, i.e. the reticular, compact, and lateral parts. The ventral tegmental area was also delineated based on its appearance in the sMRI/DTI data and its relative position to the substantia nigra.
Quantification of parvalbumin neurons using the QUINT workflow with WHSv4
To demonstrate the practical value of the complete rat brain atlas for analysis, interpretation, visualization and communication of data, we re-analysed a published data set with the new atlas using the QUINT workflow. The raw and derived data used for this analysis24 are available through the EBRAINS Knowledge Graph, and were previously interpreted using WHS rat brain atlas v2. We here re-analysed these data for one of the subjects (rat 25205), using the derived data set50 and version 4 of the atlas. To optimize the spatial registration to WHSv4, we adjusted the .json file (ext-d000008_PVRat_25205_nonlinear.json) provided with the data set42 using VisuAlign (v0.9; RRID:SCR_017978). We then exported the new customized atlas maps and combined these with the existing segmentation maps using Nutil Quantifier. We used the same Nutil parameters as previously24 except that we created new custom regions tailored to WHS rat brain atlas v4. The custom regions largely correspond to the highest level of detail in the atlas, but all white matter tracts were merged, as well as very fine grey matter areas such as cortical layers. Upon running Nutil Quantifier, we post-processed the reports as described previously24 to arrive at total number and density estimates for each custom brain region.