K-t GRAPPA accelerated phase contrast MRI: Improved assessment of blood flow and 3-directional myocardial motion during breath-hold

The role of 4-dimensional flow in the assessment of bicuspid aortic valve and its valvulo-aortopathies

Bicuspid aortic valve is the most common congenital cardiac malformation and the leading cause of aortopathy and aortic stenosis in younger patients. Aortic wall remodelling secondary to altered haemodynamic flow patterns, changes in peak velocity, and wall shear stress may be implicated in the development of aortopathy in the presence of bicuspid aortic valve and dysfunction. Assessment of these parameters as potential predictors of disease severity and progression is thus desirable. The anatomic and functional information acquired from 4D flow MRI can allow simultaneous visualisation and quantification of the pathological geometric and haemodynamic changes of the aorta. We review the current clinical utility of haemodynamic quantities including velocity, wall sheer stress and energy losses, as well as visual descriptors such as vorticity and helicity, and flow direction in assessing the aortic valve and associated aortopathies.

Myocardial Strain Evaluation with Cardiovascular MRI: Physics, Principles, and Clinical Applications

Myocardial strain is a measure of myocardial deformation, which is a more sensitive imaging biomarker of myocardial disease than the commonly used ventricular ejection fraction. Although myocardial strain is commonly evaluated by using speckle-tracking echocardiography, cardiovascular MRI (CMR) is increasingly performed for this purpose. The most common CMR technique is feature tracking (FT), which involves postprocessing of routinely acquired cine MR images. Other CMR strain techniques require dedicated sequences, including myocardial tagging, strain-encoded imaging, displacement encoding with stimulated echoes, and tissue phase mapping. The complex systolic motion of the heart can be resolved into longitudinal strain, circumferential strain, radial strain, and torsion. Myocardial strain metrics include strain, strain rate, displacement, velocity, torsion, and torsion rate. Wide variability exists in the reference ranges for strain dependent on the imaging technique, analysis software, operator, patient demographics, and hemodynamic factors. In anticancer therapy cardiotoxicity, CMR myocardial strain can help identify left ventricular dysfunction before the decline of ejection fraction. CMR myocardial strain is also valuable for identifying patients with left ventricle dyssynchrony who will benefit from cardiac resynchronization therapy. CMR myocardial strain is also useful in ischemic heart disease, cardiomyopathies, pulmonary hypertension, and congenital heart disease. The authors review the physics, principles, and clinical applications of CMR strain techniques. Online supplemental material is available for this article. ^©RSNA, 2022.

Automated segmentation of biventricular contours in tissue phase mapping using deep learning

Tissue phase mapping (TPM) is an MRI technique for quantification of regional biventricular myocardial velocities. Despite its potential, clinical use is limited due to the requisite labor-intensive manual segmentation of cardiac contours for all time frames. The purpose of this study was to develop a deep learning (DL) network for automated segmentation of TPM images, without significant loss in segmentation and myocardial velocity quantification accuracy compared with manual segmentation. We implemented a multi-channel 3D (three dimensional; 2D + time) dense U-Net that trained on magnitude and phase images and combined cross-entropy, Dice, and Hausdorff distance loss terms to improve the segmentation accuracy and suppress unnatural boundaries. The dense U-Net was trained and tested with 150 multi-slice, multi-phase TPM scans (114 scans for training, 36 for testing) from 99 heart transplant patients (44 females, 1-4 scans/patient), where the magnitude and velocity-encoded (V_x , V_y , V_z ) images were used as input and the corresponding manual segmentation masks were used as reference. The accuracy of DL segmentation was evaluated using quantitative metrics (Dice scores, Hausdorff distance) and linear regression and Bland-Altman analyses on the resulting peak radial and longitudinal velocities (V_r and V_z ). The mean segmentation time was about 2 h per patient for manual and 1.9 ± 0.3 s for DL. Our network produced good accuracy (median Dice = 0.85 for left ventricle (LV), 0.64 for right ventricle (RV), Hausdorff distance = 3.17 pixels) compared with manual segmentation. Peak V_r and V_z measured from manual and DL segmentations were strongly correlated (R ≥ 0.88) and in good agreement with manual analysis (mean difference and limits of agreement for V_z and V_r were -0.05 ± 0.98 cm/s and -0.06 ± 1.18 cm/s for LV, and -0.21 ± 2.33 cm/s and 0.46 ± 4.00 cm/s for RV, respectively). The proposed multi-channel 3D dense U-Net was capable of reducing the segmentation time by 3,600-fold, without significant loss in accuracy in tissue velocity measurements.

