Repeat it without me: Crowdsourcing the T1 mapping common ground via the ISMRM reproducibility challenge.

PURPOSE: T1 mapping is a widely used quantitative MRI technique, but its tissue-specific values remain inconsistent across protocols, sites, and vendors. The ISMRM Reproducible Research and Quantitative MR study groups jointly launched a challenge to assess the reproducibility of a well-established inversion-recovery T1 mapping technique, using acquisition details from a seminal T1 mapping paper on a standardized phantom and in human brains. METHODS: The challenge used the acquisition protocol from Barral et al. (2010). Researchers collected T1 mapping data on the ISMRM/NIST phantom and/or in human brains. Data submission, pipeline development, and analysis were conducted using open-source platforms. Intersubmission and intrasubmission comparisons were performed. RESULTS: Eighteen submissions (39 phantom and 56 human datasets) on scanners by three MRI vendors were collected at 3 T (except one, at 0.35 T). The mean coefficient of variation was 6.1% for intersubmission phantom measurements, and 2.9% for intrasubmission measurements. For humans, the intersubmission/intrasubmission coefficient of variation was 5.9/3.2% in the genu and 16/6.9% in the cortex. An interactive dashboard for data visualization was also developed: https://rrsg2020.dashboards.neurolibre.org. CONCLUSION: The T1 intersubmission variability was twice as high as the intrasubmission variability in both phantoms and human brains, indicating that the acquisition details in the original paper were insufficient to reproduce a quantitative MRI protocol. This study reports the inherent uncertainty in T1 measures across independent research groups, bringing us one step closer to a practical clinical baseline of T1 variations in vivo.

Motion-induced phase-corrected homodyne reconstruction for partial Fourier single-shot diffusion-weighted echo planar imaging of the liver.

Partial Fourier encoding is popular in single-shot (ss) diffusion-weighted (DW) echo planar imaging (EPI) because it enables a shorter echo time (TE) and, hence, improves the signal-to-noise-ratio. Motion during diffusion encoding causes k-space shifting and dispersion, which compromises the quality of the homodyne reconstruction. This work provides a comprehensive understanding of the artifacts in homodyne reconstruction of partial Fourier ss-DW-EPI data in the presence of motion-induced phase and proposes the motion-induced phase-corrected homodyne (mpc-hdyne) reconstruction method to ameliorate these artifacts. Simulations with different types of motion-induced phase were performed to provide an understanding of the potential artifacts that occur in the homodyne reconstruction of partial Fourier ss-DW-EPI data. To correct for the artifacts, the mpc-hdyne reconstruction is proposed. The algorithm recenters k-space, updates the partial Fourier factor according to detected global k-space shifts, and removes low-resolution nonlinear phase before the conventional homodyne reconstruction. The mpc-hdyne reconstruction is tested on both simulation and in vivo data. Motion-induced phase can cause signal overestimation, worm artifacts, and signal loss in partial Fourier ss-DW-EPI data with the conventional homodyne reconstruction. Simulation and in vivo data showed that the proposed mpc-hdyne reconstruction ameliorated artifacts, yielding higher quality DW images compared with conventional homodyne reconstruction. Based on the understanding of the artifacts in homodyne reconstruction of partial Fourier ss-DW-EPI data, the mpc-hdyne reconstruction was proposed and showed superior performance compared with the conventional homodyne reconstruction on both simulation and in vivo data.

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.

Partial Fourier in the presence of respiratory motion in prostate diffusion-weighted echo planar imaging.

