Quantitative mapping of cerebral metabolic rate of oxygen (CMRO2 ) using quantitative susceptibility mapping (QSM)

Vessel-specific quantification of cerebral venous oxygenation with velocity-encoding preparation and rapid acquisition

PURPOSE

Non-invasive measurement of cerebral venous oxygenation (Yv) is of critical importance in brain diseases. The present work proposed a fast method to quantify regional Yv map for both large and small veins.

METHODS

A new sequence was developed, referred to as TRU-VERA (T2 relaxation under velocity encoding and rapid acquisition, which isolates blood spins from static tissue with velocity-encoding preparation, modulates the T2 weighting of venous signal with T2-preparation and utilizes a bSSFP readout to achieve fast acquisition with high resolution. The sequence was first optimized to achieve best sensitivity for both large and small veins, and then validated with TRUST (T2 relaxation under spin tagging), TRUPC (T2 relaxation under phase contrast), and accelerated TRUPC MRI. Regional difference of Yv was evaluated, and test-retest reproducibility was examined.

RESULTS

Optimal Venc was determined to be 3 cm/s, while recovery time and balanced SSFP flip angle within reasonable range had minimal effect on SNR efficiency. Venous T2 measured with TRU-VERA was highly correlated with T2 from TRUST (R2 = 0.90), and a conversion equation was established for further calibration to Yv. TRU-VERA sequences showed consistent Yv estimation with TRUPC (R2 = 0.64) and accelerated TRUPC (R2 = 0.79). Coefficient of variation was 0.84% for large veins and 2.49% for small veins, suggesting an excellent test-retest reproducibility.

CONCLUSION

The proposed TRU-VERA sequence is a promising method for vessel-specific oxygenation assessment.

Recent Advances in MR Imaging-based Quantification of Brain Oxygen Metabolism

The metabolic rate of oxygen (MRO2) is fundamental to tissue metabolism. Determination of MRO2 demands knowledge of the arterio-venous difference in hemoglobin-bound oxygen concentration, typically expressed as oxygen extraction fraction (OEF), and blood flow rate (BFR). MRI is uniquely suited for measurement of both these quantities, yielding MRO2 in absolute physiologic units of µmol O2 min-1/100 g tissue. Two approaches are discussed, both relying on hemoglobin magnetism. Emphasis will be on cerebral oxygen metabolism expressed in terms of the cerebral MRO2 (CMRO2), but translation of the relevant technologies to other organs, including kidney and placenta will be touched upon as well. The first class of methods exploits the blood's bulk magnetic susceptibility, which can be derived from field maps. The second is based on measurement of blood water T2, which is modulated by diffusion and exchange in the local-induced fields within and surrounding erythrocytes. Some whole-organ methods achieve temporal resolution adequate to permit time-series studies of brain energetics, for instance, during sleep in the scanner with concurrent electroencephalogram (EEG) sleep stage monitoring. Conversely, trading temporal for spatial resolution has led to techniques for spatially resolved approaches based on quantitative blood oxygen level dependent (BOLD) or calibrated BOLD models, allowing regional assessment of vascular-metabolic parameters, both also exploiting deoxyhemoglobin paramagnetism like their whole-organ counterparts.

Unveiling metabolic patterns in dementia: Insights from high-resolution quantitative blood-oxygenation-level-dependent MRI

BACKGROUND

Oxygen extraction fraction (OEF) and deoxyhemoglobin (DoHb) levels reflect variations in cerebral oxygen metabolism in demented patients.

PURPOSE

Delineating the metabolic profiles evident throughout different phases of dementia necessitates an integrated analysis of OEF and DoHb levels. This is enabled by leveraging high-resolution quantitative blood oxygenation level dependent (qBOLD) analysis of magnitude images obtained from a multi-echo gradient-echo MRI (mGRE) scan performed on a 3.0 Tesla scanner.

METHODS

Achieving superior spatial resolution in qBOLD necessitates the utilization of an mGRE scan with only four echoes, which in turn limits the number of measurements compared to the parameters within the qBOLD model. Consequently, it becomes imperative to discard non-essential parameters to facilitate further analysis. This process entails transforming the qBOLD model into a format suitable for fitting the log-magnitude difference (L-MDif) profiles of the four echo magnitudes present in each brain voxel. In order to bolster spatial specificity, the log-difference qBOLD model undergoes refinement into a representative form, termed as r-qBOLD, particularly when applied to class-averaged L-MDif signals derived through k-means clustering of L-MDif signals from all brain voxels into a predetermined number of clusters. The agreement between parameters estimated using r-qBOLD for different cluster sizes is validated using Bland-Altman analysis, and the model's goodness-of-fit is evaluated using aχ 2 ${\chi ^2}$ -test. Retrospective MRI data of Alzheimer's disease (AD), mild cognitive impairment (MCI), and non-demented patients without neuropathological disorders, pacemakers, other implants, or psychiatric disorders, who completed a minimum of three visits prior to MRI enrolment, are utilized for the study.

RESULTS

Utilizing a cohort comprising 30 demented patients aged 65-83 years in stages 4-6 representing mild, moderate, and severe stages according to the clinical dementia rating (CDR), matched with an age-matched non-demented control group of 18 individuals, we conducted joint observations of OEF and DoHb levels estimated using r-qBOLD. The observations elucidate metabolic signatures in dementia based on OEF and DoHb levels in each voxel. Our principal findings highlight the significance of spatial patterns of metabolic profiles (metabolic patterns) within two distinct regimes: OEF levels exceeding the normal range (S1-regime), and OEF levels below the normal range (S2-regime). The S1-regime, accompanied by low DoHb levels, predominantly manifests in fronto-parietal and perivascular regions with increase in dementia severity. Conversely, the S2-regime, accompanied by low DoHb levels, is observed in medial temporal (MTL) regions. Other regions with abnormal metabolic patterns included the orbitofrontal cortex (OFC), medial-orbital prefrontal cortex (MOPFC), hypothalamus, ventro-medial prefrontal cortex (VMPFC), and retrosplenial cortex (RSP). Dysfunction in the OFC and MOPFC indicated cognitive and emotional impairment, while hypothalamic involvement potentially indicated preclinical dementia. Reduced metabolic activity in the RSP suggested early-stage AD related functional abnormalities.

CONCLUSIONS

Integrated analysis of OEF and DoHb levels using r-qBOLD reveals distinct metabolic signatures across dementia phases, highlighting regions susceptible to neuronal loss, vascular involvement, and preclinical indicators.

Reduced oxygen extraction fraction in deep cerebral veins associated with cognitive impairment in multiple sclerosis

Studying the relationship between cerebral oxygen utilization and cognitive impairment is essential to understanding neuronal functional changes in the disease progression of multiple sclerosis (MS). This study explores the potential of using venous susceptibility in internal cerebral veins (ICVs) as an imaging biomarker for cognitive impairment in relapsing-remitting MS (RRMS) patients. Quantitative susceptibility mapping derived from fully flow-compensated MRI phase data was employed to directly measure venous blood oxygen saturation levels (SvO2) in the ICVs. Results revealed a significant reduction in the susceptibility of ICVs (212.4 ± 30.8 ppb vs 239.4 ± 25.9 ppb) and a significant increase of SvO2 (74.5 ± 1.89% vs 72.4 ± 2.23%) in patients with RRMS compared with age- and sex-matched healthy controls. Both the susceptibility of ICVs (r = 0.508, p = 0.031) and the SvO2 (r = -0.498, p = 0.036) exhibited a moderate correlation with cognitive decline in these patients assessed by the Paced Auditory Serial Addition Test, while no significant correlation was observed with clinical disability measured by the Expanded Disability Status Scale. The findings suggest that venous susceptibility in ICVs has the potential to serve as a specific indicator of oxygen metabolism and cognitive function in RRMS. .

Plug-and-Play latent feature editing for orientation-adaptive quantitative susceptibility mapping neural networks

Quantitative susceptibility mapping (QSM) is a post-processing technique for deriving tissue magnetic susceptibility distribution from MRI phase measurements. Deep learning (DL) algorithms hold great potential for solving the ill-posed QSM reconstruction problem. However, a significant challenge facing current DL-QSM approaches is their limited adaptability to magnetic dipole field orientation variations during training and testing. In this work, we propose a novel Orientation-Adaptive Latent Feature Editing (OA-LFE) module to learn the encoding of acquisition orientation vectors and seamlessly integrate them into the latent features of deep networks. Importantly, it can be directly Plug-and-Play (PnP) into various existing DL-QSM architectures, enabling reconstructions of QSM from arbitrary magnetic dipole orientations. Its effectiveness is demonstrated by combining the OA-LFE module into our previously proposed phase-to-susceptibility single-step instant QSM (iQSM) network, which was initially tailored for pure-axial acquisitions. The proposed OA-LFE-empowered iQSM, which we refer to as iQSM+, is trained in a simulated-supervised manner on a specially-designed simulation brain dataset. Comprehensive experiments are conducted on simulated and in vivo human brain datasets, encompassing subjects ranging from healthy individuals to those with pathological conditions. These experiments involve various MRI platforms (3T and 7T) and aim to compare our proposed iQSM+ against several established QSM reconstruction frameworks, including the original iQSM. The iQSM+ yields QSM images with significantly improved accuracies and mitigates artifacts, surpassing other state-of-the-art DL-QSM algorithms. The PnP OA-LFE module's versatility was further demonstrated by its successful application to xQSM, a distinct DL-QSM network for dipole inversion. In conclusion, this work introduces a new DL paradigm, allowing researchers to develop innovative QSM methods without requiring a complete overhaul of their existing architectures.

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.

Reduced Oxygen Extraction Fraction in Deep Cerebral Veins Associated with Cognitive Impairment in Multiple Sclerosis

Studying the relationship between cerebral oxygen utilization and cognitive impairment is essential to understanding neuronal functional changes in the disease progression of multiple sclerosis (MS). This study explores the potential of using venous susceptibility in internal cerebral veins (ICVs) as an imaging biomarker for cognitive impairment in relapsing-remitting MS (RRMS) patients. Quantitative susceptibility mapping derived from fully flow-compensated MRI phase data was employed to directly measure venous blood oxygen saturation levels (SvO2) in the ICVs. Results revealed a significant reduction in the susceptibility of ICVs (212.4 ± 30.8 ppb vs 239.4 ± 25.9 ppb) and a significant increase of SvO2 (74.5 ± 1.89 % vs 72.4 ± 2.23 %) in patients with RRMS compared with age- and sex-matched healthy controls. Both the susceptibility of ICVs (r = 0.646, p = 0.004) and the SvO2 (r = -0.603, p = 0.008) exhibited a strong correlation with cognitive decline in these patients assessed by the Paced Auditory Serial Addition Test, while no significant correlation was observed with clinical disability measured by the Expanded Disability Status Scale. The findings suggest that venous susceptibility in ICVs has the potential to serve as a specific indicator of oxygen metabolism and cognitive function in RRMS.

Transient susceptibility imaging as a measure of hemodynamic compromise: A pilot study

BACKGROUND AND PURPOSE

The purpose of this study was to test the hypothesis that hemodynamically compromised brains exhibit transient changes in magnetic susceptibility throughout the cardiac cycle, and to model these changes using Linear System Theory to extract an index that reflects cerebrovascular reserve.

MATERIALS AND METHODS

Eleven patients with angiographically-confirmed intracranial atherosclerotic disease with >50% stenosis were imaged with susceptibility weighted, cardiac-gated single shot images of cerebral Oxygen Extraction Fraction (OEF) at different timepoints of the cardiac cycle. Cardiac gating of the OEF acquisition allowed interrogation of oxygenated blood and the detection of changes throughout the cardiac cycle. Independent component analysis (ICA) of raw k-space data across the cardiac phase allowed MRI signal decomposition into dynamic and static components for image reconstruction. An asymmetry index score of the resultant parametric images were compared to test the hypothesis that variation in hemoglobin-induced susceptibility across the cardiac cycle indeed reflects pathophysiology of cerebrovascular disease. A mathematical model was derived to parameterize physiologic changes induced by the presence of a hemodynamically significant stenosis in the brain as a tissue impulse response parameter (β).

RESULTS

OEF was elevated in the affected hemisphere (50.34 ± 12.13% vs 46.93 ± 12.34%), but failed to reach statistical significance (p < .0796). Transient changes in the OEF signal showed significant distinction between healthy and compromised tissue (0.56 ± 0.067 vs 0.44 ± 0.067, p < .019)). The derived tissue impulse response function was found to be significant as well (10.72 ± 3.48 10-3 ms-1, 9.69 ± 3.51 10-3 ms-1; p < .037).

