Venous oxygenation mapping using velocity-selective excitation and arterial nulling

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.

Direct angiographic comparison of different velocity-selective saturation, inversion, and DANTE labeling modules on cerebral arteries

PURPOSE

This study evaluated the velocity-selective (VS) MRA with different VS labeling modules, including double refocused hyperbolic tangent, eight-segment B1-insensitive rotation, delay alternating with nutation for tailored excitation, Fourier transform-based VS saturation, and Fourier transform-based inversion.

METHODS

These five VS labeling modules were evaluated first through Bloch simulations, and then using VSMRA directly on various cerebral arteries of healthy subjects. The relative signal ratios from arterial ROIs and surrounding tissues as well as relative arteria-tissue contrast ratios of different methods were compared.

RESULTS

Double refocused hyperbolic tangent and eight-segment B1-insensitive rotation showed very similar labeling effects. Delay alternating with nutation for tailored excitation yielded high arterial signal but with residual tissue signal due to the spatial banding effect. Fourier transform-based VS saturation with half the time of other techniques serves as an efficient nonsubtractive VSMRA method, but the remaining tissue signal still obscured some small distal arteries that were delineated by other subtraction-based VSMRA, allowing more complete cancelation of static tissue. Fourier transform-based inversion produced the highest arterial signal in VSMRA with minimal tissue background.

CONCLUSION

This is the first study that angiographically compared five different VS labeling modules. Their labeling characteristics on arteries and tissue and implications for VSMRA and VS arterial spin labeling are discussed.

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.

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.

Quantification of regional CMRO2 in human brain using dynamic 17O-MRI at 3T

OBJECTIVE

To investigate the feasibility of cerebral metabolic rate of oxygen consumption (CMRO2) measurements with MRI at 3 Tesla in different brain regions.

METHODS

CMRO2 represents a key indicator of the physiological state of brain tissue. Dynamic 17O-MRI with inhalation of isotopically enriched 17O gas has been used to quantify global CMRO2 in brain white (WM) and gray matter (GM). However, global CMRO2 can only reflect the overall oxygen metabolism of the brain and cannot provide enough information on local tissue oxygen metabolism. To investigate the feasibility of determination of regional CMRO2 at a clinical 3 T MRI system, CMRO2 values in frontal, parietal and occipital WM and GM were determined in 5 healthy volunteers and compared to evaluate the regional differences of oxygen metabolism in WM and GM. Additionally, regional CMRO2 values were determined in deep brain structures including thalamus, dorsal striatum, caudate nucleus and insula cortex and in the cerebella, and compared with literature values from 15O-PET studies.

RESULTS

In cortical GM the determined CMRO2 values were in good agreement with the literature, whereas values in WM were about 32-48% higher than literature values. Regional analysis revealed a significantly higher CMRO2 in the occipital GM compared to the frontal and parietal GM. By contrast, no significant difference of CMRO2 was observed across the WM. In addition, CMRO2 in deep brain structures was lower compared to literature values and in the cerebella a good hemispheric symmetry of the tissue oxygen metabolism was found.

CONCLUSION

Dynamic 17O-MRI enables direct, non-invasive determination of regional CMRO2 in brain structures in healthy volunteers at 3T.

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.

Robust dual-module velocity-selective arterial spin labeling (dm-VSASL) with velocity-selective saturation and inversion

PURPOSE

Compared to conventional arterial spin labeling (ASL) methods, velocity-selective ASL (VSASL) is more sensitive to artifacts from eddy currents, diffusion attenuation, and motion. Background suppression is typically suboptimal in VSASL, especially of CSF. As a result, the temporal SNR and quantification accuracy of VSASL are compromised, hindering its application despite its advantage of being delay-insensitive.

METHODS

A novel dual-module VSASL (dm-VSASL) strategy is developed to improve the SNR efficiency and the temporal SNR with a more balanced gradient configuration in the label/control image acquisition. This strategy applies for both VS saturation (VSS) and VS inversion (VSI) labeling. The dm-VSASL schemes were compared with single-module labeling and a previously developed multi-module schemes for the SNR performance, background suppression efficacy, and sensitivity to artifacts in simulation and in vivo experiments, using pulsed ASL as the reference.

RESULTS

Dm-VSASL enabled more robust labeling and efficient backgroud suppre across brain tissues, especially of CSF, resulting in significantly reduced artifacts and improved temporal SNR. Compared to single-module labeling, dm-VSASL significantly improved the temporal SNR in gray (by 90.8% and 94.9% for dm-VSS and dm-VSI, respectively; P < 0.001) and white (by 41.5% and 55.1% for dm-VSS and dm-VSI, respectively; P < 0.002) matter. Dm-VSI also improved the SNR of VSI by 5.4% (P = 0.018).

CONCLUSION

Dm-VSASL can significantly improve the robustness of VS labeling, reduce artifacts, and allow efficient background suppression. When implemented with VSI, it provides the highest SNR efficiency among VSASL methods. Dm-VSASL is a powerful ASL method for robust, accurate, and delay-insensitive perfusion mapping.

Velocity-selective excitation: Principles and applications

Velocity-selective (VS) excitation is a relatively new type of excitation that can be useful for generating image contrast based on spin's motion. This review aims to explain the principles of VS excitation and their utilization for clinical applications. We first review the generalized excitation k-space formalism, which reveals a Fourier relationship between sequence parameters and excitation profiles for spins with arbitrary spatial location, off-resonance, and velocity. Based on the k-space framework, we analyze practical VS excitation pulse sequences that yield sinusoidal or sinc-shaped velocity profiles. Then we demonstrate how these two types of VS excitation can be used as magnetization preparation for clinical applications, including saturation- or inversion-based arterial spin labeling and black- or bright-blood angiography. We also discuss practical considerations and issues for each application, including the determination of design parameters and the effects of MR system errors, such as magnetic field offsets and eddy currents.

Recent Technical Developments in ASL: A Review of the State of the Art

