Generalized parameter estimation in multi-echo gradient-echo-based chemical species separation

Accelerated liver water T1 mapping using single-shot continuous inversion-recovery spiral imaging

PURPOSE

Liver T1 mapping techniques typically require long breath holds or long scan time in free-breathing, need correction for B 1 + inhomogeneities and process composite (water and fat) signals. The purpose of this work is to accelerate the multi-slice acquisition of liver water selective T1 (wT1) mapping in a single breath hold, improving the k-space sampling efficiency.

METHODS

The proposed continuous inversion-recovery (IR) Look-Locker methodology combines a single-shot gradient echo spiral readout, Dixon processing and a dictionary-based analysis for liver wT1 mapping at 3 T. The sequence parameters were adapted to obtain short scan times. The influence of fat, B 1 + inhomogeneities and TE on the estimation of T1 was first assessed using simulations. The proposed method was then validated in a phantom and in 10 volunteers, comparing it with MRS and the modified Look-Locker inversion-recovery (MOLLI) method. Finally, the clinical feasibility was investigated by comparing wT1 maps with clinical scans in nine patients.

RESULTS

The phantom results are in good agreement with MRS. The proposed method encodes the IR-curve for the liver wT1 estimation, is minimally sensitive to B 1 + inhomogeneities and acquires one slice in 1.2 s. The volunteer results confirmed the multi-slice capability of the proposed method, acquiring nine slices in a breath hold of 11 s. The present work shows robustness to B 1 + inhomogeneities (wT 1 , No B 1 + = 1.07 wT 1 , B 1 + - 45.63 , R 2 = 0.99 ) , good repeatability (wT 1 , 2 ° = 1 . 0 wT 1 , 1 ° - 2.14 , R 2 = 0.96 ) and is in better agreement with MRS (wT 1 = 0.92 wT 1 MRS + 103.28 , R 2 = 0.38 ) than is MOLLI (wT 1 MOLLI = 0.76 wT 1 MRS + 254.43 , R 2 = 0.44 ) . The wT1 maps in patients captured diverse lesions, thus showing their clinical feasibility.

CONCLUSION

A single-shot spiral acquisition can be combined with a continuous IR Look-Locker method to perform rapid repeatable multi-slice liver water T1 mapping at a rate of 1.2 s per slice without a B 1 + map. The proposed method is suitable for nine-slice liver clinical applications acquired in a single breath hold of 11 s.

This Is What Metabolic Dysfunction-Associated Steatotic Liver Disease Looks Like: Potential of a Multiparametric MRI Protocol

Metabolic dysfunction-associated steatotic liver disease (MASLD) is a prevalent condition with a broad spectrum defined by liver biopsy. This gold standard method evaluates three features: steatosis, activity (ballooning and lobular inflammation), and fibrosis, attributing them to certain grades or stages using a semiquantitative scoring system. However, liver biopsy is subject to numerous restrictions, creating an unmet need for a reliable and reproducible method for MASLD assessment, grading, and staging. Noninvasive imaging modalities, such as magnetic resonance imaging (MRI), offer the potential to assess quantitative liver parameters. This review aims to provide an overview of the available MRI techniques for the three criteria evaluated individually by liver histology. Here, we discuss the possibility of combining multiple MRI parameters to replace liver biopsy with a holistic, multiparametric MRI protocol. In conclusion, the development and implementation of such an approach could significantly improve the diagnosis and management of MASLD, reducing the need for invasive procedures and paving the way for more personalized treatment strategies.

Comparative review of algorithms and methods for chemical-shift-encoded quantitative fat-water imaging

PURPOSE

To propose a standardized comparison between state-of-the-art open-source fat-water separation algorithms for proton density fat fraction (PDFF) andR 2 * $$ {R}_2^{\ast } $$ quantification using an open-source multi-language toolbox.

METHODS

