Robust image-based estimation of cardiac tissue parameters and their uncertainty from noisy data

Fast Posterior Estimation of Cardiac Electrophysiological Model Parameters via Bayesian Active Learning

Probabilistic estimation of cardiac electrophysiological model parameters serves an important step toward model personalization and uncertain quantification. The expensive computation associated with these model simulations, however, makes direct Markov Chain Monte Carlo (MCMC) sampling of the posterior probability density function (pdf) of model parameters computationally intensive. Approximated posterior pdfs resulting from replacing the simulation model with a computationally efficient surrogate, on the other hand, have seen limited accuracy. In this study, we present a Bayesian active learning method to directly approximate the posterior pdf function of cardiac model parameters, in which we intelligently select training points to query the simulation model in order to learn the posterior pdf using a small number of samples. We integrate a generative model into Bayesian active learning to allow approximating posterior pdf of high-dimensional model parameters at the resolution of the cardiac mesh. We further introduce new acquisition functions to focus the selection of training points on better approximating the shape rather than the modes of the posterior pdf of interest. We evaluated the presented method in estimating tissue excitability in a 3D cardiac electrophysiological model in a range of synthetic and real-data experiments. We demonstrated its improved accuracy in approximating the posterior pdf compared to Bayesian active learning using regular acquisition functions, and substantially reduced computational cost in comparison to existing standard or accelerated MCMC sampling.

Prediction of Left Ventricular Mechanics Using Machine Learning

The goal of this paper was to provide a real-time left ventricular (LV) mechanics simulator using machine learning (ML). Finite element (FE) simulations were conducted for the LV with different material properties to obtain a training set. A hyperelastic fiber-reinforced material model was used to describe the passive behavior of the myocardium during diastole. The active behavior of the heart resulting from myofiber contractions was added to the passive tissue during systole. The active and passive properties govern the LV constitutive equation. These mechanical properties were altered using optimal Latin hypercube design of experiments to obtain training FE models with varied active properties (volume and pressure predictions) and varied passive properties (stress predictions). For prediction of LV pressures, we used eXtreme Gradient Boosting (XGboost) and Cubist, and XGBoost was used for predictions of LV pressures, volumes as well as LV stresses. The LV pressure and volume results obtained from ML were similar to FE computations. The ML results could capture the shape of LV pressure as well as LV pressure-volume loops. The results predicted by Cubist were smoother than those from XGBoost. The mean absolute errors were as follows: XGBoost volume: 1.734 ± 0.584 ml, XGBoost pressure: 1.544 ± 0.298 mmHg, Cubist volume: 1.495 ± 0.260 ml, Cubist pressure: 1.623 ± 0.191 mmHg, myofiber stress: 0.334 ± 0.228 kPa, cross myofiber stress: 0.075 ± 0.024 kPa, and shear stress: 0.050 ± 0.032 kPa. The simulation results show ML can predict LV mechanics much faster than the FE method. The ML model can be used as a tool to predict LV behavior. Training of our ML model based on a large group of subjects can improve its predictability for real world applications.

Population-based priors in cardiac model personalisation for consistent parameter estimation in heterogeneous databases

Personalised cardiac models are a virtual representation of the patient heart, with parameter values for which the simulation fits the available clinical measurements. Models usually have a large number of parameters while the available data for a given patient are typically limited to a small set of measurements; thus, the parameters cannot be estimated uniquely. This is a practical obstacle for clinical applications, where accurate parameter values can be important. Here, we explore an original approach based on an algorithm called Iteratively Updated Priors (IUP), in which we perform successive personalisations of a full database through maximum a posteriori (MAP) estimation, where the prior probability at an iteration is set from the distribution of personalised parameters in the database at the previous iteration. At the convergence of the algorithm, estimated parameters of the population lie on a linear subspace of reduced (and possibly sufficient) dimension in which for each case of the database, there is a (possibly unique) parameter value for which the simulation fits the measurements. We first show how this property can help the modeller select a relevant parameter subspace for personalisation. In addition, since the resulting priors in this subspace represent the population statistics in this subspace, they can be used to perform consistent parameter estimation for cases where measurements are possibly different or missing in the database, which we illustrate with the personalisation of a heterogeneous database of 811 cases.

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.

High-resolution data assimilation of cardiac mechanics applied to a dyssynchronous ventricle

