4108

Susceptibility tensor imaging of ex vivo human hemibrain using 7T human MR scanner
Yuto Uchida1, Hyeong-Geol Shin2, Javier Redding-Ochoa3, Kengo Onda1, Alexander Barrett3, Adnan Bibic1,2, Juan C. Troncoso3, Peter van Zijl1,2, Kenichi Oishi1,4, and Xu Li1,2
1Department of Radiology, Johns Hopkins University, Baltimore, MD, United States, 2F.M. Kirby Research Center for Functional Brain Imaging, Kennedy Krieger Institute, Baltimore, MD, United States, 3Department of Pathology, Johns Hopkins University, Baltimore, MD, United States, 4The Richman Family Precision Medicine Center of Excellence in Alzheimer's Disease, Baltimore, MD, United States

Synopsis

Keywords: Preclinical Image Analysis, Microstructure

Motivation: High-resolution ex vivo susceptibility tensor imaging (STI) shows promise in visualizing detailed microstructural neuroanatomy related to tissue magnetic susceptibility contrast but is usually challenging with large human brain samples.

Goal(s): This work aimed to demonstrate a feasible protocol for STI on postmortem human hemibrain using a 7T human MR scanner.

Approach: De-identified human brain samples of the left hemisphere were collected, prepared, and scanned with a 3D multi-echo gradient echo sequence with 0.5 mm isotropic resolution.

Results: Multi-orientation quantitative susceptibility mapping and STI with a maximum of 12 orientations were obtained. Myelinated fibers and iron deposition in cortical substructures were visualized.

Impact: The proposed ex vivo MRI protocol is expected to be helpful to researchers interested in STI. The ex vivo STI acquired through this protocol may provide anatomical references for in vivo STI studies.

INTRODUCTION

Susceptibility tensor imaging (STI) aims at quantifying anisotropic tissue magnetic susceptibility from MR phase measures acquired at multiple sample orientations relative to the main magnetic field (B0) using a second-order tensor model1-3. High-resolution STI metrics, such as mean magnetic susceptibility (MMS), magnetic susceptibility anisotropy (MSA), and MSA-weighted principal eigenvector (PEV) maps, show promise in visualizing detailed microstructural neuroanatomy related to tissue magnetic susceptibility contrast3-5. Here, we aimed to test a STI scan protocol with 0.5 mm isotropic resolution using postmortem human hemibrains on a human 7T MRI scanner (Philips Healthcare, Best, the Netherlands).

MATERIALS AND METHODS

Brain sample
De-identified postmortem hemibrains of the left hemisphere were provided by the Johns Hopkins Brain Resource Center. At autopsy, the fresh brains weighed within the normal ranges and gross examinations revealed no abnormal findings. The hemibrains were fixed in 10% formaldehyde for more than 3 months before the preparation for ex vivo MRI scans to ensure complete fixation of the entire sample as performed in previous studies5, 6.
Designing and creating the rotatable ellipsoid-shaped container
To facilitate MRI scans for the hemibrain at multi-orientations along the x (right-left) and y (anterior-posterior) axes, a rotatable ellipsoid-shaped container was designed to hold the sample tightly. The container was 3D printed for an ellipsoidal shape with two equal radii of 87 mm and one larger radius of 112 mm (Figure 1). Four bases were also 3D printed to hold the container at different angles around the x-axis.
Pre-imaging steps
A set of procedures for the sample preparation from formalin fixation to MRI scan is summarized (Figure 2).
MRI scans
3D multi-echo GRE sequence was acquired at 7T: repetition time (TR)/echo time (TE)/ΔTE = 47/3/4 ms, 5 unipolar echoes, flip angle = 12°, field of view (FOV) = 224 × 180 × 90 mm3, 0.5 mm isotropic resolution, bandwidth = 502 Hz/voxel, SENSE factor of 2 × 2, scan time = 17 min 46 s per orientation. For STI, GRE scans using the same protocol were collected at 12 different brain orientations relative to the B0 field. For comparison with diffusion tensor imaging (DTI), diffusion MRI was acquired at 3T with a multi-shot EPI sequence: TR/TE=12500/81 ms, FOV =224 × 128 × 72 mm3, 1.5 mm isotropic resolution, EPI factor of 21, 2 shells with b-values to be 5000 and 10000 s/mm2 with 24 diffusion weighting directions per shell, scan time = 4 h 3 min. Diffusion MRI with reversed phase encoding (FH/HF) was also acquired.
Data processing
Phase images were preprocessed using best-path based phase unwrapping7, weighted averaging for echo combination8, and VSHARP for background field removal with maximum SMV kernel size of 6 mm9, 10. Preprocessed phase images at different orientations were nonlinearly co-registered to the reference position using ANTs11. Qualitative susceptibility mapping (QSM) using COSMOS12 was then calculated with reference to the whole brain susceptibitliy. MMS, MSA, and MSA-weighted PEV maps were calculated using asymmetric STI approach without explicit regularization13. Diffusion MRI was processed using Tortoise 4.014 and MRtrix315 to calculate mean diffusivity (MD), fractional anisotropy (FA), and FA-weighted PEV maps and co-registered to STI space.