PURPOSE: To investigate the effect of respiratory motion in terms of signal loss in prostate diffusion-weighted imaging (DWI), and to evaluate the usage of partial Fourier in a free-breathing protocol in a clinically relevant b-value range using both single-shot and multi-shot acquisitions. METHODS: A controlled breathing DWI acquisition was first employed at 3 T to measure signal loss from deep breathing patterns. Single-shot and multi-shot (2-shot) acquisitions without partial Fourier (no pF) and with partial Fourier (pF) factors of 0.75 and 0.65 were employed in a free-breathing protocol. The apparent SNR and ADC values were evaluated in 10 healthy subjects to measure if low pF factors caused low apparent SNR or overestimated ADC. RESULTS: Controlled breathing experiments showed a difference in signal coefficient of variation between shallow and deep breathing. In free-breathing single-shot acquisitions, the pF 0.65 scan showed a significantly (p < 0.05) higher apparent SNR than pF 0.75 and no pF in the peripheral zone (PZ) of the prostate. In the multi-shot acquisitions in the PZ, pF 0.75 had a significantly higher apparent SNR than 0.65 pF and no pF. The single-shot pF 0.65 scan had a significantly lower ADC than single-shot no pF. CONCLUSION: Deep breathing patterns can cause intravoxel dephasing in prostate DWI. For single-shot acquisitions at a b-value of 800 s/mm2, any potential risks of motion-related artefacts at low pF factors (pF 0.65) were outweighed by the increase in signal from a lower TE, as shown by the increase in apparent SNR. In multi-shot acquisitions however, the minimum pF factor should be larger, as shown by the lower apparent SNR at low pF factors.

Silicone implant and fibrous capsule assessment based on water-fat-silicone chemical shift encoding-based species separation in breast MRI.

BACKGROUND: With rising breast augmentations worldwide, there is an increasing clinical need for an early and accurate detection of implant complications. PURPOSE: To compare the quality of chemical shift encoding-based (CSE) water-fat-silicone separation compared to double inversion recovery (DIR) silicone-only imaging in breast magnetic resonance imaging (MRI). MATERIAL AND METHODS: This retrospective, single-center study included women with silicone implants subjected to 3-T MRI between January 2021 and March 2022. MRI included (i) two-dimensional silicone-only T2-weighted turbo spin echo DIR acquisition and (ii) three-dimensional CSE imaging based on multi-echo gradient-echo sequence enabling water-, fat-, and silicone-image separation. Images were evaluated and compared by three independent radiologists using a clinically established rating including differentiability of the silicone implant, visibility and contouring of the adjacent fibrous capsule, and accuracy of intralesional folds in a ranking of 1-5. The apparent contrast-to-noise (aCNR) was calculated. RESULTS: In 71 women, the average quality of water-fat-silicone images from CSE imaging was assessed as "good" (assessment 4 ± 0.9). In 68 (96%) patients, CSE imaging achieved a concise delineation of the silicone implant and precise visualization of the fibrous capsule that was not distinguishable in DIR imaging. Implant ruptures were more easily detected in CSE imaging. The aCNR was higher in CSE compared to DIR imaging (18.43 ± 9.8 vs. 14.73 ± 2.5; P = 0.002). CONCLUSION: Intrinsically co-registered water-fat-silicone-separated CSE-based images enable a reliable assessment of silicone implants. The simultaneously improved differentiability of the implant and fibrous capsule may provide clinicians with a valuable tool for an accurate evaluation of implant integrity and early detection of potential complications.

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.

Motion compensated renal diffusion weighted imaging.

PURPOSE: To assess the effect of respiratory motion and cardiac driven pulsation in renal DWI and to examine asymmetrical velocity-compensated diffusion encoding waveforms for robust ADC mapping in the kidneys. METHODS: The standard monopolar Stejskal-Tanner pulsed gradient spin echo (pgse) and the asymmetric bipolar velocity-compensated (asym-vc) diffusion encoding waveforms were used for coronal renal DWI at 3T. The robustness of the ADC quantification in the kidneys was tested with the aforementioned waveforms in respiratory-triggered and breath-held cardiac-triggered scans at different trigger delays in 10 healthy subjects. RESULTS: The pgse waveform showed higher ADC values in the right kidney at short trigger delays in comparison to longer trigger delays in the respiratory triggered scans when the diffusion gradient was applied in the feet-head (FH) direction. The coefficient of variation over all respiratory trigger delays, averaged over all subjects was 0.15 for the pgse waveform in the right kidney when diffusion was measured in the FH direction; the corresponding coefficient of variation for the asym-vc waveform was 0.06. The effect of cardiac driven pulsation was found to be small in comparison to the effect of respiratory motion. CONCLUSION: Short trigger delays in respiratory-triggered scans can cause higher ADC values in comparison to longer trigger delays in renal DWI, especially in the right kidney when diffusion is measured in the FH direction. The asym-vc waveform can reduce ADC variation due to respiratory motion in respiratory-triggered scans at the cost of reduced SNR compared to the pgse waveform.

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.

