The effect of and correction for through-slice dephasing on 2D gradient-echo double angle B 1 + mapping.

PURPOSE: To show that B 0 $$ {\mathrm{B}}_0 $$ variations through slice and slice profile effects are two major confounders affecting 2D dual angle B 1 + $$ {\mathrm{B}}_1^{+} $$ maps using gradient-echo signals and thus need to be corrected to obtain accurate B 1 + $$ {\mathrm{B}}_1^{+} $$ maps. METHODS: The 2D gradient-echo transverse complex signal was Bloch-simulated and integrated across the slice dimension including nonlinear variations in B 0 $$ {\mathrm{B}}_0 $$ inhomogeneities through slice. A nonlinear least squares fit was used to find the B 1 + $$ {\mathrm{B}}_1^{+} $$ factor corresponding to the best match between the two gradient-echo signals experimental ratio and the Bloch-simulated ratio. The correction was validated in phantom and in vivo at 3T. RESULTS: For our RF excitation pulse, the error in the B 1 + $$ {\mathrm{B}}_1^{+} $$ factor scales by approximately 3.8% for every 10 Hz/cm variation in B 0 $$ {\mathrm{B}}_0 $$ along the slice direction. Higher accuracy phantom B 1 + $$ {\mathrm{B}}_1^{+} $$ maps were obtained after applying the proposed correction; the root mean square B 1 + $$ {\mathrm{B}}_1^{+} $$ error relative to the gold standard B 1 + $$ {\mathrm{B}}_1^{+} $$ decreased from 6.4% to 2.6%. In vivo whole-liver T 1 $$ {\mathrm{T}}_1 $$ maps using the corrected B 1 + $$ {\mathrm{B}}_1^{+} $$ map registered a significant decrease in T 1 $$ {\mathrm{T}}_1 $$ gradient through slice. CONCLUSION: B 0 $$ {\mathrm{B}}_0 $$ inhomogeneities varying through slice were seen to have an impact on the accuracy of 2D double angle B 1 + $$ {\mathrm{B}}_1^{+} $$ maps using gradient-echo sequences. Consideration of this confounder is crucial for research relying on accurate knowledge of the true excitation flip angles, as is the case of T 1 $$ {\mathrm{T}}_1 $$ mapping using a spoiled gradient recalled echo sequence.

Optimal flip angles for in vivo liver 3D T1 mapping and B1+ mapping at 3T.

PURPOSE: The spoiled gradient recalled echo (SPGR) sequence with variable flip angles (FAs) enables whole liver T 1 $$ {T}_1 $$ mapping at high spatial resolutions but is strongly affected by B 1 + $$ {B}_1^{+} $$ inhomogeneities. The aim of this work was to study how the precision of acquired T 1 $$ {T}_1 $$ maps is affected by the T 1 $$ {T}_1 $$ and B 1 + $$ {B}_1^{+} $$ ranges observed in the liver at 3T, as well as how noise propagates from the acquired signals into the resulting T 1 $$ {T}_1 $$ map. THEORY: The T 1 $$ {T}_1 $$ variance was estimated through the Fisher information matrix with a total noise variance including, for the first time, the B 1 + $$ {B}_1^{+} $$ map noise as well as contributions from the SPGR noise. METHODS: Simulations were used to find the optimal FAs for both the B 1 + $$ {B}_1^{+} $$ mapping and T 1 $$ {T}_1 $$ mapping. The simulations results were validated in 10 volunteers. RESULTS: Four optimized SPGR FAs of 2°, 2°, 15°, and 15° (TR = 4.1 ms) and B 1 + $$ {B}_1^{+} $$ map FAs of 65° and 130° achieved a T 1 $$ {T}_1 $$ coefficient of variation of 6.2 ± 1.7% across 10 volunteers and validated our theoretical model. Four optimal FAs outperformed five uniformly spaced FAs, saving the patient one breath-hold. For the liver B 1 + $$ {B}_1^{+} $$ and T 1 $$ {T}_1 $$ parameter space at 3T, a higher return in T 1 $$ {T}_1 $$ precision was obtained by investing FAs in the SPGR acquisition rather than in the B 1 + $$ {B}_1^{+} $$ map. CONCLUSION: A novel framework was developed and validated to calculate the SPGR T 1 $$ {T}_1 $$ variance. This framework efficiently identifies optimal FA values and determines the total number of SPGR and B 1 + $$ {B}_1^{+} $$ measurements needed to achieve a desired T 1 $$ {T}_1 $$ precision.

