Quantitative susceptibility mapping of the rat brain after traumatic brain injury

ComBating inter-site differences in field strength: harmonizing preclinical traumatic brain injury MRI data

Integrating datasets from multiple sites and scanners can increase statistical power for neuroimaging studies but can also introduce significant inter-site confounds. We evaluated the effectiveness of ComBat, an empirical Bayes approach, to combine longitudinal preclinical MRI data acquired at 4.7 or 9.4 T at two different sites in Australia. Male Sprague Dawley rats underwent MRI on Days 2, 9, 28, and 150 following moderate/severe traumatic brain injury (TBI) or sham injury as part of Project 1 of the NIH/NINDS-funded Centre Without Walls EpiBioS4Rx project. Diffusion-weighted and multiple-gradient-echo images were acquired, and outcomes included QSM, FA, and ADC. Acute injury measures including apnea and self-righting reflex were consistent between sites. Mixed-effect analysis of ipsilateral and contralateral corpus callosum (CC) summary values revealed a significant effect of site on FA and ADC values, which was removed following ComBat harmonization. Bland-Altman plots for each metric showed reduced variability across sites following ComBat harmonization, including for QSM, despite appearing to be largely unaffected by inter-site differences and no effect of site observed. Following harmonization, the combined inter-site data revealed significant differences in the imaging metrics consistent with previously reported outcomes. TBI resulted in significantly reduced FA and increased susceptibility in the ipsilateral CC, and significantly reduced FA in the contralateral CC compared with sham-injured rats. Additionally, TBI rats also exhibited a reversal in ipsilateral CC ADC values over time with significantly reduced ADC at Day 9, followed by increased ADC 150 days after injury. Our findings demonstrate the need for harmonizing multi-site preclinical MRI data and show that this can be successfully achieved using ComBat while preserving phenotypical changes due to TBI.

Myo-inositol Levels in the Dorsal Anterior Cingulate Cortex Predicts Anxiety-to-Eat in Anorexia Nervosa

Background

Anorexia nervosa (AN) is a mental and behavioral health condition characterized by an intense fear of weight or fat gain, severe restriction of food intake resulting in low body weight, and distorted self-perception of body shape or weight. While substantial research has focused on general anxiety in AN, less is known about eating-related anxiety and its underlying neural mechanisms. Therefore, we sought to characterize anxiety-to-eat in AN and examine the neurometabolic profile within the dorsal anterior cingulate cortex (dACC), a brain region putatively involved in magnifying the threat response.

Methods

Women seeking inpatient treatment for AN and women of healthy weight without a lifetime history of an eating disorder (healthy controls; HC) completed a computer-based behavioral task assessing anxiety-to-eat in response to images of higher (HED) and lower (LED) energy density foods. Participants also underwent magnetic resonance spectroscopy of the dACC in a 3 Tesla scanner.

Results

The AN group reported greater anxiety to eat HED and LED foods relative to the HC group. Both groups reported greater anxiety to eat HED foods relative to LED foods. The neurometabolite myo-inositol (mI) was lower in the dACC in AN relative to HC, and mI levels negatively predicted anxiety to eat HED but not LED foods in the AN group only. mI levels in the dACC were independent of body weight, body mass, and general anxiety.

Conclusions

These findings provide critical new insight into the clinically challenging feature and underlying neural mechanisms of eating-related anxiety and indicate mI levels in the dACC could serve as a novel biomarker of illness severity that is independent of body weight to identify individuals vulnerable to disordered eating or eating pathology as well as a potential therapeutic target.

Improving quantitative susceptibility mapping for the identification of traumatic brain injury neurodegeneration at the individual level

BACKGROUND