Assessing breast density using the chemical-shift encoding-based proton density fat fraction in 3-T MRI.

OBJECTIVES: There is a clinical need for a non-ionizing, quantitative assessment of breast density, as one of the strongest independent risk factors for breast cancer. This study aims to establish proton density fat fraction (PDFF) as a quantitative biomarker for fat tissue concentration in breast MRI and correlate mean breast PDFF to mammography. METHODS: In this retrospective study, 193 women were routinely subjected to 3-T MRI using a six-echo chemical shift encoding-based water-fat sequence. Water-fat separation was based on a signal model accounting for a single T2* decay and a pre-calibrated 7-peak fat spectrum resulting in volumetric fat-only, water-only images, PDFF- and T2*-values. After semi-automated breast segmentation, PDFF and T2* values were determined for the entire breast and fibroglandular tissue. The mammographic and MRI-based breast density was classified by visual estimation using the American College of Radiology Breast Imaging Reporting and Data System categories (ACR A-D). RESULTS: The PDFF negatively correlated with mammographic and MRI breast density measurements (Spearman rho: -0.74, p < .001) and revealed a significant distinction between all four ACR categories. Mean T2* of the fibroglandular tissue correlated with increasing ACR categories (Spearman rho: 0.34, p < .001). The PDFF of the fibroglandular tissue showed a correlation with age (Pearson rho: 0.56, p = .03). CONCLUSION: The proposed breast PDFF as an automated tissue fat concentration measurement is comparable with mammographic breast density estimations. Therefore, it is a promising approach to an accurate, user-independent, and non-ionizing breast density assessment that could be easily incorporated into clinical routine breast MRI exams. KEY POINTS: • The proposed PDFF strongly negatively correlates with visually determined mammographic and MRI-based breast density estimations and therefore allows for an accurate, non-ionizing, and user-independent breast density measurement. • In combination with T2*, the PDFF can be used to track structural alterations in the composition of breast tissue for an individualized risk assessment for breast cancer.

Evaluation of a deep learning-based reconstruction method for denoising and image enhancement of shoulder MRI in patients with shoulder pain.

OBJECTIVES: To evaluate the diagnostic performance of an automated reconstruction algorithm combining MR imaging acquired using compressed SENSE (CS) with deep learning (DL) in order to reconstruct denoised high-quality images from undersampled MR images in patients with shoulder pain. METHODS: Prospectively, thirty-eight patients (14 women, mean age 40.0 ± 15.2 years) with shoulder pain underwent morphological MRI using a pseudo-random, density-weighted k-space scheme with an acceleration factor of 2.5 using CS only. An automated DL-based algorithm (CS DL) was used to create reconstructions of the same k-space data as used for CS reconstructions. Images were analyzed by two radiologists and assessed for pathologies, image quality, and visibility of anatomical landmarks using a 4-point Likert scale. RESULTS: Overall agreement for the detection of pathologies between the CS DL reconstructions and CS images was substantial to almost perfect (κ 0.95 (95% confidence interval 0.82-1.00)). Image quality and the visibility of the rotator cuff, articular cartilage, and axillary recess were overall rated significantly higher for CS DL images compared to CS (p < 0.03). Contrast-to-noise ratios were significantly higher for cartilage/fluid (CS DL 198 ± 24.3, CS 130 ± 32.2, p = 0.02) and ligament/fluid (CS DL 184 ± 17.3, CS 141 ± 23.5, p = 0.03) and SNR values were significantly higher for ligaments and muscle of the CS DL reconstructions (p < 0.04). CONCLUSION: Evaluation of shoulder pathologies was feasible using a DL-based algorithm for MRI reconstruction and denoising. In clinical routine, CS DL may be beneficial in particular for reducing image noise and may be useful for the detection and better discrimination of discrete pathologies. Assessment of shoulder pathologies was feasible with improved image quality as well as higher SNR using a compressed sensing deep learning-based framework for image reconstructions and denoising. KEY POINTS: • Automated deep learning-based reconstructions showed a significant increase in signal-to-noise ratio and contrast-to-noise ratio (p < 0.04) with only a slight increase of reconstruction time of 40 s compared to CS. • All pathologies were accurately detected with no loss of diagnostic information or prolongation of the scan time. • Significant improvements of the image quality as well as the visibility of the rotator cuff, articular cartilage, and axillary recess were detected.

