Fast measurement of blood T1 in the human jugular vein at 3 Tesla

Recommendations for quantitative cerebral perfusion MRI using multi-timepoint arterial spin labeling: Acquisition, quantification, and clinical applications

Accurate assessment of cerebral perfusion is vital for understanding the hemodynamic processes involved in various neurological disorders and guiding clinical decision-making. This guidelines article provides a comprehensive overview of quantitative perfusion imaging of the brain using multi-timepoint arterial spin labeling (ASL), along with recommendations for its acquisition and quantification. A major benefit of acquiring ASL data with multiple label durations and/or post-labeling delays (PLDs) is being able to account for the effect of variable arterial transit time (ATT) on quantitative perfusion values and additionally visualize the spatial pattern of ATT itself, providing valuable clinical insights. Although multi-timepoint data can be acquired in the same scan time as single-PLD data with comparable perfusion measurement precision, its acquisition and postprocessing presents challenges beyond single-PLD ASL, impeding widespread adoption. Building upon the 2015 ASL consensus article, this work highlights the protocol distinctions specific to multi-timepoint ASL and provides robust recommendations for acquiring high-quality data. Additionally, we propose an extended quantification model based on the 2015 consensus model and discuss relevant postprocessing options to enhance the analysis of multi-timepoint ASL data. Furthermore, we review the potential clinical applications where multi-timepoint ASL is expected to offer significant benefits. This article is part of a series published by the International Society for Magnetic Resonance in Medicine (ISMRM) Perfusion Study Group, aiming to guide and inspire the advancement and utilization of ASL beyond the scope of the 2015 consensus article.

VICTR: Venous transit time imaging by changes in T1 relaxation

PURPOSE

Abnormalities in cerebral veins are a common finding in many neurological diseases, yet there is a scarcity of MRI techniques to assess venous hemodynamic function. The present study aims to develop a noncontrast technique to measure a novel blood flow circulatory measure, venous transit time (VTT), which denotes the time it takes for water to travel from capillary to major veins.

METHODS

The proposed sequence, venous transit time imaging by changes in T1 relaxation (VICTR), is based on the notion that as water molecules transition from the tissue into the veins, they undergo a change in T1 relaxation time. The validity of the measured VTT was tested by studying the VTT along the anatomically known flow trajectory of venous vessels as well as using a physiological vasoconstrictive challenge of caffeine ingestion. Finally, we compared the VTT measured with VICTR MRI to a bolus-tracking method using gadolinium-based contrast agent.

RESULTS

VTT was measured to be 3116.3 ± 326.0 ms in the posterior superior sagittal sinus (SSS), which was significantly longer than 2865.0 ± 390.8 ms at the anterior superior sagittal sinus (p = 0.004). The test-retest assessment showed an interclass correlation coefficient of 0.964. VTT was significantly increased by 513.8 ± 239.3 ms after caffeine ingestion (p < 0.001). VTT measured with VICTR MRI revealed a strong correlation (R = 0.84, p = 0.002) with that measured with the contrast-based approach. VTT was found inversely correlated to cerebral blood flow and venous oxygenation across individuals.

CONCLUSION

A noncontrast MRI technique, VICTR MRI, was developed to measure the VTT of the brain.

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.

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.

Comparative Analysis of Blood T2 Values Measured by T2 -TRIR and TRUST

BACKGROUND

Venous blood oxygenation (Yv), which can be derived from venous blood T₂ (T₂ b), combined with oxygen-extraction fraction (OEF) and cerebral metabolic rate of oxygen, is considered indicative for tissue viability and brain functioning and frequently assessed in patients with sickle cell disease. Recently, T₂ -Prepared-Blood-Relaxation-Imaging-with-Inversion-Recovery (T₂ -TRIR) was introduced allowing for simultaneous measurements of blood T₂ and T₁ (T₁ b), potentially improving Yv estimation by overcoming the need to estimate hematocrit.

PURPOSE

To optimize and compare T₂ -TRIR with T₂ -relaxation-under-spin-tagging (TRUST) sequence.

STUDY TYPE

Prospective.

POPULATION

A total of 12 healthy volunteers (six female, 27 ± 3 years old) and 7 patients with sickle cell disease (five female, 32 ± 12 years old).

FIELD STRENGTH/SEQUENCE

3 T; turbo field echo planar imaging (TFEPI), echo planar imaging (EPI), and fast field echo (FFE).

ASSESSMENT

T₂ b, Yv, and OEF from TRUST and T₂ -TRIR were compared and T₂ -TRIR-derived T₁ b was assessed. Within- and between-session repeatability was quantified in the controls, whereas sensitivity to hemodynamic changes after acetazolamide (ACZ) administration was assessed in the patients.

STATISTICAL TESTS

Shapiro-Wilk, one-sample and paired-sample t-test, repeated measures ANOVA, mixed linear model, Bland-Altman analysis and correlation analysis. Sidak multiple-comparison correction was performed. Significance level was 0.05.

RESULTS

In controls, T₂ b from T₂ -TRIR (70 ± 11 msec) was higher compared to TRUST (60 ± 8 msec). In patients, T₂ b values were lower pre- compared to post-ACZ administration (TRUST: 80 ± 15 msec and 106 ± 23 msec and T₂ -TRIR: 95 ± 21 msec and 125 ± 36 msec). Consequently, Yv and OEF were lower and higher pre- compared to post-ACZ administration (TRUST Yv: 68% ± 7% and 77% ± 8%, T₂ -TRIR Yv: 74% ± 8% and 80% ± 6%, TRUST OEF: 30% ± 7% and 21% ± 8%, and T₂ -TRIR OEF: 25% ± 8% and 18% ± 6%).

