Sensitivity assessment of QSM+qBOLD (or QQ) in detecting elevated oxygen extraction fraction (OEF) in physiological change

Abstract

The study investigated the sensitivity of a novel MRI-based OEF mapping, quantitative susceptibility mapping plus quantitative blood oxygen level-dependent imaging (QSM+qBOLD or QQ), to physiological changes, particularly increased oxygen extraction fraction (OEF) by using hyperventilation as a vasoconstrictive stimulus. While QQ’s sensitivity to decreased OEF during hypercapnia has been demonstrated, its sensitivity to increased OEF levels, crucial for cerebrovascular disorders like vascular dementia and Parkinson’s disease, remains unexplored. In comparison with a previous QSM-based OEF, we evaluated QQ’s sensitivity to high OEF values. MRI data were obtained from 11 healthy subjects during resting state (RS) and hyperventilation state (HV) using a 3 T MRI with a three-dimensional multi-echo gradient echo sequence (mGRE) and arterial spin labeling (ASL). Region of interest (ROI) analysis and paired t-tests were used to compare OEF, CMRO₂ and CBF between QQ and QSM. Similar to QSM, QQ showed higher OEF during HV compared to RS: in cortical gray matter, QQ-OEF and QSM-OEF was 36.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 4.7% and 35.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 12.5% at RS and 45.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 11.6% and 45.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 14.8% in HV, respectively. These findings demonstrate QQ’s ability to detect physiological changes and suggest its potential in studying brain metabolism in neurological disorders.

Introduction

Cerebral metabolic rate of oxygen (CMRO₂) and oxygen extraction fraction (OEF) have been considered as vital biomarkers for assessing tissue viability and activity.^1,2 Recent advancements in MRI techniques have made their mapping possible for diagnosing and evaluating neurological conditions such as ischemic stroke, ³ Alzheimer disease, ⁴ and tumor. ⁵ Primarily, MRI methods have focused on the magnitude signal modeling for OEF estimation, including calibrated fMRI,^{6 –9} quantitative imaging of extraction of oxygen and tissue consumption (QUIXOTIC), ¹⁰ quantitative blood oxygen level dependent imaging (qBOLD). ¹¹ However, novel OEF techniques that model the often-neglected MRI phase signal have been introduced, including susceptometry-based oximetry,^12,13 and quantitative susceptibility mapping (QSM)-based OEF.^14,15

QSM-based OEF methods initially focused on estimating macrovascular OEF based on the relationship between venous oxygenation and susceptibility to infer regional OEF information.^{16 –19} However, issues such as low spatial specificities (i.e., difficulties in identifying accurate brain regions draining into the corresponding veins) and partial volume effects have promoted the development of microvascular QSM-based OEF techniques.^{16 –20} Microvascular QSM-based OEF mapping, here referred to as QSM-OEF, directly measures regional tissue oxygen metabolism by separating susceptibility distributions of microvasculature (e.g., capillaries) from non-blood tissue within the cerebral parenchyma. ¹⁸ Hence, QSM-OEF enables OEF quantification at the voxel level.^{16 –20} For the susceptibility separation, QSM-OEF utilizes vascular challenges, such as hyperventilation. ¹⁹ Hyperventilation reduces the partial pressure of carbon dioxide (CO₂) in arterial blood leading to vasoconstriction and a subsequent decrease in cerebral blood flow (CBF), while cerebral metabolic rate of oxygen (CMRO₂) remains relatively stable.^21,22 However, this vascular challenge procedure is impractical in clinical settings. Moreover, QSM-OEF assumes that venous blood volume (v) is known from its linear relationship with CBF, which was derived from the whole-brain linear regression analysis in healthy subjects. ²³ This assumption may not be valid at individual voxel level or in disease.

Recently, a novel OEF mapping technique, termed QQ, has been introduced, which removes the need for the vascular challenges and the unrealistic v-CBF assumption in QSM-OEF. QQ integrates QSM-OEF technique with qBOLD for comprehensive MRI signal modeling, utilizing both phase and magnitude signal (QSM + qBOLD = QQ). ²⁴ QQ estimates OEF maps from multi-echo gradient (mGRE) data alone in a single brain state (e.g., resting state), whereas QSM-OEF requires both mGRE and arterial spin labeling (ASL) in two brain states (e.g., resting state and hyperventilation). Furthermore, the robustness of QQ against noise has been significantly improved by using unsupervised machine learning algorithms.^{25 –29} QQ has been validated against ¹⁵O-PET showing a great scan-rescan producibility in healthy subjects, ³⁰ tested in other imaging sequences, ³¹ and applied in neurologic disorders.^{32 –42} QQ utilizes a practical data acquisition scheme, i.e., it uses a single routine sequence, mGRE, available in every MRI scanner and does not require any vascular challenges, unlike other OEF techniques that require specialized sequences and vascular challenges.^{6 –10} Consequently, QQ is readily applicable in clinical settings.

