Patient-Specific Identification of Atrial Flutter Vulnerability-A Computational Approach to Reveal Latent Reentry Pathways

Digital twins in medicine

Medical digital twins, which are potentially vital for personalized medicine, have become a recent focus in medical research. Here we present an overview of the state of the art in medical digital twin development, especially in oncology and cardiology, where it is most advanced. We discuss major challenges, such as data integration and privacy, and provide an outlook on future advancements. Emphasizing the importance of this technology in healthcare, we highlight the potential for substantial improvements in patient-specific treatments and diagnostics.

Transition mechanisms from atrial flutter to atrial fibrillation during anti-tachycardia pacing therapy

BACKGROUND

Atrial anti-tachycardia pacing (aATP) has been shown to be effective for the termination of atrial tachyarrhythmias, but its success rate is still not high enough.

OBJECTIVE

The main objective of this study was to investigate the mechanisms of atrial flutter (AFL) termination by aATP and the transition from AFL to atrial fibrillation (AF) during aATP.

METHODS

We developed a multi-scale model of the human atrium based on magnetic resonance images and examined the atrial electrophysiology of AFL during aATP with a ramp protocol.

RESULTS

In successful cases of aATP, paced excitation entered the excitable gap and collided with the leading edge of the reentrant wave front. Furthermore, the excitation propagating in the opposite direction collided with the trailing edge of the reentrant wave to terminate AFL. The second collision was made possible by the distribution of the wave propagation velocity in the atria. The detailed analysis revealed that the slowing of propagation velocity occurred at the exit of the sub-Eustachian isthmus, probably due to source-sink mismatch. During the transition from AFL to AF, the excitation collided with the refractory zone of the preceding wave and broke into multiple wave fronts to induce AF. A similar observation was made for the transition from AF to sinus rhythm. In both cases, the complex anatomy of the atria played an essential role.

CONCLUSION

The complex anatomy of atria plays an essential role in the maintenance of stable AFL and its termination by aATP, which were revealed by the realistic three-dimensional simulation model.

lifex-ep: a robust and efficient software for cardiac electrophysiology simulations

BACKGROUND

Simulating the cardiac function requires the numerical solution of multi-physics and multi-scale mathematical models. This underscores the need for streamlined, accurate, and high-performance computational tools. Despite the dedicated endeavors of various research teams, comprehensive and user-friendly software programs for cardiac simulations, capable of accurately replicating both normal and pathological conditions, are still in the process of achieving full maturity within the scientific community.

RESULTS

This work introduces [Formula: see text]-ep, a publicly available software for numerical simulations of the electrophysiology activity of the cardiac muscle, under both normal and pathological conditions. [Formula: see text]-ep employs the monodomain equation to model the heart's electrical activity. It incorporates both phenomenological and second-generation ionic models. These models are discretized using the Finite Element method on tetrahedral or hexahedral meshes. Additionally, [Formula: see text]-ep integrates the generation of myocardial fibers based on Laplace-Dirichlet Rule-Based Methods, previously released in Africa et al., 2023, within [Formula: see text]-fiber. As an alternative, users can also choose to import myofibers from a file. This paper provides a concise overview of the mathematical models and numerical methods underlying [Formula: see text]-ep, along with comprehensive implementation details and instructions for users. [Formula: see text]-ep features exceptional parallel speedup, scaling efficiently when using up to thousands of cores, and its implementation has been verified against an established benchmark problem for computational electrophysiology. We showcase the key features of [Formula: see text]-ep through various idealized and realistic simulations conducted in both normal and pathological scenarios. Furthermore, the software offers a user-friendly and flexible interface, simplifying the setup of simulations using self-documenting parameter files.

CONCLUSIONS

[Formula: see text]-ep provides easy access to cardiac electrophysiology simulations for a wide user community. It offers a computational tool that integrates models and accurate methods for simulating cardiac electrophysiology within a high-performance framework, while maintaining a user-friendly interface. [Formula: see text]-ep represents a valuable tool for conducting in silico patient-specific simulations.

MedalCare-XL: 16,900 healthy and pathological synthetic 12 lead ECGs from electrophysiological simulations