Eight recent open-source fat-water separation algorithms were compared in silico, in vitro, and in vivo. Multi-echo data were synthesized with varying fat-fractions, B0 off-resonance, SNR and TEs. Experimental evaluation was conducted using calibrated fat-water phantoms acquired at 3T and multi-site open-source phantoms data. Algorithms' performances were observed on challenging in vivo datasets at 3T. Finally, reconstruction algorithms were investigated with different fat spectra to evaluate the importance of the fat model.

RESULTS

In silico and in vitro results proved most algorithms to be not sensitive to fat-water swaps andB 0 $$ {\mathrm{B}}_0 $$ offsets with five or more echoes. However, two methods remained inaccurate even with seven echoes and SNR = 50, and two other algorithms' precision depended on the echo spacing scheme (p < 0.05). The remaining four algorithms provided reliable performances with limits of agreement under 2% for PDFF and 6 s-1 forR 2 * $$ {R}_2^{\ast } $$ . The choice of fat spectrum model influenced quantification of PDFF mildly (<2% bias) and ofR 2 * $$ {R}_2^{\ast } $$ more severely, with errors up to 20 s-1 .

CONCLUSION

In promoting standardized comparisons of MRI-based fat and iron quantification using chemical-shift encoded multi-echo methods, this benchmark work has revealed some discrepancies between recent approaches for PDFF andR 2 * $$ {R}_2^{\ast } $$ mapping. Explicit choices and parameterization of the fat-water algorithm appear necessary for reproducibility. This open-source toolbox further enables the user to optimize acquisition parameters by predicting algorithms' margins of errors.

Robust breast quantitative susceptibility mapping in the presence of silicone

PURPOSE

To (a) develop a preconditioned water-fat-silicone total field inversion (wfsTFI) algorithm that directly estimates the susceptibility map from complex multi-echo data in the breast in the presence of silicone and to (b) evaluate the performance of wfsTFI for breast quantitative susceptibility mapping (QSM) in silico and in vivo in comparison with formerly proposed methods.

METHODS

Numerical simulations and in vivo multi-echo gradient echo breast measurements were performed to compare wfsTFI to a previously proposed field map-based linear total field inversion algorithm (lTFI) with and without the consideration of the chemical shift of silicone in the field map estimation step. Specifically, a simulation based on an in vivo scan and data from five patients were included in the analysis.

RESULTS

In the simulation, wfsTFI is able to significantly decrease the normalized root mean square error from lTFI without (4.46) and with (1.77) the consideration of the chemical shift of silicone to 0.68. Both the in silico and in vivo wfsTFI susceptibility maps show reduced shadowing artifacts in local tissue adjacent to silicone, reduced streaking artifacts and no erroneous single voxels of diamagnetic susceptibility in proximity to silicone.

CONCLUSION

The proposed wfsTFI method can automatically distinguish between subjects with and without silicone. Furthermore wfsTFI accounts for the presence of silicone in the QSM dipole inversion and allows for the robust estimation of susceptibility in proximity to silicone breast implants and hence allows the visualization of structures that would otherwise be dominated by artifacts on susceptibility maps.

On the water-fat in-phase assumption for quantitative susceptibility mapping

PURPOSE

To (a) define multi-peak fat model-based effective in-phase echo times for quantitative susceptibility mapping (QSM) in water-fat regions, (b) analyze the relationship between fat fraction, field map quantification bias and susceptibility bias, and (c) evaluate the susceptibility mapping performance of the proposed effective in-phase echoes in comparison to single-peak in-phase echoes and water-fat separation for regions where both water and fat are present.

METHODS

Effective multipeak in-phase echo times for a bone marrow and a liver fat spectral model were derived from a single voxel simulation. A Monte Carlo simulation was performed to assess the field map estimation error as a function of fat fraction for the different in-phase echoes. Additionally, a phantom scan and in vivo scans in the liver, spine, and breast were performed and evaluated with respect to quantification accuracy.

RESULTS