In this study, we aim to investigate the sensitivity of QQ to physiological changes, particularly to increased OEF. While a comparison study with calibrated fMRI demonstrated QQ’s sensitivity to decreased OEF driven by hypercapnia, ⁴³ its sensitivity to increased OEF values has yet to be explored, which is crucial for investigating cerebrovascular disorders associated elevated OEF, including vascular dementia and Parkinson’s disease, where higher OEF is expected due to impaired CBF.^{44 –47} Comparing with QSM-OEF, we evaluated the sensitivity of QQ to increased OEF values.

Materials and methods

Data acquisition

The dataset utilized in this study was retrospectively obtained from previous QSM-OEF study which utilized resting state and hyperventilation. ¹⁹ This retrospective image analysis study received approval from the Weill Cornell Medicine local Institutional Review Board, and all participants provided written informed consent in accordance with the ethical standards of the 1975 Helsinki Declaration and its subsequent amendments. Healthy volunteers (n = 11, 1 female, mean age 34 ± 12 years) underwent brain MRI scanning twice: once while at resting state (RS) and again while the person was hyperventilating (HV). MRI scan was perform using a 3 Tesla scanner (HDxt, GE Healthcare, Waukesha, Wisconsin) utilizing an eight-channel brain receive coil.

Hyperventilation experiment

Participants were guided to hyperventilate at a specific rate (one inhale and exhale for every 2 to 3 seconds) prior to their scans. MRI scans were conducted both before and during hyperventilation using a protocol that included a 3D FSE ASL sequence and a 3D SPGR sequence. The scans during the hyperventilation state began 30 seconds after the start of hyperventilation. The hyperventilation period lasts for about 7 minutes. To minimize head motion, restrictive padding was used inside the head coil. Consequently, no noticeable motion artifacts were observed in either the ASL or SPGR datasets.

The imaging protocol includes an anatomical T1w BRAVO sequence, a 3D fast spin echo (FSE) ASL for CBF measurement, and a 3D mGRE sequence for OEF estimation. The parameters for the 3D FSE ASL sequence were as follows: a labeling period of 1500-ms, post label delay of 1525-ms, spiral sampling of eight interleaves with 512 readout points per leaf, readout bandwidth of 62.5-kHz, a 20-cm field of view (FOV), 3-mm isotropic resolution, 35 axial slices, 10.5-ms echo time (TE), 4796-ms pulse repetition time (TR), three signal averages, and 6-min scan time. The parameters for the 3D SPGR sequence includes seven equally spaced echoes, with a first TE of 4.3-ms, echo spacing of 7.9-ms, TR of 56.6-ms, 0.52-mm in plane resolution, 2-mm slice thickness, identical volume coverage as the 3D FSE ASL sequence, readout bandwidth of 62.5-kHz, 15 $°$ flip angle, and a scan time of 7-min. The pulse sequence was flow compensated in all three directions. ⁴⁸ The parameters of T1 BRAVO sequence included a TE of 2.92-ms, TR of 7.68-ms, preparation time 450-ms, 1.2-mm isotropic resolution, and identical volume coverage as the 3D FSE ASL sequence, 19.5-kHz readout bandwidth, 15 $°$ flip angle and a scan time of 2-min.

Data processing

CBF and QSM processing

CBF maps (ml/100 g/min) are created from the ASL data using the FuncTool software package (GE Healthcare, Waukesha, Wisconsin). QSM maps were calculated from mGRE data. Initially, susceptibility estimation involves obtaining the total field through a linear fit of the mGRE phase. ⁴⁹ Then, the local field is estimated via the projection onto dipole field (PDF) method. ⁵⁰ Subsequently, susceptibility is computed using the Morphology Enabled Dipole Inversion with automatic uniform cerebrospinal fluid zero reference (MEDI + 0) algorithm.^{15,51 –53}

QQ-OEF and CMRO₂ mapping from mGRE data at baseline and at hyperventilation separately

The QQ model estimates OEF based on the venous deoxyheme-dependent signal in the mGRE phase signal using QSM and the magnitude signal using qBOLD. ²⁵ To obtain the OEF, QQ solves for venous oxygenation ( $\mathrm{Y}$ ) using OEF $=1-\mathrm{Y}/{\mathrm{Y}}_{\mathrm{a}}$ where ${\mathrm{Y}}_{\mathrm{a}}$ is arterial oxygenation ( ${\mathrm{Y}}_{\mathrm{a}}$  = 0.98):

