Alternating low-rank tensor reconstruction for improved multiparametric mapping with cardiovascular MR Multitasking

PURPOSE

To develop a novel low-rank tensor reconstruction approach leveraging the complete acquired data set to improve precision and repeatability of multiparametric mapping within the cardiovascular MR Multitasking framework.

METHODS

A novel approach that alternated between estimation of temporal components and spatial components using the entire data set acquired (i.e., including navigator data and imaging data) was developed to improve reconstruction. The precision and repeatability of the proposed approach were evaluated on numerical simulations, 10 healthy subjects, and 10 cardiomyopathy patients at multiple scan times for 2D myocardial T1/T2 mapping with MR Multitasking and were compared with those of the previous navigator-derived fixed-basis approach.

RESULTS

In numerical simulations, the proposed approach outperformed the previous fixed-basis approach with lower T1 and T2 error against the ground truth at all scan times studied and showed better motion fidelity. In human subjects, the proposed approach showed no significantly different sharpness or T1/T2 measurement and significantly improved T1 precision by 20%-25%, T2 precision by 10%-15%, T1 repeatability by about 30%, and T2 repeatability by 25%-35% at 90-s and 50-s scan times The proposed approach at the 50-s scan time also showed comparable results with that of the previous fixed-basis approach at the 90-s scan time.

CONCLUSION

The proposed approach improved precision and repeatability for quantitative imaging with MR Multitasking while maintaining comparable motion fidelity, T1/T2 measurement, and septum sharpness and had the potential for further reducing scan time from 90 s to 50 s.

Dynamic Regularized Adaptive Cluster Optimization (DRACO) for Quantitative Cardiac Cine MRI in Complex Arrhythmias

BACKGROUND

Irregular cardiac motion can render conventional segmented cine MRI nondiagnostic. Clustering has been proposed for cardiac motion binning and may be optimized for complex arrhythmias.

PURPOSE

To develop an adaptive cluster optimization method for irregular cardiac motion, and to generate the corresponding time-resolved cine images.

STUDY TYPE

Prospective.

SUBJECTS

Thirteen with atrial fibrillation, four with premature ventricular contractions, and one patient in sinus rhythm.

FIELD STRENGTH/SEQUENCE

Free-running balanced steady state free precession (bSSFP) with sorted golden-step, reference real-time sequence.

ASSESSMENT

Each subject underwent both the sorted golden-step bSSFP and the reference Cartesian real-time imaging. Golden-step bSSFP images were reconstructed using the dynamic regularized adaptive cluster optimization (DRACO) method and k-means clustering. Image quality (4-point Likert scale), signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), edge sharpness, and ventricular function were assessed.

STATISTICAL TESTS

Paired t-tests, Friedman test, regression analysis, Fleiss' Kappa, Bland-Altman analysis. Significance level P < 0.05.

RESULTS

The DRACO method had the highest percent of images with scores ≥3 (96% for diastolic frame, 93% for systolic frame, and 93% for multiphase cine) and the percentages were significantly higher compared with both the k-means and real-time methods. Image quality scores, SNR, and CNR were significantly different between DRACO vs. k-means and between DRACO vs. real-time. Cardiac function analysis showed no significant differences between DRACO vs. the reference real-time.

CONCLUSION

DRACO with time-resolved reconstruction generated high quality images and has early promise for quantitative cine cardiac MRI in patients with complex arrhythmias including atrial fibrillation.

LEVEL OF EVIDENCE: 1

TECHNICAL EFFICACY

Stage 2.

Non-electrocardiogram-gated, free-breathing, off-resonance reduced, high-resolution, whole-heart myocardial T2 * mapping at 3 T within 5 min

PURPOSE

Widely used conventional 2D T2 * approaches that are based on breath-held, electrocardiogram (ECG)-gated, multi-gradient-echo sequences are prone to motion artifacts in the presence of incomplete breath holding or arrhythmias, which is common in cardiac patients. To address these limitations, a 3D, non-ECG-gated, free-breathing T2 * technique that enables rapid whole-heart coverage was developed and validated.

METHODS

A continuous random Gaussian 3D k-space sampling was implemented using a low-rank tensor framework for motion-resolved 3D T2 * imaging. This approach was tested in healthy human volunteers and in swine before and after intravenous administration of ferumoxytol.

RESULTS

Spatial-resolution matched T2 * images were acquired with 2-3-fold reduction in scan time using the proposed T2 * mapping approach relative to conventional T2 * mapping. Compared with the conventional approach, T2 * images acquired with the proposed method demonstrated reduced off-resonance and flow artifacts, leading to higher image quality and lower coefficient of variation in T2 *-weighted images of the myocardium of swine and humans. Mean myocardial T2 * values determined using the proposed and conventional approaches were highly correlated and showed minimal bias.

CONCLUSION

The proposed non-ECG-gated, free-breathing, 3D T2 * imaging approach can be performed within 5 min or less. It can overcome critical image artifacts from undesirable cardiac and respiratory motion and bulk off-resonance shifts at the heart-lung interface. The proposed approach is expected to facilitate faster and improved cardiac T2 * mapping in those with limited breath-holding capacity or arrhythmias.

ECG-free cine MRI with data-driven clustering of cardiac motion for quantification of ventricular function

BACKGROUND

Despite the widespread use of cine MRI for evaluation of cardiac function, existing real-time methods do not easily enable quantification of ventricular function. Moreover, segmented cine MRI assumes periodicity of cardiac motion. We aim to develop a self-gated, cine MRI acquisition scheme with data-driven cluster-based binning of cardiac motion.

METHODS

A Cartesian golden-step balanced steady-state free precession sequence with sorted k-space ordering was designed. Image data were acquired with breath-holding. Principal component analysis and k-means clustering were used for binning of cardiac phases. Cluster compactness in the time dimension was assessed using temporal variability, and dispersion in the spatial dimension was assessed using the Caliński-Harabasz index. The proposed and the reference electrocardiogram (ECG)-gated cine methods were compared using a four-point image quality score, SNR and CNR values, and Bland-Altman analyses of ventricular function.

RESULTS

A total of 10 subjects with sinus rhythm and 8 subjects with arrhythmias underwent cardiac MRI at 3.0 T. The temporal variability was 45.6 ms (cluster) versus 24.6 ms (ECG-based) (p < 0.001), and the Caliński-Harabasz index was 59.1 ± 9.1 (cluster) versus 22.0 ± 7.1 (ECG based) (p < 0.001). In subjects with sinus rhythm, 100% of the end-systolic and end-diastolic images from both the cluster and reference approach received the highest image quality score of 4. Relative to the reference cine images, the cluster-based multiphase (cine) image quality consistently received a one-point lower score (p < 0.05), whereas the SNR and CNR values were not significantly different (p = 0.20). In cases with arrhythmias, 97.9% of the end-systolic and end-diastolic images from the cluster approach received an image quality score of 3 or more. The mean bias values for biventricular ejection fraction and volumes derived from the cluster approach versus reference cine were negligible.

CONCLUSION

ECG-free cine cardiac MRI with data-driven clustering for binning of cardiac motion is feasible and enables quantification of cardiac function.

Mitigation of Fetal Irradiation Injury from Mid-Gestation Total Body Radiation with Mitochondrial-Targeted GS-Nitroxide JP4-039

Victims of a radiation terrorist event will include pregnant women and unborn fetuses. Mitochondrial dysfunction and oxidative stress are key pathogenic factors of fetal irradiation injury. The goal of this preclinical study is to investigate the efficacy of mitigating fetal irradiation injury by maternal administration of the mitochondrial-targeted gramicidin S (GS)- nitroxide radiation mitigator, JP4-039. Pregnant female C57BL/6NTac mice received 3 Gy total body ionizing irradiation (TBI) at mid-gestation embryonic day 13.5 (E13.5). Using novel time- and-motion-resolved 4D in utero magnetic resonance imaging (4D-uMRI), we found TBI caused extensive injury to the fetal brain that included cerebral hemorrhage, loss of cerebral tissue, and hydrocephalus with excessive accumulation of cerebrospinal fluid (CSF). Histopathology of the fetal mouse brain showed broken cerebral vessels and elevated apoptosis. Further use of novel 4D Oxy-wavelet MRI capable of probing in vivo mitochondrial function in intact brain revealed significant reduction of mitochondrial function in the fetal brain after 3Gy TBI. This was validated by ex vivo Oroboros mitochondrial respirometry. Maternal administration JP4-039 one day after TBI (E14.5), which can pass through the placental barrier, significantly reduced fetal brain radiation injury and improved fetal brain mitochondrial respiration. This also preserved cerebral brain tissue integrity and reduced cerebral hemorrhage and cell death. As JP4-039 administration did not change litter sizes or fetus viability, together these findings indicate JP4-039 can be deployed as a safe and effective mitigator of fetal radiation injury from mid-gestational in utero ionizing radiation exposure.

One Sentence Summary

Mitochondrial-targeted gramicidin S (GS)-nitroxide JP4-039 is safe and effective radiation mitigator for mid-gestational fetal irradiation injury.

The future of cardiovascular magnetic resonance: All-in-one vs. real-time (Part 1)

Cardiovascular magnetic resonance (CMR) protocols can be lengthy and complex, which has driven the research community to develop new technologies to make these protocols more efficient and patient-friendly. Two different approaches to improving CMR have been proposed, specifically "all-in-one" CMR, where several contrasts and/or motion states are acquired simultaneously, and "real-time" CMR, in which the examination is accelerated to avoid the need for breathholding and/or cardiac gating. The goal of this two-part manuscript is to describe these two different types of emerging rapid CMR. To this end, the vision of each is described, along with techniques which have been devised and tested along the pathway of clinical implementation. The pros and cons of the different methods are presented, and the remaining open needs of each are detailed. Part 1 will tackle the "all-in-one" approaches, and Part 2 the "real-time" approaches along with an overall summary of these emerging methods.

