Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 07 May 2021
Sec. Medical Physics and Imaging

Feasibility of Hepatic Fat Quantification Using Proton Density Fat Fraction by Multi-Echo Chemical-Shift-Encoded MRI at 7T

\nRadim Ko&#x;ínekRadim Kořínek1Lorenz PflegerLorenz Pfleger2Korbinian EcksteinKorbinian Eckstein3Hannes BeiglbckHannes Beiglböck2Simon Daniel RobinsonSimon Daniel Robinson3Michael KrebsMichael Krebs2Siegfried Trattnig,Siegfried Trattnig3,4Zenon Star
uk Jr.Zenon Starčuk Jr.1Martin Kr&#x;&#x;k,,
Martin Krššák2,3,4*
  • 1Magnetic Resonance group, Institute of Scientific Instruments of the Czech Academy of Sciences, Brno, Czechia
  • 2Division of Endocrinology and Metabolism, Department of Medicine III, Medical University of Vienna, Vienna, Austria
  • 3Department of Biomedical Imaging and Image-Guided Therapy, High-Field Magnetic Resonance Centre, Medical University of Vienna, Vienna, Austria
  • 4Christian Doppler Laboratory for Clinical Molecular Imaging, CD Laboratory for Clinical Molecular MR Imaging (MOLIMA), Medical University of Vienna, Vienna, Austria

Fat fraction quantification and assessment of its distribution in the hepatic tissue become more important with the growing epidemic of obesity, and the increasing prevalence of diabetes mellitus type 2 and non-alcoholic fatty liver disease. At 3Tesla, the multi-echo, chemical-shift-encoded magnetic resonance imaging (CSE-MRI)-based acquisition allows the measurement of proton density fat-fraction (PDFF) even in clinical protocols. Further improvements in SNR can be achieved by the use of phased array coils and increased static magnetic field. The purpose of the study is to evaluate the feasibility of PDFF imaging using a multi-echo CSE-MRI technique at ultra-high magnetic field (7Tesla). Thirteen volunteers (M/F) with a broad range of age, body mass index, and hepatic PDFF were measured at 3 and 7T by multi-gradient-echo MRI and single-voxel spectroscopy MRS. All measurements were performed in breath-hold (exhalation); the MRI protocols were optimized for a short measurement time, thus minimizing motion-related problems. 7T data were processed off-line using Matlab® (MRI:multi-gradient-echo) and jMRUI (MRS), respectively. For quantitative validation of the PDFF results, a similar protocol was performed at 3T, including on-line data processing provided by the system manufacturer, and correlation analyses between 7 and 3T data were performed off-line. The multi-echo CSE-MRI measurements at 7T with a phased-array coil configuration and an optimal post-processing yielded liver volume coverage ranging from 30 to 90% for high- and low-BMI subjects, respectively. PDFFs ranged between 1 and 20%. We found significant correlations between 7T MRI and -MRS measurements (R2 ≅ 0.97; p < 0.005), and between MRI-PDFF at 7T and 3T fields (R2 ≅ 0.94; p < 0.005) in the evaluated volumes. Based on the measurements and analyses performed, the multi-echo CSE-MRI method using a 32-channel coil at 7T showed its aptitude for MRI-based quantitation of PDFF in the investigated volumes. The results are the first step toward qMRI of the whole liver at 7T with further improvements in hardware.

Introduction

The importance of fat fraction quantification and distribution in the human tissue is growing with the epidemic of obesity, and the increasing prevalence of diabetes mellitus type 2 (T2DM) [1, 2] and non-alcoholic fatty liver disease (NAFLD) [35].

Modern 3T clinical MR systems, providing high signal-to-noise ratio (SNR) and high resolution, allow the measurement of tissue proton density fat fraction (PDFF) [6, 7] even in routine clinical protocols. Multi-echo data is required to achieve accurate fat quantification based on Chemical Shift Encoded Magnetic Resonance Imaging (CSE-MRI) [813]. To determine of accurate PDFF distribution in the inner organs, high image resolution is required. Then, however, the SNR is reduced and subsequent data processing, in which the PDFF is calculated for each image voxel/pixel, is adversely affected. In principle, the SNR can be increased by the use of phased array RF (multi-channel) coils [1416] and/or increased static magnetic field [1719], both having some practical limits and incurring higher costs of hardware [20]. Recent installations of ultra-high-field (UHF) MR systems (7T or more), along with improvements in RF hardware and acquisition methods, have clearly demonstrated superior data quality for neuro- and musculoskeletal imaging [21, 22]. Nevertheless, problems remain; in particular, the abdominal region is affected by water-fat displacement and stronger susceptibility artifacts. Moreover, at high fields, the chemical shifts between the water peak and the multiple spectral peaks in the fat signal are increased. To minimize chemical shift displacement, strong gradients and large acquisition bandwidths must be used, causing SNR reduction and potentially substantial eddy currents [23, 24]. The other UHF effects are a prolongation of relaxation time T1 and a shortening of T2 and T2* relaxation times [2527]. An important parameter in UHF MR for patient safety is the specific absorption rate (SAR), which, in principle, increases quadratically with B0 field, but the specific spatial pattern depends in a complex way on the B1 frequency and a variety of factors, including the subject [28, 29]. These factors are a problem mainly with whole-body or large-volume coverage coils, such as birdcage, saddle, and TEM coils. A related UHF-MR problem is inhomogeneity of the excitation RF-field (B1+), which leads to inhomogeneous excitation and errors in quantitative imaging if not properly addressed. To overcome the problems of heterogeneous SAR, excitation, and detection sensitivity, dedicated multi-channel RF coils in transmit and receive modes [3033] have been used and field-specific adaptations to MR acquisition protocols have been implemented [34, 35].