$$\begin{array}{c}\mathrm{argmin}\\ \mathrm{Y},v,{\mathrm{R}}_2,{\mathrm{S}}_0,{\mathrm{\chi }}_{\mathrm{n}}\end{array}\left\{\mathrm{w}{‖\mathrm{\chi }-{\mathrm{F}}_{\mathrm{QSM}}(\mathrm{Y},v,{\mathrm{\chi }}_{\mathrm{n}})‖}_2^2+\sum _{\mathrm{j}}{‖\left|{\mathrm{S}}_{\mathrm{j}}\right|-{\mathrm{F}}_{\mathrm{qBOLD}}({\mathrm{S}}_0,\mathrm{Y},v,{\mathrm{R}}_2,{\mathrm{\chi }}_{\mathrm{n}},{\mathrm{t}}_{\mathrm{j}})‖}_2^2+\mathrm{R}\left(\mathrm{Y}\right)\right\}$$

(1)

where w is weighting on the first QSM term.

In the first term of equation (1), the QSM approach separates voxel wise susceptibility ( $\mathrm{\chi }$ ) into blood susceptibility (which includes $\mathrm{Y}$ ) and non-blood neural tissue susceptibility ( ${\mathrm{\chi }}_{\mathrm{n}}$ ).

$${\mathrm{F}}_{\mathrm{QSM}}\left(\mathrm{Y},v,{\mathrm{\chi }}_{\mathrm{n}}\right)=\left[\frac{{\mathrm{\chi }}_{\mathrm{ba}}}{\mathrm{\alpha }}+{\mathrm{\psi }}_{\mathrm{Hb}}\cdot \mathrm{\Delta }{\mathrm{\chi }}_{\mathrm{Hb}}\cdot \left(-\mathrm{Y}+\frac{1-\left(1-\mathrm{\alpha }\right)\cdot {\mathrm{Y}}_{\mathrm{a}}}{\mathrm{\alpha }}\right)\right]\cdot v+\left(1-\frac{v}{\mathrm{\alpha }}\right)\cdot {\mathrm{\chi }}_{\mathrm{n}}$$

(2)

where ${\mathrm{\chi }}_{\mathrm{ba}}=-108.3\hspace{0.17em}$ ppb is the fully oxygenated blood susceptibility, ¹⁹ ${\mathrm{\psi }}_{\mathrm{Hb}}=0.0909$ is the hemoglobin volume fraction assuming Hct = 0.357,^54,55 $\mathrm{\Delta }{\mathrm{\chi }}_{\mathrm{Hb}}$  = 12522 ppb is the susceptibility difference between deoxyhemoglobin and oxyhemoglobin,^18,56 $\mathrm{\alpha }=0.77$ is the ratio between the venous and total blood volumes. ⁵⁷

In the second term of equation (1), qBOLD models the OEF effect on the magnitude signal at the echo time, ${\mathrm{t}}_{\mathrm{j}}$ .

$${\mathrm{F}}_{\mathrm{qBOLD}}\left({\mathrm{S}}_0,{\mathrm{R}}_2,\mathrm{Y},v,{\mathrm{\chi }}_{\mathrm{n}},{\mathrm{t}}_{\mathrm{j}}\right)={{\mathrm{S}}_0\cdot \mathrm{e}}^{-{\mathrm{R}}_2\cdot {\mathrm{t}}_{\mathrm{j}}}\cdot \mathrm{F}\mathrm{BOLD}\mathrm{ }\left(\mathrm{Y},v,{\mathrm{\chi }}_{\mathrm{n}},{\mathrm{t}}_{\mathrm{j}}\right)\cdot \mathrm{G}\left({\mathrm{t}}_{\mathrm{j}}\right)$$

(3)

where ${\mathrm{S}}_0$ is the signal intensity at ${\mathrm{t}}_{\mathrm{j}}=0$ , ${\mathrm{R}}_2$ is the transverse relaxation rate, $\mathrm{F}\mathrm{BOLD}\hspace{0.17em}\left(\mathrm{Y},v,{\mathrm{\chi }}_{\mathrm{n}},\mathrm{t}\right)=\exp \left[-v\cdot {\mathrm{f}}_{\mathrm{s}}\left(\mathrm{\delta }\mathrm{\omega }\cdot \mathrm{t}\right)\right]$ ⁵⁸ where ${\mathrm{f}}_{\mathrm{s}}$ is the signal decay due to presence of blood vessel network, $\mathrm{\delta }\mathrm{\omega }\left(\mathrm{Y},{\mathrm{\chi }}_{\mathrm{n}}\right)=\frac{1}{3}\cdot \mathrm{\gamma }\cdot {\mathrm{B}}_0\cdot \left[{\mathrm{\psi }}_{\mathrm{Hb}}\cdot \mathrm{\Delta }{\mathrm{\chi }}_{\mathrm{Hb}}\cdot \left(1-\mathrm{Y}\right)+{\mathrm{\chi }}_{\mathrm{ba}}-{\mathrm{\chi }}_{\mathrm{n}}\right]$ is the characteristic frequency by susceptibility difference between deoxygenated blood and surrounding tissue ²⁴ where $\mathrm{\gamma }$  = 267.51 (rad $/$ s)/T, ${\mathrm{B}}_0$ is the main magnetic field strength, and $\mathrm{G}\left({\mathrm{t}}_{\mathrm{j}}\right)$ is the macroscopic field inhomogeneity to mGRE signal estimated using voxel spread function (see Appendix of Cho et al. ²⁴ ) In the third term of equation (1), $\mathrm{R}\left(\mathrm{Y}\right)$ is the regularization on $\mathrm{Y}$ , such as total variation regularization which can suppress noise propagation in $\mathrm{Y}$ .