This review article provides an overview of a range of recent technical developments in advanced arterial spin labeling (ASL) methods that have been developed or adopted by the community since the publication of a previous ASL consensus paper by Alsop et al. It is part of a series of review/recommendation papers from the International Society for Magnetic Resonance in Medicine Perfusion Study Group. Here, we focus on advancements in readouts and trajectories, image reconstruction, noise reduction, partial volume correction, quantification of nonperfusion parameters, fMRI, fingerprinting, vessel selective ASL, angiography, deep learning, and ultrahigh field ASL. We aim to provide a high level of understanding of these new approaches and some guidance for their implementation, with the goal of facilitating the adoption of such advances by research groups and by MRI vendors. Topics outside the scope of this article that are reviewed at length in separate articles include velocity selective ASL, multiple-timepoint ASL, body ASL, and clinical ASL recommendations.

Velocity-selective arterial spin labeling perfusion MRI: A review of the state of the art and recommendations for clinical implementation

This review article provides an overview of the current status of velocity-selective arterial spin labeling (VSASL) perfusion MRI and is part of a wider effort arising from the International Society for Magnetic Resonance in Medicine (ISMRM) Perfusion Study Group. Since publication of the 2015 consensus paper on arterial spin labeling (ASL) for cerebral perfusion imaging, important advancements have been made in the field. The ASL community has, therefore, decided to provide an extended perspective on various aspects of technical development and application. Because VSASL has the potential to become a principal ASL method because of its unique advantages over traditional approaches, an in-depth discussion was warranted. VSASL labels blood based on its velocity and creates a magnetic bolus immediately proximal to the microvasculature within the imaging volume. VSASL is, therefore, insensitive to transit delay effects, in contrast to spatially selective pulsed and (pseudo-) continuous ASL approaches. Recent technical developments have improved the robustness and the labeling efficiency of VSASL, making it a potentially more favorable ASL approach in a wide range of applications where transit delay effects are of concern. In this review article, we (1) describe the concepts and theoretical basis of VSASL; (2) describe different variants of VSASL and their implementation; (3) provide recommended parameters and practices for clinical adoption; (4) describe challenges in developing and implementing VSASL; and (5) describe its current applications. As VSASL continues to undergo rapid development, the focus of this review is to summarize the fundamental concepts of VSASL, describe existing VSASL techniques and applications, and provide recommendations to help the clinical community adopt VSASL.

Anatomy imaging and hemodynamics research on the cerebral vein and venous sinus among individuals without cranial sinus and jugular vein diseases

In recent years, imaging technology has allowed the visualization of intracranial and extracranial vascular systems. However, compared with the cerebral arterial system, the relative lack of image information, individual differences in the anatomy of the cerebral veins and venous sinuses, and several unique structures often cause neurologists and radiologists to miss or over-diagnose. This increases the difficulty of the clinical diagnosis and treatment of cerebral venous system diseases. This review focuses on applying different imaging methods to the normal anatomical morphology of the cerebral venous system and special structural and physiological parameters, such as hemodynamics, in people without cranial sinus and jugular vein diseases and explores its clinical significance. We hope this study will reinforce the importance of studying the cerebral venous system anatomy and imaging data and will help diagnose and treat systemic diseases.

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.

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.

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.

Magnetic resonance angiography and perfusion mapping by arterial spin labeling using Fourier transform-based velocity-selective pulse trains: Examination on a commercial perfusion phantom

PURPOSE

Benchmarking of flow and perfusion MR techniques on standardized phantoms can facilitate the use of advanced angiography and perfusion-mapping techniques across multiple sites, field strength, and vendors. Here, MRA and perfusion mapping by arterial spin labeling (ASL) using Fourier transform (FT)-based velocity-selective saturation and inversion pulse trains were evaluated on a commercial perfusion phantom.

METHODS

The FT velocity-selective saturation-based MRA and FT velocity-selective inversion-based ASL perfusion imaging were compared with time-of-flight and pseudo-continuous ASL at 3 T on the perfusion phantom at two controlled flow rates, 175 mL/min and 350 mL/min. Velocity-selective MRA (VSMRA) and velocity-selective ASL (VSASL) were each performed with three velocity-encoding directions: foot-head, left-right, and oblique 45°. The contrast-to-noise ratio for MRA scans and perfusion-weighted signal, as well as labeling efficiency for ASL methods, were quantified.

RESULTS

On this phantom with feeding tubes having only vertical and transverse flow directions, VSMRA and VSASL exhibited the dependence of velocity-encoding directions. The foot-head-encoded VSMRA and VSASL generated similar signal contrasts as time of flight and pseudo-continuous ASL for the two flow rates, respectively. The oblique 45°-encoded VSMRA yielded more uniform contrast-to-noise ratio across slices than foot-head and left-right-encoded VSMRA scans. The oblique 45°-encoded VSASL elevated labeling efficiency from 0.22-0.68 to 0.82-0.90 through more uniform labeling of the entire feeding tubes.

CONCLUSION

Both FT velocity-selective saturation-based VSMRA and FT velocity-selective inversion-based VSASL were characterized on a commercial perfusion phantom. Careful selection of velocity-encoding directions along the major vessels is recommended for their applications in various organs.

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.

Three-dimensional whole-brain mapping of cerebral blood volume and venous cerebral blood volume using Fourier transform-based velocity-selective pulse trains

PURPOSE

To develop 3D MRI methods for cerebral blood volume (CBV) and venous cerebral blood volume (vCBV) estimation with whole-brain coverage using Fourier transform-based velocity-selective (FT-VS) pulse trains.

METHODS

For CBV measurement, FT-VS saturation pulse trains were used to suppress static tissue, whereas CSF contamination was corrected voxel-by-voxel using a multi-readout acquisition and a fast CSF T₂ scan. The vCBV mapping was achieved by inserting an arterial-nulling module that included a FT-VS inversion pulse train. Using these methods, CBV and vCBV maps were obtained on 6 healthy volunteers at 3 T.

RESULTS

The mean CBV and vCBV values in gray matter and white matter in different areas of the brain showed high correlation (r = 0.95 and P < .0001). The averaged CBV and vCBV values of the whole brain were 5.4 ± 0.6 mL/100 g and 2.5 ± 0.3 mL/100 g in gray matter, and 2.6 ± 0.5 mL/100 g and 1.5 ± 0.2 mL/100 g in white matter, respectively, comparable to the literature.