Rapid three-dimensional quantification of high-intensity plaques from coronary atherosclerosis T1-weighted characterization to predict periprocedural myocardial injury

BACKGROUND

High-intensity plaque (HIP) on magnetic resonance imaging (MRI) has been documented as a powerful predictor of periprocedural myocardial injury (PMI) following percutaneous coronary intervention (PCI). Despite the recent proposal of three-dimensional HIP quantification to enhance the predictive capability, the conventional pulse sequence, which necessitates the separate acquisition of anatomical reference images, hinders accurate three-dimensional segmentation along the coronary vasculature. Coronary atherosclerosis T1-weighted characterization (CATCH) enables the simultaneous acquisition of inherently coregistered dark-blood plaque and bright-blood coronary artery images. We aimed to develop a novel HIP quantification approach using CATCH and to ascertain its superior predictive performance compared to the conventional two-dimensional assessment based on plaque-to-myocardium signal intensity ratio (PMR).

METHODS

In this prospective study, CATCH MRI was conducted before elective stent implantation in 137 lesions from 125 patients. On CATCH images, dedicated software automatically generated tubular three-dimensional volumes of interest on the dark-blood plaque images along the coronary vasculature, based on the precisely matched bright-blood coronary artery images, and subsequently computed PMR and HIP volume (HIPvol). Specifically, HIPvol was calculated as the volume of voxels with signal intensity exceeding that of the myocardium, weighted by their respective signal intensities. PMI was defined as post-PCI cardiac troponin-T > 5 × the upper reference limit.

RESULTS

The entire analysis process was completed within 3 min per lesion. PMI occurred in 44 lesions. Based on the receiver operating characteristic curve analysis, HIPvol outperformed PMR for predicting PMI (C-statistics, 0.870 [95% CI, 0.805-0.936] vs. 0.787 [95% CI, 0.706-0.868]; p = 0.001). This result was primarily driven by the higher sensitivity HIPvol offered: 0.886 (95% CI, 0.754-0.962) vs. 0.750 for PMR (95% CI, 0.597-0.868; p = 0.034). Multivariable analysis identified HIPvol as an independent predictor of PMI (odds ratio, 1.15 per 10-μL increase; 95% CI, 1.01-1.30, p = 0.035).

CONCLUSIONS

Our semi-automated method of analyzing coronary plaque using CATCH MRI provided rapid HIP quantification. Three-dimensional assessment using this approach had a better ability to predict PMI than conventional two-dimensional assessment.

Physics-informed deep learning for T2-deblurred superresolution turbo spin echo MRI

PURPOSE

Deep learning superresolution (SR) is a promising approach to reduce MRI scan time without requiring custom sequences or iterative reconstruction. Previous deep learning SR approaches have generated low-resolution training images by simple k-space truncation, but this does not properly model in-plane turbo spin echo (TSE) MRI resolution degradation, which has variable T2 relaxation effects in different k-space regions. To fill this gap, we developed a T2 -deblurred deep learning SR method for the SR of 3D-TSE images.

METHODS

A SR generative adversarial network was trained using physically realistic resolution degradation (asymmetric T2 weighting of raw high-resolution k-space data). For comparison, we trained the same network structure on previous degradation models without TSE physics modeling. We tested all models for both retrospective and prospective SR with 3 × 3 acceleration factor (in the two phase-encoding directions) of genetically engineered mouse embryo model TSE-MR images.

RESULTS

The proposed method can produce high-quality 3 × 3 SR images for a typical 500-slice volume with 6-7 mouse embryos. Because 3 × 3 SR was performed, the image acquisition time can be reduced from 15 h to 1.7 h. Compared to previous SR methods without TSE modeling, the proposed method achieved the best quantitative imaging metrics for both retrospective and prospective evaluations and achieved the best imaging-quality expert scores for prospective evaluation.

CONCLUSION

The proposed T2 -deblurring method improved accuracy and image quality of deep learning-based SR of TSE MRI. This method has the potential to accelerate TSE image acquisition by a factor of up to 9.

Retrospective T2 quantification from conventional weighted MRI of the prostate based on deep learning

Purpose

To develop a deep learning-based method to retrospectively quantify T2 from conventional T1- and T2-weighted images.

Methods

Twenty-five subjects were imaged using a multi-echo spin-echo sequence to estimate reference prostate T2 maps. Conventional T1- and T2-weighted images were acquired as the input images. A U-Net based neural network was developed to directly estimate T2 maps from the weighted images using a four-fold cross-validation training strategy. The structural similarity index (SSIM), peak signal-to-noise ratio (PSNR), mean percentage error (MPE), and Pearson correlation coefficient were calculated to evaluate the quality of network-estimated T2 maps. To explore the potential of this approach in clinical practice, a retrospective T2 quantification was performed on a high-risk prostate cancer cohort (Group 1) and a low-risk active surveillance cohort (Group 2). Tumor and non-tumor T2 values were evaluated by an experienced radiologist based on region of interest (ROI) analysis.

Results

The T2 maps generated by the trained network were consistent with the corresponding reference. Prostate tissue structures and contrast were well preserved, with a PSNR of 26.41 ± 1.17 dB, an SSIM of 0.85 ± 0.02, and a Pearson correlation coefficient of 0.86. Quantitative ROI analyses performed on 38 prostate cancer patients revealed estimated T2 values of 80.4 ± 14.4 ms and 106.8 ± 16.3 ms for tumor and non-tumor regions, respectively. ROI measurements showed a significant difference between tumor and non-tumor regions of the estimated T2 maps (P < 0.001). In the two-timepoints active surveillance cohort, patients defined as progressors exhibited lower estimated T2 values of the tumor ROIs at the second time point compared to the first time point. Additionally, the T2 difference between two time points for progressors was significantly greater than that for non-progressors (P = 0.010).

Conclusion

A deep learning method was developed to estimate prostate T2 maps retrospectively from clinically acquired T1- and T2-weighted images, which has the potential to improve prostate cancer diagnosis and characterization without requiring extra scans.

Direct synthesis of multi-contrast brain MR images from MR multitasking spatial factors using deep learning

PURPOSE

To develop a deep learning method to synthesize conventional contrast-weighted images in the brain from MR multitasking spatial factors.

METHODS

Eighteen subjects were imaged using a whole-brain quantitative T1 -T2 -T1ρ MR multitasking sequence. Conventional contrast-weighted images consisting of T1 MPRAGE, T1 gradient echo, and T2 fluid-attenuated inversion recovery were acquired as target images. A 2D U-Net-based neural network was trained to synthesize conventional weighted images from MR multitasking spatial factors. Quantitative assessment and image quality rating by two radiologists were performed to evaluate the quality of deep-learning-based synthesis, in comparison with Bloch-equation-based synthesis from MR multitasking quantitative maps.

RESULTS

The deep-learning synthetic images showed comparable contrasts of brain tissues with the reference images from true acquisitions and were substantially better than the Bloch-equation-based synthesis results. Averaging on the three contrasts, the deep learning synthesis achieved normalized root mean square error = 0.184 ± 0.075, peak SNR = 28.14 ± 2.51, and structural-similarity index = 0.918 ± 0.034, which were significantly better than Bloch-equation-based synthesis (p < 0.05). Radiologists' rating results show that compared with true acquisitions, deep learning synthesis had no notable quality degradation and was better than Bloch-equation-based synthesis.

CONCLUSION

A deep learning technique was developed to synthesize conventional weighted images from MR multitasking spatial factors in the brain, enabling the simultaneous acquisition of multiparametric quantitative maps and clinical contrast-weighted images in a single scan.

Motion-compensated T1 mapping in cardiovascular magnetic resonance imaging: a technical review

T 1 mapping is becoming a staple magnetic resonance imaging method for diagnosing myocardial diseases such as ischemic cardiomyopathy, hypertrophic cardiomyopathy, myocarditis, and more. Clinically, most T 1 mapping sequences acquire a single slice at a single cardiac phase across a 10 to 15-heartbeat breath-hold, with one to three slices acquired in total. This leaves opportunities for improving patient comfort and information density by acquiring data across multiple cardiac phases in free-running acquisitions and across multiple respiratory phases in free-breathing acquisitions. Scanning in the presence of cardiac and respiratory motion requires more complex motion characterization and compensation. Most clinical mapping sequences use 2D single-slice acquisitions; however newer techniques allow for motion-compensated reconstructions in three dimensions and beyond. To further address confounding factors and improve measurement accuracy, T 1 maps can be acquired jointly with other quantitative parameters such as T 2 , T 2 ∗ , fat fraction, and more. These multiparametric acquisitions allow for constrained reconstruction approaches that isolate contributions to T 1 from other motion and relaxation mechanisms. In this review, we examine the state of the literature in motion-corrected and motion-resolved T 1 mapping, with potential future directions for further technical development and clinical translation.

Retrospective quantification of clinical abdominal DCE-MRI using pharmacokinetics-informed deep learning: a proof-of-concept study

Introduction

Dynamic contrast-enhanced (DCE) MRI has important clinical value for early detection, accurate staging, and therapeutic monitoring of cancers. However, conventional multi-phasic abdominal DCE-MRI has limited temporal resolution and provides qualitative or semi-quantitative assessments of tissue vascularity. In this study, the feasibility of retrospectively quantifying multi-phasic abdominal DCE-MRI by using pharmacokinetics-informed deep learning to improve temporal resolution was investigated.

Method

