Quantification of whole-organ individual and bilateral renal metabolic rate of oxygen.

PURPOSE: Renal metabolic rate of oxygen (rMRO2 ) is a potentially important biomarker of kidney function. The key parameters for rMRO2 quantification include blood flow rate (BFR) and venous oxygen saturation (SvO2 ) in a draining vessel. Previous approaches to quantify renal metabolism have focused on the single organ. Here, both kidneys are considered as one unit to quantify bilateral rMRO2 . A pulse sequence to facilitate bilateral rMRO2 quantification is introduced. METHODS: To quantify bilateral rMRO2 , measurements of BFR and SvO2 are made along the inferior vena cava (IVC) at suprarenal and infrarenal locations. From the continuity equation, these four parameters can be related to derive an expression for bilateral rMRO2 . The recently reported K-MOTIVE pulse sequence was implemented at four locations: left kidney, right kidney, suprarenal IVC, and infrarenal IVC. A dual-band variant of K-MOTIVE (db-K-MOTIVE) was developed by incorporating simultaneous-multi-slice imaging principles. The sequence simultaneously measures BFR and SvO2 at suprarenal and infrarenal locations in a single pass of 21 s, yielding bilateral rMRO2 . RESULTS: SvO2 and BFR are higher in suprarenal versus infrarenal IVC, and the renal veins are highly oxygenated (SvO2 >90%). Bilateral rMRO2 quantified in 10 healthy subjects (8 M, 30 ± 8 y) was found to be 291 ± 247 and 349 ± 300 (µmolO2 /min)/100 g, derived from K-MOTIVE and db-K-MOTIVE, respectively. In comparison, total rMRO2 from combining left and right was 329 ± 273 (µmolO2 /min)/100 g. CONCLUSION: The present work demonstrates that bilateral rMRO2 quantification is feasible with fair reproducibility and physiological plausibility. The indirect method is a promising approach to compute bilateral rMRO2 when individual rMRO2 quantification is difficult.

Cranial bone imaging using ultrashort echo-time bone-selective MRI as an alternative to gradient-echo based "black-bone" techniques.

OBJECTIVES: CT is the clinical standard for surgical planning of craniofacial abnormalities in pediatric patients. This study evaluated three MRI cranial bone imaging techniques for their strengths and limitations as a radiation-free alternative to CT. METHODS: Ten healthy adults were scanned at 3 T with three MRI sequences: dual-radiofrequency and dual-echo ultrashort echo time sequence (DURANDE), zero echo time (ZTE), and gradient-echo (GRE). DURANDE bright-bone images were generated by exploiting bone signal intensity dependence on RF pulse duration and echo time, while ZTE bright-bone images were obtained via logarithmic inversion. Three skull segmentations were derived, and the overlap of the binary masks was quantified using dice similarity coefficient. Craniometric distances were measured, and their agreement was quantified. RESULTS: There was good overlap of the three masks and excellent agreement among craniometric distances. DURANDE and ZTE showed superior air-bone contrast (i.e., sinuses) and soft-tissue suppression compared to GRE. DISCUSSIONS: ZTE has low levels of acoustic noise, however, ZTE images had lower contrast near facial bones (e.g., zygomatic) and require effective bias-field correction to separate bone from air and soft-tissue. DURANDE utilizes a dual-echo subtraction post-processing approach to yield bone-specific images, but the sequence is not currently manufacturer-supported and requires scanner-specific gradient-delay corrections.

Cerebral oxygen metabolism from MRI susceptibility.

This article provides an overview of MRI methods exploiting magnetic susceptibility properties of blood to assess cerebral oxygen metabolism, including the tissue oxygen extraction fraction (OEF) and the cerebral metabolic rate of oxygen (CMRO2). The first section is devoted to describing blood magnetic susceptibility and its effect on the MRI signal. Blood circulating in the vasculature can have diamagnetic (oxyhemoglobin) or paramagnetic properties (deoxyhemoglobin). The overall balance between oxygenated and deoxygenated hemoglobin determines the induced magnetic field which, in turn, modulates the transverse relaxation decay of the MRI signal via additional phase accumulation. The following sections of this review then illustrate the principles underpinning susceptibility-based techniques for quantifying OEF and CMRO2. Here, it is detailed whether these techniques provide global (OxFlow) or local (Quantitative Susceptibility Mapping - QSM, calibrated BOLD - cBOLD, quantitative BOLD - qBOLD, QSM+qBOLD) measurements of OEF or CMRO2, and what signal components (magnitude or phase) and tissue pools they consider (intravascular or extravascular). Validations studies and potential limitations of each method are also described. The latter include (but are not limited to) challenges in the experimental setup, the accuracy of signal modeling, and assumptions on the measured signal. The last section outlines the clinical uses of these techniques in healthy aging and neurodegenerative diseases and contextualizes these reports relative to results from gold-standard PET.

