Three-dimensional regional bi-ventricular shape remodeling is associated with exercise capacity in endurance athletes

AIMS

Endurance athletes develop cardiac remodeling to cope with increased cardiac output during exercise. This remodeling is both anatomical and functional and shows large interindividual variability. In this study, we quantify local geometric ventricular remodeling related to long-standing endurance training and assess its relationship with cardiovascular performance during exercise.

METHODS

We extracted 3D models of the biventricular shape from end-diastolic cine magnetic resonance images acquired from a cohort of 89 triathlon athletes and 77 healthy sedentary subjects. Additionally, the athletes underwent cardio-pulmonary exercise testing, together with an echocardiographic study at baseline and few minutes after maximal exercise. We used statistical shape analysis to identify regional bi-ventricular shape differences between athletes and non-athletes.

RESULTS

The ventricular shape was significantly different between athletes and controls (p < 1e-6). The observed regional remodeling in the right heart was mainly a shift of the right ventricle (RV) volume distribution towards the right ventricular infundibulum, increasing the overall right ventricular volume. In the left heart, there was an increment of left ventricular mass and a dilation of the left ventricle. Within athletes, the amount of such remodeling was independently associated to higher peak oxygen pulse (p < 0.001) and weakly with greater post-exercise RV free wall longitudinal strain (p = 0.03).

CONCLUSIONS

We were able to identify specific bi-ventricular regional remodeling induced by long-lasting endurance training. The amount of remodeling was associated with better cardiopulmonary performance during an exercise test.

Benchmarking framework for myocardial tracking and deformation algorithms: an open access database

In this paper we present a benchmarking framework for the validation of cardiac motion analysis algorithms. The reported methods are the response to an open challenge that was issued to the medical imaging community through a MICCAI workshop. The database included magnetic resonance (MR) and 3D ultrasound (3DUS) datasets from a dynamic phantom and 15 healthy volunteers. Participants processed 3D tagged MR datasets (3DTAG), cine steady state free precession MR datasets (SSFP) and 3DUS datasets, amounting to 1158 image volumes. Ground-truth for motion tracking was based on 12 landmarks (4 walls at 3 ventricular levels). They were manually tracked by two observers in the 3DTAG data over the whole cardiac cycle, using an in-house application with 4D visualization capabilities. The median of the inter-observer variability was computed for the phantom dataset (0.77 mm) and for the volunteer datasets (0.84 mm). The ground-truth was registered to 3DUS coordinates using a point based similarity transform. Four institutions responded to the challenge by providing motion estimates for the data: Fraunhofer MEVIS (MEVIS), Bremen, Germany; Imperial College London - University College London (IUCL), UK; Universitat Pompeu Fabra (UPF), Barcelona, Spain; Inria-Asclepios project (INRIA), France. Details on the implementation and evaluation of the four methodologies are presented in this manuscript. The manually tracked landmarks were used to evaluate tracking accuracy of all methodologies. For 3DTAG, median values were computed over all time frames for the phantom dataset (MEVIS=1.20mm, IUCL=0.73 mm, UPF=1.10mm, INRIA=1.09 mm) and for the volunteer datasets (MEVIS=1.33 mm, IUCL=1.52 mm, UPF=1.09 mm, INRIA=1.32 mm). For 3DUS, median values were computed at end diastole and end systole for the phantom dataset (MEVIS=4.40 mm, UPF=3.48 mm, INRIA=4.78 mm) and for the volunteer datasets (MEVIS=3.51 mm, UPF=3.71 mm, INRIA=4.07 mm). For SSFP, median values were computed at end diastole and end systole for the phantom dataset(UPF=6.18 mm, INRIA=3.93 mm) and for the volunteer datasets (UPF=3.09 mm, INRIA=4.78 mm). Finally, strain curves were generated and qualitatively compared. Good agreement was found between the different modalities and methodologies, except for radial strain that showed a high variability in cases of lower image quality.

Simulation of late gadolinium enhancement cardiac magnetic resonance studies

In this study we propose a pipeline for simulation of late gadolinium enhancement images. We used a modified version of the XCAT phantom to improve simulation realism. Modifications included the modeling of trabeculae and papillary muscles, and the increase of sublabels to resemble tissue intensity variability. Magnetic properties for each body tissue were sampled in three settings: from Gaussian distributions, combining Rayleigh-Gaussian distributions, and from Rayleigh distributions. Thirty-two simulated datasets were compared with 32 clinical datasets from infarcted patients. Histograms were obtained for five tissues: lung, pericardium, myocardium, blood and hyper-enhanced area. Real and simulated histograms were compared with the Chi-square dissimilarity metric (χ(2)) and Kullback-Leibler divergence (KL). The generated simulated images look similar to real images according to both metrics. Rayleigh and the Rayleigh-Gaussian models obtained comparable average results (respectively: χ(2)= 0.16 ± 0.12 and 0.18 ± 0.11; KL=0.15 ± 0.17 and 0.16 ± 0.18).

3D mesh based wall thickness measurement: identification of left ventricular hypertrophy phenotypes

Left ventricular hypertrophy (LVH) is a complex cardiac condition mainly identified by the thickening of the myocardial wall. Although most of the contemporary cardiac imaging modalities provide high resolution 3D images, the wall thickness (WT) is still measured within the acquired planes. This way of measurement may introduce an error as cardiac wall is not necessarily orthogonal to the plane. In this study we analyze how different approaches to measure WT can affect an automatic identification of hypertrophy. The compared approaches are: WT measured along surface normal and the one provided by a medial surface. For both approaches we evaluated their ability to identify LVH phenotypes by testing with two classifiers: Transductive Confidence Machine-k Nearest Neighbor (TCM-kNN) and Linear Discriminant Analysis (LDA). Fifty three subjects were included in this study: 18 patients with hypertrophic cardiomyopathy (HCM), 13 patients with hypertensive heart disease (HDD) and 22 sedentary subjects (CG). Medial surface based approach allowed obtaining higher classification accuracy in HDD patients, while normal based approach allowed for higher classification accuracy in HCM patients.