The relationship between myocardial microstructure and strain in chronic infarction using cardiovascular magnetic resonance diffusion tensor imaging and feature tracking

BACKGROUND

Cardiac diffusion tensor imaging (cDTI) using cardiovascular magnetic resonance (CMR) is a novel technique for the non-invasive assessment of myocardial microstructure. Previous studies have shown myocardial infarction to result in loss of sheetlet angularity, derived by reduced secondary eigenvector (E2A) and reduction in subendocardial cardiomyocytes, evidenced by loss of myocytes with right-handed orientation (RHM) on helix angle (HA) maps. Myocardial strain assessed using feature tracking-CMR (FT-CMR) is a sensitive marker of sub-clinical myocardial dysfunction. We sought to explore the relationship between these two techniques (strain and cDTI) in patients at 3 months following ST-elevation MI (STEMI).

METHODS

32 patients (F = 28, 60 ± 10 years) underwent 3T CMR three months after STEMI (mean interval 105 ± 17 days) with second order motion compensated (M2), free-breathing spin echo cDTI, cine gradient echo and late gadolinium enhancement (LGE) imaging. HA maps divided into left-handed HA (LHM, - 90 < HA < - 30), circumferential HA (CM, - 30° < HA < 30°), and right-handed HA (RHM, 30° < HA < 90°) were reported as relative proportions. Global and segmental analysis was undertaken.

RESULTS

Mean left ventricular ejection fraction (LVEF) was 44 ± 10% with a mean infarct size of 18 ± 12 g and a mean infarct segment LGE enhancement of 66 ± 21%. Mean global radial strain was 19 ± 6, mean global circumferential strain was - 13 ± - 3 and mean global longitudinal strain was - 10 ± - 3. Global and segmental radial strain correlated significantly with E2A in infarcted segments (p = 0.002, p = 0.011). Both global and segmental longitudinal strain correlated with RHM of infarcted segments on HA maps (p < 0.001, p = 0.003). Mean Diffusivity (MD) correlated significantly with the global infarct size (p < 0.008). When patients were categorised according to LVEF (reduced, mid-range and preserved), all cDTI parameters differed significantly between the three groups.

CONCLUSION

Change in sheetlet orientation assessed using E2A from cDTI correlates with impaired radial strain. Segments with fewer subendocardial cardiomyocytes, evidenced by a lower proportion of myocytes with right-handed orientation on HA maps, show impaired longitudinal strain. Infarct segment enhancement correlates significantly with E2A and RHM. Our data has demonstrated a link between myocardial microstructure and contractility following myocardial infarction, suggesting a potential role for CMR cDTI to clinically relevant functional impact.

Regional analysis of inflammation and contractile function in reperfused acute myocardial infarction by in vivo 19F cardiovascular magnetic resonance in pigs

Inflammatory cell infiltration is central to healing after acute myocardial infarction (AMI). The relation of regional inflammation to edema, infarct size (IS), microvascular obstruction (MVO), intramyocardial hemorrhage (IMH), and regional and global LV function is not clear. Here we noninvasively characterized regional inflammation and contractile function in reperfused AMI in pigs using fluorine (¹⁹F) cardiovascular magnetic resonance (CMR). Adult anesthetized pigs underwent left anterior descending coronary artery instrumentation with either 90 min occlusion (n = 17) or without occlusion (sham, n = 5). After 3 days, in surviving animals a perfluorooctyl bromide nanoemulsion was infused intravenously to label monocytes/macrophages. At day 6, in vivo ¹H-CMR was performed with cine, T2 and T2* weighted imaging, T2 and T1 mapping, perfusion and late gadolinium enhancement followed by ¹⁹F-CMR. Pigs were sacrificed for subsequent ex vivo scans and histology. Edema extent was 35 ± 8% and IS was 22 ± 6% of LV mass. Six of ten surviving AMI animals displayed both MVO and IMH (3.3 ± 1.6% and 1.9 ± 0.8% of LV mass). The ¹⁹F signal, reflecting the presence and density of monocytes/macrophages, was consistently smaller than edema volume or IS and not apparent in remote areas. The ¹⁹F signal-to-noise ratio (SNR) > 8 in the infarct border zone was associated with impaired remote systolic wall thickening. A whole heart value of ¹⁹F integral (¹⁹F SNR × milliliter) > 200 was related to initial LV remodeling independently of edema, IS, MVO, and IMH. Thus, ¹⁹F-CMR quantitatively characterizes regional inflammation after AMI and its relation to edema, IS, MVO, IMH and regional and global LV function and remodeling.

Performance of a machine-learning algorithm for fully automatic LGE scar quantification in the large multi-national derivate registry

Funding Acknowledgements

Type of funding sources: Other. Main funding source(s): J. Schwitter receives research support by “ Bayer Schweiz AG “. C.N.C. received grant by Siemens. Gianluca Pontone received institutional fees by General Electric, Bracco, Heartflow, Medtronic, and Bayer. U.J.S received grand by Astellas, Bayer, General Electric. This work was supported by Italian Ministry of Health, Rome, Italy (RC 2017 R659/17-CCM698). This work was supported by Gyrotools, Zurich, Switzerland.

Background

 Late Gadolinium enhancement (LGE) scar quantification is generally recognized as an accurate and reproducible technique, but it is observer-dependent and time consuming. Machine learning (ML) potentially offers to solve this problem. 

Purpose

 to develop and validate a ML-algorithm to allow for scar quantification thereby fully avoiding observer variability, and to apply this algorithm to the prospective international multicentre Derivate cohort.

Method

 The Derivate Registry collected heart failure patients with LV ejection fraction &lt;50% in 20 European and US centres. In the post-myocardial infarction patients (n = 689) quality of the LGE short-axis breath-hold images was determined (good, acceptable, sufficient, borderline, poor, excluded) and ground truth (GT) was produced (endo-epicardial contours, 2 remote reference regions, artefact elimination) to determine mass of non-infarcted myocardium and of dense (≥5SD above mean-remote) and non-dense scar (&gt;2SD to &lt;5SD above mean-remote). Data were divided into the learning (total n = 573; training: n = 289; testing: n = 284) and validation set (n = 116). A Ternaus-network (loss function = average of dice and binary-cross-entropy) produced 4 outputs (initial prediction, test time augmentation (TTA), threshold-based prediction (TB), and TTA + TB) representing normal myocardium, non-dense, and dense scar (Figure 1).Outputs were evaluated by dice metrics, Bland-Altman, and correlations. 

