Neurite density from magnetic resonance diffusion measurements at ultrahigh field: Comparison with light microscopy and electron microscopy
Introduction
The utility of diffusion weighted imaging, in particular in acute stroke (Hjort et al., 2005, Moseley et al., 1990), continues to drive research in the underlying biophysical mechanisms. During a typical diffusion experiment, water molecules in living tissue probe length scales on the order of 5 to 20 μm, making diffusion sensitive to a plethora of microstructural and physiological detail. The exquisite sensitivity of diffusion weighted MRI to tissue microstructure fuels the hope that the MRI scanner may eventually serve as an in vivo microscope of the brain and body. This would allow scientists and clinicians to obtain information not only about fiber tract trajectories (Ciccarelli et al., 2008), but also more detailed traits of tissue morphology such as, e.g., cell and fiber density, compartment sizes, and various cytoarchitectural attributes. Many studies have shown that parameters routinely obtained from diffusion weighted MRI, e.g. fractional anisotropy and mean diffusivity, act as biomarkers for different tissue properties (Basser, 1995, Basser, 1997, Gupta et al., 2008, Kozlowski et al., 2008, Tofts, 2003). Such MRI-derived physiological stains are very useful, but often suffer the fundamental disadvantage of lacking specificity—thus, more direct measures of histological characteristics would be advantageous. In order to achieve this goal, careful biophysical modeling of diffusion and independent validation using existing histological methods is a prerequisite. An early example of this approach demonstrated that membrane permeability and cellular surface to volume fraction could be obtained from modeling diffusion in packed erythrocytes (Latour et al., 1994). More recently, Yablonskiy et al. modeled diffusion in the airways of the lung, and showed that quantitative geometric information about the alveoli could be extracted from diffusion MRI (Sukstanskii and Yablonskiy, 2008, Woods et al., 2006, Yablonskiy et al., 2002). Two groups have recently obtained good agreement to histology in measuring axon diameters with q-space imaging, (Bar-Shir and Cohen, 2008, Ong et al., 2008). However, the general application of q-space methods for inferring microstructure requires careful interpretation, for example in the presence of permeable barriers (Sukstanskii et al., 2004). Assaf et al. (2008) accomplished good agreement between histological measurements of axonal diameter distributions and diffusion MRI by utilizing their modeling framework “AxCaliber.”
Recently, we introduced a biophysical model of diffusion in the brain (Jespersen et al., 2007). The model was based on a biophysical description of brain microstructure, and was shown to describe high b-factor diffusion weighted MRI experiments of a baboon brain better than alternative models. The fundamental assumption of the model is that diffusion can be described in terms of two non-exchanging components. One component is associated with diffusion in cylindrically symmetric structures, such as cell processes with exchange of water being sufficiently slow to be considered impermeable on the time scale of the diffusion experiment. Dendrites and axons, collectively termed neurites, were assumed to fulfill these assumptions. The net signal from this component then arises as a sum of the signal from all neurites weighted by an orientation distribution function, i.e. a probability density function specifying the number of neurites in every direction. The second component of the diffusion signal accounts for diffusion everywhere else, in particular in cell bodies, extracellular space, and glia cells. Here diffusion is hindered and described in terms of Gaussian isotropic diffusion with an effective diffusion constant. Several cytoarchitectural parameters can be extracted from this framework, one of the more interesting being the neurite density. Noninvasive, and eventually in vivo access to this parameter could be useful in several research areas, e.g. brain mapping (Schleicher et al., 1999), postnatal ontogeny (Amunts et al., 1997), aging (Stark and Pakkenberg, 2004), comparative neuroanatomy (Sherwood et al., 2004, Zilles et al., 1986), multiple sclerosis (Peterson et al., 2001, Trapp et al., 1998), schizophrenia (Colon, 1972), Alzheimer's disease (Stark et al., 2005), and alcoholism (Tang et al., 2004).
The purpose of the present work was to compare model parameters to measurements obtained from light microscopy and quantitative electron microscopy (stereology) in order to evaluate the interpretation of the model parameters. In particular, MR diffusion measurements of neurite density obtained by applying the diffusion model were compared to measures of neurite density and cell density at the light and electron microscopic level. The results were contrasted to an identical analysis using parameters from the DTI and cumulant expansions. Furthermore, we assessed the accuracy of the estimated neurite orientation functions qualitatively.
Section snippets
Theory
We analyzed the data in terms of the model presented previously (Jespersen et al., 2007), as well as an extension to account more directly for anisotropic diffusion in the extra-cellular space, and compared to two other diffusion models in common use, the diffusion tensor model and a fourth order cumulant expansion. In the model introduced in (Jespersen et al., 2007) (Model I), the signal S consisted of two non-exchanging components, one with cylindrical symmetry describing diffusion in
Results and discussion
In Fig. 2, images of Akaike's information criterion (AIC) are shown for models I, II, DTI, and CumExp. AIC provides a means of comparing models with differing numbers of free parameters by penalizing chi-square by two times the number p of free parameters in the model: AIC = χ2 + 2p (Akaike, 1974). From the AIC images, we conclude that DTI accounts only poorly for the variation of the diffusion signal with high diffusion weighting, as expected when bD≫1. The most faithful description of the
Conclusion
In this work, we have compared output from various diffusion models to a number of different measures obtained from light microscopy and stereological methods applied to electron microscopy. Models I and II are two variations of a biophysical diffusion model introduced in (Jespersen et al., 2007), whereas DTI and kurtosis models essentially are cumulant expansions of the spin phase distribution (i.e., the MR signal) terminated at second (DTI) or fourth order (CumExp). The cumulant expansion is
Acknowledgments
The authors wish to thank the Danish National Research foundation (CFIN and inSpin), and the Danish Council for Strategic Research, for funding. SNJ thanks Martin Carlsen and Daniel Otykier for technical assistance. The histological lab assistance of Ms. Dorete Jensen is greatly acknowledged. The authors thank two anonymous reviewers whose comments led to improvements of the paper.
References (48)
- et al.
Does diffusion kurtosis imaging lead to better neural tissue characterization? A rodent brain maturation study
NeuroImage
(2009) - et al.
Diffusion-based tractography in neurological disorders: concepts, applications, and future developments
Lancet Neurol.
(2008) - et al.
Effect of impermeable boundaries on diffusion-attenuated MR signal
J. Magn. Reson.
(2006) - et al.
The effect of impermeable boundaries of arbitrary geometry on the apparent diffusion coefficient
J. Magn. Reson.
(2008) - et al.
DTI derived indices correlate with immunohistochemistry obtained matrix metalloproteinase (MMP-9) expression in cellular fraction of brain tuberculoma
J. Neurol. Sci.
(2008) - et al.
Mapping the parameter space of a T2-dependent model of water diffusion MR in brain tissue
Magn. Reson. Imaging
(2006) - et al.
Modeling dendrite density from magnetic resonance diffusion measurements
NeuroImage
(2007) - et al.
In vivo visualization of displacement-distribution-derived parameters in q-space imaging
Magn. Reson. Imaging
(2008) - et al.
Indirect measurement of regional axon diameter in excised mouse spinal cord with q-space imaging: simulation and experimental studies
NeuroImage
(2008) - et al.
Observer-independent method for microstructural parcellation of cerebral cortex: a quantitative approach to cytoarchitectonics
NeuroImage
(1999)
In vivo lung morphometry with hyperpolarized 3He diffusion MRI: theoretical background
J. Magn. Reson.
Effects of permeable boundaries on the diffusion-attenuated MR signal: insights from a one-dimensional model
J. Magn. Reson.
Myelinated nerve fibres in the subcortical white matter of cerebral hemispheres are preserved in alcoholic subjects
Brain Res.
New look at statistical-model identification
IEEE Trans. Automat. Contr. AC19
Postnatal development of interhemispheric asymmetry in the cytoarchitecture of human area 4
Anat. Embryol.(Berl)
AxCaliber: a method for measuring axon diameter distribution from diffusion MRI
Magn. Reson. Med.
High b-value q-space diffusion MRS of nerves: structural information and comparison with histological evidence
NMR Biomed.
Inferring microstructural features and the physiological state of tissues from diffusion-weighted images
NMR Biomed.
New histological and physiological stains derived from diffusion-tensor MR images
Imaging Brain Struct. Funct.
Principles of nuclear magnetic resonance microscopy
Biexponential diffusion attenuation in the rat spinal cord: computer simulations based on anatomic images of axonal architecture
Magn. Reson. Med.
Quantitative cytoarchitectonics of the human cerebral cortex in schizophrenic dementia
Acta Neuropathol. (Berl)
Age-related non-Gaussian diffusion patterns in the prefrontal brain
J. Magn. Reson. Imaging
McLaren's improved snub cube and other new spherical designs in three dimensions
Discrete Comput. Geom.
Cited by (246)
Neonatal Nucleus Accumbens Microstructure Modulates Individual Susceptibility to Preconception Maternal Stress in Relation to Externalizing Behaviors
2024, Journal of the American Academy of Child and Adolescent PsychiatryAxial and radial axonal diffusivities and radii from single encoding strongly diffusion-weighted MRI
2023, Medical Image AnalysisComparison of different neurite density metrics with brain asymmetry evaluation
2023, Zeitschrift fur Medizinische Physik