Carotid Geometry and Wall Shear Stress Independently Predict Increased Wall Thickness-A Longitudinal 3D MRI Study in High-Risk Patients

Introduction: Carotid geometry and wall shear stress (WSS) have been proposed as independent risk factors for the progression of carotid atherosclerosis, but this has not yet been demonstrated in larger longitudinal studies. Therefore, we investigated the impact of these biomarkers on carotid wall thickness in patients with high cardiovascular risk.

Methods: Ninety-seven consecutive patients with hypertension, at least one additional cardiovascular risk factor and internal carotid artery (ICA) plaques (wall thickness ≥ 1.5 mm and degree of stenosis ≤ 50%) were prospectively included. They underwent high-resolution 3D multi-contrast and 4D flow MRI at 3 Tesla both at baseline and follow-up. Geometry (ICA/common carotid artery (CCA)-diameter ratio, bifurcation angle, tortuosity and wall thickness) and hemodynamics [WSS, oscillatory shear index (OSI)] of both carotid bifurcations were measured at baseline. Their predictive value for changes of wall thickness 12 months later was calculated using linear regression analysis for the entire study cohort (group 1, 97 patients) and after excluding patients with ICA stenosis ≥10% to rule out relevant inward remodeling (group 2, 61 patients).

Results: In group 1, only tortuosity at baseline was independently associated with carotid wall thickness at follow-up (regression coefficient = −0.52, p < 0.001). However, after excluding patients with ICA stenosis ≥10% in group 2, both ICA/CCA-ratio (0.49, p < 0.001), bifurcation angle (0.04, p = 0.001), tortuosity (−0.30, p = 0.040), and WSS (−0.03, p = 0.010) at baseline were independently associated with changes of carotid wall thickness at follow-up.

Conclusions: A large ICA bulb and bifurcation angle and low WSS seem to be independent risk factors for the progression of carotid atherosclerosis in the absence of ICA stenosis. By contrast, a high carotid tortuosity seems to be protective both in patients without and with ICA stenosis. These biomarkers may be helpful for the identification of patients who are at particular risk of wall thickness progression and who may benefit from intensified monitoring and treatment.

Quantitative assessment of velocity and flow using compressed SENSE in children and young adults with adequate acquired temporal resolution

BACKGROUND

Phase contrast (PC) cardiovascular magnetic resonance (CMR) imaging with parallel imaging acceleration is established and validated for measuring velocity and flow. However, additional acceleration to further shorten acquisition times would be beneficial in patients with complex vasculature who need multiple PC-CMR measurements, especially pediatric patients with higher heart rates.

METHODS

PC-CMR images acquired with compressed sensitivity encoding (C-SENSE) factors of 3 to 6 and standard of care PC-CMR with sensitivity encoding (SENSE) factor of 2 (S2) acquired as part of clinical CMR examinations performed between November 2020 and January 2021 were analyzed retrospectively. The velocity and flow through the ascending aorta (AAo), descending aorta (DAo), and superior vena cava (SVC) in a transverse plane at the level of pulmonary artery bifurcation were compared. Additionally, frequency power distribution and dynamic time warp distance were calculated for these acquisitions. To further validate the adequate temporal resolution requirement, patients with S2 PC-CMR in the same acquisition plane were added in frequency power distribution analysis.

RESULTS

Twenty-eight patients (25 males; 15.9 ± 1.9 years; body surface area (BSA) 1.7 ± 0.2 m2; heart rate 81 ± 16 bpm) underwent all five PC-CMR acquisitions during the study period. An additional 22 patients (16 males; 17.5 ± 7.7 years; BSA 1.6 ± 0.5 m2; heart rate 91 ± 16 bpm) were included for frequency power spectrum analysis. As expected, scan time decreased with increasing C-SENSE acceleration factor = 3 (37.5 ± 6.5 s, 26.4 ± 7.6%), 4 (28.1 ± 4.9 s, 44.7 ± 5.6%), 5 (21.6 ± 3.6 s, 57.6 ± 4.4%), and 6 (19.1 ± 3.2 s, 62.3 ± 4.2%) relative to SENSE = 2 (51.3 ± 10.1 s) PC-CMR acquisition. Mean peak velocity, net flow, and cardiac output were comparable (p > 0.87) between the five PC-CMR acquisitions with mean differences less than < 4%, < 2%, and < 3% respectively. All individual blood vessels showed a non-significant dependence of difference in fmax99 (< 4 Hz, p > 0.2), and dynamic time warp distance (p > 0.3) on the C-SENSE acceleration factor used. There was a strongly correlated (r = 0.74) increase in fmax99 (10.5 ± 2.2, range: 7.1-16.4 Hz) with increasing heart rate. The computed minimum required cardiac phase number was 15 ± 2.0 (range: 11-20) over the heart rate of 86 ± 15 bpm (range: 58-113 bpm).

