Estimating prognosis in patients with acute myocardial infarction using personalized computational heart models

A multiscale finite element model of left ventricular mechanics incorporating baroreflex regulation

Cardiovascular function is regulated by a short-term hemodynamic baroreflex loop, which tries to maintain arterial pressure at a normal level. In this study, we present a new multiscale model of the cardiovascular system named MyoFE. This framework integrates a mechanistic model of contraction at the myosin level into a finite-element-based model of the left ventricle pumping blood through the systemic circulation. The model is coupled with a closed-loop feedback control of arterial pressure inspired by a baroreflex algorithm previously published by our team. The reflex loop mimics the afferent neuron pathway via a normalized signal derived from arterial pressure. The efferent pathway is represented by a kinetic model that simulates the net result of neural processing in the medulla and cell-level responses to autonomic drive. The baroreflex control algorithm modulates parameters such as heart rate and vascular tone of vessels in the lumped-parameter model of systemic circulation. In addition, it spatially modulates intracellular Ca2+ dynamics and molecular-level function of both the thick and the thin myofilaments in the left ventricle. Our study demonstrates that the baroreflex algorithm can maintain arterial pressure in the presence of perturbations such as acute cases of altered aortic resistance, mitral regurgitation, and myocardial infarction. The capabilities of this new multiscale model will be utilized in future research related to computational investigations of growth and remodeling.

Image-based estimation of the left ventricular cavity volume using deep learning and Gaussian process with cardio-mechanical applications

In this investigation, an image-based method has been developed to estimate the volume of the left ventricular cavity using cardiac magnetic resonance (CMR) imaging data. Deep learning and Gaussian processes have been applied to bring the estimations closer to the cavity volumes manually extracted. CMR data from 339 patients and healthy volunteers have been used to train a stepwise regression model that can estimate the volume of the left ventricular cavity at the beginning and end of diastole. We have decreased the root mean square error (RMSE) of cavity volume estimation approximately from 13 to 8 ml compared to the common practice in the literature. Considering the RMSE of manual measurements is approximately 4 ml on the same dataset, 8 ml of error is notable for a fully automated estimation method, which needs no supervision or user-hours once it has been trained. Additionally, to demonstrate a clinically relevant application of automatically estimated volumes, we inferred the passive material properties of the myocardium given the volume estimates using a well-validated cardiac model. These material properties can be further used for patient treatment planning and diagnosis.

Identification and quantification of chimeric sequencing reads in a highly multiplexed RAD-seq protocol

Highly multiplexed approaches have become common in genomic studies. They have improved the cost-effectiveness of genotyping hundreds of individuals using combinatorially barcoded adapters. These strategies, however, can potentially misassigned reads to incorrect samples. Here, we used a modified quaddRAD protocol to analyse the occurrence of index hopping and PCR chimeras in a series of experiments with up to 100 multiplexed samples per sequencing lane (639 samples in total). We created two types of sequencing libraries: four libraries of type A, where PCRs were run on individual samples before multiplexing, and three libraries of type B, where PCRs were run on pooled samples. We used fixed pairs of inner barcodes to identify chimeric reads. Type B libraries show a higher percentage of misassigned reads (1.15%) than type A libraries (0.65%). We also quantify the commonly undetectable chimeric sequences that occur whenever multiplexed groups of samples with different outer barcodes are sequenced together on a single flow cell. Our results suggest that these types of chimeric sequences represent up to 1.56% and 1.29% of reads in type A and B libraries, respectively. We also show that increasing the number of mismatches allowed for barcode rescue to above 2 dramatically increases the number of recovered chimeric reads. We provide recommendations for developing highly multiplexed RAD-seq protocols and analysing the resulting data to minimize the generation of chimeric sequences, allowing their quantification and a finer control on the number of PCR cycles necessary to generate enough input DNA for library preparation.

Dissecting the relationship between plasma and tissue metabolome in a cohort of women with obesity: Analysis of subcutaneous and visceral adipose, muscle, and liver

