Strain analysis in CRT candidates using the novel segment length in cine (SLICE) post-processing technique on standard CMR cine images.

OBJECTIVES: Although myocardial strain analysis is a potential tool to improve patient selection for cardiac resynchronization therapy (CRT), there is currently no validated clinical approach to derive segmental strains. We evaluated the novel segment length in cine (SLICE) technique to derive segmental strains from standard cardiovascular MR (CMR) cine images in CRT candidates. METHODS: Twenty-seven patients with left bundle branch block underwent CMR examination including cine imaging and myocardial tagging (CMR-TAG). SLICE was performed by measuring segment length between anatomical landmarks throughout all phases on short-axis cines. This measure of frame-to-frame segment length change was compared to CMR-TAG circumferential strain measurements. Subsequently, conventional markers of CRT response were calculated. RESULTS: Segmental strains showed good to excellent agreement between SLICE and CMR-TAG (septum strain, intraclass correlation coefficient (ICC) 0.76; lateral wall strain, ICC 0.66). Conventional markers of CRT response also showed close agreement between both methods (ICC 0.61-0.78). Reproducibility of SLICE was excellent for intra-observer testing (all ICC ≥0.76) and good for interobserver testing (all ICC ≥0.61). CONCLUSIONS: The novel SLICE post-processing technique on standard CMR cine images offers both accurate and robust segmental strain measures compared to the 'gold standard' CMR-TAG technique, and has the advantage of being widely available. KEY POINTS: • Myocardial strain analysis could potentially improve patient selection for CRT. • Currently a well validated clinical approach to derive segmental strains is lacking. • The novel SLICE technique derives segmental strains from standard CMR cine images. • SLICE-derived strain markers of CRT response showed close agreement with CMR-TAG. • Future studies will focus on the prognostic value of SLICE in CRT candidates.

Characterization of normal regional myocardial function by MRI cardiac tagging.

PURPOSE: The aim of this study was to establish reference values for segmental myocardial strain measured by magnetic resonance (MR) cardiac tagging in order to characterize the regional function of the heart. MATERIAL AND METHODS: We characterized the left ventricular (LV) systolic deformation in 39 subjects (26 women and 13 men, age 58.8 ± 11.6 years) whose cardiological study had not revealed any significant abnormality. The deformation was measured from MR-tagged (Siemens 1.5T MR) images using an algorithm based on sine wave modeling. Circumferential and radial peak systolic strain values along with the torsion angle and circumferential-longitudinal (CL) shear were determined in 16 LV segments in order to settle the reference values for these parameters. RESULTS: Circumferential strain was highest at the anterior and lateral walls (-20.2 ± 4.0% and -21.8 ± 4.3%, respectively; P < 0.05) and was lowest at the base level (-17.2 ± 3.1% vs. -20.1 ± 3.1% "mid level," P < 0.05; -17.2 ± 3.1% vs. -20.3 ± 3.0% "apical level," P < 0.05). Radial strain highest values were from inferior and lateral walls (13.7 ± 7.4% and 12.8 ± 7.8%, respectively; P < 0.05) and it was lowest medially (9.4 ± 4.1% vs. 13.1 ± 4.1% "base level," P < 0.05; 9.4 ± 4.1% vs. 12.1 ± 4.4% "apical level," P < 0.05). Torsion angle (counterclockwise when viewed from the apex) increased with the distance from the base (7.9 ± 2.4° vs. 16.8 ± 4.4°, P < 0.05), and the highest and lowest values were found at lateral (medial lateral: 12.0 ± 4.4°, apical lateral: 25.1 ± 6.4°, P < 0.05) and septal wall (medial septal: 3.6 ± 2.1°, apical septal: 8.3 ± 5.3°, P < 0.05), respectively. These differences were found again in CL shear values, around the LV circumference. However, CL shear remained constant with increasing distance from the base (9.1 ± 2.6°, medium and 9.8 ± 2.4°, apex). CONCLUSION: In summary, this study provides reference values for the assessment of regional myocardial function by MR cardiac tagging. Comparison of patient deformation parameters with normal deformation patterns may permit early detection of regional systolic dysfunction.