Forty-five subjects consisting of healthy controls, pancreatic ductal adenocarcinoma (PDAC), and chronic pancreatitis (CP) were imaged with a 2-s temporal-resolution quantitative DCE sequence, from which 30-s temporal-resolution multi-phasic DCE-MRI was synthesized based on clinical protocol. A pharmacokinetics-informed neural network was trained to improve the temporal resolution of the multi-phasic DCE before the quantification of pharmacokinetic parameters. Through ten-fold cross-validation, the agreement between pharmacokinetic parameters estimated from synthesized multi-phasic DCE after deep learning inference was assessed against reference parameters from the corresponding quantitative DCE-MRI images. The ability of the deep learning estimated parameters to differentiate abnormal from normal tissues was assessed as well.

Results

The pharmacokinetic parameters estimated after deep learning have a high level of agreement with the reference values. In the cross-validation, all three pharmacokinetic parameters (transfer constant K trans , fractional extravascular extracellular volume v e , and rate constant k ep ) achieved intraclass correlation coefficient and R2 between 0.84-0.94, and low coefficients of variation (10.1%, 12.3%, and 5.6%, respectively) relative to the reference values. Significant differences were found between healthy pancreas, PDAC tumor and non-tumor, and CP pancreas.

Discussion

Retrospective quantification (RoQ) of clinical multi-phasic DCE-MRI is possible by deep learning. This technique has the potential to derive quantitative pharmacokinetic parameters from clinical multi-phasic DCE data for a more objective and precise assessment of cancer.

MR Multitasking-based multi-dimensional assessment of cardiovascular system (MT-MACS) with extended spatial coverage and water-fat separation

PURPOSE

To extend the MR MultiTasking-based Multidimensional Assessment of Cardiovascular System (MT-MACS) technique with larger spatial coverage and water-fat separation for comprehensive aortocardiac assessment.

METHODS

MT-MACS adopts a low-rank tensor image model for 7D imaging, with three spatial dimensions for volumetric imaging, one cardiac motion dimension for cine imaging, one respiratory motion dimension for free-breathing imaging, one T2-prepared inversion recovery time dimension for multi-contrast assessment, and one T2*-decay time dimension for water-fat separation. Nine healthy subjects were recruited for the 3T study. Overall image quality was scored on bright-blood (BB), dark-blood (DB), and gray-blood (GB) contrasts using a 4-point scale (0-poor to 3-excellent) by two independent readers, and their interreader agreement was evaluated. Myocardial wall thickness and left ventricular ejection fraction (LVEF) were quantified on DB and BB contrasts, respectively. The agreement in these metrics between MT-MACS and conventional breath-held, electrocardiography-triggered 2D sequences were evaluated.

RESULTS

MT-MACS provides both water-only and fat-only images with excellent image quality (average score = 3.725/3.780/3.835/3.890 for BB/DB/GB/fat-only images) and moderate to high interreader agreement (weighted Cohen's kappa value = 0.727/0.668/1.000/1.000 for BB/DB/GB/fat-only images). There were good to excellent agreements in myocardial wall thickness measurements (intraclass correlation coefficients [ICC] = 0.781/0.929/0.680/0.878 for left atria/left ventricle/right atria/right ventricle) and LVEF quantification (ICC = 0.716) between MT-MACS and 2D references. All measurements were within the literature range of healthy subjects.

CONCLUSION

The refined MT-MACS technique provides multi-contrast, phase-resolved, and water-fat imaging of the aortocardiac systems and allows evaluation of anatomy and function. Clinical validation is warranted.

Free-breathing 3D CEST MRI of human liver at 3.0 T

PURPOSE

To develop a novel 3D abdominal CEST MRI technique at 3 T using MR multitasking, which enables entire-liver coverage with free-breathing acquisition.

METHODS

k-Space data were continuously acquired with repetitive steady-state CEST (ss-CEST) modules. The stack-of-stars acquisition pattern was used for k-space sampling. MR multitasking was used to reconstruct motion-resolved 3D CEST images of 53 frequency offsets with entire-liver coverage and 2.0 × 2.0 × 6.0 mm3 spatial resolution. The total scan time was 9 min. The sensitivity of amide proton transfer (APT)-CEST (magnetization transfer asymmetry [MTRasym ] at 3.5 ppm) and glycogen CEST (glycoCEST) (mean MTRasym around 1.0 ppm) signals generated with the proposed method were tested with fasting experiments.

RESULTS

Both APT-CEST and glycoCEST signals showed high sensitivity between post-fasting and post-meal acquisitions. APT-CEST and glycoCEST MTRasym signals from post-mean scans were significantly increased (APT-CEST: -0.019 ± 0.017 in post-fasting scans, 0.014 ± 0.021 in post-meal scans, p < 0.01; glycoCEST: 0.003 ± 0.009 in post-fasting scans, 0.027 ± 0.021 in post-meal scans, p < 0.01).

CONCLUSION

The proposed 3D abdominal steady-state CEST method using MR multitasking can generate CEST images of the entire liver during free breathing.

Multitasking dynamic contrast enhanced magnetic resonance imaging can accurately differentiate chronic pancreatitis from pancreatic ductal adenocarcinoma

Background and aims

Accurate differentiation of chronic pancreatitis (CP) and pancreatic ductal adenocarcinoma (PDAC) is an area of unmet clinical need. In this study, a novel Multitasking dynamic contrast enhanced (DCE) magnetic resonance imaging (MRI) technique was used to quantitatively evaluate the microcirculation properties of pancreas in CP and PDAC and differentiate between them.

Methods

The Multitasking DCE technique was able to acquire one 3D image per second during the passage of MRI contrast agent, allowing the quantitative estimation of microcirculation properties of tissue, including blood flow F_p, plasma volume fraction v_p, transfer constant K^trans, and extravascular extracellular volume fraction v_e. Receiver operating characteristic (ROC) analysis was performed to differentiate the CP pancreas, PDAC pancreas, normal control pancreas, PDAC tumor, PDAC upstream, and PDAC downstream. ROCs from quantitative analysis and conventional analysis were compared.

Results

Fourteen PDAC patients, 8 CP patients and 20 healthy subjects were prospectively recruited. The combination of F_p, v_p, K^trans, and v_e can differentiate CP versus PDAC pancreas with good AUC (AUC [95% CI] = 0.821 [0.654 - 0.988]), CP versus normal pancreas with excellent AUC (1.000 [1.000 - 1.000]), PDAC pancreas versus normal pancreas with excellent AUC (1.000 [1.000 - 1.000]), CP versus PDAC tumor with excellent AUC (1.000 [1.000 - 1.000]), CP versus PDAC downstream with excellent AUC (0.917 [0.795 - 1.000]), and CP versus PDAC upstream with fair AUC (0.722 [0.465 - 0.980]). This quantitative analysis outperformed conventional analysis in differentiation of each pair.

Conclusion

Multitasking DCE MRI is a promising clinical tool that is capable of unbiased quantitative differentiation between CP from PDAC.

MR multitasking-based dynamic imaging for cerebrovascular evaluation (MT-DICE): Simultaneous quantification of permeability and leakage-insensitive perfusion by dynamic T 1 / T 2 * mapping

PURPOSE

To develop an MR multitasking-based dynamic imaging for cerebrovascular evaluation (MT-DICE) technique for simultaneous quantification of permeability and leakage-insensitive perfusion with a single-dose contrast injection.

METHODS

MT-DICE builds on a saturation-recovery prepared multi-echo fast low-angle shot sequence. The k-space is randomly sampled for 7.6 min, with single-dose contrast agent injected 1.5 min into the scan. MR multitasking is used to model the data into six dimensions, including three spatial dimensions for whole-brain coverage, a saturation-recovery time dimension, and a TE dimension for dynamicT 1 $$ {\mathrm{T}}_1 $$ andT 2 * $$ {\mathrm{T}}_2^{\ast } $$ quantification, respectively, and a contrast dynamics dimension for capturing contrast kinetics. The derived pixel-wiseT 1 / T 2 * $$ {\mathrm{T}}_1/{\mathrm{T}}_2^{\ast } $$ time series are converted into contrast concentration-time curves for calculation of kinetic metrics. The technique was assessed for its agreement with reference methods inT 1 $$ {\mathrm{T}}_1 $$ andT 2 * $$ {\mathrm{T}}_2^{\ast } $$ measurements in eight healthy subjects and, in three of them, inter-session repeatability of permeability and leakage-insensitive perfusion parameters. Its feasibility was also demonstrated in four patients with brain tumors.

RESULTS

MT-DICET 1 / T 2 * $$ {\mathrm{T}}_1/{\mathrm{T}}_2^{\ast } $$ values of normal gray matter and white matter were in excellent agreement with reference values (intraclass correlation coefficients = 0.860/0.962 for gray matter and 0.925/0.975 for white matter ). Both permeability and perfusion parameters demonstrated good to excellent intersession agreement with the lowest intraclass correlation coefficients at 0.694. Contrast kinetic parameters in all healthy subjects and patients were within the literature range.

CONCLUSION

Based on dynamicT 1 / T 2 * $$ {\mathrm{T}}_1/{\mathrm{T}}_2^{\ast } $$ mapping, MT-DICE allows for simultaneous quantification of permeability and leakage-insensitive perfusion metrics with a single-dose contrast injection.

Free-breathing, non-ECG, simultaneous myocardial T1 , T2 , T2 *, and fat-fraction mapping with motion-resolved cardiovascular MR multitasking

PURPOSE

To develop a free-breathing, non-electrocardiogram technique for simultaneous myocardial T1 , T2 , T2 *, and fat-fraction (FF) mapping in a single scan.

METHODS

The MR Multitasking framework is adapted to quantify T1 , T2 , T2 *, and FF simultaneously. A variable TR scheme is developed to preserve temporal resolution and imaging efficiency. The underlying high-dimensional image is modeled as a low-rank tensor, which allows accelerated acquisition and efficient reconstruction. The accuracy and/or repeatability of the technique were evaluated on static and motion phantoms, 12 healthy volunteers, and 3 patients by comparing to the reference techniques.

RESULTS

In static and motion phantoms, T1 /T2 /T2 */FF measurements showed substantial consistency (R > 0.98) and excellent agreement (intraclass correlation coefficient > 0.93) with reference measurements. In human subjects, the proposed technique yielded repeatable T1 , T2 , T2 *, and FF measurements that agreed with those from references.

CONCLUSIONS

The proposed free-breathing, non-electrocardiogram, motion-resolved Multitasking technique allows simultaneous quantification of myocardial T1 , T2 , T2 *, and FF in a single 2.5-min scan.

OSCILLATE: A low-rank approach for accelerated magnetic resonance elastography

PURPOSE

MR elastography (MRE) is a technique to characterize brain mechanical properties in vivo. Due to the need to capture tissue deformation in multiple directions over time, MRE is an inherently long acquisition, which limits achievable resolution and use in challenging populations. The purpose of this work is to develop a method for accelerating MRE acquisition by using low-rank image reconstruction to exploit inherent spatiotemporal correlations in MRE data.

THEORY AND METHODS

The proposed MRE sampling and reconstruction method, OSCILLATE (Observing Spatiotemporal Correlations for Imaging with Low-rank Leveraged Acceleration in Turbo Elastography), involves alternating which k-space points are sampled between each repetition by a reduction factor, ROSC. Using a predetermined temporal basis from a low-resolution navigator in a joint low-rank image reconstruction, all images can be accurately reconstructed from a reduced amount of k-space data.

RESULTS

Decomposition of MRE displacement data demonstrated that, on average, 96.1% of all energy from an MRE dataset is captured at rank L = 12 (reduced from a full rank of 24). Retrospectively undersampling data with ROSC  = 2 and reconstructing at low-rank (L = 12) yields highly accurate stiffness maps with voxel-wise error of 5.8% ± 0.7%. Prospectively undersampled data at ROSC  = 2 were successfully reconstructed without loss of material property map fidelity, with average global stiffness error of 1.0% ± 0.7% compared to fully sampled data.

CONCLUSIONS

OSCILLATE produces whole-brain MRE data at 2 mm isotropic resolution in 1 min 48 s.

Single projection driven real-time multi-contrast (SPIDERM) MR imaging using pre-learned spatial subspace and linear transformation

Objective. To develop and test the feasibility of a novel Single ProjectIon DrivEn Real-time Multi-contrast (SPIDERM) MR imaging technique that can generate real-time 3D images on-the-fly with flexible contrast weightings and a low latency. Approach. In SPIDERM, a ‘prep’ scan is first performed, with sparse k-space sampling periodically interleaved with the central k-space line (navigator data), to learn a subject-specific model, incorporating a spatial subspace and a linear transformation between navigator data and subspace coordinates. A ‘live’ scan is then performed by repeatedly acquiring the central k-space line only to dynamically determine subspace coordinates. With the ‘prep’-learned subspace and ‘live’ coordinates, real-time 3D images are generated on-the-fly with computationally efficient matrix multiplication. When implemented based on a multi-contrast pulse sequence, SPIDERM further allows for data-driven image contrast regeneration to convert real-time contrast-varying images into contrast-frozen images at user’s discretion while maintaining motion states. Both digital phantom and in-vivo experiments were performed to evaluate the technical feasibility of SPIDERM. Main results. The elapsed time from the input of the central k-space line to the generation of real-time contrast-frozen 3D images was approximately 45 ms, permitting a latency of 55 ms or less. Motion displacement measured from SPIDERM and reference images showed excellent correlation ( R 2 ≥ 0.983 ). Geometric variation from the ground truth in the digital phantom was acceptable as demonstrated by pancreas contour analysis (Dice ≥ 0.84, mean surface distance ≤ 0.95 mm). Quantitative image quality metrics showed good consistency between reference images and contrast-varying SPIDREM images in in-vivo studies (mean NMRSE = 0.141 , PSNR = 3 0.12 , SSIM = 0.88 ). Significance. SPIDERM is capable of generating real-time multi-contrast 3D images with a low latency. An imaging framework based on SPIDERM has the potential to serve as a standalone package for MR-guided radiation therapy by offering adaptive simulation through a ‘prep’ scan and real-time image guidance through a ‘live’ scan.

Simultaneous Multi-Slice Cardiac MR Multitasking for Motion-Resolved, Non-ECG, Free-Breathing T1–T2 Mapping

The aim of this study is to simultaneously quantify T1/T2 across three slices of the left-ventricular myocardium without breath-holds or ECG monitoring, all within a 3 min scan. Radial simultaneous multi-slice (SMS) encoding, self-gating, and image reconstruction was incorporated into the cardiovascular magnetic resonance (CMR) Multitasking framework to simultaneously image three short-axis slices. A T2prep-IR FLASH sequence with two flip angles was designed and implemented to allow B1+-robust T1 and T2 mapping. The proposed Multitasking-SMS method was validated in a standardized phantom and 10 healthy volunteers, comparing T1 and T2 measurements and scan-rescan repeatability against corresponding reference methods in one layer of phantom vials and in 16 American Heart Association (AHA) myocardial segments. In phantom, Multitasking-SMS T1/T2 measurements showed substantial correlation (R² > 0.996) and excellent agreement [intraclass correlation coefficients (ICC) ≥ 0.999)] with reference measurements. In healthy volunteers, Multitasking-SMS T1/T2 maps reported similar myocardial T1/T2 values (1,215 ± 91.0/41.5 ± 6.3 ms) to the reference myocardial T1/T2 values (1,239 ± 67.5/42.7 ± 4.1 ms), with P = 0.347 and P = 0.296, respectively. Bland–Altman analyses also demonstrated good in vivo repeatability in both the multitasking and references, with segment-wise coefficients of variation of 4.7% (multitasking T1), 8.9% (multitasking T2), 2.4% [modified look-locker inversion recovery (MOLLI)], and 4.6% (T2-prep FLASH), respectively. In summary, multitasking-SMS is feasible for free-breathing, non-ECG, myocardial T1/T2 quantification in 16 AHA segments over 3 short-axis slices in 3 min. The method shows the great potential for reducing exam time for quantitative CMR without ECG or breath-holds.