Mechanistic cardiac electrophysiology models allow for personalized simulations of the electrical activity in the heart and the ensuing electrocardiogram (ECG) on the body surface. As such, synthetic signals possess known ground truth labels of the underlying disease and can be employed for validation of machine learning ECG analysis tools in addition to clinical signals. Recently, synthetic ECGs were used to enrich sparse clinical data or even replace them completely during training leading to improved performance on real-world clinical test data. We thus generated a novel synthetic database comprising a total of 16,900 12 lead ECGs based on electrophysiological simulations equally distributed into healthy control and 7 pathology classes. The pathological case of myocardial infraction had 6 sub-classes. A comparison of extracted features between the virtual cohort and a publicly available clinical ECG database demonstrated that the synthetic signals represent clinical ECGs for healthy and pathological subpopulations with high fidelity. The ECG database is split into training, validation, and test folds for development and objective assessment of novel machine learning algorithms.

Recommender system for ablation lines to treat complex atrial tachycardia

BACKGROUND AND OBJECTIVE

Planning the optimal ablation strategy for the treatment of complex atrial tachycardia (CAT) is a time consuming task and is error-prone. Recently, directed network mapping, a technology based on graph theory, proved to efficiently identify CAT based solely on data of clinical interventions. Briefly, a directed network was used to model the atrial electrical propagation and reentrant activities were identified by looking for closed-loop paths in the network. In this study, we propose a recommender system, built as an optimization problem, able to suggest the optimal ablation strategy for the treatment of CAT.

METHODS

The optimization problem modeled the optimal ablation strategy as that one interrupting all reentrant mechanisms while minimizing the ablated atrial surface. The problem was designed on top of directed network mapping. Considering the exponential complexity of finding the optimal solution of the problem, we introduced a heuristic algorithm with polynomial complexity. The proposed algorithm was applied to the data of i) 6 simulated scenarios including both left and right atrial flutter; and ii) 10 subjects that underwent a clinical routine.

RESULTS

The recommender system suggested the optimal strategy in 4 out of 6 simulated scenarios. On clinical data, the recommended ablation lines were found satisfactory on 67% of the cases according to the clinician's opinion, while they were correctly located in 89%. The algorithm made use of only data collected during mapping and was able to process them nearly real-time.

CONCLUSIONS

The first recommender system for the identification of the optimal ablation lines for CAT, based solely on the data collected during the intervention, is presented. The study may open up interesting scenarios for the application of graph theory for the treatment of CAT.

Personalized ablation vs. conventional ablation strategies to terminate atrial fibrillation and prevent recurrence

AIMS

The long-term success rate of ablation therapy is still sub-optimal in patients with persistent atrial fibrillation (AF), mostly due to arrhythmia recurrence originating from arrhythmogenic sites outside the pulmonary veins. Computational modelling provides a framework to integrate and augment clinical data, potentially enabling the patient-specific identification of AF mechanisms and of the optimal ablation sites. We developed a technology to tailor ablations in anatomical and functional digital atrial twins of patients with persistent AF aiming to identify the most successful ablation strategy.

METHODS AND RESULTS

Twenty-nine patient-specific computational models integrating clinical information from tomographic imaging and electro-anatomical activation time and voltage maps were generated. Areas sustaining AF were identified by a personalized induction protocol at multiple locations. State-of-the-art anatomical and substrate ablation strategies were compared with our proposed Personalized Ablation Lines (PersonAL) plan, which consists of iteratively targeting emergent high dominant frequency (HDF) regions, to identify the optimal ablation strategy. Localized ablations were connected to the closest non-conductive barrier to prevent recurrence of AF or atrial tachycardia. The first application of the HDF strategy had a success of >98% and isolated only 5-6% of the left atrial myocardium. In contrast, conventional ablation strategies targeting anatomical or structural substrate resulted in isolation of up to 20% of left atrial myocardium. After a second iteration of the HDF strategy, no further arrhythmia episode could be induced in any of the patient-specific models.

CONCLUSION

The novel PersonAL in silico technology allows to unveil all AF-perpetuating areas and personalize ablation by leveraging atrial digital twins.