CONCLUSION

In this pilot study, we found transient OEF and β to be significant predictors of hemispheric compromise.

Current therapeutic targets and multifaceted physiological impacts of caffeine

Caffeine, which shares consubstantial structural similarity with purine adenosine, has been demonstrated as a nonselective adenosine receptor antagonist for eliciting most of the biological functions at physiologically relevant dosages. Accumulating evidence supports caffeine's beneficial effects against different disorders, such as total cardiovascular diseases and type 2 diabetes. Conversely, paradoxical effects are also linked to caffeine ingestion in humans including hypertension-hypotension and tachycardia-bradycardia. These observations suggest the association of caffeine action with its ingested concentration and/or concurrent interaction with preferential molecular targets to direct explicit events in the human body. Thus, a coherent analysis of the functional targets of caffeine, relevant to normal physiology, and disease pathophysiology, is required to understand the pharmacology of caffeine. This review provides a broad overview of the experimentally validated targets of caffeine, particularly those of therapeutic interest, and the impacts of caffeine on organ-specific physiology and pathophysiology. Overall, the available empirical and epidemiological evidence supports the dose-dependent functional activities of caffeine and advocates for further studies to get insights into the caffeine-induced changes under specific conditions, such as asthma, DNA repair, and cancer, in view of its therapeutic applications.

Evaluation of whole-brain oxygen metabolism in Alzheimer's disease using QSM and quantitative BOLD

OBJECTIVE

The objective of this study was to evaluate the whole-brain pattern of oxygen extraction fraction (OEF), cerebral blood flow (CBF), and cerebral metabolic rate of oxygen consumption (CMRO2) perturbation in Alzheimer's disease (AD) and investigate the relationship between regional cerebral oxygen metabolism and global cognition.

METHODS

Twenty-six AD patients and 25 age-matched healthy controls (HC) were prospectively recruited in this study. Mini-Mental State Examination (MMSE) was used to evaluate cognitive status. We applied the QQ-CCTV algorithm which combines quantitative susceptibility mapping and quantitative blood oxygen level-dependent models (QQ) for OEF calculation. CBF map was computed from arterial spin labeling and CMRO2 was generated based on Fick's principle. Whole-brain and regional OEF, CBF, and CMRO2 analyses were performed. The associations between these measures in substructures of deep brain gray matter and MMSE scores were assessed.

RESULTS

