Effect of respiratory hyperoxic challenge on magnetic susceptibility in human brain assessed by quantitative susceptibility mapping (QSM)

BFRnet: A deep learning-based MR background field removal method for QSM of the brain containing significant pathological susceptibility sources

INTRODUCTION

Background field removal (BFR) is a critical step required for successful quantitative susceptibility mapping (QSM). However, eliminating the background field in brains containing significant susceptibility sources, such as intracranial hemorrhages, is challenging due to the relatively large scale of the field induced by these pathological susceptibility sources.

METHOD

This study proposes a new deep learning-based method, BFRnet, to remove the background field in healthy and hemorrhagic subjects. The network is built with the dual-frequency octave convolutions on the U-net architecture, trained with synthetic field maps containing significant susceptibility sources. The BFRnet method is compared with three conventional BFR methods and one previous deep learning method using simulated and in vivo brains from 4 healthy and 2 hemorrhagic subjects. Robustness against acquisition field-of-view (FOV) orientation and brain masking are also investigated.

RESULTS

For both simulation and in vivo experiments, BFRnet led to the best visually appealing results in the local field and QSM results with the minimum contrast loss and the most accurate hemorrhage susceptibility measurements among all five methods. In addition, BFRnet produced the most consistent local field and susceptibility maps between different sizes of brain masks, while conventional methods depend drastically on precise brain extraction and further brain edge erosions. It is also observed that BFRnet performed the best among all BFR methods for acquisition FOVs oblique to the main magnetic field.

CONCLUSION

The proposed BFRnet improved the accuracy of local field reconstruction in the hemorrhagic subjects compared with conventional BFR algorithms. The BFRnet method was effective for acquisitions of tilted orientations and retained whole brains without edge erosion as often required by traditional BFR methods.

Optimization of quantitative susceptibility mapping for regional estimation of oxygen extraction fraction in the brain

PURPOSE

We sought to determine the degree to which oxygen extraction fraction (OEF) estimated using quantitative susceptibility mapping (QSM) depends on two critical acquisition parameters that have a significant impact on acquisition time: voxel size and final echo time.

METHODS

Four healthy volunteers were imaged using a range of isotropic voxel sizes and final echo times. The 0.7 mm data were downsampled at different stages of QSM processing by a factor of 2 (to 1.4 mm), 3 (2.1 mm), or 4 (2.8 mm) to determine the impact of voxel size on each analysis step. OEF was estimated from 11 veins of varying diameter. Inter- and intra-session repeatability were estimated for the optimal protocol by repeat scanning in 10 participants.

RESULTS

Final echo time was found to have no significant effect on OEF. The effect of voxel size was significant, with larger voxel sizes underestimating OEF, depending on the proximity of the vein to the superficial surface of the brain and on vein diameter. The last analysis step of estimating vein OEF values from susceptibility images had the largest dependency on voxel size. Inter-session coefficients of variation on OEF estimates of between 5.2% and 8.7% are reported, depending on the vein.

CONCLUSION

QSM acquisition times can be minimized by reducing the final echo time but an isotropic voxel size no larger than 1 mm is needed to accurately estimate OEF in most medium/large veins in the brain. Such acquisitions can be achieved in under 4 min.

Noninvasive quantification of oxygen saturation in the portal and hepatic veins in healthy mice and those with colorectal liver metastases using QSM MRI

PURPOSE

This preclinical study investigated the use of QSM MRI to noninvasively measure venous oxygen saturation (SvO2) in the hepatic and portal veins.

METHODS

QSM data were acquired from a cohort of healthy mice (n = 10) on a 9.4 Tesla MRI scanner under normoxic and hyperoxic conditions. Susceptibility was measured in the portal and hepatic veins and used to calculate SvO2 in each vessel under each condition. Blood was extracted from the inferior vena cava of 3 of the mice under each condition, and SvO2 was measured with a blood gas analyzer for comparison. QSM data were also acquired from a cohort of mice bearing liver tumors under normoxic conditions. Susceptibility was measured, and SvO2 calculated in the portal and hepatic veins and compared to the healthy mice. Statistical significance was assessed using a Wilcoxon matched-pairs signed rank test (normoxic vs. hyperoxic) or a Mann-Whitney test (healthy vs. tumor bearing).

RESULTS

SvO2 calculated from QSM measurements in healthy mice under hyperoxia showed significant increases of 15% in the portal vein (P < 0.05) and 21% in the hepatic vein (P < 0.01) versus normoxia. These values agreed with inferior vena cava measurements from the blood gas analyzer (26% increase). SvO2 in the hepatic vein was significantly lower in the colorectal liver metastases cohort (30% ± 11%) than the healthy mice (53% ± 17%) (P < 0.05); differences in the portal vein were not significant.

CONCLUSION