In this study, we applied a machine learning algorithm, termed temporal clustering, tissue composition, and total variation (CCTV), to solve QQ robustly. First, CCTV utilized temporal clustering to group voxels with similar mGRE signals into the same clusters by using x-means algorithm, an improved version of k-means that automatically determines the optimal number of clusters. This approach is based on that voxels exhibiting similar mGRE signal decay have similar model parameter values. Then, each cluster was further divided into three sub-clusters (gray matter (GM), white matter (WM), and cerebrospinal fluid (CSF)) through the tissue composition process in CCTV, based on that the same type of tissue has similar tissue properties (including v and R₂). This tissue type integration allows comprehensive clustering, which leads to enhanced optimization initialization. ²⁶ Since the number of clusters is typically much smaller than the number of voxels (∼50 vs.10⁶), averaging within these clusters can significantly enhance the signal-to-noise ratio (SNR). ²⁵ Note that, above a certain number of clusters (e.g., 10 clusters), the number of clusters does not significantly affect the resultant OEF map, ²⁵ which suggests that uniform OEF map from QQ-CCTV in this study may not be caused by the clustering, but it primarily represents the expected uniform oxygen extraction throughout the whole brain as seen in previous ¹⁵O-PET studies.^30,59,60 Additionally, CCTV applies total variation regularization to minimize noise propagation into the OEF map. ⁶¹

Using the OEF, CMRO₂ was estimated from the equation CMRO₂ = OEF × CBF × [H]_a, where [H]_a = 7.377 µmol/mL is the oxygenated heme molar concentration in an arteriole with haematocrit (Hct) = 0.357. ¹⁷

QSM-OEF mapping at baseline and hyperventilation utilizing mGRE and ASL data in both resting state and hyperventilation together

The QSM-OEF model was implemented according to Zhang et al. ¹⁹ QSM-OEF solves for three unknowns (resting state OEF, hyperventilation OEF, and ${\chi }_n$ ) using three equations: two susceptibility equations (which corresponds to RS and HV, respectively) and one regularization on CMRO₂. In the susceptibility equations, the voxel-wise susceptibility is set as the sum of ${\chi }_n$ and blood susceptibility for each RS and HV. ${\chi }_n$ is assumed to be invariant between the two states. In the CMRO₂ regularization equation, the ratio of CMRO₂ between RS and HV in each individual voxel should match the global CMRO₂ estimated from venous oxygenation in a large draining vein, straight sinus.

ROI analysis

Cortical gray matter (CGM), deep gray matters (DGMs; thalamus, caudate, putamen, pallidum), and white matter (WM) masks were constructed by using FreeSurfer algorithm ⁶² on T1w BRAVO images, which were visually confirmed by a neuroradiologist (A.G.). CSF voxels were then excluded from the ROIs using a CSF mask, which was created by thresholding R2* values below 10 Hz and expanding it by approximately 2 voxels at the boundary to include partial-volumed edge voxels. For each ROI, the mean and standard deviation of OEF, CBF, and CMRO₂ from QQ and QSM technique were measured in both before RS and during HV. Paired sample t-tests and Bland-Altman (BA) plots were performed to compare the differences in OEF, CBF, and CMRO_2, in CGM during RS and HV across the subjects. The Shapiro-Wilk test confirmed the normality of each ROI data distribution (all p-values >0.05). A P-value less than 0.05 was considered statistically significant. P-values are adjusted using False Discovery Rate (FDR) correction to address multiple comparison. All data were expressed as mean ± standard deviation.

Results

Figure 1 shows OEF maps in two representative subjects at RS and HV using QSM and QQ. Compared to QSM, QQ showed a more uniform OEF. In CGM, QQ showed similar average OEF values compared to QSM at both RS and HV. Specifically, the average OEF values for the 11 healthy participants at RS were 36.4 ± 4.7% for QQ and 35.3 ± 12.5% for QSM (P = 0.51), and at HV were 45.0 ± 11.6% for QQ and 45.0 ± 14.8% for QSM (P = 0.98), as detailed in Table 1. Both QSM and QQ showed significantly increased OEF at HV compared to RS: QQ (P = 0.005) and QSM (P < 0.001).

