Ahmad Raza Khan1, Andery Chuhutin1, Ove Wiborg2, Christopher D Kroenke3, Jens R Nyengaard4, Brian Hansen1, and Sune Nørhøj Jespersen1,5
1Center of Functionally Integrative Neuroscience, Aarhus University, Aarhus, Denmark, 2Centre for Psychiatric Research, Aarhus University, Aarhus, Denmark, 3Advanced Imaging Research Center, Portland, OR, United States, 4Stereology and Electron Microscopy Laboratory, Aarhus University, Aarhus, Denmark, 5Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark
Synopsis
Biophysical
modelling of diffusion MRI data allows detection of specific tissue
microstructures such as neurite density. However, histological validation of MR-derived
indication of microstructural alteration is limited due extensive time labour
and invasive character, even though histological validation is crucial because
it remains the gold standard. The present study applies Matlab based image
processing and analysis tools to compute histological neurite density to
validate diffusion MRI based neurite density changes in the amygdala of chronic
mild stress rat brains. The image processing and analyses provides novel tools
to validate diffusion data robustly. Introduction:
Alteration
in shape, size and packing density of cellular and sub cellular microstructure
affect the diffusion of water, so the diffusion MRI (dMRI) signal ; however
traditional dMRI parameters are not sensitive enough to probe specific
microstructural components. Modeling of dMRI data has successfully detected
changes in cellular and sub-cellular microstructure, for example neurite
density (ν), axonal
water fraction, intra and extra axonal diffusivity, to name a few (1, 2, 3). Neurite
density parameters are sensitive and potentially useful in disease, diagnosis
and monitoring treatment outcome (1, 2). However, histological
validation of these findings is limited, but highly valuable, and may provide more
pertinent explanations of parameter alteration as well as to establish the usefulness
of model parameters (4). Although histology is considered a gold standard,
correlating MRI to histology is often unfeasible because it requires extensive
time and labor, and is thus a low
throughput technique. Here we employed
morphological operation tools included with the image processing toolbox of Matlab
to compute histology based neurite, nuclear and glial volume fractions from the
amygdala region of control and chronic mild stress (CMS) rat brains. The novel
image processing and analysis framework can process and analyze the images robustly
with high throughput.
Material and Methods
Twenty four
adult wistar rats (Taconic, Denmark) were used in the present study. The CMS paradigm
was implemented (5) and distinguished as anhedonic and resilient on the basis
of sucrose consumption test, described elsewhere (5). Brains were perfusion
fixed and post-fixed in buffered formalin. After ex-vivo dMRI (TR/TE =6500 /26
ms, Δ/δ 15/5 ms, 12 directions, 14 b-values (0-8000 s/mm2) and 250 µm isotropic resolution), the
same hemispheres were sectioned into 60 µm thick sections on a vibratome. One
set of tissue sections were double stained with lipophilic dye (DiI), and
nuclear stain (Hoechst). Another set of tissue sections were immunostained with
ALDH1L1, an astrocyte marker. Image stacks
in 3D from the amygdala region were acquired from DiI and Hoechst stained
sections on a confocal microscope (Figure 1a, 1d). Multiple images were
captured on a light microscope from the amygdala region of immunostained
(ALDH1L1) sections (Figure 1g). Images in 3D underwent depth dependent
intensity correction and were subsequently imported as ‘tiff’ images in Matlab,
followed by contrast enhancement (Figure 2b, e, h). Contrast enhanced images
were thresholded to binary images and morphologically opened with rectangular structuring elements. Binary
images were compared with contrast enhanced image to visually confirm
successful delineation of targeted microstructure (Figure 2c, f, i). Pixels in the binary images were
counted and normalized to approximate neurite, nuclear and astrocyte volume
fractions. MRI and histological data were separately fit in Matlab to a linear
mixed effect models with animals as random effects and group as a fixed effect.
Figures 2 and 3 show estimated means with confidence intervals (95%) providing
an estimate of fixed effect size and variability. F tests were used to assess
the level of significance.
Results
Biophysical
modeling of dMRI data based neurite density model showed significant increase
in neurite density only in amygdala in anhedonic (p = 0.008) and resilient groups (p = 0.011) (Figure 2).
Histology based neurite density was also significantly increased in anhedonic
and resilient groups (p = 0.018 and p = 0.048) respectively, in comparison to control,
while nuclear and astrocyte densities showed modest non-significant decrease in
both stressed group (Figure 3).
Discussion
The present study demonstrates
histological validation of dMRI based neurite density increase in the amygdala
of CMS rat brains. Histological validation of MR based findings establishes the
microstructural underpinnings of the model parameters and provides a clear
understanding of microstructural alterations driving observed parameter.
Despite the fundamental difference between the MRI and microscopy modalities, previous
quantitative analysis has shown a significant correlation between the two (2, 4).
Here, we present a novel image processing and analysis framework of
histological data to validate dMRI derived neurite density. Our analysis
comparing amygdala from stressed rat to normal controls provides results in
agreement with a previous neuronal tracing study which found dendritic
hypertrophy in amygdala after chronic stress (6). The modest decrease in
nuclear and astrocyte volume fraction reported here is also supported by
previous stereology based findings in major depressive disorder (7).
Conclusion
Biophysical
modelling combined with dMRI revealed increases in neurite density in the
amygdala of stressed rats. We developed novel
robust and high throughput image processing and analysis tools validating these
findings with histology and immunohistochemistry. These techniques can be
applied to validate specificity of MRI model parameters in an unbiased and
efficient manner.
Acknowledgements
Lundbeck Foundation grant
R83–A7548 and Simon Fougner Hartmanns Familiefond. AC and BH acknowledge supported from
NIH 1R01EB012874-01. The authors wish to thank Lippert’s Foundation and
Korning’s Foundation for financial support. The 9.4T lab was made possible by
funding from the Danish Research Counsil's Infrastructure program, the Velux
Foundations, and the Department of Clinical Medicine, AU. References
1. Jespersen,
S. N., et al. Modeling dendrite density from magnetic resonance
diffusion measurements. NeuroImage (2007); 34(4): 1473-1486.
2. Jespersen, S., et al. Neurite
density from magnetic resonance diffusion measurements at ultrahigh field:
comparison with light microscopy and electron microscopy. NeuroImage (2010); 49:
205-216.
3. Assaf,
Y., et al. New modeling and experimental framework to characterize hindered and
restricted water diffusion in brain white matter. Magnetic Resonance in
Medicine (2004); 52(5): 965-978.
4. Jespersen, S. N. et al. Determination of
axonal and dendritic orientation distributions within the developing cerebral
cortex by diffusion tensor imaging. Medical
Imaging, IEEE Transactions on (2012); 31(1):
16-32.
5. Wiborg, O., Chronic
mild stress for modeling anhedonia. (2013); Cell and tissue research 354(1):
155-169.
6. Vyas, A., et al. Chronic
stress induces contrasting patterns of dendritic remodeling in hippocampal and
amygdaloid neurons. The Journal of Neuroscience (2002); 22(15): 6810-6818.
7. Bowley, M. P., et al. Low glial
numbers in the amygdala in major depressive disorder. Biological psychiatry (2002);
52(5): 404-412.