Imaging of osteoarthritis from the ankle through the midfoot

Ankle, hindfoot, and midfoot osteoarthritis (OA) is most commonly posttraumatic and tends to become symptomatic in younger patients. It often results from instability due to insufficiency of supportive soft tissue structures, such as ligaments and tendons. Diagnostic imaging can be helpful to detect and characterize the distribution of OA, and to assess the integrity of these supportive structures, which helps determine prognosis and guide treatment. However, the imaging findings associated with OA and instability may be subtle and unrecognized until the process is advanced, which may ultimately limit therapeutic options to salvage procedures. It is important to understand the abilities and limitations of various imaging modalities used to assess ankle, hindfoot, and midfoot OA, and to be familiar with the imaging findings of OA and instability patterns.

Cancellation of streak artifacts in radial abdominal imaging using interference null space projection

PURPOSE

In radial abdominal imaging, it has been commonly observed that signal from the arms cause streaks due to system imperfections. We previously introduced a streak removal technique (B-STAR), which is inherently spatially variant and limited to work in image space. In this work, we propose a spatially invariant streak cancellation technique (CACTUS), which can be applied in either image space or k-space and is compatible with iterative reconstructions.

THEORY AND METHODS

Streak sources are typically spatially localized and can be represented using a low-dimensional subspace. CACTUS identifies the streak subspace by leveraging the spatial redundancy of receiver coils and projects the data onto the streak null space to eliminate the streaks. When applied in k-space, CACTUS can be combined with iterative reconstructions. CACTUS was tested in phantoms and in vivo abdominal imaging using a radial turbo spin-echo pulse sequence.

RESULTS

In phantoms, CACTUS improved T2 estimation in comparison to previous de-streaking methods. In vivo experiments showed that CACTUS reduced streaks and yielded T2 estimation, in regions affected by streaks, closer to a streak-free reference. Evaluation using a clinical abdominal dataset (n = 20) showed that CACTUS is comparable to B-STAR and yields significantly better signal preservation and streak cancellation than coil removal and suppression methods.

CONCLUSION

CACTUS provides superior signal preservation and streak reduction performance compared to coil removal and suppression methods. As a clear advantage over B-STAR, CACTUS can be integrated with iterative reconstruction methods. In abdominal T2 mapping, CACTUS improves the accuracy of parameter estimation in areas affected by streaks.

Abdominal T2-Weighted Imaging and T2 Mapping Using a Variable Flip Angle Radial Turbo Spin-Echo Technique

BACKGROUND

T2 mapping is of great interest in abdominal imaging but current methods are limited by low resolution, slice coverage, motion sensitivity, or lengthy acquisitions.

PURPOSE

Develop a radial turbo spin-echo technique with refocusing variable flip angles (RADTSE-VFA) for high spatiotemporal T2 mapping and efficient slice coverage within a breath-hold and compare to the constant flip angle counterpart (RADTSE-CFA).

STUDY TYPE

Prospective technical efficacy.

SUBJECTS

Testing performed on agarose phantoms and 12 patients. Focal liver lesion classification tested on malignant (N = 24) and benign (N = 11) lesions.

FIELD STRENGTH/SEQUENCE

1.5 T/RADTSE-VFA, RADTSE-CFA.

ASSESSMENT

A constrained objective function was used to optimize the refocusing flip angles. Phantom and/or in vivo data were used to assess relative contrast, T2 estimation, specific absorption rate (SAR), and focal liver lesion classification. STATISTICAL TESTS: t-Tests or Mann-Whitney Rank Sum tests were used.

RESULTS

Phantom data did not show significant differences in mean relative contrast (P = 0.10) and T2 accuracy (P = 0.99) between RADTSE-VFA and RADTSE-CFA. Adding noise caused T2 overestimation predominantly for RADTSE-CFA and low T2 values. In vivo results did not show significant differences in mean spleen-to-liver (P = 0.62) and kidney-to-liver (P = 0.49) relative contrast between RADTSE-VFA and RADTSE-CFA. Mean T2 values were not significantly different between the two techniques for spleen (T2VFA  = 109.2 ± 12.3 msec; T2CFA  = 110.7 ± 11.1 msec; P = 0.78) and kidney-medulla (T2VFA  = 113.0 ± 8.7 msec; T2CFA  = 114.0 ± 8.6 msec; P = 0.79). Liver T2 was significantly higher for RADTSE-CFA (T2VFA  = 52.6 ± 6.6 msec; T2CFA  = 60.4 ± 8.0 msec) consistent with T2 overestimation in the phantom study. Focal liver lesion classification had comparable T2 distributions for RADTSE-VFA and RADTSE-CFA for malignancies (P = 1.0) and benign lesions (P = 0.39). RADTSE-VFA had significantly lower SAR than RADTSE-CFA increasing slice coverage by 1.5.

DATA CONCLUSION

RADTSE-VFA provided noise-robust T2 estimation compared to the constant flip angle counterpart while generating T2-weighted images with comparable contrast. The VFA scheme minimized SAR improving slice efficiency for breath-hold imaging.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 1.

Rapid high-resolution volumetric T1 mapping using a highly accelerated stack-of-stars Look Locker technique

PURPOSE

To develop a fast volumetric T1 mapping technique.

MATERIALS AND METHODS

A stack-of-stars (SOS) Look Locker technique based on the acquisition of undersampled radial data (>30× relative to Nyquist) and an efficient multi-slab excitation scheme is presented. A principal-component based reconstruction is used to reconstruct T1 maps. Computer simulations were performed to determine the best choice of partitions per slab and degree of undersampling. The technique was validated in phantoms against reference T1 values measured with a 2D Cartesian inversion-recovery spin-echo technique. The SOS Look Locker technique was tested in brain (n = 4) and prostate (n = 5). Brain T1 mapping was carried out with and without kz acceleration and results between the two approaches were compared. Prostate T1 mapping was compared to standard techniques. A reproducibility study was conducted in brain and prostate. Statistical analyses were performed using linear regression and Bland Altman analysis.

RESULTS

Phantom T1 values showed excellent correlations between SOS Look Locker and the inversion-recovery spin-echo reference (r2 = 0.9965; p < 0.0001) and between SOS Look Locker with slab-selective and non-slab selective inversion pulses (r2 = 0.9999; p < 0.0001). In vivo results showed that full brain T1 mapping (1 mm3) with kz acceleration is achieved in 4 min 21 s. Full prostate T1 mapping (0.9 × 0.9 × 4 mm3) is achieved in 2 min 43 s. T1 values for brain and prostate were in agreement with literature values. A reproducibility study showed coefficients of variation in the range of 0.18-0.2% (brain) and 0.15-0.18% (prostate).

CONCLUSION

A rapid volumetric T1 mapping technique was developed. The technique enables high-resolution T1 mapping with adequate anatomical coverage in a clinically acceptable time.

A Stacked Generalization of 3D Orthogonal Deep Learning Convolutional Neural Networks for Improved Detection of White Matter Hyperintensities in 3D FLAIR Images

BACKGROUND AND PURPOSE

Accurate and reliable detection of white matter hyperintensities and their volume quantification can provide valuable clinical information to assess neurologic disease progression. In this work, a stacked generalization ensemble of orthogonal 3D convolutional neural networks, StackGen-Net, is explored for improving automated detection of white matter hyperintensities in 3D T2-FLAIR images.

MATERIALS AND METHODS

Individual convolutional neural networks in StackGen-Net were trained on 2.5D patches from orthogonal reformatting of 3D-FLAIR (n = 21) to yield white matter hyperintensity posteriors. A meta convolutional neural network was trained to learn the functional mapping from orthogonal white matter hyperintensity posteriors to the final white matter hyperintensity prediction. The impact of training data and architecture choices on white matter hyperintensity segmentation performance was systematically evaluated on a test cohort (n = 9). The segmentation performance of StackGen-Net was compared with state-of-the-art convolutional neural network techniques on an independent test cohort from the Alzheimer's Disease Neuroimaging Initiative-3 (n = 20).

RESULTS