Hydration and glycogen affect T1 relaxation times of liver tissue.

T1 mapping is a useful tool for the assessment of patients with nonalcoholic fatty liver disease but still suffers from a large unexplained variance in healthy subjects. This study aims to characterize the potential effects of liver glycogen concentration and body hydration status on liver shortened modified Look-Locker inversion recovery (shMOLLI) T1 measurements. Eleven glycogen phantoms and 12 healthy volunteers (mean age: 31 years, three females) were scanned at 3 T using inversion recovery spin echo, multiple contrast spin echo (in phantoms), shMOLLI T1 mapping, multiple-echo spoiled gradient recalled echo and 13 C spectroscopy (in healthy volunteers). Phantom r1 and r2 relaxivities were determined from measured T1 and T2 values. Participants underwent a series of five metabolic experiments to vary their glycogen concentration and hydration levels: feeding, food fasting, exercising, underhydration, and rehydration. Descriptive statistics were calculated for shMOLLI T1 , inferior vena cava to aorta cross-sectional area ratio (IVC/Ao) as a marker of body hydration status, glycogen concentration, T2 * and proton density fat fraction values. A linear mixed model for shMOLLI R1 was constructed to determine the effects of glycogen concentration and IVC/Ao ratio. The mean shMOLLI T1 after fasting was 737 ± 67 ms. The mean within-subject change was 80 ± 45 ms. The linear mixed model revealed a glycogen r1 relaxivity in volunteers (0.18 M-1 s-1 , p = 0.03) close to that determined in phantoms (0.28 M-1 s-1 ). A unit change in IVC/Ao ratio was associated with a drop of -0.113 s-1 in R1 (p < 0.001). This study demonstrated a dependence of liver shMOLLI T1 values on liver glycogen concentration and overall body hydration status. Interparticipant variation of hydration status should be minimized in future liver MRI studies. Additionally, caution is advised when interpreting liver T1 measurements in participants with excess liver glycogen.

Navigator-based reacquisition and estimation of motion-corrupted data: Application to multi-echo spin echo for carotid wall MRI.

PURPOSE: To assess whether artifacts in multi-slice multi-echo spin echo neck imaging, thought to be caused by brief motion events such as swallowing, can be corrected by reacquiring corrupted central k-space data and estimating the remainder with parallel imaging. METHODS: A single phase-encode line (ky = 0, phase-encode direction anteroposterior) navigator echo was used to identify motion-corrupted data and guide the online reacquisition. If motion corruption was detected in the 7 central k-space lines, they were replaced with reacquired data. Subsequently, GRAPPA reconstruction was trained on the updated central portion of k-space and then used to estimate the remaining motion-corrupted k-space data from surrounding uncorrupted data. Similar compressed sensing-based approaches have been used previously to compensate for respiration in cardiac imaging. The g-factor noise amplification was calculated for the parallel imaging reconstruction of data acquired with a 10-channel neck coil. The method was assessed in scans with 9 volunteers and 12 patients. RESULTS: The g-factor analysis showed that GRAPPA reconstruction of 2 adjacent motion-corrupted lines causes high noise amplification; therefore, the number of 2-line estimations should be limited. In volunteer scans, median ghosting reduction of 24% was achieved with 2 adjacent motion-corrupted lines correction, and image quality was improved in 2 patient scans that had motion corruption close to the center of k-space. CONCLUSION: Motion-corrupted echo-trains can be identified with a navigator echo. Combined reacquisition and parallel imaging estimation reduced motion artifacts in multi-slice MESE when there were brief motion events, especially when motion corruption was close to the center of k-space.

Scattering matrix imaging pulse design for real-time respiration and cardiac motion monitoring.