Emerging evidence suggests that traumatic brain injury (TBI) is a major risk factor for developing neurodegenerative disease later in life. Quantitative susceptibility mapping (QSM) has been used by an increasing number of studies in investigations of pathophysiological changes in TBI. However, generating artefact-free quantitative susceptibility maps in brains with large focal lesions, as in the case of moderate-to-severe TBI (ms-TBI), is particularly challenging. To address this issue, we utilized a novel two-pass masking technique and reconstruction procedure (two-pass QSM) to generate quantitative susceptibility maps (QSMxT; Stewart et al., 2022, Magn Reson Med.) in combination with the recently developed virtual brain grafting (VBG) procedure for brain repair (Radwan et al., 2021, NeuroImage) to improve automated delineation of brain areas. We used QSMxT and VBG to generate personalised QSM profiles of individual patients with reference to a sample of healthy controls.

METHODS

Chronic ms-TBI patients (N = 8) and healthy controls (N = 12) underwent (multi-echo) GRE, and anatomical MRI (MPRAGE) on a 3T Siemens PRISMA scanner. We reconstructed the magnetic susceptibility maps using two-pass QSM from QSMxT. We then extracted values of magnetic susceptibility in grey matter (GM) regions (following brain repair via VBG) across the whole brain and determined if they deviate from a reference healthy control group [Z-score < -3.43 or > 3.43, relative to the control mean], with the aim of obtaining personalised QSM profiles.

RESULTS

Using two-pass QSM, we achieved susceptibility maps with a substantial increase in quality and reduction in artefacts irrespective of the presence of large focal lesions, compared to single-pass QSM. In addition, VBG minimised the loss of GM regions and exclusion of patients due to failures in the region delineation step. Our findings revealed deviations in magnetic susceptibility measures from the HC group that differed across individual TBI patients. These changes included both increases and decreases in magnetic susceptibility values in multiple GM regions across the brain.

CONCLUSIONS

We illustrate how to obtain magnetic susceptibility values at the individual level and to build personalised QSM profiles in ms-TBI patients. Our approach opens the door for QSM investigations in more severely injured patients. Such profiles are also critical to overcome the inherent heterogeneity of clinical populations, such as ms-TBI, and to characterize the underlying mechanisms of neurodegeneration at the individual level more precisely. Moreover, this new personalised QSM profiling could in the future assist clinicians in assessing recovery and formulating a neuroscience-guided integrative rehabilitation program tailored to individual TBI patients.

Spatiospectral image processing workflow considerations for advanced MR spectroscopy of the brain

Magnetic resonance spectroscopy (MRS) is one of the few non-invasive imaging modalities capable of making neurochemical and metabolic measurements in vivo. Traditionally, the clinical utility of MRS has been narrow. The most common use has been the "single-voxel spectroscopy" variant to discern the presence of a lactate peak in the spectra in one location in the brain, typically to evaluate for ischemia in neonates. Thus, the reduction of rich spectral data to a binary variable has not classically necessitated much signal processing. However, scanners have become more powerful and MRS sequences more advanced, increasing data complexity and adding 2 to 3 spatial dimensions in addition to the spectral one. The result is a spatially- and spectrally-variant MRS image ripe for image processing innovation. Despite this potential, the logistics for robustly accessing and manipulating MRS data across different scanners, data formats, and software standards remain unclear. Thus, as research into MRS advances, there is a clear need to better characterize its image processing considerations to facilitate innovation from scientists and engineers. Building on established neuroimaging standards, we describe a framework for manipulating these images that generalizes to the voxel, spectral, and metabolite level across space and multiple imaging sites while integrating with LCModel, a widely used quantitative MRS peak-fitting platform. In doing so, we provide examples to demonstrate the advantages of such a workflow in relation to recent publications and with new data. Overall, we hope our characterizations will lower the barrier of entry to MRS processing for neuroimaging researchers.

Meta-analysis and Open-source Database for In Vivo Brain Magnetic Resonance Spectroscopy in Health and Disease