CSE-MRI-based quantification of PDFF is a fast and reliable way to determine the distribution of fat in a tissue. This approach, proposed by Dixon [36] in 1984, has undergone considerable changes and evolution [8, 9, 13, 37, 38]. These various changes are frequently called “Dixon” methods even for multi-echo approaches; however, the term “CSE-MRI” is usually used for advanced methods that include a multi-echo acquisition scheme. Generally, the signal model representing the behavior of a vector of magnetization during the measurement sequence is crucial for the calculation of correct MRI-PDFF values. Due to the complex lipid spectrum [6, 39, 40], which contains several spectral lines in the frequency range of ~4.5 ppm (main peak at 1.3 ppm -CH2-), the signal model must contain prior knowledge about the spectral position and the relative intensity of each peak. A further common problem of all “Dixon”/CSE-MRI water-fat separation methods is field map (B0) estimation [11, 41]. At higher static magnetic fields, larger local magnetic field gradients result in more phase overflows of 2π in complex images compared to those at lower fields. The occurrence of these phase wraps may lead to water-fat swaps in the calculated MRI-PDFF maps, necessitating the application of phase-unwrapping algorithms. Several approaches have addressed this problem successfully [37, 41]. In hepatic iron overload, especially at excessive concentration, the tissue signal is significantly dephased, and hence, T2* relaxation is distinctly shortened; thus, the R2* = 1/T2* must be included in the signal model [8]; otherwise, substantial errors in PDFF estimation [42, 43] may occur.

In fact, the current clinical 3T protocols for PDFF measurements collect data in low resolution due to the requirement for a short measurement time (breath-hold) and a good SNR. The full coverage of the abdominal space is standard at 3T; however, the lower resolution smooths the information about the fat or iron distribution in the liver. This is a potential problem for the detailed study of fat distribution in the liver. Higher resolution that prolongs measurement time influences other measuring parameters (TE, TR, BW, and many others), and mainly deteriorates SNR in acquired data, which leads to noisy parametric maps (PDFF, R2*). Generally, quantitative magnetic resonance imaging (qMRI) is a crucial component of the many therapies and diagnostic [44], and provides the relatively stable and reproducible results [45]. The using of high field potentially yield benefit in the form of higher SNR compare to a low field, and it can improve the qMRI [46]. We have to note at the outset that due to incomparable coil configuration, the quality of images is not compared (SNR) in the study (it would not be objective); and the main focus is on the comparison of the quantitative results from both magnetic field, 3 and 7T. This study is the first step in exploring the possibilities of abdominal quantitative MRI (qMRI) that could provide improvement in diagnostic accuracy for a wide range of chronic liver diseases due to higher sensitivity at 7T in combination with the appropriate hardware equipment.

The purpose of this study was to assess the feasibility of proton density fat fraction (PDFF) quantification using multi-echo MRI at 7T with a 32-channel phase-array coil without B1 shimming. To demonstrate the potential value of the method at UHF with the best possible available hardware configuration at our institution (at the time of the study), PDFF derived from 7T MRI measurement (MRI-PDFF) was compared to PDFF determined by 7T MRS (MRS-PDFF) and to gold-standard [47, 48] multi-echo MRI-based 3T measurements [49, 50] on the same group of subjects.

Materials and Methods

Subjects

Thirteen subjects (4f/9m; age, 44.7 ± 14.7 years; body mass index (BMI) 25.6 ± 4.7 kg.m−2; mean ± SD) participated in this study. Volunteers were recruited based on the hepatic PDFF values obtained in previous studies [5153] to cover a broad range of PDFFs (0–20%) without a focus on their respective health status or diagnosis. The group comprised six lean volunteers (BMI = 21.7 ± 1.9 kg.m−2) and seven volunteers with high BMI (BMI = 29.6 ± 2.7 kg.m−2).

Ethics Statement

The study was approved by the local ethics committee. Informed consent was obtained from all individual participants included in the study.

Acquisition and Reconstruction

7T Measurements

At 7T (MAGNETOM, Siemens Healthineers, Erlangen, Germany), a phased array receive/transmit surface (32-channel) coil (Cardiac Transceiver Array RF Coil, MRI.TOOLS GmbH) was used. MR image data were acquired with an accelerated 3D-SPGR sequence with bipolar readout gradients with the following parameters: field of view (FOV) = 38.0 cm × 33.2 cm; acquisition bandwidth (BW) = 1,395 Hz/pixel; repetition time (TR) = 9.6 ms; flip angle (FA) = 4°; 32 slices (slice gap of 20%); acquisition matrix size in-plane = 256 × 224 pixels; voxel size 0.74 × 0.74 × 4 mm3; and six echoes with an equidistant echo spacing of ΔTE = 1.81 ms (the shortest possible), where the first TE = 1.45 ms, and acquisition time TA = 11.8 s, with a GRAPPA [54] acceleration factor of 8. Multi-channel multi-echo data were combined with the scanner image reconstructor using ASPIRE [55] and the PDFF maps were generated using the Graph-Cut approach [37] in a MATLAB® toolbox [56, 57], including the prior knowledge of the multi-frequency fat spectrum [39]. Single-voxel proton spectroscopic measurements were performed using a modified STEAM sequence [27] with TR = 5 s, with echo times TE = 6, 12, and 20 ms, TM = 10 ms, and a voxel size of 30 × 30 × 30 mm3. Due to a relatively narrow frequency bandwidth of the excitation pulse, the measurement was repeated with the same parameters, but with the excitation frequency (delta frequency) offset by −3.4 ppm from the water frequency (4.7 ppm) to fully cover the frequency band around the main fat resonances at 1.3 ppm (CH3). The acquired spectra were evaluated in jMRUI [58] with the AMARES fitting algorithm [59, 60] with a prior knowledge of the fat spectral components [39]. The MRS-based PDFF was calculated from T2-corrected spectra as the ratio of the estimated relative proton density of mobile lipids to the sum of the estimated relative proton densities of mobile water and mobile lipids.