CONCLUSIONS

Stroke volume, cardiac output, and mean peak velocity measurements using PC-CMR with C-SENSE of up to 6 agree with measurements by standard of care PC-CMR with SENSE = 2 and resulted in up to a 65% reduction in acquisition time. Adequate temporal sampling can be ensured by acquiring 20 cardiac phases throughout the entire cardiac cycle over a wide range of pediatric and young adult heart rates.

Carotid geometry is an independent predictor of wall thickness - a 3D cardiovascular magnetic resonance study in patients with high cardiovascular risk

BACKGROUND

The posterior wall of the proximal internal carotid artery (ICA) is the predilection site for the development of stenosis. To optimally prevent stroke, identification of new risk factors for plaque progression is of high interest. Therefore, we studied the impact of carotid geometry and wall shear stress on cardiovascular magnetic resonance (CMR)-depicted wall thickness in the ICA of patients with high cardiovascular disease risk.

METHODS

One hundred twenty-one consecutive patients ≥50 years with hypertension, ≥1 additional cardiovascular risk factor and ICA plaque ≥1.5 mm thickness and < 50% stenosis were prospectively included. High-resolution 3D-multi-contrast (time of flight, T1, T2, proton density) and 4D flow CMR were performed for the assessment of morphological (bifurcation angle, ICA/common carotid artery (CCA) diameter ratio, tortuosity, and wall thickness) and hemodynamic parameters (absolute/systolic wall shear stress (WSS), oscillatory shear index (OSI)) in 242 carotid bifurcations.

RESULTS

We found lower absolute/systolic WSS, higher OSI and increased wall thickness in the posterior compared to the anterior wall of the ICA bulb (p < 0.001), whereas this correlation disappeared in ≥10% stenosis. Higher carotid tortuosity (regression coefficient = 0.764; p < 0.001) and lower ICA/CCA diameter ratio (regression coefficient = - 0.302; p < 0.001) were independent predictors of increased wall thickness even after adjustment for cardiovascular risk factors. This association was not found for bifurcation angle, WSS or OSI in multivariate regression analysis.

CONCLUSIONS

High carotid tortuosity and low ICA diameter were independent predictors for wall thickness of the ICA bulb in this cross-sectional study, whereas this association was not present for WSS or OSI. Thus, consideration of geometric parameters of the carotid bifurcation could be helpful to identify patients at increased risk of carotid plaque generation. However, this association and the potential benefit of WSS measurement need to be further explored in a longitudinal study.

4D Flow with MRI

Magnetic resonance imaging (MRI) has become an important tool for the clinical evaluation of patients with cardiac and vascular diseases. Since its introduction in the late 1980s, quantitative flow imaging with MRI has become a routine part of standard-of-care cardiothoracic and vascular MRI for the assessment of pathological changes in blood flow in patients with cardiovascular disease. More recently, time-resolved flow imaging with velocity encoding along all three flow directions and three-dimensional (3D) anatomic coverage (4D flow MRI) has been developed and applied to enable comprehensive 3D visualization and quantification of hemodynamics throughout the human circulatory system. This article provides an overview of the use of 4D flow applications in different cardiac and vascular regions in the human circulatory system, with a focus on using 4D flow MRI in cardiothoracic and cerebrovascular diseases.

Autocalibrated cardiac tissue phase mapping with multiband imaging and k-t acceleration

PURPOSE

To develop an autocalibrated multiband (MB) CAIPIRINHA acquisition scheme with in-plane k-t acceleration enabling multislice three-directional tissue phase mapping in one breath-hold.

METHODS

A k-t undersampling scheme was integrated into a time-resolved electrocardiographic-triggered autocalibrated MB gradient-echo sequence. The sequence was used to acquire data on 4 healthy volunteers with MB factors of two (MB2) and three (MB3), which were reconstructed using a joint reconstruction algorithm that tackles both k-t and MB acceleration. Forward simulations of the imaging process were used to tune the reconstruction model hyperparameters. Direct comparisons between MB and single-band tissue phase-mapping measurements were performed.

RESULTS

Simulations showed that the velocities could be accurately reproduced with MB2 k-t (average ± twice the SD of the RMS error of 0.08 ± 0.22 cm/s and velocity peak reduction of 1.03% ± 6.47% compared with fully sampled velocities), whereas acceptable results were obtained with MB3 k-t (RMS error of 0.13 ± 0.58 cm/s and peak reduction of 2.21% ± 13.45%). When applied to tissue phase-mapping data, the proposed technique allowed three-directional velocity encoding to be simultaneously acquired at two/three slices in a single breath-hold of 18 heartbeats. No statistically significant differences were detected between MB2/MB3 k-t and single-band k-t motion traces averaged over the myocardium. Regional differences were found, however, when using the American Heart Association model for segmentation.

