Baseline oxygenation in the brain: Correlation between respiratory-calibration and susceptibility methods

Validation of a new 3D quantitative BOLD based cerebral oxygen extraction mapping

Quantitative BOLD (qBOLD) MRI allows evaluation of oxidative metabolism of the brain based purely on an endogenous contrast mechanism. The method quantifies deoxygenated blood volume (DBV) and hemoglobin oxygen saturation level of venous blood (Yv), yielding oxygen extraction fraction (OEF), and along with a separate measurement of cerebral blood flow, cerebral metabolic rate of oxygen (CMRO2) maps. Here, we evaluated our recently reported 3D qBOLD method that rectifies a number of deficiencies in prior qBOLD approaches in terms of repeat reproducibility and sensitivity to hypercapnia on the metabolic parameters, and in comparison to dual-gas calibrated BOLD (cBOLD) MRI for determining resting-state oxygen metabolism. Results suggested no significant difference between test-retest qBOLD scans in either DBV and OEF. Exposure to hypercapnia yielded group averages of 38 and 28% for OEF and 151 and 146 µmol/min/100 g for CMRO2 in gray matter at baseline and hypercapnia, respectively. The decrease of OEF during hypercapnia was significant (p ≪ 0.01), whereas CMRO2 did not change significantly (p = 0.25). Finally, baseline OEF (37 vs. 39%) and CMRO2 (153 vs. 145 µmol/min/100 g) in gray matter using qBOLD and dual-gas cBOLD were found to be in good agreement with literature values, and were not significantly different from each other (p > 0.1).

Cerebral oxygen metabolism from MRI susceptibility

This article provides an overview of MRI methods exploiting magnetic susceptibility properties of blood to assess cerebral oxygen metabolism, including the tissue oxygen extraction fraction (OEF) and the cerebral metabolic rate of oxygen (CMRO2). The first section is devoted to describing blood magnetic susceptibility and its effect on the MRI signal. Blood circulating in the vasculature can have diamagnetic (oxyhemoglobin) or paramagnetic properties (deoxyhemoglobin). The overall balance between oxygenated and deoxygenated hemoglobin determines the induced magnetic field which, in turn, modulates the transverse relaxation decay of the MRI signal via additional phase accumulation. The following sections of this review then illustrate the principles underpinning susceptibility-based techniques for quantifying OEF and CMRO2. Here, it is detailed whether these techniques provide global (OxFlow) or local (Quantitative Susceptibility Mapping - QSM, calibrated BOLD - cBOLD, quantitative BOLD - qBOLD, QSM+qBOLD) measurements of OEF or CMRO2, and what signal components (magnitude or phase) and tissue pools they consider (intravascular or extravascular). Validations studies and potential limitations of each method are also described. The latter include (but are not limited to) challenges in the experimental setup, the accuracy of signal modeling, and assumptions on the measured signal. The last section outlines the clinical uses of these techniques in healthy aging and neurodegenerative diseases and contextualizes these reports relative to results from gold-standard PET.

Mapping oxidative metabolism in the human brain with calibrated fMRI in health and disease

Conventional functional MRI (fMRI) with blood-oxygenation level dependent (BOLD) contrast is an important tool for mapping human brain activity non-invasively. Recent interest in quantitative fMRI has renewed the importance of oxidative neuroenergetics as reflected by cerebral metabolic rate of oxygen consumption (CMR_O2) to support brain function. Dynamic CMR_O2 mapping by calibrated fMRI require multi-modal measurements of BOLD signal along with cerebral blood flow (CBF) and/or volume (CBV). In human subjects this "calibration" is typically performed using a gas mixture containing small amounts of carbon dioxide and/or oxygen-enriched medical air, which are thought to produce changes in CBF (and CBV) and BOLD signal with minimal or no CMR_O2 changes. However non-human studies have demonstrated that the "calibration" can also be achieved without gases, revealing good agreement between CMR_O2 changes and underlying neuronal activity (e.g., multi-unit activity and local field potential). Given the simpler set-up of gas-free calibrated fMRI, there is evidence of recent clinical applications for this less intrusive direction. This up-to-date review emphasizes technological advances for such translational gas-free calibrated fMRI experiments, also covering historical progression of the calibrated fMRI field that is impacting neurological and neurodegenerative investigations of the human brain.