Quantifying the impact of shape uncertainty on predicted arrhythmias

BACKGROUND

Personalised computer models are increasingly used to diagnose cardiac arrhythmias and tailor treatment. Patient-specific models of the left atrium are often derived from pre-procedural imaging of anatomy and fibrosis. These images contain noise that can affect simulation predictions. There are few computationally tractable methods for propagating uncertainties from images to clinical predictions.

METHOD

We describe the left atrium anatomy using our Bayesian shape model that captures anatomical uncertainty in medical images and has been validated on 63 independent clinical images. This algorithm describes the left atrium anatomy using Nmodes=15 principal components, capturing 95% of the shape variance and calculated from 70 clinical cardiac magnetic resonance (CMR) images. Latent variables encode shape uncertainty: we evaluate their posterior distribution for each new anatomy. We assume a normally distributed prior. We use the unscented transform to sample from the posterior shape distribution. For each sample, we assign the local material properties of the tissue using the projection of late gadolinium enhancement CMR (LGE-CMR) onto the anatomy to estimate local fibrosis. To test which activation patterns an atrium can sustain, we perform an arrhythmia simulation for each sample. We consider 34 possible outcomes (31 macro-re-entries, functional re-entry, atrial fibrillation, and non-sustained arrhythmia). For each sample, we determine the outcome by comparing pre- and post-ablation activation patterns following a cross-field stimulus.

RESULTS

We create patient-specific atrial electrophysiology models of ten patients. We validate the mean and standard deviation maps from the unscented transform with the same statistics obtained with 12,000 Monte Carlo (ground truth) samples. We found discrepancies <3% and <2% for the mean and standard deviation for fibrosis burden and activation time, respectively. For each patient case, we then compare the predicted outcome from a model built on the clinical data (deterministic approach) with the probability distribution obtained from the simulated samples. We found that the deterministic approach did not predict the most likely outcome in 80% of the cases. Finally, we estimate the influence of each source of uncertainty independently. Fixing the anatomy to the posterior mean and maintaining uncertainty in fibrosis reduced the prediction of self-terminating arrhythmias from ≃14% to ≃7%. Keeping the fibrosis fixed to the sample mean while retaining uncertainty in shape decreased the prediction of substrate-driven arrhythmias from ≃33% to ≃18% and increased the prediction of macro-re-entries from ≃54% to ≃68%.

CONCLUSIONS

We presented a novel method for propagating shape uncertainty in atrial models through to uncertainty in numerical simulations. The algorithm takes advantage of the unscented transform to compute the output distribution of the outcomes. We validated the unscented transform as a viable sampling strategy to deal with anatomy uncertainty. We then showed that the prediction computed with a deterministic model does not always coincide with the most likely outcome. Finally, we found that shape uncertainty affects the predictions of macro-re-entries, while fibrosis uncertainty affects the predictions of functional re-entries.

Calibrating cardiac electrophysiology models using latent Gaussian processes on atrial manifolds

Models of electrical excitation and recovery in the heart have become increasingly detailed, but have yet to be used routinely in the clinical setting to guide personalized intervention in patients. One of the main challenges is calibrating models from the limited measurements that can be made in a patient during a standard clinical procedure. In this work, we propose a novel framework for the probabilistic calibration of electrophysiology parameters on the left atrium of the heart using local measurements of cardiac excitability. Parameter fields are represented as Gaussian processes on manifolds and are linked to measurements via surrogate functions that map from local parameter values to measurements. The posterior distribution of parameter fields is then obtained. We show that our method can recover parameter fields used to generate localised synthetic measurements of effective refractory period. Our methodology is applicable to other measurement types collected with clinical protocols, and more generally for calibration where model parameters vary over a manifold.

Case report: Personalized computational model guided ablation for left atrial flutter