DATA CONCLUSION

TRUST and T₂ -TRIR are reproducible, but T₂ -TRIR-derived T₂ b values are significantly higher compared to TRUST, resulting in higher Yv and lower OEF estimates. This bias might be considered when evaluating cerebral oxygen homeostasis.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 2.

Effect of hematocrit on cerebral blood flow measured by pseudo-continuous arterial spin labeling MRI: A comparative study with 15O-water positron emission tomography

INTRODUCTION

In cerebral blood flow (CBF) quantification with pseudo-continuous arterial spin labeling (pCASL) MRI, arterial blood T1 (T1a) is usually fixed to a typical value (e.g., 1650 ms). However, individual T1a depends strongly on hematocrit (Hct) level. To investigate the utility of Hct-based T1a as an alternative to the fixed T1a method, we performed a comparative study with ¹⁵O-water positron emission tomography (PET).

METHODS

For patients with unilateral occlusion or stenosis of major arteries, hemispheric CBF on the healthy side was measured using pCASL and ¹⁵O-water PET. The pCASL CBFs were calculated with both (a) fixed T1a (1650 ms) and (b) individual T1a estimated from blood-sampled Hct (Hct-based T1a). Correlation coefficients of Hct-CBF were calculated and compared between pCASL and PET.

RESULTS

In pCASL, CBF with fixed T1a showed a strong negative correlation with Hct (r = -0.568), which was reduced with individual Hct-based T1a (r = -0.341 to -0.190), consistent with the Hct-CBF relation measured with PET (r = -0.349).

DISCUSSION AND CONCLUSION

We demonstrated that Hct-based T1a resulted in smaller inter-individual variations in pCASL CBF and an inverse Hct-CBF relationship more similar to that of PET. Care must be taken in the interpretation of pCASL CBF imaging in relation to Hct level even in subjects without anemia. Further comparative studies are needed to investigate whether advanced techniques improve pCASL CBF quantification at the individual level.

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.

Ensuring both velocity and spatial responses robust to B0/B1+ field inhomogeneities for velocity-selective arterial spin labeling through dynamic phase-cycling

PURPOSE

To evaluate both velocity and spatial responses of velocity-selective arterial spin labeling (VS-ASL), using velocity-insensitive and velocity-compensated waveforms for control modules, as well as a novel dynamic phase-cycling approach, at different B0 / B 1 + field inhomogeneities.

METHODS

In the presence of imperfect refocusing, the mechanism of phase-cycling the refocusing pulses through four dynamics was first theoretically analyzed with the conventional velocity-selective saturation (VSS) pulse train. Numerical simulations were then deployed to compare the performance of the Fourier-transform based velocity-selective inversion (FT-VSI) with these three different schemes in terms of both velocity and spatial responses under various B0 / B 1 + conditions. Phantom and human brain scans were performed to evaluate the three methods at B 1 + scales of 0.8, 1.0, and 1.2.

RESULTS

The simulations of FT-VSI showed that, under nonuniform B0 / B 1 + conditions, the scheme with velocity-insensitive control was susceptible to DC bias of the static spins as systematic error, while the scheme with velocity-compensated control had deteriorated velocity-selective labeling profiles and, thus, reduced labeling efficiency. Through numerical simulation, phantom scans, and brain perfusion measurements, the dynamic phase-cycling method demonstrated considerable improvements over these issues.

CONCLUSION

The proposed dynamic phase-cycling approach was demonstrated for the velocity-selective label and control modules with both velocity and spatial responses robust to a wide range of B0 and B 1 + field inhomogeneities.

Developmental patterns of CBF and BOLD responses to visual stimulus

To investigate the developmental changes of cerebral blood flow (CBF) and hemodynamic responses to changing neural activity, we used the arterial spin label (ASL) technique to measure resting CBF and simultaneous CBF / blood-oxygen-level dependent (BOLD) signal changes during visual stimulation in 97 typically developing children and young adults (age 13.35 [6.02, 25.25] (median [min, max]) years old at the first time point). The longitudinal study protocol included three MRIs (2.7 ± 0.06 obtained), one year apart, for each participant. Mixed-effect linear and non-linear statistical models were used to analyze age effects on CBF and BOLD signals. Resting CBF decreased exponentially with age (p = 0.0001) throughout the brain, and developmental trajectories differed across brain lobes. The absolute CBF increase in visual cortex during stimulation was constant over the age range, but the fractional CBF change increased with age (p = 0.0001) and the fractional BOLD signal increased with age (p = 0.0001) correspondingly. These findings suggest that the apparent neural hemodynamic coupling in visual cortex does not change after age six years, but age-related BOLD signal changes continue through adolescence primarily due to the changes with age in resting CBF.

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.

Non-contrast-enhanced abdominal MRA at 3 T using velocity-selective pulse trains

PURPOSE