CONCLUSION

The feasibility of FT-VS-based CBV and vCBV estimation was demonstrated for 3D acquisition with large spatial coverage.

Validation of T2 -based oxygen extraction fraction measurement with 15 O positron emission tomography

PURPOSE

To evaluate the accuracy of T₂ -based whole-brain oxygen extraction fraction (OEF) estimation by comparing it with gold standard ¹⁵ O-PET measurements.

METHODS

Sixteen healthy adult subjects underwent MRI and ¹⁵ O-PET OEF measurements on the same day. On MRI, whole-brain OEF was quantified by T₂ -relaxation-under-spin-tagging (TRUST) MRI, based on subject-specific hematocrit. The TRUST OEF was compared to the whole-brain averaged OEF produced by ¹⁵ O-PET. Agreement between TRUST and ¹⁵ O-PET whole-brain OEF measurements was examined in terms of intraclass correlation coefficient (ICC) and in absolute OEF values. In a subset of 10 subjects, test-retest reproducibility of whole-brain OEF was also evaluated and compared between the two modalities.

RESULTS

Across the 16 subjects, the mean whole-brain OEF of TRUST and ¹⁵ O-PET were 36.44 ± 4.07% and 36.45 ± 3.65%, respectively, showing no difference between the two modalities (P = .99). TRUST whole-brain OEF strongly correlated with that of ¹⁵ O-PET (N = 16, ICC = 0.90, P = 4 × 10^-7 ). The coefficient-of-variation of TRUST and ¹⁵ O-PET whole-brain OEF measurements were 1.79 ± 0.67% and 2.06 ± 1.55%, respectively, showing no difference between the two modalities (N = 10, P = .64). Further analyses on the effect of hematocrit revealed that correlation between PET OEF and TRUST OEF with assumed hematocrit remained significant (ICC = 0.8, P < 2 × 10^-5 ).

CONCLUSION

Whole-brain OEF measured by TRUST was in excellent agreement with gold standard ¹⁵ O-PET, with highly comparable accuracy and reproducibility. These findings suggest that TRUST MRI can provide accurate quantification of whole-brain OEF noninvasively.

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.

Normal variations in brain oxygen extraction fraction are partly attributed to differences in end-tidal CO2

Cerebral oxygen extraction fraction is an important physiological index of the brain's oxygen consumption and supply and has been suggested to be a potential biomarker for a number of diseases such as stroke, Alzheimer's disease, multiple sclerosis, sickle cell disease, and metabolic disorders. However, in order for oxygen extraction fraction to be a sensitive biomarker for personalized disease diagnosis, inter-subject variations in normal subjects must be minimized or accounted for, which will otherwise obscure its interpretation. Therefore, it is essential to investigate the physiological underpinnings of normal differences in oxygen extraction fraction. This work used two studies, one discovery study and one verification study, to examine the extent to which an individual's end-tidal CO₂ can explain variations in oxygen extraction fraction. It was found that, across normal subjects, oxygen extraction fraction is inversely correlated with end-tidal CO₂. Approximately 50% of the inter-subject variations in oxygen extraction fraction can be attributed to end-tidal CO₂ differences. In addition, oxygen extraction fraction was found to be positively associated with age and systolic blood pressure. By accounting for end-tidal CO₂, age, and systolic blood pressure of the subjects, normal variations in oxygen extraction fraction can be reduced by 73%, which is expected to substantially enhance the utility of oxygen extraction fraction as a disease biomarker.

The potential for gas-free measurements of absolute oxygen metabolism during both baseline and activation states in the human brain

Quantitative functional magnetic resonance imaging methods make it possible to measure cerebral oxygen metabolism (CMRO₂) in the human brain. Current methods require the subject to breathe special gas mixtures (hypercapnia and hyperoxia). We tested a noninvasive suite of methods to measure absolute CMRO₂ in both baseline and dynamic activation states without the use of special gases: arterial spin labeling (ASL) to measure baseline and activation cerebral blood flow (CBF), with concurrent measurement of the blood oxygenation level dependent (BOLD) signal as a dynamic change in tissue R₂*; VSEAN to estimate baseline O₂ extraction fraction (OEF) from a measurement of venous blood R₂, which in combination with the baseline CBF measurement yields an estimate of baseline CMRO₂; and FLAIR-GESSE to measure tissue R₂' to estimate the scaling parameter needed for calculating the change in CMRO₂ in response to a stimulus with the calibrated BOLD method. Here we describe results for a study sample of 17 subjects (8 female, mean age = 25.3 years, range 21-31 years). The primary findings were that OEF values measured with the VSEAN method were in good agreement with previous PET findings, while estimates of the dynamic change in CMRO₂ in response to a visual stimulus were in good agreement between the traditional hypercapnia calibration and calibration based on R₂'. These results support the potential of gas-free methods for quantitative physiological measurements.

Quantification of whole-brain oxygenation extraction fraction and cerebral metabolic rate of oxygen consumption in adults with sickle cell anemia using individual T2 -based oxygenation calibrations

PURPOSE

To evaluate different T₂ -oxygenation calibrations for estimating venous oxygenation in people with sickle cell anemia (SCA).

METHODS

Blood T₂ values were measured at 3 T in the internal jugular veins of 12 healthy volunteers and 11 SCA participants with no history of stroke, recent transfusion, or renal impairment. T₂ -oxygenation relationships of both sickled and normal blood samples were calibrated individually and compared with values generated from published models. After converting venous T₂ values to venous oxygenation, whole-brain oxygen extraction fraction and cerebral metabolic rate of oxygen were calculated.

RESULTS

Sickle blood samples' oxygenation values calculated from our individual calibrations agreed well with measurements using a blood analyzer, whereas previous T₂ calibrations based on normal blood samples showed 13%-19% underestimation. Meanwhile, oxygenation values calculated from previous grouped T₂ calibration for sickle blood agreed well with experimental measurement on averaged values, but showed up to 20% variation for several individual samples. Using individual T₂ calibrations, the whole-brain oxygen extraction fraction and cerebral metabolic rate of oxygen of SCA participants were 0.38 ± 0.08 and 172 ± 42 µmol/min/100 g, respectively, which were comparable to those values measured on healthy volunteers.