RESULTS

Orientation dependence of GRE tissue frequency map (Figure 3)
For co-registered GRE scans, unit vectors of the B0 field in the subject frame of reference spanned uniformly across a half-spherical surface (Figure 3A). Tissue frequency maps demonstrated large orientation dependence across the brain, especially around the basal ganglia and the major white matter fiber bundles (Figure 3B).
Anatomical identification on QSM (Figure 4)
The lines of Gennari within the visual cortex around the calcarine fissure (Figure 4A), cortical substructures in the hippocampus and entorhinal cortex (Figure 4B), and detailed structures around the basal ganglia could be visualized (Figure 4C).
Comparison between STI and DTI (Figure 5)
The cortical substructures in the medial temporal lobe could be seen in the high-resolution MMS (Figure 5A). The MSA-weighted PEV in STI showed limited consistency with the FA-weighted PEV in DTI (Figure 5B).

DISCUSSION

Our ex vivo human brain STI protocol was useful for the visualization of microstructural neuroanatomy as well as the accumulation of postmortem human GRE data for future training and validation of data-driven STI reconstructions with sparse orientation sampling16, 17. Discrepancies between the PEV derived from STI and that from DTI might partly come from the enhanced chemical exchange effect and reduced susceptibility anisotropy due to fixation18.

CONCLUSION

We have demonstrated a feasible ex vivo STI protocol to scan human hemibrain at multiple orientations. Fine microstructural details in the neuroanatomy of the human brain can be visualized using the obtained QSM and STI data.

Acknowledgements

This work is supported by NIH NIBIB (P41EB031771), the Richman Family Precision Medicine Center of Excellence in Alzheimer's Disease including significant contributions from the Richman Family Foundation, the Rick Sharp Alzheimer’s Foundation, the Sharp Family Foundation and others. Kenichi Oishi is a consultant for “AnatomyWorks” and “Corporate-M.” Peter van Zijl has research support from and technology licensed to Philips Healthcare and has also been a paid speaker. This arrangement is being managed by the Johns Hopkins University in accordance with its conflict-of-interest policies. The authors thank the Brain Resource Center for providing the brain specimens and Mary McAllister for English-language editing.

References

References

1. Liu C. Susceptibility tensor imaging. Magn Reson Med 2010;63:1471-1477.

2. Li W, Wu B, Avram AV, Liu C. Magnetic susceptibility anisotropy of human brain in vivo and its molecular underpinnings. Neuroimage 2012;59:2088-2097.

3. Li X, Vikram DS, Lim IA, et al. Mapping magnetic susceptibility anisotropies of white matter in vivo in the human brain at 7 T. Neuroimage 2012;62:314-330.

4. Liu C, Li W, Wu B, et al. 3D fiber tractography with susceptibility tensor imaging. Neuroimage 2012;59:1290-1298.

5. Li W, Liu C, Duong TQ, et al. Susceptibility tensor imaging (STI) of the brain. NMR Biomed 2017;30.

6. Dawe RJ, Bennett DA, Schneider JA, et al. Postmortem MRI of human brain hemispheres: T2 relaxation times during formaldehyde fixation. Magn Reson Med 2009;61:810-818.

7. Shatil AS, Uddin MN, Matsuda KM, Figley CR. Quantitative Ex Vivo MRI Changes due to Progressive Formalin Fixation in Whole Human Brain Specimens: Longitudinal Characterization of Diffusion, Relaxometry, and Myelin Water Fraction Measurements at 3T. Front Med (Lausanne) 2018;5:31.

8. Abdul-Rahman HS, Gdeisat MA, Burton DR, et al. Fast and robust three-dimensional best path phase unwrapping algorithm. Appl Opt 2007;46:6623-6635.

