Diffusion imaging of whole, post-mortem human brains on a clinical MRI scanner

Advanced Magnetic Resonance Imaging Biomarkers of the Injured Spinal Cord: A Comparative Study of Imaging and Histology in Human Traumatic Spinal Cord Injury

A significant problem in the diagnosis and management of traumatic spinal cord injury (tSCI) is the heterogeneity of secondary injury and the prediction of neurological outcome. Imaging biomarkers specific to myelin loss and inflammation after tSCI would enable detailed assessment of the pathophysiological processes underpinning secondary damage to the cord. Such biomarkers could be used to biologically stratify injury severity and better inform prognosis for neurological recovery. While much work has been done to establish magnetic resonance imaging (MRI) biomarkers for SCI in animal models, the relationship between imaging findings and the underlying pathology has been difficult to discern in human tSCI because of the paucity of human spinal cord tissue. We utilized post-mortem spinal cords from individuals who had a tSCI to examine this relationship by performing ex vivo MRI scans before histological analysis. We investigated the correlation between the histological distribution of myelin loss and inflammatory cells in the injured spinal cord and a number of myelin and inflammation-sensitive MRI measures: myelin water fraction (MWF), inhomogeneous magnetization transfer ratio (ihMTR), and diffusion tensor and diffusion kurtosis imaging-derived fractional anisotropy (FA) and axial, radial, and mean diffusivity (AD, RD, MD). The histological features were analyzed by staining with Luxol Fast Blue (LFB) for myelin lipids and Class II major histocompatibility complex (Class II MHC) and CD68 for microglia and macrophages. Both MWF and ihMTR were strongly correlated with LFB staining for myelin, supporting the use of both as biomarkers for myelin loss after SCI. A decrease in ihMTR was also correlated with the presence of Class II MHC positive immune cells. FA and RD correlated with both Class II MHC and CD68 and may therefore be useful biomarkers for inflammation after tSCI. Our work demonstrates the utility of advanced MRI techniques sensitive to biological tissue damage after tSCI, which is an important step toward using these MRI techniques in the clinic to aid in decision-making.

Toward an optimal cadaveric brain model for neurosurgical education: assessment of preservation, parenchyma, vascular injection, and imaging

OBJECTIVE

We assessed types of cadaveric head and brain tissue specimen preparations that are used in a high throughput neurosurgical research laboratory to determine optimal preparation methods for neurosurgical anatomical research, education, and training.

METHODS

Cadaveric specimens (N = 112) prepared using different preservation and vascular injection methods were imaged, dissected, and graded by 11 neurosurgeons using a 21-point scale. We assessed the quality of tissue and preservation in both the anterior and posterior circulations. Tissue quality was evaluated using a 9-point magnetic resonance imaging (MRI) scale.

RESULTS

Formalin-fixed specimens yielded the highest scores for assessment (mean ± SD [17.0 ± 2.8]) vs. formalin-flushed (17.0 ± 3.6) and MRI (6.9 ± 2.0). Cadaver assessment and MRI scores were positively correlated (P < 0.001, R2 0.60). Analysis showed significant associations between cadaver assessment scores and specific variables: nonformalin fixation (β = -3.3), preservation within ≤72 h of death (β = 1.8), and MRI quality score (β = 0.7). Formalin-fixed specimens exhibited greater hardness than formalin-flushed and nonformalin-fixed specimens (P ≤ 0.006). Neurosurgeons preferred formalin-flushed specimens injected with colored latex.

CONCLUSION

For better-quality specimens for neurosurgical education and training, formalin preservation within ≤72 h of death was preferable, as was injection with colored latex. Formalin-flushed specimens more closely resembled live brain parenchyma. Assessment scores were lower for preparation techniques performed > 72 h postmortem and for nonformalin preservation solutions. The positive correlation between cadaver assessment scores and our novel MRI score indicates that donation organizations and institutional buyers should incorporate MRI as a screening tool for the selection of high-quality specimens.

Romer-EPTI: rotating-view motion-robust super-resolution EPTI for SNR-efficient distortion-free in-vivo mesoscale dMRI and microstructure imaging

Purpose

To overcome the major challenges in dMRI acquisition, including low SNR, distortion/blurring, and motion vulnerability.

Methods

A novel Romer-EPTI technique is developed to provide distortion-free dMRI with significant SNR gain, high motion-robustness, sharp spatial resolution, and simultaneous multi-TE imaging. It introduces a ROtating-view Motion-robust supEr-Resolution technique (Romer) combined with a distortion/blurring-free EPTI encoding. Romer enhances SNR by a simultaneous multi-thick-slice acquisition with rotating-view encoding, while providing high motion-robustness through a motion-aware super-resolution reconstruction, which also incorporates slice-profile and real-value diffusion, to resolve high-isotropic-resolution volumes. The in-plane encoding is performed using distortion/blurring-free EPTI, which further improves effective spatial resolution and motion robustness by preventing not only T2/T2*-blurring but also additional blurring resulting from combining encoded volumes with inconsistent geometries caused by dynamic distortions. Self-navigation was incorporated to enable efficient phase correction. Additional developments include strategies to address slab-boundary artifacts, achieve minimal TE for SNR gain at 7T, and achieve high robustness to strong phase variations at high b-values.

Results

Using Romer-EPTI, we demonstrate distortion-free whole-brain mesoscale in-vivo dMRI at both 3T (500-μm-iso) and 7T (485-μm-iso) for the first time, with high SNR efficiency (e.g., 25 × ), and high image quality free from distortion and slab-boundary artifacts with minimal blurring. Motion experiments demonstrate Romer-EPTI's high motion-robustness and ability to recover sharp images in the presence of motion. Romer-EPTI also demonstrates significant SNR gain and robustness in high b-value (b=5000s/mm2) and time-dependent dMRI.

Conclusion

Romer-EPTI significantly improves SNR, motion-robustness, and image quality, providing a highly efficient acquisition for high-resolution dMRI and microstructure imaging.

A cost-effective and efficient ex vivo, ex situ human whole brain perfusion protocol for immunohistochemistry

BACKGROUND

Chemical fixation of the brain can be executed through either the immersion method or the perfusion method. Perfusion fixation allows for better preservation of the brain tissue's ultrastructure, as it provides rapid and uniform delivery of the fixative to the tissue. Still, not all facilities have the expertise to perform perfusion fixation, with initial high cost and complexity of perfusion systems as the main factors limiting its widespread usage.

NEW METHOD

Here we present our low-cost approach of whole brain ex situ perfusion fixation to overcome the aforementioned limitations. Our self-made perfusion system, constructed utilising commercially accessible and affordable medical resources alongside laboratory and everyday items, demonstrates the capability to generate superior histological stainings of brain tissue. The perfused tissue can be stored prior to proceeding with IHC for at least one year.

RESULTS

Our method yielded high-quality results in histological stainings using both free-floating cryosections and paraffin-embedded tissue sections. The system is fully reusable and complies with the principles of sustainable management.

COMPARISON WITH EXISTING METHODS

Our whole brain perfusion system has been assembled from simple components and is able to achieve a linear flow with a pressure of 70 mmHg corresponding to the perfusion pressure of the brain.

CONCLUSIONS

Our ex situ method can be especially useful in research settings where expensive perfusion systems are not affordable or in any field with high time pressure, making it suitable for the field of forensic medicine or pathology in general.

Eddy current-induced artifact correction in high b-value ex vivo human brain diffusion MRI with dynamic field monitoring

PURPOSE

To investigate whether spatiotemporal magnetic field monitoring can correct pronounced eddy current-induced artifacts incurred by strong diffusion-sensitizing gradients up to 300 mT/m used in high b-value diffusion-weighted (DW) EPI.

METHODS

A dynamic field camera equipped with 16 1 H NMR field probes was first used to characterize field perturbations caused by residual eddy currents from diffusion gradients waveforms in a 3D multi-shot EPI sequence on a 3T Connectom scanner for different gradient strengths (up to 300 mT/m), diffusion directions, and shots. The efficacy of dynamic field monitoring-based image reconstruction was demonstrated on high-gradient strength, submillimeter resolution whole-brain ex vivo diffusion MRI. A 3D multi-shot image reconstruction framework was developed that incorporated the nonlinear phase evolution measured with the dynamic field camera.

RESULTS

Phase perturbations in the readout induced by residual eddy currents from strong diffusion gradients are highly nonlinear in space and time, vary among diffusion directions, and interfere significantly with the image encoding gradients, changing the k-space trajectory. During the readout, phase modulations between odd and even EPI echoes become non-static and diffusion encoding direction-dependent. Superior reduction of ghosting and geometric distortion was achieved with dynamic field monitoring compared to ghosting reduction approaches such as navigator- and structured low-rank-based methods or MUSE followed by image-based distortion correction with the FSL tool "eddy."

CONCLUSION

Strong eddy current artifacts characteristic of high-gradient strength DW-EPI can be well corrected with dynamic field monitoring-based image reconstruction.

Fluid preservation in brain banking: a review

