Quantification of brain oxygen extraction fraction using QSM and a hyperoxic challenge

Quantitative susceptibility mapping through model-based deep image prior (MoDIP)

The data-driven approach of supervised learning methods has limited applicability in solving dipole inversion in Quantitative Susceptibility Mapping (QSM) with varying scan parameters across different objects. To address this generalization issue in supervised QSM methods, we propose a novel training-free model-based unsupervised method called MoDIP (Model-based Deep Image Prior). MoDIP comprises a small, untrained network and a Data Fidelity Optimization (DFO) module. The network converges to an interim state, acting as an implicit prior for image regularization, while the optimization process enforces the physical model of QSM dipole inversion. Experimental results demonstrate MoDIP's excellent generalizability in solving QSM dipole inversion across different scan parameters. It exhibits robustness against pathological brain QSM, achieving over 32 % accuracy improvement than supervised deep learning methods. It is also 33 % more computationally efficient and runs 4 times faster than conventional DIP-based approaches, enabling 3D high-resolution image reconstruction in under 4.5 min.

Multi-Echo Complex Quantitative Susceptibility Mapping and Quantitative Blood Oxygen Level-Dependent Magnitude (mcQSM + qBOLD or mcQQ) for Oxygen Extraction Fraction (OEF) Mapping

Oxygen extraction fraction (OEF), the fraction of oxygen that tissue extracts from blood, is an essential biomarker used to directly assess tissue viability and function in neurologic disorders. In ischemic stroke, for example, increased OEF can indicate the presence of penumbra-tissue with low perfusion yet intact cellular integrity-making it a primary therapeutic target. However, practical OEF mapping methods are not currently available in clinical settings, owing to the impractical data acquisitions in positron emission tomography (PET) and the limitations of existing MRI techniques. Recently, a novel MRI-based OEF mapping technique, termed QQ, was proposed. It shows high potential for clinical use by utilizing a routine sequence and removing the need for impractical multiple gas inhalations. However, QQ relies on the assumptions of Gaussian noise in susceptibility and multi-echo gradient echo (mGRE) magnitude signals for OEF estimation. This assumption is unreliable in low signal-to-noise ratio (SNR) regions like disease-related lesions, risking inaccurate OEF estimation and potentially impacting clinical decisions. Addressing this, our study presents a novel multi-echo complex QQ (mcQQ) that models realistic Gaussian noise in mGRE complex signals. We implemented mcQQ using a deep learning framework (mcQQ-NET) and compared it with the existing QQ-NET in simulations, ischemic stroke patients, and healthy subjects, using identical training and testing datasets and schemes. In simulations, mcQQ-NET provided more accurate OEF than QQ-NET. In the subacute stroke patients, mcQQ-NET showed a lower average OEF ratio in lesions relative to unaffected contralateral normal tissue than QQ-NET. In the healthy subjects, mcQQ-NET provided uniform OEF maps, similar to QQ-NET, but without unrealistically high OEF outliers in areas of low SNR, such as SNR ≤ 15 (dB). Therefore, mcQQ-NET improves OEF accuracy by more accurately reflecting realistic Gaussian noise in complex mGRE signals. Its enhanced sensitivity to OEF abnormalities, based on more realistic biophysics modeling, suggests that mcQQ-NET has potential for investigating tissue variability in neurologic disorders.

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.

Tale of Two Cities: narrative review of oxygen

The human brain contributes 2% of the body weight yet receives 15% of cardiac output and demands a constant supply of oxygen (O ₂) and nutrients to meet its metabolic needs. Cerebral autoregulation is responsible for maintaining a constant cerebral blood flow that provides the supply of oxygen and maintains the energy storage capacity. We selected oxygen administration-related studies published between 1975-2021 that included meta-analysis, original research, commentaries, editorial, and review articles. In the present narrative review, several important aspects of the oxygen effects on brain tissues and cerebral autoregulation are discussed, as well the role of exogenous O ₂ administration in patients with chronic ischemic cerebrovascular disease: We aimed to revisit the utility of O ₂ administration in pathophysiological situations whether or not being advantageous. Indeed, a compelling clinical and experimental body of evidence questions the utility of routine oxygen administration in acute and post-recovery brain ischemia, as evident by studies in neurophysiology imaging. While O ₂ is still part of common clinical practice, it remains unclear whether its routine use is safe.

Development of a novel 10-echo multi-contrast sequence based on EPIK to deliver simultaneous quantification of T2 and T2 * with application to oxygen extraction fraction