Cerebral oxygen extraction fraction MRI: Techniques and applications

The human brain constitutes 2% of the body's total mass but uses 20% of the oxygen. The rate of the brain's oxygen utilization can be derived from a knowledge of cerebral blood flow and the oxygen extraction fraction (OEF). Therefore, OEF is a key physiological parameter of the brain's function and metabolism. OEF has been suggested to be a useful biomarker in a number of brain diseases. With recent advances in MRI techniques, several MRI-based methods have been developed to measure OEF in the human brain. These MRI OEF techniques are based on the T₂ of blood, the blood signal phase, the magnetic susceptibility of blood-containing voxels, the effect of deoxyhemoglobin on signal behavior in extravascular tissue, and the calibration of the BOLD signal using gas inhalation. Compared to ¹⁵ O PET, which is considered the "gold standard" for OEF measurement, MRI-based techniques are non-invasive, radiation-free, and are more widely available. This article provides a review of these emerging MRI-based OEF techniques. We first briefly introduce the role of OEF in brain oxygen homeostasis. We then review the methodological aspects of different categories of MRI OEF techniques, including their signal mechanisms, acquisition methods, and data analyses. The strengths and limitations of the techniques are discussed. Finally, we review key applications of these techniques in physiological and pathological conditions.

Comparison of MRI methods for measuring whole‐brain oxygen extraction fraction under different geometric conditions at 7T

Background and Purpose

Methods

Results

Conclusion

Investigating the effect of flow compensation and quantitative susceptibility mapping method on the accuracy of venous susceptibility measurement

Quantitative susceptibility mapping (QSM) is a promising non-invasive method for obtaining information relating to oxygen metabolism. However, the optimal acquisition sequence and QSM reconstruction method for reliable venous susceptibility measurements are unknown. Full flow compensation is generally recommended to correct for the influence of venous blood flow, although the effect of flow compensation on the accuracy of venous susceptibility values has not been systematically evaluated. In this study, we investigated the effect of different acquisition sequences, including different flow compensation schemes, and different QSM reconstruction methods on venous susceptibilities. Ten healthy subjects were scanned with five or six distinct QSM sequence designs using monopolar readout gradients and different flow compensation schemes. All data sets were processed using six different QSM pipelines and venous blood susceptibility was evaluated in whole-brain segmentations of the venous vasculature and single veins. The quality of vein segmentations and the accuracy of venous susceptibility values were analyzed and compared between all combinations of sequences and reconstruction methods. The influence of the QSM reconstruction method on average venous susceptibility values was found to be 2.7-11.6 times greater than the influence of the acquisition sequence, including flow compensation. The majority of the investigated QSM reconstruction methods tended to underestimate venous susceptibility values in the vein segmentations that were obtained. In summary, we found that multi-echo gradient-echo acquisition sequences without full flow compensation yielded venous susceptibility values comparable to sequences with full flow compensation. However, the QSM reconstruction method had a great influence on susceptibility values and thus needs to be selected carefully for accurate venous QSM.

Cerebral oxygen extraction fraction (OEF): Comparison of challenge-free gradient echo QSM+qBOLD (QQ) with 15O PET in healthy adults

