Journal of Cardiovascular Magnetic Resonance最新文献

Background: Automatic myocardial scar segmentation from late gadolinium enhancement (LGE) images using neural networks promises an alternative to time-consuming and observer-dependent semi-automatic approaches. However, alterations in data acquisition, reconstruction as well as post-processing may compromise network performance. The objective of the present work was to systematically assess network performance degradation due to a mismatch of point-spread function between training and testing data.

Methods: Thirty-six high-resolution (0.7×0.7×2.0 mm³) LGE k-space datasets were acquired post-mortem in porcine models of myocardial infarction. The in-plane point-spread function and hence in-plane resolution Δx was retrospectively degraded using k-space lowpass filtering, while field-of-view and matrix size were kept constant. Manual segmentation of the left ventricle (LV) and healthy remote myocardium was performed to quantify location and area (% of myocardium) of scar by thresholding (≥ SD5 above remote). Three standard U-Nets were trained on training resolutions Δx_train = 0.7, 1.2 and 1.7 mm to predict endo- and epicardial borders of LV myocardium and scar. The scar prediction of the three networks for varying test resolutions (Δx_test = 0.7 to 1.7 mm) was compared against the reference SD5 thresholding at 0.7 mm. Finally, a fourth network trained on a combination of resolutions (Δx_train = 0.7 to 1.7 mm) was tested.

Results: The prediction of relative scar areas showed the highest precision when the resolution of the test data was identical to or close to the resolution used during training. The median fractional scar errors and precisions (IQR) from networks trained and tested on the same resolution were 0.0 percentage points (p.p.) (1.24 - 1.45), and - 0.5 - 0.0 p.p. (2.00 - 3.25) for networks trained and tested on the most differing resolutions, respectively. Deploying the network trained on multiple resolutions resulted in reduced resolution dependency with median scar errors and IQRs of 0.0 p.p. (1.24 - 1.69) for all investigated test resolutions.

Conclusion: A mismatch of the imaging point-spread function between training and test data can lead to degradation of scar segmentation when using current U-Net architectures as demonstrated on LGE porcine myocardial infarction data. Training networks on multi-resolution data can alleviate the resolution dependency.

Background: Cardiovascular magnetic resonance (CMR) is an important imaging modality for the assessment and management of adult patients with congenital heart disease (CHD). However, conventional techniques for three-dimensional (3D) whole-heart acquisition involve long and unpredictable scan times and methods that accelerate scans via k-space undersampling often rely on long iterative reconstructions. Deep-learning-based reconstruction methods have recently attracted much interest due to their capacity to provide fast reconstructions while often outperforming existing state-of-the-art methods. In this study, we sought to adapt and validate a non-rigid motion-corrected model-based deep learning (MoCo-MoDL) reconstruction framework for 3D whole-heart MRI in a CHD patient cohort.

Methods: The previously proposed deep-learning reconstruction framework MoCo-MoDL, which incorporates a non-rigid motion-estimation network and a denoising regularization network within an unrolled iterative reconstruction, was trained in an end-to-end manner using 39 CHD patient datasets. Once trained, the framework was evaluated in eight CHD patient datasets acquired with seven-fold prospective undersampling. Reconstruction quality was compared with the state-of-the-art non-rigid motion-corrected patch-based low-rank reconstruction method (NR-PROST) and against reference images (acquired with three-or-four-fold undersampling and reconstructed with NR-PROST).

Results: Seven-fold undersampled scan times were 2.1 ± 0.3 minutes and reconstruction times were ∼30 seconds, approximately 240 times faster than an NR-PROST reconstruction. Image quality comparable to the reference images was achieved using the proposed MoCo-MoDL framework, with no statistically significant differences found in any of the assessed quantitative or qualitative image quality measures. Additionally, expert image quality scores indicated the MoCo-MoDL reconstructions were consistently of a higher quality than the NR-PROST reconstructions of the same data, with the differences in 12 of the 22 scores measured for individual vascular structures found to be statistically significant.

Conclusion: The MoCo-MoDL framework was applied to an adult CHD patient cohort, achieving good quality 3D whole-heart images from ∼2-minute scans with reconstruction times of ∼30 seconds.

