Cobiveco: Consistent biventricular coordinates for precise and intuitive description of position in the heart - with MATLAB implementation

Improved algorithm for generating evenly-spaced streamlines from an orientation field on a triangulated surface

BACKGROUND

Vector fields such as cardiac fiber orientation can be visualized on a surface using streamlines. The application of evenly-spaced streamline generation to the construction of interconnected cable structure for cardiac propagation models has more stringent requirements imperfectly fulfilled by current algorithms.

METHOD

We developed an open-source C++/python package for the placement of evenly-spaced streamlines on a triangulated surface. The new algorithm improves upon previous works by more accurately handling streamline extremities, U-turns and limit cycles, by providing stronger geometrical guarantees on inter-streamline minimal distance, particularly when a high density of streamlines (up to 10μm spacing) is desired, and by making a more efficient parallel implementation available. The approach requires finding intersections between geometrical capsules and triangles to update an occupancy mask defined on the triangles. This enables fast streamline integration from thousands of seed points to identify optimal streamline placement.

RESULTS

The algorithm was assessed qualitatively on different left atrial models of fiber orientation with varying mesh resolutions (up to 375k triangles) and quantitatively by measuring streamline lengths and distribution of inter-streamline minimal distance. The complexity and the computational performance of the algorithm were studied as a function of streamline spacing in relation to triangular mesh resolution.

CONCLUSION

More accurate geometrical computations, attention to details and fine-tuning led to an algorithm more amenable to applications that require precise positioning of streamlines.

Digital twinning of the human ventricular activation sequence to Clinical 12-lead ECGs and magnetic resonance imaging using realistic Purkinje networks for in silico clinical trials