Proton ( 1 H) Magnetic Resonance Spectroscopy (MRS) is a non-invasive tool capable of quantifying brain metabolite concentrations in vivo . Prioritization of standardization and accessibility in the field has led to the development of universal pulse sequences, methodological consensus recommendations, and the development of open-source analysis software packages. One on-going challenge is methodological validation with ground-truth data. As ground-truths are rarely available for in vivo measurements, data simulations have become an important tool. The diverse literature of metabolite measurements has made it challenging to define ranges to be used within simulations. Especially for the development of deep learning and machine learning algorithms, simulations must be able to produce accurate spectra capturing all the nuances of in vivo data. Therefore, we sought to determine the physiological ranges and relaxation rates of brain metabolites which can be used both in data simulations and as reference estimates. Using the Preferred Reporting Items for Systematic reviews and Meta-Analyses (PRISMA) guidelines, we've identified relevant MRS research articles and created an open-source database containing methods, results, and other article information as a resource. Using this database, expectation values and ranges for metabolite concentrations and T 2 relaxation times are established based upon a meta-analyses of healthy and diseased brains.

Axon fiber orientation as the source of T1 relaxation anisotropy in white matter: A study on corpus callosum in vivo and ex vivo

PURPOSE

Recent studies indicate that T₁ in white matter (WM) is influenced by fiber orientation in B₀ . The purpose of the study was to investigate the interrelationships between axon fiber orientation in corpus callosum (CC) and T₁ relaxation time in humans in vivo as well as in rat brain ex vivo.

METHODS

Volunteers were scanned for relaxometric and diffusion MRI at 3 T and 7 T. Angular T₁ plots from WM were computed using fractional anisotropy and fiber-to-field-angle maps. T₁ and fiber-to-field angle were measured in five sections of CC to estimate the effects of inherently varying fiber orientations on T₁ within the same tracts in vivo. Ex vivo rat-brain preparation encompassing posterior CC was rotated in B₀ and T₁ , and diffusion MRI images acquired at 9.4 T. T₁ angular plots were determined at several rotation angles in B₀ .

RESULTS

Angular T₁ plots from global WM provided reference for estimated fiber orientation-linked T₁ changes within CC. In anterior midbody of CC in vivo, where small axons are dominantly present, a shift in axon orientation is accompanied by a change in T₁ , matching that estimated from WM T₁ data. In CC, where large and giant axons are numerous, the measured T₁ change is about 2-fold greater than the estimated one. Ex vivo rotation of the same midsagittal CC region of interest produced angular T₁ plots at 9.4 T, matching those observed at 7 T in vivo.

CONCLUSION

These data causally link axon fiber orientation in B₀ to the T₁ relaxation anisotropy in WM.

Investigating the effect of oblique image acquisition on the accuracy of QSM and a robust tilt correction method

PURPOSE

Quantitative susceptibility mapping (QSM) is used increasingly for clinical research where oblique image acquisition is commonplace, but its effects on QSM accuracy are not well understood.

THEORY AND METHODS

The QSM processing pipeline involves defining the unit magnetic dipole kernel, which requires knowledge of the direction of the main magnetic fieldB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ with respect to the acquired image volume axes. The direction ofB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ is dependent on the axis and angle of rotation in oblique acquisition. Using both a numerical brain phantom and in vivo acquisitions in 5 healthy volunteers, we analyzed the effects of oblique acquisition on magnetic susceptibility maps. We compared three tilt-correction schemes at each step in the QSM pipeline: phase unwrapping, background field removal and susceptibility calculation, using the RMS error and QSM-tuned structural similarity index.

RESULTS

Rotation of wrapped phase images gave severe artifacts. Background field removal with projection onto dipole fields gave the most accurate susceptibilities when the field map was first rotated into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ . Laplacian boundary value and variable-kernel sophisticated harmonic artifact reduction for phase data background field removal methods gave accurate results without tilt correction. For susceptibility calculation, thresholded k-space division, iterative Tikhonov regularization, and weighted linear total variation regularization, all performed most accurately when local field maps were rotated into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ before susceptibility calculation.

CONCLUSION

For accurate QSM, oblique acquisition must be taken into account. Rotation of images into alignment withB ^ 0 $$ {\hat{\boldsymbol{B}}}_{\mathbf{0}} $$ should be carried out after phase unwrapping and before background-field removal. We provide open-source tilt-correction code to incorporate easily into existing pipelines: https://github.com/o-snow/QSM_TiltCorrection.git.