Whole brain voxel-wise analysis showed that CBF and CMRO2 values significantly decreased in AD predominantly in the bilateral angular gyrus, precuneus gyrus and parieto-temporal regions. Regional analysis showed that CBF value decreased in the bilateral caudal hippocampus and left rostral hippocampus and CMRO2 value decreased in left caudal and rostral hippocampus in AD patients. Considering all subjects in the AD and HC groups combined, the mean CBF and CMRO2 values in the bilateral hippocampus positively correlated with the MMSE score.

CONCLUSION

CMRO2 mapping with the QQ-CCTV method - which is readily available in MR systems for clinical practice - can be a potential biomarker for AD. In addition, CMRO2 in the hippocampus may be a useful tool for monitoring cognitive impairment.

Toward vendor-independent measurement of cerebral venous oxygenation: Comparison of TRUST MRI across three major MRI manufacturers and association with end-tidal CO2

Cerebral venous oxygenation (Yv ) is a valuable biomarker for a variety of brain diseases. T2 relaxation under spin tagging (TRUST) MRI is a widely used method for Yv quantification. In this work, there were two main objectives. The first was to evaluate the reproducibility of TRUST Yv measurements across MRI scanners from different vendors. The second was to examine the correlation between Yv and end-tidal CO2 (EtCO2 ) in a multisite, multivendor setting and determine the usefulness of this correlation to account for variations in Yv caused by normal variations and physiological fluctuations. Standardized TRUST pulse sequences were implemented on three scanners from major MRI vendors (GE, Siemens, Philips). These scanners were located at two research institutions. Ten healthy subjects were scanned. On each scanner, the subject underwent two scan sessions, each of which included three TRUST scans, to evaluate the intrasession and intersession reproducibility of Yv . Each scanner was also equipped with a capnograph device to record the EtCO2 of the subject during the MRI scan. We found no significant bias in Yv measurements across the three scanners (P = 0.18). The measured Yv values on the three scanners were also strongly correlated with each other (intraclass correlation coefficients > 0.85, P < 0.001). The intrasession and intersession coefficients of variation of Yv were less than 4% and showed no significant difference among the scanners. In addition, our results revealed that (1) within the same subject, Yv increased with EtCO2 at a rate of 1.24 ± 0.17%/mmHg (P < 0.0001), and (2) across different subjects, individuals with a higher EtCO2 had a higher Yv , at a rate of 0.94 ± 0.36%/mmHg (P = 0.01). These results suggest that (1) the standardized TRUST sequences had similar accuracies and reproducibilities for the quantification of Yv across the scanners, and (2) recording of EtCO2 may be a useful complement to Yv measurement to account for CO2 -related physiological fluctuations in Yv in multisite, multivendor studies.

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.

Capillary transit time heterogeneity inhibits cerebral oxygen metabolism in patients with reduced cerebrovascular reserve capacity from steno-occlusive disease

The healthy cerebral perfusion demonstrates a homogenous distribution of capillary transit times. A disruption of this homogeneity may inhibit the extraction of oxygen. A high degree of capillary transit time heterogeneity (CTH) describes that some capillaries have very low blood flows, while others have excessively high blood flows and consequently short transit times. Very short transit times could hinder the oxygen extraction due to insufficient time for diffusion of oxygen into the tissue. CTH could be a consequence of cerebral vessel disease. We examined whether patients with cerebral steno-occlusive vessel disease demonstrate high CTH and if elevation of cerebral blood flow (CBF) by administration of acetazolamide (ACZ) increases the cerebral metabolic rate of oxygen (CMRO₂), or if some patients demonstrate reduced CMRO₂ related to detrimental CTH. Thirty-four patients and thirty-one healthy controls participated. Global CBF and CMRO₂ were acquired using phase-contrast MRI. Regional brain maps of CTH were acquired using dynamic contrast-enhanced MRI. Patients with impaired cerebrovascular reserve capacity demonstrated elevated CTH and a significant reduction of CMRO₂ after administration of ACZ, which could be related to high CTH. Impaired oxygen extraction from CTH could be a contributing part of the declining brain health observed in patients with cerebral vessel disease.

Measurement of Cerebral Oxygen Extraction Fraction Using Quantitative BOLD Approach: A Review

Quantification of brain oxygenation and metabolism, both of which are indicators of the level of brain activity, plays a vital role in understanding the cerebral perfusion and the pathophysiology of brain disorders. Magnetic resonance imaging (MRI), a widely used clinical imaging technique, which is very sensitive to magnetic susceptibility, has the possibility of substituting positron emission tomography (PET) in measuring oxygen metabolism. This review mainly focuses on the quantitative blood oxygenation level-dependent (qBOLD) method for the evaluation of oxygen extraction fraction (OEF) in the brain. Here, we review the theoretic basis of qBOLD, as well as existing acquisition and quantification methods. Some published clinical studies are also presented, and the pros and cons of qBOLD method are discussed as well.

Neurodegenerative disease of the brain: a survey of interdisciplinary approaches

Neurodegenerative diseases of the brain pose a major and increasing global health challenge, with only limited progress made in developing effective therapies over the last decade. Interdisciplinary research is improving understanding of these diseases and this article reviews such approaches, with particular emphasis on tools and techniques drawn from physics, chemistry, artificial intelligence and psychology.

Quantitative susceptibility mapping using plug-and-play alternating direction method of multipliers

Quantitative susceptibility mapping employs regularization to reduce artifacts, yet many recent denoisers are unavailable for reconstruction. We developed a plug-and-play approach to QSM reconstruction (PnP QSM) and show its flexibility using several patch-based denoisers. We developed PnP QSM using alternating direction method of multiplier framework and applied collaborative filtering denoisers. We apply the technique to the 2016 QSM Challenge and in 10 glioblastoma multiforme datasets. We compared its performance with four published QSM techniques and a multi-orientation QSM method. We analyzed magnetic susceptibility accuracy using brain region-of-interest measurements, and image quality using global error metrics. Reconstructions on glioblastoma data were analyzed using ranked and semiquantitative image grading by three neuroradiologist observers to assess image quality (IQ) and sharpness (IS). PnP-BM4D QSM showed good correlation (β = 0.84, R² = 0.98, p < 0.05) with COSMOS and no significant bias (bias = 0.007 ± 0.012). PnP-BM4D QSM achieved excellent quality when assessed using structural similarity index metric (SSIM = 0.860), high frequency error norm (HFEN = 58.5), cross correlation (CC = 0.804), and mutual information (MI = 0.475) and also maintained good conspicuity of fine features. In glioblastoma datasets, PnP-BM4D QSM showed higher performance (IQ_Grade = 2.4 ± 0.4, IS_Grade = 2.7 ± 0.3, IQ_Rank = 3.7 ± 0.3, IS_Rank = 3.9 ± 0.3) compared to MEDI (IQ_Grade = 2.1 ± 0.5, IS_Grade = 2.1 ± 0.6, IQ_Rank = 2.4 ± 0.6, IS_Rank = 2.9 ± 0.2) and FANSI-TGV (IQ_Grade = 2.2 ± 0.6, IS_Grade = 2.1 ± 0.6, IQ_Rank = 2.7 ± 0.3, IS_Rank = 2.2 ± 0.2). We illustrated the modularity of PnP QSM by interchanging two additional patch-based denoisers. PnP QSM reconstruction was feasible, and its flexibility was shown using several patch-based denoisers. This technique may allow rapid prototyping and validation of new denoisers for QSM reconstruction for an array of useful clinical applications.

Blood-brain barrier permeability in response to caffeine challenge

PURPOSE

Caffeine is known to alter brain perfusion by acting as an adenosine antagonist, but its effect on blood-brain barrier (BBB) permeability is not fully elucidated. This study aimed to dynamically monitor BBB permeability to water after a single dose of caffeine tablet using a non-contrast MRI technique.