Figure 1.

OEF maps from QSM and QQ between RS and HV in two subjects. QQ shows more uniform OEF than QSM. Both QQ and QSM shows decreased OEF in RS when compared to HV. The unit of color scale is %. The corresponding T1 BRAVO images are shown below for each subject.

Table 1.

Average OEF, CMRO₂ and CBF in CGM.

	OEF				CMRO2
	QSM		QQ		QSM		QQ		CBF
Subject	RS	HV	RS	HV	RS	HV	RS	HV	RS	HV
1	39.2 $\hspace{0.17em}\pm \hspace{0.17em}$ 15.8	44.6 $\hspace{0.17em}\pm \hspace{0.17em}$ 17.4	38.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 4.5	40.6 ± 11.8	164.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 75.4	155.2 $\hspace{0.17em}\pm \hspace{0.17em}$ 76.5	163.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 41.7	140.1 $\hspace{0.17em}\pm \hspace{0.17em}$ 54.5	57.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 13.1	45.6 $\hspace{0.17em}\pm \hspace{0.17em}$ 14.2
2	37.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 14.7	42.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 17.1	36.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 5.0	49.2 ± 15.8	186.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 88.0	156.1 $\hspace{0.17em}\pm \hspace{0.17em}$ 80.2	180.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 45.8	176.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 71.3	67.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 15.0	47.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 14.9
3	34.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 10.3	35.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 10.1	30.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 1.9	36.6 ± 3.1	162.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 68.1	134.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 55.5	142.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 34.4	133.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 32.4	63.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 15.1	49.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 11.4
4	35.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 15.0	48.6 $\hspace{0.17em}\pm \hspace{0.17em}$ 17.5	42.2 $\hspace{0.17em}\pm \hspace{0.17em}$ 6.1	46.9 ± 10.9	161.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 105.4	145.1 $\hspace{0.17em}\pm \hspace{0.17em}$ 102.6	180.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 59.9	122.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 63.8	58.2 $\hspace{0.17em}\pm \hspace{0.17em}$ 17.5	35.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 17.2
5	31.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 13.7	45.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 16.5	39.6 $\hspace{0.17em}\pm \hspace{0.17em}$ 5.5	62.0 ± 18.7	128.2 $\hspace{0.17em}\pm \hspace{0.17em}$ 72.5	125.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 78.9	156.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 50.1	157.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 73.5	53.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 15.9	34.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 13.7
6	37.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 13.7	56.1 $\hspace{0.17em}\pm \hspace{0.17em}$ 20.3	35.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 4.7	32.0 ± 14.2	160.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 71.3	152.6 $\hspace{0.17em}\pm \hspace{0.17em}$ 76.6	149.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 36.8	83.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 43.5	57.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 11.8	35.1 $\hspace{0.17em}\pm \hspace{0.17em}$ 10.3
7	35.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 10.7	44.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 12.6	36.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 4.1	45.1 ± 9.1	139.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 64.8	129.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 70.1	139.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 40.4	123.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 48.0	51.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 14.0	37.1 $\hspace{0.17em}\pm \hspace{0.17em}$ 12.7
8	30.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 11.4	37.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 12.0	38.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 4.2	44.4 ± 10.0	163.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 83.9	122.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 69.7	196.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 55.4	139.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 57.5	70.1 $\hspace{0.17em}\pm \hspace{0.17em}$ 18.4	42.6 $\hspace{0.17em}\pm \hspace{0.17em}$ 15.5
9	40.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 11.3	46.2 $\hspace{0.17em}\pm \hspace{0.17em}$ 12.8	40.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 8.6	49.5 ± 13.3	175.1 $\hspace{0.17em}\pm \hspace{0.17em}$ 78.4	154.2 $\hspace{0.17em}\pm \hspace{0.17em}$ 74.1	169.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 56.2	160.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 61.1	56.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 14.1	42.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 12.9
10	32.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 10.1	52.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 14.3	32.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 3.5	51.5 ± 12.5	130.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 61.7	139.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 73.4	127.2 $\hspace{0.17em}\pm \hspace{0.17em}$ 37.6	128.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 51.4	52.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 14.5	33.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 11.3
11	32.6 $\hspace{0.17em}\pm \hspace{0.17em}$ 10.3	42.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 11.6	29.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 3.4	36.7 ± 7.8	164.1 $\hspace{0.17em}\pm \hspace{0.17em}$ 72.3	162.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 67.4	146.7 $\hspace{0.17em}\pm \hspace{0.17em}$ 38.6	138.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 44.4	66.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 15.9	50.9 $\hspace{0.17em}\pm \hspace{0.17em}$ 12.9
Avg $\hspace{0.17em}\pm \hspace{0.17em}$ SD	35.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 12.5	45.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 14.8	36.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 4.7	45.0 ± 11.6	157.8 $\hspace{0.17em}\pm \hspace{0.17em}$ 76.5	143.5 $\hspace{0.17em}\pm \hspace{0.17em}$ 75.0	159.2 $\hspace{0.17em}\pm \hspace{0.17em}$ 45.2	136.6 $\hspace{0.17em}\pm \hspace{0.17em}$ 54.7	59.6 $\hspace{0.17em}\pm \hspace{0.17em}$ 15.0	41.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 13.4

