In atrial fibrillation (AF), atrial biomechanics are altered, reducing atrial movement. It remains unclear whether these changes are due to altered anatomy, myocardial stiffness, or constraints from surrounding structures. Understanding the causes of changed atrial deformation in AF could enhance tissue characterization and inform AF diagnosis, stratification, and treatment. We created patient-specific anatomical models of the left atrium (LA) from CT images. Passive LA biomechanics were simulated using finite deformation continuum mechanics equations. LA stiffness was represented by the Guccione material law, where α scaled the anisotropic stiffness parameters. Regional passive stiffness parameters were calibrated to peak regional deformation during the reservoir phase and validated against deformation transients derived from retrospective gated CT images during the reservoir and conduit phase. Physiological LA deformation varies regionally, with the roof deforming significantly less than other regions during the reservoir phase. The fitted model matched peak patient deformations globally and regionally with an average error of [Formula: see text] mm over our cohort. We compared deformation transients through the reservoir and conduit phases and found that the simulated deformation transients were within an average of [Formula: see text] mm per unit time of the CT-derived deformation transients. Regional stiffness varied across the atria with average α values of 1.8, 1.6, 2.2, 1.6 and 2.1 across the cohort in the anterior, posterior, septum, lateral and roof regions respectively. Using mixed effect models, we found no correlation between regional patient LA deformation and regional estimates of wall thickness or regional volumes of epicardial adipose tissue. We found a significant correlation between regionally calibrated stiffness and CT-derived LA biomechanics (p = 0.023). We have shown that regional heterogeneity in stiffness contributes to regional LA biomechanics, while anatomical features appeared less important. These findings provide insight into the underlying causes of altered LA biomechanics in AF.

Quantifying uncertainty of predictions has been identified as one way to develop more trustworthy artificial intelligence (AI) models beyond conventional reporting of performance metrics. When considering their role in a clinical decision support setting, AI classification models should ideally avoid confident wrong predictions and maximise the confidence of correct predictions. Models that do this are said to be well calibrated with regard to confidence. However, relatively little attention has been paid to how to improve calibration when training these models, i.e. to make the training strategy uncertainty-aware. In this work we: (i) evaluate three novel uncertainty-aware training strategies with regard to a range of accuracy and calibration performance measures, comparing against two state-of-the-art approaches, (ii) quantify the data (aleatoric) and model (epistemic) uncertainty of all models and (iii) evaluate the impact of using a model calibration measure for model selection in uncertainty-aware training, in contrast to the normal accuracy-based measures. We perform our analysis using two different clinical applications: cardiac resynchronisation therapy (CRT) response prediction and coronary artery disease (CAD) diagnosis from cardiac magnetic resonance (CMR) images. The best-performing model in terms of both classification accuracy and the most common calibration measure, expected calibration error (ECE) was the Confidence Weight method, a novel approach that weights the loss of samples to explicitly penalise confident incorrect predictions. The method reduced the ECE by 17% for CRT response prediction and by 22% for CAD diagnosis when compared to a baseline classifier in which no uncertainty-aware strategy was included. In both applications, as well as reducing the ECE there was a slight increase in accuracy from 69% to 70% and 70% to 72% for CRT response prediction and CAD diagnosis respectively. However, our analysis showed a lack of consistency in terms of optimal models when using different calibration measures. This indicates the need for careful consideration of performance metrics when training and selecting models for complex high risk applications in healthcare.

Cardiac pump function arises from a series of highly orchestrated events across multiple scales. Computational electromechanics can encode these events in physics-constrained models. However, the large number of parameters in these models has made the systematic study of the link between cellular, tissue, and organ scale parameters to whole heart physiology challenging. A patient-specific anatomical heart model, or digital twin, was created. Cellular ionic dynamics and contraction were simulated with the Courtemanche-Land and the ToR-ORd-Land models for the atria and the ventricles, respectively. Whole heart contraction was coupled with the circulatory system, simulated with CircAdapt, while accounting for the effect of the pericardium on cardiac motion. The four-chamber electromechanics framework resulted in 117 parameters of interest. The model was broken into five hierarchical sub-models: tissue electrophysiology, ToR-ORd-Land model, Courtemanche-Land model, passive mechanics and CircAdapt. For each sub-model, we trained Gaussian processes emulators (GPEs) that were then used to perform a global sensitivity analysis (GSA) to retain parameters explaining 90% of the total sensitivity for subsequent analysis. We identified 45 out of 117 parameters that were important for whole heart function. We performed a GSA over these 45 parameters and identified the systemic and pulmonary peripheral resistance as being critical parameters for a wide range of volumetric and hemodynamic cardiac indexes across all four chambers. We have shown that GPEs provide a robust method for mapping between cellular properties and clinical measurements. This could be applied to identify parameters that can be calibrated in patient-specific models or digital twins, and to link cellular function to clinical indexes.