METHODS

Ten young healthy volunteers who were not regular coffee drinkers were studied. The experiment began with a pre-caffeine measurement, followed by four measurements at the post-caffeine stage. Water-extraction-with-phase-contrast-arterial-spin-tagging (WEPCAST) MRI was used to assess the time dependence of BBB permeability to water following the ingestion of 200 mg caffeine. Other cerebral physiological parameters including cerebral blood flow (CBF), venous oxygenation (Yv ), and cerebral metabolic rate of oxygen (CMRO2 ) were also examined. The relationships between cerebral physiological parameters and time were studied with mixed-effect models.

RESULTS

It was found that, after caffeine ingestion, CBF and Yv showed a time-dependent decrease (p < 0.001), while CMRO2 did not change significantly. The fraction of arterial water crossing the BBB (E) showed a significant increase (p < 0.001). In contrast, the permeability-surface-area product (PS), i.e., BBB permeability to water, remained constant (p = 0.94). Additionally, it was observed that changes in physiological parameters were non-linear with regard to time and occurred at as early as 9 min after caffeine tablet ingestion.

CONCLUSION

These results suggest an unchanged BBB permeability despite alterations in perfusion during a vasoconstrictive caffeine challenge.

Oxygen extraction fraction (OEF) assesses cerebral oxygen metabolism of deep gray matter in patients with pre-eclampsia

OBJECTIVES

The objective of this study was to compare oxygen extraction fraction (OEF) values in the deep gray matter (GM) of pre-eclampsia (PE) patients, pregnant healthy controls (PHCs), and non-pregnant healthy controls (NPHCs) to explore their brain oxygen metabolism differences in GM.

METHODS

Forty-seven PE patients, forty NPHCs, and twenty-one PHCs were included. Brain OEF values were computed from quantitative susceptibility mapping (QSM) plus quantitative blood oxygen level-dependent magnitude (QSM + qBOLD = QQ)-based mapping. One-way ANOVA was used to compare mean OEF values in the three groups. The area under the curve of the mean OEF value in each region of interest was estimated using a receiver operating characteristic curve analysis.

RESULTS

We found that the mean OEF values in the thalamus, putamen, caudate nucleus, pallidum, and substantia nigra were significantly different in these three groups (F = 5.867, p = 0.004; F = 5.142, p = 0007; F = 6.158, p = 0.003; F = 6.319, p = 0.003; F = 5.491, p = 0.005). The mean OEF values for these 5 regions were higher in PE patients than in NPHCs and in PHCs (p < 0.05). The AUC of these ROIs ranged from 0.673 to 0.692 (p < 0.01) and cutoff values varied from 35.1 to 36.6%, indicating that the OEF values could discriminate patients with and without PE. Stepwise multivariate analysis revealed that the OEF values correlated with hematocrit in pregnant women (r = 0.353, p = 0.003).

CONCLUSION

OEF values in the brains of pregnant women can be measured in clinical practice using QQ-based OEF mapping for noninvasive assessment of hypertensive disorders.

KEY POINTS

• Pre-eclampsia is a hypertensive disorder associated with abnormalities in brain oxygen extraction. • Oxygen extraction fraction (OEF) is an indicator of brain tissue viability and function. QQ-based mapping of OEF is a new MRI technique that can noninvasively quantify brain oxygen metabolism. • OEF values in the brains of pregnant women can be measured for noninvasive assessment of hypertensive disorders in clinical practice.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Quantitative susceptibility mapping as an imaging biomarker for Alzheimer’s disease: The expectations and limitations

Alzheimer’s disease (AD) is the most common type of dementia and a distressing diagnosis for individuals and caregivers. Researchers and clinical trials have mainly focused on β-amyloid plaques, which are hypothesized to be one of the most important factors for neurodegeneration in AD. Meanwhile, recent clinicopathological and radiological studies have shown closer associations of tau pathology rather than β-amyloid pathology with the onset and progression of Alzheimer’s symptoms. Toward a biological definition of biomarker-based research framework for AD, the 2018 National Institute on Aging–Alzheimer’s Association working group has updated the ATN classification system for stratifying disease status in accordance with relevant pathological biomarker profiles, such as cerebral β-amyloid deposition, hyperphosphorylated tau, and neurodegeneration. In addition, altered iron metabolism has been considered to interact with abnormal proteins related to AD pathology thorough generating oxidative stress, as some prior histochemical and histopathological studies supported this iron-mediated pathomechanism. Quantitative susceptibility mapping (QSM) has recently become more popular as a non-invasive magnetic resonance technique to quantify local tissue susceptibility with high spatial resolution, which is sensitive to the presence of iron. The association of cerebral susceptibility values with other pathological biomarkers for AD has been investigated using various QSM techniques; however, direct evidence of these associations remains elusive. In this review, we first briefly describe the principles of QSM. Second, we focus on a large variety of QSM applications, ranging from common applications, such as cerebral iron deposition, to more recent applications, such as the assessment of impaired myelination, quantification of venous oxygen saturation, and measurement of blood– brain barrier function in clinical settings for AD. Third, we mention the relationships among QSM, established biomarkers, and cognitive performance in AD. Finally, we discuss the role of QSM as an imaging biomarker as well as the expectations and limitations of clinically useful diagnostic and therapeutic implications for AD.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO 2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

A flow-diffusion model of oxygen transport for quantitative mapping of cerebral metabolic rate of oxygen (CMRO2) with single gas calibrated fMRI

One promising approach for mapping CMRO₂ is dual-calibrated functional MRI (dc-fMRI). This method exploits the Fick Principle to combine estimates of CBF from ASL, and OEF derived from BOLD-ASL measurements during arterial O₂ and CO₂ modulations. Multiple gas modulations are required to decouple OEF and deoxyhemoglobin-sensitive blood volume. We propose an alternative single gas calibrated fMRI framework, integrating a model of oxygen transport, that links blood volume and CBF to OEF and creates a mapping between the maximum BOLD signal, CBF and OEF (and CMRO₂). Simulations demonstrated the method's viability within physiological ranges of mitochondrial oxygen pressure, P_mO₂, and mean capillary transit time. A dc-fMRI experiment, performed on 20 healthy subjects using O₂ and CO₂ challenges, was used to validate the approach. The validation conveyed expected estimates of model parameters (e.g., low P_mO₂), with spatially uniform OEF maps (grey matter, GM, OEF spatial standard deviation ≈ 0.13). GM OEF estimates obtained with hypercapnia calibrated fMRI correlated with dc-fMRI (r = 0.65, p = 2·10^-3). For 12 subjects, OEF measured with dc-fMRI and the single gas calibration method were correlated with whole-brain OEF derived from phase measures in the superior sagittal sinus (r = 0.58, p = 0.048; r = 0.64, p = 0.025 respectively). Simplified calibrated fMRI using hypercapnia holds promise for clinical application.

T2 -oximetry-based cerebral venous oxygenation mapping using Fourier-transform-based velocity-selective pulse trains

Purpose

To develop a T₂‐oximetry method for quantitative mapping of cerebral venous oxygenation fraction (Y_v) using Fourier‐transform–based velocity‐selective (FT‐VS) pulse trains.

Methods

The venous isolation preparation was achieved by using an FT‐VS inversion plus a nonselective inversion (NSI) pulse to null the arterial blood signal while minimally affected capillary blood flows out into the venular vasculature during the outflow time (TO), and then applying an Fourier transform based velocity selective saturation (FT‐VSS) pulse to suppress the tissue signal. A multi‐echo readout was employed to obtain venous T₂ (T_2,v) efficiently with the last echo used to detect the residual CSF signal and correct its contamination in the fitting. Here we compared the performance of this FT‐VS–based venous isolation preparations with a traditional velocity‐selective saturation (VSS)–based approach (quantitative imaging of extraction of oxygen and tissue consumption [QUIXOTIC]) with different cutoff velocities for Y_v mapping on 6 healthy volunteers at 3 Tesla.