The Significance of Hypothalamic Inflammation and Gliosis for the Pathogenesis of Obesity in Humans

Accumulated preclinical literature demonstrates that hypothalamic inflammation and gliosis are underlying causal components of diet-induced obesity in rodent models. This review summarizes and synthesizes available translational data to better understand the applicability of preclinical findings to human obesity and its comorbidities. The published literature in humans includes histopathologic analyses performed postmortem and in vivo neuroimaging studies measuring indirect markers of hypothalamic tissue microstructure. Both support the presence of hypothalamic inflammation and gliosis in children and adults with obesity. Findings predominantly point to tissue changes in the region of the arcuate nucleus of the hypothalamus, although findings of altered tissue characteristics in whole hypothalamus or other hypothalamic regions also emerged. Moreover, the severity of hypothalamic inflammation and gliosis has been related to comorbid conditions, including glucose intolerance, insulin resistance, type 2 diabetes, and low testosterone levels in men, independent of elevated body adiposity. Cross-sectional findings are augmented by a small number of prospective studies suggesting that a greater degree of hypothalamic inflammation and gliosis may predict adiposity gain and worsening insulin sensitivity in susceptible individuals. In conclusion, existing human studies corroborate a large preclinical literature demonstrating that hypothalamic neuroinflammatory responses play a role in obesity pathogenesis. Extensive or permanent hypothalamic tissue remodeling may negatively affect the function of neuroendocrine regulatory circuits and promote the development and maintenance of elevated body weight in obesity and/or comorbid endocrine disorders.

Diffusion MRI approaches for investigating microstructural complexity in a rat model of traumatic brain injury

Our study explores the potential of conventional and advanced diffusion MRI techniques including diffusion tensor imaging (DTI), and single-shell 3-tissue constrained spherical deconvolution (SS3T-CSD) to investigate complex microstructural changes following severe traumatic brain injury in rats at a chronic phase. Rat brains after sham-operation or lateral fluid percussion (LFP) injury were scanned ex vivo in a 9.4 T scanner. Our region-of-interest-based approach of tensor-, and SS3T-CSD derived fixel-, 3-tissue signal fraction maps were sensitive to changes in both white matter (WM) and grey matter (GM) areas. Tensor-based measures, such as fractional anisotropy (FA) and radial diffusivity (RD), detected more changes in WM and GM areas as compared to fixel-based measures including apparent fiber density (AFD), peak FOD amplitude and primary fiber bundle density, while 3-tissue signal fraction maps revealed distinct changes in WM, GM, and phosphate-buffered saline (PBS) fractions highlighting the complex tissue microstructural alterations post-trauma. Track-weighted imaging demonstrated changes in track morphology including reduced curvature and average pathlength distal from the primary lesion in severe TBI rats. In histological analysis, changes in the diffusion MRI measures could be associated to decreased myelin density, loss of myelinated axons, and increased cellularity, revealing progressive microstructural alterations in these brain areas five months after injury. Overall, this study highlights the use of combined conventional and advanced diffusion MRI measures to obtain more precise insights into the complex tissue microstructural alterations in chronic phase of severe brain injury.

A Novel Rodent Model of Hypertensive Cerebral Small Vessel Disease with White Matter Hyperintensities and Peripheral Oxidative Stress

Cerebral small vessel disease (CSVD) is the second most common cause of stroke and a major contributor to dementia. Manifestations of CSVD include cerebral microbleeds, intracerebral hemorrhages (ICH), lacunar infarcts, white matter hyperintensities (WMH) and enlarged perivascular spaces. Chronic hypertensive models have been found to reproduce most key features of the disease. Nevertheless, no animal models have been identified to reflect all different aspects of the human disease. Here, we described a novel model for CSVD using salt-sensitive ‘Sabra’ hypertension-prone rats (SBH/y), which display chronic hypertension and enhanced peripheral oxidative stress. SBH/y rats were either administered deoxycorticosteroid acetate (DOCA) (referred to as SBH/y-DOCA rats) or sham-operated and provided with 1% NaCl in drinking water. Rats underwent neurological assessment and behavioral testing, followed by ex vivo MRI and biochemical and histological analyses. SBH/y-DOCA rats show a neurological decline and cognitive impairment and present multiple cerebrovascular pathologies associated with CSVD, such as ICH, lacunes, enlarged perivascular spaces, blood vessel stenosis, BBB permeability and inflammation. Remarkably, SBH/y-DOCA rats show severe white matter pathology as well as WMH, which are rarely reported in commonly used models. Our model may serve as a novel platform for further understanding the mechanisms underlying CSVD and for testing novel therapeutics.