Intertechnique agreement and interstudy reproducibility of strain and diastolic strain rate at 1.5 and 3 Tesla: a comparison of feature-tracking and tagging in patients with aortic stenosis.

PURPOSE: To determine the interstudy reproducibility of myocardial strain and peak early-diastolic strain rate (PEDSR) measurement on cardiovascular magnetic resonance imaging (MRI) assessed with feature tracking (FT) and tagging, in patients with aortic stenosis (AS). MATERIALS AND METHODS: Cardiac MRI was performed twice (1-14 days apart) in 18 patients (8 at 1.5 Tesla [T], 10 at 3T) with moderate-severe AS. Circumferential peak systolic strain (PSS) and PEDSR were measured in all patients. Longitudinal PSS and PEDSR were assessed using FT in all patients, and tagging in the 3T sub-group. RESULTS: PSS was higher with FT than tagging (21.0 ± 1.9% versus 17.0 ± 3.4% at 1.5T, 21.4 ± 4.0% versus 17.7 ± 3.0% at 3T, P < 0.05), as was PEDSR (1.3 ± 0.3 s(-1) versus 1.0 ± 0.3 s(-1) , P = 0.10 at 1.5T and 1.3 ± 0.4 s(-1) versus 0.8 ± 0.3 s(-1) , P < 0.05 at 3T). The reproducibility of PSS was excellent with FT (coefficient of variation [CoV] 9-10%) and good with tagging at 1.5T (13-19%). Reproducibility of circumferential PEDSR was best at 1.5T when only basal/mid slices were included (CoV 12%), but moderate to poor at 3T (29-35%). Reproducibility of longitudinal strain was good with FT (10-16%) but moderate for PEDSR (∼30%). CONCLUSION: In patients with AS, FT consistently produces higher values compared with tagging. The interstudy reproducibility of PSS is excellent with FT and good with tagging. The reproducibility of circumferential PEDSR at 1.5T is good when only basal and mid slices are used.

Assessment of myocardial partition coefficient of gadolinium (λ) in dilated cardiomyopathy and its impact on segmental and global systolic function.

PURPOSE: 1) To assess the myocardial partition coefficient (λ) of gadolinium quantified using T1 mapping in dilated cardiomyopathy (DCM); and 2) to assess the impact of increased λ on left ventricular (LV) circumferential strain and ejection fraction in DCM. MATERIALS AND METHODS: Seventeen patients with DCM and 11 controls were prospectively included. All patients and controls underwent a 1.5 T MRI using: 1) cine to quantify LV volumes and function; 2) tagging to quantify circumferential strain in mid-LV; 3) T1 mapping before and 9 minutes after contrast injection to quantify R1, ΔR1, and λ; and 4) inversion recovery 3D Flash was used to assess late gadolinium enhancement (LGE) 10 minutes after Gd DOTA injection (0.2 mmol/kg). We used Student's t-test to compare means, Pearson's test for correlation assessment, and a mixed linear model to integrate the dependency between myocardial segments. RESULTS: No difference in median λ was measured between patients with (0.52 [interquartile range = 0.48-0.56]) and without enhancement on LGE (0.51 [0.47-0.54]; P = 0.07). Circumferential strain value measured in each segment was correlated with the λ measured in the corresponding segment (r = 0.55; P < 0.0001). Multivariate analysis revealed a significant link between the λ in each segment and circumferential strain (0.002 ± 0.001; P = 0.009) and also with ejection fraction (-0.001 ± 0.0008; P = 0.04). CONCLUSION: In DCM, λ correlates independently with circumferential strain and ejection fraction, suggesting that there is a link between λ and systolic function.