Atypical atrial flutter is seen post-ablation in patients, and it can be challenging to map. These flutters are typically set up around areas of scar in the left atrium. MRI can reliably identify left atrial scar. We propose a personalized computational model using patient specific scar information, to generate a monodomain model. In the model conductivities are adjusted for different tissue regions and flutter was induced with a premature pacing protocol. The model was tested prospectively in patients undergoing atypical flutter ablation. The simulation-predicted flutters were visualized and presented to clinicians. Validation of the computational model was motivated by recording from electroanatomical mapping. These personalized models successfully predicted clinically observed atypical flutter circuits and at times even better than invasive maps leading to flutter termination at isthmus sites predicted by the model.

An Automata-Based Cardiac Electrophysiology Simulator to Assess Arrhythmia Inducibility

Personalized cardiac electrophysiology simulations have demonstrated great potential to study cardiac arrhythmias and help in therapy planning of radio-frequency ablation. Its application to analyze vulnerability to ventricular tachycardia and sudden cardiac death in infarcted patients has been recently explored. However, the detailed multi-scale biophysical simulations used in these studies are very demanding in terms of memory and computational resources, which prevents their clinical translation. In this work, we present a fast phenomenological system based on cellular automata (CA) to simulate personalized cardiac electrophysiology. The system is trained on biophysical simulations to reproduce cellular and tissue dynamics in healthy and pathological conditions, including action potential restitution, conduction velocity restitution and cell safety factor. We show that a full ventricular simulation can be performed in the order of seconds, emulate the results of a biophysical simulation and reproduce a patient’s ventricular tachycardia in a model that includes a heterogeneous scar region. The system could be used to study the risk of arrhythmia in infarcted patients for a large number of scenarios.

Fast Characterization of Inducible Regions of Atrial Fibrillation Models With Multi-Fidelity Gaussian Process Classification

Computational models of atrial fibrillation have successfully been used to predict optimal ablation sites. A critical step to assess the effect of an ablation pattern is to pace the model from different, potentially random, locations to determine whether arrhythmias can be induced in the atria. In this work, we propose to use multi-fidelity Gaussian process classification on Riemannian manifolds to efficiently determine the regions in the atria where arrhythmias are inducible. We build a probabilistic classifier that operates directly on the atrial surface. We take advantage of lower resolution models to explore the atrial surface and combine seamlessly with high-resolution models to identify regions of inducibility. We test our methodology in 9 different cases, with different levels of fibrosis and ablation treatments, totalling 1,800 high resolution and 900 low resolution simulations of atrial fibrillation. When trained with 40 samples, our multi-fidelity classifier that combines low and high resolution models, shows a balanced accuracy that is, on average, 5.7% higher than a nearest neighbor classifier. We hope that this new technique will allow faster and more precise clinical applications of computational models for atrial fibrillation. All data and code accompanying this manuscript will be made publicly available at: https://github.com/fsahli/AtrialMFclass.

Ablation of Left Atrial Tachycardia following Catheter Ablation of Atrial Fibrillation: 12-Month Success Rates

The treatment of atrial tachycardia following catheter ablation of atrial fibrillation is often challenging. Electrophysiological studies using high-resolution 3D mapping systems have contributed significantly to their understanding, and new ablation approaches have shown high rates of acute terminations with low recurrences for the clinical AT. However, patient populations are very heterogeneous, and long-term data of the freedom from any atrial tachycardia or any arrhythmia are still sparse. To evaluate long-term success, a unified patient population and predefined ablation strategies are preferred. In this study, we present 12-month success and mean 30 month follow-up data of catheter ablation of left atrial tachycardia. All 35 patients had a history of pulmonary vein isolation (PVI), 71% of which had a previous substrate modification. A total of 54 ATs, with a mean cycle length 297 ± 86 ms, 31 macro-reentries, and 4 localized reentries, were targeted. The ablation strategy to be used was given by the study protocol, depending on the type of reentry and the number of critical isthmuses. All available ablation strategies were included: standard (anatomical) lines, individual lines, critical isthmuses, and focal ablation. All ATs were terminated by ablation. A total of 91% terminated upon the first ablation strategy. Freedom from any AT after 12 months was 82%, and from any arrhythmia, it was 77%. The multi-procedure success after 30 months was 65% for any AT and 55% for any arrhythmia. In conclusion, individual ablation strategies based on the reentry mechanism and the number of critical isthmuses seems promising and demonstrates a high long-term clinical success. Tachycardia comprising a single critical isthmus can be ablated by critical isthmus ablation only. These patients present with the highest 12-month and long-term success rates.

