PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Neuroimage. Author manuscript; available in PMC Jan 1, 2012.
Published in final edited form as:
PMCID: PMC2962762
NIHMSID: NIHMS232798
An MRI-based Atlas and Database of the Developing Mouse Brain
Nelson Chuang,1 Susumu Mori,1,5 Akira Yamamoto,1 Hangyi Jiang,1 Xin Ye,2 Xin Xu,1 Linda J. Richards,6,7 Jeremy Nathans,2,4 Michael I. Miller,3,8 Arthur W.Toga,9 Richard L. Sidman,10 and Jiangyang Zhang1
1Department of Radiology, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
2Department of Molecular Biology and Genetics, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
3Department of Biomedical Engineering, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
4Howard Hughes Medical Institution, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
5F.M. Kirby Functional Imaging Center, Kennedy Krieger Institute, Baltimore, Maryland, USA
6Department of Anatomy and Neurobiology, University of Maryland School of Medicine, Baltimore, Maryland, USA
7The Queensland Brain Institute and the School of Biomedical Sciences, University of Queensland, Brisbane, Australia
8Center of Imaging Science, Johns Hopkins University, Baltimore, Maryland, USA
9Laboratory of Neuro Imaging, Department of Neurology, University of California, Los Angeles, California, USA
10Department of Neurology, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, Massachusetts, USA
Correspondence to: Jiangyang Zhang, Ph.D. Department of Radiology, Johns Hopkins University School of Medicine, 217 Traylor Bldg, 720 Rutland Ave., Baltimore, MD, 21205, USA ; jzhang3/at/jhmi.edu Phone: (410)-502-9856
The advent of mammalian gene engineering and genetically modified mouse models has led to renewed interest in developing resources for referencing and quantitative analysis of mouse brain anatomy. In this study, we used diffusion tensor imaging (DTI) for quantitative characterization of anatomical phenotypes in the developing mouse brain. As an anatomical reference for neuroscience research using mouse models, this paper presents DTI based atlases of ex vivo C57BL/6 mouse brains at several developmental stages. The atlas complements existing histology and MRI-based atlases by providing users access to three-dimensional, high-resolution images of the developing mouse brain, with distinct tissue contrasts and segmentations of major gray matter and white matter structures. The usefulness of the atlas and database was demonstrated by quantitative measurements of the development of major gray matter and white matter structures. Population average images of the mouse brain at several postnatal stages were created using large deformation diffeomorphic metric mapping and their anatomical variations were quantitatively characterized. The atlas and database enhance our ability to examine the neuroanatomy in normal or genetically engineered mouse strains and mouse models of neurological diseases.
The mammalian brain is one of the most complicated structures within living organisms. The three-dimensional organization (comprising the spatial dimensions), which consists of more than 100 billion neurons in the adult brain, combined with the developmental processes (comprising a temporal dimension), constitute a vast four-dimensional (4D) entity. Although there is much literature that describes the anatomy of the rodent brain and its development, for example, the development of the cortex (Bayer and Altman, 1991; Sidman and Rakic, 1973), thalamus (Altman and Bayer, 1988), and cerebellum (Altman and Bayer, 1985), the spatiotemporal anatomy of the developing rodent brain remains poorly documented. The need for such information has become even more urgent because of increasingly sophisticated gene and cell engineering technologies. For example, mapping gene and protein expression, detecting anatomical changes caused by gene mutations, and delivery and extraction of cells in the brain all require precise 4D anatomical information of the rodent brain. Existing atlases contain only limited coverage of developmental stages. Most of these atlases are histology-based and two-dimensional, and therefore, have limited slice locations and orientations. As a result, several key aspects of the spatiotemporal development of the mouse brain have not been well documented. Questions about the volume and coordinates of specific brain structures in the developing C57BL/6 mice, and the range of their anatomical variation can only be answered with a three-dimensional atlas of the developing mouse brain.
With modern magnetic resonance imaging (MRI), we can now generate three-dimensional images of the mouse brain with high spatial resolution (Benveniste et al., 2000; Jacobs et al., 1999; Johnson et al., 2002; Kovacevic et al., 2005; Natt et al., 2002; Petiet et al., 2008). This opens up new possibilities for establishing 4D atlases of the developing mouse brain. The advantages of the MRI-based atlases are clear. The image is inherently three-dimensional (3D) without distortions due to embedding and sectioning, and the entire mouse brain can be imaged with a one-to-one correspondence between image pixels and tissue regions. The data are collected and stored in digitized formats, which makes operations such as rotation, segmentation, quantification, and standardization straightforward. The accuracy (less distortion) and the time efficiency (the entire 3D image reconstruction requires only hours) also allow investigations of anatomical variation among individuals and the creation of a probabilistic atlas. MRI-based atlases, however, have a limited resolution (up to 20 ~ 30 μm) and a limited range of image contrasts, meaning many brain structures, such as various nuclei in the thalamus, are difficult to delineate. Nonetheless, MRI-based atlases carry valuable information that is complementary to histology-based atlases and could be an important tool for future brain research.
Several MRI-based atlases have been established (Aggarwal et al., 2009; Badea et al., 2007; Dorr et al., 2008; Jacobs et al., 1999; Ma et al., 2005; MacKenzie-Graham et al., 2004; Petiet et al., 2008). Most recently, Chan et al. introduced a combined MRI-CT 3D atlas of the adult mouse brain (Chan et al., 2007), demonstrating exciting opportunities for three-dimensional navigation of brain anatomy. Our group recently also published a stereotaxic mouse brain atlas that combines CT and diffusion tensor imaging (DTI) data (Aggarwal et al., 2009). Promising applications for MRI-based atlases, such as registration of histology-based information (e.g., gene expression maps), have also been demonstrated (MacKenzie-Graham et al., 2004).
The atlas described here is an extension of these previous studies with several important advancements. It was created from a database of 3D high resolution T2-weighted and diffusion tensor images of the developing C57BL/6 mouse brains, beginning at embryonic day 12 (E12). The inclusion of DTI data, of which the contrasts mainly reflect tissue microstructural organization (Basser and Pierpaoli, 1996; Le Bihan, 2003; Mori and Zhang, 2006), are key for delineation of various white matter structures and certain gray matter structures in the developing mouse brains (Mori et al., 2001; Zhang et al., 2003). Based on this database, major gray and white matter structures were identified, labeled, and segmented in MR images based on existing histology-based atlases. Quantitative characterization of normal brain development, including the sizes of the whole brain, neocortex, cerebellum, hippocampus, and several white matter tracts, was performed. Finally, we characterized anatomical variability at several developmental stages. The atlas and database will be made publicly available at http://lbam.med.jhmi.edu to facilitate research using this 4D mouse brain MRI database.
Experimental animals
All experiments and procedures were approved by the Animal Research Committee of the Johns Hopkins University School of Medicine and the University of Maryland School Of Medicine. C57BL/6 mice (The Jackson Laboratory, Bar Harbor, ME, USA) at various developmental stages (Table 1) were perfusion-fixed with 4% paraformaldehyde (PFA) in 0.1 M phosphate buffer saline (PBS) either by immersion fixation for embryonic and neonatal samples, or by trans-cardiac perfusion fixation for postnatal and adult (P60 and older) mice. Mouse heads were removed and immersed in 4% PFA for more than 48 hours. To preserve the shape of the brains during imaging, the brains were left inside the skull. Before imaging, the samples were placed in PB with 0.1 mM gadopentetate dimeglumine (Magnevist, Berlex Imaging, Wayne, NJ, USA) for more than 48 hours. The addition of gadolinium enhanced MR signals and reduced imaging time (Johnson et al., 2002). We then transferred the samples into custom-built, MR-compatible tubes. The tubes were filled with fomblin (Fomblin Profludropolyether, Ausimont, Thorofare, New Jersey, USA), a non-hydrogen containing liquid for susceptibility matching and to prevent dehydration.
Table 1
Table 1
Number of animals and imaging resolution for data included in the developing mouse brain atlas.
Data Acquisition
Ex vivo imaging was performed on an 11.7 Tesla NMR spectrometer (Bruker Biospin, Billerica, MA, USA). Bird-cage RF coils with diameters of 10 ~ 20 mm were used for transmission and reception. A 3D diffusion-weighted multiple spin echo sequence was used to acquire images with four spin echoes (six spin echoes for embryonic and neonatal brains), a repetition time (TR) of 0.7s, an echo time (TE) of 35 ms, and two signal averages. The field of view ranged from 9 mm × 6 mm × 6 mm for the smallest sample (E12) to 20 mm × 14 mm × 10 mm for the adult samples. The native imaging resolution ranged from 0.08 mm × 0.08 mm × 0.08 mm for embryonic samples to 0.125 mm × 0.125 mm × 0.125 mm for adult samples. Images reconstructed from each echo were added together to form one image to enhance the signal-to-noise-ratio (SNR). For each sample, eight diffusion-weighted images were acquired with a diffusion gradient duration (δ) of 6 ms and a separation between the pair of diffusion gradients (Δ) of 14 ms for postnatal and adult brains and δ = 5 ms / Δ = 11 ms for embryonic and neonatal brains (P7 and younger). Two images were acquired with a minimum b value (50 s/mm2) and six images with a maximum b value (1,000~1,200 s/mm2 for postnatal and adult samples and 700 ~ 800 s/mm2 for embryonic and neonatal samples). Diffusion sensitizing gradients were applied along six different orientations: [0.707, 0.707, 0], [0.707, 0, 0.707], [0, 0.707, 0.707], [−0.707, 0.707, 0], [0.707, 0, −0.707], and [0, −0.707, 0.707]. The total imaging time was approximately 24 hours. Co-registered 3D T2-weighted images were acquired with the same field of view and resolution using a fast spin echo sequence, with a TE/TR of 50/900ms, a flip angle of 40 degrees, four signal averages, an echo train length of four (the first echo was located at the center of the k-space), and an imaging time of less than two hours.
Image processing
The ex vivo 3D T2-weighted and diffusion-weighted images were reconstructed on an off-line PC workstation using IDL (ITT Visual Information Solutions, Boulder, CO, USA). Diffusion tensor fitting was performed using DtiStudio (Jiang et al., 2006). The six elements of diffusion tensor were determined by log-linear fitting. The tensor was diagonalized to obtain three eigenvalues (λ1–3) and corresponding eigenvectors (ν1–3). Fractional anisotropy (FA) (Basser and Pierpaoli, 1996) was calculated.
equation M1
Direction-encoded color map (DEC) images were generated by combining the images of the primary eigenvector (V1) and FA into RGB images (V1·FA). In DEC images, the ratio among the G(reen), R(ed), and B(lue) components of each pixel was defined by the ratio of the absolute values of the x, y, and z components of the primary eigenvector, and the intensity was proportional to the FA. In our atlas, red was assigned to the medial-lateral axis, green to the rostral-caudal axis, and blue to the dorsal-ventral axis. Average diffusion-weighted (aDW) images are the average of all diffusion-weighted images from each specimen. Signals from skull tissues in these images were manually removed and the total brain volume of each specimen was measured. At each developmental stage, a representative brain with the median brain volume was selected as the template and manually re-oriented to the standard orientations as defined in two histology-based atlases ((Foster, 1998) for embryonic brains and (Paxinos and Franklin, 2003) for postnatal and adult brains). In the re-oriented images, the left-right and rostral-caudal axes of the brain were along the x and y axes, respectively. For postnatal and adult brains, the dorsal, flat part of the cortex was perpendicular to the z axis. For embryonic brains, the line that connected the bases of the hypothalamus and brainstem in the mid-sagittal plane was perpendicular to the z axis. The aDW images from other specimens at the same developmental stage were rigidly aligned to the aDW image of template using the six-parameter rigid registration function in the automated image registration (AIR) package (Woods et al., 1998a; Woods et al., 1998b). The derived rigid transformation was then applied to the T2-weighted and tensor data. The orientation of each tensor was aligned using the methods described by Xu et al. (Xu et al., 2003).
Segmentation, volume, and cross-sectional area measurements and fiber tract reconstruction
Major brain compartments, e.g., neocortex, hippocampus, and cerebellum, were manually segmented in the rigidly aligned MR data using Amira (Visage Imaging, Inc. San Diego, CA, USA), and their volumes were calculated based on the segmentation results. We followed the definition of compartmental boundaries in the Paxinos' atlases (Paxinos and Franklin, 2003; Paxinos et al., 2006). The cross-sectional area of major white matter tracts was measured at selected locations along each white matter tract in the DEC images. Using DtiStudio, we first located the tract in the DEC image, then rotated the image so that the viewing plane was orthogonal to the tract, and finally defined the area by tracing the boundary of the tract in the rotated DEC image. This operation was repeated three times for each tract, and the average area was saved as the cross-sectional area of the tract. Three-dimensional tract reconstruction was performed using the FACT method (DtiStudio) (Mori et al., 1999; Xue et al., 1999) with an FA threshold of 0.2 and fiber angles less than 40° between two connected pixels. A multi-ROI approach was used to reconstruct tracts of interest (Conturo et al., 1999; Huang et al., 2004).
Quantification of anatomical structures during development
Changes in the volumes of major brain compartments and cross-sectional areas of major white matter tracts during development were fitted to a sigmoidal model, originally proposed by Wingert (Wingert, 1969).
equation M2
The sigmoidal model characterizes the growth of brain structures using three parameters: a, x0, and b. The first parameter, a, represents the estimated volume or cross-sectional area of a structure in the adult mouse brain. The second parameter, x0, is the time to reach half the estimated adult volume (a). The third parameter, b, describes the growth curve; the smaller the b value, the faster the growth. The variable x is the gestational age (days since conception). Volume data was fitted to the above sigmoidal model using nonlinear regression (SigmaPlot 10.0, Systat Software, San Jose, CA, USA) with a threshold for significance of p < 0.0001.
Generation of population averaged brains and individual variation maps
Ex vivo images acquired from mouse brains at 7 postnatal stages (P7, P11, P15, P21, P28, P42, and P60, n=6 at each stage) were used to generate population averaged mouse brain images at each stage. In order to minimize differences in aDW image intensity caused by different receiver gain values used during acquisition, the intensity values of the rigidly aligned aDW images from each specimen were adjusted using a multiplicative scaling factor to the representative brain with median brain volume at each stage. The multiplicative scaling factor was obtained by comparing the intensity values of several gray matter structures (the cortex, thalamus, and striatum) and white matter structures (the genu and splenium of the corpus callosum, anterior commissure, and cerebral peduncle). The intensity values of the aDW images were then scaled uniformly to a 0 to 255 range, while the intensity values in the FA images, which range from 0 to 1, were multiplied by 255. For generation of population averaged mouse brain images, the iterative procedure previously described in (Kovacevic et al., 2005) was adopted. Image registration was performed using dual contrast (aDW + FA) Large Deformation Diffeomorphic Metric Mapping (LDDMM) with equal weighting for aDW and FA signals, which can accurately register images from two subjects based on the tissue contrasts in the aDW and FA images (Ceritoglu et al., 2009). The transformation of the tensor field was carried out by the method described in (Xu et al., 2003). The results included population averaged T2-weighted, diffusion weighted, and diffusion tensor data and the deformation maps between individual mouse brain images and the population averaged images, from which the covariance tensor at each pixel, x0, in the population averaged images was computed as
equation M3
where Dcov is the covariance tensor, n is the number of subjects, xk is the coordinates in the image of subject k that corresponds to pixel x0 in the population averaged image according to the deformation map, and equation M4 is the average of all xk, k = 1, 2, …, n. The trace value of the covariance tensor characterizes the magnitude of anatomical variation at each pixel.
The MRI / DTI-based atlas of the developing mouse brain
Figure 1 shows coronal T2-weighted, FA, and direction-encoded color-map (DEC) images of an E16, a P0 (postnatal day 0), and an adult (P60) mouse brains, overlaid with segmented gray and white matter structures. The atlas focuses on the major white matter tracts in the mouse brains at embryonic, postnatal and adult stages. Figure 2 and figure 3 show three-dimensional views of the segmented structures in the E16, P0, and adult mouse brains. The spatial relationships among various neural structures and anatomical changes during development are clearly demonstrated. White matter structures, with their intermediate to high fractional anisotropy (FA) and unique orientations (as distinctive colors in the DEC image) dictated by the trajectory of each tract, can be identified in the embryonic and neonatal mouse brains. For example, the trigeminal nerve (5n), vestibulococlear nerve (8n), optic nerve (2n), and optic tract are already visible in the images of an E16 mouse brain (Fig. 2A). The three major commissural structures, the anterior commissure, the hippocampal commissure, and the corpus callosum, are barely perceptible in our MR data at E16 (not shown), but become visible by P0 in the mid-sagittal plane (Fig. 2B). In the thalamic region, the cerebral peduncles, the stria terminalis, and the fasciculus retroflexus are visible at E16 (Fig. 3B). Tracts related to the cerebellum are difficult to detect in our MR data at E16. For example, the middle cerebellar peduncle first becomes visible at P0 (Fig. 3A). Examination of the cross-sectional areas of white matter tracts at comparable locations during development suggests that they have different temporal profiles of growth. For example, the cross-sectional area of the fasciculus retroflexus (fr inFig. 3B) only shows a small increase from P0 to adult, while the cerebral peduncle (cp in Fig. 3A), middle cerebellar peduncle (mcp in Fig. 3A), and corpus callosum (indicated by white arrows in Fig. 2B) show dramatic growth during the same period.
Figure 1
Figure 1
Single subject coronal images of the E16, P0, and adult mouse brains from the mouse brain database. At each stage, the brain specimen with the median total brain volume was shown. The images are overlaid with reconstructed white and gray matter structures (more ...)
Figure 2
Figure 2
Development of the cranial nerves and commissural tracts in the mouse brain. At each stage,single subject images from the brain specimen with the median total brain volume was shown. A) The trigeminal nerve (5n, green), facial nerve (7n, yellow), vestibulocochlear (more ...)
Figure 3
Figure 3
Development of the white matter tracts related to the cerebellum and thalamus. At each stage,single subject images from the brain specimen with the median total brain volume was shown. A) The middle cerebellar peduncle (mcp, green), the internal capsule (more ...)
Volumetric Development
In Figure 4, the volume of the whole brain, neocortex, hippocampus, and cerebellum were examined from E14 (gestational day 14) to adult. We observed sigmoidal growth patterns in our MR-based volume measurements and the estimated parameters are summarized in Table 2. These measurements establish the normal growth curve of C57BL/6 mice and the normal range of individual variability measured by ex vivo MRI (via the 95% prediction bands). Close inspection of the growth of individual structures reveals delayed (greater x0) growth of the cerebellum, which coincides with the known delayed anatomical development of the cerebellum (Table 2).
Figure 4
Figure 4
Volumes of the whole brain, neocortex, hippocampus, and cerebellum (left vertical axis) and their percentage of adult volumes (right vertical axis). The volume measurements were fitted to a sigmoidal model with 95% confidence band (long dash lines) and (more ...)
Table 2
Table 2
Estimated parameters and standard errors of the sigmoidal model that characterize the volumetric development of the whole brain, neocortex, cerebellum, and hippocampus in the C57BL/6 mouse. The parameters for the C57BL/6 mice were estimated, based on (more ...)
White Matter Fiber Tract Cross-Sectional Areas
Figure 5 shows the changes in cross-sectional areas from birth to adult in twelve white matter tracts. The measurements were fitted to the sigmoidal model, and the estimated parameters and standard errors for each tract are tabulated in Table 3. At birth (18.5 days after conception, P0), the cross-sectional areas of most tracts were less than 40% of their estimated adult values (a), except the hippocampal commissure, whose cross-sectional area at P0 was approximately 60% of its estimated adult value (a = 21.11 ± 0.36 mm2). The estimated x0s (time to reach 50% of the estimated adult values), based on the cross-sectional areas, suggested that these white matter tracts can be classified into three groups: an early group (x0 < 18.5 days (P0), before birth) that includes the hipppocampal commissure (hc); a middle group (x0 between 18.5 days (P0) and 26.5 days (P8), the first week after birth) that includes the anterior commissure (ac), corpus callosum (cp), fornix (fx), stria medularis (sm), fasciculus retroflexus (fr), and the trigeminal nerve (5n); and a late group (x0 larger than 26.5 days or P8) for the rest. In terms of the growth rate of the cross-sectional area, the fasciculus retroflexus showed the highest growth rate (the smallest b).
Figure 5
Figure 5
Changes in white matter cross-sectional areas in twelve major white matter tracts. The area data were fitted to a sigmoidal model and the 95% confidence intervals are shown in the figures, with figure background indicative of the percentile of the estimated (more ...)
Table 3
Table 3
Estimated parameters and standard erros of the sigmoidal model that characterize the postnatal development of white matter tract cross-sectional areas in the C57BL/6 mouse brains. Structural abbreviations are: ac: anterior commissure; hc: hippocampal (more ...)
Population-Averaged Mouse Brain Atlas
Population averaged T2-weighted and Direction-encoded colormap (DEC) images of the mouse brain at several stages are shown in Figure 6 with segmentation of major gray and white matter structures. The average mouse brain images have high signal-to-noise ratio and sharply defined boundaries of major white matter tracts. Furthermore, with the transformations generated by LDDMM, the morphological variation in adult mouse brains was quantitatively characterized using a covariance tensor at each pixel. Figure 6 shows the amount of anatomical variation in the P7, P21, and P60 mouse brain with respect to the population averaged mouse brain images using the trace of the covariance tensor (Fig. 6, Variation). The distribution of anatomical variation in the brain is not uniform and the amount of variation varies at different developmental stages. Figure 6 also illustrates the amounts of variation on the surfaces of the cortex, hippocampus, lateral ventricles, and striatum in the adult mouse brain. The frontal and posterior dorsal regions of the cortex display variation at the level of 0.1 mm. The anterior part of the hippocampus, which is adjacent to the lateral ventricles, displayed a high degree of variation (~ 0.3 mm). The striatum has a relatively low degree of variation.
Figure 6
Figure 6
Top: Population average T2, directional encoded colormap (DEC) images of the P7, P21, and P60 mouse brain and maps of segmentation and maps of anatomical variation overlaid with outlines of structural segmenations. The unit in the variation maps is mm (more ...)
Compared to existing histology or MR-based atlases, the atlas presented here is unique in its characterization of developing mouse brains using high resolution DTI technique. DTI has become a useful phenotyping tool in neuroscience research, especially at early developmental stages. While conventional T1 and T2 based MRI show poor tissue contrast in immature mouse brains, probably due to lack of myelin, DTI was able to provide satisfactory and relatively consistent tissue contrasts in the developing mouse brains. Ex vivo mouse brain DTI has been used to examine the white matter structures and their development in wild-type and mutant strains (Andrews et al., 2006; Tyszka et al., 2006; Verma et al., 2005; Wang et al., 2006; Zhang et al., 2005; Zhang et al., 2003). With these developments, this atlas will be a useful guide for locating anatomical structures in mouse DTI data, and our database may serve as a baseline for the study of abnormal mouse brain development.
The information in the atlas is presented in 2D and 3D formats, and the atlas is simple to use. The 2D atlas contains images from developing brains with structural labels, and can be used in the same ways as current histology-based atlases. We followed the structural nomenclature used in the two atlases as closely as possible. The 3D atlas provides visualization of the segmented white and gray matter structures of interest, thus allowing exploration of the 3D spatial relationships between these structures, e.g., the 3D trajectories of a particular white matter tract. By coregistering the atlas and MRI/DTI data to CT data for the external skull features, the 3D atlas may also be useful for surgical planning among other interventions (Aggarwal et al., 2009). It is necessary to note that even though the images in the atlas have been rotated to the orientation defined in the Paxinos atlas (Paxinos and Franklin, 2003) and Foster atlas (Foster, 1998) for postnatal and embryonic mouse brains, respectively, small discrepancies in slice orientation may still present because accurate matching of our 3D MRI based atlases to the 2D section based atlases is difficult.
To avoid tissue deformation or damage during dissection and imaging, the images were acquired with the brains preserved in the skull. The disadvantages associated with imaging brains within the skull are in-skull susceptibility artifacts near the ear drum, and reduced SNR and resolution due to the use of larger, less sensitive RF coils. Y. Ma et al. (Ma et al., 2005) reported significant differences in brain shapes between images obtained from brain specimens preserved inside and outside the skull, which was attributed to tissue relaxation and stretching after dissection. The volume measurements reported here are comparable with data from previous reports based on histology and/or MRI. A previous histology-based study on albino mice showed that the volumetric growth of the brain, neocortex, and cerebellum follow a sigmoidal model (Wingert, 1969), with relatively slow growth in the late embryonic period (E16 ~ E18), followed by a rapid increase during the neonatal and early postnatal period (P0 ~ P20), subsequently slowly approaching their adult values. The data showed that the development of different brain compartments is not uniform. For example, while the volumetric growth of the neocortex and hippocampus followed the volumetric growth of the whole brain closely, the growth of the cerebellum is slower than the growth of the whole brain (Fig. 4). The estimated x0 for the cerebellum is around P8 to P9, compared to an estimated x0s of P5 to P6 for the whole brain, neocortex, and hippocampus. The sigmoidal model, however, is a general model used here to characterize the brain growth throughout the embryonic and postnatal periods. For more detailed analysis on volume development, it is necessary to focus on a particular period with more sophisticated models (Uylings et al., 1990) than the model adopted here.
While volume changes are prevalent in various types of human brain diseases, e.g., Huntington's disease (Nopoulos et al., 2007; Rosas et al., 2004) and Alzheimer's disease (Csernansky et al., 2005; Thompson et al., 2003) brain morphology has rarely been used as a phenotypic marker in mouse disease models. With the increasing availability of small animal imaging instruments, this situation is rapidly changing. The normal growth curves of C57BL/6 with 95% prediction bands will be a useful reference for future morphological phenotype studies.
In this study, population-averaged mouse brain images at several postnatal stages were generated by normalizing individual mouse brain image to a template using linear (rigid) and non-linear (LDDMM) transformations. After the transformations, we could obtain information about anatomical variability among individual mice. Recently, several groups have reported measurements of anatomical variation in mouse brains based on segmentation of MRI data (Badea et al., 2007; Dorr et al., 2008; Kovacevic et al., 2005; Ma et al., 2005). In this study, anatomical variation in adult mouse brains was characterized without segmentation. Instead, we relied on LDDMM, driven by the rich anatomical information in DTI data, to generate measurements of anatomical variation for each pixel. Figure 6 demonstrates that certain brain regions have greater variation than other regions. In addition, if the atlases are used for stereotactic operations, the variability data provide a level of confidence for the delivery of cells, chemicals, or electrodes, depending on the experiment. This type of information is difficult to obtain from histology-based morphology studies.
In this study, the maturation process of various white matter tracts was also quantitatively characterized in terms of their sizes. DTI can monitor only the macroscopic sizes of axonal bundles, which by no means provides precise cellular information, such as axonal counts, axonal calibers, and myelination. However, macroscopic axonal degeneration or misguidance is prevalent in various types of mutants, which are sometimes difficult to comprehend from serial histology sections. We believe that our atlases and quantitative growth curves of each white matter tract provide anatomical information in 0.1 – 1 mm scales, that is complementary to histological information (which is typically on a μm scale) and will serve as important references for axonal tract development studies.
Currently, the main limitation of our atlas is a lack of precise gray matter assignment based on correlations with histology. Although it is rather straightforward to identify major white matter tracts based on their trajectory information provided by DTI, assignment of gray matter structures is not always straightforward. MR-histology correlation studies will be an important future effort to establish more comprehensive atlases. Other limitations of our atlas are mainly DTI related. First, DTI is inherently a technique with low-signal-to-noise ratio, and therefore, it is more difficult to achieve a high spatial resolution image and is more time-consuming than conventional MRI and other 3D imaging modalities. The resolution achieved in this atlas ranges from 80 μm to 125 μm isotropic resolution. Previous reports based on T1- and T2/T2*-weighted MRI were able to reach up to 20 ~ 30 μm resolution (Badea et al., 2007; Dorr et al., 2008; Ma et al., 2005), while episcopic and microCT can achieve even higher resolution (Johnson et al., 2006; Weninger and Mohun, 2002). Second, DTI is based on a simplified model of water diffusion. When the anatomical structures are complex, for example in regions where two or more populations of fibers cross each other, DTI oversimplifies the underlying anatomy. Contrasts for such problematic areas could be improved by more sophisticated models, such as diffusion spectrum imaging, Q-ball imaging, and HARD (Frank, 2001; Tuch et al., 2002; Tuch et al., 2003; Wedeen et al., 2000).
In summary, an MRI-based atlas and a database of the developing mouse brain have been established. The atlas and database will be useful for the quantitative examination of brain phenotypes and their development in wildtype and mutant mouse strains. The atlas and database can also be used to study abnormal anatomy and lesions in mouse models of neurological diseases, and as a dataset for the development of novel analysis tools. In the future, we will continue to collect high-resolution data to add to our current database and pursue correlations between information provided by histology and MRI.
Research Highlights
  • Mouse brain development studied using MRI/DTI
  • MRI/DTI atlas and database of developing mouse brain
  • Characterization of anatomical changes in mouse brain development