We aimed to validate oxygen extraction fraction (OEF) estimations by quantitative susceptibility mapping plus quantitative blood oxygen-level dependence (QSM+qBOLD, or QQ) using ¹⁵O-PET. In ten healthy adult brains, PET and MRI were acquired simultaneously on a PET/MR scanner. PET was acquired using C[¹⁵O], O[¹⁵O], and H₂[¹⁵O]. Image-derived arterial input functions and standard models of oxygen metabolism provided quantification of PET. MRI included T1-weighted imaging, time-of-flight angiography, and multi-echo gradient-echo imaging that was processed for QQ. Region of interest (ROI) analyses compared PET OEF and QQ OEF. In ROI analyses, the averaged OEF differences between PET and QQ were generally small and statistically insignificant. For whole brains, the average and standard deviation of OEF was 32.8 ± 6.7% for PET; OEF was 34.2 ± 2.6% for QQ. Bland-Altman plots quantified agreement between PET OEF and QQ OEF. The interval between the 95% limits of agreement was 16.9 ± 4.0% for whole brains. Our validation study suggests that respiratory challenge-free QQ-OEF mapping may be useful for non-invasive clinical assessment of regional OEF impairment.

Temporal lobe epilepsy with or without hippocampal sclerosis: Structural and functional connectivity using advanced MRI techniques

BACKGROUND AND PURPOSE

The aim of this study was to investigate the differences in structural connectivity based on diffusion tensor imaging (DTI) and functional connectivity based on arterial spin labeling (ASL) MRI between temporal lobe epilepsy (TLE) patients with and without hippocampal sclerosis (HS).

METHODS

We enrolled 50 patients with TLE, including 25 patients with HS and 25 patients without HS, who underwent brain MRI, including DTI and ASL. We calculated the network parameters of structural connectivity based on DTI and functional connectivity based on ASL using a graph theoretical analysis. The parameters included global network measures (radius, diameter, characteristic path length, global efficiency, local efficiency, mean clustering coefficient, transitivity, assortative coefficient, and small-worldness index) and a local network measure (betweenness centrality).

RESULTS

The global and local network measures of structural connectivity were not different between TLE patients with and without HS. However, significant differences in functional connectivity existed between the two groups. The radius and diameter of the global network measures in the TLE patients with HS were significantly increased compared with those without HS (4.140 vs. 3.140, p = 0.045; 6.812 vs. 5.132, p = 0.049; respectively). No differences were detected between other global network measures of functional connectivity and local network measure.

CONCLUSIONS

Significant differences in global network measures of functional connectivity based on ASL existed between TLE patients with and without HS. These findings suggest that TLE patients with HS exhibit a more disconnected functional brain network than those without HS.

Temporal clustering, tissue composition, and total variation for mapping oxygen extraction fraction using QSM and quantitative BOLD

PURPOSE

To improve the accuracy of quantitative susceptibility mapping plus quantitative blood oxygen level-dependent magnitude (QSM+qBOLD or QQ) based mapping of oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO₂ ) using temporal clustering, tissue composition, and total variation (CCTV).

METHODS

Three-dimensional multi-echo gradient echo and arterial spin labeling images were acquired from 11 healthy subjects and 33 ischemic stroke patients. Diffusion-weighted imaging (DWI) was also obtained from patients. The CCTV mapping was developed for incorporating tissue-type information into clustering of the previous cluster analysis of time evolution (CAT) and applying total variation (TV). The QQ-based OEF and CMRO₂ were reconstructed with CAT, CAT+TV (CATV), and the proposed CCTV, and results were compared using region-of-interest analysis, Kruskal-Wallis test, and post hoc Wilcoxson rank sum test.

RESULTS

In simulation, CCTV provided more accurate and precise OEF than CAT or CATV. In healthy subjects, QQ-based OEF was less noisy and more uniform with CCTV than CAT. In subacute stroke patients, OEF with CCTV had a greater contrast-to-noise ratio between DWI-defined lesions and the unaffected contralateral side than with CAT or CATV: 1.9 ± 1.3 versus 1.1 ± 0.7 (P = .01) versus 0.7 ± 0.5 (P < .001).

CONCLUSION

The CCTV mapping significantly improves the robustness of QQ-based OEF against noise.

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.