Results

The FT‐VS–based methods yielded higher venous blood signal and temporal SNR with less CSF contamination than the velocity‐selective saturation–based results. The averaged Y_v values across the whole slice measured in different experiments were close to the global Y_v measured from the individual internal jugular vein.

Conclusion

The feasibility of the FT‐VS–based Y_v estimation was demonstrated on healthy volunteers. The obtained high venous signal as well as the mitigation of CSF contamination led to a good agreement between the T_2,v and Y_v measured in the proposed method with the values in the literature.

Click here for author‐reader discussions

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.

Brain pathological changes during neurodegenerative diseases and their identification methods: How does QSM perform in detecting this process?

The presence of iron is essential for many biological processes in the body. But sometimes, for various reasons, the amount of iron deposition in different areas of the brain increases, which leads to problems related to the nervous system. Quantitative susceptibility mapping (QSM) is one of the newest magnetic resonance imaging (MRI)-based methods for assessing iron accumulation in target areas. This Narrative Review article aims to evaluate the performance of QSM compared to other methods of assessing iron deposition in the clinical field. Based on the results, we introduced related basic definitions, some neurodegenerative diseases, methods of examining iron deposition in these diseases, and their advantages and disadvantages. This article states that the QSM method can be introduced as a new, reliable, and non-invasive technique for clinical evaluations.

Assessing Mild Cognitive Impairment in Parkinson's Disease by Magnetic Resonance Quantitative Susceptibility Mapping Combined Voxel-Wise and Radiomic Analysis

BACKGROUND

The relationship between iron accumulation in the central nervous system and cognitive decline in Parkinson's disease (PD) has not been fully elucidated. This study aimed to explore the value of quantitative susceptibility mapping in assessment of mild cognitive impairment (MCI) in PD.

METHODS

Sixteen PD patients with MCI (PD-MCI), sixteen normal cognition PD patients (PD-NC), and 28 healthy controls (HCs) were included. The differences in the magnetic susceptibility and Radiomic indicators among groups and their correlations with Montreal Cognitive Assessment-Basic (MoCA-B) scores and Unified Parkinson's Disease Rating Scale Part III (UPDRS-III) were analyzed. Receiver operating characteristic curves were used to evaluate the diagnostic performance.

RESULTS

Higher iron deposition was observed in the cortical and subcortical structures of the PD patients compared with HCs, including limbic system, orbitofrontal cortex, cuneus, red nucleus, and substantia nigra. Combined magnetic susceptibility and texture index in hippocampus achieved the best diagnostic performance (area under curves: 0.828) in differentiating PD-MCI from PD-NC. The magnetic susceptibilities of the substantia nigra, red nucleus, putamen, globus pallidus, hippocampus, and thalamus were negatively correlated with the MoCA-B scores (all p < 0.05), and of the putamen and amygdala were positively correlated with the UPDRS-III scores (both p < 0.05).

CONCLUSION

Higher iron deposition was observed in the cortical and subcortical structures of the PD-MCI and PD-NC groups. The susceptibility values of vulnerable brain subregions shown significant correlation with MoCA-B and UPDRS-III. Together with the texture index, magnetic susceptibility values could provide robust performance in distinguishing PD-MCI patients from PD-NC.

Brain oxygen extraction and neural tissue susceptibility are associated with cognitive impairment in older individuals

BACKGROUND AND PURPOSE

We investigated the effects of aging, white matter hyperintensities (WMH), and cognitive impairment on brain iron levels and cerebral oxygen metabolism, known to be altered in Alzheimer's disease (AD), using quantitative susceptibility mapping and MR-based cerebral oxygen extraction fraction (OEF).

METHODS

In 100 individuals over the age of 50 (68/32 cognitively impaired/intact), OEF and neural tissue susceptibility (χ_n ) were computed retrospectively from MRI multi-echo gradient echo data, obtained on a 3 Tesla MRI scanner. The effects of age and WMH on OEF and χ_n were assessed within groups, and OEF and χ_n were assessed between groups, using multivariate regression analyses.

RESULTS

Cognitively impaired subjects were found to have 19% higher OEF and 34% higher χ_n than cognitively intact subjects in the cortical gray matter and several frontal, temporal, and parietal regions (p < .05). Increased WMH burden was significantly associated with decreased OEF in the cognitively impaired, but not in the cognitively intact. Older age had a stronger association with decreased OEF in the cognitively intact group. Both older age and increased WMH burden were significantly associated with increased χ_n in temporoparietal regions in the cognitively impaired.

CONCLUSIONS

Higher brain OEF and χ_n in cognitively impaired older individuals may reflect altered oxygen metabolism and iron in areas with underlying AD pathology. Both age and WMH have associations with OEF and χ_n but are modified by the presence of cognitive impairment.

Brain oxygen extraction fraction mapping in patients with multiple sclerosis

We aimed to demonstrate the feasibility of whole brain oxygen extraction fraction (OEF) mapping for measuring lesion specific and regional OEF abnormalities in multiple sclerosis (MS) patients. In 22 MS patients and 11 healthy controls (HC), OEF and neural tissue susceptibility (χn) maps were computed from MRI multi-echo gradient echo data. In MS patients, 80 chronic active lesions with hyperintense rim on quantitative susceptibility mapping were identified, and the mean OEF and χn within the rim and core were compared using linear mixed-effect model analysis. The rim showed higher OEF and χn than the core: relative to their adjacent normal appearing white matter, OEF contrast = -6.6 ± 7.0% vs. -9.8 ± 7.8% (p < 0.001) and χn contrast = 33.9 ± 20.3 ppb vs. 25.7 ± 20.5 ppb (p = 0.017). Between MS and HC, OEF and χn were compared using a linear regression model in subject-based regions of interest. In the whole brain, compared to HC, MS had lower OEF, 30.4 ± 3.3% vs. 21.4 ± 4.4% (p < 0.001), and higher χn, -23.7 ± 7.0 ppb vs. -11.3 ± 7.7 ppb (p = 0.018). Our feasibility study suggests that OEF may serve as a useful quantitative marker of tissue oxygen utilization in MS.

Whole-brain 3D mapping of oxygen metabolism using constrained quantitative BOLD

Quantitative BOLD (qBOLD) MRI permits noninvasive evaluation of hemodynamic and metabolic states of the brain by quantifying parametric maps of deoxygenated blood volume (DBV) and hemoglobin oxygen saturation level of venous blood (Y_v), and along with a measurement of cerebral blood flow (CBF), the cerebral metabolic rate of oxygen (CMRO₂). The method, thus should have potential to provide important information on many neurological disorders as well as normal cerebral physiology. One major challenge in qBOLD is to separate de-oxyhemoglobin’s contribution to R₂′ from other sources modulating the voxel signal, for instance, R₂, R₂′ from non-heme iron (R′_2,nh), and macroscopic magnetic field variations. Further, even with successful separation of the several confounders, it is still challenging to extract DBV and Y_v from the heme-originated R₂′ because of limited sensitivity of the qBOLD model. These issues, which have not been fully addressed in currently practiced qBOLD methods, have so far precluded 3D whole-brain implementation of qBOLD. Thus, the purpose of this work was to develop a new 3D MRI oximetry technique that enables robust qBOLD parameter mapping across the entire brain. To achieve this goal, we employed a rapid, R₂′-sensitive, steady-state 3D pulse sequence (termed ‘AUSFIDE’) for data acquisition, and implemented a prior-constrained qBOLD processing pipeline that exploits a plurality of preliminary parameters obtained via AUSFIDE, along with additionally measured cerebral venous blood volume. Numerical simulations and in vivo studies at 3 T were performed to evaluate the performance of the proposed, constrained qBOLD mapping in comparison to the parent qBOLD method. Measured parameters (Y_v, DBV, R′_2,nh, nonblood magnetic susceptibility) in ten healthy subjects demonstrate the expected contrast across brain territories, while yielding group-averages of 64.0 ± 2.3 % and 62.2 ± 3.1 % for Y_v and 2.8 ± 0.5 % and 1.8 ± 0.4 % for DBV in cortical gray and white matter, respectively. Given the Y_v measurements, additionally quantified CBF in seven of the ten study subjects enabled whole-brain 3D CMRO₂ mapping, yielding group averages of 134.2 ± 21.1 and 79.4 ± 12.6 µmol/100 g/min for cortical gray and white matter, in good agreement with literature values. The results suggest feasibility of the proposed method as a practical and reliable means for measuring neurometabolic parameters over an extended brain coverage.