Data-Consistent non-Cartesian deep subspace learning for efficient dynamic MR image reconstruction

Non-Cartesian sampling with subspace-constrained image reconstruction is a popular approach to dynamic MRI, but slow iterative reconstruction limits its clinical application. Data-consistent (DC) deep learning can accelerate reconstruction with good image quality, but has not been formulated for non-Cartesian subspace imaging. In this study, we propose a DC non-Cartesian deep subspace learning framework for fast, accurate dynamic MR image reconstruction. Four novel DC formulations are developed and evaluated: two gradient decent approaches, a directly solved approach, and a conjugate gradient approach. We applied a U-Net model with and without DC layers to reconstruct T1-weighted images for cardiac MR Multitasking (an advanced multidimensional imaging method), comparing our results to the iteratively reconstructed reference. Experimental results show that the proposed framework significantly improves reconstruction accuracy over the U-Net model without DC, while significantly accelerating the reconstruction over conventional iterative reconstruction.

Three-dimensional simultaneous brain mapping of T1, T2, T2∗ and magnetic susceptibility with MR Multitasking

PURPOSE

To develop a new technique that enables simultaneous quantification of whole-brain T1 , T2 , T 2 ∗ , as well as susceptibility and synthesis of six contrast-weighted images in a single 9.1-minute scan.

METHODS

The technique uses hybrid T2 -prepared inversion-recovery pulse modules and multi-echo gradient-echo readouts to collect k-space data with various T1, T2, and T 2 ∗ weightings. The underlying image is represented as a six-dimensional low-rank tensor consisting of three spatial dimensions and three temporal dimensions corresponding to T1 recovery, T2 decay, and multi-echo behaviors, respectively. Multiparametric maps were fitted from reconstructed image series. The proposed method was validated on phantoms and healthy volunteers, by comparing quantitative measurements against corresponding reference methods. The feasibility of generating six contrast-weighted images was also examined.

RESULTS

High quality, co-registered T1 , T2 , and T 2 ∗ susceptibility maps were generated that closely resembled the reference maps. Phantom measurements showed substantial consistency (R2 > 0.98) with the reference measurements. Despite the significant differences of T1 (p < .001), T2 (p = .002), and T 2 ∗ (p = 0.008) between our method and the references for in vivo studies, excellent agreement was achieved with all intraclass correlation coefficients greater than 0.75. No significant difference was found for susceptibility (p = .900). The framework is also capable of synthesizing six contrast-weighted images.

CONCLUSION

The MR Multitasking-based 3D brain mapping of T1 , T2 , T 2 ∗ , and susceptibility agrees well with the reference and is a promising technique for multicontrast and quantitative imaging.

Motion-robust quantitative multiparametric brain MRI with motion-resolved MR multitasking

PURPOSE

To address head motion in brain MRI with a novel motion-resolved imaging framework, with application to motion-robust quantitative multiparametric mapping.

METHODS

MR multitasking conceptualizes the underlying multiparametric image in the presence of motion as a multidimensional low-rank tensor. By incorporating a motion-state dimension into the parameter dimensions and introducing unsupervised motion-state binning and outlier motion reweighting mechanisms, the brain motion can be readily resolved for motion-robust quantitative imaging. A numerical motion phantom was used to simulate different discrete and periodic motion patterns under various translational and rotational scenarios, as well as investigate the sensitivity to exceptionally small and large displacements. In vivo brain MRI was performed to also evaluate different real discrete and periodic motion patterns. The effectiveness of motion-resolved imaging for simultaneous T₁ /T₂ /T_1ρ mapping in the brain was investigated.