Cerebral oxygen extraction fraction: Comparison of dual-gas challenge calibrated BOLD with CBF and challenge-free gradient echo QSM+qBOLD

PURPOSE

To compare cortical gray matter oxygen extraction fraction (OEF) estimated from 2 MRI methods: (1) the quantitative susceptibility mapping (QSM) plus quantitative blood oxygen level dependent imaging (qBOLD) (QSM+qBOLD or QQ), and (2) the dual-gas calibrated-BOLD (DGCB) in healthy subjects; and to investigate the validity of iso-cerebral metabolic rate of oxygen consumption assumption during hypercapnia using QQ.

METHODS

In 10 healthy subjects, 3 tesla MRI including a multi-echo gradient echo sequence at baseline and hypercapnia for QQ, as well as an EPI dual-echo pseudo-continuous arterial spin labeling for DGCB, were performed under a hypercapnic and a hyperoxic condition. OEFs from QQ and DGCB were compared using region of interest analysis and paired t test. For QQ, cerebral metabolic rate of oxygen consumption = cerebral blood flow*OEF*arterial oxygen content was generated for both baseline and hypercapnia, which were compared.

RESULTS

Average OEF in cortical gray matter across 10 subjects from QQ versus DGCB was 35.5 ± 6.7% versus 38.0 ± 9.1% (P = .49) at baseline and 20.7 ± 4.4% versus 28.4 ± 7.6% (P = .02) in hypercapnia: OEF in cortical gray matter was significantly reduced as measured in QQ (P < .01) and in DGCB (P < .01). Cerebral metabolic rate of oxygen consumption (in μmol O₂ /min/100 g) was 168.2 ± 54.1 at baseline from DGCB and was 153.1 ± 33.8 at baseline and 126.4 ± 34.2 (P < .01) in hypercapnia from QQ.

CONCLUSION

The differences in OEF obtained from QQ and DGCB are small and nonsignificant at baseline but are statistically significant during hypercapnia. In addition, QQ shows a cerebral metabolic rate of oxygen consumption decrease (17.4%) during hypercapnia.

Detection of hemodynamic impairment on 15O gas PET using visual assessment of arterial spin-labeling MR imaging in patients with moyamoya disease

It is unclear whether the visual assessment of noninvasive arterial spin labeling magnetic resonance imaging (ASL) can identify instances of hemodynamic compromise including an elevated oxygen extraction fraction (OEF) measured by ¹⁵O-gas positron emission tomography (PET). Here we evaluated the relationship between a four-point visual assessment system referred to as 'ASL scores' using ASL with two postlabeling delays (PLDs; 1525 ms and 2525 ms) and some quantitative hemodynamic parameters measured by PET. We retrospectively evaluated the cases of 18 Japanese patients with moyamoya disease who underwent ASL and PET. We compared the patients' regional ASL scores on two ASL images to the regional values of PET parameters, and we observed a significant trend in accord with the presumed clinical severity among all PET parameters and ASL scores (p < .003). The ASL score of the long PLD (2525 ms) showed the highest specificity (98.5%) for elevated OEF. Our results suggest that hemodynamic impairment (including elevated OEF) in patients with moyamoya disease may be grossly assessed by a visual assessment of noninvasive ASL images, which can be easily obtained in clinical settings.

Cerebral oxygenation during locomotion is modulated by respiration

In the brain, increased neural activity is correlated with increases of cerebral blood flow and tissue oxygenation. However, how cerebral oxygen dynamics are controlled in the behaving animal remains unclear. We investigated to what extent cerebral oxygenation varies during locomotion. We measured oxygen levels in the cortex of awake, head-fixed mice during locomotion using polarography, spectroscopy, and two-photon phosphorescence lifetime measurements of oxygen sensors. We find that locomotion significantly and globally increases cerebral oxygenation, specifically in areas involved in locomotion, as well as in the frontal cortex and the olfactory bulb. The oxygenation increase persists when neural activity and functional hyperemia are blocked, occurred both in the tissue and in arteries feeding the brain, and is tightly correlated with respiration rate and the phase of respiration cycle. Thus, breathing rate is a key modulator of cerebral oxygenation and should be monitored during hemodynamic imaging, such as in BOLD fMRI.

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.

