An Asymmetric Wall-Thickening Pattern Predicts Response to Cardiac Resynchronization Therapy

Advancing clinical translation of cardiac biomechanics models: a comprehensive review, applications and future pathways

Cardiac mechanics models are developed to represent a high level of detail, including refined anatomies, accurate cell mechanics models, and platforms to link microscale physiology to whole-organ function. However, cardiac biomechanics models still have limited clinical translation. In this review, we provide a picture of cardiac mechanics models, focusing on their clinical translation. We review the main experimental and clinical data used in cardiac models, as well as the steps followed in the literature to generate anatomical meshes ready for simulations. We describe the main models in active and passive mechanics and the different lumped parameter models to represent the circulatory system. Lastly, we provide a summary of the state-of-the-art in terms of ventricular, atrial, and four-chamber cardiac biomechanics models. We discuss the steps that may facilitate clinical translation of the biomechanics models we describe. A well-established software to simulate cardiac biomechanics is lacking, with all available platforms involving different levels of documentation, learning curves, accessibility, and cost. Furthermore, there is no regulatory framework that clearly outlines the verification and validation requirements a model has to satisfy in order to be reliably used in applications. Finally, better integration with increasingly rich clinical and/or experimental datasets as well as machine learning techniques to reduce computational costs might increase model reliability at feasible resources. Cardiac biomechanics models provide excellent opportunities to be integrated into clinical workflows, but more refinement and careful validation against clinical data are needed to improve their credibility. In addition, in each context of use, model complexity must be balanced with the associated high computational cost of running these models.

Predicting one-year left ventricular mass index regression following transcatheter aortic valve replacement in patients with severe aortic stenosis: A new era is coming

Aortic stenosis (AS) is the most common valvular heart disease in the western world, particularly worrisome with an ever-aging population wherein postoperative outcome for aortic valve replacement is strongly related to the timing of surgery in the natural course of disease. Yet, guidelines for therapy planning overlook insightful, quantified measures from medical imaging to educate clinical decisions. Herein, we leverage statistical shape analysis (SSA) techniques combined with customized machine learning methods to extract latent information from segmented left ventricle (LV) shapes. This enabled us to predict left ventricular mass index (LVMI) regression a year after transcatheter aortic valve replacement (TAVR). LVMI regression is an expected phenomena in patients undergone aortic valve replacement reported to be tightly correlated with survival one and five year after the intervention. In brief, LV geometries were extracted from medical images of a cohort of AS patients using deep learning tools, and then analyzed to create a set of statistical shape models (SSMs). Then, the supervised shape features were extracted to feed a support vector regression (SVR) model to predict the LVMI regression. The average accuracy of the predictions was validated against clinical measurements calculating root mean square error and R 2 score which yielded the satisfactory values of 0.28 and 0.67, respectively, on test data. Our work reveals the promising capability of advanced mathematical and bioinformatics approaches such as SSA and machine learning to improve medical output prediction and treatment planning.

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.

The role of cardiac magnetic resonance in identifying appropriate candidates for cardiac resynchronization therapy - a systematic review of the literature

Despite the strict indications for cardiac resynchronization therapy (CRT) implantation, a significant proportion of patients will fail to adequately respond to the treatment. This systematic review aims to present the existing evidence about the role of cardiac magnetic resonance (CMR) in identifying patients who are likely to respond better to the CRT. A systematic search in the MedLine database and Cochrane Library from their inception to August 2021 was performed, without any limitations, by two independent investigators. We considered eligible observational studies or randomized clinical trials (RCTs) that enrolled patients > 18 years old with heart failure (HF) of ischaemic or non-ischaemic aetiology and provided data about the association of baseline CMR variables with clinical or echocardiographic response to CRT for at least 3 months. This systematic review was performed in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA Statement). Following our search strategy, 47 studies were finally included in our review. CMR appears to have an additive role in identifying the subgroup of patients who will respond better to CRT. Specifically, the presence and the extent of myocardial scar were associated with increased non-response rates, while those with no scar respond better. Furthermore, existing data show that scar location can be associated with CRT response rates. CMR-derived markers of mechanical desynchrony can also be used as predictors of CRT response. CMR data can be used to optimize the position of the left ventricular lead during the CRT implantation procedure. Specifically, positioning the left ventricular lead in a branch of the coronary sinus that feeds an area with transmural scar was associated with poorer response to CRT. CMR can be used as a non-invasive optimization tool to identify patients who are more likely to achieve better clinical and echocardiographic response following CRT implantation.

A Poroelastic Approach for Modelling Myocardial Oedema in Acute Myocarditis