Hybrid machine learning to localize atrial flutter substrates using the surface 12-lead electrocardiogram

Aims

Atrial flutter (AFlut) is a common re-entrant atrial tachycardia driven by self-sustainable mechanisms that cause excitations to propagate along pathways different from sinus rhythm. Intra-cardiac electrophysiological mapping and catheter ablation are often performed without detailed prior knowledge of the mechanism perpetuating AFlut, likely prolonging the procedure time of these invasive interventions. We sought to discriminate the AFlut location [cavotricuspid isthmus-dependent (CTI), peri-mitral, and other left atrium (LA) AFlut classes] with a machine learning-based algorithm using only the non-invasive signals from the 12-lead electrocardiogram (ECG).

Methods and results

Hybrid 12-lead ECG dataset of 1769 signals was used (1424 in silico ECGs, and 345 clinical ECGs from 115 patients—three different ECG segments over time were extracted from each patient corresponding to single AFlut cycles). Seventy-seven features were extracted. A decision tree classifier with a hold-out classification approach was trained, validated, and tested on the dataset randomly split after selecting the most informative features. The clinical test set comprised 38 patients (114 clinical ECGs). The classifier yielded 76.3% accuracy on the clinical test set with a sensitivity of 89.7%, 75.0%, and 64.1% and a positive predictive value of 71.4%, 75.0%, and 86.2% for CTI, peri-mitral, and other LA class, respectively. Considering majority vote of the three segments taken from each patient, the CTI class was correctly classified at 92%.

Conclusion

Our results show that a machine learning classifier relying only on non-invasive signals can potentially identify the location of AFlut mechanisms. This method could aid in planning and tailoring patient-specific AFlut treatments.

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.

A bi-atrial statistical shape model for large-scale in silico studies of human atria: Model development and application to ECG simulations

Large-scale electrophysiological simulations to obtain electrocardiograms (ECG) carry the potential to produce extensive datasets for training of machine learning classifiers to, e.g., discriminate between different cardiac pathologies. The adoption of simulations for these purposes is limited due to a lack of ready-to-use models covering atrial anatomical variability. We built a bi-atrial statistical shape model (SSM) of the endocardial wall based on 47 segmented human CT and MRI datasets using Gaussian process morphable models. Generalization, specificity, and compactness metrics were evaluated. The SSM was applied to simulate atrial ECGs in 100 random volumetric instances. The first eigenmode of our SSM reflects a change of the total volume of both atria, the second the asymmetry between left vs. right atrial volume, the third a change in the prominence of the atrial appendages. The SSM is capable of generalizing well to unseen geometries and 95% of the total shape variance is covered by its first 24 eigenvectors. The P waves in the 12-lead ECG of 100 random instances showed a duration of 109.7±12.2 ms in accordance with large cohort studies. The novel bi-atrial SSM itself as well as 100 exemplary instances with rule-based augmentation of atrial wall thickness, fiber orientation, inter-atrial bridges and tags for anatomical structures have been made publicly available. This novel, openly available bi-atrial SSM can in future be employed to generate large sets of realistic atrial geometries as a basis for in silico big data approaches.

Enabling forward uncertainty quantification and sensitivity analysis in cardiac electrophysiology by reduced order modeling and machine learning

We present a new, computationally efficient framework to perform forward uncertainty quantification (UQ) in cardiac electrophysiology. We consider the monodomain model to describe the electrical activity in the cardiac tissue, coupled with the Aliev-Panfilov model to characterize the ionic activity through the cell membrane. We address a complete forward UQ pipeline, including both: (i) a variance-based global sensitivity analysis for the selection of the most relevant input parameters, and (ii) a way to perform uncertainty propagation to investigate the impact of intra-subject variability on outputs of interest depending on the cardiac potential. Both tasks exploit stochastic sampling techniques, thus implying overwhelming computational costs because of the huge amount of queries to the high-fidelity, full-order computational model obtained by approximating the coupled monodomain/Aliev-Panfilov system through the finite element method. To mitigate this computational burden, we replace the full-order model with computationally inexpensive projection-based reduced-order models (ROMs) aimed at reducing the state-space dimensionality. Resulting approximation errors on the outputs of interest are finally taken into account through artificial neural network (ANN)-based models, enhancing the accuracy of the whole UQ pipeline. Numerical results show that the proposed physics-based ROMs outperform regression-based emulators relying on ANNs built with the same amount of training data, in terms of both numerical accuracy and overall computational efficiency.