RESULTS

For all 14 simulation scenarios of small, intermediate, and large motion displacements, the motion-resolved approach produced T₁ /T₂ /T_1ρ maps with less absolute difference errors against the ground truth, lower RMSE, and higher structural similarity index measure of T₁ /T₂ /T1_ρ measurements compared to motion removal, and no motion handling. For in vivo experiments, the motion-resolved approach produced T₁ /T₂ /T_1ρ maps with the best image quality free from motion artifacts under random discrete motion, tremor, periodic shaking, and nodding patterns compared to motion removal and no motion handling. The proposed method also yielded T₁ /T₂ /T_1ρ measurement distributions closest to the motion-free reference, with minimal measurement bias and variance.

CONCLUSION

Motion-resolved quantitative brain imaging is achieved with multitasking, which is generalizable to various head motion patterns without explicit need for registration-based motion correction.

Multiparametric mapping in the brain from conventional contrast-weighted images using deep learning

PURPOSE

To develop a deep-learning-based method to quantify multiple parameters in the brain from conventional contrast-weighted images.

METHODS

Eighteen subjects were imaged using an MR Multitasking sequence to generate reference T₁ and T₂ maps in the brain. Conventional contrast-weighted images consisting of T₁ MPRAGE, T₁ GRE, and T₂ FLAIR were acquired as input images. A U-Net-based neural network was trained to estimate T₁ and T₂ maps simultaneously from the contrast-weighted images. Six-fold cross-validation was performed to compare the network outputs with the MR Multitasking references.

RESULTS

The deep-learning T₁ /T₂ maps were comparable with the references, and brain tissue structures and image contrasts were well preserved. A peak signal-to-noise ratio >32 dB and a structural similarity index >0.97 were achieved for both parameter maps. Calculated on brain parenchyma (excluding CSF), the mean absolute errors (and mean percentage errors) for T₁ and T₂ maps were 52.7 ms (5.1%) and 5.4 ms (7.1%), respectively. ROI measurements on four tissue compartments (cortical gray matter, white matter, putamen, and thalamus) showed that T₁ and T₂ values provided by the network outputs were in agreement with the MR Multitasking reference maps. The mean differences were smaller than ± 1%, and limits of agreement were within ± 5% for T₁ and within ± 10% for T₂ after taking the mean differences into account.

CONCLUSION

A deep-learning-based technique was developed to estimate T₁ and T₂ maps from conventional contrast-weighted images in the brain, enabling simultaneous qualitative and quantitative MRI without modifying clinical protocols.

Free-breathing multitasking multi-echo MRI for whole-liver water-specific T1 , proton density fat fraction, and R2∗ quantification

PURPOSE

To develop a 3D multitasking multi-echo (MT-ME) technique for the comprehensive characterization of liver tissues with 5-min free-breathing acquisition; whole-liver coverage; a spatial resolution of 1.5 × 1.5 × 6 mm³ ; and simultaneous quantification of T₁ , water-specific T₁ (T_1w ), proton density fat fraction (PDFF), and R2∗ .

METHODS

Six-echo bipolar spoiled gradient echo readouts following inversion recovery preparation was performed to generate T₁ , water/fat, and R2∗ contrast. MR multitasking was used to reconstruct the MT-ME images with 3 spatial dimensions: 1 T₁ recovery dimension, 1 multi-echo dimension, and 1 respiratory dimension. A basis function-based approach was developed for T_1w quantification, followed by the estimation of R2∗ and T₁ -corrected PDFF. The intrasession repeatability and agreement against references of MT-ME measurements were tested on a phantom and 15 clinically healthy subjects. In addition, 4 patients with confirmed liver diseases were recruited, and the agreement between MT-ME measurements and references was assessed.

RESULTS

MT-ME produced high-quality, coregistered T₁ , T_1w , PDFF, and R2∗ maps with good intrasession repeatability and substantial agreement with references on phantom and human studies. The intra-class coefficients of T₁ , T_1w , PDFF, and R2∗ from the repeat MT-ME measurements on clinically healthy subjects were 0.989, 0.990, 0.999, and 0.988, respectively. The intra-class coefficients of T₁ , PDFF, and R2∗ between the MT-ME and reference measurements were 0.924, 0.987, and 0.975 in healthy subjects and 0.980, 0.999, and 0.998 in patients. The T_1w was independent to PDFF (R = -0.029, P = .904).

CONCLUSION

The proposed MT-ME technique quantifies T₁ , T_1w , PDFF, and R2∗ simultaneously and is clinically promising for the comprehensive characterization of liver tissue properties.

Coronary High-Intensity Plaques at T1-weighted MRI in Stable Coronary Artery Disease: Comparison with Near-Infrared Spectroscopy Intravascular US

Background The histologic nature of coronary high-intensity plaques (HIPs) at T1-weighted MRI in patients with stable coronary artery disease remains to be fully understood. Coronary atherosclerosis T1-weighted characterization (CATCH) enables HIP detection by simultaneously acquiring dark-blood plaque and bright-blood anatomic reference images. Purpose To determine if intraplaque hemorrhage (IPH) or lipid is the predominant substrate of HIPs on T1-weighted images by comparing CATCH MRI scans with findings on near-infrared spectroscopy (NIRS) intravascular US (IVUS) images. Materials and Methods This study retrospectively included consecutive patients who underwent CATCH MRI before NIRS IVUS between December 2019 and February 2021 at two facilities. At MRI, HIP was defined as plaque-to-myocardium signal intensity ratio of at least 1.4. The presence of an echolucent zone at IVUS (reported to represent IPH) was recorded. NIRS was used to determine the lipid component of atherosclerotic plaque. Lipid core burden index (LCBI) was calculated as the fraction of pixels with a probability of lipid-core plaque greater than 0.6 within a region of interest. Plaque with maximum LCBI within any 4-mm-long segment (maxLCBI_{4 mm}) greater than 400 was regarded as lipid rich. Multivariable analysis was performed to evaluate NIRS IVUS-derived parameters associated with HIPs. Results There were 205 plaques analyzed in 95 patients (median age, 74 years; interquartile range [IQR], 67-78 years; 75 men). HIPs (n = 42) at MRI were predominantly associated with an echolucent zone at IVUS (79% [33 of 42] vs 8.0% [13 of 163], respectively; P < .001) and a higher maxLCBI_{4 mm} at NIRS (477 [IQR, 258-738] vs 232 [IQR, 59-422], respectively; P < .001) than non-HIPs. In the multivariable model, HIPs were independently associated with an echolucent zone (odds ratio, 24.5; 95% CI: 9.3, 64.7; P < .001), but not with lipid-rich plaque (odds ratio, 2.0; 95% CI: 0.7, 5.4; P = .20). Conclusion The predominant substrate of T1-weighed MRI-defined high-intensity plaques in stable coronary artery disease was intraplaque hemorrhage, not lipid. © RSNA, 2021 Online supplemental material is available for this article. See also the editorial by Stuber in this issue.

Dual flip-angle IR-FLASH with spin history mapping for B1+ corrected T1 mapping: Application to T1 cardiovascular magnetic resonance multitasking

PURPOSE

To develop a single-scan method for B 1 + -corrected T1 mapping and apply it for free-breathing (FB) cardiac MR multitasking without electrocardiogram (ECG) triggering.

METHODS

One dual flip-angle (2FA) inversion recovery (IR)-FLASH scan provides two observations of T 1 ∗ (apparent T1 ) corresponding to two distinct combinations of the nominal FA α and B 1 + . Spatiotemporally coregistered T1 and B 1 + spin history maps are obtained by fitting the 2FA signal model. T1 estimate accuracy and repeatability for single flip-angle (1FA) and 2FA IR-FLASH sequence MR multitasking were evaluated at 3T. A T1 phantom was first imaged on the scanner table, then on two human subjects' thoraxes in both breath-hold (BH) and FB conditions. IR-turbo spin echo (IR-TSE) static phantom T1 measurements served as reference. In 10 healthy subjects, myocardial T1 was evaluated with ECG-free, FB multitasking sequences alongside ECG-triggered BH MOLLI.

RESULTS

For phantom-on-table T1 estimates, 2FA agreed better with IR-TSE (intraclass correlation coefficient [ICC] = 0.996, mean error ± SD = -1.6% ± 1.9%) than did 1FA (ICC = 0.922; mean error ± SD = -4.3% ± 12%). For phantom-on-thorax, 2FA was more repeatable and robust to respiration than 1FA (coefficient of variation [CoV] = 1.2% 2FA, = 11.3% 1FA). In vivo, in intrasession T1 repeatability, 2FA (septal CoV = 2.4%, six-segment CoV = 4.4%) outperformed 1FA (septal CoV = 3.1%, six-segment CoV = 5.5%). In six-segment T1 homogeneity, 2FA (CoV = 7.9%) also outperformed 1FA (CoV = 11.1%).

CONCLUSION

The 2FA IR-FLASH improves T1 estimate accuracy and repeatability over 1FA IR-FLASH, and enables single-scan B 1 + -corrected T1 mapping without BHs or ECG when used with MR multitasking.

Whole-brain steady-state CEST at 3 T using MR Multitasking

PURPOSE

To perform fast 3D steady-state CEST (ss-CEST) imaging using MR Multitasking.

METHODS

A continuous acquisition sequence with repetitive ss-CEST modules was developed. Each ss-CEST module contains a single-lobe Gaussian saturation pulse, followed by a spoiler gradient and eight FLASH readouts (one "training line" + seven "imaging lines"). Three-dimensional Cartesian encoding was used for k-space acquisition. Reconstructed CEST images were quantified with four-pool Lorentzian fitting.

