Time-resolved MRI oximetry for quantifying CMRO(2) and vascular reactivity

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.

Repeatability assessment for simultaneous measurement of arterial blood flow, venous oxygen saturation, and muscle perfusion following dynamic exercise

The purpose of the present study was to demonstrate a new sequence and determine the repeatability of simultaneous dynamic measurements of blood flow, venous oxygen saturation (SvO₂ ), and relative perfusion (change from resting perfusion) in calf muscle during recovery from plantar flexion exercise. The feasibility of near simultaneous measurement of bio-energetic parameters was also demonstrated. A sequence was developed to simultaneously measure arterial blood flow using flow-encoded projection, SvO₂ using susceptibility-based oximetry, and relative perfusion using arterial spin labeling in combination with dynamic plantar flexion exercise. The parameters were determined at rest and during recovery from single leg plantar flexion exercise. Test-retest repeatability was analyzed using Bland-Altman analysis and intraclass correlation coefficients (ICC). The mitochondrial capacity of skeletal muscle was also measured immediately afterwards with dynamic phosphorus magnetic resonance spectroscopy. Eight healthy subjects participated in the study for test-retest repeatability. Popliteal artery blood flow at rest was 1.79 ± 0.58 ml/s and increased to 11.18 ± 3.02 ml/s immediately after exercise. Popliteal vein SvO₂ decreased to 45.93% ± 6.5% from a resting value of 70.46% ± 4.76% following exercise. Relative perfusion (change from rest value) was 51.83 ± 15.00 ml/100 g/min at the cessation of exercise. The recovery of blood flow and SvO₂ was modeled as a single exponential with time constants of 38.03 ± 6.91 and 71.19 ± 14.53 s, respectively. All the measured parameters exhibited good repeatability with ICC ranging from 0.8 to 0.95. Bioenergetics measurements were within normal range, demonstrating the feasibility of near simultaneous measurement of hemodynamic and energetic parameters. Clinical feasibility was assessed with Barth syndrome patients, demonstrating reduced oxygen extraction from the blood and reduced mitochondrial oxidative capacity compared with healthy controls. The proposed protocol allows rapid imaging of multiple parameters in skeletal muscle that might be affected in disease.

Noninvasive Quantification of Cerebral Blood Flow Using Hybrid PET/MR Imaging to Extract the [15 O]H2 O Image-Derived Input Function Free of Partial Volume Errors

BACKGROUND

Quantification of cerebral blood flow (CBF) with [¹⁵ O]H₂ O-positron emission tomography (PET) requires arterial sampling to measure the input function. This invasive procedure can be avoided by extracting an image-derived input function (IDIF); however, IDIFs are sensitive to partial volume errors due to the limited spatial resolution of PET.

PURPOSE

To present an alternative hybrid PET/MR imaging of CBF (PMRFlow_IDIF ) that uses phase-contrast (PC) MRI measurements of whole-brain (WB) CBF to calibrate an IDIF extracted from a WB [¹⁵ O]H₂ O time-activity curve.

STUDY TYPE

Technical development and validation.

ANIMAL MODEL

Twelve juvenile Duroc pigs (83% female).

POPULATION

Thirteen healthy individuals (38% female).

FIELD STRENGTH/SEQUENCES

3 T; gradient-echo PC-MRI.

ASSESSMENT

PMRFlow_IDIF was validated against PET-only in a porcine model that included arterial sampling. CBF maps were generated by applying PMRFlow_IDIF and two previous PMRFlow methods (PC-PET and double integration method [DIM]) to [¹⁵ O]H₂ O-PET data acquired from healthy individuals.

STATISTICAL TESTS

PMRFlow and PET CBF measurements were compared with regression and correlation analyses. Paired t-tests were performed to evaluate differences. Potential biases were assessed using one-sample t-tests. Reliability was assessed by intraclass correlation coefficients. Statistical significance: α  = 0.05.

RESULTS

In the animal study, strong agreement was observed between PMRFlow_IDIF (average voxel-wise CBF, 58.0 ± 16.9 mL/100 g/min) and PET (63.0 ± 18.9 mL/100 g/min). In the human study, PMRFlow_DIM (y = 1.11x - 5.16, R²  = 0.99 ± 0.01) and PMRFlow_PC-PET (y = 0.87x + 3.82, R²  = 0.97 ± 0.02) performed similarly to PMRFlow_IDIF, and CBF was within the expected range (eg, 49.7 ± 7.2 mL/100 g/min for gray matter).