The use of single-peak in-phase echoes can introduce a worst-case susceptibility bias of 0.43 $$ 0.43 $$  ppm. The use of effective multipeak in-phase echoes shows a similar quantitative performance in the numerical simulation, the phantom and in all in vivo anatomies when compared to water-fat separation-based QSM.

CONCLUSION

QSM based on the proposed effective multipeak in-phase echoes can alleviate the quantification bias present in QSM based on single-peak in-phase echoes. When compared to water-fat separation-based QSM the proposed effective in-phase echo times achieve a similar quantitative performance while drastically reducing the computational expense for field map estimation.

On quantification errors of R 2 * $$ {R}_2^{\ast } $$ and proton density fat fraction mapping in trabecularized bone marrow in the static dephasing regime

PURPOSE

To study the effect of field inhomogeneity distributions in trabecularized bone regions on the gradient echo (GRE) signal with short TEs and to characterize quantification errors on R 2 * $$ {R}_2^{\ast } $$ and proton density fat fraction (PDFF) maps when using a water-fat model with an exponential R 2 * $$ {R}_2^{\ast } $$ decay model at short TEs.

METHODS

Field distortions were simulated based on a trabecular bone micro CT dataset. Simulations were performed for different bone volume fractions (BV/TV) and for different bone-fat composition values. A multi-TE UTE acquisition was developed to acquire multiple UTEs with random order to minimize eddy currents. The acquisition was validated in phantoms and applied in vivo in a volunteer's ankle and knee. Chemical shift encoded MRI (CSE-MRI) based on a Cartesian multi-TE GRE scan was acquired in the spine of patients with metastatic bone disease.

RESULTS

Simulations showed that signal deviations from the exponential signal decay at short TEs were more prominent for a higher BV/TV. UTE multi-TE measurements reproduced in vivo the simulation-based predicted behavior. In regions with high BV/TV, the presence of field inhomogeneities induced an R 2 * $$ {R}_2^{\ast } $$ underestimation in trabecularized bone marrow when using CSE-MRI at 3T with a short TE.

CONCLUSION

R 2 * $$ {R}_2^{\ast } $$ can be underestimated when using short TEs (<2 ms at 3 T) and a water-fat model with an exponential R 2 * $$ {R}_2^{\ast } $$ decay model in multi-echo GRE acquisitions of trabecularized bone marrow.

MR-based proton density fat fraction (PDFF) of the vertebral bone marrow differentiates between patients with and without osteoporotic vertebral fractures

The bone marrow proton density fat fraction (PDFF) assessed with MRI enables the differentiation between osteoporotic/osteopenic patients with and without vertebral fractures. Therefore, PDFF may be a potentially useful biomarker for bone fragility assessment.

INTRODUCTION

To evaluate whether magnetic resonance imaging (MRI)-based proton density fat fraction (PDFF) of vertebral bone marrow can differentiate between osteoporotic/osteopenic patients with and without vertebral fractures.

METHODS

Of the 52 study patients, 32 presented with vertebral fractures of the lumbar spine (66.4 ± 14.4 years, 62.5% women; acute low-energy osteoporotic/osteopenic vertebral fractures, N = 25; acute high-energy traumatic vertebral fractures, N = 7). These patients were frequency matched for age and sex to patients without vertebral fractures (N = 20, 69.3 ± 10.1 years, 70.0% women). Trabecular bone mineral density (BMD) values were derived from quantitative computed tomography. Chemical shift encoding-based water-fat MRI of the lumbar spine was performed, and PDFF maps were calculated. Associations between fracture status and PDFF were assessed using multivariable linear regression models.

RESULTS

Over all patients, mean PDFF and trabecular BMD correlated significantly (r =  - 0.51, P < 0.001). In the osteoporotic/osteopenic group, those patients with osteoporotic/osteopenic fractures had a significantly higher PDFF than those without osteoporotic fractures after adjusting for age, sex, weight, height, and trabecular BMD (adjusted mean difference [95% confidence interval], 20.8% [10.4%, 30.7%]; P < 0.001), although trabecular BMD values showed no significant difference between the subgroups (P = 0.63). For the differentiation of patients with and without vertebral fractures in the osteoporotic/osteopenic subgroup using mean PDFF, an area under the receiver operating characteristic (ROC) curve (AUC) of 0.88 (P = 0.006) was assessed. When evaluating all patients with vertebral fractures, those with high-energy traumatic fractures had a significantly lower PDFF than those with low-energy osteoporotic/osteopenic vertebral fractures (P < 0.001).