CONCLUSION

An autocalibrated MB k-t acquisition/reconstruction framework is presented that allows three-directional velocity encoding of the myocardial velocities at multiple slices in one breath-hold.

Evaluating Biventricular Myocardial Velocity and Interventricular Dyssynchrony in Adult Patients During the First Year After Heart Transplantation

BACKGROUND

Magnetic resonance tissue phase mapping (TPM) measures three-directional myocardial velocities of the left and right ventricle (LV, RV). This noninvasive technique may supplement endomyocardial biopsy (EMB) in monitoring grafts post-heart transplantation (HTx).

PURPOSE

To assess biventricular myocardial velocity alterations in grafts and investigate the relationship between velocities and acute cellular rejection (ACR) episodes.

STUDY TYPE

Prospective.

SUBJECTS

Twenty-seven patients within 1 year post-HTx (49 ± 13 years, 19 M) and 18 age-matched controls (49 ± 15 years, 12 M).

FIELD STRENGTH/SEQUENCE

1.5T, 2D balanced steady-state free precession, and TPM.

ASSESSMENT

Ventricular function: end-diastolic and end-systolic volumes, stroke volumes, ejection fraction (EF), and myocardial mass. TPM velocities: peak-systolic and peak-diastolic velocities, cardiac twist, and interventricular dyssynchrony. ACR rejection episodes: International Society for Heart and Lung Transplantation grading of EMB specimens.

STATISTICAL TESTS

The Lilliefors test for normality, unpaired t-tests, and Wilcoxon rank-sum tests for normally and nonnormally distributed data, respectively, were used, as well as multivariate regression for confounding variables and Pearson's correlation for associations between TPM velocities and global function.

RESULTS

Compared to controls, HTx patients demonstrated reduced biventricular systolic longitudinal velocities (LV: 5.2 ± 2.1 vs. 4.0 ± 1.5 cm/s, P < 0.05; RV: 4.2 ± 1.3 vs. 3.1 ± 1.2 cm/s, P < 0.01). Correlation analysis revealed significant positive relationships for biventricular EF with radial peak velocities of the same ventricle in both systole and diastole (LV systole: r = 0.48, P < 0.01; LV diastole: r = 0.28, P < 0.05; RV systole: r = 0.35, P < 0.01; RV diastole: r = 0.36, P < 0.01). Segmentally, longitudinal velocities were impaired in 7/16 LV segments and 5/10 RV segments in systole and 7/10 RV segments in diastole. TPM analysis in studies with >4 preceding ACR episodes showed globally reduced RV and LV systolic radial velocity, and segmentally reduced radial and longitudinal systolic velocities.

DATA CONCLUSION

Biventricular global and segmental velocities were reduced in HTx patients. Patients with >4 rejection episodes showed reduced myocardial velocities. The TPM sequence may add functional information for monitoring graft dysfunction.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY STAGE: 2 J. Magn. Reson. Imaging 2020;52:920-929.

Impact of age and cardiac disease on regional left and right ventricular myocardial motion in healthy controls and patients with repaired tetralogy of fallot

The assessment of both left (LV) and right ventricular (RV) motion is important to understand the impact of heart disease on cardiac function. The MRI technique of tissue phase mapping (TPM) allows for the quantification of regional biventricular three-directional myocardial velocities. The goal of this study was to establish normal LV and RV velocity parameters across a wide range of pediatric to adult ages and to investigate the feasibility of TPM for detecting impaired regional biventricular function in patients with repaired tetralogy of Fallot (TOF). Thirty-six healthy controls (age = 1-75 years) and 12 TOF patients (age = 5-23 years) underwent cardiac MRI including TPM in short-axis locations (base, mid, apex). For ten adults, a second TPM scan was used to assess test-retest reproducibility. Data analysis included the calculation of biventricular radial, circumferential, and long-axis velocity components, quantification of systolic and diastolic peak velocities in an extended 16 + 10 LV + RV segment model, and assessment of inter-ventricular dyssynchrony. Biventricular velocities showed good test-retest reproducibility (mean bias ≤ 0.23 cm/s). Diastolic radial and long-axis peak velocities for LV and RV were significantly reduced in adults compared to children (19-61%, p < 0.001-0.02). In TOF patients, TPM identified significantly reduced systolic and diastolic LV and RV long-axis peak velocities (20-50%, p < 0.001-0.05) compared to age-matched controls. In conclusion, tissue phase mapping enables comprehensive analysis of global and regional biventricular myocardial motion. Changes in myocardial velocities associated with age underline the importance of age-matched controls. This pilot study in TOF patients shows the feasibility to detect regionally abnormal LV and RV motion.