DATA CONCLUSION

Accuracy of PMRFlow_IDIF was confirmed in the animal study with the primary source of error attributed to differences in WB CBF measured by PC MRI and PET. In the human study, differences in CBF from PMRFlow_IDIF , PMRFlow_DIM , and PMRFlow_PC-PET were due to the latter two not accounting for blood-borne activity.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1.

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.

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.

Cerebral metabolic rate of oxygen during transition from wakefulness to sleep measured with high temporal resolution OxFlow MRI with concurrent EEG

During slow-wave sleep, synaptic transmissions are reduced with a concomitant reduction in brain energy consumption. We used 3 Tesla MRI to noninvasively quantify changes in the cerebral metabolic rate of O2 (CMRO2) during wakefulness and sleep, leveraging the 'OxFlow' method, which provides venous O2 saturation (SvO2) along with cerebral blood flow (CBF). Twelve healthy subjects (31.3 ± 5.6 years, eight males) underwent 45-60 min of continuous scanning during wakefulness and sleep, yielding one image set every 3.4 s. Concurrent electroencephalography (EEG) data were available in eight subjects. Mean values of the metabolic parameters measured during wakefulness were stable, with coefficients of variation below 7% (average values: CMRO2 = 118 ± 12 µmol O2/min/100 g, SvO2 = 67.0 ± 3.7% HbO2, CBF = 50.6 ±4.3 ml/min/100 g). During sleep, on average, CMRO2 decreased 21% (range: 14%-32%; average nadir = 98 ± 16 µmol O2/min/100 g), while EEG slow-wave activity, expressed in terms of δ -power, increased commensurately. Following sleep onset, CMRO2 was found to correlate negatively with relative δ -power (r = -0.6 to -0.8, P < 0.005), and positively with heart rate (r = 0.5 to 0.8, P < 0.0005). The data demonstrate that OxFlow MRI can noninvasively measure dynamic changes in cerebral metabolism associated with sleep, which should open new opportunities to study sleep physiology in health and disease.

A Noninvasive Method for Quantifying Cerebral Metabolic Rate of Oxygen by Hybrid PET/MRI: Validation in a Porcine Model

The gold standard for imaging the cerebral metabolic rate of oxygen (CMRO2) is positron emission tomography (PET); however, it is an invasive and complex procedure that also requires correction for recirculating 15O-H2O and the blood-borne activity. We propose a noninvasive reference-based hybrid PET/magnetic resonance imaging (MRI) method that uses functional MRI techniques to calibrate 15O-O2-PET data. Here, PET/MR imaging of oxidative metabolism (PMROx) was validated in an animal model by comparison to PET-alone measurements. Additionally, we investigated if the MRI-perfusion technique arterial spin labelling (ASL) could be used to further simplify PMROx by replacing 15O-H2O-PET, and if the PMROx was sensitive to anesthetics-induced changes in metabolism. Methods: 15O-H2O and 15O-O2 PET data were acquired in a hybrid PET/MR scanner (3 T Siemens Biograph mMR), together with simultaneous functional MRI (OxFlow and ASL), from juvenile pigs (n = 9). Animals were anesthetized with 3% isoflurane and 6 mL/kg/h propofol for the validation experiments and arterial sampling was performed for PET-alone measurements. PMROx estimates were obtained using whole-brain (WB) CMRO2 from OxFlow and local cerebral blood flow (CBF) from either noninvasive 15O-H2O-PET or ASL (PMROxASL). Changes in metabolism were investigated by increasing the propofol infusion to 20 mL/kg/h. Results: Good agreement and correlation were observed between regional CMRO2 measurements from PMROx and PET-alone. No significant differences were found between OxFlow and PET-only measurements of WB oxygen extraction fraction (0.30 ± 0.09 and 0.31 ± 0.09) and CBF (54.1 ± 16.7 and 56.6 ± 21.0 mL/100 g/min), or between PMROx and PET-only CMRO2 estimates (1.89 ± 0.16 and 1.81 ± 0.10 mLO2/100 g/min). Moreover, PMROx and PMROxASL were sensitive to propofol-induced reduction in CMRO2 Conclusion: This study provides initial validation of a noninvasive PET/MRI technique that circumvents many of the complexities of PET CMRO2 imaging. PMROx does not require arterial sampling and has the potential to reduce PET imaging to 15O-O2 only; however, future validation involving human participants are required.