Electro-Mechanical Whole-Heart Digital Twins: A Fully Coupled Multi-Physics Approach

Mathematical models of the human heart are evolving to become a cornerstone of precision medicine and support clinical decision making by providing a powerful tool to understand the mechanisms underlying pathophysiological conditions. In this study, we present a detailed mathematical description of a fully coupled multi-scale model of the human heart, including electrophysiology, mechanics, and a closed-loop model of circulation. State-of-the-art models based on human physiology are used to describe membrane kinetics, excitation-contraction coupling and active tension generation in the atria and the ventricles. Furthermore, we highlight ways to adapt this framework to patient specific measurements to build digital twins. The validity of the model is demonstrated through simulations on a personalized whole heart geometry based on magnetic resonance imaging data of a healthy volunteer. Additionally, the fully coupled model was employed to evaluate the effects of a typical atrial ablation scar on the cardiovascular system. With this work, we provide an adaptable multi-scale model that allows a comprehensive personalization from ion channels to the organ level enabling digital twin modeling.

A Reproducible Protocol to Assess Arrhythmia Vulnerability in silico: Pacing at the End of the Effective Refractory Period

In both clinical and computational studies, different pacing protocols are used to induce arrhythmia and non-inducibility is often considered as the endpoint of treatment. The need for a standardized methodology is urgent since the choice of the protocol used to induce arrhythmia could lead to contrasting results, e.g., in assessing atrial fibrillation (AF) vulnerabilty. Therefore, we propose a novel method-pacing at the end of the effective refractory period (PEERP)-and compare it to state-of-the-art protocols, such as phase singularity distribution (PSD) and rapid pacing (RP) in a computational study. All methods were tested by pacing from evenly distributed endocardial points at 1 cm inter-point distance in two bi-atrial geometries. Seven different atrial models were implemented: five cases without specific AF-induced remodeling but with decreasing global conduction velocity and two persistent AF cases with an increasing amount of fibrosis resembling different substrate remodeling stages. Compared with PSD and RP, PEERP induced a larger variety of arrhythmia complexity requiring, on average, only 2.7 extra-stimuli and 3 s of simulation time to initiate reentry. Moreover, PEERP and PSD were the protocols which unveiled a larger number of areas vulnerable to sustain stable long living reentries compared to RP. Finally, PEERP can foster standardization and reproducibility, since, in contrast to the other protocols, it is a parameter-free method. Furthermore, we discuss its clinical applicability. We conclude that the choice of the inducing protocol has an influence on both initiation and maintenance of AF and we propose and provide PEERP as a reproducible method to assess arrhythmia vulnerability.

Electrophysiological and anatomical factors determine arrhythmic risk in acute myocardial ischaemia and its modulation by sodium current availability

Acute myocardial ischaemia caused by coronary artery disease is one of the main causes of sudden cardiac death. Even though sodium current blockers are used as anti-arrhythmic drugs, decreased sodium current availability, also caused by mutations, has been shown to increase arrhythmic risk in ischaemic patients. The mechanisms are still unclear. Our goal is to exploit perfect control and data transparency of over 300 high-performance computing simulations to investigate arrhythmia mechanisms in acute myocardial ischaemia with variable sodium current availability. The human anatomically based torso-biventricular electrophysiological model used includes representation of realistic ventricular anatomy and fibre architecture, as well as ionic to electrocardiographic properties. Simulations show that reduced sodium current availability increased arrhythmic risk in acute regional ischaemia due to both electrophysiological (increased dispersion of refractoriness across the ischaemic border zone) and anatomical factors (conduction block from the thin right ventricle to thick left ventricle). The asymmetric ventricular anatomy caused high arrhythmic risk specifically for ectopic stimuli originating from the right ventricle and ventricular base. Increased sodium current availability was ineffective in reducing arrhythmic risk for septo-basal ectopic excitation. Human-based multiscale modelling and simulations reveal key electrophysiological and anatomical factors determining arrhythmic risk in acute ischaemia with variable sodium current availability.