Fluid preservation is nearly universally used in brain banking to store fixed tissue specimens for future research applications. However, the effects of long-term immersion on neural circuitry and biomolecules are not well characterized. As a result, there is a need to synthesize studies investigating fluid preservation of brain tissue. We searched PubMed and other databases to identify studies measuring the effects of fluid preservation in nervous system tissue. We categorized studies based on the fluid preservative used: formaldehyde solutions, buffer solutions, alcohol solutions, storage after tissue clearing, and cryoprotectant solutions. We identified 91 studies containing 197 independent observations of the effects of long-term storage on cellular morphology. Most studies did not report any significant alterations due to long-term storage. When present, the most frequent alteration was decreased antigenicity, commonly attributed to progressive crosslinking by aldehydes that renders biomolecules increasingly inaccessible over time. To build a mechanistic understanding, we discuss biochemical aspects of long-term fluid preservation. A subset of lipids appears to be chemical altered or extracted over time due to incomplete retention in the crosslinked gel. Alternative storage fluids mitigate the problem of antigen masking but have not been extensively characterized and may have other downsides. We also compare fluid preservation to cryopreservation, paraffin embedding, and resin embedding. Overall, existing evidence suggests that fluid preservation provides maintenance of neural architecture for decades, including precise structural details. However, to avoid the well-established problem of overfixation caused by storage in high concentration formaldehyde solutions, fluid preservation procedures can use an initial fixation step followed by an alternative long-term storage fluid. Further research is warranted on optimizing protocols and characterizing the generalizability of the storage artifacts that have been identified.

The Brainbox-a tool to facilitate correlation of brain magnetic resonance imaging features to histopathology

Magnetic resonance imaging (MRI) has limitations in identifying underlying tissue pathology, which is relevant for neurological diseases such as multiple sclerosis, stroke or brain tumours. However, there are no standardized methods for correlating MRI features with histopathology. Thus, here we aimed to develop and validate a tool that can facilitate the correlation of brain MRI features to corresponding histopathology. For this, we designed the Brainbox, a waterproof and MRI-compatible 3D printed container with an integrated 3D coordinate system. We used the Brainbox to acquire post-mortem ex vivo MRI of eight human brains, fresh and formalin-fixed, and correlated focal imaging features to histopathology using the built-in 3D coordinate system. With its built-in 3D coordinate system, the Brainbox allowed correlation of MRI features to corresponding tissue substrates. The Brainbox was used to correlate different MR image features of interest to the respective tissue substrate, including normal anatomical structures such as the hippocampus or perivascular spaces, as well as a lacunar stroke. Brain volume decreased upon fixation by 7% (P = 0.01). The Brainbox enabled degassing of specimens before scanning, reducing susceptibility artefacts and minimizing bulk motion during scanning. In conclusion, our proof-of-principle experiments demonstrate the usability of the Brainbox, which can contribute to improving the specificity of MRI and the standardization of the correlation between post-mortem ex vivo human brain MRI and histopathology. Brainboxes are available upon request from our institution.

A cellular resolution atlas of Broca's area

Brain cells are arranged in laminar, nuclear, or columnar structures, spanning a range of scales. Here, we construct a reliable cell census in the frontal lobe of human cerebral cortex at micrometer resolution in a magnetic resonance imaging (MRI)-referenced system using innovative imaging and analysis methodologies. MRI establishes a macroscopic reference coordinate system of laminar and cytoarchitectural boundaries. Cell counting is obtained with a digital stereological approach on the 3D reconstruction at cellular resolution from a custom-made inverted confocal light-sheet fluorescence microscope (LSFM). Mesoscale optical coherence tomography enables the registration of the distorted histological cell typing obtained with LSFM to the MRI-based atlas coordinate system. The outcome is an integrated high-resolution cellular census of Broca's area in a human postmortem specimen, within a whole-brain reference space atlas.

Sampling strategies and integrated reconstruction for reducing distortion and boundary slice aliasing in high-resolution 3D diffusion MRI

PURPOSE

To develop a new method for high-fidelity, high-resolution 3D multi-slab diffusion MRI with minimal distortion and boundary slice aliasing.

METHODS

Our method modifies 3D multi-slab imaging to integrate blip-reversed acquisitions for distortion correction and oversampling in the slice direction (kz ) for reducing boundary slice aliasing. Our aim is to achieve robust acceleration to keep the scan time the same as conventional 3D multi-slab acquisitions, in which data are acquired with a single direction of blip traversal and without kz -oversampling. We employ a two-stage reconstruction. In the first stage, the blip-up/down images are respectively reconstructed and analyzed to produce a field map for each diffusion direction. In the second stage, the blip-reversed data and the field map are incorporated into a joint reconstruction to produce images that are corrected for distortion and boundary slice aliasing.

RESULTS

We conducted experiments at 7T in six healthy subjects. Stage 1 reconstruction produces images from highly under-sampled data (R = 7.2) with sufficient quality to provide accurate field map estimation. Stage 2 joint reconstruction substantially reduces distortion artifacts with comparable quality to fully-sampled blip-reversed results (2.4× scan time). Whole-brain in-vivo results acquired at 1.22 mm and 1.05 mm isotropic resolutions demonstrate improved anatomical fidelity compared to conventional 3D multi-slab imaging. Data demonstrate good reliability and reproducibility of the proposed method over multiple subjects.

CONCLUSION

The proposed acquisition and reconstruction framework provide major reductions in distortion and boundary slice aliasing for 3D multi-slab diffusion MRI without increasing the scan time, which can potentially produce high-quality, high-resolution diffusion MRI.

Three-dimensional diffusion MRI using simultaneous multislab with blipped-CAIPI in a 4D k-space framework

PURPOSE

To develop an efficient simultaneous multislab imaging method with blipped-controlled aliasing in parallel imaging (blipped-SMSlab) in a 4D k-space framework, and to demonstrate its efficacy in high-resolution diffusion MRI (dMRI).

THEORY AND METHODS

First, the SMSlab 4D k-space signal expression is formulated, and the phase interferences from intraslab and interslab encodings on the same physical z-axis are analyzed. Then, the blipped-SMSlab dMRI sequence is designed, with blipped-controlled aliasing in parallel imaging (blipped-CAIPI) gradients for interslab encoding, and a 2D multiband accelerated navigator for inter-kz-shot phase correction. Third, strategies are developed to remove the phase interferences, by RF phase modulation and/or phase correction during reconstruction, thus decoupling intraslab and interslab encodings that are otherwise entangled. In vivo experiments are performed to validate the blipped-SMSlab method and preliminarily evaluate its performance in high-resolution dMRI compared with traditional 2D imaging.

RESULTS

In the 4D k-space framework, interslab and intraslab phase interferences of blipped-SMSlab are successfully removed using the proposed strategies. Compared with non-CAIPI sampling, the blipped-SMSlab acquisition reduces the g-factor and g-factor-related SNR penalty by about 12%. In addition, in vivo experiments show the SNR advantage of blipped-SMSlab dMRI over traditional 2D dMRI for 1.3-mm and 1.0-mm isotropic resolution imaging with matched acquisition time.

CONCLUSION

Removing interslab and intraslab phase interferences enables SMSlab dMRI with blipped-CAIPI in a 4D k-space framework. The proposed blipped-SMSlab dMRI is demonstrated to be more SNR-efficient than 2D dMRI and thus capable of high-quality, high-resolution fiber orientation detection.

Nigral Pathology Contributes to Microstructural Integrity of Striatal and Frontal Tracts in Parkinson's Disease

BACKGROUND

Motor and cognitive impairment in Parkinson's disease (PD) is associated with dopaminergic dysfunction that stems from substantia nigra (SN) degeneration and concomitant α-synuclein accumulation. Diffusion magnetic resonance imaging (MRI) can detect microstructural alterations of the SN and its tracts to (sub)cortical regions, but their pathological sensitivity is still poorly understood.

OBJECTIVE

To unravel the pathological substrate(s) underlying microstructural alterations of SN, and its tracts to the dorsal striatum and dorsolateral prefrontal cortex (DLPFC) in PD.

METHODS

Combining post-mortem in situ MRI and histopathology, T1-weighted and diffusion MRI, and neuropathological samples of nine PD, six PD with dementia (PDD), five dementia with Lewy bodies (DLB), and 10 control donors were collected. From diffusion MRI, mean diffusivity (MD) and fractional anisotropy (FA) were derived from the SN, and tracts between the SN and caudate nucleus, putamen, and DLPFC. Phosphorylated-Ser129-α-synuclein and tyrosine hydroxylase immunohistochemistry was included to quantify nigral Lewy pathology and dopaminergic degeneration, respectively.

RESULTS

Compared to controls, PD and PDD/DLB showed increased MD of the SN and SN-DLPFC tract, as well as increased FA of the SN-caudate nucleus tract. Both PD and PDD/DLB showed nigral Lewy pathology and dopaminergic loss compared to controls. Increased MD of the SN and FA of SN-caudate nucleus tract were associated with SN dopaminergic loss. Whereas increased MD of the SN-DLPFC tract was associated with increased SN Lewy neurite load.

CONCLUSIONS

In PD and PDD/DLB, diffusion MRI captures microstructural alterations of the SN and tracts to the dorsal striatum and DLPFC, which differentially associates with SN dopaminergic degeneration and Lewy neurite pathology. © 2023 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

The acute effects of repetitive transcranial magnetic stimulation on laminar diffusion anisotropy of neocortical gray matter