Note: The unit of OEF, CMRO₂ and CBF is %, µmol/min/100 g and mL/100 g/min respectively.

RS: resting state; HV: hyperventilation; OEF: oxygen extraction fraction; CMRO₂: cerebral metabolic rate of oxygen; CBF: cerebral blood flow; QSM: quantitative susceptibility mapping; QQ: quantitative susceptibility mapping + quantitative blood oxygen level dependent imaging.

Figure 2 shows CMRO₂ and CBF maps in the same two subjects at RS and HV for QSM and QQ. QQ showed clear CMRO₂ contrast CGM and WM. In CGM, QQ showed lower CMRO₂ value to that of QSM: during RS, the average CMRO₂ values for the 11 healthy participants were 159.2 ± 45.2 µmol/100 g/min for QQ and 157.8 ± 76.5 µmol/100 g/min for QSM (P = 0.86), and during HV were 136.6 ± 54.7 µmol/100 g/min for QQ and 143.5 ± 75.0 µmol/100 g/min for QSM (P = 0.52). Both QSM and QQ showed a decrease in CMRO₂ during HV compared to RS: QQ (P = 0.0199) and QSM (P = 0.0139). CBF was higher during RS compared to HV, with the average CBF values for the 11 subjects being 59.6 ± 15.0 mL/100 g/min for RS and 41.4 ± 13.4 mL/100 g/min for HV (P < 0.05) in CGM.

Figure 2.

CMRO₂ and CBF maps from QSM and QQ between RS and HV in the same two subjects as in Figure 1. QQ shows good CMRO₂ contrast between CGM and WM. The CBF was obtained by registering to OEF resolution. The unit of color scale is mol/100 g/min for CMRO₂ and mL/min/100 g for CBF.

Figure 3 shows boxplots of OEF, CMRO₂ and CBF in the ROIs. For OEF, both QQ and QSM show significantly higher OEF in HV compared to RS across all ROIs. Compared to QSM, QQ shows similar OEF values in CGM and Thalamus, but higher values in WM and lower values in Caudate, Putamen, and Pallium. For CMRO₂, no significant difference was found between RS and HV for both QSM and QQ in all ROIs, except CGM, where QQ shows lower CMRO₂ at HV. Compared to QSM, QQ shows similar CMRO₂ values in CGM and Thalamus, but higher values in WM and lower values in Caudate, Putamen, and Pallidum. CBF is lower in RS compared to HV in all the ROIs. Detailed ROI analysis results (i.e., mean ± standard deviation in each ROI) are provided in Supporting Information Tables S1, S2, and S3.

Figure 3.

Box plots of OEF, CMRO₂ and CBF in CGM were plotted for QQ and QSM. Red line, blue box, black whisker, and red cross indicates the median, interquartile range, and the range extending to 1.5 of the interquartile range, outlier beyond the whisker range. The significant differences (P < 0.05) between RS and HV approach are marked with an asterisk.

Figure 4 shows the Bland-Altman plots comparing the QSM and QQ approach for OEF and CMRO₂ across all subjects. All four comparisons show good agreement in CGM and Thalamus. For example, in CGM, the mean difference for OEF was not significant between QSM and QQ, 7% (P = 0.51) at RS and 5% (P = 0.98) at HV, and the mean difference for CMRO₂ were 29.29 mol/100 g/min (P = 0.86) at RS and 20.59 mol/100 g/min (P = 0.52) at HV. On the other hand, in Caudate, Putamen, and Pallidum, QQ and QSM showed significant biases. For example, in the Caudate, the mean difference for OEF was 15.9% (P < 0.001) at RS and 16.4% (P = 0.002) at HV. Similarly, in the Pallidum, the mean difference for OEF was 23.6% (P < 0.001) at RS and 15.4% (P = 0.003) at HV. The biases were also prominent for CMRO₂, with the Caudate showing a mean difference of 64.2% (P < 0.001) at RS and 52.7% (P < 0.001) at HV.

Figure 4.

Bland-Altman plots comparing OEF and CMRO₂ values in CGM between QQ and QSM. (a) QQ-OEF vs. QSM-OEF at RS, (b) QQ-OEF vs. QSM-OEF during, (c) HV QQ-CMRO₂ vs. QSM-CMRO₂ at RS, and (d) QQ-CMRO₂ vs. QSM-CMRO₂ during HV.

Discussion

