Elsevier

NeuroImage

Volume 49, Issue 1, 1 January 2010, Pages 205-216
NeuroImage

Neurite density from magnetic resonance diffusion measurements at ultrahigh field: Comparison with light microscopy and electron microscopy

https://doi.org/10.1016/j.neuroimage.2009.08.053Get rights and content

Abstract

Due to its unique sensitivity to tissue microstructure, diffusion-weighted magnetic resonance imaging (MRI) has found many applications in clinical and fundamental science. With few exceptions, a more precise correspondence between physiological or biophysical properties and the obtained diffusion parameters remain uncertain due to lack of specificity. In this work, we address this problem by comparing diffusion parameters of a recently introduced model for water diffusion in brain matter to light microscopy and quantitative electron microscopy. Specifically, we compare diffusion model predictions of neurite density in rats to optical myelin staining intensity and stereological estimation of neurite volume fraction using electron microscopy. We find that the diffusion model describes data better and that its parameters show stronger correlation with optical and electron microscopy, and thus reflect myelinated neurite density better than the more frequently used diffusion tensor imaging (DTI) and cumulant expansion methods. Furthermore, the estimated neurite orientations capture dendritic architecture more faithfully than DTI diffusion ellipsoids.

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)

  • SukstanskiiA.L. et al.

    In vivo lung morphometry with hyperpolarized 3He diffusion MRI: theoretical background

    J. Magn. Reson.

    (2008)
  • SukstanskiiA.L. et al.

    Effects of permeable boundaries on the diffusion-attenuated MR signal: insights from a one-dimensional model

    J. Magn. Reson.

    (2004)
  • TangY. et al.

    Myelinated nerve fibres in the subcortical white matter of cerebral hemispheres are preserved in alcoholic subjects

    Brain Res.

    (2004)
  • AkaikeH.

    New look at statistical-model identification

    IEEE Trans. Automat. Contr. AC19

    (1974)
  • AmuntsK. et al.

    Postnatal development of interhemispheric asymmetry in the cytoarchitecture of human area 4

    Anat. Embryol.(Berl)

    (1997)
  • AssafY. et al.

    AxCaliber: a method for measuring axon diameter distribution from diffusion MRI

    Magn. Reson. Med.

    (2008)
  • Bar-ShirA. et al.

    High b-value q-space diffusion MRS of nerves: structural information and comparison with histological evidence

    NMR Biomed.

    (2008)
  • BasserP.J.

    Inferring microstructural features and the physiological state of tissues from diffusion-weighted images

    NMR Biomed.

    (1995)
  • BasserP.J.

    New histological and physiological stains derived from diffusion-tensor MR images

    Imaging Brain Struct. Funct.

    (1997)
  • CallaghanP.T.

    Principles of nuclear magnetic resonance microscopy

    (1991)
  • ChinC.L. et al.

    Biexponential diffusion attenuation in the rat spinal cord: computer simulations based on anatomic images of axonal architecture

    Magn. Reson. Med.

    (2002)
  • ColonE.J.

    Quantitative cytoarchitectonics of the human cerebral cortex in schizophrenic dementia

    Acta Neuropathol. (Berl)

    (1972)
  • FalangolaM.F. et al.

    Age-related non-Gaussian diffusion patterns in the prefrontal brain

    J. Magn. Reson. Imaging

    (2008)
  • HardinR.H. et al.

    McLaren's improved snub cube and other new spherical designs in three dimensions

    Discrete Comput. Geom.

    (1996)
  • Cited by (246)

    View all citing articles on Scopus
    View full text