Repetitive transcranial magnetic stimulation (rTMS) is increasingly used to treat neuropsychiatric disorders. Inhibitory and excitatory regimens have been both adopted but the exact mechanism of action remains unclear, and investigating their differential effects on laminar diffusion profiles of neocortex may add important evidence. Twenty healthy participants were randomly assigned to receive a low-frequency/inhibitory or high-frequency/excitatory rTMS targeting the left dorsolateral prefrontal cortex (DLPFC). With the brand-new submillimeter diffusion tensor imaging of whole brain and specialized surface-based laminar analysis, fractional anisotropy (FA) and mean diffusion (MD) profiles of cortical layers at different cortical depths were characterized before/after rTMS. Inhibitory and excitatory rTMS both showed impacts on diffusion metrics of somatosensory, limbic, and sensory regions, but different patterns of changes were observed-increased FA with inhibitory rTMS, whereas decreased FA with excitatory rTMS. More importantly, laminar analysis indicated laminar specificity of changes in somatosensory regions during different rTMS patterns-inhibitory rTMS affected the superficial layers contralateral to the DLPFC, while excitatory rTMS led to changes in the intermediate/deep layers bilateral to the DLPFC. These findings provide novel insights into acute neurobiological effects on diffusion profiles of rTMS that may add critical evidence relevant to different protocols of rTMS on neocortex.

High angular resolution susceptibility imaging and estimation of fiber orientation distribution functions in primate brain

Uncovering brain-tissue microstructure including axonal characteristics is a major neuroimaging research focus. Within this scope, anisotropic properties of magnetic susceptibility in white matter have been successfully employed to estimate primary axonal trajectories using mono-tensorial models. However, anisotropic susceptibility has not yet been considered for modeling more complex fiber structures within a voxel, such as intersecting bundles, or an estimation of orientation distribution functions (ODFs). This information is routinely obtained by high angular resolution diffusion imaging (HARDI) techniques. In applications to fixed tissue, however, diffusion-weighted imaging suffers from an inherently low signal-to-noise ratio and limited spatial resolution, leading to high demands on the performance of the gradient system in order to mitigate these limitations. In the current work, high angular resolution susceptibility imaging (HARSI) is proposed as a novel, phase-based methodology to estimate ODFs. A multiple gradient-echo dataset was acquired in an entire fixed chimpanzee brain at 61 orientations by reorienting the specimen in the magnetic field. The constant solid angle method was adapted for estimating phase-based ODFs. HARDI data were also acquired for comparison. HARSI yielded information on whole-brain fiber architecture, including identification of peaks of multiple bundles that resembled features of the HARDI results. Distinct differences between both methods suggest that susceptibility properties may offer complementary microstructural information. These proof-of-concept results indicate a potential to study the axonal organization in post-mortem primate and human brain at high resolution.

A Postmortem MRI Study of Cerebrovascular Disease and Iron Content at End-Stage of Fragile X-Associated Tremor/Ataxia Syndrome

Brain changes at the end-stage of fragile X-associated tremor/ataxia syndrome (FXTAS) are largely unknown due to mobility impairment. We conducted a postmortem MRI study of FXTAS to quantify cerebrovascular disease, brain atrophy and iron content, and examined their relationships using principal component analysis (PCA). Intracranial hemorrhage (ICH) was observed in 4/17 FXTAS cases, among which one was confirmed by histologic staining. Compared with seven control brains, FXTAS cases showed higher ratings of T2-hyperintensities (indicating cerebral small vessel disease) in the cerebellum, globus pallidus and frontoparietal white matter, and significant atrophy in the cerebellar white matter, red nucleus and dentate nucleus. PCA of FXTAS cases revealed negative associations of T2-hyperintensity ratings with anatomic volumes and iron content in the white matter, hippocampus and amygdala, that were independent from a highly correlated number of regions with ICH and iron content in subcortical nuclei. Post-hoc analysis confirmed PCA findings and further revealed increased iron content in the white matter, hippocampus and amygdala in FXTAS cases compared to controls, after adjusting for T2-hyperintensity ratings. These findings indicate that both ischemic and hemorrhagic brain damage may occur in FXTAS, with the former being marked by demyelination/iron depletion and atrophy, and the latter by ICH and iron accumulation in basal ganglia.

High-resolution multi-shot diffusion imaging of structural networks in healthy neurocognitive aging

Healthy neurocognitive aging has been associated with the microstructural degradation of white matter pathways that connect distributed gray matter regions, assessed by diffusion-weighted imaging (DWI). However, the relatively low spatial resolution of standard DWI has limited the examination of age-related differences in the properties of smaller, tightly curved white matter fibers, as well as the relatively more complex microstructure of gray matter. Here, we capitalize on high-resolution multi-shot DWI, which allows spatial resolutions < 1 mm3 to be achieved on clinical 3T MRI scanners. We assessed whether traditional diffusion tensor-based measures of gray matter microstructure and graph theoretical measures of white matter structural connectivity assessed by standard (1.5 mm3 voxels, 3.375 μl volume) and high-resolution (1 mm3 voxels, 1μl volume) DWI were differentially related to age and cognitive performance in 61 healthy adults 18-78 years of age. Cognitive performance was assessed using an extensive battery comprising 12 separate tests of fluid (speed-dependent) cognition. Results indicated that the high-resolution data had larger correlations between age and gray matter mean diffusivity, but smaller correlations between age and structural connectivity. Moreover, parallel mediation models including both standard and high-resolution measures revealed that only the high-resolution measures mediated age-related differences in fluid cognition. These results lay the groundwork for future studies planning to apply high-resolution DWI methodology to further assess the mechanisms of both healthy aging and cognitive impairment.

Mixed methodology in human brain research: integrating MRI and histology

Postmortem magnetic resonance imaging (MRI) can provide a bridge between histological observations and the in vivo anatomy of the human brain. Approaches aimed at the co-registration of data derived from the two techniques are gaining interest. Optimal integration of the two research fields requires detailed knowledge of the tissue property requirements for individual research techniques, as well as a detailed understanding of the consequences of tissue fixation steps on the imaging quality outcomes for both MRI and histology. Here, we provide an overview of existing studies that bridge between state-of-the-art imaging modalities, and discuss the background knowledge incorporated into the design, execution and interpretation of postmortem studies. A subset of the discussed challenges transfer to animal studies as well. This insight can contribute to furthering our understanding of the normal and diseased human brain, and to facilitate discussions between researchers from the individual disciplines.

Multi-tensor diffusion abnormalities of gray matter in an animal model of cortical dysplasia

Focal cortical dysplasias are a type of malformations of cortical development that are a common cause of drug-resistant focal epilepsy. Surgical treatment is a viable option for some of these patients, with their outcome being highly related to complete surgical resection of lesions visible in magnetic resonance imaging (MRI). However, subtle lesions often go undetected on conventional imaging. Several methods to analyze MRI have been proposed, with the common goal of rendering subtle cortical lesions visible. However, most image-processing methods are targeted to detect the macroscopic characteristics of cortical dysplasias, which do not always correspond to the microstructural disarrangement of these cortical malformations. Quantitative analysis of diffusion-weighted MRI (dMRI) enables the inference of tissue characteristics, and novel methods provide valuable microstructural features of complex tissue, including gray matter. We investigated the ability of advanced dMRI descriptors to detect diffusion abnormalities in an animal model of cortical dysplasia. For this purpose, we induced cortical dysplasia in 18 animals that were scanned at 30 postnatal days (along with 19 control animals). We obtained multi-shell dMRI, to which we fitted single and multi-tensor representations. Quantitative dMRI parameters derived from these methods were queried using a curvilinear coordinate system to sample the cortical mantle, providing inter-subject anatomical correspondence. We found region- and layer-specific diffusion abnormalities in experimental animals. Moreover, we were able to distinguish diffusion abnormalities related to altered intra-cortical tangential fibers from those associated with radial cortical fibers. Histological examinations revealed myelo-architectural abnormalities that explain the alterations observed through dMRI. The methods for dMRI acquisition and analysis used here are available in clinical settings and our work shows their clinical relevance to detect subtle cortical dysplasias through analysis of their microstructural properties.

Magnetic resonance imaging at 9.4 T: the Maastricht journey

The 9.4 T scanner in Maastricht is a whole-body magnet with head gradients and parallel RF transmit capability. At the time of the design, it was conceptualized to be one of the best fMRI scanners in the world, but it has also been used for anatomical and diffusion imaging. 9.4 T offers increases in sensitivity and contrast, but the technical ultra-high field (UHF) challenges, such as field inhomogeneities and constraints set by RF power deposition, are exacerbated compared to 7 T. This article reviews some of the 9.4 T work done in Maastricht. Functional imaging experiments included blood oxygenation level-dependent (BOLD) and blood-volume weighted (VASO) fMRI using different readouts. BOLD benefits from shorter T₂* at 9.4 T while VASO from longer T₁. We show examples of both ex vivo and in vivo anatomical imaging. For many applications, pTx and optimized coils are essential to harness the full potential of 9.4 T. Our experience shows that, while considerable effort was required compared to our 7 T scanner, we could obtain high-quality anatomical and functional data, which illustrates the potential of MR acquisitions at even higher field strengths. The practical challenges of working with a relatively unique system are also discussed.

Tissue optimization strategies for high-quality ex vivo diffusion imaging