This study investigated the sensitivity of QQ to physiologically increased OEF during hyperventilation, in comparison with the microvascular QSM-based OEF technique (QSM-OEF). As abnormally high OEF is expected in a variety of cerebrovascular disorders, including vascular dementia,^{44 –46} Parkinson’s disease, ⁴⁷ and moyamoya disease,^{63 –65} this sensitivity study demonstrates the QQ’s feasibility for such disorders. Our results indicate that QQ provided similar OEF values in CGM compared to QSM-OEF at resting state (RS) and significantly higher OEF during hyperventilation (HV) compared to RS. QQ has higher potential in clinical settings than QSM-OEF in three aspects. First, QQ is based on more accurate biophysics modeling than QSM-OEF, as it eliminates the unreliable assumption in QSM-OEF that venous blood volume (v) is linearly proportional to CBF in every voxel. Our QQ provided significantly decreased v during HV compared to RS, e.g., 3.94 ± 0.27% vs. 3.71 ± 0.63% (P = 0.0029) in CGM (Supporting Information Figure S1), which demonstrates that QQ can measure vasoconstriction induced by hyperventilation. Second, QQ does not require a vascular challenge (i.e., it only needs a single brain state, e.g. RS), whereas QSM-OEF requires a vascular challenge, e.g. both RS and HV, to estimate OEF during RS. Such vascular challenge is not practical in clinical settings. Third, the reconstruction of QQ may be more robust due to its utilization of both phase and magnitude signals and its implementation of a machine learning algorithm, CCTV. The sensitivity of QQ to high OEF abnormalities, in addition to its previously observed sensitivity to low OEF abnormalities, ⁴³ guarantees QQ’s wide feasibility in various clinical settings.

OEF during resting state in CGM estimated from the two methods were similar and statistically not significantly different (Table 1): 36.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 4.7% and 35.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 12.5% (P = 0.51) for QSM-OEF and QQ-OEF, respectively. Both the OEF values fall within the range previously reported using QSM: 29 $\hspace{0.17em}\pm \hspace{0.17em}$ 3%∼50 $\hspace{0.17em}\pm \hspace{0.17em}$ 5% (Table 4 in literature ²⁰ ), using PET: 35 $\hspace{0.17em}\pm \hspace{0.17em}$ 7% ⁶⁶ and 40 $\hspace{0.17em}\pm \hspace{0.17em}$ 9%, ⁶⁷ and using calibrated BOLD: 35 $\hspace{0.17em}\pm \hspace{0.17em}$ 4% ∼44 $\hspace{0.17em}\pm \hspace{0.17em}$ 14% (Table 3 in literature ²⁰ ). On the other hand, in deep gray matters, including caudate, putamen, and pallidum, as well as WM, QSM-OEF showed different OEF values from QQ, e.g., abnormally high OEF values in deep gray matters (e.g., OEF > 60%). These abnormal OEF in QSM-OEF may be due to low SNR in ASL measurements (i.e., inaccurate CBF estimation) in those regions, ¹⁹ which leads to inaccurate v estimation and subsequently inaccurate OEF estimation. Additionally, QSM values have high uncertainties in brain regions with high susceptibility, including caudate and putamen. ⁵² These uncertainties may propagate to OEF estimation in the QSM-OEF method, whereas in QQ, the uncertainty propagation was suppressed by the inclusion of qBOLD and the application of machine learning algorithm, CCTV, which effectively enhances SNR. ²⁶ Similarly, greater intersubject OEF variation in QSM-OEF (coefficient of variation in CGM = 13% for QQ vs 35% for QSM) may be attributed to the involvement of low SNR CBF signals in QSM-OEF processing, ¹⁹ whereas the smaller OEF variation in QQ could have benefited from the utilization of CCTV. ²⁶

In addition, QQ showed relatively uniform OEF between CGM and WM (Figure 1), which agrees well with previous PET studies,^60,68,69 whereas QSM-OEF showed OEF contrast between CGM and WM. This contrast may be caused by the inclusion of CBF in QSM-OEF. QQ computes OEF independently for RS and HV without utilizing CBF information. In contrast, QSM-OEF incorporates CBF information for OEF estimation. It treats v as a constant, not as a model parameter, by assuming that v is linearly related to CBF. However, CBF estimation in WM may not be accurate owing to low SNR in ASL measurements. ⁷⁰ This leads to inaccurate v estimation in WM and subsequently inaccurate OEF estimation in WM using QSM-OEF approach. This inaccuracy in OEF estimation in WM may cause the observed OEF contrast between CGM and WM in QSM-OEF.