Background: The prognostic value of follow-up cardiovascular magnetic resonance (CMR) in dilated cardiomyopathy (DCM) patients is unclear. We aimed to investigate the prognostic value of cardiac function, structure, and tissue characteristics at mid-term CMR follow-up.

Methods: The study population was a prospectively enrolled cohort of DCM patients who underwent guideline-directed medical therapy with baseline and follow-up CMR, which included measurement of biventricular volume and ejection fraction, late gadolinium enhancement, native T1, native T2, and extracellular volume. During follow-up, major adverse cardiac events (MACE) were defined as a composite endpoint of cardiovascular death, heart transplantation, and heart-failure readmission.

Results: Among 235 DCM patients (median CMR interval: 15.3 months; interquartile range: 12.5-19.2 months), 54 (23.0%) experienced MACE during follow-up (median: 31.2 months; interquartile range: 20.8-50.0 months). In multivariable Cox regression, follow-up CMR models showed significantly superior predictive value than baseline CMR models. Stepwise multivariate Cox regression showed that follow-up left ventricular ejection fraction (LVEF; hazard ratio [HR], 0.93; 95% confidence interval [CI], 0.91-0.96; p < 0.001) and native T1 (HR, 1.01; 95% CI, 1.00-1.01; p = 0.030) were independent predictors of MACE. Follow-up LVEF ≥ 40% or stable LVEF < 40% with T1 ≤ 1273 ms indicated low risk (annual event rate < 4%), while stable LVEF < 40% and T1 > 1273 ms or LVEF < 40% with deterioration indicated high risk (annual event rate > 15%).

Conclusions: Follow-up CMR provided better risk stratification than baseline CMR. Improvements in the LVEF and T1 mapping are associated with a lower risk of MACE.

Background: Ascending thoracic aortic aneurysm (ATAA) is a silent and threatening dilation of the ascending aorta (AscAo). Maximal aortic diameter which is currently used for ATAA patients management and surgery planning has been shown to inadequately characterize risk of dissection in a large proportion of patients. Our aim was to propose a comprehensive quantitative evaluation of aortic morphology and pressure-flow-wall associations from four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) data in healthy aging and in patients with ATAA.

Methods: We studied 17 ATAA patients (64.7 ± 14.3 years, 5 females) along with 17 age- and sex-matched healthy controls (59.7 ± 13.3 years, 5 females) and 13 younger healthy subjects (33.5 ± 11.1 years, 4 females). All subjects underwent a CMR exam, including 4D flow and three-dimensional anatomical images of the aorta. This latter dataset was used for aortic morphology measurements, including AscAo maximal diameter (iD_MAX) and volume, indexed to body surface area. 4D flow MRI data were used to estimate 1) cross-sectional local AscAo spatial (∆P_S) and temporal (∆P_T) pressure changes as well as the distance (∆D_PS) and time duration (∆T_PT) between local pressure peaks, 2) AscAo maximal wall shear stress (WSS_MAX) at peak systole, and 3) AscAo flow vorticity amplitude (V_MAX), duration (V_FWHM), and eccentricity (V_ECC).

Results: Consistency of flow and pressure indices was demonstrated through their significant associations with AscAo iD_MAX (WSS_MAX:r = -0.49, p < 0.001; V_ECC:r = -0.29, p = 0.045; V_FWHM:r = 0.48, p < 0.001; ∆D_PS:r = 0.37, p = 0.010; ∆T_PT:r = -0.52, p < 0.001) and indexed volume (WSS_MAX:r = -0.63, V_ECC:r = -0.51, V_FWHM:r = 0.53, ∆D_PS:r = 0.54, ∆T_PT:r = -0.63, p < 0.001 for all). Intra-AscAo cross-sectional pressure difference, ∆P_S, was significantly and positively associated with both V_MAX (r = 0.55, p = 0.002) and WSS_MAX (r = 0.59, p < 0.001) in the 30 healthy subjects (48.3 ± 18.0 years). Associations remained significant after adjustment for iD_MAX, age, and systolic blood pressure. Superimposition of ATAA patients to normal aging trends between ∆P_S and WSS_MAX as well as V_MAX allowed identifying patients with substantially high pressure differences concomitant with AscAo dilation.