Ex vivo diffusion imaging can be used to study healthy and pathological tissue microstructure in the rodent brain with high resolution, providing a link between in vivo MRI and ex vivo microscopy techniques. Major challenges for the successful acquisition of ex vivo diffusion imaging data however are changes in the relaxivity and diffusivity of brain tissue following perfusion fixation. In this study we address this question by examining the combined effects of tissue preparation factors that influence signal-to-noise ratio (SNR) and consequently image quality, including fixative concentration, contrast agent concentration and tissue rehydration time. We present an optimization strategy combining these factors to manipulate theT 1 andT 2 of fixed tissue and maximize SNR efficiency. We apply this strategy in the rat brain, for a diffusion-weighted spin echo protocol with TE = 27 ms on a 9.4 T scanner with a 39 mm volume coil and 660 mT/m 114 mm gradient insert. We used a reduced fixative concentration of 2% paraformaldehyde (PFA), rehydration time more than 20 days, 15 mM Gd-DTPA in perfusate and TR 250 ms. This resulted in a doubling of SNR and an increase in SNR per unit time of 135% in cortical grey matter and 88% in white matter compared with 4% PFA and no contrast agent. This improved SNR efficiency enabled the acquisition of excellent-quality high-resolution (78 μ m isotropic voxel size) diffusion data with b = 4000 s/mm2 , 30 diffusion directions and a field of view of 40 × 13 × 18 mm3 in less than 4 days. It was also possible to achieve comparable data quality for a standard resolution (150 μ m) diffusion dataset in 2 1 4 h. In conclusion, the tissue optimization strategy presented here may be used to improve SNR, increase spatial resolution and/or allow faster acquisitions in preclinical ex vivo diffusion MRI experiments.

Eddy current-induced artifacts correction in high gradient strength diffusion MRI with dynamic field monitoring: demonstration in ex vivo human brain imaging

Purpose

To demonstrate the advantages of spatiotemporal magnetic field monitoring to correct eddy current-induced artifacts (ghosting and geometric distortions) in high gradient strength diffusion MRI (dMRI).

Methods

A dynamic field camera with 16 NMR field probes was used to characterize eddy current fields induced from diffusion gradients for different gradients strengths (up to 300 mT/m), diffusion directions, and shots in a 3D multi-shot EPI sequence on a 3T Connectom scanner. The efficacy of dynamic field monitoring-based image reconstruction was demonstrated on high-resolution whole brain ex vivo dMRI. A 3D multi-shot image reconstruction framework was informed with the actual nonlinear phase evolution measured with the dynamic field camera, thereby accounting for high-order eddy currents fields on top of the image encoding gradients in the image formation model.

Results

Eddy current fields from diffusion gradients at high gradient strength in a 3T Connectom scanner are highly nonlinear in space and time, inducing high-order spatial phase modulations between odd/even echoes and shots that are not static during the readout. Superior reduction of ghosting and geometric distortion was achieved with dynamic field monitoring compared to ghosting approaches such as navigator- and structured low-rank-based methods or MUSE, followed by image-based distortion correction with eddy. Improved dMRI analysis is demonstrated with diffusion tensor imaging and high-angular resolution diffusion imaging.

Conclusion

Strong eddy current artifacts characteristic of high gradient strength dMRI can be well corrected with dynamic field monitoring-based image reconstruction, unlike the two-step approach consisting of ghosting correction followed by geometric distortion reduction with eddy.

High-fidelity, high-spatial-resolution diffusion magnetic resonance imaging of ex vivo whole human brain at ultra-high gradient strength with structured low-rank echo-planar imaging ghost correction

Diffusion magnetic resonance imaging (dMRI) of whole ex vivo human brain specimens enables three-dimensional (3D) mapping of structural connectivity at the mesoscopic scale, providing detailed evaluation of fiber architecture and tissue microstructure at a spatial resolution that is difficult to access in vivo. To account for the short T2 and low diffusivity of fixed tissue, ex vivo dMRI is often acquired using strong diffusion-sensitizing gradients and multishot/segmented 3D echo-planar imaging (EPI) sequences to achieve high spatial resolution. However, the combination of strong diffusion-sensitizing gradients and multishot/segmented EPI readout can result in pronounced ghosting artifacts incurred by nonlinear spatiotemporal variations in the magnetic field produced by eddy currents. Such ghosting artifacts cannot be corrected with conventional correction solutions and pose a significant roadblock to leveraging human MRI scanners with ultrahigh gradients for ex vivo whole-brain dMRI. Here, we show that ghosting-correction approaches that correct for either polarity-related ghosting or shot-to-shot variations in a separate manner are suboptimal for 3D multishot diffusion-weighted EPI experiments in fixed human brain specimens using strong diffusion-sensitizing gradients on the 3-T Connectom MRI scanner, resulting in orientationally biased dMRI estimates. We apply a recently developed advanced k-space reconstruction method based on structured low-rank matrix (SLM) modeling that handles both polarity-related ghosting and shot-to-shot variation simultaneously, to mitigate artifacts in high-angular resolution multishot dMRI data acquired in several fixed human brain specimens at 0.7-0.8-mm isotropic spatial resolution using b-values up to 10,000 s/mm2 and gradient strengths up to 280 mT/m. We demonstrate the improved mapping of diffusion tensor imaging and fiber orientation distribution functions in key neuroanatomical areas distributed across the whole brain using SLM-based EPI ghost correction compared with alternative techniques.

A postmortem MRI study of cerebrovascular disease and iron content at end-stage of fragile X-associated tremor/ataxia syndrome

Brain changes at end-stage of fragile X-associated tremor/ataxia syndrome (FXTAS) are largely unknown due to mobility impairment. We conducted a postmortem MRI study of FXTAS to quantify cerebrovascular disease, brain atrophy, and iron content and examined their relationships using principal component analysis (PCA). Intracranial hemorrhage (ICH) was observed in 4/17 FXTAS cases among which one was confirmed by histologic staining. Compared with seven control brains, FXTAS cases showed higher ratings of T2-hyperintensities (indicating cerebral small vessel disease) in the cerebellum, globus pallidus, and frontoparietal white matter and significant atrophy in cerebellar white matter, red nucleus, and dentate nucleus. PCA of FXTAS cases revealed negative associations of T2-hyperintensity ratings with anatomic volumes and iron content in the white matter, hippocampus, and amygdala, that were independent from highly correlated number of regions with ICH and iron content in subcortical nuclei. Post hoc analysis confirmed PCA findings and further revealed increased iron content in the white matter, hippocampus, and amygdala in FXTAS cases than controls after adjusting for T2-hyperintensity ratings. These findings indicate that both ischemic and hemorrhagic brain damage may occur in FXTAS, with the former marked by demyelination/iron depletion and atrophy and the latter, ICH and iron accumulation in basal ganglia.

Tensor image registration library: Deformable registration of stand-alone histology images to whole-brain post-mortem MRI data

BACKGROUND

Accurate registration between microscopy and MRI data is necessary for validating imaging biomarkers against neuropathology, and to disentangle complex signal dependencies in microstructural MRI. Existing registration methods often rely on serial histological sampling or significant manual input, providing limited scope to work with a large number of stand-alone histology sections. Here we present a customisable pipeline to assist the registration of stand-alone histology sections to whole-brain MRI data.

METHODS

Our pipeline registers stained histology sections to whole-brain post-mortem MRI in 4 stages, with the help of two photographic intermediaries: a block face image (to undistort histology sections) and coronal brain slab photographs (to insert them into MRI space). Each registration stage is implemented as a configurable stand-alone Python script using our novel platform, Tensor Image Registration Library (TIRL), which provides flexibility for wider adaptation. We report our experience of registering 87 PLP-stained histology sections from 14 subjects and perform various experiments to assess the accuracy and robustness of each stage of the pipeline.

RESULTS

All 87 histology sections were successfully registered to MRI. Histology-to-block registration (Stage 1) achieved 0.2-0.4 mm accuracy, better than commonly used existing methods. Block-to-slice matching (Stage 2) showed great robustness in automatically identifying and inserting small tissue blocks into whole brain slices with 0.2 mm accuracy. Simulations demonstrated sub-voxel level accuracy (0.13 mm) of the slice-to-volume registration (Stage 3) algorithm, which was observed in over 200 actual brain slice registrations, compensating 3D slice deformations up to 6.5 mm. Stage 4 combined the previous stages and generated refined pixelwise aligned multi-modal histology-MRI stacks.

CONCLUSIONS

Our open-source pipeline provides robust automation tools for registering stand-alone histology sections to MRI data with sub-voxel level precision, and the underlying framework makes it readily adaptable to a diverse range of microscopy-MRI studies.

COnstrained Reference frame diffusion TEnsor Correlation Spectroscopic (CORTECS) MRI: A practical framework for high-resolution diffusion tensor distribution imaging