Arterial Spin Labeling Reveals Disrupted Brain Networks and Functional Connectivity in Drug-Resistant Temporal Epilepsy

Resting-state networks (RSNs) and functional connectivity (FC) have been increasingly exploited for mapping brain activity and identifying abnormalities in pathologies, including epilepsy. The majority of studies currently available are based on blood-oxygenation-level-dependent (BOLD) contrast in combination with either independent component analysis (ICA) or pairwise region of interest (ROI) correlations. Despite its success, this approach has several shortcomings as BOLD is only an indirect and non-quantitative measure of brain activity. Conversely, promising results have recently been achieved by arterial spin labeling (ASL) MRI, primarily developed to quantify brain perfusion. However, the wide application of ASL-based FC has been hampered by its complexity and relatively low robustness to noise, leaving several aspects of this approach still largely unexplored. In this study, we firstly aimed at evaluating the effect of noise reduction on spatio-temporal ASL analyses and quantifying the impact of two ad-hoc processing pipelines (basic and advanced) on connectivity measures. Once the optimal strategy had been defined, we investigated the applicability of ASL for connectivity mapping in patients with drug-resistant temporal epilepsy vs. controls (10 per group), aiming at revealing between-group voxel-wise differences in each RSN and ROI-wise FC changes. We first found ASL was able to identify the main network (DMN) along with all the others generally detected with BOLD but never previously reported from ASL. For all RSNs, ICA-based denoising (advanced pipeline) allowed to increase their similarity with the corresponding BOLD template. ASL-based RSNs were visibly consistent with literature findings; however, group differences could be identified in the structure of some networks. Indeed, statistics revealed areas of significant FC decrease in patients within different RSNs, such as DMN and cerebellum (CER), while significant increases were found in some cases, such as the visual networks. Finally, the ROI-based analyses identified several inter-hemispheric dysfunctional links (controls > patients) mainly between areas belonging to the DMN, right-left thalamus and right-left temporal lobe. Conversely, fewer connections, predominantly intra-hemispheric, showed the opposite pattern (controls < patients). All these elements provide novel insights into the pathological modulations characterizing a "network disease" as epilepsy, shading light on the importance of perfusion-based approaches for identifying the disrupted areas and communications between brain regions.

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.

Pros and cons of ultra-high-field MRI/MRS for human application

Magnetic resonance imaging and spectroscopic techniques are widely used in humans both for clinical diagnostic applications and in basic research areas such as cognitive neuroimaging. In recent years, new human MR systems have become available operating at static magnetic fields of 7 T or higher (≥300 MHz proton frequency). Imaging human-sized objects at such high frequencies presents several challenges including non-uniform radiofrequency fields, enhanced susceptibility artifacts, and higher radiofrequency energy deposition in the tissue. On the other side of the scale are gains in signal-to-noise or contrast-to-noise ratio that allow finer structures to be visualized and smaller physiological effects to be detected. This review presents an overview of some of the latest methodological developments in human ultra-high field MRI/MRS as well as associated clinical and scientific applications. Emphasis is given to techniques that particularly benefit from the changing physical characteristics at high magnetic fields, including susceptibility-weighted imaging and phase-contrast techniques, imaging with X-nuclei, MR spectroscopy, CEST imaging, as well as functional MRI. In addition, more general methodological developments such as parallel transmission and motion correction will be discussed that are required to leverage the full potential of higher magnetic fields, and an overview of relevant physiological considerations of human high magnetic field exposure is provided.

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.

Assessing the repeatability of absolute CMRO2, OEF and haemodynamic measurements from calibrated fMRI