CONCLUSION

Our results confirm that sickle blood T₂ values not only depend on the hematocrit and oxygenation values, but also on other hematological factors. The individual T₂ calibrations minimized the effect of heterogeneity of sickle blood between different SCA populations and improved the accuracy of T₂ -based oximetry. The measured oxygen extraction fraction and cerebral metabolic rate of oxygen of this group of SCA participants were found to not differ significantly from those of healthy individuals.

Characterization of MRI techniques to assess neonatal brain oxygenation and blood flow

There is increasing interest in applying physiological MRI in neonates, based on the premise that physiological parameters may provide an early biomarker of neonatal brain health and injury. Two commonly used techniques are oxygen extraction fraction (OEF) measurement using T₂ -relaxation-under-spin-tagging (TRUST) MRI and cerebral blood flow measurement using phase-contrast (PC) quantitative flow MRI, which collectively provide an assessment of the brain's oxygen consumption. However, prior research has only demonstrated proof of principle of these methods in neonates, without characterization or benchmarking of the techniques. This is because available time is limited in neonatal subjects, especially when scans are performed as add-ons to clinical scans (typically less than 5 min). The work presented aims to examine the TRUST and PC MRI sequences systematically in normal neonates, through research-dedicated scan sessions. A series of characterization and optimization studies were conducted in a total of 26 radiographically normal neonates on 3 T systems. Our results show that TRUST MRI at the superior sagittal sinus (SSS) provides an OEF measurement equivalent to that at the internal jugular vein (r = 0.80, n = 10), yet with shorter scan time. Lower resolution provided better precision in the TRUST measurement (p = 0.001, n = 9). Therefore, the preferred OEF measurement is to apply TRUST MRI at the SSS using a spatial resolution of 2.5 mm. For PC MRI, our results showed that non-gated PC MRI yielded blood flow measurements comparable to those from the more time-consuming gated approach in neonates (r = 0.89, n = 7). It was also found that blood flow could be overestimated by 18% when imaging resolution is larger than 0.3 mm (n = 7). Therefore, non-gated PC MRI with a spatial resolution of 0.3 mm is recommended for neonatal applications. In conclusion, this study verifies consistency of neonatal brain oxygenation and flow measurements across acquisition schemes and points to optimal strategies in parameter selection when using these sequences.

Vessel-specific quantification of neonatal cerebral venous oxygenation

PURPOSE

Noninvasive measurement of cerebral venous oxygenation (Y_v ) in neonates is important in the assessment of brain oxygen extraction and consumption, and may be useful in characterizing brain development and neonatal brain diseases. This study aims to develop a rapid method for vessel-specific measurement of Y_v in neonates.

METHODS

We developed a pulse sequence, named accelerated T₂ -relaxation-under-phase-contrast (aTRUPC), which consists of velocity-encoding phase-contrast module to isolate pure blood signal, flow-insensitive T₂ -preparation to quantify blood T₂ , and turbo-field-echo (TFE) scheme for rapid image acquisition, which is critical for neonatal MRI. A series of studies were conducted. First, the pulse sequence was optimized in terms of TFE factor, velocity encoding (VENC), and slice thickness for best sensitivity. Second, to account for the influence of TFE acquisition on T₂ quantification, simulation and experiments were conducted to establish the relationship between TFE-T₂ and standard T₂ . Finally, the complete aTRUPC sequence was applied on a group of healthy neonates and normative Y_v values were determined.

RESULTS

Optimal parameters of aTRUPC in neonates were found to be a TFE factor of 15, VENC of 5 cm/s, and slice thickness of 10 mm. The TFE-T₂ was on average 3.9% lower than standard T₂ . These two measures were strongly correlated (R² = 0.86); thus their difference can be accounted for by a correction equation, T_2,standard = 1.2002 × T_2,TFE - 10.6276. Neonatal Y_v values in veins draining cortical brain and those draining central brain were 64.8 ± 2.9% and 70.2 ± 3.3%, respectively, with a significant difference (P =.02).

CONCLUSION

The aTRUPC MRI has the potential to provide vessel-specific quantification of cerebral Y_v in neonates.

Oxygen metabolism in acute ischemic stroke

Gaining insights into brain oxygen metabolism has been one of the key areas of research in neurosciences. Extensive efforts have been devoted to developing approaches capable of providing measures of brain oxygen metabolism not only under normal physiological conditions but, more importantly, in various pathophysiological conditions such as cerebral ischemia. In particular, quantitative measures of cerebral metabolic rate of oxygen using positron emission tomography (PET) have been shown to be capable of discerning brain tissue viability during ischemic insults. However, the complex logistics associated with oxygen-15 PET have substantially hampered its wide clinical applicability. In contrast, magnetic resonance imaging (MRI)-based approaches have provided quantitative measures of cerebral oxygen metabolism similar to that obtained using PET. Given the wide availability, MRI-based approaches may have broader clinical impacts, particularly in cerebral ischemia, when time is a critical factor in deciding treatment selection. In this article, we review the pathophysiological basis of altered cerebral hemodynamics and oxygen metabolism in cerebral ischemia, how quantitative measures of cerebral metabolism were obtained using the Kety-Schmidt approach, the physical concepts of non-invasive oxygen metabolism imaging approaches, and, finally, clinical applications of the discussed imaging approaches.

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

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

Calibrated fMRI for mapping absolute CMRO2: Practicalities and prospects

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

Interleaved quantitative BOLD: Combining extravascular R2' - and intravascular R2-measurements for estimation of deoxygenated blood volume and hemoglobin oxygen saturation