Assessment of glenoid bone loss and other osseous shoulder pathologies comparing MR-based CT-like images with conventional CT.

OBJECTIVES: To evaluate and compare the diagnostic performance of CT-like images based on a 3D T1-weighted spoiled gradient-echo sequence (T1 GRE), an ultra-short echo time sequence (UTE), and a 3D T1-weighted spoiled multi-echo gradient-echo sequence (FRACTURE) with conventional CT in patients with suspected osseous shoulder pathologies. MATERIALS AND METHODS: Patients with suspected traumatic dislocation of the shoulder (n = 46, mean age 40 ± 14.5 years, 19 women) were prospectively recruited and received 3-T MR imaging including 3D T1 GRE, UTE, and 3D FRACTURE sequences. CT was performed in patients with acute fractures and served as standard of reference (n = 25). Agreement of morphological features between the modalities was analyzed including the glenoid bone loss, Hill-Sachs interval, glenoid track, and the anterior straight-line length. Agreement between the modalities was assessed using Bland-Altman plots, Student's t-test, and Pearson's correlation coefficient. Inter- and intrareader assessment was evaluated with weighted Cohen's κ and intraclass correlation coefficient. RESULTS: All osseous pathologies were detected accurately on all three CT-like sequences (n = 25, κ = 1.00). No significant difference in the percentage of glenoid bone loss was found between CT (mean ± standard deviation, 20.3% ± 8.0) and CT-like MR images (FRACTURE 20.6% ± 7.9, T1 GRE 20.4% ± 7.6, UTE 20.3% ± 7.7, p > 0.05). When comparing the different measurements on CT-like images, measurements performed using the UTE images correlated best with CT. CONCLUSION: Assessment of bony Bankart lesions and other osseous pathologies was feasible and accurate using CT-like images based on 3-T MRI compared with conventional CT. Compared to the T1 GRE and FRACTURE sequence, the UTE measurements correlated best with CT. CLINICAL RELEVANCE STATEMENT: In an acute trauma setting, CT-like images based on a T1 GRE, UTE, or FRACTURE sequence might be a useful alternative to conventional CT scan sparing associated costs as well as radiation exposure. KEY POINTS: • No significant differences were found for the assessment of the glenoid bone loss when comparing measurements of CT-like MR images with measurements of conventional CT images. • Compared to the T1 GRE and FRACTURE sequence, the UTE measurements correlated best with CT whereas the FRACTURE sequence appeared to be the most robust regarding motion artifacts. • The T1 GRE sequence had the highest resolution with high bone contrast and detailed delineation of even small fractures but was more susceptible to motion artifacts.

Diagnostic value of water-fat-separated images and CT-like susceptibility-weighted images extracted from a single ultrashort echo time sequence for the evaluation of vertebral fractures and degenerative changes of the spine.