OntoVIP: an ontology for the annotation of object models used for medical image simulation.

This paper describes the creation of a comprehensive conceptualization of object models used in medical image simulation, suitable for major imaging modalities and simulators. The goal is to create an application ontology that can be used to annotate the models in a repository integrated in the Virtual Imaging Platform (VIP), to facilitate their sharing and reuse. Annotations make the anatomical, physiological and pathophysiological content of the object models explicit. In such an interdisciplinary context we chose to rely on a common integration framework provided by a foundational ontology, that facilitates the consistent integration of the various modules extracted from several existing ontologies, i.e. FMA, PATO, MPATH, RadLex and ChEBI. Emphasis is put on methodology for achieving this extraction and integration. The most salient aspects of the ontology are presented, especially the organization in model layers, as well as its use to browse and query the model repository.

Comparison of local sine wave modeling with harmonic phase analysis for the assessment of myocardial strain.

PURPOSE: To compare local sine-wave modeling (SinMod) with harmonic phase analysis (HARP), for assessment of left ventricular (LV) circumferential strain (εcc) from tagged cardiovascular magnetic resonance images. MATERIALS AND METHODS: Mid-ventricular spatial modulation of magnetization was performed in 60 participants (15 each with hypertrophic, dilated or ischemic cardiomyopathy and 15 healthy controls) at 1.5 Tesla. Global and segmental peak transmural εcc were measured using HARP and SinMod. Repeated measurements were performed on 25% of examinations to assess observer variability. Effect of contrast was assessed in 10 additional patients. RESULTS: SinMod showed a high level of agreement with HARP for global εcc (mean difference -0.02, 95% limits of agreement -6.46 to 6.43%). Agreement was much lower for segmental εcc. Both methods showed excellent observer agreement for global εcc (intraclass correlation coefficient >0.75). Observer agreement for segmental εcc was also excellent with SinMod, but was significantly lower with HARP. Analysis time was significantly shorter using SinMod. Pre- and postcontrast εcc measurements were not significantly different using either technique, although postcontrast measurements showed greater variability with HARP. CONCLUSION: SinMod and HARP-based measurements of global εcc have a high level of agreement, but segmental agreement is substantially lower. SinMod has generally lower observer variability, is faster and is less affected by contrast, but requires further validation.

Pixel-wise statistical analysis of myocardial injury in STEMI patients with delayed enhancement MRI.

Objectives: Myocardial injury assessment from delayed enhancement magnetic resonance images is routinely limited to global descriptors such as size and transmurality. Statistical tools from computational anatomy can drastically improve this characterization, and refine the assessment of therapeutic procedures aiming at infarct size reduction. Based on these techniques, we propose a new characterization of myocardial injury up to the pixel resolution. We demonstrate it on the imaging data from the Minimalist Immediate Mechanical Intervention randomized clinical trial (MIMI: NCT01360242), which aimed at comparing immediate and delayed stenting in acute ST-Elevation Myocardial Infarction (STEMI) patients. Methods: We analyzed 123 patients from the MIMI trial (62 ± 12 years, 98 male, 65 immediate 58 delayed stenting). Early and late enhancement images were transported onto a common geometry using techniques inspired by statistical atlases, allowing pixel-wise comparisons across population subgroups. A practical visualization of lesion patterns against specific clinical and therapeutic characteristics was also proposed using state-of-the-art dimensionality reduction. Results: Infarct patterns were roughly comparable between the two treatments across the whole myocardium. Subtle but significant local differences were observed for the LCX and RCA territories with higher transmurality for delayed stenting at lateral and inferior/inferoseptal locations, respectively (15% and 23% of myocardial locations with a p-value <0.05, mainly in these regions). In contrast, global measurements were comparable for all territories (no statistically significant differences for all-except-one measurements before standardization / for all after standardization), although immediate stenting resulted in more subjects without reperfusion injury. Conclusion: Our approach substantially empowers the analysis of lesion patterns with standardized comparisons up to the pixel resolution, and may reveal subtle differences not accessible with global observations. On the MIMI trial data as illustrative case, it confirmed its general conclusions regarding the lack of benefit of delayed stenting, but revealed subgroups differences thanks to the standardized and finer analysis scale.