StackGen-Net outperformed individual convolutional neural networks in the ensemble and their combination using averaging or majority voting. In a comparison with state-of-the-art white matter hyperintensity segmentation techniques, StackGen-Net achieved a significantly higher Dice score (0.76 [SD, 0.08], F1-lesion (0.74 [SD, 0.13]), and area under precision-recall curve (0.84 [SD, 0.09]), and the lowest absolute volume difference (13.3% [SD, 9.1%]). StackGen-Net performance in Dice scores (median = 0.74) did not significantly differ (P = .22) from interobserver (median = 0.73) variability between 2 experienced neuroradiologists. We found no significant difference (P = .15) in white matter hyperintensity lesion volumes from StackGen-Net predictions and ground truth annotations.

CONCLUSIONS

A stacked generalization of convolutional neural networks, utilizing multiplanar lesion information using 2.5D spatial context, greatly improved the segmentation performance of StackGen-Net compared with traditional ensemble techniques and some state-of-the-art deep learning models for 3D-FLAIR.

Improving subspace constrained radial fast spin echo MRI using block matching driven non-local low rank regularization

Subspace-constrained reconstruction methods restrict the relaxation signals (of size M) in the scene to a pre-determined subspace (of size K≪M) and allow multi-contrast imaging and parameter mapping from accelerated acquisitions. However, these constraints yield poor image quality at some imaging contrasts, which can impact the parameter mapping performance. Additional regularization such as the use of joint-sparse (JS) or locally-low-rank (LLR) constraints can help improve the recovery of these images but are not sufficient when operating at high acceleration rates. We propose a method, non-local rank 3D (NLR3D), that is built on block matching and transform domain low rank constraints to allow high quality recovery of subspace-coefficient images (SCI) and subsequent multi-contrast imaging and parameter mapping. The performance of NLR3D was evaluated using Monte-Carlo (MC) simulations and compared against the JS and LLR methods. In vivo T 2 mapping results are presented on brain and knee datasets. MC results demonstrate improved bias, variance, and MSE behavior in both the multi-contrast images and parameter maps when compared to the JS and LLR methods. In vivo brain and knee results at moderate and high acceleration rates demonstrate improved recovery of high SNR early TE images as well as parameter maps. No significant difference was found in the T2 values measured in ROIs between the NLR3D reconstructions and the reference images (Wilcoxon signed rank test). The proposed method, NLR3D, enables recovery of high-quality SCI and, consequently, the associated multi-contrast images and parameter maps.

A multi-scale residual network for accelerated radial MR parameter mapping

A deep learning MR parameter mapping framework which combines accelerated radial data acquisition with a multi-scale residual network (MS-ResNet) for image reconstruction is proposed. The proposed supervised learning strategy uses input image patches from multi-contrast images with radial undersampling artifacts and target image patches from artifact-free multi-contrast images. Subspace filtering is used during pre-processing to denoise input patches. For each anatomy and relaxation parameter, an individual network is trained. in vivo T1 mapping results are obtained on brain and abdomen datasets and in vivo T2 mapping results are obtained on brain and knee datasets. Quantitative results for the T2 mapping of the knee show that MS-ResNet trained using either fully sampled or undersampled data outperforms conventional model-based compressed sensing methods. This is significant because obtaining fully sampled training data is not possible in many applications. in vivo brain and abdomen results for T1 mapping and in vivo brain results for T2 mapping demonstrate that MS-ResNet yields contrast-weighted images and parameter maps that are comparable to those achieved by model-based iterative methods while offering two orders of magnitude reduction in reconstruction times. The proposed approach enables recovery of high-quality contrast-weighted images and parameter maps from highly accelerated radial data acquisitions. The rapid image reconstructions enabled by the proposed approach makes it a good candidate for routine clinical use.

Fully Automated Segmentation of Globes for Volume Quantification in CT Images of Orbits using Deep Learning

BACKGROUND AND PURPOSE

Fast and accurate quantification of globe volumes in the event of an ocular trauma can provide clinicians with valuable diagnostic information. In this work, an automated workflow using a deep learning-based convolutional neural network is proposed for prediction of globe contours and their subsequent volume quantification in CT images of the orbits.

MATERIALS AND METHODS

An automated workflow using a deep learning -based convolutional neural network is proposed for prediction of globe contours in CT images of the orbits. The network, 2D Modified Residual UNET (MRes-UNET2D), was trained on axial CT images from 80 subjects with no imaging or clinical findings of globe injuries. The predicted globe contours and volume estimates were compared with manual annotations by experienced observers on 2 different test cohorts.

RESULTS

On the first test cohort (n = 18), the average Dice, precision, and recall scores were 0.95, 96%, and 95%, respectively. The average 95% Hausdorff distance was only 1.5 mm, with a 5.3% error in globe volume estimates. No statistically significant differences (P = .72) were observed in the median globe volume estimates from our model and the ground truth. On the second test cohort (n = 9) in which a neuroradiologist and 2 residents independently marked the globe contours, MRes-UNET2D (Dice = 0.95) approached human interobserver variability (Dice = 0.94). We also demonstrated the utility of inter-globe volume difference as a quantitative marker for trauma in 3 subjects with known globe injuries.

CONCLUSIONS

We showed that with fast prediction times, we can reliably detect and quantify globe volumes in CT images of the orbits across a variety of acquisition parameters.

An efficient 3D stack-of-stars turbo spin echo pulse sequence for simultaneous T2-weighted imaging and T2 mapping

PURPOSE

To design a pulse sequence for efficient 3D T2-weighted imaging and T2 mapping.

METHODS

A stack-of-stars turbo spin echo pulse sequence with variable refocusing flip angles and a flexible pseudorandom view ordering is proposed for simultaneous T2-weighted imaging and T2 mapping. An analytical framework is introduced for the selection of refocusing flip angles to maximize relative tissue contrast while minimizing T2 estimation errors and maintaining low specific absorption rate. Images at different echo times are generated using a subspace constrained iterative reconstruction algorithm. T2 maps are obtained by modeling the signal evolution using the extended phase graph model. The technique is evaluated using phantoms and demonstrated in vivo for brain, knee, and carotid imaging.

RESULTS

Numerical simulations demonstrate an improved point spread function with the proposed pseudorandom view ordering compared to golden angle view ordering. Phantom experiments show that T2 values estimated from the stack-of-stars turbo spin echo pulse sequence with variable refocusing flip angles have good concordance with spin echo reference values. In vivo results show the proposed pulse sequence can generate qualitatively comparable T2-weighted images as conventional Cartesian 3D SPACE in addition to simultaneously generating 3D T2 maps.

CONCLUSION

The proposed stack-of-stars turbo spin echo pulse sequence with pseudorandom view ordering and variable refocusing flip angles allows high resolution isotropic T2 mapping in clinically acceptable scan times. The optimization framework for the selection of refocusing flip angles improves T2 estimation accuracy while generating T2-weighted contrast comparable to conventional Cartesian imaging.

Radial streak artifact reduction using phased array beamforming

PURPOSE

A new method for streak artifact reduction in radial MRI based on phased array filtering.

THEORY

Radial imaging in applications that require large fields-of-view can be susceptible to streaking artifacts due to gradient nonlinearities. Coil removal methods prune the coils contributing the most to streaking artifacts at the expense of signal loss. Phased array beamforming is a form of spatial filtering used to suppress unwanted signals. The proposed method uses interference covariance generated from the streaking artifact samples which are manually extracted with phased array beamforming to suppress streaking in the images.

METHODS

The performance of the proposed method was evaluated on abdomen radial fast spin echo images acquired on a 1.5T Siemens scanner and compared with previously proposed methods.

RESULTS

Our results demonstrate that the proposed method can effectively suppress streaking artifacts without any noticeable loss in signal levels. Coil removal methods can suppress streaks as well but they may incur significant signal loss due to coil pruning. Quantitative metrics also demonstrate the superiority of the proposed method over earlier methods.

CONCLUSION

The use of interference covariance with phased array beamforming can help reduce streaking artifacts.

Validation of Peripheral Quantitative Computed Tomography-Derived Thigh Adipose Tissue Subcompartments in Young Girls Using a 3 T MRI Scanner

The ability to assess skeletal muscle adipose tissue is important given the negative clinical implications associated with greater fat infiltration of the muscle. Computed tomography and magnetic resonance imaging (MRI) are highly accurate for measuring appendicular soft tissue and muscle composition, but have limitations. Peripheral quantitative computed tomography (pQCT) is an alternative that investigators find valuable because of its low radiation, fast scan time, and comparatively lower costs. The present investigation sought to assess the accuracy of pQCT-derived estimates of total, subcutaneous, skeletal muscle, intermuscular, and calculated intramuscular adipose tissue areas, and muscle density in the midthigh of young girls using the gold standard, 3 T MRI, as the criterion. Cross-sectional data were analyzed for 26 healthy girls aged 9-12 years. Midthigh soft tissue composition was assessed by both pQCT and 3 T MRI. Mean tissue area for corresponding adipose compartments by pQCT and MRI was compared using t tests, regression analysis, and Bland-Altman plots. Muscle density was regressed on MRI skeletal muscle adipose tissue, intermuscular adipose tissue, and intramuscular adipose tissue, each expressed as a percentage of total muscle area. Correlations were high between MRI and pQCT for total adipose tissue (r² = 0.98), subcutaneous adipose tissue (r² = 0.95), skeletal muscle adipose tissue (r² = 0.83), and intermuscular adipose tissue (r² = 0.82), and pQCT muscle density correlated well with both MRI skeletal muscle adipose tissue (r² = 0.70) and MRI intermuscular adipose tissue (r² = 0.70). There was a slight, but statistically significant underestimation by pQCT for total and subcutaneous adipose tissue, whereas no significant difference was observed for skeletal muscle adipose tissue. Both pQCT-estimated intramuscular adipose tissue and muscle density were weakly correlated with MRI-intramuscular adipose tissue. We conclude that pQCT is a valid measurement technique for estimating all adipose subcompartments, except for intramuscular adipose tissue, for the midthigh region in young/adolescent girls.

Rapid high-resolution T1 mapping using a highly accelerated radial steady-state free-precession technique

BACKGROUND

T₁ mapping is often used in some clinical protocols. Existing techniques are limited in slice coverage, and/or spatial-temporal resolution, or require long acquisitions. Here we present a multi-slice inversion-recovery (IR) radial steady-state free precession (radSSFP) pulse sequence combined with a principal component (PC) based reconstruction that overcomes these limitations.

PURPOSE

To develop a fast technique for multi-slice high-resolution T₁ mapping.

STUDY TYPE

Technical efficacy study done prospectively.

PHANTOM/SUBJECTS

IR-radSSFP was tested in phantoms, five healthy volunteers, and four patients with abdominal lesions.

FIELD STRENGTH/SEQUENCE

IR-radSSFP was implemented at 3T.

ASSESSMENT

Computer simulations were performed to optimize the flip angle for T₁ estimation; testing was done in phantoms using as reference an IR spin-echo pulse sequence. T₁ mapping with IR-radSSFP was also assessed in vivo (brain and abdomen) and T₁ values were compared with literature. T₁ maps were also compared with a radial IR-FLASH technique.

STATISTICAL TESTS

A two-tailed t-test was used to compare T₁ values in phantoms. A repeatability study was carried out in vivo using Bland-Altman analysis.

RESULTS

Simulations and phantom experiments showed that a flip angle of 20˚ was optimal for T₁ mapping. When comparing single to multi-slice experiments in phantoms there were no significant differences between the means T₁ values (P = 0.0475). In vivo results show that T₁ maps with spatial resolution as high as 0.69 mm × 0.69 mm × 2.00 mm (brain) and 0.83 mm × 0.83 mm × 3.00 mm (abdomen) can be generated for 84 brain slices in 3 min and 10 abdominal slices in a breath-hold; T₁ values were comparable to those reported in literature. The coefficients of variation from the repeatability study were 1.7% for brain and 2.5-2.7% in the abdomen.

DATA CONCLUSION

A multi-slice IR-radSSFP technique combined with a PC-based reconstruction was demonstrated for higher resolution T₁ mapping. This technique is fast, motion-insensitive and yields repeatable T₁ values comparable to those in literature.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;49:239-252.

A multi-band double-inversion radial fast spin-echo technique for T2 cardiovascular magnetic resonance mapping of the heart

BACKGROUND

Double inversion recovery (DIR) fast spin-echo (FSE) cardiovascular magnetic resonance (CMR) sequences are used clinically for black-blood T2-weighted imaging. However, these sequences suffer from slice inefficiency due to the non-selective inversion pulses. We propose a multi-band (MB) encoded DIR radial FSE (MB-DIR-RADFSE) technique to simultaneously excite two slices. This sequence has improved signal-to-noise ratio per unit time compared to a single slice excitation. It is also motion robust and enables the reconstruction of high-resolution black-blood T2-weighted images and T2 maps for the excited slices.

METHODS

Hadamard encoded MB pulses were used in MB-DIR-RADFSE to simultaneously excite two slices. A principal component based iterative reconstruction was used to jointly reconstruct black-blood T2-weighted images and T2 maps. Phantom and in vivo experiments were performed to evaluate T2 mapping performance and results were compared to a T2-prepared balanced steady state free precession (bSSFP) method. The inter-segment variability of the T2 maps were assessed using data acquired on healthy subjects. A reproducibility study was performed to evaluate reproducibility of the proposed technique.

RESULTS

Phantom experiments show that the T2 values estimated from MB-DIR-RADFSE are comparable to the spin-echo based reference, while T2-prepared bSSFP over-estimated T2 values. The relative contrast of the black-blood images from the multi-band scheme was comparable to those from a single slice acquisition. The myocardial segment analysis on 8 healthy subjects indicated a significant difference (p-value < 0.01) in the T2 estimates from the apical slice when compared to the mid-ventricular slice. The mean T2 estimate from 12 subjects obtained using T2-prepared bSSFP was significantly higher (p-value = 0.012) compared to MB-DIR-RADFSE, consistent with the phantom results. The Bland-Altman analysis showed excellent reproducibility between the MB-DIR-RADFSE measurements, with a mean T2 difference of 0.12 ms and coefficient of reproducibility of 2.07 in 15 clinical subjects. The utility of this technique is demonstrated in two subjects where the T2 maps show elevated values in regions of pathology.

CONCLUSIONS

The use of multi-band pulses for excitation improves the slice efficiency of the double inversion fast spin-echo pulse sequence. The use of a radial trajectory and a joint reconstruction framework allows reconstruction of TE images and T2 maps for the excited slices.

Accelerated MR parameter mapping with a union of local subspaces constraint

PURPOSE

A new reconstruction method for multi-contrast imaging and parameter mapping based on a union of local subspaces constraint is presented.

THEORY

Subspace constrained reconstructions use a predetermined subspace to explicitly constrain the relaxation signals. The choice of subspace size ( K ) impacts the approximation error vs noise-amplification tradeoff associated with these methods. A different approach is used in the model consistency constraint (MOCCO) framework to leverage the subspace model to enforce a softer penalty. Our proposed method, MOCCO-LS, augments the MOCCO model with a union of local subspaces (LS) approach. The union of local subspaces model is coupled with spatial support constraints and incorporated into the MOCCO framework to regularize the contrast signals in the scene.

METHODS

The performance of the MOCCO-LS method was evaluated in vivo on T₁ and T₂ mapping of the human brain and with Monte-Carlo simulations and compared against MOCCO and the explicit subspace constrained models.

RESULTS

The results demonstrate a clear improvement in the multi-contrast images and parameter maps. We sweep across the model order space ( K ) to compare the different reconstructions and demonstrate that the reconstructions have different preferential operating points. Experiments on T₂ mapping show that the proposed method yields substantial improvements in performance even when operating at very high acceleration rates.

CONCLUSIONS

The use of a union of local subspace constraints coupled with a sparsity promoting penalty leads to improved reconstruction quality of multi-contrast images and parameter maps.

Clinical Utility of a Novel Ultrafast T2-Weighted Sequence for Spine Imaging

BACKGROUND AND PURPOSE

TSE-based T2-weighted imaging of the spine has long scan times. This work proposes a fast imaging protocol using variable refocusing flip angles, optimized for blurring and specific absorption rate.

MATERIALS AND METHODS

A variable refocusing flip angle echo-train was optimized for the spine to improve the point spread function and minimize the specific absorption rate, yielding images with improved spatial resolution and SNR compared with the constant flip angle sequence. Data were acquired from 51 patients (35 lumbar, 16 whole-spine) using conventional TSE and the proposed sequence, with a single-shot variant for whole-spine. Noninferiority analysis was performed to evaluate the efficiency of the proposed technique.

RESULTS

The proposed multishot sequence resulted in a 2× shorter scan time with a >1.5× lower specific absorption rate. The variable flip angle sequence was noninferior to the conventional TSE (P < .025) for all image-quality and clinical criteria except signal-to-noise ratio for the lumbar spine protocol. However, mean image scores for the TSE-variable refocusing flip angle were ≥4.3 for all criteria, and concordance analysis showed high agreement (>90%) with the TSE, indicating clinical equivalence. The single-shot sequence resulted in 4× shorter whole-spine scans, and image scores were ≥4.4 for all criteria, attesting to its clinical utility.

CONCLUSIONS

We present a fast T2-weighted spine protocol using variable refocusing flip angles, including a single-shot variant. The sequences have better point spread function behavior than their constant flip angle counterparts and, being faster, should be less sensitive to patient motion, often seen in the longer TSE scans.

Reproducible automated breast density measure with no ionizing radiation using fat-water decomposition MRI

BACKGROUND

Increased breast density is a significant independent risk factor for breast cancer, and recent studies show that this risk is modifiable. Hence, breast density measures sensitive to small changes are desired.

PURPOSE

Utilizing fat-water decomposition MRI, we propose an automated, reproducible breast density measurement, which is nonionizing and directly comparable to mammographic density (MD).

STUDY TYPE

Retrospective study.

POPULATION

The study included two sample sets of breast cancer patients enrolled in a clinical trial, for concordance analysis with MD (40 patients) and reproducibility analysis (10 patients).

FIELD STRENGTH/SEQUENCE

The majority of MRI scans (59 scans) were performed with a 1.5T GE Signa scanner using radial IDEAL-GRASE sequence, while the remaining (seven scans) were performed with a 3T Siemens Skyra using 3D Cartesian 6-echo GRE sequence with a similar fat-water separation technique.

ASSESSMENT

After automated breast segmentation, breast density was calculated using FraGW, a new measure developed to reliably reflect the amount of fibroglandular tissue and total water content in the entire breast. Based on its concordance with MD, FraGW was calibrated to MR-based breast density (MRD) to be comparable to MD. A previous breast density measurement, Fra80-the ratio of breast voxels with <80% fat fraction-was also calculated for comparison with FraGW.

STATISTICAL TESTS

Pearson correlation was performed between MD (reference standard) and FraGW (and Fra80). Test-retest reproducibility of MRD was evaluated using the difference between test-retest measures (Δ_1-2 ) and intraclass correlation coefficient (ICC).

RESULTS

Both FraGW and Fra80 were strongly correlated with MD (Pearson ρ: 0.96 vs. 0.90, both P < 0.0001). MRD converted from FraGW showed higher test-retest reproducibility (Δ_1-2 variation: 1.1% ± 1.2%; ICC: 0.99) compared to MD itself (literature intrareader ICC ≤0.96) and Fra80.

DATA CONCLUSION

The proposed MRD is directly comparable with MD and highly reproducible, which enables the early detection of small breast density changes and treatment response.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018;48:971-981.

Segmentation of the right ventricle in four chamber cine cardiac MR images using polar dynamic programming

The four chamber plane is currently underutilized in the right ventricular segmentation community. Four chamber information can be useful to determine ventricular short axis stacks and provide a rough estimate of the right ventricle in short axis stacks. In this study, we develop and test a semi-automated technique for segmenting the right ventricle in four chamber cine cardiac magnetic resonance images. The three techniques that use minimum cost path algorithms were used. The algorithms are: Dijkstra's shortest path algorithm (Dijkstra), an A* algorithm that uses length, curvature and torsion into an active contour model (ALCT), and a variation of polar dynamic programming (PDP). The techniques are evaluated against the expert traces using 175 cardiac images from 7 patients. The evaluation first looks at mutual overlap metrics and then focuses on clinical measures such as fractional area change (FAC). The mean mutual overlap between the physician's traces ranged from 0.85 to 0.88. Using as reference physician 1's landmarks and traces (i.e., comparing the traces from physician 1 to the semi-automated segmentation using physician 1's landmarks), the PDP algorithm has a mean mutual overlap of 0.8970 compared to 0.8912 for ALCT and 0.8879 for Dijkstra. The mean mutual overlap between the BP regions generated by physician 1 and physician 2 landmarks are 0.9674, 0.9605 and 0.9531 for PDP, ALCT and Dijkstra, respectively. The FAC correlation coefficient between the physician's traces ranged from 0.73 to 0.93.