Optimal timing for measuring cerebral blood flow after acetazolamide administration to detect preexisting cerebral hemodynamics and metabolism in patients with bilateral major cerebral artery steno-occlusive diseases: 15O positron emission tomography studies

The present study determined the optimal timing of scanning for measurement of cerebral blood flow (CBF) after acetazolamide (ACZ) administration for detection of preexisting cerebral hemodynamics and metabolism in bilateral major cerebral artery steno-occlusive diseases. Thirty three patients underwent ¹⁵O gas positron emission tomography (PET) and each parameter was obtained in the bilateral middle cerebral artery (MCA) territories. CBF was also obtained using H₂ ¹⁵O PET scanning performed at baseline and at 5, 15, and 30 min after ACZ administration. Relative CBF at each time point after ACZ administration to baseline CBF was calculated. For MCA territories with normal cerebral blood volume (CBV) and cerebral metabolic rate of oxygen (CMRO₂), CBF continued increasing until 15 min after ACZ administration. For MCA territories with abnormally increased CBV, CBF decreased 5 min after ACZ administration. After that, CBF continued increasing until 30 min after ACZ administration. For MCA territories with abnormally decreased CMRO₂, CBF did not change 5 min after ACZ administration. Ten min later, CBF increased. The accuracy to detect abnormally increased CBV was significantly greater for relative CBF₅ than for relative CBF₁₅. The accuracy to detect abnormally decreased CMRO₂ was significantly greater for relative CBF₅ or CBF₁₅ than for relative CBF₃₀. For detecting abnormally increased oxygen extraction fraction, the accuracy did not differ among each relative CBF. These findings suggested that CBF measurement at 5 min after ACZ administration is the optimal timing for detection of preexisting cerebral hemodynamics and metabolism in bilateral major cerebral artery steno-occlusive diseases.

QQ-NET - using deep learning to solve quantitative susceptibility mapping and quantitative blood oxygen level dependent magnitude (QSM+qBOLD or QQ) based oxygen extraction fraction (OEF) mapping

PURPOSE

To improve accuracy and speed of quantitative susceptibility mapping plus quantitative blood oxygen level-dependent magnitude (QSM+qBOLD or QQ) -based oxygen extraction fraction (OEF) mapping using a deep neural network (QQ-NET).

METHODS

The 3D multi-echo gradient echo images were acquired in 34 ischemic stroke patients and 4 healthy subjects. Arterial spin labeling and diffusion weighted imaging (DWI) were also performed in the patients. NET was developed to solve the QQ model inversion problem based on Unet. QQ-based OEF maps were reconstructed with previously introduced temporal clustering, tissue composition, and total variation (CCTV) and NET. The results were compared in simulation, ischemic stroke patients, and healthy subjects using a two-sample Kolmogorov-Smirnov test.

RESULTS

In the simulation, QQ-NET provided more accurate and precise OEF maps than QQ-CCTV with 150 times faster reconstruction speed. In the subacute stroke patients, OEF from QQ-NET had greater contrast-to-noise ratio (CNR) between DWI-defined lesions and their unaffected contralateral normal tissue than with QQ-CCTV: 1.9 ± 1.3 vs 6.6 ± 10.7 (p = 0.03). In healthy subjects, both QQ-CCTV and QQ-NET provided uniform OEF maps.

CONCLUSION

QQ-NET improves the accuracy of QQ-based OEF with faster reconstruction.

Impaired cerebral vascular and metabolic responses to parametric N-back tasks in subjective cognitive decline

Previous studies reported abnormally increased and/or decreased blood oxygen level-dependent (BOLD) activations during functional tasks in subjective cognitive decline (SCD). The neurophysiological basis underlying these functional aberrations remains debated. This study aims to investigate vascular and metabolic responses and their dependence on cognitive processing loads during functional tasks in SCD. Twenty-one SCD and 18 control subjects performed parametric N-back working-memory tasks during MRI scans. Task-evoked percentage changes (denoted as δ) in cerebral blood volume (δCBV), cerebral blood flow (δCBF), BOLD signal (δBOLD) and cerebral metabolic rate of oxygen (δCMRO₂) were evaluated. In the frontal lobe, trends of decreased δCBV, δCBF and δCMRO₂ and increased δBOLD were observed in SCD compared with control subjects under lower loads, and these trends increased to significant differences under the 3-back load. δCBF was significantly correlated with δCMRO₂ in controls, but not in SCD subjects. As N-back loads increased, the differences between SCD and control subjects in δCBF and δCMRO₂ tended to enlarge. In the parietal lobe, no significant between-group difference was observed. Our findings suggested that impaired vascular and metabolic responses to functional tasks occurred in the frontal lobe of SCD, which contributed to unusual BOLD hyperactivation and was modulated by cognitive processing loads.

Hyperpolarized MRI, functional MRI, MR spectroscopy and CEST to provide metabolic information in vivo

Access to metabolic information in vivo using magnetic resonance (MR) technologies has generally been the niche of MR spectroscopy (MRS) and spectroscopic imaging (MRSI). Metabolic fluxes can be studied using the infusion of substrates labeled with magnetic isotopes, with the use of hyperpolarization especially powerful. Unfortunately, these promising methods are not yet accepted clinically, where fast, simple, and reliable measurement and diagnosis are key. Recent advances in functional MRI and chemical exchange saturation transfer (CEST) MRI allow the use of water imaging to study oxygen metabolism and tissue metabolite levels. These, together with the use of novel data analysis approaches such as machine learning for all of these metabolic MR approaches, are increasing the likelihood of their clinical translation.

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 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.

Multiecho complex total field inversion method (mcTFI) for improved signal modeling in quantitative susceptibility mapping

PURPOSE

Typical quantitative susceptibility mapping (QSM) reconstruction steps consist of first estimating the magnetization field from the gradient-echo images, and then reconstructing the susceptibility map from the estimated field. The errors from the field-estimation steps may propagate into the final QSM map, and the noise in the estimated field map may no longer be zero-mean Gaussian noise, thus, causing streaking artifacts in the resulting QSM. A multiecho complex total field inversion (mcTFI) method was developed to compute the susceptibility map directly from the multiecho gradient echo images using an improved signal model that retains the Gaussian noise property in the complex domain. It showed improvements in QSM reconstruction over the conventional field-to-source inversion.

METHODS

The proposed mcTFI method was compared with the nonlinear total field inversion (nTFI) method in a numerical brain with hemorrhage and calcification, the numerical brains provided by the QSM Challenge 2.0, 18 brains with intracerebral hemorrhage scanned at 3T, and 6 healthy brains scanned at 7T.