MRI Quantification of Cortical Bone Porosity, Mineralization, and Morphologic Structure in Postmenopausal Osteoporosis.

Background Preclinical studies have suggested that solid-state MRI markers of cortical bone porosity, morphologic structure, mineralization, and osteoid density are useful measures of bone health. Purpose To explore whether MRI markers of cortical bone porosity, morphologic structure, mineralization, and osteoid density are affected in postmenopausal osteoporosis (OP) and to examine associations between MRI markers and bone mineral density (BMD) in postmenopausal women. Materials and Methods In this single-center study, postmenopausal women were prospectively recruited from January 2019 to October 2020 into two groups: participants with OP who had not undergone treatment, defined as having any dual-energy x-ray absorptiometry (DXA) T-score of -2.5 or less, and age-matched control participants without OP (hereafter, non-OP). Participants underwent MRI in the midtibia, along with DXA in the hip and spine, and peripheral quantitative CT in the midtibia. Specifically, MRI measures of cortical bone porosity (pore water and total water), osteoid density (bound water [BW]), morphologic structure (cortical bone thickness), and mineralization (phosphorous [P] density [31P] and 31P-to-BW concentration ratio) were quantified at 3.0 T. MRI measures were compared between OP and non-OP groups and correlations with BMD were assessed. Results Fifteen participants with OP (mean age, 63 years ± 5 [SD]) and 19 participants without OP (mean age, 65 years ± 6) were evaluated. The OP group had elevated pore water (11.6 mol/L vs 9.5 mol/L; P = .007) and total water densities (21.2 mol/L vs 19.7 mol/L; P = .03), and had lower cortical bone thickness (4.8 mm vs 5.6 mm; P < .001) and 31P density (6.4 mol/L vs 7.5 mol/L; P = .01) than the non-OP group, respectively, although there was no evidence of a difference in BW or 31P-to-BW concentration ratio. Pore and total water densities were inversely associated with DXA and peripheral quantitative CT BMD (P < .001), whereas cortical bone thickness and 31P density were positively associated with DXA and peripheral quantitative CT BMD (P = .01). BW, 31P density, and 31P-to-BW concentration ratio were positively associated with DXA (P < .05), but not with peripheral quantitative CT. Conclusion Solid-state MRI of cortical bone was able to help detect potential impairments in parameters reflecting porosity, morphologic structure, and mineralization in postmenopausal osteoporosis. © RSNA, 2023 Supplemental material is available for this article. See also the editorial by Bae in this issue.

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

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

Alternating unbalanced SSFP for 3D R2' mapping of the human brain.