A non-invasive reference-based method for imaging the cerebral metabolic rate of oxygen by PET/MR: theory and error analysis

Positron emission tomography (PET) remains the gold standard for quantitative imaging of the cerebral metabolic rate of oxygen (CMRO₂); however, it is an invasive and complex procedure that requires accounting for recirculating [¹⁵O]H₂O (RW) and the cerebral blood volume (CBV). This study presents a non-invasive reference-based technique for imaging CMRO₂ that was developed for PET/magnetic resonance imaging (MRI) with the goal of simplifying the PET procedure while maintaining its ability to quantify metabolism. The approach is to use whole-brain (WB) measurements of oxygen extraction fraction (OEF) and cerebral blood flow (CBF) to calibrate [¹⁵O]O₂-PET data, thereby avoiding the need for invasive arterial sampling. Here we present the theoretical framework, along with error analyses, sensitivity to PET noise and inaccuracies in input parameters, and initial assessment on PET data acquired from healthy participants. Simulations showed that neglecting RW and CBV corrections caused errors in CMRO₂ of less than ±10% for changes in regional OEF of ±25%. These predictions were supported by applying the reference-based approach to PET data, which resulted in remarkably similar CMRO₂ images to those generated by analyzing the same data using a modeling approach that incorporated the arterial input functions and corrected for CBV contributions. Significant correlations were observed between regional CMRO₂ values from the two techniques (slope = 1.00 ± 0.04, R ² > 0.98) with no significant differences found for integration times of 3 and 5 min. In summary, results demonstrate the feasibility of producing quantitative CMRO₂ images by PET/MRI without the need for invasive blood sampling.

Quantitative theory for the transverse relaxation time of blood water

An integrative model is proposed to describe the dependence of the transverse relaxation rate of blood water protons (R_2blood = 1/T_2blood ) on hematocrit fraction and oxygenation fraction (Y). This unified model takes into account (a) the diamagnetic effects of albumin, hemoglobin and the cell membrane; (b) the paramagnetic effect of hemoglobin; (c) the effect of compartmental exchange between plasma and erythrocytes under both fast and slow exchange conditions that vary depending on field strength and compartmental relaxation rates and (d) the effect of diffusion through field gradients near the erythrocyte membrane. To validate the model, whole-blood and lysed-blood R₂ data acquired previously using Carr-Purcell-Meiboom-Gill measurements as a function of inter-echo spacing τ_cp at magnetic fields of 3.0, 7.0, 9.4 and 11.7 T were fitted to determine the lifetimes (field-independent physiological constants) for water diffusion and exchange, as well as several physical constants, some of which are field-independent (magnetic susceptibilities) and some are field-dependent (relaxation rates for water protons in solutions of albumin and oxygenated and deoxygenated hemoglobin, ie, blood plasma and erythrocytes, respectively). This combined exchange-diffusion model allowed excellent fitting of the curve of the τ_cp -dependent relaxation rate dispersion at all four fields using a single average erythrocyte water lifetime, τ_ery = 9.1 ± 1.4 ms, and an averaged diffusional correlation time, τ_D = 3.15 ± 0.43 ms. Using this model and the determined physiological time constants and relaxation parameters, blood T₂ values published by multiple groups based on measurements at magnetic field strengths of 1.5 T and higher could be predicted correctly within error. Establishment of this theory is a fundamental step for quantitative modeling of the BOLD effect underlying functional MRI.

Cluster analysis of time evolution (CAT) for quantitative susceptibility mapping (QSM) and quantitative blood oxygen level-dependent magnitude (qBOLD)-based oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO2 ) mapping

PURPOSE

To improve the accuracy of QSM plus quantitative blood oxygen level-dependent magnitude (QSM + qBOLD or QQ)-based mapping of the oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO2 ) using cluster analysis of time evolution (CAT).