RESULTS

Compared with nTFI, the proposed mcTFI showed more accurate QSM reconstruction around the lesions in the numerical simulations. The mcTFI reconstructed QSM also showed the best image quality with the least artifacts in the brains with intracerebral hemorrhage scanned at 3T and healthy brains scanned at 7T.

CONCLUSION

The proposed multiecho complex total field inversion improved QSM reconstruction over traditional field-to-source inversion through better signal modeling.

Non-Invasive measurement of the cerebral metabolic rate of oxygen using MRI in rodents

Malfunctions of oxygen metabolism are suspected to play a key role in a number of neurological and psychiatric disorders, but this hypothesis cannot be properly investigated without an in-vivo non-invasive measurement of brain oxygen consumption. We present a new way to measure the Cerebral Metabolic Rate of Oxygen (CMRO2) by combining two existing magnetic resonance imaging techniques, namely arterial spin-labelling and oxygen extraction fraction mapping. This method was validated by imaging rats under different anaesthetic regimes and was strongly correlated to glucose consumption measured by autoradiography.

Prediction of an oxygen extraction fraction map by convolutional neural network: validation of input data among MR and PET images

Purpose

Oxygen extraction fraction (OEF) is a biomarker for the viability of brain tissue in ischemic stroke. However, acquisition of the OEF map using positron emission tomography (PET) with oxygen-15 gas is uncomfortable for patients because of the long fixation time, invasive arterial sampling, and radiation exposure. We aimed to predict the OEF map from magnetic resonance (MR) and PET images using a deep convolutional neural network (CNN) and to demonstrate which PET and MR images are optimal as inputs for the prediction of OEF maps.

Methods

Cerebral blood flow at rest (CBF) and during stress (sCBF), cerebral blood volume (CBV) maps acquired from oxygen-15 PET, and routine MR images (T1-, T2-, and T2*-weighted images) for 113 patients with steno-occlusive disease were learned with U-Net. MR and PET images acquired from the other 25 patients were used as test data. We compared the predicted OEF maps and intraclass correlation (ICC) with the real OEF values among combinations of MRI, CBF, CBV, and sCBF.

Results

Among the combinations of input images, OEF maps predicted by the model learned with MRI, CBF, CBV, and sCBF maps were the most similar to the real OEF maps (ICC: 0.597 ± 0.082). However, the contrast of predicted OEF maps was lower than that of real OEF maps.

Conclusion

These results suggest that the deep CNN learned useful features from CBF, sCBF, CBV, and MR images and predict qualitatively realistic OEF maps. These findings suggest that the deep CNN model can shorten the fixation time for ¹⁵O PET by skipping ¹⁵O₂ scans. Further training with a larger data set is required to predict accurate OEF maps quantitatively.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11548-021-02356-7.

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.

Quantitative susceptibility-weighted imaging may be an accurate method for determining stroke hypoperfusion and hypoxia of penumbra

OBJECTIVES

To quantitatively evaluate the volume of the ischemic penumbra using susceptibility-weighted imaging and mapping (SWIM) of asymmetrical prominent cortical veins (APCVs) in patients with acute ischemic stroke.

METHODS

Eighty-five eligible patients with acute ischemic stroke on admission within 12 h from symptom onset were studied. The APCVs on SWIM were quantitatively (SWI-volume) and semi-quantitatively (SWI-Alberta Stroke Program Early CT Score, SWI-ASPECTS) evaluated to calculate mismatch. To assess the diagnostic efficacy of APCVs on SWIM, comparative analyses were performed between SWI_volume-DWI mismatch and SWI_ASPECTS-DWI mismatch, using PWI-DWI mismatch as a reference. Correlations were calculated between the mismatches, as well as between SWI-volume and time-to-maximum (Tmax) > 6 s volume. Additionally, each of these mismatches was correlated with the National Institute of Health Stroke Scale (NIHSS).

RESULTS

The sensitivity, negative predictive value, and accuracy of SWI_volume-DWI mismatch were demonstrably higher than SWI_ASPECTS-DWI mismatch (100% vs. 53.7%, 100% vs. 9.5%, 97.7% vs. 54.5%, respectively). A significant positive correlation was found between SWI_volume-DWI and PWI-DWI mismatch (r = 0.691, p < 0.01), as well as between SWI-volume and Tmax > 6 s volume (r = 0.786, p < 0.001). A significant negative correlation was found between SWI_volume-DWI mismatch and NIHSS (r = - 0.360, p = 0.022), as well as between SWI_ASPECTS-DWI mismatch and NIHSS (r = - 0.499, p = 0.001).

CONCLUSIONS

SWI_volume-DWI mismatch had higher diagnostic efficacy than SWI_ASPECTS-DWI mismatch in defining the ischemic penumbra and showed good consistency with PWI-DWI mismatch in acute ischemic stroke. Quantitation of APCVs using SWIM provided an accurate method for determining hypoperfusion and provided a reliable method to reflect the hypoxia of penumbra.

KEY POINTS

• SWI_volume-DWI mismatch has higher diagnostic efficacy than SWI_ASPECTS-DWI mismatch in defining the ischemic penumbra. • SWI_volume-DWI mismatch shows good consistency with PWI-DWI mismatch in managing penumbra in acute ischemic stroke. • Quantitation of APCV volume using SWIM provided an accurate method for determining the hypoperfusion area and provided a reliable method to reflect the hypoxia of penumbra.

Quantitative transport mapping (QTM) of the kidney with an approximate microvascular network

Purpose

Proof‐of‐concept study of mapping renal blood flow vector field according to the inverse solution to a mass transport model of time resolved tracer‐labeled MRI data.

Theory and Methods

To determine tissue perfusion according to the underlying physics of spatiotemporal tracer concentration variation, the mass transport equation is integrated over a voxel with an approximate microvascular network for fitting time‐resolved tracer imaging data. The inverse solution to the voxelized transport equation provides the blood flow vector field, which is referred to as quantitative transport mapping (QTM). A numerical microvascular network modeling the kidney with computational fluid dynamics reference was used to verify the accuracy of QTM and the current Kety’s method that uses a global arterial input function. Multiple post‐label delay arterial spin labeling (ASL) of the kidney on seven subjects was used to assess QTM in vivo feasibility.

Results

Against the ground truth in the numerical model, the error in flow estimated by QTM (18.6%) was smaller than that in Kety’s method (45.7%, 2.5‐fold reduction). The in vivo kidney perfusion quantification by QTM (cortex: 443 ± 58 mL/100 g/min and medulla: 190 ± 90 mL/100 g/min) was in the range of that by Kety’s method (482 ± 51 mL/100 g/min in the cortex and 242 ± 73 mL/100 g/min in the medulla), and QTM provided better flow homogeneity in the cortex region.

Conclusions

QTM flow velocity mapping is feasible from multi‐delay ASL MRI data based on inverting the transport equation. In a numerical simulation, QTM with deconvolution in space and time provided more accurate perfusion quantification than Kety’s method with deconvolution in time only.

Measurement of CMRO2 and its relationship with CBF in hypoxia with an extended calibrated BOLD method

Cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO₂) are physiological parameters that not only reflect brain health and disease but also jointly contribute to blood oxygen level-dependent (BOLD) signals. Nevertheless, unsolved issues remain concerning the CBF-CMRO₂ relationship in the working brain under various oxygen conditions. In particular, the CMRO₂ responses to functional tasks in hypoxia are less studied. We extended the calibrated BOLD model to incorporate CMRO₂ measurements in hypoxia. The extended model, which was cross-validated with a multicompartment BOLD model, considers the influences of the reduced arterial saturation level and increased baseline cerebral blood volume (CBV) and deoxyhemoglobin concentration on the changes of BOLD signals in hypoxia. By implementing a pulse sequence to simultaneously acquire the CBV-, CBF- and BOLD-weighted signals, we investigated the effects of mild hypoxia on the CBF and CMRO₂ responses to graded visual stimuli. Compared with normoxia, mild hypoxia caused significant alterations in both the amplitude and the trend of the CMRO₂ responses but did not impact the corresponding CBF responses. Our observations suggested that the flow-metabolism coupling strategies in the brain during mild hypoxia were different from those during normoxia.

Venous cerebral blood volume mapping in the whole brain using venous-spin-labeled 3D turbo spin echo

PURPOSE

Venous cerebral blood volume (CBVv ) is a major contributor to BOLD contrast, and therefore is an important parameter for understanding the underlying mechanism. Here, we propose a velocity-selective venous spin labeling (VS-VSL)-prepared 3D turbo spin echo pulse sequence for whole-brain baseline CBVv mapping.

METHODS

Unlike previous CBVv measurement techniques that exploit the interrelationship between BOLD signals and CBVv , in the proposed VS-VSL technique both arterial blood and cerebrospinal fluid (CSF) signals were suppressed before the VS pulse train for exclusive labeling of venous blood, while a single-slab 3D turbo spin echo readout was used because of its relative immunity to magnetic field variations. Furthermore, two approximations were made to the VS-VSL signal model for simplified derivation of CBVv . In vivo studies were performed at 3T field strength in 8 healthy subjects. The performance of the proposed VS-VSL method in baseline CBVv estimation was first evaluated in comparison to the existing, hyperoxia-based method. Then, data were also acquired using VS-VSL under hypercapnic and hyperoxic gas breathing challenges for further validation of the technique.

RESULTS

The proposed technique yielded physiologically plausible baseline CBVv values, and when compared with the hyperoxia-based method, showed no statistical difference. Furthermore, data acquired using VS-VSL yielded average CBVv of 2.89%/1.78%, 3.71%/2.29%, and 2.88%/1.76% for baseline, hypercapnia, and hyperoxia, respectively, in gray/white matter regions. As expected, hyperoxia had negligible effect (P > .8), whereas hypercapnia demonstrated vasodilation (P << .01).

CONCLUSION

Upon further validation of the quantification model, the method is expected to have merit for 3D CBVv measurements across the entire brain.

Initial Experience of Challenge-Free MRI-Based Oxygen Extraction Fraction Mapping of Ischemic Stroke at Various Stages: Comparison With Perfusion and Diffusion Mapping

MRI-based oxygen extraction fraction imaging has a great potential benefit in the selection of clinical strategies for ischemic stroke patients. This study aimed to evaluate the performance of a challenge-free oxygen extraction fraction (OEF) mapping in a cohort of acute and subacute ischemic stroke patients. Consecutive ischemic stroke patients (a total of 30 with 5 in the acute stage, 19 in the early subacute stage, and 6 in the late subacute stage) were recruited. All subjects underwent MRI including multi-echo gradient echo (mGRE), diffusion weighted imaging (DWI), and 3D-arterial spin labeling (ASL). OEF maps were generated from mGRE phase + magnitude data, which were processed using quantitative susceptibility mapping (QSM) + quantitative blood oxygen level-dependent (qBOLD) imaging with cluster analysis of time evolution. Cerebral blood flow (CBF) and apparent diffusion coefficient (ADC) maps were reconstructed from 3D-ASL and DWI, respectively. Further, cerebral metabolic rate of oxygen (CMRO2) was calculated as the product of CBF and OEF. OEF, CMRO2, CBF, and ADC values in the ischemic cores (absolute values) and their contrasts to the contralateral regions (relative values) were evaluated. One-way analysis of variance (ANOVA) was used to compare OEF, CMRO2, CBF, and ADC values and their relative values among different stroke stages. The OEF value of infarct core showed a trend of decrease from acute, to early subacute, and to late subacute stages of ischemic stroke. Significant differences among the three stroke stages were only observed in the absolute OEF (F = 6.046, p = 0.005) and relative OEF (F = 5.699, p = 0.009) values of the ischemic core, but not in other measurements (absolute and relative CMRO2, CBF, ADC values, all values of p > 0.05). In conclusion, the challenge-free QSM + qBOLD-generated OEF mapping can be performed on stroke patients. It can provide more information on tissue viability that was not available with CBF and ADC and, thus, may help to better manage ischemic stroke patients.

Simultaneous QSM and metabolic imaging of the brain using SPICE: Further improvements in data acquisition and processing

PURPOSE

To achieve high-resolution mapping of brain tissue susceptibility in simultaneous QSM and metabolic imaging.

METHODS

Simultaneous QSM and metabolic imaging was first achieved using SPICE (spectroscopic imaging by exploiting spatiospectral correlation), but the QSM maps thus obtained were at relatively low-resolution (2.0 × 3.0 × 3.0 mm³ ). We overcome this limitation using an improved SPICE data acquisition method with the following novel features: 1) sampling (k, t)-space in dual densities, 2) sampling central k-space fully to achieve nominal spatial resolution of 3.0 × 3.0 × 3.0 mm³ for metabolic imaging, and 3) sampling outer k-space sparsely to achieve spatial resolution of 1.0 × 1.0 × 1.9 mm³ for QSM. To keep the scan time short, we acquired spatiospectral encodings in echo-planar spectroscopic imaging trajectories in central k-space but in CAIPIRINHA (controlled aliasing in parallel imaging results in higher acceleration) trajectories in outer k-space using blipped phase encodings. For data processing and image reconstruction, a union-of-subspaces model was used, effectively incorporating sensitivity encoding, spatial priors, and spectral priors of individual molecules.

RESULTS

In vivo experiments were carried out to evaluate the feasibility and potential of the proposed method. In a 6-min scan, QSM maps at 1.0 × 1.0 × 1.9 mm³ resolution and metabolic maps at 3.0 × 3.0 × 3.0 mm³ nominal resolution were obtained simultaneously. Compared with the original method, the QSM maps obtained using the new method reveal fine-scale brain structures more clearly.

CONCLUSION

We demonstrated the feasibility of achieving high-resolution QSM simultaneously with metabolic imaging using a modified SPICE acquisition method. The improved capability of SPICE may further enhance its practical utility in brain mapping.

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.

Quantitative Magnetic Resonance Imaging of Cortical and Trabecular Bone

Bone is a composite material consisting of mineral, organic matrix, and water. Water in bone can be categorized as bound water (BW), which is bound to bone mineral and organic matrix, or as pore water (PW), which resides in Haversian canals as well as in lacunae and canaliculi. Bone is generally classified into two types: cortical bone and trabecular bone. Cortical bone is much denser than trabecular bone that is surrounded by marrow and fat. Magnetic resonance (MR) imaging has been increasingly used for noninvasive assessment of both cortical bone and trabecular bone. Bone typically appears as a signal void with conventional MR sequences because of its short T2*. Ultrashort echo time (UTE) sequences with echo times 100 to 1,000 times shorter than those of conventional sequences allow direct imaging of BW and PW in bone. This article summarizes several quantitative MR techniques recently developed for bone evaluation. Specifically, we discuss the use of UTE and adiabatic inversion recovery prepared UTE sequences to quantify BW and PW, UTE magnetization transfer sequences to quantify collagen backbone protons, UTE quantitative susceptibility mapping sequences to assess bone mineral, and conventional sequences for high-resolution imaging of PW as well as the evaluation of trabecular bone architecture.