Most existing non-contrast-enhanced methods for abdominal MR arteriography rely on a spatially selective inversion (SSI) pulse with a delay to null both static tissue and venous blood, and are limited to small spatial coverage due to the sensitivity to slow arterial inflow. Velocity-selective inversion (VSI) based approach has been shown to preserve the arterial blood inside the imaging volume at 1.5 T. Recently, velocity-selective saturation (VSS) pulse trains were applied to suppress the static tissue and have been combined with SSI pulses for cerebral MR arteriography at 3 T. The aim of this study is to construct an abdominal MRA protocol with large spatial coverage at 3 T using advanced velocity-selective pulse trains.

METHODS

Multiple velocity-selective MRA protocols with different sequence modules and 3D acquisition methods were evaluated. Sequences using VSS only as well as SSI+VSS and VSI+VSS preparations were then compared among a group of healthy young and middle-aged volunteers. Using MRA without any preparations as reference, relative signal ratios and relative contrast ratios of different vascular segments were quantitatively analyzed.

RESULTS

Both SSI+VSS and VSI+VSS arteriograms achieved high artery-to-tissue and artery-to-vein relative contrast ratios above aortic bifurcation. The SSI+VSS sequence yielded lower signal at the bilateral iliac arteries than VSI+VSS, reflecting the benefit of the VSI preparation for imaging the distal branches.

CONCLUSION

The feasibility of noncontrast 3D MR abdominal arteriography was demonstrated on healthy volunteers using a combination of VSS pulse trains and SSI or VSI pulse.

Improved velocity-selective-inversion arterial spin labeling for cerebral blood flow mapping with 3D acquisition

PURPOSE

To further optimize the velocity-selective arterial spin labeling (VSASL) sequence utilizing a Fourier-transform based velocity-selective inversion (FT-VSI) pulse train, and to evaluate its utility for 3D mapping of cerebral blood flow (CBF) with a gradient- and spin-echo (GRASE) readout.

METHODS

First, numerical simulations and phantom experiments were done to test the susceptibility to eddy currents and B₁ field inhomogeneities for FT-VSI pulse trains with block and composite refocusing pulses. Second, the choices of the post-labeling delay (PLD) for FT-VSI prepared 3D VSASL were evaluated for the sensitivity to perfusion signal. The study was conducted among a young-age and a middle-age group at 3T. Both signal-to-noise ratio (SNR) and CBF were quantitatively compared with pseudo-continuous ASL (PCASL). The optimized 3D VSI-ASL was also qualitatively compared with PCASL in a whole-brain coverage among two healthy volunteers and a brain tumor patient.

RESULTS

The simulations and phantom test showed that composite refocusing pulses are more robust to both eddy-currents and B₁ field inhomogeneities than block pulses. 3D VSASL images with FT-VSI preparation were acquired over a range of PLDs and PLD = 1.2 s was selected for its higher perfusion signal. FT-VSI labeling produced quantitative CBF maps with 27% higher SNR in gray matter compared to PCASL. 3D whole-brain CBF mapping using VSI-ASL were comparable to the corresponding PCASL results.

CONCLUSION

FT-VSI with 3D-GRASE readout was successfully implemented and showed higher sensitivity to perfusion signal than PCASL for both young and middle-aged healthy volunteers.

Towards measuring the effect of flow in blood T 1 assessed in a flow phantom and in vivo

Measurement of the blood T ₁ time using conventional myocardial T ₁ mapping methods has gained clinical significance in the context of extracellular volume (ECV) mapping and synthetic hematocrit (Hct). However, its accuracy is potentially compromised by in-flow of non-inverted/non-saturated spins and in-flow of spins which are not partially saturated from previous imaging pulses. Bloch simulations were used to analyze various flow effects separately. T ₁ measurements of gadolinium doped water were performed using a flow phantom with adjustable flow velocities at 3 T. Additionally, in vivo blood T ₁ measurements were performed in 6 healthy subjects (26 ± 5 years, 2 female). To study the T ₁ time as a function of the instantaneous flow velocity, T ₁ times were evaluated in an axial imaging slice of the descending aorta. Velocity encoded cine measurements were performed to quantify the flow velocity throughout the cardiac cycle. Simulation results show more than 30% loss in accuracy for 10% non-prepared in-flowing spins. However, in- and out-flow to the imaging plane only demonstrated minor impact on the T ₁ time. Phantom T ₁ times were decreased by up to 200 ms in the flow phantom, due to in-flow of non-prepared spins. High flow velocities cause in-flow of spins that lack partial saturation from the imaging pulses but only lead to negligible T ₁ time deviation (less than 30 ms). In vivo measurements confirm a substantial variation of the T ₁ time depending on the flow velocity. The highest aortic T ₁ times are observed at the time point of minimal flow with increased flow velocity leading to reduction of the measured T ₁ time by up to [Formula: see text] at peak velocity. In this work we attempt to dissect the effects of flow on T ₁ times, by using simulations, well-controlled, simplified phantom setup and the linear flow pattern in the descending aorta in vivo.

Cerebral blood volume mapping using Fourier-transform-based velocity-selective saturation pulse trains

PURPOSE

Velocity-selective saturation (VSS) pulse trains provide a viable alternative to the spatially selective methods for measuring cerebral blood volume (CBV) by reducing the sensitivity to arterial transit time. This study is to compare the Fourier-transform-based velocity-selective saturation (FT-VSS) pulse trains with the conventional flow-dephasing VSS techniques for CBV quantification.

METHODS