PURPOSE

The simultaneous quantification of T₂ and T₂ * maps based on fast sequences for combined GE and SE acquisition has rekindled research and clinical interest by offering a wide range of attractive applications, e.g., dynamic tracking of oxygen extraction fraction (OEF). However, previously published methods based on EPI-readouts have been hindered by resolution and the number of acquired echoes.

METHODS

This work presents a novel 10-echo GE-SE EPIK (EPI with keyhole) sequence for the rapid quantification of T₂ '. T₂ /T₂ * maps from the GE-SE EPIK sequence were validated using three phantoms and 15 volunteers at 3T. The incorporation of a sliding window approach, combined with the full sampling of the k-space center inherent to EPIK, enables a high effective temporal resolution. That is, for an eight-slice breath-hold experiment, a temporal sampling rate of eight reconstructed slices per 1.1 s.

RESULTS

In comparison with repeated single-echo SE, multi-echo GE, and spectroscopy methods, the GE-SE EPIK sequence shows good agreement in quantifying T₂ /T₂ * values, while the gray matter/white matter separation yielded the expected contrast differentiation. The OEF was calculated with a view to an initial application with clinical relevance, producing results comparable to those in the literature and with good sensitivity in breath-hold experiments.

CONCLUSIONS

GE-SE EPIK provides increased resolution and more echoes, including two SEs, than comparable sequences. Moreover, GE-SE EPIK achieves this within an acquisition time of 57 s for 20 slices (matrix size = 128×128; FOV = 24 cm) and with a reasonably short TE for the final echo (114 ms). The sequence can dynamically track OEF changes in a breath-hold experiment.

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.

A novel gradient echo data based vein segmentation algorithm and its application for the detection of regional cerebral differences in venous susceptibility

Accurate segmentation of cerebral venous vasculature from gradient echo data is of central importance in several areas of neuroimaging such as for the susceptibility-based assessment of brain oxygenation or planning of electrode placement in deep brain stimulation. In this study, a vein segmentation algorithm for single- and multi-echo gradient echo data is proposed. First, susceptibility maps, true susceptibility-weighted images, and, in the multi-echo case, R₂^* maps were generated from the gradient echo data. These maps were filtered with an inverted Hamming filter to suppress background contrast as well as artifacts from field inhomogeneities at the brain boundaries. A shearlet-based scale-wise representation was generated to calculate a vesselness function and to generate segmentations based on local thresholding. The accuracy of the proposed algorithm was evaluated for different echo times and image resolutions using a manually generated reference segmentation and two vein segmentation algorithms (Frangi vesselness-based, recursive vesselness filter) as a reference with the Dice and Cohen's coefficients as well as the modified Hausdorff distance. The Frangi-based and recursive vesselness filter methods were significantly outperformed with regard to all error metrics. Applying the algorithm, susceptibility differences likely related to differences in blood oxygenation between superficial and deep venous territories could be demonstrated.

Recovering SWI-filtered phase data using deep learning

PURPOSE

To develop a deep neural network to recover filtered phase from clinical MR phase images to enable the computation of QSMs.

METHODS

Eighteen deep learning networks were trained to recover combinations of 13 SWI phase-filtering pipelines. SWI-filtered data were computed offline from five multiorientation, multiecho MRI scans yielding 132 3D volumes (118/7/7 training/validation/testing). Two experiments were conducted to show the efficacy of the networks. First, using QSM processing, local fields were computed from the raw phase and subsequently filtered using the SWI-filtering pipelines. The networks were then trained to invert the filtering operation. Second, the trained networks were fine-tuned to recover unfiltered local fields from filtered local fields computed by applying QSM processing to the SWI-filtered phase. Susceptibility maps were computed from the recovered fields and compared with gold standard multiple orientation sampling reconstructions.

RESULTS

Susceptibility maps computed from the raw phase using standard QSM processing have a normalized root mean square error (NRMSE) of 0.732 ± 0.095. Susceptibility maps computed from the recovered phase obtained NRMSEs of 0.725 ± 0.095. The network trained using all 13 processing methods generalized well, obtaining NRMSEs of 0.725 ± 0.89 on filters it has not seen, while matching the reconstruction accuracy of networks trained to recover a single filter.

CONCLUSION

It is feasible to recover SWI-filtered phase using deep learning. QSM can be computed from the recovered phase from SWI acquisition with comparable accuracy to standard QSM processing.

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.