OBJECTIVES: To evaluate the performance of single-echo Dixon water-fat imaging and computed tomography (CT)-like imaging based on a single ultrashort echo time (sUTE) MR sequence for imaging of vertebral fractures as well as degenerative bone changes of the spine in comparison to conventional CT and MR sequences. METHODS: Thirty patients with suspected acute vertebral fractures were examined using a 3-T MRI, including an sUTE sequence as well as short-tau inversion recovery (STIR) and T1-weighted sequences. During postprocessing, water-fat separation was performed by solving the smoothness-constrained inverse water-fat problem based on a single-complex UTE image. By removing the unwanted low-frequency phase terms, additional MR-based susceptibility-weighted-like (SW-like) images with CT-like contrast were created. Two radiologists evaluated semi-quantitative and quantitative features of fractures and degenerative changes independently and separately on CT and MR images. RESULTS: In total, all 58 fractures were accurately detected of whom 24 were correctly classified as acute fractures with an edema detected on the water-fat-separated UTE images, using STIR and T1w sequences as standard of reference. For the morphological assessment of fractures and degenerative changes, the overall agreement between SW-like images and CT was substantial to excellent (e.g., Genant: κ 0.90 (95% confidence interval 0.54-1.00); AO/Magerl: κ 0.75 (95% confidence interval 0.43-1.00)). Overall inter-reader agreement for water-fat-separated UTE images and SW-like images was substantial to almost perfect. CONCLUSION: Detection and assessment of vertebral fractures and degenerative bone changes of the spine were feasible and accurate using water-fat-separated images as well as SW-like images, both derived from the same sUTE-Dixon sequence. KEY POINTS: • The detection of acute vertebral fractures was feasible using water-fat-separated images and CT-like images reconstructed from one sUTE sequence. • Assessment of the vertebral fractures using SW-like images with CT-like contrast was found to be comparable to conventional CT. • sUTE imaging of the spine can help reduce examination times and radiation exposure.

Computed high-b-value high-resolution DWI improves solid lesion detection in IPMN of the pancreas.

OBJECTIVES: To examine the effect of high-b-value computed diffusion-weighted imaging (cDWI) on solid lesion detection and classification in pancreatic intraductal papillary mucinous neoplasm (IPMN), using endoscopic ultrasound (EUS) and histopathology as a standard of reference. METHODS: Eighty-two patients with known or suspected IPMN were retrospectively enrolled. Computed high-b-value images at b = 1000 s/mm2 were calculated from standard (b = 0, 50, 300, and 600 s/mm2) DWI images for conventional full field-of-view (fFOV, 3 × 3 × 4 mm3 voxel size) DWI. A subset of 39 patients received additional high-resolution reduced-field-of-view (rFOV, 2.5 × 2.5 × 3 mm3 voxel size) DWI. In this cohort, rFOV cDWI was compared against fFOV cDWI additionally. Two experienced radiologists evaluated (Likert scale 1-4) image quality (overall image quality, lesion detection and delineation, fluid suppression within the lesion). In addition, quantitative image parameters (apparent signal-to-noise ratio (aSNR), apparent contrast-to-noise ratio (aCNR), contrast ratio (CR)) were assessed. Diagnostic confidence regarding the presence/absence of diffusion-restricted solid nodules was assessed in an additional reader study. RESULTS: High-b-value cDWI at b = 1000 s/mm2 outperformed acquired DWI at b = 600 s/mm2 regarding lesion detection, fluid suppression, aCNR, CR, and lesion classification (p = < .001-.002). Comparing cDWI from fFOV and rFOV revealed higher image quality in high-resolution rFOV-DWI compared to conventional fFOV-DWI (p ≤ .001-.018). High-b-value cDWI images were rated non-inferior to directly acquired high-b-value DWI images (p = .095-.655). CONCLUSIONS: High-b-value cDWI may improve the detection and classification of solid lesions in IPMN. Combining high-resolution imaging and high-b-value cDWI may further increase diagnostic precision. CLINICAL RELEVANCE STATEMENT: This study shows the potential of computed high-resolution high-sensitivity diffusion-weighted magnetic resonance imaging for solid lesion detection in pancreatic intraductal papillary mucinous neoplasia (IPMN). The technique may enable early cancer detection in patients under surveillance. KEY POINTS: • Computed high-b-value diffusion-weighted imaging (cDWI) may improve the detection and classification of intraductal papillary mucinous neoplasms (IPMN) of the pancreas. • cDWI calculated from high-resolution imaging increases diagnostic precision compared to cDWI calculated from conventional-resolution imaging. • cDWI has the potential to strengthen the role of MRI for screening and surveillance of IPMN, particularly in view of the rising incidence of IPMNs combined with now more conservative therapeutic approaches.