The proposed FT-VSS label and control modules were compared with VSS pulse trains utilizing double refocused hyperbolic tangent (DRHT) and 8-segment B1-insensitive rotation (BIR-8). This was done using both numerical simulations and phantom studies to evaluate their sensitivities to gradient imperfections such as eddy currents. DRHT, BIR-8, and FT-VSS prepared CBV mapping was further compared for velocity-encoding gradients along 3 orthogonal directions in healthy subjects at 3T.

RESULTS

The phantom studies exhibited more consistent immunity to gradient imperfections for the utilized FT-VSS pulse trains. Compared to DRHT and BIR-8, FT-VSS delivered more robust CBV results across the 3 VS encoding directions with significantly reduced artifacts along the superior-inferior direction and improved temporal signal-to-noise ratio (SNR) values. Average CBV values obtained from FT-VSS based sequences were 5.3 mL/100 g for gray matter and 2.3 mL/100 g for white matter, comparable to literature expectations.

CONCLUSION

Absolute CBV quantification utilizing advanced FT-VSS pulse trains had several advantages over the existing approaches using flow-dephasing VSS modules. A greater immunity to gradient imperfections and the concurrent tissue background suppression of FT-VSS pulse trains enabled more robust CBV measurements and higher SNR than the conventional VSS pulse trains.

Volume-localized measurement of oxygen extraction fraction in the brain using MRI

Purpose

T₂‐relaxation‐under‐spin‐tagging (TRUST) is an MR technique for the non‐invasive assessment of whole‐brain cerebral oxygen extraction fraction (OEF), through measurement of the venous blood T₂ relaxation time in the sagittal sinus. A key limitation of TRUST, however, is the lack of spatial specificity of the measurement. We sought to develop a modified TRUST sequence, selective localized TRUST (SL‐TRUST), having sensitivity to venous blood T₂ within a targeted brain region, and therefore achieving spatially localized measurements of cerebral tissue OEF, while still retaining acquisition in the sagittal sinus.

Methods

A method for selective localization of TRUST sequence was developed, and the reproducibility of the technique was evaluated in healthy participants. Regional measurements were achieved for a single hemisphere and for a 3D‐localized 70 × 70 × 80 mm³ tissue region using SL‐TRUST and compared to a global TRUST measure. An additional measure of venous blood T₁ in the sagittal sinus was used to estimate subject‐specific hematocrit. Six subjects were scanned over 4 sessions, including intra‐session repeat measurements.

Results

The average T₂ in the sagittal sinus was found to be 60.8 ± 8.9, 62.7 ± 7.9, 64.6 ± 8.4, and 66.3 ± 10.3 ms (mean ± SD) for conventional TRUST, global SL‐TRUST, hemispheric SL‐TRUST, and 3D‐localized SL‐TRUST, respectively. Intra‐, inter‐session, and inter‐subject coefficients of variation for OEF using SL‐TRUST were found to be comparable and in some cases superior to those obtained using TRUST.

Conclusion

OEF comparison of 2 contralateral regions was achievable in under 5 min suggesting SL‐TRUST offers potential for quantifying regional OEF differences in both healthy and clinical populations.

Precision of T1-relaxation time measurements in the hepatic portal vein: influence of measurement technique and sequence parameters

OBJECTIVE

To investigate the effects of a range of parameter settings on T₁ measurement stability in the portal vein using the T1-mapping sequences Look-Locker (LL) and Modified Look-Locker inversion recovery (MOLLI).

MATERIALS AND METHODS

Ten different versions of LL and MOLLI sequences were tested and compared to a reference sequence provided by the MR manufacturer. Ten healthy volunteers were imaged multiple times on two separate scan days at 3T. The mean T₁ values and coefficient of variation (CoV) were calculated for each of the ten sequences and compared to the reference sequence.

RESULTS

Six of the tested sequences had T₁ values close to the reference sequence; among those, three sequences achieved lower CoV than the reference sequence. Lowest CoV was achieved using a non-triggered LL sequence with 5 beat readout and a 45^o flip angle (mean T₁ 1733 ms ± 89 ms, CoV 1.3% ± 0.58%).

CONCLUSION

T₁-measurements in the hepatic portal vein can be performed with high precision using either MOLLI or LL sequences provided that LL sampling duration is sufficiently long and flip angle sufficiently high. The advantage of constant timing outweighed the advantage of ECG-triggering.

Cerebral perfusion characteristics show differences in younger versus older children with sickle cell anaemia: Results from a multiple-inflow-time arterial spin labelling study

Sickle cell anaemia (SCA) is associated with chronic anaemia and oxygen desaturation, which elevate cerebral blood flow (CBF) and increase the risk of neurocognitive complications. Arterial spin labelling (ASL) provides a methodology for measuring CBF non-invasively; however, ASL techniques using only a single inflow time are not sufficient to fully characterize abnormal haemodynamic behaviour in SCA. This study investigated haemodynamic parameters from a multi-inflow-time ASL acquisition in younger (8-12 years) and older (13-18 years) children with SCA with and without silent cerebral infarction (SCI+/-) (n = 20 and 19 respectively, 6 and 4 SCI+ respectively) and healthy controls (n = 9 and 7 respectively). Compared with controls, CBF was elevated globally in both groups of patients. In the younger SCA patients, blood oxygen content was negatively correlated with CBF in the middle and posterior cerebral artery territories and significantly positively correlated with bolus arrival time (BAT) in the anterior and middle cerebral artery territories. In older children, SCA patients had significantly shorter BAT than healthy controls and there was a significant negative correlation between CBF and oxygen content only in the territory of the posterior cerebral artery, with a trend for a correlation in the anterior cerebral artery but no relationship for the middle cerebral artery territory. In the younger group, SCI+ patients had significantly higher CBF in the posterior cerebral artery territory (SCI+ mean = 92.78 ml/100 g/min; SCI- mean = 72.71 ml/100 g/min; F = 4.28, p = 0.04), but this no longer reached significance when two children with abnormal transcranial Doppler and one with haemoglobin SC disease were excluded, and there were no significant differences between patients with and without SCI in the older children. With age, there appears to be increasing disparity between patients and controls in terms of the relationship between CBF and oxygen content in the anterior circulation, potentially predicting the risk of acute and chronic compromise of brain tissue.

