An automatic service for the personalization of ventricular cardiac meshes

Left Atrial Roof Enlargement Is a Distinct Feature of Heart Failure With Preserved Ejection Fraction

BACKGROUND

It remains unknown to what extent intrinsic atrial cardiomyopathy or left ventricular diastolic dysfunction drive atrial remodeling and functional failure in heart failure with preserved ejection fraction (HFpEF). Computational 3-dimensional (3D) models fitted to cardiovascular magnetic resonance allow state-of-the-art anatomic and functional assessment, and we hypothesized to identify a phenotype linked to HFpEF.

METHODS

Patients with exertional dyspnea and diastolic dysfunction on echocardiography (E/e', >8) were prospectively recruited and classified as HFpEF or noncardiac dyspnea based on right heart catheterization. All patients underwent rest and exercise-stress right heart catheterization and cardiovascular magnetic resonance. Computational 3D anatomic left atrial (LA) models were generated based on short-axis cine sequences. A fully automated pipeline was developed to segment cardiovascular magnetic resonance images and build 3D statistical models of LA shape and find the 3D patterns discriminant between HFpEF and noncardiac dyspnea. In addition, atrial morphology and function were quantified by conventional volumetric analyses and deformation imaging. A clinical follow-up was conducted after 24 months for the evaluation of cardiovascular hospitalization.

RESULTS

Beyond atrial size, the 3D LA models revealed roof dilation as the main feature found in masked HFpEF (diagnosed during exercise-stress only) preceding a pattern shift to overall atrial size in overt HFpEF (diagnosed at rest). Characteristics of the 3D model were integrated into the LA HFpEF shape score, a biomarker to characterize the gradual remodeling between noncardiac dyspnea and HFpEF. The LA HFpEF shape score was able to discriminate HFpEF (n=34) to noncardiac dyspnea (n=34; area under the curve, 0.81) and was associated with a risk for atrial fibrillation occurrence (hazard ratio, 1.02 [95% CI, 1.01-1.04]; P=0.003), as well as cardiovascular hospitalization (hazard ratio, 1.02 [95% CI, 1.00-1.04]; P=0.043).

CONCLUSIONS

LA roof dilation is an early remodeling pattern in masked HFpEF advancing to overall LA enlargement in overt HFpEF. These distinct features predict the occurrence of atrial fibrillation and cardiovascular hospitalization.

REGISTRATION

URL: https://www.clinicaltrials.gov; Unique identifier: NCT03260621.

Multi-class point cloud completion networks for 3D cardiac anatomy reconstruction from cine magnetic resonance images

Cine magnetic resonance imaging (MRI) is the current gold standard for the assessment of cardiac anatomy and function. However, it typically only acquires a set of two-dimensional (2D) slices of the underlying three-dimensional (3D) anatomy of the heart, thus limiting the understanding and analysis of both healthy and pathological cardiac morphology and physiology. In this paper, we propose a novel fully automatic surface reconstruction pipeline capable of reconstructing multi-class 3D cardiac anatomy meshes from raw cine MRI acquisitions. Its key component is a multi-class point cloud completion network (PCCN) capable of correcting both the sparsity and misalignment issues of the 3D reconstruction task in a unified model. We first evaluate the PCCN on a large synthetic dataset of biventricular anatomies and observe Chamfer distances between reconstructed and gold standard anatomies below or similar to the underlying image resolution for multiple levels of slice misalignment. Furthermore, we find a reduction in reconstruction error compared to a benchmark 3D U-Net by 32% and 24% in terms of Hausdorff distance and mean surface distance, respectively. We then apply the PCCN as part of our automated reconstruction pipeline to 1000 subjects from the UK Biobank study in a cross-domain transfer setting and demonstrate its ability to reconstruct accurate and topologically plausible biventricular heart meshes with clinical metrics comparable to the previous literature. Finally, we investigate the robustness of our proposed approach and observe its capacity to successfully handle multiple common outlier conditions.

Reconstruction and completion of high-resolution 3D cardiac shapes using anisotropic CMRI segmentations and continuous implicit neural representations

Since the onset of computer-aided diagnosis in medical imaging, voxel-based segmentation has emerged as the primary methodology for automatic analysis of left ventricle (LV) function and morphology in cardiac magnetic resonance images (CMRI). In standard clinical practice, simultaneous multi-slice 2D cine short-axis MR imaging is performed under multiple breath-holds resulting in highly anisotropic 3D images. Furthermore, sparse-view CMRI often lacks whole heart coverage caused by large slice thickness and often suffers from inter-slice misalignment induced by respiratory motion. Therefore, these volumes only provide limited information about the true 3D cardiac anatomy which may hamper highly accurate assessment of functional and anatomical abnormalities. To address this, we propose a method that learns a continuous implicit function representing 3D LV shapes by training an auto-decoder. For training, high-resolution segmentations from cardiac CT angiography are used. The ability of our approach to reconstruct and complete high-resolution shapes from manually or automatically obtained sparse-view cardiac shape information is evaluated by using paired high- and low-resolution CMRI LV segmentations. The results show that the reconstructed LV shapes have an unconstrained subvoxel resolution and appear smooth and plausible in through-plane direction. Furthermore, Bland-Altman analysis reveals that reconstructed high-resolution ventricle volumes are closer to the corresponding reference volumes than reference low-resolution volumes with bias of [limits of agreement] -3.51 [-18.87, 11.85] mL, and 12.96 [-10.01, 35.92] mL respectively. Finally, the results demonstrate that the proposed approach allows recovering missing shape information and can indirectly correct for limited motion-induced artifacts.

Left ventricular anatomy in obstructive hypertrophic cardiomyopathy: beyond basal septal hypertrophy

AIMS

Obstructive hypertrophic cardiomyopathy (oHCM) is characterized by dynamic obstruction of the left ventricular (LV) outflow tract (LVOT). Although this may be mediated by interplay between the hypertrophied septal wall, systolic anterior motion of the mitral valve, and papillary muscle abnormalities, the mechanistic role of LV shape is still not fully understood. This study sought to identify the LV end-diastolic morphology underpinning oHCM.

METHODS AND RESULTS

Cardiovascular magnetic resonance images from 2398 HCM individuals were obtained as part of the NHLBI HCM Registry. Three-dimensional LV models were constructed and used, together with a principal component analysis, to build a statistical shape model capturing shape variations. A set of linear discriminant axes were built to define and quantify (Z-scores) the characteristic LV morphology associated with LVOT obstruction (LVOTO) under different physiological conditions and the relationship between LV phenotype and genotype. The LV remodelling pattern in oHCM consisted not only of basal septal hypertrophy but a combination with LV lengthening, apical dilatation, and LVOT inward remodelling. Salient differences were observed between obstructive cases at rest and stress. Genotype negative cases showed a tendency towards more obstructive phenotypes both at rest and stress.

CONCLUSIONS

LV anatomy underpinning oHCM consists of basal septal hypertrophy, apical dilatation, LV lengthening, and LVOT inward remodelling. Differences between oHCM cases at rest and stress, as well as the relationship between LV phenotype and genotype, suggest different mechanisms for LVOTO. Proposed Z-scores render an opportunity of redefining management strategies based on the relationship between LV anatomy and LVOTO.

Estimation of anisotropic properties of CMR patient-specific left ventricle using the virtual field method

Left ventricle (LV) myocardial dysfunction has been recently investigated using the estimation of isotropic myocardial stiffness from magnetic resonance imaging (MRI). However, Myocardium is known to have a 3D complex geometry with anisotropic stiffness. The assessment of the anisotropy properties characterizes structural changes in myocardium as a consequence of heart failure (HF). From image data, the virtual field method (VFM) can determine material stiffness in a non-invasive manner. In the present work, the objective is to compare two inverse identification methods, given the isotropic and anisotropic models in the characterization of properties of myocardium in acute lymphoblastic leukemia (ALL) survivors using VFM and MRI. Two types of VFM approach are presented. Using the first, the virtual displacements (VFs) allow whole-field LV to be imposed into VFM formulation and caused to directly estimate two independent parameters from isotropic constitutive relation. With the second, anisotropic parameters are estimated using piece-wise (Finite element-based) VFM. The resulting values showed significant differences between the subjects in comparative study of leukemia survivors, and variance in estimated parameters by two different VFM approach. This approach would be an efficient tool to characterize early cardiac dysfunction. This work elucidates the benefits and shortcomings of using VFM to determine anisotropic parameters of LV myocardium in linear elastic and of using the FEM application to generate meshes of patient-specific LVs from MRI images.

Interpretable cardiac anatomy modeling using variational mesh autoencoders

Cardiac anatomy and function vary considerably across the human population with important implications for clinical diagnosis and treatment planning. Consequently, many computer-based approaches have been developed to capture this variability for a wide range of applications, including explainable cardiac disease detection and prediction, dimensionality reduction, cardiac shape analysis, and the generation of virtual heart populations. In this work, we propose a variational mesh autoencoder (mesh VAE) as a novel geometric deep learning approach to model such population-wide variations in cardiac shapes. It embeds multi-scale graph convolutions and mesh pooling layers in a hierarchical VAE framework to enable direct processing of surface mesh representations of the cardiac anatomy in an efficient manner. The proposed mesh VAE achieves low reconstruction errors on a dataset of 3D cardiac meshes from over 1,000 patients with acute myocardial infarction, with mean surface distances between input and reconstructed meshes below the underlying image resolution. We also find that it outperforms a voxelgrid-based deep learning benchmark in terms of both mean surface distance and Hausdorff distance while requiring considerably less memory. Furthermore, we explore the quality and interpretability of the mesh VAE's latent space and showcase its ability to improve the prediction of major adverse cardiac events over a clinical benchmark. Finally, we investigate the method's ability to generate realistic virtual populations of cardiac anatomies and find good alignment between the synthesized and gold standard mesh populations in terms of multiple clinical metrics.

Beyond Simpson's Rule: Accounting for Orientation and Ellipticity Assumptions

Simpson's biplane rule (SBR) is considered the gold standard method for left ventricle (LV) volume quantification from echocardiography but relies on a summation-of-disks approach that makes assumptions about LV orientation and cross-sectional shape. We aim to identify key limiting factors in SBR and to develop a new robust standard for volume quantification. Three methods for computing LV volume were studied: (i) SBR, (ii) addition of a truncated basal cone (TBC) to SBR and (iii) a novel method of basal-oriented disks (BODs). Three retrospective cohorts representative of the young, adult healthy and heart failure populations were used to study the impact of anatomical variations in volume computations. Results reveal how basal slanting can cause over- and underestimation of volume, with errors by SBR and TBC >10 mL for slanting angles >6°. Only the BOD method correctly accounted for basal slanting, reducing relative volume errors by SBR from -2.23 ± 2.21% to -0.70 ± 1.91% in the adult population and similar qualitative performance in the other two cohorts. In conclusion, the summation of basal oriented disks, a novel interpretation of SBR, is a more accurate and precise method for estimating LV volume.

A three-dimensional atlas of child's cardiac anatomy and the unique morphological alterations associated with obesity

AIMS

Statistical shape models (SSMs) of cardiac anatomy provide a new approach for analysis of cardiac anatomy. In adults, specific cardiac morphologies associate with cardiovascular risk factors and early disease stages. However, the relationships between morphology and risk factors in children remain unknown. We propose an SSM of the paediatric left ventricle to describe its morphological variability, examine its relationship with biometric parameters and identify adverse anatomical remodelling associated with obesity.

METHODS AND RESULTS

This cohort includes 2631 children (age 10.2 ± 0.6 years), mostly Western European (68.3%) with a balanced sex distribution (51.3% girls) from Generation R study. Cardiac magnetic resonance short-axis cine scans were segmented. Three-dimensional left ventricular (LV) meshes are automatically fitted to the segmentations to reconstruct the anatomies. We analyse the relationships between the LV anatomical features and participants' body surface area (BSA), age, and sex, and search for features uniquely related to obesity based on body mass index (BMI). In the SSM, 19 modes described over 90% of the population's LV shape variability. Main modes of variation were related to cardiac size, sphericity, and apical tilting. BSA, age, and sex were mostly correlated with modes describing LV size and sphericity. The modes correlated uniquely with BMI suggested that obese children present with septo-lateral tilting (R2 = 4.0%), compression in the antero-posterior direction (R2 = 3.3%), and decreased eccentricity (R2 = 2.0%).

CONCLUSIONS

We describe the variability of the paediatric heart morphology and identify anatomical features related to childhood obesity that could aid in risk stratification. Web service is released to provide access to the new shape parameters.

Understanding and Improving Risk Assessment After Myocardial Infarction Using Automated Left Ventricular Shape Analysis

BACKGROUND

Left ventricular ejection fraction (LVEF) and end-systolic volume (ESV) remain the main imaging biomarkers for post-acute myocardial infarction (AMI) risk stratification. However, they are limited to global systolic function and fail to capture functional and anatomical regional abnormalities, hindering their performance in risk stratification.

OBJECTIVES

This study aimed to identify novel 3-dimensional (3D) imaging end-systolic (ES) shape and contraction descriptors toward risk-related features and superior prognosis in AMI.

METHODS

A multicenter cohort of AMI survivors (n = 1,021; median age 63 years; 74.5% male) who underwent cardiac magnetic resonance (CMR) at a median of 3 days after infarction were considered for this study. The clinical endpoint was the 12-month rate of major adverse cardiac events (MACE; n = 73), consisting of all-cause death, reinfarction, and new congestive heart failure. A fully automated pipeline was developed to segment CMR images, build 3D statistical models of shape and contraction in AMI, and find the 3D patterns related to MACE occurrence.

RESULTS

The novel ES shape markers proved to be superior to ESV (median cross-validated area under the receiver-operating characteristic curve 0.681 [IQR: 0.679-0.684] vs 0.600 [IQR: 0.598-0.602]; P < 0.001); and 3D contraction to LVEF (0.716 [IQR: 0.714-0.718] vs 0.681 [IQR: 0.679-0.684]; P < 0.001) in MACE occurrence prediction. They also contributed to a significant improvement in a multivariable setting including CMR markers, cardiovascular risk factors, and basic patient characteristics (0.747 [IQR: 0.745-0.749]; P < 0.001). Based on these novel 3D descriptors, 3 impairments caused by AMI were identified: global, anterior, and basal, the latter being the most complementary signature to already known predictors.

CONCLUSIONS

The quantification of 3D differences in ES shape and contraction, enabled by a fully automated pipeline, improves post-AMI risk prediction and identifies shape and contraction patterns related to MACE occurrence.

Le Cœur en Sabot: shape associations with adverse events in repaired tetralogy of Fallot

BACKGROUND

Maladaptive remodelling mechanisms occur in patients with repaired tetralogy of Fallot (rToF) resulting in a cycle of metabolic and structural changes. Biventricular shape analysis may indicate mechanisms associated with adverse events independent of pulmonary regurgitant volume index (PRVI). We aimed to determine novel remodelling patterns associated with adverse events in patients with rToF using shape and function analysis.

METHODS

Biventricular shape and function were studied in 192 patients with rToF (median time from TOF repair to baseline evaluation 13.5 years). Linear discriminant analysis (LDA) and principal component analysis (PCA) were used to identify shape differences between patients with and without adverse events. Adverse events included death, arrhythmias, and cardiac arrest with median follow-up of 10 years.

RESULTS

LDA and PCA showed that shape characteristics pertaining to adverse events included a more circular left ventricle (LV) (decreased eccentricity), dilated (increased sphericity) LV base, increased right ventricular (RV) apical sphericity, and decreased RV basal sphericity. Multivariate LDA showed that the optimal discriminative model included only RV apical ejection fraction and one PCA mode associated with a more circular and dilated LV base (AUC = 0.77). PRVI did not add value, and shape changes associated with increased PRVI were not predictive of adverse outcomes.

CONCLUSION

Pathological remodelling patterns in patients with rToF are significantly associated with adverse events, independent of PRVI. Mechanisms related to incident events include LV basal dilation with a reduced RV apical ejection fraction.

Development and Evaluation of a Rib Statistical Shape Model for Thoracic Surgery

Patients with advanced cancer undergoing chest wall resection may require reconstruction. Currently, rib prostheses are created by segmenting computed tomography images, which is time-consuming and labour intensive. The aim was to optimise the production of digital rib models based on a patient's age, weight, height and gender. A statistical shape model of human ribs was created and used to synthetise rib models, which were compared to the ones produced by segmentation and mirroring. The segmentation took 11.56±1.60 min compared to 0.027 ±0.009 min using the new technique. The average mesh error between the mirroring technique and segmentation was 0.58±0.25 mm (right ribs), and 0.87±0.18 mm (left ribs), compared to 1.37±0.66 mm ( ) and 1.68 ±0.77 mm ( ), respectively, for the new technique. The new technique is promising for the efficiency and ease-of-use in the clinical environment. Clinical Relevance- This is an optimised 3D modelling method providing clinicians with a time-efficient technique to create patient-specific rib prostheses, without any expertise or software knowledge required.

A computationally efficient physiologically comprehensive 3D-0D closed-loop model of the heart and circulation

Computer models of cardiac electro-mechanics (EM) show promise as an effective means for the quantitative analysis of clinical data and, potentially, for predicting therapeutic responses. To realize such advanced applications methodological key challenges must be addressed. Enhanced computational efficiency and robustness is crucial to facilitate, within tractable time frames, model personalization, the simulation of prolonged observation periods under a broad range of conditions, and physiological completeness encompassing therapy-relevant mechanisms is needed to endow models with predictive capabilities beyond the mere replication of observations. Here, we introduce a universal feature-complete cardiac EM modeling framework that builds on a flexible method for coupling a 3D model of bi-ventricular EM to the physiologically comprehensive 0D CircAdapt model representing atrial mechanics and closed-loop circulation. A detailed mathematical description is given and efficiency, robustness, and accuracy of numerical scheme and solver implementation are evaluated. After parameterization and stabilization of the coupled 3D-0D model to a limit cycle under baseline conditions, the model's ability to replicate physiological behaviors is demonstrated, by simulating the transient response to alterations in loading conditions and contractility, as induced by experimental protocols used for assessing systolic and diastolic ventricular properties. Mechanistic completeness and computational efficiency of this novel model render advanced applications geared towards predicting acute outcomes of EM therapies feasible.

A computationally efficient physiologically comprehensive 3D-0D closed-loop model of the heart and circulation

Computer models of cardiac electro-mechanics (EM) show promise as an effective means for the quantitative analysis of clinical data and, potentially, for predicting therapeutic responses. To realize such advanced applications methodological key challenges must be addressed. Enhanced computational efficiency and robustness is crucial to facilitate, within tractable time frames, model personalization, the simulation of prolonged observation periods under a broad range of conditions, and physiological completeness encompassing therapy-relevant mechanisms is needed to endow models with predictive capabilities beyond the mere replication of observations. Here, we introduce a universal feature-complete cardiac EM modeling framework that builds on a flexible method for coupling a 3D model of bi-ventricular EM to the physiologically comprehensive 0D CircAdapt model representing atrial mechanics and closed-loop circulation. A detailed mathematical description is given and efficiency, robustness, and accuracy of numerical scheme and solver implementation are evaluated. After parameterization and stabilization of the coupled 3D-0D model to a limit cycle under baseline conditions, the model's ability to replicate physiological behaviors is demonstrated, by simulating the transient response to alterations in loading conditions and contractility, as induced by experimental protocols used for assessing systolic and diastolic ventricular properties. Mechanistic completeness and computational efficiency of this novel model render advanced applications geared towards predicting acute outcomes of EM therapies feasible.

Neural network-based left ventricle geometry prediction from CMR images with application in biomechanics

Combining biomechanical modelling of left ventricular (LV) function and dysfunction with cardiac magnetic resonance (CMR) imaging has the potential to improve the prognosis of patient-specific cardiovascular disease risks. Biomechanical studies of LV function in three dimensions usually rely on a computerized representation of the LV geometry based on finite element discretization, which is essential for numerically simulating in vivo cardiac dynamics. Detailed knowledge of the LV geometry is also relevant for various other clinical applications, such as assessing the LV cavity volume and wall thickness. Accurately and automatically reconstructing personalized LV geometries from conventional CMR images with minimal manual intervention is still a challenging task, which is a pre-requisite for any subsequent automated biomechanical analysis. We propose a deep learning-based automatic pipeline for predicting the three-dimensional LV geometry directly from routinely-available CMR cine images, without the need to manually annotate the ventricular wall. Our framework takes advantage of a low-dimensional representation of the high-dimensional LV geometry based on principal component analysis. We analyze how the inference of myocardial passive stiffness is affected by using our automatically generated LV geometries instead of manually generated ones. These insights will inform the development of statistical emulators of LV dynamics to avoid computationally expensive biomechanical simulations. Our proposed framework enables accurate LV geometry reconstruction, outperforming previous approaches by delivering a reconstruction error 50% lower than reported in the literature. We further demonstrate that for a nonlinear cardiac mechanics model, using our reconstructed LV geometries instead of manually extracted ones only moderately affects the inference of passive myocardial stiffness described by an anisotropic hyperelastic constitutive law. The developed methodological framework has the potential to make an important step towards personalized medicine by eliminating the need for time consuming and costly manual operations. In addition, our method automatically maps the CMR scan into a low-dimensional representation of the LV geometry, which constitutes an important stepping stone towards the development of an LV geometry-heterogeneous emulator.

High-Order Cardiomyopathy Human Heart Model and Mesh Generation

Faithful, accurate, and successful cardiac biomechanics and electrophysiological simulations require patient-specific geometric models of the heart. Since the cardiac geometry consists of highly-curved boundaries, the use of high-order meshes with curved elements would ensure that the various curves and features present in the cardiac geometry are well-captured and preserved in the corresponding mesh. Most other existing mesh generation techniques require computer-aided design files to represent the geometric boundary, which are often not available for biomedical applications. Unlike such methods, our technique takes a high-order surface mesh, generated from patient medical images, as input and generates a high-order volume mesh directly from the curved surface mesh. In this paper, we use our direct high-order curvilinear tetrahedral mesh generation method [1] to generate several second-order cardiac meshes. Our meshes include the left ventricle myocardia of a healthy heart and hearts with dilated and hypertrophic cardiomyopathy. We show that our high-order cardiac meshes do not contain inverted elements and are of sufficiently high quality for use in cardiac finite element simulations.

Learning-Based Regularization for Cardiac Strain Analysis via Domain Adaptation

Reliable motion estimation and strain analysis using 3D+ time echocardiography (4DE) for localization and characterization of myocardial injury is valuable for early detection and targeted interventions. However, motion estimation is difficult due to the low-SNR that stems from the inherent image properties of 4DE, and intelligent regularization is critical for producing reliable motion estimates. In this work, we incorporated the notion of domain adaptation into a supervised neural network regularization framework. We first propose a semi-supervised Multi-Layered Perceptron (MLP) network with biomechanical constraints for learning a latent representation that is shown to have more physiologically plausible displacements. We extended this framework to include a supervised loss term on synthetic data and showed the effects of biomechanical constraints on the network's ability for domain adaptation. We validated the semi-supervised regularization method on in vivo data with implanted sonomicrometers. Finally, we showed the ability of our semi-supervised learning regularization approach to identify infarct regions using estimated regional strain maps with good agreement to manually traced infarct regions from postmortem excised hearts.

A publicly available virtual cohort of four-chamber heart meshes for cardiac electro-mechanics simulations

Computational models of the heart are increasingly being used in the development of devices, patient diagnosis and therapy guidance. While software techniques have been developed for simulating single hearts, there remain significant challenges in simulating cohorts of virtual hearts from multiple patients. To facilitate the development of new simulation and model analysis techniques by groups without direct access to medical data, image analysis techniques and meshing tools, we have created the first publicly available virtual cohort of twenty-four four-chamber hearts. Our cohort was built from heart failure patients, age 67±14 years. We segmented four-chamber heart geometries from end-diastolic (ED) CT images and generated linear tetrahedral meshes with an average edge length of 1.1±0.2mm. Ventricular fibres were added in the ventricles with a rule-based method with an orientation of -60° and 80° at the epicardium and endocardium, respectively. We additionally refined the meshes to an average edge length of 0.39±0.10mm to show that all given meshes can be resampled to achieve an arbitrary desired resolution. We ran simulations for ventricular electrical activation and free mechanical contraction on all 1.1mm-resolution meshes to ensure that our meshes are suitable for electro-mechanical simulations. Simulations for electrical activation resulted in a total activation time of 149±16ms. Free mechanical contractions gave an average left ventricular (LV) and right ventricular (RV) ejection fraction (EF) of 35±1% and 30±2%, respectively, and a LV and RV stroke volume (SV) of 95±28mL and 65±11mL, respectively. By making the cohort publicly available, we hope to facilitate large cohort computational studies and to promote the development of cardiac computational electro-mechanics for clinical applications.

Preterm Birth Is a Novel, Independent Risk Factor for Altered Cardiac Remodeling and Early Heart Failure: Is it Time for a New Cardiomyopathy?

PURPOSE OF REVIEW

Around 10% of the global population is born preterm (< 37 weeks' gestation). Preterm birth is associated with an increased risk of cardiovascular events, with preterm-born individuals demonstrating a distinct cardiac phenotype. This review aims to summarize the main phenotypic features of the preterm heart and directions for future research to develop novel intervention strategies.

RECENT FINDINGS

Being born between 28 and 31 weeks' gestation results in a 4-fold higher risk of heart failure in childhood and adolescence and 17-fold increased risk when born less than 28 weeks' gestation. In support of this being due to a reduction in myocardial functional reserve, preterm-born young adults have an impaired left ventricular cardiac systolic response to moderate and high intensity physiological stress, despite having a preserved resting left ventricular ejection fraction. Similar impairments under physiological stress were also recently reported regarding the right ventricle in young adults born preterm. Preterm birth relates to a unique cardiac phenotype with an impaired response to stress conditions. These data, combined with the work in animal models, suggest that being born preterm may lead to a novel form of cardiomyopathy. Understanding the driving mechanisms leading to this unique cardiac phenotype is important to reduce risk of future heart failure and cardiovascular events.

Fast uncertainty quantification of activation sequences in patient-specific cardiac electrophysiology meeting clinical time constraints

We present a fast, patient-specific methodology for uncertainty quantification in electrophysiology, aimed at meeting the time constraints of clinical practitioners. We focus on computing the statistics of the activation map, given the uncertainties associated with the conductivity tensor modeling the fiber orientation in the heart. We use a fast parallel solution method implemented on a graphics processing unit for the eikonal approximation, in order to compute the activation map and to sample the random fiber field with correlation on the basis of geodesic distances. While this enables to perform uncertainty quantification studies with a manageable computational effort, the required time frame still exceeds clinically suitable time expectations. In order to reduce it further by 2 orders of magnitude, we rely on Bayesian multifidelity methods. In particular, we propose a low-fidelity model that is patient-specific and free from the additional training cost associated with reduced models. This is achieved by a sound physics-based simplification of the full eikonal model. The low-fidelity output is then corrected by the standard multifidelity framework. In practice, the complete procedure only requires approximately 100 new runs of our eikonal graphics processing unit solver for producing the sought estimates and their associated credible intervals, enabling a full online analysis in less than 5 minutes.

Patient-Specific Cardiovascular Computational Modeling: Diversity of Personalization and Challenges

Patient-specific computer models have been developed representing a variety of aspects of the cardiovascular system spanning the disciplines of electrophysiology, electromechanics, solid mechanics, and fluid dynamics. These physiological mechanistic models predict macroscopic phenomena such as electrical impulse propagation and contraction throughout the entire heart as well as flow and pressure dynamics occurring in the ventricular chambers, aorta, and coronary arteries during each heartbeat. Such models have been used to study a variety of clinical scenarios including aortic aneurysms, coronary stenosis, cardiac valvular disease, left ventricular assist devices, cardiac resynchronization therapy, ablation therapy, and risk stratification. After decades of research, these models are beginning to be incorporated into clinical practice directly via marketed devices and indirectly by improving our understanding of the underlying mechanisms of health and disease within a clinical context.

Influence of atrial contraction dynamics on cardiac function

In recent years, there has been a move from monoventricular or biventricular models of the heart, to more complex models that incorporate the electromechanical function in all 4 chambers. However, the biophysical foundation is still underdeveloped, with most work in atrial cellular models having focused on electrophysiological properties. Here, we present a biophysical model of human atrial contraction at body temperature and use it to study the effects of atrial contraction on whole organ function and a study of the effects of remodelling due to atrial fibrillation on atrial and ventricular function.

Decreasing Compensatory Ability of Concentric Ventricular Hypertrophy in Aortic-Banded Rat Hearts

The cardiac system compensates for variations in physiological and pathophysiological conditions through a dynamic remodeling at the organ, tissue, and intracellular levels in order to maintain function. However, on longer time scales following the onset of ventricular pressure overload, such remodeling may begin to inhibit physiological function and ultimately lead to heart failure. This progression from compensatory to decompensatory behavior is poorly understood, in particular owing to the absence of a unified perspective of the concomitantly remodeling subsystems. To address this issue, the present study investigates the evolution of compensatory mechanisms, in response to overload, by integrating diffusion-tensor MRI, echocardiography, and intracellular and hemodynamic measurements within consistent computational simulations of aortic-banded rat hearts. This approach allows a comparison of the relative leverage of different cardiac properties (geometry, passive mechanical stiffness, fiber configuration, diastolic and peak calcium concentrations, calcium-binding affinity, and aortic impedance) to affect cardiac contraction. Measurements indicate that, following aortic banding, an ejection fraction (EF) of 75% was maintained, relative to control rats, despite significant remodeling of the left-ventricular wall thickness (increasing by ~90% over 4 weeks). Applying our framework, we identified the left-ventricular wall thickness (concentric hypertrophy) and the intracellular calcium dynamics as playing the dominant roles in preserving EF acutely, whereas the significance of hypertrophy decreased subsequently. This trend suggests an increasing reliance on intracellular mechanisms (average increase ~50%), rather than on anatomical features (average decrease ~60%), to achieve compensation of pump function in the early phase of heart failure.

Advances in computational modelling for personalised medicine after myocardial infarction

Myocardial infarction (MI) is a leading cause of premature morbidity and mortality worldwide. Determining which patients will experience heart failure and sudden cardiac death after an acute MI is notoriously difficult for clinicians. The extent of heart damage after an acute MI is informed by cardiac imaging, typically using echocardiography or sometimes, cardiac magnetic resonance (CMR). These scans provide complex data sets that are only partially exploited by clinicians in daily practice, implying potential for improved risk assessment. Computational modelling of left ventricular (LV) function can bridge the gap towards personalised medicine using cardiac imaging in patients with post-MI. Several novel biomechanical parameters have theoretical prognostic value and may be useful to reflect the biomechanical effects of novel preventive therapy for adverse remodelling post-MI. These parameters include myocardial contractility (regional and global), stiffness and stress. Further, the parameters can be delineated spatially to correspond with infarct pathology and the remote zone. While these parameters hold promise, there are challenges for translating MI modelling into clinical practice, including model uncertainty, validation and verification, as well as time-efficient processing. More research is needed to (1) simplify imaging with CMR in patients with post-MI, while preserving diagnostic accuracy and patient tolerance (2) to assess and validate novel biomechanical parameters against established prognostic biomarkers, such as LV ejection fraction and infarct size. Accessible software packages with minimal user interaction are also needed. Translating benefits to patients will be achieved through a multidisciplinary approach including clinicians, mathematicians, statisticians and industry partners.

Disproportionate cardiac hypertrophy during early postnatal development in infants born preterm

BackgroundAdults born very preterm have increased cardiac mass and reduced function. We investigated whether a hypertrophic phenomenon occurs in later preterm infants and when this occurs during early development.MethodsCardiac ultrasound was performed on 392 infants (33% preterm at mean gestation 34±2 weeks). Scans were performed during fetal development in 137, at birth and 3 months of postnatal age in 200, and during both fetal and postnatal development in 55. Cardiac morphology and function was quantified and computational models created to identify geometric changes.ResultsAt birth, preterm offspring had reduced cardiac mass and volume relative to body size with a more globular heart. By 3 months, ventricular shape had normalized but both left and right ventricular mass relative to body size were significantly higher than expected for postmenstrual age (left 57.8±41.9 vs. 27.3±29.4%, P<0.001; right 39.3±38.1 vs. 16.6±40.8, P=0.002). Greater changes were associated with lower gestational age at birth (left P<0.001; right P=0.001).ConclusionPreterm offspring, including those born in late gestation, have a disproportionate increase in ventricular mass from birth up to 3 months of postnatal age. These differences were not present before birth. Early postnatal development may provide a window for interventions relevant to long-term cardiovascular health.

A model of cardiac contraction based on novel measurements of tension development in human cardiomyocytes

Experimental data from human cardiac myocytes at body temperature is crucial for a quantitative understanding of clinically relevant cardiac function and development of whole-organ computational models. However, such experimental data is currently very limited. Specifically, important measurements to characterize changes in tension development in human cardiomyocytes that occur with perturbations in cell length are not available. To address this deficiency, in this study we present an experimental data set collected from skinned human cardiac myocytes, including the passive and viscoelastic properties of isolated myocytes, the steady-state force calcium relationship at different sarcomere lengths, and changes in tension following a rapid increase or decrease in length, and after constant velocity shortening. This data set is, to our knowledge, the first characterization of length and velocity-dependence of tension generation in human skinned cardiac myocytes at body temperature. We use this data to develop a computational model of contraction and passive viscoelasticity in human myocytes. Our model includes troponin C kinetics, tropomyosin kinetics, a three-state crossbridge model that accounts for the distortion of crossbridges, and the cellular viscoelastic response. Each component is parametrized using our experimental data collected in human cardiomyocytes at body temperature. Furthermore we are able to confirm that properties of length-dependent activation at 37°C are similar to other species, with a shift in calcium sensitivity and increase in maximum tension. We revise our model of tension generation in the skinned isolated myocyte to replicate reported tension traces generated in intact muscle during isometric tension, to provide a model of human tension generation for multi-scale simulations. This process requires changes to calcium sensitivity, cooperativity, and crossbridge transition rates. We apply this model within multi-scale simulations of biventricular cardiac function and further refine the parametrization within the whole organ context, based on obtaining a healthy ejection fraction. This process reveals that crossbridge cycling rates differ between skinned myocytes and intact myocytes.

Novel Computational Analysis of Left Atrial Anatomy Improves Prediction of Atrial Fibrillation Recurrence after Ablation

The left atrium (LA) can change in size and shape due to atrial fibrillation (AF)-induced remodeling. These alterations can be linked to poorer outcomes of AF ablation. In this study, we propose a novel comprehensive computational analysis of LA anatomy to identify what features of LA shape can optimally predict post-ablation AF recurrence. To this end, we construct smooth 3D geometrical models from the segmentation of the LA blood pool captured in pre-procedural MR images. We first apply this methodology to characterize the LA anatomy of 144 AF patients and build a statistical shape model that includes the most salient variations in shape across this cohort. We then perform a discriminant analysis to optimally distinguish between recurrent and non-recurrent patients. From this analysis, we propose a new shape metric called vertical asymmetry, which measures the imbalance of size along the anterior to posterior direction between the superior and inferior left atrial hemispheres. Vertical asymmetry was found, in combination with LA sphericity, to be the best predictor of post-ablation recurrence at both 12 and 24 months (area under the ROC curve: 0.71 and 0.68, respectively) outperforming other shape markers and any of their combinations. We also found that model-derived shape metrics, such as the anterior-posterior radius, were better predictors than equivalent metrics taken directly from MRI or echocardiography, suggesting that the proposed approach leads to a reduction of the impact of data artifacts and noise. This novel methodology contributes to an improved characterization of LA organ remodeling and the reported findings have the potential to improve patient selection and risk stratification for catheter ablations in AF.

Improved identifiability of myocardial material parameters by an energy-based cost function

Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novel CF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the HF cases and [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the healthy case.

Improved identifiability of myocardial material parameters by an energy-based cost function

Myocardial stiffness is a valuable clinical biomarker for the monitoring and stratification of heart failure (HF). Cardiac finite element models provide a biomechanical framework for the assessment of stiffness through the determination of the myocardial constitutive model parameters. The reported parameter intercorrelations in popular constitutive relations, however, obstruct the unique estimation of material parameters and limit the reliable translation of this stiffness metric to clinical practice. Focusing on the role of the cost function (CF) in parameter identifiability, we investigate the performance of a set of geometric indices (based on displacements, strains, cavity volume, wall thickness and apicobasal dimension of the ventricle) and a novel CF derived from energy conservation. Our results, with a commonly used transversely isotropic material model (proposed by Guccione et al.), demonstrate that a single geometry-based CF is unable to uniquely constrain the parameter space. The energy-based CF, conversely, isolates one of the parameters and in conjunction with one of the geometric metrics provides a unique estimation of the parameter set. This gives rise to a new methodology for estimating myocardial material parameters based on the combination of deformation and energetics analysis. The accuracy of the pipeline is demonstrated in silico, and its robustness in vivo, in a total of 8 clinical data sets (7 HF and one control). The mean identified parameters of the Guccione material law were [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the HF cases and [Formula: see text] and [Formula: see text] ([Formula: see text], [Formula: see text], [Formula: see text]) for the healthy case.

Determinants of biventricular cardiac function: a mathematical model study on geometry and myofiber orientation

In patient-specific mathematical models of cardiac electromechanics, usually a patient-specific geometry and a generic myofiber orientation field are used as input, upon which myocardial tissue properties are tuned to clinical data. It remains unclear to what extent deviations in myofiber orientation and geometry between model and patient influence model predictions on cardiac function. Therefore, we evaluated the sensitivity of cardiac function for geometry and myofiber orientation in a biventricular (BiV) finite element model of cardiac mechanics. Starting out from a reference geometry in which myofiber orientation had no transmural component, two new geometries were defined with either a 27 % decrease in LV short- to long-axis ratio, or a 16 % decrease of RV length, but identical LV and RV cavity and wall volumes. These variations in geometry caused differences in both local myofiber and global pump work below 6 %. Variation of fiber orientation was induced through adaptive myofiber reorientation that caused an average change in fiber orientation of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\sim }8^\circ $$\end{document}∼8∘ predominantly through the formation of a component in transmural direction. Reorientation caused a considerable increase in local myofiber work \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$({\sim }18\,\%)$$\end{document}(∼18%) and in global pump work \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$({\sim }17\,\%)$$\end{document}(∼17%) in all three geometries, while differences between geometries were below 5 %. The findings suggest that implementing a realistic myofiber orientation is at least as important as defining a patient-specific geometry. The model for remodeling of myofiber orientation seems a useful approach to estimate myofiber orientation in the absence of accurate patient-specific information.

Right ventricular morphology and function following stage I palliation with a modified Blalock-Taussig shunt versus a right ventricle-to-pulmonary artery conduit

OBJECTIVES

The Norwood procedure for hypoplastic left heart syndrome (HLHS) is performed either via a right ventricle-to-pulmonary artery (RVPA) conduit or a modified Blalock–Taussig (MBT) shunt. Cardiac magnetic resonance (CMR) data was used to assess the effects of the RVPA conduit on ventricular shape and function through a computational analysis of anatomy and assessment of indices of strain.

METHODS

A retrospective analysis of 93 CMR scans of subjects with HLHS was performed (59 with MBT shunt, 34 with RVPA conduit), incorporating data at varying stages of surgery from two congenital centres. Longitudinal and short-axis cine images were used to create a computational cardiac atlas and assess global strain.

RESULTS

Those receiving an RVPA conduit had significant differences (P< 0.0001) in the shape of the RV corresponding to increased ventricular dilatation (P = 0.001) and increased sphericity (P = 0.006). Differences were evident only following completion of stage II surgery. Despite preserved ejection fraction in both groups, functional strain in the RVPA conduit group compared with that in the MBT shunt group was reduced across multiple ventricular axes, including a reduced systolic longitudinal strain rate (P< 0.0001), reduced diastolic longitudinal strain rate (P = 0.0001) and reduced midventricular systolic circumferential strain (P < 0.0001).

CONCLUSIONS

Computational modelling analysis reveals differences in ventricular remodelling in patients with HLHS undergoing an RVPA conduit insertion with focal scarring and volume loading leading to decreased functional markers of strain. The need for continued surveillance is warranted, as deleterious effects may not become apparent until later years.

Breast Milk Consumption in Preterm Neonates and Cardiac Shape in Adulthood

BACKGROUND AND OBJECTIVES

Preterm birth relates to long-term alterations in cardiac morphology and function. Understanding whether preterm postnatal life is a tractable period of cardiovascular development that can be positively altered by nutrition is relevant to long-term outcomes. We hypothesized that being fed human breast milk during early postnatal life is beneficial to long-term cardiac structure and function in preterm-born individuals compared with infant formulas.

METHODS

A total of 926 preterm-born infants originally took part in a randomized controlled trial of postnatal milk-feeding regimens between 1982 and 1985 across 5 different UK centers. Preterm-born individuals were randomly assigned to either breast milk donated by unrelated lactating women or nutrient-enriched formulas. We followed 102 individuals from this cohort: 30 of whom had been randomized to being fed exclusively human milk and 16 to being fed exclusively formula. As a comparison group, we recruited an additional 102 individuals born term to uncomplicated pregnancies. Cardiac morphology and function were assessed by MRI.

RESULTS

Preterm-born individuals fed exclusively human milk as infants had increased left and right ventricular end-diastolic volume index (+9.73%, P = .04 and +18.2%, P < .001) and stroke volume index (+9.79%, P = .05 and +22.1%, P = .01) compared with preterm-born individuals who were exclusively formula fed as infants.

CONCLUSIONS

This study provides the first evidence of a beneficial association between breast milk and cardiac morphology and function in adult life in those born preterm and supports promotion of human milk for the care of preterm infants to reduce long-term cardiovascular risk.

The relative role of patient physiology and device optimisation in cardiac resynchronisation therapy: A computational modelling study

Cardiac resynchronisation therapy (CRT) is an established treatment for heart failure, however the effective selection of patients and optimisation of therapy remain controversial. While extensive research is ongoing, it remains unclear whether improvements in patient selection or therapy planning offers a greater opportunity for the improvement of clinical outcomes. This computational study investigates the impact of both physiological conditions that guide patient selection and the optimisation of pacing lead placement on CRT outcomes. A multi-scale biophysical model of cardiac electromechanics was developed and personalised to patient data in three patients. These models were separated into components representing cardiac anatomy, pacing lead location, myocardial conductivity and stiffness, afterload, active contraction and conduction block for each individual, and recombined to generate a cohort of 648 virtual patients. The effect of these components on the change in total activation time of the ventricles (ΔTAT) and acute haemodynamic response (AHR) was analysed. The pacing site location was found to have the largest effect on ΔTAT and AHR. Secondary effects on ΔTAT and AHR were found for functional conduction block and cardiac anatomy. The simulation results highlight a need for a greater emphasis on therapy optimisation in order to achieve the best outcomes for patients.

•

Ventricular conduction block indicates patient response to CRT.

•

Placement of CRT pacing leads strongly affects response to therapy.

•

Improved treatment planning should be prioritised in order to maximise CRT outcomes.

Clinical Diagnostic Biomarkers from the Personalization of Computational Models of Cardiac Physiology

Computational modelling of the heart is rapidly advancing to the point of clinical utility. However, the difficulty of parameterizing and validating models from clinical data indicates that the routine application of truly predictive models remains a significant challenge. We argue there is significant value in an intermediate step towards prediction. This step is the use of biophysically based models to extract clinically useful information from existing patient data. Specifically in this paper we review methodologies for applying modelling frameworks for this goal in the areas of quantifying cardiac anatomy, estimating myocardial stiffness and optimizing measurements of coronary perfusion. Using these indicative examples of the general overarching approach, we finally discuss the value, ongoing challenges and future potential for applying biophysically based modelling in the clinical context.

Improving the stability of cardiac mechanical simulations

In the field of cardiac modeling, the mechanical action of the heart is often simulated using finite element methods. These simulations are becoming increasingly challenging as the computational domain is customized to a patient's anatomy, within which large heterogeneous tension gradients are generated via biophysical cell models which drive simulations of the cardiac pump cycle. The convergence of nonlinear solvers in simulations of large deformation mechanics depends on many factors. When extreme stress or irregular deformations are modeled, commonly used numerical methods can often fail to find a solution, which can prevent investigation of interesting parameter variations or use of models in a clinical context with high standards for robustness. This paper outlines a novel numerical method that is straightforward to implement and which significantly improves the stability of these simulations. The method involves adding a compressibility penalty to the standard incompressible formulation of large deformation mechanics. We compare the method's performance when used with both a direct discretization of the equations for incompressible solid mechanics, as well as the formulation based on an isochoric/deviatoric split of the deformation gradient. The addition of this penalty decreases the tendency for solutions to deviate from the incompressibility constraint, and significantly improves the ability of the Newton solver to find a solution. Additionally, our method maintains the expected order of convergence under mesh refinement, has nearly identical solutions for the pressure-volume relations, and stabilizes the solver to allow challenging simulations of both diastolic and systolic function on personalized patient geometries.

Big heart data: advancing health informatics through data sharing in cardiovascular imaging

The burden of heart disease is rapidly worsening due to the increasing prevalence of obesity and diabetes. Data sharing and open database resources for heart health informatics are important for advancing our understanding of cardiovascular function, disease progression and therapeutics. Data sharing enables valuable information, often obtained at considerable expense and effort, to be reused beyond the specific objectives of the original study. Many government funding agencies and journal publishers are requiring data reuse, and are providing mechanisms for data curation and archival. Tools and infrastructure are available to archive anonymous data from a wide range of studies, from descriptive epidemiological data to gigabytes of imaging data. Meta-analyses can be performed to combine raw data from disparate studies to obtain unique comparisons or to enhance statistical power. Open benchmark datasets are invaluable for validating data analysis algorithms and objectively comparing results. This review provides a rationale for increased data sharing and surveys recent progress in the cardiovascular domain. We also highlight the potential of recent large cardiovascular epidemiological studies enabling collaborative efforts to facilitate data sharing, algorithms benchmarking, disease modeling and statistical atlases.

Quantification of left ventricular torsion and diastolic recoil using cardiovascular magnetic resonance myocardial feature tracking

OBJECTIVES

Cardiovascular magnetic resonance feature tracking (CMR-FT) offers quantification of myocardial deformation from routine cine images. However, data using CMR-FT to quantify left ventricular (LV) torsion and diastolic recoil are not yet available. We therefore sought to evaluate the feasibility and reproducibility of CMR-FT to quantify LV torsion and peak recoil rate using an optimal anatomical approach.

METHODS

Short-axis cine stacks were acquired at rest and during dobutamine stimulation (10 and 20 µg · kg(-1) · min(-1)) in 10 healthy volunteers. Rotational displacement was analysed for all slices. A complete 3D-LV rotational model was developed using linear interpolation between adjacent slices. Torsion was defined as the difference between apical and basal rotation, divided by slice distance. Depending on the distance between the most apical (defined as 0% LV distance) and basal (defined as 100% LV distance) slices, four different models for the calculation of torsion were examined: Model-1 (25-75%), Model-2 (0-100%), Model-3 (25-100%) and Model-4 (0-75%). Analysis included subendocardial, subepicardial and global torsion and recoil rate (mean of subendocardial and subepicardial values).

RESULTS

Quantification of torsion and recoil rate was feasible in all subjects. There was no significant difference between the different models at rest. However, only Model-1 (25-75%) discriminated between rest and stress (Global Torsion: 2.7 ± 1.5° cm(-1), 3.6 ± 2.0° cm(-1), 5.1 ± 2.2° cm(-1), p<0.01; Global Recoil Rate: -30.1 ± 11.1° cm(-1) s(-1),-46.9 ± 15.0° cm(-1) s(-1),-68.9 ± 32.3° cm(-1) s(-1), p<0.01; for rest, 10 and 20 µg · kg(-)1 · min(-1) of dobutamine, respectively). Reproducibility was sufficient for all parameters as determined by Bland-Altman analysis, intraclass correlation coefficients and coefficient of variation.

CONCLUSIONS

CMR-FT based derivation of myocardial torsion and recoil rate is feasible and reproducible at rest and with dobutamine stress. Using an optimal anatomical approach measuring rotation at 25% and 75% apical and basal LV locations allows effective quantification of torsion and recoil dynamics. Application of these new measures of deformation by CMR-FT should next be explored in disease states.

Methodology for image-based reconstruction of ventricular geometry for patient-specific modeling of cardiac electrophysiology

Patient-specific modeling of ventricular electrophysiology requires an interpolated reconstruction of the 3-dimensional (3D) geometry of the patient ventricles from the low-resolution (Lo-res) clinical images. The goal of this study was to implement a processing pipeline for obtaining the interpolated reconstruction, and thoroughly evaluate the efficacy of this pipeline in comparison with alternative methods. The pipeline implemented here involves contouring the epi- and endocardial boundaries in Lo-res images, interpolating the contours using the variational implicit functions method, and merging the interpolation results to obtain the ventricular reconstruction. Five alternative interpolation methods, namely linear, cubic spline, spherical harmonics, cylindrical harmonics, and shape-based interpolation were implemented for comparison. In the thorough evaluation of the processing pipeline, Hi-res magnetic resonance (MR), computed tomography (CT), and diffusion tensor (DT) MR images from numerous hearts were used. Reconstructions obtained from the Hi-res images were compared with the reconstructions computed by each of the interpolation methods from a sparse sample of the Hi-res contours, which mimicked Lo-res clinical images. Qualitative and quantitative comparison of these ventricular geometry reconstructions showed that the variational implicit functions approach performed better than others. Additionally, the outcomes of electrophysiological simulations (sinus rhythm activation maps and pseudo-ECGs) conducted using models based on the various reconstructions were compared. These electrophysiological simulations demonstrated that our implementation of the variational implicit functions-based method had the best accuracy.

Images as drivers of progress in cardiac computational modelling

Computational models have become a fundamental tool in cardiac research. Models are evolving to cover multiple scales and physical mechanisms. They are moving towards mechanistic descriptions of personalised structure and function, including effects of natural variability. These developments are underpinned to a large extent by advances in imaging technologies. This article reviews how novel imaging technologies, or the innovative use and extension of established ones, integrate with computational models and drive novel insights into cardiac biophysics. In terms of structural characterization, we discuss how imaging is allowing a wide range of scales to be considered, from cellular levels to whole organs. We analyse how the evolution from structural to functional imaging is opening new avenues for computational models, and in this respect we review methods for measurement of electrical activity, mechanics and flow. Finally, we consider ways in which combined imaging and modelling research is likely to continue advancing cardiac research, and identify some of the main challenges that remain to be solved.