Quantitative BOLD (qBOLD), a non-invasive MRI method for assessment of hemodynamic and metabolic properties of the brain in the baseline state, provides spatial maps of deoxygenated blood volume fraction (DBV) and hemoglobin oxygen saturation (HbO₂) by means of an analytical model for the temporal evolution of free-induction-decay signals in the extravascular compartment. However, mutual coupling between DBV and HbO₂ in the signal model results in considerable estimation uncertainty precluding achievement of a unique set of solutions. To address this problem, we developed an interleaved qBOLD method (iqBOLD) that combines extravascular R₂' and intravascular R₂ mapping techniques so as to obtain prior knowledge for the two unknown parameters. To achieve these goals, asymmetric spin echo and velocity-selective spin-labeling (VSSL) modules were interleaved in a single pulse sequence. Prior to VSSL, arterial blood and CSF signals were suppressed to produce reliable estimates for cerebral venous blood volume fraction (CBV_v) as well as venous blood R₂ (to yield HbO₂). Parameter maps derived from the VSSL module were employed to initialize DBV and HbO₂ in the qBOLD processing. Numerical simulations and in vivo experiments at 3 T were performed to evaluate the performance of iqBOLD in comparison to the parent qBOLD method. Data obtained in eight healthy subjects yielded plausible values averaging 60.1 ± 3.3% for HbO₂ and 3.1 ± 0.5 and 2.0 ± 0.4% for DBV in gray and white matter, respectively. Furthermore, the results show that prior estimates of CBV_v and HbO₂ from the VSSL component enhance the solution stability in the qBOLD processing, and thus suggest the feasibility of iqBOLD as a promising alternative to the conventional technique for quantifying neurometabolic parameters.

Hemodynamic and Metabolic Assessment of Neonates With Punctate White Matter Lesions Using Phase-Contrast MRI and T2-Relaxation-Under-Spin-Tagging (TRUST) MRI

The brain's hemodynamic and metabolism of punctate white matter lesions (PWML) is poorly understood due to a scarcity of non-invasive imaging techniques. The aim of this study was to apply new MRI techniques to quantify cerebral metabolic rate of oxygen (CMRO₂), global cerebral blood flow (CBF), oxygen saturation fractions in venous blood (Yv) and oxygen extraction fraction (OEF) in neonates with PWML, for better understanding of the pathophysiology of PWML. Fifty-one newborns were recruited continuously, including 23 neonatal patients with PWML and 28 normal control neonates. Phase-contrast (PC) MRI and T2-Relaxation-Under-Spin-Tagging (TRUST) MRI were performed for the measurement of CBF and Yv. OEF and CMRO₂ were calculated from the CBF and Yv values. The total maturation score (TMS) was assessed for each neonate on standard T1, 2-weighted images to evaluate cerebral maturation. The CMRO₂, CBF, Yv, and OEF values were compared between groups, and their associations with age and TMS were evaluated. Significant differences between PWML group and control group were found in CMRO₂ (P = 0.020), CBF (P = 0.027), Yv (P = 0.012), OEF (P = 0.018). After age/maturation is accounted for, Yv and OEF showed significant dependence on the groups (P < 0.05). Newborns with PWML had lower OEF and higher Yv. CMRO₂, CBF and brain volume were correlated with age (P < 0.001) and TMS (P < 0.05). It is feasible to use non-invasive MRI methods to measure cerebral oxygen supply and consumption in neonates with PWML. Newborns with PWML have lower oxygen consumption. Yv and OEF may be helpful for the diagnosis of PWML. The positive correlation between CBF and TMS, and between CMRO₂ and TMS suggested that as myelination progresses, the blood supply and oxygen metabolism in the brain increase to meet the escalating energy demand.

Quantitative assessment of cerebral venous blood T2 in mouse at 11.7T: Implementation, optimization, and age effect

PURPOSE

To develop a non-contrast-agent MRI technique to quantify cerebral venous T₂ in mice.

METHODS

We implemented and optimized a T₂ -relaxation-under-spin-tagging (TRUST) sequence on an 11.7 Tesla animal imaging system. A flow-sensitive-alternating-inversion-recovery (FAIR) module was used to generate control and label images, pair-wise subtraction of which yielded blood signals. Then, a T₂ -preparation module was applied to produce T₂ -weighted images, from which blood T₂ was quantified. We conducted a series of technical studies to optimize the imaging slice position, inversion slab thickness, post-labeling delay (PLD), and repetition time. We also performed three physiological studies to examine the venous T₂ dependence on hyperoxia (N = 4), anesthesia (N = 3), and brain aging (N = 5).

RESULTS

Our technical studies suggested that, for efficient data acquisition with minimal bias in estimated T₂ , a preferred TRUST protocol was to place the imaging slice at the confluence of sagittal sinuses with an inversion-slab thickness of 2.5-mm, a PLD of 1000 ms and a repetition time of 3.5 s. Venous T₂ values under normoxia and hyperoxia (inhaling pure oxygen) were 26.9 ± 1.7 and 32.3 ± 2.2 ms, respectively. Moreover, standard isoflurane anesthesia resulted in a higher venous T₂ compared with dexmedetomidine anesthesia (N = 3; P = 0.01) which is more commonly used in animal functional MRI studies to preserve brain function. Venous T₂ exhibited a decrease with age (N = 5; P < 0.001).

CONCLUSION

We have developed and optimized a noninvasive method to quantify cerebral venous blood T₂ in mouse at 11.7 T. This method may prove useful in studies of brain physiology and pathophysiology in animal models. Magn Reson Med 80:521-528, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Three-dimensional mapping of brain venous oxygenation using R2* oximetry

PURPOSE

Cerebral venous oxygenation (Y_v ) is an important biomarker for brain diseases. This study aims to develop an R2*-based MR oximetry that can measure cerebral Y_v in 3D.

METHODS

This technique separates blood signal from tissue by velocity-encoding phase contrast and measures the R2* of pure blood by multi-gradient-echo acquisition. The blood R2* was converted to Y_v using an R2*-versus-oxygenation (Y) calibration curve, which was obtained by in vitro bovine blood experiments. Reproducibility, sensitivity, validity, and resolution dependence of the technique were evaluated.

RESULTS