3T Measurements

At 3T (Trio/PrismaFit, Siemens Healthineers, Erlangen, Germany), a combination of phased-array abdominal (18 channels) and spinal (32-channel) receiver coils and a whole-body transmit coil supplied by the MR-system manufacturer was used for data acquisition. MR image data were acquired by an accelerated 3D-SPGR (Spoiled Gradient Echo) [61] sequence with unipolar readout gradients with the following parameters: FOV = 38.0 cm × 31.4 cm; BW = 1.040 Hz/pixel; TR = 9.32 ms; FA = 3° (to minimize T1 effects); 48 slices (slice gap of 20%); acquisition matrix size in-plane = 160 × 104 pixels (interpolated to 320 × 264 pixels); voxel size 1.2 × 1.2 × 3.5 mm3; and six echoes with an equidistant echo spacing of ΔTE = 1.31 ms, where the first TE = 1.23 ms, TA = 6.9 s, and a CAIPIRINHA [62, 63] acceleration factor of 4 (2 × 2). The MRI protocol contains the online water-fat separation provided by the system manufacturer, which allows direct visualization of PDFF maps immediately after the acquisition. Those maps were used in our evaluation. The implemented online water-fat separation is certified for clinical use; therefore, such reconstructed data (PDFF maps) were used as the 3T MRI reference. Moreover, supplementary PDFF maps were reconstructed from 3T MRI data using the same approach as in 7T to allow a comparison of the same processing algorithm applied to data from different magnetic fields.

Volume of Interest Selection

The VOI for the 7T MRI-PDFF (similar or almost identical volumes compare to MRS) data analysis and 7T MRS data acquisition was placed in a homogeneous hepatic tissue carefully chosen to avoid contamination from liver vessels and subcutaneous tissue, but in the vicinity of the multi-channel coil, thus ensuring the best possible signal-to-noise. Selection of a VOI (similar size as 7T case) in the same position was attempted on the 3T PDFF maps.

Effective Liver Volume Coverage

To evaluate the range of coil combinations used at 7T, the term “effective liver volume coverage” was introduced where the liver volume at 7T was compared with the “true liver volume.” The “true liver volume” was estimated from 3T images for each subject where the full coverage of the abdominal region (not only the liver) was expected (100% coverage); the liver segmentation was performed manually. At 7T, the image noise thresholding segmentations on measured data (echo images of each subject for the longest TE) were performed to identify the background noise regions in the images, and these masks were applied to the measured data to achieve masked images. Subsequently, the liver was manually segmented from the masked images.

Statistical Analyses

The reconstructed 3T- and 7T-MRI-PDFF data where displayed in a box-and-whisker diagram to show the distribution of fat within the investigated VOIs for each subject. To prove the relations between 7T-MRI, 3T-MRI and 7T-MRS PDFF measurements, linear regression and Bland-Altman analyses were performed. All statistical tests were performed in MATLAB (MathWorks, Natick, MA, USA).

Results

Examples of the acquired MRI-PDFF volumes for two subjects with different body sizes and composition, and with high (S6) and low (S3) BMI, at 3 and 7T, are shown in Figure 1 (and related figures of PDFFs with corresponding 7T anatomical images in Supplementary Material 1). The position of the spectroscopic volume of interest (VOI) and the ROI used for the comparison is delineated for each subject by the white box in the respective PDFF map; the calculated field maps are shown in Supplementary Material 2. The effective liver volume coverage for MRI-PDFF showed 100% of the liver volume, as well as the whole abdomen for all subjects at 3T. In the case of the 32-channel RF coil used at 7T, the effective liver volume coverage varied from ~30 to 90% based on body size and composition. Liver coverage was highest in the subjects with low BMI and smaller body (torso) size (example in Figure 1 S13, Supplementary Figures 1, 3 in Supplementary Material 1). A case of a subject with very high BMI and non-optimal body size is shown in Figure 1 S6 and Supplementary Figures 1, 2 in Supplementary Material 1. The whiskers diagrams (Figure 2) show the distributions of MRI-PDFF within the VOI and the related MRS-PDFF values acquired from similar or almost identical volumes.

FIGURE 1
www.frontiersin.org

Figure 1. Positions of VOIs in the cranial section of the liver of subjects S6 and S13 with lowest and highest BMI (within the experimental group of our subjects), respectively. The investigated region is depicted by white solid boxes across the volume. The 7T data are noisy due to higher resolution and bandwidth per pixel compare to 3T.

FIGURE 2
www.frontiersin.org

Figure 2. The distributions of MRI-PDFF for individual subjects at 3T (yellow boxes) and 7T (blue boxes). The related 7T-MRS-PDFF values are depicted by black solid horizontal crosses. The red horizontal lines in the bars are the medians and the dark blue solid oblique crosses are the mean values of these distributions. The bottom and top of the boxes represent the 25th and 75th percentile (Q1 and Q3 quartiles) of MRI-PDFF distributions. The ends (black horizontal lines) of the extended whiskers (dashed black lines) define upper and lower extremes (± 2.7σ ≈ 99.3% coverage).