RESULTS

Steady-state CEST with whole-brain coverage was performed in 5.6 s per saturation frequency offset at the spatial resolution of 1.7 × 1.7 × 3.0 mm³ . The total scan time was 5.5 min for 55 different frequency offsets. Quantitative CEST maps from multipool fitting showed consistent image quality across the volume.

CONCLUSION

Three-dimensional ss-CEST with whole-brain coverage can be done at 3 T within 5.5 min using MR Multitasking.

Publisher Correction to: Electrocardiogram-less, free-breathing myocardial extracellular volume fraction mapping in small animals at high heart rates using motion-resolved cardiovascular magnetic resonance multitasking: a feasibility study in a heart failure with preserved ejection fraction rat model

An amendment to this paper has been published and can be accessed via the original article.

Electrocardiogram-less, free-breathing myocardial extracellular volume fraction mapping in small animals at high heart rates using motion-resolved cardiovascular magnetic reesonance multitasking: a feasibility study in a heart failure with preserved ejection fraction rat model

BACKGROUND

Extracellular volume fraction (ECV) quantification with cardiovascular magnetic resonance (CMR) T1 mapping is a powerful tool for the characterization of focal or diffuse myocardial fibrosis. However, it is technically challenging to acquire high-quality T1 and ECV maps in small animals for preclinical research because of high heart rates and high respiration rates. In this work, we developed an electrocardiogram (ECG)-less, free-breathing ECV mapping method using motion-resolved CMR Multitasking on a 9.4 T small animal CMR system. The feasibility of characterizing diffuse myocardial fibrosis was tested in a rat heart failure model with preserved ejection fraction (HFpEF).

METHODS

High-salt fed rats diagnosed with HFpEF (n = 9) and control rats (n = 9) were imaged with the proposed ECV Multitasking technique. A 25-min exam, including two 4-min T1 Multitasking scans before and after gadolinium injection, were performed on each rat. It allows a cardiac temporal resolution of 20 ms for a heart rate of ~ 300 bpm. Myocardial ECV was calculated from the hematocrit (HCT) and fitted T1 values of the myocardium and the blood pool. Masson's trichrome stain was used to measure the extent of fibrosis. Welch's t-test was performed between control and HFpEF groups.

RESULTS

ECV was significantly higher in the HFpEF group (22.4% ± 2.5% vs. 18.0% ± 2.1%, P = 0.0010). A moderate correlation between the ECV and the extent of fibrosis was found (R = 0.59, P = 0.0098).

CONCLUSIONS

Motion-resolved ECV Multitasking CMR can quantify ECV in the rat myocardium at high heart rates without ECG triggering or respiratory gating. Elevated ECV found in the HFpEF group is consistent with previous human studies and well correlated with histological data. This technique has the potential to be a viable imaging tool for myocardial tissue characterization in small animal models.

Five-dimensional quantitative low-dose Multitasking dynamic contrast- enhanced MRI: Preliminary study on breast cancer

PURPOSE

To develop a low-dose Multitasking DCE technique (LD-MT-DCE) for breast imaging, enabling dynamic T₁ mapping-based quantitative characterization of tumor blood flow and vascular properties with whole-breast coverage, a spatial resolution of 0.9 × 0.9 × 1.1 mm³ , and a temporal resolution of 1.4 seconds using a 20% gadolinium dose (0.02 mmol/kg).

METHODS

Magnetic resonance Multitasking was used to reconstruct 5D images with three spatial dimensions, one T₁ recovery dimension for dynamic T₁ quantification, and one DCE dimension for contrast kinetics. Kinetic parameters F p , v p , K trans , and v e were estimated from dynamic T₁ maps using the two-compartment exchange model. The LD-MT-DCE repeatability and agreement against standard-dose MT-DCE were evaluated in 20 healthy subjects. In 7 patients with triple-negative breast cancer, LD-MT-DCE image quality and diagnostic results were compared with that of standard-dose clinical DCE in the same imaging session. One-way unbalanced analysis of variance with Tukey test was performed to evaluate the statistical significance of the kinetic parameters between control and patient groups.

RESULTS

The LD-MT-DCE technique was repeatable, agreed with standard-dose MT-DCE, and showed excellent image quality. The diagnosis using LD-MT-DCE matched well with clinical results. The values of F p , v p , and K trans were significantly different between malignant tumors and normal breast tissue (P < .001, < .001, and < .001, respectively), and between malignant and benign tumors (P = .020, .003, and < .001, respectively).

CONCLUSION

The LD-MT-DCE technique was repeatable and showed excellent image quality and equivalent diagnosis compared with standard-dose clinical DCE. The estimated kinetic parameters were capable of differentiating between normal breast tissue and benign and malignant tumors.

Free-breathing diffusion tensor MRI of the whole left ventricle using second-order motion compensation and multitasking respiratory motion correction

PURPOSE

We aimed to develop a novel free-breathing cardiac diffusion tensor MRI (DT-MRI) approach, M2-MT-MOCO, capable of whole left ventricular coverage that leverages second-order motion compensation (M2) diffusion encoding and multitasking (MT) framework to efficiently correct for respiratory motion (MOCO).

METHODS

Imaging was performed in 16 healthy volunteers and 3 heart failure patients with symptomatic dyspnea. The healthy volunteers were scanned to compare the accuracy of interleaved multislice coverage of the entire left ventricle with a single-slice acquisition and the accuracy of the free-breathing conventional MOCO and MT-MOCO approaches with reference breath-hold DT-MRI. Mean diffusivity (MD), fractional anisotropy (FA), helix angle transmurality (HAT), and intrascan repeatability were quantified and compared.

RESULTS

In all subjects, free-breathing M2-MT-MOCO DT-MRI yielded DWI of the entire left ventricle without bulk motion-induced signal loss. No significant differences were seen in the global values of MD, FA, and HAT in the multislice and single-slice acquisitions. Furthermore, global quantification of MD, FA, and HAT were also not significantly different between the MT-MOCO and breath-hold, whereas conventional MOCO yielded significant differences in MD, FA, and HAT with MT-MOCO and FA with breath-hold. In heart failure patients, M2-MT-MOCO DT-MRI was feasible yielding higher MD, lower FA, and lower HAT compared with healthy volunteers. Substantial agreement was found between repeated scans across all subjects for MT-MOCO.

CONCLUSION

M2-MT-MOCO enables free-breathing DT-MRI of the entire left ventricle in 10 min, while preserving quantification of myocardial microstructure compared to breath-held and single-slice acquisitions and is feasible in heart failure patients.

Three-dimensional whole-brain simultaneous T1, T2, and T1ρ quantification using MR Multitasking: Method and initial clinical experience in tissue characterization of multiple sclerosis

PURPOSE

To develop a 3D whole-brain simultaneous T1/T2/T1ρ quantification method with MR Multitasking that provides high quality, co-registered multiparametric maps in 9 min.

METHODS

MR Multitasking conceptualizes T1/T2/T1ρ relaxations as different time dimensions, simultaneously resolving all three dimensions with a low-rank tensor image model. The proposed method was validated on a phantom and in healthy volunteers, comparing quantitative measurements against corresponding reference methods and evaluating the scan-rescan repeatability. Initial clinical validation was performed in age-matched relapsing-remitting multiple sclerosis (RRMS) patients to examine the feasibility of quantitative tissue characterization and to compare with the healthy control cohort. The feasibility of synthesizing six contrast-weighted images was also examined.

RESULTS

Our framework produced high quality, co-registered T1/T2/T1ρ maps that closely resemble the reference maps. Multitasking T1/T2/T1ρ measurements showed substantial agreement with reference measurements on the phantom and in healthy controls. Bland-Altman analysis indicated good in vivo repeatability of all three parameters. In RRMS patients, lesions were conspicuously delineated on all three maps and on four synthetic weighted images (T2-weighted, T2-FLAIR, double inversion recovery, and a novel "T1ρ-FLAIR" contrast). T1 and T2 showed significant differences for normal appearing white matter between patients and controls, while T1ρ showed significant differences for normal appearing white matter, cortical gray matter, and deep gray matter. The combination of three parameters significantly improved the differentiation between RRMS patients and healthy controls, compared to using any single parameter alone.

CONCLUSION

MR Multitasking simultaneously quantifies whole-brain T1/T2/T1ρ and is clinically promising for quantitative tissue characterization of neurological diseases, such as MS.

Magnetic resonance multitasking for multidimensional assessment of cardiovascular system: Development and feasibility study on the thoracic aorta

PURPOSE

To develop an MR multitasking-based multidimensional assessment of cardiovascular system (MT-MACS) with electrocardiography-free and navigator-free data acquisition for a comprehensive evaluation of thoracic aortic diseases.

METHODS

The MT-MACS technique adopts a low-rank tensor image model with a cardiac time dimension for phase-resolved cine imaging and a T₂ -prepared inversion-recovery dimension for multicontrast assessment. Twelve healthy subjects and 2 patients with thoracic aortic diseases were recruited for the study at 3 T, and both qualitative (image quality score) and quantitative (contrast-to-noise ratio between lumen and wall, lumen and wall area, and aortic strain index) analyses were performed in all healthy subjects. The overall image quality was scored based on a 4-point scale: 3, excellent; 2, good; 1, fair; and 0, poor. Statistical analysis was used to test the measurement agreement between MT-MACS and its corresponding 2D references.

RESULTS