High-resolution imaging studies have consistently shown that in cortical tissue water diffuses preferentially along radial and tangential orientations with respect to the cortical surface, in agreement with histology. These dominant orientations do not change significantly even if the relative contributions from microscopic water pools to the net voxel signal vary across experiments that use different diffusion times, b-values, TEs, and TRs. With this in mind, we propose a practical new framework for imaging non-parametric diffusion tensor distributions (DTDs) by constraining the microscopic diffusion tensors of the DTD to be diagonalized using the same orthonormal reference frame of the mesoscopic voxel. In each voxel, the constrained DTD (cDTD) is completely determined by the correlation spectrum of the microscopic principal diffusivities associated with the axes of the voxel reference frame. Consequently, all cDTDs are inherently limited to the domain of positive definite tensors and can be reconstructed efficiently using Inverse Laplace Transform methods. Moreover, the cDTD reconstruction can be performed using only data acquired efficiently with single diffusion encoding, although it also supports datasets with multiple diffusion encoding. In tissues with a well-defined architecture, such as the cortex, we can further constrain the cDTD to contain only cylindrically symmetric diffusion tensors and measure the 2D correlation spectra of principal diffusivities along the radial and tangential orientation with respect to the cortical surface. To demonstrate this framework, we perform numerical simulations and analyze high-resolution dMRI data from a fixed macaque monkey brain. We estimate 2D cDTDs in the cortex and derive, in each voxel, the marginal distributions of the microscopic principal diffusivities, the corresponding distributions of the microscopic fractional anisotropies and mean diffusivities along with their 2D correlation spectra to quantify the cDTD shape-size characteristics. Signal components corresponding to specific bands in these cDTD-derived spectra show high specificity to cortical laminar structures observed with histology. Our framework drastically simplifies the measurement of non-parametric DTDs in high-resolution datasets with mesoscopic voxel sizes much smaller than the radius of curvature of the underlying anatomy, e.g., cortical surface, and can be applied retrospectively to analyze existing diffusion MRI data from fixed cortical tissues.

High-resolution cortical MAP-MRI reveals areal borders and laminar substructures observed with histological staining

The variations in cellular composition and tissue architecture measured with histology provide the biological basis for partitioning the brain into distinct cytoarchitectonic areas and for characterizing neuropathological tissue alterations. Clearly, there is an urgent need to develop whole-brain neuroradiological methods that can assess cortical cyto- and myeloarchitectonic features non-invasively. Mean apparent propagator (MAP) MRI is a clinically feasible diffusion MRI method that quantifies efficiently and comprehensively the net microscopic displacements of water molecules diffusing in tissues. We investigate the sensitivity of high-resolution MAP-MRI to detecting areal and laminar variations in cortical cytoarchitecture and compare our results with observations from corresponding histological sections in the entire brain of a rhesus macaque monkey. High-resolution images of MAP-derived parameters, in particular the propagator anisotropy (PA), non-gaussianity (NG), and the return-to-axis probability (RTAP) reveal cortical area-specific lamination patterns in good agreement with the corresponding histological stained sections. In a few regions, the MAP parameters provide superior contrast to the five histological stains used in this study, delineating more clearly boundaries and transition regions between cortical areas and laminar substructures. Throughout the cortex, various MAP parameters can be used to delineate transition regions between specific cortical areas observed with histology and to refine areal boundaries estimated using atlas registration-based cortical parcellation. Using surface-based analysis of MAP parameters we quantify the cortical depth dependence of diffusion propagators in multiple regions-of-interest in a consistent and rigorous manner that is largely independent of the cortical folding geometry. The ability to assess cortical cytoarchitectonic features efficiently and non-invasively, its clinical feasibility, and translatability make high-resolution MAP-MRI a promising 3D imaging tool for studying whole-brain cortical organization, characterizing abnormal cortical development, improving early diagnosis of neurodegenerative diseases, identifying targets for biopsies, and complementing neuropathological investigations.

Diffusion MRI anisotropy in the cerebral cortex is determined by unmyelinated tissue features

Diffusion magnetic resonance imaging (dMRI) is commonly used to assess the tissue and cellular substructure of the human brain. In the white matter, myelinated axons are the principal neural elements that shape dMRI through the restriction of water diffusion; however, in the gray matter the relative contributions of myelinated axons and other tissue features to dMRI are poorly understood. Here we investigate the determinants of diffusion in the cerebral cortex. Specifically, we ask whether myelinated axons significantly shape dMRI fractional anisotropy (dMRI-FA), a measure commonly used to characterize tissue properties in humans. We compared ultra-high resolution ex vivo dMRI data from the brain of a marmoset monkey with both myelin- and Nissl-stained histological sections obtained from the same brain after scanning. We found that the dMRI-FA did not match the spatial distribution of myelin in the gray matter. Instead dMRI-FA was more closely related to the anisotropy of stained tissue features, most prominently those revealed by Nissl staining and to a lesser extent those revealed by myelin staining. Our results suggest that unmyelinated neurites such as large caliber apical dendrites are the primary features shaping dMRI measures in the cerebral cortex.

Post mortem mapping of connectional anatomy for the validation of diffusion MRI

Diffusion MRI (dMRI) is a unique tool for the study of brain circuitry, as it allows us to image both the macroscopic trajectories and the microstructural properties of axon bundles in vivo. The Human Connectome Project ushered in an era of impressive advances in dMRI acquisition and analysis. As a result of these efforts, the quality of dMRI data that could be acquired in vivo improved substantially, and large collections of such data became widely available. Despite this progress, the main limitation of dMRI remains: it does not image axons directly, but only provides indirect measurements based on the diffusion of water molecules. Thus, it must be validated by methods that allow direct visualization of axons but that can only be performed in post mortem brain tissue. In this review, we discuss methods for validating the various features of connectional anatomy that are extracted from dMRI, both at the macro-scale (trajectories of axon bundles), and at micro-scale (axonal orientations and other microstructural properties). We present a range of validation tools, including anatomic tracer studies, Klingler's dissection, myelin stains, label-free optical imaging techniques, and others. We provide an overview of the basic principles of each technique, its limitations, and what it has taught us so far about the accuracy of different dMRI acquisition and analysis approaches.

High-resolution neural network-driven mapping of multiple diffusion metrics leveraging asymmetries in the balanced steady-state free precession frequency profile

We propose to utilize the rich information content about microstructural tissue properties entangled in asymmetric balanced steady-state free precession (bSSFP) profiles to estimate multiple diffusion metrics simultaneously by neural network (NN) parameter quantification. A 12-point bSSFP phase-cycling scheme with high-resolution whole-brain coverage is employed at 3 and 9.4 T for NN input. Low-resolution target diffusion data are derived based on diffusion-weighted spin-echo echo-planar-imaging (SE-EPI) scans, that is, mean, axial, and radial diffusivity (MD, AD, and RD), fractional anisotropy (FA), as well as the spherical coordinates (azimuth Φ and inclination ϴ) of the principal diffusion eigenvector. A feedforward NN is trained with incorporated probabilistic uncertainty estimation. The NN predictions yielded highly reliable results in white matter (WM) and gray matter structures for MD. The quantification of FA, AD, and RD was overall in good agreement with the reference but the dependence of these parameters on WM anisotropy was somewhat biased (e.g. in corpus callosum). The inclination ϴ was well predicted for anisotropic WM structures, while the azimuth Φ was overall poorly predicted. The findings were highly consistent across both field strengths. Application of the optimized NN to high-resolution input data provided whole-brain maps with rich structural details. In conclusion, the proposed NN-driven approach showed potential to provide distortion-free high-resolution whole-brain maps of multiple diffusion metrics at high to ultrahigh field strengths in clinically relevant scan times.

A novel male Japanese quail structural connectivity atlas using ultra-high field diffusion MRI at 11.7 T

The structural connectivity of animal brains can be revealed using post-mortem diffusion-weighted magnetic resonance imaging (MRI). Despite the existence of several structural atlases of avian brains, few of them address the bird’s structural connectivity. In this study, a novel atlas of the structural connectivity is proposed for the male Japanese quail (Coturnix japonica), aiming at investigating two lines divergent on their emotionality trait: the short tonic immobility (STI) and the long tonic immobility (LTI) lines. The STI line presents a low emotionality trait, while the LTI line expresses a high emotionality trait. 21 male Japanese quail brains from both lines were scanned post-mortem for this study, using a preclinical Bruker 11.7 T MRI scanner. Diffusion-weighted MRI was performed using a 3D segmented echo planar imaging (EPI) pulsed gradient spin-echo (PGSE) sequence with a 200 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document}μm isotropic resolution, 75 diffusion-encoding directions and a b-value fixed at 4500 s/mm². Anatomical MRI was likewise performed using a 2D anatomical T₂-weighted spin-echo (SE) sequence with a 150 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document}μm isotropic resolution. This very first anatomical connectivity atlas of the male Japanese quail reveals 34 labeled fiber tracts and the existence of structural differences between the connectivity patterns characterizing the two lines. Thus, the link between the male Japanese quail’s connectivity and its underlying anatomical structures has reached a better understanding.

The Digital Brain Bank, an open access platform for post-mortem imaging datasets

Post-mortem magnetic resonance imaging (MRI) provides the opportunity to acquire high-resolution datasets to investigate neuroanatomy and validate the origins of image contrast through microscopy comparisons. We introduce the Digital Brain Bank (open.win.ox.ac.uk/DigitalBrainBank), a data release platform providing open access to curated, multimodal post-mortem neuroimaging datasets. Datasets span three themes-Digital Neuroanatomist: datasets for detailed neuroanatomical investigations; Digital Brain Zoo: datasets for comparative neuroanatomy; and Digital Pathologist: datasets for neuropathology investigations. The first Digital Brain Bank data release includes 21 distinctive whole-brain diffusion MRI datasets for structural connectivity investigations, alongside microscopy and complementary MRI modalities. This includes one of the highest-resolution whole-brain human diffusion MRI datasets ever acquired, whole-brain diffusion MRI in fourteen nonhuman primate species, and one of the largest post-mortem whole-brain cohort imaging studies in neurodegeneration. The Digital Brain Bank is the culmination of our lab's investment into post-mortem MRI methodology and MRI-microscopy analysis techniques. This manuscript provides a detailed overview of our work with post-mortem imaging to date, including the development of diffusion MRI methods to image large post-mortem samples, including whole, human brains. Taken together, the Digital Brain Bank provides cross-scale, cross-species datasets facilitating the incorporation of post-mortem data into neuroimaging studies.