Conclusion: Local variations in pressures within ascending aortic cross-sections derived from 4D flow MRI were associated with flow changes, as quantified by vorticity, and with stress exerted by blood on the aortic wall, as quantified by wall shear stress. Such flow-wall and pressure interactions might help for the identification

Background: Non-Cartesian magnetic resonance imaging trajectories at golden angle increments have the advantage of allowing motion correction and gating using intermediate real-time reconstructions. However, when the acquired data are cardiac binned for cine imaging, trajectories can cluster together at certain heart rates (HR) causing image artifacts. Here, we demonstrate an approach to reduce clustering by inserting additional angular increments within the trajectory, and optimizing them while still allowing for intermediate reconstructions.

Methods: Three acquisition models were simulated under constant and variable HR: golden angle (M_trd), random additional angles (M_rnd), and optimized additional angles (M_opt). The standard deviations of trajectory angular differences (STAD) were compared through their interquartile ranges (IQR) and the Kolmogorov-Smirnov test (significance level: p = 0.05). Agreement between an image reconstructed with uniform sampling and images from M_trd, M_rnd, and M_opt was analyzed using the structural similarity index measure (SSIM). M_trd and M_opt were compared in three adults at high, low, and no HR variability.

Results: STADs from M_trd were significantly different (p < 0.05) from M_opt and M_rnd. STAD (IQR × 10^-2 rad) showed that M_opt (0.5) and M_rnd (0.5) reduced clustering relative to M_trd (1.9) at constant HR. For variable HR, M_opt (0.5) and M_rnd (0.5) outperformed M_trd (0.9). The SSIM (IQR) showed that M_opt (0.011) produced the best image quality, followed by M_rnd (0.014), and M_trd (0.030). M_opt outperformed M_trd at reduced HR variability in in-vivo studies. At high HR variability, both models performed well.

Conclusion: This approach reduces clustering in k-space and improves image quality.

Background: 4D flow MRI enables assessment of cardiac function and intra-cardiac blood flow dynamics from a single acquisition. However, due to the poor contrast between the chambers and surrounding tissue, quantitative analysis relies on the segmentation derived from a registered cine MRI acquisition. This requires an additional acquisition and is prone to imperfect spatial and temporal inter-scan alignment. Therefore, in this work we developed and evaluated deep learning-based methods to segment the left ventricle (LV) from 4D flow MRI directly.

Methods: We compared five deep learning-based approaches with different network structures, data pre-processing and feature fusion methods. For the data pre-processing, the 4D flow MRI data was reformatted into a stack of short-axis view slices. Two feature fusion approaches were proposed to integrate the features from magnitude and velocity images. The networks were trained and evaluated on an in-house dataset of 101 subjects with 67,567 2D images and 3030 3D volumes. The performance was evaluated using various metrics including Dice, average surface distance (ASD), end-diastolic volume (EDV), end-systolic volume (ESV), LV ejection fraction (LVEF), LV blood flow kinetic energy (KE) and LV flow components. The Monte Carlo dropout method was used to assess the confidence and to describe the uncertainty area in the segmentation results.

Results: Among the five models, the model combining 2D U-Net with late fusion method operating on short-axis reformatted 4D flow volumes achieved the best results with Dice of 84.52% and ASD of 3.14 mm. The best averaged absolute and relative error between manual and automated segmentation for EDV, ESV, LVEF and KE was 19.93 ml (10.39%), 17.38 ml (22.22%), 7.37% (13.93%) and 0.07 mJ (5.61%), respectively. Flow component results derived from automated segmentation showed high correlation and small average error compared to results derived from manual segmentation.

Conclusions: Deep learning-based methods can achieve accurate automated LV segmentation and subsequent quantification of volumetric and hemodynamic LV parameters from 4D flow MRI without requiring an additional cine MRI acquisition.

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 T₁-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 (HIP_vol). Specifically, HIP_vol 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, HIP_vol 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 HIP_vol 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 HIP_vol 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.

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.

Background: Cardiovascular magnetic resonance (CMR) imaging is an important tool for evaluating the severity of aortic stenosis (AS), co-existing aortic disease, and concurrent myocardial abnormalities. Acquiring this additional information requires protocol adaptations and additional scanner time, but is not necessary for the majority of patients who do not have AS. We observed that the relative signal intensity of blood in the ascending aorta on a balanced steady state free precession (bSSFP) 3-chamber cine was often reduced in those with significant aortic stenosis. We investigated whether this effect could be quantified and used to predict AS severity in comparison to existing gold-standard measurements.