In vitro R2*-Y calibration plot revealed a strong dependence of blood R2* on oxygenation, with additional dependence on hematocrit. In vivo results demonstrated that the technique can provide a 3D venous oxygenation map that depicts both large sinuses and smaller cortical veins, with venous oxygenation ranging from 57 to 72%. Intrasession coefficient of variation of the measurement was 3.0%. The technique detected an average Y_v increase of 10.8% as a result of hyperoxia, which was validated by global oxygenation measurement from T₂ -Relaxation-Under-Spin-Tagging (TRUST) MRI. Two spatial resolutions, one with an isotropic voxel dimension and the other with a nonisotropic dimension, were tested for full brain coverage.

CONCLUSIONS

This study demonstrated the feasibility of 3D brain oxygenation mapping without using contrast agent. Magn Reson Med 79:1304-1313, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Clinical quantitative susceptibility mapping (QSM): Biometal imaging and its emerging roles in patient care

Quantitative susceptibility mapping (QSM) has enabled magnetic resonance imaging (MRI) of tissue magnetic susceptibility to advance from simple qualitative detection of hypointense blooming artifacts to precise quantitative measurement of spatial biodistributions. QSM technology may be regarded to be sufficiently developed and validated to warrant wide dissemination for clinical applications of imaging isotropic susceptibility, which is dominated by metals in tissue, including iron and calcium. These biometals are highly regulated as vital participants in normal cellular biochemistry, and their dysregulations are manifested in a variety of pathologic processes. Therefore, QSM can be used to assess important tissue functions and disease. To facilitate QSM clinical translation, this review aims to organize pertinent information for implementing a robust automated QSM technique in routine MRI practice and to summarize available knowledge on diseases for which QSM can be used to improve patient care. In brief, QSM can be generated with postprocessing whenever gradient echo MRI is performed. QSM can be useful for diseases that involve neurodegeneration, inflammation, hemorrhage, abnormal oxygen consumption, substantial alterations in highly paramagnetic cellular iron, bone mineralization, or pathologic calcification; and for all disorders in which MRI diagnosis or surveillance requires contrast agent injection. Clinicians may consider integrating QSM into their routine imaging practices by including gradient echo sequences in all relevant MRI protocols.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 5 J. Magn. Reson. Imaging 2017;46:951-971.

Quantitative measurement of cerebral blood volume using velocity-selective pulse trains

PURPOSE

To develop a non-contrast-enhanced MRI method for cerebral blood volume (CBV) mapping using velocity-selective (VS) pulse trains.

METHODS

The new pulse sequence applied velocity-sensitive gradient waveforms in the VS label modules and velocity-compensated ones in the control scans. Sensitivities to the gradient imperfections (e.g., eddy currents) were evaluated through phantom studies. CBV quantification procedures based on simulated labeling efficiencies for arteriolar, capillary, and venular blood as a function of cutoff velocity (Vc) are presented. Experiments were conducted on healthy volunteers at 3T to examine the effects of unbalanced diffusion weighting, cerebrospinal (CSF) contamination and variation of Vc.

RESULTS

Phantom results of the used VS pulse trains demonstrated robustness to eddy currents. The mean CBV values of gray matter and white matter for the experiments using Vc = 3.5 mm/s and velocity-compensated control with CSF-nulling were 5.1 ± 0.6 mL/100 g and 2.4 ± 0.2 mL/100 g, respectively, which were 23% and 32% lower than results from the experiment with velocity-insensitive control, corresponding to 29% and 25% lower in averaged temporal signal-to-noise ratio values.

CONCLUSION

A novel technique using VS pulse trains was demonstrated for CBV mapping. The results were both qualitatively and quantitatively close to those from existing methods. Magn Reson Med 77:92-101, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Empirical model of human blood transverse relaxation at 3 T improves MRI T2 oximetry

PURPOSE

We sought a human blood T₂ -oximetery calibration curve over the wide range of hematocrits commonly found in anemic patients applicable with T₂ relaxation under spin tagging (TRUST).

METHODS

Blood was drawn from five healthy control subjects. Ninety-three in vitro blood transverse relaxation (T_2b ) measurements were performed at 37°C over a broad range of hematocrits (10-55%) and oxygen saturations (14-100%) at 3 Tesla (T). In vivo TRUST was performed on 35 healthy African American control subjects and 11 patients with chronic anemia syndromes.

RESULTS

1/T₂ rose linearly with hematocrit (r²  = 0.96), for fully saturated blood. Upon desaturation, 1/T₂ rose linearly with the square of the oxygen extraction, (1-Y)² , and the slope was linearly proportional to hematocrit (r²  = 0.88). The resulting bilinear model between 1/T₂ , (1-Y)² , and hematocrit had a combined r² of 0.96 and a coefficient of variation of 6.1%. Using the in vivo data, the bilinear model had significantly lower bias and variability than existing calibrations, particularly for low hematocrits. In vivo Bland Altman analysis demonstrated clinically relevant bias that was -6% (absolute saturation) for hematocrits near 30% and rose to + 6% for hematocrits near 45%.

CONCLUSION

This work introduces a robust bilinear calibration model that should be used for MRI oximetry. Magn Reson Med 77:2364-2371, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

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

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

Non-invasive imaging of oxygen extraction fraction in adults with sickle cell anaemia

Sickle cell anaemia is a monogenetic disorder with a high incidence of stroke. While stroke screening procedures exist for children with sickle cell anaemia, no accepted screening procedures exist for assessing stroke risk in adults. The purpose of this study is to use novel magnetic resonance imaging methods to evaluate physiological relationships between oxygen extraction fraction, cerebral blood flow, and clinical markers of cerebrovascular impairment in adults with sickle cell anaemia. The specific goal is to determine to what extent elevated oxygen extraction fraction may be uniquely present in patients with higher levels of clinical impairment and therefore may represent a candidate biomarker of stroke risk. Neurological evaluation, structural imaging, and the non-invasive T2-relaxation-under-spin-tagging magnetic resonance imaging method were applied in sickle cell anaemia (n = 34) and healthy race-matched control (n = 11) volunteers without sickle cell trait to assess whole-brain oxygen extraction fraction, cerebral blood flow, degree of vasculopathy, severity of anaemia, and presence of prior infarct; findings were interpreted in the context of physiological models. Cerebral blood flow and oxygen extraction fraction were elevated (P < 0.05) in participants with sickle cell anaemia (n = 27) not receiving monthly blood transfusions (interquartile range cerebral blood flow = 46.2-56.8 ml/100 g/min; oxygen extraction fraction = 0.39-0.50) relative to controls (interquartile range cerebral blood flow = 40.8-46.3 ml/100 g/min; oxygen extraction fraction = 0.33-0.38). Oxygen extraction fraction (P < 0.0001) but not cerebral blood flow was increased in participants with higher levels of clinical impairment. These data provide support for T2-relaxation-under-spin-tagging being able to quickly and non-invasively detect elevated oxygen extraction fraction in individuals with sickle cell anaemia with higher levels of clinical impairment. Our results support the premise that magnetic resonance imaging-based assessment of elevated oxygen extraction fraction might be a viable screening tool for evaluating stroke risk in adults with sickle cell anaemia.