Segmental biventricular analysis of myocardial function using high temporal and spatial resolution tissue phase mapping

OBJECTIVE

Myocardial dysfunction of the right ventricle (RV) is an important indicator of RV diseases, e.g. RV infarction or pulmonary hypertension. Tissue phase mapping (TPM) has been widely used to determine function of the left ventricle (LV) by analyzing myocardial velocities. The analysis of RV motion is more complicated due to the different geometry and smaller wall thickness. The aim of this work was to adapt and optimize TPM to the demands of the RV.

MATERIALS AND METHODS

TPM measurements were acquired in 25 healthy volunteers using a velocity-encoded phase-contrast sequence and kt-accelerated parallel imaging in combination with optimized navigator strategy and blood saturation. Post processing was extended by a 10-segment RV model and a detailed biventricular analysis of myocardial velocities was performed.

RESULTS

High spatio-temporal resolution (1.0 × 1.0 × 6 mm³, 21.3 ms) and the optimized blood saturation enabled good delineation of the RV and its velocities. Global and segmental velocities, as well as time to peak velocities showed significant differences between the LV and RV. Furthermore, complex timing of the RV could be demonstrated by segmental time to peak analysis.

CONCLUSION

High spatio-temporal resolution TPM enables a detailed biventricular analysis of myocardial motion and might provide a reliable tool for description and detection of diseases affecting left and right ventricular function.

Improved k-t PCA Algorithm Using Artificial Sparsity in Dynamic MRI

The k-t principal component analysis (k-t PCA) is an effective approach for high spatiotemporal resolution dynamic magnetic resonance (MR) imaging. However, it suffers from larger residual aliasing artifacts and noise amplification when the reduction factor goes higher. To further enhance the performance of this technique, we propose a new method called sparse k-t PCA that combines the k-t PCA algorithm with an artificial sparsity constraint. It is a self-calibrated procedure that is based on the traditional k-t PCA method by further eliminating the reconstruction error derived from complex subtraction of the sampled k-t space from the original reconstructed k-t space. The proposed method is tested through both simulations and in vivo datasets with different reduction factors. Compared to the standard k-t PCA algorithm, the sparse k-t PCA can improve the normalized root-mean-square error performance and the accuracy of temporal resolution. It is thus useful for rapid dynamic MR imaging.

Magnetic resonance tissue phase mapping demonstrates altered left ventricular diastolic function in children with chronic kidney disease

BACKGROUND

Echocardiographic examinations have revealed functional cardiac abnormalities in children with chronic kidney disease.

OBJECTIVE

To assess the feasibility of MRI tissue phase mapping in children and to assess regional left ventricular wall movements in children with chronic kidney disease.

MATERIALS AND METHODS

Twenty pediatric patients with chronic kidney disease (before or after renal transplantation) and 12 healthy controls underwent tissue phase mapping (TPM) to quantify regional left ventricular function through myocardial long (Vz) and short-axis (Vr) velocities at all 3 levels of the left ventricle.

RESULTS

Patients and controls (age: 8 years-20 years) were matched for age, height, weight, gender and heart rate. Patients had higher systolic blood pressure. No patient had left ventricular hypertrophy on MRI or diastolic dysfunction on echocardiography. Fifteen patients underwent tissue Doppler echocardiography, with normal z-scores for mitral early diastolic (V_E), late diastolic (V_A) and peak systolic (V_S) velocities. Throughout all left ventricular levels, peak diastolic Vz and Vr (cm/s) were reduced in patients: Vz_base -10.6 ± 1.9 vs. -13.4 ± 2.0 (P < 0.0003), Vz_mid -7.8 ± 1.6 vs. -11 ± 1.5 (P < 0.0001), Vz_apex -3.8 ± 1.6 vs. -5.3 ± 1.6 (P = 0.01), Vr_base -4.2 ± 0.8 vs. -4.9 ± 0.7 (P = 0.01), Vr_mid -4.7 ± 0.7 vs. -5.4 ± 0.7 (P = 0.01), Vr_apex -4.7 ± 1.4 vs. -5.6 ± 1.1 (P = 0.05).

CONCLUSION