Computational models of cardiac mechanics, personalized to a patient, offer access to mechanical information above and beyond direct medical imaging. Additionally, such models can be used to optimize and plan therapies in-silico, thereby reducing risks and improving patient outcome. Model personalization has traditionally been achieved by data assimilation, which is the tuning or optimization of model parameters to match patient observations. Current data assimilation procedures for cardiac mechanics are limited in their ability to efficiently handle high-dimensional parameters. This restricts parameter spatial resolution, and thereby the ability of a personalized model to account for heterogeneities that are often present in a diseased or injured heart. In this paper, we address this limitation by proposing an adjoint gradient-based data assimilation method that can efficiently handle high-dimensional parameters. We test this procedure on a synthetic data set and provide a clinical example with a dyssynchronous left ventricle with highly irregular motion. Our results show that the method efficiently handles a high-dimensional optimization parameter and produces an excellent agreement for personalized models to both synthetic and clinical data.

Infarct Localization From Myocardial Deformation: Prediction and Uncertainty Quantification by Regression From a Low-Dimensional Space

Diagnosing and localizing myocardial infarct is crucial for early patient management and therapy planning. We propose a new method for predicting the location of myocardial infarct from local wall deformation, which has value for risk stratification from routine examinations such as (3D) echocardiography. The pipeline combines non-linear dimensionality reduction of deformation patterns and two multi-scale kernel regressions. Confidence in the diagnosis is assessed by a map of local uncertainties, which integrates plausible infarct locations generated from the space of reduced dimensionality. These concepts were tested on 500 synthetic cases generated from a realistic cardiac electromechanical model, and 108 pairs of 3D echocardiographic sequences and delayed-enhancement magnetic resonance images from real cases. Infarct prediction is made at a spatial resolution around 4 mm, more than 10 times smaller than the current diagnosis, made regionally. Our method is accurate, and significantly outperforms the clinically-used thresholding of the deformation patterns (on real data: sensitivity/specificity of 0.828/0.804, area under the curve: 0.909 versus 0.742 for the most predictive strain component). Uncertainty adds value to refine the diagnosis and eventually re-examine suspicious cases.

Multiphysics and multiscale modelling, data-model fusion and integration of organ physiology in the clinic: ventricular cardiac mechanics

With heart and cardiovascular diseases continually challenging healthcare systems worldwide, translating basic research on cardiac (patho)physiology into clinical care is essential. Exacerbating this already extensive challenge is the complexity of the heart, relying on its hierarchical structure and function to maintain cardiovascular flow. Computational modelling has been proposed and actively pursued as a tool for accelerating research and translation. Allowing exploration of the relationships between physics, multiscale mechanisms and function, computational modelling provides a platform for improving our understanding of the heart. Further integration of experimental and clinical data through data assimilation and parameter estimation techniques is bringing computational models closer to use in routine clinical practice. This article reviews developments in computational cardiac modelling and how their integration with medical imaging data is providing new pathways for translational cardiac modelling.

Towards Personalized Cardiology: Multi-Scale Modeling of the Failing Heart

BACKGROUND

Despite modern pharmacotherapy and advanced implantable cardiac devices, overall prognosis and quality of life of HF patients remain poor. This is in part due to insufficient patient stratification and lack of individualized therapy planning, resulting in less effective treatments and a significant number of non-responders.

METHODS AND RESULTS

State-of-the-art clinical phenotyping was acquired, including magnetic resonance imaging (MRI) and biomarker assessment. An individualized, multi-scale model of heart function covering cardiac anatomy, electrophysiology, biomechanics and hemodynamics was estimated using a robust framework. The model was computed on n=46 HF patients, showing for the first time that advanced multi-scale models can be fitted consistently on large cohorts. Novel multi-scale parameters derived from the model of all cases were analyzed and compared against clinical parameters, cardiac imaging, lab tests and survival scores to evaluate the explicative power of the model and its potential for better patient stratification. Model validation was pursued by comparing clinical parameters that were not used in the fitting process against model parameters.

CONCLUSION

This paper illustrates how advanced multi-scale models can complement cardiovascular imaging and how they could be applied in patient care. Based on obtained results, it becomes conceivable that, after thorough validation, such heart failure models could be applied for patient management and therapy planning in the future, as we illustrate in one patient of our cohort who received CRT-D implantation.