Automated segmentation of a motion mask to preserve sliding motion in deformable registration of thoracic CT.

PURPOSE: Deformable registration generally relies on the assumption that the sought spatial transformation is smooth. Yet, breathing motion involves sliding of the lung with respect to the chest wall, causing a discontinuity in the motion field, and the smoothness assumption can lead to poor matching accuracy. In response, alternative registration methods have been proposed, several of which rely on prior segmentations. We propose an original method for automatically extracting a particular segmentation, called a motion mask, from a CT image of the thorax. METHODS: The motion mask separates moving from less-moving regions, conveniently allowing simultaneous estimation of their motion, while providing an interface where sliding occurs. The sought segmentation is subanatomical and based on physiological considerations, rather than organ boundaries. We therefore first extract clear anatomical features from the image, with respect to which the mask is defined. Level sets are then used to obtain smooth surfaces interpolating these features. The resulting procedure comes down to a monitored level set segmentation of binary label images. The method was applied to sixteen inhale-exhale image pairs. To illustrate the suitability of the motion masks, they were used during deformable registration of the thorax. RESULTS: For all patients, the obtained motion masks complied with the physiological requirements and were consistent with respect to patient anatomy between inhale and exhale. Registration using the motion mask resulted in higher matching accuracy for all patients, and the improvement was statistically significant. Registration performance was comparable to that obtained using lung masks when considering the entire lung region, but the use of motion masks led to significantly better matching near the diaphragm and mediastinum, for the bony anatomy and for the trachea. The use of the masks was shown to facilitate the registration, allowing to reduce the complexity of the spatial transformation considerably, while maintaining matching accuracy. CONCLUSIONS: We proposed an automated segmentation method for obtaining motion masks, capable of facilitating deformable registration of the thorax. The use of motion masks during registration leads to matching accuracies comparable to the use of lung masks for the lung region but motion masks are more suitable when registering the entire thorax.

Measurement of local orientation of cardiomyocyte aggregates in human left ventricle free wall samples using X-ray phase-contrast microtomography.

Most cardiomyocytes in the left ventricle wall are grouped in aggregates of four to five units that are quasi-parallel to each other. When one or more "cardiomyocyte aggregates" are delimited by two cleavage planes, this defines a "sheetlet" that can be considered as a "work unit" that contributes to the thickening of the wall during the cardiac cycle. In this paper, we introduce the skeleton method to measure the local three-dimensional (3D) orientation of cardiomyocyte aggregates in the sheetlets in three steps: data segmentation; extraction of the skeleton of the sheetlets; and calculation of the local orientation of the cardiomyocyte aggregates inside the sheetlets. These data include a series of virtual tissue volumes and five transmural human left ventricle free wall samples, imaged with 3D synchrotron radiation phase-contrast microtomography, and reconstructed with a 3.5×3.5×3.5µm3 voxel size. We computed the local orientation of the cardiomyocyte aggregates inside the sheetlets with a working window of 112×112×112µm3 in size. These data demonstrate that the skeleton method can provide accurate 3D measurements and reliable screening of the 3D evolution of the orientation of cardiomyocyte aggregates within the sheetlets. We showed that in regions that contain one population of quasi-parallel sheetlets, the orientation of the cardiomyocyte aggregates undergo "oscillations" along the perpendicular direction of the sheetlets. In regions that contain two populations of sheetlets with a different angular range, we demonstrate some discontinuity of the helix angle of the cardiomyocyte aggregates at the interface between the two populations.

CMRSegTools: An open-source software enabling reproducible research in segmentation of acute myocardial infarct in CMR images.