Cardiac in silico clinical trials can virtually assess the safety and efficacy of therapies using human-based modelling and simulation. These technologies can provide mechanistic explanations for clinically observed pathological behaviour. Designing virtual cohorts for in silico trials requires exploiting clinical data to capture the physiological variability in the human population. The clinical characterisation of ventricular activation and the Purkinje network is challenging, especially non-invasively. Our study aims to present a novel digital twinning pipeline that can efficiently generate and integrate Purkinje networks into human multiscale biventricular models based on subject-specific clinical 12-lead electrocardiogram and magnetic resonance recordings. Essential novel features of the pipeline are the human-based Purkinje network generation method, personalisation considering ECG R wave progression as well as QRS morphology, and translation from reduced-order Eikonal models to equivalent biophysically-detailed monodomain ones. We demonstrate ECG simulations in line with clinical data with clinical image-based multiscale models with Purkinje in four control subjects and two hypertrophic cardiomyopathy patients (simulated and clinical QRS complexes with Pearson's correlation coefficients > 0.7). Our methods also considered possible differences in the density of Purkinje myocardial junctions in the Eikonal-based inference as regional conduction velocities. These differences translated into regional coupling effects between Purkinje and myocardial models in the monodomain formulation. In summary, we demonstrate a digital twin pipeline enabling simulations yielding clinically consistent ECGs with clinical CMR image-based biventricular multiscale models, including personalised Purkinje in healthy and cardiac disease conditions.

A universal biventricular coordinate system incorporating valve annuli: Validation in congenital heart disease

Universal coordinate systems have been proposed to facilitate anatomic registration between three-dimensional images, data and models of the ventricles of the heart. However, current universal ventricular coordinate systems do not account for the outflow tracts and valve annuli where the anatomy is complex. Here we propose an extension to the 'Cobiveco' biventricular coordinate system that also accounts for the intervalvular bridges of the base and provides a tool for anatomically consistent registration between widely varying biventricular shapes. CobivecoX uses a novel algorithm to separate intervalvular bridges and assign new coordinates, including an inflow-outflow coordinate, to describe local positions in these regions uniquely and consistently. Anatomic consistency of registration was validated using curated three-dimensional biventricular shape models derived from cardiac MRI measurements in normal hearts and hearts from patients with congenital heart diseases. This new method allows the advantages of universal cardiac coordinates to be used for three-dimensional ventricular imaging data and models that include the left and right ventricular outflow tracts and valve annuli.

Non-invasive localization of the ventricular excitation origin without patient-specific geometries using deep learning

Cardiovascular diseases account for 17 million deaths per year worldwide. Of these, 25% are categorized as sudden cardiac death, which can be related to ventricular tachycardia (VT). This type of arrhythmia can be caused by focal activation sources outside the sinus node. Catheter ablation of these foci is a curative treatment in order to inactivate the abnormal triggering activity. However, the localization procedure is usually time-consuming and requires an invasive procedure in the catheter lab. To facilitate and expedite the treatment, we present two novel localization support techniques based on convolutional neural networks (CNNs) that address these clinical needs. In contrast to existing methods, our approaches were designed to be independent of the patient-specific geometry and directly applicable to surface ECG signals, while also delivering a binary transmural position. Moreover, one of the method's outputs can be interpreted as several ranked solutions. The CNNs were trained on a dataset containing only simulated data and evaluated both on simulated test data and clinical data. On a novel large and open simulated dataset, the median test error was below 3 mm. The median localization error on the unseen clinical data ranged from 32 mm to 41 mm without optimizing the pre-processing and CNN to the clinical data. Interpreting the output of one of the approaches as ranked solutions, the best median error of the top-3 solutions decreased to 20 mm on the clinical data. The transmural position was correctly detected in up to 82% of all clinical cases. These results demonstrate a proof of principle to utilize CNNs to localize the activation source without the intrinsic need for patient-specific geometrical information. Furthermore, providing multiple solutions can assist physicians in identifying the true activation source amongst more than one possible location. With further optimization to clinical data, these methods have high potential to accelerate clinical interventions, replace certain steps within these procedures and consequently reduce procedural risk and improve VT patient outcomes.

Multi-modality cardiac image computing: A survey

Multi-modality cardiac imaging plays a key role in the management of patients with cardiovascular diseases. It allows a combination of complementary anatomical, morphological and functional information, increases diagnosis accuracy, and improves the efficacy of cardiovascular interventions and clinical outcomes. Fully-automated processing and quantitative analysis of multi-modality cardiac images could have a direct impact on clinical research and evidence-based patient management. However, these require overcoming significant challenges including inter-modality misalignment and finding optimal methods to integrate information from different modalities. This paper aims to provide a comprehensive review of multi-modality imaging in cardiology, the computing methods, the validation strategies, the related clinical workflows and future perspectives. For the computing methodologies, we have a favored focus on the three tasks, i.e., registration, fusion and segmentation, which generally involve multi-modality imaging data, either combining information from different modalities or transferring information across modalities. The review highlights that multi-modality cardiac imaging data has the potential of wide applicability in the clinic, such as trans-aortic valve implantation guidance, myocardial viability assessment, and catheter ablation therapy and its patient selection. Nevertheless, many challenges remain unsolved, such as missing modality, modality selection, combination of imaging and non-imaging data, and uniform analysis and representation of different modalities. There is also work to do in defining how the well-developed techniques fit in clinical workflows and how much additional and relevant information they introduce. These problems are likely to continue to be an active field of research and the questions to be answered in the future.

Comparison of novel ventricular pacing strategies using an electro-mechanical simulation platform

AIMS

Focus of pacemaker therapy is shifting from right ventricular (RV) apex pacing (RVAP) and biventricular pacing (BiVP) to conduction system pacing. Direct comparison between the different pacing modalities and their consequences to cardiac pump function is difficult, due to the practical implications and confounding variables. Computational modelling and simulation provide the opportunity to compare electrical, mechanical, and haemodynamic consequences in the same virtual heart.

METHODS AND RESULTS

Using the same single cardiac geometry, electrical activation maps following the different pacing strategies were calculated using an Eikonal model on a three-dimensional geometry, which were then used as input for a lumped mechanical and haemodynamic model (CircAdapt). We then compared simulated strain, regional myocardial work, and haemodynamic function for each pacing strategy. Selective His-bundle pacing (HBP) best replicated physiological electrical activation and led to the most homogeneous mechanical behaviour. Selective left bundle branch (LBB) pacing led to good left ventricular (LV) function but significantly increased RV load. RV activation times were reduced in non-selective LBB pacing (nsLBBP), reducing RV load but increasing heterogeneity in LV contraction. LV septal pacing led to a slower LV and more heterogeneous LV activation than nsLBBP, while RV activation was similar. BiVP led to a synchronous LV-RV, but resulted in a heterogeneous contraction. RVAP led to the slowest and most heterogeneous contraction. Haemodynamic differences were small compared to differences in local wall behaviour.

CONCLUSION

Using a computational modelling framework, we investigated the mechanical and haemodynamic outcome of the prevailing pacing strategies in hearts with normal electrical and mechanical function. For this class of patients, nsLBBP was the best compromise between LV and RV function if HBP is not possible.

The Impact of Standard Ablation Strategies for Atrial Fibrillation on Cardiovascular Performance in a Four-Chamber Heart Model

PURPOSE

Atrial fibrillation is one of the most frequent cardiac arrhythmias in the industrialized world and ablation therapy is the method of choice for many patients. However, ablation scars alter the electrophysiological activation and the mechanical behavior of the affected atria. Different ablation strategies with the aim to terminate atrial fibrillation and prevent its recurrence exist but their impact on the performance of the heart is often neglected.

METHODS

In this work, we present a simulation study analyzing five commonly used ablation scar patterns and their combinations in the left atrium regarding their impact on the pumping function of the heart using an electromechanical whole-heart model. We analyzed how the altered atrial activation and increased stiffness due to the ablation scars affect atrial as well as ventricular contraction and relaxation.

RESULTS

We found that systolic and diastolic function of the left atrium is impaired by ablation scars and that the reduction of atrial stroke volume of up to 11.43% depends linearly on the amount of inactivated tissue. Consequently, the end-diastolic volume of the left ventricle, and thus stroke volume, was reduced by up to 1.4 and 1.8%, respectively. During ventricular systole, left atrial pressure was increased by up to 20% due to changes in the atrial activation sequence and the stiffening of scar tissue.

CONCLUSION

This study provides biomechanical evidence that atrial ablation has acute effects not only on atrial contraction but also on ventricular performance. Therefore, the position and extent of ablation scars is not only important for the termination of arrhythmias but is also determining long-term pumping efficiency. If confirmed in larger cohorts, these results have the potential to help tailoring ablation strategies towards minimal global cardiovascular impairment.

Functional Characterization of a Spectrum of Novel Romano-Ward Syndrome KCNQ1 Variants

The KCNQ1 gene encodes the α-subunit of the cardiac voltage-gated potassium (Kv) channel KCNQ1, also denoted as Kv7.1 or KvLQT1. The channel assembles with the ß-subunit KCNE1, also known as minK, to generate the slowly activating cardiac delayed rectifier current I_Ks, a key regulator of the heart rate dependent adaptation of the cardiac action potential duration (APD). Loss-of-function variants in KCNQ1 cause the congenital Long QT1 (LQT1) syndrome, characterized by delayed cardiac repolarization and a QT interval prolongation in the surface electrocardiogram (ECG). Autosomal dominant loss-of-function variants in KCNQ1 result in the LQT syndrome called Romano-Ward syndrome (RWS), while autosomal recessive variants affecting function, lead to Jervell and Lange-Nielsen syndrome (JLNS), associated with deafness. The aim of this study was the characterization of novel KCNQ1 variants identified in patients with RWS to widen the spectrum of known LQT1 variants, and improve the interpretation of the clinical relevance of variants in the KCNQ1 gene. We functionally characterized nine human KCNQ1 variants using the voltage-clamp technique in Xenopus laevis oocytes, from which we report seven novel variants. The functional data was taken as input to model surface ECGs, to subsequently compare the functional changes with the clinically observed QTc times, allowing a further interpretation of the severity of the different LQTS variants. We found that the electrophysiological properties of the variants correlate with the severity of the clinically diagnosed phenotype in most cases, however, not in all. Electrophysiological studies combined with in silico modelling approaches are valuable components for the interpretation of the pathogenicity of KCNQ1 variants, but assessing the clinical severity demands the consideration of other factors that are included, for example in the Schwartz score.

Dyssynchronous Left Ventricular Activation is Insufficient for the Breakdown of Wringing Rotation

Cardiac resynchronization therapy is a valuable tool to restore left ventricular function in patients experiencing dyssynchronous ventricular activation. However, the non-responder rate is still as high as 40%. Recent studies suggest that left ventricular torsion or specifically the lack thereof might be a good predictor for the response of cardiac resynchronization therapy. Since left ventricular torsion is governed by the muscle fiber orientation and the heterogeneous electromechanical activation of the myocardium, understanding the relation between these components and the ability to measure them is vital. To analyze if locally altered electromechanical activation in heart failure patients affects left ventricular torsion, we conducted a simulation study on 27 personalized left ventricular models. Electroanatomical maps and late gadolinium enhanced magnetic resonance imaging data informed our in-silico model cohort. The angle of rotation was evaluated in every material point of the model and averaged values were used to classify the rotation as clockwise or counterclockwise in each segment and sector of the left ventricle. 88% of the patient models (n = 24) were classified as a wringing rotation and 12% (n = 3) as a rigid-body-type rotation. Comparison to classification based on in vivo rotational NOGA XP maps showed no correlation. Thus, isolated changes of the electromechanical activation sequence in the left ventricle are not sufficient to reproduce the rotation pattern changes observed in vivo and suggest that further patho-mechanisms are involved.

A Fully-Coupled Electro-Mechanical Whole-Heart Computational Model: Influence of Cardiac Contraction on the ECG

The ECG is one of the most commonly used non-invasive tools to gain insights into the electrical functioning of the heart. It has been crucial as a foundation in the creation and validation of in silico models describing the underlying electrophysiological processes. However, so far, the contraction of the heart and its influences on the ECG have mainly been overlooked in in silico models. As the heart contracts and moves, so do the electrical sources within the heart responsible for the signal on the body surface, thus potentially altering the ECG. To illuminate these aspects, we developed a human 4-chamber electro-mechanically coupled whole heart in silico model and embedded it within a torso model. Our model faithfully reproduces measured 12-lead ECG traces, circulatory characteristics, as well as physiological ventricular rotation and atrioventricular valve plane displacement. We compare our dynamic model to three non-deforming ones in terms of standard clinically used ECG leads (Einthoven and Wilson) and body surface potential maps (BSPM). The non-deforming models consider the heart at its ventricular end-diastatic, end-diastolic and end-systolic states. The standard leads show negligible differences during P-Wave and QRS-Complex, yet during T-Wave the leads closest to the heart show prominent differences in amplitude. When looking at the BSPM, there are no notable differences during the P-Wave, but effects of cardiac motion can be observed already during the QRS-Complex, increasing further during the T-Wave. We conclude that for the modeling of activation (P-Wave/QRS-Complex), the associated effort of simulating a complete electro-mechanical approach is not worth the computational cost. But when looking at ventricular repolarization (T-Wave) in standard leads as well as BSPM, there are areas where the signal can be influenced by cardiac motion of the heart to an extent that should not be ignored.

The openCARP simulation environment for cardiac electrophysiology

BACKGROUND AND OBJECTIVE

Cardiac electrophysiology is a medical specialty with a long and rich tradition of computational modeling. Nevertheless, no community standard for cardiac electrophysiology simulation software has evolved yet. Here, we present the openCARP simulation environment as one solution that could foster the needs of large parts of this community.

METHODS AND RESULTS

openCARP and the Python-based carputils framework allow developing and sharing simulation pipelines which automate in silico experiments including all modeling and simulation steps to increase reproducibility and productivity. The continuously expanding openCARP user community is supported by tailored infrastructure. Documentation and training material facilitate access to this complementary research tool for new users. After a brief historic review, this paper summarizes requirements for a high-usability electrophysiology simulator and describes how openCARP fulfills them. We introduce the openCARP modeling workflow in a multi-scale example of atrial fibrillation simulations on single cell, tissue, organ and body level and finally outline future development potential.

CONCLUSION

As an open simulator, openCARP can advance the computational cardiac electrophysiology field by making state-of-the-art simulations accessible. In combination with the carputils framework, it offers a tailored software solution for the scientific community and contributes towards increasing use, transparency, standardization and reproducibility of in silico experiments.