PURPOSE: Measuring the transverse-relaxation rate R2' provides valuable information in quantitative evaluation of tissue microstructure, for example, in terms of oxygenation levels. Here, we propose an alternating unbalanced SSFP pulse sequence for rapid whole-brain 3D R2' mapping. METHODS: Unlike currently practiced, spin echo-based R2' measurement techniques, the proposed method alternates between SSFP-FID and SSFP-ECHO modes for rapid 3D encoding of transverse relaxation rates expressed as R2 + R2' and R2-R2' . Z-shimming gradients embedded into multi-echo trains of each SSFP module are designed to achieve relative immunity to large-scale magnetic-field variations (ΔB0 ). Appropriate models for the temporal evolution of the two groups of SSFP signals were constructed with ΔB0 -induced modulations accounted for, leading to ΔB0 -corrected estimation of R2 , R2' , and R2∗ (= R2 + R2' ). Additionally, relative magnetic susceptibility (Δχ) maps were obtained by quantitative susceptibility mapping of the phase data. Numerical simulations were performed to optimize scan parameters, followed by in vivo studies at 3 T in 7 healthy subjects. Measured parameters were evaluated in six brain regions, and subjected to interparameter correlation analysis. RESULTS: The resultant maps of R2' and additionally derived R2 , R2∗ , and Δχ all demonstrated the expected contrast across brain territories (eg, deep brain structures versus cortex), with the measured values in good agreement with previous reports. Furthermore, regression analyses yielded strong linear relationships for the transverse relaxation parameters ( R2' , R2 , and R2∗ ) against Δχ. CONCLUSION: Results suggest feasibility of the proposed method as a practical and reliable means for measuring R2' , R2 , R2∗ , and Δχ across the entire brain.

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.

Impact of gradient imperfections on bone water quantification with UTE MRI.

PURPOSE: The impact of gradient imperfections on UTE images and UTE image-derived bone water quantification was investigated at 3 T field strength. METHODS: The effects of simple gradient time delays and eddy currents on UTE images, as well as the effects of gradient error corrections, were studied with simulation and phantom experiments. The k-space trajectory was mapped with a 2D sequence with phase encoding on both spatial axes by measuring the phase of the signal in small time increments during ramp-up of the read gradient. In vivo 3D UTE images were reconstructed with and without gradient error compensation to determine the bias in bone water quantification. Finally, imaging was performed on 2 equally configured Siemens TIM Trio systems (Siemens Medical Solutions, Erlangen, Germany) to investigate the impact of such gradient imperfections on inter-scanner measurement bias. RESULTS: Compared to values derived from UTE images with full gradient error compensation, total bone water was found to deviate substantially with no (up to 17%) or partial (delay-only) compensation (up to 10.8%). Bound water, obtained with inversion recovery-prepared UTE, was somewhat less susceptible to gradient errors (up to 2.2% for both correction strategies). Inter-scanner comparison indicated a statistically significant bias between measurements from the 2 MR systems for both total and bound water, which either vanished or was substantially reduced following gradient error correction. CONCLUSION: Gradient imperfections impose spatially dependent artifacts on UTE images, which compromise not only bone water quantification accuracy but also inter-scanner measurement agreement if left uncompensated.

Rapid dual-RF, dual-echo, 3D ultrashort echo time craniofacial imaging: A feasibility study.

PURPOSE: To develop a dual-radiofrequency (RF), dual-echo, 3D ultrashort echo-time (UTE) pulse sequence and bone-selective image reconstruction for rapid high-resolution craniofacial MRI. METHODS: The proposed pulse sequence builds on recently introduced dual-RF UTE imaging. While yielding enhanced bone specificity by exploiting high sensitivity of short T2 signals to variable RF pulse widths, the parent technique exacts a 2-fold scan time penalty relative to standard dual-echo UTE. In the proposed method, the parent sequence's dual-RF scheme was incorporated into dual-echo acquisitions while radial view angles are varied every pulse-to-pulse repetition period. The resulting 4 echoes (2 for each RF) were combined by view-sharing to construct 2 sets of k-space data sets, corresponding to short and long TEs, respectively, leading to a 2-fold increase in imaging efficiency. Furthermore, by exploiting the sparsity of bone signals in echo-difference images, acceleration was achieved by solving a bone-sparsity constrained image reconstruction problem. In vivo studies were performed to evaluate the effectiveness of the proposed acceleration approaches in comparison to the parent method. RESULTS: The proposed technique achieves 1.1-mm isotropic skull imaging in 3 minutes without visual loss of image quality, compared to the parent technique (scan time = 12 minutes). Bone-specific images and corresponding 3D renderings of the skull were found to depict the expected craniofacial anatomy over the entire head. CONCLUSION: The proposed method is able to achieve high-resolution volumetric craniofacial images in a clinically practical imaging time, and thus may prove useful as a potential alternative to computed tomography.

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

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

T2 -prepared balanced steady-state free precession (bSSFP) for quantifying whole-blood oxygen saturation at 1.5T.

PURPOSE: To establish a calibration equation to convert human blood T2 to the full range of oxygen saturation levels (HbO2 ) and physiologic hematocrit (Hct) values using a T2 -prepared balanced steady-state free precession sequence (T2 -SSFP) at 1.5T. METHODS: Blood drawn from 10 healthy donors (29.1 ± 3.9 years old) was prepared into samples of varying HbO2 and Hct (n = 79), and imaged using T2 -SSFP sequence at 37°C and interrefocusing interval τ180 = 12 ms. The relationship between blood T2 , HbO2 , and Hct was established based on the model R2=R2,plasma+Hct (R2,RBC-R2,plasma)+k·Hct·(1-Hct)·(1-HbO2)2. Measured R2 and HbO2 levels were fit by the model yielding values of R2,plasma, R2,RBC, and k. T2 -SSFP and the established calibration equation were applied to extract HbO2 at the superior sagittal sinus (SSS) in vivo and were compared with susceptometry-based oximetry. RESULTS: Constants derived from the fit were: k = 74.2 [s-1 ], R2,plasma = 1.5 [s-1 ], R2,RBC = 11.6 [s-1 ], the R2 of the fit was 0.95. Average HbO2 at the SSS in seven healthy volunteers was 65% ± 7% and 66% ± 7% via T2 - and susceptometry-based oximetry, respectively. Bland-Altman analysis indicated agreement between the two oximetric methods with no significant bias. CONCLUSION: The calibration constants presented here should ensure improved accuracy for whole-blood oximetry based on T2 -SSFP at 1.5T. Magn Reson Med 79:1893-1900, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Multiplexed MRI methods for rapid estimation of global cerebral metabolic rate of oxygen consumption.

The global cerebral metabolic rate of oxygen (CMRO2), which reflects metabolic activity of the brain under various physiologic conditions, can be quantified using a method, referred to as 'OxFlow', which simultaneously measures hemoglobin oxygen saturation in a draining vein (Yv) and total cerebral blood flow (tCBF). Conventional OxFlow (Conv-OxFlow) entails four interleaves incorporated in a single pulse sequence - two for phase-contrast based measurement of tCBF in the supplying arteries of the neck, and two to measure the intra- to extravascular phase difference in the superior sagittal sinus to derive Yv [Jain et al., JCBFM 2010]. However, this approach limits achievable temporal resolution thus precluding capture of rapid changes of brain metabolic states such as the response to apneic stimuli. Here, we developed a time-efficient, multiplexed OxFlow method and evaluated its potential for measuring dynamic alterations in global CMRO2 during a breath-hold challenge. Two different implementations of multiplexed OxFlow were investigated: 1) simultaneous-echo-refocusing based OxFlow (SER-OxFlow) and 2) simultaneous-multi-slice imaging-based dual-band OxFlow (DB-OxFlow). The two sequences were implemented on 3T scanners (Siemens TIM Trio and Prisma) and their performance was evaluated in comparison to Conv-OxFlow in ten healthy subjects for baseline CMRO2 quantification. Comparison of measured parameters (Yv, tCBF, CMRO2) revealed no significant bias of SER-OxFlow and DB-OxFlow, with respect to the reference Conv-OxFlow while improving temporal resolution two-fold (12.5 versus 25s). Further acceleration shortened scan time to 8 and 6s for SER and DB-OxFlow, respectively, for time-resolved CMRO2 measurement. The two sequences were able of capturing smooth transitions of Yv, tCBF, and CMRO2 over the time course consisting of 30s of normal breathing, 30s of volitional apnea, and 90s of recovery. While both SER- and DB-OxFlow techniques provide significantly improved temporal resolution (by a factor of 3 - 4 relative to Conv-OxFlow), DB-OxFlow was found to be superior for the study of short physiologic stimuli.