Acknowledgments
This study was supported by National Institutes of Health AG20012, EB003543, NS059529, and ES012665. Support was also provided, in part, by P41 RR013642 and U24 RR021760 MouseBIRN to A. W. Toga and NIH R21 NS04584 to S. Mori, L. J. Richards, and Frank L. Margolis (University of Maryland).
Footnotes
This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
  • Aggarwal M, Zhang J, Miller MI, Sidman RL, Mori S. Magnetic resonance imaging and micro-computed tomography combined atlas of developing and adult mouse brains for stereotaxic surgery. Neuroscience. 2009;162(4):1339–50. [PMC free article] [PubMed]
  • Altman J, Bayer SA. Embryonic development of the rat cerebellum. I. Delineation of the cerebellar primordium and early cell movements. J Comp Neurol. 1985;231:1–26. [PubMed]
  • Altman J, Bayer SA. Development of the rat thalamus: I. Mosaic organization of the thalamic neuroepithelium. J Comp Neurol. 1988;275:346–377. [PubMed]
  • Andrews W, Liapi A, Plachez C, Camurri L, Zhang J, Mori S, Murakami F, Parnavelas JG, Sundaresan V, Richards LJ. Robo1 regulates the development of major axon tracts and interneuron migration in the forebrain. Development. 2006;133:2243–2252. [PubMed]
  • Badea A, Ali-Sharief AA, Johnson GA. Morphometric analysis of the C57BL/6J mouse brain. Neuroimage. 2007;37(3):683–93. [PMC free article] [PubMed]
  • Basser P, Pierpaoli C. Microstructural and physiological features of tissues elucidated by quantitative-dfifusion-tensor MRI. J Magn Reson B. 1996;111:209–219. [PubMed]
  • Bayer SA, Altman J. Neocortical Development. Raven Press; New York: 1991.
  • Benveniste H, Kim K, Zhang L, Johnson GA. Magnetic resonance microscopy of the C57BL mouse brain. Neuroimage. 2000;11:601–611. [PubMed]
  • Ceritoglu C, Oishi K, Li X, Chou MC, Younes L, Albert M, Lyketsos C, van Zijl PC, Miller MI, Mori S. Multi-contrast large deformation diffeomorphic metric mapping for diffusion tensor imaging. Neuroimage. 2009;47:618–627. [PMC free article] [PubMed]
  • Chan E, Kovacevic N, Ho SK, Henkelman RM, Henderson JT. Development of a high resolution three-dimensional surgical atlas of the murine head for strains 129S1/SvImJ and C57Bl/6J using magnetic resonance imaging and microcomputed tomography. Neuroscience. 2007;144:604–615. [PubMed]
  • Conturo TE, Lori NF, Cull TS, Akbudak E, Snyder AZ, Shimony JS, McKinstry RC, Burton H, Raichle ME. Tracking neuronal fiber pathways in the living human brain. Proc. Natl. Acad. Sci. USA. 1999;96:10422–10427. [PubMed]
  • Csernansky JG, Wang L, Swank J, Miller JP, Gado M, McKeel D, Miller MI, Morris JC. Preclinical detection of Alzheimer's disease: hippocampal shape and volume predict dementia onset in the elderly. Neuroimage. 2005;25:783–792. [PubMed]
  • Dorr AE, Lerch JP, Spring S, Kabani N, Henkelman RM. High resolution three-dimensional brain atlas using an average magnetic resonance image of 40 adult C57Bl/6J mice. Neuroimage. 2008;42:60–69. [PubMed]
  • Foster GA. Chemical Neuroanatomy of the Prenatal Rat Brain: A Developmental Atlas. 1st ed. Oxford University Press; Oxford: 1998.
  • Frank LR. Anisotropy in high angular resolution diffusion-weighted MRI. Magn Reson Med. 2001;45:935–939. [PubMed]
  • Huang H, Zhang J, van Zijl PC, Mori S. Analysis of noise effects on DTI-based tractography using the brute-force and multi-ROI approach. Magn Reson Med. 2004;52:559–565. [PubMed]
  • Jacobs RE, Ahrens ET, Dickinson ME, Laidlaw D. Towards a microMRI atlas of mouse development. Comput Med Imaging Graph. 1999;23:15–24. [PubMed]
  • Jiang H, van Zijl PC, Kim J, Pearlson GD, Mori S. DtiStudio: resource program for diffusion tensor computation and fiber bundle tracking. Comput Methods Programs Biomed. 2006;81:106–116. [PubMed]
  • Johnson GA, Cofer GP, Gewalt SL, Hedlund LW. Morphologic phenotyping with MR microscopy: the visible mouse. Radiology. 2002;222:789–793. [PubMed]
  • Johnson JT, Hansen MS, Wu I, Healy LJ, Johnson CR, Jones GM, Capecchi MR, Keller C. Virtual histology of transgenic mouse embryos for high-throughput phenotyping. PLoS Genet. 2006;2:e61. [PMC free article] [PubMed]
  • Kovacevic N, Henderson JT, Chan E, Lifshitz N, Bishop J, Evans AC, Henkelman RM, Chen XJ. A three-dimensional MRI atlas of the mouse brain with estimates of the average and variability. Cereb Cortex. 2005;15:639–645. [PubMed]
  • Le Bihan D. Looking into the functional architecture of the brain with diffusion MRI. Nat Rev Neurosci. 2003;4:469–480. [PubMed]
  • Ma Y, Hof PR, Grant SC, Blackband SJ, Bennett R, Slatest L, McGuigan MD, Benveniste H. A three-dimensional digital atlas database of the adult C57BL/6J mouse brain by magnetic resonance microscopy. Neuroscience. 2005;135:1203–1215. [PubMed]
  • MacKenzie-Graham A, Lee EF, Dinov ID, Bota M, Shattuck DW, Ruffins S, Yuan H, Konstantinidis F, Pitiot A, Ding Y, Hu G, Jacobs RE, Toga AW. A multimodal, multidimensional atlas of the C57BL/6J mouse brain. J Anat. 2004;204:93–102. [PubMed]
  • Mori S, Crain BJ, Chacko VP, van Zijl PCM. Three dimensional tracking of axonal projections in the brain by magnetic resonance imaging. Annal. Neurol. 1999;45:265–269. [PubMed]
  • Mori S, Itoh R, Zhang J, Kaufmann WE, van Zijl PCM, Solaiyappan M, Yarowsky P. Diffusion tensor imaging of the developing mouse brain. Magn Reson Med. 2001;46:18–23. [PubMed]
  • Mori S, Zhang J. Principles of diffusion tensor imaging and its applications to basic neuroscience research. Neuron. 2006;51:527–539. [PubMed]
  • Natt O, Watanabe T, Boretius S, Radulovic J, Frahm J, Michaelis T. High-resolution 3D MRI of mouse brain reveals small cerebral structures in vivo. J Neurosci Methods. 2002;120:203–209. [PubMed]
  • Nopoulos P, Magnotta VA, Mikos A, Paulson H, Andreasen NC, Paulsen JS. Morphology of the cerebral cortex in preclinical Huntington's disease. Am J Psychiatry. 2007;164:1428–1434. [PubMed]
  • Paxinos G, Franklin KBJ. The Mouse Brain in Stereotaxic Coordinates. 2nd ed. Academic Press; Scan Diego: 2003.
  • Paxinos G, Halliday GM, Watson C, Koutcherov Y, Wang H. Atlas of the Developing Mouse Brain at E17.5, P0 and P6. Academic Press; New York: 2006.
  • Petiet AE, Kaufman MH, Goddeeris MM, Brandenburg J, Elmore SA, Johnson GA. High-resolution magnetic resonance histology of the embryonic and neonatal mouse: a 4D atlas and morphologic database. Proc Natl Acad Sci U S A. 2008;105:12331–12336. [PubMed]
  • Rosas HD, Feigin AS, Hersch SM. Using advances in neuroimaging to detect, understand, and monitor disease progression in Huntington's disease. NeuroRx. 2004;1:263–272. [PubMed]
  • Sidman RL, Rakic P. Neuronal migration, with special reference to developing human brain: a review. Brain Res. 1973;62:1–35. [PubMed]
  • Thompson PM, Hayashi KM, de Zubicaray G, Janke AL, Rose SE, Semple J, Herman D, Hong MS, Dittmer SS, Doddrell DM, Toga AW. Dynamics of gray matter loss in Alzheimer's disease. J Neurosci. 2003;23:994–1005. [PubMed]
  • Tuch DS, Reese TG, Wiegell MR, Makris N, Belliveau JW, Wedeen VJ. High angular resolution diffusion imaging reveals intravoxel white matter fiber heterogeneity. Magn Reson Med. 2002;48:577–582. [PubMed]
  • Tuch DS, Reese TG, Wiegell MR, Wedeen VJ. Diffusion MRI of complex neural architecture. Neuron. 2003;40:885–895. [PubMed]
  • Tyszka JM, Readhead C, Bearer EL, Pautler RG, Jacobs RE. Statistical diffusion tensor histology reveals regional dysmyelination effects in the shiverer mouse mutant. Neuroimage. 2006;29:1058–1065. [PMC free article] [PubMed]
  • Uylings H, Van Eden C, Parnavelas J, Kalsbeek A. The prenatal and Postnatal Development of Rat Cerebral Cortex. In: Kolb B, Tees R, editors. The Cerebral Cortex of the Rat. MIT Press; Boston, MA: 1990. pp. 35–76.
  • Verma R, Mori S, Shen D, Yarowsky P, Zhang J, Davatzikos C. Spatiotemporal maturation patterns of murine brain quantified by diffusion tensor MRI and deformation-based morphometry. Proc Natl Acad Sci U S A. 2005;102:6978–6983. [PubMed]
  • Wang Y, Zhang J, Mori S, Nathans J. Axonal growth and guidance defects in Frizzled3 knock-out mice: a comparison of diffusion tensor magnetic resonance imaging, neurofilament staining, and genetically directed cell labeling. J Neurosci. 2006;26:355–364. [PubMed]
  • Wedeen V, Reese T, Tuch DS, Wieigel MR, Dou JG, Weisskoff R, Chessler D. Mapping fiber orientation spectra in cerebral white matter with Fourier-transform diffusion MRI. Annual meeting of the international society of Magnetic Resonance in Medicine (ISMRM); Denver, Colorodo, USA. 2000. p. 82.
  • Weninger WJ, Mohun T. Phenotyping transgenic embryos: a rapid 3-D screening method based on episcopic fluorescence image capturing. Nat Genet. 2002;30:59–65. [PubMed]
  • Wingert F. Biometrical analysis of growth functions of cerebral component parts and body weight. J Hirnforsch. 1969;11:133–197. [PubMed]
  • Woods RP, Grafton ST, Holmes CJ, Cherry SR, Mazziotta JC. Automated image registration: I. General methods and intrasubject, intramodality validation. J Comput Assist Tomogr. 1998a;22:139–152. [PubMed]
  • Woods RP, Grafton ST, Watson JD, Sicotte NL, Mazziotta JC. Automated image registration: II. Intersubject validation of linear and nonlinear models. J Comput Assist Tomogr. 1998b;22:153–165. [PubMed]
  • Xu D, Mori S, Shen D, van Zijl PC, Davatzikos C. Spatial normalization of diffusion tensor fields. Magn Reson Med. 2003;50:175–182. [PubMed]
  • Xue R, van Zijl PCM, Crain BJ, Solaiyappan M, Mori S. In vivo three-dimensional reconstruction of rat brain axonal projections by diffusion tensor imaging. Magn. Reson. Med. 1999;42:1123–1127. [PubMed]
  • Zhang J, Chen YB, Hardwick JM, Miller MI, Plachez C, Richards LJ, Yarowsky P, van Zijl P, Mori S. Magnetic resonance diffusion tensor microimaging reveals a role for Bcl-x in brain development and homeostasis. J Neurosci. 2005;25:1881–1888. [PubMed]
  • Zhang J, Richards LJ, Yarowsky P, Huang H, van Zijl PC, Mori S. Three-dimensional anatomical characterization of the developing mouse brain by diffusion tensor microimaging. Neuroimage. 2003;20:1639–1648. [PubMed]