Detection of glioblastoma multiforme using quantitative molecular magnetic resonance imaging based on 5-aminolevulinic acid: in vitro and in vivo studies

OBJECTIVES

We demonstrated a novel metabolic method based on sequential administration of 5-aminolevulinic acid (ALA) and iron supplement, and ferric ammonium citrate (FAC), for glioblastoma multiforme (GBM) detection using R₂' and quantitative susceptibility mapping (QSM).

MATERIALS AND METHODS

Intra-cellular iron accumulation in glioblastoma cells treated with ALA and/or FAC was measured. Cell phantoms containing glioblastoma cells and Wistar rats bearing C₆ glioblastoma were imaged using a 3 T MRI scanner after sequential administration of ALA and FAC. The relaxivity and QSM analysis were performed on the images.

RESULTS

The intra-cellular iron deposition was significantly higher in the glioma cells with sequential treatment of ALA and FAC for 6 h compared to those treated with the controls. The relaxivity and magnetic susceptibility values of the glioblastoma cells and rat brain tumors treated with ALA + FAC (115 ± 5 s-1 for R₂', and 0.1 ± 0.02 ppm for magnetic susceptibility) were significantly higher than those treated with the controls (55 ± 18 (FAC), 45 ± 15 (ALA) s-1 for R₂', p < 0.05, and 0.03 ± 0.03 (FAC), 0.02 ± 0.02 (ALA) ppm for magnetic susceptibility, p < 0.05).

DISCUSSION

Sequential administration of ALA and iron supplements increases the iron deposition in glioblastoma cells, enabling clinical 3 T MRI to detect GBM using R₂' or QSM.

Microstructural Tissue Changes in a Rat Model of Mild Traumatic Brain Injury

Our study investigates the potential of diffusion MRI (dMRI), including diffusion tensor imaging (DTI), fixel-based analysis (FBA) and neurite orientation dispersion and density imaging (NODDI), to detect microstructural tissue abnormalities in rats after mild traumatic brain injury (mTBI). The brains of sham-operated and mTBI rats 35 days after lateral fluid percussion injury were imaged ex vivo in a 11.7-T scanner. Voxel-based analyses of DTI-, fixel- and NODDI-based metrics detected extensive tissue changes in directly affected brain areas close to the primary injury, and more importantly, also in distal areas connected to primary injury and indirectly affected by the secondary injury mechanisms. Histology revealed ongoing axonal abnormalities and inflammation, 35 days after the injury, in the brain areas highlighted in the group analyses. Fractional anisotropy (FA), fiber density (FD) and fiber density and fiber bundle cross-section (FDC) showed similar pattern of significant areas throughout the brain; however, FA showed more significant voxels in gray matter areas, while FD and FDC in white matter areas, and orientation dispersion index (ODI) in areas most damage based on histology. Region-of-interest (ROI)-based analyses on dMRI maps and histology in selected brain regions revealed that the changes in MRI parameters could be attributed to both alterations in myelinated fiber bundles and increased cellularity. This study demonstrates that the combination of dMRI methods can provide a more complete insight into the microstructural alterations in white and gray matter after mTBI, which may aid diagnosis and prognosis following a mild brain injury.

Quantitative MRI susceptibility mapping reveals cortical signatures of changes in iron, calcium and zinc in malformations of cortical development in children with drug-resistant epilepsy