9. Wu B, Li W, Avram AV, et al. Fast and tissue-optimized mapping of magnetic susceptibility and T2* with multi-echo and multi-shot spirals. Neuroimage 2012;59:297-305.

10. Schweser F, Deistung A, Lehr BW, Reichenbach JR. Quantitative imaging of intrinsic magnetic tissue properties using MRI signal phase: an approach to in vivo brain iron metabolism? Neuroimage 2011;54:2789-2807.

11. Wu B, Li W, Guidon A, Liu C. Whole brain susceptibility mapping using compressed sensing. Magn Reson Med 2012;67:137-147.

12. Tustison NJ, Cook PA, Holbrook AJ, et al. The ANTsX ecosystem for quantitative biological and medical imaging. Sci Rep 2021;11:9068.

13. Liu T, Spincemaille P, de Rochefort L, et al. Calculation of susceptibility through multiple orientation sampling (COSMOS): a method for conditioning the inverse problem from measured magnetic field map to susceptibility source image in MRI. Magn Reson Med 2009;61:196-204.

14. Cao S, Wei H, Chen J, Liu C. Asymmetric susceptibility tensor imaging. Magn Reson Med 2021;86:2266-2275.

15. Irfanoglu MO, Modi P, Nayak A, et al. DR-BUDDI (Diffeomorphic Registration for Blip-Up blip-Down Diffusion Imaging) method for correcting echo planar imaging distortions. Neuroimage 2015;106:284-299.

16. Tournier JD, Smith R, Raffelt D, et al. MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation. Neuroimage 2019;202:116137.

17. Bao L, Xiong C, Wei W, et al. Diffusion-regularized susceptibility tensor imaging (DRSTI) of tissue microstructures in the human brain. Med Image Anal 2021;67:101827.

18. Fang Z, Lai KW, van Zijl P, et al. DeepSTI: Towards tensor reconstruction using fewer orientations in susceptibility tensor imaging. Med Image Anal 2023;87:102829.

19. Chan KS, Hédouin R, Mollink J, et al. Imaging white matter microstructure with gradient-echo phase imaging: Is ex vivo imaging with formalin-fixed tissue a good approximation of the in vivo brain? Magn Reson Med 2022;88:380-390.

Figures

Figure 1: The blueprint for the container and the base. (A) The whole container. (B) The inner container. (C) The base of 0 degrees. (D) The base of 30 degrees. (E) The base of 60 degrees. (F) The base of 90 degrees.

Figure 2: A set of procedures for brain sample preparation from formalin fixation to MRI scanning. (A) The brain sample in phosphate-buffered saline two days before ex vivo MRI scans. (B) Black arrows indicate the top holes from which Fomblin was poured to the container. (C) The inner lid with plastic pipes was attached to the containe. (D) Fomblin was gently poured from one of the top holes. (E) The outer lid was fitted with spray adhesives. (F) Vacuum chamber for more than 12 h. (G) The degree of the angles was marked on the container. (H) The rotatable container with a base of 0 degrees.

Figure 3. Orientation dependence of GRE tissue frequency map. (A) GRE scans were acquired at 12 different brain orientations relative to the B0 field, with unit vectors of the main field (color-coded lines) in the subject frame of reference spanning uniformly across a half-spherical surface. (B) Tissue frequency maps were largely influenced by the sample orientation with respect to the field showing orientation-dependent contrasts across the brain, e.g., the basal ganglia (Yellow arrow) and optic radiation (Yellow arrowhead).

Detailed anatomical structures on QSM. (A) Yellow arrowhead indicates iron deposition on the line of Gennari within the visual cortex around the calcarine fissure. (B) Yellow arrowhead indicates the dentate gyrus of the hippocampus and yellow arrow indicates the entorhinal layer II islands in the entorhinal cortex. (C) Detailed structures around the basal ganglia are visualized. The red bounding boxes are 2.5 × magnified to clearly visualize these microstructures.

Figure 5. Comparison between STI and DTI. MMS, MSA, and PEV maps are displayed as STI metrics, whereas MD, FA, and FA-weighted PEV maps are displayed as DTI metrics. (A) The entorhinal cortex in the medial temporal lobe can be seen in MMS (Yellow arrowhead). (B) The optic radiation fiber appears smooth in FA-weighted PEV (Yellow arrow).

Proc. Intl. Soc. Mag. Reson. Med. 32 (2024)
4108
DOI: https://doi.org/10.58530/2024/4108