PURPOSE: The scattering matrix (S-matrix) of a parallel transmit (pTx) coil is sensitive to physiological motion but requires additional monitoring RF pulses to be measured. In this work, we present and evaluate pTx RF pulse designs that simultaneously excite for imaging and measure the S-matrix to generate real-time motion signals without prolonging the image sequence. THEORY AND METHODS: Three pTx waveforms for measuring the S-matrix were identified and superimposed onto the imaging excitation RF pulses: (1) time division multiplexing, (2) frequency division multiplexing, and (3) code division multiplexing. These 3 methods were evaluated in healthy volunteers for scattering sensitivity and image artefacts. The S-matrix and real-time motion signals were calculated on the image calculation environment of the MR scanner. Prospective cardiac triggers were identified in early systole as a high rate of change of the cardiac motion signal. Monitoring accuracy was compared against electrocardiogram or the imaged diaphragm position. RESULTS: All 3 monitoring approaches measure the S-matrix during image excitation with quality correlated to input power. No image artefacts were observed for frequency multiplexing, and low energy artefacts were observed in the other methods. The accuracy of the achieved prospective cardiac gating was 15 ± 16 ms for breath hold and 24 ± 17 ms during free breathing. The diaphragm position prediction accuracy was 1.3 ± 0.9 mm. In all volunteers, good quality cine images were acquired for breath hold scans and dual gated CINEs were demonstrated. CONCLUSION: The S-matrix can be measured during image excitation to generate real-time cardiac and respiratory motion signals for prospective gating. No artefacts are introduced when frequency division multiplexing is used.

Feasibility of absolute quantification for 31 P MRS at 7 T.

PURPOSE: Phosphorus spectroscopy can differentiate among liver disease stages and types. To quantify absolute concentrations of phosphorus metabolites, sensitivity calibration and transmit field ( B1+ ) correction are required. The trend toward ultrahigh fields (7 T) and the use of multichannel RF coils makes this ever more challenging. We investigated the constraints on reference phantoms, and implemented techniques for the absolute quantification of human liver phosphorus spectra acquired using a 10-cm loop and a 16-channel array at 7 T. METHODS: The effect of phantom conductivity was assessed at 25.8 MHz (1.5 T), 49.9 MHz (3 T), and 120.3 MHz (7 T) by electromagnetic modeling. Radiofrequency field maps ( B1± ) were measured in phosphate phantoms (18 mM and 40 mM) at 7 T. These maps were used to assess the correction of 4 phantom 3D-CSI data sets using 3 techniques: phantom replacement, explicit normalization, and simplified normalization. In vivo liver spectra acquired with a 10-cm loop were corrected with all 3 methods. Simplified normalization was applied to in vivo 16-channel array data sets. RESULTS: Simulations show that quantification errors of less than 3% are achievable using a uniform electrolyte phantom with a conductivity of 0.23-0.86 S.m-1 at 1.5 T, 0.39-0.58 S.m-1 at 3 T, and 0.34-0.42 S.m-1 (16-19 mM KH2 PO4(aq) ) at 7 T. The mean γ-ATP concentration quantified in vivo at 7 T was 1.39 ± 0.30 mmol.L-1 to 1.71 ± 0.35 mmol.L-1 wet tissue for the 10-cm loop and 1.88 ± 0.25 mmol.L-1 wet tissue for the array. CONCLUSION: It is essential to select a calibration phantom with appropriate conductivity for quantitative phosphorus spectroscopy at 7 T. Using an 18-mM phosphate phantom and simplified normalization, human liver phosphate metabolite concentrations were successfully quantified at 7 T.

Magnitude-intrinsic water-fat ambiguity can be resolved with multipeak fat modeling and a multipoint search method.

PURPOSE: To develop a postprocessing algorithm for multiecho chemical-shift encoded water-fat separation that estimates proton density fat fraction (PDFF) maps over the full dynamic range (0-100%) using multipeak fat modeling and multipoint search optimization. To assess its accuracy, reproducibility, and agreement with state-of-the-art complex-based methods, and to evaluate its robustness to artefacts in abdominal PDFF maps. METHODS: We introduce MAGO (MAGnitude-Only), a magnitude-based reconstruction that embodies multipeak liver fat spectral modeling and multipoint optimization, and which is compatible with asymmetric echo acquisitions. MAGO is assessed first for accuracy and reproducibility on publicly available phantom data. Then, MAGO is applied to N = 178 UK Biobank cases, in which its liver PDFF measures are compared using Bland-Altman analysis with those from a version of the hybrid iterative decomposition of water and fat with echo asymmetry and least squares estimation (IDEAL) algorithm, LiverMultiScan IDEAL (LMS IDEAL, Perspectum Diagnostics Ltd, Oxford, UK). Finally, MAGO is tested on a succession of high field challenging cases for which LMS IDEAL generated artefacts in the PDFF maps. RESULTS: Phantom data showed accurate, reproducible MAGO PDFF values across manufacturers, field strengths, and acquisition protocols. Moreover, we report excellent agreement between MAGO and LMS IDEAL for 6-echo, 1.5 tesla human acquisitions (bias = -0.02% PDFF, 95% confidence interval = ±0.13% PDFF). When tested on 12-echo, 3 tesla cases from different manufacturers, MAGO was shown to be more robust to artefacts compared to LMS IDEAL. CONCLUSION: MAGO resolves the water-fat ambiguity over the entire fat fraction dynamic range without compromising accuracy, therefore enabling robust PDFF estimation where phase data is inaccessible or unreliable and complex-based and hybrid methods fail.