METHODS

3D multi-echo gradient echo and arterial spin labeling images were acquired in 11 healthy subjects and 5 ischemic stroke patients. DWI was also carried out on patients. CAT was developed for analyzing signal evolution over TE. QQ-based OEF and CMRO2 were reconstructed with and without CAT, and results were compared using region of interest analysis and a paired t-test.

RESULTS

Simulations demonstrated that CAT substantially reduced noise error in QQ-based OEF. In healthy subjects, QQ-based OEF appeared less noisy and more uniform with CAT than without CAT; average OEF with and without CAT in cortical gray matter was 32.7 ± 4.0% and 37.9 ± 4.5%, with corresponding CMRO2 of 148.4 ± 23.8 and 171.4 ± 22.4 μmol/100 g/min, respectively. In patients, regions of low OEF were confined within the ischemic lesions defined on DWI when using CAT, which was not observed without CAT.

CONCLUSION

The cluster analysis of time evolution (CAT) significantly improves the robustness of QQ-based OEF against noise.

Relationship between Abnormal Hyperintensity on T2-Weighted Images Around Developmental Venous Anomalies and Magnetic Susceptibility of Their Collecting Veins: In-Vivo Quantitative Susceptibility Mapping Study

Objective

A developmental venous anomaly (DVA) is a vascular malformation of ambiguous clinical significance. We aimed to quantify the susceptibility of draining veins (χ_vein) in DVA and determine its significance with respect to oxygen metabolism using quantitative susceptibility mapping (QSM).

Materials and Methods

Brain magnetic resonance imaging of 27 consecutive patients with incidentally detected DVAs were retrospectively reviewed. Based on the presence of abnormal hyperintensity on T2-weighted images (T2WI) in the brain parenchyma adjacent to DVA, the patients were grouped into edema (E+, n = 9) and non-edema (E−, n = 18) groups. A 3T MR scanner was used to obtain fully flow-compensated gradient echo images for susceptibility-weighted imaging with source images used for QSM processing. The χ_vein was measured semi-automatically using QSM. The normalized χ_vein was also estimated. Clinical and MR measurements were compared between the E+ and E− groups using Student's t-test or Mann-Whitney U test. Correlations between the χ_vein and area of hyperintensity on T2WI and between χ_vein and diameter of the collecting veins were assessed. The correlation coefficient was also calculated using normalized veins.

Results

The DVAs of the E+ group had significantly higher χ_vein (196.5 ± 27.9 vs. 167.7 ± 33.6, p = 0.036) and larger diameter of the draining veins (p = 0.006), and patients were older (p = 0.006) than those in the E− group. The χ_vein was also linearly correlated with the hyperintense area on T2WI (r = 0.633, 95% confidence interval 0.333–0.817, p < 0.001).

Conclusion

DVAs with abnormal hyperintensity on T2WI have higher susceptibility values for draining veins, indicating an increased oxygen extraction fraction that might be associated with venous congestion.

A computational model integrating brain electrophysiology and metabolism highlights the key role of extracellular potassium and oxygen

The human brain is a small organ which uses a disproportionate amount of the total metabolic energy production in the body. While it is well understood that the most significant energy sink is the maintenance of the neuronal membrane potential during the brain signaling activity, the role of astrocytes in the energy balance continues to be the topic of a lot of research. A key function of astrocytes, besides clearing glutamate from the synaptic clefts, is the potassium clearing after neuronal activation. Extracellular potassium plays a significant role in triggering neuronal firing, and elevated concentration of potassium may lead to abnormal firing patterns, e.g., seizures, thus emphasizing the importance of the glial K⁺ buffering role. The predictive mathematical model proposed in this paper elucidates the role of glial potassium clearing in brain energy metabolism, integrating a detailed model of the ion dynamics which regulates neuronal firing with a four compartment metabolic model. Because of the very different characteristic time scales of electrophysiology and metabolism, care must be taken when coupling the two models to ensure that the predictions, e.g., neuronal firing frequencies and the oxygen-glucose index (OGI) of the brain during activation and rest, are in agreement with empirical observations. The temporal multi-scale nature of the problem requires the design of new computational tools to ensure a stable and accurate numerical treatment. The model predictions for different protocols, including combinations of elevated activation and ischemic episodes, are in good agreement with experimental observations reported in the literature.