The MT-MACS images reconstructed from acquisitions as short as 6 minutes demonstrated good or excellent image quality for bright-blood (2.58 ± 0.46), dark-blood (2.58 ± 0.50), and gray-blood (2.17 ± 0.53) contrast weightings, respectively. The contrast-to-noise ratios for the three weightings were 49.2 ± 12.8, 20.0 ± 5.8 and 2.8 ± 1.8, respectively. There were good agreements in the lumen and wall area (intraclass correlation coefficient = 0.993, P < .001 for lumen; intraclass correlation coefficient = 0.969, P < .001 for wall area) and strain (intraclass correlation coefficient = 0.947, P < .001) between MT-MACS and conventional 2D sequences.

CONCLUSION

The MT-MACS technique provides high-quality, multidimensional images for a comprehensive assessment of the thoracic aorta. Technical feasibility was demonstrated in healthy subjects and patients with thoracic aortic diseases. Further clinical validation is warranted.

Multiparametric Mapping Magnetic Resonance Imaging of Pancreatic Disease

Background

Current magnetic resonance imaging (MRI) of pancreatic disease is qualitative in nature. Quantitative imaging offers several advantages, including increased reproducibility and sensitivity to detect mild or diffuse disease. The role of multiparametric mapping MRI in characterizing various tissue types in pancreatic disease such as chronic pancreatitis (CP) and pancreatic ductal adenocarcinoma (PDAC) has rarely been evaluated.

Purpose

To evaluate the feasibility of multiparametric mapping [T1, T2, and apparent diffusion coefficient (ADC)] in defining tissue characteristics that occur in CP and PDAC to improve disease diagnosis.

Materials and Methods:

Pancreatic MRI was performed in 17 patients with PDAC undergoing therapy, 7 patients with CP, and 29 healthy volunteers with no pancreatic disease. T1 modified Look-Locker Inversion Recovery (T1 MOLLI), T2-prepared gradient-echo, and multi-slice single-shot echo-planar diffusion weighted imaging (SS-EPI DWI) sequences were used for data acquisition. Regions of interest (ROIs) of pancreas in PDAC, CP, and control subjects were outlined by an experienced radiologist. One-way analysis of variance (ANOVA) was used to compare the difference between groups and regions of the pancreas, and Tukey tests were used for multiple comparison testing within groups. Receiver operator characteristic (ROC) curves were analyzed, and the areas under the curves (AUCs) were calculated using single parameter and combined parameters, respectively.

Results

T1, T2, and ADC values of the entire pancreas among PDAC, CP, and control subjects; and between upstream and downstream portions of the pancreas in PDAC patients were all significantly different (p < 0.05). The AUC values were 0.90 for T1, 0.55 for T2, and 0.71 for ADC for independent prediction of PDAC. By combining T1, T2, and ADC, the AUC value was 0.94 (sensitivity 91.54%, specificity 85.81%, 95% CI: 0.92–0.96), which yielded higher accuracy than any one parameter only (p < 0.001).

Conclusion

Multiparametric mapping MRI is feasible for the evaluation of the differences between PDAC, CP, and normal pancreas tissues. The combination of multiple parameters of T1, T2, and ADC provides a higher accuracy than any single parameter alone in tissue characterization of the pancreas.

Six-dimensional quantitative DCE MR Multitasking of the entire abdomen: Method and application to pancreatic ductal adenocarcinoma

PURPOSE

To develop a quantitative DCE MRI technique enabling entire-abdomen coverage, free-breathing acquisition, 1-second temporal resolution, and T₁ -based quantification of contrast agent concentration and kinetic modeling for the characterization of pancreatic ductal adenocarcinoma (PDAC).

METHODS

Segmented FLASH readouts following saturation-recovery preparation with randomized 3D Cartesian undersampling was used for incoherent data acquisition. MR Multitasking was used to reconstruct 6-dimensional images with 3 spatial dimensions, 1 T₁ recovery dimension for dynamic T₁ quantification, 1 respiratory dimension to resolve respiratory motion, and 1 DCE time dimension to capture the contrast kinetics. Sixteen healthy subjects and 14 patients with pathologically confirmed PDAC were recruited for the in vivo studies, and kinetic parameters v_p , K^trans , v_e , and K_ep were evaluated for each subject. Intersession repeatability of Multitasking DCE was assessed in 8 repeat healthy subjects. One-way unbalanced analysis of variance was performed between control and patient groups.

RESULTS

In vivo studies demonstrated that v_p , K^trans , and K_ep of PDAC were significantly lower compared with nontumoral regions in the patient group (P = .002, .003, .004, respectively) and normal pancreas in the control group (P = .011, <.001, <.001, respectively), while v_e was significantly higher than nontumoral regions (P < .001) and healthy pancreas (P < .001). The kinetic parameters showed good in vivo repeatability (interclass correlation coefficient: v_p , 0.95; K^trans , 0.98; v_e , 0.96; K_ep , 0.99).

CONCLUSION

The proposed Multitasking DCE is promising for the quantification of vascular properties of PDAC. Quantitative DCE parameters were repeatable in vivo and showed significant differences between normal pancreas and both tumor and nontumoral regions in patients with PDAC.

Three-dimensional simultaneous brain T1 , T2 , and ADC mapping with MR Multitasking

PURPOSE

To develop a simultaneous T₁ , T₂ , and ADC mapping method that provides co-registered, distortion-free images and enables multiparametric quantification of 3D brain coverage in a clinically feasible scan time with the MR Multitasking framework.

METHODS

The T₁ /T₂ /diffusion weighting was generated by a series of T₂ preparations and diffusion preparations. The underlying multidimensional image containing 3 spatial dimensions, 1 T₁ weighting dimension, 1 T₂ -preparation duration dimension, 1 b-value dimension, and 1 diffusion direction dimension was modeled as a 5-way low-rank tensor. A separate real-time low-rank model incorporating time-resolved phase correction was also used to compensate for both inter-shot and intra-shot phase inconsistency induced by physiological motion. The proposed method was validated on both phantom and 16 healthy subjects. The quantification of T₁ /T₂ /ADC was evaluated for each case. Three post-surgery brain tumor patients were scanned for demonstration of clinical feasibility.

RESULTS

Multitasking T₁ /T₂ /ADC maps were perfectly co-registered and free from image distortion. Phantom studies showed substantial quantitative agreement ( R 2 = 0.999 ) with reference protocols for T₁ /T₂ /ADC. In vivo studies showed nonsignificant T₁ (P = .248), T₂ (P = .97), ADC (P = .328) differences among the frontal, parietal, and occipital regions. Although Multitasking showed significant differences of T₁ (P = .03), T₂ (P < .001), and ADC (P = .001) biases against the references, the mean bias estimates were small (ΔT₁ % < 5%, ΔT₂ % < 7%, ΔADC% < 5%), with all intraclass correlation coefficients greater than 0.82 indicating "excellent" agreement. Patient studies showed that Multitasking T₁ /T₂ /ADC maps were consistent with the clinical qualitative images.

CONCLUSION

The Multitasking approach simultaneously quantifies T₁ /T₂ /ADC with substantial agreement with the references and is promising for clinical applications.

Deep learning within a priori temporal feature spaces for large-scale dynamic MR image reconstruction: Application to 5-D cardiac MR Multitasking

High spatiotemporal resolution dynamic magnetic resonance imaging (MRI) is a powerful clinical tool for imaging moving structures as well as to reveal and quantify other physical and physiological dynamics. The low speed of MRI necessitates acceleration methods such as deep learning reconstruction from under-sampled data. However, the massive size of many dynamic MRI problems prevents deep learning networks from directly exploiting global temporal relationships. In this work, we show that by applying deep neural networks inside a priori calculated temporal feature spaces, we enable deep learning reconstruction with global temporal modeling even for image sequences with >40,000 frames. One proposed variation of our approach using dilated multi-level Densely Connected Network (mDCN) speeds up feature space coordinate calculation by 3000x compared to conventional iterative methods, from 20 minutes to 0.39 seconds. Thus, the combination of low-rank tensor and deep learning models not only makes large-scale dynamic MRI feasible but also practical for routine clinical application.

Free-breathing, non-ECG, continuous myocardial T1 mapping with cardiovascular magnetic resonance multitasking

PURPOSE

To evaluate the accuracy and repeatability of a free-breathing, non-electrocardiogram (ECG), continuous myocardial T₁ and extracellular volume (ECV) mapping technique adapted from the Multitasking framework.

METHODS

The Multitasking framework is adapted to quantify both myocardial native T₁ and ECV with a free-breathing, non-ECG, continuous acquisition T₁ mapping method. We acquire interleaved high-spatial resolution image data and high-temporal resolution auxiliary data following inversion-recovery pulses at set intervals and perform low-rank tensor imaging to reconstruct images at 344 inversion times, 20 cardiac phases, and 6 respiratory phases. The accuracy and repeatability of Multitasking T₁ mapping in generating native T₁ and ECV maps are compared with conventional techniques in a phantom, a simulation, 12 healthy subjects, and 10 acute myocardial infarction patients.

RESULTS

In phantoms, Multitasking T₁ mapping correlated strongly with the gold-standard spin-echo inversion recovery (R² = 0.99). A simulation study demonstrated that Multitasking T₁ mapping has similar myocardial sharpness to the fully sampled ground truth. In vivo native T₁ and ECV values from Multitasking T₁ mapping agree well with conventional MOLLI values and show good repeatability for native T₁ and ECV mapping for 60 seconds, 30 seconds, or 15 seconds of data. Multitasking native T₁ and ECV in myocardial infarction patients correlate positively with values from MOLLI.

CONCLUSION

Multitasking T₁ mapping can quantify native T₁ and ECV in the myocardium with free-breathing, non-ECG, continuous scans with good image quality and good repeatability in vivo in healthy subjects, and correlation with MOLLI T₁ and ECV in acute myocardial infarction patients.

Quantitative 3D dynamic contrast-enhanced (DCE) MR imaging of carotid vessel wall by fast T1 mapping using Multitasking