Measurement of oxygen extraction fraction (OEF): An optimized BOLD signal model for use with hypercapnic and hyperoxic calibration

Several techniques have been proposed to estimate relative changes in cerebral metabolic rate of oxygen consumption (CMRO2) by exploiting combined BOLD fMRI and cerebral blood flow data in conjunction with hypercapnic or hyperoxic respiratory challenges. More recently, methods based on respiratory challenges that include both hypercapnia and hyperoxia have been developed to assess absolute CMRO2, an important parameter for understanding brain energetics. In this paper, we empirically optimize a previously presented "original calibration model" relating BOLD and blood flow signals specifically for the estimation of oxygen extraction fraction (OEF) and absolute CMRO2. To do so, we have created a set of synthetic BOLD signals using a detailed BOLD signal model to reproduce experiments incorporating hypercapnic and hyperoxic respiratory challenges at 3T. A wide range of physiological conditions was simulated by varying input parameter values (baseline cerebral blood volume (CBV0), baseline cerebral blood flow (CBF0), baseline oxygen extraction fraction (OEF0) and hematocrit (Hct)). From the optimization of the calibration model for estimation of OEF and practical considerations of hypercapnic and hyperoxic respiratory challenges, a new "simplified calibration model" is established which reduces the complexity of the original calibration model by substituting the standard parameters α and β with a single parameter θ. The optimal value of θ is determined (θ=0.06) across a range of experimental respiratory challenges. The simplified calibration model gives estimates of OEF0 and absolute CMRO2 closer to the true values used to simulate the experimental data compared to those estimated using the original model incorporating literature values of α and β. Finally, an error propagation analysis demonstrates the susceptibility of the original and simplified calibration models to measurement errors and potential violations in the underlying assumptions of isometabolism. We conclude that using the simplified calibration model results in a reduced bias in OEF0 estimates across a wide range of potential respiratory challenge experimental designs.

A novel Bayesian approach to accounting for uncertainty in fMRI-derived estimates of cerebral oxygen metabolism fluctuations

Calibrated blood oxygenation level dependent (BOLD) imaging is a multimodal functional MRI technique designed to estimate changes in cerebral oxygen metabolism from measured changes in cerebral blood flow and the BOLD signal. This technique addresses fundamental ambiguities associated with quantitative BOLD signal analysis; however, its dependence on biophysical modeling creates uncertainty in the resulting oxygen metabolism estimates. In this work, we developed a Bayesian approach to estimating the oxygen metabolism response to a neural stimulus and used it to examine the uncertainty that arises in calibrated BOLD estimation due to the presence of unmeasured model parameters. We applied our approach to estimate the CMRO2 response to a visual task using the traditional hypercapnia calibration experiment as well as to estimate the metabolic response to both a visual task and hypercapnia using the measurement of baseline apparent R2' as a calibration technique. Further, in order to examine the effects of cerebral spinal fluid (CSF) signal contamination on the measurement of apparent R2', we examined the effects of measuring this parameter with and without CSF-nulling. We found that the two calibration techniques provided consistent estimates of the metabolic response on average, with a median R2'-based estimate of the metabolic response to CO2 of 1.4%, and R2'- and hypercapnia-calibrated estimates of the visual response of 27% and 24%, respectively. However, these estimates were sensitive to different sources of estimation uncertainty. The R2'-calibrated estimate was highly sensitive to CSF contamination and to uncertainty in unmeasured model parameters describing flow-volume coupling, capillary bed characteristics, and the iso-susceptibility saturation of blood. The hypercapnia-calibrated estimate was relatively insensitive to these parameters but highly sensitive to the assumed metabolic response to CO2.

Vessel-specific quantification of blood oxygenation with T2-relaxation-under-phase-contrast MRI

PURPOSE

Measurement of venous oxygenation (Yv) is a critical step toward quantitative assessment of brain oxygen metabolism, a key index in many brain disorders. The present study aims to develop a noninvasive, rapid, and reproducible method to measure Yv in a vessel-specific manner.

THEORY

The method, T2-Relaxation-Under-Phase-Contrast MRI, utilizes complex subtraction of phase-contrast to isolate pure blood signal, applies nonslice-selective T2-preparation to measure T2, and converts T2 to oxygenation using a calibration plot.

METHODS

Following feasibility demonstration, several technical aspects were examined, including validation with an established global Yv technique, test-retest reproducibility, sensitivity to detect oxygenation changes due to hypoxia and caffeine challenges, applicability of echo-planar-imaging (EPI) acquisition to shorten scan duration, and ability to study veins with a caliber of 1-2 mm.

RESULTS

T2-Relaxation-Under-Phase-Contrast was able to simultaneously measure Yv in all major veins in the brain, including sagittal sinus, straight sinus, great vein, and internal cerebral vein. T2-Relaxation-Under-Phase-Contrast results showed an excellent agreement with the reference technique, high sensitivity to oxygenation changes, and test-retest variability of 3.5 ± 1.0%. The use of segmented-EPI was able to reduce the scan duration to 1.5 minutes. It was also feasible to study pial veins and deep veins.

CONCLUSION

T2-Relaxation-Under-Phase-Contrast MRI is a promising technique for vessel-specific oxygenation measurement.

Multisite evaluations of a T2 -relaxation-under-spin-tagging (TRUST) MRI technique to measure brain oxygenation

PURPOSE