QSM is a feasible tool for noninvasively measuring SvO2 in the liver and can detect differences due to increased oxygen consumption in livers bearing colorectal metastases.

Patents on Quantitative Susceptibility Mapping (QSM) of Tissue Magnetism

BACKGROUND

Quantitative susceptibility mapping (QSM) depicts biodistributions of tissue magnetic susceptibility sources, including endogenous iron and calcifications, as well as exogenous paramagnetic contrast agents and probes. When comparing QSM with simple susceptibility weighted MRI, QSM eliminates blooming artifacts and shows reproducible tissue susceptibility maps independent of field strength and scanner manufacturer over a broad range of image acquisition parameters. For patient care, QSM promises to inform diagnosis, guide surgery, gauge medication, and monitor drug delivery. The Bayesian framework using MRI phase data and structural prior knowledge has made QSM sufficiently robust and accurate for routine clinical practice.

OBJECTIVE

To address the lack of a summary of US patents that is valuable for QSM product development and dissemination into the MRI community.

METHOD

We searched the USPTO Full-Text and Image Database for patents relevant to QSM technology innovation. We analyzed the claims of each patent to characterize the main invented method and we investigated data on clinical utility.

RESULTS

We identified 17 QSM patents; 13 were implemented clinically, covering various aspects of QSM technology, including the Bayesian framework, background field removal, numerical optimization solver, zero filling, and zero-TE phase.

CONCLUSION

Our patent search identified patents that enable QSM technology for imaging the brain and other tissues. QSM can be applied to study a wide range of diseases including neurological diseases, liver iron disorders, tissue ischemia, and osteoporosis. MRI manufacturers can develop QSM products for more seamless integration into existing MRI scanners to improve medical care.

A robust multi-scale approach to quantitative susceptibility mapping

Quantitative Susceptibility Mapping (QSM), best known as a surrogate for tissue iron content, is becoming a highly relevant MRI contrast for monitoring cellular and vascular status in aging, addiction, traumatic brain injury and, in general, a wide range of neurological disorders. In this study we present a new Bayesian QSM algorithm, named Multi-Scale Dipole Inversion (MSDI), which builds on the nonlinear Morphology-Enabled Dipole Inversion (nMEDI) framework, incorporating three additional features: (i) improved implementation of Laplace's equation to reduce the influence of background fields through variable harmonic filtering and subsequent deconvolution, (ii) improved error control through dynamic phase-reliability compensation across spatial scales, and (iii) scalewise use of the morphological prior. More generally, this new pre-conditioned QSM formalism aims to reduce the impact of dipole-incompatible fields and measurement errors such as flow effects, poor signal-to-noise ratio or other data inconsistencies that can lead to streaking and shadowing artefacts. In terms of performance, MSDI is the first algorithm to rank in the top-10 for all metrics evaluated in the 2016 QSM Reconstruction Challenge. It also demonstrated lower variance than nMEDI and more stable behaviour in scan-rescan reproducibility experiments for different MRI acquisitions at 3 and 7 Tesla. In the present work, we also explored new forms of susceptibility MRI contrast making explicit use of the differential information across spatial scales. Specifically, we show MSDI-derived examples of: (i) enhanced anatomical detail with susceptibility inversions from short-range dipole fields (hereby referred to as High-Pass Susceptibility Mapping or HPSM), (ii) high specificity to venous-blood susceptibilities for highly regularised HPSM (making a case for MSDI-based Venography or VenoMSDI), (iii) improved tissue specificity (and possibly statistical conditioning) for Macroscopic-Vessel Suppressed Susceptibility Mapping (MVSSM), and (iv) high spatial specificity and definition for HPSM-based Susceptibility-Weighted Imaging (HPSM-SWI) and related intensity projections.

Lesion magnetic susceptibility response to hyperoxic challenge: A biomarker for malignant brain tumor microenvironment?

BACKGROUND AND PURPOSE

Quantitative susceptibility mapping has been previously used to differentiate lesions in patients with brain tumors. The aim of this work was to characterize the response of magnetic susceptibility differences in malignant brain tumors and surrounding edema to hyperoxic and hypercapnic respiratory challenges.

METHODS

Images of malignant brain tumor patients (2 glioblastoma multiforme, 2 anaplastic astrocytoma, 1 brain metastasis) with clinical MRI exams (contrast-enhanced T1w) were acquired at 3T. 3D multi-gradient-echo data sets were acquired while the patients inhaled medical-air (21% O2), oxygen (100% O2), and carbogen (95% O2, 5% CO2). Susceptibility maps were generated from real and imaginary data. Regions of interest were analyzed with respect to respiration-gas-induced susceptibility changes.

RESULTS