Mapping the Human Connectome using Diffusion MRI at 300 mT/m Gradient Strength: Methodological Advances and Scientific Impact

Tremendous efforts have been made in the last decade to advance cutting-edge MRI technology in pursuit of mapping structural connectivity in the living human brain with unprecedented sensitivity and speed. The first Connectom 3T MRI scanner equipped with a 300 mT/m whole-body gradient system was installed at the Massachusetts General Hospital in 2011 and was specifically constructed as part of the Human Connectome Project. Since that time, numerous technological advances have been made to enable the broader use of the Connectom high gradient system for diffusion tractography and tissue microstructure studies and leverage its unique advantages and sensitivity to resolving macroscopic and microscopic structural information in neural tissue for clinical and neuroscientific studies. The goal of this review article is to summarize the technical developments that have emerged in the last decade to support and promote large-scale and scientific studies of the human brain using the Connectom scanner. We provide a brief historical perspective on the development of Connectom gradient technology and the efforts that led to the installation of three other Connectom 3T MRI scanners worldwide - one in the United Kingdom in Cardiff, Wales, another in Continental Europe in Leipzig, Germany, and the latest in Asia in Shanghai, China. We summarize the key developments in gradient hardware and image acquisition technology that have formed the backbone of Connectom-related research efforts, including the rich array of high-sensitivity receiver coils, pulse sequences, image artifact correction strategies and data preprocessing methods needed to optimize the quality of high-gradient strength dMRI data for subsequent analyses. Finally, we review the scientific impact of the Connectom MRI scanner, including advances in diffusion tractography, tissue microstructural imaging, ex vivo validation, and clinical investigations that have been enabled by Connectom technology. We conclude with brief insights into the unique value of strong gradients for dMRI and where the field is headed in the coming years.

A time-course study of actively stained mouse brains: Diffusion tensor imaging parameters and connectomic stability over 1 year

While the application of diffusion tensor imaging (DTI), tractography, and connectomics to fixed tissue is a common practice today, there have been limited studies examining the effects of fixation on brain microstructure over extended periods. This mouse model time-course study reports the changes of regional brain volumes and diffusion scalar parameters, such as fractional anisotropy, across 12 representative brain regions as measures of brain structural stability. The scalar DTI parameters and regional volumes were highly variable over the first 2 weeks after fixation. The same parameters were consistent over a 2-8-week window after fixation, which means confounds from tissue stability over that scanning window were minimal. Quantitative connectomes were analyzed over the same time with extension out to 1 year. While there was some change in the scalar metrics at 1 year after fixation, these changes were sufficiently small, particularly in white matter, to support reproducible connectomes over a period ranging from 2-weeks to 1-year post-fixation. These findings delineate a scanning period, during which brain volumes, diffusion scalar metrics, and connectomes are remarkably consistent.

A method to remove the influence of fixative concentration on postmortem T2 maps using a kinetic tensor model

Formalin fixation has been shown to substantially reduce T₂ estimates, primarily driven by the presence of fixative in tissue. Prior to scanning, post‐mortem samples are often placed into a fluid that has more favourable imaging properties. This study investigates whether there is evidence for a change in T₂ in regions close to the tissue surface due to fixative outflux into this surrounding fluid. Furthermore, we investigate whether a simulated spatial map of fixative concentration can be used as a confound regressor to reduce T₂ inhomogeneity. To achieve this, T₂ maps and diffusion tensor estimates were obtained in 14 whole, formalin‐fixed post‐mortem brains placed in Fluorinert approximately 48 hr prior to scanning. Seven brains were fixed with 10% formalin and seven brains were fixed with 10% neutral buffered formalin (NBF). Fixative outflux was modelled using a proposed kinetic tensor (KT) model, which incorporates voxelwise diffusion tensor estimates to account for diffusion anisotropy and tissue‐specific diffusion coefficients. Brains fixed with 10% NBF revealed a spatial T₂ pattern consistent with modelled fixative outflux. Confound regression of fixative concentration reduced T₂ inhomogeneity across both white and grey matter, with the greatest reduction attributed to the KT model versus simpler models of fixative outflux. No such effect was observed in brains fixed with 10% formalin. Correlations between the transverse relaxation rate R₂ and ferritin/myelin proteolipid protein (PLP) histology lead to an increased similarity for the relationship between R₂ and PLP for the two fixative types after KT correction.

Deep learning based multiplexed sensitivity-encoding (DL-MUSE) for high-resolution multi-shot DWI

PURPOSE

A phase correction method for high-resolution multi-shot (MSH) diffusion weighted imaging (DWI) is proposed. The efficacy and generalization capability of the method were validated on both healthy volunteers and patients.

THEORY AND METHODS

Conventionally, inter-shot phase variations for MSH echo-planar imaging (EPI) DWI are corrected by model-based algorithms. However, many acquisition imperfections are hard to measure accurately for conventional model-based methods, making the phase estimation and artifacts suppression unreliable. We propose a deep learning multiplexed sensitivity-encoding (DL-MUSE) framework to improve the phase estimations based on convolutional neural network (CNN) reconstruction. Aliasing-free single-shot (SSH) DW images, which have been used routinely in clinical settings, were used for training before the aliasing correction of MSH-DWI images. A dual-channel U-net comprising multiple convolutional layers was used for the phase estimation of MSH-EPI. The network was trained on a dataset containing 30 healthy volunteers and tested on another dataset of 52 healthy subjects and 15 patients with lesions or tumors with different shot numbers (4, 6 and 8). To further validate the generalization capability of our network, we acquired a dataset with different numbers of shots, TEs, partial Fourier factors, resolutions, ETLs, FOVs, coil numbers, and image orientations from two sites. We also compared the reconstruction performance of our proposed method with that of the conventional MUSE and SSH-EPI qualitatively and quantitatively.

RESULTS

Our results show that DL-MUSE is capable of correcting inter-shot phase errors with high and robust performance. Compared to conventional model-based MUSE, our method, by applying deep learning-based phase corrections, showed reduced distortion, noise level, and signal loss in high b-value DWIs. The improvements of image quality become more evident as the shot number increases from 4 to 8, especially in those central regions of the images, where g-factor artifacts are severe. Furthermore, the proposed method could provide the information about the orientation of the white matter with better consistency and achieve finer fibers delineation compared to the SSH-EPI method. Besides, the experiments on volunteers and patients from two different sites demonstrated the generalizability of our proposed method preliminarily.

CONCLUSION

A deep learning-based reconstruction algorithm for MSH-EPI images, which helps improve image quality greatly, was proposed. Results from healthy volunteers and tumor patients demonstrated the feasibility and generalization performances of our method for high-resolution MSH-EPI DWI, which can be used for routine clinical applications as well as neuroimaging research.

Connectome 2.0: Developing the next-generation ultra-high gradient strength human MRI scanner for bridging studies of the micro-, meso- and macro-connectome

The first phase of the Human Connectome Project pioneered advances in MRI technology for mapping the macroscopic structural connections of the living human brain through the engineering of a whole-body human MRI scanner equipped with maximum gradient strength of 300 mT/m, the highest ever achieved for human imaging. While this instrument has made important contributions to the understanding of macroscale connectional topology, it has also demonstrated the potential of dedicated high-gradient performance scanners to provide unparalleled in vivo assessment of neural tissue microstructure. Building on the initial groundwork laid by the original Connectome scanner, we have now embarked on an international, multi-site effort to build the next-generation human 3T Connectome scanner (Connectome 2.0) optimized for the study of neural tissue microstructure and connectional anatomy across multiple length scales. In order to maximize the resolution of this in vivo microscope for studies of the living human brain, we will push the diffusion resolution limit to unprecedented levels by (1) nearly doubling the current maximum gradient strength from 300 mT/m to 500 mT/m and tripling the maximum slew rate from 200 T/m/s to 600 T/m/s through the design of a one-of-a-kind head gradient coil optimized to minimize peripheral nerve stimulation; (2) developing high-sensitivity multi-channel radiofrequency receive coils for in vivo and ex vivo human brain imaging; (3) incorporating dynamic field monitoring to minimize image distortions and artifacts; (4) developing new pulse sequences to integrate the strongest diffusion encoding and highest spatial resolution ever achieved in the living human brain; and (5) calibrating the measurements obtained from this next-generation instrument through systematic validation of diffusion microstructural metrics in high-fidelity phantoms and ex vivo brain tissue at progressively finer scales with accompanying diffusion simulations in histology-based micro-geometries. We envision creating the ultimate diffusion MRI instrument capable of capturing the complex multi-scale organization of the living human brain - from the microscopic scale needed to probe cellular geometry, heterogeneity and plasticity, to the mesoscopic scale for quantifying the distinctions in cortical structure and connectivity that define cyto- and myeloarchitectonic boundaries, to improvements in estimates of macroscopic connectivity.

Allografts for Skin Closure during In Utero Spina Bifida Repair in a Sheep Model