Preconditioned water-fat total field inversion: Application to spine quantitative susceptibility mapping.

PURPOSE: To (a) develop a preconditioned water-fat total field inversion (wfTFI) algorithm that directly estimates the susceptibility map from complex multi-echo gradient echo data for water-fat regions and to (b) evaluate the performance of the proposed wfTFI quantitative susceptibility mapping (QSM) method in comparison with a local field inversion (LFI) method and a linear total field inversion (TFI) method in the spine. METHODS: Numerical simulations and in vivo spine multi-echo gradient echo measurements were performed to compare wfTFI to an algorithm based on disjoint background field removal (BFR) and LFI and to a formerly proposed TFI algorithm. The data from 1 healthy volunteer and 10 patients with metastatic bone disease were included in the analysis. Clinical routine computed tomography (CT) images were used as a reference standard to distinguish osteoblastic from osteolytic changes. The ability of the QSM methods to distinguish osteoblastic from osteolytic changes was evaluated. RESULTS: The proposed wfTFI method was able to decrease the normalized root mean square error compared to the LFI and TFI methods in the simulation. The in vivo wfTFI susceptibility maps showed reduced BFR artifacts, noise amplification, and streaking artifacts compared to the LFI and TFI maps. wfTFI provided a significantly higher diagnostic confidence in differentiating osteolytic and osteoblastic lesions in the spine compared to the LFI method (p = .012). CONCLUSION: The proposed wfTFI method can minimize BFR artifacts, noise amplification, and streaking artifacts in water-fat regions and can thus better differentiate between osteoblastic and osteolytic changes in patients with metastatic disease compared to LFI and the original TFI method.

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 B1 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.

Influence of aging and gadolinium exposure on T1, T2, and T2*-relaxation in healthy women with an increased risk of breast cancer with and without prior exposure to gadoterate meglumine at 3.0-T brain MR imaging.

OBJECTIVES: We examined the effects of aging and of gadolinium-based contrast agent (GBCA) exposure on MRI measurements in brain nuclei of healthy women. METHODS: This prospective, IRB-approved single-center case-control study enrolled 100 healthy participants of our high-risk screening center for hereditary breast cancer, who had received at least six doses of macrocyclic GBCA (exposed group) or were newly entering the program (GBCA-naïve group). The cutoff "at least six doses" was chosen to be able to include a sufficient number of highly exposed participants. All participants underwent unenhanced 3.0-T brain MRI including quantitative T1, T2, and R2* mapping and T1- and T2-weighted imaging. The relaxation times/signal intensities were derived from region of interest measurements in the brain nuclei performed by a radiologist and a neuroradiologist, both board certified. Statistical analysis was based on descriptive evaluations and uni-/multivariable analyses. RESULTS: The participants (exposed group: 49, control group: 51) were aged 42 ± 9 years. In a multivariable model, age had a clear impact on R2* (p < 0.001-0.012), T2 (p = 0.003-0.048), and T1 relaxation times/signal intensities (p < 0.004-0.046) for the majority of deep brain nuclei, mostly affecting the substantia nigra, globus pallidus (GP), nucleus ruber, thalamus, and dentate nucleus (DN). The effect of prior GBCA administration on T1 relaxation times was statistically significant for the DN, GP, and pons (p = 0.019-0.037). CONCLUSIONS: In a homogeneous group of young to middle-aged healthy females aging had an effect on T2 and R2* relaxation times and former GBCA applications influenced the measured T1 relaxation times. KEY POINTS: The quantitative T1, T2, and R2* relaxation times measured in women at high risk of developing breast cancer showed characteristic bandwidth for all brain nuclei examined at 3.0-T MRI. The effect of participant age had a comparatively strong impact on R2*, T2, and T1 relaxation times for the majority of brain nuclei examined. The effect of prior GBCA administrations on T1 relaxation times rates was comparatively less pronounced, yielding statistically significant results for the dentate nucleus, globus pallidus, and pons. Healthy women with and without previous GBCA-enhanced breast MRI exhibited age-related T2* and T2 relaxation alterations at 3.0 T-brain MRI. T1 relaxation alterations due to prior GBCA administration were comparatively less pronounced.