Results

 In the validation and test data sets, both not used for training, the dense scar GT was 20.8 ± 9.6% and 21.9 ± 13.3% of LV mass, respectively. The TTA-network yielded the best results with small biases vs GT (-2.2 ± 6.1%, p &lt; 0.02; -1.7 ± 6.0%, p &lt; 0.003, respectively) and 95%CI vs GT in the range of inter-human comparisons, i.e. TTA yielded SD of the differences vs GT in the validation and test data of 6.1 and 6.0 percentage points (%p), respectively (Fig 2), which was comparable to the 7.7%p for the inter-observer comparison (n = 40). For non-dense scar, TTA performance was similar with small biases (-1.9 ± 8.6%, p &lt; 0.0005, -1.4 ± 8.2%, p &lt; 0.0001, in the validation and test sets, respectively, GT 39.2 ± 13.8% and 42.1 ± 14.2%) and acceptable 95%CI with SD of the differences of 8.6 and 8.2%p for TTA vs GT, respectively, and 9.3%p for inter-observer. 

Conclusions

 In the large Derivate cohort from 20 centres, performance of the presented ML-algorithm to quantify dense and non-dense scar fully automatically is comparable to that of experienced humans with small bias and acceptable 95%-CI. Such a tool could facilitate scar quantification in clinical routine as it eliminates human observer variability and can handle large data sets.

Diffuse myocardial fibrosis precedes impairment of myocardial strain in patients with systemic sclerosis

Funding Acknowledgements

Type of funding sources: None.

Background - Myocardial involvement is common in patients with systemic sclerosis (SSc) and causes myocardial fibrosis and subtle ventricular dysfunction. However, the temporal onset of myocardial involvement during the progression of the disease is yet unknown.

Purpose - To investigated the presence of subclinical functional impairment and diffuse myocardial fibrosis in patients with very early diagnosis of SSc (VEDOSS) and to compared the findings to patients with established SSc and healthy controls.

Methods - 110 SSc patients (86 with established SSc and 24 with VEDOSS) and 15 healthy controls were prospectively recruited. The study subjects underwent cardiovascular magnetic resonance on a clinical 1.5T system. Pre- and post-contrast T1 mapping was performed using a MOLLI (Modified Look-Locker Inversion Recovery) sequence. For extracellular volume (ECV) measurements, a single bolus protocol with image acquisition 15-20 min. post-contrast injection was used. For the assessment of subtle functional impairment, global longitudinal (GLS) and circumferential (GCS) myocardial strain were evaluated.

Results - Native T1 values and ECV were elevated in VEDOSS and in patients with established SSc compared to controls (p &lt; 0.001; Figure 1 A & B). GLS was similar in VEDOSS and controls but significantly reduced in patients with established SSc (p &lt; 0.001; Figure 1 C). GCS was similar over all groups (p = 0.88). Patients with clinical evidence of pulmonary or gastrointestinal involvement had higher ECV or T1 values, respectively. Patients with clinical signs of cardiac involvement had lower absolute GLS. SSc subtype, classification or disease duration were not associated with the extent of myocardial fibrosis or impaired strain.

Conclusion - Subclinical myocardial involvement first manifests as diffuse myocardial fibrosis identified by expansion of ECV and increased native T1 in VEDOSS patients while subtle functional impairment as measured by GLS only occurs in established SSc. No single clinical feature of SSc shows a strong association with subtle myocardial involvement.

Microstructural characteristics of chronic infarct segments assessed using diffusion tensor imaging

Funding Acknowledgements

Type of funding sources: Foundation. Main funding source(s): British Heart Foundation

Background

The microstructural changes following myocardial infarction (MI) can be characterised in-vivo with cardiac diffusion tensor imaging (cDTI) imaging, using mean diffusivity (MD), fractional anisotropy (FA), secondary eigenvector angle (E2A) and helix angle (HA) maps. In this study, we use cDTI to explore the microstructural differences between subendocardial and transmural chronic infarct segments.

Method