Mapping tissue water T1 in the liver using the MOLLI T1 method in the presence of fat, iron and B0 inhomogeneity.

Modified Look-Locker inversion recovery (MOLLI) T1 mapping sequences can be useful in cardiac and liver tissue characterization, but determining underlying water T1 is confounded by iron, fat and frequency offsets. This article proposes an algorithm that provides an independent water MOLLI T1 (referred to as on-resonance water T1 ) that would have been measured if a subject had no fat and normal iron, and imaging had been done on resonance. Fifteen NiCl2 -doped agar phantoms with different peanut oil concentrations and 30 adults with various liver diseases, nineteen (63.3%) with liver steatosis, were scanned at 3 T using the shortened MOLLI (shMOLLI) T1 mapping, multiple-echo spoiled gradient-recalled echo and 1 H MR spectroscopy sequences. An algorithm based on Bloch equations was built in MATLAB, and water shMOLLI T1 values of both phantoms and human participants were determined. The quality of the algorithm's result was assessed by Pearson's correlation coefficient between shMOLLI T1 values and spectroscopically determined T1 values of the water, and by linear regression analysis. Correlation between shMOLLI and spectroscopy-based T1 values increased, from r = 0.910 (P < 0.001) to r = 0.998 (P < 0.001) in phantoms and from r = 0.493 (for iron-only correction; P = 0.005) to r = 0.771 (for iron, fat and off-resonance correction; P < 0.001) in patients. Linear regression analysis revealed that the determined water shMOLLI T1 values in patients were independent of fat and iron. It can be concluded that determination of on-resonance water (sh)MOLLI T1 independent of fat, iron and macroscopic field inhomogeneities was possible in phantoms and human subjects.

Measuring inorganic phosphate and intracellular pH in the healthy and hypertrophic cardiomyopathy hearts by in vivo 7T 31P-cardiovascular magnetic resonance spectroscopy.

BACKGROUND: Cardiovascular phosphorus MR spectroscopy (31P-CMRS) is a powerful tool for probing energetics in the human heart, through quantification of phosphocreatine (PCr) to adenosine triphosphate (ATP) ratio. In principle, 31P-CMRS can also measure cardiac intracellular pH (pHi) and the free energy of ATP hydrolysis (ΔGATP). However, these require determination of the inorganic phosphate (Pi) signal frequency and amplitude that are currently not robustly accessible because blood signals often obscure the Pi resonance. Typical cardiac 31P-CMRS protocols use low (e.g. 30°) flip-angles and short repetition time (TR) to maximise signal-to-noise ratio (SNR) within hardware limits. Unfortunately, this causes saturation of Pi with negligible saturation of the flowing blood pool. We aimed to show that an adiabatic 90° excitation, long-TR, 7T 31P-CMRS protocol will reverse this balance, allowing robust cardiac pHi measurements in healthy subjects and patients with hypertrophic cardiomyopathy (HCM). METHODS: The cardiac Pi T1 was first measured by the dual TR technique in seven healthy subjects. Next, ten healthy subjects and three HCM patients were scanned with 7T 31P-MRS using long (6 s) TR protocol and adiabatic excitation. Spectra were fitted for cardiac metabolites including Pi. RESULTS: The measured Pi T1 was 5.0 ± 0.3 s in myocardium and 6.4 ± 0.6 s in skeletal muscle. Myocardial pH was 7.12 ± 0.04 and Pi/PCr ratio was 0.11 ± 0.02. The coefficients of repeatability were 0.052 for pH and 0.027 for Pi/PCr quantification. The pH in HCM patients did not differ (p = 0.508) from volunteers. However, Pi/PCr was higher (0.24 ± 0.09 vs. 0.11 ± 0.02; p = 0.001); Pi/ATP was higher (0.44 ± 0.14 vs. 0.24 ± 0.05; p = 0.002); and PCr/ATP was lower (1.78 ± 0.07 vs. 2.10 ± 0.20; p = 0.020), in HCM patients, which is in agreement with previous reports. CONCLUSION: A 7T 31P-CMRS protocol with adiabatic 90° excitation and long (6 s) TR gives sufficient SNR for Pi and low enough blood signal to permit robust quantification of cardiac Pi and hence pHi. Pi was detectable in every subject scanned for this study, both in healthy subjects and HCM patients. Cardiac pHi was unchanged in HCM patients, but both Pi/PCr and Pi/ATP increased that indicate an energetic impairment in HCM. This work provides a robust technique to quantify cardiac Pi and pHi.