Gaussian process manifold interpolation for probabilistic atrial activation maps and uncertain conduction velocity

In patients with atrial fibrillation, local activation time (LAT) maps are routinely used for characterizing patient pathophysiology. The gradient of LAT maps can be used to calculate conduction velocity (CV), which directly relates to material conductivity and may provide an important measure of atrial substrate properties. Including uncertainty in CV calculations would help with interpreting the reliability of these measurements. Here, we build upon a recent insight into reduced-rank Gaussian processes (GPs) to perform probabilistic interpolation of uncertain LAT directly on human atrial manifolds. Our Gaussian process manifold interpolation (GPMI) method accounts for the topology of the atrium, and allows for calculation of statistics for predicted CV. We demonstrate our method on two clinical cases, and perform validation against a simulated ground truth. CV uncertainty depends on data density, wave propagation direction and CV magnitude. GPMI is suitable for probabilistic interpolation of other uncertain quantities on non-Euclidean manifolds. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Challenges Associated with Interpreting Mechanisms of AF

Determining optimal treatment strategies for complex arrhythmogenesis in AF is confounded by the lack of consensus regarding the mechanisms causing AF. Studies report different mechanisms for AF, ranging from hierarchical drivers to anarchical multiple activation wavelets. Differences in the assessment of AF mechanisms are likely due to AF being recorded across diverse models using different investigational tools, spatial scales and clinical populations. The authors review different AF mechanisms, including anatomical and functional re-entry, hierarchical drivers and anarchical multiple wavelets. They then describe different cardiac mapping techniques and analysis tools, including activation mapping, phase mapping and fibrosis identification. They explain and review different data challenges, including differences between recording devices in spatial and temporal resolutions, spatial coverage and recording surface, and report clinical outcomes using different data modalities. They suggest future research directions for investigating the mechanisms underlying human AF.

Sensitivity of Ablation Targets Prediction to Electrophysiological Parameter Variability in Image-Based Computational Models of Ventricular Tachycardia in Post-infarction Patients

Ventricular tachycardia (VT), which could lead to sudden cardiac death, occurs frequently in patients with myocardial infarction. Computational modeling has emerged as a powerful platform for the non-invasive investigation of lethal heart rhythm disorders in post-infarction patients and for guiding patient VT ablation. However, it remains unclear how VT dynamics and predicted ablation targets are influenced by inter-patient variability in action potential duration (APD) and conduction velocity (CV). The goal of this study was to systematically assess the effect of changes in the electrophysiological parameters on the induced VTs and predicted ablation targets in personalized models of post-infarction hearts. Simulations were conducted in 5 patient-specific left ventricular models reconstructed from late gadolinium-enhanced magnetic resonance imaging scans. We comprehensively characterized all possible pre-ablation and post-ablation VTs in simulations conducted with either an “average human VT”-based electrophysiological representation (i.e., EP_avg) or with ±10% APD or CV (i.e., EP_var); additional simulations were also executed in some models for an extended range of these parameters. The results showed that: (1) a subset of reentries (76.2–100%, depending on EP parameter set) conducted with ±10% APD/CV was observed in approximately the same locations as reentries observed in EP_avg cases; (2) emergent VTs could be induced sometimes after ablation in EP_avg models, and these emergent VTs often corresponded to the pre-ablation reentries in simulations with EP_var parameter sets. These findings demonstrate that the VT ablation target uncertainty in patient-specific ventricular models with an average representation of VT-remodeled electrophysiology is relatively low and the ablation targets stable, as the localization of the induced VTs was primarily driven by the remodeled structural substrate. Thus, personalized ventricular modeling with an average representation of infarct-remodeled electrophysiology may uncover most targets for VT ablation.