Contrast-enhancing tumor regions with high baseline magnetic susceptibility exhibited a marked susceptibility reduction under hyperoxic challenges, with a stronger effect (-0.040 to -0.100ppm) under hypercapnia compared to hyperoxia (-0.010 to -0.067ppm). In contrast, regions attributed to necrotic tissue and to edema showed smaller changes of opposite sign, i.e. paramagnetic shift. There was a correlation between malignant tumor tissue magnetic susceptibility at baseline under normoxia and the corresponding susceptibility reduction under hypercapnia and - to a lesser degree - under hyperoxia.

CONCLUSION

In this small cohort of analysis, quantification of susceptibility changes in response to respiratory challenges allowed a complementary, functional differentiation of tumorous sub-regions. Those changes, together with the correlations observed between baseline susceptibility under normoxia and susceptibility reduction with challenges, could prove helpful for a non-invasive characterization of local tumor microenvironment.

Changes in brain iron concentration after exposure to high-altitude hypoxia measured by quantitative susceptibility mapping

Hypoxia can induce physiological changes. This study aims to explore effects of high-altitude (HA) hypoxia on cerebral iron concentration. Twenty-nine healthy sea-level participants were tested shortly before and after approximately 4-week adaptation to the HA environment at fQinghai-Tibet Plateau (4200m), and were re-investigated after re-adaptation to the sea-level environment one year later. Iron concentration was quantified with quantitative susceptibility mapping (QSM), and the results were compared with transverse relaxation rate (R^*₂) measurements. The variations of magnetic susceptibility indicate that the iron concentration in gray matter regions, especially in basal ganglia, including caudate nucleus, putamen, globus pallidus and substantia nigra, increases significantly after HA exposure. This increase appears consistent with the conclusion from R^*₂ value variations. However, unlike QSM, the R^*₂ value fails to demonstrate the statistical difference of iron content in red nucleus. The re-investigation results show that most variations are recovered after sea-level re-adaptation for one year. Additionally, hemisphere- and gender-related differences in iron concentration changes were analyzed among cerebral regions. The results show greater possibilities in the right hemisphere and females. Further studies based on diffusion tensor imaging (DTI) suggest that the fractional anisotropy increases and the mean diffusivity decreases after HA exposure in six deep gray matter nuclei, with linear dependence on iron concentration only in putamen. In conclusion, the magnetic susceptibility value can serve as a quantitative marker of brain iron, and variations of regional susceptibility reported herein indicate that HA hypoxia can result in significant iron deposition in most deep gray matter regions. Additionally, the linear dependence of DTI metrics on iron concentration in putamen indicates a potential relationship between ferritin and water diffusion.

Probing neuronal activation by functional quantitative susceptibility mapping under a visual paradigm: A group level comparison with BOLD fMRI and PET

Dynamic changes of brain-tissue magnetic susceptibility provide the basis for functional MR imaging (fMRI) via T2*-weighted signal-intensity modulations. Promising initial work on a detection of neuronal activity via quantitative susceptibility mapping (fQSM) has been published but consistently reported on ill-understood positive and negative activation patterns (Balla et al., 2014; Chen and Calhoun, 2015a). We set out to (i) demonstrate that fQSM can exploit established fMRI data acquisition and processing methods and to (ii) better describe aspects of the apparent activation patterns using fMRI and PET as standards of reference. Under a standardized visual-stimulation paradigm PET and 3-T gradient-echo EPI-based fQSM, fMRI data from 9 healthy volunteers were acquired and analyzed by means of Independent Component Analysis (ICA) at subject level and, for the first time, at group level. Numbers of activated (z-score>2.0) voxels were counted and their mean z-scores calculated in volumes of interest (occipital lobe (Nocc_lobe), segmented occipital gray-matter (NGM_occ_lobe), large veins (Nveins)), and in occipital-lobe voxels commonly activated in fQSM and fMRI component maps. Common but not entirely congruent regions of apparent activation were found in the occipital lobe in z-score maps from all modalities, fQSM, fMRI and PET, with distinct BOLD-negatively correlated regions in fQSM data. At subject-level, Nocc_lobe, NGM_occ_lobe and their mean z-scores were significantly smaller in fQSM than in fMRI, but their ratio, NGM_occ_lobe/Nocc_lobe, was comparable. Nveins did not statistically differ and the ratio Nveins/NGM_occ_lobe as well as the mean z-scores were higher for fQSM than for fMRI. In veins and immediate vicinity, z-score maps derived from both phase and fQSM-data showed positive and negative lobes resembling dipole shapes in simulated field and phase maps with no correlate in fMRI or PET data. Our results show that standard fMRI tools can directly be used for fQSM processing, and suggest that fQSM may have the potential to detect gray-matter activation distant from large veins, to improve detection of veins with stimulus-induced venous oxygen saturation (SvO2) variations, and to better localize areas of activation. However, our results seem to clearly expose issues that phenomenologically resemble an incomplete dipolar inversion and that need to be subject to further investigation.