Cardiac gating using scattering of an 8-channel parallel transmit coil at 7T.

PURPOSE: To establish a cardiac signal from scattering matrix or scattering coefficient measurements made on a 7T 8-channel parallel transmit (pTx) system, and to evaluate its use for cardiac gating. METHODS: Measurements of the scattering matrix and scattering coefficients were acquired using a monitoring pulse sequence and during a standard cine acquisition, respectively. Postprocessing used an independent component analysis and gating feature identification. The effect of the phase of the excitation radiofrequency (RF) field ( B1+ shim) on the cardiac signal was simulated for multiple B1+ shim configurations, and cine images were reconstructed from both the scattering coefficients and electrocardiogram (ECG). RESULTS: The cardiac motion signal was successfully identified in all subjects with a mean signal-to-noise ratio of 33.1 and 5.7 using the scattering matrix and scattering coefficient measurements, respectively. The dominant gating feature in the cardiac signal was a peak aligned with end-systole that occurred on average at 311 and 391 ms after the ECG trigger, with a mean standard deviation of 13.4 and 18.1 ms relative to ECG when using the scattering matrix and scattering coefficients measurements, respectively. The scattering coefficients showed a dependence on B1+ shim with some shim configurations not showing any cardiac signal. Cine images were successfully reconstructed using the scattering coefficients with minimal differences compared to those using ECG. CONCLUSION: We have shown that the scattering of a pTx RF coil can be used to estimate a cardiac signal, and that scattering matrix and coefficients can be used to cardiac gate MRI acquisitions with the scattering matrix providing a superior cardiac signal. Magn Reson Med 80:633-640, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Diaphragm position can be accurately estimated from the scattering of a parallel transmit RF coil at 7 T.

PURPOSE: To evaluate the use of radiofrequency scattering of a parallel transmit coil to track diaphragm motion. METHODS: Measurements made during radiofrequency excitation on an 8-channel parallel transmit coil by the directional couplers of the radiofrequency safety monitor were combined and converted into diaphragm position. A 30-s subject-specific calibration with an MRI navigator was used to determine a diaphragm estimate from each directional-coupler measure. Seven healthy volunteers were scanned at 7 T, in which images of the diaphragm were continuously acquired and directional couplers were monitored during excitation radiofrequency pulses. The ability to detect coughing was evaluated in one subject. The method was implemented on the scanner and evaluated for diaphragm gating of a free-breathing cardiac cine. RESULTS: Six of the seven scans were successful. In these subjects, the root mean square difference between MRI and scattering estimation of the superior-inferior diaphragm position was 1.4 ± 0.5 mm. On the scanner, the position was calculated less than 2 ms after every radiofrequency pulse. A prospectively gated (echocardiogram and respiration) high-resolution free-breathing cine showed no respiratory artifact and sharp blood-myocardium definition. CONCLUSIONS: Transmit coil scattering is sensitive to diaphragm motion and provides rapid, quantitative, and accurate monitoring of respiration. Magn Reson Med 79:2164-2169, 2018. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Hexagonal gradient scheme with RF spoiling improves spoiling performance for high-flip-angle fast gradient echo imaging.