The correlation analyses for the data presented in Figure 2 are shown in Figure 3, where the 7T-MRI-PDFFs are compared with (Figure 3A) 3T-MRI, and (Figure 3B) 7T-MRS PDFF values. In the both cases A and B (Figure 3), high R2-values with low p-values, indicating strong agreement between the 3T and 7T measurements, were observed. The calculated R2-values (with p-values) of A and B (Figure 3) cases were 0.936 (p ≈ 2.75·10−7; 7T MRI vs. 3T MRI) and 0.970 (p ≈ 1.32·10−6; 7T MRI vs. 7T MRS), respectively; in addition, the slopes were calculated. Then, the Bland-Altman (BA) analyses of the previous (A) and (B) cases were also performed, which can be seen in Figure 4. In the first case, the BA plot shows a mean difference (bias) of ≈-0.60% and upper and lower confidence intervals (CI of 95%) of 0.58 and −1.79%, respectively. In the second case (Figure 4 on the right), there was a mean difference (bias) of ≈-0.54% and upper and lower confidence intervals (CI of 95%) of 0.52 and −1.60%, respectively. The extended boxplots and related analyses of PDFF maps that were calculated by the offline toolbox using Graph-Cut algorithm/approach from 3T-MRI data are given in Supplementary Material 3.

FIGURE 3
www.frontiersin.org

Figure 3. The correlation analysis of (A) 7T-MRI-PDFF vs. 3T-MRI-PDFF, and (B) 7T-MRI-PDFF vs. 7T-MRS-PDFF. The red line represents the linear fit of the values, and 95% confidence intervals for the slope of the lines are depicted by red dotted lines (upper and lower bounds). The black dashed line corresponds to a perfect match. *Its data point.

FIGURE 4
www.frontiersin.org

Figure 4. The Bland Altman analysis for (left) 3T and 7T MRI, and (right) 7T MRI and MRS measurements. The red line depicts the mean difference (bias) of all measurements, and the blue dotted lines represent the 95% confidence intervals for the mean.

In subject S2, the 7T MRS measurement was not performed because of health problems (cramps) not related to the measurement. The reconstructed (7T) MRI-PDFF maps (mainly the volume of interest) from S10 were affected by strong water-fat swaps that could not be removed by changing the parameters of the algorithm used, and therefore, the results were excluded. The 7T spectra from subjects S3 and S9 were not included in the analysis due to movement artifacts in the data.

Discussion

In our study, we tested the feasibility of MRI-PDFF assessment at 7T. Even-though MRS was considered for the “reference method” in hepatic fat content quantification, recent developments in MRI-PDFF at 3T have substantially improved the accuracy of the approach and made it the method of choice for hepatic fat content quantification [47, 48] in clinical settings.

The results depicted in Figure 2 demonstrate an agreement between the MRI-PDFF distributions from 3T and 7T measurements in the investigated volumes. In most cases, the spectroscopic values were within the 25th and the 75th percentile of MRI-PDFF. In only one case was the MRS value not in this interval (Figure 2), probably due to subtle motion of the subject between the MRS and MRI measurements. Another possible deteriorating effect is intra-voxel inhomogeneity, due to movement-related changes in local B0, which can lead to an improper estimation of PDFF. In 7T, MRI-DFF distributions show a larger range of values, and the number of outliers (red crosses, Figure 2) is higher than at 3T. This can be caused by the inhomogeneity of B1- over the FOV if the investigated volume is relatively distant from the coil surface, and hence, the water-fat separation process is more prone to errors due to lower SNR.

Correlation analyses confirm good agreement between MRI-PDFF 7T and 3T measurements in the investigated volumes. The level of agreement is attributable to measurement conditions, such as the measurement sequences (minimizing the acquisition time), patient measurement management (short time period between the 3T and 7T measurements–not more than 1 h), and exhalation breath-hold. From our practical experience, the bias in BA-plots appears relatively low, given that we are comparing measurements at different magnetic field. Although initial preliminary studies concerning hepatic liver fat quantitation by MRI at 7T have already been published [64, 65], in this study, a larger number of subjects was investigated, the distribution of MRI-PDFF in the investigated volumes was analyzed and compared to already established 7T MRS, and 3T MRI measurements were performed.

Limitations

At the 7T, due to the unavailability of whole-body 1H transmit RF hardware and B1 shimming equipment, precise positioning of the multichannel Tx/Rx surface coil was necessary, and the quality of data and liver volume coverage also depended on the subject/patient body size and composition. In subjects with a low BMI, the liver volume coverage was sufficient and water-fat separation provided PDFF maps without water/fat swaps within the liver tissue. In subjects with a high BMI and body size, the PDFF maps were affected by water-fat swaps and determining the optimal parameters for the water-fat separation process was more complex. Nevertheless, the Graph-Cut approach applied here is flexible enough with regard to the input parameters (such as the regularization parameter, spatial subsampling for field map estimation, the range of field map, and many others), and, in many cases, a sufficient solution was found (i.e., the PDFF maps of liver cross-section or at least the VOI not affected by water/fat swaps). Furthermore, it is important to note that the use of ASPIRE for the combination of the multi-channel phase data ensured the correct input data for the used GC approach.