Venous oxygenation (Yv ) is an important index of brain physiology and may be indicative of brain diseases. A T2 -relaxation-under-spin-tagging (TRUST) MRI technique was recently developed to measure Yv . A multisite evaluation of this technique would be an important step toward broader availability and potential clinical utilizations of Yv measures.

METHODS

TRUST MRI was performed on a total of 250 healthy subjects, 125 from the developer's site and 25 each from five other sites. All sites were equipped with a 3 Tesla (T) MRI of the same vendor. The estimated Yv and the standard error (SE) of the estimation εYv were compared across sites.

RESULTS

The averaged Yv and εYv across six sites were 61.1% ± 1.4% and 1.3% ± 0.2%, respectively. Multivariate regression analysis showed that the estimated Yv was dependent on age (P = 0.009) but not on performance site. In contrast, the SE of the Yv estimation was site-dependent (P = 0.024) but was less than 1.5%. Further analysis revealed that εYv was positively associated with the amount of subject motion (P < 0.001) but negatively associated with blood signal intensity (P < 0.001).

CONCLUSION

This work suggests that TRUST MRI can yield equivalent results of Yv estimation across different sites.

Correction and optimization of a T2-based approach to map blood oxygenation in small cerebral veins

PURPOSE

Cerebral venous blood oxygenation (Yv ) is an important biomarker in brain physiology and function. The present study proposes a procedure to provide a quantitative map of the brain's intravascular Yv.

THEORY AND METHODS

The method is based on a pulse sequence, T2 -Relaxation-Under-Phase-Contrast (TRU-PC) MRI, with postprocessing approaches to correct eddy-current effects. A complete scan protocol consists of four TRU-PC scans sensitized to large and small vessels with anterior-posterior and foot-head flow-encoding directions, and the data are analyzed conjunctively. Eddy-current correction was performed by fitting the tissue phase to a hyperplane, and then subtracting the eddy-current phase from the measured vessel phase. The reproducibility of the Yv-maps was examined in five participants. Sensitivity of the Yv map to a caffeine challenge was studied in another five participants.

RESULTS

Removal of eddy-current induced artifact allowed for the correction of T2 measurements, as demonstrated in vivo and with simulation. A Yv-map depicting all vessels in the slice can be obtained with the proposed protocol. Test-retest variability of the Yv -map was 3.7 ± 1.2%. Yv reduction can be reliably detected (P < 0.001) following the caffeine ingestion.

CONCLUSION

With the proposed TRU-PC protocol and eddy-current correction procedure, an accurate, vessel-specific Yv map of the human brain can be obtained.

Rapid T2- and susceptometry-based CMRO2 quantification with interleaved TRUST (iTRUST)

Susceptometry-based oximetry (SBO) and T2-relaxation-under-spin-tagging (TRUST) are two promising methods for quantifying the cerebral metabolic rate of oxygen (CMRO2), a critical parameter of brain function. We present a combined method, interleaved TRUST (iTRUST), which achieves rapid, simultaneous quantification of both susceptometry- and T2-based CMRO2 via insertion of a flow-encoded, dual-echo gradient-recalled echo (OxFlow) module within the T1 recovery portion of the TRUST sequence. In addition to allowing direct comparison between SBO- and TRUST-derived venous oxygen saturation (Yv) values, iTRUST substantially improves TRUST temporal resolution for CMRO2 quantification and obviates the need for a separate blood flow measurement following TRUST acquisition. iTRUST was compared directly to TRUST and OxFlow alone in three resting subjects at baseline, exhibiting close agreement with the separate techniques and comparable precision. These baseline data as well as simulation results support the use of two instead of the traditional four T2 preparation times for T2 fitting, allowing simultaneous quantification of susceptometry- and T2-based Yv (and CMRO2) with three- and six-second temporal resolution, respectively. In 10 young healthy subjects, iTRUST was applied during a 5% CO2 gas mixture-breathing paradigm. T2-based Yv values were lower at baseline relative to susceptometry (62.3 ± 3.1 vs. 66.7 ± 5.1 %HbO2, P<0.05), but increased more in response to hypercapnia. As a result, T2-based CMRO2 decreased from 140.4 ± 9.7 to 120.0 ± 9.5 μMol/100g/min, a significant -14.6 ± 3.6% response (P < 0.0001), whereas susceptometry-based CMRO2 changed insignificantly from 123.4 ± 18.7 to 127.9 ± 25.7, a 3.3 ± 9.7% response (P = 0.31). These differing results are in accord with previous studies applying the parent OxFlow or TRUST sequences individually, thus supporting the reliability of iTRUST but also strongly suggesting that a systematic bias exists between the susceptometry- and T2-based Yv quantification techniques.

Global intravascular and local hyperoxia contrast phase-based blood oxygenation measurements

The measurement of venous cerebral blood oxygenation (Yv) has potential applications in the study of patient groups where oxygen extraction and/or metabolism are compromised. It is also useful for fMRI studies to assess the stimulus-induced changes in Yv, particularly since basal Yv partially accounts for inter-subject variation in the haemodynamic response to a stimulus. A range of MRI-based methods of measuring Yv have been developed recently. Here, we use a method based on the change in phase in the MR image arising from the field perturbation caused by deoxygenated haemoglobin in veins. We build on the existing phase based approach (Method I), where Yv is measured in a large vein (such as the superior sagittal sinus) based on the field shift inside the vein with assumptions as to the vein's shape and orientation. We demonstrate two novel modifications which address limitations of this method. The first modification (Method II), maps the actual form of the vein, rather than assume a given shape and orientation. The second modification (Method III) uses the intra and perivascular phase change in response to a known change in Yv on hyperoxia to measure normoxic Yv in smaller veins. Method III can be applied to veins whose shape, size and orientation are not accurately known, thus allowing more localised measures of venous oxygenation. Results demonstrate that the use of an overly fine spatial filter caused an overestimation in Yv for Method I, whilst the measurement of Yv using Method II was less sensitive to this bias, giving Yv = 0.62 ± 0.03. Method III was applied to mapping of Yv in local veins across the brain, yielding a distribution of values with a mode of Yv = 0.661 ± 0.008.