PURPOSE: To present a framework in which time-varying gradients are applied with RF spoiling to reduce unwanted signal, particularly at high flip angles. METHODS: A time-varying gradient spoiler scheme compatible with RF spoiling is defined, in which spoiler gradients cycle through the vertices of a hexagon, which we call hexagonal spoiling. The method is compared with a traditional constant spoiling gradient both in the transition to and in the steady state. Extended phase graph (EPG) simulations, phantom acquisitions, and in vivo images were used to assess the method. RESULTS: Simulations, phantom and in vivo experiments showed that unwanted signal was markedly reduced by employing hexagonal spoiling, both in the transition to and in the steady state. For adipose tissue at 1.5 Tesla, the unwanted signal in the steady state with a 60 ° flip angle was reduced from 22% with constant spoiling to 2% with hexagonal spoiling. CONCLUSIONS: A time-varying gradient spoiler scheme that works with RF spoiling, called "hexagonal spoiling," has been presented and found to offer improved spoiling over the traditional constant spoiling gradient. Magn Reson Med 77:1231-1237, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Creatine kinase rate constant in the human heart measured with 3D-localization at 7 tesla.

PURPOSE: We present a new Bloch-Siegert four Angle Saturation Transfer (BOAST) method for measuring the creatine kinase (CK) first-order effective rate constant kf in human myocardium at 7 tesla (T). BOAST combines a variant of the four-angle saturation transfer (FAST) method using amplitude-modulated radiofrequency pulses, phosphorus Bloch-Siegert B1+-mapping to determine the per-voxel flip angles, and nonlinear fitting to Bloch simulations for postprocessing. METHODS: Optimal flip angles and repetition time parameters were determined from Monte Carlo simulations. BOAST was validated in the calf muscle of two volunteers at 3T and 7T. The myocardial CK forward rate constant was then measured in 10 volunteers at 7T in 82 min (after 1 H localization). RESULTS: BOAST kfCK values were 0.281 ± 0.002 s-1 in the calf and 0.35 ± 0.05 s-1 in myocardium. These are consistent with literature values from lower fields. Using a literature values for adenosine triphosphate concentration, we computed CK flux values of 4.55 ± 1.52 mmol kg-1 s-1 . The sensitive volume for BOAST depends on the B1 inhomogeneity of the transmit coil. CONCLUSION: BOAST enables measurement of the CK rate constant in the human heart at 7T, with spatial localization in three dimensions to 5.6 mL voxels, using a 10-cm loop coil. Magn Reson Med 78:20-32, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Simultaneous assessment of cardiac metabolism and perfusion using copolarized [1-13 C]pyruvate and 13 C-urea.

PURPOSE: To demonstrate the feasibility of imaging a bolus of co-polarized [1-13 C]pyruvate and 13 C-urea to simultaneously assess both metabolism and perfusion in the rodent heart. METHODS: Copolarized [1-13 C]pyruvate and 13 C-urea was imaged using a multi-echo, flow-sensitized spiral pulse sequence. Healthy rats were scanned in a two-factor factorial design (n = 12 total; metabolism: overnight fasting versus fed with dichloroacetate injection; perfusion: rest versus adenosine stress-induced hyperemia). RESULTS: Alterations in metabolism were detected by changes in pyruvate metabolism into 13 C-bicarbonate. Statistically independent alterations in perfusion were detected by changes in myocardial pyruvate and urea signals. CONCLUSION: The new pulse sequence was used to obtain maps of metabolism and perfusion in the rodent heart in a single acquisition. This hyperpolarized 13 C imaging test is expected to enable new studies in which the cardiac metabolism/perfusion mismatch can be studied in the acute environment. Magn Reson Med 77:151-158, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

A look-locker acquisition scheme for quantitative myocardial perfusion imaging with FAIR arterial spin labeling in humans at 3 tesla.