Myocarditis is a general set of mechanisms that manifest themselves into the inflammation of the heart muscle. In 2017, more than 3 million people were affected by this disease worldwide, causing about 47,000 deaths. Many aspects of the origin of this disease are well known, but several important questions regarding the disease remain open. One of them is why some patients develop a significantly localised inflammation while others develop a much more diffuse inflammation, reaching across large portions of the heart. Furthermore, the specific role of the pathogenic agent that causes inflammation as well as the interaction with the immune system in the progression of the disease are still under discussion. Providing answers to these crucial questions can have an important impact on patient treatment. In this scenario, computational methods can aid specialists to understand better the relationships between pathogens and the immune system and elucidate why some patients develop diffuse myocarditis. This paper alters a recently developed model to study the myocardial oedema formation in acute infectious myocarditis. The model describes the finite deformation regime using partial differential equations to represent tissue displacement, fluid pressure, fluid phase, and the concentrations of pathogens and leukocytes. A sensitivity analysis was performed to understand better the influence of the most relevant model parameters on the disease dynamics. The results showed that the poroelastic model could reproduce local and diffuse myocarditis dynamics in simplified and complex geometrical domains.

Septal curvature as a robust and reproducible marker for basal septal hypertrophy

BACKGROUND

Basal septal hypertrophy (BSH) is an asymmetric, localized thickening of the upper interventricular septum and constitutes a marker of an early remodelling in patients with hypertension. This morphological trait has been extensively researched because of its prevalence in hypertension, yet its clinical and prognostic value for individual patients remains undetermined. One of the reasons is the lack of a reliable and reproducible metric to quantify the presence and the extent of BSH. This article proposes the use of the curvature of the left ventricular endocardium as a robust feature for BSH characterization, and as an objective criterion to quantify current subjective 'visual assessment' of the presence of sigmoidal septum. The proposed marker, called average septal curvature, is defined as the inverse of the radius adjacent to each point of the endocardial contour along the basal and mid inferoseptal segments of the left ventricle.

METHOD

Robustness and reproducibility were assessed on a cohort of 220 patients, including 161 hypertensive patients (32 with BSH) and 59 healthy controls.

RESULTS

The results show that compared with the conventionally used wall thickness metrics, the new marker is more reproducible (relative standard deviation of errors of 7 vs. 13%, and 8 vs. 38% for intra-observer and inter-observer variability, respectively) and better correlates to the functional parameters related to BSH, with main difference (absolute rank correlation 0.417 vs. 0.341) in local deformation changes assessed by longitudinal strain.

CONCLUSION

Average septal curvature is a more precisely defined and reproducible metric than thickness ratios, it can be fully automated, and better infers the functional remodelling related to hypertension.

Linking statistical shape models and simulated function in the healthy adult human heart

Cardiac anatomy plays a crucial role in determining cardiac function. However, there is a poor understanding of how specific and localised anatomical changes affect different cardiac functional outputs. In this work, we test the hypothesis that in a statistical shape model (SSM), the modes that are most relevant for describing anatomy are also most important for determining the output of cardiac electromechanics simulations. We made patient-specific four-chamber heart meshes (n = 20) from cardiac CT images in asymptomatic subjects and created a SSM from 19 cases. Nine modes captured 90% of the anatomical variation in the SSM. Functional simulation outputs correlated best with modes 2, 3 and 9 on average (R = 0.49 ± 0.17, 0.37 ± 0.23 and 0.34 ± 0.17 respectively). We performed a global sensitivity analysis to identify the different modes responsible for different simulated electrical and mechanical measures of cardiac function. Modes 2 and 9 were the most important for determining simulated left ventricular mechanics and pressure-derived phenotypes. Mode 2 explained 28.56 ± 16.48% and 25.5 ± 20.85, and mode 9 explained 12.1 ± 8.74% and 13.54 ± 16.91% of the variances of mechanics and pressure-derived phenotypes, respectively. Electrophysiological biomarkers were explained by the interaction of 3 ± 1 modes. In the healthy adult human heart, shape modes that explain large portions of anatomical variance do not explain equivalent levels of electromechanical functional variation. As a result, in cardiac models, representing patient anatomy using a limited number of modes of anatomical variation can cause a loss in accuracy of simulated electromechanical function.

The 'Digital Twin' to enable the vision of precision cardiology