Accounting for the role of hematocrit in between-subject variations of MRI-derived baseline cerebral hemodynamic parameters and functional BOLD responses

Baseline hematocrit fraction (Hct) is a determinant for baseline cerebral blood flow (CBF) and between-subject variation of Hct thus causes variation in task-based BOLD fMRI signal changes. We first verified in healthy volunteers (n = 12) that Hct values can be derived reliably from venous blood T₁ values by comparison with the conventional lab test. Together with CBF measured using phase-contrast MRI, this noninvasive estimation of Hct, instead of using a population-averaged Hct value, enabled more individual determination of oxygen delivery (DO₂ ), oxygen extraction fraction (OEF), and cerebral metabolic rate of oxygen (CMRO₂ ). The inverse correlation of CBF and Hct explained about 80% of between-subject variation of CBF in this relatively uniform cohort of subjects, as expected based on the regulation of DO₂ to maintain constant CMRO₂ . Furthermore, we compared the relationships of visual task-evoked BOLD response with Hct and CBF. We showed that Hct and CBF contributed 22%-33% of variance in BOLD signal and removing the positive correlation with Hct and negative correlation with CBF allowed normalization of BOLD signal with 16%-22% lower variability. The results of this study suggest that adjustment for Hct effects is useful for studies of MRI perfusion and BOLD fMRI. Hum Brain Mapp 39:344-353, 2018. © 2017 Wiley Periodicals, Inc.

Whole-brain arteriography and venography: Using improved velocity-selective saturation pulse trains

PURPOSE

To develop velocity-selective (VS) MR angiography (MRA) protocols for arteriography and venography with whole-brain coverage.

METHODS

Tissue suppression using velocity-selective saturation (VSS) pulse trains is sensitive to radiofrequency field (B₁ +) inhomogeneity. To reduce its sensitivity, we replaced the low-flip-angle hard pulses in the VSS pulse train with optimal composite (OCP) pulses. Additionally, new pulse sequences for arteriography and venography were developed by placing spatially selective inversion pulses with a delay to null signals from either venous or arterial blood. The VS MRA techniques were compared to the time-of-flight (TOF) MRA in six healthy subjects and two patients at 3T.

RESULTS

More uniform suppression of stationary tissue was observed when the hard pulses were replaced by OCP pulses in the VSS pulse trains, which improved contrast ratios between blood vessels and tissue background for both arteries (0.87 vs. 0.77) and veins (0.80 vs. 0.59). Both arteriograms and venograms depicted all major cervical and intracranial arteries and veins, respectively. Compared to TOF MRA, VS MRA not only offers larger spatial coverage but also depicts more small vessels. Initial clinical feasibility was shown in two patients with comparisons to TOF protocols.

CONCLUSION

Noncontrast-enhanced whole-brain arteriography and venography can be obtained without losing sensitivity to small vessel detection. Magn Reson Med 79:2014-2023, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Fast measurement of blood T1 in the human carotid artery at 3T: Accuracy, precision, and reproducibility

PURPOSE

To develop a fast protocol for measuring T1 values in the internal carotid artery (ICA), to validate this technique with in vitro measurements, and to evaluate its reproducibility.

METHODS

A modified Look-Locker sequence was optimized to enable rapid determination of T1 in the ICA at 3T. T1 values from the ICA were compared with in vitro measurements on individually sampled venous blood oxygenated to arterial levels. A test-retest reproducibility study was also conducted.

RESULTS

The group-averaged arterial blood T1 value was 1908 ± 77 ms for six women (hematocrit = 0.39 ± 0.03) and 1785 ± 55 ms for seven men (hematocrit = 0.45 ± 0.02), which is 100-200 ms longer than the widely adopted value obtained from bovine blood experiments. The arterial T1 value per subject correlated significantly with individual hematocrit values. The intrasession and intersession coefficients of variation were 1.1% and 2.1%, respectively, indicating good precision and reproducibility of our method. Reasonable agreement was observed between the in vivo and in vitro results with a correlation coefficient of 0.78.

CONCLUSION

The proposed method can provide fast arterial T1 measurement on individual subjects. When not performing such a subject-specific measurement, we recommend the use of 1908 ms and 1785 ms for healthy women and men, respectively, or 1841 ms for adults in general. Magn Reson Med 77:2296-2302, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

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.

Velocity-selective-inversion prepared arterial spin labeling

PURPOSE

To develop a Fourier-transform based velocity-selective inversion (FT-VSI) pulse train for velocity-selective arterial spin labeling (VSASL).

METHODS