PURPOSE

To develop a dynamic contrast-enhanced (DCE) MRI method capable of high spatiotemporal resolution, 3D carotid coverage, and T1-based quantification of contrast agent concentration for the assessment of carotid atherosclerosis using a newly developed Multitasking technique.

METHODS

5D imaging with 3 spatial dimensions, 1 T1 recovery dimension, and 1 DCE time dimension was performed using MR Multitasking based on low-rank tensor modeling, which allows direct T1 quantification with high spatiotemporal resolution (0.7 mm isotropic and 595 ms, respectively). Saturation recovery preparations followed by 3D segmented fast low angle shot readouts were implemented with Gaussian-density random 3D Cartesian sampling. A bulk motion removal scheme was developed to improve image quality. The proposed protocol was tested in phantom and human studies. In vivo scans were performed on 14 healthy subjects and 7 patients with carotid atherosclerosis. Kinetic parameters including area under the concentration versus time curve (AUC), v_p , K^trans , and v_e were evaluated for each case.

RESULTS

Phantom experiments showed that T1 measurements using the proposed protocol were in good agreement with reference value ( R 2 = 0.96 ). In vivo studies demonstrated that AUC, v_p , and K^trans in the patient group were significantly higher than in the control group (0.63 ± 0.13 versus 0.42 ± 0.12, P < 0.001; 0.14 ± 0.05 versus 0.11 ± 0.03, P = 0.034; and 0.13 ± 0.04 versus 0.08 ± 0.02, P < 0.001, respectively). Results from repeated subjects showed good interscan reproducibility (intraclass correlation coefficient: v_p , 0.83; K^trans , 0.87; v_e , 0.92; AUC, 0.94).

CONCLUSION

Multitasking DCE is a promising approach for quantitatively assessing the vascularity properties of the carotid vessel wall.

Chemical exchange saturation transfer fingerprinting for exchange rate quantification

PURPOSE

There is an increased interest to determine the exchange rate using CEST to provide pH information. However, current CEST quantification methods require lengthy scan times and do not address magnetization transfer effects. The purpose of this work was to apply the magnetic resonance fingerprinting (MRF) concept to CEST to achieve more efficient and accurate exchange rate quantification.

METHODS

The proposed CEST fingerprinting method used varying saturation powers and saturation times to create unique signal evolutions for different exchange rates. The acquired signal was matched to a predefined dictionary to determine the exchange rate. The magnetization transfer effects were also addressed in the framework of CEST fingerprinting: The simulated dictionary could predict the signal curves without magnetization transfer effects, and comparing the dictionary to the acquired signals allowed the correction of the magnetization transfer effects. The CEST fingerprinting method was compared with the conventional pulsed quantitative CEST method using omega plots in the creatine phantom study.

RESULTS

The CEST fingerprinting method has a significantly reduced scan time (10 minutes versus 50 minutes) while providing more accurate exchange rate quantification using literature values as the reference.

CONCLUSION

In this study, we demonstrate that CEST fingerprinting is more efficient (5 times faster) compared with pulsed quantitative CEST. It is also shown that the results of the proposed CEST fingerprinting technique are much closer to the literature values than pulsed quantitative CEST at 3 T.

Accelerated Cardiac Diffusion Tensor Imaging Using Joint Low-Rank and Sparsity Constraints

OBJECTIVE

The purpose of this paper is to accelerate cardiac diffusion tensor imaging (CDTI) by integrating low-rankness and compressed sensing.

METHODS

Diffusion-weighted images exhibit both transform sparsity and low-rankness. These properties can jointly be exploited to accelerate CDTI, especially when a phase map is applied to correct for the phase inconsistency across diffusion directions, thereby enhancing low-rankness. The proposed method is evaluated both ex vivo and in vivo, and is compared to methods using either a low-rank or sparsity constraint alone.

RESULTS

Compared to using a low-rank or sparsity constraint alone, the proposed method preserves more accurate helix angle features, the transmural continuum across the myocardium wall, and mean diffusivity at higher acceleration, while yielding significantly lower bias and higher intraclass correlation coefficient.

CONCLUSION

Low-rankness and compressed sensing together facilitate acceleration for both ex vivo and in vivo CDTI, improving reconstruction accuracy compared to employing either constraint alone.

SIGNIFICANCE

Compared to previous methods for accelerating CDTI, the proposed method has the potential to reach higher acceleration while preserving myofiber architecture features, which may allow more spatial coverage, higher spatial resolution, and shorter temporal footprint in the future.

Accelerated High-Dimensional MR Imaging With Sparse Sampling Using Low-Rank Tensors

High-dimensional MR imaging often requires long data acquisition time, thereby limiting its practical applications. This paper presents a low-rank tensor based method for accelerated high-dimensional MR imaging using sparse sampling. This method represents high-dimensional images as low-rank tensors (or partially separable functions) and uses this mathematical structure for sparse sampling of the data space and for image reconstruction from highly undersampled data. More specifically, the proposed method acquires two datasets with complementary sampling patterns, one for subspace estimation and the other for image reconstruction; image reconstruction from highly undersampled data is accomplished by fitting the measured data with a sparsity constraint on the core tensor and a group sparsity constraint on the spatial coefficients jointly using the alternating direction method of multipliers. The usefulness of the proposed method is demonstrated in MRI applications; it may also have applications beyond MRI.

Fast dynamic electron paramagnetic resonance (EPR) oxygen imaging using low-rank tensors

Hypoxic tumors are resistant to radiotherapy, motivating the development of tools to image local oxygen concentrations. It is generally believed that stable or chronic hypoxia is the source of resistance, but more recent work suggests a role for transient hypoxia. Conventional EPR imaging (EPRI) is capable of imaging tissue pO2in vivo, with high pO2 resolution and 1mm spatial resolution but low imaging speed (10min temporal resolution for T1-based pO2 mapping), which makes it difficult to investigate the oxygen changes, e.g., transient hypoxia. Here we describe a new imaging method which accelerates dynamic EPR oxygen imaging, allowing 3D imaging at 2 frames per minute, fast enough to image transient hypoxia at the "speed limit" of observed pO2 change. The method centers on a low-rank tensor model that decouples the tradeoff between imaging speed, spatial coverage/resolution, and number of inversion times (pO2 accuracy). We present a specialized sparse sampling strategy and image reconstruction algorithm for use with this model. The quality and utility of the method is demonstrated in simulations and in vivo experiments in tumor bearing mice.

Improved subspace estimation for low-rank model-based accelerated cardiac imaging

Sparse sampling methods have emerged as effective tools to accelerate cardiac magnetic resonance imaging (MRI).Low-rank model-based cardiac imaging uses a predetermined temporal subspace for image reconstruction from highly undersampled (k, t)-space data and has been demonstrated effective for high-speed cardiac MRI. The accuracy of the temporal subspace isa key factor in these methods, yet little work has been published on data acquisition strategies to improve subspace estimation. This paper investigates the use of non-Cartesian k-space trajectories to replace the Cartesian trajectories that are omnipresent but are highly sensitive to readout direction. We also propose "self-navigated" pulse sequences that collect both navigator data (for determining the temporal subspace) and imaging data after every RF pulse, allowing for even greater acceleration. We investigate subspace estimation strategies through analysis of phantom images and demonstrate in vivo cardiac imaging in rats and mice without the use of ECG or respiratory gating. The proposed methods achieved 3-D imaging of wall motion, first-pass myocardial perfusion, and late gadolinium enhancement in rats at 74 frames/s,as well as 2-D imaging of wall motion in mice at 97 frames/s.

Self-navigated low-rank MRI for MPIO-labeled immune cell imaging of the heart

Super-paramagnetic iron oxide (SPIO) particles can magnetically label immune cells in circulation; the accumulation of labeled cells can then be detected by magnetic resonance imaging (MRI). This has enormous potential for imaging inflammatory responses in the heart, but it has been difficult to do in vivo using conventional free-breathing, ungated cardiac imaging. Subspace imaging with temporal navigation and sparse sampling of (k, t)-space has previously been used to accelerate several cardiac imaging applications, conventionally alternating between acquiring navigator data and sparse data every other TR. Here we describe a more efficient self-navigated pulse sequence to acquire both navigator and sparse (k, t)-space data in the space of a single TR, doubling imaging speed to approach 100 frames per second (fps). We show the feasibility of using the resulting method to assess myocardial inflammation in a pre-clinical rodent ischemic reperfusion injury (IRI) model using micron-sized paramagnetic iron oxide (MPIO) particles to label immune cells in situ.

High-resolution cardiovascular MRI by integrating parallel imaging with low-rank and sparse modeling

Magnetic resonance imaging (MRI) has long been recognized as a powerful tool for cardiovascular imaging because of its unique potential to measure blood flow, cardiac wall motion, and tissue properties jointly. However, many clinical applications of cardiac MRI have been limited by low imaging speed. In this paper, we present a novel method to accelerate cardiovascular MRI through the integration of parallel imaging, low-rank modeling, and sparse modeling. This method consists of a novel image model and specialized data acquisition. Of particular novelty is the proposed low-rank model component, which is specially adapted to the particular low-rank structure of cardiovascular signals. Simulations and in vivo experiments were performed to evaluate the method, as well as an analysis of the low-rank structure of a numerical cardiovascular phantom. Cardiac imaging experiments were carried out on both human and rat subjects without the use of ECG or respiratory gating and without breath holds. The proposed method reconstructed 2-D human cardiac images up to 22 fps and 1.0 mm × 1.0 mm spatial resolution and 3-D rat cardiac images at 67 fps and 0.65 mm × 0.65 mm × 0.31 mm spatial resolution. These capabilities will enhance the practical utility of cardiovascular MRI.