PURPOSE: A novel method for quantitative measurement of myocardial blood flow (MBF) using arterial spin labeling (ASL) in a single breath-hold is presented, evaluated by simulations, phantom studies and in vivo studies and tested for reproducibility and variability. METHODS: A flow-sensitive alternating inversion recovery (FAIR) ASL method with Look-Locker readout (LL-FAIR-ASL) was implemented at 3 tesla. Scans were performed on 10 healthy volunteers and MBF measured in three slices. The method was investigated for reproducibility by Bland-Altman analysis and statistical measures, the coefficients of reproducibility (CR) and variation (CV) are reported. RESULTS: The MBF values for the basal, mid, and apical slices were 1.04 ± 0.40, 1.06 ± 0.46, and 1.06 ± 0.38 ml/g/min, respectively (mean ± SD), which compare well with literature values. The CV across all scans, 43%, was greater than the between-session and within-session values, at 16 and 13%, respectively, for the mid-ventricular slice. The change in MBF required for detection, from the CR, was 61% between-session and 53% within-session for the mid-ventricle. CONCLUSION: This study shows the feasibility of the LL-FAIR-ASL method for the quantification of MBF. The statistical measures reported will allow the planning of future clinical research studies involving rest and stress measurements. Magn Reson Med 78:541-549, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Adiabatic excitation for 31 P MR spectroscopy in the human heart at 7 T: A feasibility study.

PURPOSE: Phosphorus magnetic resonance spectroscopy (31 P-MRS) provides a unique tool for assessing cardiac energy metabolism, often quantified using the phosphocreatine (PCr)/adenosine triphosphate (ATP) ratio. Surface coils are typically used for excitation for 31 P-MRS, but they create an inhomogeneous excitation field across the myocardium, producing undesirable, spatially varying partial saturation. Therefore, we implemented adiabatic excitation in a 3D chemical shift imaging (CSI) sequence for cardiac 31 P-MRS at 7 Tesla (T). METHODS: We optimized an adiabatic half passage pulse with bandwidth sufficient to excite PCr and γ-ATP together. In addition, the CSI sequence was modified to allow interleaved excitation of PCr and γ-ATP, then 2,3-DPG, to enable PCr/ATP determination with blood correction. Nine volunteers were scanned at 2 transmit voltages to confirm that measured PCr/ATP was independent of B1+ (i.e. over the adiabatic threshold). Six septal voxels were evaluated for each volunteer. RESULTS: Phantom experiments showed that adiabatic excitation can be reached at the depth of the heart using our pulse. The mean evaluated cardiac PCr/ATP ratio from all 9 volunteers corrected for blood signal was 2.14 ± 0.16. Comparing the two acquisitions with different voltages resulted in a minimal mean difference of -0.005. CONCLUSION: Adiabatic excitation is possible in the human heart at 7 T, and gives consistent PCr/ATP ratios. Magn Reson Med 78:1667-1673, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Phosphodiester content measured in human liver by in vivo 31 P MR spectroscopy at 7 tesla.

PURPOSE: Phosphorus (31 P) metabolites are emerging liver disease biomarkers. Of particular interest are phosphomonoester and phosphodiester (PDE) "peaks" that comprise multiple overlapping resonances in 31 P spectra. This study investigates the effect of improved spectral resolution at 7 Tesla (T) on quantifying hepatic metabolites in cirrhosis. METHODS: Five volunteers were scanned to determine metabolite T1 s. Ten volunteers and 11 patients with liver cirrhosis were scanned at 7T. Liver spectra were acquired in 28 min using a 16-channel 31 P array and 3D chemical shift imaging. Concentrations were calculated using γ-adenosine-triphosphate (γ-ATP) = 2.65 mmol/L wet tissue. RESULTS: T1 means ± standard deviations: phosphatidylcholine 1.05 ± 0.28 s, nicotinamide-adenine-dinucleotide (NAD+ ) 2.0 ± 1.0 s, uridine-diphosphoglucose (UDPG) 3.3 ± 1.4 s. Concentrations in healthy volunteers: α-ATP 2.74 ± 0.11 mmol/L wet tissue, inorganic phosphate 2.23 ± 0.20 mmol/L wet tissue, glycerophosphocholine 2.34 ± 0.46 mmol/L wet tissue, glycerophosphoethanolamine 1.50 ± 0.28 mmol/L wet tissue, phosphocholine 1.06 ± 0.16 mmol/L wet tissue, phosphoethanolamine 0.77 ± 0.14 mmol/L wet tissue, NAD+ 2.37 ± 0.14 mmol/L wet tissue, UDPG 2.00 ± 0.22 mmol/L wet tissue, phosphatidylcholine 1.38 ±â€Š0.31 mmol/L wet tissue. Inorganic phosphate and phosphatidylcholine concentrations were significantly lower in patients; glycerophosphoethanolamine concentrations were significantly higher (P < 0.05). CONCLUSION: We report human in vivo hepatic T1 s for phosphatidylcholine, NAD+ , and UDPG for the first time at 7T. Our protocol allows high signal-to-noise, repeatable measurement of metabolite concentrations in human liver. The splitting of PDE into its constituent peaks at 7T may allow more insight into changes in metabolism. Magn Reson Med 78:2095-2105, 2017. © 2017 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Quantitative ultrashort echo time imaging for assessment of massive iron overload at 1.5 and 3 Tesla.