This new pulse contains paired and phase cycled refocusing pulses. Its sensitivities to B0/B1 inhomogeneity and gradient imperfections such as eddy currents were evaluated through simulation and phantom studies. Cerebral blood flow (CBF) quantification using FT-VSI prepared VSASL was compared with conventional VSASL and pseudocontinuous ASL (PCASL) at 3 Tesla.

RESULTS

Simulation and phantom results of the proposed FT-VSI pulse train demonstrated excellent robustness to B0/B1 field inhomogeneity and eddy currents. The estimated CBF of gray matter and white matter for the FT-VSI prepared VSASL, averaged among eight healthy volunteers, were 49.5 ± 7.5 mL/100 g/min and 14.8 ± 2.4 mL/100 g/min, respectively. Excellent correlation and agreement between the FT-VSI method and conventional VSASL and PCASL were found. The averaged signal-to-noise ratio (SNR) value in gray matter of the FT-VSI method was 39% higher than VSASL using conventional double refocused hyperbolic tangent pulses and 9% lower than PCASL.

CONCLUSION

A novel FT-VSI pulse train was demonstrated to be a suitable labeling module for VSASL with robustness of velocity-selective profile to B0/B1 field inhomogeneity and gradient imperfections. Compared with conventional VSASL, FT-VSI prepared VSASL produced consistent CBF maps with higher SNR values. Magn Reson Med 76:1136-1148, 2016. © 2015 Wiley Periodicals, Inc.

Accuracy of Parenchymal Cerebral Blood Flow Measurements Using Pseudocontinuous Arterial Spin-Labeling in Healthy Volunteers

BACKGROUND AND PURPOSE

The arterial spin-labeling method for CBF assessment is widely available, but its accuracy is not fully established. We investigated the accuracy of a whole-brain arterial spin-labeling technique for assessing the mean parenchymal CBF and the effect of aging in healthy volunteers. Phase-contrast MR imaging was used as the reference method.

MATERIALS AND METHODS

Ninety-two healthy volunteers were included: 49 young (age range, 20-30 years) and 43 elderly (age range, 65-80 years). Arterial spin-labeling parenchymal CBF values were averaged over the whole brain to quantify the mean pCBF(ASL) value. Total CBF was assessed with phase-contrast MR imaging as the sum of flows in the internal carotid and vertebral arteries, and subsequent division by brain volume returned the pCBF(PCMRI) value. Accuracy was considered as good as that of the reference method if the systematic difference was less than 5 mL/min/100 g of brain tissue and if the 95% confidence intervals were equal to or better than ±10 mL/min/100 g.

RESULTS

pCBF(ASL) correlated to pCBF(PCMRI) (r = 0.73; P < .001). Significant differences were observed between the pCBF(ASL) and pCBF(PCMRI) values in the young (P = .001) and the elderly (P < .001) volunteers. The systematic differences (mean ± 2 standard deviations) were -4 ± 14 mL/min/100 g in the young subjects and 6 ± 12 mL/min/100 g in the elderly subjects. Young subjects showed higher values than the elderly subjects for pCBF(PCMRI) (young, 57 ± 8 mL/min/100 g; elderly, 54 ± 7 mL/min/100 g; P = .05) and pCBF(ASL) (young, 61 ± 10 mL/min/100 g; elderly, 48 ± 10 mL/min/100 g; P < .001).

CONCLUSIONS

The limits of agreement were too wide for the arterial spin-labeling method to be considered satisfactorily accurate, whereas the systematic overestimation in the young subjects and underestimation in the elderly subjects were close to acceptable. The age-related decrease in parenchymal CBF was augmented in arterial spin-labeling compared with phase-contrast MR imaging.

Quantitative theory for the longitudinal relaxation time of blood water

PURPOSE

To propose and evaluate a model for the blood water T1 that takes into account the effects of hematocrit fraction, oxygenation fraction, erythrocyte hemoglobin concentration, methemoglobin fraction, and plasma albumin concentration.

METHODS

Whole blood and lysed blood T1 data were acquired at magnetic fields of 3 Tesla (T), 7T, 9.4T, and 11.7T using inversion-recovery measurements and a home-built blood circulation system for maintaining physiological conditions. A quantitative model was derived based on multivariable fitting of this data.

RESULTS

Fitting of the model to the data allowed determination of the different parameters describing the blood water T1 such as those for the diamagnetic and paramagnetic effects of albumin and hemoglobin, and the contribution of methemoglobin. The model correctly predicts blood T1 at multiple fields, as verified by comparison with existing literature.

CONCLUSION

The model provides physical and physiological parameters describing the effects of hematocrit fraction, oxygenation, hemoglobin concentration, methemoglobin fraction, and albumin concentration on blood water T1 . It can be used to predict blood T1 at multiple fields. Magn Reson Med 76:270-281, 2016. © 2015 Wiley Periodicals, Inc.

Aortic 4D flow: quantification of signal-to-noise ratio as a function of field strength and contrast enhancement for 1.5T, 3T, and 7T

PURPOSE

To investigate for the first time the feasibility of aortic four-dimensional (4D) flow at 7T, both contrast enhanced (CE) and non-CE. To quantify the signal-to-noise ratio (SNR) in aortic 4D flow as a function of field strength and CE with gadobenate dimeglumine (MultiHance).

METHODS

Six healthy male volunteers were scanned at 1.5T, 3T, and 7T with both non-CE and CE acquisitions. Temporal SNR was calculated. Flip angle optimization for CE 4D flow was carried out using Bloch simulations that were validated against in vivo measurements.