Untargeted metabolomics of blood samples has become widely applied to study metabolic alterations underpinning disease and to identify biomarkers. However, understanding the relevance of a blood metabolite marker can be challenging if it is unknown whether it reflects the concentration in relevant tissues. To explore this field, metabolomic and lipidomic profiles of plasma, four sites of adipose tissues (ATs) from peripheral or central depot, two sites of muscle tissue, and liver tissue from a group of nondiabetic women with obesity who were scheduled to undergo bariatric surgery (n = 21) or other upper GI surgery (n = 5), were measured by liquid chromatography coupled with mass spectrometry. Relationships between plasma and tissue profiles were examined using Pearson correlation analysis subject to Benjamini-Hochberg correction. Plasma metabolites and lipids showed the highest number of significantly positive correlations with their corresponding concentrations in liver tissue, including lipid species of ceramide, mono- and di-hexosylceramide, sphingomyelin, phosphatidylcholine (PC), phosphatidylethanolamine (PE), lysophosphatidylethanolamine, dimethyl phosphatidylethanolamine, ether-linked PC, ether-linked PE, free fatty acid, cholesteryl ester, diacylglycerol and triacylglycerol, and polar metabolites linked to several metabolic functions and gut microbial metabolism. Plasma also showed significantly positive correlations with muscle for several phospholipid species and polar metabolites linked to metabolic functions and gut microbial metabolism, and with AT for several triacylglycerol species. In conclusion, plasma metabolomic and lipidomic profiles were reflective more of the liver profile than any of the muscle or AT sites examined in the present study. Our findings highlighted the importance of taking into consideration the metabolomic relationship of various tissues with plasma when postulating plasma metabolites marker to underlying mechanisms occurring in a specific tissue.

Bayesian optimisation for efficient parameter inference in a cardiac mechanics model of the left ventricle

We consider parameter inference in cardio-mechanic models of the left ventricle, in particular the one based on the Holtzapfel-Ogden (HO) constitutive law, using clinical in vivo data. The equations underlying these models do not admit closed form solutions and hence need to be solved numerically. These numerical procedures are computationally expensive making computational run times associated with numerical optimisation or sampling excessive for the uptake of the models in the clinical practice. To address this issue, we adopt the framework of Bayesian optimisation (BO), which is an efficient statistical technique of global optimisation. BO seeks the optimum of an unknown black-box function by sequentially training a statistical surrogate-model and using it to select the next query point by leveraging the associated exploration-exploitation trade-off. To guarantee that the estimates based on the in vivo data are realistic also for high-pressures, unobservable in vivo, we include a penalty term based on a previously published empirical law developed using ex vivo data. Two case studies based on real data demonstrate that the proposed BO procedure outperforms the state-of-the-art inference algorithm for the HO constitutive law.

Effects of dispersed fibres in myocardial mechanics, Part I: passive response

It is widely acknowledged that an imbalanced biomechanical environment can have significant effects on myocardial pathology, leading to adverse remodelling of cardiac function if it persists. Accurate stress prediction essentially depends on the strain energy function which should have competent descriptive and predictive capabilities. Previous studies have focused on myofibre dispersion, but not on fibres along other directions. In this study, we will investigate how fibre dispersion affects myocardial biomechanical behaviours by taking into account both the myofibre dispersion and the sheet fibre dispersion, with a focus on the sheet fibre dispersion. Fibre dispersion is incorporated into a widely-used myocardial strain energy function using the discrete fibre bundle approach. We first study how different dispersion affects the descriptive capability of the strain energy function when fitting to ex vivo experimental data, and then the predictive capability in a human left ventricle during diastole. Our results show that the chosen strain energy function can achieve the best goodness-of-fit to the experimental data by including both fibre dispersion. Furthermore, noticeable differences in stress can be found in the LV model. Our results may suggest that it is necessary to include both dispersion for myofibres and the sheet fibres for the improved descriptive capability to the ex vivo experimental data and potentially more accurate stress prediction in cardiac mechanics.

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.

The role of regional myocardial topography post-myocardial infarction on infarct extension

Infarct extension involves necrosis of healthy myocardium in the border zone (BZ), progressively enlarging the infarct zone (IZ) and recruiting the remote zone (RZ) into the BZ, eventually leading to heart failure. The mechanisms underlying infarct extension remain unclear, but myocyte stretching has been suggested as the most likely cause. Using human patient-specific left-ventricular (LV) numerical simulations established from cardiac magnetic resonance imaging (MRI) of myocardial infarction (MI) patients, the correlation between infarct extension and regional mechanics abnormality was investigated by analysing the fibre stress-strain loops (FSSLs). FSSL abnormality was characterised using the directional regional external work (DREW) index, which measures FSSL area and loop direction. Sensitivity studies were also performed to investigate the effect of infarct stiffness on regional myocardial mechanics and potential for infarct extension. We found that infarct extension was correlated to severely abnormal FSSL in the form of counter-clockwise loop at the RZ close to the infarct, as indicated by negative DREW values. In regions demonstrating negative DREW values, we observed substantial fibre stretching in the isovolumic relaxation (IVR) phase accompanied by a reduced rate of systolic shortening. Such stretching in IVR phase in part of the RZ was due to its inability to withstand the high LV pressure that was still present and possibly caused by regional myocardial stiffness inhomogeneity. Further analysis revealed that the occurrence of severely abnormal FSSL due to IVR fibre stretching near the RZ-BZ boundary was due to a large amount of surrounding infarcted tissue, or an excessively stiff IZ.

Big data and new information technology: what cardiologists need to know