As energy metabolism in the brain is largely oxidative, the measurement of cerebral metabolic rate of oxygen consumption (CMRO2) is a desirable biomarker for quantifying brain activity and tissue viability. Currently, PET techniques based on oxygen isotopes are the gold standard for obtaining whole brain CMRO2 maps. Among MRI techniques that have been developed as an alternative are dual calibrated fMRI (dcFMRI) methods, which exploit simultaneous measurements of BOLD and ASL signals during a hypercapnic-hyperoxic experiment to modulate brain blood flow and oxygenation. In this study we quantified the repeatability of a dcFMRI approach developed in our lab, evaluating its limits and informing its application in studies aimed at characterising the metabolic state of human brain tissue over time. Our analysis focussed on the estimates of oxygen extraction fraction (OEF), cerebral blood flow (CBF), CBF-related cerebrovascular reactivity (CVR) and CMRO2 based on a forward model that describes analytically the acquired dual echo GRE signal. Indices of within- and between-session repeatability are calculated from two different datasets both at a bulk grey matter and at a voxel-wise resolution and finally compared with similar indices obtained from previous MRI and PET measurements. Within- and between-session values of intra-subject coefficient of variation (CVintra) calculated from bulk grey matter estimates 6.7 ± 6.6% (mean ± std.) and 10.5 ± 9.7% for OEF, 6.9 ± 6% and 5.5 ± 4.7% for CBF, 12 ± 9.7% and 12.3 ± 10% for CMRO2. Coefficient of variation (CV) and intraclass correlation coefficient (ICC) maps showed the spatial distribution of the repeatability metrics, informing on the feasibility limits of the method. In conclusion, results show an overall consistency of the estimated physiological parameters with literature reports and a satisfactory level of repeatability considering the higher spatial sensitivity compared to other MRI methods, with varied performance depending on the specific parameter under analysis, on the spatial resolution considered and on the study design.

Calibrated fMRI for mapping absolute CMRO2: Practicalities and prospects

Functional magnetic resonance imaging (fMRI) is an essential workhorse of modern neuroscience, providing valuable insight into the functional organisation of the brain. The physiological mechanisms underlying the blood oxygenation level dependent (BOLD) effect are complex and preclude a straightforward interpretation of the signal. However, by employing appropriate calibration of the BOLD signal, quantitative measurements can be made of important physiological parameters including the absolute rate of cerebral metabolic oxygen consumption or oxygen metabolism (CMRO₂) and oxygen extraction (OEF). The ability to map such fundamental parameters has the potential to greatly expand the utility of fMRI and to broaden its scope of application in clinical research and clinical practice. In this review article we discuss some of the practical issues related to the calibrated-fMRI approach to the measurement of CMRO₂. We give an overview of the necessary precautions to ensure high quality data acquisition, and explore some of the pitfalls and challenges that must be considered as it is applied and interpreted in a widening array of diseases and research questions.

Vascular and Tissue Changes of Magnetic Susceptibility in the Mouse Brain After Transient Cerebral Ischemia

Quantitative susceptibility mapping (QSM) has been recently introduced as a novel MRI post-processing technique of gradient recalled echo (GRE) data. QSM is useful in depicting both brain anatomy and for detecting abnormalities. Its utility in the context of ischemic stroke has, however, not been extensively characterized so far. In this study, we explored the potential of QSM to characterize vascular and tissue changes in the transient middle cerebral artery occlusion (tMCAO) mouse model of cerebral ischemia. We acquired GRE data of mice brains at different time points after tMCAO, from which we computed QSM and MR frequency maps, and compared these maps with diffusion imaging and multi-slice multi-echo imaging data acquired in the same animals. Prominent vessels with increased magnetic susceptibility were visible surrounding the lesion on both frequency and magnetic susceptibility maps at all time points (mostly visible at > 12 h after reperfusion). Immunohistochemistry revealed the presence of compressed capillaries and dilated larger vessels, suggesting that the appearance of prominent vessels after reestablishment of reperfusion may serve compensatory purposes. In addition, on both contrast maps, tissue regions of decreased magnetic susceptibility were observed at 24 and 48 h after reperfusion that were distinctly different from the lesions seen on maps of the apparent diffusion coefficient and T₂ relaxation time constant. Since QSM can be extracted as an add-on from GRE data and thus requires no additional acquisition time in the course of acute stroke MRI examination, it may provide unique and complementary information during the course of acute stroke MRI examinations.