Catheter ablation is currently the only curative treatment for scar-related ventricular tachycardias (VTs). However, not only are ablation procedures long, with relatively high risk, but success rates are punitively low, with frequent VT recurrence. Personalized in-silico approaches have the opportunity to address these limitations. However, state-of-the-art reaction diffusion (R-D) simulations of VT induction and subsequent circuits used for in-silico ablation target identification require long execution times, along with vast computational resources, which are incompatible with the clinical workflow. Here, we present the Virtual Induction and Treatment of Arrhythmias (VITA), a novel, rapid and fully automated computational approach that uses reaction-Eikonal methodology to induce VT and identify subsequent ablation targets. The rationale for VITA is based on finding isosurfaces associated with an activation wavefront that splits in the ventricles due to the presence of an isolated isthmus of conduction within the scar; once identified, each isthmus may be assessed for their vulnerability to sustain a reentrant circuit, and the corresponding exit site automatically identified for potential ablation targeting. VITA was tested on a virtual cohort of 7 post-infarcted porcine hearts and the results compared to R-D simulations. Using only a standard desktop machine, VITA could detect all scar-related VTs, simulating activation time maps and ECGs (for clinical comparison) as well as computing ablation targets in 48 minutes. The comparable VTs probed by the R-D simulations took 68.5 hours on 256 cores of high-performance computing infrastructure. The set of lesions computed by VITA was shown to render the ventricular model VT-free. VITA could be used in near real-time as a complementary modality aiding in clinical decision-making in the treatment of post-infarction VTs.

We present a novel multimodal deep learning framework for cardiac resynchronisation therapy (CRT) response prediction from 2D echocardiography and cardiac magnetic resonance (CMR) data. The proposed method first uses the 'nnU-Net' segmentation model to extract segmentations of the heart over the full cardiac cycle from the two modalities. Next, a multimodal deep learning classifier is used for CRT response prediction, which combines the latent spaces of the segmentation models of the two modalities. At test time, this framework can be used with 2D echocardiography data only, whilst taking advantage of the implicit relationship between CMR and echocardiography features learnt from the model. We evaluate our pipeline on a cohort of 50 CRT patients for whom paired echocardiography/CMR data were available, and results show that the proposed multimodal classifier results in a statistically significant improvement in accuracy compared to the baseline approach that uses only 2D echocardiography data. The combination of multimodal data enables CRT response to be predicted with 77.38% accuracy (83.33% sensitivity and 71.43% specificity), which is comparable with the current state-of-the-art in machine learning-based CRT response prediction. Our work represents the first multimodal deep learning approach for CRT response prediction.

BACKGROUND: Cardiac Resynchronization Therapy (CRT) is one of the few effective treatments for heart failure patients with ventricular dyssynchrony. The pacing location of the left ventricle is indicated as a determinant of CRT outcome. OBJECTIVE: Patient specific computational models allow the activation pattern following CRT implant to be predicted and this may be used to optimize CRT lead placement. METHODS: In this study, the effects of heterogeneous cardiac substrate (scar, fast endocardial conduction, slow septal conduction, functional block) on accurately predicting the electrical activation of the LV epicardium were tested to determine the minimal detail required to create a rule based model of cardiac electrophysiology. Non-invasive clinical data (CT or CMR images and 12 lead ECG) from eighteen patients from two centers were used to investigate the models. RESULTS: Validation with invasive electro-anatomical mapping data identified that computer models with fast endocardial conduction were able to predict the electrical activation with a mean distance errors of 9.2 ± 0.5 mm (CMR data) or (CT data) 7.5 ± 0.7 mm. CONCLUSION: This study identified a simple rule-based fast endocardial conduction model, built using non-invasive clinical data that can be used to rapidly and robustly predict the electrical activation of the heart. Pre-procedural prediction of the latest electrically activating region to identify the optimal LV pacing site could potentially be a useful clinical planning tool for CRT procedures.

This paper presents a morphological analysis of fibrotic scarring in non-ischemic dilated cardiomyopathy, and its relationship to electrical instabilities which underlie reentrant arrhythmias.Two dimensional electrophysiological simulation models were constructed from a set of 699 late gadolinium enhanced cardiac magnetic resonance images originating from 157 patients.Areas of late gadolinium enhancement (LGE) in each image were assigned one of 10 possible microstructures, which modelled the details of fibrotic scarring an order of magnitude below the MRI scan resolution.A simulated programmed electrical stimulation protocol tested each model for the possibility of generating either a transmural block or a transmural reentry.The outcomes of the simulations were compared against morphological LGE features extracted from the images.Models which blocked or reentered, grouped by microstructure, were significantly different from one another in myocardial-LGE interface length, number of components and entropy, but not in relative area and transmurality.With an unknown microstructure, transmurality alone was the best predictor of block, whereas a combination of interface length, transmurality and number of components was the best predictor of reentry in linear discriminant analysis. Author summaryNon-ischemic dilated cardiomyopathy is a disease in which the lower left chamber of the heart is abnormally large.The cause of the disease can be anything that is not a loss of blood supply to the heart.Many patients with non-ischemic dilated cardiomyopathy have scars in their hearts which can be detected with magnetic resonance imaging.These scars are thought to disrupt the flow of electricity through the heart and cause deadly rhythm