Providing therapies tailored to each patient is the vision of precision medicine, enabled by the increasing ability to capture extensive data about individual patients. In this position paper, we argue that the second enabling pillar towards this vision is the increasing power of computers and algorithms to learn, reason, and build the 'digital twin' of a patient. Computational models are boosting the capacity to draw diagnosis and prognosis, and future treatments will be tailored not only to current health status and data, but also to an accurate projection of the pathways to restore health by model predictions. The early steps of the digital twin in the area of cardiovascular medicine are reviewed in this article, together with a discussion of the challenges and opportunities ahead. We emphasize the synergies between mechanistic and statistical models in accelerating cardiovascular research and enabling the vision of precision medicine.

An orthotropic electro-viscoelastic model for the heart with stress-assisted diffusion

We propose and analyse the properties of a new class of models for the electromechanics of cardiac tissue. The set of governing equations consists of nonlinear elasticity using a viscoelastic and orthotropic exponential constitutive law, for both active stress and active strain formulations of active mechanics, coupled with a four-variable phenomenological model for human cardiac cell electrophysiology, which produces an accurate description of the action potential. The conductivities in the model of electric propagation are modified according to stress, inducing an additional degree of nonlinearity and anisotropy in the coupling mechanisms, and the activation model assumes a simplified stretch-calcium interaction generating active tension or active strain. The influence of the new terms in the electromechanical model is evaluated through a sensitivity analysis, and we provide numerical validation through a set of computational tests using a novel mixed-primal finite element scheme.

Sex-Dependent QRS Guidelines for Cardiac Resynchronization Therapy Using Computer Model Predictions

Cardiac resynchronization therapy (CRT) is an important treatment for heart failure. Low female enrollment in clinical trials means that current CRT guidelines may be biased toward males. However, females have higher response rates at lower QRS duration (QRSd) thresholds. Sex differences in the left ventricle (LV) size could provide an explanation for the improved female response at lower QRSd. We aimed to test if sex differences in CRT response at lower QRSd thresholds are explained by differences in LV size and hence predict sex-specific guidelines for CRT. We investigated the effect that LV size sex difference has on QRSd between male and females in 1093 healthy individuals and 50 CRT patients using electrophysiological computer models of the heart. Simulations on the healthy mean shape models show that LV size sex difference can account for 50–100% of the sex difference in baseline QRSd in healthy individuals. In the CRT patient cohort, model simulations predicted female-specific guidelines for CRT, which were 9–13 ms lower than current guidelines. Sex differences in the LV size are able to account for a significant proportion of the sex difference in QRSd and provide a mechanistic explanation for the sex difference in CRT response. Simulations accounting for the smaller LV size in female CRT patients predict 9–13 ms lower QRSd thresholds for female CRT guidelines.

Cardiac Resynchronization Therapy Optimization: A Comprehensive Approach

Since the first report on biventricular pacing in 1994, cardiac resynchronization therapy (CRT) has become standard for patients with advanced heart failure (HF) and ventricular conduction delay. CRT improves myocardial function by resynchronizing myocardial contraction, which results in reverse left ventricular remodeling and improves symptoms and clinical outcomes. Despite the accelerated development of CRT device technology and its increased application in treating HF patients, almost one-third of these patients do not respond to the therapy or gain any clinical benefit from device implantation. Over the last decade, multiple cardiac imaging modalities have provided a deeper understanding of myocardial pathophysiology, thereby improving HF treatment management. However, the optimal strategy for improving the CRT response remains debatable. This article provides an updated overview of the electropathophysiology of myocardial dysfunction in ventricular conduction delay and the diagnostic approaches involving the use of multiple modalities.

Independent Left Ventricular Morphometric Atlases Show Consistent Relationships with Cardiovascular Risk Factors: A UK Biobank Study

Left ventricular (LV) mass and volume are important indicators of clinical and pre-clinical disease processes. However, much of the shape information present in modern imaging examinations is currently ignored. Morphometric atlases enable precise quantification of shape and function, but there has been no objective comparison of different atlases in the same cohort. We compared two independent LV atlases using MRI scans of 4547 UK Biobank participants: (i) a volume atlas derived by automatic non-rigid registration of image volumes to a common template, and (ii) a surface atlas derived from manually drawn epicardial and endocardial surface contours. The strength of associations between atlas principal components and cardiovascular risk factors (smoking, diabetes, high blood pressure, high cholesterol and angina) were quantified with logistic regression models and five-fold cross validation, using area under the ROC curve (AUC) and Akaike Information Criterion (AIC) metrics. Both atlases exhibited similar principal components, showed similar relationships with risk factors, and had stronger associations (higher AUC and lower AIC) than a reference model based on LV mass and volume, for all risk factors (DeLong p < 0.05). Morphometric variations associated with each risk factor could be quantified and visualized and were similar between atlases. UK Biobank LV shape atlases are robust to construction method and show stronger relationships with cardiovascular risk factors than mass and volume.