- 1Department of Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, USA
- 2The Brain Observatory, San Diego, CA, USA
- 3Department of Radiology, University of California San Diego, San Diego, CA, USA
Diffusion weighted imaging (DWI) has provided unparalleled insight into the microscopic structure and organization of the central nervous system. Diffusion tensor imaging (DTI) and other models of the diffusion MRI signal extract microstructural properties of tissues with relevance to the normal and injured brain. Despite the prevalence of such techniques and applications, accurate and large-scale validation has proven difficult, particularly in the human brain. In this report, human brain sections obtained from a digital public brain bank were employed to quantify anisotropy and fiber orientation using structure tensor analysis. The derived maps depict the intricate complexity of white matter fibers at a resolution not attainable with current DWI experiments. Moreover, the effects of multiple fiber bundles (i.e., crossing fibers) and intravoxel fiber dispersion were demonstrated. Examination of the cortex and hippocampal regions validated-specific features of previous in vivo and ex vivo DTI studies of the human brain. Despite the limitation to two dimensions, the resulting images provide a unique depiction of white matter organization at resolutions currently unattainable with DWI. The method of analysis may be used to validate tissue properties derived from DTI and alternative models of the diffusion signal.
Introduction
The brain is a network with multiple levels of connectivity. The connections of each neuron determine its functional properties, and the functionality of the brain as a whole is determined by both these local and large-scale connections. Since anatomical and functional connectivity are intrinsically linked, a major goal of neuroscience has been to construct comprehensive maps of brain connectivity, or “connectomes” (Sporns et al., 2005; Sporns, 2011). Given the many levels of connectivity, from the microscopic to the macroscopic, the tools used to create such maps vary (Leergaard et al., 2012). No single modality is currently capable of reconstructing all connections in the human brain (Kleinfeld et al., 2011). Thus, a number of methods are needed to bridge the gap from the connectivity of individual neurons to those that measure large-scale patterns of anatomical connectivity.
Diffusion weighted imaging (DWI) has become an indispensible tool to investigate the large-scale anatomical connectivity in the brain. As a completely non-invasive modality, DWI has enabled unprecedented insights into brain microstructure and connectivity. Technological advances in the acquisition, analysis, and modeling of DWI data have propelled investigations of the normal brain architecture and the consequences of disease and injury. Despite these achievements, cross-validation of diffusion imaging and microscopy has been somewhat limited. Since DWI measures the mobility of water molecules, it is an indirect measure of structure. Inferences and mathematical models form the basis to derive structural information from the DWI signal. Diffusion tensor imaging (DTI) employs a single tensor model to obtain measures of primary fiber orientation and anisotropy. However, the DTI model has well-known limitations in brain regions that are composed of more than a single coherent fiber bundle, which occurs at high proportions in the human brain (Jeurissen et al., 2012). To remedy these limitations, other models of the DWI signal have been developed to resolve complex fiber orientations and more accurately map intersecting fiber pathways (Wedeen et al., 2012). Additionally, other interesting properties of brain tissue that may be relevant and important in both health and disease have been extracted from models of the DWI signal, such as neurite densities (Jespersen et al., 2010; Zhang et al., 2012), fiber dispersion (Sotiropoulos et al., 2012), and axonal diameters (Barazany et al., 2009; Zhang et al., 2011). However, robust methods of validating these features using other modalities have been limited.
Microscopy remains the standard for definitive examination of tissue microstructure. Post-mortem assessments of tissue organization using optical or electron microscopy are essential to directly visualize cellular and subcellular structure. Often, these methods are limited to small tissue regions and rely on selective sampling of brain tissues. In recent years, advances have been made in large-scale microscopic imaging of whole-brain sections (Denk and Horstmann, 2004; Dodt et al., 2007), including human brains (Annese, 2012). The vast increases in both resolution and coverage have necessitated quantitative methods to extract pertinent information and reduce the data to a manageable format for dissemination. With these advances, microscopy techniques are converging with DWI to enable a comprehensive map of brain microstructural anatomy and connectivity. This convergence will enable cross-validation and will provide insight into the enormous structural complexity of the brain. These insights will advance the understanding of nervous system disease and injury and therefore has considerable importance.
In this report, a method to quantify microscopic tissue properties was applied to human brain sections. Specifically, high-resolution images of whole-human brain sections obtained from a publically available collection underwent structure tensor analysis (Budde and Frank, 2012) to extract microscopic features of anisotropy and orientation, since these are the most common parameters derived from DTI. The work highlights specific features of the brain that are commonly investigated in DTI and other DWI studies, including the effect of crossing fibers and fiber dispersion. Insights to specific anatomical regions, including the cortex and hippocampal structures (Augustinack et al., 2010; Yassa et al., 2010) is also presented. The objective was to introduce the technique as a complement to previous and future DWI studies of the human brain, since accurate validation has been notoriously lacking (Hagmann et al., 2010). The results and methodology outlined here will be useful in validating DTI and advanced models of the diffusion MRI signal and reveal features of the microstructure unattainable with resolution of current DWI examinations.
Materials and Methods
Data Analysis
Human brain sections stained were obtained from a public histological resource for the human brain that is available on The Brain Observatory's web site (2012). Frozen sections were cut at an interval of 70 μm, mounted on glass slides, and stained using a modified silver stain for myelin (Gallyas, 1979). Staining was performed in large-batches, so all materials analyzed in the context of this communication underwent the same protocol. The penetration of the silver salts during impregnation and homogeneous development of the salts into visible colloidal silver was tested empirically on the finished slides using a high-magnification, and high numerical aperture oil immersion objective (Olympus PlanApo N 60×, NA: 1.42). The digitization of sections was performed using a custom-engineered, large-format, slide-scanning microscope using a 20× objective. Images were acquired at a resolution of 0.37 μm2 per pixel. The analysis described below used images from 15 whole-brain coronal sections subsampled at a working resolution of 2.96 μm2 per pixel. After fixation in 4% formaldehyde phosphate-buffered solution at 4°C, the brain specimen was cryoprotected in increasing concentration of sucrose and embedded using stereotaxic landmarks (Talairach and Tournoux, 1988). Accordingly, sections were obtained approximately perpendicular to a line intersecting the anterior and posterior commissures (Annese, 2012). Further details of tissue processing, including sectioning, staining, and digitization are provided in Bartsch et al. (2012).
Data Analysis
Structure tensor analysis of histological sections was previously described for the rat brain (Budde and Frank, 2012) and was derived from original methods proposed by Bigun and Granlund (1987). All procedures were implemented as custom routines in Matlab. Briefly, images were convolved with a discrete approximation of 2D Gaussian derivative filters (Figure 1) separately along and to obtain the image partial derivatives, fx and fy, respectively, where σ is the width of the Gaussian window in pixels that defines the inner, or local scale. The structure tensor, J, of an image is a matrix representation of the image partial derivatives for each pixel
where w is a Gaussian convolution kernel of width σ that defines the outer, or integration scale. The primary orientation, corresponding to the direction of the largest eigenvector of the tensor, was obtained by
and anisotropy index was obtained by
where λ1 and λ2 are the largest and smallest eigenvalues of J, respectively. For visualization of images in their native resolution, the values of orientation, anisotropy, and staining intensity were reflected as hue, saturation, and brightness (HSB) images, respectively. Alternatively, downsampling was employed to mimic the resolutions achievable by MRI by summing f2x, fy2, and fxy from regularly spaced square subregions of a predetermined size (i.e., “voxels”) to derive J as shown above by omitting the outer scale. In this representation, the primary orientation and anisotropy were displayed as hue and brightness, respectively that mimic directionally encoded color (DEC) images derived from DTI. Additionally, fiber orientation distributions (FODs) were obtained for each voxel for visualization and quantification of the complete in-plane angular profile. The FODs were computed as histograms of pixel orientations, θ, within each voxel using 64 equally spaced bins. To quantify fiber dispersion, FODs were fit to a single Gaussian function with a mean orientation of μ and standard deviation of σ using a non-linear least squares routine in the Matlab Curve Fitting toolbox. The standard deviation was displayed as a scalar map to depict regional variations in intravoxel dispersion within the corpus callosum (CC).
Figure 1. Structure Tensor Analysis. Images from a region of high anisotropy in the corpus callosum (A) and a region with low anisotropy from the cortical gray matter (B) are shown as examples. Original images were convolved with Gaussian derivative filters, and maps of anisotropy and orientation were derived from the structure tensor matrix at each pixel in the image. Hue-saturation-brightness (HSB) images visually display the orientation, anisotropy, and intensity (inverted for these sections), respectively, in a single image. The angular histogram of orientation over the image are quantified and displayed as fiber orientation distributions (FODs).
Results
White Matter Microstructure
Microscopic tissue orientation and anisotropy was derived using structure tensor analysis applied to images of human brain sections acquired at sub-micron resolution. Myelin staining intensity was largely uniform across the cerebral white matter (Figure 2A) and did not distinguish individual tracts. In contrast, structure tensor analysis revealed individual tracts and the microscopic heterogeneity of the white matter. The derived orientation and anisotropy reflects the microscopic composition of white matter tracts (Figure 2B). Specifically, the CC is the most prominent white matter structure in the brain. As expected, the anisotropy of the CC was greater than the other white matter tracts in the brain. Along with regional variations in orientation indicated by hue, local variations in anisotropy were evident throughout the white matter. Since anisotropy in the current study reflects fibers oriented within the section, low anisotropy could be caused by through-plane fibers. Alternatively, regions that contain multiple fiber populations within a single voxel would also have an anisotropy less than that of coherent fiber bundles. These two effects could not be disambiguated by the current approach. Nonetheless, many large-tracts such as the optic radiations and pathways in the pons were readily visualized in the color-coded images.
Figure 2. Whole-section maps of fiber microstructure. Images of the (inverted) staining intensity (A) and directionally encoded color maps (DEC; B) were reconstructed at a resolution of 68 μm2. Staining intensity was nearly uniform in the white matter, but the DEC images highlight the intricate microstructure of white matter tracts. Numbers indicate the regions shown in greater detail in subsequent figures.
The effect of crossing fibers was further demonstrated in a region within the centrum semi-ovale, where fibers of the longitudinal fasciculus and corticospinal tract intersect with transcallosal fibers (Wedeen et al., 2012). DEC images were reconstructed with square pixels of 400 μm2 (Figure 3A), substantially smaller than those normally used for whole-brain studies at field strengths of 3T and below. Anisotropy was lower in the region of multiple fiber intersections than adjacent fiber bundles with a single, coherent orientation such as the CC. Reconstructions at high-resolution, 16 μm2, resulted in greater anisotropy than the lower resolution, which can be attributed to the decreased effect of partial volumes. The interdigitated fiber pathways were also more readily visualized in the high-resolution reconstructions (Figure 3B). Compared to the unmodified images at their original resolution, 0.37 μm2, the same images color-coded with the local fiber orientation clearly revealed that the intersection of major fiber pathways is composed of small fiber bundles. Likewise, the two separable peaks in the FOD depicted the existence of two orthogonal (in-plane) fiber populations.
Figure 3. Partial volume effects and crossing fibers. Anisotropy in regions containing multiple intersecting fiber tracts was lower compared to single, coherent fiber pathways (A). Reconstructions at higher resolution (B) improve the visualization of small, interdigitated fiber bundles in the region containing crossing fibers (B). Compared to the unmodified stained sections (C), crossing fibers were more apparent in the same images color-coded with orientation (D). The two peaks in the FOD highlight the presence of two discrete fiber populations within this subregion of the centrum semi-ovale (E).
In addition to crossing fibers, bending and dispersion of fibers within a single voxel can affect the estimated fiber orientation in diffusion MRI data. The CC is the most prominent structure in the human brain and is considered to consist of a single, coherent fiber pathway. Although the primary direction of fibers in the CC is largely consistent and follows the macroscopic curvature (Figure 4), there are regional variations in orientation that are visualized using structure tensor analysis of high-resolution microscopy data. Smaller fiber bundles that exhibit slightly differences in orientation are likely related to the foliated structure of the CC (Wiggins et al., 2010). The microscopic dispersion was quantified by fitting the FODs from each reconstructed voxel to a Gaussian function (Figure 4), where the standard deviation of the Gaussian fit curve reflects the intravoxel dispersion of fibers. The increased dispersion within the medial and lateral regions of the CC was evident as a broadening of the FODs and the corresponding map of dispersion. Across the whole CC in this section, the mean intravoxel dispersion was 14.6 ± 3.6°.
Figure 4. Angular dispersion of the corpus callosum. In the body of the CC (A), the microscopic heterogeneity is apparent on images color coded by orientation reconstructed at high resolution (B and C). The FODs (D) were fit to a Gaussian function to quantify the microscopic orientation dispersion (E). The medial and lateral aspects of the CC had greater dispersion in this section.
Cortical Microstructure
The human cerebral cortex has a highly intricate and complex architecture that has been investigated with DWI. High-resolution in vivo or ex vivo DTI has consistently identified the cortex as anisotropic with its predominant orientation perpendicular to the cortical surface (Heidemann et al., 2012). In an example from the frontal cortex (Figure 5), structure tensor analysis confirms the radial organization in the cortex having high anisotropy. While the white matter fibers within the gyri were perpendicular to the gyral surface throughout, white matter fibers had sharp bends into the cortex (Figure 5B) along the sulcal surface as shown previously with DWI (Takahashi et al., 2011). In the depths of the sulcus, the cortical myeloarchitecture consisted of fibers arranged predominantly parallel to the white- and gray matter interface (Figure 5D). The white matter fibers just superficial to the interface were largely composed of U-shaped association fibers that follow the trajectory of the sulcal depths (Catani et al., 2012).
Figure 5. Cortical Structure. The myelinated fibers of the cortex (A) are predominantly oriented perpendicular to the cortical surface. Within the gyrus, white matter fibers are oriented along the main axis of the gyrus. However, (B) fibers entering the cortex from the white matter turn at nearly right angles as evidenced by the abrupt changes in orientation (D) along the white-gray matter interface (dotted line). In cortex at the deepest portion of the sulci, the cortex has a mixture of radial and tangential fibers (C). The white matter fibers just below the deepest part of the sulci are oriented tangentially, consistent with the U-shaped association fibers.
The visual, or striate cortex (Figure 6), displays an architecture different than most other cortical regions. Again, the gray matter was predominantly anisotropic with a radial organization, but different cortical layers were exceptions. Fibers in layer 1 were arranged parallel to the cortical surface. Layers 4 and 6 consisted of fibers arranged both parallel and perpendicular to the cortical surface. This was evident as decreased anisotropy as well as the two opposing peaks in the FOD (Figure 6D), and may be the cause of decreased anisotropy measured with DTI at the gray and white matter interface (Guilfoyle et al., 2003). Although layers 2–3 and 5 have different myelin contents, the anisotropy and orientation information derived from these two layers was similar. The lamina-specific orientations are consistent with ex vivo DTI findings by Leuze et al. (2012).
Figure 6. Laminar-specific anisotropy in striate cortex. A section from the primary visual cortex (A) reveals the highly anisotropic myeloarchitecture of area V1 with a clear laminar subdivision. Layer 1 is predominantly tangential (parallel) to the cortical surface as evidenced by the high anisotropy (B) and orientation (red; C and D). In contrast, layers 2, 3, and 5 have a radial orientation (green; C and D) with high anisotropy (B), despite their different myelin content. Layers 4 and 6 are composed of both radial and tangential fibers as evidenced by the low anisotropy (B) in these regions and the orthogonal peaks in the in the overlaid FODs (D). The microscopic organization of the myelin-stained cortex is shown in (E).
Perforant Pathway
The identification of small structures using in vivo or ex vivo diffusion MRI could have implications as biomarkers of disease activity. The perforant pathway (PP) has been identified in the human brain in vivo (Yassa et al., 2010) and ex vivo (Augustinack et al., 2010) and has implications to aging and dementia. The PP is a fiber bundle within the medial temporal lobe that connects the hippocampus to the entorhinal cortex. Structure tensor analysis was applied to a human brain section obtained at the level of the hippocampus. Although the PP was evident as a diffuse, myelinated structure within the unmyelinated subiculum (Figure 7), its projection through the myelinated angular bundle was obscured on the myelin-stained section (Figure 7A). However, the map of anisotropy depicted the PP as a diffuse region of high anisotropy (Figure 7B). The orientation of the PP differentiated it from the angular bundle and was consistent with its projection from the entorhinal cortex to the hippocampus (Figure 7C).
Figure 7. Histological evidence of the perforant pathway. The perforant pathway is a set of fibers projecting from the entorhinal cortex through the angular bundle and subiculum to the hippocampal formation (A). In the myelin-stained section, the perforant pathway is seen as it penetrates the subiculum, but is not discriminated from the angular bundle. In the maps of anisotropy (B) and orientation (C), the perforant pathway is evident as a region of high anisotropy (arrows) projecting from the entorhinal cortex through the subiculum. AB, angular bundle; EC, entorhinal cortex; Hipp, hippocampal formation; PP, perforant pathway; Sub, subiculum.
Discussion
Diffusion MRI is unique in its ability to examine the microstructure and large-scale connectivity of the brain non-invasively. However, since measurements of water mobility are an indirect measure of structure, validation with other techniques that directly assess structure are important for robust interpretations of the DWI signal and parameters derived from it. In this report, we derived anisotropy and fiber orientations directly from stained histological sections of the human brain. The results demonstrate that structure tensor analysis is a robust method to extract microstructural features of the brain. The results directly demonstrate the limitations of DTI in regions of complex fiber crossings, and the effect of resolution of the measured anisotropy. The limitations of DTI are analyzed thoroughly in a recent paper by Jones et al. (2012). Furthermore, the method reveals the microstructural properties of anatomical regions that have been previously studied with DWI. For example, in addition to showing the anisotropy and radial organization of the cerebral cortex, the crossing fibers in specific cortical layers have implications to DTI and multi-fiber models. Visualization of other structures of interest, such as the PP, was facilitated by quantitative measures of anisotropy and direction, as these myelinated fibers are difficult to appreciate in myelin-stained sections. The method and results described herein may prove useful in relating microscopy-based measures of structure and connectivity with in vivo and ex vivo DWI (McNab et al., 2009; Miller et al., 2011, 2012) of the human brain.
The current work has several limitations. Among them, the lack of a third dimension is the most important consideration that limits the potential of the work and has been previously discussed (Budde and Frank, 2012). While the crossing fiber problem could be directly visualized at high-resolution (Figure 3), the lack of three-dimensional information precluded reconstruction of fiber tracts and connectivity maps that reveal unique properties of brain organization (Wedeen et al., 2012). Microscopy techniques to image the brain in three-dimensions have been demonstrated for small tissue samples, such as rodent brains (Dodt et al., 2007; Keller and Dodt, 2011; Mikula et al., 2012). Large-scale application to the human brain is challenging, but not impossible (Bartsch et al., 2012). DTI and other models of characterizing brain structure with diffusion MRI have shown tremendous potential to study the normal and diseased microscopic anatomy. However, validation remains an important, but often overlooked, aspect of studies. As we have shown, the analysis of stained histological sections provides valuable corroborating information that complements diffusion MRI findings from previous in vivo and ex vivo experiments. A second limitation of the current study is the use of myelin-stained sections. In the white matter, the myeloarchitecture accurately reflects the microstructure as a whole. This feature is exploited in polarized light microscopy to quantify the microstructural organization of tissue in thick sections, including through-plane fibers (Axer et al., 2011a,b). However, in the cortex and subcortical structures, the relatively greater proportions of other cells and subcellular elements, including neuronal somas, dendrites, and glial cells, all contribute to the measured diffusion signal. Since biological membranes are the primary impediments to water mobility in the brain, the use of lipophilic dyes (Budde and Frank, 2012) or serial electron microscopy (Denk and Horstmann, 2004) may also prove essential to validate the tissue properties being probed by advanced diffusion models.
The current work both supports observations highlighted in previous diffusion MRI studies and adds greater insight into the microanatomy of the human brain. In the visual cortex, Leuze et al. recently demonstrated the laminar-specific orientation of the human visual cortex and high-resolution using ex vivo DTI and PLM (Leuze et al., 2012). The current work supports the finding that the cortex is predominantly arranged radial to the cortical surface with notable exceptions. Brodman layer 1 is oriented tangentially to the cortical surface and has high anisotropy. In contrast, layer 4 exhibits both radial and tangential fibers and therefore has a much lower anisotropy. At the interface of the white and gray matter, low anisotropy is observed in this study and previous diffusion MRI studies (Miller et al., 2011). The radial and tangential organization of fibers in layer 6 is demonstrated. However, within the visual sulcus, the cortical fibers are orientated radial to the cortex, whereas the underlying white matter fibers are largely tangential to the surface. In the majority of brain regions, fibers that project to the cortex within the sulci have sharp bends at the white-gray matter interface (Takahashi et al., 2011), while fibers projecting to the gyri maintain their straight trajectories into the gyri. These features are highlighted in Figure 5 and may underlie the diffusion MRI findings that gyri and sulci have differing cortical fiber orientations (McNab et al., 2012) and axonal termination densities (Nie et al., 2012).
The quantitative findings of intravoxel dispersion measured in the current study (14.6°) agree with previous studies of the rat brain (14.4; ) (Leergaard et al., 2010) or the owl monkey (12°) (Choe et al., 2012). These histological-derived measures of dispersion were generally consistent with estimates of the same parameter obtained with diffusion MRI using a model that incorporated intravoxel bending or fanning (Sotiropoulos et al., 2012) although variations in the regional pattern within the CC may be different. Since the CC is routinely considered to be composed of a single, coherent fiber population, it is often employed to derive a deconvolution kernel for multiple-fiber models (Tournier et al., 2004). Our high-resolution orientation maps of the CC (Figure 4C) highlight the microscopic variation within the CC that could potentially influence the accuracy of multiple-fiber models (Parker et al., 2013). The inclusion of multiple subjects and additional sections, including those obtained from different orientations, would provide a more comprehensive analysis of the effects of intravoxel dispersion. Importantly, since fractional anisotropy measured with DTI has a direct relationship with fiber density (Choe et al., 2012) but an inverse relationship with fiber dispersion (Zhang et al., 2012), these features might also be important physical phenomenon underlying FA variability in different regions of the brain and across different groups.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors would like to thank the staff of The Brain Observatory for their hard work creating the digital histological resource that was used in this study.
Funding
This project was funded through the Research and Education Initiative Fund, a component of the Advancing a Healthier Wisconsin endowment at the Medical College of Wisconsin (5520207 to Matthew D. Budde) and research grants from the National Eye Institute (R01 EY01835902 and ARRA R01 EY01835902S1 to Jacopo Annese) and the National Institute for Mental Health (R01MH084756 to Jacopo Annese).
References
The Brain Observatory. (2012). Concise Digital Atlas of the Human Brain [Online]. Available online at: http://thebrainobservatory.ucsd.edu/cortex-1 [Accessed 2012].
Annese, J. (2012). The importance of combining MRI and large-scale digital histology in neuroimaging studies of brain connectivity and disease. Front. Neuroinform. 6:13. doi: 10.3389/fninf.2012.00013
Augustinack, J. C., Helmer, K., Huber, K. E., Kakunoori, S., Zollei, L., and Fischl, B. (2010). Direct visualization of the perforant pathway in the human brain with ex vivo diffusion tensor imaging. Front. Hum. Neurosci. 4:42. doi: 10.3389/fnhum.2010.00042
Axer, H., Beck, S., Axer, M., Schuchardt, F., Heepe, J., Flucken, A., et al. (2011a). Microstructural analysis of human white matter architecture using polarized light imaging: views from neuroanatomy. Front. Neuroinform. 5:28. doi: 10.3389/fninf.2011.00028
Axer, M., Amunts, K., Grassel, D., Palm, C., Dammers, J., Axer, H., et al. (2011b). A novel approach to the human connectome: ultra-high resolution mapping of fiber tracts in the brain. Neuroimage 54, 1091–1101.
Barazany, D., Basser, P. J., and Assaf, Y. (2009). In vivo measurement of axon diameter distribution in the corpus callosum of rat brain. Brain 132, 1210–1220.
Bartsch, H., Maechler, P., and Annese, J. (2012). Automated determination of axonal orientation in the deep white matter of the human brain. Brain Connect. 2, 284–290.
Bigun, J., and Granlund, G. (1987). “Optimal orientation detection of linear symmetry,” in Proceedings of the IEEE First International Conference on Computer Vision (London), 433–438.
Budde, M. D., and Frank, J. A. (2012). Examining brain microstructure using structure tensor analysis of histological sections. Neuroimage 63, 1–10.
Catani, M., Dell'Acqua, F., Vergani, F., Malik, F., Hodge, H., Roy, P., et al. (2012). Short frontal lobe connections of the human brain. Cortex 48, 273–291.
Choe, A. S., Stepniewska, I., Colvin, D. C., Ding, Z., and Anderson, A. W. (2012). Validation of diffusion tensor MRI in the central nervous system using light microscopy: quantitative comparison of fiber properties. NMR Biomed. 25, 900–908.
Denk, W., and Horstmann, H. (2004). Serial block-face scanning electron microscopy to reconstruct three-dimensional tissue nanostructure. PLoS Biol. 2:e329. doi: 10.1371/journal.pbio.0020329
Dodt, H. U., Leischner, U., Schierloh, A., Jahrling, N., Mauch, C. P., Deininger, K., et al. (2007). Ultramicroscopy: three-dimensional visualization of neuronal networks in the whole mouse brain. Nat. Methods 4, 331–336.
Gallyas, F. (1979). Silver staining of myelin by means of physical development. Neurol. Res. 1, 203–209.
Guilfoyle, D. N., Helpern, J. A., and Lim, K. O. (2003). Diffusion tensor imaging in fixed brain tissue at 7.0 T. NMR Biomed. 16, 77–81.
Hagmann, P., Cammoun, L., Gigandet, X., Gerhard, S., Grant, P. E., Wedeen, V., et al. (2010). MR connectomics: principles and challenges. J. Neurosci. Methods 194, 34–45.
Heidemann, R. M., Anwander, A., Feiweier, T., Knosche, T. R., and Turner, R. (2012). k-space and q-space: combining ultra-high spatial and angular resolution in diffusion imaging using ZOOPPA at 7 T. Neuroimage 60, 967–978.
Jespersen, S. N., Bjarkam, C. R., Nyengaard, J. R., Chakravarty, M. M., Hansen, B., Vosegaard, T., et al. (2010). Neurite density from magnetic resonance diffusion measurements at ultrahigh field: comparison with light microscopy and electron microscopy. Neuroimage 49, 205–216.
Jeurissen, B., Leemans, A., Tournier, J. D., Jones, D. K., and Sijbers, J. (2012). Investigating the prevalence of complex fiber configurations in white matter tissue with diffusion magnetic resonance imaging. Hum. Brain Mapp. doi: 10.1002/hbm.22099. [Epub ahead of print].
Jones, D. K., Knosche, T. R., and Turner, R. (2012). White matter integrity, fiber count, and other fallacies: the do's and don'ts of diffusion MRI. Neuroimage. doi: 10.1016/j.neuroimage.2012.06.081. [Epub ahead of print].
Keller, P. J., and Dodt, H. U. (2011). Light sheet microscopy of living or cleared specimens. Curr. Opin. Neurobiol. 22, 138–143.
Kleinfeld, D., Bharioke, A., Blinder, P., Bock, D. D., Briggman, K. L., Chklovskii, D. B., et al. (2011). Large-scale automated histology in the pursuit of connectomes. J. Neurosci. 31, 16125–16138.
Leergaard, T. B., Hilgetag, C. C., and Sporns, O. (2012). Mapping the connectome: multi-level analysis of brain connectivity. Front. Neuroinform. 6:14. doi: 10.3389/fninf.2012.00014
Leergaard, T. B., White, N. S., de Crespigny, A., Bolstad, I., D'Arceuil, H., Bjaalie, J. G., et al. (2010). Quantitative histological validation of diffusion MRI fiber orientation distributions in the rat brain. PloS ONE 5:e8595. doi: 10.1371/journal.pone.0008595
Leuze, C. W., Anwander, A., Bazin, P. L., Dhital, B., Stuber, C., Reimann, K., et al. (2012). Layer-specific intracortical connectivity revealed with diffusion MRI. Cereb. Cortex. doi: 10.1093/cercor/bhs311. [Epub ahead of print].
McNab, J. A., Jbabdi, S., Deoni, S. C., Douaud, G., Behrens, T. E., and Miller, K. L. (2009). High resolution diffusion-weighted imaging in fixed human brain using diffusion-weighted steady state free precession. Neuroimage 46, 775–785.
McNab, J. A., Polimeni, J. R., Wang, R., Augustinack, J. C., Fujimoto, K., Player, A., et al. (2012). Surface based analysis of diffusion orientation for identifying architectonic domains in the in vivo human cortex. Neuroimage 69C, 87–100.
Mikula, S., Binding, J., and Denk, W. (2012). Staining and embedding the whole mouse brain for electron microscopy. Nat. Methods 9, 1198–1201.
Miller, K. L., McNab, J. A., Jbabdi, S., and Douaud, G. (2012). Diffusion tractography of post-mortem human brains: optimization and comparison of spin echo and steady-state free precession techniques. Neuroimage 59, 2284–2297.
Miller, K. L., Stagg, C. J., Douaud, G., Jbabdi, S., Smith, S. M., Behrens, T. E., et al. (2011). Diffusion imaging of whole, post-mortem human brains on a clinical MRI scanner. Neuroimage 57, 167–181.
Nie, J., Guo, L., Li, K., Wang, Y., Chen, G., Li, L., et al. (2012). Axonal fiber terminations concentrate on gyri. Cereb. Cortex 22, 2831–2839.
Parker, G. D., Marshall, D., Rosin, P. L., Drage, N., Richmond, S., and Jones, D. K. (2013). A pitfall in the reconstruction of fibre ODFs using spherical deconvolution of diffusion MRI data. Neuroimage 65, 433–448.
Sotiropoulos, S. N., Behrens, T. E., and Jbabdi, S. (2012). Ball and rackets: inferring fiber fanning from diffusion-weighted MRI. Neuroimage 60, 1412–1425.
Sporns, O., Tononi, G., and Kotter, R. (2005). The human connectome: a structural description of the human brain. PLoS Comput. Biol. 1:e42. doi: 10.1371/journal.pcbi.0010042
Takahashi, E., Dai, G., Rosen, G. D., Wang, R., Ohki, K., Folkerth, R. D., et al. (2011). Developing neocortex organization and connectivity in cats revealed by direct correlation of diffusion tractography and histology. Cereb. Cortex 21, 200–211.
Talairach, J., and Tournoux, P. (1988). Co-planar Stereotaxic Atlas of the Human Brain: 3-Dimensional Proportional System - an Approach to Cerebral Imaging. New York, NY: Thieme Medical Publishers.
Tournier, J. D., Calamante, F., Gadian, D. G., and Connelly, A. (2004). Direct estimation of the fiber orientation density function from diffusion-weighted MRI data using spherical deconvolution. Neuroimage 23, 1176–1185.
Wedeen, V. J., Rosene, D. L., Wang, R., Dai, G., Mortazavi, F., Hagmann, P., et al. (2012). The geometric structure of the brain fiber pathways. Science 335, 1628–1634.
Wiggins, C. J., Schäfer, A., Dhital, B., Le Bihan, D., and Turner, R. (2010). Imaging at 7T Reveals New Septated Fine Structure in the Human Corpus Callosum. Stockholm: Poster presented at ISMRM Annual Meeting.
Yassa, M. A., Muftuler, L. T., and Stark, C. E. (2010). Ultrahigh-resolution microstructural diffusion tensor imaging reveals perforant path degradation in aged humans in vivo. Proc. Natl. Acad. Sci. U.S.A. 107, 12687–12691.
Zhang, H., Hubbard, P. L., Parker, G. J., and Alexander, D. C. (2011). Axon diameter mapping in the presence of orientation dispersion with diffusion MRI. Neuroimage 56, 1301–1315.
Keywords: anisotropy, diffusion tensor imaging, structure tensor, crossing fibers, histology, fiber orientation
Citation: Budde MD and Annese J (2013) Quantification of anisotropy and fiber orientation in human brain histological sections. Front. Integr. Neurosci. 7:3. doi: 10.3389/fnint.2013.00003
Received: 25 October 2012; Paper pending published: 03 December 2012;
Accepted: 11 January 2013; Published online: 01 February 2013.
Edited by:
Alpay Özcan, Virginia Polytechnic and State University, USAReviewed by:
Robert Turner, Max Planck Institute for Human Cognitive and Brain Sciences, GermanyHao Huang, University of Texas Southwestern Medical Center, USA
Peng Sun, Washington University in St. Louis, USA
Copyright © 2013 Budde and Annese. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Matthew D. Budde, VA Medical Center-Research 151, 5000 West National Ave., Milwaukee, WI 53295, USA. e-mail: mdbudde@mcw.edu