Abnormal cardiac motion can indicate different forms of disease, which can manifest at different spatial scales in the myocardium. Many studies have sought to characterise particular motion abnormalities associated with specific diseases, and to utilise motion information to improve diagnoses. However, the importance of spatial scale in the analysis of cardiac deformation has not been extensively investigated. We build on recent work on the analysis of myocardial strains at different spatial scales using a cardiac motion atlas to find the optimal scales for estimating different cardiac biomarkers. We apply a multi-scale strain analysis to a 43 patient cohort of cardiac resynchronisation therapy (CRT) patients using tagged magnetic resonance imaging data for (1) predicting response to CRT, (2) identifying septal flash, (3) estimating QRS duration, and (4) identifying the presence of ischaemia. A repeated, stratified cross-validation is used to demonstrate the importance of spatial scale in our analysis, revealing different optimal spatial scales for the estimation of different biomarkers.

We present a framework for combining a cardiac motion atlas with non-motion data. The atlas represents cardiac cycle motion across a number of subjects in a common space based on rich motion descriptors capturing 3D displacement, velocity, strain and strain rate. The non-motion data are derived from a variety of sources such as imaging, electrocardiogram (ECG) and clinical reports. Once in the atlas space, we apply a novel supervised learning approach based on random projections and ensemble learning to learn the relationship between the atlas data and some desired clinical output. We apply our framework to the problem of predicting response to Cardiac Resynchronisation Therapy (CRT). Using a cohort of 34 patients selected for CRT using conventional criteria, results show that the combination of motion and non-motion data enables CRT response to be predicted with 91.2% accuracy (100% sensitivity and 62.5% specificity), which compares favourably with the current state-of-the-art in CRT response prediction.

Patients with drug-refractory heart failure can greatly benefit from cardiac resynchronization therapy (CRT). A CRT device can resynchronize the contractions of the left ventricle (LV) leading to reduced mortality. Unfortunately, 30%-50% of patients do not respond to treatment when assessed by objective criteria such as cardiac remodeling. A significant contributing factor is the suboptimal placement of the LV lead. It has been shown that placing this lead away from scar and at the point of latest mechanical activation can improve response rates. This paper presents a comprehensive and highly automated system that uses scar and mechanical activation to plan and guide CRT procedures. Standard clinical preoperative magnetic resonance imaging is used to extract scar and mechanical activation information. The data are registered to a single 3-D coordinate system and visualized in novel 2-D and 3-D American Heart Association plots enabling the clinician to select target segments. During the procedure, the planning information is overlaid onto live fluoroscopic images to guide lead deployment. The proposed platform has been used during 14 CRT procedures and validated on synthetic, phantom, volunteer, and patient data.

An unresolved issue in patients with diastolic dysfunction is that the estimation of myocardial stiffness cannot be decoupled from diastolic residual active tension (AT) because of the impaired ventricular relaxation during diastole. To address this problem, this paper presents a method for estimating diastolic mechanical parameters of the left ventricle (LV) from cine and tagged MRI measurements and LV cavity pressure recordings, separating the passive myocardial constitutive properties and diastolic residual AT. Dynamic C1-continuous meshes are automatically built from the anatomy and deformation captured from dynamic MRI sequences. Diastolic deformation is simulated using a mechanical model that combines passive and active material properties. The problem of non-uniqueness of constitutive parameter estimation using the well known Guccione law is characterized by reformulation of this law. Using this reformulated form, and by constraining the constitutive parameters to be constant across time points during diastole, we separate the effects of passive constitutive properties and the residual AT during diastolic relaxation. Finally, the method is applied to two clinical cases and one control, demonstrating that increased residual AT during diastole provides a potential novel index for delineating healthy and pathological cases.

We present an adaptive algorithm which uses a fast electrophysiological (EP) model to estimate apparent electrical conductivity and local conduction velocity from noncontact mapping of the endocardial surface potential. Development of such functional imaging revealing hidden parameters of the heart can be instrumental for improved diagnosis and planning of therapy for cardiac arrhythmia and heart failure, for example during procedures such as radio-frequency ablation and cardiac resynchronisation therapy. The proposed model is validated on synthetic data and applied to clinical data derived using hybrid X-ray/magnetic resonance imaging. We demonstrate a qualitative match between the estimated conductivity parameter and pathology locations in the human left ventricle. We also present a proof of concept for an electrophysiological model which utilizes the estimated apparent conductivity parameter to simulate the effect of pacing different ventricular sites. This approach opens up possibilities to directly integrate modelling in the cardiac EP laboratory.