Abstract P6-09-19: Breast density change at 6 months is associated with change at 12 months as measured by fat-water decomposition MRI in women on tamoxifen

Objective: Tamoxifen (TAM) lowers breast cancer recurrence by 40-50% with evidence of individual variability in responsiveness. A ≥10% decrease in mammography-determined breast density (BD) after 12–18 months of TAM use has been associated with clinical benefit. Early determination of changes in BD may offer a strategy to tailor hormone therapy in non-responders; for responders, it may encourage adherence. Fat-water decomposition MRI (FWD-MRI) is an accurate and fast (&lt; 5 minutes) method for measuring BD without ionizing radiation or contrast agent. Here, we examined whether change in FWD-MRI-derived BD predicts decrease in BD at earlier time points than observable with a 12-month measure of BD.

Methods: The study population included a subset of 44 pre- and post-menopausal women receiving TAM for treatment of early-stage breast cancer or prevention who were enrolled in a randomized, placebo-controlled trial of diindolylmethane. Eligibility for this analysis included participants with FWD-MRI scans at baseline, 6 and 12 months. Median time on TAM at baseline was 13 months (IQR, 5–26 months). All MRI images were acquired on a 1.5T GE Signa NV-CV/i scanner. Automated breast segmentation was performed using MATLAB software and validated against manual ROI drawings. MRI-based BD was calculated as the ratio of breast voxels with &lt;80% apparent fat fraction (Fra80) over the entire breast, a measure previously shown by our group to be highly correlated with mammography-derived BD. For 40 participants, the unaffected, contralateral breast was analyzed. For 4 patients with two unaffected breasts, BD data from the left breast were analyzed. Change in BD was conservatively defined as &gt; 2 times the test-retest variability of Fra80 (0.032). McNemar's test was used to test the association between change from baseline to 6 months and change from baseline to 12 months.