Tissue phase mapping is feasible in children and adolescents. Children with chronic kidney disease show significantly reduced peak diastolic long- and short-axis left ventricular wall velocities, reflecting impaired early diastolic filling. Thus, tissue phase mapping detects chronic kidney disease-related functional myocardial changes before overt left ventricular hypertrophy or echocardiographic diastolic dysfunction occurs.

Evaluation of Aortic Blood Flow and Wall Shear Stress in Aortic Stenosis and Its Association With Left Ventricular Remodeling

BACKGROUND

Aortic stenosis (AS) leads to variable stress for the left ventricle (LV) and consequently a broad range of LV remodeling. The aim of this study was to describe blood flow patterns in the ascending aorta of patients with AS and determine their association with remodeling.

METHODS AND RESULTS

Thirty-seven patients with AS (14 mild, 8 moderate, 15 severe; age, 63±13 years) and 37 healthy controls (age, 60±10 years) underwent 4-dimensional-flow magnetic resonance imaging. Helical and vortical flow formations and flow eccentricity were assessed in the ascending aorta. Normalized flow displacement from the vessel center and peak systolic wall shear stress in the ascending aorta were quantified. LV remodeling was assessed based on LV mass index and the ratio of LV mass:end-diastolic volume (relative wall mass). Marked helical and vortical flow formation and eccentricity were more prevalent in patients with AS than in healthy subjects, and patients with AS exhibited an asymmetrical and elevated distribution of peak systolic wall shear stress. In AS, aortic orifice area was strongly negatively associated with vortical flow formation (P=0.0274), eccentricity (P=0.0070), and flow displacement (P=0.0021). Bicuspid aortic valve was associated with more intense helical (P=0.0098) and vortical flow formation (P=0.0536), higher flow displacement (P=0.11), and higher peak systolic wall shear stress (P=0.0926). LV mass index and relative wall mass were significantly associated with aortic orifice area (P=0.0611, P=0.0058) and flow displacement (P=0.0058, P=0.0283).

CONCLUSIONS

In this pilot study, AS leads to abnormal blood flow pattern and peak systolic wall shear stress in the ascending aorta. In addition to aortic orifice area, normalized flow displacement was significantly associated with LV remodeling.

Efficient method for volumetric assessment of peak blood flow velocity using 4D flow MRI

PURPOSE

To test the feasibility and effectiveness of using maximum intensity plots (MIPs) based on 4D flow magnetic resonance imaging (MRI) velocity data to assess systolic peak velocities in a cohort of bicuspid aortic valve (BAV) patients.

MATERIALS AND METHODS

4D flow MRI at 1.5T was performed on 51 BAV patients. MIPs were generated from the 4D flow MRI velocity data and used by two users to determine peak velocities in three regions of interest (ROIs): ascending aorta (AAo), aortic arch, and descending aorta. 4D flow MRI peak velocities in the AAo were compared to peak velocities recorded by 2D phase contrast MRI (2D PCMRI) in a subcohort of 36 patients and by Doppler echocardiography in a subcohort of 34 patients. 4D flow MRI peak velocities recorded by each observer were compared for all ROIs to test for interobserver variability.

RESULTS

4D flow MRI recorded significantly higher velocities compared to 2D PCMRI (2.04 ± 0.71 m/s vs. 1.69 ± 0.79 m/s, 17.2% difference, P < 0.001) and similar velocities compared to Doppler echocardiography. There was excellent agreement between the observers, with a mean difference of 0.005 m/s and an intraclass correlation coefficient of 0.98.

CONCLUSION

4D flow MRI velocity MIPs allow for efficient measurement of peak velocities in BAV patients with higher accuracy than 2D PCMRI and similar accuracy to Doppler echocardiography. J. Magn. Reson. Imaging 2016;44:1673-1682.

Analyzing myocardial torsion based on tissue phase mapping cardiovascular magnetic resonance

BACKGROUND

The purpose of this work is to analyze differences in left ventricular torsion between volunteers and patients with non-ischemic cardiomyopathy based on tissue phase mapping (TPM) cardiovascular magnetic resonance (CMR).

METHODS

TPM was performed on 27 patients with non-ischemic cardiomyopathy and 14 normal volunteers. Patients underwent a standard CMR including late gadolinium enhancement (LGE) for the assessment of myocardial scar and ECG-gated cine CMR for global cardiac function. TPM was acquired in short-axis orientation at base, mid, and apex for all subjects. After evaluation by experienced observers, the patients were divided in subgroups according to the presence or absence of LGE (LGE+/LGE-), local wall motion abnormalities (WM+/WM-), and having a preserved (≥50%) or reduced (<50%) ejection fraction (EF+/EF-). TPM data was semi-automatically segmented and global LV torsion was computed for each cardiac time frame for endocardial and epicardial layers, and for the entire myocardium.