In the last decade, a large number of clinical trials have been deployed using Cardiac Magnetic Resonance (CMR) to evaluate cardioprotective strategies aiming at reducing the irreversible myocardial damage at the time of reperfusion. In these studies, segmentation and quantification of myocardial infarct lesion are often performed with a commercial software or an in-house closed-source code development thus creating a barrier for reproducible research. This paper introduces CMRSegTools: an open-source application software designed for the segmentation and quantification of myocardial infarct lesion enabling full access to state-of-the-art segmentation methods and parameters, easy integration of new algorithms and standardised results sharing. This post-processing tool has been implemented as a plug-in for the OsiriX/Horos DICOM viewer leveraging its database management functionalities and user interaction features to provide a bespoke tool for the analysis of cardiac MR images on large clinical cohorts. CMRSegTools includes, among others, user-assisted segmentation of the left-ventricle, semi- and automatic lesion segmentation methods, advanced statistical analysis and visualisation based on the American Heart Association 17-segment model. New segmentation methods can be integrated into the plug-in by developing components based on image processing and visualisation libraries such as ITK and VTK in C++ programming language. CMRSegTools allows the creation of training and testing data sets (labeled features such as lesion, microvascular obstruction and remote ROI) for supervised Machine Learning methods, and enables the comparative assessment of lesion segmentation methods via a single and integrated platform. The plug-in has been successfully used by several CMR imaging studies.

Spatiotemporal motion estimation for respiratory-correlated imaging of the lungs.

PURPOSE: Four-dimensional computed tomography (4D CT) can provide patient-specific motion information for radiotherapy planning and delivery. Motion estimation in 4D CT is challenging due to the reduced image quality and the presence of artifacts. We aim to improve the robustness of deformable registration applied to respiratory-correlated imaging of the lungs, by using a global problem formulation and pursuing a restrictive parametrization for the spatiotemporal deformation model. METHODS: A spatial transformation based on free-form deformations was extended to the temporal domain, by explicitly modeling the trajectory using a cyclic temporal model based on B-splines. A global registration criterion allowed to consider the entire image sequence simultaneously and enforce the temporal coherence of the deformation throughout the respiratory cycle. To ensure a parametrization capable of capturing the dynamics of respiratory motion, a prestudy was performed on the temporal dimension separately. The temporal parameters were tuned by fitting them to diaphragm motion data acquired for a large patient group. Suitable properties were retained and applied to spatiotemporal registration of 4D CT data. Registration results were validated using large sets of landmarks and compared to consecutive spatial registrations. To illustrate the benefit of the spatiotemporal approach, we also assessed the performance in the presence of motion-induced artifacts. RESULTS: Cubic B-splines gave better or similar fitting results as lower orders and were selected because of their inherently stronger regularization. The fitting and registration errors increased gradually with the temporal control point spacing, representing a trade-off between achievable accuracy and sensitivity to noise and artifacts. A piecewise smooth trajectory model, allowing for a discontinuous change of speed at end-inhale, was found most suitable to account for the sudden changes of motion at this breathing phase. The spatiotemporal modeling allowed a reduction of the number of parameters of 45%, while maintaining registration accuracy within 0.1 mm. The approach reduced the sensitivity to artifacts. CONCLUSIONS: Spatiotemporal registration can provide accurate motion estimation for 4D CT and improves the robustness to artifacts.

Regional myocardial function at preclinical disease stage of hypertrophic cardiomyopathy in female gene variant carriers.