Current-induced alternating reversed dual-echo-steady-state for joint estimation of tissue relaxation and electrical properties.

PURPOSE: To develop a current-induced, alternating reversed dual-echo-steady-state-based magnetic resonance electrical impedance tomography for joint estimation of tissue relaxation and electrical properties. METHODS: The proposed method reverses the readout gradient configuration of conventional, in which steady-state-free-precession (SSFP)-ECHO is produced earlier than SSFP-free-induction-decay (FID) while alternating current pulses are applied in between the two SSFPs to secure high sensitivity of SSFP-FID to injection current. Additionally, alternating reversed dual-echo-steady-state signals are modulated by employing variable flip angles over two orthogonal injections of current pulses. Ratiometric signal models are analytically constructed, from which T1 , T2 , and current-induced Bz are jointly estimated by solving a nonlinear inverse problem for conductivity reconstruction. Numerical simulations and experimental studies are performed to investigate the feasibility of the proposed method in estimating relaxation parameters and conductivity. RESULTS: The proposed method, if compared with conventional magnetic resonance electrical impedance tomography, enables rapid data acquisition and simultaneous estimation of T1 , T2 , and current-induced Bz , yielding a comparable level of signal-to-noise ratio in the parameter estimates while retaining a relative conductivity contrast. CONCLUSION: We successfully demonstrated the feasibility of the proposed method in jointly estimating tissue relaxation parameters as well as conductivity distributions. It can be a promising, rapid imaging strategy for quantitative conductivity estimation. Magn Reson Med 78:107-120, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Rapid whole-brain gray matter imaging using single-slab three-dimensional dual-echo fast spin echo: A feasibility study.