PURPOSE: Hepatic iron content (HIC) quantification via transverse relaxation rate (R2*)-MRI using multi-gradient echo (mGRE) imaging is compromised toward high HIC or at higher fields due to the rapid signal decay. Our study aims at presenting an optimized 2D ultrashort echo time (UTE) sequence for R2* quantification to overcome these limitations. METHODS: Two-dimensional UTE imaging was realized via half-pulse excitation and radial center-out sampling. The sequence includes chemically selective saturation pulses to reduce streaking artifacts from subcutaneous fat, and spatial saturation (sSAT) bands to suppress out-of-slice signals. The sequence employs interleaved multi-echo readout trains to achieve dense temporal sampling of rapid signal decays. Evaluation was done at 1.5 Tesla (T) and 3T in phantoms, and clinical applicability was demonstrated in five patients with biopsy-confirmed massively high HIC levels (>25 mg Fe/g dry weight liver tissue). RESULTS: In phantoms, the sSAT pulses were found to remove out-of-slice contamination, and R2* results were in excellent agreement to reference mGRE R2* results (slope of linear regression: 1.02/1.00 for 1.5/3T). UTE-based R2* quantification in patients with massive iron overload proved successful at both field strengths and was consistent with biopsy HIC values. CONCLUSION: The UTE sequence provides a means to measure R2* in patients with massive iron overload, both at 1.5T and 3T. Magn Reson Med 78:1839-1851, 2017. © 2017 Wiley Periodicals, Inc.

A model for hepatic fibrosis: the competing effects of cell loss and iron on shortened modified Look-Locker inversion recovery T1 (shMOLLI-T1 ) in the liver.

PURPOSE: To propose a simple multicompartment model of the liver and use Bloch-McConnell simulations to demonstrate the effects of iron and fibrosis on shortened-MOLLI (shMOLLI) T1 measurements. Liver T1 values have shown sensitivity to inflammation and fibrosis, but are also affected by hepatic iron content. Modified Look-Locker inversion recovery (MOLLI) T1 measurements are biased by the lower T2 associated with high iron. MATERIALS AND METHODS: A tissue model was generated consisting of liver cells and extracellular fluid (ECF), with iron-dependent relaxation rates. Fibrosis was imitated by increasing the ECF proportion. Simulations of the shMOLLI sequence produced a look-up table (LUT) of shMOLLI-T1 for a given ECF fraction and iron content. The LUT was used to calculate ECF(shMOLLI-T1 ), assuming normal hepatic iron content (HIC), and ECF(shMOLLI- T1,T2*), accounting for HIC determined by T2*, for 77 patients and compared to fibrosis assessed by liver biopsy. RESULTS: Simulations showed that increasing HIC decreases shMOLLI-T1 , with an increase in HIC from 1.0 to 2.5 mg/g at normal ECF fraction decreasing shMOLLI-T1 by 160 msec, while increasing ECF increased ShMOLLI-T1 , with an increase of 20% ECF at normal iron increasing shMOLLI-T1 by 200 msec. Calculated patient ECF(shMOLLI-T1 ) showed a strong dependence on Ishak score (3.3 ± 0.8 %ECF/Ishak stage) and 1/T2* (-0.23 ± 0.04 %ECF/Hz). However, when iron was accounted for to produce ECF(shMOLLI- T1,T2*), it was independent of HIC but retained sensitivity to Ishak score. CONCLUSION: Use of this multicompartment model of the liver with Bloch-McConnell simulation should enable compensation of iron effects when using shMOLLI-T1 to assess fibrosis. LEVEL OF EVIDENCE: 1 J. Magn. Reson. Imaging 2017;45:450-462.