In comparison between RS and HV, the two methods showed similar OEFs in CGM, e.g., 36.4 $\hspace{0.17em}\pm \hspace{0.17em}$ 4.7% and 35.3 $\hspace{0.17em}\pm \hspace{0.17em}$ 12.5% (P = 0.51) at RS and 45.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 11.6% and 45.0 $\hspace{0.17em}\pm \hspace{0.17em}$ 14.8% (P = 0.98) at HV (Table 1). Furthermore, similar to QSM, QQ demonstrated higher OEF at HV than at RS. This suggests that QQ is sensitive to high OEF, which is crucial for investigating cerebrovascular disorders, including vascular dementia where high OEF is expected.^{44 –47} The higher OEF observed in hyperventilation agrees with previous vascular challenge studies that have shown an increase in OEF during vasoconstriction.^19,21,71 Hyperventilation induces vasoconstriction ¹⁹ (i.e., decrease in blood flow), which leads to increase in OEF. Together with our previous study that showed QQ’s sensitivity to low OEF, ⁴³ these results suggest that QQ can be effectively used to investigate both increased and decreased OEF abnormalities in various neurologic disorders.^{32 –36,43 –47} The OEF in the DGMs also shows increased OEF during hyperventilation compared to resting state.

In CGM, QQ and QSM showed a significant CMRO₂ decrease during HV compared to RS. The CGM CMRO₂ decrease during HV is consistent with previous studies. For instance, one study observed a substantial decrease in CMRO₂ with the immediate initiation of hyperventilation when CBF decreases. ⁷² Similar studies have also shown reduction in CMRO₂ with the onset of hypocarbia.^{73 –77} In this study, QQ showed a 14.3% CMRO₂ reduction in CGM, which agrees with a CMRO₂ reduction of 17.4% in hypercapnia in a similar study. ⁴³ This is in line with recent studies suggesting that blood flow changes, such as hyper-hypoxia and hyper-hypocapnia, can alter global CMRO₂ by 5 to 20%.^{78 –82} However, further detailed studies are warranted, as there has not been a clear consensus regarding decreases, increases, and no significant changes in CMRO₂ responding to increased CO₂ level in blood.^71,83,84 Also, the CMRO₂ response to hyperventilation may depend on participant characteristics (e.g., healthy subjects vs. patient with neurologic disorders) and the experimental setup (e.g., duration of hyperventilation). Additionally, we found that an average of 30.5% reduction in CBF in CGM, which is similar to previously reported values.^21,22

The study has several limitations. First, the results from both QSM-OEF and QQ-OEF may be affected by the partial volume effect. QSM-OEF utilizes a low-resolution CBF both directly (by using CBF in the assumption regarding the CMRO₂ ratio between RS and HV) and indirectly (by estimating v from CBF) to estimate a higher-resolution OEF. Hence, inaccuracies in CBF driven by partial voluming can lead to inaccurate OEF estimation in the QSM-OEF approach. In QQ-OEF, the OEF estimation in large veins may be problematic, as qBOLD in QQ assumes numerous randomly oriented cylindrical blood vessels within a voxel, which may not be valid in large veins. Second, the hematocrit (Hct) was not measured and assumed to be consistent (Hct = 0.357) in this study. In QSM-OEF, the maximum CMRO₂ estimation error caused by intersubject variability of Hct was 2.4%. ¹⁹ Further studies are warranted to investigate potential error in OEF estimation due to Hct variability. Third, an improved ASL protocol which follows a recent white paper ⁸⁵ could potentially improve SNR in CBF and subsequently enhance v and OEF estimation in the QSM-OEF approach. Fourth, the qBOLD term in QQ assumes two source components within a voxel: iron in venous structure (deoxyhemoglobin) and a diffusively distributed source (medium), such as non-blood tissue protein or ferritin. However, the presence of other structured strong sources like myelin in white matter and structured ferritin in deep gray matter could lead to intravoxel field variation, causing greater MRI signal decay than deoxyhemoglobin effects alone. This may result in an overestimation of QQ-OEF in WM, ³⁰ though this issue was observed to be not significant in our previous study. ³⁰

In conclusion, our study demonstrated QQ’s sensitivity to physiologically increased OEF levels in healthy subjects, in comparison to QSM-OEF. The non-invasive, readily applicable, challenge-free QQ-OEF mapping is now well-positioned for further evaluation in patients exhibiting increased OEF due to pathophysiological abnormalities.

Supplemental Material

sj-pdf-1-jcb-10.1177_0271678X241298584 - Supplemental material for Sensitivity assessment of QSM+qBOLD (or QQ) in detecting elevated oxygen extraction fraction (OEF) in physiological change

Supplemental material, sj-pdf-1-jcb-10.1177_0271678X241298584 for Sensitivity assessment of QSM+qBOLD (or QQ) in detecting elevated oxygen extraction fraction (OEF) in physiological change by Praveena Elanghovan, Thanh Nguyen, Pascal Spincemaille, Ajay Gupta, Yi Wang and Junghun Cho in Journal of Cerebral Blood Flow & Metabolism