Methods: Multi-centre, multi-vendor retrospective analysis of patients with AS undergoing CMR and transthoracic echocardiography (TTE). Blood signal intensity was measured in a ∼1 cm² region of interest (ROI) in the aorta and left ventricle (LV) in the 3-chamber bSSFP cine. Because signal intensity varied across patients and scanner vendors, a ratio of the mean signal intensity in the aorta ROI to the LV ROI (Ao:LV) was used. This ratio was compared using Pearson correlations against TTE parameters of AS severity: aortic valve peak velocity, mean pressure gradient and the dimensionless index. The study also assessed whether field strength (1.5 T vs. 3 T) and patient characteristics (presence of bicuspid aortic valves (BAV), dilated aortic root and low flow states) altered this signal relationship.

Results: 314 patients (median age 69 [IQR 57-77], 64% male) who had undergone both CMR and TTE were studied; 84 had severe AS, 78 had moderate AS, 66 had mild AS and 86 without AS were studied as a comparator group. The median time between CMR and TTE was 12 weeks (IQR 4-26). The Ao:LV ratio at 1.5 T strongly correlated with peak velocity (r = -0.796, p = 0.001), peak gradient (r = -0.772, p = 0.001) and dimensionless index (r = 0.743, p = 0.001). An Ao:LV ratio of < 0.86 was 84% sensitive and 82% specific for detecting AS of any severity and a ratio of 0.58 was 83% sensitive and 92% specific for severe AS. The ability of Ao:LV ratio to predict AS severity remained for patients with bicuspid aortic valves, dilated aortic root or low indexed stroke volume. The relationship between Ao:LV ratio and AS severity was weaker at 3 T.

Conclusions: The Ao:LV ratio, derived from bSSFP 3-chamber cine images, shows a good correlation with existing measures of AS severity. It demonstrates utility at 1.5 T and offers an easily calculable metric that can be used at the time of scanning or automated to identify on an adaptive basis which patients benefit from dedicated imaging to assess which patients should have additional sequences to assess AS.

Background: The prognostic value of left ventricular (LV) myocardial trabecular complexity on cardiovascular magnetic resonance (CMR) in dilated cardiomyopathy (DCM) remains unknown. This study aimed to evaluate the prognostic value of LV myocardial trabecular complexity using fractal analysis in patients with DCM.

Methods: Consecutive patients with DCM who underwent CMR between March 2017 and November 2021 at two hospitals were prospectively enrolled. The primary endpoints were defined as the combination of all-cause death and heart failure hospitalization. The events of cardiac death alone were defined as the secondary endpoints.LV trabeculae complexity was quantified by measuring the fractal dimension (FD) of the endocardial border based on fractal geometry on CMR. Cox proportional hazards regression and Kaplan-Meier survival analysis were used to examine the association between variables and outcomes. The incremental prognostic value of FD was assessed in nested models.

Results: A total of 403 patients with DCM (49.31 ± 14.68 years, 69% male) were recruited. After a median follow-up of 43 months (interquartile range, 28-55 months), 87 and 24 patients reached the primary and secondary endpoints, respectively. Age, heart rate, New York Heart Association functional class >II, N-terminal pro-B-type natriuretic peptide, LV ejection fraction, LV end-diastolic volume index, LV end-systolic volume index, LV mass index, presence of late gadolinium enhancement, global FD, LV mean apical FD, and LV maximal apical FD were univariably associated with the outcomes (all P < 0.05). After multivariate adjustment, LV maximal apical FD remained a significant independent predictor of outcome [hazard ratio = 1.179 (1.116, 1.246), P < 0.001]. The addition of LV maximal apical FD in the nested models added incremental prognostic value to other common clinical and imaging risk factors (all <0.001; C-statistic: 0.84-0.88, P < 0.001).

Conclusion: LV maximal apical FD was an independent predictor of the adverse clinical outcomes in patients with DCM and provided incremental prognostic value over conventional clinical and imaging risk factors.