Quantitative susceptibility mapping-based cerebral metabolic rate of oxygen mapping with minimum local variance

PURPOSE

The objective of this study was to demonstrate the feasibility of a cerebral metabolic rate of oxygen (CMRO₂ ) mapping method based on its minimum local variance (MLV) without vascular challenge using quantitative susceptibility mapping (QSM) and cerebral blood flow (CBF).

METHODS

Three-dimensional multi-echo gradient echo imaging and arterial spin labeling were performed in 11 healthy subjects to calculate QSM and CBF. Minimum local variance was used to compute whole-brain CMRO₂ map from QSM and CBF. The MLV method was compared with a reference method using the caffeine challenge. Their agreement within the cortical gray matter (CGM) was assessed on CMRO₂ and oxygen extraction fraction (OEF) maps at both baseline and challenge states.

RESULTS

Mean CMRO₂ (in µmol/100 g/min) obtained in CGM using the caffeine challenge and MLV were 142 ± 16.5 and 139 ± 14.8 µmol/100 g/min, respectively; the corresponding baseline OEF were 33.0 ± 4.0% and 31.8 ± 3.2%, respectively. The MLV and caffeine challenge methods showed no statistically significant differences across subjects with small ( < 4%) biases in CMRO₂ and OEF values.

CONCLUSIONS

Minimum local variance-based CMRO₂ mapping without vascular challenge using QSM and arterial spin labeling is feasible in healthy subjects. Magn Reson Med 79:172-179, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Analytical Expression for the NO Concentration Profile Following NONOate Decomposition in the Presence of Oxygen

We have derived an analytical expression for the time dependent NO concentration from NONOate donors in the presence of oxygen for the process of NO release from NO donors following autoxidation. This analytical solution incorporates the kinetics of the releases with the autoxidation and is used to fit the simulated NO concentration profile to the experimental data. This allows one to determine the NO release rate constant, k 1, the NO release stoichiometric coefficient, v NO, and the NO autoxidation reaction rate constant, k 2. This analytical solution also allows us to predict the real NO concentration released from NO donors under aerobic conditions, while v NO is reportedly two under aerobic conditions, it falls to lower values in the presence of oxygen.

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

PURPOSE

To quantitatively map cerebral metabolic rate of oxygen ( CMRO2) and oxygen extraction fraction ( OEF) in human brains using quantitative susceptibility mapping (QSM) and arterial spin labeling-measured cerebral blood flow (CBF) before and after caffeine vasoconstriction.

METHODS

Using the multiecho, three-dimensional gradient echo sequence and an oral bolus of 200 mg caffeine, whole brain CMRO2 and OEF were mapped at 3-mm isotropic resolution on 13 healthy subjects. The QSM-based CMRO2 was compared with an R2*-based CMRO2 to analyze the regional consistency within cortical gray matter (CGM) with the scaling in the R2* method set to provide same total CMRO2 as the QSM method for each subject.

RESULTS

Compared to precaffeine, susceptibility increased (5.1 ± 1.1 ppb; P < 0.01) and CBF decreased (-23.6 ± 6.7 ml/100 g/min; P < 0.01) at 25-min postcaffeine in CGM. This corresponded to a CMRO2 of 153.0 ± 26.4 μmol/100 g/min with an OEF of 33.9 ± 9.6% and 54.5 ± 13.2% (P < 0.01) pre- and postcaffeine, respectively, at CGM, and a CMRO2 of 58.0 ± 26.6 μmol/100 g/min at white matter. CMRO2 from both QSM- and R2*-based methods showed good regional consistency (P > 0.05), but quantitation of R2*-based CMRO2 required an additional scaling factor.

CONCLUSION

QSM can be used with perfusion measurements pre- and postcaffeine vascoconstriction to map CMRO2 and OEF.

Pulse sequence programming in a dynamic visual environment: SequenceTree

PURPOSE

To describe SequenceTree, an open source, integrated software environment for implementing MRI pulse sequences and, ideally, exporting them to actual MRI scanners. The software is a user-friendly alternative to vendor-supplied pulse sequence design and editing tools and is suited for programmers and nonprogrammers alike.