RESULTS

Maximum myocardial torsion was significantly lower for patients with reduced EF compared to controls (0.21 ± 0.15°/mm vs. 0.36 ± 0.11°/mm, p = 0.018), but also for patients with wall motion abnormalities (0.21 ± 0.13°/mm vs. 0.36 ± 0.11°/mm, p = 0.004). Global myocardial torsion showed a positive correlation (r = 0.54, p < 0.001) with EF. Moreover, endocardial torsion was significantly higher than epicardial torsion for EF+ subjects (0.56 ± 0.33°/mm vs. 0.34 ± 0.18°/mm, p = 0.039) and for volunteers (0.46 ± 0.16°/mm vs. 0.30 ± 0.09°/mm, p = 0.004). The difference in maximum torsion between endo- and epicardial layers was positively correlated with EF (r = 0.47, p = 0.002) and age (r = 0.37, p = 0.016) for all subjects.

CONCLUSIONS

TPM can be used to detect significant differences in LV torsion in patients with reduced EF and in the presence of local wall motion abnormalities. We were able to quantify torsion differences between the endocardium and epicardium, which vary between patient subgroups and are correlated to age and EF.

Rapid breath-hold assessment of myocardial velocities using spiral UNFOLD-ed SENSE tissue phase mapping

PURPOSE

To develop and validate a rapid breath-hold tissue phase mapping (TPM) sequence.

MATERIALS AND METHODS

The sequence was based on an efficient uniform density spiral acquisition, combined with data acceleration. A novel acquisition and reconstruction strategy enabled combination of UNFOLD (2×) and SENSE (3×): UNFOLD-ed SENSE. The sequence was retrospectively cardiac-gated, and a graphics processing unit (GPU) was used for rapid "online" reconstruction. The optimal UNFOLD parameters for the data were calculated using an in silico model. The technique was validated on a 1.5T MR scanner in 15 patients with known aortic valve disease, against a respiratory self-navigated free-breathing TPM technique. Quantitative image quality measures (velocity-to-noise and edge sharpness) were made as well as calculation of longitudinal, radial, and tangential myocardial velocities in the left ventricle.

RESULTS

The proposed breath-hold TPM data took eight heartbeats to acquire. The breath-hold TPM images had significantly higher edge sharpness (P = 0.0014) than the self-navigated TPM images, but with significantly lower velocity-to-noise ratio (P < 0.0001). There was excellent agreement (r > 0.94) in the longitudinal, radial, and tangential velocities between the self-navigated data and the proposed breath-hold TPM sequence.

CONCLUSION

We demonstrate the feasibility of using spiral UNFOLD-ed SENSE to measure myocardial velocities using a rapid breath-hold spiral TPM sequence. This novel technique might enable accurate measurement of myocardial velocities, in a short scan time, which is especially important in a busy clinical workflow. J. MAGN. RESON. IMAGING 2016;44:1003-1009.

Advanced flow MRI: emerging techniques and applications

Magnetic resonance imaging (MRI) techniques provide non-invasive and non-ionising methods for the highly accurate anatomical depiction of the heart and vessels throughout the cardiac cycle. In addition, the intrinsic sensitivity of MRI to motion offers the unique ability to acquire spatially registered blood flow simultaneously with the morphological data, within a single measurement. In clinical routine, flow MRI is typically accomplished using methods that resolve two spatial dimensions in individual planes and encode the time-resolved velocity in one principal direction, typically oriented perpendicular to the two-dimensional (2D) section. This review describes recently developed advanced MRI flow techniques, which allow for more comprehensive evaluation of blood flow characteristics, such as real-time flow imaging, 2D multiple-venc phase contrast MRI, four-dimensional (4D) flow MRI, quantification of complex haemodynamic properties, and highly accelerated flow imaging. Emerging techniques and novel applications are explored. In addition, applications of these new techniques for the improved evaluation of cardiovascular (aorta, pulmonary arteries, congenital heart disease, atrial fibrillation, coronary arteries) as well as cerebrovascular disease (intra-cranial arteries and veins) are presented.

Myocardial dysfunction in patients with aortic stenosis and hypertensive heart disease assessed by MR tissue phase mapping

PURPOSE

To identify abnormalities of myocardial velocities in patients with left ventricular pressure overload using magnetic resonance tissue phase mapping (TPM).

MATERIAL AND METHODS