We recently showed more severe diastolic dysfunction at the time of myectomy in female compared to male patients with obstructive hypertrophic cardiomyopathy. Early recognition of aberrant cardiac contracility using cardiovascular magnetic resonance (CMR) imaging may identify women at risk of cardiac dysfunction. To define myocardial function at an early disease stage, we studied regional cardiac function using CMR imaging with tissue tagging in asymptomatic female gene variant carriers. CMR imaging with tissue tagging was done in 13 MYBPC3, 11 MYH7 and 6 TNNT2 gene carriers and 16 age-matched controls. Regional peak circumferential strain was derived from tissue tagging images of the basal and midventricular segments of the septum and lateral wall. Left ventricular wall thickness and global function were comparable between MYBPC3, MYH7, TNNT2 carriers and controls. MYH7 gene variant carriers showed a different strain pattern as compared to the other groups, with higher septal peak circumferential strain at the basal segments compared to the lateral wall, whereas MYBPC3, TNNT2 carriers and controls showed higher strain at the lateral wall compared to the septum. Only subtle gene-specific changes in strain pattern occur in the myocardium preceding development of cardiac hypertrophy. Overall, our study shows that there are no major contractile deficits in asymptomatic females carrying a pathogenic gene variant, which would justify the use of CMR imaging for earlier diagnosis.

In vivo estimation of normal left ventricular stiffness and contractility based on routine cine MR acquisition.

Post-myocardial infarction remodeling process is known to alter the mechanical properties of the heart. Biomechanical parameters, such as tissue stiffness and contractility, would be useful for clinicians to better assess the severity of the diseased heart. However, these parameters are difficult to obtain in the current clinical practice. In this paper, we estimated subject-specific in vivo myocardial stiffness and contractility from 21 healthy volunteers, based on left ventricle models constructed from data acquired from routine cardiac MR acquisition only. The subject-specific biomechanical parameters were quantified using an inverse finite-element modelling approach. The personalized models were evaluated against relevant clinical metrics extracted from the MR data, such as circumferential strain, wall thickness and fractional thickening. We obtained the ranges of healthy biomechanical indices of 1.60 ± 0.22 kPa for left ventricular stiffness and 95.13 ± 14.56 kPa for left ventricular contractility. These reference normal values can be used for future model-based investigation on the stiffness and contractility of ischemic myocardium.

A Spatiotemporal exploration and 3D modeling of blood flow in healthy carotid artery bifurcation from two modalities: Ultrasound-Doppler and phase contrast MRI.

In the present study, we investigated the velocity profile over the carotid bifurcation in ten healthy volunteers by combining velocity measurements from two imaging modalities (PC-MRI and US-Doppler) and hemodynamic modeling in order to determine the optimal combination for the most realistic velocity estimation. The workflow includes data acquisition, velocity profile extraction at three sites (CCA, ECA and ICA), the arterial geometrical model reconstruction, a mesh generation and a rheological modeling. The results showed that US-Doppler measurements yielded higher velocity values as compared to PC-MRI (about 26% shift in CCA, 52% in ECA and 53% in ICA). This implies higher simulated velocities based on US-Doppler inlet as compared to simulated velocities based on PC-MRI inlet. Overall, PC-MRI inlet based simulations are closer to measurements than US-Doppler inlet based simulations. Moreover, the measured velocities showed that blood flow keeps a parabolic sectional profile distal from CCA, ECA and ICA, while being quite disturbed in the carotid sinus with a significant decrease in magnitude making this site very prone to atherosclerosis.

Towards production-level cardiac image analysis with grids.

Production exploitation of cardiac image analysis tools is hampered by the lack of proper IT infrastructure in health institutions, the non trivial integration of heterogeneous codes in coherent analysis procedures, and the need to achieve complete automation of these methods. HealthGrids are promising technologies to address these difficulties. This paper details how they can be complemented by high level problem solving environments such as workflow managers to improve the performance of applications both in terms of execution time and robustness of results. Two of the most important important cardiac image analysis tasks are considered, namely myocardium segmentation and motion estimation in a 4D sequence. Results are shown on the corresponding pipelines, using two different execution environments on the EGEE grid production infrastructure.

A gradient-based optical-flow cardiac motion estimation method for cine and tagged MR images.