RESULTS

The 7T provided 2.2 times the SNR of 3T while 3T provided 1.7 times the SNR of 1.5T in non-CE acquisitions in the descending aorta. The SNR gains achieved by CE were 1.8-fold at 1.5T, 1.7-fold at 3T, and 1.4-fold at 7T, respectively.

CONCLUSION

The 7T provides a new tool to explore aortic 4D flow, yielding higher SNR that can be used to push the boundaries of acceleration and resolution. Field strength and contrast enhancement at all fields provide significant improvements in SNR.

Impact of neonate haematocrit variability on the longitudinal relaxation time of blood: Implications for arterial spin labelling MRI

BACKGROUND AND PURPOSE

The longitudinal relaxation time of blood (T 1b) is influenced by haematocrit (Hct) which is known to vary in neonates. The purpose of this study was threefold: to obtain T 1b values in neonates, to investigate how the T 1b influences quantitative arterial spin labelling (ASL), and to evaluate if known relationships between T 1b and haematocrit (Hct) hold true when Hct is measured by means of a point-of-care device.

MATERIALS AND METHODS

One hundred and four neonates with 120 MR scan sessions (3 T) were included. The T 1b was obtained from a T 1 inversion recovery sequence. T 1b-induced changes in ASL cerebral blood flow estimates were evaluated. The Hct was obtained by means of a point-of-care device. Linear regression analysis was used to investigate the relation between Hct and MRI-derived R1 of blood (the inverse of the T 1b).

RESULTS

Mean T 1b was 1.85 s (sd 0.2 s). The mean T 1b in preterm neonates was 1.77 s, 1.89 s in preterm neonates scanned at term-equivalent age (TEA) and 1.81 s in diseased neonates. The T 1b in the TEA was significantly different from the T 1b in the preterm (p < 0.05). The change in perfusion induced by the T 1b was -11% (sd 9.1%, p < 0.001). The relation between arterial-drawn Hct and R1b was R1b = 0.80 × Hct + 0.22, which falls within the confidence interval of the previously established relationships, whereas capillary-drawn Hct did not correlate with R1b.

CONCLUSION

We demonstrated a wide variability of the T 1b in neonates and the implications it could have in methods relying on the actual T 1b as for instance ASL. It was concluded that arterial-drawn Hct values obtained from a point-of-care device can be used to infer the T 1b whereas our data did not support the use of capillary-drawn Hct for T 1b correction.

Three-dimensional whole-brain perfusion quantification using pseudo-continuous arterial spin labeling MRI at multiple post-labeling delays: accounting for both arterial transit time and impulse response function

Measurement of the cerebral blood flow (CBF) with whole-brain coverage is challenging in terms of both acquisition and quantitative analysis. In order to fit arterial spin labeling-based perfusion kinetic curves, an empirical three-parameter model which characterizes the effective impulse response function (IRF) is introduced, which allows the determination of CBF, the arterial transit time (ATT) and T(1,eff). The accuracy and precision of the proposed model were compared with those of more complicated models with four or five parameters through Monte Carlo simulations. Pseudo-continuous arterial spin labeling images were acquired on a clinical 3-T scanner in 10 normal volunteers using a three-dimensional multi-shot gradient and spin echo scheme at multiple post-labeling delays to sample the kinetic curves. Voxel-wise fitting was performed using the three-parameter model and other models that contain two, four or five unknown parameters. For the two-parameter model, T(1,eff) values close to tissue and blood were assumed separately. Standard statistical analysis was conducted to compare these fitting models in various brain regions. The fitted results indicated that: (i) the estimated CBF values using the two-parameter model show appreciable dependence on the assumed T(1,eff) values; (ii) the proposed three-parameter model achieves the optimal balance between the goodness of fit and model complexity when compared among the models with explicit IRF fitting; (iii) both the two-parameter model using fixed blood T1 values for T(1,eff) and the three-parameter model provide reasonable fitting results. Using the proposed three-parameter model, the estimated CBF (46 ± 14 mL/100 g/min) and ATT (1.4 ± 0.3 s) values averaged from different brain regions are close to the literature reports; the estimated T(1,eff) values (1.9 ± 0.4 s) are higher than the tissue T1 values, possibly reflecting a contribution from the microvascular arterial blood compartment.

Noninvasive functional imaging of cerebral blood volume with vascular-space-occupancy (VASO) MRI

Functional MRI (fMRI) based on changes in cerebral blood volume (CBV) can probe directly vasodilatation and vasoconstriction during brain activation or physiologic challenges, and can provide important insights into the mechanism of blood oxygenation level-dependent (BOLD) signal changes. At present, the most widely used CBV fMRI technique in humans is called vascular-space-occupancy (VASO) MRI, and this article provides a technical review of this method. VASO MRI utilizes T1 differences between blood and tissue to distinguish between these two compartments within a voxel, and employs a blood-nulling inversion recovery sequence to yield an MR signal proportional to 1 - CBV. As such, vasodilatation will result in a VASO signal decrease and vasoconstriction will have the reverse effect. The VASO technique can be performed dynamically with a temporal resolution comparable with several other fMRI methods, such as BOLD or arterial spin labeling (ASL), and is particularly powerful when conducted in conjunction with these complementary techniques. The pulse sequence and imaging parameters of VASO can be optimized such that the signal change is predominantly of CBV origin, but careful considerations should be taken to minimize other contributions, such as those from the BOLD effect, cerebral blood flow (CBF) and cerebrospinal fluid (CSF). The sensitivity of the VASO technique is the primary disadvantage when compared with BOLD, but this technique is increasingly demonstrating its utility in neuroscientific and clinical applications.