CONCLUSION

MR-based PDFF enables the differentiation between osteoporotic/osteopenic patients with and without vertebral fractures, suggesting the use of PDFF as a potential biomarker for bone fragility.

Spatially resolved free-induction decay spectroscopy using a 3D ultra-short echo time multi-echo imaging sequence with systematic echo shifting and compensation of B0 field drifts

PURPOSE

Biologically interesting signals can exhibit fast transverse relaxation and frequency shifts compared to free water. For spectral assignment, a ultra-short echo time (UTE) imaging sequence was modified to provide pixel-wise free-induction decay (FID) acquisition.

METHODS

The UTE-FID approach presented relies on a multi-echo 3D spiral UTE sequence with six echoes per radiofrequency (RF) excitation (TE_min 0.05 ms, echo spacing 3 ms). A complex pixel-wise raw data set for FID spectroscopy is obtained by several multi-echo UTE measurements with systematic shifting of the readout by 0.25 or 0.5 ms, until the time domain is filled for 18 or 45 ms. B₀ drifts are compensated by mapping and according phase correction. Autoregressive extrapolation of the signal is performed before Gaussian filtering. This method was applied to a phantom containing collagen-water solutions of different concentrations. To calculate the collagen content, a 19-peak collagen model was extracted from a non-selective FID spectrum (50% collagen solution). Proton-density-collagen-fraction (PDCF) was calculated for 10 collagen solutions (2%-50%). Furthermore, an in vivo UTE-FID spectrum of adipose tissue was recorded.

RESULTS

UTE-FID signal patterns agreed well with the non-spatially selective pulse-acquire FID spectrum from a sphere filled with 50% collagen. Differentiation of collagen solution from distilled water in the PDCF map was possible from 4% collagen concentration for a UTE-FID sequence with 128 × 128 × 64 matrix (voxel size 1 × 1 × 2.85 mm³ ). The mean values of the PDCF correlate linearly with collagen concentration.

CONCLUSION

The presented UTE-FID approach allows pixel-wise raw data acquisition similar to non-spatially selective pulse-acquire spectroscopy. Spatially resolved applications for assessment of spectra of rapidly decaying signals seem feasible.

Assessment of vertebral fractures and edema of the thoracolumbar spine based on water-fat and susceptibility-weighted images derived from a single ultra-short echo time scan

PURPOSE

To develop a methodology to simultaneously perform single echo Dixon water-fat imaging and susceptibility-weighted imaging (SWI) based on a single echo time (TE) ultra-short echo time (UTE) (sUTE) scan to assess vertebral fractures and degenerative bone changes in the thoracolumbar spine.

METHODS

A methodology was developed to solve the smoothness-constrained inverse water-fat problem to separate water and fat while removing unwanted low-frequency phase terms. Additionally, the corrected UTE phase was used for SWI. UTE imaging (TE: 0.14 ms, 3T MRI) was performed in the lumbar spine of nine patients with vertebral fractures and bone marrow edema (BME). All images were reviewed by two radiologists. Water- and fat-separated images were analyzed in comparison with short-tau inversion recovery (STIR) and with respect to BME visibility. The visibility of fracture lines and cortical outlining of the UTE magnitude images were analyzed in comparison with computed tomography.

RESULTS

Unwanted phase components, dominated by the B₁ phase, were removed from the UTE phase images. The rating of the diagnostic quality of BME visualization showed a high preference for the sUTE-Dixon water- and fat-separated images in comparison with STIR. The UTE magnitude images enabled better visualizing fracture lines compared with STIR and slightly better visibility of cortical outlining. With increasing SWI weighting osseous structures and fatty tissues were enhanced.