A new method is proposed to quantify the myocardial motion from both 2D C(ine)-MRI and T(agged)-MRI sequences. The tag pattern offers natural landmarks within the image that makes it possible to accurately quantify the motion within the myocardial wall. Therefore, several methods have been proposed for T-MRI. However, the lack of salient features within the cardiac wall in C-MRI hampers local motion estimation. Our method aims to ensure the local intensity and shape features invariance during motion through the iterative minimization of a cost function via a random walk scheme. The proposed approach is evaluated on realistic simulated C-MRI and T-MRI sequences. The results show more than 53% improvements on displacement estimation, and more than 24% on strain estimation for both C-MRI and T-MRI sequences, as compared to state-of-the-art cardiac motion estimators. Preliminary experiments on clinical data have shown a good ability of the proposed method to detect abnormal motion patterns related to pathology. If those results are confirmed on large databases, this would open up the possibility for more accurate diagnosis of cardiac function from standard C-MRI examinations and also the retrospective study of prior studies.

Myocardial adaptation after surgical therapy differs for aortic valve stenosis and hypertrophic obstructive cardiomyopathy.

Surgical therapies in aortic valve stenosis (AVS) and hypertrophic obstructive cardiomyopathy (HOCM) aim to relief intraventricular pressure overload and improve clinical outcome. It is currently unknown to what extent myocardial adaptation concurs with restoration of intraventricular pressures, and whether this is similar in both patient groups. The aim of this study was to investigate changes in myocardial adaptation after surgical therapies for AVS and HOCM. Ten AVS and ten HOCM patients were enrolled and underwent cardiac magnetic resonance cine imaging and myocardial tagging prior to, and 4 months after aortic valve replacement (AVR) and septal myectomy, respectively. Global left ventricular (LV) analyses were derived from cine images. Circumferential strain was assessed from myocardial tagging images at the septal and lateral wall of the mid ventricle. Pressure gradients significantly decreased in both AVS and HOCM after surgery (p < 0.01), with a concomitant decrease in left atrial volume (p < 0.05) suggesting lower diastolic filling pressures. Also, LV volumes, mass and septal wall thickness decreased in both, but to a larger extent in AVS than in HOCM patients. AVR improved wall thickening (p < 0.05) and did not change systolic strain rate. Myectomy did not affect wall thickening and reduced septal systolic strain rate (p = 0.03). Both AVR and myectomy induced positive structural remodeling in line with a reduction of pressure overload. A concomitant recovery in systolic function however was found in AVR only. The systolic functional deterioration in HOCM patients seems to be inherent to myectomy and the ongoing and irreversible disease.

Quantitative comparison of human myocardial fiber orientations derived from DTI and polarized light imaging.

Diffusion tensor imaging (DTI) is a non-invasive technique used to obtain the 3D fiber structure of whole human hearts, for both in vivo and ex vivo cases. However, by essence, DTI does not measure directly the orientations of myocardial fibers. In contrast, polarized light imaging (PLI) allows for physical measurements of fiber orientations, but only for ex vivo case. This work aims at quantitatively comparing the myocardial fiber orientations of whole human hearts obtained from cardiac DTI with those measured by PLI. Whole human neonatal and infant hearts were first imaged using DTI. The same whole hearts were then imaged using PLI. After DTI and PLI data are registered, the orientations of fibers from the two imaging modalities were finally quantitatively compared. The results show that DTI and PLI have similar variation patterns of elevation and azimuth angles, with some differences in transmural elevation angle range. DTI itself induces an underestimation of the range of transmural elevation angles by a factor of about 25° at the basal and equatorial slices and the reduction of spatial resolution further decreases this range. PLI data exhibit a 15° ± 5° (P < 0.01) narrower transmural elevation angle range at apical slices than in basal or equatorial slices. This phenomenon is not observed in DTI data. In both modalities, the azimuth angle maps exhibit curved or twisting boundaries at the basal and apical slices. The experimental results globally enforce DTI as a valid imaging technique to reasonably characterize fiber orientations of the human heart noninvasively.