Functional oxygen extraction fraction (OEF) imaging with turbo gradient spin echo QUIXOTIC (Turbo QUIXOTIC)

PURPOSE

QUantitative Imaging of eXtraction of Oxygen and TIssue Consumption (QUIXOTIC) is a recent technique that measures voxel-wise oxygen extraction fraction (OEF) but suffers from long scan times, limiting its application. We implemented multiecho QUIXOTIC dubbed turbo QUIXOTIC (tQUIXOTIC) that reduces scan time eightfold and then applied it in functional MRI.

METHODS

tQUIXOTIC utilizes a novel turbo gradient spin echo readout enabling measurement of venular blood transverse relaxation rate in a single tag-control acquisition. Using tQUIXOTIC, we estimated cortical gray matter (GM) OEF, created voxel-by-voxel GM OEF maps, and quantified changes in visual cortex OEF during a blocked design flashing checkerboard visual stimulus. Contamination from cerebrospinal fluid partial volume averaging was estimated and corrected.

RESULTS

The average cortical GM OEF was estimated as 0.38 ± 0.06 (n = 8) using a 3.4-min acquisition. The average OEF in the visual cortex was estimated as 0.43 ± 0.04 at baseline and 0.35 ± 0.05 during activation, with an average %ΔOEF of -20%. These values are consistent with those of past studies.

CONCLUSION

tQUIXOTIC successfully estimated cortical GM OEF in clinical scan times and detected changes in OEF during blocked design visual stimulation. tQUIXOTIC will be useful to monitor regional OEF clinically and in blocked design or event-related functional MRI experiments. Magn Reson Med 79:2713-2723, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Susceptibility-Based Neuroimaging: Standard Methods, Clinical Applications, and Future Directions

The evaluation of neuropathologies using MRI methods that leverage tissue susceptibility have become standard practice, especially to detect blood products or mineralization. Additionally, emerging MRI techniques have the ability to provide new information based on tissue susceptibility properties in a robust and quantitative manner. This paper discusses these advanced susceptibility imaging techniques and their clinical applications.

MRI-based methods for quantification of the cerebral metabolic rate of oxygen

The brain depends almost entirely on oxidative metabolism to meet its significant energy requirements. As such, the cerebral metabolic rate of oxygen (CMRO2) represents a key measure of brain function. Quantification of CMRO2 has helped elucidate brain functional physiology and holds potential as a clinical tool for evaluating neurological disorders including stroke, brain tumors, Alzheimer's disease, and obstructive sleep apnea. In recent years, a variety of magnetic resonance imaging (MRI)-based CMRO2 quantification methods have emerged. Unlike positron emission tomography - the current "gold standard" for measurement and mapping of CMRO2 - MRI is non-invasive, relatively inexpensive, and ubiquitously available in modern medical centers. All MRI-based CMRO2 methods are based on modeling the effect of paramagnetic deoxyhemoglobin on the magnetic resonance signal. The various methods can be classified in terms of the MRI contrast mechanism used to quantify CMRO2: T2*, T2', T2, or magnetic susceptibility. This review article provides an overview of MRI-based CMRO2 quantification techniques. After a brief historical discussion motivating the need for improved CMRO2 methodology, current state-of-the-art MRI-based methods are critically appraised in terms of their respective tradeoffs between spatial resolution, temporal resolution, and robustness, all of critical importance given the spatially heterogeneous and temporally dynamic nature of brain energy requirements.

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.