Objectives: Use of off-label tissue graft materials, such as acellular dermal matrix (ADM), for in utero repair of severe spina bifida (SB), where primary skin layer closure is not possible, is associated with poor neurological outcomes. The cryopreserved human umbilical cord (HUC) patch has regenerative, anti-inflammatory, and anti-scarring properties, and provides watertight SB repair. We tested the hypothesis that the HUC is a superior skin patch to ADM for reducing inflammation at the repair site and preserving spinal cord function. Methods: In timed-pregnant ewes with twins, on gestational day (GD) 75, spina bifida was created without a myelotomy (functional model). On GD 95, repair was performed using HUC vs. ADM patches (randomly assigned) by suturing them to the skin edges. Additionally, full thickness skin closure as a primary skin closure (PSC) served as a positive control. Delivery was performed on GD 140, followed by blinded to treatment neurological assessments of the lambs using the Texas Spinal Cord Injury Scale (TSCIS) for gait, proprioception, and nociception. Lambs without spina bifida were used as controls (CTL). Ex vivo magnetic resonance imaging of spines at the repair site were performed, followed by quantitative pathological assessments. Histological assessments (blinded) included Masson’s trichrome, and immunofluorescence for myeloperoxidase (MPO; neutrophils) and for reactive astrocytes (inflammation) by co-staining vimentin and GFAP. Results: The combined hind limbs’ TSCIS was significantly higher in the HUC group than in ADM and PSC groups, p = 0.007. Both ADM and PSC groups exhibited loss of proprioception and mild to moderate ataxia compared to controls. MRI showed increased pathological findings in the PSC group when compared to the HUC group, p = 0.045. Histologically, the meningeal layer was thickened (inflammation) by 2–3 fold in ADM and PSC groups when compared to HUC and CTL groups, p = 0.01. There was lower MPO positive cells in the HUC group than in the ADM group, p = 0.018. Posterior column astrocyte activation was increased in ADM and PSC lambs compared to HUC lambs, p = 0.03. Conclusion: The HUC as a skin patch for in utero spina bifida repair preserves spinal cord function by reducing underlying inflammation when compared to ADM.

Slab boundary artifact correction in multislab imaging using convolutional-neural-network-enabled inversion for slab profile encoding

PURPOSE

This study aims to propose a novel algorithm for slab boundary artifact correction in both single-band multislab imaging and simultaneous multislab (SMSlab) imaging.

THEORY AND METHODS

In image domain, the formation of slab boundary artifacts can be regarded as modulating the artifact-free images using the slab profiles and introducing aliasing along the slice direction. Slab boundary artifact correction is the inverse problem of this process. An iterative algorithm based on convolutional neural networks (CNNs) is proposed to solve the problem, termed CNN-enabled inversion for slab profile encoding (CPEN). Diffusion-weighted SMSlab images and reference images without slab boundary artifacts were acquired in 7 healthy subjects for training. Images of 5 healthy subjects were acquired for testing, including single-band multislab and SMSlab images with 1.3-mm or 1-mm isotropic resolution. CNN-enabled inversion for slab profile encoding was compared with a previously reported method (i.e., nonlinear inversion for slab profile encoding [NPEN]).

RESULTS

CNN-enabled inversion for slab profile encoding reduces the slab boundary artifacts in both single-band multislab and SMSlab images. It also suppresses the slab boundary artifacts in the diffusion metric maps. Compared with NPEN, CPEN shows fewer residual artifacts in different acquisition protocols and more significant improvements in quantitative assessment, and it also accelerates the computation by more than 35 times.

CONCLUSION

CNN-enabled inversion for slab profile encoding can reduce the slab boundary artifacts in multislab acquisitions. It shows better slab boundary artifact correction capacity, higher robustness, and computation efficiency when compared with NPEN. It has the potential to improve the accuracy of multislab acquisitions in high-resolution DWI and functional MRI.

Diffusion-Weighted Imaging: Recent Advances and Applications

Quantitative diffusion imaging techniques enable the characterization of tissue microstructural properties of the human brain "in vivo", and are widely used in neuroscientific and clinical contexts. In this review, we present the basic physical principles behind diffusion imaging and provide an overview of the current diffusion techniques, including standard and advanced techniques as well as their main clinical applications. Standard diffusion tensor imaging (DTI) offers sensitivity to changes in microstructure due to diseases and enables the characterization of single fiber distributions within a voxel as well as diffusion anisotropy. Nonetheless, its inability to represent complex intravoxel fiber topologies and the limited biological specificity of its metrics motivated the development of several advanced diffusion MRI techniques. For example, high-angular resolution diffusion imaging (HARDI) techniques enabled the characterization of fiber crossing areas and other complex fiber topologies in a single voxel and supported the development of higher-order signal representations aiming to decompose the diffusion MRI signal into distinct microstructure compartments. Biophysical models, often known by their acronym (e.g., CHARMED, WMTI, NODDI, DBSI, DIAMOND) contributed to capture the diffusion properties from each of such tissue compartments, enabling the computation of voxel-wise maps of axonal density and/or morphology that hold promise as clinically viable biomarkers in several neurological and neuroscientific applications; for example, to quantify tissue alterations due to disease or healthy processes. Current challenges and limitations of state-of-the-art models are discussed, including validation efforts. Finally, novel diffusion encoding approaches (e.g., b-tensor or double diffusion encoding) may increase the biological specificity of diffusion metrics towards intra-voxel diffusion heterogeneity in clinical settings, holding promise in neurological applications.

Imaging multiple sclerosis pathology at 160 μm isotropic resolution by human whole-brain ex vivo magnetic resonance imaging at 3 T

Postmortem magnetic resonance imaging (MRI) of the fixed healthy and diseased human brain facilitates spatial resolutions and image quality that is not achievable with in vivo MRI scans. Though challenging-and almost exclusively performed at 7 T field strength-depicting the tissue architecture of the entire brain in fine detail is invaluable since it enables the study of neuroanatomy and uncovers important pathological features in neurological disorders. The objectives of the present work were (1) to develop a 3D isotropic ultra-high-resolution imaging approach for human whole-brain ex vivo acquisitions working on a standard clinical 3 T MRI system; and (2) to explore the sensitivity and specificity of this concept for specific pathoanatomical features of multiple sclerosis. The reconstructed images demonstrate unprecedented resolution and soft tissue contrast of the diseased human brain at 3 T, thus allowing visualization of sub-millimetric lesions in the different cortical layers and in the cerebellar cortex, as well as unique cortical lesion characteristics such as the presence of incomplete/complete iron rims, and patterns of iron accumulation. Further details such as the subpial molecular layer, the line of Gennari, and some intrathalamic nuclei are also well distinguishable.

Improving distortion correction for isotropic high-resolution 3D diffusion MRI by optimizing Jacobian modulation

PURPOSE

To improve distortion correction for isotropic high-resolution whole-brain 3D diffusion MRI when in a time-saving acquisition scenario.

THEORY AND METHODS

Data were acquired using simultaneous multi-slab (SMSlab) acquisitions, with a b = 0 image pair encoded by reversed polarity gradients (RPG) for phase encoding (PE) and diffusion weighted images encoded by a single PE direction. Eddy current-induced distortions were corrected first. During the following susceptibility distortion correction, image deformation was first corrected by the field map estimated from the b = 0 image pair. Intensity variation was subsequently corrected by Jacobian modulation. Two Jacobian modulation methods were compared. They calculated the Jacobian modulation map from the field map, or from the deformation corrected b = 0 image pair, termed as J_Field and J_RPG , respectively. A modified version of the J_RPG method, with proper smoothing, was further proposed for improved correction performance, termed as J_RPG-smooth .

RESULTS

Compared to J_Field modulation, less remaining distortions are observed when using the J_RPG and J_RPG-smooth methods, especially in areas with large B0 field inhomogeneity. The original J_RPG method causes signal-to-noise ratio (SNR) deficiency problem, which manifests as degraded SNR of the diffusion weighted images, while the J_RPG-smooth method maintains the original image SNR. Less estimation errors of diffusion metrics are observed when using the J_RPG-smooth method.

CONCLUSION

This study improves the distortion correction for isotropic high-resolution whole-brain 3D diffusion MRI by optimizing Jacobian modulation. The optimized method outperforms the conventional J_Field method regarding intensity variation correction and accuracy of diffusion metrics estimation, and outperforms the original J_RPG method regarding SNR performance.

Multi-modal imaging of a single mouse brain over five orders of magnitude of resolution

Mammalian neurons operate at length scales spanning six orders of magnitude; they project millimeters to centimeters across brain regions, are composed of micrometer-scale-diameter myelinated axons, and ultimately form nanometer scale synapses. Capturing these anatomical features across that breadth of scale has required imaging samples with multiple independent imaging modalities. Translating between the different modalities, however, requires imaging the same brain with each. Here, we imaged the same postmortem mouse brain over five orders of spatial resolution using MRI, whole brain micrometer-scale synchrotron x-ray tomography (μCT), and large volume automated serial electron microscopy. Using this pipeline, we can track individual myelinated axons previously relegated to axon bundles in diffusion tensor MRI or arbitrarily trace neurons and their processes brain-wide and identify individual synapses on them. This pipeline provides both an unprecedented look across a single brain’s multi-scaled organization as well as a vehicle for studying the brain’s multi-scale pathologies.

A 48-channel receive array coil for mesoscopic diffusion-weighted MRI of ex vivo human brain on the 3 T connectome scanner