METHODS

The integrated user interface was programmed using the Qt4/C++ toolkit. As parameters and code are modified, the pulse sequence diagram is automatically updated within the user interface. Several aspects of pulse programming are handled automatically, allowing users to focus on higher-level aspects of sequence design. Sequences can be simulated using a built-in Bloch equation solver and then exported for use on a Siemens MRI scanner. Ideally, other types of scanners will be supported in the future.

RESULTS

SequenceTree has been used for 8 years in our laboratory and elsewhere and has contributed to more than 50 peer-reviewed publications in areas such as cardiovascular imaging, solid state and nonproton NMR, MR elastography, and high-resolution structural imaging.

CONCLUSION

SequenceTree is an innovative, open source, visual pulse sequence environment for MRI combining simplicity with flexibility and is ideal both for advanced users and users with limited programming experience.

Feasibility and reproducibility of measurement of whole muscle blood flow, oxygen extraction, and VO2 with dynamic exercise using MRI

PURPOSE

Develop an MRI method to estimate skeletal muscle oxygen consumption (VO2 ) with dynamic exercise using simultaneous measurement of venous blood flow (VBF) and venous oxygen saturation (SvO2 ).

METHODS

Real-time imaging of femoral VBF using a complex-difference method was interleaved with imaging of venous hemoglobin oxygen saturation (SvO2 ) using magnetic susceptometry to estimate muscle VO2 (Fick principle). Nine healthy subjects performed repeated 5-watt knee-extension (quadriceps) exercise within the bore of a 1.5 Tesla MRI scanner, for test/re-test comparison. VBF, SvO2 , and derived VO2 were estimated at baseline and immediately (<1 s) postexercise and every 2.4 s for 4 min.

RESULTS

Quadriceps muscle mass was 2.43 ± 0.31 kg. Mean baseline values were VBF = 0.13 ± 0.06 L/min/kg, SvO2 = 69.4 ± 10.1%, and VO2 = 6.8 ± 4.1 mL/min/kg. VBF, SvO2 , and VO2 values from peak exercise had good agreement between trials (VBF = 0.9 ± 0.1 versus 1.0 ± 0.1 L/min/kg, R(2) = 0.83, CV = 7.6%; SvO2 = 43.2 ± 13.5 versus 40.9 ± 13.1%, R(2) = 0.88, CV = 15.6%; VO2 = 95.7 ± 18.0 versus 108.9 ± 17.3 mL/min/kg, R(2) = 0.88, CV = 12.3%), as did the VO2 recovery time constant (26.1 ± 3.5 versus 26.0 ± 4.0 s, R(2) = 0.85, CV = 6.0%). CV = coefficient of variation.

CONCLUSION

Rapid imaging of VBF and SvO2 for the estimation of whole muscle VO2 is compatible with dynamic exercise for the estimation of peak values and recovery dynamics following exercise with good reproducibility.

Investigation of the neurovascular coupling in positive and negative BOLD responses in human brain at 7 T

Decreases in stimulus-dependent blood oxygenation level dependent (BOLD) signal and their underlying neurovascular origins have recently gained considerable interest. In this study a multi-echo, BOLD-corrected vascular space occupancy (VASO) functional magnetic resonance imaging (fMRI) technique was used to investigate neurovascular responses during stimuli that elicit positive and negative BOLD responses in human brain at 7 T. Stimulus-induced BOLD, cerebral blood volume (CBV), and cerebral blood flow (CBF) changes were measured and analyzed in 'arterial' and 'venous' blood compartments in macro- and microvasculature. We found that the overall interplay of mean CBV, CBF and BOLD responses is similar for tasks inducing positive and negative BOLD responses. Some aspects of the neurovascular coupling however, such as the temporal response, cortical depth dependence, and the weighting between 'arterial' and 'venous' contributions, are significantly different for the different task conditions. Namely, while for excitatory tasks the BOLD response peaks at the cortical surface, and the CBV change is similar in cortex and pial vasculature, inhibitory tasks are associated with a maximum negative BOLD response in deeper layers, with CBV showing strong constriction of surface arteries and a faster return to baseline. The different interplays of CBV, CBF and BOLD during excitatory and inhibitory responses suggests different underlying hemodynamic mechanisms.