Technological progress in medicine is constantly garnering pace, requiring that physicians constantly update their knowledge. The new wave of technologies breaking through into clinical practice includes the following: a) mHealth, which allows constant monitoring of biological parameters, anytime, anyplace, of hundreds of patients at the same time; b) artificial intelligence, which, powered by new deep learning techniques, are starting to beat human experts at their own game: diagnosis by imaging or electrocardiography; c) 3-dimensional printing, which may lead to patient-specific prostheses; d) systems medicine, which has arisen from big data, and which will open the way to personalized medicine by bringing together genetic, epigenetic, environmental, clinical and social data into complex integral mathematical models to design highly personalized therapies. This state-of-the-art review aims to summarize in a single document the most recent and most important technological trends that are being applied to cardiology, and to provide an overall view that will allow readers to discern at a glance the direction of cardiology in the next few years.

Surrogate models based on machine learning methods for parameter estimation of left ventricular myocardium

A long-standing problem at the frontier of biomechanical studies is to develop fast methods capable of estimating material properties from clinical data. In this paper, we have studied three surrogate models based on machine learning (ML) methods for fast parameter estimation of left ventricular (LV) myocardium. We use three ML methods named K-nearest neighbour (KNN), XGBoost and multi-layer perceptron (MLP) to emulate the relationships between pressure and volume strains during the diastolic filling. Firstly, to train the surrogate models, a forward finite-element simulator of LV diastolic filling is used. Then the training data are projected in a low-dimensional parametrized space. Next, three ML models are trained to learn the relationships of pressure-volume and pressure-strain. Finally, an inverse parameter estimation problem is formulated by using those trained surrogate models. Our results show that the three ML models can learn the relationships of pressure-volume and pressure-strain very well, and the parameter inference using the surrogate models can be carried out in minutes. Estimated parameters from both the XGBoost and MLP models have much less uncertainties compared with the KNN model. Our results further suggest that the XGBoost model is better for predicting the LV diastolic dynamics and estimating passive parameters than other two surrogate models. Further studies are warranted to investigate how XGBoost can be used for emulating cardiac pump function in a multi-physics and multi-scale framework.

Biomechanics of infarcted left ventricle: a review of modelling

Mathematical modelling in biomechanics of infarcted left ventricle (LV) serves as an indispensable tool for remodelling mechanism exploration, LV biomechanical property estimation and therapy assessment after myocardial infarction (MI). However, a review of mathematical modelling after MI has not been seen in the literature so far. In the paper, a systematic review of mathematical models in biomechanics of infarcted LV was established. The models include comprehensive cardiovascular system model, essential LV pressure-volume and stress-stretch models, constitutive laws for passive myocardium and scars, tension models for active myocardium, collagen fibre orientation optimization models, fibroblast and collagen fibre growth/degradation models and integrated growth-electro-mechanical model after MI. The primary idea, unique characteristics and key equations of each model were identified and extracted. Discussions on the models were provided and followed research issues on them were addressed. Considerable improvements in the cardiovascular system model, LV aneurysm model, coupled agent-based models and integrated electro-mechanical-growth LV model are encouraged. Substantial attention should be paid to new constitutive laws with respect to stress-stretch curve and strain energy function for infarcted passive myocardium, collagen fibre orientation optimization in scar, cardiac rupture and tissue damage and viscoelastic effect post-MI in the future.

Fast parameter inference in a biomechanical model of the left ventricle by using statistical emulation

A central problem in biomechanical studies of personalized human left ventricular modelling is estimating the material properties and biophysical parameters from in vivo clinical measurements in a timeframe that is suitable for use within a clinic. Understanding these properties can provide insight into heart function or dysfunction and help to inform personalized medicine. However, finding a solution to the differential equations which mathematically describe the kinematics and dynamics of the myocardium through numerical integration can be computationally expensive. To circumvent this issue, we use the concept of emulation to infer the myocardial properties of a healthy volunteer in a viable clinical timeframe by using in vivo magnetic resonance image data. Emulation methods avoid computationally expensive simulations from the left ventricular model by replacing the biomechanical model, which is defined in terms of explicit partial differential equations, with a surrogate model inferred from simulations generated before the arrival of a patient, vastly improving computational efficiency at the clinic. We compare and contrast two emulation strategies: emulation of the computational model outputs and emulation of the loss between the observed patient data and the computational model outputs. These strategies are tested with two interpolation methods, as well as two loss functions. The best combination of methods is found by comparing the accuracy of parameter inference on simulated data for each combination. This combination, using the output emulation method, with local Gaussian process interpolation and the Euclidean loss function, provides accurate parameter inference in both simulated and clinical data, with a reduction in the computational cost of about three orders of magnitude compared with numerical integration of the differential equations by using finite element discretization techniques.