At this moment, PDFF-MRI at 7T with the phased-array Tx/Rx coil used does not provide satisfactory coverage of the whole liver for patients with a higher BMI (Figure 1, Supplementary Materials 1, 2), and does not provide any significant advantage over PDFF-MRI at 3T. The effective liver coverage is based on the patient body composition and the resulting electrical properties of a measured subject. We have to admit that the effective liver coverage dropping to 30% in some cases is a significant drawback of our 7T configuration compared to current 3T measurements where full liver coverage is not an issue, and that appears to be the main limitation of the study. The other potential limitations are minimal echo time and echo spacing, but in fact it can be a problem especially in approaches where are estimated individual fat components [66, 67]. In our case, we have the prior knowledge of fat spectral model. However, there is strong potential to improve the quality of whole-liver imaging at 7T. Subjects with a high BMI frequently accompanied by oversized abdominal organs results in practical problems, such as B0 shimming, RF power settings, and the coverage of the whole liver volume. The quality of 7T data could be improved using interactive B1 shimming [68, 69], and liver coverage could be increased by the use of a volume body coil [70] in combination with an array coil [71], but neither of these was available for this study. Nevertheless, there are several possible pulse sequence options that may improve data quality in the future: Shorter echo spacing (use of shorter excitation pulses) and implementation of CAIPIRINHA acceleration can provide a significant reduction of acquisition time compared to GRAPPA, which was the only available parallel imaging option on our 7T system.

Conclusion

Our results confirm the feasibility of hepatic fat content quantification by MRI-PDFF based on a multi-gradient-echo acquisition method at ultra-high field (7T) using a 32-channel Tx/Rx array coil. In addition to a coil configuration and other hardware equipment, the success of fat quantification using MRI-PDFF is based on the water-fat separation algorithm, including prior knowledge of the fat spectral model. In the presence of rapid field changes at 7T, a robust solution that avoids or at least minimizes water/fat swaps in the reconstructed MRI-PDFF maps (most importantly in the regions of interest) is required. The advanced coil configuration with the further envisioned hardware improvement will provide the opening for further improvement of whole-abdomen imaging and liver fat quantification for patients with a higher BMI.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee, Medical University of Vienna. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

RK and MKrš designed the study, processed data, drafted the initial manuscript, reviewed and revised the manuscript, and coordinated and supervised data collection. RK, LP, HB, and MKre performed the measurements and ensured the patients for the study. KE and SR are responsible for ASPIRE reconstruction. MKrš, ZS, ST, and SR critically reviewed the manuscript for important intellectual content. All authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

Funding