Thirty-three patients (nine with hypertensive heart disease [HYP], 24 with aortic stenosis [AS]) and 41 healthy controls were enrolled. To assess left ventricular motion, a basal, midventricular, and apical slice were acquired using three-directional velocity-encoded phase-contrast MR with a 3T system. Target parameters were peak longitudinal (Vz ) and radial (Vr ) velocity in systole and diastole (Peaksys , Peakdias ). Analysis was done on each myocardial segment. In a subgroup (n = 7 HYP, n = 12 AS, n = 24 controls), measurement was repeated during handgrip exercise.

RESULTS

AS had significantly lower Vz -Peaksys in the inferolateral and inferoseptal wall (P = 0.003-0.029) and Vr -Peaksys in the septum and anterior wall (P = 0.001-0.013) than controls. Vz -Peakdias and Vr -Peakdias were lower in AS than in controls in almost all segments (P < 0.001-0.028). HYP showed reduced Vz -Peakdias compared to controls in all basal segments as well as in the lateral midventricular wall (P < 0.001-0.045), and reduced Vr -Peakdias compared to controls predominantly in the midventricular and apical segments (P = 0.004-0.042). AS patients with focal fibrosis had significantly reduced myocardial velocities (P = 0.001-0.047) in segments without late enhancement. During exercise, Vz -Peaksys , Vr -Peaksys , and Vz -Peakdias remained unchanged in AS and HYP, but decreased in the lateral wall in controls (P < 0.001-0.043).

CONCLUSION

Even with preserved left ventricle (LV) ejection fraction, peak longitudinal and radial velocities of the LV are reduced in AS and HYP, indicating early functional impairment. J. Magn. Reson. Imaging 2016;44:168-177.

Accelerated three-dimensional cine phase contrast imaging using randomly undersampled echo planar imaging with compressed sensing reconstruction

The aim of this study was to implement and evaluate an accelerated three-dimensional (3D) cine phase contrast MRI sequence by combining a randomly sampled 3D k-space acquisition sequence with an echo planar imaging (EPI) readout. An accelerated 3D cine phase contrast MRI sequence was implemented by combining EPI readout with randomly undersampled 3D k-space data suitable for compressed sensing (CS) reconstruction. The undersampled data were then reconstructed using low-dimensional structural self-learning and thresholding (LOST). 3D phase contrast MRI was acquired in 11 healthy adults using an overall acceleration of 7 (EPI factor of 3 and CS rate of 3). For comparison, a single two-dimensional (2D) cine phase contrast scan was also performed with sensitivity encoding (SENSE) rate 2 and approximately at the level of the pulmonary artery bifurcation. The stroke volume and mean velocity in both the ascending and descending aorta were measured and compared between two sequences using Bland-Altman plots. An average scan time of 3 min and 30 s, corresponding to an acceleration rate of 7, was achieved for 3D cine phase contrast scan with one direction flow encoding, voxel size of 2 × 2 × 3 mm(3) , foot-head coverage of 6 cm and temporal resolution of 30 ms. The mean velocity and stroke volume in both the ascending and descending aorta were statistically equivalent between the proposed 3D sequence and the standard 2D cine phase contrast sequence. The combination of EPI with a randomly undersampled 3D k-space sampling sequence using LOST reconstruction allows a seven-fold reduction in scan time of 3D cine phase contrast MRI without compromising blood flow quantification.

Longitudinal myocardial peak velocities using high temporal resolution phase-contrast and simple averaging are comparable to tissue Doppler echocardiography

OBJECT

Phase contrast imaging is widely used to measure blood velocity. However tissue Doppler imaging (TDI) echocardiography is the reference for myocardial velocity assessment. This study aims at validating the ability of phase contrast (PC) sequences to correctly assess myocardial velocities and to compare these velocities to TDI. The phase contrast sequence was performed with breath-hold parameters and with parameters tuned to increase temporal resolution in free breathing.

MATERIALS AND METHODS

Left and Right auriculo-ventricular annuluses longitudinal velocities were recorded on six healthy volunteers with different temporal resolutions (TDI: 5 ms, breath-hold PC: 94 ms and free-breathing PC: 19 ms). Free-breathing PC was obtained by averaging of three excitations. Amplitudes of four standard echocardiographic and clinically relevant myocardial longitudinal velocity waves were compared: Early filling and auricular, systolic and isovolumic contractions.

RESULTS

Isovolumic contraction waves were only visible with free-breathing PC and TDI. The differences with the reference TDI wave velocities were lower (p = 0.02) for free-breathing PC (19.2 ± 2.6%) than for breath-hold PC (28.1 ± 2.9%). These differences for free-breathing PC were close to (p = 0.21) the coefficient of variation of the measurements provided by TDI (14.8 ± 1.2%).

CONCLUSION

Myocardial longitudinal peak velocities can be assessed with a PC sequence tuned to optimize temporal resolution.