In vivo diffusion-weighted magnetic resonance imaging is limited in signal-to-noise-ratio (SNR) and acquisition time, which constrains spatial resolution to the macroscale regime. Ex vivo imaging, which allows for arbitrarily long scan times, is critical for exploring human brain structure in the mesoscale regime without loss of SNR. Standard head array coils designed for patients are sub-optimal for imaging ex vivo whole brain specimens. The goal of this work was to design and construct a 48-channel ex vivo whole brain array coil for high-resolution and high b-value diffusion-weighted imaging on a 3T Connectome scanner. The coil was validated with bench measurements and characterized by imaging metrics on an agar brain phantom and an ex vivo human brain sample. The two-segment coil former was constructed for a close fit to a whole human brain, with small receive elements distributed over the entire brain. Imaging tests including SNR and G-factor maps were compared to a 64-channel head coil designed for in vivo use. There was a 2.9-fold increase in SNR in the peripheral cortex and a 1.3-fold gain in the center when compared to the 64-channel head coil. The 48-channel ex vivo whole brain coil also decreases noise amplification in highly parallel imaging, allowing acceleration factors of approximately one unit higher for a given noise amplification level. The acquired diffusion-weighted images in a whole ex vivo brain specimen demonstrate the applicability and advantage of the developed coil for high-resolution and high b-value diffusion-weighted ex vivo brain MRI studies.

High-resolution whole-brain diffusion MRI at 3T using simultaneous multi-slab (SMSlab) acquisition

High-resolution diffusion MRI (dMRI) is a crucial tool in neuroscience studies to detect fine fiber structure, depict complex fiber architecture and analyze cortical anisotropy. However, high-resolution dMRI is limited by its intrinsically low SNR due to diffusion attenuation. A series of techniques have been proposed to improve the SNR performance, but most of them are at the cost of long scan time, which in turn sacrifice the SNR efficiency, especially for large FOV imaging, such as whole-brain imaging. Recently, a combination of 3D multi-slab acquisition and simultaneous multi-slice (SMS) excitation, namely simultaneous multi-slab (SMSlab), has been demonstrated to have potential for high-resolution diffusion imaging with high SNR and SNR efficiency. In our previous work, we have proposed a 3D Fourier encoding and reconstruction framework for SMSlab acquisition. In this study, we extend this 3D k-space framework to diffusion imaging, by developing a novel navigator acquisition strategy and exploring a k-space-based phase correction method. In vivo brain data are acquired using the proposed SMSlab method and compared with a series of different acquisitions, including the traditional 3D multi-slab, 2D SMS and 2D single-shot EPI (ss-EPI) acquisitions. The results demonstrate that SMSlab has a better SNR performance compared with 3D multi-slab and 2D SMS. The detection capacity of fine fiber structures is improved using SMSlab, compared with the low-resolution diffusion imaging using conventional 2D ss-EPI.

In vivo human whole-brain Connectom diffusion MRI dataset at 760 µm isotropic resolution

We present a whole-brain in vivo diffusion MRI (dMRI) dataset acquired at 760 μm isotropic resolution and sampled at 1260 q-space points across 9 two-hour sessions on a single healthy participant. The creation of this benchmark dataset is possible through the synergistic use of advanced acquisition hardware and software including the high-gradient-strength Connectom scanner, a custom-built 64-channel phased-array coil, a personalized motion-robust head stabilizer, a recently developed SNR-efficient dMRI acquisition method, and parallel imaging reconstruction with advanced ghost reduction algorithm. With its unprecedented resolution, SNR and image quality, we envision that this dataset will have a broad range of investigational, educational, and clinical applications that will advance the understanding of human brain structures and connectivity. This comprehensive dataset can also be used as a test bed for new modeling, sub-sampling strategies, denoising and processing algorithms, potentially providing a common testing platform for further development of in vivo high resolution dMRI techniques. Whole brain anatomical T1-weighted and T2-weighted images at submillimeter scale along with field maps are also made available.

Automated cerebral cortex segmentation based solely on diffusion tensor imaging for investigating cortical anisotropy

To extract Diffusion Tensor Imaging (DTI) parameters from the human cortex, the inner and outer boundaries of the cortex are usually defined on 3D-T1-weighted images and then applied to the co-registered DTI. However, this analysis requires the acquisition of an additional high-resolution structural image that may not be practical in various imaging studies. Here an automatic cortical boundary segmentation method was developed to work directly only on the native DTI images by using fractional anisotropy (FA) maps and mean diffusion weighted images (DWI), the latter with acceptable gray-white matter image contrast. This new method was compared to the conventional cortical segmentations generated from high-resolution T1 structural images in 5 participants. In addition, the proposed method was applied to 15 healthy young adults (10 cross-sectional, 5 test-retest) to measure FA, MD, and radiality of the primary eigenvector across the cortex on whole-brain 1.5 mm isotropic images acquired in 3.5 min at 3T. The proposed method generated reasonable segmentations of the cortical boundaries for all individuals and large proportions of the proposed method segmentations (more than 85%) were within ±1 mm from those generated with the conventional approach on higher resolution T1 structural images. Both FA (~0.15) and MD (~0.77 × 10^-3 mm²/s) extracted halfway between the cortical boundaries were relatively stable across the cortex, although focal regions such as the posterior bank of the central sulcus, anterior insula, and medial temporal lobe showed higher FA. The primary eigenvectors were primarily oriented radially to the middle cortical surface, but there were tangential orientations in the sulcal fundi as well as in the posterior bank of the central sulcus. The proposed method demonstrates the feasibility and accuracy of cortical analysis in native DTI space while avoiding the acquisition of other imaging contrasts like 3D T1-weighted scans.

WIKIBrainStem: An online atlas to manually segment the human brainstem at the mesoscopic scale from ultrahigh field MRI

The brainstem is one of the most densely packed areas of the central nervous system in terms of gray, but also white, matter structures and, therefore, is a highly functional hub. It has mainly been studied by the means of histological techniques, which requires several hundreds of slices with a loss of the 3D coherence of the whole specimen. Access to the inner structure of the brainstem is possible using Magnetic Resonance Imaging (MRI), but this method has a limited spatial resolution and contrast in vivo. Here, we scanned an ex vivo specimen using an ultra-high field (11.7T) preclinical MRI scanner providing data at a mesoscopic scale for anatomical T2-weighted (100 µm and 185 µm isotropic) and diffusion-weighted imaging (300 µm isotropic). We then proposed a hierarchical segmentation of the inner gray matter of the brainstem and defined a set of rules for each segmented anatomical class. These rules were gathered in a freely accessible web-based application, WIKIBrainStem (https://fibratlas.univ-tours.fr/brainstems/index.html), for 99 structures, from which 13 were subdivided into 29 substructures. This segmentation is, to date, the most detailed one developed from ex vivo MRI of the brainstem. This should be regarded as a tool that will be complemented by future results of alternative methods, such as Optical Coherence Tomography, Polarized Light Imaging or histology… This is a mandatory step prior to segmenting multiple specimens, which will be used to create a probabilistic automated segmentation method of ex vivo, but also in vivo, brainstem and may be used for targeting anatomical structures of interest in managing some degenerative or psychiatric disorders.

Extreme capsule is a bottleneck for ventral pathway

As neuroscience literature suggests, extreme capsule is considered a whiter matter tract. Nevertheless, it is not clear whether extreme capsule itself is an association fiber pathway or only a bottleneck for other association fibers to pass. Via our review, investigating anatomical position, connectivity and cognitive role of the bundles in extreme capsule, and by analyzing data from the dissection, it can be argued that extreme capsule is probably a bottleneck for the passage of uncinated fasciculus (UF) and inferior fronto-occipital fasciculus (IFOF), and these fasciculi are responsible for the respective roles in language processing.

Dedicated container for postmortem human brain ultra-high field magnetic resonance imaging

BACKGROUND

The emerging field of ultra-high field MRI (UHF-MRI, 7 Tesla and higher) provides the opportunity to image human brains at a higher resolution and with higher signal-to-noise ratios compared to the more widely available 1.5 and 3T scanners. Scanning postmortem tissue additionally allows for greatly increased scan times and fewer movement issues leading to improvements in image quality. However, typical postmortem neuroimaging routines involve placing the tissue within plastic bags that leave room for susceptibility artifacts from tissue-air interfaces, inadequate submersion, and leakage issues. To address these challenges in postmortem imaging, a custom-built nonferromagnetic container was developed that allows whole brain hemispheres to be scanned at sub-millimeter resolution within typical head-coils.

METHOD

The custom-built polymethylmethacrylaat container consists of a cylinder with a hemispheric side and a lid with valves on the adjacent side. This shape fits within common MR head-coils and allows whole hemispheres to be submerged and vacuum sealed within it reducing imaging artifacts that would otherwise arise at air-tissue boundaries. Two hemisphere samples were scanned on a Siemens 9.4T Magnetom MRI scanner. High resolution T2* weighted data was obtained with a custom 3D gradient echo (GRE) sequence and diffusion-weighted imaging (DWI) scans were obtained with a 3D kT-dSTEAM sequence along 48 directions.

RESULTS

The custom-built container proved to submerge and contain tissue samples effectively and showed no interferences with MR scanning acquisition. The 3D GRE sequence provided high resolution isotropic T2* weighted data at 250 μm which showed a clear visualization of gray and white matter structures. DWI scans allowed for dense reconstruction of structural white matter connections via tractography.

CONCLUSION

Using this custom-built container worked towards achieving high quality MR images of postmortem brain material. This procedure can have advantages over traditional schemes including utilization of a standardized protocol and the reduced likelihood of leakage. This methodology could be adjusted and used to improve typical postmortem imaging routines.