CONCLUSION

The proposed sUTE-Dixon-SWI methodology allows the removal of unwanted low-frequency phases and enables water-fat separation and SWI processing from a single complex UTE image. The methodology can be used for the simultaneous assessment of vertebral fractures and BME of the thoracolumbar spine.

Vertebral bone marrow T2* mapping using chemical shift encoding-based water-fat separation in the quantitative analysis of lumbar osteoporosis and osteoporotic fractures

BACKGROUND

Chemical shift encoding-based water-fat separation techniques have been used for fat quantification [proton density fat fraction (PDFF)], but they also enable the assessment of bone marrow T2*, which has previously been reported to be a potential biomarker for osteoporosis and may give insight into the cause of vertebral fractures (i.e., osteoporotic vs. traumatic) and the microstructure of the bone when applied to vertebral bone marrow.

METHODS

The 32 patients (78.1% with low-energy osteopenic/osteoporotic fractures, mean age 72.3±9.8 years, 76% women; 21.9% with high-energy traumatic fractures, 47.3±12.8 years, no women) were frequency-matched for age and sex to subjects without vertebral fractures (n=20). All study patients underwent 3T-MRI of the lumbar spine including sagittally acquired spoiled gradient echo sequences for chemical shift encoding-based water-fat separation, from which T2* values were obtained. Volumetric trabecular bone mineral density (BMD) and trabecular bone parameters describing the three-dimensional structural integrity of trabecular bone were derived from quantitative CT. Associations between T2* measurements, fracture status and trabecular bone parameters were assessed using multivariable linear regression models.

RESULTS

Mean T2* values of non fractured vertebrae in all patients showed a significant correlation with BMD (r=-0.65, P<0.001), trabecular number (TbN) (r=-0.56, P<0.001) and trabecular spacing (TbSp) (r=0.61, P<0.001); patients with low-energy osteoporotic vertebral fractures showed significantly higher mean T2* values than those with traumatic fractures (13.6±4.3 vs. 8.4±2.2 ms, P=0.01) as well as a significantly lower TbN (0.69±0.08 vs. 0.93±0.03 mm-1, P<0.01) and a significantly larger trabecular spacing (1.06±0.16 vs. 0.56±0.08 mm, P<0.01). Mean T2* values of osteoporotic patients with and without vertebral fracture showed no significant difference (13.5±3.4 vs. 15.6±3.5 ms, P=0.40). When comparing the mean T2* of the fractured vertebrae, no significant difference could be detected between low-energy osteoporotic fractures and high-energy traumatic fractures (12.6±5.4 vs. 8.1±2.4 ms, P=0.10).

CONCLUSIONS

T2* mapping of vertebral bone marrow using using chemical shift encoding-based water-fat separation allows for assessing osteoporosis as well as the trabecular microstructure and enables a radiation-free differentiation between patients with low-energy osteoporotic and high-energy traumatic vertebral fractures, suggesting its potential as a biomarker for bone fragility.

Magnetic Resonance Imaging Techniques for Brown Adipose Tissue Detection

Magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) methods can non-invasively assess brown adipose tissue (BAT) structure and function. Recently, MRI and MRS have been proposed as a means to differentiate BAT from white adipose tissue (WAT) and to extract morphological and functional information on BAT inaccessible by other means. Specifically, proton MR (¹H) techniques, such as proton density fat fraction mapping, diffusion imaging, and intermolecular multiple quantum coherence imaging, have been employed to access BAT microstructure; MR thermometry, relaxometry, and MRI and MRS with ³¹P, ²H, ¹³C, and ¹²⁹Xe have shown to provide complementary information on BAT function. The purpose of the present review is to provide a comprehensive overview of MR imaging and spectroscopy techniques used to detect BAT in rodents and in humans. The present work discusses common challenges of current methods and provides an outlook on possible future directions of using MRI and MRS in BAT studies.