PURPOSE: To achieve rapid, high resolution whole-brain gray matter (GM) imaging by developing a novel, single-slab three-dimensional dual-echo fast-spin-echo pulse sequence and GM-selective reconstruction. METHODS: Unlike conventional GM imaging that uses time-consuming double-inversion-recovery preparation, the proposed pulse sequence was designed to have two split portions along the echo train, in which the first half was dedicated to yield short inversion recovery (IR)-induced white matter suppression and variable-flip-angle-induced two-step GM signal evolution while the second half cerebrospinal fluid-only signals. Multi-step variable-flip-angle schedules and sampling reordering were optimized to yield high GM signals while balancing cerebrospinal fluid signals between ECHOes. GM-selective images were then reconstructed directly from the weighted subtraction between ECHOes by solving a sparse signal recovery problem. In vivo studies were performed to validate the effectiveness of the proposed method over conventional double-inversion-recovery. RESULTS: The proposed method, while achieving one millimeter isotropic, whole-brain GM imaging within 5.5 min, showed superior performance than conventional double-inversion-recovery in producing GM-only images without apparent artifacts and noise. CONCLUSION: We successfully demonstrated the feasibility of the proposed method in achieving whole-brain GM imaging in a clinically acceptable imaging time. The proposed method is expected to be a promising alternative to conventional double-inversion-recovery in clinical applications. Magn Reson Med 78:1691-1699, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Alternating steady state free precession for estimation of current-induced magnetic flux density: A feasibility study.

PURPOSE: To develop a novel, current-controlled alternating steady-state free precession (SSFP)-based conductivity imaging method and corresponding MR signal models to estimate current-induced magnetic flux density (Bz ) and conductivity distribution. METHODS: In the proposed method, an SSFP pulse sequence, which is in sync with alternating current pulses, produces dual oscillating steady states while yielding nonlinear relation between signal phase and Bz . A ratiometric signal model between the states was analytically derived using the Bloch equation, wherein Bz was estimated by solving a nonlinear inverse problem for conductivity estimation. A theoretical analysis on the signal-to-noise ratio of Bz was given. Numerical and experimental studies were performed using SSFP-FID and SSFP-ECHO with current pulses positioned either before or after signal encoding to investigate the feasibility of the proposed method in conductivity estimation. RESULTS: Given all SSFP variants herein, SSFP-FID with alternating current pulses applied before signal encoding exhibits the highest Bz signal-to-noise ratio and conductivity contrast. Additionally, compared with conventional conductivity imaging, the proposed method benefits from rapid SSFP acquisition without apparent loss of conductivity contrast. CONCLUSION: We successfully demonstrated the feasibility of the proposed method in estimating current-induced Bz and conductivity distribution. It can be a promising, rapid imaging strategy for quantitative conductivity imaging.

Validation of a new 3D quantitative BOLD based cerebral oxygen extraction mapping.