The use of dynamic tracer concentration in veins for quantitative DCE-MRI kinetic analysis in head and neck

BACKGROUND

Head and neck Magnetic Resonance (MR) Images are vulnerable to the arterial blood in-flow effect. To compensate for this effect and enhance accuracy and reproducibility, dynamic tracer concentration in veins was proposed and investigated for quantitative dynamic contrast-enhanced (DCE) MRI analysis in head and neck.

METHODOLOGY

21 patients with head and neck tumors underwent DCE-MRI at 3T. An automated method was developed for blood vessel selection and separation. Dynamic concentration-time-curves (CTCs) in arteries and veins were used for the Tofts model parameter estimations. The estimation differences by using CTCs in arteries and veins were compared. Artery and vein voxels were accurately separated by the automated method. Remarkable inter-slice tracer concentration differences were found in arteries while the inter-slice concentration differences in veins were moderate. Tofts model fitting by using the CTCs in arteries and veins produced significantly different parameter estimations. The individual artery CTCs resulted in large (>50% generally) inter-slice parameter estimation variations. Better inter-slice consistency was achieved by using the vein CTCs.

CONCLUSIONS

The use of vein CTCs helps to compensate for arterial in-flow effect and reduce kinetic parameter estimation error and inconsistency for head and neck DCE-MRI.

Hematocrit and oxygenation dependence of blood (1)H(2)O T(1) at 7 Tesla

Knowledge of blood (1)H2O T1 is critical for perfusion-based quantification experiments such as arterial spin labeling and cerebral blood volume-weighted MRI using vascular space occupancy. The dependence of blood (1)H2O T1 on hematocrit fraction (Hct) and oxygen saturation fraction (Y) was determined at 7 T using in vitro bovine blood in a circulating system under physiological conditions. Blood (1)H2O R1 values for different conditions could be readily fitted using a two-compartment (erythrocyte and plasma) model, which are described by a monoexponential longitudinal relaxation rate constant dependence. It was found that T1 = 2171 ± 39 ms for Y = 1 (arterial blood) and 2010 ± 41 ms for Y = 0.6 (venous blood), for a typical Hct of 0.42. The blood (1)H2O T1 values in the normal physiological range (Hct from 0.35 to 0.45, and Y from 0.6 to 1.0) were determined to range from 1900 to 2300 ms. The influence of oxygen partial pressure (pO2) and the effect of plasma osmolality for different anticoagulants were also investigated. It is discussed why blood (1)H2O T1 values measured in vivo for human blood may be about 10-20% larger than found in vitro for bovine blood at the same field strength.

Measurement of T1 of human arterial and venous blood at 7T

Techniques for measuring cerebral perfusion require accurate longitudinal relaxation (T1) of blood, an MRI parameter that is field dependent. T1 of arterial and venous human blood was measured at 7T using three different sources - pathology laboratory, blood bank and in vivo. The T1 of venous blood was measured from sealed samples from a pathology lab and in vivo. Samples from a blood bank were oxygenated and mixed to obtain different physiological concentrations of hematocrit and oxygenation. T1 relaxation times were estimated using a three-point fit to a simple inversion recovery equation. At 37°C, the T1 of blood at arterial pO2 was 2.29±0.1s and 2.07±0.12 at venous pO2. The in vivo T1 of venous blood, in three subjects, was slightly longer at 2.45±0.11s. T1 of arterial and venous blood at 7T was measured and found to be significantly different. The T1 values were longer in vivo than in vitro. While the exact cause for the discrepancy is unknown, the additives in the blood samples, degradation during experiment, oxygenation differences, and the non-stagnant nature of blood in vivo could be potential contributors to the lower values of T1 in the venous samples.

Noninvasive Measurements of Cerebral Blood Flow, Oxygen Extraction Fraction, and Oxygen Metabolic Index in Human with Inhalation of Air and Carbogen using Magnetic Resonance Imaging

Noninvasive magnetic resonance (MR) methods have been explored to provide quantitative measurements of cerebral blood flow (CBF), oxygen extraction fraction (OEF), and oxygen metabolic index (OMI = CBF × OEF). In this study, we sought to evaluate whether MR measured OEF, CBF, and OMI can consistently detect the expected physiological changes in humans under normal and hyperoxic hypercapnic conditions. Nine healthy human subjects were scanned while breathing through a mask, alternating inhaled gas in a sequential order as room air, carbogen (3% CO2 mixed with 97% O2), room air, carbogen, and room air. OEF, CBF, and OMI were obtained from the whole brain, gray matter (GM), and white matter (WM) at each gas inhalation state. Similar to previous positron emission tomography findings, our study consistently demonstrated a 10-12% decrease in OEF with a 10% increase of CBF and a stable OMI during carbogen inhalation. Moreover, GM/WM ratio in CBF and OMI remained constant during air and carbogen breathing. In addition, OEF, CBF, and OMI were highly reproducible if the same inhaled gas was used. In summary, our results demonstrate that noninvasive MR measurements can provide reproducible measurements of OEF, CBF, and OMI in normal subjects under normal and altered physiological conditions.