Deep learning-based acceleration of Compressed Sense MR imaging of the ankle.

OBJECTIVES: To evaluate a compressed sensing artificial intelligence framework (CSAI) to accelerate MRI acquisition of the ankle. METHODS: Thirty patients were scanned at 3T. Axial T2-w, coronal T1-w, and coronal/sagittal intermediate-w scans with fat saturation were acquired using compressed sensing only (12:44 min, CS), CSAI with an acceleration factor of 4.6-5.3 (6:45 min, CSAI2x), and CSAI with an acceleration factor of 6.9-7.7 (4:46 min, CSAI3x). Moreover, a high-resolution axial T2-w scan was obtained using CSAI with a similar scan duration compared to CS. Depiction and presence of abnormalities were graded. Signal-to-noise and contrast-to-noise were calculated. Wilcoxon signed-rank test and Cohen's kappa were used to compare CSAI with CS sequences. RESULTS: The correlation was perfect between CS and CSAI2x (κ = 1.0) and excellent for CS and CSAI3x (κ = 0.86-1.0). No significant differences were found for the depiction of structures between CS and CSAI2x and the same abnormalities were detected in both protocols. For CSAI3x the depiction was graded lower (p ≤ 0.001), though most abnormalities were also detected. For CSAI2x contrast-to-noise fluid/muscle was higher compared to CS (p ≤ 0.05), while no differences were found for other tissues. Signal-to-noise and contrast-to-noise were higher for CSAI3x compared to CS (p ≤ 0.05). The high - resolution axial T2-w sequence specifically improved the depiction of tendons and the tibial nerve (p ≤ 0.005). CONCLUSIONS: Acquisition times can be reduced by 47% using CSAI compared to CS without decreasing diagnostic image quality. Reducing acquisition times by 63% is feasible but should be reserved for specific patients. The depiction of specific structures is improved using a high-resolution axial T2-w CSAI scan. KEY POINTS: • Prospective study showed that CSAI enables reduction in acquisition times by 47% without decreasing diagnostic image quality. • Reducing acquisition times by 63% still produces images with an acceptable diagnostic accuracy but should be reserved for specific patients. • CSAI may be implemented to scan at a higher resolution compared to standard CS images without increasing acquisition times.

Gradient nonlinearity correction in liver DWI using motion-compensated diffusion encoding waveforms.

OBJECTIVE : To experimentally characterize the effectiveness of a gradient nonlinearity correction method in removing ADC bias for different motion-compensated diffusion encoding waveforms. METHODS: The diffusion encoding waveforms used were the standard monopolar Stejskal-Tanner pulsed gradient spin echo (pgse) waveform, the symmetric bipolar velocity-compensated waveform (sym-vc), the asymmetric bipolar velocity-compensated waveform (asym-vc) and the asymmetric bipolar partial velocity-compensated waveform (asym-pvc). The effectiveness of the gradient nonlinearity correction method using the spherical harmonic expansion of the gradient coil field was tested with the aforementioned waveforms in a phantom and in four healthy subjects. RESULTS: The gradient nonlinearity correction method reduced the ADC bias in the phantom experiments for all used waveforms. The range of the ADC values over a distance of ± 67.2 mm from isocenter reduced from 1.29 × 10-4 to 0.32 × 10-4 mm2/s for pgse, 1.04 × 10-4 to 0.22 × 10-4 mm2/s for sym-vc, 1.22 × 10-4 to 0.24 × 10-4 mm2/s for asym-vc and 1.07 × 10-4 to 0.11 × 10-4 mm2/s for asym-pvc. The in vivo results showed that ADC overestimation due to motion or bright vessels can be increased even further by the gradient nonlinearity correction. CONCLUSION: The investigated gradient nonlinearity correction method can be used effectively with various motion-compensated diffusion encoding waveforms. In coronal liver DWI, ADC errors caused by motion and residual vessel signal can be increased even further by the gradient nonlinearity correction.