Quantitative BOLD (qBOLD) MRI allows evaluation of oxidative metabolism of the brain based purely on an endogenous contrast mechanism. The method quantifies deoxygenated blood volume (DBV) and hemoglobin oxygen saturation level of venous blood (Yv), yielding oxygen extraction fraction (OEF), and along with a separate measurement of cerebral blood flow, cerebral metabolic rate of oxygen (CMRO2) maps. Here, we evaluated our recently reported 3D qBOLD method that rectifies a number of deficiencies in prior qBOLD approaches in terms of repeat reproducibility and sensitivity to hypercapnia on the metabolic parameters, and in comparison to dual-gas calibrated BOLD (cBOLD) MRI for determining resting-state oxygen metabolism. Results suggested no significant difference between test-retest qBOLD scans in either DBV and OEF. Exposure to hypercapnia yielded group averages of 38 and 28% for OEF and 151 and 146 µmol/min/100 g for CMRO2 in gray matter at baseline and hypercapnia, respectively. The decrease of OEF during hypercapnia was significant (p ≪ 0.01), whereas CMRO2 did not change significantly (p = 0.25). Finally, baseline OEF (37 vs. 39%) and CMRO2 (153 vs. 145 µmol/min/100 g) in gray matter using qBOLD and dual-gas cBOLD were found to be in good agreement with literature values, and were not significantly different from each other (p > 0.1).

Automated Quantification of Total Cerebral Blood Flow from Phase-Contrast MRI and Deep Learning.

Knowledge of input blood to the brain, which is represented as total cerebral blood flow (tCBF), is important in evaluating brain health. Phase-contrast (PC) magnetic resonance imaging (MRI) enables blood velocity mapping, allowing for noninvasive measurements of tCBF. In the procedure, manual selection of brain-feeding arteries is an essential step, but is time-consuming and often subjective. Thus, the purpose of this work was to develop and validate a deep learning (DL)-based technique for automated tCBF quantifications. To enhance the DL segmentation performance on arterial blood vessels, in the preprocessing step magnitude and phase images of PC MRI were multiplied several times. Thereafter, a U-Net was trained on 218 images for three-class segmentation. Network performance was evaluated in terms of the Dice coefficient and the intersection-over-union (IoU) on 40 test images, and additionally, on externally acquired 20 datasets. Finally, tCBF was calculated from the DL-predicted vessel segmentation maps, and its accuracy was statistically assessed with the correlation of determination (R2), the intraclass correlation coefficient (ICC), paired t-tests, and Bland-Altman analysis, in comparison to manually derived values. Overall, the DL segmentation network provided accurate labeling of arterial blood vessels for both internal (Dice=0.92, IoU=0.86) and external (Dice=0.90, IoU=0.82) tests. Furthermore, statistical analyses for tCBF estimates revealed good agreement between automated versus manual quantifications in both internal (R2=0.85, ICC=0.91, p=0.52) and external (R2=0.88, ICC=0.93, p=0.88) test groups. The results suggest feasibility of a simple and automated protocol for quantifying tCBF from neck PC MRI and deep learning.

SNR-optimized phase-sensitive dual-acquisition turbo spin echo imaging: a fast alternative to FLAIR.

Phase-sensitive dual-acquisition single-slab three-dimensional turbo spin echo imaging was recently introduced, producing high-resolution isotropic cerebrospinal fluid attenuated brain images without long inversion recovery preparation. Despite the advantages, the weighted-averaging-based technique suffers from noise amplification resulting from different levels of cerebrospinal fluid signal modulations over the two acquisitions. The purpose of this work is to develop a signal-to-noise ratio-optimized version of the phase-sensitive dual-acquisition single-slab three-dimensional turbo spin echo. Variable refocusing flip angles in the first acquisition are calculated using a three-step prescribed signal evolution while those in the second acquisition are calculated using a two-step pseudo-steady state signal transition with a high flip-angle pseudo-steady state at a later portion of the echo train, balancing the levels of cerebrospinal fluid signals in both the acquisitions. Low spatial frequency signals are sampled during the high flip-angle pseudo-steady state to further suppress noise. Numerical simulations of the Bloch equations were performed to evaluate signal evolutions of brain tissues along the echo train and optimize imaging parameters. In vivo studies demonstrate that compared with conventional phase-sensitive dual-acquisition single-slab three-dimensional turbo spin echo, the proposed optimization yields 74% increase in apparent signal-to-noise ratio for gray matter and 32% decrease in imaging time. The proposed method can be a potential alternative to conventional fluid-attenuated imaging.