Strain imaging to predict response to cardiac resynchronization therapy: a systematic comparison of strain parameters using multiple imaging techniques.

AIMS: Various strain parameters and multiple imaging techniques are presently available including cardiovascular magnetic resonance (CMR) tagging (CMR-TAG), CMR feature tracking (CMR-FT), and speckle tracking echocardiography (STE). This study aims to compare predictive performance of different strain parameters and evaluate results per imaging technique to predict cardiac resynchronization therapy (CRT) response. METHODS AND RESULTS: Twenty-seven patients were prospectively enrolled and underwent CMR and echocardiographic examination before CRT implantation. Strain analysis was performed in circumferential (CMR-TAG, CMR-FT, and STE-circ) and longitudinal (STE-long) orientations. Regional strain values, parameters of dyssynchrony, and discoordination were calculated. After 12 months, CRT response was measured by the echocardiographic change in left ventricular (LV) end-systolic volume (LVESV). Twenty-six patients completed follow-up; mean LVESV change was -29 ± 27% with 17 (65%) patients showing ≥15% LVESV reduction. Measures of dyssynchrony (SD-TTPLV ) and discoordination (ISFLV ) were strongly related to CRT response when using CMR-TAG (R2 0.61 and R2 0.57, respectively), but showed poor correlations for CMR-FT and STE (all R2  ≤ 0.32). In contrast, the end-systolic septal strain (ESSsep ) parameter showed a consistent high correlation with LVESV change for all techniques (CMR-TAG R2 0.60; CMR-FT R2 0.50; STE-circ R2 0.43; and STE-long R2 0.43). After adjustment for QRS duration and QRS morphology, ESSsep remained an independent predictor of response per technique. CONCLUSIONS: End-systolic septal strain was the only parameter with a consistent good relation to reverse remodelling after CRT, irrespective of assessment technique. In clinical practice, this measure can be obtained by any available strain imaging technique and provides predictive value on top of current guideline criteria.

Comparison of strain imaging techniques in CRT candidates: CMR tagging, CMR feature tracking and speckle tracking echocardiography.

Parameters using myocardial strain analysis may predict response to cardiac resynchronization therapy (CRT). As the agreement between currently available strain imaging modalities is unknown, three different modalities were compared. Twenty-seven CRT-candidates, prospectively included in the MARC study, underwent cardiac magnetic resonance (CMR) imaging and echocardiographic examination. Left ventricular (LV) circumferential strain was analysed with CMR tagging (CMR-TAG), CMR feature tracking (CMR-FT), and speckle tracking echocardiography (STE). Basic strain values and parameters of dyssynchrony and discoordination obtained with CMR-FT and STE were compared to CMR-TAG. Agreement of CMR-FT and CMR-TAG was overall fair, while agreement between STE and CMR-TAG was often poor. For both comparisons, agreement on discoordination parameters was highest, followed by dyssynchrony and basic strain parameters. For discoordination parameters, agreement on systolic stretch index was highest, with fair intra-class correlation coefficients (ICC) (CMR-FT: 0.58, STE: 0.55). ICC of septal systolic rebound stretch (SRSsept) was poor (CMR-FT: 0.41, STE: 0.30). Internal stretch factor of septal and lateral wall (ISFsep-lat) showed fair ICC values (CMR-FT: 0.53, STE: 0.46), while the ICC of the total LV (ISFLV) was fair for CMR-FT (0.55) and poor for STE (ICC: 0.32). The CURE index had a fair ICC for both comparisons (CMR-FT: 0.49, STE 0.41). Although comparison of STE to CMR-TAG was limited by methodological differences, agreement between CMR-FT and CMR-TAG was overall higher compared to STE and CMR-TAG. CMR-FT is a potential clinical alternative for CMR-TAG and STE, especially in the detection of discoordination in CRT-candidates.