Twenty STEMI patients (15 men, 5 women, mean age 59) underwent 3T CMR scan at 3 months following presentation (mean interval 107 ± 18 days). Scan protocol included: second order motion compensated (M012) free-breathing spin echo DTI (3 slices, 18 diffusion directions at b-values 100s/mm2[3], 200s/mm2[3] and 500s/mm2[12], acquired resolution was 2.20x2.27x8mm3; cine gradient echo and LGE imaging. Average MD, FA, E2A and HA parameters were calculated on a 16-AHA-segmental level. HA maps were described by dividing values into left-handed HA (LHM, -90° &lt; HA &lt; -30°), circumferential HA (CM, -30° &lt; HA &lt; 30°), and right-handed HA (RHM, 30° &lt; HA &lt; 90°) and reported as relative proportions. Infarct segments were identified using LGE; patients were categorised according to the maximal transmurality of their infarct segments, into subendocardial (&lt;50% LGE) or transmural (&gt;50% LGE) MI.

Results

DTI acquisition was successful in all patients (acquisition time 13 ± 5mins). Ten patients had transmural MI. The results are shown in table 1. Transmurally infarcted segments had significantly lower FA (FA subendocardial MI = 0.27 ± 0.04, FA transmural MI = 0.23 ± 0.02, p &lt; 0.01), lower E2A (E2A subendocardial MI = 47 ± 7°, E2A transmural MI = 38 ± 6°, p &lt; 0.01) and lower proportions of right-handed cardiomyocytes (RHM subendocardial MI = 21 ± 5%, RHM transmural MI = 14 ± 5%, p &lt; 0.01) than subendocardial infarct segments. 

Conclusion

Compared to subendocardial MI segments, the diffusion of water molecules is more isotropic in transmurally infarcted myocardium as evidenced by lower FA values, signifying increased structural disarray. The significantly lower E2A values suggest that laminar sheetlets of transmural infarct segments remain fixed at shallower angles during systole and are unable to reach their usual contractile configuration. The lower proportions of RHM on HA maps highlight the significantly greater loss of subendocardial cardiomyocytes in transmural infarct segments. Further studies are required to assess if these segmental changes can be predictive of long-term LV remodelling.

The relationship between myocardial microstructure and strain in chronic infarcts, assessed using diffusion tensor imaging and feature tracking

Funding Acknowledgements

Type of funding sources: Foundation. Main funding source(s): British Heart Foundation

Background

Cardiac diffusion tensor imaging (cDTI) is a novel technique for the non-invasive assessment of myocardial microstructure.  It allows in-vivo characterisation of microstructural changes post myocardial infarction (MI). Previously published evidence shows significant loss of sheetlet orientation as derived by cDTI secondary eigenvector (E2A), and loss of subendocardial cardiomyocytes derived by reductions in the proportions of myocytes with right-handed orientation (RHM) on helix angle (HA) maps. The assessment of myocardial strain by feature tracking (FT) allows the measurement of radial strain (RS), thought to be driven by the dynamic reorientation of laminar sheetlets, and longitudinal strain (LS), which is thought to relate to subendocardial function. We sought to explore the relationship between the strain and cDTI parameters in patients at 3 months following ST-elevation MI (STEMI).

Methods

Twenty five STEMI patients (M:F = 18:7, mean age 58 ± 9) underwent 3T CMR scan (mean interval 106 ± 17 days) with the following protocol: second order motion compensated (M2), free-breathing spin echo DTI (3 slices, 18 diffusion directions at b-values 100s/mm2, 200s/mm2 and 500s/mm2, acquired resolution was 2.20*2.27*8mm3; cine gradient echo and Late Gadolinium Enhancement (LGE) imaging. HA maps were described by dividing values into left-handed HA (LHM, -90&lt; HA &lt; -30), circumferential HA (CM, -30° &lt; HA &lt; 30°), and right-handed HA (RHM, 30° &lt; HA &lt; 90°) and reported as relative proportions. Segmental analysis were undertaken to derive: HA proportions, E2A, longitudinal strain and LGE%. Segments positive for LGE were classed as infarct segments.

Results

cDTI acquisition was successful in all patients (acquisition time 13 ± 5mins). Mean ejection fraction was 47 ± 8% with mean LGE in the infarcted segment of 57 ± 27%. Mean radial strain was 21 (95% confidence interval, 15-26). The mean E2A was 44 (95% confidence interval 41-47). There was a significant correlation between segmental radial strain and segmental E2A in infarcted segments (p &lt; 0.001, figure 1). In addition, segmental longitudinal strain correlated with the proportion of RHM on HA maps (p &lt; 0.02, figure 2).

Conclusion

Through the combined use of cDTI and FT in patients with chronic infarcts, our results show that the loss of sheetlet orientation assessed using E2A, correlates with worsening radial strain. Segments with less subendocardial cardiomyocytes, evidenced by a lower proportion of myocytes with right-handed orientation on HA maps, correlated with worse longitudinal strain. While this could potentially elucidate the complex association between myocardial microstructure and regional function, further studies are needed to define the incremental clinical value of cDTI.

Segmental strain analysis for the detection of chronic ischemic scars in non-contrast cardiac MRI cine images

Cardiac magnetic resonance imaging (MRI) with late gadolinium enhancement (LGE) is considered the gold standard for scar detection after myocardial infarction. In times of increasing skepticism about gadolinium depositions in brain tissue and contraindications of gadolinium administration in some patient groups, tissue strain-based techniques for detecting ischemic scars should be further developed as part of clinical protocols. Therefore, the objective of the present work was to investigate whether segmental strain is noticeably affected in chronic infarcts and thus can be potentially used for infarct detection based on routinely acquired non-contrast cine images in patients with known coronary artery disease (CAD). Forty-six patients with known CAD and chronic scars in LGE images (5 female, mean age 52 ± 19 years) and 24 gender- and age-matched controls with normal cardiac MRI (2 female, mean age 47 ± 13 years) were retrospectively enrolled. Global (global peak circumferential [GPCS], global peak longitudinal [GPLS], global peak radial strain [GPRS]) and segmental (segmental peak circumferential [SPCS], segmental peak longitudinal [SPLS], segmental peak radial strain [SPRS]) strain parameters were calculated from standard non-contrast balanced SSFP cine sequences using commercially available software (Segment CMR, Medviso, Sweden). Visual wall motion assessment of short axis cine images as well as segmental circumferential strain calculations (endo-/epicardially contoured short axis cine and resulting polar plot strain map) of every patient and control were presented in random order to two independent blinded readers, which should localize potentially infarcted segments in those datasets blinded to LGE images and patient information. Global strain values were impaired in patients compared to controls (GPCS p = 0.02; GPLS p = 0.04; GPRS p = 0.01). Patients with preserved ejection fraction showed also impeded GPCS compared to healthy individuals (p = 0.04). In patients, mean SPCS was significantly impaired in subendocardially (-  5.4% ± 2) and in transmurally infarcted segments (- 1.2% ± 3) compared to remote myocardium (- 12.9% ± 3, p = 0.02 and 0.03, respectively). ROC analysis revealed an optimal cut-off value for SPCS for discriminating infarcted from remote myocardium of - 7.2% with a sensitivity of 89.4% and specificity of 85.7%. Mean SPRS was impeded in transmurally infarcted segments (15.9% ± 6) compared to SPRS of remote myocardium (31.4% ± 5; p = 0.02). The optimal cut-off value for SPRS for discriminating scar tissue from remote myocardium was 16.6% with a sensitivity of 83.3% and specificity of 76.5%. 80.3% of all in LGE infarcted segments (118/147) were correctly localized in segmental circumferential strain calculations based on non-contrast cine images compared to 53.7% (79/147) of infarcted segments detected by visual wall motion assessment (p > 0.01). Global strain parameters are impaired in patients with chronic infarcts compared to controls. Mean SPCS and SPRS in scar tissue is impeded compared to remote myocardium in infarcts patients. Blinded to LGE images, two readers correctly localized 80% of infarcted segments in segmental circumferential strain calculations based on non-contrast cine images, in contrast to only 54% of infarcted segments detected due to wall motion abnormalities in visual wall motion assessment. Analysis of segmental circumferential strain shows a promising method for detection of chronic scars in routinely acquired, non-contrast cine images for patients who cannot receive or decline gadolinium.

The effect of microvascular obstruction on the myocardial microstructure: a diffusion tensor imaging study

Funding Acknowledgements

Type of funding sources: Foundation. Main funding source(s): Heart Research UK

Background

Diffusion tensor cardiac magnetic resonance (DT-CMR) imaging allows for characterising myocardial microstructure in-vivo using mean diffusivity (MD), fractional anisotropy (FA), secondary eigenvector angle (E2A) and helix angle (HA) maps. Following myocardial infarction (MI), alterations in MD, FA and HA proportions have previously been reported. E2A depicts the contractile state of myocardial sheetlets, however the behaviour of E2A in infarct segments, and all DTI markers in areas of microvascular obstruction (MVO) is also not fully understood. 

Purpose

We performed spin echo DTI in patients following ST-elevation MI (STEMI) in order to investigate acute changes in DTI parameters in remote and infarct segments both with and without MVO.

Method

Twenty STEMI patients (16 men, 4 women, mean age 59) had acute (5 ± 2d) 3T CMR scans. CMR protocol included: second order motion compensated (M012) free-breathing spin echo DTI (3 slices, 18 diffusion directions at b-values 100s/mm2[3], 200s/mm2[3] and 500s/mm2[12], reconstructed resolution was 1.66x1.66x8mm); cine and late gadolinium enhancement (LGE) imaging. Average MD, FA, E2A HA parameters were calculated on a  16 AHA segmental level. HA maps were described by dividing values into left-handed HA (LHM, -90° &lt; HA &lt; -30°), circumferential HA (CM, -30° &lt; HA &lt; 30°), and right-handed HA (RHM, 30° &lt; HA &lt; 90°) and reported as relative proportions. Segments were defined as infarct (positive for LGE) and remote (opposite to the infarct). 

Results

DTI acquisition was successful in all patients (acquisition time 13 ± 5mins). Ten patients had evidence of MVO on LGE images. MD was significantly higher in infarct regions in comparison to remote; MVO-ve infarct segments had significantly higher MD than MVO + ve infarct segments (MD remote= 1.46 ± 0.12x10-3mm2/s, MD MVO + ve = 1.59 ± 0.12x10-3mm2/s, MD MVO-ve  = 1.75 ± 0.12x10-3mm2/s, ANOVA p &lt; 0.01). FA was reduced in infarct segments in comparison to remote; MVO-ve infarct segments had significantly lower FA than MVO + ve infarct segments (FAremote= 0.37 ± 0.02, FA MVO + ve = 0.31 ± 0.02 x 10-3mm2/s, MD MVO-ve =0.25 ± 0.02, ANOVA p &lt; 0.01).

E2A values were significantly lower in infarct segments compared to remote; MVO + ve infarct segments had significantly lower values than MVO-ve. (E2A remote= 57.4 ± 5.2°, E2A MVO-ve = 46.8 ± 2.5°, E2A MVO + ve = 36.8 ± 3.1°, ANOVA p &lt; 0.001). RHM% (corresponding to subendocardium) was significantly lower in infarct segments compared to remote; MVO + ve infarct segments had significantly lower RHM% than MVO-ve. (RHM remote= 37 ± 3%, RHM RHM MVO-ve= 28 ± 7%, MVO + ve= 8 ± 5%, ANOVA p &lt; 0.001).

Conclusion

The presence of MVO results in a decrease in MD and increase in FA in comparison to surrounding infarct segments. However, the reduction in E2A and right-handed myocytes on HA in infarct segments is further exacerbated by the presence of MVO. Further study is required to investigate the underlying mechanisms for such alterations in signal intensity.

Abstract Figure. A case of transmural septal MI with MVO

Age-dependent cardiac remodelling – role of sex hormones

Background

While cardiovascular mortality in women has exceeded those in men, women continue to be underrepresented in cardiovascular clinical trials. Further, preclinical experiments are predominantly conducted in male animals, rendering sex-specific variables contributing to cardiovascular disease largely unknown. As age and menopause remain to be key risk factors for cardiovascular disease in women, the aim of this study was to identify key variables of cardiac remodelling in the aging female and male heart, as well as to assess effects of sex hormone deprivation on left ventricular (LV) morphology, LV function and cardiac sympathetic activity.

Materials and methods

Gonadectomized and sham-operated FVB/N mice of both sexes were subjected to positron emission tomography (PET) and cardiac magnetic resonance (CMR) imaging at the age of 4 (young cohort) and 20 (aged cohort) months (total n=123, 55% females). Following tail-vein injection of [11C]meta-hydroxynorephedrine ([11C]mHED), a widely used PET probe in preclinical and clinical assessment of cardiac sympathetic integrity, animals were scanned and cardiac sympathetic outflow was derived from myocardial [11C]mHED uptake. Cardiac parameters including LV volumes and left ventricular ejection fraction (LVEF) were obtained from electrocardiogram (ECG)-gated CMR imaging.

Results and discussion

A significant increase of LVEF was observed in aging females (p=0.012, Figure 1), but not in males. The latter was not associated with a higher cardiac output, and was a consequence of reduced LV end-systolic volumes (p=0.008), unveiling a substantial reduction of size in the aging female heart. As this age-dependent observation was not present in gonadectomized animals (p=0.414), the lack of growth-stimulating estrogen might account for reduction of cardiac size in aging females. Thus, despite a significantly heightened body weight, female heart size is reduced with age. Accordingly, sufficient cardiac output was maintained via increased heart rate (p=0.005) and cardiac sympathetic activity (p=0.040, Figure 1). Gonadectomy accelerated age-dependent changes in LV morphology and function in female mice. While sex hormone deprivation blunted cardiac sympathetic activity and norepinephrine levels in male mice, an opposite trend was observed in females.

Conclusion

Despite increasing body weight with age, aged female and male hearts maintain a stable circulatory blood supply, however, by distinct mechanisms. While the “shrinking” female heart requires an increased heart rate and cardiac sympathetic activity to compensate for smaller ventricular volumes, aging males maintain cardiac size. Importantly, sex hormone deprivation at a young age accelerates age-dependent changes in LV morphology and function in female mice, but not in male mice. The increased sympathetic activity reflects a higher stress level in aged females that might expose them to a higher cardiac vulnerability at postmenopausal age.

Funding Acknowledgement

Type of funding source: Public grant(s) – National budget only. Main funding source(s): Swiss National Science Foundation; Swissheart Foundation

Spatial Resolution and the Magnitude of Infarct Volume Measurement Error in DWI in Acute Ischemic Stroke

BACKGROUND AND PURPOSE

Infarct volume in acute ischemic stroke is an important prognostic marker and determines endovascular treatment decisions. This study evaluates the magnitude and potential clinical impact of the error related to partial volume effects in infarct volume measurement on diffusion-weighted MR imaging in acute stroke and explores how increasing spatial resolution could reduce this error.

MATERIALS AND METHODS

Diffusion-weighted imaging of 393 patients with acute stroke, of whom 56 had anterior circulation large-vessel occlusion, was coregistered to standard space. Lesion boundaries were manually segmented. A 3D lesion-volume model was resampled for voxel sizes from 4 × 4 × 8 to 1 × 1 × 2 mm, and the surface-volume, corresponding to the partial volume error, was calculated. The number of cases with anterior circulation large-vessel occlusion, in which the endovascular therapy core threshold of 70 mL was contained within the margin of error, was calculated as a function of imaging resolution.

RESULTS

The mean infarct core volume was 27.2 ± 49.9 mL. The mean surface volume was 14.7 ± 20.8 mL for 2 × 2 × 4 mm resolution and 7.4 ± 10.7 mL for 1 × 1 × 2 mm resolution. With a resolution of 2 × 2 × 4 mm, 70 mL was contained within the margin of error in 7/56 cases (12.5%) with large-vessel occlusion, while with a 1 × 1 × 2 mm voxel size, the margin of error was 3/56 (5%). The lesion-volume range of potentially misclassified lesions dropped from 46.5-94.1 mL for a 2 × 2 × 4 mm resolution to 64.4-80.1 mL for a 1 × 1 × 2 mm resolution.

CONCLUSIONS

Partial volume effect is an important source of error in infarct volume measurement in acute stroke. Increasing spatial resolution substantially decreases the mean error. Standard use of high-resolution DWI should be considered to increase the reliability of infarct volume measurements.

Equilibrated warping: Finite element image registration with finite strain equilibrium gap regularization

In this paper, we propose a novel continuum finite strain formulation of the equilibrium gap regularization for image registration. The equilibrium gap regularization essentially penalizes any deviation from the solution of a hyperelastic body in equilibrium with arbitrary loads prescribed at the boundary. It thus represents a regularization with strong mechanical basis, especially suited for cardiac image analysis. We describe the consistent linearization and discretization of the regularized image registration problem, in the framework of the finite elements method. The method is implemented using FEniCS & VTK, and distributed as a freely available python library. We show that the equilibrated warping method is effective and robust: regularization strength and image noise have minimal impact on motion tracking, especially when compared to strain-based regularization methods such as hyperelastic warping. We also show that equilibrated warping is able to extract main deformation features on both tagged and untagged cardiac magnetic resonance images.

Heterogeneous growth-induced prestrain in the heart

Even when entirely unloaded, biological structures are not stress-free, as shown by Y.C. Fung׳s seminal opening angle experiment on arteries and the left ventricle. As a result of this prestrain, subject-specific geometries extracted from medical imaging do not represent an unloaded reference configuration necessary for mechanical analysis, even if the structure is externally unloaded. Here we propose a new computational method to create physiological residual stress fields in subject-specific left ventricular geometries using the continuum theory of fictitious configurations combined with a fixed-point iteration. We also reproduced the opening angle experiment on four swine models, to characterize the range of normal opening angle values. The proposed method generates residual stress fields which can reliably reproduce the range of opening angles between 8.7±1.8 and 16.6±13.7 as measured experimentally. We demonstrate that including the effects of prestrain reduces the left ventricular stiffness by up to 40%, thus facilitating the ventricular filling, which has a significant impact on cardiac function. This method can improve the fidelity of subject-specific models to improve our understanding of cardiac diseases and to optimize treatment options.

[Epicardial fat: Imaging and implications for diseases of the cardiovascular system]

Since the discovery of the obese (ob) gene product leptin, fat has been considered an endocrine organ. Especially epicardial fat has gained increasing attention in recent years. The epicardial fat plays a major role in fat metabolism; however, harmful properties have also been reported. Echocardiography, computed tomography and cardiac magnetic resonance imaging are the non-invasive tools used to measure epicardial fat volume. This review briefly introduces the basic physiological and pathophysiological considerations concerning epicardial fat. The main issue of this review is the presentation of non-invasive measurement techniques of epicardial fat using various imaging modalities and a literature overview of associations between epicardial fat and common cardiovascular diseases.

Improved UTE-based attenuation correction for cranial PET-MR using dynamic magnetic field monitoring

PURPOSE

Ultrashort echo time (UTE) MRI has been proposed as a way to produce segmented attenuation maps for PET, as it provides contrast between bone, air, and soft tissue. However, UTE sequences require samples to be acquired during rapidly changing gradient fields, which makes the resulting images prone to eddy current artifacts. In this work it is demonstrated that this can lead to misclassification of tissues in segmented attenuation maps (AC maps) and that these effects can be corrected for by measuring the true k-space trajectories using a magnetic field camera.

METHODS

The k-space trajectories during a dual echo UTE sequence were measured using a dynamic magnetic field camera. UTE images were reconstructed using nominal trajectories and again using the measured trajectories. A numerical phantom was used to demonstrate the effect of reconstructing with incorrect trajectories. Images of an ovine leg phantom were reconstructed and segmented and the resulting attenuation maps were compared to a segmented map derived from a CT scan of the same phantom, using the Dice similarity measure. The feasibility of the proposed method was demonstrated in in vivo cranial imaging in five healthy volunteers. Simulated PET data were generated for one volunteer to show the impact of misclassifications on the PET reconstruction.

RESULTS

Images of the numerical phantom exhibited blurring and edge artifacts on the bone-tissue and air-tissue interfaces when nominal k-space trajectories were used, leading to misclassification of soft tissue as bone and misclassification of bone as air. Images of the tissue phantom and the in vivo cranial images exhibited the same artifacts. The artifacts were greatly reduced when the measured trajectories were used. For the tissue phantom, the Dice coefficient for bone in MR relative to CT was 0.616 using the nominal trajectories and 0.814 using the measured trajectories. The Dice coefficients for soft tissue were 0.933 and 0.934 for the nominal and measured cases, respectively. For air the corresponding figures were 0.991 and 0.993. Compared to an unattenuated reference image, the mean error in simulated PET uptake in the brain was 9.16% when AC maps derived from nominal trajectories was used, with errors in the SUV max for simulated lesions in the range of 7.17%-12.19%. Corresponding figures when AC maps derived from measured trajectories were used were 0.34% (mean error) and -0.21% to +1.81% (lesions).

CONCLUSIONS

Eddy current artifacts in UTE imaging can be corrected for by measuring the true k-space trajectories during a calibration scan and using them in subsequent image reconstructions. This improves the accuracy of segmented PET attenuation maps derived from UTE sequences and subsequent PET reconstruction.

Group sparse reconstruction using intensity-based clustering

Compressed sensing has been of great interest to speed up the acquisition of MR images. The k-t group sparse (k-t GS) method has recently been introduced for dynamic MR images to exploit not just the sparsity, as in compressed sensing, but also the spatial group structure in the sparse representation. k-t GS achieves higher acceleration factors compared to the conventional compressed sensing method. However, it assumes a spatial structure in the sparse representation and it requires a time consuming hard-thresholding reconstruction scheme. In this work, we propose to modify k-t GS by incorporating prior information about the sorted intensity of the signal in the sparse representation, for a more general and robust group assignment. This approach is referred to as group sparse reconstruction using intensity-based clustering. The feasibility of the proposed method is demonstrated for static 3D hyperpolarized lung images and applications with both dynamic and intensity changes, such as 2D cine and perfusion cardiac MRI, with retrospective undersampling. For all reported acceleration factors the proposed method outperforms the original compressed sensing method. Improved reconstruction over k-t GS method is demonstrated when k-t GS assumptions are not satisfied. The proposed method was also applied to cardiac cine images with a prospective sevenfold acceleration, outperforming the standard compressed sensing reconstruction.

Toward true 3D visualization of active catheters using compressed sensing

A crucial requirement in MR-guided interventions is the visualization of catheter devices in real time. However, true 3D visualization of the full length of catheters has hitherto been impossible given scan time constraints. Compressed sensing (CS) has recently been proposed as a method to accelerate MR imaging of sparse objects. Images acquired with active interventional devices exhibit a high CNR and are inherently sparse, therefore rendering CS ideally suited for accelerating data acquisition. A framework for true visualization of active catheters in 3D is proposed employing CS to gain high undersampling factors making real-time applications feasible. Constraints are introduced taking into account prior knowledge of catheter geometry and catheter motion over time to improve and accelerate image reconstruction. The potential of the method is demonstrated using computer simulations and phantom experiments and in vivo feasibility is demonstrated in a pig experiment.

High wall shear stress measured by magnetic resonance is a predictor of restenosis in the femoral-artery after balloon angioplasty

AIM

Wall shear stress (WSS) has been implied in the pathogenesis of restenosis after percutaneous transluminal angioplasty (PTA). Aims of the present study were to calculate WSS in the superficial femoral artery (SFA) from magnetic resonance imaging (MRI) and from duplex sonography in healthy controls and in patients after PTA of the SFA to evaluate the predictive value of WSS for restenosis.

METHODS

WSS was assessed by calculating the slope of velocity profiles at the vessel wall from data obtained with velocity encoded cine MR and with duplex using the formula: Shear stress=4xblood viscosityxpeak blood velocity/internal diameter. Seventeen patients were studied 1 day after successful PTA of the SFA. Restenosis was determined by duplex ultrasound at the 6-months follow-up visit.

RESULTS

In healthy controls WSS values calculated from MRI and from duplex were similar (1.86+/-0.35 N/m2 vs 1.88+/-0.34 N/m2, n.s.). In patients the values obtained with duplex were higher than those obtained with MRI (4.1+/-2.3 N/m2 vs 2.4+/-1.2 N/m2, p=0.002). With both methods post-interventional WSS was higher in patients developing restenosis (duplex 5.4+/-2.2 N/m2, MRI 3.1+/-0.9 N/m2) than in those without restenosis (duplex 2.7+/-1.4 N/m2, MRI 1.5+/-0.7 N/m2) and was revealed to be an independent predictor of restenosis (p=0.03).

CONCLUSIONS

This is the first study demonstrating that increased post-interventional WSS in the SFA is predictive for restenosis. WSS values obtained with MRI and duplex were different in patients, however with both methods higher WSS was associated with restenosis.

Coronary flow quantification by Fourier velocity encoded MRI

Recently, the feasibility of measuring coronary blood flow using fast magnetic resonance (MR) techniques was reported. Thus, MR holds potential to non-invasively assess significance of coronary stenosis. However, the accurate determination of flow and vessel area still remains challenging. High spatial and temporal resolution is required to assess reliable flow profiles within the coronary arteries. For this purpose, Fourier velocity encoding (FVE) was implemented with a small number of encoding steps. Simulations and in-vitro experiments have been performed to demonstrate the benefit of FVE for flow quantification. Further, initial volunteer measurements have shown its potential for invivo application.

Aortic and mitral regurgitation: quantification using moving slice velocity mapping

Comprehensive assessment of the severity of valvular insufficiency includes quantification of regurgitant volumes. Previous methods lack reliable slice positioning with respect to the valve and are prone to velocity offsets due to through-plane motion of the valvular plane of the heart. Recently, the moving slice velocity mapping technique was proposed. In this study, the technique was applied for quantification of mitral and aortic regurgitation. Time-efficient navigator-based respiratory artifact suppression was achieved by implementing a prospective k-space reordering scheme in conjunction with slice position correction. Twelve patients with aortic insufficiency and three patients with mitral insufficiency were studied. Aortic regurgitant volumes were calculated from diastolic velocities mapped with a moving slice 5 mm distal to the aortic valve annulus. Mitral regurgitant flow was indirectly assessed by measuring mitral inflow at the level of the mitral annulus and net aortic outflow. Regurgitant fractions, derived from velocity data corrected for through-plane motion, were compared to data without correction for through-plane motion. In patients with mild and moderate aortic regurgitation, regurgitant fractions differed by 60% and 15%, on average, when comparing corrected and uncorrected data, respectively. Differences in severe aortic regurgitation were less (7%). Due to the large orifice area of the mitral valve, differences were still substantial in moderate-to-severe mitral regurgitation (19%). The moving slice velocity mapping technique was successfully applied in patients with aortic and mitral regurgitation. The importance of correction for valvular through-plane motion is demonstrated.

Visualization of flow patterns distal to aortic valve prostheses in humans using a fast approach for cine 3D velocity mapping

The fluid dynamic performance of mechanical heart valves differs from normal valves and thus is considered related to late clinical complications in patients. Since flow patterns evolving around heart valves are complex in space and time, flow visualization based on time-resolved 3D velocity data might add important information regarding the performance of specific valve designs in vivo. However, previous cine 3D techniques for three-directional phase-contrast velocity mapping suffer from long scan duration and therefore might hamper assessment in patients. A hybrid 3D phase-contrast sequence combining segmented k-space acquisition with short EPI readout trains is presented with its validation in vitro. The technique was applied to study flow patterns downstream from a bileaflet aortic prosthesis in six patients. Navigator echoes were incorporated for respiratory motion compensation. Before flow visualization, spurious phase errors due to concomitant gradient fields and eddy currents were corrected. Flow visualization was based on particle paths and animated velocity vector plots. Dedicated algorithms for particle path integration were implemented to account for the considerable motion of the ascending aorta during the cardiac cycle. A distinct flow pattern reflecting the valve design was observed closest to the valve during early flow acceleration. Reverse flow occurred adjacent to high velocity jets and above the hinge housings. Later in systole, flow became confined to the central vessel area and reverse flow along the inner aortic curvature developed. Further downstream from the valve, flow patterns varied considerably among patients, indicating the impact of varying aortic anatomy in vivo. It is concluded that MR velocity mapping is a potential tool for studying 3D flow patterns evolving around heart valve prostheses in humans. J. Magn. Reson. Imaging 2001;13:690-698.

Heart motion-adapted MR velocity mapping of blood velocity distribution downstream of aortic valve prostheses: initial experience

PURPOSE

To investigate blood flow velocities and shear rates at two distances downstream of an artificial aortic valve in patients.

MATERIALS AND METHODS

Blood velocity was quantified downstream of the valve prosthesis (for replacement after aortic valve stenosis or combined stenosis and regurgitation) in 10 patients by using a magnetic resonance (MR) cine velocity mapping method in which the imaging section position is adapted according to the excursion of the valvular plane of the heart. Two acquisitions were performed to display the blood velocity distributions one-fourth valve diameter and one valve diameter downstream of the valve and to quantify blood volumes and shear rates.

RESULTS

The velocity profiles measured during flow acceleration one-fourth valve diameter downstream were characterized by a distinct pattern of two lateral jets and one central jet of antegrade flow. High shear rates were found along the leaflet tips. The profiles obtained one valve diameter downstream were skewed, with varying velocity patterns among patients. Peak shear rates were found close to the vessel wall. With correction for through-plane motion of the valve, the mean apparent regurgitant fraction (+/- SD) was 14% +/- 6; the mean regurgitant fraction without correction was 9% +/- 5.

CONCLUSION

The described noninvasive procedure for velocity mapping enables measurements close to the valve and thus evaluation of blood flow patterns with respect to valve design in humans.

Heart motion adapted cine phase-contrast flow measurements through the aortic valve

A method for magnetic resonance cine velocity mapping through heart valves with adaptation of both slice offset and angulation according to the motion of the valvular plane of the heart is presented. By means of a subtractive labeling technique, basal myocardial markers are obtained and automatically extracted for quantification of heart motion at the valvular level. The captured excursion of the basal plane is used to calculate the slice offset and angulation of each required time frame for cine velocity mapping. Through-plane velocity offsets are corrected by subtracting velocities introduced by basal plane motion from the measured velocities. For evaluation of the method, flow measurements downstream from the aortic valve were performed both with and without slice adaptation in 11 healthy volunteers and in four patients with aortic regurgitation. Maximum through-plane motion at the aortic root level as calculated from the labeled markers averaged 8.9 mm in the volunteers and 6.5 mm in the patients. The left coronary root was visible in 2-4 (mean: 2.2) time frames during early diastole when imaging with a spatially fixed slice. Time frames obtained with slice adaptation did not contain the coronary roots. Motion correction increased the apparent regurgitant volume by 5.7 +/- 0.4 ml for patients with clinical aortic regurgitation, for an increase of approximately 50%. The proposed method provides flow measurements with correction for through-plane motion perpendicular to the aortic root between the valvular annulus and the coronary ostia throughout the cardiac cycle. Magn Reson Med 42:970-978, 1999.

Automatic vessel segmentation using active contours in cine phase contrast flow measurements

The segmentation of images obtained by cine magnetic resonance (MR) phase contrast velocity mapping using manual or semi-automated methods is a time consuming and observer-dependent process that still hampers the use of flow quantification in a clinical setting. A fully automatic segmentation method based on active contour model algorithms for defining vessel boundaries has been developed. For segmentation, the phase image, in addition to the magnitude image, is used to address image distortions frequently seen in the magnitude image of disturbed flow fields. A modified definition for the active contour model is introduced to reduce the influence of missing or spurious edge information of the vessel wall. The method was evaluated on flow phantom data and on in vivo images acquired in the ascending aorta of humans. Phantom experiments resulted in an error of 0.8% in assessing the luminal area of a flow phantom equipped with an artificial heart valve. Blinded evaluation of the volume flow rates from automatic vs. manual segmentation of gradient echo (FFE) phase contrast images obtained in vivo resulted in a mean difference of -0.9 +/- 3%. The mean difference from automatic vs. manual segmentation of images acquired with a hybrid phase contrast sequence (TFEPI) within a single breath-hold was -0.9 +/- 6%.

Quantitative abdominal aortic flow measurements at controlled levels of ergometer exercise

Measuring the exercise-induced flow changes in the arteries of the body is a major challenge. The use of quantitative MR flow measurements for this purpose is hampered by movement artifacts and ECG triggering problems. To quantify exercise-induced flow changes in the abdominal aorta, we applied a fast hybrid phase contrast sequence with K-space segmentation and echo planar imaging readouts during a 12 heart beat, single breathhold post exercise scanning window after ergometer exercise in nine volunteers. Central k-space was acquired first. The changes in heart rate throughout the scanning window were quantified. The mean decrease in heart rate after six heart beats post exercise was less than 4% and less than 14% after 11 heart beats indicating that the exercise state was very well represented during the acquisition of central k-space. Abdominal aortic flow increased from 1.4+/-0.3 l/min at rest to 7.9+/-1.1 l/min at 131 watt. Retrograde flow reached a maximum value of 1.2 l/min at rest, and lasted 140 ms on average. Only for one out of the nine volunteers was there any retrograde flow present during exercise (at 33 watt and 65 watt exercise). It was concluded that retrograde flow patterns in the abdominal aorta associated with oscillating wall shear stresses and development of atherosclerosis disappeared with increasing levels of exercise. The feasibility of using fast quantitative phase contrast measurements during a post exercise scanning window to represent controlled exercise levels was demonstrated.

Automatic accurate non-invasive quantitation of blood flow, cross-sectional vessel area, and wall shear stress by modelling of magnetic resonance velocity data

OBJECTIVES

To apply a new, automatic and non-invasive method for quantification of blood flow, dynamic cross-sectional vessel area, and wall shear stress (WSS) by in vivo magnetic resonance velocity mapping of normal subjects.

DESIGN

Prospective, open study.

MATERIALS

Six young volunteers.

METHODS

A three-dimensional paraboloid model enabling automatic determination of blood flow, vessel distensibility and WSS was applied to blood velocity determinations in the common carotid artery. Blood flow was also determined by a manual edge detection method.

RESULTS

Using the new method, the common carotid mean blood flow was 7.28 (5.61-9.63) (mean (range)) ml/s. By the manual-method blood flow was 7.21 (5.55-9.60) ml/s. Mean luminal vessel area was 26% larger in peak systole than in diastole. Mean/peak WSS was 0.82/2.28 N/m2. Manually and automatically determined flows correlated (r2 = 0.998, p < 0.0001). WSS and peak centre velocity were associated (r2 = 0.805, p < 0.0001).

CONCLUSIONS

Blood flow, luminal vessel area dilatation, and WSS can be determined by the automatic three-dimensional paraboloid method. The hypothesis of association between peak centre velocity and WSS was not contradicted by the results of the present study.

Quantitation of circumferential subpixel vessel wall position and wall shear stress by multiple sectored three-dimensional paraboloid modeling of velocity encoded cine MR

Methods are lacking for accurate, noninvasive circumferential edge detection and wall shear stress calculation. Using standard MR phase contrast sequences, parts of the velocity profiles were fitted to a multiple sectored three-dimensional paraboloid model enabling exact calculation of vessel wall position and wall shear stress in 24 locations evenly distributed around the luminal vessel wall. The model was evaluated by in vitro scans and computer simulations and applied to the common carotid artery of humans. In vitro, the luminal area of a glass tube was assessed with an error of 0.9%. Computer simulations of peak systolic data revealed errors of +/-0.9% (vessel area) and +/-3.25% (wall shear stress). The in vivo results showed substantial difference between anterior and posterior wall shear stress values due to skewed velocity profiles. A new noninvasive method for highly accurate measurement of circumferential subpixel vessel wall position and wall shear stress has been developed.

Accurate noninvasive quantitation of blood flow, cross-sectional lumen vessel area and wall shear stress by three-dimensional paraboloid modeling of magnetic resonance imaging velocity data

OBJECTIVES

We present a new method in which a priori knowledge of the blood velocity fields within the boundary layer at the vessel wall, combined with acquisition of high resolution magnetic resonance imaging (MRI) blood velocity data, allow exact modeling at the subpixel level.

BACKGROUND

Methods are lacking for accurate, noninvasive estimation of blood flow, dynamic cross-sectional lumen vessel area and wall shear stress.

METHODS

Using standard acquisition of MRI blood flow velocity data, we fitted all data points (n = 69) within the boundary layer of the velocity profile to a three-dimensional paraboloid, which enabled calculation of absolute volume blood flow, circumferential vessel wall position, lumen vessel area and wall shear stress. The method was tested in a 8.00 +/ 0.01-mm diameter glass tube model and applied in vivo to the common carotid artery of seven volunteers.

RESULTS

In vitro the lumen area was assessed with a mean error of 0.6%. The 95% confidence interval included the specified tube dimensions. Common carotid mean blood flow was 7.42 ml/s, and mean (standard error) diastolic/systolic vessel area was 33.25 (0.72 [2.2%])/43.46 (0.65 [1.5%]) mm2. Mean/peak wall shear stress was 0.95 (0.04 [4.2%])/2.56 (0.08 [3.1%]) N/m2.

CONCLUSIONS

We describe a new noninvasive method for highly accurate estimation of blood flow, cross-sectional lumen vessel area and wall shear stress. In vitro results and statistical analysis demonstrate the feasibility of the method, and the first in vivo results are comparable to published data.