This work was partially supported by the European Commission and the Ministry of Education, Youth, and Sports of the Czech Republic (projects Nos. CZ.02.1.01/0.0/0.0/16_013/0001775, LM2018129; 7AMB18AT023; 8J18AT023; AKTION AUT-CZE # 74p6), the Czech Academy of Sciences (project No. MSM100651801), Austrian Federal Ministry of Education, Science and Research; Contract Grant No. BMWFW WTZ Mobility, CZ09-2019 and the Austrian Science Fund FWF31452. SR was also supported by the Marie Skłodowska-Curie Action MS-fMRI-QSM 794298.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2021.665562/full#supplementary-material

References

1. Anderwald C, Bernroider E, Krssak M, Stingl H, Brehm A, Bischof MG, et al. Effects of insulin treatment in type 2 diabetic patients on intracellular lipid content in liver and skeletal muscle. Diabetes. (2002) 51:3025–32. doi: 10.2337/diabetes.51.10.3025

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Krssak M, Falk PK, Dresner A, DiPietro L, Vogel SM, Rothman DL, et al. Intramyocellular lipid concentrations are correlated with insulin sensitivity in humans: a 1H NMR spectroscopy study. Diabetologia. (1999) 42:113–6. doi: 10.1007/s001250051123

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Stepanova M, Younossi Z. Independent association between nonalcoholic fatty liver disease and cardiovascular disease in the US population. Clin Gastroenterol Hepatol. (2012) 10:646–50. doi: 10.1016/j.cgh.2011.12.039

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Seo SW, Gottesman RF, Clark JM, Hernaez R, Chang Y, Kim C, et al. Nonalcoholic fatty liver disease is associated with cognitive function in adults. Neurology. (2016) 86:1136–42. doi: 10.1212/WNL.0000000000002498

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Younossim Z. Non-alcoholic fatty liver disease - a global public health perspective. J Hepatol. (2019) 70:531–44. doi: 10.1016/j.jhep.2018.10.033

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Reeder S, Hu H, Sirlin C. Proton density fat-fraction: a standardized MR-based biomarker of tissue fat concentration. J Magn Reson Imaging. (2012) 36:1011–4. doi: 10.1002/jmri.23741

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Idilman IS, Aniktar H, Idilman R, Kabacam G, Savas B, Elhan A, et al. Hepatic steatosis: quantification by proton density fat fraction with MR imaging versus liver biopsy. Radiology. (2013) 267:767–75. doi: 10.1148/radiol.13121360

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Yu H, Shimakawa A, McKenzie C, Brodsky E, Brittain J, Reeder S. Multiecho water-fat separation and simultaneous R2* estimation with multifrequency fat spectrum modeling. Magn Reson Med. (2008) 60:1122–34. doi: 10.1002/mrm.21737

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Tsao J, Jiang Y. Hierarchical IDEAL: fast, robust, and multiresolution separation of multiple chemical species from multiple echo times. Magn Reson Med. (2013) 70:155–9. doi: 10.1002/mrm.24441

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Reeder S, Wen Z, Yu H, Pineda AR, Gold GE, Markl M, et al. Multicoil dixon chemical species separation with an iterative least-squares estimation method. Magn Reson Med. (2004) 51:35–45. doi: 10.1002/mrm.10675

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Glover G. Multipoint dixon technique for water and fat proton and susceptibility imaging. J Magn Reson Imaging. (1991) 1:521–30. doi: 10.1002/jmri.1880010504

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yu H, McKenzie C, Shimakawa A, Vu AT, Brau AC, Beatty PJ, et al. Multiecho reconstruction for simultaneous water-fat decomposition and T2* estimation. J Magn Reson Imaging. (2007) 26:1153–61. doi: 10.1002/jmri.21090

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Yu H, Shimakawa A, Hines CD, McKenzie CA, Hamilton G, Sirlin CB, et al. Combination of complex-based and magnitude-based multiecho water-fat separation for accurate quantification of fat-fraction. Magn Reson Med. (2011) 66:199–206. doi: 10.1002/mrm.22840

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Dietrich O, Raya J, Reeder S, Reiser M, Schoenberg S. Measurement of signal-to-noise ratios in MR images: influence of multichannel coils, parallel imaging, and reconstruction filters. J Magn Reson Imaging. (2007) 26:375–85. doi: 10.1002/jmri.20969

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Stockmann J, Witzel T, Keil B, Polimeni JR, Mareyam A, LaPierre C, et al. A 32-channel combined RF and B0 shim array for 3T brain imaging. Magn Reson Med. (2016) 75:441–51. doi: 10.1002/mrm.25587

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Navarro de Lara L, Windischberger C, Kuehne A, Woletz M, Sieg J, Bestmann S, et al. A novel coil array for combined TMS/fMRI experiments at 3 T. Magn Reson Med. (2015) 74:1492–501. doi: 10.1002/mrm.25535

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wen H, Denison T, Singerman R, Balaban R. The intrinsic signal-to-noise ratio in human cardiac imaging at 1.5, 3, and 4 T. J Magn Reson. (1997) 125:65–71. doi: 10.1006/jmre.1996.1072

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Vaughan JT, Garwood M, Collins CM, Liu W, DelaBarre L, Adriany G, et al. 7T vs. 4T: RF power, homogeneity, and signal-to-noise comparison in head images. Magn Reson Med. (2001) 46:24–30. doi: 10.1002/mrm.1156

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Pohmann R, Speck O, Scheffler K. Signal-to-noise ratio and MR tissue parameters in human brain imaging at 3, 7, and 9.4 tesla using current receive coil arrays. Magn Reson Med. (2016) 75:801–9. doi: 10.1002/mrm.25677

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Balchandani P, Naidich T. Ultra-High-Field MR Neuroimaging. AJNR Am J Neuroradiol. (2015) 36:1204–15. doi: 10.3174/ajnr.A4180

CrossRef Full Text | Google Scholar

21. Trattnig S, Springer E, Bogner W, Hangel G, Strasser B, Dymerska B, et al. Key clinical benefits of neuroimaging at 7T. Neuroimage. (2018) 168:477–89. doi: 10.1016/j.neuroimage.2016.11.031

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Trattnig S, Zbýn S, Schmitt B, Friedrich K, Juras V, Szomolanyi P, et al. Advanced MR methods at ultra-high field (7 Tesla) for clinical musculoskeletal applications. Eur Radiol. (2012) 22:2338–46. doi: 10.1007/s00330-012-2508-0

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Korinek R, Bartusek B, Starcuk Z. Water-Fat separation in rat by MRI at high field (9.4T). In: Measurement 2015 Proceedings of the 10th International Conference on Measurement. Smolenice (2015). p. 133–6.

Google Scholar

24. Kickler N, van der Zwaag W, Mekle R, Kober T, Marques JP, Krueger G, et al. Eddy current effects on a clinical 7T-68 cm bore scanner. MAGMA. (2010) 23:39–43. doi: 10.1007/s10334-009-0192-0

PubMed Abstract | CrossRef Full Text | Google Scholar

25. de Graaf RA, Brown PB, McIntyre S, Nixon TW, Behar KL, Rothman DL. High magnetic field water and metabolite proton T1 and T2 relaxation in rat brain in vivo. Magn Reson Med. (2006) 56:386–94. doi: 10.1002/mrm.20946

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Gajdošík M, Chmelík M, Just-Kukurová I, Bogner W, Valkovič L, Trattnig S, et al. In vivo relaxation behavior of liver compounds at 7 Tesla, measured by single-voxel proton MR spectroscopy. J Magn Reson Imaging. (2014) 40:1365–74. doi: 10.1002/jmri.24489

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Gajdošík M, Chadzynski G, Hangel G, Mlynárik V, Chmelík M, Valkovič L, et al. Ultrashort-TE stimulated echo acquisition mode (STEAM) improves the quantification of lipids and fatty acid chain unsaturation in the human liver at 7 T. NMR Biomed. (2015) 28:1283–93. doi: 10.1002/nbm.3382

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Hoult D, Phil D. Sensitivity and power deposition in a high-field imaging experiment. J Magn Reson Imaging. (2000) 12:46–67. doi: 10.1002/1522-2586(200007)12:1<46::AID-JMRI6>3.0.CO;2-D

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Collins C, Wang Z. Calculation of radiofrequency electromagnetic fields and their effects in MRI of human subjects. Magn Reson Med. (2011) 65:1470–82. doi: 10.1002/mrm.22845

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Hardy CJ, Cline HE, Giaquinto RO, Niendorf T, Grant AK, Sodickson DK. 32-element receiver-coil array for cardiac imaging. Magn Reson Med. (2006) 55:1142–9. doi: 10.1002/mrm.20870

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Keil B, Blau JN, Biber S, Hoecht P, Tountcheva V, Setsompop K, et al. A 64-channel 3T array coil for accelerated brain MRI. Magn Reson Med. (2013) 70:248–58. doi: 10.1002/mrm.24427

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zhu Y. Parallel excitation with an array of transmit coils. Magn Reson Med. (2004) 51:775–784. doi: 10.1002/mrm.20011

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Lattanzi R, Sodickson DK, Grant AK, Zhu Y. Electrodynamic constraints on homogeneity and radiofrequency power deposition in multiple coil excitations. Magn Reson Med. (2009) 61:315–334. doi: 10.1002/mrm.21782

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Zhang B, Seifert AC, Kim JW, Borrello J, Xu J. 7 Tesla 22-channel wrap-around coil array for cervical spinal cord and brainstem imaging. Magn Reson Med. (2017) 78:1623–34. doi: 10.1002/mrm.26538

CrossRef Full Text | Google Scholar

35. Rietsch SHG, Orzada S, Maderwald S, Brunheim S, Philips BWJ, Scheenen TWJ, et al. 7T ultra-high field body MR imaging with an 8-channel transmit/32-channel receive radiofrequency coil array. Med Phys. (2018) 45:2978–90. doi: 10.1002/mp.12931

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Dixon W. Simple proton spectroscopic imaging. Radiology. (1984) 153:189–94. doi: 10.1148/radiology.153.1.6089263

CrossRef Full Text | Google Scholar

37. Hernando D, Kellman P, Haldar J, Liang Z. Robust water/fat separation in the presence of large field inhomogeneities using a graph cut algorithm. Magn Reson Med. (2010) 63:79–90. doi: 10.1002/mrm.22177

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Sharma SD, Hu HH, Nayak KS. Accelerated water-fat imaging using restricted subspace field map estimation and compressed sensing. Magn Reson Med. (2012) 67:650–9. doi: 10.1002/mrm.23052

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Ren J, Dimitrov I, Sherry AD, Malloy CR. Composition of adipose tissue and marrow fat in humans by 1H NMR at 7 tesla. J Lipid Res. (2008) 49:2055–62. doi: 10.1194/jlr.D800010-JLR200

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Hamilton G, Smith DL, Bydder M, Nayak KS, Hu HH. MR properties of brown and white adipose tissues. J Magn Reson Imaging. (2011) 34:468–73. doi: 10.1002/jmri.22623

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Yu H, Reeder SB, Shimakawa A, Brittain JH, Pelc NJ. Field map estimation with a region growing scheme for iterative 3-point water-fat decomposition. Magn Reson Med. (2005) 54:1032–9. doi: 10.1002/mrm.20654

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Hussain HK, Chenevert TL, Londy FJ, Gulani V, Swanson SD, McKenna BJ, et al. Hepatic fat fraction: MR imaging for quantitative measurement and display–early experience. Radiology. (2005) 237:1048–55. doi: 10.1148/radiol.2373041639

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Bydder M, Yokoo T, Hamilton G, Middleton MS, Chavez AD, Schwimmer JB, et al. Relaxation effects in the quantification of fat using gradient echo imaging. Magn Reson Imaging. (2008) 26:347–59. doi: 10.1016/j.mri.2007.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Curtis WA, Fraum TJ, An H, Chen Y, Shetty AS, Fowler KJ. Quantitative MRI of diffuse liver disease: current applications and future directions. Radiology. (2019) 290:23–30. doi: 10.1148/radiol.2018172765

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Gracien RM, Maiworm M, Brüche N, Shrestha M, Nöth U, Hattingen E, et al. How stable is quantitative MRI? - Assessment of intra- and inter-scanner-model reproducibility using identical acquisition sequences and data analysis programs. Neuroimage. (2020) 207:116364. doi: 10.1016/j.neuroimage.2019.116364

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Phal PM, Usmanov A, Nesbit GM, Anderson JC, Spencer D, Wang P, et al. Qualitative comparison of 3-T and 1.5-T MRI in the evaluation of epilepsy. AJR Am J Roentgenol. (2008) 191:890–5. doi: 10.2214/AJR.07.3933

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Bohte AE, van Werven JR, Bipat S, Stoker J. The diagnostic accuracy of US, CT, MRI and 1H-MRS for the evaluation of hepatic steatosis compared with liver biopsy: a meta-analysis. Eur Radiol. (2011) 21:87–97. doi: 10.1007/s00330-010-1905-5

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Korínek R, Gajdošík M, Trattnig S, Starčuk Z Jr, Krššák M. Low-level fat fraction quantification at 3 T: comparative study of different tools for water-fat reconstruction and MR spectroscopy. MAGMA. (2020) 33:455–68. doi: 10.1007/s10334-020-00825-9

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Krssák M, Hofer H, Wrba F, Meyerspeer M, Brehm A, Lohninger A, et al. Non-invasive assessment of hepatic fat accumulation in chronic hepatitis C by 1H magnetic resonance spectroscopy. Eur J Radiol. (2010) 74:e60–6. doi: 10.1016/j.ejrad.2009.03.062

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Hájek M, Dezortová M, Wagnerová D, Skoch A, Voska L, Hejlová I, et al. MR spectroscopy as a tool for in vivo determination of steatosis in liver transplant recipients. MAGMA. (2011) 24:297–304. doi: 10.1007/s10334-011-0264-9

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Pfleger L, Gajdošík M, Wolf P, Smajis S, Fellinger P, Kuehne A, et al. Absolute quantification of phosphor-containing metabolites in the liver using 31 P MRSI and hepatic lipid volume correction at 7T suggests no dependence on body mass index or age. J Magn Reson Imaging. (2019) 49:597–607. doi: 10.1002/jmri.26225

CrossRef Full Text | Google Scholar

52. Wolf P, Fellinger P, Pfleger L, Smajis S, Beiglböck H, Gajdošík M, et al. Reduced hepatocellular lipid accumulation and energy metabolism in patients with long standing type 1 diabetes mellitus. Sci Rep. (2019) 9:2576. doi: 10.1038/s41598-019-39362-4

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Smajis S, Gajdošík M, Pfleger L, Traussnigg S, Kienbacher C, Halilbasic E, et al. Metabolic effects of a prolonged, very-high-dose dietary fructose challenge in healthy subjects. Am J Clin Nutr. (2020) 9:369–77. doi: 10.1093/ajcn/nqz271

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, Wang J, et al. Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn Reson Med. (2002) 47:1202–10. doi: 10.1002/mrm.10171

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Eckstein K, Dymerska B, Bachrata B, Bogner W, Poljanc K, Trattnig S, et al. Computationally efficient combination of multi-channel phase data from multi-echo acquisitions (ASPIRE). Magn Reson Med. (2018) 79:2996–3006. doi: 10.1002/mrm.26963

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Hu HH, Bornert P, Hernando D, Kellman P, Ma J, Reeder S, et al. ISMRM workshop on fat-water separation: insights, applications and progress in MRI. Magn Reson Med. (2012) 68:378–88. doi: 10.1002/mrm.24369

PubMed Abstract | CrossRef Full Text

57. WaterFat12 ISMRM toolbox. (2012). Available online at: https://www.ismrm.org/workshops/FatWater12/data.htm

58. Naressi A, Couturier C, Devos JM, Janssen M, Mangeat C, de Beer R, et al. Java-based graphical user interface for the MRUI quantitation package. MAGMA. (2001) 12:141–52. doi: 10.1007/BF02668096

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Vanhamme L, van den Boogaart A, Van Huffel S. Improved method for accurate and efficient quantification of MRS data with use of prior knowledge. J Magn Reson. (1997) 129:35–43. doi: 10.1006/jmre.1997.1244

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Graveron-Demilly D. Quantification in magnetic resonance spectroscopy based on semi-parametric approaches. MAGMA. (2014) 27:113–30. doi: 10.1007/s10334-013-0393-4

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Haase A, Frahm J, Matthaei K. FLASH imaging: rapid NMR imaging using low flip angle pulses. J Magn Reson. (1986) 67:258–66. doi: 10.1016/0022-2364(86)90433-6

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Breuer FA, Blaimer M, Heidemann RM, Mueller MF, Griswold MA, Jakob PM. Controlled aliasing in parallel imaging results in higher acceleration (CAIPIRINHA) for multi-slice imaging. Magn Reson Med. (2005) 53:684–91. doi: 10.1002/mrm.20401

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Breuer FA, Blaimer M, Mueller MF, Seiberlich N, Heidemann RM, Griswold MA, et al. Controlled aliasing in volumetric parallel imaging (2D CAIPIRINHA). Magn Reson Med. (2006) 55:549–56. doi: 10.1002/mrm.20787

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Korínek R, Eckstein K, Stračuk Z Jr, Trattnig S, Krššák M. Feasibility of abdominal quantitative imaging at 7T: pilot study. In: Proceeding International Society Magnetic Resonance Medicine. Montreal, QC (2019).

65. Williams H, Doran E, Bawden S, Mirfin C, Gowland P. Liver mDixon at 7T. In: Proceeding International Society Magnetic Resonance Medicine. Montreal, QC (2019).

66. Leporq B, Lambert SA, Ronot M, Boucenna I, Colinart P, Cauchy F, et al. Hepatic fat fraction and visceral adipose tissue fatty acid composition in mice: quantification with 7.0T MRI. Magn Reson Med. (2016) 76:510–8. doi: 10.1002/mrm.25895

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Berglund J, Ahlström H, Kullberg J. Model-based mapping of fat unsaturation and chain length by chemical shift imaging–phantom validation and in vivo feasibility. Magn Reson Med. (2012) 68:1815–27. doi: 10.1002/mrm.24196

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Metzger GJ, Auerbach EJ, Akgun C, Simonson J, Bi X, Ugurbil K, et al. Dynamically applied B1+ shimming solutions for non-contrast enhanced renal angiography at 7.0 Tesla. Magn Reson Med. (2013) 69:114–26. doi: 10.1002/mrm.24237

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Damen M, van Houtum Q, van Leeuwen M, Luijten P, Webb A, Klomp, et al. Quantitative T1 and T2 measurements of pancreas at 7 Tesla using a multi-transmit system. In: Proceeding International Society Magnetic Resonance Medicine (Honolulu, HI) (2017).

Google Scholar

70. Vaughan JT, Snyder CJ, DelaBarre LJ, Bolan PJ, Tian J, Bolinger L, et al. Whole-body imaging at 7T: preliminary results. Magn Reson Med. (2009) 61:244–8. doi: 10.1002/mrm.21751

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Vaughan J, Snyder C, Delabarre L, Tian J, Adriany G, Andersen P, et al. Clinical imaging at 7T with a 16 channel whole body coil and 32 receive channels. In: Proceeding International Society Magnetic Resonance Medicine. Honolulu (2009).

Keywords: CSE-MRI, ultra-high magnetic field, 7T, feasibility, liver, PDFF

Citation: Kořínek R, Pfleger L, Eckstein K, Beiglböck H, Robinson SD, Krebs M, Trattnig S, Starčuk Z Jr and Krššák M (2021) Feasibility of Hepatic Fat Quantification Using Proton Density Fat Fraction by Multi-Echo Chemical-Shift-Encoded MRI at 7T. Front. Phys. 9:665562. doi: 10.3389/fphy.2021.665562

Received: 08 February 2021; Accepted: 29 March 2021;
Published: 07 May 2021.

Edited by:

Nirbhay Yadav, Johns Hopkins University, United States

Reviewed by:

Simon Auguste Lambert, Université Claude Bernard Lyon 1, France
Christian Herbert Ziener, German Cancer Research Center (DKFZ), Germany

Copyright © 2021 Kořínek, Pfleger, Eckstein, Beiglböck, Robinson, Krebs, Trattnig, Starčuk and Krššák. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Martin Krššák, martin.krssak@meduniwien.ac.at

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.