Cerebrospinal fluid pulse wave velocity measurements: In vitro and in vivo evaluation of a novel multiband cine phase-contrast MRI sequence.

PURPOSE: Intracranial and intraspinal compliance are parameters of interest in the diagnosis and prediction of treatment outcome in patients with normal pressure hydrocephalus and other forms of communicating hydrocephalus. A noninvasive method to estimate the spinal cerebrospinal fluid (CSF) pulse wave velocity (PWV) as a measure of compliance was developed using a multiband cine phase-contrast MRI sequence and a foot-to-foot algorithm. METHODS: We used computational simulations to estimate the accuracy of the MRI acquisition and transit-time algorithm. In vitro measurements were performed to investigate the reproducibility and accuracy of the measurements under controlled conditions. In vivo measurements in 20 healthy subjects and 2 patients with normal pressure hydrocephalus were acquired to show the technical feasibility in a clinical setting. RESULTS: Simulations showed a mean deviation of the calculated CSF PWV of 3.41% ± 2.68%. In vitro results were in line with theory, showing a square-root relation between PWV and transmural pressure and a good reproducibility with SDs of repeated measurements below 5%. Mean CSF PWV over all healthy subjects was 5.83 ± 3.36 m/s. The CSF PWV measurements in the patients with normal pressure hydrocephalus were distinctly higher before CSF shunt surgery (33.80 ± 6.75 m/s and 31.31 ± 7.82 m/s), with a decrease 5 days after CSF shunt surgery (15.69 ± 3.37 m/s). CONCLUSION: This study evaluates the feasibility of CSF PWV measurements using a multiband cine phase-contrast MRI sequence. In vitro and in vivo measurements showed that this method is a potential tool for the noninvasive estimation of intraspinal compliance.

Respiratory motion corrected 4D flow using golden radial phase encoding.

PURPOSE: To minimize respiratory motion artifacts while achieving predictable scan times with 100% scan efficiency for thoracic 4D flow MRI. METHODS: A 4D flow sequence with golden radial phase encoding (GRPE) was acquired in 9 healthy volunteers covering the heart, aorta, and venae cavae. Scan time was 15 min, and data were acquired without motion gating during acquisition. Data were retrospectively re-binned into respiratory and cardiac phases based on respiratory self-navigation and the electrocardiograph signals, respectively. Nonrigid respiratory motion fields were extracted and corrected for during the k-t SENSE reconstruction. A respiratory-motion corrected (GRPE-MOCO) and a free-breathing (GRPE-UNCORR) 4D flow dataset was reconstructed using 100% of the acquired data. For comparison, a respiratory gated Cartesian 4D flow acquisition (CART-REF) covering the aorta was acquired. Stroke volumes and peak flows were compared. Additionally, an internal flow validation based on mass conservation was performed on the GRPE-MOCO and GRPE-UNCORR. Statistically significant differences were analyzed using a paired Wilcoxon test. RESULTS: Stroke volumes and peak flows in the aorta between GRPE-MOCO and the CART-REF showed a mean difference of -1.5 ± 10.3 mL (P > 0.05) and 25.2 ± 55.9 mL/s (P > 0.05), respectively. Peak flow in the GRPE-UNCORR data was significantly different compared with CART-REF (P < 0.05). GRPE-MOCO showed higher accuracy for internal consistency analysis than GRPE-UNCORR. CONCLUSION: The proposed 4D flow sequence allows a straight-forward planning by covering the entire thorax and ensures a predictable scan time independent of cardiac cycle variations and breathing patterns.