Results and Discussion: At 12 months, 15 (34%) participants had a decrease in BD, whereas 29 (66%) remained unchanged or increased. Of these 29, 28 also had no decrease at 6 months (specificity = 97%), and 9 of the 15 women who showed a decrease at 12 months had a decrease at 6 months (sensitivity = 60%; McNemar's test, P = 0.06). Conversely, for those women with a measured decrease in BD from baseline to 6 months, 9 of 10 had a measured decrease at 12 months. A study limitation is inclusion of participants on TAM for varying duration as the greatest change in BD likely would have occurred earlier. Ongoing efforts will focus on FWD-MRI for measures of change in BD in patients initiating TAM.

Conclusion: Use of the specified cut point would fail to detect a decrease in BD at 12 months in 40% of women. However, a decrease in BD from baseline to 6 months was highly associated with decrease from baseline to 12 months and in some women may be useful as an early biomarker of effect. Ongoing effort is needed to determine the impact of factors such as baseline BD, menopausal status, and time on TAM in misclassification of BD change using the 6-month measure.

Acknowledgement: NIH grants CA149417, CA161534.Objective: Tamoxifen (TAM) lowers breast cancer recurrence by 40-50% with evidence of individual variability in responsiveness. A ≥10% decrease in mammography-determined breast density (BD) after 12–18 months of TAM use has been associated with clinical benefit. Early determination of changes in BD may offer a strategy to tailor hormone therapy in non-responders; for responders, it may encourage adherence. Fat-water decomposition MRI (FWD-MRI) is an accurate and fast (&lt; 5 minutes) method for measuring BD without ionizing radiation or contrast agent. Here, we examined whether change in FWD-MRI-derived BD predicts decrease in BD at earlier time points than observable with a 12-month measure of BD.

Methods: The study population included a subset of 44 pre- and post-menopausal women receiving TAM for treatment of early-stage breast cancer or prevention who were enrolled in a randomized, placebo-controlled trial of diindolylmethane. Eligibility for this analysis included participants with FWD-MRI scans at baseline, 6 and 12 months. Median time on TAM at baseline was 13 months (IQR, 5–26 months). All MRI images were acquired on a 1.5T GE Signa NV-CV/i scanner. Automated breast segmentation was performed using MATLAB software and validated against manual ROI drawings. MRI-based BD was calculated as the ratio of breast voxels with &lt;80% apparent fat fraction (Fra80) over the entire breast, a measure previously shown by our group to be highly correlated with mammography-derived BD. For 40 participants, the unaffected, contralateral breast was analyzed. For 4 patients with two unaffected breasts, BD data from the left breast were analyzed. Change in BD was conservatively defined as &gt; 2 times the test-retest variability of Fra80 (0.032). McNemar's test was used to test the association between change from baseline to 6 months and change from baseline to 12 months.

Results and Discussion: At 12 months, 15 (34%) participants had a decrease in BD, whereas 29 (66%) remained unchanged or increased. Of these 29, 28 also had no decrease at 6 months (specificity = 97%), and 9 of the 15 women who showed a decrease at 12 months had a decrease at 6 months (sensitivity = 60%; McNemar's test, P = 0.06). Conversely, for those women with a measured decrease in BD from baseline to 6 months, 9 of 10 had a measured decrease at 12 months. A study limitation is inclusion of participants on TAM for varying duration as the greatest change in BD likely would have occurred earlier. Ongoing efforts will focus on FWD-MRI for measures of change in BD in patients initiating TAM.

Conclusion: Use of the specified cut point would fail to detect a decrease in BD at 12 months in 40% of women. However, a decrease in BD from baseline to 6 months was highly associated with decrease from baseline to 12 months and in some women may be useful as an early biomarker of effect. Ongoing effort is needed to determine the impact of factors such as baseline BD, menopausal status, and time on TAM in misclassification of BD change using the 6-month measure.

Acknowledgement: NIH grants CA149417, CA161534.

Citation Format: Ding J, Thompson PA, Wertheim BC, Roe DJ, Marron MT, Altbach MI, Galons J-P, Wang F, Thomson CA, Huang C, Stopeck A. Breast density change at 6 months is associated with change at 12 months as measured by fat-water decomposition MRI in women on tamoxifen [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr P6-09-19.

Abstract P3-02-03: Accurate and reliable automated breast density measurements with no ionizing radiation using fat-water decomposition MRI

Objective Breast density(BD) is a measure of the distribution of variable tissue types within the breast and higher BD has been shown to positively correlate with breast cancer risk. As such, the accurate measurement of BD has become a priority for risk assessment and for evaluating the effects of prevention strategies aimed at reducing BD. Mammography(MG) is the most common method of BD determination but is limited by the exposure to ionizing radiation, particularly for studies requiring repeated measures. BD derived from fat-water decomposition magnetic resonance imaging(FWMRI-BD) has been proposed as an alternative, safe, and quantitative method for BD. To optimize its use, we developed a new FWMRI-BD that is automated, more accurate and reliable. In this study, we compare our automated method to digital MG and a previous reported algorithm for MRI derived BD.

Methods From a completed prevention trial, 42 pre- and post-menopausal patients receiving tamoxifen therapy for early stage breast cancer or as primary chemoprevention were identified. Patients had undergone prior digital MG within 6 months from the date of MRI scan and MG-BD was calculated using a well-established method(Cumulus). MRI scans were performed on a 1.5T GE Signa NV-CV/i scanner using an axial radial IDEAL-GRASE sequence to generate quantitative fat fraction maps of the entire breast. Total acquisition time was &lt; 5 min and automated breast segmentation was applied to all scans. Only the contralateral, unaffected breast was analyzed. Pearson correlation analysis compared BD as measured by MG(range 0-100%) and FWMRI based methods. BD by FWMRI was initially calculated as the ratio of breast voxels with&lt;80% apparent fat fraction(Fra80). Fra80 had been previously shown by our group to correlate with MG-BD(Spearman ρ=0.86, p&lt;0.001). Here, BD was calculated using a new algorithm(FraG+W) that accounts for the total amount of fibroglandular tissue and water content in the breast after correction for fat-water signal intensity bias and fat-water signal shine-through. Reliability of FWMRI measurements was tested in 24 repeated scans from 9 patients and evaluated using intra-class correlation(ICC) analysis.

Results Table 1 shows the correlation and reliability analysis results between MG-BD and FWMRI-BD. Both FWMRI-BD measures(Fra80 and FraG+W) were strongly correlated with MG-BD. More importantly, they exhibit superior test-retest reliability(ICC&gt;0.98) compared to MG-BD values from the literature(reported ICC range 0.91-0.95). FraG+W showed improvement over Fra80 in all measures tested including correlation to MG-BD, dynamic range, standard errors and ICC.

Table 1. Accuracy and Reliability of the FWMRI-BD measuresFWMRI-BDFra80FraG+WPearson correlation coefficient* with MG-BDR=0.86R=0.94Test-retest reliabilitystandard error0.02300.0134dynamic range0.0902 – 0.65370.0736 – 0.6588standard error/ dynamic range4.1%2.3%ICC [95% confidence interval]0.985 [0.966,0.993]0.990 [0.976,0.995]* All P-values &lt; 1e-10

Conclusion The refined and automated FWMRI-BD that quantifies the entire fibroglandular and water content of the breast(FraG+W) strongly correlates with MG-BD and is more accurate and reliable than previous FWMRI-BD method.

Acknowledgement NIH grants CA149417, CA161534.

Citation Format: Ding J, Thompson PA, Gao Y, Marron MT, Wertheim BC, Altbach MI, Galons J-P, Roe DJ, Wang F, Maskarinec G, Thomson CA, Stopeck A, Huang C. Accurate and reliable automated breast density measurements with no ionizing radiation using fat-water decomposition MRI [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr P3-02-03.

Multiresolution imaging using golden angle stack-of-stars and compressed sensing for dynamic MR urography

PURPOSE

To develop a novel multiresolution MRI methodology for accurate estimation of glomerular filtration rate (GFR) in vivo.

MATERIALS AND METHODS

A three-dimensional golden-angle radial stack-of-stars (SoS) trajectory was used for data acquisition on a 3 Tesla MRI scanner. Multiresolution reconstruction and analysis was performed using arterial input function reconstructed at 1-s. temporal resolution and renal dynamic data reconstructed using compressed sensing (CS) with 4-s temporal resolution. The method was first validated using simulations and the clinical utility of the technique was evaluated by comparing the GFR estimates from the proposed method to the estimated GFR (eGFR) obtained from serum creatinine for 10 subjects.

RESULTS

The 4-s temporal resolution CS images minimized streaking artifacts and noise while the 1-s temporal resolution AIF minimized errors in GFR estimates. A paired t-test showed that there was no statistically significant difference between MRI based total GFR values and serum creatinine based eGFR estimates (P = 0.92).

CONCLUSION

We have demonstrated the feasibility of multiresolution MRI using a golden angle radial stack-of-stars scheme to accurately estimate GFR as well as produce diagnostic quality dynamic images in vivo.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 3 J. MAGN. RESON. IMAGING 2017;46:303-311.

Task-based optimization of flip angle for fibrosis detection in T1-weighted MRI of liver

Chronic liver disease is a worldwide health problem, and hepatic fibrosis (HF) is one of the hallmarks of the disease. The current reference standard for diagnosing HF is biopsy followed by pathologist examination; however, this is limited by sampling error and carries a risk of complications. Pathology diagnosis of HF is based on textural change in the liver as a lobular collagen network that develops within portal triads. The scale of collagen lobules is characteristically in the order of 1 to 5 mm, which approximates the resolution limit of in vivo gadolinium-enhanced magnetic resonance imaging in the delayed phase. We use MRI of formalin-fixed human ex vivo liver samples as phantoms that mimic the textural contrast of in vivo Gd-MRI. We have developed a local texture analysis that is applied to phantom images, and the results are used to train model observers to detect HF. The performance of the observer is assessed with the area-under-the-receiver-operator-characteristic curve (AUROC) as the figure-of-merit. To optimize the MRI pulse sequence, phantoms were scanned with multiple times at a range of flip angles. The flip angle that was associated with the highest AUROC was chosen as optimal for the task of detecting HF.

T2 mapping of the heart with a double-inversion radial fast spin-echo method with indirect echo compensation

BACKGROUND

The abnormal signal intensity in cardiac T2-weighted images is associated with various pathologies including myocardial edema. However, the assessment of pathologies based on signal intensity is affected by the acquisition parameters and the sensitivities of the receiver coils. T2 mapping has been proposed to overcome limitations of T2-weighted imaging, but most methods are limited in spatial and/or temporal resolution. Here we present and evaluate a double inversion recovery radial fast spin-echo (DIR-RADFSE) technique that yields data with high spatiotemporal resolution for cardiac T2 mapping.

METHODS

DIR-RADFSE data were collected at 1.5 T on phantoms and subjects with echo train length (ETL) = 16, receiver bandwidth (BW) = ±32 kHz, TR = 1RR, matrix size = 256 × 256. Since only 16 views per echo time (TE) are collected, two algorithms designed to reconstruct highly undersampled radial data were used to generate images for 16 time points: the Echo-Sharing (ES) and the CUrve Reconstruction via pca-based Linearization with Indirect Echo compensation (CURLIE) algorithm. T2 maps were generated via least-squares fitting or the Slice-resolved Extended Phase Graph (SEPG) model fitting. The CURLIE-SEPG algorithm accounts for the effect of indirect echoes. The algorithms were compared based on reproducibility, using Bland-Altman analysis on data from 7 healthy volunteers, and T2 accuracy (against a single-echo spin-echo technique) using phantoms.

RESULTS

Both reconstruction algorithms generated in vivo images with high spatiotemporal resolution and showed good reproducibility. Mean T2 difference between repeated measures and the coefficient of repeatability were 0.58 ms and 2.97 for ES and 0.09 ms and 4.85 for CURLIE-SEPG. In vivo T2 estimates from ES were higher than those from CURLIE-SEPG. In phantoms, CURLIE-SEPG yielded more accurate T2s compared to reference values (error was 7.5-13.9% for ES and 0.6-2.1% for CURLIE-SEPG), consistent with the fact that CURLIE-SEPG compensates for the effects of indirect echoes. The potential of T2 mapping with CURLIE-SEPG is demonstrated in two subjects with known heart disease. Elevated T2 values were observed in areas of suspected pathology.

CONCLUSIONS

DIR-RADFSE yielded TE images with high spatiotemporal resolution. Two algorithms for generating T2 maps from highly undersampled data were evaluated in terms of accuracy and reproducibility. Results showed that CURLIE-SEPG yields T2 estimates that are reproducible and more accurate than ES.

Correcting partial volume effects in biexponential T2 estimation of small lesions

PURPOSE

T2 mapping provides a quantitative approach for focal liver lesion characterization. For small lesions, a biexponential model should be used to account for partial volume effects (PVE). However, conventional biexponential fitting suffers from large uncertainty of the fitted parameters when noise is present. The purpose of this work is to develop a more robust method to correct for PVE affecting small lesions.

METHODS

We developed a region of interest-based joint biexponential fitting (JBF) algorithm to estimate the T2 of lesions affected by PVE. JBF takes advantage of the lesion fraction variation among voxels within a region of interest. JBF is compared to conventional approaches using Cramér-Rao lower bound analysis, numerical simulations, phantom, and in vivo data.

RESULTS

JBF provides more accurate and precise T2 estimates in the presence of PVE. Furthermore, JBF is less sensitive to region of interest drawing. Phantom and in vivo results show that JBF can be combined with a reconstruction method for highly undersampled data, enabling the characterization of small abdominal lesions from data acquired in a single breath hold.

CONCLUSION

The JBF algorithm provides more accurate and stable T2 estimates for small structures than conventional techniques when PVE is present. It should be particularly useful for the characterization of small abdominal lesions.

T2 relaxometry with indirect echo compensation from highly undersampled data

PURPOSE

To develop an algorithm for fast and accurate T2 estimation from highly undersampled multi-echo spin-echo data.

METHODS

The algorithm combines a model-based reconstruction with a signal decay based on the slice-resolved extended phase graph (SEPG) model with the goal of reconstructing T2 maps from highly undersampled radial multi-echo spin-echo data with indirect echo compensation. To avoid problems associated with the nonlinearity of the SEPG model, principal component decomposition is used to linearize the signal model. The proposed CUrve Reconstruction via principal component-based Linearization with Indirect Echo compensation (CURLIE) algorithm is used to estimate T2 curves from highly undersampled data. T2 maps are obtained by fitting the curves to the SEPG model.

RESULTS

Results on phantoms showed T2 biases (1.9% to 18.4%) when indirect echoes are not taken into account. The T2 biases were reduced (< 3.2%) when the CURLIE reconstruction was performed along with SEPG fitting even for high degrees of undersampling (4% sampled). Experiments in vivo for brain, liver, and heart followed the same trend as the phantoms.

CONCLUSION

The CURLIE reconstruction combined with SEPG fitting enables accurate T2 estimation from highly undersampled multi-echo spin-echo radial data thus, yielding a fast T2 mapping method without errors caused by indirect echoes.

T2 mapping from highly undersampled data by reconstruction of principal component coefficient maps using compressed sensing

Recently, there has been an increased interest in quantitative MR parameters to improve diagnosis and treatment. Parameter mapping requires multiple images acquired with different timings usually resulting in long acquisition times. While acquisition time can be reduced by acquiring undersampled data, obtaining accurate estimates of parameters from undersampled data is a challenging problem, in particular for structures with high spatial frequency content. In this work, principal component analysis is combined with a model-based algorithm to reconstruct maps of selected principal component coefficients from highly undersampled radial MRI data. This novel approach linearizes the cost function of the optimization problem yielding a more accurate and reliable estimation of MR parameter maps. The proposed algorithm--reconstruction of principal component coefficient maps using compressed sensing--is demonstrated in phantoms and in vivo and compared with two other algorithms previously developed for undersampled data.

Automated registration of sequential breath-hold dynamic contrast-enhanced MR images: a comparison of three techniques

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) is increasingly in use as an investigational biomarker of response in cancer clinical studies. Proper registration of images acquired at different time points is essential for deriving diagnostic information from quantitative pharmacokinetic analysis of these data. Motion artifacts in the presence of time-varying intensity due to contrast enhancement make this registration problem challenging. DCE-MRI of chest and abdominal lesions is typically performed during sequential breath-holds, which introduces misregistration due to inconsistent diaphragm positions and also places constraints on temporal resolution vis-à-vis free-breathing. In this work, we have employed a computer-generated DCE-MRI phantom to compare the performance of two published methods, Progressive Principal Component Registration and Pharmacokinetic Model-Driven Registration, with Sequential Elastic Registration (SER) to register adjacent time-sample images using a published general-purpose elastic registration algorithm. In all three methods, a 3D rigid-body registration scheme with a mutual information similarity measure was used as a preprocessing step. The DCE-MRI phantom images were mathematically deformed to simulate misregistration, which was corrected using the three schemes. All three schemes were comparably successful in registering large regions of interest (ROIs) such as muscle, liver, and spleen. SER was superior in retaining tumor volume and shape, and in registering smaller but important ROIs such as tumor core and tumor rim. The performance of SER on clinical DCE-MRI data sets is also presented.

Mild hypothermia delays the development of stone heart from untreated sustained ventricular fibrillation--a cardiovascular magnetic resonance study

BACKGROUND

'Stone heart' resulting from ischemic contracture of the myocardium, precludes successful resuscitation from ventricular fibrillation (VF). We hypothesized that mild hypothermia might slow the progression to stone heart.

METHODS

Fourteen swine (27 ± 1 kg) were randomized to normothermia (group I; n=6) or hypothermia groups (group II; n=8). Mild hypothermia (34 ± 2 °C) was induced with ice packs prior to VF induction. The LV and right ventricular (RV) cross-sectional areas were followed by cardiovascular magnetic resonance until the development of stone heart. A commercial 1.5T GE Signa NV-CV/i scanner was used. Complete anatomic coverage of the heart was acquired using a steady-state free precession (SSFP) pulse sequence gated at baseline prior to VF onset. Un-gated SSFP images were obtained serially after VF induction. The ventricular endocardium was manually traced and LV and RV volumes were calculated at each time point.

RESULTS

In group I, the LV was dilated compared to baseline at 5 minutes after VF and this remained for 20 minutes. Stone heart, arbitrarily defined as LV volume <1/3 of baseline at the onset of VF, occurred at 29 ± 3 minutes. In group II, there was less early dilation of the LV (p<0.05) and the development of stone heart was delayed to 52 ± 4 minutes after onset of VF (P<0.001).

CONCLUSIONS

In this closed-chest swine model of prolonged untreated VF, hypothermia reduced the early LV dilatation and importantly, delayed the onset of stone heart thereby extending a known, morphologic limit of resuscitability.

Rapid water and lipid imaging with T2 mapping using a radial IDEAL-GRASE technique

Three-point Dixon methods have been investigated as a means to generate water and fat images without the effects of field inhomogeneities. Recently, an iterative algorithm (IDEAL, iterative decomposition of water and fat with echo asymmetry and least squares estimation) was combined with a gradient and spin-echo acquisition strategy (IDEAL-GRASE) to provide a time-efficient method for lipid-water imaging with correction for the effects of field inhomogeneities. The method presented in this work combines IDEAL-GRASE with radial data acquisition. Radial data sampling offers robustness to motion over Cartesian trajectories as well as the possibility of generating high-resolution T(2) maps in addition to the water and fat images. The radial IDEAL-GRASE technique is demonstrated in phantoms and in vivo for various applications including abdominal, pelvic, and cardiac imaging.

Cardiac MRI is Complementary to Echocardiography in the Assessment of Cardiac Masses

Despite the fact that the incidence of cardiac tumors is low, the prompt evaluation and adequate intervention of these is highly important. Although most tumors of the heart are considered histologically benign, there are significant risks associated with these "benign" tumors. These are associated with significant morbidity and mortality due to obstruction of blood flow, alterations of conduction, propagation of arrhythmias, and thromboembolism, depending on their size, location, and nature. With the advent of noninvasive imaging modalities--traditionally echocardiography; but more recently using cross-sectional imaging with cardiac computed tomography and magnetic resonance imaging--cardiac tumors can be optimally assessed providing a greater opportunity for curative treatments by cardiothoracic surgery.

Cardiac MRI Is an Important Complementary Tool to Doppler Echocardiography in the Management of Patients with Pulmonary Regurgitation

Cardiac MRI (CMR) is a noninvasive diagnostic tool with comprehensive capabilities similar to that of two-dimensional echocardiography with Doppler. In addition to the ability to evaluate the etiology and severity of pulmonary valve regurgitation (PR), CMR is well designed to serially monitor the impact of the PR on the right ventricle (RV). Importantly, RV dilation and dysfunction is a critical determinate to time surgical intervention. CMR gives the silent RV, suffering from PR, a voice.

Fast decomposition of water and lipid using a GRASE technique with the IDEAL algorithm

Three-point Dixon techniques achieve good lipid-water separation by estimating the phase due to field inhomogeneities. Recently it was demonstrated that the combination of an iterative algorithm (iterative decomposition of water and fat with echo asymmetry and least-squares estimation (IDEAL)) with a fast spin-echo (FSE) three-point Dixon method yielded robust lipid-water decomposition. As an alternative to FSE, the gradient- and spin-echo (GRASE) technique has been developed for efficient data collection. In this work we present a method for lipid-water separation by combining IDEAL with the GRASE technique. An approach to correct for errors in the lipid-water decomposition caused by phase distortions due to the switching of the readout gradient polarities inherent to GRASE is presented. The IDEAL-GRASE technique is demonstrated in phantoms and in vivo for various applications, including pelvic, musculoskeletal, and (breath-hold) cardiac imaging.

Isotropic diffusion weighting in radial fast spin-echo magnetic resonance imaging

Radial fast spin-echo (radial-FSE) methods enable multishot diffusion-weighted MRI (DWMRI) to be carried out without significant artifacts due to motion and/or susceptibility and can be used to generate DWMRI images with high spatial resolution. In this work, a novel method that allows isotropic diffusion weighting to be obtained in a single radial k-space data set is presented. This is accomplished by altering the direction of diffusion weighting gradients between groups of TR periods, which yield sets of radial lines that possess diffusion weighting sensitive to motion in different directions. By altering the diffusion weighting directions and controlling the view ordering appropriately within the sequence, an effectively isotropic diffusion-weighted image can be obtained within one radial-FSE scan. The order in which radial lines are acquired can also be controlled to yield data sets without significant artifacts due to motion, T(2) decay, and/or diffusion anisotropy.

Magnetic resonance imaging during untreated ventricular fibrillation reveals prompt right ventricular overdistention without left ventricular volume loss

BACKGROUND

Most out-of-hospital ventricular fibrillation (VF) is prolonged (>5 minutes), and defibrillation from prolonged VF typically results in asystole or pulseless electrical activity. Recent visual epicardial observations in an open-chest, open-pericardium model of swine VF indicate that blood flows from the high-pressure arterial system to the lower-pressure venous system during untreated VF, thereby overdistending the right ventricle and apparently decreasing left ventricular size. Therefore, inadequate left ventricular stroke volume after defibrillation from prolonged VF has been postulated as a major contributor to the development of pulseless rhythms.

METHODS AND RESULTS

Ventricular dimensions were determined by MRI for 30 minutes of untreated VF in a closed-chest, closed-pericardium model in 6 swine. Within 1 minute of untreated VF, mean right ventricular volume increased by 29% but did not increase thereafter. During the first 5 minutes of untreated VF, mean left ventricular volume increased by 34%. Between 20 and 30 minutes of VF, stone heart occurred as manifested by dramatic thickening of the myocardium and concomitant substantial decreases in left ventricular volume.

CONCLUSIONS

In this closed-chest swine model of VF, substantial right ventricular volume changes occurred early and did not result in smaller left ventricular volumes. The changes in ventricular volumes before the late development of stone heart do not explain why defibrillation from brief duration VF (<5 minutes) typically results in a pulsatile rhythm with return of spontaneous circulation, whereas defibrillation from prolonged VF (5 to 15 minutes) does not.

Radial GRASE: Implementation and applications

RAD-GRASE is an MRI sequence that combines radial (RAD) k-space scanning with the gradient and spin-echo (GRASE) technique. RAD-GRASE has the advantages of all radial data acquisition methods in that it can reduce motion sensitivity and correct motion-induced data errors, which can be exploited to achieve high-resolution diffusion-weighted imaging (DWI). One can obtain different types of image contrast, including DWI, T(1), T(2), and T(2)*, in RAD-GRASE by controlling the magnetization preparation and sequence timing. Moreover, because there is oversampling of the low spatial frequencies inherent to radial sequences, partial data reconstruction can be used to achieve multiple forms of image contrast from a single acquired data set, and to generate parametric image maps of equilibrium magnetization, T(2), and T(2) (dagger). The RAD-GRASE technique can also be used to achieve fat-suppressed and/or separated fat and water images by choosing the appropriate timing parameters.

Interleaved acquisition of lipid and water images of the heart using a double-inversion fast spin-echo method

In this work we present a new method for the improved detection of lipid infiltration in the heart. The method employs a double-inversion fast spin-echo technique where the acquisition of water- and lipid-suppressed k-space data is alternated between TR periods to produce co-registered lipid and water images from data acquired in a breath hold. The lipid and water images can then be combined to generate a lipid/water image with reduced artifacts due to flow and excellent contrast between lipid and myocardium. The method is demonstrated in ex vivo tissue and in vivo. This novel method may improve the detection of lipid infiltration in the heart in pathologies such as arrhythmogenic right ventricular dysplasia.

Processing of radial fast spin-echo data for obtaining T2 estimates from a single k-space data set

Radially acquired fast spin-echo data can be processed to obtain T2-weighted images and a T2 map from a single k-space data set. The general approach is to use data at a specific TE (or narrow TE range) in the center of k-space and data at other TE values in the outer part of k-space. With this method high-resolution T2-weighted images and T2 maps are obtained in a time efficient manner. The mixing of TE data, however, introduces errors in the T2-weighted images and T2 maps that affect the accuracy of the T2 estimates. In this work, various k-space data processing methods for reconstructing T2-weighted images and T2 maps from a single radial fast spin-echo k-space data set are analyzed in terms of the accuracy of T2 estimates. The analysis is focused on the effect of image artifacts, object dependency, and noise on the T2 estimates. Results are presented in computer-generated phantoms and in vivo.

View-ordering in radial fast spin-echo imaging

Radial MRI sequences are frequently used to obtain images with reduced sensitivity to motion. To decrease imaging time, multiple spin-echo acquisitions can be incorporated into radial sequences. In this case, different radial lines of Fourier data have different TE times and the resulting images can contain streaking artifacts due to T(2) decay. The streaking is not only dependent on the T(2) of the object and the timing of the data acquisition, but also on the order in which radial lines are collected (view order). The view ordering can easily be controlled to minimize artifacts due to T(2) decay as well as motion. Four view-ordering techniques are presented and evaluated for the radial FSE sequence.

Radial fast spin-echo method for T2-weighted imaging and T2 mapping of the liver

PURPOSE

To evaluate a multishot radial fast-spin echo (RAD-FSE) method developed to improve the quality of abdominal T2-weighted imaging as well as the characterization of focal liver lesions.

MATERIALS AND METHODS

The RAD-FSE sequence used in this work consisted of a preparatory period followed by a short echo train (ETL = 16). A novel radial k-space trajectory was used to minimize streaking artifacts due to T2 variations and motion. Small diffusion gradients (b = 1.2 mm/s(2)) were used to improve flow suppression. The quality of images obtained with RAD-FSE was compared to multishot 2DFT fast spin-echo (2DFT-FSE) and half-Fourier acquisition single-shot turbo-spin-echo (HASTE) images using data from 16 patients. A postprocessing algorithm was used to generate multiple high-resolution images (at different effective TE values) as well as a T2 map from a single RAD-FSE data set. The T2 maps were used to differentiate malignant from benign lesions for a set of 33 lesions ranging from 0.8-194 cm(3).

RESULTS

RAD-FSE produces high-resolution images of the liver in a breath-hold without the motion artifacts of 2DFT-FSE methods, and without the blurriness and loss of small lesion detectability of HASTE. The inclusion of diffusion weighting in RAD-FSE decreases the signal from blood in hepatic vessels, which improves lesion visualization. The T2 values obtained by postprocessing a single RAD-FSE data set can differentiate malignant from benign lesions. The mean T2 values obtained for malignancies, hemangiomas, and cysts are 108 +/- 30 msec, 240 +/- 14 msec, and 572 +/- 334 msec, respectively.

CONCLUSION

These results indicate that RAD-FSE produces abdominal images of higher quality than 2DFT-FSE and HASTE. In addition, lesions can be characterized using T2 maps generated from a single RAD-FSE data set.

Chemical-shift imaging utilizing the positional shifts along the readout gradient direction

In this work, we describe a method that uses the linear phase acquired during the readout period due to chemical shift to generate individual magnetic resonance (MR) images of chemically shifted species. The method utilizes sets of Fourier (or k-space) data acquired with different directions of the readout gradient and a postprocessing algorithm to generate chemical shift images. The methodology is developed for both Cartesian data acquisition and for radial data acquisition. The method is presented here for two chemically shifted species but it can be extended to more species. In this work, we present the theory, show the results in phantoms and in human images, and discuss the artifacts and signal-to-noise ratio of the images obtained with the technique.

Reconstruction of MR images from data acquired on a general nonregular grid by pseudoinverse calculation

A minimum-norm least-squares image-reconstruction method for the reconstruction of magnetic resonance images from non-Cartesian sampled data is proposed. The method is based on a general formalism for continuous-to-discrete mapping and pseudoinverse calculation. It does not involve any regridding or interpolation of the data and therefore the methodology differs fundamentally from existing regridding-based methods. Moreover, the method uses a continuous representation of objects in the image domain instead of a discretized representation. Simulations and experiments show the possibilities of the method in both radial and spiral imaging. Simulations revealed that minimum-norm least-squares image reconstruction can result in a drastic decrease of artifacts compared with regridding-based reconstruction. Besides, both in vivo and phantom experiments showed that minimum-norm least-squares image reconstruction leads to contrast improvement and increased signal-to-noise ratio compared with image reconstruction based on regridding. As an appendix, an analytical calculation of the raw data corresponding to the well-known Shepp and Logan software head phantom is presented.

High-resolution diffusion imaging with DIFRAD-FSE (diffusion-weighted radial acquisition of data with fast spin-echo) MRI

A novel MRI method, DIFRAD-FSE (diffusion-weighted radial acquisition of data with fast spin-echo), is demonstrated that enables rapid, high-resolution multi-shot diffusion-weighted MRI without significant artifacts due to motion. Following a diffusion-weighting spin-echo preparation period, multiple radial lines of Fourier data are acquired using spin-echo refocusing. Images can be reconstructed from the radial data set using a magnitude-only filtered back-projection reconstruction algorithm that removes phase errors due to motion. Results from human brain imaging demonstrate the ability of DIFRAD-FSE to acquire multiple radial lines of Fourier data each TR period without significant artifacts due to relaxation and to produce high-resolution diffusion-weighted MRI images without significant artifacts from motion.

Early increases in breast tumor xenograft water mobility in response to paclitaxel therapy detected by non-invasive diffusion magnetic resonance imaging

An important goal in cancer chemotherapy is to sensitively and quantitatively monitor the response of individual patients' tumors to successful, or unsuccessful, therapy so that regimens can be altered iteratively. Currently, tumor response is monitored by frank changes in tumor morphology, yet these markers take long to manifest and are not quantitative. Recent studies suggest that the apparent diffusion coefficient of water (ADCw), measured noninvasively with magnetic resonance imaging, is sensitively and reliably increased in response to successful CTx. In the present study, we investigate the combination chemotherapy response of human breast cancer tumor xenografts sensitive or resistant to Paclitaxel by monitoring changes in the ADCw. Our results indicate that there is a clear, substantial, and early increase in the ADCw after successful therapy in drug sensitive tumors and that there is no change in the ADCw in p-glycoprotein-positive tumors, which are resistant to Paclitaxel. The mechanism underlying these changes is unknown yet is consistent with apoptotic cell shrinkage and a concomitant increase in the extracellular water fraction.

Plasmalemmal pH-gradients in drug-sensitive and drug-resistant MCF-7 human breast carcinoma xenografts measured by 31P magnetic resonance spectroscopy

31p Magnetic resonance spectroscopy (MRS) was employed to investigate tumor pH in xenografts of drug-sensitive and drug-resistant MCF-7 human breast carcinoma cells. Measured extracellular pH values were found to be lower than the intracellular pH in all three tumor types investigated. The magnitude of this acid-outside plasmalemmal pH gradient increased with increasing tumor size in tumors of two drug-resistant variants of MCF-7 cells, but not in tumors of the parent (drug-sensitive) cells. The partitioning of weak-base or weak-acid drug molecules across the plasma membrane of a tumor cell is dependent upon the acid-dissociation constant (pKa) of the drug as well as the plasmalemmal pH gradient. A large acid-outside pH gradient, such as those seen in MCF-7 xenografts, can exert a protective effect on the cell from weak-base drugs such as anthracyclines and Vinca alkaloids, which have pKa values of 7.5 to 9.5. The possibility of enhancing the therapeutic efficacy of weak-base drugs by dietary or metabolic manipulation of the extracellular pH, in order to reduce or reverse the plasmalemmal pH gradient, deserves investigation.

MRI and NMR spectroscopy of the lipids of atherosclerotic plaque in rabbits and humans

The early stages of atherosclerosis are characterized by the deposition of cholesteryl esters and triglycerides into the arterial wall. In the excised human atherosclerotic plaque these lipids are in a liquid-like state at body temperature and observable via MRI and NMR spectroscopy. To assess the ability of MRI to quantitatively image the lipids of atherosclerotic plaque in vivo, we have investigated eight New Zealand White rabbits fed atherogenic diets (2 weight (wt)% cholesterol, 1 wt% cholesterol + 6 wt% peanut oil, and 1 wt% cholesterol + 6 wt% com oil). Postmortem examination indicated that all rabbits developed atherosclerosis in the aorta. Except for one animal, magnetic resonance angiography showed no noticeable obstruction in the aorta. MRI was carried out in an attempt to image atherosclerotic plaque lipids directly, but no signal was detected in vivo. However, a plaque lipid signal was observed from excised tissue using a small diameter RF coil. 1H NMR spectroscopy of the atherosclerotic plaque from excised aortas indicated that the major fraction of plaque lipids in rabbits is not in a liquid state at physiological temperature and are only marginally MRI-visible compared to human plaque lipid. The differences in the MRI characteristics of rabbit and human plaque are due to differences in the fatty acid profile of the cholesteryl esters, chiefly a decrease of linoleic acid in rabbit lesions.