OBJECTIVE

Malformations of cortical development (MCD), including focal cortical dysplasia (FCD), are the most common cause of drug-resistant focal epilepsy in children. Histopathological lesion characterisation demonstrates abnormal cell types and lamination, alterations in myelin (typically co-localised with iron), and sometimes calcification. Quantitative susceptibility mapping (QSM) is an emerging MRI technique that measures tissue magnetic susceptibility (χ) reflecting it's mineral composition. We used QSM to investigate abnormal tissue composition in a group of children with focal epilepsy with comparison to effective transverse relaxation rate (R2*) and Synchrotron radiation X-ray fluorescence (SRXRF) elemental maps. Our primary hypothesis was that reductions in χ would be found in FCD lesions, resulting from alterations in their iron and calcium content. We also evaluated deep grey matter nuclei for changes in χ with age.

METHODS

QSM and R2* maps were calculated for 40 paediatric patients with suspected MCD (18 histologically confirmed) and 17 age-matched controls. Patients' sub-groups were defined based on concordant electro-clinical or histopathology data. Quantitative investigation of QSM and R2* was performed within lesions, using a surface-based approach with comparison to homologous regions, and within deep brain regions using a voxel-based approach with regional values modelled with age and epilepsy as covariates. Synchrotron radiation X-ray fluorescence (SRXRF) was performed on brain tissue resected from 4 patients to map changes in iron, calcium and zinc and relate them to MRI parameters.

RESULTS

Compared to fluid-attenuated inversion recovery (FLAIR) or T1-weighted imaging, QSM improved lesion conspicuity in 5% of patients. In patients with well-localised lesions, quantitative profiling demonstrated decreased χ, but not R2*, across cortical depth with respect to the homologous regions. Contra-lateral homologous regions additionally exhibited increased χ at 2-3 mm cortical depth that was absent in lesions. The iron decrease measured by the SRXRF in FCDIIb lesions was in agreement with myelin reduction observed by Luxol Fast Blue histochemical staining. SRXRF analysis in two FCDIIb tissue samples showed increased zinc and calcium in one patient, and decreased iron in the brain region exhibiting low χ and high R2* in both patients. QSM revealed expected age-related changes in the striatum nuclei, substantia nigra, sub-thalamic and red nucleus.

CONCLUSION

QSM non-invasively revealed cortical/sub-cortical tissue alterations in MCD lesions and in particular that χ changes in FCDIIb lesions were consistent with reduced iron, co-localised with low myelin and increased calcium and zinc content. These findings suggest that measurements of cortical χ could be used to characterise tissue properties non-invasively in epilepsy lesions.

Emerging Applications for Quantitative Susceptibility Mapping in the Detection of Traumatic Brain Injury Pathology

Traumatic brain injury (TBI) is a common but heterogeneous injury underpinned by numerous complex and interrelated pathophysiological mechanisms. An essential trace element, iron is abundant within the brain and involved in many fundamental neurobiological processes, including oxygen transportation, oxidative phosphorylation, myelin production and maintenance, as well as neurotransmitter synthesis and metabolism. Excessive levels of iron are neurotoxic and thus iron homeostasis is tightly regulated in the brain, however, many details about the mechanisms by which this is achieved are yet to be elucidated. A key mediator of oxidative stress, mitochondrial dysfunction and neuroinflammatory response, iron dysregulation is an important contributor to secondary injury in TBI. Advances in neuroimaging that leverage magnetic susceptibility properties have enabled increasingly comprehensive investigations into the distribution and behaviour of iron in the brain amongst healthy individuals as well as disease states such as TBI. Quantitative Susceptibility Mapping (QSM) is an advanced neuroimaging technique that promises quantitative estimation of local magnetic susceptibility at the voxel level. In this review, we provide an overview of brain iron and its homeostasis, describe recent advances enabling applications of QSM within the context of TBI and summarise the current state of the literature. Although limited, the emergent research suggests that QSM is a promising neuroimaging technique that can be used to investigate a host of pathophysiological changes that are associated with TBI.