Six-minute, in vivo MRI quantification of proximal femur trabecular bone 3D microstructure.

BACKGROUND: Assessment of proximal femur trabecular bone microstructure in vivo by magnetic resonance imaging has recently been validated for acquiring information independent of bone mineral density in osteoporotic patients. However, the requisite signal-to-noise ratio (SNR) and resolution for interrogation of the trabecular microstructure at this anatomical location prolongs the scan duration and renders the imaging protocol clinically infeasible. Parallel imaging and compressed sensing (PICS) techniques can reduce the scan duration of the imaging protocol without substantially compromising image quality. The present work investigates the limits of acceleration for a commonly used PICS technique, ℓ1-ESPIRiT, for the purpose of quantifying measures of trabecular bone microarchitecture. Based on a desired error tolerance, a six-minute, prospectively accelerated variant of the imaging protocol was developed and assessed for intersession reproducibility and agreement with the longer reference scan. PURPOSE: To investigate the limits of acceleration for MRI-based trabecular bone quantification by parallel imaging and compressed sensing reconstruction, and to develop a prototypical imaging protocol for assessing the proximal femur microstructure in a clinically practical scan time. METHODS: Healthy participants (n = 11) were scanned by a 3D balanced steady-state free precession (bSSFP) sequence satisfying the Nyquist criterion with a scan duration of about 18 min. The raw data were retrospectively undersampled and reconstructed to mimic various acceleration factors ranging from 2 to 6. Trabecular volumes-of-interest in four major femoral regions (greater trochanter, intertrochanteric region, femoral neck, and femoral head) were analyzed and six relevant measures of trabecular bone microarchitecture (bone volume fraction, surface-to-curve ratio, erosion index, elastic modulus, trabecular thickness, plates-to-rods ratio) were obtained for images of all accelerations. To assess agreement, median percent error and intraclass correlation coefficients (ICCs) were computed using the fully-sampled data as reference. Based on this analysis, a prospectively 3-fold accelerated sequence with a duration of about 6 min was developed and the analysis was repeated. RESULTS: A prospective acceleration factor of 3 demonstrated comparable performance in reproducibility and absolute agreement to the fully-sampled scan. The median CoV over all image-derived metrics was generally <6 % and ICCs >0.70. Also, measurements from prospectively 3-fold accelerated scans demonstrated in general median percent errors of <7 % and ICCs >0.70. CONCLUSION: The present work proposes a method to make in vivo quantitative assessment of proximal femur trabecular microstructure with a clinically practical scan duration of about 6 min.

Superior sagittal sinus flow as a proxy for tracking global cerebral blood flow dynamics during wakefulness and sleep.

Sleep, a state of reduced consciousness, affects brain oxygen metabolism and lowers cerebral metabolic rate of oxygen (CMRO2). Previously, we quantified CMRO2 during sleep via Fick's Principle, with a single-band MRI sequence measuring both hemoglobin O2 saturation (SvO2) and superior sagittal sinus (SSS) blood flow, which was upscaled to obtain total cerebral blood flow (tCBF). The procedure involves a brief initial calibration scan to determine the upscaling factor (fc), assumed state-invariant. Here, we used a dual-band sequence to simultaneously provide SvO2 in SSS and tCBF in the neck every 16 seconds, allowing quantification of fc dynamically. Ten healthy subjects were scanned by MRI with simultaneous EEG for 80 minutes, yielding 300 temporal image frames per subject. Four volunteers achieved slow-wave sleep (SWS), as evidenced by increased δ-wave activity (per American Academy of Sleep Medicine criteria). SWS was maintained for 13.5 ± 7.0 minutes, with CMRO2 28.6 ± 5.5% lower than pre-sleep wakefulness. Importantly, there was negligible bias between tCBF obtained by upscaling SSS-blood flow, and tCBF measured directly in the inflowing arteries of the neck (intra-class correlation 0.95 ± 0.04, averaged across all subjects), showing that the single-band approach is a valid substitute for quantifying tCBF, simplifying image data collection and analysis without sacrificing accuracy.